Molecular and Structural Aspects of Clinically Relevant Mutations of SARS-CoV-2 RNA-Dependent RNA Polymerase in Remdesivir-Treated Patients

(1) Background: SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) is a promising therapeutic target to fight COVID-19, and many RdRp inhibitors nucleotide/nucleoside analogs, such as remdesivir, have been identified or are in clinical studies. However, the appearance of resistant mutations could reduce their efficacy. In the present work, we structurally evaluated the impact of RdRp mutations found at baseline in 39 patients treated with remdesivir and associated with a different degree of antiviral response in vivo. (2) Methods: A refined bioinformatics approach was applied to assign SARS-CoV-2 clade and lineage, and to define RdRp mutational profiles. In line with such a method, the same mutations were built and analyzed by combining docking and thermodynamics evaluations with both molecular dynamics and representative pharmacophore models. (3) Results: Clinical studies revealed that patients bearing the most prevalent triple mutant P323L+671S+M899I, which was present in 41% of patients, or the more complex mutational profile P323L+G671S+L838I+D738Y+K91E, which was found with a prevalence of 2.6%, showed a delayed reduced response to remdesivir, as confirmed by the increase in SARS-CoV-2 viral load and by a reduced theoretical binding affinity versus RdRp (ΔGbindWT = -122.70 kcal/mol; ΔGbindP323L+671S+M899I = -84.78 kcal/mol; ΔGbindP323L+G671S+L838I+D738Y+K91E = -96.74 kcal/mol). Combined computational approaches helped to rationalize such clinical observations, offering a mechanistic understanding of the allosteric effects of mutants on the global motions of the viral RNA synthesis machine and in the changes of the interactions patterns of remdesivir during its binding.

Galidesivir Triphosphate Promotes Stalling of Dengue-2 Virus Polymerase Immediately Prior to Incorporation

Millions of people are infected by the dengue and Zika viruses each year, resulting in significant morbidity and mortality. Galidesivir is an adenosine nucleoside analog that can attenuate flavivirus replication in cell-based assays and animal models of infection. Galidesivir is converted to the triphosphorylated form by host kinases and subsequently incorporated into viral RNA by viral RNA polymerases. This has been proposed to lead to the delayed termination of RNA synthesis. Here, we report direct in vitro testing of the effects of Galidesivir triphosphate on dengue-2 and Zika virus polymerase activity. Galidesivir triphosphate was chemically synthesized, and inhibition of RNA synthesis followed using a dinucleotide-primed assay with a homopolymeric poly(U) template. Galidesivir triphosphate was equipotent against dengue-2 and Zika polymerases, with IC50 values of 42 ± 12 μM and 47 ± 5 μM, respectively, at an ATP concentration of 20 μM. RNA primer extension assays show that the dengue-2 polymerase stalls while attempting to add a Galidesivir nucleotide to the nascent RNA chain, evidenced by the accumulation of RNA products truncated immediately upstream of Galidesivir incorporation sites. Nevertheless, Galidesivir is incorporated at isolated sites with low efficiency, leading to the subsequent synthesis of full-length RNA with no evidence of delayed chain termination. The incorporation of Galidesivir at consecutive sites is strongly disfavored, highlighting the potential for modulation of inhibitory effects of nucleoside analogs by the template sequence. Our results suggest that attenuation of dengue replication by Galidesivir may not derive from the early termination of RNA synthesis following Galidesivir incorporation.

Celastrol attenuates hepatitis C virus translation and inflammatory response in mice by suppressing heat shock protein 90β

Hepatitis C virus (HCV) infection is one of the major factors to trigger a sustained hepatic inflammatory response and hence hepatocellular carcinoma (HCC), but direct-acting-antiviral (DAAs) was not efficient to suppress HCC development. Heat shock protein 90 kDa (HSP90) is highly abundant in different types of cancers, and especially controls protein translation, endoplasmic reticulum stress, and viral replication. In this study we investigated the correlation between the expression levels of HSP90 isoforms and inflammatory response marker NLRP3 in different types of HCC patients as well as the effect of a natural product celastrol in suppression of HCV translation and associated inflammatory response in vivo. We identified that the expression level of HSP90β isoform was correlated with that of NLRP3 in the liver tissues of HCV positive HCC patients (R2 = 0.3867, P < 0.0101), but not in hepatitis B virus-associated HCC or cirrhosis patients. We demonstrated that celastrol (3, 10, 30 μM) dose-dependently suppressed the ATPase activity of both HSP90α and HSP90β, while its anti-HCV activity was dependent on the Ala47 residue in the ATPase pocket of HSP90β. Celastrol (200 nM) halted HCV internal ribosomal entry site (IRES)-mediated translation at the initial step by disrupting the association between HSP90β and 4EBP1. The inhibitory activity of celastrol on HCV RNA-dependent RNA polymerase (RdRp)-triggered inflammatory response also depended on the Ala47 residue of HSP90β. Intravenous injection of adenovirus expressing HCV NS5B (pAde-NS5B) in mice induced severe hepatic inflammatory response characterized by significantly increased infiltration of immune cells and hepatic expression level of Nlrp3, which was dose-dependently ameliorated by pretreatment with celastrol (0.2, 0.5 mg/kg, i.p.). This study reveals a fundamental role of HSP90β in governing HCV IRES-mediated translation as well as hepatic inflammation, and celastrol as a novel inhibitor of HCV translation and associated inflammation by specifically targeting HSP90β, which could be developed as a lead for the treatment of HSP90β positive HCV-associated HCC.

