Phase-field approach to chemotactic driving of neutrophil morphodynamics

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.

Generation of a Virtual Cell using a Phase Field Approach to Model Amoeboid Crawling

The process of chemotaxis of living cells is complex. Cells follow gradients of an external signal because the interior of the cells gets polarized. The description of the exterior and the interior of the cell together with its motion for the convenient realization of the computational modeling of the whole process is a complex technical problem. Here, we employ a phase field model to characterize the interior of the cell, permitting the integration of stochastic partial differential equations, responsible for the polarization in the interior of the cell, and simultaneously, the calculation of the shape deformations of the cell, including its locomotion. We detail the mathematical description of the process and the procedure to calculate numerically the phase field with a simple reaction-diffusion equation for a single concentration.

The ECM and tissue architecture are major determinants of early invasion mediated by E-cadherin dysfunction

Germline mutations of E-cadherin cause Hereditary Diffuse Gastric Cancer (HDGC), a highly invasive cancer syndrome characterised by the occurrence of diffuse-type gastric carcinoma and lobular breast cancer. In this disease, E-cadherin-defective cells are detected invading the adjacent stroma since very early stages. Although E-cadherin loss is well established as a triggering event, other determinants of the invasive process persist largely unknown. Herein, we develop an experimental strategy that comprises in vitro extrusion assays using E-cadherin mutants associated to HDGC, as well as mathematical models epitomising epithelial dynamics and its interaction with the extracellular matrix (ECM). In vitro, we verify that E-cadherin dysfunctional cells detach from the epithelial monolayer and extrude basally into the ECM. Through phase-field modelling we demonstrate that, aside from loss of cell-cell adhesion, increased ECM attachment further raises basal extrusion efficiency. Importantly, by combining phase-field and vertex model simulations, we show that the cylindrical structure of gastric glands strongly promotes the cell's invasive ability. Moreover, we validate our findings using a dissipative particle dynamics simulation of epithelial extrusion. Overall, we provide the first evidence that cancer cell invasion is the outcome of defective cell-cell linkages, abnormal interplay with the ECM, and a favourable 3D tissue structure.

Spontaneous transitions between amoeboid and keratocyte-like modes of migration

The motility of adherent eukaryotic cells is driven by the dynamics of the actin cytoskeleton. Despite the common force-generating actin machinery, different cell types often show diverse modes of locomotion that differ in their shape dynamics, speed, and persistence of motion. Recently, experiments in Dictyostelium discoideum have revealed that different motility modes can be induced in this model organism, depending on genetic modifications, developmental conditions, and synthetic changes of intracellular signaling. Here, we report experimental evidence that in a mutated D. discoideum cell line with increased Ras activity, switches between two distinct migratory modes, the amoeboid and fan-shaped type of locomotion, can even spontaneously occur within the same cell. We observed and characterized repeated and reversible switchings between the two modes of locomotion, suggesting that they are distinct behavioral traits that coexist within the same cell. We adapted an established phenomenological motility model that combines a reaction-diffusion system for the intracellular dynamics with a dynamic phase field to account for our experimental findings.

A computational study of amoeboid motility in 3D: the role of extracellular matrix geometry, cell deformability, and cell-matrix adhesion

Amoeboid cells often migrate using pseudopods, which are membrane protrusions that grow, bifurcate, and retract dynamically, resulting in a net cell displacement. Many cells within the human body, such as immune cells, epithelial cells, and even metastatic cancer cells, can migrate using the amoeboid phenotype. Amoeboid motility is a complex and multiscale process, where cell deformation, biochemistry, and cytosolic and extracellular fluid motions are coupled. Furthermore, the extracellular matrix (ECM) provides a confined, complex, and heterogeneous environment for the cells to navigate through. Amoeboid cells can migrate without significantly remodeling the ECM using weak or no adhesion, instead utilizing their deformability and the microstructure of the ECM to gain enough traction. While a large volume of work exists on cell motility on 2D substrates, amoeboid motility is 3D in nature. Despite recent progress in modeling cellular motility in 3D, there is a lack of systematic evaluations of the role of ECM microstructure, cell deformability, and adhesion on 3D motility. To fill this knowledge gap, here we present a multiscale, multiphysics modeling study of amoeboid motility through 3D-idealized ECM. The model is a coupled fluid‒structure and coarse-grain biochemistry interaction model that accounts for large deformation of cells, pseudopod dynamics, cytoplasmic and extracellular fluid motion, stochastic dynamics of cell-ECM adhesion, and microstructural (pore-scale) geometric details of the ECM. The key finding of the study is that cell deformation and matrix porosity strongly influence amoeboid motility, while weak adhesion and microscale structural details of the ECM have secondary but subtle effects.

Adhesion modulates cell morphology and migration within dense fibrous networks

One of the most fundamental abilities required for the sustainability of complex life forms is active cell migration, since it is essential in diverse processes from morphogenesis to leukocyte chemotaxis in immune response. The movement of a cell is the result of intricate mechanisms, that involve the coordination between mechanical forces, biochemical regulatory pathways and environmental cues. In particular, epithelial cancer cells have to employ mechanical strategies in order to migrate through the tissue's basement membrane and infiltrate the bloodstream during the invasion stage of metastasis. In this work we explore how mechanical interactions such as spatial restriction and adhesion affect migration of a self-propelled droplet in dense fibrous media. We have performed a systematic analysis using a phase-field model and we propose a novel approach to simulate cell migration with dissipative particle dynamics modelling. With this purpose we have measured in our simulation the cell's velocity and quantified its morphology as a function of the fibre density and of its adhesiveness to the matrix fibres. Furthermore, we have compared our results to a previousin vitromigration assay of fibrosarcoma cells in fibrous matrices. The results show good agreement between the two methodologies and experiments in the literature, which indicates that these minimalist descriptions are able to capture the main features of the system. Our results indicate that adhesiveness is critical for cell migration, by modulating cell morphology in crowded environments and by enhancing cell velocity. In addition, our analysis suggests that matrix metalloproteinases (MMPs) play an important role as adhesiveness modulators. We propose that new assays should be carried out to address the role of adhesion and the effect of different MMPs in cell migration under confined conditions.

How cortical waves drive fission of motile cells

Cytokinesis-the division of a cell into two daughter cells-is a key step in cell growth and proliferation. It typically occurs in synchrony with the cell cycle to ensure that a complete copy of the genetic information is passed on to the next generation of daughter cells. In animal cells, cytokinesis commonly relies on an actomyosin contractile ring that drives equatorial furrowing and separation into the two daughter cells. However, also contractile ring-independent forms of cell division are known that depend on substrate-mediated traction forces. Here, we report evidence of an as yet unknown type of contractile ring-independent cytokinesis that we termed wave-mediated cytofission. It is driven by self-organized cortical actin waves that travel across the ventral membrane of oversized, multinucleated Dictyostelium discoideum cells. Upon collision with the cell border, waves may initiate the formation of protrusions that elongate and eventually pinch off to form separate daughter cells. They are composed of a stable elongated wave segment that is enclosed by a cell membrane and moves in a highly persistent fashion. We rationalize our observations based on a noisy excitable reaction-diffusion model in combination with a dynamic phase field to account for the cell shape and demonstrate that daughter cells emerging from wave-mediated cytofission exhibit a well-controlled size.

A minimal computational model for three-dimensional cell migration

During migration, eukaryotic cells can continuously change their three-dimensional morphology, resulting in a highly dynamic and complex process. Further complicating this process is the observation that the same cell type can rapidly switch between different modes of migration. Modelling this complexity necessitates models that are able to track deforming membranes and that can capture the intracellular dynamics responsible for changes in migration modes. Here we develop an efficient three-dimensional computational model for cell migration, which couples cell mechanics to a simple intracellular activator-inhibitor signalling system. We compare the computational results to quantitative experiments using the social amoeba Dictyostelium discoideum. The model can reproduce the observed migration modes generated by varying either mechanical or biochemical model parameters and suggests a coupling between the substrate and the biomechanics of the cell.

Plasticity of cell migration resulting from mechanochemical coupling

Eukaryotic cells can migrate using different modes, ranging from amoeboid-like, during which actin filled protrusions come and go, to keratocyte-like, characterized by a stable morphology and persistent motion. How cells can switch between these modes is not well understood but waves of signaling events are thought to play an important role in these transitions. Here we present a simple two-component biochemical reaction-diffusion model based on relaxation oscillators and couple this to a model for the mechanics of cell deformations. Different migration modes, including amoeboid-like and keratocyte-like, naturally emerge through transitions determined by interactions between biochemical traveling waves, cell mechanics and morphology. The model predictions are explicitly verified by systematically reducing the protrusive force of the actin network in experiments using Dictyostelium discoideum cells. Our results indicate the importance of coupling signaling events to cell mechanics and morphology and may be applicable in a wide variety of cell motility systems.

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.

Modeling random crawling, membrane deformation and intracellular polarity of motile amoeboid cells

Amoeboid movement is one of the most widespread forms of cell motility that plays a key role in numerous biological contexts. While many aspects of this process are well investigated, the large cell-to-cell variability in the motile characteristics of an otherwise uniform population remains an open question that was largely ignored by previous models. In this article, we present a mathematical model of amoeboid motility that combines noisy bistable kinetics with a dynamic phase field for the cell shape. To capture cell-to-cell variability, we introduce a single parameter for tuning the balance between polarity formation and intracellular noise. We compare numerical simulations of our model to experiments with the social amoeba Dictyostelium discoideum. Despite the simple structure of our model, we found close agreement with the experimental results for the center-of-mass motion as well as for the evolution of the cell shape and the overall intracellular patterns. We thus conjecture that the building blocks of our model capture essential features of amoeboid motility and may serve as a starting point for more detailed descriptions of cell motion in chemical gradients and confined environments.

A computational model of amoeboid cell motility in the presence of obstacles

Locomotion of amoeboid cells is mediated by finger-like protrusions of the cell body, known as pseudopods, which grow, bifurcate, and retract in a dynamic fashion. Pseudopods are the primary mode of locomotion for many cells within the human body, such as leukocytes, embryonic cells, and metastatic cancer cells. Amoeboid motility is a complex and multiscale process, which involves bio-molecular reactions, cell deformation, and cytoplasmic and extracellular fluid motion. Additionally, cells within the human body are subject to a confined 3D environment known as the extra-cellular matrix (ECM), which resembles a fluid-filled porous medium. In this article, we present a 3D, multiphysics computational approach coupling fluid mechanics, solid mechanics, and a pattern formation model to simulate locomotion of amoeboid cells through a porous matrix composed of a viscous fluid and an array of finite-sized spherical obstacles. The model combines reaction-diffusion of activator/inhibitors, extreme deformation of the cell, pseudopod dynamics, cytoplasmic and extracellular fluid motion, and fully resolved extracellular matrix. A surface finite-element method is used to obtain the cell deformation and activator/inhibitor concentrations, while the fluid motion is solved using a combined finite-volume and spectral method. The immersed-boundary methods are used to couple the cell deformation, obstacles, and fluid. The model is able to recreate squeezing and weaving motion of cells through the matrix. We study the influence of matrix porosity, obstacle size, and cell deformability on the motility behavior. It is found that below certain values of these parameters, cell motion is completely inhibited. Phase diagrams are presented depicting such motility limits. Interesting dynamics seen in the presence of obstacles but absent in unconfined medium, such as freezing or cell arrest, probing, doubling-back, and tug-of-war are predicted. Furthermore, persistent unidirectional motion of cells that is often observed in an unconfined medium is shown to be lost in presence of obstacles, and is attributed to an alteration of the pseudopod dynamics. The same mechanism, however, allows the cell to find a new direction to penetrate further into the matrix without being stuck in one place. The results and analysis presented here show a strong coupling between cell deformability and ECM properties, and provide new fluid mechanical insights on amoeboid motility in confined medium.

Critical behavior of subcellular density organization during neutrophil activation and migration

Physical theories of active matter continue to provide a quantitative understanding of dynamic cellular phenomena, including cell locomotion. Although various investigations of the rheology of cells have identified important viscoelastic and traction force parameters for use in these theoretical approaches, a key variable has remained elusive both in theoretical and experimental approaches: the spatiotemporal behavior of the subcellular density. The evolution of the subcellular density has been qualitatively observed for decades as it provides the source of image contrast in label-free imaging modalities (e.g., differential interference contrast, phase contrast) used to investigate cellular specimens. While these modalities directly visualize cell structure, they do not provide quantitative access to the structures being visualized. We present an established quantitative imaging approach, non-interferometric quantitative phase microscopy, to elucidate the subcellular density dynamics in neutrophils undergoing chemokinesis following uniform bacterial peptide stimulation. Through this approach, we identify a power law dependence of the neutrophil mean density on time with a critical point, suggesting a critical density is required for motility on 2D substrates. Next we elucidate a continuum law relating mean cell density, area, and total mass that is conserved during neutrophil polarization and migration. Together, our approach and quantitative findings will enable investigators to define the physics coupling cytoskeletal dynamics with subcellular density dynamics during cell migration.

The chaser and the chased: a phase-field model of an immune response

In this paper we present a model for an immune response to an invading pathogen. Particularly, we follow the motion of a neutrophil as it migrates to the site of infection guided by chemical cues, a mechanism termed chemotaxis, with the ability to reorient itself as the pathogen changes its position. In the process, the cell undergoes morphological alterations, in addition to the structural changes observed at its leading edge. Also, we derive a condition for a successful immune reaction by relating the speed of the neutrophil to that of the pathogen and to the diffusion coefficient of the chemical attractant.