Optimization of the shape (and topology) of the initial conditions for diffusion parameter identification

The design of an experiment, e.g., the setting of initial conditions, strongly influences the accuracy of the whole process of determining model parameters from data. We impose a sensitivity-based approach for choosing optimal design variables and study the optimization of the shape (and topology) of the initial conditions for an inverse problem of a diffusion parameter identification. Our approach, although case independent, is illustrated at the FRAP (Fluorescence Recovery After Photobleaching) experimental technique. The core idea resides in the maximization of a sensitivity measure, which depends on a specific experimental setting of initial conditions. By a numerical optimization, we find an interesting pattern of increasingly complicated (with respect to connectivity) optimal initial shapes. The proposed modification of the FRAP experimental protocol is rather surprising but entirely realistic and the resulting enhancement of the parameter estimate accuracy is significant.

[1]  Levi A. Gheber,et al.  Universal Approach to FRAP Analysis of Arbitrary Bleaching Patterns , 2015 .

[2]  Benn Macdonald,et al.  Closed-form formulas vs. PDE based numerical solution for the FRAP data processing: Theoretical and practical comparison , 2017, Comput. Math. Appl..

[3]  Stefan Kindermann,et al.  On Data Space Selection and Data Processing for Parameter Identification in a Reaction-Diffusion Model Based on FRAP Experiments , 2015 .

[4]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[5]  W. Webb,et al.  Mobility measurement by analysis of fluorescence photobleaching recovery kinetics. , 1976, Biophysical journal.

[6]  Ivo F. Sbalzarini,et al.  Analysis, modeling, and simulation of diffusion processes in cell biology , 2006 .

[7]  Joel A. Tropp,et al.  Topics in sparse approximation , 2004 .

[8]  C. Mullineaux,et al.  Mobility of photosynthetic complexes in thylakoid membranes , 1997, Nature.

[9]  B. Vandeginste Nonlinear regression analysis: Its applications, D. M. Bates and D. G. Watts, Wiley, New York, 1988. ISBN 0471‐816434. Price: £34.50 , 1989 .

[10]  Stefan Kindermann,et al.  FRAP & FLIP: Two Sides of the Same Coin? , 2015, IWBBIO.

[11]  Davide Mazza,et al.  FRAP and kinetic modeling in the analysis of nuclear protein dynamics: what do we really know? , 2010, Current opinion in cell biology.

[12]  J. Rohn Computing the norm ∥A∥∞,1 is NP-hard , 2000 .

[13]  Ctirad Matonoha,et al.  Phycobilisome Mobility and Its Role in the Regulation of Light Harvesting in Red Algae1[W] , 2014, Plant Physiology.

[14]  Howard J. Worman,et al.  Nuclear Membrane Dynamics and Reassembly in Living Cells: Targeting of an Inner Nuclear Membrane Protein in Interphase and Mitosis , 1997, The Journal of cell biology.

[15]  Sabine Attinger,et al.  Are assumptions about the model type necessary in reaction-diffusion modeling? A FRAP application. , 2011, Biophysical journal.

[16]  Stepán Papácek,et al.  Estimation of diffusivity of phycobilisomes on thylakoid membrane based on spatio-temporal FRAP images , 2013, Math. Comput. Model..

[17]  Douglas M. Bates,et al.  Nonlinear Regression Analysis and Its Applications , 1988 .