An adhesion-dependent switch between mechanisms that determine motile cell shape

Increased vesicular dynamics and nanoscale clustering of IL-2 after T cell activation

T cells coordinate intercellular communication through the meticulous regulation of cytokine secretion. Direct visualization of vesicular transport and intracellular distribution of cytokines provides valuable insights into the temporal and spatial mechanisms involved in regulation. Employing Jurkat E6-1 T cells and interleukin-2 (IL-2) as a model system, we investigated vesicular dynamics using single-particle tracking and the nanoscale distribution of intracellular IL-2 in fixed T cells using superresolution microscopy. Live-cell imaging revealed that in vitro activation resulted in increased vesicular dynamics. Direct stochastic optical reconstruction microscopy and 3D structured illumination microscopy revealed nanoscale clustering of IL-2. In vitro activation correlated with spatial accumulation of IL-2 nanoclusters into more pronounced and elongated clusters. These observations provide visual evidence that accelerated vesicular transport and spatial concatenation of IL-2 clusters at the nanoscale may constitute a potential mechanism for modulating cytokine release by Jurkat T cells.

Mesenchymal cell migration on one-dimensional micropatterns

Quantitative studies of mesenchymal cell motion are important to elucidate cytoskeleton function and mechanisms of cell migration. To this end, confinement of cell motion to one dimension (1D) significantly simplifies the problem of cell shape in experimental and theoretical investigations. Here we review 1D migration assays employing micro-fabricated lanes and reflect on the advantages of such platforms. Data are analyzed using biophysical models of cell migration that reproduce the rich scenario of morphodynamic behavior found in 1D. We describe basic model assumptions and model behavior. It appears that mechanical models explain the occurrence of universal relations conserved across different cell lines such as the adhesion-velocity relation and the universal correlation between speed and persistence (UCSP). We highlight the unique opportunity of reproducible and standardized 1D assays to validate theory based on statistical measures from large data of trajectories and discuss the potential of experimental settings embedding controlled perturbations to probe response in migratory behavior.

Influence of Curvature on Cell Motility and Morphology during Cancer Migration in Confined Microchannels

Microchannels often serve as highways for cancer migration, and their topology largely determines the migration efficiency. Curvature, a topological parameter in biological systems, has recently been reported to be efficient in guiding cell polarization and migration. Curvature varies widely along curved microchannels, while its influence on cell migration remains elusive. Here, we recapitulated the curved microchannels, as observed in clinical tumor tissues with hydrogels, and studied how cancer cells respond to curvature. We found that cells bend more significantly in a larger curvature and exhibit less spreading as well as lower motility. The underlying mechanism is probably based on the hindrance of the movement of cytoskeletal molecules at the curved microchannel walls. Collectively, our results demonstrated that the accelerated actin retrograde flow rate under local curvature has an effective negative regulation on cell motility and morphology, leading to shortened and bent cell morphologies as well as hampered cell migration efficiency.

Three-component contour dynamics model to simulate and analyze amoeboid cell motility in two dimensions

Amoeboid cell motility is relevant in a wide variety of biomedical processes such as wound healing, cancer metastasis, and embryonic morphogenesis. It is characterized by pronounced changes of the cell shape associated with expansions and retractions of the cell membrane, which result in a crawling kind of locomotion. Despite existing computational models of amoeboid motion, the inference of expansion and retraction components of individual cells, the corresponding classification of cells, and the a priori specification of the parameter regime to achieve a specific motility behavior remain challenging open problems. We propose a novel model of the spatio-temporal evolution of two-dimensional cell contours comprising three biophysiologically motivated components: a stochastic term accounting for membrane protrusions and two deterministic terms accounting for membrane retractions by regularizing the shape and area of the contour. Mathematically, these correspond to the intensity of a self-exciting Poisson point process, the area-preserving curve-shortening flow, and an area adjustment flow. The model is used to generate contour data for a variety of qualitatively different, e.g., polarized and non-polarized, cell tracks that visually resemble experimental data very closely. In application to experimental cell tracks, we inferred the protrusion component and examined its correlation to common biomarkers: the F-actin density close to the membrane and its local motion. Due to the low model complexity, parameter estimation is fast, straightforward, and offers a simple way to classify contour dynamics based on two locomotion types: the amoeboid and a so-called fan-shaped type. For both types, we use cell tracks segmented from fluorescence imaging data of the model organism Dictyostelium discoideum. An implementation of the model is provided within the open-source software package AmoePy, a Python-based toolbox for analyzing and simulating amoeboid cell motility.

Myosin-independent stiffness sensing by fibroblasts is regulated by the viscoelasticity of flowing actin

The stiffness of the extracellular matrix induces differential tension within integrin-based adhesions, triggering differential mechanoresponses. However, it has been unclear if the stiffness-dependent differential tension is induced solely by myosin activity. Here, we report that in the absence of myosin contractility, 3T3 fibroblasts still transmit stiffness-dependent differential levels of traction. This myosin-independent differential traction is regulated by polymerizing actin assisted by actin nucleators Arp2/3 and formin where formin has a stronger contribution than Arp2/3 to both traction and actin flow. Intriguingly, despite only slight changes in F-actin flow speed observed in cells with the combined inhibition of Arp2/3 and myosin compared to cells with sole myosin inhibition, they show a 4-times reduction in traction than cells with myosin-only inhibition. Our analyses indicate that traditional models based on rigid F-actin are inadequate for capturing such dramatic force reduction with similar actin flow. Instead, incorporating the F-actin network's viscoelastic properties is crucial. Our new model including the F-actin viscoelasticity reveals that Arp2/3 and formin enhance stiffness sensitivity by mechanically reinforcing the F-actin network, thereby facilitating more effective transmission of flow-induced forces. This model is validated by cell stiffness measurement with atomic force microscopy and experimental observation of model-predicted stiffness-dependent actin flow fluctuation.

The excitable nature of polymerizing actin and the Belousov-Zhabotinsky reaction

The intricate regulatory processes behind actin polymerization play a crucial role in cellular biology, including essential mechanisms such as cell migration or cell division. However, the self-organizing principles governing actin polymerization are still poorly understood. In this perspective article, we compare the Belousov-Zhabotinsky (BZ) reaction, a classic and well understood chemical oscillator known for its self-organizing spatiotemporal dynamics, with the excitable dynamics of polymerizing actin. While the BZ reaction originates from the domain of inorganic chemistry, it shares remarkable similarities with actin polymerization, including the characteristic propagating waves, which are influenced by geometry and external fields, and the emergent collective behavior. Starting with a general description of emerging patterns, we elaborate on single droplets or cell-level dynamics, the influence of geometric confinements and conclude with collective interactions. Comparing these two systems sheds light on the universal nature of self-organization principles in both living and inanimate systems.

Motility of Synthetic Cells from Engineered Lipids

