Enhancing the dissolution of phenylbutazone using Syloid®based mesoporous silicas for oral equine applications

Lur J.Wters,John P.Hnrhn,Joseph M.Toin,Ctherine V.Finch,Greth M.B.Prkes,Shmsuddeen A.Ahmd,Frj Mohmmd,Mri Sleem

aSchool of Applied Sciences,University of Hudders field,Queensgate,Hudders field HD1 3DH,UK

bGlantreo Ltd,ERI Building,Lee Road,Cork City T23 XE10,Ireland

1.Introduction

Mesoporous silica has bee n shown to exhibit a great potential to aid in the formulation of pharmaceutical compounds with poor aqueous solubility,as reviewed by Choudhari et al.[1].As a drug carrier system,mesoporous silica can accommodate drugs that have been introduced through organic solvent immersion,incipient wetness impregnation or melted in[2].Specific advantages of using excipients such as mesoporous silicas are their nanoporous structures,high surface areas,clinical safety and large pore volumes[3].Current opinion is that substantial progress has been made in recent years in the characterisation and development of mesoporous drug delivery systems although more work is needed regarding dissolution enhancement potential and related physicochemical properties[4].There are several reasons for this need to continue exploring the possible use of mesoporous silica such as practical considerations such as manufacturability to large scale quantities(e.g.tonne)and regulation,as well as physicochemical considerations such as the possibility of re-adsorption onto the silica surface[4].Adsorption of small drug particles on the surface of large excipients has been a successful strategy for low-dose drugs,poorly water soluble drugs,targeted drug release[5],sustained drug delivery[6]and stability enhancement.This is mainly a result of improving the dissolution pro file by increasing drug surface area or transformation of the drug from a crystalline to amorphous form[7],and its ability to be retained within the silica pores[8].In many cases,the method of formulation can be critical in defining the properties of the resultant formulation.For example,silica-based drug delivery vehicleshavebeen investigated to avoid hydrolisation of the active compound using supercritical CO2[9,10],a formulation method known for its high drug-loading ability[11]amongst other advantages[12].Several other formulation methods have also been attempted,for example,to create liquid(also known as liquisolid)formulations[13]and pediatric(solvent free)formulations[14].The work within our group that has previously confirmed the application of microwave irradiation for mesoporous silicas[15].Furthermore,there is clearly an interest in developing mesoporous silica formulations as evidenced by recent work to predict in vivo performance,for example,using in silico techniques[16],to overcome multidrug resistance[17]as well as to ameliorate toxic side effects[18].

One particular category of mesoporous silicas where only very limited studies have been conducted to date is regarding Syloid®silica based formulations.These forms of silica have a highly developed network of mesopores that provide access to the large surface area,i.e.a combination of a high adsorption capacity,along with a desirable pore size and surface morphology.For these reasons,these silicas tend to be used to improve the flow properties of pharmaceuticals where liquid ingredients can be converted into free- flowing powders.Although these properties are beneficial,their suitability to enhance dissolution has only briefly been considered(by publication)for two forms of Syloid®silica(244 and AL-1)with two model drugs,namely,indomethacin[19]and itraconazole[20].Interestingly,for both compounds,an enhancement in the rate and extent of dissolution was observed in both studies.Yet surprisingly,other forms of Syloid®silica have not yet been considered even though they may provide a plethora of advantages for drug-loading formulations.

One specific drug renowned for having poor aqueous solubility∶0.05 mg/mL[21],and therefore problematic dissolution with potentially low bioavailability,is phenylbutazone.This particular compound is commonly used in equine environments as a nonsteroidal anti-inflammatory drug(NSAID)[22,23],often prescribed for pain control[24,25].Although drug solubility is significantly greater in ethanol and 1-octanol[26],the low level of aqueous solubility results in complications for formulators.One study has successfully enhanced the dissolution through the creation of a solid dispersion with polyethylene glycol 8000[27]and another with SBA-15[28],yet there is still a clear need for developing alternative formulations that can achieve an even greater enhancement in release of the active compound.Phenylbutazone is an excellent candidate for exploring the potential to enhance solubility through formulation with a dissolution-limiting low solubility yet incredibly significant usage within the equine community.This is because many of the present formulations available on the market tend to be unfavourable with issues surrounding drug delivery and poor palatability[29].Thus ways to enhance phenylbutazone-based formulations are highly desirable.

