Slow folding kinetics of RNase P RNA.

Understanding the folding mechanisms of large, highly structured RNAs is important for understanding how these molecules carry out their function. Although models for the three-dimensional architecture of several large RNAs have been constructed, the process by which these structures are formed is only now beginning to be explored. The kinetic folding pathway of the Tetrahymena ribozyme involves multiple intermediates and both Mg2+-dependent and Mg2+-independent steps. To determine whether this general mechanism is representative of folding of other large RNAs, a study of RNase P RNA folding was undertaken. We show, using a kinetic oligonucleotide hybridization assay, that there is at least one slow step on the folding pathway of RNase P RNA, resulting in conformational changes in the P7 helix region on the minute timescale. Although this folding event requires the presence of Mg2+, the slow step itself does not involve Mg2+ binding. The P7 and P2 helix regions exhibit distinctly different folding behavior and ion dependence, implying that RNase P folding is likely to be a complex process. Furthermore, there are distinct similarities in the folding of RNase P RNA from both Bacillus subtilis and Escherichia coli, indicating that the folding pathway may also be conserved along with the final structure. The slow folding kinetics, Mg2+-independence of the rate, and existence of intermediates are basic features of the folding mechanism of the Tetrahymena group I intron that are also found in RNase P RNA, suggesting these may be general features of the folding of large RNAs.     

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.     

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.     

Domain structure of the ribozyme from eubacterial ribonuclease P.

Large RNAs can be composed of discrete domains that fold independently. One such "folding domain" has been identified previously in the ribozyme from Bacillus subtilis ribonuclease P (denoted P RNA). This domain contains roughly one-third of all residues. Folding of an RNA construct consisting of the remaining two-thirds of B. subtilis P RNA was examined by Fe(II)-EDTA hydroxyl radical protection. This molecule folds into the proper higher-order structure under identical conditions as the full-length P RNA, suggesting the presence of a second folding domain in B. subtilis P RNA. Folding analysis of the Escherichia coli P RNA by hydroxyl radical protection shows that this P RNA is completely folded at 5-6 mM Mg2+. In order to analyze the structural organization of folding domains in E. coli P RNA, constructs were designed based on the domain structure of B. subtilis P RNA. Fe(II)-EDTA protection indicates that E. coli P RNA also contains two folding domains. Despite the significant differences at the secondary structure level, both P RNAs appear to converge structurally at the folding domain level. The pre-tRNA substrate, localized in previous studies, may bind across the folding domains with the acceptor stem/3'CCA contacting the domain including the active site and the T stem-loop contacting the other. Because all eubacterial P RNAs share considerable homology in secondary structure to either B. subtilis or E. coli P RNA, these results suggest that this domain structure may be applicable for most, if not all, eubacterial P RNAs. Identification of folding domains should be valuable in dissecting structure-function relationship of large RNAs.     

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.     

Reduced contact order and RNA folding rates

We investigated the relationship between RNA structure and folding rates accounting for hierarchical structural formation. Folding rates of two-state folding proteins correlate well with relative contact order, a quantitative measure of the number and sequence distance between tertiary contacts. These proteins do not form stable structures prior to the rate-limiting step. In contrast, most secondary structures are stably formed prior to the rate-limiting step in RNA folding. Accordingly, we introduce "reduced contact order", a metric that reflects only the number of residues available to participate in the conformational search after the formation of secondary structure. Plotting the folding rates and the reduced contact order from ten different RNAs suggests that RNA folding can be divided into two classes. To examine this division, folding rates of circularly permutated isomers are compared for two RNAs, one from each class. Folding rates vary by tenfold for circularly permuted Bacillus subtilis RNase P RNA isomers, whereas folding rates vary by only 1.2-fold for circularly permuted catalytic domains. This difference is likely related to the dissimilar natures of their rate-limiting steps.     

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.     

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.     

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.     

Understanding the Role of Three-Dimensional Topology in Determining the?Folding Intermediates of Group I Introns

