'Dicty dynamics': Dictyostelium motility as persistent random motion

Data-driven classification of individual cells by their non-Markovian motion

We present a method to differentiate organisms solely by their motion based on the generalized Langevin equation (GLE) and use it to distinguish two different swimming modes of strongly confined unicellular microalgae Chlamydomonas reinhardtii. The GLE is a general model for active or passive motion of organisms and particles that can be derived from a time-dependent general many-body Hamiltonian and in particular includes non-Markovian effects (i.e., the trajectory memory of its past). We extract all GLE parameters from individual cell trajectories and perform an unbiased cluster analysis to group them into different classes. For the specific cell population employed in the experiments, the GLE-based assignment into the two different swimming modes works perfectly, as checked by control experiments. The classification and sorting of single cells and organisms is important in different areas; our method, which is based on motion trajectories, offers wide-ranging applications in biology and medicine.

Learning dynamical models of single and collective cell migration: a review

Single and collective cell migration are fundamental processes critical for physiological phenomena ranging from embryonic development and immune response to wound healing and cancer metastasis. To understand cell migration from a physical perspective, a broad variety of models for the underlying physical mechanisms that govern cell motility have been developed. A key challenge in the development of such models is how to connect them to experimental observations, which often exhibit complex stochastic behaviours. In this review, we discuss recent advances in data-driven theoretical approaches that directly connect with experimental data to infer dynamical models of stochastic cell migration. Leveraging advances in nanofabrication, image analysis, and tracking technology, experimental studies now provide unprecedented large datasets on cellular dynamics. In parallel, theoretical efforts have been directed towards integrating such datasets into physical models from the single cell to the tissue scale with the aim of conceptualising the emergent behaviour of cells. We first review how this inference problem has been addressed in both freely migrating and confined cells. Next, we discuss why these dynamics typically take the form of underdamped stochastic equations of motion, and how such equations can be inferred from data. We then review applications of data-driven inference and machine learning approaches to heterogeneity in cell behaviour, subcellular degrees of freedom, and to the collective dynamics of multicellular systems. Across these applications, we emphasise how data-driven methods can be integrated with physical active matter models of migrating cells, and help reveal how underlying molecular mechanisms control cell behaviour. Together, these data-driven approaches are a promising avenue for building physical models of cell migration directly from experimental data, and for providing conceptual links between different length-scales of description.

Nonlinear Growth Dynamics of Neuronal Cells Cultured on Directional Surfaces

During the development of the nervous system, neuronal cells extend axons and dendrites that form complex neuronal networks, which are essential for transmitting and processing information. Understanding the physical processes that underlie the formation of neuronal networks is essential for gaining a deeper insight into higher-order brain functions such as sensory processing, learning, and memory. In the process of creating networks, axons travel towards other recipient neurons, directed by a combination of internal and external cues that include genetic instructions, biochemical signals, as well as external mechanical and geometrical stimuli. Although there have been significant recent advances, the basic principles governing axonal growth, collective dynamics, and the development of neuronal networks remain poorly understood. In this paper, we present a detailed analysis of nonlinear dynamics for axonal growth on surfaces with periodic geometrical patterns. We show that axonal growth on these surfaces is described by nonlinear Langevin equations with speed-dependent deterministic terms and gaussian stochastic noise. This theoretical model yields a comprehensive description of axonal growth at both intermediate and long time scales (tens of hours after cell plating), and predicts key dynamical parameters, such as speed and angular correlation functions, axonal mean squared lengths, and diffusion (cell motility) coefficients. We use this model to perform simulations of axonal trajectories on the growth surfaces, in turn demonstrating very good agreement between simulated growth and the experimental results. These results provide important insights into the current understanding of the dynamical behavior of neurons, the self-wiring of the nervous system, as well as for designing innovative biomimetic neural network models.

Size-dependent self-avoidance enables superdiffusive migration in macroscopic unicellulars

Many cells face search problems, such as finding food, mates, or shelter, where their success depends on their search strategy. In contrast to other unicellular organisms, the slime mold Physarum polycephalum forms a giant network-shaped plasmodium while foraging for food. What is the advantage of the giant cell on the verge of multicellularity? We experimentally study and quantify the migration behavior of P. polycephalum plasmodia on the time scale of days in the absence and presence of food. We develop a model which successfully describes its migration in terms of ten data-derived parameters. Using the mechanistic insights provided by our data-driven model, we find that regardless of the absence or presence of food, P. polycephalum achieves superdiffusive migration by performing a self-avoiding run-and-tumble movement. In the presence of food, the run duration statistics change, only controlling the short-term migration dynamics. However, varying organism size, we find that the long-term superdiffusion arises from self-avoidance determined by cell size, highlighting the potential evolutionary advantage that this macroscopically large cell may have.

An entropy-based approach for assessing the directional persistence of cell migration

Cell migration, which is primarily characterized by directional persistence, is essential for the development of normal tissues and organs, as well as for numerous pathological processes. However, there is a lack of simple and efficient tools to analyze the systematic properties of persistence based on cellular trajectory data. Here, we present a novel approach, the entropy of angular distribution , which combines cellular turning dynamics and Shannon entropy to explore the statistical and time-varying properties of persistence that strongly correlate with cellular migration modes. Our results reveal the changes in the persistence of multiple cell lines that are tightly regulated by both intra- and extracellular cues, including Arpin protein, collagen gel/substrate, and physical constraints. Significantly, some previously unreported distinctive details of persistence have also been captured, helping to elucidate how directional persistence is distributed and evolves in different cell populations. The analysis suggests that the entropy of angular distribution-based approach provides a powerful metric for evaluating directional persistence and enables us to better understand the relationships between cellular behaviors and multiscale cues, which also provides some insights into the migration dynamics of cell populations, such as collective cell invasion.

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.

Random walk and cell morphology dynamics in Naegleria gruberi

