Search

RNA molecules can show high levels of cooperativity in their global folding and interactions with divalent ions. However, cooperativity at individual ligand-RNA interaction sites remains poorly understood. Here, we investigated the binding of thiamine and methylene diphosphonic acid (MDP, a soluble structural analogue of pyrophosphate) to the thiamine pyrophosphate riboswitch. These ligands each bind weakly at proximal subsites, with 10 μM and 1 mM affinities, respectively. The affinity of MDP moderately improves when thiamine or thiamine-like fragments are pre-bound to the RNA. Covalent linking of thiamine and MDP substantially increases riboswitch binding to a notable high affinity of 20 nM. Crystal structures and single-molecule correlated chemical probing revealed favorable induced fit effects upon binding of individual ligands and, unexpectedly, a substantial thermodynamically unfavorable RNA structural rearrangement upon binding of the linked thiamine-MDP ligand. Thus, linking of two ligands of modest affinity, accompanied by an unfavorable structural rearrangement, still yields a potent linked RNA-binding compound. Since complex ligands often bind riboswitches and other RNAs at proximal subsites, principles derived from this work inform and support fragment-linking strategies for identifying small molecules that interact with RNA specifically and with high affinity.

Higher-order structure governs function for many RNAs. However, discerning this structure for large RNA molecules in solution is an unresolved challenge. Here, we present SHAPE-JuMP (selective 2'-hydroxyl acylation analyzed by primer extension and juxtaposed merged pairs) to interrogate through-space RNA tertiary interactions. A bifunctional small molecule is used to chemically link proximal nucleotides in an RNA structure. The RNA cross-link site is then encoded into complementary DNA (cDNA) in a single, direct step using an engineered reverse transcriptase that "jumps" across cross-linked nucleotides. The resulting cDNAs contain a deletion relative to the native RNA sequence, which can be detected by sequencing, that indicates the sites of cross-linked nucleotides. SHAPE-JuMP measures RNA tertiary structure proximity concisely across large RNA molecules at nanometer resolution. SHAPE-JuMP is especially effective at measuring interactions in multihelix junctions and loop-to-helix packing, enables modeling of the global fold for RNAs up to several hundred nucleotides in length, facilitates ranking of structural models by consistency with through-space restraints, and is poised to enable solution-phase structural interrogation and modeling of complex RNAs.

7SK is a conserved non-coding RNA that regulates transcription by sequestering the transcription factor P-TEFb. 7SK function entails complex changes in RNA structure, but characterizing RNA dynamics in cells remains an unsolved challenge. We developed a single-molecule chemical probing strategy, DANCE-MaP (deconvolution and annotation of ribonucleic conformational ensembles), that defines per-nucleotide reactivity, direct base pairing interactions, tertiary interactions, and thermodynamic populations for each state in RNA structural ensembles from a single experiment. DANCE-MaP reveals that 7SK RNA encodes a large-scale structural switch that couples dissolution of the P-TEFb binding site to structural remodeling at distal release factor binding sites. The 7SK structural equilibrium shifts in response to cell growth and stress and can be targeted to modulate expression of P-TEFb responsive genes. Our study reveals that RNA structural dynamics underlie 7SK function as an integrator of diverse cellular signals to control transcription, and establishes the power of DANCE-MaP to define RNA dynamics in cells.

RNA molecules adopt three-dimensional structures that are critical to their function and of interest in drug discovery. Few RNA structures are known, however, and predicting them computationally has proven challenging. We introduce a machine learning approach that enables identification of accurate structural models without assumptions about their defining characteristics, despite being trained with only 18 known RNA structures. The resulting scoring function, the Atomic Rotationally Equivariant Scorer (ARES), substantially outperforms previous methods and consistently produces the best results in community-wide blind RNA structure prediction challenges. By learning effectively even from a small amount of data, our approach overcomes a major limitation of standard deep neural networks. Because it uses only atomic coordinates as inputs and incorporates no RNA-specific information, this approach is applicable to diverse problems in structural biology, chemistry, materials science, and beyond.

Nucleocapsid protein (N-protein) is required for multiple steps in betacoronaviruses replication. SARS-CoV-2-N-protein condenses with specific viral RNAs at particular temperatures making it a powerful model for deciphering RNA sequence specificity in condensates. We identify two separate and distinct double-stranded, RNA motifs (dsRNA stickers) that promote N-protein condensation. These dsRNA stickers are separately recognized by N-protein's two RNA binding domains (RBDs). RBD1 prefers structured RNA with sequences like the transcription-regulatory sequence (TRS). RBD2 prefers long stretches of dsRNA, independent of sequence. Thus, the two N-protein RBDs interact with distinct dsRNA stickers, and these interactions impart specific droplet physical properties that could support varied viral functions. Specifically, we find that addition of dsRNA lowers the condensation temperature dependent on RBD2 interactions and tunes translational repression. In contrast RBD1 sites are sequences critical for sub-genomic (sg) RNA generation and promote gRNA compression. The density of RBD1 binding motifs in proximity to TRS-L/B sequences is associated with levels of sub-genomic RNA generation. The switch to packaging is likely mediated by RBD1 interactions which generate particles that recapitulate the packaging unit of the virion. Thus, SARS-CoV-2 can achieve biochemical complexity, performing multiple functions in the same cytoplasm, with minimal protein components based on utilizing multiple distinct RNA motifs that control N-protein interactions.

Chemical probing experiments, coupled with empirically determined free energy change relationships, can enable accurate modeling of the secondary structures of diverse and complex RNAs. A current frontier lies in modeling large and structurally heterogeneous transcripts, including complex eukaryotic RNAs. To validate and improve on experimentally driven approaches for modeling large transcripts, we obtained high-quality SHAPE data for the protein-free human 18S and 28S ribosomal RNAs (rRNAs). To our surprise, SHAPE-directed structure models for the human rRNAs poorly matched accepted structures. Analysis of predicted rRNA structures based on low-SHAPE and low-entropy (lowSS) metrics revealed that, whereas ∼75% of Escherichia coli rRNA sequences form well-determined lowSS secondary structure, only ∼40% of the human rRNAs do. Critically, regions of the human rRNAs that specifically fold into well-determined lowSS structures were modeled to high accuracy using SHAPE data. This work reveals that eukaryotic rRNAs are more unfolded than are those of prokaryotic rRNAs and indeed are largely unfolded overall, likely reflecting increased protein dependence for eukaryotic ribosome structure. In addition, those regions and substructures that are well-determined can be identified de novo and successfully modeled by SHAPE-directed folding.

Chemical probing is an important tool for characterizing the complex folded structures of RNA molecules, many of which play key cellular roles. Electrophilic SHAPE reagents create adducts at the 2'-hydroxyl position on the RNA backbone of flexible ribonucleotides with relatively little dependence on nucleotide identity. Strategies for adduct detection such as mutational profiling (MaP) allow accurate, automated calculation of relative adduct frequencies for each nucleotide in a given RNA or group of RNAs. A number of alternative reagents and adduct detection strategies have been proposed, especially for use in living cells. Here we evaluate five SHAPE reagents: three previously well-validated reagents 1M7 (1-methyl-7-nitroisatoic anhydride), 1M6 (1-methyl-6-nitroisatoic anhydride), and NMIA ( N-methylisatoic anhydride), one more recently proposed NAI (2-methylnicotinic acid imidazolide), and one novel reagent 5NIA (5-nitroisatoic anhydride). We clarify the importance of carefully designed software in reading out SHAPE experiments using massively parallel sequencing approaches. We examine SHAPE modification in living cells in diverse cell lines, compare MaP and reverse transcription-truncation as SHAPE adduct detection strategies, make recommendations for SHAPE reagent choice, and outline areas for future development.

Every class of RNA forms base-paired structures that impact biological functions. Chemical probing of RNA structure, especially with the advent of strategies such as SHAPE-MaP, vastly expands the scale and quantitative accuracy over which RNA structure can be examined. These methods have enabled large-scale structural studies of mRNAs and lncRNAs, but the length and complexity of these RNAs makes interpretation of the data challenging. We have created modules available through the open-source Integrative Genomics Viewer (IGV) for straightforward visualization of RNA structures and associated and complementary experimental data. Here we present detailed and stepwise strategies for exploring and visualizing complex RNA structures in IGV. Individuals can use these instructions and supplied sample data to become adept at using IGV to visualize RNA structure models in conjunction with useful allied information.

We report that the SARS-CoV-2 nucleocapsid protein (N-protein) undergoes liquid-liquid phase separation (LLPS) with viral RNA. N-protein condenses with specific RNA genomic elements under physiological buffer conditions and condensation is enhanced at human body temperatures (33°C and 37°C) and reduced at room temperature (22°C). RNA sequence and structure in specific genomic regions regulate N-protein condensation while other genomic regions promote condensate dissolution, potentially preventing aggregation of the large genome. At low concentrations, N-protein preferentially crosslinks to specific regions with single-stranded RNA flanked by structure and these features specify the location, number, and strength of N-protein binding sites (valency). Liquid-like N-protein condensates form in mammalian cells in a concentration-dependent manner and can be altered by small molecules. Condensation of N-protein is RNA sequence and structure specific, sensitive to human body temperature, and manipulatable with small molecules, and presents a screenable process for identifying antiviral compounds effective against SARS-CoV-2., Graphical Abstract, Highlights Phase separation occurs with the viral genome (gRNA) and at human body temperature. Phase separation is driven by specific elements in gRNA. RBD mutant N-protein fails to undergo LLPS; exhibits altered RNA crosslinking. N-protein forms liquid-like droplets in cells., Iserman and Roden et al. demonstrate phase separation (LLPS) of SARS-CoV-2 nucleocapsid (N-protein) with viral RNA. Viral RNA sequences promote or oppose phase separation depending on binding patterns of N-protein and genomic RNA. LLPS-promoting sequences occur at 5′ and 3′-Ends of the genome, suggestive of a genome packaging role.

A mechanistic understanding of the SARS-CoV-2 viral replication cycle is essential to develop new therapies for the COVID-19 global health crisis. In this study, we show that the SARS-CoV-2 nucleocapsid protein (N-protein) undergoes liquid-liquid phase separation (LLPS) with the viral genome, and propose a model of viral packaging through LLPS. N-protein condenses with specific RNA sequences in the first 1000 nts (5’-End) under physiological conditions and is enhanced at human upper airway temperatures. N-protein condensates exclude non-packaged RNA sequences. We comprehensively map sites bound by N-protein in the 5’-End and find preferences for single-stranded RNA flanked by stable structured elements. Liquid-like N-protein condensates form in mammalian cells in a concentration-dependent manner and can be altered by small molecules. Condensation of N-protein is sequence and structure specific, sensitive to human body temperature, and manipulatable with small molecules thus presenting screenable processes for identifying antiviral compounds effective against SARS-CoV-2.

RNA–protein interaction networks govern many biological processes but are difficult to examine comprehensively. We devised ribonucleoprotein networks analyzed by mutational profiling (RNP-MaP), a live-cell chemical probing strategy that maps cooperative interactions among multiple proteins bound to single RNA molecules at nucleotide resolution. RNP-MaP uses a hetero-bifunctional crosslinker to freeze interacting proteins in place on RNA and then maps multiple bound proteins on single RNA strands by read-through reverse transcription and DNA sequencing. RNP-MaP revealed that RNase P and RMRP, two sequence-divergent but structurally related non-coding RNAs, share RNP networks and that network hubs define functional sites in these RNAs. RNP-MaP also identified protein interaction networks conserved between mouse and human XIST long non-coding RNAs and defined protein communities whose binding sites colocalize and form networks in functional regions of XIST. RNP-MaP enables discovery and efficient validation of functional protein interaction networks on long RNAs in living cells. Networks of proteins bound to single RNAs are identified by correlated chemical crosslinking.

The functions of most long non-coding RNAs (lncRNAs) are unknown. In contrast to proteins, lncRNAs with similar functions often lack linear sequence homology; thus, the identification of function in one lncRNA rarely informs the identification of function in others. We developed a sequence comparison method to deconstruct linear sequence relationships in lncRNAs and evaluate similarity based on the abundance of short motifs called k-mers. We found that lncRNAs of related function often had similar k-mer profiles despite lacking linear homology, and that k-mer profiles correlated with protein binding to lncRNAs and with their subcellular localization. Using a novel assay to quantify Xist-like regulatory potential, we directly demonstrated that evolutionarily unrelated lncRNAs can encode similar function through different spatial arrangements of related sequence motifs. K-mer-based classification is a powerful approach to detect recurrent relationships between sequence and function in lncRNAs.

This protocol is an extension to: Nat. Protoc. 10, 1643-1669 (2015); doi:10.1038/nprot.2015.103; published online 01 October 2015RNAs play key roles in many cellular processes. The underlying structure of RNA is an important determinant of how transcripts function, are processed, and interact with RNA-binding proteins and ligands. RNA structure analysis by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) takes advantage of the reactivity of small electrophilic chemical probes that react with the 2'-hydroxyl group to assess RNA structure at nucleotide resolution. When coupled with mutational profiling (MaP), in which modified nucleotides are detected as internal miscodings during reverse transcription and then read out by massively parallel sequencing, SHAPE yields quantitative per-nucleotide measurements of RNA structure. Here, we provide an extension to our previous in vitro SHAPE-MaP protocol with detailed guidance for undertaking and analyzing SHAPE-MaP probing experiments in live cells. The MaP strategy works for both abundant-transcriptome experiments and for cellular RNAs of low to moderate abundance, which are not well examined by whole-transcriptome methods. In-cell SHAPE-MaP, performed in roughly 3 d, can be applied in cell types ranging from bacteria to cultured mammalian cells and is compatible with a variety of structure-probing reagents. We detail several strategies by which in-cell SHAPE-MaP can inform new biological hypotheses and emphasize downstream analyses that reveal sequence or structure motifs important for RNA interactions in cells.

