Abstract We explore efficient estimation of statistical quantities, particularly rare event probabilities, for stochastic reaction networks. Consequently, we propose an importance sampling (IS) approach to improve the Monte Carlo (MC) estimator efficiency based on approximate tau-leap scheme. The crucial step in IS framework is choosing appropriate change probability measure achieve substantial...