Experimental and computational methods for studying the dynamics of RNA-RNA interactions in SARS-COV2 genomes

Long-range ribonucleic acid (RNA)-RNA interactions (RRI) are prevalent in positive-strand RNA viruses, including Beta-coronaviruses, and these take part in regulatory roles, including the regulation of sub-genomic RNA production rates. Crosslinking of interacting RNAs and short read-based deep sequencing of resulting RNA-RNA hybrids have shown that these long-range structures exist in severe acute respiratory syndrome coronavirus (SARS-CoV)-2 on both genomic and sub-genomic levels and in dynamic topologies. Furthermore, co-evolution of coronaviruses with their hosts is navigated by genetic variations made possible by its large genome, high recombination frequency and a high mutation rate. SARS-CoV-2's mutations are known to occur spontaneously during replication, and thousands of aggregate mutations have been reported since the emergence of the virus. Although many long-range RRIs have been experimentally identified using high-throughput methods for the wild-type SARS-CoV-2 strain, evolutionary trajectory of these RRIs across variants, impact of mutations on RRIs and interaction of SARS-CoV-2 RNAs with the host have been largely open questions in the field. In this review, we summarize recent computational tools and experimental methods that have been enabling the mapping of RRIs in viral genomes, with a specific focus on SARS-CoV-2. We also present available informatics resources to navigate the RRI maps and shed light on the impact of mutations on the RRI space in viral genomes. Investigating the evolution of long-range RNA interactions and that of virus-host interactions can contribute to the understanding of new and emerging variants as well as aid in developing improved RNA therapeutics critical for combating future outbreaks.

Improved prime editing allows for routine predictable gene editing in Physcomitrium patens

Efficient and precise gene editing is the gold standard of any reverse genetic study. The recently developed prime editing approach, a modified CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein] editing method, has reached the precision goal but its editing rate can be improved. We present an improved methodology that allows for routine prime editing in the model plant Physcomitrium patens, whilst exploring potential new prime editing improvements. Using a standardized protoplast transfection procedure, multiple prime editing guide RNA (pegRNA) structural and prime editor variants were evaluated targeting the APT reporter gene through direct plant selection. Together, enhancements of expression of the prime editor, modifications of the 3' extension of the pegRNA, and the addition of synonymous mutation in the reverse transcriptase template sequence of the pegRNA dramatically improve the editing rate without affecting the quality of the edits. Furthermore, we show that prime editing is amenable to edit a gene of interest through indirect selection, as demonstrated by the generation of a Ppdek10 mutant. Additionally, we determine that a plant retrotransposon reverse transcriptase enables prime editing. Finally, we show for the first time the possibility of performing prime editing with two independently coded peptides.

Structure and function analysis of a type III preQ1-I riboswitch from Escherichia coli reveals direct metabolite sensing by the Shine-Dalgarno sequence

Riboswitches are small noncoding RNAs found primarily in the 5' leader regions of bacterial messenger RNAs where they regulate expression of downstream genes in response to binding one or more cellular metabolites. Such noncoding RNAs are often regulated at the translation level, which is thought to be mediated by the accessibility of the Shine-Dalgarno sequence (SDS) ribosome-binding site. Three classes (I-III) of prequeuosine1 (preQ1)-sensing riboswitches are known that control translation. Class I is divided into three subtypes (types I-III) that have diverse mechanisms of sensing preQ1, which is involved in queuosine biosynthesis. To provide insight into translation control, we determined a 2.30 Å-resolution cocrystal structure of a class I type III preQ1-sensing riboswitch identified in Escherichia coli (Eco) by bioinformatic searches. The Eco riboswitch structure differs from previous preQ1 riboswitch structures because it has the smallest naturally occurring aptamer and the SDS directly contacts the preQ1 metabolite. We validated structural observations using surface plasmon resonance and in vivo gene-expression assays, which showed strong switching in live E. coli. Our results demonstrate that the Eco riboswitch is relatively sensitive to mutations that disrupt noncanonical interactions that form the pseudoknot. In contrast to type II preQ1 riboswitches, a kinetic analysis showed that the type III Eco riboswitch strongly prefers preQ1 over the chemically similar metabolic precursor preQ0. Our results reveal the importance of noncanonical interactions in riboswitch-driven gene regulation and the versatility of the class I preQ1 riboswitch pseudoknot as a metabolite-sensing platform that supports SDS sequestration.

Structural and Functional Insights into Viral Programmed Ribosomal Frameshifting

Protein synthesis by the ribosome is the final stage of biological information transfer and represents an irreversible commitment to gene expression. Accurate translation of messenger RNA is therefore essential to all life, and spontaneous errors by the translational machinery are highly infrequent (∼1/100,000 codons). Programmed -1 ribosomal frameshifting (-1PRF) is a mechanism in which the elongating ribosome is induced at high frequency to slip backward by one nucleotide at a defined position and to continue translation in the new reading frame. This is exploited as a translational regulation strategy by hundreds of RNA viruses, which rely on -1PRF during genome translation to control the stoichiometry of viral proteins. While early investigations of -1PRF focused on virological and biochemical aspects, the application of X-ray crystallography and cryo-electron microscopy (cryo-EM), and the advent of deep sequencing and single-molecule approaches have revealed unexpected structural diversity and mechanistic complexity. Molecular players from several model systems have now been characterized in detail, both in isolation and, more recently, in the context of the elongating ribosome. Here we provide a summary of recent advances and discuss to what extent a general model for -1PRF remains a useful way of thinking.

Structure of LARP7 Protein p65-telomerase RNA Complex in Telomerase Revealed by Cryo-EM and NMR

La-related protein 7 (LARP7) are a family of RNA chaperones that protect the 3'-end of RNA and are components of specific ribonucleoprotein complexes (RNP). In Tetrahymena thermophila telomerase, LARP7 protein p65 together with telomerase reverse transcriptase (TERT) and telomerase RNA (TER) form the core RNP. p65 has four known domains-N-terminal domain (NTD), La motif (LaM), RNA recognition motif 1 (RRM1), and C-terminal xRRM2. To date, only the xRRM2 and LaM and their interactions with TER have been structurally characterized. Conformational dynamics leading to low resolution in cryo-EM density maps have limited our understanding of how full-length p65 specifically recognizes and remodels TER for telomerase assembly. Here, we combined focused classification of Tetrahymena telomerase cryo-EM maps with NMR spectroscopy to determine the structure of p65-TER. Three previously unknown helices are identified, one in the otherwise intrinsically disordered NTD that binds the La module, one that extends RRM1, and another preceding xRRM2, that stabilize p65-TER interactions. The extended La module (αN, LaM and RRM1) interacts with the four 3' terminal U nucleotides, while LaM and αN additionally interact with TER pseudoknot, and LaM with stem 1 and 5' end. Our results reveal the extensive p65-TER interactions that promote TER 3'-end protection, TER folding, and core RNP assembly and stabilization. The structure of full-length p65 with TER also sheds light on the biological roles of genuine La and LARP7 proteins as RNA chaperones and core RNP components.

Computational design and experimental verification of pseudoknotted ribozymes

The design of new RNA sequences that retain the function of a model RNA structure is a challenge in bioinformatics because of the structural complexity of these molecules. RNA can fold into its secondary and tertiary structures by forming stem-loops and pseudoknots. A pseudoknot is a set of base pairs between a region within a stem-loop and nucleotides outside of this stem-loop; this motif is very important for numerous functional structures. It is important for any computational design algorithm to take into account these interactions to give a reliable result for any structures that include pseudoknots. In our study, we experimentally validated synthetic ribozymes designed by Enzymer, which implements algorithms allowing for the design of pseudoknots. Enzymer is a program that uses an inverse folding approach to design pseudoknotted RNAs; we used it in this study to design two types of ribozymes. The ribozymes tested were the hammerhead and the glmS, which have a self-cleaving activity that allows them to liberate the new RNA genome copy during rolling-circle replication or to control the expression of the downstream genes, respectively. We demonstrated the efficiency of Enzymer by showing that the pseudoknotted hammerhead and glmS ribozymes sequences it designed were extensively modified compared to wild-type sequences and were still active.

