Involvement of a cytosine side chain in proton transfer in the rate-determining step of ribozyme self-cleavage

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

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

Conjugation of RNA via 2'-OH acylation: Mechanisms determining nucleotide reactivity

The acylation reactivity of RNA 2'-OH groups has proven broadly useful for labeling and mapping RNA. Here we perform kinetics studies to test the mechanisms governing this reaction, and we find strong steric and inductive effects modulating reactivity. The results shed light on new strategies for improved conjugation and mapping.

Calculations of pKa Values for a Series of Naturally Occurring Modified Nucleobases

Modified nucleobases are found in functionally important regions of RNA and are often responsible for essential structural roles. Many of these nucleobase modifications are dynamically regulated in nature, with each modification having a different biological role in RNA. Despite the high abundance of modifications, many of their characteristics are still poorly understood. One important property of a nucleobase is its pK_a value, which has been widely studied for unmodified nucleobases, but not for the modified versions. In this study, the pK_a values of modified nucleobases were determined by performing ab initio quantum mechanical calculations using a B3LYP density functional with the 6-31+G(d,p) basis set and a combination of implicit-explicit solvation systems. This method, which was previously employed to determine the pK_a values of unmodified nucleobases, is applicable to a variety of modified nucleobases. Comparisons of the pK_a values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.

Small Molecule Rescue and Glycosidic Conformational Analysis of the Twister Ribozyme

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

Cellular Small Molecules Contribute to Twister Ribozyme Catalysis

The number of self-cleaving small ribozymes has increased sharply in recent years. Advances have been made in describing these ribozymes in terms of four catalytic strategies: α describes in-line attack, β describes neutralization of the nonbridging oxygens, γ describes activation of the nucleophile, and δ describes stabilization of the leaving group. Current literature presents the rapid self-cleavage of the twister ribozyme in terms of all four of these classic catalytic strategies. Herein, we describe the nonspecific contribution of small molecules to ribozyme catalysis. At biological pH, the rate of the wild-type twister ribozyme is enhanced up to 5-fold in the presence of moderate buffer concentrations, similar to the 3-5-fold effects reported previously for buffer catalysis for protein enzymes. We observe this catalytic enhancement not only with standard laboratory buffers, but also with diverse biological small molecules, including imidazole, amino acids, and amino sugars. Brønsted plots suggest that small molecules assist in proton transfer, most likely with δ catalysis. Cellular small molecules provide a simple way to overcome the limited functional diversity of RNA and have the potential to participate in the catalytic mechanisms of many ribozymes in vivo.

Single-Molecule Titration in a Protein Nanoreactor Reveals the Protonation/Deprotonation Mechanism of a C:C Mismatch in DNA

Measurement of single-molecule reactions can elucidate microscopic mechanisms that are often hidden from ensemble analysis. Herein, we report the acid-base titration of a single DNA duplex confined within the wild-type α-hemolysin (α-HL) nanopore for up to 3 h, while monitoring the ionic current through the nanopore. Modulation between two states in the current-time trace for duplexes containing the C:C mismatch in proximity to the latch constriction of α-HL is attributed to the base flipping of the C:C mismatch. As the pH is lowered, the rate for the C:C mismatch to flip from the intra-helical state to the extra-helical state ( k_intra-extra) decreases, while the rate for base flipping from the extra-helical state to the intra-helical state ( k_extra-intra) remains unchanged. Both k_intra-extra and k_extra-intra are on the order of 1 × 10^-2 s^-1 to 1 × 10^-1 s^-1 and remain stable over the time scale of the measurement (several hours). Analysis of the pH-dependent kinetics of base flipping using a hidden Markov kinetic model demonstrates that protonation/deprotonation occurs while the base pair is in the intra-helical state. We also demonstrate that the rate of protonation is limited by transport of H⁺ into the α-HL nanopore. Single-molecule kinetic isotope experiments exhibit a large kinetic isotope effect (KIE) for k_intra-extra ( k_H/ k_D ≈ 5) but a limited KIE for k_extra-intra ( k_H/ k_D ≈ 1.3), supporting our model. Our experiments correspond to the longest single-molecule measurements performed using a nanopore, and demonstrate its application in interrogating mechanisms of single-molecule reactions in confined geometries.

The Novel Chemical Mechanism of the Twister Ribozyme

