A comparison of vanadate to a 2'-5' linkage at the active site of a small ribozyme suggests a role for water in transition-state stabilization

Simulation-Guided Conformational Space Exploration to Assess Reactive Conformations of a Ribozyme

Self-splicing ribozymes are small ribonucleic acid (RNA) enzymes that catalyze their own cleavage through a transphosphoesterification reaction. While this process is involved in some specific steps of viral RNA replication and splicing, it is also of importance in the context of the (putative) first autocatalytic RNA-based systems that could have preceded the emergence of modern life. The uncatalyzed phosphoester bond formation is thermodynamically very unfavorable, and many experimental studies have focused on understanding the molecular features of catalysis in these ribozymes. However, chemical reaction paths are short-lived and not easily characterized by experimental approaches, so molecular simulation approaches appear as an ideal tool to unveil the molecular details of the reaction. Here, we focus on the model hairpin ribozyme. We show that identifying a relevant initial conformation for reactivity studies, which is frequently overlooked in mixed quantum-classical studies that predominantly concentrate on the chemical reaction itself, can be highly challenging. These challenges stem from limitations in both available experimental structures (which are chemically altered to prevent self-cleavage) and the accuracy of force fields, together with the necessity for comprehensive sampling. We show that molecular dynamics simulations, combined with extensive conformational phase space exploration with Hamiltonian replica-exchange simulations, enable us to characterize the relevant conformational basins of the minimal hairpin ribozyme in the ligated state prior to self-cleavage. We find that what is usually considered a canonical reactive conformation with active site geometries and hydrogen-bond patterns that are optimal for the addition-elimination reaction with general acid/general base catalysis is metastable and only marginally populated. The thermodynamically stable conformation appears to be consistent with the expectations of a mechanism that does not require the direct participation of ribozyme residues in the reaction. While these observations may suffer from forcefield inaccuracies, all investigated forcefields lead to the same conclusions upon proper sampling, contrasting with previous investigations on shorter timescales suggesting that at least one reparametrization of the Amber99 forcefield allowed to stabilize aligned active site conformations. Our study demonstrates that identifying the most pertinent reactant state conformation holds equal importance alongside the accurate determination of the thermodynamics and kinetics of the chemical steps of the reaction.

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Investigation of the pKa of the Nucleophilic O2' of the Hairpin Ribozyme

Small ribozymes cleave their RNA phosphodiester backbone by catalyzing a transphosphorylation reaction wherein a specific O2' functions as the nucleophile. While deprotonation of this alcohol through its acidification would increase its nucleophilicity, little is known about the pK_a of this O2' in small ribozymes, in part because high pK_a's are not readily accessible experimentally. Herein, we turn to molecular dynamics to calculate the pK_a of the nucleophilic O2' in the hairpin ribozyme and to study interactions within the active site that may impact its value. We estimate the pK_a of the nucleophilic O2' in the wild-type hairpin ribozyme to be 18.5 ± 0.8, which is higher than the reference compound, and identify a correlation between proper positioning of the O2' for nucleophilic attack and elevation of its pK_a. We find that monovalent ions may play a role in depression of the O2' pK_a, while the exocyclic amine appears to be important for organizing the ribozyme active site. Overall, this study suggests that the pK_a of the O2' is raised in the ground state and lowers during the course of the reaction owing to positioning and metal ion interactions.

RR3DD: an RNA global structure-based RNA three-dimensional structural classification database

The three-dimensional (3D) structure of RNA usually plays an important role in the recognition with RNA-binding protein. Along with the discovering of RNAs, several RNA databases are developed to study the functions of RNA based on sequence, secondary structure, local 3D structural motif and global structure. Based on RNA function and structure, different RNAs are classified and stored in SCOR and DARTS, respectively. The classification of RNA structures is useful in RNA structure prediction and function annotation. However, the SCOR and DARTS are not updated any more. In this study, we present an RNA classification database RR3DD based on RNA fold with the global 3D structural similarity. The RR3DD includes 13,601 RNA chains from PDB and mmCIF format structures which are classified into 780 RNA folds. The RNA chains from PDB and mmCIF format structures are aligned and clustered into 675 and 220 RNA folds, respectively. By analysing the RNA structure in RR3DD, we find that there are 11 clusters with more than 50 members. These clusters include rRNAs, riboswitches, tRNAs and so on. By mapping RR3DD into Rfam, we found that some RNAs without annotation by Rfam can be annotated through structural alignment. For example, we analysed tRNAs and found that tRNA were successfully grouped in RR3DD for which Rfam did not classify them into one family. Finally, we provide a web interface of RR3DD offering functions of browsing RR3DD, annotating RNA 3D structure and finding templates for RNA homology modelling.

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.

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.

Intrinsic Base-Pair Rearrangement in the Hairpin Ribozyme Directs RNA Conformational Sampling and Tertiary Interface Formation