Right Coronary Artery Button Pseudoaneurysm After the Modified Bentall Procedure

Anastomoses of the coronary buttons are the Achilles' heel of the modified Bentall procedure (MBP) for the repair of the aortic root and ascending aorta. We present a rare case of post-MBP right coronary artery button pseudoaneurysm in a 30-year-old man. The contained leak, attributed to a pseudoknot in the polypropylene suture, was visualized via computed tomography angiography and transesophageal echocardiogram and repaired under deep hypothermic circulatory arrest.

Pseudoknots in RNA Structure Prediction

RNA molecules play active roles in the cell and are important for numerous applications in biotechnology and medicine. The function of an RNA molecule stems from its structure. RNA structure determination is time consuming, challenging, and expensive using experimental methods. Thus, much research has been directed at RNA structure prediction through computational means. Many of these methods focus primarily on the secondary structure of the molecule, ignoring the possibility of pseudoknotted structures. However, pseudoknots are known to play functional roles in many RNA molecules or in their method of interaction with other molecules. Improving the accuracy and efficiency of computational methods that predict pseudoknots is an ongoing challenge for single RNA molecules, RNA-RNA interactions, and RNA-protein interactions. To improve the accuracy of prediction, many methods focus on specific applications while restricting the length and the class of the pseudoknotted structures they can identify. In recent years, computational methods for structure prediction have begun to catch up with the impressive developments seen in biotechnology. Here, we provide a non-comprehensive overview of available pseudoknot prediction methods and their best-use cases. © 2023 Wiley Periodicals LLC.

Genomic Analysis of Non-B Nucleic Acids Structures in SARS-CoV-2: Potential Key Roles for These Structures in Mutability, Translation, and Replication?

Non-B nucleic acids structures have arisen as key contributors to genetic variation in SARS-CoV-2. Herein, we investigated the presence of defining spike protein mutations falling within inverted repeats (IRs) for 18 SARS-CoV-2 variants, discussed the potential roles of G-quadruplexes (G4s) in SARS-CoV-2 biology, and identified potential pseudoknots within the SARS-CoV-2 genome. Surprisingly, there was a large variation in the number of defining spike protein mutations arising within IRs between variants and these were more likely to occur in the stem region of the predicted hairpin stem-loop secondary structure. Notably, mutations implicated in ACE2 binding and propagation (e.g., ΔH69/V70, N501Y, and D614G) were likely to occur within IRs, whilst mutations involved in antibody neutralization and reduced vaccine efficacy (e.g., T19R, ΔE156, ΔF157, R158G, and G446S) were rarely found within IRs. We also predicted that RNA pseudoknots could predominantly be found within, or next to, 29 mutations found in the SARS-CoV-2 spike protein. Finally, the Omicron variants BA.2, BA.4, BA.5, BA.2.12.1, and BA.2.75 appear to have lost two of the predicted G4-forming sequences found in other variants. These were found in nsp2 and the sequence complementary to the conserved stem-loop II-like motif (S2M) in the 3' untranslated region (UTR). Taken together, non-B nucleic acids structures likely play an integral role in SARS-CoV-2 evolution and genetic diversity.

Replication of alphaviruses requires a pseudoknot that involves the poly(A) tail

Alphaviruses, such as the Sindbis virus and the Chikungunya virus, are RNA viruses with a positive sense single-stranded RNA genome that infect various vertebrates, including humans. A conserved sequence element (CSE) of ∼19 nt in the 3' noncoding region is important for replication. Despite extensive mutational analysis of the CSE, no comprehensive model of this element exists to date. Here, it is shown that the CSE can form an RNA pseudoknot with part of the poly(A) tail and is similar to the human telomerase pseudoknot with which it shares 17 nt. Mutants that alter the stability of the pseudoknot were investigated in the context of a replicon of the Sindbis virus and by native gel electrophoresis. These studies reveal that the pseudoknot is required for virus replication and is stabilized by UAU base triples. The new model is discussed in relation to previous data on Sindbis virus mutants and revertants lacking (part of) the CSE.

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.

Comparative Analysis of Novel Strains of Porcine Astrovirus Type 3 in the USA

Porcine astrovirus type 3 (PoAstV3) has been previously identified as a cause of polioencephalomyelitis in swine and continues to cause disease in the US swine industry. Herein, we describe the characterization of both untranslated regions, frameshifting signal, putative genome-linked virus protein (VPg) and conserved antigenic epitopes of several novel PoAstV3 genomes. Twenty complete coding sequences (CDS) were obtained from 32 diagnostic cases originating from 11 individual farms/systems sharing a nucleotide (amino acid) percent identity of 89.74-100% (94.79-100%), 91.9-100% (96.3-100%) and 90.71-100% (93.51-100%) for ORF1a, ORF1ab and ORF2, respectively. Our results indicate that the 5'UTR of PoAstV3 is highly conserved highlighting the importance of this region in translation initiation while their 3'UTR is moderately conserved among strains, presenting alternative configurations including multiple putative protein binding sites and pseudoknots. Moreover, two predicted conserved antigenic epitopes were identified matching the 3' termini of VP27 of PoAstV3 USA strains. These epitopes may aid in the design and development of vaccine components and diagnostic assays useful to control outbreaks of PoAstV3-associated CNS disease. In conclusion, this is the first analysis predicting the structure of important regulatory motifs of neurotropic mamastroviruses, which differ from those previously described in human astroviruses.

Small molecule targeting of biologically relevant RNA tertiary and quaternary structures

Initial successes in developing small molecule ligands for non-coding RNAs have underscored their potential as therapeutic targets. More recently, these successes have been aided by advances in biophysical and structural techniques for identification and characterization of more complex RNA structures; these higher-level folds present protein-like binding pockets that offer opportunities to design small molecules that could achieve a degree of selectivity often hard to obtain at the primary and secondary structure level. More specifically, identification and small molecule targeting of RNA tertiary and quaternary structures have allowed researchers to probe several human diseases and have resulted in promising clinical candidates. In this review we highlight a selection of diverse and exciting successes and the experimental approaches that led to their discovery. These studies include examples of recent developments in RNA-centric assays and ligands that provide insight into the features responsible for the affinity and biological outcome of RNA-targeted chemical probes. This report highlights the potential and emerging opportunities to selectively target RNA tertiary and quaternary structures as a route to better understand and, ultimately, treat many diseases.

The RNA Architecture of the SARS-CoV-2 3'-Untranslated Region

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the current COVID-19 pandemic. The 3' untranslated region (UTR) of this β-CoV contains essential cis-acting RNA elements for the viral genome transcription and replication. These elements include an equilibrium between an extended bulged stem-loop (BSL) and a pseudoknot. The existence of such an equilibrium is supported by reverse genetic studies and phylogenetic covariation analysis and is further proposed as a molecular switch essential for the control of the viral RNA polymerase binding. Here, we report the SARS-CoV-2 3' UTR structures in cells that transcribe the viral UTRs harbored in a minigene plasmid and isolated infectious virions using a chemical probing technique, namely dimethyl sulfate (DMS)-mutational profiling with sequencing (MaPseq). Interestingly, the putative pseudoknotted conformation was not observed, indicating that its abundance in our systems is low in the absence of the viral nonstructural proteins (nsps). Similarly, our results also suggest that another functional cis-acting element, the three-helix junction, cannot stably form. The overall architectures of the viral 3' UTRs in the infectious virions and the minigene-transfected cells are almost identical.

