Ancestral reconstruction of duplicated signaling proteins reveals the evolution of signaling specificity

Evolutionary paths that link orthogonal pairs of binding proteins

Some protein-binding pairs exhibit extreme specificities that functionally insulate them from homologs. Such pairs evolve mostly by accumulating single-point mutations, and mutants are selected if they exhibit sufficient affinity. Until now, finding a fully functional single-mutation path connecting orthogonal pairs could only be achieved by full enumeration of intermediates and was restricted to pairs that were mutationally close. We present a computational framework for discovering single-mutation paths with low molecular strain and apply it to two orthogonal bacterial endonuclease-immunity pairs separated by 17 interfacial mutations. By including mutations that bridge identities that could not be exchanged by single-nucleotide mutations, we discovered a strain-free 19-mutation path that was fully functional in vivo. The change in binding preference occurred remarkably abruptly, resulting from only one radical mutation in each partner. Furthermore, each of the specificity-switch mutations increased fitness, demonstrating that functional divergence could be driven by positive Darwinian selection.

Transcriptomic and proteomic ramifications of segmental amplification in Escherichia coli

Gene amplification can drive adaptation by rapidly increasing the cellular dosage of critical gene products. Segmental amplifications often encompass large genomic regions surrounding the gene(s) under selection for higher dosage. Overexpression of coamplified neighboring genes imposes a substantial metabolic burden. While compensatory mutations can decrease inappropriate overexpression of coamplified genes, it takes time for such mutations to arise. The extent to which intrinsic regulatory mechanisms modulate expression of coamplified genes in the immediate aftermath of segmental amplification is largely unknown. To address the collateral effects of segmental amplification, we evolved replicate cultures of an Escherichia coli mutant under conditions that select for higher dosage of an inefficient enzyme whose weak activity limits growth rate. Segmental amplifications encompassing the gene encoding the weak-link enzyme arose in all populations. Amplified regions ranged in size (33 to 125 kb) and copy number (2 to ≥14 copies). We performed RNA-seq and label-free proteomics to quantify expression of amplified genes present at 2, 6, and 14 copies. mRNA expression generally scales with gene copy number, but protein expression scales less well with both gene copy number and mRNA expression. We characterize the molecular mechanisms underlying discrepancies between gene copy number and expression for several cases. We also show that segmental amplifications can have system-wide consequences by indirectly altering expression of nonamplified genes. Our findings indicate that the fitness benefit derived from segmental amplification depends on the combined effects of amplicon size, gene content, and copy number as well as collateral effects on nonamplified genes.

High-throughput amino acid-level characterization of the interactions of plasminogen activator inhibitor-1 with variably divergent proteases

While members of large paralogous protein families share structural features, their functional niches often diverge significantly. Serine protease inhibitors (SERPINs), whose members typically function as covalent inhibitors of serine proteases, are one such family. Plasminogen activator inhibitor-1 (PAI-1) is a prototypic SERPIN, which canonically inhibits tissue- and urokinase-type plasminogen activators (tPA and uPA) to regulate fibrinolysis. PAI-1 has been shown to also inhibit other serine proteases, including coagulation factor XIIa (FXIIa) and transmembrane serine protease 2 (TMPRSS2). The structural determinants of PAI-1 inhibitory function toward these non-canonical protease targets, and the biological significance of these functions, are unknown. We applied deep mutational scanning (DMS) to assess the effects of ~80% of all possible single-amino acid substitutions in PAI-1 on its ability to inhibit three putative serine protease targets (uPA, FXIIa, and TMPRSS2). Selection with each target protease generated a unique PAI-1 mutational landscape, with the determinants of protease specificity distributed throughout PAI-1's primary sequence. Next, we conducted a comparative analysis of extant orthologous sequences, demonstrating that key residues modulating PAI-1 inhibition of uPA and FXIIa, but not TMPRSS2, are maintained by purifying selection (also referred to as "negative selection"). PAI-1's activity toward FXIIa may reflect how protease evolutionary relationships predict SERPIN functional divergence, which we support via a cophylogenetic analysis of secreted SERPINs and their cognate serine proteases. This work provides insight into the functional diversification of SERPINs and lays the framework for extending these studies to other proteases and their regulators.

Enantiocomplementary Gut Bacterial Enzymes Metabolize Dietary Polyphenols

