MDSCAN: RMSD-based HDBSCAN clustering of long molecular dynamics
Autor
González-Alemán, Roy
Platero-Rochart, Daniel
Rodríguez-Serradet, Alejandro
Hernández-Rodríguez, Erix W.
Caballero, Julio
Leclerc, Fabrice
Montero-Cabrera, Luis
Fecha
2022Resumen
Motivation
The term clustering designates a comprehensive family of unsupervised learning methods allowing to group similar elements into sets called clusters. Geometrical clustering of molecular dynamics (MD) trajectories is a well-established analysis to gain insights into the conformational behavior of simulated systems. However, popular variants collapse when processing relatively long trajectories because of their quadratic memory or time complexity. From the arsenal of clustering algorithms, HDBSCAN stands out as a hierarchical density-based alternative that provides robust differentiation of intimately related elements from noise data. Although a very efficient implementation of this algorithm is available for programming-skilled users (HDBSCAN*), it cannot treat long trajectories under the de facto molecular similarity metric RMSD.
Results
Here, we propose MDSCAN, an HDBSCAN-inspired software specifically conceived for non-programmers users to perform memory-efficient RMSD-based clustering of long MD trajectories. Methodological improvements over the original version include the encoding of trajectories as a particular class of vantage-point tree (decreasing time complexity), and a dual-heap approach to construct a quasi-minimum spanning tree (reducing memory complexity). MDSCAN was able to process a trajectory of 1 million frames using the RMSD metric in about 21 h with <8 GB of RAM, a task that would have taken a similar time but more than 32 TB of RAM with the accelerated HDBSCAN* implementation generally used.
Fuente
Bioinformatics, 38(23), 5191-5198Link de Acceso
Click aquí para ver el documentoIdentificador DOI
doi.org/10.1093/bioinformatics/btac666Colecciones
La publicación tiene asociados los siguientes ficheros de licencia: