Next Article in Journal
CRISPR-Cas9 Mediated Mutation in OsPUB43 Improves Grain Length and Weight in Rice by Promoting Cell Proliferation in Spikelet Hull
Next Article in Special Issue
Melissa officinalis: Composition, Pharmacological Effects and Derived Release Systems—A Review
Previous Article in Journal
Diffracted X-ray Tracking Method for Measuring Intramolecular Dynamics of Membrane Proteins
Previous Article in Special Issue
An Engineered Nanocomplex with Photodynamic and Photothermal Synergistic Properties for Cancer Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ab Initio Insight into the Interaction of Metal-Decorated Fluorinated Carbon Fullerenes with Anti-COVID Drugs

by
Konstantin P. Katin
1,2,*,
Alexey I. Kochaev
1,3,
Savas Kaya
4,
Fadoua El-Hajjaji
5 and
Mikhail M. Maslov
1,2
1
Laboratory of Computational Design of Nanostructures, Nanodevices, and Nanotechnologies, Research Institute for the Development of Scientific and Educational Potential of Youth, Aviatorov Str. 14/55, 119620 Moscow, Russia
2
Institute of Nanotechnologies in Electronics, Spintronics and Photonics, National Research Nuclear University “MEPhI”, Kashirskoe Shosse 31, 115409 Moscow, Russia
3
Research and Education Center “Silicon and Carbon Nanotechnologies”, Ulyanovsk State University, 42 Leo Tolstoy Str., 432017 Ulyanovsk, Russia
4
Department of Chemistry, Faculty of Science, Cumhuriyet University, Sivas 58140, Turkey
5
Engineering Laboratory of Organometallic, Molecular Materials, and Environment, Faculty of Sciences, University Sidi Mohamed Ben Abdellah, Fez 1796, Morocco
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(4), 2345; https://doi.org/10.3390/ijms23042345
Submission received: 20 January 2022 / Revised: 14 February 2022 / Accepted: 17 February 2022 / Published: 21 February 2022
(This article belongs to the Special Issue Nanotechnology in Targeted Drug Delivery)

Abstract

:
We theoretically investigated the adsorption of two common anti-COVID drugs, favipiravir and chloroquine, on fluorinated C60 fullerene, decorated with metal ions Cr3+, Fe2+, Fe3+, Ni2+. We focused on the effect of fluoridation on the interaction of fullerene with metal ions and drugs in an aqueous solution. We considered three model systems, C60, C60F2 and C60F48, and represented pristine, low-fluorinated and high-fluorinated fullerenes, respectively. Adsorption energies, deformation of fullerene and drug molecules, frontier molecular orbitals and vibrational spectra were investigated in detail. We found that different drugs and different ions interacted differently with fluorinated fullerenes. Cr3+ and Fe2+ ions lead to the defluorination of low-fluorinated fullerenes. Favipiravir also leads to their defluorination with the formation of HF molecules. Therefore, fluorinated fullerenes are not suitable for the delivery of favipiravir and similar drugs molecules. In contrast, we found that fluorine enhances the adsorption of Ni2+ and Fe3+ ions on fullerene and their activity to chloroquine. Ni2+-decorated fluorinated fullerenes were found to be stable and suitable carriers for the loading of chloroquine. Clear shifts of infrared, ultraviolet and visible spectra can provide control over the loading of chloroquine on Ni2+-doped fluorinated fullerenes.

1. Introduction

Carbon-based nanomaterials interact with a variety of pharmaceutical drugs [1]. Fullerenes are the smallest carbon nanostructures. They are hollow spherical cages, and sp2-hybridized carbon atoms form pentagons and hexagons on their surface [2,3]. Fullerenes are quite stable and soluble in organic solvents. The functionalization of fullerenes by various radicals provides a wide variation in their solubility. Fullerenes consist of an even number of atoms. The smallest fullerene, C20, has a strained skeleton consisting of 12 pentagons [4]. However, it has a high kinetic stability [5]. Giant fullerenes can contain several hundred of atoms [6,7]. However, fullerenes C60 and C70 (especially C60) are the most stable and are obtained with the highest yield [7].
Carbon fullerenes are considered promising nanoparticles for drug delivery [8]. The most common fullerene, C60, has a well-reproducible stable carbon skeleton. In contrast to other sp2-hybridized carbon nanostructures, such as nanotubes and graphene flakes, the shape and size of C60 can be easily reproduced from the experimental point of view. The chirality of nanotubes and the exact size of nanoflakes are difficult to control, whereas C60 fullerene has a well-defined structure. In addition, the curvature of the fullerene surface provides its increased chemical activity [9,10,11]. Like other carbon nanostructures, fullerenes are not toxic to the human body [12,13,14]. The doping of fullerenes results in a wide family of derivative compounds, which possess a wide variety of physical and chemical characteristics, including adsorption properties, electronic structure and solubility [15,16,17].
There are three common ways of fullerenes doping. Endohedral doping involves the formation of guest-host complexes, in which the doping atom or molecule cannot leave the fullerene cage due to spatial difficulties. Substitution doping involves the embedding of foreign atoms into the fullerene skeleton. This type of doping is very common for drug delivery applications since the embedded atom becomes the reaction centre to which the drugs can attach. In particular, fullerenes doped with various elements, including silicon [18], germanium, boron and nitrogen [19], metals [20,21] and other elements [22,23,24], were considered for drug delivery. External doping implies the formation of the chemical bond between the carbon atom and the doping functional group. A significant advantage of such kind of doping is that it does not require the destruction of the carbon skeleton of fullerene. Despite the strained structure of fullerene, the breaking of carbon–carbon bonds in its skeleton requires high activation energy of about 7 eV [25] and, therefore, high temperatures or very active reagents. Thus, external doping seems to be the “gentlest” type of doping, providing good reproducibility of the structure of doped fullerenes.
Fluorine is a very suitable dopant for carbon nanoparticles used as the drug carriers. This very electronegative element forms strong hydrogen bonds with drug molecules. Recently, synthesized and approbated fluorinated graphene nanoflakes have demonstrated other advantages [26,27,28]. The carbon–fluorine bonds provide greater activity in the near-infrared range, providing heat delivery to the system (the heat is necessary for the release of drugs and additional therapeutic effects [29]). In addition, fluorine helps to track nanoparticles inside the human body using nuclear magnetic resonance [30]. For these reasons, fluorinated fullerenes appear to be the preferred carriers of drugs compared to pristine ones.
Fluorinated fullerenes do not have free chemical bonds; therefore, they can interact with drugs molecules only through hydrogen or van-der-Waals bonds. Stronger bonds are possible via the introduction of metal ions. These ions act as linkers, forming complexes with drug molecules and fullerenes. In Ref. [31], the interaction of metal-decorated fullerenes with several drugs that presumably exhibit anti-COVID19 activity was investigated. The authors considered five promising transition metals—Ti, Cr, Fe, Ni, and Zn. Among them, three metals (Cr, Fe, Ni) demonstrated the best interaction with drug molecules [31]. Note that the presence of negatively charged fluorine atoms in the system can contribute to the attraction between metal ions and fullerenes.
In this paper, we examined the effect of fluorination on the interaction of fullerenes with metal ions and two anti-COVID19 drugs, favipiravir and chloroquine. Their structures are presented in Figure 1. According to recent studies, “favipiravir has shown rapid viral clearance and faster clinical improvement” [32]. It was approved in many countries, including Russia and India, for COVID-19 treatment. Recent studies also suggest the high therapeutic efficiency of chloroquine [33,34]. Note that the adsorption of both drugs on pristine and doped fullerenes was intensively investigated [19,21,35,36], but the effect of fullerene fluorination was not considered so far. Here we applied the electron density functional theory to clarify the structural and spectroscopic features of fluorinated fullerenes decorated with previously selected metal ions (Cr, Fe, Ni) and investigate their interaction with two mentioned drugs. The rest of the article is organized as follows. Section 2 describes the results of the study. Section 3 is devoted to the description of computational details and methods. Section 4 contains the conclusion.

