Wrapped Up: The Motility of Polarly Flagellated Bacteria.

A huge number of bacterial species are motile by flagella, which allow them to actively move toward favorable environments and away from hazardous areas and to conquer new habitats. The general perception of flagellum-mediated movement and chemotaxis is dominated by the Escherichia coli paradigm, with its peritrichous flagellation and its famous run-and-tumble navigation pattern, which has shaped the view on how bacteria swim and navigate in chemical gradients. However, a significant amount-more likely the majority-of bacterial species exhibit a (bi)polar flagellar localization pattern instead of lateral flagella. Accordingly, these species have evolved very different mechanisms for navigation and chemotaxis. Here, we review the earlier and recent findings on the various modes of motility mediated by polar flagella.

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

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

Cargo Size Limits and Forces of Cell-Driven Microtransport.

The integration of motile cells into biohybrid microrobots offers unique properties such as sensitive responses to external stimuli, resilience, and intrinsic energy supply. Here, biohybrid cell-cargo systems that are driven by amoeboid Dictyostelium discoideum cells are studied and how the cargo speed and the resulting viscous drag force scales with increasing radius of the spherical cargo particle are explored. Using a simplified geometrical model of the cell-cargo interaction, the findings toward larger cargo sizes, which are not accessible with the experimental setup, are extrapolated and a maximal cargo size is predicted beyond which active cell-driven movements will stall. The active forces exerted by the cells to move a cargo show mechanoresponsive adaptation and increase dramatically when challenged by an external pulling force, a mechanism that may become relevant when navigating cargo through complex heterogeneous environments.

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.

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.

Versatile Patterns in the Actin Cortex of Motile Cells: Self-Organized Pulses Can Coexist with Macropinocytic Ring-Shaped Waves.

Self-organized patterns in the actin cytoskeleton are essential for eukaryotic cellular life. They are the building blocks of many functional structures that often operate simultaneously to facilitate, for example, nutrient uptake and movement of cells. However, identifying how qualitatively distinct actin patterns can coexist remains a challenge. Using bifurcation theory of a mass conserved activator-inhibitor system, we uncover a generic mechanism of how different actin waves-traveling waves and excitable pulses-organize and simultaneously emerge. Live-cell imaging experiments indeed reveal that narrow, planar, and fast-moving excitable pulses may coexist with ring-shaped macropinocytic actin waves in the cortex of motile amoeboid cells.

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.

Inferring the Chemotactic Strategy of P. putida and E. coli Using Modified Kramers-Moyal Coefficients.

Many bacteria perform a run-and-tumble random walk to explore their surrounding and to perform chemotaxis. In this article we present a novel method to infer the relevant parameters of bacterial motion from experimental trajectories including the tumbling events. We introduce a stochastic model for the orientation angle, where a shot-noise process initiates tumbles, and analytically calculate conditional moments, reminiscent of Kramers-Moyal coefficients. Matching them with the moments calculated from experimental trajectories of the bacteria E. coli and Pseudomonas putida, we are able to infer their respective tumble rates, the rotational diffusion constants, and the distributions of tumble angles in good agreement with results from conventional tumble recognizers. We also define a novel tumble recognizer, which explicitly quantifies the error in recognizing tumbles. In the presence of a chemical gradient we condition the moments on the bacterial direction of motion and thereby explore the chemotaxis strategy. For both bacteria we recover and quantify the classical chemotactic strategy, where the tumble rate is smallest along the chemical gradient. In addition, for E. coli we detect some cells, which bias their mean tumble angle towards smaller values. Our findings are supported by a scaling analysis of appropriate ratios of conditional moments, which are directly calculated from experimental data.

Non-Gaussianity, population heterogeneity, and transient superdiffusion in the spreading dynamics of amoeboid cells.

