An alternative route for the folding of large RNAs: apparent two-state folding by a group II intron ribozyme

Despite a growing literature on the folding of RNA, our understanding of tertiary folding in large RNAs derives from studies on a small set of molecular examples, with primary focus on group I introns and RNase P RNA. To broaden the scope of RNA folding models and to better understand group II intron function, we have examined the tertiary folding of a ribozyme (D135) that is derived from the self-splicing ai5gamma intron from yeast mitochondria. The D135 ribozyme folds homogeneously and cooperatively into a compact, well-defined tertiary structure that includes all regions critical for active-site organization and substrate recognition. When D135 was treated with increasing concentrations of Mg(2+) and then subjected to hydroxyl radical footprinting, similar Mg(2+) dependencies were seen for internalization of all regions of the molecule, suggesting a highly cooperative folding behavior. In this work, we show that global folding and compaction of the molecule have the same magnesium dependence as the local folding previously observed. Furthermore, urea denaturation studies indicate highly cooperative unfolding of the ribozyme that is governed by thermodynamic parameters similar to those for forward folding. In fact, D135 folds homogeneously and cooperatively from the unfolded state to its native, active structure, thereby demonstrating functional reversibility in RNA folding. Taken together, the data are consistent with two-state folding of the D135 ribozyme, which is surprising given the size and multi-domain structure of the RNA. The findings establish that the accumulation of stable intermediates prior to formation of the native state is not a universal feature of RNA folding and that there is an alternative paradigm in which the folding landscape is relatively smooth, lacking rugged features that obstruct folding to the native state.     

An obligate intermediate along the slow folding pathway of a group II intron ribozyme

Most RNA molecules collapse rapidly and reach the native state through a pathway that contains numerous traps and unproductive intermediates. The D135 group II intron ribozyme is unusual in that it can fold slowly and directly to the native state, despite its large size and structural complexity. Here we use hydroxyl radical footprinting and native gel analysis to monitor the timescale of tertiary structure collapse and to detect the presence of obligate intermediates along the folding pathway of D135. We find that structural collapse and native folding of Domain 1 precede assembly of the entire ribozyme, indicating that D1 contains an on-pathway intermediate to folding of the D135 ribozyme. Subsequent docking of Domains 3 and 5, for which D1 provides a preorganized scaffold, appears to be very fast and independent of one another. In contrast to other RNAs, the D135 ribozyme undergoes slow tertiary collapse to a compacted state, with a rate constant that is also limited by the formation D1. These findings provide a new paradigm for RNA folding and they underscore the diversity of RNA biophysical behaviors.     

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.     

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.     

Compaction of a bacterial group I ribozyme coincides with the assembly of core helices

Counterions are critical to the self-assembly of RNA tertiary structure because they neutralize the large electrostatic forces which oppose the folding process. Changes in the size and shape of the Azoarcus group I ribozyme as a function of Mg(2+) and Na(+) concentration were followed by small angle neutron scattering. In low salt buffer, the RNA was expanded, with an average radius of gyration (R(g)) of 53 +/- 1 A. A highly cooperative transition to a compact form (R(g) = 31.5 +/- 0.5 A) was observed between 1.6 and 1.7 mM MgCl(2). The collapse transition, which is unusually sharp in Mg(2+), has the characteristics of a first-order phase transition. Partial digestion with ribonuclease T1 under identical conditions showed that this transition correlated with the assembly of double helices in the ribozyme core. Fivefold higher Mg(2+) concentrations were required for self-splicing, indicating that compaction occurs before native tertiary interactions are fully stabilized. No further decrease in R(g) was observed between 1.7 and 20 mM MgCl(2), indicating that the intermediates have the same dimensions as the native ribozyme, within the uncertainty of the data (+/-1 A). A more gradual transition to a final R(g) of approximately 33.5 A was observed between 0.45 and 2 M NaCl. This confirms the expectation that monovalent ions not only are less efficient in charge neutralization but also contract the RNA less efficiently than multivalent ions.     

Small angle X-ray scattering reveals a compact intermediate in RNA folding

We have used small angle X-ray scattering (SAXS) to monitor changes in the overall size and shape of the Tetrahymena ribozyme as it folds. The native ribozyme, formed in the presence of Mg2+, is much more compact and globular than the ensemble of unfolded conformations. Time-resolved measurements show that most of the compaction occurs at least 20-fold faster than the overall folding to the native state, suggesting that a compact intermediate or family of intermediates is formed early and then rearranges in the slow steps that limit the overall folding rate. These results lead to a kinetic folding model in which an initial 'electrostatic collapse' of the RNA is followed by slower rearrangements of elements that are initially mispositioned.     

Role of counterion condensation in folding of the Tetrahymena ribozyme. II. Counterion-dependence of folding kinetics

