Mechanism of reaction of RNA-dependent RNA polymerase from SARS-CoV-2

Viral RNA Replication Suppression of SARS-CoV-2: Atomistic Insights into Inhibition Mechanisms of RdRp Machinery by ddhCTP

The nonstructural protein 12, known as RNA-dependent RNA polymerase (RdRp), is essential for both replication and repair of the viral genome. The RdRp of SARS-CoV-2 has been used as a promising candidate for drug development since the inception of the COVID-19 spread. In this work, we performed an in silico investigation on the insertion of the naturally modified pyrimidine nucleobase ddhCTP into the SARS-CoV-2 RdRp active site, in a comparative analysis with the natural one (CTP). The modification in ddhCTP involves the removal of the 3'-hydroxyl group that prevents the addition of subsequent nucleotides into the nascent strand, acting as an RNA chain terminator inhibitor. Quantum mechanical investigations helped to shed light on the mechanistic source of RdRp activity on the selected nucleobases, and comprehensive all-atom simulations provided insights about the structural rearrangements occurring in the active-site region when inorganic pyrophosphate (PPi) is formed. Subsequently, the intricate pathways for the release of PPi, the catalytic product of RdRp, were investigated using Umbrella Sampling simulations. The results are in line with the available experimental data and contribute to a more comprehensive point of view on such an important viral enzyme.

Targeting RdRp of SARS-CoV-2 with De Novo Molecule Generation

Viruses are known for their extremely high mutation rates, allowing them to evade both the human immune system and many forms of standard medicine. Despite this, the RNA dependent RNA polymerase (RdRp) of the RNA viruses has been largely conserved, and any significant mutation of this protein is unlikely. The recent COVID-19 pandemic presents a need for therapeutics. We have designed a de novo drug design algorithm that generates strong binding ligands from scratch, based on only the structure of the target protein's receptor. In this paper, we applied our method to target SARS-CoV-2 RdRp and generated several de novo molecules. We then chose some drug molecules based on the structural similarity to some of our strongest binding de novo molecules. Subsequently, we showed, using rigorous all-atom explicit-water free energy calculations in near-microsecond time scales using state-of-the-art well-tempered metadynamics simulations, that some of our de novo generated ligands bind more strongly to RdRp than the recent FDA approved drug remdesivir in its active form, remdesivir triphosphate (RTP). We elucidated the binding mechanism for some of the top binders and compared it with RTP. We believe that this work will be useful both by presenting lead structures for RdRp inhibition and by delivering key insights into the residues of the protein potentially involved in the binding/unbinding of these small molecule drugs, leading to more targeted studies in the future.

The Functional Implications of Broad Spectrum Bioactive Compounds Targeting RNA-Dependent RNA Polymerase (RdRp) in the Context of the COVID-19 Pandemic

BACKGROUND

As long as COVID-19 endures, viral surface proteins will keep changing and new viral strains will emerge, rendering prior vaccines and treatments decreasingly effective. To provide durable targets for preventive and therapeutic agents, there is increasing interest in slowly mutating viral proteins, including non-surface proteins like RdRp.

METHODS

A scoping review of studies was conducted describing RdRp in the context of COVID-19 through MEDLINE/PubMed and EMBASE. An iterative approach was used with input from content experts and three independent reviewers, focused on studies related to either RdRp activity inhibition or RdRp mechanisms against SARS-CoV-2.

RESULTS

Of the 205 records screened, 43 studies were included in the review. Twenty-five evaluated RdRp activity inhibition, and eighteen described RdRp mechanisms of existing drugs or compounds against SARS-CoV-2. In silico experiments suggested that RdRp inhibitors developed for other RNA viruses may be effective in disrupting SARS-CoV-2 replication, indicating a possible reduction of disease progression from current and future variants. In vitro, in vivo, and human clinical trial studies were largely consistent with these findings.

CONCLUSIONS

Future risk mitigation and treatment strategies against forthcoming SARS-CoV-2 variants should consider targeting RdRp proteins instead of surface proteins.

SARS-CoV-2 RNA-Dependent RNA Polymerase Follows Asynchronous Translocation Pathway for Viral Transcription and Replication

Translocation is one essential step for the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) to exert viral replication and transcription. Although cryo-EM structures of SARS-CoV-2 RdRp are available, the molecular mechanisms of dynamic translocation remain elusive. Herein, we constructed a Markov state model based on extensive molecular dynamics simulations to elucidate the translocation dynamics of the SARS-CoV-2 RdRp. We identified two intermediates that pinpoint the rate-limiting step of translocation and characterize the asynchronous movement of the template-primer duplex. The 3'-terminal nucleotide in the primer strand lags behind due to the uneven distribution of protein-RNA interactions, while the translocation of the template strand is delayed by the hurdle residue K500. Even so, the two strands share the same "ratchet" to stabilize the polymerase in the post-translocation state, suggesting a Brownian-ratchet model. Overall, our study provides intriguing insights into SARS-CoV-2 replication and transcription, which would open a new avenue for drug discoveries.