Synthetic cells are artificial systems that resemble natural cells. Significant efforts have been made over the years to construct synthetic protocells that can mimic biological mechanisms and perform various complex processes. These include compartmentalization, metabolism, energy supply, communication, and gene reproduction. Cell motility is also of great importance, as nature uses elegant mechanisms for intracellular trafficking, immune response, and embryogenesis. In this review, we discuss the motility of synthetic cells made from lipid vesicles and relevant molecular mechanisms. Synthetic cell motion may be classified into surface-based or solution-based depending on whether it involves interactions with surfaces or movement in fluids. Collective migration behaviors have also been demonstrated. The swarm motion requires additional mechanisms for intercellular signaling and directional motility that enable communication and coordination among the synthetic vesicles. In addition, intracellular trafficking for molecular transport has been reconstituted in minimal cells with the help of DNA nanotechnology. These efforts demonstrate synthetic cells that can move, detect, respond, and interact. We envision that new developments in protocell motility will enhance our understanding of biological processes and be instrumental in bioengineering and therapeutic applications.

Friction patterns guide actin network contraction

The shape of cells is the outcome of the balance of inner forces produced by the actomyosin network and the resistive forces produced by cell adhesion to their environment. The specific contributions of contractile, anchoring and friction forces to network deformation rate and orientation are difficult to disentangle in living cells where they influence each other. Here, we reconstituted contractile actomyosin networks in vitro to study specifically the role of the friction forces between the network and its anchoring substrate. To modulate the magnitude and spatial distribution of friction forces, we used glass or lipids surface micropatterning to control the initial shape of the network. We adapted the concentration of Nucleating Promoting Factor on each surface to induce the assembly of actin networks of similar densities and compare the deformation of the network toward the centroid of the pattern shape upon myosin-induced contraction. We found that actin network deformation was faster and more coordinated on lipid bilayers than on glass, showing the resistance of friction to network contraction. To further study the role of the spatial distribution of these friction forces, we designed heterogeneous micropatterns made of glass and lipids. The deformation upon contraction was no longer symmetric but biased toward the region of higher friction. Furthermore, we showed that the pattern of friction could robustly drive network contraction and dominate the contribution of asymmetric distributions of myosins. Therefore, we demonstrate that during contraction, both the active and resistive forces are essential to direct the actin network deformation.

G-actin diffusion is insufficient to achieve F-actin assembly in fast-treadmilling protrusions

To generate forces that drive migration of a eukaryotic cell, arrays of actin filaments (F-actin) are assembled at the cell's leading membrane edge. To maintain cell propulsion and respond to dynamic external cues, actin filaments must be disassembled to regenerate the actin monomers (G-actin), and transport of G-actin from sites of disassembly back to the leading edge completes the treadmilling cycle and limits the flux of F-actin assembly. Whether or not molecular diffusion is sufficient for G-actin transport has been a long-standing topic of debate, in part because the dynamic nature of cell motility and migration hinders the estimation of transport parameters. In this work, we applied an experimental system in which cells adopt an approximately constant and symmetrical shape; they cannot migrate but exhibit fast, steady treadmilling in the thin region protruding from the cell. Using fluorescence recovery after photobleaching, we quantified the relative concentrations and corresponding fluxes of F- and G-actin in this system. In conjunction with mathematical modeling, constrained by measured features of each region of interest, this approach revealed that diffusion alone cannot account for the transport of G-actin to the leading edge. Although G-actin diffusion and vectorial transport might vary with position in the protruding region, good agreement with the fluorescence recovery after photobleaching measurements was achieved by a model with constant G-actin diffusivity ∼2 μm2/s and anterograde G-actin velocity less than 1 μm/s.

Crawling, waving, inch worming, dilating, and pivoting mechanics of migrating cells: Lessons from Ken Jacobson

Research on the locomotion of single cells on hard, flat surfaces brought insight into the mechanisms of leading-edge protrusion, spatially graded adhesion, front-rear coordination, and how intracellular and traction forces are harnessed to execute various maneuvers. Here, we highlight how, by studying a variety of cell types, shapes, and movements, Ken Jacobson and his collaborators made several discoveries that triggered the mechanistic understanding of cell motility. We then review the recent advancements and current perspectives in this field.

Whole-genome screens reveal regulators of differentiation state and context-dependent migration in human neutrophils

Neutrophils are the most abundant leukocyte in humans and provide a critical early line of defense as part of our innate immune system. We perform a comprehensive, genome-wide assessment of the molecular factors critical to proliferation, differentiation, and cell migration in a neutrophil-like cell line. Through the development of multiple migration screen strategies, we specifically probe directed (chemotaxis), undirected (chemokinesis), and 3D amoeboid cell migration in these fast-moving cells. We identify a role for mTORC1 signaling in cell differentiation, which influences neutrophil abundance, survival, and migratory behavior. Across our individual migration screens, we identify genes involved in adhesion-dependent and adhesion-independent cell migration, protein trafficking, and regulation of the actomyosin cytoskeleton. This genome-wide screening strategy, therefore, provides an invaluable approach to the study of neutrophils and provides a resource that will inform future studies of cell migration in these and other rapidly migrating cells.

Synchronization in collectively moving inanimate and living active matter

Whether one considers swarming insects, flocking birds, or bacterial colonies, collective motion arises from the coordination of individuals and entails the adjustment of their respective velocities. In particular, in close confinements, such as those encountered by dense cell populations during development or regeneration, collective migration can only arise coordinately. Yet, how individuals unify their velocities is often not understood. Focusing on a finite number of cells in circular confinements, we identify waves of polymerizing actin that function as a pacemaker governing the speed of individual cells. We show that the onset of collective motion coincides with the synchronization of the wave nucleation frequencies across the population. Employing a simpler and more readily accessible mechanical model system of active spheres, we identify the synchronization of the individuals' internal oscillators as one of the essential requirements to reach the corresponding collective state. The mechanical 'toy' experiment illustrates that the global synchronous state is achieved by nearest neighbor coupling. We suggest by analogy that local coupling and the synchronization of actin waves are essential for the emergent, self-organized motion of cell collectives.

Adhesion tunes speed and persistence by coordinating protrusions and extracellular matrix remodeling

Cell migration through 3D environments is essential to development, disease, and regeneration processes. Conceptual models of migration have been developed primarily on the basis of 2D cell behaviors, but a general understanding of 3D cell migration is still lacking due to the added complexity of the extracellular matrix. Here, using a multiplexed biophysical imaging approach for single-cell analysis of human cell lines, we show how the subprocesses of adhesion, contractility, actin cytoskeletal dynamics, and matrix remodeling integrate to produce heterogeneous migration behaviors. This single-cell analysis identifies three modes of cell speed and persistence coupling, driven by distinct modes of coordination between matrix remodeling and protrusive activity. The framework that emerges establishes a predictive model linking cell trajectories to distinct subprocess coordination states.

From actin waves to mechanism and back: How theory aids biological understanding

Actin dynamics in cell motility, division, and phagocytosis is regulated by complex factors with multiple feedback loops, often leading to emergent dynamic patterns in the form of propagating waves of actin polymerization activity that are poorly understood. Many in the actin wave community have attempted to discern the underlying mechanisms using experiments and/or mathematical models and theory. Here, we survey methods and hypotheses for actin waves based on signaling networks, mechano-chemical effects, and transport characteristics, with examples drawn from Dictyostelium discoideum, human neutrophils, Caenorhabditis elegans, and Xenopus laevis oocytes. While experimentalists focus on the details of molecular components, theorists pose a central question of universality: Are there generic, model-independent, underlying principles, or just boundless cell-specific details? We argue that mathematical methods are equally important for understanding the emergence, evolution, and persistence of actin waves and conclude with a few challenges for future studies.

