A continuous cable method for determining the transient potential in passive dendritic trees of known geometry

Models using cable equations are increasingly employed in neurophysiological analyses, but the amount of computer time and memory required for their implementation are prohibitively large for many purposes and many laboratories. A mathematical procedure for determining the transient voltage response to injected current or synaptic input in a passive dendritic tree of known geometry is presented that is simple to implement since it is based on one equation. It proved to be highly accurate when results were compared to those obtained analytically for dendritic trees satisfying equivalent cylinder constraints. In this method the passive cable equation is used to express the potential for each interbranch segment of the dendritic tree. After applying boundary conditions at branch points and terminations, a system of equations for the Laplace transform of the potential at the ends of the segments can be readily obtained by inspection of the dendritic tree. Except for the starting equation, all of the equations have a simple format that varies only with the number of branches meeting at a branch point. The system of equations is solved in the Laplace domain, and the result is numerically inverted back to the time domain for each specified time point (the method is independent of any time increment Δt). The potential at any selected location in the dendritic tree can be obtained using this method. Since only one equation is required for each interbranch segment, this procedure uses far fewer equations than comparable compartmental approaches. By using significantly less computer memory and time than other methods to attain similar accuracy, this method permits extensive analyses to be performed rapidly on small computers. One hopes that this will involve more investigators in modeling studies and will facilitate their motivation to undertake realistically complex dendritic models.

[1]  T. Poggio,et al.  Nonlinear interactions in a dendritic tree: localization, timing, and role in information processing. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Christof Koch,et al.  A simple algorithm for solving the cable equation in dendritic trees of arbitrary geometry , 1985, Journal of Neuroscience Methods.

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

[4]  V. Torre,et al.  A Volterra representation for some neuron models , 1977, Biological Cybernetics.

[5]  W. Rall Electrophysiology of a dendritic neuron model. , 1962, Biophysical journal.

[6]  B. Davies,et al.  Numerical Inversion of the Laplace Transform: A Survey and Comparison of Methods , 1979 .

[7]  J. Miller,et al.  Factors that control amplitude of EPSPs in dendritic neurons. , 1983, Journal of neurophysiology.

[8]  J. P. Miller,et al.  Modeling the electrical behavior of anatomically complex neurons using a network analysis program: Passive membrane , 1985, Biological Cybernetics.

[9]  T. Poggio,et al.  A New Approach to Synaptic Interactions , 1978 .

[10]  B Horwitz,et al.  Unequal diameters and their effects on time-varying voltages in branched neurons. , 1983, Biophysical journal.

[11]  J. Jack,et al.  Electric current flow in excitable cells , 1975 .

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

[13]  D A Turner,et al.  Segmental cable evaluation of somatic transients in hippocampal neurons (CA1, CA3, and dentate). , 1984, Biophysical journal.

[14]  J Rinzel,et al.  Transient response in a dendritic neuron model for current injected at one branch. , 1974, Biophysical journal.

[15]  I. Segev,et al.  Modeling the electrical behavior of anatomically complex neurons using a network analysis program: Excitable membrane , 2004, Biological Cybernetics.

[16]  T. Poggio,et al.  Retinal ganglion cells: a functional interpretation of dendritic morphology. , 1982, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[17]  J J Jack,et al.  The propagation of transient potentials in some linear cable structures , 1971, The Journal of physiology.

[18]  J. Cowan,et al.  Transient potentials in dendritic systems of arbitrary geometry. , 1974, Biophysical journal.

[19]  G. V. Narayanan,et al.  Numerical operational methods for time‐dependent linear problems , 1982 .

[20]  D. Perkel,et al.  Electrotonic properties of neurons: steady-state compartmental model. , 1978, Journal of neurophysiology.

[21]  W. Rall Branching dendritic trees and motoneuron membrane resistivity. , 1959, Experimental neurology.

[22]  Andrew H. Sherman,et al.  Algorithm 533: NSPIV, a Fortran subroutine for sparse Gaussian elimination with partial pivoting [F4] , 1978, TOMS.

[23]  D. Perkel,et al.  Quantitative methods for predicting neuronal behavior , 1981, Neuroscience.

[24]  S. Redman The attenuation of passively propagating dendritic potentials in a motoneurone cable model , 1973, The Journal of physiology.

[25]  B Horwitz,et al.  An analytical method for investigating transient potentials in neurons with branching dendritic trees. , 1981, Biophysical journal.