Engineering orthogonal signalling pathways reveals the sparse occupancy of sequence space

The genetic architecture of protein interaction affinity and specificity

The encoding and evolution of specificity and affinity in protein-protein interactions is poorly understood. Here, we address this question by quantifying how all mutations in one protein, JUN, alter binding to all other members of a protein family, the 54 human basic leucine zipper transcription factors. We fit a global thermodynamic model to the data to reveal that most affinity changing mutations equally affect JUN's affinity to all its interaction partners. Mutations that alter binding specificity are relatively rare but distributed throughout the interaction interface. Specificity is determined both by features that promote on-target interactions and by those that prevent off-target interactions. Approximately half of the specificity-defining residues in JUN contribute both to promoting on-target binding and preventing off-target binding. Nearly all specificity-altering mutations in the interaction interface are pleiotropic, also altering affinity to all partners. In contrast, mutations outside the interface can tune global affinity without affecting specificity. Our results reveal the distributed encoding of specificity and affinity in an interaction interface and how coiled-coils provide an elegant solution to the challenge of optimizing both specificity and affinity in a large protein 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. 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.

Deep-learning-based design of synthetic orthologs of SH3 signaling domains

Evolution-based deep generative models represent an exciting direction in understanding and designing proteins. An open question is whether such models can learn specialized functional constraints that control fitness in specific biological contexts. Here, we examine the ability of generative models to produce synthetic versions of Src-homology 3 (SH3) domains that mediate signaling in the Sho1 osmotic stress response pathway of yeast. We show that a variational autoencoder (VAE) model produces artificial sequences that experimentally recapitulate the function of natural SH3 domains. More generally, the model organizes all fungal SH3 domains such that locality in the model latent space (but not simply locality in sequence space) enriches the design of synthetic orthologs and exposes non-obvious amino acid constraints distributed near and far from the SH3 ligand-binding site. The ability of generative models to design ortholog-like functions in vivo opens new avenues for engineering protein function in specific cellular contexts and environments.

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.

Engineering Chimeric Chemoreceptors and Two-Component Systems for Orthogonal and Leakless Biosensing of Extracellular γ-Aminobutyric Acid

Two-component systems (TCSs) sensing and responding to various stimuli outside and inside cells are valuable resources for developing biosensors with synthetic biology applications. However, the use of TCS-based biosensors suffers from a limited effector spectrum, hypersensitivity, low dynamic range, and unwanted signal crosstalk. Here, we developed a tailor-made Escherichia coli whole-cell γ-aminobutyric acid (GABA) biosensor by engineering a chimeric GABA chemoreceptor PctC and TCS. By testing different TCSs, the chimeric PctC/PhoQ showed the response to GABA. Chimera-directed evolution and introduction of the insulated chimeric pair PctC/PhoQ*PhoP* produced biosensors with up to 3.50-fold dynamic range and good orthogonality. To further enhance the dynamic range and lower the basal leakage, three strategies, engineering of PhoP DNA binding sites, fine-tuning reporter expression by optimizing transcription/translation components, and a tobacco etch virus protease-controlled protein degradation, were integrated. This chimeric biosensor displayed a low basal leakage, a large dynamic range (15.8-fold), and a high threshold level (22.7 g L-1). Finally, the optimized biosensor was successfully applied in the high-throughput microdroplet screening of GABA-overproducing Corynebacterium glutamicum, demonstrating its desired properties for extracellular signal biosensing.

Accelerating Genetic Sensor Development, Scale-up, and Deployment Using Synthetic Biology

Living cells are exquisitely tuned to sense and respond to changes in their environment. Repurposing these systems to create engineered biosensors has seen growing interest in the field of synthetic biology and provides a foundation for many innovative applications spanning environmental monitoring to improved biobased production. In this review, we present a detailed overview of currently available biosensors and the methods that have supported their development, scale-up, and deployment. We focus on genetic sensors in living cells whose outputs affect gene expression. We find that emerging high-throughput experimental assays and evolutionary approaches combined with advanced bioinformatics and machine learning are establishing pipelines to produce genetic sensors for virtually any small molecule, protein, or nucleic acid. However, more complex sensing tasks based on classifying compositions of many stimuli and the reliable deployment of these systems into real-world settings remain challenges. We suggest that recent advances in our ability to precisely modify nonmodel organisms and the integration of proven control engineering principles (e.g., feedback) into the broader design of genetic sensing systems will be necessary to overcome these hurdles and realize the immense potential of the field.

Multiscale regulation of nutrient stress responses in Escherichia coli from chromatin structure to small regulatory RNAs

