Email updates

Keep up to date with the latest news and content from BMC Neuroscience and BioMed Central.

This article is part of the supplement: Seventeenth Annual Computational Neuroscience Meeting: CNS*2008

Open Access Poster presentation

Reaction-diffusion in complex 3D geometries: mesh construction and stochastic simulation with STEPS

Stefan Wils12* and Erik De Schutter12

Author Affiliations

1 Computational Neuroscience Unit, Okinawa Institute of Science and Technology, Onna, Okinawa 904-0411, Japan

2 Theoretical Neurobiology Laboratory, University of Antwerp, Antwerp B-2610, Belgium

For all author emails, please log on.

BMC Neuroscience 2008, 9(Suppl 1):P85  doi:10.1186/1471-2202-9-S1-P85

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2202/9/S1/P85


Published:11 July 2008

© 2008 Wils and De Schutter; licensee BioMed Central Ltd.

Poster presentation

Computational modeling and simulation of signaling pathways has become a valuable and established tool for studying the molecular aspects of biological systems. Aside from the reaction kinetics themselves, the spatial organization and morphology of the cell is known to play an active role in many of these pathways through chemical compartmentalization and by various transport processes. The first step towards simulating these morphological aspects, therefore, is to construct a computational mesh that can adequately represent such features. This is especially challenging for neurons, where functionally relevant details span a wide range of shapes and scales, from the dendritic arborization all the way down to the ultrastructural details surrounding and supporting the synapses.

We show how key features of neural morphology can be described procedurally using various geometric primitives, constructive solid geometry and transformations. Using an adaptation of the elegant algorithm described in [1], this procedural description can then be transformed into tetrahedral meshes with good computational properties. Such meshes can accurately trace complex boundary shapes and are capable of locally adapting their resolution. This adaptation allows a wide range of spatial detail to be represented within a single mesh, while keeping the number of voxels required minimal. In contrast to many third-party packages that are commonly used for tetrahedral mesh generation, this approach is very flexible and can be extended easily to deal with new problems.

We also show how these meshes can be used for simulation, by extending Gillespie's Direct Method (DM) [2] to deal with spatial processes. Our method treats each tetrahedral voxel as a well-mixed volume in which reaction processes can occur. Voxels can be coupled to their neighbouring voxels by diffusive fluxes or to neighbouring boundary triangles to allow simulation of volume-surface interactions, such as ligand-receptor binding and channel currents.

All this work has been implemented in a modular simulator called STEPS. The basic platform of STEPS is provided by Python, which is a general and widely adopted scripting language. Using Python scripts, modelers can describe the kinetics of the reaction-diffusion model; generate and annotate the 3D tetrahedral mesh; simulate and control the system; and process, visualize or export the simulation results. Efficiency is not compromised, since the computational core is still implemented in C/C++. Multiple solvers are implemented and can be used interchangeably; this includes the extension of DM mentioned above, but also a well-mixed solver for simulations in which diffusion and concentration gradients are not required. Our software runs on various platforms, including Un*x, Mac OS X and Windows (via Cygwin).

References

  1. Persson PO, Strang G: A Simple Mesh Generator in MATLAB.

    SIAM Rev 2004, 46(2):329-345. Publisher Full Text OpenURL

  2. Gillespie DT: A general method for numerically simulating the stochastic evolution of coupled chemical reactions.

    J Comp Phys 1976, 22:403-434. Publisher Full Text OpenURL