Dynamic fluctuations in RNA structure enable conformational changes that are required for catalysis and recognition. In the hairpin ribozyme, the catalytically active structure is formed as an intricate tertiary interface between two RNA internal loops. Substantial alterations in the structure of each loop are observed upon interface formation, or docking. The very slow on-rate for this relatively tight interaction has led us to hypothesize a double conformational capture mechanism for RNA-RNA recognition. We used extensive molecular dynamics simulations to assess conformational sampling in the undocked form of the loop domain containing the scissile phosphate (loop A). We observed several major accessible conformations with distinctive patterns of hydrogen bonding and base stacking interactions in the active-site internal loop. Several important conformational features characteristic of the docked state were observed in well-populated substates, consistent with the kinetic sampling of docking-competent states by isolated loop A. Our observations suggest a hybrid or multistage binding mechanism, in which initial conformational selection of a docking-competent state is followed by induced-fit adjustment to an in-line, chemically reactive state only after formation of the initial complex with loop B.

Multiscale methods for computational RNA enzymology

RNA catalysis is of fundamental importance to biology and yet remains ill-understood due to its complex nature. The multidimensional "problem space" of RNA catalysis includes both local and global conformational rearrangements, changes in the ion atmosphere around nucleic acids and metal ion binding, dependence on potentially correlated protonation states of key residues, and bond breaking/forming in the chemical steps of the reaction. The goal of this chapter is to summarize and apply multiscale modeling methods in an effort to target the different parts of the RNA catalysis problem space while also addressing the limitations and pitfalls of these methods. Classical molecular dynamics simulations, reference interaction site model calculations, constant pH molecular dynamics (CpHMD) simulations, Hamiltonian replica exchange molecular dynamics, and quantum mechanical/molecular mechanical simulations will be discussed in the context of the study of RNA backbone cleavage transesterification. This reaction is catalyzed by both RNA and protein enzymes, and here we examine the different mechanistic strategies taken by the hepatitis delta virus ribozyme and RNase A.

Strand-specific (asymmetric) contribution of phosphodiester linkages on RNA polymerase II transcriptional efficiency and fidelity

Nonenzymatic RNA polymerization in early life is likely to introduce backbone heterogeneity with a mixture of 2'-5' and 3'-5' linkages. On the other hand, modern nucleic acids are dominantly composed of 3'-5' linkages. RNA polymerase II (pol II) is a key modern enzyme responsible for synthesizing 3'-5'-linked RNA with high fidelity. It is not clear how modern enzymes, such as pol II, selectively recognize 3'-5' linkages over 2'-5' linkages of nucleic acids. In this work, we systematically investigated how phosphodiester linkages of nucleic acids govern pol II transcriptional efficiency and fidelity. Through dissecting the impacts of 2'-5' linkage mutants in the pol II catalytic site, we revealed that the presence of 2'-5' linkage in RNA primer only modestly reduces pol II transcriptional efficiency without affecting pol II transcriptional fidelity. In sharp contrast, the presence of 2'-5' linkage in DNA template leads to dramatic decreases in both transcriptional efficiency and fidelity. These distinct effects reveal that pol II has an asymmetric (strand-specific) recognition of phosphodiester linkage. Our results provided important insights into pol II transcriptional fidelity, suggesting essential contributions of phosphodiester linkage to pol II transcription. Finally, our results also provided important understanding on the molecular basis of nucleic acid recognition and genetic information transfer during molecular evolution. We suggest that the asymmetric recognition of phosphodiester linkage by modern nucleic acid enzymes likely stems from the distinct evolutionary pressures of template and primer strand in genetic information transfer during molecular evolution.

Evidence for the role of active site residues in the hairpin ribozyme from molecular simulations along the reaction path

The hairpin ribozyme accelerates a phosphoryl transfer reaction without catalytic participation of divalent metal ions. Residues A38 and G8 have been implicated as playing roles in general acid and base catalysis, respectively. Here we explore the structure and dynamics of key active site residues using more than 1 μs of molecular dynamics simulations of the hairpin ribozyme at different stages along the catalytic pathway. Analysis of results indicates hydrogen bond interactions between the nucleophile and proR nonbridging oxygen are correlated with active inline attack conformations. Further, the simulation results suggest a possible alternative role for G8 to promote inline fitness and facilitate activation of the nucleophile by hydrogen bonding, although this does not necessarily exclude an additional role as a general base. Finally, we suggest that substitution of G8 with N7- or N3-deazaguanosine which have elevated pK_a values, both with and without thio modifications at the 5′ leaving group position, would provide valuable insight into the specific role of G8 in catalysis.

Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme

We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.

Intermolecular domain docking in the hairpin ribozyme: metal dependence, binding kinetics and catalysis

The hairpin ribozyme is a prototype small, self-cleaving RNA motif. It exists naturally as a four-way RNA junction containing two internal loops on adjoining arms. These two loops interact in a cation-driven docking step prior to chemical catalysis to form a tightly integrated structure, with dramatic changes occurring in the conformation of each loop upon docking. We investigate the thermodynamics and kinetics of the docking process using constructs in which loop A and loop B reside on separate molecules. Using a novel CD difference assay to isolate the effects of metal ions linked to domain docking, we find the intermolecular docking process to be driven by sub-millimolar concentrations of the exchange-inert Co(NH₃)₆³⁺. RNA self-cleavage requires binding of lower-affinity ions with greater apparent cooperativity than the docking process itself, implying that, even in the absence of direct coordination to RNA, metal ions play a catalytic role in hairpin ribozyme function beyond simply driving loop-loop docking. Surface plasmon resonance assays reveal remarkably slow molecular association, given the relatively tight loop-loop interaction. This observation is consistent with a “double conformational capture” model in which only collisions between loop A and loop B molecules that are simultaneously in minor, docking-competent conformations are productive for binding.

