Impact of mutation on proton transfer reactions in ketosteroid isomerase: insights from molecular dynamics simulations.

The two proton transfer reactions catalyzed by ketosteroid isomerase (KSI) involve a dienolate intermediate stabilized by hydrogen bonds with Tyr14 and Asp99. Molecular dynamics simulations based on an empirical valence bond model are used to examine the impact of mutating these residues on the hydrogen-bonding patterns, conformational changes, and van der Waals and electrostatic interactions during the proton transfer reactions. While the rate constants for the two proton transfer steps are similar for wild-type (WT) KSI, the simulations suggest that the rate constant for the first proton transfer step is smaller in the mutants due to the significantly higher free energy of the dienolate intermediate relative to the reactant. The calculated rate constants for the mutants D99L, Y14F, and Y14F/D99L relative to WT KSI are qualitatively consistent with the kinetic experiments indicating a significant reduction in the catalytic rates along the series of mutants. In the simulations, WT KSI retained two hydrogen-bonding interactions between the substrate and the active site, while the mutants typically retained only one hydrogen-bonding interaction. A new hydrogen-bonding interaction between the substrate and Tyr55 was observed in the double mutant, leading to the prediction that mutation of Tyr55 will have a greater impact on the proton transfer rate constants for the double mutant than for WT KSI. The electrostatic stabilization of the dienolate intermediate relative to the reactant was greater for WT KSI than for the mutants, providing a qualitative explanation for the significantly reduced rates of the mutants. The active site exhibited restricted motion during the proton transfer reactions, but small conformational changes occurred to facilitate the proton transfer reactions by strengthening the hydrogen-bonding interactions and by bringing the proton donor and acceptor closer to each other with the proper orientation for proton transfer. Thus, these calculations suggest that KSI forms a preorganized active site but that the structure of this preorganized active site is altered upon mutation. Moreover, small conformational changes due to stochastic thermal motions are required within this preorganized active site to facilitate the proton transfer reactions.

[1]  A. Mildvan,et al.  Hydrogen Bonding at the Active Site of Δ5-3-Ketosteroid Isomerase† , 1997 .

[2]  B. Oh,et al.  Asp-99 donates a hydrogen bond not to Tyr-14 but to the steroid directly in the catalytic mechanism of Delta 5-3-ketosteroid isomerase from Pseudomonas putida biotype B. , 2000, Biochemistry.

[3]  Theoretical investigation of the role of hydrogen bonding during ketosteroid isomerase catalysis , 2000 .

[4]  J. Åqvist,et al.  The catalytic power of ketosteroid isomerase investigated by computer simulation. , 2002, Biochemistry.

[5]  Kim F. Wong,et al.  Calculation of the transition state theory rate constant for a general reaction coordinate: Application to hydride transfer in an enzyme , 2006 .

[6]  Arieh Warshel,et al.  Electrostatic contributions to binding of transition state analogues can be very different from the corresponding contributions to catalysis: phenolates binding to the oxyanion hole of ketosteroid isomerase. , 2007, Biochemistry.

[7]  Richard Wolfenden,et al.  Analog approaches to the structure of the transition state in enzyme reactions , 1972 .

[8]  D. Shortle,et al.  Kinetic and ultraviolet spectroscopic studies of active-site mutants of delta 5-3-ketosteroid isomerase. , 1989, Biochemistry.

[9]  A. Mildvan,et al.  NMR evidence for the participation of a low-barrier hydrogen bond in the mechanism of delta 5-3-ketosteroid isomerase. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[10]  M. Massiah,et al.  Solution structure of Delta 5-3-ketosteroid isomerase complexed with the steroid 19-nortestosterone hemisuccinate. , 1999, Biochemistry.

[11]  David Chandler,et al.  Barrier crossings:. classical theory of rare but important events , 1998 .

[12]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[13]  S. Benkovic,et al.  Free-energy landscape of enzyme catalysis. , 2008, Biochemistry.

[14]  Gregory K. Schenter,et al.  Generalized transition state theory in terms of the potential of mean force , 2003 .

[15]  David J Weber,et al.  Quantitative interpretations of double mutations of enzymes. , 1992, Archives of biochemistry and biophysics.

[16]  Daniel Herschlag,et al.  Testing Electrostatic Complementarity in Enzyme Catalysis: Hydrogen Bonding in the Ketosteroid Isomerase Oxyanion Hole , 2006, PLoS biology.

[17]  Sharon Hammes-Schiffer,et al.  Hybrid quantum/classical path integral approach for simulation of hydrogen transfer reactions in enzymes. , 2006, The Journal of chemical physics.

[18]  Nam-Chul Ha,et al.  Structural double-mutant cycle analysis of a hydrogen bond network in ketosteroid isomerase from Pseudomonas putida biotype B. , 2004, The Biochemical journal.

[19]  S. Benkovic,et al.  Preorganization and protein dynamics in enzyme catalysis. , 2002, Chemical record.

[20]  M. Summers,et al.  Mechanistic Insights from the Three-Dimensional Structure of 3-Oxo-Δ5-steroid Isomerase , 1999 .

[21]  Barry Honig,et al.  Extending the accuracy limits of prediction for side-chain conformations. , 2001 .

[22]  G. Petsko,et al.  Hydrogen bond coupling in the ketosteroid isomerase active site. , 2009, Biochemistry.

[23]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[24]  R. Pollack,et al.  Electrophilic Assistance by Asp-99 of 3-Oxo-Δ5-steroid Isomerase† , 1998 .

[25]  R. Pollack,et al.  Energetics of 3-oxo-delta 5-steroid isomerase: source of the catalytic power of the enzyme. , 1991, Biochemistry.

[26]  James B. Anderson,et al.  Statistical theories of chemical reactions. Distributions in the transition region , 1973 .

[27]  Computational modeling of enzymatic keto-enol isomerization reactions , 2002 .

[28]  Do ligand binding and solvent exclusion alter the electrostatic character within the oxyanion hole of an enzymatic active site? , 2007, Journal of the American Chemical Society.

[29]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[30]  A. Mildvan,et al.  Ultraviolet spectroscopic evidence for decreased motion of the active site tyrosine residue of delta 5-3-ketosteroid isomerase by steroid binding. , 1995, Biochemistry.

[31]  Ralph E. Christoffersen,et al.  Algorithms for Chemical Computations , 1977 .

[32]  B. Oh,et al.  Role of catalytic residues in enzymatic mechanisms of homologous ketosteroid isomerases. , 2000, Biochemistry.

[33]  G. T. Marks,et al.  Short, strong hydrogen bonds on enzymes: NMR and mechanistic studies , 2002 .

[34]  M. S. Kim,et al.  Contribution of the hydrogen-bond network involving a tyrosine triad in the active site to the structure and function of a highly proficient ketosteroid isomerase from Pseudomonas putida biotype B. , 2000, Biochemistry.

[35]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[36]  Dhruva K. Chakravorty,et al.  Hybrid quantum/classical molecular dynamics simulations of the proton transfer reactions catalyzed by ketosteroid isomerase: analysis of hydrogen bonding, conformational motions, and electrostatics. , 2009, Biochemistry.

[37]  A. Mildvan,et al.  Studies of the catalytic mechanism of an active-site mutant (Y14F) of delta 5-3-ketosteroid isomerase by kinetic deuterium isotope effects. , 1991, Biochemistry.

[38]  A. Mildvan,et al.  13C NMR relaxation studies of backbone and side chain motion of the catalytic tyrosine residue in free and steroid-bound delta 5-3-ketosteroid isomerase. , 1996, Biochemistry.

[39]  Arieh Warshel,et al.  Computer Modeling of Chemical Reactions in Enzymes and Solutions , 1991 .

[40]  Arieh Warshel,et al.  Ketosteroid isomerase provides further support for the idea that enzymes work by electrostatic preorganization , 2010, Proceedings of the National Academy of Sciences.

