Comparative molecular docking and simulation analysis of molnupiravir and remdesivir with SARS-CoV-2 RNA dependent RNA polymerase (RdRp)

Treatment of SARS-CoV-2 targeting its RNA dependent RNA polymerase (RdRp) is of current interest. Remdesivir has been approved for the treatment of COVID-19 around the world. However, the drug has been linked with pharmacological limitations like adverse effects and reduced efficiency. Nevertheless, recent advancements have depicted molnupiravir as an effective therapeutic agent to target the SARS-CoV-2 RdRp. The drug has cleared both in vitro and in vivo screening. It is in phase-III clinical trial. Nonetheless, there are no data on themolecular binding interaction of molnupiravir with RdRp. Therefore, it is of interest to report the binding interaction of molnupiravir using molecular docking. It is also of interest to show its stability during interaction using molecular dynamics and binding free energy calculations along with drug likeliness and pharmacokinetic properties in comparison with remdesivir.     

Lycorine,a non-nucleoside RNA dependent RNA polymerase inhibitor,as potential treatment for emerging coronavirus infections

Background: Highly effective novel treatments need to be developed to suppress emerging coronavirus (CoV) infections such as COVID-19. The RNA dependent RNA polymerase (RdRp) among the viral proteins is known as an effective antiviral target. Lycorine is a phenanthridine Amaryllidaceae alkaloid isolated from the bulbs of Lycoris radiata (L'Hér.) Herb. and has various pharmacological bioactivities including antiviral function.Purpose: We investigated the direct-inhibiting action of lycorine on CoV's RdRp, as potential treatment for emerging CoV infections.Methods: We examined the inhibitory effect of lycorine on MERS-CoV, SARS-CoV, and SARS-CoV-2 infections, and then quantitatively measured the inhibitory effect of lycorine on MERS-CoV RdRp activity using a cell-based reporter assay. Finally, we performed the docking simulation with lycorine and SARS-CoV-2 RdRp.Results: Lycorine efficiently inhibited these CoVs with IC₅₀ values of 2.123 ± 0.053, 1.021 ± 0.025, and 0.878 ± 0.022 μM, respectively, comparable with anti-CoV effects of remdesivir. Lycorine directly inhibited MERS-CoV RdRp activity with an IC₅₀ of 1.406 ± 0.260 μM, compared with remdesivir's IC₅₀ value of 6.335 ± 0.731 μM. In addition, docking simulation showed that lycorine interacts with SARS-CoV-2 RdRp at the Asp623, Asn691, and Ser759 residues through hydrogen bonding, at which the binding affinities of lycorine (−6.2 kcal/mol) were higher than those of remdesivir (−4.7 kcal/mol).Conclusions: Lycorine is a potent non-nucleoside direct-acting antiviral against emerging coronavirus infections and acts by inhibiting viral RdRp activity; therefore, lycorine may be a candidate against the current COVID-19 pandemic.     

Characterization of the NiRAN domain from RNA-dependent RNA polymerase provides insights into a potential therapeutic target against SARS-CoV-2

Apart from the canonical fingers, palm and thumb domains, the RNA dependent RNA polymerases (RdRp) from the viral order Nidovirales possess two additional domains. Of these, the function of the Nidovirus RdRp associated nucleotidyl transferase domain (NiRAN) remains unanswered. The elucidation of the 3D structure of RdRp from the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), provided the first ever insights into the domain organisation and possible functional characteristics of the NiRAN domain. Using in silico tools, we predict that the NiRAN domain assumes a kinase or phosphotransferase like fold and binds nucleoside triphosphates at its proposed active site. Additionally, using molecular docking we have predicted the binding of three widely used kinase inhibitors and five well characterized anti-microbial compounds at the NiRAN domain active site along with their drug-likeliness. For the first time ever, using basic biochemical tools, this study shows the presence of a kinase like activity exhibited by the SARS-CoV-2 RdRp. Interestingly, a well-known kinase inhibitor- Sorafenib showed a significant inhibition and dampened viral load in SARS-CoV-2 infected cells. In line with the current global COVID-19 pandemic urgency and the emergence of newer strains with significantly higher infectivity, this study provides a new anti-SARS-CoV-2 drug target and potential lead compounds for drug repurposing against SARS-CoV-2.     

Genetic interactions among the West Nile virus methyltransferase, the RNA-dependent RNA polymerase, and the 5' stem-loop of genomic RNA