The HIV-1 RNA genome contains complex structures with many structural elements playing regulatory roles during viral replication. A recent study has identified multiple RNA structures with unknown functions that are conserved among HIV-1 and two simian immunodeficiency viruses. To explore the roles of these conserved RNA structures, we introduced synonymous mutations into the HIV-1 genome to disrupt each structure. These mutants exhibited similar particle production, viral infectivity, and replication kinetics relative to the parent NL4-3 virus. However, when replicating in direct competition with the wild-type NL4-3 virus, mutations of RNA structures at inter-protein domain junctions can cause fitness defects. These findings reveal the ability of HIV-1 to tolerate changes in its sequences, even in apparently highly conserved structures, which permits high genetic diversity in HIV-1 population. Our results also suggest that some conserved RNA structures may function to fine-tune viral replication.

RNA viruses encode the information required to usurp cellular metabolism and gene regulation and to enable their own replication in two ways: in the linear sequence of their RNA genomes and in higher-order structures that form when the genomic RNA strand folds back on itself. Application of high-resolution SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) structure probing to viral RNA genomes has identified numerous new regulatory elements, defined new principles by which viral RNAs interact with the cellular host and evade host immune responses, and revealed relationships between virus evolution and RNA structure. This review summarizes our current understanding of genome structure-function interrelationships for RNA viruses, as informed by SHAPE structure probing, and outlines opportunities for future studies.

Understanding RNA-mediated functions requires a detailed characterization of the underlying RNA structure. In many cases, structure probing experiments are performed on RNA that has been "refolded" in some way, which may cause the conformation to differ from that of the native RNA. We used SHAPE-MaP (selective 2'-hydroxyl acylation analyzed by primer extension, read out by mutational profiling), to probe the structure of the Dengue virus (DENV) RNA genome after gentle extraction of the native RNA from intact virions (ex virion) and after heat denaturation and refolding. Comparison of mutiple SHAPE-informed structural features revealed that refolded RNA is more highly structured and samples fewer conformations than the ex virion RNA. Regions with similar structural features are generally those with low SHAPE reactivity and low Shannon entropy (lowSS regions), which correspond to elements with high levels of well-determined structure. This high-structure and low-entropy analysis framework, previously shown to make possible discovery of functional RNA structures, is thus now shown to allow de novo identification of structural elements in a refolded RNA that are likely to recapitulate RNA structures in the ex virion RNA state. Regions with less well-defined structures, which occurred more frequently in the more native-like ex virion RNA and may contain RNA switches, are challenging to recapitulate using refolded RNA.

Summary: Long non-coding RNAs (lncRNAs) are critical regulators of numerous physiological processes and diseases, especially cancers. However, development of lncRNA-based therapies is limited because the mechanisms of many lncRNAs are obscure, and interactions with functional partners, including proteins, remain uncharacterized. The lncRNA SLNCR1 binds to and regulates the androgen receptor (AR) to mediate melanoma invasion and proliferation in an androgen-independent manner. Here, we use biochemical analyses coupled with selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) RNA structure probing to show that the N-terminal domain of AR binds a pyrimidine-rich motif in an unstructured region of SLNCR1. This motif is predictive of AR binding, as we identify an AR-binding motif in lncRNA HOXA11-AS-203. Oligonucleotides that bind either the AR N-terminal domain or the AR RNA motif block the SLNCR1-AR interaction and reduce SLNCR1-mediated melanoma invasion. Delivery of oligos that block SLNCR1-AR interaction thus represent a plausible therapeutic strategy. : Androgen receptor (AR)-RNA complexes have been implicated in cancer, including melanoma. Schmidt et al. demonstrate that AR binds a single strand sequence in the long non-coding RNA (lncRNA) SLNCR. Point mutations or oligonucleotides that abrogate AR binding to SLNCR block melanoma invasion, suggesting that targeting lncRNA-protein complexes holds therapeutic promise. Keywords: long non-coding RNA, hormone receptor, androgen receptor, HOXA11AS, LINC00673, antisense oligos, melanoma, invasion, RNA scaffold, SHAPE

[Image: see text]

RNA molecules are essential for cellular information transfer and gene regulation, and RNAs have been implicated in many human diseases. Messenger and non-coding RNAs contain highly structured elements, and evidence suggests that many of these structures are important for function. Targeting these RNAs with small molecules offers opportunities to therapeutically modulate numerous cellular processes, including those linked to 'undruggable' protein targets. Despite this promise, there is currently only a single class of human-designed small molecules that target RNA used clinically - the linezolid antibiotics. However, a growing number of small-molecule RNA ligands are being identified, leading to burgeoning interest in the field. Here, we discuss principles for discovering small-molecule drugs that target RNA and argue that the overarching challenge is to identify appropriate target structures - namely, in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets. If focus is placed on such druggable binding sites in RNA, extensive knowledge of the typical physicochemical properties of drug-like small molecules could then enable small-molecule drug discovery for RNA targets to become (only) roughly as difficult as for protein targets.

Messenger RNAs (mRNAs) can fold into complex structures that regulate gene expression. Resolving such structures de novo has remained challenging and has limited understanding of the prevalence and functions of mRNA structure. We use SHAPE-MaP experiments in living E. coli cells to derive quantitative, nucleotide-resolution structure models for 194 endogenous transcripts encompassing approximately 400 genes. Individual mRNAs have exceptionally diverse architectures, and most contain well-defined structures. Active translation destabilizes mRNA structure in cells. Nevertheless, mRNA structure remains similar between in-cell and cell-free environments, indicating broad potential for structure-mediated gene regulation. We find that translation efficiency of endogenous genes is regulated by unfolding kinetics of structures overlapping the ribosome binding site. We discover conserved structured elements in 35% of untranslated regions, several of which we validate as novel protein binding motifs. RNA structure regulates every gene studied here in a meaningful way, implying that most functional structures remain to be discovered.

RNA and membraneless organelles Membraneless compartments can form in cells through liquidliquid phase separation (see the Perspective by Polymenidou). But what prevents these cellular condensates from randomly fusing together? Using the RNA-binding protein (RBP) Whi3, Langdon et al. demonstrated that the secondary structure of different RNA components determines the distinct biophysical and biological properties of the two types of condensates that Whi3 forms. Several RBPs, such as FUS and TDP43, contain prion-like domains and are linked to neurodegenerative diseases. These RBPs are usually soluble in the nucleus but can form pathological aggregates in the cytoplasm. Maharana et al. showed that local RNA concentrations determine distinct phase separation behaviors in different subcellular locations. The higher RNA concentrations in the nucleus act as a buffer to prevent phase separation of RBPs; when mislocalized to the cytoplasm, lower RNA concentrations trigger aggregation. Science , this issue p. 922 , p. 918 ; see also p. 859

There are large differences between the cellular environment and the conditions widely used to study RNA in vitro. SHAPE RNA structure probing in Escherichia coli cells has shown that the cellular environment stabilizes both long-range and local tertiary interactions in the adenine riboswitch aptamer domain. Synthetic crowding agents are widely used to understand the forces that stabilize RNA structure and in efforts to recapitulate the cellular environment under simplified experimental conditions. Here, we studied the structure and ligand binding ability of the adenine riboswitch in the presence of the macromolecular crowding agent, polyethylene glycol (PEG). Ethylene glycol and low-molecular mass PEGs destabilized RNA structure and caused the riboswitch to sample secondary structures different from those observed in simple buffered solutions or in cells. In the presence of larger PEGs, longer-range loop-loop interactions were more similar to those in cells than in buffer alone, consistent with prior work showing that larger PEGs stabilize compact RNA states. Ligand affinity was weakened by low-molecular mass PEGs but increased with high-molecular mass PEGs, indicating that PEG cosolvents exert complex chemical and steric effects on RNA structure. Regardless of polymer size, however, nucleotide-resolution structural characteristics observed in cells were not recapitulated in PEG solutions. Our results reveal that the cellular environment is difficult to recapitulate in vitro; mimicking the cellular state will likely require a combination of crowding agents and other chemical species.

SHAPE chemistries exploit small electrophilic reagents that react with the 2′-hydroxyl group to interrogate RNA structure at single-nucleotide resolution. Mutational profiling (MaP) identifies modified residues based on the ability of reverse transcriptase to misread a SHAPE-modified nucleotide and then counting the resulting mutations by massively parallel sequencing. The SHAPE-MaP approach measures the structure of large and transcriptome-wide systems as accurately as for simple model RNAs. This protocol describes the experimental steps, implemented over three days, required to perform SHAPE probing and construct multiplexed SHAPE-MaP libraries suitable for deep sequencing. These steps include RNA folding and SHAPE structure probing, mutational profiling by reverse transcription, library construction, and sequencing. Automated processing of MaP sequencing data is accomplished using two software packages. ShapeMapper converts raw sequencing files into mutational profiles, creates SHAPE reactivity plots, and provides useful troubleshooting information, often within an hour. SuperFold uses these data to model RNA secondary structures, identify regions with well-defined structures, and visualize probable and alternative helices, often in under a day. We illustrate these algorithms with the E. coli thiamine pyrophosphate riboswitch, E. coli 16S rRNA, and HIV-1 genomic RNAs. SHAPE-MaP can be used to make nucleotide-resolution biophysical measurements of individual RNA motifs, rare components of complex RNA ensembles, and entire transcriptomes. The straightforward MaP strategy greatly expands the number, length, and complexity of analyzable RNA structures.

RNA-Puzzles is a collective experiment in blind 3D RNA structure prediction. We report here a third round of RNA-Puzzles. Five puzzles, 4, 8, 12, 13, 14, all structures of riboswitch aptamers and puzzle 7, a ribozyme structure, are included in this round of the experiment. The riboswitch structures include biological binding sites for small molecules (S-adenosyl methionine, cyclic diadenosine monophosphate, 5-amino 4-imidazole carboxamide riboside 5′-triphosphate, glutamine) and proteins (YbxF), and one set describes large conformational changes between ligand-free and ligand-bound states. The Varkud satellite ribozyme is the most recently solved structure of a known large ribozyme. All puzzles have established biological functions and require structural understanding to appreciate their molecular mechanisms. Through the use of fast-track experimental data, including multidimensional chemical mapping, and accurate prediction of RNA secondary structure, a large portion of the contacts in 3D have been predicted correctly leading to similar topologies for the top ranking predictions. Template-based and homology-derived predictions could predict structures to particularly high accuracies. However, achieving biological insights from de novo prediction of RNA 3D structures still depends on the size and complexity of the RNA. Blind computational predictions of RNA structures already appear to provide useful structural information in many cases. Similar to the previous RNA-Puzzles Round II experiment, the prediction of non-Watson–Crick interactions and the observed high atomic clash scores reveal a notable need for an algorithm of improvement. All prediction models and assessment results are available at http://ahsoka.u-strasbg.fr/rnapuzzles/.

Many biological processes are RNA-mediated, but higher-order structures for most RNAs are unknown, which makes it difficult to understand how RNA structure governs function. Here we describe selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) that makes possible de novo and large-scale identification of RNA functional motifs. Sites of 2'-hydroxyl acylation by SHAPE are encoded as noncomplementary nucleotides during cDNA synthesis, as measured by massively parallel sequencing. SHAPE-MaP-guided modeling identified greater than 90% of accepted base pairs in complex RNAs of known structure, and we used it to define a new model for the HIV-1 RNA genome. The HIV-1 model contains all known structured motifs and previously unknown elements, including experimentally validated pseudoknots. SHAPE-MaP yields accurate and high-resolution secondary-structure models, enables analysis of low-abundance RNAs, disentangles sequence polymorphisms in single experiments and will ultimately democratize RNA-structure analysis.

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

There are large differences between the intracellular environment and the conditions widely used to study RNA structure and function in vitro. To assess the effects of the crowded cellular environment on RNA, we examined the structure and ligand binding function of the adenine riboswitch aptamer domain in healthy, growing Escherichia coli cells at single-nucleotide resolution on the minute time scale using SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension). The ligand-bound aptamer structure is essentially the same in cells and in buffer at 1 mM Mg(2+), the approximate Mg(2+) concentration we measured in cells. In contrast, the in-cell conformation of the ligand-free aptamer is much more similar to the fully folded ligand-bound state. Even adding high Mg(2+) concentrations to the buffer used for in vitro analyses did not yield the conformation observed for the free aptamer in cells. The cellular environment thus stabilizes the aptamer significantly more than does Mg(2+) alone. Our results show that the intracellular environment has a large effect on RNA structure that ultimately favors highly organized conformations.

The length of the CAG repeat region in the huntingtin messenger RNA is predictive of Huntington’s disease. Structural studies of CAG repeat-containing RNAs suggest that these sequences form simple hairpin structures; however, in the context of the full-length huntingtin mRNA, CAG repeats may form complex structures that could be targeted for therapeutic intervention. We examined the structures of transcripts spanning the first exon of the huntingtin mRNA with both healthy and disease-prone repeat lengths. In transcripts with 17 to 70 repeats, the CAG sequences base paired extensively with bases in the 5′ UTR and with a conserved region downstream of the CCG repeat region. In huntingtin transcripts with healthy numbers of repeats, the previously observed CAG hairpin was either absent or short. In contrast, in transcripts with disease-associated numbers of repeats, a CAG hairpin was present and extended from a three-helix junction. Our findings demonstrate the profound importance of sequence context in RNA folding and identify specific structural differences between healthy and disease-inducing huntingtin alleles that may be targets for therapeutic intervention.

RNA chaperones are ubiquitous, heterogeneous proteins essential for RNA structural biogenesis and function. We investigated the mechanism of chaperone-mediated RNA folding by following the time-resolved dimerization of the packaging domain of a retroviral RNA at nucleotide resolution. In the absence of the nucleocapsid (NC) chaperone, dimerization proceeded through multiple, slow-folding intermediates. In the presence of NC, dimerization occurred rapidly through a single structural intermediate. The RNA binding domain of heterogeneous nuclear ribonucleoprotein A1 protein, a structurally unrelated chaperone, also accelerated dimerization. Both chaperones interacted primarily with guanosine residues. Replacing guanosine with more weakly pairing inosine yielded an RNA that folded rapidly without a facilitating chaperone. These results show that RNA chaperones can simplify RNA folding landscapes by weakening intramolecular interactions involving guanosine and explain many RNA chaperone activities.

Structured RNA elements within messenger RNA often direct or modulate the cellular production of active proteins. As reviewed here, RNA structures have been discovered that govern nearly every step in protein production: mRNA production and stability; translation initiation, elongation, and termination; protein folding; and cellular localization. Regulatory RNA elements are common within RNAs from every domain of life. This growing body of RNA-mediated mechanisms continues to reveal new ways in which mRNA structure regulates translation. We integrate examples from several different classes of RNA structure-mediated regulation to present a global perspective that suggests that the secondary and tertiary structure of RNA ultimately constitutes an additional level of the genetic code that both guides and regulates protein biosynthesis.

Single-nucleotide-resolution chemical mapping for structured RNA is being rapidly advanced by new chemistries, faster readouts, and coupling to computational algorithms. Recent tests have shown that selective 2´-hydroxyl acylation by primer extension (SHAPE) can give near-zero error rates (0–2%) in modeling the helices of RNA secondary structure. Here, we benchmark the method on six molecules for which crystallographic data are available: tRNA(phe) and 5S rRNA from E. coli; the P4-P6 domain of the Tetrahymena group I ribozyme; and ligand-bound domains from riboswitches for adenine, cyclic di-GMP, and glycine. SHAPE-directed modeling of these highly structured RNAs gave an overall false negative rate (FNR) of 17% and a false discovery rate (FDR) of 21%, with at least one helix prediction error in five of the six cases. Extensive variations of data processing, normalization, and modeling parameters did not significantly mitigate modeling errors. Only one varation, filtering out data collected with deoxyinosine triphosphate during primer extension, gave a modest improvement (FNR=12% and FDR=14%). The residual structure modeling errors are explained by insufficient information content of these RNAs’ SHAPE data, as evaluated by a nonparametric bootstrapping analysis inspired by approaches in phylogenetic inference. Beyond these benchmark cases, bootstrapping analysis suggests low confidence (

Chemical probing of RNA and DNA structure is a widely used and highly informative approach for examining nucleic acid structure and for evaluating interactions with protein and small-molecule ligands. Use of capillary electrophoresis to analyze chemical probing experiments yields hundreds of nucleotides of information per experiment and can be performed on automated instruments. Extraction of the information from capillary electrophoresis electropherograms is a computationally intensive multistep analytical process, and no current software provides rapid, automated, and accurate data analysis. To overcome this bottleneck, we developed a platform-independent, user-friendly software package, QuShape, that yields quantitatively accurate nucleotide reactivity information with minimal user supervision. QuShape incorporates newly developed algorithms for signal decay correction, alignment of time-varying signals within and across capillaries and relative to the RNA nucleotide sequence, and signal scaling across channels or experiments. An analysis-by-reference option enables multiple, related experiments to be fully analyzed in minutes. We illustrate the usefulness and robustness of QuShape by analysis of RNA SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) experiments.

RNA has the intrinsic property to base pair, forming complex structures fundamental to its diverse functions. Here, we develop PARIS, a method based on reversible psoralen crosslinking for global mapping of RNA duplexes with near base-pair resolution in living cells. PARIS analysis in three human and mouse cell types reveals frequent long-range structures, higher-order architectures, and RNA-RNA interactions in trans across the transcriptome. PARIS determines base-pairing interactions on an individual-molecule level, revealing pervasive alternative conformations. We used PARIS-determined helices to guide phylogenetic analysis of RNA structures and discovered conserved long-range and alternative structures. XIST, a long noncoding RNA (lncRNA) essential for X chromosome inactivation, folds into evolutionarily conserved RNA structural domains that span many kilobases. XIST A-repeat forms complex inter-repeat duplexes that nucleate higher-order assembly of the key epigenetic silencing protein SPEN. PARIS is a generally applicable and versatile method that provides novel insights into the RNA structurome and interactome. VIDEO ABSTRACT.

Many RNA structures are comprised of simple secondary structure elements linked by a few, critical, tertiary interactions. SHAPE chemistry has made interrogation of RNA dynamics at single-nucleotide resolution straightforward. However, de novo identification of nucleotides involved in tertiary interactions remains a challenge. Here we show that nucleotides that form non-canonical or tertiary contacts are detected by comparing information obtained using two SHAPE reagents, N-methylisatoic anhydride (NMIA) and 1-methyl-6-nitroisatoic anhydride (1M6). Nucleotides that react preferentially with NMIA exhibit slow local nucleotide dynamics and preferentially adopt the less common C2′-endo ribose conformation. Experiments and first-principle calculations show 1M6 reacts preferentially with nucleotides in which one face of the nucleobase allows an unhindered stacking interaction with the reagent. Differential SHAPE reactivities were used to detect non-canonical and tertiary interactions in four RNAs with diverse structures and to identify pre-formed non-canonical interactions in partially folded RNAs. Differential SHAPE reactivity analysis will enable experimentally concise, large-scale identification of tertiary structure elements and ligand binding sites in complex RNAs and in diverse biological environments.

Polymorphisms near the interferon lambda 3 (IFNL3) gene strongly predict clearance of hepatitis C virus (HCV) infection. We analyzed a variant (rs4803217 G/T) located within the IFNL3 mRNA 3′ untranslated region (UTR); the G allele (protective allele) is associated with elevated therapeutic HCV clearance. We show that the IFNL3 3′ UTR represses mRNA translation and the rs4803217 allele modulates the extent of translational regulation. We analyzed the structures of IFNL3 variant mRNAs at nucleotide resolution by SHAPE-MaP. The rs4803217 G allele mRNA forms well-defined 3′ UTR structure while the T allele mRNA is more dynamic. The observed differences between alleles are among the largest possible RNA structural alterations that can be induced by a single nucleotide change and transform the UTR from a single well-defined conformation to one with multiple dynamic interconverting structures. These data illustrate that non-coding genetic variants can have significant functional effects by impacting RNA structure.

RNA is the central conduit for gene expression. This role depends on an ability to encode information at two levels: in its linear sequence and in the complex structures RNA can form by folding back on itself. Understanding the global structure-function interrelationships mediated by RNA remains a great challenge in molecular and structural biology. In this Account, we discuss evolving work in our laboratory focused on creating facile, generic, quantitative, accurate, and highly informative approaches for understanding RNA structure in biologically important environments. The core innovation derives from our discovery that the nucleophilic reactivity of the ribose 2'-hydroxyl in RNA is gated by local nucleotide flexibility. The 2'-hydroxyl is reactive at conformationally flexible positions but is unreactive at nucleotides constrained by base pairing. Sites of modification in RNA can be detected efficiently either using primer extension or by protection from exoribonucleolytic degradation. This technology is now called SHAPE, for selective 2'-hydroxyl acylation analyzed by primer extension (or protection from exoribonuclease). SHAPE reactivities are largely independent of nucleotide identity but correlate closely with model-free measurements of molecular order. The simple SHAPE reaction is thus a robust, nucleotide-resolution, biophysical measurement of RNA structure. SHAPE can be used to provide an experimental correction to RNA folding algorithms and, in favorable cases, yield kilobase-scale secondary structure predictions with high accuracies. SHAPE chemistry is based on very simple reactive carbonyl centers that can be varied to yield slow- and fast-reacting reagents. Differential SHAPE reactivities can be used to detect specific RNA positions with slow local nucleotide dynamics. These positions, which are often in the C2'-endo conformation, have the potential to function as molecular timers that regulate RNA folding and function. In addition, fast-reacting SHAPE reagents can be used to visualize RNA structural biogenesis and RNA-protein assembly reactions in one second snapshots in very straightforward experiments. The application of SHAPE to challenging problems in biology has revealed surprises in well-studied systems. New regions have been identified that are likely to have critical functional roles on the basis of their high levels of RNA structure. For example, SHAPE analysis of large RNAs, such as authentic viral RNA genomes, suggests that RNA structure organizes regulatory motifs and regulates splicing, protein folding, genome recombination, and ribonucleoprotein assembly. SHAPE has also revealed limitations to the hierarchical model for RNA folding. Continued development and application of SHAPE technologies will advance our understanding of the many ways in which the genetic code is expressed through the underlying structure of RNA.

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) is a powerful approach for characterizing RNA structure and dynamics at single-nucleotide resolution. However, SHAPE technology is limited, sometimes severely, because primer extension detection obscures structural information for approximately 15 nts at the 5' end and 40-60 nts at the 3' end of the RNA. Moreover, detection by primer extension is more complex than the actual structure-selective chemical interrogation step. Here we quantify covalent adducts in RNA directly by adduct-inhibited exoribonuclease degradation. RNA 2'-O-adducts block processivity of a 3'-->5' exoribonuclease, RNase R, to produce fragments that terminate three nucleotides 3' of the modification site. We analyzed the structure of the native thiamine pyrophosphate (TPP) riboswitch aptamer domain and identified large changes in local nucleotide dynamics and global RNA structure upon ligand binding. In addition to numerous changes that can be attributed to ligand recognition, we identify a single nucleotide bulge register shift, distant from the binding site, that stabilizes the ligand-bound structure. Selective 2'-hydroxyl acylation analyzed by protection from exoribonuclease (RNase-detected SHAPE) should prove broadly useful for facile structural analysis of small noncoding RNAs and for RNAs that have functionally critical structures at their 5' and 3' ends.

Proteins play diverse and critical roles in cellular ribonucleoproteins (RNPs) including promoting formation of and stabilizing active RNA conformations. Yet, the conformational changes required to convert large RNAs into an active RNPs have proven difficult to characterize fully. Here we use high-resolution approaches to monitor both local nucleotide flexibility and solvent accessibility for nearly all nucleotides in the bI3 group I intron RNP in four assembly states: the free RNA, maturase-bound RNA, Mrs1-bound RNA, and the complete six-component holocomplex. The free RNA is misfolded relative to the secondary structure required for splicing. The maturase and Mrs1 proteins each stabilized long-range tertiary interactions but neither protein alone induced folding into the functional secondary structure. In contrast, simultaneous binding by both proteins results in large secondary structure rearrangements in the RNA and yielded the catalytically active group I intron structure. Secondary and tertiary folding of the RNA component of the bI3 RNP are thus not independent: RNA folding is strongly non-hierarchical. These results emphasize that protein-mediated stabilization of RNA tertiary interactions functions to pull the secondary structure into an energetically disfavored, but functional, conformation and emphasize a new role for facilitator proteins in RNP assembly.

Tertiary structure prediction is important for understanding structure–function relationships for RNAs whose structures are unknown and for characterizing RNA states recalcitrant to direct analysis. However, it is unknown what root-mean-square deviation (RMSD) corresponds to a statistically significant RNA tertiary structure prediction. We use discrete molecular dynamics to generate RNA-like folds for structures up to 161 nucleotides (nt) that have complex tertiary interactions and then determine the RMSD distribution between these decoys. These distributions are Gaussian-like. The mean RMSD increases with RNA length and is smaller if secondary structure constraints are imposed while generating decoys. The compactness of RNA molecules with true tertiary folds is intermediate between closely packed spheres and a freely jointed chain. We use this scaling relationship to define an expression relating RMSD with the confidence that a structure prediction is better than that expected by chance. This is the prediction significance, and corresponds to a P-value. For a 100-nt RNA, the RMSD of predicted structures should be within 25 Å of the accepted structure to reach the P ≤ 0.01 level if the secondary structure is predicted de novo and within 14 Å if secondary structure information is used as a constraint. This significance approach should be useful for evaluating diverse RNA structure prediction and molecular modeling algorithms.

RNA selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry exploits the discovery that conformationally dynamic nucleotides preferentially adopt configurations that facilitate reaction between the 2'-OH group and a hydroxyl-selective electrophile, such as benzoyl cyanide (BzCN), to form a 2'-O-adduct. BzCN is ideally suited for quantitative, time-resolved analysis of RNA folding and ribonucleoprotein (RNP) assembly mechanisms because this reagent both reacts with flexible RNA nucleotides and also undergoes auto-inactivating hydrolysis with a half-life of 0.25 s at 37 degrees C. RNA folding is initiated by addition of Mg(2+) or protein, or other change in solution conditions, and nucleotide resolution structural images are obtained by adding aliquots of the evolving reaction to BzCN and then 'waiting' for 1 second. Sites of the 2'-O-adduct formation are subsequently scored as stops to primer extension using reverse transcriptase. This time-resolved SHAPE protocol makes it possible to obtain 1-second structural snapshots in time-resolved kinetic studies for RNAs of arbitrary length and complexity in a straightforward and concise experiment.

Hydroxyl-selective electrophiles, including N-methylisatoic anhydride (NMIA) and 1-methyl-7-nitroisatoic anhydride (1M7), are broadly useful for RNA structure analysis because they react preferentially with the ribose 2′-OH group at conformationally unconstrained or flexible nucleotides. Each nucleotide in an RNA has the potential to form an adduct with these reagents to yield a comprehensive, nucleotide-resolution, view of RNA structure. However, it is possible that factors other than local structure modulate reactivity. To evaluate the influence of base identity on the intrinsic reactivity of each nucleotide, we analyze NMIA and 1M7 reactivity using four distinct RNAs, under both native and denaturing conditions. We show that guanosine and adenosine residues have identical intrinsic 2′-hydroxyl reactivities at pH 8.0 and are 1.4 and 1.7 times more reactive than uridine and cytidine, respectively. These subtle, but statistically significant, differences do not impact the ability of selective 2′-hydroxyl acylation analyzed by primer extension-based (SHAPE) methods to establish an RNA secondary structure or monitor RNA folding in solution because base-specific influences are much smaller than the reactivity differences between paired and unpaired nucleotides.

Most functional RNAs require proteins to facilitate formation of their active structures. In the case of the yeast bI3 group I intron, splicing requires binding by two proteins, the intron-encoded bI3 maturase and the nuclear encoded Mrs1. Here, we use selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry coupled with analysis of point mutants to map long-range interactions in this RNA. This analysis reveals two critical features of the free RNA state. First, the catalytic intron is separated from the flanking exons via a stable anchoring helix. This anchoring helix creates an autonomous structural domain for the intron and functions to prevent misfolding with the flanking exons. Second, the thermodynamically most stable structure for the free RNA is not consistent with the catalytically active conformation as phylogenetically conserved elements form stable, non-native structures. These results highlight a fragile bI3 RNA for which binding by protein cofactors functions to promote extensive secondary structure rearrangements that are an obligatory prerequisite for forming the catalytically active tertiary structure.

The mammalian mitochondrial genome encodes 13 proteins, which are synthesized at the direction of nine monocistronic and two dicistronic mRNAs. These mRNAs lack both 5′ and 3′ untranslated regions. The mechanism by which the specialized mitochondrial translational apparatus locates start codons and initiates translation of these leaderless mRNAs is currently unknown. To better understand this mechanism, the secondary structures near the start codons of all 13 open reading frames have been analyzed using RNA SHAPE chemistry. The extent of structure in these mRNAs as assessed experimentally is distinctly lower than would be predicted by current algorithms based on free energy minimization alone. We find that the 5′ ends of all mitochondrial mRNAs are highly unstructured. The first 35 nucleotides for all mitochondrial mRNAs form structures with free energies less favorable than −3 kcal/mol, equal to or less than a single typical base pair. The start codons, which lie at the very 5′ ends of these mRNAs, are accessible within single stranded motifs in all cases, making them potentially poised for ribosome binding. These data are consistent with a model in which the specialized mitochondrial ribosome preferentially allows passage of unstructured 5′ sequences into the mRNA entrance site to participate in translation initiation.

The emergence of catalysis in early genetic polymers such as RNA is considered a key transition in the origin of life, pre-dating the appearance of protein enzymes. DNA also demonstrates the capacity to fold into three-dimensional structures and form catalysts in vitro. However, to what degree these natural biopolymers comprise functionally privileged chemical scaffolds for folding or the evolution of catalysis is not known. The ability of synthetic genetic polymers (XNAs) with alternative backbone chemistries not found in nature to fold into defined structures and bind ligands raises the possibility that these too might be capable of forming catalysts (XNAzymes). Here we report the discovery of such XNAzymes, elaborated in four different chemistries (arabino nucleic acids, ANA; 2'-fluoroarabino nucleic acids, FANA; hexitol nucleic acids, HNA; and cyclohexene nucleic acids, CeNA) directly from random XNA oligomer pools, exhibiting in trans RNA endonuclease and ligase activities. We also describe an XNA-XNA ligase metalloenzyme in the FANA framework, establishing catalysis in an entirely synthetic system and enabling the synthesis of FANA oligomers and an active RNA endonuclease FANAzyme from its constituent parts. These results extend catalysis beyond biopolymers and establish technologies for the discovery of catalysts in a wide range of polymer scaffolds not found in nature. Evolution of catalysis independent of any natural polymer has implications for the definition of chemical boundary conditions for the emergence of life on Earth and elsewhere in the Universe. ispartof: Nature vol:518 issue:7539 pages:427-430 ispartof: location:England status: published

SHAPE-MaP is unique among RNA structure probing strategies in that it both measures flexibility at single-nucleotide resolution and quantifies the uncertainties in these measurements. We report a straightforward analytical framework that incorporates these uncertainties to enable detection of RNA structural differences between any two states, and we use it here to detect RNA-protein interactions in healthy mouse trophoblast stem cells. We validate this approach by analysis of three model cytoplasmic and nuclear ribonucleoprotein complexes, in 2-minute in-cell probing experiments. In contrast, data produced by alternative in-cell SHAPE probing methods correlate poorly (r = 0.2) with those generated by SHAPE-MaP and do not yield accurate signals for RNA-protein interactions. We then examine RNA-protein and RNA-substrate interactions in the RNase MRP complex and, by comparing in-cell interaction sites with disease-associated mutations, characterize these non-coding mutations in terms of molecular phenotype. Together, these results reveal that SHAPE-MaP can define true interaction sites and infer RNA functions under native cellular conditions with limited pre-existing knowledge of the proteins or RNAs involved.

Profound insights regarding nucleic acid structure and function can be gleaned from very simple, direct, and chemistry-based strategies. Our approach strives to incorporate the elegant physical insights that Don Crothers instilled in those who trained in his laboratory. Don emphasized the advantages of focusing on direct and concise experiments even when the final objective was to understand something complex-potentially including the large-scale architectures of the genomes of RNA viruses and the transcriptomes of cells. Here, the author reviews the intellectual path, and a few detours, that led to the development of the SHAPE-MaP and RING-MaP technologies for interrogating RNA structure and function at large scales. The author also argues that greater attention to creating direct, less inferential experiments will convert "omics" investigations into lasting and definitive contributions to our understanding of biological function.

In cells, RNAs likely adopt numerous intermediate conformations prior to formation of functional RNA-protein complexes. We used single-nucleotide resolution selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) to probe the structure of Escherichia coli 16S rRNA in healthy growing bacteria. SHAPE-directed modeling indicated that the predominant steady-state RNA conformational ensemble in dividing cells had a base-paired structure different from that expected on the basis of comparative sequence analysis and high-resolution studies of the 30S ribosomal subunit. We identified the major cause of these differences by stopping ongoing in-cell transcription (in essence, an in-cell RNA structure pulse-chase experiment) which caused the RNA to chase into a structure that closely resembled the expected one. Most helices that formed alternate RNA conformations under growth conditions interact directly with tertiary-binding ribosomal proteins and form a C-shape that surrounds the mRNA channel and decoding site. These in-cell experiments lead to a model in which ribosome assembly factors function as molecular struts to preorganize this intermediate and emphasize that the final stages of ribonucleoprotein assembly involve extensive protein-facilitated RNA conformational changes.

We introduce a melded chemical and computational approach for probing and modeling higher-order intramolecular tertiary interactions in RNA. 2'-Hydroxyl molecular interference (HMX) identifies nucleotides in highly packed regions of an RNA by exploiting the ability of bulky adducts at the 2'-hydroxyl position to disrupt overall RNA structure. HMX was found to be exceptionally selective for quantitative detection of higher-order and tertiary interactions. When incorporated as experimental constraints in discrete molecular dynamics (DMD) simulations, HMX information yielded accurate three-dimensional models, emphasizing the power of molecular interference to guide RNA tertiary structure analysis and fold refinement. In the case of a large, multi-domain RNA, the Tetrahymena group I intron, HMX identified multiple distinct sets of tertiary structure interaction groups in a single, concise experiment.

We describe structural analysis of RNAs by SHAPE chemical probing. RNAs are treated with 1-methyl-7-nitroisatoic anhydride (1M7), a reagent that detects local nucleotides flexibility, and N-methylisatoic anhydride (NMIA) and 1-methyl-6-nitroisatoic anhydride (1M6), reagents which together detect higher-order and non-canonical interactions. Chemical adducts are detected as stops during reverse transcriptase-mediated primer extension. Probing information can be used to infer conformational changes and ligand binding, and to develop highly accurate models of RNA secondary structures.