C2'-endo nucleotides as molecular timers suggested by the folding of an RNA domain

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures

RNA structures play a fundamental role in nearly every aspect of cellular physiology and pathology. Gaining insights into the functions of RNA molecules requires accurate predictions of RNA secondary structures. However, the existing thermodynamic folding models remain less accurate than desired, even when chemical probing data, such as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivities, are used as restraints. Unlike most SHAPE-directed algorithms that only consider SHAPE restraints for base pairing, we extract two-dimensional structural features encoded in SHAPE data and establish robust relationships between characteristic SHAPE patterns and loop motifs of various types (hairpin, internal, and bulge) and lengths (2-11 nucleotides). Such characteristic SHAPE patterns are closely related to the sugar pucker conformations of loop residues. Based on these patterns, we propose a computational method, SHAPELoop, which refines the predicted results of the existing methods, thereby further improving their prediction accuracy. In addition, SHAPELoop can provide information about local or global structural rearrangements (including pseudoknots) and help researchers to easily test their hypothesized secondary structures.

SHAPE Directed Discovery of New Functions in Large RNAs

RNA lies upstream of nearly all biology and functions as the central conduit of information exchange in all cells. RNA molecules encode information both in their primary sequences and in complex structures that form when an RNA folds back on itself. From the time of discovery of mRNA in the late 1950s until quite recently, we had only a rudimentary understanding of RNA structure across vast regions of most messenger and noncoding RNAs. This deficit is now rapidly being addressed, especially by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry, mutational profiling (MaP), and closely related platform technologies that, collectively, create chemical microscopes for RNA. These technologies make it possible to interrogate RNA structure, quantitatively, at nucleotide resolution, and at large scales, for entire mRNAs, noncoding RNAs, and viral RNA genomes. By applying comprehensive structure probing to diverse problems, we and others are showing that control of biological function mediated by RNA structure is ubiquitous across prokaryotic and eukaryotic organisms.Work over the past decade using SHAPE-based analyses has clarified key principles. First, the method of RNA structure probing matters. SHAPE-MaP, with its direct and one-step readout that probes nearly every nucleotide by reaction at the 2'-hydroxyl, gives a more detailed and accurate readout than alternatives. Second, comprehensive chemical probing is essential. Focusing on fragments of large RNAs or using meta-gene or statistical analyses to compensate for sparse data sets misses critical features and often yields structure models with poor predictive power. Finally, every RNA has its own internal structural personality. There are myriad ways in which RNA structure modulates sequence accessibility, protein binding, translation, splice-site choice, phase separation, and other fundamental biological processes. In essentially every instance where we have applied rigorous and quantitative SHAPE technologies to study RNA structure-function interrelationships, new insights regarding biological regulatory mechanisms have emerged. RNA elements with more complex higher-order structures appear more likely to contain high-information-content clefts and pockets that bind small molecules, broadly informing a vigorous field of RNA-targeted drug discovery.The broad implications of this collective work are twofold. First, it is long past time to abandon depiction of large RNAs as simple noodle-like or gently flowing molecules. Instead, we need to emphasize that nearly all RNAs are punctuated with distinctive internal structures, a subset of which modulate function in profound ways. Second, structure probing should be an integral component of any effort that seeks to understand the functional nexuses and biological roles of large RNAs.

Time-Resolved, Single-Molecule, Correlated Chemical Probing of RNA

Capturing the folding dynamics of large, functionally important RNAs has relied primarily on global measurements of structure or on per-nucleotide chemical probing. These approaches infer, but do not directly measure, through-space structural interactions. Here we introduce trimethyloxonium (TMO) as a chemical probe for RNA. TMO alkylates RNA at high levels in seconds, and thereby enables time-resolved, single-molecule, through-space probing of RNA folding using the RING-MaP correlated chemical probing framework. Time-resolved correlations in the RNase P RNA-a functional RNA with a complex structure stabilized by multiple noncanonical interactions-revealed that a long-range tertiary interaction guides native RNA folding for both secondary and tertiary structure. This unanticipated nonhierarchical folding mechanism was directly validated by examining the consequences of concise disruption of the through-space interaction. Single-molecule, time-resolved RNA structure probing with TMO is poised to reveal a wide range of dynamic RNA folding processes and principles of RNA folding.

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.

Deconstructing internal ribosome entry site elements: an update of structural motifs and functional divergences

Beyond the general cap-dependent translation initiation, eukaryotic organisms use alternative mechanisms to initiate protein synthesis. Internal ribosome entry site (IRES) elements are cis-acting RNA regions that promote internal initiation of translation using a cap-independent mechanism. However, their lack of primary sequence and secondary RNA structure conservation, as well as the diversity of host factor requirement to recruit the ribosomal subunits, suggest distinct types of IRES elements. In spite of this heterogeneity, conserved motifs preserve sequences impacting on RNA structure and RNA–protein interactions important for IRES-driven translation. This conservation brings the question of whether IRES elements could consist of basic building blocks, which upon evolutionary selection result in functional elements with different properties. Although RNA-binding proteins (RBPs) perform a crucial role in the assembly of ribonucleoprotein complexes, the versatility and plasticity of RNA molecules, together with their high flexibility and dynamism, determines formation of macromolecular complexes in response to different signals. These properties rely on the presence of short RNA motifs, which operate as modular entities, and suggest that decomposition of IRES elements in short modules could help to understand the different mechanisms driven by these regulatory elements. Here we will review evidence suggesting that model IRES elements consist of the combination of short modules, providing sites of interaction for ribosome subunits, eIFs and RBPs, with implications for definition of criteria to identify novel IRES-like elements genome wide.

Ribosome-dependent conformational flexibility changes and RNA dynamics of IRES domains revealed by differential SHAPE

Internal ribosome entry site (IRES) elements are RNA regions that recruit the translation machinery internally. Here we investigated the conformational changes and RNA dynamics of a picornavirus IRES upon incubation with distinct ribosomal fractions. Differential SHAPE analysis of the free RNA showed that nucleotides reaching the final conformation on long timescales were placed at domains 4 and 5, while candidates for long-range interactions were located in domain 3. Salt-washed ribosomes induced a fast RNA local flexibility modification of domains 2 and 3, while ribosome-associated factors changed domains 4 and 5. Consistent with this, modeling of the three-dimensional RNA structure indicated that incubation of the IRES with native ribosomes induced a local rearrangement of the apical region of domain 3, and a reorientation of domains 4 and 5. Furthermore, specific motifs within domains 2 and 3 showed a decreased flexibility upon incubation with ribosomal subunits in vitro, and presence of the IRES enhanced mRNA association to the ribosomal subunits in whole cell lysates. The finding that RNA modules can provide direct IRES-ribosome interaction suggests that linking these motifs to additional sequences able to recruit trans-acting factors could be useful to design synthetic IRESs with novel activities.

SHAPE analysis of the htrA RNA thermometer from Salmonella enterica

RNA thermometers regulate expression of some genes involved in virulence of pathogenic bacteria such as Yersinia, Neisseria, and Salmonella They often function through temperature-dependent conformational changes that alter accessibility of the ribosome-binding site. The 5'-untranslated region (UTR) of the htrA mRNA from Salmonella enterica contains a very short RNA thermometer. We have systematically characterized the structure and dynamics of this thermometer at single-nucleotide resolution using SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) assays. Our results confirm that the htrA thermometer adopts the predicted hairpin conformation at low temperatures, with conformational change occurring over a physiological temperature regime. Detailed SHAPE melting curves for individual nucleotides suggest that the thermometer unfolds in a cooperative fashion, with nucleotides from both upper and lower portions of the stem gaining flexibility at a common transition temperature. Intriguingly, analysis of an extended htrA 5' UTR sequence revealed not only the presence of the RNA thermometer, but also an additional, stable upstream structure. We generated and analyzed point mutants of the htrA thermometer, revealing elements that modulate its stability, allowing the hairpin to melt under the slightly elevated temperatures experienced during the infection of a warm-blooded host. This work sheds light on structure-function relationships in htrA and related thermometers, and it also illustrates the utility of SHAPE assays for detailed study of RNA thermometer systems.

