Influence of emulsification mode, stirring speed and volume on ibuprofen-loaded PLGA microparticles

Overview

Microparticles based on biodegradable synthetic lactic acid and glycolic acid copolymer (PLGA) were successfully prepared by the solvent evaporation method. Ibuprofen was chosen as the model drug. Various formulation and process parameters have been used to prepare each sample with emphasis on size reduction. The effect of the emulsification method (direct emulsification or emulsification using an ULTRA-TURRAX or a NE-1000 dispenser), the volume of the aqueous phase (200, 800 ml) and the stirring speed of the emulsion system (600, 1000 rpm) on the characteristic properties of microparticles, such as encapsulation efficiency, drug loading and particle morphology, was observed. The resulting microparticles were evaluated by optical microscopy or laser diffraction and the dissolution test was performed. It was found that the sample prepared by direct emulsification with 800 ml of an aqueous phase at 600 rpm provided the most favorable results, meanwhile the emulsification pre-step using a homogenizer caused promising particle size reduction. Gradual emulsification was evaluated as inapplicable due to great losses.

Key words: microparticles • solvent evaporation • PLGA • ibuprofen • size reduction

Keywords:

microparticles – solvent evaporation – PLGA – ibuprofen – size reduction

Introduction

Microdispersion and nanodispersion dosage forms have a prominent place in contemporary research of controlled release^1–4). Microparticle drug delivery systems are generally used to prolong drug release, to enhance their stability, to target a specific site and to improve bioavailability⁵⁾. With a growing need of drug incorporation into polymeric materials, several chemical and physical-mechanical methods have been gradually developed. The chemical methods include coacervation, interphase polymerization and polycondensation and cross-linking methods. Physical-mechanical methods include, for example, extrusion and spheronization, spray drying and cooling, or molding and coating of microparticles in the fluidized bed. One of the most frequently used physical-mechanical methods is the solvent evaporation method, which is very important for the preparation of poly(lactic-co-glycolic acid) (PLGA) microparticles⁶⁾. PLGA is a biodegradable polymer with very well-known degradation properties⁷⁾ which has been extensively utilized for controlled drug delivery systems⁸⁾ and belongs to the best-defined biodegradable materials available for similar systems⁹⁾.

The basic method principle is the evaporation of that emulsion portion which contains a dissolved or dispersed drug and polymer. During evaporation, the dissolved polymer solidifies and forms a matrix which entraps the drug in the particle structure. Microparticle preparation by a simple way of this method can be divided into several steps¹⁰⁾. In the first stage, the polymer is dissolved in a volatile organic solvent which is immiscible with water¹¹⁾. Subsequently, the drug is dissolved or dispersed in this polymer solution. The resulting lipophilic phase is then emulsified into an aqueous continuous phase, usually combined with emulsifiers, to form a fine O/W emulsion (O is the reference for the oil phase and W is the water phase). To achieve the required emulsification level, for example, a propeller stirrer with adjustable speed, a homogenization device or an ultrasonic bath are used¹⁰⁾. The organic solvent subsequently diffuses into the aqueous phase and vaporizes in the last step on the water/air interface. Lack of solvent and stirring causes a formation of microparticles from the polymer. The resulting microparticles, suspended in the continuous phase, are then filtered, washed and dried¹²⁾.

Preparation process is influenced by a great number of both the process and formulation parameters, including stirring speed, emulsification approach, or ratios of used excipients. Altering these variables can seriously affect encapsulation efficiency, yield and particle size, resulting in different dissolution profiles^{12, 13)}. The aim of this study was to prepare ibuprofen-loaded PLGA microparticles by the solvent evaporation method and to evaluate the influence of the outer aqueous phase volume, stirring speed and the mode of the emulsification on the measured particle parameters, mainly encapsulation efficiency, drug loading and particle size, and on dissolution profiles. Drug release profiles were then analyzed using kinetic model equations to approximate the drug release mechanism.

Experimental part

Materials

Ibuprofen (Zentiva, Czech Republic) served as the model drug, poly(lactide-co-glycolide) acid PLGA; L-lactide/glycolide = 50/50 (Resomer® RG 504 H, Boehringer Ingelheim Pharma GmbH & Co., Germany) was used for the formation of the polymer matrix in the oil phase. Dichloromethane – DM (Penta, Czech Republic) was as the organic solvent used for the oil phase and polyvinyl alcohol – PVA (Mw 31.000–50.000) (Sigma Aldrich, USA) served as the emulsifier. Phosphate buffer of pH 6.8 for dissolution test was prepared from sodium phosphate dodecahydrate and potassium dihydrogen phosphate (both Merck KGaA, Germany). All materials were of Ph. Eur. quality.

Microparticle preparation

Microparticles were prepared by the single emulsion (O/W) solvent evaporation technique. The samples preparation differed in the formulation parameter of the outer phase volume (200 or 800 ml), and the process parameters, namely the stirring speed (600 or 1000 revolutions per minute) and the mode of emulsification – direct emulsification, direct emulsification with a pre-emulsification step and emulsification using a NE-1000 dispenser (New Era Pump Systems, USA) with a dispensing rate of 1.12 ml/min. The samples characteristics and their designation are listed in Table 1.

For the formation of the oil phase, 200 mg of ibuprofen and 700 mg of PLGA were dissolved in 5 ml of dichloromethane. In the next step, the thus formed O phase was emulsified into the aqueous phase of a 0.1% PVA solution by one of the three above-mentioned emulsification modes. The organic solvent was evaporated under a mechanical stirrer (Heidolph RZR 2021, Sigma Aldrich, USA) for 90 minutes at 450 rpm. After evaporation, the prepared microparticles were collected using an 80 µm mesh sieve, washed three times with purified water and then dried at 25 °C in a cabinet drier for 24 hours (HORO – 048B, Dr. Hofmann GmbH, Germany). To compare the results with our previous studies¹⁴⁾, sample A was prepared as a reference under previously tested conditions. The preparation of the last sample (I) was similar to other samples; however, the internal phase was directly emulsified by a rotor-stator homogenizer ULTRA-TURRAX (T25 basic, IKA-Werke, Germany) for 60 seconds at 10 000 rpm to ensure a formation of a fine micro emulsion (pre-emulsion step). As a result of this modification, at the end of the preparation, the microparticles were not detectable by the naked eye. Isolation of the resulting micro suspension by sieve was not possible, but it was accomplished by centrifugation (EBA 20 Hettich, Germany) at 6000 rpm for 2 minutes. Particles were re-suspended in a small amount of water and collected by filtration on a Buchner funnel with membrane filter paper. Each sample was prepared in quadruplicate, combined after the collection to represent one sample. A total of 9 PLGA samples were prepared.

Microparticle characteristics

Drug content analysis

The ibuprofen content in PLGA microparticles was determined using an UV/Vis spectrometer (Lambda 25, Perkin Elmer Instruments, USA). Each sample was prepared by weighing an exact quantity of dried microparticles (50 mg) into a 50 ml volumetric flask with dichloromethane added to the mark. Each sample was analyzed in triplicate, absorbance was measured at 264 nm (the absorption maximum for ibuprofen) and the obtained results were expressed as mean values and their standard deviations. According to the data obtained from calibration curve, the drug content in the microparticles was determined.

The obtained values also served to determine encapsulation efficiency (EE) [1], drug load (DL) [2] and practical yield [3] by using the equations below^15–17).                                                

 

                                                                                                     

where w₁ represents the actual weight of the drug in microparticles, c_t is the theoretical amount of the drug, w₂ is the total weight of prepared microparticles and w_t is the theoretical yield (total amount of the drug and polymer used for the microparticle preparation).

Optical microscope analysis

Morphological properties of the prepared micro-particles such as sphericity factor and equivalent size distribution were evaluated by a NIKON SMZ 1500 optical stereomicroscope (Nikon, Japan) and a 72AUC02 USB camera (The Imaging Source, Germany). Randomly selected 200 microparticles were evaluated by the computer software NIS-Elements AR 4.0 (Nikon, Japan). Equivalent diameter and sphericity were calculated from the measured values according to equations [4, 5] and expressed as arithmetic mean ± standard deviation^{18, 19)}. Due to a smaller size, optical analysis of sample I was performed using a NIKON ECLIPSE E200 (Nikon, Japan) optical microscope. A picture of each sample was taken.                                 

 

where S means sphericity, A is the area of particles and p their perimeter.

Scanning electron microscopy

