Synthesis and characterization
MCM-41 MSNs were synthesized through sol-gel method using TEOS and CTAB as silica resource and template surfactant, respectively. The structure of MSNs was evaluated by X-ray diffraction in low angles. According to the diffraction pattern of Fig. 1, three distinct characteristic peaks at 2.2°, 3.9°, 4.6° and a mild peak at 6.1° are observed for MSNs which are attributed to the reflections 100, 110, 200 and 210 planes of the hexagonal ordering [35]. TEM imaging was also used to depict the mesoporous structure of MSNs. Figure 2a and 2b show the formed mesopores of MSNs in two different magnifications. It can be clearly observed the well-ordered porous channels with honeycomb-like hexagonal structures, that verifies the formation of MSNs and X-ray diffraction analysis results.
In order to identify the pores of MSNs and determine the surface area, nitrogen gas absorption and desorption isotherm was obtained based on Brunauer–Emmett–Teller (BET) theory. Figure 3a shows that the nitrogen adsorption-desorption isotherm follows the type IV according to the classification by IUPAC. Also, there is a hysteresis loop in high ratios of P/P0. This result indicates the formation of mesoporous structure and the presence of very narrow and capillary pores of MSNs. According to BET method, the surface area of MSNs was calculated to be about 1023 m2g− 1. In addition, the average diameter of the pores was obtained 3.50 nm, from BJH study (Fig. 3b).
To covering MSNs with chitosan, GPTMS was used as a crosslinker. In the acidic condition, trimethoxysilane groups of GPTMS hydrolyzed to form silanol groups which condensed in the presence of silica network of MSNs to form a cross-linked structure through Si-O-Si covalent bonds [36, 37]. On the other hand, oxirane group of GPTMS reacted with amino groups of chitosan by the acid-catalyzed addition reaction and chitosan coated MSNs were prepared [38]. To load Que in nanopores, the pH-responsive property of chitosan was utilized. In this way, reducing the pH to the acidic ranges, Cs-MSNs swelled due to protonation of amino groups in chitosan, and Que entered in the pores of MSNs. Readjusting the pH to neutral ranges caused shrinkage in Cs-MSNs structure, capping of the pores due to deprotonating of amino groups of chitosan and entrapment of Que in the pores of MSNs (Scheme. 1).
In addition, the morphology and size distribution of the Que@Cs-MSNs were investigated by FESEM analysis (Fig. 4). The results indicated that particles have nearly spherical shapes and a narrow size distribution with average diameter of 185.09 ± 4.21 nm. Previous studies have shown that nanoparticles with dimensions in the range of 100 to 200 nm can accumulate in the tumor site due to passive targeting caused by enhanced permeability and retention effect (EPR) [39]. Therefore, the dimensions obtained for nanoparticles here can be useful for passive targeting of cancer cells.
To confirm the loading of Que into Cs-MSNs, FTIR analysis was performed on pure Que, Cs-MSNs and Que@Cs-MSNs samples. According to FTIR spectrum of Que in Fig. 5, the characteristic band around 3406 cm− 1 was attributed to –OH groups, 1666 cm− 1 to C = O stretching, 1611 cm− 1 and 1521 cm− 1 to aromatic C = C stretching vibrations. The C–O stretching of oxygen in the ring was found at 1263 cm− 1 beside other C–O bands at 1168 cm− 1 and 1013 cm− 1. Also the bands at 1450 cm− 1, 1382 cm− 1, 1319 cm− 1, 864 cm− 1 and 823 cm− 1 belong to C-H bending vibrations [40–42]. Comparing the spectrum of Que@Cs-MSNs with the pure Que implies that the bands at 1658 cm− 1, 1512 cm− 1, 1357 cm− 1, 1317 cm− 1 and 819 cm− 1 belong to Que which confirms the successful loading of Que in Cs-MSNs carriers. The existence of some shifts in the spectrum of Que in nanoparticles can be due to its trapping in the Cs-MSNs pores and overlapping of some peaks by chitosan bands.
In-vitro release behavior
Fig. 6 shows the cumulative release of Que from the Que@Cs-MSNs at pH 5 and 7.4. In addition, to study the effect of pH-sensitive chitosan coating on the release behavior, the release of Que from MSNs without chitosan coating was also investigated and shown in the Fig. 6. According to the results, the release of Que from nanoparticles is burst in the first 5 h, followed by a controlled manner at both pH, but the amount of released drug depends on presence of chitosan coating and pH value. The amount of Que released from Que@Cs-MSNs at pH 7.4 reached about 23.85 % after 36 h, while this amount for uncoated MSNs was about 49.60 %. Indeed, at pH range above 6.5, chitosan caps the pores of MSNs as a gatekeeper due to deprotonation of amine groups and shrinkage, so prevents the release of more Que from MSNs. In acidic pH, although the amount of Que released from uncoated MSNs is equal to 57.38 % after 36 h, this amount reached to 46.02 % for Que@Cs-MSNs. Owing to the pH-sensitive property of chitosan, at lower pH, amine groups are protonated, cause the structure of chitosan to swell and the gate of the pores are opened for the release of Que. This result has also been reported in similar studies [43, 33]. This property is beneficial for drug delivery to the tumor sites. Since the rate of glycolysis is higher in the tumor environment, it has more acidity and the pH value in tumor tissue reaches less than 6.8, while the pH of normal tissues and blood is about 7.4. Therefore, the drug cargo inside the MSNs is protected without leakage into the blood and damage to normal cells, and when it reaches the acidic tumor site, the drug is released by uncapping the holes of MSNs and dissolution of chitosan.
Release kinetic mechanisms
Mathematical kinetic models can be applied to better understand the mechanism and predict the release behavior of drugs in in-vivo conditions using the data obtained from in-vitro studies. Herein, the kinetic mechanism of Que release was evaluated by the Zero-order, First-order, Higuchi and Kormeyer-Peppas mathematical models.
Zero order model
The Zero-order kinetic model is an ideal model for drug release. For drug delivery systems that follow this model, the release rate is independent of drug concentration. In other words, drug release is constant throughout the delivery time in the blood flow [44]. The Zero-order kinetic can be represented by the following equation.
$${Q}_{t}={Q}_{0}+{k}_{0}t$$
1
Where \({Q}_{0}\) is the initial concentration of drug in release medium which is usually zero, \({Q}_{t}\) is released drug concentration at time \(t\) and \({k}_{0}\) is the constant of Zero-order equation.
First-order model
This model is presented based on Fick's law for diffusion and the release rate is dependent on drug concentration. Accordingly, due concentration gradient inside the carrier and outside of it, the drug diffuses to the release medium and as a result, its concentration inside the carrier decreases logarithmically with time [45]. The simplified equation of this model is as follows
$$\text{log}{C}_{t}=\text{log}{C}_{0}-({k}_{1}t/2.303)$$
2
Where, \({C}_{0}\) is the initial amount of drug in the carrier, \({C}_{t}\) is the amount of drug remaining in the carrier at time t and \({k}_{1}\) is the First-order constant.
Higuchi model
Higuchi developed a widely used model for the release of hydrophilic and hydrophobic compounds loaded in solid and semi-solid matrices with various geometries based on Fick's first law [46]. In Higuchi's simplified model, the amount of drug released () is related to the root of time and is Higuchi release constant.
$${Q}_{t}={k}_{H}\times \sqrt{t}$$
3
Korsmeyer-Peppas model
Korsmeyer and Peppas presented an empirical kinetic model for release mechanisms, both for Fickian and non-Fickian diffusion phenomena [47]. In this model according to Eq. 4, is the fraction of drug released at time, is Korsmeyer-Peppas constant which depends on structural and geometric characteristics of the system and is release exponent. The values obtained for predict the release mechanism, when the release mechanism corresponds to the Fickian diffusion and non-Fickian or anomalous transport.
$$\text{log}{Q}_{t}/{Q}_{\infty }=\text{log}{k}_{KP}+n\text{log}t$$
4
The fitting of the experimental data of Que release from Que@MSNs and Que@Cs-MSNs on these kinetic models are shown in the graph of Fig. 7, and the values obtained for the constants along with the coefficient of determination (R2) are given in Table. 1. The R2 value being closer to 1 means a better fit of the data with the model, and these values are higher for the Higuchi model in the examined samples. Therefore, diffusion phenomena are involved in the release of Que. On the other hand, the comparison of \(n\) values obtained from the Korsmeyer-Peppas model show that for Que@MSNs, the \(n\) values are less than 0.5 at both pH of 5 and 7.4, that means the release kinetic follow the Fickian diffusion. However, \(n\) values are 0.8850 and 0.6863 for Que@Cs-MSNs at pH of 5 and 7.4, respectively. These results indicate that in the presence of chitosan coating, the release kinetics of Que from MSNs changes and follows the non-Fickian or anomalous transport model. In this case, Que molecules are released both due to the diffusion from the MSNs matrix and the swelling or dissolution of the chitosan coating.
Table 1
R2 values and kinetic constants obtained from data fitting with mathematical models
Model | coefficients | Que@MSNs pH = 5 | Que@MSNs pH = 7.4 | QUE@Cs-MSNs pH = 5 | QUE@Cs-MSNs pH = 7.4 |
Zero order | k0 R2 | 1.291 0.9070 | 1.102 0.8688 | 1.290 0.8509 | 0.5805 0.8902 |
First order | k1 R2 | 0.0213 0.9499 | 0.0166 0.9140 | 0.0181 0.8904 | 0.0068 0.9103 |
Korsmeyer-Peppas | kKP n R2 | 1.080 0.4646 0.9613 | 1.002 0.4860 0.9274 | 0.4661 0.8850 0.8811 | 0.4107 0.6863 0.8490 |
Higuchi | kH R2 | 9.413 0.9745 | 8.148 0.9594 | 9.580 0.9484 | 4.255 0.9659 |
Cytotoxicity study
To study the cytotoxicity, G-292 osteosarcoma cells were treated with different concentrations of Que@Cs-MSNs (0-450 µM) and free Que for 24 and 72 h and evaluated by MTT assay. As shown in Fig. 8, after 24 h of treatment, in the concentration range of 50–150 µM, cell viability is above 88% for both Que@Cs-MSNs and free Que, and cytotoxicity is not significant. However, at higher concentrations, significant cytotoxicity is observed for both formulations compared to the control group. Therefore, cytotoxicity is concentration dependent but free Que and Que@Cs-MSNs does not exhibit a significant cytotoxicity difference up to 450 µM. In 72 h treatments, however, at the concentration of 100 µM and above, Que@Cs-MSNs show significant higher cytotoxicity. In addition, the dependence of cytotoxicity on the concentration after 72 h is more significant so that Que@Cs-MSNs and Que show IC50 values of 280.30 µM and 494.7 µM, respectively. These results demonstrate that firstly, the cytotoxicity of both formulations increase with enhancing incubation time; Secondly, Que loading within Cs-MSNs improves its anticancer property against G-292 osteosarcoma cells. The increased cytotoxicity of Que@Cs-MSNs after 72 h of treatment can be due to the fact that nanoparticles, by improving the limitations of Que caused by low solubility and instability, using specific cell uptake pathways, discharge the drug cargo in a controlled manner into the cells. Meanwhile, due to its hydrophobic nature, Que molecules may undergo agglomeration after 72 h and their cellular uptake will decrease compared to Que@Cs-MSNs.
In addition to the cytotoxicity assessment, optical microscope images of untreated G-292 cells and the cells treated with Que and Que@Cs-MSNs at IC50 concentration of Que@Cs-MSNs were prepared. According to the images shown in Fig. 9 untreated cells have healthy morphology, intact membrane, high density and adherent together on the surface of the plate. However, in cells treated with formulations, cell morphology is completely changed and membrane shrinkage due to the condensation of the cytoplasm is visible. Moreover, the density of the cells is decreased, the adjacent cells are not adhesive to each other and to the surface of the plate. These changes can be caused by the apoptosis induction in the treated cells. Also, the above changes are more significant in cells treated with Que@Cs-MSNs than free Que. Therefore, Que@Cs-MSNs inhibits G-292 cancer cells more effectively probably through apoptosis induction pathways, which is consistent with the results obtained from MTT assay.