Bioinformatics Vol. 19 no. 18 2003
pages 2359-2368
© 2003 Oxford University Press
Efficient estimation of emission probabilities in profile hidden Markov models
1 Department of Statistics, FIN-20014 University of Turku, Finland, 2 Turku Centre for Computer Science, FIN-20520 Turku, Finland, 3 Institute of Medical Technology, FIN-33014 University of Tampere, Finland and 4 Research Unit, Tampere University Hospital, FIN-33520 Tampere, Finland
Received on November 19, 2002
; revised on March 17, 2003
; accepted on June 13, 2003
Motivation: 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.
Availability: The algorithms written in C language are available on request from the authors.
Contact: virahola{at}utu.fi