Evolution of protein complexes by duplication of homomeric interactions

From duplication to fusion: Expanding Dayhoff's model of protein evolution

Dayhoff's hypothesis suggests that complex proteins emerged from simpler peptides or domains, which duplicated and fused to create more complex proteins and novel functions. These processes expanded and diversified the protein repertoire within organisms. Extensive studies and reviews over the past two decades have highlighted the impact of gene duplication on protein evolution. However, the role of fusion in this evolutionary narrative remains less understood. This perspective seeks to address this gap by emphasizing the role of fusion in evolution. Fusion is critical in determining the evolutionary fate of duplicated protomers, either preserving their ancestral function or evolving entirely new functions. It complements mutations, insertions, and deletions as evolutionary steps to enhance protein evolvability by expanding the capacity of the protein to explore new structural and functional space.

Symmetry facilitated the evolution of heterospecificity and high-order stoichiometry in vertebrate hemoglobin

Many proteins form paralogous multimers-molecular complexes in which evolutionarily related proteins are arranged into specific quaternary structures. Little is known about the mechanisms by which they acquired their stoichiometry (the number of total subunits in the complex) and heterospecificity (the preference of subunits for their paralogs rather than other copies of the same protein). Here, we use ancestral protein reconstruction and biochemical experiments to study historical increases in stoichiometry and specificity during the evolution of vertebrate hemoglobin (Hb), an α2β2 heterotetramer that evolved from a homodimeric ancestor after a gene duplication. We show that the mechanisms for this evolutionary transition were simple. One hydrophobic substitution in subunit β after the gene duplication was sufficient to cause the ancestral dimer to homotetramerize with high affinity across a new interface. During this same interval, a single-residue deletion in subunit α at the older interface conferred specificity for the heterotetrameric form and the trans-orientation of subunits within it. These sudden transitions in stoichiometry and specificity were possible because the interfaces in Hb are isologous, binding via the same surface patch on interacting subunits, but rotated 180° relative to each other. This architecture amplifies the impacts of individual mutations on stoichiometry and specificity, especially in higher-order complexes, and allows single substitutions to differentially affect heteromeric and homomeric interactions. Our findings suggest that elaborate and specific symmetrical molecular complexes may often evolve via simple genetic and physical mechanisms.

Symmetry facilitated the evolution of heterospecificity and high-order stoichiometry in vertebrate hemoglobin

Many proteins form paralogous multimers - molecular complexes in which evolutionarily related proteins are arranged into specific quaternary structures. Little is known about the mechanisms by which they acquired their stoichiometry (the number of total subunits in the complex) and heterospecificity (the preference of subunits for their paralogs rather than other copies of the same protein). Here we use ancestral protein reconstruction and biochemical experiments to study historical increases in stoichiometry and specificity during the evolution of vertebrate hemoglobin (Hb), aα 2 β 2 heterotetramer that evolved from a homodimeric ancestor after a gene duplication. We show that the mechanisms for this evolutionary transition were simple. One hydrophobic substitution in subunit β after the gene duplication was sufficient to cause the ancestral dimer to homotetramerize with high affinity across a new interface. During this same interval, a single-residue deletion in subunit α at the older interface conferred specificity for the heterotetrameric form and the trans-orientation of subunits within it. These sudden transitions in stoichiometry and specificity were possible because the interfaces in Hb are isologous - involving the same surface patch on interacting subunits, rotated 180° relative to each other. This architecture amplifies the impacts of individual mutations on stoichiometry and specificity, especially in higher-order complexes, and allows single substitutions to differentially affect heteromeric vs homomeric interactions. Our findings suggest that elaborate and specific symmetrical molecular complexes may often evolve via simple genetic and physical mechanisms.

Altered assembly paths mitigate interference among paralogous complexes

Protein complexes are fundamental to all cellular processes, so understanding their evolutionary history and assembly processes is important. Gene duplication followed by divergence is considered a primary mechanism for diversifying protein complexes. Nonetheless, to what extent assembly of present-day paralogous complexes has been constrained by their long evolutionary pathways and how cross-complex interference is avoided remain unanswered questions. Subunits of protein complexes are often stabilized upon complex formation, whereas unincorporated subunits are degraded. How such cooperative stability influences protein complex assembly also remains unclear. Here, we demonstrate that subcomplexes determined by cooperative stabilization interactions serve as building blocks for protein complex assembly. We further develop a protein stability-guided method to compare the assembly processes of paralogous complexes in cellulo. Our findings support that oligomeric state and the structural organization of paralogous complexes can be maintained even if their assembly processes are rearranged. Our results indicate that divergent assembly processes by paralogous complexes not only enable the complexes to evolve new functions, but also reinforce their segregation by establishing incompatibility against deleterious hybrid assemblies.

Compensatory mutations potentiate constructive neutral evolution by gene duplication

The functions of proteins generally depend on their assembly into complexes. During evolution, some complexes have transitioned from homomers encoded by a single gene to heteromers encoded by duplicate genes. This transition could occur without adaptive evolution through intermolecular compensatory mutations. Here, we experimentally duplicated and evolved a homodimeric enzyme to determine whether and how this could happen. We identified hundreds of deleterious mutations that inactivate individual homodimers but produce functional enzymes when coexpressed as duplicated proteins that heterodimerize. The structure of one such heteromer reveals how both losses of function are buffered through the introduction of asymmetry in the complex that allows them to subfunctionalize. Constructive neutral evolution can thus occur by gene duplication followed by only one deleterious mutation per duplicate.

Following the Evolutionary Paths of Dscam1 Proteins toward Highly Specific Homophilic Interactions

