Feed-forward stimulation of CAMK2 by the oncogenic pseudokinase PEAK1 generates a therapeutically 'actionable' signalling axis in triple negative breast cancer.

The PEAK family of pseudokinases, comprising PEAK1-3, are signalling scaffolds that play oncogenic roles in several poor prognosis human cancers, including triple negative breast cancer (TNBC). However, therapeutic targeting of pseudokinases is challenging due to their lack of catalytic activity. To address this, we screened for PEAK1 effectors by affinity purification and mass spectrometry, identifying calcium/calmodulin-dependent protein kinase 2 (CAMK2)D and CAMK2G. PEAK1 promoted CAMK2D/G activation in TNBC cells via a novel feed-forward mechanism involving PEAK1/PLCγ1/Ca 2+ signalling and direct binding via a consensus CAMK2 interaction motif in the PEAK1 N-terminus. In turn, CAMK2 phosphorylated PEAK1 to enhance association with PEAK2, which is critical for PEAK1 oncogenic signalling. To achieve pharmacologic targeting of PEAK1/CAMK2, we repurposed RA306, a second generation CAMK2 inhibitor under pre-clinical development for treatment of cardiovascular disease. RA306 demonstrated on-target activity against CAMK2 in TNBC cells and inhibited PEAK1-enhanced migration and invasion in vitro . Moreover, RA306 significantly attenuated TNBC xenograft growth and blocked metastasis in a manner mirrored by CRISPR-mediated PEAK1 ablation. Overall, these studies establish PEAK1 as a critical cell signalling nexus, identify a novel mechanism for regulation of Ca 2+ signalling and its integration with tyrosine kinase signals, and identify CAMK2 as a therapeutically 'actionable' target downstream of PEAK1.

Integrative modeling and analysis of signaling crosstalk reveal molecular switches coordinating Yes-associated protein transcriptional activities.

The transcriptional co-activator YAP forms complexes with distinct transcription factors, controlling cell fate decisions, such as proliferation and apoptosis. However, the mechanisms underlying its context-dependent function are poorly defined. This study explores the interplay between the TGF-ß and Hippo pathways and their influence on YAP's association with specific transcription factors. By integrating iterative mathematical modeling with experimental validation, we uncover molecular switches, predominantly controlled by RASSF1A and ITCH, which dictate the formation of YAP-SMAD (proliferative) and YAP-p73 (apoptotic) complexes. Our results show that RASSF1A enhances the formation of apoptotic complexes, whereas ITCH promotes the formation of proliferative complexes. Notably, higher levels of ITCH transform YAP-SMAD activity from a transient to a sustained state, impacting cellular behaviors. Extending these findings to various breast cancer cell lines highlights the role of cellular context in YAP regulation. Our study provides new insights into the mechanisms of YAP transcriptional activities and their therapeutic implications.

Integrative modeling uncovers p21-driven drug resistance and prioritizes therapies for PIK3CA-mutant breast cancer.

Utility of PI3Kα inhibitors like BYL719 is limited by the acquisition of genetic and non-genetic mechanisms of resistance which cause disease recurrence. Several combination therapies based on PI3K inhibition have been proposed, but a way to systematically prioritize them for breast cancer treatment is still missing. By integrating published and in-house studies, we have developed in silico models that quantitatively capture dynamics of PI3K signaling at the network-level under a BYL719-sensitive versus BYL719 resistant-cell state. Computational predictions show that signal rewiring to alternative components of the PI3K pathway promote resistance to BYL719 and identify PDK1 as the most effective co-target with PI3Kα rescuing sensitivity of resistant cells to BYL719. To explore whether PI3K pathway-independent mechanisms further contribute to BYL719 resistance, we performed phosphoproteomics and found that selection of high levels of the cell cycle regulator p21 unexpectedly promoted drug resistance in T47D cells. Functionally, high p21 levels favored repair of BYL719-induced DNA damage and bypass of the associated cellular senescence. Importantly, targeted inhibition of the check-point inhibitor CHK1 with MK-8776 effectively caused death of p21-high T47D cells, thus establishing a new vulnerability of BYL719-resistant breast cancer cells. Together, our integrated studies uncover hidden molecular mediators causing resistance to PI3Kα inhibition and provide a framework to prioritize combination therapies for PI3K-mutant breast cancer.

SynDISCO: A Mechanistic Modeling-Based Framework for Predictive Prioritization of Synergistic Drug Combinations Targeting Cell Signalling Networks.

The widespread development of resistance to cancer monotherapies has prompted the need to identify combinatorial treatment approaches that circumvent drug resistance and achieve more durable clinical benefit. However, given the vast space of possible combinations of existing drugs, the inaccessibility of drug screens to candidate targets with no available drugs, and the significant heterogeneity of cancers, exhaustive experimental testing of combination treatments remains highly impractical. There is thus an urgent need to develop computational approaches that complement experimental efforts and aid the identification and prioritization of effective drug combinations. Here, we provide a practical guide to SynDISCO, a computational framework that leverages mechanistic ODE modeling to predict and prioritize synergistic combination treatments directed at signaling networks. We demonstrate the key steps of SynDISCO and its application to the EGFR-MET signaling network in triple negative breast cancer as an illustrative example. SynDISCO is, however, a network- and cancer-independent framework, and given a suitable ODE model of the network of interest, it could be leveraged to discover cancer-specific combination treatments.

A Boolean-based machine learning framework identifies predictive biomarkers of HSP90-targeted therapy response in prostate cancer.

Precision medicine has emerged as an important paradigm in oncology, driven by the significant heterogeneity of individual patients' tumour. A key prerequisite for effective implementation of precision oncology is the development of companion biomarkers that can predict response to anti-cancer therapies and guide patient selection for clinical trials and/or treatment. However, reliable predictive biomarkers are currently lacking for many anti-cancer therapies, hampering their clinical application. Here, we developed a novel machine learning-based framework to derive predictive multi-gene biomarker panels and associated expression signatures that accurately predict cancer drug sensitivity. We demonstrated the power of the approach by applying it to identify response biomarker panels for an Hsp90-based therapy in prostate cancer, using proteomic data profiled from prostate cancer patient-derived explants. Our approach employs a rational feature section strategy to maximise model performance, and innovatively utilizes Boolean algebra methods to derive specific expression signatures of the marker proteins. Given suitable data for model training, the approach is also applicable to other cancer drug agents in different tumour settings.

The pseudokinase NRBP1 activates Rac1/Cdc42 via P-Rex1 to drive oncogenic signalling in triple-negative breast cancer.

We have determined that expression of the pseudokinase NRBP1 positively associates with poor prognosis in triple negative breast cancer (TNBC) and is required for efficient migration, invasion and proliferation of TNBC cells in culture as well as growth of TNBC orthotopic xenografts and experimental metastasis. Application of BioID/MS profiling identified P-Rex1, a known guanine nucleotide exchange factor for Rac1, as a NRBP1 binding partner. Importantly, NRBP1 overexpression enhanced levels of GTP-bound Rac1 and Cdc42 in a P-Rex1-dependent manner, while NRBP1 knockdown reduced their activation. In addition, NRBP1 associated with P-Rex1, Rac1 and Cdc42, suggesting a scaffolding function for this pseudokinase. NRBP1-mediated promotion of cell migration and invasion was P-Rex1-dependent, while constitutively-active Rac1 rescued the effect of NRBP1 knockdown on cell proliferation and invasion. Generation of reactive oxygen species via a NRBP1/P-Rex1 pathway was implicated in these oncogenic roles of NRBP1. Overall, these findings define a new function for NRBP1 and a novel oncogenic signalling pathway in TNBC that may be amenable to therapeutic intervention.

Dynamic modelling of the PI3K/MTOR signalling network uncovers biphasic dependence of mTORC1 activity on the mTORC2 subunit SIN1.

The PI3K/MTOR signalling network regulates a broad array of critical cellular processes, including cell growth, metabolism and autophagy. The mechanistic target of rapamycin (MTOR) kinase functions as a core catalytic subunit in two physically and functionally distinct complexes mTORC1 and mTORC2, which also share other common components including MLST8 (also known as GßL) and DEPTOR. Despite intensive research, how mTORC1 and 2 assembly and activity are coordinated, and how they are functionally linked remain to be fully characterized. This is due in part to the complex network wiring, featuring multiple feedback loops and intricate post-translational modifications. Here, we integrate predictive network modelling, in vitro experiments and -omics data analysis to elucidate the emergent dynamic behaviour of the PI3K/MTOR network. We construct new mechanistic models that encapsulate critical mechanistic details, including mTORC1/2 coordination by MLST8 (de)ubiquitination and the Akt-to-mTORC2 positive feedback loop. Model simulations validated by experimental studies revealed a previously unknown biphasic, threshold-gated dependence of mTORC1 activity on the key mTORC2 subunit SIN1, which is robust against cell-to-cell variation in protein expression. In addition, our integrative analysis demonstrates that ubiquitination of MLST8, which is reversed by OTUD7B, is regulated by IRS1/2. Our results further support the essential role of MLST8 in enabling both mTORC1 and 2's activity and suggest MLST8 as a viable therapeutic target in breast cancer. Overall, our study reports a new mechanistic model of PI3K/MTOR signalling incorporating MLST8-mediated mTORC1/2 formation and unveils a novel regulatory linkage between mTORC1 and mTORC2.

Evaluation of FGFR targeting in breast cancer through interrogation of patient-derived models.

BACKGROUND: Particular breast cancer subtypes pose a clinical challenge due to limited targeted therapeutic options and/or poor responses to the existing targeted therapies. While cell lines provide useful pre-clinical models, patient-derived xenografts (PDX) and organoids (PDO) provide significant advantages, including maintenance of genetic and phenotypic heterogeneity, 3D architecture and for PDX, tumor-stroma interactions. In this study, we applied an integrated multi-omic approach across panels of breast cancer PDXs and PDOs in order to identify candidate therapeutic targets, with a major focus on specific FGFRs. METHODS: MS-based phosphoproteomics, RNAseq, WES and Western blotting were used to characterize aberrantly activated protein kinases and effects of specific FGFR inhibitors. PDX and PDO were treated with the selective tyrosine kinase inhibitors AZD4547 (FGFR1-3) and BLU9931 (FGFR4). FGFR4 expression in cancer tissue samples and PDOs was assessed by immunohistochemistry. METABRIC and TCGA datasets were interrogated to identify specific FGFR alterations and their association with breast cancer subtype and patient survival. RESULTS: Phosphoproteomic profiling across 18 triple-negative breast cancers (TNBC) and 1 luminal B PDX revealed considerable heterogeneity in kinase activation, but 1/3 of PDX exhibited enhanced phosphorylation of FGFR1, FGFR2 or FGFR4. One TNBC PDX with high FGFR2 activation was exquisitely sensitive to AZD4547. Integrated 'omic analysis revealed a novel FGFR2-SKI fusion that comprised the majority of FGFR2 joined to the C-terminal region of SKI containing the coiled-coil domains. High FGFR4 phosphorylation characterized a luminal B PDX model and treatment with BLU9931 significantly decreased tumor growth. Phosphoproteomic and transcriptomic analyses confirmed on-target action of the two anti-FGFR drugs and also revealed novel effects on the spliceosome, metabolism and extracellular matrix (AZD4547) and RIG-I-like and NOD-like receptor signaling (BLU9931). Interrogation of public datasets revealed FGFR2 amplification, fusion or mutation in TNBC and other breast cancer subtypes, while FGFR4 overexpression and amplification occurred in all breast cancer subtypes and were associated with poor prognosis. Characterization of a PDO panel identified a luminal A PDO with high FGFR4 expression that was sensitive to BLU9931 treatment, further highlighting FGFR4 as a potential therapeutic target. CONCLUSIONS: This work highlights how patient-derived models of human breast cancer provide powerful platforms for therapeutic target identification and analysis of drug action, and also the potential of specific FGFRs, including FGFR4, as targets for precision treatment.

Akt phosphorylates insulin receptor substrate to limit PI3K-mediated PIP3 synthesis.

The phosphoinositide 3-kinase (PI3K)-Akt network is tightly controlled by feedback mechanisms that regulate signal flow and ensure signal fidelity. A rapid overshoot in insulin-stimulated recruitment of Akt to the plasma membrane has previously been reported, which is indicative of negative feedback operating on acute timescales. Here, we show that Akt itself engages this negative feedback by phosphorylating insulin receptor substrate (IRS) 1 and 2 on a number of residues. Phosphorylation results in the depletion of plasma membrane-localised IRS1/2, reducing the pool available for interaction with the insulin receptor. Together these events limit plasma membrane-associated PI3K and phosphatidylinositol (3,4,5)-trisphosphate (PIP3) synthesis. We identified two Akt-dependent phosphorylation sites in IRS2 at S306 (S303 in mouse) and S577 (S573 in mouse) that are key drivers of this negative feedback. These findings establish a novel mechanism by which the kinase Akt acutely controls PIP3 abundance, through post-translational modification of the IRS scaffold.

For the body to work properly, cells must constantly 'talk' to each other using signalling molecules. Receiving a chemical signal triggers a series of molecular events in a cell, a so-called 'signal transduction pathway' that connects a signal with a precise outcome. Disturbing cell signalling can trigger disease, and strict control mechanisms are therefore in place to ensure that communication does not break down or become erratic. For instance, just as a thermostat turns off the heater once the right temperature is reached, negative feedback mechanisms in cells switch off signal transduction pathways when the desired outcome has been achieved. The hormone insulin is a signal for growth that increases in the body following a meal to promote the storage of excess blood glucose (sugar) in muscle and fat cells. The hormone binds to insulin receptors at the cell surface and switches on a signal transduction pathway that makes the cell take up glucose from the bloodstream. If the signal is not engaged diseases such as diabetes develop. Conversely, if the signal cannot be adequately switched of cancer can develop. Determining exactly how insulin works would help to understand these diseases better and to develop new treatments. Kearney et al. therefore set out to examine the biochemical 'fail-safes' that control insulin signalling. Experiments using computer simulations of the insulin signalling pathway revealed a potential new mechanism for negative feedback, which centred on a molecule known as Akt. The models predicted that if the negative feedback were removed, then Akt would become hyperactive and accumulate at the cell's surface after stimulation with insulin. Further manipulation of the 'virtual' insulin signalling pathway and studies of live cells in culture confirmed that this was indeed the case. The cell biology experiments also showed how Akt, once at the cell surface, was able to engage the negative feedback and shut down further insulin signalling. Akt did this by inactivating a protein required to pass the signal from the insulin receptor to the rest of the cell. Overall, this work helps to understand cell communication by revealing a previously unknown, and critical component of the insulin signalling pathway.

Feedback, Crosstalk and Competition: Ingredients for Emergent Non-Linear Behaviour in the PI3K/mTOR Signalling Network.

The PI3K/mTOR signalling pathway plays a central role in the governing of cell growth, survival and metabolism. As such, it must integrate and decode information from both external and internal sources to guide efficient decision-making by the cell. To facilitate this, the pathway has evolved an intricate web of complex regulatory mechanisms and elaborate crosstalk with neighbouring signalling pathways, making it a highly non-linear system. Here, we describe the mechanistic biological details that underpin these regulatory mechanisms, covering a multitude of negative and positive feedback loops, feed-forward loops, competing protein interactions, and crosstalk with major signalling pathways. Further, we highlight the non-linear and dynamic network behaviours that arise from these regulations, uncovered through computational and experimental studies. Given the pivotal role of the PI3K/mTOR network in cellular homeostasis and its frequent dysregulation in pathologies including cancer and diabetes, a coherent and systems-level understanding of the complex regulation and consequential dynamic signalling behaviours within this network is imperative for advancing biology and development of new therapeutic approaches.

Signaling Heterogeneity is Defined by Pathway Architecture and Intercellular Variability in Protein Expression.

Insulin's activation of PI3K/Akt signaling, stimulates glucose uptake by enhancing delivery of GLUT4 to the cell surface. Here we examined the origins of intercellular heterogeneity in insulin signaling. Akt activation alone accounted for ~25% of the variance in GLUT4, indicating that additional sources of variance exist. The Akt and GLUT4 responses were highly reproducible within the same cell, suggesting the variance is between cells (extrinsic) and not within cells (intrinsic). Generalized mechanistic models (supported by experimental observations) demonstrated that the correlation between the steady-state levels of two measured signaling processes decreases with increasing distance from each other and that intercellular variation in protein expression (as an example of extrinsic variance) is sufficient to account for the variance in and between Akt and GLUT4. Thus, the response of a population to insulin signaling is underpinned by considerable single-cell heterogeneity that is largely driven by variance in gene/protein expression between cells.

Control of Glucocorticoid Receptor Levels by PTEN Establishes a Failsafe Mechanism for Tumor Suppression.

The PTEN tumor suppressor controls cell death and survival by regulating functions of various molecular targets. While the role of PTEN lipid-phosphatase activity on PtdIns(3,4,5)P3 and inhibition of PI3K pathway is well characterized, the biological relevance of PTEN protein-phosphatase activity remains undefined. Here, using knockin (KI) mice harboring cancer-associated and functionally relevant missense mutations, we show that although loss of PTEN lipid-phosphatase function cooperates with oncogenic PI3K to promote rapid mammary tumorigenesis, the additional loss of PTEN protein-phosphatase activity triggered an extensive cell death response evident in early and advanced mammary tumors. Omics and drug-targeting studies revealed that PI3Ks act to reduce glucocorticoid receptor (GR) levels, which are rescued by loss of PTEN protein-phosphatase activity to restrain cell survival. Thus, we find that the dual regulation of GR by PI3K and PTEN functions as a rheostat that can be exploited for the treatment of PTEN loss-driven cancers.

Global redox proteome and phosphoproteome analysis reveals redox switch in Akt.

Protein oxidation sits at the intersection of multiple signalling pathways, yet the magnitude and extent of crosstalk between oxidation and other post-translational modifications remains unclear. Here, we delineate global changes in adipocyte signalling networks following acute oxidative stress and reveal considerable crosstalk between cysteine oxidation and phosphorylation-based signalling. Oxidation of key regulatory kinases, including Akt, mTOR and AMPK influences the fidelity rather than their absolute activation state, highlighting an unappreciated interplay between these modifications. Mechanistic analysis of the redox regulation of Akt identified two cysteine residues in the pleckstrin homology domain (C60 and C77) to be reversibly oxidized. Oxidation at these sites affected Akt recruitment to the plasma membrane by stabilizing the PIP3 binding pocket. Our data provide insights into the interplay between oxidative stress-derived redox signalling and protein phosphorylation networks and serve as a resource for understanding the contribution of cellular oxidation to a range of diseases.

Coupled feedback regulation of nuclear factor of activated T-cells (NFAT) modulates activation-induced cell death of T cells.

A properly functioning immune system is vital for an organism's wellbeing. Immune tolerance is a critical feature of the immune system that allows immune cells to mount effective responses against exogenous pathogens such as viruses and bacteria, while preventing attack to self-tissues. Activation-induced cell death (AICD) in T lymphocytes, in which repeated stimulations of the T-cell receptor (TCR) lead to activation and then apoptosis of T cells, is a major mechanism for T cell homeostasis and helps maintain peripheral immune tolerance. Defects in AICD can lead to development of autoimmune diseases. Despite its importance, the regulatory mechanisms that underlie AICD remain poorly understood, particularly at an integrative network level. Here, we develop a dynamic multi-pathway model of the integrated TCR signalling network and perform model-based analysis to characterize the network-level properties of AICD. Model simulation and analysis show that amplified activation of the transcriptional factor NFAT in response to repeated TCR stimulations, a phenomenon central to AICD, is tightly modulated by a coupled positive-negative feedback mechanism. NFAT amplification is predominantly enabled by a positive feedback self-regulated by NFAT, while opposed by a NFAT-induced negative feedback via Carabin. Furthermore, model analysis predicts an optimal therapeutic window for drugs that help minimize proliferation while maximize AICD of T cells. Overall, our study provides a comprehensive mathematical model of TCR signalling and model-based analysis offers new network-level insights into the regulation of activation-induced cell death in T cells.

Systems modelling of the EGFR-PYK2-c-Met interaction network predicts and prioritizes synergistic drug combinations for triple-negative breast cancer.

Prediction of drug combinations that effectively target cancer cells is a critical challenge for cancer therapy, in particular for triple-negative breast cancer (TNBC), a highly aggressive breast cancer subtype with no effective targeted treatment. As signalling pathway networks critically control cancer cell behaviour, analysis of signalling network activity and crosstalk can help predict potent drug combinations and rational stratification of patients, thus bringing therapeutic and prognostic values. We have previously showed that the non-receptor tyrosine kinase PYK2 is a downstream effector of EGFR and c-Met and demonstrated their crosstalk signalling in basal-like TNBC. Here we applied a systems modelling approach and developed a mechanistic model of the integrated EGFR-PYK2-c-Met signalling network to identify and prioritize potent drug combinations for TNBC. Model predictions validated by experimental data revealed that among six potential combinations of drug pairs targeting the central nodes of the network, including EGFR, c-Met, PYK2 and STAT3, co-targeting of EGFR and PYK2 and to a lesser extent of EGFR and c-Met yielded strongest synergistic effect. Importantly, the synergy in co-targeting EGFR and PYK2 was linked to switch-like cell proliferation-associated responses. Moreover, simulations of patient-specific models using public gene expression data of TNBC patients led to predictive stratification of patients into subgroups displaying distinct susceptibility to specific drug combinations. These results suggest that mechanistic systems modelling is a powerful approach for the rational design, prediction and prioritization of potent combination therapies for individual patients, thus providing a concrete step towards personalized treatment for TNBC and other tumour types.

Dissecting Cell-Fate Determination Through Integrated Mathematical Modeling of the ERK/MAPK Signaling Pathway.

The past three decades have witnessed an enormous progress in the elucidation of the ERK/MAPK signaling pathway and its involvement in various cellular processes. Because of its importance and complex wiring, the ERK pathway has been an intensive subject for mathematical modeling, which facilitates the unraveling of key dynamic properties and behaviors of the pathway. Recently, however, it became evident that the pathway does not act in isolation but closely interacts with many other pathways to coordinate various cellular outcomes under different pathophysiological contexts. This has led to an increasing number of integrated, large-scale models that link the ERK pathway to other functionally important pathways. In this chapter, we first discuss the essential steps in model development and notable models of the ERK pathway. We then use three examples of integrated, multipathway models to investigate how crosstalk of ERK signaling with other pathways regulates cell-fate decision-making in various physiological and disease contexts. Specifically, we focus on ERK interactions with the phosphoinositide-3 kinase (PI3K), c-Jun N-terminal kinase (JNK), and ß-adrenergic receptor (ß-AR) signaling pathways. We conclude that integrated modeling in combination with wet-lab experimentation have been and will be instrumental in gaining an in-depth understanding of ERK signaling in multiple biological contexts.

Targeting of PYK2 Synergizes with EGFR Antagonists in Basal-like TNBC and Circumvents HER3-Associated Resistance via the NEDD4-NDRG1 Axis.

Triple-negative breast cancer (TNBC) is a highly aggressive, heterogeneous disease with poor prognosis and no effective targeted therapies. EGFR is highly expressed in basal-like TNBC and is considered as a potential therapeutic target. However, EGFR targeting exerts only marginal clinical benefits, possibly due to activation of compensatory signaling pathways, which are frequently associated with HER3 upregulation. Here we show that concomitant targeting of EGFR and the nonreceptor tyrosine kinases PYK2/FAK synergistically inhibits the proliferation of basal-like TNBC cells in vitro and attenuates tumor growth in a mouse xenograft model. Dual targeting of EGFR and PYK2/FAK inhibited complementary key growth and survival pathways mediated by AKT, S6K, STAT3, and ERK1/2 activation. PYK2 inhibition also abrogated HER3 upregulation in response to EGFR antagonists, thereby circumventing HER3-associated drug resistance. Mechanistically, PYK2 inhibition facilitated the proteasomal degradation of HER3 while inducing upregulation of NDRG1 (N-myc downstream regulated 1 gene). NDRG1 enhanced the interaction of HER3 with the ubiquitin ligase NEDD4, while PYK2, which interacts with NEDD4 and HER3, interfered with NEDD4-HER3 binding, suggesting that the PYK2-NDRG1-NEDD4 circuit has a critical role in receptor degradation, drug response, and resistance mechanism. Our studies offer a preclinical proof of concept for a strategy of cotargeting the EGFR and PYK2/FAK kinases to improve TNBC therapy. Cancer Res; 77(1); 86-99. ©2016 AACR.

A Regulated Double-Negative Feedback Decodes the Temporal Gradient of Input Stimulation in a Cell Signaling Network.

Revealing the hidden mechanism of how cells sense and react to environmental signals has been a central question in cell biology. We focused on the rate of increase of stimulation, or temporal gradient, known to cause different responses of cells. We have investigated all possible three-node enzymatic networks and identified a network motif that robustly generates a transient or sustained response by acute or gradual stimulation, respectively. We also found that a regulated double-negative feedback within the motif is essential for the temporal gradient-sensitive switching. Our analysis highlights the essential structure and mechanism enabling cells to properly respond to dynamic environmental changes.

Unveiling Hidden Dynamics of Hippo Signalling: A Systems Analysis.

The Hippo signalling pathway has recently emerged as an important regulator of cell apoptosis and proliferation with significant implications in human diseases. In mammals, the pathway contains the core kinases MST1/2, which phosphorylate and activate LATS1/2 kinases. The pro-apoptotic function of the MST/LATS signalling axis was previously linked to the Akt and ERK MAPK pathways, demonstrating that the Hippo pathway does not act alone but crosstalks with other signalling pathways to coordinate network dynamics and cellular outcomes. These crosstalks were characterised by a multitude of complex regulatory mechanisms involving competitive protein-protein interactions and phosphorylation mediated feedback loops. However, how these different mechanisms interplay in different cellular contexts to drive the context-specific network dynamics of Hippo-ERK signalling remains elusive. Using mathematical modelling and computational analysis, we uncovered that the Hippo-ERK network can generate highly diverse dynamical profiles that can be clustered into distinct dose-response patterns. For each pattern, we offered mechanistic explanation that defines when and how the observed phenomenon can arise. We demonstrated that Akt displays opposing, dose-dependent functions towards ERK, which are mediated by the balance between the Raf-1/MST2 protein interaction module and the LATS1 mediated feedback regulation. Moreover, Ras displays a multi-functional role and drives biphasic responses of both MST2 and ERK activities; which are critically governed by the competitive protein interaction between MST2 and Raf-1. Our study represents the first in-depth and systematic analysis of the Hippo-ERK network dynamics and provides a concrete foundation for future studies.