Marginal specificity in protein interactions constrains evolution of a paralogous family

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.

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.

Intrinsic Disorder and Other Malleable Arsenals of Evolved Protein Multifunctionality

Microscopic evolution at the functional biomolecular level is an ongoing process. Leveraging functional and high-throughput assays, along with computational data mining, has led to a remarkable expansion of our understanding of multifunctional protein (and gene) families over the past few decades. Various molecular and intermolecular mechanisms are now known that collectively meet the cumulative multifunctional demands in higher organisms along an evolutionary path. This multitasking ability is attributed to a certain degree of intrinsic or adapted flexibility at the structure-function level. Evolutionary diversification of structure-function relationships in proteins highlights the functional importance of intrinsically disordered proteins/regions (IDPs/IDRs) which are highly dynamic biological soft matter. Multifunctionality is favorably supported by the fluid-like shapes of IDPs/IDRs, enabling them to undergo disorder-to-order transitions upon binding to different molecular partners. Other new malleable members of the protein superfamily, such as those involved in fold-switching, also undergo structural transitions. This new insight diverges from all traditional notions of functional singularity in enzyme classes and emphasizes a far more complex, multi-layered diversification of protein functionality. However, a thorough review in this line, focusing on flexibility and function-driven structural transitions related to evolved multifunctionality in proteins, is currently missing. This review attempts to address this gap while broadening the scope of multifunctionality beyond single protein sequences. It argues that protein intrinsic disorder is likely the most striking mechanism for expressing multifunctionality in proteins. A phenomenological analogy has also been drawn to illustrate the increasingly complex nature of modern digital life, driven by the need for multitasking, particularly involving media.

Constitutive activation of two-component systems reveals regulatory network interactions in Streptococcus agalactiae

Bacterial two-component systems (TCSs) are signaling modules that control physiology, adaptation, and host interactions. A typical TCS consists of a histidine kinase (HK) that activates a response regulator via phosphorylation in response to environmental signals. Here, we systematically test the effect of inactivating the conserved phosphatase activity of HKs to activate TCS signaling pathways. Transcriptome analyses of 14 HK mutants in Streptococcus agalactiae, the leading cause of neonatal meningitis, validate the conserved HK phosphatase mechanism and its role in the inhibition of TCS activity in vivo. Constitutive TCS activation, independent of environmental signals, enables high-resolution mapping of the regulons for several TCSs (e.g., SaeRS, BceRS, VncRS, DltRS, HK11030, HK02290) and reveals the functional diversity of TCS signaling pathways, ranging from highly specialized to interconnected global regulatory networks. Targeted analysis shows that the SaeRS-regulated PbsP adhesin acts as a signaling molecule to activate CovRS signaling, thereby linking the major regulators of host-pathogen interactions. Furthermore, constitutive BceRS activation reveals drug-independent activity, suggesting a role in cell envelope homeostasis beyond antimicrobial resistance. This study highlights the versatility of constitutive TCS activation, via phosphatase-deficient HKs, to uncover regulatory networks and biological processes.

Engineering two-component systems for advanced biosensing: From architecture to applications in biotechnology

Two-component systems (TCSs) are prevalent signaling pathways in bacteria. These systems mediate phosphotransfer between histidine kinase and a response regulator, facilitating responses to diverse physical, chemical, and biological stimuli. Advancements in synthetic and structural biology have repurposed TCSs for applications in monitoring heavy metals, disease-associated biomarkers, and the production of bioproducts. However, the utility of many TCS biosensors is hindered by undesired performance due to the lack of effective engineering methods. Here, we briefly discuss the architectures and regulatory mechanisms of TCSs. We also summarize the recent advancements in TCS engineering by experimental or computational-based methods to fine-tune the biosensor functional parameters, such as response curve and specificity. Engineered TCSs have great potential in the medical, environmental, and biorefinery fields, demonstrating a crucial role in a wide area of biotechnology.

Predicting the involvement of polyQ- and polyA in protein-protein interactions by their amino acid context

Homorepeats, specifically polyglutamine (polyQ) and polyalanine (polyA), are often implicated in protein-protein interactions (PPIs). So far, a method to predict the participation of homorepeats in protein interactions is lacking. We propose a machine learning approach to identify PPI-involved polyQ and polyA regions within the human proteome based on known interacting regions. Using the dataset of human homorepeats, we identified 157 polyQ and 745 polyA regions potentially involved in PPIs. Machine learning models, trained on amino acid context and homorepeat length, demonstrated high precision (0.90-0.98) but variable recall (0.42-0.85). Random forest outperformed other models (AUC polyQ = 0.686, AUC polyA = 0.732) using the positions surrounding the homorepeat -10 to +10. Integrating paralog information marginally improved predictions but was excluded for model simplicity. Further optimization revealed that for polyQ, using amino acid surrounding positions from -6 to +6 increased AUC to 0.715. For polyA, no improvement was found. Incorporating coiled coil overlap information enhanced polyA predictions (AUC = 0.745) but not polyQ. Finally, we applied these models to predict PPI involvement across all polyQ and polyA regions, identifying potential interactions. Case studies illustrated the method's predictive capacity, highlighting known interacting regions with high scores and elucidating potential false negatives.

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.

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.

