Cell stiffness is a key determinant of how cells deform, migrate, and adapt to mechanically restrictive environments, yet existing single-cell stiffness assays remain difficult to combine with molecular analysis and downstream functional studies. To address these limitations, we introduce a microfluidic platform, stiffness-based ferrohydrodynamic cell sorting (Stiff-FCS), designed for high-throughput quantification of single-cell stiffness, on-chip molecular analysis, and post-assay cell recovery. Stiff-FCS combines ferrofluid-driven actuation with graded confinement channels to control cell movement, induce deformation, and spatially separate cells based on stiffness. An inverse computational model converts cell position and morphology into quantitative Young's modulus values. We demonstrate stiffness profiling of hundreds to thousands of cells per chip within minutes, same-cell fluorescence-based protein analysis, and recovery of stiffness-defined cells for downstream assays. Across diverse human and mouse cell lines, Lamin A/C showed the most consistent association with stiffness, whereas softer cells exhibited greater migratory capacity than stiffer cells. In a series of human
VTX is an open-source molecular visualization software designed to overcome the scaling limitations of existing real-time molecular visualization software when handling massive molecular datasets. VTX employs a meshless molecular graphics engine utilizing impostor-based techniques and adaptive level-of-detail (LOD) rendering. This approach significantly reduces memory usage and enables real-time visualization and manipulation of large molecular systems. Performance benchmarks against VMD, PyMOL, and ChimeraX using a 114-million-bead Martini minimal whole-cell model demonstrate VTX's efficiency, maintaining consistent frame rates even under interactive manipulation on standard computer hardware. VTX incorporates features such as screen-space ambient occlusion (SSAO) for enhanced depth perception and free-fly navigation for intuitive exploration of large molecular systems. VTX is open-source and free for non commercial use. Binaries for Windows and Ubuntu Linux are available at http://vtx.drugdesign.fr. VTX source code is available at https://github.com/VTX-Molecular-Visualization.
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
Activation of naive CD8 T-cells can lead to the generation of multiple effector and memory subsets. Multiple parameters associated with activation conditions are involved in generating this diversity that is associated with heterogeneous molecular contents of activated cells. Although naive cell polarisation upon antigenic stimulation and the resulting asymmetric division are known to be a major source of heterogeneity and cell fate regulation, the consequences of stochastic uneven partitioning of molecular content upon subsequent divisions remain unclear yet. Here we aim at studying the impact of uneven partitioning on molecular-content heterogeneity and then on the immune response dynamics at the cellular level. To do so, we introduce a multiscale mathematical model of the CD8 T-cell immune response in the lymph node. In the model, cells are described as agents evolving and interacting in a 2D environment while a set of differential equations, embedded in each cell, models the regulation of intra and extracellular proteins involved in cell differentiation. Based on the analysis of in silico data at the single cell level, we 1 show that immune response dynamics can be explained by
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
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 division, aging, and stress recovery triggers spatial reorganization of cellular components in the cytoplasm, including membrane bound organelles, with molecular changes in their compositions and structures. However, it is not clear how these events are coordinated and how they integrate with regulation of molecular crowding. We use the budding yeast Saccharomyces cerevisiae as a model system to study these questions using recent progress in optical fluorescence microscopy and crowding sensing probe technology. We used a Förster Resonance Energy Transfer (FRET) based sensor, illuminated by confocal microscopy for high throughput analyses and Slimfield microscopy for single-molecule resolution, to quantify molecular crowding. We determine crowding in response to cellular growth of both mother and daughter cells, in addition to osmotic stress, and reveal hot spots of crowding across the bud neck in the burgeoning daughter cell. This crowding might be rationalized by the packing of inherited material, like the vacuole, from mother cells. We discuss recent advances in understanding the role of crowding in cellular regulation and key current challenges and conclude by presenting ou
The physical and chemical environment inside cells is of fundamental importance to all life but has traditionally been difficult to determine on a subcellular basis. Here we combine cutting-edge genomically integrated FRET biosensing to readout localized molecular crowding in single live yeast cells. Confocal microscopy allows us to build subcellular crowding heatmaps using ratiometric FRET, while whole-cell analysis demonstrates crowding is reduced when yeast is grown in elevated glucose concentrations. Simulations indicate that the cell membrane is largely inaccessible to these sensors and that cytosolic crowding is broadly uniform across each cell over a timescale of seconds. Millisecond single-molecule optical microscopy was used to track molecules and obtain brightness estimates that enabled calculation of crowding sensor copy numbers. The quantification of diffusing molecule trajectories paves the way for correlating subcellular processes and the physicochemical environment of cells under stress.
According to the National Cancer Institute, there were 9.5 million cancer-related deaths in 2018. A challenge in improving treatment is resistance in genetically unstable cells. The purpose of this study is to evaluate unsupervised machine learning on classifying treatment-resistant phenotypes in heterogeneous tumors through analysis of single cell RNA sequencing(scRNAseq) data with a pipeline and evaluation metrics. scRNAseq quantifies mRNA in cells and characterizes cell phenotypes. One scRNAseq dataset was analyzed (tumor/non-tumor cells of different molecular subtypes and patient identifications). The pipeline consisted of data filtering, dimensionality reduction with Principal Component Analysis, projection with Uniform Manifold Approximation and Projection, clustering with nine approaches (Ward, BIRCH, Gaussian Mixture Model, DBSCAN, Spectral, Affinity Propagation, Agglomerative Clustering, Mean Shift, and K-Means), and evaluation. Seven models divided tumor versus non-tumor cells and molecular subtype while six models classified different patient identification (13 of which were presented in the dataset); K-Means, Ward, and BIRCH often ranked highest with ~80% accuracy on th
The dynamic interplay between collective cell movement and the various molecules involved in the accompanying cell signalling mechanisms plays a crucial role in many biological processes including normal tissue development and pathological scenarios such as wound healing and cancer. Information about the various structures embedded within these processes allows a detailed exploration of the binding of molecular species to cell-surface receptors within the evolving cell population. In this paper we establish a general spatio-temporal-structural framework that enables the description of molecular binding to cell membranes coupled with the cell population dynamics. We first provide a general theoretical description for this approach and then illustrate it with two examples arising from cancer invasion.
A rigorous understanding of how multicellular behaviors arise from the actions of single cells requires quantitative frameworks that bridge the gap between genetic circuits, the arrangement of cells in space, and population-level behaviors. Here, we provide such a framework for a ubiquitous class of multicellular systems - namely, "secrete-and-sense cells" that communicate by secreting and sensing a signaling molecule. By using formal, mathematical arguments and introducing the concept of a phenotype diagram, we show how these cells tune their degrees of autonomous and collective behavior to realize distinct single-cell and population-level phenotypes; these phenomena have biological analogs, such as quorum sensing or paracrine signaling. We also define the "entropy of population," a measurement of the number of arrangements that a population of cells can assume, and demonstrate how a decrease in the entropy of population accompanies the formation of ordered spatial patterns. Our conceptual framework ties together diverse systems, including tissues and microbes, with common principles.
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
Recent findings show that single, non-neuronal cells are also able to learn signalling responses developing cellular memory. In cellular learning nodes of signalling networks strengthen their interactions e.g. by the conformational memory of intrinsically disordered proteins, protein translocation, miRNAs, lncRNAs, chromatin memory and signalling cascades. This can be described by a generalized, unicellular Hebbian learning process, where those signalling connections, which participate in learning, become stronger. Here we review those scenarios, where cellular signalling is not only repeated in a few times (when learning occurs), but becomes too frequent, too large, or too complex and overloads the cell. This leads to desensitisation of signalling networks by decoupling signalling components, receptor internalization, and consequent downregulation. These molecular processes are examples of anti-Hebbian learning and forgetting of signalling networks. Stress can be perceived as signalling overload inducing the desensitisation of signalling pathways. Aging occurs by the summative effects of cumulative stress downregulating signalling. We propose that cellular learning desensitisation
Myxococcus xanthus cells self-organize into aligned groups, clusters, at various stages of their lifecycle. Formation of these clusters is crucial for the complex dynamic multi-cellular behavior of these bacteria. However, the mechanism underlying the cell alignment and clustering is not fully understood. Motivated by studies of clustering in self-propelled rods, we hypothesized that M. xanthus cells can align and form clusters through pure mechanical interactions among cells and between cells and substrate. We test this hypothesis using an agent-based simulation framework in which each agent is based on the biophysical model of an individual M. xanthus cell. We show that model agents, under realistic cell flexibility values, can align and form cell clusters but only when periodic reversals of cell directions are suppressed. However, by extending our model to introduce the observed ability of cells to deposit and follow slime trails, we show that effective trail-following leads to clusters in reversing cells. Furthermore, we conclude that mechanical cell alignment combined with slime-trail-following is sufficient to explain the distinct clustering behaviors observed for wild-type a
We define a model for random (abstract) cell complexes (CCs), similiar to the well-known Erdős-Rényi model for graphs and its extensions for simplicial complexes. To build a random cell complex, we first draw from an Erdős-Rényi graph, and consecutively augment the graph with cells for each dimension with a specified probability. As the number of possible cells increases combinatorially -- e.g., 2-cells can be represented as cycles, or permutations -- we derive an approximate sampling algorithm for this model limited to two-dimensional abstract cell complexes. As a basis for this algorithm, we first introduce a spanning-tree-based method that samples simple cycles and allows the efficient approximation of various properties, most notably the probability of occurence of a given cycle. This approximation is of independent interest as it enables the approximation of a wide variety of cycle-related graph statistics using importance sampling. We use this to approximate the number of cycles of a given length on a graph, allowing us to calculate the sampling probability to arrive at a desired expected number of sampled 2-cells. The probability approximation also trivially leads to a sampl
Understanding mechanosensitivity, i.e. how cells sense the stiffness of their environment is very important, yet there is a fundamental difficulty in understanding its mechanism: to measure an elastic modulus one requires two points of application of force - a measuring and a reference point. The cell in contact with substrate has only one (adhesion) point to work with, and thus a new method of measurement needs to be invented. The aim of this theoretical work is to develop a self-consistent physical model for mechanosensitivity, a process by which a cell detects the mechanical stiffness of its environment (e.g. a substrate it is attached to via adhesion points) and generates an appropriate chemical signaling to remodel itself in response to this environment. The model uses the molecular mechanosensing complex of latent TGF-$β$ attached to the adhesion point as the biomarker. We show that the underlying Brownian motion in the substrate is the reference element in the measuring process. The model produces the closed expression for the rate of release of active TGF-$β$, which depends on the substrate stiffness and the pulling force coming from the cell in a subtle and non-trivial way
Most microorganisms regulate their cell size. We review here some of the mathematical formulations of the problem of cell size regulation. We focus on coarse-grained stochastic models and the statistics they generate. We review the biologically relevant insights obtained from these models. We then describe cell cycle regulation and their molecular implementations, protein number regulation, and population growth, all in relation to size regulation. Finally, we discuss several future directions for developing understanding beyond phenomenological models of cell size regulation.
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
Organisms across all domains of life regulate the size of their cells. However, the means by which this is done is poorly understood. We study two abstracted "molecular" models for size regulation: inhibitor dilution and initiator accumulation. We apply the models to two settings: bacteria like Escherichia coli, that grow fully before they set a division plane and divide into two equally sized cells, and cells that form a bud early in the cell division cycle, confine new growth to that bud, and divide at the connection between that bud and the mother cell, like the budding yeast Saccharomyces cerevisiae. In budding cells, delaying cell division until buds reach the same size as their mother leads to very weak size control, with average cell size and standard deviation of cell size increasing over time and saturating up to 100-fold higher than those values for cells that divide when the bud is still substantially smaller than its mother. In budding yeast, both inhibitor dilution or initiator accumulation models are consistent with the observation that the daughters of diploid cells add a constant volume before they divide. This adder behavior has also been observed in bacteria. We f