Boris N. Kholodenko

Cell-signalling dynamics in time and space

The specificity of cellular responses to receptor stimulation is encoded by the spatial and temporal dynamics of downstream signalling networks. Temporal dynamics are coupled to spatial gradients of signalling activities, which guide pivotal intracellular processes and tightly regulate signal propagation across a cell. Computational models provide insights into the complex relationships between the stimuli and the cellular responses, and reveal the mechanisms that are responsible for signal amplification, noise reduction and generation of discontinuous bistable dynamics or oscillations.

Signaling switches and bistability arising from multisite phosphorylation in protein kinase cascades

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.

Quantification of Short Term Signaling by the Epidermal Growth Factor Receptor

During the past decade, our knowledge of molecular mechanisms involved in growth factor signaling has proliferated almost explosively. However, the kinetics and control of information transfer through signaling networks remain poorly understood. This paper combines experimental kinetic analysis and computational modeling of the short term pattern of cellular responses to epidermal growth factor (EGF) in isolated hepatocytes. The experimental data show transient tyrosine phosphorylation of the EGF receptor (EGFR) and transient or sustained response patterns in multiple signaling proteins targeted by EGFR. Transient responses exhibit pronounced maxima, reached within 15–30 s of EGF stimulation and followed by a decline to relatively low (quasi-steady-state) levels. In contrast to earlier suggestions, we demonstrate that the experimentally observed transients can be accounted for without requiring receptor-mediated activation of specific tyrosine phosphatases, following EGF stimulation. The kinetic model predicts how the cellular response is controlled by the relative levels and activity states of signaling proteins and under what conditions activation patterns are transient or sustained. EGFR signaling patterns appear to be robust with respect to variations in many elemental rate constants within the range of experimentally measured values. On the other hand, we specify which changes in the kinetic scheme, rate constants, and total amounts of molecular factors involved are incompatible with the experimentally observed kinetics of signal transfer. Quantitation of signaling network responses to growth factors allows us to assess how cells process information controlling their growth and differentiation.

Signalling ballet in space and time

Although we have amassed extensive catalogues of signalling network components, our understanding of the spatiotemporal control of emergent network structures has lagged behind. Dynamic behaviour is starting to be explored throughout the genome, but analysis of spatial behaviours is still confined to individual proteins. The challenge is to reveal how cells integrate temporal and spatial information to determine specific biological functions. Key findings are the discovery of molecular signalling machines such as Ras nanoclusters, spatial activity gradients and flexible network circuitries that involve transcriptional feedback. They reveal design principles of spatiotemporal organization that are crucial for network function and cell fate decisions.

Untangling the wires: A strategy to trace functional interactions in signaling and gene networks

Emerging technologies have enabled the acquisition of large genomics and proteomics data sets. However, current methodologies for analysis do not permit interpretation of the data in ways that unravel cellular networking. We propose a quantitative method for determining functional interactions in cellular signaling and gene networks. It can be used to explore cell systems at a mechanistic level or applied within a “modular” framework, which dramatically decreases the number of variables to be assayed. This method is based on a mathematical derivation that demonstrates how the topology and strength of network connections can be retrieved from experimentally measured network responses to successive perturbations of all modules. Importantly, our analysis can reveal functional interactions even when the components of the system are not all known. Under these circumstances, some connections retrieved by the analysis will not be direct but correspond to the interaction routes through unidentified elements. The method is tested and illustrated by using computer-generated responses of a modeled mitogen-activated protein kinase cascade and gene network.

Cross-talk between mitogenic Ras/MAPK and survival PI3K/Akt pathways: a fine balance

In the present paper, we describe multiple levels of cross-talk between the PI3K (phosphoinositide 3-kinase)/Akt and Ras/MAPK (mitogen-activated protein kinase) signalling pathways. Experimental data and computer simulations demonstrate that cross-talk is context-dependent and that both pathways can activate or inhibit each other. Positive influence of the PI3K pathway on the MAPK pathway is most effective at sufficiently low doses of growth factors, whereas negative influence of the MAPK pathway on the PI3K pathway is mostly pronounced at high doses of growth factors. Pathway cross-talk endows a cell with emerging capabilities for processing and decoding signals from multiple receptors activated by different combinations of extracellular cues.

