A comparative analysis of clustering algorithms: O2 migration in truncated hemoglobin I from transition networks

Cazade, Pierre-André and Zheng, Wenwei and Prada-Gracia, Diego and Berezovska, Ganna and Rao, Francesco and Clementi, Cecilia and Meuwly, Markus. (2015) A comparative analysis of clustering algorithms: O2 migration in truncated hemoglobin I from transition networks. Journal of Chemical Physics, 142 (2). 025103.

PDF - Published Version

Official URL: http://edoc.unibas.ch/41802/

Downloads: Statistics Overview


The ligand migration network for O2-diffusion in truncated Hemoglobin N is analyzed based on three different clustering schemes. For coordinate-based clustering, the conventional k-means and the kinetics-based Markov Clustering (MCL) methods are employed, whereas the locally scaled diffusion map (LSDMap) method is a collective-variable-based approach. It is found that all three methods agree well in their geometrical definition of the most important docking site, and all experimentally known docking sites are recovered by all three methods. Also, for most of the states, their population coincides quite favourably, whereas the kinetics of and between the states differs. One of the major differences between k-means and MCL clustering on the one hand and LSDMap on the other is that the latter finds one large primary cluster containing the Xe1a, IS1, and ENT states. This is related to the fact that the motion within the state occurs on similar time scales, whereas structurally the state is found to be quite diverse. In agreement with previous explicit atomistic simulations, the Xe3 pocket is found to be a highly dynamical site which points to its potential role as a hub in the network. This is also highlighted in the fact that LSDMap cannot identify this state. First passage time distributions from MCL clusterings using a one- (ligand-position) and two-dimensional (ligand-position and protein-structure) descriptor suggest that ligand- and protein-motions are coupled. The benefits and drawbacks of the three methods are discussed in a comparative fashion and highlight that depending on the questions at hand the best-performing method for a particular data set may differ.
Faculties and Departments:05 Faculty of Science > Departement Chemie > Chemie > Physikalische Chemie (Meuwly)
UniBasel Contributors:Meuwly, Markus
Item Type:Article, refereed
Article Subtype:Research Article
Publisher:AIP Publishing
Note:Publication type according to Uni Basel Research Database: Journal article -- The final publication is available at AIP, see DOI link.
Identification Number:
edoc DOI:
Last Modified:07 Dec 2016 15:08
Deposited On:25 Apr 2016 11:35

Repository Staff Only: item control page