2. Results and Discussion

2.1. Atomistic Models of Fluorinated Fullerenes

The concentration of fluorine strongly affects the interaction of fluorinated fullerene with ions and drugs. At a low fluorine concentration, the fullerene surface remains partially free and can participate in the interaction with drugs. On the other hand, all fullerene surfaces are coated with fluorine radicals at a higher fluorine concentration. To cover both cases, we considered two fluorinated isomers, C60F2 and C60F48. They correspond to the lowest and the highest fluorinated fullerenes, respectively, experimentally observed in significant yields [37]. The structures of these models are described in computational details section.

2.2. Interaction of Drugs Molecules with the Pristine and Fluorinated Fullerenes

According to previous studies, the drugs interact poorly with the pristine C60. The binding energy between C60 and favipiravir in the gas phase is as low as 0.3 eV, and the presence of solvents changes it only by about 10% [35]. Chloroquine does not contain active oxygen atoms and, therefore, interacts even weaker with C60, possessing the binding energy of about 0.06 eV [21]. Therefore, any modification of pristine C60 is necessary for efficient drugs loading.
First, we investigated the interaction of fluorinated fullerenes C60F2 and C60F48 with favipiravir. During geometry optimizing, hydroxyl groups of favipiravir become near fluorine atoms to form O-H..F hydrogen bonds. However, we observed the detachment of fluorine from fullerene and the forming of the HF acid after that. We tried many different initial positions of fullerene and the drug, but we always observed the formation of HF as a result of optimization. Similar defluorination of fullerene by OH groups was observed experimentally [38]. Thus, favipiravir results in defluorination of fullerene. Therefore, fluorinated fullerenes are not a suitable drug carrier for favipiravir and similar drugs with active OH groups. In contrast, we observed the successful loading of chloroquine on both C60F2 and C60F48 (see Figure 2).The calculated results are collected in Table 1. One can see that the presence of fluorine leads to an increase in the binding energy by more than ten times. For low-fluorinated fullerene C60F2, electrostatic attraction between nitrogen and fluorine combines with π-π interaction of non-coated fullerene surface with aromatic rings of chloroquine (the distance between parallel aromatic rings belong to fullerene and drug is about 3.3 Å). As a result, low-fluorinated fullerene C60F2 possesses almost as strong binding to the drug as high-fluorinated fullerene C60F48. Note that slight interaction of fluorine with carbon aromatic rings was reported in Ref. [39]. The stronger distortion of the chloroquine in the C60F48 + drug complex can be explained by the attraction of the non-aromatic part of the drug to fluorine. However, distortions of both fullerenes and drugs are quite slight, and their total contribution to interaction energy does not exceed 10%. Electron localization function for both complexes presented in Figure 2 are plotted in Supplementary materials (Figure S1a,b).

2.3. Interaction of Metals Ions Molecules with the Pristine and Fluorinated Fullerenes

Binding energies between chloroquine drug and fluorinated fullerenes considered above are lower than 0.5 eV. Such energies are insufficient to provide durable drug loading. Therefore, the introduction of metal ions is required. Cr, Fe, and Ni were recognized as the best metals for doping fullerene-based drug carriers [31]. Therefore, we considered the interaction of Cr3+, Fe2+, Fe2+, and Ni2+ ions with pristine and fluorinated C60. We tried different configurations of doping ions (top, bridge, hollow) and selected ones with the lowest energies. Resulted geometries are shown in Figure 3, whereas corresponding characteristics are presented in Table 2.
We found that all ions strongly interacted with the carbon cage. The ion formed three bonds with carbon atoms, whereas others ions preferred a bridge configuration under the centre of the C-C bond. Note that the observed behaviour of charged ions in solution differs from the behaviour of neutral atoms in a vacuum (for example, neutral Ni atoms prefer a hollow position on fullerene [40]). One can see that ions with a charge of +3 possess much higher Eb values than ions with the charge of +2. It is a remarkable fact that the low concentration of fluorine results in the increase of the Eb values for all metal ions. Preferable positions of metal ions on low-fluorinated fullerenes are located near the fluorine atoms. However, Cr3+ and Fe2+ ions lead to the detachment of fluorine from the carbon cage (see Figure 3). In highly fluorinated fullerene, the carbon surface is completely coated. Therefore, metal ions bind only to fluorine. The nickel ion forms two bonds with fluorine atoms, whereas all other ions form three bonds each. The binding energy, in this case, is significantly lower than for low-fluorinated fullerenes (see Table 2).
Deformation energies of fullerenes C60Fn (n = 0, 2, 48) were calculated as energy differences between distorted and relaxed fluorinated fullerenes, as was described above. We were not surprised that the deformation energies are significantly higher for fluorinated fullerenes, in which fluorine atoms can rotate to adjust to doping ions.
It should be noted that the nickel ion Ni2+ seems to be the most suitable for doping fullerene. This ion does not threaten the stability of fluorinated fullerene and demonstrates moderate binding energy, which is sufficient for reliable adsorption on fullerene.

