SIMULATION-DRIVEN DYNAMIC CLAMPING OF NEURONS

We describe an experimentation environment that enables large-scale numerical simulations of neural microphysiology to be fed back onto living neurons in-vitro via dynamic wholecell patch clamping - in effect making living neurons and simulated neurons part of the same neural circuit. Owing to high computational demands, the experimental testbed will be dispersed over a local area network comprising several high performance computing resources. Parallel execution, including feedback between the simulation components, will be managed by the Tarragon, a programming model and run time library that supports asynchronous data driven execution. Tarragon's execution model matches the underlying dynamics of Monte Carlo simulation of diffusive processes and it masks the long network latencies entailed in coupled dispersed simulations. We discuss Tarragon and show how its data driven execution model can be used to dynamically feed back the results of a neural circuit simulation onto living cells in order to better understand the underlying signaling pathways between and within living cells.

[1]  David Potter Computational physics , 1973 .

[2]  Ian Foster,et al.  The Grid 2 - Blueprint for a New Computing Infrastructure, Second Edition , 1998, The Grid 2, 2nd Edition.

[3]  Larry Carter,et al.  Bandwidth-centric allocation of independent tasks on heterogeneous platforms , 2002, Proceedings 16th International Parallel and Distributed Processing Symposium.

[4]  Sanjay Ranka,et al.  Parallel Incremental Graph Partitioning , 1997, IEEE Trans. Parallel Distributed Syst..

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

[6]  K. Mani Chandy,et al.  Asynchronous distributed simulation via a sequence of parallel computations , 1981, CACM.

[7]  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.

[8]  Scott B. Baden,et al.  Performance Tradeoffs in Multi-tier Formulation of a Finite Difference Method , 2001, International Conference on Computational Science.

[9]  Eric L. Schwartz,et al.  Computational Neuroscience , 1993, Neuromethods.

[10]  Steven G. Parker,et al.  Large-scale Computational Science Applications using the SCIRun Problem Solving Environment , 2000 .

[11]  Bradley C. Kuszmaul,et al.  Cilk: an efficient multithreaded runtime system , 1995, PPOPP '95.

[12]  Bradley C. Kuszmaul,et al.  Cilk: an efficient multithreaded runtime system , 1995, PPOPP '95.

[13]  S.R. Kohn,et al.  A Parallel Software Infrastructure for Structured Adaptive Mesh Methods , 1995, Proceedings of the IEEE/ACM SC95 Conference.

[14]  Arvind,et al.  Executing a Program on the MIT Tagged-Token Dataflow Architecture , 1990, IEEE Trans. Computers.

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

[16]  Andrew S. Grimshaw,et al.  Portable run-time support for dynamic object-oriented parallel processing , 1996, TOCS.

[17]  Evangelos P. Markatos,et al.  Load Balancing vs. Locality Management in Shared-Memory Multiprocessors , 1992, ICPP.

[18]  R. Fujimoto Parallel and distributed simulation , 1995, Winter Simulation Conference Proceedings, 1995..

[19]  Francine Berman,et al.  High-performance schedulers , 1998 .

[20]  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.

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

[22]  Scott B. Baden,et al.  A Programming Methodology for Dual-Tier Multicomputers , 2000, IEEE Trans. Software Eng..

[23]  Scott B. Baden,et al.  A large scale Monte Carlo simulator for cellular microphysiology , 2004, 18th International Parallel and Distributed Processing Symposium, 2004. Proceedings..

[24]  Ami Marowka,et al.  The GRID: Blueprint for a New Computing Infrastructure , 2000, Scalable Comput. Pract. Exp..

[25]  W. E Nagel 1988 International conference on supercomputing , 1988 .

[26]  Robert G. Babb,et al.  Parallel Processing with Large-Grain Data Flow Techniques , 1984, Computer.

[27]  Scott B. Baden,et al.  THEMIS: Component Dependence Metadata in Adaptive Parallel Applications , 2001, Parallel Process. Lett..

[28]  Keshav Pingali,et al.  Date movement and control substrate for parallel adaptive applications , 2002, Concurr. Comput. Pract. Exp..

[29]  Scott B. Baden,et al.  Programming Abstractions for Dynamically Partitioning and Coordinating Localized Scientific Calculations Running on Multiprocessors , 1991, SIAM J. Sci. Comput..

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

[31]  Anoop Gupta,et al.  Data locality and load balancing in COOL , 1993, PPOPP '93.

[32]  Scott B. Baden,et al.  Communication overlap in multi-tier parallel algorithms , 1998, Proceedings of the IEEE/ACM SC98 Conference.

[33]  Allen D. Malony,et al.  SMARTS: exploiting temporal locality and parallelism through vertical execution , 1999, ICS '99.

[34]  Hironori Kasahara,et al.  A Data-Localization Compilation Scheme Using Partial-Static Task Assignment for Fortran Coarse-Grain Parallel Processing , 1998, Parallel Comput..

[35]  Philip Heidelberger,et al.  Message proxies for efficient, protected communication on SMP clusters , 1997, Proceedings Third International Symposium on High-Performance Computer Architecture.

[36]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

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

[38]  Scott B. Baden,et al.  Dynamic Partitioning of Non-Uniform Structured Workloads with Spacefilling Curves , 1996, IEEE Trans. Parallel Distributed Syst..

[39]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

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

[41]  Ian Watson,et al.  The Manchester prototype dataflow computer , 1985, CACM.

[42]  Zhiling Lan,et al.  Dynamic Load Balancing of SAMR Applications on Distributed Systems , 2001, ACM/IEEE SC 2001 Conference (SC'01).

[43]  R.M. Fujimoto,et al.  Parallel and distributed simulation systems , 2001, Proceeding of the 2001 Winter Simulation Conference (Cat. No.01CH37304).

[44]  Boris D. Lubachevsky,et al.  Efficient Parallel Simulations of Asynchronous Cellular Arrays , 2005, Complex Syst..

[45]  Jack B. Dennis,et al.  Data Flow Supercomputers , 1980, Computer.

[46]  Nicholas T. Carnevale,et al.  The NEURON Simulation Environment , 1997, Neural Computation.

[47]  Joseph E. Flaherty,et al.  Parallel Partitioning Strategies for the Adaptive Solution of Conservation Laws , 1995 .