Recent research has indicated the presence of heterochromatin-like regions of extended protein occupancy and transcriptional silencing of bacterial genomes. We utilized an integrative approach to track chromatin structure and transcription in E. coli K-12 across a wide range of nutrient conditions. In the process, we identified multiple loci which act similarly to facultative heterochromatin in eukaryotes, normally silenced but permitting expression of genes under specific conditions. We also found a strong enrichment of small regulatory RNAs (sRNAs) among the set of differentially expressed transcripts during nutrient stress. Using a newly developed bioinformatic pipeline, the transcription factors regulating sRNA expression were bioinformatically predicted, with experimental follow-up revealing novel relationships for 36 sRNA-transcription factors candidates. Direct regulation of sRNA expression was confirmed by mutational analysis for five sRNAs of metabolic interest: IsrB, CsrB and CsrC, GcvB, and GadY. Our integrative analysis thus reveals additional layers of complexity in the nutrient stress response in E. coli and provides a framework for revealing similar poorly understood regulatory logic in other organisms.

Epistasis facilitates functional evolution in an ancient transcription factor

A protein's genetic architecture - the set of causal rules by which its sequence produces its functions - also determines its possible evolutionary trajectories. Prior research has proposed that the genetic architecture of proteins is very complex, with pervasive epistatic interactions that constrain evolution and make function difficult to predict from sequence. Most of this work has analyzed only the direct paths between two proteins of interest - excluding the vast majority of possible genotypes and evolutionary trajectories - and has considered only a single protein function, leaving unaddressed the genetic architecture of functional specificity and its impact on the evolution of new functions. Here, we develop a new method based on ordinal logistic regression to directly characterize the global genetic determinants of multiple protein functions from 20-state combinatorial deep mutational scanning (DMS) experiments. We use it to dissect the genetic architecture and evolution of a transcription factor's specificity for DNA, using data from a combinatorial DMS of an ancient steroid hormone receptor's capacity to activate transcription from two biologically relevant DNA elements. We show that the genetic architecture of DNA recognition consists of a dense set of main and pairwise effects that involve virtually every possible amino acid state in the protein-DNA interface, but higher-order epistasis plays only a tiny role. Pairwise interactions enlarge the set of functional sequences and are the primary determinants of specificity for different DNA elements. They also massively expand the number of opportunities for single-residue mutations to switch specificity from one DNA target to another. By bringing variants with different functions close together in sequence space, pairwise epistasis therefore facilitates rather than constrains the evolution of new functions.

The genetic landscape of a metabolic interaction

While much prior work has explored the constraints on protein sequence and evolution induced by physical protein-protein interactions, the sequence-level constraints emerging from non-binding functional interactions in metabolism remain unclear. To quantify how variation in the activity of one enzyme constrains the biochemical parameters and sequence of another, we focus on dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), a pair of enzymes catalyzing consecutive reactions in folate metabolism. We use deep mutational scanning to quantify the growth rate effect of 2696 DHFR single mutations in 3 TYMS backgrounds under conditions selected to emphasize biochemical epistasis. Our data are well-described by a relatively simple enzyme velocity to growth rate model that quantifies how metabolic context tunes enzyme mutational tolerance. Together our results reveal the structural distribution of epistasis in a metabolic enzyme and establish a foundation for the design of multi-enzyme systems.

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.

REC protein family expansion by the emergence of a new signaling pathway

IMPORTANCE

We explore when and why large classes of proteins expand into new sequence space. We used an unsupervised machine learning approach to observe the sequence landscape of REC domains of bacterial response regulator proteins. We find that within-gene recombination can switch effector domains and, consequently, change the regulatory context of the duplicated protein.

Customizing cellular signal processing by synthetic multi-level regulatory circuits

As synthetic biology permeates society, the signal processing circuits in engineered living systems must be customized to meet practical demands. Towards this mission, novel regulatory mechanisms and genetic circuits with unprecedented complexity have been implemented over the past decade. These regulatory mechanisms, such as transcription and translation control, could be integrated into hybrid circuits termed "multi-level circuits". The multi-level circuit design will tremendously benefit the current genetic circuit design paradigm, from modifying basic circuit dynamics to facilitating real-world applications, unleashing our capabilities to customize cellular signal processing and address global challenges through synthetic biology.

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 Genetic Landscape of a Metabolic Interaction

Enzyme abundance, catalytic activity, and ultimately sequence are all shaped by the need of growing cells to maintain metabolic flux while minimizing accumulation of deleterious intermediates. While much prior work has explored the constraints on protein sequence and evolution induced by physical protein-protein interactions, the sequence-level constraints emerging from non-binding functional interactions in metabolism remain unclear. To quantify how variation in the activity of one enzyme constrains the biochemical parameters and sequence of another, we focused on dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), a pair of enzymes catalyzing consecutive reactions in folate metabolism. We used deep mutational scanning to quantify the growth rate effect of 2,696 DHFR single mutations in 3 TYMS backgrounds under conditions selected to emphasize biochemical epistasis. Our data are well-described by a relatively simple enzyme velocity to growth rate model that quantifies how metabolic context tunes enzyme mutational tolerance. Together our results reveal the structural distribution of epistasis in a metabolic enzyme and establish a foundation for the design of multi-enzyme systems.

