Two Active Site Divalent Ions in the Crystal Structure of the Hammerhead Ribozyme Bound to a Transition State Analogue

Identification of Hammerhead-variant ribozyme sequences in SARS-CoV-2

The SARS-CoV-2 RNA virus and variants, responsible for the COVID-19 pandemic has become endemic, raised a need for further understanding of the viral genome and biology. Despite vast research on SARS-CoV-2, no ribozymes have been found in the virus genome. Here we report the identification of 39 Hammerhead-variant ribozyme sequences (CoV-HHRz) in SARS-CoV-2. These sequences are highly conserved within SARS-CoV-2 variants but show large diversity among other coronaviruses. In vitro CoV-HHRz sequences possess the characteristics of typical ribozymes; cleavage is pH and ion dependent, although their activity is relatively low and Mn2+ is required for cleavage. The cleavage sites of four CoV-HHRz coincide with the breakpoint of expressed subgenomic RNA (sgRNAs) in SARS-CoV-2 transcriptome data suggesting in vivo activity. The CoV-HHRz are involved in processing sgRNAs for ORF7b, ORF 10 and ORF1ab nsp13 which are essential for viral packaging and life cycle.

Interaction of Vanadium Complexes with Proteins: Revisiting the Reported Structures in the Protein Data Bank (PDB) since 2015

The structural determination and characterization of molecules, namely proteins and enzymes, is crucial to gaining a better understanding of their role in different chemical and biological processes. The continuous technical developments in the experimental and computational resources of X-ray diffraction (XRD) and, more recently, cryogenic Electron Microscopy (cryo-EM) led to an enormous growth in the number of structures deposited in the Protein Data Bank (PDB). Bioinorganic chemistry arose as a relevant discipline in biology and therapeutics, with a massive number of studies reporting the effects of metal complexes on biological systems, with vanadium complexes being one of the relevant systems addressed. In this review, we focus on the interactions of vanadium compounds (VCs) with proteins. Several types of binding are established between VCs and proteins/enzymes. Considering that the V-species that bind may differ from those initially added, the mentioned structural techniques are pivotal to clarifying the nature and variety of interactions of VCs with proteins and to proposing the mechanisms involved either in enzymatic inhibition or catalysis. As such, we provide an account of the available structural information of VCs bound to proteins obtained by both XRD and/or cryo-EM, mainly exploring the more recent structures, particularly those containing organic-based vanadium complexes.

Biological Catalysis and Information Storage Have Relied on N-Glycosyl Derivatives of β-D-Ribofuranose since the Origins of Life

Most naturally occurring nucleotides and nucleosides are N-glycosyl derivatives of β-d-ribose. These N-ribosides are involved in most metabolic processes that occur in cells. They are essential components of nucleic acids, forming the basis for genetic information storage and flow. Moreover, these compounds are involved in numerous catalytic processes, including chemical energy production and storage, in which they serve as cofactors or coribozymes. From a chemical point of view, the overall structure of nucleotides and nucleosides is very similar and simple. However, their unique chemical and structural features render these compounds versatile building blocks that are crucial for life processes in all known organisms. Notably, the universal function of these compounds in encoding genetic information and cellular catalysis strongly suggests their essential role in the origins of life. In this review, we summarize major issues related to the role of N-ribosides in biological systems, especially in the context of the origin of life and its further evolution, through the RNA-based World(s), toward the life we observe today. We also discuss possible reasons why life has arisen from derivatives of β-d-ribofuranose instead of compounds based on other sugar moieties.

Viroids: Non-Coding Circular RNAs Able to Autonomously Replicate and Infect Higher Plants

Viroids are a unique type of infectious agent, exclusively composed of a relatively small (246-430 nt), highly base-paired, circular, non-coding RNA. Despite the small size and non-coding nature, the more-than-thirty currently known viroid species infectious of higher plants are able to autonomously replicate and move systemically through the host, thereby inducing disease in some plants. After recalling viroid discovery back in the late 60s and early 70s of last century and discussing current hypotheses about their evolutionary origin, this article reviews our current knowledge about these peculiar infectious agents. We describe the highly base-paired viroid molecules that fold in rod-like or branched structures and viroid taxonomic classification in two families, Pospiviroidae and Avsunviroidae, likely gathering nuclear and chloroplastic viroids, respectively. We review current knowledge about viroid replication through RNA-to-RNA rolling-circle mechanisms in which host factors, notably RNA transporters, RNA polymerases, RNases, and RNA ligases, are involved. Systemic movement through the infected plant, plant-to-plant transmission and host range are also discussed. Finally, we focus on the mechanisms of viroid pathogenesis, in which RNA silencing has acquired remarkable importance, and also for the initiation of potential biotechnological applications of viroid molecules.

In vitro evolution of ribonucleases from expanded genetic alphabets

