PLoS OnePLoS ONEplosplosonePLoS ONE1932-6203Public Library of ScienceSan Francisco, CA USA26161510449889910.1371/journal.pone.0130253PONE-D-15-16158Research ArticleSynthesis Mechanism and Thermal Optimization of an Economical Mesoporous Material Using Silica: Implications for the Effective Removal or Delivery of IbuprofenLow-Cost Mesoporous Material for the Removal or Delivery of IbuprofenKittappaShanmuga1CuiMingcan2RamalingamMalarvili3IbrahimShaliza1KhimJeehyeong2YoonYeomin4SnyderShane A.5JangMin16*Department of Civil Engineering, Faculty of Engineering, University of Malaya, Kuala Lumpur, MalaysiaSchool of Civil, Environmental, and Architectural Engineering, Korea University, Seoul, Republic of KoreaDepartment of Chemistry, Jalan Sultan Petaling Jaya, Selangor, MalaysiaDepartment of Civil and Environmental Engineering, University of South Carolina, Columbia, South Carolina, United States of AmericaDepartment of Chemical and Environmental Engineering, University of Arizona, Tucson, Arizona, United States of AmericaNanotechnology and Catalysis Research Centre (NANOCAT), University of Malaya, Kuala Lumpur, MalaysiaTajmir-RiahiHeidar-AliEditorUniversity of Quebec at Trois-Rivieres, CANADA

Competing Interests: The authors have declared that no competing interests exist.

Conceived and designed the experiments: MJ. Performed the experiments: SK. Analyzed the data: MC MR SI JK YY SAS MJ. Contributed reagents/materials/analysis tools: MJ. Wrote the paper: SK YY MJ.

* E-mail: minjang@um.edu.my
10720152015107e0130253144201519520152015Kittappa et alThis is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Mesoporous silica materials (MSMs) were synthesized economically using silica (SiO2) as a precursor via a modified alkaline fusion method. The MSM prepared at 500°C (MSM–500) had the highest surface area, pore size, and volume, and the results of isotherms and the kinetics of ibuprofen (IBP) removal indicated that MSM–500 had the highest sorption capacity and fastest removal speed vs. SBA–15 and zeolite. Compared with commercial granular activated carbon (GAC), MSM–500 had a ~100 times higher sorption rate at neutral pH. IBP uptake by MSM–500 was thermodynamically favorable at room temperature, which was interpreted as indicating relatively weak bonding because the entropy (∆adsS, –0.07 J mol–1 K–1) was much smaller. Five times recycling tests revealed that MSM–500 had 83–87% recovery efficiencies and slower uptake speeds due to slight deformation of the outer pore structure. In the IBP delivery test, MSM–500 drug loading was 41%, higher than the reported value of SBA–15 (31%). The in vitro release of IBP was faster, almost 100%, reaching equilibrium within a few hours, indicating its effective loading and unloading characteristics. A cost analysis study revealed that the MSM was ~10–70 times cheaper than any other mesoporous silica material for the removal or delivery of IBP.

This research was funded by the Malaysian Government Ministry of Higher Education through the High Impact Research Grant (D000062-16001), and partly supported by the Basic Science Research Program through a National Research Foundation of Korea (NRF) grant funded by the Korea government (MEST, No. KRF-2009-0092799). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Data AvailabilityAll relevant data are within the paper.
Data Availability

All relevant data are within the paper.

Introduction

Mesoporous silica materials (MSMs) have large, ordered pores, ranging from 2 to 50 nm, as classified by the International Union of Pure and Applied Chemistry (IUPAC) [1, 2]. MSMs have thick framework walls with interconnected channel structures and a high surface area that can provide superior properties for adsorption, catalysis, and sensing [1, 3, 4]. Because of evenly distributed pores (6–10 nm in size), MSMs readily allow organic molecules, including pharmaceuticals, to penetrate through their pores, resulting in less pore diffusion [1, 5].

In addition to environmental applications (e.g., water treatment), mesoporous materials have emerged as suitable candidates for drug delivery systems and other biomedical applications [6, 7]. Reference materials such as SBA–15, MCM–41, and MCM–48 have been studied extensively for the loading and unloading of pharmaceuticals [3, 4, 79]. However, these media face a major hurdle in actual application due to the use of expensive organometallic chemicals, such as tetraethyl orthosilicate (TEOS, Si(OC2H5)4) [3, 10], which is used as a cross-linking agent in the synthesis process [2, 11]. Thus, a more economical approach is essential for synthesizing MSMs if they are to be used in pharmaceutical removal or delivery processes.

Ibuprofen (IBP) was selected as a model drug for this study because it is one of the most widely consumed drugs worldwide; it is classified as a nonsteroidal anti-inflammatory drug (NSAID) [12, 13]. Indeed, due to the high concentration of IBP in water sources, it has been classified as a pharmaceutical pollutant by the World Health Organization (WHO). Moreover, IBP has been studied widely for loading and unloading purposes as a model drug [14].

To synthesize MSMs for IBP removal or delivery economically, we developed a modified alkaline fusion method, which broadens the choice of framework precursors. For the first time, we synthesized MSMs using SiO2 as a framework precursor, replacing TEOS. The main objectives of this investigation were to (i) assess the effects of thermal treatment (calcination) through characterization using X-ray diffraction (XRD), transmission electron microscopy (TEM), Fourier transform infrared spectroscopy (FTIR), and N2 gas isotherm; (ii) assess the sorption capacities and rates of IBP in multi-recycle runs; (iii) study IBP loading and unloading by MSMs; and finally (iii) determine the adsorption mechanism.

Materials and Methods

Pluronic P123 (EO20PO70EO20) and SiO2 were obtained from Sigma-Aldrich (St. Louis, MO, USA) and R&M Chemicals (Edmonton, AB, Canada), respectively. Sodium hydroxide (NaOH) and hydrochloric acid (HCl) (37%) were purchased from Merck (Darmstadt, Germany). Analytical grade IBP was supplied by Alfa Aesar (Royston, Hertfordshire, UK).

