Principal component and clustering analysis on molecular dynamics data of the ribosomal L11·23S subdomain

Abstract

With improvements in computer speed and algorithm efficiency, MD simulations are sampling larger amounts of molecular and biomolecular conformations. Being able to qualitatively and quantitatively sift these conformations into meaningful groups is a difficult and important task, especially when considering the structure-activity paradigm. Here we present a study that combines two popular techniques, principal component (PC) analysis and clustering, for revealing major conformational changes that occur in molecular dynamics (MD) simulations. Specifically, we explored how clustering different PC subspaces effects the resulting clusters versus clustering the complete trajectory data. As a case example, we used the trajectory data from an explicitly solvated simulation of a bacteria’s L11·23S ribosomal subdomain, which is a target of thiopeptide antibiotics. Clustering was performed, using K-means and average-linkage algorithms, on data involving the first two to the first five PC subspace dimensions. For the average-linkage algorithm we found that data-point membership, cluster shape, and cluster size depended on the selected PC subspace data. In contrast, K-means provided very consistent results regardless of the selected subspace. Since we present results on a single model system, generalization concerning the clustering of different PC subspaces of other molecular systems is currently premature. However, our hope is that this study illustrates a) the complexities in selecting the appropriate clustering algorithm, b) the complexities in interpreting and validating their results, and c) by combining PC analysis with subsequent clustering valuable dynamic and conformational information can be obtained.Electronic supplementary materialThe online version of this article (doi:10.1007/s00894-012-1563-4) contains supplementary material, which is available to authorized users.

Dynamics of Bcl-xL in Water and Membrane: Molecular Simulations

“…All the trajectories were oriented with respect to the uncomplexed Bcl-xl that was averaged over the last 10ns of the simulations in water, with the tail region excluded (owing to the large flexibilities). In the first set, more than 80% of the variance was covered by the first three PCs (see Figure S12), which is in agreement with previous reports on other systems [49,50]. Distributions of the conformations projected along the first three PCs in Figure 10 shows separate clusters originating from the complexed and uncomplexed systems in water and in membrane; the complex in the membrane is well separated from the other systems (Movie S7).…”

Use of restrained molecular dynamics to predict the conformations of phosphorylated receiver domains in two‐component signaling systems

“…To reduce the dimensionality of the data, we performed PCA on the combined trajectories for each rec domain (data not shown). Performing PCA prior to clustering can save significantly on computational cost and time, as well as filter out high‐frequency conformational variance (noise) from the data …”

Solving the riddle: Unraveling the mechanisms of blocking the binding of leukotoxin by therapeutic antagonists in periodontal diseases

“…Principal component analysis (PCA) is a method that is used to measure the changes in motions, dynamics, and conformation of proteins. 68,69 This technique helps to quantify the direction and magnitude (amplitude) of protein motions, which are described as eigenvectors and eigenvalues, respectively. 68,70 In other words, the eigenvector describes the directions of motion with a corresponding eigenvalue which represents the energetic contribution of a particular component to the protein motion.…”