Detecting coevolving amino acid sites using Bayesian mutational mapping
1Department of Biological Statistics and Computational Biology, Cornell University Ithaca, NY 13101, USA
2Center for Bioinformatics, University of Copenhagen Copenhagen 2100 Kbh Ø, Denmark
*To whom correspondence should be addressed at Divergence, Inc., St. Louis, MO 63141, USA.
Motivation: The evolution of protein sequences is constrained by complex interactions between amino acid residues. Because harmful substitutions may be compensated for by other substitutions at neighboring sites, residues can coevolve. We describe a Bayesian phylogenetic approach to the detection of coevolving residues in protein families. This method, Bayesian mutational mapping (BMM), assigns mutations to the branches of the evolutionary tree stochastically, and then test statistics are calculated to determine whether a coevolutionary signal exists in the mapping. Posterior predictive P-values provide an estimate of significance, and specificity is maintained by integrating over uncertainty in the estimation of the tree topology, branch lengths and substitution rates. A coevolutionary Markov model for codon substitution is also described, and this model is used as the basis of several test statistics.
Results: Results on simulated coevolutionary data indicate that the BMM method can successfully detect nearly all coevolving sites when the model has been correctly specified, and that non-parametric statistics such as mutual information are generally less powerful than parametric statistics. On a dataset of eukaryotic proteins from the phosphoglycerate kinase (PGK) family, interdomain site contacts yield a significantly greater coevolutionary signal than interdomain non-contacts, an indication that the method provides information about interacting sites. Failure to account for the heterogeneity in rates across sites in PGK resulted in a less discriminating test, yielding a marked increase in the number of reported positives at both contact and non-contact sites.
Contact: matt{at}dimmic.net
Supplementary information: http://www.dimmic.net/supplement/
Received on January 15, 2005; accepted on March 27, 2005
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  M. Anisimova and C. Kosiol Investigating Protein-Coding Sequence Evolution with Probabilistic Codon Substitution Models Mol. Biol. Evol., February 1, 2009; 26(2): 255 - 271. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 V. N Minin and M. A Suchard Fast, accurate and simulation-free stochastic mapping Phil Trans R Soc B, December 27, 2008; 363(1512): 3985 - 3995. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 K. Y. Yip, P. Patel, P. M. Kim, D. M. Engelman, D. McDermott, and M. Gerstein An integrated system for studying residue coevolution in proteins Bioinformatics, January 15, 2008; 24(2): 290 - 292. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 D. S. Horner, W. Pirovano, and G. Pesole Correlated substitution analysis and the prediction of amino acid structural contacts Brief Bioinform, January 1, 2008; 9(1): 46 - 56. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. A. Fares and S. A. A. Travers A Novel Method for Detecting Intramolecular Coevolution: Adding a Further Dimension to Selective Constraints Analyses Genetics, May 1, 2006; 173(1): 9 - 23. [Abstract] [Full Text] [PDF]
|
|