Genetic variation and RNA structure regulate microRNA biogenesis

MiRNA biogenesis is highly regulated at the post-transcriptional level; however, the role of sequence and secondary RNA structure in this process has not been extensively studied. A single G to A substitution present in the terminal loop of pri-mir-30c-1 in breast and gastric cancer patients had been previously described to result in increased levels of mature miRNA. Here, we report that this genetic variant directly affects Drosha-mediated processing of pri-mir-30c-1 in vitro and in cultured cells. Structural analysis of this variant revealed an altered RNA structure that facilitates the interaction with SRSF3, an SR protein family member that promotes pri-miRNA processing. Our results are compatible with a model whereby a genetic variant in pri-mir-30c-1 leads to a secondary RNA structure rearrangement that facilitates binding of SRSF3 resulting in increased levels of miR-30c. These data highlight that primary sequence determinants and RNA structure are key regulators of miRNA biogenesis.

'Z-DNA like' fragments in RNA: a recurring structural motif with implications for folding, RNA/protein recognition and immune response

Since the work of Alexander Rich, who solved the first Z-DNA crystal structure, we have known that d(CpG) steps can adopt a particular structure that leads to forming left-handed helices. However, it is still largely unrecognized that other sequences can adopt ‘left-handed’ conformations in DNA and RNA, in double as well as single stranded contexts. These ‘Z-like’ steps involve the coexistence of several rare structural features: a C2’-endo puckering, a syn nucleotide and a lone pair–π stacking between a ribose O4’ atom and a nucleobase. This particular arrangement induces a conformational stress in the RNA backbone, which limits the occurrence of Z-like steps to ≈0.1% of all dinucleotide steps in the PDB. Here, we report over 600 instances of Z-like steps, which are located within r(UNCG) tetraloops but also in small and large RNAs including riboswitches, ribozymes and ribosomes. Given their complexity, Z-like steps are probably associated with slow folding kinetics and once formed could lock a fold through the formation of unique long-range contacts. Proteins involved in immunologic response also specifically recognize/induce these peculiar folds. Thus, characterizing the conformational features of these motifs could be a key to understanding the immune response at a structural level.

The role of RNA conformation in RNA-protein recognition

Interactions between protein and RNA play a key role in many biological processes in the gene expression pathway. Those interactions are mediated through a variety of RNA-binding protein domains, among them the highly abundant RNA recognition motif (RRM). Here we studied protein-RNA complexes from different RNA binding domain families solved by NMR and x-ray crystallography. Characterizing the structural properties of the RNA at the binding interfaces revealed an unexpected number of nucleotides with unusual RNA conformations, specifically found in RNA-RRM complexes. Moreover, we observed that the RNA nucleotides that are directly involved in interactions with the RRM domains, via hydrogen bonds and hydrophobic contacts, are significantly enriched with unique RNA conformations. Further examination of the sequences binding the RRM domain showed a preference for G nucleotides in syn conformation to precede or to follow U nucleotides in the anti-conformation, and U nucleotides in C2' endo conformation to precede U and G nucleotides possessing the more common C3' endo conformation. These findings imply a possible mode of RNA recognition by the RRM domains which enables the recognition of a wide variety of different RNA sequences and shapes. Overall, this study suggests an additional way by which the RRM domain recognizes its RNA target, involving a conformational readout.

Quantitative and predictive model of kinetic regulation by E. coli TPP riboswitches

Riboswitches are non-coding elements upstream or downstream of mRNAs that, upon binding of a specific ligand, regulate transcription and/or translation initiation in bacteria, or alternative splicing in plants and fungi. We have studied thiamine pyrophosphate (TPP) riboswitches regulating translation of thiM operon and transcription and translation of thiC operon in E. coli, and that of THIC in the plant A. thaliana. For all, we ascertained an induced-fit mechanism involving initial binding of the TPP followed by a conformational change leading to a higher-affinity complex. The experimental values obtained for all kinetic and thermodynamic parameters of TPP binding imply that the regulation by A. thaliana riboswitch is governed by mass-action law, whereas it is of kinetic nature for the two bacterial riboswitches. Kinetic regulation requires that the RNA polymerase pauses after synthesis of each riboswitch aptamer to leave time for TPP binding, but only when its concentration is sufficient. A quantitative model of regulation highlighted how the pausing time has to be linked to the kinetic rates of initial TPP binding to obtain an ON/OFF switch in the correct concentration range of TPP. We verified the existence of these pauses and the model prediction on their duration. Our analysis also led to quantitative estimates of the respective efficiency of kinetic and thermodynamic regulations, which shows that kinetically regulated riboswitches react more sharply to concentration variation of their ligand than thermodynamically regulated riboswitches. This rationalizes the interest of kinetic regulation and confirms empirical observations that were obtained by numerical simulations.

Fingerprinting the junctions of RNA structure by an open-paddlewheel diruthenium compound

RNA function is determined by its structural organization. The RNA structure consists of the combination of distinct secondary structure motifs connected by junctions that play an essential role in RNA folding. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) probing is an established methodology to analyze the secondary structure of long RNA molecules in solution, which provides accurate data about unpaired nucleotides. However, the residues located at the junctions of RNA structures usually remain undetected. Here we report an RNA probing method based on the use of a novel open-paddlewheel diruthenium (OPW-Ru) compound [Ru2Cl2(µ-DPhF)3(DMSO)] (DPhF = N,N'-diphenylformamidinate). This compound has four potential coordination sites in a singular disposition to establish covalent bonds with substrates. As a proof of concept, we have analyzed the reactivity of OPW-Ru toward RNA using two viral internal ribosome entry site (IRES) elements whose function depends on the structural organization of the molecule. Our study suggests that the compound OPW-Ru preferentially attacks at positions located one or two nucleotides away from junctions or bulges of the RNA structure. The OPW-Ru fingerprinting data differ from that obtained by other chemical reagents and provides new information about RNA structure features.

RNA approaches the B-form in stacked single strand dinucleotide contexts

Duplex RNA adopts an A-form structure, while duplex DNA interconverts between the A- and B-forms depending on the environment. The C2'-endo sugar pucker seen in B-form DNA can occur infrequently in ribose sugars as well, but RNA is not understood to assume B-form conformations. Through analysis of over 45,000 stacked single strand dinucleotide (SSD) crystal structure conformations, this study demonstrates that RNA is capable of adopting a wide conformational range between the canonical A- and B-forms at the localized SSD level, including many B-form-like conformations. It does so through C2'-endo ribose conformations in one or both nucleotides, and B-form-like neighboring base stacking patterns. As chemical reactions on nucleic acids involve localized changes in chemical bonds, the understanding of how enzymes distinguish between DNA and RNA nucleotides is altered by the energetic accessibility of these rare B-form-like RNA SSD conformations. The existence of these conformations also has direct implications in parametrization of molecular mechanics energy functions used extensively to model nucleic acid behavior., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 65-82, 2016.

Global shape mimicry of tRNA within a viral internal ribosome entry site mediates translational reading frame selection

The dicistrovirus intergenic region internal ribosome entry site (IRES) adopts a triple-pseudoknotted RNA structure and occupies the core ribosomal E, P, and A sites to directly recruit the ribosome and initiate translation at a non-AUG codon. A subset of dicistrovirus IRESs directs translation in the 0 and +1 frames to produce the viral structural proteins and a +1 overlapping open reading frame called ORFx, respectively. Here we show that specific mutations of two unpaired adenosines located at the core of the three-helical junction of the honey bee dicistrovirus Israeli acute paralysis virus (IAPV) IRES PKI domain can uncouple 0 and +1 frame translation, suggesting that the structure adopts distinct conformations that contribute to 0 or +1 frame translation. Using a reconstituted translation system, we show that ribosomes assembled on mutant IRESs that direct exclusive 0 or +1 frame translation lack reading frame fidelity. Finally, a nuclear magnetic resonance/small-angle X-ray scattering hybrid approach reveals that the PKI domain of the IAPV IRES adopts an RNA structure that resembles a complete tRNA. The tRNA shape-mimicry enables the viral IRES to gain access to the ribosome tRNA-binding sites and form intermolecular contacts with the ribosome that are necessary for initiating IRES translation in a specific reading frame.

Selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) for direct, versatile and accurate RNA structure analysis

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistries exploit small electrophilic reagents that react with 2'-hydroxyl groups to interrogate RNA structure at single-nucleotide resolution. Mutational profiling (MaP) identifies modified residues by using 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 can be done for simple model RNAs. This protocol describes the experimental steps, implemented over 3 d, that are required to perform SHAPE probing and to construct multiplexed SHAPE-MaP libraries suitable for deep 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. SuperFold uses these data to model RNA secondary structures, identify regions with well-defined structures and visualize probable and alternative helices, often in under 1 d. SHAPE-MaP can be used to make nucleotide-resolution biophysical measurements of individual RNA motifs, rare components of complex RNA ensembles and entire transcriptomes.

Review toward all RNA structures, concisely

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.

Magnesium-dependent folding of a picornavirus IRES element modulates RNA conformation and eIF4G interaction

Internal ribosome entry site (IRES) elements are high-order RNA structures that promote internal initiation of translation to allow protein synthesis under situations that compromise the general cap-dependent translation mechanism. Picornavirus IRES elements are highly efficient elements with a modular RNA structure organization. Here we investigated the effect of Mg(2+) concentration on the local flexibility and solvent accessibility of the foot-and-mouth disease virus (FMDV) IRES element measured on the basis of selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivity and hydroxyl radical cleavage. We have found that Mg(2+) concentration affects the organization of discrete IRES regions, mainly the apical region of domain 3, the 10 nt loop of domain 4, and the pyrimidine tract of domain 5. In support of the effect of RNA structure on IRES activity, substitution or deletion mutants of the 10 nt loop of domain 4 impair internal initiation. In addition, divalent cations affect the binding of eIF4G, a eukaryotic initiation factor that is essential for IRES-dependent translation that interacts with domain 4. Binding of eIF4G is favored by the local RNA flexibility adopted at low Mg(2+) concentration, while eIF4B interacts with the IRES independently of the compactness of the RNA structure. Our study shows that the IRES element adopts a near-native structure in the absence of proteins, shedding light on the influence of Mg(2+) ions on the local flexibility and binding of eIF4G in a model IRES element.

RNA secondary structure modeling at consistent high accuracy using differential SHAPE

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

Hierarchy of RNA functional dynamics

RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.

The cellular environment stabilizes adenine riboswitch RNA structure

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.

Conformational readout of RNA by small ligands

RNA molecules have highly versatile structures that can fold into myriad conformations, providing many potential pockets for binding small molecules. The increasing number of available RNA structures, in complex with proteins, small ligands and in free form, enables the design of new therapeutically useful RNA-binding ligands. Here we studied RNA ligand complexes from 10 RNA groups extracted from the protein data bank (PDB), including adaptive and non-adaptive complexes. We analyzed the chemical, physical, structural and conformational properties of binding pockets around the ligand. Comparing the properties of ligand-binding pockets to the properties of computed pockets extracted from all available RNA structures and RNA-protein interfaces, revealed that ligand-binding pockets, mainly the adaptive pockets, are characterized by unique properties, specifically enriched in rare conformations of the nucleobase and the sugar pucker. Further, we demonstrate that nucleotides possessing the rare conformations are preferentially involved in direct interactions with the ligand. Overall, based on our comprehensive analysis of RNA-ligand complexes, we suggest that the unique conformations adopted by RNA nucleotides play an important role in RNA recognition by small ligands. We term the recognition of a binding site by a ligand via the unique RNA conformations "RNA conformational readout." We propose that "conformational readout" is a general way by which RNA binding pockets are recognized and selected from an ensemble of different RNA states.

