Correction for phylogeny, small number of observations and data redundancy improves the identification of coevolving amino acid pairs using mutual information

被引:66
作者
Marino Buslje, Cristina [1 ,2 ]
Santos, Javier [1 ,2 ]
Maria Delfino, Jose [1 ,2 ]
Nielsen, Morten [3 ]
机构
[1] Univ Buenos Aires, Dept Biol Chem, RA-1113 Buenos Aires, DF, Argentina
[2] Univ Buenos Aires, Inst Biochem & Biophys, IQUIFIB, Sch Pharm & Biochem, RA-1113 Buenos Aires, DF, Argentina
[3] Tech Univ Denmark, Ctr Biol Sequence Anal, Dept Syst Biol, DK-2800 Lyngby, Denmark
基金
美国国家卫生研究院;
关键词
PROTEIN; COEVOLUTION; PREDICTION; RESIDUES;
D O I
10.1093/bioinformatics/btp135
中图分类号
Q5 [生物化学];
学科分类号
071010 ; 081704 ;
摘要
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 functionally important in a given fold or protein family. Accurate identification of coevolving positions in protein sequences is difficult due to the high background signal imposed by phylogeny and noise. Several 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 redundancy and low-count corrections to account for small number of observations in limited size sequence families, can significantly improve the predictability of MI. The evaluation is made on large sets of both in silico-generated alignments as well as on biological sequence 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 sequence-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 sequences <62% identical 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.
引用
收藏
页码:1125 / 1131
页数:7
相关论文
共 24 条
[1]   Gapped BLAST and PSI-BLAST: a new generation of protein database search programs [J].
Altschul, SF ;
Madden, TL ;
Schaffer, AA ;
Zhang, JH ;
Zhang, Z ;
Miller, W ;
Lipman, DJ .
NUCLEIC ACIDS RESEARCH, 1997, 25 (17) :3389-3402
[2]  
[Anonymous], 1991, ELEMENTS INFORM THEO
[3]  
CHIU DKY, 1991, COMPUT APPL BIOSCI, V7, P347
[4]   Missense meanderings in sequence space: A biophysical view of protein evolution [J].
DePristo, MA ;
Weinreich, DM ;
Hartl, DL .
NATURE REVIEWS GENETICS, 2005, 6 (09) :678-687
[5]   Mutual information without the influence of phylogeny or entropy dramatically improves residue contact prediction [J].
Dunn, S. D. ;
Wahl, L. M. ;
Gloor, G. B. .
BIOINFORMATICS, 2008, 24 (03) :333-340
[6]   A novel method for detecting intramolecular coevolution: Adding a further dimension to selective constraints analyses [J].
Fares, Mario A. ;
Travers, Simon A. A. .
GENETICS, 2006, 173 (01) :9-23
[7]   Influence of conservation on calculations of amino acid covariance in multiple sequence alignments [J].
Fodor, AA ;
Aldrich, RW .
PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS, 2004, 56 (02) :211-221
[8]   Mutual information in protein multiple sequence alignments reveals two classes of coevolving positions [J].
Gloor, GB ;
Martin, LC ;
Wahl, LM ;
Dunn, SD .
BIOCHEMISTRY, 2005, 44 (19) :7156-7165
[9]   Finding coevolving amino acid residues using row and column weighting of mutual information and multi-dimensional amino acid representation [J].
Gouveia-Oliveira, Rodrigo ;
Pedersen, Anders G. .
ALGORITHMS FOR MOLECULAR BIOLOGY, 2007, 2 (1)
[10]  
HOBOHM U, 1992, PROTEIN SCI, V1, P409