The enzymes for genome size increase and maintenance of large (+)RNA viruses

A yeast-based reverse genetics system to generate HCoV-OC43 reporter viruses encoding an eighth subgenomic RNA

Coronaviruses have large, positive-sense single-stranded RNA genomes that challenge conventional strategies for mutagenesis. Yeast genetics has been used to manipulate large viral genomes, including those of herpesviruses and coronaviruses. This method, known as transformation-associated recombination (TAR), involves assembling complete viral genomes from dsDNA copies of viral genome fragments via homologous recombination in Saccharomyces cerevisiae. Here, we report our development of a TAR assembly and mutagenesis system for the endemic, seasonal human coronavirus (HCoV) strain OC43. HCoV-OC43 generally causes mild respiratory symptoms and is classified as a biosafety level 2 agent, making it useful for studying fundamental aspects of coronavirus biology and for comparative studies of more highly pathogenic betacoronaviruses. Following cDNA synthesis from HCoV-OC43 viral RNA, we generated five plasmids encompassing ~7.2 kb portions of the ORF1ab gene, the NS2 to M segment, or the N gene and structured to facilitate reporter gene insertions in the M-to-N intergenic region. Using these plasmids, we completed independent assemblies of yeast centromeric plasmids encoding ORF1ab, NS2a to N, as well as full-length HCoV-OC43 plasmids. A wild-type virus (OC43YA), as well as mClover3-H2B (OC43-mCloYA), mRuby3-H2B (OC43-mRubyYA), and mCardinal (OC43-mCardYA) reporter viruses, were rescued. The OC43-mCloYA reporter virus replicated comparably to an OC43 reference strain and produced the mClover3-H2B protein from a novel subgenomic RNA through insertion of an eighth body transcription regulatory sequence, preventing the need to delete or mutate viral genes. This updated HCoV-OC43 reverse genetics system will contribute to a better understanding of betacoronavirus host-pathogen interactions and can accelerate studies of novel antivirals.

IMPORTANCE

Coronaviruses are ubiquitous pathogens that infect humans resulting in both mild and severe respiratory infections. Human coronavirus strain OC43 (HCoV-OC43) is one of many viruses responsible for common colds and is a useful model of more severe coronavirus infections. In this study, we describe an updated HCoV-OC43 mutagenesis system that uses yeast to capture six DNA fragments of the viral RNA genome and assemble them into full-length genomes in yeast/bacterial plasmids. The design of this system allowed for the rapid assembly and rescue of functional HCoV-OC43 viruses, including fluorescent reporter viruses with expanded genetic capacity. This updated reverse genetics system will enhance our ability to monitor viral replication, through building new reporter viruses, while also enhancing the study of betacoronavirus biology through the generation of mutant HCoV-OC43 viruses.

Giant RNA genomes: Roles of host, translation elongation, genome architecture, and proteome in nidoviruses

Positive-strand RNA viruses of the order Nidovirales have the largest known RNA genomes of vertebrate and invertebrate viruses with 36.7 and 41.1 kb, respectively. The acquisition of a proofreading exoribonuclease (ExoN) by an ancestral nidovirus enabled crossing of the 20 kb barrier. Other factors constraining genome size variations in nidoviruses remain poorly defined. We assemble 76 genome sequences of invertebrate nidoviruses from >500.000 published transcriptome experiments and triple the number of known nidoviruses with >36 kb genomes, including a 64 kb RNA genome. Many of the identified viral lineages acquired putative enzymatic and other protein domains linked to genome size, host phyla, or virus families. The inserted domains may regulate viral replication and virion formation, or modulate infection otherwise. We classify ExoN-encoding nidoviruses into seven groups and four subgroups, according to canonical and noncanonical modes of viral replicase expression by ribosomes and genomic organization (reModes). The most-represented group employing the canonical reMode comprises invertebrate and vertebrate nidoviruses, including coronaviruses. Six groups with noncanonical reModes include invertebrate nidoviruses with 31-to-64 kb genomes. Among them are viruses with segmented genomes and viruses utilizing dual ribosomal frameshifting that we validate experimentally. Moreover, largest polyprotein length and genome size in nidoviruses show reMode- and host phylum-dependent relationships. We hypothesize that the polyprotein length increase in nidoviruses may be limited by the host-inherent translation fidelity, ultimately setting a nidovirus genome size limit. Thus, expansion of ExoN-encoding RNA virus genomes, the vertebrate/invertebrate host division, the control of viral replicase expression, and translation fidelity are interconnected.

