RACER-m leverages structural features for sparse T cell specificity prediction.

Reliable prediction of T cell specificity against antigenic signatures is a formidable task, complicated by the immense diversity of T cell receptor and antigen sequence space and the resulting limited availability of training sets for inferential models. Recent modeling efforts have demonstrated the advantage of incorporating structural information to overcome the need for extensive training sequence data, yet disentangling the heterogeneous TCR-antigen interface to accurately predict MHC-allele-restricted TCR-peptide interactions has remained challenging. Here, we present RACER-m, a coarse-grained structural model leveraging key biophysical information from the diversity of publicly available TCR-antigen crystal structures. Explicit inclusion of structural content substantially reduces the required number of training examples and maintains reliable predictions of TCR-recognition specificity and sensitivity across diverse biological contexts. Our model capably identifies biophysically meaningful point-mutant peptides that affect binding affinity, distinguishing its ability in predicting TCR specificity of point-mutants from alternative sequence-based methods. Its application is broadly applicable to studies involving both closely related and structurally diverse TCR-peptide pairs.

Dual role of CASP8AP2/FLASH in regulating epithelial-to-mesenchymal transition plasticity (EMP).

BACKGROUND: Epithelial-to-mesenchymal transition (EMT) is a developmental program that consists of the loss of epithelial features concomitant with the acquisition of mesenchymal features. Activation of EMT in cancer facilitates the acquisition of aggressive traits and cancer invasion. EMT plasticity (EMP), the dynamic transition between multiple hybrid states in which cancer cells display both epithelial and mesenchymal markers, confers survival advantages for cancer cells in constantly changing environments during metastasis. METHODS: RNAseq analysis was performed to assess genome-wide transcriptional changes in cancer cells depleted for histone regulators FLASH, NPAT, and SLBP. Quantitative PCR and Western blot were used for the detection of mRNA and protein levels. Computational analysis was performed on distinct sets of genes to determine the epithelial and mesenchymal score in cancer cells and to correlate FLASH expression with EMT markers in the CCLE collection. RESULTS: We demonstrate that loss of FLASH in cancer cells gives rise to a hybrid E/M phenotype with high epithelial scores even in the presence of TGFß, as determined by computational methods using expression of predetermined sets of epithelial and mesenchymal genes. Multiple genes involved in cell-cell junction formation are similarly specifically upregulated in FLASH-depleted cells, suggesting that FLASH acts as a repressor of the epithelial phenotype. Further, FLASH expression in cancer lines is inversely correlated with the epithelial score. Nonetheless, subsets of mesenchymal markers were distinctly up-regulated in FLASH, NPAT, or SLBP-depleted cells. CONCLUSIONS: The ZEB1low/SNAILhigh/E-cadherinhigh phenotype described in FLASH-depleted cancer cells is driving a hybrid E/M phenotype in which epithelial and mesenchymal markers coexist.

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.

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.

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.

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.

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.

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.

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.