Computational study of the mechanism of Bcl-2 apoptotic switch

Programmed cell death - apoptosis is one of the most studied biological phenomenon of recent years. Apoptotic regulatory network contains several significant control points, including probably the most important one - Bcl--2 apoptotic switch. There are two proposed hypotheses regarding its internal working - the indirect activation and direct activation models. Since these hypotheses form extreme poles of full continuum of intermediate models, we have constructed more general model with these two models as extreme cases. By studying relationship between model parameters and steady-state response ultrasensitivity we have found optimal interaction pattern which reproduces behavior of Bcl-2 apoptotic switch. Our results show, that stimulus-response ultrasensitivity is negatively related to spontaneous activation of Bcl-2 effectors - subgroup of Bcl-2 proteins. We found that ultrasensitivity requires effector's activation, mediated by another subgroup of Bcl-2 proteins - activators. We have shown that the auto-activation of effectors forms ultrasensitivity enhancing feedback loop, only if mediated by monomers, but not by oligomers. Robustness analysis revealed that interaction pattern proposed by direct activation hypothesis is able to conserve stimulus-response dependence and preserve ultrasensitivity despite large changes of its internal parameters. This ability is strongly reduced as for the intermediate to indirect side of the models. Computer simulation of the more general model presented here suggest, that stimulus-response ultrasensitivity is an emergent property of the direct activation model, that cannot originate within model of indirect activation. Introduction of indirect-model-specific interactions does not provide better explanation of Bcl-2 functioning compared to direct model

Enzyme sequestration as a tuning point in controlling response dynamics of signalling networks

Signalling networks result from combinatorial interactions among many enzymes and scaffolding proteins. These complex systems generate response dynamics that are often essential for correct decision-making in cells. Uncovering biochemical design principles that underpin such response dynamics is a prerequisite to understand evolved signalling networks and to design synthetic ones. Here, we use in silico evolution to explore the possible biochemical design space for signalling networks displaying ultrasensitive and adaptive response dynamics. By running evolutionary simulations mimicking different biochemical scenarios, we find that enzyme sequestration emerges as a key mechanism for enabling such dynamics. Inspired by these findings, and to test the role of sequestration, we design a generic, minimalist model of a signalling cycle, featuring two enzymes and a single scaffolding protein. We show that this simple system is capable of displaying both ultrasensitive and adaptive response dynamics. Furthermore, we find that tuning the concentration or kinetics of the sequestering protein can shift system dynamics between these two response types. These empirical results suggest that enzyme sequestration through scaffolding proteins is exploited by evolution to generate diverse response dynamics in signalling networks and could provide an engineering point in synthetic biology applications

Impact of upstream and downstream constraints on a signaling module's ultrasensitivity

Much work has been done on the study of the biochemical mechanisms that result in ultrasensitive behavior of simple biochemical modules. However, in a living cell, such modules are embedded in a bigger network that constrains the range of inputs that the module will receive as well as the range of the module’s outputs that network will be able to detect. Here, we studied how the effective ultrasensitivity of a modular system is affected by these restrictions. We use a simple setup to explore to what extent the dynamic range spanned by upstream and downstream components of an ultrasensitive module impact on the effective sensitivity of the system. Interestingly, we found for some ultrasensitive motifs that dynamic range limitations imposed by downstream components can produce effective sensitivities much larger than that of the original module when considered in isolation.Fil: Altszyler Lemcovich, Edgar Jaim. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Fisiología, Biología Molecular y Neurociencias. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Fisiología, Biología Molecular y Neurociencias; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Ventura, Alejandra. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Fisiología, Biología Molecular y Neurociencias. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Fisiología, Biología Molecular y Neurociencias; ArgentinaFil: Colman Lerner, Alejandro Ariel. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Fisiología, Biología Molecular y Neurociencias. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Fisiología, Biología Molecular y Neurociencias; ArgentinaFil: Chernomoretz, Ariel. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina. Fundación Instituto Leloir; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; Argentin

Response dynamics of phosphorelays suggest their potential utility in cell signalling