This work investigates the suitability of using three types of Syloid®silica based excipients to quantify their potential to enhance the rate of dissolution of phenylbutazone and determine the causes of any enhancements observed.

2.Materials and methods

2.1.Materials

Phenylbutazone,potassium phosphate dibasic,and potassium phosphate monobasic(all≥99%)were purchased from Sigma Aldrich(Dorset,UK)and used as received.Syloid®silicas(AL-1 FP,XDP 3050 and XDP 3150)were kindly donated by Glantreo Ltd,Cork,Ireland and W.R.Grace&Co,Maryland,USA.Table 1 provides a summary of the physicochemical properties of the Syloid®silicas,and the data presented were determined using nitrogen gas sorption isotherms.These were measured at 77K using a Micromeritics TriStar II surface area analyzer(Micromeritics,Norcross,GA,USA).Samples were pre-treated by heating at 200°C undernitrogen for 12 h.The surface area was measured using the Brunauer-Emmett-Teller(BET)method.The pore volume and pore diameter data was calculated using the Barrett,Joyner and Halenda(BJH)method[2].Specific surface areas were calculated from the measured relative pressure in the range of P/P0=0.01 to P/P0=0.3.Mesoporous volumes were estimated from the volume of nitrogen adsorbed after the micropores have been filled until after condensation into the mesopores was complete.Of particular interest is the range of surface areas and pore volumes exhibited by the three Syloid®silicas as based on previous research,such properties may influence dissolution.For example,pore size has been known to effect drug release pro files for other mesoporous systems[30].

2.2.Methods

2.2.1.Formulation methods

200mg of Syloid®silica XDP 3050 was placed in a beaker whereupon 40mL of deionised water was gradually added,followed by 200mg of phenylbutazone to achieve a total drug and silica mass of 400mg.Over a period of 60min the solution was stirred and heated to a maximum of 90°C,cooled to room temperature,vacuum filtered and dried overnight at 60°C,and then sieved to remove agglomerates larger than 250μm.This process was repeated in triplicate and then with the replacement of XDP 3050 with XDP 3150 and AL-1 FP to produce a total of three unique drug-Syloid®silica formulations.A series of variable ratios of drug∶Syloid drug∶Syloid®silica formulations were also formulated but based on dissolution pro file data(not shown),no significant differences in release pro files were observed between the formulations;thus this paper only presents formulations at a 1∶1 ratio.A f i nal formulation was produced that involved phenylbutazone undergoing the formulation process(but without the presence of Syloid®silica)to determine if it was the processing that affected dissolution or the presence of each Syloid®silica itself.Water was used as a ‘carrier’to disperse the drug within the mixture,rather than dissolving the drug with organic solvent,followed by heating to help achieve maximum dispersion within the mixture.

2.2.2.Characterisation methods

Powder X-ray diffraction(XRD)data were collected on a Bruker D2-Phaser equipped with a Cu Kα1radiation source at 30 kV and 10 mA current.Particle size distribution of the formulated products was analysed using a Malvern Mastersizer 2000(Worcestershire,UK)using 5–10 mg of powder per sample with one drop of surfactant(IGEPAL®CA-630)at a stirring speed of 2000rpm.Triplicate data was subsequently analysed using Mastersizer 2000 software(V5.61).Drug loading was verified to be 100%in all formulated samples by UV analysis of the filtrates(λ=282 nm)with no residual drug detected(＜1%),thus confirming all of the drug remained within the formulation(rather than washed away with the filtrate during the formulation process).For stability conf i rmation the infrared spectrum for the pure samples and their formulations was recorded using a Nicolet-380 Fourier Transform Infrared spectrometer(FT-IR)with an ATR crystal.Powder samples were placed directly onto the diamond crystal and the anvil was lowered to ensure that sample was in full contact with the diamond.Each spectrum was obtainedintherange of500–4000 cm-1with 2 cm-1resolution.In this study,the morphology of the prepared samples was characterised using scanning electron microscopy(SEM)(JEOL JSM-6060LV,Japan)with gold-plating using a sputter coater(SC7620)prior to imaging.

2.2.3.In vitro phenylbutazone release