Protective Effects of Caffeine on Chikungunya and Zika Virus Infections: An in Vitro and in Silico Study

Infection by viruses Chikungunya (CHIKV) and Zika (ZIKV) continue to be serious problems in tropical and subtropical areas of the world. Here, we evaluated the antiviral and virucidal activity of caffeine against CHIKV and ZIKV in Vero, A549, and Huh-7 cell lines. Results showed that caffeine displays antiviral properties against both viruses. By pre-and post-infection treatment, caffeine significantly inhibited CHIKV and ZIKV replication in a dose-dependent manner. Furthermore, caffeine showed a virucidal effect against ZIKV. Molecular docking suggests the possible binding of caffeine with envelope protein and RNA-dependent RNA polymerase of CHIKV and ZIKV. This is the first study that showed an antiviral effect of caffeine against CHIKV and ZIKV. Although further studies are needed to better understand the mechanism of caffeine-mediated repression of viral replication, caffeine appears to be a promising compound that could be used for in vivo studies, perhaps in synergy with other compounds present in daily beverages.

Investigating the Mechanism of Action of Anti-Dengue Compounds as Potential Binders of Zika Virus RNA-Dependent RNA Polymerase

The World Health Organization (WHO) has designated the Zika virus (ZIKV) as a significant risk to the general public's health. Currently, there are no vaccinations or medications available to treat or prevent infection with the Zika virus. Thus, it is urgently required to develop a highly efficient therapeutic molecule. In the presented study, a computationally intensive search was carried out to identify potent compounds that have the potential to bind and block the activity of ZIKV NS5 RNA-dependent RNA polymerase (RdRp). The anti-dengue chemical library was subjected to high-throughput virtual screening and MM/GBSA analysis in order to rate the potential candidates. The top three compounds were then chosen. According to the MM/GBSA analysis, compound 127042987 from the database had the highest binding affinity to the protein with a minimum binding free energy of -77.16 kcal/mole. Compound 127042987 had the most stable RMSD trend and the greatest number of hydrogen bond interactions when these chemical complexes were evaluated further under a 100 ns molecular dynamics simulation. Compound 127042987 displayed the best binding free energy (GBind) of -96.50 kcal/mol, surpassing the native ligand binding energy (-66.17 kcal/mole). Thereafter, an MM/GBSA binding free energy study was conducted to validate the stability of selected chemical complexes. Overall, this study illustrated that compound 127042987 showed preferred binding free energies, suggesting a possible inhibitory mechanism against ZIKV-RdRp. As per this study, it was proposed that compound 127042987 could be used as a therapeutic option to prevent Zika virus infection. These compounds need to be tested in experiments for further validation.

Mutual information networks reveal evolutionary relationships within the influenza A virus polymerase

The influenza A virus (IAV) RNA polymerase is an essential driver of IAV evolution. Mutations that the polymerase introduces into viral genome segments during replication are the ultimate source of genetic variation, including within the three subunits of the IAV polymerase (polymerase basic protein 2, polymerase basic protein 1, and polymerase acidic protein). Evolutionary analysis of the IAV polymerase is complicated, because changes in mutation rate, replication speed, and drug resistance involve epistatic interactions among its subunits. In order to study the evolution of the human seasonal H3N2 polymerase since the 1968 pandemic, we identified pairwise evolutionary relationships among ∼7000 H3N2 polymerase sequences using mutual information (MI), which measures the information gained about the identity of one residue when a second residue is known. To account for uneven sampling of viral sequences over time, we developed a weighted MI (wMI) metric and demonstrate that wMI outperforms raw MI through simulations using a well-sampled severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) dataset. We then constructed wMI networks of the H3N2 polymerase to extend the inherently pairwise wMI statistic to encompass relationships among larger groups of residues. We included hemagglutinin (HA) in the wMI network to distinguish between functional wMI relationships within the polymerase and those potentially due to hitch-hiking on antigenic changes in HA. The wMI networks reveal coevolutionary relationships among residues with roles in replication and encapsidation. Inclusion of HA highlighted polymerase-only subgraphs containing residues with roles in the enzymatic functions of the polymerase and host adaptability. This work provides insight into the factors that drive and constrain the rapid evolution of influenza viruses.

A potential role for Giardia chaperone protein GdDnaJ in regulating Giardia proliferation and Giardiavirus replication

BACKGROUND