Synthesis method of the MSMs

MSMs were synthesized by reacting Pluronic P123 with SiO2 in the presence of NaOH and HCl. SiO2 (1 M) was dissolved in a 1 M NaOH solution and stirred using a magnetic stirrer at 45°C for 20 h. The pore-templating agent was prepared separately using 4 g of Pluronic P123, dissolved in 120 mL of a 2 M HCl solution with continuous stirring at 45°C. The two solutions were then mixed and stirred for 3 h at 45°C, and then for an additional 12 h at room temperature. The molar ratio of chemicals used in the synthesis of the MSM was 1 SiO2:1 NaOH: 5.28 HCl: 0.015 Pluronic: 200 H2O. The solution was then aged in a Teflon bottle at 90°C for 20 h. The precipitated solid product was recovered by filtering it using a 0.45 μm pore size cellulose acetate membrane filter. Finally, the material was washed with deionized (DI) water and ethanol (50%) and dried at 60°C for 24 h. A Western furnace was used to calcine the dried samples at temperatures ranging from 500 to 900°C for 4–6 h.

Characterization of the MSMs

XRD analysis was performed using an X-ray diffractometer (Empyrean; PANalytical, Almelo, The Netherlands). The samples were scanned at 2θ from 0.5 to 2.5° using a step size of 0.0070° and a scanning time of 19.9260 s. Nitrogen adsorption and desorption isotherms were measured using a TriStar II 3020 system (Micromeritics, Norcross, GA, USA). The samples were analyzed at 77.35 K. The adsorption and desorption data were calculated using Brunauer–Emmett–Teller (BET) theory to determine the specific surface area of the samples, and the pore-size distributions and pore volumes were determined using the Barrett–Joyner–Halenda (BJH) theory. Infrared (IR) spectra were obtained using a NICOLET IS 10 spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). Microscopic images of nanoscale pore structures were taken using a transmission electron microscope (HT 7700 TEM; Hitachi, Tokyo, Japan) at 120 kV.

Preparation of IBP stock solution

The required amount of IBP powder (200 mg) was dissolved in 10% methanol (CH3OH; Thermo Fisher Scientific) and dissolved in DI water in a 250 mL volumetric flask. It was then stirred for 12 h, sonicated for 2 h, and filtered using a 0.45 μm pore-size cellulose acetate membrane filter. A sodium chloride (NaCl) solution was added to the filtrate to create ionic strength (0.01 M) in the solution. The pH was adjusted to 7 using a sodium phosphate (Na3PO4) solution.

Isotherms and kinetics of IBP removal by the MSMs

Adsorption isotherms and kinetics were determined for four synthesized MSMs calcined at 500, 600, 700, and 900°C, which were denoted as MSM–500,–600,–700 and–900, respectively.

Isotherms

In 20 mL vials, 10–500 mg of MSM and 10 mL of IBP solution (150 mg L–1) were added together. The vials were placed in an electric shaker and agitated at 250 rpm and 25°C. The suspension was filtered using a 0.45 μm pore-size cellulose acetate membrane filter, and the filtrate was analyzed for IBP.

The data were fitted with the Langmuir and Freundlich isotherms. When the adsorption is retained on a uniform monolayer surface, the Langmuir model fits the data of the isotherm. The maximum adsorption capacity is achieved when all sorption sites are saturated. The linear form of the Langmuir model can be written as qeq=QmaxKLCeq1+KLCeq where qeq is the amount of solute adsorbed per unit weight of adsorbent (mg g–1) at equilibrium, Ceq is the equilibrium concentration of the solute in the bulk solution (mg L–1), Qmax is the maximum adsorption capacity (mg g–1), and KL is the Langmuir constant related to the energy of adsorption.

Assuming that the adsorption is retained on the heterogeneous surface of the adsorbent, the Freundlich model fits the isotherm data better. In the Freundlich model, chemisorption and physisorption are pertinent to monolayer and multilayer adsorption, respectively. The linear form of the Freundlich equation can be expressed as logqeq=logKF+1nlogCeq where KF and n are Freundlich isotherm constants related to the adsorption capacity and adsorption intensity, respectively.

Kinetics

IBP solutions (200 mL) with concentrations of 100 mg L–1 were added to four separate beakers in which 1,000 mg of MSM had been pre-added. These four beakers were stirred constantly at 200 rpm for 200 min at room temperature. At timed intervals, samples from each of the four beakers were collected, filtered, and analyzed for IBP. After obtaining the kinetic data, a pseudo-second-order kinetic model equation was used to obtain the kinetic constants of the reaction, tqt=1K2qeq2+tqeq where K2 is the rate constant of pseudo-second-order adsorption (g mg–1 min–1), K2qeq2 or v0 (mg g–1 min–1) is the initial adsorption rate, and qt is the amount of adsorbate adsorbed at time t (min).

Thermodynamics

The kinetics of IBP uptake by MSM-500 were examined at 299, 309, and 319 K. Thermodynamic values, such as the Gibbs free energy (ΔadsG0), enthalpy (ΔadsH0), and entropy (ΔadsS0), were determined using the kinetic data. This adsorption process is represented by a reversible heterogeneous equilibrium: IBP in solution ↔ IBP adsorbed.

The thermodynamic calculation was adapted from a mechanism and calculation reported previously [15]. The equilibrium constant is defined using the following equation [16, 17] K0=qeCs where qe and Cs are concentrations of adsorbed IBP and in solution at equilibrium, respectively. The K value was substituted to acquire ΔadsG0, ΔadsH0, and ΔadsS0 using the following equations, ΔadsGo=RTlnK0 lnK0=ΔadsS0RΔadsH0RT where R and T are the gas constant and absolute temperature, respectively. adsH0 and adsS0 were obtained from the slope and intercept of the line of ln K vs. 1/T.

Regeneration

