Inter-domain cross-linking and molecular modelling of the hairpin ribozyme

Identification of over 200-fold more hairpin ribozymes than previously known in diverse circular RNAs

Self-cleaving ribozymes are catalytic RNAs that cut themselves at a specific inter-nucleotide linkage. They serve as a model of RNA catalysis, and as an important tool in biotechnology. For most of the nine known structural classes of self-cleaving ribozymes, at least hundreds of examples are known, and some are present in multiple domains of life. By contrast, only four unique examples of the hairpin ribozyme class are known, despite its discovery in 1986. We bioinformatically predicted 941 unique hairpin ribozymes of a different permuted form from the four previously known hairpin ribozymes, and experimentally confirmed several diverse predictions. These results profoundly expand the number of natural hairpin ribozymes, enabling biochemical analysis based on natural sequences, and suggest that a distinct permuted form is more biologically relevant. Moreover, all novel hairpins were discovered in metatranscriptomes. They apparently reside in RNA molecules that vary both in size—from 381 to 5170 nucleotides—and in protein content. The RNA molecules likely replicate as circular single-stranded RNAs, and potentially provide a dramatic increase in diversity of such RNAs. Moreover, these organisms have eluded previous attempts to isolate RNA viruses from metatranscriptomes—suggesting a significant untapped universe of viruses or other organisms hidden within metatranscriptome sequences.

RNA structure and dynamics: a base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson-Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson-Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.

Predicting and modeling RNA architecture

A general approach for modeling the architecture of large and structured RNA molecules is described. The method exploits the modularity and the hierarchical folding of RNA architecture that is viewed as the assembly of preformed double-stranded helices defined by Watson-Crick base pairs and RNA modules maintained by non-Watson-Crick base pairs. Despite the extensive molecular neutrality observed in RNA structures, specificity in RNA folding is achieved through global constraints like lengths of helices, coaxiality of helical stacks, and structures adopted at the junctions of helices. The Assemble integrated suite of computer tools allows for sequence and structure analysis as well as interactive modeling by homology or ab initio assembly with possibilities for fitting within electronic density maps. The local key role of non-Watson-Crick pairs guides RNA architecture formation and offers metrics for assessing the accuracy of three-dimensional models in a more useful way than usual root mean square deviation (RMSD) values.

Exploring RNA structure by integrative molecular modelling

RNA molecular modelling is adequate to rapidly tackle the structure of RNA molecules. With new structured RNAs constituting a central class of cellular regulators discovered every year, the need for swift and reliable modelling methods is more crucial than ever. The pragmatic method based on interactive all-atom molecular modelling relies on the observation that specific structural motifs are recurrently found in RNA sequences. Once identified by a combination of comparative sequence analysis and biochemical data, the motifs composing the secondary structure of a given RNA can be extruded in three dimensions (3D) and used as building blocks assembled manually during a bioinformatic interactive process. Comparing the models to the corresponding crystal structures has validated the method as being powerful to predict the RNA topology and architecture while being less accurate regarding the prediction of base-base interactions. These aspects as well as the necessary steps towards automation will be discussed.

Synthesis of guanosine and deoxyguanosine phosphoramidites with cross-linkable thioalkyl tethers for direct incorporation into RNA and DNA

We describe the synthesis of protected phosphoramidites of deoxyriboguanosine and guanosine derivatives containing a thiopropyl tether at the guanine N2 (7a,b) for site-specific crosslinking from the minor groove of either DNA or RNA to a thiol of a protein or another nucleic acid. The thiol is initially protected as a tert-butyl disulfide that is stable during oligonucleotide synthesis. While the completed oligonucleotide is still attached to the support, or after purification, the tert-butyl thiol can readily be removed or replaced by thioethylamine or 5-thio-2-nitrobenzoic acid, which have more favorable crosslinking rates.

DNA aptamer-mediated regulation of the hairpin ribozyme by human alpha-thrombin

The combination of specific ligand-binding aptamers with hairpin ribozyme catalysis generates molecules that can be controlled by external factors. Here we have generated hairpin ribozymes that can be regulated by a short DNA aptamer specific for human alpha-thrombin. This was achieved by constructing a ribozyme variant harboring an RNA sequence complementary to the aptamer, to which the aptamer can hybridize forming a heteroduplex. In this way, the DNA aptamer completely abolishes the catalytic activity of the ribozyme, due to the formation of an inactive ribozyme conformation. However, in the presence of the aptamer's target protein human alpha-thrombin, the inhibitory effect of the DNA aptamer is competitively neutralized and the ribozyme is activated in a highly specific fashion. Protein-responsive allosteric ribozymes are proposed to act as tools with potential applications in medicine where fast detection of clinically relevant targets is required.

Topological rearrangement yields structural stabilization and interhelical distance constraints in the Kin.46 self-phosphorylating ribozyme

