Exploring a phylogenetic approach for the detection of correlated substitutions in proteins

The remarkable conservation of protein structure, compared with that of sequences, suggests that in the course of evolution, residue substitutions which tend to destabilize a particular structure must be compensated by other substitutions that confer greater stability on that structure. Several approaches have been designed to detect correlated changes in a set of homologous sequences. However, most of them do not take into account the phylogeny of the sequences, and it has been shown that their detection power is weak. It remains unclear whether coevolution could be a general process at the level of amino acids of proteins. In the present study, we analyze the phylogenetic reconstruction of 15 sets of homologous proteins to assess, under different conditions, whether a significant amount of coevolving sites can be detected. Two criteria are used to detect significantly cosubstituting sites. One criterion corresponds to that of Shindyalov, Kolchanov, and Sander. The second one is based on intensive simulations of evolution of protein sequences along a phylogeny to estimate the significance of the number of observed cosubstitutions for pairs of sites. Our results show an important sensitivity of the detection of cosubstituting sites to the model used for the phylogenetic reconstruction. Not considering the uncertainty associated with the reconstructed data might lead to detecting numerous false-positive pairs of sites. Finally, significant amounts of coevolving pairs could be found only when substitutions affecting the physicochemical properties of the amino acids were considered. Such results suggest evidence of a cosubstitution mechanism in protein evolution. However, the identification of nonambiguous coevolving sites is still unresolved.