Is an intuitive convergence definition of molecular dynamics simulations solely based on the root mean square deviation possible?

J Comput Biol. 2011 Aug;18(8):997-1005. doi: 10.1089/cmb.2010.0237. Epub 2011 Jun 24.

Abstract

Molecular dynamics is a commonly used technique in computational biology. One key issue of each molecular dynamics simulation is: When does this simulation reach equilibrium state? A widely used way to determine this is the visual and intuitive inspection of root mean square deviation (RMSD) plots of the simulation. Although this technique has been criticized several times, it is still often used. Therefore, we present a study proving that this method is not reliable at all. We conducted a survey with participants from the field in which we illustrated different RMSD plots to scientists in the field of molecular dynamics. These plots were randomized and repeated, using a statistical model and different variants of the plots. We show that there is no mutual consent about the point of equilibrium. The decisions are severely biased by different parameters. Therefore, we conclude that scientists should not discuss the equilibration of a molecular dynamics simulation on the basis of a RMSD plot.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms
  • Bias
  • Computational Biology* / methods
  • Computational Biology* / statistics & numerical data
  • Data Collection
  • Humans
  • Models, Statistical*
  • Molecular Dynamics Simulation / standards*
  • Molecular Dynamics Simulation / statistics & numerical data
  • Observer Variation