Spatial Stochastic Modeling with MCell and CellBlender

This chapter provides a brief introduction to the theory and practice of spatial stochastic simulations. It begins with an overview of different methods available for biochemical simulations highlighting their strengths and limitations. Spatial stochastic modeling approaches are indicated when diffusion is relatively slow and spatial inhomogeneities involve relatively small numbers of particles. The popular software package MCell allows particle-based stochastic simulations of biochemical systems in complex three dimensional (3D) geometries, which are important for many cell biology applications. Here, we provide an overview of the simulation algorithms used by MCell and the underlying theory. We then give a tutorial on building and simulating MCell models using the CellBlender graphical user interface, that is built as a plug-in to Blender, a widely-used and freely available software platform for 3D modeling. The tutorial starts with simple models that demonstrate basic MCell functionality and then advances to a number of more complex examples that demonstrate a range of features and provide examples of important biophysical effects that require spatially-resolved stochastic dynamics to capture.

[1]  W. Jäger,et al.  On explosions of solutions to a system of partial differential equations modelling chemotaxis , 1992 .

[2]  Katherine C. Chen,et al.  Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. , 2003, Current opinion in cell biology.

[3]  Alexandra Jilkine,et al.  A Density-Dependent Switch Drives Stochastic Clustering and Polarization of Signaling Molecules , 2011, PLoS Comput. Biol..

[4]  David Fange,et al.  Noise-Induced Min Phenotypes in E. coli , 2006, PLoS Comput. Biol..

[5]  Andreas Hellander,et al.  URDME: a modular framework for stochastic simulation of reaction-transport processes in complex geometries , 2012, BMC Systems Biology.

[6]  D. Bray,et al.  Stochastic simulation of chemical reactions with spatial resolution and single molecule detail , 2004, Physical biology.

[7]  T. Bartol,et al.  Monte Carlo simulation of miniature endplate current generation in the vertebrate neuromuscular junction. , 1991, Biophysical journal.

[8]  Scott B. Baden,et al.  Fast Monte Carlo Simulation Methods for Biological Reaction-Diffusion Systems in Solution and on Surfaces , 2008, SIAM J. Sci. Comput..

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

[10]  Karen Lipkow,et al.  Model for Protein Concentration Gradients in the Cytoplasm , 2008, Cellular and molecular bioengineering.

[11]  T. Bartol,et al.  Miniature endplate current rise times less than 100 microseconds from improved dual recordings can be modeled with passive acetylcholine diffusion from a synaptic vesicle. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[12]  B. Kholodenko Cell-signalling dynamics in time and space , 2006, Nature Reviews Molecular Cell Biology.

[13]  Jason M Haugh,et al.  A unified model for signal transduction reactions in cellular membranes. , 2002, Biophysical journal.

[14]  Linda R Petzold,et al.  Validity conditions for stochastic chemical kinetics in diffusion-limited systems. , 2014, The Journal of chemical physics.

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

[16]  Andrew D Rutenberg,et al.  Pattern formation inside bacteria: fluctuations due to the low copy number of proteins. , 2003, Physical review letters.

[17]  H V Westerhoff,et al.  Why cytoplasmic signalling proteins should be recruited to cell membranes. , 2000, Trends in cell biology.

[18]  T. Bartol,et al.  Monte Carlo Methods for Simulating Realistic Synaptic Microphysiology Using MCell , 2000 .

[19]  Johan Hattne,et al.  Stochastic reaction-diffusion simulation with MesoRD , 2005, Bioinform..

[20]  L. Loew,et al.  The Virtual Cell: a software environment for computational cell biology. , 2001, Trends in biotechnology.

[21]  Olaf Wolkenhauer,et al.  Stochastic approaches in systems biology , 2011, Wiley interdisciplinary reviews. Systems biology and medicine.

[22]  James R Faeder,et al.  Efficient modeling, simulation and coarse-graining of biological complexity with NFsim , 2011, Nature Methods.

[23]  D. Lauffenburger,et al.  Physicochemical modelling of cell signalling pathways , 2006, Nature Cell Biology.

[24]  J. Elf,et al.  Spontaneous separation of bi-stable biochemical systems into spatial domains of opposite phases. , 2004, Systems biology.

[25]  Andreas Hellander,et al.  Perspective: Stochastic algorithms for chemical kinetics. , 2013, The Journal of chemical physics.

[26]  James R Faeder,et al.  Rule‐based modeling: a computational approach for studying biomolecular site dynamics in cell signaling systems , 2014, Wiley interdisciplinary reviews. Systems biology and medicine.

[27]  Frank Noé,et al.  Simulation tools for particle-based reaction-diffusion dynamics in continuous space , 2014, BMC biophysics.

[28]  P. K. Maini,et al.  Overview of Mathematical Approaches Used to Model Bacterial Chemotaxis II: Bacterial Populations , 2008, Bulletin of mathematical biology.

[29]  D A Lauffenburger,et al.  Physical modulation of intracellular signaling processes by locational regulation. , 1997, Biophysical journal.

[30]  Roger Brent,et al.  Detailed Simulations of Cell Biology with Smoldyn 2.1 , 2010, PLoS Comput. Biol..

[31]  S. Schnell,et al.  Modelling reaction kinetics inside cells. , 2008, Essays in biochemistry.