Towards control of cellular decision-making networks in the epithelial-to-mesenchymal transition

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.

Modeling stress-induced responses: plasticity in continuous state space and gradual clonal evolution

Mathematical models of cancer and bacterial evolution have generally stemmed from a gene-centric framework, assuming clonal evolution via acquisition of resistance-conferring mutations and selection of their corresponding subpopulations. More recently, the role of phenotypic plasticity has been recognized and models accounting for phenotypic switching between discrete cell states (e.g., epithelial and mesenchymal) have been developed. However, seldom do models incorporate both plasticity and mutationally driven resistance, particularly when the state space is continuous and resistance evolves in a continuous fashion. In this paper, we develop a framework to model plastic and mutational mechanisms of acquiring resistance in a continuous gradual fashion. We use this framework to examine ways in which cancer and bacterial populations can respond to stress and consider implications for therapeutic strategies. Although we primarily discuss our framework in the context of cancer and bacteria, it applies broadly to any system capable of evolving via plasticity and genetic evolution.

MaxCal can infer models from coupled stochastic trajectories of gene expression and cell division

Gene expression is inherently noisy due to small numbers of proteins and nucleic acids inside a cell. Likewise, cell division is stochastic, particularly when tracking at the level of a single cell. The two can be coupled when gene expression affects the rate of cell division. Single-cell time-lapse experiments can measure both fluctuations by simultaneously recording protein levels inside a cell and its stochastic division. These information-rich noisy trajectory data sets can be harnessed to learn about the underlying molecular and cellular details that are often not known a priori. A critical question is: How can we infer a model given data where fluctuations at two levels-gene expression and cell division-are intricately convoluted? We show the principle of maximum caliber (MaxCal)-integrated within a Bayesian framework-can be used to infer several cellular and molecular details (division rates, protein production, and degradation rates) from these coupled stochastic trajectories (CSTs). We demonstrate this proof of concept using synthetic data generated from a known model. An additional challenge in data analysis is that trajectories are often not in protein numbers, but in noisy fluorescence that depends on protein number in a probabilistic manner. We again show that MaxCal can infer important molecular and cellular rates even when data are in fluorescence, another example of CST with three confounding factors-gene expression noise, cell division noise, and fluorescence distortion-all coupled. Our approach will provide guidance to build models in synthetic biology experiments as well as general biological systems where examples of CSTs are abundant.

Benchmarking of protein interaction databases for integration with manually reconstructed signalling network models

Protein interaction databases are critical resources for network bioinformatics and integrating molecular experimental data. Interaction databases may also enable construction of predictive computational models of biological networks, although their fidelity for this purpose is not clear. Here, we benchmark protein interaction databases X2K, Reactome, Pathway Commons, Omnipath and Signor for their ability to recover manually curated edges from three logic-based network models of cardiac hypertrophy, mechano-signalling and fibrosis. Pathway Commons performed best at recovering interactions from manually reconstructed hypertrophy (137 of 193 interactions, 71%), mechano-signalling (85 of 125 interactions, 68%) and fibroblast networks (98 of 142 interactions, 69%). While protein interaction databases successfully recovered central, well-conserved pathways, they performed worse at recovering tissue-specific and transcriptional regulation. This highlights a knowledge gap where manual curation is critical. Finally, we tested the ability of Signor and Pathway Commons to identify new edges that improve model predictions, revealing important roles of protein kinase C autophosphorylation and Ca2+ /calmodulin-dependent protein kinase II phosphorylation of CREB in cardiomyocyte hypertrophy. This study provides a platform for benchmarking protein interaction databases for their utility in network model construction, as well as providing new insights into cardiac hypertrophy signalling. KEY POINTS: Protein interaction databases are used to recover signalling interactions from previously developed network models. The five protein interaction databases benchmarked recovered well-conserved pathways, but did poorly at recovering tissue-specific pathways and transcriptional regulation, indicating the importance of manual curation. We identify new signalling interactions not previously used in the network models, including a role for Ca2+ /calmodulin-dependent protein kinase II phosphorylation of CREB in cardiomyocyte hypertrophy.

