going beyond clustering in md trajectory analysis an application to villin headpiece folding超越在md轨迹聚类分析应用程序villin帽子折叠.pdfVIP

Going beyond Clustering in MD Trajectory Analysis: An Application to Villin Headpiece Folding 1 2¤ 1,2 Aruna Rajan , Peter L. Freddolino , Klaus Schulten * 1 Department of Physics, University of Illinois, Urbana-Champaign, Illinois, United States of America, 2 Beckman Institute for Advanced Science and Technology, University of Illinois, Urbana-Champaign, Illinois, United States of America Abstract Recent advances in computing technology have enabled microsecond long all-atom molecular dynamics (MD) simulations of biological systems. Methods that can distill the salient features of such large trajectories are now urgently needed. Conventional clustering methods used to analyze MD trajectories suffer from various setbacks, namely (i) they are not data driven, (ii) they are unstable to noise and changes in cut-off parameters such as cluster radius and cluster number, and (iii) they do not reduce the dimensionality of the trajectories, and hence are unsuitable for finding collective coordinates. We advocate the application of principal component analysis (PCA) and a non-metric multidimensional scaling (nMDS) method to reduce MD trajectories and overcome the drawbacks of clustering. To illustrate the superiority of nMDS over other methods in reducing data and reproducing salient features, we analyze three complete villin headpiece folding trajectories. Our analysis suggests that the folding process of the villin headpiece is structurally heterogeneous. Citation: Rajan A, Freddolino PL, Schulten K (2010) Going beyond Clustering in MD Trajectory Analysis: An Application to Villin Headpiece Folding. PLoS ONE 5(4): e9890. doi:10.1371/journal.pone.0009890 Editor: Ashley M. Buckle, Monash University, Australia Received December 16, 2009; Accept