Flavivirus methyltransferase catalyzes both guanine N7 and ribose 2'-OH methylations of the viral RNA cap (GpppA-RNA-->m(7)GpppAm-RNA). The methyltransferase is physically linked to an RNA-dependent RNA polymerase (RdRp) in the flaviviral NS5 protein. Here, we report genetic interactions of West Nile virus (WNV) methyltransferase with the RdRp and the 5'-terminal stem-loop of viral genomic RNA. Genome-length RNAs, containing amino acid substitutions of D146 (a residue essential for both cap methylations) in the methyltransferase, were transfected into BHK-21 cells. Among the four mutant RNAs (D146L, D146P, D146R, and D146S), only D146S RNA generated viruses in transfected cells. Sequencing of the recovered viruses revealed that, besides the D146S change in the methyltransferase, two classes of compensatory mutations had reproducibly emerged. Class 1 mutations were located in the 5'-terminal stem-loop of the genomic RNA (a G35U substitution or U38 insertion). Class 2 mutations resided in NS5 (K61Q in methyltransferase and W751R in RdRp). Mutagenesis analysis, using a genome-length RNA and a replicon of WNV, demonstrated that the D146S substitution alone was lethal for viral replication; however, the compensatory mutations rescued replication, with the highest rescuing efficiency occurring when both classes of mutations were present. Biochemical analysis showed that a low level of N7 methylation of the D146S methyltransferase is essential for the recovery of adaptive viruses. The methyltransferase K61Q mutation facilitates viral replication through improved N7 methylation activity. The RdRp W751R mutation improves viral replication through an enhanced polymerase activity. Our results have clearly established genetic interactions among flaviviral methyltransferase, RdRp, and the 5' stem-loop of the genomic RNA.     

“Bucket brigade” using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2

The RNA-dependent RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a promising drug target for coronavirus disease 2019 (COVID-19) because it plays the most important role in the replication of the RNA genome. Nucleotide analogs such as remdesivir and favipiravir are thought to interfere with the RNA replication by RdRp. More specifically, they are expected to compete with nucleoside triphosphates, such as ATP. However, the process in which these drug molecules and nucleoside triphosphates are taken up by RdRp remains unknown. In this study, we performed all-atom molecular dynamics simulations to clarify the recognition mechanism of RdRp for these drug molecules and ATP that were at a distance. The ligand recognition ability of RdRp decreased in the order of remdesivir, favipiravir, and ATP. We also identified six recognition paths. Three of them were commonly found in all ligands, and the remaining three paths were ligand-dependent ones. In the common two paths, it was observed that the multiple lysine residues of RdRp carried the ligands to the binding site like a “bucket brigade.” In the remaining common path, the ligands directly reached the binding site. Our findings contribute to the understanding of the efficient ligand recognition by RdRp at the atomic level.     

“Bucket brigade” using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2

The RNA-dependent RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a promising drug target for coronavirus disease 2019 (COVID-19) because it plays the most important role in the replication of the RNA genome. Nucleotide analogs such as remdesivir and favipiravir are thought to interfere with the RNA replication by RdRp. More specifically, they are expected to compete with nucleoside triphosphates, such as ATP. However, the process in which these drug molecules and nucleoside triphosphates are taken up by RdRp remains unknown. In this study, we performed all-atom molecular dynamics simulations to clarify the recognition mechanism of RdRp for these drug molecules and ATP that were at a distance. The ligand recognition ability of RdRp decreased in the order of remdesivir, favipiravir, and ATP. We also identified six recognition paths. Three of them were commonly found in all ligands, and the remaining three paths were ligand-dependent ones. In the common two paths, it was observed that the multiple lysine residues of RdRp carried the ligands to the binding site like a “bucket brigade.” In the remaining common path, the ligands directly reached the binding site. Our findings contribute to the understanding of the efficient ligand recognition by RdRp at the atomic level.     

Cell Penetrable Humanized-VH/VHH That Inhibit RNA Dependent RNA Polymerase (NS5B) of HCV

