Quality by design (QbD) based fabrication of atazanavir loaded nanostructured lipid carriers for lymph targeting: Bioavailability enhancement using chylomicron flow block model and toxicity study


 Atazanavir (ATV) is widely used as anti-HIV agent with poor aqueous solubility which requires fabrication of novel drug delivery system to enhance therapeutic activity and safety. For this purpose, the quality by design (QbD) based ATV loaded nanostructured lipid carriers (NLCs) to address the challenges of bioavailability and its safety on oral administration. Herein, the main objective was to identify the influencing variables for the production of quality product. Considering this objective, quality target product profile (QTPP) was assigned and a systematic risk assessment study was performed to identify the critical material attributes (CMAs) and critical process parameter (CPP) having an influence on critical quality attributes (CQAs). Lipid concentrations, surfactant concentrations, and pressure of high-pressure homogenizer were identified as CMAs and CPP. ATV-NLCs were prepared by emulsification-high pressure homogenization method and further lyophilized to obtain solid-state NLCs. The effect of formulation variables (CMAs and CPP) on responses like particle size (Y1), polydispersity index (Y2), and zeta potential (Y3) was observed by central composite rotatable design (CCRD). The data were statistically evaluated by ANOVA for confirmation of a significant level (P<0.05). The optimal conditions of NLCs were obtained by generating design space and desirability value. The lyophilized ATV-NLCs were characterized by DSC, PXRD, and FT-IR analysis. The morphology of NLCs was revealed by TEM and FESEM. In vitro study suggested a sustained release pattern of drug (92.37±1.03 %) with a mechanism of Korsmeyer-Peppas model (r2 =0.925, and n=0.63). In vivo evaluation in Wistar rats showed significantly higher (p<0.001) plasma drug concentration of ATV-NLCs as compared to ATV-suspension using chylomicron flow block model. The relative bioavailability of ATV-NLCs was obtained to be 2.54 folds. Thus, a safe and promising drug targeting system was successfully developed to improve bioavailability and avoiding first-pass effect ensures to circumvent the acute-toxicity of liver.


Introduction
To date, antiretroviral (ARV) drugs have been made remarkable progress in the treatment and control of human immunodeficiency virus (HIV) infection [1]. The combined therapy of ARV drugs so-called highly active antiretroviral therapy (HAART) is effective to extend the lifespan of individuals suffered from HIV infection [2,3]. However, they are inefficient to eradicate HIV completely from the host owing to their poor biodistribution in certain inaccessible compartments of the body such as CNS, genitals, and lymphatics [2]. Therefore, it is noteworthy to develop a novel strategy using efficient modalities for the treatment of HIV. Many scientists have believed that targeting lymphatics can be improved bioavailability of drug molecules [4].
Atazanavir (ATV), which belongs to the category of azapeptide protease inhibitor was selected as a model drug for the present investigation. ATV is categorized as a BCS-II drug exhibiting poor aqueous solubility and high permeability due to high log p value (log p=4. 11). It is the most recommended drug of choice in HAART regimen for treatment of HIV [5]. ATV is marketed as conventional therapeutic dose of 400 mg for adults. Regardless of promising molecule, ATV shows several biopharmaceutical and toxicological issues including rapid first-pass metabolism and P-gp efflux leads to a marked reduction in oral bioavailability to 60 % [6]. Additionally, low absorption and frequent administration causes variation in plasma concentration and produces serious hepatic toxicity that limits clinical applications [7]. To overcome the aforementioned issues, there is a strong need to develop a novel strategy of ATV delivery to improve oral bioavailability and lowers hepatic toxicity.
Few works have been accounted for ATV to improve bioavailability including solid lipid nanoparticles (SNLs) [8,9], self-nanoemulsifying drug delivery systems (SNEDDS) [5], solid dispersion [10], and nanoparticles of eudragit RL 100 [6]. However, no report has been 5 addressed to target the lymphatics to improve bioavailability and subsequently safety. In this context, lipid-based delivery system holds the assurance to access the drugs into the intestinal lymphatics in spite of blood portal avoiding first-pass metabolism contributing to bioavailability improvement [11][12][13][14][15].
Nanostructured lipid carriers (NLCs) have been introduced as an option for the traditional carriers, such as liposomes, emulsions, microparticles, solid lipid nanoparticles (SLNs), and their polymeric counterparts due to numerous advantages [16,17]. NLCs are referred to as nextgeneration and better than SLNs because of some advantages including higher drug entrapment, avoidance of drug expulsion, excellent biocompatibility, non-toxicity, shielding the drug from enzymatic degradation, controlled drug release, and industrial scalability [18]. NLCs are the combination of solid lipid and liquid lipid (oil) of GRAS (generally regarded as safe) status that makes the matrix system imperfect to entrap a higher amount of drugs [17,19,20]. NLCs having a smaller particle size < 200 nm reduces the hepatic uptake and extends blood circulation time [21]. More specifically smaller particle size and minimum polydispersity index (PDI) are the common attributes of NLCs for performing the targeting activity. Several researchers have been developed nanocarriers to target the lymphatics by anti-HIV drugs such as ritonavir, lopinavir, and efavirenz [1,15,22,23].
The development of NLCs involves multistage production and multiple variables including lipid concentrations, surfactant concentrations, amount of organic solvents, speed of high-speed homogenizer, and pressure of high pressure homogenizer. Controlling all the variables using several experiments is difficult. Therefore, a systematic approach is required to achieve the desired product. Thus, the application of quality by design (QbD) paradigm helps to identify, understand, and control the variables influencing the results [24]. QbD approach consigns the 7 of 10 °C higher than the melting point of solid lipid. Various solid lipids such as Precirol ® ATO5, Compritol ® ATO 888, and Geleol™ Mono and Diglyceride (GMDG) and Compritol ® HD5 ATO, Gellucire ® 44/14, Palmitic acid, and Myristic acid were screened visually according to the solubilization behavior of ATV.
Likewise, the liquid lipid (oil) was screened by the continuous addition method. Briefly, an excess amount of ATV was added in a fixed amount (2 mL) of oil to achieve a saturation point.
Each sample was centrifuged and analyzed by UV spectrophotometrically at 249 nm [28] The compatibility ratio of lipid mixture was selected based on phase separation parameter. The different mixture of solid lipid and liquid lipid in the proportion of 9:1, 8:2, 7:3, 6:4, 5:5, 6:4, and 7:3 was taken and heated at 5 °C above the melting point of solid lipid and was evaluated for the phase separation after 1 h, 12 h and 24 h of solidification. No sign of phase separation was observed in the ratio of lipid mixture that was selected throughout the study [29]. The surfactant and co-surfactant were selected based on solubility and HLB values which lead to required particle size, PDI, and stability [30].