Amoeboid cell movement and migration are wide-spread across various cell types and species. Microscopy-based analysis of the model systems Dictyostelium and neutrophils over the years have uncovered generality in their overall cell movement pattern. Under no directional cues, the centroid movement can be quantitatively characterized by their persistence to move in a straight line and the frequency of re-orientation. Mathematically, the cells essentially behave as a persistent random walker with memory of two characteristic time-scale. Such quantitative characterization is important from a cellular-level ethology point of view as it has direct connotation to their exploratory and foraging strategies. Interestingly, outside the amoebozoa and metazoa, there are largely uncharacterized species in the excavate taxon Heterolobosea including amoeboflagellate Naegleria. While classical works have shown that these cells indeed show typical amoeboid locomotion on an attached surface, their quantitative features are so far unexplored. Here, we analyzed the cell movement of Naegleria gruberi by employing long-time phase contrast imaging that automatically tracks individual cells. We show that the cells move as a persistent random walker with two time-scales that are close to those known in Dictyostelium and neutrophils. Similarities were also found in the shape dynamics which are characterized by the appearance, splitting and annihilation of the curvature waves along the cell edge. Our analysis based on the Fourier descriptor and a neural network classifier point to importance of morphology features unique to Naegleria including complex protrusions and the transient bipolar dumbbell morphologies.

Fast detection of slender bodies in high density microscopy data

Computer-aided analysis of biological microscopy data has seen a massive improvement with the utilization of general-purpose deep learning techniques. Yet, in microscopy studies of multi-organism systems, the problem of collision and overlap remains challenging. This is particularly true for systems composed of slender bodies such as swimming nematodes, swimming spermatozoa, or the beating of eukaryotic or prokaryotic flagella. Here, we develop a end-to-end deep learning approach to extract precise shape trajectories of generally motile and overlapping slender bodies. Our method works in low resolution settings where feature keypoints are hard to define and detect. Detection is fast and we demonstrate the ability to track thousands of overlapping organisms simultaneously. While our approach is agnostic to area of application, we present it in the setting of and exemplify its usability on dense experiments of swimming Caenorhabditis elegans. The model training is achieved purely on synthetic data, utilizing a physics-based model for nematode motility, and we demonstrate the model's ability to generalize from simulations to experimental videos.

Biased Random Walk Model of Neuronal Dynamics on Substrates with Periodic Geometrical Patterns

Neuronal networks are complex systems of interconnected neurons responsible for transmitting and processing information throughout the nervous system. The building blocks of neuronal networks consist of individual neurons, specialized cells that receive, process, and transmit electrical and chemical signals throughout the body. The formation of neuronal networks in the developing nervous system is a process of fundamental importance for understanding brain activity, including perception, memory, and cognition. To form networks, neuronal cells extend long processes called axons, which navigate toward other target neurons guided by both intrinsic and extrinsic factors, including genetic programming, chemical signaling, intercellular interactions, and mechanical and geometrical cues. Despite important recent advances, the basic mechanisms underlying collective neuron behavior and the formation of functional neuronal networks are not entirely understood. In this paper, we present a combined experimental and theoretical analysis of neuronal growth on surfaces with micropatterned periodic geometrical features. We demonstrate that the extension of axons on these surfaces is described by a biased random walk model, in which the surface geometry imparts a constant drift term to the axon, and the stochastic cues produce a random walk around the average growth direction. We show that the model predicts key parameters that describe axonal dynamics: diffusion (cell motility) coefficient, average growth velocity, and axonal mean squared length, and we compare these parameters with the results of experimental measurements. Our findings indicate that neuronal growth is governed by a contact-guidance mechanism, in which the axons respond to external geometrical cues by aligning their motion along the surface micropatterns. These results have a significant impact on developing novel neural network models, as well as biomimetic substrates, to stimulate nerve regeneration and repair after injury.

Heterogeneous individual motility biases group composition in a model of aggregating cells

Aggregative life cycles are characterized by alternating phases of unicellular growth and multicellular development. Their multiple, independent evolutionary emergence suggests that they may have coopted pervasive properties of single-celled ancestors. Primitive multicellular aggregates, where coordination mechanisms were less efficient than in extant aggregative microbes, must have faced high levels of conflict between different co-aggregating populations. Such conflicts within a multicellular body manifest in the differential reproductive output of cells of different types. Here, we study how heterogeneity in cell motility affects the aggregation process and creates a mismatch between the composition of the population and that of self-organized groups of active adhesive particles. We model cells as self-propelled particles and describe aggregation in a plane starting from a dispersed configuration. Inspired by the life cycle of aggregative model organisms such as Dictyostelium discoideum or Myxococcus xanthus, whose cells interact for a fixed duration before the onset of chimeric multicellular development, we study finite-time configurations for identical particles and in binary mixes. We show that co-aggregation results in three different types of frequency-dependent biases, one of which is associated to evolutionarily stable coexistence of particles with different motility. We propose a heuristic explanation of such observations, based on the competition between delayed aggregation of slower particles and detachment of faster particles. Unexpectedly, despite the complexity and non-linearity of the system, biases can be largely predicted from the behavior of the two corresponding homogenous populations. This model points to differential motility as a possibly important factor in driving the evolutionary emergence of facultatively multicellular life-cycles.

Feedback-controlled dynamics of neuronal cells on directional surfaces

The formation of neuronal networks is a complex phenomenon of fundamental importance for understanding the development of the nervous system. The basic process underlying the network formation is axonal growth, a process involving the extension of axons from the cell body and axonal navigation toward target neurons. Axonal growth is guided by the interactions between the tip of the axon (growth cone) and its extracellular environmental cues, which include intercellular interactions, the biochemical landscape around the neuron, and the mechanical and geometrical features of the growth substrate. Here, we present a comprehensive experimental and theoretical analysis of axonal growth for neurons cultured on micropatterned polydimethylsiloxane (PDMS) surfaces. We demonstrate that closed-loop feedback is an essential component of axonal dynamics on these surfaces: the growth cone continuously measures environmental cues and adjusts its motion in response to external geometrical features. We show that this model captures all the characteristics of axonal dynamics on PDMS surfaces for both untreated and chemically modified neurons. We combine experimental data with theoretical analysis to measure key parameters that describe axonal dynamics: diffusion (cell motility) coefficients, speed and angular distributions, and cell-substrate interactions. The experiments performed on neurons treated with Taxol (inhibitor of microtubule dynamics) and Y-27632 (disruptor of actin filaments) indicate that the internal dynamics of microtubules and actin filaments plays a critical role for the proper function of the feedback mechanism. Our results demonstrate that axons follow geometrical patterns through a contact-guidance mechanism, in which high-curvature geometrical features impart high traction forces to the growth cone. These results have important implications for our fundamental understanding of axonal growth as well as for bioengineering novel substrate to guide neuronal growth and promote nerve repair.

Cluster size distribution of cells disseminating from a primary tumor

The first stage of the metastatic cascade often involves motile cells emerging from a primary tumor either as single cells or as clusters. These cells enter the circulation, transit to other parts of the body and finally are responsible for growth of secondary tumors in distant organs. The mode of dissemination is believed to depend on the EMT nature (epithelial, hybrid or mesenchymal) of the cells. Here, we calculate the cluster size distribution of these migrating cells, using a mechanistic computational model, in presence of different degree of EMT-ness of the cells; EMT is treated as given rise to changes in their active motile forces (μ) and cell-medium surface tension (Γ). We find that, for (μ > μ_min, Γ > 1), when the cells are hybrid in nature, the mean cluster size, N¯∼Γ2.0/μ2.8, where μ_min increases with increase in Γ. For Γ ≤ 0, N¯=1, the cells behave as completely mesenchymal. In presence of spectrum of hybrid states with different degree of EMT-ness (motility) in primary tumor, the cells which are relatively more mesenchymal (higher μ) in nature, form larger clusters, whereas the smaller clusters are relatively more epithelial (lower μ). Moreover, the heterogeneity in μ is comparatively higher for smaller clusters with respect to that for larger clusters. We also observe that more extended cell shapes promote the formation of smaller clusters. Overall, this study establishes a framework which connects the nature and size of migrating clusters disseminating from a primary tumor with the phenotypic composition of the tumor, and can lead to the better understanding of metastasis.

In the process of metastasis, tumor cells disseminate from the primary tumor either as single cells or multicellular clusters. These clusters are potential contributor to the initiation of secondary tumor in distant organs. Our computational model captures the size distribution of migrating clusters depending on the adhesion and motility of the cells (which determine the degree of their EMT nature). Furthermore, we investigate the effect of heterogeneity of cell types in the primary tumor on the resultant heterogeneity of cell types in clusters of different sizes. We believe that the understanding the formation and nature of these clusters, dangerous actors in the deadly aspect of cancer progression, will be useful for improving prognostic methods and eventually better treatments.

Axonal growth on surfaces with periodic geometrical patterns

The formation of neuron networks is a complex phenomenon of fundamental importance for understanding the development of the nervous system, and for creating novel bioinspired materials for tissue engineering and neuronal repair. The basic process underlying the network formation is axonal growth, a process involving the extension of axons from the cell body towards target neurons. Axonal growth is guided by environmental stimuli that include intercellular interactions, biochemical cues, and the mechanical and geometrical features of the growth substrate. The dynamics of the growing axon and its biomechanical interactions with the growing substrate remains poorly understood. In this paper, we develop a model of axonal motility which incorporates mechanical interactions between the axon and the growth substrate. We combine experimental data with theoretical analysis to measure the parameters that describe axonal growth on surfaces with micropatterned periodic geometrical features: diffusion (cell motility) coefficients, speed and angular distributions, and axon bending rigidities. Experiments performed on neurons treated Taxol (inhibitor of microtubule dynamics) and Blebbistatin (disruptor of actin filaments) show that the dynamics of the cytoskeleton plays a critical role in the axon steering mechanism. Our results demonstrate that axons follow geometrical patterns through a contact-guidance mechanism, in which high-curvature geometrical features impart high traction forces to the growth cone. These results have important implications for our fundamental understanding of axonal growth as well as for bioengineering novel substrates that promote neuronal growth and nerve repair.

Cell migration guided by long-lived spatial memory

Living cells actively migrate in their environment to perform key biological functions-from unicellular organisms looking for food to single cells such as fibroblasts, leukocytes or cancer cells that can shape, patrol or invade tissues. Cell migration results from complex intracellular processes that enable cell self-propulsion, and has been shown to also integrate various chemical or physical extracellular signals. While it is established that cells can modify their environment by depositing biochemical signals or mechanically remodelling the extracellular matrix, the impact of such self-induced environmental perturbations on cell trajectories at various scales remains unexplored. Here, we show that cells can retrieve their path: by confining motile cells on 1D and 2D micropatterned surfaces, we demonstrate that they leave long-lived physicochemical footprints along their way, which determine their future path. On this basis, we argue that cell trajectories belong to the general class of self-interacting random walks, and show that self-interactions can rule large scale exploration by inducing long-lived ageing, subdiffusion and anomalous first-passage statistics. Altogether, our joint experimental and theoretical approach points to a generic coupling between motile cells and their environment, which endows cells with a spatial memory of their path and can dramatically change their space exploration.

Neuronal Growth and Formation of Neuron Networks on Directional Surfaces

The formation of neuron networks is a process of fundamental importance for understanding the development of the nervous system and for creating biomimetic devices for tissue engineering and neural repair. The basic process that controls the network formation is the growth of an axon from the cell body and its extension towards target neurons. Axonal growth is directed by environmental stimuli that include intercellular interactions, biochemical cues, and the mechanical and geometrical properties of the growth substrate. Despite significant recent progress, the steering of the growing axon remains poorly understood. In this paper, we develop a model of axonal motility, which incorporates substrate-geometry sensing. We combine experimental data with theoretical analysis to measure the parameters that describe axonal growth on micropatterned surfaces: diffusion (cell motility) coefficients, speed and angular distributions, and cell-substrate interactions. Experiments performed on neurons treated with inhibitors for microtubules (Taxol) and actin filaments (Y-27632) indicate that cytoskeletal dynamics play a critical role in the steering mechanism. Our results demonstrate that axons follow geometrical patterns through a contact-guidance mechanism, in which geometrical patterns impart high traction forces to the growth cone. These results have important implications for bioengineering novel substrates to guide neuronal growth and promote nerve repair.

Deriving time-varying cellular motility parameters via wavelet analysis

