Modeling of 2D diffusion processes based on microscopy data: parameter estimation and practical identifiability analysis

BackgroundDiffusion is a key component of many biological processes such as chemotaxis, developmental differentiation and tissue morphogenesis. Since recently, the spatial gradients caused by diffusion can be assessed in-vitro and in-vivo using microscopy based imaging techniques. The resulting time-series of two dimensional, high-resolutions images in combination with mechanistic models enable the quantitative analysis of the underlying mechanisms. However, such a model-based analysis is still challenging due to measurement noise and sparse observations, which result in uncertainties of the model parameters.MethodsWe introduce a likelihood function for image-based measurements with log-normal distributed noise. Based upon this likelihood function we formulate the maximum likelihood estimation problem, which is solved using PDE-constrained optimization methods. To assess the uncertainty and practical identifiability of the parameters we introduce profile likelihoods for diffusion processes.Results and conclusionAs proof of concept, we model certain aspects of the guidance of dendritic cells towards lymphatic vessels, an example for haptotaxis. Using a realistic set of artificial measurement data, we estimate the five kinetic parameters of this model and compute profile likelihoods. Our novel approach for the estimation of model parameters from image data as well as the proposed identifiability analysis approach is widely applicable to diffusion processes. The profile likelihood based method provides more rigorous uncertainty bounds in contrast to local approximation methods.

1. 1
   The chemical basis of morphogenesis

   Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences195211,630 citationsDataRank 1.4
2. 2
   Uncertainty Analysis for Non-identifiable Dynamical Systems: Profile Likelihoods, Bootstrapping and More

   Lecture Notes in Computer Science201439 citationsDataRank 0.553
3. 3
   Tailored parameter optimization methods for ordinary differential equation models with steady-state constraints

   BMC Systems Biology201628 citationsDataRank 0.505
4. 4
   Radial Basis Function Approximations of Bayesian Parameter Posterior Densities for Uncertainty Analysis

   Lecture Notes in Computer Science20143 citationsDataRank 0.208
5. 5
   Quantitative comparison of competing PDE models for Pomlp dynamics in fission yeast * *The authors acknowledge financial support from the Postdoctoral Fellowship Program (PFP) of the Helmholtz Zentrum Munchen.

   IFAC-PapersOnLine20162 citationsDataRank 0.165