Protein engineering via sequence-performance mapping

Discovery and evolution of new and improved proteins has empowered molecular therapeutics, diagnostics, and industrial biotechnology. Discovery and evolution both require efficient screens and effective libraries, although they differ in their challenges because of the absence or presence, respectively, of an initial protein variant with the desired function. A host of high-throughput technologies-experimental and computational-enable efficient screens to identify performant protein variants. In partnership, an informed search of sequence space is needed to overcome the immensity, sparsity, and complexity of the sequence-performance landscape. Early in the historical trajectory of protein engineering, these elements aligned with distinct approaches to identify the most performant sequence: selection from large, randomized combinatorial libraries versus rational computational design. Substantial advances have now emerged from the synergy of these perspectives. Rational design of combinatorial libraries aids the experimental search of sequence space, and high-throughput, high-integrity experimental data inform computational design. At the core of the collaborative interface, efficient protein characterization (rather than mere selection of optimal variants) maps sequence-performance landscapes. Such quantitative maps elucidate the complex relationships between protein sequence and performance-e.g., binding, catalytic efficiency, biological activity, and developability-thereby advancing fundamental protein science and facilitating protein discovery and evolution.

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.

Data-driven design of orthogonal protein-protein interactions

Engineering protein-protein interactions to generate new functions presents a challenge with great potential for many applications, ranging from therapeutics to synthetic biology. To avoid unwanted cross-talk with preexisting protein interaction networks in a cell, the specificity and selectivity of newly engineered proteins must be controlled. Here, we developed a computational strategy that mimics gene duplication and the divergence of preexisting interacting protein pairs to design new interactions. We used the bacterial PhoQ-PhoP two-component system as a model system to demonstrate the feasibility of this strategy and validated the approach with known experimental results. The designed protein pairs are predicted to exclusively interact with each other and to be insulated from potential cross-talk with their native partners. Thus, our approach enables exploration of uncharted regions of the protein sequence space and the design of new interacting protein pairs.

Structural features of sensory two component systems: a synthetic biology perspective

All living organisms include a set of signaling devices that confer the ability to dynamically perceive and adapt to the fluctuating environment. Two-component systems are part of this sensory machinery that regulates the execution of different genetic and/or biochemical programs in response to specific physical or chemical signals. In the last two decades, there has been tremendous progress in our molecular understanding on how signals are detected, the allosteric mechanisms that control intramolecular information transmission and the specificity determinants that guarantee correct wiring. All this information is starting to be exploited in the development of new synthetic networks. Connecting multiple molecular players, analogous to programming lines of code, can provide the resources to build new sophisticated biocomputing systems. The Synthetic Biology field is starting to revolutionize several scientific fields, such as biomedicine and agriculture, propelling the development of new solutions. Expanding the spectrum of available nanodevices in the toolbox is key to unleash its full potential. This review aims to discuss, from a structural perspective, how to take advantage of the vast array of sensor and effector protein modules involved in two-component systems for the construction of new synthetic circuits.

An Evolutionary Paradigm Favoring Cross Talk between Bacterial Two-Component Signaling Systems

The prevalent paradigm governing bacterial two-component signaling systems (TCSs) is specificity, wherein the histidine kinase (HK) of a TCS exclusively activates its cognate response regulator (RR). Cross talk, where HKs activate noncognate RRs, is considered evolutionarily disadvantageous because it can compromise adaptive responses by leaking signals. Yet cross talk is observed in several bacteria. Here, to resolve this paradox, we propose an alternative paradigm where cross talk can be advantageous. We envisioned programmed environments, wherein signals appear in predefined sequences. In such environments, cross talk that primes bacteria to upcoming signals may improve adaptive responses and confer evolutionary benefits. To test this hypothesis, we employed mathematical modeling of TCS signaling networks and stochastic evolutionary dynamics simulations. We considered the comprehensive set of bacterial phenotypes, comprising thousands of distinct cross talk patterns competing in varied signaling environments. Our simulations predicted that in programmed environments phenotypes with cross talk facilitating priming would outcompete phenotypes without cross talk. In environments where signals appear randomly, bacteria without cross talk would dominate, explaining the specificity widely seen. Additionally, a testable prediction was that the phenotypes selected in programmed environments would display one-way cross talk, ensuring priming to future signals. Interestingly, the cross talk networks we deduced from available data on TCSs of Mycobacterium tuberculosis all displayed one-way cross talk, which was consistent with our predictions. Our study thus identifies potential evolutionary underpinnings of cross talk in bacterial TCSs, suggests a reconciliation of specificity and cross talk, makes testable predictions of the nature of cross talk patterns selected, and has implications for understanding bacterial adaptation and the response to interventions. IMPORTANCE Bacteria use two-component signaling systems (TCSs) to sense and respond to environmental changes. The prevalent paradigm governing TCSs is specificity, where signal flow through TCSs is insulated; leakage to other TCSs is considered evolutionarily disadvantageous. Yet cross talk between TCSs is observed in many bacteria. Here, we present a potential resolution of this paradox. We envision programmed environments, wherein stimuli appear in predefined sequences. Cross talk that primes bacteria to upcoming stimuli could then confer evolutionary benefits. We demonstrate this benefit using mathematical modeling and evolutionary simulations. Interestingly, we found signatures of predicted cross talk patterns in Mycobacterium tuberculosis. Furthermore, specificity was selected in environments where stimuli occurred randomly, thus reconciling specificity and cross talk. Implications follow for understanding bacterial evolution and for interventions.

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.