Assignment of QTPP and CQA
Implementation of QbD based concept was initiated by assigning the QTPP which can be achieved with the view of patient-centric approach considering safety and efficacy of final formulation [31]. It was based on the recommendation of ICH Q8 guidelines, experimental judgment, literature, and prior knowledge [32]. The essential parameters concerned with the 8 route of administration, strength, dosage form, pharmacokinetics, and release profile were especially focused on the production of targeted product.
Moreover, the significant CQAs were identified which have a greater impact on the quality of product. The identified CQAs was monitored and controlled to ensure the desired quality of the product [33]. These attributes were further explored based on previous experience, and knowledge gained through literature [34].

Risk assessment
Risk assessment was performed to identify critical formulation variables like critical material attributes (CMAs) and critical process parameters (CPPs) which have a decisive effect on proposed product performance [34]. The Ishikawa diagram ( Fig. 1) was constructed to enlist potential high-risk variables with their cause and effect to achieve the quality of the final product [35]. Furthermore, the risk assessment matrix (RAM) method was employed to find the effective and potential formulation variables. Thus, risk assessment of the proposed formulation was executed by combining the Ishikawa diagram and RAM method [36].

HPLC analysis
The ATV content present in formulation, plasma, spleen, and thymus were analyzed by the modified RP-HPLC method previously reported by Singh and Pai, 2016 [6]. The mobile phase was prepared using potassium dihydrogen phosphate (pH 3.4) and acetonitrile in a ratio of (60:40, v/v). The sample (20 µL) for ATV analysis and bioanalytical sample (20 µL) was injected into Nucleosil C18 (4.6 mm I. D. × 250 mm) column at ambient temperature. The flow of 1.2 mL/min was adjusted for eluting ATV and was monitored at 249 nm by a UV detector.
Furthermore, the developed method was validated as per the International Conference on 9 Harmonization (ICH) guidelines, Q2 (R1) (2005) [37]. The standard curve was found to be linear (r 2 = 0.991) over the range of 0.050-2.0 µg/mL. Before analysis of sample, standard solutions of ATV and saquinavir (SAQ) Internal standard (IS) of strength 1000 μg/mL were prepared separately. A freshly collected rat plasma (blank) of approximately 0.5 mL was spiked into ATV and SAQ (IS) solutions of strength 1000 μg/mL. To this mixture, ethyl acetate was added and vortexed for 2 min in order to extract the drug content.
Subsequently, the supernatant was obtained by centrifugation at 3000 rpm for 10 min and filtered through the 0.22 μm membrane filter. After evaporation, the dry residue was reconstituted with the mobile phase and the resulting solution of a 20 µL sample was injected into HPLC system.

