A three-dimensional mathematical model for the signal propagation on a neuron's membrane

In order to be able to examine the extracellular potential's influence on network activity and to better understand dipole properties of the extracellular potential, we present and analyze a three-dimensional formulation of the cable equation which facilitates numeric simulations. When the neuron's intra- and extracellular space is assumed to be purely resistive (i.e., no free charges), the balance law of electric fluxes leads to the Laplace equation for the distribution of the intra- and extracellular potential. Moreover, the flux across the neuron's membrane is continuous. This observation already delivers the three dimensional cable equation. The coupling of the intra- and extracellular potential over the membrane is not trivial. Here, we present a continuous extension of the extracellular potential to the intracellular space and combine the resulting equation with the intracellular problem. This approach makes the system numerically accessible. On the basis of the assumed pure resistive intra- and extracellular spaces, we conclude that a cell's out-flux balances out completely. As a consequence neurons do not own any current monopoles. We present a rigorous analysis with spherical harmonics for the extracellular potential by approximating the neuron's geometry to a sphere. Furthermore, we show with first numeric simulations on idealized circumstances that the extracellular potential can have a decisive effect on network activity through ephaptic interactions.

[1]  Klas H. Pettersen,et al.  Modeling the Spatial Reach of the LFP , 2011, Neuron.

[2]  C. Koch,et al.  On the origin of the extracellular action potential waveform: A modeling study. , 2006, Journal of neurophysiology.

[3]  Gabriel Wittum,et al.  UG 4: A novel flexible software system for simulating PDE based models on high performance computers , 2013, Comput. Vis. Sci..

[4]  G. Holt A critical reexamination of some assumptions and implications of cable theory in neurobiology , 1998 .

[5]  A. Hodgkin,et al.  Propagation of electrical signals along giant nerve fibres , 1952, Proceedings of the Royal Society of London. Series B - Biological Sciences.

[6]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[7]  Wilfrid Rall,et al.  Theoretical significance of dendritic trees for neuronal input-output relations , 1964 .

[8]  Gabriel Wittum,et al.  Modeling and simulation for three-dimensional signal propagation in passive dendrites , 2006 .

[9]  Alwyn C. Scott,et al.  The electrophysics of a nerve fiber , 1975 .

[10]  W. Rall Theory of Physiological Properties of Dendrites , 1962, Annals of the New York Academy of Sciences.

[11]  James P. Evans The Origin , 2009, Genetics in Medicine.

[12]  Christof Koch,et al.  Ephaptic coupling of cortical neurons , 2011, Nature Neuroscience.

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

[14]  Gabriel Wittum,et al.  A three-dimensional mathematical model of active signal processing in axons , 2010, Comput. Vis. Sci..

[15]  Christof Koch,et al.  Electrical Interactions via the Extracellular Potential Near Cell Bodies , 1999, Journal of Computational Neuroscience.

[16]  O Herreras,et al.  Activity-dependent changes of tissue resistivity in the CA1 region in vivo are layer-specific: modulation of evoked potentials , 2001, Neuroscience.

[17]  S. Moran,et al.  On the Origin , 2012 .

[18]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .

[19]  C. Koch,et al.  The origin of extracellular fields and currents — EEG, ECoG, LFP and spikes , 2012, Nature Reviews Neuroscience.

[20]  Hartmut Bossel,et al.  Modeling and simulation , 1994 .