Visualizing the solvent-inaccessible core of a group II intron ribozyme

Group II introns are well recognized for their remarkable catalytic capabilities, but little is known about their three-dimensional structures. In order to obtain a global view of an active enzyme, hydroxyl radical cleavage was used to define the solvent accessibility along the backbone of a ribozyme derived from group II intron ai5gamma. These studies show that a highly homogeneous ribozyme population folds into a catalytically compact structure with an extensively internalized catalytic core. In parallel, a model of the intron core was built based on known tertiary contacts. Although constructed independently of the footprinting data, the model implicates the same elements for involvement in the catalytic core of the intron.     

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.     

Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have?a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming?a unique, stable fold.     

A folding control element for tertiary collapse of a group II intron ribozyme

Ribozymes derived from the group II intron ai5gamma collapse to a compact intermediate, folding to the native state through a slow, direct pathway that is unperturbed by kinetic traps. Molecular collapse of ribozyme D135 requires high magnesium concentrations and is thought to involve a structural element in domain 1 (D1). We used nucleotide analog interference mapping, in combination with nondenaturing gel electrophoresis, to identify RNA substructures and functional groups that are essential for D135 tertiary collapse. This revealed that the most crucial atoms for compaction are located within a small section of D1 that includes the kappa and zeta elements. This small substructure controls specific collapse of the molecule and, in later steps of the folding pathway, it forms the docking site for catalytic D5. In this way, the stage is set for proper active site formation during the earliest steps of ribozyme folding.     

Domains 2 and 3 interact to form critical elements of the group II intron active site

Group II introns are self-splicing RNA molecules that also behave as mobile genetic elements. The secondary structure of group II intron RNAs is typically described as a series of six domains that project from a central wheel. Most structural and mechanistic analyses of the intron have focused on domains 1 and 5, which contain the residues essential for catalysis, and on domain 6, which contains the branch-point adenosine. Domains 2 and 3 (D2, D3) have been shown to make important contributions to intronic activity; however, information about their function is quite limited. To elucidate the role of D2 and D3 in group II ribozyme catalysis, we built a series of multi-piece ribozyme constructs based on the ai5gamma group II intron. These constructs are designed to shed light on the roles of D2 and D3 in some of the major reactions catalyzed by the intron: 5'-exon cleavage, branching, and substrate hydrolysis. Reactions with these constructs demonstrate that D3 stimulates the chemical rate constant of group II intron reactions, and that it behaves as a form of catalytic effector. However, D3 is unable to associate independently with the ribozyme core. Docking of D3 is mediated by a short duplex that is found at the base of D2. In addition to recruiting D3 into the core, the D2 stem directs the folding of the adjacent j(2/3) linker, which is among the most conserved elements in the group II intron active site. In turn, the D2 stem contributes to 5'-splice site docking and ribozyme conformational change. Nucleotide analog interference mapping suggests an interaction between the D2 stem and D3 that builds on the known theta-theta' interaction and extends it into D3. These results establish that D3 and the base of D2 are key elements of the group II intron core and they suggest a hierarchy for active-site assembly.     

Monitoring intermediate folding states of the td group I intron in vivo

Group I introns consist of two major structural domains, the P4-P6 and P3-P9 domains, which assemble through interactions with peripheral extensions to fold into an active ribozyme. To assess group I intron folding in vivo, we probed the structure of td wild-type and mutant introns using dimethyl sulfate. The results suggest that the majority of the intron population is in the native state in accordance with the current structural model, which was refined to include two novel tertiary contacts. The importance of the loop E motif in the P7.1-P7.2 extension in assisting ribozyme folding was deduced from modeling and mutational analyses. Destabilization of stem P6 results in a deficiency in tertiary structure formation in both major domains, while weakening of stem P7 only interferes with folding of the P3-P9 domain. The different impact of mutations on the tertiary structure suggests that they interfere with folding at different stages. These results provide a first insight into the structure of folding intermediates and suggest a putative order of events in a hierarchical folding pathway in vivo.     

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.     

Structural requirement for Mg2+ binding in the group I intron core

Divalent metal ions are required for splicing of group I introns, but their role in maintaining the structure of the active site is still under investigation. Ribonuclease and hydroxyl radical footprinting of a small group I intron from Azoarcus pre-tRNA(Ile) showed that tertiary interactions between helical domains are stable in a variety of cations. Only Mg(2+), however, induced a conformational change in the intron core that correlates with self-splicing activity. Three metal ion binding sites in the catalytic core were identified by Tb(III)-dependent cleavage. Two of these are near bound substrates in a three-dimensional model of the ribozyme. A third metal ion site is near an A minor motif in P3. In the pre-tRNA, Tb(3+) cleavage was redirected to the 5' and 3' splice sites, consistent with metal-dependent activation of splice site phosphodiesters. The results show that many counterions induce global folding, but organization of the group I active site is specifically linked to Mg(2+) binding at a few sites.     

