Concerted folding of a Candida ribozyme into the catalytically active structure posterior to a rapid RNA compaction

Folding of the major population of Tetrahymena intron RNA into the catalytically active structure is trapped in a slow pathway. In this report, folding of Candida albicans intron was investigated using the trans-acting Ca.L-11 ribozyme as a model. We demonstrated that both the catalytic activity (k_obs) and compact folding equilibrium of Ca.L-11 are strongly dependent on Mg²⁺ at physiological concentrations, with both showing an Mg²⁺ Hill coefficient of 3. Formation of the compact structure of Ca.L-11 is shown to occur very rapidly, on a subsecond time scale similar to that of RNase T1 cleavage. Most of the ribozyme RNA population folds into the catalytically active structure with a rate constant of 2 min^–1 at 10 mM Mg²⁺; neither slower kinetics nor obvious Mg²⁺ inhibition is observed. These results suggest that folding of the Ca.L-11 ribozyme is initiated by a rapid magnesium-dependent RNA compaction, which is followed by a slower searching for the native contacts to form the catalytically active structure without interference from the long-lived trapped states. This model thus provides an ideal system to address a range of interesting aspects of RNA folding, such as conformational searching, ion binding and the role of productive intermediates.     

Metal Ion Dependence of Cooperative Collapse Transitions in RNA

Positively charged counterions drive RNA molecules into compact configurations that lead to their biologically active structures. To understand how the valence and size of the cations influences the collapse transition in RNA, small-angle X-ray scattering was used to follow the decrease in the radius of gyration (R_g) of the Azoarcus and Tetrahymena ribozymes in different cations. Small, multivalent cations induced the collapse of both ribozymes more efficiently than did monovalent ions. Thus, the cooperativity of the collapse transition depends on the counterion charge density. Singular value decomposition of the scattering curves showed that folding of the smaller and more thermostable Azoarcus ribozyme is well described by two components, whereas collapse of the larger Tetrahymena ribozyme involves at least one intermediate. The ion-dependent persistence length, extracted from the distance distribution of the scattering vectors, shows that the Azoarcus ribozyme is less flexible at the midpoint of transition in low-charge-density ions than in high-charge-density ions. We conclude that the formation of sequence-specific tertiary interactions in the Azoarcus ribozyme overlaps with neutralization of the phosphate charge, while tertiary folding of the Tetrahymena ribozyme requires additional counterions. Thus, the stability of the RNA structure determines its sensitivity to the valence and size of the counterions.     

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 (τ < 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 Mg²⁺, 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 °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 I_C 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.     

Folding of the group I intron ribozyme from the 26S rRNA gene of Candida albicans

Preincubation of the group I intron Ca.LSU from Candida albicans at 37°C in the absence of divalent cations results in partial folding of this intron. This is indicated by increased resistance to T1 ribonuclease cleavage of many G residues in most local helices, including P4-P6, as well as the non-local helix P7, where the G binding site is located. These changes correlate with increased gel mobility and activation of catalysis by precursor RNA containing this intron after preincubation. The presence of divalent cations or spermidine during preincubation results in formation of the predicted helices, as indicated by protection of additional G residues. However, addition of these cations during preincubation of the precursor RNA alters its gel mobility and eliminates the preincubation activation of precursor RNA seen in the absence of cations. These results suggest that, in the presence of divalent cations or spermidine, Ca.LSU folds into a more ordered, stable but misfolded conformation that is less able to convert into the catalytically active form than the ribozyme preincubated without cations. These results indicate that, like the group I intron of Tetrahymena, multiple folding pathways exist for Ca.LSU. However, it appears that the role cations play in the multiple folding pathways leading to the catalytically active form may differ between folding of these two group I introns.     

An active site rearrangement within the Tetrahymena group I ribozyme releases nonproductive interactions and allows formation of catalytic interactions

Biological catalysis hinges on the precise structural integrity of an active site that binds and transforms its substrates and meeting this requirement presents a unique challenge for RNA enzymes. Functional RNAs, including ribozymes, fold into their active conformations within rugged energy landscapes that often contain misfolded conformers. Here we uncover and characterize one such “off-pathway” species within an active site after overall folding of the ribozyme is complete. The Tetrahymena group I ribozyme (E) catalyzes cleavage of an oligonucleotide substrate (S) by an exogenous guanosine (G) cofactor. We tested whether specific catalytic interactions with G are present in the preceding E•S•G and E•G ground-state complexes. We monitored interactions with G via the effects of 2′- and 3′-deoxy (–H) and −amino (–NH₂) substitutions on G binding. These and prior results reveal that G is bound in an inactive configuration within E•G, with the nucleophilic 3′-OH making a nonproductive interaction with an active site metal ion termed M_A and with the adjacent 2′-OH making no interaction. Upon S binding, a rearrangement occurs that allows both –OH groups to contact a different active site metal ion, termed M_C, to make what are likely to be their catalytic interactions. The reactive phosphoryl group on S promotes this change, presumably by repositioning the metal ions with respect to G. This conformational transition demonstrates local rearrangements within an otherwise folded RNA, underscoring RNA''s difficulty in specifying a unique conformation and highlighting Nature''s potential to use local transitions of RNA in complex function.     

Metal ion coordination to 2′ functionality of guanosine mediates substrate-guanosine coupling in group I ribozymes: implications for conserved role of metal ions and for variability in RNA folding in ribozyme catalysis

Divalent metal ions are necessary in the self splicing reaction of group I introns, and we report that metal interaction to the 2′ position of guanosine for the Azoarcus ribozyme is required for catalysis. Moreover, this metal coordination promotes the guanosine-substrate coupled binding to the ribozyme, which is another conserved feature seen across phylogenetic boundaries. Typically there is a 4-9-fold difference in binding of G to E_free versus E · S. In the Tetrahymena ribozyme’s case this substrate-guanosine communication was attributed to conformational change(s) that lead to cooperative binding of the two cofactors which is almost nonexistent at low temperatures (4 °C). In the prokaryotic Azoarcus ribozyme we also see a 4-5-fold difference in binding of the guanosine/substrate to E_free versus E · G or E · S at 10 °C that is attributed to guanosine-substrate coupling. This coupling is diminished when the metal (Mg²⁺) coordination to the 2′ is disrupted with use of 2′-amino-2′-deoxyguanosine. The coupling is restored when softer Mn²⁺ ions are added to the buffer. This evidence generalizes a model for group I ribozyme catalysis that involves metal coordination to the 2′ position of guanosine. However, we see one striking difference in that the guanosine-substrate coupling is reversed. In the Azoarcus system (10 °C) the guanosine/substrate binds 5-fold more tightly to E_free than to E · S or E · G, which is the opposite for Tetrahymena even when the later is run at 4 °C. One implication for this difference in coupling is that the Azoarcus is in a folded state well accommodated for guanosine or substrate binding. This initial binding actually causes a conformational change that retards the subsequent binding of the second cofactor, which contrasts what was found for the Tetrahymena ribozyme. These results indicate that while the role for the metal ions in the chemical catalysis is conserved across phylogenetic boundaries, there is variability in the folding pattern of the ribozyme that leads to phosphoryl transfer.     

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 k_obs ≈ 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 (I_trap). Here we use catalytic activity to probe the folding transitions from I_trap 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 I_trap 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.     

Monitoring equilibrium changes in RNA structure by 'peroxidative' and 'oxidative' hydroxyl radical footprinting

RNA molecules play an essential role in biology. In addition to transmitting genetic information, RNA can fold into unique tertiary structures fulfilling a specific biologic role as regulator, binder or catalyst. Information about tertiary contact formation is essential to understand the function of RNA molecules. Hydroxyl radicals (•OH) are unique probes of the structure of nucleic acids due to their high reactivity and small size.¹ When used as a footprinting probe, hydroxyl radicals map the solvent accessible surface of the phosphodiester backbone of DNA¹ and RNA² with as fine as single nucleotide resolution. Hydroxyl radical footprinting can be used to identify the nucleotides within an intermolecular contact surface, e.g. in DNA-protein¹ and RNA-protein complexes. Equilibrium³ and kinetic⁴ transitions can be determined by conducting hydroxyl radical footprinting as a function of a solution variable or time, respectively. A key feature of footprinting is that limited exposure to the probe (e.g., ''single-hit kinetics'') results in the uniform sampling of each nucleotide of the polymer.⁵In this video article, we use the P4-P6 domain of the Tetrahymena ribozyme to illustrate RNA sample preparation and the determination of a Mg(II)-mediated folding isotherms. We describe the use of the well known hydroxyl radical footprinting protocol that requires H₂O₂ (we call this the ''peroxidative'' protocol) and a valuable, but not widely known, alternative that uses naturally dissolved O₂ (we call this the ''oxidative'' protocol). An overview of the data reduction, transformation and analysis procedures is presented.Download video file.^{(82M, mov)}     

Time-resolved hydroxyl-radical footprinting of RNA using Fe(II)-EDTA

Chemical footprinting methods have been used extensively to probe the structures of biologically important RNAs at nucleotide resolution. One of these methods, hydroxyl-radical footprinting, has recently been employed to study the kinetics of RNA folding. Hydroxyl radicals can be generated by a number of different methods, including Fe(II)-EDTA complexes, synchrotron radiation, and peroxynitrous acid disproportionation. The latter two methods have been used for kinetic studies of RNA folding. We have taken advantage of rapid hydroxyl-radical generation by Fe(II)-EDTA-hydrogen peroxide solutions to develop a benchtop method to study folding kinetics of RNA complexes. This technique can be performed using commercially available chemicals, and can be used to accurately define RNA folding rate constants slower than 6 min(-1). Here we report the method and an example of time-resolved footprinting on the hairpin ribozyme, a small endoribonuclease and RNA ligase.     

Molecular beacons as probes of RNA unfolding under native conditions

Hybridization of fluorescent molecular beacons provides real-time detection of RNA secondary structure with high specificity. We used molecular beacons to measure folding and unfolding rates of the Tetrahymena group I ribozyme under native conditions. A molecular beacon targeted against 15 nt in the 5′ strand of the P3 helix specifically hybridized with misfolded forms of the ribozyme, without invading the native tertiary structure. The beacon associated with the misfolded ribozyme 300 times more slowly than with an unstructured oligonucleotide containing the same target sequence, suggesting that the misfolded ribozyme core remains structured in the absence of Mg²⁺. The rate of beacon hybridization under native conditions revealed a linear relationship between the free energy of unfolding and Mg²⁺ concentration. A small fraction of the RNA population unfolded very rapidly, suggesting parallel unfolding in one step or through misfolded intermediates.     

Increased Ribozyme Activity in Crowded Solutions

Noncoding RNAs must function in the crowded environment of the cell. Previous small-angle x-ray scattering experiments showed that molecular crowders stabilize the structure of the Azoarcus group I ribozyme, allowing the ribozyme to fold at low physiological Mg²⁺ concentrations. Here, we used an RNA cleavage assay to show that the PEG and Ficoll crowder molecules increased the biochemical activity of the ribozyme, whereas sucrose did not. Crowding lowered the Mg²⁺ threshold at which activity was detected and increased total RNA cleavage at high Mg²⁺ concentrations sufficient to fold the RNA in crowded or dilute solution. After correcting for solution viscosity, the observed reaction rate was proportional to the fraction of active ribozyme. We conclude that molecular crowders stabilize the native ribozyme and favor the active structure relative to compact inactive folding intermediates.     

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.     

Productive folding to the native state by a group II intron ribozyme.

Group II introns are large catalytic RNA molecules that fold into compact structures essential for the catalysis of splicing and intron mobility reactions. Despite a growing body of information on the folded state of group II introns at equilibrium, there is currently no information on the folding pathway and little information on the ionic requirements for folding. Folding isotherms were determined by hydroxyl radical footprinting for the 32 individual protections that are distributed throughout a group II intron ribozyme derived from intron ai5gamma. The isotherms span a similar range of Mg(2+) concentrations and share a similar index of cooperativity. Time-resolved hydroxyl radical footprinting studies show that all regions of the ribozyme fold slowly and with remarkable synchrony into a single catalytically active structure at a rate comparable to those of other ribozymes studied thus far. The rate constants for the formation of tertiary contacts and recovery of catalytic activity are identical within experimental error. Catalytic activity analyses in the presence of urea provide no evidence that the slow folding of the ai5gamma intron is attributable to the presence of unproductive kinetic traps along the folding pathway. Taken together, the data suggest that the rate-limiting step for folding of group II intron ai5gamma occurs early along the reaction pathway. We propose that this behavior resembles protein folding that is limited in rate by high contact order, or the need to form key tertiary interactions from partners that are located far apart in the primary or secondary structure.     

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.     

Hinge stiffness is a barrier to RNA folding

Cation-mediated RNA folding from extended to compact, biologically active conformations relies on a temporal balance of forces. The Mg^2 +-mediated folding of the Tetrahymena thermophila ribozyme is characterized by rapid nonspecific collapse followed by tertiary-contact-induced compaction. This article focuses on an autonomously folding portion of the Tetrahymena ribozyme, its P4-P6 domain, in order to probe one facet of the rapid collapse: chain flexibility. The time evolution of P4-P6 folding was followed by global and local measures as a function of Mg^2 + concentration. While all concentrations of Mg^2 + studied are sufficient to screen the charge on the helices, the rates of compaction and tertiary contact formation diverge as the concentration of Mg^2 + increases; collapse is greatly accelerated by Mg^2 +, while tertiary contact formation is not. These studies highlight the importance of chain stiffness to RNA folding; at 10 mM Mg^2 +, a stiff hinge limits the rate of P4-P6 folding. At higher magnesium concentrations, the rate-limiting step shifts from hinge bending to tertiary contact formation.     

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active “docked” conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.     

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.