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.     

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.     

Conformational flexibility of four-way junctions in RNA

Helical junctions are common architectural features in RNA. They are particularly important in autonomously folding molecules, as exemplified by the hairpin ribozyme. We have used single-molecule fluorescence spectroscopy to study the dynamic properties of the perfect (4H) four-way helical junction derived from the hairpin ribozyme. In the presence of Mg(2+), the junction samples parallel and antiparallel conformations and both stacking conformers, with a bias towards one antiparallel stacking conformer. There is continual interconversion between the forms, such that there are several transitions per second under physiological conditions. Our data suggest that interconversion proceeds via an open intermediate with reduced cation binding in which coaxial stacking between helices is disrupted. The rate of interconversion becomes slower at higher Mg(2+) concentrations, yet the activation barrier decreases under these conditions, indicating that entropic effects are important. Transitions also occur in the presence of Na(+) only; however, the coaxial stacking appears incomplete under these conditions. The polymorphic and dynamic character of the four-way RNA junction provides a source of structural diversity, from which particular conformations required for biological function might be stabilised by additional RNA interactions or protein binding.     

Folding cooperativity in RNA and DNA is dependent on position in the helix

Secondary structural motifs play essential roles in the folding and function of RNA and DNA molecules. Previous work from our lab compared the folding of small DNA and RNA hairpin loops containing a sheared GA pair [Moody, E. M., Feerar, J. C., and Bevilacqua, P. C. (2004) Biochemistry 43, 7992-7998]. We found that the small DNA hairpins fold in a highly cooperative manner with indirect coupling, while their RNA counterparts fold in a much less cooperative fashion and display direct coupling. Herein, we extend this study to the double-stranded helix. We carried out double mutant cycles on base pairs having identical nearest-neighbor contexts but located in either external or internal helical registers. In the external register, both RNA and DNA exhibit extensive folding cooperativity between the penultimate and terminal base pair, which is independent of mismatch identity. In contrast, DNA exhibits virtually no folding cooperativity in the center of the helix, while RNA maintains substantial coupling, which is dependent on mismatch identity. Two models account for these non-nearest-neighbor effects: one involves the unfavorable entropy of helix initiation common to DNA and RNA, and the other involves steric and electrostatic strain peculiar to RNA. These data show that RNA can display cooperativity less than, greater than, or equal to that of DNA depending on context and position.     

Pyrite footprinting of RNA

In RNA, function follows form. Mapping the surface of RNA molecules with chemical and enzymatic probes has revealed invaluable information about structure and folding. Hydroxyl radicals (()OH) map the surface of nucleic acids by cutting the backbone where it is accessible to solvent. Recent studies showed that a microfluidic chip containing pyrite (FeS(2)) can produce sufficient ()OH to footprint DNA. The 49-nt Diels-Alder RNA enzyme catalyzes the C-C bond formation between a diene and a dienophile. A crystal structure, molecular dynamics simulation and atomic mutagenesis studies suggest that nucleotides of an asymmetric bulge participate in the dynamic architecture of the ribozyme's active center. Of note is that residue U42 directly interacts with the product in the crystallized RNA/product complex. Here, we use powdered pyrite held in a commercially available cartridge to footprint the Diels-Alderase ribozyme with single nucleotide resolution. Residues C39 to U42 are more reactive to ()OH than predicted by the solvent accessibility calculated from the crystal structure suggesting that this loop is dynamic in solution. The loop's flexibility may contribute to substrate recruitment and product release. Our implementation of pyrite-mediated ()OH footprinting is a readily accessible approach to gleaning information about the architecture of small RNA molecules.     

Efficient RNA ligation by reverse-joined hairpin ribozymes and engineering of twin ribozymes consisting of conventional and reverse-joined hairpin ribozyme units

