The paradoxical behavior of a highly structured misfolded intermediate in RNA folding

Tracking Native Tetrahymena Ribozyme Folding with Fluorescence

Folding of the Tetrahymena group I intron ribozyme and other structured RNAs has been measured using a catalytic activity assay to monitor the native state formation by cleavage of a radiolabeled oligonucleotide substrate. While highly effective, the assay has inherent limitations present in any radioactivity- and gel-based assay. Administrative and safety considerations arise from the radioisotope, and data collection is laborious due to the use of polyacrylamide gels. Here we describe a fluorescence-based, solution assay that allows for more efficient data acquisition. The substrate is labeled with 6-carboxyfluorescein (6FAM) fluorophore and black hole quencher (BHQ1) at the 5' and 3' ends, respectively. Substrate cleavage results in release of the quencher, increasing the fluorescence signal by an average of 30-fold. A side-by-side comparison with the radioactivity-based assay shows good agreement in monitoring Tetrahymena ribozyme folding from a misfolded conformation to the native state, albeit with increased uncertainty. The lower precision of the fluorescence assay is compensated for by the relative ease and efficiency of the workflow. In addition, this assay will allow institutions that do not use radioactive materials to monitor native folding of the Tetrahymena ribozyme, and the same strategy should be amenable to native folding of other ribozymes.

Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time

Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process.

Topological crossing in the misfolded Tetrahymena ribozyme resolved by cryo-EM

Taking advantage of single-particle cryogenic electron microscopy (cryo-EM) to analyze highly heterogeneous or flexible samples, we obtained long-awaited three-dimensional (3D) structures of the misfolded Tetrahymena ribozyme. These structures provide clear evidence for a previously proposed topological isomer model, in which the stereochemically impossible crossing of two core RNA strands prevents rapid rearrangement of the misfolded state to the native state. Topological isomers may be widespread in misfolding of complex RNA, and these cryo-EM structures set a foundation for dissecting their detailed kinetic mechanisms and functional consequences in a paradigmatic model system.

The Tetrahymena group I intron has been a key system in the understanding of RNA folding and misfolding. The molecule folds into a long-lived misfolded intermediate (M) in vitro, which has been known to form extensive native-like secondary and tertiary structures but is separated by an unknown kinetic barrier from the native state (N). Here, we used cryogenic electron microscopy (cryo-EM) to resolve misfolded structures of the Tetrahymena L-21 ScaI ribozyme. Maps of three M substates (M1, M2, M3) and one N state were achieved from a single specimen with overall resolutions of 3.5 Å, 3.8 Å, 4.0 Å, and 3.0 Å, respectively. Comparisons of the structures reveal that all the M substates are highly similar to N, except for rotation of a core helix P7 that harbors the ribozyme’s guanosine binding site and the crossing of the strands J7/3 and J8/7 that connect P7 to the other elements in the ribozyme core. This topological difference between the M substates and N state explains the failure of 5′-splice site substrate docking in M, supports a topological isomer model for the slow refolding of M to N due to a trapped strand crossing, and suggests pathways for M-to-N refolding.

Cryo-EM reveals an entangled kinetic trap in the folding of a catalytic RNA

Functional RNAs fold through complex pathways that can contain misfolded "kinetic traps." A complete model of RNA folding requires understanding the formation of these misfolded states, but they are difficult to characterize because of their transient and potentially conformationally dynamic nature. We used cryo-electron microscopy (cryo-EM) to visualize a long-lived misfolded state in the folding pathway of the Tetrahymena thermophila group I intron, a paradigmatic RNA structure-function model system. The structure revealed how this state forms native-like secondary structure and tertiary contacts but contains two incorrectly crossed strands, consistent with a previous model. This incorrect topology mispositions a critical catalytic domain and cannot be resolved locally as extensive refolding is required. This work provides a structural framework for interpreting decades of biochemical and functional studies and demonstrates the power of cryo-EM for the exploration of RNA folding pathways.

Thoughts on how to think (and talk) about RNA structure

Recent events have pushed RNA research into the spotlight. Continued discoveries of RNA with unexpected diverse functions in healthy and diseased cells, such as the role of RNA as both the source and countermeasure to a severe acute respiratory syndrome coronavirus 2 infection, are igniting a new passion for understanding this functionally and structurally versatile molecule. Although RNA structure is key to function, many foundational characteristics of RNA structure are misunderstood, and the default state of RNA is often thought of and depicted as a single floppy strand. The purpose of this perspective is to help adjust mental models, equipping the community to better use the fundamental aspects of RNA structural information in new mechanistic models, enhance experimental design to test these models, and refine data interpretation. We discuss six core observations focused on the inherent nature of RNA structure and how to incorporate these characteristics to better understand RNA structure. We also offer some ideas for future efforts to make validated RNA structural information available and readily used by all researchers.

Measurement of ATP utilization in RNA unwinding and RNA chaperone activities by DEAD-box helicase proteins

RNA helicase proteins perform coupled reactions in which cycles of ATP binding and hydrolysis are used to drive local unwinding of double-stranded RNA (dsRNA). For some helicases in the ubiquitous DEAD-box family, these local unwinding events are integral to folding transitions in structured RNAs, and thus these helicases function as RNA chaperones. An important measure of the efficiency of the helicase-catalyzed reaction is the ATP utilization value, which represents the average number of ATP molecules hydrolyzed during RNA unwinding or a chaperone-assisted RNA structural rearrangement. Here we outline procedures that can be used to measure the ATP utilization value in RNA unwinding or folding transitions. As an example of an RNA folding transition, we focus on the refolding of the Tetrahymena thermophila group I intron ribozyme from a long-lived misfolded structure to its native structure, and we outline strategies for adapting this assay to other RNA folding transitions. For a simple dsRNA unwinding event, the ATP utilization value provides a measure of the coupling between the ATPase and RNA unwinding activities, and for a complex RNA structural transition it can give insight into the scope of the rearrangement and the efficiency with which the helicase uses the energy from ATPase cycles to promote the rearrangement.

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. We isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

How to Kinetically Dissect an RNA Machine

RNA-based machines are ubiquitous in Nature and increasingly important for medicines. They fold into complex, dynamic structures that process information and catalyze reactions, including reactions that generate new RNAs and proteins across biology. What are the experimental strategies and steps that are necessary to understand how these complex machines work? Two 1990 papers from Herschlag and Cech on "Catalysis of RNA Cleavage by the Tetrahymena thermophila Ribozyme" provide a master class in dissecting an RNA machine through kinetics approaches. By showing how to propose a kinetic framework, fill in the numbers, do cross-checks, and make comparisons across mutants and different RNA systems, the papers illustrate how to take a mechanistic approach and distill the results into general insights that are difficult to attain through other means.

