title: Maximising the size of non-redundant protein datasets using graph theory.
abstract: Analysis of protein data sets often requires prior removal of redundancy, so that data is not biased by containing similar proteins. This is usually achieved by pairwise comparison of sequences, followed by purging so that no two pairs have similarities above a chosen threshold. From a starting set, such as the PDB or a genome, one should remove as few sequences as possible, to give the largest possible non-redundant set for subsequent analysis. Protein redundancy can be represented as a graph, with proteins as nodes connected by undirected edges, if they have a pairwise similarity above the chosen threshold. The problem is then equivalent to finding the maximum independent set (MIS), where as few nodes are removed as possible to remove all edges. We tested seven MIS algorithms, three of which are new. We applied the methods to the PDB, subsets of the PDB, various genomes and the BHOLSIB benchmark datasets. For PDB subsets of up to 1000 proteins, we could compare to the exact MIS, found by the Cliquer algorithm. The best algorithm was the new method, Leaf. This works by adding clique members that have no edges to nodes outside the clique to the MIS, starting with the smallest cliques. For PDB subsets of up to 1000 members, it usually finds the MIS and is fast enough to apply to data sets of tens of thousands of proteins. Leaf gives sets that are around 10% larger than the commonly used PISCES algorithm, that are of identical quality. We therefore suggest that Leaf should be the method of choice for generating non-redundant protein data sets, though it is ineffective on dense graphs, such as the BHOLSIB benchmarks. The Leaf algorithm is available at: https://github.com/SimonCB765/Leaf, and sets from genomes and the PDB are available at: http://www.bioinf.manchester.ac.uk/leaf/.
authors: Manchester Institute of Biotechnology, Faculty of Life Sciences, The University of Manchester, Manchester, United Kingdom.
authors: Simon C
authors: SC
authors: Bull
authors: Mark R
authors: MR
authors: Muldoon
authors: Andrew J
authors: AJ
authors: Doig
mesh_terms: D000465:Algorithms; D030562:Databases, Protein; D008958:Models, Molecular; D011506:Proteins; D016415:Sequence Alignment
chemical_list: D011506:Proteins
references: Bioinformatics. 2003 Aug 12;19(12):1589-91
references: 12912846
references: J Mol Biol. 1990 Oct 5;215(3):403-10
references: 2231712
references: Protein Sci. 1994 Mar;3(3):522-4
references: 8019422
references: Nucleic Acids Res. 1997 Sep 1;25(17):3389-402
references: 9254694