Cell migration, which is regulated by intracellular signaling pathways (ICSP) and extracellular matrix (ECM), plays an indispensable role in many physiological and pathological process such as normal tissue development and cancer metastasis. However, there is a lack of rigorous and quantitative tools for analyzing the time-varying characteristics of cell migration in heterogeneous microenvironment, resulted from, e.g. the time-dependent local stiffness due to microstructural remodeling by migrating cells. Here, we develop a wavelet-analysis approach to derive the time-dependent motility parameters from cell migration trajectories, based on the time-varying persistent random walk model. In particular, the wavelet denoising and wavelet transform are employed to analyze migration velocities and obtain the wavelet power spectrum. Subsequently, the time-dependent motility parameters are derived via Lorentzian power spectrum. Our results based on synthetic data indicate the superiority of the method for estimating the intrinsic transient motility parameters, robust against a variety of stochastic noises. We also carry out a systematic parameter study and elaborate the effects of parameter selection on the performance of the method. Moreover, we demonstrate the utility of our approach via analyzing experimental data ofin vitrocell migration in distinct microenvironments, including the migration of MDA-MB-231 cells in confined micro-channel arrays and correlated migration of MCF-10A cells due to ECM-mediated mechanical coupling. Our analysis shows that our approach can be as a powerful tool to accurately derive the time-dependent motility parameters, and further analyze the time-dependent characteristics of cell migration regulated by complex microenvironment.

Shannon entropy for time-varying persistence of cell migration

Cell migration, which can be significantly affected by intracellular signaling pathways and extracellular matrix, plays a crucial role in many physiological and pathological processes. Cell migration is typically modeled as a persistent random walk, which depends on two critical motility parameters, i.e., migration speed and persistence time. It is generally very challenging to efficiently and accurately quantify the migration dynamics from noisy experimental data. Here, we introduce the normalized Shannon entropy (SE) based on the FPS of cellular velocity autocovariance function to quantify migration dynamics. The SE introduced here possesses a similar physical interpretation as the Gibbs entropy for thermal systems in that SE naturally reflects the degree of order or randomness of cellular migration, attaining the maximal value of unity for purely diffusive migration (i.e., SE = 1 for the most "random" dynamics) and the minimal value of 0 for purely ballistic dynamics (i.e., SE = 0 for the most "ordered" dynamics). We also find that SE is strongly correlated with the migration persistence but is less sensitive to the migration speed. Moreover, we introduce the time-varying SE based on the WPS of cellular dynamics and demonstrate its superior utility to characterize the time-dependent persistence of cell migration, which typically results from complex and time-varying intra- or extracellular mechanisms. We employ our approach to analyze experimental data of in vitro cell migration regulated by distinct intracellular and extracellular mechanisms, exhibiting a rich spectrum of dynamic characteristics. Our analysis indicates that the SE and wavelet transform (i.e., SE-based approach) offers a simple and efficient tool to quantify cell migration dynamics in complex microenvironment.

Learning the dynamics of cell-cell interactions in confined cell migration

The migratory dynamics of cells in physiological processes, ranging from wound healing to cancer metastasis, rely on contact-mediated cell-cell interactions. These interactions play a key role in shaping the stochastic trajectories of migrating cells. While data-driven physical formalisms for the stochastic migration dynamics of single cells have been developed, such a framework for the behavioral dynamics of interacting cells still remains elusive. Here, we monitor stochastic cell trajectories in a minimal experimental cell collider: a dumbbell-shaped micropattern on which pairs of cells perform repeated cellular collisions. We observe different characteristic behaviors, including cells reversing, following, and sliding past each other upon collision. Capitalizing on this large experimental dataset of coupled cell trajectories, we infer an interacting stochastic equation of motion that accurately predicts the observed interaction behaviors. Our approach reveals that interacting noncancerous MCF10A cells can be described by repulsion and friction interactions. In contrast, cancerous MDA-MB-231 cells exhibit attraction and antifriction interactions, promoting the predominant relative sliding behavior observed for these cells. Based on these experimentally inferred interactions, we show how this framework may generalize to provide a unifying theoretical description of the diverse cellular interaction behaviors of distinct cell types.

Non-Markovian data-driven modeling of single-cell motility

Trajectories of human breast cancer cells moving on one-dimensional circular tracks are modeled by the non-Markovian version of the Langevin equation that includes an arbitrary memory function. When averaged over cells, the velocity distribution exhibits spurious non-Gaussian behavior, while single cells are characterized by Gaussian velocity distributions. Accordingly, the data are described by a linear memory model which includes different random walk models that were previously used to account for various aspects of cell motility such as migratory persistence, non-Markovian effects, colored noise, and anomalous diffusion. The memory function is extracted from the trajectory data without restrictions or assumptions, thus making our approach truly data driven, and is used for unbiased single-cell comparison. The cell memory displays time-delayed single-exponential negative friction, which clearly distinguishes cell motion from the simple persistent random walk model and suggests a regulatory feedback mechanism that controls cell migration. Based on the extracted memory function we formulate a generalized exactly solvable cell migration model which indicates that negative friction generates cell persistence over long timescales. The nonequilibrium character of cell motion is investigated by mapping the non-Markovian Langevin equation with memory onto a Markovian model that involves a hidden degree of freedom and is equivalent to the underdamped active Ornstein-Uhlenbeck process.

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

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

Inferring the Dynamics of Underdamped Stochastic Systems

Many complex systems, ranging from migrating cells to animal groups, exhibit stochastic dynamics described by the underdamped Langevin equation. Inferring such an equation of motion from experimental data can provide profound insight into the physical laws governing the system. Here, we derive a principled framework to infer the dynamics of underdamped stochastic systems from realistic experimental trajectories, sampled at discrete times and subject to measurement errors. This framework yields an operational method, Underdamped Langevin Inference, which performs well on experimental trajectories of single migrating cells and in complex high-dimensional systems, including flocks with Viscek-like alignment interactions. Our method is robust to experimental measurement errors, and includes a self-consistent estimate of the inference error.

Topotaxis of active Brownian particles

Recent experimental studies have demonstrated that cellular motion can be directed by topographical gradients, such as those resulting from spatial variations in the features of a micropatterned substrate. This phenomenon, known as topotaxis, has been extensively studied for topographical gradients at the subcellular scale, but can also occur in the presence of a spatially varying density of cell-sized features. Such a large-scale topotaxis has recently been observed in highly motile cells that persistently crawl within an array of obstacles with smoothly varying lattice spacing. We introduce a toy model of large-scale topotaxis, based on active Brownian particles. Using numerical simulations and analytical arguments, we demonstrate that topographical gradients introduce a spatial modulation of the particles' persistence, leading to directed motion toward regions of higher persistence. Our results demonstrate that persistent motion alone is sufficient to drive large-scale topotaxis and could serve as a starting point for more detailed studies on self-propelled particles and cells.

