Tissue mechanical properties such as rigidity and fluidity, and changes in these properties driven by jamming-unjamming transitions (UJT), have come under recent highlight as mechanical markers of health and disease in various biological processes including cancer. However, most analysis of these mechanical properties and UJT have sidestepped the effect of cellular death and division in these systems. Cellular apoptosis (programmed cell death) and mitosis (cell division) can drive significant changes in tissue properties. The balance between the two is crucial in maintaining tissue function, and an imbalance between the two is seen in situations such as cancer progression, wound healing and necrosis. In this work we investigate the impact of cell death and division on tissue mechanical properties, by incorporating specific mechanosensitive triggers of cell death and division based on the size and geometry of the cell within in silico models of tissue dynamics. Specifically, we look at cell migration, tissue response to external stress, tissue extrusion propensity and self-organization of different cell types within the tissue, as a function of cell death and division and the rules
Understanding deaths and life-death boundaries of cells is a fundamental challenge in biological sciences. In this study, we present a theoretical framework for investigating cell death. We conceptualize cell death as a controllability problem within dynamical systems, and compute the life-death boundary through the development of "stoichiometric rays". This method utilizes enzyme activity as control parameters, exploiting the inherent property of enzymes to enhance reaction rates without shifting equilibrium states. This approach facilitates the efficient evaluation of the global controllability of models. We demonstrate the utility of our framework using its application to a toy metabolic model, where we delineate the life-death boundary. The formulation of cell death through mathematical principles provides a foundation for the theoretical study of cellular mortality.
Cellular responses to death-promoting stimuli typically proceed through a differentiated multistage process, involving a lag phase, extensive death, and potential adaptation. Deregulation of this chain of events is at the root of many diseases. Improper adaptation is particularly important because it allows cell sub-populations to survive even in the continuous presence of death conditions, which results, among others, in the eventual failure of many targeted anticancer therapies. Here, I show that these typical responses arise naturally from the interplay of intracellular variability with a threshold-based control mechanism that detects cellular changes in addition to just the cellular state itself. Implementation of this mechanism in a quantitative model for T-cell apoptosis, a prototypical example of programmed cell death, captures with exceptional accuracy experimental observations for different expression levels of the oncogene Bcl-xL and directly links adaptation with noise in an ATP threshold below which cells die. These results indicate that oncogenes like Bcl-xL, besides regulating absolute death values, can have a novel role as active controllers of cell-cell variability
Death ligand mediated apoptotic activation is a mode of programmed cell death that is widely used in cellular and physiological situations. Interest in studying death ligand induced apoptosis has increased due to the promising role of recombinant soluble forms of death ligands (mainly recombinant TRAIL) in anti-cancer therapy. A clear elucidation of how death ligands activate the type 1 and type 2 apoptotic pathways in healthy and cancer cells may help develop better chemotherapeutic strategies. In this work, we use kinetic Monte Carlo simulations to address the problem of type 1/ type 2 choice in death ligand mediated apoptosis of cancer cells. Our study provides insights into the activation of membrane proximal death module that results from complex interplay between death and decoy receptors. Relative abundance of death and decoy receptors was shown to be a key parameter for activation of the initiator caspases in the membrane module. Increased concentration of death ligands frequently increased the type 1 activation fraction in cancer cells, and, in certain cases changes the signaling phenotype from type 2 to type 1. Results of this study also indicate that inherent differences
BH3 mimetics have been proposed as new anticancer therapeutics. They target anti-apoptotic Bcl-2 proteins, up-regulation of which has been implicated in the resistance of many cancer cells, particularly leukemia and lymphoma cells, to apoptosis. Using probabilistic computational modeling of the mitochondrial pathway of apoptosis, verified by single-cell experimental observations, we develop a model of Bcl-2 inhibition of apoptosis. Our results clarify how Bcl-2 imparts its anti-apoptotic role by increasing the time-to-death and cell-to-cell variability. We also show that although the commitment to death is highly impacted by differences in protein levels at the time of stimulation, inherent stochastic fluctuations in apoptotic signaling are sufficient to induce cell-to-cell variability and to allow single cells to escape death. This study suggests that intrinsic cell-to-cell stochastic variability in apoptotic signaling is sufficient to cause fractional killing of cancer cells after exposure to BH3 mimetics. This is an unanticipated facet of cancer chemoresistance.
Nearly all cell models explicitly or implicitly deal with the biophysical constraints that must be respected for life to persist. Despite this, there is almost no systematicity in how these constraints are implemented, and we lack a principled understanding of how cellular dynamics interact with them and how they originate in actual biology. Computational cell biology will only overcome these concerns once it treats the life-death boundary as a central concept, creating a theory of cellular viability. We lay the foundation for such a development by demonstrating how specific geometric structures can separate regions of qualitatively similar survival outcomes in our models, offering new global organizing principles for cell fate. We also argue that idealized models of emergent individuals offer a tractable way to begin understanding life's intrinsically generated limits.
The crawling motility of many eukaryotic cells is driven by filamentous actin (F-actin), and regulated by a network of signaling proteins and lipids (including small GTPases). The tangle of positive and negative feedback loops gives rise to various experimentally observed dynamic patterns (``actin waves''). Here we consider a recent prototypical model for actin waves in which F-actin exerts negative feedback onto a GTPase. Guided by recent numerical PDE bifurcation analysis in Hughes (2025) and Hughes et al (2026), we explore cell shapes and motility associated with polar, oscillatory, and traveling waves solutions of a mass-conserved partial differential equation (PDE) model. We use Morpheus (cellular Potts) simulations to investigate the implications of such regimes of behavior on the shapes and motion of cells, and on transitions between modes of behavior. The model demonstrates various cell states, including resting (spatially uniform GTPase), polar cells (static ``zones'' of GTPase), and traveling waves along the cell edge. In some parameter regimes, such states can coexist, so that cells can transition from one behavior to another in response to noisy stimuli.
Cell-cell communication is essential for tissue development, regeneration and function, and its disruption can lead to diseases and developmental abnormalities. The revolution of single-cell genomics technologies offers unprecedented insights into cellular identities, opening new avenues to resolve the intricate cellular interactions present in tissue niches. CellPhoneDB is a bioinformatics toolkit designed to infer cell-cell communication by combining a curated repository of bona fide ligand-receptor interactions with a set of computational and statistical methods to integrate them with single-cell genomics data. Importantly, CellPhoneDB captures the multimeric nature of molecular complexes, thus representing cell-cell communication biology faithfully. Here we present CellPhoneDB v5, an updated version of the tool, which offers several new features. Firstly, the repository has been expanded by one-third with the addition of new interactions. These encompass interactions mediated by non-protein ligands such as endocrine hormones and GPCR ligands. Secondly, it includes a differentially expression-based methodology for more tailored interaction queries. Thirdly, it incorporates novel
The process of programmed cell death, namely apoptosis, is a natural mechanism that regulates healthy tissue, multicellular structures, and homeostasis. An improved understanding of apoptosis can significantly enhance our knowledge of biological processes and systems. For instance, pathogens can manipulate the apoptotic process to either evade immune detection or to facilitate their spread. Furthermore, of particular clinical interest is the ability of cancer cells to evade apoptosis, hence allowing them to survive and proliferate uncontrollably. Thus, in this work, we propose a phase-field framework for simulating intrinsic or extrinsic apoptosis induced by an activation field, including deriving the configurational mechanics underlying such phenomena. Along with exploring varying conditions needed to initiate or reduce apoptosis, this can serve as a starting point for computational therapeutic testing. To showcase model capabilities, we present simulations exhibiting different types of cellular dynamics produced when varying the mechanisms underlying apoptosis. The model is subsequently applied to probe different morphological transitions, such as cell shrinkage, membrane blebbin
Tumor development is characterized by a compromised balance between cell life and death decision mechanisms, which are tighly regulated in normal cells. Understanding this process provides insights for developing new treatments for fighting with cancer. We present a study of a mathematical model describing cellular choice between survival and two alternative cell death modalities: apoptosis and necrosis. The model is implemented in discrete modeling formalism and allows to predict probabilities of having a particular cellular phenotype in response to engagement of cell death receptors. Using an original parameter sensitivity analysis developed for discrete dynamic systems, we determine the critical parameters affecting cellular fate decision variables that appear to be critical in the cellular fate decision and discuss how they are exploited by existing cancer therapies.
Regulation of cell proliferation is a crucial aspect of tissue development and homeostasis and plays a major role in morphogenesis, wound healing, and tumor invasion. A phenomenon of such regulation is contact inhibition, which describes the dramatic slowing of proliferation, cell migration and individual cell growth when multiple cells are in contact with each other. While many physiological, molecular and genetic factors are known, the mechanism of contact inhibition is still not fully understood. In particular, the relevance of cellular signaling due to interfacial contact for contact inhibition is still debated. Cellular automata (CA) have been employed in the past as numerically efficient mathematical models to study the dynamics of cell ensembles, but they are not suitable to explore the origins of contact inhibition as such agent-based models assume fixed cell sizes. We develop a minimal, data-driven model to simulate the dynamics of planar cell cultures by extending a probabilistic CA to incorporate size changes of individual cells during growth and cell division. We successfully apply this model to previous in-vitro experiments on contact inhibition in epithelial tissue: A
The diversity of cell populations is regulated by extracellular and intracellular variability. The latter includes genetic, epigenetic and stochastic variability, all contributing to the experimentally observed heterogeneity in response to external death-inducing stimuli. Studies of sources and regulation of variability in commitment to apoptotic cancer cell death are likely to identify the fundamental features of apoptotic protein networks that are responsible for determining the ultimate cell fate. Systems biology approaches, involving computer simulations of the biochemical reactions accompanied, if possible, by experimental verification of selected components of the model, are proving useful in determining the origins of cell-to-cell variability in response to external stress stimuli. Here we summarize our current understanding of the origins of stochastic variability in cells' commitment to apoptosis, and its implications in the field on cancer therapy.
It has been shown recently that changing the fluidic properties of a drug can improve its efficacy in ablating solid tumors. We develop a modeling framework for tumor ablation, and present the simplest possible model for drug diffusion in a spherical tumor with leaky boundaries and assuming cell death eventually leads to ablation of that cell effectively making the two quantities numerically equivalent. The death of a cell after a given exposure time depends on both the concentration of the drug and the amount of oxygen available to the cell. Higher oxygen availability leads to cell death at lower drug concentrations. It can be assumed that a minimum concentration is required for a cell to die, effectively connecting diffusion with efficacy. The concentration threshold decreases as exposure time increases, which allows us to compute dose-response curves. Furthermore, these curves can be plotted at much finer time intervals compared to that of experiments, which is used to produce a dose-threshold-response surface giving an observer a complete picture of the drug's efficacy for an individual. In addition, since the diffusion, leak coefficients, and the availability of oxygen is diff
The purpose of this paper is to implement a random death process into a persistent random walk model which produces subballistic superdiffusion (Lévy walk). We develop a Markovian model of cell motility with the extra residence variable $τ.$ The model involves a switching mechanism for cell velocity with dependence of switching rates on $τ$. This dependence generates intermediate subballistic superdiffusion. We derive master equations for the cell densities with the generalized switching terms involving the tempered fractional material derivatives. We show that the random death of cells has an important implication for the transport process through tempering of superdiffusive process. In the long-time limit we write stationary master equations in terms of exponentially truncated fractional derivatives in which the rate of death plays the role of tempering of a Lévy jump distribution. We find the upper and lower bounds for the stationary profiles corresponding to the ballistic transport and diffusion with the death rate dependent diffusion coefficient. Monte Carlo simulations confirm these bounds.
Lymphocyte populations, stimulated in vitro or in vivo, grow as cells divide. Stochastic models are appropriate because some cells undergo multiple rounds of division, some die, and others of the same type in the same conditions do not divide at all. If individual cells behave independently, each can be imagined as sampling from a probability density of times to division. The most convenient choice of density in mathematical and computational work, the exponential density, overestimates the probability of short division times. We consider a multi-stage model that produces an Erlang distribution of times to division, and an exponential distribution of times to die. The resulting dynamics of competing fates is a type of cyton model. Using Approximate Bayesian Computation, we compare our model to published cell counts, obtained after CFSE-labelled OT-I and F5 T cells were transferred to lymphopenic mice. The death rate is assumed to scale linearly with the generation (number of divisions) and the number of stages of undivided cells (generation $0$) is allowed to differ from that of cells that have divided at least once (generation greater than zero). Multiple stages are preferred in p
The transition from single-cell to multicellular behavior is important in early development but rarely studied. The starvation-induced aggregation of the social amoeba Dictyostelium discoideum into a multicellular slug is known to result from single-cell chemotaxis towards emitted pulses of cyclic adenosine monophosphate (cAMP). However, how exactly do transient short-range chemical gradients lead to coherent collective movement at a macroscopic scale? Here, we use a multiscale model verified by quantitative microscopy to describe wide-ranging behaviors from chemotaxis and excitability of individual cells to aggregation of thousands of cells. To better understand the mechanism of long-range cell-cell communication and hence aggregation, we analyze cell-cell correlations, showing evidence for self-organization at the onset of aggregation (as opposed to following a leader cell). Surprisingly, cell collectives, despite their finite size, show features of criticality known from phase transitions in physical systems. Application of external cAMP perturbations in our simulations near the sensitive critical point allows steering cells into early aggregation and towards certain locations b
In this work, the capability of primary, monocyte-derived dendritic cells (DCs) to uptake iron oxide magnetic nanoparticles (MNPs) is assessed and a strategy to induce selective cell death in these MNP-loaded DCs using external alternating magnetic fields (AMFs) is reported. No significant decrease in the cell viability of MNP-loaded DCs, compared to the control samples, was observed after five days of culture. The amount of MNPs incorporated into the cytoplasm was measured by magnetometry, which confirmed that 1 to 5 pg of the particles were uploaded per cell. The intracellular distribution of these MNPs, assessed by transmission electron microscopy, was found to be primarily inside the endosomic structures. These cells were then subjected to an AMF for 30 min, and the viability of the blank DCs (i.e., without MNPs), which were used as control samples, remained essentially unaffected. However, a remarkable decrease of viability from approximately 90% to 2-5% of DCs previously loaded with MNPs was observed after the same 30 min exposure to an AMF. The same results were obtained using MNPs having either positive (NH2+) or negative (COOH-) surface functional groups. In spite of the m
This paper demonstrates that the pattern of lipid spiculesthat emerge on the surface of red blood cells in the classic 'Discocyte to Echinocyte' shape change is a generative spiral, and presents a qualitative, fluid- driven mechanism for their production, compatible with the work of Douady and Couder. Implications for the dynamics of cell growth, plant cell phyllotaxy, programmed cell death and gravity sensitivity are explained in terms of a new qualitative model of cellular fluid dynamics.
We generated a computational approach to analyze the biomechanics of epithelial cell aggregates, either island or stripes or entire monolayers, that combines both vertex and contact-inhibition-of-locomotion models to include both cell-cell and cell-substrate adhesion. Examination of the distribution of cell protrusions (adhesion to the substrate) in the model predicted high order profiles of cell organization that agree with those previously seen experimentally. Cells acquired an asymmetric distribution of basal protrusions, traction forces and apical aspect ratios that decreased when moving from the edge to the island center. Our in silico analysis also showed that tension on cell-cell junctions and apical stress is not homogeneous across the island. Instead, these parameters are higher at the island center and scales up with island size, which we confirmed experimentally using laser ablation assays and immunofluorescence. Without formally being a 3-dimensional model, our approach has the minimal elements necessary to reproduce the distribution of cellular forces and mechanical crosstalk as well as distribution of principal stress in cells within epithelial cell aggregates. By mak
We study the dynamics of a thick polar epithelium subjected to the action of both an electric and a flow field in a planar geometry. We develop a generalized continuum hydrodynamic description and describe the tissue as a two component fluid system. The cells and the interstitial fluid are the two components and we keep all terms allowed by symmetry. In particular we keep track of the cell pumping activity for both solvent flow and electric current and discuss the corresponding orders of magnitude. We study the growth dynamics of tissue slabs, their steady states and obtain the dependence of the cell velocity, net cell division rate, and cell stress on the flow strength and the applied electric field. We find that finite thickness tissue slabs exist only in a restricted region of phase space and that relatively modest electric fields or imposed external flows can induce either proliferation or death.