Preparation of ATV loaded NLCs (ATV-NLCs)
ATV-NLCs were prepared by the emulsification-hot high pressure homogenization (HPH) technique [38]. Briefly, a measured quantity of GMDG (M. P. 54.0 -64.0 °C), Capryol TM 90 was added in 10 mL chloroform and melted on a water bath (Equitron water bath). To this phase, a measured amount of ATV (0.04 %) and Span ® 20 (co-surfactant) were dissolved. The organic solvent was completely evaporated by heating on a water bath at 70-75 ºC. On the other hand, Tween ® 80 (surfactant) was dissolved in distilled water and heated up to the same temperature as that of the lipid phase. After a clear and homogeneous lipid phase, the hot surfactant solution was dispersed into a melted lipid phase to obtain a pre-emulsion. This pre-emulsion was homogenized by high-speed homogenizer (Ultra-Turrax, T25 Basic, Ika Werke, Stanfer, Germany) at 10,000 rpm for 10 min. It was then passed through a hot high-pressure homogenizer (HPH) (PANDA 2K, Niro Soavi, Italy) at 600-800 Bar pressure and 5-7 cycles to produce the dispersion. The resultant dispersion was allowed to cool at room temperature to obtain NLCs and was stored in a cool place until characterizations [39].

Central composite rotatable design (CCRD) and statistical optimization
The CCRD was employed using Design Expert® version 12 (DX12) software for the influence of different variables on responses. As summarized in Table 1, the formulation variables such as lipid concentrations (A), surfactant concentrations (B), and HPH pressure (C) were studied at five different levels as low (−1), medium (0), high (+1), relatively low (−α) and relatively high (+α). The responses were selected as particle size (Y1), PDI (Y2), and zeta potential (Y3). Thus, 20 experimental runs including 2 axial points (±α), 8 factorial design (±1), and 10 central points (0) were studied for evaluation of variables as shown in Table 2. Further, the responses were examined for interactions and 3D response plots were generated.
To achieve optimal formulation, the graphical and numerical analysis was performed based on constraints like minimum particle size, minimum PDI, and optimum zeta potential value [40]. A compromised solution for optimization was obtained using a desirability value that lies between 0 and 1 [41]. The observed result was compared with the predicted values and bias (%) was determined in order to validate the selected experimental domain.

Size, PDI, and zeta potential analysis
A zeta sizer (Nano ZS90, Malvern, Worcestershire, UK) was used to measure particle size and PDI of NLCs [42]. The formulation was dispersed in Milli-Q water (aqueous medium) in the ratio of 1:100 and analyzed at 25 °C in triplicate. Zeta potential was measured by electrophoretic mobility of particles using zeta sizer [43]. NLCs were diluted 100 times using distilled Milli-Q water. The measurements were repeated in triplicate at 25 °C to obtain the average zeta potential value.
Determination of percent entrapment efficiency (EE) and drug loading 11 The percent drug entrapment of NLCs was used to determined by a premier tabletop ultracentrifuge (Optima Max-XP, Beckman Coulter Inc., Brea, USA). Briefly, a measured amount of dispersion of ATV-NLCs was centrifuged at 10,000 × g at 4 °C for 15 min. to separate the supernatant and nanocarrier pellets. The supernatant (1 mL) were suitably diluted analyzed by spectrophotometrically (UV-1900, Shimadzu, Tokyo, Japan) at 249 nm. The percentage of drug entrapment was obtained in triplicates and results were reported accordingly. EE (%) and drug loading (% DL) was calculated by mathematical equations:

Lyophilization process
The optimized dispersion was lyophilized (Virtis-Bench Top Lyophilizer, Spinco Biotech Pvt. Ltd.) at a condition of −74 °C temperature, 35 mT vacuum pressure for 42 h [44]. The different quantities of mannitol (2 %, 2.5 %, 3%, and 3.5%) were tested to find the optimal level to obtain a free-flowing solid powder and prevent lysis of NLCs.

Drug content (%)
An accurate 100 µg/mL sample of ATV-NLC was prepared and diluted suitably using the mobile phase. The sample was sonicated followed by filtering through a 0.22 µ membrane filter.
Subsequently, the sample was injected into the HPLC system. The content was calculated by the following mathematical equation.

Differential scanning calorimetry (DSC)
Thermal behavior of the ATV (pure drug), GMDG (solid lipid), physical mixture (PM) of ATV+GMDG, and optimized NLC formulation was recorded by DSC (Mettler Toledo, Italy) enabled with Mettler STARe V8.10 software. The test sample (2 mg) was sealed in an aluminum pan whereas an empty pan was used as a reference [45]. The onset temperature, end set, and peak temperature were recorded at a temperature range of 40 to 200 °C using a heating rate of 10 °C per min in a nitrogen environment.