The role of sensory kinase proteins in two-component signal transduction

Two-component systems (TCSs) are modular signaling circuits that regulate diverse aspects of microbial physiology in response to environmental cues. These molecular circuits comprise a sensor histidine kinase (HK) protein that contains a conserved histidine residue, and an effector response regulator (RR) protein with a conserved aspartate residue. HKs play a major role in bacterial signaling, since they perceive specific stimuli, transmit the message across the cytoplasmic membrane, and catalyze their own phosphorylation, and the trans-phosphorylation and dephosphorylation of their cognate response regulator. The molecular mechanisms by which HKs co-ordinate these functions have been extensively analyzed by genetic, biochemical, and structural approaches. Here, we describe the most common modular architectures found in bacterial HKs, and address the operation mode of the individual functional domains. Finally, we discuss the use of these signaling proteins as drug targets or as sensing devices in whole-cell biosensors with medical and biotechnological applications.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Directed evolution of colE1 plasmid replication compatibility: a fast tractable tunable model for investigating biological orthogonality

Plasmids of the ColE1 family are among the most frequently used in molecular biology. They were adopted early for many biotechnology applications, and as models to study plasmid biology. Their mechanism of replication is well understood, involving specific interactions between a plasmid encoded sense-antisense gene pair (RNAI and RNAII). Due to such mechanism, two plasmids with the same origin cannot be stably maintained in cells-a process known as incompatibility. While mutations in RNAI and RNAII can make colE1 more compatible, there has been no systematic effort to engineer new compatible colE1 origins, which could bypass technical design constraints for multi-plasmid applications. Here, we show that by diversifying loop regions in RNAI (and RNAII), it is possible to select new viable colE1 origins compatible with the wild-type one. We demonstrate that sequence divergence is not sufficient to enable compatibility and pairwise interactions are not an accurate guide for higher order interactions. We identify potential principles to engineer plasmid copy number independently from other regulatory strategies and we propose plasmid compatibility as a tractable model to study biological orthogonality.

Therapeutic cell engineering: designing programmable synthetic genetic circuits in mammalian cells

Cell therapy approaches that employ engineered mammalian cells for on-demand production of therapeutic agents in the patient's body are moving beyond proof-of-concept in translational medicine. The therapeutic cells can be customized to sense user-defined signals, process them, and respond in a programmable and predictable way. In this paper, we introduce the available tools and strategies employed to design therapeutic cells. Then, various approaches to control cell behaviors, including open-loop and closed-loop systems, are discussed. We also highlight therapeutic applications of engineered cells for early diagnosis and treatment of various diseases in the clinic and in experimental disease models. Finally, we consider emerging technologies such as digital devices and their potential for incorporation into future cell-based therapies.

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

Gene duplication is crucial to generating novel signaling pathways during evolution. However, it remains unclear how the redundant proteins produced by gene duplication ultimately acquire new interaction specificities to establish insulated paralogous signaling pathways. Here, we used ancestral sequence reconstruction to resurrect and characterize a bacterial two-component signaling system that duplicated in α-proteobacteria. We determined the interaction specificities of the signaling proteins that existed before and immediately after this duplication event and then identified key mutations responsible for establishing specificity in the two systems. Just three mutations, in only two of the four interacting proteins, were sufficient to establish specificity of the extant systems. Some of these mutations weakened interactions between paralogous systems to limit crosstalk. However, others strengthened interactions within a system, indicating that the ancestral interaction, although functional, had the potential to be strengthened. Our work suggests that protein-protein interactions with such latent potential may be highly amenable to duplication and divergence.

Co-evolution of interacting proteins through non-contacting and non-specific mutations