ATTfold: RNA Secondary Structure Prediction With Pseudoknots Based on Attention Mechanism

Accurate RNA secondary structure information is the cornerstone of gene function research and RNA tertiary structure prediction. However, most traditional RNA secondary structure prediction algorithms are based on the dynamic programming (DP) algorithm, according to the minimum free energy theory, with both hard and soft constraints. The accuracy is particularly dependent on the accuracy of soft constraints (from experimental data like chemical and enzyme detection). With the elongation of the RNA sequence, the time complexity of DP-based algorithms will increase geometrically, as a result, they are not good at coping with relatively long sequences. Furthermore, due to the complexity of the pseudoknots structure, the secondary structure prediction method, based on traditional algorithms, has great defects which cannot predict the secondary structure with pseudoknots well. Therefore, few algorithms have been available for pseudoknots prediction in the past. The ATTfold algorithm proposed in this article is a deep learning algorithm based on an attention mechanism. It analyzes the global information of the RNA sequence via the characteristics of the attention mechanism, focuses on the correlation between paired bases, and solves the problem of long sequence prediction. Moreover, this algorithm also extracts the effective multi-dimensional features from a great number of RNA sequences and structure information, by combining the exclusive hard constraints of RNA secondary structure. Hence, it accurately determines the pairing position of each base, and obtains the real and effective RNA secondary structure, including pseudoknots. Finally, after training the ATTfold algorithm model through tens of thousands of RNA sequences and their real secondary structures, this algorithm was compared with four classic RNA secondary structure prediction algorithms. The results show that our algorithm significantly outperforms others and more accurately showed the secondary structure of RNA. As the data in RNA sequence databases increase, our deep learning-based algorithm will have superior performance. In the future, this kind of algorithm will be more indispensable.

The Halastavi árva Virus Intergenic Region IRES Promotes Translation by the Simplest Possible Initiation Mechanism

Dicistrovirus intergenic region internal ribosomal entry sites (IGR IRESs) do not require initiator tRNA, an AUG codon, or initiation factors and jumpstart translation from the middle of the elongation cycle via formation of IRES/80S complexes resembling the pre-translocation state. eEF2 then translocates the [codon-anticodon]-mimicking pseudoknot I (PKI) from ribosomal A sites to P sites, bringing the first sense codon into the decoding center. Halastavi árva virus (HalV) contains an IGR that is related to previously described IGR IRESs but lacks domain 2, which enables these IRESs to bind to individual 40S ribosomal subunits. By using in vitro reconstitution and cryoelectron microscopy (cryo-EM), we now report that the HalV IGR IRES functions by the simplest initiation mechanism that involves binding to 80S ribosomes such that PKI is placed in the P site, so that the A site contains the first codon that is directly accessible for decoding without prior eEF2-mediated translocation of PKI.

Mechanisms of Theta Plasmid Replication in Enterobacteria and Implications for Adaptation to Its Host

Plasmids are autonomously replicating sequences that help cells adapt to diverse stresses. Theta plasmids are the most frequent plasmid class in enterobacteria. They co-opt two host replication mechanisms: replication at oriC, a DnaA-dependent pathway leading to replisome assembly (theta class A), and replication fork restart, a PriA-dependent pathway leading to primosome assembly through primer extension and D-loop formation (theta classes B, C, and D). To ensure autonomy from the host's replication and to facilitate copy number regulation, theta plasmids have unique mechanisms of replication initiation at the plasmid origin of replication (ori). Tight plasmid copy number regulation is essential because of the major and direct impact plasmid gene dosage has on gene expression. The timing of plasmid replication and segregation are also critical for optimizing plasmid gene expression. Therefore, we propose that plasmid replication needs to be understood in its biological context, where complex origins of replication (redundant origins, mosaic and cointegrated replicons), plasmid segregation, and toxin-antitoxin systems are often present. Highlighting their tight functional integration with ori function, we show that both partition and toxin-antitoxin systems tend to be encoded in close physical proximity to the ori in a large collection of Escherichia coli plasmids. We also propose that adaptation of plasmids to their host optimizes their contribution to the host's fitness while restricting access to broad genetic diversity, and we argue that this trade-off between adaptation to host and access to genetic diversity is likely a determinant factor shaping the distribution of replicons in populations of enterobacteria.

The RNA of Maize Chlorotic Mottle Virus, an Obligatory Component of Maize Lethal Necrosis Disease, Is Translated via a Variant Panicum Mosaic Virus-Like Cap-Independent Translation Element

Maize chlorotic mottle virus (MCMV) combines with a potyvirus in maize lethal necrosis disease (MLND), a serious emerging disease worldwide. To inform resistance strategies, we characterized the translation initiation mechanism of MCMV. We report that MCMV RNA contains a cap-independent translation element (CITE) in its 3' untranslated region (UTR). The MCMV 3' CITE (MTE) was mapped to nucleotides 4164 to 4333 in the genomic RNA. 2'-Hydroxyl acylation analyzed by primer extension (SHAPE) probing revealed that the MTE is a distinct variant of the panicum mosaic virus-like 3' CITE (PTE). Like the PTE, electrophoretic mobility shift assays (EMSAs) indicated that eukaryotic translation initiation factor 4E (eIF4E) binds the MTE despite the absence of an m7GpppN cap structure, which is normally required for eIF4E to bind RNA. Using a luciferase reporter system, mutagenesis to disrupt and restore base pairing revealed that the MTE interacts with the 5' UTRs of both genomic RNA and subgenomic RNA1 via long-distance kissing stem-loop interaction to facilitate translation. The MTE stimulates a relatively low level of translation and has a weak, if any, pseudoknot, which is present in the most active PTEs, mainly because the MTE lacks the pyrimidine-rich tract that base pairs to a G-rich bulge to form the pseudoknot. However, most mutations designed to form a pseudoknot decreased translation activity. Mutations in the viral genome that reduced or restored translation prevented and restored virus replication, respectively, in maize protoplasts and in plants. In summary, the MTE differs from the canonical PTE but falls into a structurally related class of 3' CITEs.IMPORTANCE In the past decade, maize lethal necrosis disease has caused massive crop losses in East Africa. It has also emerged in China and parts of South America. Maize chlorotic mottle virus (MCMV) infection is required for this disease. While some tolerant maize lines have been identified, there are no known resistance genes that confer immunity to MCMV. In order to improve resistance strategies against MCMV, we focused on how the MCMV genome is translated, the first step of gene expression by all positive-strand RNA viruses. We identified a structure (cap-independent translation element) in the 3' untranslated region of the viral RNA genome that allows the virus to usurp a host translation initiation factor, eIF4E, in a way that differs from host mRNA interactions with the translational machinery. This difference indicates eIF4E may be a soft target for engineering of-or breeding for-resistance to MCMV.

Nucleobase mutants of a bacterial preQ1-II riboswitch that uncouple metabolite sensing from gene regulation

Riboswitches are a class of nonprotein-coding RNAs that directly sense cellular metabolites to regulate gene expression. They are model systems for analyzing RNA-ligand interactions and are established targets for antibacterial agents. Many studies have analyzed the ligand-binding properties of riboswitches, but this work has outpaced our understanding of the underlying chemical pathways that govern riboswitch-controlled gene expression. To address this knowledge gap, we prepared 15 mutants of the preQ1-II riboswitch-a structurally and biochemically well-characterized HLout pseudoknot that recognizes the metabolite prequeuosine1 (preQ1). The mutants span the preQ1-binding pocket through the adjoining Shine-Dalgarno sequence (SDS) and include A-minor motifs, pseudoknot-insertion helix P4, U·A-U base triples, and canonical G-C pairs in the anti-SDS. As predicted-and confirmed by in vitro isothermal titration calorimetry measurements-specific mutations ablated preQ1 binding, but most aberrant binding effects were corrected by compensatory mutations. In contrast, functional analysis in live bacteria using a riboswitch-controlled GFPuv-reporter assay revealed that each mutant had a deleterious effect on gene regulation, even when compensatory changes were included. Our results indicate that effector binding can be uncoupled from gene regulation. We attribute loss of function to defects in a chemical interaction network that links effector binding to distal regions of the fold that support the gene-off RNA conformation. Our findings differentiate effector binding from biological function, which has ramifications for riboswitch characterization. Our results are considered in the context of synthetic ligands and drugs that bind tightly to riboswitches without eliciting a biological response.