Disentangling the behavioural variability of confined cell migration

Cell-to-cell variability is inherent to numerous biological processes, including cell migration. Quantifying and characterizing the variability of migrating cells is challenging, as it requires monitoring many cells for long time windows under identical conditions. Here, we observe the migration of single human breast cancer cells (MDA-MB-231) in confining two-state micropatterns. To describe the stochastic dynamics of this confined migration, we employ a dynamical systems approach. We identify statistics to measure the behavioural variance of the migration, which significantly exceeds that predicted by a population-averaged stochastic model. This additional variance can be explained by the combination of an ‘ageing’ process and population heterogeneity. To quantify population heterogeneity, we decompose the cells into subpopulations of slow and fast cells, revealing the presence of distinct classes of dynamical systems describing the migration, ranging from bistable to limit cycle behaviour. Our findings highlight the breadth of migration behaviours present in cell populations.

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.

Neuron dynamics on directional surfaces

Geometrical features play a very important role in neuronal growth and the formation of functional connections between neuronal cells. Here, we analyze the dynamics of axonal growth for neuronal cells cultured on micro-patterned polydimethylsiloxane surfaces. We utilize fluorescence microscopy to image axons, quantify their dynamics, and demonstrate that periodic geometrical patterns impart strong directional bias to neuronal growth. We quantify axonal alignment and present a general stochastic approach that quantitatively describes the dynamics of the growth cones. Neuronal growth is described by a general phenomenological model, based on a simple automatic controller with a closed-loop feedback system. We demonstrate that axonal alignment on these substrates is determined by the surface geometry, and it is quantified by the deterministic part of the stochastic (Langevin and Fokker-Planck) equations. We also show that the axonal alignment with the surface patterns is greatly suppressed by the neuron treatment with Blebbistatin, a chemical compound that inhibits the activity of myosin II. These results give new insight into the role played by the molecular motors and external geometrical cues in guiding axonal growth, and could lead to novel approaches for bioengineering neuronal regeneration platforms.

The nucleus does not significantly affect the migratory trajectories of amoeba in two-dimensional environments

For a wide range of cells, from bacteria to mammals, locomotion movements are a crucial systemic behavior for cellular life. Despite its importance in a plethora of fundamental physiological processes and human pathologies, how unicellular organisms efficiently regulate their locomotion system is an unresolved question. Here, to understand the dynamic characteristics of the locomotion movements and to quantitatively study the role of the nucleus in the migration of Amoeba proteus we have analyzed the movement trajectories of enucleated and non-enucleated amoebas on flat two-dimensional (2D) surfaces using advanced non-linear physical-mathematical tools and computational methods. Our analysis shows that both non-enucleated and enucleated amoebas display the same kind of dynamic migration structure characterized by highly organized data sequences, super-diffusion, non-trivial long-range positive correlations, persistent dynamics with trend-reinforcing behavior, and move-step fluctuations with scale invariant properties. Our results suggest that the presence of the nucleus does not significantly affect the locomotion of amoeba in 2D environments.

Stochastic and Heterogeneous Cancer Cell Migration: Experiment and Theory

Cell migration, an essential process for normal cell development and cancer metastasis, differs from a simple random walk: the mean-square displacement (〈(Δr)²(t)〉) of cells sometimes shows non-Fickian behavior, and the spatiotemporal correlation function (G(r, t)) of cells is often non-Gaussian. We find that this intriguing cell migration should be attributed to heterogeneity in a cell population, even one with a homogeneous genetic background. There are two limiting types of heterogeneity in a cell population: cellular heterogeneity and temporal heterogeneity. Cellular heterogeneity accounts for the cell-to-cell variation in migration capacity, while temporal heterogeneity arises from the temporal noise in the migration capacity of single cells. We illustrate that both cellular and temporal heterogeneity need to be taken into account simultaneously to elucidate cell migration. We investigate the two-dimensional migration of A549 lung cancer cells using time-lapse microscopy and find that the migration of A549 cells is Fickian but has a non-Gaussian spatiotemporal correlation. We find that when a theoretical model considers both cellular and temporal heterogeneity, the model reproduces all of the anomalous behaviors of cancer cell migration.

Collective regulation of cell motility using an accurate density-sensing system

The capacity of living cells to sense their population density and to migrate accordingly is essential for the regulation of many physiological processes. However, the mechanisms used to achieve such functions are poorly known. Here, based on the analysis of multiple trajectories of vegetative Dictyostelium discoideum cells, we investigate such a system extensively. We show that the cells secrete a high-molecular-weight quorum-sensing factor (QSF) in their medium. This extracellular signal induces, in turn, a reduction of the cell movements, in particular, through the downregulation of a mode of motility with high persistence time. This response appears independent of cAMP and involves a G-protein-dependent pathway. Using a mathematical analysis of the cells' response function, we evidence a negative feedback on the QSF secretion, which unveils a powerful generic mechanism for the cells to detect when they exceed a density threshold. Altogether, our results provide a comprehensive and dynamical view of this system enabling cells in a scattered population to adapt their motion to their neighbours without physical contact.

Anomalous diffusion for neuronal growth on surfaces with controlled geometries

Geometrical cues are known to play a very important role in neuronal growth and the formation of neuronal networks. Here, we present a detailed analysis of axonal growth and dynamics for neuronal cells cultured on patterned polydimethylsiloxane surfaces. We use fluorescence microscopy to image neurons, quantify their dynamics, and demonstrate that the substrate geometrical patterns cause strong directional alignment of axons. We quantify axonal growth and report a general stochastic approach that quantitatively describes the motion of growth cones. The growth cone dynamics is described by Langevin and Fokker-Planck equations with both deterministic and stochastic contributions. We show that the deterministic terms contain both the angular and speed dependence of axonal growth, and that these two contributions can be separated. Growth alignment is determined by surface geometry, and it is quantified by the deterministic part of the Langevin equation. We combine experimental data with theoretical analysis to measure the key parameters of the growth cone motion: speed and angular distributions, correlation functions, diffusion coefficients, characteristics speeds and damping coefficients. We demonstrate that axonal dynamics displays a cross-over from Brownian motion (Ornstein-Uhlenbeck process) at earlier times to anomalous dynamics (superdiffusion) at later times. The superdiffusive regime is characterized by non-Gaussian speed distributions and power law dependence of the axonal mean square length and the velocity correlation functions. These results demonstrate the importance of geometrical cues in guiding axonal growth, and could lead to new methods for bioengineering novel substrates for controlling neuronal growth and regeneration.

