Actin-myosin viscoelastic flow in the keratocyte lamellipod

A simple active fluid model unites cytokinesis, cell crawling, and axonal outgrowth

Axonal outgrowth, cell crawling, and cytokinesis utilize actomyosin, microtubule-based motors, cytoskeletal dynamics, and substrate adhesions to produce traction forces and bulk cellular motion. While it has long been appreciated that growth cones resemble crawling cells and that the mechanisms that drive cytokinesis help power cell crawling, they are typically viewed as unique processes. To better understand the relationship between these modes of motility, here, we developed a unified active fluid model of cytokinesis, amoeboid migration, mesenchymal migration, neuronal migration, and axonal outgrowth in terms of cytoskeletal flow, adhesions, viscosity, and force generation. Using numerical modeling, we fit subcellular velocity profiles of the motions of cytoskeletal structures and docked organelles from previously published studies to infer underlying patterns of force generation and adhesion. Our results indicate that, during cytokinesis, there is a primary converge zone at the cleavage furrow that drives flow towards it; adhesions are symmetric across the cell, and as a result, cells are stationary. In mesenchymal, amoeboid, and neuronal migration, the site of the converge zone shifts, and differences in adhesion between the front and back of the cell drive crawling. During neuronal migration and axonal outgrowth, the primary convergence zone lies within the growth cone, which drives actin retrograde flow in the P-domain and bulk anterograde flow of the axonal shaft. They differ in that during neuronal migration, the cell body is weakly attached to the substrate and thus moves forward at the same velocity as the axon. In contrast, during axonal outgrowth, the cell body strongly adheres to the substrate and remains stationary, resulting in a decrease in flow velocity away from the growth cone. The simplicity with which cytokinesis, cell crawling, and axonal outgrowth can be modeled by varying coefficients in a simple model suggests a deep connection between them.

Biomechanics of parasite migration within hosts

The dissemination of protozoan and metazoan parasites through host tissues is hindered by cellular barriers, dense extracellular matrices, and fluid forces in the bloodstream. To overcome these diverse biophysical impediments, parasites implement versatile migratory strategies. Parasite-exerted mechanical forces and upregulation of the host's cellular contractile machinery are the motors for these strategies, and these are comparably better characterized for protozoa than for helminths. Using the examples of the protozoans, Toxoplasma gondii and Plasmodium, and the metazoan, Schistosoma mansoni, we highlight how quantitative tools such as traction force and reflection interference contrast microscopies have improved our understanding of how parasites alter host mechanobiology to promote their migration.

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 laminin-keratin link shields the nucleus from mechanical deformation and signalling

The mechanical properties of the extracellular matrix dictate tissue behaviour. In epithelial tissues, laminin is a very abundant extracellular matrix component and a key supporting element. Here we show that laminin hinders the mechanoresponses of breast epithelial cells by shielding the nucleus from mechanical deformation. Coating substrates with laminin-111-unlike fibronectin or collagen I-impairs cell response to substrate rigidity and YAP nuclear localization. Blocking the laminin-specific integrin β4 increases nuclear YAP ratios in a rigidity-dependent manner without affecting the cell forces or focal adhesions. By combining mechanical perturbations and mathematical modelling, we show that β4 integrins establish a mechanical linkage between the substrate and keratin cytoskeleton, which stiffens the network and shields the nucleus from actomyosin-mediated mechanical deformation. In turn, this affects the nuclear YAP mechanoresponses, chromatin methylation and cell invasion in three dimensions. Our results demonstrate a mechanism by which tissues can regulate their sensitivity to mechanical signals.

The Influence of Nucleus Mechanics in Modelling Adhesion-independent Cell Migration in Structured and Confined Environments

Recent biological experiments (Lämmermann et al. in Nature 453(7191):51-55, 2008; Reversat et al. in Nature 7813:582-585, 2020; Balzer et al. in ASEB J Off Publ Fed Am Soc Exp Biol 26(10):4045-4056, 2012) have shown that certain types of cells are able to move in structured and confined environments even without the activation of focal adhesion. Focusing on this particular phenomenon and based on previous works (Jankowiak et al. in Math Models Methods Appl Sci 30(03):513-537, 2020), we derive a novel two-dimensional mechanical model, which relies on the following physical ingredients: the asymmetrical renewal of the actin cortex supporting the membrane, resulting in a backward flow of material; the mechanical description of the nuclear membrane and the inner nuclear material; the microtubule network guiding nucleus location; the contact interactions between the cell and the external environment. The resulting fourth order system of partial differential equations is then solved numerically to conduct a study of the qualitative effects of the model parameters, mainly those governing the mechanical properties of the nucleus and the geometry of the confining structure. Coherently with biological observations, we find that cells characterized by a stiff nucleus are unable to migrate in channels that can be crossed by cells with a softer nucleus. Regarding the geometry, cell velocity and ability to migrate are influenced by the width of the channel and the wavelength of the external structure. Even though still preliminary, these results may be potentially useful in determining the physical limit of cell migration in confined environments and in designing scaffolds for tissue engineering.

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 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.