Dissolution pro files were determined using a USP Type II(paddle method)PharmaTest DT70 system,with manual sampling for a period of 30 min.Formulated samples with a total drug content of 22.5 mg were placed in 900 mL of pH 7.0 phosphate buffer,stirred at 75 rpm and maintained at 37.0 ± 0.5°C,maintaining sink conditions throughout the duration of the experiment.Filtered samples were removed every 5 min,replaced with phosphate buffer,and analysed using UV spectroscopy(Cary 60,Agilent)set at a wavelength(λ)of 282 nm with conversion to percentage drug release using a standard calibration plot.Samples were analysed in triplicate to determine mean drug release percentages and associated error limits.

3.Results and discussion

3.1.Characterisation of formulations

XRD patterns for samples of the three mesoporous silicas both with and without the presence of phenylbutazone are shown in Fig.1.Previous XRD studies using naproxen noted that the diminishment of peak intensities confirmed that the drug had loaded into channels of a mesoporous material[31],resulting in an amorphous formulation with an absence of characteristic peaks[15].A similar result was observed in this work whereby the purely crystalline phenylbutazone that could be seen in Fig.1A(and after processing in Fig.1B)was converted to the amorphous form following formulation with the three Syloid®silicas(Figs.1C-E).As discussed earlier,from analysing the filtrates and confirming all of the drug had remained within the formulation,the absence of peaks cannot be explained by a reduced concentration of drug and can only be explained by a transformation to the amorphous form.

Particle size analysis confirmed that phenylbutazone(prior to formulation)exhibited an average particle sizes of 65–70μm.The sizes of the three Syloid®silicas prior to formulation are presented in Table 1 and were confirmed in this study to have average values of 10,50 and 110μm for AL-1 FP,XDP 3050 and XDP 3150,respectively.These three Syloid®silicas display an interesting range of particle sizes prior to formulation.Yet their subsequent dissolution pro files may actually be more dependent upon their size after formulation(through the formation of aggregates);therefore,it is this parameter that is of interest in this work.Firstly,AL-1 FP displayed an increase in the average particle size and a slightly broader distribution of sizes with the majority of particles between 5 and 100 μm after formulation.Secondly,a similar result was seen for XDP 3050 with the majority of particles between 40 and 100μm.An explanation for this increase in size and diversity of sizes for both Syloid®silicas is most likely a consequence of particle agglomeration as a result of drug incorporation and/or processing effects.Thirdly,Syloid®silica 3150 did not exhibit any significant increase in average particle size following formulation although there was an increase in the polydispersity of particle size.Again,this indicates agglomeration may have occurred to a limited extent but not to the same degree as that seen for the other Syloid®silicas.

Fig.1.XRD patterns for(a)phenylbutazone,(b)processed phenylbutazone,(c)phenylbutazone and AL-1 FP,(d)phenylbutazone and XDP 3150,and(e)phenylbutazone and XDP 3050.

Fig.2.FT-IR analysis of(a)phenylbutazone,(b)AL-1 FP,(c)XDP 3150,and(d)XDP 3050 silicas prior to formulation.

FT-IR spectroscopy was used to monitor the presence of phenylbutazone and determine interactions with the three silicas(Figs.2 and 3).Analysis of spectra for phenylbutazone showed the expected absorption bands at wavenumbers(with corresponding functional groups)of 754 and 1483cm-1(C-H),1270cm-1(C-N)and 1720cm-1(C=O).Analysis of the spectra for phenylbutazone subjected to the processing method did not reveal any changes in the specific absorption bands for the drug,suggesting a lack of degradation as a result of the formulation process.The three Syloid®silicas were analysed using FT-IR spectroscopy and all displayed the expected intense Si-O absorption band at 1060–1070cm-1[32].For the three phenylbutazone-silica formulated products,the results indicated a significant disappearance of the drug,mainly displaying spectra corresponding to just each type of silica present.Furthermore,the spectra did not display any obvious additional peaks,thus indicating there had been no significant changes in the chemical structure or drug-silica interactions.

Fig.3.FT-IR analysis of(a)phenylbutazone,(b)processed phenylbutazone,(c)phenylbutazone and AL-1 FP,(d)phenylbutazone and XDP 3150,and(e)phenylbutazone and XDP 3050 silicas after formulation.

Fig.4.Scanning electron micrographs for particles of(A)phenylbutazone,(B)AL-1 FP,(C)XDP 3150,and(D)XDP 3050 silicas prior to formulation.