The ability of nucleic acids to catalyze reactions (as well as store and transmit information) is important for both basic and applied science, the first in the context of molecular evolution and the origin of life and the second for biomedical applications. However, the catalytic power of standard nucleic acids (NAs) assembled from just four nucleotide building blocks is limited when compared with that of proteins. Here, we assess the evolutionary potential of libraries of nucleic acids with six nucleotide building blocks as reservoirs for catalysis. We compare the outcomes of in vitro selection experiments toward RNA-cleavage activity of two nucleic acid libraries: one built from the standard four independently replicable nucleotides and the other from six, with the two added nucleotides coming from an artificially expanded genetic information system (AEGIS). Results from comparative experiments suggest that DNA libraries with increased chemical diversity, higher information density, and larger searchable sequence spaces are one order of magnitude richer reservoirs of molecules that catalyze the cleavage of a phosphodiester bond in RNA than DNA libraries built from a standard four-nucleotide alphabet. Evolved AEGISzymes with nitro-carrying nucleobase Z appear to exploit a general acid–base catalytic mechanism to cleave that bond, analogous to the mechanism of the ribonuclease A family of protein enzymes and heavily modified DNAzymes. The AEGISzyme described here represents a new type of catalysts evolved from libraries built from expanded genetic alphabets.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Solution Structure of NPSL2, A Regulatory Element in the oncomiR-1 RNA

The miR-17 ∼ 92a polycistron, also known as oncomiR-1, is commonly overexpressed in multiple cancers and has several oncogenic properties. OncomiR-1 encodes six constituent microRNAs (miRs), each enzymatically processed with different efficiencies. However, the structural mechanism that regulates this differential processing remains unclear. Chemical probing of oncomiR-1 revealed that the Drosha cleavage sites of pri-miR-92a are sequestered in a four-way junction. NPSL2, an independent stem loop element, is positioned just upstream of pri-miR-92a and sequesters a crucial part of the sequence that constitutes the basal helix of pri-miR-92a. Disruption of the NPSL2 hairpin structure could promote the formation of a pri-miR-92a structure that is primed for processing by Drosha. Thus, NPSL2 is predicted to function as a structural switch, regulating pri-miR-92a processing. Here, we determined the solution structure of NPSL2 using solution NMR spectroscopy. This is the first high-resolution structure of an oncomiR-1 element. NPSL2 adopts a hairpin structure with a large, but highly structured, apical and internal loops. The 10-bp apical loop contains a pH-sensitive A+·C mismatch. Additionally, several adenosines within the apical and internal loops have elevated pKa values. The protonation of these adenosines can stabilize the NPSL2 structure through electrostatic interactions. Our study provides fundamental insights into the secondary and tertiary structure of an important RNA hairpin proposed to regulate miR biogenesis.

Structural variants and modifications of hammerhead ribozymes targeting influenza A virus conserved structural motifs

The naturally occurring structure and biological functions of RNA are correlated, which includes hammerhead ribozymes. We proposed new variants of hammerhead ribozymes targeting conserved structural motifs of segment 5 of influenza A virus (IAV) (+)RNA. The variants carry structural and chemical modifications aiming to improve the RNA cleavage activity of ribozymes. We introduced an additional hairpin motif and attempted to select ribozyme-target pairs with sequence features that enable the potential formation of the trans-Hoogsteen interactions that are present in full-length, highly active hammerhead ribozymes. We placed structurally defined guanosine analogs into the ribozyme catalytic core. Herein, the significantly improved synthesis of 2′-deoxy-2′-fluoroarabinoguanosine derivatives is described. The most potent hammerhead ribozymes were applied to chimeric short hairpin RNA (shRNA)-ribozyme plasmid constructs to improve the antiviral activity of the two components. The modified hammerhead ribozymes showed moderate cleavage activity. Treatment of IAV-infected Madin-Darby canine kidney (MDCK) cells with the plasmid constructs resulted in significant inhibition of virus replication. Real-time PCR analysis revealed a significant (80%–88%) reduction in viral RNA when plasmids carriers were used. A focus formation assay (FFA) for chimeric plasmids showed inhibition of virus replication by 1.6–1.7 log₁₀ units, whereas the use of plasmids carrying ribozymes or shRNAs alone resulted in lower inhibition.

RNA Electrostatics: How Ribozymes Engineer Active Sites to Enable Catalysis

Electrostatic interactions are fundamental to RNA structure and function, and intimately influenced by solvation and the ion atmosphere. RNA enzymes, or ribozymes, are catalytic RNAs that are able to enhance reaction rates over a million-fold, despite having only a limited repertoire of building blocks and available set of chemical functional groups. Ribozyme active sites usually occur at junctions where negatively charged helices come together, and in many cases leverage this strained electrostatic environment to recruit metal ions in solution that can assist in catalysis. Similar strategies have been implicated in related artificially engineered DNA enzymes. Herein, we apply Poisson-Boltzmann, 3D-RISM, and molecular simulations to study a set of metal-dependent small self-cleaving ribozymes (hammerhead, pistol, and Varkud satellite) as well as an artificially engineered DNAzyme (8-17) to examine electrostatic features and their relation to the recruitment of monovalent and divalent metal ions important for activity. We examine several fundamental roles for these ions that include: (1) structural integrity of the catalytically active state, (2) pKa tuning of residues involved in acid-base catalysis, and (3) direct electrostatic stabilization of the transition state via Lewis acid catalysis. Taken together, these examples demonstrate how RNA electrostatics orchestrates the site-specific and territorial binding of metal ions to play important roles in catalysis.

The role of Na+ in catalysis by the 8-17 DNAzyme

