How Cells Stay Together: A Mechanism for Maintenance of a Robust Cluster Explored by Local and Non-local Continuum Models.

Formation of organs and specialized tissues in embryonic development requires migration of cells to specific targets. In some instances, such cells migrate as a robust cluster. We here explore a recent local approximation of non-local continuum models by Falcó et al. (SIAM J Appl Math 84:17-42, 2023). We apply their theoretical results by specifying biologically-based cell-cell interactions, showing how such cell communication results in an effective attraction-repulsion Morse potential. We then explore the clustering instability, the existence and size of the cluster, and its stability. For attractant-repellent chemotaxis, we derive an explicit condition on cell and chemical properties that guarantee the existence of robust clusters. We also extend their work by investigating the accuracy of the local approximation relative to the full non-local model.

Cell Repolarization: A Bifurcation Study of Spatio-Temporal Perturbations of Polar Cells.

The intrinsic polarity of migrating cells is regulated by spatial distributions of protein activity. Those proteins (Rho-family GTPases, such as Rac and Rho) redistribute in response to stimuli, determining the cell front and back. Reaction-diffusion equations with mass conservation and positive feedback have been used to explain initial polarization of a cell. However, the sensitivity of a polar cell to a reversal stimulus has not yet been fully understood. We carry out a PDE bifurcation analysis of two polarity models to investigate routes to repolarization: (1) a single-GTPase ("wave-pinning") model and (2) a mutually antagonistic Rac-Rho model. We find distinct routes to reversal in (1) vs. (2). We show numerical simulations of full PDE solutions for the RD equations, demonstrating agreement with predictions of the bifurcation results. Finally, we show that simulations of the polarity models in deforming 1D model cells are consistent with biological experiments.

Bridging from single to collective cell migration: A review of models and links to experiments.

Mathematical and computational models can assist in gaining an understanding of cell behavior at many levels of organization. Here, we review models in the literature that focus on eukaryotic cell motility at 3 size scales: intracellular signaling that regulates cell shape and movement, single cell motility, and collective cell behavior from a few cells to tissues. We survey recent literature to summarize distinct computational methods (phase-field, polygonal, Cellular Potts, and spherical cells). We discuss models that bridge between levels of organization, and describe levels of detail, both biochemical and geometric, included in the models. We also highlight links between models and experiments. We find that models that span the 3 levels are still in the minority.

Cell Size, Mechanical Tension, and GTPase Signaling in the Single Cell.

Cell polarization requires redistribution of specific proteins to the nascent front and back of a eukaryotic cell. Among these proteins are Rac and Rho, members of the small GTPase family that regulate the actin cytoskeleton. Rac promotes actin assembly and protrusion of the front edge, whereas Rho activates myosin-driven contraction at the back. Mathematical models of cell polarization at many levels of detail have appeared. One of the simplest based on "wave-pinning" consists of a pair of reaction-diffusion equations for a single GTPase. Mathematical analysis of wave-pinning so far is largely restricted to static domains in one spatial dimension. Here, we extend the analysis to cells that change in size, showing that both shrinking and growing cells can lose polarity. We further consider the feedback between mechanical tension, GTPase activation, and cell deformation in both static, growing, shrinking, and moving cells. Special cases (spatially uniform cell chemistry, the absence or presence of mechanical feedback) are analyzed, and the full model is explored by simulations in 1D. We find a variety of novel behaviors, including "dilution-induced" oscillations of Rac activity and cell size, as well as gain or loss of polarization and motility in the model cell.

Correlated random walks inside a cell: actin branching and microtubule dynamics.

Correlated random walks (CRW) have been explored in many settings, most notably in the motion of individuals in a swarm or flock. But some subcellular systems such as growth or disassembly of bio-polymers can also be described with similar models and understood using related mathematical methods. Here we consider two examples of growing cytoskeletal elements, actin and microtubules. We use CRW or generalized CRW-like PDEs to model their spatial distributions. In each case, the linear models can be reduced to a Telegrapher's equation. A combination of explicit solutions (in one case) and numerical solutions (in the other) demonstrates that the approach to steady state can be accompanied by (decaying) waves.

A space-jump derivation for non-local models of cell-cell adhesion and non-local chemotaxis.

Cellular adhesion provides one of the fundamental forms of biological interaction between cells and their surroundings, yet the continuum modelling of cellular adhesion has remained mathematically challenging. In 2006, Armstrong et al. proposed a mathematical model in the form of an integro-partial differential equation. Although successful in applications, a derivation from an underlying stochastic random walk has remained elusive. In this work we develop a framework by which non-local models can be derived from a space-jump process. We show how the notions of motility and a cell polarization vector can be naturally included. With this derivation we are able to include microscopic biological properties into the model. We show that particular choices yield the original Armstrong model, while others lead to more general models, including a doubly non-local adhesion model and non-local chemotaxis models. Finally, we use random walk simulations to confirm that the corresponding continuum model represents the mean field behaviour of the stochastic random walk.