Mechanism of the single-step nanoprecipitation method for the production of the hybrid nanoparticles
The nanoprecipitation method is one of the fast and repeatable ones for preparing the
HLPNPs (
Figure 2) (
25). This method was mainly based on the polymer precipitation from the lipophilic solution which was a combination of polar solvent and water. In particular, the used polymer (
PLGA) was precipitated as a hydrophobic core to encapsulate the less water-soluble drug. Notably, adding a lipid layer between the
PLGA polymer core and
PEG shell results in i) restriction of the release of small drug molecules from the polymer core which improves the encapsulation and loading efficiencies, and ii) reduction of the water penetration into the polymer core, and as a result, reduction in the hydrolysis rate of the
PLGA polymer which causes a slow drug release from the
NPS.
Validation parameters
Based on the linear equation taken from the melphalan calibration curve the values of σ and S were (0.0346 and 65.02 respectively). Consequently, according to Equation 3, the amount of LOQ was obtained at 5.32 µg/mL.
Analysis and optimization of the central composite design
Statistical analysis based on
RSM was used to predict the most appropriate model to describe the response surfaces (nanoparticle size and
EE%). Each experiment was repeated three times, and the response surface was determined in each experiment.
Table 2 presents the outputs. The results of the experimental design indicated that the designed system was affected by the amount of lipid, polymer, and
PVA, resulting in various drug
EE% and
NPs sizes. As shown in
Table 3 and
Table 4, the best state for each quadratic model response compared to the linear model was the quadratic two-factor model, which had the highest correlation coefficient (R
2). Thus the quadratic model was selected to describe binary interactions of independent variables on each response.
Effect on the size of the nanoparticles
The
NPs size is imperative to determine the efficiency of loaded therapeutics where
NPs with a smaller size have a higher chance to internalize into cells, resulting in the higher intracellular concentration of the loaded therapeutics (
26).
NPs with a size smaller than 300 nm are effectively up taken by target cells and exert their pharmaceutical activity (
27).
In the present study, the
NPs size for the 20 synthesized formulations was in the range of 92.01 – 194.03 nm. According to
Table 2, the smallest and largest
NPs were produced in experiments 20 and 10, respectively. The proposed quadratic model to describe the effect of the independent variables on the particle size is presented in Equation 5:
Equation 5.
Statistical results demonstrated that by increasing the amount of polymer (A) and lipid (B), the particle size with positive coefficients also increased, as shown in Equation 5.
Analysis of variance in the statistical calculations does not express the results with 100% certainty, thus the results are expressed with a percentage of probability. Therefore, by increasing the
F-value and decreasing the
P-value, the importance of the
NPs size would be significant and could be considered as a parameter affecting the process. Accordingly, concerning the results from data variance in
Table 5, all three independent variables were found to affect the
NPs size.
Figure.3 depicts the three-dimensional (3D) curve of the particle size response surface for a better understanding of the binary interaction of the independent variables on the response surface. According to
Figure.3a, at the high polymer concentrations, the particle size was risen by increasing the lipid concentration. This increase in the particle size might be due to increasing the viscosity of the solution, which in turn caused an increase in the liquid phase resistance of the particle dispersion. Consequently, the particle size can be increased by increasing the interconnection of particles among each other. Therefore, by increasing the number of particles, the interconnection rate between particles was increased, resulting in the production of larger
NPs (
28,
29). Also, by increasing the viscosity, the evaporation rate of the organic solvent was decreased, and particles with larger sizes were produced (
30,
31).
Figure 3b depicts that by increasing the initial concentration of the surfactant, the particle size decreased. This decrease in the particle size could be explained by the fact that at high concentrations of surfactant, the surfactant molecules tend to be accumulated and were sufficient to coat the
NPs. Hence, the surfactant activity increased and exert a significant effect on the
NPs size (
32,
33).
Effect on the drug entrapment efficiency
Encapsulation efficiency is a critical factor that affects the efficacy of the drug delivery carrier (
34). For example, low drug entrapment efficiency in polymeric carriers causes poor treatment outcomes and drug release properties and, as a result, insufficient efficacy of drug delivery systems (
19).
In the present study, EE% was calculated for all the prepared formulations to evaluate the effects of the independent variables (polymer, lipid, and PVA) on the drug EE% at different concentrations where the drug concentration was constant. Variation in the EE% at different concentrations was described by Equation 6:
Equation 6.
According to
Table 2, Experiments number 10 and 14 exhibited the highest and lowest
EE%, respectively. Regarding the importance of the
F-value parameter and statistical results reported in
Table 6, the two independent variables of polymer concentration (A) and amount of surfactant (C) had significant effects on the
EE% (
p-value < 0.05).
Figure 4 demonstrates the impacts of the independent variables on the drug
EE%. As can be seen in
Figure 4a, increasing the polymer concentration, the
EE% also increased. This could be a result of the fact that by increasing the polymer concentration, the encapsulation spaces for the drug also increased; hence, a relatively compressed matrix was created. Also, the hydrophobicity of
melphalan helped to achieve a high
EE%. According to the previous research, polymeric core and the drug hydrophobicity, as the two significant factors, promoted the
EE%. Furthermore, changes in the lipid content had no significant effect on
EE% and only affected the lipid thickness which resulted in a partial drug release from the nanoparticle (
6). As
Figure.4b demonstrated that by decreasing the polymer concentration and increasing the surfactant content, the
EE% decreased. These results from the fact that with increasing the surfactant concentration, the solubility of the drug from the organic phase to the aqueous phase increased, which caused a reduction in the viscosity, and consequently, a reduction in permeation during the process. This issue, in turn, led to a decrease in
EE% (
35).
Optimization
After evaluating the response surfaces by analyzing variance, numerical optimization was performed by applying optimal limitations and specifications for independent variables. Then optimal response surfaces were obtained with the values predicted by the software (
Figure 5). Then, the formulation of the
HLPNPs with the predicted values was prepared, and the
NPs size,
PDI (
Figure 6), and
EE% were measured and calculated.
Table 7 reports values predicted by the software and actual values of the response surfaces. According to the data presented in
Table 7, there was no significant difference between the actual values of the response surfaces and the predicted values. Therefore, the efficiency of encapsulation and loading of
melphalan was obtained to be 84.43 ± 3.67% and 2.98 ± 1.3%, respectively.
Moreover, the desirability of the optimized values was 0.947. In general, it could be concluded that those NPs used as a carrier for anticancer drugs, can easily reach the cell membrane and increase the drug concentration at the cell surface compared to the standard drug. consequently, such a condition increased the therapeutic effect of the anticancer drug. Hence, a little encapsulation efficiency would be of high importance.
Morphology and zeta potential of the optimal formulation
TEM was used to identify the morphology of the optimized hybrid
NPs.
Figure 7 shows images obtained from
TEM, indicating a small ring of lipid covering around the polymer core. Moreover, it could be seen that the prepared
NPs had a smooth surface, uniform, and integrated pattern with the spherical structure, suggesting a slow release of the drug. There was a difference, but not significant, between the particle size obtained from
TEM (94.41 nm) and those obtained by zeta sizer (96.24 nm) (
Figure 6). These results were somewhat consistent with the results of similar research, reported previously (
36,
37). Zeta potential determines the stability of the colloidal nanosuspensions.
Figure 8 depicts the zeta potential of the optimized value attributed to the nonionic nature of
PVA compared to the anionic nature of phosphatidylcholine. Additionally, negative zeta potential created a large repulsive force between
NPs, prevented their aggregation, and resulted in their stability (
38).
Drug release study (in-vitro)
Controlled-release drug delivery systems have remarkable advantages in comparison with conventional dosage forms. These systems 1) cause a significant decrease in the dosing frequency and provide more convenience for patients, 2) cause a minimum in the fluctuation of drug concentration
in-vivo and preserve the concentration of drugs within the proper range, 3) can deliver drugs site-specifically, and 4) can reduce the drug side effects (
19). Also, drug release from nanocarriers is a critical factor affecting the therapeutic outcome (
39).
Figure 9 demonstrates a pattern of sustained
melphalan release (for both the standard and encapsulated forms) at any time point. According to the Figure, the rate of drug release from the
NPs was much lower than that of the free drug release, indicating that the
NPs were able to encapsulate the drug, and release it in a controlled manner in that only 17.39% of the encapsulated drug was released after 48 h. The drug release from the
NPs was initiated with a burst release, in which 29% of the total release occurred in the first hour of the study. This could result from the release of the adsorbed drug to the
NPs surface. Releasing continued at a reduced rate until the end of the study. Lack of the rapid release of the drug from the
NPs suggested the proper interaction of
melphalan with
HLPNPs. Overall, the pattern of drug release from the
NPs indicated the potency of the particles as a controlled drug delivery system.
In-vitro cell viability and IC50 of nanoformulation
NPs can increase the therapeutic effects of anticancer drugs because of their capability to enhance the drug concentration in tumor cells.
NPs perform this by increasing drug circulation time. Moreover,
NPs can deliver drugs site-specifically
in-vivo, resulting in the restriction of the drug side effects (
26). In the present study, the cytotoxicity effects of
melphalan and
melphalan-loaded
HLPNPs against human ovarian cancer
A2780CP and
SKOV3 cells were evaluated. As
melphalan is used for the treatment of ovarian cancer; therefore,
A2780CP and
SKOV3 cells were used as
in-vitro models of the disease. The results demonstrated that the cytotoxicity of both formulations (
melphalan and
melphalan-loaded
NPs) was increased in a dose-dependent manner (
Figures 10 and
11). Also, the cytotoxicity was found to be cell type-dependent, as both formulations caused higher cytotoxicity in
A2780CP cells compared to
SKOV3 cells. However,
melphalan-loaded
HLPNPs were more potent compared to the standard drug (at the same drug concentration) to inhibit the growth of cancer cells, indicating the potency of the
HLPNPs to increase the cytotoxicity effects of
melphalan. Increasing the cytotoxicity effects of the
melphalan-loaded
NPs compared to
melphalan, resulted from the controlled drug release from the
HLPNPs. Moreover, the nanoformulation at high concentrations inhibited the toxic effects of
melphalan. This feature increases the maximum tolerable dose, and thus higher concentrations of the drug can be used, which in turn decreases the risk of tumor drug resistance. The cytotoxicity effects were also evaluated by calculating the half-maximal inhibitory concentration (
IC50): the drug concentration required to kill 50% of the cells incubated over the determined period).
Figure 12 shows the
IC50 values calculated for both cell lines over the intended time in comparison to the free drug. As seen in the Figure, the
IC50 values of the
HLPNPs were lower than that of the free drug over time because the free drug quickly passed through the cell membrane while the encapsulated form of the drug chose a specific pathway and released the drug in a controlled way.
Standard curve of melphalan
Scheme of the preparation process of the HLPNPs using the single-step nanoprecipitation
Three-dimensional curve of the effect of the independent variables on the response surface Y1. (a) Interaction of polymer and lipid concentration and (b) Interaction of polymer and polyvinyl alcohol concentration
Three-dimensional curve of the effect of the independent variables on the response surface Y2. (a) Interaction of polymer and lipid concentration and (b) Interaction of polymer and polyvinyl alcohol concentration
Optimization plots of the response surfaces and desirability of the optimal values
Size distribution of the HLPNPs with the values predicted by software
TEM micrograph of the HLPHNPs containing melphalan
Surface charge potential of the optimal Nano formulation
Release profile of Melphalan from HLPNPs free drug using the dialysis method within 48 h at 37 °C
Cytotoxicity effects of free melphalan and HLPNPs loaded with melphalan on the A2780CP cell line after 24, 48, and 72 h of incubation. Data is expressed as mean ± SD (n = 3). *(p < 0.05), ** (p < 0.01), and *** (p < 0.001) indicate a significant difference with HLPNPs
Cytotoxicity effects of free melphalan and HLPNPs loaded with melphalan on the SKOV3 cell line after 24, 48, and 72 h of incubation. Data is expressed as mean ± SD (n = 3). *(p < 0.05), **(p < 0.01), and ***(p < 0.001) indicate a significant difference with HLPNPs
IC50 of free melphalan, HLPNPs loaded with melphalan on the A2780CP and SKOV3 cell line after 72 h. Data is expressed as mean ± SD (n = 3)
| Independent variables Levels |
|---|
| | -1.68 | -1 | 0 | 1 | 1.68 |
| A: Polymer (mg.mL) | 3.3 | 4 | 5 | 6 | 6.7 |
| B: Lipid (mg.mL) | 2.3 | 3 | 4 | 5 | 5.7 |
| C: PVA (%) | 0.3 | 1 | 2 | 3 | 3.7 |
| Dependent variables | Constrains |
| Y1: Particle size (nm) | | | | Minimize | |
| Y2: Entrapment efficiency (%) | | | Maximize | |
| Std | Run | Polymer(mg.mL) | lipid(mg.mL) | PVA (%) | Size(nm) | EE (%) |
| 10 | 5 | 1.682 | 0 | 0 | 162.85 ± 5.5 | 94.5 ± 2.2 |
| 4 | 10 | 1 | 1 | -1 | 194.03 ± 2.9 | 94.45 ± 2.5 |
| 8 | 12 | 1 | 1 | 1 | 144.07 ± 3.5 | 93.85 ± 1.6 |
| 2 | 13 | 1 | -1 | -1 | 156.1 ± 5.1 | 91.5 ± 1.1 |
| 6 | 17 | 1 | -1 | 1 | 107 ± 2.4 | 90.55 ± 1.8 |
| 20 | 1 | 0 | 0 | 0 | 134.29 ± 2.6 | 88.45 ± 2.4 |
| 11 | 2 | 0 | -1.68 | 0 | 108.69 ± 1.9 | 88.45 ± 1.3 |
| 17 | 4 | 0 | 0 | 0 | 136.13 ± 1.1 | 88.65 ± 1.7 |
| 13 | 7 | 0 | 0 | -1.68 | 165.74 ± 5.1 | 90.01 ± 2.1 |
| 16 | 8 | 0 | 0 | 0 | 134.29 ± 2.5 | 87.89 ± 1.4 |
| 18 | 11 | 0 | 0 | 0 | 133.77 ± 3.3 | 88.5 ± 1.12 |
| 15 | 16 | 0 | 0 | 0 | 134.26 ± 2.8 | 88.4 ± 2.5 |
| 12 | 18 | 0 | 1.68 | 0 | 164.46 ± 3.1 | 89.5 ± 1.6 |
| 19 | 19 | 0 | 0 | 0 | 132.32 ± 3 | 89 ± 3.1 |
| 14 | 20 | 0 | 0 | 1.68 | 92.01 ± 2.4 | 88.58 ± 1.2 |
| 3 | 3 | -1 | 1 | -1 | 154.75 ± 4.4 | 81.48 ± 2.3 |
| 7 | 6 | -1 | 1 | -1 | 120.19 ± 3.3 | 80 ± 2.5 |
| 1 | 9 | -1 | -1 | -1 | 124.88 ± 4.1 | 86.5 ± 1.2 |
| 5 | 15 | -1 | -1 | 1 | 94.24 ± 2.6 | 84.45 ± 1.2 |
| 9 | 14 | -1.68 | 0 | 0 | 114.36 ± 1.8 | 76.45 ± 2.3 |
| Y1( size(nm)) |
|---|
| Source | Std.Dev. | R-Squared | Adjusted R-Squared | Predicted R-Squared | p-value | |
|---|
| Linear | 4.65 | 0.9728 | 0.9677 | 0.9506 | 0.0018 | |
| 2FI | 3.45 | 0.9879 | 0.9823 | 0.9679 | 0.0067 | |
| Quadratic | 1.84 | 0.9973 | 0.9949 | 0.9832 | 0.0953 | Suggested |
| Cubic | 1.95 | 0.9982 | 0.9943 | 0.7333 | 0.0239 | Aliased |
| Y2 (EE%) |
|---|
| Source | Std.Dev. | R-Squared | Adjusted R-Squared | Predicted R-Squared | p-value | |
|---|
| Linear | 1.85 | 0.8631 | 0.8375 | 0.7498 | 100.16 | |
| 2FI | 1.34 | 0.9418 | 0.9149 | 0.8722 | 51.17 | |
| Quadratic | 0.62 | 0.9905 | 0.982 | 0.9379 | 24.87 | Suggested |
| Cubic | 0.33 | 0.9984 | 0.9948 | 0.9957 | 1.72 | Aliased |
| Source | Sum ofSquares | df | MeanSquares | FValue | p-valueProb> F | |
|---|
| Model | 12698.91 | 9 | 1410.99 | 414.75 | < 0.0001 | significant |
| A-polymer | 2607.24 | 1 | 2607.24 | 766.38 | < 0.0001 | |
| B-lipid | 3694.58 | 1 | 3694.58 | 1086.00 | < 0.0001 | |
| C-PVA | 6084.58 | 1 | 6084.58 | 1788.53 | < 0.0001 | |
| AB | 46.04 | 1 | 46.04 | 13.53 | 0.0043 | |
| AC | 143.31 | 1 | 143.31 | 42.12 | < 0.0001 | |
| BC | 2.86 | 1 | 2.86 | 0.84 | 0.3811 | |
| A2 | 59.30 | 1 | 59.30 | 17.43 | 0.0019 | |
| B2 | 24.79 | 1 | 24.79 | 7.29 | 0.0223 | |
| C2 | 28.71 | 1 | 28.71 | 8.44 | 0.0157 | |
| Residual | 34.02 | 10 | 3.40 | | | |
| Lack of Fit | 26.55 | 5 | 5.31 | 3.55 | 0.0953 | not significant |
| Pure Error | 7.47 | 5 | 1.49 | | | |
| Cor Total | 12732.93 | 19 | | | | |
| Source | Sum ofSquares | df | MeanSquares | FValue | p-valueProb>F | |
|---|
| Model | 396.60 | 9 | 44.07 | 115.88 | < 0.0001 | significant |
| A-polymer | 341.34 | 1 | 341.34 | 897.60 | < 0.0001 | |
| B-lipid | 0.15 | 1 | 0.15 | 0.41 | 0.5378 | |
| C-PVA | 4.10 | 1 | 4.10 | 10.79 | 0.0082 | |
| AB | 30.89 | 1 | 30.89 | 81.23 | < 0.0001 | |
| AC | 0.49 | 1 | 0.49 | 1.29 | 0.2828 | |
| BC | 0.11 | 1 | 0.11 | 0.28 | 0.6094 | |
| A2 | 16.50 | 1 | 16.50 | 43.39 | < 0.0001 | |
| B2 | 0.40 | 1 | 0.40 | 1.06 | 0.3270 | |
| C2 | 1.13 | 1 | 1.13 | 2.98 | 0.1148 | |
| Residual | 3.80 | 10 | 0.38 | | | |
| Lack of Fit | 3.15 | 5 | 0.63 | 4.81 | 0.0550 | not significant |
| Pure Error | 0.66 | 5 | 0.13 | | | |
| Cor Total | 400.40 | 19 | | | | |
| Polymer (mg.mL) | Lipid (mg.mL) | PVA (%) | Size (nm) | EE (%) | Desirability |
|---|
| predicted formulation | 4.8 | 3 | 3 | 95.77 | 87.73 | 0.94 |
| Actual Optimized formulation | 4.8 | 3 | 3 | 96.24 ± 2.62 | 83.43 ± 3.67 | …………… |