Enzyme polarization of substrates of dihydrofolate reductase by different theoretical methods

We have investigated the importance of polarization by the enzyme dihydrofolate reductase (DHFR) on its substrates, folate and dihydrofolate, using a series of quantum mechanical (QM) techniques (Hartree‐Fock (HF), Møller‐Plesset second‐order perturbation theory (MP2), local density approximation (LDA) and generalized gradient approximation (GGA) density functional theory (DFT) calculations) in which the bulk enzyme is included in the calculations as point charges. Polarization, in terms of both charges on components (residues) of the folate and dihydrofolate molecules and changes in the electron density, particularly of the pterin ring of the substrates, and the implications for the catalytic reduction are discussed. Significant differences in polarization behavior are observed for the different theoretical methods employed. The consequences of this, particularly for choosing an appropriate model for quantum mechanical/molecular mechanical (QM/MM) calculations, are pointed out. The HF and MP2 QM methods show small polarizations (∼0.04 electrons) of the pterin ring but quite large polarizations with both LDA and GGA DFT methods (0.3–0.5 electrons). This large difference in polarization for both folate and dihydrofolate arises as a result of substantial differences between the charge distributions for the gas‐phase DFT and HF calculations, specifically the charges on the dianionic glutamate side chain. Some recent literature reports of incorrect representation of anionic systems by DFT methods are noted. The DFT results are similar to the previously reported LDA DFT results of Bajorath et al. predicting a large polarization of the pterin ring of folate (Proteins 9:217–224,1991) and dihydrofolate (PNAS 88:6423–6426, 1991) of ∼0.5–0.6 electrons. Proteins 1999;37:157–165. ©1999 Wiley‐Liss, Inc.

[1]  J. Kraut,et al.  The electrostatic potential of Escherichia coli dihydrofolate reductase , 1991, Proteins.

[2]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[3]  S. Benkovic,et al.  Computational studies on pterins and speculations on the mechanism of action of dihydrofolate reductase. , 1989, Biochemical and biophysical research communications.

[4]  S. Benkovic,et al.  Consideration of the pH-dependent inhibition of dihydrofolate reductase by methotrexate. , 1997, Journal of molecular biology.

[5]  G C Roberts,et al.  Domain motions in dihydrofolate reductase: a molecular dynamics study. , 1997, Journal of molecular biology.

[6]  K. Morokuma,et al.  ONIOM: A Multilayered Integrated MO + MM Method for Geometry Optimizations and Single Point Energy Predictions. A Test for Diels−Alder Reactions and Pt(P(t-Bu)3)2 + H2 Oxidative Addition , 1996 .

[7]  J. Gready,et al.  Mechanistic aspects of biological redox reactions involving NADH 2: A combined semiempirical and ab initio study of hydride‐ion transfer between the NADH analogue, 1‐methyl‐dihydronicotinamide, and folate and dihydrofolate analogue substrates of dihydrofolate reductase , 1990 .

[8]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

[9]  J. Gready,et al.  Ionization state and pKa of pterin-analogue ligands bound to dihydrofolate reductase. , 1994, European journal of biochemistry.

[10]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[11]  P. Kraulis.  A program to produce both detailed and schematic plots of protein structures , 1991 .

[12]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[13]  Barbara J. Garrison,et al.  ELECTROSTATIC CHARACTERIZATION OF ENZYME COMPLEXES : EVALUATION OF THE MECHANISM OF CATALYSIS OF DIHYDROFOLATE REDUCTASE , 1997 .

[14]  Application of density functional theory to calculation of in-crystal anionic polarizability , 1999 .

[15]  Raymond L. Blakley,et al.  Eukaryotic dihydrofolate reductase. , 2006, Advances in enzymology and related areas of molecular biology.

[16]  J. Gready,et al.  Enzymic properties of a new mechanism-based substrate for dihydrofolate reductase , 1989 .

[17]  J. Gready,et al.  Molecular dynamics and free energy perturbation study of hydride‐ion transfer step in dihydrofolate reductase using combined quantum and molecular mechanical model , 1998 .

[18]  J Kraut,et al.  Theoretical studies on the dihydrofolate reductase mechanism: electronic polarization of bound substrates. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[19]  J. Kraut,et al.  Loop and subdomain movements in the mechanism of Escherichia coli dihydrofolate reductase: crystallographic evidence. , 1997, Biochemistry.

[20]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[21]  P. Kollman,et al.  A Simulation of the Catalytic Mechanism of Aspartylglucosaminidase Using ab Initio Quantum Mechanics and Molecular Dynamics , 1997 .

[22]  J. Kraut,et al.  Exploring the molecular mechanism of dihydrofolate reductase. , 1992, Faraday discussions.

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

[24]  Jürgen Bajorath,et al.  Electron redistribution on binding of a substrate to an enzyme: Folate and dihydrofolate reductase , 1991, Proteins.

[25]  David R. Lowis,et al.  Pterin 1H–3H tautomerism and its possible relevance to the binding of folate to dihydrofolate reductase , 1993 .

[26]  J. Kraut,et al.  Changes in the electron density of the cofactor NADPH on binding to E. coli dihydrofolate reductase , 1991, Proteins.

[27]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[28]  Bernard R. Brooks,et al.  ENZYME MECHANISMS WITH HYBRID QUANTUM AND MOLECULAR MECHANICAL POTENTIALS.I. THEORETICAL CONSIDERATIONS , 1996 .

[29]  G. Wagner,et al.  Detection of long-lived bound water molecules in complexes of human dihydrofolate reductase with methotrexate and NADPH. , 1995, Journal of molecular biology.

[30]  E A Merritt,et al.  Raster3D: photorealistic molecular graphics. , 1997, Methods in enzymology.

[31]  J. Gready.  Theoretical studies on the activation of the pterin cofactor in the catalytic mechanism of dihydrofolate reductase. , 1985, Biochemistry.

[32]  S. Benkovic,et al.  Insights into enzyme function from studies on mutants of dihydrofolate reductase. , 1988, Science.

[33]  John C. Slater,et al.  The self-consistent field for molecules and solids , 1974 .

[34]  Wang,et al.  Accurate and simple analytic representation of the electron-gas correlation energy. , 1992, Physical review. B, Condensed matter.

[35]  Gregory S. Tschumper,et al.  PREDICTING ELECTRON AFFINITIES WITH DENSITY FUNCTIONAL THEORY: SOME POSITIVE RESULTS FOR NEGATIVE IONS , 1997 .

[36]  Krishnan Raghavachari,et al.  Assessment of Gaussian-2 and density functional theories for the computation of ionization potentials and electron affinities , 1998 .

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