Ligand‐dependent responses of the ErbB signaling network: experimental and modeling analyses

Deregulation of ErbB signaling plays a key role in the progression of multiple human cancers. To help understand ErbB signaling quantitatively, in this work we combine traditional experiments with computational modeling, building a model that describes how stimulation of all four ErbB receptors with epidermal growth factor (EGF) and heregulin (HRG) leads to activation of two critical downstream proteins, extracellular‐signal‐regulated kinase (ERK) and Akt. Model analysis and experimental validation show that (i) ErbB2 overexpression, which occurs in approximately 25% of all breast cancers, transforms transient EGF‐induced signaling into sustained signaling, (ii) HRG‐induced ERK activity is much more robust to the ERK cascade inhibitor U0126 than EGF‐induced ERK activity, and (iii) phosphoinositol‐3 kinase is a major regulator of post‐peak but not pre‐peak EGF‐induced ERK activity. Sensitivity analysis leads to the hypothesis that ERK activation is robust to parameter perturbation at high ligand doses, while Akt activation is not.

The Mammalian MAPK/ERK Pathway Exhibits Properties of a Negative Feedback Amplifier

Analysis of ERK pathway circuitry suggests appropriate targets for inhibition, providing a guide for drug development. Biological Circuits Inform Drug Development The mitogen-activated protein kinase (MAPK) pathway involves a three-tiered kinase module, which amplifies the signal. Many cells also have negative feedback loops from the last kinase in the module to various points upstream in the pathway. Sturm et al. showed that, with negative feedback loops, the MAPK module results in a system like that of a negative feedback amplifier (NFA), which is an engineering design that smoothens the output to changes in input and makes a system robust to change. These NFA-like properties may explain why some cells are sensitive to inhibition of the second kinase in the cascade (they lack feedback loops), whereas other cells are resistant to inhibition at this point (their feedback loops are intact). These results also have implications for drug development, because inhibitors that target components that are outside the NFA are more effective at inhibiting the pathway. Three-tiered kinase modules, such as the Raf–MEK (mitogen-activated or extracellular signal–regulated protein kinase kinase)–ERK (extracellular signal–regulated kinase) mitogen-activated protein kinase pathway, are widespread in biology, suggesting that this structure conveys evolutionarily advantageous properties. We show that the three-tiered kinase amplifier module combined with negative feedback recapitulates the design principles of a negative feedback amplifier (NFA), which is used in electronic circuits to confer robustness, output stabilization, and linearization of nonlinear signal amplification. We used mathematical modeling and experimental validation to demonstrate that the ERK pathway has properties of an NFA that (i) converts intrinsic switch-like activation kinetics into graded linear responses, (ii) conveys robustness to changes in rates of reactions within the NFA module, and (iii) stabilizes outputs in response to drug-induced perturbations of the amplifier. These properties determine biological behavior, including activation kinetics and the response to drugs.

Systems-level interactions between insulin-EGF networks amplify mitogenic signaling.

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 (PIP3) concentration into enhanced Ras/ERK activity. The model predicts and experiments confirm that insulin‐induced amplification of mitogenic signaling is abolished by disrupting PIP3‐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.

SPATIAL GRADIENTS OF CELLULAR PHOSPHO-PROTEINS

If a protein is rapidly phosphorylated and dephosphorylated at separate cellular locations and protein diffusion is slow, then a spatial gradient of the phosphorylated form of the protein may develop within the cell. We have estimated the potential size of such gradients using measured values of protein diffusion coefficients and protein kinase and phosphatase activities. We analysed two different cellular geometries: (1) where the kinases is located on the plasma membrane of a spherical cell and the phospatase is distributed homogenously in the cytoplasm and (2) where the kinase is located on one planar membrane and the phosphatase on a second parallel planar membrane. The estimated gradients of phospho‐proteins were potentially very large, which has important implications for cellular signalling.