High-throughput mutate-map-rescue evaluates SHAPE-directed RNA structure and uncovers excited states

Structural and biophysical dissection of RNA conformational ensembles

RNA's ability to form and interconvert between multiple secondary and tertiary structures is critical to its functional versatility and the traditional view of RNA structures as static entities has shifted towards understanding them as dynamic conformational ensembles. In this review we discuss RNA structural ensembles and their dynamics, highlighting the concept of conformational energy landscapes as a unifying framework for understanding RNA processes such as folding, misfolding, conformational changes, and complex formation. Ongoing advancements in cryo-electron microscopy and chemical probing techniques are significantly enhancing our ability to investigate multiple structures adopted by conformationally dynamic RNAs, while traditional methods such as nuclear magnetic resonance spectroscopy continue to play a crucial role in providing high-resolution, quantitative spatial and temporal information. We discuss how these methods, when used synergistically, can provide a comprehensive understanding of RNA conformational ensembles, offering new insights into their regulatory functions.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

With hundreds of copies of rDNA, it is unknown whether they possess sequence variations that form different types of ribosomes. Here, we developed an algorithm for long-read variant calling, termed RGA, which revealed that variations in human rDNA loci are predominantly insertion-deletion (indel) variants. We developed full-length rRNA sequencing (RIBO-RT) and in situ sequencing (SWITCH-seq), which showed that translating ribosomes possess variation in rRNA. Over 1,000 variants are lowly expressed. However, tens of variants are abundant and form distinct rRNA subtypes with different structures near indels as revealed by long-read rRNA structure probing coupled to dimethyl sulfate sequencing. rRNA subtypes show differential expression in endoderm/ectoderm-derived tissues, and in cancer, low-abundance rRNA variants can become highly expressed. Together, this study identifies the diversity of ribosomes at the level of rRNA variants, their chromosomal location, and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms

Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.

Structural and computational studies of HIV-1 RNA

Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.

Dissecting the energetic architecture within an RNA tertiary structural motif via high-throughput thermodynamic measurements

Structured RNAs and RNA/protein complexes perform critical cellular functions. They often contain structurally conserved tertiary contact "motifs," whose occurrence simplifies the RNA folding landscape. Prior studies have focused on the conformational and energetic modularity of intact motifs. Here, we turn to the dissection of one common motif, the 11nt receptor (11ntR), using quantitative analysis of RNA on a massively parallel array to measure the binding of all single and double 11ntR mutants to GAAA and GUAA tetraloops, thereby probing the energetic architecture of the motif. While the 11ntR behaves as a motif, its cooperativity is not absolute. Instead, we uncovered a gradient from high cooperativity amongst base-paired and neighboring residues to additivity between distant residues. As expected, substitutions at residues in direct contact with the GAAA tetraloop resulted in the largest decreases to binding, and energetic penalties of mutations were substantially smaller for binding to the alternate GUAA tetraloop, which lacks tertiary contacts present with the canonical GAAA tetraloop. However, we found that the energetic consequences of base partner substitutions are not, in general, simply described by base pair type or isostericity. We also found exceptions to the previously established stability-abundance relationship for 11ntR sequence variants. These findings of "exceptions to the rule" highlight the power of systematic high-throughput approaches to uncover novel variants for future study in addition to providing an energetic map of a functional RNA.

Programmable antivirals targeting critical conserved viral RNA secondary structures from influenza A virus and SARS-CoV-2

Influenza A virus's (IAV's) frequent genetic changes challenge vaccine strategies and engender resistance to current drugs. We sought to identify conserved and essential RNA secondary structures within IAV's genome that are predicted to have greater constraints on mutation in response to therapeutic targeting. We identified and genetically validated an RNA structure (packaging stem-loop 2 (PSL2)) that mediates in vitro packaging and in vivo disease and is conserved across all known IAV isolates. A PSL2-targeting locked nucleic acid (LNA), administered 3 d after, or 14 d before, a lethal IAV inoculum provided 100% survival in mice, led to the development of strong immunity to rechallenge with a tenfold lethal inoculum, evaded attempts to select for resistance and retained full potency against neuraminidase inhibitor-resistant virus. Use of an analogous approach to target SARS-CoV-2, prophylactic administration of LNAs specific for highly conserved RNA structures in the viral genome, protected hamsters from efficient transmission of the SARS-CoV-2 USA_WA1/2020 variant. These findings highlight the potential applicability of this approach to any virus of interest via a process we term 'programmable antivirals', with implications for antiviral prophylaxis and post-exposure therapy.

Computational approaches for RNA structure ensemble deconvolution from structure probing data

RNA structure probing experiments have emerged over the last decade as a straightforward way to determine the structure of RNA molecules in a number of different contexts. Although powerful, the ability of RNA to dynamically interconvert between, and to simultaneously populate, alternative structural configurations, poses a nontrivial challenge to the interpretation of data derived from these experiments. Recent efforts aimed at developing computational methods for the reconstruction of coexisting alternative RNA conformations from structure probing data are paving the way to the study of RNA structure ensembles, even in the context of living cells. In this review, we critically discuss these methods, their limitations and possible future improvements.

Mutate-and-chemical-shift-fingerprint (MCSF) to characterize excited states in RNA using NMR spectroscopy

It is important to understand the dynamics and higher energy structures of RNA, called excited states, to achieve better understanding of RNA function. R1ρ relaxation dispersion NMR spectroscopy (RD) determines chemical shift differences between the most stable, ground state and the short-lived, low-populated excited states. We describe a procedure for deducing the excited state structure from these chemical shift differences using the mutate-and-chemical-shift-fingerprint (MCSF) method, which requires ~2-6 weeks and moderate understanding of NMR and RNA structure. We recently applied the MCSF methodology to elucidate the excited state of microRNA 34a targeting the SIRT1 mRNA and use this example to demonstrate the analysis. The protocol comprises the following steps: (i) determination of the secondary structure of the excited state from RD chemical shift data, (ii) design of trapped excited state RNA, (iii) validation of the excited state structure by NMR, and (iv) MCSF analysis comparing the chemical shifts of the trapped excited state with the RD-derived chemical shift differences. MCSF enables observation of the short-lived RNA structures, which can be functionally and structurally characterized by entrapment.

