Constraints on the expansion of paralogous protein families

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.

Elaboration of the Homer1 recognition landscape reveals incomplete divergence of paralogous EVH1 domains

Short sequences that mediate interactions with modular binding domains are ubiquitous throughout eukaryotic proteomes. Networks of short linear motifs (SLiMs) and their corresponding binding domains orchestrate many cellular processes, and the low mutational barrier to evolving novel interactions provides a way for biological systems to rapidly sample selectable phenotypes. Mapping SLiM binding specificity and the rules that govern SLiM evolution is fundamental to uncovering the pathways regulated by these networks and developing the tools to manipulate them. We used high-throughput screening of the human proteome to identify sequences that bind to the Enabled/VASP homology 1 (EVH1) domain of the postsynaptic density scaffolding protein Homer1. This expanded our understanding of the determinants of Homer EVH1 binding preferences and defined a new motif that can facilitate the discovery of additional Homer-mediated interactions. Interestingly, the Homer1 EVH1 domain preferentially binds to sequences containing an N-terminally overlapping motif that is bound by the paralogous family of Ena/VASP actin polymerases, and many of these sequences can bind to EVH1 domains from both protein families. We provide evidence from orthologous EVH1 domains in pre-metazoan organisms that the overlap in human Ena/VASP and Homer binding preferences corresponds to an incomplete divergence from a common Ena/VASP ancestor. Given this overlap in binding profiles, promiscuous sequences that can be recognized by both families either achieve specificity through extrinsic regulatory strategies or may provide functional benefits via multi-specificity. This may explain why these paralogs incompletely diverged despite the accessibility of further diverged isoforms.

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.

Elaboration of the Homer1 Recognition Landscape Reveals Incomplete Divergence of Paralogous EVH1 Domains

Short sequences that mediate interactions with modular binding domains are ubiquitous throughout eukaryotic proteomes. Networks of Short Linear Motifs (SLiMs) and their corresponding binding domains orchestrate many cellular processes, and the low mutational barrier to evolving novel interactions provides a way for biological systems to rapidly sample selectable phenotypes. Mapping SLiM binding specificity and the rules that govern SLiM evolution is fundamental to uncovering the pathways regulated by these networks and developing the tools to manipulate them. We used high-throughput screening of the human proteome to identify sequences that bind to the Enabled/VASP homology 1 (EVH1) domain of the postsynaptic density scaffolding protein Homer1. In doing so, we expanded current understanding of the determinants of Homer EVH1 binding preferences and defined a new motif that can facilitate the discovery of additional Homer-mediated interactions. Interestingly, the Homer1 EVH1 domain preferentially binds to sequences containing an N-terminally overlapping motif that is bound by the paralogous family of Ena/VASP actin polymerases, and many of these sequences can bind to EVH1 domains from both protein families. We provide evidence from orthologous EVH1 domains in pre-metazoan organisms that the overlap in human Ena/VASP and Homer binding preferences corresponds to an incomplete divergence from a common Ena/VASP ancestor. Given this overlap in binding profiles, promiscuous sequences that can be recognized by both families either achieve specificity through extrinsic regulatory strategies or may provide functional benefits via multi-specificity. This may explain why these paralogs incompletely diverged despite the accessibility of further diverged isoforms.

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.

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.

Extended disordered regions of ribosome-associated NAC proteins paralogs belong only to the germline in Drosophila melanogaster

The nascent polypeptide-associated complex (NAC) consisting of α- and β-subunits is an essential ribosome-associated protein conserved in eukaryotes. NAC is a ubiquitously expressed co-translational regulator of nascent protein folding and sorting providing for homeostasis of cellular proteins. Here we report on discovering the germline-specific NACαβ paralogs (gNACs), whose β-subunits, non-distinguishable by ordinary immunodetection, are encoded by five highly homologous gene copies, while the α-subunit is encoded by a single αNAC gene. The gNAC expression is detected in the primordial embryonic and adult gonads via immunostaining. The germline-specific α and β subunits differ from the ubiquitously expressed paralogs by the extended intrinsically disordered regions (IDRs) acquired at the N- and C-termini of the coding regions, predicted to be phosphorylated. The presence of distinct phosphorylated isoforms of gNAC-β subunits is confirmed by comparing of their profiles by 2D-isoeletrofocusing resolution before and after phosphatase treatment of testis ribosomes. We revealed that the predicted S/T sites of phosphorylation in the individual orthologous IDRs of gNAC-β sequences of Drosophila species are positionally conserved despite these disordered regions are drastically different. We propose the IDR-dependent molecular crowding and specific coordination of NAC and other proteostasis regulatory factors at the ribosomes of germinal cells. Our findings imply that there may be a functional crosstalk between the germinal and ubiquitous α- and β-subunits based on assessing their depletion effects on the fly viability and gonad development.