Conserved patterns of sequence diversification provide insight into the evolution of two-component systems in Enterobacteriaceae

Two-component regulatory systems (TCSs) are a major mechanism used by bacteria to sense and respond to their environments. Many of the same TCSs are used by biologically diverse organisms with different regulatory needs, suggesting that the functions of TCS must evolve. To explore this topic, we analysed the amino acid sequence divergence patterns of a large set of broadly conserved TCS across different branches of Enterobacteriaceae, a family of Gram-negative bacteria that includes biomedically important genera such as Salmonella, Escherichia, Klebsiella and others. Our analysis revealed trends in how TCS sequences change across different proteins or functional domains of the TCS, and across different lineages. Based on these trends, we identified individual TCS that exhibit atypical evolutionary patterns. We observed that the relative extent to which the sequence of a given TCS varies across different lineages is generally well conserved, unveiling a hierarchy of TCS sequence conservation with EnvZ/OmpR as the most conserved TCS. We provide evidence that, for the most divergent of the TCS analysed, PmrA/PmrB, different alleles were horizontally acquired by different branches of this family, and that different PmrA/PmrB sequence variants have highly divergent signal-sensing domains. Collectively, this study sheds light on how TCS evolve, and serves as a compendium for how the sequences of the TCS in this family have diverged over the course of evolution.

Evolutionary history of the cytochrome P450s from Colletotrichum species and prediction of their putative functional roles during host-pathogen interactions

The genomes of species belonging to the genus Colletotrichum harbor a substantial number of cytochrome P450 monooxygenases (CYPs) encoded by a broad diversity of gene families. However, the biological role of their CYP complement (CYPome) has not been elucidated. Here, we investigated the putative evolutionary scenarios that occurred during the evolution of the CYPome belonging to the Colletotrichum Graminicola species complex (s.c.) and their biological implications. The study revealed that most of the CYPome gene families belonging to the Graminicola s.c. experienced gene contractions. The reductive evolution resulted in species restricted CYPs are predominant in each CYPome of members from the Graminicola s.c., whereas only 18 families are absolutely conserved among these species. However, members of CYP families displayed a notably different phylogenetic relationship at the tertiary structure level, suggesting a putative convergent evolution scenario. Most of the CYP enzymes of the Graminicola s.c. share redundant functions in secondary metabolite biosynthesis and xenobiotic metabolism. Hence, this current work suggests that the presence of a broad CYPome in the genus Colletotrichum plays a critical role in the optimization of the colonization capability and virulence.

A balancing act: investigations on the impact of altered signal sensitivity in bacterial quorum sensing

IMPORTANCE

Quorum sensing (QS) is a widespread form of cell-cell signaling that regulates group behaviors important for competition and cooperation within bacterial communities. The QS systems from different bacterial species have diverse properties, but the functional consequences of this diversity are largely unknown. Taking advantage of hyper- and hypo-sensitive QS receptor variants in the opportunistic pathogen Pseudomonas aeruginosa, we examine the costs and benefits of altered signal sensitivity. We find that the sensitivity of a model QS receptor, LasR, impacts the timing and level of quorum gene expression, and fitness during intra- and interspecies competition. These findings suggest competition with kin and with other bacterial species work together to tune signal sensitivity.

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.

Evolutionary innovation through transcription factor rewiring in microbes is shaped by levels of transcription factor activity, expression, and existing connectivity

The survival of a population during environmental shifts depends on whether the rate of phenotypic adaptation keeps up with the rate of changing conditions. A common way to achieve this is via change to gene regulatory network (GRN) connections-known as rewiring-that facilitate novel interactions and innovation of transcription factors. To understand the success of rapidly adapting organisms, we therefore need to determine the rules that create and constrain opportunities for GRN rewiring. Here, using an experimental microbial model system with the soil bacterium Pseudomonas fluorescens, we reveal a hierarchy among transcription factors that are rewired to rescue lost function, with alternative rewiring pathways only unmasked after the preferred pathway is eliminated. We identify 3 key properties-high activation, high expression, and preexisting low-level affinity for novel target genes-that facilitate transcription factor innovation. Ease of acquiring these properties is constrained by preexisting GRN architecture, which was overcome in our experimental system by both targeted and global network alterations. This work reveals the key properties that determine transcription factor evolvability, and as such, the evolution of GRNs.