A transition-state interaction shifts nucleobase ionization toward neutrality to facilitate small ribozyme catalysis

One mechanism by which ribozymes can accelerate biological reactions is by adopting folds that favorably perturb nucleobase ionization. Herein we used Raman crystallography to directly measure pK(a) values for the Ade38 N1 imino group of a hairpin ribozyme in distinct conformational states. A transition-state analogue gave a pK(a) value of 6.27 ± 0.05, which agrees strikingly well with values measured by pH-rate analyses. To identify the chemical attributes that contribute to the shifted pK(a), we determined crystal structures of hairpin ribozyme variants containing single-atom substitutions at the active site and measured their respective Ade38 N1 pK(a) values. This approach led to the identification of a single interaction in the transition-state conformation that elevates the base pK(a) > 0.8 log unit relative to the precatalytic state. The agreement of the microscopic and macroscopic pK(a) values and the accompanying structural analysis supports a mechanism in which Ade38 N1(H)+ functions as a general acid in phosphodiester bond cleavage. Overall the results quantify the contribution of a single electrostatic interaction to base ionization, which has broad relevance for understanding how RNA structure can control chemical reactivity.

Crystallographic analysis of small ribozymes and riboswitches

Ribozymes and riboswitches are RNA motifs that accelerate biological reactions and regulate gene expression in response to metabolite recognition, respectively. These RNA molecules gain functionality via complex folding that cannot be predicted a priori, and thus requires high-resolution three-dimensional structure determination to locate key functional attributes. Herein, we present an overview of the methods used to determine small RNA structures with an emphasis on RNA preparation, crystallization, and structure refinement. We draw upon examples from our own research in the analysis of the leadzyme ribozyme, the hairpin ribozyme, a class I preQ(1) riboswitch, and variants of a larger class II preQ(1) riboswitch. The methods presented provide a guide for comparable investigations of noncoding RNA molecules including a 48-solution, "first choice" RNA crystal screen compiled from our prior successes with commercially available screens.

Multistrand RNA secondary structure prediction and nanostructure design including pseudoknots

We are presenting NanoFolder, a method for the prediction of the base pairing of potentially pseudoknotted multistrand RNA nanostructures. We show that the method outperforms several other structure prediction methods when applied to RNA complexes with non-nested base pairs. We extended this secondary structure prediction capability to allow RNA sequence design. Using native PAGE, we experimentally confirm that four in silico designed RNA strands corresponding to a triangular RNA structure form the expected stable complex.

QM/MM studies of hairpin ribozyme self-cleavage suggest the feasibility of multiple competing reaction mechanisms