What is the underlying diffusion process governing the spreading dynamics and search strategies employed by amoeboid cells? Based on the statistical analysis of experimental single-cell tracking data of the two-dimensional motion of the Dictyostelium discoideum amoeboid cells, we quantify their diffusive behaviour based on a number of standard and complementary statistical indicators. We compute the ensemble- and time-averaged mean-squared displacements (MSDs) of the diffusing amoebae cells and observe a pronounced spread of short-time diffusion coefficients and anomalous MSD-scaling exponents for individual cells. The distribution functions of the cell displacements, the long-tailed distribution of instantaneous speeds, and the velocity autocorrelations are also computed. In particular, we observe a systematic superdiffusive short-time behaviour for the ensemble- and time-averaged MSDs of the amoeboid cells. Also, a clear anti-correlation of scaling exponents and generalised diffusivity values for different cells is detected. Most significantly, we demonstrate that the distribution function of the cell displacements has a strongly non-Gaussian shape and-using a rescaled spatio-temporal variable-the cell-displacement data collapse onto a universal master curve. The current analysis of single-cell motions can be implemented for quantifying diffusive behaviours in other living-matter systems, in particular, when effects of active transport, non-Gaussian displacements, and heterogeneity of the population are involved in the dynamics.

Variability and Order in Cytoskeletal Dynamics of Motile Amoeboid Cells.

The chemotactic motion of eukaryotic cells such as leukocytes or metastatic cancer cells relies on membrane protrusions driven by the polymerization and depolymerization of actin. Here we show that the response of the actin system to a receptor stimulus is subject to a threshold value that varies strongly from cell to cell. Above the threshold, we observe pronounced cell-to-cell variability in the response amplitude. The polymerization time, however, is almost constant over the entire range of response amplitudes, while the depolymerization time increases with increasing amplitude. We show that cell-to-cell variability in the response amplitude correlates with the amount of Arp2/3, a protein that enhances actin polymerization. A time-delayed feedback model for the cortical actin concentration is consistent with all our observations and confirms the role of Arp2/3 in the observed cell-to-cell variability. Taken together, our observations highlight robust regulation of the actin response that enables a reliable timing of cell movement.

Signaling in chemotactic amoebae remains spatially confined to stimulated membrane regions.

Recent work has demonstrated that the receptor-mediated signaling system in chemotactic amoeboid cells shows typical properties of an excitable system. Here, we delivered spatially confined stimuli of the chemoattractant cAMP to the membrane of differentiated Dictyostelium discoideum cells to investigate whether localized receptor stimuli can induce the spreading of excitable waves in the G-protein-dependent signal transduction system. By imaging the spatiotemporal dynamics of fluorescent markers for phosphatidylinositol (3,4,5)-trisphosphate (PIP3), PTEN and filamentous actin, we observed that the activity of the signaling pathway remained spatially confined to the stimulated membrane region. Neighboring parts of the membrane were not excited and no receptor-initiated spatial spreading of excitation waves was observed. To generate localized cAMP stimuli, either particles that carried covalently bound cAMP molecules on their surface were brought into contact with the cell or a patch of the cell membrane was aspirated into a glass micropipette to shield this patch against freely diffusing cAMP molecules in the surrounding medium. Additionally, the binding site of the cAMP receptor was probed with different surface-immobilized cAMP molecules, confirming results from earlier ligand-binding studies.

Actin and PIP3 waves in giant cells reveal the inherent length scale of an excited state.

The membrane and actin cortex of a motile cell can autonomously differentiate into two states, one typical of the front, the other of the tail. On the substrate-attached surface of Dictyostelium discoideum cells, dynamic patterns of front-like and tail-like states are generated that are well suited to monitor transitions between these states. To image large-scale pattern dynamics independently of boundary effects, we produced giant cells by electric-pulse-induced cell fusion. In these cells, actin waves are coupled to the front and back of phosphatidylinositol (3,4,5)-trisphosphate (PIP3)-rich bands that have a finite width. These composite waves propagate across the plasma membrane of the giant cells with undiminished velocity. After any disturbance, the bands of PIP3 return to their intrinsic width. Upon collision, the waves locally annihilate each other and change direction; at the cell border they are either extinguished or reflected. Accordingly, expanding areas of progressing PIP3 synthesis become unstable beyond a critical radius, their center switching from a front-like to a tail-like state. Our data suggest that PIP3 patterns in normal-sized cells are segments of the self-organizing patterns that evolve in giant cells.

Noisy Oscillations in the Actin Cytoskeleton of Chemotactic Amoeba.

