Search

We describe the application of a popular and widely available electrical circuit simulation program called SPICE to modeling the electrical behavior of neurons with passive membrane properties and arbitrarily complex dendritic trees. Transient responses may be calculated at any location in the cell model following current, voltage or conductance perturbations at any point. A numbering method is described for binary trees which is helpful in transforming complex dendritic structures into a coded list of short cylindrical dendritic segments suitable for input to SPICE. Individual segments are modeled as isopotential compartments comprised of a parallel resistor and capacitor, representing the transmembrane impedance, in series with one or two core resistors. Synaptic current is modeled by a current source controlled by the local membrane potential and an "alpha-shaped" voltage, thus simulating a conductance change in series with a driving potential. Extensively branched test cell circuits were constructed which satisfied the equivalent cylinder constraints (Rall 1959). These model neurons were perturbed by independent current sources and by synaptic currents. Responses calculated by SPICE are compared with analytical results. With appropriately chosen model parameters, extremely accurate transient calculations may be obtained. Details of the SPICE circuit elements are presented, along with illustrative examples sufficient to allow implementation of passive nerve cell models on a number of common computers. Methods for modeling excitable membrane are presented in the companion paper (Bunow et al. 1985).

We present methods for using the general-purpose network analysis program, SPICE, to construct computer models of excitable membrane displaying Hodgkin-Huxley-like kinetics. The four non-linear partial differential equations of Hodgkin and Huxley (H-H; 1952) are implemented using electrical circuit elements. The H-H rate constants, alpha and beta, are approximated by polynomial functions rather than exponential functions, since the former are handled more efficiently by SPICE. The process of developing code to implement the H-H sodium conductance is described in detail. The Appendix contains a complete listing of the code required to simulate an H-H action potential. The behavior of models so constructed is validated by comparison with the space-clamped and propagating action potentials of Hodgkin and Huxley. SPICE models of multiply branched axons were tested and found to behave as predicted by previous numerical solutions for propagation in inhomogeneous axons. New results are presented for two cases. First, a detailed, anatomically based model is constructed of group Ia input to an alpha-motoneuron with an excitable soma, a myelinated axon and passive dendrites. Second, we simulate interactions among clusters of mixed excitable and passive dendritic spines on an idealized neuron. The methods presented in this paper and its companion (Segev et al. 1985) should permit neurobiologists to construct and explore models which simulate much more closely the real morphological and physiological characteristics of nerve cells.