Evaluation of type-B RR dimerization in poplar: A mechanism to preserve signaling specificity?

Plants possess specific signaling pathways, such as the MultiStep Phosphorelay (MSP), which is involved in cytokinin and ethylene sensing, and light, drought or osmotic stress sensing. These MSP comprise histidine-aspartate kinases (HKs) as receptors, histidine phosphotransfer (HPts) proteins acting as phosphorelay proteins, and response regulators (RRs), some of which act as transcription factors (type-B RRs). In previous studies, we identified partners of the poplar osmosensing signaling pathway, composed of two HKs, three main HPts, and six type-B RRs. To date, it is unresolved as to how cytokinin or osmotic stress signal specificity is achieved in the MSP in order to generate specific responses. Here, we present a large-scale interaction study of poplar type-B RR dimerization. Using the two-hybrid assay, we were able to show the homodimerization of type-B RRs, the heterodimerization of duplicated type-B RRs, and surprisingly, a lack of interaction between some type-B RRs belonging to different duplicates. The lack of interaction of the duplicates RR12-14 and RR18-19, which are involved in the osmosensing pathway has been confirmed by BiFC experiments. This study reveals, for the first time, an overview of type-B RR dimerization in poplar and makes way for the hypothesis that signal specificity for cytokinin or osmotic stress could be in part due to the fact that it is impossible for specific type-B RRs to heterodimerize.

Adaptive gene misregulation

Because gene expression is important for evolutionary adaptation, its misregulation is an important cause of maladaptation. A misregulated gene can be incorrectly silent ("off") when a transcription factor (TF) that is required for its activation does not binds its regulatory region. Conversely, a misregulated gene can be incorrectly active ("on") when a TF not normally involved in its activation binds its regulatory region, a phenomenon also known as regulatory crosstalk. DNA mutations that destroy or create TF binding sites on DNA are an important source of misregulation and crosstalk. Although misregulation reduces fitness in an environment to which an organism is well-adapted, it may become adaptive in a new environment. Here, I derive simple yet general mathematical expressions that delimit the conditions under which misregulation can be adaptive. These expressions depend on the strength of selection against misregulation, on the fraction of DNA sequence space filled with TF binding sites, and on the fraction of genes that must be expressed for optimal adaptation. I then use empirical data from RNA sequencing, protein-binding microarrays, and genome evolution, together with population genetic simulations to ask when these conditions are likely to be met. I show that they can be met under realistic circumstances, but these circumstances may vary among organisms and environments. My analysis provides a framework in which improved theory and data collection can help us demonstrate the role of misregulation in adaptation. It also shows that misregulation, like DNA mutation, is one of life's many imperfections that can help propel Darwinian evolution.

Programmable protein circuit design

A fundamental challenge in synthetic biology is to create molecular circuits that can program complex cellular functions. Because proteins can bind, cleave, and chemically modify one another and interface directly and rapidly with endogenous pathways, they could extend the capabilities of synthetic circuits beyond what is possible with gene regulation alone. However, the very diversity that makes proteins so powerful also complicates efforts to harness them as well-controlled synthetic circuit components. Recent work has begun to address this challenge, focusing on principles such as orthogonality and composability that permit construction of diverse circuit-level functions from a limited set of engineered protein components. These approaches are now enabling the engineering of circuits that can sense, transmit, and process information; dynamically control cellular behaviors; and enable new therapeutic strategies, establishing a powerful paradigm for programming biology.

Uncovering the basis of protein-protein interaction specificity with a combinatorially complete library

Protein-protein interaction specificity is often encoded at the primary sequence level. However, the contributions of individual residues to specificity are usually poorly understood and often obscured by mutational robustness, sequence degeneracy, and epistasis. Using bacterial toxin-antitoxin systems as a model, we screened a combinatorially complete library of antitoxin variants at three key positions against two toxins. This library enabled us to measure the effect of individual substitutions on specificity in hundreds of genetic backgrounds. These distributions allow inferences about the general nature of interface residues in promoting specificity. We find that positive and negative contributions to specificity are neither inherently coupled nor mutually exclusive. Further, a wild-type antitoxin appears optimized for specificity as no substitutions improve discrimination between cognate and non-cognate partners. By comparing crystal structures of paralogous complexes, we provide a rationale for our observations. Collectively, this work provides a generalizable approach to understanding the logic of molecular recognition.