Pathway modulation, circular permutation and rapid RNA folding under kinetic control

The thermodynamics and folding kinetics of a circularly permuted construct of the ribozyme from Bacillus subtilis RNase P are analyzed and compared with the folding properties of the wild-type ribozyme using optical spectroscopy and catalytic activity. The folding of the wild-type ribozyme is slow due to the rearrangement of kinetically trapped species containing misfolded structures. To test whether any misfolded structure arises from interactions between the two independently folding domains of the RNase P RNA, a circular permuted form was created where one of the two phosphodiester bonds connecting these domains is broken. This construct folds approximately 15-fold faster (t1/2 approximately nine seconds) than the wild-type ribozyme at 37 degreesC. While the complete folding of both domains is kinetically indistinguishable in the wild-type ribozyme, one domain folds much faster than the other domain in the circularly permuted construct. Hence, the major kinetic trap in the folding of the wild-type RNase P RNA involves interdomain interactions. This kinetic trap is avoidable at 37 degreesC in the circularly permuted RNA. However, at temperatures below 30 degreesC or when refolding begins from an equilibrium intermediate stabilized by submillimolar concentrations of Mg2+, a subpopulation containing an interdomain misfold still forms. These results indicate that the folding pathway of this large RNA is highly malleable and can be under kinetic control.     

Perturbed folding kinetics of circularly permuted RNAs with altered topology

The folding pathway of the Tetrahymena ribozyme correlates inversely with the sequence distance between native interactions, or contact order. The rapidly folding P4-P6 domain has a low contact order, while the slowly folding P3-P7 region has a high contact order. To examine the role of topology and contact order in RNA folding, we screened for circular permutants of the ribozyme that retain catalytic activity. Permutants beginning in the P4-P6 domain fold 5 to 20 times more slowly than the wild-type ribozyme. By contrast, 50% of a permuted RNA that disjoins a non-native interaction in P3 folds tenfold faster than the wild-type ribozyme. Hence, the probability of rapidly folding to the native state depends on the topology of tertiary domains.     

Communication between RNA folding domains revealed by folding of circularly permuted ribozymes

To study the role of sequence and topology in RNA folding, we determined the kinetic folding pathways of two circularly permuted variants of the Tetrahymena group I ribozyme, using time-resolved hydroxyl radical footprinting. Circular permutation changes the distance between interacting residues in the primary sequence, without changing the native structure of the RNA. In the natural ribozyme, tertiary interactions in the P4-P6 domain form in 1 s, while interactions in the P3-P9 form in 1-3 min at 42 degrees C. Permutation of the 5' end to G111 in the P4 helix allowed the stable P4-P6 domain to fold in 200 ms at 30 degrees C, five times faster than in the wild-type RNA, while the other domains folded five times more slowly (5-8 min). By contrast, circular permutation of the 5' end to G303 in J8/7 decreased the folding rate of the P4-P6 domain. In this permuted RNA, regions joining P2, P3 and P4 were protected in 500 ms, while the P3-P9 domain was 60-80% folded within 30 s. RNase T(1) digestion and FMN photocleavage showed that circular permutation of the RNA sequence alters the initial ensemble of secondary structures, thereby changing the tertiary folding pathways. Our results show that the natural 5'-to-3' order of the structural domains in group I ribozymes optimizes structural communication between tertiary domains and promotes self-assembly of the catalytic center.     

Fast folding of a ribozyme by stabilizing core interactions: evidence for multiple folding pathways in RNA

Folding of the Tetrahymena ribozyme under physiological conditions in vitro is limited by slow conversion of long-lived intermediates to the active structure. These intermediates arise because the most stable domain of the ribozyme folds 10-50 times more rapidly than the core region containing helix P3. Native gel electrophoresis and time-resolved X-ray-dependent hydroxyl radical cleavage revealed that mutations that weaken peripheral interactions between domains accelerated folding fivefold, while a point mutation that stabilizes P3 enabled 80 % of the mutant RNA to reach the native conformation within 30 seconds at 22 degrees C. The P3 mutation increased the folding rate of the catalytic core as much as 50-fold, so that both domains of the ribozyme were formed at approximately the same rate. The results show that the ribozyme folds rapidly without significantly populating metastable intermediates when native interactions in the ribozyme core are stabilized relative to peripheral structural elements.     

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.     

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.     

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^ΔP5abc). 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^ΔP5abc 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.     

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.     

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.     

Folding mechanism of the Tetrahymena ribozyme P4-P6 domain

Synchrotron X-ray-dependent hydroxyl radical footprinting was used to probe the folding kinetics of the P4-P6 domain of the Tetrahymena group I ribozyme, which forms a stable, closely packed tertiary structure. The 160-nt domain folds independently at a similar rate (approximately 2 s(-1)) as it does in the ribozyme, when folding is measured in 10 mM sodium cacodylate and 10 mM MgCl(2). Surprisingly, tertiary interactions around a three-helix junction (P5abc) within the P4-P6 domain fold at least 25 times more rapidly (k >/= 50 s(-1)) in isolation, than when part of the wild-type P4-P6 RNA. This difference implies that long-range interactions in the P4-P6 domain can interfere with folding of P5abc. P4-P6 was observed to fold much faster at higher ionic strength than in 10 mM sodium cacodylate. Analytical centrifugation was used to measure the sedimentation and diffusion coefficients of the unfolded RNA. The hydrodynamic radius of the RNA decreased from 58 to 46 A over the range of 0-100 mM NaCl. We propose that at low ionic strength, the addition of Mg(2+) causes the domain to collapse to a compact intermediate where P5abc is trapped in a non-native structure. At high ionic strength, the RNA rapidly collapses to the native structure. Faster folding most likely results from a different average initial conformation of the RNA in higher salt conditions.     

Folding path of P5abc RNA involves direct coupling of secondary and tertiary structures

Folding mechanisms in which secondary structures are stabilized through the formation of tertiary interactions are well documented in protein folding but challenge the folding hierarchy normally assumed for RNA. However, it is increasingly clear that RNA could fold by a similar mechanism. P5abc, a small independently folding tertiary domain of the Tetrahymena thermophila group I ribozyme, is known to fold by a secondary structure rearrangement involving helix P5c. However, the extent of this rearrangement and the precise stage of folding that triggers it are unknown. We use experiments and simulations to show that the P5c helix switches to the native secondary structure late in the folding pathway and is directly coupled to the formation of tertiary interactions in the A-rich bulge. P5c mutations show that the switch in P5c is not rate-determining and suggest that non-native interactions in P5c aid folding rather than impede it. Our study illustrates that despite significant differences in the building blocks of proteins and RNA, there may be common ways in which they self-assemble.     

Transient RNA-protein interactions in RNA folding

The RNA folding trajectory features numerous off-pathway folding traps, which represent conformations that are often equally as stable as the native functional ones. Therefore, the conversion between these off-pathway structures and the native correctly folded ones is the critical step in RNA folding. This process, referred to as RNA refolding, is slow, and is represented by a transition state that has a characteristic high free energy. Because this kinetically limiting process occurs in vivo, proteins (called RNA chaperones) have evolved that facilitate the (re)folding of RNA molecules. Here, we present an overview of how proteins interact with RNA molecules in order to achieve properly folded states. In this respect, the discrimination between static and transient interactions is crucial, as different proteins have evolved a multitude of mechanisms for RNA remodeling. For RNA chaperones that act in a sequence-unspecific manner and without the use of external sources of energy, such as ATP, transient RNA-protein interactions represent the basis of the mode of action. By presenting stretches of positively charged amino acids that are positioned in defined spatial configurations, RNA chaperones enable the RNA backbone, via transient electrostatic interactions, to sample a wider conformational space that opens the route for efficient refolding reactions.     

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 major tertiary folding transitions of the Tetrahymena catalytic RNA.

The L-21 Tetrahymena ribozyme, an RNA molecule with sequence-specific endoribonuclease activity derived from a self-splicing group I intron, provides a model system for studying the RNA folding problem. A 160 nucleotide, independently folding domain of tertiary structure (the P4-P6 domain) comprises about half of the ribozyme. We now apply Fe(II)-EDTA cleavage to mutants of the ribozyme to explore the role of individual structural elements in tertiary folding of the RNA at equilibrium. Deletion of peripheral elements near the 3' end of the ribozyme destabilizes a region of the catalytic core (P3-P7) without altering the folding of the P4-P6 domain. Three different mutations within the P4-P6 domain that destabilize its folding also shift the folding of the P3-P7 region of the catalytic core to higher MgCl2 concentrations. We conclude that the role of the extended P4-P6 domain and of the 3'-terminal peripheral elements is at least in part to stabilize the catalytic core. The organization of RNA into independently folding domains of tertiary structure may be common in large RNAs, including ribosomal RNAs. Furthermore, the observation of domain-domain interactions in a catalytic RNA supports the feasibility of a primitive spliceosome without any proteins.     

Kinetic pathway for folding of the Tetrahymena ribozyme revealed by three UV-inducible crosslinks.

The kinetics of RNA folding were examined in the L-21 ribozyme, an RNA enzyme derived from the self-splicing Tetrahymena intron. Three UV-inducible crosslinks were mapped, characterized, and used as indicators for the folded state of the ribozyme. Together these data suggest that final structures are adopted first by the P4-P6 independently folding domain and only later in a region that positions the P1 helix (including the 5' splice site), a region whose folding is linked to that of a portion of the catalytic core. At intermediate times, a non-native structure forms in the region of the triple helical scaffold, which connects the major folding domains. At 30 degrees C, the unfolded ribozyme passes through these stages with a half-life of 2 min from the time magnesium cations are provided. At higher temperatures, the half-life is shortened but the order of events is unchanged. Thermal melting of the fully folded ribozyme also revealed a multi-stage process in which the steps of folding are reversed: the kinetically slowest structure is the least stable and melts first. These structures of the ribozyme also bind Mg2+ cooperatively and their relative affinity for binding seems to be a major determinant in the order of events during folding. Na+ can also substitute for Mg2+ to give rise to the same crosslinkable structures, but only at much higher concentrations. Specific binding sites for Mg2+ may make this cation particularly efficient at electrostatic stabilization during folding of these ribozyme structures.