Collective behaviors in cellular systems are regulated not only by biochemical signalling pathways but also by intercellular mechanical forces, whose quantification in contractile monolayers remains poorly understood. Here, by integrating traction force microscopy and numerical simulations, we reconstruct the stress distribution in C2C12 myoblast monolayers to reveal the roles of local mechanical forces in determining the collective cellular structures. We find that contractile monolayers exhibit positive maximum and negative minimum principal stresses, reflecting the intrinsic anisotropy of active tension. Distinct stress patterns emerge around topological defects, coinciding with singularities in cell alignment, density, and morphology, indicating a strong coupling between mechanical forces and structural organization. Moreover, tensile stresses are preferentially transmitted along the cell elongation axis and compressive stresses transversely, demonstrating that local stress guides cell arrangement. This mechanical guidance appears to be universal among contractile systems, as observed also in bone marrow-derived mesenchymal stem cells. Together, our work establishes a quantitat
The mechanics of animal cells is strongly determined by stress fibers, which are contractile filament bundles that form dynamically in response to extracellular cues. Stress fibers allow the cell to adapt its mechanics to environmental conditions and to protect it from structural damage. While the physical description of single stress fibers is well-developed, much less is known about their spatial distribution on the level of whole cells. Here, we combine a finite element method for one-dimensional fibers embedded in an elastic bulk medium with dynamical rules for stress fiber formation based on genetic algorithms. We postulate that their main goal is to achieve minimal mechanical stress in the bulk material with as few fibers as possible. The fiber positions and configurations resulting from this optimization task alone are in good agreement with those found in experiments where cells in 3D-scaffolds were mechanically strained at one attachment point. For optimized configurations, we find that stress fibers typically run through the cell in a diagonal fashion, similar to reinforcement strategies used for composite material.
To understand cell migration, it is crucial to gain knowledge on how cells exert and integrate forces on/from their environment. A quantity of prime interest for biophysicists interested in cell movements modeling is the intracellular stresses. Up to now, three different methods have been proposed to calculate it, they are all in the regime of the thin plate approximation. Two are based on solving the mechanical equilibrium equation inside the cell material (Monolayer Stress Microscopy, and Bayesian Inference Stress Microscopy) and one is based on the continuity of displacement at the cell/substrate interface (Intracellular Stress Microscopy). We show here using 3D FEM modeling that these techniques do not calculate the same quantities (as was previously assumed), the first techniques calculate the sum of the active and resistive stresses within the cell, whereas the last one only calculate the resistive component. Combining these techniques should in principle permit to get access to the active stress alone.
We investigate a continuum mechanical model for an adherent cell on two dimensional adhesive micropatterned substrates. The cell is modeled as an isotropic and homogeneous elastic material subject to uniform internal contractile stresses. The build-up of tension from cortical actin bundles at the cell periphery is incorporated by introducing an energy cost for bending of the cell boundary, resulting to a resistance to changes in local curvature. Integrin-based adhesions are modeled as harmonic springs, that pin the cell to adhesive patches of a predefined geometry. Using Monte Carlo simulations and analytical techniques we investigate the competing effects of bulk contractility and cortical bending rigidity in regulating cell shapes on non-adherent regions. We show that the crossover from convex to concave cell edges is controlled by the interplay between contractile stresses and boundary bending rigidity. In particular, the cell boundary becomes concave beyond a critical value of the contractile stress that is proportional to the cortical bending rigidity. Furthermore, the intracellular stresses are found largely concentrated at the concave edge of the cell. The model can be used
We introduce AFSTRESS, the first multi-label corpus of self-reported stress narratives in Dari (Eastern Persian), comprising 737 responses collected from Afghan individuals during an ongoing humanitarian crisis. Participants describe experienced stress and select emotion and stressor labels via Dari checklists. The dataset enables analysis at three levels: computational (multi-label classification), social (structural drivers and gender disparities), and psychological (learned helplessness, chronic stress, and emotional cascade patterns). It includes 12 binary labels (5 emotions, 7 stressors), with high label cardinality (5.54) and density (0.462), reflecting complex, multi-dimensional stress. Structural stressors dominate: uncertain future (62.6 percent) and education closure (60.0 percent) exceed emotional states, indicating stress is primarily structurally driven. The strongest co-occurrence is between hopelessness and uncertain future (J = 0.388). Baseline experiments show that character TF-IDF with Linear SVM achieves Micro-F1 = 0.663 and Macro-F1 = 0.651, outperforming ParsBERT and XLM-RoBERTa, while threshold tuning improves Micro-F1 by 10.3 points. AFSTRESS provides the fir
Survival of living tumor cells underlies many influences such as nutrient saturation, oxygen level, drug concentrations or mechanical forces. Data-supported mathematical modeling can be a powerful tool to get a better understanding of cell behavior in different settings. However, under consideration of numerous environmental factors mathematical modeling can get challenging. We present an approach to model the separate influences of each environmental quantity on the cells in a collective manner by introducing the "environmental stress level". It is an immeasurable auxiliary variable, which quantifies to what extent viable cells would get in a stressed state, if exposed to certain conditions. A high stress level can inhibit cell growth, promote cell death and influence cell movement. As a proof of concept, we compare two systems of ordinary differential equations, which model tumor cell dynamics under various nutrient saturations respectively with and without considering an environmental stress level. Particle-based Bayesian inversion methods are used to quantify uncertainties and calibrate unknown model parameters with time resolved measurements of in vitro populations of liver ca
Collective cell migration is a highly regulated process involved in wound healing, cancer metastasis and morphogenesis. Mechanical interactions among cells provide an important regulatory mechanism to coordinate such collective motion. Using a Self-Propelled Voronoi (SPV) model that links cell mechanics to cell shape and cell motility, we formulate a generalized mechanical inference method to obtain the spatio-temporal distribution of cellular stresses from measured traction forces in motile tissues and show that such traction-based stresses match those calculated from instantaneous cell shapes. We additionally use stress information to characterize the rheological properties of the tissue. We identify a motility-induced swim stress that adds to the interaction stress to determine the global contractility or extensibility of epithelia. We further show that the temporal correlation of the interaction shear stress determines an effective viscosity of the tissue that diverges at the liquid-solid transition, suggesting the possibility of extracting rheological information directly from traction data.
Many morphogenetic processes involve mechanical rearrangement of epithelial tissues that is driven by precisely regulated cytoskeletal forces and cell adhesion. The mechanical state of the cell and intercellular adhesion are not only the targets of regulation, but are themselves likely signals that coordinate developmental process. Yet, because it is difficult to directly measure mechanical stress {\it in vivo} on sub-cellular scale, little is understood about the role of mechanics of development. Here we present an alternative approach which takes advantage of the recent progress in live imaging of morphogenetic processes and uses computational analysis of high resolution images of epithelial tissues to infer relative magnitude of forces acting within and between cells. We model intracellular stress in terms of bulk pressure and interfacial tension, allowing these parameters to vary from cell to cell and from interface to interface. Assuming that epithelial cell layers are close to mechanical equilibrium, we use the observed geometry of the two dimensional cell array to infer interfacial tensions and intracellular pressures. Here we present the mathematical formulation of the prop
In Escherichia coli the response to DNA damage shows strong cell-to-cell-heterogenity. This results in a random delay in cell division and asymmetrical binary fission of single cells, which can compromise the size homeostasis of the population. To quantify the effect of the heterogeneous response to genotoxic stress (called SOS response in E. coli) on the growth of the bacterial population, we propose a flexible time-continuous parametric model of individual-based population dynamics. We construct a stochastic model based on the "adder" size-control mechanism, extended to incorporate the dynamics of the SOS response and its effect on cell division. The model is fitted to individual lineage data obtained in a 'mother machine' microfluidic device. We show that the heterogeneity of the SOS response can bias the observed division rate. In particular, we show that the adder division rate is decreased by SOS induction and that this perturbative effect is stronger in fast-growing conditions.
The warming climate is causing livestock to experience heat stress at an increasing frequency. Holstein cows are particularly susceptible to heat stress because of their high metabolic rate. Heat stress negatively affects immune function, particularly with respect to the cell-mediated immune response, which leads to increased susceptibility to disease. Cattle identified as having enhanced immune response have lower incidence of disease. Therefore, the objective of this study was to evaluate the impact of in vitro heat challenge on blood mononuclear cells from dairy cattle, that had previously been ranked for immune response, in terms of heat shock protein 70 concentration, nitric oxide production, and cell proliferation. Bovine blood mononuclear cells, from Holstein dairy cattle previously ranked for immune response based on their estimated breeding values, were subjected to three heat treatments: thermoneutral, heat stress 1 and heat stress 2. Cells of each treatment were evaluated for heat shock protein 70, cell proliferation and nitric oxide production. Blood mononuclear cells from dairy cattle classified as high immune responders, based on their estimated breeding values for an
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
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
Cellular mechanics plays an important role in epithelial morphogenesis, a process wherein cells reshape and rearrange to produce tissue-scale deformations. However, the study of tissue-scale mechanics is impaired by the difficulty of direct measurement of stress in-vivo. Alternative, image-based inference schemes aim to estimate stress from snapshots of cellular geometry but are challenged by sensitivity to fluctuations and measurement noise as well as the dependence on boundary conditions. Here we overcome these difficulties by introducing a new variational approach - the Geometrical Variation Method (GVM) - which exploits the fundamental duality between stress and cellular geometry that exists in the state of mechanical equilibrium of discrete mechanical networks that approximate cellular tissues. In the Geometrical Variation Method, the two dimensional apical geometry of an epithelial tissue is approximated by a 2D tiling with Circular Arc Polygons (CAP) in which the arcs represent intercellular interfaces defined by the balance of local line tension and pressure differentials between adjacent cells. We take advantage of local constraints that mechanical equilibrium imposes on C
Using a popular vertex-based model to describe a spatially disordered planar epithelial monolayer, we examine the relationship between cell shape and mechanical stress at the cell and tissue level. Deriving expressions for stress tensors starting from an energetic formulation of the model, we show that the principal axes of stress for an individual cell align with the principal axes of shape, and we determine the bulk effective tissue pressure when the monolayer is isotropic at the tissue level. Using simulations for a monolayer that is not under peripheral stress, we fit parameters of the model to experimental data for Xenopus embryonic tissue. The model predicts that mechanical interactions can generate mesoscopic patterns within the monolayer that exhibit long-range correlations in cell shape. The model also suggests that the orientation of mechanical and geometric cues for processes such as cell division are likely to be strongly correlated in real epithelia. Some limitations of the model in capturing geometric features of Xenopus epithelial cells are highlighted.
The motility of a cell can be triggered or inhibited not only by an applied force but also by a mechanically neutral force couple. This type of loading, represented by an applied stress and commonly interpreted as either squeezing or stretching, can originate from extrinsic interaction of a cell with its neighbors. To quantify the effect of applied stresses on cell motility we use an analytically transparent one-dimensional model accounting for active myosin contraction and induced actin turnover. We show that stretching can polarize static cells and initiate cell motility while squeezing can symmetrize and arrest moving cells. We show further that sufficiently strong squeezing can lead to the loss of cell integrity. The overall behavior of the system depends on the two dimensionless parameters characterizing internal driving (chemical activity) and external loading (applied stress). We construct a phase diagram in this parameter space distinguishing between, static, motile and collapsed states. The obtained results are relevant for the mechanical understanding of contact inhibition and the epithelial-to-mesenchymal transition.
Collective cell migration in cohesive units is vital for tissue morphogenesis, wound repair, and immune response. While the fundamental driving forces for collective cell motion stem from contractile and protrusive activities of individual cells, it remains unknown how their balance is optimized to maintain tissue cohesiveness and the fluidity for motion. Here we present a cell-based computational model for collective cell migration during wound healing that incorporates mechanochemical coupling of cell motion and adhesion kinetics with stochastic transformation of active motility forces. We show that a balance of protrusive motility and actomyosin contractility is optimized for accelerating the rate of wound repair, which is robust to variations in cell and substrate mechanical properties. This balance underlies rapid collective cell motion during wound healing, resulting from a tradeoff between tension mediated collective cell guidance and active stress relaxation in the tissue.
The complexity of the cells can be described and understood by a number of networks such as protein-protein interaction, cytoskeletal, organelle, signalling, gene transcription and metabolic networks. All these networks are highly dynamic producing continuous rearrangements in their links, hubs, network-skeleton and modules. Here we describe the adaptation of cellular networks after various forms of stress causing perturbations, congestions and network damage. Chronic stress decreases link-density, decouples or even quarantines modules, and induces an increased competition between network hubs and bridges. Extremely long or strong stress may induce a topological phase transition in the respective cellular networks, which switches the cell to a completely different mode of cellular function. We summarize our initial knowledge on network restoration after stress including the role of molecular chaperones in this process. Finally, we discuss the implications of stress-induced network rearrangements in diseases and ageing, and propose therapeutic approaches both to increase the robustness and help the repair of cellular networks.
Reduced blood flow, as occurring in ischemia or resulting from exposure to microgravity such as encountered in space flights, induces a decrease in the level of shear stress sensed by the endothelial cells forming the inner part of blood vessels. In the present study, we use a microvasculature-on-a-chip device in order to investigate in vitro the effect of such a reduction in shear stress on shear-adapted endothelial cells. We find that, within one hour of exposition to reduced wall shear stress, human umbilical vein endothelial cells undergo a reorganization of their actin skeleton, with a decrease in the number of stress fibers and actin being recruited into the cells' peripheral band, indicating a fairly fast change in cells' phenotype due to altered flow.
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