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.     

Helical junctions as determinants for RNA folding: origin of tertiary structure stability of the hairpin ribozyme

Helical junctions are ubiquitous structural elements that govern the folding and tertiary structure of RNAs. The tobacco ringspot virus hairpin ribozyme consists of two helix-loop-helix elements that lie on adjacent arms of a four-way junction. In the active form of the hairpin ribozyme, the loops are in proximity. The nature of the helical junction determines the stability of the hairpin ribozyme tertiary structure [Walter, N. G., Burke, J. M., and Millar, D. P. (1999) Nat. Struct. Biol. 6, 544-549] and thus its catalytic activity. We used two-, three-, and four-way junction hairpin ribozymes as model systems to investigate the thermodynamic basis for the different tertiary structure stabilities. The equilibrium between docked and extended conformers was analyzed as a function of temperature using time-resolved fluorescence resonance energy transfer (trFRET). As the secondary and tertiary structure transitions overlap, information from UV melting curves and trFRET had to be combined to gain insight into the thermodynamics of both structural transitions. It turned out that the higher tertiary structure stability observed in the context of a four-way junction is the result of a lower entropic cost for the docking process. In the two- and three-way junction ribozymes, a high entropic cost counteracts the favorable enthalpic term, rendering the docked conformer only marginally stable. Thus, two- and three-way junction tertiary structures are more sensitive toward regulation by ligands, whereas four-way junctions provide a stable scaffold. Altogether, RNA folding and stability appear to be governed by principles similar to those for the folding of proteins.     

Two major tertiary folding transitions of the Tetrahymena catalytic RNA.

The L-21 Tetrahymena ribozyme, an RNA molecule with sequence-specific endoribonuclease activity derived from a self-splicing group I intron, provides a model system for studying the RNA folding problem. A 160 nucleotide, independently folding domain of tertiary structure (the P4-P6 domain) comprises about half of the ribozyme. We now apply Fe(II)-EDTA cleavage to mutants of the ribozyme to explore the role of individual structural elements in tertiary folding of the RNA at equilibrium. Deletion of peripheral elements near the 3' end of the ribozyme destabilizes a region of the catalytic core (P3-P7) without altering the folding of the P4-P6 domain. Three different mutations within the P4-P6 domain that destabilize its folding also shift the folding of the P3-P7 region of the catalytic core to higher MgCl2 concentrations. We conclude that the role of the extended P4-P6 domain and of the 3'-terminal peripheral elements is at least in part to stabilize the catalytic core. The organization of RNA into independently folding domains of tertiary structure may be common in large RNAs, including ribosomal RNAs. Furthermore, the observation of domain-domain interactions in a catalytic RNA supports the feasibility of a primitive spliceosome without any proteins.     

Tertiary structure stabilization promotes hairpin ribozyme ligation.

The hairpin ribozyme catalyzes a reversible RNA cleavage reaction that participates in processing intermediates of viral satellite RNA replication in plants. A minimal hairpin ribozyme consists of two helix-loop-helix segments. These segments associate noncoaxially in the active folded structure in a way that brings catalytically important loop nucleotides into close proximity. The hairpin ribozyme in the satellite RNA of Tobacco Ringspot Virus assembles in the context of a four-way helical junction. Recent physical characterization of hairpin ribozyme structures using fluorescence resonance energy transfer demonstrated enhanced stability of the folded structure in the context of a four-way helical junction compared to minimal hairpin ribozyme variants. Analysis of the functional consequences of this modification of the helical junction has revealed two changes in the hairpin ribozyme kinetic mechanism. First, ribozymes with a four-way helical junction bind 3' cleavage products with much higher affinity than minimal hairpin ribozymes, evidence that tertiary interactions within the folded structure contribute to product binding energy. Second, the balance between ligation and cleavage shifts in favor of ligation. The enhanced ligation activity of hairpin ribozymes that contain a four-way helical junction supports the notion that tertiary structure stability is a major determinant of the hairpin ribozyme proficiency as a ligase and illustrates the link between RNA structure and biological function.     

Enhanced specificity against misfolding in a thermostable mutant of the Tetrahymena ribozyme