RNA molecules with conserved catalytic cores but variable peripheries fold along unique energetically optimized pathways

Functional and kinetic constraints must be efficiently balanced during the folding process of all biopolymers. To understand how homologous RNA molecules with different global architectures fold into a common core structure we determined, under identical conditions, the folding mechanisms of three phylogenetically divergent group I intron ribozymes. These ribozymes share a conserved functional core defined by topologically equivalent tertiary motifs but differ in their primary sequence, size, and structural complexity. Time-resolved hydroxyl radical probing of the backbone solvent accessible surface and catalytic activity measurements integrated with structural-kinetic modeling reveal that each ribozyme adopts a unique strategy to attain the conserved functional fold. The folding rates are not dictated by the size or the overall structural complexity, but rather by the strength of the constituent tertiary motifs which, in turn, govern the structure, stability, and lifetime of the folding intermediates. A fundamental general principle of RNA folding emerges from this study: The dominant folding flux always proceeds through an optimally structured kinetic intermediate that has sufficient stability to act as a nucleating scaffold while retaining enough conformational freedom to avoid kinetic trapping. Our results also suggest a potential role of naturally selected peripheral A-minor interactions in balancing RNA structural stability with folding efficiency.     

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.     

Group II introns: highly specific endonucleases with modular structures and diverse catalytic functions

Group II introns are large catalytic RNAs with a remarkable repertoire of reactions. Here we present construct designs and protocols that were used to develop a set of kinetic frameworks for studying the structure and reaction mechanisms of group II introns and ribozymes derived from them. In addition, we discuss application of these systems to structure/function analysis of the ai5gamma group II intron.     

Stability and cooperativity of individual tertiary contacts in RNA revealed through chemical denaturation

For proteins, understanding tertiary interactions involved in local versus global unfolding has become increasingly important for understanding the nature of the native state ensemble, the mechanisms of unfolding, and the stability of both the native and intermediate states in folding. In this work we have addressed related questions with respect to RNA structure by combining chemical denaturation and hydroxyl radical footprinting methods. We have determined unfolding isotherms for each of 26 discrete sites of protection located throughout the Tetrahymena thermophila group I ribozyme. The cooperativity of folding, m-value, and the free energy, DeltaG degrees N-U, associated with formation of each tertiary contact was determined by analysis of the isotherms. The DeltaG degrees N-U values measured in this study vary from 1.7 +/- 0.2 to 7. 6 +/- 1.2 kcal mol-1. Thus, the stability of these discrete tertiary contacts vary by almost 104. In addition, an intradomain contact and three interdomain contacts show high cooperativity (m-values of 1.1 +/- 0.2 to 1.7 +/- 0.3 kcal mol-1 M-1) indicating that these contacts exhibit global cooperatively in their folding behavior. This new approach to examining RNA stability provides an exciting comparison to our understanding of protein structure and folding mechanisms.     

Metal ion binding sites in a group II intron core

Group II introns are catalytic RNA molecules that require divalent metal ions for folding, substrate binding, and chemical catalysis. Metal ion binding sites in the group II core have now been elucidated by monitoring the site-specific RNA hydrolysis patterns of bound ions such as Tb(3+) and Mg(2+). Major sites are localized near active site elements such as domain 5 and its surrounding tertiary interaction partners. Numerous sites are also observed at intron substructures that are involved in binding and potentially activating the splice sites. These results highlight the locations of specific metal ions that are likely to play a role in ribozyme catalysis.     

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.     

Studies of RNA cleavage by photolysis of N-hydroxypyridine-2(1H)-thione. A new photochemical footprinting method

N-Hydroxypyridine-2(1H)-thione (N-HPT) has been studied as a photochemical source of hydroxyl radicals for use in photoinitiated nucleic acid footprinting experiments. Steady-state photolysis of dilute aqueous solutions of N-HPT at 350 nm in the presence of a 385 nucleotide (32)P-labeled RNA, the Tetrahymena L-21 ribozyme, resulted in cleavage of the RNA at nucleotide resolution. No cleavage of the RNA occurred in the absence of light or in the absence of N-HPT. Photolysis of the analogous pyridine lacking the N-hydroxyl group did not result in detectable amounts of RNA cleavage. The addition of RNA to preirradiated solutions of N-HPT gave no apparent RNA cleavage products, suggesting that the photoproducts of N-HPT do not result in RNA modification. Cleavage of RNA, upon photolysis in the presence of N-HPT, occurred in a sequence-independent fashion with double-stranded RNA being cleaved as efficiently as single-stranded RNA. Based on these observations, we conclude that photochemically generated diffusable hydroxyl radicals are responsible for the RNA cleavage. Experiments involving the photolysis of N-HPT in the presence of the Tetrahymena ribozyme and magnesium showed a magnesium-dependent protection from RNA cleavage due to formation of a folded RNA tertiary structure. The locations and amount of protection were identical to those observed in footprinting experiments performed with other hydroxyl radical sources. The presence of N-HPT had no effect on either the rate of folding or the catalytic activity of the folded RNA, indicating that this reagent does not disrupt RNA tertiary structure or otherwise affect activity. Thus, N-HPT is established as a new reagent for use in photoinitiated RNA footprinting experiments.     

