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

Chemo-Enzymatic Generation of Highly Fluorescent Nucleoside Analogs Using Purine-Nucleoside Phosphorylase

Chemo-enzymatic syntheses of strongly fluorescent nucleoside analogs, potentially applicable in analytical biochemistry and cell biology are reviewed. The syntheses and properties of fluorescent ribofuranosides of several purine, 8-azapurine, and etheno-purine derivatives, obtained using various types of purine nucleoside phosphorylase (PNP) as catalysts, as well as α-ribose-1-phosphate (r1P) as a second substrate, are described. In several instances, the ribosylation sites are different to the canonical purine N9. Some of the obtained ribosides show fluorescence yields close to 100%. Possible applications of the new analogs include assays of PNP, nucleoside hydrolases, and other enzyme activities both in vitro and within living cells using fluorescence microscopy.

Calculated pKa Values for a Series of Aza- and Deaza-Modified Nucleobases

A variety of synthetic modified nucleobases have been used to investigate the structure and function of RNA and DNA or act as enzyme inhibitors. A set of these modifications involves the addition or removal of a nitrogen atom in the ring. These aza and deaza modifications have garnered interest as useful biochemical tools, but information on some of their physical characteristics is lacking. In this study, the B3LYP density functional with the 6-31+G(d,p) basis set and an implicit-explicit solvent model was used to perform ab initio quantum mechanical studies to estimate pK_a values of aza- and deaza-modified nucleobases. A comparison between theoretical and known experimental pK_a values was carried out, and adjustment factors were applied to 57 pK_a values in the purine and pyrimidine data sets.

Atypical structures of GAA/TTC trinucleotide repeats underlying Friedreich's ataxia: DNA triplexes and RNA/DNA hybrids

Expansion of the GAA/TTC repeats in the first intron of the FXN gene causes Friedreich's ataxia. Non-canonical structures are linked to this expansion. DNA triplexes and R-loops are believed to arrest transcription, which results in frataxin deficiency and eventual neurodegeneration. We present a systematic in silico characterization of the possible DNA triplexes that could be assembled with GAA and TTC strands; the two hybrid duplexes [r(GAA):d(TTC) and d(GAA):r(UUC)] in an R-loop; and three hybrid triplexes that could form during bidirectional transcription when the non-template DNA strand bonds with the hybrid duplex (collapsed R-loops, where the two DNA strands remain antiparallel). For both Y·R:Y and R·R:Y DNA triplexes, the parallel third strand orientation is more stable; both parallel and antiparallel protonated d(GA+A)·d(GAA):d(TTC) triplexes are stable. Apparent contradictions in the literature about the R·R:Y triplex stability is probably due to lack of molecular resolution, since shifting the third strand by a single nucleotide alters the stability ranking. In the collapsed R-loops, antiparallel d(TTC+)·d(GAA):r(UUC) is unstable, while parallel d(GAA)·r(GAA):d(TTC) and d(GA+A)·r(GAA):d(TTC) are stable. In addition to providing new structural perspectives for specific therapeutic aims, our results contribute to a systematic structural basis for the emerging field of quantitative R-loop biology.

Excited-State Proton Transfer in 8-Azapurines I: A Kinetic Analysis of 8-Azaxanthine Fluorescence

A super-continuum white laser with a half-pulse width of ~75 ps was used to observe the kinetics of a postulated excited-state proton transfer in 8-azaxanthine and its 8-methyl derivative. Both compounds exhibited dual emissions in weakly acidified alcoholic media, but only one band was present in aqueous solutions, exhibiting an abnormal Stokes shift (>12,000 cm⁻¹). It was shown that long-wavelength emissions were delayed relative to the excitation pulse within alcoholic media. The rise time was calculated to be 0.4–0.5 ns in both methanol and deuterated methanol. This is equal to the main component of the fluorescence decay in the short-wavelength band (340 nm). Time-resolved emission spectra (TRES) indicated a two-state photo-transformation model in both compounds. Global analysis of the time dependence revealed three exponential components in each compound, one of which had an identical rise-time, with the second attributed to a long-wavelength band decay (6.4 ns for aza-xanthine and 8.3 ns for its 8-methyl derivative). The origin of the third, intermediate decay time (1.41 ns for aza-xanthine and 0.87 ns for 8-methyl-azaxanthine) is uncertain, but decay-associated spectra (DAS) containing both bands suggest the participation of a contact ion pair. These results confirm the model of phototautomerism proposed earlier, but the question of the anomalous isotope effect remains unsolved.