Structured RNAs encode native conformations that are more stable than the vast ensembles of alternative conformations, but how this specificity is evolved is incompletely understood. Here we show that a variant of the Tetrahymena group I intron ribozyme that was generated previously by in vitro selection for enhanced thermostability also displays modestly enhanced specificity against a stable misfolded structure that is globally similar to the native state, despite the absence of selective pressure to increase the energy gap between these structures. The enhanced specificity for native folding arises from mutations in two nucleotides that are close together in space in the native structure, and additional experiments show that these two mutations do not affect the stability of the misfolded conformation relative to the largely unstructured transition state ensemble for interconversion between the native and misfolded conformers. Thus, they selectively stabilize the native state, presumably by strengthening a local tertiary contact network that cannot form in the misfolded conformation. The stabilization is larger in the presence of the peripheral element P5abc, suggesting that cooperative tertiary structure formation plays a key role in the enhanced stability. The increased specificity in the absence of explicit selection suggests that the large energy gap in the wild-type RNA may have arisen analogously, a consequence of selective pressure for stability of the functional structure. More generally, the structural rigidity and intricate networks of contacts in structured RNAs may allow them to evolve substantial structural specificity without explicit negative selection, even against closely related alternative structures.     

A conserved base pair within helix P4 of the Tetrahymena ribozyme helps to form the tertiary structure required for self-splicing.

Site-specific mutagenesis of the self-splicing Tetrahymena intron has been used to investigate the function of C109-G212, a conserved base pair in the P4 stem of group I introns. Mutation of C109 to G affects splicing only slightly, whereas mutation of G212 to A or C reduces the rate of splicing substantially (500-fold reduction in kcat/Km under standard in vitro splicing conditions for the G212C mutant). Splicing activity of the compensatory double mutant (C109G:G212C) is intermediate between those of the two single mutants. Thus, the stability of the P4 stem as well as the identity of the base at position 212 are important for self-splicing. Single and double mutants containing the G212C substitution have a decreased temperature optimum for self-splicing and are partially Mg2+ suppressible, both indicative of structural destabilization. Chemical structure mapping indicates that the mutations do not redirect the global folding of the RNA, but affect the structure locally and at one other site (A183) that is distant in the secondary structure. We propose that, in addition to its pairing in P4, G212 is involved in a base triplet or an alternate base pair that contributes to the catalytically active tertiary structure of the ribozyme.     

Optimized electrostatic surfaces parallel increased thermostability: a structural bioinformatic analysis

It has been known for some time that thermophilic proteins generally have increased numbers of non-covalent interactions (salt bridges, hydrogen bonds, etc.) compared with their mesophilic orthologs. Recently, anecdotal structural comparisons suggest that non-specific acid-base ion pairs on the protein surface can be an evolutionary efficient mechanism to increase thermostability. In this comprehensive structural analysis, we confirm this to be the case. Comparison of 127 orthologous mesophilic- thermophilic protein groups indicates a clear preference for stabilizing acid-base pairs on the surface of thermophilic proteins. Compared with positions in the core, stabilizing surface mutations are less likely to disrupt the tertiary structure, and thus more likely to be evolutionarily selected. Therefore, we believe that our results, in addition to being theoretically interesting, will facilitate identification of charge-altering mutations likely to increase the stability of a particular protein structure.     

Design of hairpin ribozyme variants with improved activity for poorly processed substrates

Application of ribozymes for knockdown of RNA targets requires the identification of suitable target sites according to the consensus sequence. For the hairpin ribozyme, this was originally defined as Y?2 N?1 *G+1 U+2 Y+3 B+?, with Y = U or C, and B = U, C or G, and C being the preferred nucleobase at positions -2 and +4. In the context of development of ribozymes for destruction of an oncogenic mRNA, we have designed ribozyme variants that efficiently process RNA substrates at U?2 G?1 *G+1 U+2 A+3 A+? sites. Substrates with G?1 *G+1 U+2 A+3 sites were previously shown to be processed by the wild-type hairpin ribozyme. However, our study demonstrates that, in the specific sequence context of the substrate studied herein, compensatory base changes in the ribozyme improve activity for cleavage (eight-fold) and ligation (100-fold). In particular, we show that A+3 and A+? are well tolerated if compensatory mutations are made at positions 6 and 7 of the ribozyme strand. Adenine at position +4 is neutralized by G? →U, owing to restoration of a Watson-Crick base pair in helix 1. In this ribozyme-substrate complex, adenine at position +3 is also tolerated, with a slightly decreased cleavage rate. Additional substitution of A? with uracil doubled the cleavage rate and restored ligation, which was lost in variants with A?, C? and G?. The ability to cleave, in conjunction with the inability to ligate RNA, makes these ribozyme variants particularly suitable candidates for RNA destruction.     

