1
|
El-Hendawy MM, Garate JA, English NJ, O'Reilly S, Mooney DA. Diffusion and interactions of carbon dioxide and oxygen in the vicinity of the active site of Rubisco: Molecular dynamics and quantum chemical studies. J Chem Phys 2012; 137:145103. [DOI: 10.1063/1.4757021] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Morad M El-Hendawy
- SFI Strategic Research Cluster in Solar Energy Conversion, University College Dublin, Belfield, Dublin 4, Ireland
| | | | | | | | | |
Collapse
|
2
|
Doron D, Major DT, Kohen A, Thiel W, Wu X. Hybrid Quantum and Classical Simulations of the Dihydrofolate Reductase Catalyzed Hydride Transfer Reaction on an Accurate Semi-Empirical Potential Energy Surface. J Chem Theory Comput 2011; 7:3420-37. [PMID: 26598171 DOI: 10.1021/ct2004808] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Dihydrofolate reductase (DHFR) catalyzes the reduction of 7,8-dihydrofolate by nicotinamide adenine dinucleotide phosphate hydride (NADPH) to form 5,6,7,8-tetrahydrofolate and oxidized nicotinamide. DHFR is a small, flexible, monomeric protein with no metals or SS bonds and serves as one of the enzymes commonly used to examine basic aspects in enzymology. In the current work, we present extensive benchmark calculations for several model reactions in the gas phase that are relevant to the DHFR catalyzed hydride transfer. To this end, we employ G4MP2 and CBS-QB3 ab initio calculations as well as numerous density functional theory methods. Using these results, we develop two specific reaction parameter (SRP) Hamiltonians based on the semiempirical AM1 method. The first generation SRP Hamiltonian does not account for dispersion, while the second generation SRP accounts for dispersion implicitly via the AM1 core-repulsion functions. These SRP semiempirical Hamiltonians are subsequently used in hybrid quantum mechanics/molecular mechanics simulations of the DHFR catalyzed reaction. Finally, kinetic isotope effects are computed using a mass-perturbation-based path-integral approach.
Collapse
Affiliation(s)
- Dvir Doron
- Department of Chemistry, The Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Dan Thomas Major
- Department of Chemistry, The Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Amnon Kohen
- Department of Chemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Xin Wu
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
3
|
Jung J, Park HJ, Uhm KN, Kim D, Kim HK. Asymmetric synthesis of (S)-ethyl-4-chloro-3-hydroxy butanoate using a Saccharomyces cerevisiae reductase: enantioselectivity and enzyme-substrate docking studies. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2010; 1804:1841-9. [PMID: 20601218 DOI: 10.1016/j.bbapap.2010.06.011] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/02/2010] [Revised: 06/05/2010] [Accepted: 06/14/2010] [Indexed: 11/25/2022]
Abstract
Ethyl (S)-4-chloro-3-hydroxy butanoate (ECHB) is a building block for the synthesis of hypercholesterolemia drugs. In this study, various microbial reductases have been cloned and expressed in Escherichia coli. Their reductase activities toward ethyl-4-chloro oxobutanoate (ECOB) have been assayed. Amidst them, Baker's yeast YDL124W, YOR120W, and YOL151W reductases showed high activities. YDL124W produced (S)-ECHB exclusively, whereas YOR120W and YOL151W made (R)-form alcohol. The homology models and docking models with ECOB and NADPH elucidated their substrate specificities and enantioselectivities. A glucose dehydrogenase-coupling reaction was used as NADPH recycling system to perform continuously the reduction reaction. Recombinant E. coli cell co-expressing YDL124W and Bacillus subtilis glucose dehydrogenase produced (S)-ECHB exclusively.
Collapse
Affiliation(s)
- Jihye Jung
- Division of Biotechnology, The Catholic University of Korea, Bucheon 420-743, Republic of Korea
| | | | | | | | | |
Collapse
|
4
|
Giese TJ, York DM. Charge-dependent model for many-body polarization, exchange, and dispersion interactions in hybrid quantum mechanical/molecular mechanical calculations. J Chem Phys 2008; 127:194101. [PMID: 18035873 DOI: 10.1063/1.2778428] [Citation(s) in RCA: 54] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This work explores a new charge-dependent energy model consisting of van der Waals and polarization interactions between the quantum mechanical (QM) and molecular mechanical (MM) regions in a combined QMMM calculation. van der Waals interactions are commonly treated using empirical Lennard-Jones potentials, whose parameters are often chosen based on the QM atom type (e.g., based on hybridization or specific covalent bonding environment). This strategy for determination of QMMM nonbonding interactions becomes tedious to parametrize and lacks robust transferability. Problems occur in the study of chemical reactions where the "atom type" is a complex function of the reaction coordinate. This is particularly problematic for reactions, where atoms or localized functional groups undergo changes in charge state and hybridization. In the present work we propose a new model for nonelectrostatic nonbonded interactions in QMMM calculations that overcomes many of these problems. The model is based on a scaled overlap model for repulsive exchange and attractive dispersion interactions that is a function of atomic charge. The model is chemically significant since it properly correlates atomic size, softness, polarizability, and dispersion terms with minimal one-body parameters that are functions of the atomic charge. Tests of the model are examined for rare-gas interactions with neutral and charged atoms in order to demonstrate improved transferability. The present work provides a new framework for modeling QMMM interactions with improved accuracy and transferability.
Collapse
Affiliation(s)
- Timothy J Giese
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA
| | | |
Collapse
|
5
|
Khavrutskii IV, Price DJ, Lee J, Brooks CL. Conformational change of the methionine 20 loop of Escherichia coli dihydrofolate reductase modulates pKa of the bound dihydrofolate. Protein Sci 2007; 16:1087-100. [PMID: 17473015 PMCID: PMC2206655 DOI: 10.1110/ps.062724307] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2006] [Revised: 03/05/2007] [Accepted: 03/06/2007] [Indexed: 10/23/2022]
Abstract
We evaluate the pK(a) of dihydrofolate (H(2)F) at the N(5) position in three ternary complexes with Escherichia coli dihydrofolate reductase (ecDHFR), namely ecDHFR(NADP(+):H(2)F) in the closed form (1), and the Michaelis complexes ecDHFR(NADPH:H(2)F) in the closed (2) and occluded (3) forms, by performing free energy perturbation with molecular dynamics simulations (FEP/MD). Our simulations suggest that in the Michaelis complex the pK(a) is modulated by the Met20 loop fluctuations, providing the largest pK(a) shift in substates with a "tightly closed" loop conformation; in the "partially closed/open" substates, the pK(a) is similar to that in the occluded complex. Conducive to the protonation, tightly closing the Met20 loop enhances the interactions of the cofactor and the substrate with the Met20 side chain and aligns the nicotinamide ring of the cofactor coplanar with the pterin ring of the substrate. Overall, the present study favors the hypothesis that N(5) is protonated directly from solution and provides further insights into the mechanism of the substrate protonation.
Collapse
Affiliation(s)
- Ilja V Khavrutskii
- The Scripps Research Institute, Department of Molecular Biology, TPC6, La Jolla, California 92037, USA
| | | | | | | |
Collapse
|
6
|
Cummins PL, Rostov IV, Gready JE. Calculation of a Complete Enzymic Reaction Surface: Reaction and Activation Free Energies for Hydride-Ion Transfer in Dihydrofolate Reductase. J Chem Theory Comput 2007; 3:1203-11. [DOI: 10.1021/ct600313b] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Peter L. Cummins
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Ivan V. Rostov
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Jill E. Gready
- Computational Proteomics Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| |
Collapse
|
7
|
Balog E, Smith JC, Perahia D. Conformational heterogeneity and low-frequency vibrational modes of proteins. Phys Chem Chem Phys 2006; 8:5543-8. [PMID: 17136269 DOI: 10.1039/b610075a] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Molecular dynamics simulation and normal mode analysis are used to calculate the vibrational density of states of dihydrofolate reductase complexed with nicotinamide adenine dinucleotide phosphate at 120 K and the results are compared with the experimental spectrum derived from inelastic neutron scattering. The simulation results indicate that the experimental spectrum arises from an average over proteins trapped in different conformations with structural differences mainly in the loop regions, and that these conformations have significantly different low-frequency (<20 cm(-1)) spectra. Thus, the experimentally measured spectrum is an average over the vibrational modes of different protein conformations and is thus inhomogeneously broadened. The implications of this broadening for future neutron scattering experiments and ligand binding calculations are discussed.
Collapse
Affiliation(s)
- Erika Balog
- Institut des Hautes Etudes Scientifiques, 35 route de Chartres, 91440, Bures sur Yvette, France
| | | | | |
Collapse
|
8
|
Sergi A, Watney JB, Wong KF, Hammes-Schiffer S. Freezing a Single Distal Motion in Dihydrofolate Reductase. J Phys Chem B 2006; 110:2435-41. [PMID: 16471835 DOI: 10.1021/jp056939u] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Constraining a single motion between distal residues separated by approximately 28 A in hybrid quantum/classical molecular dynamics simulations is found to increase the free energy barrier for hydride transfer in dihydrofolate reductase by approximately 3 kcal/mol. Our analysis indicates that a single distal constraint alters equilibrium motions throughout the enzyme on a wide range of time scales. This alteration of the conformational sampling of the entire system is sufficient to significantly increase the free energy barrier and decrease the rate of hydride transfer. Despite the changes in conformational sampling introduced by the constraint, the system assumes a similar transition state conformation with a donor-acceptor distance of approximately 2.72 A to enable the hydride transfer reaction. The modified thermal sampling leads to a substantial increase in the average donor-acceptor distance for the reactant state, however, thereby decreasing the probability of sampling the transition state conformations with the shorter distances required for hydride transfer. These simulations indicate that fast thermal fluctuations of the enzyme, substrate, and cofactor lead to conformational sampling of configurations that facilitate hydride transfer. The fast thermal motions are in equilibrium as the reaction progresses along the collective reaction coordinate, and the overall average equilibrium conformational changes occur on the slower time scale measured experimentally. Recent single molecule experiments suggest that at least some of these thermally averaged equilibrium conformational changes occur on the millisecond time scale of the hydride transfer reaction. Thus, introducing a constraint that modifies the conformational sampling of an enzyme could significantly impact its catalytic activity.
Collapse
Affiliation(s)
- Alessandro Sergi
- Department of Chemistry, 104 Chemistry Building, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | | | | | | |
Collapse
|
9
|
Garcia-Viloca M, Truhlar DG, Gao J. Reaction-path energetics and kinetics of the hydride transfer reaction catalyzed by dihydrofolate reductase. Biochemistry 2004; 42:13558-75. [PMID: 14622003 DOI: 10.1021/bi034824f] [Citation(s) in RCA: 179] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have studied the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) and the coenzyme nicotinamide adenine dinucleotide phosphate (NADPH); the substrate is 5-protonated 7,8-dihydrofolate, and the product is tetrahydrofolate. The potential energy surface is modeled by a combined quantum mechanical-molecular mechanical (QM/MM) method employing Austin model 1 (AM1) and a simple valence bond potential for 69 QM atoms and employing the CHARMM22 and TIP3P molecular mechanics force fields for the other 21 399 atoms; the QM and MM regions are joined by two boundary atoms treated by the generalized hybrid orbital (GHO) method. All simulations are carried out using periodic boundary conditions at neutral pH and 298 K. In stage 1, a reaction coordinate is defined as the difference between the breaking and forming bond distances to the hydride ion, and a quasithermodynamic free energy profile is calculated along this reaction coordinate. This calculation includes quantization effects on bound vibrations but not on the reaction coordinate, and it is used to locate the variational transition state that defines a transition state ensemble. Then, the key interactions at the reactant, variational transition state, and product are analyzed in terms of both bond distances and electrostatic energies. The results of both analyses support the conclusion derived from previous mutational studies that the M20 loop of DHFR makes an important contribution to the electrostatic stabilization of the hydride transfer transition state. Third, transmission coefficients (including recrossing factors and multidimensional tunneling) are calculated and averaged over the transition state ensemble. These averaged transmission coefficients, combined with the quasithermodynamic free energy profile determined in stage 1, allow us to calculate rate constants, phenomenological free energies of activation, and primary and secondary kinetic isotope effects. A primary kinetic isotope effect (KIE) of 2.8 has been obtained, in good agreement with the experimentally determined value of 3.0 and with the value 3.2 calculated previously. The primary KIE is mainly a consequence of the quantization of bound vibrations. In contrast, the secondary KIE, with a value of 1.13, is almost entirely due to dynamical effects on the reaction coordinate, especially tunneling.
Collapse
Affiliation(s)
- Mireia Garcia-Viloca
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | | | | |
Collapse
|
10
|
Thorpe IF, Brooks CL. Barriers to Hydride Transfer in Wild Type and Mutant Dihydrofolate Reductase from E. coli. J Phys Chem B 2003. [DOI: 10.1021/jp035734n] [Citation(s) in RCA: 53] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ian F. Thorpe
- Department of Molecular Biology (TPC6), Center for Theoretical Biological Physics, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037
| | - Charles L. Brooks
- Department of Molecular Biology (TPC6), Center for Theoretical Biological Physics, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037
| |
Collapse
|
11
|
|
12
|
Computational methods for the study of enzymic reaction mechanisms. II. An overlapping mechanically embedded method for hybrid semi-empirical-QM/MM calculations. ACTA ACUST UNITED AC 2003. [DOI: 10.1016/s0166-1280(03)00303-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
13
|
Watney JB, Agarwal PK, Hammes-Schiffer S. Effect of mutation on enzyme motion in dihydrofolate reductase. J Am Chem Soc 2003; 125:3745-50. [PMID: 12656604 DOI: 10.1021/ja028487u] [Citation(s) in RCA: 97] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Hybrid quantum-classical molecular dynamics simulations of a mutant Escherichia coli dihydrofolate reductase enzyme are presented. Although residue 121 is on the exterior of the enzyme, experimental studies have shown that the mutation of Gly-121 to valine reduces the rate of hydride transfer by a factor of 163. The simulations indicate that the decrease in the hydride transfer rate for the G121V mutant is due to an increase in the free energy barrier. The calculated free energy barrier is higher for the mutant than for the wild-type enzyme by an amount that is consistent with the experimentally observed rate reduction. The calculated transmission coefficients are comparable for the wild-type and mutant enzymes. The simulations suggest that this mutation may interrupt a network of coupled promoting motions proposed to play an important role in DHFR catalysis. This phenomenon has broad implications for protein engineering and drug design.
Collapse
Affiliation(s)
- James B Watney
- Department of Chemistry, 152 Davey Laboratory, Pennsylvania State University, University Park 16802, USA
| | | | | |
Collapse
|
14
|
Garcia-Viloca M, Truhlar DG, Gao J. Importance of substrate and cofactor polarization in the active site of dihydrofolate reductase. J Mol Biol 2003; 327:549-60. [PMID: 12628257 DOI: 10.1016/s0022-2836(03)00123-2] [Citation(s) in RCA: 73] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
By using a combined quantum-mechanical and molecular-mechanical potential in molecular dynamics simulations, we have investigated the effects of the enzyme electric field of dihydrofolate reductase on the electronic polarization of its 5-protonated dihydrofolate substrate at various stages of the catalyzed hydride transfer reaction. Energy decomposition of the total electrostatic interaction energy between the ligands and the enzyme shows that the polarization effect is 4% of the total electrostatic interaction energy, and, significantly, it accounts for 9kcal/mol of transition state stabilization relative to the reactant state. Therefore it is essential to take account of substrate polarization for quantitative interpretation of enzymatic function and for calculation of binding free energies of inhibitors to a protein. Atomic polarizations are calculated as the differences in the average atomic charges on the atoms in gas phase and in molecular simulations of the enzyme; this analysis shows that the glutamate tail and the pterin ring are the highly polarized regions of the substrate. Electron density difference plots of the reactant and product complexes at instantaneous configurations in the enzyme active center confirm the inferences made on the basis of partial atomic charges.
Collapse
Affiliation(s)
- Mireia Garcia-Viloca
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street, SE Minneapolis 55455-0431, USA.
| | | | | |
Collapse
|
15
|
Titmuss SJ, Cummins PL, Rendell AP, Bliznyuk AA, Gready JE. Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis. J Comput Chem 2002; 23:1314-22. [PMID: 12214314 DOI: 10.1002/jcc.10122] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
QM/MM methods have been developed as a computationally feasible solution to QM simulation of chemical processes, such as enzyme-catalyzed reactions, within a more approximate MM representation of the condensed-phase environment. However, there has been no independent method for checking the quality of this representation, especially for highly nonisotropic protein environments such as those surrounding enzyme active sites. Hence, the validity of QM/MM methods is largely untested. Here we use the possibility of performing all-QM calculations at the semiempirical PM3 level with a linear-scaling method (MOZYME) to assess the performance of a QM/MM method (PM3/AMBER94 force field). Using two model pathways for the hydride-ion transfer reaction of the enzyme dihydrofolate reductase studied previously (Titmuss et al., Chem Phys Lett 2000, 320, 169-176), we have analyzed the reaction energy contributions (QM, QM/MM, and MM) from the QM/MM results and compared them with analogous-region components calculated via an energy partitioning scheme implemented into MOZYME. This analysis further divided the MOZYME components into Coulomb, resonance and exchange energy terms. For the model in which the MM coordinates are kept fixed during the reaction, we find that the MOZYME and QM/MM total energy profiles agree very well, but that there are significant differences in the energy components. Most significantly there is a large change (approximately 16 kcal/mol) in the MOZYME MM component due to polarization of the MM region surrounding the active site, and which arises mostly from MM atoms close to (<10 A) the active-site QM region, which is not modelled explicitly by our QM/MM method. However, for the model where the MM coordinates are allowed to vary during the reaction, we find large differences in the MOZYME and QM/MM total energy profiles, with a discrepancy of 52 kcal/mol between the relative reaction (product-reactant) energies. This is largely due to a difference in the MM energies of 58 kcal/mol, of which we can attribute approximately 40 kcal/mol to geometry effects in the MM region and the remainder, as before, to MM region polarization. Contrary to the fixed-geometry model, there is no correlation of the MM energy changes with distance from the QM region, nor are they contributed by only a few residues. Overall, the results suggest that merely extending the size of the QM region in the QM/MM calculation is not a universal solution to the MOZYME- and QM/MM-method differences. They also suggest that attaching physical significance to MOZYME Coulomb, resonance and exchange components is problematic. Although we conclude that it would be possible to reparameterize the QM/MM force field to reproduce MOZYME energies, a better way to account for both the effects of the protein environment and known deficiencies in semiempirical methods would be to parameterize the force field based on data from DFT or ab initio QM linear-scaling calculations. Such a force field could be used efficiently in MD simulations to calculate free energies.
Collapse
Affiliation(s)
- Stephen J Titmuss
- Computational Proteomics and Therapy Design Group, Division of Biochemistry and Molecular Biology, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra, ACT 2601, Australia
| | | | | | | | | |
Collapse
|
16
|
Cummins PL, Greatbanks SP, Rendell AP, Gready JE. Computational Methods for the Study of Enzymic Reaction Mechanisms. 1. Application to the Hydride Transfer Step in the Catalysis of Dihydrofolate Reductase. J Phys Chem B 2002. [DOI: 10.1021/jp021070q] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Peter L. Cummins
- Computational Proteomics and Therapy Design Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia, and Department of Computer Science, Faculty of Engineering and Information Technology, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Stephen P. Greatbanks
- Computational Proteomics and Therapy Design Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia, and Department of Computer Science, Faculty of Engineering and Information Technology, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Alistair P. Rendell
- Computational Proteomics and Therapy Design Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia, and Department of Computer Science, Faculty of Engineering and Information Technology, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| | - Jill E. Gready
- Computational Proteomics and Therapy Design Group, John Curtin School of Medical Research, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia, and Department of Computer Science, Faculty of Engineering and Information Technology, Australian National University, P.O. Box 334, Canberra ACT 2601, Australia
| |
Collapse
|
17
|
Agarwal PK, Billeter SR, Hammes-Schiffer S. Nuclear Quantum Effects and Enzyme Dynamics in Dihydrofolate Reductase Catalysis. J Phys Chem B 2002. [DOI: 10.1021/jp020190v] [Citation(s) in RCA: 156] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Pratul K. Agarwal
- Department of Chemistry, 152 Davey Laboratory, Pennsylvania State University, University Park, Pennsylvania 16802
| | - Salomon R. Billeter
- Department of Chemistry, 152 Davey Laboratory, Pennsylvania State University, University Park, Pennsylvania 16802
| | - Sharon Hammes-Schiffer
- Department of Chemistry, 152 Davey Laboratory, Pennsylvania State University, University Park, Pennsylvania 16802
| |
Collapse
|
18
|
Dybala-Defratyka A, Paneth P. Theoretical evaluation of the hydrogen kinetic isotope effect on the first step of the methylmalonyl-CoA mutase reaction. J Inorg Biochem 2001; 86:681-9. [PMID: 11583786 DOI: 10.1016/s0162-0134(01)00230-6] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
We have calculated hydrogen kinetic isotope effects (KIEs) for the first step of the methylmalonyl-CoA mutase reaction, including multidimensional tunneling correction at the zero curvature (ZCT) level, and compared them with the experimental values. Both alternative mechanisms of this step, concerted and stepwise, can be accommodated. It turned out to be essential to include Arg207 hydrogen-bonded to the reactant in the mechanism predicting simultaneous breaking of the Co-C bond of AdoCbl and hydrogen atom transfer. The consequence of the stepwise mechanism is a much larger facilitation of the homolytic dissociation of the carbon-cobalt bond by the enzyme than currently appreciated; our results suggest lowering of the activation energy by about 23 kcal mol(-1). We have also shown that large hydrogen KIEs of tunneling origin do not necessarily break the Swain-Schaad equation. Furthermore, when this equation does not hold, the exponent may be smaller in the presence of tunneling than it is at the semi-classical limit, indicating that nonclassical behavior may be a more common phenomenon than expected.
Collapse
Affiliation(s)
- A Dybala-Defratyka
- Institute of Applied Radiation Chemistry, Department of Chemistry, Technical University of Lodz, Zeromskiego 116, 90-924, Lodz, Poland
| | | |
Collapse
|
19
|
Wang W, Donini O, Reyes CM, Kollman PA. Biomolecular simulations: recent developments in force fields, simulations of enzyme catalysis, protein-ligand, protein-protein, and protein-nucleic acid noncovalent interactions. ANNUAL REVIEW OF BIOPHYSICS AND BIOMOLECULAR STRUCTURE 2001; 30:211-43. [PMID: 11340059 DOI: 10.1146/annurev.biophys.30.1.211] [Citation(s) in RCA: 389] [Impact Index Per Article: 16.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Computer modeling has been developed and widely applied in studying molecules of biological interest. The force field is the cornerstone of computer simulations, and many force fields have been developed and successfully applied in these simulations. Two interesting areas are (a) studying enzyme catalytic mechanisms using a combination of quantum mechanics and molecular mechanics, and (b) studying macromolecular dynamics and interactions using molecular dynamics (MD) and free energy (FE) calculation methods. Enzyme catalysis involves forming and breaking of covalent bonds and requires the use of quantum mechanics. Noncovalent interactions appear ubiquitously in biology, but here we confine ourselves to review only noncovalent interactions between protein and protein, protein and ligand, and protein and nucleic acids.
Collapse
Affiliation(s)
- W Wang
- Graduate Group in Biophysics, University of California San Francisco, California 94143, USA.
| | | | | | | |
Collapse
|
20
|
Lewandowicz A, Rudziński J, Tronstad L, Widersten M, Ryberg P, Matsson O, Paneth P. Chlorine kinetic isotope effects on the haloalkane dehalogenase reaction. J Am Chem Soc 2001; 123:4550-5. [PMID: 11457241 DOI: 10.1021/ja003503d] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have found chlorine kinetic isotope effects on the dehalogenation catalyzed by haloalkane dehalogenase from Xanthobacter autotrophicus GJ10 to be 1.0045 +/- 0.0004 for 1,2-dichloroethane and 1.0066 +/- 0.0004 for 1-chlorobutane. The latter isotope effect approaches the intrinsic chlorine kinetic isotope effect for the dehalogenation step. The intrinsic isotope effect has been modeled using semiempirical and DFT theory levels using the ONIOM QM/QM scheme. Our results indicate that the dehalogenation step is reversible; the overall irreversibility of the enzyme-catalyzed reaction is brought about by a step following the dehalogenation.
Collapse
Affiliation(s)
- A Lewandowicz
- Department of Chemistry, Technical University of Lodz, Zeromskiego 116, 90-924 Lodz, Poland
| | | | | | | | | | | | | |
Collapse
|
21
|
Castillo R, Andrés J, Moliner V. Quantum Mechanical/Molecular Mechanical Study on the Favorskii Rearrangement in Aqueous Media. J Phys Chem B 2001. [DOI: 10.1021/jp003264g] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- R. Castillo
- Departament de Ciències Experimentals, Universitat Jaume I, Box 224, 12080 Castelló, Spain
| | - J. Andrés
- Departament de Ciències Experimentals, Universitat Jaume I, Box 224, 12080 Castelló, Spain
| | - V. Moliner
- Departament de Ciències Experimentals, Universitat Jaume I, Box 224, 12080 Castelló, Spain
| |
Collapse
|
22
|
|
23
|
Orozco M, Luque FJ. Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. Chem Rev 2000; 100:4187-4226. [PMID: 11749344 DOI: 10.1021/cr990052a] [Citation(s) in RCA: 454] [Impact Index Per Article: 18.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Modesto Orozco
- Departament de Bioquímica i Biologia Molecular, Facultat de Química, Universitat de Barcelona, Martí i Franqués 1, E-08028 Barcelona, Spain, and Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona, Avgda. Diagonal s/n, E-08028 Barcelona, Spain
| | | |
Collapse
|
24
|
Greatbanks SP, Gready JE, Limaye AC, Rendell AP. Comparison of enzyme polarization of ligands and charge-transfer effects for dihydrofolate reductase using point-charge embeddedab initio quantum mechanical and linear-scaling semiempirical quantum mechanical methods. J Comput Chem 2000. [DOI: 10.1002/(sici)1096-987x(20000715)21:9<788::aid-jcc7>3.0.co;2-q] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
25
|
Svanberg M, Pettersson JBC, Bolton K. Coupled QM/MM Molecular Dynamics Simulations of HCl Interacting with Ice Surfaces and Water Clusters − Evidence of Rapid Ionization. J Phys Chem A 2000. [DOI: 10.1021/jp0012698] [Citation(s) in RCA: 96] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marcus Svanberg
- Department of Chemistry, Physical Chemistry, Göteborg University, SE-412 96 Göteborg, Sweden, School of Engineering, University of Borås, SE-501 90 Borås, Sweden, and School of Environmental Sciences, Göteborg University, SE-412 96 Göteborg, Sweden
| | - Jan B. C. Pettersson
- Department of Chemistry, Physical Chemistry, Göteborg University, SE-412 96 Göteborg, Sweden, School of Engineering, University of Borås, SE-501 90 Borås, Sweden, and School of Environmental Sciences, Göteborg University, SE-412 96 Göteborg, Sweden
| | - Kim Bolton
- Department of Chemistry, Physical Chemistry, Göteborg University, SE-412 96 Göteborg, Sweden, School of Engineering, University of Borås, SE-501 90 Borås, Sweden, and School of Environmental Sciences, Göteborg University, SE-412 96 Göteborg, Sweden
| |
Collapse
|
26
|
Cummins PL, Gready JE. Combined Quantum and Molecular Mechanics (QM/MM) Study of the Ionization State of 8-Methylpterin Substrate Bound to Dihydrofolate Reductase. J Phys Chem B 2000. [DOI: 10.1021/jp993153l] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Peter L. Cummins
- Computational Molecular Biology and Drug Design Group, John Curtin School of Medical Research, Australian National University, P.O. BOX 334, Canberra ACT, 2601 Australia
| | - Jill E. Gready
- Computational Molecular Biology and Drug Design Group, John Curtin School of Medical Research, Australian National University, P.O. BOX 334, Canberra ACT, 2601 Australia
| |
Collapse
|
27
|
Titmuss SJ, Cummins PL, Bliznyuk AA, Rendell AP, Gready JE. Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods applied to enzyme reactions. Chem Phys Lett 2000. [DOI: 10.1016/s0009-2614(00)00215-3] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
28
|
Cummins PL, Gready JE. QM/MM and SCRF studies of the ionization state of 8-methylpterin substrate bound to dihydrofolate reductase: existence of a low-barrier hydrogen bond. J Mol Graph Model 2000; 18:42-9. [PMID: 10935206 DOI: 10.1016/s1093-3263(00)00034-6] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Abstract
Using combined semiempirical quantum mechanics and molecular mechanics (QM/MM) and ab initio self-consistent reaction field (SCRF) calculations, we determined that a low-barrier hydrogen bond (LBHB) is formed when the mechanism-based substrate 8-methylpterin binds to dihydrofolate reductase (DHFR). The substrate initially was assumed bound either in the ion-pair form corresponding to N3-protonated substrate hydrogen (H) bonded to the unprotonated (carboxylate) of the conserved Glu30 residue in the active site, or in the neutral-pair form corresponding to unprotonated substrate H bonded to the neutral (carboxylic acid) from of Glu30. The free energy of interaction of these H-bonded systems with the protein/solvent surroundings was computed using a coordinate-coupled free energy perturbation (FEP) method implemented within the molecular dynamics (MD) simulation scheme and using a semiempirical (PM3) QM/MM force field. The free energy obtained from the QM/MM force-field simulations corresponds most closely with the corresponding free energy component obtained from HF/6-31G* SCRF calculations using a value of 2 for the dielectric constant (epsilon) for the solvated protein. Calculations were performed at levels ranging from HF/6-31G to MP2/6-31G* to B3LYP/6-31 + G**, with varying dielectric constants. The energy-minimized path for motion of the proton in the H bond along a one-dimensional reaction coordinate was calculated at HF/6-31G, HF/6-31G* (epsilon = 1) and B3LYP/6-31G* (epsilon = 2) levels. These calculations identified a second neutral-pair complex, involving the 2-amino group of substrate, which also interacts with Glu30, which is lower in energy than the ion-pair form. A harmonic vibrational analysis shows that the first vibrational state appears to lie near or above the TS connecting potential energy minima corresponding to the two neutral-pair configurations, thus indicating an LBHB. Consequently, the H-bonded system will have a significant probability of being found in the ion-pair form, in agreement with experimental spectral studies indicating an enzyme-bound cation and suggesting that the LBHB would activate substrate towards hydride-ion transfer from NADPH.
Collapse
Affiliation(s)
- P L Cummins
- Computational Molecular Biology and Drug Design Group, John Curtin School of Medical Research, Australian National University, Canberra, Australia
| | | |
Collapse
|
29
|
Castillo R, Andrés J, Moliner V. Catalytic Mechanism of Dihydrofolate Reductase Enzyme. A Combined Quantum-Mechanical/Molecular-Mechanical Characterization of Transition State Structure for the Hydride Transfer Step. J Am Chem Soc 1999. [DOI: 10.1021/ja9843019] [Citation(s) in RCA: 67] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- R. Castillo
- Contribution from the Departament de Ciències Experimentals, Universitat Jaume I, Castelló, Spain
| | - J. Andrés
- Contribution from the Departament de Ciències Experimentals, Universitat Jaume I, Castelló, Spain
| | - V. Moliner
- Contribution from the Departament de Ciències Experimentals, Universitat Jaume I, Castelló, Spain
| |
Collapse
|
30
|
Greatbanks SP, Gready JE, Limaye AC, Rendell AP. Enzyme polarization of substrates of dihydrofolate reductase by different theoretical methods. Proteins 1999. [DOI: 10.1002/(sici)1097-0134(19991101)37:2<157::aid-prot2>3.0.co;2-j] [Citation(s) in RCA: 20] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
31
|
Monard G, Merz KM. Combined Quantum Mechanical/Molecular Mechanical Methodologies Applied to Biomolecular Systems. Acc Chem Res 1999. [DOI: 10.1021/ar970218z] [Citation(s) in RCA: 293] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Gérald Monard
- 152 Davey Laboratory, Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802
| | - Kenneth M. Merz
- 152 Davey Laboratory, Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802
| |
Collapse
|
32
|
Cummins PL, Gready JE. Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules. J Comput Chem 1999. [DOI: 10.1002/(sici)1096-987x(19990730)20:10<1028::aid-jcc5>3.0.co;2-7] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
33
|
Boresch S, Ringhofer S, Höchtl P, Steinhauser O. Towards a better description and understanding of biomolecular solvation. Biophys Chem 1999; 78:43-68. [PMID: 17030304 DOI: 10.1016/s0301-4622(98)00235-x] [Citation(s) in RCA: 41] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/1998] [Revised: 11/12/1998] [Accepted: 11/23/1998] [Indexed: 11/26/2022]
Abstract
We introduce a flexible framework for the correct description of the solvation of biological macromolecules, the dielectric field equation (DFE). The formalism permits the use of any combination of quantum mechanical (QM), molecular mechanical (MM) and continuum electrostatic (CE) based techniques. For the CE region a method that yields the electric field rather than the potential is outlined. The DFE formalism makes clear the need to consider and to calibrate a dielectric boundary region surrounding the simulation system. The details of how to do this are presented for the case of the Ewald summation method; the effects are demonstrated by calculations of the dielectric properties and the spatially resolved Kirkwood G-factor, G(K)(r), of TIP3P water. Computing the dielectric properties of a multi-component system provides a sensitive method to better understand the solvation of biological macromolecules. Towards this goal a rigorous analysis of the dielectric properties of solvated biomolecules based on a decomposition of the frequency-dependent dielectric constant (or susceptibility) of the full system is presented. The meaning of our approach is investigated, and the results of a first application are reported. Using the method of Voronoi polyhedra, the dielectric properties of the first two solvation shells and bulk water are compared by re-analyzing a 12-ns trajectory of a zinc finger peptide in water [Löffler et al. J. Mol. Biol. 270 (1997) 520]. It is found that the first shell behaves considerably different; in addition, there is a non-negligible contribution to the total susceptibility of the system from coupling between the protein and the bulk water phase, i.e. the water molecules not in the immediate vicinity of the solute.
Collapse
Affiliation(s)
- S Boresch
- Institut für Theoretische Chemie und Molekulare Strukturbiologie, Universität Wien, Währingerstrasse 17, A-1090 Vienna, Austria
| | | | | | | |
Collapse
|
34
|
|
35
|
King WA, Gready JE, Andrews TJ. Quantum chemical analysis of the enolization of ribulose bisphosphate: the first hurdle in the fixation of CO2 by Rubisco. Biochemistry 1998; 37:15414-22. [PMID: 9799503 DOI: 10.1021/bi981598e] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
A study, using ab initio quantum chemical methods, of the first step in the reaction mechanism of Rubisco, the enolization of the substrate, ribulose bisphosphate, is reported. This is the first such study that takes into account the likely roles of critical features within the active site. On the basis of molecular dynamics relaxation of the complex between activated enzyme and substrate using X-ray crystallographic structures as starting coordinates, a 29-atom fragment that mimicked the active site was constructed. States along a proposed reaction pathway were calculated using density functional theory and Moller-Plesset second-order perturbation theory. The results are consistent with the postulate that the base that abstracts the C3 proton of ribulose bisphosphate is the metal-stabilized carbamate of Lys-201 formed during the activation process. The calculations suggest that the active-site residue, Lys-175, is charged before enolization commences and they indicate a possible means by which the enzyme directs the incoming CO2 to attack the C2 carbon atom of the enediol, rather than the chemically very similar C3 atom.
Collapse
Affiliation(s)
- W A King
- Computational Molecular Biology and Drug Design Group, John Curtin School of Medical Research, Australian National University, Canberra
| | | | | |
Collapse
|