i c J o

Sensitivity analysis for stochastic chemical reaction networks with multiple time-scales Ankit Gupta ∗ Mustafa Khammash∗ Stochastic models for chemical reaction networks have become very popular in re-cent years. For such models, the estimation of parameter sensitivities is an impor-tant and challenging problem. Sensitivity values help in analyzing the network, un-derstanding its robustness properties and also in identifying the key reactions for a given outcome. Most of the methods that exist in the literature for the estimation of parameter sensitivities, rely on Monte Carlo simulations using Gillespie’s stochastic simulation algorithm or its variants. It is well-known that such simulation methods can be prohibitively expensive when the network contains reactions firing at differ-ent time-scales, which is a feature of many important biochemical networks. For such networks, it is often possible to exploit the time-scale separation and approximately capture the original dynamics by simulating a “reduced " model, which is obtained by eliminating the fast reactions in a certain way. The aim of this paper is to tie these model reduction techniques with sensitivity analysis. We prove that under some con-ditions, the sensitivity values for the reduced model can be used to approximately re-cover the sensitivity values for the original model. Through an example we illustrate how our result can help in sharply reducing the computational costs for the estima-tion of parameter sensitivities for reaction networks with multiple time-scales. To prove our result, we use coupling arguments based on the random time change rep-resentation of Kurtz. We also exploit certain connections between the distributions of the occupation times of Markov chains and multi-dimensional wave equations