The 8-17 DNAzyme is the most studied deoxyribozyme in terms of its molecular mechanism; hence it has become a model system to understand the basis behind DNA catalysis. New functional studies and the recent attainment of high-resolution X-ray structures, in addition to theoretical calculations have offered a great opportunity to gain a broader comprehension of its mechanism; however many aspects are unclear yet, especially regarding the precise role of metal ions in catalysis. Recently, molecular dynamics simulations have suggested for the first time a specific and dynamical participation of Na⁺ in the mechanism through the reaction pathway, besides the roles proposed for divalent metal cofactors. Herein, we present experimental evidence of a cooperative role of the monovalent cation Na⁺ in catalysis that is in line with these theoretical suggestions. Our findings show a clear influence of the concentration of Na⁺ on the activity of the 8-17 DNAzyme when Pb²⁺ is used as the cofactor. Interestingly, this effect is not noticed with Mg²⁺, indicating a particular contribution of the monovalent ion to catalysis that would operate preferentially with Pb²⁺. We have also found that Na⁺ affects the pK_a of the general base and the general acid, indicating its influence on general acid-base catalysis, already identified as part of the mechanism of the 8-17 DNAzyme. Finally, our results emphasize the need to consider Na⁺ carefully in the design and analysis of functional studies of catalytic DNAs and its possible specific role in their mechanisms.

Structural Basis for Fluorescence Activation by Pepper RNA

Pepper is a fluorogenic RNA aptamer tag that binds to a variety of benzylidene-cyanophenyl (HBC) derivatives with tight affinity and activates their fluorescence. To investigate how Pepper RNA folds to create a binding site for HBC, we used antibody-assisted crystallography to determine the structures of Pepper bound to HBC530 and HBC599 to 2.3 and 2.7 Å resolutions, respectively. The structural data show that Pepper folds into an elongated structure and organizes nucleotides within an internal bulge to create the ligand binding site, assisted by an out-of-plane platform created by tertiary interactions with an adjacent bulge. As predicted from a lack of K⁺ dependence, Pepper does not use a G-quadruplex to form a binding pocket for HBC. Instead, Pepper uses a unique base-quadruple·base-triple stack to sandwich the ligand with a U·G wobble pair. Site-bound Mg²⁺ ions support ligand binding structurally and energetically. This research provides insight into the structural features that allow the Pepper aptamer to bind HBC and show how Pepper's function may expand to allow the in vivo detection of other small molecules and metals.

Structural modifications as tools in mechanistic studies of the cleavage of RNA phosphodiester linkages

The cleavage of RNA phosphodiester bonds by RNase A and hammerhead ribozyme at neutral pH fundamentally differs from the spontaneous reactions of these bonds under the same conditions. While the predominant spontaneous reaction is isomerization of the 3',5'-phosphodiester linkages to their 2',5'-counterparts, this reaction has never been reported to compete with the enzymatic cleavage reaction, not even as a minor side reaction. Comparative kinetic measurements with structurally modified di-nucleoside monophosphates and oligomeric phosphodiesters have played an important role in clarification of mechanistic details of the buffer-independent and buffer-catalyzed reactions. More recently, heavy atom isotope effects and theoretical calculations have refined the picture. The primary aim of all these studies has been to form a solid basis for mechanistic analyses of the action of more complicated catalytic machineries. In other words, to contribute to conception of a plausible unified picture of RNA cleavage by biocatalysts, such as RNAse A, hammerhead ribozyme and DNAzymes. In addition, structurally modified trinucleoside monophosphates as transition state models for Group I and II introns have clarified some features of the action of large ribozymes.

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.

Theoretical Studies of the Self Cleavage Pistol Ribozyme Mechanism

Ribozymes are huge complex biological catalysts composed of a combination of RNA and proteins. Nevertheless, there is a reduced number of small ribozymes, the self-cleavage ribozymes, that are formed just by RNA and, apparently, they existed in cells of primitive biological systems. Unveiling the details of these “fossils” enzymes can contribute not only to the understanding of the origins of life but also to the development of new simplified artificial enzymes. A computational study of the reactivity of the pistol ribozyme carried out by means of classical MD simulations and QM/MM hybrid calculations is herein presented to clarify its catalytic mechanism. Analysis of the geometries along independent MD simulations with different protonation states of the active site basic species reveals that only the canonical system, with no additional protonation changes, renders reactive conformations. A change in the coordination sphere of the Mg2+ ion has been observed during the simulations, which allows proposing a mechanism to explain the unique mode of action of the pistol ribozyme by comparison with other ribozymes. The present results are at the center of the debate originated from recent experimental and theoretical studies on pistol ribozyme.

Hydrated metal ion as a general acid in the catalytic mechanism of the 8-17 DNAzyme

The RNA-cleaving 8-17 DNAzyme, which is a metalloenzyme that depends on divalent metal ions for its function, is the most studied catalytic DNA in terms of its mechanism. By the end of 2017, a report of the crystal structure of the enzyme-substrate complex in the presence of Pb²⁺ probed some of the previous findings and opened new questions, especially around the participation of the metal ion in the catalytic mechanism and the promiscuity exhibited by the enzyme in terms of the metal cofactor required for catalysis. In this article we explore the role of the divalent metal ion in the mechanism of the 8-17 DNAzyme as a general acid, by measuring the influence of pH over the activity of a slower variant of the enzyme in the presence of Pb²⁺. We replaced G14, which has been identified as a general base in the mechanism of the enzyme, by the unnatural analog 2-aminopurine, with a lower pK_a value of the N1 group. With this approach, we obtained a bell-shaped pH-rate profile with experimental pK_a values of 5.4 and 7.0. Comparing these results with previous pH-rate profiles in the presence of Mg²⁺, our findings suggest the stabilization of the 5'-O leaving group by the hydrated metal ion acting as a general acid, in addition to the activation of the 2'-OH nucleophile by the general base G14.

The Varkud Satellite Ribozyme: A Thirty-Year Journey through Biochemistry, Crystallography, and Computation

The discovery of catalytic RNAs or ribozymes introduced a new class of enzymes to biology. In addition to their increasingly important roles in modern life, ribozymes are key players in the RNA World hypothesis, which posits that life started or flourished with RNA supporting both genetic and enzymatic functions. Therefore, investigations into the mechanisms of ribozyme function provide an exciting opportunity to examine the foundational principles of biological catalysis. Ribozymes are also attractive model systems to investigate the relationship between structure and function in RNA. Endonucleolytic ribozymes represent the largest class of catalytic RNA, of which the Varkud satellite (VS) ribozyme is structurally the most complex. The last ribozyme to be discovered by accident, the VS ribozyme had eluded structural determination for over two decades. When we solved the first crystal structures of the VS ribozyme, an extensive body of biochemical and biophysical data had accumulated over the years with which we could evaluate the functional relevance of the structure. Conversely, the structures provided a new perspective from which to reexamine the functional data and test new hypotheses. The VS ribozyme is organized in a modular fashion where independently folding domains assemble into the active conformation of the ribozyme via three-way junctions. Structures of the VS ribozyme in complex with its substrate at different stages of activation enabled us to map the structural reorganization of the substrate that must precede catalysis. In addition to defining the global architecture of the RNA, the essential interactions between the substrate and catalytic domains, and the rearrangements in the substrate prior to catalysis, these structures provided detailed snapshots of the ribozyme active site, revealing potential catalytic interactions. High resolution structures of the active site bolstered the view that the catalytic mechanism involved nucleobase-mediated general acid-base catalysis and uncovered additional catalytic interactions between the cleavage site and catalytic residues. Informed by the crystal structures of the VS ribozyme, an integrated experimental and computational approach identified the key players and essential interactions that define the active site of the ribozyme. This confluence of biochemical, structural, and computational studies revealed the catalytic mechanism of the ribozyme at unprecedented detail. Additionally, comparative analyses of the active site structures of the VS ribozyme and other nucleic acid-based endoribonucleases revealed common architectural motifs and strikingly similar catalytic strategies. In this Account, we document the progress of VS ribozyme research starting from its discovery and extending to the elucidation of its detailed catalytic mechanism 30 years later.

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.

RNA-Puzzles Round IV: 3D structure predictions of four ribozymes and two aptamers

RNA-Puzzles is a collective endeavor dedicated to the advancement and improvement of RNA 3D structure prediction. With agreement from crystallographers, the RNA structures are predicted by various groups before the publication of the crystal structures. We now report the prediction of 3D structures for six RNA sequences: four nucleolytic ribozymes and two riboswitches. Systematic protocols for comparing models and crystal structures are described and analyzed. In these six puzzles, we discuss (i) the comparison between the automated web servers and human experts; (ii) the prediction of coaxial stacking; (iii) the prediction of structural details and ligand binding; (iv) the development of novel prediction methods; and (v) the potential improvements to be made. We show that correct prediction of coaxial stacking and tertiary contacts is essential for the prediction of RNA architecture, while ligand binding modes can only be predicted with low resolution and simultaneous prediction of RNA structure with accurate ligand binding still remains out of reach. All the predicted models are available for the future development of force field parameters and the improvement of comparison and assessment tools.

Single Mutation in Hammerhead Ribozyme Favors Cleavage Activity with Manganese over Magnesium

Hammerhead ribozymes are one of the most studied classes of ribozymes so far, from both the structural and biochemical point of views. The activity of most hammerhead ribozymes is cation-dependent. Mg2+ is one of the most abundant divalent cations in the cell and therefore plays a major role in cleavage activity for most hammerhead ribozymes. Besides Mg2+, cleavage can also occur in the presence of other cations such as Mn2+. The catalytic core of hammerhead ribozymes is highly conserved, which could contribute to a preference of hammerhead ribozymes toward certain cations. Here, we show a naturally occurring variation in the catalytic core of hammerhead ribozymes, A6C, that can favor one metallic ion, Mn2+, over several other cations.

The L-platform/L-scaffold framework: a blueprint for RNA-cleaving nucleic acid enzyme design

We develop an L-platform/L-scaffold framework we hypothesize may serve as a blueprint to facilitate site-specific RNA-cleaving nucleic acid enzyme design. Building on the L-platform motif originally described by Suslov and coworkers, we identify new critical scaffolding elements required to anchor a conserved general base guanine ("L-anchor") and bind functionally important metal ions at the active site ("L-pocket"). Molecular simulations, together with a broad range of experimental structural and functional data, connect the L-platform/L-scaffold elements to necessary and sufficient conditions for catalytic activity. We demonstrate that the L-platform/L-scaffold framework is common to five of the nine currently known naturally occurring ribozyme classes (Twr, HPr, VSr, HHr, Psr), and intriguingly from a design perspective, the framework also appears in an artificially engineered DNAzyme (8-17dz). The flexibility of the L-platform/L-scaffold framework is illustrated on these systems, highlighting modularity and trends in the variety of known general acid moieties that are supported. These trends give rise to two distinct catalytic paradigms, building on the classifications proposed by Wilson and coworkers and named for the implicated general base and acid. The "G + A" paradigm (Twr, HPr, VSr) exclusively utilizes nucleobase residues for chemistry, and the "G + M + " paradigm (HHr, 8-17dz, Psr) involves structuring of the "L-pocket" metal ion binding site for recruitment of a divalent metal ion that plays an active role in the chemical steps of the reaction. Finally, the modularity of the L-platform/L-scaffold framework is illustrated in the VS ribozyme where the "L-pocket" assumes the functional role of the "L-anchor" element, highlighting a distinct mechanism, but one that is functionally linked with the hammerhead ribozyme.

Confluence of theory and experiment reveals the catalytic mechanism of the Varkud satellite ribozyme

The Varkud satellite ribozyme catalyses site-specific RNA cleavage and ligation, and serves as an important model system to understand RNA catalysis. Here, we combine stereospecific phosphorothioate substitution, precision nucleobase mutation and linear free-energy relationship measurements with molecular dynamics, molecular solvation theory and ab initio quantum mechanical/molecular mechanical free-energy simulations to gain insight into the catalysis. Through this confluence of theory and experiment, we unify the existing body of structural and functional data to unveil the catalytic mechanism in unprecedented detail, including the degree of proton transfer in the transition state. Further, we provide evidence for a critical Mg²⁺ in the active site that interacts with the scissile phosphate and anchors the general base guanine in position for nucleophile activation. This novel role for Mg²⁺ adds to the diversity of known catalytic RNA strategies and unifies functional features observed in the Varkud satellite, hairpin and hammerhead ribozyme classes.

Crucial Roles of Two Hydrated Mg2+ Ions in Reaction Catalysis of the Pistol Ribozyme

Pistol ribozymes constitute a new class of small self‐cleaving RNAs. Crystal structures have been solved, providing three‐dimensional snapshots along the reaction coordinate of pistol phosphodiester cleavage, corresponding to the pre‐catalytic state, a vanadate mimic of the transition state, and the product. The results led to the proposed underlying chemical mechanism. Importantly, a hydrated Mg²⁺ ion remains innersphere‐coordinated to N7 of G33 in all three states, and is consistent with its likely role as acid in general acid base catalysis (δ and β catalysis). Strikingly, the new structures shed light on a second hydrated Mg²⁺ ion that approaches the scissile phosphate from its binding site in the pre‐cleavage state to reach out for water‐mediated hydrogen bonding in the cyclophosphate product. The major role of the second Mg²⁺ ion appears to be the stabilization of product conformation. This study delivers a mechanistic understanding of ribozyme‐catalyzed backbone cleavage.

Molecular crowding and RNA catalysis

Molecular crowding promotes RNA folding and catalysis and could have played vital roles in the evolution of primordial ribozymes and protocells.

From general base to general acid catalysis in a sodium-specific DNAzyme by a guanine-to-adenine mutation

Recently, a few Na+-specific RNA-cleaving DNAzymes were reported, where nucleobases are likely to play critical roles in catalysis. The NaA43 and NaH1 DNAzymes share the same 16-nt Na+-binding motif, but differ in one or two nucleotides in a small catalytic loop. Nevertheless, they display an opposite pH-dependency, implicating distinct catalytic mechanisms. In this work, rational mutation studies locate a catalytic adenine residue, A22, in NaH1, while previous studies found a guanine (G23) to be important for the catalysis of NaA43. Mutation with pKa-perturbed analogs, such as 2-aminopurine (∼3.8), 2,6-diaminopurine (∼5.6) and hypoxanthine (∼8.7) affected the overall reaction rate. Therefore, we propose that the N1 position of G23 (pKa ∼6.6) in NaA43 functions as a general base, while that of A22 (pKa ∼6.3) in NaH1 as a general acid. Further experiments with base analogs and a phosphorothioate-modified substrate suggest that the exocyclic amine in A22 and both of the non-bridging oxygens at the scissile phosphate are important for catalysis for NaH1. This is an interesting example where single point mutations can change the mechanism of cleavage from general base to general acid, and it can also explain this Na+-dependent DNAzyme scaffold being sensitive to a broad range of metal ions and molecules.

Dynamical ensemble of the active state and transition state mimic for the RNA-cleaving 8-17 DNAzyme in solution

We perform molecular dynamics simulations, based on recent crystallographic data, on the 8-17 DNAzyme at four states along the reaction pathway to determine the dynamical ensemble for the active state and transition state mimic in solution. A striking finding is the diverse roles played by Na+ and Pb2+ ions in the electrostatically strained active site that impact all four fundamental catalytic strategies, and share commonality with some features recently inferred for naturally occurring hammerhead and pistol ribozymes. The active site Pb2+ ion helps to stabilize in-line nucleophilic attack, provides direct electrostatic transition state stabilization, and facilitates leaving group departure. A conserved guanine residue is positioned to act as the general base, and is assisted by a bridging Na+ ion that tunes the pKa and facilitates in-line fitness. The present work provides insight into how DNA molecules are able to solve the RNA-cleavage problem, and establishes functional relationships between the mechanism of these engineered DNA enzymes with their naturally evolved RNA counterparts. This adds valuable information to our growing body of knowledge on general mechanisms of phosphoryl transfer reactions catalyzed by RNA, proteins and DNA.

Molecular simulations of the pistol ribozyme: unifying the interpretation of experimental data and establishing functional links with the hammerhead ribozyme

