Data-driven modeling reconciles kinetics of ERK phosphorylation, localization, and activity states

A guide to ERK dynamics, part 1: mechanisms and models

Extracellular signal-regulated kinase (ERK) has long been studied as a key driver of both essential cellular processes and disease. A persistent question has been how this single pathway is able to direct multiple cell behaviors, including growth, proliferation, and death. Modern biosensor studies have revealed that the temporal pattern of ERK activity is highly variable and heterogeneous, and critically, that these dynamic differences modulate cell fate. This two-part review discusses the current understanding of dynamic activity in the ERK pathway, how it regulates cellular decisions, and how these cell fates lead to tissue regulation and pathology. In part 1, we cover the optogenetic and live-cell imaging technologies that first revealed the dynamic nature of ERK, as well as current challenges in biosensor data analysis. We also discuss advances in mathematical models for the mechanisms of ERK dynamics, including receptor-level regulation, negative feedback, cooperativity, and paracrine signaling. While hurdles still remain, it is clear that higher temporal and spatial resolution provide mechanistic insights into pathway circuitry. Exciting new algorithms and advanced computational tools enable quantitative measurements of single-cell ERK activation, which in turn inform better models of pathway behavior. However, the fact that current models still cannot fully recapitulate the diversity of ERK responses calls for a deeper understanding of network structure and signal transduction in general.

Subcellular partitioning of protein kinase activity revealed by functional kinome profiling

Protein kinases and their substrates form signaling networks partitioned across subcellular compartments to facilitate critical biological processes. While the subcellular roles of many individual kinases have been elucidated, a comprehensive assessment of the synaptic subkinome is lacking. Further, most studies of kinases focus on transcript, protein, and/or phospho-protein expression levels, providing an indirect measure of protein kinase activity. Prior work suggests that gene expression levels are not a good predictor of protein function. Thus, we assessed global serine/threonine protein kinase activity profiles in synaptosomal, nuclear, and cytosolic fractions from rat frontal cortex homogenate using peptide arrays. Comparisons made between fractions demonstrated differences in overall protein kinase activity. Upstream kinase analysis revealed a list of cognate kinases that were enriched in the synaptosomal fraction compared to the nuclear fraction. We identified many kinases in the synaptic fraction previously implicated in this compartment, while also identifying other kinases with little or no evidence for synaptic localization. Our results show the feasibility of assessing subcellular fractions with peptide activity arrays, as well as suggesting compartment specific activity profiles associated with established and novel kinases.

Regulation of STAT3 signaling in IFNγ and IL10 pathways and in their cross-talk

The pro-inflammatory IFNγ-STAT1 pathway and anti-inflammatory IL10-STAT3 pathway elicit cellular responses primarily utilizing their canonical STATs. However IL10 mediated STAT1 and IFNγ mediated STAT3 activation is also observed, suggesting crosstalk of these functionally opposing signaling pathways can potentially reshape the canonical dynamics both STATs and alter the expression of their target genes. Herein, we measured the dynamics of STATs in response to different doses of IL10 or IFNγ and in their co-stimulation and employed quantitative modeling to understand the regulatory mechanisms controlling signal responses in individual and co-simulation scenarios. Our experiments show, STAT3 in particular, exhibits a bell-shaped dose-response while treated with IFNγ or IL10 and our model quantiatively captured the dose-dependent dynamics of both the STATs in both pathways. The model next predicted and subsequent experiments validated that STAT3 dynamics would robustly remain IL10 specific when subjected to a co-stimulation of both IFNγ and IL10. Genes common to both pathways also exhibited IL10 specific expression during the co-stimulation. The findings thus uncover anovel feature of the IL10-STAT3 signaling axis during pathway crosstalk. Finally, parameter sampling coupled to information theory based analysis showed that bell-shaped signal-response of STAT3 in both pathways is primarily dependent on receptor concentration whereas robustness of IL10-STAT3 signaling axis in co-stimulation results from the negative regulation of the IFNγ pathway.

Cell geometry, signal dampening, and a bimodal transcriptional response underlie the spatial precision of an ERK-mediated embryonic induction

Precise control of lineage segregation is critical for the development of multicellular organisms, but our quantitative understanding of how variable signaling inputs are integrated to activate lineage-specific gene programs remains limited. Here, we show how precisely two out of eight ectoderm cells adopt neural fates in response to ephrin and FGF signals during ascidian neural induction. In each ectoderm cell, FGF signals activate ERK to a level that mirrors its cell contact surface with FGF-expressing mesendoderm cells. This gradual interpretation of FGF inputs is followed by a bimodal transcriptional response of the immediate early gene, Otx, resulting in its activation specifically in the neural precursors. At low levels of ERK, Otx is repressed by an ETS family transcriptional repressor, ERF2. Ephrin signals are critical for dampening ERK activation levels across ectoderm cells so that only neural precursors exhibit above-threshold levels, evade ERF repression, and "switch on" Otx transcription.

Signaling Pathways in Proton and Non-proton ASIC1a Activation

Acid-sensing ion channels (ASICs) regulate synaptic activities and play important roles in neurodegenerative diseases as well as pain conditions. Classically, ASICs are described as transiently activated by a reduced pH, followed by desensitization; the activation allows sodium influx, and in the case of ASIC1a-composed channels, also calcium to some degree. Several factors are emerging and extensively analyzed as modulators, activating, inhibiting, and potentiating specific channel subunits. However, the signaling pathways triggered by channel activation are only starting to be revealed.The channel has been recently shown to be activated through a mechanism other than proton-mediated. Indeed, the large extracellular loop of these channels opens the possibility that other non-proton ligands might exist. One such molecule discovered was a toxin present in the Texas coral snake venom. The finding was associated with the activation of the channel at neutral pH via the toxin and causing intense and unremitting pain.By using different pharmacological tools, we analyzed the downstream signaling pathway triggered either by the proton and non-proton activation for human, mouse, and rat ASIC1a-composed channels in in vitro models. We show that for all species analyzed, the non-protonic mode of activation determines the activation of the ERK signaling cascade at a higher level and duration compared to the proton mode.This study adds to the growing evidence of the important role ASIC1a channels play in different physiological and pathological conditions and also hints at a possible pathological mechanism for a sustained effect.

