Edinburgh Research Explorer

Parameter estimation for biochemical reaction networks using Wasserstein distances

Research output: Contribution to journalArticle

Original languageEnglish
Number of pages24
JournalJournal of Physics A: Mathematical and Theoretical
DOIs
Publication statusAccepted/In press - 18 Nov 2019

Abstract

We present a method for estimating parameters in stochastic models of biochemical reaction networks by fitting steady-state distributions using Wasserstein distances. We simulate a reaction network at different parameter settings and train a Gaussian process to learn the Wasserstein distance between observations and the simulator output for all parameters. We then use Bayesian optimization to find parameters minimizing this distance based on the trained Gaussian process. The effectiveness of our method is demonstrated on the three-stage model of gene expression and a genetic feedback loop for which moment-based methods are known to perform poorly. Our method is applicable to any simulator model of stochastic reaction networks, including Brownian Dynamics.

    Research areas

  • Wasserstein distance, Bayesian optimization, Chemical Master Equation, parameter estimation

ID: 121352898