Working Paper Article Version 1 This version is not peer-reviewed

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

Version 1 : Received: 22 January 2021 / Approved: 25 January 2021 / Online: 25 January 2021 (12:26:31 CET)

A peer-reviewed article of this Preprint also exists.

Godoy, B.I.; Vickers, N.A.; Andersson, S.B. An Estimation Algorithm for General Linear Single Particle Tracking Models with Time-Varying Parameters. Molecules 2021, 26, 886. Godoy, B.I.; Vickers, N.A.; Andersson, S.B. An Estimation Algorithm for General Linear Single Particle Tracking Models with Time-Varying Parameters. Molecules 2021, 26, 886.

Journal reference: Molecules 2021, 26, 886
DOI: 10.3390/molecules26040886

Abstract

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.

Subject Areas

single particle tracking; single molecule biophysics; fluorescence

Comments (0)

We encourage comments and feedback from a broad range of readers. See criteria for comments and our diversity statement.

Leave a public comment
Send a private comment to the author(s)
Views 0
Downloads 0
Comments 0
Metrics 0


×
Alerts
Notify me about updates to this article or when a peer-reviewed version is published.
We use cookies on our website to ensure you get the best experience.
Read more about our cookies here.