A distributed software package for global sensitivity analysis of biological models

Mathematical modeling is a vital tool for studying biological systems. Due to the system complexity and technical challenges in molecular level measurement, it is commonly the case that a large number of model parameters have uncertain values. Analyzing model dynamics from a single estimated parameter set is insufficient and liable for misleading results. Global sensitivity analysis (GSA) has been recommended as a must-have step in the process of developing reliable models. However, the technique comes at high computation costs as it is based on Monte Carlo simulation which requires a large number of model evaluations and manipulating on massive data for sensitivity estimation. In this work, we develop a software package for global sensitivity estimation of biological models. The software is deployed on KISTI high performance computing (HPC) cluster environment to provide web-service to system biology modelers.

[1]  Sriram Krishnan,et al.  Opal web services for biomedical applications , 2010, Nucleic Acids Res..

[2]  A. Rundell,et al.  Comparative study of parameter sensitivity analyses of the TCR-activated Erk-MAPK signalling pathway. , 2006, Systems biology.

[3]  I. Sobol.  Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates , 2001 .

[4]  John G. Albeck,et al.  Systematic calibration of a cell signaling network model , 2010, BMC Bioinformatics.

[5]  Jacky L. Snoep,et al.  BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems , 2005, Nucleic Acids Res..

[6]  Zhike Zi,et al.  SBML-SAT: a systems biology markup language (SBML) based sensitivity analysis tool , 2008, BMC Bioinformatics.

[7]  H. Kitano,et al.  Computational systems biology , 2002, Nature.

[8]  M. Jansen.  Analysis of variance designs for model output , 1999 .

[9]  Satoru Shiono,et al.  Control mechanism of JAK/STAT signal transduction pathway , 2003, FEBS letters.

[10]  A. Hindmarsh,et al.  CVODE, a stiff/nonstiff ODE solver in C , 1996 .

[11]  Michael Hucka,et al.  LibSBML: an API Library for SBML , 2008, Bioinform..

[12]  Max D. Morris,et al.  Factorial sampling plans for preliminary computational experiments , 1991 .

[13]  D-M Hu,et al.  Time-dependent sensitivity analysis of biological networks: coupled MAPK and PI3K signal transduction pathways. , 2006, The journal of physical chemistry. A.

[14]  Kwang-Hyun Cho,et al.  Experimental Design in Systems Biology, Based on Parameter Sensitivity Analysis Using a Monte Carlo Method: A Case Study for the TNFα-Mediated NF-κ B Signal Transduction Pathway , 2003, Simul..

[15]  Stefano Tarantola,et al.  Winding Stairs: A sampling tool to compute sensitivity indices , 2000, Stat. Comput..

[16]  Mats Jirstrand,et al.  Systems biology Systems Biology Toolbox for MATLAB : a computational platform for research in systems biology , 2006 .

[17]  Z Zi,et al.  Sensitivity analysis approaches applied to systems biology models. , 2011, IET systems biology.

[18]  Saltelli Andrea,et al.  Global Sensitivity Analysis: The Primer , 2008 .