Efficient modeling, simulation and coarse-graining of biological complexity with NFsim

Managing the overwhelming numbers of molecular states and interactions is a fundamental obstacle to building predictive models of biological systems. Here we introduce the Network-Free Stochastic Simulator (NFsim), a general-purpose modeling platform that overcomes the combinatorial nature of molecular interactions. Unlike standard simulators that represent molecular species as variables in equations, NFsim uses a biologically intuitive representation: objects with binding and modification sites acted on by reaction rules. During simulations, rules operate directly on molecular objects to produce exact stochastic results with performance that scales independently of the reaction network size. Reaction rates can be defined as arbitrary functions of molecular states to provide powerful coarse-graining capabilities, for example to merge Boolean and kinetic representations of biological networks. NFsim enables researchers to simulate many biological systems that were previously inaccessible to general-purpose software, as we illustrate with models of immune system signaling, microbial signaling, cytoskeletal assembly and oscillating gene expression.

[1]  C. Rao,et al.  Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm , 2003 .

[2]  David P. Nickerson,et al.  An overview of the CellML API and its implementation , 2010, BMC Bioinformatics.

[3]  B. Palsson,et al.  Towards genome-scale signalling-network reconstructions , 2010, Nature Reviews Genetics.

[4]  J. Tyson,et al.  Computational Cell Biology , 2010 .

[5]  Dietmar Wendt,et al.  The Markoff Automaton: A New Algorithm For Simulating The Time-Evolution Of Large Stochastic Dynamic Systems , 1995 .

[6]  T D Pollard,et al.  Molecular mechanisms controlling actin filament dynamics in nonmuscle cells. , 2000, Annual review of biophysics and biomolecular structure.

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

[8]  J. Gunawardena,et al.  Unlimited multistability in multisite phosphorylation systems , 2009, Nature.

[9]  S. Leibler,et al.  An ultrasensitive bacterial motor revealed by monitoring signaling proteins in single cells. , 2000, Science.

[10]  Richard G Posner,et al.  Trivalent antigens for degranulation of mast cells. , 2007, Organic letters.

[11]  William S. Hlavacek,et al.  RuleMonkey: software for stochastic simulation of rule-based models , 2010, BMC Bioinformatics.

[12]  Ellen S. Vitetta,et al.  An allosteric model for heterogeneous receptor complexes : Understanding bacterial chemotaxis responses to multiple stimuli , 2006 .

[13]  Vincent Danos,et al.  Rule-Based Modelling of Cellular Signalling , 2007, CONCUR.

[14]  William S. Hlavacek,et al.  Simulation of large-scale rule-based models , 2009, Bioinform..

[15]  William S. Hlavacek,et al.  Rejection-free kinetic Monte Carlo simulation of multivalent biomolecular interactions , 2008, 0812.4619.

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

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

[18]  J. S. Parkinson,et al.  Bacterial chemoreceptors: high-performance signaling in networked arrays. , 2008, Trends in biochemical sciences.

[19]  Yuichiro Maéda,et al.  Capping protein binding to actin in yeast , 2004, The Journal of cell biology.

[20]  T. Pollard,et al.  Cellular Motility Driven by Assembly and Disassembly of Actin Filaments , 2003, Cell.

[21]  Kelly R. Spiece,et al.  Redefining the “Scientific Method” , 2000 .

[22]  James R Faeder,et al.  Rule-based modeling of biochemical systems with BioNetGen. , 2009, Methods in molecular biology.

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

[24]  Xiongwu Wu,et al.  Chemoreceptors in Caulobacter crescentus: Trimers of Receptor Dimers in a Partially Ordered Hexagonally Packed Array , 2008, Journal of bacteriology.

[25]  D. Lauffenburger,et al.  Input–output behavior of ErbB signaling pathways as revealed by a mass action model trained against dynamic data , 2009, Molecular systems biology.

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

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

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

[29]  R. B. Jensen,et al.  Location and architecture of the Caulobacter crescentus chemoreceptor array , 2008, Molecular microbiology.

[30]  Vincent Danos,et al.  Scalable Simulation of Cellular Signaling Networks , 2007, APLAS.

[31]  Suresh Kumar Poovathingal,et al.  Global parameter estimation methods for stochastic biochemical systems , 2010, BMC Bioinformatics.

[32]  Thomas D Pollard,et al.  Pathway of Actin Filament Branch Formation by Arp2/3 Complex* , 2008, Journal of Biological Chemistry.

[33]  William S. Hlavacek,et al.  Rule-based modeling of biochemical networks , 2005, Complex..

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

[35]  Michael Hucka,et al.  A Correction to the Review Titled "Rules for Modeling Signal-Transduction Systems" by W. S. Hlavacek et al. , 2006, Science's STKE.

[36]  G. L. Hazelbauer,et al.  Cellular Stoichiometry of the Components of the Chemotaxis Signaling Complex , 2004, Journal of bacteriology.

[37]  Ned S. Wingreen,et al.  Chemotaxis in Escherichia coli: A Molecular Model for Robust Precise Adaptation , 2007, PLoS Comput. Biol..

[38]  T D Pollard,et al.  Rate constants for the reactions of ATP- and ADP-actin with the ends of actin filaments , 1986, The Journal of cell biology.

[39]  W. S. Hlavacek,et al.  Mathematical and computational models of immune-receptor signalling , 2004, Nature Reviews Immunology.

[40]  T D Pollard,et al.  Direct real-time observation of actin filament branching mediated by Arp2/3 complex using total internal reflection fluorescence microscopy , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[41]  J. Tyson,et al.  Design principles of biochemical oscillators , 2008, Nature Reviews Molecular Cell Biology.

[42]  W. S. Hlavacek,et al.  Modeling multivalent ligand-receptor interactions with steric constraints on configurations of cell-surface receptor aggregates. , 2010, Biophysical journal.

[43]  T. Henzinger,et al.  Executable cell biology , 2007, Nature Biotechnology.

[44]  James R Faeder,et al.  Kinetic Monte Carlo method for rule-based modeling of biochemical networks. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[45]  Laurent Blanchoin,et al.  Stochastic severing of actin filaments by actin depolymerizing factor/cofilin controls the emergence of a steady dynamical regime. , 2008, Biophysical journal.

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

[47]  Yuhai Tu,et al.  How white noise generates power-law switching in bacterial flagellar motors. , 2005, Physical review letters.

[48]  D. Bray,et al.  A free-energy-based stochastic simulation of the Tar receptor complex. , 1999, Journal of molecular biology.

[49]  Tuomas P. J. Knowles,et al.  An Analytical Solution to the Kinetics of Breakable Filament Assembly , 2009, Science.

[50]  Jonathan A. Cooper,et al.  Dynamics of capping protein and actin assembly in vitro: uncapping barbed ends by polyphosphoinositides , 1996, The Journal of cell biology.

[51]  A S Perelson,et al.  Equilibrium theory for the clustering of bivalent cell surface receptors by trivalent ligands. Application to histamine release from basophils. , 1984, Biophysical journal.

[52]  W. S. Hlavacek,et al.  Investigation of Early Events in FcεRI-Mediated Signaling Using a Detailed Mathematical Model1 , 2003, The Journal of Immunology.

[53]  Ned S. Wingreen,et al.  Precise adaptation in bacterial chemotaxis through “assistance neighborhoods” , 2006, Proceedings of the National Academy of Sciences.

[54]  Jeff Hasty,et al.  Delay-induced degrade-and-fire oscillations in small genetic circuits. , 2009, Physical review letters.

[55]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..