Cryo-EM structures of full-length Tetrahymena ribozyme at 3.1 Å resolution

Single-particle cryogenic electron microscopy (cryo-EM) has become a standard technique for determining protein structures at atomic resolution^1-3. However, cryo-EM studies of protein-free RNA are in their early days. The Tetrahymena thermophila group I self-splicing intron was the first ribozyme to be discovered and has been a prominent model system for the study of RNA catalysis and structure-function relationships⁴, but its full structure remains unknown. Here we report cryo-EM structures of the full-length Tetrahymena ribozyme in substrate-free and bound states at a resolution of 3.1 Å. Newly resolved peripheral regions form two coaxially stacked helices; these are interconnected by two kissing loop pseudoknots that wrap around the catalytic core and include two previously unforeseen (to our knowledge) tertiary interactions. The global architecture is nearly identical in both states; only the internal guide sequence and guanosine binding site undergo a large conformational change and a localized shift, respectively, upon binding of RNA substrates. These results provide a long-sought structural view of a paradigmatic RNA enzyme and signal a new era for the cryo-EM-based study of structure-function relationships in ribozymes.

Sisyphus observed: Unraveling the high ATP usage of an RNA chaperone

DEAD-box proteins are nonprocessive RNA helicases that can function as RNA chaperones by coupling ATP binding and hydrolysis to structural reorganization of RNA. Here, Jarmoskaite et al. quantify the ATP utilization of an RNA chaperone during refolding of a misfolded ribozyme substrate. Strikingly, 100 ATP hydrolysis events are needed per successfully refolded ribozyme, suggesting that each round of unfolding requires ten ATP molecules, since 90% of substrate unfolding cycles only lead back to the kinetically favored misfolded state. This near-Sisyphean effort reveals a potentially conserved model for RNA reorganization by RNA chaperones.

ATP utilization by a DEAD-box protein during refolding of a misfolded group I intron ribozyme

DEAD-box helicase proteins perform ATP-dependent rearrangements of structured RNAs throughout RNA biology. Short RNA helices are unwound in a single ATPase cycle, but the ATP requirement for more complex RNA structural rearrangements is unknown. Here we measure the amount of ATP used for native refolding of a misfolded group I intron ribozyme by CYT-19, a Neurospora crassa DEAD-box protein that functions as a general chaperone for mitochondrial group I introns. By comparing the rates of ATP hydrolysis and ribozyme refolding, we find that several hundred ATP molecules are hydrolyzed during refolding of each ribozyme molecule. After subtracting nonproductive ATP hydrolysis that occurs in the absence of ribozyme refolding, we find that approximately 100 ATPs are hydrolyzed per refolded RNA as a consequence of interactions specific to the misfolded ribozyme. This value is insensitive to changes in ATP and CYT-19 concentration and decreases with decreasing ribozyme stability. Because of earlier findings that ∼90% of global ribozyme unfolding cycles lead back to the kinetically preferred misfolded conformation and are not observed, we estimate that each global unfolding cycle consumes ∼10 ATPs. Our results indicate that CYT-19 functions as a general RNA chaperone by using a stochastic, energy-intensive mechanism to promote RNA unfolding and refolding, suggesting an evolutionary convergence with protein chaperones.

Antisense probing of dynamic RNA structures

RNA regulation is influenced by the dynamic changes in conformational accessibility on the transcript. Here we discuss the initial validation of a cell-free antisense probing method for structured RNAs, using the Tetrahymena group I intron as a control target. We observe changes in signal that qualitatively match prior traditional DMS footprinting experiments. Importantly, we have shown that application of this technique can elucidate new RNA information given its sensitivity for detecting rare intermediates that are not as readily observed by single-hit kinetics chemical probing techniques. Observing changes in RNA accessibility has broad applications in determining the effect that regulatory elements have on regional structures. We speculate that this method could be useful in quickly observing those interactions, along with other phenomena that influence RNA accessibility including RNA-RNA interactions and small molecules.

Probing RNA-Protein Interactions and RNA Compaction by Sedimentation Velocity Analytical Ultracentrifugation

Recent advances in multi-wavelength analytical ultracentrifugation (MWL-AUC) combine the power of an exquisitely sensitive hydrodynamic-based separation technique with the added dimension of spectral separation. This added dimension has opened up new doors to much improved characterization of multiple, interacting species in solution. When applied to structural investigations of RNA, MWL-AUC can precisely report on the hydrodynamic radius and the overall shape of an RNA molecule by enabling precise measurements of its sedimentation and diffusion coefficients and identify the stoichiometry of interacting components based on spectral decomposition. Information provided in this chapter will allow an investigator to design experiments for probing ion and/or protein-induced global conformational changes of an RNA molecule and exploit spectral differences between proteins and RNA to characterize their interactions in a physiological solution environment.

Ribosomal Protein L11 Selectively Stabilizes a Tertiary Structure of the GTPase Center rRNA Domain

The GTPase Center (GAC) RNA domain in bacterial 23S rRNA is directly bound by ribosomal protein L11, and this complex is essential to ribosome function. Previous cocrystal structures of the 58-nucleotide GAC RNA bound to L11 revealed the intricate tertiary fold of the RNA domain, with one monovalent and several divalent ions located in specific sites within the structure. Here, we report a new crystal structure of the free GAC that is essentially identical to the L11-bound structure, which retains many common sites of divalent ion occupation. This new structure demonstrates that RNA alone folds into its tertiary structure with bound divalent ions. In solution, we find that this tertiary structure is not static, but rather is best described as an ensemble of states. While L11 protein cannot bind to the GAC until the RNA has adopted its tertiary structure, new experimental data show that L11 binds to Mg²⁺-dependent folded states, which we suggest lie along the folding pathway of the RNA. We propose that L11 stabilizes a specific GAC RNA tertiary state, corresponding to the crystal structure, and that this structure reflects the functionally critical conformation of the rRNA domain in the fully assembled ribosome.

The Story of RNA Folding, as Told in Epochs

The past decades have witnessed tremendous developments in our understanding of RNA biology. At the core of these advances have been studies aimed at discerning RNA structure and at understanding the forces that influence the RNA folding process. It is easy to take the present state of understanding for granted, but there is much to be learned by considering the path to our current understanding, which has been tortuous, with the birth and death of models, the adaptation of experimental tools originally developed for characterization of protein structure and catalysis, and the development of novel tools for probing RNA. In this review we tour the stages of RNA folding studies, considering them as "epochs" that can be generalized across scientific disciplines. These epochs span from the discovery of catalytic RNA, through biophysical insights into the putative primordial RNA World, to characterization of structured RNAs, the building and testing of models, and, finally, to the development of models with the potential to yield generalizable predictive and quantitative models for RNA conformational, thermodynamic, and kinetic behavior. We hope that this accounting will aid others as they navigate the many fascinating questions about RNA and its roles in biology, in the past, present, and future.

