Cell State Transition Models Stratify Breast Cancer Cell Phenotypes and Reveal New Therapeutic Targets.

Understanding signaling patterns of transformation and controlling cell phenotypes is a challenge of current biology. Here we applied a cell State Transition Assessment and Regulation (cSTAR) approach to a perturbation dataset of single cell phosphoproteomic patterns of multiple breast cancer (BC) and normal breast tissue-derived cell lines. Following a separation of luminal, basal, and normal cell states, we identified signaling nodes within core control networks, delineated causal connections, and determined the primary drivers underlying oncogenic transformation and transitions across distinct BC subtypes. Whereas cell lines within the same BC subtype have different mutational and expression profiles, the architecture of the core network was similar for all luminal BC cells, and mTOR was a main oncogenic driver. In contrast, core networks of basal BC were heterogeneous and segregated into roughly four major subclasses with distinct oncogenic and BC subtype drivers. Likewise, normal breast tissue cells were separated into two different subclasses. Based on the data and quantified network topologies, we derived mechanistic cSTAR models that serve as digital cell twins and allow the deliberate control of cell movements within a Waddington landscape across different cell states. These cSTAR models suggested strategies of normalizing phosphorylation networks of BC cell lines using small molecule inhibitors.

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.

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.

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.

Cell state transition analysis identifies interventions that improve control of M. 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 TB-susceptible and TB-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 the TB-susceptible macrophages towards that of TB-resistant cells. Finally, we demonstrate that the compounds identified with this approach enhance resistance of the TB-susceptible mouse macrophages to virulent M. tuberculosis .

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.

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.

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.

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.

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.

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 KD 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.

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.

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.

Dissecting RAF Inhibitor Resistance by Structure-based Modeling Reveals Ways to Overcome Oncogenic RAS Signaling.

Clinically used RAF inhibitors are ineffective in RAS mutant tumors because they enhance homo- and heterodimerization of RAF kinases, leading to paradoxical activation of ERK signaling. Overcoming enhanced RAF dimerization and the resulting resistance is a challenge for drug design. Combining multiple inhibitors could be more effective, but it is unclear how the best combinations can be chosen. We built a next-generation mechanistic dynamic model to analyze combinations of structurally different RAF inhibitors, which can efficiently suppress MEK/ERK signaling. This rule-based model of the RAS/ERK pathway integrates thermodynamics and kinetics of drug-protein interactions, structural elements, posttranslational modifications, and cell mutational status as model rules to predict RAF inhibitor combinations for inhibiting ERK activity in oncogenic RAS and/or BRAFV600E backgrounds. Predicted synergistic inhibition of ERK signaling was corroborated by experiments in mutant NRAS, HRAS, and BRAFV600E cells, and inhibition of oncogenic RAS signaling was associated with reduced cell proliferation and colony formation.

On the chemotherapeutic agents localization in tissue by means of snake venoms.

The efficiency of anti-tumour drug strongly depends on its dose. Higher drug doses and exposure times usually result in better treatment. It is why the implementation of high-dose treatment is always attractive. However, most of the drug delivery techniques meet essential limitations. In isolated regional perfusion a tumour can be exposed to high-dose therapeutic influence but the target organ may be isolated from the rest of circulatory system only for a relatively short period of time. During systemic injection of anti-tumour agents dose limitations are dictated by side toxicity danger. Viperidae venoms are known to cause local stagnation of blood flow and blood-tissue exchange processes in the place of snakebite. In present paper we suggest to use Viperidae snake venoms in addition to anti-tumour drugs for regional anti-cancer therapy. We suppose that Viperidae venoms will assist in drug localization. We state that their usage will help in high-dosage therapy implementation.

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.

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.