We describe the multifactorial origins of catalysis by the twister ribozyme. We provide evidence that the adenine immediately 3' to the scissile phosphate (A1) acts as a general acid. Substitution of ring nitrogen atoms indicates that very unusually the N3 of A1 is the proton donor to the oxyanion leaving group. A1 is accommodated in a specific binding pocket that raises its pKa toward neutrality, juxtaposes its N3 with the O5' to be protonated, and helps create the in-line trajectory required for nucleophilic attack. A1 performs general acid catalysis while G33 acts as a general base. A 100-fold stereospecific phosphorothioate effect at the scissile phosphate is consistent with a significant stabilization of the transition state by the ribozyme, and functional group substitution at G33 indicates that its exocyclic N2 interacts directly with the scissile phosphate. A model of the ribozyme active site is proposed that accommodates these catalytic strategies.

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

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

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.

Transition State Features in the Hepatitis Delta Virus Ribozyme Reaction Revealed by Atomic Perturbations

Endonucleolytic ribozymes constitute a class of non-coding RNAs that catalyze single-strand RNA scission. With crystal structures available for all of the known ribozymes, a major challenge involves relating functional data to the physically observed RNA architecture. In the case of the hepatitis delta virus (HDV) ribozyme, there are three high-resolution crystal structures, the product state of the reaction and two precursor variants, with distinct mechanistic implications. Here, we develop new strategies to probe the structure and catalytic mechanism of a ribozyme. First, we use double-mutant cycles to distinguish differences in functional group proximity implicated by the crystal structures. Second, we use a corrected form of the Brønsted equation to assess the functional significance of general acid catalysis in the system. Our results delineate the functional relevance of atomic interactions inferred from structure, and suggest that the HDV ribozyme transition state resembles the cleavage product in the degree of proton transfer to the leaving group.

Ammonia as an in situ sanitizer: inactivation kinetics and mechanisms of the ssRNA virus MS2 by NH3

Sanitizing human and animal waste (e.g., urine, fecal sludge, or grey water) is a critical step in reducing the spread of disease and ensuring microbially safe reuse of waste materials. Viruses are particularly persistent pathogens and can be transmitted through inadequately sanitized waste. However, adequate storage or digestion of waste can strongly reduce the number of viruses due to increases in pH and uncharged aqueous ammonia (NH3), a known biocide. In this study we investigated the kinetics and mechanisms of inactivation of the single-stranded RNA virus MS2 under temperature, pH and NH3 conditions representative of waste storage. MS2 inactivation was mainly controlled by the activity of NH3 over a pH range of 7.0–9.5 and temperatures lower than 40 °C. Other bases (e.g., hydroxide, carbonate) additionally contributed to the observed reduction of infective MS2. The loss in MS2 infectivity could be rationalized by a loss in genome integrity, which was attributed to genome cleavage via alkaline transesterification. The contribution of each base to genome transesterification, and hence inactivation, could be related to the base pKa by means of a Bronsted relationship. The Bronsted relationship in conjunction with the activity of bases in solution enabled an accurate prediction of MS2 inactivation rates.

New tools provide a second look at HDV ribozyme structure, dynamics and cleavage

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme essential for processing viral transcripts during rolling circle viral replication. The first crystal structure of the cleaved ribozyme was solved in 1998, followed by structures of uncleaved, mutant-inhibited and ion-complexed forms. Recently, methods have been developed that make the task of modeling RNA structure and dynamics significantly easier and more reliable. We have used ERRASER and PHENIX to rebuild and re-refine the cleaved and cis-acting C75U-inhibited structures of the HDV ribozyme. The results correct local conformations and identify alternates for RNA residues, many in functionally important regions, leading to improved R values and model validation statistics for both structures. We compare the rebuilt structures to a higher resolution, trans-acting deoxy-inhibited structure of the ribozyme, and conclude that although both inhibited structures are consistent with the currently accepted hammerhead-like mechanism of cleavage, they do not add direct structural evidence to the biochemical and modeling data. However, the rebuilt structures (PDBs: 4PR6, 4PRF) provide a more robust starting point for research on the dynamics and catalytic mechanism of the HDV ribozyme and demonstrate the power of new techniques to make significant improvements in RNA structures that impact biologically relevant conclusions.

HDV family of self-cleaving ribozymes

