A minimal physical model for curvotaxis driven by curved protein complexes at the cell's leading edge.

Cells often migrate on curved surfaces inside the body, such as curved tissues, blood vessels, or highly curved protrusions of other cells. Recent in vitro experiments provide clear evidence that motile cells are affected by the curvature of the substrate on which they migrate, preferring certain curvatures to others, termed "curvotaxis." The origin and underlying mechanism that gives rise to this curvature sensitivity are not well understood. Here, we employ a "minimal cell" model which is composed of a vesicle that contains curved membrane protein complexes, that exert protrusive forces on the membrane (representing the pressure due to actin polymerization). This minimal-cell model gives rise to spontaneous emergence of a motile phenotype, driven by a lamellipodia-like leading edge. By systematically screening the behavior of this model on different types of curved substrates (sinusoidal, cylinder, and tube), we show that minimal ingredients and energy terms capture the experimental data. The model recovers the observed migration on the sinusoidal substrate, where cells move along the grooves (minima), while avoiding motion along the ridges. In addition, the model predicts the tendency of cells to migrate circumferentially on convex substrates and axially on concave ones. Both of these predictions are verified experimentally, on several cell types. Altogether, our results identify the minimization of membrane-substrate adhesion energy and binding energy between the membrane protein complexes as key players of curvotaxis in cell migration.

Mitotic outcomes and errors in fibrous environments.

During mitosis, cells round up and utilize the interphase adhesion sites within the fibrous extracellular matrix (ECM) as guidance cues to orient the mitotic spindles. Here, using suspended ECM-mimicking nanofiber networks, we explore mitotic outcomes and error distribution for various interphase cell shapes. Elongated cells attached to single fibers through two focal adhesion clusters (FACs) at their extremities result in perfect spherical mitotic cell bodies that undergo significant 3-dimensional (3D) displacement while being held by retraction fibers (RFs). Increasing the number of parallel fibers increases FACs and retraction fiber-driven stability, leading to reduced 3D cell body movement, metaphase plate rotations, increased interkinetochore distances, and significantly faster division times. Interestingly, interphase kite shapes on a crosshatch pattern of four fibers undergo mitosis resembling single-fiber outcomes due to rounded bodies being primarily held in position by RFs from two perpendicular suspended fibers. We develop a cortex-astral microtubule analytical model to capture the retraction fiber dependence of the metaphase plate rotations. We observe that reduced orientational stability, on single fibers, results in increased monopolar mitotic defects, while multipolar defects become dominant as the number of adhered fibers increases. We use a stochastic Monte Carlo simulation of centrosome, chromosome, and membrane interactions to explain the relationship between the observed propensity of monopolar and multipolar defects and the geometry of RFs. Overall, we establish that while bipolar mitosis is robust in fibrous environments, the nature of division errors in fibrous microenvironments is governed by interphase cell shapes and adhesion geometries.

A minimal cell model for lamellipodia-based cellular dynamics and migration.

One ubiquitous cellular structure for performing various tasks, such as spreading and migration over external surfaces, is the sheet-like protrusion called a lamellipodium, which propels the leading edge of the cell. Despite the detailed knowledge about the many components of this cellular structure, it is not yet fully understood how these components self-organize spatiotemporally to form lamellipodia. We review here recent theoretical works where we have demonstrated that membrane-bound protein complexes that have intrinsic curvature and recruit the protrusive forces of the cytoskeleton result in a simple, yet highly robust, organizing feedback mechanism that organizes the cytoskeleton and the membrane. This self-organization mechanism accounts for the formation of flat lamellipodia at the leading edge of cells spreading over adhesive substrates, allowing for the emergence of a polarized, motile 'minimal cell' model. The same mechanism describes how lamellipodia organize to drive robust engulfment of particles during phagocytosis and explains in simple physical terms the spreading and migration of cells over fibers and other curved surfaces. This Review highlights that despite the complexity of cellular composition, there might be simple general physical principles that are utilized by the cell to drive cellular shape dynamics.

Spatiotemporal dynamics of animal contests arise from effective forces between contestants.

Competition among animals for resources, notably food, territories, and mates, is ubiquitous at all scales of life. This competition is often resolved through contests among individuals, which are commonly understood according to their outcomes and in particular, how these outcomes depend on decision-making by the contestants. Because they are restricted to end-point predictions, these approaches cannot predict real-time or real-space dynamics of animal contest behavior. This limitation can be overcome by studying systems that feature typical contest behavior while being simple enough to track and model. Here, we propose to use such systems to construct a theoretical framework that describes real-time movements and behaviors of animal contestants. We study the spatiotemporal dynamics of contests in an orb-weaving spider, in which all the common elements of animal contests play out. The confined arena of the web, on which interactions are dominated by vibratory cues in a two-dimensional space, simplifies the analysis of interagent interactions. We ask whether these seemingly complex decision-makers can be modeled as interacting active particles responding only to effective forces of attraction and repulsion due to their interactions. By analyzing the emergent dynamics of "contestant particles," we provide mechanistic explanations for real-time dynamical aspects of animal contests, thereby explaining competitive advantages of larger competitors and demonstrating that complex decision-making need not be invoked in animal contests to achieve adaptive outcomes. Our results demonstrate that physics-based classification and modeling, in terms of effective rules of interaction, provide a powerful framework for understanding animal contest behaviors.

The geometry of decision-making in individuals and collectives.

Choosing among spatially distributed options is a central challenge for animals, from deciding among alternative potential food sources or refuges to choosing with whom to associate. Using an integrated theoretical and experimental approach (employing immersive virtual reality), we consider the interplay between movement and vectorial integration during decision-making regarding two, or more, options in space. In computational models of this process, we reveal the occurrence of spontaneous and abrupt "critical" transitions (associated with specific geometrical relationships) whereby organisms spontaneously switch from averaging vectorial information among, to suddenly excluding one among, the remaining options. This bifurcation process repeats until only one option-the one ultimately selected-remains. Thus, we predict that the brain repeatedly breaks multichoice decisions into a series of binary decisions in space-time. Experiments with fruit flies, desert locusts, and larval zebrafish reveal that they exhibit these same bifurcations, demonstrating that across taxa and ecological contexts, there exist fundamental geometric principles that are essential to explain how, and why, animals move the way they do.

Chemokine-biased robust self-organizing polarization of migrating cells in vivo.

To study the mechanisms controlling front-rear polarity in migrating cells, we used zebrafish primordial germ cells (PGCs) as an in vivo model. We find that polarity of bleb-driven migrating cells can be initiated at the cell front, as manifested by actin accumulation at the future leading edge and myosin-dependent retrograde actin flow toward the other side of the cell. In such cases, the definition of the cell front, from which bleb-inhibiting proteins such as Ezrin are depleted, precedes the establishment of the cell rear, where those proteins accumulate. Conversely, following cell division, the accumulation of Ezrin at the cleavage plane is the first sign for cell polarity and this aspect of the cell becomes the cell back. Together, the antagonistic interactions between the cell front and back lead to a robust polarization of the cell. Furthermore, we show that chemokine signaling can bias the establishment of the front-rear axis of the cell, thereby guiding the migrating cells toward sites of higher levels of the attractant. We compare these results to a theoretical model according to which a critical value of actin treadmilling flow can initiate a positive feedback loop that leads to the generation of the front-rear axis and to stable cell polarization. Together, our in vivo findings and the mathematical model, provide an explanation for the observed nonoriented migration of primordial germ cells in the absence of the guidance cue, as well as for the directed migration toward the region where the gonad develops.

Polarization and motility of one-dimensional multi-cellular trains.

Collective cell migration, whereby cells adhere to form multi-cellular clusters that move as a single entity, play an important role in numerous biological processes, such as during development and cancer progression. Recent experimental work focused on migration of one-dimensional cellular clusters, confined to move along adhesive lanes, as a simple geometry in which to systematically study this complex system. One-dimensional migration also arises in the body when cells migrate along blood vessels, axonal projections, and narrow cavities between tissues. We explore here the modes of one-dimensional migration of cellular clusters ("trains") by implementing cell-cell interactions in a model of cell migration that contains a mechanism for spontaneous cell polarization. We go beyond simple phenomenological models of the cells as self-propelled particles by having the internal polarization of each cell depend on its interactions with the neighboring cells that directly affect the actin polymerization activity at the cell's leading edges. Both contact inhibition of locomotion and cryptic lamellipodia interactions between neighboring cells are introduced. We find that this model predicts multiple motility modes of the cell trains, which can have several different speeds for the same polarization pattern. Compared to experimental data, we find that Madin-Darby canine kidney cells are poised along the transition region where contact inhibition of locomotion and cryptic lamellipodia roughly balance each other, where collective migration speed is most sensitive to the values of the cell-cell interaction strength.

A simple cognitive model explains movement decisions in zebrafish while following leaders.

