Dissecting nucleotide selectivity in viral RNA polymerases

Trapping a non-cognate nucleotide upon initial binding for replication fidelity control in SARS-CoV-2 RNA dependent RNA polymerase

The RNA dependent RNA polymerase (RdRp) in SARS-CoV-2 is a highly conserved enzyme responsible for viral genome replication/transcription. To understand how the viral RdRp achieves fidelity control during such processes, here we computationally investigate the natural non-cognate vs. cognate nucleotide addition and selectivity during viral RdRp elongation. We focus on the nucleotide substrate initial binding (RdRp active site open) to the prechemical insertion (active site closed) of the RdRp. The current studies were first carried out using microsecond ensemble equilibrium all-atom molecular dynamics (MD) simulations. Due to the slow conformational changes (from open to closed) during nucleotide insertion and selection, enhanced or umbrella sampling methods have been further employed to calculate the free energy profiles of the nucleotide insertion. Our studies find notable stability of noncognate dATP and GTP upon initial binding in the active-site open state. The results indicate that while natural cognate ATP and Remdesivir drug analogue (RDV-TP) are biased toward stabilization in the closed state to facilitate insertion, the natural non-cognate dATP and GTP can be well trapped in off-path initial binding configurations and prevented from insertion so that to be further rejected. The current work thus presents the intrinsic nucleotide selectivity of SARS-CoV-2 RdRp for natural substrate fidelity control, which should be considered in antiviral drug design.

Kill or corrupt: Mechanisms of action and drug-resistance of nucleotide analogues against SARS-CoV-2

Nucleoside/tide analogues (NAs) have long been used in the fight against viral diseases, and now present a promising option for the treatment of COVID-19. Once activated to the 5'-triphosphate state, NAs act by targeting the viral RNA-dependent RNA-polymerase for incorporation into the viral RNA genome. Incorporated analogues can either 'kill' (terminate) synthesis, or 'corrupt' (genetically or chemically) the RNA. Against coronaviruses, the use of NAs has been further complicated by the presence of a virally encoded exonuclease domain (nsp14) with proofreading and repair capacities. Here, we describe the mechanism of action of four promising anti-COVID-19 NAs; remdesivir, molnupiravir, favipiravir and bemnifosbuvir. Their distinct mechanisms of action best exemplify the concept of 'killers' and 'corruptors'. We review available data regarding their ability to be incorporated and excised, and discuss the specific structural features that dictate their overall potency, toxicity, and mutagenic potential. This should guide the synthesis of novel analogues, lend insight into the potential for resistance mutations, and provide a rational basis for upcoming combinations therapies.

Atypical Mutational Spectrum of SARS-CoV-2 Replicating in the Presence of Ribavirin

We report that ribavirin exerts an inhibitory and mutagenic activity on SARS-CoV-2-infecting Vero cells, with a therapeutic index higher than 10. Deep sequencing analysis of the mutant spectrum of SARS-CoV-2 replicating in the absence or presence of ribavirin indicated an increase in the number of mutations, but not in deletions, and modification of diversity indices, expected from a mutagenic activity. Notably, the major mutation types enhanced by replication in the presence of ribavirin were A→G and U→C transitions, a pattern which is opposite to the dominance of G→A and C→U transitions previously described for most RNA viruses. Implications of the inhibitory activity of ribavirin, and the atypical mutational bias produced on SARS-CoV-2, for the search for synergistic anti-COVID-19 lethal mutagen combinations are discussed.

Unveiling the "Template-Dependent" Inhibition on the Viral Transcription of SARS-CoV-2

Remdesivir is one nucleotide analogue prodrug capable to terminate RNA synthesis in SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) by two distinct mechanisms. Although the "delayed chain termination" mechanism has been extensively investigated, the "template-dependent" inhibitory mechanism remains elusive. In this study, we have demonstrated that remdesivir embedded in the template strand seldom directly disrupted the complementary NTP incorporation at the active site. Instead, the translocation of remdesivir from the +2 to the +1 site was hindered due to the steric clash with V557. Moreover, we have elucidated the molecular mechanism characterizing the drug resistance upon V557L mutation. Overall, our studies have provided valuable insight into the "template-dependent" inhibitory mechanism exerted by remdesivir on SARS-CoV-2 RdRp and paved venues for an alternative antiviral strategy for the COVID-19 pandemic. As the "template-dependent" inhibition occurs across diverse viral RdRps, our findings may also shed light on a common acting mechanism of inhibitors.

