Quantification and functional analysis of modular protein evolution in a dense phylogenetic tree

Modularity is a hallmark of molecular evolution. Whether considering gene regulation, the components of metabolic pathways or signaling cascades, the ability to reuse autonomous modules in different molecular contexts can expedite evolutionary innovation. Similarly, protein domains are the modules of proteins, and modular domain rearrangements can create diversity with seemingly few operations in turn allowing for swift changes to an organism's functional repertoire. Here, we assess the patterns and functional effects of modular rearrangements at high resolution. Using a well resolved and diverse group of pancrustaceans, we illustrate arrangement diversity within closely related organisms, estimate arrangement turnover frequency and establish, for the first time, branch-specific rate estimates for fusion, fission, domain addition and terminal loss. Our results show that roughly 16 new arrangements arise per million years and that between 64% and 81% of these can be explained by simple, single-step modular rearrangement events. We find evidence that the frequencies of fission and terminal deletion events increase over time, and that modular rearrangements impact all levels of the cellular signaling apparatus and thus may have strong adaptive potential. Novel arrangements that cannot be explained by simple modular rearrangements contain a significant amount of repeat domains that occur in complex patterns which we term “supra-repeats”. Furthermore, these arrangements are significantly longer than those with a single-step rearrangement solution, suggesting that such arrangements may result from multi-step events. In summary, our analysis provides an integrated view and initial quantification of the patterns and functional impact of modular protein evolution in a well resolved phylogenetic tree. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

Emergence of bimodal cell population responses from the interplay between analog single-cell signaling and protein expression noise

ABSTRACT: BACKGROUND: Cell-to-cell variability in protein expression can be large, and its propagation through signaling networks affects biological outcomes. Here, we apply deterministic and probabilistic models and biochemical measurements to study how network topologies and cell-to-cell protein abundance variations interact to shape signaling responses. RESULTS: We observe bimodal distributions of extracellular signal-regulated kinase (ERK) responses to epidermal growth factor (EGF) stimulation, which are generally thought to indicate bistable or ultrasensitive signaling behavior in single cells. Surprisingly, we find that a simple MAPK/ERK-cascade model with negative feedback that displays graded, analog ERK responses at a single cell level can explain the experimentally observed bimodality at the cell population level. Model analysis suggests that a conversion of graded input--output responses in single cells to digital responses at the population level is caused by a broad distribution of ERK pathway activation thresholds brought about by cell-to-cell variability in protein expression. CONCLUSIONS: Our results show that bimodal signaling response distributions do not necessarily imply digital (ultrasensitive or bistable) single cell signaling, and the interplay between protein expression noise and network topologies can bring about digital population responses from analog single cell dose responses. Thus, cells can retain the benefits of robustness arising from negative feedback, while simultaneously generating population-level on/off responses that are thought to be critical for regulating cell fate decisions.     

Reassembly of JIP1 Scaffold Complex in JNK MAP Kinase Pathway Using Heterologous Protein Interactions

Formation of signaling protein complexes is crucial for proper signal transduction. Scaffold proteins in MAP kinase pathways are thought to facilitate complex assembly, thereby promoting efficient and specific signaling. To elucidate the assembly mechanism of scaffold complexes in mammals, we attempted to rationally rewire JIP1-dependent JNK MAP kinase pathway via alternative assembly of JIP1 complex. When JIP1-JNK docking interaction in the complex was replaced with heterologous protein interaction domains, such as PDZ domains and JNK-binding domains, a functional scaffold complex was reconstituted, and JNK signaling was rescued. Reassembly of JIP1 complex using heterologous protein interactions was sufficient for restoring of JNK MAP kinase pathway to induce signaling responses, including JNK activation and cell death. These results suggest a simple yet modular mechanism for JIP1 scaffold assembly in mammals.     

Substrate competition as a source of ultrasensitivity in the inactivation of Wee1

The mitotic regulators Wee1 and Cdk1 can inactivate each other through inhibitory phosphorylations. This double-negative feedback loop is part of a bistable trigger that makes the transition into mitosis abrupt and decisive. To generate a bistable response, some component of a double-negative feedback loop must exhibit an ultrasensitive response to its upstream regulator. Here, we experimentally demonstrate that Wee1 exhibits a highly ultrasensitive response to Cdk1. Several mechanisms can, in principle, give rise to ultrasensitivity, including zero-order effects, multisite phosphorylation, and competition mechanisms. We found that the ultrasensitivity in the inactivation of Wee1 arises mainly through two competition mechanisms: competition between two sets of phosphorylation sites in Wee1 and between Wee1 and other high-affinity Cdk1 targets. Based on these findings, we were able to reconstitute a highly ultrasensitive Wee1 response with purified components. Competition provides a simple way of generating the equivalent of a highly cooperative allosteric response.     

