ISC, DOAJ, CAS, Google Scholar......

Molecular Dynamics Simulation and Adsorption Study of Different Anticoagulants on Fullerene

Document Type : Research Paper

Authors

1 Department of Chemistry, Payame Noor University, Tehran, Iran.

2 Department of Chemistry, Shahid Bahonar University of Kerman, Kerman, Iran

3 Medical Biotechnology Research Center, Ashkezar Branch, Islamic Azad University, Yazd, Iran

10.22034/ijnc.2024.720190
Abstract
The different properties of fullerenes in drug delivery are the subject of various discussions. In this research, we investigated the structure of anticoagulants (Apixaban, Dabigartan, Edoxaban, Rivaroxaban, and Warfarin with fullerenes to provide a drug delivery system. The structure of the anticoagulants and the two-dimensional structure of C60 fullerene were obtained from the PubChem database. With the Spartan software series, their structure was optimized and three-dimensional with the B3LYP 6-311G* DFT function. Then we used the ChemOffice collection to obtain the structure with the lowest energy and to calculate the HOMO and LOMO energy and the total energy of the structure in question; we used the SPARTAN software and the DFT function B3LYP 6-311G*. Fullerene drug complex calculations showed that all reactions are exothermic, and experimentally, there is the possibility of forming a drug bond with fullerene. In all drugs, after the formation of the drug-fullerene complex, the amount of electrophilicity and the maximum transfer charge to the system increased, which means that the drug-fullerene complex is more electrophilic. Also, after the bonding between C60 and drugs in all compounds, their hardness decreased, and their chemical softness increased, proving the reaction and bonding in experimental conditions. The HOMO and LUMO calculations done in the MD study indicate that there is a possibility to examine the electrostatics between the compounds and C60. The reactivity of the drugs is lower than that of C60, making the system stable if drugs interact with C60. The total and adsorption energy results showed high interaction energy between Db-C60, Wr-C60, Rv-C60, Ed-C60, and Ap-C60, respectively. The interactions of Ap, Db, Ed, Rv, and Wr towards C60 are of the Van der Waals type, making C60 an effective carrier for these drugs.

Keywords

Subjects

[1]        B. H. Al-Tamimi, S. B.H. Farid, Fullerenes and Nanodiamonds for Medical Drug Delivery, in Nanocrystals [Working Title], (2021).
[2]        I. V. Mikheev, M. M. Sozarukova, D. Y. Izmailov, I. E. Kareev, E. V. Proskurnina, M. A. Proskurnin, Antioxidant potential of aqueous dispersions of fullerenes c60, c70, and gd@c82, Int. J. Mol. Sci., 22, 11, (2021).
[3]        M. Kumar and K. Raza, C60-fullerenes as Drug Delivery Carriers for Anticancer Agents: Promises and Hurdles, Pharm. Nanotechnol., 5, 3, (2017).
[4]        A. Montellano, T. Da Ros, A. Bianco, M. Prato, Fullerene C60 as a multifunctional system for drug and gene delivery, Nanoscale, 3, 10, (2011).
[5]        C. H. Yeh, K. Hogg, J. I. Weitz, Overview of the new oral anticoagulants: Opportunities and challenges,  Arterioscler. Thromb. Vasc. Biol., 35, 5, (2015).
[6]        Y. H. Mekaj, A. Y. Mekaj, S. B. Duci, E. I. Miftari, New oral anticoagulants: Their advantages and disadvantages compared with vitamin K antagonists in the prevention and treatment of patients with thromboembolic events, Ther. Clin. Risk Manag., 11, (2015).
[7]        R. Bauersachs, Edoxaban: a new oral direct factor Xa inhibitor for the prevention and treatment of thromboembolic disorders, Clin. Investig. (Lond)., vol. 4, no. 7, (2014).
[8]        E. Shakerzadeh, Theoretical investigations of interactions between boron nitride nanotubes and drugs, in Boron Nitride Nanotubes in Nanomedicine, (2016).
[9]        C. J. Cramer, Essentials of Computational Chemistry Theories and Models, 42, 2. (2004).
[10]      G. B. Goh, N. O. Hodas, and A. Vishnu, Deep learning for computational chemistry, Journal of Computational Chemistry, 38, 16. (2017).
[11]      R. Kumari, R. Rathi, S. R. Pathak, and V. Dalal, Computational investigation of potent inhibitors against YsxC: structure-based pharmacophore modeling, molecular docking, molecular dynamics, and binding free energy, J. Biomol. Struct. Dyn., (2021).
[12]      S. McArdle, S. Endo, A. Aspuru-Guzik, S. C. Benjamin, X. Yuan, Quantum computational chemistry, Rev. Mod. Phys., 92, 1, (2020).
[13]      W. Li, H. Ma, S. Li, J. Ma, Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning, Chemical Science, 12, 45. (2021).
[14]      E. G. Lewars, Computational chemistry: Introduction to the theory and applications of molecular and quantum mechanics. (2011).
[15]      H. Bukheet Hassan Ȧ, Density Function Theory B3LYP/6-31G**Calculation of Geometry Optimization and Energies of Donor-Bridge-Acceptor Molecular System, Res. Artic. Int. J. Curr. Eng. Technol. Accept., 4, 4, (2014).
[16]      R. Pelalak et al., Molecular dynamics simulation of novel diamino-functionalized hollow mesosilica spheres for adsorption of dyes from synthetic wastewater,” J. Mol. Liq., 322, (2021).
[17]      R. Eshaghi Malekshah, B. Fahimirad, M. Aallaei, A. Khaleghian, Synthesis and toxicity assessment of Fe3O4 NPs grafted by ∼ NH2-Schiff base as anticancer drug: modeling and proposed molecular mechanism through docking and molecular dynamic simulation,  Drug Deliv., 27, 1, (2020).
[18]      P. Pourhakkak, A. Taghizadeh, M. Taghizadeh, M. Ghaedi, S. Haghdoust, Fundamentals of adsorption technology, in Interface Science and Technology, 33, (2021).
[19]      T. Lemaoui et al., Prediction of Electrical Conductivity of Deep Eutectic Solvents Using COSMO-RS Sigma Profiles as Molecular Descriptors: A Quantitative Structure-Property Relationship Study, Ind. Eng. Chem. Res., 59, 29, (2020).
[20]      R. Pelalak et al., “Synthesis, molecular dynamics simulation and adsorption study of different pollutants on functionalized mesosilica,” Sci. Rep., 11, 1, (2021).
[21]      A. Ghasemi, M. Asgarpour Khansary, A. Marjani, S. Shirazian, Using quantum chemical modeling and calculations for evaluation of cellulose potential for estrogen micropollutants removal from water effluents, Chemosphere, 178, (2017).
[22]      W. Gao, Y. Chen, B. Li, S. P. Liu, X. Liu, Q. Jiang, Determining the adsorption energies of small molecules with the intrinsic properties of adsorbates and substrates, Nat. Commun., 11, 1, (2020).
[23]      J. R. De Lile, S. G. Kang, Y. A. Son, S. G. Lee, Do HOMO-LUMO energy levels and band gaps provide sufficient understanding of Dye-sensitizer activity trends for water purification?,  ACS Omega, 5, 25, (2020).
[24]      M. K. Harbola, P. K. Chattaraj, R. G. Parr, Aspects of the Softness and Hardness Concepts of Density‐Functional Theory, Isr. J. Chem., 31, 4, (1991).
[25]      J. L. Gázquez, Hardness and softness in density functional theory, in Chemical Hardness, (2006).
[26]      İ. Muz, F. Göktaş, M. Kurban, Size dependence in the electronic and optical properties of a BN analogue of two-dimensional graphdiyne: A theoretical study, Chem. Phys., 539, (2020).
[27]      R. G. Parr, L. V. Szentpály, S. Liu, Electrophilicity index, J. Am. Chem. Soc., 121, 9, (1999).