Int J App Pharm, Vol 10, Issue 2, 2018, 6269Original Article
APPLICATION OF RESPONSE SURFACE METHODOLOGY IN OPTIMIZATION OF PACLITAXEL LIPOSOMES PREPARED BY THIN FILM HYDRATION TECHNIQUE
AMOL A. TATODE^{1*}, ARUN T. PATIL^{2}, MILIND J. UMEKAR^{2}
^{1}Department of Pharmaceutics, University Department of Pharmaceutical Sciences, Rashtrasant Tukadoji Maharaj Nagpur University, Nagpur, Maharashtra 440033, India, ^{2}Smt Kishoritai Bhoyar College of Pharmacy, Kamptee, Nagpur, Maharashtra 441002, India
Email: aatatode@rediffmail.com
Received: 12 Dec 2017, Revised and Accepted: 23 Jan 2018
ABSTRACT
Objective: The present investigation was aimed to optimize the formula of paclitaxelloaded liposomes (PTL) by using the application of response surface methodology (RSM).
Methods: Paclitaxelloaded liposome (PTL) was optimized by response surface methodology based on two parameters, namely, percent entrapment efficiency (% EE) and percent in vitro drug release at 12 h (% DR). The liposome formula was prepared using 3^{2} factorial design, and the selected independent variables were, phospholipid (phospholipon 90G) and cholesterol (CH) concentrations. Nine formulas of paclitaxelloaded liposome were prepared by thin film hydration technique (THF). The entrapment efficiency, in vitro release studies and drug content, were evaluated using on UVvisible spectrophotometer at λmax230 nm. The developed PTL formulation vesicle morphology, particle size, polydispersity index (PDI) and zeta potential (ζ) were evaluated by Motic digital microscope and Malvern zetasizer respectively.
Results: Using response surface methodology the estimated coefficient values obtained for independent variables in the regression equations, exhibited that the phospholipid (PL90G) and cholesterol (CH) molar concentration was observed to be highly influencing variables in optimizing % EE (86.67±0.67) and % DR (63.49±1.21). In the prediction of % EE and % DR values, the percent relative errors (PRE) was found to be low (–0.290%) and (0.058%) respectively. This suggests that designdeveloped model was found to be suitable for PTL formulations and thus, validate the model.
Conclusion: Experimental results show that the observed responses were in close agreement with the predicted values and this demonstrates the reliability of the RSM in an optimization of % EE and % DR in paclitaxel liposomal (PTL) formulations.
Keywords: Response surface methodology, Paclitaxel, Liposomes, Thin film hydration technique
© 2017 The Authors. Published by Innovare Academic Sciences Pvt Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/)
DOI: http://dx.doi.org/10.22159/ijap.2018v10i2.24238
INTRODUCTION
Paclitaxel (PT) is a chemical compound isolated from the bark of Taxus brevifolia (northwest Pacific Yew Tree), empirical formula (C_{47}H_{51}NO_{14}) and on the basis of characterization named it as Taxol [1]. It has the potential anticancer drug, based on previous reports, PT shows anticancer activities towards breast cancer [2], ovarian cancer [3], lung cancer [4] and pancreatic cancer [5]. However, PT exhibits poor aqueous solubility and permeability owing to biopharmaceutical classification system (BCS) class IV drug, which directs it to low bioavailability. Therefore, by consideration of these problems, there is a need to develop a novel formulation of such effective and efficient anticancer drug.
For improving poor aqueous solubility and permeability of PT, a number of formulation strategies have been developed and used. Some of them were modified, due to some excipient–drug interactions. For improving the solubility, PT dissolved in a mixture of polyoxyethylated castor oil (Cremophore EL): dehydrated ethanol (1:1) ratio as a delivery vehicle. The formulation produced hypersensitivity and nonlinear pharmacokinetic behavior after intravenous administration. The hypersensitivity reaction at the site of administration could be due to an inclusion of Cremophore EL [67]. After that, the delivery vehicle was replaced with the addition of tween 80 alone or combination of tween 80: dehydrated alcohol, and diluted with aqueous media. The diluted formulation showed the precipitation of PT from solution due to low solubility [8]. These attempted techniques, with persistent low solubility problem, has been overcome by creating novel formulation with the aim of improving aqueous solubility, permeability, and bioavailability of PT. It includes novel oral formulation [9], novel PT selfemulsifying drug delivery system (SEDDS) [10], novel ligands based PT targeting formulation [11], micellar formulation [12], liposomal formulation [13], bioconjugates [14], dendrimers [15] and nanocarrier systems [16]. In all these formulation techniques the problem associated with PT was shown to be improved significantly.
The liposome is emerging techniques for specialized drug delivery [17] and best suitable for lipophilic drug due to its biocompatibility and reducing drug toxicity, with maintaining efficacy of the anticancer drug for a maximum period of time. Some previous studies include asulacrine [18], docetaxel [19] and tamoxifen [20] with these approaches, their poor aqueous solubility and bioavailability were found to be improved. So, need to develop and optimize the paclitaxelloaded liposomes (PTL) for effective anticancer treatment. In pharmaceutical technology, in the development and optimization of different pharmaceutical dosage forms, there are a high number of factors which influence the product characteristics. Therefore, complex, expensive and timeconsuming formulation studies are often necessary for the development of a product with required and desired properties. The experimental design methodology is a strategy to use a smaller number of experiments and to avoid unnecessary experiments [2122].
Experiments were designed to determine the effect of the independent variables (factor) on the dependent variable (parameter/response) of a process or formulation. RSM, one of the designs of experiments, is a powerful tool for determining the relationship between a response and a set of quantitative involved factors. RSM is a technique used to find the optimum response by using the quadratic polynomial model [23]. The advantage of RSM is the reduced amount of experiments required, thereby reducing the cost of expensive analysis methods. The application of RSM is useful for understanding or mapping a region of the response surface, finding the variable level of optimum response, and selecting the process condition or formula to meet the specifications [24]. This research was carried out to optimize PTL formula with independent variables such as phospholipid (phospholipon 90G) concentration and cholesterol (CH) concentration. The optimum formula was obtained from RSM using 3^{2} full factorial design. The optimization approach was applied to obtain desired % EE and % DR for PTL.
MATERIALS AND METHODS
Materials
Paclitaxel (PT), (purity>90%) was received as a gift sample from MACCHEM Products (India) Pvt. Ltd. Bhoisar, Thane, India. The phospholipids samples viz., Phospholipon 90G^{®} (PL90G), Phospholipon 80H^{®} (PL80H) and Phospholipon 90H^{®}(PL90H) with purity >90%, was obtained as a free gift sample from Lipoid GmbH, Ludwigshafen, Germany. The solvents namely chloroform and methanol were purchased from Merck Ltd. Mumbai, India. Cholesterol (CH), potassium dihydrogen phosphate and sodium hydroxide pellets were obtained from Sigma Chemicals, SigmaAldrich Corporation, St. Louis, MO. Chemical used in this work were of analytical grade (AR).
Experimental design (3^{2} full factorial design)
To reduce the number of trials and attain the highest amount of information on product properties, the screening was done by applying full factorial design (3^{2}), systematically study the joint influence of the independent variables on the dependent variables. So, in this study two factors were evaluated, each at three levels and experimental trials was performed at all nine possible combinations. Amount of phospholipid (PL90G1, 2 and 3 moles) was taken as the first independent variable (X_{1}, w) and amount of cholesterol (CH1, 2 and 3 moles) was selected as the second independent variables (X_{2}, w) for liposomes. These variables varied at three levels, low level (1), medium level (0), and high level (+1). All the calculations were done at milligram level. Amount of PT (10 µM) and final formulation volume 15 ml was kept constant. Percent entrapment efficiency (% EE) (Y_{1}) and percent in vitro drug release at 12 h (% DR) (Y_{2}) were selected as dependent variables. Values of variables and batch codes are shown in table 1 and 2. Design Expert^{®} DX 10.0.7.0 (StatEase Inc., MN) license version software was used for the generation and evaluation of statistical experimental design [2526].
Preparation of paclitaxelloaded liposomes (PTL)
Liposomes were prepared by the thin film hydration method (TFH). The 10 µM PT (mol. wt., 853.9) constant for all batches and the required quantities of phospholipid (PL90G) (mol. wt., 758.07) and cholesterol (CH) (mol. wt., 386.67) were taken in a 100 ml round bottom flask and dissolved in 10 ml chloroform. All the batches were prepared according to the experimental design in table 1. Chloroform was evaporated using rotary vacuum evaporator (Model: PBV–7D, Vertical condenser, rotavap, superfit^{TM} continental Pvt. Ltd., Mumbai, India) and kept overnight under vacuum. Then it was hydrated by 15 ml of phosphate buffer pH 7.4 for 1 h with 10 min of extensive vortexing. The suspension of liposomes was sonicated in the water bath at 60 °C to reduce the size of liposomes. Nonincorporated PT was separated by ultracentrifuge at 10,000 rpm for 30 min at 4 °C. The supernatant was then discarded and PT loaded liposomes in the precipitate were redispersed in required volume of phosphate buffer pH 7.4. This was transferred to vials and stored at 4 °C [27].
Table 1: 3^{2}full factorial design: factors, factor levels and responses for PTL formulation
Factors (Independent variables)  Factor levels used 

