New pathways in folding of the Tetrahymena group I RNA enzyme

A high-throughput search for intracellular factors that affect RNA folding identifies E. coli proteins PepA and YagL as RNA chaperones that promote RNA remodelling

General RNA chaperones are RNA-binding proteins (RBPs) that interact transiently and non-specifically with RNA substrates and assist in their folding into their native state. In bacteria, these chaperones impact both coding and non-coding RNAs and are particularly important for large, structured RNAs which are prone to becoming kinetically trapped in misfolded states. Currently, due to the limited number of well-characterized examples and the lack of a consensus structural or sequence motif, it is difficult to identify general RNA chaperones in bacteria. Here, we adapted a previously published in vivo RNA regional accessibility probing assay to screen genome wide for intracellular factors in E. coli affecting RNA folding, among which we aimed to uncover novel RNA chaperones. Through this method, we identified eight proteins whose deletion gives changes in regional accessibility within the exogenously expressed Tetrahymena group I intron ribozyme. Furthermore, we purified and measured in vitro properties of two of these proteins, YagL and PepA, which were especially attractive as general chaperone candidates. We showed that both proteins bind RNA and that YagL accelerates native refolding of the ribozyme from a long-lived misfolded state. Further dissection of YagL showed that a putative helix-turn-helix (HTH) domain is responsible for most of its RNA-binding activity, but only the full protein shows chaperone activity. Altogether, this work expands the current repertoire of known general RNA chaperones in bacteria.

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.

Structural Expansion of Catalytic RNA Nanostructures through Oligomerization of a Cyclic Trimer of Engineered Ribozymes

The multimolecular assembly of three-dimensionally structured proteins forms their quaternary structures, some of which have high geometric symmetry. The size and complexity of protein quaternary structures often increase in a hierarchical manner, with simpler, smaller structures serving as units for larger quaternary structures. In this study, we exploited oligomerization of a ribozyme cyclic trimer to achieve larger ribozyme-based RNA assembly. By installing kissing loop (KL) interacting units to one-, two-, or three-unit RNA molecules in the ribozyme trimer, we constructed dimers, open-chain oligomers, and branched oligomers of ribozyme trimer units. One type of open-chain oligomer preferentially formed a closed tetramer containing 12 component RNAs to provide 12 ribozyme units. We also observed large assembly of ribozyme trimers, which reached 1000 nm in size.

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.

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.

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.

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.

Iterative annealing mechanism explains the functions of the GroEL and RNA chaperones

Molecular chaperones are ATP-consuming machines, which facilitate the folding of proteins and RNA molecules that are kinetically trapped in misfolded states. Unassisted folding occurs by the kinetic partitioning mechanism according to which folding to the native state, with low probability as well as misfolding to one of the many metastable states, with high probability, occur rapidly. GroEL is an all-purpose stochastic machine that assists misfolded substrate proteins to fold. The RNA chaperones such as CYT-19, which are ATP-consuming enzymes, help the folding of ribozymes that get trapped in metastable states for long times. GroEL does not interact with the folded proteins but CYT-19 disrupts both the folded and misfolded ribozymes. The structures of GroEL and RNA chaperones are strikingly different. Despite these differences, the iterative annealing mechanism (IAM) quantitatively explains all the available experimental data for assisted folding of proteins and ribozymes. Driven by ATP binding and hydrolysis and GroES binding, GroEL undergoes a catalytic cycle during which it samples three allosteric states, T (apo), R (ATP bound), and R^″ (ADP bound). Analyses of the experimental data show that the efficiency of the GroEL-GroES machinery and mutants is determined by the resetting rate k _{R ″  → T} , which is largest for the wild-type (WT) GroEL. Generalized IAM accurately predicts the folding kinetics of Tetrahymena ribozyme and its variants. Chaperones maximize the product of the folding rate and the steady-state native state fold by driving the substrates out of equilibrium. Neither the absolute yield nor the folding rate is optimized.

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.

Quantitative tests of a reconstitution model for RNA folding thermodynamics and kinetics

Decades of study of the architecture and function of structured RNAs have led to the perspective that RNA tertiary structure is modular, made of locally stable domains that retain their structure across RNAs. We formalize a hypothesis inspired by this modularity-that RNA folding thermodynamics and kinetics can be quantitatively predicted from separable energetic contributions of the individual components of a complex RNA. This reconstitution hypothesis considers RNA tertiary folding in terms of ΔG_align, the probability of aligning tertiary contact partners, and ΔG_tert, the favorable energetic contribution from the formation of tertiary contacts in an aligned state. This hypothesis predicts that changes in the alignment of tertiary contacts from different connecting helices and junctions (ΔG_HJH) or from changes in the electrostatic environment (ΔG_+/-) will not affect the energetic perturbation from a mutation in a tertiary contact (ΔΔG_tert). Consistent with these predictions, single-molecule FRET measurements of folding of model RNAs revealed constant ΔΔG_tert values for mutations in a tertiary contact embedded in different structural contexts and under different electrostatic conditions. The kinetic effects of these mutations provide further support for modular behavior of RNA elements and suggest that tertiary mutations may be used to identify rate-limiting steps and dissect folding and assembly pathways for complex RNAs. Overall, our model and results are foundational for a predictive understanding of RNA folding that will allow manipulation of RNA folding thermodynamics and kinetics. Conversely, the approaches herein can identify cases where an independent, additive model cannot be applied and so require additional investigation.

RNA Structural Modules Control the Rate and Pathway of RNA Folding and Assembly

Structured RNAs fold through multiple pathways, but we have little understanding of the molecular features that dictate folding pathways and determine rates along a given pathway. Here, we asked whether folding of a complex RNA can be understood from its structural modules. In a two-piece version of the Tetrahymena group I ribozyme, the separated P5abc subdomain folds to local native secondary and tertiary structure in a linked transition and assembles with the ribozyme core via three tertiary contacts: a kissing loop (P14), a metal core-receptor interaction, and a tetraloop-receptor interaction, the first two of which are expected to depend on native P5abc structure from the local transition. Native gel, NMR, and chemical footprinting experiments showed that mutations that destabilize the native P5abc structure slowed assembly up to 100-fold, indicating that P5abc folds first and then assembles with the core by conformational selection. However, rate decreases beyond 100-fold were not observed because an alternative pathway becomes dominant, with nonnative P5abc binding the core and then undergoing an induced-fit rearrangement. P14 is formed in the rate-limiting step along the conformational selection pathway but after the rate-limiting step along the induced-fit pathway. Strikingly, the assembly rate along the conformational selection pathway resembles that of an isolated kissing loop similar to P14, and the rate along the induced-fit pathway resembles that of an isolated tetraloop-receptor interaction. Our results indicate substantial modularity in RNA folding and assembly and suggest that these processes can be understood in terms of underlying structural modules.

Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch

Intermediates play important roles in RNA folding but can be difficult to characterize when short-lived or not significantly populated. By combining (15)N relaxation dispersion NMR with chemical probing, we visualized a fast (kex=k1+k-1≈423 s(-1)) secondary structural switch directed towards a low-populated (∼3%) partially folded intermediate in tertiary folding of the P5abc subdomain of the 'Tetrahymena' group I intron ribozyme. The secondary structure switch changes the base-pairing register across the P5c hairpin, creating a native-like structure, and occurs at rates of more than two orders of magnitude faster than tertiary folding. The switch occurs robustly in the absence of tertiary interactions, Mg(2+) or even when the hairpin is excised from the three-way junction. Fast, highly modular secondary structural switches may be quite common during RNA tertiary folding where they may help smoothen the folding landscape by allowing folding to proceed efficiently via additional pathways.

A kinetic and thermodynamic framework for the Azoarcus group I ribozyme reaction

Determination of quantitative thermodynamic and kinetic frameworks for ribozymes derived from the Azoarcus group I intron and comparisons to their well-studied analogs from the Tetrahymena group I intron reveal similarities and differences between these RNAs. The guanosine (G) substrate binds to the Azoarcus and Tetrahymena ribozymes with similar equilibrium binding constants and similar very slow association rate constants. These and additional literature observations support a model in which the free ribozyme is not conformationally competent to bind G and in which the probability of assuming the binding-competent state is determined by tertiary interactions of peripheral elements. As proposed previously, the slow binding of guanosine may play a role in the specificity of group I intron self-splicing, and slow binding may be used analogously in other biological processes. The internal equilibrium between ribozyme-bound substrates and products is similar for these ribozymes, but the Azoarcus ribozyme does not display the coupling in the binding of substrates that is observed with the Tetrahymena ribozyme, suggesting that local preorganization of the active site and rearrangements within the active site upon substrate binding are different for these ribozymes. Our results also confirm the much greater tertiary binding energy of the 5'-splice site analog with the Azoarcus ribozyme, binding energy that presumably compensates for the fewer base-pairing interactions to allow the 5'-exon intermediate in self splicing to remain bound subsequent to 5'-exon cleavage and prior to exon ligation. Most generally, these frameworks provide a foundation for design and interpretation of experiments investigating fundamental properties of these and other 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.

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.

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.

The kinetics of ribozyme cleavage: a tool to analyze RNA folding as a function of catalysis