Giardia duodenalis (referred to as Giardia) is a flagellated binucleate protozoan parasite, which causes one of the most common diarrheal diseases, giardiasis, worldwide. Giardia can be infected by Giardiavirus (GLV), a small endosymbiotic dsRNA virus belongs to the Totiviridae family. However, the regulation of GLV and a positive correlation between GLV and Giardia virulence is yet to be elucidated.

METHODS

To identify potential regulators of GLV, we performed a yeast two-hybrid (Y2H) screen to search for interacting proteins of RdRp. GST pull-down, co-immunoprecipitation and bimolecular fluorescence complementation (BiFC) assay were used to verify the direct physical interaction between GLV RdRp and its new binding partner. In addition, their in vivo interaction and colocalization in Giardia trophozoites were examined by using Duolink proximal ligation assay (Duolink PLA).

RESULTS

From Y2H screen, the Giardia chaperone protein, Giardia DnaJ (GdDnaJ), was identified as a new binding partner for GLV RdRp. The direct interaction between GdDnaJ and GLV RdRp was verified via GST pull-down, co-immunoprecipitation and BiFC. In addition, colocalization and in vivo interaction between GdDnaJ and RdRp in Giardia trophozoites were confirmed by Duolink PLA. Further analysis revealed that KNK437, the inhibitor of GdDnaJ, can significantly reduce the replication of GLVs and the proliferation of Giardia.

CONCLUSION

Taken together, our results suggested a potential role of GdDnaJ in regulating Giardia proliferation and GLV replication through interaction with GLV RdRp.

First report of Capsicum chlorosis virus naturally infecting Ageratum conyzoides in China