In recent years major progress has been made in elucidating the mechanism and structure of catalytic RNA molecules, and we are now beginning to understand ribozymes well enough to turn them into useful tools. Work in our laboratory has focused on the development of twin ribozymes for site-specific RNA sequence alteration. To this end, we followed a strategy that relies on the combination of two ribozyme units into one molecule (hence dubbed twin ribozyme). Here, we present reverse-joined hairpin ribozymes that are structurally optimized and which, in addition to cleavage, catalyse efficient RNA ligation. The most efficient variant ligated its appropriate RNA substrate with a single turnover rate constant of 1.1 min(-1) and a final yield of 70%. We combined a reverse-joined hairpin ribozyme with a conventional hairpin ribozyme to create a twin ribozyme that mediates the insertion of four additional nucleotides into a predetermined position of a substrate RNA, and thus mimics, at the RNA level, the repair of a short deletion mutation; 17% of the initial substrate was converted to the insertion product.     

Spermine supports catalysis of hairpin ribozyme variants to differing extents

Because of the ability to cleave RNA substrates in trans, the hairpin ribozyme has great potential for therapeutic application. Activity of a three-stranded version of the minimal truncated form is enhanced by the presence of the polyamine spermine. Since spermine is the most abundant polyamine in eucariots, improved prospects for the hairpin ribozyme as therapeutic agent were predicted. We have found that not all hairpin ribozyme variants accept spermine equally well as counter-ion. Particularly the two-stranded versions commonly used for therapeutic studies show rather decreased activity when spermine is present. We have investigated a number of hairpin ribozyme derivatives regarding their ability to carry out spermine supported catalysis. Among the studied structures a two-stranded reverse-joined hairpin ribozyme displayed the highest cleavage rates in a synergistic mixture of magnesium ions and spermine. The specific features of this ribozyme along with its potential for in vivo application are discussed.     

Soft Interactions with Model Crowders and Non-canonical Interactions with Cellular Proteins Stabilize RNA Folding

Living cells contain diverse biopolymers, creating a heterogeneous crowding environment, the impact of which on RNA folding is poorly understood. Here, we have used single-molecule fluorescence resonance energy transfer to monitor tertiary structure formation of the hairpin ribozyme as a model to probe the effects of polyethylene glycol and yeast cell extract as crowding agents. As expected, polyethylene glycol stabilizes the docked, catalytically active state of the ribozyme, in part through excluded volume effects; unexpectedly, we found evidence that it additionally displays soft, non-specific interactions with the ribozyme. Yeast extract has a profound effect on folding at protein concentrations 1000-fold lower than found intracellularly, suggesting the dominance of specific interactions over volume exclusion. Gel shift assays and affinity pull-down followed by mass spectrometry identified numerous non-canonical RNA-binding proteins that stabilize ribozyme folding; the apparent chaperoning activity of these ubiquitous proteins significantly compensates for the low-counterion environment of the cell.     

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.     

In vitro selection of adenine-dependent hairpin ribozymes

Adenine-dependent hairpin ribozymes were isolated by in vitro selection from a degenerated hairpin ribozyme population. Two new adenine-dependent ribozymes catalyze their own reversible cleavage in the presence of free adenine. Both aptamers have Mg(2+) requirements for adenine-assisted cleavage similar to the wild-type hairpin ribozyme. Cleavage kinetics studies in the presence of various other small molecules were compared. The data suggest that adenine does not induce RNA self-cleavage in the same manner for both aptamers. In addition, investigations of pH effects on catalytic rates show that both adenine-dependent aptamers are more active in basic conditions, suggesting that they use new acid/base catalytic strategies in which adenine could be involved directly. The discovery of hairpin ribozymes dependent on adenine for their reversible self-cleavage presents considerable biochemical and evolutionary interests because we show that RNA is able to use exogenous reactive molecules to enhance its own catalytic activity. Such a mechanism may have been a means by which the ribozymes of the RNA world enlarged their chemical repertoire.     

Substrate shape specificity of E coli RNase P ribozyme is dependent on the concentration of magnesium ion

