Gapped BLAST and PSI-BLAST: a new generation of protein database search programs

Abstract

The BLAST programs are widely used tools for searching protein and DNA databases for sequence similarities. For protein comparisons, a variety of definitional, algorithmic and statistical refinements described here permits the execution time of the BLAST programs to be decreased substantially while enhancing their sensitivity to weak similarities. A new criterion for triggering the extension of word hits, combined with a new heuristic for generating gapped alignments, yields a gapped BLAST program that runs at approximately three times the speed of the original. In addition, a method is introduced for automatically combining statistically significant alignments produced by BLAST into a position-specific score matrix, and searching the database using this matrix. The resulting Position-Specific Iterated BLAST (PSI-BLAST) program runs at approximately the same speed per iteration as gapped BLAST, but in many cases is much more sensitive to weak but biologically relevant sequence similarities. PSI-BLAST is used to uncover several new and interesting members of the BRCT superfamily.

Cite this paper

Related papers

1. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice1994
2. STAR: ultrafast universal RNA-seq aligner2012
3. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability2013
4. MUSCLE: multiple sequence alignment with high accuracy and high throughput2004
5. Cutadapt removes adapter sequences from high-throughput sequencing reads2011