Phosphorelays are extended two-component signalling systems found in diverse bacteria, lower eukaryotes and plants. Only few of these systems are characterized, and we still lack a full understanding of their signalling abilities. Here, we aim to achieve a global understanding of phosphorelay signalling and its dynamical properties. We develop a generic model, allowing us to systematically analyse response dynamics under different assumptions. Using this model, we find that the steady-state concentration of phosphorylated protein at the final layer of a phosphorelay is a linearly increasing, but eventually saturating function of the input. In contrast, the intermediate layers can display ultrasensitivity. We find that such ultrasensitivity is a direct result of the phosphorelay biochemistry; shuttling of a single phosphate group from the first to the last layer. The response dynamics of the phosphorelay results in tolerance of cross-talk, especially when it occurs as cross-deactivation. Further, it leads to a high signal-to-noise ratio for the final layer. We find that a relay length of four, which is most commonly observed, acts as a saturating point for these dynamic properties. These findings suggest that phosphorelays could act as a mechanism to reduce noise and effects of cross-talk on the final layer of the relay and enforce its input–response relation to be linear. In addition, our analysis suggests that middle layers of phosphorelays could embed thresholds. We discuss the consequence of these findings in relation to why cells might use phosphorelays along with enzymatic kinase cascades

BioJazz : In silico evolution of cellular networks with unbounded complexity using rule-based modeling

Systems biologists aim to decipher the structure and dynamics of signaling and regulatory networks underpinning cellular responses; synthetic biologists can use this insight to alter existing networks or engineer de novo ones. Both tasks will benefit from an understanding of which structural and dynamic features of networks can emerge from evolutionary processes, through which intermediary steps these arise, and whether they embody general design principles. As natural evolution at the level of network dynamics is difficult to study, in silico evolution of network models can provide important insights. However, current tools used for in silico evolution of network dynamics are limited to ad hoc computer simulations and models. Here we introduce BioJazz, an extendable, user-friendly tool for simulating the evolution of dynamic biochemical networks. Unlike previous tools for in silico evolution, BioJazz allows for the evolution of cellular networks with unbounded complexity by combining rule-based modeling with an encoding of networks that is akin to a genome. We show that BioJazz can be used to implement biologically realistic selective pressures and allows exploration of the space of network architectures and dynamics that implement prescribed physiological functions. BioJazz is provided as an open-source tool to facilitate its further development and use. Source code and user manuals are available at: http://oss-lab.github.io/biojazz and http://osslab.lifesci.warwick.ac.uk/BioJazz.aspx

Biosensor Architectures for High-Fidelity Reporting of Cellular Signaling

Understanding mechanisms of information processing in cellular signaling networks requires quantitative measurements of protein activities in living cells. Biosensors are molecular probes that have been developed to directly track the activity of specific signaling proteins and their use is revolutionizing our understanding of signal transduction. The use of biosensors relies on the assumption that their activity is linearly proportional to the activity of the signaling protein they have been engineered to track. We use mechanistic mathematical models of common biosensor architectures (single-chain FRET-based biosensors), which include both intramolecular and intermolecular reactions, to study the validity of the linearity assumption. As a result of the classic mechanism of zero-order ultrasensitivity, we find that biosensor activity can be highly nonlinear so that small changes in signaling protein activity can give rise to large changes in biosensor activity and vice versa. This nonlinearity is abolished in architectures that favor the formation of biosensor oligomers, but oligomeric biosensors produce complicated FRET states. Based on this finding, we show that high-fidelity reporting is possible when a single-chain intermolecular biosensor is used that cannot undergo intramolecular reactions and is restricted to forming dimers. We provide phase diagrams that compare various trade-offs, including observer effects, which further highlight the utility of biosensor architectures that favor intermolecular over intramolecular binding. We discuss challenges in calibrating and constructing biosensors and highlight the utility of mathematical models in designing novel probes for cellular signaling

Spatio-temporal correlations can drastically change the response of a MAPK pathway