Hexokinase 2 Promotes Cell Growth and Tumor Formation Through the Raf/MEK/ERK Signaling Pathway in Cervical Cancer

Hexokinase 2 (HK2) is a member of the hexokinases (HK) that has been reported to be a key regulator during glucose metabolism linked to malignant growth in many types of cancers. In this study, stimulation of HK2 expression was observed in squamous cervical cancer (SCC) tissues, and HK2 expression promoted the proliferation of cervical cancer cells in vitro and tumor formation in vivo by accelerating cell cycle progression, upregulating cyclin A1, and downregulating p27 expression. Moreover, transcriptome sequencing analysis revealed that MAPK3 (ERK1) was upregulated in HK2-overexpressing HeLa cells. Further experiments found that the protein levels of p-Raf, p-MEK1/2, ERK1/2, and p-ERK1/2 were increased in HK2 over-expressing SiHa and HeLa cells. When ERK1/2 and p-ERK1/2 expression was blocked by an inhibitor (FR180204), reduced cyclin A1 expression was observed in HK2 over-expressing cells, with induced p27 expression and inhibited cell growth. Therefore, our data demonstrated that HK2 promoted the proliferation of cervical cancer cells by upregulating cyclin A1 and down-regulating p27 expression through the Raf/MEK/ERK signaling pathway.

Optical control of MAP kinase kinase 6 (MKK6) reveals that it has divergent roles in pro-apoptotic and anti-proliferative signaling

The selective pressure imposed by extrinsic death signals and stressors adds to the challenge of isolating and interpreting the roles of proteins in stress-activated signaling networks. By expressing a kinase with activating mutations and a caged lysine blocking the active site, we can rapidly switch on catalytic activity with light and monitor the ensuing dynamics. Applying this approach to MAP kinase 6 (MKK6), which activates the p38 subfamily of MAPKs, we found that decaging active MKK6 in fibroblasts is sufficient to trigger apoptosis in a p38-dependent manner. Both in fibroblasts and in a murine melanoma cell line expressing mutant B-Raf, MKK6 activation rapidly and potently inhibited the pro-proliferative extracellular signal-regulated kinase (ERK) pathway; to our surprise, this negative cross-regulation was equally robust when all p38 isoforms were inhibited. These results position MKK6 as a new pleiotropic signal transducer that promotes both pro-apoptotic and anti-proliferative signaling, and they highlight the utility of caged, light-activated kinases for dissecting stress-activated signaling networks.

Optical control of protein phosphatase function

Protein phosphatases are involved in embryonic development, metabolic homeostasis, stress response, cell cycle transitions, and many other essential biological mechanisms. Unlike kinases, protein phosphatases remain understudied and less characterized. Traditional genetic and biochemical methods have contributed significantly to our understanding; however, these methodologies lack precise and acute spatiotemporal control. Here, we report the development of a light-activated protein phosphatase, the dual specificity phosphatase 6 (DUSP6 or MKP3). Through genetic code expansion, MKP3 is placed under optical control via two different approaches: (i) incorporation of a caged cysteine into the active site for controlling catalytic activity and (ii) incorporation of a caged lysine into the kinase interaction motif for controlling the protein-protein interaction between the phosphatase and its substrate. Both strategies are expected to be applicable to the engineering of a wide range of light-activated phosphatases. Applying the optogenetically controlled MKP3 in conjunction with live cell reporters, we discover that ERK nuclear translocation is regulated in a graded manner in response to increasing MKP3 activity.

Genetically Encoded Fluorescent Biosensors Illuminate the Spatiotemporal Regulation of Signaling Networks

Cellular signaling networks are the foundation which determines the fate and function of cells as they respond to various cues and stimuli. The discovery of fluorescent proteins over 25 years ago enabled the development of a diverse array of genetically encodable fluorescent biosensors that are capable of measuring the spatiotemporal dynamics of signal transduction pathways in live cells. In an effort to encapsulate the breadth over which fluorescent biosensors have expanded, we endeavored to assemble a comprehensive list of published engineered biosensors, and we discuss many of the molecular designs utilized in their development. Then, we review how the high temporal and spatial resolution afforded by fluorescent biosensors has aided our understanding of the spatiotemporal regulation of signaling networks at the cellular and subcellular level. Finally, we highlight some emerging areas of research in both biosensor design and applications that are on the forefront of biosensor development.

Generation of recombinant affinity reagents against a two-phosphosite epitope of ATF2

Activating Transcription Factor 2 (ATF2) plays an important role in mammalian cell proliferation, apoptosis and DNA repair. Its activation is dependent on the sequential phosphorylation of residue threonine 71 (T71) followed by threonine 69 (T69) in its transactivation domain. While these modifications can be directed by a variety of kinases, the time to reach full phosphorylation is dependent on which signaling pathway has been activated, which is thought to be important for proper temporal regulation. To explore this phenomenon further, there have been ongoing efforts to generate affinity reagents for monitoring phosphorylation events in cellular assays. While phospho-specific antibodies have been valuable tools for monitoring cell signaling events, those raised against a peptide containing two or more adjacent phosphosites tend to cross-react with that peptide's various phospho-states, rendering such reagents unusable for studying sequential phosphorylation. As an alternative, we have employed the N-terminal Forkhead-associated 1 (FHA1) domain of yeast Rad53p as a scaffold to generate recombinant affinity reagents via phage display and were successful in generating a set of reagents that can distinguish between the dual-phosphorylated epitope, 63-IVADQpTPpTPTRFLK-77, and the mono-phosphorylated epitope, 63-IVADQpTPTPTRFLK-77, in the human ATF2 transactivation domain.

Importance of Feedback and Feedforward Loops to Adaptive Immune Response Modeling