Self-organization in amoeboid motility

Amoeboid motility has come to refer to a spectrum of cell migration modes enabling a cell to move in the absence of strong, specific adhesion. To do so, cells have evolved a range of motile surface movements whose physical principles are now coming into view. In response to external cues, many cells—and some single-celled-organisms—have the capacity to turn off their default migration mode. and switch to an amoeboid mode. This implies a restructuring of the migration machinery at the cell scale and suggests a close link between cell polarization and migration mediated by self-organizing mechanisms. Here, I review recent theoretical models with the aim of providing an integrative, physical picture of amoeboid migration.

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.

Two Complementary Signaling Pathways Depict Eukaryotic Chemotaxis: A Mechanochemical Coupling Model

Many eukaryotic cells, including neutrophils and Dictyostelium cells, are able to undergo correlated random migration in the absence of directional cues while reacting to shallow gradients of chemoattractants with exquisite precision. Although progress has been made with regard to molecular identities, it remains elusive how molecular mechanics are integrated with cell mechanics to initiate and manipulate cell motility. Here, we propose a two dimensional (2D) cell migration model wherein a multilayered dynamic seesaw mechanism is accompanied by a mechanical strain-based inhibition mechanism. In biology, these two mechanisms can be mapped onto the biochemical feedback between phosphoinositides (PIs) and Rho GTPase and the mechanical interplay between filamin A (FLNa) and FilGAP. Cell migration and the accompanying morphological changes are demonstrated in numerical simulations using a particle-spring model, and the diffusion in the cell membrane are simulations using a one dimensional (1D) finite differences method (FDM). The fine balance established between endogenous signaling and a mechanically governed inactivation scheme ensures the endogenous cycle of self-organizing pseudopods, accounting for the correlated random migration. Furthermore, this model cell manifests directional and adaptable responses to shallow graded signaling, depending on the overwhelming effect of the graded stimuli guidance on strain-based inhibition. Finally, the model cell becomes trapped within an obstacle-ridden spatial region, manifesting a shuttle run for local explorations and can chemotactically “escape”, illustrating again the balance required in the complementary signaling pathways.

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.

Discrete mechanical model of lamellipodial actin network implements molecular clutch mechanism and generates arcs and microspikes

Mechanical forces, actin filament turnover, and adhesion to the extracellular environment regulate lamellipodial protrusions. Computational and mathematical models at the continuum level have been used to investigate the molecular clutch mechanism, calculating the stress profile through the lamellipodium and around focal adhesions. However, the forces and deformations of individual actin filaments have not been considered while interactions between actin networks and actin bundles is not easily accounted with such methods. We develop a filament-level model of a lamellipodial actin network undergoing retrograde flow using 3D Brownian dynamics. Retrograde flow is promoted in simulations by pushing forces from the leading edge (due to actin polymerization), pulling forces (due to molecular motors), and opposed by viscous drag in cytoplasm and focal adhesions. Simulated networks have densities similar to measurements in prior electron micrographs. Connectivity between individual actin segments is maintained by permanent and dynamic crosslinkers. Remodeling of the network occurs via the addition of single actin filaments near the leading edge and via filament bond severing. We investigated how several parameters affect the stress distribution, network deformation and retrograde flow speed. The model captures the decrease in retrograde flow upon increase of focal adhesion strength. The stress profile changes from compression to extension across the leading edge, with regions of filament bending around focal adhesions. The model reproduces the observed reduction in retrograde flow speed upon exposure to cytochalasin D, which halts actin polymerization. Changes in crosslinker concentration and dynamics, as well as in the orientation pattern of newly added filaments demonstrate the model's ability to generate bundles of filaments perpendicular (actin arcs) or parallel (microspikes) to the protruding direction.

Estimation Algorithm for a Hybrid PDE–ODE Model Inspired by Immunocompetent Cancer-on-Chip Experiment