The hepatitis delta virus (HDV) ribozymes are catalytic RNAs capable of cleaving their own sugar-phosphate backbone. The HDV virus possesses the ribozymes in both sense and antisense genomic transcripts, where they are essential for processing during replication. These ribozymes have been the subject of intense biochemical scrutiny and have yielded a wealth of mechanistic insights. In recent years, many HDV-like ribozymes have been identified in nearly all branches of life. The ribozymes are implicated in a variety of biological events, including episodic memory in mammals and retrotransposition in many eukaryotes. Detailed analysis of additional HDV-like ribozyme isolates will likely reveal many more biological functions and provide information about the evolution of this unique RNA.

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein-folding, protein-protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH-dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi-site λ-dynamics (CPHMD(MSλD)). In the CPHMD(MSλD) framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi-site λ-dynamics, and designed novel biasing potentials to ensure that the physical end-states are predominantly sampled. We show that explicit solvent CPHMD(MSλD) simulations model realistic pH-dependent properties of proteins such as the Hen-Egg White Lysozyme (HEWL), binding domain of 2-oxoglutarate dehydrogenase (BBL) and N-terminal domain of ribosomal protein L9 (NTL9), and the pKa predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pKa units. With the recent development of the explicit solvent CPHMD(MSλD) framework for nucleic acids, accurate modeling of pH-dependent properties of both major class of biomolecules-proteins and nucleic acids is now possible.

Thio effects and an unconventional metal ion rescue in the genomic hepatitis delta virus ribozyme

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg(2+) ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies of the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn(2+). Here, we combine experiment and theory to explore phosphorothioate substitutions at the scissile phosphate. We report significant thio effects at the scissile phosphate and metal ion rescue with Cd(2+). Reaction profiles with an SP-phosphorothioate substitution are indistinguishable from those of the unmodified substrate in the presence of Mg(2+) or Cd(2+), supporting the idea that the pro-SP oxygen does not coordinate metal ions. The RP-phosphorothioate substitution, however, exhibits biphasic kinetics, with the fast-reacting phase displaying a thio effect of up to 5-fold and the slow-reacting phase displaying a thio effect of ~1000-fold. Moreover, the fast- and slow-reacting phases give metal ion rescues in Cd(2+) of up to 10- and 330-fold, respectively. The metal ion rescues are unconventional in that they arise from Cd(2+) inhibiting the oxo substrate but not the RP substrate. This metal ion rescue suggests a direct interaction of the catalytic metal ion with the pro-RP oxygen, in line with experiments with the antigenomic HDV ribozyme. Experiments without divalent ions, with a double mutant that interferes with Mg(2+) binding, or with C75 deleted suggest that the pro-RP oxygen plays at most a redundant role in positioning C75. Quantum mechanical/molecular mechanical (QM/MM) studies indicate that the metal ion contributes to catalysis by interacting with both the pro-RP oxygen and the nucleophilic 2'-hydroxyl, supporting the experimental findings.

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.

The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies

Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.

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.

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.

Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.

Deciphering the role of glucosamine-6-phosphate in the riboswitch action of glmS ribozyme

The GlmS ribozyme is believed to exploit a general acid-base catalytic mechanism in the presence of glucosamine-6-phosphate (GlcN6P) to accelerate self-cleavage by approximately six orders of magnitude. The general acid and general base are not known, and the role of the GlcN6P cofactor is even less well understood. The amine group of GlcN6P has the ability to either accept or donate a proton and could therefore potentially act as an acid or a base. In order to decipher the role of GlcN6P in the self-cleavage of glmS, we have determined the preferred protonation state of the amine group in the wild-type and an inactive G40A mutant using molecular dynamics simulations and free energy calculations. Here we show that, upon binding of GlcN6P to wild-type glmS, the pK(a) of the amine moiety is altered by the active site environment, decreasing by about 2.2 from a solution pK(a) of about 8.2. On the other hand, we show that the pK(a) of the amine group slightly increases to about 8.4 upon binding to the G40A inactive mutant of glmS. These results suggest that GlcN6P acts as a general acid in the self-cleavage of glmS. Upon binding to glmS, GlcN6P can easily release a proton to the 5'-oxygen of G1 during self-cleavage of the backbone phosphodiester bond. However, in the G40A inactive mutant of glmS, the results suggest that the ability of GlcN6P to easily release its proton is diminished, in addition to the possible lack of G40 as an effective base.

Identification of catalytic metal ion ligands in ribozymes