Multisite covalent modification of proteins is omnipresent in eukaryotic cells. A well-known example is the mitogen-activated protein kinase (MAPK) cascade, where in each layer of the cascade a protein is phosphorylated at two sites. It has long been known that the response of a MAPK pathway strongly depends on whether the enzymes that modify the protein act processively or distributively: distributive mechanism, in which the enzyme molecules have to release the substrate molecules in between the modification of the two sites, can generate an ultrasensitive response and lead to hysteresis and bistability. We study by Green's Function Reaction Dynamics, a stochastic scheme that makes it possible to simulate biochemical networks at the particle level and in time and space, a dual phosphorylation cycle in which the enzymes act according to a distributive mechanism. We find that the response of this network can differ dramatically from that predicted by a mean-field analysis based on the chemical rate equations. In particular, rapid rebindings of the enzyme molecules to the substrate molecules after modification of the first site can markedly speed up the response, and lead to loss of ultrasensitivity and bistability. In essence, rapid enzyme-substrate rebindings can turn a distributive mechanism into a processive mechanism. We argue that slow ADP release by the enzymes can protect the system against these rapid rebindings, thus enabling ultrasensitivity and bistability

Threshold responses to morphogen gradients by zero-order ultrasensitivity

Translating a graded morphogen distribution into tight response borders is central to all developmental processes. Yet, the molecular mechanisms generating such behavior are poorly understood. During patterning of the Drosophila embryonic ventral ectoderm, a graded mitogen-activated protein kinase (MAPK) activation is converted into an all-or-none degradation switch of the Yan transcriptional repressor. Replacing the cardinal phosphorylated amino acid of Yan by a phosphomimetic residue allowed its degradation in a MAPK-independent manner, consistent with Yan phosphorylation being the critical event in generating the switch. Several alternative threshold mechanisms that could, in principle, be realized by this phosphorylation, including first order, cooperativity, positive feedback and zero-order ultrasensitivity, were analyzed. We found that they can be distinguished by their kinetics and steady-state responses to Yan overexpression. In agreement with the predictions for zero-order kinetics, an increase in Yan levels did not shift the degradation border, but significantly elevated the time required to reach steady state. We propose that a reversible loop of Yan phosphorylation implements a zero-order ultrasensitivity-like threshold mechanism, with the capacity to form sharp thresholds that are independent of the level of Yan

Forward and Reverse Engineering of Cellular Decision-Making

Cells reside in highly dynamic environments to which they must adapt. Throughout its lifetime, an individual cell receives numerous chemical and mechanical signals, communicated through dense molecular networks, and eliciting a diverse array of responses. A large number of these signals necessitate discrete, all-or-none responses. For instance, a cell receiving proliferation signals must respond by committing to the cell-cycle and dividing, or by not initiating the process at all; that is, the cell must not adopt an intermediate route. Analogously, a stem-cell receiving signals for different lineages must commit exclusively to one of these lineages. How individual cells integrate multiple, possibly conflicting, noisy inputs, and make discrete decisions is poorly understood. Detailed insight into cellular decision-making can enable cell-based therapies, shed light on diseases arising out of dysregulation of control, and suggest practical design strategies for implementing this behavior in synthetic systems for research and industrial use. In this thesis, we have employed both mathematical modeling and experiments to further elucidate the mechanistic underpinnings of decision-making in cells. First, we describe a computational study that assesses the entire space of minimal networks to identify topologies that can not only make decisions but can do so robustly in the dynamic and noisy cellular environment. Second, via model-driven, quantitative experiments in a megakaryocyte erythroid progenitor line, we demonstrate that a simple network with mutual antagonism and autoregulation captures the dynamics of the master transcription factors at the level of individual cells. Expansion of this model to account for extrinsic cues reconciles the competing stochastic and instructive theories of hematopoietic lineage commitment, and implicates cytokine receptors in broader regulatory roles. Third, to assess the impact of specific genetic perturbations on the distribution of the population, and on commitment trajectories of individual cells, we implemented the core mutual antagonism and autoregulation topology synthetically in yeast cells. Our approach of using orthogonal variants of a single core protein represents a general, modular design strategy for building synthetic circuits, and model-driven experiments elucidate how gene dosage, repression strength, and promoter architecture can modulate decision-making behavior