Combinatorial interventions inhibit TGFβ-driven epithelial-to-mesenchymal transition and support hybrid cellular phenotypes

Output Stabilizing Control of Complex Biological Networks Based on Boolean Algebra Analysis

Output stabilizing control of biological systems is of utmost importance in systems biology since key phenotypes of biological networks are often encoded by a small subset of their phenotypic marker nodes. This study addresses the challenge of output stabilizing control for complex biological systems modeled by Boolean networks (BNs). The objective is to identify a set of constant control inputs capable of driving the BN toward a desirable long-term behavior with respect to specified output nodes. Leveraging the algebraic properties of Boolean logic, we develop a novel control algorithm that reformulates the output stabilizing control problem into a simple graph theoretic problem involving auxiliary BNs, the scale of which significantly decreases compared to the original BN. The proposed method ensures superiority over previous results in terms of both the number of control inputs and computational loads, since it searches for the solution within the reduced BNs while retaining essential structures needed for output stabilization. The efficacy of the proposed control scheme is demonstrated through extensive numerical experiments with complex random BNs and real biological networks. To support the reproducible research initiative, detailed results of numerical experiments are provided in the supplementary material, and all the implementation codes are made accessible at https://github.com/choonlog/OutputStabilization.

Systems Biology of the Cancer Cell

Questions in cancer have engaged systems biologists for decades. During that time, the quantity of molecular data has exploded, but the need for abstractions, formal models, and simplifying insights has remained the same. This review brings together classic breakthroughs and recent findings in the field of cancer systems biology, focusing on cancer cell pathways for tumorigenesis and therapeutic response. Cancer cells mutate and transduce information from their environment to alter gene expression, metabolism, and phenotypic states. Understanding the molecular architectures that make each of these steps possible is a long-term goal of cancer systems biology pursued by iterating between quantitative models and experiments. We argue that such iteration is the best path to deploying targeted therapies intelligently so that each patient receives the maximum benefit for their cancer.

CanSeer: a translational methodology for developing personalized cancer models and therapeutics

Computational modeling and analysis of biomolecular network models annotated with omics data are emerging as a versatile tool for designing personalized therapies. Current endeavors aimed at employing in silico models towards personalized cancer therapeutics remain limited in providing all-in-one approach that ascertains actionable targets, re-positions FDA (Food and Drug Administration) approved drugs, furnishes quantitative cues on therapy responses such as efficacy and cytotoxic effect, and identifies novel drug combinations. Here we propose "CanSeer"-a methodology for developing personalized therapeutics. CanSeer employs patient-specific genetic alterations and RNA-seq data to annotate in silico models followed by dynamical network analyses towards assessment of treatment responses. To exemplify, three use cases involving paired samples, unpaired samples, and cancer samples only, of lung squamous cell carcinoma (LUSC) patients are provided. CanSeer reveals the effectiveness of repositioned drugs along with the identification of several novel LUSC treatment combinations including Afuresertib + Palbociclib, Dinaciclib + Trametinib, Afatinib + Oxaliplatin, Ulixertinib + Olaparib, etc.

A Boolean network model of hypoxia, mechanosensing and TGF-β signaling captures the role of phenotypic plasticity and mutations in tumor metastasis

The tumor microenvironment aids cancer progression by promoting several cancer hallmarks, independent of cancer-related mutations. Biophysical properties of this environment, such as the stiffness of the matrix cells adhere to and local cell density, impact proliferation, apoptosis, and the epithelial to mesenchymal transition (EMT). The latter is a rate-limiting step for invasion and metastasis, enhanced in hypoxic tumor environments but hindered by soft matrices and/or high cell densities. As these influences are often studied in isolation, the crosstalk between hypoxia, biomechanical signals, and the classic EMT driver TGF-β is not well mapped, limiting our ability to predict and anticipate cancer cell behaviors in changing tumor environments. To address this, we built a Boolean regulatory network model that integrates hypoxic signaling with a mechanosensitive model of EMT, which includes the EMT-promoting crosstalk of mitogens and biomechanical signals, cell cycle control, and apoptosis. Our model reproduces the requirement of Hif-1α for proliferation, the anti-proliferative effects of strong Hif-1α stabilization during hypoxia, hypoxic protection from anoikis, and hypoxia-driven mechanosensitive EMT. We offer experimentally testable predictions about the effect of VHL loss on cancer hallmarks, with or without secondary oncogene activation. Taken together, our model serves as a predictive framework to synthesize the signaling responses associated with tumor progression and metastasis in healthy vs. mutant cells. Our single-cell model is a key step towards more extensive regulatory network models that cover damage-response and senescence, integrating most cell-autonomous cancer hallmarks into a single model that can, in turn, control the behavior of in silico cells within a tissue model of epithelial homeostasis and carcinoma.

An integrative phenotype-structured partial differential equation model for the population dynamics of epithelial-mesenchymal transition

Phenotypic heterogeneity along the epithelial-mesenchymal (E-M) axis contributes to cancer metastasis and drug resistance. Recent experimental efforts have collated detailed time-course data on the emergence and dynamics of E-M heterogeneity in a cell population. However, it remains unclear how different intra- and inter-cellular processes shape the dynamics of E-M heterogeneity. Here, using Cell Population Balance model, we capture the dynamics of cell density along E-M phenotypic axis resulting from interplay between-(a) intracellular regulatory interaction among biomolecules, (b) cell division and death and (c) stochastic cell-state transition. We find that while the existence of E-M heterogeneity depends on intracellular regulation, heterogeneity gets enhanced with stochastic cell-state transitions and diminished by growth rate differences. Further, resource competition among E-M cells can lead to both bi-phasic growth of the total population and/or bi-stability in the phenotypic composition. Overall, our model highlights complex interplay between cellular processes shaping dynamic patterns of E-M heterogeneity.

Gene regulatory network inference during cell fate decisions by perturbation strategies

With rapid advances in biological technology and computational approaches, inferring specific gene regulatory networks from data alone during cell fate decisions, including determining direct regulations and their intensities between biomolecules, remains one of the most significant challenges. In this study, we propose a general computational approach based on systematic perturbation, statistical, and differential analyses to infer network topologies and identify network differences during cell fate decisions. For each cell fate state, we first theoretically show how to calculate local response matrices based on perturbation data under systematic perturbation analysis, and we also derive the wild-type (WT) local response matrix for specific ordinary differential equations. To make the inferred network more accurate and eliminate the impact of perturbation degrees, the confidence interval (CI) of local response matrices under multiple perturbations is applied, and the redefined local response matrix is proposed in statistical analysis to determine network topologies across all cell fates. Then in differential analysis, we introduce the concept of relative local response matrix, which enables us to identify critical regulations governing each cell state and dominant cell states associated with specific regulations. The epithelial to mesenchymal transition (EMT) network is chosen as an illustrative example to verify the feasibility of the approach. Largely consistent with experimental observations, the differences of inferred networks at the three cell states can be quantitatively identified. The approach presented here can be also applied to infer other regulatory networks related to cell fate decisions.