The Kin.46 ribozyme catalyzes transfer of the gamma (thio)phosphoryl group of ATP (or ATPgammaS) to the ribozyme's 5' hydroxyl. Single-turnover catalytic activities of topologically rearranged versions of Kin.46 were studied to gain insight into its overall tertiary architecture. The distal ends of stems P3 and P4 were tethered through a single-stranded connection domain that altered the interhelical connectivity. The shortest linkers interfered with catalysis, while seven or more nucleotides (nt) in the linker allowed near-normal catalytic rates, suggesting that a distance of roughly 25-35 A optimally separates the termini of these helices. Activity was maximal when the tether contained 15 nt, at which point the k(cat) (0.016 min(-1)) and Km (1.2 mM) values were identical to those of a nontethered control. The presence of the tether alters Mg(2+) dependence, in that Mg2+ binding appears to be more cooperative in the tethered ribozyme (Hill coefficient 1.4-1.8 versus 0.8 for the nontethered ribozyme). Binding affinity for the ATPgammaS substrate increases at elevated concentrations of Mg2+, particularly for the tethered ribozyme. The tethered ribozyme displays significantly enhanced thermal stability, with a maximum initial velocity (0.126 min(-1)) at 60 degrees C, whereas the nontethered ribozyme has a lower maximum initial velocity (0.051 min(-1)) at 50 degrees C. The tether also significantly reduces the apparent entropy of activation. Both of these effects can be understood in terms of stabilization of the ribozyme in a conformation that is on-path with respect to catalysis, and in terms of facilitating formation of the allosteric activation helix P4.

The interaction networks of structured RNAs

All pairwise interactions occurring between bases which could be detected in three-dimensional structures of crystallized RNA molecules are annotated on new planar diagrams. The diagrams attempt to map the underlying complex networks of base–base interactions and, especially, they aim at conveying key relationships between helical domains: co-axial stacking, bending and all Watson–Crick as well as non-Watson–Crick base pairs. Although such wiring diagrams cannot replace full stereographic images for correct spatial understanding and representation, they reveal structural similarities as well as the conserved patterns and distances between motifs which are present within the interaction networks of folded RNAs of similar or unrelated functions. Finally, the diagrams could help devising methods for meaningfully transforming RNA structures into graphs amenable to network analysis.

Conformational heterogeneity at position U37 of an all-RNA hairpin ribozyme with implications for metal binding and the catalytic structure of the S-turn

The hairpin ribozyme is an RNA enzyme that performs site-specific phosphodiester bond cleavage between nucleotides A-1 and G+1 within its cognate substrate. Previous functional studies revealed that the minimal hairpin ribozyme exhibited "gain-of-function" cleavage properties resulting from U39C or U39 to propyl linker (C3) modifications. Furthermore, each "mutant" displayed different magnesium-dependence in its activity. To investigate the molecular basis for these gain-of-function variants, crystal structures of minimal, junctionless hairpin ribozymes were solved in native (U39), and mutant U39C and U39(C3) forms. The results revealed an overall molecular architecture comprising two docked internal loop domains folded into a wishbone shape, whose tertiary interface forms a sequestered active site. All three minimal hairpin ribozymes bound Co(NH(3))(6)(3+) at G21/A40, the E-loop/S-turn boundary. The native structure also showed that U37 of the S-turn adopts both sequestered and exposed conformations that differ by a maximum displacement of 13 A. In the sequestered form, the U37 base packs against G36, and its 2'-hydroxyl group forms a water mediated hydrogen bond to O4' of G+1. These interactions were not observed in previous four-way-junction hairpin ribozyme structures due to crystal contacts with the U1A splicing protein. Interestingly, the U39C and U39(C3) mutations shifted the equilibrium conformation of U37 into the sequestered form through formation of new hydrogen bonds in the S-turn, proximal to the essential nucleotide A38. A comparison of all three new structures has implications for the catalytically relevant conformation of the S-turn and suggests a rationale for the distinctive metal dependence of each mutant.

Competitive regulation of modular allosteric aptazymes by a small molecule and oligonucleotide effector

The hairpin ribozyme can catalyze the cleavage of RNA substrates by employing its conformational flexibility. To form a catalytic complex, the two domains A and B of the hairpin-ribozyme complex must interact with one another in a folding step called docking. We have constructed hairpin ribozyme variants harboring an aptamer sequence that can be allosterically induced by flavin mononucleotide (FMN). Domains A and B are separated by distinct bridge sequences that communicate the formation of the FMN-aptamer complex to domains A and B, facilitating their docking. In the presence of a short oligonucleotide that is complementary to the aptamer, catalytic activity of the ribozyme is completely abolished, due to the formation of an extended conformer that cannot perform catalysis. However, in the presence of the small molecule effector FMN, the inhibitory effect of the oligonucleotide is competitively neutralized and the ribozyme is activated 150-fold. We thus have established a new principle for the regulation of ribozyme catalysis in which two regulatory factors (an oligonucleotide and a small molecule) that switch the ribozyme's activity in opposite directions compete for the same binding site in the aptamer domain.

Investigation of Mg2+- and temperature-dependent folding of the hairpin ribozyme by photo-crosslinking: effects of photo-crosslinker tether length and chemistry