Molybdenum-dependent catechol dehydroxylases in gut Actinobacteria catalyze the removal of para-hydroxyl groups from catechols, a central reaction in the microbial metabolism of polyphenol compounds. However, the substrates of most putative catechol dehydroxylases remain unidentified due to the challenges of obtaining these enzymes from standard heterologous expression systems. In this work, we establish Gordonibacter urolithinfaciens as a versatile bacterial host to express active catechol dehydroxylases. Using this system, we rapidly deorphanize eight previously uncharacterized gut bacterial catechol dehydroxylases that selectively dehydroxylate intermediates in the gut bacterial metabolism of plant-derived catechins and lignans. Unexpectedly, we discover multiple instances of distinct catechol dehydroxylases that have evolved to selectively metabolize individual substrate enantiomers, setting the stage for future efforts to elucidate their mechanisms and evolution. Altogether, these findings greatly increase our knowledge of these metalloenzymes, illustrating the power of bacterial genetics to accelerate enzyme discovery and providing a more complete understanding of transformations relevant to the health benefits of phytochemicals.

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.

Zebrafish: unraveling genetic complexity through duplicated genes

The zebrafish is an invaluable model organism for genetic, developmental, and disease research. Although its high conservation with humans is often cited as justification for its use, the zebrafish harbors oft-ignored genetic characteristics that may provide unique insights into gene structure and function. Zebrafish, along with other teleost fish, underwent an additional round of whole genome duplication after their split from tetrapods-resulting in an abundance of duplicated genes when compared to other vertebrates. These duplicated genes have evolved in distinct ways over the ensuing 350 million years. Thus, each gene within a duplicated gene pair has nuanced differences that create a unique identity. By investigating both members of the gene pair together, we can elucidate the mechanisms that underly protein structure and function and drive the complex interplay within biological systems, such as signal transduction cascades, genetic regulatory networks, and evolution of tissue and organ function. It is crucial to leverage such studies to explore these molecular dynamics, which could have far-reaching implications for both basic science and therapeutic development. Here, we will review the role of gene duplications and the existing models for gene divergence and retention following these events. We will also highlight examples within each of these models where studies comparing duplicated genes in the zebrafish have yielded key insights into protein structure, function, and regulation.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known. Here, we use Saccharomyces cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, the resulting tyrosine phosphorylation is biologically spurious. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3500 proteins. The number of spurious pY sites generated correlates strongly with decreased growth, and we predict over 1000 pY events to be deleterious. However, we also find that many of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with tyrosine kinases. Our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

High-throughput amino acid-level characterization of the interactions of plasminogen activator inhibitor-1 with variably divergent proteases

While members of large paralogous protein families share structural features, their functional niches often diverge significantly. Serine protease inhibitors (SERPINs), whose members typically function as covalent inhibitors of serine proteases, are one such family. Plasminogen activator inhibitor-1 (PAI-1) is a prototypic SERPIN, which canonically inhibits tissue-and urokinase-type plasminogen activators (tPA and uPA) to regulate fibrinolysis. PAI-1 has been shown to also inhibit other serine proteases, including coagulation factor XIIa (FXIIa) and transmembrane serine protease 2 (TMPRSS2). The structural determinants of PAI-1 inhibitory function toward these non-canonical protease targets, and the biological significance of these functions, are unknown. We applied deep mutational scanning (DMS) to assess the effects of ∼80% of all possible single amino acid substitutions in PAI-1 on its ability to inhibit three putative serine protease targets (uPA, FXIIa, and TMPRSS2). Selection with each target protease generated a unique PAI-1 mutational landscape, with the determinants of protease specificity distributed throughout PAI-1's primary sequence. Next, we conducted a comparative analysis of extant orthologous sequences, demonstrating that key residues modulating PAI-1 inhibition of uPA and FXIIa, but not TMPRSS2, are maintained by purifying selection. PAI-1's activity toward FXIIa may reflect how protease evolutionary relationships predict SERPIN functional divergence, which we support via a cophylogenetic analysis of all secreted SERPINs and their cognate serine proteases. This work provides insight into the functional diversification of SERPINs and lays the framework for extending these studies to other proteases and their regulators.

Crosstalk involving two-component systems in Staphylococcus aureus signaling networks

