Repository logo
 

Limits of Accuracy for Parameter Estimation and Localization in Single-Molecule Microscopy via Sequential Monte Carlo Methods

Accepted version
Peer-reviewed

Type

Article

Change log

Authors

Marie D’Avigneau, A 
Singh, SS 
Ober, RJ 

Abstract

Assessing the quality of parameter estimates for models describing the motion of single molecules in cellular environments is an important problem in fluorescence microscopy. In this work, we consider the fundamental data model, where molecules emit photons at random time instances and these photons arrive at random locations on the detector ac- cording to complex point spread functions (PSFs). The randomness and non-Gaussian PSF of the detection process, and the random trajectory of the molecule, makes inference challenging. Moreover, the presence of other closely spaced molecules causes further un- certainty in the origin of the measurements, which impacts the statistical precision of the estimates. We quantify the limits of accuracy of model parameter estimates and separation distance between closely spaced molecules (known as the resolution problem) by computing the Cram ́er-Rao lower bound (CRLB), or equivalently the inverse of the Fisher informa- tion matrix (FIM), for the variance of estimates. Results on the CRLB obtained from the fundamental model are crucial, in that they provide a lower bound for more practical scenar- ios. While analytic expressions for the FIM can be derived for static and deterministically moving molecules, the analytical tools to evaluate the FIM for molecules whose trajecto- ries follow stochastic differential equations (SDEs) are still for the most part missing. We address this by presenting a general sequential Monte Carlo (SMC) based methodology for both parameter inference and computing the desired accuracy limits for non-static molecules and a non-Gaussian fundamental detection model. For the first time, we are able to estimate the FIM for stochastically moving molecules observed through the Airy and Born and Wolf detection models. This is achieved by estimating the score and observed information matrix via SMC. We summarise the outcome of our numerical work by delineating the qualitative behaviours for the accuracy limits as functions of various experimental settings like collected photon count, molecule diffusion, etc. We also verify that we can recover known results from the static molecule case.

Description

Keywords

&nbsp, single-molecule microscopy, sequential Monte Carlo, SMC, Fisher information matrix, fluorescence, microscopy, particle smoothing, stochastic differential equations, SDEs, particle filtering

Journal Title

SIAM Journal on Imaging Sciences

Conference Name

Journal ISSN

1936-4954
1936-4954

Volume Title

Publisher

Society for Industrial & Applied Mathematics (SIAM)

Rights

All rights reserved