The effects of Remdesivir's functional groups on its antiviral potency and resistance against the SARS-CoV-2 polymerase

Remdesivir (RDV, Veklury®) is the first FDA-approved antiviral treatment for COVID-19. It is a nucleotide analogue (NA) carrying a 1'-cyano (1'-CN) group on the ribose and a pseudo-adenine nucleobase whose contributions to the mode of action (MoA) are not clear. Here, we dissect these independent contributions by employing RDV-TP analogues. We show that while the 1'-CN group is directly responsible for transient stalling of the SARS-CoV-2 replication/transcription complex (RTC), the nucleobase plays a role in the strength of this stalling. Conversely, RNA extension assays show that the 1'-CN group plays a role in fidelity and that RDV-TP can be incorporated as a GTP analogue, albeit with lower efficiency. However, a mutagenic effect by the viral polymerase is not ascertained by deep sequencing of viral RNA from cells treated with RDV. We observe that once added to the 3' end of RNA, RDV-MP is sensitive to excision and its 1'-CN group does not impact its nsp14-mediated removal. A >14-fold RDV-resistant SARS-CoV-2 isolate can be selected carrying two mutations in the nsp12 sequence, S759A and A777S. They confer both RDV-TP discrimination over ATP by nsp12 and stalling during RNA synthesis, leaving more time for excision-repair and potentially dampening RDV efficiency. We conclude that RDV presents a multi-faced MoA. It slows down or stalls overall RNA synthesis but is efficiently repaired from the primer strand, whereas once in the template, read-through inhibition adds to this effect. Its efficient incorporation may corrupt proviral RNA, likely disturbing downstream functions in the virus life cycle.

Pioneering role of RNA in the early evolution of life

The catalytic, regulatory and structural properties of RNA, combined with their extraordinary ubiquity in cellular processes, are consistent with the proposal that this molecule played a much more conspicuous role in heredity and metabolism during the early stages of biological evolution. This review explores the pivotal role of RNA in the earliest life forms and its relevance in modern biological systems. It examines current models that study the early evolution of life, providing insights into the primordial RNA world and its legacy in contemporary biology.

Mapping glycoprotein structure reveals Flaviviridae evolutionary history

Viral glycoproteins drive membrane fusion in enveloped viruses and determine host range, tissue tropism and pathogenesis1. Despite their importance, there is a fragmentary understanding of glycoproteins within the Flaviviridae2, a large virus family that include pathogens such as hepatitis C, dengue and Zika viruses, and numerous other human, animal and emergent viruses. For many flaviviruses the glycoproteins have not yet been identified, for others, such as the hepaciviruses, the molecular mechanisms of membrane fusion remain uncharacterized3. Here we combine phylogenetic analyses with protein structure prediction to survey glycoproteins across the entire Flaviviridae. We find class II fusion systems, homologous to the Orthoflavivirus E glycoprotein in most species, including highly divergent jingmenviruses and large genome flaviviruses. However, the E1E2 glycoproteins of the hepaciviruses, pegiviruses and pestiviruses are structurally distinct, may represent a novel class of fusion mechanism, and are strictly associated with infection of vertebrate hosts. By mapping glycoprotein distribution onto the underlying phylogeny, we reveal a complex evolutionary history marked by the capture of bacterial genes and potentially inter-genus recombination. These insights, made possible through protein structure prediction, refine our understanding of viral fusion mechanisms and reveal the events that have shaped the diverse virology and ecology of the Flaviviridae.

A ~40-kb flavi-like virus does not encode a known error-correcting mechanism