Selection of a novel class of RNA-RNA interaction motifs based on the ligase ribozyme with defined modular architecture

To develop molecular tools for the detection and control of RNA molecules whose functions rely on their 3D structures, we have devised a selection system to isolate novel RNA motifs that interact with a target RNA structure within a given structural context. In this system, a GAAA tetraloop and its specific receptor motif (11-ntR) from an artificial RNA ligase ribozyme with modular architecture (the DSL ribozyme) were replaced with a target structure and random sequence, respectively. Motifs recognizing the target structure can be identified by in vitro selection based on ribozyme activity. A model selection targeting GAAA-loop successfully identified motifs previously known as GAAA-loop receptors. In addition, a new selection targeting a C-loop motif also generated novel motifs that interact with this structure. Biochemical analysis of one of the C-loop receptor motifs revealed that it could also function as an independent structural unit.     

Evolution of Tetrahymena ribozyme mutants with increased structural stability

Determining how large RNA molecules stabilize their tertiary structures is critical for understanding how they perform their biological functions. Here we use in vitro selection to identify active variants of the Tetrahymena ribozyme with increased stability. The mutant pool converged to a single family that shared nine mutations; an RNA representing the consensus sequence was structurally more stable by 10.5 degrees C and catalytically active at elevated temperatures. Remarkably, of the nine altered sites, most are already known to be involved in tertiary interactions, and the stabilizing mutations primarily improve the packing interactions in the molecular interior. The wild type ribozyme and the selected mutants provide pairs of mesophilic and thermophilic homologs for studying the origin of their thermal stability.     

Molecular dynamics investigations of hammerhead ribozyme RNA

The hammerhead ribozyme, a small catalytic RNA molecule, cleaves, in the presence of magnesium ions, a specific phosphodiester bond within its own backbone, leading to 23-cyclic phosphate and 5-OH extremities. In order to study the dynamical flexibility of the hammerhead RNA, we performed molecular dynamics simulations of the solvated crystal structure of an active hammerhead ribozyme, obtained after flash-freezing crystals soaked with magnesium. Because of a careful equilibration protocol and the use of the Ewald summation in calculating the electrostatic interactions, the RNA structure remained close to the crystal structure, as attested by a root-mean-square deviation below 2.5 A after 750 ps of simulation. All Watson-Crick base pairs were intact at the end of the simulations. The tertiary interactions, such as the sheared G.A pairs and the U-turn, important for the stabilisation of the three-dimensional RNA fold, were also retained. The results demonstrate that molecular dynamics simulations can be successfully used to investigate the dynamical behaviour of a ribozyme, thus, opening a road to study the role of transient structural changes involved in ribozyme catalysis.     

The tolerance to exchanges of the Watson Crick base pair in the hammerhead ribozyme core is determined by surrounding elements

Tertiary interacting elements are important features of functional RNA molecules, for example, in all small nucleolytic ribozymes. The recent crystal structure of a tertiary stabilized type I hammerhead ribozyme revealed a conventional Watson-Crick base pair in the catalytic core, formed between nucleotides C3 and G8. We show that any Watson-Crick base pair between these positions retains cleavage competence in two type III ribozymes. In the Arabidopsis thaliana sequence, only moderate differences in cleavage rates are observed for the different base pairs, while the peach latent mosaic viroid (PLMVd) ribozyme exhibits a preference for a pyrimidine at position 3 and a purine at position 8. To understand these differences, we created a series of chimeric ribozymes in which we swapped sequence elements that surround the catalytic core. The kinetic characterization of the resulting ribozymes revealed that the tertiary interacting loop sequences of the PLMVd ribozyme are sufficient to induce the preference for Y3-R8 base pairs in the A. thaliana hammerhead ribozyme. In contrast to this, only when the entire stem-loops I and II of the A. thaliana sequences are grafted on the PLMVd ribozyme is any Watson-Crick base pair similarly tolerated. The data provide evidence for a complex interplay of secondary and tertiary structure elements that lead, mediated by long-range effects, to an individual modulation of the local structure in the catalytic core of different hammerhead ribozymes.     

Solution structure and backbone dynamics of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin mutant: a molecular analysis of its reduced thermal stability

