Simple adjustment of the sequence weight algorithm remarkably enhances PSI-BLAST performance
Open Access
- 2 June 2017
- journal article
- software
- Published by Springer Nature in BMC Bioinformatics
- Vol. 18 (1) , 1-8
Abstract
PSI-BLAST, an extremely popular tool for sequence similarity search, features the utilization of Position-Specific Scoring Matrix (PSSM) constructed from a multiple sequence alignment (MSA). PSSM allows the detection of more distant homologs than a general amino acid substitution matrix does. An accurate estimation of the weights for sequences in an MSA is crucially important for PSSM construction. PSI-BLAST divides a given MSA into multiple blocks, for which sequence weights are calculated. When the block width becomes very narrow, the sequence weight calculation can be odd. We demonstrate that PSI-BLAST indeed generates a significant fraction of blocks having width less than 5, thereby degrading the PSI-BLAST performance. We revised the code of PSI-BLAST to prevent the blocks from being narrower than a given minimum block width (MBW). We designate the modified application of PSI-BLAST as PSI-BLASTexB. When MBW is 25, PSI-BLASTexB notably outperforms PSI-BLAST consistently for three independent benchmark sets. The performance boost is even more drastic when an MSA, instead of a sequence, is used as a query. Our results demonstrate that the generation of narrow-width blocks during the sequence weight calculation is a critically important factor that restricts the PSI-BLAST search performance. By preventing narrow blocks, PSI-BLASTexB upgrades the PSI-BLAST performance remarkably. Binaries and source codes of PSI-BLASTexB (MBW = 25) are available at https://github.com/kyungtaekLIM/PSI-BLASTexB .Keywords
All Related Versions
Funding Information
- AMED
This publication has 24 references indexed in Scilit:
- UniProt Protein Knowledgebase2017
- The Pfam protein families database: towards a more sustainable futureNucleic Acids Research, 2015
- SCOPe: Structural Classification of Proteins—extended, integrating SCOP and ASTRAL data and classification of new structuresNucleic Acids Research, 2013
- Revisiting amino acid substitution matrices for identifying distantly related proteinsBioinformatics, 2013
- MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and UsabilityMolecular Biology and Evolution, 2013
- HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignmentNature Methods, 2011
- Code Optimization of the Subroutine to Remove Near Identical Matches in the Sequence Database Homology Search Tool PSI-BLASTJournal of Computational Biology, 2010
- Sequence context-specific profiles for homology searchingProceedings of the National Academy of Sciences, 2009
- PSI-BLAST pseudocounts and the minimum description length principleNucleic Acids Research, 2008
- Amino acid substitution matrices from protein blocks.Proceedings of the National Academy of Sciences, 1992