Site-Selective Ribosylation of Fluorescent Nucleobase Analogs Using Purine-Nucleoside Phosphorylase as a Catalyst: Effects of Point Mutations

Enzymatic ribosylation of fluorescent 8-azapurine derivatives, like 8-azaguanine and 2,6-diamino-8-azapurine, with purine-nucleoside phosphorylase (PNP) as a catalyst, leads to N9, N8, and N7-ribosides. The final proportion of the products may be modulated by point mutations in the enzyme active site. As an example, ribosylation of the latter substrate by wild-type calf PNP gives N7- and N8-ribosides, while the N243D mutant directs the ribosyl substitution at N9- and N7-positions. The same mutant allows synthesis of the fluorescent N7-β-d-ribosyl-8-azaguanine. The mutated form of the E. coli PNP, D204N, can be utilized to obtain non-typical ribosides of 8-azaadenine and 2,6-diamino-8-azapurine as well. The N7- and N8-ribosides of the 8-azapurines can be analytically useful, as illustrated by N7-β-d-ribosyl-2,6-diamino-8-azapurine, which is a good fluorogenic substrate for mammalian forms of PNP, including human blood PNP, while the N8-riboside is selective to the E. coli enzyme.

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.

In vitro evolution of distinct self-cleaving ribozymes in diverse environments

In vitro evolution experiments have long been used to evaluate the roles of RNA in both modern and ancient biology, and as a tool for biotechnology applications. The conditions under which these experiments have been conducted, however, do not reflect the range of cellular environments in modern biology or our understanding of chemical environments on the early earth, when the atmosphere and oceans were largely anoxic and soluble Fe²⁺ was abundant. To test the impact of environmental factors relevant to RNA's potential role in the earliest forms of life, we evolved populations of self-cleaving ribozymes in an anoxic atmosphere with varying pH in the presence of either Fe²⁺ or Mg²⁺. Populations evolved under these different conditions are dominated by different sequences and secondary structures, demonstrating global differences in the underlying fitness landscapes. Comparisons between evolutionary outcomes and catalytic activities also indicate that Mg²⁺ can readily take the place of Fe²⁺ in supporting the catalysis of RNA cleavage at neutral pH, but not at lower pH. These results highlight the importance of considering the specific environments in which functional biopolymers evolve when evaluating their potential roles in the origin of life, extant biology, or biotechnology.

8-Azapurines as isosteric purine fluorescent probes for nucleic acid and enzymatic research

The 8-azapurines, and their 7-deaza and 9-deaza congeners, represent a unique class of isosteric (isomorphic) analogues of the natural purines, frequently capable of substituting for the latter in many biochemical processes. Particularly interesting is their propensity to exhibit pH-dependent room-temperature fluorescence in aqueous medium, and in non-polar media. We herein review the physico-chemical properties of this class of compounds, with particular emphasis on the fluorescence emission properties of their neutral and/or ionic species, which has led to their widespread use as fluorescent probes in enzymology, including enzymes involved in purine metabolism, agonists/antagonists of adenosine receptors, mechanisms of catalytic RNAs, RNA editing, etc. They are also exceptionally useful fluorescent probes for analytical and clinical applications in crude cell homogenates.

Excited-state proton transfer and phototautomerism in nucleobase and nucleoside analogs: a mini-review

Intermolecular excited-state proton transfer (ESPT) has been observed in several fluorescent nucleobase and/or nucleoside analogs. In the present work, some new examples of ESPT in this class of compounds are presented together with a brief recapitulation of the previously published data. The nucleobases, nucleosides, and their analogs contain many basic and acidic centers and therefore their ESPT behavior may be complex. To interpret the complex data, it is usually necessary to determine the microscopic pK* values for each (or most) of the possible ESPT centers. Typical approach to solve this problem is by analysis of the alkyl derivatives, in which the possibility of the ESPT is reduced. Of particular interest are examples of "phototautomerization via the cation," observed in several systems, which in the neutral media do not undergo ESPT. Protonation of the molecule in the ground state facilitates the two-step phototautomerism in several systems, including formycin A and 2-amino-8-azadenine. Fluorescence of the nucleobase and nucleoside analogs undergoing ESPT is usually solvent-, isotope-, and buffer-ion sensitive, and in some systems the ESPT can be promoted by environmental factors, e.g., the presence of buffer ions. This sensitivity to the microenvironment parameters makes the ESPT systems potentially useful for biological applications.

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.

Study of pK values and effective dielectric constants of ionizable residues in pentapeptides and in staphylococcal nuclease (SNase) using a mean-field approach

The determination of pK values of amino acid residues as a function of temperature and ionic concentration is crucial to understanding the dynamics of various biological processes such as adsorption of peptides and their interactions with active sites of enzymes. In this study we developed a mean-field model to calculate the position-dependent dielectric constants of ionizable groups and the mean electrostatic potential on the surface. Such potential, which takes into account the contributions exerted by neighboring groups and ions in solution, is responsible for the fine-tuning of the pK value of each residue. The proposed model was applied to the amino acids Asp, Glu, Lys, His, Tyr, and Cys, and since the results were consistent with experimentally obtained values, the model was extended and applied to computation of pK values of Gly and Ala pentapeptides and of ionizable residues of the enzyme staphylococcal nuclease (SNase). In this latter case, we used an approach similar to a first-neighbors approximation, and the results turned out to be in good agreement with previously reported data when considering only the interactions of charged groups located at distances of maximally 20 Å. These considerations and the little computational cost involved turn the suggested approach into a promising tool for the modeling of force fields in computational simulations.

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.

Experimental approaches for measuring pKa's in RNA and DNA

RNA and DNA carry out diverse functions in biology including catalysis, splicing, gene regulation, and storage of genetic information. Interest has grown in understanding how nucleic acids perform such sophisticated functions given their limited molecular repertoire. RNA can fold into diverse shapes that often perturb pKa values and allow it to ionize appreciably under biological conditions, thereby extending its molecular diversity. The goal of this chapter is to enable experimental measurement of pKa's in RNA and DNA. A number of experimental methods for measuring pKa values in RNA and DNA have been developed over the last 10 years, including RNA cleavage kinetics; UV-, fluorescence-, and NMR-detected pH titrations; and Raman crystallography. We begin with general considerations for choosing a pKa assay and then describe experimental conditions, advantages, and disadvantages for these assays. Potential pitfalls in measuring a pKa are provided including the presence of apparent pKa's due to a kinetic pKa or coupled acid- and alkali-promoted RNA unfolding, as well as degradation of RNA, precipitation of metal hydroxides and poor baselines. Use of multiple data fitting procedures and the study of appropriate mutants are described as ways to avoid some of these pitfalls. Application of these experimental methods to RNA and DNA will increase the number of available nucleic acid pKa values in the literature, which should deepen insight into biology and provide benchmarks for pKa calculations. Future directions for measuring pKa's in nucleic acids are discussed.

Enzymatic synthesis of highly fluorescent 8-azapurine ribosides using a purine nucleoside phosphorylase reverse reaction: variable ribosylation sites

Various forms of purine-nucleoside phosphorylase (PNP) were used as catalysts of enzymatic ribosylation of selected fluorescent 8-azapurines. It was found that the recombinant calf PNP catalyzes ribosylation of 2,6-diamino-8-azapurine in a phosphate-free medium, with ribose-1-phosphate as ribose donor, but the ribosylation site is predominantly N7 and N8, with the proportion of N8/N7 ribosylated products markedly dependent on the reaction conditions. Both products are fluorescent. Application of the E. coli PNP gave a mixture of N8 and N9-substituted ribosides. Fluorescence of the ribosylated 2,6-diamino-8-azapurine has been briefly characterized. The highest quantum yield, ~0.9, was obtained for N9-β-d-riboside (λ_max 365 nm), while for N8-β-d-riboside, emitting at ~430 nm, the fluorescence quantum yield was found to be close to 0.4. Ribosylation of 8-azaguanine with calf PNP as a catalyst goes exclusively to N9. By contrast, the E. coli PNP ribosylates 8-azaGua predominantly at N9, with minor, but highly fluorescent products ribosylated at N8/N7.

Towards Accurate Prediction of Protonation Equilibrium of Nucleic Acids

The role of protonated nucleotides in modulating the pH-dependent properties of nucleic acids is one of the emerging frontiers in the field of nucleic acid biology. The recent development of a constant pH molecular dynamics simulation (CPHMD^MSλD) framework for simulating nucleic acids has provided a tool for realistic simulations of pH-dependent dynamics. We enhanced the CPHMD^MSλD framework with pH-based replica exchange (pH-REX), which significantly improves the sampling of both titration and spatial coordinates. The results from our pK_a calculations for the GAAA tetraloop, which was predicted with lower accuracy previously due to sampling challenges, demonstrates that pH-REX reduces the average unsigned error (AUE) to 0.7 pK_a units, and the error of the most poorly predicted residue A17 was drastically reduced from 2.9 to 1.2 pK_a unit. Lastly, we show that pH-REX CPHMD^MSλD simulations can be used to identify the dominant conformation of nucleic acid structures in alternate pH environments. This work suggests that pH-REX CPHMD^MSλD simulations provide a practical tool for predicting nucleic acid protonation equilibrium from first-principles, and offering structural and mechanistic insight into the study of pH-dependent properties of nucleic acids.