MSM–500 was used to test the regeneration of MSM for IBP adsorption. After being used for the first adsorption test, MSM–500 was filtered, washed in methanol for 4–6 h, and dried overnight in an oven. The same medium was reused for four subsequent adsorption tests.

Drug loading and in vitro drug release

To load IBP, 0.2 g of MSM–500 was added to a hexane solution containing 30 mg mL–1 IBP. The powder-containing solution was soaked for 2 days and stirred at 150 rpm at room temperature until the concentration in the solution remained constant. The loaded amount of IBP was determined using the difference in the initial and final IBP concentration. The powder was then washed with the hexane solution and dried under vacuum.

The in vitro drug release test was conducted by soaking 0.2 g of IBP containing MSM–500 in 100 mL of phosphate-buffered saline (PBS) at pH 7.4. The solution was stirred at 150 rpm and kept at 36.9°C, similar to body temperature. At each time point, 3 mL of solution was extracted and analyzed, and an equal volume of fresh solution was replaced. The IBP was analyzed to assess the release of IBP. The test was also repeated with the PBS solution at a lower pH, i.e., 4–6.

Analysis of IBP

IBP was analyzed using a UV 2600 spectrophotometer (Shimadzu, Otsu, Japan). A standard IBP solution was subjected to a full-spectrum scan, ranging from 900 to 0 nm wavelengths, and maximum absorbance (λmax) was detected at 264 nm. A calibration curve, which was obtained by analyzing various concentrations of the prepared IBP standard solution, showed good linearity, with a determination coefficient (R2) of 0.9997.

Results and DiscussionThermal effect analysis and characterization of the MSMs

Fig 1(A) shows a low-angle XRD diffraction pattern of selected MSMs. Among the synthesized media, as-synthesized MSM, MSM–500, and MSM–600 showed well-resolved peaks indexed at the d(100) reflection [1, 7, 18]. Additionally, MSM–500 had two more weak peaks at 2θ angles of 1.08 and 1.58°, indexed at d(110) and d(200), respectively. These types of XRD reflection patterns can be assigned to a well-ordered mesoporous structure with P6mm hexagonal symmetry [1, 2, 19].

10.1371/journal.pone.0130253.g001

(A) XRD patterns of as-synthesized mesoporous material and MSMs calcined at different temperatures, (B) the N2 sorption isotherm, and (C) pore-size distribution of the MSMs (A: MSM-as synthesized, B: MSM-500, C: MSM-600, D: MSM-700, E: MSM-800, and F: MSM-900).

The XRD results indicated that MSM–500 had the highest peak intensity compared with the other MSMs, reflecting the formation of more uniform mesopores [4, 19]. The shifting of peaks to lower 2θ values with higher calcination temperature shows an increase in the pore sizes of the MSMs, indicating the possibility that MSM–500 had larger pores than MSM–600 and –700. However, MSMs calcined at >700°C showed no XRD diffraction peak, demonstrating the collapse of the pore structure [19]. This phenomenon was investigated further using N2 adsorption–desorption analysis.

Park et al. reported that synthesis of mesoporous materials using silica glass at pH >10 yielded a white gel that transformed into a mesoporous material due to a hydrothermal process [20]. However, in the current method, we synthesized the MSMs under acidic conditions (pH <1), whereby a clear white precipitate was formed immediately after the silicate (HSiO3) and acidified Pluronic P123 solutions were mixed together. The precipitate was treated using a hydrothermal process, in which it was heated and condensed to form the MSM. This method allows the excess hydrogen (H+) and chlorine (Cl) to coordinate electrostatically between the nonionic pore-templating agent (So) and the hydrogen-associated silicic acid (Io, Si(OH)4), which may finally lead to the formation of SoH+ClIo. Given these observations, the reaction mechanism of this synthesis can be represented using reactions 1–3.

First, the mixing of NaOH and SiO2 may lead to the formation of silicate ions (HSiO3) (R1).

SiO2+NaOHNa++HSiO3

Second, when both the silicate and acidic Pluronic solution are mixed together, silicic acid (Si(OH)4) is formed under acidic conditions (R2).

HSiO3+HCl+H2OSi(OH)4+Cl

Finally, through a hydrothermal process, two silicic acids are condensed and form a silicate dimer, which further polymerizes to form the framework of the MSM (R3).

Si(OH)4+Si(OH)4Δ90°CSi2O(OH)6+H2O

Fig 1(B) shows the N2 adsorption–desorption isotherms of the various MSMs. According to the IUPAC classification, the as-synthesized MSM, MSM–500, and MSM–600 had a type-IV sorption pattern [21, 22], whereas MSMs calcined at >700°C showed no sorption profile. At a relative pressure of P/Po < 0.2, MSM–500 had more primary micropores (<2 nm) than MSM–600. Reversible sorption of N2 gas for MSM–500 and–600 occurred up to a P/Po of 0.6, whereas that for as-synthesized MSM happened much earlier, at a P/Po of 0.2. These results indicate that MSM–500 and–600 had more uniform pores with a hysteresis loop pattern type H1, similar to SBA–15 [23, 24]. However, this MSM behavior differed from that of MCM–41, in which the adsorption and desorption isotherm process were fully reversible [25].

Using the BJH method, the pore-size distribution was obtained for as-synthesized MSM, MSM–500, and–600 (Fig 1(C)). These results showed that MSM–500 had the narrowest pore-size distribution, with a primary peak at 6.2 nm, while MSM–600 displayed a broader distribution and a primary peak at 5.6 nm. This indicates that MSM–500 had a larger pore structure and pore volume than MSM––600, consistent with the XRD results in Fig 1(A).

The pore structure characteristics of MSMs and other referenced media are summarized in Table 1. Among the calcined MSMs, MSM–500 had the highest BET specific surface area (310 m2 g–1) and pore volume (0.541 cm3 g–1). The BET specific surface area and pore volume decreased substantially when the calcination temperature increased from 500 to 900°C. A significant decrease was noted at 700°C, with the BET specific surface area reduced to <10 m2 g–1. Generally, the MSMs synthesized had smaller BET specific surface areas and pore volumes than reported for SBA–15 [9, 18].

