A large scale Monte Carlo simulator for cellular microphysiology

Summary form only given. Biological structures are extremely complex at the cellular level. The MCell project has been highly successful in simulating the microphysiology of systems of modest size, but many larger problems require too much storage and computation time to be simulated on a single workstation. MCell-K, a new parallel variant of MCell, has been implemented using the KeLP framework and is running on NPACl's Blue Horizon. MCell-K not only produces validated results consistent with the serial version of MCell but does so with unprecedented scalability. We have thus found a level of description and a way to simulate cellular systems that can approach the complexity of nature on its own terms. At the heart of MCell is a 3D random walk that models diffusion using a Monte Carlo method. We discuss two challenging issues that arose in parallelizing the diffusion process - detecting time-step termination efficiently and performing parallel diffusion of particles in a biophysically accurate way. We explore the scalability limits of the present parallel algorithm and discuss ways to improve upon these limits.

[1]  Scott B. Baden,et al.  Efficient Run-Time Support for Irregular Block-Structured Applications , 1998, J. Parallel Distributed Comput..

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

[3]  J A McCammon,et al.  Analysis of synaptic transmission in the neuromuscular junction using a continuum finite element model. , 1998, Biophysical journal.

[4]  Joseph E. Flaherty,et al.  A hierarchical partition model for adaptive finite element computation , 2000 .

[5]  R. Roitzsch,et al.  KASKADE 3.0 - An Object Oriented Adaptive Finite Element Code , 1995 .

[6]  M Hines,et al.  A program for simulation of nerve equations with branching geometries. , 1989, International journal of bio-medical computing.

[7]  Scott B. Baden,et al.  Abstract KeLP , 2001, PPSC.

[8]  Shahid H. Bokhari,et al.  A Partitioning Strategy for Nonuniform Problems on Multiprocessors , 1987, IEEE Transactions on Computers.

[9]  Joel R. Stiles,et al.  Monte Carlo simulation of neuro-transmitter release using MCell, a general simulator of cellular physiological processes , 1998 .

[10]  D. Gillespie.  A rigorous derivation of the chemical master equation , 1992 .

[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]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[13]  David M. Nicol,et al.  Rectilinear Partitioning of Irregular Data Parallel Computations , 1994, J. Parallel Distributed Comput..

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

[15]  Michael L. Hines,et al.  NEURON — A Program for Simulation of Nerve Equations , 1993 .

[16]  Francine Berman,et al.  Distributing MCell Simulations on the Grid , 2001, Int. J. High Perform. Comput. Appl..

[17]  James M. Bower,et al.  The book of GENESIS - exploring realistic neural models with the GEneral NEural SImulation System (2. ed.) , 1994 .

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

[19]  L M Loew,et al.  A general computational framework for modeling cellular structure and function. , 1997, Biophysical journal.

[20]  A. Arkin,et al.  Simulation of prokaryotic genetic circuits. , 1998, Annual review of biophysics and biomolecular structure.

[21]  Laxmikant V. Kalé,et al.  NAMD: Biomolecular Simulation on Thousands of Processors , 2002, ACM/IEEE SC 2002 Conference (SC'02).

[22]  Bodo Erdmann,et al.  An Object-Oriented Adaptive Finite Element Code , 1997 .

[23]  Masaru Tomita,et al.  E-CELL: software environment for whole-cell simulation , 1999, Bioinform..

[24]  J. D. Teresco,et al.  Parallel structures and dynamic load balancing for adaptive finite element computation , 1998 .

[25]  J. Stiles,et al.  The temperature sensitivity of miniature endplate currents is mostly governed by channel gating: evidence from optimized recordings and Monte Carlo simulations. , 1999, Biophysical journal.

[26]  Joel R. Stiles,et al.  Acetylcholinesterase density and turnover number at frog neuromuscular junctions, with modeling of their role in synaptic function , 1994, Neuron.