Skip to main content
Article thumbnail
Location of Repository

Convergence of Molecular Dynamics Simulations of Membrane Proteins

By Alan Grossfield, Scott E. Feller and Michael C. Pitman


ABSTRACT The central question in evaluating almost any result from a molecular dynamics simulation is whether the calculation has converged. Unfortunately, assessing the ergodicity of a single trajectory is very difficult to do. In this work, we assess the sampling of molecular dynamics simulations of the membrane protein rhodopsin by comparing the results from 26 independent trajectories, each run for 100 ns. By examining principal components and cluster populations, we show that rhodopsin’s fluctuations are not well described by 100 ns of dynamics, and that the sampling is not fully converged even for individual loops. The results serve as a reminder of the caution required when interpreting molecular dynamics simulations of macromolecules. Proteins 2007;67:31–40. VC 2007 Wiley-Liss, Inc. Key words: molecular dynamics; rhodopsin; ergodicity; principal component; cluster analysis; convergenc

Year: 2009
OAI identifier: oai:CiteSeerX.psu:
Provided by: CiteSeerX
Download PDF:
Sorry, we are unable to provide the full text but you may find it at the following location(s):
  • (external link)
  • http://membrane.urmc.rochester... (external link)
  • Suggested articles

    To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.