The present work is motivated by the development of a mathematical model mimicking the mechanisms observed in lab-on-chip experiments, made to reproduce on microfluidic chips the in vivo reality. Here we consider the Cancer-on-Chip experiment where tumor cells are treated with chemotherapy drug and secrete chemical signals in the environment attracting multiple immune cell species. The in silico model here proposed goes towards the construction of a “digital twin” of the experimental immune cells in the chip environment to better understand the complex mechanisms of immunosurveillance. To this aim, we develop a tumor-immune microfluidic hybrid PDE–ODE model to describe the concentration of chemicals in the Cancer-on-Chip environment and immune cells migration. The development of a trustable simulation algorithm, able to reproduce the immunocompetent dynamics observed in the chip, requires an efficient tool for the calibration of the model parameters. In this respect, the present paper represents a first methodological work to test the feasibility and the soundness of the calibration technique here proposed, based on a multidimensional spline interpolation technique for the time-varying velocity field surfaces obtained from cell trajectories.

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.

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.

Tools for computational analysis of moving boundary problems in cellular mechanobiology

A cell's ability to change shape is one of the most fundamental biological processes and is essential for maintaining healthy organisms. When the ability to control shape goes awry, it often results in a diseased system. As such, it is important to understand the mechanisms that allow a cell to sense and respond to its environment so as to maintain cellular shape homeostasis. Because of the inherent complexity of the system, computational models that are based on sound theoretical understanding of the biochemistry and biomechanics and that use experimentally measured parameters are an essential tool. These models involve an inherent feedback, whereby shape is determined by the action of regulatory signals whose spatial distribution depends on the shape. To carry out computational simulations of these moving boundary problems requires special computational techniques. A variety of alternative approaches, depending on the type and scale of question being asked, have been used to simulate various biological processes, including cell motility, division, mechanosensation, and cell engulfment. In general, these models consider the forces that act on the system (both internally generated, or externally imposed) and the mechanical properties of the cell that resist these forces. Moving forward, making these techniques more accessible to the non-expert will help improve interdisciplinary research thereby providing new insight into important biological processes that affect human health. This article is categorized under: Cancer > Cancer>Computational Models Cancer > Cancer>Molecular and Cellular Physiology.

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.

Active random forces can drive differential cellular positioning and enhance motor-driven transport

Cells are remarkable machines capable of performing an exquisite range of functions, many of which depend crucially on the activity of molecular motors that generate forces. Recent experiments have shown that intracellular random movements are not solely thermal in nature but also arise from stochasticity in the forces from these molecular motors. Here we consider the effects of these nonthermal random forces. We show that stochastic motor force not only enhances diffusion but also leads to size-dependent transport of objects that depends on the local density of the cytoskeletal filaments on which motors operate. As a consequence, we find that objects that are larger than the mesh size of the cytoskeleton should be attracted to regions of high cytoskeletal density, while objects that are smaller than the mesh size will preferentially avoid these regions. These results suggest a mechanism for size-based organelle positioning and also suggest that motor-driven random forces can additionally enhance motor-driven transport.

Influence of myosin activity and mechanical impact on keratocyte polarization

In cell migration, polarization is the process by which a stationary cell breaks symmetry and initiates motion. Although a lot is known about the mechanisms involved in cell polarization, the role played by myosin contraction remains unclear. In addition, cell polarization by mechanical impact has received little attention. Here, we study the influence of myosin activity on cell polarization and the initiation of motion induced by mechanical cues using a computational model for keratocytes. The model accounts for cell deformation, the dynamics of myosin and the signaling protein RhoA (a member of the Rho GTPases family), as well as the forces acting on the actomyosin network. Our results show that the attainment of a steady polarized state depends on the strength of myosin down- or up-regulation and that myosin upregulation favors cell polarization. Our results also confirm the existence of a threshold level for cell polarization, which is determined by the level of polarization of the Rho GTPases at the time the external stimuli vanish. In all, this paper shows that capturing the interactions between the signaling proteins (Rho GTPases for keratocytes) and the compounds of the motile machinery in a moving cell is crucial to study cell polarization.

A minimal mechanosensing model predicts keratocyte evolution on flexible substrates

A mathematical model is proposed for shape evolution and locomotion of fish epidermal keratocytes on elastic substrates. The model is based on mechanosensing concepts: cells apply contractile forces onto the elastic substrate, while cell shape evolution depends locally on the substrate stress generated by themselves or external mechanical stimuli acting on the substrate. We use the level set method to study the behaviour of the model numerically, and predict a number of distinct phenomena observed in experiments, such as (i) symmetry breaking from the stationary centrosymmetric to the well-known steadily propagating crescent shape, (ii) asymmetric bipedal oscillations and travelling waves in the lamellipodium leading edge, (iii) response to remote mechanical stress externally applied to the substrate (tensotaxis) and (iv) changing direction of motion towards an interface with a rigid substrate (durotaxis).

Modeling the Shape of Synaptic Spines by Their Actin Dynamics

Dendritic spines are the morphological basis of excitatory synapses in the cortex and their size and shape correlates with functional synaptic properties. Recent experiments show that spines exhibit large shape fluctuations that are not related to activity-dependent plasticity but nonetheless might influence memory storage at their synapses. To investigate the determinants of such spontaneous fluctuations, we propose a mathematical model for the dynamics of the spine shape and analyze it in 2D-related to experimental microscopic imagery-and in 3D. We show that the spine shape is governed by a local imbalance between membrane tension and the expansive force from actin bundles that originates from discrete actin polymerization foci. Experiments have shown that only few such polymerization foci co-exist at any time in a spine, each having limited life time. The model shows that the momentarily existing set of such foci pushes the membrane along certain directions until foci are replaced and other directions may now be affected. We explore these relations in depth and use our model to predict shape and temporal characteristics of spines from the different biophysical parameters involved in actin polymerization. Approximating the model by a single recursive equation we finally demonstrate that the temporal evolution of the number of active foci is sufficient to predict the size of the model-spines. Thus, our model provides the first platform to study the relation between molecular and morphological properties of the spine with a high degree of biophysical detail.

Active gel segment behaving as an active particle

We reduce a one-dimensional model of an active segment (AS), which is used, for instance, in the description of contraction-driven cell motility, to a zero-dimensional model of an active particle (AP) characterized by two internal degrees of freedom: position and polarity. Both models give rise to hysteretic force-velocity relations showing that an active agent can support two opposite polarities under the same external force and that it can maintain the same polarity while being dragged by external forces with opposite orientations. This double bistability results in a rich dynamic repertoire which we illustrate by studying static, stalled, motile, and periodically repolarizing regimes displayed by an active agent confined in a viscoelastic environment. We show that the AS and AP models can be calibrated to generate quantitatively similar dynamic responses.

Flower-like patterns in multi-species bacterial colonies

Diverse interactions among species within bacterial colonies lead to intricate spatiotemporal dynamics, which can affect their growth and survival. Here, we describe the emergence of complex structures in a colony grown from mixtures of motile and non-motile bacterial species on a soft agar surface. Time-lapse imaging shows that non-motile bacteria 'hitchhike' on the motile bacteria as the latter migrate outward. The non-motile bacteria accumulate at the boundary of the colony and trigger an instability that leaves behind striking flower-like patterns. The mechanism of the front instability governing this pattern formation is elucidated by a mathematical model for the frictional motion of the colony interface, with friction depending on the local concentration of the non-motile species. A more elaborate two-dimensional phase-field model that explicitly accounts for the interplay between growth, mechanical stress from the motile species, and friction provided by the non-motile species, fully reproduces the observed flower-like patterns.

Dual role of the nucleus in cell migration on planar substrates

Cell migration is essential to sustain life. There have been significant advances in the understanding of the mechanisms that control cell crawling, but the role of the nucleus remains poorly understood. The nucleus exerts a tight control of cell migration in 3D environments, but its influence in 2D migration on planar substrates remains unclear. Here, we study the role of the cell nucleus in 2D cell migration using a computational model of fish keratocytes. Our results indicate that the apparently minor role played by the nucleus emerges from two antagonist effects: While the nucleus modifies the spatial distributions of actin and myosin in a way that reduces cell velocity (e.g., the nucleus displaces myosin to the sides and front of the cell), its mechanical connection with the cytoskeleton alters the intracellular stresses promoting cell migration. Overall, the favorable effect of the nucleus-cytoskeleton connection prevails, which may explain why regular cells usually move faster than enucleated cells.

Crawling in a Fluid

There is increasing evidence that mammalian cells not only crawl on substrates but can also swim in fluids. To elucidate the mechanisms of the onset of motility of cells in suspension, a model which couples actin and myosin kinetics to fluid flow is proposed and solved for a spherical shape. The swimming speed is extracted in terms of key parameters. We analytically find super- and subcritical bifurcations from a nonmotile to a motile state and also spontaneous polarity oscillations that arise from a Hopf bifurcation. Relaxing the spherical assumption, the obtained shapes show appealing trends.

Stick-slip dynamics of migrating cells on viscoelastic substrates

Stick-slip motion, a common phenomenon observed during crawling of cells, is found to be strongly sensitive to the substrate stiffness. Stick-slip behaviors have previously been investigated typically using purely elastic substrates. For a more realistic understanding of this phenomenon, we propose a theoretical model to study the dynamics on a viscoelastic substrate. Our model, based on a reaction-diffusion framework, incorporates known important interactions such as retrograde flow of actin, myosin contractility, force-dependent assembly, and disassembly of focal adhesions coupled with cell-substrate interaction. We show that consideration of a viscoelastic substrate not only captures the usually observed stick-slip jumps but also predicts the existence of an optimal substrate viscosity corresponding to maximum traction force and minimum retrograde flow which was hitherto unexplored. Moreover, our theory predicts the time evolution of individual bond force that characterizes the stick-slip patterns on soft versus stiff substrates. Our analysis also elucidates how the duration of the stick-slip cycles are affected by various cellular parameters.

