Evolution of coronavirus frameshifting elements: Competing stem networks explain conservation and variability

The frameshifting RNA element (FSE) in coronaviruses (CoVs) regulates the programmed -1 ribosomal frameshift (-1 PRF) mechanism common to many viruses. The FSE is of particular interest as a promising drug candidate. Its associated pseudoknot or stem loop structure is thought to play a large role in frameshifting and thus viral protein production. To investigate the FSE structural evolution, we use our graph theory-based methods for representing RNA secondary structures in the RNA-As-Graphs (RAG) framework to calculate conformational landscapes of viral FSEs with increasing sequence lengths for representative 10 Alpha and 13 Beta-CoVs. By following length-dependent conformational changes, we show that FSE sequences encode many possible competing stems which in turn favor certain FSE topologies, including a variety of pseudoknots, stem loops, and junctions. We explain alternative competing stems and topological FSE changes by recurring patterns of mutations. At the same time, FSE topology robustness can be understood by shifted stems within different sequence contexts and base pair coevolution. We further propose that the topology changes reflected by length-dependent conformations contribute to tuning the frameshifting efficiency. Our work provides tools to analyze virus sequence/structure correlations, explains how sequence and FSE structure have evolved for CoVs, and provides insights into potential mutations for therapeutic applications against a broad spectrum of CoV FSEs by targeting key sequence/structural transitions.

A potential long-range RNA-RNA interaction in the HIV-1 RNA

It is well-established that viral and cellular mRNAs alike harbour functional long-range intra-molecular RNA-RNA interactions. Despite the biological importance of such interactions, their identification and characterization remain challenging. Here we present a computational method for the identification of certain kinds of long-range intra-molecular RNA-RNA interactions involving the loop nucleotides of a hairpin loop. Using the computational method, we analysed 4272 HIV-1 genomic mRNAs. A potential long-range intra-molecular RNA-RNA interaction within the HIV-1 genomic RNA was identified. The long-range interaction is mediated by a kissing loop structure between two stem-loops of the previously reported SHAPE-based secondary structure of the entire HIV-1 genome. Structural modelling studies were carried out to show that the kissing loop structure not only is sterically feasible, but also contains a conserved RNA structural motif often found in compact RNA pseudoknots. The computational method should be generally applicable to the identification of potential long-range intra-molecular RNA-RNA interactions in any viral or cellular mRNA sequence.Communicated by Ramaswamy H. Sarma.

HIV co-opts a cellular antiviral mechanism, activation of stress kinase PKR by its RNA, to enable splicing of rev/tat mRNA

BACKGROUND

Activation of RNA-dependent stress kinase PKR, especially by viral double-stranded RNA, induces eukaryotic initiation factor 2 α-chain (eIF2α) phosphorylation, attenuating thereby translation. We report that this RNA-mediated negative control mechanism, considered a cornerstone of the cell's antiviral response, positively regulates splicing of a viral mRNA.

RESULTS

Excision of the large human immunodeficiency virus (HIV) rev/tat intron depends strictly on activation of PKR by the viral RNA and on eIF2α phosphorylation. Rev/tat mRNA splicing was blocked by viral PKR antagonists Vaccinia E3L and Ebola VP35, as well as by a trans-dominant negative mutant of PKR, yet enhanced by overexpressing PKR. Expression of non-phosphorylatable mutant eIF2αS51A, but not of wild type eIF2α, abrogated efficient splicing of rev/tat mRNA. By contrast, expression of eIF2αS51D, a phosphomimetic mutant of eIF2α, left rev/tat mRNA splicing intact. Unlike eIF2αS51A, eIF2αS51D does not inhibit eIF2α phosphorylation by activated PKR. All HIV mRNA species contain terminal trans-activation response (TAR) stem-loop sequences that potentially could activate PKR, yet even upon TAR deletion, HIV mRNA production remained sensitive to inhibitors of PKR activation. Bioinformatic and mutational analyses revealed a compact RNA pseudoknot upstream of 3'-terminal TAR that promotes splicing by activating PKR. Supporting its essential role in control of splicing, this pseudoknot is conserved among diverse HIV and nonhuman primate SIVcpz isolates. The pseudoknot and 3'-terminal TAR collaborate in mediating PKR-regulated splicing of rev/tat intron, the pseudoknot being dominant.

CONCLUSIONS

