Catalytic mechanism of RNA backbone cleavage by ribonuclease H from quantum mechanics/molecular mechanics simulations.

We use quantum mechanics/molecular mechanics simulations to study the cleavage of the ribonucleic acid (RNA) backbone catalyzed by ribonuclease H. This protein is a prototypical member of a large family of enzymes that use two-metal catalysis to process nucleic acids. By combining Hamiltonian replica exchange with a finite-temperature string method, we calculate the free energy surface underlying the RNA-cleavage reaction and characterize its mechanism. We find that the reaction proceeds in two steps. In a first step, catalyzed primarily by magnesium ion A and its ligands, a water molecule attacks the scissile phosphate. Consistent with thiol-substitution experiments, a water proton is transferred to the downstream phosphate group. The transient phosphorane formed as a result of this nucleophilic attack decays by breaking the bond between the phosphate and the ribose oxygen. In the resulting intermediate, the dissociated but unprotonated leaving group forms an alkoxide coordinated to magnesium ion B. In a second step, the reaction is completed by protonation of the leaving group, with a neutral Asp132 as a likely proton donor. The overall reaction barrier of ∼15 kcal mol(-1), encountered in the first step, together with the cost of protonating Asp132, is consistent with the slow measured rate of ∼1-100/min. The two-step mechanism is also consistent with the bell-shaped pH dependence of the reaction rate. The nonmonotonic relative motion of the magnesium ions along the reaction pathway agrees with X-ray crystal structures. Proton-transfer reactions and changes in the metal ion coordination emerge as central factors in the RNA-cleavage reaction.

[1]  K. Morikawa,et al.  Proposal for new catalytic roles for two invariant residues in Escherichia coli ribonuclease HI. , 1996, Protein engineering.

[2]  W. E,et al.  Finite temperature string method for the study of rare events. , 2002, Journal of Physical Chemistry B.

[3]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

[4]  E. Goldsmith,et al.  Docking interactions induce exposure of activation loop in the MAP kinase ERK2. , 2006, Structure.

[5]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[6]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[7]  P. V. von Hippel,et al.  A two-metal ion mechanism operates in the hammerhead ribozyme-mediated cleavage of an RNA substrate. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[8]  J. Champoux,et al.  RNase H activity: structure, specificity, and function in reverse transcription. , 2008, Virus research.

[9]  Bernard R. Brooks,et al.  Interfacing Q‐Chem and CHARMM to perform QM/MM reaction path calculations , 2007, J. Comput. Chem..

[10]  Stanley T Crooke,et al.  Human RNase H1 Uses One Tryptophan and Two Lysines to Position the Enzyme at the 3′-DNA/5′-RNA Terminus of the Heteroduplex Substrate* , 2003, Journal of Biological Chemistry.

[11]  S. Kanaya,et al.  Identification of the amino acid residues involved in an active site of Escherichia coli ribonuclease H by site-directed mutagenesis. , 1990, The Journal of biological chemistry.

[12]  T. Darden,et al.  Reaction mechanism of the epsilon subunit of E. coli DNA polymerase III: insights into active site metal coordination and catalytically significant residues. , 2009, Journal of the American Chemical Society.

[13]  Suse Broyde,et al.  A water-mediated and substrate-assisted catalytic mechanism for Sulfolobus solfataricus DNA polymerase IV. , 2007, Journal of the American Chemical Society.

[14]  Matteo Dal Peraro,et al.  Phosphodiester cleavage in ribonuclease H occurs via an associative two-metal-aided catalytic mechanism. , 2008, Journal of the American Chemical Society.

[15]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[16]  S. Halford,et al.  DNA cleavage by the EcoRV restriction endonuclease: pH dependence and proton transfers in catalysis. , 1999, Journal of molecular biology.

[17]  M. Klein,et al.  Modeling the charge distribution at metal sites in proteins for molecular dynamics simulations. , 2007, Journal of structural biology.

[18]  S. Crooke,et al.  Investigating the Structure of Human RNase H1 by Site-directed Mutagenesis* , 2001, The Journal of Biological Chemistry.