Condensed counterions contribute to the stability of compact structures in RNA, largely by reducing electrostatic repulsion among phosphate groups. Varieties of cations induce a collapsed state in the Tetrahymena ribozyme that is readily transformed to the catalytically active structure in the presence of Mg2+. Native gel electrophoresis was used to compare the effects of the valence and size of the counterion on the kinetics of this transition. The rate of folding was found to decrease with the charge of the counterion. Transitions in monovalent ions occur 20- to 40-fold faster than transitions induced by multivalent metal ions. These results suggest that multivalent cations yield stable compact structures, which are slower to reorganize to the native conformation than those induced by monovalent ions. The folding kinetics are 12-fold faster in the presence of spermidine3+ than [Co(NH3)6]3+, consistent with less effective stabilization of long-range RNA interactions by polyamines. Under most conditions, the observed folding rate decreases with increasing counterion concentration. In saturating amounts of counterion, folding is accelerated by addition of urea. These observations indicate that reorganization of compact intermediates involves partial unfolding of the RNA. We find that folding of the ribozyme is most efficient in a mixture of monovalent salt and Mg2+. This is attributed to competition among counterions for binding to the RNA. The counterion dependence of the folding kinetics is discussed in terms of the ability of condensed ions to stabilize compact structures in RNA.     

Folding of the Tetrahymena ribozyme by polyamines: importance of counterion valence and size

Polyamines are abundant metabolites that directly influence gene expression. Although the role of polyamines in DNA condensation is well known, their role in RNA folding is less understood. Non-denaturing gel electrophoresis was used to monitor the equilibrium folding transitions of the Tetrahymena ribozyme in the presence of polyamines. All of the polyamines tested induce near-native structures that readily convert to the native conformation in Mg(2+). The stability of the folded structure increases with the charge of the polyamine and decreases with the size of the polyamine. When the counterion excluded volume becomes large, the transition to the native state does not go to completion even under favorable folding conditions. Brownian dynamics simulations of a model polyelectrolyte suggest that the kinetics of counterion-mediated collapse and the dimensions of the collapsed RNA chains depend on the structure of the counterion. The results are consistent with delocalized condensation of polyamines around the RNA. However, the effective charge of the counterions is lowered by their excluded volume. The stability of the folded RNA is enhanced when the spacing between amino groups matches the distance between adjacent phosphate groups. These results show how changes in intracellular polyamine concentrations could alter RNA folding pathways.     

Counterion charge density determines the position and plasticity of RNA folding transition states

The self-assembly of RNA structure depends on the interactions of counterions with the RNA and with each other. Comparison of various polyamines showed that the tertiary structure of the Tetrahymena ribozyme is more stable when the counterions are small and highly charged. By monitoring the folding kinetics of the ribozyme as a function of polyamine concentration, we now find that the charge density of the counterions determines the positions of the folding transition states. The transition state ensemble (TSE) between U and N moves away from the native state as the counterion valence and charge density increase, as predicted by the Hammond postulate. The TSE is broader and less structured when the RNA is refolded in polyamines rather than Mg2+. That the charge density of the counterions determines the plasticity of the TSE demonstrates the importance of interactions among condensed counterions for the self-assembly of RNA structures. We propose that the major barrier to RNA folding is dominated by entropy changes when counterion charge density is low and enthalpy differences when it is high.     

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.     

The fastest global events in RNA folding: electrostatic relaxation and tertiary collapse of the Tetrahymena ribozyme

Large RNAs can collapse into compact conformations well before the stable formation of the tertiary contacts that define their final folds. This study identifies likely physical mechanisms driving these early compaction events in RNA folding. We have employed time-resolved small-angle X-ray scattering to monitor the fastest global shape changes of the Tetrahymena ribozyme under different ionic conditions and with RNA mutations that remove long-range tertiary contacts. A partial collapse in each of the folding time-courses occurs within tens of milliseconds with either monovalent or divalent cations. Combined with comparison to predictions from structural models, this observation suggests a relaxation of the RNA to a more compact but denatured conformational ensemble in response to enhanced electrostatic screening at higher ionic concentrations. Further, the results provide evidence against counterion-correlation-mediated attraction between RNA double helices, a recently proposed model for early collapse. A previous study revealed a second 100 ms phase of collapse to a globular state. Surprisingly, we find that progression to this second early folding intermediate requires RNA sequence motifs that eventually mediate native long-range tertiary interactions, even though these regions of the RNA were observed to be solvent-accessible in previous footprinting studies under similar conditions. These results help delineate an analogy between the early conformational changes in RNA folding and the "burst phase" changes and molten globule formation in protein folding.     

A thermodynamic framework and cooperativity in the tertiary folding of a Mg2+-dependent ribozyme