The human adaptive immune system is a very complex network of different types of cells, cytokines, and signaling molecules. This complex network makes it difficult to understand the system level regulations. To properly explain the immune system, it is necessary to explicitly investigate the presence of different feedback and feedforward loops (FFLs) and their crosstalks. Considering that these loops increase the complexity of the system, the mathematical modeling has been proved to be an important tool to explain such complex biological systems. This review focuses on these regulatory loops and discusses their importance on systems modeling of the immune system.

Signaling pathways as linear transmitters

One challenge in biology is to make sense of the complexity of biological networks. A good system to approach this is signaling pathways, whose well-characterized molecular details allow us to relate the internal processes of each pathway to their input-output behavior. In this study, we analyzed mathematical models of three metazoan signaling pathways: the canonical Wnt, MAPK/ERK, and Tgfβ pathways. We find an unexpected convergence: the three pathways behave in some physiological contexts as linear signal transmitters. Testing the results experimentally, we present direct measurements of linear input-output behavior in the Wnt and ERK pathways. Analytics from each model further reveal that linearity arises through different means in each pathway, which we tested experimentally in the Wnt and ERK pathways. Linearity is a desired property in engineering where it facilitates fidelity and superposition in signal transmission. Our findings illustrate how cells tune different complex networks to converge on the same behavior.

A novel arylbenzofuran induces cervical cancer cell apoptosis and G1/S arrest through ERK-mediated Cdk2/cyclin-A signaling pathway

7-hydroxy-5,4'-dimethoxy-2-arylbenzofuran (Ary) is purified from Livistona. It has been demonstrated to have anticancer activity to various tumors in including cervical cancer, but its mechanism is still unclear. In the present, we show that Ary induces cervical cancer cells apoptosis through mitochondria degradation and mediates cervical cancer cell arrest. Further, Ary-inducing cell cycle G1/S-phase arrest is associated with increased cyclin A2 and cyclin dependent kinase 2 (Cdk2) proteins. Knockdown of cyclin A2 using small interfering RNA (siRNA), and inhibiting Cdk2 activity with flavopiridol, strikingly reduced G1/S-phase arrest. Moreover, Ary sustainedly induced phosphorylation of extracellular signal-regulated kinase1/2 (ERK1/2). And ERK1/2 phosphorylation inhibition using specific inhibitor U0126 effectively suppressed cyclin A2 expression, and reduced G1/S-phase arrest induced by Ary. All the experiments in vitro and in vivo verified that Ary has an anticancer effect on cervical cancer. These data provide novel evidences that Ary induces cervical cancer cells apoptosis through mitochondria degradation and cell G1/S-phase arrest. These findings also suggest that ERK-mediated Cdk2/cyclin A signaling pathway is involved in Ary-induced G1/S-phase arrest.

Use of mathematics to guide target selection in systems pharmacology; application to receptor tyrosine kinase (RTK) pathways

A key element of the drug discovery process is target selection. Although the topic is subject to much discussion and experimental effort, there are no defined quantitative rules around optimal selection. Often 'rules of thumb', that have not been subject to rigorous exploration, are used. In this paper we explore the 'rule of thumb' notion that the molecule that initiates a pathway signal is the optimal target. Given the multi-factorial and complex nature of this question, we have simplified an example pathway to its logical minimum of two steps and used a mathematical model of this to explore the different options in the context of typical small and large molecule drugs. In this paper, we report the conclusions of our analysis and describe the analysis tool and methods used. These provide a platform to enable a more extensive enquiry into this important topic.

Identification of PTPN1 as a novel negative regulator of the JNK MAPK pathway using a synthetic screening for pathway-specific phosphatases

The mitogen activated protein kinase (MAPK) signaling cascades transmit extracellular stimulations to generate various cellular responses via the sequential and reversible phosphorylation of kinases. Since the strength and duration of kinase phosphorylation within the pathway determine the cellular response, both kinases and phosphatases play an essential role in the precise control of MAPK pathway activation and attenuation. Thus, the identification of pathway-specific phosphatases is critical for understanding the functional mechanisms by which the MAPK pathway is regulated. To identify phosphatases specific to the c-Jun N-terminal kinase (JNK) MAPK pathway, a synthetic screening approach was utilized in which phosphatases were individually tethered to the JNK pathway specific-JIP1 scaffold protein. Of 77 mammalian phosphatases tested, PTPN1 led to the inhibition of JNK pathway activation. Further analyses revealed that of three pathway member kinases, PTPN1 directly dephosphorylates JNK, the terminal kinase of the pathway, and negatively regulates the JNK MAPK pathway. Specifically, PTPN1 appears to regulate the overall signaling magnitude, rather than the adaptation timing, suggesting that PTPN1 might be involved in the control and maintenance of signaling noise. Finally, the negative regulation of the JNK MAPK pathway by PTPN1 was found to reduce the tumor necrosis factor α (TNFα)-dependent cell death response.

Two faces of competition: target-mediated reverse signalling in microRNA and mitogen-activated protein kinase regulatory networks

Biomolecular regulatory networks are organised around hubs, which can interact with a large number of targets. These targets compete with each other for access to their common hubs, but whether the effect of this competition would be significant in magnitude and in function is not clear. In this review, the authors discuss recent in vivo studies that analysed the system level retroactive effects induced by target competition in microRNA and mitogen-activated protein kinase regulatory networks. The results of these studies suggest that downstream targets can regulate the overall state of their upstream regulators, and thus cannot be ignored in analysing biomolecular networks.

Modulation of Glutathione Hemostasis by Inhibition of 12/15-Lipoxygenase Prevents ROS-Mediated Cell Death after Hepatic Ischemia and Reperfusion

Background

Reactive oxygen species- (ROS-) mediated ischemia-reperfusion injury (IRI) detrimentally impacts liver transplantation and resection. 12/15-Lipoxygenase (12/15-LOX), an antagonistic protein of the glutathione peroxidase 4 (GPX4) signaling cascade, was proven to mediate cell death in postischemic cerebral and myocardial tissue. The aim of this study was to investigate the impact of 12/15-LOX inhibition on hepatic IRI.

Methods

Livers of C57BL/6 mice were exposed to 60 minutes of partial warm ischemia and 90 minutes of reperfusion after previous Baicalein administration, an inhibitor of 12/15-LOX. Tissue samples were analyzed by TUNEL assay, Western blot, and spectral photometry.