The pistol ribozyme (Psr) is among the most recently discovered RNA enzymes and has been the subject of experiments aimed at elucidating the mechanism. Recent biochemical studies have revealed exciting clues about catalytic interactions in the active site not apparent from available crystallographic data. The present work unifies the interpretation of the existing body of structural and functional data on Psr by providing a dynamical model for the catalytically active state in solution from molecular simulation. Our results suggest that a catalytic Mg2+ ion makes inner-sphere contact with G33:N7 and outer-sphere coordination to the pro-RP of the scissile phosphate, promoting electrostatic stabilization of the dianionic transition state and neutralization of the developing charge of the leaving group through a metal-coordinated water molecule that is made more acidic by a hydrogen bond donated from the 2'OH of P32. This model is consistent with experimental activity-pH and mutagenesis data, including sensitivity to G33(7cG) and phosphorothioate substitution/metal ion rescue. The model suggests several experimentally testable predictions, including the response of cleavage activity to mutations at G42 and P32 positions in the ribozyme, and thio substitutions of the substrate in the presence of different divalent metal ions. Further, the model identifies striking similarities of Psr to the hammerhead ribozyme (HHr), including similar global fold, organization of secondary structure around an active site three-way junction, catalytic metal ion binding mode, and guanine general base. However, the specific binding mode and role of the Mg2+ ion, as well as a conserved 2'-OH in the active site, are interrelated but subtly different between the ribozymes.

RiboSketch: versatile visualization of multi-stranded RNA and DNA secondary structure

Summary

Creating clear, visually pleasing 2D depictions of RNA and DNA strands and their interactions is important to facilitate and communicate insights related to nucleic acid structure. Here we present RiboSketch, a secondary structure image production application that enables the visualization of multistranded structures via layout algorithms, comprehensive editing capabilities, and a multitude of simulation modes. These interactive features allow RiboSketch to create publication quality diagrams for structures with a wide range of composition, size and complexity. The program may be run in any web browser without the need for installation, or as a standalone Java application.

Availability and implementation

https://rnastructure.cancer.gov/ribosketch.

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.

Hatchet ribozyme structure and implications for cleavage mechanism

Self-cleaving ribozymes are RNAs that catalyze position-specific cleavage of their phosphodiester backbone. The cleavage site of the newly discovered hatchet ribozyme is located at the very 5′ end of its consensus secondary structure motif. Here we report on the 2.1-Å crystal structure of the hatchet ribozyme in the product state, which defines its intricate tertiary fold and identifies key residues lining the catalytic pocket. This in turn has allowed us to propose a model of the precatalytic state structure and a role in catalysis for a conserved guanine. This study therefore provides a structure-based platform toward an improved understanding of the catalytic mechanism of hatchet ribozymes.

Small self-cleaving ribozymes catalyze site-specific cleavage of their own phosphodiester backbone with implications for viral genome replication, pre-mRNA processing, and alternative splicing. We report on the 2.1-Å crystal structure of the hatchet ribozyme product, which adopts a compact pseudosymmetric dimeric scaffold, with each monomer stabilized by long-range interactions involving highly conserved nucleotides brought into close proximity of the scissile phosphate. Strikingly, the catalytic pocket contains a cavity capable of accommodating both the modeled scissile phosphate and its flanking 5′ nucleoside. The resulting modeled precatalytic conformation incorporates a splayed-apart alignment at the scissile phosphate, thereby providing structure-based insights into the in-line cleavage mechanism. We identify a guanine lining the catalytic pocket positioned to contribute to cleavage chemistry. The functional relevance of structure-based insights into hatchet ribozyme catalysis is strongly supported by cleavage assays monitoring the impact of selected nucleobase and atom-specific mutations on ribozyme activity.

RNA Back and Forth: Looking through Ribozyme and Viroid Motifs

Current cellular facts allow us to follow the link from chemical to biochemical metabolites, from the ancient to the modern world. In this context, the "RNA world" hypothesis proposes that early in the evolution of life, the ribozyme was responsible for the storage and transfer of genetic information and for the catalysis of biochemical reactions. Accordingly, the hammerhead ribozyme (HHR) and the hairpin ribozyme belong to a family of endonucleolytic RNAs performing self-cleavage that might occur during replication. Furthermore, regarding the widespread occurrence of HHRs in several genomes of modern organisms (from mammals to small parasites and elsewhere), these small ribozymes have been regarded as living fossils of a primitive RNA world. They fold into 3D structures that generally require long-range intramolecular interactions to adopt the catalytically active conformation under specific physicochemical conditions. By studying viroids as plausible remains of ancient RNA, we recently demonstrated that they replicate in non-specific hosts, emphasizing their adaptability to different environments, which enhanced their survival probability over the ages. All these results exemplify ubiquitous features of life. Those are the structural and functional versatility of small RNAs, ribozymes, and viroids, as well as their diversity and adaptability to various extreme conditions. All these traits must have originated in early life to generate novel RNA populations.

Bulky cations greatly increase the turnover of a native hammerhead ribozyme

