Structural metals in the group I intron: a ribozyme with a multiple metal ion core

Metal ions play key roles in the folding and function for many structured RNAs, including group I introns. We determined the X-ray crystal structure of the Azoarcus bacterial group I intron in complex with its 5' and 3' exons. In addition to 222 nucleotides of RNA, the model includes 18 Mg(2+) and K(+) ions. Five of the metals bind within 12 A of the scissile phosphate and coordinate the majority of the oxygen atoms biochemically implicated in conserved metal-RNA interactions. The metals are buried deep within the structure and form a multiple metal ion core that is critical to group I intron structure and function. Eight metal ions bind in other conserved regions of the intron structure, and the remaining five interact with peripheral structural elements. Each of the 18 metals mediates tertiary interactions, facilitates local bends in the sugar-phosphate backbone or binds in the major groove of helices. The group I intron has a rich history of biochemical efforts aimed to identify RNA-metal ion interactions. The structural data are correlated to the biochemical results to further understand the role of metal ions in group I intron structure and function.     

RNA tertiary interactions mediate native collapse of a bacterial group I ribozyme

Large RNAs collapse into compact intermediates in the presence of counterions before folding to the native state. We previously found that collapse of a bacterial group I ribozyme correlates with the formation of helices within the ribozyme core, but occurs at Mg2+ concentrations too low to support stable tertiary structure and catalytic activity. Here, using small-angle X-ray scattering, we show that Mg2+-induced collapse is a cooperative folding transition that can be fit by a two-state model. The Mg2+ dependence of collapse is similar to the Mg2+ dependence of helix assembly measured by partial ribonuclease T1 digestion and of an unfolding transition measured by UV hypochromicity. The correspondence between multiple probes of RNA structure further supports a two-state model. A mutation that disrupts tertiary contacts between the L9 tetraloop and its helical receptor destabilized the compact state by 0.8 kcal/mol, while mutations in the central triplex were less destabilizing. These results show that native tertiary interactions stabilize the compact folding intermediates under conditions in which the RNA backbone remains accessible to solvent.     

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

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

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.     

Functionally distinct nucleic acid binding sites for a group I intron encoded RNA maturase/DNA homing endonuclease

A large number of group I introns encode a family of homologous proteins that either promote intron splicing (maturases) or are site-specific DNA endonucleases that function in intron mobility (a process called "homing"). Genetic studies have shown that some of these proteins have both activities, yet how a single protein carries out both functions remains obscure. The similarity between respective DNA-binding sites and the RNA structure near the 5' and 3' splice sites has fueled speculation that such proteins may use analogous interactions to perform both functions. The Aspergillus nidulans mitochondrial COB group I intron encodes a bi-functional protein, I-AniI, that has both RNA maturase and site-specific DNA endonuclease activities in vitro. Here, we show that I-AniI shows distinctive features of the endonuclease family to which it belongs, including highly specific, tight binding and sequential DNA strand cleavage. Competition experiments demonstrate that I-AniI binds the COB intron RNA even in saturating concentrations of its DNA target site substrate, suggesting that the protein has a separate binding site for RNA. In addition, we provide evidence that two different DNA-binding site mutants of I-AniI have little effect on the protein's RNA maturation activity. Since RNA splicing is likely a secondary adaptation of the protein, these observations support a model in which homing endonucleases may have developed maturase function by utilizing a hitherto "non-functional" protein surface.     

Hidden Structural Modules in a Cooperative RNA Folding Transition

1. Download high-res image (325KB)
2. Download full-size image

     

Toward predicting self-splicing and protein-facilitated splicing of group I introns

In the current era of massive discoveries of noncoding RNAs within genomes, being able to infer a function from a nucleotide sequence is of paramount interest. Although studies of individual group I introns have identified self-splicing and nonself-splicing examples, there is no overall understanding of the prevalence of self-splicing or the factors that determine it among the >2300 group I introns sequenced to date. Here, the self-splicing activities of 12 group I introns from various organisms were assayed under six reaction conditions that had been shown previously to promote RNA catalysis for different RNAs. Besides revealing that assessing self-splicing under only one condition can be misleading, this survey emphasizes that in vitro self-splicing efficiency is correlated with the GC content of the intron (>35% GC was generally conductive to self-splicing), and with the ability of the introns to form particular tertiary interactions. Addition of the Neurospora crassa CYT-18 protein activated splicing of two nonself-splicing introns, but inhibited the second step of self-splicing for two others. Together, correlations between sequence, predicted structure and splicing begin to establish rules that should facilitate our ability to predict the self-splicing activity of any group I intron from its sequence.     

Analysis of Ca2+/Mg2+ selectivity in α‐lactalbumin and Ca2+‐binding lysozyme reveals a distinct Mg2+‐specific site in lysozyme

The triggering of Ca²⁺ signaling pathways relies on Ca²⁺/Mg²⁺ specificity of proteins mediating these pathways. Two homologous milk Ca²⁺‐binding proteins, bovine α‐lactalbumin (bLA) and equine lysozyme (EQL), were analyzed using the simplest “four‐state” scheme of metal‐ and temperature‐induced structural changes in a protein. The association of Ca²⁺/Mg²⁺ by native proteins is entropy‐driven. Both proteins exhibit strong temperature dependences of apparent affinities to Ca²⁺ and Mg²⁺, due to low thermal stabilities of their apo‐forms and relatively high unfavorable enthalpies of Mg²⁺ association. The ratios of their apparent affinities to Ca²⁺ and Mg²⁺, being unusually high at low temperatures (5.3–6.5 orders of magnitude), reach the values inherent to classical EF‐hand motifs at physiological temperatures. The comparison of phase diagrams predicted within the model of competitive Ca²⁺ and Mg²⁺ binding with experimental data strongly suggests that the association of Ca²⁺ and Mg²⁺ ions with bLA is a competitive process, whereas the primary Mg²⁺ site of EQL is different from its Ca²⁺‐binding site. The later conclusion is corroborated by qualitatively different molar ellipticity changes in near‐UV region accompanying Mg²⁺ and Ca²⁺ association. The Ca²⁺/Mg²⁺ selectivity of Mg²⁺‐site of EQL is below an order of magnitude. EQL exhibits a distinct Mg²⁺‐specific site, probably arising as an adaptation to the extracellular environment. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Metal ion specificities for folding and cleavage activity in the Schistosoma hammerhead ribozyme

The effects of various metal ions on cleavage activity and global folding have been studied in the extended Schistosoma hammerhead ribozyme. Fluorescence resonance energy transfer was used to probe global folding as a function of various monovalent and divalent metal ions in this ribozyme. The divalent metals ions Ca²⁺, Mg²⁺, Mn²⁺, and Sr²⁺ have a relatively small variation (less than sixfold) in their ability to globally fold the hammerhead ribozyme, which contrasts with the very large difference (>10,000-fold) in apparent rate constants for cleavage for these divalent metal ions in single-turnover kinetic experiments. There is still a very large range (>4600-fold) in the apparent rate constants for cleavage for these divalent metal ions measured in high salt (2 M NaCl) conditions where the ribozyme is globally folded. These results demonstrate that the identity of the divalent metal ion has little effect on global folding of the Schistosoma hammerhead ribozyme, whereas it has a very large effect on the cleavage kinetics. Mechanisms by which the identity of the divalent metal ion can have such a large effect on cleavage activity in the Schistosoma hammerhead ribozyme are discussed.     

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.     

Mg2+-induced variations in the conformation and dynamics of HIV-1 TAR RNA probed using NMR residual dipolar couplings