Low (1)  Medium (0)  High (+1)  
Amount (moles) of phospholipid (PL90G) (X_{1}, w)  1  2  3 
Amount (moles) of cholesterol (CH) (X_{2}, w)  1  2  3 
Responses (Dependent variable)  
Y_{1} = Percent entrapment efficiency (% EE)  
Y_{2} = Percent in vitro drug release (% DR) 
Statistical analysis and optimization of formulation using RSM
Response surface modeling and evaluation of the quality of fit of the model for the current study were performed employing Design Expert^{®} DX 10.0.7.0 license version software [2326, 30]. Polynomial models including linear, interaction and quadratic terms were generated for all the response variables using multiple linear regression analysis (MLRA). A secondorder polynomial equation that describes the effect of independent factors on the response is expressed in the following forms:
Linear model: Y = β _{0}+b_{1}X_{1}+b_{2}X_{2} (1)
2FI (interaction model): Y = β _{0}+β _{1}X_{1}+β _{2}X_{2}+β _{12}X_{1}X_{2}(2)
Quadratic model = β _{0}+β _{1}X_{1}+β _{2}X_{2}+β _{12}X_{1}X_{2}+β _{11}X_{1}^{2}+β _{22}X_{2}^{2}(3)
Where Y is the dependent variable; β _{0} is the arithmetic mean response of the nine runs and β_{i} (β_{1}; β_{2}; β_{12}; β_{11} and β_{22}) is the estimated coefficient for the corresponding factor X_{i} (X_{1}, X_{2}, X_{1}X_{2}, X_{1}X_{1}, and X_{2}X_{2}). The main effects (X_{1} and X_{2}) represent the average result of changing one factor at a time from its low to high value. The interaction terms (X_{1}X_{2}) show how the response changes when two factors are simultaneously changed. The polynomial terms (X_{1}^{2} and X_{2}^{2}) are included to investigate nonlinearity. The equations enable the study of the effects of each factor and their interaction over the considered responses. The polynomial equations were used to draw conclusions after considering the magnitude of coefficients and the mathematical sign they carry, i.e. positive or negative. A positive sign signifies a synergistic effect, whereas a negative sign stands for an antagonistic effect. The best fitting mathematical model was selected based on the comparisons of several statistical parameters, including the coefficient of variation (CV), the coefficient of determination (R^{2}), adjusted coefficient of determination (Adjusted R^{2}) and the predicted residual sum of square (PRESS), provided by Design Expert software. Among them, PRESS indicates how well the model fits the data and for the chosen model it should be small relative to the other models under consideration. Level of significance was considered at p<0.05. Mathematical relationships in the form of polynomial equations are generated using multiple linear regression analysis (MLRA) and used to find out the relative influence of each factor on the response. Analysis of variance (ANOVA) for the responses was performed to identify a significant effect of factors on responses and the model parameters were obtained. The relationship between the dependent and independent variables was further elucidated using contour and response surface plots. Twodimensional contour plots and threedimensional response surface plots resulting from equations were obtained by the Design Expert software. These plots are very useful in a study of the effects of two factors on the response at one time and predict the responses of dependent variables at the intermediate levels of independent variables. Subsequently, a numerical optimization technique by the desirability and graphical optimization technique by the overlay plot approach were used to generate the new formulation with the desired responses. An optimized formulation was developed by setting constraints (goals) on the dependent and independent variables. To validate the chosen experimental design, the resultant experimental values of the responses were quantitatively compared with those of the predicted values and calculated the percent relative error (PRE) by the following equation 5.
% Relative error = … (4)
Determination of percent entrapment efficiency (% EE)
Purification of PTL formulation was done by the ultracentrifugation method [28]. To quantify the amount of entrapped PT, 2 ml of the vesicular dispersion was centrifuged at 10,000 rpm for 1 h at the controlled temperature of 4 °C (Remi cooling centrifuge, Remi Elektrotechnik limited, India). Supernatant contains unentrapped drug was withdrawn and measured UV spectrophotometrically (at λmax230 nm) (Model: SPECTRO 2060 PLUS, Analytical Technologies Ltd., Gujarat, India) against 30:70 ratio of methanol: phosphate buffer solution (PBS) (pH 7.4). All the determinations were made in triplicate. A calibration plot was produced by diluting stock solutions of PT with 30:70 ratio of methanol and PBS (pH 7.4). % EE was calculated and expressed as a percent of the available dissolved solute actually encapsulated. The amount of drug entrapped in liposomes was determined by equation 1.
% entrapment efficiency (% EE) = (5)
Percent in vitro drug release study (% DR) at 12 h.
The in vitro drug release study for PT from different PTL formulation was carried out as per the procedure described by Utreja [29] with little modifications. In brief, the Franz diffusion cell apparatus was employed for this study. The apparatus is consisted of donor and receptor compartment, with an effective surface area for dissolution was (2.303 cm^{2}). The dialysis membrane (LA395, Dialysis Membrane–110 AV, flat width ~ 31.12 mm, Average diameter ~ 21.5 mm, and approximate capacity is ~ 3.63 ml/cm; HI media laboratories, Mumbai, India) was employed and pretreated as per the directions were given by the manufacturer. After proper pretreatment, the membrane was cut into desired size and shape, then mounted between the effective surface area of donor and receptor compartment. The PTL dispersion (2 ml) was placed over the membrane, accompanied by an addition of PBS (20 ml, pH 7.4) contain 0.1% tween 80 as dissolution media in the receptor compartment. The contents of receptor compartment were stirred at 100 rpm using magnetic stirrer at 37±1.0 °C. At specified time intervals, 2 ml aliquots were withdrawn from sampling port of apparatus, diluted suitably with fresh media and the absorbance of the resulting solution was read at 230 nm using UVvisible spectrophotometer (Model: SPECTRO 2060 PLUS, Analytical Technologies Ltd., Gujarat, India).
Vesicle morphology study of liposomes
The liposomes were mounted on glass slides and viewed under a Motic Digital Microscope (type DM1802) for morphological observation after suitable dilution. The size analysis of PTL was examined at the magnification of (×40) using calibrated eyepiece micrometer. The images were recorded using Motic Image plus 2.0 ML software, accompanying with the instrument.
Determination of percent drug content
One milliliter of dispersion was pipette from the PTL formulation and lysed with methanol. It was further diluted with 30:70 ratio of methanol: phosphate buffer solution (PBS) (pH 7.4) and sample were analyzed spectrophotometrically at λmax 230 nm for PT. [26]
Determination of particle size, polydispersity index (PDI) and zeta potential (ζ)
The size of liposomes was measured by dynamic light scattering with a Malvern zetasizer. Diluted (1:100) PTL dispersion was added to the sample cuvette and then cuvette is placed in zetasizer. The sample is stabilized for two minutes and reading was measured. The average particle size was measured after performing the experiment in triplicate. The zeta potential of developed PTL formulation was determined using Malvern zetasizer (Malvern zetasizer ver. 6.20, UK). The zeta potential was calculated by HelmholtzSmoluchowski's equation from the electrophoretic mobility of liposomes at 25 °C [26].
RESULTS AND DISCUSSION
Experimental design and data acquiring (3^{2}full factorial design)
Full factorial design (3^{2}) was applied to optimize the PTL formulation. All nine batches of PTL were prepared according to the formulation variables as shown in table 2. Liposomes were obtained by the TFH method. RSM was exploited to estimate the influence of the molar ratio of PL90G and CH as independent variables and their interactions on the investigated responses (dependent variables; % EE and % DR). This experiment was aimed to identify considerable factor effect influencing the formulation performance and to set up to their excellent levels for the desirability of responses shown in table 2.
Statistical analysis and optimization of formulation using RSM
To evaluate the quantitative effects of factors (X_{1} and X_{2}) and their levels low (1), middle (0), and high (+1) on the preferred responses, the experimental values of the flux were analyzed by Design Expert^{®} DX 10.0.7.0 license version software and mathematical models obtained for each response [2526, 3031]. The mathematical relationship generated using multiple linear regression analysis (MLRA) for the studied response variables (% EE and % DR at 12 h.) that were relating different response and independent variables are expressed as following polynomial equations (quadratic model).
Y_{1} (% EE) = 76.63+3.10X_{1–}9.44X_{2}+3.18 X_{1}X_{2–}9.91X_{1}^{2}+0.35X_{2}^{2} (6)
Y_{2} (% DR12 h.) = 64.74+6.79X_{1}–4.63X_{2}–1.38X_{1}X_{2}–15.29X_{1}^{2}–5.82X_{2}^{2} (7)
The above equations expose the quantifiable effect of the independent variables, a molar ratio of PL90G and CH, on the responses such as % EE (Y_{1}) and in vitro % DR at 12 h (Y_{2}) as dependent variables. The fitted polynomial equation (quadratic model) related to % EE and percent in vitro % DR used to draw a conclusion after considering the coefficient and the mathematical sign it carries. i.e. positive and negative. The correlation coefficient (r^{2}) of the quadratic model (0.9736) for response Y_{1} (% EE) and (0.9779) for response Y_{2} (% DR) was found to be significant.
Response 1 (Percent entrapment efficiency) (% EE)
Regression analysis of above equation (6) of response Y_{1} (% EE) revealed that the coefficient of β_{1} was positive and β_{2} was negative, this indicated that as PL90G (X_{1}) increased the % EE increased but as we further increased the PL90G (X_{1}) to higher level the % EE decreased and on increasing cholesterol (X_{2}) the % EE decreased. The higher concentration of cholesterol leads to rigidity in the vesicles [26] which in turn decreased the % EE. The % EE of different liposomal batches was in a range of 51.68 to 86.67%. The maximum entrapment was observed in batch L4 (table 2) with the composition of PL90G: CH (2:1 molar ratio) (0,1).
Table 2: Composition 3^{2} full factorial design with measured responses of PTL formulation
Batches  Variable level in coaded form  Variable level in actual form  Response variables  
X_{1}  X_{2}  Phospholipid (PL90G) in moles (X_{1}, W)  Cholesterol (CH) in moles (X_{2}, W) 
Percent entrapment efficiency* (% EE)±SD  Percent In vitro drug release* (12 h) (% DR)±SD  
L1  1  1  1  1  75.65±1.48  40.28±1.58 
L2  1  0  1  2  64.23±0.82  43.68±0.98 
L3  1  +1  1  3  51.68±0.75  32.39±1.02 
L4  0  1  2  1  86.67±0.67  63.49±1.21 
L5  0  0  2  2  78.69±1.61  62.12±0.63 
L6  0  +1  2  3  65.23±1.17  56.98±1.33 
L7  +1  1  3  1  77.12±1.53  56.32±1.42 
L8  +1  0  3  2  67.16±0.91  57.85±0.87 
L9  +1  +1  3  3  65.89±1.32  42.92±1.19 
*Values represented as mean±SD, n = 3, All baches contain drug 10 µM and 15 ml phosphate buffer (pH 7.4) for hydration.
Table 3: Analysis of variance (ANOVA) table of % EE
Source  Sum of squares  df  Mean squares  F Value  pvalue Prob>F  
Model  829.57  5  165.91  22.14  0.0142  Significant 
X_{1}–Phospholipid (PL90G)  57.72  1  57.72  7.70  0.0692  
X_{2}–Cholesterol (CH)  534.68  1  534.68  71.36  0.0035  
X_{1}X_{2}  40.58  1  40.58  5.42  0.1024  
X_{1}^{2}  196.35  1  196.35  26.20  0.0144  
X_{2}^{2}  0.24  1  0.24  0.032  0.8693  
Residual  22.48  3  7.49  
Cor–total  852.05  8 
For estimation of the significance of the model, the analysis of variance (ANOVA) was executed, from the ANOVA data, the model Fvalue of response (Y_{1}) (22.14) indicated that the model is significant shown in table 3. There is only a 1.42% chance that an Fvalue this large could occur due to noise. Values of “prob>F” less than 0.0500 indicate model terms are significant. In these case X_{2}, X_{1}^{2} are significant model terms. Values greater than 0.1000 indicate that model terms are not significant.
Table 4: Parameter of selected quadratic model of % EE
Std. dev.  2.74  Rsquared (r^{2})  0.9736 
Mean  70.26  Adjusted RSquared  0.9296 
C. V. %  3.90  Predicted RSquared  0.7315 
PRESS  228.77  Adequate Precision  15.693 
The predicted Rsquared value of 0.7315 is in reasonable agreement with the adjusted Rsquared of 0.9296; i.e. the difference is less than 0.2. Adequate precision measures the signal to noise ratio. A ratio greater than 4 is desirable and the result of adequate precision was 15.693 indicates an adequate signal. So, this model can use to navigate the design space.
The relationship between the dependent and independent variables was further elucidated using contour and response surface plots. The contour (fig. 1) and 3D response surface plots (fig. 2) of % EE clearly indicated that X_{1} and X_{2} highly influenced the response 1 (% EE). The change in % EE as a function of X_{1} and X_{2} was depicted in the form of contour and response surface plots based on full factorial design (3^{2}). So, middle level of X_{1} and low level of X_{2} was found to be favorable conditions for obtaining higher % EE.
Fig. 1: Counter plot showing the effect of phospholipid (PL90G) (X_{1}) and cholesterol (CH) (X_{2}) on % EE (Y_{1}) of PTL
Response 2 (Percent drug release at 12 h) (% DR)
The effect on drug release at 12 h (% DR) (Y_{2}) was observed to be significant (P<0.05) by ANOVA and the polynomial equation (7) revealed that the coefficient of β_{1} was positive and β_{2} was negative, this indicated that as PL90G (X_{1}) increased the % DR increased and on increasing cholesterol (X_{2}) the % DR decreased. The % DR increased with increased concentration of lipid and at a certain level the percent release is retarded above that and the release was decreased at higher levels of cholesterol. This is because cholesterol at higher levels makes the lipid bilayers more rigid and retards the release of the drug. This was evident by the higher cholesterol concentration of vesicles showed around 50 % of the release except for (L6) formulations. The L4 formulation found to have 63.43 % DR at 12 h (table 2) with the composition of PL90G: CH (2:1 molar ratio) (0,1). At lower concentration of phospholipid and cholesterol, the drug release was very less due to the formation of stagnant layer [26].
Fig. 2: Response surface plot showing the effect of phospholipid (PL90G) (X_{1}) and cholesterol (CH) (X_{2}) on % EE (Y_{1}) of PTL
Table 5: Analysis of variance (ANOVA) table of % DR at 12 h
Source  Sum of squares  df  Mean squares  F Value  pvalue prob>F  
Model  948.33  5  189.67  26.53  0.0109  significant 
APhospholipid  276.62  1  276.62  38.70  0.0084  
BCholesterol  128.81  1  128.81  18.02  0.0239  
AB  7.59  1  7.59  1.06  0.3786  
A^{2}  467.57  1  467.57  65.41  0.0040  
B^{2}  67.74  1  67.74  9.48  0.0542  
Residual  21.44  3  7.15  
Cor Total  969.78  8 
From the ANOVA data, the model Fvalue of response (Y_{2}) (26.53) indicated that the model is significantly shown in table 5. There is only a 1.09% chance that a model Fvalue this large could occur due to noise. Values of “prob>F" less than 0.0500 indicate model terms are significant. In these case, X_{1}, X_{2}, and X_{1}^{2} are significant model terms. Values greater than 0.1000 indicate that model terms are not significant.
Table 6: Parameter of selected quadratic model of % DR at 12h
Std. Dev.  2.67  Rsquared (r^{2})  0.9779 
Mean  50.67  Adjusted RSquared  0.9410 
C. V. %  5.28  Predicted RSquared  0.8063 
PRESS  187.85  Adequate Precision  14.272 
The predicted Rsquared value of 0.8063 is in reasonable agreement with the adjusted Rsquared of 0.9410; i.e. the difference is less than 0.2. Adequate precision measures the signal to noise ratio. A ratio greater than 4 is desirable and result of adequate precision was 14.272 indicates an adequate signal. So, this model can be used to navigate the design space.
Fig. 3: Counter plot showing the effect of phospholipid (PL90G) (X_{1}) and cholesterol (CH) (X_{2}) on % DR at 12 h (Y_{2}) of PTL
Fig. 4: Response surface plot showing the effect of phospholipid (PL90G) (X_{1}) and cholesterol (CH) (X_{2}) on % DR at 12 h (Y_{2}) of PTL
To envisage the effect of an independent factor on the response (Y_{2}) the contour plot (fig. 3) and 3Dresponse surface plots (fig. 4) of % DR at 12 h shows the curvature with a change in the factor (X_{1} and X_{2}). The plot was found to be curvilinear and indicated that a high value of % DR (63.43%) can be obtained for a combination middle level of X_{1} and low level of X_{2} factors.
Desirability and overlay plot
The aim of pharmaceutical formulation optimization is generally to find the levels of the variable that affect the chosen responses and determine the levels of the variable from which a robust product with highquality characteristics may be produced. All the measured responses that may affect the quality of the product were taken into consideration during the optimization procedure. The % EE and % DR at 12 h were set out in the maximum criteria. Each response criterion was combined (overlay plot) to obtain the optimum value (fig. 5). The optimization results of this research can be seen in table 7.
Validation of RSM results
In order to evaluate the optimization capability of models generated according to the results of the RSM (3^{2}full factorial design), PTL formulation was prepared using the optimal process variables settings that X_{1} and X_{2} were equal to 2:1. The response Y_{1} (% EE) and Y_{2} (% DR at 12 h) obtained with predicted models and the experimental model were shown in table 8. The percent relative error was obtained using equation 4. The percent relative error (PRE) for response Y_{1} (% EE) and Y_{2} (% DR at 12 h) were found to be (–0.290) and (0.058) respectively. The maximum PRE value was (–0.290). However, the values were found to be<2 % and hence it confirmed the suitability of experimental design. The results showed good agreement on preparation properties with theoretical properties.
Table 7: Characteristics of optimum formula
Objects  Phospholipid (PL90G) in moles (X_{1}, W) 
Cholesterol (CH) in moles (X_{2}, W) 
% EE (Y_{1}, %) 
In vitro % DR at 12 h (Y_{2}, %) 
Desirability  
Predicted  1.996  1.000  86.419  63.527  0.996  Selected 
Actual (L4)  2  1  86.67±0.67  63.49±1.21 
The optimization parameter of desirability was determined by regulating the optimum input variables to obtain one or more optimal parameters. The desirability value ranged between 0 and 1, where a value of 1 is perfect, i.e., the ideal parameter value [30]. The PTL desirability plot was shown in fig. 6. The optimizing desirability of PTL formulation was 0.996. This value was near to ideal value (1), meaning that the predicted parameters were desired parameter values. The composition of the predicted formulations was matching with L4 liposomes (table 7).
Fig. 5: Overlay plot of % EE and % DR at 12 h
Fig. 6: PTL desirability plot
Table 8: Validation of predicted and experimental PTL batch
Response  Experimental values  Predicted value  % Relative error (PRE) 
Y_{1} (% EE)  86.67±0.67  86.419  –0.290 
Y_{2} (% DR at 12 h)  63.49±1.21  63.527  0.058 
Vesicle morphology, percent drug content, particle size, PDI and zeta potential (ζ) of developed PTL formulation
Vesicle morphology of developed PTL formulation was observed by Motic Digital Microscope (type DM1802). The liposomes were spherical in shape with a smooth surface shown in fig. 7.
The developed liposomal percent drug content was found to be 98±1.0 %. (mean±SD, n = 3). Percent drug content indicated that the PT was uniformly distributed in vesicular dispersions and percent drug content near to 100 % indicated no loss of the material during the preparation.
Fig. 7: Microphotograph of PTL by motic image plus 2.0 ML software
Fig. 8: The particle size (A) and Zeta potential (B) of developed PTL formulation
The particle size (fig. 8A) and (ζ) zeta potential (fig. 8B) of developed PTL formulation was found to be 144.4 nm and () 22.6 mV. The polydispersity index was calculated as 0.224. The low polydispersity index indicates a narrow range of particle size distribution. The zeta potential was a reliable indicator in the prediction of stability of particles in a liquid medium and the possible interactions with other materials. The behavior (size and size distribution) of vesicles completely depended on the amount of selected variables were also reported [32].
CONCLUSION
In the present study, the response surface methodology i.e.3^{2} full factorial design was successfully employed for the optimization of PTL formulations. The PTL was prepared by thin film hydration method. The results of above optimization study displayed that the phospholipids (PL90G) and cholesterol (CH) with a molar ratio (2:1) showed an enhancement in rate and extent of in vitro release of PT from designoptimized PTL formulations. Thus, we conclude that the proposed RSM could be useful for the preparation and optimization of paclitaxelbased liposomal formulations.
ACKNOWLEDGEMENT
The authors wish to express their thanks and gratitude to University Dept. of Pharmaceutical Sciences, RTM Nagpur University, Nagpur and SKB College of Pharmacy, Kamptee, Nagpur for providing all necessary facilities for completion of this project work.
AUTHORS CONTRIBUTION
This work was carried out in collaboration between all the four authors in the concept and design of the work, collection, assembly, analysis and interpretation of data, writing, critical revision and approval of the final manuscript.
CONFLICT OF INTERESTS
All authors have none to declare
REFERENCES
Louage B, De Wever O, Hennink WE, De Geest BG. Developments and future clinical outlook of taxane nanomedicines. J Controlled Release 2017;253:13752.
Ye J, Liu Y, Xia X, Meng L, Dong W, Wang R, et al. Improved safety and efficacy of a lipid emulsion loaded with a paclitaxelcholesterol complex for the treatment of breast tumors. Oncol Rep 2016;36:399409.
Kampan NC, Madondo MT, McNally OM, Quinn M, Plebanski M. Paclitaxel and its evolving role in the management of ovarian cancer. BioMed Res Int 2015:121. http://dx.doi.org/10.1155/2015/413076
Wang L, Yu RS, Yang WL, Luan SJ, Qin BK, Pang XB, et al. Effects of paclitaxel loadeddrug micelles on cell proliferation and apoptosis of human lung cancer A549 cells. Acta Pharma Sin 2015;50:12405.
Lemstrova R, Melichar B, MohelnikovaDuchonova B. Therapeutic potential of taxanes in the treatment of metastatic pancreatic cancer. Cancer Chemother Pharmacol 2016;78:110111.
Webster LK, Woodcock DM, Rischin D, Millward MJ. Review: cremophor: the pharmacological activity of an" inert" solubiliser. J Oncol Pharm Pract 1997;3:18692.
Cragun JM, Baggs JH, Rollins C, Chambers SK. Case report hypersensitivity reaction to parenteral nutrition after severe hypersensitivity reaction to paclitaxel. Am J Clin Exp Obstet Gynecol 2013;1:6975.
Singla AK, Garg A, Aggarwal D. Paclitaxel and its formulations. Int J Pharm 2002;235:17992.
Moes J, Koolen S, Huitema A, Schellens J, Beijnen J, Nuijen B. Development of an oral solid dispersion formulation for use in lowdose metronomic chemotherapy of paclitaxel. Eur J Pharm Biopharm 2013;83:8794.
Sandhu PS, Beg S, Mehta F, Singh B, Trivedi P. Novel dietary lipidbased selfnano emulsifying drug delivery systems of paclitaxel with pgp inhibitor: implications on cytotoxicity and biopharmaceutical performance. Expert Opin Drug Delivery 2015;12:180922.
Li J, Wang F, Sun D, Wang R. A review of the ligands and related targeting strategies for active targeting of paclitaxel to tumours. J Drug Target 2016;24:590602.
Xu J, Zhang X, Chen Y, Huang Y, Wang P, Wei Y, et al. Improved micellar formulation for enhanced delivery for paclitaxel. Mol Pharm 2016;14:3141.
Bonde S, Nair S. Advances in liposomal drug delivery system: fascinating types and potential applications. Int J Appl Pharm 2017;9:17.
Li J, Huang P, Chang L, Long X, Dong A, Liu J, et al. Tumor targeting and pHresponsive polyelectrolyte complex nanoparticles based on hyaluronic acidpaclitaxel conjugates and chitosan for oral delivery of paclitaxel. Macromol Res 2013;21:13317.
Teow HM, Zhou Z, Najlah M, Yusof SR, Abbott NJ, D’Emanuele A. Delivery of paclitaxel across cellular barriers using a dendrimerbased nanocarrier. Int J Pharm 2013;441:70111.
AlNajjar BY, Hussain SA. Chitosan microspheres for the delivery of chemotherapeutic agents: paclitaxel as a model. Asian J Pharm Clin Res 2017;10:15.
Kulkarni PR, Yadav JD, Vaidya KA. Liposomes: a novel drug delivery system. Int J Curr Pharm Res 2011;3:108.
Zhang W, Wang G, Falconer JR, Baguley BC, Shaw JP, Liu J, et al. Strategies to maximize liposomal drug loading for a poorly watersoluble anticancer drug. Pharm Res 2015;32:145161.
Zhang H, Gong W, Wang ZY, Yuan SJ, Xie XY, Yang YF, et al. Preparation, characterization, and pharmacodynamics of thermosensitive liposomes containing docetaxel. J Pharm Sci 2014;103:217783.
Bhatia A, Singh B, Raza K, Shukla A, Amarji B, Katare OP. Tamoxifenloaded novel liposomal formulations: evaluation of anticancer activity on DMBATPA induced mouse skin carcinogenesis. J Drug Target 2012;20:54450.
Carbone C, Tomasello B, Ruozi B, Renis M, Puglisi G. Preparation and optimization of PIT solid lipid nanoparticles via statistical factorial design. Eur J Med Chem 2010;49:1107.
Hermans K, Van den Plas D, Everaert A, Weyenberg W, Ludwig A. Full factorial design, physicochemical characterisation and biological assessment of cyclosporine loaded cationic nanoparticles. Eur J Pharm Biopharm 2012;82:2735.
Krishnaiah D, Bono A, Sarbatly R, Nithyanandam R, Anisuzzaman SM. Optimisation of spray drying operating conditions of Morinda citrifolia L. Fruit extracts using response surface methodology. J King Saud Univ 2012;27:2636.
Minitab. What is a Response Surface Design; 2016. Available from: http://www.support.minitab.com/enus/minitab/17/ topiclibrary/modelingstatistics/doe/responsesurfacedesigns/whatisaresponsesurfacedesign. [Last accessed on 10 Nov 2017].
Ghanbarzadeh S, Valizadeh H, Zakeri Milani P. Application of response surface methodology in the development of sirolimus liposomes prepared by thin film hydration technique. BioImpacts 2013;3:7581.
Sudhakar B, Krishna MC, Murthy KV. Factorial design studies of antiretroviral drugloaded stealth liposomal injectable: PEGylation, lyophilization and pharmacokinetic studies. Appl Nanosci 2016;6:4360.
Dua JS, Rana AC, Bhandari AK. Liposome: methods of preparation and applications. Int J Pharm Stud Res 2012;3:1420.
Yang T, Cui FD, Choi MK, Lin H, Chung SJ, Shim CK, et al. Liposome formulation of paclitaxel with enhanced solubility and stability. Drug Delivery 2007;14:3018.
Utreja P, Jain S, Tiwary AK. Localized delivery of paclitaxel using elastic liposomes: Formulation development and evaluation. Drug Delivery 2011;18:36776.
Damai RS, Silvia S, Etik M. Optimization of luteolinloaded transfersome using response surface methodology. Int J Appl Pharm 2017;9(Suppl):10711.
Majumdar S, Debnath R, Bhattacharjee A, Banerjee A, Patro CS. Statistical optimization and characterization of prepared Fluconazole topical liposomal gel for improved skin permeation. Pelagia Res Lib 2014;5:4255.
Tatode AA, Patil AT, Umekar MJ, Telange DR. Investigation of effect of phospholipids on the physical and functional characterization of paclitaxel liposomes. Int J Pharm Pharm Sci 2017;9:1416.
About this article




Refbacks
 There are currently no refbacks.