No general strategy for thermostability has been yet established, because the extra stability of thermophiles appears to be the sum of different cumulative stabilizing interactions. In addition, the increase of conformational rigidity observed in many thermophilic proteins, which in some cases disappears when mesophilic and thermophilic proteins are compared at their respective physiological temperatures, suggests that evolutionary adaptation tends to maintain corresponding states with respect to conformational flexibility. In this study, we accomplished a structural analysis of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin (BacTrx) mutant, which has reduced heat resistance with respect to the thermostable wild-type. Furthermore, we have also achieved a detailed study, carried out at 25, 45, and 65 degrees C, of the backbone dynamics of both the BacTrx and its K18G/R82E mutant. Our findings clearly indicate that the insertion of the two mutations causes a loss of energetically favorable long-range interactions and renders the secondary structure elements of the double mutants more similar to those of the mesophilic Escherichia coli thioredoxin. Moreover, protein dynamics analysis shows that at room temperature the BacTrx, as well as the double mutant, are globally as rigid as the mesophilic thioredoxins; differently, at 65 degrees C, which is in the optimal growth temperature range of A. acidocaldarius, the wild-type retains its rigidity while the double mutant is characterized by a large increase of the amplitude of the internal motions. Finally, our research interestingly shows that fast motions on the pico- to nanosecond time scale are not detrimental to protein stability and provide an entropic stabilization of the native state. This study further confirms that protein thermostability is reached through diverse stabilizing interactions, which have the key role to maintain the structural folding stable and functional at the working temperature.     

Self-Replication Reactions Dependent on Tertiary Interaction Motifs in an RNA Ligase Ribozyme

RNA can function both as an informational molecule and as a catalyst in living organisms. This duality is the premise of the RNA world hypothesis. However, one flaw in the hypothesis that RNA was the most essential molecule in primitive life is that no RNA self-replicating system has been found in nature. To verify whether RNA has the potential for self-replication, we constructed a new RNA self-assembling ribozyme that could have conducted an evolvable RNA self-replication reaction. The artificially designed, in vitro selected ligase ribozyme was employed as a prototype for a self-assembling ribozyme. The ribozyme is composed of two RNA fragments (form R1·Z1) that recognize another R1·Z1 molecule as their substrate and perform the high turnover ligation reaction via two RNA tertiary interaction motifs. Furthermore, the substrate recognition of R1·Z1 is tolerant of mutations, generating diversity in the corresponding RNA self-replicating network. Thus, we propose that our system implies the significance of RNA tertiary motifs in the early RNA molecular evolution of the RNA world.     

Secondary structure content of the HDV ribozyme in 95% formamide.

The Hepatitis Delta Virus (HDV) ribozyme self-cleaving activity in 20 M formamide solutions is unique. Does this catalytic activity result from the conservation of its tertiary structure in 20 M formamide? We followed the ribozyme structure in formamide solutions by monitoring the amount of bound Ethidium Bromide (EB). We were able to measure the quantity of dye bound using time-resolved fluorescence spectroscopy, as an estimate of the ribozyme double helical content. This method, calibrated by using oligonucleotides with defined tertiary structure and denaturing solvents, parallels NMR and UV measurements as a function of temperature. Measurements with the HDV ribozyme lead to three conclusions: (a) both the precursor and product RNAs are structured to 24 M (95% w/w) formamide or 4 M H2O solutions which is equivalent to 4 M H2O; (b) the HDV ribozyme is the only RNA sequence investigated in this study that retains so much structure in formamide; and (c) DNA analogs of formamide resistant HDV ribozyme sequences lose their structure at less than 15 M formamide. Thus, the structural integrity of the HDV ribozyme is an intrinsic property of the RNA molecule and its sequence.     

Evidence for processivity and two-step binding of the RNA substrate from studies of J1/2 mutants of the Tetrahymena ribozyme.