Role of geometrical cues in neuronal growth

Geometrical cues play an essential role in neuronal growth. Here, we quantify axonal growth on surfaces with controlled geometries and report a general stochastic approach that quantitatively describes the motion of growth cones. We show that axons display a strong directional alignment on micropatterned surfaces when the periodicity of the patterns matches the dimension of the growth cone. The growth cone dynamics on surfaces with uniform geometry is described by a linear Langevin equation with both deterministic and stochastic contributions. In contrast, axonal growth on surfaces with periodic patterns is characterized by a system of two generalized Langevin equations with both linear and quadratic velocity dependence and stochastic noise. We combine experimental data with theoretical analysis to measure the key parameters of the growth cone motion: angular distributions, correlation functions, diffusion coefficients, characteristics speeds, and damping coefficients. We demonstrate that axonal dynamics displays a crossover from an Ornstein-Uhlenbeck process to a nonlinear stochastic regime when the geometrical periodicity of the pattern approaches the linear dimension of the growth cone. Growth alignment is determined by surface geometry, which is fully quantified by the deterministic part of the Langevin equation. These results provide insight into the role of curvature sensing proteins and their interactions with geometrical cues.

Harnessing Motile Amoeboid Cells as Trucks for Microtransport and -Assembly

Cell-driven microtransport is one of the most prominent applications in the emerging field of biohybrid systems. While bacterial cells have been successfully employed to drive the swimming motion of micrometer-sized cargo particles, the transport capacities of motile adherent cells remain largely unexplored. Here, it is demonstrated that motile amoeboid cells can act as efficient and versatile trucks to transport microcargo. When incubated together with microparticles, cells of the social amoeba Dictyostelium discoideum readily pick up and move the cargo particles. Relying on the unspecific adhesive properties of the amoeba, a wide range of different cargo materials can be used. The cell-driven transport can be directionally guided based on the chemotactic responses of amoeba to chemoattractant gradients. On the one hand, the cargo can be assembled into clusters in a self-organized fashion, relying on the developmentally induced chemotactic aggregation of cells. On the other hand, chemoattractant gradients can be externally imposed to guide the cellular microtrucks to a desired location. Finally, larger cargo particles of different shapes that exceed the size of a single cell by more than an order of magnitude, can also be transported by the collective effort of large numbers of motile 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.

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?

Modelling of Dictyostelium discoideum movement in a linear gradient of chemoattractant

Chemotaxis is a ubiquitous biological phenomenon in which cells detect a spatial gradient of chemoattractant, and then move towards the source. Here we present a position-dependent advection-diffusion model that quantitatively describes the statistical features of the chemotactic motion of the social amoeba Dictyostelium discoideum in a linear gradient of cAMP (cyclic adenosine monophosphate). We fit the model to experimental trajectories that are recorded in a microfluidic setup with stationary cAMP gradients and extract the diffusion and drift coefficients in the gradient direction. Our analysis shows that for the majority of gradients, both coefficients decrease over time and become negative as the cells crawl up the gradient. The extracted model parameters also show that besides the expected drift in the direction of the chemoattractant gradient, we observe a nonlinear dependency of the corresponding variance on time, which can be explained by the model. Furthermore, the results of the model show that the non-linear term in the mean squared displacement of the cell trajectories can dominate the linear term on large time scales.

Discrete Modeling of Amoeboid Locomotion and Chemotaxis in Dictyostelium discoideum by Tracking Pseudopodium Growth Direction

Dictyostelium discoideum amoeba is a well-established model organism for studying the crawling locomotion of eukaryotic cells. These amoebae extend pseudopodium - a temporary actin-based protrusion of their body membrane to probe the medium and crawl through it. Experiments show highly-ordered patterns in the growth direction of these pseudopodia, which results in persistence cell motility. Here, we propose a discrete model for studying and investigating the cell locomotion based on the experimental evidences. According to our model, Dictyostelium selects its pseudopodium growth direction based on a second-order Markov chain process, in the absence of external cues. Consequently, compared to a random walk process, our model indicates stronger growth in the mean-square displacement of cells, which is consistent with empirical findings. In the presence of external chemical stimulants, cells tend to align with the gradient of chemoattractant molecules. To quantify this tendency, we define a coupling coefficient between the pseudopodium extension direction and the gradient of an external stimulant, which depends on the local stimulant concentration and its gradient. Additionally, we generalize the model to weak-coupling regime by utilizing perturbation methods.

Anomalous diffusion and q-Weibull velocity distributions in epithelial cell migration

In multicellular organisms, cell motility is central in all morphogenetic processes, tissue maintenance, wound healing and immune surveillance. Hence, the control of cell motion is a major demand in the creation of artificial tissues and organs. Here, cell migration assays on plastic 2D surfaces involving normal (MDCK) and tumoral (B16F10) epithelial cell lines were performed varying the initial density of plated cells. Through time-lapse microscopy quantities such as speed distributions, velocity autocorrelations and spatial correlations, as well as the scaling of mean-squared displacements were determined. We find that these cells exhibit anomalous diffusion with q-Weibull speed distributions that evolves non-monotonically to a Maxwellian distribution as the initial density of plated cells increases. Although short-ranged spatial velocity correlations mark the formation of small cell clusters, the emergence of collective motion was not observed. Finally, simulational results from a correlated random walk and the Vicsek model of collective dynamics evidence that fluctuations in cell velocity orientations are sufficient to produce q-Weibull speed distributions seen in our migration assays.

Determining the impact of cell mixing on signaling during development