Site-bound metal ions participate in the catalytic mechanisms of many ribozymes. Understanding these mechanisms therefore requires knowledge of the specific ligands on both substrate and ribozyme that coordinate these catalytic metal ions. A number of different structural and biochemical strategies have been developed and refined for identifying metal ion binding sites within ribozymes, and for assessing the catalytic contributions of the metal ions bound at those sites. We review these approaches and provide examples of their application, focusing in particular on metal ion rescue experiments and their roles in the construction of the transition state models for the Tetrahymena group I and RNase P ribozymes.

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.

A catalytic metal ion interacts with the cleavage Site G.U wobble in the HDV ribozyme

The HDV ribozyme self-cleaves by a chemical mechanism involving general acid-base catalysis to generate 2',3'-cyclic phosphate and 5'-hydroxyl termini. Biochemical studies from several laboratories have implicated C75 as the general acid and hydrated magnesium as the general base. We have previously shown that C75 has a pK(a) shifted >2 pH units toward neutrality [Gong, B., Chen, J. H., Chase, E., Chadalavada, D. M., Yajima, R., Golden, B. L., Bevilacqua, P. C., and Carey, P. R. (2007) J. Am. Chem. Soc. 129, 13335-13342], while in crystal structures, it is well-positioned for proton transfer. However, no evidence for a hydrated magnesium poised to serve as a general base in the reaction has been observed in high-resolution crystal structures of various reaction states and mutants. Herein, we use solution kinetic experiments and parallel Raman crystallographic studies to examine the effects of pH on the rate and Mg(2+) binding properties of wild-type and 7-deazaguanosine mutants of the HDV ribozyme. These data suggest that a previously unobserved hydrated magnesium ion interacts with N7 of the cleavage site G.U wobble base pair. Integrating this metal ion binding site with the available crystal structures provides a new three-dimensional model for the active site of the ribozyme that accommodates all available biochemical data and appears competent for catalysis. The position of this metal is consistent with a role of a magnesium-bound hydroxide as a general base as dictated by biochemical data.

Mechanistic characterization of the HDV genomic ribozyme: solvent isotope effects and proton inventories in the absence of divalent metal ions support C75 as the general acid

The hepatitis delta virus (HDV) ribozyme uses the nucleobase C75 and a hydrated Mg(2+) ion as the general acid-base catalysts in phosphodiester bond cleavage at physiological salt. A mechanistic framework has been advanced that involves one Mg(2+)-independent and two Mg(2+)-dependent channels. The rate-pH profile for wild-type (WT) ribozyme in the Mg(2+)-free channel is inverted relative to the fully Mg(2+)-dependent channel, with each having a near-neutral pKa. Inversion of the rate-pH profile was used as the crux of a mechanistic argument that C75 serves as general acid both in the presence and absence of Mg(2+). However, subsequent studies on a double mutant (DM) ribozyme suggested that the pKa observed for WT in the absence of Mg(2+) arises from ionization of C41, a structural nucleobase. To investigate this further, we acquired rate-pH/pD profiles and proton inventories for WT and DM in the absence of Mg(2+). Corrections were made for effects of ionic strength on hydrogen ion activity and pH meter readings. Results are accommodated by a model wherein the Mg(2+)-free pKa observed for WT arises from ionization of C75, and DM reactivity is compromised by protonation of C41. The Brønsted base appears to be water or hydroxide ion depending on pH. The observed pKa's are related to salt-dependent pH titrations of a model oligonucleotide, as well as electrostatic calculations, which support the local environment for C75 in the absence of Mg(2+) being similar to that in the presence of Mg(2+) and impervious to bulk ions. Accordingly, the catalytic role of C75 as the general acid does not appear to depend on divalent ions or the identity of the Brønsted base.

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid-base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK (a) values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK (a) values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid-base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid-base catalysis indicate that k (1), the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK (a). The inference that k (1) is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK (a) is that of the general acid and the higher pK (a) is that of the general base.

Catalytic strategies of self-cleaving ribozymes