Fibrotic extracellular matrix preferentially induces a partial Epithelial-Mesenchymal Transition phenotype in a 3-D agent based model of fibrosis

One of the main drivers of fibrotic diseases is epithelial-mesenchymal transition (EMT): a transdifferentiation process in which cells undergo a phenotypic change from an epithelial state to a pro-migratory state. The cytokine transforming growth factor-β1 (TGF-β1) has been previously shown to regulate EMT. TGF-β1 binds to fibronectin (FN) fibrils, which are the primary extracellular matrix (ECM) component in renal fibrosis. We have previously demonstrated experimentally that inhibition of FN fibrillogenesis and/or TGF-β1 tethering to FN inhibits EMT. However, these studies have only been conducted on 2-D cell monolayers, and the role of TGF-β1-FN tethering in 3-D cellular environments is not clear. As such, we sought to develop a 3-D computational model of epithelial spheroids that captured both EMT signaling dynamics and TGF-β1-FN tethering dynamics. We have incorporated the bi-stable EMT switch model developed by Tian et al. (2013) into a 3-D multicellular model to capture both temporal and spatial TGF-β1 signaling dynamics. We showed that the addition of increasing concentrations of exogeneous TGF-β1 led to faster EMT progression, indicated by increased expression of mesenchymal markers, decreased cell proliferation and increased migration. We then incorporated TGF-β1-FN fibril tethering by locally reducing the TGF-β1 diffusion coefficient as a function of EMT to simulate the reduced movement of TGF-β1 when tethered to FN fibrils during fibrosis. We showed that incorporation of TGF-β1 tethering to FN fibrils promoted a partial EMT state, independent of exogenous TGF-β1 concentration, indicating a mechanism by which fibrotic ECM can promote a partial EMT state.

Low dimensionality of phenotypic space as an emergent property of coordinated teams in biological regulatory networks

Cell-fate decisions involve coordinated genome-wide expression changes, typically leading to a limited number of phenotypes. Although often modeled as simple toggle switches, these rather simplistic representations often disregard the complexity of regulatory networks governing these changes. Here, we unravel design principles underlying complex cell decision-making networks in multiple contexts. We show that the emergent dynamics of these networks and corresponding transcriptomic data are consistently low-dimensional, as quantified by the variance explained by principal component 1 (PC1). This low dimensionality in phenotypic space arises from extensive feedback loops in these networks arranged to effectively enable the formation of two teams of mutually inhibiting nodes. We use team strength as a metric to quantify these feedback interactions and show its strong correlation with PC1 variance. Using artificial networks of varied topologies, we also establish the conditions for generating canalized cell-fate landscapes, offering insights into diverse binary cellular decision-making networks.

AI-driven automated discovery tools reveal diverse behavioral competencies of biological networks

Many applications in biomedicine and synthetic bioengineering rely on understanding, mapping, predicting, and controlling the complex behavior of chemical and genetic networks. The emerging field of diverse intelligence investigates the problem-solving capacities of unconventional agents. However, few quantitative tools exist for exploring the competencies of non-conventional systems. Here, we view gene regulatory networks (GRNs) as agents navigating a problem space and develop automated tools to map the robust goal states GRNs can reach despite perturbations. Our contributions include: (1) Adapting curiosity-driven exploration algorithms from AI to discover the range of reachable goal states of GRNs, and (2) Proposing empirical tests inspired by behaviorist approaches to assess their navigation competencies. Our data shows that models inferred from biological data can reach a wide spectrum of steady states, exhibiting various competencies in physiological network dynamics without requiring structural changes in network properties or connectivity. We also explore the applicability of these 'behavioral catalogs' for comparing evolved competencies across biological networks, for designing drug interventions in biomedical contexts and synthetic gene networks for bioengineering. These tools and the emphasis on behavior-shaping open new paths for efficiently exploring the complex behavior of biological networks. For the interactive version of this paper, please visit https://developmentalsystems.org/curious-exploration-of-grn-competencies.

Attractor-Transition Control of Complex Biological Networks: A Constant Control Approach

This article presents attractor-transition control of complex biological networks represented by Boolean networks (BNs) wherein the BN is steered from a prescribed initial attractor toward a desired one. The proposed approach leverages the similarity between attractors and Boolean algebraic properties embedded in the underlying state transition equations. To enhance the clarity of expression regarding stabilization toward the desired attractor, a simple coordinate transformation is performed on the considered BN. Based on the characteristics of transformed state equations, self-stabilizing state variables requiring no control efforts are derived in the first. Next, by applying the feedback vertex set (FVS) control scheme, control inputs stabilizing the remaining state variables are determined. The proposed control scheme exhibits versatility by accommodating both fixed-point and cyclic attractors. We validate the effectiveness of the proposed strategy through extensive numerical experiments conducted on random BNs as well as complex biological systems. In adherence to the reproducible research initiative, detailed results of numerical experiments and all the implementation codes are provided on the authors' website: https://github.com/choonlog/AttractorTransition.

A Boolean model explains phenotypic plasticity changes underlying hepatic cancer stem cells emergence

In several carcinomas, including hepatocellular carcinoma, it has been demonstrated that cancer stem cells (CSCs) have enhanced invasiveness and therapy resistance compared to differentiated cancer cells. Mathematical-computational tools could be valuable for integrating experimental results and understanding the phenotypic plasticity mechanisms for CSCs emergence. Based on the literature review, we constructed a Boolean model that recovers eight stable states (attractors) corresponding to the gene expression profile of hepatocytes and mesenchymal cells in senescent, quiescent, proliferative, and stem-like states. The epigenetic landscape associated with the regulatory network was analyzed. We observed that the loss of p53, p16, RB, or the constitutive activation of β-catenin and YAP1 increases the robustness of the proliferative stem-like phenotypes. Additionally, we found that p53 inactivation facilitates the transition of proliferative hepatocytes into stem-like mesenchymal phenotype. Thus, phenotypic plasticity may be altered, and stem-like phenotypes related to CSCs may be easier to attain following the mutation acquisition.

Comprehensive molecular interaction map of TGFβ induced epithelial to mesenchymal transition in breast cancer

Breast cancer is one of the prevailing cancers globally, with a high mortality rate. Metastatic breast cancer (MBC) is an advanced stage of cancer, characterised by a highly nonlinear, heterogeneous process involving numerous singling pathways and regulatory interactions. Epithelial-mesenchymal transition (EMT) emerges as a key mechanism exploited by cancer cells. Transforming Growth Factor-β (TGFβ)-dependent signalling is attributed to promote EMT in advanced stages of breast cancer. A comprehensive regulatory map of TGFβ induced EMT was developed through an extensive literature survey. The network assembled comprises of 312 distinct species (proteins, genes, RNAs, complexes), and 426 reactions (state transitions, nuclear translocations, complex associations, and dissociations). The map was developed by following Systems Biology Graphical Notation (SBGN) using Cell Designer and made publicly available using MINERVA ( http://35.174.227.105:8080/minerva/?id=Metastatic_Breast_Cancer_1 ). While the complete molecular mechanism of MBC is still not known, the map captures the elaborate signalling interplay of TGFβ induced EMT-promoting MBC. Subsequently, the disease map assembled was translated into a Boolean model utilising CaSQ and analysed using Cell Collective. Simulations of these have captured the known experimental outcomes of TGFβ induced EMT in MBC. Hub regulators of the assembled map were identified, and their transcriptome-based analysis confirmed their role in cancer metastasis. Elaborate analysis of this map may help in gaining additional insights into the development and progression of metastatic breast cancer.

Data- and theory-driven approaches for understanding paths of epithelial-mesenchymal transition

Reversible transitions between epithelial and mesenchymal cell states are a crucial form of epithelial plasticity for development and disease progression. Recent experimental data and mechanistic models showed multiple intermediate epithelial-mesenchymal transition (EMT) states as well as trajectories of EMT underpinned by complex gene regulatory networks. In this review, we summarize recent progress in quantifying EMT and characterizing EMT paths with computational methods and quantitative experiments including omics-level measurements. We provide perspectives on how these studies can help relating fundamental cell biology to physiological and pathological outcomes of EMT.

Proneural-mesenchymal antagonism dominates the patterns of phenotypic heterogeneity in glioblastoma

The aggressive nature of glioblastoma (GBM) - one of the deadliest forms of brain tumors - is majorly attributed to underlying phenotypic heterogeneity. Early attempts to classify this heterogeneity at a transcriptomic level in TCGA GBM cohort proposed the existence of four distinct molecular subtypes: Proneural, Neural, Classical, and Mesenchymal. Further, a single-cell RNA sequencing (scRNA-seq) analysis of primary tumors also reported similar four subtypes mimicking neurodevelopmental lineages. However, it remains unclear whether these four subtypes identified via bulk and single-cell transcriptomics are mutually exclusive or not. Here, we perform pairwise correlations among individual genes and gene signatures corresponding to these proposed subtypes and show that the subtypes are not distinctly mutually antagonistic in either TCGA or scRNA-seq data. We observed that the proneural (or neural progenitor-like)-mesenchymal axis is the most prominent antagonistic pair, with the other two subtypes lying on this spectrum. These results are reinforced through a meta-analysis of over 100 single-cell and bulk transcriptomic datasets as well as in terms of functional association with metabolic switching, cell cycle, and immune evasion pathways. Finally, this proneural-mesenchymal antagonistic trend percolates to the association of relevant transcription factors with patient survival. These results suggest rethinking GBM phenotypic characterization for more effective therapeutic targeting efforts.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics

Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.

Theoretical and computational tools to model multistable gene regulatory networks

The last decade has witnessed a surge of theoretical and computational models to describe the dynamics of complex gene regulatory networks, and how these interactions can give rise to multistable and heterogeneous cell populations. As the use of theoretical modeling to describe genetic and biochemical circuits becomes more widespread, theoreticians with mathematical and physical backgrounds routinely apply concepts from statistical physics, non-linear dynamics, and network theory to biological systems. This review aims at providing a clear overview of the most important methodologies applied in the field while highlighting current and future challenges. It also includes hands-on tutorials to solve and simulate some of the archetypical biological system models used in the field. Furthermore, we provide concrete examples from the existing literature for theoreticians that wish to explore this fast-developing field. Whenever possible, we highlight the similarities and differences between biochemical and regulatory networks and 'classical' systems typically studied in non-equilibrium statistical and quantum mechanics.

Theoretical and computational tools to model multistable gene regulatory networks

The last decade has witnessed a surge of theoretical and computational models to describe the dynamics of complex gene regulatory networks, and how these interactions can give rise to multistable and heterogeneous cell populations. As the use of theoretical modeling to describe genetic and biochemical circuits becomes more widespread, theoreticians with mathematical and physical backgrounds routinely apply concepts from statistical physics, non-linear dynamics, and network theory to biological systems. This review aims at providing a clear overview of the most important methodologies applied in the field while highlighting current and future challenges. It also includes hands-on tutorials to solve and simulate some of the archetypical biological system models used in the field. Furthermore, we provide concrete examples from the existing literature for theoreticians that wish to explore this fast-developing field. Whenever possible, we highlight the similarities and differences between biochemical and regulatory networks and 'classical' systems typically studied in non-equilibrium statistical and quantum mechanics.

Dynamics of Epithelial-Mesenchymal Plasticity: What Have Single-Cell Investigations Elucidated So Far?

Epithelial-mesenchymal plasticity (EMP) is a key driver of cancer metastasis and therapeutic resistance, through which cancer cells can reversibly and dynamically alter their molecular and functional traits along the epithelial-mesenchymal spectrum. While cells in the epithelial phenotype are usually tightly adherent, less metastatic, and drug-sensitive, those in the hybrid epithelial/mesenchymal and/or mesenchymal state are more invasive, migratory, drug-resistant, and immune-evasive. Single-cell studies have emerged as a powerful tool in gaining new insights into the dynamics of EMP across various cancer types. Here, we review many recent studies that employ single-cell analysis techniques to better understand the dynamics of EMP in cancer both in vitro and in vivo. These single-cell studies have underlined the plurality of trajectories cells can traverse during EMP and the consequent heterogeneity of hybrid epithelial/mesenchymal phenotypes seen at both preclinical and clinical levels. They also demonstrate how diverse EMP trajectories may exhibit hysteretic behavior and how the rate of such cell-state transitions depends on the genetic/epigenetic background of recipient cells, as well as the dose and/or duration of EMP-inducing growth factors. Finally, we discuss the relationship between EMP and patient survival across many cancer types. We also present a next set of questions related to EMP that could benefit much from single-cell observations and pave the way to better tackle phenotypic switching and heterogeneity in clinic.

Learning in Transcriptional Network Models: Computational Discovery of Pathway-Level Memory and Effective Interventions

Trainability, in any substrate, refers to the ability to change future behavior based on past experiences. An understanding of such capacity within biological cells and tissues would enable a particularly powerful set of methods for prediction and control of their behavior through specific patterns of stimuli. This top-down mode of control (as an alternative to bottom-up modification of hardware) has been extensively exploited by computer science and the behavioral sciences; in biology however, it is usually reserved for organism-level behavior in animals with brains, such as training animals towards a desired response. Exciting work in the field of basal cognition has begun to reveal degrees and forms of unconventional memory in non-neural tissues and even in subcellular biochemical dynamics. Here, we characterize biological gene regulatory circuit models and protein pathways and find them capable of several different kinds of memory. We extend prior results on learning in binary transcriptional networks to continuous models and identify specific interventions (regimes of stimulation, as opposed to network rewiring) that abolish undesirable network behavior such as drug pharmacoresistance and drug sensitization. We also explore the stability of created memories by assessing their long-term behavior and find that most memories do not decay over long time periods. Additionally, we find that the memory properties are quite robust to noise; surprisingly, in many cases noise actually increases memory potential. We examine various network properties associated with these behaviors and find that no one network property is indicative of memory. Random networks do not show similar memory behavior as models of biological processes, indicating that generic network dynamics are not solely responsible for trainability. Rational control of dynamic pathway function using stimuli derived from computational models opens the door to empirical studies of proto-cognitive capacities in unconventional embodiments and suggests numerous possible applications in biomedicine, where behavior shaping of pathway responses stand as a potential alternative to gene therapy.

"In medio stat virtus": Insights into hybrid E/M phenotype attitudes

Epithelial-mesenchymal plasticity (EMP) refers to the ability of cells to dynamically interconvert between epithelial (E) and mesenchymal (M) phenotypes, thus generating an array of hybrid E/M intermediates with mixed E and M features. Recent findings have demonstrated how these hybrid E/M rather than fully M cells play key roles in most of physiological and pathological processes involving EMT. To this regard, the onset of hybrid E/M state coincides with the highest stemness gene expression and is involved in differentiation of either normal and cancer stem cells. Moreover, hybrid E/M cells are responsible for wound healing and create a favorable immunosuppressive environment for tissue regeneration. Nevertheless, hybrid state is responsible of metastatic process and of the increasing of survival, apoptosis and therapy resistance in cancer cells. The present review aims to describe the main features and the emerging concepts regulating EMP and the formation of E/M hybrid intermediates by describing differences and similarities between cancer and normal hybrid stem cells. In particular, the comprehension of hybrid E/M cells biology will surely advance our understanding of their features and how they could be exploited to improve tissue regeneration and repair.

Network topology metrics explaining enrichment of hybrid epithelial mesenchymal phenotypes in metastasis

Epithelial to Mesenchymal Transition (EMT) and its reverse—Mesenchymal to Epithelial Transition (MET) are hallmarks of metastasis. Cancer cells use this reversible cellular programming to switch among Epithelial (E), Mesenchymal (M), and hybrid Epithelial/Mesenchymal (hybrid E/M) state(s) and seed tumors at distant sites. Hybrid E/M cells are often more aggressive and metastatic than the “pure” E and M cells. Thus, identifying mechanisms to inhibit hybrid E/M cells can be promising in curtailing metastasis. While multiple gene regulatory networks (GRNs) based mathematical models for EMT/MET have been developed recently, identifying topological signatures enriching hybrid E/M phenotypes remains to be done. Here, we investigate the dynamics of 13 different GRNs and report an interesting association between “hybridness” and the number of negative/positive feedback loops across the networks. While networks having more negative feedback loops favor hybrid phenotype(s), networks having more positive feedback loops (PFLs) or many HiLoops–specific combinations of PFLs, support terminal (E and M) phenotypes. We also establish a connection between “hybridness” and network-frustration by showing that hybrid phenotypes likely result from non-reinforcing interactions among network nodes (genes) and therefore tend to be more frustrated (less stable). Our analysis, thus, identifies network topology-based signatures that can give rise to, as well as prevent, the emergence of hybrid E/M phenotype in GRNs underlying EMP. Our results can have implications in terms of targeting specific interactions in GRNs as a potent way to restrict switching to the hybrid E/M phenotype(s) to curtail metastasis.

Landscape of epithelial-mesenchymal plasticity as an emergent property of coordinated teams in regulatory networks

Elucidating the design principles of regulatory networks driving cellular decision-making has fundamental implications in mapping and eventually controlling cell-fate decisions. Despite being complex, these regulatory networks often only give rise to a few phenotypes. Previously, we identified two 'teams' of nodes in a small cell lung cancer regulatory network that constrained the phenotypic repertoire and aligned strongly with the dominant phenotypes obtained from network simulations (Chauhan et al., 2021). However, it remained elusive whether these 'teams' exist in other networks, and how do they shape the phenotypic landscape. Here, we demonstrate that five different networks of varying sizes governing epithelial-mesenchymal plasticity comprised of two 'teams' of players - one comprised of canonical drivers of epithelial phenotype and the other containing the mesenchymal inducers. These 'teams' are specific to the topology of these regulatory networks and orchestrate a bimodal phenotypic landscape with the epithelial and mesenchymal phenotypes being more frequent and dynamically robust to perturbations, relative to the intermediary/hybrid epithelial/mesenchymal ones. Our analysis reveals that network topology alone can contain information about corresponding phenotypic distributions, thus obviating the need to simulate them. We propose 'teams' of nodes as a network design principle that can drive cell-fate canalization in diverse decision-making processes.

Cooperative RNA degradation stabilizes intermediate epithelial-mesenchymal states and supports a phenotypic continuum

Multiple intermediate epithelial-mesenchymal transition (EMT) states reflecting hybrid epithelial and mesenchymal phenotypes were observed in physiological and pathological conditions. Previous theoretical models explaining multiple EMT states rely on regulatory loops involving transcriptional feedback, which produce three or four attractors. This is incompatible with the observed continuum-like EMT spectrum. Here, we used mass-action-based models to describe post-transcriptional regulations, finding that cooperative RNA degradation via multiple microRNA binding sites can generate four-attractor systems without transcriptional feedback. Furthermore, the newly identified intermediates-enabling circuits are common in the EMT regulatory network, and they can synergize with transcriptional feedback to support phenotypic continuum. Finally, our model predicted a role of miR-101 in multistate EMT, and we identified evidence from single-cell RNA-sequencing data that support the prediction. Our work reveals a previously unknown role of cooperative RNA degradation and microRNAs in EMT, providing a framework that can bridge the gap between mechanistic models and single-cell experiments.

Probabilistic edge weights fine-tune Boolean network dynamics

Biological systems are noisy by nature. This aspect is reflected in our experimental measurements and should be reflected in the models we build to better understand these systems. Noise can be especially consequential when trying to interpret specific regulatory interactions, i.e. regulatory network edges. In this paper, we propose a method to explicitly encode edge-noise in Boolean dynamical systems by probabilistic edge-weight (PEW) operators. PEW operators have two important features: first, they introduce a form of edge-weight into Boolean models through the noise, second, the noise is dependent on the dynamical state of the system, which enables more biologically meaningful modeling choices. Moreover, we offer a simple-to-use implementation in the already well-established BooleanNet framework. In two application cases, we show how the introduction of just a few PEW operators in Boolean models can fine-tune the emergent dynamics and increase the accuracy of qualitative predictions. This includes fine-tuning interactions which cause non-biological behaviors when switching between asynchronous and synchronous update schemes in dynamical simulations. Moreover, PEW operators also open the way to encode more exotic cellular dynamics, such as cellular learning, and to implementing edge-weights for regulatory networks inferred from omics data.

SNAIL driven by a feed forward loop motif promotes TGFβinduced epithelial to mesenchymal transition

Epithelial to Mesenchymal Transition (EMT) plays an important role in tissue regeneration, embryonic development, and cancer metastasis. Several signaling pathways are known to regulate EMT, among which the modulation of TGFβ(Transforming Growth Factor-β) induced EMT is crucial in several cancer types. Several mathematical models were built to explore the role of core regulatory circuit of ZEB/miR-200, SNAIL/miR-34 double negative feedback loops in modulating TGFβinduced EMT. Different emergent behavior including tristability, irreversible switching, existence of hybrid EMT states were inferred though these models. Some studies have explored the role of TGFβreceptor activation, SMADs nucleocytoplasmic shuttling and complex formation. Recent experiments have revealed that MDM2 along with SMAD complex regulates SNAIL expression driven EMT. Encouraged by this, in the present study we developed a mathematical model for p53/MDM2 dependent TGFβinduced EMT regulation. Inclusion of p53 brings in an additional mechanistic perspective in exploring the EM transition. The network formulated comprises a C1FFL moderating SNAIL expression involving MDM2 and SMAD complex, which functions as a noise filter and persistent detector. The C1FFL was also observed to operate as a coincidence detector driving the SNAIL dependent downstream signaling into phenotypic switching decision. Systems modelling and analysis of the devised network, displayed interesting dynamic behavior, systems response to various inputs stimulus, providing a better understanding of p53/MDM2 dependent TGF-βinduced Epithelial to Mesenchymal Transition.

The Wnt pathway can stabilize hybrid phenotypes in the epithelial-mesenchymal transition: A logical modeling approach

The Wnt pathway is important to regulate a variety of biochemical functions and can contribute to cancer development through its influence on the epithelial-mesenchymal transition (EMT). Multiple circuits have been reported to participate in the regulation of the Wnt signaling, however, the way these circuits coordinately regulate this signaling is still unclear. Moreover, the mechanisms responsible for the appearance of hybrid phenotypes (cells presenting both E and M features) are not well determined. The hybrid phenotype can present much higher metastatic potential than the mesenchymal phenotype. In this study, we propose a Boolean model of the Wnt pathway signaling contemplating recent published biochemical information on hepatocarcinoma. The model presents good coherence with experimental data for perturbed and wild-type cases. With the model, we propose two new molecular circuits involving several molecules that can stabilize hybrid states during the EMT. Moreover, we found that the two well studied circuits, AKT1/β-catenin and SNAIL1/miR-34, can cooperate with the predicted ones to favor the stabilization of the hybrid states. These findings highlight some possible unrecognized mechanisms during Wnt signaling and may provide alternative therapeutic strategies to control cancer metastatization.

Translocating proteins compartment-specifically alter the fate of epithelial-mesenchymal transition in a compartmentalized Boolean network model

Regulation of translocating proteins is crucial in defining cellular behaviour. Epithelial-mesenchymal transition (EMT) is important in cellular processes, such as cancer progression. Several orchestrators of EMT, such as key transcription factors, are known to translocate. We show that translocating proteins become enriched in EMT-signalling. To simulate the compartment-specific functions of translocating proteins we created a compartmentalized Boolean network model. This model successfully reproduced known biological traits of EMT and as a novel feature it also captured organelle-specific functions of proteins. Our results predicted that glycogen synthase kinase-3 beta (GSK3B) compartment-specifically alters the fate of EMT, amongst others the activation of nuclear GSK3B halts transforming growth factor beta-1 (TGFB) induced EMT. Moreover, our results recapitulated that the nuclear activation of glioma associated oncogene transcription factors (GLI) is needed to achieve a complete EMT. Compartmentalized network models will be useful to uncover novel control mechanisms of biological processes. Our algorithmic procedures can be automatically rerun on the https://translocaboole.linkgroup.hu website, which provides a framework for similar future studies.

Heterogeneity and dynamic of EMT through the plasticity of ribosome and mRNA translation

Growing evidence exposes translation and its translational machinery as key players in establishing and maintaining physiological and pathological biological processes. Examining translation may not only provide new biological insight but also identify novel innovative therapeutic targets in several fields of biology, including that of epithelial-to-mesenchymal transition (EMT). EMT is currently considered as a dynamic and reversible transdifferentiation process sustaining the transition from an epithelial to mesenchymal phenotype, known to be mainly driven by transcriptional reprogramming. However, it seems that the characterization of EMT plasticity is challenging, relying exclusively on transcriptomic and epigenetic approaches. Indeed, heterogeneity in EMT programs was reported to depend on the biological context. Here, by reviewing the involvement of translational control, translational machinery and ribosome biogenesis characterizing the different types of EMT, from embryonic and adult physiological to pathological contexts, we discuss the added value of integrating translational control and its machinery to depict the heterogeneity and dynamics of EMT programs.

Population Dynamics of Epithelial-Mesenchymal Heterogeneity in Cancer Cells

Phenotypic heterogeneity is a hallmark of aggressive cancer behaviour and a clinical challenge. Despite much characterisation of this heterogeneity at a multi-omics level in many cancers, we have a limited understanding of how this heterogeneity emerges spontaneously in an isogenic cell population. Some longitudinal observations of dynamics in epithelial-mesenchymal heterogeneity, a canonical example of phenotypic heterogeneity, have offered us opportunities to quantify the rates of phenotypic switching that may drive such heterogeneity. Here, we offer a mathematical modeling framework that explains the salient features of population dynamics noted in PMC42-LA cells: (a) predominance of EpCAM^high subpopulation, (b) re-establishment of parental distributions from the EpCAM^high and EpCAM^low subpopulations, and (c) enhanced heterogeneity in clonal populations established from individual cells. Our framework proposes that fluctuations or noise in content duplication and partitioning of SNAIL—an EMT-inducing transcription factor—during cell division can explain spontaneous phenotypic switching and consequent dynamic heterogeneity in PMC42-LA cells observed experimentally at both single-cell and bulk level analysis. Together, we propose that asymmetric cell division can be a potential mechanism for phenotypic heterogeneity.

Capturing Biomarkers and Molecular Targets in Cellular Landscapes From Dynamic Reaction Network Models and Machine Learning

Computational dynamic ODE models of cell function describing biochemical reactions have been created for decades, but on a small scale. Still, they have been highly effective in describing and predicting behaviors. For example, oscillatory phospho-ERK levels were predicted and confirmed in MAPK signaling encompassing both positive and negative feedback loops. These models typically were limited and not adapted to large datasets so commonly found today. But importantly, ODE models describe reaction networks in well-mixed systems representing the cell and can be simulated with ordinary differential equations that are solved deterministically. Stochastic solutions, which can account for noisy reaction networks, in some cases, also improve predictions. Today, dynamic ODE models rarely encompass an entire cell even though it might be expected that an upload of the large genomic, transcriptomic, and proteomic datasets may allow whole cell models. It is proposed here to combine output from simulated dynamic ODE models, completed with omics data, to discover both biomarkers in cancer a priori and molecular targets in the Machine Learning setting.

An Integrated View of Virus-Triggered Cellular Plasticity Using Boolean Networks

Virus-related mortality and morbidity are due to cell/tissue damage caused by replicative pressure and resource exhaustion, e.g., HBV or HIV; exaggerated immune responses, e.g., SARS-CoV-2; and cancer, e.g., EBV or HPV. In this context, oncogenic and other types of viruses drive genetic and epigenetic changes that expand the tumorigenic program, including modifications to the ability of cancer cells to migrate. The best-characterized group of changes is collectively known as the epithelial–mesenchymal transition, or EMT. This is a complex phenomenon classically described using biochemistry, cell biology and genetics. However, these methods require enormous, often slow, efforts to identify and validate novel therapeutic targets. Systems biology can complement and accelerate discoveries in this field. One example of such an approach is Boolean networks, which make complex biological problems tractable by modeling data (“nodes”) connected by logical operators. Here, we focus on virus-induced cellular plasticity and cell reprogramming in mammals, and how Boolean networks could provide novel insights into the ability of some viruses to trigger uncontrolled cell proliferation and EMT, two key hallmarks of cancer.

Landscape and kinetic path quantify critical transitions in epithelial-mesenchymal transition

Epithelial-mesenchymal transition (EMT), a basic developmental process that might promote cancer metastasis, has been studied from various perspectives. Recently, the early warning theory has been used to anticipate critical transitions in EMT from mathematical modeling. However, the underlying mechanisms of EMT involving complex molecular networks remain to be clarified. Especially, how to quantify the global stability and stochastic transition dynamics of EMT and what the underlying mechanism for early warning theory in EMT is remain to be fully clarified. To address these issues, we constructed a comprehensive gene regulatory network model for EMT and quantified the corresponding potential landscape. The landscape for EMT displays multiple stable attractors, which correspond to E, M, and some other intermediate states. Based on the path-integral approach, we identified the most probable transition paths of EMT, which are supported by experimental data. Correspondingly, the results of transition actions demonstrated that intermediate states can accelerate EMT, consistent with recent studies. By integrating the landscape and path with early warning concept, we identified the potential barrier height from the landscape as a global and more accurate measure for early warning signals to predict critical transitions in EMT. The landscape results also provide an intuitive and quantitative explanation for the early warning theory. Overall, the landscape and path results advance our mechanistic understanding of dynamical transitions and roles of intermediate states in EMT, and the potential barrier height provides a new, to our knowledge, measure for critical transitions and quantitative explanations for the early warning theory.

EP300/CBP is crucial for cAMP-PKA pathway to alleviate podocyte dedifferentiation via targeting Notch3 signaling

Podocyte injury is the hallmark of proteinuric glomerular diseases. Notch3 is neo-activated simultaneously in damaged podocytes and podocyte's progenitor cells of FSGS, indicating a unique role of Notch3. We previously showed that activation of cAMP-PKA pathway alleviated podocyte injury possibly via inhibiting Notch3 expression. However, the mechanisms are unknown. In the present study, Notch3 signaling was significantly activated in ADR-induced podocytes in vitro and in PAN nephrosis rats and patients with idiopathic FSGS in vivo, concomitantly with podocyte dedifferentiation. In cultured podocytes, pCPT-cAMP, a selective cAMP-PKA activator, dramatically blocked ADR-induced activation of Notch3 signaling as well as inhibition of cAMP-PKA pathway, thus alleviating the decreased cell viability and podocyte dedifferentiation. Bioinformatics analysis revealed EP300/CBP, a transcriptional co-activator, as a central hub for the crosstalk between these two signaling pathways. Additionally, CREB/KLF15 in cAMP-PKA pathway competed with RBP-J the major transcriptional factor of Notch3 signaling for binding to EP300/CBP. EP300/CBP siRNA significantly inhibited these two signaling transduction pathways and disrupted the interactions between the above major transcriptional factors. These data indicate a crucial role of EP300/CBP in regulating the crosstalk between cAMP-PKA pathway and Notch3 signaling and modulating the phenotypic change of podocytes, and enrich the reno-protective mechanisms of cAMP-PKA pathway.

A Personalized Therapeutics Approach Using an In Silico Drosophila Patient Model Reveals Optimal Chemo- and Targeted Therapy Combinations for Colorectal Cancer

In silico models of biomolecular regulation in cancer, annotated with patient-specific gene expression data, can aid in the development of novel personalized cancer therapeutic strategies. Drosophila melanogaster is a well-established animal model that is increasingly being employed to evaluate such preclinical personalized cancer therapies. Here, we report five Boolean network models of biomolecular regulation in cells lining the Drosophila midgut epithelium and annotate them with colorectal cancer patient-specific mutation data to develop an in silico Drosophila Patient Model (DPM). We employed cell-type-specific RNA-seq gene expression data from the FlyGut-seq database to annotate and then validate these networks. Next, we developed three literature-based colorectal cancer case studies to evaluate cell fate outcomes from the model. Results obtained from analyses of the proposed DPM help: (i) elucidate cell fate evolution in colorectal tumorigenesis, (ii) validate cytotoxicity of nine FDA-approved CRC drugs, and (iii) devise optimal personalized treatment combinations. The personalized network models helped identify synergistic combinations of paclitaxel-regorafenib, paclitaxel-bortezomib, docetaxel-bortezomib, and paclitaxel-imatinib for treating different colorectal cancer patients. Follow-on therapeutic screening of six colorectal cancer patients from cBioPortal using this drug combination demonstrated a 100% increase in apoptosis and a 100% decrease in proliferation. In conclusion, this work outlines a novel roadmap for decoding colorectal tumorigenesis along with the development of personalized combinatorial therapeutics for preclinical translational studies.

Cell Line-Specific Network Models of ER+ Breast Cancer Identify Potential PI3Kα Inhibitor Resistance Mechanisms and Drug Combinations

Durable control of invasive solid tumors necessitates identifying therapeutic resistance mechanisms and effective drug combinations. In this work, we used a network-based mathematical model to identify sensitivity regulators and drug combinations for the PI3Kα inhibitor alpelisib in estrogen receptor positive (ER⁺) PIK3CA-mutant breast cancer. The model-predicted efficacious combination of alpelisib and BH3 mimetics, for example, MCL1 inhibitors, was experimentally validated in ER⁺ breast cancer cell lines. Consistent with the model, FOXO3 downregulation reduced sensitivity to alpelisib, revealing a novel potential resistance mechanism. Cell line-specific sensitivity to combinations of alpelisib and BH3 mimetics depended on which BCL2 family members were highly expressed. On the basis of these results, newly developed cell line-specific network models were able to recapitulate the observed differential response to alpelisib and BH3 mimetics. This approach illustrates how network-based mathematical models can contribute to overcoming the challenge of cancer drug resistance. SIGNIFICANCE: Network-based mathematical models of oncogenic signaling and experimental validation of its predictions can identify resistance mechanisms for targeted therapies, as this study demonstrates for PI3Kα-specific inhibitors in breast cancer.

Parity and time reversal elucidate both decision-making in empirical models and attractor scaling in critical Boolean networks

We present new applications of parity inversion and time reversal to the emergence of complex behavior from simple dynamical rules in stochastic discrete models. Our parity-based encoding of causal relationships and time-reversal construction efficiently reveal discrete analogs of stable and unstable manifolds. We demonstrate their predictive power by studying decision-making in systems biology and statistical physics models. These applications underpin a novel attractor identification algorithm implemented for Boolean networks under stochastic dynamics. Its speed enables resolving a long-standing open question of how attractor count in critical random Boolean networks scales with network size and whether the scaling matches biological observations. Via 80-fold improvement in probed network size (N = 16,384), we find the unexpectedly low scaling exponent of 0.12 ± 0.05, approximately one-tenth the analytical upper bound. We demonstrate a general principle: A system's relationship to its time reversal and state-space inversion constrains its repertoire of emergent behaviors.

Computational systems-biology approaches for modeling gene networks driving epithelial-mesenchymal transitions

Epithelial-mesenchymal transition (EMT) is an important biological process through which epithelial cells undergo phenotypic transitions to mesenchymal cells by losing cell-cell adhesion and gaining migratory properties that cells use in embryogenesis, wound healing, and cancer metastasis. An important research topic is to identify the underlying gene regulatory networks (GRNs) governing the decision making of EMT and develop predictive models based on the GRNs. The advent of recent genomic technology, such as single-cell RNA sequencing, has opened new opportunities to improve our understanding about the dynamical controls of EMT. In this article, we review three major types of computational and mathematical approaches and methods for inferring and modeling GRNs driving EMT. We emphasize (1) the bottom-up approaches, where GRNs are constructed through literature search; (2) the top-down approaches, where GRNs are derived from genome-wide sequencing data; (3) the combined top-down and bottom-up approaches, where EMT GRNs are constructed and simulated by integrating bioinformatics and mathematical modeling. We discuss the methodologies and applications of each approach and the available resources for these studies.

Classification of triple-negative breast cancers through a Boolean network model of the epithelial-mesenchymal transition

Predicting the metastasis risk in patients with a primary breast cancer tumor is of fundamental importance to decide the best therapeutic strategy in the framework of personalized medicine. Here, we present ARIADNE, a general algorithmic strategy to assess the risk of metastasis from transcriptomic data of patients with triple-negative breast cancer, a subtype of breast cancer with poorer prognosis with respect to the other subtypes. ARIADNE identifies hybrid epithelial/mesenchymal phenotypes by mapping gene expression data into the states of a Boolean network model of the epithelial-mesenchymal pathway. Using this mapping, it is possible to stratify patients according to their prognosis, as we show by validating the strategy with three independent cohorts of triple-negative breast cancer patients. Our strategy provides a prognostic tool that could be applied to other biologically relevant pathways, in order to estimate the metastatic risk for other breast cancer subtypes or other tumor types. A record of this paper's transparent peer review process is included in the supplemental information.

In Silico Logical Modelling to Uncover Cooperative Interactions in Cancer

The multistep development of cancer involves the cooperation between multiple molecular lesions, as well as complex interactions between cancer cells and the surrounding tumour microenvironment. The search for these synergistic interactions using experimental models made tremendous contributions to our understanding of oncogenesis. Yet, these approaches remain labour-intensive and challenging. To tackle such a hurdle, an integrative, multidisciplinary effort is required. In this article, we highlight the use of logical computational models, combined with experimental validations, as an effective approach to identify cooperative mechanisms and therapeutic strategies in the context of cancer biology. In silico models overcome limitations of reductionist approaches by capturing tumour complexity and by generating powerful testable hypotheses. We review representative examples of logical models reported in the literature and their validation. We then provide further analyses of our logical model of Epithelium to Mesenchymal Transition (EMT), searching for additional cooperative interactions involving inputs from the tumour microenvironment and gain of function mutations in NOTCH.

Measuring and Modelling the Epithelial- Mesenchymal Hybrid State in Cancer: Clinical Implications

The epithelial-mesenchymal (E/M) hybrid state has emerged as an important mediator of elements of cancer progression, facilitated by epithelial mesenchymal plasticity (EMP). We review here evidence for the presence, prognostic significance, and therapeutic potential of the E/M hybrid state in carcinoma. We further assess modelling predictions and validation studies to demonstrate stabilised E/M hybrid states along the spectrum of EMP, as well as computational approaches for characterising and quantifying EMP phenotypes, with particular attention to the emerging realm of single-cell approaches through RNA sequencing and protein-based techniques.

Mathematical modeling of the Candida albicans yeast to hyphal transition reveals novel control strategies

Candida albicans, an opportunistic fungal pathogen, is a significant cause of human infections, particularly in immunocompromised individuals. Phenotypic plasticity between two morphological phenotypes, yeast and hyphae, is a key mechanism by which C. albicans can thrive in many microenvironments and cause disease in the host. Understanding the decision points and key driver genes controlling this important transition and how these genes respond to different environmental signals is critical to understanding how C. albicans causes infections in the host. Here we build and analyze a Boolean dynamical model of the C. albicans yeast to hyphal transition, integrating multiple environmental factors and regulatory mechanisms. We validate the model by a systematic comparison to prior experiments, which led to agreement in 17 out of 22 cases. The discrepancies motivate alternative hypotheses that are testable by follow-up experiments. Analysis of this model revealed two time-constrained windows of opportunity that must be met for the complete transition from the yeast to hyphal phenotype, as well as control strategies that can robustly prevent this transition. We experimentally validate two of these control predictions in C. albicans strains lacking the transcription factor UME6 and the histone deacetylase HDA1, respectively. This model will serve as a strong base from which to develop a systems biology understanding of C. albicans morphogenesis.

Multi-Stability and Consequent Phenotypic Plasticity in AMPK-Akt Double Negative Feedback Loop in Cancer Cells

Adaptation and survival of cancer cells to various stress and growth factor conditions is crucial for successful metastasis. A double-negative feedback loop between two serine/threonine kinases AMPK (AMP-activated protein kinase) and Akt can regulate the adaptation of breast cancer cells to matrix-deprivation stress. This feedback loop can significantly generate two phenotypes or cell states: matrix detachment-triggered pAMPK^high/ pAkt^low state, and matrix (re)attachment-triggered pAkt^high/ pAMPK^low state. However, whether these two cell states can exhibit phenotypic plasticity and heterogeneity in a given cell population, i.e., whether they can co-exist and undergo spontaneous switching to generate the other subpopulation, remains unclear. Here, we develop a mechanism-based mathematical model that captures the set of experimentally reported interactions among AMPK and Akt. Our simulations suggest that the AMPK-Akt feedback loop can give rise to two co-existing phenotypes (pAkt^high/ pAMPK^low and pAMPK^high/pAkt^low) in specific parameter regimes. Next, to test the model predictions, we segregated these two subpopulations in MDA-MB-231 cells and observed that each of them was capable of switching to another in adherent conditions. Finally, the predicted trends are supported by clinical data analysis of The Cancer Genome Atlas (TCGA) breast cancer and pan-cancer cohorts that revealed negatively correlated pAMPK and pAkt protein levels. Overall, our integrated computational-experimental approach unravels that AMPK-Akt feedback loop can generate multi-stability and drive phenotypic switching and heterogeneity in a cancer cell population.

Context Matters: NOTCH Signatures and Pathway in Cancer Progression and Metastasis

The Notch signaling pathway is a critical player in embryogenesis but also plays various roles in tumorigenesis, with both tumor suppressor and oncogenic activities. Mutations, deletions, amplifications, or over-expression of Notch receptors, ligands, and a growing list of downstream Notch-activated genes have by now been described for most human cancer types. Yet, it often remains unclear what may be the functional impact of these changes for tumor biology, initiation, and progression, for cancer therapy, and for personalized medicine. Emerging data indicate that Notch signaling can also contribute to increased aggressive properties such as invasion, tumor heterogeneity, angiogenesis, or tumor cell dormancy within solid cancer tissues; especially in epithelial cancers, which are in the center of this review. Notch further supports the “stemness” of cancer cells and helps define the stem cell niche for their long-term survival, by integrating the interaction between cancer cells and the cells of the tumor microenvironment (TME). The complexity of Notch crosstalk with other signaling pathways and its roles in cell fate and trans-differentiation processes such as epithelial-to-mesenchymal transition (EMT) point to this pathway as a decisive player that may tip the balance between tumor suppression and promotion, differentiation and invasion. Here we not only review the literature, but also explore genomic databases with a specific focus on Notch signatures, and how they relate to different stages in tumor development. Altered Notch signaling hereby plays a key role for tumor cell survival and coping with a broad spectrum of vital issues, contributing to failed therapies, poor patient outcome, and loss of lives.

Hybrid E/M Phenotype(s) and Stemness: A Mechanistic Connection Embedded in Network Topology

Metastasis remains an unsolved clinical challenge. Two crucial features of metastasizing cancer cells are (a) their ability to dynamically move along the epithelial-hybrid-mesenchymal spectrum and (b) their tumor initiation potential or stemness. With increasing functional characterization of hybrid epithelial/mesenchymal (E/M) phenotypes along the spectrum, recent in vitro and in vivo studies have suggested an increasing association of hybrid E/M phenotypes with stemness. However, the mechanistic underpinnings enabling this association remain unclear. Here, we develop a mechanism-based mathematical modeling framework that interrogates the emergent nonlinear dynamics of the coupled network modules regulating E/M plasticity (miR-200/ZEB) and stemness (LIN28/let-7). Simulating the dynamics of this coupled network across a large ensemble of parameter sets, we observe that hybrid E/M phenotype(s) are more likely to acquire stemness relative to "pure" epithelial or mesenchymal states. We also integrate multiple "phenotypic stability factors" (PSFs) that have been shown to stabilize hybrid E/M phenotypes both in silico and in vitro-such as OVOL1/2, GRHL2, and NRF2-with this network, and demonstrate that the enrichment of hybrid E/M phenotype(s) with stemness is largely conserved in the presence of these PSFs. Thus, our results offer mechanistic insights into recent experimental observations of hybrid E/M phenotype(s) that are essential for tumor initiation and highlight how this feature is embedded in the underlying topology of interconnected EMT (Epithelial-Mesenchymal Transition) and stemness networks.

Hypoxia, partial EMT and collective migration: Emerging culprits in metastasis

Epithelial-mesenchymal transition (EMT) is a cellular biological process involved in migration of primary cancer cells to secondary sites facilitating metastasis. Besides, EMT also confers properties such as stemness, drug resistance and immune evasion which can aid a successful colonization at the distant site. EMT is not a binary process; recent evidence suggests that cells in partial EMT or hybrid E/M phenotype(s) can have enhanced stemness and drug resistance as compared to those undergoing a complete EMT. Moreover, partial EMT enables collective migration of cells as clusters of circulating tumor cells or emboli, further endorsing that cells in hybrid E/M phenotypes may be the 'fittest' for metastasis. Here, we review mechanisms and implications of hybrid E/M phenotypes, including their reported association with hypoxia. Hypoxia-driven activation of HIF-1α can drive EMT. In addition, cyclic hypoxia, as compared to acute or chronic hypoxia, shows the highest levels of active HIF-1α and can augment cancer aggressiveness to a greater extent, including enriching for a partial EMT phenotype. We also discuss how metastasis is influenced by hypoxia, partial EMT and collective cell migration, and call for a better understanding of interconnections among these mechanisms. We discuss the known regulators of hypoxia, hybrid EMT and collective cell migration and highlight the gaps which needs to be filled for connecting these three axes which will increase our understanding of dynamics of metastasis and help control it more effectively.

Simulation of multiple microenvironments shows a pivot role of RPTPs on the control of Epithelial-to-Mesenchymal Transition

Epithelial-to-Mesenchymal Transition (EMT) is a natural and reversible process involved in embryogenesis, wound healing and thought to participate in the process of metastasis. Multiple signals from the microenvironment have been reported to drive EMT. However, the tight control of this process on physiological scenarios and how it is disrupted during cancer progression is not fully understood. Here, we analysed a regulatory network of EMT accounting for 10 key microenvironment signals focusing on the impact of two cell contact signals on the reversibility of EMT and the stability of resulting phenotypes. The analysis showed that the microenvironment is not enough for stabilizing Hybrid and Amoeboid-like phenotypes, requiring intracellular de-regulations as reported during cancer progression. Our simulations demonstrated that RPTP activation by cell contacts have the potential to inhibit the process of EMT and trigger its reversibility under tissue growth and chronic inflammation scenarios. Simulations also showed that hypoxia inhibits the capacity of RPTPs to control EMT. Our analysis further provided a theoretical explanation for the observed correlation between hypoxia and metastasis under chronic inflammation, and predicted that de-regulations in FAT4 signalling may promote Hybrid stabilization. Taken together, we propose a natural control mechanism of EMT that supports the idea that EMT is tightly regulated by the microenvironment.