[19]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[20]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[21]  Y. Uchiyama,et al.  Studies of the interactions between Escherichia coli ribonuclease HI and its substrate. , 1994, Journal of molecular biology.

[22]  M. Oobatake,et al.  Dispensability of glutamic acid 48 and aspartic acid 134 for Mn2+-dependent activity of Escherichia coli ribonuclease HI. , 2003, Biochemistry.

[23]  E. Ohtsuka,et al.  Recognition of 2′‐hydroxyl groups by Escherichia coli ribonuclease HI , 1995, FEBS letters.

[24]  Bernard R. Brooks,et al.  Artificial reaction coordinate “tunneling” in free‐energy calculations: The catalytic reaction of RNase H , 2009, J. Comput. Chem..

[25]  Wei Yang,et al.  Structure of human RNase H1 complexed with an RNA/DNA hybrid: insight into HIV reverse transcription. , 2007, Molecular cell.

[26]  R. Ghirlando,et al.  Specific recognition of RNA/DNA hybrid and enhancement of human RNase H1 activity by HBD , 2008, The EMBO journal.

[27]  Wei Yang,et al.  An equivalent metal ion in one- and two-metal-ion catalysis , 2008, Nature Structural &Molecular Biology.

[28]  Wei Yang,et al.  Crystal Structures of RNase H Bound to an RNA/DNA Hybrid: Substrate Specificity and Metal-Dependent Catalysis , 2005, Cell.

[29]  Carlo Adamo,et al.  Predicting proton transfer barriers with density functional methods , 1999 .

[30]  Wei Yang,et al.  Stepwise analyses of metal ions in RNase H catalysis from substrate destabilization to product release , 2006, The EMBO journal.

[31]  C. M. Dupureur.  Roles of metal ions in nucleases. , 2008, Current opinion in chemical biology.

[32]  Christoph Dellago,et al.  Potential energy landscape for proton transfer in (H2O)3H+: comparison of density functional theory and wavefunction-based methods , 2000 .

[33]  Brigitta Elsässer,et al.  Atomistic details of the associative phosphodiester cleavage in human ribonuclease H. , 2010, Physical chemistry chemical physics : PCCP.

[34]  J. Grasby,et al.  The pH-dependence of the Escherichia coli RNase HII-catalysed reaction suggests that an active site carboxylate group participates directly in catalysis. , 2007, Journal of molecular biology.

[35]  V. Barone,et al.  Proton transfer in model hydrogen-bonded systems by a density functional approach , 1994 .

[36]  E. Kanaya,et al.  Kinetic analysis of Escherichia coli ribonuclease HI using oligomeric DNA/RNA substrates suggests an alternative mechanism for the interaction between the enzyme and the substrate. , 1995, European journal of biochemistry.

[37]  W. Stec,et al.  Stereochemical Course of Escherichia coli RNase H , 2002, Chembiochem : a European journal of chemical biology.

[38]  A. Jeltsch,et al.  On the divalent metal ion dependence of DNA cleavage by restriction endonucleases of the EcoRI family. , 2009, Journal of molecular biology.

[39]  R. Radhakrishnan,et al.  Computational delineation of the catalytic step of a high‐fidelity DNA polymerase , 2010, Protein science : a publication of the Protein Society.

[40]  Gerhard Hummer,et al.  Pore opening and closing of a pentameric ligand-gated ion channel , 2010, Proceedings of the National Academy of Sciences.

[41]  Carmay Lim,et al.  A combined experimental and theoretical study of divalent metal ion selectivity and function in proteins: application to E. coli ribonuclease H1. , 2003, Journal of the American Chemical Society.

[42]  H. Shinagawa,et al.  Two basic residues, Lys‐107 and Lys‐118, of RuvC resolvase are involved in critical contacts with the Holliday junction for its resolution , 2000, Genes to cells : devoted to molecular & cellular mechanisms.

[43]  R. Radhakrishnan.  Coupling of fast and slow modes in the reaction pathway of the minimal hammerhead ribozyme cleavage. , 2007, Biophysical journal.

[44]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .