3.1. X-ray Diffraction (XRD)
The synthesized materials were investigated using the X-ray diffraction technique. In order to evaluate possible impacts on the formation of hydroxyapatite, caused by co-doping with europium and gadolinium, results of materials synthesized only with europium (HA-Eu(2%) and MSN/HA-Eu(2%)) will also be presented, for purposes of comparison.
The results evidenced reflection peaks in all the diffractograms that can be assigned to the hexagonal hydroxyapatite crystalline phase with spatial group P63/m, according to the card number 9-432 from Powder Diffraction File database (PDF2, International Centre of Diffraction Data—ICDD) as shown in
Figure 3. The characteristic peaks at 32°, 33°, and 34° Bragg angles (2θ) can be respectively attributed to the (211), (300), and (202) crystallographic planes, in agreement with the literature [
25]. Similar work described in the literature also showed the formation of HA phase, confirmed due to the appearance of characteristic patterns in the Bragg angles [
21,
26].
The MSN/HA-Eu(2%), MSN/HA-Eu(2%)-Gd(1%), MSN/HA-Eu(2%)/P(MAA) and MSN/HA-Eu(2%)-Gd(1%)/P(MAA) samples revealed a diffuse halo in the 15° to 35° 2θ range (
Figure 3, the broad peak centered around 22°), which can be assigned to the characteristic short-range periodicity of the (SiO4)
4−, a tetrahedral unit of amorphous silica nanoparticles suggesting that the synthesis of the nanocomposites was successfully achieved. This finding is in agreement with previous works in the literature [
19].
Figure 4 compares XRD results of the HA, HA-Eu(2%) and HA-Eu(2%)-Gd(1%) samples. The results show that the HA-Eu(2%) and HA-Eu(2%)-Gd(1%) samples showed a small peak at the Bragg angle 2θ of approximately 29º (
Figure 4a). This peak can be attributed to the presence of a second phase after the doping processes, confirmed by the Rietveld refinement, as shown in
Figure 4b–d and
Table 1. For HA-Eu(2%) and HA-Eu(2%)-Gd(1%) samples, it is observed that there was the formation of a secondary phase of Ca
2P
2O
7, which corresponds, respectively, to 8.3% and 25.5% of the synthesized material. Despite this, these materials are mainly composed of the hydroxyapatite phase, as shown in
Table 1. Therefore, the XRD result indicates the synthesis of crystalline materials with HA as the primary phase in the dopant materials.
The results show good agreement with the reference values (ICSD16742). A decrease in the unit cell volume is observed in all synthesized materials in relation to the reference value (530.139 Å
3). Different values of this parameter are found in the literature ranging from 525.580 Å
3 to 532.48 Å
3 [
27,
28]. Aldén and collaborators suggest that these variations can be attributed mainly to the presence of contaminants, such as carbonates and fluorides, which can occupy the positions of Ca and the phosphate group in the unitary cell of hydroxyapatite [
27]. Furthermore, the addition of rare earths into the hydroxyapatite structure resulted in reduced mean grain size (
Table 1). This may occur due to the decreasing of Ca-O distances and smaller ionic radius for rare earths resulting in a reduction of cell parameters as shown in
Table 1 [
29,
30,
31].
3.2. Fourier Transform Infrared Spectroscopy (FTIR)
Figure 5 shows the FTIR spectra for HA, HA-Eu(2%)-Gd(1%), MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/P(MAA). The spectra evidence the presence of hydroxyapatite in all the samples, due to the characteristic transmittance of the phosphate (1098–1037, 960, 603, 560, and 470 cm
−1) and hydroxyl stretching (3571 and 638 cm
−1) [
19,
26,
32]. The bands observed around 1454– 1412 cm
−1 in the HA and HA-Eu(2%)-Gd(1%) spectra can be attributed to the carbonate (CO
32−) vibrational modes. These ions may have originated due to the exposure of these samples to atmospheric CO
2 before the hydrothermal treatment or due to the decomposition of the CTAB during the calcination process [
21,
26].
FTIR spectra of MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)P(MAA) samples present the typical MSN vibrational modes. The characteristic transmittance bands from the angular deformation mode of Si-O-Si bonds can be found around 464 cm
−1 [
33,
34], the Si-O symmetrical stretching mode is located at 810 cm
−1 [
34], the Si-OH stretching mode at 961 cm
−1 [
35], and the a symmetric stretching of Si-O-Si bonds around 1075 a 1224 cm
-1 [
36]. The presence of the characteristic bands of MSN and HA in the FTIR spectra of the nanocomposites suggests that mesoporous silica can be formed in the hydroxyapatite nanoparticles, in agreement with the XRD results.
The presence of functional groups related to the polymeric phase is noted. The peaks around 3020 and 2450 cm
−1 (
Figure 4a) are attributed to the symmetrical and asymmetrical stretching modes of the -CH
2 and CH
3 groups, which arose due to the incorporation of the polymer phases [
10,
37]. In the deconvolution of region of 1750-1400 cm
-1 (
Figure 4b), it becomes evident a band around 1650 cm
-1 that can be attributed to the stretching modes of the carbonyl groups (C=O) present in the monomer [
38,
39]. The band between 1700 and 1715 cm
-1 can be related to the characteristic stretching modes of the carboxylate group, which is part of the MAA monomer used during polymerization, reinforcing the evidence of the presence of the polymeric phase in the hybrid [
19,
40]. These bands have already been observed by other works in studies involving nanoparticles to which polymeric phases have been incorporated [
10,
19].
Figure 6 shows the FTIR spectra obtained for the samples MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/P(MAA) to which the antitumor drug DOX were incorporated. The bands already observed for these materials prior to the incorporation of the drugs continued to be observed, indicating the maintenance of the phases of the materials even after the entire incorporation process. Furthermore, typical doxorubicin bands were identified in the sample spectra indicating that drug incorporation process was achieved. The bands at 3340 cm
-1 N–H (3525 cm
-1 ), O-H (2987 cm
-1), C–H (2897 cm
-1), C=O (1720 cm
-1), C=C ring (1690, 11581 and 1540 cm
-1), C–O–C (1071 cm
-1), C=H bend, C=C ring bend (800 cm
-1) can be interpreted as indicative of the presence of the drug in these samples [
41,
42].
3.4. Potential Zeta (ζ)
The zeta potential is used to characterize the surface charges of the particles and the values obtained are presented in
Table 2. The results obtained indicate a value around -17 mV for HA, which is attributed to the ionization of the hydroxyl groups (OH) present in the surface of the material, a result that agrees with the load presented by most of the hydroxyapatite ceramics in the literature [
46,
47]. After the doping process with Eu and Gd, there is no significant difference in the Zeta potential. This fact can be explained by the HA unit cell that has calcium in two positions, a more internal and another more external atom in relation to their structure. Therefore, the doping process may have occurred in the most internal positions of the unit cell, not being detected by this technique. The presence of silica in the system caused a decrease in the Zeta potential (-23.2 ± 0.23). Probably, this result may be related to the presence of negative charge on the surface of this material due to the presence of silanol groups characteristic of this material [
48]. After polymerization with methacrylic acid, an increase in the negative value of zeta potential is observed, probably due to the presence of carboxyl groups [
37]. Furthermore, the Zeta potential of the hybrid sample changed from -27.2 mV to -21.0 mV after the drug incorporation process. Possibly, the amino group present in doxorubicin may have contributed to the decrease in the negative charge of the MSN/HA-Eu(2%)-Gd(1%)/P(MAA)-DOX system, allowing us to suggest that DOXO was incorporated into the material. This result is in line with what was observed in the FTIR.
The zeta potential helps evaluate the stability of colloidal dispersions. In drug delivery systems, the zeta potential measurements help optimize the stability and performance of drug delivery systems. It assists in determining the surface charge and stability of these systems, which can affect their circulation time, interaction with cells, and drug release behavior. The ideal zeta potential value for a controlled drug delivery system can vary depending on various factors, including the specific formulation, the intended mode of action, and the target site within the body [
6]. However, in general, a moderately high zeta potential (either positive or negative) is often preferred for controlled drug delivery systems. In general, the margin that defines the stability of colloidal suspensions is -30 mV and +30 mV, where the higher the zeta potential values, in modulus, the more stable the suspensions. It is observed that HA and HA/Eu (2%)-Gd(1%) materials present values a little different from these in modulus, suggesting that these materials do not present high stability in colloidal suspensions. It is observed that the addition of silica increased the value in modulus of the zeta potential in the MSN/HA-Eu(2%)-Gd(1%) sample. Nevertheless, the value presented by this material is still a little far from the ideal value. After the polymerization process with P(MAA), the zeta potential value obtained was -27.2. This appears as a promising result and shows that the presence of the polymer in the system can minimize the tendency to agglomeration and, consequently, improve the stability of the material. In this sense, to determine the stability of nanoparticles in aqueous dispersion, measurements of zeta potential were carried out over 15 days.
Figure 8 presents the results of this study through the distribution curves of zeta potential measurements for the different synthesized systems.
Hydroxyapatite, when in aqueous dispersion, showed an altered surface charge throughout the studied period, indicating that the dispersion presents high instability over time, as expected for this system. In the results obtained for the HA-Eu(2%)-Gd(1%) system, it is observed that the nanoparticles doped with europium and gadolinium presented greater stability over 15 days in aqueous dispersion. The material remained stable with the same zeta potential distribution indicating that there was no significant change in the surface charge of the material during the analyzed period. In contrast, the MSN/HA-Eu(2%)-Gd(1%) nanocomposite showed some instability in the zeta potential value. Despite this, it is observed that after the polymerization process on the surface of this material, there was an improvement in the stability of these particles over 15 days. The hybrid system MSN/HA-Eu(2%)-Gd(1%)/P(MAA) maintained the zeta potential distribution over the 15 days. These data reinforce the need for a polymerization process to obtain a more stable material. The stability of this material is essential for a better application performance in biological systems.
3.6. N2 Adsorption
Nitrogen adsorption and desorption isotherms of the HA-Eu(2%)-Gd(1%), MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/P(MAA) samples are shown in
Figure 9. This technique was used to evaluate the pore structure of samples, since the performance of the nanomaterial in controlled drug delivery is largely influenced by its pore structure.
It is possible to observe that sample MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/P(MAA) samples exhibited type IV isotherms typical for mesoporous structures with constant cross-section porosity (cylindrical or hexagonal) and three-dimensionally ordered pores with lower or no blockage [
19,
46]. Four well-defined regions can be seen in these isotherms. Initially, in region I, a linear increase of the adsorbed volume happened due to the monolayer and multilayer adsorption on the surface of the material, including the surface of the mesoporous. In region II, the increase of the adsorbed volume at intermediate relative pressures was attributed to the capillary condensation phenomenon in the mesopores. In this point, the amount of adsorbed gas increased sharply with a small change in pressure due to the condensation of the adsorbate molecules below their vapor pressure. The adsorbate in the liquid state filled up the primary mesopores until the point in which the slope of the curve change, almost forming plateau. In region III, the increase in the adsorbed volume under high relative pressures was associated with the multilayer adsorption and/or to the condensation of the liquid in the secondary mesopores. Finally, in IV, the increase of the adsorbed volume corresponded to the filling of the voids between the particles that can be considered as porosity.
In sample HA-Eu(2%)-Gd(1%) it was possible to observe a sudden increase in the curves for high relative pressures. The abrupt increase in the adsorbed volume in regions of elevated relative pressure is characteristic of pores formed between the particle agglomerates (macroporosity).
The HA-Eu(2%)-Gd(1%) presented type-H3 hysteresis loop which is often observed with aggregates of plate-like particles that give rise to slit-shape pores [
32].The MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/P(MAA) samples exhibited H4 hysteresis loops, typical for mesoporous structures with constant cross-section porosity (cylindrical or hexagonal) and three- dimensionally ordered pores with lower or no blockage [
19,
49].
Furthermore, it is possible to observe that MSN/ HA/Eu(2%)-Gd(1%)/P(MAA) sample exhibited significant reductions in the adsorbed N2 volume at all relative pressures; this result indicates that the polymerization processes led to significant changes in the pore structure, indicating the presence of polymeric phase in the system, as evidenced by the previous characterization results.
The calculated texture and porosity parameters are presented in
Table 4. It is possible to observe a significant increase in surface area after the addition of silica to the system, which can be attributed to the intrinsic characteristics of this material, such as high surface area and pore volume. This feature is essential for nanosystems with an emphasis on drug delivery applications. The surface area obtained for the nanocomposite MSN/HA-Eu(2%)-Gd(1%) of 289.13 m
2.g
−1 is in agreement with a similar nanocomposite described in the literature, with application in drug delivery [
19,
22,
50].
It is observed that after polymerization with P(MAA), there was a reduction in the surface area and the pore volume of the MSN/HA-Eu(2%)-Gd(1%)/P(MAA). This can be ascribed to the hydrogel polymerization occurring both on the surface and inside the pore structure of the MSN/HA-Eu(2%)-Gd(1%) phase [
19,
51]. The pore size diameters remained unchanged after the hydrogel polymerization, since the BJH method calculates the average value of the nanoparticle diameter without distinguishing the amount of large or small pores present in the sample [
19,
37].
3.7. Photoluminescence Analysis
The room temperature PL spectra of the samples were measured in a wide wavelength range (300–850 nm) by employing an excitation laser line at , where no emissions were observed between 457 and 600 nm for all samples. Furthermore, as expected, the HA sample does not display significant luminescence in its whole spectrum. Thus, to present more detailed spectra, they were separated into the following two regions: (i) high-energy emissions (300-470 nm) with a total integrated PL intensity named I1, and (ii) low-energy emissions (550-870 nm) with a total integrated PL intensity named I2.
Figure 10 presents the PL spectra of all samples measured in the high-energy region (300-470 nm), where, except for the HA sample, the several observed peak emissions could be ascribed according to the indicated radiative transitions between the energy levels of Eu
3+ ions. No emissions from the levels of the Gd
3+ ions were identified. Since the laser line (λ
exc= 457 nm) has no enough energy to excite the high energy levels of the Eu
3+ ions, the PL emissions involving these states can only be explained by a complex phonon assisted and cooperative upconversion mechanism of nonradiative energy transfer involving dipole-dipole and dipole-quadrupole interactions between rare earth ions in the nearest neighbor positions, which can be enhanced by exchange coupling between these ions [
52,
53,
54,
55].
Specifically for the samples with two unlike trivalent rare earth ions, Eu
3+ and Gd
3+, (
HA-Eu(2%)-Gd(1%), MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/P(MAA) samples), this energy transfer mechanism can predominantly be governed by the dipole-quadrupole interaction [
53], mainly for the states of the larger rare earth ions that lie high in energy [
52].
Figure 11 presents the PL spectra of all samples measured in the low-energy region (550-850 nm), showing several peak emissions excluding the spectrum of the HA sample. As specified in this figure, these emissions were attributed according to the radiative transitions involving the levels of the
Eu3+ ions:
5D
1→
7F
J (J = 4, 5, 6) and
5D
0→
7F
J (J = 2, 3, 4, 5, 6). It is interesting to note that the characteristic emissions
5D
0→
7F
0 (~570–585 nm) and
5D
0→
7F
1 (~585–600 nm) of
Eu3+ ions [
56] were not observed. This effect can be attributed to the relatively high wavelength of the excitation laser line (
) which is activating the transitions indicated in the
Figure 11 and deactivating the transitions
5D
0→
7F
J (J = 0, 2).
It is also possible to notice in
Figure 11 that the addition of gadolinium in the HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%) samples has caused an enhancement in their PL intensities in comparison to the HA-Eu(2%) and MSN/HA-Eu(2%) samples, respectively. Xie and collaborators obtained a similar result, where the increase in Eu
3+ luminescence intensity could be observed by co-doping with Gd
3+. These authors obtained an increase of 60% to 120% in the emission intensity depending on the molar ratio of the europium and gadolinium content R(Eu/Gd), by explaining this luminescent enhancement phenomenon based on the already mentioned cooperative upconversion mechanism [
55] involving the adjacent
Eu3+ and Gd3+ ions.
Figure 12A presents an illustration of this cooperative upconversion mechanism in a Dieke's energy level scheme, where initially the laser line
can excite by photon absorption, involving phonons due to the nonresonance between the energies, the
5D
2 level of the
Eu3+ ions.
Then, a nonradiative thermalization assisted by phonons occurs from the
5D
2 to both the
5D
1 and
5D
0 levels of
Eu3+ ions, as indicated by the doted wave lines. In sequency, the cooperative upconversion mechanism of nonradiative energy transfer involving dipole-dipole and dipole-quadrupole interactions of the adjacent rare earth ions (
Eu3+ and Gd3+) takes place promoting excitation of the upper energy levels of the Eu3+ ions with consequent nonradiative energy transfer to the 5P
J levels of the
Gd3+ ions. Since no emissions from the gadolinium’s levels were identified in the high-energy region (300-470 nm) of PL spectra of samples (
Figure 10), it is concluded that a back nonradiative energy transfer occurs from the
Gd3+ to the Eu3+ ions. The nonradiative excitation of the upper levels of Gd3+ ions and the subsequent back energy transfer to the upper levels of the Eu3+ ions are indicated in the Figure 12A by the curved arrows. From the upper energy levels of the
Eu3+ ions, several radiative emissions were identified in the PL spectra shown in the Figure 10. Finally, a part of energy from the upper energy levels of the
Eu3+ ions is also returned by nonradiative transitions to their lower energy levels (
5D
2,
5D
1 and
5D
0), causing the enhancement in PL intensities observed in
Figure 11 ascribed to the transitions
5D
1→
7F
J (J = 4, 5, 6) and
5D
0→
7F
J (J = 2, 3, 4, 5, 6).
To quantify this complex cooperative upconversion mechanism of nonradiative energy transfer in the all samples with rare earth ions (
Eu3+ and Gd3+), it was carried out an analysis of the total integrated PL intensities I
1 and I
2 for the two regions: (i) I
1: high-energy emissions (300-470 nm) (
Figure 10), and (ii) I
2: low-energy emissions (550-870 nm) (
Figure 11).
Figure 12B(a) shows a comparison of the percentage of these total integrated PL intensities, E
1 and E
2, for these samples with the following values: (i) Eu
2O
3: I
1= 49.29% and I
2=50.71%; (ii) HA-Eu(2%): I
1= 71.68% and I
2=28.32%; (iii) HA-Eu(2%)-Gd(1%): I
1= 54.96% and I
2=45.04%; (iv) MSN/HA-Eu(2%): I
1= 61.95% and I
2=38.05%; (v) MSN/HA-Eu(2%)-Gd(1%): I
1= 45.09% and I
2=54.91%; (vi)
MSN/HA-Eu(2%)-Gd(1%)/P(MAA): I
1= 56.85% and I
2=43.15%. In Figure 14a the HA-Eu(2%) presents a substantial increase of I
1 and decrease of I
2 in comparison with the Eu
2O
3, evidencing the strong cooperative upconversion mechanism between the Eu
3+ ions in the nearest neighbor positions. The substantial increase of I
1 and decrease of I
2 the HA-Eu(2%) evidences its strong cooperative upconversion mechanism in comparison with the Eu
2O
3 but with a less efficient return of energy, by nonradiative transitions, from the upper energy levels to the lower energy levels of the
Eu3+ ions. In contrast with the HA-Eu(2%), the growth of I
2 and reduction of I
1 for the HA-Eu(2%)-Gd(1%) confirms that the inclusion of the Gd
3+ ions has favored the nonradiative back energy transfer to the lower energy levels of the
Eu3+ ions by enhancing their PL intensities in low-energy emissions (
Figure 11).
In comparison with the HA-Eu(2%), the MSN/HA-Eu(2%) samples shows increase of I
2 and decrease of I
1, evidencing its better efficiency nonradiative back energy transfer to the lower energy levels of the
Eu3+ ions. The MSN/HA-Eu(2%)-Gd(1%) sample, in contrast with the MSN/HA-Eu(2%), displays a strong increase of I
2 and decrease of I
1, again evidencing the improvement of the cooperative upconversion mechanism caused the presence of the Gd
3+ ions with consequent intensification of the PL intensities in low-energy emissions. Finally, in comparison with the MSN/HA-Eu(2%)-Gd(1%), the
MSN/HA-Eu(2%)-Gd(1%)/P(MAA) samples has decrease of I2 and increase of I1, which can be attributed the functionalization of the nanocomposite with the P(MAA).
In agreement with the results of the
Figure 12B(a), the
Figure 12B(b) presents a comparison of the ratio I
2/I
1, normalized for HA-Eu(2%) sample, in order to provide a undoubted understanding about the enhancement of the PL intensities in low-energy emissions (550-870 nm) caused by the gadolinium. The HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%) shows strong growth of the ratio I
2/I
1 , from 1.00 to 2.07 and from 1.55 to 3.08, in comparison with HA-Eu(2%) and MSN/HA-Eu(2%), respectively. Once more, this result confirms the high efficiency cooperative upconversion mechanism of nonradiative energy transfer involving multipolar interactions between Gd
3+ and Gd
3+ ions, with the subsequent partial nonradiative back energy transfer from the upper energy levels of the
Eu3+ ions to their lower energy levels (
5D
2,
5D
1 and
5D
0). Furthermore, the functionalization of the nanocomposite with the P(MAA) in the MSN/HA-Eu(2%)-Gd(1%)/P(MAA) sample has caused the weakening of the cooperative upconversion mechanism, reducing the ratio I2/I1 from 3.08 to 1.92 compared to the MSN/HA-Eu(2%)-Gd(1%), which can be attributed to an energy loss involving the vibrational states of the P(MAA). Nonetheless, this value 1.92 is still very higher (almost double) of the value 1.00 observed for the HA-Eu(2%) sample, evidencing that the nonradiative multipolar energy transfer between Eu3+ and Gd3+ ions is maintained very efficient.
3.9. Drug Loading and Releasing Results
In this study, we evaluated the incorporation and release of the antitumor drug doxorubicin in the luminescent materials MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/P (MAA), using UV-VIS spectroscopy. The incorporation efficiency of the antitumor drug doxorubicin in these materials was esti mated. The results obtained were 71% for the MSN/HA-Eu(2%)-Gd(1%) sample and 95% for the MSN/HA-Eu(2%)-Gd(1%) / P(MAA) ). It is possible to notice a higher incorporation efficiency (E.I.) of the luminescent hybrid when compared to the MSN/HA-Eu(2%)-Gd(1%) sample. This behavior may be related to the cationic character presented by the doxorubicin molecule (pKa 8.6) which is different from the anionic character presented by the carboxylic group present in methacrylic acid (pKa 5.5-6.0). In this sense, in incorporation conditions with pH around 7.0, these characteristics favor the electrostatic interaction between DOX and the MAA hydrogel, leading to a greater efficiency of drug incorporation [
59].
Drug release assays were performed as a function of time for 168 h at two different pH's, 5 and 7 (
Figure 14). This test is essential to evaluate the potential of this material for application in controlled drug release, based on the difference in pH that exists between normal and tumor cells. The MSN/HA-Eu(2%)-Gd(1%) sample showed a total release percentage of 22% at pH 5 and 16% at pH 7, indicating that the difference in the medium does not affect significatively the release profile. In the other hand, for the luminescent hybrid system MSN/HA-Eu(2% )- Gd(1%)/P(MAA), it is possible to notice significant changes in the total DOX release values, 27% at pH 5 and 1.5% at pH 7. These results indicate that the release of DOX from MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/P(MAA) materials is pH dependent. of the medium, being favored at pH 5. Despite this, a greater responsiveness to the pH of the hybrid system is evident, as expected.
Possibly, the greater amount of DOX released in the tests performed at pH 5 occurs due to the protonation of the carboxyl groups present in the hydrogel, which results in lower electrostatic interactions between it and the drug. At pH 7, smaller amounts of the drug were released and could probably be related to the strong electrostatic interaction between the DOX molecules and the MAA segments present in the hybrid materials [
60]. On the other hand, it is observed that not all incorporated material was released for the two studied samples. Possibly there may be the intermolecular interactions between the materials and the drug occurs through.
The study of dynamic swelling and ionic gel equilibrium in fluids used in controlled drug delivery studies are important to understand the diffusion process [
61]. To perform the kinetic studies, the data obtained were adjusted to the Korsmeyer-Peppas equation, represented by the following expression [
62,
63,
64].
where Mt corresponds to drug release at time t, Mo represents the amount of drug initially incorporated into the material, and n is the release index, indicative of the mechanism of drug release, and k is a kinetic constant characteristic of the drug-carrier system.
Peppas (1985) used the value of n to characterize the different release mechanisms, having reached values of n = 0.5 for Fickian diffusion and values of n between 0.5 and 1.0, for mass transfer according to a non-Fickian model. This model is generally used to analyze the release of polymeric dosage forms, when the release mechanism is not well known or when more than one type of release may be involved [
65].
Table 5 presents the results of the release kinetics of the drug doxorubicin obtained for the samples MSN/HA-Eu(2%)-Gd(1%) and MSN/HA-Eu(2%)-Gd(1%)/ P (MAA), as well as the correlation coefficients (R
2), used to assess the accuracy of the fit.
The regressions that originated the release kinetics data are illustrated in
Figure 15. Analyzing the presented curves, it is possible to perceive different inclinations for all samples. Therefore, the graphs were separated by regions and the values (R
2) obtained for all samples were greater than 0.97, indicating that the linear correlation coefficient was ideal.
It can be inferred, for the hybrid system MSN/HA-Eu(2%)-Gd(1%)-P(MAA), that the release model predominantly follows the Fickian diffusion, with n ≤ 0.5 [
66]. Note, for this system, that the test conducted at pH 5 (region I), initially presented a release model that does not follow the Fickian diffusion. This behavior may be related to the burst release, attributed to the desorption of doxorubicin molecules located on the surface of the synthesized materials [
32]. Sousa et al. observed that the MCM-41-HA system showed a fast delivery rate during the first 2 hours of the test, releasing about 18% of the incorporated model drug. After that, the system showed a slower rate, with a cumulative release of approximately 45% after 160 h of testing [
22].
Furthermore, when analyzing the drug release kinetics for the MSN/HA-Eu(2%)-Gd(1%) nanocomposite, it is observed that for most of the test time, this system does not follow Fickian diffusion. This result demonstrates a more uncontrolled release of the system that does not contain the polymeric phase and reinforces the previously discussed theory that the hybrid system presents a more controlled release of the drug. In addition, it is noted that the MSN/HA-Eu(2%)-Gd(1%)/P(MAA) hybrid system has a higher rate of doxorubicin drug release than MSN/HA-Eu(2%) - Gd(1 %) in pH 5, associated to the visibly greater kinetic constant K, that is related with the drug release rate. In other words, this result shows that the presence of the polymer in the system affects the release kinetics. The influence of the pH of the medium on the release kinetics is also observed, corroborating what was previously discussed.
This result is in line with what was expected in this work, whose objective was to develop a system capable of delivering antitumor drugs at a more acidic pH (tumor region), and not reaching healthy cells, which have a more alkaline pH. This fact demonstrates that the system developed in this work is promising for targeted delivery of antitumoral drug in a targeted and controlled manner from responsive polymers.
The greater control of drug release by the MSN/HA-Eu(2%)-Gd(1%)/P(MAA) luminescent hybrid system, when subjected to pH changes, may indicate that the material developed in this work has the potential to targeted and controlled drug delivery. This release behavior is desirable for cancer treatment, as most of the drug remains encapsulated in the nanocarrier during circulation and when it reaches the acidic tumor tissue, the low pH triggers drug release.