Biological systems with their complex biochemical networks are known to be intrinsically noisy. Here we investigate the dynamics of actin polymerization of amoeboid cells, which are close to the onset of oscillations. We show that the large phenotypic variability in the polymerization dynamics can be accurately captured by a generic nonlinear oscillator model in the presence of noise. We determine the relative role of the noise with a single dimensionless, experimentally accessible parameter, thus providing a quantitative description of the variability in a population of cells. Our approach, which rests on a generic description of a system close to a Hopf bifurcation and includes the effect of noise, can characterize the dynamics of a large class of noisy systems close to an oscillatory instability.

How imperfect mixing and differential diffusion accelerate the rate of nonlinear reactions in microfluidic channels.

In this paper, we show experimentally that inside a microfluidic device, where the reactants are segregated, the reaction rate of an autocatalytic clock reaction is accelerated in comparison to the case where all the reactants are well mixed. We also find that, when mixing is enhanced inside the microfluidic device by introducing obstacles into the flow, the clock reaction becomes slower in comparison to the device where mixing is less efficient. Based on numerical simulations, we show that this effect can be explained by the interplay of nonlinear reaction kinetics (cubic autocatalysis) and differential diffusion, where the autocatalytic species diffuses slower than the substrate.

Actin cytoskeleton of chemotactic amoebae operates close to the onset of oscillations.

The rapid reorganization of the actin cytoskeleton in response to external stimuli is an essential property of many motile eukaryotic cells. Here, we report evidence that the actin machinery of chemotactic Dictyostelium cells operates close to an oscillatory instability. When averaging the actin response of many cells to a short pulse of the chemoattractant cAMP, we observed a transient accumulation of cortical actin reminiscent of a damped oscillation. At the single-cell level, however, the response dynamics ranged from short, strongly damped responses to slowly decaying, weakly damped oscillations. Furthermore, in a small subpopulation, we observed self-sustained oscillations in the cortical F-actin concentration. To substantiate that an oscillatory mechanism governs the actin dynamics in these cells, we systematically exposed a large number of cells to periodic pulse trains of different frequencies. Our results indicate a resonance peak at a stimulation period of around 20 s. We propose a delayed feedback model that explains our experimental findings based on a time-delay in the regulatory network of the actin system. To test the model, we performed stimulation experiments with cells that express GFP-tagged fusion proteins of Coronin and actin-interacting protein 1, as well as knockout mutants that lack Coronin and actin-interacting protein 1. These actin-binding proteins enhance the disassembly of actin filaments and thus allow us to estimate the delay time in the regulatory feedback loop. Based on this independent estimate, our model predicts an intrinsic period of 20 s, which agrees with the resonance observed in our periodic stimulation experiments.

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.

A bacterial swimmer with two alternating speeds of propagation.

We recorded large data sets of swimming trajectories of the soil bacterium Pseudomonas putida. Like other prokaryotic swimmers, P. putida exhibits a motion pattern dominated by persistent runs that are interrupted by turning events. An in-depth analysis of their swimming trajectories revealed that the majority of the turning events is characterized by an angle of ϕ1 = 180° (reversals). To a lesser extent, turning angles of ϕ2 = 0° are also found. Remarkably, we observed that, upon a reversal, the swimming speed changes by a factor of two on average-a prominent feature of the motion pattern that, to our knowledge, has not been reported before. A theoretical model, based on the experimental values for the average run time and the rotational diffusion, recovers the mean-square displacement of P. putida if the two distinct swimming speeds are taken into account. Compared to a swimmer that moves with a constant intermediate speed, the mean-square displacement is strongly enhanced. We furthermore observed a negative dip in the directional autocorrelation at intermediate times, a feature that is only recovered in an extended model, where the nonexponential shape of the run-time distribution is taken into account.

Cross-diffusion in the two-variable Oregonator model.

We explore the effect of cross-diffusion on pattern formation in the two-variable Oregonator model of the Belousov-Zhabotinsky reaction. For high negative cross-diffusion of the activator (the activator being attracted towards regions of increased inhibitor concentration) we find, depending on the values of the parameters, Turing patterns, standing waves, oscillatory Turing patterns, and quasi-standing waves. For the inhibitor, we find that positive cross-diffusion (the inhibitor being repelled by increasing concentrations of the activator) can induce Turing patterns, jumping waves and spatially modulated bulk oscillations. We qualitatively explain the formation of these patterns. With one model we can explain Turing patterns, standing waves and jumping waves, which previously was done with three different models.

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.

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.