Kinetic Basis for DNA Target Specificity of CRISPR-Cas12a

Class 2 CRISPR-Cas nucleases are programmable genome editing tools with promising applications in human health and disease. However, DNA cleavage at off-target sites that resemble the target sequence is a pervasive problem that remains poorly understood mechanistically. Here, we use quantitative kinetics to dissect the reaction steps of DNA targeting by Acidaminococcus sp Cas12a (also known as Cpf1). We show that Cas12a binds DNA tightly in two kinetically separable steps. Protospacer-adjacent motif (PAM) recognition is followed by rate-limiting R-loop propagation, leading to inevitable DNA cleavage of both strands. Despite functionally irreversible binding, Cas12a discriminates strongly against mismatches along most of the DNA target sequence. This result implies substantial reversibility during R-loop formation-a late transition state-and defies common descriptions of a "seed" region. Our results provide a quantitative basis for the DNA cleavage patterns measured in vivo and observations of greater reported target specificity for Cas12a than for the Cas9 nuclease.

Distinct RNA-unwinding mechanisms of DEAD-box and DEAH-box RNA helicase proteins in remodeling structured RNAs and RNPs

Structured RNAs and RNA-protein complexes (RNPs) fold through complex pathways that are replete with misfolded traps, and many RNAs and RNPs undergo extensive conformational changes during their functional cycles. These folding steps and conformational transitions are frequently promoted by RNA chaperone proteins, notably by superfamily 2 (SF2) RNA helicase proteins. The two largest families of SF2 helicases, DEAD-box and DEAH-box proteins, share evolutionarily conserved helicase cores, but unwind RNA helices through distinct mechanisms. Recent studies have advanced our understanding of how their distinct mechanisms enable DEAD-box proteins to disrupt RNA base pairs on the surfaces of structured RNAs and RNPs, while some DEAH-box proteins are adept at disrupting base pairs in the interior of RNPs. Proteins from these families use these mechanisms to chaperone folding and promote rearrangements of structured RNAs and RNPs, including the spliceosome, and may use related mechanisms to maintain cellular messenger RNAs in unfolded or partially unfolded conformations.

Optimization of a novel biophysical model using large scale in vivo antisense hybridization data displays improved prediction capabilities of structurally accessible RNA regions

Current approaches to design efficient antisense RNAs (asRNAs) rely primarily on a thermodynamic understanding of RNA–RNA interactions. However, these approaches depend on structure predictions and have limited accuracy, arguably due to overlooking important cellular environment factors. In this work, we develop a biophysical model to describe asRNA–RNA hybridization that incorporates in vivo factors using large-scale experimental hybridization data for three model RNAs: a group I intron, CsrB and a tRNA. A unique element of our model is the estimation of the availability of the target region to interact with a given asRNA using a differential entropic consideration of suboptimal structures. We showcase the utility of this model by evaluating its prediction capabilities in four additional RNAs: a group II intron, Spinach II, 2-MS2 binding domain and glgC 5΄ UTR. Additionally, we demonstrate the applicability of this approach to other bacterial species by predicting sRNA–mRNA binding regions in two newly discovered, though uncharacterized, regulatory RNAs.

The DEAD-Box Protein CYT-19 Uses Arginine Residues in Its C-Tail To Tether RNA Substrates

DEAD-box proteins are nonprocessive RNA helicases that play diverse roles in cellular processes. The Neurospora crassa DEAD-box protein CYT-19 promotes mitochondrial group I intron splicing and functions as a general RNA chaperone. CYT-19 includes a disordered, arginine-rich "C-tail" that binds RNA, positioning the helicase core to capture and unwind nearby RNA helices. Here we probed the C-tail further by varying the number and positions of arginines within it. We found that removing sets of as few as four of the 11 arginines reduced RNA unwinding activity (k_cat/K_M) to a degree equivalent to that seen upon removal of the C-tail, suggesting that a minimum or "threshold" number of arginines is required. In addition, a mutant with 16 arginines displayed RNA unwinding activity greater than that of wild-type CYT-19. The C-tail modifications impacted unwinding only of RNA helices within constructs that included an adjacent helix or structured RNA element that would allow C-tail binding, indicating that the helicase core remained active in the mutants. In addition, changes in RNA unwinding efficiency of the mutants were mirrored by changes in functional RNA affinity, as determined from the RNA concentration dependence of ATPase activity, suggesting that the C-tail functions primarily to increase RNA affinity. Interestingly, the salt concentration dependence of RNA unwinding activity is unaffected by C-tail composition, suggesting that the C-tail uses primarily hydrogen bonding, not electrostatic interactions, to bind double-stranded RNA. Our results provide insights into how an unstructured C-tail contributes to DEAD-box protein activity and suggest parallels with other families of RNA- and DNA-binding proteins.

Capturing RNA Folding Free Energy with Coarse-Grained Molecular Dynamics Simulations

We introduce a coarse-grained RNA model for molecular dynamics simulations, RACER (RnA CoarsE-gRained). RACER achieves accurate native structure prediction for a number of RNAs (average RMSD of 2.93 Å) and the sequence-specific variation of free energy is in excellent agreement with experimentally measured stabilities (R2 = 0.93). Using RACER, we identified hydrogen-bonding (or base pairing), base stacking, and electrostatic interactions as essential driving forces for RNA folding. Also, we found that separating pairing vs. stacking interactions allowed RACER to distinguish folded vs. unfolded states. In RACER, base pairing and stacking interactions each provide an approximate stability of 3-4 kcal/mol for an A-form helix. RACER was developed based on PDB structural statistics and experimental thermodynamic data. In contrast with previous work, RACER implements a novel effective vdW potential energy function, which led us to re-parameterize hydrogen bond and electrostatic potential energy functions. Further, RACER is validated and optimized using a simulated annealing protocol to generate potential energy vs. RMSD landscapes. Finally, RACER is tested using extensive equilibrium pulling simulations (0.86 ms total) on eleven RNA sequences (hairpins and duplexes).

Nucleobases Undergo Dynamic Rearrangements during RNA Tertiary Folding