The morphology and dimension of sample I was also evaluated by scanning electron microscopy (SEM; MIRA3, Tescan Orsay Holding, Czech Republic) equipped with a secondary electron detector (SED). The sample was mounted on a SEM specimen stub using carbon conductive double-faced adhesive tape (Agar Scientific, United Kingdom). To eliminate charging artefacts, microparticles were coated with a 20 nm layer of gold using the metal sputtering coating method with argon atmosphere (Q150R ES Rotary-Pumped Sputter Coater/Carbon Coater, Quorum Technologies, United Kingdom). SEM images were obtained at accelerating voltage of 3 kV.

Laser diffraction

To determine sample I particle size distribution, laser diffraction was performed using a Mastersizer 2000 (Malvern, United Kingdom). To obtain a clear signal, approximately 0.5 g of sample was poured into a beaker filled with 300 ml of degassed purified water. The sample was measured immediately and it was analyzed for volume-weighted size distribution. The measurement was performed in triplicate.

Multivariate data analysis

For evaluation of process/formulation variables (stirring speed, aqueous phase volume and emulsification method) and their interaction with response variables (EE, DL, yield, diameter), factor analysis was used, including rotation of factors (Varimax normalized). The data matrix for multivariate analysis was composed of 8 measurements/objects. Prior to modelling, variables were automatically adjusted by autoscaling. Statistical evaluation of data was performed using the program Statistica 12 (StatSoft, USA).

In vitro release studies

The samples were weighed to obtain 15 mg of the drug, placed in 500 ml of 6.8 pH buffer and tested in an automated dissolution device (SOTAX AT 7 SMART On-Line System, Donau Lab, Switzerland) equipped with baskets. The temperature was kept at 37.0 ± 0.5 °C and the stirring speed at 75 rpm. Sampling was performed at pre-set time intervals: every 15 minutes for 1 hour, every 60 minutes for the rest of the first 24 hours and every 120 minutes for next two days to obtain 72-hour profiles. Analysis of the samples was performed by a UV spectrophotometer (Lambda 25, PerkinElmer, USA) at 222 nm. The measurement result was a dissolution curve expressing the cumulative drug release in time. Sustained drug release rate was expected since ibuprofen is practically insoluble in water. The obtained data were correlated with the equations of the drug release mathematical models²⁰⁾:

 

where M₀ is the initial amount of the drug; M_t is the amount of the drug released in time t; M_∞ is the absolute cumulative amount of the drug released at an infinite time; K₀, K₁, K_H, K_KP, K_S and K_BL are the zero order, first order, Higuchi, Korsmeyer-Peppas, Hixson-Crowell and Baker-Lonsdale release constants. Release exponent n of the Korsmeyer-Peppas model describes the mechanism of the drug release: n = 0.5 corresponds with the Fickian diffusion, 0.5 < n < 1.0 suggests an anomalous transport, n = 1.0 means the zero-order release kinetics and n > 1.0 suggests the super Case II transport.

The similarity factor f₂ [11] indicates the similarity percentage between two dissolution curves. The range of this factor can vary from 0 to 100. If a similarity factor is 100, it means that the dissolution profiles are identical²⁰⁾. If f₂ values range between 50 and 100, the dissolution profiles are evaluated as similar²¹⁾. If the value is 50 or higher (f₂ ≥ 50), then it can be said that the drug release profiles are more than 90% similar. While if this value is lower than 50 (f₂ < 50) then the profiles are not similar and the observed influence of the process variable can be considered as significant.

where T_i is the drug amount (%) released at time interval i in the tested sample, R_i is the drug amount (%) released at time interval i in the reference sample, n is the total number of samplings and w_i is optional weight factor at time interval i.

Results and discussion

Encapsulation process, drug loading and yield

The encapsulation efficiency, drug loading and yield results are shown in Table 2. The encapsulation efficiency ranged between 15.69 ± 0.12 % and 54.07 ± 1.22 %, the drug loading took values from 15.19 ± 0.12 to 20.15 ± 0.32 and the yield was from 18.60 % to 61.67 %. From the results it is evident that the emulsification using a dispenser provided values vastly inferior to the direct emulsification method, primarily the yield values from which the poor encapsulation efficiency also resulted. It is possible that the gradual addition of the O phase disrupts already formed half-made particles, which leads to higher material losses through polymer merging. From Table 2 it can be also seen that at 600 rpm the samples prepared with a higher volume of the W phase provided much better values of encapsulation efficiency and yield compared to their corresponding samples prepared with 200 ml outer phase.

