Moment fitting for parameter inference in repeatedly and partially observed stochastic biological models.
The inference of reaction rate parameters in biochemical network models from time series concentration data is a central task in computational systems biology. Under the assumption of well mixed conditions the network dynamics are typically described by the chemical master equation, the Fokker Planck equation, the linear noise approximation or the macroscopic rate equation. The inverse problem of estimating the parameters of the underlying network model can be approached in deterministic and stochastic ways, and available methods often compare individual or mean concentration traces obtained from experiments with theoretical model predictions when maximizing likelihoods, minimizing regularized least squares functionals, approximating posterior distributions or sequentially processing the data. In this article we assume that the biological reaction network can be observed at least partially and repeatedly over time such that sample moments of species molecule numbers for various time points can be calculated from the data. Based on the chemical master equation we furthermore derive closed systems of parameter dependent nonlinear ordinary differential equations that predict the time evolution of the statistical moments. For inferring the reaction rate parameters we suggest to not only compare the sample mean with the theoretical mean prediction but also to take the residual of higher order moments explicitly into account. Cost functions that involve residuals of higher order moments may form landscapes in the parameter space that have more pronounced curvatures at the minimizer and hence may weaken or even overcome parameter sloppiness and uncertainty. As a consequence both deterministic and stochastic parameter inference algorithms may be improved with respect to accuracy and efficiency. We demonstrate the potential of moment fitting for parameter inference by means of illustrative stochastic biological models from the literature and address topics for future research.
Mathematical Methods in Molecular and Systems Biology, Johann Radon Institute for Computational and Applied Mathematics, Linz, Austria. philipp.kuegler@oeaw.ac.at
Philipp
P
Kügler
D000465:Algorithms; D003198:Computer Simulation; D019281:Dimerization; D007700:Kinetics; D008954:Models, Biological; D015398:Signal Transduction; D013269:Stochastic Processes; D049490:Systems Biology; D016159:Tumor Suppressor Protein p53 / Q000378:metabolism
D016159:Tumor Suppressor Protein p53