pH-dependent dynamics of complex RNA macromolecules

The role of pH-dependent protonation equilibrium in modulating RNA dynamics and function is one of the key unanswered questions in RNA biology. Molecular dynamics (MD) simulations can provide insight into the mechanistic roles of protonated nucleotides, but it is only capable of modeling fixed protonation states and requires prior knowledge of the key residue's protonation state. Recently, we developed a framework for constant pH molecular dynamics simulations (CPHMDMSλD) of nucleic acids, where the nucleotides' protonation states are modeled as dynamic variables that are coupled to the structural dynamics of the RNA. In the present study, we demonstrate the application of CPHMDMSλD to the lead-dependent ribozyme; establishing the validity of this approach for modeling complex RNA structures. We show that CPHMDMSλD accurately predicts the direction of the pKa shifts and reproduces experimentally-measured microscopic pKa values with an average unsigned error of 1.3 pKa units. The effects of coupled titration states in RNA structures are modeled, and the importance of conformation sampling is highlighted. The general accuracy of CPHMDMSλD simulations in reproducing pH-dependent observables reported in this work demonstrates that constant pH simulations provides a powerful tool to investigate pH-dependent processes in nucleic acids.

Lewis acid catalysis of phosphoryl transfer from a copper(II)-NTP complex in a kinase ribozyme

The chemical strategies used by ribozymes to enhance reaction rates are revealed in part from their metal ion and pH requirements. We find that kinase ribozyme K28(1-77)C, in contrast with previously characterized kinase ribozymes, requires Cu²⁺ for optimal catalysis of thiophosphoryl transfer from GTPγS. Phosphoryl transfer from GTP is greatly reduced in the absence of Cu²⁺, indicating a specific catalytic role independent of any potential interactions with the GTPγS thiophosphoryl group. In-line probing and ATPγS competition both argue against direct Cu²⁺ binding by RNA; rather, these data establish that Cu²⁺ enters the active site within a Cu²⁺•GTPγS or Cu²⁺•GTP chelation complex, and that Cu²⁺•nucleobase interactions further enforce Cu²⁺ selectivity and position the metal ion for Lewis acid catalysis. Replacing Mg²⁺ with [Co(NH₃)₆]³⁺ significantly reduced product yield, but not k_obs, indicating that the role of inner-sphere Mg²⁺ coordination is structural rather than catalytic. Replacing Mg²⁺ with alkaline earths of increasing ionic radii (Ca²⁺, Sr²⁺ and Ba²⁺) gave lower yields and approximately linear rates of product accumulation. Finally, we observe that reaction rates increased with pH in log-linear fashion with an apparent pKa = 8.0 ± 0.1, indicating deprotonation in the rate-limiting step.

The glmS ribozyme cofactor is a general acid-base catalyst

The glmS ribozyme is the first natural self-cleaving ribozyme known to require a cofactor. The d-glucosamine-6-phosphate (GlcN6P) cofactor has been proposed to serve as a general acid, but its role in the catalytic mechanism has not been established conclusively. We surveyed GlcN6P-like molecules for their ability to support self-cleavage of the glmS ribozyme and found a strong correlation between the pH dependence of the cleavage reaction and the intrinsic acidity of the cofactors. For cofactors with low binding affinities, the contribution to rate enhancement was proportional to their intrinsic acidity. This linear free-energy relationship between cofactor efficiency and acid dissociation constants is consistent with a mechanism in which the cofactors participate directly in the reaction as general acid-base catalysts. A high value for the Brønsted coefficient (β ~ 0.7) indicates that a significant amount of proton transfer has already occurred in the transition state. The glmS ribozyme is the first self-cleaving RNA to use an exogenous acid-base catalyst.

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.

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.

Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent

The nucleosides of adenine and cytosine have pKa values of 3.50 and 4.08, respectively, and are assumed to be unprotonated under physiological conditions. However, evidence from recent NMR and X-Ray crystallography studies has revealed the prevalence of protonated adenine and cytosine in RNA macromolecules. Such nucleotides with elevated pKa values may play a role in stabilizing RNA structure and participate in the mechanism of ribozyme catalysis. With the work presented here, we establish the framework and demonstrate the first constant pH MD simulations (CPHMD) for nucleic acids in explicit solvent in which the protonation state is coupled to the dynamical evolution of the RNA system via λ-dynamics. We adopt the new functional form λ(Nexp) for λ that was recently developed for Multi-Site λ-Dynamics (MSλD) and demonstrate good sampling characteristics in which rapid and frequent transitions between the protonated and unprotonated states at pH = pKa are achieved. Our calculated pKa values of simple nucleotides are in a good agreement with experimentally measured values, with a mean absolute error of 0.24 pKa units. This work demonstrates that CPHMD can be used as a powerful tool to investigate pH-dependent biological properties of RNA macromolecules.

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.

An active-site guanine participates in glmS ribozyme catalysis in its protonated state

Active-site guanines that occupy similar positions have been proposed to serve as general base catalysts in hammerhead, hairpin, and glmS ribozymes, but no specific roles for these guanines have been demonstrated conclusively. Structural studies place G33(N1) of the glmS ribozyme of Bacillus anthracis within hydrogen-bonding distance of the 2'-OH nucleophile. Apparent pK(a) values determined from the pH dependence of cleavage kinetics for wild-type and mutant glmS ribozymes do not support a role for G33, or any other active-site guanine, in general base catalysis. Furthermore, discrepancies between apparent pK(a) values obtained from functional assays and microscopic pK(a) values obtained from pH-fluorescence profiles with ribozymes containing a fluorescent guanosine analogue, 8-azaguanosine, at position 33 suggest that the pH-dependent step in catalysis does not involve G33 deprotonation. These results point to an alternative model in which G33(N1) in its neutral, protonated form donates a hydrogen bond to stabilize the transition state.

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.

Arginine cofactors on the polymerase ribozyme

The RNA world hypothesis states that the early evolution of life went through a stage in which RNA served both as genome and as catalyst. The central catalyst in an RNA world organism would have been a ribozyme that catalyzed RNA polymerization to facilitate self-replication. An RNA polymerase ribozyme was developed previously in the lab but it is not efficient enough for self-replication. The factor that limits its polymerization efficiency is its weak sequence-independent binding of the primer/template substrate. Here we tested whether RNA polymerization could be improved by a cationic arginine cofactor, to improve the interaction with the substrate. In an RNA world, amino acid-nucleic acid conjugates could have facilitated the emergence of the translation apparatus and the transition to an RNP world. We chose the amino acid arginine for our study because this is the amino acid most adept to interact with RNA. An arginine cofactor was positioned at ten different sites on the ribozyme, using conjugates of arginine with short DNA or RNA oligonucleotides. However, polymerization efficiency was not increased in any of the ten positions. In five of the ten positions the arginine reduced or modulated polymerization efficiency, which gives insight into the substrate-binding site on the ribozyme. These results suggest that the existing polymerase ribozyme is not well suited to using an arginine cofactor.

The glmS ribozyme tunes the catalytically critical pK(a) of its coenzyme glucosamine-6-phosphate

The glmS ribozyme riboswitch is the first known natural catalytic RNA that employs a small-molecule cofactor. Binding of glucosamine-6-phosphate (GlcN6P) uncovers the latent self-cleavage activity of the RNA, which adopts a catalytically competent conformation that is nonetheless inactive in the absence of GlcN6P. Structural and analogue studies suggest that the amine of GlcN6P functions as a general acid-base catalyst, while its phosphate is important for binding affinity. However, the solution pK(a) of the amine is 8.06 ± 0.05, which is not optimal for proton transfer. Here we used Raman crystallography directly to determine the pK(a)'s of GlcN6P bound to the glmS ribozyme. Binding to the RNA lowers the pK(a) of the amine of GlcN6P to 7.26 ± 0.09 and raises the pK(a) of its phosphate to 6.35 ± 0.09. Remarkably, the pK(a)'s of these two functional groups are unchanged from their values for free GlcN6P (8.06 ± 0.05 and 5.98 ± 0.05, respectively) when GlcN6P binds to the catalytically inactive but structurally unperturbed G40A mutant of the ribozyme, thus implicating the ribozyme active site guanine in pK(a) tuning. This is the first demonstration that a ribozyme can tune the pK(a) of a small-molecule ligand. Moreover, the anionic glmS ribozyme in effect stabilizes the neutral amine of GlcN6P by lowering its pK(a). This is unprecedented and illustrates the chemical sophistication of ribozyme active sites.