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

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.

Mechanistic role of nucleobases in self-cleavage catalysis of hairpin ribozyme at ambient versus high-pressure conditions

Ribozymes catalyze the site-specific self-cleavage of intramolecular phosphodiester bonds. Initially thought to act as metalloenzymes, they are now known to be functional even in the absence of divalent metal ions and specific nucleobases directly participate in the self-cleavage reaction. Here, we use extensive replica exchange molecular dynamics simulations to probe the precise mechanistic role of nucleobases by simulating precatalytic reactant and active precursor states of a hairpin ribozyme along its reaction path at ambient as well as high-pressure conditions. The results provide novel key insights into the self-cleavage of ribozymes. We find that deprotonation of the hydroxyl group is crucial and might be the penultimate step to the self-cleavage. The G8 nucleobase is found to stabilize the activated precursor into inline arrangement for facile nucleophilic attack of the scissile phosphate only after deprotonation of the hydroxyl group. The protonated A38 nucleobase, in contrast, mainly acts a proton donor to the O5'-oxygen leaving group that eventually leads to the self-cleavage. Indeed, systematic high-pressure simulations of catalytically relevant states confirm these findings and, moreover, provide support to the role of ribozymes as piezophilic biocatalysts with regard to their relevance in early life under extreme conditions in the realm of RNA world hypothesis.

The importance of pKa in an analysis of the interaction of compounds with DNA

pKa of a compound is crucial for determining the contributions of different forms of it towards overall binding with DNA. Hence it is important to use correct pKa values in DNA interaction studies. This study takes a look at the importance of pKa values to realize binding of compounds with DNA. Since pKa of a compound determined in the presence of DNA is quite different from that determined in its absence hence, presence of different forms of a compound during interaction with DNA is different from that realized if the determination of pKa is done in normal aqueous solution in absence of DNA. Hence, calculations determining contributions of different forms of a compound interacting with DNA are affected accordingly. Two simple analogues of anthracyclines, alizarin and purpurin, were used to investigate the influence DNA has on pKa values. Indeed, they were different in presence of DNA than when determined in normal aqueous solution. pKa1 for alizarin and purpurin determined in the absence and presence of calf thymus DNA were used in equations that determine contributions of two forms (neutral and anionic) towards overall binding with DNA. The study concludes that correct pKa values, determined correctly i.e. under appropriate conditions, must be used for DNA binding experiments to evaluate contributions of individual forms.

Facile synthesis of a 3-deazaadenosine phosphoramidite for RNA solid-phase synthesis

Access to 3-deazaadenosine (c³A) building blocks for RNA solid-phase synthesis represents a severe bottleneck in modern RNA research, in particular for atomic mutagenesis experiments to explore mechanistic aspects of ribozyme catalysis. Here, we report the 5-step synthesis of a c³A phosphoramidite from cost-affordable starting materials. The key reaction is a silyl-Hilbert-Johnson nucleosidation using unprotected 6-amino-3-deazapurine and benzoyl-protected 1-O-acetylribose. The novel path is superior to previously described syntheses in terms of efficacy and ease of laboratory handling.

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

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

Exploring the mode of binding of the bioflavonoid kaempferol with B and protonated forms of DNA using spectroscopic and molecular docking studies

Mode of binding of kaempferol with B and protonated form of calf thymus DNA.

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.

General acid-base catalysis mediated by nucleobases in the hairpin ribozyme

The catalytic mechanism by which the hairpin ribozyme accelerates cleavage or ligation of the phosphodiester backbone of RNA has been incompletely understood. There is experimental evidence for an important role for an adenine (A38) and a guanine (G8), and it has been proposed that these act in general acid-base catalysis. In this work we show that a large reduction in cleavage rate on substitution of A38 by purine (A38P) can be reversed by replacement of the 5'-oxygen atom at the scissile phosphate by sulfur (5'-PS), which is a much better leaving group. This is consistent with A38 acting as the general acid in the unmodified ribozyme. The rate of cleavage of the 5'-PS substrate by the A38P ribozyme increases with pH log-linearly, indicative of a requirement for a deprotonated base with a relatively high pK(a). On substitution of G8 by diaminopurine, the 5'-PS substrate cleavage rate at first increases with pH and then remains at a plateau, exhibiting an apparent pK(a) consistent with this nucleotide acting in general base catalysis. Alternative explanations for the pH dependence of hairpin ribozyme reactivity are discussed, from which we conclude that general acid-base catalysis by A38 and G8 is the simplest and most probable explanation consistent with all the experimental data.

The chemical versatility of RNA

The ability of RNA to both store genetic information and catalyse chemical reactions has led to the hypothesis that it predates DNA and proteins. While there is no doubt that RNA is capable of storing the genetic information of a primitive organism, only two classes of reactions-phosphoryl transfer and peptide bond formation-have been observed to be catalysed by RNA in nature. However, these naturally occurring ribozymes use a wide range of catalytic strategies that could be applied to other reactions. Furthermore, RNA can bind several cofactors that are used by protein enzymes to facilitate a wide variety of chemical processes. Despite its limited functional groups, these observations indicate RNA is a versatile molecule that could, in principle, catalyse the myriad reactions necessary to sustain life.

Charged nucleobases and their potential for RNA catalysis