2.4. Loading of Chloroquine on Metal-Decorated Fluorinated Fullerenes

We examined the interaction of the chloroquine drug with the pristine and fluorinated fullerenes doped with metal ions. We observed the formation of the covalent M–N bond between the metal ion M and the drug. Corresponding bonds lengths and other calculated characteristics of fluorinated fullerenes doped with metal ions and loaded with chloroquine are presented in Table 3. We were not able to consider complexes containing C60F2Cr3+ and C60F2Fe2+ compounds, because Cr3+ and Fe2+ ions induce fullerene defluorination, as described in above.
Note that a high concentration of fluorine significantly enhances the binding between the metal-doped fullerenes and the drug. Complexes containing Ni2+ ion demonstrated a moderate binding with the drug, suitable for reliable loading. The geometries of these complexes are shown in Figure 4. This figure confirms that both covalent and non-covalent interactions occur between the drug and fullerene.
Deformation energies of the drug and curriers are also presented in Table 3. One can see that they make up a significant proportion of the binding energy. Deformation energies of high-fluorinated fullerenes are substantially higher. This fact indicates that high-fluorinated fullerenes can “adjust” their shapes to the drug. In most cases, fullerenes doped with Ni2+ ion provide the lowest deformation of both drug and carrier, as indicated in Table 3.

2.5. Spectral Fingerprints of the C60F2Ni2+-Chloroquine and C60F48Ni2+-Chloroquine Complexes

Here we considered spectral features of Ni-doped fluorinated fullerenes loaded with chloroquine. We chose nickel as the best doping element, providing moderate binding energies and low drug deformations and not leading to carrier defluorination. The formation of a new M-N bond and the deformation of molecules associated with the drug loading should affect the optical spectra of the considered systems. Thus, one can use spectroscopic methods to control the process of loading. Figure 5 shows the ultraviolet and visible spectra of the drug, carriers and corresponding complexes. It can be seen that the activity of the drug in the high-energy region (200–300 nm) disappears when it loads on fullerene. In addition, the formation of the “C60F48Ni2+ + chloroquine” complex is accompanied by the appearance of a peak at 418 nm.
Figure 6 shows the active infrared frequencies and the corresponding intensities for the systems under consideration. It can be seen from Figure 6 that the spectrum of the “carrier + drug” system is not a simple sum of the spectra of its parts. Therefore, one can control the loading process by spectral measurements. Integral infrared absorption intensity of the “C60F2Ni2+ + drug” complex is much higher than absorption intensities of its parts (see Figure 6a). Therefore, drug loading on C60F2Ni2+ can be easily indicated. Loading the drug onto the C60F48Ni2+ cluster does not lead to such a significant change in the infrared spectrum. However, drug loading on C60F48Ni2+ results to an appearance of a series of active frequencies in the range of 590–700 cm−1.

3. Computational Details

For C60F2, we constructed the most stable isomer, in which two fluorine atoms occupy the same carbon hexagon in the para position [41]. For C60F48, we selected the low-energy isomer with D3 symmetry, which provides the best fit of the X-ray diffraction spectra of C60F48 [42]. For geometry optimizing, we used B3LYP exchange–correlation functional, which is widely used for investigations of drugs on fullerenes (see, for example, [43,44]). Two Pople’s electronic basis sets 6–31 G*_ldz [45] and 6–31 G** [46], were consistently used for metal ions and other atoms, respectively. Grimme’s D3 corrections [47] were introduced to take into account non-covalent dispersion interactions. The solvent effect was simulated with the conductor-like polarization model COSMO [48]. We assumed that the solvent is water with a dielectric constant of 78.4. Note that all calculations, including geometry optimizing, were carried out by taking the solvent into account. To accelerate the calculations rate, we combined the power of GPU-based TeraChem software [49] with the efficiency of the geomeTRIC optimizer [50]. Vibrational frequencies and infrared (IR) spectra were calculated with the same B3LYP functional. Ultraviolet and visible spectra (UV–Vis) were calculated with CAM-B3LYP functional, which is more suitable for such calculation due to taking into account the long-range Coulomb correlations [51]. Twenty excited states were considered with the Tamm-Dancoff time-dependent density functional approach [52].
The binging energies were calculated as follows. To characterize chloroquine–fullerene interaction quantitatively (Section 2.2, Table 1), we calculated binding energies Eb as:
Eb = E(C60Fn) + E(chloroquine) − E(C60Fn + chloroquine) + BSSE.
Here n = 2 or 48 for low- and high-fluorinated fullerene, respectively. We also calculated deformation energies Edef as the energy differences between distorted fullerene/drugs included into the complex and relaxed fullerene/drugs. The total interaction energy between distorted fullerene and drugs can be calculated as a sum of three components:
Eint = Eb + Edef(C60Fn) + Edef(chloroquine).
To characterize the interaction of pristine or fluorinated fullerene with a metal ion (Section 2.3, Table 2), binding energies Eb were calculated as
Eb = E(C60Fn) + E(M ion) − E(C60FnM ion) + BSSE.
Here n = 0, 2 or 48 for pristine, low- and high-fluorinated fullerene, respectively. For metal-decorated fluorinated fullerenes loaded with drugs (Section 2.4, Table 3), the binding energies Eb were calculated as
Eb = E(C60FnM ion) + E(drug) − E(complex) + BSSE.
Here n = 0, 2 or 48 for pristine, low- and high-fluorinated fullerene, respectively; M = Cr3+, Fe2+, Fe3+, Ni2+. In Formulas (1), (3), and (4), basis set superposition error (BSSE) were taken into account with the ghost atoms method implemented in TeraChem.

