Mitogen-activated protein kinase (MAPK) cascades can operate as bistable switches residing in either of two different stable states. MAPK cascades are often embedded in positive feedback loops, which are considered to be a prerequisite for bistable behavior. Here we demonstrate that in the absence of any imposed feedback regulation, bistability and hysteresis can arise solely from a distributive kinetic mechanism of the two-site MAPK phosphorylation and dephosphorylation. Importantly, the reported kinetic properties of the kinase (MEK) and phosphatase (MKP3) of extracellular signal–regulated kinase (ERK) fulfill the essential requirements for generating a bistable switch at a single MAPK cascade level. Likewise, a cycle where multisite phosphorylations are performed by different kinases, but dephosphorylation reactions are catalyzed by the same phosphatase, can also exhibit bistability and hysteresis. Hence, bistability induced by multisite covalent modification may be a widespread mechanism of the control of protein activity.

Grb2-associated binder 1 (GAB1) is a scaffold protein involved in numerous interactions that propagate signaling by growth factor and cytokine receptors. Here we explore in silico and validate in vivo the role of GAB1 in the control of mitogenic (Ras/MAPK) and survival (phosphatidylinositol 3-kinase (PI3K)/Akt) signaling stimulated by epidermal growth factor (EGF). We built a comprehensive mechanistic model that allows for reliable predictions of temporal patterns of cellular responses to EGF under diverse perturbations, including different EGF doses, GAB1 suppression, expression of mutant proteins, and pharmacological inhibitors. We show that the temporal dynamics of GAB1 tyrosine phosphorylation is significantly controlled by positive GAB1-PI3K feedback and negative MAPK-GAB1 feedback. Our experimental and computational results demonstrate that the essential function of GAB1 is to enhance PI3K/Akt activation and extend the duration of Ras/MAPK signaling. By amplifying positive interactions between survival and mitogenic pathways, GAB1 plays the critical role in cell proliferation and tumorigenesis.

After activation, many receptors and their adaptor proteins act as scaffolds displaying numerous docking sites and engaging multiple targets. The consequent assemblage of a variety of protein complexes results in a combinatorial increase in the number of feasible molecular species presenting different states of a receptor-scaffold signaling module. Tens of thousands of such microstates emerge even for the initial signal propagation events, greatly impeding a quantitative analysis of networks. Here, we demonstrate that the assumption of independence of molecular events occurring at distinct sites enables us to approximate a mechanistic picture of all possible microstates by a macrodescription of states of separate domains, i.e., macrostates that correspond to experimentally verifiable variables. This analysis dissects a highly branched network into interacting pathways originated by protein complexes assembled on different sites of receptors and scaffolds. We specify when the temporal dynamics of any given microstate can be expressed using the product of the relative concentrations of individual sites. The methods presented here are equally applicable to deterministic and stochastic calculations of the temporal dynamics. Our domain-oriented approach drastically reduces the number of states, processes, and kinetic parameters to be considered for quantification of complex signaling networks that propagate distinct physiological responses.

A hallmark of protein kinase/phosphatase cascades, including mitogen-activated protein kinase (MAPK) pathways, is the spatial separation of their components within cells. The top-level kinase, MAP3K, is phosphorylated at the cell membrane, and cytoplasmic kinases at sequential downstream levels (MAP2K and MAPK) spread the signal to distant targets. Given measured protein diffusivity and phosphatase activities, signal propagation by diffusion would result in a steep decline of MAP2K activity and low bisphosphorylated MAPK (ppMAPK) levels near the nucleus, especially in large cells, such as oocytes. Here, we show that bistability in a two-site MAPK (de)phosphorylation cycle generates a novel type of phosphoprotein wave that propagates from the surface deep into the cell interior. Positive feedback from ppMAPK to cytoplasmic MAP2K enhances the propagation span of the ppMAPK wave, making it possible to convey phosphorylation signals over exceedingly long distances. The finding of phosphorylation waves traveling with constant amplitude and high velocity may solve a long-standing enigma of survival signaling in developing neurons.

A systems biology approach is applied to gain a quantitative understanding of the integration of signalling by the small GTPase Ras. The Ras protein acts as a critical switch in response to signals that determine the cell's fate. In unstimulated cells, Ras switching between an inactive GDP-binding and active GTP-binding state is controlled by the intrinsic catalytic activities of Ras. The calculated high sensitivity of the basal Ras-GTP fraction to changes in the rate constant of GTP-hydrolysis by Ras can account for the carcinogenic potential of Ras mutants with decreased GTPase activities. Extracelluar stimuli initiate Ras interactions with GDP/GTP exchange factors such as SOS, and GTP-hydrolysis activating proteins such as RasGAP. Our data on freshly isolated hepatocytes stimulated with epidermal growth factor (EGF) show transient SOS activation and sustained Ras-GTP patterns. We demonstrate that these dose-response data can only be explained by transient RasGAP activitation, and not by merely switching off the SOS signal, e.g. by inhibitory phosphorylation of SOS. A transient RasGAP activity can be brought about by a number of mechanisms. A comprehensive kinetic model of the EGF receptor (EGFR) network was developed to explore feasible molecular scenarios, including the receptor-mediated recruitment of SOS and RasGAP to the plasma membrane, phosphorylation of RasGAP and p190 RhoGAP by soluble tyrosine kinases, and RasGAP interactions with phosphoinositides and p190 RhoGAP. We show that a transient RasGAP association with EGFR followed by the capture of RasGAP through the formation of complexes with p190 RhoGAP can account for data on hepatocytes. In summary, our results demonstrate that a combination of experimental monitoring and integrated dynamic analysis is capable of dissecting regulatory mechanisms that govern cellular signal transduction.