The intrinsic, or mitochondrial, pathway of caspase activation is essential for apoptosis induction by various stimuli including cytotoxic stress. It depends on the cellular context, whether cytochrome c released from mitochondria induces caspase activation gradually or in an all-or-none fashion, and whether caspase activation irreversibly commits cells to apoptosis. By analyzing a quantitative kinetic model, we show that inhibition of caspase-3 (Casp3) and Casp9 by inhibitors of apoptosis (IAPs) results in an implicit positive feedback, since cleaved Casp3 augments its own activation by sequestering IAPs away from Casp9. We demonstrate that this positive feedback brings about bistability (i.e., all-or-none behaviour), and that it cooperates with Casp3-mediated feedback cleavage of Casp9 to generate irreversibility in caspase activation. Our calculations also unravel how cell-specific protein expression brings about the observed qualitative differences in caspase activation (gradual versus all-or-none and reversible versus irreversible). Finally, known regulators of the pathway are shown to efficiently shift the apoptotic threshold stimulus, suggesting that the bistable caspase cascade computes multiple inputs into an all-or-none caspase output. As cellular inhibitory proteins (e.g., IAPs) frequently inhibit consecutive intermediates in cellular signaling cascades (e.g., Casp3 and Casp9), the feedback mechanism described in this paper is likely to be a widespread principle on how cells achieve ultrasensitivity, bistability, and irreversibility.

Crosstalk mechanisms have not been studied as thoroughly as individual signaling pathways. We exploit experimental and computational approaches to reveal how a concordant interplay between the insulin and epidermal growth factor (EGF) signaling networks can potentiate mitogenic signaling. In HEK293 cells, insulin is a poor activator of the Ras/ERK (extracellular signal-regulated kinase) cascade, yet it enhances ERK activation by low EGF doses. We find that major crosstalk mechanisms that amplify ERK signaling are localized upstream of Ras and at the Ras/Raf level. Computational modeling unveils how critical network nodes, the adaptor proteins GAB1 and insulin receptor substrate (IRS), Src kinase, and phosphatase SHP2, convert insulin-induced increase in the phosphatidylinositol-3,4,5-triphosphate (PIP 3 ) concentration into enhanced Ras/ERK activity. The model predicts and experiments confirm that insulin-induced amplification of mitogenic signaling is abolished by disrupting PIP 3 -mediated positive feedback via GAB1 and IRS. We demonstrate that GAB1 behaves as a non-linear amplifier of mitogenic responses and insulin endows EGF signaling with robustness to GAB1 suppression. Our results show the feasibility of using computational models to identify key target combinations and predict complex cellular responses to a mixture of external cues.

In most biological organisms intracellular signal processing is carried out by networks composed of enzymes that activate and inactivate each other by covalent modification. Signals received at the cell membrane ripple through signalling networks via covalent modification events to reach various locations in the cell and ultimately cause cellular responses. The biochemical building blocks of these networks are frequently enzyme pairs, such as a kinase and a phosphatase, that form covalent modification cycles in which the target enzyme is covalently modified at single or multiple sites in a reversible manner.In some experiments, the stimulus-response curves display strong sigmoidal dependencies in vivo, for example, in the activation of the mitogen-activated protein kinase (MAPK) cascade [2] and Sic1 [3], and in vitro, for example, in the phosphorylation of isocitrate dehydrogenase [4], muscle glycolysis [5] and in postsynaptic calcium signalling [6]. Sigmoidal stimulus-response curves imply that the responses are highly sensitive to changes in signals around the threshold level. Thus it is more sensitive than a typical Michaelis-Menten-like response, a property that has been termed ultrasensitivity [1].

The cytokine TGFβ provides important information during embryonic development, adult tissue homeostasis, and regeneration. Alterations in the cellular response to TGFβ are involved in severe human diseases. To understand how cells encode the extracellular input and transmit its information to elicit appropriate responses, we acquired quantitative time‐resolved measurements of pathway activation at the single‐cell level. We established dynamic time warping to quantitatively compare signaling dynamics of thousands of individual cells and described heterogeneous single‐cell responses by mathematical modeling. Our combined experimental and theoretical study revealed that the response to a given dose of TGFβ is determined cell specifically by the levels of defined signaling proteins. This heterogeneity in signaling protein expression leads to decomposition of cells into classes with qualitatively distinct signaling dynamics and phenotypic outcome. Negative feedback regulators promote heterogeneous signaling, as a SMAD7 knock‐out specifically affected the signal duration in a subpopulation of cells. Taken together, we propose a quantitative framework that allows predicting and testing sources of cellular signaling heterogeneity.

Cellular decision‐making and environmental adaptation are dependent upon a heterogeneous response of gene expression to external cues. Heterogeneity arises in transcription from random switching between transcriptionally active and inactive states, resulting in bursts of RNA synthesis. Furthermore, the cellular state influences the competency of transcription, thereby globally affecting gene expression in a cell‐specific manner. We determined how external stimuli interplay with cellular state to modulate the kinetics of bursting. To this end, single‐cell dynamics of nascent transcripts were monitored at the endogenous estrogen‐responsive GREB1 locus. Stochastic modeling of gene expression implicated a two‐state promoter model in which the estrogen stimulus modulates the frequency of transcriptional bursting. The cellular state affects transcriptional dynamics by altering initiation and elongation kinetics and acts globally, as GREB1 alleles in the same cell correlate in their transcriptional output. Our results suggest that cellular state strongly affects the first step of the central dogma of gene expression, to promote heterogeneity in the transcriptional output of isogenic cells.