Two-metal-ion mechanism for hammerhead-ribozyme catalysis

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

Theoretical Studies of the Self Cleavage Pistol Ribozyme Mechanism

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

Kinetic Study of the Avocado Sunblotch Viroid Self-Cleavage Reaction Reveals Compensatory Effects between High-Pressure and High-Temperature: Implications for Origins of Life on Earth

A high pressure apparatus allowing one to study enzyme kinetics under pressure was used to study the self-cleavage activity of the avocado sunblotch viroid. The kinetics of this reaction were determined under pressure over a range up to 300 MPa (1-3000 bar). It appears that the initial rate of this reaction decreases when pressure increases, revealing a positive ΔV≠ of activation, which correlates with the domain closure accompanying the reaction and the decrease of the surface of the viroid exposed to the solvent. Although, as expected, temperature increases the rate of the reaction whose energy of activation was determined, it appeared that it does not significantly influence the ΔV≠ of activation and that pressure does not influence the energy of activation. These results provide information about the structural aspects or this self-cleavage reaction, which is involved in the process of maturation of this viroid. The behavior of ASBVd results from the involvement of the hammerhead ribozyme present at its catalytic domain, indeed a structural motif is very widespread in the ancient and current RNA world.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Divalent Metal Ion Activation of a Guanine General Base in the Hammerhead Ribozyme: Insights from Molecular Simulations

The hammerhead ribozyme is a well-studied nucleolytic ribozyme that catalyzes the self-cleavage of the RNA phosphodiester backbone. Despite experimental and theoretical efforts, key questions remain about details of the mechanism with regard to the activation of the nucleophile by the putative general base guanine (G12). Straightforward interpretation of the measured activity-pH data implies the pK_a value of the N1 position in the G12 nucleobase is significantly shifted by the ribozyme environment. Recent crystallographic and biochemical work has identified pH-dependent divalent metal ion binding at the N7/O6 position of G12, leading to the hypothesis that this binding mode could induce a pK_a shift of G12 toward neutrality. We present computational results that support this hypothesis and provide a model that unifies the interpretation of available structural and biochemical data and paints a detailed mechanistic picture of the general base step of the reaction. Experimentally testable predictions are made for mutational and rescue effects on G12, which will give further insights into the catalytic mechanism. These results contribute to our growing knowledge of the potential roles of divalent metal ions in RNA catalysis.

Catalytic metal ions and enzymatic processing of DNA and RNA

CONSPECTUS: Two-metal-ion-dependent nucleases cleave the phosphodiester bonds of nucleic acids via the two-metal-ion (2M) mechanism. Several high-resolution X-ray structures portraying the two-metal-aided catalytic site, together with mutagenesis and kinetics studies, have demonstrated a functional role of the ions for catalysis in numerous metallonucleases. Overall, the experimental data confirm the general mechanistic hypothesis for 2M-aided phosphoryl transfer originally reported by Steitz and Steitz ( Proc. Natl. Acad. Sci. U.S.A. 1993 , 90 ( 14 ), 6498 - 6502 ). This seminal paper proposed that one metal ion favors the formation of the nucleophile, while the nearby second metal ion facilitates leaving group departure during RNA hydrolysis. Both metals were suggested to stabilize the enzymatic transition state. Nevertheless, static X-ray structures alone cannot exhaustively unravel how the two ions execute their functional role along the enzymatic reaction during processing of DNA or RNA strands when moving from reactants to products, passing through metastable intermediates and high-energy transition states. In this Account, we discuss the role of multiscale molecular simulations in further disclosing mechanistic insights of 2M-aided catalysis for two prototypical enzymatic targets for drug discovery, namely, ribonuclease H (RNase H) and type II topoisomerase (topoII). In both examples, first-principles molecular simulations, integrated with structural data, emphasize a cooperative motion of the bimetal motif during catalysis. The coordinated motion of both ions is crucial for maintaining a flexible metal-centered structural architecture exquisitely tailored to accommodate the DNA or RNA sugar-phosphate backbone during phosphodiester bond cleavage. Furthermore, our analysis of RNase H and the N-terminal domain (PAN) of influenza polymerase shows that classical molecular dynamics simulations coupled with enhanced sampling techniques have contributed to describe the modulatory effect of metal ion concentration and metal uptake on the 2M mechanism and efficiency. These aspects all point to the emerging and intriguing role of additional adjacent ions potentially involved in the modulation of phosphoryl transfer reactions and enzymatic turnover in 2M-catalysis, as recently observed experimentally in polymerase η and homing endonuclease I-DmoI. These computational results, integrated with experimental findings, describe and reinforce the nascent concept of a functional and cooperative dynamics of the catalytic metal ions during the 2M-dependent enzymatic processing of DNA and RNA. Encouraged by the insights provided by computational approaches, we foresee further experiments that will feature the functional and joint dynamics of the catalytic metal ions for nucleic acid processing. This could impact the de novo design of artificial metallonucleases and the rational design of potent metal-chelating inhibitors of pharmaceutically relevant enzymes.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