Our results on HIV provide the first example of a virus co-opting activation of PKR by its RNA, a cellular antiviral mechanism, to promote splicing. They raise the question whether other viruses may use local activation of host kinase PKR through RNA elements within their genome to achieve efficient splicing of their mRNA. Our experiments reveal an indispensable role for eIF2α phosphorylation in HIV rev/tat mRNA splicing that accounts for the need for PKR activation.

Investigating the correlation between Xrn1-resistant RNAs and frameshifter pseudoknots

Xrn1-resistant RNA structures are multifunctional elements employed by an increasing number of RNA viruses. One of such elements is the coremin motif, discovered in plant virus RNAs, of which the structure has been hypothesized to form a yet unelucidated pseudoknot. Recently, the coremin motif was shown to be capable of stalling not only Xrn1, but scanning ribosomes as well. Following that observation, in this study we demonstrate that the coremin motif can promote -1 ribosomal frameshifting, similar to better-characterized viral frameshifting pseudoknots. Since this function was lost in concert with substitutions that were known to disturb Xrn1-resistance, we developed a frameshifting screen for finding novel Xrn1-resistant RNAs by randomizing parts of the coremin motif. This yielded new insights into the coremin motif structure, as Xrn1-resistant variations were identified that more clearly indicate a pseudoknot interaction. In addition, we show that the Xrn1-resistant RNA of Zika virus promotes frameshifting as well, while known -1 programmed ribosomal frameshifting pseudoknots do not stall Xrn1, suggesting that promoting frameshifting is a universal characteristic of Xrn1-resistant RNAs, but that Xrn1-resistance requires more than just a frameshifting pseudoknot.

VfoldMCPX: predicting multistrand RNA complexes

Multistrand RNA complexes play a critical role in RNA-related biological processes. The understanding of RNA functions and the rational design of RNA nanostructures require accurate prediction of the structure and folding stability of the complexes, including those containing pseudoknots. Here, we present VfoldMCPX, a new model for predicting two-dimensional (2D) structures and folding stabilities of multistrand RNA complexes. Based on a partition function-based algorithm combined with physical loop free energy parameters, the VfoldMCPX model predicts not only the native structure but also the folding stability of the complex. An important advantage of the model is the ability to treat pseudoknotted structures. Extensive tests on structure predictions show the VfoldMCPX model provides improved accuracy for multistranded RNA complexes, especially for RNA complexes with three or more strands and/or containing pseudoknots. We have developed a freely accessible VfoldMCPX web server at http://rna.physics.missouri.edu/vfoldMCPX2.

Restriction of SARS-CoV-2 replication by targeting programmed -1 ribosomal frameshifting

A large variety of RNA viruses, including the novel coronavirus SARS-CoV-2, contain specific RNA structures that promote programmed ribosomal frameshifting (PRF) to regulate viral gene expression. From a high-throughput compound screen, we identified a PRF inhibitor for SARS-CoV-2 and found that it substantially impeded viral replication in cultured cells. Interestingly, the compound could target not only SARS-CoV-2 but also other coronaviruses that use similar RNA structures to promote frameshifting. These results suggest targeting PRF is a plausible, effective, and broad-spectrum antiviral strategy for SARS-CoV-2 and other coronaviruses.

Translation of open reading frame 1b (ORF1b) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) requires a programmed −1 ribosomal frameshift (−1 PRF) promoted by an RNA pseudoknot. The extent to which SARS-CoV-2 replication may be sensitive to changes in −1 PRF efficiency is currently unknown. Through an unbiased, reporter-based high-throughput compound screen, we identified merafloxacin, a fluoroquinolone antibacterial, as a −1 PRF inhibitor for SARS-CoV-2. Frameshift inhibition by merafloxacin is robust to mutations within the pseudoknot region and is similarly effective on −1 PRF of other betacoronaviruses. Consistent with the essential role of −1 PRF in viral gene expression, merafloxacin impedes SARS-CoV-2 replication in Vero E6 cells, thereby providing proof-of-principle for targeting −1 PRF as a plausible and effective antiviral strategy for SARS-CoV-2 and other coronaviruses.

Xrn1-resistant RNA structures are well-conserved within the genus flavivirus

