skip to main content

Explore Research Products in the NSF-PAR

Title: Probing remdesivir nucleotide analogue insertion to SARS-CoV-2 RNA dependent RNA polymerase in viral replication
Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. In this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (∼1.5 kcal mol −1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick base pairing with the template nt, and inserts into the active site with a medium low free energy barrier (∼2.6 kcal mol −1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.  more » « less
Award ID(s):
2028935 1763272
NSF-PAR ID:
10321349
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Molecular Systems Design & Engineering
Volume:
6
Issue:
11
ISSN:
2058-9689
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Scientific Reports)
    Abstract

    RNA dependent RNA polymerase (RdRp), is an essential in the RNA replication within the life cycle of the severely acute respiratory coronavirus-2 (SARS-CoV-2), causing the deadly respiratory induced sickness COVID-19. Remdesivir is a prodrug that has seen some success in inhibiting this enzyme, however there is still the pressing need for effective alternatives. In this study, we present the discovery of four non-nucleoside small molecules that bind favorably to SARS-CoV-2 RdRp over the active form of the popular drug remdesivir (RTP) and adenosine triphosphate (ATP) by utilizing high-throughput virtual screening (HTVS) against the vast ZINC compound database coupled with extensive molecular dynamics (MD) simulations. After post-trajectory analysis, we found that the simulations of complexes containing both ATP and RTP remained stable for the duration of their trajectories. Additionally, it was revealed that the phosphate tail of RTP was stabilized by both the positive amino acid pocket and magnesium ions near the entry channel of RdRp which includes residues K551, R553, R555 and K621. It was also found that residues D623, D760, and N691 further stabilized the ribose portion of RTP with U10 on the template RNA strand forming hydrogen pairs with the adenosine motif. Using these models of RdRp, we employed them to screen the ZINC database of ~ 17 million molecules. Using docking and drug properties scoring, we narrowed down our selection to fourteen candidates. These were subjected to 200 ns simulations each underwent free energy calculations. We identified four hit compounds from the ZINC database that have similar binding poses to RTP while possessing lower overall binding free energies, with ZINC097971592 having a binding free energy two times lower than RTP.

     
    more » « less
  2. ; ; ; ; ( , International Journal of Molecular Sciences)
    A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a −1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2. 
    more » « less
  3. ; ; ; ( , Nature Communications)
    Abstract

    Many viruses utilize trimeric spikes to gain entry into host cells. However, without in situ structures of these trimeric spikes, a full understanding of this dynamic and essential process of viral infections is not possible. Here we present four in situ and one isolated cryoEM structures of the trimeric spike of the cytoplasmic polyhedrosis virus, a member of the non-envelopedReoviridaefamily and a virus historically used as a model in the discoveries of RNA transcription and capping. These structures adopt two drastically different conformations, closed spike and opened spike, which respectively represent the penetration-inactive and penetration-active states. Each spike monomer has four domains: N-terminal, body, claw, and C-terminal. From closed to opened state, the RGD motif-containing C-terminal domain is freed to bind integrins, and the claw domain rotates to expose and project its membrane insertion loops into the cellular membrane. Comparison between turret vertices before and after detachment of the trimeric spike shows that thetrimericspike anchors its N-terminal domain in the iris of thepentamericRNA-capping turret. Sensing of cytosolic S-adenosylmethionine (SAM) and adenosine triphosphate (ATP) by the turret triggers a cascade of events: opening of the iris, detachment of the spike, and initiation of endogenous transcription.

     
    more » « less
  4. ; ; ; ; ; ; ( , Chemistry & Biodiversity)
    Abstract

    As COVID‐19 infection caused severe public health concerns recently, the development of novel antivirals has become the need of the hour. Main protease (Mpro) has been an attractive target for antiviral drugs since it plays a vital role in polyprotein processing and virus maturation. Herein we report the discovery of a novel class of inhibitors against the SARS‐CoV‐2, bearing histidineα‐nitrile motif embedded on a simple dipeptide framework.In‐vitroandin‐silicostudies revealed that the histidineα‐nitrile motif envisioned to target the Mprocontributes to the inhibitory activity. Among a series of dipeptides synthesized featuring this novel structural motif, some dipeptides displayed strong viral reduction (EC50=0.48 μM) with a high selectivity index, SI>454.54. These compounds also exhibit strong binding energies in the range of −28.7 to −34.2 Kcal/mol. The simple dipeptide structural framework, amenable to quick structural variations, coupled with ease of synthesis from readily available commercial starting materials are the major attractive features of this novel class of SARS‐CoV‐2 inhibitors. The histidineα‐nitrile dipeptides raise the hope of discovering potent drug candidates based on this motif to fight the dreaded SARS‐CoV‐2.

     
    more » « less
  5. ; ; ; ; ; ; ( , Nucleic Acids Research)
    Abstract

    The enteroviral 2C protein is a therapeutic target, but the absence of a mechanistic framework for this enzyme limits our understanding of inhibitor mechanisms. Here, we use poliovirus 2C and a derivative thereof to elucidate the first biochemical mechanism for this enzyme and confirm the applicability of this mechanism to other members of the enterovirus genus. Our biochemical data are consistent with a dimer forming in solution, binding to RNA, which stimulates ATPase activity by increasing the rate of hydrolysis without impacting affinity for ATP substantially. Both RNA and DNA bind to the same or overlapping site on 2C, driven by the phosphodiester backbone, but only RNA stimulates ATP hydrolysis. We propose that RNA binds to 2C driven by the backbone, with reorientation of the ribose hydroxyls occurring in a second step to form the catalytically competent state. 2C also uses a two-step mechanism for binding to ATP. Initial binding is driven by the α and β phosphates of ATP. In the second step, the adenine base and other substituents of ATP are used to organize the active site for catalysis. These studies provide the first biochemical description of determinants driving specificity and catalytic efficiency of a picornaviral 2C ATPase.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email