Salt contribution to RNA tertiary structure folding stability

Accurate quantification of the ionic contribution to RNA folding stability could greatly enhance our ability to understand and predict RNA functions. Recently, motivated by the potential importance of ion correlation and fluctuation in RNA folding, we developed the tightly bound ion (TBI) model. Extensive experimental tests showed that the TBI model can lead to better treatment of multivalent ions than the Poisson-Boltzmann equation. In this study, we use the model to quantify the contribution of salt (Na⁺ and Mg²⁺) to the RNA tertiary structure folding free energy. Folding of the RNA tertiary structure often involves intermediates. We focus on the folding transition from an intermediate state to the native state, and compute the electrostatic folding free energy of the RNA. Based on systematic calculations for a variety of RNA molecules, we derive a set of formulas for the electrostatic free energy for tertiary structural folding as a function of the sequence length and compactness of the RNA and the Na⁺ and Mg²⁺ concentrations. Extensive comparisons with experimental data suggest that our model and the extracted empirical formulas are quite reliable.     

Energetics of hydrogen bond networks in RNA: hydrogen bonds surrounding G+1 and U42 are the major determinants for the tertiary structure stability of the hairpin ribozyme

The hairpin ribozyme, a small catalytic RNA consisting of two helix-loop-helix motifs, serves as a paradigm for RNA folding. In the active conformer, the ribozyme is docked into a compact structure via loop-loop interactions. The crystal structure of the docked hairpin ribozyme shows an intricate network of hydrogen bonding interactions at the docking interface, mediated by the base, sugar, and phosphate groups of U42 and G+1 [Rupert, P. B., and Ferre-D'Amare, A. R. (2001) Nature 410, 780-786]. To elucidate the determinants for tertiary structure stability in the hairpin ribozyme, we evaluated the energetic contributions of hydrogen bonds surrounding U42 and G+1 by time-resolved fluorescence resonance energy transfer using modified ribozymes that lack one or more of the individual interactions. Elimination of a single tertiary hydrogen bond consistently resulted in a net destabilization of approximately 2 kJ/mol. The results of double- and triple-mutant cycles suggest that individual hydrogen bonds surrounding G+1 or U42 act cooperatively and form extended hydrogen bond networks that stabilize the docked ribozyme. These results demonstrate that RNAs, similar to proteins, can exploit coupled hydrogen bond networks to stabilize the docking of distant structural domains.     

Stability of hairpin ribozyme tertiary structure is governed by the interdomain junction.

The equilibrium distributions of hairpin ribozyme conformational isomers have been examined by time-resolved fluorescence resonance energy transfer. Ribozymes partition between active (docked) and inactive (extended) conformers, characterized by unique interdomain distance distributions, which define differences in folding free energy. The active tertiary structure is stabilized both by specific interactions between the catalytic and the substrate-binding domains and by the structure of the intervening helical junction. Under physiological conditions, the docking equilibrium of the natural four-way junction dramatically favors the active conformer, while those of a three-way and the two-way junction used in gene therapy applications favor the inactive conformer.     

Nucleotide analog interference mapping of the hairpin ribozyme: implications for secondary and tertiary structure formation.

The hairpin ribozyme is a small, naturally occurring RNA capable of folding into a distinct three-dimensional structure and catalyzing a specific phosphodiester transfer reaction. We have adapted a high throughput screening procedure entitled nucleotide analog interference mapping (NAIM) to identify functional groups important for proper folding and catalysis of this ribozyme. A total of 18 phosphorothioate-tagged nucleotide analogs were used to determine the contribution made by individual ribose 2'-OH and purine functional groups to the hairpin ribozyme ligation reaction. Substitution with 2'-deoxy-nucleotide analogs disrupted activity at six sites within the ribozyme, and a unique interference pattern was observed at each of the 11 conserved purine nucleotides. In most cases where such information is available, the NAIM data agree with the previously reported single-site substitution results. The interference patterns are interpreted in comparison to the isolated loop A and loop B NMR structures and a model of the intact ribozyme. These data provide biochemical evidence in support of many, but not all, of the non-canonical base-pairs observed by NMR in each loop, and identify the functional groups most likely to participate in the tertiary interface between loop A and loop B. These groups include the 2'-OH groups of A10, G11, U12, C25, and A38, the exocyclic amine of G11, and the minor groove edge of A9 and A24. The data also predict non-A form sugar pucker geometry at U39 and U41. Based upon these results, a revised model for the loop A tertiary interaction with loop B is proposed. This work defines the chemical basis of purine nucleotide conservation in the hairpin ribozyme, and provides a basis for the design and interpretation of interference suppression experiments.     

Probing RNA tertiary structure: interhelical crosslinking of the hammerhead ribozyme.

Distinct structural models for the hammerhead ribozyme derived from single-crystal X-ray diffraction and fluorescence resonance energy transfer (FRET) measurements have been compared. Both models predict the same overall geometry, a wishbone shape with helices II and III nearly colinear and helix I positioned close to helix II. However, the relative orientations of helices I and II are different. To establish whether one of the models represents a kinetically active structure, a new crosslinking procedure was developed in which helices I and II of hammerhead ribozymes were disulfide-crosslinked via the 2' positions of specific sugar residues. Crosslinking residues on helices I and II that are close according to the X-ray structure did not appreciably reduce the catalytic efficiency. In contrast, crosslinking residues closely situated according to the FRET model dramatically reduced the cleavage rate by at least three orders of magnitude. These correlations between catalytic efficiencies and spatial proximities are consistent with the X-ray structure.     

Metal ion binding and the folding of the hairpin ribozyme

The hairpin ribozyme comprises two formally unpaired loops carried on two arms of a four-way helical RNA junction. Addition of divalent metal ions brings about a conformational transition into an antiparallel structure in which there is an intimate association between the loops to generate the active form of the ribozyme. In this study, we have used fluorescence resonance energy transfer to analyze the global folding of the complete ribozyme, and the simple four-way junction derived from it, over a wide concentration range of divalent and monovalent metal ions. The simple junction undergoes an ion-induced rotation into an antiparallel form. In the presence of a constant background concentration of sodium ions, the magnesium-ion-induced transition is characterized by noncooperative binding with a Hill coefficient n = 1. By contrast, the magnesium-ion-induced folding of the complete ribozyme is more complex, involving two distinct binding phases. The first phase occurs in the micromolar range, and involves the cooperative binding of at least three magnesium ions. This can also be achieved by high concentrations of sodium ions, and is therefore likely to be due to diffuse binding of cations at the junction and the interface of the loop-loop interaction. The second phase occurs in the millimolar range, and can only be induced by divalent metal ions. This transition occurs in response to the noncooperative, site-specific binding of magnesium ions. We observe a good correlation between the extent of ion-induced folding and cleavage activity.     

The tertiary structure of the hairpin ribozyme is formed through a slow conformational search

The hairpin ribozyme is a small catalytic RNA comprised of two internal loops carried on two adjacent arms of a four-way helical junction (4WJ). To achieve catalytic activity, the ribozyme folds into a compact conformation that facilitates the formation of tertiary interactions between the two loops. We have investigated the folding kinetics of the natural 4WJ form of the hairpin ribozyme, as well as a minimal construct consisting of just the two loop-containing duplexes, by means of stopped-flow fluorescence resonance energy transfer between donor and acceptor probes attached to the ends of the loop-bearing arms. Folding was initiated by the addition of Mg(2+) ions or a pseudosubstrate strand to the ribozyme, and the ensuing changes in the emission of both donor and acceptor were monitored over time. Both ribozyme constructs exhibited slow, biphasic kinetic behavior, attributed to two parallel folding pathways leading to compact, docked structures. Two distinct folding rates were observed across a range of Mg(2+) concentrations, and increasing amounts of Mg(2+) accelerated both rates. Notably, both rates were essentially independent of temperature, indicating that the corresponding activation enthalpies were negligible, in contrast to the large activation enthalpies generally observed for RNA folding processes. Instead, the slow folding was due to unfavorable entropy changes in reaching the transition state, indicating that the ribozyme tertiary structure forms through a slow conformational search. These features were observed in both forms of the ribozyme, indicating that the conformational search is confined to the two loop regions and is largely independent of the overall ribozyme architecture. Conformational search may be a general mechanism of tertiary structure formation in RNA.     

The fastest global events in RNA folding: electrostatic relaxation and tertiary collapse of the Tetrahymena ribozyme

Large RNAs can collapse into compact conformations well before the stable formation of the tertiary contacts that define their final folds. This study identifies likely physical mechanisms driving these early compaction events in RNA folding. We have employed time-resolved small-angle X-ray scattering to monitor the fastest global shape changes of the Tetrahymena ribozyme under different ionic conditions and with RNA mutations that remove long-range tertiary contacts. A partial collapse in each of the folding time-courses occurs within tens of milliseconds with either monovalent or divalent cations. Combined with comparison to predictions from structural models, this observation suggests a relaxation of the RNA to a more compact but denatured conformational ensemble in response to enhanced electrostatic screening at higher ionic concentrations. Further, the results provide evidence against counterion-correlation-mediated attraction between RNA double helices, a recently proposed model for early collapse. A previous study revealed a second 100 ms phase of collapse to a globular state. Surprisingly, we find that progression to this second early folding intermediate requires RNA sequence motifs that eventually mediate native long-range tertiary interactions, even though these regions of the RNA were observed to be solvent-accessible in previous footprinting studies under similar conditions. These results help delineate an analogy between the early conformational changes in RNA folding and the "burst phase" changes and molten globule formation in protein folding.     

The folding of the hairpin ribozyme: dependence on the loops and the junction

In its natural context, the hairpin ribozyme is constructed around a four-way helical junction. This presents the two loops that interact to form the active site on adjacent arms, requiring rotation into an antiparallel structure to bring them into proximity. In the present study we have compared the folding of this form of the ribozyme and subspecies lacking either the loops or the helical junction using fluorescence resonance energy transfer. The complete ribozyme as a four-way junction folds into an antiparallel structure by the cooperative binding of magnesium ions, requiring 20-40 microM for half-maximal extent of folding ([Mg2+]1/2) and a Hill coefficient n = 2. The isolated junction (lacking the loops) also folds into a corresponding antiparallel structure, but does so noncooperatively (n = 1) at a higher magnesium ion concentration ([Mg2+]1/2 = 3 mM). Introduction of a G + 1A mutation into loop A of the ribozyme results in a species with very similar folding to the simple junction, and complete loss of ribozyme activity. Removal of the junction from the ribozyme, replacing it either with a strand break (serving as a hinge) or a GC5 bulge, results in greatly impaired folding, with [Mg2+]1/2 > 20 mM. The results indicate that the natural form of the ribozyme undergoes ion-induced folding by the cooperative formation of an antiparallel junction and loop-loop interaction to generate the active form of the ribozyme. The four-way junction thus provides a scaffold in the natural RNA that facilitates the folding of the ribozyme into the active form.     

RNA folding and misfolding of the hammerhead ribozyme

The hammerhead ribozyme undergoes a well-defined two-stage folding process induced by the sequential binding of two magnesium ions. These probably correspond to the formation of domain 2 (0-500 microM magnesium ions) and domain 1 (1-20 mM magnesium ions), respectively. In this study we have used fluorescence resonance energy transfer (FRET) to analyze the ion-induced folding of a number of variants of the hammerhead ribozyme. We find that both A14G and G8U mutations are highly destabilizing, such that these species are essentially unfolded under all conditions. Thus they appear to be blocked in the first stage of the folding process, and using uranyl-induced photocleavage we show that the core is completely accessible to this probe under these conditions. Changes at G5 do not affect the first transition but appear to provide a blockage at the second stage of folding; this is true of changes in the sugar (removal of the 2'-hydroxyl group) and base (G5C mutation, previously studied by comparative gel electrophoresis). Arrest of folding at this intermediate stage leads to a pattern of uranyl-induced photocleavage that is changed from the wild-type, but suggests a structure less open than the A14G mutant. Specific photocleavage at G5 is found only in the wild-type sequence, suggesting that this ion-binding site is formed late in the folding process. In addition to folding that is blocked at selected stages, we have also observed misfolding. Thus the A13G mutation appears to result in the ion-induced formation of a novel tertiary structure.     

Functional equivalence of the uridine turn and the hairpin as building blocks of tertiary structure in the Neurospora VS ribozyme.

Mutational, kinetic, and chemical modification experiments show that one of the three-way helical junctions in the Neurospora VS ribozyme contains a uridine turn that is important for organizing the functional three-dimensional structure of this junction. Disruption of the uridine turn disrupts the structure of the junction and decreases the self-cleavage activity of the ribozyme; however, substitution of the uridine turn with a variety of different hairpins, thereby transforming the three-way junction into a four-way junction, maintains catalytic activity. Chemical modification structure probing reveals that both the native junction and the hairpin-containing junction support the same tertiary interactions required elsewhere in the ribozyme for catalysis. These observations show that functionally equivalent three-dimensional RNA structures can be built from different secondary structure elements.     

The kinetic mechanism of the hairpin ribozyme in vivo: influence of RNA helix stability on intracellular cleavage kinetics

The relationship between hairpin ribozyme structure, and cleavage and ligation kinetics, and equilibria has been characterized extensively under a variety of reaction conditions in vitro. We developed a quantitative assay of hairpin ribozyme cleavage activity in yeast to learn how structure-function relationships defined for RNA enzymes in vitro relate to RNA-mediated reactions in cells. Here, we report the effects of variation in the stability of an essential secondary structure element, H1, on intracellular cleavage kinetics. H1 is the base-paired helix formed between ribozyme and 3' cleavage product RNAs. H1 sequences with fewer than three base-pairs fail to support full activity in vitro or in vivo, arguing against any significant difference in the stability of short RNA helices under in vitro and intracellular conditions. Under standard conditions in vitro that include 10 mM MgCl(2), the internal equilibrium between cleavage and ligation of ribozyme-bound products favors ligation. Consequently, ribozymes with stable H1 sequences display sharply reduced self-cleavage rates, because cleavage is reversed by rapid re-ligation of bound products. In contrast, ribozymes with as many as 26 base-pairs in H1 continue to self-cleave at maximum rates in vivo. The failure of large products to inhibit cleavage could be explained if intracellular conditions promote rapid product dissociation or shift the internal equilibrium to favor cleavage. Model experiments in vitro suggest that the internal equilibrium between cleavage and ligation of bound products is likely to favor cleavage under intracellular ionic conditions.     

Topological constraints are major determinants of tRNA tertiary structure and dynamics and provide basis for tertiary folding cooperativity

Recent studies have shown that basic steric and connectivity constraints encoded at the secondary structure level are key determinants of 3D structure and dynamics in simple two-way RNA junctions. However, the role of these topological constraints in higher order RNA junctions remains poorly understood. Here, we use a specialized coarse-grained molecular dynamics model to directly probe the thermodynamic contributions of topological constraints in defining the 3D architecture and dynamics of transfer RNA (tRNA). Topological constraints alone restrict tRNA''s allowed conformational space by over an order of magnitude and strongly discriminate against formation of non-native tertiary contacts, providing a sequence independent source of folding specificity. Topological constraints also give rise to long-range correlations between the relative orientation of tRNA''s helices, which in turn provides a mechanism for encoding thermodynamic cooperativity between distinct tertiary interactions. These aspects of topological constraints make it such that only several tertiary interactions are needed to confine tRNA to its native global structure and specify functionally important 3D dynamics. We further show that topological constraints are conserved across tRNA''s different naturally occurring secondary structures. Taken together, our results emphasize the central role of secondary-structure-encoded topological constraints in defining RNA 3D structure, dynamics and folding.     

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.     

Mechanistic considerations for general acid-base catalysis by RNA: revisiting the mechanism of the hairpin ribozyme

Several small ribozymes carry out self-cleavage at a specific phosphodiester bond to yield 2',3'-cyclic phosphate and 5'-hydroxyl termini. Prior mechanistic and structural studies on the HDV ribozymes led to the proposal that the pK(a) of C75 is shifted toward neutrality, making it an effective general acid. Recent mechanistic studies on the hairpin ribozyme have led to models in which protonation of G8 is required for phosphodiester cleavage, either for general acid catalysis or for electrostatic stabilization. Inspection of recent crystal structures of the hairpin ribozyme, including a complex with a vanadate transition state mimic, suggests an alternative model involving general acid-base catalysis with G8 serving as the general base and A38 as the general acid. This model is consistent with the literature on the hairpin ribozyme, including pH-rate profiles of wild-type and mutant ribozymes and solvent isotope effects. General mechanistic considerations for RNA catalysis suggest that the penalty for having general acids and bases with pK(a)s removed from neutrality is not as severe as expected. These considerations suggest that general acid-base catalysis may be a common mechanistic strategy of RNA enzymes.     

Importance of specific nucleotides in the folding of the natural form of the hairpin ribozyme

The hairpin ribozyme in its natural context consists of two loops in RNA duplexes that are connected as arms of a four-way helical junction. Magnesium ions induce folding into the active conformation in which the two loops are in proximity. In this study, we have investigated nucleotides that are important to this folding process. We have analyzed the folding in terms of the cooperativity and apparent affinity for magnesium ions as a function of changes in base sequence and functional groups, using fluorescence resonance energy transfer. Our results suggest that the interaction between the loops is the sum of a number of component interactions. Some sequence variants such as A10U, G+1A, and C25U exhibit loss of cooperativity and reduced affinity of apparent magnesium ion binding. These variants are also very impaired in ribozyme cleavage activity. Nucleotides A10, G+1, and C25 thus appear to be essential in creating the conformational environment necessary for ion binding. The double variant G+1A/C25U exhibits a marked recovery of both folding and catalytic activity compared to either individual variant, consistent with the proposal of a triple-base interaction among A9, G+1, and C25 [Pinard, R., Lambert, D., Walter, N. G., Heckman, J. E., Major, F., and Burke, J. M. (1999) Biochemistry 38, 16035-16039]. However, substitution of A9 leads to relatively small changes in folding properties and cleavage activity, and the double variant G+1DAP/C25U (DAP is 2,6-diaminopurine), which could form an isosteric triple-base interaction, exhibits folding and cleavage activities that are both very impaired compared to those of the natural sequence. Our results indicate an important role for a Watson--Crick base pair between G+1 and C25; this may be buttressed by an interaction with A9, but the loss of this has less significant consequences for folding. 2'-Deoxyribose substitution leads to folding with reduced magnesium ion affinity in the following order: unmodified RNA > dA9 > dA10 > dC25 approximately dA10 plus dC25. The results are interpreted in terms of an interaction between the ribose ring of C25 and the ribose and base of A10, in agreement with the proposal of Ryder and Strobel [Ryder, S. P., and Strobel, S. A. (1999) J. Mol. Biol. 291, 295-311]. In general, there is a correlation between the ability to undergo ion-induced folding and the rate of ribozyme cleavage. An exception to this is provided by G8, for which substitution with uridine leads to severe impairment of cleavage but folding characteristics that are virtually unaltered from those of the natural species. This is consistent with a direct role for the nucleobase of G8 in the chemistry of cleavage.     

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

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

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.     

Solution structure of the loop B domain from the hairpin ribozyme

The hairpin ribozyme is a small catalytic RNA with a unique two-domain structure. Here we present the solution structure of the loop B domain of the hairpin ribozyme, which contains most of the catalytically essential nucleotides. The 38-nucleotide domain contains a 16-nucleotide internal loop that forms one of the largest non-Watson-Crick segments of base pairing thus far determined by either NMR or crystallography. Since the solution structure of the smaller loop A domain has been previously solved, an NMR structure-based model of the 22,000 Mr hairpin ribozyme-substrate open complex emerges by joining the two domain structures. Strikingly, catalytically essential functional groups for the loop B domain are concentrated within an expanded minor groove, presenting a clear docking surface for the loop A domain.