The tertiary structure of the GTPase center (GAC) of 23S ribosomal RNA (rRNA) as seen in cocrystals is extremely compact. It is stabilized by long-range hydrogen bonds and nucleobase stacking and by a triloop that forms within its three-way junction. Its folding pathway from secondary structure to tertiary structure has not been previously observed, but it was shown to require Mg²⁺ ions in equilibrium experiments. The fluorescent nucleotide 2-aminopurine was substituted at selected sites within the 60-nt GAC. Fluorescence intensity changes upon addition of MgCl₂ were monitored over a time-course from 1ms to 100s as the RNA folds. The folding pathway is revealed here to be hierarchical through several intermediates. Observation of the nucleobases during folding provides a new perspective on the process and the pathway, revealing the dynamics of nucleobase conformational exchange during the folding transitions.

Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway

The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.

Formation of Tertiary Interactions during rRNA GTPase Center Folding

The 60-nt GTPase center (GAC) of 23S rRNA has a phylogenetically conserved secondary structure with two hairpin loops and a 3-way junction. It folds into an intricate tertiary structure upon addition of Mg(2+) ions, which is stabilized by the L11 protein in cocrystal structures. Here, we monitor the kinetics of its tertiary folding and Mg(2+)-dependent intermediate states by observing selected nucleobases that contribute specific interactions to the GAC tertiary structure in the cocrystals. The fluorescent nucleobase 2-aminopurine replaced three individual adenines, two of which make long-range stacking interactions and one that also forms hydrogen bonds. Each site reveals a unique response to Mg(2+) addition and temperature, reflecting its environmental change from secondary to tertiary structure. Stopped-flow fluorescence experiments revealed that kinetics of tertiary structure formation upon addition of MgCl2 are also site specific, with local conformational changes occurring from 5 ms to 4s and with global folding from 1 to 5s. Site-specific substitution with (15)N-nucleobases allowed observation of stable hydrogen bond formation by NMR experiments. Equilibrium titration experiments indicate that a stable folding intermediate is present at stoichiometric concentrations of Mg(2+) and suggest that there are two initial sites of Mg(2+) ion association.

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

MOTIVATION

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

RESULTS

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

AVAILABILITY AND IMPLEMENTATION

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

CONTACT

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

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exploiting post-transcriptional regulation to probe RNA structures in vivo via fluorescence

While RNA structures have been extensively characterized in vitro, very few techniques exist to probe RNA structures inside cells. Here, we have exploited mechanisms of post-transcriptional regulation to synthesize fluorescence-based probes that assay RNA structures in vivo. Our probing system involves the co-expression of two constructs: (i) a target RNA and (ii) a reporter containing a probe complementary to a region in the target RNA attached to an RBS-sequestering hairpin and fused to a sequence encoding the green fluorescent protein (GFP). When a region of the target RNA is accessible, the area can interact with its complementary probe, resulting in fluorescence. By using this system, we observed varied patterns of structural accessibility along the length of the Tetrahymena group I intron. We performed in vivo DMS footprinting which, along with previous footprinting studies, helped to explain our probing results. Additionally, this novel approach represents a valuable tool to differentiate between RNA variants and to detect structural changes caused by subtle mutations. Our results capture some differences from traditional footprinting assays that could suggest that probing in vivo via oligonucleotide hybridization facilitates the detection of folding intermediates. Importantly, our data indicate that intracellular oligonucleotide probing can be a powerful complement to existing RNA structural probing methods.

DEAD-box helicase proteins disrupt RNA tertiary structure through helix capture

Single-molecule fluorescence experiments reveal how DEAD-box proteins unfold structured RNAs to promote conformational transitions and refolding to the native functional state.

DEAD-box helicase proteins accelerate folding and rearrangements of highly structured RNAs and RNA–protein complexes (RNPs) in many essential cellular processes. Although DEAD-box proteins have been shown to use ATP to unwind short RNA helices, it is not known how they disrupt RNA tertiary structure. Here, we use single molecule fluorescence to show that the DEAD-box protein CYT-19 disrupts tertiary structure in a group I intron using a helix capture mechanism. CYT-19 binds to a helix within the structured RNA only after the helix spontaneously loses its tertiary contacts, and then CYT-19 uses ATP to unwind the helix, liberating the product strands. Ded1, a multifunctional yeast DEAD-box protein, gives analogous results with small but reproducible differences that may reflect its in vivo roles. The requirement for spontaneous dynamics likely targets DEAD-box proteins toward less stable RNA structures, which are likely to experience greater dynamic fluctuations, and provides a satisfying explanation for previous correlations between RNA stability and CYT-19 unfolding efficiency. Biologically, the ability to sense RNA stability probably biases DEAD-box proteins to act preferentially on less stable misfolded structures and thereby to promote native folding while minimizing spurious interactions with stable, natively folded RNAs. In addition, this straightforward mechanism for RNA remodeling does not require any specific structural environment of the helicase core and is likely to be relevant for DEAD-box proteins that promote RNA rearrangements of RNP complexes including the spliceosome and ribosome.

In addition to carrying genetic information from DNA to protein, RNAs function in many essential cellular processes. This often requires the RNA to form a specific three-dimensional structure, and some functions require cycling between multiple structures. However, RNAs have a strong propensity to become trapped in nonfunctional, misfolded structures. Due to the intrinsic stability of local structure for RNA, these misfolded species can be long-lived and therefore accumulate. ATP-dependent RNA chaperone proteins called DEAD-box proteins are known to promote native RNA folding by disrupting misfolded structures. Although these proteins can unwind short RNA helices, the mechanism by which they act upon higher order tertiary contacts is unknown. Our current work shows that DEAD-box proteins capture transiently exposed RNA helices, preventing any tertiary contacts from reforming and potentially destabilizing the global RNA architecture. Helix unwinding by the DEAD-box protein then allows the product RNA strands to form new contacts. This helix capture mechanism for manipulation of RNA tertiary structure does not require a specific binding motif or structural environment and may be general for DEAD-box helicase proteins that act on structured RNAs.

DEAD-box protein CYT-19 is activated by exposed helices in a group I intron RNA

DEAD-box proteins are nonprocessive RNA helicases and can function as RNA chaperones, but the mechanisms of their chaperone activity remain incompletely understood. The Neurospora crassa DEAD-box protein CYT-19 is a mitochondrial RNA chaperone that promotes group I intron splicing and has been shown to resolve misfolded group I intron structures, allowing them to refold. Building on previous results, here we use a series of tertiary contact mutants of the Tetrahymena group I intron ribozyme to demonstrate that the efficiency of CYT-19-mediated unfolding of the ribozyme is tightly linked to global RNA tertiary stability. Efficient unfolding of destabilized ribozyme variants is accompanied by increased ATPase activity of CYT-19, suggesting that destabilized ribozymes provide more productive interaction opportunities. The strongest ATPase stimulation occurs with a ribozyme that lacks all five tertiary contacts and does not form a compact structure, and small-angle X-ray scattering indicates that ATPase activity tracks with ribozyme compactness. Further, deletion of three helices that are prominently exposed in the folded structure decreases the ATPase stimulation by the folded ribozyme. Together, these results lead to a model in which CYT-19, and likely related DEAD-box proteins, rearranges complex RNA structures by preferentially interacting with and unwinding exposed RNA secondary structure. Importantly, this mechanism could bias DEAD-box proteins to act on misfolded RNAs and ribonucleoproteins, which are likely to be less compact and more dynamic than their native counterparts.

Detecting RNA tertiary folding by sedimentation velocity analytical ultracentrifugation

Analytical Ultracentrifugation (AUC) is a highly sensitive technique for detecting global conformational features of biological molecules and molecular interactions in solution. When operated in a sedimentation velocity (SV) recording mode, it reports precisely on the hydrodynamic properties of a molecule, including its sedimentation and diffusion coefficients, which can be used to calculate its hydrated radius, as well as, to estimate its global shape. This chapter describes the application of SV-AUC to the detection of global conformational changes accompanying equilibrium counterion induced tertiary folding of structured RNA molecules. A brief theoretical background is provided at the beginning, aimed at familiarizing the readers with the operational principle of the technique; then, a detailed set of instructions is provided on how to design, conduct, and analyze the data from an equilibrium RNA folding experiment, using SV-AUC.

Folding pathways of the Tetrahymena ribozyme

Like many structured RNAs, the Tetrahymena group I intron ribozyme folds through multiple pathways and intermediates. Under standard conditions in vitro, a small fraction reaches the native state (N) with kobs ≈ 0.6 min(-1), while the remainder forms a long-lived misfolded conformation (M) thought to differ in topology. These alternative outcomes reflect a pathway that branches late in folding, after disruption of a trapped intermediate (Itrap). Here we use catalytic activity to probe the folding transitions from Itrap to the native and misfolded states. We show that mutations predicted to weaken the core helix P3 do not increase the rate of folding from Itrap but they increase the fraction that reaches the native state rather than forming the misfolded state. Thus, P3 is disrupted during folding to the native state but not to the misfolded state, and P3 disruption occurs after the rate-limiting step. Interestingly, P3-strengthening mutants also increase native folding. Additional experiments show that these mutants are rapidly committed to folding to the native state, although they reach the native state with approximately the same rate constant as the wild-type ribozyme (~1 min(-1)). Thus, the P3-strengthening mutants populate a distinct pathway that includes at least one intermediate but avoids the M state, most likely because P3 and the correct topology are formed early. Our results highlight multiple pathways in RNA folding and illustrate how kinetic competitions between rapid events can have long-lasting effects because the "choice" is enforced by energy barriers that grow larger as folding progresses.

Single-molecule studies of riboswitch folding

The folding dynamics of riboswitches are central to their ability to modulate gene expression in response to environmental cues. In most cases, a structural competition between the formation of a ligand-binding aptamer and an expression platform (or some other competing off-state) determines the regulatory outcome. Here, we review single-molecule studies of riboswitch folding and function, predominantly carried out using single-molecule FRET or optical trapping approaches. Recent results have supplied new insights into riboswitch folding energy landscapes, the mechanisms of ligand binding, the roles played by divalent ions, the applicability of hierarchical folding models, and kinetic vs. thermodynamic control schemes. We anticipate that future work, based on improved data sets and potentially combining multiple experimental techniques, will enable the development of more complete models for complex RNA folding processes. This article is part of a Special Issue entitled: Riboswitches.

RNA helicase proteins as chaperones and remodelers

Superfamily 2 helicase proteins are ubiquitous in RNA biology and have an extraordinarily broad set of functional roles. Central among these roles are the promotion of rearrangements of structured RNAs and the remodeling of ribonucleoprotein complexes (RNPs), allowing formation of native RNA structure or progression through a functional cycle of structures. Although all superfamily 2 helicases share a conserved helicase core, they are divided evolutionarily into several families, and it is principally proteins from three families, the DEAD-box, DEAH/RHA, and Ski2-like families, that function to manipulate structured RNAs and RNPs. Strikingly, there are emerging differences in the mechanisms of these proteins, both between families and within the largest family (DEAD-box), and these differences appear to be tuned to their RNA or RNP substrates and their specific roles. This review outlines basic mechanistic features of the three families and surveys individual proteins and the current understanding of their biological substrates and mechanisms.

Structural dynamics of a mitochondrial tRNA possessing weak thermodynamic stability

Folding dynamics are ubiquitously involved in controlling the multivariate functions of RNAs. While the high thermodynamic stabilities of some RNAs favor purely native states at equilibrium, it is unclear whether weakly stable RNAs exist in random, partially folded states or sample well-defined, globally folded conformations. Using a folding assay that precisely tracks the formation of native aminoacylable tRNA, we show that the folding of a weakly stable human mitochondrial (hmt) leucine tRNA is hierarchical with a distinct kinetic folding intermediate. The stabilities of the native and intermediate conformers are separated by only about 1.2 kcal/mol, and the species are readily interconvertible. Comparison of folding dynamics between unmodified and fully modified tRNAs reveals that post-transcriptional modifications produce a more constrained native structure that does not sample intermediate conformations. These structural dynamics may thus be crucial for recognition by some modifying enzymes in vivo, especially those targeting the globular core region, by allowing access to pretransition state conformers. Reduced conformational sampling of the native, modified tRNAs could then permit improved performance in downstream processes of translation. More generally, weak stabilities of small RNAs that fold in the absence of chaperone proteins may facilitate conformational switching that is central to biological function.

Mapping the kinetic barriers of a Large RNA molecule's folding landscape

The folding of linear polymers into discrete three-dimensional structures is often required for biological function. The formation of long-lived intermediates is a hallmark of the folding of large RNA molecules due to the ruggedness of their energy landscapes. The precise thermodynamic nature of the barriers (whether enthalpic or entropic) that leads to intermediate formation is still poorly characterized in large structured RNA molecules. A classic approach to analyzing kinetic barriers are temperature dependent studies analyzed with Eyring's transition state theory. We applied Eyring's theory to time-resolved hydroxyl radical (•OH) footprinting kinetics progress curves collected at eight temperature from 21.5 °C to 51 °C to characterize the thermodynamic nature of folding intermediate formation for the Mg(2+)-mediated folding of the Tetrahymena thermophila group I ribozyme. A common kinetic model configuration describes this RNA folding reaction over the entire temperature range studied consisting of primary (fast) transitions to misfolded intermediates followed by much slower secondary transitions, consistent with previous studies. Eyring analysis reveals that the primary transitions are moderate in magnitude and primarily enthalpic in nature. In contrast, the secondary transitions are daunting in magnitude and entropic in nature. The entropic character of the secondary transitions is consistent with structural rearrangement of the intermediate species to the final folded form. This segregation of kinetic control reveals distinctly different molecular mechanisms during the two stages of RNA folding and documents the importance of entropic barriers to defining rugged RNA folding landscapes.

