Knapp B, Frantal S, Cibena M, Schreiner W, Bauer P
Department for Biosimulation and Bioinformatics, Medical University of Vienna, Vienna, Austria.
J Comput Biol. 2011 Aug;18(8):997-1005. doi: 10.1089/cmb.2010.0237. Epub 2011 Jun 24.
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.
分子动力学是计算生物学中常用的技术。每次分子动力学模拟的一个关键问题是:该模拟何时达到平衡状态?一种广泛使用的确定方法是对模拟的均方根偏差(RMSD)图进行直观的视觉检查。尽管该技术已多次受到批评,但仍经常被使用。因此,我们进行了一项研究,证明该方法根本不可靠。我们对该领域的参与者进行了一项调查,向分子动力学领域的科学家展示了不同的RMSD图。这些图是随机的且重复展示,使用了统计模型和不同的图变体。我们表明,对于平衡点没有达成一致意见。决策受到不同参数的严重影响。因此,我们得出结论,科学家不应基于RMSD图来讨论分子动力学模拟的平衡情况。