Modeling the RNA 2'OH activation: possible roles of metal ion and nucleobase as catalysts in self-cleaving ribozymes

The RNA 2'OH activation as taking place in the first chemical step of self-cleaving ribozymes is studied theoretically by DFT and MP2 methods using a continuum solvation model (CPCM). The reaction of proton transfer is studied in the presence of two kinds of catalysts: a fully hydrated metal ion (Mg(2+)) or partially hydrated nucleobase (guanine), taken separately or together leading to three different modes of activation. The metal ion is either directly bound (inner-sphere) or indirectly bound (outer-sphere) to the 2'OH group and a hydroxide ion acts as a general or specific base; the nucleobase is taken in anionic or in neutral enol-tautomeric forms playing itself the role of general base. The presence of a close metal ion (outer-sphere) lowers the pK(a) value of the 2'OH group by several log units in both metal-ion and nuleobase catalysis. The direct metal coordination to the 2'OH group (inner-sphere) further stabilizes the developing negative charge on the nucleophile. The switching from the inner-sphere to the outer-sphere coordination appears to be driven by the energy cost for reorganizing the first coordination shell rather than by the electrostatic repulsion between the ligands. The metal-ion catalysis is more effective with a specific base in the dianionic mechanism. On the other hand, the nucleobase catalysis is more effective in the monoanionic mechanism and in the presence of a metal ion acting as a cofactor through nonspecific electrostatic interactions. The results establish a baseline to study the possible roles of metal and nucleobase catalysts and their environment in more realistic models for self-cleaving ribozymes.

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

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

Editing mechanism of aminoacyl-tRNA synthetases operates by a hybrid ribozyme/protein catalyst

Aminoacyl-tRNA synthetases (aaRSs) are critical for the translational process, catalyzing the attachment of specific amino acids to their cognate tRNAs. To ensure formation of the correct aminoacyl-tRNA, and thereby enhance the reliability of translation, several aaRSs have an editing capability that hinders formation of misaminoacylated tRNAs. We investigated theoretically the mechanism of the editing reaction for a class I enzyme, leucyl-tRNA synthetase (LeuRS), complexed with a misaminoacylated tRNA(Leu), employing ab initio hybrid quantum mechanical/molecular mechanical potentials in conjunction with molecular dynamics simulations. It is shown that the water molecule that acts as the nucleophile in the editing reaction is activated by a 3'-hydroxyl group at the 3'-end of tRNA(Leu) and that the O2' atom of the leaving group of the substrate is capped by one of the water's hydrogen atoms. Thus, it is shown that editing is a self-cleavage reaction of the tRNA and so it is the tRNA, and not the protein, that drives the reaction. The protein does, however, have an important stabilizing effect on some high-energy intermediates along the reaction path, which is more efficient than the ribozyme would be alone. This indicates that editing is achieved by a novel "hybrid ribozyme/protein catalyst". Analysis of existing experimental data and additional modeling shows that this ribozymal mechanism appears to be widespread, occurring in the ribosome as well as in other aaRSs. It also suggests transitional forms that could have played an important role in the RNA world hypothesis for the origin of life.

On the role of Hoogsteen:Hoogsteen interactions in RNA: ab initio investigations of structures and energies

We use a combination of database analysis and quantum chemical studies to investigate the role of cis and trans Hoogsteen:Hoogsteen (H:H) base pairs and associated higher-order structures in RNA. We add three new examples to the list of previously identified base-pair combinations belonging to these families and, in addition to contextual classification and characterization of their structural and energetic features, we compare their interbase interaction energies and propensities toward participation in triplets and quartets. We find that some base pairs, which are nonplanar in their isolated minimum energy geometries, attain planarity and stability upon triplet formation. A:A H:H trans is the most frequent H:H combination in RNA structures. This base pair occurs at many distinct positions in known rRNA structures, where it helps in the interaction of ribosomal domains in the 50S subunit. It is also present as a part of tertiary interaction in tRNA structures. Although quantum chemical studies suggest an intrinsically nonplanar geometry for this base pair in isolated form, it has the tendency to attain planar geometry in RNA crystal structures by forming higher-order tertiary interactions or in the presence of additional base-phosphate interactions. The tendency of this base pair to form such additional interactions may be helpful in bringing together different segments of RNA, thus making it suitable for the role of facilitator for RNA folding. This also explains the high occurrence frequency of this base pair among all H:H interactions.

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

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