Linking the group II intron catalytic domains: tertiary contacts and structural features of domain 3

Despite its importance for group II intron catalytic activity, structural information on conserved domain 3 (D3) is extremely limited. This domain is known to specifically stimulate the chemical rate of catalysis and to function as a 'catalytic effector'. Of all the long-range tertiary contacts that have been identified within group II introns, none has included D3 residues. Furthermore, little is known about the atoms and functional groups in D3 that contribute to catalysis. Using a nucleotide analog interference mapping assay with an extended repertoire of nucleotide analogs, we have identified functional groups in D3 that are critical for ribozyme activity. These data, together with mutational analysis, suggest the formation of noncanonical base pairs within the phylogenetically conserved internal loop at the base of D3. Finally, a related nucleotide analog interference suppression study resulted in the identification of a direct tertiary interaction between D3 and catalytic domain 5, which sheds new light on D3 function in the group II intron structure and mechanism.     

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.     

Multiple monovalent ion-dependent pathways for the folding of the L-21 Tetrahymena thermophila ribozyme

Synchrotron hydroxyl radical (*OH) footprinting is a technique that monitors the local changes in solvent accessibility of the RNA backbone on milliseconds to minutes time-scales. The Mg(2+)-dependent folding of the L-21 Sca 1 Tetrahymena thermophila ribozyme has been followed using this technique at an elevated concentration of monovalent ion (200 mM NaCl) and as a function of the initial annealing conditions and substrate. Previous studies conducted at low concentrations of monovalent ion displayed sequential folding of the P4-P6 domain, the peripheral helices and the catalytic core, with each protection displaying monophasic kinetics. For ribozyme annealed in buffer containing 200 mM NaCl and folded by the addition of 10 mM MgCl(2), multiple kinetic phases are observed for *OH protections throughout the ribozyme. The independently folding P4-P6 domain is the first to fold with its protections displaying 50-90% burst phase amplitudes. That the folding of P4-P6 within the ribozyme does not display the 100% burst phase of isolated P4-P6 at 200 mM NaCl shows that interactions with the remainder of the ribozyme impede this domain's folding. In addition, *OH protections constituting each side of a tertiary contact are not coincident in some cases, consistent with the formation of transient non-native interactions. While the peripheral contacts and triple helical scaffold exhibit substantial burst phases, the slowest protection to appear is J8/7 in the catalytic core, which displays a minimal burst amplitude and whose formation is coincident with the recovery of catalytic activity. The number of kinetic phases as well as their amplitudes and rates are different when the ribozyme is annealed in low-salt buffer and folded by the concomitant addition of monovalent and divalent cations. Annealed substrate changes the partitioning of the ribozyme among the multiple folding populations. These results provide a map of the early steps in the ribozyme's folding landscape and the degree to which the preferred pathways are dependent upon the initial reaction conditions.     

Unexpected metal ion requirements specific for catalysis of the branching reaction in a group II intron

The splicing process catalyzed by group II intron ribozymes follows the same two-step pathway as nuclear pre-mRNA splicing. In vivo, the first splicing step of wild-type introns is a transesterification reaction giving rise to a branched lariat intron-3'-exon intermediate characteristic of this splicing mode. In the wild-type introns, the ribozyme core and the substrate intron-exon junctions are carried by the same precursor molecule, making it difficult to distinguish between RNA folding and catalysis under normal splicing reactions. To characterize the catalytic step of the first transesterification reaction, we studied the reversal of this reaction, reverse branching. In this reverse reaction, the excised lariat intron and the substrate 5'-exon can be preincubated and folded separately, allowing the measure of the catalytic rate of the reaction. To measure the catalytic rate of the second splicing step, purified lariat intron-3'-exon intermediate molecules were preincubated and folded prior to the addition of 5'-exon. Conditions could be found where chemistry appeared rate limiting for both catalytic steps. Study of the metal ion requirements under these conditions resulted in the unexpected finding that, for the intron studied, substitution of magnesium ions by manganese ions enhanced the rate of the first transesterification reaction by two orders of magnitude but had virtually no effect on the second transesterification reaction or the 5' splice site cleavage by hydrolysis. Finally, the catalytic rates measured under optimal conditions for both splicing steps were faster by three orders of magnitude in the branching pathway than in the hydrolytic pathway.