Comparison of Pseudoknotted RNA Secondary Structures by Topological Centroid Identification and Tree Edit Distance

Comparison of RNA structures is one of the most crucial analysis for elucidating their individual functions and promoting medical applications. Because it is widely accepted that their functions and structures are strongly correlated, various methods for RNA secondary structure analysis have been proposed owing to the difficulty in predicting RNA three-dimensional structure directly from its sequence. However, there are few methods dealing with RNA secondary structures with a specific and complex partial structure called pseudoknot despite its significance to biological process, which is a big obstacle for analyzing their functions. In this study, we propose a novel tree representation of pseudoknotted RNA secondary structures by topological centroid identification and their comparison methods based on the tree edit distance. In the proposed method, a given graph representing an RNA secondary structure is transformed to a tree rooted at one of the vertices constituting the topological centroid that is identified by removing cycles with peeling processing for the graph. When comparing tree-represented RNA secondary structures collected from a public database using the tree edit distance and functional gene groups defined by Gene Ontology (GO), the proposed method showed better clustering results according to their GOs than canonical RNA sequence-based comparison. In addition, we also report a case that the combination of the tree edit distance and the sequence edit distance shows a better classification of the pseudoknotted RNA secondary structures.

The conserved structure of plant telomerase RNA provides the missing link for an evolutionary pathway from ciliates to humans

Telomerase is essential for maintaining telomere integrity. Although telomerase function is widely conserved, the integral telomerase RNA (TR) that provides a template for telomeric DNA synthesis has diverged dramatically. Nevertheless, TR molecules retain 2 highly conserved structural domains critical for catalysis: a template-proximal pseudoknot (PK) structure and a downstream stem-loop structure. Here we introduce the authentic TR from the plant Arabidopsis thaliana, called AtTR, identified through next-generation sequencing of RNAs copurifying with Arabidopsis TERT. This RNA is distinct from the RNA previously described as the templating telomerase RNA, AtTER1. AtTR is a 268-nt Pol III transcript necessary for telomere maintenance in vivo and sufficient with TERT to reconstitute telomerase activity in vitro. Bioinformatics analysis identified 85 AtTR orthologs from 3 major clades of plants: angiosperms, gymnosperms, and lycophytes. Through phylogenetic comparisons, a secondary structure model conserved among plant TRs was inferred and verified using in vitro and in vivo chemical probing. The conserved plant TR structure contains a template-PK core domain enclosed by a P1 stem and a 3' long-stem P4/5/6, both of which resemble a corresponding structural element in ciliate and vertebrate TRs. However, the plant TR contains additional stems and linkers within the template-PK core, allowing for expansion of PK structure from the simple PK in the smaller ciliate TR during evolution. Thus, the plant TR provides an evolutionary bridge that unites the disparate structures of previously characterized TRs from ciliates and vertebrates.

Structure and function of cis-acting RNA elements of flavivirus

The genus Flavivirus is a group of single-stranded, positive-sense RNA viruses that includes numerous human pathogens with global impact, such as dengue virus (DENV), yellow fever virus (YFV), West Nile virus (WNV), and Zika virus (ZIKV). The approximately 11-kilobase genome is flanked by highly structured untranslated regions (UTRs), which contain various cis-acting RNA elements with unique structures and functions. Moreover, local RNA elements circularize the genome non-covalently through long-range interactions. Interestingly, many flavivirus cis-acting RNA elements contain group-specific motifs or are specific for the given phylogenetic groups, suggesting their potential association with flavivirus evolution and diversification. In this review, we summarize recent advances about the structure and function of cis-acting RNA elements in flavivirus genomes and highlight the potential implications for flavivirus evolution. Finally, the scientific questions remained to be answered in the field are also discussed.

A bacterial riboswitch class for the thiamin precursor HMP-PP employs a terminator-embedded aptamer

We recently implemented a bioinformatics pipeline that can uncover novel, but rare, riboswitch candidates as well as other noncoding RNA structures in bacteria. A prominent candidate revealed by our initial search efforts was called the ‘thiS motif’ because of its frequent association with a gene coding for the ThiS protein, which delivers sulfur to form the thiazole moiety of the thiamin precursor HET-P. In the current report, we describe biochemical and genetic data demonstrating that thiS motif RNAs function as sensors of the thiamin precursor HMP-PP, which is fused with HET-P ultimately to form the final active coenzyme thiamin pyrophosphate (TPP). HMP-PP riboswitches exhibit a distinctive architecture wherein an unusually small ligand-sensing aptamer is almost entirely embedded within an otherwise classic intrinsic transcription terminator stem. This arrangement yields remarkably compact genetic switches that bacteria use to tune the levels of thiamin precursors during the biosynthesis of this universally distributed coenzyme.

Many bacteria use small genetic devices called riboswitches to sense molecules that are essential for life and regulate the genes necessary to make, break or move these molecules. Riboswitches are made of molecules of RNA and appear to have ancient origins that predate the evolution of bacteria and other lifeforms made of cells. Inside modern bacteria, chunks of DNA in the genome provide the instructions to make riboswitches and around 40 different types of riboswitch have been identified so far. However, it has been proposed that the instructions for thousands more riboswitches may remain hidden in the DNA of bacteria.

All of the currently known riboswitches contain a region called an aptamer that binds to a target molecule. This binding causes another structure in the riboswitch RNA to switch a specific gene on or off. For example, the aptamer binding might cause a hairpin-like structure called a terminator to form, which stops a gene being used to make new RNA molecules.

In 2019 a team of researchers reported using a computational approach to identify new riboswitches in bacteria. This approach identified many different chunks of DNA that might code for a riboswitch, including one known as the thiS motif. This potential new riboswitch appeared to be associated with a gene that encodes a protein required to make a vitamin called thiamin (also known as vitamin B₁).

To test the new computational approach, Atilho et al. including several of the researchers involved in the earlier work used genetic and biochemical techniques to study the thiS motif. The experiments revealed that the motif binds to a molecule called HMP-PP, which bacteria use to make thiamin. Unexpectedly, the aptamer of the riboswitch was nested within a terminator, rather than being a separate entity.

The findings of Atilho et al. reveal that riboswitches can be even more compact than previously thought. Furthermore, these findings reveal new insights into how bacteria use riboswitches to manage their vitamin levels. In the future it may be possible to develop drugs that target such riboswitches to starve bacteria of these essential molecules.

DMfold: A Novel Method to Predict RNA Secondary Structure With Pseudoknots Based on Deep Learning and Improved Base Pair Maximization Principle