Surface morphologies of the pure phenylbutazone and the three Syloid®silicas prior to formulation,processed phenylbutazone,and Syloid®silica-based formulations–XDP 3050,XDP 3150 and AL-1 FP are presented in Figs.4 and 5.The drug's crystalline state,along with the disordered irregular shapes of AL-1 FP,XDP 3150,and XDP 3050 silicas was evident by SEM(Fig.4).The SEM image confirmed the insignificant effect of processed phenylbutazone as the drug retained a crystalline structure.However,there was a uniform distribution of phenylbutazone on the surface of AL-1 FP due to a larger surface area,smaller pore volume and pore diameter.For Syloid®XDP 3150 and XDP 3050 based formulations,there was an even distribution of the former particles with phenylbutazone particles reduced in size while for the latter,more of the drug was con fined in the pores and on the surface,visible in the SEM images(Fig.5).

3.2.In vitro phenylbutazone release

Dissolution pro files of phenylbutazone loaded Syloid®silicas were investigated for a period of 30 min in pH 7.0 phosphate buffer.As can be seen in Fig.6,pure phenylbutazone that had not undergone the formulation process exhibited 7.2%(±1.4%)drug release after 5 min yet only increased to a maximum of 43.7%(±2.3%)release after 30 min.For many drugs,this low percentage of drug release after this time would be deemed unsuitably low and may limit bioavailability.Through undertaking the formulation process with the drug alone,i.e.hydration,heating, filtering,drying then sieving,the percentage of drug release,or more accurately in this case,dissolution after 30 min was 43.8%(±7.9%).Therefore,it has been confirmed that exposure of the drug to the formulation process did not affect the pro file observed,i.e.hydrating through sieving did not enhance the effects observed for phenylbutazone.AllthreeSyloid®silicabasedformulations exhibited a dramatic enhancement in percentage dissolution,confirming that the presence of Syloid®silica contributed to the increase.Firstly,XDP 3150 achieved a percentage release of 42.4%(±1.9%)after only 5 min,i.e.almost equal to that observed for drug alone after twice as long.After a period of 30 min,this value had increased to 78.3%(±2.2%),far higher than that seen for drug alone or drug that had undergone the formulation process.Secondly,Syloid®silicas AL-1 FP did not show such a promising percentage release after 5 min(30.2%(±1.2%))compared with XDP 3150,yet after a total of 30 min had exceeded the former Syloid®silica to reach a maximum percentage release of 86.0%(±4.2%).Finally,XDP 3050 was found to be the most successful Syloid®silica for enhancing percentage release with an impressive 49.4%(±0.8%)released after 5 min,i.e.greater than the total seen for pure drug after 30 min,increasing to a maximum of 99.6%(±3.0%)release after 30 min.

When determining why all three Syloid®silicas enhanced the percentageofdissolution following astandard formulation method,it would appear that the transformation from the crystalline to amorphous form(as evidenced by XRD and dissolution pro files of processed samples)plays a key role.This has been the conclusion of other researchers,when investigating alternative mesoporous materials[31],and fits well with the results from this work.However,when considering why the three Syloid®silicas did not facilitate the same increase in percentage release,it is more appropriate to consider their relative physicochemical properties,specifically those identified in Table 1.For example,AL-1 FP and XDP 3050 pore sizes are very different,in that AL-1 FP has small mesopores,i.e.a smaller pore volume and diameter compared with XDP 3050.Based on the pattern of increasing percentage release,i.e.from XDP 3150 to AL-1 FP to XDP 3050,it would appear that two properties of the Syloid®silicas may play a key role in controlling the process,namely,surface area and/or pore diameter.Interestingly,pore volume does not appear to be an influential factor for the rate and extent of dissolution,yet pore diameter is.In this work it appears that a large pore diameter,with a small surface area,maximises the extent of dissolution,which again, fits well with the findings of other studies with mesoporous microspheres[30,33].As a consequence of this,it is not only possible to dramatically enhance the rate and extent of dissolution,but also to vary the percentage depending upon the type of Syloid®silica used.Another potentially influential factor is the formation of aggregates which may affect the drug release pro file through the creation of particle aggregation.If this is the case,then it can be proposed that there are two unique structures within the formulation∶drug within pores and aggregates between particles which can both contribute to drug release.

Fig.5.Scanning electron microscope images(SEM)of(A)processed phenylbutazone,(B)phenylbutazone and AL-1 FP,(C)phenylbutazone and XDP 3150,and(D)phenylbutazone and XDP 3050 after formulations.

Fig.6.Phenylbutazone release pro files for phenylbutazone(PhB),processed phenylbutazone,and Syloid® silica based formulations–XDP 3050,XDP 3150 and AL-1 FP.Each data point represents the mean of triplicate results(±SD).

4.Conclusions

In summary,it has been confirmed that it is possible to formulate Syloid®silica based formulations to enhance the dissolution of a poorly soluble drug,in this case,phenylbutazone.Characterisation data implies that this enhancement is a result of a change in crystallinity and the ability of the drug to enter pores within the Syloid®silica structure.All three Syloid®silicas analysed demonstrated a dramatic increase in percentage release with their final percentage values linked to the Syloid®silica pore diameter and/or surface area.This finding can be of benefit for not only phenylbutazone-based equine formulations but potentially a far wider range of compounds that exhibit poor aqueous solubility,which will help alleviate bioavailability issues.To ensure that longterm stability is not a limiting factor for formulation possibilities,it is the intended subject of future sample analysis,using techniques such as XRD and SEM.

Conflicts of interest

The authors declare that there are no conflicts of interest.

[1]Y.Choudhari,H.Hoefer,C.Libanati,et al.,Mesoporous silica drug delivery systems.N.Shah,H.Sandhu,D.S.Choi,et al.(Eds.),Amorphous Solid Dispersions∶Theory and Practice,Springer,New York,2014∶665–693.

[2]W.Xu,J.Riikonen,V.P.Lehto,Mesoporous systems for poorly soluble drugs,Int.J.Pharm.453(2012)181–197.

[3]S.C.Shen,W.K.Ng,L.S.O.Chia,et al.,Applications of mesoporous materials as excipients for innovative drug delivery and formulation,Curr.Pharm.Des.19(2013)6270–6289.

[4]C.A.McCarthy,R.J.Ahern,R.Dontireddy,et al.,Mesoporous silica formulation strategies for drug dissolution enhancement∶a review,Expert Opin.Drug Deliv.13(2016)93–108.

[5]S.H.Cheng,W.N.Liao,L.M.Chen,et al.,pH-controllable release using functionalized mesoporous silica nanoparticles as an oral drug delivery system,J.Mater.Chem.21(2011)7130–7137.

[6]Y.Hu,X.Dong,L.Ke,et al.,Polysaccharides/mesoporous silica nanoparticles hybrid composite hydrogel beads for sustained drug delivery,J.Mater.Sci.52(2017)3095–3109.

[7]H.Wen,Y.Qiu,Adsorption of small drug particles at the surface of large excipients,Pharm.Technol.Eur.18(2006)39–44.

[8]A.Kiwilsza,A.Pajzderska,J.Mielcarek,et al.,Dynamical properties of nimodipine molecules con fined in SBA-15 matrix,Chem.Phys.475(2016)126–130.

[9]N.Murillo-Cremaes,A.M.López-Periago,J.Saurina,et al.,Nanostructured silica-based drug delivery vehicles for hydrophobic and moisture sensitive drugs,J.Supercrit.Fluids 73(2013)34–42.

[10]A.Patil,U.N.Chirmade,V.Trivedi,et al.,Encapsulation of water insoluble drugs in mesoporous silica nanoparticles using supercritical carbon dioxide,J.Nanomed.Nanotechnol.2(2011)111–119.

[11]R.J.Ahern,A.M.Crean,K.B.Ryan,The influence of supercritical carbon dioxide(SC-CO2)processing conditions on drug loading and physicochemical properties,Int.J.Pharm.439(2012)92–99.

[12]R.K.Kankala,Y.S.Zhang,S.B.Wang,et al.,Supercritical fluid technology∶an emphasis on drug delivery and related biomedical applications advanced healthcare,Materials 6(2017)1700433.

[13]Y.Choudhari,U.Reddy,F.Monsuur,et al.,Comparative evaluation of porous silica based carriers for lipids and liquid drug formulations,Mesoporous Biomater.1(2014)61–74.

[14]F.Monsuur,Y.Choudhari,U.Reddy,et al.,Solvent free amorphisation for pediatric formulations(minitablets)using mesoporous silica,Int.J.Pharm.511(2016)1135–1136.