Morphomigrational description as a new approach connecting cell's migration with its morphology

The examination of morphology and migration of cells plays substantial role in understanding the cellular behaviour, being described by plethora of quantitative parameters and models. These descriptions, however, treat cell migration and morphology as independent properties of temporal cell state, while not taking into account their strong interdependence in adherent cells. Here we present the new and simple mathematical parameter called signed morphomigrational angle (sMM angle) that links cell geometry with translocation of cell centroid, considering them as one morphomigrational behaviour. The sMM angle combined with pre-existing quantitative parameters enabled us to build a new tool called morphomigrational description, used to assign the numerical values to several cellular behaviours. Thus, the cellular activities that until now were characterized using verbal description or by complex mathematical models, are described here by a set of numbers. Our tool can be further used in automatic analysis of cell populations as well as in studies focused on cellular response to environmental directional signals.

Stiffness sensing by smooth muscle cells: Continuum mechanics modeling of the acto-myosin role

Aortic smooth muscle cells (SMCs) play a vital role in maintaining homeostasis in the aorta by sensing and responding to mechanical stimuli. However, the mechanisms that underlie the ability of SMCs to sense and respond to stiffness change in their environment are still partially unclear. In this study, we focus on the role of acto-myosin contractility in stiffness sensing and introduce a novel continuum mechanics approach based on the principles of thermal strains. Each stress fiber satisfies a universal stress-strain relationship driven by a Young's modulus, a contraction coefficient scaling the fictitious thermal strain, a maximum contraction stress and a softening parameter describing the sliding effects between actin and myosin filaments. To account for the inherent variability of cellular responses, large populations of SMCs are modeled with the finite-element method, each cell having a random number and a random arrangement of stress fibers. Moreover, the level of myosin activation in each stress fiber satisfies a Weibull probability density function. Model predictions are compared to traction force measurements on different SMC lineages. It is demonstrated that the model not only predicts well the effects of substrate stiffness on cellular traction, but it can also successfully approximate the statistical variations of cellular tractions induced by intercellular variability. Finally, stresses in the nuclear envelope and in the nucleus are computed with the model, showing that the variations of cytoskeletal forces induced by substrate stiffness directly induce deformations of the nucleus which can potentially alter gene expression. The predictability of the model combined to its relative simplicity are promising assets for further investigation of stiffness sensing in 3D environments. Eventually, this could contribute to decipher the effects of mechanosensitivity impairment, which are known to be at the root of aortic aneurysms.

Positive, negative and controlled durotaxis

Cell migration is a physical process central to life. Among others, it regulates embryogenesis, tissue regeneration, and tumor growth. Therefore, understanding and controlling cell migration represent fundamental challenges in science. Specifically, the ability of cells to follow stiffness gradients, known as durotaxis, is ubiquitous across most cell types. Even so, certain cells follow positive stiffness gradients while others move along negative gradients. How the physical mechanisms involved in cell migration work to enable a wide range of durotactic responses is still poorly understood. Here, we provide a mechanistic rationale for durotaxis by integrating stochastic clutch models for cell adhesion with an active gel theory of cell migration. We show that positive and negative durotaxis found across cell types are explained by asymmetries in the cell adhesion dynamics. We rationalize durotaxis by asymmetric mechanotransduction in the cell adhesion behavior that further polarizes the intracellular retrograde flow and the protruding velocity at the cell membrane. Our theoretical framework confirms previous experimental observations and explains positive and negative durotaxis. Moreover, we show how durotaxis can be engineered to manipulate cell migration, which has important implications in biology, medicine, and bioengineering.

On multistability and constitutive relations of cell motion on fibronectin lanes

Cell motility on flat substrates exhibits coexisting steady and oscillatory morphodynamics, the biphasic adhesion-velocity relation, and the universal correlation between speed and persistence (UCSP) as simultaneous observations common to many cell types. Their universality and concurrency suggest a unifying mechanism causing all three of them. Stick-slip models for cells on one-dimensional lanes suggest multistability to arise from the nonlinear friction of retrograde flow. This study suggests a mechanical mechanism controlled by integrin signaling on the basis of a biophysical model and analysis of trajectories of MDA-MB-231 cells on fibronectin lanes, which additionally explains the constitutive relations. The experiments exhibit cells with steady or oscillatory morphodynamics and either spread or moving with spontaneous transitions between the dynamic regimes, spread and moving, and spontaneous direction reversals. Our biophysical model is based on the force balance at the protrusion edge, the noisy clutch of retrograde flow, and a response function of friction and membrane drag to integrin signaling. The theory reproduces the experimentally observed cell states, characteristics of oscillations, and state probabilities. Analysis of experiments with the biophysical model establishes a stick-slip oscillation mechanism, and explains multistability of cell states and the statistics of state transitions. It suggests protrusion competition to cause direction reversal events, the statistics of which explain the UCSP. The effect of integrin signaling on drag and friction explains the adhesion-velocity relation and cell behavior at fibronectin density steps. The dynamics of our mechanism are nonlinear flow mechanics driven by F-actin polymerization and shaped by the noisy clutch of retrograde flow friction, protrusion competition via membrane tension, and drag forces. Integrin signaling controls the parameters of the mechanical system.

Coupling cell shape and velocity leads to oscillation and circling in keratocyte galvanotaxis

During wound healing, fish keratocyte cells undergo galvanotaxis where they follow a wound-induced electric field. In addition to their stereotypical persistent motion, keratocytes can develop circular motion without a field or oscillate while crawling in the field direction. We developed a coarse-grained phenomenological model that captures these keratocyte behaviors. We fit this model to experimental data on keratocyte response to an electric field being turned on. A critical element of our model is a tendency for cells to turn toward their long axis, arising from a coupling between cell shape and velocity, which gives rise to oscillatory and circular motion. Galvanotaxis is influenced not only by the field-dependent responses, but also cell speed and cell shape relaxation rate. When the cell reacts to an electric field being turned on, our model predicts that stiff, slow cells react slowly but follow the signal reliably. Cells that polarize and align to the field at a faster rate react more quickly and follow the signal more reliably. When cells are exposed to a field that switches direction rapidly, cells follow the average of field directions, while if the field is switched more slowly, cells follow a "staircase" pattern. Our study indicated that a simple phenomenological model coupling cell speed and shape is sufficient to reproduce a broad variety of different keratocyte behaviors, ranging from circling to oscillation to galvanotactic response, by only varying a few parameters.

Confined keratocytes mimic in vivo migration and reveal volume-speed relationship

