Bridging the gap between collective motility and epithelial-mesenchymal transitions through the active finite voronoi model.

We introduce an active version of the recently proposed finite Voronoi model of epithelial tissue. The resultant Active Finite Voronoi (AFV) model enables the study of both confluent and non-confluent geometries and transitions between them, in the presence of active cells. Our study identifies six distinct phases, characterized by aggregation-segregation, dynamical jamming-unjamming, and epithelial-mesenchymal transitions (EMT), thereby extending the behavior beyond that observed in previously studied vertex-based models. The AFV model with rich phase diagram provides a cohesive framework that unifies the well-observed progression to collective motility via unjamming with the intricate dynamics enabled by EMT. This approach should prove useful for challenges in developmental biology systems as well as the complex context of cancer metastasis. The simulation code is also provided.

Modeling the Role of Immune Cell Conversion in the Tumor-Immune Microenvironment.

Tumors develop in a complex physical, biochemical, and cellular milieu, referred to as the tumor microenvironment. Of special interest is the set of immune cells that reciprocally interact with the tumor, the tumor-immune microenvironment (TIME). The diversity of cell types and cell-cell interactions in the TIME has led researchers to apply concepts from ecology to describe the dynamics. However, while tumor cells are known to induce immune cells to switch from anti-tumor to pro-tumor phenotypes, this type of ecological interaction has been largely overlooked. To address this gap in cancer modeling, we develop a minimal, ecological model of the TIME with immune cell conversion, to highlight this important interaction and explore its consequences. A key finding is that immune conversion increases the range of parameters supporting a co-existence phase in which the immune system and the tumor reach a stalemate. Our results suggest that further investigation of the consequences of immune cell conversion, using detailed, data-driven models, will be critical for greater understanding of TIME dynamics.

Physics of smart active matter: integrating active matter and control to gain insights into living systems.

We offer our opinion on the benefits of integration of insights from active matter physics with principles of regulatory interactions and control to develop a field we term "smart active matter". This field can provide insight into important principles in living systems as well as aid engineering of responsive, robust and functional collectives.

Infiltration of tumor spheroids by activated immune cells.

Recent years have seen a tremendous growth of interest in understanding the role that the adaptive immune system could play in interdicting tumor progression. In this context, it has been shown that the density of adaptive immune cells inside a solid tumor serves as a favorable prognostic marker across different types of cancer. The exact mechanisms underlying the degree of immune cell infiltration is largely unknown. Here, we quantify the temporal dynamics of the density profile of activated immune cells around a solid tumor spheroid. We propose a computational model incorporating immune cells with active, persistent movement and a proliferation rate that depends on the presence of cancer cells, and show that the model able to reproduce semi-quantitatively the experimentally measured infiltration profile. Studying the density distribution of immune cells inside a solid tumor can help us better understand immune trafficking in the tumor micro-environment, hopefully leading towards novel immunotherapeutic strategies.

The alternate ligand Jagged enhances the robustness of Notch signaling patterns.

The Notch pathway, an example of juxtacrine signaling, is an evolutionary conserved cell-cell communication mechanism. It governs emergent spatiotemporal patterning in tissues during development, wound healing and tumorigenesis. Communication occurs when Notch receptors of one cell bind to either of its ligands, Delta/Jagged of the neighboring cell. In general, Delta-mediated signaling drives neighboring cells to have an opposite fate (lateral inhibition) whereas Jagged-mediated signaling drives cells to maintain similar fates (lateral induction). Here, by deriving and solving a reduced set of 12 coupled ordinary differential equations for the Notch-Delta-Jagged system on a hexagonal grid of cells, we determine the allowed states across different parameter sets. We also show that Jagged (at low dose) acts synergistically with Delta to enable more robust pattern formation by making the neighboring cell states more distinct from each other, despite its lateral induction property. Our findings extend our understanding of the possible synergistic role of Jagged with Delta which had been previously proposed through experiments and models in the context of chick inner ear development. Finally, we show how Jagged can help to expand the bistable (both uniform and hexagon phases are stable) region, where a local perturbation can spread over time in an ordered manner to create a biologically relevant, perfectly ordered lateral inhibition pattern.