RNA catalytic activity as a probe of chaperone-mediated RNA folding

For structured RNAs that possess catalytic activity, this activity provides a powerful probe for measuring the progress of folding and the effects of RNA chaperone proteins on the folding rate. The crux of this approach is that only the natively folded RNA is able to perform the catalytic reaction. This method can provide a quantitative measure of the fraction of native RNA over time, and it can readily distinguish the native state from all misfolded conformations. Here we describe an activity-based method measuring native folding of ribozymes derived from self-splicing group I introns, and we show how the assay can be used to monitor acceleration of native folding by DEAD-box RNA helicase proteins that function as general RNA chaperones. By measuring the amount of substrate that is converted to product in a rapid first turnover, we describe how to determine the fraction of the ribozyme population that is present in the native state. Further, we describe how to perform a two-stage or discontinuous assay in which folding proceeds in stage one and then solution conditions are changed in stage two to permit catalytic activity and block further folding. This protocol allows folding to be followed under a broad range of solution conditions, including those that do not support catalytic activity, and facilitates studies of chaperone proteins.

Evolutionary evidence for alternative structure in RNA sequence co-variation

Sequence conservation and co-variation of base pairs are hallmarks of structured RNAs. For certain RNAs (e.g. riboswitches), a single sequence must adopt at least two alternative secondary structures to effectively regulate the message. If alternative secondary structures are important to the function of an RNA, we expect to observe evolutionary co-variation supporting multiple conformations. We set out to characterize the evolutionary co-variation supporting alternative conformations in riboswitches to determine the extent to which alternative secondary structures are conserved. We found strong co-variation support for the terminator, P1, and anti-terminator stems in the purine riboswitch by extending alignments to include terminator sequences. When we performed Boltzmann suboptimal sampling on purine riboswitch sequences with terminators we found that these sequences appear to have evolved to favor specific alternative conformations. We extended our analysis of co-variation to classic alignments of group I/II introns, tRNA, and other classes of riboswitches. In a majority of these RNAs, we found evolutionary evidence for alternative conformations that are compatible with the Boltzmann suboptimal ensemble. Our analyses suggest that alternative conformations are selected for and thus likely play functional roles in even the most structured of RNAs.

RNA (Ribonucleic Acid) is a messenger of genetic information, master regulator, and catalyst in the cell. To carry out its function, RNA can fold into complex three-dimensional structures. Certain classes of RNAs, called riboswitches, adopt at least two alternative structures to act as a switch. We set out to detect the evolutionary signal for alternative structures in riboswitches as we hypothesize that these RNA sequences must have evolved to allow both conformations. We find that indeed such signals exist when we compare the sequences of riboswitches from multiple species. When we extend this analysis to other RNA regulators in the cell that are not thought of as switches, we detect equivalent evolutionary support for alternative structures. Viewed through the lens of evolutionary structure conservation RNA sequences appear to have adapted to adopt multiple conformations.

The long-range P3 helix of the Tetrahymena ribozyme is disrupted during folding between the native and misfolded conformations

RNAs are prone to misfolding, but how misfolded structures are formed and resolved remains incompletely understood. The Tetrahymena group I intron ribozyme folds in vitro to a long-lived misfolded conformation (M) that includes extensive native structure but is proposed to differ in topology from the native state (N). A leading model predicts that exchange of the topologies requires unwinding of the long-range, core helix P3, despite the presence of P3 in both conformations. To test this model, we constructed 16 mutations to strengthen or weaken P3. Catalytic activity and in-line probing showed that nearly all of the mutants form the M state before folding to N. The P3-weakening mutations accelerated refolding from M (3- to 30-fold) and the P3-strengthening mutations slowed refolding (6- to 1400-fold), suggesting that P3 indeed unwinds transiently. Upon depletion of Mg(2+), the mutations had analogous effects on unfolding from N to intermediates that subsequently fold to M. The magnitudes for the P3-weakening mutations were larger than in refolding from M, and small-angle X-ray scattering showed that the ribozyme expands rapidly to intermediates from which P3 is disrupted subsequently. These results are consistent with previous results indicating unfolding of native peripheral structure during refolding from M, which probably permits rearrangement of the core. Together, our results demonstrate that exchange of the native and misfolded conformations requires loss of a core helix in addition to peripheral structure. Further, the results strongly suggest that misfolding arises from a topological error within the ribozyme core, and a specific topology is proposed.

Understanding the role of three-dimensional topology in determining the folding intermediates of group I introns

Many RNA molecules exert their biological function only after folding to unique three-dimensional structures. For long, noncoding RNA molecules, the complexity of finding the native topology can be a major impediment to correct folding to the biologically active structure. An RNA molecule may fold to a near-native structure but not be able to continue to the correct structure due to a topological barrier such as crossed strands or incorrectly stacked helices. Achieving the native conformation thus requires unfolding and refolding, resulting in a long-lived intermediate. We investigate the role of topology in the folding of two phylogenetically related catalytic group I introns, the Twort and Azoarcus group I ribozymes. The kinetic models describing the Mg(2+)-mediated folding of these ribozymes were previously determined by time-resolved hydroxyl (∙OH) radical footprinting. Two intermediates formed by parallel intermediates were resolved for each RNA. These data and analytical ultracentrifugation compaction analyses are used herein to constrain coarse-grained models of these folding intermediates as we investigate the role of nonnative topology in dictating the lifetime of the intermediates. Starting from an ensemble of unfolded conformations, we folded the RNA molecules by progressively adding native constraints to subdomains of the RNA defined by the ∙OH time-progress curves to simulate folding through the different kinetic pathways. We find that nonnative topologies (arrangement of helices) occur frequently in the folding simulations despite using only native constraints to drive the reaction, and that the initial conformation, rather than the folding pathway, is the major determinant of whether the RNA adopts nonnative topology during folding. From these analyses we conclude that biases in the initial conformation likely determine the relative flux through parallel RNA folding pathways.

An in vitro-selected RNA receptor for the GAAC loop: modular receptor for non-GNRA-type tetraloop