We have used photo-crosslinking to investigate the structure and dynamics of four-way junction hairpin ribozyme constructs. Four phenylazide photo-crosslinkers were coupled to 2′-NH₂-modified U+2 in the substrate and irradiated at different Mg²⁺ concentrations and temperatures. Consistent with the role of divalent metal ions in hairpin ribozyme folding, we observed more interdomain crosslinks in the presence of Mg²⁺ than in its absence. In general, we observed intradomain crosslinks to nucleotides 2–11 and interdomain crosslinks to the U1A binding loop. Crosslinks to A26 and G36 in domain B were also observed when crosslinking was carried out at −78°C. In contrast to crosslinking results at higher temperatures (0, 25 and 37°C), similar crosslinks were obtained in the presence and absence of Mg²⁺ at −78°C, suggesting Mg²⁺ stabilizes a low-energy hairpin ribozyme conformation. We also evaluated the effects of photo-crosslinker structure and mechanism on crosslinks. First, most crosslinks were to unpaired nucleotides. Second, shorter and longer photo-crosslinkers formed crosslinks to intradomain locations nearer to and farther from photo-crosslinker modification, respectively. Finally, fluorine substitutions on the phenylazide ring did not change the locations of crosslinks, but rather decreased crosslinking efficiency. These findings have implications for the use of phenylazide photo-crosslinkers in structural studies of RNA.

Sensing complex regulatory networks by conformationally controlled hairpin ribozymes

The hairpin ribozyme catalyses RNA cleavage by a mechanism utilizing its conformational flexibility during the docking of two independently folded internal loop domains A and B. Based on this mechanism, we designed hairpin ribozyme variants that can be induced or repressed by external effector oligonucleotides influencing the docking process. We incorporated a third domain C to assimilate alternate stable RNA motifs such as a pseudo-half-knot or an internal stem-loop structure. Small sequence changes in domain C allowed targeted switching of ribozyme activity: the same effector oligonucleotide can either serve as an inducer or repressor. The ribozymes were applied to trp leader mRNA, the RNA sequence tightly bound by l-tryptophan-activated trp-RNA-binding attenuation protein (TRAP). When domain C is complementary to this mRNA, ribozyme activity can be altered by annealing trp leader mRNA, then specifically reverted by its TRAP/tryptophan-mediated sequestration. This approach allows to precisely sense the activity status of a protein controlled by its metabolite molecule.

Modifications and deletions of helices within the hairpin ribozyme-substrate complex: an active ribozyme lacking helix 1

Within the hairpin ribozyme, structural elements required for formation of the active tertiary structure are localized in two independently folding domains, each consisting of an internal loop flanked by helical elements. Here, we present results of a systematic examination of the relationship between the structure of the helical elements and the ability of the RNA to form the catalytically active tertiary structure. Deletions and mutational analyses indicate that helix 1 (H1) in domain A can be entirely eliminated, while segments of helices 2, 3, and 4 can also be deleted. From these results, we derive a new active minimal ribozyme that contains three helical elements, an internal loop, and a terminal loop. A three-dimensional model of this truncated ribozyme was generated using MC-SYM, and confirms that the catalytic core of the minimized construct can adopt a tertiary structure that is very similar to that of the nontruncated version. A new strategy is described to study the functional importance of various residues and chemical groups and to identify specific interdomain interactions. This approach uses two physically separated and truncated domains derived from the minimal motif.

In vitro selection of second site revertants analysis of the hairpin ribozyme active site

We have used in vitro genetics to evaluate the function and interactions of the conserved base G8 in the hairpin ribozyme catalytic RNA. Second site revertant selection for a G8X mutant, where X is any of the other three natural nucleobases, yielded a family of second site suppressors of the G8U mutant, but not of G8C or G8A, indicating that only G and U can be tolerated at position 8 of the ribozyme. This result is consistent with recent observations that point to the functional importance of G8 N-1 in the chemistry of catalysis by this ribozyme reaction. Suppression of the G8U mutation was observed when changes were made directly across loop A from the mutated base at substrate position +2 or positions +2 and +3 in combination. The same changes made in the context of the natural G8 sequence resulted in a very large drop in activity. Thus, the G8U mutation results in a change in specificity of the ribozyme from 5'-N / GUC-3' to 5'-N / GCU-3'. The results presented imply that G8 interacts directly with U+2 during catalysis. We propose that this interaction favors the correct positioning of the catalytic determinants of G8. The implications for the folding of the ribozyme and the catalytic mechanism are discussed.

Ribozymes: recent advances in the development of RNA tools

The discovery 20 years ago that some RNA molecules, called ribozymes, are able to catalyze chemical reactions was a breakthrough in biology. Over the last two decades numerous natural RNA motifs endowed with catalytic activity have been described. They all fit within a few well-defined types that respond to a specific RNA structure. The prototype catalytic domain of each one has been engineered to generate trans-acting ribozymes that catalyze the site-specific cleavage of other RNA molecules. On the 20th anniversary of ribozyme discovery we briefly summarize the main features of the different natural catalytic RNAs. We also describe progress towards developing strategies to ensure an efficient ribozyme-based technology, dedicating special attention to the ones aimed to achieve a new generation of therapeutic agents.

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.

Developments in RNA chemistry, a personal view

This review describes some of the contributions of chemistry to the RNA field with a personal bias towards the phosphorothioate modification and the derivatives at the ribose 2'-position. The usefulness of these modifications is discussed and documented with some examples.