As catalytically active RNAs, ribozymes can be characterized by kinetic measurements similar to classical enzyme kinetics. However, in contrast to standard protein enzymes, for which reactions can usually be started by mixing the enzyme with its substrate, ribozymes are typically self-cleaving. The reaction has to be initiated by folding the RNA into its active conformation. Thus, ribozyme kinetics are influenced by both folding and catalytic components and often enable indirect observation of RNA folding. Here, I describe how to obtain quantitative ribozyme cleavage data via denaturing polyacrylamide gel electrophoresis (PAGE) of radioactively labeled in vitro transcripts and discuss general considerations for subsequent kinetic analysis.

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.

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.

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.

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.

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.

Roles of DEAD-box proteins in RNA and RNP Folding

RNAs and RNA-protein complexes (RNPs) traverse rugged energy landscapes as they fold to their native structures, and many continue to undergo conformational rearrangements as they function. Due to the inherent stability of local RNA structure, proteins are required to assist with RNA conformational transitions during initial folding and in exchange between functional structures. DEAD-box proteins are superfamily 2 RNA helicases that are ubiquitously involved in RNA-mediated processes. Some of these proteins use an ATP-dependent cycle of conformational changes to disrupt RNA structure nonprocessively, accelerating structural transitions of RNAs and RNPs in a manner that bears a strong resemblance to the activities of certain groups of protein chaperones. This review summarizes recent work using model substrates and tractable self-splicing intron RNAs, which has given new insights into how DEAD-box proteins promote RNA folding steps and conformational transitions, and it summarizes recent progress in identifying sites and mechanisms of DEAD-box protein activity within more complex cellular targets.

Multiple unfolding events during native folding of the Tetrahymena group I ribozyme

Despite the ubiquitous nature of misfolded intermediates in RNA folding, little is known about their physical properties or the folding transitions that allow them to continue folding productively. Folding of the Tetrahymena group I ribozyme includes sequential accumulation of two intermediates, termed I(trap) and misfolded (M). Here, we probe the structure and folding transition of I(trap) and compare them to those of M. Hydroxyl radical and dimethyl sulfate footprinting show that both I(trap) and M are extensively structured and crudely resemble the native RNA. However, regions of the core P3-P8 domain are more exposed to solvent in I(trap) than in M. I(trap) rearranges to continue folding nearly 1000-fold faster than M, and urea accelerates folding of I(trap) much less than M. Thus, the rate-limiting transition from I(trap) requires a smaller increase in exposed surface. Mutations that disrupt peripheral tertiary contacts give large and nearly uniform increases in re-folding of M, whereas the same mutations give at most modest increases in folding from I(trap). Intriguingly, mutations within the peripheral element P5abc give 5- to 10-fold accelerations in escape from I(trap), whereas ablation of P13, which lies on the opposite surface in the native structure, near the P3-P8 domain, has no effect. Thus, the unfolding required from I(trap) appears to be local, whereas the unfolding of M appears to be global. Further, the modest effects from several mutations suggest that there are multiple pathways for escape from I(trap) and that escape is aided by loosening nearby native structural constraints, presumably to facilitate local movements of nucleotides or segments that have not formed native contacts. Overall, these and prior results suggest a model in which the global architecture and peripheral interactions of the RNA are achieved relatively early in folding. Multiple folding and re-folding events occur on the predominant pathway to the native state, with increasing native core interactions and cooperativity as folding progresses.

Catalytic activity as a probe of native RNA folding

As RNAs fold to functional structures, they traverse complex energy landscapes that include many partially folded and misfolded intermediates. For structured RNAs that possess catalytic activity, this activity can provide a powerful means of monitoring folding that is complementary to biophysical approaches. RNA catalysis can be used to track accumulation of the native RNA specifically and quantitatively, readily distinguishing the native structure from intermediates that resemble it and may not be differentiated by other approaches. Here, we outline how to design and interpret experiments using catalytic activity to monitor RNA folding, and we summarize adaptations of the method that have been used to probe aspects of folding well beyond determination of the folding rates.

Protein roles in group I intron RNA folding: the tyrosyl-tRNA synthetase CYT-18 stabilizes the native state relative to a long-lived misfolded structure without compromising folding kinetics

The Neurospora crassa CYT-18 protein is a mitochondrial tyrosyl-tRNA synthetase that also promotes self-splicing of group I intron RNAs by stabilizing the functional structure in the conserved core. CYT-18 binds the core along the same surface as a common peripheral element, P5abc, suggesting that CYT-18 can replace P5abc functionally. In addition to stabilizing structure generally, P5abc stabilizes the native conformation of the Tetrahymena group I intron relative to a globally similar misfolded conformation that has only local differences within the core and is populated significantly at equilibrium by a ribozyme variant lacking P5abc (E(DeltaP5abc)). Here, we show that CYT-18 specifically promotes formation of the native group I intron core from this misfolded conformation. Catalytic activity assays demonstrate that CYT-18 shifts the equilibrium of E(DeltaP5abc) toward the native state by at least 35-fold, and binding assays suggest an even larger effect. Thus, similar to P5abc, CYT-18 preferentially recognizes the native core, despite the global similarity of the misfolded core and despite forming crudely similar complexes, as revealed by dimethyl sulfate footprinting. Interestingly, the effects of CYT-18 and P5abc on folding kinetics differ. Whereas P5abc inhibits refolding of the misfolded conformation by forming peripheral contacts that must break during refolding, CYT-18 does not display analogous inhibition, most likely because it relies to a greater extent on direct interactions with the core. Although CYT-18 does not encounter this RNA in vivo, our results suggest that it stabilizes its cognate group I introns relative to analogous misfolded intermediates. By specifically recognizing native structural features, CYT-18 may also interact with earlier folding intermediates to avoid RNA misfolding or to trap native contacts as they form. More generally, our results highlight the ability of a protein cofactor to stabilize a functional RNA structure specifically without incurring associated costs in RNA folding kinetics.