While predicting the secondary structure of RNA is vital for researching its function, determining RNA secondary structure is challenging, especially for that with pseudoknots. Typically, several excellent computational methods can be utilized to predict the secondary structure (with or without pseudoknots), but they have their own merits and demerits. These methods can be classified into two categories: the multi-sequence method and the single-sequence method. The main advantage of the multi-sequence method lies in its use of the auxiliary sequences to assist in predicting the secondary structure, but it can only successfully predict in the presence of multiple highly homologous sequences. The single-sequence method is associated with the major merit of easy operation (only need the target sequence to predict secondary structure), but its folding parameters are the common features of diversity RNA, which cannot describe the unique characteristics of RNA, thus potentially resulting in the low prediction accuracy in some RNA. In this paper, "DMfold," a method based on the Deep Learning and Improved Base Pair Maximization Principle, is proposed to predict the secondary structure with pseudoknots, which fully absorbs the advantages and avoids some disadvantages of those two methods. Notably, DMfold could predict the secondary structure of RNA by learning similar RNA in the known structures, which uses the similar RNA sequences instead of the highly homogeneous sequences in the multi-sequence method, thereby reducing the requirement for auxiliary sequences. In DMfold, it only needs to input the target sequence to predict the secondary structure. Its folding parameters are fully extracted automatically by deep learning, which could avoid the lack of folding parameters in the single-sequence method. Experiments show that our method is not only simple to operate, but also improves the prediction accuracy compared to multiple excellent prediction methods. A repository containing our code can be found at https://github.com/linyuwangPHD/RNA-Secondary-Structure-Database.

The 3' Untranslated Region of a Plant Viral RNA Directs Efficient Cap-Independent Translation in Plant and Mammalian Systems

Many plant viral RNA genomes lack a 5′ cap, and instead are translated via a cap-independent translation element (CITE) in the 3′ untranslated region (UTR). The panicum mosaic virus-like CITE (PTE), found in many plant viral RNAs, binds and requires the cap-binding translation initiation factor eIF4E to facilitate translation. eIF4E is structurally conserved between plants and animals, so we tested cap-independent translation efficiency of PTEs of nine plant viruses in plant and mammalian systems. The PTE from thin paspalum asymptomatic virus (TPAV) facilitated efficient cap-independent translation in wheat germ extract, rabbit reticulocyte lysate, HeLa cell lysate, and in oat and mammalian (BHK) cells. Human eIF4E bound the TPAV PTE but not a PTE that did not stimulate cap-independent translation in mammalian extracts or cells. Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) footprinting revealed that both human and wheat eIF4E protected the conserved guanosine (G)-rich domain in the TPAV PTE pseudoknot. The central G plays a key role, as it was found to be required for translation and protection from SHAPE modification by eIF4E. These results provide insight on how plant viruses gain access to the host’s translational machinery, an essential step in infection, and raise the possibility that similar PTE-like mechanisms may exist in mRNAs of mammals or their viruses.

A Boltzmann Sampler for 1-Pairs with Double Filtration

Recently, a framework considering RNA sequences and their RNA secondary structures as pairs led to some information-theoretic perspectives on how the semantics encoded in RNA sequences can be inferred. This pairing arises naturally from the energy model of RNA secondary structures. Fixing the sequence in the pairing produces the RNA energy landscape, whose partition function was discovered by McCaskill. Dually, fixing the structure induces the energy landscape of sequences. The latter has been considered originally for designing more efficient inverse folding algorithms and subsequently enhanced by facilitating the sampling of sequences. We present here a partition function of sequence/structure pairs, with endowed Hamming distance and base pair distance filtration. This partition function is an augmentation of the previous mentioned (dual) partition function. We develop an efficient dynamic programming routine to recursively compute the partition function with this double filtration. Our framework is capable of dealing with RNA secondary structures as well as 1-structures, where a 1-structure is an RNA pseudoknot structure consisting of "building blocks" of genus 0 or 1. In particular, 0-structures, consisting of only "building blocks" of genus 0, are exactly RNA secondary structures. The time complexity for calculating the partition function of 1-pairs, that is, sequence/structure pairs where the structures are 1-structures, is O(h³b³n⁶), where h, b, n denote the Hamming distance, base pair distance, and sequence length, respectively. The time complexity for the partition function of 0-pairs is O(h²b²n³).

Effects of flanking regions on HDV cotranscriptional folding kinetics

Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.

Identification of the RNA Pseudoknot within the 3' End of the Porcine Reproductive and Respiratory Syndrome Virus Genome as a Pathogen-Associated Molecular Pattern To Activate Antiviral Signaling via RIG-I and Toll-Like Receptor 3

Once infected by viruses, cells can detect pathogen-associated molecular patterns (PAMPs) on viral nucleic acid by host pattern recognition receptors (PRRs) to initiate the antiviral response. Porcine reproductive and respiratory syndrome virus (PRRSV) is the causative agent of porcine reproductive and respiratory syndrome (PRRS), characterized by reproductive failure in sows and respiratory diseases in pigs of different ages. To date, the sensing mechanism of PRRSV has not been elucidated. Here, we reported that the pseudoknot region residing in the 3' untranslated regions (UTR) of the PRRSV genome, which has been proposed to regulate RNA synthesis and virus replication, was sensed as nonself by retinoic acid-inducible gene I (RIG-I) and Toll-like receptor 3 (TLR3) and strongly induced type I interferons (IFNs) and interferon-stimulated genes (ISGs) in porcine alveolar macrophages (PAMs). The interaction between the two stem-loops inside the pseudoknot structure was sufficient for IFN induction, since disruption of the pseudoknot interaction powerfully dampened the IFN induction. Furthermore, transfection of the 3' UTR pseudoknot transcripts in PAMs inhibited PRRSV replication in vitro Importantly, the predicted similar structures of other arterivirus members, including equine arteritis virus (EAV), lactate dehydrogenase-elevating virus (LDV), and simian hemorrhagic fever virus (SHFV), also displayed strong IFN induction activities. Together, in this work we identified an innate recognition mechanism by which the PRRSV 3' UTR pseudoknot region served as PAMPs of arteriviruses and activated innate immune signaling to produce IFNs that inhibit virus replication. All of these results provide novel insights into innate immune recognition during virus infection.IMPORTANCE PRRS is the most common viral disease in the pork industry. It is caused by PRRSV, a positive single-stranded RNA virus, whose infection often leads to persistent infection. To date, it is not yet clear how PRRSV is recognized by the host and what is the exact mechanism of IFN induction. Here, we investigated the nature of PAMPs on PRRSV and the associated PRRs. We found that the 3' UTR pseudoknot region of PRRSV, which has been proposed to regulate viral RNA synthesis, could act as PAMPs recognized by RIG-I and TLR3 to induce type I IFN production to suppress PRRSV infection. This report is the first detailed description of pattern recognition for PRRSV, which is important in understanding the antiviral response of arteriviruses, especially PRRSV, and extends our knowledge on virus recognition.

New tRNA contacts facilitate ligand binding in a Mycobacterium smegmatis T box riboswitch

T box riboswitches are RNA regulatory elements widely used by organisms in the phyla Firmicutes and Actinobacteria to regulate expression of amino acid-related genes. Expression of T box family genes is down-regulated by transcription attenuation or inhibition of translation initiation in response to increased charging of the cognate tRNA. Three direct contacts with tRNA have been described; however, one of these contacts is absent in a subclass of T box RNAs and the roles of several structural domains conserved in most T box RNAs are unknown. In this study, structural elements of a Mycobacterium smegmatis ileS T box riboswitch variant with an Ultrashort (US) Stem I were sequentially deleted, which resulted in a progressive decrease in binding affinity for the tRNA^Ile ligand. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) revealed structural changes in conserved riboswitch domains upon interaction with the tRNA ligand. Cross-linking and mutational analyses identified two interaction sites, one between the S-turn element in Stem II and the T arm of tRNA^Ile and the other between the Stem IIA/B pseudoknot and the D loop of tRNA^Ile These newly identified RNA contacts add information about tRNA recognition by the T box riboswitch and demonstrate a role for the S-turn and pseudoknot elements, which resemble structural elements that are common in many cellular RNAs.

Evaluating the quality of SHAPE data simulated by k-mers for RNA structure prediction

