51:Journal of Molecular Evolution, 65:485-95 (2007)

Coevolutionary patterns in cytochrome c oxidase subunit I depend on structure and functional context

Wang ZO and Pollock DD

The strength and pattern of coevolution between amino acid residues varies depending on their structural and functional environment. This context dependence, along with differences in analytical technique, is responsible for different results among coevolutionary analyses of different proteins. It is thus important to perform detailed study of individual proteins to gain better insight into how context dependence can affect coevolutionary patterns even within individual proteins, and to unravel the details of context dependence with respect to structure and function. Here, we extend our previous study by presenting further analysis of residue coevolution in cytochrome c oxidase subunit I sequences from 231 vertebrates using a statistically robust phylogeny-based maximum likelihood ratio method. As in previous studies, a strong overall coevolutionary signal was detected, and coevolution within structural regions was significantly related to the Ca distances between residues. While the strong selection for adjacent residues among predicted coevolving pairs in the surface region indicates that the statistical method is highly selective for biologically relevant interactions, the coevolutionary signal was strongest in the transmembrane region, although the distances between coevolving residues were greater. This indicates that coevolution may act to maintain more global structural and functional constraints in the transmembrane region. In the transmembrane region, sites that coevolved according to polarity and hydrophobicity rather than volume had a greater tendency to co-localize with just one of the predicted proton channels (channel H). Thus, the details of coevolution in cytochrome c oxidase subunit I depend greatly on domain structure and residue physicochemical characteristics, but proximity to function appears to play a critical role. We hypothesize that the association of coevolutionary sites with channel H was caused by adaptive coevolution, and is indicative of a more important functional role for this channel.

Figure 1 (ppt pdf)

Figure 2 (ppt pdf)

Figure 3 (ppt pdf)

Figure 4 (pdf)

Supplementary Data (pdf)

Supplementary Figure S1 (png)

Supplementary Figure S2A (gif movie)

Supplementary Figure S2B (gif movie)