Molecular titration and ultrasensitivity in regulatory networks

Protein sequestration occurs when an active protein is sequestered by a repressor into an inactive complex. Using mathematical and computational modeling, we show how this regulatory mechanism (called “molecular titration”) can generate ultrasensitive or “all-or-none” responses that are equivalent to highly cooperative processes. The ultrasensitive nature of the input-output response is mainly determined by two parameters: the dimer dissociation constant and the repressor concentration. Because in vivo concentrations are tunable through a variety of mechanisms, molecular titration represents a flexible mechanism for generating ultrasensitivity. Using physiological parameters, we report how details of in vivo protein degradation affect the strength of the ultrasensitivity at steady state. Given that developmental systems often transduce signals into cell-fate decisions on timescales incompatible with steady state, we further examine whether molecular titration can produce ultrasensitive responses within physiologically relevant time intervals. Using Drosophila somatic sex determination as a developmental paradigm, we demonstrate that molecular titration can generate ultrasensitivity on timescales compatible with most cell-fate decisions. Gene duplication followed by loss-of-function mutations can create dominant negatives that titrate and compete with the original protein. Dominant negatives are abundant in gene regulatory circuits, and our results suggest that molecular titration might be generating an ultrasensitive response in these networks.     

A three-step proteolytic cascade mediates the activation of the peptidoglycan-induced toll pathway in an insect

The recognition of lysine-type peptidoglycans (PG) by the PG recognition complex has been suggested to cause activation of the serine protease cascade leading to the processing of Sp?tzle and subsequent activation of the Toll signaling pathway. So far, two serine proteases involved in the lysine-type PG Toll signaling pathway have been identified. One is a modular serine protease functioning as an initial enzyme to be recruited into the lysine-type PG recognition complex. The other is the Drosophila Sp?tzle processing enzyme (SPE), a terminal enzyme that converts Sp?tzle pro-protein to its processed form capable of binding to the Toll receptor. However, it remains unclear how the initial PG recognition signal is transferred to Sp?tzle resulting in Toll pathway activation. Also, the biochemical characteristics and mechanism of action of a serine protease linking the modular serine protease and SPE have not been investigated. Here, we purified and cloned a novel upstream serine protease of SPE that we named SAE, SPE-activating enzyme, from the hemolymph of a large beetle, Tenebrio molitor larvae. This enzyme was activated by Tenebrio modular serine protease and in turn activated the Tenebrio SPE. The biochemical ordered functions of these three serine proteases were determined in vitro, suggesting that the activation of a three-step proteolytic cascade is necessary and sufficient for lysine-type PG recognition signaling. The processed Sp?tzle by this cascade induced antibacterial activity in vivo. These results demonstrate that the three-step proteolytic cascade linking the PG recognition complex and Sp?tzle processing is essential for the PG-dependent Toll signaling pathway.     

Study of Circular Dichroism Modes Through Decomposition of Planar Nanostructures

We extend the plasmon hybridization method from a single nanoparticle to a complex planar nanostructure, decomposing the complex nanostructure into fundamental nanoparticle building blocks. Using gammadion nanostructure as an example, we validated the theory by comparing the field profile in the gammadion’s arms under the influence of an incident circularly polarized wave. This allows us to address the origin of the plasmonics modes in the circular dichroism (CD) spectrum. The use of this hybridization method provides a simple and intuitive explanation on how conductive and inductive coupling may result from complex planar nanostructures, allowing us to study its optical properties. Using our approach, top down hybridization studies can be applied to other complex planar structures to gain further insight on the origin of the CD modes and enhance ultrasensitive sensing of chiral micro and macro molecules.     

Equilibrium analysis of allosteric interactions shows zero-order effects

The binding of effector to an allosteric protein exhibits a non-Michaelis-Menten behavior, resulting in either ultrasensitive or subsensitive response. In the present work, a modular approach has been developed to determine the response curve for allosteric systems at higher concentration of allosteric enzyme than that of effector (zero-order sensitivity, as observed in enzyme cascades) by equilibrium analysis. The analysis shows that, in an allosteric system, the zero-order effect can make the response ultrasensitive or subsensitive with respect to the enzyme concentration. The response is dependent on the number of binding sites, cooperativity, and the total effector concentration. The framework was further applied to a well studied allosteric protein, the Escherichia coli aspartate transcarbamoylase. The predictions are found to be consistent with the reported experimental data.     