Powder X-ray diffraction (PXRD)
Diffraction patterns of crystalline substances were determined using PXRD study. This method was used to detect alterations in crystalline forms during development. In this study, the sample of ATV, GMDG (solid lipid), ATV+GMDG (PM), and optimized NLCs were analyzed through the X-ray diffractometer (Brucker AXS D8 advance ® , Karlsruhe, Germany). Each sample was scanned at a 2θ scale with scanning range from 2 to 50°C using voltage 40.0 kV and 35 mA current intensity.

Fourier-transform infrared (FTIR) spectroscopy
The IR spectra of pure ATV, GMDG (solid lipid), ATV+GMDG (PM), and optimized NLCs were recorded using FTIR spectrometer-430 (Shimadzu 8400S, Tokyo, Japan). Samples were physically mixed with dried potassium bromide (KBr) and applied a pressure to obtain a thin pellet. The wavelength for the scanning of KBr pellets was varied from 4000-500 cm −1 .
Transmission electron microscopy (TEM) and field emission scanning electron microscope (FESEM) 13 For surface morphology, the sample was prepared using negatively stained by phosphotungstic acid (2 wt %) and positively stained by uranyl acetate (1 wt %) with NLC dispersion on a copper grid. Following the air dry of sample, the image was captured using a TEM (TEM, Jeol/JEM 2100) instrument at a voltage 200 kV [46].
The surface structure of optimized nanocarriers was studied by FESEM (Quanta 250 FEG; FEI, Switzerland). The aqueous suspension of the NLCs (20 µL) was spread on clean glass slides and allowed to air-dry overnight. The dried glass slides were sputtered with 2 nm thick gold and palladium layer at 20 mA before imaging [47].

Accelerated stability studies
The optimized ATV-NLCs were subjected to accelerated stability study for 90 days as per ICH guidelines at 4 °C/75 % RH and 25 °C ± 2 °C/60 % RH ± 5 % RH. The lyophilized sample was sealed in a glass vial and stored in a stability chamber (CHM-10S, REMI Instruments Ltd., Mumbai, India) and examined for particle size, PDI, zeta potential, and % EE at 0 th , 45 th , and 90 th days respectively [48].

In vitro study and mechanism of release kinetics
In vitro drug release study for optimized ATV-NLCs and ATV-suspension was performed using the dialysis bag (soaked overnight in PBS solution to activate before use). A fixed volume of ATV-NLCs and drug suspension (containing 0.5 % carboxymethyl cellulose) were sealed in the dialysis bag separately. The bag was hanged in 200 mL simulated gastric fluid (pH 1.2) for the first 2 h followed by simulated intestinal gastric fluid (pH 6.8) for a period of 24 h with a stirrer at 50 rpm at 37 °C maintaining a sink condition. Samples were withdrawn at various point of interval and analyzed spectrophotometrically at 249 nm after suitable dilution [49].

Chylomicron flow block model for bioavailability and tissue distribution study
The experiment was designed as per the protocol approved by the Institutional Animal Ethics Healthy Wistar male rats weighing in the range of 220-230 ± 5 g were randomly divided into three cohorts (each having 5 rats in a cage) as follows.
Cohort-I was given oral aqueous ATV-suspension,

Cohort-II was administered ATV-NLC, and
Cohort-III was administered ATV-NLC after the cycloheximide treatment.
The animals were administered a single oral dose of 7.2 mg/kg of their body weight using oral feeding needle [6]. The cohort-III was pre-treated with an intraperitoneal injection of cycloheximide (3 mg/kg) dissolved in saline solution (3 mg/mL) one hour before the experiment.
The cycloheximide, a chylomicron flow blocker, prevents the entry of a drug into the lymphatic system. Rats were anesthetized by intraperitoneal injection of 1 mL/kg of ketamine-xylazine 15 (9%:1%) solution respectively [51]. Afterward, plasma was separated by centrifugation at 3500 rpm for 15 min and analyzed for ATV content using saquinavir (SAQ) as an internal standard (IS) by RP-HPLC with modifications in the reported method [6].
To investigate the tissue biodistribution, the rats were sacrificed by cervical dislocation technique at each time point. Surgical procedures were undertaken as described by Ahammed et al., [23]. The organs of our interest such as the spleen and thymus (containing a maximum amount of lymphocyte) were isolated and cut into pieces. The homogenate was prepared in methanol using a high-speed homogenizer (Ultra-Turrax, T25 Basic, Ika Werke, Stanfer, Germany). These samples were stored at −20 °C until analysis.
For organ samples, the homogenate sample was subjected to the extraction process (described above for the plasma sample). Further this sample was centrifuged at 15000×g for 15 min at 4ºC to obtain a supernatant. The obtained supernatant was evaporated to get the dry residue and analyzed for drug content by estimated RP-HPLC method.

