Now showing 1 - 10 of 18
  • Publication
    Crosstalk and Signaling Switches in Mitogen-Activated Protein Kinase Cascades
    Mitogen-activated protein kinase (MAPK) cascades control cell fate decisions, such as proliferation, differentiation, and apoptosis by integrating and processing intra- and extracellular cues. However, similar MAPK kinetic profiles can be associated with opposing cellular decisions depending on cell type, signal strength, and dynamics. This implies that signaling by each individual MAPK cascade has to be considered in the context of the entire MAPK network. Here, we develop a dynamic model of feedback and crosstalk for the three major MAPK cascades; extracellular signal-regulated kinase (ERK), p38 mitogen-activated protein kinase (p38), c-Jun N-terminal kinase (JNK), and also include input from protein kinase B (AKT) signaling. Focusing on the bistable activation characteristics of the JNK pathway, this model explains how pathway crosstalk harmonizes different MAPK responses resulting in pivotal cell fate decisions. We show that JNK can switch from a transient to sustained activity due to multiple positive feedback loops. Once activated, positive feedback locks JNK in a highly active state and promotes cell death. The switch is modulated by the ERK, p38, and AKT pathways. ERK activation enhances the dual specificity phosphatase (DUSP) mediated dephosphorylation of JNK and shifts the threshold of the apoptotic switch to higher inputs. Activation of p38 restores the threshold by inhibiting ERK activity via the PP1 or PP2A phosphatases. Finally, AKT activation inhibits the JNK positive feedback, thus abrogating the apoptotic switch and allowing only proliferative signaling. Our model facilitates understanding of how cancerous deregulations disturb MAPK signal processing and provides explanations for certain drug resistances. We highlight a critical role of DUSP1 and DUSP2 expression patterns in facilitating the switching of JNK activity and show how oncogene induced ERK hyperactivity prevents the normal apoptotic switch explaining the failure of certain drugs to induce apoptosis.
    Scopus© Citations 135  388
  • Publication
    BAX and SMAC regulate bistable properties of the apoptotic caspase system
    The initiation of apoptosis is a core mechanism in cellular biology by which organisms control the removal of damaged or unnecessary cells. The irreversible activation of caspases is essential for apoptosis, and mathematical models have demonstrated that the process is tightly regulated by positive feedback and a bistable switch. BAX and SMAC are often dysregulated in diseases such as cancer or neurodegeneration and are two key regulators that interact with the caspase system generating the apoptotic switch. Here we present a mathematical model of how BAX and SMAC control the apoptotic switch. Formulated as a system of ordinary differential equations, the model summarises experimental and computational evidence from the literature and incorporates the biochemical mechanisms of how BAX and SMAC interact with the components of the caspase system. Using simulations and bifurcation analysis, we find that both BAX and SMAC regulate the time-delay and activation threshold of the apoptotic switch. Interestingly, the model predicted that BAX (not SMAC) controls the amplitude of the apoptotic switch. Cell culture experiments using siRNA mediated BAX and SMAC knockdowns validated this model prediction. We further validated the model using data of the NCI-60 cell line panel using BAX protein expression as a cell-line specific parameter and show that model simulations correlated with the cellular response to DNA damaging drugs and established a defined threshold for caspase activation that could distinguish between sensitive and resistant melanoma cells. In summary, we present an experimentally validated dynamic model that summarises our current knowledge of how BAX and SMAC regulate the bistable properties of irreversible caspase activation during apoptosis.
      16Scopus© Citations 11
  • Publication
    A dynamic model of the MYCN regulated DNA damage response in Neuroblastoma
    Neuroblastoma is the most common the most common cancer in infancy with an extremely heterogeneous phenotype that is mainly driven by the MYCN oncogene. The MYCN transcription factor and its amplification is commonly associated with poor prognosis in patients, although it has also been shown that elevated MYCN levels correlates with apoptosis sensitization in cells. HMGA1 is one of MYCN target genes and is involved in triggering apoptosis through a DNA Damage Response (DDR) by inducing ataxia-telangiectasia-mutated (ATM) gene expression. But HMGA1 is also involved in preventing apoptosis by directly binding HIPK2 and decreasing its presence in the nucleus, therefore decreasing phosphorylation of p53 at serine 46 which is required for the activation of p53 apoptotic targets. In this article, we propose a model in which MYCN protein regulates the HMGA1-ATM-p53 and HMGA1-HIPK2-p53 subsystems. Because the molecular details concerning the HMGA1-HMGA1 interaction are uncertain several possibilities were explored in simulations. Our model points towards an important role of MYCN-dependent regulation of HMGA1 expression levels and the subsequent HIPK2 nuclear/cytoplasmic re-localization and led to experimentally testable predictions that can discern between alternative model structures.  
      187
  • Publication
    Signaling pathway models as biomarkers: Patient-specific simulations of JNK activity predict the survival of neuroblastoma patients
    Signaling pathways control cell fate decisions that ultimately determine the behavior of cancer cells. Therefore, the dynamics of pathway activity may contain prognostically relevant information different from that contained in the static nature of other types of biomarkers. To investigate this hypothesis, we characterized the network that regulated stress signaling by the c-Jun N-terminal kinase (JNK) pathway in neuroblastoma cells. We generated an experimentally calibrated and validated computational model of this network and used the model to extract prognostic information from neuroblastoma patient–specific simulations of JNK activation. Switch-like JNK activation mediates cell death by apoptosis. An inability to initiate switch-like JNK activation in the simulations was significantly associated with poor overall survival for patients with neuroblastoma with or without MYCN amplification, indicating that patient-specific simulations of JNK activation could stratify patients. Furthermore, our analysis demonstrated that extracting information about a signaling pathway to develop a prognostically useful model requires understanding of not only components and disease-associated changes in the abundance or activity of the components but also how those changes affect pathway dynamics.
      1349Scopus© Citations 109
  • Publication
    Applications of personalised signalling network models in precision oncology
    As our ability to provide in-depth, patient-specific characterisation of the molecular alterations within tumours rapidly improves, it is becoming apparent that new approaches will be required to leverage the power of this data and derive the full benefit for each individual patient. Systems biology approaches are beginning to emerge within this field as a potential method of incorporating large volumes of network level data and distilling a coherent, clinically-relevant prediction of drug response. However, the initial promise of this developing field is yet to be realised. Here we argue that in order to develop these precise models of individual drug response and tailor treatment accordingly, we will need to develop mathematical models capable of capturing both the dynamic nature of drug-response signalling networks and key patient-specific information such as mutation status or expression profiles. We also review the modelling approaches commonly utilised within this field, and outline recent examples of their use in furthering the application of systems biology for a precision medicine approach to cancer treatment.
    Scopus© Citations 13  457
  • Publication
    Frequency modulation of ERK activation dynamics rewires cell fate
    Transient versus sustained ERK MAP kinase (MAPK) activation dynamics induce proliferation versus differentiation in response to epidermal (EGF) or nerve (NGF) growth factors in PC-12 cells. Duration of ERK activation has therefore been proposed to specify cell fate decisions. Using a biosensor to measure ERK activation dynamics in single living cells reveals that sustained EGF/NGF application leads to a heterogeneous mix of transient and sustained ERK activation dynamics in distinct cells of the population, different than the population average. EGF biases toward transient, while NGF biases toward sustained ERK activation responses. In contrast, pulsed growth factor application can repeatedly and homogeneously trigger ERK activity transients across the cell population. These datasets enable mathematical modeling to reveal salient features inherent to the MAPK network. Ultimately, this predicts pulsed growth factor stimulation regimes that can bypass the typical feedback activation to rewire the system toward cell differentiation irrespective of growth factor identity.
    Scopus© Citations 144  364
  • Publication
    Impact of measurement noise, experimental design, and estimation methods on Modular Response Analysis based network reconstruction
    Modular Response Analysis (MRA) is a method to reconstruct signalling networks from steady-state perturbation data which has frequently been used in different settings. Since these data are usually noisy due to multi-step measurement procedures and biological variability, it is important to investigate the effect of this noise onto network reconstruction. Here we present a systematic study to investigate propagation of noise from concentration measurements to network structures. Therefore, we design an in silico study of the MAPK and the p53 signalling pathways with realistic noise settings. We make use of statistical concepts and measures to evaluate accuracy and precision of individual inferred interactions and resulting network structures. Our results allow to derive clear recommendations to optimize the performance of MRA based network reconstruction: First, large perturbations are favorable in terms of accuracy even for models with non-linear steady-state response curves. Second, a single control measurement for different perturbation experiments seems to be sufficient for network reconstruction, and third, we recommend to execute the MRA workflow with the mean of different replicates for concentration measurements rather than using computationally more involved regression strategies.
      387Scopus© Citations 9
  • Publication
    On the personalised modelling of cancer signalling
    Dynamic modelling has long been used to understand fundamental principles of cell signalling and its dysregulation in cancer. More recently these models have also been used to understand the individual risks of cancer patients, and predict their survival probabilities. However, the current methodologies for integrating tumour data and generating patient-specific simulations suffer from the lack of general applicability; they only work for cell signalling models in which only posttranslational protein modifications are considered, so that the total protein concentrations are conserved. Here, we present novel, generally applicable method. The method is based on a simple theoretical framework for modelling gene-regulation, and the indirect estimation of patient-specific parameters from tumour data. Because our method does not require time-invariance of the total-protein concentrations, it can be applied to models of any nature, including the many cancer signalling models involving gene-regulation.
      325Scopus© Citations 2
  • Publication
    Polyubiquitin chain assembly and organization determine the dynamics of protein activation and degradation
    Protein degradation via ubiquitination is a major proteolytic mechanism in cells. Once a protein is destined for degradation, it is tagged by multiple ubiquitin (Ub) molecules. The synthesized polyubiquitin chains can be recognized by the 26S proteosome where proteins are degraded. These chains form through multiple ubiquitination cycles that are similar to multi-site phosphorylation cycles. As kinases and phosphatases, two opposing enzymes (E3 ligases and deubiquitinases DUBs) catalyze (de)ubiquitination cycles. Although multi-ubiquitination cycles are fundamental mechanisms of controlling protein concentrations within a cell, their dynamics have never been explored. Here, we fill this knowledge gap. We show that under permissive physiological conditions, the formation of polyubiquitin chain of length greater than two and subsequent degradation of the ubiquitinated protein, which is balanced by protein synthesis, can display bistable, switch-like responses. Interestingly, the occurrence of bistability becomes pronounced, as the chain grows, giving rise to "all-or-none" regulation at the protein levels. We give predictions of protein distributions under bistable regime awaiting experimental verification. Importantly, we show for the first time that sustained oscillations can robustly arise in the process of formation of ubiquitin chain, largely due to the degradation of the target protein. This new feature is opposite to the properties of multi-site phosphorylation cycles, which are incapable of generating oscillation if the total abundance of interconverted protein forms is conserved. We derive structural and kinetic constraints for the emergence of oscillations, indicating that a competition between different substrate forms and the E3 and DUB is critical for oscillation. Our work provides the first detailed elucidation of the dynamical features brought about by different molecular setups of the polyubiquitin chain assembly process responsible for protein degradation.
      323Scopus© Citations 24
  • Publication
    Hippocampal contribution to vector model hypothesis during cue-dependent navigation
    (Cold Spring Harbor Laboratory Press, 2013-06-17) ; ;
    Learning to navigate toward a goal is an essential skill. Place learning is thought to rely on the ability of animals to associate the location of a goal with surrounding environmental cues. Using the Morris water maze, a task popularly used to examine place learning, we demonstrate that distal cues provide animals with distance and directional information. We show how animals use the cues in a visually dependent guidance manner to find the goal. Further, we demonstrate how hippocampal lesions disrupt this learning mechanism. Our results can be explained through the vector model of navigation built on associative learning principles rather than evoking a cognitive map.
    Scopus© Citations 10  340