Catalysis in living cells is carried out by both proteins and RNA. Protein enzymes have been known for over 200 years, but RNA enzymes, or "ribozymes", were discovered only 30 years ago. Developing insight into RNA enzyme mechanisms is invaluable for better understanding both extant biological catalysis as well as the primitive catalysis envisioned in an early RNA-catalyzed life. Natural ribozymes include large RNAs such as the group I and II introns; small RNAs such as the hepatitis delta virus and the hairpin, hammerhead, VS, and glmS ribozymes; and the RNA portion of the ribosome and spliceosome. RNA enzymes use many of the same catalytic strategies as protein enzymes, but do so with much simpler side chains. Among these strategies are metal ion, general acid-base, and electrostatic catalysis. In this Account, we examine evidence for participation of charged nucleobases in RNA catalysis. Our overall approach is to integrate direct measurements on catalytic RNAs with thermodynamic studies on oligonucleotide model systems. The charged amino acids make critical contributions to the mechanisms of nearly all protein enzymes. Ionized nucleobases should be critical for RNA catalysis as well. Indeed, charged nucleobases have been implicated in RNA catalysis as general acid-bases and oxyanion holes. We provide an overview of ribozyme studies involving nucleobase catalysis and the complications involved in developing these mechanisms. We also consider driving forces for perturbation of the pK(a) values of the bases. Mechanisms for pK(a) values shifting toward neutrality involve electrostatic stabilization and the addition of hydrogen bonding. Both mechanisms couple protonation with RNA folding, which we treat with a thermodynamic formalism and conceptual models. Furthermore, ribozyme reaction mechanisms can be multichannel, which demonstrates the versatility of ribozymes but makes analysis of experimental data challenging. We examine advances in measuring and analyzing perturbed pK(a) values in RNA. Raman crystallography and fluorescence spectroscopy have been especially important for pK(a) measurement. These methods reveal pK(a) values for the nucleobases A or C equal to or greater than neutrality, conferring potential histidine- and lysine/arginine-like behavior on them. Structural support for ionization of the nucleobases also exists: an analysis of RNA structures in the databases conducted herein suggests that charging of the bases is neither especially uncommon nor difficult to achieve under cellular conditions. Our major conclusions are that cationic and anionic charge states of the nucleobases occur in RNA enzymes and that these states make important catalytic contributions to ribozyme activity. We conclude by considering outstanding questions and possible experimental and theoretical approaches for further advances.

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.

Probing DNA interstrand cross-link formation by an oxidized abasic site using nonnative nucleotides

The C4'-oxidized abasic site (C4-AP) forms two types of interstrand cross-links with the adjacent nucleotides in DNA. Previous experiments revealed that dG does not react with the lesion and that formation of one type of cross-link is catalyzed by the opposing dA. iso-Guanosine·dC and 2-aminopurine·dT base pairs were used to determine why dG does not cross-link with C4-AP despite its well known reactivity with other bis-electrophiles. 7-Deaza-2'-deoxyadenosine was used to probe the role of the nucleotide opposite C4-AP in the catalysis of interstrand cross-link formation.

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.

Do the hairpin and VS ribozymes share a common catalytic mechanism based on general acid-base catalysis? A critical assessment of available experimental data

The active centers of the hairpin and VS ribozymes are both generated by the interaction of two internal loops, and both ribozymes use guanine and adenine nucleobases to accelerate cleavage and ligation reactions. The centers are topologically equivalent and the relative positioning of key elements the same. There is good evidence that the cleavage reaction of the VS ribozyme is catalyzed by the guanine (G638) acting as general base and the adenine (A756) as general acid. We now critically evaluate the experimental mechanistic evidence for the hairpin ribozyme. We conclude that all the available data are fully consistent with a major contribution to catalysis by general acid-base catalysis involving the adenine (A38) and guanine (G8). It appears that the two ribozymes are mechanistically equivalent.

Nucleotide docking: prediction of reactant state complexes for ribonuclease enzymes

Ribonuclease enzymes (RNases) play key roles in the maturation and metabolism of all RNA molecules. Computational simulations of the processes involved can help to elucidate the underlying enzymatic mechanism and is often employed in a synergistic approach together with biochemical experiments. Theoretical calculations require atomistic details regarding the starting geometries of the molecules involved, which, in the absence of crystallographic data, can only be achieved from computational docking studies. Fortunately, docking algorithms have improved tremendously in recent years, so that reliable structures of enzyme-ligand complexes can now be successfully obtained from computation. However, most docking programs are not particularly optimized for nucleotide docking. In order to assist our studies on the cleavage of RNA by the two most important ribonuclease enzymes, RNase A and RNase H, we evaluated four docking tools-MOE2009, Glide 5.5, QXP-Flo+0802, and Autodock 4.0-for their ability to simulate complexes between these enzymes and RNA oligomers. To validate our results, we analyzed the docking results with respect to the known key interactions between the protein and the nucleotide. In addition, we compared the predicted complexes with X-ray structures of the mutated enzyme as well as with structures obtained from previous calculations. In this manner, we were able to prepare the desired reaction state complex so that it could be used as the starting structure for further DFT/B3LYP QM/MM reaction mechanism studies.

Hammerhead ribozymes: true metal or nucleobase catalysis? Where is the catalytic power from?

The hammerhead ribozyme was first considered as a metalloenzyme despite persistent inconsistencies between structural and functional data. In the last decade, metal ions were confirmed as catalysts in self-splicing ribozymes but displaced by nucleobases in self-cleaving ribozymes. However, a model of catalysis just relying on nucleobases as catalysts does not fully fit some recent data. Gathering and comparing data on metal ions in self-cleaving and self-splicing ribozymes, the roles of divalent metal ions and nucleobases are revisited. Hypothetical models based on cooperation between metal ions and nucleobases are proposed for the catalysis and evolution of this prototype in RNA catalysis.