Myc Dysregulation in Activated Macrophages Initiates Iron-Mediated Lipid Peroxidation that Fuels Type I Interferon and Compromises TB Resistance

A quarter of human population is infected with Mycobacterium tuberculosis, but less than 10% of those infected develop clinical, mostly pulmonary, TB. To dissect mechanisms of susceptibility in immunocompetent individuals, we developed a genetically defined sst1-susceptible mouse model that uniquely reproduces a defining feature of human TB: development of necrotic lung lesions after infection with virulent Mtb. In this study, we explored the connectivity of the sst1-regulated pathways during prolonged macrophage activation with TNF. We determined that the aberrant response of the sst1-susceptible macrophages to TNF was primarily driven by conflicting Myc and antioxidant response pathways that resulted in a coordinated failure to properly sequester intracellular iron and activate ferroptosis inhibitor enzymes. Consequently, iron-mediated lipid peroxidation fueled IFNβ superinduction and sustained the Type I Interferon (IFN-I) pathway hyperactivity that locked the sst1-susceptible macrophages in a state of unresolving stress and compromised their resistance to Mtb. The accumulation of the aberrantly activated, stressed, macrophages within granuloma microenvironment led to the local failure of anti-tuberculosis immunity and tissue necrosis. Our findings suggest a novel link between metabolic dysregulation in macrophages and susceptibility to TB, offering insights into potential therapeutic targets aimed at modulating macrophage function and improving TB control.

VEGF counteracts shear stress-determined arterial fate specification during capillary remodeling

A key feature of arteriogenesis is capillary-to-arterial endothelial cell fate transition. Although a number of studies in the past two decades suggested this process is driven by VEGF activation of Notch signaling, how arteriogenesis is regulated remains poorly understood. Here we report that arterial specification is mediated by fluid shear stress (FSS) independent of VEGFR2 signaling and that a decline in VEGFR2 signaling is required for arteriogenesis to fully take place. VEGF does not induce arterial fate in capillary ECs and, instead, counteracts FSS-driven capillary-to-arterial cell fate transition. Mechanistically, FSS-driven arterial program involves both Notch-dependent and Notch-independent events. Sox17 is the key mediator of the FSS-induced arterial specification and a target of VEGF-FSS competition. These findings suggest a new paradigm of VEGF-FSS crosstalk coordinating angiogenesis, arteriogenesis and capillary maintenance.

Reversing pathological cell states: the road less travelled can extend the therapeutic horizon

Acquisition of omics data advances at a formidable pace. Yet, our ability to utilize these data to control cell phenotypes and design interventions that reverse pathological states lags behind. Here, we posit that cell states are determined by core networks that control cell-wide networks. To steer cell fate decisions, core networks connecting genotype to phenotype must be reconstructed and understood. A recent method, cell state transition assessment and regulation (cSTAR), applies perturbation biology to quantify causal connections and mechanistically models how core networks influence cell phenotypes. cSTAR models are akin to digital cell twins enabling us to purposefully convert pathological states back to physiologically normal states. While this capability has a range of applications, here we discuss reverting oncogenic transformation.

cSTAR analysis identifies endothelial cell cycle as a key regulator of flow-dependent artery remodeling

Fluid shear stress (FSS) from blood flow is sensed by vascular endothelial cells (ECs) to determine vessel stability, remodeling and susceptibility to atherosclerosis and other inflammatory diseases but the regulatory networks that govern these behaviors are only partially understood. We used cSTAR, a powerful new computational method, to define EC transcriptomic states under low shear stress (LSS) that triggers vessel inward remodeling, physiological shear stress (PSS) that stabilizes vessels, high shear stress (HSS) that triggers outward remodeling, and oscillatory shear stress (OSS) that confers disease susceptibility, all in comparison to cells under static conditions (STAT). We combined these results with the LINCS database where EC transcriptomic responses to drug treatments to define a preliminary regulatory network in which the cyclin-dependent kinases CDK1/2 play a central role in promoting vessel stability. Experimental analysis showed that PSS induced a strong late G1 cell cycle arrest in which CDK2 was activated. EC deletion of CDK2 in mice resulted in inward artery remodeling and both pulmonary and systemic hypertension. These results validate use of cSTAR to determine EC state and in vivo vessel behavior, reveal unexpected features of EC phenotype under different FSS conditions, and identify CDK2 as a key element within the EC regulatory network that governs artery remodeling.

Cell state transition analysis identifies interventions that improve control of Mycobacterium tuberculosis infection by susceptible macrophages

Understanding cell state transitions and purposefully controlling them to improve therapies is a longstanding challenge in biological research and medicine. Here, we identify a transcriptional signature that distinguishes activated macrophages from the tuberculosis (TB) susceptible and resistant mice. We then apply the cSTAR (cell state transition assessment and regulation) approach to data from screening-by-RNA sequencing to identify chemical perturbations that shift the transcriptional state of tumor necrosis factor (TNF)-activated TB-susceptible macrophages toward that of TB-resistant cells, i.e., prevents their aberrant activation without suppressing beneficial TNF responses. Last, we demonstrate that the compounds identified with this approach enhance the resistance of the TB-susceptible mouse macrophages to virulent Mycobacterium tuberculosis.

A Combination of Conformation-Specific RAF Inhibitors Overcome Drug Resistance Brought about by RAF Overexpression

Cancer cells often adapt to targeted therapies, yet the molecular mechanisms underlying adaptive resistance remain only partially understood. Here, we explore a mechanism of RAS/RAF/MEK/ERK (MAPK) pathway reactivation through the upregulation of RAF isoform (RAFs) abundance. Using computational modeling and in vitro experiments, we show that the upregulation of RAFs changes the concentration range of paradoxical pathway activation upon treatment with conformation-specific RAF inhibitors. Additionally, our data indicate that the signaling output upon loss or downregulation of one RAF isoform can be compensated by overexpression of other RAF isoforms. We furthermore demonstrate that, while single RAF inhibitors cannot efficiently inhibit ERK reactivation caused by RAF overexpression, a combination of two structurally distinct RAF inhibitors synergizes to robustly suppress pathway reactivation.

Control of cell state transitions

Understanding cell state transitions and purposefully controlling them is a longstanding challenge in biology. Here we present cell state transition assessment and regulation (cSTAR), an approach for mapping cell states, modelling transitions between them and predicting targeted interventions to convert cell fate decisions. cSTAR uses omics data as input, classifies cell states, and develops a workflow that transforms the input data into mechanistic models that identify a core signalling network, which controls cell fate transitions by influencing whole-cell networks. By integrating signalling and phenotypic data, cSTAR models how cells manoeuvre in Waddington's landscape1 and make decisions about which cell fate to adopt. Notably, cSTAR devises interventions to control the movement of cells in Waddington's landscape. Testing cSTAR in a cellular model of differentiation and proliferation shows a high correlation between quantitative predictions and experimental data. Applying cSTAR to different types of perturbation and omics datasets, including single-cell data, demonstrates its flexibility and scalability and provides new biological insights. The ability of cSTAR to identify targeted perturbations that interconvert cell fates will enable designer approaches for manipulating cellular development pathways and mechanistically underpinned therapeutic interventions.

Can Systems Biology Advance Clinical Precision Oncology?

Precision oncology is perceived as a way forward to treat individual cancer patients. However, knowing particular cancer mutations is not enough for optimal therapeutic treatment, because cancer genotype-phenotype relationships are nonlinear and dynamic. Systems biology studies the biological processes at the systems' level, using an array of techniques, ranging from statistical methods to network reconstruction and analysis, to mathematical modeling. Its goal is to reconstruct the complex and often counterintuitive dynamic behavior of biological systems and quantitatively predict their responses to environmental perturbations. In this paper, we review the impact of systems biology on precision oncology. We show examples of how the analysis of signal transduction networks allows to dissect resistance to targeted therapies and inform the choice of combinations of targeted drugs based on tumor molecular alterations. Patient-specific biomarkers based on dynamical models of signaling networks can have a greater prognostic value than conventional biomarkers. These examples support systems biology models as valuable tools to advance clinical and translational oncological research.

Channeling macrophage polarization by rocaglates increases macrophage resistance to Mycobacterium tuberculosis

Macrophages contribute to host immunity and tissue homeostasis via alternative activation programs. M1-like macrophages control intracellular bacterial pathogens and tumor progression. In contrast, M2-like macrophages shape reparative microenvironments that can be conducive for pathogen survival or tumor growth. An imbalance of these macrophages phenotypes may perpetuate sites of chronic unresolved inflammation, such as infectious granulomas and solid tumors. We have found that plant-derived and synthetic rocaglates sensitize macrophages to low concentrations of the M1-inducing cytokine IFN-gamma and inhibit their responsiveness to IL-4, a prototypical activator of the M2-like phenotype. Treatment of primary macrophages with rocaglates enhanced phagosome-lysosome fusion and control of intracellular mycobacteria. Thus, rocaglates represent a novel class of immunomodulators that can direct macrophage polarization toward the M1-like phenotype in complex microenvironments associated with hypofunction of type 1 and/or hyperactivation of type 2 immunity, e.g., chronic bacterial infections, allergies, and, possibly, certain tumors.

Modeling the Nonlinear Dynamics of Intracellular Signaling Networks

This protocol illustrates a pipeline for modeling the nonlinear behavior of intracellular signaling pathways. At fixed spatial points, nonlinear signaling dynamics are described by ordinary differential equations (ODEs). At constant parameters, these ODEs may have multiple attractors, such as multiple steady states or limit cycles. Standard optimization procedures fine-tune the parameters for the system trajectories localized within the basin of attraction of only one attractor, usually a stable steady state. The suggested protocol samples the parameter space and captures the overall dynamic behavior by analyzing the number and stability of steady states and the shapes of the assembly of nullclines, which are determined as projections of quasi-steady-state trajectories into different 2D spaces of system variables. Our pipeline allows identifying main qualitative features of the model behavior, perform bifurcation analysis, and determine the borders separating the different dynamical regimes within the assembly of 2D parametric planes. Partial differential equation (PDE) systems describing the nonlinear spatiotemporal behavior are derived by coupling fixed point dynamics with species diffusion.

Reengineering protein-phosphorylation switches

Phosphorylation circuits operate as logic gates that rapidly toggle a system between two stable states

A systematic analysis of signaling reactivation and drug resistance

Increasing evidence suggests that the reactivation of initially inhibited signaling pathways causes drug resistance. Here, we analyze how network topologies affect signaling responses to drug treatment. Network-dependent drug resistance is commonly attributed to negative and positive feedback loops. However, feedback loops by themselves cannot completely reactivate steady-state signaling. Newly synthesized negative feedback regulators can induce a transient overshoot but cannot fully restore output signaling. Complete signaling reactivation can only occur when at least two routes, an activating and inhibitory, connect an inhibited upstream protein to a downstream output. Irrespective of the network topology, drug-induced overexpression or increase in target dimerization can restore or even paradoxically increase downstream pathway activity. Kinase dimerization cooperates with inhibitor-mediated alleviation of negative feedback. Our findings inform drug development by considering network context and optimizing the design drug combinations. As an example, we predict and experimentally confirm specific combinations of RAF inhibitors that block mutant NRAS signaling.

Kholodenko et al. uncover signaling network circuitries and molecular mechanisms necessary and sufficient for complete reactivation or overshoot of steady-state signaling after kinase inhibitor treatment. The two means to revive signaling output fully are through network topology or reactivation of the kinase activity of the primary drug target. Blocking RAF dimer activity by a combination of type I½ and type II RAF inhibitors efficiently blocks mutant NRAS-driven ERK signaling.

Systems biology approaches to macromolecules: the role of dynamic protein assemblies in information processing

Macromolecular protein assemblies govern many cellular processes and are disturbed in many diseases including cancer. Often seen as static molecular machines, protein complexes involved in signal transduction networks exhibit intricate dynamics that are critical for their function. Using the RAS-RAF-MEK-ERK pathway as example we discuss recent progress in our understanding of protein complex dynamics achieved through mathematical modelling, computational simulations and structural studies. The emerging picture highlights that both spatial and temporal dynamics cooperate to enable correct signal processing and the fine tuning of timing, duration and strengths of signalling. These dynamic processes are subverted by oncogenic mutations and contribute to tumorigenesis and drug resistance.

Periodic propagating waves coordinate RhoGTPase network dynamics at the leading and trailing edges during cell migration

Migrating cells need to coordinate distinct leading and trailing edge dynamics but the underlying mechanisms are unclear. Here, we combine experiments and mathematical modeling to elaborate the minimal autonomous biochemical machinery necessary and sufficient for this dynamic coordination and cell movement. RhoA activates Rac1 via DIA and inhibits Rac1 via ROCK, while Rac1 inhibits RhoA through PAK. Our data suggest that in motile, polarized cells, RhoA–ROCK interactions prevail at the rear, whereas RhoA-DIA interactions dominate at the front where Rac1/Rho oscillations drive protrusions and retractions. At the rear, high RhoA and low Rac1 activities are maintained until a wave of oscillatory GTPase activities from the cell front reaches the rear, inducing transient GTPase oscillations and RhoA activity spikes. After the rear retracts, the initial GTPase pattern resumes. Our findings show how periodic, propagating GTPase waves coordinate distinct GTPase patterns at the leading and trailing edge dynamics in moving cells.

Extensive rewiring of the EGFR network in colorectal cancer cells expressing transforming levels of KRASG13D

Protein-protein-interaction networks (PPINs) organize fundamental biological processes, but how oncogenic mutations impact these interactions and their functions at a network-level scale is poorly understood. Here, we analyze how a common oncogenic KRAS mutation (KRASG13D) affects PPIN structure and function of the Epidermal Growth Factor Receptor (EGFR) network in colorectal cancer (CRC) cells. Mapping >6000 PPIs shows that this network is extensively rewired in cells expressing transforming levels of KRASG13D (mtKRAS). The factors driving PPIN rewiring are multifactorial including changes in protein expression and phosphorylation. Mathematical modelling also suggests that the binding dynamics of low and high affinity KRAS interactors contribute to rewiring. PPIN rewiring substantially alters the composition of protein complexes, signal flow, transcriptional regulation, and cellular phenotype. These changes are validated by targeted and global experimental analysis. Importantly, genetic alterations in the most extensively rewired PPIN nodes occur frequently in CRC and are prognostic of poor patient outcomes.

An integrative computational approach for a prioritization of key transcription regulators associated with nanomaterial-induced toxicity

A rapid increase of new nanomaterial products poses new challenges for their risk assessment. Current traditional methods for estimating potential adverse health effect of nanomaterials (NMs) are complex, time consuming and expensive. In order to develop new prediction tests for nanotoxicity evaluation, a systems biology approach and data from high-throughput omics experiments can be used. We present a computational approach that combines reverse engineering techniques, network analysis and pathway enrichment analysis for inferring the transcriptional regulation landscape and its functional interpretation. To illustrate this approach, we used published transcriptomic data derived from mice lung tissue exposed to carbon nanotubes (NM-401 and NRCWE-26). Because fibrosis is the most common adverse effect of these NMs, we included in our analysis the data for bleomycin (BLM) treatment, which is a well-known fibrosis inducer. We inferred gene regulatory networks for each NM and BLM to capture functional hierarchical regulatory structures between genes and their regulators. Despite the different nature of the lung injury caused by nanoparticles and BLM, we identified several conserved core regulators for all agents. We reason that these regulators can be considered as early predictors of toxic responses after NMs exposure. This integrative approach, which refines traditional methods of transcriptomic analysis, can be useful for prioritization of potential core regulators and generation of new hypothesis about mechanisms of nanoparticles toxicity.

Mapping connections in signaling networks with ambiguous modularity

Modular Response Analysis (MRA) is a suite of methods that under certain assumptions permits the precise reconstruction of both the directions and strengths of connections between network modules from network responses to perturbations. Standard MRA assumes that modules are insulated, thereby neglecting the existence of inter-modular protein complexes. Such complexes sequester proteins from different modules and propagate perturbations to the protein abundance of a downstream module retroactively to an upstream module. MRA-based network reconstruction detects retroactive, sequestration-induced connections when an enzyme from one module is substantially sequestered by its substrate that belongs to a different module. Moreover, inferred networks may surprisingly depend on the choice of protein abundances that are experimentally perturbed, and also some inferred connections might be false. Here, we extend MRA by introducing a combined computational and experimental approach, which allows for a computational restoration of modular insulation, unmistakable network reconstruction and discrimination between solely regulatory and sequestration-induced connections for a range of signaling pathways. Although not universal, our approach extends MRA methods to signaling networks with retroactive interactions between modules arising from enzyme sequestration effects.

Modeling cell line-specific recruitment of signaling proteins to the insulin-like growth factor 1 receptor

Receptor tyrosine kinases (RTKs) typically contain multiple autophosphorylation sites in their cytoplasmic domains. Once activated, these autophosphorylation sites can recruit downstream signaling proteins containing Src homology 2 (SH2) and phosphotyrosine-binding (PTB) domains, which recognize phosphotyrosine-containing short linear motifs (SLiMs). These domains and SLiMs have polyspecific or promiscuous binding activities. Thus, multiple signaling proteins may compete for binding to a common SLiM and vice versa. To investigate the effects of competition on RTK signaling, we used a rule-based modeling approach to develop and analyze models for ligand-induced recruitment of SH2/PTB domain-containing proteins to autophosphorylation sites in the insulin-like growth factor 1 (IGF1) receptor (IGF1R). Models were parameterized using published datasets reporting protein copy numbers and site-specific binding affinities. Simulations were facilitated by a novel application of model restructuration, to reduce redundancy in rule-derived equations. We compare predictions obtained via numerical simulation of the model to those obtained through simple prediction methods, such as through an analytical approximation, or ranking by copy number and/or K_D value, and find that the simple methods are unable to recapitulate the predictions of numerical simulations. We created 45 cell line-specific models that demonstrate how early events in IGF1R signaling depend on the protein abundance profile of a cell. Simulations, facilitated by model restructuration, identified pairs of IGF1R binding partners that are recruited in anti-correlated and correlated fashions, despite no inclusion of cooperativity in our models. This work shows that the outcome of competition depends on the physicochemical parameters that characterize pairwise interactions, as well as network properties, including network connectivity and the relative abundances of competitors.

Cells rely on networks of interacting biomolecules to sense and respond to environmental perturbations and signals. However, it is unclear how information is processed to generate appropriate and specific responses to signals, especially given that these networks tend to share many components. For example, receptors that detect distinct ligands and regulate distinct cellular activities commonly interact with overlapping sets of downstream signaling proteins. Here, to investigate the downstream signaling of a well-studied receptor tyrosine kinase (RTK), the insulin-like growth factor 1 (IGF1) receptor (IGF1R), we formulated and analyzed 45 cell line-specific mathematical models, which account for recruitment of 18 different binding partners to six sites of receptor autophosphorylation in IGF1R. The models were parameterized using available protein copy number and site-specific affinity measurements, and restructured to allow for network generation. We find that recruitment is influenced by the protein abundance profile of a cell, with different patterns of recruitment in different cell lines. Furthermore, in a given cell line, we find that pairs of IGF1R binding partners may be recruited in a correlated or anti-correlated fashion. We demonstrate that the simulations of the model have greater predictive power than protein copy number and/or binding affinity data, and that even a simple analytical model cannot reproduce the predicted recruitment ranking obtained via simulations. These findings represent testable predictions and indicate that the outputs of IGF1R signaling depend on cell line-specific properties in addition to the properties that are intrinsic to the biomolecules involved.

Impact of measurement noise, experimental design, and estimation methods on Modular Response Analysis based network reconstruction

Modular Response Analysis (MRA) is a method to reconstruct signalling networks from steady-state perturbation data which has frequently been used in different settings. Since these data are usually noisy due to multi-step measurement procedures and biological variability, it is important to investigate the effect of this noise onto network reconstruction. Here we present a systematic study to investigate propagation of noise from concentration measurements to network structures. Therefore, we design an in silico study of the MAPK and the p53 signalling pathways with realistic noise settings. We make use of statistical concepts and measures to evaluate accuracy and precision of individual inferred interactions and resulting network structures. Our results allow to derive clear recommendations to optimize the performance of MRA based network reconstruction: First, large perturbations are favorable in terms of accuracy even for models with non-linear steady-state response curves. Second, a single control measurement for different perturbation experiments seems to be sufficient for network reconstruction, and third, we recommend to execute the MRA workflow with the mean of different replicates for concentration measurements rather than using computationally more involved regression strategies.

New insights into RAS biology reinvigorate interest in mathematical modeling of RAS signaling

RAS is the most frequently mutated gene across human cancers, but developing inhibitors of mutant RAS has proven to be challenging. Given the difficulties of targeting RAS directly, drugs that impact the other components of pathways where mutant RAS operates may potentially be effective. However, the system-level features, including different localizations of RAS isoforms, competition between downstream effectors, and interlocking feedback and feed-forward loops, must be understood to fully grasp the opportunities and limitations of inhibiting specific targets. Mathematical modeling can help us discern the system-level impacts of these features in normal and cancer cells. New technologies enable the acquisition of experimental data that will facilitate development of realistic models of oncogenic RAS behavior. In light of the wealth of empirical data accumulated over decades of study and the advancement of experimental methods for gathering new data, modelers now have the opportunity to advance progress toward realization of targeted treatment for mutant RAS-driven cancers.

Transcriptionally inducible Pleckstrin homology-like domain, family A, member 1, attenuates ErbB receptor activity by inhibiting receptor oligomerization

Feedback control is a key mechanism in signal transduction, intimately involved in regulating the outcome of the cellular response. Here, we report a novel mechanism by which PHLDA1, Pleckstrin homology-like domain, family A, member 1, negatively regulates ErbB receptor signaling by inhibition of receptor oligomerization. We have found that the ErbB3 ligand, heregulin, induces PHILDA1 expression in MCF-7 cells. Transcriptionally-induced PHLDA1 protein directly binds to ErbB3, whereas knockdown of PHLDA1 increases complex formation between ErbB3 and ErbB2. To provide insight into the mechanism for our time-course and single-cell experimental observations, we performed a systematic computational search of network topologies of the mathematical models based on receptor dimer-tetramer formation in the ErbB activation processes. Our results indicate that only a model in which PHLDA1 inhibits formation of both dimers and tetramer can explain the experimental data. Predictions made from this model were further validated by single-molecule imaging experiments. Our studies suggest a unique regulatory feature of PHLDA1 to inhibit the ErbB receptor oligomerization process and thereby control the activity of receptor signaling network.

Three-factor models versus time series models: quantifying time-dependencies of interactions between stimuli in cell biology and psychobiology for short longitudinal data

Signal integration determines cell fate on the cellular level, affects cognitive processes and affective responses on the behavioural level, and is likely to be involved in psychoneurobiological processes underlying mood disorders. Interactions between stimuli may subjected to time effects. Time-dependencies of interactions between stimuli typically lead to complex cell responses and complex responses on the behavioural level. We show that both three-factor models and time series models can be used to uncover such time-dependencies. However, we argue that for short longitudinal data the three factor modelling approach is more suitable. In order to illustrate both approaches, we re-analysed previously published short longitudinal data sets. We found that in human embryonic kidney 293 cells cells the interaction effect in the regulation of extracellular signal-regulated kinase (ERK) 1 signalling activation by insulin and epidermal growth factor is subjected to a time effect and dramatically decays at peak values of ERK activation. In contrast, we found that the interaction effect induced by hypoxia and tumour necrosis factor-alpha for the transcriptional activity of the human cyclo-oxygenase-2 promoter in HEK293 cells is time invariant at least in the first 12-h time window after stimulation. Furthermore, we applied the three-factor model to previously reported animal studies. In these studies, memory storage was found to be subjected to an interaction effect of the beta-adrenoceptor agonist clenbuterol and certain antagonists acting on the alpha-1-adrenoceptor / glucocorticoid-receptor system. Our model-based analysis suggests that only if the antagonist drug is administer in a critical time window, then the interaction effect is relevant.

Probing the Heterogeneity of Protein Kinase Activation in Cells by Super-resolution Microscopy

Heterogeneity of mitogen-activated protein kinase (MAPK) activation in genetically identical cells, which occurs in response to epidermal growth factor receptor (EGFR) signaling, remains poorly understood. MAPK cascades integrate signals emanating from different EGFR spatial locations, including the plasma membrane and endocytic compartment. We previously hypothesized that in EGF-stimulated cells the MAPK phosphorylation (pMAPK) level and activity are largely determined by the spatial organization of the EGFR clusters within the cell. For experimental testing of this hypothesis, we used super-resolution microscopy to define EGFR clusters by receptor numbers (N) and average intracluster distances (d). From these data, we predicted the extent of pMAPK with 85% accuracy on a cell-to-cell basis with control data returning 54% accuracy (P < 0.001). For comparison, the prediction accuracy was only 61% (P = 0.382) when the diffraction-limited averaged fluorescence intensity/cluster was used. Large clusters (N ≥ 3) with d > 50 nm were most predictive for pMAPK level in cells. Electron microscopy revealed that these large clusters were primarily localized to the limiting membrane of multivesicular bodies (MVB). Many tighter packed dimers/multimers (d < 50 nm) were found on intraluminal vesicles within MVBs, where they were unlikely to activate MAPK because of the physical separation. Our results suggest that cell-to-cell differences in N and d contain crucial information to predict EGFR-activated cellular pMAPK levels and explain pMAPK heterogeneity in isogenic cells.

Integrating network reconstruction with mechanistic modeling to predict cancer therapies

Signal transduction networks are often rewired in cancer cells. Identifying these alterations will enable more effective cancer treatment. We developed a computational framework that can identify, reconstruct, and mechanistically model these rewired networks from noisy and incomplete perturbation response data and then predict potential targets for intervention. As a proof of principle, we analyzed a perturbation data set targeting epidermal growth factor receptor (EGFR) and insulin-like growth factor 1 receptor (IGF1R) pathways in a panel of colorectal cancer cells. Our computational approach predicted cell line-specific network rewiring. In particular, feedback inhibition of insulin receptor substrate 1 (IRS1) by the kinase p70S6K was predicted to confer resistance to EGFR inhibition, suggesting that disrupting this feedback may restore sensitivity to EGFR inhibitors in colorectal cancer cells. We experimentally validated this prediction with colorectal cancer cell lines in culture and in a zebrafish (Danio rerio) xenograft model.

MAPK kinase signalling dynamics regulate cell fate decisions and drug resistance

The RAS/RAF/MEK/MAPK kinase pathway has been extensively studied for more than 25 years, yet we continue to be puzzled by its intricate dynamic control and plasticity. Different spatiotemporal MAPK dynamics bring about distinct cell fate decisions in normal vs cancer cells and developing organisms. Recent modelling and experimental studies provided novel insights in the versatile MAPK dynamics concerted by a plethora of feedforward/feedback regulations and crosstalk on multiple timescales. Multiple cancer types and various developmental disorders arise from persistent alterations of the MAPK dynamics caused by RAS/RAF/MEK mutations. While a key role of the MAPK pathway in multiple diseases made the development of novel RAF/MEK inhibitors a hot topic of drug development, these drugs have unexpected side-effects and resistance inevitably occurs. We review how RAF dimerization conveys drug resistance and recent breakthroughs to overcome this resistance.

HER2-HER3 dimer quantification by FLIM-FRET predicts breast cancer metastatic relapse independently of HER2 IHC status

Overexpression of HER2 is an important prognostic marker, and the only predictive biomarker of response to HER2-targeted therapies in invasive breast cancer. HER2-HER3 dimer has been shown to drive proliferation and tumor progression, and targeting of this dimer with pertuzumab alongside chemotherapy and trastuzumab, has shown significant clinical utility. The purpose of this study was to accurately quantify HER2-HER3 dimerisation in formalin fixed paraffin embedded (FFPE) breast cancer tissue as a novel prognostic biomarker.FFPE tissues were obtained from patients included in the METABRIC (Molecular Taxonomy of Breast Cancer International Consortium) study. HER2-HER3 dimerisation was quantified using an improved fluorescence lifetime imaging microscopy (FLIM) histology-based analysis. Analysis of 131 tissue microarray cores demonstrated that the extent of HER2-HER3 dimer formation as measured by Förster Resonance Energy Transfer (FRET) determined through FLIM predicts the likelihood of metastatic relapse up to 10 years after surgery (hazard ratio 3.91 (1.61-9.5), p = 0.003) independently of HER2 expression, in a multivariate model. Interestingly there was no correlation between the level of HER2 protein expressed and HER2-HER3 heterodimer formation. We used a mathematical model that takes into account the complex interactions in a network of all four HER proteins to explain this counterintuitive finding.Future utility of this technique may highlight a group of patients who do not overexpress HER2 protein but are nevertheless dependent on the HER2-HER3 heterodimer as driver of proliferation. This assay could, if validated in a group of patients treated with, for instance pertuzumab, be used as a predictive biomarker to predict for response to such targeted therapies.

Protein-protein interactions generate hidden feedback and feed-forward loops to trigger bistable switches, oscillations and biphasic dose-responses

Protein-protein interactions (PPIs) defined as reversible association of two proteins to form a complex, are undoubtedly among the most common interaction motifs featured in cells. Recent large-scale proteomic studies have revealed an enormously complex interactome of the cell, consisting of tens of thousands of PPIs with numerous signalling hubs. PPIs have functional roles in regulating a wide range of cellular processes including signal transduction and post-translational modifications, and de-regulation of PPIs is implicated in many diseases including cancers and neuro-degenerative disorders. Despite the ubiquitous appearance and physiological significance of PPIs, our understanding of the dynamic and functional consequences of these simple motifs remains incomplete, particularly when PPIs occur within large biochemical networks. We employ quantitative, dynamic modelling to computationally analyse salient dynamic features of the PPI motifs and PPI-containing signalling networks varying in topological architecture. Our analyses surprisingly reveal that simple reversible PPI motifs, when being embedded into signalling cascades, could give rise to extremely rich and complex regulatory dynamics in the absence of explicit positive and negative feedback loops. Our work represents a systematic investigation of the dynamic properties of PPIs in signalling networks, and the results shed light on how this simple event may potentiate diverse and intricate behaviours in vivo.

The complexities and versatility of the RAS-to-ERK signalling system in normal and cancer cells

The intricate dynamic control and plasticity of RAS to ERK mitogenic, survival and apoptotic signalling has mystified researches for more than 30 years. Therapeutics targeting the oncogenic aberrations within this pathway often yield unsatisfactory, even undesired results, as in the case of paradoxical ERK activation in response to RAF inhibition. A direct approach of inhibiting single oncogenic proteins misses the dynamic network context governing the network signal processing. In this review, we discuss the signalling behaviour of RAS and RAF proteins in normal and in cancer cells, and the emerging systems-level properties of the RAS-to-ERK signalling network. We argue that to understand the dynamic complexities of this control system, mathematical models including mechanistic detail are required. Looking into the future, these dynamic models will build the foundation upon which more effective, rational approaches to cancer therapy will be developed.

Phosphorylation of RAF Kinase Dimers Drives Conformational Changes that Facilitate Transactivation

RAF kinases are key players in the MAPK signaling pathway and are important targets for personalized cancer therapy. RAF dimerization is part of the physiological activation mechanism, together with phosphorylation, and is known to convey resistance to RAF inhibitors. Herein, molecular dynamics simulations are used to show that phosphorylation of a key N-terminal acidic (NtA) motif facilitates RAF dimerization by introducing several interprotomer salt bridges between the αC-helix and charged residues upstream of the NtA motif. Additionally, we show that the R-spine of RAF interacts with a conserved Trp residue in the vicinity of the NtA motif, connecting the active sites of two protomers and thereby modulating the cooperative interactions in the RAF dimer. Our findings provide a first structure-based mechanism for the auto-transactivation of RAF and could be generally applicable to other kinases, opening new pathways for overcoming dimerization-related drug resistance.

Abstract CN05-03: Personalized cancer diagnostics and therapeutics based on the computational modeling of signal transduction networks

Signalling pathways exert finely tuned control over cell fate decisions that ultimately determine the behaviour of cancer cells. It could therefore be expected that the dynamics of key pathway activation may contain prognostically relevant information over and above that which is contained in the static nature of traditional biomarkers. To investigate this hypothesis we have characterised the network architecture regulating JNK stress signalling in neuroblastoma cells and applied an experimentally calibrated and validated computational model of this network to extract prognostic information from neuroblastoma patient-specific simulations of JNK activation. Survival analysis based upon the dynamics of these simulations revealed that an inability to initiate switch-like JNK signalling in silico was significantly associated with poor overall survival for both MYCN amplified and non-amplified neuroblastoma patients. Furthermore, our analysis demonstrated that in order to extract prognostic information from a signalling pathway, deciphering the extant network structure is a vital consideration in model development. We also show that network based analysis can lead to the discovery of new therapeutic targets.

Citation Format: Walter Kolch, Dirk Fey, Melinda Halasz, Nora Rauch, Amaya Garcia Munoz, Ruth Pilkington, Boris N. Kholodenko, David R. Croucher, Sean P. Kennedy, Jordan F. Hastings, Frank Westermann, Daniel Dreidax, Matthias Fischer, David Duffy, Aleksandar Krstic, Thomas Schwarzl. Personalized cancer diagnostics and therapeutics based on the computational modeling of signal transduction networks. [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2015 Nov 5-9; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2015;14(12 Suppl 2):Abstract nr CN05-03.

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.

MST2-RASSF protein-protein interactions through SARAH domains

The detailed, atomistic-level understanding of molecular signaling along the tumor-suppressive Hippo signaling pathway that controls tissue homeostasis by balancing cell proliferation and death through apoptosis is a promising avenue for the discovery of novel anticancer drug targets. The activation of kinases such as Mammalian STE20-Like Protein Kinases 1 and 2 (MST1 and MST2)-modulated through both homo- and heterodimerization (e.g. interactions with Ras association domain family, RASSF, enzymes)-is a key upstream event in this pathway and remains poorly understood. On the other hand, RASSFs (such as RASSF1A or RASSF5) act as important apoptosis activators and tumor suppressors, although their exact regulatory roles are also unclear. We present recent molecular studies of signaling along the Ras-RASSF-MST pathway, which controls growth and apoptosis in eukaryotic cells, including a variety of modern molecular modeling and simulation techniques. Using recently available structural information, we discuss the complex regulatory scenario according to which RASSFs perform dual signaling functions, either preventing or promoting MST2 activation, and thus control cell apoptosis. Here, we focus on recent studies highlighting the special role being played by the specific interactions between the helical Salvador/RASSF/Hippo (SARAH) domains of MST2 and RASSF1a or RASSF5 enzymes. These studies are crucial for integrating atomistic-level mechanistic information about the structures and conformational dynamics of interacting proteins, with information available on their system-level functions in cellular signaling.

Competing to coordinate cell fate decisions: the MST2-Raf-1 signaling device

How do biochemical signaling pathways generate biological specificity? This question is fundamental to modern biology, and its enigma has been accentuated by the discovery that most proteins in signaling networks serve multifunctional roles. An answer to this question may lie in analyzing network properties rather than individual traits of proteins in order to elucidate design principles of biochemical networks that enable biological decision-making. We discuss how this is achieved in the MST2/Hippo-Raf-1 signaling network with the help of mathematical modeling and model-based analysis, which showed that competing protein interactions with affinities controlled by dynamic protein modifications can function as Boolean computing devices that determine cell fate decisions. In addition, we discuss areas of interest for future research and highlight how systems approaches would be of benefit.

Drug Resistance Resulting from Kinase Dimerization Is Rationalized by Thermodynamic Factors Describing Allosteric Inhibitor Effects

Treatment of cancer patients with ATP-competitive inhibitors of BRAF/CRAF kinases surprisingly increases total kinase activity, especially in wild-type BRAF cells, subverting the desired clinical outcome. Similar inhibition resistance is observed for numerous kinases involving homo/heterodimerization in their activation cycles. Here, I demonstrate that drug resistance resulting from kinase dimerization can be explained using thermodynamic principles. I show that allosteric regulation by inhibitors is described by thermodynamic factors that quantify inhibitor-induced changes in kinase dimerization and the difference in the drug affinity for a free monomer versus a dimer harboring one drug molecule. The analysis extends to kinase homo- and heterodimers, allows for their symmetric and asymmetric conformations, and predicts how thermodynamic factors influence dose-response dependencies. I show how two inhibitors, ineffective on their own, when combined can abolish drug resistance at lower doses than either inhibitor applied alone. Thus, the mechanistic models suggest ways to overcome resistance to kinase inhibitors.

The dynamic control of signal transduction networks in cancer cells

Cancer is often considered a genetic disease. However, much of the enormous plasticity of cancer cells to evolve different phenotypes, to adapt to challenging microenvironments and to withstand therapeutic assaults is encoded by the structure and spatiotemporal dynamics of signal transduction networks. In this Review, we discuss recent concepts concerning how the rich signalling dynamics afforded by these networks are regulated and how they impinge on cancer cell proliferation, survival, invasiveness and drug resistance. Understanding this dynamic circuitry by mathematical modelling could pave the way to new therapeutic approaches and personalized treatments.

G Protein-Coupled Receptor Signaling Networks from a Systems Perspective

The signal-transduction network of a mammalian cell integrates internal and external cues to initiate adaptive responses. Among the cell-surface receptors are the G protein-coupled receptors (GPCRs), which have remarkable signal-integrating capabilities. Binding of extracellular signals stabilizes intracellular-domain conformations that selectively activate intracellular proteins. Hereby, multiple signaling routes are activated simultaneously to degrees that are signal-combination dependent. Systems-biology studies indicate that signaling networks have emergent processing capabilities that go far beyond those of single proteins. Such networks are spatiotemporally organized and capable of gradual, oscillatory, all-or-none, and subpopulation-generating responses. Protein-protein interactions, generating feedback and feedforward circuitry, are generally required for these spatiotemporal phenomena. Understanding of information processing by signaling networks therefore requires network theories in addition to biochemical and biophysical concepts. Here we review some of the key signaling systems behaviors that have been discovered recurrently across signaling networks. We emphasize the role of GPCRs, so far underappreciated receptors in systems-biology research.

Mitogen-Inducible Gene-6 Mediates Feedback Inhibition from Mutated BRAF towards the Epidermal Growth Factor Receptor and Thereby Limits Malignant Transformation

BRAF functions in the RAS-extracellular signal-regulated kinase (ERK) signaling cascade. Activation of this pathway is necessary to mediate the transforming potential of oncogenic BRAF, however, it may also cause a negative feedback that inhibits the epidermal growth factor receptor (EGFR). Mitogen-inducible gene-6 (MIG-6) is a potent inhibitor of the EGFR and has been demonstrated to function as a tumor suppressor. As MIG-6 can be induced via RAS-ERK signaling, we investigated its potential involvement in this negative regulatory loop. Focus formation assays were performed and demonstrated that MIG-6 significantly reduces malignant transformation induced by oncogenic BRAF. Although this genetic interaction was mirrored by a physical interaction between MIG-6 and BRAF, we did not observe a direct regulation of BRAF kinase activity by MIG-6. Interestingly, a selective chemical EGFR inhibitor suppressed transformation to a similar degree as MIG-6, whereas combining these approaches had no synergistic effect. By analyzing a range of BRAF mutated and wildtype cell line models, we could show that BRAF V600E causes a strong upregulation of MIG-6, which was mediated at the transcriptional level via the RAS-ERK pathway and resulted in downregulation of EGFR activation. This feedback loop is operational in tumors, as shown by the analysis of almost 400 patients with papillary thyroid cancer (PTC). Presence of BRAF V600E correlated with increased MIG-6 expression on the one hand, and with inactivation of the EGFR and of PI3K/AKT signaling on the other hand. Importantly, we also observed a more aggressive disease phenotype when BRAF V600E coexisted with low MIG-6 expression. Finally, analysis of methylation data was performed and revealed that higher methylation of MIG-6 correlated to its decreased expression. Taken together, we demonstrate that MIG-6 efficiently reduces cellular transformation driven by oncogenic BRAF by orchestrating a negative feedback circuit directed towards the EGFR.

Ubiquitin chain specific auto-ubiquitination triggers sustained oscillation, bistable switches and excitable firing

Ubiquitin modification of cellular proteins commonly targets them for proteosomal degradation, but can also convey non‐proteolytic functions. Over the past years, advances in experimental approaches have helped uncover the extensive involvement of ubiquitination in protein regulation. However, our understanding of the dynamics of the ubiquitination‐related networks have lagged behind. A common regulatory theme for many E3 ligases is the ability to self‐catalyse their own ubiquitination without involving external E3 ligating enzymes. Here, the authors have explored computational models of both proteolytic and non‐proteolytic auto‐ubiquitination of E3 ligases and characterised the dynamic properties of these regulatory motifs. Remarkably, in both cases auto‐ubiquitination coupled with multi‐step de‐ubiquitination process can bring about sustained oscillatory behaviour. In addition, the same basic wiring structures can trigger bistable switches of activity and excitable firing of the dynamic responses of the ubiquitinated E3 ligase. Bifurcation analysis allows one to derive parametric conditions that govern these dynamics. They also show that these complex non‐linear behaviours persist for a more detailed mechanistic description that involves the E1 and E2 enzymes. Their work therefore provides new insights into the dynamic features of auto‐ubiquitination in different cellular contexts.

Advances in dynamic modeling of colorectal cancer signaling-network regions, a path toward targeted therapies

The interconnected network of pathways downstream of the TGFβ, WNT and EGF-families of receptor ligands play an important role in colorectal cancer pathogenesis.We studied and implemented dynamic simulations of multiple downstream pathways and described the section of the signaling network considered as a Molecular Interaction Map (MIM). Our simulations used Ordinary Differential Equations (ODEs), which involved 447 reactants and their interactions.Starting from an initial "physiologic condition", the model can be adapted to simulate individual pathologic cancer conditions implementing alterations/mutations in relevant onco-proteins. We verified some salient model predictions using the mutated colorectal cancer lines HCT116 and HT29. We measured the amount of MYC and CCND1 mRNAs and AKT and ERK phosphorylated proteins, in response to individual or combination onco-protein inhibitor treatments. Experimental and simulation results were well correlated. Recent independently published results were also predicted by our model.Even in the presence of an approximate and incomplete signaling network information, a predictive dynamic modeling seems already possible. An important long term road seems to be open and can be pursued further, by incremental steps, toward even larger and better parameterized MIMs. Personalized treatment strategies with rational associations of signaling-proteins inhibitors, could become a realistic goal.

Nonlinear signalling networks and cell-to-cell variability transform external signals into broadly distributed or bimodal responses

We show theoretically and experimentally a mechanism behind the emergence of wide or bimodal protein distributions in biochemical networks with nonlinear input-output characteristics (the dose-response curve) and variability in protein abundance. Large cell-to-cell variation in the nonlinear dose-response characteristics can be beneficial to facilitate two distinct groups of response levels as opposed to a graded response. Under the circumstances that we quantify mathematically, the two distinct responses can coexist within a cellular population, leading to the emergence of a bimodal protein distribution. Using flow cytometry, we demonstrate the appearance of wide distributions in the hypoxia-inducible factor-mediated response network in HCT116 cells. With help of our theoretical framework, we perform a novel calculation of the magnitude of cell-to-cell heterogeneity in the dose-response obtained experimentally.

Species differential regulation of COX2 can be described by an NFκB-dependent logic AND gate

Cyclooxygenase 2 (COX2), a key regulatory enzyme of the prostaglandin/eicosanoid pathway, is an important target for anti-inflammatory therapy. It is highly induced by pro-inflammatory cytokines in a Nuclear factor kappa B (NFκB)-dependent manner. However, the mechanisms determining the amplitude and dynamics of this important pro-inflammatory event are poorly understood. Furthermore, there is significant difference between human and mouse COX2 expression in response to the inflammatory stimulus tumor necrosis factor alpha (TNFα). Here, we report the presence of a molecular logic AND gate composed of two NFκB response elements (NREs) which controls the expression of human COX2 in a switch-like manner. Combining quantitative kinetic modeling and thermostatistical analysis followed by experimental validation in iterative cycles, we show that the human COX2 expression machinery regulated by NFκB displays features of a logic AND gate. We propose that this provides a digital, noise-filtering mechanism for a tighter control of expression in response to TNFα, such that a threshold level of NFκB activation is required before the promoter becomes active and initiates transcription. This NFκB-regulated AND gate is absent in the mouse COX2 promoter, most likely contributing to its differential graded response in promoter activity and protein expression to TNFα. Our data suggest that the NFκB-regulated AND gate acts as a novel mechanism for controlling the expression of human COX2 to TNFα, and its absence in the mouse COX2 provides the foundation for further studies on understanding species-specific differential gene regulation.

Feedforward regulation of mRNA stability by prolonged extracellular signal-regulated kinase activity

Extracellular signal-regulated kinase (ERK) plays a central role in signal transduction networks and cell fate decisions. Sustained ERK activation induces cell differentiation, whereas transient ERK results in the proliferation of several types of cells. Sustained ERK activity stabilizes the proteins of early-response gene products. However, the effect of ERK activity duration on mRNA stability is unknown. We analyzed the quantitative relationship between the duration of four ERK activity kinetics and the mRNA expression profile in growth factor-treated cells. Time-course transcriptome analysis revealed that the cells with prolonged ERK activity generally showed sustained mRNA expression of late response genes but not early or mid genes. Selected late response genes decayed more rapidly in the presence of a specific ERK inhibitor than a general transcription inhibitor and the decay rate was not related to the number of AU-rich elements. Our results suggest that sustained ERK activity plays an important role in the lifespan of the mRNA encoded by late response genes, in addition to the previously demonstrated role in protein stabilization of early-response genes, including transcription factors regulating the transcription of mid and late genes. This double-positive regulation of ligand-induced genes, also termed feedforward regulation, is critical in cell fate decisions.

The ErbB4 CYT2 variant protects EGFR from ligand-induced degradation to enhance cancer cell motility

The epidermal growth factor receptor (EGFR) is a member of the ErbB family that can promote the migration and proliferation of breast cancer cells. Therapies that target EGFR can promote the dimerization of EGFR with other ErbB receptors, which is associated with the development of drug resistance. Understanding how interactions among ErbB receptors alter EGFR biology could provide avenues for improving cancer therapy. We found that EGFR interacted directly with the CYT1 and CYT2 variants of ErbB4 and the membrane-anchored intracellular domain (mICD). The CYT2 variant, but not the CYT1 variant, protected EGFR from ligand-induced degradation by competing with EGFR for binding to a complex containing the E3 ubiquitin ligase c-Cbl and the adaptor Grb2. Cultured breast cancer cells overexpressing both EGFR and ErbB4 CYT2 mICD exhibited increased migration. With molecular modeling, we identified residues involved in stabilizing the EGFR dimer. Mutation of these residues in the dimer interface destabilized the complex in cells and abrogated growth factor-stimulated cell migration. An exon array analysis of 155 breast tumors revealed that the relative mRNA abundance of the ErbB4 CYT2 variant was increased in ER+ HER2- breast cancer patients, suggesting that our findings could be clinically relevant. We propose a mechanism whereby competition for binding to c-Cbl in an ErbB signaling heterodimer promotes migration in response to a growth factor gradient.