Finding an effective measure to predict a more accurate RNA secondary structure is a challenging problem. In the last decade, an experimental method, known as selective [Formula: see text]-hydroxyl acylation analyzed by primer extension (SHAPE), was proposed to measure the tendency of forming a base pair for almost all nucleotides in an RNA sequence. These SHAPE reactivities are then utilized to improve the accuracy of RNA structure prediction. Due to a significant impact of SHAPE reactivity and in order to reduce the experimental costs, we propose a new model called HL-k-mer. This model simulates the SHAPE reactivity for each nucleotide in an RNA sequence. This is done by fetching the SHAPE reactivities for all sub-sequences of length k (k-mers) appearing in helix and loop regions. For evaluating the quality of simulated SHAPE data, ESD-Fold method is used based on the SHAPE data simulated by the HL-k-mer model ([Formula: see text]). Also, for further evaluation of simulated SHAPE data, three different methods are employed. We also extend this model to simulate the SHAPE data for the RNA pseudoknotted structure. The results indicate that the average accuracies of prediction using the SHAPE data simulated by our models (for [Formula: see text]) are higher compared to the experimental SHAPE data.

Programmed Ribosomal Frameshifting Generates a Copper Transporter and a Copper Chaperone from the Same Gene

Metal efflux pumps maintain ion homeostasis in the cell. The functions of the transporters are often supported by chaperone proteins, which scavenge the metal ions from the cytoplasm. Although the copper ion transporter CopA has been known in Escherichia coli, no gene for its chaperone had been identified. We show that the CopA chaperone is expressed in E. coli from the same gene that encodes the transporter. Some ribosomes translating copA undergo programmed frameshifting, terminate translation in the -1 frame, and generate the 70 aa-long polypeptide CopA(Z), which helps cells survive toxic copper concentrations. The high efficiency of frameshifting is achieved by the combined stimulatory action of a "slippery" sequence, an mRNA pseudoknot, and the CopA nascent chain. Similar mRNA elements are not only found in the copA genes of other bacteria but are also present in ATP7B, the human homolog of copA, and direct ribosomal frameshifting in vivo.

Structure-Function Model for Kissing Loop Interactions That Initiate Dimerization of Ty1 RNA

The genomic RNA of the retrotransposon Ty1 is packaged as a dimer into virus-like particles. The 5′ terminus of Ty1 RNA harbors cis-acting sequences required for translation initiation, packaging and initiation of reverse transcription (TIPIRT). To identify RNA motifs involved in dimerization and packaging, a structural model of the TIPIRT domain in vitro was developed from single-nucleotide resolution RNA structural data. In general agreement with previous models, the first 326 nucleotides of Ty1 RNA form a pseudoknot with a 7-bp stem (S1), a 1-nucleotide interhelical loop and an 8-bp stem (S2) that delineate two long, structured loops. Nucleotide substitutions that disrupt either pseudoknot stem greatly reduced helper-Ty1-mediated retrotransposition of a mini-Ty1, but only mutations in S2 destabilized mini-Ty1 RNA in cis and helper-Ty1 RNA in trans. Nested in different loops of the pseudoknot are two hairpins with complementary 7-nucleotide motifs at their apices. Nucleotide substitutions in either motif also reduced retrotransposition and destabilized mini- and helper-Ty1 RNA. Compensatory mutations that restore base-pairing in the S2 stem or between the hairpins rescued retrotransposition and RNA stability in cis and trans. These data inform a model whereby a Ty1 RNA kissing complex with two intermolecular kissing-loop interactions initiates dimerization and packaging.

Single-molecule FRET-Rosetta reveals RNA structural rearrangements during human telomerase catalysis

Maintenance of telomeres by telomerase permits continuous proliferation of rapidly dividing cells, including the majority of human cancers. Despite its direct biomedical significance, the architecture of the human telomerase complex remains unknown. Generating homogeneous telomerase samples has presented a significant barrier to developing improved structural models. Here we pair single-molecule Förster resonance energy transfer (smFRET) measurements with Rosetta modeling to map the conformations of the essential telomerase RNA core domain within the active ribonucleoprotein. FRET-guided modeling places the essential pseudoknot fold distal to the active site on a protein surface comprising the C-terminal element, a domain that shares structural homology with canonical polymerase thumb domains. An independently solved medium-resolution structure of Tetrahymena telomerase provides a blind test of our modeling methodology and sheds light on the structural homology of this domain across diverse organisms. Our smFRET-Rosetta models reveal nanometer-scale rearrangements within the RNA core domain during catalysis. Taken together, our FRET data and pseudoatomic molecular models permit us to propose a possible mechanism for how RNA core domain rearrangement is coupled to template hybrid elongation.

An Adaptive Defect Weighted Sampling Algorithm to Design Pseudoknotted RNA Secondary Structures

Computational design of RNA sequences that fold into targeted secondary structures has many applications in biomedicine, nanotechnology and synthetic biology. An RNA molecule is made of different types of secondary structure elements and an important RNA element named pseudoknot plays a key role in stabilizing the functional form of the molecule. However, due to the computational complexities associated with characterizing pseudoknotted RNA structures, most of the existing RNA sequence designer algorithms generally ignore this important structural element and therefore limit their applications. In this paper we present a new algorithm to design RNA sequences for pseudoknotted secondary structures. We use NUPACK as the folding algorithm to compute the equilibrium characteristics of the pseudoknotted RNAs, and describe a new adaptive defect weighted sampling algorithm named Enzymer to design low ensemble defect RNA sequences for targeted secondary structures including pseudoknots. We used a biological data set of 201 pseudoknotted structures from the Pseudobase library to benchmark the performance of our algorithm. We compared the quality characteristics of the RNA sequences we designed by Enzymer with the results obtained from the state of the art MODENA and antaRNA. Our results show our method succeeds more frequently than MODENA and antaRNA do, and generates sequences that have lower ensemble defect, lower probability defect and higher thermostability. Finally by using Enzymer and by constraining the design to a naturally occurring and highly conserved Hammerhead motif, we designed 8 sequences for a pseudoknotted cis-acting Hammerhead ribozyme. Enzymer is available for download at https://bitbucket.org/casraz/enzymer.

Structural and Functional Characterization of Programmed Ribosomal Frameshift Signals in West Nile Virus Strains Reveals High Structural Plasticity Among cis-Acting RNA Elements

West Nile virus (WNV) is a prototypical emerging virus for which no effective therapeutics currently exist. WNV uses programmed -1 ribosomal frameshifting (-1 PRF) to synthesize the NS1' protein, a C terminally extended version of its non-structural protein 1, the expression of which enhances neuro-invasiveness and viral RNA abundance. Here, the NS1' frameshift signals derived from four WNV strains were investigated to better understand -1 PRF in this quasispecies. Sequences previously predicted to promote -1 PRF strongly promote this activity, but frameshifting was significantly more efficient upon inclusion of additional 3' sequence information. The observation of different rates of -1 PRF, and by inference differences in the expression of NS1', may account for the greater degrees of pathogenesis associated with specific WNV strains. Chemical modification and mutational analyses of the longer and shorter forms of the -1 PRF signals suggests dynamic structural rearrangements between tandem stem-loop and mRNA pseudoknot structures in two of the strains. A model is suggested in which this is employed as a molecular switch to fine tune the relative expression of structural to non-structural proteins during different phases of the viral replication cycle.

Secondary structure confirmation and localization of Mg2+ ions in the mammalian CPEB3 ribozyme

Most of today's knowledge of the CPEB3 ribozyme, one of the few small self-cleaving ribozymes known to occur in humans, is based on comparative studies with the hepatitis delta virus (HDV) ribozyme, which is highly similar in cleavage mechanism and probably also in structure. Here we present detailed NMR studies of the CPEB3 ribozyme in order to verify the formation of the predicted nested double pseudoknot in solution. In particular, the influence of Mg(2+), the ribozyme's crucial cofactor, on the CPEB3 structure is investigated. NMR titrations, Tb(3+)-induced cleavage, as well as stoichiometry determination by hydroxyquinoline sulfonic acid fluorescence and equilibrium dialysis, are used to evaluate the number, location, and binding mode of Mg(2+)ions. Up to eight Mg(2+)ions interact site-specifically with the ribozyme, four of which are bound with high affinity. The global fold of the CPEB3 ribozyme, encompassing 80%-90% of the predicted base pairs, is formed in the presence of monovalent ions alone. Low millimolar concentrations of Mg(2+)promote a more compact fold and lead to the formation of additional structures in the core of the ribozyme, which contains the inner small pseudoknot and the active site. Several Mg(2+)binding sites, which are important for the functional fold, appear to be located in corresponding locations in the HDV and CPEB3 ribozyme, demonstrating the particular relevance of Mg(2+)for the nested double pseudoknot structure.

Reprogramming the genetic code: The emerging role of ribosomal frameshifting in regulating cellular gene expression

Reading frame maintenance is a critical property of ribosomes. However, a number of genetic elements have been described that can induce ribosomes to shift on mRNAs, the most well understood of which are a class that directs ribosomal slippage by one base in 5' (‐1) direction. This is referred to as programmed ‐1 ribosomal frameshifting (‐1 PRF). Recently, a new ‐1 PRF promoting element was serendipitously discovered in a study examining the effects of stretches of adenosines in the coding sequences of mRNAs. Here, we discuss this finding, recent studies describing how ‐1 PRF is used to control gene expression in eukaryotes, and how ‐1 PRF is itself regulated. The implications of dysregulation of ‐1 PRF on human health are examined, as are possible new areas in which novel ‐1 PRF promoting elements might be discovered.

Also watch the https://youtu.be/1mPXIINCRcY.

A dynamic RNA loop in an IRES affects multiple steps of elongation factor-mediated translation initiation

Internal ribosome entry sites (IRESs) are powerful model systems to understand how the translation machinery can be manipulated by structured RNAs and for exploring inherent features of ribosome function. The intergenic region (IGR) IRESs from the Dicistroviridae family of viruses are structured RNAs that bind directly to the ribosome and initiate translation by co-opting the translation elongation cycle. These IRESs require an RNA pseudoknot that mimics a codon-anticodon interaction and contains a conformationally dynamic loop. We explored the role of this loop and found that both the length and sequence are essential for translation in different types of IGR IRESs and from diverse viruses. We found that loop 3 affects two discrete elongation factor-dependent steps in the IRES initiation mechanism. Our results show how the IRES directs multiple steps after 80S ribosome placement and highlights the often underappreciated significance of discrete conformationally dynamic elements within the context of structured RNAs.

DOI:http://dx.doi.org/10.7554/eLife.08146.001

Many viruses store their genetic information in the form of strands of ribonucleic acid (RNA), which contain building blocks called nucleotides. Once inside an infected cell, the virus hijacks the cellular structures that build proteins (called ribosomes), which forces the cell to start making viral proteins.

Many RNA viruses manipulate the cell’s ribosomes using RNA elements called Internal Ribosome Entry Sites, or IRESs. In a family of viruses called Dicistroviridae, which infect a number of insects, a section of the IRES RNA binds directly to the ribosome. Proteins called elongation factors then trigger a series of events that lead to the cell starting to make the viral proteins.

By mutating the RNA of many different Dicistroviridae viruses that infect a variety of invertebrates, Ruehle et al. have now investigated how a particular loop in the structure of the IRES helps to make cells build the viral proteins. This loop is flexible, and interacts with the ribosome to enable the IRES to move through the ribosome. Mutations that shorten the loop or alter the sequence of nucleotides in the loop prevent the occurrence of two of the steps that need to occur for the cell to make viral proteins. Both of these steps depend on elongation factors. Determining how the entire IRES might change shape as it moves through the ribosome is an important next step, since the ribosome is exquisitely sensitive to the shape and motions of its binding partners.

DOI:http://dx.doi.org/10.7554/eLife.08146.002

Cell cycle control (and more) by programmed -1 ribosomal frameshifting: implications for disease and therapeutics

Like most basic molecular mechanisms, programmed –1 ribosomal frameshifting (−1 PRF) was first identified in viruses. Early observations that global dysregulation of −1 PRF had deleterious effects on yeast cell growth suggested that −1 PRF may be used to control cellular gene expression, and the cell cycle in particular. Collection of sufficient numbers of viral −1 PRF signals coupled with advances in computer sciences enabled 2 complementary computational approaches to identify −1 PRF signals in free living organisms. The unexpected observation that almost all −1 PRF events on eukaryotic mRNAs direct ribosomes to premature termination codons engendered the hypothesis that −1 PRF signals post-transcriptionally regulate gene expression by functioning as mRNA destabilizing elements. Emerging research suggests that some human diseases are associated with global defects in −1 PRF. The recent discovery of −1 PRF signal-specific trans-acting regulators may provide insight into novel therapeutic strategies aimed at treating diseases caused by changes in gene expression patterns.

A Highly Coupled Network of Tertiary Interactions in the SAM-I Riboswitch and Their Role in Regulatory Tuning

RNA folding in vivo is significantly influenced by transcription, which is not necessarily recapitulated by Mg(2+)-induced folding of the corresponding full-length RNA in vitro. Riboswitches that regulate gene expression at the transcriptional level are an ideal system for investigating this aspect of RNA folding as ligand-dependent termination is obligatorily co-transcriptional, providing a clear readout of the folding outcome. The folding of representative members of the SAM-I family of riboswitches has been extensively analyzed using approaches focusing almost exclusively upon Mg(2+) and/or S-adenosylmethionine (SAM)-induced folding of full-length transcripts of the ligand binding domain. To relate these findings to co-transcriptional regulatory activity, we have investigated a set of structure-guided mutations of conserved tertiary architectural elements of the ligand binding domain using an in vitro single-turnover transcriptional termination assay, complemented with phylogenetic analysis and isothermal titration calorimetry data. This analysis revealed a conserved internal loop adjacent to the SAM binding site that significantly affects ligand binding and regulatory activity. Conversely, most single point mutations throughout key conserved features in peripheral tertiary architecture supporting the SAM binding pocket have relatively little impact on riboswitch activity. Instead, a secondary structural element in the peripheral subdomain appears to be the key determinant in observed differences in regulatory properties across the SAM-I family. These data reveal a highly coupled network of tertiary interactions that promote high-fidelity co-transcriptional folding of the riboswitch but are only indirectly linked to regulatory tuning.

Influence of RNA structural elements on Ty1 retrotransposition

The long-terminal repeat (LTR)-retrotransposon Ty1 is a mobile genetic element that replicates through an RNA intermediate. Retroelement genomic transcripts contain internal structures fundamental to gene expression and propagation. In addition, long non-coding antisense RNAs overlap the 5′-terminal region of the genomic RNA and confer post-translational copy number control. Although LTR- retrotransposons are functionally related to retroviruses, little is known about the structural determinants required for genomic RNA packaging or reverse transcription. This commentary summarizes two recent papers that provide the first snapshot of genomic RNA structures from the retrotransposon Ty1 involved in transposition. We combined structural approaches with functional and genetic assays to determine if antisense RNAs anneal with the genomic RNA. Analysis of various steps in the Ty1 life cycle showed that a novel RNA pseudoknot contributes to retrotransposon function. Comparing different RNA states provides additional information about regions potentially involved in Ty1 RNA dimerization or packaging.

The structural basis of transfer RNA mimicry and conformational plasticity by a viral RNA

