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

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.

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.

Optimal calibration of optical tweezers with arbitrary integration time and sampling frequencies: a general framework [Invited]

Optical tweezers (OT) have become an essential technique in several fields of physics, chemistry, and biology as precise micromanipulation tools and microscopic force transducers. Quantitative measurements require the accurate calibration of the trap stiffness of the optical trap and the diffusion constant of the optically trapped particle. This is typically done by statistical estimators constructed from the position signal of the particle, which is recorded by a digital camera or a quadrant photodiode. The finite integration time and sampling frequency of the detector need to be properly taken into account. Here, we present a general approach based on the joint probability density function of the sampled trajectory that corrects exactly the biases due to the detector's finite integration time and limited sampling frequency, providing theoretical formulas for the most widely employed calibration methods: equipartition, mean squared displacement, autocorrelation, power spectral density, and force reconstruction via maximum-likelihood-estimator analysis (FORMA). Our results, tested with experiments and Monte Carlo simulations, will permit users of OT to confidently estimate the trap stiffness and diffusion constant, extending their use to a broader set of experimental conditions.

Spectral fingerprints of nonequilibrium dynamics: The case of a Brownian gyrator

The same system can exhibit a completely different dynamical behavior when it evolves in equilibrium conditions or when it is driven out-of-equilibrium by, e.g., connecting some of its components to heat baths kept at different temperatures. Here we concentrate on an analytically solvable and experimentally relevant model of such a system-the so-called Brownian gyrator-a two-dimensional nanomachine that performs a systematic, on average, rotation around the origin under nonequilibrium conditions, while no net rotation takes place under equilibrium ones. On this example, we discuss a question whether it is possible to distinguish between two types of a behavior judging not upon the statistical properties of the trajectories of components but rather upon their respective spectral densities. The latter are widely used to characterize diverse dynamical systems and are routinely calculated from the data using standard built-in packages. From such a perspective, we inquire whether the power spectral densities possess some "fingerprint" properties specific to the behavior in nonequilibrium. We show that indeed one can conclusively distinguish between equilibrium and nonequilibrium dynamics by analyzing the cross-correlations between the spectral densities of both components in the short frequency limit, or from the spectral densities of both components evaluated at zero frequency. Our analytical predictions, corroborated by experimental and numerical results, open a new direction for the analysis of a nonequilibrium dynamics.

Renormalization group approach to connect discrete- and continuous-time descriptions of Gaussian processes

Discretization of continuous stochastic processes is needed to numerically simulate them or to infer models from experimental time series. However, depending on the nature of the process, the same discretization scheme may perform very differently for the two tasks, if it is not accurate enough. Exact discretizations, which work equally well at any scale, are characterized by the property of invariance under coarse-graining. Motivated by this observation, we build an explicit renormalization group (RG) approach for Gaussian time series generated by autoregressive models. We show that the RG fixed points correspond to discretizations of linear SDEs, and only come in the form of first order Markov processes or non-Markovian ones. This fact provides an alternative explanation of why standard delay-vector embedding procedures fail in reconstructing partially observed noise-driven systems. We also suggest a possible effective Markovian discretization for the inference of partially observed underdamped equilibrium processes based on the exploitation of the Einstein relation.

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.

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.

CompuCell3D Simulations Reproduce Mesenchymal Cell Migration on Flat Substrates

Mesenchymal cell crawling is a critical process in normal development, in tissue function, and in many diseases. Quantitatively predictive numerical simulations of cell crawling thus have multiple scientific, medical, and technological applications. However, we still lack a low-computational-cost approach to simulate mesenchymal three-dimensional (3D) cell crawling. Here, we develop a computationally tractable 3D model (implemented as a simulation in the CompuCell3D simulation environment) of mesenchymal cells crawling on a two-dimensional substrate. The Fürth equation, the usual characterization of mean-squared displacement (MSD) curves for migrating cells, describes a motion in which, for increasing time intervals, cell movement transitions from a ballistic to a diffusive regime. Recent experiments have shown that for very short time intervals, cells exhibit an additional fast diffusive regime. Our simulations' MSD curves reproduce the three experimentally observed temporal regimes, with fast diffusion for short time intervals, slow diffusion for long time intervals, and intermediate time -interval-ballistic motion. The resulting parameterization of the trajectories for both experiments and simulations allows the definition of time- and length scales that translate between computational and laboratory units. Rescaling by these scales allows direct quantitative comparisons among MSD curves and between velocity autocorrelation functions from experiments and simulations. Although our simulations replicate experimentally observed spontaneous symmetry breaking, short-timescale diffusive motion, and spontaneous cell-motion reorientation, their computational cost is low, allowing their use in multiscale virtual-tissue simulations. Comparisons between experimental and simulated cell motion support the hypothesis that short-time actomyosin dynamics affects longer-time cell motility. The success of the base cell-migration simulation model suggests its future application in more complex situations, including chemotaxis, migration through complex 3D matrices, and collective cell motion.

Heterogeneous T cell motility behaviors emerge from a coupling between speed and turning in vivo

T cells in vivo migrate primarily via undirected random walks, but it remains unresolved how these random walks generate an efficient search. Here, we use light sheet microscopy of T cells in the larval zebrafish as a model system to study motility across large populations of cells over hours in their native context. We show that cells do not perform Levy flight; rather, there is substantial cell-to-cell variability in speed, which persists over timespans of a few hours. This variability is amplified by a correlation between speed and directional persistence, generating a characteristic cell behavioral manifold that is preserved under a perturbation to cell speeds, and seen in Mouse T cells and Dictyostelium. Together, these effects generate a broad range of length scales over which cells explore in vivo.

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.

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?