[Structure: see text]. Five naturally occurring nucleolytic ribozymes have been identified: the hammerhead, hairpin, glmS, hepatitis delta virus (HDV), and Varkud satellite (VS) ribozymes. All of these RNA enzymes catalyze self-scission of the RNA backbone using a chemical mechanism equivalent to that of RNase A. RNase A uses four basic strategies to promote this reaction: geometric constraints, activation of the nucleophile, transition-state stabilization, and leaving group protonation. In this Account, we discuss the current thinking on how nucleolytic ribozymes harness RNase A's four sources of catalytic power. The geometry of the phosphodiester cleavage reaction constrains the nucleotides flanking the scissile phosphate so that they are unstacked from a canonical A-form helix and thus require alternative stabilization. Crystal structures and mutational analysis reveal that cross-strand base pairing, along with unconventional stacking and tertiary hydrogen-bonding interactions, work to stabilize the splayed conformation in nucleolytic ribozymes. Deprotonation of the 2'-OH nucleophile greatly increases its nucleophilicity in the strand scission reaction. Crystal structures of the hammerhead, hairpin, and glmS ribozymes reveal the N1 of a G residue within hydrogen-bonding distance of the 2'-OH. In each case, this residue has also been shown to be important for catalysis. In the HDV ribozyme, a hydrated magnesium has been implicated as the general base. Catalysis by the VS ribozyme requires both an A and a G, but the precise role of either has not been elucidated. Enzymes can lower the energy of a chemical reaction by binding more tightly to the transition state than to the ground states. Comparison of the hairpin ground- and transition-state mimic structures reveal greater hydrogen bonding to the transition-state mimic structure, suggesting transition-state stabilization as a possible catalytic strategy. However, the hydrogen-bonding pattern in the glmS ribozyme transition-state mimic structure and the ground-state structures are equivalent. Protonation of the 5'-O leaving group by a variety of functional groups can promote the cleavage reaction. In the HDV ribozyme, the general acid is a conserved C residue. In the hairpin ribozyme, a G residue has been implicated in protonation of the leaving group. An A in the hammerhead ribozyme probably plays a similar role. In the glmS ribozyme, an exogenous cofactor may provide the general acid. This diversity is in contrast to the relatively small number of functional groups that serve as a general base, where at least three of the nucleolytic ribozymes may use the N1 of a G.

An important role of G638 in the cis-cleavage reaction of the Neurospora VS ribozyme revealed by a novel nucleotide analog incorporation method

We describe a chemical coupling procedure that allows joining of two RNAs, one of which contains a site-specific base analog substitution, in the absence of divalent ions. This method allows incorporation of nucleotide analogs at specific positions even into large, cis-cleaving ribozymes. Using this method we have studied the effects of substitution of G638 in the cleavage site loop of the VS ribozyme with a variety of purine analogs having different functional groups and pK(a) values. Cleavage rate versus pH profiles combined with kinetic solvent isotope experiments indicate an important role for G638 in proton transfer during the rate-limiting step of the cis-cleavage reaction.

Evidence for proton transfer in the rate-limiting step of a fast-cleaving Varkud satellite ribozyme

A fast-cleaving version of the Varkud satellite ribozyme, called RG, shows an apparent cis-cleavage rate constant of 5 sec(-1), similar to the rates of protein enzymes that catalyze similar reactions. Here, we describe mutational, pH-rate, and kinetic solvent isotope experiments that investigate the identity and rate constant of the rate-limiting step in this reaction. Self-cleavage of RG exhibits a bell-shaped rate vs. pH profile with apparent pK(a)s of 5.8 and 8.3, consistent with the protonation state of two nucleotides being important for the rate of cleavage. Cleavage experiments in heavy water (D(2)O) revealed a kinetic solvent isotope effect consistent with proton transfer in the rate-limiting step. A mutant RNA that disrupts a peripheral loop-loop interaction involved in RNA folding exhibits pH- and D(2)O-independent cleavage approximately 10(3)-fold slower than wild type, suggesting that this mutant is limited by a different step than wild type. Substitution of adenosine 756 in the putative active-site loop with cytosine also decreases the cleavage rate approximately 10(3)-fold, but the A756C mutant retains pH- and D(2)O-sensitivity similar to wild type, consistent with this mutant and wild type being limited by the chemical step of the reaction. These results suggest that the RG ribozyme provides a good experimental system to investigate the nature of fast, rate-limiting steps in a ribozyme cleavage reaction.

Structural Roles of Monovalent Cations in the HDV Ribozyme