NS5B is pivotal RNA dependent RNA polymerase (RdRp) of HCV and NS5B function interfering halts the virus infective cycle. This work aimed to produce cell penetrable humanized single domain antibodies (SdAb; VH/V_HH) that interfere with the RdRp activity. Recombinant NS5BΔ55 of genotype 3a HCV with de novo RNA synthetic activity was produced and used in phage biopanning for selecting phage clones that displayed NS5BΔ55 bound VH/V_HH from a humanized-camel VH/V_HH display library. VH/V_HH from E. coli transfected with four selected phage clones inhibited RdRp activity when tested by ELISA inhibition using 3′di-cytidylate 25 nucleotide directed in vitro RNA synthesis. Deduced amino acid sequences of two clones showed V_HH hallmark and were designated V_HH6 and V_HH24; other clones were conventional VH, designated VH9 and VH13. All VH/V_HH were linked molecularly to a cell penetrating peptide, penetratin. The cell penetrable VH9, VH13, V_HH6 and V_HH24 added to culture of Huh7 cells transfected with JHF-1 RNA of genotype 2a HCV reduced the amounts of RNA intracellularly and in culture medium implying that they inhibited the virus replication. VH/V_HH mimotopes matched with residues scattered on the polymerase fingers, palm and thumb which were likely juxtaposed to form conformational epitopes. Molecular docking revealed that the antibodies covered the RdRp catalytic groove. The transbodies await further studies for in vivo role in inhibiting HCV replication.     

VPg-primed RNA synthesis of norovirus RNA-dependent RNA polymerases by using a novel cell-based assay

Molecular studies of human noroviruses (NoV) have been hampered by the lack of a permissive cell culture system. We have developed a sensitive and reliable mammalian cell-based assay for the human NoV GII.4 strain RNA-dependent RNA polymerase (RdRp). The assay is based on the finding that RNAs synthesized by transiently expressed RdRp can stimulate retinoic acid-inducible gene I (RIG-I)-dependent reporter luciferase production via the beta interferon promoter. Comparable activities were observed for the murine norovirus (MNV) RdRp. RdRps with mutations at divalent metal ion binding residues did not activate RIG-I signaling. Furthermore, both NoV and MNV RdRp activities were stimulated by the coexpression of their respective VPg proteins, while mutations in the putative site of nucleotide linkage on VPg abolished most of their stimulatory effects. Sequencing of the RNAs linked to VPg revealed that the cellular trans-Golgi network protein 2 (TGOLN2) mRNA was the template for VPg-primed RNA synthesis. Small interfering RNA knockdown of RNase L abolished the enhancement of signaling that occurred in the presence of VPg. Finally, the coexpression of each of the other NoV proteins revealed that p48 (also known as NS1-2) and VP1 enhanced and that VP2 reduced the RdRp activity. The assay should be useful for the dissection of the requirements for NoV RNA synthesis as well as the identification of inhibitors of the NoV RdRp.     

Fidelity variants of RNA dependent RNA polymerases uncover an indirect, mutagenic activity of amiloride compounds

In a screen for RNA mutagen resistance, we isolated a high fidelity RNA dependent RNA polymerase (RdRp) variant of Coxsackie virus B3 (CVB3). Curiously, this variant A372V is also resistant to amiloride. We hypothesize that amiloride has a previously undescribed mutagenic activity. Indeed, amiloride compounds increase the mutation frequencies of CVB3 and poliovirus and high fidelity variants of both viruses are more resistant to this effect. We hypothesize that this mutagenic activity is mediated through alterations in intracellular ions such as Mg2+ and Mn2+, which in turn increase virus mutation frequency by affecting RdRp fidelity. Furthermore, we show that another amiloride-resistant RdRp variant, S299T, is completely resistant to this mutagenic activity and unaffected by changes in ion concentrations. We show that RdRp variants resist the mutagenic activity of amiloride via two different mechanisms: 1) increased fidelity that generates virus populations presenting lower basal mutation frequencies or 2) resisting changes in divalent cation concentrations that affect polymerase fidelity. Our results uncover a new antiviral approach based on mutagenesis.     

Poly(A)- and primer-independent RNA polymerase of Norovirus

Replication of positive-strand caliciviruses is mediated by a virus-encoded RNA-dependent RNA polymerase (RdRp). To study the replication of Norovirus (NV), a member of the family Caliciviridae, we used a recombinant baculovirus system to express an enzymatically active RdRp protein from the 3D region of the NV genome and defined conditions for optimum enzymatic activity. Using an RNA template from the NV 3' genomic region, we observed similar levels of enzymatic activity in assays with and without a poly(A) tail. RdRp activity was not significantly affected by the addition of an RNA primer to the reaction mixture. Thus, the NV RdRp exhibited primer- and poly(A)-independent RNA polymerase activity. While the RdRp inhibitor phosphonoacetic acid inhibited NV RdRp activity, another gliotoxin did not. The active recombinant NV RdRp will be of benefit to studies of NV replication and will facilitate the development of specific inhibitors of NV proliferation.     