Cell motility dependence on adhesive wetting

Adhesive cell-substrate interactions are crucial for cell motility and are responsible for the necessary traction that propels cells. These interactions can also change the shape of the cell, analogous to liquid droplet wetting on adhesive substrates. To address how these shape changes affect cell migration and cell speed we model motility using deformable, 2D cross-sections of cells in which adhesion and frictional forces between cell and substrate can be varied separately. Our simulations show that increasing the adhesion results in increased spreading of cells and larger cell speeds. We propose an analytical model which shows that the cell speed is inversely proportional to an effective height of the cell and that increasing this height results in increased internal shear stress. The numerical and analytical results are confirmed in experiments on motile eukaryotic cells.

A discrete in continuous mathematical model of cardiac progenitor cells formation and growth as spheroid clusters (Cardiospheres)

We propose a discrete in continuous mathematical model describing the in vitro growth process of biophsy-derived mammalian cardiac progenitor cells growing as clusters in the form of spheres (Cardiospheres). The approach is hybrid: discrete at cellular scale and continuous at molecular level. In the present model, cells are subject to the self-organizing collective dynamics mechanism and, additionally, they can proliferate and differentiate, also depending on stochastic processes. The two latter processes are triggered and regulated by chemical signals present in the environment. Numerical simulations show the structure and the development of the clustered progenitors and are in a good agreement with the results obtained from in vitro experiments.

Stability, Convergence, and Sensitivity Analysis of the FBLM and the Corresponding FEM

We study in this paper the filament-based lamellipodium model (FBLM) and the corresponding finite element method (FEM) used to solve it. We investigate fundamental numerical properties of the FEM and justify its further use with the FBLM. We show that the FEM satisfies a time step stability condition that is consistent with the nature of the problem and propose a particular strategy to automatically adapt the time step of the method. We show that the FEM converges with respect to the (two-dimensional) space discretization in a series of characteristic and representative chemotaxis and haptotaxis experiments. We embed and couple the FBLM with a complex and adaptive extracellular environment comprised of chemical and adhesion components that are described by their macroscopic density and study their combined time evolution. With this combination, we study the sensitivity of the FBLM on several of its controlling parameters and discuss their influence in the dynamics of the model and its future evolution. We finally perform a number of numerical experiments that reproduce biological cases and compare the results with the ones reported in the literature.

Strain energy storage and dissipation rate in active cell mechanics

When living cells are observed at rest on a flat substrate, they can typically exhibit a rounded (symmetric) or an elongated (polarized) shape. Although the cells are apparently at rest, the active stress generated by the molecular motors continuously stretches and drifts the actin network, the cytoskeleton of the cell. In this paper we theoretically compare the energy stored and dissipated in this active system in two geometric configurations of interest: symmetric and polarized. We find that the stored energy is larger for a radially symmetric cell at low activation regime, while the polar configuration has larger strain energy when the active stress is beyond a critical threshold. Conversely, the dissipation of energy in a symmetric cell is always larger than that of a nonsymmetric one. By a combination of symmetry arguments and competition between surface and bulk stress, we argue that radial symmetry is an energetically expensive metastable state that provides access to an infinite number of lower-energy states, the polarized configurations.

Computational modeling of single-cell mechanics and cytoskeletal mechanobiology

Cellular cytoskeletal mechanics plays a major role in many aspects of human health from organ development to wound healing, tissue homeostasis and cancer metastasis. We summarize the state-of-the-art techniques for mathematically modeling cellular stiffness and mechanics and the cytoskeletal components and factors that regulate them. We highlight key experiments that have assisted model parameterization and compare the advantages of different models that have been used to recapitulate these experiments. An overview of feed-forward mechanisms from signaling to cytoskeleton remodeling is provided, followed by a discussion of the rapidly growing niche of encapsulating feedback mechanisms from cytoskeletal and cell mechanics to signaling. We discuss broad areas of advancement that could accelerate research and understanding of cellular mechanobiology. A precise understanding of the molecular mechanisms that affect cell and tissue mechanics and function will underpin innovations in medical device technologies of the future. WIREs Syst Biol Med 2018, 10:e1407. doi: 10.1002/wsbm.1407 This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.

