corner
corner

Phys. Rev. E 65, 031910 (2002) [10 pages]

Convergence of sampling in protein simulations

Download: PDF (410 kB) Buy this article Export: BibTeX or EndNote (RIS)

Berk Hess
Department of Biophysical Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands

Received 3 August 2001; published 1 March 2002

With molecular dynamics protein dynamics can be simulated in atomic detail. Current computers are not fast enough to probe all available conformations, but fluctuations around one conformation can be sampled to a reasonable extent. The motions with the largest fluctuations can be filtered out of a simulation using covariance or principal component analysis. A problem with this analysis is that random diffusion can appear as correlated motion. An analysis is presented of how long a simulation should be to obtain relevant results for global motions. The analysis reveals that the cosine content of the principal components is a good indicator for bad sampling.

© 2002 The American Physical Society

URL:
http://link.aps.org/doi/10.1103/PhysRevE.65.031910
DOI:
10.1103/PhysRevE.65.031910
PACS:
87.14.Ee, 87.15.Aa, 87.15.He