While moving, animals must frequently make decisions about their future travel direction, whether they are alone or in a group. Here we investigate this process for zebrafish (Danio rerio), which naturally move in cohesive groups. Employing state-of-the-art virtual reality, we study how real fish (RF) follow one or several moving, virtual conspecifics (leaders). These data are used to inform, and test, a model of social response that includes a process of explicit decision-making, whereby the fish can decide which of the virtual conspecifics to follow, or to follow in some average direction. This approach is in contrast with previous models where the direction of motion was based on a continuous computation, such as directional averaging. Building upon a simplified version of this model (Sridharet al2021Proc. Natl Acad. Sci.118e2102157118), which was limited to a one-dimensional projection of the fish motion, we present here a model that describes the motion of the RF as it swims freely in two-dimensions. Motivated by experimental observations, the swim speed of the fish in this model uses a burst-and-coast swimming pattern, with the burst frequency being dependent on the distance of the fish from the followed conspecific(s). We demonstrate that this model is able to explain the observed spatial distribution of the RF behind the virtual conspecifics in the experiments, as a function of their average speed and number. In particular, the model naturally explains the observed critical bifurcations for a freely swimming fish, which appear in the spatial distributions whenever the fish makes a decision to follow only one of the virtual conspecifics, instead of following them as an averaged group. This model can provide the foundation for modeling a cohesive shoal of swimming fish, while explicitly describing their directional decision-making process at the individual level.

Cell confinement reveals a branched-actin independent circuit for neutrophil polarity.

Migratory cells use distinct motility modes to navigate different microenvironments, but it is unclear whether these modes rely on the same core set of polarity components. To investigate this, we disrupted actin-related protein 2/3 (Arp2/3) and the WASP-family verprolin homologous protein (WAVE) complex, which assemble branched actin networks that are essential for neutrophil polarity and motility in standard adherent conditions. Surprisingly, confinement rescues polarity and movement of neutrophils lacking these components, revealing a processive bleb-based protrusion program that is mechanistically distinct from the branched actin-based protrusion program but shares some of the same core components and underlying molecular logic. We further find that the restriction of protrusion growth to one site does not always respond to membrane tension directly, as previously thought, but may rely on closely linked properties such as local membrane curvature. Our work reveals a hidden circuit for neutrophil polarity and indicates that cells have distinct molecular mechanisms for polarization that dominate in different microenvironments.

A theoretical model of efficient phagocytosis driven by curved membrane proteins and active cytoskeleton forces.

Phagocytosis is the process of engulfment and internalization of comparatively large particles by cells, and plays a central role in the functioning of our immune system. We study the process of phagocytosis by considering a simplified coarse grained model of a three-dimensional vesicle, having a uniform adhesion interaction with a rigid particle, and containing curved membrane-bound protein complexes or curved membrane nano-domains, which in turn recruit active cytoskeletal forces. Complete engulfment is achieved when the bending energy cost of the vesicle is balanced by the gain in the adhesion energy. The presence of curved (convex) proteins reduces the bending energy cost by self-organizing with a higher density at the highly curved leading edge of the engulfing membrane, which forms the circular rim of the phagocytic cup that wraps around the particle. This allows the engulfment to occur at much smaller adhesion strength. When the curved membrane-bound protein complexes locally recruit actin polymerization machinery, which leads to outward forces being exerted on the membrane, we found that engulfment is achieved more quickly and at a lower protein density. We consider spherical and non-spherical particles and found that non-spherical particles are more difficult to engulf in comparison to the spherical particles of the same surface area. For non-spherical particles, the engulfment time crucially depends on the initial orientation of the particles with respect to the vesicle. Our model offers a mechanism for the spontaneous self-organization of the actin cytoskeleton at the phagocytic cup, in good agreement with recent high-resolution experimental observations.

Local actin dynamics couple speed and persistence in a cellular Potts model of cell migration.

Cell migration is astoundingly diverse. Molecular signatures, cell-cell interactions, and environmental structures each play their part in shaping cell motion, yielding numerous morphologies and migration modes. Nevertheless, in recent years, a simple unifying law was found to describe cell migration across many different cell types and contexts: faster cells turn less frequently. This universal coupling between speed and persistence (UCSP) was explained by retrograde actin flow from front to back, but it remains unclear how this mechanism generalizes to cells with complex shapes and cells migrating in structured environments, which may not have a well-defined front-to-back orientation. Here, we present an in-depth characterization of an existing cellular Potts model, in which cells polarize dynamically from a combination of local actin dynamics (stimulating protrusions) and global membrane tension along the perimeter (inhibiting protrusions). We first show that the UCSP emerges spontaneously in this model through a cross talk of intracellular mechanisms, cell shape, and environmental constraints, resembling the dynamic nature of cell migration in vivo. Importantly, we find that local protrusion dynamics suffice to reproduce the UCSP-even in cases in which no clear global, front-to-back polarity exists. We then harness the spatial nature of the cellular Potts model to show how cell shape dynamics limit both the speed and persistence a cell can reach and how a rigid environment such as the skin can restrict cell motility even further. Our results broaden the range of potential mechanisms underlying the speed-persistence coupling that has emerged as a fundamental property of migrating cells.

Cell cluster migration: Connecting experiments with physical models.

In multi-cellular organisms, the migration of cohesive clusters of cells containing many individual cells is a common occurrence. Examples include the migration of cells during processes such as the development of the embryo, wound healing, immune response, and the spread of cancer. The migration process depends not only on the traction forces applied by the cluster on its surroundings, in order to move, but also on the viscoelastic properties of both the surrounding matrix and the migrating cellular cluster. Characterizing the viscoelastic properties of the cluster, its environment and the forces within the cluster, in great detail, is difficult both in-vitro and certainly in-vivo. We review here several examples where theoretical studies using simplified models can be used to gain insights into the basic underlying mechanisms that control the cellular migration patterns.

A random first-order transition theory for an active glass.

How does nonequilibrium activity modify the approach to a glass? This is an important question, since many experiments reveal the near-glassy nature of the cell interior, remodeled by activity. However, different simulations of dense assemblies of active particles, parametrized by a self-propulsion force, [Formula: see text], and persistence time, [Formula: see text], appear to make contradictory predictions about the influence of activity on characteristic features of glass, such as fragility. This calls for a broad conceptual framework to understand active glasses; here, we extend the random first-order transition (RFOT) theory to a dense assembly of self-propelled particles. We compute the active contribution to the configurational entropy through an effective model of a single particle in a caging potential. This simple active extension of RFOT provides excellent quantitative fits to existing simulation results. We find that whereas [Formula: see text] always inhibits glassiness, the effect of [Formula: see text] is more subtle and depends on the microscopic details of activity. In doing so, the theory automatically resolves the apparent contradiction between the simulation models. The theory also makes several testable predictions, which we verify by both existing and new simulation data, and should be viewed as a step toward a more rigorous analytical treatment of active glass.

Active Trap Model.

Motivated by the dynamics of particles embedded in active gels, both in vitro and inside the cytoskeleton of living cells, we study an active generalization of the classical trap model. We demonstrate that activity leads to dramatic modifications in the diffusion compared to the thermal case: the mean square displacement becomes subdiffusive, spreading as a power law in time, when the trap depth distribution is a Gaussian and is slower than any power law when it is drawn from an exponential distribution. The results are derived for a simple, exactly solvable, case of harmonic traps. We then argue that the results are robust for more realistic trap shapes when the activity is strong.

Cellular Blebs and Membrane Invaginations Are Coupled through Membrane Tension Buffering.

Bleb-type cellular protrusions play key roles in a range of biological processes. It was recently found that bleb growth is facilitated by a local supply of membrane from tubular invaginations, but the interplay between the expanding bleb and the membrane tubes remains poorly understood. On the one hand, the membrane area stored in tubes may serve as a reservoir for bleb expansion. On the other hand, the sequestering of excess membrane in stabilized invaginations may effectively increase the cell membrane tension, which suppresses spontaneous protrusions. Here, we investigate this duality through physical modeling and in vivo experiments. In agreement with observations, our model describes the transition into a tube-flattening mode of bleb expansion while also predicting that the blebbing rate is impaired by elevating the concentration of the curved membrane proteins that form the tubes. We show both theoretically and experimentally that the stabilizing effect of tubes could be counterbalanced by the cortical myosin contractility. Our results largely suggest that proteins able to induce membrane tubulation, such as those containing N-BAR domains, can buffer the effective membrane tension-a master regulator of all cell deformations.

Bi-stability in cooperative transport by ants in the presence of obstacles.

To cooperatively carry large food items to the nest, individual ants conform their efforts and coordinate their motion. Throughout this expedition, collective motion is driven both by internal interactions between the carrying ants and a response to newly arrived informed ants that orient the cargo towards the nest. During the transport process, the carrying group must overcome obstacles that block their path to the nest. Here, we investigate the dynamics of cooperative transport, when the motion of the ants is frustrated by a linear obstacle that obstructs the motion of the cargo. The obstacle contains a narrow opening that serves as the only available passage to the nest, and through which single ants can pass but not with the cargo. We provide an analytical model for the ant-cargo system in the constrained environment that predicts a bi-stable dynamic behavior between an oscillatory mode of motion along the obstacle and a convergent mode of motion near the opening. Using both experiments and simulations, we show how for small cargo sizes, the system exhibits spontaneous transitions between these two modes of motion due to fluctuations in the applied force on the cargo. The bi-stability provides two possible problem solving strategies for overcoming the obstacle, either by attempting to pass through the opening, or take large excursions to circumvent the obstacle.

Theoretical study of vesicle shapes driven by coupling curved proteins and active cytoskeletal forces.

Eukaryote cells have a flexible shape, which dynamically changes according to the function performed by the cell. One mechanism for deforming the cell membrane into the desired shape is through the expression of curved membrane proteins. Furthermore, these curved membrane proteins are often associated with the recruitment of the cytoskeleton, which then applies active forces that deform the membrane. This coupling between curvature and activity was previously explored theoretically in the linear limit of small deformations, and low dimensionality. Here we explore the unrestricted shapes of vesicles that contain active curved membrane proteins, in three-dimensions, using Monte-Carlo numerical simulations. The activity of the proteins is in the form of protrusive forces that push the membrane outwards, as may arise from the cytoskeleton of the cell due to actin or microtubule polymerization occurring near the membrane. For proteins that have an isotropic convex shape, the additional protrusive force enhances their tendency to aggregate and form membrane protrusions (buds). In addition, we find another transition from deformed spheres with necklace type aggregates, to flat pancake-shaped vesicles, where the curved proteins line the outer rim. This second transition is driven by the active forces, coupled to the spontaneous curvature, and the resulting configurations may shed light on the formation of sheet-like protrusions and lamellipodia of adhered and motile cells.

Emergent oscillations assist obstacle negotiation during ant cooperative transport.

Collective motion by animal groups is affected by internal interactions, external constraints, and the influx of information. A quantitative understanding of how these different factors give rise to different modes of collective motion is, at present, lacking. Here, we study how ants that cooperatively transport a large food item react to an obstacle blocking their path. Combining experiments with a statistical physics model of mechanically coupled active agents, we show that the constraint induces a deterministic collective oscillatory mode that facilitates obstacle circumvention. We provide direct experimental evidence, backed by theory, that this motion is an emergent group effect that does not require any behavioral changes at the individual level. We trace these relaxation oscillations to the interplay between two forces; informed ants pull the load toward the nest whereas uninformed ants contribute to the motion's persistence along the tangential direction. The model's predictions that oscillations appear above a critical system size, that the group can spontaneously transition into its ordered phase, and that the system can exhibit complete rotations are all verified experimentally. We expect that similar oscillatory modes emerge in collective motion scenarios where the structure of the environment imposes conflicts between individually held information and the group's tendency for cohesiveness.

Modeling and analysis of collective cell migration in an in vivo three-dimensional environment.

A long-standing question in collective cell migration has been what might be the relative advantage of forming a cluster over migrating individually. Does an increase in the size of a collectively migrating group of cells enable them to sample the chemical gradient over a greater distance because the difference between front and rear of a cluster would be greater than for single cells? We combined theoretical modeling with experiments to study collective migration of the border cells in-between nurse cells in the Drosophila egg chamber. We discovered that cluster size is positively correlated with migration speed, up to a particular point above which speed plummets. This may be due to the effect of viscous drag from surrounding nurse cells together with confinement of all of the cells within a stiff extracellular matrix. The model predicts no relationship between cluster size and velocity for cells moving on a flat surface, in contrast to movement within a 3D environment. Our analyses also suggest that the overall chemoattractant profile in the egg chamber is likely to be exponential, with the highest concentration in the oocyte. These findings provide insights into collective chemotaxis by combining theoretical modeling with experimentation.

Theory of Epithelial Cell Shape Transitions Induced by Mechanoactive Chemical Gradients.

Cell shape is determined by a balance of intrinsic properties of the cell as well as its mechanochemical environment. Inhomogeneous shape changes underlie many morphogenetic events and involve spatial gradients in active cellular forces induced by complex chemical signaling. Here, we introduce a mechanochemical model based on the notion that cell shape changes may be induced by external diffusible biomolecules that influence cellular contractility (or equivalently, adhesions) in a concentration-dependent manner-and whose spatial profile in turn is affected by cell shape. We map out theoretically the possible interplay between chemical concentration and cellular structure. Besides providing a direct route to spatial gradients in cell shape profiles in tissues, we show that the dependence on cell shape helps create robust mechanochemical gradients.