Markovian domain fingerprinting: statistical segmentation of protein sequences.
BACKGROUND
Characterization of a protein family by its distinct sequence domains is crucial for functional annotation and correct classification of newly discovered proteins. Conventional Multiple Sequence Alignment (MSA) based methods find difficulties when faced with heterogeneous groups of proteins. However, even many families of proteins that do share a common domain contain instances of several other domains, without any common underlying linear ordering. Ignoring this modularity may lead to poor or even false classification results. An automated method that can analyze a group of proteins into the sequence domains it contains is therefore highly desirable.


RESULTS
We apply a novel method to the problem of protein domain detection. The method takes as input an unaligned group of protein sequences. It segments them and clusters the segments into groups sharing the same underlying statistics. A Variable Memory Markov (VMM) model is built using a Prediction Suffix Tree (PST) data structure for each group of segments. Refinement is achieved by letting the PSTs compete over the segments, and a deterministic annealing framework infers the number of underlying PST models while avoiding many inferior solutions. We show that regions of similar statistics correlate well with protein sequence domains, by matching a unique signature to each domain. This is done in a fully automated manner, and does not require or attempt an MSA. Several representative cases are analyzed. We identify a protein fusion event, refine an HMM superfamily classification into the underlying families the HMM cannot separate, and detect all 12 instances of a short domain in a group of 396 sequences.


BACKGROUND
jill@cs.huji.ac.il; tishby@cs.huji.ac.il.
Bejerano
G
G
School of Computer Science & Engineering, The Hebrew University, Jerusalem 91904, Israel.
Seldin
Y
Y
Margalit
H
H
Tishby
N
N
D000465:Algorithms; D019295:Computational Biology; D004250:DNA Topoisomerases, Type II / Q000235:genetics; D005982:Glutathione Transferase / Q000235:genetics; D018398:Homeodomain Proteins / Q000235:genetics; D008390:Markov Chains; D010449:Peptide Mapping / Q000706:statistics & numerical data*; D017434:Protein Structure, Tertiary; D011506:Proteins / Q000235:genetics*; D016415:Sequence Alignment / Q000706:statistics & numerical data; D014157:Transcription Factors / Q000235:genetics
D018398:Homeodomain Proteins; D011506:Proteins; D014157:Transcription Factors; D005982:Glutathione Transferase; D004250:DNA Topoisomerases, Type II