The folding thermodynamics of the catalytic domain from the Bacillus subtilis RNase P RNA is analyzed using circular dichroism and fluorescence spectroscopies, hydroxyl radical protection, and catalytic activity. Folding of this 255-nucleotide ribozyme can be described with three populated species: unfolded (U), intermediate (I), and native (N) states. The U-to-I transition primarily involves secondary structure formation, whereas the I-to-N transition is dominated by tertiary structure formation. The I-to-N transition is highly cooperative as indicated by the coincidence of the four probes applied here. Two isothermal methods are used to determine the stability of the N state relative to the I state at 10 and 37 degrees C. The first method measures the extent of Mg(2+)-induced folding without urea or at constant urea concentrations. The second method measures the extent of urea-induced unfolding at constant Mg(2+) concentrations. Via application of a cooperative binding analysis, the Mg(2+) transition midpoint (K(Mg)), the Hill constant (n), and the urea-dependent surface burial parameter (m value) determined by both methods are identical, indicating that they report the same, reversible folding event. Three conclusions can be drawn from these results. (i) The folding free energy of a Mg(2+)-dependent tertiary RNA structure can be described by the K(Mg) and n parameters according to a cooperative Mg(2+) binding model. (ii) The Hill constant for this tertiary RNA structure probably represents the differential number of Mg(2+) ions bound in the I-to-N transition. (iii) Under physiological conditions, the stability of this large ribozyme is similar to that of small globular proteins.     

The tertiary structure of the hairpin ribozyme is formed through a slow conformational search

The hairpin ribozyme is a small catalytic RNA comprised of two internal loops carried on two adjacent arms of a four-way helical junction (4WJ). To achieve catalytic activity, the ribozyme folds into a compact conformation that facilitates the formation of tertiary interactions between the two loops. We have investigated the folding kinetics of the natural 4WJ form of the hairpin ribozyme, as well as a minimal construct consisting of just the two loop-containing duplexes, by means of stopped-flow fluorescence resonance energy transfer between donor and acceptor probes attached to the ends of the loop-bearing arms. Folding was initiated by the addition of Mg(2+) ions or a pseudosubstrate strand to the ribozyme, and the ensuing changes in the emission of both donor and acceptor were monitored over time. Both ribozyme constructs exhibited slow, biphasic kinetic behavior, attributed to two parallel folding pathways leading to compact, docked structures. Two distinct folding rates were observed across a range of Mg(2+) concentrations, and increasing amounts of Mg(2+) accelerated both rates. Notably, both rates were essentially independent of temperature, indicating that the corresponding activation enthalpies were negligible, in contrast to the large activation enthalpies generally observed for RNA folding processes. Instead, the slow folding was due to unfavorable entropy changes in reaching the transition state, indicating that the ribozyme tertiary structure forms through a slow conformational search. These features were observed in both forms of the ribozyme, indicating that the conformational search is confined to the two loop regions and is largely independent of the overall ribozyme architecture. Conformational search may be a general mechanism of tertiary structure formation in RNA.     

Role of counterion condensation in folding of the Tetrahymena ribozyme. I. Equilibrium stabilization by cations

Folding of RNA into an ordered, compact structure requires substantial neutralization of the negatively charged backbone by positively charged counterions. Using a native gel electrophoresis assay, we have examined the effects of counterion condensation upon the equilibrium folding of the Tetrahymena ribozyme. Incubation of the ribozyme in the presence of mono-, di- and trivalent ions induces a conformational state that is capable of rapidly forming the native structure upon brief exposure to Mg2+. The cation concentration dependence of this transition is directly correlated with the charge of the counterion used to induce folding. Substrate cleavage assays confirm the rapid onset of catalytic activity under these conditions. These results are discussed in terms of classical counterion condensation theory. A model for folding is proposed which predicts effects of charge, ionic radius and temperature on counterion-induced RNA folding transitions.     

Monovalent ion-mediated folding of the Tetrahymena thermophila ribozyme

The time-course of monovalent cation-induced folding of the L-21 Sca1 Tetrahymena thermophila ribozyme and a selected mutant was quantitatively followed using synchrotron X-ray (.OH) footprinting. Initiating folding by increasing the concentration of either Na+ or K+ to 1.5M from an initial condition of approximately 0.008 M Na+ at 42 degrees C resulted in the complete formation of tertiary contacts within the P5abc subdomain and between the peripheral helices within the dead time of our measurements (k>50 s(-1)). These results contrast with folding rates of 2-0.2 s(-1) previously observed for formation of these contacts in 10mM Mg2+ from the same initial condition. Thus, the initial formation of native tertiary contacts is inhibited by divalent but not monovalent cations. The native contacts within the catalytic core form without a detectable burst phase at rates of 0.4-1.0 s(-1) in a manner reminiscent of the Mg2+-dependent folding behavior, although tenfold faster. The tertiary interactions stabilizing the catalytic core interaction with P4-P6 and P2.1, as well as one of the protections internal for the P4-P6 domain, display progress curves with appreciable burst amplitudes and a phase comparable in rate to that of the catalytic core. That the slow folding of the ribozyme's core is a consequence of the alt-P3 secondary structure is shown by the 100% burst phase amplitudes that are observed for folding of the U273A mutant ribozyme within which the native secondary structure (P3) is strengthened. Thus, formation of a misfolded intermediate(s) resulting from the alt-P3 secondary structure is independent of ion valency while the rate at which the respective intermediates are resolved is sensitive to ion valency. The overall portrait painted by these results is that ion valency differentially affects steps in the folding process and that folding in monovalent ion alone for the U273A mutant Tetrahymena ribozyme is fast and direct.     

Principles of RNA compaction: insights from the equilibrium folding pathway of the P4-P6 RNA domain in monovalent cations

Counterions are required for RNA folding, and divalent metal ions such as Mg(2+) are often critical. To dissect the role of counterions, we have compared global and local folding of wild-type and mutant variants of P4-P6 RNA derived from the Tetrahymena group I ribozyme in monovalent and in divalent metal ions. A remarkably simple picture of the folding thermodynamics emerges. The equilibrium folding pathway in monovalent ions displays two phases. In the first phase, RNA molecules that are initially in an extended conformation enforced by charge-charge repulsion are relaxed by electrostatic screening to a state with increased flexibility but without formation of long-range tertiary contacts. At higher concentrations of monovalent ions, a state that is nearly identical to the native folded state in the presence of Mg(2+) is formed, with tertiary contacts that involve base and backbone interactions but without the subset of interactions that involve specific divalent metal ion-binding sites. The folding model derived from these and previous results provides a robust framework for understanding the equilibrium and kinetic folding of RNA.     

Perturbation of the hierarchical folding of a large RNA by the destabilization of its Scaffold's tertiary structure

The P4-P6 domain serves as a scaffold against which the periphery and catalytic core organize and fold during Mg2+-mediated folding of the Tetrahymena thermophila ribozyme. The most prominent structural motif of the P4-P6 domain is the tetraloop-tetraloop receptor interaction which "clamps" the distal parts of its hairpin-like structure. Destabilization of the tertiary structure of the P4-P6 domain by perturbation of the tetraloop-tetraloop receptor interaction alters the Mg2+-mediated folding pathway. The folding hierarchy of P5c approximately P4-P6 > periphery > catalytic core that is a striking attribute of the folding of the wild-type RNA is abolished. The initial steps in folding of the mutant RNA are > or =50-fold faster than those of the wild-type ribozyme with the earliest observed tertiary contacts forming around regions known to specifically bind Mg2+. The interaction between the mutant tetraloop and the tetraloop receptor appears coincidently with slowly forming catalytic core tertiary contacts. Thus, the stability conferred upon the P4-P6 domain by the tetraloop-tetraloop receptor interaction dictates the preferred folding pathway by stabilizing an early intermediate. A sub-denaturing concentration of urea diminishes the early barrier to folding the wild-type ribozyme along with complex effects on the subsequent steps of folding the wild-type and mutant RNA.     

An early transition state for folding of the P4-P6 RNA domain

Tertiary folding of the 160-nt P4-P6 domain of the Tetrahymena group I intron RNA involves burying of substantial surface area, providing a model for the folding of other large RNA domains involved in catalysis. Stopped-flow fluorescence was used to monitor the Mg2+-induced tertiary folding of pyrene-labeled P4-P6. At 35 degrees C with [Mg2+] approximately 10 mM, P4-P6 folds on the tens of milliseconds timescale with k(obs) = 15-31 s(-1). From these values, an activation free energy deltaG(double dagger) of approximately 8-16 kcal/mol is calculated, where the large range for deltaG(double dagger) arises from uncertainty in the pre-exponential factor relating k(obs) and delta G(double dagger). The folding rates of six mutant P4-P6 RNAs were measured and found to be similar to that of the wild-type RNA, in spite of significant thermodynamic destabilization or stabilization. The ratios of the kinetic and thermodynamic free energy changes phi = delta deltaG(double dagger)/delta deltaG(o') are approximately 0, implying a folding transition state in which most of the native-state tertiary contacts are not yet formed (an early folding transition state). The k(obs) depends on the Mg2+ concentration, and the initial slope of k(obs) versus [Mg2+] suggests that only approximately 1 Mg2+ ion is bound in the rate-limiting folding step. This is consistent with an early folding transition state, because folded P4-P6 binds many Mg2+ ions. The observation of a substantial deltaG(double dagger) despite an early folding transition state suggests that a simple two-state folding diagram for Mg2+-induced P4-P6 folding is incomplete. Our kinetic data are some of the first to provide quantitative values for an activation barrier and location of a transition state for tertiary folding of an RNA domain.