skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Virtual Screening of the Indonesian Medicinal Plant and Zinc Databases for Potential Inhibitors of the RNA-Dependent RNA Polymerase (RdRp) of 2019 Novel Coronavirus
The novel coronavirus disease 19 (Covid-19) which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a pandemic across the world, which necessitate the need for the antiviral drug discovery. One of the potential protein targets for coronavirus treatment is RNA-dependent RNA polymerase. It is the key enzyme in the viral replication machinery, and it does not exist in human beings, therefore its targeting has been considered as a strategic approach. Here we describe the identification of potential hits from Indonesian Herbal and ZINC databases. The pharmacophore modeling was employed followed by molecular docking and dynamics simulation for 40 ns. 151 and 14480 hit molecules were retrieved from Indonesian herbal and ZINC databases, respectively. Three hits that were selected based on the structural analysis were stable during 40 ns, while binding energy prediction further implied that ZINC1529045114, ZINC169730811, and 9-Ribosyl-trans-zeatin had tighter binding affinities compared to Remdesivir. The ZINC169730811 had the strongest affinity toward RdRp compared to the other two hits including Remdesivir and its binding was corroborated by electrostatic, van der Waals, and nonpolar contribution for solvation energies. The present study offers three hits showing tighter binding to RdRp based on MM-PBSA binding energy prediction for further experimental verification.  more » « less
Award ID(s):
1904797
PAR ID:
10234026
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Indonesian Journal of Chemistry
Volume:
20
Issue:
6
ISSN:
1411-9420
Page Range / eLocation ID:
1430; 1440
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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. 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
  3. Hepatitis C Virus (HCV) is reported to develop GS-461203 resistance because of multiple mutations within the RNA-dependent RNA Polymerase (RdRp) of HCV. The lack of a high-resolution structure of these RdRp mutants in complex with GS-461203 hinders efforts to understand the drug resistance. Here we decipher the binding differences of GS-461203 in the wild type and mutated systems T179A or M289L of HCV RdRp Genotype 2a using homology modeling, molecular docking, and molecular dynamics simulation. Key residues responsible for GS-461203 binding were identified to be Arg48, Arg158, Asp318, Asp319, and Asp220, and that mutations T179A or M289L have caused conformational changes of GS-461203 in the RdRp active site. The affinities of GS-461203 were reduced in T179A system, but it became slightly stronger in the M289L system. Furthermore, we designed two new analogues of GS-461203 which encouragingly induced more stable interactions than GS-461203, and thus resulted in much better binding energies. This present study reveals how a single mutation, T179A or M289L, will modulate GS-461203 binding in HCV RdRp Genotype 2a, while introducing two novel analogues to overcome the drug resistance which may be good candidate for further experimental verification. 
    more » « less
  4. Abstract When it is in the template RNA, the naturally occurring m1A epitranscriptomic RNA modification was recently reported to be able to stop the RNA polymerization reaction catalyzed by the RNA dependent RNA polymerase (RdRp) of SARS-CoV-2. In this report, we report that m1A via its triphosphate form (m1ATP) can be incorporated into RNA by the same RdRp. These two findings point a new direction for antiviral drug development based on m1A for combatting COVID-19. More broadly, it is possible that the large pool of epigenetic RNA as well as DNA modifications could serve as a treasury for drug discovery aimed at combating various infectious and other diseases. 
    more » « less
  5. Daròs, José-Antonio (Ed.)
    Viroids, a fascinating group of plant pathogens, are subviral agents composed of single-stranded circular noncoding RNAs. It is well-known that nuclear-replicating viroids exploit host DNA-dependent RNA polymerase II (Pol II) activity for transcription from circular RNA genome to minus-strand intermediates, a classic example illustrating the intrinsic RNA-dependent RNA polymerase activity of Pol II. The mechanism for Pol II to accept single-stranded RNAs as templates remains poorly understood. Here, we reconstituted a robust in vitro transcription system and demonstrated that Pol II also accepts minus-strand viroid RNA template to generate plus-strand RNAs. Further, we purified the Pol II complex on RNA templates for nano-liquid chromatography-tandem mass spectrometry analysis and identified a remodeled Pol II missing Rpb4, Rpb5, Rpb6, Rpb7, and Rpb9, contrasting to the canonical 12-subunit Pol II or the 10-subunit Pol II core on DNA templates. Interestingly, the absence of Rpb9, which is responsible for Pol II fidelity, explains the higher mutation rate of viroids in comparison to cellular transcripts. This remodeled Pol II is active for transcription with the aid of TFIIIA-7ZF and appears not to require other canonical general transcription factors (such as TFIIA, TFIIB, TFIID, TFIIE, TFIIF, TFIIH, and TFIIS), suggesting a distinct mechanism/machinery for viroid RNA-templated transcription. Transcription elongation factors, such as FACT complex, PAF1 complex, and SPT6, were also absent in the reconstituted transcription complex. Further analyses of the critical zinc finger domains in TFIIIA-7ZF revealed the first three zinc finger domains pivotal for RNA template binding. Collectively, our data illustrated a distinct organization of Pol II complex on viroid RNA templates, providing new insights into viroid replication, the evolution of transcription machinery, as well as the mechanism of RNA-templated transcription. 
    more » « less