The hairpin ribozyme is a prominent member of small ribozymes since it does not require metal ions to achieve catalysis. Guanine 8 (G8) and adenine 38 (A38) have been identified as key participants in self-cleavage and -ligation. We have carried out hybrid quantum-mechanical/molecular mechanical (QM/MM) calculations to evaluate the energy along several putative reaction pathways. The error of our DFT description of the QM region was tested and shown to be ~1 kcal/mol. We find that self-cleavage of the hairpin ribozyme may follow several competing microscopic reaction mechanisms, all with calculated activation barriers in good agreement with those from experiment (20-21 kcal/mol). The initial nucleophilic attack of the A-1(2'-OH) group on the scissile phosphate is predicted to be rate-limiting in all these mechanisms. An unprotonated G8(-) (together with A38H(+)) yields a feasible activation barrier (20.4 kcal/mol). Proton transfer to a nonbridging phosphate oxygen also leads to feasible reaction pathways. Finally, our calculations consider thio-substitutions of one or both nonbridging oxygens of the scissile phosphate and predict that they have only a negligible effect on the reaction barrier, as observed experimentally.

Two distinct catalytic strategies in the hepatitis δ virus ribozyme cleavage reaction

The hepatitis delta virus (HDV) ribozyme and related RNAs are widely dispersed in nature. This RNA is a small nucleolytic ribozyme that self-cleaves to generate products with a 2',3'-cyclic phosphate and a free 5'-hydroxyl. Although small ribozymes are dependent on divalent metal ions under biologically relevant buffer conditions, they function in the absence of divalent metal ions at high ionic strengths. This characteristic suggests that a functional group within the covalent structure of small ribozymes is facilitating catalysis. Structural and mechanistic analyses have demonstrated that the HDV ribozyme active site contains a cytosine with a perturbed pK(a) that serves as a general acid to protonate the leaving group. The reaction of the HDV ribozyme in monovalent cations alone never approaches the velocity of the Mg(2+)-dependent reaction, and there is significant biochemical evidence that a Mg(2+) ion participates directly in catalysis. A recent crystal structure of the HDV ribozyme revealed that there is a metal binding pocket in the HDV ribozyme active site. Modeling of the cleavage site into the structure suggested that this metal ion can interact directly with the scissile phosphate and the nucleophile. In this manner, the Mg(2+) ion can serve as a Lewis acid, facilitating deprotonation of the nucleophile and stabilizing the conformation of the cleavage site for in-line attack of the nucleophile at the scissile phosphate. This catalytic strategy had previously been observed only in much larger ribozymes. Thus, in contrast to most large and small ribozymes, the HDV ribozyme uses two distinct catalytic strategies in its cleavage reaction.

Noncanonical hydrogen bonding in nucleic acids. Benchmark evaluation of key base-phosphate interactions in folded RNA molecules using quantum-chemical calculations and molecular dynamics simulations

RNA molecules are stabilized by a wide range of noncanonical interactions that are not present in DNA. Among them, the recently classified base-phosphate (BPh) interactions belong to the most important ones. Twelve percent of nucleotides in the ribosomal crystal structures are involved in BPh interactions. BPh interactions are highly conserved and provide major constraints on RNA sequence evolution. Here we provide assessment of the energetics of BPh interactions using MP2 computations extrapolated to the complete basis set of atomic orbitals and corrected for higher-order electron correlation effects. The reference computations are compared with DFT-D and DFT-D3 approaches, the SAPT method, and the molecular mechanics force field. The computations, besides providing the basic benchmark for the BPh interactions, allow some refinements of the original classification, including identification of some potential doubly bonded BPh patterns. The reference computations are followed by analysis of some larger RNA fragments that consider the context of the BPh interactions. The computations demonstrate the complexity of interaction patterns utilizing the BPh interactions in real RNA structures. The BPh interactions are often involved in intricate interaction networks. We studied BPh interactions of protonated adenine that can contribute to catalysis of hairpin ribozyme, the key BPh interaction in the S-turn motif of the sarcin-ricin loop, which may predetermine the S-turn topology and complex BPh patterns from the glmS riboswitch. Finally, the structural stability of BPh interactions in explicit solvent molecular dynamics simulations is assessed. The simulations well preserve key BPh interactions and allow dissection of structurally/functionally important water-meditated BPh bridges, which could not be considered in earlier bioinformatics classification of BPh interactions.

Base ionization and ligand binding: how small ribozymes and riboswitches gain a foothold in a protein world

Genome sequencing has produced thousands of nonprotein coding (nc)RNA sequences including new ribozymes and riboswitches. Such RNAs are notable for their extraordinary functionality, which entails exquisite folding that culminates in biocatalytic or ligand-binding capabilities. Here we discuss advances in relating ncRNA form to function with an emphasis on base pK(a) shifting by the hairpin and hepatitis delta virus ribozymes. We then describe ligand binding by the two smallest riboswitches, which target preQ(1) and S-adenosyl-(l)-homocysteine, followed by an analysis of a second-messenger riboswitch that binds cyclic-di-GMP. Each riboswitch is then compared to a protein that binds the same ligand to contrast binding properties. The results showcase the breadth of functionality attainable from ncRNAs, as well as molecular features notable for antibacterial design.

The pH dependence of hairpin ribozyme catalysis reflects ionization of an active site adenine

Understanding how self-cleaving ribozymes mediate catalysis is crucial in light of compelling evidence that human and bacterial gene expression can be regulated through RNA self-cleavage. The hairpin ribozyme catalyzes reversible phosphodiester bond cleavage through a mechanism that does not require divalent metal cations. Previous structural and biochemical evidence implicated the amidine group of an active site adenosine, A38, in a pH-dependent step in catalysis. We developed a way to determine microscopic pK(a) values in active ribozymes based on the pH-dependent fluorescence of 8-azaadenosine (8azaA). We compared the microscopic pK(a) for ionization of 8azaA at position 38 with the apparent pK(a) for the self-cleavage reaction in a fully functional hairpin ribozyme with a unique 8azaA at position 38. Microscopic and apparent pK(a) values were virtually the same, evidence that A38 protonation accounts for the decrease in catalytic activity with decreasing pH. These results implicate the neutral unprotonated form of A38 in a transition state that involves formation of the 5'-oxygen-phosphorus bond.

The glmS ribozyme: use of a small molecule coenzyme by a gene-regulatory RNA

The glmS ribozyme is the first known example of a natural ribozyme that has evolved to require binding of an exogenous small molecule for activity. In Gram-positive bacteria, this RNA domain is part of the messenger RNA (mRNA) encoding the essential enzyme that synthesizes glucosamine-6-phosphate (GlcN6P). When present at physiologic concentration, this small molecule binds to the glmS ribozyme and uncovers a latent self-cleavage activity that ultimately leads to degradation of the mRNA. Biochemical and structural studies reveal that the RNA adopts a rigid fold stabilized by three pseudoknots and the packing of a peripheral domain against the ribozyme core. GlcN6P binding to this pre-organized RNA does not induce conformational changes; rather, the small molecule functions as a coenzyme, providing a catalytically essential amine group to the active site. The ribozyme is not a passive player, however. Active site functional groups are essential for catalysis, even in the presence of GlcN6P. In addition to being a superb experimental system with which to analyze how RNA catalysts can exploit small molecule coenzymes to broaden their chemical versatility, the presence of the glmS ribozyme in numerous pathogenic bacteria make this RNA an attractive target for the development of new antibiotics and antibacterial strategies.

Use of the spliceosomal protein U1A to facilitate crystallization and structure determination of complex RNAs

The structure determination of complex RNA molecules such as ribozymes, riboswitches and aptamers by X-ray crystallography hinges on the preparation of well-ordered crystals. Success usually results from molecular engineering to facilitate crystallization. An approach that has resulted in 10 new RNA structures in the past decade is the use of the U1A crystallization module. In this approach, the cognate site for the U1A spliceosomal protein is introduced into a functionally dispensable location in the RNA of interest, and the RNA is cocrystallized with the basic RNA-binding protein. In addition to facilitating crystallization, the presence of U1A can be useful for de novo phase determination. In this paper, some general considerations for the use of this approach to RNA crystallization are presented, and specifics of the application of the U1A module to the crystallization of the hairpin ribozyme and the tetracycline aptamer are reviewed.

Investigation of the catalytic mechanism of a synthetic DNAzyme with protein-like functionality: an RNaseA mimic?

The protein enzyme ribonuclease A (RNaseA) cleaves RNA with catalytic perfection, although with little sequence specificity, by a divalent metal ion (M(2+))-independent mechanism in which a pair of imidazoles provides general acid and base catalysis, while a cationic amine provides electrostatic stabilization of the transition state. Synthetic imitation of this remarkable organo-catalyst ("RNaseA mimicry") has been a longstanding goal in biomimetic chemistry. The 9(25)-11 DNAzyme contains synthetically modified nucleotides presenting both imidazole and cationic amine side chains, and catalyzes RNA cleavage with turnover in the absence of M(2+) similarly to RNaseA. Nevertheless, the catalytic roles, if any, of the "protein-like" functional groups have not been defined, and hence the question remains whether 9(25)-11 engages any of these functionalities to mimic aspects of the mechanism of RNaseA. To address this question, we report a mechanistic investigation of 9(25)-11 catalysis wherein we have employed a variety of experiments, such as DNAzyme functional group deletion, mechanism-based affinity labeling, and bridging and nonbridging phosphorothioate substitution of the scissile phosphate. Several striking parallels exist between the results presented here for 9(25)-11 and the results of analogous experiments applied previously to RNaseA. Specifically, our results implicate two particular imidazoles in general acid and base catalysis and suggest that a specific cationic amine stabilizes the transition state via diastereoselective interaction with the scissile phosphate. Overall, 9(25)-11 appears to meet the minimal criteria of an RNaseA mimic; this demonstrates how added synthetic functionality can expand the mechanistic repertoire available to a synthetic DNA-based catalyst.

Catalytic importance of a protonated adenosine in the hairpin ribozyme active site

The hairpin ribozyme accelerates the rate of phosphodiester transfer reactions by at least 5 orders of magnitude. To achieve this rate enhancement, the active site forms via a substrate helix docking event that constrains the scissile phosphate linkage and positions G8 and A38 for catalysis, both of which have been implicated as sites of proton transfer in general acid-base catalysis. To investigate the functional groups required for hairpin activity, we previously reported a series of nucleotide analogue interference mapping experiments [Ryder, S. P., et al. (2001) RNA 7, 1454-1463]. The critical functional groups implicated in those studies were largely consistent with subsequent X-ray crystal structures, but the lack of A38 interference with 8-azaadenosine (n(8)A), a pK(a) perturbed nucleotide analogue, argued against functional base ionization at this site. This is inconsistent with a transition state crystal structure and other biochemical studies. To address this discrepancy, we investigated the hairpin ribozyme with an expanded set of pK(a) perturbed adenosine analogues containing fluorine. A38 was the only site that showed persistent and strong interference with low pK(a) analogues across a variety of construct/substrate pairs. This interference pattern suggests that A38 base ionization is required for catalytic activity. The lack of n(8)A interference at A38, in spite of its reduced pK(a), likely results from n(8)A stabilization of the docked state, which requires an unusual syn glycosidic base conformation at A38 for active site assembly. The fluorinated adenosine analogues are better suited to identify sites of functional ionization in systems where structural rearrangements are closely coupled to catalytic steps. All pK(a) reduced analogues, including those of the previous study, produce selective interference at A38 when substrates are stably bound and docked, consistent with the importance of base ionization at this site.

Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme

The hairpin ribozyme is a prominent member of the group of small catalytic RNAs (RNA enzymes or ribozymes) because it does not require metal ions to achieve catalysis. Biochemical and structural data have implicated guanine 8 (G8) and adenine 38 (A38) as catalytic participants in cleavage and ligation catalyzed by the hairpin ribozyme, yet their exact role in catalysis remains disputed. To gain insight into dynamics in the active site of a minimal self-cleaving hairpin ribozyme, we have performed extensive classical, explicit-solvent molecular dynamics (MD) simulations on time scales of 50-150 ns. Starting from the available X-ray crystal structures, we investigated the structural impact of the protonation states of G8 and A38, and the inactivating A-1(2'-methoxy) substitution employed in crystallography. Our simulations reveal that a canonical G8 agrees well with the crystal structures while a deprotonated G8 profoundly distorts the active site. Thus MD simulations do not support a straightforward participation of the deprotonated G8 in catalysis. By comparison, the G8 enol tautomer is structurally well tolerated, causing only local rearrangements in the active site. Furthermore, a protonated A38H(+) is more consistent with the crystallography data than a canonical A38. The simulations thus support the notion that A38H(+) is the dominant form in the crystals, grown at pH 6. In most simulations, the canonical A38 departs from the scissile phosphate and substantially perturbs the structures of the active site and S-turn. Yet, we occasionally also observe formation of a stable A-1(2'-OH)...A38(N1) hydrogen bond, which documents the ability of the ribozyme to form this hydrogen bond, consistent with a potential role of A38 as general base catalyst. The presence of this hydrogen bond is, however, incompatible with the expected in-line attack angle necessary for self-cleavage, requiring a rapid transition of the deprotonated 2'-oxyanion to a position more favorable for in-line attack after proton transfer from A-1(2'-OH) to A38(N1). The simulations revealed a potential force field artifact, occasional but irreversible formation of "ladder-like", underwound A-RNA structure in one of the external helices. Although it does not affect the catalytic center of the hairpin ribozyme, further studies are under way to better assess possible influence of such force field behavior on long RNA simulations.

Direct Raman measurement of an elevated base pKa in the active site of a small ribozyme in a precatalytic conformation

Catalytic RNA molecules can achieve rate acceleration by shifting base pK(a) values toward neutrality. Prior evidence has suggested that base A38 of the hairpin ribozyme plays an important role in phosphoryl transfer, possibly functioning as a general acid, or by orienting a specific water molecule for proton transfer. To address the role of A38, we used Raman spectroscopy to measure directly the pK(a) of the N1-imino moiety in the context of hairpin ribozyme crystals representative of a "precatalytic" conformation. The results revealed that the pK(a) of A38 is shifted to 5.46 +/- 0.05 relative to 3.68 +/- 0.06 derived from a reference solution of the nucleotide AMP. The elevated pK(a) correlates well with the first titration point of the macroscopic pH-rate profile of the hairpin ribozyme in solution and strongly supports A38 as a general acid catalyst in bond scission. The results confirm that A38 is protonated before the transition state, which would promote phosphorane development. Overall, the results establish a cogent structure-function paradigm that expands our understanding of how RNA structure can enhance nucleobase reactivity to catalyze biological reactions.

Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.

Exploring ribozyme conformational changes with X-ray crystallography

Relating three-dimensional fold to function is a central challenge in RNA structural biology. Toward this goal, X-ray crystallography has long been considered the "gold standard" for structure determinations at atomic resolution, although NMR spectroscopy has become a powerhouse in this arena as well. In the area of dynamics, NMR remains the dominant technique to probe the magnitude and timescales of molecular motion. Although the latter area remains largely unassailable by conventional crystallographic methods, inroads have been made on proteins using Laue radiation on timescales of ms to ns. Proposed 'fourth generation' radiation sources, such as free-electron X-ray lasers, promise ps- to fs-timescale resolution, and credible evidence is emerging that supports the feasibility of single molecule imaging. At present however, the preponderance of RNA structural information has been derived from timescale and motion insensitive crystallographic techniques. Importantly, developments in computing, automation and high-flux synchrotron sources have propelled the rapidity of 'conventional' RNA crystal structure determinations to timeframes of hours once a suitable set of phases is obtained. With a sufficient number of crystal structures, it is possible to create a structural ensemble that can provide insight into global and local molecular motion characteristics that are relevant to biological function. Here we describe techniques to explore conformational changes in the hairpin ribozyme, a representative non-protein-coding RNA catalyst. The approaches discussed include: (i) construct choice and design using prior knowledge to improve X-ray diffraction; (ii) recognition of long-range conformational changes and (iii) use of single-base or single-atom changes to create ensembles. The methods are broadly applicable to other RNA systems.