[41]  G. Ciccotti,et al.  Constrained reaction coordinate dynamics for the simulation of rare events , 1989 .

[42]  K. S. Kim,et al.  Catalytic role of enzymes: short strong H-bond-induced partial proton shuttles and charge redistributions. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[43]  K. Y. Choi,et al.  Identification of active site residues by site-directed mutagenesis of delta 5-3-ketosteroid isomerase from Pseudomonas putida biotype B , 1995, Journal of bacteriology.

[44]  T. C. Bruice,et al.  Computational study of ketosteroid isomerase: insights from molecular dynamics simulation of enzyme bound substrate and intermediate. , 2003, Journal of the American Chemical Society.

[45]  Michael F. Summers,et al.  Solution Structure of 3-Oxo-Δ5-Steroid Isomerase , 1997 .

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

[47]  Walter Thiel,et al.  Analysis of the statistical error in umbrella sampling simulations by umbrella integration. , 2006, The Journal of chemical physics.

[48]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[49]  R. Pollack,et al.  Evaluation of the internal equilibrium constant for 3-oxo-delta 5-steroid isomerase using the D38E and D38N mutants: the energetic basis for catalysis. , 1994, Biochemistry.

[50]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[51]  D. Herschlag,et al.  Evaluating the potential for halogen bonding in the oxyanion hole of ketosteroid isomerase using unnatural amino acid mutagenesis. , 2009, ACS chemical biology.

[52]  Arieh Warshel,et al.  A Quantized Classical Path Approach for Calculations of Quantum Mechanical Rate Constants , 1993 .

[53]  Giovanni Ciccotti,et al.  Book Review: Classical and Quantum Dynamics in Condensed Phase Simulations , 1998 .

[54]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

[55]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[56]  Daniel Herschlag,et al.  Testing geometrical discrimination within an enzyme active site: constrained hydrogen bonding in the ketosteroid isomerase oxyanion hole. , 2008, Journal of the American Chemical Society.

[57]  Arieh Warshel,et al.  How Important Are Quantum Mechanical Nuclear Motions in Enzyme Catalysis , 1996 .

[58]  Arieh Warshel,et al.  Simulations of quantum mechanical corrections for rate constants of hydride-transfer reactions in enzymes and solutions , 1991 .

[59]  Michael W. Mahoney,et al.  A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions , 2000 .

[60]  Martin Karplus,et al.  A POSITION DEPENDENT FRICTION MODEL FOR SOLUTION REACTIONS IN THE HIGH FRICTION REGIME : PROTON TRANSFER IN TRIOSEPHOSPHATE ISOMERASE (TIM) , 1996 .

[61]  A. Mildvan,et al.  Substrate polarization by residues in Δ5‐3‐ketosteroid isomerase probed by site‐directed mutagenesis and UV resonance Raman spectroscopy , 1992, Protein science : a publication of the Protein Society.

[62]  B. Oh,et al.  Crystal Structure of Δ5-3-Ketosteroid Isomerase from Pseudomonas testosteroni in Complex with Equilenin Settles the Correct Hydrogen Bonding Scheme for Transition State Stabilization* , 1999, The Journal of Biological Chemistry.

[63]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[64]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[65]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[66]  Walter Thiel,et al.  Bridging the gap between thermodynamic integration and umbrella sampling provides a novel analysis method: "Umbrella integration". , 2005, The Journal of chemical physics.

[67]  Dhruva K. Chakravorty,et al.  Implementation of umbrella integration within the framework of the empirical valence bond approach. , 2008, Journal of chemical theory and computation.

[68]  D. Herschlag,et al.  Dissecting the paradoxical effects of hydrogen bond mutations in the ketosteroid isomerase oxyanion hole , 2010, Proceedings of the National Academy of Sciences.

[69]  James C. Keck,et al.  Variational Theory of Chemical Reaction Rates Applied to Three‐Body Recombinations , 1960 .