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

Divalent Metal Ion Activation of a Guanine General Base in the Hammerhead Ribozyme: Insights from Molecular Simulations

The hammerhead ribozyme is a well-studied nucleolytic ribozyme that catalyzes the self-cleavage of the RNA phosphodiester backbone. Despite experimental and theoretical efforts, key questions remain about details of the mechanism with regard to the activation of the nucleophile by the putative general base guanine (G12). Straightforward interpretation of the measured activity-pH data implies the pK_a value of the N1 position in the G12 nucleobase is significantly shifted by the ribozyme environment. Recent crystallographic and biochemical work has identified pH-dependent divalent metal ion binding at the N7/O6 position of G12, leading to the hypothesis that this binding mode could induce a pK_a shift of G12 toward neutrality. We present computational results that support this hypothesis and provide a model that unifies the interpretation of available structural and biochemical data and paints a detailed mechanistic picture of the general base step of the reaction. Experimentally testable predictions are made for mutational and rescue effects on G12, which will give further insights into the catalytic mechanism. These results contribute to our growing knowledge of the potential roles of divalent metal ions in RNA catalysis.

Polymerase-Mediated Site-Specific Incorporation of a Synthetic Fluorescent Isomorphic G Surrogate into RNA

An enzyme-mediated approach for the assembly of singly modified RNA constructs in which specific G residues are replaced with th G, an emissive isomorphic G surrogate, is reported. Transcription in the presence of th G and native nucleoside triphosphates enforces initiation with the unnatural analogue, yielding 5'-end modified transcripts that can be mono-phosphorylated and ligated to provide longer site-specifically modified RNA constructs. The scope of this unprecedented enzymatic approach to non-canonical purine-containing RNAs is explored via the assembly of several altered hammerhead (HH) ribozymes and a singly modified HH substrate. By strategically modifying key positions, a mechanistic insight into the ribozyme-mediated cleavage is gained. Additionally, the emissive features of the modified nucleoside and its responsiveness to environmental changes can be used to monitor cleavage in real time by steady state fluorescence spectroscopy.

RNA sociology: group behavioral motifs of RNA consortia

RNA sociology investigates the behavioral motifs of RNA consortia from the social science perspective. Besides the self-folding of RNAs into single stem loop structures, group building of such stem loops results in a variety of essential agents that are highly active in regulatory processes in cellular and non-cellular life. RNA stem loop self-folding and group building do not depend solely on sequence syntax; more important are their contextual (functional) needs. Also, evolutionary processes seem to occur through RNA stem loop consortia that may act as a complement. This means the whole entity functions only if all participating parts are coordinated, although the complementary building parts originally evolved for different functions. If complementary groups, such as rRNAs and tRNAs, are placed together in selective pressure contexts, new evolutionary features may emerge. Evolution initiated by competent agents in natural genome editing clearly contrasts with statistical error replication narratives.

Finding instances of riboswitches and ribozymes by homology search of structured RNA with Infernal

In the genomics era, computational tools are essential to extract information from sequences and annotate them to allow easy access to genes. Fortunately, many of these tools are now part of standard pipelines. As a consequence, a cornucopia of genomic features is available in multiple databases. Nevertheless, as novel genomes are sequenced and new structured RNAs are discovered, homology searches and additional analyses need to be performed. In this chapter, we propose simple ways of finding instances of riboswitches and ribozymes in databases or in unannotated genomes, as well as ways of finding variants that deviate from the typical consensus.

Characterization of hairpin ribozyme reactions

Hairpin ribozymes are small RNA catalytic motifs naturally found in the satellite RNAs of tobacco ringspot virus (TRsV), chicory yellow mottle virus (CYMoV), and arabis mosaic virus (ArMV). The catalytic activity of the hairpin ribozyme extends to both cleavage and ligation reactions. Here we describe methods for the kinetic analysis of the self-cleavage reaction under transcription reaction conditions. We also describe methods for the generation of DNA templates for subsequent in vitro transcription reaction of hairpin ribozymes. This is followed by a description of the preparation of the suitable RNA molecules for ligation reaction and their kinetic analysis.

G17-modified hammerhead ribozymes are active in vitro and in vivo

Natural hammerhead ribozymes (HHRz) feature tertiary interactions between hairpin loops or bulges in two of three helices that surround the catalytic core of conserved nucleotides. Their conservation was originally established on minimal versions lacking the tertiary interactions. While those sequence requirements in general also apply to natural versions, we show here differences for the HHRz cleavage site N17. A guanosine at this position strongly impairs cleavage activity in minimal versions, whereas we observe for the G17 variants of four tertiary stabilized HHRz significant cleavage and ligation activity in vitro. Kinetic analyses of these variants revealed a reduced rate and extent of cleavage, compared with wild-type sequences, while variants with distorted tertiary interactions cleaved at a reduced rate, but to the same extent. Contrary to this, G17 variants exhibit similar in vitro ligation activity as compared with the respective wild-type motif. To also address the catalytic performance of these motifs in vivo, we have inserted HHRz cassettes in the lacZ gene and tested this β-galactosidase reporter in Dictyostelium discoideum. In colorimetric assays, we observe differences in the enzymatic activity of β-galactosidase, which correlate well with the activity of the different HHRz variants in vitro and which can be unambiguously attributed to ribozyme cleavage by primer extension analysis.

Prediction of hammerhead ribozyme intracellular activity with the catalytic core fingerprint

Hammerhead ribozyme is a versatile tool for down-regulation of gene expression in vivo. Owing to its small size and high activity, it is used as a model for RNA structure–function relationship studies. In the present paper we describe a new extended hammerhead ribozyme HH-2 with a tertiary stabilizing motif constructed on the basis of the tetraloop receptor sequence. This ribozyme is very active in living cells, but shows low activity in vitro. To understand it, we analysed tertiary structure models of substrate–ribozyme complexes. We calculated six unique catalytic core geometry parameters as distances and angles between particular atoms that we call the ribozyme fingerprint. A flanking sequence and tertiary motif change the geometry of the general base, general acid, nucleophile and leaving group. We found almost complete correlation between these parameters and the decrease of target gene expression in the cells. The tertiary structure model calculations allow us to predict ribozyme intracellular activity. Our approach could be widely adapted to characterize catalytic properties of other RNAs.

The DNA Habitat and its RNA Inhabitants: At the Dawn of RNA Sociology

Most molecular biological concepts derive from physical chemical assumptions about the genetic code that are basically more than 40 years old. Additionally, systems biology, another quantitative approach, investigates the sum of interrelations to obtain a more holistic picture of nucleotide sequence order. Recent empirical data on genetic code compositions and rearrangements by mobile genetic elements and noncoding RNAs, together with results of virus research and their role in evolution, does not really fit into these concepts and compel a reexamination. In this review, we try to find an alternate hypothesis. It seems plausible now that if we look at the abundance of regulatory RNAs and persistent viruses in host genomes, we will find more and more evidence that the key players that edit the genetic codes of host genomes are consortia of RNA agents and viruses that drive evolutionary novelty and regulation of cellular processes in all steps of development. This agent-based approach may lead to a qualitative RNA sociology that investigates and identifies relevant behavioral motifs of cooperative RNA consortia. In addition to molecular biological perspectives, this may lead to a better understanding of genetic code evolution and dynamics.

Bridging the gap between theory and experiment to derive a detailed understanding of hammerhead ribozyme catalysis

Herein we summarize our progress toward the understanding of hammerhead ribozyme (HHR) catalysis through a multiscale simulation strategy. Simulation results collectively paint a picture of HHR catalysis: HHR first folds to form an electronegative active site pocket to recruit a threshold occupation of cationic charges, either a Mg(2+) ion or multiple monovalent cations. Catalytically active conformations that have good in-line fitness are supported by specific metal ion coordination patterns that involve either a bridging Mg(2+) ion or multiple Na(+) ions, one of which is also in a bridging coordination pattern. In the case of a single Mg(2+) ion bound in the active site, the Mg(2+) ion undergoes a migration that is coupled with deprotonation of the nucleophile (C17:O2'). As the reaction proceeds, the Mg(2+) ion stabilizes the accumulating charge of the leaving group and significantly increases the general acid ability of G8:O2'. Further computational mutagenesis simulations suggest that the disruptions due to mutations may severely impact HHR catalysis at different stages of the reaction. Catalytic mechanisms supported by the simulation results are consistent with available structural and biochemical experiments, and together they advance our understanding of HHR catalysis.

Viroid-specific small RNA in plant disease

Viroids are the smallest autonomous infectious nucleic acids known today. They are non-coding, unencapsidated, circular RNAs with sizes ranging from 250 to 400 nucleotides and infect certain plants. These RNAs are transcribed by rolling-circle mechanisms in the plant host's nuclei (Pospiviroidae) or chloroplasts (Avsunviroidae). Since viroids lack any open reading frame, their pathogenicity has for a long time been a conundrum. Recent findings, however, show that viroid infection is associated with the appearance of viroid-specific small RNA (vsRNA). These have sizes similar to endogenous small interfering RNA and microRNA and thus might alter the normal gene expression in the host plant. In this review we will summarize the current knowledge on vsRNA and discuss the current hypotheses how they connect to the induced symptoms, which vary dramatically, depending on both the plant cultivar and the viroid strain.

The ubiquitous hammerhead ribozyme

The hammerhead ribozyme is a small catalytic RNA motif capable of endonucleolytic (self-) cleavage. It is composed of a catalytic core of conserved nucleotides flanked by three helices, two of which form essential tertiary interactions for fast self-scission under physiological conditions. Originally discovered in subviral plant pathogens, its presence in several eukaryotic genomes has been reported since. More recently, this catalytic RNA motif has been shown to reside in a large number of genomes. We review the different approaches in discovering these new hammerhead ribozyme sequences and discuss possible biological functions of the genomic motifs.

Folding of the hammerhead ribozyme: pyrrolo-cytosine fluorescence separates core folding from global folding and reveals a pH-dependent conformational change

The catalytic activity of the hammerhead ribozyme is limited by its ability to fold into the native tertiary structure. Analysis of folding has been hampered by a lack of assays that can independently monitor the environment of nucleobases throughout the ribozyme-substrate complex in real time. Here, we report the development and application of a new folding assay in which we use pyrrolo-cytosine (pyC) fluorescence to (1) probe active-site formation, (2) examine the ability of peripheral ribozyme domains to support native folding, (3) identify a pH-dependent conformational change within the ribozyme, and (4) explore its influence on the equilibrium between the folded and unfolded core of the hammerhead ribozyme. We conclude that the natural ribozyme folds in two distinct noncooperative steps and the pH-dependent correlation between core folding and activity is linked to formation of the G8-C3 base pair.

Characterization of hammerhead ribozyme reactions

Hammerhead ribozymes are small catalytic RNA motifs ubiquitously present in a large variety of genomes. The reactions catalyzed by these motifs are both their self-scission and the reverse ligation reaction. Here, we describe methods for the generation of DNA templates for the subsequent in vitro transcription of hammerhead ribozymes. This is followed by a description of the preparation of suitable RNA molecules for both reaction types, and their kinetic analysis.

Sequence elements outside the catalytic core of natural hairpin ribozymes modulate the reactions differentially

Abstract Hairpin ribozymes occur naturally only in the satellite RNAs of tobacco ringspot virus (TRsV), chicory yellow mottle virus (CYMoV) and arabis mosaic virus (ArMV). The catalytic centre of the predominantly studied sTRsV hairpin ribozyme, and of sArMV is organised around a four-way helical junction. We show here that sCYMoV features a five-way helical junction instead. Mutational analysis indicates that the fifth stem does not influence kinetic parameters of the sCYMoV hairpin ribozyme in vitro reactions, and therefore seems an appendix to that junction in the other ribozymes. We report further that all three ribozymes feature a three-way helical junction outside the catalytic core in stem A, with Watson-Crick complementarity to loop nucleotides in stem B. Kinetic analyses of cleavage and ligation reactions of several variants of the sTRsV and sCYMoV hairpin ribozymes in vitro show that the presence of this junction interferes with their reactions, particularly the ligation. We provide evidence that this is not due to a presumed interaction of the afore-mentioned elements in stems A and B. The evolutionary survival of this cis-inhibiting element seems rather to be caused by the coincidence of its position with that of the hammerhead ribozyme in the other RNA polarity.

Identification of hammerhead ribozymes in all domains of life reveals novel structural variations

Hammerhead ribozymes are small self-cleaving RNAs that promote strand scission by internal phosphoester transfer. Comparative sequence analysis was used to identify numerous additional representatives of this ribozyme class than were previously known, including the first representatives in fungi and archaea. Moreover, we have uncovered the first natural examples of "type II" hammerheads, and our findings reveal that this permuted form occurs in bacteria as frequently as type I and III architectures. We also identified a commonly occurring pseudoknot that forms a tertiary interaction critical for high-speed ribozyme activity. Genomic contexts of many hammerhead ribozymes indicate that they perform biological functions different from their known role in generating unit-length RNA transcripts of multimeric viroid and satellite virus genomes. In rare instances, nucleotide variation occurs at positions within the catalytic core that are otherwise strictly conserved, suggesting that core mutations are occasionally tolerated or preferred.

From alpaca to zebrafish: hammerhead ribozymes wherever you look

The hammerhead ribozyme was originally discovered in subviral plant pathogens and was subsequently also found in a few other genomic locations. Using a secondary structure-based descriptor, we have searched publicly accessible sequence databases for new examples of type III hammerhead ribozymes. The more than 60,000 entries fulfilling the descriptor were filtered with respect to folding and stability parameters that were experimentally validated. This resulted in a set of 284 unique motifs, of which 124 represent database entries of known hammerhead ribozymes from subviral plant pathogens and A. thaliana. The remainder are 160 novel ribozyme candidates in 50 different eukaryotic genomes. With a few exceptions, the ribozymes were found either in repetitive DNA sequences or in introns of protein coding genes. Our data, which is complementary to a study by De la Peña and García-Robles in 2010, indicate that the hammerhead is the most abundant small endonucleolytic ribozyme, which, in view of no sequence conservation beyond the essential nucleotides, likely has evolved independently in different organisms.

Computational mutagenesis studies of hammerhead ribozyme catalysis

Computational studies of the mutational effects at the C3, G8, and G5 positions of the hammerhead ribozyme (HHR) are reported, based on a series of twenty-four 100-ns molecular dynamics simulations of the native and mutated HHR in the reactant state and in an activated precursor state (G8:2'OH deprotonated). Invoking the assumptions that G12 acts as the general base while the 2'OH of G8 acts as a general acid, the simulations are able to explain the origins of experimentally observed mutational effects, including several that are not easily inferred from the crystal structure. Simulations suggest that the Watson-Crick base-pairing between G8 and C3, the hydrogen bond network between C17 and G5, and the base stacking interactions between G8 and C1.1, collectively, are key to maintaining an active site structure conducive for catalytic activity. Mutation-induced disruption of any of these interactions will adversely affect activity. The simulation results predict that the C3U/G8D double mutant, where D is 2,6-diaminopurine, will have a rescue effect relative to the corresponding single mutations. Two general conclusions about the simulations emerge from this work. First, mutation simulations may require 30 ns or more to suitably relax such that the mutational effects become apparent. Second, in some cases, it is necessary to look beyond the reactant state in order to interpret mutational effects in terms of catalytically active structure. The present simulation results lead to better understanding of the origin of experimental mutational effects and provide insight into the key conserved features necessary to maintain the integrity of the active site architecture.

Phylogenetic footprinting of non-coding RNA: hammerhead ribozyme sequences in a satellite DNA family of Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae)

BACKGROUND

The great variety in sequence, length, complexity, and abundance of satellite DNA has made it difficult to ascribe any function to this genome component. Recent studies have shown that satellite DNA can be transcribed and be involved in regulation of chromatin structure and gene expression. Some satellite DNAs, such as the pDo500 sequence family in Dolichopoda cave crickets, have a catalytic hammerhead (HH) ribozyme structure and activity embedded within each repeat.

RESULTS

We assessed the phylogenetic footprints of the HH ribozyme within the pDo500 sequences from 38 different populations representing 12 species of Dolichopoda. The HH region was significantly more conserved than the non-hammerhead (NHH) region of the pDo500 repeat. In addition, stems were more conserved than loops. In stems, several compensatory mutations were detected that maintain base pairing. The core region of the HH ribozyme was affected by very few nucleotide substitutions and the cleavage position was altered only once among 198 sequences. RNA folding of the HH sequences revealed that a potentially active HH ribozyme can be found in most of the Dolichopoda populations and species.

CONCLUSIONS

The phylogenetic footprints suggest that the HH region of the pDo500 sequence family is selected for function in Dolichopoda cave crickets. However, the functional role of HH ribozymes in eukaryotic organisms is unclear. The possible functions have been related to trans cleavage of an RNA target by a ribonucleoprotein and regulation of gene expression. Whether the HH ribozyme in Dolichopoda is involved in similar functions remains to be investigated. Future studies need to demonstrate how the observed nucleotide changes and evolutionary constraint have affected the catalytic efficiency of the hammerhead.

Comparative enzymology and structural biology of RNA self-cleavage

Self-cleaving hammerhead, hairpin, hepatitis delta virus, and glmS ribozymes comprise a family of small catalytic RNA motifs that catalyze the same reversible phosphodiester cleavage reaction, but each motif adopts a unique structure and displays a unique array of biochemical properties. Recent structural, biochemical, and biophysical studies of these self-cleaving RNAs have begun to reveal how active site nucleotides exploit general acid-base catalysis, electrostatic stabilization, substrate destabilization, and positioning and orientation to reduce the free energy barrier to catalysis. Insights into the variety of catalytic strategies available to these model RNA enzymes are likely to have important implications for understanding more complex RNA-catalyzed reactions fundamental to RNA processing and protein synthesis.

The unforeseeable hammerhead ribozyme

Despite its small size, the complex behavior of the hammerhead ribozyme keeps surprising us, even more than 20 years after its discovery. Here, we summarize recent developments in the field, in particular the discovery of the first split hammerhead ribozyme.

Catalytic strategies of self-cleaving ribozymes

[Structure: see text]. Five naturally occurring nucleolytic ribozymes have been identified: the hammerhead, hairpin, glmS, hepatitis delta virus (HDV), and Varkud satellite (VS) ribozymes. All of these RNA enzymes catalyze self-scission of the RNA backbone using a chemical mechanism equivalent to that of RNase A. RNase A uses four basic strategies to promote this reaction: geometric constraints, activation of the nucleophile, transition-state stabilization, and leaving group protonation. In this Account, we discuss the current thinking on how nucleolytic ribozymes harness RNase A's four sources of catalytic power. The geometry of the phosphodiester cleavage reaction constrains the nucleotides flanking the scissile phosphate so that they are unstacked from a canonical A-form helix and thus require alternative stabilization. Crystal structures and mutational analysis reveal that cross-strand base pairing, along with unconventional stacking and tertiary hydrogen-bonding interactions, work to stabilize the splayed conformation in nucleolytic ribozymes. Deprotonation of the 2'-OH nucleophile greatly increases its nucleophilicity in the strand scission reaction. Crystal structures of the hammerhead, hairpin, and glmS ribozymes reveal the N1 of a G residue within hydrogen-bonding distance of the 2'-OH. In each case, this residue has also been shown to be important for catalysis. In the HDV ribozyme, a hydrated magnesium has been implicated as the general base. Catalysis by the VS ribozyme requires both an A and a G, but the precise role of either has not been elucidated. Enzymes can lower the energy of a chemical reaction by binding more tightly to the transition state than to the ground states. Comparison of the hairpin ground- and transition-state mimic structures reveal greater hydrogen bonding to the transition-state mimic structure, suggesting transition-state stabilization as a possible catalytic strategy. However, the hydrogen-bonding pattern in the glmS ribozyme transition-state mimic structure and the ground-state structures are equivalent. Protonation of the 5'-O leaving group by a variety of functional groups can promote the cleavage reaction. In the HDV ribozyme, the general acid is a conserved C residue. In the hairpin ribozyme, a G residue has been implicated in protonation of the leaving group. An A in the hammerhead ribozyme probably plays a similar role. In the glmS ribozyme, an exogenous cofactor may provide the general acid. This diversity is in contrast to the relatively small number of functional groups that serve as a general base, where at least three of the nucleolytic ribozymes may use the N1 of a G.

Solvent structure and hammerhead ribozyme catalysis

Although the hammerhead ribozyme is regarded as a prototype for understanding RNA catalysis, the mechanistic roles of associated metal ions and water molecules in the cleavage reaction remain controversial. We have investigated the catalytic potential of observed divalent metal ions and water molecules bound to a 2 A structure of the full-length hammerhead ribozyme by using X-ray crystallography in combination with molecular dynamics simulations. A single Mn(2+) is observed to bind directly to the A9 phosphate in the active site, accompanying a hydrogen-bond network involving a well-ordered water molecule spanning N1 of G12 (the general base) and 2'-O of G8 (previously implicated in general acid catalysis) that we propose, based on molecular dynamics calculations, facilitates proton transfer in the cleavage reaction. Phosphate-bridging metal interactions and other mechanistic hypotheses are also tested with this approach.

Origin of mutational effects at the C3 and G8 positions on hammerhead ribozyme catalysis from molecular dynamics simulations

A series of ten 60 ns molecular dynamics (MD) simulations of the native and mutated full length hammerhead ribozymes in the reactant state and in an activated precursor state (G8:2'OH deprotonated) are reported. Mutant simulations include the C3U, G8A, and G8I single mutants and a C3U/G8A double mutant that exhibits an experimental rescue effect. The results provide critical details into the origin of the observed mutation effects and support a mechanism where the 2'OH of G8 acts as a general acid catalyst that is held in position through Watson-Crick hydrogen bonding between G8 and C3.

Hammerhead redux: does the new structure fit the old biochemical data?

The cleavage rates of 78 hammerhead ribozymes containing structurally conservative chemical modifications were collected from the literature and compared to the recently determined crystal structure of the Schistosoma mansoni hammerhead. With only a few exceptions, the biochemical data were consistent with the structure, indicating that the new structure closely resembles the transition state of the reaction. Since all the biochemical data were collected on minimal hammerheads that have a very different structure, the minimal hammerhead must be dynamic and occasionally adopt the quite different extended structure in order to cleave.

Minimal and extended hammerheads utilize a similar dynamic reaction mechanism for catalysis

Analysis of the catalytic activity of identical mutations in the catalytic cores of nHH8, a very active "extended" hammerhead, and HH16, a less active "minimal" hammerhead, reveal that the tertiary Watson-Crick base pair between C3 and G8 seen in the recent structure of the Schistosoma mansoni extended hammerhead can be replaced by other base pairs in both backgrounds. This supports the model that both hammerheads utilize a similar catalytic mechanism but HH16 is slower because it infrequently samples the active conformation. The relative effect of different mutations at positions 3 and 8 also depends on the identity of residue 17 in both nHH8 and HH16. This synergistic effect can best be explained by transient pairing between residues 3 and 17 and 8 and 13, which stabilize an inactive conformation. Thus, mutants of nHH8 and possibly nHH8 itself are also in dynamic equilibrium with an inactive conformation that may resemble the X-ray structure of a minimal hammerhead. Therefore, both minimal and extended hammerhead structures must be considered to fully understand hammerhead catalysis.