Many RNA molecules exert their biological function only after folding to unique three-dimensional structures. For long, noncoding RNA molecules, the complexity of finding the native topology can be a major impediment to correct folding to the biologically active structure. An RNA molecule may fold to a near-native structure but not be able to continue to the correct structure due to a topological barrier such as crossed strands or incorrectly stacked helices. Achieving the native conformation thus requires unfolding and refolding, resulting in a long-lived intermediate. We investigate the role of topology in the folding of two phylogenetically related catalytic group I introns, the Twort and Azoarcus group I ribozymes. The kinetic models describing the Mg²⁺-mediated folding of these ribozymes were previously determined by time-resolved hydroxyl (⋅OH) radical footprinting. Two intermediates formed by parallel intermediates were resolved for each RNA. These data and analytical ultracentrifugation compaction analyses are used herein to constrain coarse-grained models of these folding intermediates as we investigate the role of nonnative topology in dictating the lifetime of the intermediates. Starting from an ensemble of unfolded conformations, we folded the RNA molecules by progressively adding native constraints to subdomains of the RNA defined by the ⋅OH time-progress curves to simulate folding through the different kinetic pathways. We find that nonnative topologies (arrangement of helices) occur frequently in the folding simulations despite using only native constraints to drive the reaction, and that the initial conformation, rather than the folding pathway, is the major determinant of whether the RNA adopts nonnative topology during folding. From these analyses we conclude that biases in the initial conformation likely determine the relative flux through parallel RNA folding pathways.     

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.     

Assembly of an exceptionally stable RNA tertiary interface in a group I ribozyme

Group I intron RNAs contain a core of highly conserved helices flanked by peripheral domains that stabilize the core structure. In the Tetrahymena group I ribozyme, the P4, P5, and P6 helices of the core pack tightly against a three-helix subdomain called P5abc. Chemical footprinting and the crystal structure of the Tetrahymena intron P4-P6 domain revealed that tertiary interactions between these two parts of the domain create an extensive solvent-inaccessible interface. We have examined the formation and stability of this tertiary interface by providing the P5abc segment in trans to a Tetrahymena ribozyme construct that lacks P5abc (EDeltaP5abc). Equilibrium gel shift experiments show that the affinity of the P5abc and EDeltaP5abc RNAs is exceptionally strong, with a Kd of approximately 100 pM at 10 mM MgCl2 (at 37 degrees C). Chemical and enzymatic footprinting shows that the RNAs are substantially folded prior to assembly of the complex. Solvent accessibility mapping reveals that, in the absence of P5abc, the intron RNA maintains a nativelike fold but its active-site helices are not tightly packed. Upon binding of P5abc, the catalytic core becomes more tightly packed through indirect effects of the tertiary interface formation. This two-component system facilitates quantitative examination of individual tertiary contacts that stabilize the folded intron.     

Mycoplasma fermentans simplifies our view of the catalytic core of ribonuclease P RNA.

The catalytic RNA moiety of (eu)bacterial RNase P is responsible for cleavage of the 5' leader sequence from precursor tRNAs. We report the sequence, the catalytic properties, and a phylogenetic-comparative structural analysis of the RNase P RNA from Mycoplasma fermentans, at 276 nt the smallest known RNase P RNA. This RNA is noteworthy in that it lacks a stem-loop structure (helix P12) that was thought previously to be universally present in bacterial RNase P RNAs. This finding suggests that helix P12 is not required for catalytic activity in vivo. In order to test this possibility in vitro, the kinetic properties of M. fermentans RNase P RNA and a mutant Escherichia coli RNase P RNA that was engineered to lack helix P12 were determined. These RNase P RNAs are catalytically active with efficiencies (Kcat/Km) comparable to that of native E. coli RNase P RNA. These results show that helix P12 is dispensable in vivo in some organisms, and therefore is unlikely to be essential for the mechanism of RNase P action. The notion that all phylogenetically volatile structures in RNase P RNA are dispensable for the catalytic mechanism was tested. A synthetic RNA representing the phylogenetic minimum RNase P RNA was constructed by deleting all evolutionarily variable structures from the M. fermentans RNA. This simplified RNA (Micro P RNA) was catalytically active in vitro with approximately 600-fold decrease in catalytic efficiency relative to the native RNA.     