Many adhesion proteins, evolutionarily related through gene duplication, exhibit distinct and precise interaction preferences and affinities crucial for cell patterning. Yet, the evolutionary paths by which these proteins acquire new specificities and prevent cross-interactions within their family members remain unknown. To bridge this gap, this study focuses on Drosophila Down syndrome cell adhesion molecule-1 (Dscam1) proteins, which are cell adhesion proteins that have undergone extensive gene duplication. Dscam1 evolved under strong selective pressure to achieve strict homophilic recognition, essential for neuronal self-avoidance and patterning. Through a combination of phylogenetic analyses, ancestral sequence reconstruction, and cell aggregation assays, we studied the evolutionary trajectory of Dscam1 exon 4 across various insect lineages. We demonstrated that recent Dscam1 duplications in the mosquito lineage bind with strict homophilic specificities without any cross-interactions. We found that ancestral and intermediate Dscam1 isoforms maintained their homophilic binding capabilities, with some intermediate isoforms also engaging in promiscuous interactions with other paralogs. Our results highlight the robust selective pressure for homophilic specificity integral to the Dscam1 function within the process of neuronal self-avoidance. Importantly, our study suggests that the path to achieving such selective specificity does not introduce disruptive mutations that prevent self-binding but includes evolutionary intermediates that demonstrate promiscuous heterophilic interactions. Overall, these results offer insights into evolutionary strategies that underlie adhesion protein interaction specificities.

Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

An atlas of protein homo-oligomerization across domains of life

Protein structures are essential to understanding cellular processes in molecular detail. While advances in artificial intelligence revealed the tertiary structure of proteins at scale, their quaternary structure remains mostly unknown. We devise a scalable strategy based on AlphaFold2 to predict homo-oligomeric assemblies across four proteomes spanning the tree of life. Our results suggest that approximately 45% of an archaeal proteome and a bacterial proteome and 20% of two eukaryotic proteomes form homomers. Our predictions accurately capture protein homo-oligomerization, recapitulate megadalton complexes, and unveil hundreds of homo-oligomer types, including three confirmed experimentally by structure determination. Integrating these datasets with omics information suggests that a majority of known protein complexes are symmetric. Finally, these datasets provide a structural context for interpreting disease mutations and reveal coiled-coil regions as major enablers of quaternary structure evolution in human. Our strategy is applicable to any organism and provides a comprehensive view of homo-oligomerization in proteomes.

Compensatory mutations potentiate constructive neutral evolution by gene duplication

Protein functions generally depend on their assembly into complexes. During evolution, some complexes have transitioned from homomers encoded by a single gene to heteromers encoded by duplicate genes. This transition could occur without adaptive evolution through intermolecular compensatory mutations. Here, we experimentally duplicate and evolve an homodimeric enzyme to examine if and how this could happen. We identify hundreds of deleterious mutations that inactivate individual homodimers but produce functional enzymes when co-expressed as duplicated proteins that heterodimerize. The structure of one such heteromer reveals how both losses of function are buffered through the introduction of asymmetry in the complex that allows them to subfunctionalize. Constructive neutral evolution can thus occur by gene duplication followed by only one deleterious mutation per duplicate.

Newfound Coding Potential of Transcripts Unveils Missing Members of Human Protein Communities

Recent proteogenomic approaches have led to the discovery that regions of the transcriptome previously annotated as non-coding regions [i.e., untranslated regions (UTRs), open reading frames overlapping annotated coding sequences in a different reading frame, and non-coding RNAs] frequently encode proteins, termed alternative proteins (altProts). This suggests that previously identified protein-protein interaction (PPI) networks are partially incomplete because altProts are not present in conventional protein databases. Here, we used the proteogenomic resource OpenProt and a combined spectrum- and peptide-centric analysis for the re-analysis of a high-throughput human network proteomics dataset, thereby revealing the presence of 261 altProts in the network. We found 19 genes encoding both an annotated (reference) and an alternative protein interacting with each other. Of the 117 altProts encoded by pseudogenes, 38 are direct interactors of reference proteins encoded by their respective parental genes. Finally, we experimentally validate several interactions involving altProts. These data improve the blueprints of the human PPI network and suggest functional roles for hundreds of altProts.

Evolution of homology: From archetype towards a holistic concept of cell type

The concept of homology lies in the heart of comparative biological science. The distinction between homology as structure and analogy as function has shaped the evolutionary paradigm for a century and formed the axis of comparative anatomy and embryology, which accept the identity of structure as a ground measure of relatedness. The advent of single-cell genomics overturned the classical view of cell homology by establishing a backbone regulatory identity of cell types, the basic biological units bridging the molecular and phenotypic dimensions, to reveal that the cell is the most flexible unit of living matter and that many approaches of classical biology need to be revised to understand evolution and diversity at the cellular level. The emerging theory of cell types explicitly decouples cell identity from phenotype, essentially allowing for the divergence of evolutionarily related morphotypes beyond recognition, as well as it decouples ontogenetic cell lineage from cell-type phylogeny, whereby explicating that cell types can share common descent regardless of their structure, function or developmental origin. The article succinctly summarizes current progress and opinion in this field and formulates a more generalistic view of biological cell types as avatars, transient or terminal cell states deployed in a continuum of states by the developmental programme of one and the same omnipotent cell, capable of changing or combining identities with distinct evolutionary histories or inventing ad hoc identities that never existed in evolution or development. It highlights how the new logic grounded in the regulatory nature of cell identity transforms the concepts of cell homology and phenotypic stability, suggesting that cellular evolution is inherently and massively network-like, with one-to-one homologies being rather uncommon and restricted to shallower levels of the animal tree of life.

Paralog Explorer: A resource for mining information about paralogs in common research organisms