Capsicum chlorosis virus (CaCV; family Tospoviridae, genus Orthotospovirus) was first reported to infect capsicum (Capsicum annuum) and tomato (Solanum lycopersicum) in Australia in 2002 (McMichael et al., 2002). Subsequently, its infection was detected in different plants including waxflower (Hoya calycina Schlecter) in the United States (Melzer et al. 2014), peanut (Arachis hypogaea) in India (Vijayalakshmi et al. 2016), and spider lily (Hymenocallis americana) (Huang et al. 2017), Chilli pepper (Capsicum annuum) (Zheng et al. 2020), and Feiji cao (Chromolaena odorata) (Chen et al. 2022) in China. Ageratum conyzoides L. (commonly known as goat weed, family Asteraceae) is a natural weed in crop fields distributed in subtropical and tropical areas and a reservoir host of numerous plant pathogens (She et al. 2013). In April 2022, we observed that 90% of plants of A. conyzoides in maize fields in Sanya, Hainan province, China, exhibited typical virus-like symptoms of vein yellowing, leaf chlorosis, and distortion (Fig. S1 A-C). Total RNA was extracted from one symptomatic leaf of A. conyzoides. Small RNA libraries were constructed using the small RNA Sample Pre Kit (Illumina, San Diego, USA) for sequencing with an Illumina Novaseq 6000 platform (Biomarker Technologies Corporation, Beijing, China). A total 15,848,189 clean reads were obtained after removing low-quality reads. Quality-controlled qualified reads were assembled into contigs using Velvet 1.0.5 software with a k-mer value of 17. One hundred contigs shared nucleotide identity ranging from 85.7% to 100% with the CaCV using BLASTn searches online (https://blast.ncbi.nlm.nih.gov/Blast.cgi?). Numerous contigs (45, 34, and 21) obtained in this study were mapped to the L, M, and S RNA segments of the CaCV-Hainan isolate (GenBank accession no. KX078565- KX078567) from spider lily (Hymenocallis americana) in Hainan province, China, respectively. The full-length of L, M, and S RNA segments of CaCV-AC were determined to be 8,913, 4,841, and 3,629 bp, respectively (GenBank accession no. OQ597167- OQ597169). Furthermore, five symptomatic leaf samples were tested to be positive for CaCV using a CaCV enzyme-linked immunosorbent assay (ELISA) kit (MEIMIAN, Jiangsu, China) (Fig. S1-D). Total RNA from these leaves was amplified by RT-PCR with two sets of primer pairs. Primers CaCV-F (5'-ACTTTCCATCAACCTCTGT-3') and CaCV-R (5'-GTTATGGCCATATTTCCCT-3') were used for the amplification of 828 bp fragment from nucleocapsid protein (NP) on CaCV S RNA. While another, primers gL3637 (5'-CCTTTAACAGTDGAAACAT-3') and gL4435c (5'-CATDGCRCAAGARTGRTARACAGA-3') were used for the amplification of 816 bp fragment from RNA-dependent RNA polymerase (RdRP) on CaCV L RNA (Fig. S1-E and -F) (Basavaraj et al. 2020). These amplicons were cloned into the pCE2 TA/Blunt-Zero vector (Vazyme, Nanjing, China) and three independent positive colonies of Escherichia coli DH5α carrying each viral amplicon were sequenced. These sequences were deposited in the GenBank database under accession nos. OP616700-OP616709. Pairwise sequence comparison revealed that nucleotide sequences of NP and RdRP genes of the five CaCV isolates shared 99.5% (812 bp out of 828 bp) and 99.4% (799 bp out of 816 bp) nucleotide identities, respectively. They showed 86.2-99.2% and 86.5-99.1% nucleotide identities with corresponding nucleotide sequences of other CaCV isolates derived from GenBank database, respectively. The highest nucleotide sequence identity (99%) of the CaCV isolates obtained in the study was observed with the CaCV-Hainan isolate. Phylogenetic analysis based on NP amino acid demonstrated that six CaCV isolates (this study = 5 and NCBI database = 1) clustered into one distinct clade (Fig. S2). Our data confirmed for the first time the presence of CaCV naturally infecting A. conyzoides plant in China, which enriches information on the host range and will be helpful for disease management.

Current trends in RNA virus detection through metatranscriptome sequencing data

With advances in sequencing technology, metatranscriptome sequencing from a variety of environmental and biological sources has revealed the existence of various previously unknown RNA viruses. This review presents recent major RNA virome studies sampled from invertebrate and vertebrate species as well as aquatic environments. In particular, we focus on severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and related RNA virus identification through metatranscriptome sequencing analyses. Recently developed bioinformatics software and databases for RNA virus identification are introduced. A relationship between newly identified RNA viruses and endogenous viral elements in host genomes is also discussed.

A Single Nonsynonymous Substitution in the RNA-Dependent RNA Polymerase of Potato virus Y Allows the Simultaneous Breakdown of Two Different Forms of Antiviral Resistance in Capsicum annuum

The dominant Pvr4 gene in pepper (Capsicum annuum) confers resistance to members of six potyvirus species, all of which belong to the Potato virus Y (PVY) phylogenetic group. The corresponding avirulence factor in the PVY genome is the NIb cistron (i.e., RNA-dependent RNA polymerase). Here, we describe a new source of potyvirus resistance in the Guatemalan accession C. annuum cv. PM949. PM949 is resistant to members of at least three potyvirus species, a subset of those controlled by Pvr4. The F₁ progeny between PM949 and the susceptible cultivar Yolo Wonder was susceptible to PVY, indicating that the resistance is recessive. The segregation ratio between resistant and susceptible plants observed in the F₂ progeny matched preferably with resistance being determined by two unlinked recessive genes independently conferring resistance to PVY. Inoculations by grafting resulted in the selection of PVY mutants breaking PM949 resistance and, less efficiently, Pvr4-mediated resistance. The codon substitution E₄₇₂K in the NIb cistron of PVY, which was shown previously to be sufficient to break Pvr4 resistance, was also sufficient to break PM949 resistance, a rare example of cross-pathogenicity effect. In contrast, the other selected NIb mutants showed specific infectivity in PM949 or Pvr4 plants. Comparison of Pvr4 and PM949 resistance, which share the same target in PVY, provides interesting insights into the determinants of resistance durability.

Molecular Epidemiology of SARS-CoV-2 Omicron Sub-Lineages Isolated from Turkish Patients Infected with COVID-19

Early detection and characterization of new variants and their impacts enable improved genomic surveillance. This study aims to evaluate the subvariant distribution of Omicron strains isolated from Turkish cases to determine the rate of antiviral resistance of RdRp and 3CLpro inhibitors. The Stanford University Coronavirus Antiviral & Resistance Database online tool was used for variant analyses of the strains uploaded to GISAID as Omicron (n = 20.959) between January 2021 and February,2023. Out of 288 different Omicron subvariants, B.1, BA.1, BA.2, BA.4, BE.1, BF.1, BM.1, BN.1, BQ.1, CK.1, CL.1, and XBB.1 were the main determined subvariants, and BA.1 (34.7%), BA.2 (30.8%), and BA.5 (23.6%) were reported most frequently. RdRp and 3CLPro-related resistance mutations were determined in n = 150, 0.72% sequences, while the rates of resistance against RdRp and 3CLpro inhibitors were reported at 0.1% and 0.6%, respectively. Mutations that were previously associated with a reduced susceptibility to remdesivir, nirmatrelvir/r, and ensitrelvir were most frequently detected in BA.2 (51.3%). The mutations detected at the highest rate were A449A/D/G/V (10.5%), T21I (10%), and L50L/F/I/V (6%). Our findings suggest that continuous monitoring of variants, due to the diversity of Omicron lineages, is necessary for global risk assessment. Although drug-resistant mutations do not pose a threat, the tracking of drug mutations will be necessary due to variant heterogenicity.

Profiling Plant Proteome and Transcriptome Changes during Grapevine Fanleaf Virus Infection

Viruses can elicit varying types and severities of symptoms during plant host infection. We investigated changes in the proteome and transcriptome of Nicotiana benthamiana plants infected by grapevine fanleaf virus (GFLV) with an emphasis on vein clearing symptom development. Comparative, time-course liquid chromatography tandem mass spectrometry and 3' ribonucleic acid sequencing analyses of plants infected by two wildtype GFLV strains, one symptomatic and one asymptomatic, and their asymptomatic mutant strains carrying a single amino acid change in the RNA-dependent RNA polymerase (RdRP) were conducted to identify host biochemical pathways involved in viral symptom development. During peak vein clearing symptom display at 7 days post-inoculation (dpi), protein and gene ontologies related to immune response, gene regulation, and secondary metabolite production were overrepresented when contrasting wildtype GFLV strain GHu and mutant GHu-1E_K802G^Pol. Prior to the onset of symptom development at 4 dpi and when symptoms faded away at 12 dpi, protein and gene ontologies related to chitinase activity, hypersensitive response, and transcriptional regulation were identified. This systems biology approach highlighted how a single amino acid of a plant viral RdRP mediates changes to the host proteome (∼1%) and transcriptome (∼8.5%) related to transient vein clearing symptoms and the network of pathways involved in the virus-host arms race.

Effects of Epitranscriptomic RNA Modifications on the Catalytic Activity of the SARS-CoV-2 Replication Complex

SARS-CoV-2 causes individualized symptoms. Many reasons have been given. We propose that an individual's epitranscriptomic system could be responsible as well. The viral RNA genome can be subject to epitranscriptomic modifications, which can be different for different individuals, and thus epitranscriptomics can affect many events including RNA replication differently. In this context, we studied the effects of modifications including pseudouridine (Ψ), 5-methylcytosine (m5 C), N6-methyladenosine (m6 A), N1-methyladenosine (m1 A) and N3-methylcytosine (m3 C) on the activity of SARS-CoV-2 replication complex (SC2RC). We found that Ψ, m5 C, m6 A and m3 C had little effect, whereas m1 A inhibited the enzyme. Both m1 A and m3 C disrupt canonical base pairing, but they had different effects. The fact that m1 A inhibits SC2RC implies that the modification can be difficult to detect. This fact also implies that individuals with upregulated m1 A including cancer, obesity and diabetes patients might have milder symptoms. However, this contradicts clinical observations. Relevant discussions are provided.

The Unusual Architecture of RNA-Dependent RNA Polymerase (RdRp)'s Catalytic Chamber Provides a Potential Strategy for Combination Therapy against COVID-19

The unusual and interesting architecture of the catalytic chamber of the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) was recently explored using Cryogenic Electron Microscopy (Cryo-EM), which revealed the presence of two distinctive binding cavities within the catalytic chamber. In this report, first, we mapped out and fully characterized the variations between the two binding sites, BS1 and BS2, for significant differences in their amino acid architecture, size, volume, and hydrophobicity. This was followed by investigating the preferential binding of eight antiviral agents to each of the two binding sites, BS1 and BS2, to understand the fundamental factors that govern the preferential binding of each drug to each binding site. Results showed that, in general, hydrophobic drugs, such as remdesivir and sofosbuvir, bind better to both binding sites than relatively less hydrophobic drugs, such as alovudine, molnupiravir, zidovudine, favilavir, and ribavirin. However, suramin, which is a highly hydrophobic drug, unexpectedly showed overall weaker binding affinities in both binding sites when compared to other drugs. This unexpected observation may be attributed to its high binding solvation energy, which disfavors overall binding of suramin in both binding sites. On the other hand, hydrophobic drugs displayed higher binding affinities towards BS1 due to its higher hydrophobic architecture when compared to BS2, while less hydrophobic drugs did not show a significant difference in binding affinities in both binding sites. Analysis of binding energy contributions revealed that the most favorable components are the ΔE_ele, ΔE_vdw, and ΔG_gas, whereas ΔG_sol was unfavorable. The ΔE_ele and ΔG_gas for hydrophobic drugs were enough to balance the unfavorable ΔG_sol, leaving the ΔE_vdw to be the most determining factor of the total binding energy. The information presented in this report will provide guidelines for tailoring SARS-CoV-2 inhibitors with enhanced binding profiles.