[15]L.J.Waters,T.Hussain,G.Parkes,et al.,Inclusion of fenofibrate in a series of mesoporous silicas using microwave irradiation,Eur.J.Pharm.Biopharm.85(2013)936–941.

[16]C.A.McCarthy,W.Faisal,J.P.O'Shea,et al.,In vitro dissolution models for the prediction of in vivo performance of an oral mesoporous silica formulation,J.Control.Release 250(2017)86–95.

[17]R.K.Kankala,C.G.Liu,A.Z.Chen,et al.,Overcoming multidrug resistance through the synergistic effects of hierarchical pH-sensitive,ROS-generating nanoreactors,ACS Biomater.Sci.Eng.3(2017)2431–2442.

[18]R.K.Kankala,Y.Kuthati,C.L.Liu,et al.,Killing cancer cells by delivering a nanoreactor for inhibition of catalase and catalytically enhancing intracellular levels of ROS,RSC Adv.5(2015)86072–86081.

[19]T.Limnell,H.A.Santos,E.Mäkilä,et al.,Drug delivery formulations of ordered and nonordered mesoporous silica∶comparison of three drug loading methods,J.Pharm.Sci.100(2011)3294–3306.

[20]P.Kinnari,E.Mäkilä,T.Heikkilä,et al.,Comparison of mesoporous silicon and non-ordered mesoporous silica materials as drug carriers for itraconazole,Int.J.Pharm.414(2011)148–156.

[21]W.Xu,J.Riikonen,V.P.Lehto,Mesoporous systems for poorly soluble drugs,Int.J.Pharm.453(2013)181–197.

[22]L.R.Soma,C.E.Uboh,G.M.Maylin,The use of phenylbutazone in the horse,J.Vet.Pharmacol.Ther.35(2012)1–12.

[23]C.Castagnetti,J.Mariella,Anti-inflammatory drugs in equine neonatal medicine.Part I∶nonsteroidal anti-inflammatory drugs,J.Equine Vet.Sci.35(2015)475–480.

[24]L.C.Sanchez,S.A.Robertson,Pain control in horses∶what do we really know?Equine Vet.J.46(2014)517–523.

[25]J.C.De Grauw,J.P.A.M.van Loon,C.H.A.van de Lest,et al.,In vivo effects of phenylbutazone on inflammation and cartilage-derived biomarkers in equine joints with acute synovitis,Vet.J.201(2014)51–56.

[26]U.Domańska,A.Pobudkowska,A.Pelczarska,et al.,Modelling,solubility and pKa of five sparingly soluble drugs,Int.J.Pharm.403(2011)115–122.

[27]S.Khan,H.Batchelor,P.Hanson,et al.,Physicochemical characterisation,drug polymer dissolution and in vitro evaluation of phenacetin and phenylbutazone solid dispersions with polyethylene glycol 8000,J.Pharm.Sci.100(2011)4281–4294.

[28]M.Van Speybroeck,V.Barillaro,T.D.Thi,et al.,Ordered mesoporous silica material SBA-15∶a broad-spectrum formulation platform for poorly soluble drugs,J.Pharm.Sci.98(2009)2648–2658.

[29]S.L.Longhofer,C.R.Reinemeyer,S.V.Radecki,Evaluation of the palatability of three nonsteroidal anti inflammatory top-dress formulations in horses,Vet.Ther.9(2008)122–127.

[30]Y.Hu,J.Wang,Z.Zhi,et al.,Facile synthesis of 3D cubic mesoporous silica microspheres with a controllable pore size and their application for improved delivery of a water-insoluble drug,J.Colloid Interface Sci.363(2011)410–417.

[31]Z.Guo,X.M.Liu,L.Ma,et al.,Effects of particle morphology,pore size and surface coating of mesoporous silica on Naproxen dissolution rate enhancement,Colloids Surf.B∶Biointerfaces 101(2013)228–235.

[32]R.Al-Oweini,H.El-Rassy,Synthesis and characterization by FTIR spectroscopy of silica aerogels prepared using several Si(OR)4 and R′′Si(OR′)3 precursors,J.Mol.Struct.919(2009)140–145.

[33]A.Martín,R.A.García,D.S.Karaman,et al.,Polyethyleneimine-functionalized large pore ordered silica materials for poorly water-soluble drug delivery,J.Mater.Sci.49(2014)1437–1447.