A comparison of approximation techniques for variance-based sensitivity analysis of biochemical reaction systems.
BACKGROUND
Sensitivity analysis is an indispensable tool for the analysis of complex systems. In a recent paper, we have introduced a thermodynamically consistent variance-based sensitivity analysis approach for studying the robustness and fragility properties of biochemical reaction systems under uncertainty in the standard chemical potentials of the activated complexes of the reactions and the standard chemical potentials of the molecular species. In that approach, key sensitivity indices were estimated by Monte Carlo sampling, which is computationally very demanding and impractical for large biochemical reaction systems. Computationally efficient algorithms are needed to make variance-based sensitivity analysis applicable to realistic cellular networks, modeled by biochemical reaction systems that consist of a large number of reactions and molecular species.


RESULTS
We present four techniques, derivative approximation (DA), polynomial approximation (PA), Gauss-Hermite integration (GHI), and orthonormal Hermite approximation (OHA), for analytically approximating the variance-based sensitivity indices associated with a biochemical reaction system. By using a well-known model of the mitogen-activated protein kinase signaling cascade as a case study, we numerically compare the approximation quality of these techniques against traditional Monte Carlo sampling. Our results indicate that, although DA is computationally the most attractive technique, special care should be exercised when using it for sensitivity analysis, since it may only be accurate at low levels of uncertainty. On the other hand, PA, GHI, and OHA are computationally more demanding than DA but can work well at high levels of uncertainty. GHI results in a slightly better accuracy than PA, but it is more difficult to implement. OHA produces the most accurate approximation results and can be implemented in a straightforward manner. It turns out that the computational cost of the four approximation techniques considered in this paper is orders of magnitude smaller than traditional Monte Carlo estimation. Software, coded in MATLAB, which implements all sensitivity analysis techniques discussed in this paper, is available free of charge.


CONCLUSIONS
Estimating variance-based sensitivity indices of a large biochemical reaction system is a computationally challenging task that can only be addressed via approximations. Among the methods presented in this paper, a technique based on orthonormal Hermite polynomials seems to be an acceptable candidate for the job, producing very good approximation results for a wide range of uncertainty levels in a fraction of the time required by traditional Monte Carlo sampling.
Zhang
Hong-Xuan
HX
Whitaker Biomedical Engineering Institute, The Johns Hopkins University, Baltimore, MD 21218, USA.
Goutsias
John
J
D000465:Algorithms; D000704:Analysis of Variance; D001669:Biochemical Phenomena; D009010:Monte Carlo Method; D015398:Signal Transduction
Oncogene. 2005 Aug 25;24(36):5533-42
16007170
J Phys Chem A. 2006 Feb 23;110(7):2474-85
16480307
Bioinformatics. 2009 Oct 1;25(19):2466-72
19648136
Chem Rev. 2005 Jul;105(7):2811-28
16011325
J Biopharm Stat. 1995 Nov;5(3):245-71
8580927
Proc Natl Acad Sci U S A. 2004 Sep 7;101(36):13210-5
15340155
Cell. 1995 Jan 27;80(2):179-85
7834738
Clin Pharmacokinet. 2003;42(10):883-908
12885263
PLoS Comput Biol. 2007 Jul;3(7):e130
17616983
J Phys Chem A. 2006 Apr 27;110(16):5361-70
16623463
Biophys J. 2004 Oct;87(4):2195-202
15454422
Nat Cell Biol. 2002 Aug;4(8):556-64
12134156
Neural Comput. 2007 Jan;19(1):231-57
17134324
Biophys J. 2004 Nov;87(5):L01-2
15454405
J Theor Biol. 2004 Oct 21;230(4):541-62
15363675
J Cell Sci. 2000 Dec;113 Pt 24:4499-510
11082043
Biophys J. 2007 Mar 15;92(6):1846-57
17208985
J Cell Sci. 2005 Jul 15;118(Pt 14):2997-3002
16014377
Mol Cell Biol. 2004 Jan;24(1):144-53
14673150
J Chem Phys. 2009 Sep 7;131(9):094101
19739843
Biochem J. 1998 Mar 15;330 ( Pt 3):1451-60
9494119