Coevolution analyses illuminate the dependencies between amino acid sites in the chaperonin system GroES-L

Abstract

BackgroundGroESL is a heat-shock protein ubiquitous in bacteria and eukaryotic organelles. This evolutionarily conserved protein is involved in the folding of a wide variety of other proteins in the cytosol, being essential to the cell. The folding activity proceeds through strong conformational changes mediated by the co-chaperonin GroES and ATP. Functions alternative to folding have been previously described for GroEL in different bacterial groups, supporting enormous functional and structural plasticity for this molecule and the existence of a hidden combinatorial code in the protein sequence enabling such functions. Describing this plasticity can shed light on the functional diversity of GroEL. We hypothesize that different overlapping sets of amino acids coevolve within GroEL, GroES and between both these proteins. Shifts in these coevolutionary relationships may inevitably lead to evolution of alternative functions.ResultsWe conducted the first coevolution analyses in an extensive bacterial phylogeny, revealing complex networks of evolutionary dependencies between residues in GroESL. These networks differed among bacterial groups and involved amino acid sites with functional importance and others with previously unsuspected functional potential. Coevolutionary networks formed statistically independent units among bacterial groups and map to structurally continuous regions in the protein, suggesting their functional link. Sites involved in coevolution fell within narrow structural regions, supporting dynamic combinatorial functional links involving similar protein domains. Moreover, coevolving sites within a bacterial group mapped to regions previously identified as involved in folding-unrelated functions, and thus, coevolution may mediate alternative functions.ConclusionsOur results highlight the evolutionary plasticity of GroEL across the entire bacterial phylogeny. Evidence on the functional importance of coevolving sites illuminates the as yet unappreciated functional diversity of proteins.

The Phenotypic Plasticity of Duplicated Genes inSaccharomyces cerevisiaeand the Origin of Adaptations

“…We aligned the 600 nucleotide sequence regions upstream of duplicated genes and their orthologs, as these are likely to include most if not all the regulatory elements of the genes (Ohler and Niemann 2001). We then measured the coefficient of conservation (CC) for each nucleotide site using the entropy equation (Cover and Thomas 2006; Halabi et al 2009; Ruiz-Gonzalez and Fares 2013):$CC=f_k^{\left(a\right)}\ln \frac{f_k^{\left(a\right)}}{q^{\left(a\right)}}+\left(1-f_k^{\left(a\right)}\right)\ln \frac{1-f_k^{\left(a\right)}}{1-q^{\left(a\right)}};\forall a\in \left[A,T,G,C\right]$In this equation, CC of a nucleotide ( a ) at position ( k ) in an alignment is defined as the entropy of the observed frequency of a at k ($f_k^{\left(a\right)}$) relative to the background frequency of a in all sequences of the alignment ( q (a) ). Therefore, the more conserved the site the higher is its CC value.…”

A comprehensive analysis of the Omp85/TpsB protein superfamily structural diversity, taxonomic occurrence, and evolution

“…The diversity observed in the Omp85 family could reflect adaptations to different substrate (“client”) proteins, as has been observed in molecular chaperone protein families. Detailed studies on molecular chaperones found in the cytoplasm show high levels of variation with respect to their copy numbers; in order to cope with the assembly of their evolving range of substrate proteins, as well as to acquire novel (sub)functions themselves ( Henderson et al, 2013 ; Ruiz-Gonzalez and Fares, 2013 ).…”

The evolution of protein moonlighting: adaptive traps and promiscuity in the chaperonins

“…In a recent study, the statistical independence of such coevolving groups was examined, finding that when each group of bacteria was analysed independently, the regions affected by signatures of co-evolution changed substantially between bacterial groups [55]. In particular, many of the amino acids identified in bacteria specific groups of co-evolution mapped to regions known to be responsible for folding-unrelated Cpn60 functions.…”