title: More than 1,001 problems with protein domain databases: transmembrane regions, signal peptides and the issue of sequence homology.
abstract: Large-scale genome sequencing gained general importance for life science because functional annotation of otherwise experimentally uncharacterized sequences is made possible by the theory of biomolecular sequence homology. Historically, the paradigm of similarity of protein sequences implying common structure, function and ancestry was generalized based on studies of globular domains. Having the same fold imposes strict conditions over the packing in the hydrophobic core requiring similarity of hydrophobic patterns. The implications of sequence similarity among non-globular protein segments have not been studied to the same extent; nevertheless, homology considerations are silently extended for them. This appears especially detrimental in the case of transmembrane helices (TMs) and signal peptides (SPs) where sequence similarity is necessarily a consequence of physical requirements rather than common ancestry. Thus, matching of SPs/TMs creates the illusion of matching hydrophobic cores. Therefore, inclusion of SPs/TMs into domain models can give rise to wrong annotations. More than 1001 domains among the 10,340 models of Pfam release 23 and 18 domains of SMART version 6 (out of 809) contain SP/TM regions. As expected, fragment-mode HMM searches generate promiscuous hits limited to solely the SP/TM part among clearly unrelated proteins. More worryingly, we show explicit examples that the scores of clearly false-positive hits, even in global-mode searches, can be elevated into the significance range just by matching the hydrophobic runs. In the PIR iProClass database v3.74 using conservative criteria, we find that at least between 2.1% and 13.6% of its annotated Pfam hits appear unjustified for a set of validated domain models. Thus, false-positive domain hits enforced by SP/TM regions can lead to dramatic annotation errors where the hit has nothing in common with the problematic domain model except the SP/TM region itself. We suggest a workflow of flagging problematic hits arising from SP/TM-containing models for critical reconsideration by annotation users.
authors.affiliation: Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), Singapore. wongwc@bii-sg.org
authors.forename: Wing-Cheong
authors.initials: WC
authors.lastname: Wong
authors.forename: Sebastian
authors.initials: S
authors.lastname: Maurer-Stroh
authors.forename: Frank
authors.initials: F
authors.lastname: Eisenhaber
mesh_terms: D000818:Animals; D019295:Computational Biology / Q000379:methods*; D030562:Databases, Protein; D006801:Humans; D008390:Markov Chains; D008565:Membrane Proteins / Q000737:chemistry / Q000145:classification; D010363:Pattern Recognition, Automated; D017510:Protein Folding; D021382:Protein Sorting Signals; D017434:Protein Structure, Tertiary; D011506:Proteins / Q000737:chemistry* / Q000145:classification; D015203:Reproducibility of Results; D017386:Sequence Homology, Amino Acid
chemical_list: D008565:Membrane Proteins; D021382:Protein Sorting Signals; D011506:Proteins