A Hands-on Guide to AmoePy - a Python-Based Software Package to Analyze Cell Migration Data.

Amoeboid cell motility is fundamental for a multitude of biological processes such as embryogenesis, immune responses, wound healing, and cancer metastasis. It is characterized by specific cell shape changes: the extension and retraction of membrane protrusions, known as pseudopodia. A common approach to investigate the mechanisms underlying this type of cell motility is to study phenotypic differences in the locomotion of mutant cell lines. To characterize such differences, methods are required to quantify the contour dynamics of migrating cells. AmoePy is a Python-based software package that provides tools for cell segmentation, contour detection as well as analyzing and simulating contour dynamics. First, a digital representation of the cell contour as a chain of nodes is extracted from each frame of a time-lapse microscopy recording of a moving cell. Then, the dynamics of these nodes-referred to as virtual markers-are tracked as the cell contour evolves over time. From these data, various quantities can be calculated that characterize the contour dynamics, such as the displacement of the virtual markers or the local stretching rate of the marker chain. Their dynamics is typically visualized in space-time plots, the so-called kymographs, where the temporal evolution is displayed for the different locations along the cell contour. Using AmoePy, you can straightforwardly create kymograph plots and videos from stacks of experimental bright-field or fluorescent images of motile cells. A hands-on guide on how to install and use AmoePy is provided in this chapter.

Non-Gaussian Displacements in Active Transport on a Carpet of Motile Cells.

We study the dynamics of micron-sized particles on a layer of motile cells. This cell carpet acts as an active bath that propels passive tracer particles via direct mechanical contact. The resulting nonequilibrium transport shows a crossover from superdiffusive to normal-diffusive dynamics. The particle displacement distribution is distinctly non-Gaussian even at macroscopic timescales exceeding the measurement time. We obtain the distribution of diffusion coefficients from the experimental data and introduce a model for the displacement distribution that matches the experimentally observed non-Gaussian statistics. We argue why similar transport properties are expected for many composite active matter systems.

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.

Black and gray box learning of amplitude equations: Application to phase field systems.

We present a data-driven approach to learning surrogate models for amplitude equations and illustrate its application to interfacial dynamics of phase field systems. In particular, we demonstrate learning effective partial differential equations describing the evolution of phase field interfaces from full phase field data. We illustrate this on a model phase field system, where analytical approximate equations for the dynamics of the phase field interface (a higher-order eikonal equation and its approximation, the Kardar-Parisi-Zhang equation) are known. For this system, we discuss data-driven approaches for the identification of equations that accurately describe the front interface dynamics. When the analytical approximate models mentioned above become inaccurate, as we move beyond the region of validity of the underlying assumptions, the data-driven equations outperform them. In these regimes, going beyond black box identification, we explore different approaches to learning data-driven corrections to the analytically approximate models, leading to effective gray box partial differential equations.

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.

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.