Molecular Dynamics Insights into the Adsorption of the COVID-19 Antiviral Remdesivir on Silica-Functionalized Graphene Oxide: Enthalpic Prevalence and Comparative Evaluation for Targeted Antiviral Delivery
نویسندگان
1 Department of Biological and Chemical Sciences, School of Arts and Sciences, Lebanese International University, Beirut, Lebanon- PO Box: 146404 Mazraa, Lebanon
2 Department of Biological and Chemical Sciences, School of Arts and Sciences, Lebanese International University, Beirut, Lebanon- PO Box: 146404 Mazraa, Lebanon
3 Department of Nutrition and Food Sciences, School of Arts and Sciences, Lebanese International University, Mouseitbeh, P.O. Box 146404, Mazraa, Beirut, Lebanon
4 Department of Biological and Chemical Sciences, School of Arts and Sciences, Lebanese International University, Beirut, Lebanon- PO Box: 146404 Mazraa, Lebanon
5 Department of Biological and Chemical Sciences, School of Arts and Sciences, Lebanese International University, Beirut, Lebanon- PO Box: 146404 Mazraa, Lebanon
6 Platform for Research and Analysis in Environmental Sciences (EDST-PRASE), Beirut P.O. Box 6573/14, Lebanon
7 Department of Biological and Chemical Sciences, School of Arts and Sciences, Lebanese International University, Beirut, Lebanon- PO Box: 146404 Mazraa, Lebanon
doi
10.48309/chemm.2025.523370.1956چکیده
Remdesivir, a broad-spectrum antiviral agent, has gained widespread attention for treating SARS-CoV-2; however, its therapeutic efficacy remains limited by poor solubility, rapid metabolism, and inefficient delivery. In this study, we employed all-atom molecular dynamics simulations to investigate the adsorption behavior of Remdesivir on a silica-functionalized graphene oxide (GO–SiNP) nanocarrier. The GO–SiNP hybrid integrates the π-conjugated structure and high surface area of graphene oxide with the surface reactivity and aqueous dispersibility of silica nanoparticles, enabling robust non-covalent interactions through hydrogen bonding, π–π stacking, and van der Waals forces. Radial distribution functions indicated a higher molecular affinity in the GO–SiNP system, with a closer peak proximity (~0.28 nm), while RMSD and RMSF analyses showed enhanced conformational stability of the drug at the hybrid interface. Key descriptors of drug–carrier interactions—including RDFs, structural fluctuation, and MM-PBSA binding energies—demonstrated increased interaction persistence and enthalpic prevalence in the GO–SiNP system. The selection of an appropriate partial charge model was found to critically affect the electrostatic interactions and overall binding behavior of Remdesivir with the GO–SiNP nanocarrier. A comparative evaluation with other nanocarriers revealed that GO–SiNP exhibited the most favorable binding affinity (−41.8 ± 1.6 kcal/mol), highlighting its superior potential for antiviral drug delivery. These findings emphasize the importance of tailored surface functionalization in optimizing drug–carrier interactions and support the GO–SiNP hybrid as a promising platform for stable, high-affinity antiviral delivery systems.
