Diffusion analysis of single particle trajectories in a Bayesian nonparametrics framework

Directional change during active diffusion of viral ribonucleoprotein particles through cytoplasm

A mesh of cytoskeletal fibers, consisting of microtubules, intermediate filaments, and fibrous actin, prevents the Brownian diffusion of particles with a diameter larger than 0.10 μm, such as vesicular stomatitis virus ribonucleoprotein (RNP) particles, in mammalian cells. Nevertheless, RNP particles do move in random directions but at a lower rate than Brownian diffusion, which is thermally driven. This nonthermal biological transport process is called "active diffusion" because it is driven by ATP. The ATP powers motor proteins such as myosin II. The motor proteins bend and cross-link actin fibers, causing the mesh to jiggle. Until recently, little was known about how RNP particles get through the mesh. It has been customary to analyze the tracks of particles like RNPs by computing the slope of the ensemble-averaged mean-squared displacement of the particles as a signature of mechanism. Although widely used, this approach "loses information" about the timing of the switches between physical mechanisms. It has been recently shown that machine learning composed of variational Bayesian analysis, Gaussian mixture models, and hidden Markov models can use "all the information" in a single track to reveal that that the positions of RNP particles are spatially clustered. Machine learning assigns a number, called a state, to each cluster. RNP particles remain in one state for 0.2-1.0 s before switching (hopping) to a different state. This earlier work is here extended to analyze the movements of a particle within a state and to determine particle directionality within and between states.

Trajectory Analysis in Single-Particle Tracking: From Mean Squared Displacement to Machine Learning Approaches

Single-particle tracking is a powerful technique to investigate the motion of molecules or particles. Here, we review the methods for analyzing the reconstructed trajectories, a fundamental step for deciphering the underlying mechanisms driving the motion. First, we review the traditional analysis based on the mean squared displacement (MSD), highlighting the sometimes-neglected factors potentially affecting the accuracy of the results. We then report methods that exploit the distribution of parameters other than displacements, e.g., angles, velocities, and times and probabilities of reaching a target, discussing how they are more sensitive in characterizing heterogeneities and transient behaviors masked in the MSD analysis. Hidden Markov Models are also used for this purpose, and these allow for the identification of different states, their populations and the switching kinetics. Finally, we discuss a rapidly expanding field-trajectory analysis based on machine learning. Various approaches, from random forest to deep learning, are used to classify trajectory motions, which can be identified by motion models or by model-free sets of trajectory features, either previously defined or automatically identified by the algorithms. We also review free software available for some of the analysis methods. We emphasize that approaches based on a combination of the different methods, including classical statistics and machine learning, may be the way to obtain the most informative and accurate results.

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions

Many imaging techniques for biological systems-like fixation of cells coupled with fluorescence microscopy-provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. These snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

Hidden Markov modeling of single-particle diffusion with stochastic tethering

The statistics of the diffusive motion of particles often serve as an experimental proxy for their interaction with the environment. However, inferring the physical properties from the observed trajectories is challenging. Inspired by a recent experiment, here we analyze the problem of particles undergoing two-dimensional Brownian motion with transient tethering to the surface. We model the problem as a hidden Markov model where the physical position is observed and the tethering state is hidden. We develop an alternating maximization algorithm to infer the hidden state of the particle and estimate the physical parameters of the system. The crux of our method is a saddle-point-like approximation, which involves finding the most likely sequence of hidden states and estimating the physical parameters from it. Extensive numerical tests demonstrate that our algorithm reliably finds the model parameters and is insensitive to the initial guess. We discuss the different regimes of physical parameters and the algorithm's performance in these regimes. We also provide a free software implementation of our algorithm.

ExTrack characterizes transition kinetics and diffusion in noisy single-particle tracks