Single-Molecule Fluorescence Studies of RNA: A Decade's Progress

Over the past decade, single-molecule fluorescence studies have elucidated the structure-function relationship of RNA molecules. The real-time observation of individual RNAs by single-molecule fluorescence has unveiled the dynamic behavior of complex RNA systems in unprecedented detail, revealing the presence of transient intermediate states and their kinetic pathways. This review provides an overview of how single-molecule fluorescence has been used to explore the dynamics of RNA folding and catalysis.

Structural rearrangements linked to global folding pathways of the Azoarcus group I ribozyme

Stable RNAs must fold into specific three-dimensional structures to be biologically active, yet many RNAs form metastable structures that compete with the native state. Our previous time-resolved footprinting experiments showed that Azoarcus group I ribozyme forms its tertiary structure rapidly (tau < 30 ms) without becoming significantly trapped in kinetic intermediates. Here, we use stopped-flow fluorescence spectroscopy to probe the global folding kinetics of a ribozyme containing 2-aminopurine in the loop of P9. The modified ribozyme was catalytically active and exhibited two equilibrium folding transitions centered at 0.3 and 1.6 mM Mg2+, consistent with previous results. Stopped-flow fluorescence revealed four kinetic folding transitions with observed rate constants of 100, 34, 1, and 0.1 s-1 at 37 degrees C. From comparison with time-resolved Fe(II)-ethylenediaminetetraacetic acid footprinting of the modified ribozyme under the same conditions, these folding transitions were assigned to formation of the IC intermediate, tertiary folding and docking of the nicked P9 tetraloop, reorganization of the P3 pseudoknot, and refolding of nonnative conformers, respectively. The footprinting results show that 50-60% of the modified ribozyme folds in less than 30 ms, while the rest of the RNA population undergoes slow structural rearrangements that control the global folding rate. The results show how small perturbations to the structure of the RNA, such as a nick in P9, populate kinetic folding intermediates that are not observed in the natural ribozyme.

Reconstructing three-dimensional shape envelopes from time-resolved small-angle X-ray scattering data

Modern computing power has made it possible to reconstruct low-resolution, three-dimensional shapes from solution small-angle X-ray scattering (SAXS) data on biomolecules without a priori knowledge of the structure. In conjunction with rapid mixing techniques, SAXS has been applied to time resolve conformational changes accompanying important biological processes, such as biomolecular folding. In response to the widespread interest in SAXS reconstructions, their value in conjunction with such time-resolved data has been examined. The group I intron from Tetrahymena thermophila and its P4-P6 subdomain are ideal model systems for investigation owing to extensive previous studies, including crystal structures. The goal of this paper is to assay the quality of reconstructions from time-resolved data given the sacrifice in signal-to-noise required to obtain sharp time resolution.

DMS footprinting of structured RNAs and RNA-protein complexes

We describe a protocol in which dimethyl sulfate (DMS) modification of the base-pairing faces of unpaired adenosine and cytidine nucleotides is used for structural analysis of RNAs and RNA-protein complexes (RNPs). The protocol is optimized for RNAs of small to moderate size (< or = 500 nt). The RNA or RNP is first exposed to DMS under conditions that promote formation of the folded structure or complex, as well as 'control' conditions that do not allow folding or complex formation. The positions and extents of modification are then determined by primer extension, polyacrylamide gel electrophoresis and quantitative analysis. From changes in the extent of modification upon folding or protein binding (appearance of a 'footprint'), it is possible to detect local changes in the secondary and tertiary structure of RNA, as well as the formation of RNA-protein contacts. This protocol takes 1.5-3 d to complete, depending on the type of analysis used.

RNA misfolding and the action of chaperones