Mechanical stress as a regulator of cell motility

The motility of a cell can be triggered or inhibited not only by an applied force but also by a mechanically neutral force couple. This type of loading, represented by an applied stress and commonly interpreted as either squeezing or stretching, can originate from extrinsic interaction of a cell with its neighbors. To quantify the effect of applied stresses on cell motility we use an analytically transparent one-dimensional model accounting for active myosin contraction and induced actin turnover. We show that stretching can polarize static cells and initiate cell motility while squeezing can symmetrize and arrest moving cells. We show further that sufficiently strong squeezing can lead to the loss of cell integrity. The overall behavior of the system depends on the two dimensionless parameters characterizing internal driving (chemical activity) and external loading (applied stress). We construct a phase diagram in this parameter space distinguishing between static, motile, and collapsed states. The obtained results are relevant for the mechanical understanding of contact inhibition and the epithelial-to-mesenchymal transition.

A free-boundary model of a motile cell explains turning behavior

To understand shapes and movements of cells undergoing lamellipodial motility, we systematically explore minimal free-boundary models of actin-myosin contractility consisting of the force-balance and myosin transport equations. The models account for isotropic contraction proportional to myosin density, viscous stresses in the actin network, and constant-strength viscous-like adhesion. The contraction generates a spatially graded centripetal actin flow, which in turn reinforces the contraction via myosin redistribution and causes retraction of the lamellipodial boundary. Actin protrusion at the boundary counters the retraction, and the balance of the protrusion and retraction shapes the lamellipodium. The model analysis shows that initiation of motility critically depends on three dimensionless parameter combinations, which represent myosin-dependent contractility, a characteristic viscosity-adhesion length, and a rate of actin protrusion. When the contractility is sufficiently strong, cells break symmetry and move steadily along either straight or circular trajectories, and the motile behavior is sensitive to conditions at the cell boundary. Scanning of a model parameter space shows that the contractile mechanism of motility supports robust cell turning in conditions where short viscosity-adhesion lengths and fast protrusion cause an accumulation of myosin in a small region at the cell rear, destabilizing the axial symmetry of a moving cell.

To understand shapes and movements of simple motile cells, we systematically explore minimal models describing a cell as a two-dimensional actin-myosin gel with a free boundary. The models account for actin-myosin contraction balanced by viscous stresses in the actin gel and uniform adhesion. The myosin contraction causes the lamellipodial boundary to retract. Actin protrusion at the boundary counters the retraction, and the balance of protrusion and retraction shapes the cell. The models reproduce a variety of motile shapes observed experimentally. The analysis shows that the mechanical state of a cell depends on a small number of parameters. We find that when the contractility is sufficiently strong, cells break symmetry and move steadily along either straight or circular trajectory. Scanning model parameters shows that the contractile mechanism of motility supports robust cell turning behavior in conditions where deformable actin gel and fast protrusion destabilize the axial symmetry of a moving cell.

Complex self-propelled rings: a minimal model for cell motility

Collective behavior of active matter is observed for self-propelled particles, such as vibrated disks and active Brownian particles, as well as for cytoskeletal filaments in motile cells. Here, a system of quasi two-dimensional penetrable self-propelled rods inside rigid rings is used to construct a complex self-propelled particle. The rods interact sterically with each other and with a stationary or mobile ring via a separation-shifted Lennard-Jones potential. They either have a sliding attachment to the inside of the ring at one of their ends, or can move freely within the ring confinement. We study the inner structure and dynamics of the mobile self-propelled rings. We find that these complex particles cannot only be characterized as active Brownian particles, but can also exhibit cell-like motility: random walks, persistent motion, circling, and run-and-circle motion.

Antagonistic Behaviors of NMY-1 and NMY-2 Maintain Ring Channels in the C. elegans Gonad

Contractile rings play critical roles in a number of biological processes, including oogenesis, wound healing, and cytokinesis. In many cases, the activity of motor proteins such as nonmuscle myosins is required for appropriate constriction of these contractile rings. In the gonad of the nematode worm Caenorhabditis elegans, ring channels are a specialized form of contractile ring that are maintained at a constant diameter before oogenesis. We propose a model of ring channel maintenance that explicitly incorporates force generation by motor proteins that can act normally or tangentially to the ring channel opening. We find that both modes of force generation are needed to maintain the ring channels. We demonstrate experimentally that the type II myosins NMY-1 and NMY-2 antagonize each other in the ring channels by producing force in perpendicular directions: the experimental depletion of NMY-1/theoretical decrease in orthogonal force allows premature ring constriction and cellularization, whereas the experimental depletion of NMY-2/theoretical decrease in tangential force opens the ring channels and prevents cellularization. Together, our experimental and theoretical results show that both forces, mediated by NMY-1 and NMY-2, are crucial for maintaining the appropriate ring channel diameter and dynamics throughout the gonad.