Single-atom imino substitutions at A9 and A10 reveal distinct effects on the fold and function of the hairpin ribozyme catalytic core

The hairpin ribozyme cleaves a phosphodiester bond within a cognate substrate. Structural and biochemical data indicate the conserved A9 and A10 bases reside close to the scissile bond but make distinct contributions to catalysis. To investigate these residues, we replaced the imino moiety of each base with N1-deazaadenosine. This single-atom change resulted in an 8-fold loss in k(obs) for A9 and displacement of the base from the active site; no effects were observed for A10. We propose that the imino moiety of A9 promotes a key water-mediated contact that favors transition-state formation, which suggests an enhanced chemical repertoire for RNA.

Comparative enzymology and structural biology of RNA self-cleavage

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

Identification of an imino group indispensable for cleavage by a small ribozyme

The hairpin ribozyme is a small, noncoding RNA (ncRNA) that catalyzes a site-specific phosphodiester bond cleavage reaction. Prior biochemical and structural analyses pinpointed the amidine moiety of base Ade38 as a key functional group in catalysis, but base changes designed to probe function resulted in localized misfolding of the active site. To define the requirements for chemical activity using a conservative modification, we synthesized and incorporated N1-deazaadenosine into the full-length ribozyme construct. This single-atom variant severely impairs activity, although the active-site fold remains intact in the accompanying crystal structures. The results demonstrate the essentiality of the imino moiety as well as the importance of its interaction with the substrate in the precatalytic and transition-state conformations. This work demonstrates the efficacy of single-atom approaches in the analysis of ncRNA structure-function relationships.