In the fluorescent spotlight: global and local conformational changes of small catalytic RNAs

RNA is a ubiquitous biopolymer that performs a multitude of essential cellular functions involving the maintenance, transfer, and processing of genetic information. RNA is unique in that it can carry both genetic information and catalytic function. Its secondary structure domains, which fold stably and independently, assemble hierarchically into modular tertiary structures. Studies of these folding events are key to understanding how catalytic RNAs (ribozymes) are able to position reaction components for site-specific chemistry. We have made use of fluorescence techniques to monitor the rates and free energies of folding of the small hairpin and hepatitis delta virus (HDV) ribozymes, found in satellite RNAs of plant and the human hepatitis B viruses, respectively. In particular, fluorescence resonance energy transfer (FRET) has been employed to monitor global conformational changes, and 2-aminopurine fluorescence quenching to probe for local structural rearrangements. In this review we illuminate what we have learned about the reaction pathways of the hairpin and HDV ribozymes, and how our results have complemented other biochemical and biophysical investigations. The structural transitions observed in these two small catalytic RNAs are likely to be found in many other biological RNAs, and the described fluorescence techniques promise to be broadly applicable.

Correlating structural dynamics and function in single ribozyme molecules

We have studied the correlation between structural dynamics and function of the hairpin ribozyme. The enzyme-substrate complex exists in either docked (active) or undocked (inactive) conformations. Using single-molecule fluorescence methods, we found complex structural dynamics with four docked states of distinct stabilities and a strong memory effect where each molecule rarely switches between different docked states. We also found substrate cleavage to be rate-limited by a combination of conformational transitions and reversible chemistry equilibrium. The complex structural dynamics quantitatively explain the heterogeneous cleavage kinetics common to many catalytic RNAs. The intimate coupling of structural dynamics and function is likely a general phenomenon for RNA.

Rescue of an abasic hairpin ribozyme by cationic nucleobases: evidence for a novel mechanism of RNA catalysis

The hairpin ribozyme catalyzes a reversible phosphodiester cleavage reaction. We examined the roles of conserved nucleobases in catalysis using an abasic ribozyme rescue strategy. Loss of the active site G8 nucleobase reduced the cleavage rate constant by 350-fold while loss of A9 and A10 nucleobases reduced activity less than 10-fold. Certain heterocyclic amines restored partial activity when provided in solution to the variant lacking G8. Heterocyclic amines that were capable of rescue shared the exocyclic amine and cyclic amide in common with the Watson-Crick hydrogen bonding face of guanine. In contrast to the shallow pH dependence of unmodified ribozyme activity, rescue activity increased sharply with decreasing pH. These results support a novel model for RNA catalysis in which a cationic nucleobase contributes electrostatic stabilization to negative charge developing in the transition state.

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.

The contribution of 2'-hydroxyls to the cleavage activity of the Neurospora VS ribozyme

We have used nucleotide analog interference mapping and site-specific substitution to determine the effect of 2'-deoxynucleotide substitution of each nucleotide in the VS ribozyme on the self-cleavage reaction. A large number of 2'-hydroxyls (2'-OHs) that contribute to cleavage activity of the VS ribozyme were found distributed throughout the core of the ribozyme. The locations of these 2'-OHs in the context of a recently developed helical orientation model of the VS ribozyme suggest roles in multi-stem junction structure, helix packing, internal loop structure and catalysis. The functional importance of three separate 2'-OHs supports the proposal that three uridine turns contribute to local and long-range tertiary structure formation. A cluster of important 2'-OHs near the loop that is the candidate region for the active site and one very important 2'-OH in the loop that contains the cleavage site confirm the functional importance of these two loops. A cluster of important 2'-OHs lining the minor groove of stem-loop I and helix II suggests that these regions of the backbone may play an important role in positioning helices in the active structure of the ribozyme.

Novel spermine-amino acid conjugates and basic tripeptides enhance cleavage of the hairpin ribozyme at low magnesium ion concentration

Combinations of the polyamine spermine and magnesium ions synergize to dramatically enhance cleavage of the hairpin ribozyme. Certain synthetic basic tripeptides stimulate hairpin cleavage significantly at limiting magnesium ion concentration, notably the tripeptide of L-diaminobutyric acid (Dab). Of a range of novel synthetic spermine-amino acid conjugates, L-Dab-spermine (but not D-Dab nor other amino acid conjugates) was more effective than spermine itself.

Ribozyme structures and mechanisms

The past few years have seen exciting advances in understanding the structure and function of catalytic RNA. Crystal structures of several ribozymes have provided detailed insight into the folds of RNA molecules. Models of other biologically important RNAs have been constructed based on structural, phylogenetic, and biochemical data. However, many questions regarding the catalytic mechanisms of ribozymes remain. This review compares the structures and possible catalytic mechanisms of four small self-cleaving RNAs: the hammerhead, hairpin, hepatitis delta virus, and in vitro-selected lead-dependent ribozymes. The organization of these small catalysts is contrasted to that of larger ribozymes, such as the group I intron.

Engineering disulfide cross-links in RNA via air oxidation

This unit presents protocols for the synthesis of alkylthiol-modified ribonucleosides, their incorporation into synthetic RNA, and the formation of intramolecular disulfide bonds in RNA by air oxidation. The disulfide bonds can be formed in quantitative yields between thiols positioned in close proximity by virtue of either the secondary or tertiary structure of the RNA. Disulfide cross-links are useful tools to probe solution structures of RNA, to monitor dynamic motion, to stabilize folded RNAs, and to study the process of tertiary structure folding.

Crystal structure of a hairpin ribozyme-inhibitor complex with implications for catalysis

The hairpin ribozyme catalyses sequence-specific cleavage of RNA. The active site of this natural RNA results from the docking of two irregular helices: stems A and B. One strand of stem A harbours the scissile bond. The 2.4 A resolution structure of a hairpin ribozyme-inhibitor complex reveals that the ribozyme aligns the 2'-OH nucleophile and the 5'-oxo leaving group by twisting apart the nucleotides that flank the scissile phosphate. The base of the nucleotide preceding the cleavage site is stacked within stem A; the next nucleotide, a conserved guanine, is extruded from stem A and accommodated by a highly complementary pocket in the minor groove of stem B. Metal ions are absent from the active site. The bases of four conserved purines are positioned potentially to serve as acid-base catalysts. This is the first structure determination of a fully assembled ribozyme active site that catalyses a phosphodiester cleavage without recourse to metal ions.

The hairpin ribozyme substrate binding-domain: a highly constrained D-shaped conformation

The two domains of the hairpin ribozyme-substrate complex, usually depicted as straight structural elements, must interact with one another in order to form an active conformation. Little is known about the internal geometry of the individual domains in an active docked complex. Using various crosslinking and structural approaches in conjunction with molecular modeling (constraint-satisfaction program MC-SYM), we have investigated the conformation of the substrate-binding domain in the context of the active docked ribozyme-substrate complex. The model generated by MC-SYM showed that the domain is not straight but adopts a bent conformation (D-shaped) in the docked state of the ribozyme, indicating that the two helices bounding the internal loop are closer than was previously assumed. This arrangement rationalizes the observed ability of hairpin ribozymes with a circularized substrate-binding strand to cleave a circular substrate, and provides essential information concerning the organization of the substrate in the active conformation. The internal geometry of the substrate-binding strand places G8 of the substrate-binding strand near the cleavage site, which has allowed us to predict the crucial role played by this nucleotide in the reaction chemistry.

Ribozyme structures and mechanisms

The past few years have seen exciting advances in understanding the structure and function of catalytic RNA. Crystal structures of several ribozymes have provided detailed insight into the folds of RNA molecules. Models of other biologically important RNAs have been constructed based on structural, phylogenetic, and biochemical data. However, many questions regarding the catalytic mechanisms of ribozymes remain. This review compares the structures and possible catalytic mechanisms of four small self-cleaving RNAs: the hammerhead, hairpin, hepatitis delta virus, and in vitro-selected lead-dependent ribozymes. The organization of these small catalysts is contrasted to that of larger ribozymes, such as the group I intron.

A base change in the catalytic core of the hairpin ribozyme perturbs function but not domain docking

The hairpin ribozyme is a small endonucleolytic RNA motif with potential for targeted RNA inactivation. It optimally cleaves substrates containing the sequence 5'-GU-3' immediately 5' of G. Previously, we have shown that tertiary structure docking of its two domains is an essential step in the reaction pathway of the hairpin ribozyme. Here we show, combining biochemical and fluorescence structure and function probing techniques, that any mutation of the substrate base U leads to a docked RNA fold, yet decreases cleavage activity. The docked mutant complex shares with the wild-type complex a common interdomain distance as measured by time-resolved fluorescence resonance energy transfer (FRET) as well as the same solvent-inaccessible core as detected by hydroxyl-radical protection; hence, the mutant complex appears nativelike. FRET experiments also indicate that mutant docking is kinetically more complex, yet with an equilibrium shifted toward the docked conformation. Using 2-aminopurine as a site-specific fluorescent probe in place of the wild-type U, a local structural rearrangement in the substrate is observed. This substrate straining accompanies global domain docking and involves unstacking of the base and restriction of its conformational dynamics, as detected by time-resolved 2-aminopurine fluorescence spectroscopy. These data appear to invoke a mechanism of functional interference by a single base mutation, in which the ribozyme-substrate complex becomes trapped in a nativelike fold preceding the chemical transition state.

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.

Folding of branched RNA species

The global structures of branched RNA species are important to their function. Branched RNA species are defined as molecules in which double-helical segments are interrupted by abrupt discontinuities. These include helical junctions of different orders, and base bulges and loops. Common helical junctions are three- and four-way junctions, often interrupted by mispairs or additional nucleotides. There are many interesting examples of functional RNA junctions, including the hammerhead and hairpin ribozymes, and junctions that serve as binding sites for proteins. The junctions display some common structural properties. These include a tendency to undergo pairwise helical stacking and ion-induced conformational transitions. Helical branchpoints can act as key architectural components and as important sites for interactions with proteins. Copyright 1999 John Wiley & Sons, Inc.

Adenine protonation in domain B of the hairpin ribozyme

Protein enzymes often use ionizable side chains, such as histidine, for general acid-base catalysis because the imidazole pK(a) is near neutral pH. RNA enzymes, on the other hand, are comprised of nucleotides which do not have apparent pK(a) values near neutral pH. Nevertheless, it has been recently shown that cytidine and adenine protonation can play an important role in both nucleic acid structure and catalysis. We have employed heteronuclear NMR methods to determine the pK(a) values and time scales of chemical exchanges associated with adenine protonation within the catalytically essential B domain of the hairpin ribozyme. The large, adenine-rich internal loop of the B domain allows us to determine adenine pK(a) values for a variety of non-Watson-Crick base pairs. We find that adenines within the internal loop have pK(a) values ranging from 4.8 to 5.8, significantly higher than the free mononucleotide pK(a) of 3. 5. Adenine protonation results in potential charge stabilization, hydrogen bond formation, and stacking interactions that are expected to stabilize the internal loop structure at low pH. Fast proton exchange times of 10-50 micros were determined for the well-resolved adenines. These results suggest that shifted pK(a) values may be a common feature of adenines in non-Watson-Crick base pairs, and identify two adenines which may participate in hairpin ribozyme active site chemistry.

In vitro selection of a ligase ribozyme carrying alkylamino groups in the side chains

A novel ligase ribozyme was in vitro selected from a random-sequence RNA library including N(6)-aminohexyl-modified adenine residues in place of natural adenine residues. This ribozyme mediated the formation of a phosphodiester bond with a DNA oligonucleotide through condensation with a 5'-triphosphate moiety on the ribozyme. Among the clones isolated from this selection, one was shown to accelerate ligation about 250-fold compared to the original random-sequence RNA library. Almost no rate acceleration was observed when the N(6)-aminohexyl-groups on adenine residues were omitted. Furthermore, ligation was dependent on the presence of a template DNA oligonucleotide that bridged the two strands.

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.

Analysis of the functional role of a G.A sheared base pair by in vitro genetics

A classical genetic strategy has been combined with an in vitro selection method to search for functional interactions between the two domains of the hairpin ribozyme. G(21) is located within internal loop B; it is proposed to form a sheared base pair with A(43) across loop B and to bind a Mg(2+) ion. Both nucleotides are important for ribozyme function, and G.A sheared base pairs are a very widespread motif in structured RNA. We took advantage of its presence in the hairpin ribozyme to study its functional role. Pseudorevertants, in which the loss of G(21) was compensated by mutations at other positions, were isolated by in vitro selection. The vast majority of G(21) revertants contained substitutions within domain A, pointing to functional communication between specific sites within the two domains of the hairpin ribozyme. The possibility of a direct or redundant contacts is supported by electrophoretic mobility shift studies showing that a complex formed between domain B of the ribozyme and the substrate was disrupted and restored by base substitutions that have analogous effects on catalytic activity. The functional significance of this complex, the role of the nucleotides involved, and the basis for magnesium ion requirement is discussed.

Design, synthesis, and analysis of conformationally constrained nucleic acids

In this review I discuss straightforward and general methods to modify nucleic acid structure with disulfide cross-links. A motivating factor in developing this chemistry was the notion that disulfide bonds would be excellent tools to probe the structure, dynamics, thermodynamics, folding, and function of DNA and RNA, much in the way that cystine cross-links have been used to study proteins. The chemistry described has been used to synthesize disulfide cross-linked hairpins and duplexes, higher order structures like triplexes, nonground-state conformations, and tRNAs. Since the cross-links form quantitatively by mild air oxidation and do not perturb either secondary or tertiary structure, this modification should prove quite useful for the study of nucleic acids.

Investigation of the proposed interdomain ribose zipper in hairpin ribozyme cleavage using 2'-modified nucleosides

The hairpin ribozyme achieves catalytic cleavage through interaction of essential nucleotides located in two distinct helical domains that include internal loops. Initial docking of the two domains is ion dependent and appears to be followed by a structural rearrangement that allows the ribozyme to achieve a catalytically active state that can undergo cleavage. The proposed structural rearrangement may also be ion dependent and is now of increased importance due to recent evidence that docking is not rate limiting and that metal ions are unlikely to be involved in the chemical cleavage step. An initial structural model of the docked hairpin ribozyme included a proposal for a ribose zipper motif that involves two pairs of hydroxyl groups at A(10) and G(11) in domain A pairing with C(25) and A(24) in domain B, respectively. We have used a chemical functional group substitution technique to study whether this proposed ribose zipper is likely to be present in the active, conformationally rearranged ribozyme that is fit for cleavage. We have chemically synthesized a series of individually modified hairpin ribozymes containing 2'-analogues of nucleosides, that include 2'-deoxy and 2'-deoxy-2'-fluoro at each of the four nucleoside positions, 2'-amino-2'-deoxy, 2'-deoxy-2'-thio, and 2'-arabino at position C(25), and 2'-oxyamino at position A(10), as well as some double substitutions, and we studied their cleavage rates under both single- and multiple-turnover conditions. We conclude that at least some of the hydrogen-bonding interactions in the ribose zipper motif, either as originally proposed or in a recently suggested structural variation, are unlikely to be present in the active rearranged form of the ribozyme that undergoes cleavage. Instead, we provide strong evidence for a very precise conformational positioning for the residue C(25) in the active hairpin. A precise conformational requirement would be expected for C(25) if it rearranges to form a base-triple with A(9) and the essential residue neighboring the cleavage site G(+1), as recently proposed by another laboratory. Our results provide further support for conformational rearrangement as an important step in hairpin ribozyme cleavage.

Structure and function of the hairpin ribozyme

The hairpin ribozyme belongs to the family of small catalytic RNAs that cleave RNA substrates in a reversible reaction that generates 2',3'-cyclic phosphate and 5'-hydroxyl termini. The hairpin catalytic motif was discovered in the negative strand of the tobacco ringspot virus satellite RNA, where hairpin ribozyme-mediated self-cleavage and ligation reactions participate in processing RNA replication intermediates. The self-cleaving hairpin, hammerhead, hepatitis delta and Neurospora VS RNAs each adopt unique structures and exploit distinct kinetic and catalytic mechanisms despite catalyzing the same chemical reactions. Mechanistic studies of hairpin ribozyme reactions provided early evidence that, like protein enzymes, RNA enzymes are able to exploit a variety of catalytic strategies. In contrast to the hammerhead and Tetrahymena ribozyme reactions, hairpin-mediated cleavage and ligation proceed through a catalytic mechanism that does not require direct coordination of metal cations to phosphate or water oxygens. The hairpin ribozyme is a better ligase than it is a nuclease while the hammerhead reaction favors cleavage over ligation of bound products by nearly 200-fold. Recent structure-function studies have begun to yield insights into the molecular bases of these unique features of the hairpin ribozyme.

RNA double cleavage by a hairpin-derived twin ribozyme

The hairpin ribozyme is a small catalytic RNA that catalyses reversible sequence-specific RNA hydrolysis in trans. It consists of two domains, which interact with each other by docking in an antiparallel fashion. There is a region between the two domains acting as a flexible hinge for interdomain interactions to occur. Hairpin ribozymes with reverse-joined domains have been constructed by dissecting the domains at the hinge and rejoining them in reverse order. We have used both the conventional and reverse-joined hairpin ribozymes for the design of a hairpin-derived twin ribozyme. We show that this twin ribozyme cleaves a suitable RNA substrate at two specific sites while maintaining the target specificity of the individual monoribozymes. For characterisation of the studied ribozymes we have evaluated a quantitative assay of sequence-specific ribozyme activity using fluorescently labelled RNA substrates in conjunction with an automated DNA sequencer. This assay was found to be applicable with hairpin and hairpin-derived ribozymes. The results demonstrate the potential of hairpin ribozymes for multi-target strategies of RNA cleavage and suggest the possibility for employing hairpin-derived twin ribozymes as powerful tools for RNA manipulation in vitro and in vivo.

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.

Interstrand disulfide cross-linking of internal sugar residues in duplex RNA

Disulfide cross-linking is being used increasingly more to study the structure and dynamics of nucleic acids. We have previously developed a procedure for the formation of disulfide cross-links through the sugar-phosphate backbone of nucleic acids. Here we report the preparation and characterization of an RNA duplex containing a disulfide interstrand cross-link. A self-complementary oligoribonucleotide duplex containing an interstrand cross-link was prepared from the corresponding 2'-amino modified oligomer. Selective modification of the 2'-amino group with an aliphatic isocyanate, containing a protected disulfide, gave the corresponding 2'-urea derivative in excellent yield. An RNA duplex containing an intrahelical, interstrand disulfide cross-link was subsequently prepared by a thiol disulfide exchange reaction in nearly quantitative yield as judged by denaturing polyacrylamide gel electrophoresis (DPAGE). The cross-linked RNA was further characterized by enzymatic digestion and the Structure of the cross-link lesion was verified by comparison to an authentic sample, prepared by chemical synthesis. The effect of the chemical modifications on duplex stability was determined by UV thermal denaturation experiments. The intrahelical cross-link stabilized the duplex considerably: the disulfide cross-linked oligomer had a melting temperature that was ca. 40 degrees C higher than that of the noncross-linked oligomer.

Stitching together RNA tertiary architectures

The powerful explanatory paradigm of molecular biology requiring form to co-evolve with function has again been proven successful when, over the recent two decades, a wealth of biological functions have been uncovered for RNA. Previously considered as a mere mediator of the genetic code, RNA is now acknowledged as a key player in a wide variety of cellular processes. Along with the discovery of novel biological functions of RNA molecules, a number of RNA three-dimensional structures have been solved which beautifully demonstrate the molecular adaptability which allows RNA to participate as a key player in these functions. A distinct repertoire of molecular motifs provides a basis for the assembly of complex RNA tertiary architectures.

Structural basis for the guanosine requirement of the hairpin ribozyme

To form a catalytically active complex, the essential nucleotides of the hairpin ribozyme, embedded within the internal loops of the two domains, must interact with one another. Little is known about the nature of these essential interdomain interactions. In the work presented here, we have used recent topographical constraints and other biochemical data in conjunction with molecular modeling (constraint-satisfaction program MC-SYM) to generate testable models of interdomain interactions. Visual analysis of the generated models has revealed a potential interdomain base pair between the conserved guanosine immediately downstream of the reactive phosphodiester (G(+1)) and C(25) within the large domain. We have tested this former model through activity assays, using all 16 combinations of bases at positions +1 and 25. When the standard ribozyme was used, catalytic activity was severely suppressed with substrates containing U(+1), C(+1), or A(+1). Similarly, mutations of the putative pairing partner (C(25) to A(25) or G(25)) reduce activity by several orders of magnitude. The U(25) substitution retains a significant level of activity, consistent with the possible formation of a G.U wobble pair. Strikingly, when combinations of Watson-Crick (or wobble) base pairs were introduced in these positions, catalytic activity was restored, strongly suggesting the existence of the proposed interaction. These results provide a structural basis for the guanosine requirement of this ribozyme and indicate that the hairpin ribozyme can now be engineered to cleave a wider range of RNA sequences.

Investigation of the recognition of an important uridine in an internal loop of a hairpin ribozyme prepared using post-synthetically modified oligonucleotides

We introduced 4-thio- ((4S)U), 2-thio- ((2S)U), 4- O -methyluridine ((4Me)U) and cytidine substitutions for U+2, which is an important base for cleavage in a substrate RNA. Oligonucleotides containing 4-thio- and 4- O -methyluridine were prepared by a new convenient post-synthetic modification method using a 4- O - p -nitrophenyl-uridine derivative. A hairpin ribozyme cleaved the substrate RNA with either C+2, (4S)U+2 or (4Me)U+2 at approximately 14-, 6- and 4-fold lower rates, respectively, than that of the natural substrate. In contrast, the substrate with a (2S)U+2 was cleaved with the same activity as the natural substrate. These results suggest that the O4 of U+2 is involved in hydrogen bonding at loop A, but the O2 of U+2 is not recognized by the active residues. Circular dichroism data of the ribozymes containing (4S)U+2 and (2S)U+2, as well as the susceptibility of the thiocarbonyl group to hydrogen peroxide, suggest that a conformational change of U+2 occurs during the domain docking in the cleavage reaction. We propose here the conformational change of U+2 from the ground state to the active molecule during the reaction.

Specificity of the hairpin ribozyme. Sequence requirements surrounding the cleavage site

Substrate sequence requirements of the hairpin ribozyme have been partially defined by both mutational and in vitro selection experiments. It was considered that the best targets were those that included the N downward arrowGUC sequence surrounding the cleavage site. In contrast to previous studies that failed to evaluate all possible combinations of these nucleotides, we have performed an exhaustive analysis of the cleavage of 64 substrate variants. They represent all possible sequence combinations of the J2/1 nucleotides except the well established G(+1). No cleavage was observed with 24 sequences. C(+2) variants showed little or no cleavage, whereas U(+2) substrates were all cleavable. The maximal cleavage rate was obtained with the AGUC substrate. Cleavage rates of sequences HGUC (H = A, C, or U), GGUN, GGGR (R = A or G), AGUU, and UGUA were up to 5 times lower than the AGUC one. This shows that other sequences besides NGUC could also be considered as good targets. A second group of sequences WGGG (W = A or U), UGUK (K = G or U), MGAG (M = A or C), AGUA, and UGGA were cleaved between 6 and 10 times less efficiently. Furthermore, the UGCU sequence of a noncleavable viral target was mutated to AGUC resulting in a proficiently cleavable substrate by its cognate hairpin ribozyme. This indicates that our conclusions may be extrapolated to other hairpin ribozymes with different specificity.

Comparative kinetic analysis of structural variants of the hairpin ribozyme reveals further potential to optimize its catalytic performance

The hairpin ribozyme derived from the minus strand of the satellite RNA associated with the tobacco ringspot virus is one of the small catalytic RNAs that has been shown to catalyze trans-cleavage reactions. There is much interest in designing hairpin ribozymes with improved catalytic activity for the development of new therapeutic agents. Extensive mutagenesis studies as well as in vitro selection experiments have been performed to define the structure and optimize its catalytic activity. This communication describes a comparative kinetic analysis of four structural variants, introduced, either alone, or in combination, into the hairpin ribozyme. We have shown that extension of the helix 2 from 4 to 6 bp resulted in a significant decrease in K(M). Furthermore, the combination of this extension with the simultaneous stabilization of helix 4, led to a more than two-fold increase in the catalytic efficiency. This variant showed a 15-fold reduction in the K(M) value in respect to the wild-type ribozyme. This could be of great interest for the in vivo application of this catalytic motif. The 9-bp enlargement of helix 4 implied about a three-fold improvement in the catalytic activity. Similarly, the U39C substitution brought up the efficiency of the ribozyme slightly. However, introduction of nucleotides at the hinge region between A and B domains reduced the catalytic activity. This reduction was gradually increased with the number of nucleotides. Results obtained with variants carrying more than one modification always agreed with the ones obtained from each single variant.