Michael
Fischer
ab
aCrystallography & Geomaterials Research, Faculty of Geosciences, University of Bremen, Klagenfurter Straße 2-4, 28359 Bremen, Germany. E-mail: michael.fischer@uni-bremen.de
bBremen Center for Computational Materials Science and MAPEX Center for Materials and Processes, University of Bremen, 28359 Bremen, Germany
First published on 10th June 2024
Zeolites have been proposed as carrier materials for the encapsulation and controlled release of the anticancer drug 5-fluorouracil (5-FU). Besides, they could also find use in the adsorption-based removal of 5-FU from water, for example in the treatment of hospital effluents. In the present work, dispersion-corrected density functional theory (DFT) calculations and DFT-based ab initio molecular dynamics (AIMD) simulations are employed to study the interaction of 5-FU with faujasite-type zeolites having different Si/Al ratios. Comparing distinct local arrangements of Al atoms and charge-balancing protons, it is evaluated to what extent “multi-site” interactions, i.e., interactions of 5-FU with more than one proton, affect the adsorption energy. While the most pronounced increase in interaction strength occurs when moving from an all-silica zeolite to a protonic zeolite having a single proton in one twelve-membered ring, a significant additional stabilisation arises if a second framework proton is present in the same ring. Typically, several hydrogen bonds are formed between 5-FU and protonic zeolite frameworks, with 5-FU simultaneously acting as donor and acceptor. AIMD simulations confirm the stability of these hydrogen bonds at room temperature in the absence of water. Additionally, infrared spectra are predicted for selected low-energy configurations in order to facilitate an experimental identification of different bonding environments. AIMD simulations probing the competitive adsorption of 5-FU and water show that the high affinity of water to the framework protons causes a breaking of hydrogen bonds, framework deprotonation, and a displacement of 5-FU from its initial position. Exposure of dehydrated 5-FU@FAU composites to water might thus be a useful approach to trigger 5-FU release in drug delivery applications.
Fig. 1 a) Molecular structure and atom labelling scheme of 5-fluorouracil. The labelling follows earlier work.15 b) Schematic representation of the FAU framework. Sodalite cages are shown in blue, double six-membered rings in red. “Supercages” connected by twelve-membered rings (12MRs) are transparent. c) Labelling scheme of proton positions within a single 12MR. The shown example corresponds to a _1_3 arrangement. Yellow = Si, turquoise = Al, red = O, white = H. |
Zeolites are crystalline inorganic materials whose structures are built up by a framework of corner-sharing TO4 tetrahedra.16 Most natural zeolites and the majority of synthetic zeolites have an aluminosilicate composition (T = Si, Al), resulting in a negative framework charge that is balanced by extra-framework cations or framework protons. However, neutral-framework all-silica zeolites and zeolite-like materials having other elements on the T sites have also been reported.17,18 The presence of accessible cavities or channels is a key feature of zeolite crystal structures, and the arising intrinsic porosity is pivotal for their large-scale use in catalysis, ion exchange, and separation processes.17 Besides these established applications, zeolites have also received considerable interest as potential drug delivery materials, along with other possible applications in medicine.19 With regard to the specific task of 5-FU delivery, Datt et al. studied the storage and release of 5-FU from zeolite Y samples (FAU framework type16) having Si/Al ratios of 2.5, 15, and 30.20 They observed similar loadings for all three samples, on the order of 0.1 g 5-FU per gram of zeolite, but rather different release characteristics, with the most Al-rich sample (Si/Al = 2.5) releasing only a small fraction of the adsorbed 5-FU under physiological conditions. In contrast, the more Si-rich samples released more than 50% within a few minutes. 27Al-NMR spectroscopy provided evidence for a binding of 5-FU to extra-framework aluminium (EFAL) sites, which were most abundant in the most Al-rich sample.
Subsequently, a number of other authors studied the 5-FU delivery properties of various zeolites.21–27 Neves, Baltazar and co-workers investigated sodium-exchanged zeolite Y (NaY) and compared it to other systems including Linde type L and Linde type A zeolites (LTL and LTA framework, respectively) and ZSM-5 (MFI framework).21,24 They also demonstrated the cytotoxicity of the 5-FU@zeolite composites to different cancer cell lines in cell viability studies, whereas the pristine zeolites were shown to be non-cytotoxic. Spanakis et al. compared the 5-FU release behaviour of zeolite NaX (FAU framework) and zeolite beta (BEA framework) and employed molecular dynamics (MD) simulations using a classical force field (FF) to explain the observed differences.22 Apart from commercially available synthetic zeolites, which can be relatively expensive, zeolites synthesised from inexpensive and widespread starting materials (metakaolin)25 and natural zeolites27 were also investigated as 5-FU carrier materials. Another avenue of research is constituted by the preparation and characterisation of functionalised zeolites and zeolite-containing nanocomposites. For example, it was shown that surface functionalisation of synthetic zeolites can enhance the bioavailability of encapsulated 5-FU.28 Magnetic nanocomposites of zeolites and Fe3O4 nanoparticles are of interest as they could be directed to the location of the tumour through application of an external magnetic field.29,30 Furthermore, the use of silver-exchanged zeolites as 5-FU carriers could help to prevent microbial infections in the tumour microenvironment due to the antimicrobial activity of Ag+ cations.31 Looking beyond zeolites to other crystalline porous materials, metal–organic frameworks (MOFs) and MOF-containing nanocomposites have also been investigated as materials permitting controlled 5-FU delivery.32–34
Apart from their potential use as 5-FU carrier materials, zeolites could also be employed as adsorbents for the selective removal of 5-FU from wastewaters. While little focus has so far been placed on 5-FU in this context, a number of experimental investigations have dealt with the removal of other pharmaceutically active ingredients from wastewaters using hydrophobic high-silica zeolites.35–44 Although 5-FU is not commonly ranked among the pharmaceuticals of most concern with regard to their environmental impact,45,46 it has been identified as a priority cytostatic drug for environmental monitoring studies due to its reported persistence and ecotoxicity.47 In addition to occasional reports of low levels of 5-FU (< 10 ng L−1) in wastewater treatment plant influents,48 5-FU has been detected primarily in hospital effluents,49 in at least one case in high concentrations exceeding 100 μg L−1.50 The toxicity (including genotoxicity) of 5-FU, alone and in mixtures, has been investigated for various aquatic organisms, and estimations of the arising environmental hazard potential generally agree that its potential environmental impact is non-negligible.51–55 This highlights the need to develop efficient removal strategies from contaminated waters, especially hospital effluents. Among other technological options, adsorption-based processes using activated carbons have been proposed for 5-FU removal.56,57 While they are more expensive to produce than activated carbons, hydrophobic zeolites might constitute an interesting alternative due to the reduced pore blockage by co-adsorbed natural organic matter38 and good adsorbent regeneration properties.41 Due to their narrow pore size distribution, they are especially well suited for the highly efficient removal of molecules that fit well into the pores.
In addition to its potential relevance for different applications, highlighted above, 5-FU adsorption in zeolites is also interesting from a more fundamental point of view. Due to the presence of two negatively polarised carbonyl oxygen atoms and two amine groups (Fig. 1), 5-FU can act as hydrogen bond donor and acceptor, and a stabilisation in the zeolite pores through various hydrogen bonds can be envisaged, especially for protonic zeolites having (at least) two framework protons in relatively close proximity. Experimentally, the nature and extent of such “multi-site” interactions is rather difficult to probe, not least due to the important influence of defects on the adsorption properties, as highlighted above when discussing the work of Datt et al.20 In contrast, atomistic modelling methods are ideally suited to investigate such interactions in a systematic fashion, giving access to adsorption energies, low-energy structures, and derived quantities like vibrational spectra. Molecular dynamics simulations allow to study the dynamic behaviour of the host-guest system. It is clear that the use of idealised models of perfect zeolite crystals constitutes a certain limitation, as not all features that may play a role in real materials (defects, intrinsic mesoporosity, …) are incorporated in the modelling. So far, atomistic modelling methods have been employed to study the adsorption of pharmaceuticals in zeolites, both in the context of “screening” studies addressing a range of compounds in different zeolites with computationally efficient force field (FF) methods58,59 and in studies of individual compounds using computationally more expensive density functional theory (DFT) calculations.37,60–63 The specific topic of 5-FU adsorption in zeolites has so far received limited attention from the computational perspective, the previously mentioned study of 5-FU diffusion by Spanakis et al. being an exception.22 However, it is worth mentioning that dispersion-corrected DFT has been used to investigate the adsorption of 5-FU at a hydroxylated cristobalite surface.64 That work highlighted the important role of hydrogen bonds in stabilising the surface-adsorbed 5-FU molecule. Moreover, a DFT-based composite electronic structure method was employed to shed light on the local structure and vibrational properties of 5-FU adsorbed in HKUST-1, a copper-containing MOF.34
In the present work, dispersion-corrected DFT calculations are employed to investigate the adsorption of 5-FU in FAU-type zeolites, considering aluminosilicate models having different amounts of Al atoms and charge-balancing framework protons as well as a reference all-silica FAU model. Focus is placed on FAU for various reasons: first, FAU-type zeolites have been addressed in several previous drug delivery investigations, as highlighted above. Second, they are available across the whole range of Si/Al ratios (FAU type systems with Si/Al < 1.5 are commonly designated “zeolite X”, those with Si/Al 1.5 are dubbed “zeolite Y”), whereas many other framework types are synthetically accessible only within a certain range. Finally, the relatively simple structure, with only one non-equivalent T site, renders the generation of models having different proton arrangements relatively straightforward. The scope of the present study remains restricted to protonic, rather than cationic zeolites because the propensity of metal cations to associate with 6MRs of the FAU structure65 renders interactions of one 5-FU molecule with several cations unlikely, except, potentially, at high cation concentrations (approaching a Si/Al ratio of 1). The aims can be summarised as follows:
• The possible extent of “multi-site” interactions is probed through calculations for FAU models with different proton arrangements. In addition to comparing adsorption energies, an analysis of the low-energy configurations gives insights into the role of hydrogen bonds.
• A calculation of infrared (IR) spectra of adsorbed 5-FU is used to identify features that would allow for the experimental detection of different bonding situations using IR spectroscopy.
• Static calculations are complemented by DFT-based AIMD calculations, which allow to assess the lability/stability of the hydrogen bonds at room temperature. Further AIMD simulations are employed to investigate the competitive adsorption of 5-FU and water.
It is clear that the model systems employed in the present work are somewhat simplified, imposing specific arrangements of Al atoms and protons in the zeolite structure, rather than accounting for the coexistence of various local arrangements in real samples. However, the detailed insights obtained through a systematic comparison of different local environments, which bring about different bonding situations, do not only enhance the fundamental understanding, but they might also be relevant for the development of new zeolite applications.
Label | Unit cell composition | Si/Al ratio | Proton arrangement | a/Å |
---|---|---|---|---|
SiO2-FAU | Si192O384 | Infinity | — | 24.227 |
H-FAU_1H | HAlSi191O384 | 191 | 1 H in one 12MR | 24.237 |
H-FAU_2H_1_3 | H2Al2Si190O384 | 95 | 2 H in one 12MR, located at next-nearest O1 atoms | 24.248 |
H-FAU_2H_1_4 | H2Al2Si190O384 | 95 | 2 H in one 12MR, located at opposite O1 atoms | 24.248 |
H-FAU_3H | H3Al3Si189O384 | 63 | 3 H in one 12MR, located at next-nearest O1 atoms | 24.258 |
H-FAU_16H | H16Al16Si176O384 | 11 | 1 H in each 12MR | 24.398 |
H-FAU_32H_1_3 | H32Al32Si160O384 | 5 | 2 H in each 12MR, located at next-nearest O1 atoms | 24.579 |
H-FAU_32H_1_4 | H32Al32Si160O384 | 5 | 2 H in each 12MR, located at opposite O1 atoms | 24.566 |
The structure of 5-FU was taken from the PubChem database.1 While a number of tautomers exist, the diketo form visualised in Fig. 1a is the most stable one, which accounts for >99.9% of the tautomer population at 375 K (evaporation temperature).70 The structure of 5-FU was initially optimised using PCFF parameters.71,72 Different adsorption configurations of 5-FU in the pores of the FAU models were generated using FF-based Monte Carlo (MC) simulations, making use of the Sorption module integrated in the Materials Studio suite (DS BIOVIA).73 Details of the MC calculations are provided in the ESI,† section S1.1. For each adsorbent model, between 3 and 10 distinct, pre-optimised 5-FU@FAU configurations were used as starting points for DFT optimisations, which are described in the following subsection. For AIMD simulations including co-adsorbed water molecules, which were performed only for a subset of FAU models, a hydration shell of 16 H2O molecules was generated using a FF-based simulated annealing procedure as implemented in the Adsorption Locator module of Materials Studio. Details of these simulations are presented in the ESI,† section S1.2.
All FF-based simulations used PCFF parameters to describe the guest molecules (5-FU, water),71,72 and partial charges and Lennard-Jones parameters reported by Emami et al. to represent framework Si and O atoms.74 Since that work does not include parameters for framework protons, Al atoms, and O atoms of the AlO4 tetrahedra, partial charges for these atoms were derived from the DFT electrostatic potential of H-FAU_1H. The REPEAT method was used to compute charges,75 which were then scaled and slightly adjusted. The whole procedure is described in more detail in the ESI,† section S1.3. Lennard-Jones parameters of H, Al, and O(AlO4) atoms were estimated using parameters of similar atom types included in the work of Emami et al. All partial charges and Lennard-Jones parameters are compiled in Table S3 of the ESI.†
ΔEads = EDFT,5FU@FAU − EDFT,FAU − EDFT,5FU |
All atomic positions were relaxed in structure optimisations, which used a BFGS optimiser, imposing a maximal gradient of 5 × 10−6 Ha a0−1 and a maximal displacement between optimisation steps of 2.5 × 10−5a0 (a0 = 0.529177 Å) as convergence criteria. Calculations of the vibrational modes and IR intensities employed the finite displacement method, using a partial Hessian approach in which only displacements of the atoms of 5-FU and (where relevant) of framework protons and framework oxygen atoms participating in hydrogen bonds were considered. Structure visualisations were prepared using VESTA.82 The analysis of vibrational properties made use of MOLDEN, using a broadening with a full width at half maximum of 10 cm−1 for plots of IR spectra.83
All AIMD simulations were carried out for a temperature of 298 K in the NVT ensemble, using a Nosé–Hoover thermostat84,85 with a timestep of 0.5 fs and a time constant of 50 fs. Three types of tasks were tackled with AIMD simulations:
(1) AIMD simulations starting from optimised 5-FU@FAU configurations were performed in order to determine internal energies of adsorption ΔUads and to investigate the stability/lability of hydrogen bonds. For this purpose, simulations were run for the 5-FU@FAU configuration, the guest-free zeolite, and for 5-FU in a box, computing three independent trajectories of 12.5 ps (25000 steps) for each system (four trajectories were computed for 5-FAU@SiO2-FAU, see Results). After discarding the first 2.5 ps as equilibration phase, internal energies of adsorption as well as estimated upper and lower boundary values were computed from the production part of the trajectories. The procedure, which follows earlier work,62,63 is described in more detail in section S1.4 of the ESI.† For atom combinations of interest, radial distribution functions (RDFs) were calculated for individual trajectories using VMD,86 and then averaged over all 3 (or 4) trajectories.
(2) To investigate the behaviour of adsorbed 5-FU in the presence of water, AIMD simulations were run for models in which 16 H2O molecules were added in the vicinity of the 5-FU molecule. While the specific amount of H2O molecules included in these calculations is, to a degree, arbitrary, the simulations of 5-FU in liquid water, described below, showed that there are typically about 13 to 14 H2O molecules within a distance of 5 Å from the ring atoms of 5-FU. On this basis, 16 H2O molecules appear as a plausible choice that should correspond to an essentially complete first hydration shell of the adsorbed 5-FU molecule. The simulations of co-adsorbed 5-FU and water considered only H-FAU_1H, H-FAU_2H_1_3, and H-FAU_2H_1_4. Since these simulations were aimed primarily at qualitative insights, only a single trajectory was computed for each system, but for an extended time period of 25 ps.
(3) To analyse the behaviour of 5-FU in liquid water, three independent trajectories were computed for a single 5-FU molecule in a box of 256 H2O molecules, with the size of the box adjusted to match the experimental density of liquid water at ambient conditions (0.997 g cm−3). Starting configurations were generated using the Amorphous Cell module of Materials Studio, using PCFF parameters.
Looking at the numerical values, an adsorption energy of −73 kJ mol−1 is predicted for purely siliceous FAU. The presence of a single framework proton in H-FAU_1H results in a near-doubling of the interaction strength, with a ΔĒads,Boltz value of −144 kJ mol−1. Inclusion of a second proton in the 12MR results in a further increase by about 25%/31% for the “1_3” and “1_4” proton arrangements of H-FAU_2H. If a third proton is present in the same 12MR, however, the adsorption energy is somewhat less negative than for H-FAU_2H, an evolution that can be understood when looking at the local environment of the adsorbed 5-FU molecule (discussed in the following section). For the more Al-rich models containing 16 and 32 protons, the computed adsorption energies are within ±5% of the values obtained for the corresponding low-Al models having the same number and arrangement of protons in a single 12MR (e.g., H-FAU_2H_1_3: −179 kJ mol−1, H-FAU_32H_1_3: −183 kJ mol−1). This indicates that the adsorption energy is dominantly influenced by the local environment, and that the overall Si/Al ratio and the more distant environment have only a minuscule effect (as long as framework topology [and, consequently, framework density] remain [essentially] constant).
For those FAU models where AIMD simulations were run, internal energies of adsorption are also included in Fig. 2 (green bars). Numerical values of ΔUads, including the estimated upper and lower boundary values, are given in the ESI† (Table S5). Typical error bars associated with the AIMD results are on the order of ±10 kJ mol−1, with the larger error of about ±20 kJ mol−1 for H-FAU_2H_1_4 being an exception. For 5-FU@SiO2-FAU, where two configurations are very close in energy (Config1 and Config2, see Fig. 3), a total of four AIMD simulations were run, two starting from each configuration. The resulting internal energy of adsorption is about 9% less negative than the ΔĒads,Boltz value computed with DZVP basis sets. Relative differences obtained for the protonic FAU models show significant scatter, falling between 2% and 13%. Given the error bars of the AIMD-derived ΔUads values, it is clear that these numbers should not be overinterpreted. Nevertheless, it is worth noting that the reduction in interaction strength caused by the inclusion of thermal motion in the simulations is considerably smaller than observed in previous, methodologically analogous investigations of carbamazepine and triclosan adsorption, where differences were on the order of 20% to 30% for all-silica FAU, and amounted to 16% for triclosan in H-FAU_1H.62,63
For 5-FU@SiO2-FAU, two configurations that are very close in energy are shown in Fig. 3. In the first configuration, which is energetically slightly favoured according to DZVP results, the 5-FU molecule lies more or less flat above one 6MR and one 4MR forming the wall of the FAU supercage. As no hydrogen bonds are present, it can be concluded that the molecule interacts with the zeolite predominantly through dispersion interactions. In the second configuration, more favourable according to the calculation using TZVP basis sets, the molecule assumes a more tilted orientation with respect to the curvature of the pore wall, forming a relatively long hydrogen bond via the H1 atom to a framework oxygen atom (Ofw). AIMD simulations starting from both configurations (two from Config1, two from Config2) resulted in fairly similar RDFs. In particular, the H1⋯Ofw RDF (Fig. S2†), averaged over all four trajectories, exhibits only a shoulder in the distance range between 2 and 2.5 Å, not a pronounced maximum. This indicates that hydrogen bonds to oxygen atoms of the all-silica framework are relatively weak and transient, as observed previously for triclosan in all-silica FAU and MOR.63
In the lowest-energy configuration obtained for 5-FU@H-FAU_1H, the framework is deprotonated by the adsorbed molecule, with the proton now being bonded to the O4 atom of 5-FU. As is visible in Fig. 3, the protonated 5-FU species is stabilised by two short hydrogen bonds from the H@O4 and H3 atoms to two Ofw atoms neighbouring the framework Al site. Protonation of O2 is observed in two other configurations (Config1 and Config5 of Table S4†), which are, however, about 10 kJ mol−1 less stable. This agrees with the previous finding by Wielińska et al. that protonation of O4 is energetically favoured.15 Configurations in which a neutral 5-FU molecule is hydrogen-bonded to the framework are even less stable. The Ofw⋯Hfw and Hfw⋯O4 RDFs visualised in Fig. S3† confirm that the proton remains closer to O4 than to Ofw for the majority of the simulation time. However, an overlap of the maxima in the distance range of 1.0 to 1.3 Å points to a certain degree of dynamic exchange.
In contrast to the scenario observed for H-FAU_1H, the 5-FU molecule is only hydrogen-bonded to framework protons in the lowest-energy configurations for both models having two protons within one 12MR (Fig. 3). Among the DFT-optimised structures, only a single configuration with a protonated 5-FU molecule was found (for H-FAU_2H_1_3), which is about 25 kJ mol−1 higher in energy than the most stable configuration. In both configurations shown in Fig. 3, the hydrogen bond involving the O2 atom is somewhat shorter than that involving the O4 atom. Repulsive interactions between the F atom and the framework are likely to prevent the formation of a very short hydrogen bond to O4. Moreover, different partial charge calculation schemes consistently deliver a more negative polarisation of the O2 atom compared to O4 (Table S6†). Additional hydrogen bonds occur between the H1 atom of 5-FU and framework oxygen atoms. In H-FAU_2H_1_4, this hydrogen bond is fairly short, involving an O atom that is part of the AlO4 tetrahedron, whereas a longer bond to an O atom of a Si–O–Si linkage is present in H-FAU_2H_1_3. It is plausible to infer that the presence of a shorter hydrogen bond to a more negatively polarised Ofw atom contributes to the higher stability of the former configuration. The RDFs (Fig. S4 and S5†) indicate that the hydrogen bonds highlighted in Fig. 3 are stable on the timescale of the AIMD simulations. Although there appears to be some fluctuation between the hydrogen-bonded state and short-lived framework deprotonation in the vicinity of the O2 atom, the cumulative RDFs indicate that the proton is located within 1.2 Å of O2 for less than 10% of the total simulation time.
The lowest-energy configuration for H-FAU_3H is qualitatively similar to the one obtained for H-FAU_2H_1_3. Although the H1⋯Ofw hydrogen bond is shorter, the adsorption energy is slightly less negative. As the third proton of the 12MR lies at a distance of 3.25 Å from the H6 atom of 5-FU, repulsive electrostatic interactions between them are likely to be responsible for the somewhat reduced stability. For the more Al-rich models (Fig. S1†), the lowest-energy configurations are generally very similar to those observed for the corresponding H-FAU_1H/2H models, with framework deprotonation occurring in the system with one proton per 12MR (H-FAU_16H) and solely hydrogen-bonded 5-FU in models with two protons per 12MR (H-FAU_32H_1_3/_1_4). Among the whole set of optimised configurations, there are two cases where 5-FU is hydrogen-bonded to framework protons pointing into different 12MRs (one for H-FAU_16H, the other for H-FAU_32H_1_3, they are labelled with asterisks in Table S4†). However, these configurations are energetically less favourable than those where all significant interactions occur within a single 12MR. This indicates that the inclusion of different proton arrangements within one ring may be sufficient to cover the most relevant low-energy bonding situations for a relatively simply and rigid molecule like 5-FU.
Tables S8† compares the calculated frequencies in the range of 1200 to 1800 cm−1 to experimental results obtained for 5-FU in an Ar matrix and to DFT-B3LYP calculations, both reported by Ivanov et al.,70 and to experimental results obtained for solid 5-FU by Datt et al. (available only for the four main modes shown in Fig. 4).20 The root mean square error of the DFT frequencies from the present work with respect to experimental results for 5-FU in an Ar matrix amounts to 17 cm−1 (including only the modes falling between 1200 and 1800 cm−1 listed in Table S8†). This agreement appears satisfactory, especially when considering that the DFT frequencies were directly taken from the calculations, without the application of any scaling factor. Scaling of vibrational frequencies computed with electronic structure methods by a method-dependent scaling factor is often done to minimise systematic deviations from experimental frequencies.87
Upon adsorption in SiO2-FAU, the most significant change in the IR spectrum is a slight red-shift of the two ν(CO) modes, correlated with a modest elongation of the CO bonds. In the higher-frequency region (not included in the figure), the ν(N1–H1) mode is red-shifted due to the presence of the H1⋯Ofw hydrogen bond (Fig. 3). Much more pronounced changes are visible in the IR spectrum of 5-FU@H-FAU_1H: the deprotonation of the framework by 5-FU results in the emergence of a new ν(O4–H) stretching mode at 2000 cm−1, which is very prominently visible. The splitting of the ν(CO) modes increases drastically, with the ν(C4O4) mode being red-shifted due to the bond elongation upon protonation of O4, whereas the ν(C2O2) mode is (slightly) blue-shifted. The ν(C5C6) mode is red-shifted, and the O4–H bending mode appears as a fairly intense signal at 1360 cm−1. A closer look at Table S7† shows that several of the bond distances in the 5-FU ring change rather significantly upon protonation. These changes can explain the shift of a mode associated with ring deformations from 1383 cm−1 in free 5-FU to 1519 cm−1. Altogether, the changes in the IR spectra are so pronounced that it should be straightforward to identify protonated 5-FU species in zeolite pores, even if only a fraction of the adsorbed molecules are protonated. For the 5-FU@H-FAU_2H_1_3 case, where 5-FU is solely hydrogen-bonded, the changes are less distinctive, but still clearly visible: here, both ν(CO) modes are red-shifted, with the larger shift occurring for the ν(C2O2) mode that participates in the stronger hydrogen bond. The ν(CC) mode is also red-shifted, but not as strongly, resulting in partial overlap of the signals arising from these three modes. The Ofw–Hfw bending modes appear at 1237 and 1261 cm−1, respectively, thus overlapping with the ν(C5–F) mode.
For the four modes highlighted in Fig. 4, experimental frequencies of 5-FU adsorbed in H-FAU samples having Si/Al ratios of 2.5, 15, and 30 were reported by Datt et al.20 Since the results for all three samples were essentially identical, it is sufficient to discuss those obtained for the most Si-rich sample: taking free 5-FU in Ar matrix as reference, these authors observed a slight red-shift of the ν(C2O2) mode by −5 cm−1 and a more prominent red-shift of the ν(C4O4) mode by −25 cm−1. These changes are much smaller than those predicted by the DFT calculations for hydrogen-bonded 5-FU adsorbed in H-FAU_2H_1_3, as well as being qualitatively different from those computed for the protonated 5-FU case. Experimentally, a slight blue-shift of the ν(CC) mode by +5 cm−1 was observed, whereas the ν(C5–F) mode was completely unaffected. Neither of these results agree with the computational predictions summarised above. These apparent discrepancies can be explained when taking into account that Datt et al. concluded that 5-FU adsorbed in their samples primarily interacted with Lewis acidic EFAL sites, not with framework protons. While a direct corroboration would require further calculations, the present results support the findings by Datt et al. that adsorption of 5-FU at Brønsted acidic framework protons was not the dominant process occurring in their samples.
These considerations already indicate that the competitive adsorption of water should result in a displacement of 5-FU away from the framework protons. AIMD simulations were employed to investigate the co-adsorption of both species in more detail. First, preliminary FF-based simulations were used to add 16 H2O molecules to the lowest-energy configurations computed for 5-FU in H-FAU_1H, H-FAU_2H_1_3, and H-FAU_2H_1_4. AIMD simulations with a total duration of 25 ps were then run for each of these (5-FU + 16H2O)@FAU systems. In the following, key observations for each trajectory are presented, focussing on representative snapshots from different stages of the simulation and selected RDFs (including time-resolved RDFs for different time intervals where relevant).
Fig. 5 shows the environment of the adsorbed 5-FU molecule in (5-FU + 16H2O)@H-FAU_1H after 0.5 ps and 2.5 ps and at the end of the AIMD simulation (after 25 ps). At the beginning of the simulation, the protonated 5-FU molecule is hydrogen-bonded to the framework, as in the lowest-energy configuration of 5-FU@H-FAU_1H in the absence of water (Fig. 3). After approximately 1.5 ps, the proton moves from the 5-FU molecule to one of the surrounding water molecules, and the RDFs computed over different time intervals, collected in Fig. S11,† show that no new O4–H bonds form at later stages of the simulation. The 5-FU molecule progressively moves away from its initial location near the framework Al atom to another part of the cage, whereas the vicinity of the Al atom is occupied by a positively charged cluster formed by water molecules and one hydronium ion. The 5-FU molecule that was initially bonded to the framework now forms numerous hydrogen bonds to surrounding water molecules. Additional RDFs compiled in Fig. S10† indicate that the O2 atom of 5-FU remains unprotonated, and that there is no formation of any new Ofw–Hfw bond; in other words, the “excess” proton remains part of the water cluster throughout the simulation.
In the case of (5-FU + 16H2O)@H-FAU_2H_1_3, the following chain of events can be identified based on representative snapshots shown in Fig. 6 (after 0.5, 2.5, 25 ps): at the very beginning of the simulation, the framework proton that participated in a hydrogen bond to the O4 atom of 5-FU in the water-free structure (Fig. 3) is removed from the framework, protonating the O4 atom. While the other framework proton is still bonded to the framework in the first snapshot shown in Fig. 6, it is removed shortly afterwards, forming a covalent O2–Hfw bond. At approximately the same time, the O4 atom is deprotonated again, with the proton being transferred to a water molecule. This is the situation after 2.5 ps. After about 6 ps, the O2 atom is also deprotonated. At the end of the simulation, the vicinity of the Al atoms is occupied by the cluster consisting of H2O molecules and two hydronium ions, and the 5-FU molecule is displaced considerably from its initial location. The time-resolved RDFs shown in Fig. S13† further illustrate these findings.
Three representative snapshots of (5-FU + 16H2O)@H-FAU_2H_1_4 are shown in Fig. 7, again visualising the system after 0.5, 2.5, and 25 ps. In this case, one of the framework protons moves to the O2 atom of 5-FU in the very initial stages of the simulation, while the O4 atom remains hydrogen-bonded to the other framework proton after 0.5 ps. After 2.5 ps, however, a rather unexpected situation is observed: at this point, both the O2 and O4 atom of 5-FU are protonated, and the N3 atom is deprotonated. Carrying a net charge of +1, 5-FU is now in a different tautomeric form than observed in all previous calculations involving protonated 5-FU. Although this tautomer was found to be less stable than the form having a single protonated O atom in a previous computational study,15 it does not seem implausible to surmise that it can be stabilised, at least temporarily, through localised interactions in a specific environment. At the end of the 25 ps trajectory, the proton that was previously attached to the O4 atom has been transferred to a water molecule, whereas the O2 and N3 atoms remain protonated and deprotonated, respectively. The absence of additional peaks in the RDFs in the time interval from 20 to 25 ps confirms the stability of this scenario on the picosecond timescale. Although the 5-FU molecule still forms a hydrogen bond to a Ofw atom neighbouring one of the two Al sites, Fig. 7 clearly shows the movement away from its initial location at the centre of the 12MR.
Fig. 7 Snapshots from three different stages of the AIMD simulation of (5-FU + 16H2O)@H-FAU_2H_1_4. Light blue arrows highlight deprotonation of the ring N atom. |
Although there are non-negligible differences among the individual AIMD trajectories, some common features can be identified: in all three cases, the co-adsorption of water and 5-FU results in framework deprotonation, a finding that falls in line with previous observations made for water clusters89–91 and for co-adsorbed water and triclosan in protonic zeolites.63 Although protonated 5-FU occurs as an intermediate species during the simulations, the protons are transferred to water molecules over the course of the 25 ps simulations. In the final snapshots, 5-FU is always in the neutral form, interacting through hydrogen bonds with the positively charged water cluster and the deprotonated zeolite framework.
It can be inferred from these observations that the behaviour of co-adsorbed 5-FU and water in the zeolite pores is strongly impacted by the Brønsted acidity of the framework protons. It is therefore interesting to compare this behaviour to that of 5-FU in aqueous solution. For this purpose, additional AIMD simulations were run for a box containing one 5-FU and 256 H2O molecules, considering three different starting configurations. The RDFs computed from these trajectories, shown in Fig. S16,† provide no indications for protonation of 5-FU oxygen atoms or deprotonation of nitrogen atoms, in other words, the neutral tautomer shown in Fig. 1 is stable in aqueous solution. This corroborates that the proton(s) that the zeolite provides are responsible for the proton transfers between 5-FU and its environment observed in the AIMD simulations.
Both in terms of adsorption energies and equilibrium structures, fairly similar results were obtained for zeolite models having vastly different Si/Al ratios, but identical proton arrangements in the 12MR where 5-FU is adsorbed. This indicates that models containing only a few framework protons can be useful to represent specific local environments. However, it should be kept in mind that the presence of 2 or 3 protons in the same 12MR at Si/Al ratios of 95 and 63, respectively, appears extremely improbable when assuming a more or less statistic distribution. Thus, these structures should be seen as representatives of possible local arrangements at lower Si/Al ratios, rather than being realistic models. An additional simplification made in the present study was the assumption that the framework protons are always bonded to O1 oxygen atoms, neglecting the possibility to occupy different sites at finite temperature.
The computational prediction of IR spectra clearly showed that different bonding environments of 5-FU could be easily distinguished by means of IR spectroscopy. Although significant discrepancies between the DFT results and previously reported experimental data were apparent, these are most probably attributable to the interaction of 5-FU with EFAL sites in the experimentally studied samples. Future computational work will look into a DFT-based modelling of the interaction between 5-FU and EFAL sites, taking EFAL models proposed in the recent literature as starting points.95,96 While it might appear counterintuitive that all three samples studied experimentally in prior work show identical shifts of IR-active modes, it has to be considered that even the most Si-rich sample contained, on average, 6 Al atoms per unit cell (Si/Al ratio of 30). Thus, it still seems possible that a significant fraction of the adsorbed 5-FU molecules can interact with EFAL sites. However, these aspects cannot be elucidated with calculations alone, and a more detailed understanding of the role of different aluminium environments would require further experimental investigations of 5-FU adsorption in zeolite samples having well-quantified amounts of framework and extra-framework Al atoms.
Finally, AIMD simulations including co-adsorbed water molecules showed that the water molecules have a higher affinity towards the framework protons, displacing 5-FU from their vicinity. Thus, it seems unlikely that interactions with the framework protons would contribute significantly to the stabilisation of 5-FU in zeolite pores in the presence of water. In the context of contaminant removal, it is therefore pivotal to maximise adsorbent hydrophobicity, as already discussed in earlier work, for example, for the case of triclosan adsorption.63,97 For drug delivery, on the other hand, there are some additional implications: if 5-FU was loaded in the absence of water (e.g., from non-aqueous solvents like acetone, as reported in some previous studies21,24,26) and the 5-FU@zeolite composite stored without exposure to humidity, the drug could be released upon contact with water, for example, in the human body. Potentially, the release kinetics could be tailored through a variation of the Si/Al ratio and/or the framework topology. In the prior experimental investigation by Datt et al., the sample having the lowest Si/Al ratio showed essentially no 5-FU release in simulated body fluid, pointing to a strong binding of 5-FU to EFAL sites.20 Based on the present work, it is suggested to focus on protonic zeolites containing few EFAL sites in future experimental investigations of zeolites as 5-FU carrier materials. An additional possibility to tailor the affinity towards 5-FU (or other drug molecules of interest) is the use of cation-exchanged zeolites containing different metal cations. A systematic exploration using an analogous computational protocol is planned for future work.
Footnotes |
† Electronic supplementary information (ESI) available: The supplementary PDF contains details of force field simulations (S1.1 to S1.3) and of calculations of internal energies of adsorption (S1.4), tables of adsorption energies (S2.1) and internal energies of adsorption (S2.2), a figure showing additional low-energy configurations (S2.3), figures showing input structures used in AIMD simulations of (5 FU + 16H2O)@FAU (S2.7), RDF plots (S2.4, S2.8, S2.9) as well as tables reporting additional results (partial charges, S2.5; intramolecular bond distances and vibrational modes, S2.6). See DOI: https://doi.org/10.1039/d4ce00344f |
‡ Input files and datasets are available from https://doi.org/10.6084/m9.figshare.25541446. The deposited ZIP archives include sample CP2K input files (Sample_inputs.zip), low-energy configurations and output files from TZVP calculations (Configs_TZVP_results.zip), results of vibrational calculations (Vibrational_calculations.zip), and results of AIMD simulations (5-FU_FAU_AIMD_results.zip, 5-FU_H2O_FAU_AIMD_results.zip, 5-FU_H2O_AIMD_results.zip). |
This journal is © The Royal Society of Chemistry 2024 |