Probing general acid catalysis in the hammerhead ribozyme

Recent crystallographic and computational studies have provided fresh insights into the catalytic mechanism of the RNA-cleaving hammerhead ribozyme. Based on these findings, specific ribozyme functional groups have been hypothesized to act directly as the general acid and base catalysts, although the catalytic role of divalent metal cations (M(2+)) remains uncertain. We now report a functional characterization of the general acid catalysis mechanism and the role of an M(2+) cofactor therein, for the S. mansoni hammerhead (an "extended" hammerhead ribozyme). We have compared hammerhead cleavage of substrates with natural (ribo-phosphodiester) versus bridging-5'-phosphorothioate scissile linkages, in the contexts of active site mutations and M(2+) substitution. Cleavage of the natural substrate is inhibited by modification of the G8 2'-OH ribozyme residue and depends strongly upon the presence and identity of an M(2+) cofactor; in contrast, cleavage of the bridging-phosphorothioate substrate is conspicuously insensitive to any of these factors. These results imply that (1) both an M(2+) cofactor and the G8 2'-OH play crucial roles in hammerhead general acid catalysis and (2) the M(2+) cofactor does not contribute to general acid catalysis via Lewis acid stabilization of the leaving group. General acid pK(a) perturbation was also demonstrated for both M(2+) substitution and G8 2'-OH modification, which suggests transition state M(2+) coordination of the G8 2'-OH, to lower its pK(a) and improve its ability to transfer a proton to the leaving group. We also report a simple method for synthesizing radiolabeled bridging-5'-phosphorothioate substrates.

Insights into functional modulation of catalytic RNA activity

RNA molecules play critical roles in cell biology, and novel findings continuously broaden their functional repertoires. Apart from their well-documented participation in protein synthesis, it is now apparent that several noncoding RNAs (i.e., micro-RNAs and riboswitches) also participate in the regulation of gene expression. The discovery of catalytic RNAs had profound implications on our views concerning the evolution of life on our planet at a molecular level. A characteristic attribute of RNA, probably traced back to its ancestral origin, is the ability to interact with and be modulated by several ions and molecules of different sizes. The inhibition of ribosome activity by antibiotics has been extensively used as a therapeutical approach, while activation and substrate-specificity alteration have the potential to enhance the versatility of ribozyme-based tools in translational research. In this review, we will describe some representative examples of such modulators to illustrate the potential of catalytic RNAs as tools and targets in research and clinical approaches.

The linkage between ribosomal crystallography, metal ions, heteropolytungstates and functional flexibility

Crystallography of ribosomes, the universal cell nucleoprotein assemblies facilitating the translation of the genetic-code into proteins, met with severe problems owing to their large size, complex structure, inherent flexibility and high conformational variability. For the case of the small ribosomal subunit, which caused extreme difficulties, post crystallization treatment by minute amounts of a heteropolytungstate cluster allowed structure determination at atomic resolution. This cluster played a dual role in ribosomal crystallography: providing anomalous phasing power and dramatically increased the resolution, by stabilization of a selected functional conformation. Thus, four out of the fourteen clusters that bind to each of the crystallized small subunits are attached to a specific ribosomal protein in a fashion that may control a significant component of the subunit internal flexibility, by "gluing" symmetrical related subunits. Here we highlight basic issues in the relationship between metal ions and macromolecules and present common traits controlling in the interactions between polymetalates and various macromolecules, which may be extended towards the exploitation of polymetalates for therapeutical treatment.

General base catalysis for cleavage by the active-site cytosine of the hepatitis delta virus ribozyme: QM/MM calculations establish chemical feasibility

The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.

Solvent structure and hammerhead ribozyme catalysis

Although the hammerhead ribozyme is regarded as a prototype for understanding RNA catalysis, the mechanistic roles of associated metal ions and water molecules in the cleavage reaction remain controversial. We have investigated the catalytic potential of observed divalent metal ions and water molecules bound to a 2 A structure of the full-length hammerhead ribozyme by using X-ray crystallography in combination with molecular dynamics simulations. A single Mn(2+) is observed to bind directly to the A9 phosphate in the active site, accompanying a hydrogen-bond network involving a well-ordered water molecule spanning N1 of G12 (the general base) and 2'-O of G8 (previously implicated in general acid catalysis) that we propose, based on molecular dynamics calculations, facilitates proton transfer in the cleavage reaction. Phosphate-bridging metal interactions and other mechanistic hypotheses are also tested with this approach.

