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

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.

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.

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.

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 resolution1-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 relationships4, 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.

Single-nucleotide control of tRNA folding cooperativity under near-cellular conditions

RNA folding is often studied by renaturing full-length RNA in vitro and tracking folding transitions. However, the intracellular transcript folds as it emerges from the RNA polymerase. Here, we investigate the folding pathways and stability of numerous late-transcriptional intermediates of yeast and Escherichia coli transfer RNAs (tRNAs). Transfer RNA is a highly regulated functional RNA that undergoes multiple steps of posttranscriptional processing and is found in very different lengths during its lifetime in the cell. The precursor transcript is extended on both the 5' and 3' ends of the cloverleaf core, and these extensions get trimmed before addition of the 3'-CCA and aminoacylation. We studied the thermodynamics and structures of the precursor tRNA and of late-transcriptional intermediates of the cloverleaf structure. We examined RNA folding at both the secondary and tertiary structural levels using multiple biochemical and biophysical approaches. Our findings suggest that perhaps nature has selected for a single-base addition to control folding to the functional 3D structure. In near-cellular conditions, yeast tRNA^Phe and E. coli tRNA^Ala transcripts fold in a single, cooperative transition only when nearly all of the nucleotides in the cloverleaf are transcribed by indirectly enhancing folding cooperativity. Furthermore, native extensions on the 5' and 3' ends do not interfere with cooperative core folding. This highly controlled cooperative folding has implications for recognition of tRNA by processing and modification enzymes and quality control of tRNA in cells.

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.

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.

Cooperative RNA Folding under Cellular Conditions Arises From Both Tertiary Structure Stabilization and Secondary Structure Destabilization

RNA folding has been studied extensively in vitro, typically under dilute solution conditions and abiologically high salt concentrations of 1 M Na⁺ or 10 mM Mg²⁺. The cellular environment is very different, with 20-40% crowding and only 10-40 mM Na⁺, 140 mM K⁺, and 0.5-2.0 mM Mg²⁺. As such, RNA structures and functions can be radically altered under cellular conditions. We previously reported that tRNA^phe secondary and tertiary structures unfold together in a cooperative two-state fashion under crowded in vivo-like ionic conditions, but in a noncooperative multistate fashion under dilute in vitro ionic conditions unless in nonphysiologically high concentrations of Mg²⁺. The mechanistic basis behind these effects remains unclear, however. To address the mechanism that drives RNA folding cooperativity, we probe effects of cellular conditions on structures and stabilities of individual secondary structure fragments comprising the full-length RNA. We elucidate effects of a diverse set of crowders on tRNA secondary structural fragments and full-length tRNA at three levels: at the nucleotide level by temperature-dependent in-line probing, at the tertiary structure level by small-angle X-ray scattering, and at the global level by thermal denaturation. We conclude that cooperative RNA folding is induced by two overlapping mechanisms: increased stability and compaction of tertiary structure through effects of Mg²⁺, and decreased stability of certain secondary structure elements through the effects of molecular crowders. These findings reveal that despite having very different chemical makeups RNA and protein can both have weak secondary structures in vivo leading to cooperative folding.

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.

Bridging the gap between in vitro and in vivo RNA folding

Deciphering the folding pathways and predicting the structures of complex three-dimensional biomolecules is central to elucidating biological function. RNA is single-stranded, which gives it the freedom to fold into complex secondary and tertiary structures. These structures endow RNA with the ability to perform complex chemistries and functions ranging from enzymatic activity to gene regulation. Given that RNA is involved in many essential cellular processes, it is critical to understand how it folds and functions in vivo. Within the last few years, methods have been developed to probe RNA structures in vivo and genome-wide. These studies reveal that RNA often adopts very different structures in vivo and in vitro, and provide profound insights into RNA biology. Nonetheless, both in vitro and in vivo approaches have limitations: studies in the complex and uncontrolled cellular environment make it difficult to obtain insight into RNA folding pathways and thermodynamics, and studies in vitro often lack direct cellular relevance, leaving a gap in our knowledge of RNA folding in vivo. This gap is being bridged by biophysical and mechanistic studies of RNA structure and function under conditions that mimic the cellular environment. To date, most artificial cytoplasms have used various polymers as molecular crowding agents and a series of small molecules as cosolutes. Studies under such in vivo-like conditions are yielding fresh insights, such as cooperative folding of functional RNAs and increased activity of ribozymes. These observations are accounted for in part by molecular crowding effects and interactions with other molecules. In this review, we report milestones in RNA folding in vitro and in vivo and discuss ongoing experimental and computational efforts to bridge the gap between these two conditions in order to understand how RNA folds in the cell.

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 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.