Proteins often accumulate neutral mutations that do not affect current functions but can profoundly influence future mutational possibilities and functions. Understanding such hidden potential has major implications for protein design and evolutionary forecasting but has been limited by a lack of systematic efforts to identify potentiating mutations. Here, through the comprehensive analysis of a bacterial toxin-antitoxin system, we identified all possible single substitutions in the toxin that enable it to tolerate otherwise interface-disrupting mutations in its antitoxin. Strikingly, the majority of enabling mutations in the toxin do not contact and promote tolerance non-specifically to many different antitoxin mutations, despite covariation in homologues occurring primarily between specific pairs of contacting residues across the interface. In addition, the enabling mutations we identified expand future mutational paths that both maintain old toxin-antitoxin interactions and form new ones. These non-specific mutations are missed by widely used covariation and machine learning methods. Identifying such enabling mutations will be critical for ensuring continued binding of therapeutically relevant proteins, such as antibodies, aimed at evolving targets.

Methods for the directed evolution of biomolecular interactions

Noncovalent interactions between biomolecules such as proteins and nucleic acids coordinate all cellular processes through changes in proximity. Tools that perturb these interactions are and will continue to be highly valuable for basic and translational scientific endeavors. By taking cues from natural systems, such as the adaptive immune system, we can design directed evolution platforms that can generate proteins that bind to biomolecules of interest. In recent years, the platforms used to direct the evolution of biomolecular binders have greatly expanded the range of types of interactions one can evolve. Herein, we review recent advances in methods to evolve protein-protein, protein-RNA, and protein-DNA interactions.

Robust and tunable signal processing in mammalian cells via engineered covalent modification cycles

Engineered signaling networks can impart cells with new functionalities useful for directing differentiation and actuating cellular therapies. For such applications, the engineered networks must be tunable, precisely regulate target gene expression, and be robust to perturbations within the complex context of mammalian cells. Here, we use bacterial two-component signaling proteins to develop synthetic phosphoregulation devices that exhibit these properties in mammalian cells. First, we engineer a synthetic covalent modification cycle based on kinase and phosphatase proteins derived from the bifunctional histidine kinase EnvZ, enabling analog tuning of gene expression via its response regulator OmpR. By regulating phosphatase expression with endogenous miRNAs, we demonstrate cell-type specific signaling responses and a new strategy for accurate cell type classification. Finally, we implement a tunable negative feedback controller via a small molecule-stabilized phosphatase, reducing output expression variance and mitigating the context-dependent effects of off-target regulation and resource competition. Our work lays the foundation for establishing tunable, precise, and robust control over cell behavior with synthetic signaling networks.

Droplet-based screening of phosphate transfer catalysis reveals how epistasis shapes MAP kinase interactions with substrates

The combination of ultrahigh-throughput screening and sequencing informs on function and intragenic epistasis within combinatorial protein mutant libraries. Establishing a droplet-based, in vitro compartmentalised approach for robust expression and screening of protein kinase cascades (>10⁷ variants/day) allowed us to dissect the intrinsic molecular features of the MKK-ERK signalling pathway, without interference from endogenous cellular components. In a six-residue combinatorial library of the MKK1 docking domain, we identified 29,563 sequence permutations that allow MKK1 to efficiently phosphorylate and activate its downstream target kinase ERK2. A flexibly placed hydrophobic sequence motif emerges which is defined by higher order epistatic interactions between six residues, suggesting synergy that enables high connectivity in the sequence landscape. Through positive epistasis, MKK1 maintains function during mutagenesis, establishing the importance of co-dependent residues in mammalian protein kinase-substrate interactions, and creating a scenario for the evolution of diverse human signalling networks.

Here, the authors use a droplet-based screen for phosphate transfer catalysis, testing variants of the human protein kinase MKK1 for its ability to activate its downstream target ERK2. Data reveal a flexible motif in the MKK1 docking domain that promotes efficient activation of ERK2, and suggest epistasis between the residues within that sequence.

Bacterial two-component systems as sensors for synthetic biology applications

Two-component systems (TCSs) are a ubiquitous family of signal transduction pathways that enable bacteria to sense and respond to diverse physical, chemical, and biological stimuli outside and inside the cell. Synthetic biologists have begun to repurpose TCSs for applications in optogenetics, materials science, gut microbiome engineering, and soil nutrient biosensing, among others. New engineering methods including genetic refactoring, DNA-binding domain swapping, detection threshold tuning, and phosphorylation cross-talk insulation are being used to increase the reliability of TCS sensor performance and tailor TCS signaling properties to the requirements of specific applications. There is now potential to combine these methods with large-scale gene synthesis and laboratory screening to discover the inputs sensed by many uncharacterized TCSs and develop a large new family of genetically-encoded sensors that respond to an unrivaled breadth of stimuli.