10.1371/journal.pone.0130253.t001Structural parameters of as-synthesized and calcined MSMs.
d-MeanWallSurfaceVteVm
spacingaporethicknesscareadPore volumeMicropore
(nm)diameterb (nm)(nm)(m2 g–1)(cm3 g–1)(cm3 g–1)
MSM as-synthesized11.95.787.93260.5410.025
MSM-5009.686.264.93100.5650.011
MSM-6009.336.933.81600.3390.011
MSM-700n/a4.76n/a6.30.00890.0069
MSM-800n/a1.98n/a1.70.00390.0029
MSM-900n/an/an/a0.70.00060.0005

a (100) interplanar spacing.

bCalculated from desorption of the N2 isotherm.

cDetermined from the difference between the unit cell parameter (ao = 2d100/√3) and the frame work pore size.

dBET specific surface area calculated from adsorption of the N2 isotherm.

eVt = total pore volume calculated using BJH.

The corresponding unit cell parameter (ao) and mean BJH desorption pore diameter were used to determine the wall thickness of the MSMs. The wall thickness was calculated as the difference between ao, where ao = 2 × d(100)/3, and the mean BJH desorption pore diameter [2, 26]. The mean BJH desorption pore diameter increased during the calcination process at 500°C and 600°C (6–7 nm) for MSMs and was smaller than that of SBA–15 (~10 nm). At the same temperature, the MSMs d-spacing and wall thickness decreased from 11.9 to 9.3 nm and 7.6 to 3.3 nm, respectively. This reduction in wall thickness is a result of silicate condensation during calcination [27]. Calcination at temperatures >700°C caused a significant decrease in the wall thickness and the framework to collapse [2, 28].

In this study, MSMs were found to have lower thermal stability (up to 600°C), which could be due to the formation of weaker bonds between the condensed silica walls. In contrast, Tung et al. reported that SBA––15 and a few other mesoporous silica-based materials had high thermal stabilities of up to 850°C [1].

Adsorption isotherms of IBP by the MSMs

The IBP adsorption isotherms were determined using MSMs calcined at different temperatures (Fig 2(A)). The Qeq of IBP by MSM–500 and–600 increased with increasing Ceq, although differences were observed in the adsorption capacity. However, MSM–700 and–900 showed no significant adsorption, consistent with the physical property analyses using XRD and N2 adsorption–desorption analysis. To better understand the adsorption phenomenon of these MSMs, Langmuir and Freundlich isotherms were fitted with the experimental data of MSM–500 and–600.

10.1371/journal.pone.0130253.g002

(A) Adsorption isotherms of the MSMs calcined at different temperatures. (B) Kinetics of IBP uptake by MSMs calcined at different temperatures.

The Langmuir and Freundlich parameters, along with determination coefficients (R2) of the linear plots, were calculated. The results are presented in Table 2. The adsorption capacities (Qmax) for MSM–500 and–600 varied considerably: 64.8 mg g–1 for MSM–500 and 31.9 mg g–1 for MSM–600.

10.1371/journal.pone.0130253.t002Isotherm parameters obtained by fitting equilibrium data with the Freundlich and Langmuir isotherms for the adsorption of IBP on MSMs.
SampleLangmuir equationFreundlich EquationRef.
pHKL (L mg–1)qm (mg g–1)R2akF1/nR2
MSM-50070.01064.80.970.930.820.95This study
MSM-60070.00931.90.740.270.930.75This study
SBA-15513.3300.410.971.500.780.97[1]
MCM-410.328[7]
GACn/an/an/a0.0720.870.99[7]
CAC~40.262153.20.99939.10.380.890[29]
CPAC~40.112416.70.99856.90.480.919[29]
AC-olive waste4.10.5512.60.9716.580.340.956[30]
Commercial GAC0.611604[31]
Cork-Based AC40.356139.236.60.303[32]

a kF in mg1–1/n L1/n g–1, KL: Langmuir constant, kF and n: Freundlich constants, qmax: maximum amount of adsorbate.

As the value of KL increased, the adsorption affinity of IBP by the adsorbent increased. The results show that MSM–500 had a higher affinity than MSM–600, with a higher KL value (Table 2).

However, the adsorption tests of MSMs were conducted at neutral pH, so comparing the results with reported KL values measured under acidic conditions is difficult.

At higher pH, two important factors, ionization and surface charge of the adsorbent, can affect the adsorption process. The acid dissociation of IBP molecules occurs at a pH higher than the pKa value (4.9), which may increase IBP solubility. As a result, a weaker interaction between the MSM surface and IBP would occur, and the adsorption capacity would decrease significantly. However, the Qmax of MSM–500 was higher than the value reported for SBA–15.

As a function of adsorption strength, the constant values of 1/n obtained from the Freundlich model for IBP removal by MSM–500 and–600 were 0.82 and 0.93, respectively. These values were larger than the constant values of other media except commercial granular activated carbon (GAC; Table 2). If 1/n equals 1, then the partition between the two phases is considered independent of the concentration. A value of 1/n > 1 demonstrates “normal” adsorption, while a value of 1/n < 1 indicates cooperative adsorption. Thus, based on the evaluation of 1/n, the adsorptive removal of IBP by MSM–500 and–600 can be concluded to be favorable [33, 34].

The adsorption densities per unit surface area for MSM–500 and–600 were similar, 0.21 and 0.20 mg m–2, respectively. These values were comparable to the reported value (0.28–0.32 mg m–2) for SBA–15 [7, 8].

Adsorptive kinetics of IBP by the MSMs