Single-particle tracking microscopy is a powerful technique to investigate how proteins dynamically interact with their environment in live cells. However, the analysis of tracks is confounded by noisy molecule localization, short tracks, and rapid transitions between different motion states, notably between immobile and diffusive states. Here, we propose a probabilistic method termed ExTrack that uses the full spatio-temporal information of tracks to extract global model parameters, to calculate state probabilities at every time point, to reveal distributions of state durations, and to refine the positions of bound molecules. ExTrack works for a wide range of diffusion coefficients and transition rates, even if experimental data deviate from model assumptions. We demonstrate its capacity by applying it to slowly diffusing and rapidly transitioning bacterial envelope proteins. ExTrack greatly increases the regime of computationally analyzable noisy single-particle tracks. The ExTrack package is available in ImageJ and Python.

ExTrack characterizes transition kinetics and diffusion in noisy single-particle tracks

Single-particle tracking microscopy is a powerful technique to investigate how proteins dynamically interact with their environment in live cells. However, the analysis of tracks is confounded by noisy molecule localization, short tracks, and rapid transitions between different motion states, notably between immobile and diffusive states. Here, we propose a probabilistic method termed ExTrack that uses the full spatio-temporal information of tracks to extract global model parameters, to calculate state probabilities at every time point, to reveal distributions of state durations, and to refine the positions of bound molecules. ExTrack works for a wide range of diffusion coefficients and transition rates, even if experimental data deviate from model assumptions. We demonstrate its capacity by applying it to slowly diffusing and rapidly transitioning bacterial envelope proteins. ExTrack greatly increases the regime of computationally analyzable noisy single-particle tracks. The ExTrack package is available in ImageJ and Python.

Receptor tyrosine kinase MET ligand-interaction classified via machine learning from single-particle tracking data

Internalin B-mediated activation of the membrane-bound receptor tyrosine kinase MET is accompanied by a change in receptor mobility. Conversely, it should be possible to infer from receptor mobility whether a cell has been treated with internalin B. Here, we propose a method based on hidden Markov modeling and explainable artificial intelligence that machine-learns the key differences in MET mobility between internalin B-treated and -untreated cells from single-particle tracking data. Our method assigns receptor mobility to three diffusion modes (immobile, slow, and fast). It discriminates between internalin B-treated and -untreated cells with a balanced accuracy of >99% and identifies three parameters that are most affected by internalin B treatment: a decrease in the mobility of slow molecules (1) and a depopulation of the fast mode (2) caused by an increased transition of fast molecules to the slow mode (3). Our approach is based entirely on free software and is readily applicable to the analysis of other membrane receptors.

Biophysical Models of PAR Cluster Transport by Cortical Flow in C. elegans Early Embryogenesis

The clustering of membrane-bound proteins facilitates their transport by cortical actin flow in early Caenorhabditis elegans embryo cell polarity. PAR-3 clustering is critical for this process, yet the biophysical processes that couple protein clusters to cortical flow remain unknown. We develop a discrete, stochastic agent-based model of protein clustering and test four hypothetical models for how clusters may interact with the flow. Results show that the canonical way to assess transport characteristics from single-particle tracking data used thus far in this area, the Péclet number, is insufficient to distinguish these hypotheses and that all models can account for transport characteristics quantified by this measure. However, using this model, we demonstrate that these different cluster-cortex interactions may be distinguished using a different metric, namely the scalar projection of cluster displacement on to the flow displacement vector. Our results thus provide a testable way to use existing single-particle tracking data to test how endogenous protein clusters may interact with the cortical flow to localize during polarity establishment. To facilitate this investigation, we also develop both improved simulation and semi-analytic methodologies to quantify motion summary statistics (e.g., Péclet number and scalar projection) for these stochastic models as a function of biophysical parameters.

Delivery of Ultrasmall Nanoparticles to the Cytosolic Compartment of Pyroptotic J774A.1 Macrophages via GSDMDNterm Membrane Pores

