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

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

Elaborated pseudoknots that stimulate -1 programmed ribosomal frameshifting or stop codon readthrough in RNA viruses

Pseudoknots assume various functions including stimulation of -1 programmed ribosomal frameshifting (PRF) or stop codon readthrough (SCR) in RNA viruses. These pseudoknots vary greatly in sizes and structural complexities. Recent biochemical and structural studies confirm the three-stemmed pseudoknots as the -1 PRF stimulators in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and related coronaviruses. We reexamined previously reported -1 PRF or SCR stimulating pseudoknots, especially those containing a relatively long connecting loop between the two pseudoknot-forming stems, for their ability to form elaborated structures. Many potential elaborated pseudoknots were identified that contain one or more of the following extra structural elements: stem-loop, embedded pseudoknot, kissing hairpins, and additional loop-loop interactions. The elaborated pseudoknots are found in several different virus families that utilize either the -1 PRF or SCR recoding mechanisms. Model-building studies were performed to not only establish the structural feasibility of the elaborated pseudoknots but also reveal potential additional structural features that cannot be readily inferred from the predicted secondary structures. Some of the structures, such as embedded double pseudoknots and compact loop-loop pseudoknots mediated by the previously established common pseudoknot motif-1 (CPK-1), represent the first of its kind in the literatures. By advancing discovery of new functional RNA structures, we significantly expand the repertoire of known elaborated pseudoknots that could potentially play a role in -1 PRF and SCR regulation. These results contribute to a better understanding of RNA structures in general, facilitating the design of engineering RNA molecules with certain desired functions.Communicated by Ramaswamy H. Sarma.

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.

The In Vivo and In Vitro Architecture of the Hepatitis C Virus RNA Genome Uncovers Functional RNA Secondary and Tertiary Structures

Hepatitis C virus (HCV) is a positive-strand RNA virus that remains one of the main contributors to chronic liver disease worldwide. Studies over the last 30 years have demonstrated that HCV contains a highly structured RNA genome and many of these structures play essential roles in the HCV life cycle. Despite the importance of riboregulation in this virus, most of the HCV RNA genome remains functionally unstudied. Here, we report a complete secondary structure map of the HCV RNA genome in vivo, which was studied in parallel with the secondary structure of the same RNA obtained in vitro. Our results show that HCV is folded extensively in the cellular context. By performing comprehensive structural analyses on both in vivo data and in vitro data, we identify compact and conserved secondary and tertiary structures throughout the genome. Genetic and evolutionary functional analyses demonstrate that many of these elements play important roles in the virus life cycle. In addition to providing a comprehensive map of RNA structures and riboregulatory elements in HCV, this work provides a resource for future studies aimed at identifying therapeutic targets and conducting further mechanistic studies on this important human pathogen. IMPORTANCE HCV has one of the most highly structured RNA genomes studied to date, and it is a valuable model system for studying the role of RNA structure in protein-coding genes. While previous studies have identified individual cases of regulatory RNA structures within the HCV genome, the full-length structure of the HCV genome has not been determined in vivo. Here, we present the complete secondary structure map of HCV determined both in cells and from corresponding transcripts generated in vitro. In addition to providing a comprehensive atlas of functional secondary structural elements throughout the genomic RNA, we identified a novel set of tertiary interactions and demonstrated their functional importance. In terms of broader implications, the pipeline developed in this study can be applied to other long RNAs, such as long noncoding RNAs. In addition, the RNA structural motifs characterized in this study broaden the repertoire of known riboregulatory elements.

The translational landscape of SARS-CoV-2 and infected cells

SARS-CoV-2 utilizes a number of strategies to modulate viral and host mRNA translation. Here, we used ribosome profiling in SARS-CoV-2 infected model cell lines and primary airway cells grown at the air-liquid interface to gain a deeper understanding of the translationally regulated events in response to virus replication. We find that SARS-CoV-2 mRNAs dominate the cellular mRNA pool but are not more efficiently translated than cellular mRNAs. SARS-CoV-2 utilized a highly efficient ribosomal frameshifting strategy in comparison to HIV-1, suggesting utilization of distinct structural elements. In the highly permissive cell models, although SARS-CoV-2 infection induced the transcriptional upregulation of numerous chemokines, cytokines and interferon stimulated genes, many of these mRNAs were not translated efficiently. Impact of SARS-CoV-2 on host mRNA translation was more subtle in primary cells, with marked transcriptional and translational upregulation of inflammatory and innate immune responses and downregulation of processes involved in ciliated cell function. Together, these data reveal the key role of mRNA translation in SARS-CoV-2 replication and highlight unique mechanisms for therapeutic development.