The hepatitis delta virus (HDV) ribozyme catalyzes viral RNA self-cleavage through general acid-base chemistry in which an active-site cytidine and at least one metal ion are involved. Monovalent metal ions support slow catalysis and were proposed to substitute for structural, but not catalytic, divalent metal ions in the RNA. To investigate the role of monovalent cations in ribozyme structure and function, we determined the crystal structure of the precursor HDV ribozyme in the presence of thallium ions (Tl(+)). Two Tl(+) ions can occupy a previously observed divalent metal ion hexahydrate-binding site located near the scissile phosphate, but are easily competed away by cobalt hexammine, a magnesium hexahydrate mimic and potent reaction inhibitor. Intriguingly, a third Tl(+) ion forms direct inner-sphere contacts with the ribose 2'-OH nucleophile and the pro-S(p) scissile phosphate oxygen. We discuss possible structural and catalytic implications of monovalent cation binding for the HDV ribozyme mechanism.

Long-range impact of peripheral joining elements on structure and function of the hepatitis delta virus ribozyme

The HDV ribozyme is an RNA enzyme from the human pathogenic hepatitis delta virus (HDV) that has recently also been identified in the human genome. It folds into a compact, nested double-pseudoknot. We examined here the functional relevance of the capping loop L4 and the helical crossover J1/2, which tightly interlace the two helical stacks of the ribozyme. Peripheral structural elements such as these are present in cis-acting, but not trans-acting ribozymes, which may explain the order-of-magnitude decrease in cleavage activity observed in trans-acting ribozymes with promise in gene therapy applications. Comparison of a systematic set of cis- and trans-acting HDV ribozymes shows that the absence of either L4 or J1/2 significantly and independently impacts catalytic activity. Using terbium(III) footprinting and affinity studies, as well as distance measurements based on time-resolved fluorescence resonance energy transfer, we find that J1/2 is most important for conferring structural properties similar to those of the cis-acting ribozyme. Our results are consistent with a model in which removal of either a helical crossover or surprisingly a capping loop induces greater dynamics and expansion of the catalytic core at long range, impacting local and global folding, as well as catalytic function.

Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.

Chemical rescue, multiple ionizable groups, and general acid-base catalysis in the HDV genomic ribozyme

In the ribozyme from the hepatitis delta virus (HDV) genomic strand RNA, a cytosine side chain is proposed to facilitate proton transfer in the transition state of the reaction and, thus, act as a general acid-base catalyst. Mutation of this active-site cytosine (C75) reduced RNA cleavage rates by as much as one million-fold, but addition of exogenous cytosine and certain nucleobase or imidazole analogs can partially rescue activity in these mutants. However, pH-rate profiles for the rescued reactions were bell shaped, and only one leg of the pH-rate curve could be attributed to ionization of the exogenous nucleobase or buffer. When a second potential ionizable nucleobase (C41) was removed, one leg of the bell-shaped curve was eliminated in the chemical-rescue reaction. With this construct, the apparent pK(a) determined from the pH-rate profile correlated with the solution pK(a) of the buffer, and the contribution of the buffer to the rate enhancement could be directly evaluated in a free-energy or Brønsted plot. The free-energy relationship between the acid dissociation constant of the buffer and the rate constant for cleavage (Brønsted value, beta, = approximately 0.5) was consistent with a mechanism in which the buffer acted as a general acid-base catalyst. These data support the hypothesis that cytosine 75, in the intact ribozyme, acts as a general acid-base catalyst.

Ligand requirements for glmS ribozyme self-cleavage

Natural RNA catalysts (ribozymes) perform essential reactions in biological RNA processing and protein synthesis, whereby catalysis is intrinsic to RNA structure alone or in combination with metal ion cofactors. The recently discovered glmS ribozyme is unique in that it functions as a glucosamine-6-phosphate (GlcN6P)-dependent catalyst believed to enable "riboswitch" regulation of amino-sugar biosynthesis in certain prokaryotes. However, it is unclear whether GlcN6P functions as an effector or coenzyme to promote ribozyme self-cleavage. Herein, we demonstrate that ligand is absolutely requisite for glmS ribozyme self-cleavage activity. Furthermore, catalysis both requires and is dependent upon the acid dissociation constant (pKa) of the amine functionality of GlcN6P and related compounds. The data demonstrate that ligand is integral to catalysis, consistent with a coenzyme role for GlcN6P and illustrating an expanded capacity for biological RNA catalysis.

Chemical biology at the crossroads of molecular structure and mechanism

Chemical insight into biological function is the holy grail of structural biology. Small molecules are central players as building blocks, effectors and probes of macromolecular structure and function.

Non-identical electronic characters of the internucleotidic phosphates in RNA modulate the chemical reactivity of the phosphodiester bonds