Fig 2(B) shows the kinetics of IBP removal by the MSMs at pH 7. Rapid adsorption of IBP by MSM–500 and–600 was observed for the first 5 min, reaching a plateau after 15 min. After 60 min, the IBP adsorption increased again, suggesting the possibility of IBP molecules forming a dimer [7, 19, 35]. Thus, the contact time needed for MSM–500 and–600 was much shorter than that of activated carbon (AC), which requires ~2 h [29]. The faster adsorption of MSM–500 and–600 may relate to their pore structures. MSMs have mostly mesopores and only 3–5% micropores, while AC has predominantly micropores (as high as 70%), leading to slow diffusion. The pore size of MSMs (~7 nm) is uniform and large, and sufficient for the penetration of IBP molecules (~1 nm). The pore size of MSMs allows the IBP molecules to diffuse readily and adsorb onto the surface of MSMs within a short time, which indicates that MSM–500 and–600 may be advantageous in terms of designing a compact water treatment system, enhancing their economical use commercially. Furthermore, the homogenous pore structures may make them suitable candidates for biomedical applications as drug carrier media.

The adsorption kinetics were investigated using a pseudo-second-order kinetic model, in which t/qt vs. t was plotted to obtain the rate parameters. The results showed that the fit of the model to the data had a high determination coefficient (0.99).

Table 3 lists the important parameters obtained from the current kinetic study and other references.

10.1371/journal.pone.0130253.t003Kinetic parameters obtained after fitting to the pseudo-second-order model.
SampleC0 (mg dm–3)pHK2 (g mg–1 min–1)R2K2 qeq2 or v0 (mg g–1 min–1)qe,calc. (mg g–1)Ref.
MSM-50010070.100.99943.520.4This study
MSM-60010070.030.9999.5817.7This study
MSM-70010070.020.9650.925.6This study
MSM-90010070.010.6220.192.5This study
CAC90~40.070.999834112.4[36]
CPAC90~40.230.9992500106.4[36]
Zeolites MNCZ0.170.2140.9990.098[37]
Commercial GAC70.00110.95269.95[31]
Cork-based GAC9040.070.999834112.4[32]

The K2 values for MSM–500 and–600 were 0.10 and 0.03 g mg–1 min–1, respectively. The higher K2 value for MSM–500 indicates that the medium was much faster in adsorbing IBP than the other MSMs. Although several types of AC have higher K2 values than MSMs (Table 3), most of the tests were conducted under acidic conditions (pH ~4). At pH 7, MSM–500 had a ~100 times higher sorption speed than commercial GAC, which had a K2 value of 0.0011 g mg–1 min–1. Zeolites had a higher K2 value than the MSMs, but the Qmax was only 0.0976 mg g–1.

Thermodynamics

Fig 3 shows kinetic data of IBP adsorption at 299, 309, and 319 K, and fit lines of the pseudo-second-order kinetic model. The qeq of IBP uptake by MSM–500 was 22.8 mg g–1 at 299 K, but decreased to 18.28 and 7.7 mg g–1 at 309 and 319 K, respectively. This phenomenon has been seen in several studies on removing various pharmaceuticals with silica-based media and multi-wall carbon nanotubes (MWCNTs) [15, 3840]. Based on the calculations of thermodynamic parameters, the adsG0 values were –10.4, 1.01, and 11.3 kJ mol–1 at 299, 309, and 319 K, respectively. Thus, the IBP adsorption is a spontaneous physical process at room temperature, but becomes non-spontaneous at higher temperatures. Moreover, adsG0, adsH0, and adsS0 were also calculated for comparison with other references (Table 4).

10.1371/journal.pone.0130253.g003Temperature effect of IBP adsorption onto MSM-500.10.1371/journal.pone.0130253.t004Thermodynamics parameters (<italic>∆</italic><sub><italic>ads</italic></sub><italic>H</italic><sup><italic>0</italic></sup> and <italic>∆</italic><sub><italic>ads</italic></sub><italic>S</italic><sup><italic>0</italic></sup>) of IBP uptake by MSM-500 and comparison with other references.
SampleTemp. (K)EDCΔadsH0 (kJ mol–1)ΔadsS0 (J mol–1 K–1)Reference
UMSa290–313Molsidomine–15–27[15]
MMSb290–313Molsidomine–44–127[15]
PhMSc290–313Molsidomine–65–186[15]
β-CDd298–323IBP–52.08–98.22[38]
MWCNTe288–318Bisphenol AF–17.16–53.7[39]
MWCNT278–333Estrone–14.16–560[40]
MWCNT278–33317β-estradiol–12.08–460[40]
MWCNT278–33317a-ethinylestradiol–10.11–420[40]
MSM500299–319IBP–23.0–0.07This study

UMSa: unmodified silica

MMSb: mercaptopropyl grafted silica material

PhMSc: phenol-modified silica

β-CDd: beta cyclodextrin

MWCNTe: multi-wall carbon nanotube.

The calculated adsH0 was –23.0 kJ mol–1, indicating that the IBP adsorption onto MSM–500 is an exothermic reaction, and adsS0 was found to be –0.07 J mol–1 K–1. Accordingly, as a consequence of IBP adsorption onto the pores of MSM–500, the randomness of IBP decreased because the degree of freedom or mobility of IBP is limited compared with IBP in solution. The adsS0 value was much smaller than other references, indicating that less energy is needed to restore the system back to its original state. The negative adsH0 and adsS0 suggest that the removal process of IBP can be considered an enthalpy-driven process.

Regeneration

Fig 4(A) presents kinetic data for five cycles of the adsorption of IBP by MSM–500 and fit lines of the pseudo-second-order kinetic model. Fig 4(B–D) show the qeq, K2, and v0, respectively, at each cycle. The first cycle had the highest qeq (22.8 mg g–1). This result is similar to the kinetic tests conducted previously, suggesting that the experiments were reproducible under the same conditions. All consecutive cycles showed a qeq of 19–19.8 mg g–1, slightly lower than the qeq for the first cycle. Thus, based on the first qeq, MSM–500 had recovery efficiencies of 83–87%.

10.1371/journal.pone.0130253.g004IBP uptake by MSM-500 at several recycles: (A) kinetics, (B) <italic>q</italic><sub><italic>eq</italic></sub>, (C) <italic>K</italic><sub><italic>2</italic></sub>, and (D) <italic>v</italic><sub><italic>0</italic></sub> according to the number of recycles.

In the first cycle, 90% of the qeq for IBP uptake was achieved at 5 min, but it slowed to ~30 min in the third and fourth cycles. In detail, the values of K2 and v0 were comparable in the first and second cycles, but dropped considerably at the third and decreased further at the fourth and fifth cycles. Accordingly, the reduction trends of IBP uptake capacity and speed were different, and these differences might have been affected by the deformation of the pore structure or the existence of some IBP molecules trapped in the pores after extraction with methanol. AC is composed mostly of micropores, and IBP molecules are adsorbed at the supermicropores, because its critical dimension is 0.72 nm [41]. Thus, extraction-based regeneration is not practical for spent AC due to the relatively high adsorption energy of AC. Generally, loaded AC is regenerated ex situ using heat or steam, which is a high energy-consuming process. Also, during this process, a large proportion of the AC can be lost [36]. In terms of regeneration, however, MSM–500 has advantages over AC due to simple and fast regeneration using methanol.

Fig 5 shows TEM images of MSM–500 before and after the fifth adsorption of IBP. As shown in Fig 5(A), MSM–500 had homogeneous pore and wall structures. The measured pore size was in the range of 5.44 ± 0.4 nm, ~11% smaller than the primary pore size (6.2 nm) as determined by the N2 gas isotherm. The TEM image of MSM–500 obtained after the fifth adsorption showed slight deformation at the inlet side of the pore structure, but no major defect. Thus, this result may be linked to the kinetic results for re-adsorption.

10.1371/journal.pone.0130253.g005TEM images of MSM-500 (A) before and (B) after the fifth adsorption of IBP.
Mechanism of adsorption of IBP by the MSMs

Fig 6 shows the FTIR spectra (4,000–500 cm–1) of all the MSMs prepared at different calcination temperatures, as well as those of IBP-retained MSMs. The as-synthesized MSM had obvious bands for OH (~3,400 cm–1), H2O (1,644 cm–1), Si–O–Si (symmetric stretching vibrations (vs.) at 1,057 cm–1 and ~800 cm–1), Si–OH (vs. at ~950 cm–1), and O–Si–O (deformation vibration (δ) at 435cm–1). Si–O–Si (lines a and c in Fig 6) and Si–OH (line b) of MSM-500 and -600 show FTIR spectra typical of MSMs [9]. The intensity of the Si–O–Si peak for the MSM calcined at >600°C decreased markedly, indicating dissociation of the Si–O–Si bond during calcination.

10.1371/journal.pone.0130253.g006FTIR spectra of MSM materials before (A) and after adsorption (B) of IBP.

After the adsorption of IBP, the Si–O–Si (lines a′ and c′ in Fig 6) and Si–OH (line b′) peaks were preserved in the MSMs, showing stability of the pore structure. A new band, centered at 1,630 cm–1 (box labeled as d′), can be attributed to the asymmetric vibration of the carboxylic functional group (COO). Two more sharp peaks, at 1,435 and 1,500 cm–1, were assigned to the C–H bonding of carbon atoms. The peak at 2,920 cm–1 (box labeled as e′) was assigned to a strong C–H bond, confirming the adsorption of IBP into the channels of MSM–500 [8, 9, 42]. In contrast, MSM–600 had a weaker adsorption peak, as shown in the d′ and e′ boxes, while the other MSMs showed no peaks in this range.

These FTIR data suggest that the reaction mechanism between IBP and MSM is a hydrophilic reaction [1]. At pH 7, the hydrogen of the carboxylic group of IBP is dissociated to form a COOfunctional group. The formation of a hydrogen bonding reaction between the COOof IBP and silanol groups on the pore surface is favorable because it requires a lower activation energy [1, 7], SiOH+OOCBCHSiOHOOCBCH where ≡Si–OH and−OOC–B···CH are silanol and IBP, respectively. Fig 7 shows the overall mechanisms of MSM synthesis and IBP adsorption.

10.1371/journal.pone.0130253.g007Mechanism of MSM synthesis and IBP adsorption.
IBP loading and release by the MSMs

MSM–500 showed a 41% IBP drug loading capacity, comparable to the value reported for SBA–15 (20–30%) [6]. Fig 8 shows the in vitro drug release of the MSMs: almost 100% released within a few hours. Compared with other mesostructured materials, such as SBA–15 (~10 h), the unloading time of MSM–500 was much shorter. MSM–500 has a larger pore size than SBA–15 and MCM–41, which increases the diffusion of IBP molecules through the pores and subsequently increases the loading and unloading rates. The release behavior was also affected by the pH of the solution. At a lower pH, the IBP release was faster and higher, which is advantageous because most tumor cells in the body are slightly acidic and the drug release can be more specific to such target areas. Moreover pure SiO2 was utilized as a precursor, instead of TEOS to synthesis MSMs. This is expected to produce less toxic materials than the conventional mesoporous materials. Nakashima et al. have reported the acute and subchronic inhalation toxicity of TEOS in the synthesis of SBA–15, MCM–41 and MCM–48 [43].

10.1371/journal.pone.0130253.g008<italic>In vitro</italic> release of IBP from MSM-500 at different pH.
IBP removal and drug loading cost analysis

In this study, the initial synthesis costs of materials were calculated to estimate the removal or delivery costs for IBP, for which the IBP concentration and volume were 10 mg L–1 and 1 m3, respectively. As shown in Table 5, the MSM had the lowest cost ($0.76) vs. other silica-based mesoporous media.

10.1371/journal.pone.0130253.t005Comparison of pore characteristics and cost analyses for IBP removal.
MaterialsSpecific surfacePore vol.Conc. IBPtQmaxpH1Mass2CostsReference
area(cm3 g–1)(mg L–1)(min)(mg g–1)(g)(USD)
(m2 g–1)
Granular SBA– 157670.860.1140015766751.2[44]
3MSN4508171.7420042098.371028.22[45]
MSM3100.5411505–1064.871540.76This study

1Mass of media required for treating 1 m3 of water containing IBP (10 mg L–1)

2Costs (USD) of media for treating 1 L of IBP containing water (10 mg L–1), which does not include electricity consumption, personnel, or regeneration.

3MSN450: mesoporous silica nanoparticles synthesized with 450 W of microwave power.

For example, the MSM was 67 times cheaper than granular SBA-15. The reduction in cost was mainly due to the replacement of TEOS with the much cheaper silicon source, SiO2. The allocated price of TEOS is ~96% of the total costs for synthesis and TEOS is ~200 times more expensive than SiO2. Furthermore, the reusability of MSM will further reduce the cost significantly.

Conclusions

Adsorption isotherms and kinetics revealed that the MSMs synthesized showed relatively higher adsorption capacities and faster adsorption rates than other comparable media. The adsorption of IBP by MSM–500 was thermodynamically favorable at room temperature, but it involved relatively weak bonding because the calculated entropy was much smaller than in other references. In consecutive sorption studies, MSM–500 showed 83–87% recovery efficiencies, although it had a slower pattern of uptake. Based on the FTIR results, the reaction mechanism was identified as a hydrophilic interaction between the COOof IBP and Si–OH of the pore surface. MSM–500 was also found to be suitable for IBP drug loading and unloading purposes, with a 41% loading capacity and almost 100% unloading efficiency within a few hours. The MSM was synthesized successfully using SiO2, so that its production costs were much lower than that of SBA-15. Thus, the new synthesis route developed in this study could have high potential for cost-effective mass production, and the MSM could be a suitable medium for industrial applications (for water treatment and drug delivery), as shown in the cost analysis study [6, 46, 47]. Furthermore, the use of inert precursor SiO2 is expected to produce less toxic MSMs, which are more suitable for water treatment and biomedical application.

This project was funded by the Malaysian Government Ministry of Higher Education through the High Impact Research Grant (D000062-16001), and partly supported by the Basic Science Research Program through a National Research Foundation of Korea (NRF) grant funded by the Korea government (MEST, No. KRF-2009-0092799).