Specific reverse transcriptase slippage at the HIV ribosomal frameshift sequence: potential implications for modulation of GagPol synthesis

Synthesis of HIV GagPol involves a proportion of ribosomes translating a U6A shift site at the distal end of the gag gene performing a programmed -1 ribosomal frameshift event to enter the overlapping pol gene. In vitro studies here show that at the same shift motif HIV reverse transcriptase generates -1 and +1 indels with their ratio being sensitive to the relative concentration ratio of dNTPs specified by the RNA template slippage-prone sequence and its 5' adjacent base. The GGG sequence 3' adjacent to the U6A shift/slippage site, which is important for ribosomal frameshifting, is shown here to limit reverse transcriptase base substitution and indel 'errors' in the run of A's in the product. The indels characterized here have either 1 more or less A, than the corresponding number of template U's. cDNA with 5 A's may yield novel Gag product(s), while cDNA with an extra base, 7 A's, may only be a minor contributor to GagPol polyprotein. Synthesis of a proportion of non-ribosomal frameshift derived GagPol would be relevant in efforts to identify therapeutically useful compounds that perturb the ratio of GagPol to Gag, and pertinent to the extent in which specific polymerase slippage is utilized in gene expression.

Conformational dynamics of the frameshift stimulatory structure in HIV-1

Programmed ribosomal frameshifting (PRF) in HIV-1 is thought to be stimulated by a hairpin in the mRNA, although a pseudoknot-like triplex has also been proposed. Because the conformational dynamics of the stimulatory structure under tension applied by the ribosomal helicase during translation may play an important role in PRF, we used optical tweezers to apply tension to the HIV stimulatory structure and monitor its unfolding and refolding dynamics. The folding and unfolding kinetics and energy landscape of the hairpin were measured by ramping the force on the hairpin up and down, providing a detailed biophysical characterization. Unexpectedly, whereas unfolding reflected the simple two-state behavior typical of many hairpins, refolding was more complex, displaying significant heterogeneity. Evidence was found for multiple refolding pathways as well as previously unsuspected, partially folded intermediates. Measuring a variant mRNA containing only the sequence required to form the proposed triplex, it behaved largely in the same way. Nonetheless, very rarely, high-force unfolding events characteristic of pseudoknot-like structures were observed. The rare occurrence of the triplex suggests that the hairpin is the functional stimulatory structure. The unusual heterogeneity of the hairpin dynamics under tension suggests a possible functional role in PRF similar to the dynamics of other stimulatory structures.

Expression of a MORN repeat protein from Euplotes octocarinatus requires a +1 programmed ribosomal frameshifting

Analysis of transcriptome revealed that a membrane occupation and recognition nexus (MORN) repeat protein-encoding gene of Euplotes octocarinatus (Eo-morn-9-31) was a candidate for programmed +1 ribosomal frameshifting (+1 PRF). In this study, a dual-luciferase assay was performed to detect its expression. The result showed that the MORN repeat protein (Eo-MORN-9-31) could be produced by the +1 PRF event during the process of translation in yeast and the frameshifting efficiency was about 4-5%. We further confirmed its reality by western blot and mass spectrometry. This study provided experimental evidence indicating that the expression of the Eo-MORN-9-31 of E. octocarinatus required the +1 PRF.

Changes in the Plasticity of HIV-1 Nef RNA during the Evolution of the North American Epidemic

Because of a high mutation rate, HIV exists as a viral swarm of many sequence variants evolving under various selective pressures from the human immune system. Although the Nef gene codes for the most immunogenic of HIV accessory proteins, which alone makes it of great interest to HIV research, it also encodes an RNA structure, whose contribution to HIV virulence has been largely unexplored. Nef RNA helps HIV escape RNA interference (RNAi) through nucleotide changes and alternative folding. This study examines Historic and Modern Datasets of patient HIV-1 Nef sequences during the evolution of the North American epidemic for local changes in RNA plasticity. By definition, RNA plasticity refers to an RNA molecule’s ability to take alternative folds (i.e., alternative conformations). Our most important finding is that an evolutionarily conserved region of the HIV-1 Nef gene, which we denote by R2, recently underwent a statistically significant increase in its RNA plasticity. Thus, our results indicate that Modern Nef R2 typically accommodates an alternative fold more readily than Historic Nef R2. Moreover, the increase in RNA plasticity resides mostly in synonymous nucleotide changes, which cannot be a response to selective pressures on the Nef protein. R2 may therefore be of interest in the development of antiviral RNAi therapies.

