Identifying diffusive motions in single-particle trajectories on the plasma membrane via fractional time-series models

Confined modes of single-particle trajectories induced by stochastic resetting

Random trajectories of single particles in living cells contain information about the interaction between particles, as well as with the cellular environment. However, precise consideration of the underlying stochastic properties, beyond normal diffusion, remains a challenge as applied to each particle trajectory separately. In this paper, we show how positions of confined particles in living cells can obey not only the Laplace distribution, but the Linnik one. This feature is detected in experimental data for the motion of G proteins and coupled receptors in cells, and its origin is explained in terms of stochastic resetting. This resetting process generates power-law waiting times, giving rise to the Linnik statistics in confined motion, and also includes exponentially distributed times as a limit case leading to the Laplace one. The stochastic process, which is affected by the resetting, can be Brownian motion commonly found in cells. Other possible models producing similar effects are discussed.

Single molecule characterization of anomalous transport in a thin, anisotropic film

The diffusion of small, charged molecules incorporated in an anisotropic polyelectrolyte multilayer (PEM) was tracked in three dimensions by combining single-molecule fluorescence localization (to characterize lateral diffusion) with Förster resonance energy transfer (FRET) between diffusing molecules and the supporting surface (to measure diffusion in the surface-normal direction). Analysis of the surface-normal diffusion required model-based statistical analysis to account for the inherently noisy FRET signal. Combining these distinct single-molecule methods, which are inherently sensitive to different length-scales, permitted simultaneous characterization of severely anisotropic diffusion, which was more than three orders of magnitude slower in the surface-normal direction. We hypothesize that the anomalously slow surface-normal diffusion was related to the periodic distribution of charge in the PEM, which created electrostatic barriers. The motion was strongly subdiffusive, with anomalous temporal scaling exponents in lateral and normal directions, suggesting a connection to the transient, random fractal conformation of polymer chains in the film's matrix.

An Estimation Algorithm for General Linear Single Particle Tracking Models with Time-Varying Parameters

Single Particle Tracking (SPT) is a powerful class of methods for studying the dynamics of biomolecules inside living cells. The techniques reveal the trajectories of individual particles, with a resolution well below the diffraction limit of light, and from them the parameters defining the motion model, such as diffusion coefficients and confinement lengths. Most existing algorithms assume these parameters are constant throughout an experiment. However, it has been demonstrated that they often vary with time as the tracked particles move through different regions in the cell or as conditions inside the cell change in response to stimuli. In this work, we propose an estimation algorithm to determine time-varying parameters of systems that discretely switch between different linear models of motion with Gaussian noise statistics, covering dynamics such as diffusion, directed motion, and Ornstein-Uhlenbeck dynamics. Our algorithm consists of three stages. In the first stage, we use a sliding window approach, combined with Expectation Maximization (EM) to determine maximum likelihood estimates of the parameters as a function of time. These results are only used to roughly estimate the number of model switches that occur in the data to guide the selection of algorithm parameters in the second stage. In the second stage, we use Change Detection (CD) techniques to identify where the models switch, taking advantage of the off-line nature of the analysis of SPT data to create non-causal algorithms with better precision than a purely causal approach. Finally, we apply EM to each set of data between the change points to determine final parameter estimates. We demonstrate our approach using experimental data generated in the lab under controlled conditions.