3.1. Factorial Designs
As previously mentioned in the methodology, the selection of excipients for the NLC core and their proportions were determined in a previous study through Raman mapping and DOE [
1]. From this previous study, it was established that Crodamol
TM CP would be used as a solid lipid, SR
TM DMI, and SR
TM Lauryl Lactate as a liquid lipid (SR
TM DMI is indeed a solvent, however since it goes within with lauryl lactate, it will be denoted here ‘liquid lipid’ as well) and several drugs were evaluated, including Butamben. The proportions between Crodamol
TM CP and Super Refined
TM Lauryl Lactate were set at 1:1.6 (w/w), and the ratio of SR
TM DMI to Super Refined
TM Lauryl Lactate was 1:1 (w/w).
The 2
3 factorial design aimed at identifying the optimal percentage of total lipids, the most effective surfactant and its percentage (critical material attributes - CMAs), and the optimal sonication time (critical process attribute - CPP) that would allow the achievement of suitable particle size, PDI, and zeta potential – considered critical quality attributes (CQAS).
Table 2 outlines the selected independent and dependent variables along with the specified objectives for particle size, PDI, and zeta potential.
The selected surfactants for the study were Crodasol
TM HS HP, Synperonic
TM PE/F68, and Croduret
TM 40. Through the 2³ design for each surfactant, 11 different experimental combinations were obtained for formulations and subsequent analysis (
Table 3).
3.1.1. Interpretation of Responses Obtained for the DOE Using CrodasolTM HS HP as Surfactant
The responses obtained for the design with Crodasol
TM HS HP as a surfactant are presented in
Table 4. The regression was significant for particle size and PDI and not for zeta potential, i.e. the deliberate experimental variations did not affect this response. Linear models were well adjusted for both responses. The effect of excipients on the responses of interest is detailed in
Table 4 and
Figure 1, and the residue analysis of linear models is depicted in
Figure 2.
As can be seen in
Table 4, formulation particle size ranged from 244.93 nm to 698.93 nm, PDI from 0.176 to 0.349, and zeta potential from 25.66 to 29.26.
For the Particle size response, the R² was 0.9987. Residual analysis graphs (Fig.2.A-B) indicated that the residuals follow a normal distribution, are homoscedastic and independent. Additionally, Fig. 2C shows an excellent agreement between the predicted and actual values, thus indicating the model’s good predictive capability. As described in
Table 5, the particle size response exhibited two significant main effects and three significant interactions. Analyzing the significant effects, Crodasol
TM HS HP, % surfactant (X
2), and the total lipid (X
1) variables showed positive effects. The positive effect of the surfactant was unexpected since, typically, in such formulations, the surfactant effect is negative, i.e., as the concentration increases, the size decreases. A possible explanation is that higher concentrations of Crodasol
TM HS HP result in a predominance of hydrophobic interactions due to non-polar groups, such as the long alkyl chains of the hydroxyl stearate portion of the surfactant, the interactions between the lipids and the surfactant, causing an increase of size [
26]. This outcome was also observed by Almousallam and collaborators, who have used Kolliphor® P 188 as a surfactant [
27]. The positive effect of the total lipid was attributed to the increase in the lipid concentration causing lipid particle coalescence, and therefore an increase in particle size [
26].
For the PDI response, the R² was 0.7348. The residual analysis graphs (Fig.2.D-F) show that the model residuals follow a normal distribution, are homoscedastic and independent. As indicated in
Table 5, the PDI response showed only surfactant (X
2) as a significant effect. The increase in the concentration of Crodasol
TM HS HP is responsible for the increase in the PDI of the formulations, this increase may be because the increase in surfactant causes particle aggregation, resulting from the surfactant binding on the surface of the nanoparticle [
28].
Figure 3A shows the 3D surface and it becomes evident that particle size is influenced by the interaction between variables X
1 and X
2. Higher levels of the variables X
1, X
2, result in larger particle sizes, whereas lower levels are associated with smaller particle sizes. Consequently, to achieve smaller particle sizes (indicated by the blue region in the graph), it would be ideal to work with a lower % of surfactant (5%), a lower % of total lipids (10%), and a shorter sonication time (5min). The shorter sonication time is due to the interaction between this variable and variable X2.
Figure 3B shows the surface graph which clearly indicates that the PDI is influenced solely by variable X
2. Therefore, to achieve lower PDIs (indicated by the blue region in the graph), it is advisable to work with a lower % of surfactant (5%).
After processing the data and generating the mathematical model for each property of interest (particle size, PDI, and zeta potential), it was possible to define the optimized formulation for each surfactant used. Desirability was used with the following criteria for the responses: minimize particle size and PDI and maximize zeta potential. The desirability graph indicates the formulation that best meets the desired criteria (the closer to 1 the better).
For the formulation containing Crodasol
TM HS HP as a surfactant, the desirability was 0.892 (
Figure 3C). Therefore, for the Crodasol
TM HS HP surfactant, the optimized formulation was defined with 10% Total Lipids, 5% surfactant, and 7.5 minutes of sonication time. The formulation was prepared again and monitored in the stability study, with subsequent characterization.
3.1.2. Interpretation of Responses Obtained for the DOE Using SynperonicTM PE/F68 as Surfactant
The obtained responses for the design using Synperonic
TM PE/F68 as a surfactant are presented in
Table 6. The models were significant for particle size and zeta potential whereas the experimental variations did not affect PDI. This indicates the robustness of the CQA regarding the chosen parameters. The lack of fit was nonsignificant for particle size and zeta potential.
For the Particle size, the R² was 0.8465. The residual analysis graphs (Fig.5.A-B) show that the model residuals follow a normal distribution, are homoscedastic, and independent of each other. Additionally Fig. 5C demonstrates the good agreement between predicted and actual values, indicating the predictive capacity of the model. The particle size of the formulations ranged from 157 nm to 260 nm and, as shown in
Table 7 and
Figure 4, only the surfactant Synperonic
TM PE/F68 (X
2) has a significant (negative) effect. This negative effect for the surfactant is expected, as the surfactant reduces the surface tension between the lipid and aqueous phases, resulting in the formation of smaller nanoparticles and preventing the coalescence of larger droplets [
15,
28]. This is clear when analyzing
Figure 6A, which includes the surface graphs, Therefore, to obtain smaller particle sizes (indicated by the blue region in the graph), it would be optimal to work with a higher % of surfactant (10%).
For the Zeta Potential response, the R² was 0.9116. In the residual analysis graphs (
Figure 5.D-E) it is observed that the model residuals follow a normal distribution, are homoscedastic, and independent of each other. Fig. 5F shows a good agreement between predicted and actual values, thus indicating the good predictive capability of the model. As shown in
Table 7 and
Figure 4, the Zeta Potential response presented surfactant and sonication time as a significant effects. The zeta potential of the formulations varied between |4| mV and |19.7| mV.
When analyzing
Figure 6.B, which includes the surface graphs, it is noted that the variables that affect the zeta potential are surfactant (X
2) and sonication time (X
3). The variable X
2 is associated with an increase in zeta potential when used at its lowest level and the variable X
3 is associated with an increase in zeta potential when used at its highest level. Therefore, to obtain higher potentials, it would be optimal to work with a lower % of surfactant (5%), a lower % of total lipid (10%), and a longer sonication time (10 minutes).
For the formulation containing SynperonicTM PE/F68 as a surfactant, the desirability was 0.656 (
Figure 6.C). It is worth noting that
Figure 6.C is shown in full green, as the particle size effect presented the desirability criterion (size smaller than 260 nm) at all tested points. Therefore, for the surfactant SynperonicTM PE/F68, the optimized formulation was defined with 20% Total Lipids, 10% surfactant, and 10 minutes of sonication time. The formulation was prepared again and monitored in the stability study, with subsequent characterization.
Figure 5.
Graphs containing residue analyses for each planning response containing SynperonicTM PE/F68 as surfactant. A-C: residue analysis referring to the Particle Size response. D-F: residue analysis referring to the Zeta Potential response.
Figure 5.
Graphs containing residue analyses for each planning response containing SynperonicTM PE/F68 as surfactant. A-C: residue analysis referring to the Particle Size response. D-F: residue analysis referring to the Zeta Potential response.
Figure 6.
Surface graphs for the responses of the formulation containing SynperonicTM PE/F68 as surfactant: A – Particle size. B – Zeta Potential and C – desirability graph.
Figure 6.
Surface graphs for the responses of the formulation containing SynperonicTM PE/F68 as surfactant: A – Particle size. B – Zeta Potential and C – desirability graph.
3.1.3. Interpretation of Responses Obtained for the Experimental Design Using CroduretTM 40 as Surfactant
The initial objective was to carry out experimental designs with all the surfactants selected for the project. However, when using Croduret
TM 40 in the formulations, phase separation occurred shortly after preparation, indicating that Croduret
TM 40 could not form stable NLCs. This could have happened due to its structure, it is an ethoxylated surfactant of plant origin, non-ionic, with an HLB of 13, normally used as an excipient in self-emulsifying drug delivery systems (SEDDS) [
29]. Therefore studies with this surfactant were not continued.
3.2. Comparison between Formulations with CrodasolTM HS HP, and SynperonicTM PE/F68
Table 8 presents a summary of the significant effects found for each response in each design, allowing a comparison between the use of different types of surfactants. The variation of the factors was consistent (
Table 2 and
Table 3). Therefore, the same % total lipids, % surfactant, and sonication time were consistently utilized.
SynperonicTM PE/F68 has the same chemical structure of Poloxamer 188 and consists of triple copolymers formed by Polyethylene glycol-b-polypropylene glycoly-polyethylene glycol (PEO-PPO-PEO) – with 16% PPO and 84% PEO and has an HLB >24. On the other hand, Crodasol
TM HS HP is synonymous with Macrogol 15- Hydroxystearate and consists of a mixture of monoesters and diesters of 12-hydroxystearic acid and macrogols obtained through the ethoxylation of 12-hydroxystearic acid, having an HLB of 15. It is observed that the difference in chemical composition between the two surfactants is also reflected in the values of particle size, PDI, and zeta potential.
Figure 7 displays the size, PDI, and zeta graphs of the three surfactants used in formulations with BTB.
Analyzing
Figure 7, it was observed that formulations with the Crodasol
TM HS HP surfactant resulted in larger particle size (244 nm – 698 nm), higher PDI values (0.176 – 0.319), and elevated zeta potential values (|25 -29|), in comparison to formulations containing Synperonic
TM PE/F68 (size: 157-260 nm. PDI: 0.091 - 0.172. Zeta: |4.4 – 19| mV) as the surfactant.
Regarding the difference in particle size between Crodasol
TM HS HP and Synperonic
TM PE/F68, it can be influenced by the hydrophilic-lipophilic balance (HLB) value. Crodasol
TM HS HP has an HLB of 15, while Synperonic
TM PE/F68 has an HLB of > 24. Since the solid lipid used is Crodamol
TM CP, characterized by high hydrophobicity and molecular weight, it contributes to the viscosity of the system. Consequently, this requires a surfactant with a higher HLB [
30].
3.4. Particle Size Analysis, Polydispersity Index, and Zeta Potential Determination
The optimized formulations, as outlined in
Table 11 and
Table 12, were characterized by particle size, PDI, and zeta potential.
Initially, the formulations were analyzed both with and without SR
TM DMI in the presence of BTB. These results are presented in
Figure 8.
It can be observed that the removal of SR DMITM from the formulations changed particle size and zeta potential, except for the formulation containing CrodasolTM HS HP as a surfactant, in which the zeta potential of the system did not change after removing SRTM DMI. It can be inferred that SRTM DMI helped stabilize the system. The optimized formulations without SRTM DMI exhibited larger particle sizes and lower zeta potentials, generating potentially less stable particles. NLCs with SRTM DMI, having smaller particle sizes, may suggest that SRTM DMI preferred to remain in the lipid nucleus, alongside Super RefinedTM Lauryl Lactate, a compound of medium polarity, rather than leaving the nucleus and entering the aqueous phase of the NLC.
Studies of particle size, PDI, and zeta potential with the optimized formulations without the presence of BTB were also carried out. Results are shown in
Table 12.
Table 12 shows that the size of NLC particles containing BTB is larger than that of empty NLCs, and this increase in size may indicate drug encapsulation within the particles. The PDI of NLCs with BTB (CBTB-NLC and SBTB-NLC) decreased slightly compared to S-NLC and C-NLC, remaining below 0.3, which indicates a uniform size distribution. The zeta potential of NLCs with BTB also decreased in formulations containing the two different surfactants, suggesting that the drug contributes to stability through electrostatic repulsion.
Concerning the Zeta Potentials of NLCs with BTB, it is noteworthy that the zeta values obtained from NLCs formulated with Synperonic
TM PE/F68 did not reach zeta values below 30 mV. It is known that surfactants of the Poloxamer 188 type confer protective effects due to the hydrophilic corona present in their structure, resulting in zeta potential values that may not be sufficient for electrostatic stabilization. Through electron microscopy, Bhattacharya and co-authors that Kolliphor
® P188 provides additional steric and electrostatic stabilization of formulations described it. The presence of this surfactant enhances steric stabilization on the particle surface, even in cases where the charge may be insufficient to generate electrostatic stabilization [
31].
According to Shah and co-authors, NCLs can demonstrate enhanced stability, even with a lower zeta potential value, attributed to steric repulsion. This repulsion may result from the adsorption of polymers or non-ionic surfactants on the lipid matrix surface, effectively preventing potential particle aggregation [
32].
3.5. Nanoparticle Tracking Analysis
NTA is a technique for measuring size distribution of particles in the sample through the measurement of light scattering properties and Brownian motion. The process involves a laser beam passing through the sampling chamber, where suspended particles in the beam’s path scatter the light. Additionally, it enables the determination of nanoparticle concentration (n° particles/mL) through specific and individual particle counting concerning volume [
24].
The results of the NTA analyses are presented in
Table 13.
It is observed that all NCLs exhibit particle concentrations on the order of 10¹³ particles/mL, consistent with values found in the literature [
15,
16,
33]. An important parameter to consider in NTA is the Span, which should be lower than 1 to indicate a homogeneous particle size distribution [
24]. Formulations containing Crodasol
TM HS HP (CBTB-NLC and C-NLC) displayed a Span value greater than 1, indicating a wider range of particle size distribution when compared to Synperonic
TM formulations.
When comparing the number of particles between NLCs with BTB and empty NLCs, it is observed that the addition of BTB to the formulations increased the number of nanoparticles in suspension. The same observation was made by Guilherme and collaborators, who utilized Poloxamer 188 as a surfactant and cetyl palmitate as a solid lipid [
15].
3.6. Differential Scanning Calorimetry (DSC)
Through this technique, information about the crystallinity of the particles and the interactions between the drug and lipids in the formulation were obtained.
Figure 9 displays the DSC curves of the lyophilized NLCs, BTB, Crodamol
TM CP, and the two surfactants used. Butamben, Crodamol
TM CP, and Synperonic
TM PE/F68 exhibited endothermic events at 58.10°C, 55.58°C, and 55.81°C, respectively, corresponding to the melting point of these compounds. These results are in accordance to what is described for these compounds in literature [
22,
34,
35].
In
Figure 9, it can be seen that the C-NLC and CBTB-NLC curves present endothermic events at 52.57°C and 54.84°C, respectively. These events are possibly related to the Crodamol™ CP transition. When comparing the enthalpy of fusion and temperature values, a decrease in temperature and enthalpy is noted in both events compared to the temperature and enthalpy of fusion of bulk Crodamol™ CP. This change can be attributed to the nanometer size of the NLC, as explained by the Gibbs-Thomson effect, which postulates that the smaller the particle size, the greater the surface area, leading to a decrease in the melting point of the NLC compared to the solid lipid, and due to the presence of surfactants in the NLC. Furthermore, the decrease in the endothermic event of Crodamol™ CP in NLCs may come from an interaction between the drug and Crodamol™ CP (solid lipid), which may alter the lipid's fusion behavior [
36,
37,
38,
39,
40].
Similar events are observed in the S-NLC and SBTB-NLC curves (50.56°C and 51.9°C). However, in the case of these NLCs, the surfactant used also presents an endothermic event at 55.81°C (Synperonic™ PE/F68). In these cases, the events observed in NLCs may result from the combination of events related to the solid lipid and the surfactant [
40].
When comparing the M-NLC and CBTB-NLC formulations, a decrease in the fusion enthalpy and temperature in the CBTB-NLC is observed, which is expected, as the physical mixture is merely all the solid components present in the NLC, only mixed. The same occurs when comparing M-NLC with SBTB-NLC. This indicates that M-NLC presents more ordered crystal lattice arrangements, unlike NLCs (CBTB-NLC and SBTB-NLC), in which the solid lipid is in less ordered arrangements, contributing to the encapsulation of the drug [
37,
41,
42].
Based on the presented results, it is evident that all NLCs exhibited a decrease in fusion enthalpy, indicating lower lipid crystallinity, minimizing drug expulsion during the storage period [
41].
3.7. FTIR
Infrared spectroscopy is used to explore the structural characteristics of lipids and identify potential interactions between the drug and the excipients in the formulations [
43,
44].
Analyzing
Figure 10 it is possible to compare the spectra of the constituents of the formulations with those of the NLCs and with the M-NLC, which is the physical mixture of the developed NLCs. It can be observed that the stretch and deformation signals are consistent with those present in the spectra of Crodamol
TM CP, the solid lipid used, and the bands of the surfactants. In the case of S-NLC, and SBTB-NLC, their spectra include characteristic bands of their respective surfactant Synperonic
TM PE/F68. The same applies to M-NLC, which has Synperonic
TM PE/F68 as a surfactant. Bands in the regions of 1735 cm
-1, 1250-970 cm
-1, and 3000-2800 cm
-1 are present in both Crodamol
TM CP and the nanocarrier spectra. Regarding surfactants, bands at 1281-1239 cm
-1, and 960-942 cm
-1 are found in the spectra of both surfactants used and also in the spectra of the nanocarriers. It is noteworthy that the band in the region between 1275-1245 cm
-1 (O-C stretching) is present in the BTB spectrum, however, when observing this band in the nanocarrier spectra, it is noted that both the NLCs loaded with BTB and empty NLCs exhibit similar bands.
When comparing the BTB spectrum with that of the charged NLC (
Figure 11), it is evident that the bands at 3420-3341 cm
-1 (N-H stretching), 1675-1638 cm
-1 (C=C stretching), and 1589 cm
-1 (deformation NH2) are absent in the spectra of the charged NLCs, due to the lower concentration of BTB. However, the bands at 1517 cm
-1 (C=C stretching – in ring) and 839-700 cm
-1 (C-H bending and ring puckering deformation) are identifiable in both the drug spectrum and the BTB-NLCs spectra.
When comparing the spectra of empty NLCs with BTB-NLCs, a similarity in the obtained spectra is observed. This similarity may arise from the mechanical incorporation of the drug into de particle, with only some potential interactions between the components of the formulation and the drug being detected by the technique [
37,
45,
46,
47].
Upon analyzing the spectra of BTB-NLCs and M-NLCs (
Figure 11), a distinctive band related to O-H stretching (3363-3360 cm
-1), is observed, exclusively present in these formulations. This band may be a shifted version of a BTB band, as it is absent in the spectra of empty NLCs. There is a high similarity between the M-NLC spectra and the CBTB-NLC and SBTB-NLC spectra regarding the frequency of bands present (bands 3361, 1603, and 1275 cm
-1). However, some of these frequencies have lower intensity in the CBTB-NLC and SBTB-NLC formulations, which may indicate an interaction between the NLC and the drug.
3.8. 3-Month Stability Study of the Optimized NLCs
The evaluated stability parameters were particle size, PDI, and zeta potential.
Figure 12 presents the stability results of the formulations.
Based on the results presented in
Figure 12 it is observed that in the presence of the drug, the NLCs that have Crodasol
TM HS HP in their composition did not show stability over 2 months, presenting phase separation. The only exception is the CBTB-NLC without SR
TM DMI
. NLCs with Synperonic
TM PE/F68 even with lower zeta potentials, presented higher stability than NLCs with Crodasol
TM HS HP. This is probably because Synperonic
TM PE/F68, provides additional steric and electrostatic stabilization within the NLC and even on the surface of the particles [
31].
The study makes us wonder whether it is really important to achieve this type of NLC stability in the form of a dispersion in water since this type of formulation in the liquid state can suffer instability due to oxidation reactions, hydrolysis, phase inversion, crystallization, and polymorphism, which can result in particle aggregation, gelation, and affect drug release [
48,
49].
Studies report that the ideal approach for NLC would be to freeze-dry the system or add preservatives, so that it obtains long-term stability, preserving the NLC from bacterial contamination and maintaining the properties of the nanoparticles, such as particle size, in addition to preventing degradation reactions (hydrolysis) from occurring [
28,
48,
49,
50].