Now showing 1 - 10 of 18
  • 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.
      371Scopus© Citations 144
  • Publication
    PP1 initiates the dephosphorylation of MASTL, triggering mitotic exit and bistability in human cells
    (The Company of Biologists, 2016-04-01) ; ; ;
    Entry into mitosis is driven by the phosphorylation of thousands of substrates, under the master control of Cdk1. During entry into mitosis, Cdk1, in collaboration with MASTL kinase, represses the activity of the major mitotic protein phosphatases, PP1 and PP2A, thereby ensuring mitotic substrates remain phosphorylated. For cells to complete and exit mitosis, these phosphorylation events must be removed, and hence, phosphatase activity must be reactivated. This reactivation of phosphatase activity presumably requires the inhibition of MASTL; however, it is not currently understood what deactivates MASTL and how this is achieved. In this study, we identified that PP1 is associated with, and capable of partially dephosphorylating and deactivating, MASTL during mitotic exit. Using mathematical modelling, we were able to confirm that deactivation of MASTL is essential for mitotic exit. Furthermore, small decreases in Cdk1 activity during metaphase are sufficient to initiate the reactivation of PP1, which in turn partially deactivates MASTL to release inhibition of PP2A and, hence, create a feedback loop. This feedback loop drives complete deactivation of MASTL, ensuring a strong switch-like activation of phosphatase activity during mitotic exit.
      320Scopus© Citations 40
  • Publication
    Pseudophosphatase STYX modulates cell-fate decisions and cell migration by spatiotemporal regulation of ERK1/2
    (National Academy of Sciences, 2013-07-11) ; ; ;
    Serine/threonine/tyrosine-interacting protein (STYX) is a catalytically inactive member of the dual-specificity phosphatases (DUSPs) family. Whereas the role of DUSPs in cellular signaling is well explored, the function of STYX is still unknown. Here, we identify STYX as a spatial regulator of ERK signaling. We used predictive-model simulation to test several hypotheses for possible modes of STYX action. We show that STYX localizes to the nucleus, competes with nuclear DUSP4 for binding to ERK, and acts as a nuclear anchor that regulates ERK nuclear export. Depletion of STYX increases ERK activity in both cytosol and nucleus. Importantly, depletion of STYX causes an ERK-dependent fragmentation of the Golgi apparatus and inhibits Golgi polarization and directional cell migration. Finally, we show that overexpression of STYX reduces ERK1/2 activation, thereby blocking PC12 cell differentiation. Overall, our results identify STYX as an important regulator of ERK1/2 signaling critical for cell migration and PC12 cell differentiation.
      524Scopus© Citations 45
  • 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.
      20Scopus© Citations 11
  • 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
    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.  
  • 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.
      340Scopus© Citations 10
  • Publication
    Accurate prediction of kinase-substrate networks using knowledge graphs
    Phosphorylation of specific substrates by protein kinases is a key control mechanism for vital cell-fate decisions and other cellular processes. However, discovering specific kinasesubstrate relationships is time-consuming and often rather serendipitous. Computational predictions alleviate these challenges, but the current approaches suffer from limitations like restricted kinome coverage and inaccuracy. They also typically utilise only local features without reflecting broader interaction context. To address these limitations, we have developed an alternative predictive model. It uses statistical relational learning on top of phosphorylation networks interpreted as knowledge graphs, a simple yet robust model for representing networked knowledge. Compared to a representative selection of six existing systems, our model has the highest kinome coverage and produces biologically valid highconfidence predictions not possible with the other tools. Specifically, we have experimentally validated predictions of previously unknown phosphorylations by the LATS1, AKT1, PKA and MST2 kinases in human. Thus, our tool is useful for focusing phosphoproteomic experiments, and facilitates the discovery of new phosphorylation reactions. Our model can be accessed publicly via an easy-to-use web interface (LinkPhinder).
      20Scopus© Citations 12
  • 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.
      463Scopus© Citations 13
  • Publication
    Multiscale Model of Dynamic Neuromodulation Integrating Neuropeptide-Induced Signaling Pathway Activity with Membrane Electrophysiology
    We developed a multiscale model to bridge neuropeptide receptor-activated signaling pathway activity with membrane electrophysiology. Typically, the neuromodulation of biochemical signaling and biophysics have been investigated separately in modeling studies. We studied the effects of Angiotensin II (AngII) on neuronal excitability changes mediated by signaling dynamics and downstream phosphorylation of ion channels. Experiments have shown that AngII binding to the AngII receptor type-1 elicits baseline-dependent regulation of cytosolic Ca2+ signaling. Our model simulations revealed a baseline Ca2+-dependent response to AngII receptor type-1 activation by AngII. Consistent with experimental observations, AngII evoked a rise in Ca2+ when starting at a low baseline Ca2+ level, and a decrease in Ca2+ when starting at a higher baseline. Our analysis predicted that the kinetics of Ca2+ transport into the endoplasmic reticulum play a critical role in shaping the Ca2+ response. The Ca2+ baseline also influenced the AngII-induced excitability changes such that lower Ca2+ levels were associated with a larger firing rate increase. We examined the relative contributions of signaling kinases protein kinase C and Ca2+/Calmodulin-dependent protein kinase II to AngII-mediated excitability changes by simulating activity blockade individually and in combination. We found that protein kinase C selectively controlled firing rate adaptation whereas Ca2+/Calmodulin-dependent protein kinase II induced a delayed effect on the firing rate increase. We tested whether signaling kinetics were necessary for the dynamic effects of AngII on excitability by simulating three scenarios of AngII-mediated KDR channel phosphorylation: (1), an increased steady state; (2), a step-change increase; and (3), dynamic modulation. Our results revealed that the kinetics emerging from neuromodulatory activation of the signaling network were required to account for the dynamical changes in excitability. In summary, our integrated multiscale model provides, to our knowledge, a new approach for quantitative investigation of neuromodulatory effects on signaling and electrophysiology.
      345Scopus© Citations 3