Calculation of pharmacokinetic parameters, relative bioavailability, and statistical analysis
The pharmacokinetic parameters were established using a Microsoft excel based PK solver 2.0 add-in program. Non-compartmental analysis was performed to determine the pharmacokinetic parameters. The AUC0-t (t=24) was established by the trapezoidal method and relative bioavailability was calculated using the following equation [53,54].
GraphPad Prism 6.0 (GraphPad Software Inc., California, USA) was used for statistical data evaluation. All the results were reported in mean ± S.D., (n = 3). The statistical analysis were performed by Student's t-test and one-way ANOVA (analysis of variance) followed by Dunnett's multiple comparisons. The data were considered to be statistically significant (P<0.05).

Oral sub-acute toxicity study and serum analysis
The oral sub-acute toxicity study was performed on selected male Wistar rats. They were divided into control and ATV-NLC administration cohorts. Animals were administered a low, medium, and high dose for 28 days as per the guidelines of OECD with pre-set time. They were under continuous observation for a behavioral sign, mortality rate, and body weight. Blood sampling followed by centrifugation at 1600 rpm for 20 min was done and freezed at -40 °C until further study. The liver function tests were performed in order to analyze the changes in enzyme activity and ensure safety.

Histopathology studies
The animals were administered ATV-NLCs and standard ATV-suspension formulation over a period of 28 days for toxicity study. Animals were sacrificed and liver was isolated according to the procedure. It was placed in a 10 % formalin solution at 4 °C temperature for further analysis.
The isolated tissues were washed with alcohol, xylene, and finally embedded in paraffin. These tissues were further sliced, dewaxed, and stained with hematoxylin and eosin for microscopic examination [55].

Screening of solid lipid, liquid lipid, and selection of surfactants
Among the solid lipids screened, the GMDG was selected as potential solid lipid based on their highest solubilizing capacity for ATV i.e. GMDG > Compritol ® ATO 888 > Precirol ® ATO5 > Gellucire ® 44/14 > Compritol ® HD 5 > Palmitic Acid > Myristic Acid. GMDG is a biocompatible and biodegradable lipid in combination with other ingredients. It was consistent with our previous publication that showed the highest solubility for ATV [56].
The liquid lipid, i.e. oil is the crucial excipient for NLCs production. The presence of oil can make the imperfect lattice arrangement to entrap the highest amount of drug. Among the tested oils, Capryol TM 90 was able to dissolve ATV to an appreciable extent and hence, selected for the preparation purpose [57].
The appropriate concentration of surfactant or surfactant blend reduces the interfacial tension of formulation and they are required to optimize the size of formulation. The surfactant's HLB values lies between 8 and 20 are able to form a stable formulation. Therefore, a surfactant blend of Tween 80 and Span 20 (2:1) were used to minimize the particle size and PDI and hence this blend was selected for development of NLCs [30].

Lipid mixture compatibility
Based on the evaluation of phase separation parameter, the ratio of 7:3 (solid lipid: liquid lipid) lipid mixture was preferred throughout the study. This ratio can be constructed the imperfect lattice arrangement of carrier which results in desired particle size and maximum drug entrapment.
QbD framework and its implementation 18

QTPP, CQA and risk assessment
In the QbD framework, initially the QTPP was assigned for the fabrication of NLCs achieving safe and efficacious formulation. The production of NLCs with planning from the initial stage was significant because a robust and reproducible production technique was required to meet the desired quality attributes. The identification of CQAs was linked with QTPP of expected formulation [33]. They were particle size (Y1), PDI (Y2), and zeta potential (Y3) and were presented as critical characteristics of formulation variables based on understanding and scientific knowledge [58]. For the desired quality product, the QTPP and CQA are summarized with their justifications in Table 3.
To perform the risk assessment, we plotted the Ishikawa diagram ( Fig. 1) displaying the various root causes and effects of formulation variables [36]. Employing the risk assessment matrix (RAM) method, the formulation variables had an impact on CQAs that were isolated. In this method, we highlighted the different risk levels (high for red, medium for yellow, and low for green) of material attributes and process parameters (Table 4) [20]. The variables having a low or negligible impact on responses were neglected. RAM method suggested the high-risk variables such as CMAs i.e. lipid concentrations, surfactant blend concentrations, and CPPs i.e. HPH pressure based on the significant influence on critical attributes.

Experimental design and statistical data analysis
In the present study, the CCRD was selected based on science and literature to observe the influence of formulation variables (lipid conc., surfactant conc., and HPH pressure) on their responses. As shown in Table 2, all the measurements were replicated three times to produce mean value (mean ± S.D., n = 3). The responses such as average particle size (Y1), PDI (Y2), and zeta potential (Y3) was ranged from 82.96 ±1.23 to 127.7 ± 1.82 nm, 0.116±1.93 to 0.486±1.04 and −7.1±0.21 to −13.7±1.05 respectively ( Table 2). As shown in Table 5, ANOVA suggested the fitted model was quadratic for responses and found to be significant (p<0.05) with an insignificant lack of fit p-value (p > 0.10) ( Table 6). The influence of formulation variables on the particle size and PDI of NLCs are illustrated in equations (5) and (6)