4. Conclusions

Carbon fullerene derivatives are interesting because they combine the advantages of the pristine fullerene (controlled size, high stability and biocompatibility) with additional improvements related to dopants. Systems with two different types of dopants are particularly complex since it is necessary to consider the interaction of dopants with each other. Here we have considered fullerenes functionalized simultaneously by fluorine and metal ions. It is not surprising that not all of the considered systems turned out to be stable: Cr3+ and Fe2+ ions led to the rupture of C–F bonds. In contrast, Fe3+ and Ni2+ ions are compatible with fluorinated fullerenes. However, Fe3+ can be reduced to Fe2+ under the action of antioxidants contained in the human body [53]. Therefore, we selected the only suitable ion, Ni2+, which does not violate the structures and stability of fluorinated fullerenes. Therefore, Ni-doped fluorinated fullerenes combine high adsorption to drugs due to the presence of the Ni2+ ion with well-known advantages of fluorinated carbon nanosystems. We expect that such complexes can be considered a basis for carriers suitable for chloroquine and other drugs with similar molecular structures.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23042345/s1. Reference [54] is cited in the supplementary materials.

Author Contributions

Conceptualization, K.P.K. and M.M.M.; methodology, K.P.K., A.I.K. and M.M.M.; software, K.P.K., A.I.K., S.K. and M.M.M.; validation, K.P.K., A.I.K., S.K. and M.M.M.; formal analysis, K.P.K., A.I.K., S.K. and M.M.M.; investigation, K.P.K., F.E.-H. and M.M.M.; resources, K.P.K., A.I.K. and M.M.M.; data curation, K.P.K., F.E.-H. and M.M.M.; writing—original draft preparation, K.P.K. and M.M.M.; writing—review and editing, A.I.K., S.K. and F.E.-H.; visualization, K.P.K., A.I.K. and M.M.M.; supervision, K.P.K. and M.M.M.; project administration, K.P.K.; funding acquisition, K.P.K. All authors have read and agreed to the published version of the manuscript.

Funding

The presented study was performed with the financial support of the Russian Science Foundation (Grant No. 20-73-00245).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data is available as Supplementary Materials to this article.

Acknowledgments

Konstantin P. Katin is grateful to DSEPY-RI for the provided computing resources and software as well as comprehensive support of the presented study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Veclani, D.; Tolazzi, M.; Melchior, A. Molecular Interpretation of Pharmaceuticals’ Adsorption on Carbon Nanomaterials: Theory Meets Experiments. Processes 2020, 8, 642. [Google Scholar] [CrossRef]
  2. Kroto, H.W.; Heath, J.R.; O’Brien, S.C.; Curl, R.F.; Smalley, R.E. C60: Buckminsterfullerene. Nature 1985, 318, 162. [Google Scholar] [CrossRef]
  3. Goodarzi, S.; Da Ros, T.; Conde, J.; Sefat, F.; Mozafari, M. Fullerene: Biomedical Engineers Get to Revisit an Old Friend. Materials Today 2017, 20, 460. [Google Scholar] [CrossRef] [Green Version]
  4. Wang, Z.; Ke, X.; Zhu, Z.; Zhu, F.; Ruan, M.; Chen, H.; Huang, R.; Zheng, L. A New Carbon Solid Made of the World’s Smallest Caged Fullerene C20. Phys. Lett. A 2001, 280, 351–356. [Google Scholar] [CrossRef]
  5. Katin, K.P.; Podlivaev, A.I. Dynamic Characteristics of the Low-Temperature Decomposition of the C20 Fullerene. Phys. Solid State 2010, 52, 436. [Google Scholar] [CrossRef] [Green Version]
  6. Rysaeva, L.K.; Lobzenko, I.P.; Baimova, J.A.; Dmitriev, S.V.; Zhou, K. Modeling C540-C20 Fullerene Collisions. Rev. Adv. Mater. Sci. 2018, 57, 143. [Google Scholar] [CrossRef]
  7. Dresselhaus, M.S.; Dresselhaus, G.; Eklund, P.C. Science of Fullerenes and Carbon Nanotubes; Elsevier: Amsterdam, The Netherlands, 1996. [Google Scholar] [CrossRef] [Green Version]
  8. Kumar, M.; Raza, K. C60-Fullerenes as Drug Delivery Carriers for Anticancer Agents: Promises and Hurdles. Pharm. Nanotech. 2018, 5, 169. [Google Scholar] [CrossRef]
  9. Staykov, A.; Ooishi, Y.; Ishihara, T. Immobilizing Metal Nanoparticles on Single Wall Nanotubes. Effect of Surface Curvature. J. Phys. Chem. C 2014, 118, 8907. [Google Scholar] [CrossRef]
  10. Prudkovskiy, V.S.; Katin, K.P.; Maslov, M.M.; Puech, P.; Yakimova, R.; Deligeorgis, G. Efficient Cleaning of Graphene from Residual Lithographic Polymers by Ozone Treatment. Carbon 2016, 109, 221. [Google Scholar] [CrossRef]
  11. Mesarič, T.; Baweja, L.; Drašler, B.; Drobne, D.; Makovec, D.; Dušak, P.; Dhawan, A.; Sepčić, K. Effects of Surface Curvature and Surface Characteristics of Carbon-Based Nanomaterials on the Adsorption and Activity of Acetylcholinesterase. Carbon 2013, 62, 222. [Google Scholar] [CrossRef]
  12. Andrievsky, G.; Klochkov, V.; Derevyanchenko, L. Is the C60 Fullerene Molecule Toxic? Fuller. Nanotub. Carbon Nanostruct. 2005, 13, 363. [Google Scholar] [CrossRef]
  13. Kolosnjaj, J.; Szwarc, H.; Moussa, F. Toxicity Studies of Fullerenes and Derivatives. In Bio-Applications of Nanoparticles; Springer: Berlin/Heidelberg, Germany, 2007; pp. 168–180. [Google Scholar] [CrossRef]
  14. Spohn, P.; Hirsch, C.; Hasler, F.; Bruinink, A.; Krug, H.F.; Wick, P. C60 Fullerene: A Powerful Antioxidant or a Damaging Agent? The Importance of an in-Depth Material Characterization Prior to Toxicity Assays. Environ. Pollut. 2009, 157, 1134. [Google Scholar] [CrossRef] [PubMed]
  15. Yáñez-Sedeño, P.; Campuzano, S.; Pingarrón, J. Fullerenes in Electrochemical Catalytic and Affinity Biosensing: A Review. C 2017, 3, 21. [Google Scholar] [CrossRef] [Green Version]
  16. Sharoyko, V.V.; Ageev, S.V.; Meshcheriakov, A.A.; Akentiev, A.V.; Noskov, B.A.; Rakipov, I.T.; Charykov, N.A.; Kulenova, N.A.; Shaimardanova, B.K.; Podolsky, N.E.; et al. Physicochemical Study of Water-Soluble C60(OH)24 Fullerenol. J. Mol. Liq. 2020, 311, 113360. [Google Scholar] [CrossRef]
  17. Akentiev, A.V.; Gorniaia, S.B.; Isakov, N.A.; Lebedev, V.T.; Milyaeva, O.Y.; Sedov, V.P.; Semenov, K.N.; Timoshen, K.A.; Noskov, B.A. Surface Properties of Fullerenol C60(OH)20 Solutions. J. Mol. Liq. 2020, 306, 112904. [Google Scholar] [CrossRef]
  18. Gökpek, Y.; Bilge, M.; Bilge, D.; Alver, Ö.; Parlak, C. Adsorption Mechanism, Structural and Electronic Properties: 4-Phenylpyridine & Undoped or Doped (B or Si) C60. J. Mol. Liq. 2017, 238, 225. [Google Scholar] [CrossRef]
  19. Muz, İ.; Göktaş, F.; Kurban, M. A Density Functional Theory Study on Favipiravir Drug Interaction with BN-Doped C60 Heterofullerene. Phys. E 2022, 135, 114950. [Google Scholar] [CrossRef]
  20. Rad, A.S.; Shahavi, M.H.; Esfahani, M.R.; Darvishinia, N.; Ahmadizadeh, S. Are Nickel- and Titanium- Doped Fullerenes Suitable Adsorbents for Dopamine in an Aqueous Solution? Detailed DFT and AIM Studies. J. Mol. Liq. 2021, 322, 114942. [Google Scholar] [CrossRef]
  21. Bagheri Novir, S.; Aram, M.R. Quantum Mechanical Simulation of Chloroquine Drug Interaction with C60 Fullerene for Treatment of COVID-19. Chem. Phys. Lett. 2020, 757, 137869. [Google Scholar] [CrossRef]
  22. Parlak, C.; Alver, Ö. A Density Functional Theory Investigation on Amantadine Drug Interaction with Pristine and B, Al, Si, Ga, Ge Doped C60 Fullerenes. Chem. Phys. Lett. 2017, 678, 85. [Google Scholar] [CrossRef]
  23. Muz, İ.; Kurban, M. A First-Principles Evaluation on the Interaction of 1,3,4-Oxadiazole with Pristine and B-, Al-, Ga-Doped C60 Fullerenes. J. Mol. Liq. 2021, 335, 116181. [Google Scholar] [CrossRef]
  24. Baei, M.T.; Soltani, A.; Rajabzadeh, H.; Tazikeh-Lemeski, E. Structural and Electronic Properties of XY-Doped (AlN, AlP, GaN, GaP) C58 Fullerenes: A DFT Study. Russ. J. Inorg. Chem. 2017, 62, 1067–1076. [Google Scholar] [CrossRef]
  25. Podlivaev, A.I.; Katin, K.P.; Lobanov, D.A.; Openov, L.A. Specific Features of the Stone-Wales Transformation in the C20 and C36 Fullerenes. Phys. Solid State 2011, 53, 215. [Google Scholar] [CrossRef]
  26. Khosravanian, A.; Moslehipour, A.; Ashrafian, H. A Review on Bioimaging, Biosensing, and Drug Delivery Systems Based on Graphene Quantum Dots. Prog. Chem. Biochem. Res. 2021, 4, 44. [Google Scholar] [CrossRef]
  27. Antonova, I.; Nebogatikova, N.; Zerrouki, N.; Kurkina, I.; Ivanov, A. Flexibility of Fluorinated Graphene-Based Materials. Materials 2020, 13, 1032. [Google Scholar] [CrossRef] [Green Version]
  28. Ivanov, A.I.; Nebogatikova, N.A.; Kotin, I.A.; Smagulova, S.A.; Antonova, I.V. Resistive Switching Effects in Fluorinated Graphene Films with Graphene Quantum Dots Enhanced by Polyvinyl Alcohol. Nanotechnology 2019, 30, 255701. [Google Scholar] [CrossRef]
  29. Sun, L.; Gong, P.; Liu, X.; Pang, M.; Tian, M.; Chen, J.; Du, J.; Liu, Z. Fluorinated Carbon Fiber as a Novel Nanocarrier for Cancer Chemo-Photothermal Therapy. J. Mater. Chem. B 2017, 5, 6128. [Google Scholar] [CrossRef]
  30. Romero-Aburto, R.; Narayanan, T.N.; Nagaoka, Y.; Hasumura, T.; Mitcham, T.M.; Fukuda, T.; Cox, P.J.; Bouchard, R.R.; Maekawa, T.; Kumar, D.S.; et al. Fluorinated Graphene Oxide; a New Multimodal Material for Biological Applications. Adv. Mater. 2013, 25, 5632. [Google Scholar] [CrossRef] [Green Version]
  31. Rad, A.S.; Ardjmand, M.; Esfahani, M.R.; Khodashenas, B. DFT Calculations towards the Geometry Optimization, Electronic Structure, Infrared Spectroscopy and UV–Vis Analyses of Favipiravir Adsorption on the First-Row Transition Metals Doped Fullerenes; a New Strategy for COVID-19 Therapy. Spectrochim. Acta A 2021, 247, 119082. [Google Scholar] [CrossRef]
  32. Joshi, S.; Parkar, J.; Ansari, A.; Vora, A.; Talwar, D.; Tiwaskar, M.; Patil, S.; Barkate, H. Role of Favipiravir in the Treatment of COVID-19. Int. J. Infect. Dis. 2021, 102, 501. [Google Scholar] [CrossRef]
  33. Touret, F.; de Lamballerie, X. Of Chloroquine and COVID-19. Antivir. Res. 2020, 177, 104762. [Google Scholar] [CrossRef] [PubMed]
  34. Cortegiani, A.; Ingoglia, G.; Ippolito, M.; Giarratano, A.; Einav, S. A Systematic Review on the Efficacy and Safety of Chloroquine for the Treatment of COVID-19. J. Crit. Care 2020, 57, 279. [Google Scholar] [CrossRef] [PubMed]
  35. Parlak, C.; Alver, Ö.; Şenyel, M. Computational Study on Favipiravir Adsorption onto Undoped- and Silicon-Decorated C60 Fullerenes. J. Theor. Comput. Chem. 2017, 16, 1750011. [Google Scholar] [CrossRef]
  36. Soliman, K.A.; Aal, S.A. Theoretical Investigation of Favipiravir Antiviral Drug Based on Fullerene and Boron Nitride Nanocages. Diam. Relat. Mater. 2021, 117, 108458. [Google Scholar] [CrossRef]
  37. Troyanov, I.S.; Kemnitz, E. Synthesis and Structure of Halogenated Fullerenes. Curr. Org. Chem. 2012, 16, 1060. [Google Scholar] [CrossRef]
  38. Darwish, A.D.; Avent, A.G.; Street, J.M.; Taylor, R. Electrophilic Substitution of C60F18into Phenols: HF Elimination between OH and a 1,3-Shifted Fluorine Giving Benzofurano [2′,3′:10,26]Hexadecafluoro[60]Fullerene and Derivatives. Org. Biomol. Chem. 2003, 1, 1764. [Google Scholar] [CrossRef]
  39. Tsuzuki, S.; Uchimaru, T.; Wakisaka, A.; Ono, T. Magnitude and Directionality of Halogen Bond of Benzene with C6F5X, C6H5X, and CF3X (X = I, Br, Cl, and F). J. Phys. Chem. A 2016, 120, 7020–7029. [Google Scholar] [CrossRef]
  40. Katin, K.P.; Prudkovskiy, V.S.; Maslov, M.M. Molecular Dynamics Simulation of Nickel-coated Graphene Bending. Micro&Nano Lett. 2018, 13, 160. [Google Scholar] [CrossRef]
  41. Salem, M.A.; Katin, K.P.; Kaya, S.; Kochaev, A.I.; Maslov, M.M. Interaction of Dopants and Functional Groups Adsorbed on the Carbon Fullerenes: Computational Study. Phys. E 2020, 124, 114319. [Google Scholar] [CrossRef]
  42. Papoular, R.J.; Allouchi, H.; Dzyabchenko, A.V.; Davydov, V.A.; Rakhmanina, A.V.; Boltalina, O.V.; Seppelt, K.; Agafonov, V. High-Resolution X-Ray Powder Diffraction Structure Determination of C60F48. Fuller. Nanotub. Carb. Nanostruct. 2006, 14, 279. [Google Scholar] [CrossRef]
  43. Hossain, M.R.; Hasan, M.M.; Nishat, M.; Ahmed, F.; Ferdous, T.; Hossain, M.A. DFT and QTAIM Investigations of the Adsorption of Chlormethine Anticancer Drug on the Exterior Surface of Pristine and Transition Metal Functionalized Boron Nitride Fullerene. J. Mol. Liq. 2021, 323, 114627. [Google Scholar] [CrossRef]
  44. Li, M.; Wei, Y.; Zhang, G.; Wang, F.; Li, M.; Soleymanabadi, H. A DFT Study on the Detection of Isoniazid Drug by Pristine, Si and Al Doped C70 Fullerenes. Phys. E 2020, 118, 113878. [Google Scholar] [CrossRef]
  45. Rassolov, V.A.; Pople, J.A.; Ratner, M.A.; Windus, T.L. 6-31G* Basis Set for Atoms K through Zn. J. Chem. Phys. 1998, 109, 1223. [Google Scholar] [CrossRef]
  46. Tanaka, M.; Katouda, M.; Nagase, S. Optimization of RI-MP2 Auxiliary Basis Functions for 6-31G** and 6-311G** Basis Sets for First-, Second-, and Third-Row Elements. J. Comput. Chem. 2013, 34, 2568–2575. [Google Scholar] [CrossRef]
  47. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [Green Version]
  48. Liu, F.; Luehr, N.; Kulik, H.J.; Martínez, T.J. Quantum Chemistry for Solvated Molecules on Graphical Processing Units Using Polarizable Continuum Models. J. Chem. Theory Comput. 2015, 11, 3131–3144. [Google Scholar] [CrossRef] [Green Version]
  49. Seritan, S.; Bannwarth, C.; Fales, B.S.; Hohenstein, E.G.; Kokkila-Schumacher, S.I.L.; Luehr, N.; Snyder, J.W., Jr.; Song, C.; Titov, A.V.; Ufimtsev, I.S.; et al. TeraChem: Accelerating Electronic Structure and ab Initio Molecular Dynamics with Graphical Processing Units. J. Chem. Phys. 2020, 152, 224110. [Google Scholar] [CrossRef]
  50. Wang, L.-P.; Song, C. Geometry Optimization Made Simple with Translation and Rotation Coordinates. J. Chem. Phys. 2020, 144, 214108. [Google Scholar] [CrossRef] [Green Version]
  51. Yanai, T.; Tew, D.P.; Handy, N.C. A New Hybrid Exchange–Correlation Functional Using the Coulomb-Attenuating Method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51. [Google Scholar] [CrossRef] [Green Version]
  52. Isborn, C.M.; Luehr, N.; Ufimtsev, I.S.; Martínez, T.J. Excited-State Electronic Structure with Configuration Interaction Singles and Tamm–Dancoff Time-Dependent Density Functional Theory on Graphical Processing Units. J. Chem. Theory Comput. 2011, 7, 1814. [Google Scholar] [CrossRef]
  53. Timoshnikov, V.A.; Kobzeva, T.V.; Polyakov, N.E.; Kontoghiorghes, G.J. Redox Interactions of Vitamin C and Iron: Inhibition of the Pro-Oxidant Activity by Deferiprone. Int. J. Mol. Sci. 2020, 21, 3967. [Google Scholar] [CrossRef] [PubMed]
  54. Smidstrup, S.; Markussen, T.; Vancraeyveld, P.; Wellendorff, J.; Schneider, J.; Gunst, T.; Verstichel, B.; Stradi, D.; Khomyakov, P.A.; Vej-Hansen, U.G.; et al. Quantum ATK: An integrated platform of electronic and atomic-scale modelling tools. J. Phys.: Cond. Mat. 2020, 32, 015901. [Google Scholar] [CrossRef]
Figure 1. Chemical structures of (a) favipiravir and (b) chloroquine molecules.
Figure 1. Chemical structures of (a) favipiravir and (b) chloroquine molecules.
Ijms 23 02345 g001
Figure 2. Optimized geometry of the chloroquine molecule loaded on fluorinated fullerenes (a) C60F2 and (b) C60F48.
Figure 2. Optimized geometry of the chloroquine molecule loaded on fluorinated fullerenes (a) C60F2 and (b) C60F48.
Ijms 23 02345 g002
Figure 3. Optimized geometries of pristine, low- and high-fluorinated C60 decorated with metal ions: (a) C60Cr3+, (b) C60F2Cr3+, (c) C60F48Cr3+, (d) C60Fe2+, (e) C60F2Fe2+, (f) C60F48Fe2+, (g) C60Fe3+, (h) C60F2Fe3+, (i) C60F48Fe3+, (j) C60Ni2+, (k) C60F2Ni2+, (l) C60F48Ni2+.
Figure 3. Optimized geometries of pristine, low- and high-fluorinated C60 decorated with metal ions: (a) C60Cr3+, (b) C60F2Cr3+, (c) C60F48Cr3+, (d) C60Fe2+, (e) C60F2Fe2+, (f) C60F48Fe2+, (g) C60Fe3+, (h) C60F2Fe3+, (i) C60F48Fe3+, (j) C60Ni2+, (k) C60F2Ni2+, (l) C60F48Ni2+.
Ijms 23 02345 g003
Figure 4. Optimized geometries of chloroquine drug loaded on Ni2+-decorated fluorinated fullerenes (a) C60Ni2+, (b) C60F2Ni2+ and (c) C60F48Ni2+.
Figure 4. Optimized geometries of chloroquine drug loaded on Ni2+-decorated fluorinated fullerenes (a) C60Ni2+, (b) C60F2Ni2+ and (c) C60F48Ni2+.
Ijms 23 02345 g004aIjms 23 02345 g004b
Figure 5. Ultraviolet and visible spectra of (a) chloroquine drug (black line), C60F2Ni2+ carrier (blue line) and carrier + loaded drug complex (red dashed line). (b) The same spectra for high-fluorinated C60F48Ni2+ carrier loaded with chloroquine. The calculated wavelengths corresponding to the transitions are broadened by Gaussian curves with σ = 10 nm.
Figure 5. Ultraviolet and visible spectra of (a) chloroquine drug (black line), C60F2Ni2+ carrier (blue line) and carrier + loaded drug complex (red dashed line). (b) The same spectra for high-fluorinated C60F48Ni2+ carrier loaded with chloroquine. The calculated wavelengths corresponding to the transitions are broadened by Gaussian curves with σ = 10 nm.
Ijms 23 02345 g005
Figure 6. Infrared spectra of (a) chloroquine drug (black line), C60F2Ni2+ carrier (blue line) and carrier + loaded drug complex (red dashed line). (b) The same spectra for high-fluorinated C60F48Ni2+ carrier loaded with chloroquine. The calculated wavelengths corresponding to the transitions are broadened by Gaussian curves with σ = 10 cm−1.
Figure 6. Infrared spectra of (a) chloroquine drug (black line), C60F2Ni2+ carrier (blue line) and carrier + loaded drug complex (red dashed line). (b) The same spectra for high-fluorinated C60F48Ni2+ carrier loaded with chloroquine. The calculated wavelengths corresponding to the transitions are broadened by Gaussian curves with σ = 10 cm−1.
Ijms 23 02345 g006
Table 1. Calculated geometric and energetic characteristics of fluorinated C60 loaded by chloroquine drug. Binding energy Eb (eV), fullerene and drug deformation energies Edef (eV), interaction energy Eint (eV), dipole moments D (Debye), frontier molecular orbitals HOMO and LUMO (eV) and HOMO–LUMO gaps (eV) are presented.
Table 1. Calculated geometric and energetic characteristics of fluorinated C60 loaded by chloroquine drug. Binding energy Eb (eV), fullerene and drug deformation energies Edef (eV), interaction energy Eint (eV), dipole moments D (Debye), frontier molecular orbitals HOMO and LUMO (eV) and HOMO–LUMO gaps (eV) are presented.
EbEdef (Fullerene)Edef (Drug)EintDHOMOLUMOGap
C60F203.95−5.86−3.342.52
C60F2 + drug0.410.010.040.465.61−5.54−3.282.26
C60F4800.17−9.51−4.155.36
C60F48 + drug0.290.030.060.387.01−5.37−4.161.21
Table 2. Calculated geometric and energetic characteristics of metal-decorated pristine and fluorinated fullerenes. Binding energy Eb (eV), fullerene deformation energy Edef (eV), mean length of carbon-ion bonds lC–M (Å) and the number of such bonds n (dimensionless), mean length of fluorine-ion bonds lM–F (Å) and the number of such bonds m (dimensionless), dipole moments D (Debye), frontier molecular orbitals HOMO and LUMO (eV) and HOMO–LUMO gaps (eV) are presented.
Table 2. Calculated geometric and energetic characteristics of metal-decorated pristine and fluorinated fullerenes. Binding energy Eb (eV), fullerene deformation energy Edef (eV), mean length of carbon-ion bonds lC–M (Å) and the number of such bonds n (dimensionless), mean length of fluorine-ion bonds lM–F (Å) and the number of such bonds m (dimensionless), dipole moments D (Debye), frontier molecular orbitals HOMO and LUMO (eV) and HOMO–LUMO gaps (eV) are presented.
SystemEbEdeflC–M (n)lM–F (m)DHOMOLUMOGap
non-fluorinated fullerene
C60Cr3+17.940.182.265 (3)18.79−7.08−6.140.94
C60Fe2+6.260.392.013 (2)14.85−6.85−5.831.02
C60Fe3+21.710.352.034 (2)15.17−7.14−6.200.94
C60Ni2+8.690.441.928 (2)9.45−6.91−6.000.91
low-fluorinated fullerene
C60F2Cr3+21.972.150 (1)1.684 (2)1.87−7.13−6.900.23
C60F2Fe2+8.542.046 (1)1.836 (2)2.84−5.96−5.030.93
C60F2Fe3+22.671.161.994 (1)1.866 (1)12.18−7.14−6.290.85
C60F2Ni2+9.040.891.924 (1)1.895 (1)10.98−6.93−5.791.14
high-fluorinated fullerene
C60F48Cr3+14.022.302.023 (3)29.81−9.37−7.941.43
C60F48Fe2+4.641.061.918 (3)39.71−9.92−8.301.62
C60F48Fe3+17.012.621.921 (3)26.77−9.40−8.470.93
C60F48Ni2+3.890.551.930 (2)39.48−9.92−9.560.36
Table 3. Calculated geometric and energetic characteristics of metal-decorated pristine and fluorinated C60 loaded by chloroquine drug. Binding energy Eb (eV), fullerene and drug deformation energies Edef (eV), bond lengths between the metal ion and nitrogen atom lM-N (Å), dipole moments D (Debye), frontier molecular orbitals HOMO and LUMO (eV) and HOMO–LUMO gaps (eV) are presented.
Table 3. Calculated geometric and energetic characteristics of metal-decorated pristine and fluorinated C60 loaded by chloroquine drug. Binding energy Eb (eV), fullerene and drug deformation energies Edef (eV), bond lengths between the metal ion and nitrogen atom lM-N (Å), dipole moments D (Debye), frontier molecular orbitals HOMO and LUMO (eV) and HOMO–LUMO gaps (eV) are presented.
EbEdef (Fullerene)Edef (Drug)lM-NDHOMOLUMOGap
non-fluorinated fullerene
C60Cr3+ + drug5.210.160.791.96839.53−7.04−4.422.62
C60Fe2+ + drug3.740.250.681.90420.39−5.18−4.550.63
C60Fe3+ + drug5.200.230.841.90939.81−6.58−4.671.91
C60Ni2+ + drug2.610.230.321.85328.88−5.66−5.330.34
low-fluorinated fullerene
C60F2Fe3+ + drug4.850.320.431.90834.46−7.01−5.241.77
C60F2Ni2+ + drug3.430.040.241.85215.30−6.03−5.070.96
high-fluorinated fullerene
C60F48Cr3+ + drug8.010.870.651.89568.71−6.41−6.270.14
C60F48Fe2+ + drug4.780.870.581.87544.13−5.38−4.970.41
C60F48Fe3+ + drug9.321.060.801.86868.96−7.10−5.321.78
C60F48Ni2+ + drug5.880.470.421.84646.82−5.90−5.460.45
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Katin, K.P.; Kochaev, A.I.; Kaya, S.; El-Hajjaji, F.; Maslov, M.M. Ab Initio Insight into the Interaction of Metal-Decorated Fluorinated Carbon Fullerenes with Anti-COVID Drugs. Int. J. Mol. Sci. 2022, 23, 2345. https://doi.org/10.3390/ijms23042345

AMA Style

Katin KP, Kochaev AI, Kaya S, El-Hajjaji F, Maslov MM. Ab Initio Insight into the Interaction of Metal-Decorated Fluorinated Carbon Fullerenes with Anti-COVID Drugs. International Journal of Molecular Sciences. 2022; 23(4):2345. https://doi.org/10.3390/ijms23042345

Chicago/Turabian Style

Katin, Konstantin P., Alexey I. Kochaev, Savas Kaya, Fadoua El-Hajjaji, and Mikhail M. Maslov. 2022. "Ab Initio Insight into the Interaction of Metal-Decorated Fluorinated Carbon Fullerenes with Anti-COVID Drugs" International Journal of Molecular Sciences 23, no. 4: 2345. https://doi.org/10.3390/ijms23042345

APA Style

Katin, K. P., Kochaev, A. I., Kaya, S., El-Hajjaji, F., & Maslov, M. M. (2022). Ab Initio Insight into the Interaction of Metal-Decorated Fluorinated Carbon Fullerenes with Anti-COVID Drugs. International Journal of Molecular Sciences, 23(4), 2345. https://doi.org/10.3390/ijms23042345

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop