Facile: a command-line network compiler for systems biology

BackgroundA goal of systems biology is the quantitative modelling of biochemical networks. Yet for many biochemical systems, parameter values and even the existence of interactions between some chemical species are unknown. It is therefore important to be able to easily investigate the effects of adding or removing reactions and to easily perform a bifurcation analysis, which shows the qualitative dynamics of a model for a range of parameter values.ResultsWe present Facile, a Perl command-line tool for analysing the dynamics of a systems biology model. Facile implements the law of mass action to automatically compile a biochemical network (written as, for example, E + S <-> C) into scripts for analytical analysis (Mathematica and Maple), for simulation (XPP and Matlab), and for bifurcation analysis (AUTO). Facile automatically identifies mass conservations and generates the reduced form of a model with the minimum number of independent variables. This form is essential for bifurcation analysis, and Facile produces a C version of the reduced model for AUTO.ConclusionFacile is a simple, yet powerful, tool that greatly accelerates analysis of the dynamics of a biochemical network. By acting at the command-line and because of its intuitive, text-based input, Facile is quick to learn and can be incorporated into larger programs or into automated tasks.

[1]  Clifford A. Shaffer,et al.  The JigCell Model Builder and Run Manager , 2004, Bioinform..

[2]  Nicolas Le Novère,et al.  SBMLeditor: effective creation of models in the Systems Biology Markup Language (SBML) , 2007, BMC Bioinformatics.

[3]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[4]  Bruce E. Shapiro,et al.  MathSBML: a package for manipulating SBML-based biological models , 2004, Bioinform..

[5]  김삼묘,et al.  “Bioinformatics” 특집을 내면서 , 2000 .

[6]  H. Sauro,et al.  Conservation analysis in biochemical networks: computational issues for software writers. , 2004, Biophysical chemistry.

[7]  Roger Brent,et al.  Automatic generation of cellular reaction networks with Moleculizer 1.0 , 2005, Nature Biotechnology.

[8]  Rui Alves,et al.  Tools for kinetic modeling of biochemical networks , 2006, Nature Biotechnology.

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

[10]  Hiroaki Kitano,et al.  Biological robustness , 2008, Nature Reviews Genetics.

[11]  Leslie M Loew,et al.  The Virtual Cell project. , 2002, Novartis Foundation symposium.

[12]  T. Elston,et al.  Stochasticity in gene expression: from theories to phenotypes , 2005, Nature Reviews Genetics.

[13]  Hamid Bolouri,et al.  Dizzy: Stochastic Simulation of Large-scale Genetic Regulatory Networks , 2005, J. Bioinform. Comput. Biol..

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

[15]  Michael Hucka,et al.  SBMLToolbox: an SBML toolbox for MATLAB users , 2006, Bioinform..

[16]  Bruce E. Shapiro,et al.  Cellerator: extending a computer algebra system to include biochemical arrows for signal transduction simulations , 2003, Bioinform..

[17]  Hao Zhu,et al.  Cellware-a multi-algorithmic software for computational systems biology , 2004, Bioinform..

[18]  Edda Klipp,et al.  Systems Biology , 1994 .