Constraining the electron densities in DFT method as an effective way for ab initio studies of metal‐catalyzed reactions

The use of hybrid ab initio QM/MM methods in studies of metalloenzymes and related systems presents a major challenge to computational chemists. Methods that include the metal ion in the quantum mechanical region should also include the ligands of the metal in this region. Such a treatment, however, should be very demanding if one is interested in performing the configurational averaging needed for proper calculations of activation free energies. In the present work we examine the ability of the frozen DFT (FDFT) and the constrained DFT (CDFT) approaches to be used in ab initio studies of metal‐catalyzed reactions, while allowing for an effective QM (rather than a QM/MM) treatment of the reacting complex. These approaches allow one to treat the entire enzyme by ab initio DFT methods, while confining the SCF calculations to a relatively small subsystem and keeping the electron density of the rest of the system frozen (or constrained). It is found that the FDFT and CDFT models can reproduce the trend obtained by a full DFT calculation of a proton transfer between two water molecules in a (Im)3Zn2+(H2O)2 system. This and related test cases indicate that our approximated models should be capable of providing a reliable representation of the energetics of metalloenzymes. The reasons for the special efficiency of the FDFT approach are clarified, and the strategies that can be used in FDFT studies of metalloenzymes are outlined. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1554–1561, 2000

[1]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[2]  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.

[3]  Arieh Warshel,et al.  An empirical valence bond approach for comparing reactions in solutions and in enzymes , 1980 .

[4]  B. Thole.  Molecular polarizabilities calculated with a modified dipole interaction , 1981 .

[5]  U. Singh,et al.  A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl− exchange reaction and gas phase protonation of polyethers , 1986 .

[6]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

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

[8]  William H. Miller,et al.  An empirical valence bond model for constructing global potential energy surfaces for chemical reactions of polyatomic molecular systems , 1990 .

[9]  Arieh Warshel,et al.  Free energy relationships in metalloenzyme-catalyzed reactions. Calculations of the effects of metal ion substitutions in staphylococcal nuclease , 1990 .

[10]  H. -. Kim,et al.  Equilibrium and nonequilibrium solvation and solute electronic structure , 1990 .

[11]  D. Salahub,et al.  New algorithm for the optimization of geometries in local density functional theory , 1990 .

[12]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

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

[14]  D. Salahub,et al.  Gaussian density functional calculations on the allyl and polyene radicals: C3H5 to C11H13 , 1991 .

[15]  Cortona,et al.  Self-consistently determined properties of solids without band-structure calculations. , 1991, Physical review. B, Condensed matter.

[16]  Yang,et al.  Direct calculation of electron density in density-functional theory. , 1991, Physical review letters.

[17]  H. J. Kim,et al.  A theoretical model for SN1 ionic dissociation in solution. 1. Activation free energetics and transition-state structure , 1992 .

[18]  A. Warshel,et al.  Computer simulation of the initial proton transfer step in human carbonic anhydrase I. , 1992, Journal of molecular biology.

[19]  W. Miller,et al.  Classical trajectory studies of the molecular dissociation dynamics of formaldehyde: H2CO→H2+CO , 1992 .

[20]  A. Warshel,et al.  Quantum-mechanical calculations of solvation free energies. A combined ab initio pseudopotential free-energy perturbation approach , 1992 .

[21]  Arieh Warshel,et al.  Simulation of enzyme reactions using valence bond force fields and other hybrid quantum/classical approaches , 1993 .

[22]  A. Warshel,et al.  Frozen density functional approach for ab initio calculations of solvated molecules , 1993 .

[23]  György G. Ferenczy,et al.  Quantum mechanical computations on very large molecular systems: The local self‐consistent field method , 1994, J. Comput. Chem..

[24]  G. Voth,et al.  A path integral study of electronic polarization and nonlinear coupling effects in condensed phase proton transfer reactions , 1994 .

[25]  A. Warshel,et al.  Why have mutagenesis studies not located the general base in ras p21 , 1994, Nature Structural Biology.

[26]  M. Thompson.  Hybrid quantum mechanical/molecular mechanical force field development for large flexible molecules. A molecular dynamics study of 18-crown-6 , 1995 .