RNA is arguably the most functionally diverse biological macromolecule. In some cases a single discrete RNA sequence performs multiple roles, and this can be conferred by a complex three-dimensional structure. Such multifunctionality can also be driven or enhanced by the ability of a given RNA to assume different conformational (and therefore functional) states. Despite its biological importance, a detailed structural understanding of the paradigm of RNA structure-driven multifunctionality is lacking. To address this gap it is useful to study examples from single-stranded positive-sense RNA viruses, a prototype being the tRNA-like structure (TLS) found at the 3' end of the turnip yellow mosaic virus (TYMV). This TLS not only acts like a tRNA to drive aminoacylation of the viral genomic (g)RNA, but also interacts with other structures in the 3' untranslated region of the gRNA, contains the promoter for negative-strand synthesis, and influences several infection-critical processes. TLS RNA can provide a glimpse into the structural basis of RNA multifunctionality and plasticity, but for decades its high-resolution structure has remained elusive. Here we present the crystal structure of the complete TYMV TLS to 2.0 Å resolution. Globally, the RNA adopts a shape that mimics tRNA, but it uses a very different set of intramolecular interactions to achieve this shape. These interactions also allow the TLS to readily switch conformations. In addition, the TLS structure is 'two faced': one face closely mimics tRNA and drives aminoacylation, the other face diverges from tRNA and enables additional functionality. The TLS is thus structured to perform several functions and interact with diverse binding partners, and we demonstrate its ability to specifically bind to ribosomes.

RNA secondary structure modeling at consistent high accuracy using differential SHAPE

RNA secondary structure modeling is a challenging problem, and recent successes have raised the standards for accuracy, consistency, and tractability. Large increases in accuracy have been achieved by including data on reactivity toward chemical probes: Incorporation of 1M7 SHAPE reactivity data into an mfold-class algorithm results in median accuracies for base pair prediction that exceed 90%. However, a few RNA structures are modeled with significantly lower accuracy. Here, we show that incorporating differential reactivities from the NMIA and 1M6 reagents--which detect noncanonical and tertiary interactions--into prediction algorithms results in highly accurate secondary structure models for RNAs that were previously shown to be difficult to model. For these RNAs, 93% of accepted canonical base pairs were recovered in SHAPE-directed models. Discrepancies between accepted and modeled structures were small and appear to reflect genuine structural differences. Three-reagent SHAPE-directed modeling scales concisely to structurally complex RNAs to resolve the in-solution secondary structure analysis problem for many classes of RNA.

RNA families in Epstein-Barr virus

Epstein-Barr virus (EBV) is a tumorigenic human γ-herpesvirus, which produces several known structured RNAs with functional importance: two are implicated in latency maintenance and tumorigenic phenotypes, EBER1 and EBER2; a viral small nucleolar RNA (v-snoRNA1) that may generate a small regulatory RNA; and an internal ribosomal entry site in the EBNA1 mRNA. A recent bioinformatics and RNA-Seq study of EBV identified two novel EBV non-coding (nc)RNAs with evolutionary conservation in lymphocryptoviruses and likely functional importance. Both RNAs are transcribed from a repetitive region of the EBV genome (the W repeats) during a highly oncogenic type of viral latency. One novel ncRNA can form a massive (586 nt) hairpin, while the other RNA is generated from a short (81 nt) intron and is found in high abundance in EBV-infected cells.

Possible utilization of -1 Ribosomal frame shifting in the expression of a human SEMA6C isoform

We have used bioinformatics approaches to identify a potential case of -1 ribosomal frame shifting in the mRNAs of the three variants of human SEMA6C protein. The mRNAs contain a heptanucleotide slippery sequence followed by a compact H-type pseudoknot. Unlike -1 frameshifting signals in viral or viral-like mRNAs, the slippery sequence and downstream pseudoknot in SEMA6C mRNAs locate 423 nucleotides (encoding 141 amino acids) upstream of the stop codon. The potential -1 frameshifting event would produce a polypeptide of 238 residues encoded by the -1 reading frames. Sequence similarity searches using BLAST indicate that ~90% of the 238 residues match actual protein sequences annotated as SEMA6C proteins in the database. We propose that the mRNAs of human SEMA6C utilize a pseudoknot dependent -1 ribosomal frameshifting mechanism to express novel SEMA6C isoforms.

PKscan: a program to identify H-type RNA pseudoknots in any RNA sequence with unlimited length

A computer program written in C++ has been developed which can detect all potential H-type RNA pseudoknots within any given RNA sequence. There is no limit on the length of the input sequence. A validation run of the program using the full-length (8173 nt) genomic mRNA of simian retrovirus type-1 (SRV-1) identifies the established -1 frameshift stimulating pseudokont at the gagpro junction as the most stable pseudoknot within the genomic mRNA.

Specific contacts between protein S4 and ribosomal RNA are required at multiple stages of ribosome assembly

Assembly of bacterial 30S ribosomal subunits requires structural rearrangements to both its 16S rRNA and ribosomal protein components. Ribosomal protein S4 nucleates 30S assembly and associates rapidly with the 5' domain of the 16S rRNA. In vitro, transformation of initial S4-rRNA complexes to long-lived, mature complexes involves refolding of 16S helix 18, which forms part of the decoding center. Here we use targeted mutagenesis of Geobacillus stearothermophilus S4 to show that remodeling of S4-rRNA complexes is perturbed by ram alleles associated with reduced translational accuracy. Gel mobility shift assays, SHAPE chemical probing, and in vivo complementation show that the S4 N-terminal extension is required for RNA binding and viability. Alanine substitutions in Y47 and L51 that interact with 16S helix 18 decrease S4 affinity and destabilize the helix 18 pseudoknot. These changes to the protein-RNA interface correlate with no growth (L51A) or cold-sensitive growth, 30S assembly defects, and accumulation of 17S pre-rRNA (Y47A). A third mutation, R200A, over-stabilizes the helix 18 pseudoknot yet results in temperature-sensitive growth, indicating that complex stability is finely tuned by natural selection. Our results show that early S4-RNA interactions guide rRNA folding and impact late steps of 30S assembly.

Retrotransposon Ty1 RNA contains a 5'-terminal long-range pseudoknot required for efficient reverse transcription

Ty1 retrotransposon RNA has the potential to fold into a variety of distinct structures, mutation of which affects retrotransposition frequencies. We show here that one potential functional structure is located at the 5' end of the genome and can assume a pseudoknot conformation. Chemoenzymatic probing of wild-type and mutant mini-Ty1 RNAs supports the existence of such a structure, while molecular genetic analyses show that mutations disrupting pseudoknot formation interfere with retrotransposition, indicating that it provides a critical biological function. These defects are enhanced at higher temperatures. When these mutants are combined with compensatory changes, retrotransposition is restored, consistent with pseudoknot architecture. Analyses of mutants suggest a defect in Ty1 reverse transcription. Collectively, our data allow modeling of a three-dimensional structure for this novel critical cis-acting signal of the Ty1 genome.

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study

Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 μs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the μs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the μs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-μs time scale, which could easily lead to inappropriate interpretation of the simulation data.

The influenza A segment 7 mRNA 3' splice site pseudoknot/hairpin family

The 3′ splice site of the influenza A segment 7 transcript is utilized to produce mRNA for the critical M2 ion-channel protein. In solution a 63 nt fragment that includes this region can adopt two conformations: a pseudoknot and a hairpin. In each conformation, the splice site, a binding site for the SF2/ASF exonic splicing enhancer and a polypyrimidine tract, each exists in a different structural context. The most dramatic difference occurs for the splice site. In the hairpin the splice site is between two residues that are involved in a 2 by 2 nucleotide internal loop. In the pseudoknot, however, these bases are canonically paired within one of the pseudoknotted helices. The conformational switching observed in this region has implications for the regulation of splicing of the segment 7 mRNA. A measure of stability of the structures also shows interesting trends with respect to host specificity: avian strains tend to be the most stable, followed by swine and then human.