1
|
Laref S, Harrou F, Sun Y, Gao X, Gojobori T. Exploring Antiviral Drugs on Monolayer Black Phosphorene: Atomistic Theory and Explainable Machine Learning-Assisted Platform. Int J Mol Sci 2024; 25:4897. [PMID: 38732115 PMCID: PMC11084629 DOI: 10.3390/ijms25094897] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Revised: 04/24/2024] [Accepted: 04/26/2024] [Indexed: 05/13/2024] Open
Abstract
Favipiravir (FP) and ebselen (EB) belong to a diverse class of antiviral drugs known for their significant efficacy in treating various viral infections. Utilizing molecular dynamics (MD) simulations, machine learning, and van der Waals density functional theory, we accurately elucidate the binding properties of these antiviral drugs on a phosphorene single-layer. To further investigate these characteristics, this study employs four distinct machine learning models-Random Forest, Gradient Boosting, XGBoost, and CatBoost. The Hamiltonian of antiviral molecules within a monolayer of phosphorene is appropriately trained. The key aspect of utilizing machine learning (ML) in drug design revolves around training models that are efficient and precise in approximating density functional theory (DFT). Furthermore, the study employs SHAP (SHapley Additive exPlanations) to elucidate model predictions, providing insights into the contribution of each feature. To explore the interaction characteristics and thermodynamic properties of the hybrid drug, we employ molecular dynamics and DFT calculations in a vacuum interface. Our findings suggest that this functionalized 2D complex exhibits robust thermostability, indicating its potential as an effective and enabled entity. The observed variations in free energy at different surface charges and temperatures suggest the adsorption potential of FP and EB molecules from the surrounding environment.
Collapse
Affiliation(s)
- Slimane Laref
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia; (X.G.); (T.G.)
| | - Fouzi Harrou
- Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia;
| | - Ying Sun
- Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia;
| | - Xin Gao
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia; (X.G.); (T.G.)
| | - Takashi Gojobori
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia; (X.G.); (T.G.)
| |
Collapse
|
2
|
Laref S, Harrou F, Wang B, Sun Y, Laref A, Laleg-Kirati TM, Gojobori T, Gao X. Synergy of Small Antiviral Molecules on a Black-Phosphorus Nanocarrier: Machine Learning and Quantum Chemical Simulation Insights. Molecules 2023; 28:molecules28083521. [PMID: 37110754 PMCID: PMC10142408 DOI: 10.3390/molecules28083521] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2023] [Revised: 04/03/2023] [Accepted: 04/10/2023] [Indexed: 04/29/2023] Open
Abstract
Favipiravir (FP) and Ebselen (EB) belong to a broad range of antiviral drugs that have shown active potential as medications against many viruses. Employing molecular dynamics simulations and machine learning (ML) combined with van der Waals density functional theory, we have uncovered the binding characteristics of these two antiviral drugs on a phosphorene nanocarrier. Herein, by using four different machine learning models (i.e., Bagged Trees, Gaussian Process Regression (GPR), Support Vector Regression (SVR), and Regression Trees (RT)), the Hamiltonian and the interaction energy of antiviral molecules in a phosphorene monolayer are trained in an appropriate way. However, training efficient and accurate models for approximating the density functional theory (DFT) is the final step in using ML to aid in the design of new drugs. To improve the prediction accuracy, the Bayesian optimization approach has been employed to optimize the GPR, SVR, RT, and BT models. Results revealed that the GPR model obtained superior prediction performance with an R2 of 0.9649, indicating that it can explain 96.49% of the data's variability. Then, by means of DFT calculations, we examine the interaction characteristics and thermodynamic properties in a vacuum and a continuum solvent interface. These results illustrate that the hybrid drug is an enabled, functionalized 2D complex with vigorous thermostability. The change in Gibbs free energy at different surface charges and temperatures implies that the FP and EB molecules are allowed to adsorb from the gas phase onto the 2D monolayer at different pH conditions and high temperatures. The results reveal a valuable antiviral drug therapy loaded by 2D biomaterials that may possibly open a new way of auto-treating different diseases, such as SARS-CoV, in primary terms.
Collapse
Affiliation(s)
- Slimane Laref
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science & Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Fouzi Harrou
- A Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Bin Wang
- Center for Interfacial Reaction Engineering (CIRE), School of Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, OK 73019, USA
| | - Ying Sun
- A Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Amel Laref
- Department of Physics and Astronomy, College of Science, King Saud University, Riyadh 11451, Saudi Arabia
| | - Taous-Meriem Laleg-Kirati
- A Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Takashi Gojobori
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science & Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Xin Gao
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science & Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
3
|
Laref S, Wang B, Gao X, Gojobori T. Computational Studies of Auto-Active van der Waals Interaction Molecules on Ultra-Thin Black-Phosphorus Film. Molecules 2023; 28:molecules28020681. [PMID: 36677738 PMCID: PMC9864666 DOI: 10.3390/molecules28020681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Revised: 01/03/2023] [Accepted: 01/04/2023] [Indexed: 01/11/2023] Open
Abstract
Using the van der Waals density functional theory, we studied the binding peculiarities of favipiravir (FP) and ebselen (EB) molecules on a monolayer of black phosphorene (BP). We systematically examined the interaction characteristics and thermodynamic properties in a vacuum and a continuum, solvent interface for active drug therapy. These results illustrate that the hybrid molecules are enabled functionalized two-dimensional (2D) complex systems with a vigorous thermostability. We demonstrate in this study that these molecules remain flat on the monolayer BP system and phosphorus atoms are intact. It is inferred that the hybrid FP+EB molecules show larger adsorption energy due to the van der Waals forces and planar electrostatic interactions. The changes in Gibbs free energy at different surface charge fluctuations and temperatures imply that the FP and EB are allowed to adsorb from the gas phase onto the 2D film at high temperatures. Thereby, the results unveiled beneficial inhibitor molecules on two dimensional BP nanocarriers, potentially introducing a modern strategy to enhance the development of advanced materials, biotechnology, and nanomedicine.
Collapse
Affiliation(s)
- Slimane Laref
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science & Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
- Correspondence:
| | - Bin Wang
- School of Chemical, Biological and Materials Engineering, Center for Interfacial Reaction Engineering (CIRE), University of Oklahoma, Norman, OK 73019, USA
| | - Xin Gao
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science & Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Takashi Gojobori
- Computational Bioscience Research Center (CBRC), King Abdullah University of Science & Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
4
|
Rekhviashvili SS, Bukhurova MM. Model of Interaction of Fullerene С60 with Epitaxial Graphene. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2018. [DOI: 10.1134/s0036024418100266] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
5
|
Jaiswal R, Saha U, Goswami TH, Srivastava A, Prasad N. ‘Pillar effect’ of chemically bonded fullerene in enhancing supercapacitance performances of partially reduced fullerenol graphene oxide hybrid electrode material. Electrochim Acta 2018. [DOI: 10.1016/j.electacta.2018.06.110] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
6
|
Jha SK, Roth M, Todde G, Buchanan JP, Moser RD, Shukla MK, Subramanian G. Role of Stone-Wales defects on the interfacial interactions among graphene, carbon nanotubes, and Nylon 6: A first-principles study. J Chem Phys 2018; 149:054703. [PMID: 30089374 PMCID: PMC6910596 DOI: 10.1063/1.5032081] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2018] [Accepted: 07/12/2018] [Indexed: 12/29/2022] Open
Abstract
We investigate computationally the role of Stone-Wales (SW) defects on the interfacial interactions among graphene, carbon nanotubes (CNTs), and Nylon 6 using density functional theory (DFT) and the empirical force-field. Our first-principles DFT calculations were performed using the Quantum ESPRESSO electronic structure code with the highly accurate van der Waals functional (vdW-DF2). Both pristine and SW-defected carbon nanomaterials were investigated. The computed results show that the presence of SW defects on CNTs weakens the CNT-graphene interactions. Our result that CNT-graphene interaction is much stronger than CNT-CNT interaction indicates that graphene would be able to promote the dispersion of CNTs in the polymer matrix. Our results demonstrate that carbon nanomaterials form stable complexes with Nylon 6 and that the van der Waals interactions, as revealed by the electronic charge density difference maps, play a key stabilizing role on the interfacial interactions among graphene, CNTs, and Nylon 6. Using the density of states calculations, we observed that the bandgaps of graphene and CNTs were not significantly modified due to their interactions with Nylon 6. The Young's moduli of complexes were found to be the averages of the moduli of their individual constituents.
Collapse
Affiliation(s)
- Sanjiv K. Jha
- Authors to whom correspondence should be addressed: and
| | - Michael Roth
- School of Polymer Science and Engineering, The University of Southern Mississippi, Hattiesburg, Mississippi 39406, USA
| | - Guido Todde
- Department of Chemistry and Biomedical Sciences, Linnaeus University, Norra Vägen 49, 392 34 Kalmar, Sweden
| | - J. Paige Buchanan
- Concrete and Materials Branch, Engineer Research and Development Center, US Army Corps of Engineers, 3909 Halls Ferry Road, Vicksburg, Mississippi 39180, USA
| | - Robert D. Moser
- Concrete and Materials Branch, Engineer Research and Development Center, US Army Corps of Engineers, 3909 Halls Ferry Road, Vicksburg, Mississippi 39180, USA
| | - Manoj K. Shukla
- Environmental Laboratory, Engineer Research and Development Center, US Army Corps of Engineers, 3909 Halls Ferry Road, Vicksburg, Mississippi 39180, USA
| | - Gopinath Subramanian
- School of Polymer Science and Engineering, The University of Southern Mississippi, Hattiesburg, Mississippi 39406, USA
| |
Collapse
|
7
|
Borowik A, Prylutskyy Y, Kawelski Ł, Kyzyma O, Bulavin L, Ivankov O, Cherepanov V, Wyrzykowski D, Kaźmierkiewicz R, Gołuński G, Woziwodzka A, Evstigneev M, Ritter U, Piosik J. Does C 60 fullerene act as a transporter of small aromatic molecules? Colloids Surf B Biointerfaces 2018; 164:134-143. [PMID: 29413590 DOI: 10.1016/j.colsurfb.2018.01.026] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 12/23/2017] [Accepted: 01/15/2018] [Indexed: 12/13/2022]
Abstract
C60 fullerene is reported to directly interact with biomolecules, such as aromatic mutagens or anticancer drugs. Therefore, it is extensively studied for its potential application in the fields of drug delivery and chemoprevention. Understanding the nature of fullerene-drugs interactions might contribute to optimization and modification of the existing chemotherapy systems. Possible interactions between ICR-191, a model acridine mutagen, with well-established biophysical properties and mutagenic activity, and C60 fullerene aqueous solution were investigated by broad range of biophysical methods, such as Dynamic Light Scattering, Isothermal Titration Calorimetry, and Atomic Force Microscopy. Additionally, to determine biological activity of ICR-191-C60 fullerene mixtures, Ames mutagenicity test was employed. It was demonstrated that C60 fullerene interacts non-covalently with ICR-191 and has strong affinity to bacterial membranes. The obtained results provide practical insight into C60 fullerene interactions with aromatic compounds.
Collapse
Affiliation(s)
- Agnieszka Borowik
- Laboratory of Biophysics, Intercollegiate Faculty of Biotechnology UG-MUG, Abrahama 58, 80-307 Gdańsk, Poland
| | - Yuriy Prylutskyy
- Taras Shevchenko National University of Kyiv, Volodymyrska Str., 64, 01601 Kyiv, Ukraine
| | | | - Olena Kyzyma
- Taras Shevchenko National University of Kyiv, Volodymyrska Str., 64, 01601 Kyiv, Ukraine; Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia
| | - Leonid Bulavin
- Taras Shevchenko National University of Kyiv, Volodymyrska Str., 64, 01601 Kyiv, Ukraine
| | - Oleksandr Ivankov
- Taras Shevchenko National University of Kyiv, Volodymyrska Str., 64, 01601 Kyiv, Ukraine; Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia
| | | | | | - Rajmund Kaźmierkiewicz
- Laboratory of Biomolecular Systems Simulations, Intercollegiate Faculty of Biotechnology UG-MUG, Abrahama 58, 80-307 Gdańsk, Poland
| | - Grzegorz Gołuński
- Laboratory of Biophysics, Intercollegiate Faculty of Biotechnology UG-MUG, Abrahama 58, 80-307 Gdańsk, Poland
| | - Anna Woziwodzka
- Laboratory of Biophysics, Intercollegiate Faculty of Biotechnology UG-MUG, Abrahama 58, 80-307 Gdańsk, Poland
| | - Maxim Evstigneev
- Department of Physics, Sevastopol State University, Sevastopol, Crimea, Russia
| | - Uwe Ritter
- Technical University of Ilmenau, Institute of Chemistry and Biotechnology, 25 Weimarer Str., 98693 Ilmenau, Germany
| | - Jacek Piosik
- Laboratory of Biophysics, Intercollegiate Faculty of Biotechnology UG-MUG, Abrahama 58, 80-307 Gdańsk, Poland.
| |
Collapse
|
8
|
Correa JD, Orellana PA, Pacheco M. Optoelectronic Properties of Van Der Waals Hybrid Structures: Fullerenes on Graphene Nanoribbons. NANOMATERIALS 2017; 7:nano7030069. [PMID: 28336904 PMCID: PMC5388171 DOI: 10.3390/nano7030069] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/18/2016] [Revised: 02/16/2017] [Accepted: 03/01/2017] [Indexed: 01/26/2023]
Abstract
The search for new optical materials capable of absorbing light in the frequency range from visible to near infrared is of great importance for applications in optoelectronic devices. In this paper, we report a theoretical study of the electronic and optical properties of hybrid structures composed of fullerenes adsorbed on graphene and on graphene nanoribbons. The calculations are performed in the framework of the density functional theory including the van der Waals dispersive interactions. We found that the adsorption of the C 60 fullerenes on a graphene layer does not modify its low energy states, but it has strong consequences for its optical spectrum, introducing new absorption peaks in the visible energy region. The optical absorption of fullerenes and graphene nanoribbon composites shows a strong dependence on photon polarization and geometrical characteristics of the hybrid systems, covering a broad range of energies. We show that an external electric field across the nanoribbon edges can be used to tune different optical transitions coming from nanoribbon-fullerene hybridized states, which yields a very rich electro-absorption spectrum for longitudinally polarized photons. We have carried out a qualitative analysis on the potential of these hybrids as possible donor-acceptor systems in photovoltaic cells.
Collapse
Affiliation(s)
- Julián David Correa
- Departamento de Ciencias Básicas, Universidad de Medellín, 050026 Medellín, Colombia.
| | | | - Mónica Pacheco
- Departamento de Física, Universidad Técnica Federico Santa María, 2340000 Valparaíso, Chile.
| |
Collapse
|
9
|
Iglesias D, Guerra J, Gómez MV, Rodríguez AM, Prieto P, Vázquez E, Herrero MA. Design of Assembled Systems Based on Conjugated Polyphenylene Derivatives and Carbon Nanohorns. Chemistry 2016; 22:11643-51. [PMID: 27404562 DOI: 10.1002/chem.201601282] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2016] [Indexed: 11/11/2022]
Abstract
Promising materials have been designed and fully characterised by an effective interaction between versatile platforms such as carbon nanohorns (CNHs) and conjugated molecules based on thiophene derivatives. Easy and non-aggressive methods have been described for the synthesis and purification of the final systems. Oligothiophenephenylvinylene (OTP) systems with different geometries and electron density are coupled to the CNHs. A wide range of characterization techniques have been used to confirm the effective interaction between the donor (OTP) and the acceptor (CNH) systems. These hybrid materials show potential for integration into solar cell devices. Importantly, surface-enhanced Raman spectroscopy (SERS) effects are observed without the presence of any metal surface in the system. Theoretical calculations have been performed to study the optimised geometries of the noncovalent interaction between the surface and the organic molecule. The calculations allow information on the monoelectronic energies of HOMO-LUMO orbitals and band gap of different donor systems to be extracted.
Collapse
Affiliation(s)
- Daniel Iglesias
- Departamento de Química Orgánica, Inorgánica y Bioquímica, Facultad de Ciencias y Tecnologías Químicas, IRICA Universidad de Castilla-La Mancha, Campus Universitario, 13071, Ciudad Real, Spain.,Dipartimento di Scienze Chimiche e Farmaceutiche, Center of Excellence for Nanostructured Materials (CENMAT) & Italian Interuniversity Consortium on Materials Science and Technology (INSTM - Unit of Trieste), Università degli Studi di Trieste, Piazzale Europa 1, 34127, Trieste, Italy
| | - Javier Guerra
- Departamento de Química Orgánica, Inorgánica y Bioquímica, Facultad de Ciencias y Tecnologías Químicas, IRICA Universidad de Castilla-La Mancha, Campus Universitario, 13071, Ciudad Real, Spain.,Crystal Pharma, a division of AMRI, Parque Tecnológico de Boecillo, Valladolid, 47151, Spain
| | - M Victoria Gómez
- Departamento de Química Orgánica, Inorgánica y Bioquímica, Facultad de Ciencias y Tecnologías Químicas, IRICA Universidad de Castilla-La Mancha, Campus Universitario, 13071, Ciudad Real, Spain
| | - Antonio M Rodríguez
- Departamento de Química Orgánica, Inorgánica y Bioquímica, Facultad de Ciencias y Tecnologías Químicas, IRICA Universidad de Castilla-La Mancha, Campus Universitario, 13071, Ciudad Real, Spain.,Università degli Studi di Napoli Federico II - Dipartimento di Scienze Chimiche, Via Cintia 26, Monte Sant'Angelo, 80126, Napoli Campania, Italy
| | - Pilar Prieto
- Departamento de Química Orgánica, Inorgánica y Bioquímica, Facultad de Ciencias y Tecnologías Químicas, IRICA Universidad de Castilla-La Mancha, Campus Universitario, 13071, Ciudad Real, Spain
| | - Ester Vázquez
- Departamento de Química Orgánica, Inorgánica y Bioquímica, Facultad de Ciencias y Tecnologías Químicas, IRICA Universidad de Castilla-La Mancha, Campus Universitario, 13071, Ciudad Real, Spain
| | - M Antonia Herrero
- Departamento de Química Orgánica, Inorgánica y Bioquímica, Facultad de Ciencias y Tecnologías Químicas, IRICA Universidad de Castilla-La Mancha, Campus Universitario, 13071, Ciudad Real, Spain.
| |
Collapse
|
10
|
Buckminster fullerene adhesion on graphene flakes: Numerical accuracy of dispersion corrected DFT. Polyhedron 2016. [DOI: 10.1016/j.poly.2015.11.016] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|