Bistable front dynamics in a contractile medium: Travelling wave fronts and cortical advection define stable zones of RhoA signaling at epithelial adherens junctions

Mechanical coherence of cell layers is essential for epithelia to function as tissue barriers and to control active tissue dynamics during morphogenesis. RhoA signaling at adherens junctions plays a key role in this process by coupling cadherin-based cell-cell adhesion together with actomyosin contractility. Here we propose and analyze a mathematical model representing core interactions involved in the spatial localization of junctional RhoA signaling. We demonstrate how the interplay between biochemical signaling through positive feedback, combined with diffusion on the cell membrane and mechanical forces generated in the cortex, can determine the spatial distribution of RhoA signaling at cell-cell junctions. This dynamical mechanism relies on the balance between a propagating bistable signal that is opposed by an advective flow generated by an actomyosin stress gradient. Experimental observations on the behavior of the system when contractility is inhibited are in qualitative agreement with the predictions of the model.

Physical models of collective cell motility: from cell to tissue

In this article, we review physics-based models of collective cell motility. We discuss a range of techniques at different scales, ranging from models that represent cells as simple self-propelled particles to phase field models that can represent a cell's shape and dynamics in great detail. We also extensively review the ways in which cells within a tissue choose their direction, the statistics of cell motion, and some simple examples of how cell-cell signaling can interact with collective cell motility. This review also covers in more detail selected recent works on collective cell motion of small numbers of cells on micropatterns, in wound healing, and the chemotaxis of clusters of cells.

Crawling and turning in a minimal reaction-diffusion cell motility model: Coupling cell shape and biochemistry

We study a minimal model of a crawling eukaryotic cell with a chemical polarity controlled by a reaction-diffusion mechanism describing Rho GTPase dynamics. The size, shape, and speed of the cell emerge from the combination of the chemical polarity, which controls the locations where actin polymerization occurs, and the physical properties of the cell, including its membrane tension. We find in our model both highly persistent trajectories, in which the cell crawls in a straight line, and turning trajectories, where the cell transitions from crawling in a line to crawling in a circle. We discuss the controlling variables for this turning instability and argue that turning arises from a coupling between the reaction-diffusion mechanism and the shape of the cell. This emphasizes the surprising features that can arise from simple links between cell mechanics and biochemistry. Our results suggest that similar instabilities may be present in a broad class of biochemical descriptions of cell polarity.

Computational model for amoeboid motion: Coupling membrane and cytosol dynamics

A distinguishing feature of amoeboid motion is that the migrating cell undergoes large deformations, caused by the emergence and retraction of actin-rich protrusions, called pseudopods. Here, we propose a cell motility model that represents pseudopod dynamics, as well as its interaction with membrane signaling molecules. The model accounts for internal and external forces, such as protrusion, contraction, adhesion, surface tension, or those arising from cell-obstacle contacts. By coupling the membrane and cytosol interactions we are able to reproduce a realistic picture of amoeboid motion. The model results are in quantitative agreement with experiments and show how cells may take advantage of the geometry of their microenvironment to migrate more efficiently.

Theoretical modeling of actin-retrograde-flow passing clusters of confined T cell receptors

Through the activation process of T cells, actin filaments move from the cell periphery toward the cell center. The moving filaments engage with T cell receptors and thus contribute to transportation of the signaling molecules. To study the connection between the moving actin filaments and T cell receptors, an experiment available in the literature has measured filaments flow velocity passing over a region of confined clusters of receptors. It shows that flow velocity decreases in the proximity of the receptors, and then regains its normal value after traversing the region, suggesting a dissipative friction-like connection. In this work, we develop a minimal theoretical model to re-examine this experiment. The model brings the insight that, in contrast to the first impression that the experiment gives, the direct necessity of having a minimum in the velocity profile is not the locally high friction region, but a combined driving force of push from upstream and pull from within and downstream of the system. The predicted driving force integrates our current understanding of the spatially dependent role of the myosin motor proteins and the actin-polymerization-machinery, which make the pulling and pushing forces, respectively.

Prediction of traction forces of motile cells

When crawling on a flat substrate, living cells exert forces on it via adhesive contacts, enabling them to build up tension within their cytoskeleton and to change shape. The measurement of these forces has been made possible by traction force microscopy (TFM), a technique which has allowed us to obtain time-resolved traction force maps during cell migration. This cell 'footprint' is, however, not sufficient to understand the details of the mechanics of migration, that is how cytoskeletal elements (respectively, adhesion complexes) are put under tension and reinforce or deform (respectively, mature and/or unbind) as a result. In a recent paper, we have validated a rheological model of actomyosin linking tension, deformation and myosin activity. Here, we complement this model with tentative models of the mechanics of adhesion and explore how closely these models can predict the traction forces that we recover from experimental measurements during cell migration. The resulting mathematical problem is a PDE set on the experimentally observed domain, which we solve using a finite-element approach. The four parameters of the model can then be adjusted by comparison with experimental results on a single frame of an experiment, and then used to test the predictive power of the model for following frames and other experiments. It is found that the basic pattern of traction forces is robustly predicted by the model and fixed parameters as a function of current geometry only.

Actomyosin contraction, aggregation and traveling waves in a treadmilling actin array

We use perturbation theory to derive a continuum model for the dynamic actomyosin bundle/ring in the regime of very strong crosslinking. Actin treadmilling is essential for contraction. Linear stability analysis and numerical solutions of the model equations reveal that when the actin treadmilling is very slow, actin and myosin aggregate into equidistantly spaced peaks. When treadmilling is significant, actin filament of one polarity are distributed evenly, while filaments of the opposite polarity develop a shock wave moving with the treadmilling velocity. Myosin aggregates into a sharp peak surfing the crest of the actin wave. Any actomyosin aggregation diminishes contractile stress. The easiest way to maintain higher contraction is to upregulate the actomyosin turnover which destabilizes nontrivial patterns and stabilizes the homogeneous actomyosin distributions. We discuss the model's implications for the experiment.

Growth, collapse, and stalling in a mechanical model for neurite motility

Neurites, the long cellular protrusions that form the routes of the neuronal network, are capable of actively extending during early morphogenesis or regenerating after trauma. To perform this task, they rely on their cytoskeleton for mechanical support. In this paper, we present a three-component active gel model that describes neurites in the three robust mechanical states observed experimentally: collapsed, static, and motile. These states arise from an interplay between the physical forces driven by the growth of the microtubule-rich inner core of the neurite and the acto-myosin contractility of its surrounding cortical membrane. In particular, static states appear as a mechanical balance between traction and compression of these two parallel structures. The model predicts how the response of a neurite to a towing force depends on the force magnitude and recovers the response of neurites to several drug treatments that modulate the cytoskeleton active and passive properties.

Assembly and positioning of actomyosin rings by contractility and planar cell polarity

The actomyosin cytoskeleton is a primary force-generating mechanism in morphogenesis, thus a robust spatial control of cytoskeletal positioning is essential. In this report, we demonstrate that actomyosin contractility and planar cell polarity (PCP) interact in post-mitotic Ciona notochord cells to self-assemble and reposition actomyosin rings, which play an essential role for cell elongation. Intriguingly, rings always form at the cells' anterior edge before migrating towards the center as contractility increases, reflecting a novel dynamical property of the cortex. Our drug and genetic manipulations uncover a tug-of-war between contractility, which localizes cortical flows toward the equator and PCP, which tries to reposition them. We develop a simple model of the physical forces underlying this tug-of-war, which quantitatively reproduces our results. We thus propose a quantitative framework for dissecting the relative contribution of contractility and PCP to the self-assembly and repositioning of cytoskeletal structures, which should be applicable to other morphogenetic events.

Celebrating Soft Matter's 10th Anniversary: Cell division: a source of active stress in cellular monolayers

We introduce the notion of cell division-induced activity and show that the cell division generates extensile forces and drives dynamical patterns in cell assemblies. Extending the hydrodynamic models of lyotropic active nematics we describe turbulent-like velocity fields that are generated by the cell division in a confluent monolayer of cells. We show that the experimentally measured flow field of dividing Madin-Darby Canine Kidney (MDCK) cells is reproduced by our modeling approach. Division-induced activity acts together with intrinsic activity of the cells in extensile and contractile cell assemblies to change the flow and director patterns and the density of topological defects. Finally we model the evolution of the boundary of a cellular colony and compare the fingering instabilities induced by cell division to experimental observations on the expansion of MDCK cell cultures.

The mechanisms of spatial and temporal patterning of cell-edge dynamics

Adherent cells migrate and change their shape by means of protrusion and retraction at their edges. When and where these activities occur defines the shape of the cell and the way it moves. Despite a great deal of knowledge about the structural organization, components, and biochemical reactions involved in protrusion and retraction, the origins of their spatial and temporal patterns are still poorly understood. Chemical signaling circuitry is believed to be an important source of patterning, but recent studies highlighted mechanisms based on physical forces, motion, and mechanical feedback.