Subgenomic RNAs are produced by several RNA viruses through incomplete degradation of their genomic RNA by the exoribonuclease Xrn1, and have been shown to be essential for viral growth and pathogenicity. Within the flavivirus genus of the Flaviviridae family, two distinct classes of Xrn1-resistant RNA motifs have been proposed; one for mosquito-borne and insect-specific flaviviruses, and one for tick-borne flaviviruses and no-known-vector flaviviruses. We investigated tick-borne and no-known-vector flavivirus Xrn1-resistant RNA motifs through systematic in vitro mutational analysis and showed that both classes actually possess very similar structural configurations, including a double pseudoknot and a base-triple at identical, conserved locations. For the no-known-vector flavivirus Modoc virus, we show that in vivo generation of subgenomic flaviviral RNA was affected by mutations targeted at nucleotides involved in the structural features of flaviviral Xrn1-resistant RNA motifs that were defined in this work. Our results suggest that throughout the genus flavivirus Xrn1-resistant RNA motifs adopt the same topologically conserved structure.

All genera of Flaviviridae host a conserved Xrn1-resistant RNA motif

After infection by flaviviruses like Zika and West Nile virus, eukaryotic hosts employ the well-conserved endoribonuclease Xrn1 to degrade the viral genomic RNA. Within the 3' untranslated regions, this enzyme encounters intricate Xrn1-resistant structures. This results in the accumulation of subgenomic flaviviral RNAs, an event that improves viral growth and aggravates viral pathogenicity. Xrn1-resistant RNAs have been established throughout the flaviviral genus, but not yet throughout the entire Flaviviridae family. In this work, we use previously determined characteristics of these structures to identify homologous sequences in many members of the genera pegivirus, hepacivirus and pestivirus. We used structural alignment and mutational analyses to establish that these sequences indeed represent Xrn1-resistant RNA and that they employ the general features of the flaviviral xrRNAs, consisting of a double pseudoknot formed by five base-paired regions stitched together by a crucial triple base interaction. Furthermore, we demonstrate that the pestivirus Bungowannah virus produces subgenomic RNA in vivo. Altogether, these results indicate that viruses make use of a universal Xrn1-resistant RNA throughout the Flaviviridae family.

Salt-Dependent RNA Pseudoknot Stability: Effect of Spatial Confinement

Macromolecules, such as RNAs, reside in crowded cell environments, which could strongly affect the folded structures and stability of RNAs. The emergence of RNA-driven phase separation in biology further stresses the potential functional roles of molecular crowding. In this work, we employed the coarse-grained model that was previously developed by us to predict 3D structures and stability of the mouse mammary tumor virus (MMTV) pseudoknot under different spatial confinements over a wide range of salt concentrations. The results show that spatial confinements can not only enhance the compactness and stability of MMTV pseudoknot structures but also weaken the dependence of the RNA structure compactness and stability on salt concentration. Based on our microscopic analyses, we found that the effect of spatial confinement on the salt-dependent RNA pseudoknot stability mainly comes through the spatial suppression of extended conformations, which are prevalent in the partially/fully unfolded states, especially at low ion concentrations. Furthermore, our comprehensive analyses revealed that the thermally unfolding pathway of the pseudoknot can be significantly modulated by spatial confinements, since the intermediate states with more extended conformations would loss favor when spatial confinements are introduced.

The crystal structure of a Polerovirus exoribonuclease-resistant RNA shows how diverse sequences are integrated into a conserved fold

Exoribonuclease-resistant RNAs (xrRNAs) are discrete elements that block the progression of 5' to 3' exoribonucleases using specifically folded RNA structures. A recently discovered class of xrRNA is widespread in several genera of plant-infecting viruses, within both noncoding and protein-coding subgenomic RNAs. The structure of one such xrRNA from a dianthovirus revealed three-dimensional details of the resistant fold but did not answer all questions regarding the conservation and diversity of this xrRNA class. Here, we present the crystal structure of a representative polerovirus xrRNA that contains sequence elements that diverge from the previously solved structure. This new structure rationalizes previously unexplained sequence conservation patterns and shows interactions not present in the first structure. Together, the structures of these xrRNAs from dianthovirus and polerovirus genera support the idea that these plant virus xrRNAs fold through a defined pathway that includes a programmed intermediate conformation. This work deepens our knowledge of the structure-function relationship of xrRNAs and shows how evolution can craft similar RNA folds from divergent sequences.

The crystal structure of a Polerovirus exoribonuclease-resistant RNA shows how diverse sequences are integrated into a conserved fold

Exoribonuclease-resistant RNAs (xrRNAs) are discrete elements that block the progression of 5' to 3' exoribonucleases using specifically folded RNA structures. A recently discovered class of xrRNA is widespread in several genera of plant-infecting viruses, within both noncoding and protein-coding subgenomic RNAs. The structure of one such xrRNA from a dianthovirus revealed three-dimensional details of the resistant fold but did not answer all questions regarding the conservation and diversity of this xrRNA class. Here, we present the crystal structure of a representative polerovirus xrRNA that contains sequence elements that diverge from the previously solved structure. This new structure rationalizes previously unexplained sequence conservation patterns and shows interactions not present in the first structure. Together, the structures of these xrRNAs from dianthovirus and polerovirus genera support the idea that these plant virus xrRNAs fold through a defined pathway that includes a programmed intermediate conformation. This work deepens our knowledge of the structure-function relationship of xrRNAs and shows how evolution can craft similar RNA folds from divergent sequences.

Hierarchical Conformational Dynamics Confers Thermal Adaptability to preQ1 RNA Riboswitches

Single-stranded noncoding regulatory RNAs, as exemplified by bacterial riboswitches, are highly dynamic. The conformational dynamics allow the riboswitch to reach maximum switching efficiency under appropriate conditions. Here we characterize the conformational dynamics of preQ₁ riboswitches from mesophilic and thermophilic bacterial species at various temperatures. With the integrative use of small-angle X-ray scattering, NMR, and molecular dynamics simulations, we model the ensemble-structures of the preQ₁ riboswitch aptamers without or with a ligand bound. We show that the preQ₁ riboswitch is sufficiently dynamic and fluctuating among multiple folding intermediates only near the physiological temperature of the microorganism. The hierarchical folding dynamics of the RNA involves the docking of 3'-tail to form a second RNA helix and the helical stacking to form an H-type pseudoknot structure. Further, we show that RNA secondary and tertiary dynamics can be modulated by temperature and by the length of an internal loop. The coupled equilibria between RNA folding intermediates are essential for preQ₁ binding, and a four-state exchange model can account for the change of ligand-triggered switching efficiency with temperature. Together, we have established a relationship between the hierarchical dynamics and riboswitch function, and illustrated how the RNA adapts to high temperature.

An RNA pseudoknot stimulates HTLV-1 pro-pol programmed -1 ribosomal frameshifting

Programmed -1 ribosomal frameshifts (-1 PRFs) are commonly used by viruses to regulate their enzymatic and structural protein levels. Human T-cell leukemia virus type 1 (HTLV-1) is a carcinogenic retrovirus that uses two independent -1 PRFs to express viral enzymes critical to establishing new HTLV-1 infections. How the cis-acting RNA elements in this viral transcript function to induce frameshifting is unknown. The objective of this work was to conclusively define the 3' boundary of and the RNA elements within the HTLV-1 pro-pol frameshift site. We hypothesized that the frameshift site structure was a pseudoknot and that its 3' boundary would be defined by the pseudoknot's 3' end. To test these hypotheses, the in vitro frameshift efficiencies of three HTLV-1 pro-pol frameshift sites with different 3' boundaries were quantified. The results indicated that nucleotides included in the longest construct were essential to highly efficient frameshift stimulation. Interestingly, only this construct could form the putative frameshift site pseudoknot. Next, the secondary structure of this frameshift site was determined. The dominant structure was an H-type pseudoknot which, together with the slippery sequence, stimulated frameshifting to 19.4(±0.3)%. The pseudoknot's critical role in frameshift stimulation was directly revealed by examining the impact of structural changes on HTLV-1 pro-pol -1 PRF. As predicted, mutations that occluded pseudoknot formation drastically reduced the frameshift efficiency. These results are significant because they demonstrate that a pseudoknot is important to HTLV-1 pro-pol -1 PRF and define the frameshift site's 3' boundary.

Large-scale frequent stem pattern mining in RNA families

Functionally similar non-coding RNAs are expected to be similar in certain regions of their secondary structures. These similar regions are called common structure motifs, and are structurally conserved throughout evolution to maintain their functional roles. Common structure motif identification is one of the critical tasks in RNA secondary structure analysis. Nevertheless, current approaches suffer several limitations, and/or do not scale with both structure size and the number of input secondary structures. In this work, we present a method to transform the conserved base pair stems into transaction items and apply frequent itemset mining to identify common structure motifs existing in a majority of input structures. Our experimental results on telomerase and ribosomal RNA secondary structures report frequent stem patterns that are of biological significance. Moreover, the algorithms utilized in our method are scalable and frequent stem patterns can be identified efficiently among many large structures.

An Ancient Pseudoknot in TNF-α Pre-mRNA Activates PKR, Inducing eIF2α Phosphorylation that Potently Enhances Splicing