Optimal cancer evasion in a dynamic immune microenvironment generates diverse post-escape tumor antigenicity profiles.

The failure of cancer treatments, including immunotherapy, continues to be a major obstacle in preventing durable remission. This failure often results from tumor evolution, both genotypic and phenotypic, away from sensitive cell states. Here, we propose a mathematical framework for studying the dynamics of adaptive immune evasion that tracks the number of tumor-associated antigens available for immune targeting. We solve for the unique optimal cancer evasion strategy using stochastic dynamic programming and demonstrate that this policy results in increased cancer evasion rates compared to a passive, fixed strategy. Our foundational model relates the likelihood and temporal dynamics of cancer evasion to features of the immune microenvironment, where tumor immunogenicity reflects a balance between cancer adaptation and host recognition. In contrast with a passive strategy, optimally adaptive evaders navigating varying selective environments result in substantially heterogeneous post-escape tumor antigenicity, giving rise to immunogenically hot and cold tumors.

Modeling collective cell behavior in cancer: Perspectives from an interdisciplinary conversation.

Collective cell behavior contributes to all stages of cancer progression. Understanding how collective behavior emerges through cell-cell interactions and decision-making will advance our understanding of cancer biology and provide new therapeutic approaches. Here, we summarize an interdisciplinary discussion on multicellular behavior in cancer, draw lessons from other scientific disciplines, and identify future directions.

Addressing the genetic/nongenetic duality in cancer with systems biology.

The dogma that cancer is a genetic disease is being questioned. Recent findings suggest that genetic/nongenetic duality is necessary for cancer progression. A think tank organized by the Shraman Foundation's Institute for Theoretical Biology compiled key challenges and opportunities that theoreticians, experimentalists, and clinicians can explore from a systems biology perspective to provide a better understanding of the disease as well as help discover new treatment options and therapeutic strategies.

Epigenetic memory acquired during long-term EMT induction governs the recovery to the epithelial state.

Epithelial-mesenchymal transition (EMT) and its reverse mesenchymal-epithelial transition (MET) are critical during embryonic development, wound healing and cancer metastasis. While phenotypic changes during short-term EMT induction are reversible, long-term EMT induction has been often associated with irreversibility. Here, we show that phenotypic changes seen in MCF10A cells upon long-term EMT induction by TGFß need not be irreversible, but have relatively longer time scales of reversibility than those seen in short-term induction. Next, using a phenomenological mathematical model to account for the chromatin-mediated epigenetic silencing of the miR-200 family by ZEB family, we highlight how the epigenetic memory gained during long-term EMT induction can slow the recovery to the epithelial state post-TGFß withdrawal. Our results suggest that epigenetic modifiers can govern the extent and time scale of EMT reversibility and advise caution against labelling phenotypic changes seen in long-term EMT induction as 'irreversible'.

Minimal frustration underlies the usefulness of incomplete regulatory network models in biology.

Regulatory networks as large and complex as those implicated in cell-fate choice are expected to exhibit intricate, very high-dimensional dynamics. Cell-fate choice, however, is a macroscopically simple process. Additionally, regulatory network models are almost always incomplete and/or inexact, and do not incorporate all the regulators and interactions that may be involved in cell-fate regulation. In spite of these issues, regulatory network models have proven to be incredibly effective tools for understanding cell-fate choice across contexts and for making useful predictions. Here, we show that minimal frustration-a feature of biological networks across contexts but not of random networks-can compel simple, low-dimensional steady-state behavior even in large and complex networks. Moreover, the steady-state behavior of minimally frustrated networks can be recapitulated by simpler networks such as those lacking many of the nodes and edges and those that treat multiple regulators as one. The present study provides a theoretical explanation for the success of network models in biology and for the challenges in network inference.

Decoding the coupled decision-making of the epithelial-mesenchymal transition and metabolic reprogramming in cancer.

Cancer metastasis relies on an orchestration of traits driven by different interacting functional modules, including metabolism and epithelial-mesenchymal transition (EMT). During metastasis, cancer cells can acquire a hybrid metabolic phenotype (W/O) by increasing oxidative phosphorylation without compromising glycolysis and they can acquire a hybrid epithelial/mesenchymal (E/M) phenotype by engaging EMT. Both the W/O and E/M states are associated with high metastatic potentials, and many regulatory links coupling metabolism and EMT have been identified. Here, we investigate the coupled decision-making networks of metabolism and EMT. Their crosstalk can exhibit synergistic or antagonistic effects on the acquisition and stability of different coupled metabolism-EMT states. Strikingly, the aggressive E/M-W/O state can be enabled and stabilized by the crosstalk irrespective of these hybrid states' availability in individual metabolism or EMT modules. Our work emphasizes the mutual activation between metabolism and EMT, providing an important step toward understanding the multifaceted nature of cancer metastasis.

Epigenetic factor competition reshapes the EMT landscape.

The emergence of and transitions between distinct phenotypes in isogenic cells can be attributed to the intricate interplay of epigenetic marks, external signals, and gene-regulatory elements. These elements include chromatin remodelers, histone modifiers, transcription factors, and regulatory RNAs. Mathematical models known as gene-regulatory networks (GRNs) are an increasingly important tool to unravel the workings of such complex networks. In such models, epigenetic factors are usually proposed to act on the chromatin regions directly involved in the expression of relevant genes. However, it has been well-established that these factors operate globally and compete with each other for targets genome-wide. Therefore, a perturbation of the activity of a regulator can redistribute epigenetic marks across the genome and modulate the levels of competing regulators. In this paper, we propose a conceptual and mathematical modeling framework that incorporates both local and global competition effects between antagonistic epigenetic regulators, in addition to local transcription factors, and show the counterintuitive consequences of such interactions. We apply our approach to recent experimental findings on the epithelial-mesenchymal transition (EMT). We show that it can explain the puzzling experimental data, as well as provide verifiable predictions.

Phenomenological Approach to Cancer Cell Persistence.

Drug persistence is a phenomenon by which a small percentage of cancer cells survive the presentation of targeted therapy by transitioning to a quiescent state. Eventually some of these persister cells can transition back to an active growing state and give rise to resistant tumors. Here we introduce a quantitative genetics approach to drug-exposed populations of cancer cells in order to interpret recent experimental data regarding inheritance of persister probability. Our results indicate that alternating periods of drug treatment and drug removal may not be an effective strategy for eliminating persisters.

A physical model for dynamic assembly of human salivary stem/progenitor microstructures.

The in vitro reconstructions of human salivary glands in service of their eventual medical use represent a challenge for tissue engineering. Here, we present a theoretical approach to the dynamical formation of acinar structures from human salivary cells, focusing on observed stick-slip radial expansion as well as possible growth instabilities. Our findings demonstrate the critical importance of basement membrane remodeling in controlling the growth process.

Contact map dependence of a T-cell receptor binding repertoire.

The T-cell arm of the adaptive immune system provides the host protection against unknown pathogens by discriminating between host and foreign material. This discriminatory capability is achieved by the creation of a repertoire of cells each carrying a T-cell receptor (TCR) specific to non-self-antigens displayed as peptides bound to the major histocompatibility complex (pMHC). The understanding of the dynamics of the adaptive immune system at a repertoire level is complex, due to both the nuanced interaction of a TCR-pMHC pair and to the number of different possible TCR-pMHC pairings, making computationally exact solutions currently unfeasible. To gain some insight into this problem, we study an affinity-based model for TCR-pMHC binding in which a crystal structure is used to generate a distance-based contact map that weights the pairwise amino acid interactions. We find that the TCR-pMHC binding energy distribution strongly depends both on the number of contacts and the repeat structure allowed by the topology of the contact map of choice; this in turn influences T-cell recognition probability during negative selection, with higher variances leading to higher survival probabilities. In addition, we quantify the degree to which neoantigens with mutations in sites with higher contacts are recognized at a higher rate.

Lack of Correlation between Landscape Geometry and Transition Rates.

Biological cells can exist in a variety of distinct phenotypes, determined by the steady-state solutions of genetic networks governing their cell fate. A popular way of representing these states relies on the creation of landscape related to the relative occupation of these states. It is often assumed that this landscape offers direct information regarding the state-to-state transition rates, suggesting that these are related to barrier heights separating landscape minima. Here, we study a toggle triad network exhibiting multistability and directly demonstrate the lack of any direct correlation between properties of the landscape and corresponding transition rates.

Nrf2 Modulates the Hybrid Epithelial/Mesenchymal Phenotype and Notch Signaling During Collective Cancer Migration.

Hybrid epithelial/mesenchymal cells (E/M) are key players in aggressive cancer metastasis. It remains a challenge to understand how these cell states, which are mostly non-existent in healthy tissue, become stable phenotypes participating in collective cancer migration. The transcription factor Nrf2, which is associated with tumor progression and resistance to therapy, appears to be central to this process. Here, using a combination of immunocytochemistry, single cell biosensors, and computational modeling, we show that Nrf2 functions as a phenotypic stability factor for hybrid E/M cells by inhibiting a complete epithelial-mesenchymal transition (EMT) during collective cancer migration. We also demonstrate that Nrf2 and EMT signaling are spatially coordinated near the leading edge. In particular, computational analysis of an Nrf2-EMT-Notch network and experimental modulation of Nrf2 by pharmacological treatment or CRISPR/Cas9 gene editing reveal that Nrf2 stabilizes a hybrid E/M phenotype which is maximally observed in the interior region immediately behind the leading edge. We further demonstrate that the Nrf2-EMT-Notch network enhances Dll4 and Jagged1 expression at the leading edge, which correlates with the formation of leader cells and protruding tips. Altogether, our results provide direct evidence that Nrf2 acts as a phenotypic stability factor in restricting complete EMT and plays an important role in coordinating collective cancer migration.

Changes in Triple-Negative Breast Cancer Molecular Subtypes in Patients Without Pathologic Complete Response After Neoadjuvant Systemic Chemotherapy.

PURPOSE: Lehmann et al have identified four molecular subtypes of triple-negative breast cancer (TNBC)-basal-like (BL) 1, BL2, mesenchymal (M), and luminal androgen receptor-and an immunomodulatory (IM) gene expression signature modifier. Our group previously showed that the response of TNBC to neoadjuvant systemic chemotherapy (NST) differs by molecular subtype, but whether NST affects the subtype was unknown. Here, we tested the hypothesis that in patients without pathologic complete response, TNBC subtypes can change after NST. Moreover, in cases with the changed subtype, we determined whether epithelial-to-mesenchymal transition (EMT) had occurred. MATERIALS AND METHODS: From the Pan-Pacific TNBC Consortium data set containing TNBC patient samples from four countries, we examined 64 formalin-fixed, paraffin-embedded pairs of matched pre- and post-NST tumor samples. The TNBC subtype was determined using the TNBCtype-IM assay. We analyzed a partial EMT gene expression scoring metric using mRNA data. RESULTS: Of the 64 matched pairs, 36 (56%) showed a change in the TNBC subtype after NST. The most frequent change was from BL1 to M subtypes (38%). No tumors changed from M to BL1. The IM signature was positive in 14 (22%) patients before NST and eight (12.5%) patients after NST. The EMT score increased after NST in 28 (78%) of the 36 patients with the changed subtype (v 39% of the 28 patients without change; P = .002254). CONCLUSION: We report, to our knowledge, for the first time that the TNBC molecular subtype and IM signature frequently change after NST. Our results also suggest that EMT is promoted by NST. Our findings may lead to innovative adjuvant therapy strategies in TNBC cases with residual tumor after NST.