RNA folds to a myriad of three-dimensional structures and performs an equally diverse set of functions. The ability of RNA to fold and function in vivo is all the more remarkable because, in vitro, RNA has been shown to have a strong propensity to adopt misfolded, non-functional conformations. A principal factor underlying the dominance of RNA misfolding is that local RNA structure can be quite stable even in the absence of enforcing global tertiary structure. This property allows non-native structure to persist, and it also allows native structure to form and stabilize non-native contacts or non-native topology. In recent years it has become clear that one of the central reasons for the apparent disconnect between the capabilities of RNA in vivo and its in vitro folding properties is the presence of RNA chaperones, which facilitate conformational transitions of RNA and therefore mitigate the deleterious effects of RNA misfolding. Over the past two decades, it has been demonstrated that several classes of non-specific RNA binding proteins possess profound RNA chaperone activity in vitro and when overexpressed in vivo, and at least some of these proteins appear to function as chaperones in vivo. More recently, it has been shown that certain DExD/H-box proteins function as general chaperones to facilitate folding of group I and group II introns. These proteins are RNA-dependent ATPases and have RNA helicase activity, and are proposed to function by using energy from ATP binding and hydrolysis to disrupt RNA structure and/or to displace proteins from RNA-protein complexes. This review outlines experimental studies that have led to our current understanding of the range of misfolded RNA structures, the physical origins of RNA misfolding, and the functions and mechanisms of putative RNA chaperone proteins.

Kinetic redistribution of native and misfolded RNAs by a DEAD-box chaperone

DExD/H-box proteins are ubiquitously involved in RNA-mediated processes and use ATP to accelerate conformational changes in RNA. However, their mechanisms of action, and what determines which RNA species are targeted, are not well understood. Here we show that the DExD/H-box protein CYT-19, a general RNA chaperone, mediates ATP-dependent unfolding of both the native conformation and a long-lived misfolded conformation of a group I catalytic RNA with efficiencies that depend on the stabilities of the RNA species but not on specific structural features. CYT-19 then allows the RNA to refold, changing the distribution from equilibrium to kinetic control. Because misfolding is favoured kinetically, conditions that allow unfolding of the native RNA yield large increases in the population of misfolded species. Our results suggest that DExD/H-box proteins act with sufficient breadth and efficiency to allow structured RNAs to populate a wider range of conformations than would be present at equilibrium. Thus, RNAs may face selective pressure to stabilize their active conformations relative to inactive ones to avoid significant redistribution by DExD/H-box proteins. Conversely, RNAs whose functions depend on forming multiple conformations may rely on DExD/H-box proteins to increase the populations of less stable conformations, thereby increasing their overall efficiencies.

Deletion of the P5abc peripheral element accelerates early and late folding steps of the Tetrahymena group I ribozyme

The P5abc peripheral element stabilizes the Tetrahymena group I ribozyme and enhances its catalytic activity. Despite its beneficial effects on the native structure, prior studies have shown that early formation of P5abc structure during folding can slow later folding steps. Here we use a P5abc deletion variant E(deltaP5abc) to systematically probe the role of P5abc throughout tertiary folding. Time-resolved hydroxyl radical footprinting shows that E(deltaP5abc) forms its earliest stable tertiary structure on the millisecond time scale, approximately 5-fold faster than the wild-type ribozyme, and stable structure spreads throughout E(deltaP5abc) in seconds. Nevertheless, activity measurements show that the earliest detectable formation of native E(deltaP5abc) ribozyme is much slower (approximately 0.6 min(-1)), in a manner similar to that of the wild type. Also similar, only a small fraction of E(deltaP5abc) attains the native state on this time scale under standard conditions at 25 degrees C, whereas the remainder misfolds; footprinting experiments show that the misfolded conformer shares structural features with the long-lived misfolded conformer of the wild-type ribozyme. Thus, P5abc does not have a large overall effect on the rate-limiting step(s) along this pathway. However, once misfolded, E(deltaP5abc) refolds to the native state 80-fold faster than the wild-type ribozyme and is less accelerated by urea, indicating that P5abc stabilizes the misfolded structure relative to the less-ordered transition state for refolding. Together, the results suggest that, under these conditions, even the earliest tertiary folding intermediates of the wild-type ribozyme represent misfolded species and that P5abc is principally a liability during the tertiary folding process.

Kinetic characterization of a cis- and trans-acting M2+-independent DNAzyme that depends on synthetic RNaseA-like functionality — Burst-phase kinetics from the coalescence of two active DNAzyme folds

