1Departement of Pharmaceutics, Faculty of Pharmacy, Universitas Gadjah Mada, Sekip Utara, D. I. Yogyakarta 55281, Indonesia, 2Department of Dermatology and Venereology, Faculty of Medicine, Public Health and Nursing, Universitas Gadjah Mada, Farmako Sekip Utara, D. I. Yogyakarta 55281, Indonesia
Email: [email protected]
Received: 04 Jun 2019, Revised and Accepted: 04 Jul 2019
Objective: This study aims to optimize the timolol maleate (TM) nanoparticle prepared by ionic gelation method using the factors of pectin (PC), calcium chloride (CC), and chitosan (CS) concentrations with the responses of entrapment efficiency, particle size, and polydispersity index using 23 factorial design.
Methods: TM nanoparticle suspensions were obtained by mixing of PC (0,4-0,6% (w/v)), CC (0,2-0,4% (w/v)), and CS (0,01-0,02% (w/v)) with TM concentration of 0,02% w/v. Each mixture was then tested for entrapment efficiency, particle size, and polydispersity index. The test results were analyzed with 23 factorial design using Design-Expert software in order to determine the optimum formula.
Results: The optimization study showed that all of the factors influenced the responses significantly (p<0.05) based on the analysis of variance (ANOVA) of the suggested models. The R2value and the adequate precision value of the three models were more than 0.7 and 4, respectively. The difference between Adjusted R-Squared and Predicted R-Squared value were less than 0.200. The optimum condition of TM nanoparticle was suggested at the desirability value of 0.839 with the concentration of PC, CC, and CS of 0,4% (w/v), 0,2% (w/v), and 0,01% (w/v), respectively. The entrapment efficiency, particle size, and polydispersity index of the optimum condition were 24.791±2.84%, 274.867±14.45 nm, and 0.634±0.066, respectively.
Conclusion: The 23factorial design has been proved as the suitable method to determine the optimum condition that yields the good results of the entrapment efficiency, particle size, and polydispersity index of the TM-loaded nanoparticle prepared by ionic gelation method.
Keywords: Nanoparticle, Ionic gelation, Timolol maleate, Pectin, Calcium chloride, Chitosan, Factorial design
© 2019 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/)
Timolol maleate (TM) is a β-blocker agent which has a promising alternative in the treatment of infantile hemangioma (IH), a benign vascular tumor that frequently occurs in infants with the prevalence of approximately 4-10% in the first year of life [1, 2]. Application of topical TM for IH treatment shows better efficacy and less potency of adverse effect. A recent study exhibited a significant decrease of IH surface area treated with topical timolol maleate compared with topical corticosteroid . Nevertheless, topical use of TM is reported to reach systemic circulation, which can produce undesirable effects .
Polymeric nanoparticle can be utilized in order to provide an effective therapy of IH using TM with less potency to be absorbed through the systemic circulation. A polymeric nanoparticle is one of the drug carriers that can be used in the formulation of a topical dosage form. It possesses several abilities such as protecting the drug trapped inside the nano-sized particle from degradation or denaturation and decreasing the potency of systemic absorption by providing sustained release of drug from the dosage form . The nanoparticle can be made by the method of ionic gelation, which is a method that utilized the electrostatic interaction between cationic and anionic polymer that spontaneously from nano-sized particles. This advantageous method only requires simple stirring without the use of organic solvents that mostly possess toxic properties [5-7].
The ionic gelation system can be achieved in the event of ionic interaction between the oppositely charged polymers and cross-linking agents that subsequently able to encapsulate the drug molecules. Pectin (PC) is an anionic polymer extracted from plant cell walls. The combination of PC with chitosan (CS) as the cationic polymer and calcium chloride (CC) as the cross-linker ion yields the nanoparticulate system with the nano-sized particles and drug encapsulated inside the particles. The combination has also prolonged the release of drug from nanoparticle so that the frequency of drug administration can be reduced .
The concentration of polymers and cross-linking agents used in the ionic gelation greatly affects the nanoparticle formed. An optimum value of several parameters such as entrapment efficiency, particle size, and polydispersity index can be achieved by the appropriate concentration of PC, CC, and CS. The use of factorial design as one of the methods of nanoparticle formula optimization can be conducted in order to determine the factors affecting the experiment results and to observe the level of factors that generate the desirable responses .
This study aims to observe the effect of PC, CC, and CS concentrations on the entrapment efficiency, particle size, and polydispersity index of the timolol maleate-loaded nanoparticle manufactured with ionic gelation method. Formula optimization is conducted with 23 factorial design using three replicates to analyze the effects of the three selected factors.
Timolol maleate (TM) was purchased from Octagon Chemicals Limited (China), pectin (PC) was purchased from Sigma Aldrich (Darmstadt, Germany), calcium chloride (CC) was purchased from Merck (Germany), chitosan (CS) was purchased from local company Chimultiguna (Cirebon, Indonesia), glacial acetic acid, hydrochloric acid, and sodium hydroxide were purchased from Merck (New Jersey, USA).
Instrumentation and software
The entrapment efficiency of TM in nanoparticle was analyzed using UV-Vis spectrophotometer (Thermo Scientific Genesys 10S UV). The particle size, polydispersity index, and zeta potential of nanoparticle were analyzed by the method of dynamic light scattering (DLS) using Zetasizer Nano ZS (Malvern, UK). The optimization study was analyzed using Design-Expert software (Stat-Ease Inc., Minneapolis, MN, USA). The verification study of optimum nanoparticle formula was analyzed statistically using R free software.
Formula optimization using a factorial design
Formula optimization was conducted with 23 factorial design using Design-Expert software. Concentration of PC (X1) CC (X2), and CS (X3) were used as factors, whereas entrapment efficiency (Y1), particle size (Y2), and polydispersity index (Y3) were used as responses.
Prediction and verification of optimum nanoparticle formula
The prediction of the optimal condition of nanoparticle preparation was conducted by determining the priority value of each response (Y) in Design-Expert software. The verification of optimal condition was conducted by comparing the three replicates of experiment data with the predicted value provided by the software. Data were analyzed with one-sample t-test using R free software at 95 % of confidence interval (p>0.05) and using Design-Expert software by observing the 95% confidence interval (CI) and 95% prediction interval (PI).
Preparation of TM-loaded PC-CS nanoparticle
PC solutions (0,4-0,6% (w/v)) and CC solutions (0,2-0,4% (w/v)) were obtained by dissolving certain amounts of PC and CC separately in demineralized water for 2 h with constant stirring. CS solutions (0,01-0,02% (w/v)) were obtained by dissolving certain amounts of CS with 1% of acetic acid (v/v) solution in overnight constant stirring. TM solution (0,02% (w/v)) was obtained by dissolving certain amounts of TM in demineralized water with constant stirring. Furthermore, 2 ml of PC solution was added with 0,1 ml of TM solution with constant stirring for an hour. 1,9 ml of CC solution was then added to the mixture solution of PC and TM and stirred for 30 min to complete the pre-gelation process. Finally, 1 ml of CS was added to the pre-gel solution and stirred for 30 min .
Nanoparticle sample was centrifuged with a speed of 15000 rpm for 30 min at 4 °C to separate the unentrapped drug with the nanoparticle precipitate. The supernatant was then diluted with HCl 0,1 M and analyzed for TM content using UV-Vis spectrophotometer at 296 nm. HCl 0,1 M was used as the blank. The test was conducted in three replicates for each sample. The entrapment efficiency (EE) value was calculated using following equation:
Table 1: The condition of each experiment in the factorial design and their responses (n=24)
|X1: Concentration of PC (%)||X2: Concentration of CC (%)||X3: Concentration of CS (%)||Y1: Entrapment efficiency (%)||Y2: Particle size (nm)||Y3: Polydispersity index|
Particle size and polydispersity index
The test was conducted using a particle size analyzer instrument of Zetasizer Nano ZS (Malvern, UK). A number of timolol maleate nanoparticle sample was put into the cuvette and analyzed in three replicates .
Optimization studies using a factorial design
Factorial design is one of the methods commonly used in optimization study using the design of the experimental (DoE) approach. The method can be purposed to probe factors that influence the outcome of the experiment. Moreover, the levels of factors that generate a better response in the experiment can also be observed using factorial design . Variables that can be controlled by the investigator can be used as factors in factorial design. In this study, the concentration of PC, CC, and CS were used as factors due to the fact that the process of polymeric nanoparticle preparation is influenced by the concentration of polymers and cross-linking agents [12–14]. The results of the optimization study that was conducted using factorial were presented as 24 experiment conditions (table 1).
Table 2: The regression equations of all responses used in the design of experimental
|Entrapment Efficiency (Y1)||Y1 = 64.199–43.854X1–141.348X2–447.975X3+175.470X1X2–1170.583X1X3+1949.416X2X3|
|Particle Size (Y2)||Y2 =-2457.067+5292X1+9423.833X2+182940X3-18013.333X1X2-335217X1X3–613150X2X3+1.127670X1X2X3|
|Polydispersity Index (Y3)||Y3 = 1.205–0.826X1+0.206X2–44.550X3–1.200X1X2+76.833X1X3+5.333X2X3|
The analysis of variance (ANOVA) of the suggested models generated by Design-Expert software indicated that three factors influenced the responses used that showed by a significant model (p<0.05) of the three responses. The equation model of the three responses (Y1, Y2, and Y3) was statistically able to predict the optimum condition of nanoparticle preparation. The regression equation (table 2) showed the contribution of each factor on increase and decrease of the evaluated responses marked by positive and negative signs, respectively.
Table 3: The results of the acceptance criteria of statistical parameters of the experimental design
|Response||SD||CV||R2||Adjusted R2||Predicted R2||Adequate Precision|
|Y1: Entrapment Efficiency||2.55||19.46||0.7774||0.6989||0.5564||8.489|
|Y2: Particle Size||49.39||15.33||0.7546||0.6472||0.4479||7.176|
|Y3: Polydispersity Index||0.041||6.77||0.7251||0.6281||0.4522||7.273|
Based on the factorial method parameter on DoE (table 3), R2 measure the effect of factor to the response in the equation, the higher R2 (more than 0.7) more preferable (Y1=0.7774; Y2=0.7546; Y3=0.7251). The difference between Adjusted R-Squared and Predicted R-Squared value must be less than 0.200 (Y1=0.1425; Y2=0.1993; Y3=0.1759) and a good correlation between the actual value and prediction value of the three responses were presented by predicted vs actual plots (fig. 1A, 1B, 1C). Adequate Precision showed the signal to noise ratio and must be more than 4 (Y1=8.489; Y2=7.176; Y3=7.273); thus all response indicates the adequate equation to describe the relationship between factors to responses [15–17].
The entrapment efficiency (Y1) of the 24 experiments was varied from 13.417% to 29.797% (table 1). Based on the ANOVA analysis, the model of entrapment efficiency response (Y1) was significant with the p-value of<0.0001 (p<0.05). The factor of CC concentration (X2) and interaction between factors of concentration of PC and concentration of CC (X1X2) possessed a significant effect with the p-value of 0.0003 and 0.0036, respectively.
The factor of X2 possessed a negative effect, whilst the interaction of X1X2 possessed a positive effect on the response of entrapment efficiency. The negative effect of CC concentration indicates that the decrease of entrapment efficiency with increasing the concentration of CC (fig. 2A). The excess Ca2+ion resulted in a high concentration of CC decrease the strength of gel surface in the pre-gelation process. There is a possibility that the high concentration of CC (0,4%) used in this study generates saturation of calcium ions in the nanoparticle preparation lead to a decrease of entrapment efficiency .
The entrapment efficiency is also affected by the interaction between the concentration of PC and CC. It can be seen in the table 1 that in the most of experiments with the low concentration of CC (0.2% (w/v)), the entrapment efficiency values were lower when the high concentration of PC (0.6% (w/v)) is used than the low concentration of PC (0.4% (w/v)). The decrease of entrapment efficiency can be occurred due to the increase of the solution viscosity generated by a high concentration of PC. The higher solution viscosity subsequently causes an inadequate interaction between the PC and the cross-linking agent .
The results of experiments showed that the response of particle size (Y2) of the 24 runs was varied between 177.8 nm and 495.3 nm. The ANOVA analysis indicated that the suggested model was significant with the p-value of 0.0006 (p<0.05). All of the factors (X1, X2, X3) and interaction of factors (X1X2, X1X3, X2X3) were significant with the p-values of<0.0001 (p<0.05).
Based on the regression equation of particle size response, the factors of PC (X1) and CC (X2) concentrations showed positive effects, which indicate that the higher concentration of PC and CC, the bigger particle size will be yielded. The effects were also illustrated by the 3D surface of particle size (fig. 2B). This finding is in accordance with the previous study describing that the small particles are formed when the low concentration of PC and divalent cation . The concentration of CS (X3) was also indicated a positive effect on the particle size of TM nanoparticle. The higher chitosan concentration increases the viscosity of nanoparticle preparation and subsequently slower the gelation process which yields the higher particle size . The interaction between the factor of X1 and X3 showed negative effect indicates that small and stable colloid particles will be formed in the adequate interaction between PC’s carboxylate groups and oppositely charged amino groups of CS .
The results of 24 runs showed that the values of polydispersity index (Y3) were between 0.484 and 0.698. The ANOVA analysis indicated that the suggested model was significant with the p-value of 0.0005 (p<0.05). The factors of PC concentration (X1), CS concentration (X3), and interaction between of PC and CS (X1X3) were significant with the p-values of 0.0366, 0.0003, and 0.0003, respectively (p<0.05). Based on the produced regression equation (table 2), the negative effect of PC concentration (X1) and CS concentration (X3) indicate that the higher concentration of both polymers decreases the polydispersity index value, which is also illustrated by the 3D surface of polydispersity index (fig. 2C). This finding is in accordance with another study which finds that the increase of polymer concentrations promotes a decrease in polydispersity index value . Nevertheless, the interaction between the concentration of PC and CS (X1X3) yielded a positive effect. This interaction lead to an increase of polydispersity index values due to the aggregations in the nanoparticle sample that occurs as the result of the increase of ionic interaction between the PC’s carboxylate groups and CS’s amino groups .
Fig. 1: Actual and predicted value correlation of entrapment efficiency (A), particle size (B), and polydispersity index (C)
Fig. 2: Interaction of factors of X1, X2, and X3 on entrapment efficiency (A), particle size (B), and polydispersity index (C) illustrated in the form of 3D surface diagrams
The prediction and verification of the optimum condition of nanoparticle
The selection of the optimum condition of nanoparticle preparation can be done by considering the desirability value of the solutions provided by the Design-Expert software. Desirability is an indicator that describes the closeness between the prediction of response values and the observed response values. The desirability is presented at the value between 0 and 1, where the 0 value indicates an undesirable value and 1 indicates the desirable value or the ideal response. The desirability considered as good if it has the value close to 1 .
Based on the determination of goal and importance of each response, the highest desirability value of 0.839 was obtained and considered as the optimum condition of nanoparticle preparation with the concentration of PC (X1) 0.4% (w/v), the concentration of CC (X2) 0.2% (w/v), and the concentration of chitosan (X3) 0.01% (w/v). The prediction values of each response were 27.162% for entrapment efficiency (Y1), 267.8 nm for particle size (Y2), and 0.692 for polydispersity index (Y3). The contour plot of the desirability value of 0.839 can be seen in fig. 3.
The results of the verification study of the three responses were provided in table 4. The entrapment efficiency value of the verification study was 24.791±2.84 %, which was theoretically in the range of 95% CI verification range (24.765 to 29.559%) and 95% PI range (22.117 to 32.707%). The result of the statistical test of entrapment efficiency using one-sample t-test showed that the p-value of the test was 0.285 (>0.05), which indicated that the entrapment efficiency value of the verification study was not significantly different from the predicted value. The quite low entrapment efficiency value of 24.791±2.84 % may be attributed to the physichochemical characteristic of the drug used in this study. TM is a hydrophilic drug that has a possibility to escape to the external environment that overwhelmed by hydrophilic substances. This phenomenon causes only a few interactions occurred between the drug and the polymers and lead to poor entrapment efficiency .
Fig. 3: Desirability value of optimal condition of nanoparticle preparation
Table 4: The prediction and observation value of optimum condition with the statistical range of verification
|Response||Prediction||Observation*||One Sample t-test**||95% CI||95% PI|
|Y1: Entrapment Efficiency (%)||27.162||24.791±2.84||0.285||24.765||29.559|
|Y2: Particle Size (nm)||267.8||274.867±14.45||0.486||218.013||317.587|
|Y3: Polydispersity Index||0.692||0.634±0.066||0.269||0.653||0.730|
*results expressed in the mean of n=3±standard deviation (SD), **(p>0.05)
The result of verification study of particle size was also indicated that the values observed in the verification were not significantly different with the predicted values and in the range of 95% CI and 95% PI, while the polydispersity index verification study showed that the observed value was in the range of 95% PI. Based on the one-sample t-test results of the observed value of particle size and polydispersity index, the p-values of particle size and polydispersity index were respectively 0.486 and 0.269, which indicated that the observed values were not significantly different from the predicted values.
The particle size of 274.867±14.45 nm is considered an acceptable size for the intended use of the nanoparticle. The TM nanoparticle is intended to be administrated on the skin surface as the therapy of infantile hemangioma. The maximal depth penetration of particle size in the range between 250 and 500 nm can only reach the stratum corneum and unable to penetrate deeper through the skin . The polydispersity index value of 0.634±0.066 (less than 0.7) indicates that the nanoparticle sample is able to be analyzed using dynamic light scattering method due to a narrow particle size distribution. The optimum formula possesses monodisperse particle distribution, which indicates size uniformity of the particles inside of the nanoparticle sample with small variation between individual particle [28, 29].
The optimization study to observe the optimum condition of TM nanoparticle preparation using ionic gelation technique was successfully applied using the method of 23 factorial design. The optimum condition which produced good results of the response of entrapment efficiency, particle size, and polydispersity index was found with the concentration of PC, CC, and CS of 0.4% (w/v), 0.2% (w/v), and 0.01% (w/v), respectively.
The authors would like to thank the research grant of Penelitian Dasar Unggulan Perguruan Tinggi (PDUPT) 2018 Kemenristekdikti-Universitas Gadjah Mada for funding this research.
All of the authors listed in the manuscript has contributed equally
The authors declare there is no conflict of interest is associated with this work
Danarti R, Ariwibowo L, Radiono S, Budiyanto A. Topical timolol maleate 0.5% for infantile hemangioma: its effectiveness compared to ultrapotent topical corticosteroids-a single-center experience of 278 cases. Dermatology (Basel) 2016;232:566–71.
Zheng JW, Zhang L, Zhou Q, Mai HM, Wang YA, Fan XD, et al. A practical guide to treatment of infantile hemangiomas of the head and neck. Int J Clin Exp Med 2013;6:851–60.
Frommelt P, Juern A, Siegel D, Holland K, Seefeldt M, Yu J, et al. Adverse events in young and preterm infants receiving topical timolol for infantile hemangioma. Pediatr Dermatol 2016; 33:405–14.
Zhang Z, Tsai PC, Ramezanli T, Michniak Kohn BB. Polymeric nanoparticles-based topical delivery systems for the treatment of dermatological diseases. Wiley Interdiscip Rev: Nanomed Nanobiotechnol 2013;5:205–18.
Kunjachan S, Jose S, Lammers T. Understanding the mechanism of ionic gelation for the synthesis of chitosan nanoparticles using qualitative techniques. Asian J Pharm 2014;4:148-53.
Grabnar PA, Kristl J. The manufacturing techniques of drug-loaded polymeric nanoparticles from preformed polymers. J Microencapsul 2011;28:323-35.
Nagpal K, Singh SK, Mishra DN. Chitosan nanoparticles: a promising system in novel drug delivery. Chem Pharm Bull 2010;58:1423-30.
Grabnar PA, Kristl J. Physicochemical characterization of protein-loaded pectin-chitosan nanoparticles prepared by polyelectrolyte complexation. Pharmazie 2010;65:851-2.
de Pinho Neves AL, Milioli CC, Müller L, Riella HG, Kuhnen NC, Stulzer HK. Factorial design as a tool in chitosan nanoparticles development by ionic gelation technique. Colloids Surf A 2014;445:34-9.
Ilka R, Mohseni M, Kianirad M, Naseripour M, Ashtari K, Mehravi B. Nanogel-based natural polymers as smart carriers for the controlled delivery of timolol maleate through the cornea for glaucoma. Int J Biol Macromol 2018;109:955-62.
Siafaka PI, Titopoulou A, Koukaras EN, Kostoglou M, Koutris E, Karavas E, et al. Chitosan derivatives as effective nanocarriers for ocular release of timolol drug. Int J Pharm 2015;495:249–64.
Shaw R, Festing MFW, Peers I, Furlong L. Use of factorial designs to optimize animal experiments and reduce animal use. ILAR J 2002;43:223-32.
Bhoskar M, Patil P. Development and evaluation of paclitaxel loaded nanoparticles using 24 factorial design. Int J Curr Pharm Res 2015;7:64-72.
Jain A, Jain SK. Formulation and optimization of temozolomide nanoparticles by 3 factors 2 level factorial design. Biomatter 2013;3:1-13.
Pham DT, Saelim N, Tiyaboonchai W. Design of experiments model for the optimization of silk fibroin based nanoparticle. Int J Appl Pharm 2018;10:195-201.
Shiyan S, Hertiani T, Martien R, Nugroho AK. Optimization of a novel kinetic-assisted infundation of white tea (Camellia sinensis) using central composite design. Int J Appl Pharm 2018;10:259-67.
Setyawan EI, Setyowati EP, Rohman A, Nugroho AK. Central composite design for optimizing the extraction of EGCG from green tea leaf (Camellia sinensis L.). Int J Appl Pharm 2018;10:211–6.
Pawar AP, Gadhe AR, Venkatachalam P, Sher P, Mahadik KR. Effect of core and surface cross-linking on the entrapment of metronidazole in pectin beads. Acta Pharm 2008;58:78–85.
Sharma R, Ahuja M, Kaur H. Thiolated pectin nanoparticles: preparation, characterization and ex vivo corneal permeation study. Carbohydr Polym 2012;87:1606–10.
Opanasopit P, Apirakaramwong A, Ngawhirunpat T, Rojanarata T, Ruktanonchai U. Development and characterization of pectinate micro/nanoparticles for gene delivery. AAPS PharmSciTech 2008;9:67–74.
Zaki SSO, Ibrahim MN, Katas H. Particle size affects concentration-dependent cytotoxicity of chitosan nanoparticles towards mouse hematopoietic stem cells. J Nanotechnol 2015. http://dx.doi.org/10.1155/2015/919658
Maciel VBV, Yoshida CMP, Pereira SMSS, Goycoolea FM, Franco TT. Electrostatic self-assembled chitosan-pectin nano-and microparticles for insulin delivery. Molecules 2017;22:1-21.
Melo NFS, Campos EVR, Paula ED, Rosa AH, Fraceto LF. Factorial design and characterization studies for articaine hydrochloride loaded alginate/chitosan nanoparticles. J Colloid Sci Biotechnol 2013;2:146–52.
Masarudin MJ, Cutts SM, Evison BJ, Phillips DR, Pigram PJ. Factors determining the stability, size distribution, and cellular accumulation of small, monodisperse chitosan nanoparticles as candidate vectors for anticancer drug delivery: application to the passive encapsulation of 14C-doxorubicin. Nanotechnol Sci Appl 2015;8:67-80.
Candioti LV, De Zan MM, Camara MS, Goicoechea HC. Experimental design and multiple response optimization. Using the desirability function in analytical methods development. Talanta 2014;124:123-38.
Abdalla KA, Kamoun EA, Maghraby GME. Optimization of the entrapment efficiency and release of ambroxol hydrochloride alginate beads. J Appl Pharm Sci 2015;5:13-9.
Poland CA, Read SAK, Varet J, Carse G, Christensen FM, Hankin SM. Dermal absorption of nanomaterials part of the ”Better control of nano” initiative 2012-2015. Denmark: The Danish Environmental Protection Agency; 2013.
Danaei M, Dehghankhold M, Ataei S, Hasanzadeh Davarani F, Javanmard R, Dokhani A, et al. Impact of particle size and polydispersity index on the clinical applications of lipidic nanocarrier systems. Pharmaceutics 2018;10:57.
Zayed GM, El-feky GS. Growth factor loaded functionalized gold nanoparticles as a potential targeted treatment for acute renal failure. Int J Appl Pharm 2019;11:174–85.