Classification and energetics of the base-phosphate interactions in RNA

Structured RNA molecules form complex 3D architectures stabilized by multiple interactions involving the nucleotide base, sugar and phosphate moieties. A significant percentage of the bases in structured RNA molecules in the Protein Data Bank (PDB) hydrogen-bond with phosphates of other nucleotides. By extracting and superimposing base-phosphate (BPh) interactions from a reduced-redundancy subset of 3D structures from the PDB, we identified recurrent phosphate-binding sites on the RNA bases. Quantum chemical calculations were carried out on model systems representing each BPh interaction. The calculations show that the centers of each cluster obtained from the structure superpositions correspond to energy minima on the potential energy hypersurface. The calculations also show that the most stable phosphate-binding sites occur on the Watson-Crick edge of guanine and the Hoogsteen edge of cytosine. We modified the 'Find RNA 3D' (FR3D) software suite to automatically find and classify BPh interactions. Comparison of the 3D structures of the 16S and 23S rRNAs of Escherichia coli and Thermus thermophilus revealed that most BPh interactions are phylogenetically conserved and they occur primarily in hairpin, internal or junction loops or as part of tertiary interactions. Bases that form BPh interactions, which are conserved in the rRNA 3D structures are also conserved in homologous rRNA sequence alignments.

Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme

Despite numerous structural and biochemical investigations, the catalytic mechanism of hairpin ribozyme self-cleavage remains elusive. To gain insight into the coupling of active site dynamics with activity of this small catalytic RNA, we analyzed a total of approximately 300 ns of molecular dynamics (MD) simulations. Our simulations predict improved global stability for an in vitro selected "gain of function" mutation, which is validated by native gel electrophoretic mobility shift assay. We observe that active site nucleobases and water molecules stabilize a geometry favorable to catalysis through a dynamic hydrogen bonding network. Simulations in which A38 is unprotonated show its N1 move into close proximity of the active site 2'-OH, indicating that A38 may act as a general base during cleavage, a role that has generally been discounted due to the longer distances observed in crystal structures involving inactivating substrate analogs. By contrast, simulations in which N1 of A38 is protonated place N1 in close proximity to the 5'-oxygen leaving group, which supports the proposal that A38 serves as a general acid. In analogy to protein enzymes, we discuss a plausible mechanism in which A38 acts bifunctionally and shuttles a proton directly from the 2'-OH to the 5'-oxygen. Furthermore, our simulations suggest an important role for protonation of N1 of A38 in promoting a favorable geometry similar to that observed in transition-state analog crystal structures, and support previously proposed roles of A38, G8, and long residency water molecules in transition-state stabilization.

Direct measurement of the ionization state of an essential guanine in the hairpin ribozyme

Active site guanines are critical for self-cleavage reactions of several ribozymes, but their precise functions in catalysis are unclear. To learn whether protonated or deprotonated forms of guanine predominate in the active site, microscopic pKa values were determined for ionization of 8-azaguanosine substituted for G8 in the active site of a fully functional hairpin ribozyme in order to determine microscopic pKa values for 8-azaguanine deprotonation from the pH dependence of fluorescence. Microscopic pKa values above 9 for deprotonation of 8-azaguanine in the active site were about 3 units higher than apparent pKa values determined from the pH dependence of self-cleavage kinetics. Thus, the increase in activity with increasing pH does not correlate with deprotonation of G8, and most of G8 is protonated at neutral pH. These results do not exclude a role in proton transfer, but a simple interpretation is that G8 functions in the protonated form, perhaps by donating hydrogen bonds.

The Small Ribozymes: Common and Diverse Features Observed through the FRET Lens

The hammerhead, hairpin, HDV, VS and glmS ribozymes are the five known, naturally occurring catalytic RNAs classified as the "small ribozymes". They share common reaction chemistry in cleaving their own backbone by phosphodiester transfer, but are diverse in their secondary and tertiary structures, indicating that Nature has found at least five independent solutions to a common chemical task. Fluorescence resonance energy transfer (FRET) has been extensively used to detect conformational changes in these ribozymes and dissect their reaction pathways. Common and diverse features are beginning to emerge that, by extension, highlight general biophysical properties of non-protein coding RNAs.

Structural effects of nucleobase variations at key active site residue Ade38 in the hairpin ribozyme

The hairpin ribozyme requires functional groups from Ade38 to achieve efficient bond cleavage or ligation. To identify molecular features that contribute to catalysis, structures of position 38 base variants 2,6-diaminopurine (DAP), 2-aminopurine (AP), cytosine (Cyt), and guanine (Gua) were determined between 2.2 and 2.8 A resolution. For each variant, two substrate modifications were compared: (1) a 2'-O-methyl-substituent at Ade-1 was used in lieu of the nucleophile to mimic the precatalytic state, and (2) a 3'-deoxy-2',5'-phosphodiester linkage between Ade-1 and Gua+1 was used to mimic a reaction-intermediate conformation. While the global fold of each variant remained intact, the results revealed the importance of Ade38 N1 and N6 groups. Absence of N6 resulting from AP38 coincided with failure to localize the precatalytic scissile phosphate. Cyt38 severely impaired catalysis in a prior study, and its structures here indicated an anti base conformation that sequesters the imino moiety from the scissile bond. Gua38 was shown to be even more deleterious to activity. Although the precatalytic structure was nominally affected, the reaction-intermediate conformation indicated a severe electrostatic clash between the Gua38 keto oxygen and the pro-Rp oxygen of the scissile bond. Overall, position 38 modifications solved in the presence of 2'-OMe Ade-1 deviated from in-line geometry, whereas variants with a 2',5' linkage exhibited S-turn destabilization, as well as base conformational changes from syn to anti. These findings demonstrate the importance of the Ade38 Watson-Crick face in attaining a reaction-intermediate state and the sensitivity of the RNA fold to restructuring when electrostatic and shape features fail to complement.

Electrostatic interactions in the hairpin ribozyme account for the majority of the rate acceleration without chemical participation by nucleobases

Molecular dynamics simulations using a combined quantum mechanical/molecular mechanical potential are used to determine the two-dimensional free energy profiles for the mechanism of RNA transphosphorylation in solution and catalyzed by the hairpin ribozyme. A mechanism is explored whereby the reaction proceeds without explicit chemical participation by conserved nucleobases in the active site. The ribozyme lowers the overall free energy barrier by up to 16 kcal/mol, accounting for the majority of the observed rate enhancement. The barrier reduction in this mechanism is achieved mainly by the electrostatic environment provided by the ribozyme without recruitment of active site nucleobases as acid or base catalysts. The results establish a baseline mechanism that invokes only the solvation and specific hydrogen-bonding interactions present in the ribozyme active site and provide a departure point for the exploration of alternate mechanisms where nucleobases play an active chemical role.

Shared traits on the reaction coordinates of ribonuclease and an RNA enzyme

Reaction-intermediate analogs have been used to understand how phosphoryl transfer enzymes promote catalysis. Herein we report the first structure of a small ribozyme crystallized with a 3'-OH, 2',5'-linkage in lieu of the normal phosphodiester substrate. The new structure shares features of the reaction coordinate exhibited in prior ribozyme structures including a vanadate complex that mimicked the oxyphosphorane transition state. As such, the structure exhibits reaction-intermediate traits that allow direct comparison of stabilizing interactions to the 3'-OH, 2',5'-linkage contributed by the RNA enzyme and its protein counterpart, ribonuclease. Clear similarities are observed between the respective structures including hydrogen bonds to the non-bridging oxygens of the scissile phosphate. Other commonalities include carefully poised water molecules that may alleviate charge build-up in the transition state and placement of a positive charge near the leaving group. The advantages of 2',5'-linkages to investigate phosphoryl-transfer reactions are discussed, and argue for their expanded use in structural studies.

Quantum mechanical/molecular mechanical simulation study of the mechanism of hairpin ribozyme catalysis

The molecular mechanism of hairpin ribozyme catalysis is studied with molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential with a recently developed semiempirical AM1/d-PhoT model for phosphoryl transfer reactions. Simulations are used to derive one- and two-dimensional potentials of mean force to examine specific reaction paths and assess the feasibility of proposed general acid and base mechanisms. Density-functional calculations of truncated active site models provide complementary insight to the simulation results. Key factors utilized by the hairpin ribozyme to enhance the rate of transphosphorylation are presented, and the roles of A38 and G8 as general acid and base catalysts are discussed. The computational results are consistent with available experimental data, provide support for a general acid/base mechanism played by functional groups on the nucleobases, and offer important insight into the ability of RNA to act as a catalyst without explicit participation by divalent metal ions.

Ribozyme catalysis revisited: is water involved?

Enzymatic catalysis by RNA was discovered 25 years ago, yet mechanistic insights are emerging only slowly. Thought to be metalloenzymes at first, some ribozymes proved more versatile than anticipated when shown to utilize their own functional groups for catalysis. Recent evidence suggests that some may also judiciously place structural water molecules to shuttle protons in acid-base catalyzed reactions.

Focus on function: single molecule RNA enzymology

The ability of RNA to catalyze chemical reactions was first demonstrated 25 years ago with the discovery that group I introns and RNase P function as RNA enzymes (ribozymes). Several additional ribozymes were subsequently identified, most notably the ribosome, followed by intense mechanistic studies. More recently, the introduction of single molecule tools has dissected the kinetic steps of several ribozymes in unprecedented detail and has revealed surprising heterogeneity not evident from ensemble approaches. Still, many fundamental questions of how RNA enzymes work at the molecular level remain unanswered. This review surveys the current status of our understanding of RNA catalysis at the single molecule level and discusses the existing challenges and opportunities in developing suitable assays.

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

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

Essential role of an active-site guanine in glmS ribozyme catalysis

The glmS ribozyme is a catalytic riboswitch that is activated for endonucleolytic cleavage by the coenzyme glucosamine-6-phosphate. Using kinetic assays and X-ray crystallography, we identify an active-site mutation of a conserved guanine that abolishes catalysis without perturbing coenzyme binding. Our results provide evidence that coenzyme function requires a specific nucleobase to interact with the nucleophile of the cleavage reaction.