Hybrid Quantum/Classical Molecular Dynamics Simulations of the Proton Transfer Reactions Catalyzed by Ketosteroid Isomerase: Analysis of Hydrogen Bonding, Conformational Motions, and Electrostatics
- 2 October 2009
- journal article
- research article
- Published by American Chemical Society (ACS) in Biochemistry
- Vol. 48 (44) , 10608-10619
Abstract
Hybrid quantum/classical molecular dynamics simulations of the two proton transfer reactions catalyzed by ketosteroid isomerase are presented. The potential energy surfaces for the proton transfer reactions are described with the empirical valence bond method. Nuclear quantum effects of the transferring hydrogen increase the rates by a factor of ∼8, and dynamical barrier recrossings decrease the rates by a factor of 3−4. For both proton transfer reactions, the donor−acceptor distance decreases substantially at the transition state. The carboxylate group of the Asp38 side chain, which serves as the proton acceptor and donor in the first and second steps, respectively, rotates significantly between the two proton transfer reactions. The hydrogen-bonding interactions within the active site are consistent with the hydrogen bonding of both Asp99 and Tyr14 to the substrate. The simulations suggest that a hydrogen bond between Asp99 and the substrate is present from the beginning of the first proton transfer step, whereas the hydrogen bond between Tyr14 and the substrate is virtually absent in the first part of this step but forms nearly concurrently with the formation of the transition state. Both hydrogen bonds are present throughout the second proton transfer step until partial dissociation of the product. The hydrogen bond between Tyr14 and Tyr55 is present throughout both proton transfer steps. The active site residues are more mobile during the first step than during the second step. The van der Waals interaction energy between the substrate and the enzyme remains virtually constant along the reaction pathway, but the electrostatic interaction energy is significantly stronger for the dienolate intermediate than for the reactant and product. Mobile loop regions distal to the active site exhibit significant structural rearrangements and, in some cases, qualitative changes in the electrostatic potential during the catalytic reaction. These results suggest that relatively small conformational changes of the enzyme active site and substrate strengthen the hydrogen bonds that stabilize the intermediate, thereby facilitating the proton transfer reactions. Moreover, the conformational and electrostatic changes associated with these reactions are not limited to the active site but rather extend throughout the entire enzyme.Keywords
This publication has 70 references indexed in Scilit:
- Hydrogen Bond Coupling in the Ketosteroid Isomerase Active SiteBiochemistry, 2009
- Evaluating the Potential for Halogen Bonding in the Oxyanion Hole of Ketosteroid Isomerase Using Unnatural Amino Acid MutagenesisACS Chemical Biology, 2009
- Implementation of Umbrella Integration within the Framework of the Empirical Valence Bond ApproachJournal of Chemical Theory and Computation, 2008
- Testing Geometrical Discrimination within an Enzyme Active Site: Constrained Hydrogen Bonding in the Ketosteroid Isomerase Oxyanion HoleJournal of the American Chemical Society, 2008
- Do Ligand Binding and Solvent Exclusion Alter the Electrostatic Character within the Oxyanion Hole of an Enzymatic Active Site?Journal of the American Chemical Society, 2007
- PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulationsNucleic Acids Research, 2007
- Development and testing of a general amber force fieldJournal of Computational Chemistry, 2004
- VMD: Visual molecular dynamicsJournal of Molecular Graphics, 1996
- A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP modelThe Journal of Physical Chemistry, 1993
- Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics, 1983