Three improvements to the blastp search of genome databases

The BLASTP program is a search tool for databases of protein sequences that is widely used by biolo-gists as a first step in investigating new genome se-quences. BLASTP finds high-scoring local alignments ) without gaps between a query sequence and sequences in the database. The score of an alignment is the sum of the scores of individual alignments between amino acids that make up the protein. These individual scores come from a scoring matrix modeling the rate of evolutionary mutation. Here we provide a detailed description of the original program and three separate optimisations to it. BLASTP consists of three steps, that we call neighbourhood con-struction, hit detection, and hit extension. The three opti-misations target hit extension since it accounts for 93 % of the execution time. The first optimisation alters the data representation of the query sequence and the related code for indexing the scoring matrix. The second optimisation performs extensions in step-sizes of two rather than one. The third optimisation forestalls the calling of the hit exten-sion step in cases that are unlikely to lead to a high-scoring alignment. Individually, the three optimisations show speed ups of 15%, 48%, and 63 % respectively. 1