Systematic and fully automated identification of protein sequence patterns.
We present an efficient algorithm to systematically and automatically identify patterns in protein sequence families. The procedure is based on the Splash deterministic pattern discovery algorithm and on a framework to assess the statistical significance of patterns. We demonstrate its application to the fully automated discovery of patterns in 974 PROSITE families (the complete subset of PROSITE families which are defined by patterns and contain DR records). Splash generates patterns with better specificity and undiminished sensitivity, or vice versa, in 28% of the families; identical statistics were obtained in 48% of the families, worse statistics in 15%, and mixed behavior in the remaining 9%. In about 75% of the cases, Splash patterns identify sequence sites that overlap more than 50% with the corresponding PROSITE pattern. The procedure is sufficiently rapid to enable its use for daily curation of existing motif and profile databases. Third, our results show that the statistical significance of discovered patterns correlates well with their biological significance. The trypsin subfamily of serine proteases is used to illustrate this method's ability to exhaustively discover all motifs in a family that are statistically and biologically significant. Finally, we discuss applications of sequence patterns to multiple sequence alignment and the training of more sensitive score-based motif models, akin to the procedure used by PSI-BLAST. All results are available at httpl//www.research.ibm.com/spat/.
Hart
R K
RK
IBM Computational Biology Center, T.J. Watson Research Center, Yorktown Heights, NY 10598, USA. reece@in-machina.com
Royyuru
A K
AK
Stolovitzky
G
G
Califano
A
A
D000465:Algorithms; D000595:Amino Acid Sequence; D000818:Animals; D019295:Computational Biology; D016208:Databases, Factual; D008390:Markov Chains; D008958:Models, Molecular; D010363:Pattern Recognition, Automated; D011487:Protein Conformation; D011506:Proteins / Q000737:chemistry*; D012680:Sensitivity and Specificity; D020539:Sequence Analysis, Protein / Q000706:statistics & numerical data*; D012697:Serine Endopeptidases / Q000737:chemistry; D014357:Trypsin / Q000737:chemistry
D011506:Proteins; D012697:Serine Endopeptidases; D014357:Trypsin