Endosome capture is a major physiological barrier to the successful delivery of nanomedicine. Here, we found a strategy to deliver ultrasmall nanoparticles (<10 nm) to the cytosolic compartment of pyroptotic cells with spontaneous endosomal escape. To mimic pathological pyroptotic cells, J774A.1 macrophages were stimulated with lipopolysaccharide (LPS) plus nigericin (Nig) or adenosine triphosphate (ATP) to form specific gasdermin D protein-driven membrane pores at an N-terminal domain (GSDMD^Nterm). Through GSDMD^Nterm membrane pores, both anionic and cationic nanoparticles (NPs) with diameters less than 10 nm were accessed into the cytosolic compartment of pyroptotic cells in an energy- and receptor-independent manner, while NPs larger than the size of GSDMD^Nterm membrane pores failed to enter pyroptotic cells. NPs pass through GSDMD^Nterm membrane pores via free diffusion and then access into the cytoplasm of pyroptotic cells in a microtubule-independent manner. Interestingly, we found that LPS-primed NPs may act as Trojan horse, deliver extracellular LPS into normal cells through endocytosis, and in turn induce GSDMD^Nterm membrane pores, which facilitate further internalization of NPs. This study presented a straightforward method of distinguishing normal and pyroptotic cells through GSDMD membrane pores, implicating their potential application in monitoring the delivery of desired nanomedicines in pyroptosis-related diseases and conditions.

Uncover Single Nanoparticle Dynamics on Live Cell Membrane with Data-Driven Historical Experience Analysis

Understanding the spatiotemporal dynamics of particles in a complex biological environment is crucial for the study of related biological processes. To analyze the complicated trajectories recorded from single-particle tracking (SPT), we have proposed a method named SEES based on historical experience vector analysis, which allows both the global patterns and local state continuities of a trajectory to emerge by themselves as color segments without predefined models. This method implements a data-driven strategy and thus uncovers the hidden information with less prior knowledge or subjective bias. Here, we demonstrate its efficiency by comparing its performance with the Hidden Markov model (HMM), one of the most widely used methods in time series processing. The results demonstrated that the SEES operator was more sensitive in identifying rare events and could utilize multivariable observations in the dynamic processes to uncover more details. We applied the method to analyze the dynamics of nanoparticles interacting with live cells expressing programmed death ligand 1 (PD-L1) on the membrane. The results showed that the SEES operator can successfully pinpoint the transmembrane rare events, visualize the on-membrane "Brownian searching" motion, and evaluate different dynamics among multiple trajectories. Furthermore, we found that the PD-L1 expression level on the cell membrane affected the rotation behavior of the nanoparticle as well as the cellular uptake efficiency. These findings enabled by SEES could potentially help the rational design of highly efficient nanocargoes.

Maximum likelihood estimates of diffusion coefficients from single-particle tracking experiments

Single-molecule localization microscopy allows practitioners to locate and track labeled molecules in biological systems. When extracting diffusion coefficients from the resulting trajectories, it is common practice to perform a linear fit on mean-squared-displacement curves. However, this strategy is suboptimal and prone to errors. Recently, it was shown that the increments between the observed positions provide a good estimate for the diffusion coefficient, and their statistics are well-suited for likelihood-based analysis methods. Here, we revisit the problem of extracting diffusion coefficients from single-particle tracking experiments subject to static noise and dynamic motion blur using the principle of maximum likelihood. Taking advantage of an efficient real-space formulation, we extend the model to mixtures of subpopulations differing in their diffusion coefficients, which we estimate with the help of the expectation-maximization algorithm. This formulation naturally leads to a probabilistic assignment of trajectories to subpopulations. We employ the theory to analyze experimental tracking data that cannot be explained with a single diffusion coefficient. We test how well a dataset conforms to the assumptions of a diffusion model and determine the optimal number of subpopulations with the help of a quality factor of known analytical distribution. To facilitate use by practitioners, we provide a fast open-source implementation of the theory for the efficient analysis of multiple trajectories in arbitrary dimensions simultaneously.