This work describes the kinetics of the DNAzyme 925-11, a combinatorially selected, M2+-independent ribophosphodiesterase that is covalently modified with both cationic amines and imidazoles. At 13 °C, cis- and trans-cleaving constructs of 925-11 demonstrate the highest rate constants reported to date for any M2+-independent nucleic acid catalyst, investigated at physiological ionic strength and pH 7.5 (0.3 min–1for self cleavage and 0.2 min–1for intermolecular cleavage). In contrast to the cis-cleaving species, single-turnover experiments with the trans-cleaving species exhibit biphasic cleavage data, suggesting the presence of two conformations of the catalyst–substrate complex. Pulse–chase experiments demonstrate that both complexes lead to substrate cleavage. Under multiple-turnover conditions, the higher rate constant appears in a burst phase that decays to a slower steady state exhibiting a rate constant of 0.0077 min–1, a value approximating that of the slow-cleaving phase seen in single-turnover experiments. Slow product release is excluded as the source of the burst phase. An integrated rate equation is derived to describe burst-phase kinetics based on the funneling of the initial population of fast-cleaving conformation into a steady-state population composed largely of the slow-cleaving conformation.Key words: RNase mimics, DNAzymes, ribozymes, kinetics, RNA cleavage.

Probing the role of a secondary structure element at the 5'- and 3'-splice sites in group I intron self-splicing: the tetrahymena L-16 ScaI ribozyme reveals a new role of the G.U pair in self-splicing

Several ribozyme constructs have been used to dissect aspects of the group I self-splicing reaction. The Tetrahymena L-21 ScaI ribozyme, the best studied of these intron analogues, catalyzes a reaction analogous to the first step of self-splicing, in which a 5'-splice site analogue (S) and guanosine (G) are converted into a 5'-exon analogue (P) and GA. This ribozyme preserves the active site but lacks a short 5'-terminal segment (called the IGS extension herein) that forms dynamic helices, called the P1 extension and P10 helix. The P1 extension forms at the 5'-splice site in the first step of self-splicing, and P10 forms at the 3'-splice site in the second step of self-splicing. To dissect the contributions from the IGS extension and the helices it forms, we have investigated the effects of each of these elements at each reaction step. These experiments were performed with the L-16 ScaI ribozyme, which retains the IGS extension, and with 5'- and 3'-splice site analogues that differ in their ability to form the helices. The presence of the IGS extension strengthens binding of P by 40-fold, even when no new base pairs are formed. This large effect was especially surprising, as binding of S is essentially unaffected for S analogues that do not form additional base pairs with the IGS extension. Analysis of a U.U pair immediately 3' to the cleavage site suggests that a previously identified deleterious effect from a dangling U residue on the L-21 ScaI ribozyme arises from a fortuitous active site interaction and has implications for RNA tertiary structure specificity. Comparisons of the affinities of 5'-splice site analogues that form only a subset of base pairs reveal that inclusion of the conserved G.U base pair at the cleavage site of group I introns destabilizes the P1 extension >100-fold relative to the stability of a helix with all Watson-Crick base pairs. Previous structural data with model duplexes and the recent intron structures suggest that this effect can be attributed to partial unstacking of the P1 extension at the G.U step. These results suggest a previously unrecognized role of the G.U wobble pair in self-splicing: breaking cooperativity in base pair formation between P1 and the P1 extensions. This effect may facilitate replacement of the P1 extension with P10 after the first chemical step of self-splicing and release of the ligated exons after the second step of self-splicing.

Probing the mechanisms of DEAD-box proteins as general RNA chaperones: the C-terminal domain of CYT-19 mediates general recognition of RNA

The DEAD-box protein CYT-19 functions in the folding of several group I introns in vivo and a diverse set of group I and group II RNAs in vitro. Recent work using the Tetrahymena group I ribozyme demonstrated that CYT-19 possesses a second RNA-binding site, distinct from the unwinding active site, which enhances unwinding activity by binding nonspecifically to the adjacent RNA structure. Here, we probe the region of CYT-19 responsible for that binding by constructing a C-terminal truncation variant that lacks 49 amino acids and terminates at a domain boundary, as defined by limited proteolysis. This truncated protein unwinds a six-base-pair duplex, formed between the oligonucleotide substrate of the Tetrahymena ribozyme and an oligonucleotide corresponding to the internal guide sequence of the ribozyme, with near-wild-type efficiency. However, the truncated protein is activated much less than the wild-type protein when the duplex is covalently linked to the ribozyme or single-stranded or double-stranded extensions. Thus, the active site for RNA unwinding remains functional in the truncated CYT-19, but the site that binds the adjacent RNA structure has been compromised. Equilibrium binding experiments confirmed that the truncated protein binds RNA less tightly than the wild-type protein. RNA binding by the compromised site is important for chaperone activity, because the truncated protein is less active in facilitating the folding of a group I intron that requires CYT-19 in vivo. The deleted region contains arginine-rich sequences, as found in other RNA-binding proteins, and may function by tethering CYT-19 to structured RNAs, so that it can efficiently disrupt exposed, non-native structural elements, allowing them to refold. Many other DExD/H-box proteins also contain arginine-rich ancillary domains, and some of these domains may function similarly as nonspecific RNA-binding elements that enhance general RNA chaperone activity.

Nonspecific binding to structured RNA and preferential unwinding of an exposed helix by the CYT-19 protein, a DEAD-box RNA chaperone