The effects of divalent Mg(2+) on the conformation and dynamics of the stem-loop transactivation response element (TAR) RNA from HIV-1 have been characterized using NMR residual dipolar couplings (RDCs). Order matrix analysis of one bond 13C-1H RDCs measured in TAR at [Mg(2+)]:[TAR] stoichiometric ratios of approximately 3:1 (TAR(3.0Mg)) and approximately 4.5:1 (TAR(4.5Mg)) revealed that Mg(2+) reduces the average inter-helical angle from 47(+/-5) degrees in TAR(free) to 5(+/-7) degrees in TAR(4.5Mg). In contrast to the TAR(free) state, the generalized degree of order for the two stems in TAR(4.5Mg) is found to be identical within experimental uncertainty, indicating that binding of Mg(2+) leads to an arrest of inter-helical motions in TAR(free). Results demonstrate that RDC-NMR methodology can provide new insight into the effects of Mg(2+) on both the conformation and dynamics of RNA.     

A kinetic and thermodynamic framework for the Azoarcus group I ribozyme reaction

Determination of quantitative thermodynamic and kinetic frameworks for ribozymes derived from the Azoarcus group I intron and comparisons to their well-studied analogs from the Tetrahymena group I intron reveal similarities and differences between these RNAs. The guanosine (G) substrate binds to the Azoarcus and Tetrahymena ribozymes with similar equilibrium binding constants and similar very slow association rate constants. These and additional literature observations support a model in which the free ribozyme is not conformationally competent to bind G and in which the probability of assuming the binding-competent state is determined by tertiary interactions of peripheral elements. As proposed previously, the slow binding of guanosine may play a role in the specificity of group I intron self-splicing, and slow binding may be used analogously in other biological processes. The internal equilibrium between ribozyme-bound substrates and products is similar for these ribozymes, but the Azoarcus ribozyme does not display the coupling in the binding of substrates that is observed with the Tetrahymena ribozyme, suggesting that local preorganization of the active site and rearrangements within the active site upon substrate binding are different for these ribozymes. Our results also confirm the much greater tertiary binding energy of the 5′-splice site analog with the Azoarcus ribozyme, binding energy that presumably compensates for the fewer base-pairing interactions to allow the 5′-exon intermediate in self splicing to remain bound subsequent to 5′-exon cleavage and prior to exon ligation. Most generally, these frameworks provide a foundation for design and interpretation of experiments investigating fundamental properties of these and other structured RNAs.     

A structural analysis of the group II intron active site and implications for the spliceosome

Group II introns are self-splicing, mobile genetic elements that have fundamentally influenced the organization of terrestrial genomes. These large ribozymes remain important for gene expression in almost all forms of bacteria and eukaryotes and they are believed to share a common ancestry with the eukaryotic spliceosome that is required for processing all nuclear pre-mRNAs. The three-dimensional structure of a group IIC intron was recently determined by X-ray crystallography, making it possible to visualize the active site and the elaborate network of tertiary interactions that stabilize the molecule. Here we describe the molecular features of the active site in detail and evaluate their correspondence with prior biochemical, genetic, and phylogenetic analyses on group II introns. In addition, we evaluate the structural significance of RNA motifs within the intron core, such as the major-groove triple helix and the domain 5 bulge. Having combined what is known about the group II intron core, we then compare it with known structural features of U6 snRNA in the eukaryotic spliceosome. This analysis leads to a set of predictions for the molecular structure of the spliceosomal active site.     

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.     

A class I ligase ribozyme with reduced Mg2+ dependence: Selection,sequence analysis,and identification of functional tertiary interactions

The class I ligase was among the first ribozymes to have been isolated from random sequences and represents the catalytic core of several RNA-directed RNA polymerase ribozymes. The ligase is also notable for its catalytic efficiency and structural complexity. Here, we report an improved version of this ribozyme, arising from selection that targeted the kinetics of the chemical step. Compared with the parent ribozyme, the improved ligase achieves a modest increase in rate enhancement under the selective conditions and shows a sharp reduction in [Mg²⁺] dependence. Analysis of the sequences and kinetics of successful clones suggests which mutations play the greatest part in these improvements. Moreover, backbone and nucleobase interference maps of the parent and improved ligase ribozymes complement the newly solved crystal structure of the improved ligase to identify the functionally significant interactions underlying the catalytic ability and structural complexity of the ligase ribozyme.     

Exploring purine N7 interactions via atomic mutagenesis: the group I ribozyme as a case study

Atomic mutagenesis has emerged as a powerful tool to unravel specific interactions in complex RNA molecules. An early extensive study of analogs of the exogenous guanosine nucleophile in group I intron self-splicing by Bass and Cech demonstrated structure-function relationships analogous to those seen for protein ligands and provided strong evidence for a well-formed substrate binding site made of RNA. Subsequent functional and structural studies have confirmed these interacting sites and extended our understanding of them, with one notable exception. Whereas 7-methyl guanosine did not affect reactivity in the original study, a subsequent study revealed a deleterious effect of the seemingly more conservative 7-deaza substitution. Here we investigate this paradox, studying these and other analogs with the more thoroughly characterized ribozyme derived from the Tetrahymena group I intron. We found that the 7-deaza substitution lowers binding by ~20-fold, relative to the cognate exogenous guanosine nucleophile, whereas binding and reaction with 7-methyl and 8-aza-7-deaza substitutions have no effect. These and additional results suggest that there is no functionally important contact between the N7 atom of the exogenous guanosine and the ribozyme. Rather, they are consistent with indirect effects introduced by the N7 substitution on stacking interactions and/or solvation that are important for binding. The set of analogs used herein should be valuable in deciphering nucleic acid interactions and how they change through reaction cycles for other RNAs and RNA/protein complexes.     

Architecture and folding mechanism of the Azoarcus Group I Pre-tRNA

Self-splicing RNAs must evolve to function in their specific exon context. The conformation of a group I pre-tRNA(ile) from the bacterium Azoarcus was probed by ribonuclease T(1) and hydroxyl radical cleavage, and by native gel electrophoresis. Biochemical data and three-dimensional models of the pre-tRNA showed that the tRNA is folded, and that the tRNA and intron sequences form separate tertiary domains. Models of the active site before steps 1 and 2 of the splicing reaction predict that exchange of the external G-cofactor and the 3'-terminal G is accomplished by a slight conformational change in P9.0 of the Azoarcus group I intron. Kinetic assays showed that the pre-tRNA folds in minutes, much more slowly than the intron alone. The dependence of the folding kinetics on Mg(2+) and the concentration of urea, and RNase T(1) experiments showed that formation of native pre-tRNA is delayed by misfolding of P3-P9, including mispairing between residues in P9 and the tRNA. Thus, although the intron and tRNA sequences form separate domains in the native pre-tRNA, their folding is coupled via metastable non-native base-pairs. This could help prevent premature processing of the 5' and 3' ends of unspliced pre-tRNA.     

The exocyclic amine at the RNase P cleavage site contributes to substrate binding and catalysis

Most tRNAs carry a G at their 5' termini, i.e. at position +1. This position corresponds to the position immediately downstream of the site of cleavage in tRNA precursors. Here we studied RNase P RNA-mediated cleavage of substrates carrying substitutions/modifications at position +1 in the absence of the RNase P protein, C5, to investigate the role of G at the RNase P cleavage site. We present data suggesting that the exocyclic amine (2NH2) of G+1 contributes to cleavage site recognition, ground state binding and catalysis by affecting the rate of cleavage. This is in contrast to O6, N7 and 2'OH that are suggested to affect ground state binding and rate of cleavage to significantly lesser extent. We also provide evidence that the effects caused by the absence of 2NH2 at position +1 influenced the charge distribution and conceivably Mg2+ binding at the RNase P cleavage site. These findings are consistent with models where the 2NH2 at the cleavage site (when present) interacts with RNase P RNA and/or influences the positioning of Mg2+ in the vicinity of the cleavage site. Moreover, our data suggest that the presence of the base at +1 is not essential for cleavage but its presence suppresses miscleavage and dramatically increases the rate of cleavage. Together our findings provide reasons why most tRNAs carry a guanosine at their 5' end.