An in vitro culture platform for studying the effect of collective cell migration on spatial self-organization within induced pluripotent stem cell colonies

BACKGROUND

Human induced pluripotent stem cells (hiPSCs) provide an in vitro system to identify the impact of cell behavior on the earliest stages of cell fate specification during human development. Here, we developed an hiPSC-based model to study the effect of collective cell migration in meso-endodermal lineage segregation and cell fate decisions through the control of space confinement using a detachable ring culture system.

RESULTS

The actomyosin organization of cells at the edge of undifferentiated colonies formed in a ring barrier differed from that of the cells in the center of the colony. In addition, even in the absence of exogenous supplements, ectoderm, mesoderm, endoderm, and extraembryonic cells differentiated following the induction of collective cell migration at the colony edge by removing the ring-barrier. However, when collective cell migration was inhibited by blocking E-cadherin function, this fate decision within an hiPSC colony was altered to an ectodermal fate. Furthermore, the induction of collective cell migration at the colony edge using an endodermal induction media enhanced endodermal differentiation efficiency in association with cadherin switching, which is involved in the epithelial-mesenchymal transition.

CONCLUSIONS

Our findings suggest that collective cell migration can be an effective way to drive the segregation of mesoderm and endoderm lineages, and cell fate decisions of hiPSCs.

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

NETISCE: a network-based tool for cell fate reprogramming

The search for effective therapeutic targets in fields like regenerative medicine and cancer research has generated interest in cell fate reprogramming. This cellular reprogramming paradigm can drive cells to a desired target state from any initial state. However, methods for identifying reprogramming targets remain limited for biological systems that lack large sets of experimental data or a dynamical characterization. We present NETISCE, a novel computational tool for identifying cell fate reprogramming targets in static networks. In combination with machine learning algorithms, NETISCE estimates the attractor landscape and predicts reprogramming targets using signal flow analysis and feedback vertex set control, respectively. Through validations in studies of cell fate reprogramming from developmental, stem cell, and cancer biology, we show that NETISCE can predict previously identified cell fate reprogramming targets and identify potentially novel combinations of targets. NETISCE extends cell fate reprogramming studies to larger-scale biological networks without the need for full model parameterization and can be implemented by experimental and computational biologists to identify parts of a biological system relevant to the desired reprogramming task.

Structure-based approach to identifying small sets of driver nodes in biological networks

In network control theory, driving all the nodes in the Feedback Vertex Set (FVS) by node-state override forces the network into one of its attractors (long-term dynamic behaviors). The FVS is often composed of more nodes than can be realistically manipulated in a system; for example, only up to three nodes can be controlled in intracellular networks, while their FVS may contain more than 10 nodes. Thus, we developed an approach to rank subsets of the FVS on Boolean models of intracellular networks using topological, dynamics-independent measures. We investigated the use of seven topological prediction measures sorted into three categories-centrality measures, propagation measures, and cycle-based measures. Using each measure, every subset was ranked and then evaluated against two dynamics-based metrics that measure the ability of interventions to drive the system toward or away from its attractors: To Control and Away Control. After examining an array of biological networks, we found that the FVS subsets that ranked in the top according to the propagation metrics can most effectively control the network. This result was independently corroborated on a second array of different Boolean models of biological networks. Consequently, overriding the entire FVS is not required to drive a biological network to one of its attractors, and this method provides a way to reliably identify effective FVS subsets without the knowledge of the network dynamics.

Transcriptional and post-transcriptional control of epithelial-mesenchymal plasticity: why so many regulators?

The dynamic transition between epithelial-like and mesenchymal-like cell states has been a focus for extensive investigation for decades, reflective of the importance of Epithelial-Mesenchymal Transition (EMT) through development, in the adult, and the contributing role EMT has to pathologies including metastasis and fibrosis. Not surprisingly, regulation of the complex genetic networks that underlie EMT have been attributed to multiple transcription factors and microRNAs. What is surprising, however, are the sheer number of different regulators (hundreds of transcription factors and microRNAs) for which critical roles have been described. This review seeks not to collate these studies, but to provide a perspective on the fundamental question of whether it is really feasible that so many regulators play important roles and if so, what does this tell us about EMT and more generally, the genetic machinery that controls complex biological processes.