U2 snRNA structure is influenced by SF3A and SF3B proteins but not by SF3B inhibitors

U2 snRNP is an essential component of the spliceosome. It is responsible for branch point recognition in the spliceosome A-complex via base-pairing of U2 snRNA with an intron to form the branch helix. Small molecule inhibitors target the SF3B component of the U2 snRNP and interfere with A-complex formation during spliceosome assembly. We previously found that the first SF3B inhibited-complex is less stable than A-complex and hypothesized that SF3B inhibitors interfere with U2 snRNA secondary structure changes required to form the branch helix. Using RNA chemical modifiers, we probed U2 snRNA structure in A-complex and SF3B inhibited splicing complexes. The reactivity pattern for U2 snRNA in the SF3B inhibited-complex is indistinguishable from that of A-complex, suggesting that they have the same secondary structure conformation, including the branch helix. This observation suggests SF3B inhibited-complex instability does not stem from an alternate RNA conformation and instead points to the inhibitors interfering with protein component interactions that normally stabilize U2 snRNP’s association with an intron. In addition, we probed U2 snRNA in the free U2 snRNP in the presence of SF3B inhibitor and again saw no differences. However, increased protection of nucleotides upstream of Stem I in the absence of SF3A and SF3B proteins suggests a change of secondary structure at the very 5′ end of U2 snRNA. Chemical probing of synthetic U2 snRNA in the absence of proteins results in similar protections and predicts a previously uncharacterized extension of Stem I. Because this stem must be disrupted for SF3A and SF3B proteins to stably join the snRNP, the structure has the potential to influence snRNP assembly and recycling after spliceosome disassembly.

2'-O-Methylation can increase the abundance and lifetime of alternative RNA conformational states

2'-O-Methyl (Nm) is a highly abundant post-transcriptional RNA modification that plays important biological roles through mechanisms that are not entirely understood. There is evidence that Nm can alter the biological activities of RNAs by biasing the ribose sugar pucker equilibrium toward the C3'-endo conformation formed in canonical duplexes. However, little is known about how Nm might more broadly alter the dynamic ensembles of flexible RNAs containing bulges and internal loops. Here, using NMR and the HIV-1 transactivation response (TAR) element as a model system, we show that Nm preferentially stabilizes alternative secondary structures in which the Nm-modified nucleotides are paired, increasing both the abundance and lifetime of low-populated short-lived excited states by up to 10-fold. The extent of stabilization increased with number of Nm modifications and was also dependent on Mg2+. Through phi-value analysis, the Nm modification also provided rare insights into the structure of the transition state for conformational exchange. Our results suggest that Nm could alter the biological activities of Nm-modified RNAs by modulating their secondary structural ensembles as well as establish the utility of Nm as a tool for the discovery and characterization of RNA excited state conformations.

Folding heterogeneity in the essential human telomerase RNA three-way junction

Telomeres safeguard the genome by suppressing illicit DNA damage responses at chromosome termini. To compensate for incomplete DNA replication at telomeres, most continually dividing cells, including many cancers, express the telomerase ribonucleoprotein (RNP) complex. Telomerase maintains telomere length by catalyzing de novo synthesis of short DNA repeats using an internal telomerase RNA (TR) template. TRs from diverse species harbor structurally conserved domains that contribute to RNP biogenesis and function. In vertebrate TRs, the conserved regions 4 and 5 (CR4/5) fold into a three-way junction (TWJ) that binds directly to the telomerase catalytic protein subunit and is required for telomerase function. We have analyzed the structural properties of the human TR (hTR) CR4/5 domain using a combination of in vitro chemical mapping, secondary structural modeling, and single-molecule structural analysis. Our data suggest the essential P6.1 stem-loop within CR4/5 is not stably folded in the absence of the telomerase reverse transcriptase in vitro. Rather, the hTR CR4/5 domain adopts a heterogeneous ensemble of conformations. Finally, single-molecule FRET measurements of CR4/5 and a mutant designed to stabilize the P6.1 stem demonstrate that TERT binding selects for a structural conformation of CR4/5 that is not the dominant state of the TERT-free in vitro RNA ensemble.

Accelerated cryo-EM-guided determination of three-dimensional RNA-only structures

The discovery and design of biologically important RNA molecules is outpacing three-dimensional structural characterization. Here, we demonstrate that cryo-electron microscopy can routinely resolve maps of RNA-only systems and that these maps enable subnanometer-resolution coordinate estimation when complemented with multidimensional chemical mapping and Rosetta DRRAFTER computational modeling. This hybrid 'Ribosolve' pipeline detects and falsifies homologies and conformational rearrangements in 11 previously unknown 119- to 338-nucleotide protein-free RNA structures: full-length Tetrahymena ribozyme, hc16 ligase with and without substrate, full-length Vibrio cholerae and Fusobacterium nucleatum glycine riboswitch aptamers with and without glycine, Mycobacterium SAM-IV riboswitch with and without S-adenosylmethionine, and the computer-designed ATP-TTR-3 aptamer with and without AMP. Simulation benchmarks, blind challenges, compensatory mutagenesis, cross-RNA homologies and internal controls demonstrate that Ribosolve can accurately resolve the global architectures of RNA molecules but does not resolve atomic details. These tests offer guidelines for making inferences in future RNA structural studies with similarly accelerated throughput.

