Computational modeling of three-dimensional electrodiffusion in biological systems: application to the node of Ranvier.

A computational model is presented for the simulation of three-dimensional electrodiffusion of ions. Finite volume techniques were used to solve the Poisson-Nernst-Planck equation, and a dual Delaunay-Voronoi mesh was constructed to evaluate fluxes of ions, as well as resulting electric potentials. The algorithm has been validated and applied to a generalized node of Ranvier, where numerical results for computed action potentials agree well with cable model predictions for large clusters of voltage-gated ion channels. At smaller channel clusters, however, the three-dimensional electrodiffusion predictions diverge from the cable model predictions and show a broadening of the action potential, indicating a significant effect due to each channel's own local electric field. The node of Ranvier complex is an elaborate organization of membrane-bound aqueous compartments, and the model presented here represents what we believe is a significant first step in simulating electrophysiological events with combined realistic structural and physiological data.

[1]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1952, The Journal of physiology.

[2]  A. Huxley,et al.  The action potential in the myelinated nerve fibre of Xenopus laevis as computed on the basis of voltage clamp data , 1964, The Journal of physiology.

[3]  W. Rall.  Time constants and electrotonic length of membrane cylinders and neurons. , 1969, Biophysical journal.

[4]  C. Nicholson.  Electric current flow in excitable cells J. J. B. Jack, D. Noble &R. W. Tsien Clarendon Press, Oxford (1975). 502 pp., £18.00 , 1976, Neuroscience.

[5]  M H Ellisman,et al.  Molecular specializations of the axon membrane at nodes of Ranvier are not dependent upon myelination , 1979, Journal of neurocytology.

[6]  M H Ellisman,et al.  Immunocytochemical localization of sodium channel distributions in the excitable membranes of Electrophorus electricus. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[7]  William H. Press,et al.  Numerical recipes in Fortran 77 : the art of scientificcomputing. , 1992 .

[8]  J P Wikswo,et al.  A model for axonal propagation incorporating both radial and axial ionic transport. , 1993, Biophysical journal.

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

[10]  J. Ruppersberg.  Ion Channels in Excitable Membranes , 1996 .

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

[12]  Agard,et al.  I5M: 3D widefield light microscopy with better than 100 nm axial resolution , 1999, Journal of microscopy.

[13]  M L Hines,et al.  Neuron: A Tool for Neuroscientists , 2001, The Neuroscientist : a review journal bringing neurobiology, neurology and psychiatry.

[14]  T. J. Sejnowski,et al.  An electro-diffusion model for computing membrane potentials and ionic concentrations in branching dendrites, spines and axons , 1989, Biological Cybernetics.

[15]  Terrence J. Sejnowski,et al.  Synthesis of models for excitable membranes, synaptic transmission and neuromodulation using a common kinetic formalism , 1994, Journal of Computational Neuroscience.

[16]  Hilmar Gugel,et al.  Cooperative 4Pi excitation and detection yields sevenfold sharper optical sections in live-cell microscopy. , 2004, Biophysical journal.

[17]  Mark H. Ellisman,et al.  Development of a model for microphysiological simulations , 2005, Neuroinformatics.

[18]  E. K. Kosmidis,et al.  Three fluorescent protein voltage sensors exhibit low plasma membrane expression in mammalian cells , 2007, Journal of Neuroscience Methods.