It is commonly held that there is a fundamental relationship between genome size and error rate, manifest as a notional "error threshold" that sets an upper limit on genome sizes. The genome sizes of RNA viruses, which have intrinsically high mutation rates due to a lack of mechanisms for error correction, must therefore be small to avoid accumulating an excessive number of deleterious mutations that will ultimately lead to population extinction. The proposed exceptions to this evolutionary rule are RNA viruses from the order Nidovirales (such as coronaviruses) that encode error-correcting exonucleases, enabling them to reach genome lengths greater than 40 kb. The recent discovery of large-genome flavi-like viruses (Flaviviridae), which comprise genomes up to 27 kb in length yet seemingly do not encode exonuclease domains, has led to the proposal that a proofreading mechanism is required to facilitate the expansion of nonsegmented RNA virus genomes above 30 kb. Herein, we describe a ~40 kb flavi-like virus identified in a Haliclona sponge metatranscriptome that does not encode a known exonuclease. Structural analysis revealed that this virus may have instead captured cellular domains associated with nucleic acid metabolism that have not been previously found in RNA viruses. Phylogenetic inference placed this virus as a divergent pesti-like lineage, such that we have provisionally termed it "Maximus pesti-like virus." This virus represents an instance of a flavi-like virus achieving a genome size comparable to that of the Nidovirales and demonstrates that RNA viruses have evolved multiple solutions to overcome the error threshold.

Development of Allosteric Ribozymes for ATP and l-Histidine Based on the R3C Ligase Ribozyme

During the evolution of the RNA, short RNAs are thought to have joined together to form long RNAs, enhancing their function as ribozymes. Previously, the artificial R3C ligase ribozyme (73 nucleotides) was successfully reduced to 46 nucleotides; however, its activity decreased significantly. Therefore, we aimed to develop allosteric ribozymes, whose activities could be regulated by effector compounds, based on the reduced R3C ligase ribozyme (R3C-A). Among the variants prepared by fusing an ATP-binding aptamer RNA with R3C-A, one mutant showed increased ligation activity in an ATP-dependent manner. Melting temperature measurements of the two RNA mutants suggested that the region around the aptamer site was stabilized by the addition of ATP. This resulted in a suitable conformation for the reaction at the ligation site. Another ribozyme was prepared by fusing R3C-A with a l-histidine-binding aptamer RNA, and the ligase activity increased with increasing l-histidine concentrations. Both ATP and l-histidine play prominent roles in current molecular biology and the interaction of RNAs and these molecules could be a key step in the evolution of the world of RNAs. Our results suggest promise in the development of general allosteric ribozymes that are independent of the type of effector molecule and provide important clues to the evolution of the RNA world.

AT-752 targets multiple sites and activities on the Dengue virus replication enzyme NS5

AT-752 is a guanosine analogue prodrug active against dengue virus (DENV). In infected cells, it is metabolized into 2'-methyl-2'-fluoro guanosine 5'-triphosphate (AT-9010) which inhibits RNA synthesis in acting as a RNA chain terminator. Here we show that AT-9010 has several modes of action on DENV full-length NS5. AT-9010 does not inhibit the primer pppApG synthesis step significantly. However, AT-9010 targets two NS5-associated enzyme activities, the RNA 2'-O-MTase and the RNA-dependent RNA polymerase (RdRp) at its RNA elongation step. Crystal structure and RNA methyltransferase (MTase) activities of the DENV 2 MTase domain in complex with AT-9010 at 1.97 Å resolution shows the latter bound to the GTP/RNA-cap binding site, accounting for the observed inhibition of 2'-O but not N7-methylation activity. AT-9010 is discriminated ∼10 to 14-fold against GTP at the NS5 active site of all four DENV1-4 NS5 RdRps, arguing for significant inhibition through viral RNA synthesis termination. In Huh-7 cells, DENV1-4 are equally sensitive to AT-281, the free base of AT-752 (EC50 ≈ 0.50 μM), suggesting broad spectrum antiviral properties of AT-752 against flaviviruses.