ReferencesBuiTX, ChoiH. Adsorptive removal of selected pharmaceuticals by mesoporous silica SBA-15. J Hazard Mater. 2009;168(2):6028.19327889ZhaoD, FengJ, HuoQ, MeloshN, FredricksonGH, ChmelkaBF, et al Triblock copolymer syntheses of mesoporous silica with periodic 50 to 300 angstrom pores. science. 1998;279(5350):54852. 9438845KatiyarA, YadavS, SmirniotisPG, PintoNG. Synthesis of ordered large pore SBA-15 spherical particles for adsorption of biomolecules. J ChromatogrA. 2006;1122(1):1320.ZhuY-F, ShiJ-I, LiY-S, ChenH-R, ShenW-H, DongX-P. Storage and release of ibuprofen drug molecules in hollow mesoporous silica spheres with modified pore surface. Microporous Mesoporous Mater. 2005;85(1):7581.Vallet-RegiM, RamilaA, Del RealR, Pérez-ParienteJ. A new property of MCM-41: drug delivery system. Chem Mater. 2001;13(2):30811.SongS-W, HidajatK, KawiS. Functionalized SBA-15 materials as carriers for controlled drug delivery: Influence of surface properties on matrix−drug interactions. Langmuir. 2005;21:956875. 16207037AnderssonJ, RosenholmJ, ArevaS, LindénM. Influences of material characteristics on ibuprofen drug loading and release profiles from ordered micro-and mesoporous silica matrices. Chem Mater. 2004;16(21):41607.Vallet-RegiM, BalasF, ArcosD. Mesoporous materials for drug delivery. Angew Chem Int Ed. 2007;46(40):754858.YangP, HuangS, KongD, LinJ, FuH. Luminescence functionalization of SBA-15 by YVO4: Eu3+ as a novel drug delivery system. Inorg Chem. 2007;46(8):320311. 17371013SonSJ, BaiX, NanA, GhandehariH, LeeSB. Template synthesis of multifunctional nanotubes for controlled release. J Control Release. 2006;114(2):14352. 16870299MargoleseD, MeleroJ, ChristiansenS, ChmelkaB, StuckyG. Direct syntheses of ordered SBA-15 mesoporous silica containing sulfonic acid groups. Chem Mater. 2000;12(8):244859.SadeckaJ, ČakrtM, HercegováA, PolonskýJ, SkačániI. Determination of ibuprofen and naproxen in tablets. J Pharmaceut Biomed. 2001;25(5):88191.KwonM, KimS, YoonY, JungY, HwangT-M, LeeJ, et al Comparative evaluation of ibuprofen removal by UV/H2O2 and UV/S2O82− processes for wastewater treatment. Chem Eng J. 2015;269:37990.WangX, LiuP, TianY. Ordered mesoporous carbons for ibuprofen drug loading and release behavior. Microporous Mesoporous Mater. 2011;142(1):33440.AlyoshinaNA, ParfenyukEV. Functionalized mesoporous silica materials for molsidomine adsorption: Thermodynamic study. J Solid State Chem. 2013;205:2116.LiuY, LiuY-J. Biosorption isotherms, kinetics and thermodynamics. Sep Purif Technol. 2008;61:22942.WalasSM. Phase equilibria in chemical engineering: Butterworth Boston; 1985.DoadrioA, SousaE, DoadrioJ, Pérez ParienteJ, Izquierdo-BarbaI, Vallet-RegıM. Mesoporous SBA-15 HPLC evaluation for controlled gentamicin drug delivery. J Control Release. 2004;97(1):12532. 15147810WangX, LiuP, TianY. Ordered mesoporous carbons for ibuprofen drug loading and release behavior. Microporous Mesoporous Mater. 2011;142(1):33440.ParkJ, HanY, KimH. Formation of Mesoporous Materials from Silica Dissolved in Various NaOH Concentrations: Effect of pH and Ionic Strength. J Nano Mat. 2012;2012:10 doi: 10.1155/2012/528174KresgeC, LeonowiczM, RothW, VartuliJ, BeckJ. Ordered mesoporous molecular sieves synthesized by a liquid-crystal template mechanism. Nature. 1992;359(6397):7102.KimS-S, PinnavaiaTJ. A low cost route to hexagonal mesostructured carbon molecular sieves. Chem Commun. 2001;(23):24189.ThommesM. Physical adsorption characterization of nanoporous materials. Chem Ing Tech. 2010;82(7):105973.JaroniecM. Characterization of Nanoporous Materials Access in Nanoporous Materials: Springer; 2002 p. 25572.KrukM, JaroniecM, SayariA. Application of large pore MCM-41 molecular sieves to improve pore size analysis using nitrogen adsorption measurements. Langmuir. 1997;13(23):626773.LimMH, SteinA. Comparative studies of grafting and direct syntheses of inorganic-organic hybrid mesoporous materials. Chem Mater. 1999;11(11):328595.ZhaoD, HuoQ, FengJ, ChmelkaBF, StuckyGD. Nonionic triblock and star diblock copolymer and oligomeric surfactant syntheses of highly ordered, hydrothermally stable, mesoporous silica structures. J Am Chem Soc. 1998;120(24):602436.ZhaoX, LuG, WhittakerA, MillarG, ZhuH. Comprehensive study of surface chemistry of MCM-41 using 29Si CP/MAS NMR, FTIR, pyridine-TPD, and TGA. J Phys ChemB. 1997;101(33):652531.MestreA, PiresJ, NogueiraJ, CarvalhoA. Activated carbons for the adsorption of ibuprofen. Carbon. 2007;45(10):197988.BaccarR, SarràM, BouzidJ, FekiM, BlánquezP. Removal of pharmaceutical compounds by activated carbon prepared from agricultural by-product. Chem Eng J. 2012;211:3107.GuedidiH, ReinertL, LévêqueJ-M, SonedaY, BellakhalN, DuclauxL. The effects of the surface oxidation of activated carbon, the solution pH and the temperature on adsorption of ibuprofen. Carbon. 2013;54:43243.NengN, MestreA, CarvalhoA, NogueiraJ. Cork-based activated carbons as supported adsorbent materials for trace level analysis of ibuprofen and clofibric acid in environmental and biological matrices. J Chromatogr A. 2011;1218(37):626370. doi: 10.1016/j.chroma.2011.07.025 21820664FooK, HameedB. Insights into the modeling of adsorption isotherm systems. Chem Eng J. 2010;156(1):210.TonghuanL, GuojianD, XiaojiangD, WangsuoW, YingY. Adsorptive features of polyacrylic acid hydrogel for UO2 2+. J Radioanal Nucl Chem. 2013;297(1):11925.BabonneauF, CamusL, SteunouN, RamilaA, Vallet-RegiM, editors. Encapsulation of ibuprofen in mesoporous silica: solid state NMR characterization MRS Proceedings; 2003: Cambridge Univ Press.MestreAS, PiresJ, NogueiraJMF, ParraJB, CarvalhoAP, AniaCO. Waste-derived activated carbons for removal of ibuprofen from solution: Role of surface chemistry and pore structure. Bioresource Technol. 2009;100:17206.SalemAttia TM, HuXL, YinDQ. Synthesized magnetic nanoparticles coated zeolite for the adsorption of pharmaceutical compounds from aqueous solution using batch and column studies. Chemosphere. 2013;93(9):207685. doi: 10.1016/j.chemosphere.2013.07.046 24074881ReijengaJC, IngelseBA, EveraertsFM. Thermodynamics of chiral selectivity in capillary electrophoresis: Separation of ibuprofen enantiomers with b-cyclodextrin. J Chromatogr. 1997;792:3718.ZhangL, LvJ, XuT, YangL, JiangX, LiQ. High efficiency removal and recovery of an endocrine disrupting compound–bisphenol AF from wastewaters. Sep Purif Technol. 2013;116:14553.Al-KhateebLA, ObaidAY, AsiriNA, SalamMA. Adsorption behavior of estrogenic compounds on carbon nanotubes from aqueous solutions: Kinetic and thermodynamic studies. J Indus Eng Chem. 2014;20:91624.MestreAS, PiresRA, ArosoI, FernandesEM, PintoML, ReisRL, et al Activated carbons prepared from industrial pre-treated cork: Sustainable adsorbents for pharmaceutical compounds removal. Chem Eng J. 2014;253:40817.MunozB, RamilaA, Perez-ParienteJ, DiazI, Vallet-RegiM. MCM-41 organic modification as drug delivery rate regulator. Chem Mater. 2003;15(2):5003.NakashimaH, OmaeK, SakaiT, YamazakiK, SakuraiH. Acute and subchronic inhalation toxicity of tetraethoxysilane (TEOS) in mice. Arch Toxicol. 1994;68(5):27783. 8085937KimY, BaeJ, ParkJ, SuhJ, LeeS, ParkH, et al Removal of 12 selected pharmaceuticals by granular mesoporous silica SBA-15 in aqueous phase. Chem Eng J. 2014;256:47585.KamarudinN, JalilA, TriwahyonoS, ArtikaV, SallehN, KarimA, et al Variation of the crystal growth of mesoporous silica nanoparticles and the evaluation to ibuprofen loading and release. J Colloid Interf Sci. 2014;421:613.SlowingII, TrewynBG, GiriS, LinYSY. Mesoporous silica nanoparticles for drug delivery and biosensing applications. Adv Funct Mater. 2007;17:122536.YangP, GaiS, LinJ. Functionalized mesoporous silica materials for controlled drug delivery. Chem Soc Rev. 2012;41(9):367998. doi: 10.1039/c2cs15308d 22441299