It is of primary interest for biologists to be able to visualize the dynamics of proteins within the cell. In this paper, we propose a new mapping method to robustly estimate dynamics in the entire cell from particle tracks. To obtain satisfying diffusion and drift maps, we use a spatiotemporal kernel estimator. Trajectory classification data is used as input and allows to automatically label particle movements into three classes: confined motion (or subdiffusion), Brownian motion, and directed motion (or superdiffusion). We then use this information to calculate diffusion coefficient and drift maps separately on each class of motion.

MOTIVATION: Recent advances in molecular biology and fluorescence microscopy imaging have made possible the inference of the dynamics of single molecules in living cells. Changes of dynamics can occur along a trajectory. Then, an issue is to estimate the temporal change-points that is the times at which a change of dynamics occurs. The number of points in the trajectory required to detect such changes will depend on both the magnitude and type of the motion changes. Here, the number of points per trajectory is of the order of 102, even if in practice dramatic motion changes can be detected with less points. RESULTS: We propose a non-parametric procedure based on test statistics computed on local windows along the trajectory to detect the change-points. This algorithm controls the number of false change-point detections in the case where the trajectory is fully Brownian. We also develop a strategy for aggregating the detections obtained with different window sizes so that the window size is no longer a parameter to optimize. A Monte Carlo study is proposed to demonstrate the performances of the method and also to compare the procedure to two competitive algorithms. At the end, we illustrate the efficacy of the method on real data in 2D and 3D, depicting the motion of mRNA complexes-called mRNA-binding proteins-in neuronal dendrites, Galectin-3 endocytosis and trafficking within the cell. AVAILABILITY AND IMPLEMENTATION: A user-friendly Matlab package containing examples and the code of the simulations used in the paper is available at http://serpico.rennes.inria.fr/doku.php? id=software:cpanalysis:index. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.