Homology-based method for identification of protein repeats using statistical significance estimates.
Short protein repeats, frequently with a length between 20 and 40 residues, represent a significant fraction of known proteins. Many repeats appear to possess high amino acid substitution rates and thus recognition of repeat homologues is highly problematic. Even if the presence of a certain repeat family is known, the exact locations and the number of repetitive units often cannot be determined using current methods. We have devised an iterative algorithm based on optimal and sub-optimal score distributions from profile analysis that estimates the significance of all repeats that are detected in a single sequence. This procedure allows the identification of homologues at alignment scores lower than the highest optimal alignment score for non-homologous sequences. The method has been used to investigate the occurrence of eleven families of repeats in Saccharomyces cerevisiae, Caenorhabditis elegans and Homo sapiens accounting for 1055, 2205 and 2320 repeats, respectively. For these examples, the method is both more sensitive and more selective than conventional homology search procedures. The method allowed the detection in the SwissProt database of more than 2000 previously unrecognised repeats belonging to the 11 families. In addition, the method was used to merge several repeat families that previously were supposed to be distinct, indicating common phylogenetic origins for these families.
Andrade
M A
MA
European Molecular Biology Laboratory, Meyerhofstr. 1, Heidelberg, 69012, Germany.
Ponting
C P
CP
Gibson
T J
TJ
Bork
P
P
D000465:Algorithms; D019883:Alkyl and Aryl Transferases / Q000737:chemistry; D020816:Amino Acid Motifs; D000818:Animals; D017487:Ankyrins / Q000737:chemistry; D017173:Caenorhabditis elegans / Q000737:chemistry; D029742:Caenorhabditis elegans Proteins; D002352:Carrier Proteins / Q000737:chemistry; D018797:Cell Cycle Proteins; D019295:Computational Biology / Q000379:methods* / Q000706:statistics & numerical data*; D004268:DNA-Binding Proteins / Q000737:chemistry; D016208:Databases, Factual; D051231:Farnesyltranstransferase; D016678:Genome; D020662:Guanine Nucleotide Exchange Factors; D050956:HSP40 Heat-Shock Proteins; D006360:Heat-Shock Proteins; D015801:Helminth Proteins / Q000737:chemistry; D006801:Humans; D007930:Leucine / Q000737:chemistry; D008958:Models, Molecular; D018832:Molecular Chaperones; D009687:Nuclear Proteins; D010802:Phylogeny; D017433:Protein Structure, Secondary; D011506:Proteins / Q000737:chemistry* / Q000378:metabolism*; D020449:Repetitive Sequences, Amino Acid; D012441:Saccharomyces cerevisiae / Q000737:chemistry; D012680:Sensitivity and Specificity; D016415:Sequence Alignment; D017386:Sequence Homology, Amino Acid
D017487:Ankyrins; D029742:Caenorhabditis elegans Proteins; D002352:Carrier Proteins; D018797:Cell Cycle Proteins; D004268:DNA-Binding Proteins; C485782:DNAJC7 protein, human; D020662:Guanine Nucleotide Exchange Factors; D050956:HSP40 Heat-Shock Proteins; D006360:Heat-Shock Proteins; D015801:Helminth Proteins; D018832:Molecular Chaperones; D009687:Nuclear Proteins; D011506:Proteins; C062521:RCC1 protein, human; C400683:che-2 protein, C elegans; C499303:ran-3 protein, C elegans; D019883:Alkyl and Aryl Transferases; D051231:Farnesyltranstransferase; D007930:Leucine