Grid Computing for Sensitivity Analysis of Stochastic Biological Models

Systems biology is a multidisciplinary research area aimed at investigating biological systems by developing mathematical models that approach the study and the analysis of both the structure and behaviour of a biological phenomenon from a system perspective. The dynamics described by such mathematical models can be deeply affected by many parameters, and an extensive exploration of the parameters space in order to find crucial factors is most of the time prohibitive since it requires the execution of a huge number of computer simulations. Sensitivity analysis techniques can help in understanding how much the uncertainty in the model outcome is determined by the uncertainties, or by the variations, of the model input factors (components, reactions and respective parameters). In this work we exploit the European Grid Infrastructure to manage the calculations required to perform the SA on a stochastic model of bacterial chemotaxis, using an improved version of the first order screening method of Morris. According to the results achieved in our exploratory analysis, the European Grid Infrastructure is a useful solution for distributing the stochastic simulations required to carry out the SA of a stochastic model. Considering that the more intensive the computation the more scalable the infrastructure, grid computing can be a suitable technology for large scale biological models analysis.

[1]  Kevin Burrage,et al.  Stochastic approaches for modelling in vivo reactions , 2004, Comput. Biol. Chem..

[2]  Ian T. Foster,et al.  The Anatomy of the Grid: Enabling Scalable Virtual Organizations , 2001, Int. J. High Perform. Comput. Appl..

[3]  A. Sciaba,et al.  Experience with the gLite Workload Management System in ATLAS Monte Carlo production on LCG , 2008 .

[4]  G. Wadhams,et al.  Making sense of it all: bacterial chemotaxis , 2004, Nature Reviews Molecular Cell Biology.

[5]  A. D. Meglio,et al.  Programming the Grid with gLite , 2006 .

[6]  Giancarlo Mauri,et al.  Stochastic Simulations on a Grid Framework for Parameter Sweep Applications in Biological Models , 2009, 2009 International Workshop on High Performance Computational Systems Biology.

[7]  Yang Cao,et al.  Sensitivity analysis of discrete stochastic systems. , 2005, Biophysical journal.

[8]  Stephen Gilmore,et al.  Sensitivity Analysis of Stochastic Models of Bistable Biochemical Reactions , 2008, SFM.

[9]  Paolo Milazzo,et al.  The Calculus of Looping Sequences , 2008, SFM.

[10]  A. Arkin,et al.  Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. , 1998, Genetics.

[11]  Giancarlo Mauri,et al.  Tau Leaping Stochastic Simulation Method in P Systems , 2006, Workshop on Membrane Computing.

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

[13]  B. Stoddard,et al.  Mind your B's and R's: bacterial chemotaxis, signal transduction and protein recognition. , 1998, Structure.

[14]  Giancarlo Mauri,et al.  A study on the combined interplay between stochastic fluctuations and the number of flagella in bacterial chemotaxis , 2009, COMPMOD.

[15]  P. Swain,et al.  Stochastic Gene Expression in a Single Cell , 2002, Science.

[16]  Adam P. Arkin,et al.  Efficient stochastic sensitivity analysis of discrete event systems , 2007, J. Comput. Phys..

[17]  Andrea Saltelli,et al.  An effective screening design for sensitivity analysis of large models , 2007, Environ. Model. Softw..

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