We explore the interactions of CYT-19, a DExD/H-box protein that functions in folding of group I RNAs, with a well characterized misfolded species of the Tetrahymena ribozyme. Consistent with its function, CYT-19 accelerates refolding of the misfolded RNA to its native state. Unexpectedly, CYT-19 performs another reaction much more efficiently; it unwinds the 6-bp P1 duplex formed between the ribozyme and its oligonucleotide substrate. Furthermore, CYT-19 performs this reaction 50-fold more efficiently than it unwinds the same duplex free in solution, suggesting that it forms additional interactions with the ribozyme, most likely using a distinct RNA binding site from the one responsible for unwinding. This site can apparently bind double-stranded RNA, as attachment of a simple duplex adjacent to P1 recapitulates much of the activation provided by the ribozyme. Unwinding the native P1 duplex does not accelerate refolding of the misfolded ribozyme, implying that CYT-19 can disrupt multiple contacts on the RNA, consistent with its function in folding of multiple RNAs. Further experiments showed that the P1 duplex unwinding activity is virtually the same whether the ribozyme is misfolded or native but is abrogated by formation of tertiary contacts between the P1 duplex and the body of the ribozyme. Together these results suggest a mechanism for CYT-19 and other general DExD/H-box RNA chaperones in which the proteins bind to structured RNAs and efficiently unwind loosely associated duplexes, which biases the proteins to disrupt nonnative base pairs and gives the liberated strands an opportunity to refold.

Mechanical unfolding of RNA: from hairpins to structures with internal multiloops

Mechanical unfolding of RNA structures, ranging from hairpins to ribozymes, using laser optical tweezer experiments have begun to reveal the features of the energy landscape that cannot be easily explored using conventional experiments. Upon application of constant force (f), RNA hairpins undergo cooperative transitions from folded to unfolded states whereas subdomains of ribozymes unravel one at a time. Here, we use a self-organized polymer model and Brownian dynamics simulations to probe mechanical unfolding at constant force and constant-loading rate of four RNA structures of varying complexity. For simple hairpins, such as P5GA, application of constant force or constant loading rate results in bistable cooperative transitions between folded and unfolded states without populating any intermediates. The transition state location (DeltaxFTS) changes dramatically as the loading rate is varied. At loading rates comparable to those used in laser optical tweezer experiments, the hairpin is plastic, with DeltaxFTS being midway between folded and unfolded states; whereas at high loading rates, DeltaxFTS moves close to the folded state, i.e., RNA is brittle. For the 29-nucleotide TAR RNA with the three-nucleotide bulge, unfolding occurs in a nearly two-state manner with an occasional pause in a high free energy metastable state. Forced unfolding of the 55 nucleotides of the Hepatitis IRES domain IIa, which has a distorted L-shaped structure, results in well-populated stable intermediates. The most stable force-stabilized intermediate represents straightening of the L-shaped structure. For these structures, the unfolding pathways can be predicted using the contact map of the native structures. Unfolding of a RNA motif with internal multiloop, namely, the 109-nucleotide prohead RNA that is part of the 29 DNA packaging motor, at constant value of rf occurs with three distinct rips that represent unraveling of the paired helices. The rips represent kinetic barriers to unfolding. Our work shows 1), the response of RNA to force is largely determined by the native structure; and 2), only by probing mechanical unfolding over a wide range of forces can the underlying energy landscape be fully explored.

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

Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins.

Two distinct binding modes of a protein cofactor with its target RNA

Like most cellular RNA enzymes, the bI5 group I intron requires binding by a protein cofactor to fold correctly. Here, we use single-molecule approaches to monitor the structural dynamics of the bI5 RNA in real time as it assembles with its CBP2 protein cofactor. These experiments show that CBP2 binds to the target RNA in two distinct modes with apparently opposite effects: a "non-specific" mode that forms rapidly and induces large conformational fluctuations in the RNA, and a "specific" mode that forms slowly and stabilizes the native RNA structure. The bI5 RNA folds though multiple pathways toward the native state, typically traversing dynamic intermediate states induced by non-specific binding of CBP2. These results suggest that the protein cofactor-assisted RNA folding involves sequential non-specific and specific protein-RNA interactions. The non-specific interaction potentially increases the local concentration of CBP2 and the number of conformational states accessible to the RNA, which may promote the formation of specific RNA-protein interactions.

Concordant exploration of the kinetics of RNA folding from global and local perspectives