Harnessing intercellular signals to engineer the soil microbiome

Covering: Focus on 2015 to 2020Plant and soil microbiomes consist of diverse communities of organisms from across kingdoms and can profoundly affect plant growth and health. Natural product-based intercellular signals govern important interactions between microbiome members that ultimately regulate their beneficial or harmful impacts on the plant. Exploiting these evolved signalling circuits to engineer microbiomes towards beneficial interactions with crops is an attractive goal. There are few reports thus far of engineering the intercellular signalling of microbiomes, but this article argues that it represents a tremendous opportunity for advancing the field of microbiome engineering. This could be achieved through the selection of synergistic consortia in combination with genetic engineering of signal pathways to realise an optimised microbiome.

D-Xylose Sensing in Saccharomyces cerevisiae: Insights from D-Glucose Signaling and Native D-Xylose Utilizers

Extension of the substrate range is among one of the metabolic engineering goals for microorganisms used in biotechnological processes because it enables the use of a wide range of raw materials as substrates. One of the most prominent examples is the engineering of baker's yeast Saccharomyces cerevisiae for the utilization of d-xylose, a five-carbon sugar found in high abundance in lignocellulosic biomass and a key substrate to achieve good process economy in chemical production from renewable and non-edible plant feedstocks. Despite many excellent engineering strategies that have allowed recombinant S. cerevisiae to ferment d-xylose to ethanol at high yields, the consumption rate of d-xylose is still significantly lower than that of its preferred sugar d-glucose. In mixed d-glucose/d-xylose cultivations, d-xylose is only utilized after d-glucose depletion, which leads to prolonged process times and added costs. Due to this limitation, the response on d-xylose in the native sugar signaling pathways has emerged as a promising next-level engineering target. Here we review the current status of the knowledge of the response of S. cerevisiae signaling pathways to d-xylose. To do this, we first summarize the response of the native sensing and signaling pathways in S. cerevisiae to d-glucose (the preferred sugar of the yeast). Using the d-glucose case as a point of reference, we then proceed to discuss the known signaling response to d-xylose in S. cerevisiae and current attempts of improving the response by signaling engineering using native targets and synthetic (non-native) regulatory circuits.

Programmable receptors enable bacterial biosensors to detect pathological biomarkers in clinical samples

Bacterial biosensors, or bactosensors, are promising agents for medical and environmental diagnostics. However, the lack of scalable frameworks to systematically program ligand detection limits their applications. Here we show how novel, clinically relevant sensing modalities can be introduced into bactosensors in a modular fashion. To do so, we have leveraged a synthetic receptor platform, termed EMeRALD (Engineered Modularized Receptors Activated via Ligand-induced Dimerization) which supports the modular assembly of sensing modules onto a high-performance, generic signaling scaffold controlling gene expression in E. coli. We apply EMeRALD to detect bile salts, a biomarker of liver dysfunction, by repurposing sensing modules from enteropathogenic Vibrio species. We improve the sensitivity and lower the limit-of-detection of the sensing module by directed evolution. We then engineer a colorimetric bactosensor detecting pathological bile salt levels in serum from patients having undergone liver transplant, providing an output detectable by the naked-eye. The EMeRALD technology enables functional exploration of natural sensing modules and rapid engineering of synthetic receptors for diagnostics, environmental monitoring, and control of therapeutic microbes.

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.

Constraints on the expansion of paralogous protein families

Duplication and divergence is a major mechanism by which new proteins and functions emerge in biology. Consequently, most organisms, in all domains of life, have genomes that encode large paralogous families of proteins. For recently duplicated pathways to acquire different, independent functions, the two paralogs must acquire mutations that effectively insulate them from one another. For instance, paralogous signaling proteins must acquire mutations that endow them with different interaction specificities such that they can participate in different signaling pathways without disruptive cross-talk. Although duplicated genes undoubtedly shape each other's evolution as they diverge and attain new functions, it is less clear how other paralogs impact or constrain gene duplication. Does the establishment of a new pathway by duplication and divergence require the system-wide optimization of all paralogs? The answer has profound implications for molecular evolution and our ability to engineer biological systems. Here, we discuss models, experiments, and approaches for tackling this question, and for understanding how new proteins and pathways are born.

Uncovering Transcriptional Regulators and Targets of sRNAs Using an Integrative Data-Mining Approach: H-NS-Regulated RseX as a Case Study