[27]  Feliu Maseras,et al.  IMOMM: A new integrated ab initio + molecular mechanics geometry optimization scheme of equilibrium structures and transition states , 1995, J. Comput. Chem..

[28]  Weitao Yang,et al.  A density‐matrix divide‐and‐conquer approach for electronic structure calculations of large molecules , 1995 .

[29]  W. Thiel,et al.  Hybrid Models for Combined Quantum Mechanical and Molecular Mechanical Approaches , 1996 .

[30]  William N. Lipscomb,et al.  Recent Advances in Zinc Enzymology. , 1996, Chemical reviews.

[31]  M. Thompson,et al.  QM/MMpol: A Consistent Model for Solute/Solvent Polarization. Application to the Aqueous Solvation and Spectroscopy of Formaldehyde, Acetaldehyde, and Acetone , 1996 .

[32]  Edward I. Solomon,et al.  Structural and Functional Aspects of Metal Sites in Biology. , 1996, Chemical reviews.

[33]  Jacques Weber,et al.  Kohn-Sham equations with constrained electron density: an iterative evaluation of the ground-state electron density of interacting molecules , 1996 .

[34]  J. Mccammon,et al.  Density functional based parametrization of a valence bond method and its applications in quantum‐classical molecular dynamics simulations of enzymatic reactions , 1996 .

[35]  J. A. McCammon,et al.  Quantum-Classical Molecular Dynamics Simulations of Proton Transfer Processes in Molecular Complexes and in Enzymes , 1996 .

[36]  Arieh Warshel,et al.  Ab-initio Frozen Density Functional Calculations of Proton Transfer Reactions in Solution , 1996 .

[37]  Arieh Warshel,et al.  Computer Simulations of the Action of Metalloenzymes , 1997 .

[38]  C. Kölmel,et al.  Combining ab initio techniques with analytical potential functions for structure predictions of large systems: Method and application to crystalline silica polymorphs , 1997 .

[39]  E. Neria,et al.  Molecular dynamics of an enzyme reaction: proton transfer in TIM , 1997 .

[40]  Weitao Yang.  Absolute-energy-minimum principles for linear-scaling electronic-structure calculations , 1997 .

[41]  T. Wesołowski.  Density Functional Theory with approximate kinetic energy functionals applied to hydrogen bonds , 1997 .

[42]  Jiali Gao,et al.  HYBRID AB INITIO QM/MM SIMULATION OF N-METHYLACETAMIDE IN AQUEOUS SOLUTION , 1997 .

[43]  Daniel Borgis,et al.  An extended empirical valence bond model for describing proton transfer in H+(H2O)n clusters and liquid water , 1998 .

[44]  P. Kollman,et al.  Combined ab initio and Free Energy Calculations To Study Reactions in Enzymes and Solution: Amide Hydrolysis in Trypsin and Aqueous Solution , 1998 .

[45]  Gregory A. Voth,et al.  Multistate Empirical Valence Bond Model for Proton Transport in Water , 1998 .

[46]  Dongsheng Lu,et al.  PROTON TRANSFER IN THE ENZYME CARBONIC ANHYDRASE : AN AB INITIO STUDY , 1998 .

[47]  Arieh Warshel,et al.  Hybrid ab Initio Quantum Mechanics/Molecular Mechanics Calculations of Free Energy Surfaces for Enzymatic Reactions: The Nucleophilic Attack in Subtilisin , 1998 .

[48]  Weitao Yang,et al.  FROZEN DENSITY MATRIX APPROACH FOR ELECTRONIC STRUCTURE CALCULATIONS , 1998 .

[49]  Darrin M. York,et al.  Quantum Mechanical Treatment of Biological Macromolecules in Solution Using Linear-Scaling Electronic Structure Methods , 1998 .

[50]  Incorporating Protein Environments in Density Functional Theory: A Self-Consistent Reaction Field Calculation of Redox Potentials of [2Fe2S] Clusters in Ferredoxin and Phthalate Dioxygenase Reductase , 1998 .

[51]  K. Merz,et al.  Implementation and Testing of a Frozen Density Matrix−Divide and Conquer Algorithm , 1999 .

[52]  Donald G. Truhlar,et al.  Multiconfiguration molecular mechanics algorithm for potential energy surfaces of chemical reactions , 2000 .