The identifiability of gene regulatory networks: the role of observation data

Identifying gene regulatory networks (GRN) from observation data is significant to understand biological systems. Conventional studies focus on improving the performance of identification algorithms. However, besides algorithm performance, the GRN identification is strongly depended on the observation data. In this work, for three GRN S-system models, three observation data collection schemes are used to perform the identifiability test procedure. A modified genetic algorithm-particle swarm optimization algorithm is proposed to implement this task, including the multi-level mutation operation and velocity limitation strategy. The results show that, in scheme 1 (starting from a special initial condition), the GRN systems are of identifiability using the sufficient transient observation data. In scheme 2, the observation data are short of sufficient system dynamic. The GRN systems are not of identifiability even though the state trajectories can be reproduced. As a special case of scheme 2, i.e., the steady-state observation data, the equilibrium point analysis is given to explain why it is infeasible for GRN identification. In schemes 1 and 2, the observation data are obtained from zero-input GRN systems, which will evolve to the steady state at last. The sufficient transient observation data in scheme 1 can be obtained by changing the experimental conditions. Additionally, the valid observation data can be also obtained by means of adding impulse excitation signal into GRN systems (scheme 3). Consequently, the GRN systems are identifiable using scheme 3. Owing to its universality and simplicity, these results provide a guide for biologists to collect valid observation data for identifying GRNs and to further understand GRN dynamics.

Epigenetic Reprogramming and Landscape of Transcriptomic Interactions: Impending Therapeutic Interference of Triple-Negative Breast Cancer in Molecular Medicine

The mechanisms governing the development and progression of cancers are believed to be the consequence of hereditary deformities and epigenetic modifications. Accordingly, epigenetics has become an incredible and progressively explored field of research to discover better prevention and therapy for neoplasia, especially triple-negative breast cancer (TNBC). It represents 15-20% of all invasive breast cancers and will, in general, have bellicose histological highlights and poor clinical outcomes. In the early phases of triple-negative breast carcinogenesis, epigenetic deregulation modifies chromatin structure and influences the plasticity of cells. It up-keeps the oncogenic reprogramming of malignant progenitor cells with the acquisition of unrestrained selfrenewal capacities. Genomic impulsiveness in TNBC prompts mutations, copy number variations, as well as genetic rearrangements, while epigenetic remodeling includes an amendment by DNA methylation, histone modification, and noncoding RNAs of gene expression profiles. It is currently evident that epigenetic mechanisms assume a significant part in the pathogenesis, maintenance, and therapeutic resistance of TNBC. Although TNBC is a heterogeneous malaise that is perplexing to describe and treat, the ongoing explosion of genetic and epigenetic research will help to expand these endeavors. Latest developments in transcriptome analysis have reformed our understanding of human diseases, including TNBC at the molecular medicine level. It is appealing to envision transcriptomic biomarkers to comprehend tumor behavior more readily regarding its cellular microenvironment. Understanding these essential biomarkers and molecular changes will propel our capability to treat TNBC adequately. This review will depict the different aspects of epigenetics and the landscape of transcriptomics in triple-negative breast carcinogenesis and their impending application for diagnosis, prognosis, and treatment decision with the view of molecular medicine.

WebMaBoSS: A Web Interface for Simulating Boolean Models Stochastically

WebMaBoSS is an easy-to-use web interface for conversion, storage, simulation and analysis of Boolean models that allows to get insight from these models without any specific knowledge of modeling or coding. It relies on an existing software, MaBoSS, which simulates Boolean models using a stochastic approach: it applies continuous time Markov processes over the Boolean network. It was initially built to fill the gap between Boolean and continuous formalisms, i.e., providing semi-quantitative results using a simple representation with a minimum number of parameters to fit. The goal of WebMaBoSS is to simplify the use and the analysis of Boolean models coping with two main issues: 1) the simulation of Boolean models of intracellular processes with MaBoSS, or any modeling tool, may appear as non-intuitive for non-experts; 2) the simulation of already-published models available in current model databases (e.g., Cell Collective, BioModels) may require some extra steps to ensure compatibility with modeling tools such as MaBoSS. With WebMaBoSS, new models can be created or imported directly from existing databases. They can then be simulated, modified and stored in personal folders. Model simulations are performed easily, results visualized interactively, and figures can be exported in a preferred format. Extensive model analyses such as mutant screening or parameter sensitivity can also be performed. For all these tasks, results are stored and can be subsequently filtered to look for specific outputs. This web interface can be accessed at the address: https://maboss.curie.fr/webmaboss/ and deployed locally using docker. This application is open-source under LGPL license, and available at https://github.com/sysbio-curie/WebMaBoSS.

The involvement of epithelial-to-mesenchymal transition in doxorubicin resistance: Possible molecular targets

Considering the fact that cancer cells can switch among various molecular pathways and mechanisms to ensure their progression, chemotherapy is no longer effective enough in cancer therapy. As an anti-tumor agent, doxorubicin (DOX) is derived from Streptomyces peucetius and can induce cytotoxicity by binding to topoisomerase enzymes to suppress DNA replication, leading to apoptosis and cell cycle arrest. However, efficacy of DOX in suppressing cancer progression is restricted by development of drug resistance. Cancer cells elevate their metastasis in triggering DOX resistance. The epithelial-to-mesenchymal transition (EMT) mechanism participates in transforming epithelial cells into mesenchymal cells that have fibroblast-like features. The EMT diminishes intercellular adhesion and enhances migration of cells that are necessary for carcinogenesis. Various oncogenic molecular pathways stimulate EMT in cancer. EMT can induce DOX resistance, and in this way, upstream mediators such as ZEB proteins, microRNAs, Twist1 and TGF-β play a significant role. Identification of molecular pathways involved in EMT regulation and DOX resistance has resulted in using gene therapy such as microRNA transfection and siRNA in overcoming chemoresistance. Furthermore, curcumin and formononetin, owing to their cytotoxicity against cancer cells, can suppress EMT in mediating DOX sensitivity. For promoting efficacy in DOX sensitivity, nanoparticles have been developed for boosting ability in EMT inhibition.

Gene-circuit therapy on the horizon: synthetic biology tools for engineered therapeutics

Therapeutic genome modification requires precise control over the introduced therapeutic functions. Current approaches of gene and cell therapy fail to deliver such command and rely on semi-quantitative methods with limited influence on timing, contextuality and levels of transgene expression, and hence on therapeutic function. Synthetic biology offers new opportunities for quantitative functionality in designing therapeutic systems and their components. Here, we discuss synthetic biology tools in their therapeutic context, with examples of proof-of-principle and clinical applications of engineered synthetic biomolecules and higher-order functional systems, i.e. gene circuits. We also present the prospects of future development towards advanced gene-circuit therapy.

Immunofluorescence Image Feature Analysis and Phenotype Scoring Pipeline for Distinguishing Epithelial-Mesenchymal Transition

Epithelial–mesenchymal transition (EMT) is an essential biological process, also implicated in pathological settings such as cancer metastasis, in which epithelial cells transdifferentiate into mesenchymal cells. We devised an image analysis pipeline to distinguish between tissues comprised of epithelial and mesenchymal cells, based on extracted features from immunofluorescence images of differing biochemical markers. Mammary epithelial cells were cultured with 0 (control), 2, 4, or 10 ng/mL TGF-β1, a well-established EMT-inducer. Cells were fixed, stained, and imaged for E-cadherin, actin, fibronectin, and nuclei via immunofluorescence microscopy. Feature selection was performed on different combinations of individual cell markers using a Bag-of-Features extraction. Control and high-dose images comprised the training data set, and the intermediate dose images comprised the testing data set. A feature distance analysis was performed to quantify differences between the treatment groups. The pipeline was successful in distinguishing between control (epithelial) and the high-dose (mesenchymal) groups, as well as demonstrating progress along the EMT process in the intermediate dose groups. Validation using quantitative PCR (qPCR) demonstrated that biomarker expression measurements were well-correlated with the feature distance analysis. Overall, we identified image pipeline characteristics for feature extraction and quantification of immunofluorescence images to distinguish progression of EMT.

Concomitant attenuation of HMGCR expression and activity enhances the growth inhibitory effect of atorvastatin on TGF-β-treated epithelial cancer cells

Epithelial-mesenchymal transition (EMT) in primary tumor cells is a key prerequisite for metastasis initiation. Statins, cholesterol-lowering drugs, can delay metastasis formation in vivo and attenuate the growth and proliferation of tumor cells in vitro. The latter effect is stronger in tumor cells with a mesenchymal-like phenotype than in those with an epithelial one. However, the effect of statins on epithelial cancer cells treated with EMT-inducing growth factors such as transforming growth factor-β (TGF-β) remains unclear. Here, we examined the effect of atorvastatin on two epithelial cancer cell lines following TGF-β treatment. Atorvastatin-induced growth inhibition was stronger in TGF-β-treated cells than in cells not thusly treated. Moreover, treatment of cells with atorvastatin prior to TGF-β treatment enhanced this effect, which was further potentiated by the simultaneous reduction in the expression of the statin target enzyme, 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR). Dual pharmacological targeting of HMGCR can thus strongly inhibit the growth and proliferation of epithelial cancer cells treated with TGF-β and may also improve statin therapy-mediated attenuation of metastasis formation in vivo.

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.

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.

Impact of Nischarin on EMT regulators in breast cancer cell lines

Nischarin is an integrin-binding protein, which is well known as a novel tumor suppressor. In breast cancer, Nischarin serves a critical role in breast cancer cell migration and invasion. However, the molecular mechanism underlying the role of Nischarin remains unclear. Recent findings have demonstrated that epithelial-mesenchymal transition (EMT) increases the capacity of cell migration and invasion. As a member of the integrin family, it was hypothesized that Nischarin may regulate cellular processes via various signaling pathways associated with the EMT process. The present study detected the mRNA levels of EMT regulators via reverse transcription-quantitative PCR and related protein levels via western blotting in breast cancer cells, following NISCH-overexpression and -knockdown. The results demonstrated that Nischarin inhibits cell proliferation, migration and invasion in breast cancer cells. Furthermore, when the NISCH gene was overexpressed, the relative mRNA level of E-cadherin was increased, while the relative mRNA levels of several transcription factors, such as Snail, ZEB1, N-cadherin, Slug, Twist1 and vimentin, decreased. When NISCH was silenced, these results were reversed. The present results demonstrated that Nischarin suppresses cell migration and invasion via inhibiting the EMT process.

Boolean model of anchorage dependence and contact inhibition points to coordinated inhibition but semi-independent induction of proliferation and migration

Epithelial cells respond to their physical neighborhood with mechano-sensitive behaviors required for development and tissue maintenance. These include anchorage dependence, matrix stiffness-dependent proliferation, contact inhibition of proliferation and migration, and collective migration that balances cell crawling with the maintenance of cell junctions. While required for development and tissue repair, these coordinated responses to the microenvironment also contribute to cancer metastasis. Predictive models of the signaling networks that coordinate these behaviors are critical in controlling cell behavior to halt disease. Here we propose a Boolean regulatory network model that synthesizes mechanosensitive signaling that links anchorage to a matrix of varying stiffness and cell density sensing to contact inhibition, proliferation, migration, and apoptosis. Our model can reproduce anchorage dependence and anoikis, detachment-induced cytokinesis errors, the effect of matrix stiffness on proliferation, and contact inhibition of proliferation and migration by two mechanisms that converge on the YAP transcription factor. In addition, we offer testable predictions related to cell cycle-dependent anoikis sensitivity, the molecular requirements for abolishing contact inhibition, and substrate stiffness dependent expression of the catalytic subunit of PI3K. Moreover, our model predicts heterogeneity in migratory vs. non-migratory phenotypes in sub-confluent monolayers, and co-inhibition but semi-independent induction of proliferation vs. migration as a function of cell density and mitogenic stimulation. Our model serves as a stepping-stone towards modeling mechanosensitive routes to the epithelial to mesenchymal transition, capturing the effects of the mesenchymal state on anoikis resistance, and understanding the balance between migration versus proliferation at each stage of the epithelial to mesenchymal transition.

Local and global features of genetic networks supporting a phenotypic switch

Phenotypic switches are associated with alterations in the cell's gene expression profile and are vital to many aspects of biology. Previous studies have identified local motifs of the genetic regulatory network that could underlie such switches. Recent advancements allowed the study of networks at the global, many-gene, level; however, the relationship between the local and global scales in giving rise to phenotypic switches remains elusive. In this work, we studied the epithelial-mesenchymal transition (EMT) using a gene regulatory network model. This model supports two clusters of stable steady-states identified with the epithelial and mesenchymal phenotypes, and a range of intermediate less stable hybrid states, whose importance in cancer has been recently highlighted. Using an array of network perturbations and quantifying the resulting landscape, we investigated how features of the network at different levels give rise to these landscape properties. We found that local connectivity patterns affect the landscape in a mostly incremental manner; in particular, a specific previously identified double-negative feedback motif is not required when embedded in the full network, because the landscape is maintained at a global level. Nevertheless, despite the distributed nature of the switch, it is possible to find combinations of a few local changes that disrupt it. At the level of network architecture, we identified a crucial role for peripheral genes that act as incoming signals to the network in creating clusters of states. Such incoming signals are a signature of modularity and are expected to appear also in other biological networks. Hybrid states between epithelial and mesenchymal arise in the model due to barriers in the interaction between genes, causing hysteresis at all connections. Our results suggest emergent switches can neither be pinpointed to local motifs, nor do they arise as typical properties of random network ensembles. Rather, they arise through an interplay between the nature of local interactions, and the core-periphery structure induced by the modularity of the cell.

Cell plasticity in cancer cell populations

In this review, we propose a recension of biological observations on plasticity in cancer cell populations and discuss theoretical considerations about their mechanisms.

How Do Cells Adapt? Stories Told in Landscapes

Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.

Identifying inhibitors of epithelial-mesenchymal plasticity using a network topology-based approach

Metastasis is the cause of over 90% of cancer-related deaths. Cancer cells undergoing metastasis can switch dynamically between different phenotypes, enabling them to adapt to harsh challenges, such as overcoming anoikis and evading immune response. This ability, known as phenotypic plasticity, is crucial for the survival of cancer cells during metastasis, as well as acquiring therapy resistance. Various biochemical networks have been identified to contribute to phenotypic plasticity, but how plasticity emerges from the dynamics of these networks remains elusive. Here, we investigated the dynamics of various regulatory networks implicated in Epithelial–mesenchymal plasticity (EMP)—an important arm of phenotypic plasticity—through two different mathematical modelling frameworks: a discrete, parameter-independent framework (Boolean) and a continuous, parameter-agnostic modelling framework (RACIPE). Results from either framework in terms of phenotypic distributions obtained from a given EMP network are qualitatively similar and suggest that these networks are multi-stable and can give rise to phenotypic plasticity. Neither method requires specific kinetic parameters, thus our results emphasize that EMP can emerge through these networks over a wide range of parameter sets, elucidating the importance of network topology in enabling phenotypic plasticity. Furthermore, we show that the ability to exhibit phenotypic plasticity correlates positively with the number of positive feedback loops in a given network. These results pave a way toward an unorthodox network topology-based approach to identify crucial links in a given EMP network that can reduce phenotypic plasticity and possibly inhibit metastasis—by reducing the number of positive feedback loops.

CBX7 binds the E-box to inhibit TWIST-1 function and inhibit tumorigenicity and metastatic potential

Deaths from ovarian cancer usually occur when patients succumb to overwhelmingly numerous and widespread micrometastasis. Whereas epithelial-mesenchymal transition is required for epithelial ovarian cancer cells to acquire metastatic potential, the cellular phenotype at secondary sites and the mechanisms required for the establishment of metastatic tumors are not fully determined. Using in vitro and in vivo models we show that secondary epithelial ovarian cancer cells (sEOC) do not fully reacquire the molecular signature of the primary epithelial ovarian cancer cells from which they are derived. Despite displaying an epithelial morphology, sEOC maintains a high expression of the mesenchymal effector, TWIST-1. TWIST-1 is however transcriptionally nonfunctional in these cells as it is precluded from binding its E-box by the PcG protein, CBX7. Deletion of CBX7 in sEOC was sufficient to reactivate TWIST-1-induced transcription, prompt mesenchymal transformation, and enhanced tumorigenicity in vivo. This regulation allows secondary tumors to achieve an epithelial morphology while conferring the advantage of prompt reversal to a mesenchymal phenotype upon perturbation of CBX7. We also describe a subclassification of ovarian tumors based on CBX7 and TWIST-1 expression, which predicts clinical outcomes and patient prognosis.

Multiplexed Gene Expression Tuning with Orthogonal Synthetic Gene Circuits

Understanding the individual and joint contribution of multiple protein levels toward a phenotype requires precise and tunable multigene expression control. Here we introduce a pair of mammalian synthetic gene circuits that linearly and orthogonally control the expression of two reporter genes in mammalian cells with low variability in response to chemical inducers introduced into the growth medium. These gene expression systems can be used to simultaneously probe the individual and joint effects of two gene product concentrations on a cellular phenotype in basic research or biomedical applications.

Hybrid Epithelial-Mesenchymal Phenotypes Are Controlled by Microenvironmental Factors

Epithelial-to-mesenchymal transition (EMT) has been associated with cancer cell heterogeneity, plasticity, and metastasis. However, the extrinsic signals supervising these phenotypic transitions remain elusive. To assess how selected microenvironmental signals control cancer-associated phenotypes along the EMT continuum, we defined a logical model of the EMT cellular network that yields qualitative degrees of cell adhesions by adherens junctions and focal adhesions, two features affected during EMT. The model attractors recovered epithelial, mesenchymal, and hybrid phenotypes. Simulations showed that hybrid phenotypes may arise through independent molecular paths involving stringent extrinsic signals. Of particular interest, model predictions and their experimental validations indicated that: (i) stiffening of the extracellular matrix was a prerequisite for cells overactivating FAK_SRC to upregulate SNAIL and acquire a mesenchymal phenotype and (ii) FAK_SRC inhibition of cell-cell contacts through the receptor-type tyrosine-protein phosphatases kappa led to acquisition of a full mesenchymal, rather than a hybrid, phenotype. Altogether, these computational and experimental approaches allow assessment of critical microenvironmental signals controlling hybrid EMT phenotypes and indicate that EMT involves multiple molecular programs. SIGNIFICANCE: A multidisciplinary study sheds light on microenvironmental signals controlling cancer cell plasticity along EMT and suggests that hybrid and mesenchymal phenotypes arise through independent molecular paths.

Bidirectional interplay between physical and biological approaches on studying the epithelial-to-mesenchymal transition

Physical approaches have made notable contributions to the study of epithelial-to-mesenchymal transition (EMT), and EMT serves as a model system for advancing physics theories. A collection of reviews and original research papers are included in this special issue.

Noise-reducing optogenetic negative-feedback gene circuits in human cells

Gene autorepression is widely present in nature and is also employed in synthetic biology, partly to reduce gene expression noise in cells. Optogenetic systems have recently been developed for controlling gene expression levels in mammalian cells, but most have utilized activator-based proteins, neglecting negative feedback except for in silico control. Here, we engineer optogenetic gene circuits into mammalian cells to achieve noise-reduction for precise gene expression control by genetic, in vitro negative feedback. We build a toolset of these noise-reducing Light-Inducible Tuner (LITer) gene circuits using the TetR repressor fused with a Tet-inhibiting peptide (TIP) or a degradation tag through the light-sensitive LOV2 protein domain. These LITers provide a range of nearly 4-fold gene expression control and up to 5-fold noise reduction from existing optogenetic systems. Moreover, we use the LITer gene circuit architecture to control gene expression of the cancer oncogene KRAS(G12V) and study its downstream effects through phospho-ERK levels and cellular proliferation. Overall, these novel LITer optogenetic platforms should enable precise spatiotemporal perturbations for studying multicellular phenotypes in developmental biology, oncology and other biomedical fields of research.