J1/2 of the Tetrahymena ribozyme, a sequence of three A residues, connects the RNA-binding site to the catalytic core. Addition or deletion of bases from J1/2 improves turnover and substrate specificity in the site-specific endonuclease reaction catalyzed by this ribozyme: G2CCCUCUA5 (S) + G in-equilibrium G2CCCUCU (P) + GA5. These paradoxical enhancements are caused by decreased affinity of the ribozyme for S and P [Young, B., Herschlag, D., & Cech, T.R. (1991) Cell 67, 1007]. An additional property of these mutant ribozymes, decreased fidelity of RNA cleavage, is now analyzed. (Fidelity is the ability to cleave at the correct phosphodiester bond within a particular RNA substrate.) Introduction of deoxy residues to give "chimeric" ribo/deoxyribooligonucleotides changes the positions of incorrect cleavage. Previous work indicated that S is bound to the ribozyme by both base pairing and teritary interactions involving 2'-hydroxyl groups of S. The data herein strongly suggest that the P1 duplex, which consists of S base-paired with the 5' exon binding site of the ribozyme, can dock into tertiary interactions in different registers; different 2'-hydroxyl groups of S plug into tertiary contacts with the ribozyme in the different registers. It is concluded that the mutations decrease fidelity by increasing the probability of docking out of register relative to docking in the normal register, thereby giving cleavage at different positions along S. These data also show that the contribution of J1/2 to the teritiary interactions is indirect, not direct. Thus, a structural role of the nonconserved J1/2 is indicated: this sequence positions S to optimize tertiary binding interactions and to ensure cleavage at the phosphodiester bond corresponding to the 5' splice site. Substitution of sulfur for the nonbridging pro-RP oxygen atom at the normal cleavage site has no effect on (kcat/Km)S but decreases the fraction of cleavage at the normal site in reactions catalyzed by the -3A mutant ribozyme, which has all three A residues of J1/2 removed. Thus, the ribozyme chooses where to cleave S after rate-limiting binding of S, indicating that docking can change after binding and suggesting that the ribozyme could act processively. Indeed, it is shown that the +2A ribozyme cleaves at one position along an RNA substrate and then, before releasing that RNA product, cleaves it again.(ABSTRACT TRUNCATED AT 400 WORDS)     

Maximizing RNA folding rates: a balancing act

Large ribozymes typically require very long times to refold into their active conformation in vitro, because the RNA is easily trapped in metastable misfolded structures. Theoretical models show that the probability of misfolding is reduced when local and long-range interactions in the RNA are balanced. Using the folding kinetics of the Tetrahymena ribozyme as an example, we propose that folding rates are maximized when the free energies of forming independent domains are similar to each other. A prediction is that the folding pathway of the ribozyme can be reversed by inverting the relative stability of the tertiary domains. This result suggests strategies for optimizing ribozyme sequences for therapeutics and structural studies.     

Comparative analysis of hairpin ribozyme structures and interference data

Great strides in understanding the molecular underpinnings of RNA catalysis have been achieved with advances in RNA structure determination by NMR spectroscopy and X-ray crystallography. Despite these successes the functional relevance of a given structure can only be assessed upon comparison with biochemical studies performed on functioning RNA molecules. The hairpin ribozyme presents an excellent case study for such a comparison. The active site is comprised of two stems each with an internal loop that forms a series of non-canonical base pairs. These loops dock into each other to create an active site for catalysis. Recently, three independent structures have been determined for this catalytic RNA, including two NMR structures of the isolated loop A and loop B stems and a high-resolution crystal structure of both loops in a docked conformation. These structures differ significantly both in their tertiary fold and the nature of the non-canonical base pairs formed within each loop. Several of the chemical groups required to achieve a functioning hairpin ribozyme have been determined by nucleotide analog interference mapping (NAIM). Here we compare the three hairpin structures with previously published NAIM data to assess the convergence between the structural and functional data. While there is significant disparity between the interference data and the individual NMR loop structures, there is almost complete congruity with the X-ray structure. The only significant differences cluster around an occluded pocket adjacent to the scissile phosphate. These local differences may suggest a role for these atoms in the transition state, either directly in chemistry or via a local structural rearrangement.     

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.     

Single VS ribozyme molecules reveal dynamic and hierarchical folding toward catalysis

Non-coding RNAs of complex tertiary structure are involved in numerous aspects of the replication and processing of genetic information in many organisms; however, an understanding of the complex relationship between their structural dynamics and function is only slowly emerging. The Neurospora Varkud Satellite (VS) ribozyme provides a model system to address this relationship. First, it adopts a tertiary structure assembled from common elements, a kissing loop and two three-way junctions. Second, catalytic activity of the ribozyme is essential for replication of VS RNA in vivo and can be readily assayed in vitro. Here we exploit single molecule FRET to show that the VS ribozyme exhibits previously unobserved dynamic and heterogeneous hierarchical folding into an active structure. Readily reversible kissing loop formation combined with slow cleavage of the upstream substrate helix suggests a model whereby the structural dynamics of the VS ribozyme favor cleavage of the substrate downstream of the ribozyme core instead. This preference is expected to facilitate processing of the multimeric RNA replication intermediate into circular VS RNA, which is the predominant form observed in vivo.