Non-equilibrium phase transitions in hybrid Voronoi models of cell colonies
Pith reviewed 2026-05-23 17:38 UTC · model grok-4.3
The pith
Adding short-range repulsions for nuclear compressibility to the self-propelled Voronoi model produces multiple non-equilibrium phase transitions controlled by nucleus size.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the hybrid Voronoi model, steric repulsions representing nuclear compressibility compete with vertex interactions that mimic adhesion and cytoskeletal organization; the resulting non-equilibrium dynamics yield phase transitions from motility-induced phase separation through mesenchymal-like phases to disordered confluent states, with nucleus effective size providing an independent control parameter that crosses phase boundaries in agreement with experimental observations.
What carries the argument
Hybrid self-propelled Voronoi model with added short-range repulsive forces that encode nuclear compressibility.
If this is right
- Cells with larger or stiffer nuclei remain in confluent disordered states over wider ranges of motility and adhesion parameters.
- Lowering nuclear compressibility allows the system to enter mesenchymal-like phases even at moderate adhesion strengths.
- The same repulsion term suppresses motility-induced phase separation when nucleus size is increased.
- Phase diagrams in the space of repulsion strength, motility, and adhesion strength contain at least three distinct non-equilibrium regimes.
Where Pith is reading between the lines
- Nuclear mechanics could serve as a tunable handle in tissue engineering to steer collective cell migration without altering adhesion molecules.
- The model suggests that nuclear size variation within a single tissue could create spatial heterogeneity in cell packing and motility.
- If the repulsion term is made dependent on local cell density, the transitions might sharpen into first-order lines rather than crossovers.
Load-bearing premise
The short-range repulsive forces added to the Voronoi tessellation accurately capture the quantitative effects of nuclear compressibility without needing extra subcellular details or direct experimental calibration.
What would settle it
Measure the phase boundary locations in cell colonies while varying nuclear size or stiffness (for example via lamin mutations) and check whether the observed shifts match the locations predicted by changing the repulsion strength in the model.
Figures
read the original abstract
Eukaryotic cells are characterized by a stiff nucleus whose effect in modeling the collective behavior of cell aggregates is usually underestimated. However, increasing experimental evidence links nuclear modifications with phenotypic transition, like the one between epithelial and mesenchymal states. In this work, we explore the effect of short-range repulsive forces in the non-equilibrium dynamics of the self-propelled Voronoi model. We show that the competition between steric repulsions (representing nuclear/cellular compressibility) and Vertex interactions (mimicking cell-cell adhesion/interaction and cytoskeleton organization) generate a variety of non-equilibrium phase transitions from Motility-Induced Phase Separation to mesenchymal-like phases up to disordered confluent configurations. Notably, we found that tuning the nucleus's effective size/compressibility provides an additional way to cross the boundary between the different possible phases in line with experimental observations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a hybrid self-propelled Voronoi model augmented with short-range repulsive forces intended to represent nuclear/cellular compressibility. It claims that competition between these repulsions and standard vertex interactions produces a sequence of non-equilibrium phase transitions, from motility-induced phase separation through mesenchymal-like states to disordered confluent configurations, and that varying the effective nuclear size or compressibility provides an additional control parameter that crosses these boundaries in a manner consistent with experimental observations on nuclear modifications during epithelial-mesenchymal transitions.
Significance. If the mapping from the added repulsion to nuclear mechanics can be placed on a quantitative footing and the phase boundaries shown to be robust, the work would supply a computationally tractable way to incorporate nuclear stiffness into vertex-based models of cell colonies and would identify nuclear compressibility as a tunable parameter for collective transitions. The simulation framework itself appears straightforward to implement and could be useful for exploring parameter space once calibrated.
major comments (1)
- [Model definition and results (abstract and main text)] The central claim that tuning the nucleus's effective size/compressibility crosses phase boundaries 'in line with experimental observations' rests on an uncalibrated identification of the short-range repulsive term with nuclear compressibility. No derivation relating the repulsion amplitude or range to measured nuclear Young's modulus, lamin stiffness, or nuclear volume data is supplied, nor are simulated transition loci compared quantitatively to experimental parameter ranges for EMT-related nuclear changes. This mapping is load-bearing for the biological interpretation but remains an assumption without supporting evidence or sensitivity analysis.
minor comments (1)
- [Abstract] The abstract states that 'a variety of non-equilibrium phase transitions' are generated but provides no quantitative metrics (order parameters, transition loci, error bars, or system-size dependence) that would allow a reader to assess robustness.
Simulated Author's Rebuttal
We thank the referee for their constructive comments on our manuscript. We address the major concern regarding the mapping of the repulsive term to nuclear compressibility below, clarifying the phenomenological nature of the model while agreeing that additional discussion is needed.
read point-by-point responses
-
Referee: The central claim that tuning the nucleus's effective size/compressibility crosses phase boundaries 'in line with experimental observations' rests on an uncalibrated identification of the short-range repulsive term with nuclear compressibility. No derivation relating the repulsion amplitude or range to measured nuclear Young's modulus, lamin stiffness, or nuclear volume data is supplied, nor are simulated transition loci compared quantitatively to experimental parameter ranges for EMT-related nuclear changes. This mapping is load-bearing for the biological interpretation but remains an assumption without supporting evidence or sensitivity analysis.
Authors: We agree that the identification of the short-range repulsion with nuclear compressibility is phenomenological rather than derived from a first-principles mapping to Young's modulus or lamin stiffness. The term is introduced as a minimal representation of nuclear excluded volume and stiffness within the hybrid Voronoi framework, motivated by the established role of nuclear mechanics in cell aggregates. The central result is that varying the effective nuclear size and repulsion strength provides an additional control parameter that induces the reported sequence of non-equilibrium transitions; the statement that this occurs 'in line with experimental observations' refers to the qualitative direction of EMT-related nuclear changes (e.g., softening or volume alteration facilitating mesenchymal states), not a quantitative match. We have performed extensive parameter sweeps on repulsion amplitude and range that demonstrate robustness of the phase boundaries, but we acknowledge the absence of direct quantitative comparison to experimental nuclear parameter ranges. In the revised manuscript we will add an expanded discussion section that (i) explicitly states the phenomenological character of the mapping, (ii) includes sensitivity analysis on the repulsion parameters with the resulting phase diagrams, and (iii) places the chosen parameter values in the context of available literature values for nuclear stiffness, while noting the limitations and the need for future calibration. revision: partial
Circularity Check
No circularity; phases emerge from explicit simulation of added model terms
full rationale
The paper defines a hybrid self-propelled Voronoi model by adding short-range repulsive forces as a modeling choice to represent nuclear compressibility, then runs simulations to observe emergent non-equilibrium phases from competition with vertex interactions. No equations or results reduce a claimed prediction or phase boundary to a fitted parameter by construction, nor does any load-bearing step rely on a self-citation chain that itself assumes the target result. The work is a self-contained simulation study whose outputs are not forced by redefinition of inputs.
Axiom & Free-Parameter Ledger
free parameters (1)
- nuclear repulsion strength and range
axioms (2)
- domain assumption Voronoi tessellation plus vertex interactions adequately capture cell-cell adhesion and cytoskeleton effects
- domain assumption Short-range repulsions map directly to nuclear/cellular compressibility
Reference graph
Works this paper leans on
-
[1]
4: Nucleus size as a driver for the Epithelial-to-Mesenchymal transition in cell colonies
Geometrical force The geometrical force acting on the nucleus center of cell i is obtained deriving the total Vertex energy: ⃗F i v = − dEv d⃗ ri = − dEi v d⃗ ri − X j∈N (i) dEj v d⃗ ri (5) 7 /uni00A0=/uni00A00.05 /uni00A0=/uni00A00.25 /uni00A0=/uni00A00.5 /uni00A0=/uni00A00.75 /uni00A0=/uni00A00.85 /uni00A0=/uni00A00.9 /uni00A0=/uni00A01.0 (i) (ii) (iii)...
-
[2]
Repulsive force The repulsive force acting on the nucleus center of cell i is obtained by deriving the total repulsive energy ⃗F i r = − dϕr d⃗ ri = X j∈N (i) dϕj r d⃗ ri (22) and dΦj r drα i = −12 σ12 0 r13 ij ! (rα i − rα j ) rij (23) B. Observables As order parameters describing the global proper- ties of the system, we compute (i) the average long- ti...
work page 2019
-
[3]
Physical models of col- lective cell migration
Ricard Alert and Xavier Trepat. Physical models of col- lective cell migration. Annual Review of Condensed Mat- ter Physics , 11(1):77–101, 2020
work page 2020
-
[4]
Mesoscale physical prin- ciples of collective cell organization
Xavier Trepat and Erik Sahai. Mesoscale physical prin- ciples of collective cell organization. Nature Physics , 14(7):671–682, 2018
work page 2018
-
[5]
Physical mod- els of collective cell motility: from cell to tissue
Brian A Camley and Wouter-Jan Rappel. Physical mod- els of collective cell motility: from cell to tissue. Journal of physics D: Applied physics , 50(11):113002, 2017
work page 2017
-
[6]
Active vertex model for cell-resolution description of epithelial tissue mechanics
Daniel L Barton, Silke Henkes, Cornelis J Weijer, and Rastko Sknepnek. Active vertex model for cell-resolution description of epithelial tissue mechanics. PLoS compu- tational biology, 13(6):e1005569, 2017
work page 2017
-
[7]
Emergent structures and dynamics of cell colonies by con- tact inhibition of locomotion
Bart Smeets, Ricard Alert, Jir? Pesek, Ignacio Pag- onabarraga, Herman Ramon, and Romaric Vincent. Emergent structures and dynamics of cell colonies by con- tact inhibition of locomotion. Proceedings of the National Academy of Sciences , 113(51):14621–14626, 2016
work page 2016
-
[8]
Dense active matter model of motion patterns in confluent cell monolayers
Silke Henkes, Kaja Kostanjevec, J Martin Collinson, Rastko Sknepnek, and Eric Bertin. Dense active matter model of motion patterns in confluent cell monolayers. Nature communications, 11(1):1–9, 2020
work page 2020
-
[9]
Giant fluctuations and structural effects in a flocking epithelium
Fabio Giavazzi, Chiara Malinverno, Salvatore Corallino, Francesco Ginelli, Giorgio Scita, and Roberto Cerbino. Giant fluctuations and structural effects in a flocking epithelium. Journal of Physics D: Applied Physics , 50(38):384003, 2017
work page 2017
-
[10]
From motility-induced phase-separation to glassi- ness in dense active matter
Matteo Paoluzzi, Demian Levis, and Ignacio Pagonabar- raga. From motility-induced phase-separation to glassi- ness in dense active matter. Communications Physics , 5(1), May 2022
work page 2022
-
[11]
From flocking to glassiness in dense disordered polar active matter
Matteo Paoluzzi, Demian Levis, and Ignacio Pagonabar- raga. From flocking to glassiness in dense disordered polar active matter. Communications Physics , 7(1):57, 2024
work page 2024
-
[12]
Computational model for cell morphodynamics
Danying Shao, Wouter-Jan Rappel, and Herbert Levine. Computational model for cell morphodynamics. Phys. Rev. Lett., 105:108104, Sep 2010
work page 2010
-
[13]
D. Wenzel and A. Voigt. Multiphase field models for collective cell migration. Phys. Rev. E , 104:054410, Nov 2021
work page 2021
-
[14]
Inter- cellular friction and motility drive orientational order in cell monolayers
Michael Chiang, Austin Hopkins, Benjamin Loewe, M Cristina Marchetti, and Davide Marenduzzo. Inter- cellular friction and motility drive orientational order in cell monolayers. Proceedings of the National Academy of Sciences, 121(40):e2319310121, 2024
work page 2024
-
[15]
Entropy evalua- tion sheds light on ecosystem complexity
Mattia Miotto and Lorenzo Monacelli. Entropy evalua- tion sheds light on ecosystem complexity. Phys. Rev. E , 98:042402, Oct 2018
work page 2018
-
[16]
Austin Hopkins, Michael Chiang, Benjamin Loewe, Da- vide Marenduzzo, and M. Cristina Marchetti. Local yield and compliance in active cell monolayers. Phys. Rev. Lett., 129:148101, Sep 2022
work page 2022
-
[17]
Benjamin Loewe, Michael Chiang, Davide Marenduzzo, and M. Cristina Marchetti. Solid-liquid transition of de- formable and overlapping active particles. Phys. Rev. Lett., 125:038003, Jul 2020
work page 2020
-
[18]
Yeomans, and Amin Doost- mohammadi
Romain Mueller, Julia M. Yeomans, and Amin Doost- mohammadi. Emergence of active nematic behavior in monolayers of isotropic cells. Phys. Rev. Lett. , 122:048004, Feb 2019
work page 2019
-
[19]
Fran ¸ cois Graner and James A. Glazier. Simulation of biological cell sorting using a two-dimensional extended potts model. Phys. Rev. Lett. , 69:2013–2016, Sep 1992
work page 2013
-
[20]
Souvik Sadhukhan and Saroj Kumar Nandi. Theory and simulation for equilibrium glassy dynamics in cellular potts model of confluent biological tissue. Phys. Rev. E, 103:062403, Jun 2021. 10
work page 2021
-
[21]
R. Belousov, S. Savino, P. Moghe, T. Hiiragi, L. Rondoni, and A. Erzberger. Poissonian cellular potts models reveal nonequilibrium kinetics of cell sorting. Phys. Rev. Lett. , 132:248401, Jun 2024
work page 2024
-
[22]
Glazier and Fran ¸ cois Graner
James A. Glazier and Fran ¸ cois Graner. Simulation of the differential adhesion driven rearrangement of biological cells. Phys. Rev. E , 47:2128–2154, Mar 1993
work page 1993
-
[23]
A dynamic cell model for the formation of epithelial tissues
Tatsuzo Nagai and Hisao Honda. A dynamic cell model for the formation of epithelial tissues. Philosophical Mag- azine B , 81(7):699–719, 2001
work page 2001
-
[24]
A density-independent rigidity transition in bi- ological tissues
Dapeng Bi, JH Lopez, Jennifer M Schwarz, and M Lisa Manning. A density-independent rigidity transition in bi- ological tissues. Nature Physics, 11(12):1074–1079, 2015
work page 2015
-
[25]
Dapeng Bi, Xingbo Yang, M. Cristina Marchetti, and M. Lisa Manning. Motility-driven glass and jamming transitions in biological tissues. Phys. Rev. X , 6:021011, Apr 2016
work page 2016
-
[26]
Flocking transitions in confluent tissues
Fabio Giavazzi, Matteo Paoluzzi, Marta Macchi, Dapeng Bi, Giorgio Scita, M Lisa Manning, Roberto Cerbino, and M Cristina Marchetti. Flocking transitions in confluent tissues. Soft matter , 14(18):3471–3477, 2018
work page 2018
-
[27]
A geometrically controlled rigidity transition in a model for confluent 3d tissues
Matthias Merkel and M Lisa Manning. A geometrically controlled rigidity transition in a model for confluent 3d tissues. New Journal of Physics , 20(2):022002, 2018
work page 2018
-
[28]
A minimal-length approach unifies rigidity in underconstrained materials
Matthias Merkel, Karsten Baumgarten, Brian P Tighe, and M Lisa Manning. A minimal-length approach unifies rigidity in underconstrained materials. Proceedings of the National Academy of Sciences , 116(14):6560–6568, 2019
work page 2019
-
[29]
Effect of cellular rearrangement time delays on the rheology of vertex models for confluent tissues
Gonca Erdemci-Tandogan and M Lisa Manning. Effect of cellular rearrangement time delays on the rheology of vertex models for confluent tissues. PLoS computational biology, 17(6):e1009049, 2021
work page 2021
-
[30]
Genome hetero- geneity drives the evolution of species
Mattia Miotto and Lorenzo Monacelli. Genome hetero- geneity drives the evolution of species. Phys. Rev. Res. , 2:043026, Oct 2020
work page 2020
-
[31]
Cor- relating cell shape and cellular stress in motile confluent tissues
Xingbo Yang, Dapeng Bi, Michael Czajkowski, Matthias Merkel, M Lisa Manning, and M Cristina Marchetti. Cor- relating cell shape and cellular stress in motile confluent tissues. Proceedings of the National Academy of Sciences , 114(48):12663–12668, 2017
work page 2017
-
[32]
The influence of cell mechanics, cell-cell interactions, and proliferation on ep- ithelial packing
Reza Farhadifar, Jens-Christian R¨ oper, Benoit Aigouy, Suzanne Eaton, and Frank J¨ ulicher. The influence of cell mechanics, cell-cell interactions, and proliferation on ep- ithelial packing. Current biology, 17(24):2095–2104, 2007
work page 2095
-
[33]
Collective behav- ior and self-organization in neural rosette morphogenesis
Mattia Miotto, Maria Rosito, Matteo Paoluzzi, Valeria de Turris, Viola Folli, Marco Leonetti, Giancarlo Ruocco, Alessandro Rosa, and Giorgio Gosti. Collective behav- ior and self-organization in neural rosette morphogenesis. Frontiers in Cell and Developmental Biology , 11, August 2023
work page 2023
-
[34]
Collective cell durotaxis emerges from long-range intercellular force transmission
Raimon Sunyer, Vito Conte, Jorge Escribano, Alberto Elosegui-Artola, Anna Labernadie, L´ eo Valon, Daniel Navajas, Jos´ e Manuel Garc´ ıa-Aznar, Jos´ e J Mu˜ noz, Pere Roca-Cusachs, et al. Collective cell durotaxis emerges from long-range intercellular force transmission. Science, 353(6304):1157–1161, 2016
work page 2016
-
[35]
Tuning collective cell migration by cell?cell junction regulation
Peter Friedl and Roberto Mayor. Tuning collective cell migration by cell?cell junction regulation. Cold Spring Harbor perspectives in biology , 9 (4):a029199, 2017
work page 2017
-
[36]
Mattia Miotto and Lorenzo Monacelli. Tolomeo, a novel machine learning algorithm to measure information and order in correlated networks and predict their state. En- tropy, 23(9):1138, August 2021
work page 2021
-
[37]
Jin-Ah Park, Jae Hun Kim, Dapeng Bi, Jen- nifer A. Mitchel, Nader Taheri Qazvini, Kelan Tan- tisira, Chan Young Park, Maureen McGill, Sae-Hoon Kim, Bomi Gweon, Robert Steward Jr Jacob Notbohm, Stephanie Burger, Scott H. Randell, Alvin T. Kho, Dhananjay T. Tambe, Corey Hardin, Stephanie A. Shore, Elliot Israel, David A. Weitz, Daniel J. Tschumperlin, Eliza...
work page 2015
-
[38]
Pragya Arora, Souvik Sadhukhan, Saroj Kumar Nandi, Dapeng Bi, A. K. Sood, and Rajesh Ganapathy. A shape- driven reentrant jamming transition in confluent mono- layers of synthetic cell-mimics. Nature Communications, 15(1), July 2024
work page 2024
-
[39]
How motility drives the glassy dynamics in confluent epithelial monolayers? bioRxiv, March 2024
Souvik Sadhukhan, Manoj Kumar Nandi, Satyam Pandey, Matteo Paoluzzi, Chandan Dasgupta, Nir Gov, and Saroj Kumar Nandi. How motility drives the glassy dynamics in confluent epithelial monolayers? bioRxiv, March 2024
work page 2024
-
[40]
Chiara Malinverno, Salvatore Corallino, Fabio Gi- avazzi, Martin Bergert, Qingsen Li, Marco Leoni, Andrea Disanza, Emanuela Frittoli, Amanda Oldani, Emanuele Martini, Tobias Lendenmann, Gianluca De- florian, Galina V. Beznoussenko, Dimos Poulikakos, Kok Haur Ong, Marina Uroz, Xavier Trepat, Dario Parazzoli, Paolo Maiuri, Weimiao Yu, Aldo Ferrari, Roberto ...
work page 2017
-
[41]
A novel jamming phase diagram links tumor invasion to non-equilibrium phase separation
Wenying Kang, Jacopo Ferruzzi, Catalina-Paula Spatarelu, Yu Long Han, Yasha Sharma, Stephan A Koehler, Jennifer A Mitchel, Adil Khan, James P Butler, Darren Roblyer, et al. A novel jamming phase diagram links tumor invasion to non-equilibrium phase separation. Iscience, 24(11), 2021
work page 2021
-
[42]
Jennifer A. Mitchel, Amit Das, Michael J. O’Sullivan, Ian T. Stancil, Stephen J. DeCamp, Stephan Koehler, Oscar H. Oca˜ na, James P. Butler, Jeffrey J. Fredberg, M. Angela Nieto, Dapeng Bi, and Jin-Ah Park. In pri- mary airway epithelial cells, the unjamming transition is distinct from the epithelial-to-mesenchymal transition. Nature Communications, 11(1)...
work page 2020
-
[43]
Cell deformability drives fluid-to-fluid phase transition in active cell monolayers
Nen Saito and Shuji Ishihara. Cell deformability drives fluid-to-fluid phase transition in active cell monolayers. Science Advances, 10(19), May 2024
work page 2024
-
[44]
Dealing with the nucleus during cell migra- tion
Francisco J Calero-Cuenca, C´ atia S Janota, and Edgar R Gomes. Dealing with the nucleus during cell migra- tion. Current Opinion in Cell Biology , 50:35–41, Febru- ary 2018
work page 2018
-
[45]
Arthur Boutillon, Samhita P. Banavar, and Otger Camp` as. Conserved physical mechanisms of cell and tis- sue elongation. Development, 151(10), May 2024
work page 2024
-
[46]
Katarina Wolf, Mariska te Lindert, Marina Krause, Stephanie Alexander, Joost te Riet, Amanda L. Willis, Robert M. Hoffman, Carl G. Figdor, Stephen J. Weiss, and Peter Friedl. Physical limits of cell migration: Con- trol by ecm space and nuclear deformation and tuning by proteolysis and traction force. Journal of Cell Biology , 201(7):1069–1084, June 2013
work page 2013
-
[47]
Sussman, Fr´ ed´ eric Renner, Pablo Gottheil, Erik W
Steffen Grosser, J¨ urgen Lippoldt, Linda Oswald, Matthias Merkel, Daniel M. Sussman, Fr´ ed´ eric Renner, Pablo Gottheil, Erik W. Morawetz, Thomas Fuhs, Xiao- fan Xie, Steve Pawlizak, Anatol W. Fritsch, Benjamin 11 Wolf, Lars-Christian Horn, Susanne Briest, Bahriye Ak- tas, M. Lisa Manning, and Josef A. K¨ as. Cell and nu- cleus shape as an indicator of ...
work page 2021
-
[48]
Robert Chojowski, Ulrich S. Schwarz, and Falko Ziebert. The role of the nucleus for cell mechanics: an elastic phase field approach. Soft Matter , 20(22):4488–4503, 2024
work page 2024
-
[49]
J. Tailleur and M. E. Cates. Statistical mechanics of interacting run-and-tumble bacteria. Phys. Rev. Lett. , 100:218103, May 2008
work page 2008
-
[50]
Full phase diagram of active brownian disks: from melting to motility-induced phase separation
Pasquale Digregorio, Demian Levis, Antonio Suma, Leti- cia F Cugliandolo, Giuseppe Gonnella, and Ignacio Pag- onabarraga. Full phase diagram of active brownian disks: from melting to motility-induced phase separation. Phys- ical review letters , 121(9):098003, 2018
work page 2018
-
[51]
Confinement-induced transition between wavelike col- lective cell migration modes
Vanni Petrolli, Magali Le Goff, Monika Tadrous, Kirsten Martens, C´ edric Allier, Ondrej Mandula, Lionel Herv´ e, Silke Henkes, Rastko Sknepnek, Thomas Boudou, et al. Confinement-induced transition between wavelike col- lective cell migration modes. Physical review letters , 122(16):168101, 2019
work page 2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.