Bioinformatics Advance Access originally published online on February 26, 2008
Bioinformatics 2008 24(7):908-915; doi:10.1093/bioinformatics/btn057
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Tracing evolutionary pressure
1Division of Medicinal Chemistry, Leiden/Amsterdam Center for Drug Research, Leiden University, Einsteinweg 55, 2333CC, Leiden and 2CMBI, Radboud University Nijmegen, Toernooiveld 1, 6525 ED, Nijmegen, The Netherlands
*To whom correspondence should be addressed.
 |
Abstract |
|---|
Motivation: 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.
Availability: All computer programs and datasets used in this work are available at http://nava.liacs.nl/kye/TEA-O/ for academic use
Contact: k.ye{at}lacdr.leidenuniv.nl
Associate Editor: Dmitrij Frishman
Received on September 20, 2007; revised on February 6, 2008; accepted on February 9, 2008
CiteULike 
Connotea 
Del.icio.us What's this?