Comparing geometric and kinetic cluster algorithms for molecular simulation data

Bettina Keller, Xavier Daura, Wilfred F. Van Gunsteren

Research output: Contribution to journalArticleResearchpeer-review

93 Citations (Scopus)


The identification of metastable states of a molecule plays an important role in the interpretation of molecular simulation data because the free-energy surface, the relative populations in this landscape, and ultimately also the dynamics of the molecule under study can be described in terms of these states. We compare the results of three different geometric cluster algorithms (neighbor algorithm, K-medoids algorithm, and common-nearest-neighbor algorithm) among each other and to the results of a kinetic cluster algorithm. First, we demonstrate the characteristics of each of the geometric cluster algorithms using five two-dimensional data sets. Second, we analyze the molecular dynamics data of a Β -heptapeptide in methanol-a molecule that exhibits a distinct folded state, a structurally diverse unfolded state, and a fast folding/unfolding equilibrium-using both geometric and kinetic cluster algorithms. We find that geometric clustering strongly depends on the algorithm used and that the density based common-nearest-neighbor algorithm is the most robust of the three geometric cluster algorithms with respect to variations in the input parameters and the distance metric. When comparing the geometric cluster results to the metastable states of the Β -heptapeptide as identified by kinetic clustering, we find that in most cases the folded state is identified correctly but the overlap of geometric clusters with further metastable states is often at best approximate. © 2010 American Institute of Physics.
Original languageEnglish
Article number074110
JournalJournal of Chemical Physics
Publication statusPublished - 26 Feb 2010


Dive into the research topics of 'Comparing geometric and kinetic cluster algorithms for molecular simulation data'. Together they form a unique fingerprint.

Cite this