Time-resolved small-angle X-ray scattering (SAXS) with millisecond time-resolution reveals two discrete phases of global compaction upon Mg2+-mediated folding of the Tetrahymena thermophila ribozyme. Electrostatic relaxation of the RNA occurs rapidly and dominates the first phase of compaction during which the observed radius of gyration (R(g)) decreases from 75 angstroms to 55 angstroms. A further decrease in R(g) to 45 angstroms occurs in a well-defined second phase. An analysis of mutant ribozymes shows that the latter phase depends upon the formation of long-range tertiary contacts within the P4-P6 domain of the ribozyme; disruption of the three remaining long-range contacts linking the peripheral helices has no effect on the 55-45 angstroms compaction transition. A better understanding of the role of specific tertiary contacts in compaction was obtained by concordant time-resolved hydroxyl radical (OH) analyses that report local changes in the solvent accessibility of the RNA backbone. Comparison of the global and local measures of folding shows that formation of a subset of native tertiary contacts (i.e. those defining the ribozyme core) can occur within a highly compact ensemble whose R(g) is close to that of the fully folded ribozyme. Analyses of additional ribozyme mutants and reaction conditions establish the generality of the rapid formation of a partially collapsed state with little to no detectable tertiary structure. These studies directly link global RNA compaction with formation of tertiary structure as the molecule acquires its biologically active structure, and underscore the strong dependence on salt of both local and global measures of folding kinetics.

Structural specificity conferred by a group I RNA peripheral element

Like proteins, structured RNAs must specify a native conformation that is more stable than all other possible conformations. Local structure is much more stable for RNA than for protein, so it is likely that the principal challenge for RNA is to stabilize the native structure relative to misfolded and partially folded intermediates rather than unfolded structures. Many structured RNAs contain peripheral structural elements, which surround the core elements. Although it is clear that peripheral elements stabilize structure within RNAs that contain them, it has not yet been explored whether they specifically stabilize the native states relative to alternative folds. A two-piece version of the group I intron RNA from Tetrahymena is used here to show that the peripheral element P5abc binds to the native conformation of the rest of the RNA 50,000 times more tightly than it binds to a long-lived misfolded conformation. Thus, P5abc stabilizes the native conformation by approximately 6 kcal/mol relative to this misfolded conformation. Further, activity measurements show that for the RNA lacking P5abc, the native conformation is only marginally preferred over the misfolded conformation (<0.5 kcal/mol), indicating that the peripheral structure of this RNA is required to achieve a significant thermodynamic preference for the native state. Such "structural specificity" may be a general function of RNA peripheral domains.

Single-molecule RNA folding

Single-molecule experiments significantly expand our capability to characterize complex dynamics of biological processes. This relatively new approach has contributed significantly to our understanding of the RNA folding problem. Recent single-molecule experiments, together with structural and biochemical characterizations of RNA at the ensemble level, show that RNA molecules typically fold across a highly rugged energy landscape. As a result, long-lived folding intermediates, multiple folding pathways, and heterogeneous conformational dynamics are commonly found for RNA enzymes. While initial results have suggested that stable secondary structures are partly responsible for the rugged energy landscape of RNA, a complete mechanistic understanding of the complex folding behavior has not yet been obtained. A combination of single-molecule experiments, which are well suited to analyze transient and heterogeneous dynamic behaviors, with ensemble characterizations that can provide structural information at a superior resolution will likely provide more answers.

Single-Molecule RNA Science

The development of single-molecule detection and manipulation has allowed us to monitor the behavior of individual biological molecules and molecular complexes in real time. This approach significantly expands our capability to characterize complex dynamics of biological processes, allowing transient intermediate states and parallel kinetic pathways to be directly observed. Exploring this capability to elucidate complex dynamics, recent single-molecule experiments on RNA folding and catalysis have improved our understanding of the folding energy landscape of RNA and allowed us to better dissect complex RNA catalytic reactions, including translation by the ribosome.

RNA and protein folding: common themes and variations

Visualizing the navigation of an ensemble of unfolded molecules through the bumpy energy landscape in search of the native state gives a pictorial view of biomolecular folding. This picture, when combined with concepts in polymer theory, provides a unified theory of RNA and protein folding. Just as for proteins, the major folding free energy barrier for RNA scales sublinearly with the number of nucleotides, which allows us to extract the elusive prefactor for RNA folding. Several folding scenarios can be anticipated by considering variations in the energy landscape that depend on sequence, native topology, and external conditions. RNA and protein folding mechanism can be described by the kinetic partitioning mechanism (KPM) according to which a fraction (Phi) of molecules reaches the native state directly, whereas the remaining fraction gets kinetically trapped in metastable conformations. For two-state folders Phi approximately 1. Molecular chaperones are recruited to assist protein folding whenever Phi is small. We show that the iterative annealing mechanism, introduced to describe chaperonin-mediated folding, can be generalized to understand protein-assisted RNA folding. The major differences between the folding of proteins and RNA arise in the early stages of folding. For RNA, folding can only begin after the polyelectrolyte problem is solved, whereas protein collapse requires burial of hydrophobic residues. Cross-fertilization of ideas between the two fields should lead to an understanding of how RNA and proteins solve their folding problems.