Safety and Biodistribution of Nanoligomers Targeting the SARS-CoV-2 Genome for the Treatment of COVID-19

As the world braces to enter its fourth year of the coronavirus disease 2019 (COVID-19) pandemic, the need for accessible and effective antiviral therapeutics continues to be felt globally. The recent surge of Omicron variant cases has demonstrated that vaccination and prevention alone cannot quell the spread of highly transmissible variants. A safe and nontoxic therapeutic with an adaptable design to respond to the emergence of new variants is critical for transitioning to the treatment of COVID-19 as an endemic disease. Here, we present a novel compound, called SBCoV202, that specifically and tightly binds the translation initiation site of RNA-dependent RNA polymerase within the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome, inhibiting viral replication. SBCoV202 is a Nanoligomer, a molecule that includes peptide nucleic acid sequences capable of binding viral RNA with single-base-pair specificity to accurately target the viral genome. The compound has been shown to be safe and nontoxic in mice, with favorable biodistribution, and has shown efficacy against SARS-CoV-2 in vitro. Safety and biodistribution were assessed using three separate administration methods, namely, intranasal, intravenous, and intraperitoneal. Safety studies showed the Nanoligomer caused no outward distress, immunogenicity, or organ tissue damage, measured through observation of behavior and body weight, serum levels of cytokines, and histopathology of fixed tissue, respectively. SBCoV202 was evenly biodistributed throughout the body, with most tissues measuring Nanoligomer concentrations well above the compound KD of 3.37 nM. In addition to favorable availability to organs such as the lungs, lymph nodes, liver, and spleen, the compound circulated through the blood and was rapidly cleared through the renal and urinary systems. The favorable biodistribution and lack of immunogenicity and toxicity set Nanoligomers apart from other antisense therapies, while the adaptability of the nucleic acid sequence of Nanoligomers provides a defense against future emergence of drug resistance, making these molecules an attractive potential treatment for COVID-19.

