Ion Condensation onto Ribozyme Is Site Specific and Fold Dependent

The highly charged RNA molecules, with each phosphate carrying a single negative charge, cannot fold into well-defined architectures with tertiary interactions in the absence of ions. For ribozymes, divalent cations are known to be more efficient than monovalent ions in driving them to a compact state, although Mg²⁺ ions are needed for catalytic activities. Therefore, how ions interact with RNA is relevant in understanding RNA folding. It is often thought that most of the ions are territorially and nonspecifically bound to the RNA, as predicted by the counterion condensation theory. Here, we show using simulations of Azoarcus ribozyme, based on an accurate coarse-grained three-site interaction model with explicit divalent and monovalent cations, that ion condensation is highly specific and depends on the nucleotide position. The regions with high coordination between the phosphate groups and the divalent cations are discernible even at very low Mg²⁺ concentrations when the ribozyme does not form tertiary interactions. Surprisingly, these regions also contain the secondary structural elements that nucleate subsequently in the self-assembly of RNA, implying that ion condensation is determined by the architecture of the folded state. These results are in sharp contrast to interactions of ions (monovalent and divalent) with rigid charged rods, in which ion condensation is uniform and position independent. The differences are explained in terms of the dramatic nonmonotonic shape fluctuations in the ribozyme as it folds with increasing Mg²⁺ or Ca²⁺ concentration.     

Molecular crowding overcomes the destabilizing effects of mutations in a bacterial ribozyme

The native structure of the Azoarcus group I ribozyme is stabilized by the cooperative formation of tertiary interactions between double helical domains. Thus, even single mutations that break this network of tertiary interactions reduce ribozyme activity in physiological Mg²⁺ concentrations. Here, we report that molecular crowding comparable to that in the cell compensates for destabilizing mutations in the Azoarcus ribozyme. Small angle X-ray scattering, native polyacrylamide gel electrophoresis and activity assays were used to compare folding free energies in dilute and crowded solutions containing 18% PEG1000. Crowder molecules allowed the wild-type and mutant ribozymes to fold at similarly low Mg²⁺ concentrations and stabilized the active structure of the mutant ribozymes under physiological conditions. This compensation helps explains why ribozyme mutations are often less deleterious in the cell than in the test tube. Nevertheless, crowding did not rescue the high fraction of folded but less active structures formed by double and triple mutants. We conclude that crowding broadens the fitness landscape by stabilizing compact RNA structures without improving the specificity of self-assembly.     

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.     

Hinge stiffness is a barrier to RNA folding

Cation-mediated RNA folding from extended to compact, biologically active conformations relies on a temporal balance of forces. The Mg^2 +-mediated folding of the Tetrahymena thermophila ribozyme is characterized by rapid nonspecific collapse followed by tertiary-contact-induced compaction. This article focuses on an autonomously folding portion of the Tetrahymena ribozyme, its P4-P6 domain, in order to probe one facet of the rapid collapse: chain flexibility. The time evolution of P4-P6 folding was followed by global and local measures as a function of Mg^2 + concentration. While all concentrations of Mg^2 + studied are sufficient to screen the charge on the helices, the rates of compaction and tertiary contact formation diverge as the concentration of Mg^2 + increases; collapse is greatly accelerated by Mg^2 +, while tertiary contact formation is not. These studies highlight the importance of chain stiffness to RNA folding; at 10 mM Mg^2 +, a stiff hinge limits the rate of P4-P6 folding. At higher magnesium concentrations, the rate-limiting step shifts from hinge bending to tertiary contact formation.     

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

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

Metal ion coordination to 2′ functionality of guanosine mediates substrate-guanosine coupling in group I ribozymes: implications for conserved role of metal ions and for variability in RNA folding in ribozyme catalysis