The associative nature of adenylyl transfer catalyzed by T4 DNA ligase

DNA ligase seals nicks in dsDNA using chemical energy of the phosphoanhydride bond in ATP or NAD(+) and assistance of a divalent metal cofactor Mg(2+). Molecular details of ligase catalysis are essential for understanding the mechanism of metal-promoted phosphoryl transfer reactions in the living cell responsible for a wide range of processes, e.g., DNA replication and transcription, signaling and differentiation, energy coupling and metabolism. Here we report a single-turnover (31)P solid-state NMR study of adenylyl transfer catalyzed by DNA ligase from bacteriophage T4. Formation of a high-energy covalent ligase-nucleotide complex is triggered in situ by the photo release of caged Mg(2+), and sequentially formed intermediates are monitored by NMR. Analyses of reaction kinetics and chemical-shift changes indicate that the pentacoordinated phosphorane intermediate builds up to 35% of the total reacting species after 4-5 h of reaction. This is direct experimental evidence of the associative nature of adenylyl transfer catalyzed by DNA ligase. NMR spectroscopy in rotating solids is introduced as an analytical tool for recording molecular movies of reaction processes. Presented work pioneers a promising direction in structural studies of biochemical transformations.

Roles of metal ions in nucleases

The hydrolysis of phosphodiester bonds by metallonucleases is crucial to most aspects of nucleic acid processing. In recent years, studies of the classical restriction endonucleases have given way to the characterization of metallonucleases with widely divergent active site motifs. These developments fuel debates regarding the roles of metal ions in these enzymes. It is fortuitous that the current literature also includes the increased application of a variety of computational techniques to test the roles of metal ions in nucleic acid hydrolysis by these systems. This includes recent proposals and indirect evidence that these enzymes utilize metal ion movement in these reactions.

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.

Ribozyme catalysis revisited: is water involved?

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

Hydration properties of magnesium and calcium ions from constrained first principles molecular dynamics

We studied the solvation structures of the divalent metal cations Mg(2+) and Ca(2+) in ambient water by applying a Car-Parrinello-based constrained molecular dynamics method. By employing the metal-water oxygen coordination number as a reaction coordinate, we could identify distinct aqua complexes characterized by structural variations of the first coordination shell. In particular, our estimated free-energy profile clearly shows that the global minimum for Mg(2+) is represented by a rather stable sixfold coordination in the octahedral arrangement, in agreement with experiments. Conversely, for Ca(2+) the free-energy curve shows several shallow local minima, suggesting that the hydration structure of Ca(2+) is highly variable. Implications for water exchange reactions are also discussed.

Insight into the role of Mg in hammerhead ribozyme catalysis from X-ray crystallography and molecular dynamics simulation

Results of a series of 12 ns molecular dynamics (MD) simulations of the reactant state (with and without a Mg(2+) ion), early and late transition state mimics are presented based on a recently reported crystal structure of a full-length hammerhead RNA. The simulation results support a catalytically active conformation with a Mg(2+) ion bridging the A9 and scissile phosphates. In the reactant state, the Mg(2+) spends significant time closely associated with the 2'OH of G8, but remains fairly distant from the leaving group O(5') position. In the early TS mimic simulation, where the nucleophilic O(2') and leaving group O(5') are equidistant from the phosphorus, the Mg(2+) ion remains tightly coordinated to the 2'OH of G8, but is positioned closer to the O(5') leaving group, stabilizing the accumulating charge. In the late TS mimic simulation, the coordination around the bridging Mg(2+) ion undergoes a transition whereby the coordination with the 2'OH of G8 is replace by the leaving group O(5') that has developed significant charge. At the same time, the 2'OH of G8 forms a hydrogen bond with the leaving group O(5') and is positioned to act as a general acid catalyst. This work represents the first reported simulations of the full-length hammerhead structure and TS mimics, and provides direct evidence for the possible role of a bridging Mg(2+) ion in catalysis that is consistent with both crystallographic and biochemical data.

Density Functional Theory Investigation on the Mechanism of the Hepatitis Delta Virus Ribozyme

Density functional theory methods have been used to investigate the hepatitis delta virus (HDV) ribozyme and its catalyzed phosphodiester cleavage. In particular, the effects of the environment's polarity and/or specific hydrogen-bond interactions on the proton affinity of the active site cytosine's N3 ring center have been considered. In addition, the basicities of possible hydrated Mg2+ ion species were also examined. The mechanism previously proposed for the HDV ribozyme in which the active site cytosine (C75) is protonated and thus acts as an acid while the Mg2+ species acts as the complementary base was then investigated. The possible role of tautomerization of C75 is also discussed.