Results

TUNEL labeling showed a significant reduction of hepatic cell death following baicalein pretreatment. Western Blot analysis revealed a significant downregulation of Jun-amino-terminal-kinase (JNK), caspase-3, and poly-ADP-ribose-polymerase (PARP), besides considerably lowered p44/42-MAP-kinase (ERK1/2) expression after Baicalein administration. A significant elevation of glutathione oxidation was measured in Baicalein pretreated livers.

Conclusion

Our data show that inhibition of 12/15-lipoxygenase causes significant cell death reduction after hepatic ischemia and reperfusion by enhancing glutathione metabolism. We conclude that GPX4-dependent cell death signaling cascade might play a major role in development of hepatic IRI, in which the investigated proteins JNK, caspase-3, ERK1/2, and PARP might contribute to tissue damage.

Sampling-based Bayesian approaches reveal the importance of quasi-bistable behavior in cellular decision processes on the example of the MAPK signaling pathway in PC-12 cell lines

BACKGROUND

Positive and negative feedback loops are ubiquitous motifs in biochemical signaling pathways. The mitogen-activated protein kinase (MAPK) pathway module is part of many distinct signaling networks and comprises several of these motifs, whose functioning depends on the cell line at hand and on the particular context. The maintainance of specificity of the response of the MAPK module to distinct stimuli has become a key paradigm especially in PC-12 cells, where the same module leads to different cell fates, depending on the stimulating growth factor. This cell fate is regulated by differences in the ERK (MAPK) activation profile, which shows a transient response upon stimulation with EGF, while the response is sustained in case of NGF. This behavior was explained by different effective network topologies. It is widely believed that this sustained response requires a bistable system.

RESULTS

In this study we present a sampling-based Bayesian model analysis on a dataset, in which PC-12 cells have been stimulated with different growth factors. This is combined with novel analysis methods to investigate the role of feedback interconnections to shape ERK response. Results strongly suggest that, besides bistability, an additional effect called quasi-bistability can contribute to explain the observed responses of the system to different stimuli. Quasi-bistability is the ability of a monostable system to maintain two distinct states over a long time period upon a transient signal, which is also related to positive feedback, but cannot be detected by standard steady state analysis methods.

CONCLUSIONS

Although applied on a specific example, our framework is generic enough to be also relevant for other regulatory network modeling studies that comprise positive feedback to explain cellular decision making processes. Overall, this study advices to focus not only on steady states, but also to take transient behavior into account in the analysis.

ERK signaling pathway regulates sleep duration through activity-induced gene expression during wakefulness

Wakefulness is accompanied by experience-dependent synaptic plasticity and an increase in activity-regulated gene transcription. Wake-induced genes are certainly markers of neuronal activity and may also directly regulate the duration of and need for sleep. We stimulated murine cortical cultures with the neuromodulatory signals that are known to control wakefulness in the brain and found that norepinephrine alone or a mixture of these neuromodulators induced activity-regulated gene transcription. Pharmacological inhibition of the various signaling pathways involved in the regulation of gene expression indicated that the extracellular signal-regulated kinase (ERK) pathway is the principal one mediating the effects of waking neuromodulators on gene expression. In mice, ERK phosphorylation in the cortex increased and decreased with wakefulness and sleep. Whole-body or cortical neuron-specific deletion of Erk1 or Erk2 significantly increased the duration of wakefulness in mice, and pharmacological inhibition of ERK phosphorylation decreased sleep duration and increased the duration of wakefulness bouts. Thus, this signaling pathway, which is highly conserved from Drosophila to mammals, is a key pathway that links waking experience-induced neuronal gene expression to sleep duration and quality.

Kinetic Modeling and Analysis of the Akt/Mechanistic Target of Rapamycin Complex 1 (mTORC1) Signaling Axis Reveals Cooperative, Feedforward Regulation

Mechanistic target of rapamycin complex 1 (mTORC1) controls biosynthesis and has been implicated in uncontrolled cell growth in cancer. Although many details of mTORC1 regulation are well understood, a systems-level, predictive framework synthesizing those details is currently lacking. We constructed various mathematical models of mTORC1 activation mediated by Akt and aligned the model outputs to kinetic data acquired for growth factor-stimulated cells. A model based on a putative feedforward loop orchestrated by Akt consistently predicted how the pathway was altered by depletion of key regulatory proteins. Analysis of the successful model also elucidates two dynamical motifs: neutralization of a negative regulator, which characterizes how Akt indirectly activates mTORC1, and seesaw enzyme regulation, which describes how activated and inhibited states of mTORC1 are controlled in concert to produce a nonlinear, ultrasensitive response. Such insights lend quantitative understanding of signaling networks and their precise manipulation in various contexts.

Single-Cell Imaging of ERK Signaling Using Fluorescent Biosensors

Single-cell analysis of the mitogen-activated protein kinase (MAPK) extracellular signal-regulated kinase (ERK) provides a means to perform highly detailed kinetic studies, assess heterogeneity between cells, and distinguish the subcellular localization of ERK activity. We describe here the methods needed to perform such measurements in a cell type of the investigator's choosing. We discuss the selection of appropriate reporters and provide detailed methods for stably introducing reporters, collecting live-cell data, and automatically extracting quantitative information from individual cells.

Cell Density-Dependent Increase in Tyrosine-Monophosphorylated ERK2 in MDCK Cells Expressing Active Ras or Raf

