Steady-state parameter sensitivity in stochastic modeling via trajectory reweighting

Abstract

Parameter sensitivity analysis is a powerful tool in the building and analysis of biochemical network models. For stochastic simulations, parameter sensitivity analysis can be computationally expensive, requiring multiple simulations for perturbed values of the parameters. Here, we use trajectory reweighting to derive a method for computing sensitivity coefficients in stochastic simulations without explicitly perturbing the parameter values, avoiding the need for repeated simulations. The method allows the simultaneous computation of multiple sensitivity coefficients. Our approach recovers results originally obtained by application of the Girsanov measure transform in the general theory of stochastic processes [A. Plyasunov and A. P. Arkin, J. Comp. Phys. 221, 724 (2007)]. We build on these results to show how the method can be used to compute steady-state sensitivity coefficients from a single simulation run, and we present various efficiency improvements. For models of biochemical signaling networks the method has a particularly simple implementation. We demonstrate its application to a signaling network showing stochastic focussing and to a bistable genetic switch, and present exact results for models with linear propensity functions.

17 pages (including supplementary material); 5 figures

References in corpus (1)

• Steady-state parameter sensitivity in stochastic modeling via trajectory reweighting

Cited by in corpus (12)

• Markovian Dynamics on Complex Reaction Networks
• Steady-state parameter sensitivity in stochastic modeling via trajectory reweighting
• A practical approach to the sensitivity analysis for kinetic Monte Carlo simulation of heterogeneous catalysis
• Malliavin weight sampling: a practical guide
• Efficient Finite Difference Method for Computing Sensitivities of Biochemical Reactions
• Efficient estimators for likelihood ratio sensitivity indices of complex stochastic dynamics
• Stochastic Averaging and Sensitivity Analysis for Two Scale Reaction Networks
• Trajectory reweighting for non-equilibrium steady states
• An efficient and unbiased method for sensitivity analysis of stochastic reaction networks
• Variance reduction for steady-state simulation and sensitivity analysis of stochastic chemical systems
• Parametric Sensitivity Analysis for Stochastic Molecular Systems using Information Theoretic Metrics
• Steady State Sensitivity Analysis of Continuous Time Markov Chains