Exposing the kinetic traps in RNA folding.

Large ribozymes fold on a 'glacial' timescale compared to the folding of their protein counterparts. The sluggish folding exhibited by large RNAs results from the formation of kinetically trapped, misfolded intermediates, which are nonessential features of the folding mechanism. Newly developed mutant ribozymes that avoid kinetic traps should facilitate the study of the RNA folding problem.     

The Folding Trajectory of RNase H Is Dominated by Its Topology and Not Local Stability: A Protein Engineering Study of Variants that Fold via Two-State and Three-State Mechanisms

Proteins can sample a variety of partially folded conformations during the transition between the unfolded and native states. Some proteins never significantly populate these high-energy states and fold by an apparently two-state process. However, many proteins populate detectable, partially folded forms during the folding process. The role of such intermediates is a matter of considerable debate. A single amino acid change can convert Escherichia coli ribonuclease H from a three-state folder that populates a kinetic intermediate to one that folds in an apparent two-state fashion. We have compared the folding trajectories of the three-state RNase H and the two-state RNase H, proteins with the same native-state topology but altered regional stability, using a protein engineering approach. Our data suggest that both versions of RNase H fold through a similar trajectory with similar high-energy conformations. Mutations in the core and the periphery of the protein affect similar aspects of folding for both variants, suggesting a common trajectory with folding of the core region followed by the folding of the periphery. Our results suggest that formation of specific partially folded conformations may be a general feature of protein folding that can promote, rather than hinder, efficient folding.     

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather than Electrostatic Repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg²⁺ concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.     

Linkage of monovalent and divalent ion binding in the folding of the P4-P6 domain of the Tetrahymena ribozyme

We have explored the linkage of monovalent and divalent ion binding in the folding of the P4-P6 domain of Tetrahymena thermophila ribozyme by examining the Mg2+-induced folding and the urea-induced denaturation of the folded state as a function of Na+ under equilibrium folding conditions using hydroxyl radical footprinting. These studies allowed a thermodynamic examination of eight discrete protection sites within P4-P6 that are involved in several tertiary structure contacts. Monovalent ions compete with Mg2+ ions in mediating P4-P6 folding. The urea denaturation isotherms demonstrated DeltaDeltaG values of >2 kcal x mol(-1) in experiments conducted in 10 versus 200 mM NaCl at a constant 10 mM MgCl2. However, the individual-site isotherms reported by footprinting revealed that larger than average changes in DeltaG values were localized to specific sites within the Mg2+-rich A-bulge. The competitive effects of monovalent ions were less when K+ rather than Na+ was the monovalent cation present. This result indicates the importance of the specific K+ binding sites that are associated with AA-platform structures to P4-P6 folding and stability. These site-specific footprinting data provide quantitative and site-specific measurements of the ion-linked stability for P4-P6 that are interpreted with respect to crystallographic data.     

Mapping nucleic acid structure by hydroxyl radical cleavage

Hydroxyl radical footprinting is a widely used method for following the folding of RNA molecules in solution. This method has the unique ability to provide experimental information on the solvent accessibility of each nucleotide in an RNA molecule, so that the folding of all domains of the RNA species can be followed simultaneously at single-nucleotide resolution. In recent work, hydroxyl radical footprinting has been used, often in combination with other global measures of structure, to work out detailed folding pathways and three-dimensional structures for increasingly large and complicated RNA molecules. These include synthetic ribozymes, and group I and group II ribozymes, from yeast, the Azoarcus cyanobacterium and Tetrahymena thermophila. Advances have been made in methods for analysis of hydroxyl radical data, so that the large datasets that result from kinetic folding experiments can be analyzed in a semi-automated and quantitative manner.