Paralogs are genes which arose via gene duplication, and when such paralogs retain overlapping or redundant function, this poses a challenge to functional genetics research. Recent technological advancements have made it possible to systematically probe gene function for redundant genes using dual or multiplex gene perturbation, and there is a need for a simple bioinformatic tool to identify putative paralogs of a gene(s) of interest. We have developed Paralog Explorer (https://www.flyrnai.org/tools/paralogs/), an online resource that allows researchers to quickly and accurately identify candidate paralogous genes in the genomes of the model organisms D. melanogaster, C. elegans, D. rerio, M. musculus, and H. sapiens. Paralog Explorer deploys an effective between-species ortholog prediction software, DIOPT, to analyze within-species paralogs. Paralog Explorer allows users to identify candidate paralogs, and to navigate relevant databases regarding gene co-expression, protein-protein and genetic interaction, as well as gene ontology and phenotype annotations. Altogether, this tool extends the value of current ortholog prediction resources by providing sophisticated features useful for identification and study of paralogous genes.

How gene duplication diversifies the landscape of protein oligomeric state and function

Oligomeric proteins are central to cellular life and the duplication and divergence of their genes is a key driver of evolutionary innovations. The duplication of a gene coding for an oligomeric protein has numerous possible outcomes, which motivates questions on the relationship between structural and functional divergence. How do protein oligomeric states diversify after gene duplication? In the simple case of duplication of a homo-oligomeric protein gene, what properties can influence the fate of descendant paralogs toward forming independent homomers or maintaining their interaction as a complex? Furthermore, how are functional innovations associated with the diversification of oligomeric states? Here, we review recent literature and present specific examples in an attempt to illustrate and answer these questions.

Oligomeric interactions maintain active-site structure in a noncooperative enzyme family

The evolutionary benefit accounting for widespread conservation of oligomeric structures in proteins lacking evidence of intersubunit cooperativity remains unclear. Here, crystal and cryo-EM structures, and enzymological data, demonstrate that a conserved tetramer interface maintains the active-site structure in one such class of proteins, the short-chain dehydrogenase/reductase (SDR) superfamily. Phylogenetic comparisons support a significantly longer polypeptide being required to maintain an equivalent active-site structure in the context of a single subunit. Oligomerization therefore enhances evolutionary fitness by reducing the metabolic cost of enzyme biosynthesis. The large surface area of the structure-stabilizing oligomeric interface yields a synergistic gain in fitness by increasing tolerance to activity-enhancing yet destabilizing mutations. We demonstrate that two paralogous SDR superfamily enzymes with different specificities can form mixed heterotetramers that combine their individual enzymological properties. This suggests that oligomerization can also diversify the functions generated by a given metabolic investment, enhancing the fitness advantage provided by this architectural strategy.

SIPGCN: A Novel Deep Learning Model for Predicting Self-Interacting Proteins from Sequence Information Using Graph Convolutional Networks

Protein is the basic organic substance that constitutes the cell and is the material condition for the life activity and the guarantee of the biological function activity. Elucidating the interactions and functions of proteins is a central task in exploring the mysteries of life. As an important protein interaction, self-interacting protein (SIP) has a critical role. The fast growth of high-throughput experimental techniques among biomolecules has led to a massive influx of available SIP data. How to conduct scientific research using the massive amount of SIP data has become a new challenge that is being faced in related research fields such as biology and medicine. In this work, we design an SIP prediction method SIPGCN using a deep learning graph convolutional network (GCN) based on protein sequences. First, protein sequences are characterized using a position-specific scoring matrix, which is able to describe the biological evolutionary message, then their hidden features are extracted by the deep learning method GCN, and, finally, the random forest is utilized to predict whether there are interrelationships between proteins. In the cross-validation experiment, SIPGCN achieved 93.65% accuracy and 99.64% specificity in the human data set. SIPGCN achieved 90.69% and 99.08% of these two indicators in the yeast data set, respectively. Compared with other feature models and previous methods, SIPGCN showed excellent results. These outcomes suggest that SIPGCN may be a suitable instrument for predicting SIP and may be a reliable candidate for future wet experiments.

Non-adaptive complexity and biochemical function

Intricate biochemical structures are usually thought to be useful, because natural selection preserves them from degradation by a constant hail of destructive mutations. Biochemists therefore often deliberately disrupt them to understand how complexity improves protein function or fitness. However, evolutionary theory suggests that even useless complexity that never improved fitness can become completely essential if a simple set of evolutionary conditions is fulfilled. We review evidence that stable protein complexes, protein-chaperone interactions, and complexes consisting of several paralogs all fulfill these conditions. This makes reverse genetics or destructive mutagenesis unsuitable for assigning functions to these kinds of complexity. Instead, we advocate that incorporating evolutionary approaches into biochemistry overcomes this difficulty and allows us to distinguish useless from useful biochemical complexity.

Multiple Copies of flhDC in Paraburkholderia unamae Regulate Flagellar Gene Expression, Motility, and Biofilm Formation

FlhDC is a heterohexameric complex that acts as a master regulator of flagellar biosynthesis genes in numerous bacteria. Previous studies have identified a single flhDC operon encoding this complex. However, we found that two flhDC loci are present throughout Paraburkholderia, and two additional flhC copies are also present in Paraburkholderia unamae. Systematic deletion analysis in P. unamae of the different flhDC copies showed that one of the operons, flhDC1, plays the predominant role, with deletion of its genes resulting in a severe inhibition of motility and biofilm formation. Expression analysis using promoter-lacZ fusions and real-time quantitative PCR support the primary role of flhDC1 in flagellar gene regulation, with flhDC2 a secondary contributor. Phylogenetic analysis shows the presence of the flhDC1 and flhDC2 operons throughout Paraburkholderia. In contrast, Burkholderia and other bacteria only carry the copy syntenous with flhDC2. The variations in impact each copy of flhDC has on downstream processes indicate that regulation of FlhDC in P. unamae, and likely other Paraburkholderia species, is regulated at least in part by the presence of multiple copies of these genes. IMPORTANCE Motility is important in the colonization of plant roots by beneficial and pathogenic bacteria, with flagella playing essential roles in host cell adhesion, entrance, and biofilm formation. Flagellar biosynthesis is energetically expensive. Its complex regulation by the FlhDC master regulator is well studied in peritrichous flagella expressing enterics. We report the unique presence throughout Paraburkholderia of multiple copies of flhDC. In P. unamae, the flhDC1 copy showed higher expression and a greater effect on swim motility, flagellar development, and regulation of downstream genes, than the flhDC2 copy that is syntenous to flhDC in Escherichia coli and pathogenic Burkholderia spp. The flhDC genes have evolved differently in these plant-growth-promoting bacteria, giving an additional layer of complexity in gene regulation by FlhDC.

Gene duplication and subsequent diversification strongly affect phenotypic evolvability and robustness

We study the effects of non-determinism and gene duplication on the structure of genotype-phenotype (GP) maps by introducing a non-deterministic version of the Polyomino self-assembly model. This model has previously been used in a variety of contexts to model the assembly and evolution of protein quaternary structure. Firstly, we show the limit of the current deterministic paradigm which leads to built-in anti-correlation between evolvability and robustness at the genotypic level. We develop a set of metrics to measure structural properties of GP maps in a non-deterministic setting and use them to evaluate the effects of gene duplication and subsequent diversification. Our generalized versions of evolvability and robustness exhibit positive correlation for a subset of genotypes. This positive correlation is only possible because non-deterministic phenotypes can contribute to both robustness and evolvability. Secondly, we show that duplication increases robustness and reduces evolvability initially, but that the subsequent diversification that duplication enables has a stronger, inverse effect, greatly increasing evolvability and reducing robustness relative to their original values.

Expression attenuation as a mechanism of robustness against gene duplication

Gene duplication is ubiquitous and a major driver of phenotypic diversity across the tree of life, but its immediate consequences are not fully understood. Deleterious effects would decrease the probability of retention of duplicates and prevent their contribution to long-term evolution. One possible detrimental effect of duplication is the perturbation of the stoichiometry of protein complexes. Here, we measured the fitness effects of the duplication of 899 essential genes in the budding yeast using high-resolution competition assays. At least 10% of genes caused a fitness disadvantage when duplicated. Intriguingly, the duplication of most protein complex subunits had small to nondetectable effects on fitness, with few exceptions. We selected four complexes with subunits that had an impact on fitness when duplicated and measured the impact of individual gene duplications on their protein-protein interactions. We found that very few duplications affect both fitness and interactions. Furthermore, large complexes such as the 26S proteasome are protected from gene duplication by attenuation of protein abundance. Regulatory mechanisms that maintain the stoichiometric balance of protein complexes may protect from the immediate effects of gene duplication. Our results show that a better understanding of protein regulation and assembly in complexes is required for the refinement of current models of gene duplication.

Reconstructed evolutionary history of the yeast septins Cdc11 and Shs1

Septins are GTP-binding proteins conserved across metazoans. They can polymerize into extended filaments and, hence, are considered a component of the cytoskeleton. The number of individual septins varies across the tree of life—yeast (Saccharomyces cerevisiae) has seven distinct subunits, a nematode (Caenorhabditis elegans) has two, and humans have 13. However, the overall geometric unit (an apolar hetero-octameric protomer and filaments assembled there from) has been conserved. To understand septin evolutionary variation, we focused on a related pair of yeast subunits (Cdc11 and Shs1) that appear to have arisen from gene duplication within the fungal clade. Either Cdc11 or Shs1 occupies the terminal position within a hetero-octamer, yet Cdc11 is essential for septin function and cell viability, whereas Shs1 is not. To discern the molecular basis of this divergence, we utilized ancestral gene reconstruction to predict, synthesize, and experimentally examine the most recent common ancestor (“Anc.11-S”) of Cdc11 and Shs1. Anc.11-S was able to occupy the terminal position within an octamer, just like the modern subunits. Although Anc.11-S supplied many of the known functions of Cdc11, it was unable to replace the distinct function(s) of Shs1. To further evaluate the history of Shs1, additional intermediates along a proposed trajectory from Anc.11-S to yeast Shs1 were generated and tested. We demonstrate that multiple events contributed to the current properties of Shs1: (1) loss of Shs1–Shs1 self-association early after duplication, (2) co-evolution of heterotypic Cdc11–Shs1 interaction between neighboring hetero-octamers, and (3) eventual repurposing and acquisition of novel function(s) for its C-terminal extension domain. Thus, a pair of duplicated proteins, despite constraints imposed by assembly into a highly conserved multi-subunit structure, could evolve new functionality via a complex evolutionary pathway.

Galvanization of Protein-Protein Interactions in a Dynamic Zinc Interactome

The presence of Zn²⁺ at protein-protein interfaces modulates complex function, stability, and introduces structural flexibility/complexity, chemical selectivity, and reversibility driven in a Zn²⁺-dependent manner. Recent studies have demonstrated that dynamically changing Zn²⁺ affects numerous cellular processes, including protein-protein communication and protein complex assembly. How Zn²⁺-involved protein-protein interactions (ZPPIs) are formed and dissociate and how their stability and reactivity are driven in a zinc interactome remain poorly understood, mostly due to experimental obstacles. Here, we review recent research advances on the role of Zn²⁺ in the formation of interprotein sites, their architecture, function, and stability. Moreover, we underline the importance of zinc networks in intersystemic communication and highlight bioinformatic and experimental challenges required for the identification and investigation of ZPPIs.

Determining the interaction status and evolutionary fate of duplicated homomeric proteins

Oligomeric proteins are central to life. Duplication and divergence of their genes is a key evolutionary driver, also because duplications can yield very different outcomes. Given a homomeric ancestor, duplication can yield two paralogs that form two distinct homomeric complexes, or a heteromeric complex comprising both paralogs. Alternatively, one paralog remains a homomer while the other acquires a new partner. However, so far, conflicting trends have been noted with respect to which fate dominates, primarily because different methods and criteria are being used to assign the interaction status of paralogs. Here, we systematically analyzed all Saccharomyces cerevisiae and Escherichia coli oligomeric complexes that include paralogous proteins. We found that the proportions of homo-hetero duplication fates strongly depend on a variety of factors, yet that nonetheless, rigorous filtering gives a consistent picture. In E. coli about 50%, of the paralogous pairs appear to have retained the ancestral homomeric interaction, whereas in S. cerevisiae only ~10% retained a homomeric state. This difference was also observed when unique complexes were counted instead of paralogous gene pairs. We further show that this difference is accounted for by multiple cases of heteromeric yeast complexes that share common ancestry with homomeric bacterial complexes. Our analysis settles contradicting trends and conflicting previous analyses, and provides a systematic and rigorous pipeline for delineating the fate of duplicated oligomers in any organism for which protein-protein interaction data are available.

About half of all proteins assemble as oligomers, either by self-interaction (homomers) or via interaction with another protein (heteromers). The latter can be unrelated, yet, quite commonly, the interacting proteins are paralogs, namely two genes that arose by gene duplication. Indeed, while a homomer is encoded by a single gene, heteromers demand two genes as a minimum. Duplication can therefore yield two discrete homomeric complexes or a single heteromer. Do paralogs tend to retain the ancestral homomeric interaction, or do they mostly diverge into heteromeric complexes? Despite several studies addressing this question, to date, we lack a systematic, rigorous approach for delineating the oligomeric fates of paralogs on an organism scale. To this end, we developed a new pipeline for analysis of molecular interaction databases that includes various filtering steps and unambiguous definitions of all possible oligomeric fates. Applying this method to Escherichia coli and Saccharomyces cerevisiae we noted that paralogous pairs tend to remain homomeric in the former while in the latter heteromeric complexes dominate. We consequently note a systematic trend of homomeric bacterial proteins diverging into heteromeric complexes in eukaryotes.

Tetrameric Structures of Inorganic CBS-Pyrophosphatases from Various Bacterial Species Revealed by Small-Angle X-ray Scattering in Solution

Quaternary structure of CBS-pyrophosphatases (CBS-PPases), which belong to the PPases of family II, plays an important role in their function ensuring cooperative behavior of the enzymes. Despite an intensive research, high resolution structures of the full-length CBS-PPases are not yet available making it difficult to determine the signal transmission path from the regulatory to the active center. In the present work, small-angle X-ray scattering (SAXS) combined with size-exclusion chromatography was applied to determine the solution structures of the full-length wild-type CBS-PPases from three different bacterial species. Previously, in the absence of an experimentally determined full-length CBS-PPase structure, a homodimeric model of the enzyme based on known crystal structures of the CBS domain and family II PPase without this domain has been proposed. Our SAXS analyses demonstrate, for the first time, the existence of stable tetramers in solution for all studied CBS-PPases from different sources. Our findings show that further studies are required to establish the functional properties of these enzymes. This is important not only to enhance our understanding of the relation between CBS-PPases structure and function under normal conditions but also because some human pathogens harbor this class of enzymes.

Propulsive nanomachines: the convergent evolution of archaella, flagella and cilia

Echoing the repeated convergent evolution of flight and vision in large eukaryotes, propulsive swimming motility has evolved independently in microbes in each of the three domains of life. Filamentous appendages – archaella in Archaea, flagella in Bacteria and cilia in Eukaryotes – wave, whip or rotate to propel microbes, overcoming diffusion and enabling colonization of new environments. The implementations of the three propulsive nanomachines are distinct, however: archaella and flagella rotate, while cilia beat or wave; flagella and cilia assemble at their tips, while archaella assemble at their base; archaella and cilia use ATP for motility, while flagella use ion-motive force. These underlying differences reflect the tinkering required to evolve a molecular machine, in which pre-existing machines in the appropriate contexts were iteratively co-opted for new functions and whose origins are reflected in their resultant mechanisms. Contemporary homologies suggest that archaella evolved from a non-rotary pilus, flagella from a non-rotary appendage or secretion system, and cilia from a passive sensory structure. Here, we review the structure, assembly, mechanism and homologies of the three distinct solutions as a foundation to better understand how propulsive nanomachines evolved three times independently and to highlight principles of molecular evolution.

Evolution of metabolic capabilities and molecular features of diplonemids, kinetoplastids, and euglenids

Background

The Euglenozoa are a protist group with an especially rich history of evolutionary diversity. They include diplonemids, representing arguably the most species-rich clade of marine planktonic eukaryotes; trypanosomatids, which are notorious parasites of medical and veterinary importance; and free-living euglenids. These different lifestyles, and particularly the transition from free-living to parasitic, likely require different metabolic capabilities. We carried out a comparative genomic analysis across euglenozoan diversity to see how changing repertoires of enzymes and structural features correspond to major changes in lifestyles.

Results

We find a gradual loss of genes encoding enzymes in the evolution of kinetoplastids, rather than a sudden decrease in metabolic capabilities corresponding to the origin of parasitism, while diplonemids and euglenids maintain more metabolic versatility. Distinctive characteristics of molecular machines such as kinetochores and the pre-replication complex that were previously considered specific to parasitic kinetoplastids were also identified in their free-living relatives. Therefore, we argue that they represent an ancestral rather than a derived state, as thought until the present. We also found evidence of ancient redundancy in systems such as NADPH-dependent thiol-redox. Only the genus Euglena possesses the combination of trypanothione-, glutathione-, and thioredoxin-based systems supposedly present in the euglenozoan common ancestor, while other representatives of the phylum have lost one or two of these systems. Lastly, we identified convergent losses of specific metabolic capabilities between free-living kinetoplastids and ciliates. Although this observation requires further examination, it suggests that certain eukaryotic lineages are predisposed to such convergent losses of key enzymes or whole pathways.

Conclusions

The loss of metabolic capabilities might not be associated with the switch to parasitic lifestyle in kinetoplastids, and the presence of a highly divergent (or unconventional) kinetochore machinery might not be restricted to this protist group. The data derived from the transcriptomes of free-living early branching prokinetoplastids suggests that the pre-replication complex of Trypanosomatidae is a highly divergent version of the conventional machinery. Our findings shed light on trends in the evolution of metabolism in protists in general and open multiple avenues for future research.

Paralog dependency indirectly affects the robustness of human cells

The protective redundancy of paralogous genes partly relies on the fact that they carry their functions independently. However, a significant fraction of paralogous proteins may form functionally dependent pairs, for instance, through heteromerization. As a consequence, one could expect these heteromeric paralogs to be less protective against deleterious mutations. To test this hypothesis, we examined the robustness landscape of gene loss-of-function by CRISPR-Cas9 in more than 450 human cell lines. This landscape shows regions of greater deleteriousness to gene inactivation as a function of key paralog properties. Heteromeric paralogs are more likely to occupy such regions owing to their high expression and large number of protein-protein interaction partners. Further investigation revealed that heteromers may also be under stricter dosage balance, which may also contribute to the higher deleteriousness upon gene inactivation. Finally, we suggest that physical dependency may contribute to the deleteriousness upon loss-of-function as revealed by the correlation between the strength of interactions between paralogs and their higher deleteriousness upon loss of function.

The role of structural pleiotropy and regulatory evolution in the retention of heteromers of paralogs

Gene duplication is a driver of the evolution of new functions. The duplication of genes encoding homomeric proteins leads to the formation of homomers and heteromers of paralogs, creating new complexes after a single duplication event. The loss of these heteromers may be required for the two paralogs to evolve independent functions. Using yeast as a model, we find that heteromerization is frequent among duplicated homomers and correlates with functional similarity between paralogs. Using in silico evolution, we show that for homomers and heteromers sharing binding interfaces, mutations in one paralog can have structural pleiotropic effects on both interactions, resulting in highly correlated responses of the complexes to selection. Therefore, heteromerization could be preserved indirectly due to selection for the maintenance of homomers, thus slowing down functional divergence between paralogs. We suggest that paralogs can overcome the obstacle of structural pleiotropy by regulatory evolution at the transcriptional and post-translational levels.

Evolution of interface binding strengths in simplified model of protein quaternary structure

The self-assembly of proteins into protein quaternary structures is of fundamental importance to many biological processes, and protein misassembly is responsible for a wide range of proteopathic diseases. In recent years, abstract lattice models of protein self-assembly have been used to simulate the evolution and assembly of protein quaternary structure, and to provide a tractable way to study the genotype-phenotype map of such systems. Here we generalize these models by representing the interfaces as mutable binary strings. This simple change enables us to model the evolution of interface strengths, interface symmetry, and deterministic assembly pathways. Using the generalized model we are able to reproduce two important results established for real protein complexes: The first is that protein assembly pathways are under evolutionary selection to minimize misassembly. The second is that the assembly pathway of a complex mirrors its evolutionary history, and that both can be derived from the relative strengths of interfaces. These results demonstrate that the generalized lattice model offers a powerful new idealized framework to facilitate the study of protein self-assembly processes and their evolution.

Mosaic origin of the eukaryotic kinetochore

The emergence of eukaryotes from ancient prokaryotic lineages embodied a remarkable increase in cellular complexity. While prokaryotes operate simple systems to connect DNA to the segregation machinery during cell division, eukaryotes use a highly complex protein assembly known as the kinetochore. Although conceptually similar, prokaryotic segregation systems and the eukaryotic kinetochore are not homologous. Here we investigate the origins of the kinetochore before the last eukaryotic common ancestor (LECA) using phylogenetic trees, sensitive profile-versus-profile homology detection, and structural comparisons of its protein components. We show that LECA's kinetochore proteins share deep evolutionary histories with proteins involved in a few prokaryotic systems and a multitude of eukaryotic processes, including ubiquitination, transcription, and flagellar and vesicular transport systems. We find that gene duplications played a major role in shaping the kinetochore; more than half of LECA's kinetochore proteins have other kinetochore proteins as closest homologs. Some of these have no detectable homology to any other eukaryotic protein, suggesting that they arose as kinetochore-specific folds before LECA. We propose that the primordial kinetochore evolved from proteins involved in various (pre)eukaryotic systems as well as evolutionarily novel folds, after which a subset duplicated to give rise to the complex kinetochore of LECA.

Structural principles that enable oligomeric small heat-shock protein paralogs to evolve distinct functions

Oligomeric proteins assemble with exceptional selectivity, even in the presence of closely related proteins, to perform their cellular roles. We show that most proteins related by gene duplication of an oligomeric ancestor have evolved to avoid hetero-oligomerization and that this correlates with their acquisition of distinct functions. We report how coassembly is avoided by two oligomeric small heat-shock protein paralogs. A hierarchy of assembly, involving intermediates that are populated only fleetingly at equilibrium, ensures selective oligomerization. Conformational flexibility at noninterfacial regions in the monomers prevents coassembly, allowing interfaces to remain largely conserved. Homomeric oligomers must overcome the entropic benefit of coassembly and, accordingly, homomeric paralogs comprise fewer subunits than homomers that have no paralogs.

Inferring and Using Protein Quaternary Structure Information from Crystallographic Data

A precise knowledge of the quaternary structure of proteins is essential to illuminate both their function and their evolution. The major part of our knowledge on quaternary structure is inferred from X-ray crystallography data, but this inference process is hard and error-prone. The difficulty lies in discriminating fortuitous protein contacts, which make up the lattice of protein crystals, from biological protein contacts that exist in the native cellular environment. Here, we review methods devised to discriminate between both types of contacts and describe resources for downloading protein quaternary structure information and identifying high-confidence quaternary structures. The use of high-confidence datasets of quaternary structures will be critical for the analysis of structural, functional, and evolutionary properties of proteins.

Transcription regulation by the Mediator complex

Alterations in the regulation of gene expression are frequently associated with developmental diseases or cancer. Transcription activation is a key phenomenon in the regulation of gene expression. In all eukaryotes, mediator of RNA polymerase II transcription (Mediator), a large complex with modular organization, is generally required for transcription by RNA polymerase II, and it regulates various steps of this process. The main function of Mediator is to transduce signals from the transcription activators bound to enhancer regions to the transcription machinery, which is assembled at promoters as the preinitiation complex (PIC) to control transcription initiation. Recent functional studies of Mediator with the use of structural biology approaches and functional genomics have revealed new insights into Mediator activity and its regulation during transcription initiation, including how Mediator is recruited to transcription regulatory regions and how it interacts and cooperates with PIC components to assist in PIC assembly. Novel roles of Mediator in the control of gene expression have also been revealed by showing its connection to the nuclear pore and linking Mediator to the regulation of gene positioning in the nuclear space. Clear links between Mediator subunits and disease have also encouraged studies to explore targeting of this complex as a potential therapeutic approach in cancer and fungal infections.

Sharing the SAGA

Transcription initiation is a major regulatory step in eukaryotic gene expression. Co-activators establish transcriptionally competent promoter architectures and chromatin signatures to allow the formation of the pre-initiation complex (PIC), comprising RNA polymerase II (Pol II) and general transcription factors (GTFs). Many GTFs and co-activators are multisubunit complexes, in which individual components are organized into functional modules carrying specific activities. Recent advances in affinity purification and mass spectrometry analyses have revealed that these complexes often share functional modules, rather than containing unique components. This observation appears remarkably prevalent for chromatin-modifying and remodeling complexes. Here, we use the modular organization of the evolutionary conserved Spt-Ada-Gcn5 acetyltransferase (SAGA) complex as a paradigm to illustrate how co-activators share and combine a relatively limited set of functional tools.

Gene duplication can impart fragility, not robustness, in the yeast protein interaction network

The maintenance of duplicated genes is thought to protect cells from genetic perturbations, but the molecular basis of this robustness is largely unknown. By measuring the interaction of yeast proteins with their partners in wild-type cells and in cells lacking a paralog, we found that 22 out of 56 paralog pairs compensate for the lost interactions. An equivalent number of pairs exhibit the opposite behavior and require each other's presence for maintaining their interactions. These dependent paralogs generally interact physically, regulate each other's abundance, and derive from ancestral self-interacting proteins. This reveals that gene duplication may actually increase mutational fragility instead of robustness in a large number of cases.

Signalling assemblies: the odds of symmetry

The assembly of proteins into complexes is fundamental to nearly all biological signalling processes. Symmetry is a dominant feature of the structures of experimentally determined protein complexes, observed in the vast majority of homomers and many heteromers. However, some asymmetric structures exist, and asymmetry also often forms transiently, intractable to traditional structure determination methods. Here, we explore the role of protein complex symmetry and asymmetry in cellular signalling, focusing on receptors, transcription factors and transmembrane channels, among other signalling assemblies. We highlight a recurrent tendency for asymmetry to be crucial for signalling function, often being associated with activated states. We conclude with a discussion of how consideration of protein complex symmetry and asymmetry has significant potential implications and applications for pharmacology and human disease.

Evolution of complex adaptations in molecular systems

A central challenge in evolutionary biology concerns the mechanisms by which complex adaptations arise. Such adaptations depend on the fixation of multiple, highly specific mutations, where intermediate stages of evolution seemingly provide little or no benefit. It is generally assumed that the establishment of complex adaptations is very slow in nature, as evolution of such traits demands special population genetic or environmental circumstances. However, blueprints of complex adaptations in molecular systems are pervasive, indicating that they can readily evolve. We discuss the prospects and limitations of non-adaptive scenarios, which assume multiple neutral or deleterious steps in the evolution of complex adaptations. Next, we examine how complex adaptations can evolve by natural selection in changing environment. Finally, we argue that molecular 'springboards', such as phenotypic heterogeneity and promiscuous interactions facilitate this process by providing access to new adaptive paths.

Quaternary structure is an essential component that contributes to the sophisticated allosteric regulation mechanism in a key enzyme from Mycobacterium tuberculosis

The first enzyme of the shikimate pathway, 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase (DAH7PS), adopts a range of distinct allosteric regulation mechanisms in different organisms, related to different quaternary assemblies. DAH7PS from Mycobacterium tuberculosis (MtuDAH7PS) is a homotetramer, with the allosteric sites in close proximity to the interfaces. Here we examine the importance of the quaternary structure on catalysis and regulation, by amino acid substitution targeting the tetramer interface of MtuDAH7PS. Using only single amino acid substitutions either in, or remote from the interface, two dimeric variants of MtuDAH7PS (MtuDAH7PS_F227D and MtuDAH7PS_G232P) were successfully generated. Both dimeric variants maintained activity due to the distance between the sites of amino acid substitution and the active sites, but attenuated catalytic efficiency was observed. Both dimeric variants showed significantly disrupted allosteric regulation with greatly impaired binding affinity for one of the allosteric ligands. Molecular dynamics simulations revealed changes in protein dynamics and average conformations in the dimeric variant caused by amino acid substitution remote to the tetramer interface (MtuDAH7PS_G232P), which are consistent with the observed reduction in catalytic efficiency and loss of allosteric response.

Unique Phylogenetic Distributions of the Ska and Dam1 Complexes Support Functional Analogy and Suggest Multiple Parallel Displacements of Ska by Dam1

Faithful chromosome segregation relies on kinetochores, the large protein complexes that connect chromatin to spindle microtubules. Although human and yeast kinetochores are largely homologous, they track microtubules with the unrelated protein complexes Ska (Ska-C, human) and Dam1 (Dam1-C, yeast). We here uncovered that Ska-C and Dam1-C are both widespread among eukaryotes, but in an exceptionally inverse manner, supporting their functional analogy. Within the complexes, all Ska-C and various Dam1-C subunits are ancient paralogs, showing that gene duplication shaped these complexes. We examined various evolutionary scenarios to explain the nearly mutually exclusive patterns of Ska-C and Dam1-C in present-day species. We propose that Ska-C was present in the last eukaryotic common ancestor, that subsequently Dam1-C displaced Ska-C in an early fungus and was horizontally transferred to diverse non-fungal lineages, displacing Ska-C in these lineages too.

Structural reconstruction of protein ancestry

Ancestral protein reconstruction allows the resurrection and characterization of ancient proteins based on computational analyses of sequences of modern-day proteins. Unfortunately, many protein families are highly divergent and not suitable for sequence-based reconstruction approaches. This limitation is exemplified by the antigen receptors of jawed vertebrates (B- and T-cell receptors), heterodimers formed by pairs of Ig domains. These receptors are believed to have evolved from an extinct homodimeric ancestor through a process of gene duplication and diversification; however molecular evidence has so far remained elusive. Here, we use a structural approach and laboratory evolution to reconstruct such molecules and characterize their interaction with antigen. High-resolution crystal structures of reconstructed homodimeric receptors in complex with hen-egg white lysozyme demonstrate how nanomolar affinity binding of asymmetrical antigen is enabled through selective recruitment and structural plasticity within the receptor-binding site. Our results provide structural evidence in support of long-held theories concerning the evolution of antigen receptors, and provide a blueprint for the experimental reconstruction of protein ancestry in the absence of phylogenetic evidence.

Oligomerization process of Bcl-2 associated X protein revealed from intermediate structures in solution

Linear oligomerization of ditopic BAX-dimers into tri-dimer helical units then into a rod-like structure, as revealed using integrated ESR/SAXS/MD analyses.

The origin and evolution of cell types

Cell types are the basic building blocks of multicellular organisms and are extensively diversified in animals. Despite recent advances in characterizing cell types, classification schemes remain ambiguous. We propose an evolutionary definition of a cell type that allows cell types to be delineated and compared within and between species. Key to cell type identity are evolutionary changes in the 'core regulatory complex' (CoRC) of transcription factors, that make emergent sister cell types distinct, enable their independent evolution and regulate cell type-specific traits termed apomeres. We discuss the distinction between developmental and evolutionary lineages, and present a roadmap for future research.

Protein-protein interactions can be predicted using coiled coil co-evolution patterns

Protein-protein interactions are sometimes mediated by coiled coil structures. The evolutionary conservation of interacting orthologs in different species, along with the presence or absence of coiled coils in them, may help in the prediction of interacting pairs. Here, we illustrate how the presence of coiled coils in a protein can be exploited as a potential indicator for its interaction with another protein with coiled coils. The prediction capability of our strategy improves when restricting our dataset to highly reliable, known protein-protein interactions. Our study of the co-evolution of coiled coils demonstrates that pairs of interacting proteins can be distinguished from not interacting pairs by means of their structural information. This hints at the potential of our strategy to predict new protein-protein interactions.

Comparative analysis of genome-wide Mlo gene family in Cajanus cajan and Phaseolus vulgaris

The Mlo gene was discovered in barley because the mutant 'mlo' allele conferred broad-spectrum, non-race-specific resistance to powdery mildew caused by Blumeria graminis f. sp. hordei. The Mlo genes also play important roles in growth and development of plants, and in responses to biotic and abiotic stresses. The Mlo gene family has been characterized in several crop species, but only a single legume species, soybean (Glycine max L.), has been investigated so far. The present report describes in silico identification of 18 CcMlo and 20 PvMlo genes in the important legume crops Cajanus cajan (L.) Millsp. and Phaseolus vulgaris L., respectively. In silico analysis of gene organization, protein properties and conserved domains revealed that the C. cajan and P. vulgaris Mlo gene paralogs are more divergent from each other than from their orthologous pairs. The comparative phylogenetic analysis classified CcMlo and PvMlo genes into three major clades. A comparative analysis of CcMlo and PvMlo proteins with the G. max Mlo proteins indicated close association of one CcMlo, one PvMlo with two GmMlo genes, indicating that there was no further expansion of the Mlo gene family after the separation of these species. Thus, most of the diploid species of eudicots might be expected to contain 15-20 Mlo genes. The genes CcMlo12 and 14, and PvMlo11 and 12 are predicted to participate in powdery mildew resistance. If this prediction were verified, these genes could be targeted by TILLING or CRISPR to isolate powdery mildew resistant mutants.

Paralogue Interference Affects the Dynamics after Gene Duplication

Proteins tend to form homomeric complexes of identical subunits, which act as functional units. By definition, the subunits are encoded from a single genetic locus. When such a gene is duplicated, the gene products are suggested initially to cross-interact when coexpressed, thus resulting in the phenomenon of paralogue interference. In this opinion article, we explore how paralogue interference can shape the fate of a duplicated gene. One important outcome is a prolonged time window in which both copies remain under selection increasing the chance to accumulate mutations and to develop new properties. Thereby, paralogue interference can mediate the coevolution of duplicates and here we illustrate the potential of this phenomenon in light of recent new studies.

Substrate-Induced Dimerization of Engineered Monomeric Variants of Triosephosphate Isomerase from Trichomonas vaginalis

The dimeric nature of triosephosphate isomerases (TIMs) is maintained by an extensive surface area interface of more than 1600 Ų. TIMs from Trichomonas vaginalis (TvTIM) are held in their dimeric state by two mechanisms: a ball and socket interaction of residue 45 of one subunit that fits into the hydrophobic pocket of the complementary subunit and by swapping of loop 3 between subunits. TvTIMs differ from other TIMs in their unfolding energetics. In TvTIMs the energy necessary to unfold a monomer is greater than the energy necessary to dissociate the dimer. Herein we found that the character of residue I45 controls the dimer-monomer equilibrium in TvTIMs. Unfolding experiments employing monomeric and dimeric mutants led us to conclude that dimeric TvTIMs unfold following a four state model denaturation process whereas monomeric TvTIMs follow a three state model. In contrast to other monomeric TIMs, monomeric variants of TvTIM1 are stable and unexpectedly one of them (I45A) is only 29-fold less active than wild-type TvTIM1. The high enzymatic activity of monomeric TvTIMs contrast with the marginal catalytic activity of diverse monomeric TIMs variants. The stability of the monomeric variants of TvTIM1 and the use of cross-linking and analytical ultracentrifugation experiments permit us to understand the differences between the catalytic activities of TvTIMs and other marginally active monomeric TIMs. As TvTIMs do not unfold upon dimer dissociation, herein we found that the high enzymatic activity of monomeric TvTIM variants is explained by the formation of catalytic dimeric competent species assisted by substrate binding.