Methods to facilitate the catalytic turnover of ribozymes are required for advancing oligonucleotide-based technologies. This study examined tetraalkylammonium ions for their ability to increase the efficiency of catalytic turnover of a native hammerhead ribozyme. Kinetic analysis showed that large tetraalkylammonium ions significantly increased the turnover rate of the ribozyme and was much more effective than poly(ethylene glycol) (PEG) and urea. The magnitude of the rate increase depended on the concentrations of Mg²⁺ and tetrapentylammonium ions, and the rate was enhanced by more than 180-fold at the optimal concentrations of these salts. The results provide physical insights into interactions of ribozymes with large cationic molecules through electrostatic forces and steric hindrance.

Large tetraalkylammonium ions increase the turnover rate of the ribozyme derived from an intronic ribozyme in the human genome. The rate can be enhanced by more than a hundred-fold at the optimal concentrations of Mg²⁺ and TPeA ions.

Structural Simplicity and Mechanistic Complexity in the Hammerhead Ribozyme

Natural or full-length hammerhead ribozymes are up to 1000-fold more active than their minimal counterparts that lack a complex tertiary interaction that pre-organizes and stabilizes the ribozyme active site, positioning RNA functional groups to facilitate acid-base catalysis. The recent discovery that a single tertiary contact (an AU Hoogsteen pair) between Stems I and II confers essentially all of the enhanced activity greatly simplifies our understanding of the structural requirements for hammerhead ribozyme activity. In contrast, the simplest mechanistic interpretations are challenged with the presentation of more complex alternatives. These alternatives are elucidated and critically analyzed in the context of several of the active hammerhead ribozyme structures now available.

NMR spectroscopic characterization of a model RNA duplex reflecting the core sequence of hammerhead ribozymes

Hammerhead ribozymes are a model system for studying molecular mechanism of RNA catalysis. Physicochemical data-driven mechanistic studies are an indispensable step towards understanding the catalysis of hammerhead ribozymes. Here we characterized a model RNA duplex with catalytically important sheared-type G12-A9 base pair and A9-G10.1 metal ion-binding motif in hammerhead ribozymes. By using high magnetic field NMR, all base proton signals, including catalytic residues, were unambiguously assigned. We further characterized structural features of this RNA molecule and found that it reflects the structural features of the A9-G10.1 motif of hammerhead ribozymes. Therefore, this RNA molecule is suitable for extracting an intrinsic physicochemical properties of catalytically important residues.

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.

Elucidation of Catalytic Strategies of Small Nucleolytic Ribozymes From Comparative Analysis of Active Sites

