共找到 20 条结果
The strategic induction of cell death serves as a crucial immune defense mechanism for the eradication of pathogenic infections within host cells. Investigating the molecular mechanisms underlying immunogenic cell pathways has significantly enhanced our understanding of the host's immunity. This review provides a comprehensive overview of the immunogenic cell death mechanisms triggered by pathogen infections, focusing on the critical role of pattern recognition receptors. In response to infections, host cells dictate a variety of cell death pathways, including apoptosis, pyroptosis, necrosis, and lysosomal cell death, which are essential for amplifying immune responses and controlling pathogen dissemination. Key components of these mechanisms are host cellular receptors that recognize pathogen-associated ligands. These receptors activate downstream signaling cascades, leading to the expression of immunoregulatory genes and the production of antimicrobial cytokines and chemokines. Particularly, the inflammasome, a multi-protein complex, plays a pivotal role in these responses by processing pro-inflammatory cytokines and inducing pyroptotic cell death. Pathogens, in turn, have evolve
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.
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
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
Cell type annotation is a key task in analyzing the heterogeneity of single-cell RNA sequencing data. Although recent foundation models automate this process, they typically annotate cells independently, without considering batch-level cellular context or providing explanatory reasoning. In contrast, human experts often annotate distinct cell types for different cell clusters based on their domain knowledge. To mimic this workflow, we introduce the CellPuzzles task, where the objective is to assign unique cell types to a batch of cells. This benchmark spans diverse tissues, diseases, and donor conditions, and requires reasoning across the batch-level cellular context to ensure label uniqueness. We find that off-the-shelf large language models (LLMs) struggle on CellPuzzles, with the best baseline (OpenAI's o1) achieving only 19.0% batch-level accuracy. To fill this gap, we propose Cell-o1, a 7B LLM trained via supervised fine-tuning on distilled reasoning traces, followed by reinforcement learning with batch-level rewards. Cell-o1 achieves state-of-the-art performance, outperforming o1 by over 73% and generalizing well across contexts. Further analysis of training dynamics and reas
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 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
The target theory is the most classical hypothesis explaining radiation-induced cell death, the physical or biological nature of the "target" remains ambiguous. This study hypothesizes that the distribution of DNA double-strand breaks (DSBs) within the 3D genome is a pivotal factor affecting the probability of radiation-induced cell death. We propose that clustered DSBs in DNA segments with high interaction frequencies are more susceptible to leading to cell death than isolated DSBs. Topologically associating domains (TAD) can be regarded as the reference unit for evaluating the impact of DSB clustering in the 3D genome. To quantify this correlation between the DSB distribution in 3D genome and radiation-induced effect, we developed a simplified model considering the DSB distribution across TADs. Utilizing track-structure Monte Carlo codes to simulate the electron and carbon ion irradiation, we calculated the incidence of each case across a variety of radiation doses and LETs. Our simulation results indicate that DSBs in TADs with frequent interactions (case 3) are significantly more likely to induce cell death than clustered DSBs within a single TAD (case 2). Moreover, case 2 is s
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.
We introduce a general class of branching Markov processes for the modelling of a parasite infection in a cell population. Each cell contains a quantity of parasites which evolves as a diffusion with positive jumps. The growth rate, diffusive function and positive jump rate of this quantity of parasites depend on its current value. The division rate of the cells also depends on the quantity of parasites they contain. At division, a cell gives birth to two daughter cells and shares its parasites between them. Cells may also die, at a rate which may depend on the quantity of parasites they contain. We study the long time behaviour of the parasite infection. In particular, we are interested in the quantity of parasites in a `typical' cell and on the survival of the cell population. We specifically focus on the influence of two parameters on the probability for the cell population to survive and/or contain the parasite infection: the law of the sharing of the parasites between the daughter cells at division and the form of the division and death rates of the cells as functions of the quantity of parasites they contain.
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
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
We investigate the effects of cell division and apopotosis on collective dynamics in two-dimensional epithelial tissues. Our model includes three key ingredients observed across many epithelia, namely cell-cell adhesion, cell death and a cell division process that depends on the surrounding environment. We show a rich non-equilibrium phase diagram depending on the ratio of cell death to cell division and on the adhesion strength. For large apopotosis rates, cells die out and the tissue disintegrates. As the death rate decreases, however, we show, consecutively, the existence of a gas-like phase, a gel-like phase, and a dense confluent (tissue) phase. Most striking is the observation that the tissue is self-melting through its own internal activity, ruling out the existence of any glassy phase.
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.
Apoptosis, or genetically programmed cell death, is a crucial cellular process that maintains the balance between life and death in cells. The precise molecular mechanism of apoptosis signaling and how these two pathways are differentially activated under distinct apoptotic stimuli is poorly understood. We developed a Monte Carlo-based stochastic simulation model that can characterize distinct signaling behaviors in the two major pathways of apoptotic signaling using a novel probability distribution-based approach. Specifically, we show that for a weak death signal, such as low levels of death ligand Fas (CD95) binding or under stress conditions, the type 2 mitochondrial pathway dominates apoptotic signaling. Our results also show signaling in the type 2 pathway is stochastic, where the population average over many cells does not capture the cell-to-cell fluctuations in the time course (~1 - 10 hours) of downstream caspase-3 activation. On the contrary, the probability distribution of caspase-3 activation for the mitochondrial pathway shows a distinct bimodal behavior that can be used to characterize the stochastic signaling in type 2 apoptosis. Interestingly, such stochastic fluct
Cancer stem cells are controlled by developmental networks that are often topologically indistinguishable from normal, healthy stem cells. The question is why cancer stem cells can be both phenotypically distinct and have morphological effects so different from normal stem cells. The difference between cancer stem cells and normal stem cells lies not in differences their network architecture, but rather in the spatial-temporal locality of their activation in the genome and the resulting expression in the body. The metastatic potential cancer stem cells is not based primarily on their network divergence from normal stem cells, but on non-network based genetic changes that enable the evolution of gene-based phenotypic properties of the cell that permit its escape and travel to other parts of the body. Stem cell network theory allows the precise prediction of stem cell behavioral dynamics and a mathematical description of stem cell proliferation for both normal and cancer stem cells. It indicates that the best therapeutic approach is to tackle the highest order stem cells first, otherwise spontaneous remission of so called cured cancers will always be a danger. Stem cell networks poin
For years, extensive efforts have been made to discover effective, non-invasive anti-cancer therapies. Cold atmospheric plasma (CAP), is a near room temperature ionized gas composed of reactive species, charged particles, neutral particles, and electrons. CAP also has several physical factors including thermal radiation, ultraviolet (UV) radiation, and electromagnetic (EM) waves. Most of the previously reported biological effects of CAP have relied on direct contact between bulk plasma and cells, resulting in the chemical effects generally seen after CAP treatment. In this paper, we demonstrate that the electromagnetic emission produced by CAP can lead to the death of B16F10 melanoma cancer cells via a transbarrier contactless method. When compared with the effect of reactive species, the effect of the physical factors causes much greater growth inhibition. The physical-triggered growth inhibition is due to a new type of cell death, characterized by rapid leakage of bulk water from the cells, resulting in bubbles on the cell membrane, and cytoplasm shrinkage. The results of this study introduce a new possible mechanism of CAP induced cancer cell death and build a foundation for CAP