Fast structure similarity searches among protein models: efficient clustering of protein fragments
Department of Biomedical Sciences and Technologies, University of Udine Piazzale Kolbe 4, 33100 Udine - Italy and Istituto Nazionale Biostrutture e Biosistemi Viale Medaglie d'Oro 305, 00136 Roma, Italy
Algorithms for Molecular Biology 2012, 7:16 doi:10.1186/1748-7188-7-16Published: 29 May 2012
For many predictive applications a large number of models is generated and later clustered in subsets based on structure similarity. In most clustering algorithms an all-vs-all root mean square deviation (RMSD) comparison is performed. Most of the time is typically spent on comparison of non-similar structures. For sets with more than, say, 10,000 models this procedure is very time-consuming and alternative faster algorithms, restricting comparisons only to most similar structures would be useful.
We exploit the inverse triangle inequality on the RMSD between two structures given the RMSDs with a third structure. The lower bound on RMSD may be used, when restricting the search of similarity to a reasonably low RMSD threshold value, to speed up similarity searches significantly. Tests are performed on large sets of decoys which are widely used as test cases for predictive methods, with a speed-up of up to 100 times with respect to all-vs-all comparison depending on the set and parameters used. Sample applications are shown.
The algorithm presented here allows fast comparison of large data sets of structures with limited memory requirements. As an example of application we present clustering of more than 100000 fragments of length 5 from the top500H dataset into few hundred representative fragments. A more realistic scenario is provided by the search of similarity within the very large decoy sets used for the tests. Other applications regard filtering nearly-indentical conformation in selected CASP9 datasets and clustering molecular dynamics snapshots.
A linux executable and a Perl script with examples are given in the supplementary material (Additional file 1). The source code is available upon request from the authors.