Divalent metal ions are necessary in the self splicing reaction of group I introns, and we report that metal interaction to the 2′ position of guanosine for the Azoarcus ribozyme is required for catalysis. Moreover, this metal coordination promotes the guanosine-substrate coupled binding to the ribozyme, which is another conserved feature seen across phylogenetic boundaries. Typically there is a 4-9-fold difference in binding of G to E_free versus E · S. In the Tetrahymena ribozyme’s case this substrate-guanosine communication was attributed to conformational change(s) that lead to cooperative binding of the two cofactors which is almost nonexistent at low temperatures (4 °C). In the prokaryotic Azoarcus ribozyme we also see a 4-5-fold difference in binding of the guanosine/substrate to E_free versus E · G or E · S at 10 °C that is attributed to guanosine-substrate coupling. This coupling is diminished when the metal (Mg²⁺) coordination to the 2′ is disrupted with use of 2′-amino-2′-deoxyguanosine. The coupling is restored when softer Mn²⁺ ions are added to the buffer. This evidence generalizes a model for group I ribozyme catalysis that involves metal coordination to the 2′ position of guanosine. However, we see one striking difference in that the guanosine-substrate coupling is reversed. In the Azoarcus system (10 °C) the guanosine/substrate binds 5-fold more tightly to E_free than to E · S or E · G, which is the opposite for Tetrahymena even when the later is run at 4 °C. One implication for this difference in coupling is that the Azoarcus is in a folded state well accommodated for guanosine or substrate binding. This initial binding actually causes a conformational change that retards the subsequent binding of the second cofactor, which contrasts what was found for the Tetrahymena ribozyme. These results indicate that while the role for the metal ions in the chemical catalysis is conserved across phylogenetic boundaries, there is variability in the folding pattern of the ribozyme that leads to phosphoryl transfer.     

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

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

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.     

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active “docked” conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.     

Toward resolution of ambiguity for the unfolded state

The unfolded states in proteins and nucleic acids remain weakly understood despite their importance in folding processes; misfolding diseases (Parkinson's and Alzheimer's); natively unfolded proteins (as many as 30% of eukaryotic proteins, according to Fink); and the study of ribozymes. Research has been hindered by the inability to quantify the residual (native) structure present in an unfolded protein or nucleic acid. Here, a scaling model is proposed to quantify the molar degree of folding and the unfolded state. The model takes a global view of protein structure and can be applied to a number of analytic methods and to simulations. Three examples are given of application to small-angle scattering from pressure-induced unfolding of SNase, from acid-unfolded cytochrome c, and from folding of Azoarcus ribozyme. These examples quantitatively show three characteristic unfolded states for proteins, the statistical nature of a protein folding pathway, and the relationship between extent of folding and chain size during folding for charge-driven folding in RNA.     

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid–base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK_a values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK_a values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid–base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid–base catalysis indicate that k₁, the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK_a. The inference that k₁ is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK_a is that of the general acid and the higher pK_a is that of the general base.     

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.     

Probing non-selective cation binding in the hairpin ribozyme with Tb(III)

Catalysis by the hairpin ribozyme is stimulated by a wide range of both simple and complex metallic and organic cations. This independence from divalent metal ion binding unequivocally excludes inner-sphere coordination to RNA as an obligatory role for metal ions in catalysis. Hence, the hairpin ribozyme is a unique model to study the role of outer-sphere coordinated cations in folding of a catalytically functional RNA structure. Here, we demonstrate that micromolar concentrations of a deprotonated aqueous complex of the lanthanide metal ion terbium(III), Tb(OH)(aq)(2+), reversibly inhibit the ribozyme by competing for a crucial, yet non-selective cation binding site. Tb(OH)(aq)(2+) also reports a likely location of this binding site through backbone hydrolysis, and permits the analysis of metal binding through sensitized luminescence. We propose that the critical cation-binding site is located at a position within the catalytic core that displays an appropriately-sized pocket and a high negative charge density. We show that cationic occupancy of this site is required for tertiary folding and catalysis, yet the site can be productively occupied by a wide variety of cations. It is striking that micromolar Tb(OH)(aq)(2+) concentrations are compatible with tertiary folding, yet interfere with catalysis. The motif implicated here in cation-binding has also been found to organize the structure of multi-helix loops in evolutionary ancient ribosomal RNAs. Our findings, therefore, illuminate general principles of non-selective outer-sphere cation binding in RNA structure and function that may have prevailed in primitive ribozymes of an early "RNA world".     