We here show that the electronic properties and the chemical reactivities of the internucleotidic phosphates in the heptameric ssRNAs are dissimilar in a sequence-specific manner because of their non-identical microenvironments, in contrast with the corresponding isosequential ssDNAs. This has been evidenced by monitoring the delta H8(G) shifts upon pH-dependent ionization (pK(a1)) of the central 9-guaninyl (G) to the 9-guanylate ion (G-), and its electrostatic effect on each of the internucleotidic phosphate anions, as measured from the resultant delta 31P shifts (pKa2) in the isosequential heptameric ssRNAs vis-à-vis ssDNAs: [d/r(5'-Cp1Ap2Q1p3Gp4Q2p5Ap6C-3'): Q1 = Q2 = A (5a/5b) or C (8a/8b), Q(1) = A, Q(2) = C (6a/6b), Q1 = C, Q2 = A (7a/7b)]. These oligos with single ionizable G in the centre are chosen because of the fact that the pseudoaromatic character of G can be easily modulated in a pH-dependent manner by its transformation to G- (the 2'-OH to 2-O- ionization effect is not detectable below pH 11.6 as evident from the N(1-Me)-G analog), thereby modulating/titrating the nature of the electrostatic interactions of G to G- with the phosphates, which therefore constitute simple models to interrogate how the variable pseudoaromatic characters of nucleobases under different sequence context (J. Am. Chem. Soc., 2004, 126, 8674-8681) can actually influence the reactivity of the internucleotide phosphates as a result of modulation of sequence context-specific electrostatic interactions. In order to better understand the impact of the electrostatic effect of the G to G- on the tunability of the electronic character of internucleotidic phosphates in the heptameric ssRNAs 5b, 6b, 7b and 8b, we have also performed their alkaline hydrolysis at pH 12.5 at 20 degrees C, and have identified the preferences of the cleavage sites at various phosphates, which are p2, p3 and p4 (Fig.3). The results of these alkaline hydrolysis studies have been compared with the hydrolysis of analogous N(1-Me)-G heptameric ssRNA sequences 5c, 7c and 8c under identical conditions in order to establish the role of the electrostatic effect of the 9-guanylate ion (and the 2'-OH to 2-O- ionization) on the internucleotidic phosphate. It turned out that the relative alkaline hydrolysis rate at those particular phosphates (p2, p3 and p4) in the N(1-Me)-G heptamers was reduced from 16-78% compared to those in the native counterparts [Fig. 4, and ESI 2 (Fig. S11)]. Thus, these physico-chemical studies have shown that those p2, p3 and p4 phosphates in the native heptameric RNAs, which show pKa2 as well as more deshielding (owing to weaker 31P screening) in the alkaline pH compared to those at the neutral pH, are more prone to the alkaline hydrolysis because of their relatively enhanced electrophilic character resulting from weaker 31P screening. This screening effect originates as a result of the systematic charge repulsion effect between the electron cloud in the outermost orbitals of phosphorus and the central guanylate ion, leading to delocalization of the phosphorus p(pi) charge into its dpi orbitals. It is thus likely that, just as in the non-enzymatic hydrolysis, the enzymatic hydrolysis of a specific phosphate in RNA by general base-catalysis in RNA-cleaving proteins (RNase A, RNA phosphodiesterase or nuclease) can potentially be electrostatically influenced by tuning the transient charge on the nucleobase in the steric proximity or as a result of specific sequence context owing to nearest-neighbor interactions.

The RNA-Cleaving Bipartite DNAzyme Is a Distinctive Metalloenzyme

Much interest has focused on the mechanisms of the five naturally occurring self-cleaving ribozymes, which, in spite of catalyzing the same reaction, adopt divergent strategies. These ribozymes, with the exception of the recently described glmS ribozyme, do not absolutely require divalent metal ions for their catalytic chemistries in vitro. A mechanistic investigation of an in vitro-selected, RNA-cleaving DNA enzyme, the bipartite, which catalyzes the same chemistry as the five natural self-cleaving ribozymes, found a mechanism of significant complexity. The DNAzyme showed a bell-shaped pH profile. A dissection of metal usage indicated the involvement of two catalytically relevant magnesium ions for optimal activity. The DNAzyme was able to utilize manganese(II) as well as magnesium; however, with manganese it appeared to function complexed to either one or two of those cations. Titration with hexaamminecobalt(III) chloride inhibited the activity of the bipartite; this suggests that it is a metalloenzyme that utilizes metal hydroxide as a general base for activation of its nucleophile. Overall, the bipartite DNAzyme appeared to be kinetically distinct not only from the self-cleaving ribozymes but also from other in vitro-selected, RNA-cleaving deoxyribozymes, such as the 8-17, 10-23, and 614.

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.