Comprehensive sequence-to-function mapping of cofactor-dependent RNA catalysis in the glmS ribozyme

Massively parallel, quantitative measurements of biomolecular activity across sequence space can greatly expand our understanding of RNA sequence-function relationships. We report the development of an RNA-array assay to perform such measurements and its application to a model RNA: the core glmS ribozyme riboswitch, which performs a ligand-dependent self-cleavage reaction. We measure the cleavage rates for all possible single and double mutants of this ribozyme across a series of ligand concentrations, determining kcat and KM values for active variants. These systematic measurements suggest that evolutionary conservation in the consensus sequence is driven by maintenance of the cleavage rate. Analysis of double-mutant rates and associated mutational interactions produces a structural and functional mapping of the ribozyme sequence, revealing the catalytic consequences of specific tertiary interactions, and allowing us to infer structural rearrangements that permit certain sequence variants to maintain activity.

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

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

The roles of structural dynamics in the cellular functions of RNAs

RNAs fold into 3D structures that range from simple helical elements to complex tertiary structures and quaternary ribonucleoprotein assemblies. The functions of many regulatory RNAs depend on how their 3D structure changes in response to a diverse array of cellular conditions. In this Review, we examine how the structural characterization of RNA as dynamic ensembles of conformations, which form with different probabilities and at different timescales, is improving our understanding of RNA function in cells. We discuss the mechanisms of gene regulation by microRNAs, riboswitches, ribozymes, post-transcriptional RNA modifications and RNA-binding proteins, and how the cellular environment and processes such as liquid-liquid phase separation may affect RNA folding and activity. The emerging RNA-ensemble-function paradigm is changing our perspective and understanding of RNA regulation, from in vitro to in vivo and from descriptive to predictive.

RNA base-pairing complexity in living cells visualized by correlated chemical probing

RNA structure and dynamics are critical to biological function. However, strategies for determining RNA structure in vivo are limited, with established chemical probing and newer duplex detection methods each having deficiencies. Here we convert the common reagent dimethyl sulfate into a useful probe of all 4 RNA nucleotides. Building on this advance, we introduce PAIR-MaP, which uses single-molecule correlated chemical probing to directly detect base-pairing interactions in cells. PAIR-MaP has superior resolution compared to alternative experiments, can resolve multiple sets of pairing interactions for structurally dynamic RNAs, and enables highly accurate structure modeling, including of RNAs containing multiple pseudoknots and extensively bound by proteins. Application of PAIR-MaP to human RNase MRP and 2 bacterial messenger RNA 5' untranslated regions reveals functionally important and complex structures undetected by prior analyses. PAIR-MaP is a powerful, experimentally concise, and broadly applicable strategy for directly visualizing RNA base pairs and dynamics in cells.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

Hidden Structural Modules in a Cooperative RNA Folding Transition

Large-scale, cooperative rearrangements underlie the functions of RNA in RNA-protein machines and gene regulation. To understand how such rearrangements are orchestrated, we used high-throughput chemical footprinting to dissect a seemingly concerted rearrangement in P5abc RNA, a paradigm of RNA folding studies. With mutations that systematically disrupt or restore putative structural elements, we found that this transition reflects local folding of structural modules, with modest and incremental cooperativity that results in concerted behavior. First, two distant secondary structure changes are coupled through a bridging three-way junction and Mg²⁺-dependent tertiary structure. Second, long-range contacts are formed between modules, resulting in additional cooperativity. Given the sparseness of RNA tertiary contacts after secondary structure formation, we expect that modular folding and incremental cooperativity are generally important for specifying functional structures while also providing productive kinetic paths to these structures. Additionally, we expect our approach to be useful for uncovering modularity in other complex RNAs.

In silico survey of the central conserved regions in viroids of the Pospiviroidae family for conserved asymmetric loop structures

Viroids are the smallest replicative pathogens, consisting of RNA circles (∼300 nucleotides) that require host machinery to replicate. Structural RNA elements recruit these host factors. Currently, many of these structural elements and the nature of their interactions are unknown. All Pospiviroidae have homology in the central conserved region (CCR). The CCR of potato spindle tuber viroid (PSTVd) contains a sarcin/ricin domain (SRD), the only viroid structural element with an unequivocal replication role. We assumed that every member of this family uses this region to recruit host factors, and that each CCR has an SRD-like asymmetric loop within it. Potential SRD or SRD-like motifs were sought in the CCR of each Pospiviroidae member as follows. Motif location in each CCR was predicted with MUSCLE alignment and Vienna RNAfold. Viroid-specific models of SRD-like motifs were built by superimposing noncanonical base pairs and nucleotides on a model of an SRD. The RNA geometry search engine FR3D was then used to find nucleotide groups close to the geometry suggested by this superimposition. Atomic resolution structures were assembled using the molecular visualization program Chimera, and the stability of each motif was assessed with molecular dynamics (MD). Some models required a protonated cytosine. To be stable within a cell, the pK_a of that cytosine must be shifted up. Constant pH-replica exchange MD analysis showed such a shift in the proposed structures. These data show that every Pospiviroidae member could form a motif that resembles an SRD in its CCR, and imply there could be undiscovered mimics of other RNA domains.

Dynamic ensemble of HIV-1 RRE stem IIB reveals non-native conformations that disrupt the Rev-binding site

The HIV-1 Rev response element (RRE) RNA element mediates the nuclear export of intron containing viral RNAs by forming an oligomeric complex with the viral protein Rev. Stem IIB and nearby stem II three-way junction nucleate oligomerization through cooperative binding of two Rev molecules. Conformational flexibility at this RRE region has been shown to be important for Rev binding. However, the nature of the flexibility has remained elusive. Here, using NMR relaxation dispersion, including a new strategy for directly observing transient conformational states in large RNAs, we find that stem IIB alone or when part of the larger RREII three-way junction robustly exists in dynamic equilibrium with non-native excited state (ES) conformations that have a combined population of ∼20%. The ESs disrupt the Rev-binding site by changing local secondary structure, and their stabilization via point substitution mutations decreases the binding affinity to the Rev arginine-rich motif (ARM) by 15- to 80-fold. The ensemble clarifies the conformational flexibility observed in stem IIB, reveals long-range conformational coupling between stem IIB and the three-way junction that may play roles in cooperative Rev binding, and also identifies non-native RRE conformational states as new targets for the development of anti-HIV therapeutics.

High-throughput determination of RNA structures

RNA performs and regulates a diverse range of cellular processes, with new functional roles being uncovered at a rapid pace. Interest is growing in how these functions are linked to RNA structures that form in the complex cellular environment. A growing suite of technologies that use advances in RNA structural probes, high-throughput sequencing and new computational approaches to interrogate RNA structure at unprecedented throughput are beginning to provide insights into RNA structures at new spatial, temporal and cellular scales.

The interplay between molecular flexibility and RNA chemical probing reactivities analyzed at the nucleotide level via an extensive molecular dynamics study

Determination of the tridimensional structure of ribonucleic acid molecules is fundamental for understanding their function in the cell. A common method to investigate RNA structures of large molecules is the use of chemical probes such as SHAPE (2^'-hydroxyl acylation analyzed by primer extension) reagents, DMS (dimethyl sulfate) and CMCT (1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfate), the reaction of which is dependent on the local structural properties of each nucleotide. In order to understand the interplay between local flexibility, sugar pucker, canonical pairing and chemical reactivity of the probes, we performed all-atom molecular dynamics simulations on a set of RNA molecules for which both tridimensional structure and chemical probing data are available and we analyzed the correlations between geometrical parameters and the chemical reactivity. Our study confirms that SHAPE reactivity is guided by the local flexibility of the different chemical moieties but suggests that a combination of multiple parameters is needed to better understand the implications of the reactivity at the molecular level. This is also the case for DMS and CMCT for which the reactivity appears to be more complex than commonly accepted.

The evolution of RNA structural probing methods: From gels to next-generation sequencing

RNA molecules are important players in all domains of life and the study of the relationship between their multiple flexible states and the associated biological roles has increased in recent years. For several decades, chemical and enzymatic structural probing experiments have been used to determine RNA structure. During this time, there has been a steady improvement in probing reagents and experimental methods, and today the structural biologist community has a large range of tools at its disposal to probe the secondary structure of RNAs in vitro and in cells. Early experiments used radioactive labeling and polyacrylamide gel electrophoresis as read-out methods. This was superseded by capillary electrophoresis, and more recently by next-generation sequencing. Today, powerful structural probing methods can characterize RNA structure on a genome-wide scale. In this review, we will provide an overview of RNA structural probing methodologies from a historical and technical perspective. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Methods > RNA Analyses in vitro and In Silico RNA Methods > RNA Analyses in Cells.

Classification of RNA structure change by 'gazing' at experimental data

Motivation

Mutations (or Single Nucleotide Variants) in folded RiboNucleic Acid structures that cause local or global conformational change are riboSNitches. Predicting riboSNitches is challenging, as it requires making two, albeit related, structure predictions. The data most often used to experimentally validate riboSNitch predictions is Selective 2' Hydroxyl Acylation by Primer Extension, or SHAPE. Experimentally establishing a riboSNitch requires the quantitative comparison of two SHAPE traces: wild-type (WT) and mutant. Historically, SHAPE data was collected on electropherograms and change in structure was evaluated by 'gel gazing.' SHAPE data is now routinely collected with next generation sequencing and/or capillary sequencers. We aim to establish a classifier capable of simulating human 'gazing' by identifying features of the SHAPE profile that human experts agree 'looks' like a riboSNitch.

Results

We find strong quantitative agreement between experts when RNA scientists 'gaze' at SHAPE data and identify riboSNitches. We identify dynamic time warping and seven other features predictive of the human consensus. The classSNitch classifier reported here accurately reproduces human consensus for 167 mutant/WT comparisons with an Area Under the Curve (AUC) above 0.8. When we analyze 2019 mutant traces for 17 different RNAs, we find that features of the WT SHAPE reactivity allow us to improve thermodynamic structure predictions of riboSNitches. This is significant, as accurate RNA structural analysis and prediction is likely to become an important aspect of precision medicine.

Availability and Implementation

The classSNitch R package is freely available at http://classsnitch.r-forge.r-project.org .

Contact

alain@email.unc.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Functional 5' UTR mRNA structures in eukaryotic translation regulation and how to find them

RNA molecules can fold into intricate shapes that can provide an additional layer of control of gene expression beyond that of their sequence. In this Review, we discuss the current mechanistic understanding of structures in 5' untranslated regions (UTRs) of eukaryotic mRNAs and the emerging methodologies used to explore them. These structures may regulate cap-dependent translation initiation through helicase-mediated remodelling of RNA structures and higher-order RNA interactions, as well as cap-independent translation initiation through internal ribosome entry sites (IRESs), mRNA modifications and other specialized translation pathways. We discuss known 5' UTR RNA structures and how new structure probing technologies coupled with prospective validation, particularly compensatory mutagenesis, are likely to identify classes of structured RNA elements that shape post-transcriptional control of gene expression and the development of multicellular organisms.

Allosteric mechanism of the V. vulnificus adenine riboswitch resolved by four-dimensional chemical mapping

The structural interconversions that mediate the gene regulatory functions of RNA molecules may be different from classic models of allostery, but the relevant structural correlations have remained elusive in even intensively studied systems. Here, we present a four-dimensional expansion of chemical mapping called lock-mutate-map-rescue (LM²R), which integrates multiple layers of mutation with nucleotide-resolution chemical mapping. This technique resolves the core mechanism of the adenine-responsive V. vulnificus add riboswitch, a paradigmatic system for which both Monod-Wyman-Changeux (MWC) conformational selection models and non-MWC alternatives have been proposed. To discriminate amongst these models, we locked each functionally important helix through designed mutations and assessed formation or depletion of other helices via compensatory rescue evaluated by chemical mapping. These LM²R measurements give strong support to the pre-existing correlations predicted by MWC models, disfavor alternative models, and suggest additional structural heterogeneities that may be general across ligand-free riboswitches.

Primerize-2D: automated primer design for RNA multidimensional chemical mapping

Summary

Rapid RNA synthesis of comprehensive single mutant libraries and targeted multiple mutant libraries is enabling new multidimensional chemical approaches to solve RNA structures. PCR assembly of DNA templates and in vitro transcription allow synthesis and purification of hundreds of RNA mutants in a cost-effective manner, with sharing of primers across constructs allowing significant reductions in expense. However, these protocols require organization of primer locations across numerous 96 well plates and guidance for pipetting, non-trivial tasks for which informatics and visualization tools can prevent costly errors. We report here an online tool to accelerate synthesis of large libraries of desired mutants through design and efficient organization of primers. The underlying program and graphical interface have been experimentally tested in our laboratory for RNA domains with lengths up to 300 nucleotides and libraries encompassing up to 960 variants. In addition to the freely available Primerize-2D server, the primer design code is available as a stand-alone Python package for broader applications.

Availability and Implementation

http://primerize2d.stanford.edu.

Contact

rhiju@stanford.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Molecular Dynamics Simulations Reveal an Interplay between SHAPE Reagent Binding and RNA Flexibility

The function of RNA molecules usually depends on their overall fold and on the presence of specific structural motifs. Chemical probing methods are routinely used in combination with nearest-neighbor models to determine RNA secondary structure. Among the available methods, SHAPE is relevant due to its capability to probe all RNA nucleotides and the possibility to be used in vivo. However, the structural determinants for SHAPE reactivity and its mechanism of reaction are still unclear. Here molecular dynamics simulations and enhanced sampling techniques are used to predict the accessibility of nucleotide analogs and larger RNA structural motifs to SHAPE reagents. We show that local RNA reconformations are crucial in allowing reagents to reach the 2'-OH group of a particular nucleotide and that sugar pucker is a major structural factor influencing SHAPE reactivity.

RNA structure inference through chemical mapping after accidental or intentional mutations

Despite the critical roles RNA structures play in regulating gene expression, sequencing-based methods for experimentally determining RNA base pairs have remained inaccurate. Here, we describe a multidimensional chemical-mapping method called "mutate-and-map read out through next-generation sequencing" (M2-seq) that takes advantage of sparsely mutated nucleotides to induce structural perturbations at partner nucleotides and then detects these events through dimethyl sulfate (DMS) probing and mutational profiling. In special cases, fortuitous errors introduced during DNA template preparation and RNA transcription are sufficient to give M2-seq helix signatures; these signals were previously overlooked or mistaken for correlated double-DMS events. When mutations are enhanced through error-prone PCR, in vitro M2-seq experimentally resolves 33 of 68 helices in diverse structured RNAs including ribozyme domains, riboswitch aptamers, and viral RNA domains with a single false positive. These inferences do not require energy minimization algorithms and can be made by either direct visual inspection or by a neural-network-inspired algorithm called M2-net. Measurements on the P4-P6 domain of the Tetrahymena group I ribozyme embedded in Xenopus egg extract demonstrate the ability of M2-seq to detect RNA helices in a complex biological environment.

Resolving sugar puckers in RNA excited states exposes slow modes of repuckering dynamics

Recent studies have shown that RNAs exist in dynamic equilibrium with short-lived low-abundance 'excited states' that form by reshuffling base pairs in and around non-canonical motifs. These conformational states are proposed to be rich in non-canonical motifs and to play roles in the folding and regulatory functions of non-coding RNAs but their structure proves difficult to characterize given their transient nature. Here, we describe an approach for determining sugar pucker conformation in RNA excited states through nuclear magnetic resonance measurements of C1΄ and C4΄ rotating frame spin relaxation (R1ρ) in uniformly 13C/15N labeled RNA samples. Application to HIV-1 TAR exposed slow modes of sugar repuckering dynamics at the μs and ms timescale accompanying transitions between non-helical (C2΄-endo) to helical (C3΄-endo) conformations during formation of two distinct excited states. In contrast, we did not obtain any evidence for slow sugar repuckering dynamics for nucleotides in a variety of structural contexts that do not undergo non-helical to helical transitions. Our results outline a route for significantly improving the conformational characterization of RNA excited states and suggest that slow modes of repuckering dynamics gated by transient changes in secondary structure are quite common in RNA.

RNA structure through multidimensional chemical mapping

The discoveries of myriad non-coding RNA molecules, each transiting through multiple flexible states in cells or virions, present major challenges for structure determination. Advances in high-throughput chemical mapping give new routes for characterizing entire transcriptomes in vivo, but the resulting one-dimensional data generally remain too information-poor to allow accurate de novo structure determination. Multidimensional chemical mapping (MCM) methods seek to address this challenge. Mutate-and-map (M2), RNA interaction groups by mutational profiling (RING-MaP and MaP-2D analysis) and multiplexed •OH cleavage analysis (MOHCA) measure how the chemical reactivities of every nucleotide in an RNA molecule change in response to modifications at every other nucleotide. A growing body of in vitro blind tests and compensatory mutation/rescue experiments indicate that MCM methods give consistently accurate secondary structures and global tertiary structures for ribozymes, ribosomal domains and ligand-bound riboswitch aptamers up to 200 nucleotides in length. Importantly, MCM analyses provide detailed information on structurally heterogeneous RNA states, such as ligand-free riboswitches that are functionally important but difficult to resolve with other approaches. The sequencing requirements of currently available MCM protocols scale at least quadratically with RNA length, precluding general application to transcriptomes or viral genomes at present. We propose a modify-cross-link-map (MXM) expansion to overcome this and other current limitations to resolving the in vivo 'RNA structurome'.

A critical base pair in k-turns determines the conformational class adopted, and correlates with biological function

k-turns are commonly-occurring motifs that introduce sharp kinks into duplex RNA, thereby facilitating tertiary contacts. Both the folding and conformation of k-turns are determined by their local sequence. k-turns fall into two conformational classes, called N3 and N1, that differ in the pattern of hydrogen bonding in the core. We show here that this is determined by the basepair adjacent to the critical G•A pairs. We determined crystal structures of a series of Kt-7 variants in which this 3b,3n position has been systematically varied, showing that this leads to a switch in the conformation. We have previously shown that the 3b,3n position also determines the folding characteristics of the k-turn, i.e. whether or not the k-turn can fold in the presence of metal ions alone. We have analyzed the distribution of 3b,3n sequences from four classes of k-turns from ribosomes, riboswitches and U4 snRNA, finding a strong conservation of properties for a given k-turn type. We thus demonstrate a strong association between biological function, 3b,3n sequence and k-turn folding and conformation. This has strong predictive power, and can be applied to the modeling of large RNA architectures.

Nucleic acid polymeric properties and electrostatics: Directly comparing theory and simulation with experiment

Nucleic acids are biopolymers that carry genetic information and are also involved in various gene regulation functions such as gene silencing and protein translation. Because of their negatively charged backbones, nucleic acids are polyelectrolytes. To adequately understand nucleic acid folding and function, we need to properly describe its i) polymer/polyelectrolyte properties and ii) associating ion atmosphere. While various theories and simulation models have been developed to describe nucleic acids and the ions around them, many of these theories/simulations have not been well evaluated due to complexities in comparison with experiment. In this review, I discuss some recent experiments that have been strategically designed for straightforward comparison with theories and simulation models. Such data serve as excellent benchmarks to identify limitations in prevailing theories and simulation parameters.

RNAComposer and RNA 3D structure prediction for nanotechnology

RNAs adopt specific, stable tertiary architectures to perform their activities. Knowledge of RNA tertiary structure is fundamental to understand RNA functions beginning with transcription and ending with turnover. Contrary to advanced RNA secondary structure prediction algorithms, which allow good accuracy when experimental data are integrated into the prediction, tertiary structure prediction of large RNAs still remains a significant challenge. However, the field of RNA tertiary structure prediction is rapidly developing and new computational methods based on different strategies are emerging. RNAComposer is a user-friendly and freely available server for 3D structure prediction of RNA up to 500 nucleotide residues. RNAComposer employs fully automated fragment assembly based on RNA secondary structure specified by the user. Importantly, this method allows incorporation of distance restraints derived from the experimental data to strengthen the 3D predictions. The potential and limitations of RNAComposer are discussed and an application to RNA design for nanotechnology is presented.

Automated RNA 3D Structure Prediction with RNAComposer

RNAs adopt specific structures to perform their activities and these are critical to virtually all RNA-mediated processes. Because of difficulties in experimentally assessing structures of large RNAs using NMR, X-ray crystallography, or cryo-microscopy, there is currently great demand for new high-resolution 3D structure prediction methods. Recently we reported on RNAComposer, a knowledge-based method for the fully automated RNA 3D structure prediction from a user-defined secondary structure. RNAComposer method is especially suited for structural biology users. Since our initial report in 2012, both servers, freely available at http://rnacomposer.ibch.poznan.pl and http://rnacomposer.cs.put.poznan.pl have been often visited. Therefore this chapter provides guidance for using RNAComposer and discusses points that should be considered when predicting 3D RNA structure. An application example presents current scope and limitations of RNAComposer.

Rich RNA Structure Landscapes Revealed by Mutate-and-Map Analysis

Landscapes exhibiting multiple secondary structures arise in natural RNA molecules that modulate gene expression, protein synthesis, and viral infection [corrected]. We report herein that high-throughput chemical experiments can isolate an RNA's multiple alternative secondary structures as they are stabilized by systematic mutagenesis (mutate-and-map, M2) and that a computational algorithm, REEFFIT, enables unbiased reconstruction of these states' structures and populations. In an in silico benchmark on non-coding RNAs with complex landscapes, M2-REEFFIT recovers 95% of RNA helices present with at least 25% population while maintaining a low false discovery rate (10%) and conservative error estimates. In experimental benchmarks, M2-REEFFIT recovers the structure landscapes of a 35-nt MedLoop hairpin, a 110-nt 16S rRNA four-way junction with an excited state, a 25-nt bistable hairpin, and a 112-nt three-state adenine riboswitch with its expression platform, molecules whose characterization previously required expert mutational analysis and specialized NMR or chemical mapping experiments. With this validation, M2-REEFFIT enabled tests of whether artificial RNA sequences might exhibit complex landscapes in the absence of explicit design. An artificial flavin mononucleotide riboswitch and a randomly generated RNA sequence are found to interconvert between three or more states, including structures for which there was no design, but that could be stabilized through mutations. These results highlight the likely pervasiveness of rich landscapes with multiple secondary structures in both natural and artificial RNAs and demonstrate an automated chemical/computational route for their empirical characterization.

Consistent global structures of complex RNA states through multidimensional chemical mapping

Accelerating discoveries of non-coding RNA (ncRNA) in myriad biological processes pose major challenges to structural and functional analysis. Despite progress in secondary structure modeling, high-throughput methods have generally failed to determine ncRNA tertiary structures, even at the 1-nm resolution that enables visualization of how helices and functional motifs are positioned in three dimensions. We report that integrating a new method called MOHCA-seq (Multiplexed •OHCleavage Analysis with paired-end sequencing) with mutate-and-map secondary structure inference guides Rosetta 3D modeling to consistent 1-nm accuracy for intricately folded ncRNAs with lengths up to 188 nucleotides, including a blind RNA-puzzle challenge, the lariat-capping ribozyme. This multidimensional chemical mapping (MCM) pipeline resolves unexpected tertiary proximities for cyclic-di-GMP, glycine, and adenosylcobalamin riboswitch aptamers without their ligands and a loose structure for the recently discovered human HoxA9D internal ribosome entry site regulon. MCM offers a sequencing-based route to uncovering ncRNA 3D structure, applicable to functionally important but potentially heterogeneous states.

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

Our genetic material, in the form of molecules of DNA, provides instructions for many different processes in our cells. To issue these instructions, particular sections of DNA are copied to make a type of molecule called ribonucleic acid (RNA). Some of these RNA molecules contain instructions to make proteins, but others—known as non-coding RNAs—regulate the activity of genes in cells.

The genetic information within RNA is encoded by the sequence of four different chemical parts called ‘nucleotides’. RNA can exist as a single strand of nucleotides, but the nucleotides can also pair up in specific combinations to form sections of double-stranded RNA. Therefore, a single strand of non-coding RNA can fold into a complex three-dimensional shape that contains loops, twists, and bulges.

The three-dimensional structures of non-coding RNAs are crucial for their roles in cells, but the variety and complexity of shapes that they can form makes it technically difficult to study them. In 2008, researchers developed a new method called MOHCA that can map the positions of nucleotides that are close together in the three-dimensional structure. Highly reactive chemicals are attached to the nucleotides and these can react with, and damage, other nearby nucleotides. By detecting which nucleotides have been damaged, it is possible to map the positions of these nucleotides and decipher the structure of the RNA molecule using computer algorithms.

MOHCA is a promising approach, but the initial methods to find the damaged nucleotides were tedious and required specialized equipment. Now, Cheng, Das et al.—including some of the researchers involved in the 2008 work—have developed an improved version of MOHCA that uses readily available RNA sequencing techniques to find the damaged nucleotides. The RNA sequencing data are then analyzed by a new algorithm in the Rosetta computer modeling software.

Cheng, Das et al. used this newly developed ‘MOHCA-seq’ and Rosetta to reveal the structures of a human non-coding RNA and several other non-coding RNA molecules to a much higher level of detail than before. Together, MOHCA-seq and Rosetta provide a rapid method for researchers to decipher the three-dimensional structure of non-coding RNAs. This method is likely to speed up the analysis of the complex structures of non-coding RNAs. It will be useful in future efforts to work out what roles these RNAs play in cells, including their activity in cancer, neurodegeneration, and other diseases.

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

RNA-Puzzles Round II: assessment of RNA structure prediction programs applied to three large RNA structures

This paper is a report of a second round of RNA-Puzzles, a collective and blind experiment in three-dimensional (3D) RNA structure prediction. Three puzzles, Puzzles 5, 6, and 10, represented sequences of three large RNA structures with limited or no homology with previously solved RNA molecules. A lariat-capping ribozyme, as well as riboswitches complexed to adenosylcobalamin and tRNA, were predicted by seven groups using RNAComposer, ModeRNA/SimRNA, Vfold, Rosetta, DMD, MC-Fold, 3dRNA, and AMBER refinement. Some groups derived models using data from state-of-the-art chemical-mapping methods (SHAPE, DMS, CMCT, and mutate-and-map). The comparisons between the predictions and the three subsequently released crystallographic structures, solved at diffraction resolutions of 2.5-3.2 Å, were carried out automatically using various sets of quality indicators. The comparisons clearly demonstrate the state of present-day de novo prediction abilities as well as the limitations of these state-of-the-art methods. All of the best prediction models have similar topologies to the native structures, which suggests that computational methods for RNA structure prediction can already provide useful structural information for biological problems. However, the prediction accuracy for non-Watson-Crick interactions, key to proper folding of RNAs, is low and some predicted models had high Clash Scores. These two difficulties point to some of the continuing bottlenecks in RNA structure prediction. All submitted models are available for download at http://ahsoka.u-strasbg.fr/rnapuzzles/.

Primerize: automated primer assembly for transcribing non-coding RNA domains

Customized RNA synthesis is in demand for biological and biotechnological research. While chemical synthesis and gel or chromatographic purification of RNA is costly and difficult for sequences longer than tens of nucleotides, a pipeline of primer assembly of DNA templates, in vitro transcription by T7 RNA polymerase and kit-based purification provides a cost-effective and fast alternative for preparing RNA molecules. Nevertheless, designing template primers that optimize cost and avoid mispriming during polymerase chain reaction currently requires expert inspection, downloading specialized software or both. Online servers are currently not available or maintained for the task. We report here a server named Primerize that makes available an efficient algorithm for primer design developed and experimentally tested in our laboratory for RNA domains with lengths up to 300 nucleotides. Free access: http://primerize.stanford.edu.

Automated band annotation for RNA structure probing experiments with numerous capillary electrophoresis profiles

MOTIVATION

Capillary electrophoresis (CE) is a powerful approach for structural analysis of nucleic acids, with recent high-throughput variants enabling three-dimensional RNA modeling and the discovery of new rules for RNA structure design. Among the steps composing CE analysis, the process of finding each band in an electrophoretic trace and mapping it to a position in the nucleic acid sequence has required significant manual inspection and remains the most time-consuming and error-prone step. The few available tools seeking to automate this band annotation have achieved limited accuracy and have not taken advantage of information across dozens of profiles routinely acquired in high-throughput measurements.

RESULTS

We present a dynamic-programming-based approach to automate band annotation for high-throughput capillary electrophoresis. The approach is uniquely able to define and optimize a robust target function that takes into account multiple CE profiles (sequencing ladders, different chemical probes, different mutants) collected for the RNA. Over a large benchmark of multi-profile datasets for biological RNAs and designed RNAs from the EteRNA project, the method outperforms prior tools (QuSHAPE and FAST) significantly in terms of accuracy compared with gold-standard manual annotations. The amount of computation required is reasonable at a few seconds per dataset. We also introduce an 'E-score' metric to automatically assess the reliability of the band annotation and show it to be practically useful in flagging uncertainties in band annotation for further inspection.

AVAILABILITY AND IMPLEMENTATION

The implementation of the proposed algorithm is included in the HiTRACE software, freely available as an online server and for download at http://hitrace.stanford.edu.

CONTACT

sryoon@snu.ac.kr or rhiju@stanford.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Characterizing RNA Excited States Using NMR Relaxation Dispersion

Changes in RNA secondary structure play fundamental roles in the cellular functions of a growing number of noncoding RNAs. This chapter describes NMR-based approaches for characterizing microsecond-to-millisecond changes in RNA secondary structure that are directed toward short-lived and low-populated species often referred to as "excited states." Compared to larger scale changes in RNA secondary structure, transitions toward excited states do not require assistance from chaperones, are often orders of magnitude faster, and are localized to a small number of nearby base pairs in and around noncanonical motifs. Here, we describe a procedure for characterizing RNA excited states using off-resonance R1ρ NMR relaxation dispersion utilizing low-to-high spin-lock fields (25-3000 Hz). R1ρ NMR relaxation dispersion experiments are used to measure carbon and nitrogen chemical shifts in base and sugar moieties of the excited state. The chemical shift data are then interpreted with the aid of secondary structure prediction to infer potential excited states that feature alternative secondary structures. Candidate structures are then tested by using mutations, single-atom substitutions, or by changing physiochemical conditions, such as pH and temperature, to either stabilize or destabilize the candidate excited state. The resulting chemical shifts of the mutants or under different physiochemical conditions are then compared to those of the ground and excited states. Application is illustrated with a focus on the transactivation response element from the human immune deficiency virus type 1, which exists in dynamic equilibrium with at least two distinct excited states.

Modeling complex RNA tertiary folds with Rosetta

Reliable modeling of RNA tertiary structures is key to both understanding these structures' roles in complex biological machines and to eventually facilitating their design for molecular computing and robotics. In recent years, a concerted effort to improve computational prediction of RNA structure through the RNA-Puzzles blind prediction trials has accelerated advances in the field. Among other approaches, the versatile and expanding Rosetta molecular modeling software now permits modeling of RNAs in the 100-300 nucleotide size range at consistent subhelical (~1 nm) resolution. Our laboratory's current state-of-the-art methods for RNAs in this size range involve Fragment Assembly of RNA with Full-Atom Refinement (FARFAR), which optimizes RNA conformations in the context of a physically realistic energy function, as well as hybrid techniques that leverage experimental data to inform computational modeling. In this chapter, we give a practical guide to our current workflow for modeling RNA three-dimensional structures using FARFAR, including strategies for using data from multidimensional chemical mapping experiments to focus sampling and select accurate conformations.

RNA regulons in Hox 5' UTRs confer ribosome specificity to gene regulation

Emerging evidence suggests that the ribosome has a regulatory function in directing how the genome is translated in time and space. However, how this regulation is encoded in the messenger RNA sequence remains largely unknown. Here we uncover unique RNA regulons embedded in homeobox (Hox) 5' untranslated regions (UTRs) that confer ribosome-mediated control of gene expression. These structured RNA elements, resembling viral internal ribosome entry sites (IRESs), are found in subsets of Hox mRNAs. They facilitate ribosome recruitment and require the ribosomal protein RPL38 for their activity. Despite numerous layers of Hox gene regulation, these IRES elements are essential for converting Hox transcripts into proteins to pattern the mammalian body plan. This specialized mode of IRES-dependent translation is enabled by an additional regulatory element that we term the translation inhibitory element (TIE), which blocks cap-dependent translation of transcripts. Together, these data uncover a new paradigm for ribosome-mediated control of gene expression and organismal development.

Probing viral genomic structure: alternative viewpoints and alternative structures for satellite tobacco mosaic virus RNA

Viral RNA structure prediction is a valuable tool for development of drugs against viral disease. This work discusses different approaches to predicting encapsidated viral RNA and highlights satellite tobacco mosaic virus (STMV) RNA as a model system with excellent crystallography data. Fundamentally important issues for debate include thermodynamic versus kinetic control of virus assembly and the possible consequences of quasi-species in the primary structure on RNA secondary structure prediction of a single structure or an ensemble of structures. Multiple computational tools and chemical reagents are now available for improved viral RNA structure prediction. Two different predicted structures for encapsidated STMV RNA result from differences in three main areas: a different approach and philosophy to studying encapsidated viral RNA, an emphasis on different RNA motifs, and technical differences in computational methods and chemical reagents. The experiments with traditional chemical probing and SHAPE reagents are compared in terms of chemistry, results, and interpretation for STMV RNA as well as other RNA protein assemblies, such as the 5'UTR of HIV and the ribosome. This discussion of the challenges of viral RNA structure prediction will lead to new experiments and improved future predictions for viral RNA.