Effect of geometry optimisations on QM-cluster and QM/MM studies of reaction energies in proteins
Sumner, Sophie ; Söderhjelm, Pär LU and Ryde, Ulf LU
(2013)
In Journal of Chemical Theory and Computation
9(9).
p.4205-4214
- Abstract
- We have examined the effect of geometry optimisation on energies calculated with the quantum-mechanical (QM) cluster, the combined QM and molecular-mechanics (QM/MM), the big-QM approaches (very large single-point QM calculations taken from QM/MM-optimised structures, including all atoms within 4.5 Å of the minimal active site, all buried charged groups in the protein, and truncations moved at least three residues away from the active site). We study a simple proton-transfer reaction between His-79 and Cys-546 in the active site of [Ni,Fe] hydrogenase and optimise QM systems of 50 different sizes (56–362 atoms). Geometries optimised with QM/MM are stable and reliable, whereas QM-cluster optimisations give larger changes in the structures... (More)
- We have examined the effect of geometry optimisation on energies calculated with the quantum-mechanical (QM) cluster, the combined QM and molecular-mechanics (QM/MM), the big-QM approaches (very large single-point QM calculations taken from QM/MM-optimised structures, including all atoms within 4.5 Å of the minimal active site, all buried charged groups in the protein, and truncations moved at least three residues away from the active site). We study a simple proton-transfer reaction between His-79 and Cys-546 in the active site of [Ni,Fe] hydrogenase and optimise QM systems of 50 different sizes (56–362 atoms). Geometries optimised with QM/MM are stable and reliable, whereas QM-cluster optimisations give larger changes in the structures and sometimes lead to large distortions in the active site if some hydrogen-bond partners to the metal ligands are omitted. Keeping 2–3 atoms for each truncated residue (rather than one) fixed in the optimisation improves the results, but does not solve all problems for the QM-cluster optimisations. QM-cluster energies in vacuum and a continuum solvent are insensitive to the geometry optimisations with a mean absolute change upon the optimisations of only 4–7 kJ/mol. This shows that geometry optimisations do not decrease the dependence of QM-cluster energies on how the QM system is selected – there is still a ~60 kJ/mol difference between calculations in which groups have been added to the QM system according to their distance to the active site or based on QM/MM free-energy components. QM/MM energies do not show such a difference, but they converge rather slowly with respect to the size of the QM system, although the convergence is improved by moving truncations away from the active site. The big-QM energies are stable over the 50 different optimised structures, 57±1 kJ/mol, although some smaller trends can be discerned. This shows that both QM-cluster geometries and energies should be interpreted with caution. Instead, we recommend QM/MM for geometry optimisations and energies calculated by the big-QM approach. (Less)
Please use this url to cite or link to this publication:
https://lup.lub.lu.se/record/4226429
- author
- Sumner, Sophie
; Söderhjelm, Pär
LU
and Ryde, Ulf
LU
- organization
- publishing date
- 2013
- type
- Contribution to journal
- publication status
- published
- subject
- keywords
- Quantum mechanical cluster calculations, QM/MM, geometry optimisation, density-functional theory, [Ni, Fe] hydrogenase.
- in
- Journal of Chemical Theory and Computation
- volume
- 9
- issue
- 9
- pages
- 4205 - 4214
- publisher
- The American Chemical Society (ACS)
- external identifiers
-
- wos:000330096800035
- scopus:84884185644
- pmid:26592409
- ISSN
- 1549-9618
- DOI
- 10.1021/ct400339c
- language
- English
- LU publication?
- yes
- additional info
- The information about affiliations in this record was updated in December 2015. The record was previously connected to the following departments: Theoretical Chemistry (S) (011001039)
- id
- 72f549fa-57fb-4919-b8e5-347e32b61807 (old id 4226429)
- date added to LUP
- 2016-04-01 11:14:00
- date last changed
- 2023-02-20 17:30:29
@article{72f549fa-57fb-4919-b8e5-347e32b61807,
abstract = {{We have examined the effect of geometry optimisation on energies calculated with the quantum-mechanical (QM) cluster, the combined QM and molecular-mechanics (QM/MM), the big-QM approaches (very large single-point QM calculations taken from QM/MM-optimised structures, including all atoms within 4.5 Å of the minimal active site, all buried charged groups in the protein, and truncations moved at least three residues away from the active site). We study a simple proton-transfer reaction between His-79 and Cys-546 in the active site of [Ni,Fe] hydrogenase and optimise QM systems of 50 different sizes (56–362 atoms). Geometries optimised with QM/MM are stable and reliable, whereas QM-cluster optimisations give larger changes in the structures and sometimes lead to large distortions in the active site if some hydrogen-bond partners to the metal ligands are omitted. Keeping 2–3 atoms for each truncated residue (rather than one) fixed in the optimisation improves the results, but does not solve all problems for the QM-cluster optimisations. QM-cluster energies in vacuum and a continuum solvent are insensitive to the geometry optimisations with a mean absolute change upon the optimisations of only 4–7 kJ/mol. This shows that geometry optimisations do not decrease the dependence of QM-cluster energies on how the QM system is selected – there is still a ~60 kJ/mol difference between calculations in which groups have been added to the QM system according to their distance to the active site or based on QM/MM free-energy components. QM/MM energies do not show such a difference, but they converge rather slowly with respect to the size of the QM system, although the convergence is improved by moving truncations away from the active site. The big-QM energies are stable over the 50 different optimised structures, 57±1 kJ/mol, although some smaller trends can be discerned. This shows that both QM-cluster geometries and energies should be interpreted with caution. Instead, we recommend QM/MM for geometry optimisations and energies calculated by the big-QM approach.}},
author = {{Sumner, Sophie and Söderhjelm, Pär and Ryde, Ulf}},
issn = {{1549-9618}},
keywords = {{Quantum mechanical cluster calculations; QM/MM; geometry optimisation; density-functional theory; [Ni; Fe] hydrogenase.}},
language = {{eng}},
number = {{9}},
pages = {{4205--4214}},
publisher = {{The American Chemical Society (ACS)}},
series = {{Journal of Chemical Theory and Computation}},
title = {{Effect of geometry optimisations on QM-cluster and QM/MM studies of reaction energies in proteins}},
url = {{https://lup.lub.lu.se/search/files/2491818/4226432.pdf}},
doi = {{10.1021/ct400339c}},
volume = {{9}},
year = {{2013}},
}