Efficient estimation of emission probabilities in profile hidden Markov models.
BACKGROUND
Profile hidden Markov models provide a sensitive method for performing sequence database search and aligning multiple sequences. One of the drawbacks of the hidden Markov model is that the conserved amino acids are not emphasized, but signal and noise are treated equally. For this reason, the number of estimated emission parameters is often enormous. Focusing the analysis on conserved residues only should increase the accuracy of sequence database search.


RESULTS
We address this issue with a new method for efficient emission probability (EEP) estimation, in which amino acids are divided into effective and ineffective residues at each conserved alignment position. A practical study with 20 protein families demonstrated that the EEP method is capable of detecting family members from other proteins with sensitivity of 98% and specificity of 99% on the average, even if the number of free emission parameters was decreased to 15% of the original. In the database search for TIM barrel sequences, EEP recognizes the family members nearly as accurately as HMMER or Blast, but the number of false positive sequences was significantly less than that obtained with the other methods.


BACKGROUND
The algorithms written in C language are available on request from the authors.
Ahola
Virpi
V
Department of Statistics, FIN-20014 University of Turku, Finland. virahola@utu.fi
Aittokallio
Tero
T
Uusipaikka
Esa
E
Vihinen
Mauno
M
D000465:Algorithms; D000596:Amino Acids / Q000737:chemistry / Q000145:classification; D030562:Databases, Protein; D004798:Enzymes / Q000737:chemistry / Q000145:classification; D005189:False Positive Reactions; D008390:Markov Chains; D008957:Models, Genetic; D015233:Models, Statistical; D011506:Proteins / Q000737:chemistry* / Q000145:classification; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D016415:Sequence Alignment / Q000379:methods*; D020539:Sequence Analysis, Protein / Q000379:methods*
D000596:Amino Acids; D004798:Enzymes; D011506:Proteins