The bacterial RNase P ribozyme can accept a hairpin RNA with CCA-3' tag sequence as well as a cloverleaf pre-tRNA as substrate in vitro, but the details are not known. By switching tRNA structure using an antisense guide DNA technique, we examined the Escherichia coli RNase P ribozyme specificity for substrate RNA of a given shape. Analysis of the RNase P reaction with various concentrations of magnesium ion revealed that the ribozyme cleaved only the cloverleaf RNA at below 10 mM magnesium ion. At 10 mM magnesium ion or more, the ribozyme also cleaved a hairpin RNA with a CCA-3' tag sequence. At above 20 mM magnesium ion, cleavage site wobbling by the enzyme in tRNA-derived hairpin occurred, and the substrate specificity of the enzyme became broader. Additional studies using another hairpin substrate demonstrated the same tendency. Our data strongly suggest that raising the concentration of metal ion induces a conformational change in the RNA enzyme.     

The many faces of the hairpin ribozyme: structural and functional variants of a small catalytic RNA

The hairpin ribozyme is a small catalytic RNA that has been reengineered resulting in a number of variants with extended or even new functions. Thus, manipulation of the hairpin ribozyme structure has allowed for activity control by external effectors, namely oligonucleotides, flavine mononucleotide, and adenine. Hairpin ribozyme-derived twin ribozymes that mediate RNA fragment exchange reactions as well as self-processing hairpin ribozymes were designed. Furthermore, several hairpin ribozyme variants have been engineered for knock down of specific RNA substrates by adapting the substrate-binding domain to the specific target sequence. This review will focus on hairpin ribozymes possessing structural extensions/variations and thus functionally differing from the parent hairpin ribozyme.     

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.     

Hairpins under tension: RNA versus DNA

We use optical tweezers to control the folding and unfolding of individual DNA and RNA hairpins by force. Four hairpin molecules are studied in comparison: two DNA and two RNA ones. We observe that the conformational dynamics is slower for the RNA hairpins than for their DNA counterparts. Our results indicate that structures made of RNA are dynamically more stable. This difference might contribute to the fact that DNA and RNA play fundamentally different biological roles in spite of chemical similarity.     

Imaging of single hairpin ribozymes in solution by atomic force microscopy.

The hairpin ribozyme is a short endonucleolytic RNA motif isolated from a family of related plant virus satellite RNAs. It consists of two independently folding domains, each comprising two Watson-Crick helices flanking a conserved internal loop. The domains need to physically interact (dock) for catalysis of site-specific cleavage and ligation reactions. Using tapping-mode atomic force microscopy in aqueous buffer solution, we were able to produce high quality images of individual hairpin ribozyme molecules with extended terminal helices. Three RNA constructs with either the essential cleavage site guanosine or a detrimental adenosine substitution and with or without a 6-nt insertion to confer flexibility to the interdomain hinge show structural differences that correlate with their ability to form the active docked conformation. The observed contour lengths and shapes are consistent with previous bulk-solution measurements of the transient electric dichroism decays for the same RNA constructs. The active docked construct appears as an asymmetrically docked conformation that might be an indication of a more complicated docking event than a simple collapse around the interdomain hinge.     

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 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".     

Mutational inhibition of ligation in the hairpin ribozyme: substitutions of conserved nucleobases A9 and A10 destabilize tertiary structure and selectively promote cleavage

The hairpin ribozyme acts as a reversible, site-specific endoribonuclease that ligates much more rapidly than it cleaves cognate substrate. While the reaction pathway for ligation is the reversal of cleavage, little is known about the atomic and electrostatic details of the two processes. Here, we report the functional consequences of molecular substitutions of A9 and A10, two highly conserved nucleobases located adjacent to the hairpin ribozyme active site, using G, C, U, 2-aminopurine, 2,6-diaminopurine, purine, and inosine. Cleavage and ligation kinetics were analyzed, tertiary folding was monitored by hydroxyl radical footprinting, and interdomain docking was studied by native gel electrophoresis. We determined that nucleobase substitutions that exhibit significant levels of interference with tertiary folding and interdomain docking have relatively large inhibitory effects on ligation rates while showing little inhibition of cleavage. Indeed, one variant, A10G, showed a fivefold enhancement of cleavage rate and no detectable ligation, and we suggest that this property may be uniquely well suited to intracellular targeted RNA cleavage applications. Results support a model in which formation of a kinetically stable tertiary structure is essential for ligation of the hairpin ribozyme, but is not necessary for cleavage.