Fish basal epidermal cells, known as keratocytes, are well-suited for cell migration studies. In vitro, isolated keratocytes adopt a stereotyped shape with a large fan-shaped lamellipodium and a nearly spherical cell body. However, in their native in vivo environment, these cells adopt a significantly different shape during their rapid migration toward wounds. Within the epidermis, keratocytes experience two-dimensional (2D) confinement between the outer epidermal cell layer and the basement membrane; these two deformable surfaces constrain keratocyte cell bodies to be flatter in vivo than in isolation. In vivo keratocytes also exhibit a relative elongation of the front-to-back axis and substantially more lamellipodial ruffling, as compared to isolated cells. We have explored the effects of 2D confinement, separated from other in vivo environmental cues, by overlaying isolated cells with an agarose hydrogel with occasional spacers, or with a ceiling made of polydimethylsiloxane (PDMS) elastomer. Under these conditions, isolated keratocytes more closely resemble the in vivo migratory shape phenotype, displaying a flatter apical-basal axis and a longer front-to-back axis than unconfined keratocytes. We propose that 2D confinement contributes to multiple dimensions of in vivo keratocyte shape determination. Further analysis demonstrates that confinement causes a synchronous 20% decrease in both cell speed and volume. Interestingly, we were able to replicate the 20% decrease in speed using a sorbitol hypertonic shock to shrink the cell volume, which did not affect other aspects of cell shape. Collectively, our results suggest that environmentally imposed changes in cell volume may influence cell migration speed, potentially by perturbing physical properties of the cytoplasm.

On the significance of membrane unfolding in mechanosensitive cell spreading: Its individual and synergistic effects

Mechanosensitivity of cell spread area to substrate stiffness has been established both through experiments and different types of mathematical models of varying complexity including both the mechanics and biochemical reactions in the cell. What has not been addressed in previous mathematical models is the role of cell membrane dynamics on cell spreading, and an investigation of this issue is the goal of this work. We start with a simple mechanical model of cell spreading on a deformable substrate and progressively layer mechanisms to account for the traction dependent growth of focal adhesions, focal adhesion induced actin polymerization, membrane unfolding/exocytosis and contractility. This layering approach is intended to progressively help in understanding the role each mechanism plays in reproducing experimentally observed cell spread areas. To model membrane unfolding we introduce a novel approach based on defining an active rate of membrane deformation that is dependent on membrane tension. Our modeling approach allows us to show that tension-dependent membrane unfolding plays a critical role in achieving the large cell spread areas experimentally observed on stiff substrates. We also demonstrate that coupling between membrane unfolding and focal adhesion induced polymerization works synergistically to further enhance cell spread area sensitivity to substrate stiffness. This enhancement has to do with the fact that the peripheral velocity of spreading cells is associated with contributions from the different mechanisms by either enhancing the polymerization velocity at the leading edge or slowing down of the retrograde flow of actin within the cell. The temporal evolution of this balance in the model corresponds to the three-phase behavior observed experimentally during spreading. In the initial phase membrane unfolding is found to be particularly important.

The nonlinear motion of cells subject to external forces

To develop a minimal model for a cell moving in a crowded environment such as in tissue, we investigate the response of a liquid drop of active matter moving on a flat rigid substrate to forces applied at its boundaries. We consider two different self-propulsion mechanisms, active stresses and treadmilling polymerisation, and we investigate how the active drop motion is altered by these surface forces. We find a highly non-linear response to forces that we characterise using drop velocity, drop shape, and the traction between the drop and the substrate. Each self-propulsion mechanism gives rise to two main modes of motion: a long thin drop with zero traction in the bulk, mostly occurring under strong stretching forces, and a parabolic drop with finite traction in the bulk, mostly occurring under strong squeezing forces. In each case there is a sharp transition between parabolic, and long thin drops as a function of the applied forces and indications of drop break-up where large forces stretch the drop.

Long-term molecular turnover of actin stress fibers revealed by advection-reaction analysis in fluorescence recovery after photobleaching

Fluorescence recovery after photobleaching (FRAP) is a versatile technique to evaluate the intracellular molecular exchange called turnover. Mechanochemical models of FRAP typically consider the molecular diffusion and chemical reaction that simultaneously occur on a time scale of seconds to minutes. Particularly for long-term measurements, however, a mechanical advection effect can no longer be ignored, which transports the proteins in specific directions within the cells and accordingly shifts the spatial distribution of the local chemical equilibrium. Nevertheless, existing FRAP models have not considered the spatial shift, and as such, the turnover rate is often analyzed without considering the spatiotemporally updated chemical equilibrium. Here we develop a new FRAP model aimed at long-term measurements to quantitatively determine the two distinct effects of the advection and chemical reaction, i.e., the different major sources of the change in fluorescence intensity. To validate this approach, we carried out FRAP experiments on actin in stress fibers over a time period of more than 900 s, and the advection rate was shown to be comparable in magnitude to the chemical dissociation rate. We further found that the actin–myosin interaction and actin polymerization differently affect the advection and chemical dissociation. Our results suggest that the distinction between the two effects is indispensable to extract the intrinsic chemical properties of the actin cytoskeleton from the observations of complicated turnover in cells.

Ray Optics for Gliders

Control of self-propelled particles is central to the development of many microrobotic technologies, from dynamically reconfigurable materials to advanced lab-on-a-chip systems. However, there are few physical principles by which particle trajectories can be specified and can be used to generate a wide range of behaviors. Within the field of ray optics, a single principle for controlling the trajectory of light─Snell's law─yields an intuitive framework for engineering a broad range of devices, from microscopes to cameras and telescopes. Here we show that the motion of self-propelled particles gliding across a resistance discontinuity is governed by a variant of Snell's law, and develop a corresponding ray optics for gliders. Just as the ratio of refractive indexes sets the path of a light ray, the ratio of resistance coefficients is shown to determine the trajectories of gliders. The magnitude of refraction depends on the glider's shape, in particular its aspect ratio, which serves as an analogue to the wavelength of light. This enables the demixing of a polymorphic, many-shaped, beam of gliders into distinct monomorphic, single-shaped, beams through a friction prism. In turn, beams of monomorphic gliders can be focused by spherical and gradient friction lenses. Alternatively, the critical angle for total internal reflection can be used to create shape-selective glider traps. Overall our work suggests that furthering the analogy between light and microscopic gliders may be used for sorting, concentrating, and analyzing self-propelled particles.

Integrative experimental/computational approach establishes active cellular protrusion as the primary driving force of phagocytic spreading by immune cells

The dynamic interplay between cell adhesion and protrusion is a critical determinant of many forms of cell motility. When modeling cell spreading on adhesive surfaces, traditional mathematical treatments often consider passive cell adhesion as the primary, if not exclusive, mechanistic driving force of this cellular motion. To better assess the contribution of active cytoskeletal protrusion to immune-cell spreading during phagocytosis, we here develop a computational framework that allows us to optionally investigate both purely adhesive spreading (“Brownian zipper hypothesis”) as well as protrusion-dominated spreading (“protrusive zipper hypothesis”). We model the cell as an axisymmetric body of highly viscous fluid surrounded by a cortex with uniform surface tension and incorporate as potential driving forces of cell spreading an attractive stress due to receptor-ligand binding and an outward normal stress representing cytoskeletal protrusion, both acting on the cell boundary. We leverage various model predictions against the results of a directly related experimental companion study of human neutrophil phagocytic spreading on substrates coated with different densities of antibodies. We find that the concept of adhesion-driven spreading is incompatible with experimental results such as the independence of the cell-spreading speed on the density of immobilized antibodies. In contrast, the protrusive zipper model agrees well with experimental findings and, when adapted to simulate cell spreading on discrete adhesion sites, it also reproduces the observed positive correlation between antibody density and maximum cell-substrate contact area. Together, our integrative experimental/computational approach shows that phagocytic spreading is driven by cellular protrusion, and that the extent of spreading is limited by the density of adhesion sites.

To accomplish many routine biological tasks, cells must rapidly spread over different types of surfaces. Here, we examine the biophysical underpinnings of immune cell spreading during phagocytosis, the process by which white blood cells such as neutrophils engulf pathogens or other foreign objects. Our computational framework models the case in which a human neutrophil spreads over a flat surface coated with antibodies, which we also test experimentally in a companion paper. Our primary purpose is to assess whether phagocytic spreading is actively driven by protrusive forces exerted by the cell, or passively by adhesive forces acting between receptors in the cell membrane and antibodies on the surface. By directly comparing our model predictions to experimental results, we demonstrate that phagocytic spreading is primarily driven by protrusion, but the extent of spreading is still limited by the availability of binding sites. Our findings improve the fundamental understanding of phagocytosis and may also pave the way for future investigations of the balance between adhesion and protrusion in other forms of cell spreading, such as wound healing or cancer cell migration.

3D single cell migration driven by temporal correlation between oscillating force dipoles

Directional cell locomotion requires symmetry breaking between the front and rear of the cell. In some cells, symmetry breaking manifests itself in a directional flow of actin from the front to the rear of the cell. Many cells, especially in physiological 3D matrices do not show such coherent actin dynamics and present seemingly competing protrusion/retraction dynamics at their front and back. How symmetry breaking manifests itself for such cells is therefore elusive. We take inspiration from the scallop theorem proposed by Purcell for micro-swimmers in Newtonian fluids: self-propelled objects undergoing persistent motion at low Reynolds number must follow a cycle of shape changes that breaks temporal symmetry. We report similar observations for cells crawling in 3D. We quantified cell motion using a combination of 3D live cell imaging, visualization of the matrix displacement and a minimal model with multipolar expansion. We show that our cells embedded in a 3D matrix form myosin-driven force dipoles at both sides of the nucleus, that locally and periodically pinch the matrix. The existence of a phase shift between the two dipoles is required for directed cell motion which manifests itself as cycles with finite area in the dipole-quadrupole diagram, a formal equivalence to the Purcell cycle. We confirm this mechanism by triggering local dipolar contractions with a laser. This leads to directed motion. Our study reveals that these cells control their motility by synchronizing dipolar forces distributed at front and back. This result opens new strategies to externally control cell motion as well as for the design of micro-crawlers.

Cell motility as an energy minimization process

The dynamics of active matter driven by interacting molecular motors has a nonpotential structure at the local scale. However, we show that there exists a quasipotential effectively describing the collective self-organization of the motors propelling a cell at a continuum active gel level. Such a model allows us to understand cell motility as an active phase transition problem between the static and motile steady-state configurations that minimize the quasipotential. In particular, both configurations can coexist in a metastable fashion and a small stochastic disorder in the gel is sufficient to trigger an intermittent cell dynamics where either static or motile phases are more probable, depending on which state is the global minimum of the quasipotential.

Phase field model for cell spreading dynamics

We suggest a 3D phase field model to describe 3D cell spreading on a flat substrate. The model is a simplified version of a minimal model that was developed in Winkler (Commun Phys 2:82, 2019). Our model couples the order parameter u with 3D polarization (orientation) vector field [Formula: see text] of the actin network. We derive a closed integro-differential equation governing the 3D cell spreading dynamics on a flat substrate, which includes the normal velocity of the membrane, curvature, volume relaxation rate, a function determined by the molecular effects of the subcell level, and the adhesion effect. This equation is easily solved numerically. The results are in agreement with the early fast phase observed experimentally in Dobereiner (Phys Rev Lett 93:108105, 2004). Also we find agreement with the universal power law (Cuvelier in Curr Biol 17:694-699, 2007) which suggest that cell adhesion or contact area versus time behave as [Formula: see text] in the early stage of cell spreading dynamics, and slow down at the next stages.

Modeling cell protrusion predicts how myosin II and actin turnover affect adhesion-based signaling

Orchestration of cell migration is essential for development, tissue regeneration, and the immune response. This dynamic process integrates adhesion, signaling, and cytoskeletal subprocesses across spatial and temporal scales. In mesenchymal cells, adhesion complexes bound to extracellular matrix mediate both biochemical signal transduction and physical interaction with the F-actin cytoskeleton. Here, we present a mathematical model that offers insight into both aspects, considering spatiotemporal dynamics of nascent adhesions, active signaling molecules, mechanical clutching, actin treadmilling, and nonmuscle myosin II contractility. At the core of the model is a positive feedback loop, whereby adhesion-based signaling promotes generation of barbed ends at, and protrusion of, the cell's leading edge, which in turn promotes formation and stabilization of nascent adhesions. The model predicts a switch-like transition and optimality of membrane protrusion, determined by the balance of actin polymerization and retrograde flow, with respect to extracellular matrix density. The model, together with new experimental measurements, explains how protrusion can be modulated by mechanical effects (nonmuscle myosin II contractility and adhesive bond stiffness) and F-actin turnover.

Coupling traction force patterns and actomyosin wave dynamics reveals mechanics of cell motion

Motile cells can use and switch between different modes of migration. Here, we use traction force microscopy and fluorescent labeling of actin and myosin to quantify and correlate traction force patterns and cytoskeletal distributions in Dictyostelium discoideum cells that move and switch between keratocyte-like fan-shaped, oscillatory, and amoeboid modes. We find that the wave dynamics of the cytoskeletal components critically determine the traction force pattern, cell morphology, and migration mode. Furthermore, we find that fan-shaped cells can exhibit two different propulsion mechanisms, each with a distinct traction force pattern. Finally, the traction force patterns can be recapitulated using a computational model, which uses the experimentally determined spatiotemporal distributions of actin and myosin forces and a viscous cytoskeletal network. Our results suggest that cell motion can be generated by friction between the flow of this network and the substrate.

Optogenetic control of intracellular flows and cell migration: A comprehensive mathematical analysis with a minimal active gel model

The actin cytoskeleton of cells is in continuous motion due to both polymerization of new filaments and their contraction by myosin II molecular motors. Through adhesion to the substrate, such intracellular flow can be converted into cell migration. Recently, optogenetics has emerged as a new powerful experimental method to control both actin polymerization and myosin II contraction. While optogenetic control of polymerization can initiate cell migration by generating protrusion, it is less clear if and how optogenetic control of contraction can also affect cell migration. Here we analyze the latter situation using a minimal variant of active gel theory into which we include optogenetic activation as a spatiotemporally constrained perturbation. The model can describe the symmetrical flow of the actomyosin system observed in optogenetic experiments, but not the long-lasting polarization required for cell migration. Motile solutions become possible if cytoskeletal polymerization is included through the boundary conditions. Optogenetic activation of contraction can then initiate locomotion in a symmetrically spreading cell and strengthen motility in an asymmetrically polymerizing one. If designed appropriately, it can also arrest motility even for protrusive boundaries.

Analysis of protrusion dynamics in amoeboid cell motility by means of regularized contour flows

Amoeboid cell motility is essential for a wide range of biological processes including wound healing, embryonic morphogenesis, and cancer metastasis. It relies on complex dynamical patterns of cell shape changes that pose long-standing challenges to mathematical modeling and raise a need for automated and reproducible approaches to extract quantitative morphological features from image sequences. Here, we introduce a theoretical framework and a computational method for obtaining smooth representations of the spatiotemporal contour dynamics from stacks of segmented microscopy images. Based on a Gaussian process regression we propose a one-parameter family of regularized contour flows that allows us to continuously track reference points (virtual markers) between successive cell contours. We use this approach to define a coordinate system on the moving cell boundary and to represent different local geometric quantities in this frame of reference. In particular, we introduce the local marker dispersion as a measure to identify localized membrane expansions and provide a fully automated way to extract the properties of such expansions, including their area and growth time. The methods are available as an open-source software package called AmoePy, a Python-based toolbox for analyzing amoeboid cell motility (based on time-lapse microscopy data), including a graphical user interface and detailed documentation. Due to the mathematical rigor of our framework, we envision it to be of use for the development of novel cell motility models. We mainly use experimental data of the social amoeba Dictyostelium discoideum to illustrate and validate our approach.

Amoeboid motion is a crawling-like cell migration that plays an important key role in multiple biological processes such as wound healing and cancer metastasis. This type of cell motility results from expanding and simultaneously contracting parts of the cell membrane. From fluorescence images, we obtain a sequence of points, representing the cell membrane, for each time step. By using regression analysis on these sequences, we derive smooth representations, so-called contours, of the membrane. Since the number of measurements is discrete and often limited, the question is raised of how to link consecutive contours with each other. In this work, we present a novel mathematical framework in which these links are described by regularized flows allowing a certain degree of concentration or stretching of neighboring reference points on the same contour. This stretching rate, the so-called local dispersion, is used to identify expansions and contractions of the cell membrane providing a fully automated way of extracting properties of these cell shape changes. We applied our methods to time-lapse microscopy data of the social amoeba Dictyostelium discoideum.

Modeling spatiotemporally varying protein-protein interactions in CyLaKS, the Cytoskeleton Lattice-based Kinetic Simulator

Interaction of cytoskeletal filaments, motor proteins, and crosslinking proteins drives important cellular processes such as cell division and cell movement. Cytoskeletal networks also exhibit nonequilibrium self-assembly in reconstituted systems. An emerging problem in cytoskeletal modeling and simulation is spatiotemporal alteration of the dynamics of filaments, motors, and associated proteins. This can occur due to motor crowding, obstacles along the filament, motor interactions and direction switching, and changes, defects, or heterogeneity in the filament binding lattice. How such spatiotemporally varying cytoskeletal filaments and motor interactions affect their collective properties is not fully understood. We developed the Cytoskeleton Lattice-based Kinetic Simulator (CyLaKS) to investigate such problems. The simulation model builds on previous work by incorporating motor mechanochemistry into a simulation with many interacting motors and/or associated proteins on a discretized lattice. CyLaKS also includes detailed balance in binding kinetics, movement, and lattice heterogeneity. The simulation framework is flexible and extensible for future modeling work and is available on GitHub for others to freely use or build upon. Here we illustrate the use of CyLaKS to study long-range motor interactions, microtubule lattice heterogeneity, motion of a heterodimeric motor, and how changing crosslinker number affects filament separation.

Nanotopography-Modulated Epithelial Cell Collective Migration

Matrix nanotopography plays an essential role in regulating cell behaviors including cell proliferation, differentiation, and migration. While studies on isolated single cell migration along the nanostructural orientation have been reported for various cell types, there remains a lack of understanding of how nanotopography regulates the behavior of collectively migrating cells during processes such as epithelial wound healing. We demonstrated that collective migration of epithelial cells was promoted on nanogratings perpendicular to, but not on those parallel to, the wound-healing axis. We further discovered that nanograting-modulated epithelial migration was dominated by the adhesion turnover process, which was Rho-associated protein kinase activity-dependent, and the lamellipodia protrusion at the cell leading edge, which was Rac1-GTPase activity-dependent. This work provides explanations to the distinct migration behavior of epithelial cells on nanogratings, and indicates that the effect of nanotopographic modulations on cell migration is cell-type dependent and involves complex mechanisms.

A phosphoinositide-based model of actin waves in frustrated phagocytosis

Phagocytosis is a complex process by which phagocytes such as lymphocytes or macrophages engulf and destroy foreign bodies called pathogens in a tissue. The process is triggered by the detection of antibodies that trigger signaling mechanisms that control the changes of the cellular cytoskeleton needed for engulfment of the pathogen. A mathematical model of the entire process would be extremely complicated, because the signaling and cytoskeletal changes produce large mechanical deformations of the cell. Recent experiments have used a confinement technique that leads to a process called frustrated phagocytosis, in which the membrane does not deform, but rather, signaling triggers actin waves that propagate along the boundary of the cell. This eliminates the large-scale deformations and facilitates modeling of the wave dynamics. Herein we develop a model of the actin dynamics observed in frustrated phagocytosis and show that it can replicate the experimental observations. We identify the key components that control the actin waves and make a number of experimentally-testable predictions. In particular, we predict that diffusion coefficients of membrane-bound species must be larger behind the wavefront to replicate the internal structure of the waves. Our model is a first step toward a more complete model of phagocytosis, and provides insights into circular dorsal ruffles as well.

A formalism for modelling traction forces and cell shape evolution during cell migration in various biomedical processes

The phenomenological model for cell shape deformation and cell migration Chen (BMM 17:1429–1450, 2018), Vermolen and Gefen (BMM 12:301–323, 2012), is extended with the incorporation of cell traction forces and the evolution of cell equilibrium shapes as a result of cell differentiation. Plastic deformations of the extracellular matrix are modelled using morphoelasticity theory. The resulting partial differential differential equations are solved by the use of the finite element method. The paper treats various biological scenarios that entail cell migration and cell shape evolution. The experimental observations in Mak et al. (LC 13:340–348, 2013), where transmigration of cancer cells through narrow apertures is studied, are reproduced using a Monte Carlo framework.

Dynamic Nuclear Structure Emerges from Chromatin Cross-Links and Motors

The cell nucleus houses the chromosomes, which are linked to a soft shell of lamin protein filaments. Experiments indicate that correlated chromosome dynamics and nuclear shape fluctuations arise from motor activity. To identify the physical mechanisms, we develop a model of an active, cross-linked Rouse chain bound to a polymeric shell. System-sized correlated motions occur but require both motor activity and cross-links. Contractile motors, in particular, enhance chromosome dynamics by driving anomalous density fluctuations. Nuclear shape fluctuations depend on motor strength, cross-linking, and chromosome-lamina binding. Therefore, complex chromosome dynamics and nuclear shape emerge from a minimal, active chromosome-lamina system.

Comparison of explicit and mean-field models of cytoskeletal filaments with crosslinking motors

