ENISI SDE: A novel web-based stochastic modeling tool for computational biology

Stochasticity is part of the nature of many biological processes and an important aspect in modeling and simulations of computational biology. Gillespie's algorithm was developed for modeling stochasticity of chemical reactions and has been broadly applied to stochastic modeling in computational biology. The Gillespie's algorithm accounts for the particle effect in chemical reactions. However, many biological processes including cellular and molecular immunological mechanisms have many other sources of stochasticity such as cell movement, ligand binding, or unaccounted variation in experimental settings. In this paper, we propose stochastic differential equations (SDE) being used as a generic stochastic modeling technique for systems immunology. SDE has been widely used in statistics and economics areas; but only a few isolated studies in computational biology are found using SDE to model stochastic behaviors of cells and molecules. In addition, to the best of our knowledge, there is no user-friendly SDE-based modeling tool available for computational biologists. This paper presents ENISI SDE, a web-based user-friendly stochastic modeling tool for computational biologists. This work provides three major contributions: (1) we discuss SDE and propose it as a generic approach for stochastic modeling in computational biology; (2) we develop ENISI SDE, a web-based user-friendly SDE modeling tool that only requires little extra effort beyond regular ODE-based modeling; (3) we use the model SDE modeling tool to study stochastic sources of cell heterogeneity in the context of a CD4+ T Cell differentiation process. The case study clearly shows the effectiveness of SDE as a stochastic modeling approach in biology in general and immunology in particular and the power of the SDE modeling tool we developed.

[1]  Intawat Nookaew,et al.  Prospects for systems biology and modeling of the gut microbiome. , 2011, Trends in biotechnology.

[2]  Avi Ma’ayan,et al.  Systems biology of stem cell fate and cellular reprogramming , 2009, Nature Reviews Molecular Cell Biology.

[3]  H. Weiner,et al.  Reciprocal developmental pathways for the generation of pathogenic effector TH17 and regulatory T cells , 2006, Nature.

[4]  Mike Holcombe,et al.  Development of a Three Dimensional Multiscale Computational Model of the Human Epidermis , 2010, PloS one.

[5]  Olli Yli-Harja,et al.  Stochastic Differential Equation Model for Cerebellar Granule Cell Excitability , 2008, PLoS Comput. Biol..

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

[7]  Harvey J. Gold,et al.  Mathematical modeling of biological systems. An introductory guidebook. , 1977 .

[8]  Hidde de Jong,et al.  Modeling and Simulation of Genetic Regulatory Systems: A Literature Review , 2002, J. Comput. Biol..

[9]  S. Basov Simulation and Inference for Stochastic Differential Equations: With R Examples , 2010 .

[10]  Cheng-Yan Kao,et al.  A stochastic differential equation model for quantifying transcriptional regulatory network in Saccharomyces cerevisiae , 2005, Bioinform..

[11]  J. Tyson,et al.  Modeling the control of DNA replication in fission yeast. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Elhanan Borenstein,et al.  Computational systems biology and in silico modeling of the human microbiome , 2012, Briefings Bioinform..

[13]  Ivo F Sbalzarini,et al.  A partial-propensity formulation of the stochastic simulation algorithm for chemical reaction networks with delays. , 2011, The Journal of chemical physics.

[14]  Chen Dong,et al.  Molecular antagonism and plasticity of regulatory and inflammatory T cell programs. , 2008, Immunity.

[15]  Nicolas Le Novère,et al.  STOCHSIM: modelling of stochastic biomolecular processes , 2001, Bioinform..

[16]  D. Gillespie A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[17]  Xiaodong Cai,et al.  Efficient exact and K-skip methods for stochastic simulation of coupled chemical reactions. , 2009, The Journal of chemical physics.

[18]  Mudita Singhal,et al.  COPASI - a COmplex PAthway SImulator , 2006, Bioinform..

[19]  Hiroaki Kitano,et al.  The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models , 2003, Bioinform..

[20]  Z Liu,et al.  Detailed comparison between StochSim and SSA. , 2008, IET systems biology.

[21]  J. Morgan,et al.  Mathematical modeling of plant metabolic pathways. , 2002, Metabolic engineering.

[22]  Stephanie Forrest,et al.  Infect Recognize Destroy , 1996 .

[23]  J. Tyson,et al.  Mathematical model of the cell division cycle of fission yeast. , 2001, Chaos.

[24]  Linda R Petzold,et al.  The slow-scale stochastic simulation algorithm. , 2005, The Journal of chemical physics.

[25]  T. McMillen Simulation and Inference for Stochastic Differential Equations: With R Examples , 2008 .

[26]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[27]  Madhav V. Marathe,et al.  Systems Modeling of Molecular Mechanisms Controlling Cytokine-driven CD4+ T Cell Differentiation and Phenotype Plasticity , 2013, PLoS Comput. Biol..

[28]  Phil McMinn,et al.  Modelling complex biological systems using an agent-based approach. , 2012, Integrative biology : quantitative biosciences from nano to macro.

[29]  Keijo Ruohonen,et al.  Developing Itô stochastic differential equation models for neuronal signal transduction pathways , 2006, Comput. Biol. Chem..

[30]  J. Stelling,et al.  Genome‐scale metabolic networks , 2009, Wiley interdisciplinary reviews. Systems biology and medicine.

[31]  Xiaoxiao Wan,et al.  FoxP3+RORγt+ T Helper Intermediates Display Suppressive Function against Autoimmune Diabetes , 2010, The Journal of Immunology.

[32]  Idan Segev,et al.  Ion Channel Stochasticity May Be Critical in Determining the Reliability and Precision of Spike Timing , 1998, Neural Computation.

[33]  C DeLisi Mathematical modeling in immunology. , 1983, Annual review of biophysics and bioengineering.

[34]  Yuelei Shen,et al.  TGF-β-induced Foxp3 inhibits TH17 cell differentiation by antagonizing RORγt function , 2008, Nature.