Comparison of DFT and ab initio QM/MM methods for modelling reaction in chorismate synthase

Abstract Quantum mechanics/molecular mechanics (QM/MM) methods are a popular tool in the investigation of enzyme reactions. Here, we compare B3LYP density functional theory (DFT) and ab initio QM/MM methods for modelling the conversion of 5-enolpyruvylshikimate-3-phosphate to chorismate in chorismate synthase. Good agreement with experimental data is only obtained at the SCS-MP2/CHARMM27 level for a reaction mechanism in which phosphate elimination precedes proton transfer. B3LYP predicts reaction energetics that are qualitatively wrong, stressing the need for ab initio QM/MM methods, and caution in interpretation of DFT results for this enzyme.

[1]  F. Lederer.  Extreme pKa displacements at the active sites of FMN‐dependent α‐hydroxy acid‐oxidizing enzymes , 1992 .

[2]  Jan H. Jensen,et al.  Very fast prediction and rationalization of pKa values for protein–ligand complexes , 2008, Proteins.

[3]  The transient kinetics of Escherichia coli chorismate synthase: substrate consumption, product formation, phosphate dissociation, and characterization of a flavin intermediate. , 1996, Biochemistry.

[4]  Lars Ridder,et al.  Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction. , 2004, Organic & biomolecular chemistry.

[5]  Nathalie Reuter,et al.  Frontier Bonds in QM/MM Methods: A Comparison of Different Approaches , 2000 .

[6]  S. Bornemann,et al.  Studies with Substrate and Cofactor Analogues Provide Evidence for a Radical Mechanism in the Chorismate Synthase Reaction* , 2000, The Journal of Biological Chemistry.

[7]  G Vriend,et al.  WHAT IF: a molecular modeling and drug design program. , 1990, Journal of molecular graphics.

[8]  Thom H. Dunning,et al.  Gaussian basis sets for use in correlated molecular calculations. V. Core-valence basis sets for boron through neon , 1995 .

[9]  F. Lederer.  Extreme pKa displacements at the active sites of FMN-dependent alpha-hydroxy acid-oxidizing enzymes. , 1992, Protein Science.

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

[11]  Adrian J. Mulholland,et al.  Investigations of enzyme-catalysed reactions with combined quantum mechanics/molecular mechanics (QM/MM) methods , 2010 .

[12]  S. Bornemann,et al.  Replacement of two invariant serine residues in chorismate synthase provides evidence that a proton relay system is essential for intermediate formation and catalytic activity , 2008, The FEBS journal.

[13]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[14]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[15]  J. Maclean,et al.  The structure of chorismate synthase reveals a novel flavin binding site fundamental to a unique chemical reaction. , 2003, Structure.

[16]  L. Hedstrom,et al.  Is Arg418 the catalytic base required for the hydrolysis step of the IMP dehydrogenase reaction? , 2005, Biochemistry.

[17]  Adrian J Mulholland,et al.  Atomic Description of an Enzyme Reaction Dominated by Proton Tunneling , 2006, Science.

[18]  S. Bornemann.  Flavoenzymes that catalyse reactions with no net redox change. , 2002, Natural product reports.

[19]  S. Bornemann,et al.  Studies with Flavin Analogs Provide Evidence That a Protonated Reduced FMN Is the Substrate-induced Transient Intermediate in the Reaction of Escherichia coli Chorismate Synthase* , 1996, The Journal of Biological Chemistry.

[20]  Adrian J Mulholland,et al.  Testing high-level QM/MM methods for modeling enzyme reactions: acetyl-CoA deprotonation in citrate synthase. , 2010, The journal of physical chemistry. B.

[21]  Lars Ridder,et al.  A Quantum Mechanical/Molecular Mechanical Study of the Hydroxylation of Phenol and Halogenated Derivatives by Phenol Hydroxylase , 2000 .

[22]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[23]  Adrian J Mulholland,et al.  Effects of Dispersion in Density Functional Based Quantum Mechanical/Molecular Mechanical Calculations on Cytochrome P450 Catalyzed Reactions. , 2012, Journal of chemical theory and computation.

[24]  Angela K. Wilson,et al.  Benchmark calculations with correlated molecular wave functions XII. Core correlation effects on the homonuclear diatomic molecules B2-F2 , 1997 .

[25]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[26]  Adrian J Mulholland,et al.  Analysis of chorismate mutase catalysis by QM/MM modelling of enzyme-catalysed and uncatalysed reactions. , 2011, Organic & biomolecular chemistry.

[27]  M. W. van der Kamp,et al.  Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. , 2013, Biochemistry.

[28]  M Karplus,et al.  Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparison , 1988, Proteins.

[29]  Adrian J Mulholland,et al.  High-accuracy computation of reaction barriers in enzymes. , 2006, Angewandte Chemie.

[30]  Johannes C. Hermann,et al.  Insights into enzyme catalysis from QM/MM modelling: transition state stabilization in chorismate mutase , 2003 .

[31]  S. Grimme,et al.  Is spin-component scaled second-order Møller-Plesset perturbation theory an appropriate method for the study of noncovalent interactions in molecules? , 2007, The journal of physical chemistry. A.

[32]  M. Karplus,et al.  Simulation of activation free energies in molecular systems , 1996 .

[33]  L. Hedstrom,et al.  A twisted base? The role of arginine in enzyme-catalyzed proton abstractions. , 2005, Archives of biochemistry and biophysics.

[34]  S. Grimme.  Improved second-order Møller–Plesset perturbation theory by separate scaling of parallel- and antiparallel-spin pair correlation energies , 2003 .

[35]  J. Harvey,et al.  Spin-forbidden CO ligand recombination in myoglobin. , 2004, Faraday discussions.