GeneRAGE: a robust algorithm for sequence clustering and domain detection

MOTIVATION Efficient, accurate and automatic clustering of large protein sequence datasets, such as complete proteomes, into families, according to sequence similarity. Detection and correction of false positive and negative relationships with subsequent detection and resolution of multi-domain proteins. RESULTS A new algorithm for the automatic clustering of protein sequence datasets has been developed. This algorithm represents all similarity relationships within the dataset in a binary matrix. Removal of false positives is achieved through subsequent symmetrification of the matrix using a Smith-Waterman dynamic programming alignment algorithm. Detection of multi-domain protein families and further false positive relationships within the symmetrical matrix is achieved through iterative processing of matrix elements with successive rounds of Smith-Waterman dynamic programming alignments. Recursive single-linkage clustering of the corrected matrix allows efficient and accurate family representation for each protein in the dataset. Initial clusters containing multi-domain families, are split into their constituent clusters using the information obtained by the multi-domain detection step. This algorithm can hence quickly and accurately cluster large protein datasets into families. Problems due to the presence of multi-domain proteins are minimized, allowing more precise clustering information to be obtained automatically. AVAILABILITY GeneRAGE (version 1.0) executable binaries for most platforms may be obtained from the authors on request. The system is available to academic users free of charge under license.

[1]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[2]  J. Lake,et al.  Genomic evidence for two functionally distinct gene classes. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[3]  D. Lipman,et al.  A genomic perspective on protein families. , 1997, Science.

[4]  C. Sander,et al.  Protein structure comparison by alignment of distance matrices. , 1993, Journal of molecular biology.

[5]  C R Woese,et al.  Methanococcus jannaschii genome: revisited. , 1996, Microbial & comparative genomics.

[6]  C. Chothia,et al.  Intermediate sequences increase the detection of homology between sequences. , 1997, Journal of molecular biology.

[7]  A. Bateman The structure of a domain common to archaebacteria and the homocystinuria disease protein. , 1997, Trends in biochemical sciences.

[8]  Christine A. Orengo,et al.  Protein Structure Comparison , 2002 .

[9]  E V Koonin,et al.  Evidence for a Family of Archaeal ATPases , 1997, Science.

[10]  Anton J. Enright,et al.  Protein interaction maps for complete genomes based on gene fusion events , 1999, Nature.

[11]  David C. Jones,et al.  GenTHREADER: an efficient and reliable protein fold recognition method for genomic sequences. , 1999, Journal of molecular biology.

[12]  Tim J. P. Hubbard,et al.  SCOP: a Structural Classification of Proteins database , 1999, Nucleic Acids Res..

[13]  D. Lipman,et al.  Improved tools for biological sequence comparison. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[14]  W. Pearson Effective protein sequence comparison. , 1996, Methods in enzymology.

[15]  C Sander,et al.  Mapping the Protein Universe , 1996, Science.

[16]  Richard Hughey,et al.  Hidden Markov models for detecting remote protein homologies , 1998, Bioinform..

[17]  Temple F. Smith,et al.  Comparison of biosequences , 1981 .

[18]  J. Coggins,et al.  The pentafunctional arom enzyme of Saccharomyces cerevisiae is a mosaic of monofunctional domains. , 1987, The Biochemical journal.

[19]  Thomas L. Madden,et al.  Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. , 1997, Nucleic acids research.

[20]  N. Kyrpides,et al.  Tetratrico-peptide-repeat proteins in the archaeon Methanococcus jannaschii. , 1998, Trends in biochemical sciences.