Principles for understanding the accuracy of SHAPE-directed RNA structure modeling

Accurate RNA structure modeling is an important, incompletely solved, challenge. Single-nucleotide resolution SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) yields an experimental measurement of local nucleotide flexibility that can be incorporated as pseudo-free energy change constraints to direct secondary structure predictions. Prior work from our laboratory has emphasized both the overall accuracy of this approach and the need for nuanced interpretation of modeled structures. Recent studies by Das and colleagues [Kladwang, W., et al. (2011) Biochemistry 50, 8049; Nat. Chem. 3, 954], focused on analyzing six small RNAs, yielded poorer RNA secondary structure predictions than expected on the basis of prior benchmarking efforts. To understand the features that led to these divergent results, we re-examined four RNAs yielding the poorest results in this recent work: tRNA(Phe), the adenine and cyclic-di-GMP riboswitches, and 5S rRNA. Most of the errors reported by Das and colleagues reflected nonstandard experiment and data processing choices, and selective scoring rules. For two RNAs, tRNA(Phe) and the adenine riboswitch, secondary structure predictions are nearly perfect if no experimental information is included but were rendered inaccurate by the SHAPE data of Das and colleagues. When best practices were used, single-sequence SHAPE-directed secondary structure modeling recovered ~93% of individual base pairs and >90% of helices in the four RNAs, essentially indistinguishable from the results of the mutate-and-map approach with the exception of a single helix in the 5S rRNA. The field of experimentally directed RNA secondary structure prediction is entering a phase focused on the most difficult prediction challenges. We outline five constructive principles for guiding this field forward.

Fingerprinting noncanonical and tertiary RNA structures by differential SHAPE reactivity

Many RNA structures are composed 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 noncanonical or tertiary contacts can be 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 usually adopt the less common C2'-endo ribose conformation. Experiments and first-principles calculations show that 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 noncanonical and tertiary interactions in four RNAs with diverse structures and to identify preformed noncanonical 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.

Structure of the RNA claw of the DNA packaging motor of bacteriophage Φ29

Bacteriophage DNA packaging motors translocate their genomic DNA into viral heads, compacting it to near-crystalline density. The Bacillus subtilis phage ϕ29 has a unique ring of RNA (pRNA) that is an essential component of its motor, serving as a scaffold for the packaging ATPase. Previously, deletion of a three-base bulge (18-CCA-20) in the pRNA A-helix was shown to abolish packaging activity. Here, we solved the structure of this crucial bulge by nuclear magnetic resonance (NMR) using a 27mer RNA fragment containing the bulge (27b). The bulge actually involves five nucleotides (17-UCCA-20 and A100), as U17 and A100 are not base paired as predicted. Mutational analysis showed these newly identified bulge residues are important for DNA packaging. The bulge introduces a 33–35° bend in the helical axis, and inter-helical motion around this bend appears to be restricted. A model of the functional 120b pRNA was generated using a 27b NMR structure and the crystal structure of the 66b prohead-binding domain. Fitting this model into a cryo-EM map generated a pentameric pRNA structure; five helices projecting from the pRNA ring resemble an RNA claw. Biochemical analysis suggested that this shape is important for coordinated motor action required for DNA translocation.

Csy4 relies on an unusual catalytic dyad to position and cleave CRISPR RNA

CRISPR-Cas adaptive immune systems protect prokaryotes against foreign genetic elements. crRNAs derived from CRISPR loci base pair with complementary nucleic acids, leading to their destruction. In Pseudomonas aeruginosa, crRNA biogenesis requires the endoribonuclease Csy4, which binds and cleaves the repetitive sequence of the CRISPR transcript. Biochemical assays and three co-crystal structures of wild-type and mutant Csy4/RNA complexes reveal a substrate positioning and cleavage mechanism in which a histidine deprotonates the ribosyl 2'-hydroxyl pinned in place by a serine, leading to nucleophilic attack on the scissile phosphate. The active site catalytic dyad lacks a general acid to protonate the leaving group and positively charged residues to stabilize the transition state, explaining why the observed catalytic rate constant is ∼10(4)-fold slower than that of RNase A. We show that this RNA cleavage step is essential for assembly of the Csy protein-crRNA complex that facilitates target recognition. Considering that Csy4 recognizes a single cellular substrate and sequesters the cleavage product, evolutionary pressure has likely selected for substrate specificity and high-affinity crRNA interactions at the expense of rapid cleavage kinetics.

The mechanisms of RNA SHAPE chemistry

The biological functions of RNA are ultimately governed by the local environment at each nucleotide. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry is a powerful approach for measuring nucleotide structure and dynamics in diverse biological environments. SHAPE reagents acylate the 2'-hydroxyl group at flexible nucleotides because unconstrained nucleotides preferentially sample rare conformations that enhance the nucleophilicity of the 2'-hydroxyl. The critical corollary is that some constrained nucleotides must be poised for efficient reaction at the 2'-hydroxyl group. To identify such nucleotides, we performed SHAPE on intact crystals of the Escherichia coli ribosome, monitored the reactivity of 1490 nucleotides in 16S rRNA, and examined those nucleotides that were hyper-reactive toward SHAPE and had well-defined crystallographic conformations. Analysis of these conformations revealed that 2'-hydroxyl reactivity is broadly facilitated by general base catalysis involving multiple RNA functional groups and by two specific orientations of the bridging 3'-phosphate group. Nucleotide analog studies confirmed the contributions of these mechanisms to SHAPE reactivity. These results provide a strong mechanistic explanation for the relationship between SHAPE reactivity and local RNA dynamics and will facilitate interpretation of SHAPE information in the many technologies that make use of this chemistry.

Exploring RNA structural codes with SHAPE chemistry

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.

Solution structure of RNase P RNA

The ribonucleoprotein enzyme ribonuclease P (RNase P) processes tRNAs by cleavage of precursor-tRNAs. RNase P is a ribozyme: The RNA component catalyzes tRNA maturation in vitro without proteins. Remarkable features of RNase P include multiple turnovers in vivo and ability to process diverse substrates. Structures of the bacterial RNase P, including full-length RNAs and a ternary complex with substrate, have been determined by X-ray crystallography. However, crystal structures of free RNA are significantly different from the ternary complex, and the solution structure of the RNA is unknown. Here, we report solution structures of three phylogenetically distinct bacterial RNase P RNAs from Escherichia coli, Agrobacterium tumefaciens, and Bacillus stearothermophilus, determined using small angle X-ray scattering (SAXS) and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis. A combination of homology modeling, normal mode analysis, and molecular dynamics was used to refine the structural models against the empirical data of these RNAs in solution under the high ionic strength required for catalytic activity.

The shape-shifting quasispecies of RNA: one sequence, many functional folds

E Unus pluribum, or "Of One, Many", may be at the root of decoding the RNA sequence-structure-function relationship. RNAs embody the large majority of genes in higher eukaryotes and fold in a sequence-directed fashion into three-dimensional structures that perform functions conserved across all cellular life forms, ranging from regulating to executing gene expression. While it is the most important determinant of the RNA structure, the nucleotide sequence is generally not sufficient to specify a unique set of secondary and tertiary interactions due to the highly frustrated nature of RNA folding. This frustration results in folding heterogeneity, a common phenomenon wherein a chemically homogeneous population of RNA molecules folds into multiple stable structures. Often, these alternative conformations constitute misfolds, lacking the biological activity of the natively folded RNA. Intriguingly, a number of RNAs have recently been described as capable of adopting multiple distinct conformations that all perform, or contribute to, the same function. Characteristically, these conformations interconvert slowly on the experimental timescale, suggesting that they should be regarded as distinct native states. We discuss how rugged folding free energy landscapes give rise to multiple native states in the Tetrahymena Group I intron ribozyme, hairpin ribozyme, sarcin-ricin loop, ribosome, and an in vitro selected aptamer. We further describe the varying degrees to which folding heterogeneity impacts function in these RNAs, and compare and contrast this impact with that of heterogeneities found in protein folding. Embracing that one sequence can give rise to multiple native folds, we hypothesize that this phenomenon imparts adaptive advantages on any functionally evolving RNA quasispecies.