Surprisingly, at 1000 rpm this trend is not evident and sample G prepared with 200 ml outer phase gives even better results than its corresponding sample. It could mean that stirring rate co-influences the observed characteristics together with the volume of the outer phase. This is further confirmed by interpreting the results from the perspective of the stirring rate. In the corresponding samples prepared with 800 ml outer phase the samples prepared at 600 rpm provide better results, however this trend is missing in the samples prepared with 200 ml outer phase.

Results for sample I stood apart the other samples. Emulsification pre-step ensured the smallest size, as is further discussed, and this fact influenced the observed characteristics considerably. Sample I suffered from the lowest encapsulation efficiency and drug loading values, which ware probably caused by higher drug losses during preparation. As the particles were of a smaller size, they had a greater surface area on the same weight ratio compared to the other samples, thus providing a larger area for drug leakage²²⁾.

Sphericity, equivalent diameter and morphology

Sphericity and particle size of prepared microparticles were evaluated by optical microscope analysis using specialized software. Images from the microscope are shown in Figure 1. Results of sphericity and equivalent diameter are presented in Table 2. Generally, all samples showed great sphericity as the sphericity factor ranged from 0.892 ± 0.045 to 0.951 ± 0.023 and it is difficult to find any dependency pattern. The particle size ranged from 142.6 ± 25.8 µm to 352.1 ± 80.0 µm. Results clearly show that the larger particles were formed at the stirring rate of 600 rpm. It has been reported before that, among other parameters, the size of microparticles can be affected by stirring speed^{23, 24)}. Increasing stirring speed reduces the mean diameter and yields lower microparticles polydispersity¹²⁾. Also, it could be concluded that at 600 rpm the outer phase volume can influence the size as the samples A and B had a lower size compared to samples C and D; however, because of standard deviation values overlay a multidimensional analysis was needed for the confirmation (see below). Nevertheless, similar relationship was described in Yang’s et al. article based on a greater probability of connecting microparticles with big pores in a smaller volume of water phase^{25, 26)}. In the samples prepared at 1000 rpm a possible dependency of the particle size on the volume of the aqueous phase was not shown.

Because of the significantly smaller size of the microparticles prepared by pre-emulsion step, it was not possible to establish sphericity and equivalent diameter measurement using stereomicroscopic analysis. Therefore, sample I was captured by a Nikon Eclipse optical microscope (Fig. 2). Particle size was determined by laser diffraction, results are shown in Figure 3. The frequency curve is characterized by a number-based particle size distribution. Significant part of the measured microparticles was in the size range of 3 to 30 µm, proving that emulsification via homogenizer yields a significantly lower particle size. This was further confirmed by surface morphology examination using scanning electron microscopy (SEM). SEM revealed that the microparticles of sample I had a spherical shape without visible cracks, the average size of the shown particles varied from 0.5 to 15 µm and the size distribution seemed unequable within this range. SEM also indicated smooth and creaseless (or unfolded) surface with no visible pores on it (Fig. 2).

Multivariate data analysis

To support the initial finding about possible complex co-influences, multivariate data analysis was performed. A data set of samples A to H with process variables representing the preparation method used, stirrer speed and aqueous phase volume were evaluated by factor analysis to find the relationships between the formulation/process variables and qualitative parameters of the prepared microparticles. Sphericity was not included in the factor analysis due to low inter-sample variability and the analysis did not evaluate sample I. The first two factors were used for the calculation, which together explain 79.3% of the total variability in the data. In the factor load graph (Fig. 4), the first factor is explained by the emulsification method used and the second factor is explained by the stirrer speed. In the graph of factor loads, a correlation of EE and the yield with the direct emulsification method and a negative correlation with the NE-1000 dispenser method can be observed. DL is correlated with the speed, while the particle diameter decreases with the speed. (However, changes in DL with changing speed are in the range of about 0.5% and are not of great practical importance for the microparticle preparation). According to Figure 4, the influence of the aqueous phase volume is insignificant. (This result probably comes from the fact that the qualitative parameters (EE, yield, ...) do not change in direct proportion to the volume, but the change in qualitative parameters with the change of volume is conditioned by certain stirring speed). Multivariate data analysis output further confirms the hypothesis about complex phenomenon with existing co-influences.

Drug release behavior