Staphylococcus aureus poses a serious global threat to human health due to its pathogenic nature, adaptation to environmental stress, high virulence, and the prevalence of antimicrobial resistance. The signaling network in S. aureus coordinates and integrates various internal and external inputs and stimuli to adapt and formulate a response to the environment. Two-component systems (TCSs) of S. aureus play a central role in this network where surface-expressed histidine kinases (HKs) receive and relay external signals to their cognate response regulators (RRs). Despite the purported high fidelity of signaling, crosstalk within TCSs, between HK and non-cognate RR, and between TCSs and other systems has been detected widely in bacteria. The examples of crosstalk in S. aureus are very limited, and there needs to be more understanding of its molecular recognition mechanisms, although some crosstalk can be inferred from similar bacterial systems that share structural similarities. Understanding the cellular processes mediated by this crosstalk and how it alters signaling, especially under stress conditions, may help decipher the emergence of antibiotic resistance. This review highlights examples of signaling crosstalk in bacteria in general and S. aureus in particular, as well as the effect of TCS mutations on signaling and crosstalk.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Evolutionary paths that link orthogonal pairs of binding proteins

Some protein binding pairs exhibit extreme specificities that functionally insulate them from homologs. Such pairs evolve mostly by accumulating single-point mutations, and mutants are selected if their affinity exceeds the threshold required for function1-4. Thus, homologous and high-specificity binding pairs bring to light an evolutionary conundrum: how does a new specificity evolve while maintaining the required affinity in each intermediate5,6? Until now, a fully functional single-mutation path that connects two orthogonal pairs has only been described where the pairs were mutationally close thus enabling experimental enumeration of all intermediates2. We present an atomistic and graph-theoretical framework for discovering low molecular strain single-mutation paths that connect two extant pairs, enabling enumeration beyond experimental capability. We apply it to two orthogonal bacterial colicin endonuclease-immunity pairs separated by 17 interface mutations7. We were not able to find a strain-free and functional path in the sequence space defined by the two extant pairs. But including mutations that bridge amino acids that cannot be exchanged through single-nucleotide mutations led us to a strain-free 19-mutation trajectory that is completely viable in vivo. Our experiments show that the specificity switch is remarkably abrupt, resulting from only one radical mutation on each partner. Furthermore, each of the critical specificity-switch mutations increases fitness, demonstrating that functional divergence could be driven by positive Darwinian selection. These results reveal how even radical functional changes in an epistatic fitness landscape may evolve.

Environmental adaptation and diversification of bacterial two-component systems

Bacterial two-component systems (TCS) are versatile signaling mechanisms that govern cellular responses to diverse environmental cues. These systems rely on phosphoryl-group transfers between histidine- and aspartate-containing modules of sensor histidine kinase and response regulator proteins. TCS diversity is shaped by the ecological niche of the bacterium, resulting in significant population-level variations. Consequently, orthologous TCSs can display considerable divergence throughout the signaling process. Here, we venture into the mechanisms governing the emergence of TCS variation, and explore the adaptation of orthologous TCS in bacteria with dissimilar lifestyles. The peculiar features of the bacterial adaptive response A/ultraviolet light repair Y (BarA/UvrY) and anoxic redox control B/anoxic redox control A (ArcB/ArcA) and their ortholog TCSs illustrate the remarkable capacity of TCSs to evolve and finely tune their signaling mechanisms, effectively addressing specific environmental challenges.

The permissive binding theory of cancer

The later stages of cancer, including the invasion and colonization of new tissues, are actively mysterious compared to earlier stages like primary tumor formation. While we lack many details about both, we do have an apparently successful explanatory framework for the earlier stages: one in which genetic mutations hold ultimate causal and explanatory power. By contrast, on both empirical and conceptual grounds, it is not currently clear that mutations alone can explain the later stages of cancer. Can a different type of molecular change do better? Here, I introduce the "permissive binding theory" of cancer, which proposes that novel protein binding interactions are the key causal and explanatory entity in invasion and metastasis. It posits that binding is more abundant at baseline than we observe because it is restricted in normal physiology; that any large perturbation to physiological state revives this baseline abundance, unleashing many new binding interactions; and that a subset of these cause the cellular functions at the heart of oncogenesis, especially invasion and metastasis. Significant physiological perturbations occur in cancer cells in very early stages, and generally become more extreme with progression, providing interactions that continually fuel invasion and metastasis. The theory is compatible with, but not limited to, causal roles for the diverse molecular changes observed in cancer (e.g. gene expression or epigenetic changes), as these generally act causally upstream of proteins, and so may exert their effects by changing the protein binding interactions that occur in the cell. This admits the possibility that molecular changes that appear quite different may actually converge in creating the same few protein complexes, simplifying our picture of invasion and metastasis. If correct, the theory offers a concrete therapeutic strategy: targeting the key novel complexes. The theory is straightforwardly testable by large-scale identification of protein interactions in different cancers.