A number of small, self-cleaving ribozyme classes have been identified including the hammerhead, hairpin, hepatitis delta virus (HDV), Varkud satellite (VS), glmS, twister, hatchet, pistol, and twister sister ribozymes. Within the active sites of these ribozymes, myriad functional groups contribute to catalysis. There has been extensive structure-function analysis of individual ribozymes, but the extent to which catalytic devices are shared across different ribozyme classes is unclear. As such, emergent catalytic principles for ribozymes may await discovery. Identification of conserved catalytic devices can deepen our understanding of RNA catalysis specifically and of enzymic catalysis generally. To probe similarities and differences amongst ribozyme classes, active sites from more than 80 high-resolution crystal structures of self-cleaving ribozymes were compared computationally. We identify commonalities amongst ribozyme classes pertaining to four classic catalytic devices: deprotonation of the 2'OH nucleophile (γ), neutralization of the non-bridging oxygens of the scissile phosphate (β), neutralization of the O5' leaving group (δ), and in-line nucleophilic attack (α). In addition, we uncover conservation of two catalytic devices, each of which centers on the activation of the 2'OH nucleophile by a guanine: one to acidify the 2'OH by hydrogen bond donation to it (γ') and one to acidify the 2'OH by releasing it from non-productive interactions by competitive hydrogen bonding (γ''). Our findings reveal that the amidine functionalities of G, A, and C are especially important for these strategies, and help explain absence of U at ribozyme active sites. The identified γ' and γ'' catalytic strategies help unify the catalytic strategies shared amongst catalytic RNAs and may be important for large ribozymes, as well as protein enzymes that act on nucleic acids.

Inner-Sphere Coordination of Divalent Metal Ion with Nucleobase in Catalytic RNA

Identification of the function of metal ions and the RNA moieties, particularly nucleobases, that bind metal ions is important in RNA catalysis. Here we combine single-atom and abasic substitutions to probe functions of conserved nucleobases in ribonuclease P (RNase P). Structural and biophysical studies of bacterial RNase P propose direct coordination of metal ions by the nucleobases of conserved uridine and guanosine in helix P4 of the RNA subunit (P RNA). To biochemically probe the function of metal ion interactions, we substituted the universally conserved bulged uridine (U51) in the P4 helix of circularly permuted Bacillus subtilis P RNA with 4-thiouridine, 4-deoxyuridine, and abasic modifications and G378/379 with 2-aminopurine, N7-deazaguanosine, and 6-thioguanosine. The functional group modifications of U51 decrease RNase P-catalyzed phosphodiester bond cleavage 16- to 23-fold, as measured by the single-turnover cleavage rate constant. The activity of the 4-thiouridine RNase P is partially rescued by addition of Cd(II) or Mn(II) ions. This is the first time a metal-rescue experiment provides evidence for inner-sphere divalent metal ion coordination with a nucleobase. Modifications of G379 modestly decrease the cleavage activity of RNase P, suggesting outer-sphere coordination of O6 on G379 to a metal ion. These data provide biochemical evidence for catalytically important interactions of the P4 helix of P RNA with metal ions, demonstrating that the bulged uridine coordinates at least one catalytic metal ion through an inner-sphere interaction. The combination of single-atom and abasic nucleotide substitutions provides a powerful strategy to probe functions of conserved nucleobases in large RNAs.

Structure-based mechanistic insights into catalysis by small self-cleaving ribozymes

Small self-cleaving ribozymes are widely distributed in nature and are essential for rolling-circle-based replication of satellite and pathogenic RNAs. Earlier structure-function studies on the hammerhead, hairpin, glmS, hepatitis delta virus and Varkud satellite ribozymes have provided insights into their overall architecture, their catalytic active site organization, and the role of nearby nucleobases and hydrated divalent cations in facilitating general acid-base and electrostatic-mediated catalysis. This review focuses on recent structure-function research on active site alignments and catalytic mechanisms of the Rzb hammerhead ribozyme, as well as newly-identified pistol, twister and twister-sister ribozymes. In contrast to an agreed upon mechanistic understanding of self-cleavage by Rzb hammerhead and pistol ribozymes, there exists a divergence of views as to the cleavage site alignments and catalytic mechanisms adopted by twister and twister-sister ribozymes. One approach to resolving this conundrum would be to extend the structural studies from currently available pre-catalytic conformations to their transition state mimic vanadate counterparts for both ribozymes.

Structure-based insights into self-cleavage by a four-way junctional twister-sister ribozyme

Here we report on the crystal structure and cleavage assays of a four-way junctional twister-sister self-cleaving ribozyme. Notably, 11 conserved spatially separated loop nucleotides are brought into close proximity at the ribozyme core through long-range interactions mediated by hydrated Mg²⁺ cations. The C62–A63 step at the cleavage site adopts a splayed-apart orientation, with flexible C62 directed outwards, whereas A63 is directed inwards and anchored by stacking and hydrogen-bonding interactions. Structure-guided studies of key base, sugar, and phosphate mutations in the twister-sister ribozyme, suggest contributions to the cleavage chemistry from interactions between a guanine at the active site and the non-bridging oxygen of the scissile phosphate, a feature found previously also for the related twister ribozyme. Our four-way junctional pre-catalytic structure differs significantly in the alignment at the cleavage step (splayed-apart vs. base-stacked) and surrounding residues and hydrated Mg²⁺ ions relative to a reported three-way junctional pre-catalytic structure of the twister-sister ribozyme.

Twister-sister is a self-cleaving ribozyme. Here, the authors report the 2.0 Å crystal structure of the four-way junctional twister-sister ribozyme in the pre-catalytic state and discuss mechanistic implications based on their mutagenesis experiments and comparisons with other ribozyme structures.

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.

Small Molecule Recognition Triggers Secondary and Tertiary Interactions in DNA Folding and Hammerhead Ribozyme Catalysis

We have identified tris(2-aminoethyl)amine (tren)-derived scaffolds with two (t2M) or four (t4M) melamine rings that can target oligo T/U domains in DNA/RNA. Unstructured T-rich DNAs cooperatively fold with the tren derivatives to form hairpin-like structures. Both t2M and t4M act as functional switches in a family of hammerhead ribozymes deactivated by stem or loop replacement with a U-rich sequence. Catalysis of bond scission in these hammerhead ribozymes could be restored by putative t2M/t4M refolding of stem secondary structure or tertiary bridging interactions between loop and stem. The simplicity of the t2M/t4M binding site enables programming of allostery in RNAs, recoding oligo-U domains as potential sites for secondary structure or tertiary contact. In combination with a facile and general method for installation of the t2M motif on primary amines, the method described herein streamlines design of synthetic allosteric riboswitches and small molecule-nucleic acid complexes.

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.

Cooperative Interactions in the Hammerhead Ribozyme Drive pKa Shifting of G12 and Its Stacked Base C17

General acid-base catalysis is a key mechanistic strategy in protein and RNA enzymes. Ribozymes use hydrated metal ions, nucleobases, and organic cofactors to carry this out. In most small ribozymes, a guanosine is positioned to participate in proton transfer with the nucleophilic 2'-OH. The unshifted pK_a values for nucleobases and solvated metal ions are far from neutrality, however, and thus nonideal for general acid-base catalysis. Herein, evidence is provided for cooperative interaction in the hammerhead ribozyme among the guanine that interacts with the nucleophilic 2'-OH, G12, the -1 nucleobase C17, and Mg²⁺ ions. We introduce global fitting for analyzing ribozyme rate-pH data parametric in Mg²⁺ concentration and benchmark this method on data from the hepatitis delta virus ribozyme. We then apply global fitting to new rate-pH data for the hammerhead ribozyme using a minimal three-dimensional, four-channel cooperative model. The value for the pK_a of G12 that we obtain is channel-dependent and varies from 8.1 to 9.9, shifting closest toward neutrality in the presence of two cationic species: C17H⁺ and a Mg²⁺ ion. The value for the pK_a of the -1 nucleotide, C17, is increased a remarkable 3.5-5 pK_a units toward neutrality. Shifting of the pK_a of C17 appears to be driven by an electrostatic sandwich of C17 between carbonyl groups of the 5'-neighboring U and of G12 and involves cation-π interactions. Rate-pH profiles reveal that the major reactive channel under biological Mg²⁺ and pH involves a cationic C17 rather than a second metal ion. Substitution of a cationic base for a metal underscores the versatility of RNA.

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.