Volume 16 Supplement 1

24th Annual Computational Neuroscience Meeting: CNS*2015

Open Access

A Newton-based shooting method to find synaptic threshold in active cables

BMC Neuroscience201516(Suppl 1):P116


Published: 18 December 2015

The integration of synaptic inputs in a neuron can be nonlinear not just at the axon, but also locally in the dendrites if they are imbued with active voltage-gated ion channels. For example, CA1 pyramidal neurons have high densities of sodium and potassium and currents in their dendrites, and these densities can vary substantially in the arbor [1, 2]. Such nonlinearities can lead to compartmentalized responses to inputs [3], with branches acting as individual nonlinear units in which dendritic spikes occur. A cell can thus function as a multi-layered network with the soma as final output. This motivates determining when a given set of synaptic inputs is large enough to generate a local dendritic spike, or, alternatively, determining the synaptic conductance value(s) at threshold for producing a spike.

Above- and below-threshold conditions are known to be separated by the threshold or critical surface [4]. Here it is shown that the synaptic conductance leading to a threshold solution can be found by modifying Newton methods developed to find steady-state solutions in fluid mechanics [5]. Consider a general form of the cable equation,
u t = L ( x ) u + N ( u , x ) + g s G s ( t , x ) ( u rev - u )

Here u represents the voltage and any gating variables. The first term on the right models the diffusive part of the cable equation, and the next term the nonlinearities from any active voltage-gated ion channels. The last term represents synaptic conductances at points in the dendritic tree with overall strength g s . An unstable threshold solution and accompanying synaptic strength can be found using 1) a preconditioned version of the steady-state cable equation combined with constraints requiring the difference between a shooting solution from rest, u(T,x;g s ), and the critical surface to be perpendicular to the single unstable eigenvector associated with the critical surface. The overall procedure finds the value of g s leading to a solution that asymptotes to the critical surface as t goes to ∞.

Example results from the method are shown in Figure 1. These are for the Fitzhugh-Nagumo model without recovery in a y-branched cable morphology where the daughter branches have half the diameter of the main dendrite. One observes the expected drop of threshold conductance as the synapse moves from the main to daughter branch, but this transition is not monotonic as for a passive cable because of the nonlinear active conductance in the model. Discovering phenomena such as this would be quite laborious without employing a Newton-based method.
Figure 1

Threshold Fitzhugh-Nagumo voltage profile in a branched cable for synapse position just past the branch point (left; blue is main and one daughter branch; red is the other) and critical conductance value as a function of synapse position in main and one daughter branch (right). The branch point is at position × = 40.



WLK was supported by NIH grant 1R01NS077601 & the Janelia Research Campus Visiting Scientist Program

Authors’ Affiliations

Departments of Applied Mathematics and Neurobiology, Northwestern University


  1. Spruston N: Pyramidal neurons: dendritic structure and synaptic integration. Nat. Rev. Neurosci. 2008, 9 (3): 206-21.PubMedView ArticleGoogle Scholar
  2. Johnston D, Magee JC, Colbert CM, Cristie BR: Active properties of neuronal dendrites. Annual Rev. Neurosci. 1996, 19: 165-186.View ArticleGoogle Scholar
  3. Polsky A, Mel BW, Schiller J: Computational subunits in thin dendrites of pyramidal cells. Nat. Neurosci. 2004, 7 (6): 621-7.PubMedView ArticleGoogle Scholar
  4. Mckean HP, Moll V: Stabilization to the standing wave in a simple caricature of the nerve equation. Comm. Pure Appl. Math. 1986, 39 (4): 485-529.View ArticleGoogle Scholar
  5. Tuckerman LS, Huepe C, Brachet ME: Numerical methods for bifurcation problems. Instabilities and non-equilibrium structures IX. Edited by: O. Descalzi, J. Martinez, and S. Rica. 2004, 9:Google Scholar


© Kath 2015

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.