Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis

Copyright 2002 Wiley Periodicals, Inc.

Bibliographische Detailangaben
Veröffentlicht in:Journal of computational chemistry. - 1984. - 23(2002), 14 vom: 15. Nov., Seite 1314-22
1. Verfasser: Titmuss, Stephen J (VerfasserIn)
Weitere Verfasser: Cummins, Peter L, Rendell, Alistair P, Bliznyuk, Andrey A, Gready, Jill E
Format: Aufsatz
Sprache:English
Veröffentlicht: 2002
Zugriff auf das übergeordnete Werk:Journal of computational chemistry
Schlagworte:Comparative Study Journal Article Research Support, Non-U.S. Gov't Pterins 8-methylpterin 13005-87-3 Tetrahydrofolate Dehydrogenase EC 1.5.1.3
LEADER 01000naa a22002652 4500
001 NLM120729652
003 DE-627
005 20231222191518.0
007 tu
008 231222s2002 xx ||||| 00| ||eng c
028 5 2 |a pubmed24n0403.xml 
035 |a (DE-627)NLM120729652 
035 |a (NLM)12214314 
040 |a DE-627  |b ger  |c DE-627  |e rakwb 
041 |a eng 
100 1 |a Titmuss, Stephen J  |e verfasserin  |4 aut 
245 1 0 |a Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis 
264 1 |c 2002 
336 |a Text  |b txt  |2 rdacontent 
337 |a ohne Hilfsmittel zu benutzen  |b n  |2 rdamedia 
338 |a Band  |b nc  |2 rdacarrier 
500 |a Date Completed 20.03.2003 
500 |a Date Revised 24.11.2016 
500 |a published: Print 
500 |a Citation Status MEDLINE 
520 |a Copyright 2002 Wiley Periodicals, Inc. 
520 |a 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 
650 4 |a Comparative Study 
650 4 |a Journal Article 
650 4 |a Research Support, Non-U.S. Gov't 
650 7 |a Pterins  |2 NLM 
650 7 |a 8-methylpterin  |2 NLM 
650 7 |a 13005-87-3  |2 NLM 
650 7 |a Tetrahydrofolate Dehydrogenase  |2 NLM 
650 7 |a EC 1.5.1.3  |2 NLM 
700 1 |a Cummins, Peter L  |e verfasserin  |4 aut 
700 1 |a Rendell, Alistair P  |e verfasserin  |4 aut 
700 1 |a Bliznyuk, Andrey A  |e verfasserin  |4 aut 
700 1 |a Gready, Jill E  |e verfasserin  |4 aut 
773 0 8 |i Enthalten in  |t Journal of computational chemistry  |d 1984  |g 23(2002), 14 vom: 15. Nov., Seite 1314-22  |w (DE-627)NLM098138448  |x 1096-987X  |7 nnns 
773 1 8 |g volume:23  |g year:2002  |g number:14  |g day:15  |g month:11  |g pages:1314-22 
912 |a GBV_USEFLAG_A 
912 |a SYSFLAG_A 
912 |a GBV_NLM 
912 |a GBV_ILN_350 
951 |a AR 
952 |d 23  |j 2002  |e 14  |b 15  |c 11  |h 1314-22