Tumor necrosis factor alpha (TNF-α) is expressed promptly during inflammatory responses. Efficient TNF-α mRNA splicing is achieved through a 3' UTR element that activates RNA-dependent eIF2α protein kinase (PKR). The TNF-α RNA activator, we show, folds into a pseudoknot conserved from teleost fish to humans, critical for PKR activation and mRNA splicing. The pseudoknot constrains the RNA into two double-helical stacks having parallel axes, permitting facile PKR dimerization and trans-autophosphorylation needed for kinase activation. Mutations show that the PKR activator potently enhances splicing without inhibiting translation. eIF2α phosphorylation represses translation and is essential for coping with cellular stress, yet PKR-enabled TNF mRNA splicing depends strictly on eIF2α phosphorylation. Indeed, eIF2α phosphorylation at Serine51 is necessary and sufficient to achieve highly efficient splicing, extending its role from negative control of translation to positive control of splicing. This mechanism, operational in human peripheral blood mononuclear cells (PBMCs), links stress signaling to protective immunity through TNF mRNA splicing rendered efficient upon eIF2α phosphorylation.

Viral molecular mimicry circumvents abortive infection and suppresses bacterial suicide to make hosts permissive for replication

The global interplay between bacteria and bacteriophages has generated many macromolecules useful in biotechnology, through the co-evolutionary see-saw of bacterial defense and viral counter-attack measures. Bacteria can protect themselves using abortive infection systems, which induce altruistic suicide in an infected cell and therefore protect the clonal population at the expense of the infected individual. Our recent paper describes how bacteriophage ΦTE successfully subverted the activity of a plasmid-borne abortive infection system. ΦTE evolved mimics of the small RNA antitoxin that naturally inhibits the active toxin component of this anti-viral mechanism. These mutant phages further manipulated the behavior of the host population, through transduction of the plasmid encoding the abortive infection system. Transductants thereby became enslaved by the abortive infection system, committing suicide in response to infection by the original phage population. In effect, the new host was infected by an “addictive altruism,” to the advantage of the resistant bacteriophage.

RNA structural constraints in the evolution of the influenza A virus genome NP segment

Conserved RNA secondary structures were predicted in the nucleoprotein (NP) segment of the influenza A virus genome using comparative sequence and structure analysis. A number of structural elements exhibiting nucleotide covariations were identified over the whole segment length, including protein-coding regions. Calculations of mutual information values at the paired nucleotide positions demonstrate that these structures impose considerable constraints on the virus genome evolution. Functional importance of a pseudoknot structure, predicted in the NP packaging signal region, was confirmed by plaque assays of the mutant viruses with disrupted structure and those with restored folding using compensatory substitutions. Possible functions of the conserved RNA folding patterns in the influenza A virus genome are discussed.

Highly conserved RNA pseudoknots at the Gag-Pol junction of HIV-1 suggest a novel mechanism of -1 ribosomal frameshifting

-1 programmed ribosomal frameshifting (PRF) is utilized by many viruses to synthesize their enzymatic (Pol) and structural (Gag) proteins at a defined ratio. For efficient -1 PRF, two cis-acting elements are required: a heptanucleotide frameshift site and a downstream stimulator such as a pseudoknot. We have analyzed the gag-pol junction sequences from 4254 HIV-1 strains. Approximately ninety-five percent of the sequences can form four pseudoknots PK1-PK4 (∼ 97% contain PK1, PK3, and PK4), covering ∼ 72 nt including the frameshift site. Some pseudoknots are mutually excluded due to sequence overlap. PK1 and PK3 arrange tandemly. Their stems form a quasi-continuous helix of ∼ 22 bp. We propose a novel mechanism for possible roles of these pseudoknots. Multiple alternative structures may exist at the gag-pol junction. In most strains, the PK1-PK3 tandem pseudoknots may dominate the structurally heterogeneous pool of RNA due to their greater overall stability. The tandem pseudoknots may function as a breaking system to slow down the ribosome. The ribosome unwinds PK1 and stem 1 of PK3 before it can reach the frameshift site. Then, PK4 can form rapidly because the intact stem 2 of PK3 makes up a large part of the stem 1 of PK4. The newly formed PK4 jams the entrance of the mRNA tunnel. The process then proceeds as in a typical case of -1 PRF. This mechanism incorporates several exquisite new features while still being consistent with the current paradigm of pseudoknot-dependent -1 PRF.