Prediction of amino acid positions specific for functional groups in a protein family based on local sequence similarity.

The exchange of single amino acid residue in protein can substantially affect the specificity of molecular recognition. Many protein families can be divided into the groups based on specificity to recognized ligands. Prediction of group-discriminating residues within the certain family is extremely necessary for theoretical studies, enzyme engineering, drug design, and so on. The most existing methods use the multiple sequence alignment. They have the limitations in prediction accuracy due to the family sequence divergence and ligand-based grouping. We developed a new method SPrOS (Specificity Projection On Sequence) for estimating the specificity of residues to user-defined groups. SPrOS compares the sequence segments from the test protein and training proteins. Contrary to other segment-comparison approaches extracting the string motifs, SPrOS calculates the scores for single positions by the similarity of their surroundings. The method was evaluated on the simulated sequences and real protein families. The high-prediction accuracy was achieved for simulated sequences, in which SPrOS detected specific positions not predicted with the alignment-based method. For bacterial transcription factors (LacI/GalR) clearly divided into functional groups, the predicted specific residues corresponded to the published experimental data. In a more complicated case of protein kinases classified by inhibitor specificity, the positions predicted with high significance were located in ligand-binding areas. As the ligand specificity is not necessary coincided with phylogeny, evolutionary-coupled mutations could disturb the detection of ligand-specific residues. Excluding proximate homologs of the test protein kinase from the training set, we improved the prediction of the ligand-specific residues. The SPrOS is available at Copyright (c) 2015 John Wiley & Sons, Ltd.