Tracing evolutionary pressure.
BACKGROUND
Recent advances in sequencing techniques have yielded enormous amounts of protein sequence data from various species. This large dataset allows sequence comparison between paralogous and orthologous proteins to identify motifs or functional positions that account for the differences of functional subgroups ('specificity' positions). Algorithms such as SDPpred and the two-entropies analysis (TEA) have been developed to detect such specificity positions from a multiple sequence alignment (MSA) grouped into classes according to certain biological functions. Other algorithms such as TreeDet compute a classification and then predict specificity positions associated with it. However, there are still many unresolved questions: Was the optimal subdivision of a protein family achieved? Do the definitions at different levels of the phylogenetic tree affect the prediction of specificity positions? Can the whole phylogenetic tree be used instead of only one level in it to predict specificity positions?


RESULTS
Here we present a novel method, TEA-O (Two-entropies analysis-Objective), to trace the evolutionary pressure from the root to the branches of the phylogenetic tree. At each level of the tree, a TEA plot is produced to capture the signal of the evolutionary pressure. A consensus TEA-O plot is composed from the whole series of plots to provide a condensed representation. Positions related to functions that evolved early (conserved) or later (specificity) are close to the lower-left or upper-left corner of the TEA-O plot, respectively. This novel approach allows an unbiased, user-independent, analysis of residue relevance in a protein family. We compared our TEA-O method with various algorithms using both synthetic and real protein sequences. The results show that our method is robust, sensitive to subtle differences in evolutionary pressure during evolution and comprehensive because all positions in the MSA are presented in the consensus plot.


BACKGROUND
All computer programs and datasets used in this work are available at http://nava.liacs.nl/kye/TEA-O/ for academic use.
Ye
Kai
K
Division of Medicinal Chemistry, Leiden/Amsterdam Center for Drug Research, Leiden University, Einsteinweg 55, 2333CC, Leiden, The Netherlands. k.ye@lacdr.leidenuniv.nl
Vriend
Gert
G
IJzerman
Adriaan P
AP
D000465:Algorithms; D001483:Base Sequence; D005075:Biological Evolution; D003198:Computer Simulation; D004252:DNA Mutational Analysis / Q000379:methods*; D019143:Evolution, Molecular; D008957:Models, Genetic; D008969:Molecular Sequence Data; D010802:Phylogeny; D011506:Proteins / Q000235:genetics*; D012641:Selection, Genetic; D017422:Sequence Analysis, DNA / Q000379:methods*
D011506:Proteins