Cell movement and intercellular signaling occur simultaneously to organize morphogenesis during embryonic development. Cell movement can cause relative positional changes between neighboring cells. When intercellular signals are local such cell mixing may affect signaling, changing the flow of information in developing tissues. Little is known about the effect of cell mixing on intercellular signaling in collective cellular behaviors and methods to quantify its impact are lacking. Here we discuss how to determine the impact of cell mixing on cell signaling drawing an example from vertebrate embryogenesis: the segmentation clock, a collective rhythm of interacting genetic oscillators. We argue that comparing cell mixing and signaling timescales is key to determining the influence of mixing. A signaling timescale can be estimated by combining theoretical models with cell signaling perturbation experiments. A mixing timescale can be obtained by analysis of cell trajectories from live imaging. After comparing cell movement analyses in different experimental settings, we highlight challenges in quantifying cell mixing from embryonic timelapse experiments, especially a reference frame problem due to embryonic motions and shape changes. We propose statistical observables characterizing cell mixing that do not depend on the choice of reference frames. Finally, we consider situations in which both cell mixing and signaling involve multiple timescales, precluding a direct comparison between single characteristic timescales. In such situations, physical models based on observables of cell mixing and signaling can simulate the flow of information in tissues and reveal the impact of observed cell mixing on signaling.

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.

Disentangling Random Motion and Flow in a Complex Medium

We describe a technique for deconvolving the stochastic motion of particles from large-scale fluid flow in a dynamic environment such as that found in living cells. The method leverages the separation of timescales to subtract out the persistent component of motion from single-particle trajectories. The mean-squared displacement of the resulting trajectories is rescaled so as to enable robust extraction of the diffusion coefficient and subdiffusive scaling exponent of the stochastic motion. We demonstrate the applicability of the method for characterizing both diffusive and fractional Brownian motion overlaid by flow and analytically calculate the accuracy of the method in different parameter regimes. This technique is employed to analyze the motion of lysosomes in motile neutrophil-like cells, showing that the cytoplasm of these cells behaves as a viscous fluid at the timescales examined.

Diffusion of active chiral particles

The diffusion of chiral active Brownian particles in three-dimensional space is studied analytically, by consideration of the corresponding Fokker-Planck equation for the probability density of finding a particle at position x and moving along the direction v[over ̂] at time t, and numerically, by the use of Langevin dynamics simulations. The analysis is focused on the marginal probability density of finding a particle at a given location and at a given time (independently of its direction of motion), which is found from an infinite hierarchy of differential-recurrence relations for the coefficients that appear in the multipole expansion of the probability distribution, which contains the whole kinematic information. This approach allows the explicit calculation of the time dependence of the mean-squared displacement and the time dependence of the kurtosis of the marginal probability distribution, quantities from which the effective diffusion coefficient and the "shape" of the positions distribution are examined. Oscillations between two characteristic values were found in the time evolution of the kurtosis, namely, between the value that corresponds to a Gaussian and the one that corresponds to a distribution of spherical shell shape. In the case of an ensemble of particles, each one rotating around a uniformly distributed random axis, evidence is found of the so-called effect "anomalous, yet Brownian, diffusion," for which particles follow a non-Gaussian distribution for the positions yet the mean-squared displacement is a linear function of time.

How to connect time-lapse recorded trajectories of motile microorganisms with dynamical models in continuous time

We provide a tool for data-driven modeling of motility, data being time-lapse recorded trajectories. Several mathematical properties of a model to be found can be gleaned from appropriate model-independent experimental statistics, if one understands how such statistics are distorted by the finite sampling frequency of time-lapse recording, by experimental errors on recorded positions, and by conditional averaging. We give exact analytical expressions for these effects in the simplest possible model for persistent random motion, the Ornstein-Uhlenbeck process. Then we describe those aspects of these effects that are valid for any reasonable model for persistent random motion. Our findings are illustrated with experimental data and Monte Carlo simulations.

A Model for Direction Sensing in Dictyostelium discoideum: Ras Activity and Symmetry Breaking Driven by a Gβγ-Mediated, Gα2-Ric8 -- Dependent Signal Transduction Network

Chemotaxis is a dynamic cellular process, comprised of direction sensing, polarization and locomotion, that leads to the directed movement of eukaryotic cells along extracellular gradients. As a primary step in the response of an individual cell to a spatial stimulus, direction sensing has attracted numerous theoretical treatments aimed at explaining experimental observations in a variety of cell types. Here we propose a new model of direction sensing based on experiments using Dictyostelium discoideum (Dicty). The model is built around a reaction-diffusion-translocation system that involves three main component processes: a signal detection step based on G-protein-coupled receptors (GPCR) for cyclic AMP (cAMP), a transduction step based on a heterotrimetic G protein Gα2βγ, and an activation step of a monomeric G-protein Ras. The model can predict the experimentally-observed response of cells treated with latrunculin A, which removes feedback from downstream processes, under a variety of stimulus protocols. We show that [Formula: see text] cycling modulated by Ric8, a nonreceptor guanine exchange factor for [Formula: see text] in Dicty, drives multiple phases of Ras activation and leads to direction sensing and signal amplification in cAMP gradients. The model predicts that both [Formula: see text] and Gβγ are essential for direction sensing, in that membrane-localized [Formula: see text], the activated GTP-bearing form of [Formula: see text], leads to asymmetrical recruitment of RasGEF and Ric8, while globally-diffusing Gβγ mediates their activation. We show that the predicted response at the level of Ras activation encodes sufficient 'memory' to eliminate the 'back-of-the wave' problem, and the effects of diffusion and cell shape on direction sensing are also investigated. In contrast with existing LEGI models of chemotaxis, the results do not require a disparity between the diffusion coefficients of the Ras activator GEF and the Ras inhibitor GAP. Since the signal pathways we study are highly conserved between Dicty and mammalian leukocytes, the model can serve as a generic one for direction sensing.

Patterns of cell thickness oscillations during directional migration of Physarum polycephalum

The functional relationship between the velocity of cell locomotion and intracellular spatial patterns of thickness oscillations in the acellular slime mould Physarum polycephalum was studied. The freely migrating plasmodial cells of 300-800 µm length were tadpole-shaped and displayed thickness oscillations along their longitudinal (body) axis. Two distinct patterns of intracellular thickness oscillations were observed in dependence on the locomotive velocity. The first mode consisted of a single travelling wave that propagated from the rear to the front of the cell. This pattern occurred when the plasmodium migrated slowly. The second mode was a multinodal standing wave that was found during events of fast propagation. Transitions between these two types of cell thickness oscillation patterns took place in narrow propagation velocity intervals. We discuss the possible mechanism leading to these patterns, which are conjectured to modulate both the intracellular pressure and the velocity of free locomotion of the cell.