Exploration of potent antiviral phytomedicines from Lauraceae family plants against SARS-CoV-2 RNA-dependent RNA polymerase

RNA-dependent RNA polymerase, also known as RdRp, is a possible therapeutic target that could be used to suppress the proliferation of RNA viruses such as SARS-CoV-2. This protein has two major functional sites (a) catalytic and (b) substrate entry, which regulate the natural substrate entry and its corresponding interaction with the protein. In this study, a computational drug design pipeline was applied to investigate potential inhibitors against SARS-CoV-2 RdRp from Lauraceae plants, and five top hits were selected based on the docked score (< -7 kcal/mol). The docking study suggested that the Glochidioboside had a minimum binding score of -7.8 kcal/mol. This compound showed total five hydrogen bonds while two of them were with catalytic residues Asp618 and Asp760. However, another compound, Sitogluside showed a binding score of -7.3 kcal/mol with four hydrogen bonds targeting three functional residues (Arg555, Ser759, and Asp760). Later, 100 ns explicit solvent molecular dynamics (MD) simulation was performed to evaluate the stability of the protein-ligand docked system. These compounds translocated their positions from the catalytic site to the substrate entry site, as observed in the MD simulation trajectory. However, translocation did not affect the binding strength of these compounds, and they retained the strong binding affinity (ΔG < -11.5 kcal/mol), estimated using the MM/GBSA method. In general, the findings of this study indicated the potential therapeutic compounds that may be used targeting SARS-CoV-2 RdRp. However, these compounds still need to be validated by experimentation in order to determine their inhibitory function.Communicated by Ramaswamy H. Sarma.

Fidelity of Ribonucleotide Incorporation by the SARS-CoV-2 Replication Complex

The SARS-CoV-2 coronavirus has caused a global pandemic. Despite the initial success of vaccines at preventing infection, genomic variation has led to the proliferation of variants capable of higher infectivity. Mutations in the SARS-CoV-2 genome are the consequence of replication errors, highlighting the importance of understanding the determinants of SARS-CoV-2 replication fidelity. The RNA-dependent RNA polymerase (RdRp) is the central catalytic subunit for SARS-CoV-2 RNA replication and genome transcription. Here, we report the fidelity of ribonucleotide incorporation by SARS-CoV-2 RdRp (nsp12), along with its co-factors nsp7/nsp8, using steady-state kinetic analysis. Our analysis suggests that in the absence of the proofreading subunit (nsp14), the nsp12/7/8 complex has a surprisingly low base substitution fidelity (10-1-10-3). This is orders of magnitude lower than the fidelity reported for other coronaviruses (10-6-10-7), highlighting the importance of proofreading for faithful SARS-CoV-2 replication. We performed a mutational analysis of all reported SARS-CoV-2 genomes and identified mutations in both nsp12 and nsp14 that appear likely to lower viral replication fidelity through mechanisms that include impairing the nsp14 exonuclease activity or its association with the RdRp. Our observations provide novel insight into the mechanistic basis of replication fidelity in SARS-CoV-2 and the potential effect of nsp12 and nsp14 mutations on replication fidelity, informing the development of future antiviral agents and SARS-CoV-2 vaccines.

Upregulation of galectin-3 in influenza A virus infection promotes viral RNA synthesis through its association with viral PA protein

BACKGROUND

Influenza is one of the most important viral infections globally. Viral RNA-dependent RNA polymerase (RdRp) consists of the PA, PB1, and PB2 subunits, and the amino acid residues of each subunit are highly conserved among influenza A virus (IAV) strains. Due to the high mutation rate and emergence of drug resistance, new antiviral strategies are needed. Host cell factors are involved in the transcription and replication of influenza virus. Here, we investigated the role of galectin-3, a member of the β-galactoside-binding animal lectin family, in the life cycle of IAV infection in vitro and in mice.

METHODS

We used galectin-3 knockout and wild-type mice and cells to study the intracellular role of galectin-3 in influenza pathogenesis. Body weight and survival time of IAV-infected mice were analyzed, and viral production in mouse macrophages and lung fibroblasts was examined. Overexpression and knockdown of galectin-3 in A549 human lung epithelial cells were exploited to assess viral entry, viral ribonucleoprotein (vRNP) import/export, transcription, replication, virion production, as well as interactions between galectin-3 and viral proteins by immunoblotting, immunofluorescence, co-immunoprecipitation, RT-qPCR, minireplicon, and plaque assays. We also employed recombinant galectin-3 proteins to identify specific step(s) of the viral life cycle that was affected by exogenously added galectin-3 in A549 cells.

RESULTS

Galectin-3 levels were increased in the bronchoalveolar lavage fluid and lungs of IAV-infected mice. There was a positive correlation between galectin-3 levels and viral loads. Notably, galectin-3 knockout mice were resistant to IAV infection. Knockdown of galectin-3 significantly reduced the production of viral proteins and virions in A549 cells. While intracellular galectin-3 did not affect viral entry, it increased vRNP nuclear import, RdRp activity, and viral transcription and replication, which were associated with the interaction of galectin-3 with viral PA subunit. Galectin-3 enhanced the interaction between viral PA and PB1 proteins. Moreover, exogenously added recombinant galectin-3 proteins also enhanced viral adsorption and promoted IAV infection in A549 cells.

CONCLUSION

We demonstrate that galectin-3 enhances viral infection through increases in vRNP nuclear import and RdRp activity, thereby facilitating viral transcription and replication. Our findings also identify galectin-3 as a potential therapeutic target for influenza.

Detection of RNA-dependent RNA polymerase of porcine epidemic diarrhea virus

The RNA synthesis of porcine epidemic diarrhea virus (PEDV) is a sophisticated process performed by a multilingual viral replication complex, together with cellular factors. A key enzyme of this replication complex is RNA-dependent RNA polymerase (RdRp). However, there is limited knowledge about PEDV RdRp. In our present study, a polyclonal antibody against RdRp was prepared by using a prokaryotic expression vector pET-28a-RdRp to study the function of PEDV RdRp and provide a tool to investigate PEDV pathogenesis. In addition, the enzyme activity and half-life of PEDV RdRp were investigated. The result showed that the polyclonal antibody against PEDV RdRp was successfully prepared and was able to be used to detect PEDV RdRp by immunofluorescence and western blotting. Additionally, enzyme activity of PEDV RdRp reached nearly 2 pmol/μg/h and the half-life of PEDV RdRp was 5.47 h.

Quassinoids from Eurycoma longifolia with antiviral activities by inhibiting dengue virus replication

BACKGROUND

Dengue caused by dengue virus (DENV) spreads rapidly around the world. However, there are no worldwide licensed vaccines or specific antivirals to combat DENV infection. Quassinoids are the most characteristic components of Eurycoma longifolia, which have been reported to display a variety of biological activities. However, whether quassinoids exert anti-DENV activities remains unknown.

PURPOSE

To test the quassinoids of E. longifolia for their activity against DENV and to clarify the potential mechanisms.

METHODS

The quassinoids from E. longifolia were isolated by chromatography techniques, and their chemical structures were elucidated by spectroscopic analysis. The anti-DENV activities of quassinoids on baby hamster kidney cells BHK-21 were determined by lactate dehydrogenase (LDH) assay. The synthesis of progeny virus was measured by plaque assay. The expression levels of envelope protein (E) and non-structural protein 1 (NS1) were evaluated by qRT-PCR, Western blot and immunofluorescence assays. Molecular docking was used to screen the potential targets of the most active quassinoid against DENV-2, and surface plasmon resonance analysis was employed to confirm the direct binding between the most active quassinoid and potential target.

RESULTS

Twenty-four quassinoids, including three new quassinoids (1 - 3), were isolated from the ethanol extract of E. longifolia. Quassinoids 4, 5, 9, 11, 12, 15, 16, 17, 19 and 20 significantly reduced the LDH release at the stages of viral binding and entry or intracellular replication. Among them, 19 (6α-hydroxyeurycomalactone, 6α-HEL) exhibited the best anti-DENV-2 activities with an EC₅₀ value of 0.39 ± 0.02 μM. Further experiments suggested that 6α-HEL remarkably inhibited progeny virus synthesis and mRNA and protein expression levels of E and NS1 of DENV-2. Time-of-drug-addition assay suggested that 6α-HEL inhibited intracellular replication of DENV-2 at an early stage. Moreover, 6α-HEL was shown to interact with NS5-RdRp domain at a binding affinity of -8.15 kcal/mol. SPR assay further verified 6α-HEL bound to RdRp protein with an equilibrium dissociation constant of 1.49 × 10^-7 M.

CONCLUSION

Ten quassinoids from E. longifolia showed anti-DENV activities at processes of virus binding and entry or intracellular replication. The most active quassinoid 6α-HEL exerts the anti-DENV-2 activities at intracellular replication stage by directly targeting the NS5-RdRp protein. These results suggest that 6α-HEL could be a promising candidate for the treatment of DENV-2 infection.

Direct observation of backtracking by influenza A and B polymerases upon consecutive incorporation of the nucleoside analog T1106

The antiviral pseudo-base T705 and its de-fluoro analog T1106 mimic adenine or guanine and can be competitively incorporated into nascent RNA by viral RNA-dependent RNA polymerases. Although dispersed, single pseudo-base incorporation is mutagenic, consecutive incorporation causes polymerase stalling and chain termination. Using a template encoding single and then consecutive T1106 incorporation four nucleotides later, we obtained a cryogenic electron microscopy structure of stalled influenza A/H7N9 polymerase. This shows that the entire product-template duplex backtracks by 5 nt, bringing the singly incorporated T1106 to the +1 position, where it forms an unexpected T1106:U wobble base pair. Similar structures show that influenza B polymerase also backtracks after consecutive T1106 incorporation, regardless of whether prior single incorporation has occurred. These results give insight into the unusual mechanism of chain termination by pyrazinecarboxamide base analogs. Consecutive incorporation destabilizes the proximal end of the product-template duplex, promoting irreversible backtracking to a more energetically favorable overall configuration.

Functional Importance of the Hydrophobic Residue 362 in Influenza A PB1 Subunit

PB1, acting as the catalytic subunit of the influenza polymerase, has numerous sequentially and structurally conserved regions. It has been observed that the slight modification of residues in PB1 would greatly affect the polymerase activity and even host adaptation ability. Here, we identified a critical residue, 362M, on the polymerase activity and virus replication. By means of the minireplicon assay, we assured the importance of the hydrophobicity of PB1 362, and the possibility that the size and charge of the side chain might directly interfere with the polymerase function. We also proposed a hydrophobic core between the PA-arch and the PB1 β-hairpin motifs and showed the importance of the core to the polymerase function.

Voltage-dependent anion channel 2 (VDAC2) facilitates the accumulation of rice stripe virus in the vector Laodelphax striatellus

Rice stripe virus (RSV) causes enormous losses in rice production and is transmitted by the small brown planthopper, Laodelphax striatellus, in a persistent-propagative manner. RSV accumulation within the gut lumen of the vector is indispensable for the successful transmission to rice and insects. In this study, we obtained a 1464 bp full-length cDNA of a voltage-dependent anion channel 2 from L. striatellus (LsVDAC2), which encodes a 283 amino acid protein. RSV infection increased the expression of LsVDAC2 in the midguts and ovaries of L. striatellus by 260% and 228%, respectively. Silencing of LsVDAC2 resulted in a 88% reduction of RSV loads at 24 h after RNAi, indicating that LsVDAC2 facilitates RSV accumulation in the vector. Yeast two-hybrid and GST pulldown assays demonstrated that LsVDAC2 interacted with RSV RNA-dependent RNA polymerase, RdRp. Furthermore, experiments in vivo and in vitro showed that LsVDAC2 induced the apoptotic response in RSV-infected insects and tissues. Silencing of LsVDAC2 via RNAi significantly reduced the expression of genes for apoptosis-related caspases 1a and 1c by 62% and 78%, respectively, in RSV-infected vectors. Whether LsVDAC2-induced RSV accumulation is related to RSV RdRp and LsVDAC2-induced cell apoptosis deserves further investigation.

Molecular Evolution of RNA-Dependent RNA Polymerase Region in Norovirus Genogroup I

Norovirus is the leading viral agent of gastroenteritis in humans. RNA-dependent RNA polymerase (RdRp) is essential in the replication of norovirus RNA. Here, we present a comprehensive evolutionary analysis of the norovirus GI RdRp gene. Our results show that the norovirus GI RdRp gene can be divided into three groups, and that the most recent common ancestor was 1484. The overall evolutionary rate of GI RdRp is 1.821 × 10-3 substitutions/site/year. Most of the amino acids of the GI RdRp gene were under negative selection, and only a few positively selected sites were recognized. Amino acid substitutions in the GI RdRp gene accumulated slowly over time. GI.P1, GI.P3 and GI.P6 owned the higher evolutionary rates. GI.P11 and GI.P13 had the faster accumulation rate of amino acid substitutions. GI.P2, GI.P3, GI.P4, GI.P6 and GI.P13 presented a strong linear evolution. These results reveal that the norovirus GI RdRp gene evolves conservatively, and that the molecular evolutionary characteristics of each P-genotype are diverse. Sequencing in RdRp and VP1 of norovirus should be advocated in the surveillance system to explore the effect of RdRp on norovirus activity.

Structural basis of transition from initiation to elongation in de novo viral RNA-dependent RNA polymerases

De novo viral RNA-dependent RNA polymerases (RdRPs) utilize their priming element (PE) to facilitate accurate initiation. Upon transition to elongation, the PE has to retreat from the active site to give room to the template-product RNA duplex. However, PE conformational change upon this transition and the role of PE at elongation both remain elusive. Here, we report crystal structures of RdRP elongation complex (EC) from dengue virus serotype 2 (DENV2), demonstrating a dramatic refolding of PE that allows establishment of interactions with the RNA duplex backbone approved to be essential for EC stability. Enzymology data from both DENV2 and hepatitis C virus (HCV) RdRPs suggest that critical transition of the refolding likely occurs after synthesis of a 4- to 5-nucleotide (nt) product together providing a key basis in understanding viral RdRP transition from initiation to elongation.