1
|
Vishwakarma K, Ravi S, Mittal S. Ab initio Modeling of Hydrogen Bonding of Remdesivir and Adenosine with Uridine. Chemphyschem 2024; 25:e202300552. [PMID: 37983746 DOI: 10.1002/cphc.202300552] [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: 08/03/2023] [Revised: 11/18/2023] [Accepted: 11/19/2023] [Indexed: 11/22/2023]
Abstract
Remdesivir (RDV) emerged as an effective drug against the SARS-CoV-2 virus pandemic. One of the crucial steps in the mechanism of action of RDV is its incorporation into the growing RNA strand. RDV, an adenosine analogue, forms Watson-Crick (WC) type hydrogen bonds with uridine in the complementary strand and the strength of this interaction will control efficacy of RDV. While there is a plethora of structural and energetic information available about WC H-bonds in natural base pairs, the interaction of RDV with uridine has not been studied yet at the atomic level. In this article, we aim to bridge this gap, to understand RDV and its hydrogen bonding interactions, by employing density functional theory (DFT) at the M06-2X/cc-pVDZ level. The interaction energy, QTAIM analysis, NBO and SAPT2 are performed for RDV, adenosine, and their complex with uridine to gain insights into the nature of hydrogen bonding. The computations show that RDV has similar geometry, energetic, molecular orbitals, and aromaticity as adenosine, suggesting that RDV is an effective adenosine analogue. The important geometrical parameters, such as bond distances and red-shift in the stretching vibrational modes of adenosine, RDV and uridine identify two WC-type H-bonds. The relative strength of these two H-bonds is computed using QTAIM parameters and the computed hydrogen bond energy. Finally, the SAPT2 study is performed at the minima and at non-equilibrium base pair distances to understand the dominant intermolecular physical force. This study, based on a thorough analysis of a variety of computations, suggests that both adenosine and RDV have similar structure, energetic, and hydrogen bonding behaviour.
Collapse
Affiliation(s)
- Kamini Vishwakarma
- School of Advance Science and Languages, VIT Bhopal University, Kothrikalan, Sehore, Madhya, Pradesh, 466114, India
| | - Satyam Ravi
- School of Advance Science and Languages, VIT Bhopal University, Kothrikalan, Sehore, Madhya, Pradesh, 466114, India
| | - Sumit Mittal
- School of Advance Science and Languages, VIT Bhopal University, Kothrikalan, Sehore, Madhya, Pradesh, 466114, India
| |
Collapse
|
2
|
Kulkarni S, Singh Y, Biharee A, Bhatia N, Monga V, Thareja S. Molecular docking, 3D-QSAR and simulation studies for identifying pharmacophoric features of indole derivatives as 17β-hydroxysteroid dehydrogenase type 5 (17β-HSD5) inhibitors. J Biomol Struct Dyn 2023; 41:12668-12685. [PMID: 36744535 DOI: 10.1080/07391102.2023.2175265] [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: 12/08/2022] [Accepted: 01/05/2023] [Indexed: 02/07/2023]
Abstract
Excess of androgens leads to various diseases such as Poly-Cystic Ovarian Syndrome, Prostate Cancer, Hirsutism, Obesity and Acne. 17β-Hydroxysteroid Dehydrogenase type 5 (17β-HSD5) converts androstenedione into testosterone peripherally, thereby significantly contributing to the development of these diseases. Indole-bearing scaffolds are reported as potential 17β-HSD5 inhibitors for the manifestation of diseases arising due to androgen excess. In the present work, we have extensively performed a combination of molecular docking, Gaussian field-based 3D-QSAR, Pharmacophore mapping and MD-simulation studies (100 ns) to identify the pharmacophoric features of indole-based compounds as potent 17β-HSD5 inhibitors. Molecular simulation studies of the most potent compound in the binding pocket of enzyme revealed that the compound 11 was stable in the binding pocket and showed good binding affinity through interactions with various residues of active site pocket. The Molecular mechanics Generalized Born surface area continuum solvation (MM/GBSA) and Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations revealed that the compound 11 possessed a free binding energy of -36.36 kcal/mol and -7.00 kcal/mol, respectively, which was better as compared to reference compound Desmethyl indomethacin (DES). The developed pharmacophore will be helpful to design novel indole-based molecules as potent 17β-HSD5 inhibitors for the treatment of various androgenic disorders.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Swanand Kulkarni
- Department of Pharmaceutical Sciences and Natural Products, Central University of Punjab, Bathinda, Punjab, India
| | - Yogesh Singh
- Department of Pharmaceutical Sciences and Natural Products, Central University of Punjab, Bathinda, Punjab, India
| | - Avadh Biharee
- Department of Pharmaceutical Sciences and Natural Products, Central University of Punjab, Bathinda, Punjab, India
| | - Neha Bhatia
- Department of Pharmaceutical Sciences and Natural Products, Central University of Punjab, Bathinda, Punjab, India
| | - Vikramdeep Monga
- Department of Pharmaceutical Sciences and Natural Products, Central University of Punjab, Bathinda, Punjab, India
| | - Suresh Thareja
- Department of Pharmaceutical Sciences and Natural Products, Central University of Punjab, Bathinda, Punjab, India
| |
Collapse
|
3
|
Tanimoto S, Itoh SG, Okumura H. State-of-the-Art Molecular Dynamics Simulation Studies of RNA-Dependent RNA Polymerase of SARS-CoV-2. Int J Mol Sci 2022; 23:ijms231810358. [PMID: 36142270 PMCID: PMC9499461 DOI: 10.3390/ijms231810358] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Revised: 08/26/2022] [Accepted: 08/31/2022] [Indexed: 01/18/2023] Open
Abstract
Molecular dynamics (MD) simulations are powerful theoretical methods that can reveal biomolecular properties, such as structure, fluctuations, and ligand binding, at the level of atomic detail. In this review article, recent MD simulation studies on these biomolecular properties of the RNA-dependent RNA polymerase (RdRp), which is a multidomain protein, of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are presented. Although the tertiary structures of RdRps in SARS-CoV-2 and SARS-CoV are almost identical, the RNA synthesis activity of RdRp of SARS-CoV is higher than SARS-CoV-2. Recent MD simulations observed a difference in the dynamic properties of the two RdRps, which may cause activity differences. RdRp is also a drug target for Coronavirus disease 2019 (COVID-19). Nucleotide analogs, such as remdesivir and favipiravir, are considered to be taken up by RdRp and inhibit RNA replication. Recent MD simulations revealed the recognition mechanism of RdRp for these drug molecules and adenosine triphosphate (ATP). The ligand-recognition ability of RdRp decreases in the order of remdesivir, favipiravir, and ATP. As a typical recognition process, it was found that several lysine residues of RdRp transfer these ligand molecules to the binding site such as a “bucket brigade.” This finding will contribute to understanding the mechanism of the efficient ligand recognition by RdRp. In addition, various simulation studies on the complexes of SARS-CoV-2 RdRp with several nucleotide analogs are reviewed, and the molecular mechanisms by which these compounds inhibit the function of RdRp are discussed. The simulation studies presented in this review will provide useful insights into how nucleotide analogs are recognized by RdRp and inhibit the RNA replication.
Collapse
Affiliation(s)
- Shoichi Tanimoto
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki 444-8787, Aichi, Japan
- Exploratory Research Center on Life and Living Systems, National Institutes of Natural Sciences, Okazaki 444-8787, Aichi, Japan
| | - Satoru G. Itoh
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki 444-8787, Aichi, Japan
- Exploratory Research Center on Life and Living Systems, National Institutes of Natural Sciences, Okazaki 444-8787, Aichi, Japan
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki 444-8787, Aichi, Japan
| | - Hisashi Okumura
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki 444-8787, Aichi, Japan
- Exploratory Research Center on Life and Living Systems, National Institutes of Natural Sciences, Okazaki 444-8787, Aichi, Japan
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki 444-8787, Aichi, Japan
- Correspondence:
| |
Collapse
|
4
|
Remdesivir MD Simulations Suggest a More Favourable Binding to SARS-CoV-2 RNA Dependent RNA Polymerase Mutant P323L Than Wild-Type. Biomolecules 2021; 11:biom11070919. [PMID: 34206274 PMCID: PMC8301449 DOI: 10.3390/biom11070919] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2021] [Revised: 05/18/2021] [Accepted: 05/21/2021] [Indexed: 02/02/2023] Open
Abstract
SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) protein is the target for the antiviral drug Remdesivir (RDV). With RDV clinical trials on COVID-19 patients showing a reduced hospitalisation time. During the spread of the virus, the RdRp has developed several mutations, with the most frequent being A97V and P323L. The current study sought to investigate whether A97V and P323L mutations influence the binding of RDV to the RdRp of SARS-CoV-2 compared to wild-type (WT). The interaction of RDV with WT-, A97V-, and P323L-RdRp were measured using molecular dynamic (MD) simulations, and the free binding energies were extracted. Results showed that RDV that bound to WT- and A97V-RdRp had a similar dynamic motion and internal residue fluctuations, whereas RDV interaction with P323L-RdRp exhibited a tighter molecular conformation, with a high internal motion near the active site. This was further corroborated with RDV showing a higher binding affinity to P323L-RdRp (-24.1 kcal/mol) in comparison to WT-RdRp (-17.3 kcal/mol). This study provides insight into the potential significance of administering RDV to patients carrying the SARS-CoV-2 P323L-RdRp mutation, which may have a more favourable chance of alleviating the SARS-CoV-2 illness in comparison to WT-RdRp carriers, thereby suggesting further scientific consensus for the usage of Remdesivir as clinical candidate against COVID-19.
Collapse
|
5
|
Allosteric Activation of SARS-CoV-2 RNA-Dependent RNA Polymerase by Remdesivir Triphosphate and Other Phosphorylated Nucleotides. mBio 2021; 12:e0142321. [PMID: 34154407 PMCID: PMC8262916 DOI: 10.1128/mbio.01423-21] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The catalytic subunit of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA-dependent RNA polymerase (RdRp) Nsp12 has a unique nidovirus RdRp-associated nucleotidyltransferase (NiRAN) domain that transfers nucleoside monophosphates to the Nsp9 protein and the nascent RNA. The NiRAN and RdRp modules form a dynamic interface distant from their catalytic sites, and both activities are essential for viral replication. We report that codon-optimized (for the pause-free translation in bacterial cells) Nsp12 exists in an inactive state in which NiRAN-RdRp interactions are broken, whereas translation by slow ribosomes and incubation with accessory Nsp7/8 subunits or nucleoside triphosphates (NTPs) partially rescue RdRp activity. Our data show that adenosine and remdesivir triphosphates promote the synthesis of A-less RNAs, as does ppGpp, while amino acid substitutions at the NiRAN-RdRp interface augment activation, suggesting that ligand binding to the NiRAN catalytic site modulates RdRp activity. The existence of allosterically linked nucleotidyl transferase sites that utilize the same substrates has important implications for understanding the mechanism of SARS-CoV-2 replication and the design of its inhibitors.
Collapse
|
6
|
Wang B, Svetlov V, Wolf YI, Koonin EV, Nudler E, Artsimovitch I. Allosteric activation of SARS-CoV-2 RdRp by remdesivir triphosphate and other phosphorylated nucleotides. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2021:2021.01.24.428004. [PMID: 33948598 PMCID: PMC8095223 DOI: 10.1101/2021.01.24.428004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The catalytic subunit of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp), Nsp12, has a unique NiRAN domain that transfers nucleoside monophosphates to the Nsp9 protein. The NiRAN and RdRp modules form a dynamic interface distant from their catalytic sites and both activities are essential for viral replication. We report that codon-optimized (for the pause-free translation) Nsp12 exists in inactive state in which NiRAN/RdRp interactions are broken, whereas translation by slow ribosomes and incubation with accessory Nsp7/8 subunits or NTPs partially rescue RdRp activity. Our data show that adenosine and remdesivir triphosphates promote synthesis of A-less RNAs, as does ppGpp, while amino acid substitutions at the NiRAN/RdRp interface augment activation, suggesting that ligand binding to the NiRAN catalytic site modulates RdRp activity. The existence of allosterically-linked nucleotidyl transferase sites that utilize the same substrates has important implications for understanding the mechanism of SARS-CoV-2 replication and design of its inhibitors. HIGHLIGHTS Codon-optimization of Nsp12 triggers misfolding and activity lossSlow translation, accessory Nsp7 and Nsp8 subunits, and NTPs rescue Nsp12Non-substrate nucleotides activate RNA chain synthesis, likely via NiRAN domainCrosstalk between two Nsp12 active sites that bind the same ligands.
Collapse
Affiliation(s)
- Bing Wang
- Department of Microbiology, The Ohio State University, Columbus, OH 43210, USA
- The Center for RNA Biology, The Ohio State University, Columbus, OH 43210, USA
| | - Vladimir Svetlov
- Department of Biochemistry and Molecular Pharmacology, New York University School of Medicine, New York, NY 10016, USA
| | - Yuri I Wolf
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, 20894, USA
| | - Eugene V Koonin
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, 20894, USA
| | - Evgeny Nudler
- Department of Biochemistry and Molecular Pharmacology, New York University School of Medicine, New York, NY 10016, USA
- Howard Hughes Medical Institute, New York University School of Medicine, New York, NY 10016, USA
| | - Irina Artsimovitch
- Department of Microbiology, The Ohio State University, Columbus, OH 43210, USA
- The Center for RNA Biology, The Ohio State University, Columbus, OH 43210, USA
| |
Collapse
|