Theoretical study of ArcB and its dimerization, interaction with anaerobic metabolites, and activation of ArcA

The complex metabolism of Escherichia coli has been extensively studied, including its response to oxygen availability. The ArcA/B two-component system (TCS) is the key regulator for the transition between these two environmental conditions and has been thoroughly characterized using genetic and biochemical approaches. Still, to date, limited structural data is available. The breakthrough provided by AlphaFold2 in 2021 has brought a reliable tool to the scientific community for assessing the structural features of complex proteins. In this report, we analyzed the structural aspects of the ArcA/B TCS using AlphaFold2 models. The models are consistent with the experimentally determined structures of ArcB kinase. The predicted structure of the dimeric form of ArcB is consistent with the extensive genetic and biochemical data available regarding mechanistic signal perception and regulation. The predicted interaction of the dimeric form of ArcB with its cognate response regulator (ArcA) is also consistent with both the forward and reverse phosphotransfer mechanisms. The ArcB model was used to detect putative binding cavities to anaerobic metabolites, encouraging testing of these predictions experimentally. Finally, the highly accurate models of other ArcB homologs suggest that different experimental approaches are needed to determine signal perception in kinases lacking the PAS domain. Overall, ArcB is a kinase with features that need further testing, especially in determining its crystal structure under different conditions.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known, but quantifying this is required to understand the constraints faced by cell systems as they evolve. Here, we use the model organism S. cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, most of the resulting tyrosine phosphorylation is spurious. This provides a suitable system to measure the impact of artificial protein interactions on fitness. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3,500 proteins. Examination of the fitness costs in each strain revealed a strong correlation between the number of spurious pY sites and decreased growth. Moreover, the analysis of pY effects on protein structure and on protein function revealed over 1000 pY events that we predict to be deleterious. However, we also find that a large number of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with bona fide tyrosine kinases. Taken together, our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Marginal specificity in protein interactions constrains evolution of a paralogous family

The evolution of novel functions in biology relies heavily on gene duplication and divergence, creating large paralogous protein families. Selective pressure to avoid detrimental cross-talk often results in paralogs that exhibit exquisite specificity for their interaction partners. But how robust or sensitive is this specificity to mutation? Here, using deep mutational scanning, we demonstrate that a paralogous family of bacterial signaling proteins exhibits marginal specificity, such that many individual substitutions give rise to substantial cross-talk between normally insulated pathways. Our results indicate that sequence space is locally crowded despite overall sparseness, and we provide evidence that this crowding has constrained the evolution of bacterial signaling proteins. These findings underscore how evolution selects for "good enough" rather than optimized phenotypes, leading to restrictions on the subsequent evolution of paralogs.

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

General strategies for using amino acid sequence data to guide biochemical investigation of protein function

The rapid increase of '-omics' data warrants the reconsideration of experimental strategies to investigate general protein function. Studying individual members of a protein family is likely insufficient to provide a complete mechanistic understanding of family functions, especially for diverse families with thousands of known members. Strategies that exploit large amounts of available amino acid sequence data can inspire and guide biochemical experiments, generating broadly applicable insights into a given family. Here we review several methods that utilize abundant sequence data to focus experimental efforts and identify features truly representative of a protein family or domain. First, coevolutionary relationships between residues within primary sequences can be successfully exploited to identify structurally and/or functionally important positions for experimental investigation. Second, functionally important variable residue positions typically occupy a limited sequence space, a property useful for guiding biochemical characterization of the effects of the most physiologically and evolutionarily relevant amino acids. Third, amino acid sequence variation within domains shared between different protein families can be used to sort a particular domain into multiple subtypes, inspiring further experimental designs. Although generally applicable to any kind of protein domain because they depend solely on amino acid sequences, the second and third approaches are reviewed in detail because they appear to have been used infrequently and offer immediate opportunities for new advances. Finally, we speculate that future technologies capable of analyzing and manipulating conserved and variable aspects of the three-dimensional structures of a protein family could lead to broad insights not attainable by current methods.

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.