The dissolution profiles of all samples have shown prolonged drug release and in addition, samples A–H had only a mild burst effect (5–10%, Fig. 5, Table 3). The faster initial release was observed in samples E–H when compared to their corresponding samples from A–D group²⁶⁾. The slowest release of ibuprofen was observed in sample A with lower stirring speed and higher aqueous volume, where only 27.40% of the substance was released after 72 hours. The fastest drug release was observed in sample I with the smallest particle size. During the first 30 minutes, 58.41% of drug was released. From the preparation perspective it can be concluded that the samples prepared with a higher stirring speed or with a pre-emulsification step had faster drug release, resulting from the smaller particle size²⁷⁾. From the results it is also possible to find an eventual dependency on emulsification method and outer aqueous phase volume in samples prepared at 600 rpm. In samples prepared at 1000 rpm these differences are wiped, suggesting a stirring speed as an important co-factor.

Based on the ​​obtained dissolution profiles, similarity factors were calculated (Table 4). For sample I, the values ​​of f₂ ranged from 11.66 to 20.57, meaning that the dissolution profile shows no similarity, clearly apart of other samples. Generally, the low f₂ values suggest that the observed parameters could have a significant effect on the drug release, meaning the preparation with pre-emulsification step yields a unique profile. Majority of similarities found was divided mainly among two stirring speed groups, A–D and E–H, respectively. At 600 rpm it is observable that just one change of preparation parameters still yields a similar profile. At 1000 rpm this trend is not so evident, further confirming stirring speed as a major variable during microparticle preparation.

The most important information obtained from the dissolution profile using the kinetics equation ​​is the expression of the drug release kinetics and mechanisms. The drug release from microparticles can occur by several mechanisms including physical erosion, hydrolysis of the polymer, diffusion through the surface of the pores or through the polymer, ion exchange mechanism or by and a combination of those mechanisms^{28, 29)}. A determination coefficient R² of a given equation served as a parameter for the correlation evaluation with the applied existing model. The resulting values ​​are presented in Table 5. Neither sample showed any correlation with the first order kinetic model (perhaps except for sample B, which nearly correlates at 0.987), nor the Higuchi, Baker-Lonsdale and Korsmeyer-Peppas release mechanism mathematical models. Samples E to H showed clear correlation with the zero-order kinetic model (0.991–0.994), suggesting at least partial gradual release not depending on the remaining drug concentration. Low correlation of the samples with the Higuchi model suggests that the diffusion is not a main release mechanism and the drug is probably released predominantly by erosion mechanism. General absence of the correlation with Korsmeyer-Peppas model could be explained by segmented, multi-phasic profiles of some samples. Moreover, the results were obtained using a classic linear regression; however, a more correct approach is the use of non-linear regression, as was reported recently³⁰⁾.

Conclusion

A total of 9 ibuprofen-loaded samples was successfully prepared by the O/W solvent evaporation method. As the model drug substance in these microparticles ibuprofen was picked and the PLGA polymer was used as the polymer carrier. Different preparation parameters were applied, which were evaluated in terms of particles morphological properties such as sphericity, equivalent size distribution; and further the yield, encapsulation efficiency and drug loading were monitored. All prepared samples had very high values of sphericity factor. A very negative effect on the yield and encapsulation efficiency had the utilization of a NE-1000 dispenser and an ULTRA-TURRAX which, on the other hand, provided promising reduction of the microparticles size. A change in the microparticle mean size was also evident when increasing the stirring speed from 600 rpm to 1000 rpm. Smaller particles further influenced drug release profiles, causing faster ibuprofen release, compared to corresponding 600  rpm samples, and also raising a significant burst effect of sample I. The 1000 rpm samples correlated with the drug release kinetics of the zero order, the ideal type of kinetics for this type of dosage form. From the point of view of morphology and drug content, the best result was shown in sample A prepared with 600 rpm and 800 ml of the aqueous phase, possessing the highest yield and encapsulation efficiency.

Conflict of interest: none.

Acknowledgement

This research was supported by the Ministry of Education, Youth and Sports of the Czech Republic under the project CZ.02.1.01/0.0/0.0/15_003/0000495.

PharmDr. Jakub Vysloužil, Ph.D.

Department of Pharmaceutical Technology

Faculty of Pharmacy, Masaryk University

Palackého 1, 612 42 Brno, Czech Republic

e-mail: vyslouzilj@pharm.muni.cz