A Bmp Reporter with Ultrasensitive Characteristics Reveals That High Bmp Signaling Is Not Required for Cortical Hem Fate

Insights into Bone morphogenetic protein (Bmp) functions during forebrain development have been limited by a lack of Bmp signaling readouts. Here we used a novel Bmp signaling reporter (“BRE-gal” mice) to study Bmp signaling in the dorsal telencephalon. At early stages, BRE-gal expression was restricted to the dorsal telencephalic midline. At later stages, strong BRE-gal expression occurred in neurons of the marginal zone and dentate gyrus. Comparisons to nuclear phospho-Smad1/5/8 (pSmad) and Msx1 indicated that BRE-gal expression occurred exclusively in neural cells with high-level Bmp signaling. BRE-gal responsiveness to Bmps was confirmed in reporter-negative cortical cells cultured with Bmp4, and both in vivo and in vitro, BRE-gal expression was switch-like, or ultrasensitive. In the early dorsal telencephalon, BRE-gal expression negatively correlated with the cortical selector gene Lhx2, indicating a BRE-gal expression border that coincides with the cortex-hem boundary. However, in Lhx2 null chimeras, neither BRE-gal nor nuclear pSmad increases were observed in ectopic hem cells. These findings establish BRE-gal as an ultrasensitive reporter of Bmp signaling in the dorsal telencephalon, imply that hem fate can be specified at different Bmp signaling intensities, and suggest that Lhx2 primarily regulates the responses to – rather than the intensity of – Bmp signaling in dorsal telencephalic cells.     

Modular engineering of cellular signaling proteins and networks

Living cells respond to their environment using networks of signaling molecules that act as sensors, information processors, and actuators. These signaling systems are highly modular at both the molecular and network scales, and much evidence suggests that evolution has harnessed this modularity to rewire and generate new physiological behaviors. Conversely, we are now finding that, following nature's example, signaling modules can be recombined to form synthetic tools for monitoring, interrogating, and controlling the behavior of cells. Here we highlight recent progress in the modular design of synthetic receptors, optogenetic switches, and phospho-regulated proteins and circuits, and discuss the expanding role of combinatorial design in the engineering of cellular signaling proteins and networks.     

From intracellular signaling to population oscillations: bridging size- and time-scales in collective behavior

Collective behavior in cellular populations is coordinated by biochemical signaling networks within individual cells. Connecting the dynamics of these intracellular networks to the population phenomena they control poses a considerable challenge because of network complexity and our limited knowledge of kinetic parameters. However, from physical systems, we know that behavioral changes in the individual constituents of a collectively behaving system occur in a limited number of well-defined classes, and these can be described using simple models. Here, we apply such an approach to the emergence of collective oscillations in cellular populations of the social amoeba Dictyostelium discoideum. Through direct tests of our model with quantitative in vivo measurements of single-cell and population signaling dynamics, we show how a simple model can effectively describe a complex molecular signaling network at multiple size and temporal scales. The model predicts novel noise-driven single-cell and population-level signaling phenomena that we then experimentally observe. Our results suggest that like physical systems, collective behavior in biology may be universal and described using simple mathematical models.     

A bistable model of cell polarity

Ultrasensitivity, as described by Goldbeter and Koshland, has been considered for a long time as a way to realize bistable switches in biological systems. It is not as well recognized that when ultrasensitivity and reinforcing feedback loops are present in a spatially distributed system such as the cell plasmamembrane, they may induce bistability and spatial separation of the system into distinct signaling phases. Here we suggest that bistability of ultrasensitive signaling pathways in a diffusive environment provides a basic mechanism to realize cell membrane polarity. Cell membrane polarization is a fundamental process implicated in several basic biological phenomena, such as differentiation, proliferation, migration and morphogenesis of unicellular and multicellular organisms. We describe a simple, solvable model of cell membrane polarization based on the coupling of membrane diffusion with bistable enzymatic dynamics. The model can reproduce a broad range of symmetry-breaking events, such as those observed in eukaryotic directional sensing, the apico-basal polarization of epithelium cells, the polarization of budding and mating yeast, and the formation of Ras nanoclusters in several cell types.     

Structure-based kinetic models of modular signaling protein function: focus on Shp2

We present here a computational, rule-based model to study the function of the SH2 domain-containing protein tyrosine phosphatase, Shp2, in intracellular signal transduction. The two SH2 domains of Shp2 differentially regulate the enzymatic activity by a well-characterized mechanism, but they also affect the targeting of Shp2 to signaling receptors in cells. Our kinetic model integrates these potentially competing effects by considering the intra- and intermolecular interactions of the Shp2 SH2 domains and catalytic site as well as the effect of Shp2 phosphorylation. Even for the isolated Shp2/receptor system, which may seem simple by certain standards, we find that the network of possible binding and phosphorylation states is composed of over 1000 members. To our knowledge, this is the first kinetic model to fully consider the modular, multifunctional structure of a signaling protein, and the computational approach should be generally applicable to other complex intermolecular interactions.     

Ultrasensitive synthetic protein regulatory networks using mixed decoys

Cellular protein interaction networks exhibit sigmoidal input-output relationships with thresholds and steep responses (i.e. ultrasensitivity). Although cooperativity can be a source of ultrasensitivity, we examined whether the presence of "decoy" binding sites that are not coupled to activation could also lead to this effect. To systematically vary key parameters of the system, we designed a synthetic regulatory system consisting of an autoinhibited PDZ domain coupled to an activating SH3 domain binding site. In the absence of a decoy binding site, this system is non-ultrasensitive, as predicted by modeling of this system. Addition of a high-affinity decoy site adds a threshold, but the response is not ultrasensitive. We found that sigmoidal activation profiles can be generated utilizing multiple decoys with mixtures of high and low affinities, where high affinity decoys act to set the threshold and low affinity decoys ensure a sigmoidal response. Placing the synthetic decoy system in a mitotic spindle orientation cell culture system thresholds this physiological activity. Thus, simple combinations of non-activating binding sites can lead to complex regulatory responses in protein interaction networks.     

Specificity in signal transduction: from phosphotyrosine-SH2 domain interactions to complex cellular systems

Over the last two decades, a new and unifying concept of cellular organization has emerged in which modular protein-protein interactions provide an underlying framework through which signaling pathways are assembled and controlled. In this scheme, posttranslational modifications such as phosphorylation commonly exert their biological effects by regulating molecular interactions, exemplified by the ability of phosphotyrosine sites to bind selectively to SH2 domains. Although these interactions are rather simple in isolation, they can nonetheless be exploited to generate complex cellular systems. Here, I discuss experiments that have led to this view of dynamic cellular behavior and identify some current and future areas of interest in cell signaling.     

Hormone signaling in evolution and development: a non-model system approach

Cooption and modularity are informative concepts in evolutionary developmental biology. Genes function within complex networks that act as modules in development. These modules can then be coopted in various functional and evolutionary contexts. Hormonal signaling, the main focus of this review, has a modular character. By regulating the activities of genes, proteins and other cellular molecules, a hormonal signal can have major effects on physiological and ontogenetic processes within and across tissues over a wide spatial and temporal scale. Because of this property, we argue that hormones are frequently involved in the coordination of life history transitions (LHTs) and their evolution (LHE). Finally, we promote the usefulness of a comparative, non-model system approach towards understanding how hormones function and guide development and evolution, highlighting thyroid hormone function in echinoids as an example.     

Multisite phosphorylation and network dynamics of cyclin-dependent kinase signaling in the eukaryotic cell cycle

Multisite phosphorylation of regulatory proteins has been proposed to underlie ultrasensitive responses required to generate nontrivial dynamics in complex biological signaling networks. We used a random search strategy to analyze the role of multisite phosphorylation of key proteins regulating cyclin-dependent kinase (CDK) activity in a model of the eukaryotic cell cycle. We show that multisite phosphorylation of either CDK, CDC25, wee1, or CDK-activating kinase is sufficient to generate dynamical behaviors including bistability and limit cycles. Moreover, combining multiple feedback loops based on multisite phosphorylation do not destabilize the cell cycle network by inducing complex behavior, but rather increase the overall robustness of the network. In this model we find that bistability is the major dynamical behavior of the CDK signaling network, and that negative feedback converts bistability into limit cycle behavior. We also compare the dynamical behavior of several simplified models of CDK regulation to the fully detailed model. In summary, our findings suggest that multisite phosphorylation of proteins is a critical biological mechanism in generating the essential dynamics and ensuring robust behavior of the cell cycle.