In cells, cytoskeletal filament networks are responsible for cell movement, growth, and division. Filaments in the cytoskeleton are driven and organized by crosslinking molecular motors. In reconstituted cytoskeletal systems, motor activity is responsible for far-from-equilibrium phenomena such as active stress, self-organized flow, and spontaneous nematic defect generation. How microscopic interactions between motors and filaments lead to larger-scale dynamics remains incompletely understood. To build from motor-filament interactions to predict bulk behavior of cytoskeletal systems, more computationally efficient techniques for modeling motor-filament interactions are needed. Here, we derive a coarse-graining hierarchy of explicit and continuum models for crosslinking motors that bind to and walk on filament pairs. We compare the steady-state motor distribution and motor-induced filament motion for the different models and analyze their computational cost. All three models agree well in the limit of fast motor binding kinetics. Evolving a truncated moment expansion of motor density speeds the computation by [Formula: see text]-[Formula: see text] compared to the explicit or continuous-density simulations, suggesting an approach for more efficient simulation of large networks. These tools facilitate further study of motor-filament networks on micrometer to millimeter length scales.

Collective migration during a gap closure in a two-dimensional haptotactic model

The ability of cells to respond to substrate-bound protein gradients is crucial for many physiological processes, such as immune response, neurogenesis and cancer cell migration. However, the difficulty to produce well-controlled protein gradients has long been a limitation to our understanding of collective cell migration in response to haptotaxis. Here we use a photopatterning technique to create circular, square and linear fibronectin (FN) gradients on two-dimensional (2D) culture substrates. We observed that epithelial cells spread preferentially on zones of higher FN density, creating rounded or elongated gaps within epithelial tissues over circular or linear FN gradients, respectively. Using time-lapse experiments, we demonstrated that the gap closure mechanism in a 2D haptotaxis model requires a significant increase of the leader cell area. In addition, we found that gap closures are slower on decreasing FN densities than on homogenous FN-coated substrate and that fresh closed gaps are characterized by a lower cell density. Interestingly, our results showed that cell proliferation increases in the closed gap region after maturation to restore the cell density, but that cell–cell adhesive junctions remain weaker in scarred epithelial zones. Taken together, our findings provide a better understanding of the wound healing process over protein gradients, which are reminiscent of haptotaxis.

Runx3 Induces a Cell Shape Change and Suppresses Migration and Metastasis of Melanoma Cells by Altering a Transcriptional Profile

Runt-related transcription factor-3 (Runx3) is a tumor suppressor, and its contribution to melanoma progression remains unclear. We previously demonstrated that Runx3 re-expression in B16-F10 melanoma cells changed their shape and attenuated their migration. In this study, we found that Runx3 re-expression in B16-F10 cells also suppressed their pulmonary metastasis. We performed microarray analysis and uncovered an altered transcriptional profile underlying the cell shape change and the suppression of migration and metastasis. This altered transcriptional profile was rich in Gene Ontology/Kyoto Encyclopedia of Genes and Genomes (GO/KEGG) annotations relevant to adhesion and the actin cytoskeleton and included differentially expressed genes for some major extracellular matrix (ECM) proteins as well as genes that were inversely associated with the increase in the metastatic potential of B16-F10 cells compared to B16-F0 melanoma cells. Further, we found that this altered transcriptional profile could have prognostic value, as evidenced by myelin and lymphocyte protein (MAL) and vilin-like (VILL). Finally, Mal gene expression was correlated with metastatic potential among the cells and was targeted by histone deacetylase (HDAC) inhibitors in B16-F10 cells, and the knockdown of Mal gene expression in B16-F0 cells changed their shape and enhanced the migratory and invasive traits of their metastasis. Our study suggests that self-entrapping of metastatic Runx3-negative melanoma cells via adhesion and the actin cytoskeleton could be a powerful therapeutic strategy.

Elastic wrinkling of keratocyte lamellipodia driven by myosin-induced contractile stress

During actin-based cell migration, the actin cytoskeleton in the lamellipodium both generates and responds to force, which has functional consequences for the ability of the cell to extend protrusions. However, the material properties of the lamellipodial actin network and its response to stress on the timescale of motility are incompletely understood. Here, we describe a dynamic wrinkling phenotype in the lamellipodium of fish keratocytes, in which the actin sheet buckles upward away from the ventral membrane of the cell, forming a periodic pattern of wrinkles perpendicular to the cell's leading edge. Cells maintain an approximately constant wrinkle wavelength over time despite new wrinkle formation and the lateral movement of wrinkles in the cell frame of reference, suggesting that cells have a preferred or characteristic wrinkle wavelength. Generation of wrinkles is dependent upon myosin contractility, and their wavelength scales directly with the density of the actin network and inversely with cell adhesion. These results are consistent with a simple physical model for wrinkling in an elastic sheet under compression and suggest that the lamellipodial cytoskeleton behaves as an elastic material on the timescale of cell migration despite rapid actin turnover.

On the adhesion-velocity relation and length adaptation of motile cells on stepped fibronectin lanes

The biphasic adhesion-velocity relation is a universal observation in mesenchymal cell motility. It has been explained by adhesion-promoted forces pushing the front and resisting motion at the rear. Yet, there is little quantitative understanding of how these forces control cell velocity. We study motion of MDA-MB-231 cells on microlanes with fields of alternating Fibronectin densities to address this topic and derive a mathematical model from the leading-edge force balance and the force-dependent polymerization rate. It reproduces quantitatively our measured adhesion-velocity relation and results with keratocytes, PtK1 cells, and CHO cells. Our results confirm that the force pushing the leading-edge membrane drives lamellipodial retrograde flow. Forces resisting motion originate along the whole cell length. All motion-related forces are controlled by adhesion and velocity, which allows motion, even with higher Fibronectin density at the rear than at the front. We find the pathway from Fibronectin density to adhesion structures to involve strong positive feedbacks. Suppressing myosin activity reduces the positive feedback. At transitions between different Fibronectin densities, steady motion is perturbed and leads to changes of cell length and front and rear velocity. Cells exhibit an intrinsic length set by adhesion strength, which, together with the length dynamics, suggests a spring-like front-rear interaction force. We provide a quantitative mechanistic picture of the adhesion-velocity relation and cell response to adhesion changes integrating force-dependent polymerization, retrograde flow, positive feedback from integrin to adhesion structures, and spring-like front-rear interaction.

Green algae scatter off sharp viscosity gradients

We study the behaviour of the green alga Chlamydomonas reinhardtii (CR) in the presence of neighbouring regions of different viscosity. We show that the velocity and angular diffusion of the algae decreases when the viscosity of the surrounding medium is increased. We report on a phenomenon occurring when the algae try to cross from a region of low viscosity to a highly viscous one, which causes CR to re-orient and scatter away from the interface if it is approached at a sufficiently small angle. We highlight that the effect does not occur for CR crossing from high to low viscosity regions. Lastly we show that algae do not concentrate in the region of high viscosity despite them swimming slower there. On the contrary, they concentrate in the region of low viscosity or maintain a uniform concentration profile, depending on the viscosity ratio between the two regions.

Dynamics of a cell motility model near the sharp interface limit

Phase-field models have recently had great success in describing the dynamic morphologies and motility of eukaryotic cells. In this work we investigate the minimal phase-field model introduced in Berlyand et al. (2017). Rigorous analysis of its sharp interface limit dynamics was completed in Mizuhara et al. (2016) and Mizuhara et al. (2019), where it was observed that persistent cell motion was not stable. In this work we numerically study the pre-limiting phase-field model near the sharp interface limit, to better understand this lack of persistent motion. We find that immobile, persistent, and rotating states are all exhibited in this minimal model, and investigate the loss of persistent motion in the sharp interface limit.

Deciphering the dynamics of lamellipodium in a fish keratocytes model

Motile cells depend on an intricate network of feedback loops that are essential in driving cell movement. Integrin-based focal adhesions (FAs) along with actin are the two key factors that mediate such motile behaviour. Together, they generate excitable dynamics that are essential for forming protrusions at the leading edge of the cell and, in certain cases, traveling waves along the membrane. A partial differential equation (PDE) model of a self-organizing lamellipodium in crawling keratocytes has been previously developed to understand how the three spatiotemporal patterns of activity observed in such cells, namely, stalling, waving and smooth motility, are produced. The model consisted of three key variables: the density of barbed actin filaments, newly formed FAs called nascent adhesions (NAs) and VASP, an anti-capping protein that gets sequestered by NAs during maturation. Using parameter sweeping techniques, the distinct regimes of behaviour associated with the three activity patterns were identified. In this study, we convert the PDE model into an ordinary differential equation (ODE) model to examine its excitability properties and determine all the patterns of activity exhibited by this system. Our results reveal that there are two additional regimes not previously identified, including bistability and oscillatory-like type IV excitability (generated by three steady states and their manifolds, rather than limit cycles). These regimes are also present in the PDE model. Applying slow-fast analysis on the ODE model shows that it exhibits a canard explosion through a folded-saddle and that rough motility seen in keratocytes is likely due to noise-dependent motility governed by dynamics near the interface of bistability and type IV excitability. The two parameter bifurcation suggests that the increase in the proportion of rough motion is due to a shift in activity towards the bistable and type IV excitable regimes induced by a decrease in NA maturation rate. Our results thus provide important insight into how microscopic mechanical effects are integrated to produce the observed modes of motility.

Stick-slip model for actin-driven cell protrusions, cell polarization, and crawling

Cell crawling requires the generation of intracellular forces by the cytoskeleton and their transmission to an extracellular substrate through specific adhesion molecules. Crawling cells show many features of excitable systems, such as spontaneous symmetry breaking and crawling in the absence of external cues, and periodic and propagating waves of activity. Mechanical instabilities in the active cytoskeleton network and feedback loops in the biochemical network of activators and repressors of cytoskeleton dynamics have been invoked to explain these dynamical features. Here, I show that the interplay between the dynamics of cell-substrate adhesion and linear cellular mechanics is sufficient to reproduce many nonlinear dynamical patterns observed in spreading and crawling cells. Using an analytical formalism of the molecular clutch model of cell adhesion, regulated by local mechanical forces, I show that cellular traction forces exhibit stick-slip dynamics resulting in periodic waves of protrusion/retraction and propagating waves along the cell edge. This can explain spontaneous symmetry breaking and polarization of spreading cells, leading to steady crawling or bipedal motion, and bistability, where persistent cell motion requires a sufficiently strong transient external stimulus. The model also highlights the role of membrane tension in providing the long-range mechanical communication across the cell required for symmetry breaking.

Cell-substrate mechanics guide collective cell migration through intercellular adhesion: a dynamic finite element cellular model

During the process of tissue formation and regeneration, cells migrate collectively while remaining connected through intercellular adhesions. However, the roles of cell-substrate and cell-cell mechanical interactions in regulating collective cell migration are still unclear. In this study, we employ a newly developed finite element cellular model to study collective cell migration by exploring the effects of mechanical feedback between cell and substrate and mechanical signal transmission between adjacent cells. Our viscoelastic model of cells consists many triangular elements and is of high resolution. Cadherin adhesion between cells is modeled explicitly as linear springs at subcellular level. In addition, we incorporate a mechano-chemical feedback loop between cell-substrate mechanics and Rac-mediated cell protrusion. Our model can reproduce a number of experimentally observed patterns of collective cell migration during wound healing, including cell migration persistence, separation distance between cell pairs and migration direction. Moreover, we demonstrate that cell protrusion determined by the cell-substrate mechanics plays an important role in guiding persistent and oriented collective cell migration. Furthermore, this guidance cue can be maintained and transmitted to submarginal cells of long distance through intercellular adhesions. Our study illustrates that our finite element cellular model can be employed to study broad problems of complex tissue in dynamic changes at subcellular level.

Cell Mechanics at the Rear Act to Steer the Direction of Cell Migration

Motile cells navigate complex environments by changing their direction of travel, generating left-right asymmetries in their mechanical subsystems to physically turn. Currently, little is known about how external directional cues are propagated along the length scale of the whole cell and integrated with its force-generating apparatus to steer migration mechanically. We examine the mechanics of spontaneous cell turning in fish epidermal keratocytes and find that the mechanical asymmetries responsible for turning behavior predominate at the rear of the cell, where there is asymmetric centripetal actin flow. Using experimental perturbations, we identify two linked feedback loops connecting myosin II contractility, adhesion strength and actin network flow in turning cells that are sufficient to explain the observed cell shapes and trajectories. Notably, asymmetries in actin polymerization at the cell leading edge play only a minor role in the mechanics of cell turning-that is, cells steer from the rear.

Cellular self-organization: generating order from the abyss

While the organization of inanimate systems such as gases or liquids is predominantly thermodynamically driven—a mixture of two gases will tend to mix until they reach equilibrium—biological systems frequently exhibit organization that is far from a well-mixed equilibrium. The anisotropies displayed by cells are evident in some of the dynamic processes that constitute life including cell development, movement, and division. These anisotropies operate at different length-scales, from the meso- to the nanoscale, and are proposed to reflect self-organization, a characteristic of living systems that is becoming accessible to reconstitution from purified components, and thus a more thorough understanding. Here, some examples of self-organization underlying cellular anisotropies at the cellular level are reviewed, with an emphasis on Rho-family GTPases operating at the plasma membrane. Given the technical challenges of studying these dynamic proteins, some of the successful approaches that are being employed to study their self-organization will also be considered.

Actomyosin Contraction Induces In-Bulk Motility of Cells and Droplets

Cell crawling on two-dimensional surfaces is a relatively well-understood phenomenon that is based on actin polymerization at a cell's front edge and anchoring on a substrate, allowing the cell to pull itself forward. However, some cells, such as cancer cells invading a three-dimensional matrigel, can also swim in the bulk, where surface adhesion is impossible. Although there is strong evidence that the self-organized engine that drives cells forward in the bulk involves myosin, the specific propulsion mechanism remains largely unclear. Here, we propose a minimal model for in-bulk self-motility of a droplet containing an isotropic and compressible contractile gel, representing a cell extract containing a disordered actomyosin network. In our model, contraction mediates a feedback loop between myosin-induced flow and advection-induced myosin accumulation, which leads to clustering and locally enhanced flow. The symmetry of such flow is then spontaneously broken through actomyosin-membrane interactions, leading to self-organized droplet motility relative to the underlying solvent. Depending on the balance between contraction, diffusion, detachment rate of myosin, and effective surface tension, this motion can be either straight or circular. Our simulations and analytical results shed new light on in-bulk myosin-driven cell motility in living cells and provide a framework to design a novel type of synthetic active matter droplet potentially resembling the motility mechanism of biological cells.