Nucleic Acids Research Advance Access published online on December 16, 2008
Nucleic Acids Research, doi:10.1093/nar/gkn981
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Computational Biology |
PSI-BLAST pseudocounts and the minimum description length principle
National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health Bethesda, MD 20894
*To whom correspondence should be addressed. Tel: +1 301 435 7803; Fax: +1 301 480 2288; Email: altschul{at}ncbi.nlm.nih.gov
Received October 15, 2008. Revised November 18, 2008. Accepted November 20, 2008.
Position specific score matrices (PSSMs) are derived from multiple sequence alignments to aid in the recognition of distant protein sequence relationships. The PSI-BLAST protein database search program derives the column scores of its PSSMs with the aid of pseudocounts, added to the observed amino acid counts in a multiple alignment column. In the absence of theory, the number of pseudocounts used has been a completely empirical parameter. This article argues that the minimum description length principle can motivate the choice of this parameter. Specifically, for realistic alignments, the principle supports the practice of using a number of pseudocounts essentially independent of alignment size. However, it also implies that more highly conserved columns should use fewer pseudocounts, increasing the inter-column contrast of the implied PSSMs. A new method for calculating pseudocounts that significantly improves PSI-BLAST's; retrieval accuracy is now employed by default.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
J. Mrazek Finding sequence motifs in prokaryotic genomes--a brief practical guide for a microbiologist Brief Bioinform, September 1, 2009; 10(5): 525 - 536. [Abstract] [Full Text] [PDF] |
||||