Influence of formulation variables on particle size and PDI
The particle size (Y1) and PDI (Y2) were the significant responses of formulation that determine the characteristics of NLCs. From equation 5, it was observed that a positive relationship with lipid concentration (A), since an increase in lipid concentration gradually increases the particle size and PDI. These are attributed to several experimental trials and well in agreement with our previous finding [56]. The lower concentration of lipid was significant to produce smaller particles and uniform distribution. On the other hand, larger particle size was the result of a higher concentration of lipid i.e. more than 3 % [38]. The surfactant concentration (B) was the main ingredient in the formulation to stabilize the NLC system. It was revealed that appropriate surfactant and their optimal concentrations lead to a decrease in mean particle size and PDI [59].
Moreover, owing to hydrophilic and lipophilic domains of Tween 80, they generate strong repulsion by adsorbing on particles during storage to stabilize the system [60]. The combination 20 of the surfactant (Tween 80) and co-surfactant (Span 20) was responsible for lowering particle size and PDI of the formulations. This was due to the reduction in interfacial tension and providing complete surface coverage to particles and provides steric repulsion.
The pressure of HPH significantly reduced the particle size during formulation processing. The generation of heat in homogenizer might be a contributing factor for the reduction of particle size due to the cavitational force produced in the homogenization gap. Beyond the optimal pressure, particle size increases gradually which could be due to possible collision and aggregation.
Besides, the generated heat affects the particles that could alter the entrapment efficiency [61,62].

Analysis of zeta potential (Y3)
Zeta potential (Y3) is a significant characteristic that develops on the particle surface offering stability through electrostatic repulsion [63]. In general, zeta potential of ± 30 is the desired value for good stability [64] but developed NLC showed a smaller zeta potential value as depicted in Table 2. This was due to the combined effect of non-ionic Tween 80 and Span 20 offering steric stabilization and avoiding the particle aggregation [65,66]. This indicates the particles were completely covered by the surfactants. Thus, the reduced zeta potential was an indicator of stabilized NLCs. These findings are well in agreement with other authors [25,67].

Percent entrapment efficiency (EE) and drug loading
EE (%) of ATV-NLCs was found to be from 85.5±0.24 to 96.23 ± 2.1. EE is essential for retention of drug molecules in lipid matrix. Higher solubility of the drug in lipid is the function of higher entrapment of drug. The content of oil in formulation offers an additional benefit to hold a high amount of drug in imperfect lattice arrangement. Thus, the combined effect of solid 21 lipid and oil prevents drug leakage and improves drug entrapment [56]. Moreover, drug loading capacity was found to be 1.92 %. The enhanced drug loading was due to the presence of components such as GMDG and Capryol 90 offering space to lodge the drugs due to imperfections in lipid matrix [68].

Formulation optimization and validation of experimental design
Design space was developed by obtaining a contour plot (Fig. 3A) and overlay plot (Fig. 3B) of formulation variables with given concentrations to meet the desired response. As can be seen in overlay plot (Fig. 3B), the colored region indicates all the desired criteria of responses. The working experimental optimal condition of NLCs development was selected based on highest desirability value. Using suggested formulation values, the ATV-NLC was freshly prepared and responses so obtained were compared with predicted values, and percent (%) bias was calculated as 2.05 for particle size, 16.66 for PDI, and 0.019 for zeta potential respectively (Table 7).

Lyophilization
Mannitol as a cryoprotectant offers a significant role in the lyophilization process as it covers the particles to prevent the aggregations. The optimal concentration of 3 % was found to have freeflowing characteristics of the powder [45]. The dry powder can be redispersed in distilled Milli-Q water for further analysis after short sonication [69].

Drug content
Drug content was calculated as 97.24 ± 0.42 % of optimized NLCs. The drug content indicated that the percentage of the total amount of actual ATV present in the formulation [70]. This was due to combining the solid lipid and liquid lipid domain that provides both mechanical stability and higher loading efficiency [71].

Fourier-transform infrared (FTIR) spectroscopy
We investigated the interactions among the drug and main components of NLCs using FTIR analysis presented in Fig. 4 (Fig. 4A) [72]. The spectra of GMDG (Fig. 4B), PM of GMDG, Capryol 90, and ATV

DSC study
Thermogram behavior of pure ATV, GMDG, PM, and optimized ATV-NLC are depicted in Fig.   5. Pure ATV shows a specific broad endotherm at 197 °C (Fig. 5A) which indicates the crystalline nature of drug [74]. The GMDG displayed the sharp endotherm at 60.93 °C (Fig. 5B) which is corresponding to their melting point. The disappearance of the ATV peak at its melting range in the PM due to the complete solubilization of drug in the lipid matrix (Fig. 5C). This was due to the conformational polymorphism during processing and indicates that ATV was molecularly dispersed in the lipid matrix that suggested an amorphous nature of drug present in PM [75].
Thermogram of ATV-NLC (Fig. 5D) shows the missing peak of ATV indicating the drug molecule was completely solubilized in lipid matrix and was in amorphous state. It is also noted that disappearance of peak associated with Kelvin effect, according to which particles in the nano range could melt at a temperature lower than the melting temperature of bulk material [53].
Importantly, as shown in Fig. 5, both exothermic and endothermic processes were occurred due to cold crystallization process by rapid cooling followed by crystalline melt in the heating pan.
Exothermic events can be observed due to the increase in heat flow in the sample. Moreover, the temperature experienced by the sample was higher than that of the reference pan [76].

Powder X-ray diffraction (PXRD)
The diffraction patterns of pure ATV, GMDG, PM, and optimized NLCs are depicted in Fig. 6.
The distinctive peaks of diffraction pattern were attributed to 2θ value in association with dvalue (Å). Pure ATV showed the diffraction peak at 2θ values of 3.102, 6.159, 9.23, and 12.307 during formulation processing and affected complete molecular solubility in lipid matrix [42,77,78].

TEM and FESEM photomicrograph
24 TEM (Fig. 7B) of a spherical shape with a smooth surface and FESEM (Fig. 7C) depicted the discrete nanoparticles with uniform distribution of NLC which is a close agreement with DLS results shown in Fig. 7A [46]. FESEM is an advanced imaging technique to evaluate the morphology of particles.

Stability studies
As can be seen in Table 8 the presence of cryoprotectant prevents aggregation by covering the particles significantly [79].

In vitro drug dissolution and mechanism of kinetics
The comparative drug release profiles from optimized ATV-NLC and ATV-suspension are shown in Fig. 8. The cumulative drug release 92.37 ±1.03 % and 36.13 ± 0.06 % was obtained from optimized ATV-NLC and ATV-suspension over a period of 24 h respectively (Fig. 8). The observed difference in the release pattern was due to the sustained release behavior of drugs from the encapsulated NLCs. It was found that the ATV-NLCs had an initial burst release of around 60-70 % at the end of 4 h, while the remaining amount of drug was released slowly over the period of 24 h [80]. The initial burst release from the optimized formulation for the first 4 h was attributed to the combined result of unentrapped ATV and adsorbed ATV on the surface of the nanocarriers [81]. It was also noted that the external portion of the nanocarrier was soft due to the presence of oil that contributed to the initial burst release [82]. The blend of lipids was not allowed to release the encapsulated drug particle equally due to the imperfect matrix structure. 25 Therefore, the subsequent controlled release of ATV from the lipid carrier was shown due to the diffusion of drug through the lipid matrix. In addition, the presence of Tween 80 (hydrophilic surfactant) was contributed to modulate drug release from the lipid matrix [38]. Thus, considering all evidence, it was revealed that NLCs could be the potent carriers for the sustained delivery of ATV. The Korsmeyer-Peppas release kinetics was found to be best fit model with the highest correlation coefficient (r 2 ) value of 0.925 ( Table 9). The release exponent (n) was found to be 0.63 (n=0.5-1.0) suggesting the non-fickian (anomalous) diffusion.

In vivo pharmacokinetic study and tissue distribution studies
Despite the conventional surgical procedure for evaluation of lymphatic transport, we used the method from the published article [30]. In this method, no surgical procedure was carried out for the establishment of pharmacokinetic profile. The common protein synthesis inhibitor i.e.
Cycloheximide was used for lymphatic flow blocking approach. It works by inhibiton of production of chylomicrons and inhibition of phagocytic activity of M cells thus block the pathway of lymphatic uptake of xenobiotics [1,30]. In addition, no interference of other absorption pathways and no death or other side effect was observed in cycloheximide treated rats during study. might further attempts to cause the first pass effects in liver [83]. 26 On the other hand, significantly higher (p<0.001) plasma concentration of Cmax and AUC0-t was observed in cohort-II receiving ATV-NLCs without treating cycloheximide. Because xenobiotics might be entered in lymphatic pathways avoiding reaching the liver for metabolism. This mechanism might be caused by the enhanced bioavailability of drug.
Herein, the key purpose of cycloheximide is that it only blocks the lymphatic uptake pathway and no interference of another absorption pathway. Additionally as stated earlier in the literature, the drug molecule should possess the log P value more than 5. Thus, in the present study, ATV has log p value 4.1 and have inability to reach in lymphatics by conventional mean [1].
Therefore, to achieve the main objective, we modulate novel lipid based NLC formulation to enhance the bioavailability. Several researchers were reported the lipid-based formulations that they have achieved the bioavailability using the chylomicron flow block model approach [30,[83][84][85].
From the results, ATV-suspension receiving cohort-I, it was observed that free drug suspension has not to access the lymphatic uptake, and hence the AUC level is much lower (Table 10). Thus, pharmacokinetic study showed that a significant (p<0.001) improvement in relative bioavailability of ATV in plasma 2.54 folds post-oral administration of ATV-NLCs when compared to ATV-suspension (  [87]. The concentration of ATV in suspension (Fig. 9B) and optimized NLCs (Fig. 9C) was estimated in lymphoid tissues such as thymus, spleen, and plasma ( This result was obtained due to the formulation components of the NLC that facilitated transport through the lymphatic pathway leads to rapid absorption into lymphatic organs [83,88]. The GMDG was assisted to stimulate the intestinal lymphatic absorption because the long-chain lipid is known as drug transporter in lymphatic as reported by Shete et al., 2013 [86]. The non-ionic surfactant, Tween 80, was known as a P-gp inhibitor that facilitates to stimulate the secretion of chylomicrons [84,89]. Nevertheless, Capryol 90 is effective in loosening tight junctions (TJs) and increased membrane fluidity of the intestinal epithelium that improves the lymphatic uptake [90]. As reported in Table 10, Tmax value for ATV-suspension and ATV-NLCs were found to be 2 h and 4 h respectively. At the same time, the higher MRT values of optimized ATV-NLCs indicated a longer residence time of the drug in the spleen as compared to plasma and thymus. This is well in agreement with previous findings reported in the article [23].

Oral dose toxicity and histopathology study
Animals were treated for 28 days with a repeated dose of control, and ATV-NLCs comprising of low, medium, and high dose. No significant variations in the weights of animals under study were found. However, two animals were died because of high dose and the remaining was survived throughout the study. Table 11 shows the results of the biochemical examination of 28 liver function tests like aspartate aminotransferase (AST), alanine aminotransferase (ALT), and alkaline phosphatase (ALP). Slight variations were observed in low and medium dose administered cohorts as compared to high dose administered cohorts. As can be seen in Table 11, the high dose cohorts showed an increased level of AST and ALT and ALP i.e. 2.3, 3.32, and 1.4 folds as compared to control cohorts. The low and medium dose values are normal in range for ATV-NLCs cohort when compared to the control. An increased level of AST and ALT was responsible for cell damage and the production of toxic species in the liver.
Histopathological examination was observed the normal size of vacuolation and sinusoidal space ( Fig. 10A) for control cohorts whereas suspension administered cohorts was found the liver injury that showing distortion in the liver with various degrees of sinusoidal and vacuolation spaces (Fig. 10B). In contrast, ATV-NLC administered cohort shows slight alterations (Fig. 10C) when compared to the control and suspension cohort. Thus, liver toxicity could be reduced by reduced particle size of nanocarriers and avoiding first-pass metabolism.

Conclusion
The implementation of a QbD concept throughout the product lifecycle implies the time and cost-saving factor to ensure quality products. In the present study, we successfully developed the ATV-NLCs using hot high pressure homogenization technique. The risk assessment tools such as the Ishikawa diagram and RAM method were implemented to identify the formulation variables. Optimized ATV-NLC were investigated physicochemically and revealed the existence of the amorphous state of ATV and its compatibility with lipid carriers. This was observed due to the complete molecular solubility and reduced particle of the formulation. In vitro study revealed the sustained release of ATV can be tailored to meet desired results. The enhanced bioavailability of ATV-NLCs was obtained to be 2.54 folds as compared to a pure drug suspension. The higher accessibility of ATV by the transcellular mechanism in lymphatics 29 avoiding the portal circulation. ATV-NLCs ensures lower or no toxicity associated with a low and medium dose of the drug that could be the promising carrier for the effective oral delivery of lipophilic drugs.

Compliance with ethical standards
Ethics approval and consent to participate All the animal experiments were performed according to the ARRIVE guidelines and U. K.
Animals (Scientific Procedures) Act, 1986 and associated guidelines.

Consent for publication
Not Applicable.

Availability of data and materials
Data available on request from the authors.

Declaration of interest
The authors declare that they have no competing interests.

Funding
Not Applicable.    Table 1 Different variables and their levels used in CCRD for production of ATV-NLCs.   Table 4 Risk assessment matrix (RAM) method for identifying experimental variables. Table 5 ANOVA for response measurements of particle size, PDI, and zeta potential.  Table 7 Predicted and observed responses of optimized formulations. Table 8 Stability study. Table 9 Release kinetic parameters. Table 10 Pharmacokinetic performance and tissue distribution study. Table 11 Liver function test at 0 th and 28 th day for observing effect of repeated administration of dose.