The extracellular signal-regulated kinase (ERK) is one of the principal hub proteins that transmit growth signals from upstream oncogene products including Ras and BRaf to downstream effector proteins. However, there are both reports supporting and refuting the increase in ERK activity in cancer tissues expressing the active Ras and BRaf proteins. We considered that the cell density might account for this discrepancy. To examine this possibility, we prepared Madin-Darby canine kidney (MDCK) cells that expressed an active HRas, NRas, KRas, or BRaf and an ERK biosensor based on the principle of Förster resonance energy transfer (FRET). As we anticipated, expression of the active Ras or BRaf increased ERK activity at low cell densities. However, the ERK activity was markedly suppressed at high cell densities irrespective of the expression of the active Ras or BRaf. Western blotting analysis with Phos-tag gel revealed the decrease of tyrosine and threonine-diphosphorylated active ERK and the increase of tyrosine-monophosphorylated inactive ERK at high cell density. In addition, we found that calyculin A, an inhibitor for PPP-subfamily protein serine/threonine phosphatases, decreased the tyrosine-monophosphorylated ERK. Our study suggests that PPP-subfamily phosphatases may be responsible for cell density-dependent ERK dephosphorylation in cancer cells expressing active Ras or BRaf protein.

Two New FRET Imaging Measures: Linearly Proportional to and Highly Contrasting the Fraction of Active Molecules

We developed two new FRET imaging measures for intramolecular FRET biosensors, called linearly proportional (LP) and highly contrasting (HC) measures, which can be easily calculated by the fluorescence intensities of donor and acceptor as a ratio between their weighted sums. As an alternative to the conventional ratiometric measure, which non-linearly depends on the fraction of active molecule, we first developed the LP measure, which is linearly proportional to the fraction of active molecules. The LP measure inherently unmixes bleed-through signals and is robust against fluorescence noise. By extending the LP measure, we furthermore designed the HC measure, which provides highly contrasting images of the molecular activity, more than the ratiometric measure. In addition to their advantages, these measures are insensitive to the biosensor expression level, which is a fundamental property of the ratiometric measure. Using artificial data and FRET imaging data, we showed that the LP measure effectively represents the fraction of active molecules and that the HC measure improves visual interpretability by providing high contrast images of molecular activity. Therefore, the LP and HC measures allow us to gain more quantitative and qualitative insights from FRET imaging than the ratiometric measure.

Modular bond-graph modelling and analysis of biomolecular systems

Bond graphs can be used to build thermodynamically-compliant hierarchical models of biomolecular systems. As bond graphs have been widely used to model, analyse and synthesise engineering systems, this study suggests that they can play the same rôle in the modelling, analysis and synthesis of biomolecular systems. The particular structure of bond graphs arising from biomolecular systems is established and used to elucidate the relation between thermodynamically closed and open systems. Block diagram representations of the dynamics implied by these bond graphs are used to reveal implicit feedback structures and are linearised to allow the application of control-theoretical methods. Two concepts of modularity are examined: computational modularity where physical correctness is retained and behavioural modularity where module behaviour (such as ultrasensitivity) is retained. As well as providing computational modularity, bond graphs provide a natural formulation of behavioural modularity and reveal the sources of retroactivity. A bond graph approach to reducing retroactivity, and thus inter-module interaction, is shown to require a power supply such as that provided by the ATP ⇌ ADP + Pi reaction. The mitogen-activated protein kinase cascade (Raf-MEK-ERK pathway) is used as an illustrative example.

Signal Transduction at the Single-Cell Level: Approaches to Study the Dynamic Nature of Signaling Networks

Signal transduction, or how cells interpret and react to external events, is a fundamental aspect of cellular function. Traditional study of signal transduction pathways involves mapping cellular signaling pathways at the population level. However, population-averaged readouts do not adequately illuminate the complex dynamics and heterogeneous responses found at the single-cell level. Recent technological advances that observe cellular response, computationally model signaling pathways, and experimentally manipulate cells now enable studying signal transduction at the single-cell level. These studies will enable deeper insights into the dynamic nature of signaling networks.

An approach for optimally extending mathematical models of signaling networks using omics data

Mathematical modeling is a key process in Systems Biology and the use of computational tools such as Cytoscape for omics data processing, need to be integrated in the modeling activity. In this paper we propose a new methodology for modeling signaling networks by combining ordinary differential equation models and a gene recommender system, GeneMANIA. We started from existing models, that are stored in the BioModels database, and we generated a query to use as input for the GeneMANIA algorithm. The output of the recommender system was then led back to the kinetic reactions that were finally added to the starting model. We applied the proposed methodology to EGFR-IGF1R signal transduction network, which plays an important role in translational oncology and cancer therapy of non small cell lung cancer.

A Transport Model for Estimating the Time Course of ERK Activation in the C. elegans Germline

The Caenorhabditis elegans germline is a well-studied model system for investigating the control of cell fate by signaling pathways. Cell signals at the distal tip of the germline promote cell proliferation; just before the loop, signals couple cell maturation to organism-level nutrient status; at the proximal end of the germline, signals coordinate oocyte maturation and fertilization in the presence of sperm. The latter two events require dual phosphorylation and activation of ERK, the effector molecule of the Ras/MAPK cascade. In C. elegans, ERK is known as MPK-1. At this point, none of today's methods for real-time monitoring of dually phosphorylated MPK-1 are working in the germline. Consequently, quantitative understanding of the MPK-1-dependent processes during germline development is limited. Here, we make a step toward advancing this understanding using a model-based framework that reconstructs the time course of MPK-1 activation from a snapshot of a fixed germline. Our approach builds on a number of recent studies for estimating temporal dynamics from fixed organisms, but takes advantage of the anatomy of the germline to simplify the analysis. Our model predicts that the MPK-1 signal turns on ∼30 h into germ cell progression and peaks ∼7 h later.

Cell signaling regulation by protein phosphorylation: a multivariate, heterogeneous, and context-dependent process

Proper spatiotemporal regulation of protein phosphorylation in cells and tissues is required for normal development and homeostasis, but aberrant protein phosphorylation regulation leads to various diseases. The study of signaling regulation by protein phosphorylation is complicated in part by the sheer scope of the kinome and phosphoproteome, dependence of signaling protein functionality on cellular localization, and the complex multivariate relationships that exist between protein phosphorylation dynamics and the cellular phenotypes they control. Additional complexities arise from the ability of microenvironmental factors to influence phosphorylation-dependent signaling and from the tendency for some signaling processes to occur heterogeneously among cells. These considerations should be taken into account when measuring cell signaling regulation by protein phosphorylation.

Unraveling Kinase Activation Dynamics Using Kinase-Substrate Relationships from Temporal Large-Scale Phosphoproteomics Studies

In response to stimuli, biological processes are tightly controlled by dynamic cellular signaling mechanisms. Reversible protein phosphorylation occurs on rapid time-scales (milliseconds to seconds), making it an ideal carrier of these signals. Advances in mass spectrometry-based proteomics have led to the identification of many tens of thousands of phosphorylation sites, yet for the majority of these the kinase is unknown and the underlying network topology of signaling networks therefore remains obscured. Identifying kinase substrate relationships (KSRs) is therefore an important goal in cell signaling research. Existing consensus sequence motif based prediction algorithms do not consider the biological context of KSRs, and are therefore insensitive to many other mechanisms guiding kinase-substrate recognition in cellular contexts. Here, we use temporal information to identify biologically relevant KSRs from Large-scale In Vivo Experiments (KSR-LIVE) in a data-dependent and automated fashion. First, we used available phosphorylation databases to construct a repository of existing experimentally-predicted KSRs. For each kinase in this database, we used time-resolved phosphoproteomics data to examine how its substrates changed in phosphorylation over time. Although substrates for a particular kinase clustered together, they often exhibited a different temporal pattern to the phosphorylation of the kinase. Therefore, although phosphorylation regulates kinase activity, our findings imply that substrate phosphorylation likely serve as a better proxy for kinase activity than kinase phosphorylation. KSR-LIVE can thereby infer which kinases are regulated within a biological context. Moreover, KSR-LIVE can also be used to automatically generate positive training sets for the subsequent prediction of novel KSRs using machine learning approaches. We demonstrate that this approach can distinguish between Akt and Rps6kb1, two kinases that share the same linear consensus motif, and provide evidence suggesting IRS-1 S265 as a novel Akt site. KSR-LIVE is an open-access algorithm that allows users to dissect phosphorylation signaling within a specific biological context, with the potential to be included in the standard analysis workflow for studying temporal high-throughput signal transduction data.

Gatekeeper Tyrosine Phosphorylation of SYMRK Is Essential for Synchronizing the Epidermal and Cortical Responses in Root Nodule Symbiosis

Symbiosis receptor kinase (SYMRK) is indispensable for activation of root nodule symbiosis (RNS) at both epidermal and cortical levels and is functionally conserved in legumes. Previously, we reported SYMRK to be phosphorylated on "gatekeeper" Tyr both in vitro as well as in planta. Since gatekeeper phosphorylation was not necessary for activity, the significance remained elusive. Herein, we show that substituting gatekeeper with nonphosphorylatable residues like Phe or Ala significantly affected autophosphorylation on selected targets on activation segment/αEF and β3-αC loop of SYMRK. In addition, the same gatekeeper mutants failed to restore proper symbiotic features in a symrk null mutant where rhizobial invasion of the epidermis and nodule organogenesis was unaffected but rhizobia remain restricted to the epidermis in infection threads migrating parallel to the longitudinal axis of the root, resulting in extensive infection patches at the nodule apex. Thus, gatekeeper phosphorylation is critical for synchronizing epidermal/cortical responses in RNS.

Spatio-temporal dynamics of a cell signal pathway with negative feedbacks: the MAPK/ERK pathway

We studied the spatio-temporal dynamics of a cell signal cascade with negative feedback that quantitatively emulates the regulative process that occurs in the Mitogen Activated Protein Kinase/Extracellular Regulated Kinase (MAPK/ERK) pathway. The model consists of a set of six coupled reaction-diffusion equations that describes the dynamics of the six-module pathway. In the basic module the active form of the protein transmits the signal to the next pathway’s module. As suggested by experiments, the model considers that the fifth module's kinase down-regulates the first and third modules. The feedback parameter is defined as, μ(r)( j)= k(kin)5/k(kin)(j), (j = 1, 3). We analysed the pathway's dynamics for μ(r)( j) = 0.10, 1.0, and 10 in the kinetic regimes: i) saturation of both kinases and phosphatases, ii) saturation of the phosphatases and iii) saturation of the kinases. For a regulated pathway the Total Activated Protein Profiles (TAPPs) as a function of time develop a maximum during the transient stage in the three kinetic regimes. These maxima become higher and their positions shift to longer times downstream. This scenario also applies to the TAPP's regulatory kinase that sums up its inhibitory action to that of the phosphatases leading to a maximum. Nevertheless, when μ(r)(j)= 1.0 , the TAPPs develop two maxima, with the second maximum being almost imperceptible. These results are in qualitative agreement with experimental data obtained from NIH 3T3 mouse fibroblasts. In addition, analyses of the stationary states as a function of position indicate that in the kinetic regime i) which is of physiological interest, signal transduction occurs with a relatively large propagation length for the three values of the regulative parameter. However, for μ(r)(j)= 0.10 , the sixth module concentration profile is transmitted with approximately 45% of its full value. The results obtained for μ(r)(j) = 10 , indicate that the first five concentration profiles are small with a short propagation length; nonetheless, the last concentration profile, c6, attains more than 90% of its full value with a relatively large propagation length as an indication of signal transduction. Signal transduction also occurred favourably in the kinetic regimes ii) and iii), but the signal was longer-ranged in the regime ii).

Network Architecture Predisposes an Enzyme to Either Pharmacologic or Genetic Targeting

Chemical inhibition and genetic knockdown of enzymes are not equivalent in cells, but network-level mechanisms that cause discrepancies between knockdown and inhibitor perturbations are not understood. Here we report that enzymes regulated by negative feedback are robust to knockdown but susceptible to inhibition. Using the Raf-MEK-ERK kinase cascade as a model system, we find that ERK activation is resistant to genetic knockdown of MEK but susceptible to a comparable degree of chemical MEK inhibition. We demonstrate that negative feedback from ERK to Raf causes this knockdown-versus-inhibitor discrepancy in vivo. Exhaustive mathematical modeling of three-tiered enzyme cascades suggests that this result is general: negative autoregulation or feedback favors inhibitor potency, whereas positive autoregulation or feedback favors knockdown potency. Our findings provide a rationale for selecting pharmacologic versus genetic perturbations in vivo and point out the dangers of using knockdown approaches in search of drug targets.