Bacterial small RNAs (sRNAs) play a vital role in pathogenesis by enabling rapid, efficient networks of gene attenuation during infection. In recent decades, there has been a surge in the number of proposed and biochemically-confirmed sRNAs in both Gram-positive and Gram-negative pathogens. However, limited homology, network complexity, and condition specificity of sRNA has stunted complete characterization of the activity and regulation of these RNA regulators. To streamline the discovery of the expression of sRNAs, and their post-transcriptional activities, we propose an integrative in vivo data-mining approach that couples DNA protein occupancy, RNA-seq, and RNA accessibility data with motif identification and target prediction algorithms. We benchmark the approach against a subset of well-characterized E. coli sRNAs for which a degree of in vivo transcriptional regulation and post-transcriptional activity has been previously reported, finding support for known regulation in a large proportion of this sRNA set. We showcase the abilities of our method to expand understanding of sRNA RseX, a known envelope stress-linked sRNA for which a cellular role has been elusive due to a lack of native expression detection. Using the presented approach, we identify a small set of putative RseX regulators and targets for experimental investigation. These findings have allowed us to confirm native RseX expression under conditions that eliminate H-NS repression as well as uncover a post-transcriptional role of RseX in fimbrial regulation. Beyond RseX, we uncover 163 putative regulatory DNA-binding protein sites, corresponding to regulation of 62 sRNAs, that could lead to new understanding of sRNA transcription regulation. For 32 sRNAs, we also propose a subset of top targets filtered by engagement of regions that exhibit binding site accessibility behavior in vivo. We broadly anticipate that the proposed approach will be useful for sRNA-reliant network characterization in bacteria. Such investigations under pathogenesis-relevant environmental conditions will enable us to deduce complex rapid-regulation schemes that support infection.

Directed evolution of orthogonal RNA-RBP pairs through library-vs-library in vitro selection

RNA-binding proteins (RBPs) and their RNA ligands play many critical roles in gene regulation and RNA processing in cells. They are also useful for various applications in cell biology and synthetic biology. However, re-engineering novel and orthogonal RNA-RBP pairs from natural components remains challenging while such synthetic RNA-RBP pairs could significantly expand the RNA-RBP toolbox for various applications. Here, we report a novel library-vs-library in vitro selection strategy based on Phage Display coupled with Systematic Evolution of Ligands by EXponential enrichment (PD-SELEX). Starting with pools of 1.1 × 1012 unique RNA sequences and 4.0 × 108 unique phage-displayed L7Ae-scaffold (LS) proteins, we selected RNA-RBP complexes through a two-step affinity purification process. After six rounds of library-vs-library selection, the selected RNAs and LS proteins were analyzed by next-generation sequencing (NGS). Further deconvolution of the enriched RNA and LS protein sequences revealed two synthetic and orthogonal RNA-RBP pairs that exhibit picomolar affinity and >4000-fold selectivity.

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.

The upcycled roles of pseudoenzymes in two-component signal transduction

Upon first glance at a bacterial genome, pseudoenzymes appear unremarkable due to their lack of critical motifs that facilitate catalysis. These pseudoenzymes exist within signal transduction enzymes including histidine kinases, response regulators, diguanylate cyclases, and phosphodiesterases. Here, we summarize recent studies of bacterial pseudo-histidine kinases and pseudo-response regulators that regulate cell division, capsule formation, and the circadian rhythm. These examples illuminate the mechanistic potential of catalytically dead signaling enzymes and their impact upon bacterial signal transduction. Moreover, proteins lacking characteristic catalytic features of two-component systems reveal the sophisticated underlying potential of canonical two-component systems.

Designing Biological Circuits: Synthetic Biology Within the Operon Model and Beyond

In 1961, Jacob and Monod proposed the operon model of gene regulation. At the model's core was the modular assembly of regulators, operators, and structural genes. To illustrate the composability of these elements, Jacob and Monod linked phenotypic diversity to the architectures of regulatory circuits. In this review, we examine how the circuit blueprints imagined by Jacob and Monod laid the foundation for the first synthetic gene networks that launched the field of synthetic biology in 2000. We discuss the influences of the operon model and its broader theoretical framework on the first generation of synthetic biological circuits, which were predominantly transcriptional and posttranscriptional circuits. We also describe how recent advances in molecular biology beyond the operon model-namely, programmable DNA- and RNA-binding molecules as well as models of epigenetic and posttranslational regulation-are expanding the synthetic biology toolkit and enabling the design of more complex biological circuits.

Sensing the future of bio-informational engineering

The practices of synthetic biology are being integrated into 'multiscale' designs enabling two-way communication across organic and inorganic information substrates in biological, digital and cyber-physical system integrations. Novel applications of 'bio-informational' engineering will arise in environmental monitoring, precision agriculture, precision medicine and next-generation biomanufacturing. Potential developments include sentinel plants for environmental monitoring and autonomous bioreactors that respond to biosensor signaling. As bio-informational understanding progresses, both natural and engineered biological systems will need to be reimagined as cyber-physical architectures. We propose that a multiple length scale taxonomy will assist in rationalizing and enabling this transformative development in engineering biology.

Evidence for Pentapeptide-Dependent and Independent CheB Methylesterases

Many bacteria possess multiple chemosensory pathways that are composed of homologous signaling proteins. These pathways appear to be functionally insulated from each other, but little information is available on the corresponding molecular basis. We report here a novel mechanism that contributes to pathway insulation. We show that, of the four CheB paralogs of Pseudomonas aeruginosa PAO1, only CheB₂ recognizes a pentapeptide at the C-terminal extension of the McpB (Aer2) chemoreceptor (K_D = 93 µM). McpB is the sole chemoreceptor that stimulates the Che2 pathway, and CheB₂ is the methylesterase of this pathway. Pectobacterium atrosepticum SCRI1043 has a single CheB, CheB_Pec, and 19 of its 36 chemoreceptors contain a C-terminal pentapeptide. The deletion of cheB_Pec abolished chemotaxis, but, surprisingly, none of the pentapeptides bound to CheB_Pec. To determine the corresponding structural basis, we solved the 3D structure of CheB_Pec. Its structure aligned well with that of the pentapeptide-dependent enzyme from Salmonella enterica. However, no electron density was observed in the CheB_Pec region corresponding to the pentapeptide-binding site in the Escherichia coli CheB. We hypothesize that this structural disorder is associated with the failure to bind pentapeptides. Combined data show that CheB methylesterases can be divided into pentapeptide-dependent and independent enzymes.

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.

Engineered protein switches for exogenous control of gene expression

There is an ongoing need in the synthetic biology community for novel ways to regulate gene expression. Protein switches, which sense biological inputs and respond with functional outputs, represent one way to meet this need. Despite the fact that there is already a large pool of transcription factors and signaling proteins available, the pool of existing switches lacks the substrate specificities and activities required for certain applications. Therefore, a large number of techniques have been applied to engineer switches with novel properties. Here we discuss some of these techniques by broadly organizing them into three approaches. We show how novel switches can be created through mutagenesis, domain swapping, or domain insertion. We then briefly discuss their use as biosensors and in complex genetic circuits.

De novo design of an intercellular signaling toolbox for multi-channel cell-cell communication and biological computation

Intercellular signaling is indispensable for single cells to form complex biological structures, such as biofilms, tissues and organs. The genetic tools available for engineering intercellular signaling, however, are quite limited. Here we exploit the chemical diversity of biological small molecules to de novo design a genetic toolbox for high-performance, multi-channel cell-cell communications and biological computations. By biosynthetic pathway design for signal molecules, rational engineering of sensing promoters and directed evolution of sensing transcription factors, we obtain six cell-cell signaling channels in bacteria with orthogonality far exceeding the conventional quorum sensing systems and successfully transfer some of them into yeast and human cells. For demonstration, they are applied in cell consortia to generate bacterial colony-patterns using up to four signaling channels simultaneously and to implement distributed bio-computation containing seven different strains as basic units. This intercellular signaling toolbox paves the way for engineering complex multicellularity including artificial ecosystems and smart tissues.

Engineered bacteria to report gut function: technologies and implementation

Advances in synthetic biology and microbiology have enabled the creation of engineered bacteria which can sense and report on intracellular and extracellular signals. When deployed in vivo these whole-cell bacterial biosensors can act as sentinels to monitor biomolecules of interest in human health and disease settings. This is particularly interesting in the context of the gut microbiota, which interacts extensively with the human host throughout time and transit of the gut and can be accessed from feces without requiring invasive collection. Leveraging rational engineering approaches for genetic circuits as well as an expanding catalog of disease-associated biomarkers, bacterial biosensors can act as non-invasive and easy-to-monitor reporters of the gut. Here, we summarize recent engineering approaches applied in vivo in animal models and then highlight promising technologies for designing the next generation of bacterial biosensors.

Small design from big alignment: engineering proteins with multiple sequence alignment as the starting point

Multiple sequence alignment (MSA) is a fundamental way to gain information that cannot be obtained from the analysis of any individual sequence included in the alignment. It provides ways to investigate the relationship between sequence and function from a perspective of evolution. Thus, the MSA of proteins can be employed as a reference for protein engineering. In this paper, we reviewed the recent advances to highlight how protein engineering was benefited from the MSA of proteins. These methods include (1) engineering the thermostability or solubility of proteins by making it closer to the consensus sequence of the alignment through introducing site mutations; (2) structure-based engineering proteins with comparative modeling; (3) creating paleoenzymes featured with high thermostability and promiscuity by constructing the ancestral sequences derived from multiple sequence alignment; and (4) incorporating site-mutations targeting the evolutionarily coupled sites identified from multiple sequence alignment.