Characterization of anomalous movements of spherical living cells on a silicon dioxide glassy substrate

The random walk of spherical living cells on a silicon dioxide glassy substrate was studied experimentally and numerically. This random walk trajectory exhibited erratic dancing, which seemingly obeyed anomalous diffusion (i.e., Lévy-like walk) rather than normal diffusion. Moreover, the angular distribution (-π to π) of the cells' trajectory followed a "U-shaped pattern" in comparison to the uniform distribution seen in the movements of negatively charged polystyrene microspheres. These effects could be attributable to the homeostasis-driven structural resilient character of cells and physical interactions derived from temporarily retained nonspecific binding due to weak forces between the cells and substrates. Our results provide new insights into the stochastic behavior of mesoscopic biological particles with respect to structural properties and physical interactions.

Geometry-Driven Polarity in Motile Amoeboid Cells

Motile eukaryotic cells, such as leukocytes, cancer cells, and amoeba, typically move inside the narrow interstitial spacings of tissue or soil. While most of our knowledge of actin-driven eukaryotic motility was obtained from cells that move on planar open surfaces, recent work has demonstrated that confinement can lead to strongly altered motile behavior. Here, we report experimental evidence that motile amoeboid cells undergo a spontaneous symmetry breaking in confined interstitial spaces. Inside narrow channels, the cells switch to a highly persistent, unidirectional mode of motion, moving at a constant speed along the channel. They remain in contact with the two opposing channel side walls and alternate protrusions of their leading edge near each wall. Their actin cytoskeleton exhibits a characteristic arrangement that is dominated by dense, stationary actin foci at the side walls, in conjunction with less dense dynamic regions at the leading edge. Our experimental findings can be explained based on an excitable network model that accounts for the confinement-induced symmetry breaking and correctly recovers the spatio-temporal pattern of protrusions at the leading edge. Since motile cells typically live in the narrow interstitial spacings of tissue or soil, we expect that the geometry-driven polarity we report here plays an important role for movement of cells in their natural environment.

Gaussian memory in kinematic matrix theory for self-propellers

We extend the kinematic matrix ("kinematrix") formalism [Phys. Rev. E 89, 062304 (2014)], which via simple matrix algebra accesses ensemble properties of self-propellers influenced by uncorrelated noise, to treat Gaussian correlated noises. This extension brings into reach many real-world biological and biomimetic self-propellers for which inertia is significant. Applying the formalism, we analyze in detail ensemble behaviors of a 2D self-propeller with velocity fluctuations and orientation evolution driven by an Ornstein-Uhlenbeck process. On the basis of exact results, a variety of dynamical regimes determined by the inertial, speed-fluctuation, orientational diffusion, and emergent disorientation time scales are delineated and discussed.

Quantifying the degree of persistence in random amoeboid motion based on the Hurst exponent of fractional Brownian motion

Amoebae explore their environment in a random way, unless external cues like, e.g., nutrients, bias their motion. Even in the absence of cues, however, experimental cell tracks show some degree of persistence. In this paper, we analyzed individual cell tracks in the framework of a linear mixed effects model, where each track is modeled by a fractional Brownian motion, i.e., a Gaussian process exhibiting a long-term correlation structure superposed on a linear trend. The degree of persistence was quantified by the Hurst exponent of fractional Brownian motion. Our analysis of experimental cell tracks of the amoeba Dictyostelium discoideum showed a persistent movement for the majority of tracks. Employing a sliding window approach, we estimated the variations of the Hurst exponent over time, which allowed us to identify points in time, where the correlation structure was distorted ("outliers"). Coarse graining of track data via down-sampling allowed us to identify the dependence of persistence on the spatial scale. While one would expect the (mode of the) Hurst exponent to be constant on different temporal scales due to the self-similarity property of fractional Brownian motion, we observed a trend towards stronger persistence for the down-sampled cell tracks indicating stronger persistence on larger time scales.

Normal and tumoral melanocytes exhibit q-Gaussian random search patterns

In multicellular organisms, cell motility is central in all morphogenetic processes, tissue maintenance, wound healing and immune surveillance. Hence, failures in its regulation potentiates numerous diseases. Here, cell migration assays on plastic 2D surfaces were performed using normal (Melan A) and tumoral (B16F10) murine melanocytes in random motility conditions. The trajectories of the centroids of the cell perimeters were tracked through time-lapse microscopy. The statistics of these trajectories was analyzed by building velocity and turn angle distributions, as well as velocity autocorrelations and the scaling of mean-squared displacements. We find that these cells exhibit a crossover from a normal to a super-diffusive motion without angular persistence at long time scales. Moreover, these melanocytes move with non-Gaussian velocity distributions. This major finding indicates that amongst those animal cells supposedly migrating through Lévy walks, some of them can instead perform q-Gaussian walks. Furthermore, our results reveal that B16F10 cells infected by mycoplasmas exhibit essentially the same diffusivity than their healthy counterparts. Finally, a q-Gaussian random walk model was proposed to account for these melanocytic migratory traits. Simulations based on this model correctly describe the crossover to super-diffusivity in the cell migration tracks.

Theory of diffusion of active particles that move at constant speed in two dimensions

Starting from a Langevin description of active particles that move with constant speed in infinite two-dimensional space and its corresponding Fokker-Planck equation, we develop a systematic method that allows us to obtain the coarse-grained probability density of finding a particle at a given location and at a given time in arbitrary short-time regimes. By going beyond the diffusive limit, we derive a generalization of the telegrapher equation. Such generalization preserves the hyperbolic structure of the equation and incorporates memory effects in the diffusive term. While no difference is observed for the mean-square displacement computed from the two-dimensional telegrapher equation and from our generalization, the kurtosis results in a sensible parameter that discriminates between both approximations. We carry out a comparative analysis in Fourier space that sheds light on why the standard telegrapher equation is not an appropriate model to describe the propagation of particles with constant speed in dispersive media.