Although artificial RNA motifs that can functionally replace the GNRA/receptor interaction, a class of RNA-RNA interacting motifs, were isolated from RNA libraries and used to generate designer RNA structures, receptors for non-GNRA tetraloops have not been found in nature or selected from RNA libraries. In this study, we report successful isolation of a receptor motif interacting with GAAC, a non-GNRA tetraloop, from randomized sequences embedded in a catalytic RNA. Biochemical characterization of the GAAC/receptor interacting motif within three structural contexts showed its binding affinity, selectivity and structural autonomy. The motif has binding affinity comparable with that of a GNRA/receptor, selectivity orthogonal to GNRA/receptors and structural autonomy even in a large RNA context. These features would be advantageous for usage of the motif as a building block for designer RNAs. The isolated motif can also be used as a query sequence to search for unidentified naturally occurring GANC receptor motifs.

Toward a molecular understanding of RNA remodeling by DEAD-box proteins

DEAD-box proteins are superfamily 2 helicases that function in all aspects of RNA metabolism. They employ ATP binding and hydrolysis to generate tight, yet regulated RNA binding, which is used to unwind short RNA helices non-processively and promote structural transitions of RNA and RNA-protein substrates. In the last few years, substantial progress has been made toward a detailed, quantitative understanding of the structural and biochemical properties of DEAD-box proteins. Concurrently, progress has been made toward a physical understanding of the RNA rearrangements and folding steps that are accelerated by DEAD-box proteins in model systems. Here, we review the recent progress on both of these fronts, focusing on the mitochondrial DEAD-box proteins Mss116 and CYT-19 and their mechanisms in promoting the splicing of group I and group II introns.

DEAD-box proteins as RNA helicases and chaperones

DEAD-box proteins are ubiquitous in RNA-mediated processes and function by coupling cycles of ATP binding and hydrolysis to changes in affinity for single-stranded RNA. Many DEAD-box proteins use this basic mechanism as the foundation for a version of RNA helicase activity, efficiently separating the strands of short RNA duplexes in a process that involves little or no translocation. This activity, coupled with mechanisms to direct different DEAD-box proteins to their physiological substrates, allows them to promote RNA folding steps and rearrangements and to accelerate remodeling of RNA–protein complexes. This review will describe the properties of DEAD-box proteins as RNA helicases and the current understanding of how the energy from ATPase activity is used to drive the separation of RNA duplex strands. It will then describe how the basic biochemical properties allow some DEAD-box proteins to function as chaperones by promoting RNA folding reactions, with a focus on the self-splicing group I and group II intron RNAs.

RNA catalysis as a probe for chaperone activity of DEAD-box helicases

DEAD-box proteins are vitally important to cellular processes and make up the largest class of helicases. Many DEAD-box proteins function as RNA chaperones by accelerating structural transitions of RNA, which can result in the resolution of misfolded conformers or conversion between functional structures. While the biological importance of chaperone proteins is clear, their mechanisms are incompletely understood. Here, we illustrate how the catalytic activity of certain RNAs can be used to measure RNA chaperone activity. By measuring the amount of substrate converted to product, the fraction of catalytically active molecules is measured over time, providing a quantitative measure of the formation or loss of native RNA. The assays are described with references to group I and group II introns and their ribozyme derivatives, and examples are included that illustrate potential complications and indicate how catalytic activity measurements can be combined with physical approaches to gain insights into the mechanisms of DEAD-box proteins as RNA chaperones.

Structure-function analysis from the outside in: long-range tertiary contacts in RNA exhibit distinct catalytic roles

The conserved catalytic core of the Tetrahymena group I ribozyme is encircled by peripheral elements. We have conducted a detailed structure-function study of the five long-range tertiary contacts that fasten these distal elements together. Mutational ablation of each of the tertiary contacts destabilizes the folded ribozyme, indicating a role of the peripheral elements in overall stability. Once folded, three of the five tertiary contact mutants exhibit defects in overall catalysis that range from 20- to 100-fold. These and the subsequent results indicate that the structural ring of peripheral elements does not act as a unitary element; rather, individual connections have distinct roles as further revealed by kinetic and thermodynamic dissection of the individual reaction steps. Ablation of P14 or the metal ion core/metal ion core receptor (MC/MCR) destabilizes docking of the substrate-containing P1 helix into tertiary interactions with the ribozyme's conserved core. In contrast, ablation of the L9/P5 contact weakens binding of the guanosine nucleophile by slowing its association, without affecting P1 docking. The P13 and tetraloop/tetraloop receptor (TL/TLR) mutations had little functional effect and small, local structural changes, as revealed by hydroxyl radical footprinting, whereas the P14, MC/MCR, and L9/P5 mutants show structural changes distal from the mutation site. These changes extended into regions of the catalytic core involved in docking or guanosine binding. Thus, distinct allosteric pathways couple the long-range tertiary contacts to functional sites within the conserved core. This modular functional specialization may represent a fundamental strategy in RNA structure-function interrelationships.

Effects of a protecting osmolyte on the ion atmosphere surrounding DNA duplexes

Osmolytes are small, chemically diverse, organic solutes that function as an essential component of cellular stress response. Protecting osmolytes enhance protein stability via preferential exclusion, and nonprotecting osmolytes, such as urea, destabilize protein structures. Although much is known about osmolyte effects on proteins, less is understood about osmolyte effects on nucleic acids and their counterion atmospheres. Nonprotecting osmolytes destabilize nucleic acid structures, but effects of protecting osmolytes depend on numerous factors including the type of nucleic acid and the complexity of the functional fold. To begin quantifying protecting osmolyte effects on nucleic acid interactions, we used small-angle X-ray scattering (SAXS) techniques to monitor DNA duplexes in the presence of sucrose. This protecting osmolyte is a commonly used contrast matching agent in SAXS studies of protein-nucleic acid complexes; thus, it is important to characterize interaction changes induced by sucrose. Measurements of interactions between duplexes showed no dependence on the presence of up to 30% sucrose, except under high Mg(2+) conditions where stacking interactions were disfavored. The number of excess ions associated with DNA duplexes, reported by anomalous small-angle X-ray scattering (ASAXS) experiments, was sucrose independent. Although protecting osmolytes can destabilize secondary structures, our results suggest that ion atmospheres of individual duplexes remain unperturbed by sucrose.

The Azoarcus group I intron ribozyme misfolds and is accelerated for refolding by ATP-dependent RNA chaperone proteins