Increased Ribozyme Activity in Crowded Solutions

Noncoding RNAs must function in the crowded environment of the cell. Previous small-angle x-ray scattering experiments showed that molecular crowders stabilize the structure of the Azoarcus group I ribozyme, allowing the ribozyme to fold at low physiological Mg²⁺ concentrations. Here, we used an RNA cleavage assay to show that the PEG and Ficoll crowder molecules increased the biochemical activity of the ribozyme, whereas sucrose did not. Crowding lowered the Mg²⁺ threshold at which activity was detected and increased total RNA cleavage at high Mg²⁺ concentrations sufficient to fold the RNA in crowded or dilute solution. After correcting for solution viscosity, the observed reaction rate was proportional to the fraction of active ribozyme. We conclude that molecular crowders stabilize the native ribozyme and favor the active structure relative to compact inactive folding intermediates.     

Accommodation of Ca(II) ions for catalytic activity by a group I ribozyme

The wildtype Tetrahymena ribozyme cannot catalyze detectable levels of phosphotransfer activity in vitro on an exogenous RNA substrate oligonucleotide when calcium(II) is supplied as the only available divalent ion. Nevertheless, low-error mutants of this ribozyme have been acquired through directed evolution that do have activity in 10 mM CaCl₂. The mechanisms for such Ca(II) accommodation are not known. Here, we assayed the entire molecule in an effort to identify the roles of the mutations in allowing catalytic activity in Ca(II). We used four biochemical probing techniques - native-gel electrophoresis, hydroxyl radical footprinting, terbium(III) cleavage footprinting, and phosphorothioate interference mapping - to compare the solution structure of the wildtype ribozyme with that of a Ca(II)-active five-site mutant. We compared the gross folding patterns and specific metal-binding sites in both MgCl₂ and CaCl₂ solutions. We detected no large-scale folding differences between the two RNAs in either metal. However, we did discover a limited number of local folding differences, involving regions of the RNA affected by positions 42, 188, and 270. These data support the notion that Ca(II) is accommodated by the Tetrahymena ribozyme by a slight breathing at the active site, but that alterations at, near to, and distal from the active site can all contribute to Ca(II)-based activity.     

Metal ions and RNA folding: a highly charged topic with a dynamic future

Metal ions are required to stabilize RNA tertiary structure and to begin the folding process. How different metal ions enable RNAs to fold depends on the electrostatic potential of the RNA and correlated fluctuations in the positions of the ions themselves. Theoretical models, fluorescence spectroscopy, small angle scattering and structural biology reveal that metal ions alter the RNA dynamics and folding transition states. Specifically coordinated divalent metal ions mediate conformational rearrangements within ribozyme active sites.     

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

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

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

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

Tetrahymena thermophila and Candida albicans group I intron-derived ribozymes can catalyze the trans-excision-splicing reaction

Group I intron-derived ribozymes can catalyze a variety of non-native reactions. For the trans-excision-splicing (TES) reaction, an intron-derived ribozyme from the opportunistic pathogen Pneumocystis carinii catalyzes the excision of a predefined region from within an RNA substrate with subsequent ligation of the flanking regions. To establish TES as a general ribozyme-mediated reaction, intron-derived ribozymes from Tetrahymena thermophila and Candida albicans, which are similar to but not the same as that from Pneumocystis, were investigated for their propensity to catalyze the TES reaction. We now report that the Tetrahymena and Candida ribozymes can catalyze the excision of a single nucleotide from within their ribozyme-specific substrates. Under the conditions studied, the Tetrahymena and Candida ribozymes, however, catalyze the TES reaction with lower yields and rates [Tetrahymena (k_obs) = 0.14/min and Candida (k_obs) = 0.34/min] than the Pneumocystis ribozyme (k_obs = 3.2/min). The lower yields are likely partially due to the fact that the Tetrahymena and Candida catalyze additional reactions, separate from TES. The differences in rates are likely partially due to the individual ribozymes ability to effectively bind their 3′ terminal guanosines as intramolecular nucleophiles. Nevertheless, our results demonstrate that group I intron-derived ribozymes are inherently able to catalyze the TES reaction.