RNA-dependent RNA polymerase activity associated with endogenous double-stranded RNA in rice

RNA-dependent RNA polymerase (RdRp) activity was detected in the crude microsomal fraction of rice cultured cells that contain a 14 kbp double-stranded RNA (dsRNA). RdRp activity is maximal in the presence of all four nucleotide triphosphates and Mg2+ ion and is resistant to inhibitors of DNA-dependent RNA polymerases (actinomycin D and alpha-amanitin). RdRp activity increases approximately 2.5-fold in the presence of 0.5% deoxycholate. Treatment of purified microsomal fraction with proteinase K plus deoxycholate suggests that the RdRp enzyme complex with its own 14 kb RNA template is located in vesicles. The RdRp enzyme complex was solubilized with Nonidet P-40 and purified by glycerol gradient centrifugation, then exogenous RNA templates were added. Results indicate that exogenous dsRNA reduces RNA synthesis from the endogenous 14 kb RNA template.     

NMPylation and de-NMPylation of SARS-CoV-2 nsp9 by the NiRAN domain

The catalytic subunit of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) contains two active sites that catalyze nucleotidyl-monophosphate transfer (NMPylation). Mechanistic studies and drug discovery have focused on RNA synthesis by the highly conserved RdRp. The second active site, which resides in a Nidovirus RdRp-Associated Nucleotidyl transferase (NiRAN) domain, is poorly characterized, but both catalytic reactions are essential for viral replication. One study showed that NiRAN transfers NMP to the first residue of RNA-binding protein nsp9; another reported a structure of nsp9 containing two additional N-terminal residues bound to the NiRAN active site but observed NMP transfer to RNA instead. We show that SARS-CoV-2 RdRp NMPylates the native but not the extended nsp9. Substitutions of the invariant NiRAN residues abolish NMPylation, whereas substitution of a catalytic RdRp Asp residue does not. NMPylation can utilize diverse nucleotide triphosphates, including remdesivir triphosphate, is reversible in the presence of pyrophosphate, and is inhibited by nucleotide analogs and bisphosphonates, suggesting a path for rational design of NiRAN inhibitors. We reconcile these and existing findings using a new model in which nsp9 remodels both active sites to alternately support initiation of RNA synthesis by RdRp or subsequent capping of the product RNA by the NiRAN domain.     

Identification of SARS-CoV-2 RNA-dependent RNA polymerase inhibitors from the major phytochemicals of Nigella sativa: An in silico approach

The coronavirus disease 2019 (COVID-19), which emerged in December 2019, continues to be a serious health concern worldwide. There is an urgent need to develop effective drugs and vaccines to control the spread of this disease. In the current study, the main phytochemical compounds of Nigella sativa were screened for their binding affinity for the active site of the RNA-dependent RNA polymerase (RdRp) enzyme of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The binding affinity was investigated using molecular docking methods, and the interaction of phytochemicals with the RdRp active site was analyzed and visualized using suitable software. Out of the nine phytochemicals of N. sativa screened in this study, a significant docking score was observed for four compounds, namely α-hederin, dithymoquinone, nigellicine, and nigellidine. Based on the findings of our study, we report that α-hederin, which was found to possess the lowest binding energy (–8.6 kcal/mol) and hence the best binding affinity, is the best inhibitor of RdRp of SARS-CoV-2, among all the compounds screened here. Our results prove that the top four potential phytochemical molecules of N. sativa, especially α-hederin, could be considered for ongoing drug development strategies against SARS-CoV-2. However, further in vitro and in vivo testing are required to confirm the findings of this study.     

Structure-function analysis of mutant RNA-dependent RNA polymerase complexes with VPg

The replication of the foot-and-mouth disease virus (FMDV) genome is critically dependent upon the activity of a virally encoded RNA-dependent RNA polymerase (RdRp). In this study, four mutant RdRps of FMDV were isolated from viral quasi-species treated with ribavirin, of which two were single mutants (L123F and T381A) and two were double mutants (T291I/T381I and L123F/F244L). The mutant proteins were expressed in Escherichia coli and purified by His-bind resin chromatography. In combination with real-time RT-PCR, an in vitro RNA replication system that uses genome RNA/VPg as template-primers was used to determine polymerase activity. Mutant L123F exhibited a 0.6-fold decrease (p < 0.001) in polymerase activity relative to wild-type RdRp, whereas the activity of L123F/F244L and T381A was undetectable. Surprisingly, the activity of T291I/T381I yielded a 0.7-fold increase (p < 0.001) as compared to wild-type. In order to study the structure-function relationship of RdRp, all structures of the RdRp-RNA template-primer complex were obtained through homology modeling and molecular docking. The VPg1 orientation in the RdRp-VPg1 complexes was determined and analyzed with mathematical methods. Our results reveal that the orientation of VPg after binding to the polymerase determines the FMDV RdRp catalytic activity, which provides a basis for the rational design of novel antiviral agents.     

Initiation of minus-strand RNA synthesis by the brome mosaicvirus RNA-dependent RNA polymerase: use of oligoribonucleotide primers.

Various DNA- and RNA-dependent RNA polymerases have been reported to use oligoribonucleotide primers to initiate nucleic acid synthesis. For the brome mosaic virus RNA-dependent RNA polymerase (RdRp), we determined that in reactions performed with limited GTP concentrations, minus-strand RNA synthesis can be stimulated by the inclusion of guanosine monophosphate or specific oligoribonucleotides. Furthermore, guanylyl-3',5'-guanosine (GpG) was incorporated into minus-strand RNA and increased the rate of minus-strand RNA synthesis. In the presence of GpG, RdRp's Km for GTP decreased from 50 microM to approximately 3 microM while the Kms for other nucleotides were unaffected. These results have implications for the mechanism of initiation by RdRp.     

Nodavirus RNA replication protein a induces membrane association of genomic RNA

Positive-strand RNA virus genome replication occurs in membrane-associated RNA replication complexes, whose assembly remains poorly understood. Here we show that prior to RNA replication, the multifunctional, transmembrane RNA replication protein A of the nodavirus flock house virus (FHV) recruits FHV genomic RNA1 to a membrane-associated state in both Drosophila melanogaster and Saccharomyces cerevisiae cells. Protein A has mitochondrial membrane-targeting, self-interaction, RNA-dependent RNA polymerase (RdRp), and RNA capping domains. In the absence of RdRp activity due to an active site mutation (A(D692E)), protein A stimulated RNA1 accumulation by increasing RNA1 stability. Protein A(D692E) stimulated RNA1 accumulation in wild-type cells and in xrn1(-) yeast defective in decapped RNA decay, showing that increased RNA1 stability was not due to protein A-mediated RNA1 recapping. Increased RNA1 stability was closely linked with protein A-induced membrane association of the stabilized RNA and was highly selective for RNA1. Substantial N- and C-proximal regions of protein A were dispensable for these activities. However, increased RNA1 accumulation was eliminated by deleting protein A amino acids (aa) 1 to 370 but was restored completely by adding back the transmembrane domain (aa 1 to 35) and partially by adding back peripheral membrane association sequences in aa 36 to 370. Moreover, although RNA polymerase activity was not required, even small deletions in or around the RdRp domain abolished increased RNA1 accumulation. These and other results show that prior to negative-strand RNA synthesis, multiple domains of mitochondrially targeted protein A cooperate to selectively recruit FHV genomic RNA to membranes where RNA replication complexes form.     

Molecular dynamics simulations of viral RNA polymerases link conserved and correlated motions of functional elements to fidelity

The viral RNA-dependent RNA polymerase (RdRp) is essential for multiplication of all RNA viruses. The sequence diversity of an RNA virus population contributes to its ability to infect the host. This diversity emanates from errors made by the RdRp during RNA synthesis. The physical basis for RdRp fidelity is unclear but is linked to conformational changes occurring during the nucleotide-addition cycle. To understand RdRp dynamics that might influence RdRp function, we have analyzed all-atom molecular dynamics simulations on the nanosecond timescale of four RdRps from the picornavirus family that exhibit 30-74% sequence identity. Principal component analysis showed that the major motions observed during the simulations derived from conserved structural motifs and regions of known function. The dynamics of residues participating in the same biochemical property, for example, RNA binding, nucleotide binding or catalysis, were correlated even when spatially distant on the RdRp structure. The conserved and correlated dynamics of functional structural elements suggest coevolution of dynamics with structure and function of the RdRp. Crystal structures of all picornavirus RdRps exhibit a template-nascent RNA duplex channel too small to fully accommodate duplex RNA. Simulations revealed opening and closing motions of the RNA and nucleoside triphosphate channels, which might be relevant to nucleoside triphosphate entry, inorganic pyrophosphate exit and translocation. A role for nanosecond timescale dynamics in RdRp fidelity is supported by the altered dynamics of the high-fidelity G64S derivative of PV RdRp relative to wild-type enzyme.