Versatile phosphoramidation reactions for nucleic acid conjugations with peptides, proteins, chromophores, and biotin derivatives

Chemical conjugations of nucleic acids with macromolecules or small molecules are common approaches to study nucleic acids in chemistry and biology and to exploit nucleic acids for medical applications. The conjugation of nucleic acids such as oligonucleotides with peptides is especially useful to circumvent cell delivery and specificity problems of oligonucleotides as therapeutic agents. However, current approaches are limited and inefficient in their ability to afford peptide-oligonucleotide conjugates (POCs). Here, we report an effective and reproducible approach to prepare POCs and other nucleic acid conjugates based on a newly developed nucleic acid phosphoramidation method. The development of a new nucleic acid phosphoramidation reaction was achieved by our successful synthesis of a novel amine-containing biotin derivative used to systematically optimize the reactions. The improved phosphoramidation reactions dramatically increased yields of nucleic acid-biotin conjugates up to 80% after 3 h reaction. Any nucleic acids with a terminal phosphate group are suitable reactants in phosphoramidation reactions to conjugate with amine-containing molecules such as biotin and fluorescein derivatives, proteins, and, most importantly, peptides to enable the synthesis of POCs for therapeutic applications. Polymerase chain reactions (PCRs) to study incorporation of biotin or fluorescein-tagged DNA primers into the reaction products demonstrated that appropriate controls of nucleic acid phosphoramidation reactions incur minimum adverse effects on inherited base-pairing characteristics of nucleotides in nucleic acids. The phosphoramidation approach preserves the integrity of hybridization specificity in nucleic acids when preparing POCs. By retaining integrity of the nucleic acids, their effectiveness as therapeutic reagents for gene silencing, gene therapy, and RNA interference is ensured. The potential for POC use was demonstrated by two-step phosphoramidation reactions to successfully synthesize nucleic acid-tetraglycine conjugates. In addition, phosphoramidation reactions provided a facile approach to prepare nucleic acid-BSA conjugates with good yields. In summary, the new approach to phosphoramidation reactions offers a universal method to prepare POCs and other nucleic acid conjugates with high yields in aqueous solutions. The methods can be easily adapted to typical chemistry or biology laboratory setups which will expedite the applications of POCs for basic research and medicine.

Taming free energy landscapes with RNA chaperones

Many non-coding RNAs fold into complex three-dimensional structures, yet the self-assembly of RNA structure is hampered by mispairing, weak tertiary interactions, electrostatic barriers, and the frequent requirement that the 5' and 3' ends of the transcript interact. This rugged free energy landscape for RNA folding means that some RNA molecules in a population rapidly form their native structure, while many others become kinetically trapped in misfolded conformations. Transient binding of RNA chaperone proteins destabilize misfolded intermediates and lower the transition states between conformations, producing a smoother landscape that increases the rate of folding and the probability that a molecule will find the native structure. DEAD-box proteins couple the chemical potential of ATP hydrolysis with repetitive cycles of RNA binding and release, expanding the range of conditions under which they can refold RNA structures.

Advances in RNA structure analysis by chemical probing

RNA is arguably the most versatile biological macromolecule because of its ability both to encode and to manipulate genetic information. The diverse roles of RNA depend on its ability to fold back on itself to form biologically functional structures that bind small molecule and large protein ligands, to change conformation, and to affect the cellular regulatory state. These features of RNA biology can be structurally interrogated using chemical mapping experiments. The usefulness and applications of RNA chemical probing technologies have expanded dramatically over the past five years because of several critical advances. These innovations include new sequence-independent RNA chemistries, algorithmic tools for high-throughput analysis of complex data sets composed of thousands of measurements, new approaches for interpreting chemical probing data for both secondary and tertiary structure prediction, facile methods for following time-dependent processes, and the willingness of individual research groups to tackle increasingly bold problems in RNA structural biology.