A review on architecture of the gag-pol ribosomal frameshifting RNA in human immunodeficiency virus: a variability survey of virus genotypes

Programmed '-1' ribosomal frameshifting is necessary for expressing the pol gene overlapped from a gag of human immunodeficiency virus. A viral RNA structure that requires base pairing across the overlapping sequence region suggests a mechanism of regulating ribosome and helicase traffic during expression. To get precise roles of an element around the frameshift site, a review on architecture of the frameshifting RNA is performed in combination of reported information with augments of a representative set of 19 viral samples. In spite of a different length for the viral RNAs, a canonical comparison on the element sequence allocation is performed for viewing variability associations between virus genotypes. Additionally, recent and historical insights recognized in frameshifting regulation are looked back as for indel and single nucleotide polymorphism of RNA. As specially noted, structural changes at a frameshift site, the spacer sequence, and a three-helix junction element, as well as two Watson-Crick base pairs near a bulge and a C-G pair close a loop, are the most vital strategies for the virus frameshifting regulations. All of structural changes, which are dependent upon specific sequence variations, facilitate an elucidation about the RNA element conformation-dependent mechanism for frameshifting. These facts on disrupting base pair interactions also allow solving the problem of competition between ribosome and helicase on a same RNA template, common to single-stranded RNA viruses. In a broad perspective, each new insight of frameshifting regulation in the competition systems introduced by the RNA element construct changes will offer a compelling target for antiviral therapy.

Stability of HIV Frameshift Site RNA Correlates with Frameshift Efficiency and Decreased Virus Infectivity

Human immunodeficiency virus (HIV) replication is strongly dependent upon a programmed ribosomal frameshift. Here we investigate the relationships between the thermodynamic stability of the HIV type 1 (HIV-1) RNA frameshift site stem-loop, frameshift efficiency, and infectivity, using pseudotyped HIV-1 and HEK293T cells. The data reveal a strong correlation between frameshift efficiency and local, but not overall, RNA thermodynamic stability. Mutations that modestly increase the local stability of the frameshift site RNA stem-loop structure increase frameshift efficiency 2-fold to 3-fold in cells. Thus, frameshift efficiency is determined by the strength of the thermodynamic barrier encountered by the ribosome. These data agree with previous in vitro measurements, suggesting that there are no virus- or host-specific factors that modulate frameshifting. The data also indicate that there are no sequence-specific requirements for the frameshift site stem-loop. A linear correlation between Gag-polymerase (Gag-Pol) levels in cells and levels in virions supports the idea of a stochastic virion assembly mechanism. We further demonstrate that the surrounding genomic RNA secondary structure influences frameshift efficiency and that a mutation that commonly arises in response to protease inhibitor therapy creates a functional but inefficient secondary slippery site. Finally, HIV-1 mutants with enhanced frameshift efficiencies are significantly less infectious, suggesting that compounds that increase frameshift efficiency by as little as 2-fold may be effective at suppressing HIV-1 replication.

IMPORTANCE

HIV, like many retroviruses, utilizes a -1 programmed ribosomal frameshift to generate viral enzymes in the form of a Gag-Pol polyprotein precursor. Thus, frameshifting is essential for viral replication. Here, we utilized a panel of mutant HIV strains to demonstrate that in cells, frameshifting efficiency is correlated with the stability of the local thermodynamic barrier to ribosomal translocation. Increasing the stability of the frameshift site RNA increases the frameshift efficiency 2-fold to 3-fold. Mutant viruses with increased frameshift efficiencies have significantly reduced infectivity. These data suggest that this effect might be exploited in the development of novel antiviral strategies.

The Complementarity of the Loop to the Stem in DNA Pseudoknots Gives Rise to Local TAT Base-Triplets

Pseudoknots belong to an RNA structural motif that has significant roles in the biological function of RNA. An example is ribosomal frameshifting; in this mechanism, the formation of a local triplex changes the reading frame that allows for differences in the translation of mRNAs. In this work, we have used a combination of temperature-dependent UV spectroscopy and differential scanning calorimetry (DSC) to determine the unfolding thermodynamics of a set of DNA pseudoknots with the following sequence: d(TCTCTTnAAAAAAAAGAGAT5TTTTTTT), where "Tn" is a thymine loop with n=5 (PsK-5), 7 (PsK-7), 9 (PsK-9), or 11 (PsK-11). All four oligonucleotides form intramolecular pseudoknots, and the increase in the length of this loop yielded more stable pseudoknots due to higher transition temperatures and higher unfolding enthalpies. This indicates formation of one and three TAT/TAT stacks in PsK-9 and PsK-11, respectively. We have flipped one AT for a TA base pair in the core stem of these pseudoknots, preventing in this way the formation of these base-triplet stacks. The DSC curves of these pseudoknots yielded lower unfolding enthalpies, confirming the formation of a local triplex in PsK-9 and PsK-11. Furthermore, we have investigated the reaction of PsK-5 and PsK-9 with their partially complementary strands: directly by isothermal titration calorimetry and indirectly by creating a Hess cycle with the DSC data. Relative to the PsK-5 reaction, PsK-9 reacts with its complementary strand with less favorable free energy and enthalpy contributions; this indicates PsK-9 is more stable and more compact due to the formation of a local triplex.

The role of sequence in altering the unfolding pathway of an RNA pseudoknot: a steered molecular dynamics study

This work highlights a sequence dependent unfolding pathway of an RNA pseudoknot under force-induced pulling conditions.

Characterization of a Novel Megabirnavirus from Sclerotinia sclerotiorum Reveals Horizontal Gene Transfer from Single-Stranded RNA Virus to Double-Stranded RNA Virus

Mycoviruses have been detected in all major groups of filamentous fungi, and their study represents an important branch of virology. Here, we characterized a novel double-stranded RNA (dsRNA) mycovirus, Sclerotinia sclerotiorum megabirnavirus 1 (SsMBV1), in an apparently hypovirulent strain (SX466) of Sclerotinia sclerotiorum. Two similarly sized dsRNA segments (L1- and L2-dsRNA), the genome of SsMBV1, are packaged in rigid spherical particles purified from strain SX466. The full-length cDNA sequence of L1-dsRNA/SsMBV1 comprises two large open reading frames (ORF1 and ORF2), which encode a putative coat protein and an RNA-dependent RNA polymerase (RdRp), respectively. Phylogenetic analysis of the RdRp domain clearly indicates that SsMBV1 is related to Rosellinia necatrix megabirnavirus 1 (RnMBV1). L2-dsRNA/SsMBV1 comprises two nonoverlapping ORFs (ORFA and ORFB) encoding two hypothetical proteins with unknown functions. The 5'-terminal regions of L1- and L2-dsRNA/SsMBV1 share strictly conserved sequences and form stable stem-loop structures. Although L2-dsRNA/SsMBV1 is dispensable for replication, genome packaging, and pathogenicity of SsMBV1, it enhances transcript accumulation of L1-dsRNA/SsMBV1 and stability of virus-like particles (VLPs). Interestingly, a conserved papain-like protease domain similar to a multifunctional protein (p29) of Cryphonectria hypovirus 1 was detected in the ORFA-encoded protein of L2-dsRNA/SsMBV1. Phylogenetic analysis based on the protease domain suggests that horizontal gene transfer may have occurred from a single-stranded RNA (ssRNA) virus (hypovirus) to a dsRNA virus, SsMBV1. Our results reveal that SsMBV1 has a slight impact on the fundamental biological characteristics of its host regardless of the presence or absence of L2-dsRNA/SsMBV1.

IMPORTANCE

Mycoviruses are widespread in all major fungal groups, and they possess diverse genomes of mostly ssRNA and dsRNA and, recently, circular ssDNA. Here, we have characterized a novel dsRNA virus (Sclerotinia sclerotiorum megabirnavirus 1 [SsMBV1]) that was isolated from an apparently hypovirulent strain, SX466, of Sclerotinia sclerotiorum. Although SsMBV1 is phylogenetically related to RnMBV1, SsMBV1 is markedly distinct from other reported megabirnaviruses with two features of VLPs and conserved domains. Our results convincingly showed that SsMBV1 is viable in the absence of L2-dsRNA/SsMBV1 (a potential large satellite-like RNA or genuine genomic virus component). More interestingly, we detected a conserved papain-like protease domain that commonly exists in ssRNA viruses, including members of the families Potyviridae and Hypoviridae. Phylogenetic analysis based on the protease domain suggests that horizontal gene transfer might have occurred from an ssRNA virus to a dsRNA virus, which may provide new insights into the evolutionary history of dsRNA and ssRNA viruses.

The highly conserved codon following the slippery sequence supports -1 frameshift efficiency at the HIV-1 frameshift site

HIV-1 utilises -1 programmed ribosomal frameshifting to translate structural and enzymatic domains in a defined proportion required for replication. A slippery sequence, U UUU UUA, and a stem-loop are well-defined RNA features modulating -1 frameshifting in HIV-1. The GGG glycine codon immediately following the slippery sequence (the 'intercodon') contributes structurally to the start of the stem-loop but has no defined role in current models of the frameshift mechanism, as slippage is inferred to occur before the intercodon has reached the ribosomal decoding site. This GGG codon is highly conserved in natural isolates of HIV. When the natural intercodon was replaced with a stop codon two different decoding molecules-eRF1 protein or a cognate suppressor tRNA-were able to access and decode the intercodon prior to -1 frameshifting. This implies significant slippage occurs when the intercodon is in the (perhaps distorted) ribosomal A site. We accommodate the influence of the intercodon in a model of frame maintenance versus frameshifting in HIV-1.

HIV-1 replication and the cellular eukaryotic translation apparatus

Eukaryotic translation is a complex process composed of three main steps: initiation, elongation, and termination. During infections by RNA- and DNA-viruses, the eukaryotic translation machinery is used to assure optimal viral protein synthesis. Human immunodeficiency virus type I (HIV-1) uses several non-canonical pathways to translate its own proteins, such as leaky scanning, frameshifting, shunt, and cap-independent mechanisms. Moreover, HIV-1 modulates the host translation machinery by targeting key translation factors and overcomes different cellular obstacles that affect protein translation. In this review, we describe how HIV-1 proteins target several components of the eukaryotic translation machinery, which consequently improves viral translation and replication.

Local structural and environmental factors define the efficiency of an RNA pseudoknot involved in programmed ribosomal frameshift process

In programmed -1 ribosomal frameshift, an RNA pseudoknot stalls the ribosome at specific sequence and restarts translation in a new reading frame. A precise understanding of structural characteristics of these pseudoknots and their PRF inducing ability has not been clear to date. To investigate this phenomenon, we have studied various structural aspects of a -1 PRF inducing RNA pseudoknot from BWYV using extensive molecular dynamics simulations. A set of functional and poorly functional forms, for which previous mutational data were available, were chosen for analysis. These structures differ from each other by either single base substitutions or base-pair replacements from the native structure. We have rationalized how certain mutations in RNA pseudoknot affect its function; e.g., a specific base substitution in loop 2 stabilizes the junction geometry by forming multiple noncanonical hydrogen bonds, leading to a highly rigid structure that could effectively resist ribosome-induced unfolding, thereby increasing efficiency. While, a CG to AU pair substitution in stem 1 leads to loss of noncanonical hydrogen bonds between stems and loop, resulting in a less stable structure and reduced PRF inducing ability, inversion of a pair in stem 2 alters specific base-pair geometry that might be required in ribosomal recognition of nucleobase groups, negatively affecting pseudoknot functioning. These observations illustrate that the ability of an RNA pseudoknot to induce -1 PRF with an optimal rate depends on several independent factors that contribute to either the local conformational variability or geometry.

Structure and dynamics of the HIV-1 frameshift element RNA

The HIV-1 ribosomal frameshift element is highly structured, regulates translation of all virally encoded enzymes, and is a promising therapeutic target. The prior model for this motif contains two helices separated by a three-nucleotide bulge. Modifications to this model were suggested by SHAPE chemical probing of an entire HIV-1 RNA genome. Novel features of the SHAPE-directed model include alternate helical conformations and a larger, more complex structure. These structural elements also support the presence of a secondary frameshift site within the frameshift domain. Here, we use oligonucleotide-directed structure perturbation, probing in the presence of formamide, and in-virion experiments to examine these models. Our data support a model in which the frameshift domain is anchored by a stable helix outside the conventional domain. Less stable helices within the domain can switch from the SHAPE-predicted to the two-helix conformation. Translational frameshifting assays with frameshift domain mutants support a functional role for the interactions predicted by and specific to the SHAPE-directed model. These results reveal that the HIV-1 frameshift domain is a complex, dynamic structure and underscore the importance of analyzing folding in the context of full-length RNAs.