Will molecular dynamics simulations of proteins ever reach equilibrium?
Genheden, Samuel LU and Ryde, Ulf LU
(2012)
In Physical chemistry chemical physics : PCCP
14(24).
p.8662-8677
- Abstract
- We show that conformational entropies calculated for five proteins and protein-ligand complexes with dihedral-distribution histogramming, the von Mises approach, or quasi-harmonic analysis do not converge to any useful precision even if molecular dynamics (MD) simulations of 380-500 ns length are employed (the uncertainty is 12-89 kJ mol(-1)). To explain this, we suggest a simple protein model involving dihedrals with effective barriers forming a uniform distribution and show that for such a model, the entropy increases logarithmically with time until all significantly populated dihedral states have been sampled, in agreement with the simulations (during the simulations, 52-70% of the available dihedral phase space has been visited). This... (More)
- We show that conformational entropies calculated for five proteins and protein-ligand complexes with dihedral-distribution histogramming, the von Mises approach, or quasi-harmonic analysis do not converge to any useful precision even if molecular dynamics (MD) simulations of 380-500 ns length are employed (the uncertainty is 12-89 kJ mol(-1)). To explain this, we suggest a simple protein model involving dihedrals with effective barriers forming a uniform distribution and show that for such a model, the entropy increases logarithmically with time until all significantly populated dihedral states have been sampled, in agreement with the simulations (during the simulations, 52-70% of the available dihedral phase space has been visited). This is also confirmed by the analysis of the trajectories of a 1 ms simulation of bovine pancreatic trypsin inhibitor (31 kJ mol(-1) difference in the entropy between the first and second part of the simulation). Strictly speaking, this means that it is practically impossible to equilibrate MD simulations of proteins. We discuss the implications of such a lack of strict equilibration of protein MD simulations and show that ligand-binding free energies estimated with the MM/GBSA method (molecular mechanics with generalised Born and surface-area solvation) vary by 3-15 kJ mol(-1) during a 500 ns simulation (the higher estimate is caused by rare conformational changes), although they involve a questionable but well-converged normal-mode entropy estimate, whereas free energies estimated by free-energy perturbation vary by less than 0.6 kJ mol(-1) for the same simulation. (Less)
Please use this url to cite or link to this publication:
https://lup.lub.lu.se/record/2608652
- author
- Genheden, Samuel
LU
and Ryde, Ulf
LU
- organization
- publishing date
- 2012
- type
- Contribution to journal
- publication status
- published
- subject
- in
- Physical chemistry chemical physics : PCCP
- volume
- 14
- issue
- 24
- pages
- 8662 - 8677
- publisher
- Royal Society of Chemistry
- external identifiers
-
- wos:000304605600024
- pmid:22614001
- scopus:84864261290
- pmid:22614001
- ISSN
- 1463-9084
- DOI
- 10.1039/c2cp23961b
- 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
- f09f3e75-c891-414f-acb8-22fd63e6aa94 (old id 2608652)
- date added to LUP
- 2016-04-01 14:48:05
- date last changed
- 2023-02-22 06:38:41
@article{f09f3e75-c891-414f-acb8-22fd63e6aa94,
abstract = {{We show that conformational entropies calculated for five proteins and protein-ligand complexes with dihedral-distribution histogramming, the von Mises approach, or quasi-harmonic analysis do not converge to any useful precision even if molecular dynamics (MD) simulations of 380-500 ns length are employed (the uncertainty is 12-89 kJ mol(-1)). To explain this, we suggest a simple protein model involving dihedrals with effective barriers forming a uniform distribution and show that for such a model, the entropy increases logarithmically with time until all significantly populated dihedral states have been sampled, in agreement with the simulations (during the simulations, 52-70% of the available dihedral phase space has been visited). This is also confirmed by the analysis of the trajectories of a 1 ms simulation of bovine pancreatic trypsin inhibitor (31 kJ mol(-1) difference in the entropy between the first and second part of the simulation). Strictly speaking, this means that it is practically impossible to equilibrate MD simulations of proteins. We discuss the implications of such a lack of strict equilibration of protein MD simulations and show that ligand-binding free energies estimated with the MM/GBSA method (molecular mechanics with generalised Born and surface-area solvation) vary by 3-15 kJ mol(-1) during a 500 ns simulation (the higher estimate is caused by rare conformational changes), although they involve a questionable but well-converged normal-mode entropy estimate, whereas free energies estimated by free-energy perturbation vary by less than 0.6 kJ mol(-1) for the same simulation.}},
author = {{Genheden, Samuel and Ryde, Ulf}},
issn = {{1463-9084}},
language = {{eng}},
number = {{24}},
pages = {{8662--8677}},
publisher = {{Royal Society of Chemistry}},
series = {{Physical chemistry chemical physics : PCCP}},
title = {{Will molecular dynamics simulations of proteins ever reach equilibrium?}},
url = {{https://lup.lub.lu.se/search/files/4175419/3412384.pdf}},
doi = {{10.1039/c2cp23961b}},
volume = {{14}},
year = {{2012}},
}