The catalytic diversity of RNAs

The natural RNA enzymes catalyse phosphate-group transfer and peptide-bond formation. Initially, metal ions were proposed to supply the chemical versatility that nucleotides lack. In the ensuing decades, structural and mechanistic studies have substantially altered this initial viewpoint. Whereas self-splicing ribozymes clearly rely on essential metal-ion cofactors, self-cleaving ribozymes seem to use nucleotide bases for their catalytic chemistry. Despite the overall differences in chemical features, both RNA and protein enzymes use similar catalytic strategies.

Ribozyme catalysis: not different, just worse

Evolution has resoundingly favored protein enzymes over RNA-based catalysts, yet ribozymes occupy important niches in modern cell biology that include the starring role in catalysis of protein synthesis on the ribosome. Recent results from structural and biochemical studies show that natural ribozymes use an impressive range of catalytic mechanisms, beyond metalloenzyme chemistry and analogous to more chemically diverse protein enzymes. These findings make it increasingly possible to compare details of RNA- and protein-based catalysis.

Linkage between proton binding and folding in RNA: implications for RNA catalysis

Small ribozymes use their nucleobases to catalyse phosphodiester bond cleavage. The hepatitis delta virus ribozyme employs C75 as a general acid to protonate the 5′-bridging oxygen leaving group, and to accomplish this task efficiently, it shifts its pKa towards neutrality. Simulations and thermodynamic experiments implicate linkage between folding and protonation in nucleobase pKa shifting. Even small oligonucleotides are shown to fold in a highly co-operative manner, although they do so in a context-specific fashion. Linkage between protonation and co-operativity of folding may drive pKa shifting and provide for enhanced function in RNA.

General acid catalysis by the hepatitis delta virus ribozyme

Recent crystallographic and functional analyses of RNA enzymes have raised the possibility that the purine and pyrimidine nucleobases may function as general acid-base catalysts. However, this mode of nucleobase-mediated catalysis has been difficult to establish unambiguously. Here, we used a hyperactivated RNA substrate bearing a 5'-phosphorothiolate to investigate the role of a critical cytosine residue in the hepatitis delta virus ribozyme. The hyperactivated substrate specifically suppressed the deleterious effects of cytosine mutations and pH changes, thereby linking the protonation of the nucleobase to leaving-group stabilization. We conclude that the active-site cytosine provides general acid catalysis, mediating proton transfer to the leaving group through a protonated N3-imino nitrogen. These results establish a specific role for a nucleobase in a ribozyme reaction and support the proposal that RNA nucleobases may function in a manner analogous to that of catalytic histidine residues in protein enzymes.

Role of an active site adenine in hairpin ribozyme catalysis

The hairpin ribozyme is a small catalytic RNA that accelerates reversible cleavage of a phosphodiester bond. Structural and mechanistic studies suggest that divalent metals stabilize the functional structure but do not participate directly in catalysis. Instead, two active site nucleobases, G8 and A38, appear to participate in catalytic chemistry. The features of A38 that are important for active site structure and chemistry were investigated by comparing cleavage and ligation reactions of ribozyme variants with A38 modifications. An abasic substitution of A38 reduced cleavage and ligation activity by 14,000-fold and 370,000-fold, respectively, highlighting the critical role of this nucleobase in ribozyme function. Cleavage and ligation activity of unmodified ribozymes increased with increasing pH, evidence that deprotonation of some functional group with an apparent pK(a) value near 6 is important for activity. The pH-dependent transition in activity shifted by several pH units in the basic direction when A38 was substituted with an abasic residue, or with nucleobase analogs with very high or low pK(a) values that are expected to retain the same protonation state throughout the experimental pH range. Certain exogenous nucleobases that share the amidine group of adenine restored activity to abasic ribozyme variants that lack A38. The pH dependence of chemical rescue reactions also changed according to the intrinsic basicity of the rescuing nucleobase, providing further evidence that the protonation state of the N1 position of purine analogs is important for rescue activity. These results are consistent with models of the hairpin ribozyme catalytic mechanism in which interactions with A38 provide electrostatic stabilization to the transition state.