Elucidation of the conformational dynamics and assembly of Argonaute−RNA complexes by distinct yet coordinated actions of the supplementary microRNA

Argonaute (AGO) proteins, the core of RNA-induced silencing complex, are guided by microRNAs (miRNAs) to recognize target RNA for repression. The miRNA–target RNA recognition forms initially through pairing at the seed region while the additional supplementary pairing can enhance target recognition and compensate for seed mismatch. The extension of miRNA lengths can strengthen the target affinity when pairing both in the seed and supplementary regions. However, the mechanism underlying the effect of the supplementary pairing on the conformational dynamics and the assembly of AGO–RNA complex remains poorly understood. To address this, we performed large-scale molecular dynamics simulations of AGO–RNA complexes with different pairing patterns and miRNA lengths. The results reveal that the additional supplementary pairing can not only strengthen the interaction between miRNA and target RNA, but also induce the increased plasticity of the PAZ domain and enhance the domain connectivity among the PAZ, PIWI, N domains of the AGO protein. The strong community network between these domains tightens the mouth of the supplementary chamber of AGO protein, which prevents the escape of target RNA from the complex and shields it from solvent water attack. Importantly, the inner stronger matching pairs between the miRNA and target RNA can compensate for weaker mismatches at the edge of supplementary region. These findings provide guidance for the design of miRNA mimics and anti-miRNAs for both clinical and experimental use and open the way for further engineering of AGO proteins as a new tool in the field of gene regulation.

Identification of novel SARS-CoV-2 RNA dependent RNA polymerase (RdRp) inhibitors: From in silico screening to experimentally validated inhibitory activity

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 2019 has posed a serious threat to global health and the economy for over two years, prompting the need for development of antiviral inhibitors. Due to its vital role in viral replication, RNA-dependent RNA polymerase (RdRp) is a promising therapeutic target. Herein, we analyzed amino acid sequence conservation of RdRp across coronaviruses. The conserved amino acids at the catalytic binding site served as the ligand-contacting residues for in silico screening to elucidate possible resistant mutation. Molecular docking was employed to screen inhibitors of SARS-CoV-2 from the ZINC ChemDiv database. The top-ranked compounds selected from GOLD docking were further investigated for binding modes at the conserved residues of RdRp, and ten compounds were selected for experimental validation. Of which, three compounds exhibited promising antiviral activity. The most promising candidate showed a half-maximal effective concentration (EC50) of 5.04 µM. Molecular dynamics simulations, binding free-energy calculation and hydrogen bond analysis were performed to elucidate the critical interactions providing a foundation for developing lead compounds effective against SARS-CoV-2.

Co-crystallization and structure determination: An effective direction for anti-SARS-CoV-2 drug discovery

Safer and more-effective drugs are urgently needed to counter infections with the highly pathogenic SARS-CoV-2, cause of the COVID-19 pandemic. Identification of efficient inhibitors to treat and prevent SARS-CoV-2 infection is a predominant focus. Encouragingly, using X-ray crystal structures of therapeutically relevant drug targets (PLpro, Mpro, RdRp, and S glycoprotein) offers a valuable direction for anti-SARS-CoV-2 drug discovery and lead optimization through direct visualization of interactions. Computational analyses based primarily on MMPBSA calculations have also been proposed for assessing the binding stability of biomolecular structures involving the ligand and receptor. In this study, we focused on state-of-the-art X-ray co-crystal structures of the abovementioned targets complexed with newly identified small-molecule inhibitors (natural products, FDA-approved drugs, candidate drugs, and their analogues) with the assistance of computational analyses to support the precision design and screening of anti-SARS-CoV-2 drugs.