Structured RNAs traverse complex energy landscapes that include valleys representing misfolded intermediates. In Neurospora crassa and Saccharomyces cerevisiae, efficient splicing of mitochondrial group I and II introns requires the DEAD box proteins CYT-19 and Mss116p, respectively, which promote folding transitions and function as general RNA chaperones. To test the generality of RNA misfolding and the activities of DEAD box proteins in vitro, here we measure native folding of a small group I intron ribozyme from the bacterium Azoarcus by monitoring its catalytic activity. To develop this assay, we first measure cleavage of an oligonucleotide substrate by the prefolded ribozyme. Substrate cleavage is rate-limited by binding and is readily reversible, with an internal equilibrium near unity, such that the amount of product observed is less than the amount of native ribozyme. We use this assay to show that approximately half of the ribozyme folds readily to the native state, whereas the other half forms an intermediate that transitions slowly to the native state. This folding transition is accelerated by urea and increased temperature and slowed by increased Mg(2+) concentration, suggesting that the intermediate is misfolded and must undergo transient unfolding during refolding to the native state. CYT-19 and Mss116p accelerate refolding in an ATP-dependent manner, presumably by disrupting structure in the intermediate. These results highlight the tendency of RNAs to misfold, underscore the roles of CYT-19 and Mss116p as general RNA chaperones, and identify a refolding transition for further dissection of the roles of DEAD box proteins in RNA folding.

Understanding the errors of SHAPE-directed RNA structure modeling

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 using six molecules for which crystallographic data are available: tRNA(phe) and 5S rRNA from Escherichia 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 the insufficient information content of these RNAs' SHAPE data, as evaluated by a nonparametric bootstrapping analysis. Beyond these benchmark cases, bootstrapping suggests a low level of confidence (<50%) in the majority of helices in a previously proposed SHAPE-directed model for the HIV-1 RNA genome. Thus, SHAPE-directed RNA modeling is not always unambiguous, and helix-by-helix confidence estimates, as described herein, may be critical for interpreting results from this powerful methodology.

ATP-dependent roles of the DEAD-box protein Mss116p in group II intron splicing in vitro and in vivo

The yeast DEAD-box protein Mss116p functions as a general RNA chaperone in splicing mitochondrial group I and group II introns. For most of its functions, Mss116p is thought to use ATP-dependent RNA unwinding to facilitate RNA structural transitions, but it has been suggested to assist in the folding of one group II intron (aI5γ) primarily by stabilizing a folding intermediate. Here we compare three aI5γ constructs: one with long exons, one with short exons, and a ribozyme construct lacking exons. The long exons result in slower splicing, suggesting that they misfold and/or stabilize nonnative intronic structures. Nevertheless, Mss116p acceleration of all three constructs depends on ATP and is inhibited by mutations that compromise RNA unwinding, suggesting similar mechanisms. Results of splicing assays and a new two-stage assay that separates ribozyme folding and catalysis indicate that maximal folding of all three constructs by Mss116p requires ATP-dependent RNA unwinding. ATP-independent activation is appreciable for only a subpopulation of the minimal ribozyme construct and not for constructs containing exons. As expected for a general RNA chaperone, Mss116p can also disrupt the native ribozyme, which can refold after Mss116p removal. Finally, using yeast strains with mitochondrial DNA containing only the single intron aI5γ, we show that Mss116p mutants promote splicing in vivo to degrees that correlate with their residual ATP-dependent RNA-unwinding activities. Together, our results indicate that, although DEAD-box proteins play multiple roles in RNA folding, the physiological function of Mss116p in aI5γ splicing includes a requirement for ATP-dependent local unfolding, allowing the conversion of nonfunctional RNA structure into functional RNA structure.

RNA molecules with conserved catalytic cores but variable peripheries fold along unique energetically optimized pathways

Functional and kinetic constraints must be efficiently balanced during the folding process of all biopolymers. To understand how homologous RNA molecules with different global architectures fold into a common core structure we determined, under identical conditions, the folding mechanisms of three phylogenetically divergent group I intron ribozymes. These ribozymes share a conserved functional core defined by topologically equivalent tertiary motifs but differ in their primary sequence, size, and structural complexity. Time-resolved hydroxyl radical probing of the backbone solvent accessible surface and catalytic activity measurements integrated with structural-kinetic modeling reveal that each ribozyme adopts a unique strategy to attain the conserved functional fold. The folding rates are not dictated by the size or the overall structural complexity, but rather by the strength of the constituent tertiary motifs which, in turn, govern the structure, stability, and lifetime of the folding intermediates. A fundamental general principle of RNA folding emerges from this study: The dominant folding flux always proceeds through an optimally structured kinetic intermediate that has sufficient stability to act as a nucleating scaffold while retaining enough conformational freedom to avoid kinetic trapping. Our results also suggest a potential role of naturally selected peripheral A-minor interactions in balancing RNA structural stability with folding efficiency.

Sharing and archiving nucleic acid structure mapping data

Nucleic acids are particularly amenable to structural characterization using chemical and enzymatic probes. Each individual structure mapping experiment reveals specific information about the structure and/or dynamics of the nucleic acid. Currently, there is no simple approach for making these data publically available in a standardized format. We therefore developed a standard for reporting the results of single nucleotide resolution nucleic acid structure mapping experiments, or SNRNASMs. We propose a schema for sharing nucleic acid chemical probing data that uses generic public servers for storing, retrieving, and searching the data. We have also developed a consistent nomenclature (ontology) within the Ontology of Biomedical Investigations (OBI), which provides unique identifiers (termed persistent URLs, or PURLs) for classifying the data. Links to standardized data sets shared using our proposed format along with a tutorial and links to templates can be found at http://snrnasm.bio.unc.edu.

Enhanced specificity against misfolding in a thermostable mutant of the Tetrahymena ribozyme

Structured RNAs encode native conformations that are more stable than the vast ensembles of alternative conformations, but how this specificity is evolved is incompletely understood. Here we show that a variant of the Tetrahymena group I intron ribozyme that was generated previously by in vitro selection for enhanced thermostability also displays modestly enhanced specificity against a stable misfolded structure that is globally similar to the native state, despite the absence of selective pressure to increase the energy gap between these structures. The enhanced specificity for native folding arises from mutations in two nucleotides that are close together in space in the native structure, and additional experiments show that these two mutations do not affect the stability of the misfolded conformation relative to the largely unstructured transition state ensemble for interconversion between the native and misfolded conformers. Thus, they selectively stabilize the native state, presumably by strengthening a local tertiary contact network that cannot form in the misfolded conformation. The stabilization is larger in the presence of the peripheral element P5abc, suggesting that cooperative tertiary structure formation plays a key role in the enhanced stability. The increased specificity in the absence of explicit selection suggests that the large energy gap in the wild-type RNA may have arisen analogously, a consequence of selective pressure for stability of the functional structure. More generally, the structural rigidity and intricate networks of contacts in structured RNAs may allow them to evolve substantial structural specificity without explicit negative selection, even against closely related alternative structures.