Steady-state parameter sensitivity in stochastic modeling via trajectory reweighting
- 12 March 2012
- journal article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 136 (10), 104106
- https://doi.org/10.1063/1.3690092
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. Comput. Phys. 221, 724 (2007)10.1016/j.jcp.2006.06.047]. 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.Keywords
This publication has 22 references indexed in Scilit:
- State-dependent biasing method for importance sampling in the weighted stochastic simulation algorithmThe Journal of Chemical Physics, 2010
- Enhanced Sampling of Nonequilibrium Steady StatesAnnual Review of Physical Chemistry, 2010
- Efficient computation of parameter sensitivities of discrete stochastic chemical reaction networksThe Journal of Chemical Physics, 2010
- Forward flux sampling for rare event simulationsJournal of Physics: Condensed Matter, 2009
- DNA looping provides stability and robustness to the bacteriophage λ switchProceedings of the National Academy of Sciences of the United States of America, 2009
- Spectral Methods for Parametric Sensitivity in Stochastic Dynamical SystemsBiophysical Journal, 2007
- Exact results for noise power spectra in linear biochemical reaction networksThe Journal of Chemical Physics, 2006
- COPASI—a COmplex PAthway SImulatorBioinformatics, 2006
- The finite state projection algorithm for the solution of the chemical master equationThe Journal of Chemical Physics, 2006
- Exact stochastic simulation of coupled chemical reactionsThe Journal of Physical Chemistry, 1977