ABSTRACT
In this paper, we aim to detect trapping areas (equivalently microdomains or confinement areas) within cells, corresponding to regions where molecules are trapped and thereby undergo subdiffusion. We propose an original computational approach that takes as input a set of molecule trajectories estimated by appropriate tracking methods. The core of the algorithm is based on a combination of clustering algorithms with trajectory classification procedures able to distinguish subdiffusion, superdiffusion and Brownian motion. The idea is to automatically identify trapping areas where we observe a high concentration of subdiffusive particles. We evaluate our proof of concept on artificial sequences obtained with a biophysics-based simulator (Fluosim), and we illustrate its potential on real TIRF microscopy data.
Subject(s)
Algorithms , Cells/chemistry , Cells/metabolism , Computer Simulation , Cell Movement , Cells/cytology , Diffusion , Models, MolecularABSTRACT
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.