Resolving a guanine-quadruplex structure in the SARS-CoV-2 genome through circular dichroism and multiscale molecular modeling

The genome of SARS-CoV-2 coronavirus is made up of a single-stranded RNA fragment that can assume a specific secondary structure, whose stability can influence the virus's ability to reproduce. Recent studies have identified putative guanine quadruplex sequences in SARS-CoV-2 genome fragments that are involved in coding for both structural and non-structural proteins. In this contribution, we focus on a specific G-rich sequence referred to as RG-2, which codes for the non-structural protein 10 (Nsp10) and assumes a guanine-quadruplex (G4) arrangement. We provide the secondary structure of RG-2 G4 at atomistic resolution by molecular modeling and simulation, validated by the superposition of experimental and calculated electronic circular dichroism spectra. Through both experimental and simulation approaches, we have demonstrated that pyridostatin (PDS), a widely recognized G4 binder, can bind to and stabilize RG-2 G4 more strongly than RG-1, another G4 forming sequence that was previously proposed as a potential target for antiviral drug candidates. Overall, this study highlights RG-2 as a valuable target to inhibit the translation and replication of SARS-CoV-2, paving the way towards original therapeutic approaches against emerging RNA viruses.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

Current understanding of nucleoside analogs inhibiting the SARS-CoV-2 RNA-dependent RNA polymerase

Since the outbreak of the COVID-19 pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA-dependent RNA polymerase (RdRp) has become a main target for antiviral therapeutics due to its essential role in viral replication and transcription. Thus, nucleoside analogs structurally resemble the natural RdRp substrate and hold great potential as inhibitors. Until now, extensive experimental investigations have been performed to explore nucleoside analogs to inhibit the RdRp, and concerted efforts have been made to elucidate the underlying molecular mechanisms further. This review begins by discussing the nucleoside analogs that have demonstrated inhibition in the experiments. Second, we examine the current understanding of the molecular mechanisms underlying the action of nucleoside analogs on the SARS-CoV-2 RdRp. Recent findings in structural biology and computational research are presented through the classification of inhibitory mechanisms. This review summarizes previous experimental findings and mechanistic investigations of nucleoside analogs inhibiting SARS-CoV-2 RdRp. It would guide the rational design of antiviral medications and research into viral transcriptional mechanisms.

Reflections on the Origin and Early Evolution of the Genetic Code

Examination of the genetic code (GeCo) reveals that amino acids coded by (A/U) codons display a large functional spectrum and bind RNA whereas, except for Arg, those coded by (G/C) codons do not. From a stereochemical viewpoint, the clear preference for (A/U)-rich codons to be located at the GeCo half blocks suggests they were specifically determined. Conversely, the overall lower affinity of cognate amino acids for their (G/C)-rich anticodons points to their late arrival to the GeCo. It is proposed that i) initially the code was composed of the eight (A/U) codons; ii) these codons were duplicated when G/C nucleotides were added to their wobble positions, and three new codons with G/C in their first position were incorporated; and iii) a combination of A/U and G/C nucleotides progressively generated the remaining codons.

How SARS-CoV-2 Alters the Regulation of Gene Expression in Infected Cells

Nonstructural accessory proteins in viruses play a key role in hijacking the basic cellular mechanisms, which is essential to promote the virus survival and evasion of the immune system. The immonuglobulin-like open reading frame 8 (ORF8) protein expressed by SARS-CoV-2 accumulates in the nucleus and may influence the regulation of the gene expression in infected cells. In this contribution, by using microsecond time-scale all-atom molecular dynamics simulations, we unravel the structural bases behind the epigenetic action of ORF8. In particular, we highlight how the protein is able to form stable aggregates with DNA through a histone tail-like motif, and how this interaction is influenced by post-translational modifications, such as acetylation and methylation, which are known epigenetic markers in histones. Our work not only clarifies the molecular mechanisms behind the perturbation of the epigenetic regulation caused by the viral infection but also offers an unusual perspective which may foster the development of original antivirals.

1'-Cyano Intermediate Enables Rapid and Stereoretentive Access to 1'-Modified Remdesivir Nucleosides

While biochemical, structural, and computational studies have shown the importance of remdesivir's C1'-substituent in its perturbation of SARS-CoV-2 RdRp action, we recognized the paucity of methods to stereoselectively install substituents at this position as an obstacle to rigorous explorations of SAR and mechanism. We report the utilization of an anomerically pure 1'-cyano intermediate as an entry point to a chemically diverse set of substitutions, allowing for 1'diversification while obviating the need for the tedious separation of anomeric mixtures.

On the Recognition of Natural Substrate CTP and Endogenous Inhibitor ddhCTP of SARS-CoV-2 RNA-Dependent RNA Polymerase: A Molecular Dynamics Study

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 outbreak that is affecting the entire planet. As the pandemic is still spreading worldwide, with multiple mutations of the virus, it is of interest and of help to employ computational methods for identifying potential inhibitors of the enzymes responsible for viral replication. Attractive antiviral nucleotide analogue RNA-dependent RNA polymerase (RdRp) chain terminator inhibitors are investigated with this purpose. This study, based on molecular dynamics (MD) simulations, addresses the important aspects of the incorporation of an endogenously synthesized nucleoside triphosphate, ddhCTP, in comparison with the natural nucleobase cytidine triphosphate (CTP) in RdRp. The ddhCTP species is the product of the viperin antiviral protein as part of the innate immune response. The absence of the ribose 3'-OH in ddhCTP could have important implications in its inhibitory mechanism of RdRp. We built an in silico model of the RNA strand embedded in RdRp using experimental methods, starting from the cryo-electron microscopy structure and exploiting the information obtained by spectrometry on the RNA sequence. We determined that the model was stable during the MD simulation time. The obtained results provide deeper insights into the incorporation of nucleoside triphosphates, whose molecular mechanism by the RdRp active site still remains elusive.

Modeling the Enzymatic Mechanism of the SARS-CoV-2 RNA-Dependent RNA Polymerase by DFT/MM-MD: An Unusual Active Site Leading to High Replication Rates

Viral infection relies on the hijacking of cellular machineries to enforce the reproduction of the infecting virus and its subsequent diffusion. In this context, the replication of the viral genome is a key step performed by specific enzymes, i.e., polymerases. The replication of SARS-CoV-2, the causative agent of the COVID-19 pandemics, is based on the duplication of its RNA genome, an action performed by the viral RNA-dependent RNA polymerase. In this contribution, by using highly demanding DFT/MM-MD computations coupled to 2D-umbrella sampling techniques, we have determined the chemical mechanisms leading to the inclusion of a nucleotide in the nascent viral RNA strand. These results highlight the high efficiency of the polymerase, which lowers the activation free energy to less than 10 kcal/mol. Furthermore, the SARS-CoV-2 polymerase active site is slightly different from those usually found in other similar enzymes, and in particular, it lacks the possibility to enforce a proton shuttle via a nearby histidine. Our simulations show that this absence is partially compensated by lysine whose proton assists the reaction, opening up an alternative, but highly efficient, reactive channel. Our results present the first mechanistic resolution of SARS-CoV-2 genome replication at the DFT/MM-MD level and shed light on its unusual enzymatic reactivity paving the way for the future rational design of antivirals targeting emerging RNA viruses.

State-of-the-Art Molecular Dynamics Simulation Studies of RNA-Dependent RNA Polymerase of SARS-CoV-2

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.

In Silico Study of Alkaloids: Neferine and Berbamine Potentially Inhibit the SARS-CoV-2 RNA-Dependent RNA Polymerase

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, has been a global concern. While there have been some vaccines and drugs, the rapid emergence of variants due to mutations has threatened public health. As the de novo drug development process is expensive and time-consuming, repurposing existing antiviral drugs against SARS-CoV-2 is an alternative and promising approach to mitigate the current situation. Several studies have indicated that some natural products exhibit inhibitory activities against SARS-CoV-2. This study is aimed at analyzing the potential of natural alkaloids, using various computational tools, as drug candidates against SARS-CoV-2. The molecular docking analysis predicted that naturally occurring alkaloids can bind with RNA-dependent RNA-polymerase (RdRP). The QSAR analysis was conducted by using the way2drug/PASS online web resource, and the pharmacokinetics and toxicity properties of these alkaloids were predicted using pkCSM, SwissADME, and ProTox-II webserver. Among the different alkaloids studied, neferine and berbamine were repurposed as potential drug candidates based on their binding affinity and interactions with RdRP. Further, molecular dynamics simulation of 90 ns revealed the conformational stability of the neferine-RdRP complex.

Pre-exascale HPC approaches for molecular dynamics simulations. Covid-19 research: A use case

Exascale computing has been a dream for ages and is close to becoming a reality that will impact how molecular simulations are being performed, as well as the quantity and quality of the information derived for them. We review how the biomolecular simulations field is anticipating these new architectures, making emphasis on recent work from groups in the BioExcel Center of Excellence for High Performance Computing. We exemplified the power of these simulation strategies with the work done by the HPC simulation community to fight Covid-19 pandemics. This article is categorized under:Data Science > Computer Algorithms and ProgrammingData Science > Databases and Expert SystemsMolecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods.