General acid catalysis by the hepatitis delta virus ribozyme

Chemical Synthesis of Modified RNA

Ribonucleic acids (RNAs) play a vital role in living organisms. Many of their cellular functions depend critically on chemical modification. Methods to modify RNA in a controlled manner-both in vitro and in vivo-are thus essential to evaluate and understand RNA biology at the molecular and mechanistic levels. The diversity of modifications, combined with the size and uniformity of RNA (made up of only 4 nucleotides) makes its site-specific modification a challenging task that needs to be addressed by complementary approaches. One such approach is solid-phase RNA synthesis. We discuss recent developments in this field, starting with new protection concepts in the ongoing effort to overcome current size limitations. We continue with selected modifications that have posed significant challenges for their incorporation into RNA. These include deazapurine bases required for atomic mutagenesis to elucidate mechanistic aspects of catalytic RNAs, and RNA containing xanthosine, N4-acetylcytidine, 5-hydroxymethylcytidine, 3-methylcytidine, 2'-OCF3, and 2'-N3 ribose modifications. We also discuss the all-chemical synthesis of 5'-capped mRNAs and the enzymatic ligation of chemically synthesized oligoribonucleotides to obtain long RNA with multiple distinct modifications, such as those needed for single-molecule FRET studies. Finally, we highlight promising developments in RNA-catalyzed RNA modification using cofactors that transfer bioorthogonal functionalities.

New Insights into the Dependence of CPEB3 Ribozyme Cleavage on Mn2+ and Mg2

CPEB3 ribozyme is a self-cleaving RNA that occurs naturally in mammals and requires divalent metal ions for efficient activity. Ribozymes exhibit preferences for specific metal ions, but the exact differences in the catalytic mechanisms of various metal ions on the CPEB3 ribozyme remain unclear. Our findings reveal that Mn2+ functions as a more effective cofactor for CPEB3 ribozyme catalysis compared to Mg2+, as confirmed by its stronger binding affinity to CPEB3 by EPR. Cleavage assays of CPEB3 mutants and molecular docking analyses further showed that excessive Mn2+ ions can bind to a second binding site near the catalytic site, hindering CPEB3 catalytic efficiency and contributing to the Mn2+ bell-shaped curve. These results implicate a pivotal role for the local nucleobase-Mn2+ interactions in facilitating RNA folding and modulating the directed attack of nucleophilic reagents. Our study provides new insights and experimental evidence for exploring the divalent cation dependent cleavage mechanism of the CPEB3 ribozyme.

Inhibition of Cpeb3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs and enhances object location memory

A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element-binding protein 3 (Cpeb3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the Cpeb3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.

Inhibition of CPEB3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs, and enhances object location memory

A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element binding protein 3 (CPEB3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the CPEB3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.

Dissociative Transition State in Hepatitis Delta Virus Ribozyme Catalysis

Ribonucleases and small nucleolytic ribozymes are both able to catalyze RNA strand cleavage through 2'-O-transphosphorylation, provoking the question of whether protein and RNA enzymes facilitate mechanisms that pass through the same or distinct transition states. Here, we report the primary and secondary 18O kinetic isotope effects for hepatitis delta virus ribozyme catalysis that reveal a dissociative, metaphosphate-like transition state in stark contrast to the late, associative transition states observed for reactions catalyzed by specific base, Zn2+ ions, or ribonuclease A. This new information provides evidence for a discrete ribozyme active site design that modulates the RNA cleavage pathway to pass through an altered transition state.

A Cofactor-Based Mechanism for the Origin of the Genetic Code

The origin of the genetic code is probably the central problem of the studies on the origin of life. The key question to answer is the molecular mechanism that allows the association of the amino acids with their triplet codons. We proposed that the codon-anticodon duplex located in the acceptor stem of primitive tRNAs would facilitate the chemical reactions required to synthesize cognate amino acids from simple amino acids (glycine, valine, and aspartic acid) linked to the 3' acceptor end. In our view, various nucleotide-A-derived cofactors (with reactive chemical groups) may be attached to the codon-anticodon duplex, which allows group-transferring reactions from cofactors to simple amino acids, thereby producing the final amino acid. The nucleotide-A-derived cofactors could be incorporated into the RNA duplex (helix) by docking Adenosine (cofactor) into the minor groove via an interaction similar to the A-minor motif, forming a base triple between Adenosine and one complementary base pair of the duplex. Furthermore, we propose that this codon-anticodon duplex could initially catalyze a self-aminoacylation reaction with a simple amino acid. Therefore, the sequence of bases in the codon-anticodon duplex would determine the reactions that occurred during the formation of new amino acids for selective binding of nucleotide-A-derived cofactors.

1-Deazaguanosine-Modified RNA: The Missing Piece for Functional RNA Atomic Mutagenesis

Atomic mutagenesis is the key to advance our understanding of RNA recognition and RNA catalysis. To this end, deazanucleosides are utilized to evaluate the participation of specific atoms in these processes. One of the remaining challenges is access to RNA-containing 1-deazaguanosine (c¹G). Here, we present the synthesis of this nucleoside and its phosphoramidite, allowing first time access to c¹G-modified RNA. Thermodynamic analyses revealed the base pairing parameters for c¹G-modified RNA. Furthermore, by NMR spectroscopy, a c¹G-triggered switch of Watson-Crick into Hoogsteen pairing in HIV-2 TAR RNA was identified. Additionally, using X-ray structure analysis, a guanine–phosphate backbone interaction affecting RNA fold stability was characterized, and finally, the critical impact of an active-site guanine in twister ribozyme on the phosphodiester cleavage was revealed. Taken together, our study lays the synthetic basis for c¹G-modified RNA and demonstrates the power of the completed deazanucleoside toolbox for RNA atomic mutagenesis needed to achieve in-depth understanding of RNA recognition and catalysis.

Structural basis for substrate binding and catalysis by a self-alkylating ribozyme

Ribozymes that react with small-molecule probes have important applications in transcriptomics and chemical biology, such as RNA labeling and imaging. Understanding the structural basis for these RNA-modifying reactions will enable the development of better tools for studying RNA. Nevertheless, high-resolution structures and underlying catalytic mechanisms for members of this ribozyme class remain elusive. Here, we focus on a self-alkylating ribozyme that catalyzes nitrogen-carbon bond formation between a specific guanine and a 2,3-disubstituted epoxide substrate and report the crystal structures of a self-alkylating ribozyme, including both alkylated and apo forms, at 1.71-Å and 2.49-Å resolution, respectively. The ribozyme assumes an elongated hairpin-like architecture preorganized to accommodate the epoxide substrate in a hook-shaped conformation. Observed reactivity of substrate analogs together with an inverse, log-linear pH dependence of the reaction rate suggests a requirement for epoxide protonation, possibly assisted by the ether oxygens within the substrate.

Self-cleaving ribozymes: substrate specificity and synthetic biology applications

Various self-cleaving ribozymes appearing in nature catalyze the sequence-specific intramolecular cleavage of RNA and can be engineered to catalyze cleavage of appropriate substrates in an intermolecular fashion, thus acting as true catalysts. The mechanisms of the small, self-cleaving ribozymes have been extensively studied and reviewed previously. Self-cleaving ribozymes can possess high catalytic activity and high substrate specificity; however, substrate specificity is also engineerable within the constraints of the ribozyme structure. While these ribozymes share a common fundamental catalytic mechanism, each ribozyme family has a unique overall architecture and active site organization, indicating that several distinct structures yield this chemical activity. The multitude of catalytic structures, combined with some flexibility in substrate specificity within each family, suggests that such catalytic RNAs, taken together, could access a wide variety of substrates. Here, we give an overview of 10 classes of self-cleaving ribozymes and capture what is understood about their substrate specificity and synthetic applications. Evolution of these ribozymes in an RNA world might be characterized by the emergence of a new ribozyme family followed by rapid adaptation or diversification for specific substrates.

Self-cleaving ribozymes have become important tools of synthetic biology. Here we summarize the substrate specificity and applications of the main classes of these ribozymes.

Synthesis of Oligoribonucleotides Containing a 2'-Amino-5'-S-phosphorothiolate Linkage

Oligoribonucleotides containing a photocaged 2′-amino-5′-S-phophorothiolate linkage have potential applications as therapeutic agents and biological probes to investigate the RNA structure and function. We envisioned that oligoribonucleotides containing a 2′-amino-5′-S-phosphorothiolate linkage could provide an approach to identify the general base within catalytic RNAs by chemogenetic suppression. To enable preliminary tests of this idea, we developed synthetic approaches to a dinucleotide, trinucleotide, and oligoribonucleotide containing a photocaged 2′-amino-5′-S-phosphorothiolate linkage. We incorporated the photocaged 2′-amino-5′-S-phosphorothiolate linkage into an oligoribonucleotide substrate for the hepatitis delta virus (HDV) ribozyme and investigated the pH dependence of its cleavage following UV irradiation both in the presence and absence of the ribozyme. The substrate exhibited a pH-rate profile characteristic of the modified linkage but reacted slower when bound to the ribozyme. Cleavage inhibition by the HDV ribozyme could reflect a non-productive ground-state interaction with the modified substrate’s nucleophilic 2′-NH₂ or a poor fit of the modified transition state at the ribozyme’s active site.

Beyond the Plateau: pL Dependence of Proton Inventories as a Tool for Studying Ribozyme and Ribonuclease Catalysis

Acid/base catalysis is an important catalytic strategy used by ribonucleases and ribozymes; however, understanding the number and identity of functional groups involved in proton transfer remains challenging. The proton inventory (PI) technique analyzes the dependence of the enzyme reaction rate on the ratio of D₂O to H₂O and can provide information about the number of exchangeable sites that produce isotope effects and their magnitude. The Gross-Butler (GB) equation is used to evaluate H/D fractionation factors from PI data typically collected under conditions (i.e., a "plateau" in the pH-rate profile) assuming minimal change in active site residue ionization. However, restricting PI analysis to these conditions is problematic for many ribonucleases, ribozymes, and their variants due to ambiguity in the roles of active site residues, the lack of a plateau within the accessible pL range, or cooperative interactions between active site functional groups undergoing ionization. Here, we extend the integration of species distributions for alternative enzyme states in noncooperative models of acid/base catalysis into the GB equation, first used by Bevilacqua and colleagues for the HDV ribozyme, to develop a general population-weighted GB equation that allows simulation and global fitting of the three-dimensional relationship of the D₂O ratio (n) versus pL versus k_n/k₀. Simulations using the GPW-GB equation of PI results for RNase A, HDVrz, and VSrz illustrate that data obtained at multiple selected pL values across the pL-rate profile can assist in the planning and interpreting of solvent isotope effect experiments to distinguish alternative mechanistic models.

The potential versatility of RNA catalysis

It is commonly thought that in the early development of life on this planet RNA would have acted both as a store of genetic information and as a catalyst. While a number of RNA enzymes are known in contemporary cells, they are largely confined to phosphoryl transfer reactions, whereas an RNA based metabolism would have required a much greater chemical diversity of catalysis. Here we discuss how RNA might catalyze a wider variety of chemistries, and particularly how information gleaned from riboswitches could suggest how ribozymes might recruit coenzymes to expand their chemical range. We ask how we might seek such activities in modern biology. This article is categorized under: RNA-Based Catalysis > Miscellaneous RNA-Catalyzed Reactions Regulatory RNAs/RNAi/Riboswitches > Riboswitches RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Impact of 3-deazapurine nucleobases on RNA properties

Deazapurine nucleosides such as 3-deazaadenosine (c3A) are crucial for atomic mutagenesis studies of functional RNAs. They were the key for our current mechanistic understanding of ribosomal peptide bond formation and of phosphodiester cleavage in recently discovered small ribozymes, such as twister and pistol RNAs. Here, we present a comprehensive study on the impact of c3A and the thus far underinvestigated 3-deazaguanosine (c3G) on RNA properties. We found that these nucleosides can decrease thermodynamic stability of base pairing to a significant extent. The effects are much more pronounced for 3-deazapurine nucleosides compared to their constitutional isomers of 7-deazapurine nucleosides (c7G, c7A). We furthermore investigated base pair opening dynamics by solution NMR spectroscopy and revealed significantly enhanced imino proton exchange rates. Additionally, we solved the X-ray structure of a c3A-modified RNA and visualized the hydration pattern of the minor groove. Importantly, the characteristic water molecule that is hydrogen-bonded to the purine N3 atom and always observed in a natural double helix is lacking in the 3-deazapurine-modified counterpart. Both, the findings by NMR and X-ray crystallographic methods hence provide a rationale for the reduced pairing strength. Taken together, our comparative study is a first major step towards a comprehensive understanding of this important class of nucleoside modifications.

Zn2+-dependent DNAzymes that cleave all combinations of ribonucleotides

Although several DNAzymes are known, their utility is limited by a narrow range of substrate specificity. Here, we report the isolation of two zinc-dependent DNAzymes, ZincDz1 and ZincDz2, which exhibit compact catalytic core sequences with highly versatile hydrolysis activity. They were selected through in vitro selection followed by deep sequencing analysis. Despite their sequence similarity, each DNAzyme showed different Zn²⁺-concentration and pH-dependent reaction profiles, and cleaved the target RNA sequences at different sites. Using various substrate RNA sequences, we found that the cleavage sequence specificity of ZincDz2 and its highly active mutant ZincDz2-v2 to be 5'-rN↓rNrPu-3'. Furthermore, we demonstrated that the designed ZincDz2 could cut microRNA miR-155 at three different sites. These DNAzymes could be useful in a broad range of applications in the fields of medicine and biotechnology.

Single-pass transcription by T7 RNA polymerase

RNA molecules can be conveniently synthesized in vitro by the T7 RNA polymerase (T7 RNAP). In some experiments, such as cotranscriptional biochemical analyses, continuous synthesis of RNA is not desired. Here, we propose a method for a single-pass transcription that yields a single transcript per template DNA molecule using the T7 RNAP system. We hypothesized that stalling the polymerase downstream from the promoter region and subsequent cleavage of the promoter by a restriction enzyme (to prevent promoter binding by another polymerase) would allow synchronized production of a single transcript per template. The single-pass transcription was verified in two different scenarios: a short self-cleaving ribozyme and a long mRNA. The results show that a controlled single-pass transcription using T7 RNAP allows precise measurement of cotranscriptional ribozyme activity, and this approach will facilitate the study of other kinetic events.

Fundamental studies of functional nucleic acids: aptamers, riboswitches, ribozymes and DNAzymes

This review aims at juxtaposing common versus distinct structural and functional strategies that are applied by aptamers, riboswitches, and ribozymes/DNAzymes. Focusing on recently discovered systems, we begin our analysis with small-molecule binding aptamers, with emphasis on in vitro-selected fluorogenic RNA aptamers and their different modes of ligand binding and fluorescence activation. Fundamental insights are much needed to advance RNA imaging probes for detection of exo- and endogenous RNA and for RNA process tracking. Secondly, we discuss the latest gene expression-regulating mRNA riboswitches that respond to the alarmone ppGpp, to PRPP, to NAD⁺, to adenosine and cytidine diphosphates, and to precursors of thiamine biosynthesis (HMP-PP), and we outline new subclasses of SAM and tetrahydrofolate-binding RNA regulators. Many riboswitches bind protein enzyme cofactors that, in principle, can catalyse a chemical reaction. For RNA, however, only one system (glmS ribozyme) has been identified in Nature thus far that utilizes a small molecule - glucosamine-6-phosphate - to participate directly in reaction catalysis (phosphodiester cleavage). We wonder why that is the case and what is to be done to reveal such likely existing cellular activities that could be more diverse than currently imagined. Thirdly, this brings us to the four latest small nucleolytic ribozymes termed twister, twister-sister, pistol, and hatchet as well as to in vitro selected DNA and RNA enzymes that promote new chemistry, mainly by exploiting their ability for RNA labelling and nucleoside modification recognition. Enormous progress in understanding the strategies of nucleic acids catalysts has been made by providing thorough structural fundaments (e.g. first structure of a DNAzyme, structures of ribozyme transition state mimics) in combination with functional assays and atomic mutagenesis.

The Positively Charged Active Site of the Bacterial Toxin RelE Causes a Large Shift in the General Base pKa

The bacterial toxin RelE cleaves mRNA in the ribosomal A site. Although it shares a global fold with other microbial RNases, the active site contains several positively charged residues instead of histidines and glutamates that are typical of ribonucleases. The pH dependences of wild-type and mutant RelE indicate it uses general acid-base catalysis, but either the general acid (proposed to be R81) or the general base must have a substantially downshifted pK_a. However, which group is shifted cannot be determined using available structural and biochemical data. Here, we use a phosphorothiolate at the scissile phosphate to remove the need for a general acid. We show this modification rescues nearly all of the defect of the R81A mutation, supporting R81 as the general acid. We also find that the observed pK_a of the general base is dependent on the charge of the side chain at position 81. This indicates that positive charge in the active site contributes to a general base pK_a downshifted by more than 5 units. Although this modestly reduces the effectiveness of general acid-base catalysis, it is strongly supplemented by the role of the positive charge in stabilizing the transition state for cleavage. Furthermore, we show that the ribosome is required for cleavage but not binding of mRNA by RelE. Ribosome functional groups do not directly contact the scissile phosphate, indicating that positioning and charge interactions dominate RelE catalysis. The unusual RelE active site catalyzes phosphoryl transfer at a rate comparable to those of similar enzymes, but in a ribosome-dependent fashion.

In vitro and in vivo properties of therapeutic oligonucleotides containing non-chiral 3' and 5' thiophosphate linkages

The introduction of non-bridging phosphorothioate (PS) linkages in oligonucleotides has been instrumental for the development of RNA therapeutics and antisense oligonucleotides. This modification offers significantly increased metabolic stability as well as improved pharmacokinetic properties. However, due to the chiral nature of the phosphorothioate, every PS group doubles the amount of possible stereoisomers. Thus PS oligonucleotides are generally obtained as an inseparable mixture of a multitude of diastereoisomeric compounds. Herein, we describe the introduction of non-chiral 3′ thiophosphate linkages into antisense oligonucleotides and report their in vitro as well as in vivo activity. The obtained results are carefully investigated for the individual parameters contributing to antisense activity of 3′ and 5′ thiophosphate modified oligonucleotides (target binding, RNase H recruitment, nuclease stability). We conclude that nuclease stability is the major challenge for this approach. These results highlight the importance of selecting meaningful in vitro experiments particularly when examining hitherto unexplored chemical modifications.

Insights into DNA catalysis from structural and functional studies of the 8-17 DNAzyme

The review examines functional knowledge gathered over two decades of research on the 8-17 DNAzyme, focusing on three aspects: the structural requirements for catalysis, the role of metal ions and the participation of general acid-base catalysis.

Small Molecule Rescue and Glycosidic Conformational Analysis of the Twister Ribozyme

The number of self-cleaving ribozymes has increased sharply in recent years, giving rise to elaborations of the four known ribozyme catalytic strategies, α, β, γ, and δ. One such extension is utilized by the twister ribozyme, which is hypothesized to conduct δ, or general acid catalysis, via N3 of the syn adenine +1 nucleobase indirectly via buffer catalysis at biological pH and directly at lower pH. Herein, we test the δ catalysis role of A1 via chemical rescue and the catalytic relevance of the syn orientation of the nucleobase by conformational analysis. Using inhibited twister ribozyme variants with A1(N3) deaza or A1 abasic modifications, we observe >100-fold chemical rescue effects in the presence of protonatable biological small molecules such as imidazole and histidine, similar to observed rescue values previously reported for C75U/C76Δ in the HDV ribozyme. Brønsted plots for the twister variants support a model in which small molecules rescue catalytic activity via a proton transfer mechanism, suggesting that A1 in the wild type is involved in proton transfer, most likely general acid catalysis. Additionally, through glycosidic conformational analysis in an appropriate background that accommodates the bromine atom, we observe that an 8BrA1-modified twister ribozyme is up to 10-fold faster than a nonmodified A1 ribozyme, supporting crystallographic data that show that A1 is syn when conducting proton transfer. Overall, this study provides functional evidence that the nucleotide immediately downstream of the cleavage site participates directly or indirectly in general acid-base catalysis in the twister ribozyme while occupying the syn conformation.

Novel ribozymes: discovery, catalytic mechanisms, and the quest to understand biological function

Small endonucleolytic ribozymes promote the self-cleavage of their own phosphodiester backbone at a specific linkage. The structures of and the reactions catalysed by members of individual families have been studied in great detail in the past decades. In recent years, bioinformatics studies have uncovered a considerable number of new examples of known catalytic RNA motifs. Importantly, entirely novel ribozyme classes were also discovered, for most of which both structural and biochemical information became rapidly available. However, for the majority of the new ribozymes, which are found in the genomes of a variety of species, a biological function remains elusive. Here, we concentrate on the different approaches to find catalytic RNA motifs in sequence databases. We summarize the emerging principles of RNA catalysis as observed for small endonucleolytic ribozymes. Finally, we address the biological functions of those ribozymes, where relevant information is available and common themes on their cellular activities are emerging. We conclude by speculating on the possibility that the identification and characterization of proteins that we hypothesize to be endogenously associated with catalytic RNA might help in answering the ever-present question of the biological function of the growing number of genomically encoded, small endonucleolytic ribozymes.

Classification of the nucleolytic ribozymes based upon catalytic mechanism

The nucleolytic ribozymes carry out site-specific RNA cleavage reactions by nucleophilic attack of the 2'-oxygen atom on the adjacent phosphorus with an acceleration of a million-fold or greater. A major part of this arises from concerted general acid-base catalysis. Recent identification of new ribozymes has expanded the group to a total of nine and this provides a new opportunity to identify sub-groupings according to the nature of the general base and acid. These include nucleobases, hydrated metal ions, and 2'-hydroxyl groups. Evolution has selected a number of different combinations of these elements that lead to efficient catalysis. These differences provide a new mechanistic basis for classifying these ribozymes.

Comparison of the Structures and Mechanisms of the Pistol and Hammerhead Ribozymes

Comparison of the secondary and three-dimensional structures of the hammerhead and pistol ribozymes reveals many close similarities, so in this work we have asked if they are mechanistically identical. We have determined a new crystal structure of the pistol ribozyme and have shown that G40 acts as general base in the cleavage reaction. The conformation in the active site ensures an in-line attack of the O2' nucleophile, and the conformation at the scissile phosphate and the position of the general base are closely similar to those in the hammerhead ribozyme. However, the two ribozymes differ in the nature of the general acid. 2'-Amino substitution experiments indicate that the general acid of the hammerhead ribozyme is the O2' of G8, while that of the pistol ribozyme is a hydrated metal ion. The two ribozymes are related but mechanistically distinct.

Evidence That Nucleophile Deprotonation Exceeds Bond Formation in the HDV Ribozyme Transition State

Steric constraints imposed by the active sites of protein and RNA enzymes pose major challenges to the investigation of structure-function relationships within these systems. As a strategy to circumvent such constraints in the HDV ribozyme, we have synthesized phosphoramidites from propanediol derivatives and incorporated them at the 5'-termini of RNA and DNA oligonucleotides to generate a series of novel substrates with nucleophiles perturbed electronically through geminal fluorination. In nonenzymatic, hydroxide-catalyzed intramolecular transphosphorylation of the DNA substrates, pH-rate profiles revealed that fluorine substitution reduces the maximal rate and the kinetic p K_a, consistent with the expected electron-withdrawing effect. In HDV ribozyme reactions, we observed that the RNA substrates undergo transphosphorylation relatively efficiently, suggesting that the conformational constraints imposed by a ribofuranose ring are not strictly required for ribozyme catalysis. In contrast to the nonenzymatic reactions, however, substrate fluorination modestly increases the ribozyme reaction rate, consistent with a mechanism in which (1) the 2'-hydroxyl nucleophile exists predominantly in its neutral, protonated form in the ground state and (2) the 2'-hydroxyl bears some negative charge in the rate-determining step, consistent with a transition state in which the extent of 2'-OH deprotonation exceeds the extent of P-O bond formation.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Nucleic acids: function and potential for abiogenesis

The emergence of functional cooperation between the three main classes of biomolecules - nucleic acids, peptides and lipids - defines life at the molecular level. However, how such mutually interdependent molecular systems emerged from prebiotic chemistry remains a mystery. A key hypothesis, formulated by Crick, Orgel and Woese over 40 year ago, posits that early life must have been simpler. Specifically, it proposed that an early primordial biology lacked proteins and DNA but instead relied on RNA as the key biopolymer responsible not just for genetic information storage and propagation, but also for catalysis, i.e. metabolism. Indeed, there is compelling evidence for such an 'RNA world', notably in the structure of the ribosome as a likely molecular fossil from that time. Nevertheless, one might justifiably ask whether RNA alone would be up to the task. From a purely chemical perspective, RNA is a molecule of rather uniform composition with all four bases comprising organic heterocycles of similar size and comparable polarity and pK a values. Thus, RNA molecules cover a much narrower range of steric, electronic and physicochemical properties than, e.g. the 20 amino acid side-chains of proteins. Herein we will examine the functional potential of RNA (and other nucleic acids) with respect to self-replication, catalysis and assembly into simple protocellular entities.

Chemical Bypass of General Base Catalysis in Hedgehog Protein Cholesterolysis Using a Hyper-Nucleophilic Substrate

Proteins in the hedgehog family undergo self-catalyzed endoproteolysis involving nucleophilic attack by a molecule of cholesterol. Recently, a conserved aspartate residue (D303, or D46) of hedgehog was identified as the general base that activates cholesterol during this unusual autoprocessing event; mutation of the catalyzing functional group (D303A) reduces activity by >104-fold. Here we report near total rescue of this ostensibly dead general base mutant by a synthetic substrate, 3β-hydroperoxycholestane (3HPC) in which the sterol -OH group is replaced by the hyper nucleophilic -OOH group. Other hedgehog point mutants at D303, also unreactive with cholesterol, accepted 3HPC as a substrate with the rank order: WT > D303A ≈ D303N ≫ D303R, D303E. We attribute the revived activity with 3-HPC to the α-effect, where tandem electronegative atoms exhibit exceptionally high nucleophilicity despite relatively low basicity.

Complexity in pH-Dependent Ribozyme Kinetics: Dark pKa Shifts and Wavy Rate-pH Profiles

Charged bases occur in RNA enzymes, or ribozymes, where they play key roles in catalysis. Cationic bases donate protons and perform electrostatic catalysis, while anionic bases accept protons. We previously published simulations of rate-pH profiles for ribozymes in terms of species plots for the general acid and general base that have been useful for understanding how ribozymes respond to pH. In that study, we did not consider interaction between the general acid and general base or interaction with other species on the RNA. Since that report, diverse small ribozyme classes have been discovered, many of which have charged nucleobases or metal ions in the active site that can either directly interact and participate in catalysis or indirectly interact as "influencers". Herein, we simulate experimental rate-pH profiles in terms of species plots in which reverse protonated charged nucleobases interact. These analyses uncover two surprising features of pH-dependent enzyme kinetics. (1) Cooperativity between the general acid and general base enhances population of the functional forms of a ribozyme and manifests itself as hidden or "dark" pK_a shifts, real pK_a shifts that accelerate the reaction but are not readily observed by standard experimental approaches, and (2) influencers favorably shift the pK_as of proton-transferring nucleobases and manifest themselves as "wavy" rate-pH profiles. We identify parallels with the protein enzyme literature, including reverse protonation and wavelike behavior, while pointing out that RNA is more prone to reverse protonation. The complexities uncovered, which arise from simple pairwise interactions, should aid deconvolution of complex rate-pH profiles for RNA and protein enzymes and suggest veiled catalytic devices for promoting catalysis that can be tested by experiment and calculation.

Allosteric Modulation of the Faecalibacterium prausnitzii Hepatitis Delta Virus-like Ribozyme by Glucosamine 6-Phosphate: The Substrate of the Adjacent Gene Product

Self-cleaving ribozymes were discovered 30 years ago and have been found throughout nature, from bacteria to animals, but little is known about their biological functions and regulation, particularly how cofactors and metabolites alter their activity. A hepatitis delta virus-like self-cleaving ribozyme maps upstream of a phosphoglucosamine mutase (glmM) open reading frame in the genome of the human gut bacterium Faecalibacterium prausnitzii. The presence of a ribozyme in the untranslated region of glmM suggests a regulation mechanism of gene expression. In the bacterial hexosamine biosynthesis pathway, the enzyme glmM catalyzes the isomerization of glucosamine 6-phosphate into glucosamine 1-phosphate. In this study, we investigated the effect of these metabolites on the co-transcriptional self-cleavage rate of the ribozyme. Our results suggest that glucosamine 6-phosphate, but not glucosamine 1-phosphate, is an allosteric ligand that increases the self-cleavage rate of drz-Fpra-1, providing the first known example of allosteric modulation of a self-cleaving ribozyme by the substrate of the adjacent gene product. Given that the ribozyme is activated by the glmM substrate, but not the product, this allosteric modulation may represent a potential feed-forward mechanism of gene expression regulation in bacteria.

Synthesis of 5'-Thio-3'-O-ribonucleoside Phosphoramidites

The chemical synthesis of phosphoramidite derivatives of all four 5'-deoxy-5'-thioribonucleosides is described. These phosphoramidites contained trityl (A, G, C, and U), dimethoxytrityl (A and G), or tert-butyldisulfanyl (G) as the 5'-S-protecting group. The application of several of these phosphoramidites for solid-phase synthesis of oligoribonucleotides containing a 2'-O-photocaged 5'-S-phosphorothiolate linkage or 5'-thiol-labeled RNAs is also further investigated.

Model for the Functional Active State of the TS Ribozyme from Molecular Simulation

Recently, a crystal structure has been reported of a new catalytic RNA, the TS ribozyme, that has been identified through comparative genomics and is believed to be a metalloribozyme having novel mechanistic features. Although this data provides invaluable structural information, analysis suggests a conformational change is required to arrive at a catalytically relevant state. We report results of molecular simulations that predict a spontaneous local rearrangement of the active site, leading to solution structures consistent with available functional data and providing competing mechanistic hypotheses that can be experimentally tested. The two competing hypotheses differ in the proposed identity of the catalytic general acid: either a water molecule coordinating a Mg²⁺ ion bound at the Watson-Crick edge of residue C7, or the N3 position of residue C7 itself.

How RNA acts as a nuclease: some mechanistic comparisons in the nucleolytic ribozymes

Recent structural and mechanistic studies have shed considerable light on the catalytic mechanisms of nucleolytic ribozymes. The discovery of several new ribozymes in this class has now allowed comparisons to be made, and the beginnings of mechanistic groupings to emerge.

The structure of a nucleolytic ribozyme that employs a catalytic metal ion

The TS ribozyme (originally called "twister sister") is a catalytic RNA. We present a crystal structure of the ribozyme in a pre-reactive conformation. Two co-axial helical stacks are organized by a three-way junction and two tertiary contacts. Five divalent metal ions are directly coordinated to RNA ligands, making important contributions to the RNA architecture. The scissile phosphate lies in a quasihelical loop region that is organized by a network of hydrogen bonding. A divalent metal ion is directly bound to the nucleobase 5' to the scissile phosphate, with an inner-sphere water molecule positioned to interact with the O2' nucleophile. The rate of ribozyme cleavage correlated in a log-linear manner with divalent metal ion pKa, consistent with proton transfer in the transition state, and we propose that the bound metal ion is a likely general base for the cleavage reaction. Our data indicate that the TS ribozyme functions predominantly as a metalloenzyme.

Structural and Biochemical Properties of Novel Self-Cleaving Ribozymes

Fourteen well-defined ribozyme classes have been identified to date, among which nine are site-specific self-cleaving ribozymes. Very recently, small self-cleaving ribozymes have attracted renewed interest in their structure, biochemistry, and biological function since the discovery, during the last three years, of four novel ribozymes, termed twister, twister sister, pistol, and hatchet. In this review, we mainly address the structure, biochemistry, and catalytic mechanism of the novel ribozymes. They are characterized by distinct active site architectures and divergent, but similar, biochemical properties. The cleavage activities of the ribozymes are highly dependent upon divalent cations, pH, and base-specific mutations, which can cause changes in the nucleotide arrangement and/or electrostatic potential around the cleavage site. It is most likely that a guanine and adenine in close proximity of the cleavage site are involved in general acid-base catalysis. In addition, metal ions appear to play a structural rather than catalytic role although some of their crystal structures have shown a direct metal ion coordination to a non-bridging phosphate oxygen at the cleavage site. Collectively, the structural and biochemical data of the four newest ribozymes could contribute to advance our mechanistic understanding of how self-cleaving ribozymes accomplish their efficient site-specific RNA cleavages.

TrmD: A Methyl Transferase for tRNA Methylation With m1G37

TrmD is an S-adenosyl methionine (AdoMet)-dependent methyl transferase that synthesizes the methylated m¹G37 in tRNA. TrmD is specific to and essential for bacterial growth, and it is fundamentally distinct from its eukaryotic and archaeal counterpart Trm5. TrmD is unusual by using a topological protein knot to bind AdoMet. Despite its restricted mobility, the TrmD knot has complex dynamics necessary to transmit the signal of AdoMet binding to promote tRNA binding and methyl transfer. Mutations in the TrmD knot block this intramolecular signaling and decrease the synthesis of m¹G37-tRNA, prompting ribosomes to +1-frameshifts and premature termination of protein synthesis. TrmD is unique among AdoMet-dependent methyl transferases in that it requires Mg²⁺ in the catalytic mechanism. This Mg²⁺ dependence is important for regulating Mg²⁺ transport to Salmonella for survival of the pathogen in the host cell. The strict conservation of TrmD among bacterial species suggests that a better characterization of its enzymology and biology will have a broad impact on our understanding of bacterial pathogenesis.

Crystal structure of Pistol, a class of self-cleaving ribozyme

Small self-cleaving ribozymes have been discovered in all evolutionary domains of life. They can catalyze site-specific RNA cleavage, and as a result, they have relevance in gene regulation. Comparative genomic analysis has led to the discovery of a new class of small self-cleaving ribozymes named Pistol. We report the crystal structure of Pistol at 2.97-Å resolution. Our results suggest that the Pistol ribozyme self-cleavage mechanism likely uses a guanine base in the active site pocket to carry out the phosphoester transfer reaction. The guanine G40 is in close proximity to serve as the general base for activating the nucleophile by deprotonating the 2'-hydroxyl to initiate the reaction (phosphoester transfer). Furthermore, G40 can also establish hydrogen bonding interactions with the nonbridging oxygen of the scissile phosphate. The proximity of G32 to the O5' leaving group suggests that G32 may putatively serve as the general acid. The RNA structure of Pistol also contains A-minor interactions, which seem to be important to maintain its tertiary structure and compact fold. Our findings expand the repertoire of ribozyme structures and highlight the conserved evolutionary mechanism used by ribozymes for catalysis.

Unwinding the twister ribozyme: from structure to mechanism

The twister ribozyme motif has been identified by bioinformatic means very recently. Currently, four crystal structures with ordered active sites together with a series of chemical and biochemical data provide insights into how this RNA accomplishes its efficient self‐cleavage. Of particular interest for a mechanistic proposal are structural distinctions observed in the active sites that concern the conformation of the U‐A cleavage site dinucleotide (in‐line alignment of the attacking 2′‐O nucleophile to the to‐be‐cleaved P—O5′ bond versus suboptimal alignments) as well as the presence/absence of Mg²⁺ ions at the scissile phosphate. All structures support the notion that an active site guanine and the conserved adenine at the cleavage site are important contributors to cleavage chemistry, likely being involved in general acid base catalysis. Evidence for innersphere coordination of a Mg²⁺ ion to the pro‐S nonbridging oxygen of the scissile phosphate stems from two of the four crystal structures. Together with the finding of thio/rescue effects for phosphorothioate substrates, this suggests the participation of divalent ions in the overall catalytic strategy employed by twister ribozymes. In this context, it is notable that twister retains wild‐type activity when the phylogenetically conserved stem P1 is deleted, able to cleave a single nucleotide only. WIREs RNA 2017, 8:e1402. doi: 10.1002/wrna.1402

For further resources related to this article, please visit the WIREs website.

Molecular Dynamics Study of Twister Ribozyme: Role of Mg(2+) Ions and the Hydrogen-Bonding Network in the Active Site

The recently discovered twister ribozyme is thought to utilize general acid-base catalysis in its self-cleavage mechanism, but the roles of nucleobases and metal ions in the mechanism are unclear. Herein, molecular dynamics simulations of the env22 twister ribozyme are performed to elucidate the structural and equilibrium dynamical properties, as well as to examine the role of Mg(2+) ions and possible candidates for the general base and acid in the self-cleavage mechanism. The active site region and the ends of the pseudoknots were found to be less mobile than other regions of the ribozyme, most likely providing structural stability and possibly facilitating catalysis. A purported catalytic Mg(2+) ion and the closest neighboring Mg(2+) ion remained chelated and relatively immobile throughout the microsecond trajectories, although removal of these Mg(2+) ions did not lead to any significant changes in the structure or equilibrium motions of the ribozyme on the microsecond time scale. In addition, a third metal ion, a Na(+) ion remained close to A1(O5'), the leaving group atom, during the majority of the microsecond trajectories, suggesting that it might stabilize the negative charge on A1(O5') during self-cleavage. The locations of these cations and their interactions with key nucleotides in the active site suggest that they may be catalytically relevant. The P1 stem is partially melted at its top and bottom in the crystal structure and further unwinds in the trajectories. The simulations also revealed an interconnected network comprised of hydrogen-bonding and π-stacking interactions that create a relatively rigid network around the self-cleavage site. The nucleotides involved in this network are among the highly conserved nucleotides in twister ribozymes, suggesting that this interaction network may be important to structure and function.

Topological constraints of structural elements in regulation of catalytic activity in HDV-like self-cleaving ribozymes

Self-cleaving ribozymes fold into intricate structures, which orient active site groups into catalytically competent conformations. Most ribozyme families have distinct catalytic cores stabilized by tertiary interactions between domains peripheral to those cores. We show that large hepatitis delta virus (HDV)-like ribozymes are activated by peripheral domains that bring two helical segments, P1 and P2, into proximity – a “pinch” that results in rate acceleration by almost three orders of magnitude. Kinetic analysis of ribozymes with systematically altered length and stability of the peripheral domain revealed that about one third of its free energy of formation is used to lower an activation energy barrier, likely related to a rate-limiting conformational change leading to the pre-catalytic state. These findings provide a quantitative view of enzyme regulation by peripheral domains and may shed light on the energetics of allosteric regulation.

Catalytic DNA: Scope, Applications, and Biochemistry of Deoxyribozymes

The discovery of natural RNA enzymes (ribozymes) prompted the pursuit of artificial DNA enzymes (deoxyribozymes) by in vitro selection methods. A key motivation is the conceptual and practical advantages of DNA relative to proteins and RNA. Early studies focused on RNA-cleaving deoxyribozymes, and more recent experiments have expanded the breadth of catalytic DNA to many other reactions. Including modified nucleotides has the potential to widen the scope of DNA enzymes even further. Practical applications of deoxyribozymes include their use as sensors for metal ions and small molecules. Structural studies of deoxyribozymes are only now beginning; mechanistic experiments will surely follow. Following the first report 21 years ago, the field of deoxyribozymes has promise for both fundamental and applied advances in chemistry, biology, and other disciplines.

A Two-Metal-Ion-Mediated Conformational Switching Pathway for HDV Ribozyme Activation

RNA enzymes serve as a potentially powerful platform from which to design catalysts and engineer new biotechnology. A fundamental understanding of these systems provides insight to guide design. The hepatitis delta virus ribozyme (HDVr) is a small, self-cleaving RNA motif widely distributed in nature, that has served as a paradigm for understanding basic principles of RNA catalysis. Nevertheless, questions remain regarding the precise roles of divalent metal ions and key nucleotides in catalysis. In an effort to establish a reaction mechanism model consistent with available experimental data, we utilize molecular dynamics simulations to explore different conformations and metal ion binding modes along the HDVr reaction path. Building upon recent crystallographic data, our results provide a dynamic model of the HDVr reaction mechanism involving a conformational switch between multiple non-canonical G25:U20 base pair conformations in the active site. These local nucleobase dynamics play an important role in catalysis by modulating the metal binding environments of two Mg²⁺ ions that support catalysis at different steps of the reaction pathway. The first ion plays a structural role by inducing a base pair flip necessary to obtain the catalytic fold in which C75 moves towards to the scissile phosphate in the active site. Ejection of this ion then permits a second ion to bind elsewhere in the active site and facilitate nucleophile activation. The simulations collectively describe a mechanistic scenario that is consistent with currently available experimental data from crystallography, phosphorothioate substitutions, and chemical probing studies. Avenues for further experimental verification are suggested.

A fluorescent surrogate of thymidine in duplex DNA

DMAT is a new fluorescent thymidine mimic composed of 2′-deoxyuridine fused to dimethylaniline.

A Mini-Twister Variant and Impact of Residues/Cations on the Phosphodiester Cleavage of this Ribozyme Class

Nucleolytic ribozymes catalyze site-specific cleavage of their phosphodiester backbones. A minimal version of the twister ribozyme is reported that lacks the phylogenetically conserved stem P1 while retaining wild-type activity. Atomic mutagenesis revealed that nitrogen atoms N1 and N3 of the adenine-6 at the cleavage site are indispensable for cleavage. By NMR spectroscopy, a pKa value of 5.1 was determined for a (13) C2-labeled adenine at this position in the twister ribozyme, which is significantly shifted compared to the pKa of the same adenine in the substrate alone. This finding pinpoints at a potential role for adenine-6 in the catalytic mechanism besides the previously identified invariant guanine-48 and a Mg(2+) ion, both of which are directly coordinated to the non-bridging oxygen atoms of the scissile phosphate; for the latter, additional evidence stems from the observation that Mn(2+) or Cd(2+) accelerated cleavage of phosphorothioate substrates. The relevance of this metal ion binding site is further emphasized by a new 2.6 Å X-ray structure of a 2'-OCH3 -U5 modified twister ribozyme.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.

Assessment of metal-assisted nucleophile activation in the hepatitis delta virus ribozyme from molecular simulation and 3D-RISM

The hepatitis delta virus ribozyme is an efficient catalyst of RNA 2'-O-transphosphorylation and has emerged as a key experimental system for identifying and characterizing fundamental features of RNA catalysis. Recent structural and biochemical data have led to a proposed mechanistic model whereby an active site Mg(2+) ion facilitates deprotonation of the O2' nucleophile, and a protonated cytosine residue (C75) acts as an acid to donate a proton to the O5' leaving group as noted in a previous study. This model assumes that the active site Mg(2+) ion forms an inner-sphere coordination with the O2' nucleophile and a nonbridging oxygen of the scissile phosphate. These contacts, however, are not fully resolved in the crystal structure, and biochemical data are not able to unambiguously exclude other mechanistic models. In order to explore the feasibility of this model, we exhaustively mapped the free energy surfaces with different active site ion occupancies via quantum mechanical/molecular mechanical (QM/MM) simulations. We further incorporate a three-dimensional reference interaction site model for the solvated ion atmosphere that allows these calculations to consider not only the rate associated with the chemical steps, but also the probability of observing the system in the presumed active state with the Mg(2+) ion bound. The QM/MM results predict that a pathway involving metal-assisted nucleophile activation is feasible based on the rate-controlling transition state barrier departing from the presumed metal-bound active state. However, QM/MM results for a similar pathway in the absence of Mg(2+) are not consistent with experimental data, suggesting that a structural model in which the crystallographically determined Mg(2+) is simply replaced with Na(+) is likely incorrect. It should be emphasized, however, that these results hinge upon the assumption of the validity of the presumed Mg(2+)-bound starting state, which has not yet been definitively verified experimentally, nor explored in depth computationally. Thus, further experimental and theoretical study is needed such that a consensus view of the catalytic mechanism emerges.