Conversion of graded phosphorylation into switch-like nuclear translocation via autoregulatory mechanisms in ERK signalling

The phosphorylation cascade in the extracellular signal-regulated kinase (ERK) pathway is a versatile reaction network motif that can potentially act as a switch, oscillator or memory. Nevertheless, there is accumulating evidence that the phosphorylation response is mostly linear to extracellular signals in mammalian cells. Here we find that subsequent nuclear translocation gives rise to a switch-like increase in nuclear ERK concentration in response to signal input. The switch-like response disappears in the presence of ERK inhibitor, suggesting the existence of autoregulatory mechanisms for ERK nuclear translocation involved in conversion from a graded to a switch-like response. In vitro reconstruction of ERK nuclear translocation indicates that ERK-mediated phosphorylation of nucleoporins regulates ERK translocation. A mathematical model and knockdown experiments suggest a contribution of nucleoporins to regulation of the ERK nuclear translocation response. Taken together, this study provides evidence that nuclear translocation with autoregulatory mechanisms acts as a switch in ERK signalling.

While ERK signalling can produce switch-like cell behaviour, phosphorylation of ERK increases linearly with extracellular signals. Here, the authors solve this seeming contradiction by showing that nuclear translocation of ERK behaves in a switch-like manner and is controlled by ERK activity.

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.

Phosphatase specificity and pathway insulation in signaling networks

Phosphatases play an important role in cellular signaling networks by regulating the phosphorylation state of proteins. Phosphatases are classically considered to be promiscuous, acting on tens to hundreds of different substrates. We recently demonstrated that a shared phosphatase can couple the responses of two proteins to incoming signals, even if those two substrates are from otherwise isolated areas of the network. This finding raises a potential paradox: if phosphatases are indeed highly promiscuous, how do cells insulate themselves against unwanted crosstalk? Here, we use mathematical models to explore three possible insulation mechanisms. One approach involves evolving phosphatase KM values that are large enough to prevent saturation by the phosphatase's substrates. Although this is an effective method for generating isolation, the phosphatase becomes a highly inefficient enzyme, which prevents the system from achieving switch-like responses and can result in slow response kinetics. We also explore the idea that substrate degradation can serve as an effective phosphatase. Assuming that degradation is unsaturatable, this mechanism could insulate substrates from crosstalk, but it would also preclude ultrasensitive responses and would require very high substrate turnover to achieve rapid dephosphorylation kinetics. Finally, we show that adaptor subunits, such as those found on phosphatases like PP2A, can provide effective insulation against phosphatase crosstalk, but only if their binding to substrates is uncoupled from their binding to the catalytic core. Analysis of the interaction network of PP2A's adaptor domains reveals that although its adaptors may isolate subsets of targets from one another, there is still a strong potential for phosphatase crosstalk within those subsets. Understanding how phosphatase crosstalk and the insulation mechanisms described here impact the function and evolution of signaling networks represents a major challenge for experimental and computational systems biology.

Novel Reporter for Faithful Monitoring of ERK2 Dynamics in Living Cells and Model Organisms

Uncoupling of ERK1/2 phosphorylation from subcellular localization is essential towards the understanding of molecular mechanisms that control ERK1/2-mediated cell-fate decision. ERK1/2 non-catalytic functions and discoveries of new specific anchors responsible of the subcellular compartmentalization of ERK1/2 signaling pathway have been proposed as regulation mechanisms for which dynamic monitoring of ERK1/2 localization is necessary. However, studying the spatiotemporal features of ERK2, for instance, in different cellular processes in living cells and tissues requires a tool that can faithfully report on its subcellular distribution. We developed a novel molecular tool, ERK2-LOC, based on the T2A-mediated coexpression of strictly equimolar levels of eGFP-ERK2 and MEK1, to faithfully visualize ERK2 localization patterns. MEK1 and eGFP-ERK2 were expressed reliably and functionally both in vitro and in single living cells. We then assessed the subcellular distribution and mobility of ERK2-LOC using fluorescence microscopy in non-stimulated conditions and after activation/inhibition of the MAPK/ERK1/2 signaling pathway. Finally, we used our coexpression system in Xenopus laevis embryos during the early stages of development. This is the first report on MEK1/ERK2 T2A-mediated coexpression in living embryos, and we show that there is a strong correlation between the spatiotemporal subcellular distribution of ERK2-LOC and the phosphorylation patterns of ERK1/2. Our approach can be used to study the spatiotemporal localization of ERK2 and its dynamics in a variety of processes in living cells and embryonic tissues.

Digital signaling decouples activation probability and population heterogeneity

Digital signaling enhances robustness of cellular decisions in noisy environments, but it is unclear how digital systems transmit temporal information about a stimulus. To understand how temporal input information is encoded and decoded by the NF-κB system, we studied transcription factor dynamics and gene regulation under dose- and duration-modulated inflammatory inputs. Mathematical modeling predicted and microfluidic single-cell experiments confirmed that integral of the stimulus (or area, concentration × duration) controls the fraction of cells that activate NF-κB in the population. However, stimulus temporal profile determined NF-κB dynamics, cell-to-cell variability, and gene expression phenotype. A sustained, weak stimulation lead to heterogeneous activation and delayed timing that is transmitted to gene expression. In contrast, a transient, strong stimulus with the same area caused rapid and uniform dynamics. These results show that digital NF-κB signaling enables multidimensional control of cellular phenotype via input profile, allowing parallel and independent control of single-cell activation probability and population heterogeneity.

Conditional robustness analysis for fragility discovery and target identification in biochemical networks and in cancer systems biology

Background

The study of cancer therapy is a key issue in the field of oncology research and the development of target therapies is one of the main problems currently under investigation. This is particularly relevant in different types of tumor where traditional chemotherapy approaches often fail, such as lung cancer.

Results

We started from the general definition of robustness introduced by Kitano and applied it to the analysis of dynamical biochemical networks, proposing a new algorithm based on moment independent analysis of input/output uncertainty. The framework utilizes novel computational methods which enable evaluating the model fragility with respect to quantitative performance measures and parameters such as reaction rate constants and initial conditions. The algorithm generates a small subset of parameters that can be used to act on complex networks and to obtain the desired behaviors. We have applied the proposed framework to the EGFR-IGF1R signal transduction network, a crucial pathway in lung cancer, as an example of Cancer Systems Biology application in drug discovery. Furthermore, we have tested our framework on a pulse generator network as an example of Synthetic Biology application, thus proving the suitability of our methodology to the characterization of the input/output synthetic circuits.

Conclusions

The achieved results are of immediate practical application in computational biology, and while we demonstrate their use in two specific examples, they can in fact be used to study a wider class of biological systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0216-5) contains supplementary material, which is available to authorized users.

Receptor Level Mechanisms Are Required for Epidermal Growth Factor (EGF)-stimulated Extracellular Signal-regulated Kinase (ERK) Activity Pulses

In both physiological and cell culture systems, EGF-stimulated ERK activity occurs in discrete pulses within individual cells. Many feedback loops are present in the EGF receptor (EGFR)-ERK network, but the mechanisms driving pulsatile ERK kinetics are unknown. Here, we find that in cells that respond to EGF with frequency-modulated pulsatile ERK activity, stimulation through a heterologous TrkA receptor system results in non-pulsatile, amplitude-modulated activation of ERK. We further dissect the kinetics of pulse activity using a combination of FRET- and translocation-based reporters and find that EGFR activity is required to maintain ERK activity throughout the 10-20-minute lifetime of pulses. Together, these data indicate that feedbacks operating within the core Ras-Raf-MEK-ERK cascade are insufficient to drive discrete pulses of ERK activity and instead implicate mechanisms acting at the level of EGFR.

The mouse radial spoke protein 3 is a nucleocytoplasmic shuttling protein that promotes neurogenesis

Radial spoke protein 3 (RSP3) was first identified in Chlamydomonas as a component of radial spoke, which is important for flagellar motility. The mammalian homolog of the Chlamydomonas RSP3 protein is found to be a mammalian protein kinase A-anchoring protein that binds ERK1/2. Here we show that mouse RSP3 is a nucleocytoplasmic shuttling protein. The full-length RSP3-EGFP fusion protein is mainly located in the cytoplasm of Chinese hamster ovary cells. However, by using deletion mutants of RSP3, we identified two nuclear localization signals and a nuclear export signal in RSP3. Moreover, using in utero electroporation, we found that overexpression of RSP3 in the developing cerebral cortex promotes neurogenesis. The layer II/III of the neocortex was much thicker in the RSP3-transfected region than that of the untransfected region in the neocortex. We also show that RSP3 is specifically located in the primary cilia of the radial glial cells, where it acts as a signaling mediator that regulates neurogenesis. Thus, our results suggest that RSP3 is a nucleocytoplasmic shuttling protein and plays an essential role in neurogenesis.

Using fluctuation analysis to establish causal relations between cellular events without experimental perturbation

Experimental perturbations are commonly used to establish causal relationships between the molecular components of a pathway and their cellular functions; however, this approach suffers inherent limitations. Especially in pathways with a significant level of nonlinearity and redundancy among components, such perturbations induce compensatory responses that obscure the actual function of the targeted component in the unperturbed pathway. A complementary approach uses constitutive fluctuations in component activities to identify the hierarchy of information flow through pathways. Here, we review the motivation for using perturbation-free approaches and highlight recent advances made in using perturbation-free fluctuation analysis as a means to establish causality among cellular events.

A single mechanism can explain network-wide insulin resistance in adipocytes from obese patients with type 2 diabetes

The response to insulin is impaired in type 2 diabetes. Much information is available about insulin signaling, but understanding of the cellular mechanisms causing impaired signaling and insulin resistance is hampered by fragmented data, mainly obtained from different cell lines and animals. We have collected quantitative and systems-wide dynamic data on insulin signaling in primary adipocytes and compared cells isolated from healthy and diabetic individuals. Mathematical modeling and experimental verification identified mechanisms of insulin control of the MAPKs ERK1/2. We found that in human adipocytes, insulin stimulates phosphorylation of the ribosomal protein S6 and hence protein synthesis about equally via ERK1/2 and mTORC1. Using mathematical modeling, we examined the signaling network as a whole and show that a single mechanism can explain the insulin resistance of type 2 diabetes throughout the network, involving signaling both through IRS1, PKB, and mTOR and via ERK1/2 to the nuclear transcription factor Elk1. The most important part of the insulin resistance mechanism is an attenuated feedback from the protein kinase mTORC1 to IRS1, which spreads signal attenuation to all parts of the insulin signaling network. Experimental inhibition of mTORC1 using rapamycin in adipocytes from non-diabetic individuals induced and thus confirmed the predicted network-wide insulin resistance.

Data-driven modeling reconciles kinetics of ERK phosphorylation, localization, and activity states

The extracellular signal‐regulated kinase (ERK) signaling pathway controls cell proliferation and differentiation in metazoans. Two hallmarks of its dynamics are adaptation of ERK phosphorylation, which has been linked to negative feedback, and nucleocytoplasmic shuttling, which allows active ERK to phosphorylate protein substrates in the nucleus and cytosol. To integrate these complex features, we acquired quantitative biochemical and live‐cell microscopy data to reconcile phosphorylation, localization, and activity states of ERK. While maximal growth factor stimulation elicits transient ERK phosphorylation and nuclear translocation responses, ERK activities available to phosphorylate substrates in the cytosol and nuclei show relatively little or no adaptation. Free ERK activity in the nucleus temporally lags the peak in nuclear translocation, indicating a slow process. Additional experiments, guided by kinetic modeling, show that this process is consistent with ERK's modification of and release from nuclear substrate anchors. Thus, adaptation of whole‐cell ERK phosphorylation is a by‐product of transient protection from phosphatases. Consistent with this interpretation, predictions concerning the dose‐dependence of the pathway response and its interruption by inhibition of MEK were experimentally confirmed.