Bioinformatics Advance Access published online on March 10, 2009
Bioinformatics, doi:10.1093/bioinformatics/btp135
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Correction for phylogeny, small number of observations and data redundancy improves the identification of coevolving amino acid pairs using mutual information.
1Department of Biological Chemistry and Institute of Biochemistry and Biophysics (IQUIFIB), School of Pharmacy and Biochemistry, University of Buenos Aires, Junín 956, 1113 Buenos Aires, Argentina.
2Centre for Biological Sequence Analysis, Department of Systems Biology. The Technical University of Denmark, Building 208, DK-2800 Lyngby. Denmark
*To whom correspondence should be addressed., E-mail: cmb{at}qb.ffyb.uba.ar; mniel{at}cbs.dtu.dk
| Abstract |
|---|
Motivation: Mutual information (MI) theory is often applied to predict positional correlations in a multiple sequence alignment (MSA) to make possible the analysis of those positions structurally or func-tionally important in a given fold or protein family. Accurate identifi-cation of coevolving positions in protein sequences is difficult due to the high background signal imposed by phylogeny and noise. Sev-eral methods have been proposed using MI to identify coevolving amino acids in protein families.
Results: After evaluating two current methods, we demonstrate how the use of sequence-weighting techniques to reduce sequence re-dundancy and low-count corrections to account for small number of observations in limited size sequence families, can significantly im-prove the predictability of MI. The evaluation is made on large sets of both in silico-generated alignments as well as on biological se-quence data. The methods included in the analysis are the APC (average product correction) and RCW (row-column weighting) methods. The best performing method was APC including se-quence-weighting and low-count corrections. The use of sequence permutations to calculate a MI rescaling is shown to significantly improve the prediction accuracy and allows for direct comparison of information values across protein families. Finally, we demonstrate how a lower bound of 400 unique sequences is needed in an MSA in order to achieve meaningful predictive performances. With our contribution, we achieve a noteworthy improvement on the current procedures to determine coevolution and residue contacts, and we believe that this will have potential impacts on the understanding of protein structure, function and folding.
Contact: cmb{at}qb.ffyb.uba.ar; mniel{at}cbs.dtu.dk
Associate Editor: Prof. Dmitrij Frishman
Received on October 14, 2008; revised on March 5, 2009; accepted on March 5, 2009