RNA-Dependent RNA Polymerase of the Second Human Pegivirus Exhibits a High-Fidelity Feature

The virus-encoded RNA-dependent RNA polymerase (RdRp) is responsible for viral replication, and its fidelity is closely related to viral diversity, pathogenesis, virulence, and fitness. Hepatitis C virus (HCV) and the second human pegivirus (HPgV-2) belong to the family Flaviviridae and share some features, including similar viral genome structure. Unlike HCV, HPgV-2 preserves a highly conserved genome sequence and low intrahost variation. However, the underlying mechanism remains to be elucidated. In this study, we evaluated the fidelity of HPgV-2 and HCV RdRp in an in vitro RNA polymerase reaction system. The results showed higher fidelity of HPgV-2 RdRp than HCV NS5B with respect to the misincorporation rate due to their difference in recognizing nucleoside triphosphate (NTP) substrates. Furthermore, HPgV-2 RdRp showed lower sensitivity than HCV to sofosbuvir, a nucleotide inhibitor against HCV RdRp, which explained the insusceptibility of HPgV-2 to direct-acting antiviral (DAA) therapy against HCV infection. Our results indicate that HPgV-2 could be an excellent model for studying the mechanisms involved in viral polymerase fidelity as well as RNA virus diversity and evolution. IMPORTANCE RNA viruses represent the most important pathogens for humans and animals and exhibit rapid evolution and high adaptive capacity, which is due to the high mutation rates for using the error-prone RNA-dependent RNA polymerase (RdRp) during replication. The fidelity of RdRp is closely associated with viral diversity, fitness, and pathogenesis. Previous studies have shown that the second human pegivirus (HPgV-2) exhibits a highly conserved genome sequence and low intrahost variation, which might be due to the fidelity of HPgV-2 RdRp. In this work, we used a series of in vitro RNA polymerase assays to evaluate the in vitro fidelity of HPgV-2 RdRp and compared it with that of HCV RdRp. The results indicated that HPgV-2 RdRp preserves significantly higher fidelity than HCV RdRp, which might contribute to the conservation of the HPgV-2 genome. The unique feature of HPgV-2 RdRp fidelity provides a new model for investigation of viral RdRp fidelity.

A second type of N7-guanine RNA cap methyltransferase in an unusual locus of a large RNA virus genome

The order Nidovirales is a diverse group of (+)RNA viruses, with a common genome organization and conserved set of replicative and editing enzymes. In particular, RNA methyltransferases play a central role in mRNA stability and immune escape. However, their presence and distribution in different Nidovirales families is not homogeneous. In Coronaviridae, the best characterized family, two distinct methytransferases perform methylation of the N7-guanine and 2'-OH of the RNA-cap to generate a cap-1 structure (m7GpppNm). The genes of both of these enzymes are located in the ORF1b genomic region. While 2'-O-MTases can be identified for most other families based on conservation of both sequence motifs and genetic loci, identification of the N7-guanine methyltransferase has proved more challenging. Recently, we identified a putative N7-MTase domain in the ORF1a region (N7-MT-1a) of certain members of the large genome Tobaniviridae family. Here, we demonstrate that this domain indeed harbors N7-specific methyltransferase activity. We present its structure as the first N7-specific Rossmann-fold (RF) MTase identified for (+)RNA viruses, making it remarkably different from that of the known Coronaviridae ORF1b N7-MTase gene. We discuss the evolutionary implications of such an appearance in this unexpected location in the genome, which introduces a split-off in the classification of Tobaniviridae.

Subtle structural differences of nucleotide analogs may impact SARS-CoV-2 RNA-dependent RNA polymerase and exoribonuclease activity

The rapid spread and public health impact of the novel SARS-CoV-2 variants that cause COVID-19 continue to produce major global impacts and social distress. Several vaccines were developed in record time to prevent and limit the spread of the infection, thus playing a pivotal role in controlling the pandemic. Although the repurposing of available drugs attempts to provide therapies of immediate access against COVID-19, there is still a need for developing specific treatments for this disease. Remdesivir, molnupiravir and Paxlovid remain the only evidence-supported antiviral drugs to treat COVID-19 patients, and only in severe cases. To contribute on the search of potential Covid-19 therapeutic agents, we targeted the viral RNA-dependent RNA polymerase (RdRp) and the exoribonuclease (ExoN) following two strategies. First, we modeled and analyzed nucleoside analogs sofosbuvir, remdesivir, favipiravir, ribavirin, and molnupiravir at three key binding sites on the RdRp-ExoN complex. Second, we curated and virtually screened a database containing 517 nucleotide analogs in the same binding sites. Finally, we characterized key interactions and pharmacophoric features presumably involved in viral replication halting at multiple sites. Our results highlight structural modifications that might lead to more potent SARS-CoV-2 inhibitors against an expansive range of variants and provide a collection of nucleotide analogs useful for screening campaigns.

The SARS-CoV nsp12 Polymerase Active Site Is Tuned for Large-Genome Replication

Positive-strand RNA viruses replicate their genomes using virally encoded RNA-dependent RNA polymerases (RdRP) with a common active-site structure and closure mechanism upon which replication speed and fidelity can evolve to optimize virus fitness. Coronaviruses (CoV) form large multicomponent RNA replication-transcription complexes containing a core RNA synthesis machine made of the nsp12 RdRP protein with one nsp7 and two nsp8 proteins as essential subunits required for activity. We show that assembly of this complex can be accelerated 5-fold by preincubation of nsp12 with nsp8 and further optimized with the use of a novel nsp8L7 heterodimer fusion protein construct. Using rapid kinetics methods, we measure elongation rates of up to 260 nucleotides (nt)/s for the core replicase, a rate that is unusually fast for a viral polymerase. To address the origin of this fast rate, we examined the roles of two CoV-specific residues in the RdRP active site: Ala547, which replaces a conserved glutamate above the bound NTP, and Ser759, which mutates the palm domain GDD sequence to SDD. Our data show that Ala547 allows for a doubling of replication rate, but this comes at a fidelity cost that is mitigated by using a SDD sequence in the palm domain. Our biochemical data suggest that fixation of mutations in polymerase motifs F and C played a key role in nidovirus evolution by tuning replication rate and fidelity to accommodate their large genomes. IMPORTANCE Replicating large genomes represents a challenge for RNA viruses because fast RNA synthesis is needed to escape innate immunity defenses, but faster polymerases are inherently low-fidelity enzymes. Nonetheless, the coronaviruses replicate their ≈30-kb genomes using the core polymerase structure and mechanism common to all positive-strand RNA viruses. The classic explanation for their success is that the large-genome nidoviruses have acquired an exonuclease-based repair system that compensates for the high polymerase mutation rate. In this work, we establish that the nidoviral polymerases themselves also play a key role in maintaining genome integrity via mutations at two key active-site residues that enable very fast replication rates while maintaining typical mutation rates. Our findings further demonstrate the evolutionary plasticity of the core polymerase platform by showing how it has adapted during the expansion from short-genome picornaviruses to long-genome nidoviruses.

The logic of virus evolution

Viruses are obligate intracellular parasites. Despite their dependence on host cells, viruses are evolutionarily autonomous, with their own genomes and evolutionary trajectories locked in arms races with the hosts. Here, we discuss a simple functional logic to explain virus macroevolution that appears to define the course of virus evolution. A small core of virus hallmark genes that are responsible for genome replication apparently descended from primordial replicators, whereas most virus genes, starting with those encoding capsid proteins, were subsequently acquired from hosts. The oldest of these acquisitions antedate the last universal cellular ancestor (LUCA). Host gene capture followed two major routes: convergent recruitment of genes with functions that directly benefit virus reproduction and exaptation when host proteins are repurposed for unique virus functions. These forms of host protein recruitment by viruses result in different levels of similarity between virus and host homologs, with the exapted ones often changing beyond easy recognition.