Abstract
Action potentials (APs) are initiated in the proximal axon of most neurons. In myelinated axons, a 50-times higher sodium channel density in the initial segment compared to the soma may account for this phenomenon. However, little is known about sodium channel density and gating in proximal unmyelinated axons. To study the mechanisms underlying AP initiation in unmyelinated hippocampal mossy fibers of adult mice, we recorded sodium currents in axonal and somatic membrane patches. We demonstrate that sodium channel density in the proximal axon is ∼5 times higher than in the soma. Furthermore, sodium channel activation and inactivation are ∼2 times faster. Modeling revealed that the fast activation localized the initiation site to the proximal axon even upon strong synaptic stimulation, while fast inactivation contributed to energy-efficient membrane charging during APs. Thus, sodium channel gating and density in unmyelinated mossy fiber axons appear to be specialized for robust AP initiation and propagation with minimal current flow.
Introduction
Action potential (AP) initiation is an elementary step in neuronal signal processing. The temporal and spatial pattern of AP propagation along the axon and within the dendritic tree critically depends on the location where the AP is triggered. In the mammalian CNS, APs are preferentially initiated in the proximal axon (Stuart et al., 1997). Computational modeling, sodium (Na+) imaging studies and Na+ current recordings from axonal membrane patches suggest that a 40- to 50-times higher density of functionally specialized voltage-gated Na+ (Nav) channels in the axon initial segment of neocortical pyramidal cells might be critically important for this phenomenon (Kole et al., 2008; Hu et al., 2009). In contrast to the large amount of data that has recently been published on Nav channel gating and density in the initial segment of myelinated axons, only very little is known about Nav channels in mammalian nonmyelinated axons.
In myelinated axons, APs propagate mainly passively along internodal sections. Since the specific capacitance of the myelinated membrane is very low, little energy is required to depolarize the membrane during an AP. By contrast, AP propagation in nonmyelinated axons depends largely on continuous ion flow across the axonal membrane. Since the ion gradient between extracellular space and cytoplasm needs to be maintained by the Na+-K+-ATPase, this process consumes considerably more energy than propagation in myelinated axons, and consequently, AP propagation in nonmyelinated axons of the gray matter has been assumed to be responsible for most of the energy consumption in the brain (Attwell and Laughlin, 2001; Attwell and Gibb, 2005). This view has recently been challenged by Na+ and K+ current recordings evoked by AP waveforms in presynaptic hippocampal mossy fiber boutons (MFBs), suggesting that there is much less overlap of Na+ and K+ currents during the AP than previously thought (Alle et al., 2009). However, it is unclear how the kinetics of Nav channel gating contributes to this phenomenon. This question might be of particular relevance at the AP initiation site in the proximal axon, because this is where the highest Nav channel density can be found within the axon (Kress et al., 2008) and therefore, the largest amount of Na+ ions per length is expected to flow at this location.
In the thin nonmyelinated axons of hippocampal granule cells, the so-called mossy fibers, immunohistochemical data (Kress et al., 2008, 2010) and axonal voltage recordings (Schmidt-Hieber et al., 2008) suggest high Nav channel densities. However, no direct functional analysis of these channels has been performed yet. To address this question, we analyzed the functional properties of Nav channels in outside-out patches from somata and axons of hippocampal granule cells. We found that on average, Nav channel density in the proximal axon is ∼5 times higher than in the soma. Most importantly, the functional properties of somatic and axonal Nav channels appear to be different with up to 2-times faster activation and inactivation kinetics in the axon. Compartmental cable modeling revealed that fast activation kinetics was required to provide a robust axonal initiation site with a moderately enhanced density of Nav channels, while fast axonal inactivation kinetics allowed energy-efficient initiation of APs.
Materials and Methods
Slice preparation.
Transverse 350-μm-thick slices were cut from the hippocampus of 2- to 4-month-old C57BL/6 male or female mice using a custom-made vibratome (Geiger et al., 2002). Animals were kept in an oxygenated chamber for 10 min, anesthetized with isoflurane (Forene, Abbott; 4–5% added to the inspiration air flow) and subsequently killed by decapitation, in accordance with national and institutional guidelines. Experiments were approved by the Animal Care Committee Freiburg (registry X-07/08A). For the dissection and the storage of the slices, a solution containing 87 mm NaCl, 25 mm NaHCO3, 2.5 mm KCl, 1.25 mm NaH2PO4, 25 mm glucose, 75 mm sucrose, 7 mm MgCl2, and 0.5 mm CaCl2 (equilibrated with 95% O2/5% CO2) was used. Slices were incubated at 35°C for ∼30 min and subsequently stored at room temperature.
Electrophysiology.
Slices were superfused with a physiological extracellular solution containing 125 mm NaCl, 25 mm NaHCO3, 2.5 mm KCl, 1.25 mm NaH2PO4, 25 mm glucose, 2 mm CaCl2, and 1 mm MgCl2. Mature granule cells and terminal axonal expansions resulting from the slicing procedure (bleb diameter: ∼2–3 μm) were identified using infrared differential interference contrast (IR-DIC) video microscopy (Shu et al., 2006). Whole-cell patch-clamp recordings from axonal blebs were established under IR-DIC control using a pipette solution containing 130 mm CsCl, 4 mm MgCl2, 4 mm Na2ATP, 10 mm Na2-phosphocreatine, 10 mm EGTA, 10 mm HEPES, and 100 μm AlexaFluor 594 (Invitrogen; pH adjusted to 7.2 with CsOH). Subsequently, the fluorescently labeled soma was identified using a CCD camera (EBFT 512; Princeton Instruments). The excitation light source (Polychrome II; T.I.L.L. Photonics) was coupled to the epifluorescent port of the microscope [Axioskop FS2 (Zeiss); 60× water-immersion objective (Olympus)] via a light guide. To minimize bleaching and avoid phototoxicity, the intensity of the excitation light (550 nm) was reduced to 10% of the maximal value provided by the excitation light source. Patch pipettes (6–10 MΩ for somatic outside-out patch recordings, 6–15 MΩ for axonal outside-out patch recordings, 3–4 MΩ for nucleated patch recordings) were pulled from thick-walled borosilicate glass (outer diameter, 2.0 mm; inner diameter, 0.6 mm). For estimation of Na+ current density, pipettes with similar tip resistances (7–9 MΩ) were used for somatic and axonal recordings. Series resistance in whole-cell mode before withdrawal of the pipette was 20–40 MΩ for outside-out patch recordings and 8–15 MΩ for nucleated patch recordings. Currents were measured with an Axopatch 200A amplifier (Molecular Devices), filtered at 10 kHz (internal 4-pole low-pass Bessel filter), and digitized at 50–100 kHz using a CED 1401plus or Power1401 interface (Cambridge Electronic Design). Voltage command pulses were generated using a home-made data acquisition software (FPulse, U. Fröbe, Physiologisches Institut der Universität Freiburg, Freiburg, Germany) running under Igor 5 (WaveMetrics). To isolate Na+ currents pharmacologically, 20 mm tetraethylammonium chloride (TEA), 3 mm 4-aminopyridine (4-AP), and 100 μm CdCl2 were added to the bath solution. Capacitive and leakage currents were subtracted using the pipette capacitance compensation circuit of the amplifier and two P/4 pulse sequences within each protocol. Patches were held at −80 mV throughout the experiment. Chemicals were obtained from Sigma, Merck, or Gerbu. Recordings were made at 23 ± 2°C. At this temperature, the fastest activation and deactivation time constants were ∼20 μs, close to the temporal resolution limit of the recording system.
To estimate the surface area of outside-out patches, capacitive current transients were elicited by test pulses from −60 to −110 mV. Pipettes were coated with Sylgard for this set of experiments (Sylgard 184, Dow Corning) and had resistances of 7–9 MΩ for both somatic and axonal outside-out patches. Transients were first recorded from intact outside-out patches in close proximity to a Sylgard microsphere attached to a glass tool. Then, the pipette tip was pressed against the microsphere, resulting in complete sealing of the tip (Sakmann and Neher, 1995). Care was taken to maintain the immersion depth of the pipette in the bath during this procedure. Average traces of 50–100 sweeps obtained in the two conditions were subtracted and corrected for leak current. The capacitive charge transfer was obtained by integrating over the initial 100 μs of the transient after pulse onset (Engel and Jonas, 2005). Membrane surface area A (∼4 μm2) of outside-out patches was calculated from the capacitance (∼40 fF) assuming a specific membrane capacitance of 1 μF · cm−2. Membrane surface area of a nucleated patch was calculated from its diameter on the IR-DIC image at the end of a recording. The cross-sectional area of the pipette tip was subtracted.
Data analysis.
For kinetic analysis, current traces were left-shifted in time with respect to the pulse protocol by 60 μs to account for the temporal delay introduced by the amplifier circuit, including the low-pass four-pole Bessel filter (fc = 10 kHz) and the analog–digital converter. This delay was determined from the latency between the onset of the test pulse and the time point of half-maximal rise of tail currents during the deactivation protocol (Bischofberger et al., 2002). The time course of Nav channel activation was fitted with an exponential function with delayed onset where I0 is the amplitude, τa the activation time constant, and δ the delay of Nav channel activation. The time course of deactivation was fitted with a single exponential. To account for inactivation at V = −70 mV, a monoexponential function fitted to the late part of the trace was subtracted before analyzing the trace (Oxford, 1981; Engel and Jonas, 2005). Peak activation data were fitted with a Boltzmann function of the form where Vmid is the midpoint potential and k is the slope factor. Steady-state inactivation data were fitted with a Boltzmann function of the form where f(V) refers to Equation 2. To determine peak activation and steady-state inactivation, data were digitally filtered with a Gaussian filter (fc = 1.5 kHz). This filter was not used when analyzing time constants. Gating properties are summarized in supplemental Table 1 (available at www.jneurosci.org as supplemental material).
Mean peak current–voltage (I–V) relations were fitted with the product of a Boltzmann and a linear function where f(V) is the Boltzmann function of Equation 2, s is the slope, and c is the y-intercept of the linear function. Na+ reversal potential ENa was calculated from the root of Equation 4 as ENa = −c/s. The peak Na+ conductance density ḡNa was calculated using a 30 ms test pulse to −10 mV preceded by a 50 ms prepulse to −120 mV to completely remove steady-state inactivation (see Fig. 2C,D). ḡNa was extrapolated by dividing current traces with an exponential function fitted to the decay of the current to account for the temporal overlap of activation and inactivation (supplemental Fig. 1B, available at www.jneurosci.org as supplemental material) (Keynes and Rojas, 1976). ḡNa was then calculated as ḡNa(V) = INa,c/(VA), where INa,c is the corrected Na+ peak current amplitude, V is the test pulse potential, and A is the membrane surface area of outside-out or nucleated patches, as appropriate. The axonal length was estimated from the shortest distance between the tip of the axonal recording pipette and the somatic point of origin of the axon using the IR-DIC images acquired during the experiments (Kole et al., 2007; Schmidt-Hieber et al., 2008). AP amplitudes were measured from the holding potential of −80 mV. Voltage dependence of deactivation and activation time constants (τd and τa) in Figure 3E were fitted with the function where V is voltage and pn are free parameters.
AP latencies were measured as the difference in time at half-maximal amplitudes. Negative values denote APs preceding the somatic AP. The AP initiation site was determined as the site with the minimal AP latency.
Nonparametric Wilcoxon signed-rank or Mann–Whitney rank-sum tests were used to assess statistical significance. Values are given as mean ± SEM. Error bars in the figures also indicate the SEM. Data were analyzed using C-Stimfit (C.S.-H.) and SciPy (http://www.scipy.org).
Modeling.
Differential equations for kinetic Nav channel gating models were solved using the Q-matrix approach of Colquhoun and Hawkes (1977). Simulations of compartmental cable models were performed using NEURON 7 (Carnevale and Hines, 2006). In current-clamp simulations using compartmental cable models (see Figs. 5, 6B,C, 7; supplemental Figs. 5B, 6, 8, available at www.jneurosci.org as supplemental material), voltage dependence of Nav channel gating rates was corrected for artificial shifts arising from Donnan and liquid junction potentials (∼12 mV) and from the outside-out patch configuration (∼10 mV for inactivation rates) (supplemental Fig. 3 and supplemental text, available at www.jneurosci.org as supplemental material). Experimental data and Nav channel gating models (see Figs. 1⇓⇓–4, 6A; supplemental Figs. 1–4, 5A, 7, available at www.jneurosci.org as supplemental material) are presented without correction for these shifts. Details of the modeling procedures are given in the supplemental material (available at www.jneurosci.org).
Results
A moderately higher density of Nav channels in the proximal mossy fiber axon
To directly measure the densities of Nav channels in hippocampal granule cells, we performed voltage-clamp recordings from somatic and axonal outside-out patches excised from granule cells and mossy fiber axons in acute hippocampal brain slices from adult mice. To record from similar membrane areas, two identical pipettes (range of resistances: 7–9 MΩ) were used for recording of axonal and somatic Na+ currents from the same granule cell (Fig. 1A,B, see Materials and Methods). Axonal recordings were obtained from artificial axonal endings at a distance of 5–35 μm from the soma. Using test pulses to −10 mV, peak Na+ current amplitudes were on average ∼5 times larger in the axon than in the soma (Fig. 1B; supplemental Fig. 1A, available at www.jneurosci.org as supplemental material). To estimate the patch surface area, capacitance of outside-out patches was determined from capacitive currents evoked by −50 mV voltage steps (Fig. 1C) (Sakmann and Neher, 1995). Somatic and axonal outside-out patches had very similar capacitances (44 ± 9 fF, n = 5, and 41 ± 5 fF, n = 5, respectively; p > 0.5), corresponding to a patch area of ∼4 μm2. To calculate peak Na+ conductance densities ḡNa from the peak Na+ current amplitudes and the patch area, we divided current traces with an exponential function fitted to the decay of the current to account for the temporal overlap of activation and inactivation (supplemental Fig. 1B, available at www.jneurosci.org as supplemental material) (Keynes and Rojas, 1976). ḡNa was on average 188 ± 22 pS · μm−2 (18.8 ± 2.2 mS · cm−2) in the soma and 941 ± 173 pS · μm−2 in the axon measured at a distance of 5–35 μm from the soma (n = 14 paired recordings from granule cell soma and axon, respectively; p < 0.001). As a second independent estimate of somatic Nav channel density, we performed nucleated-patch recordings (Martina and Jonas, 1997) revealing a conductance density of 203 ± 64 pS · μm−2, similar to the somatic outside-out patch conductance (Fig. 1D) (n = 6; p > 0.9). ḡNa was highest in the most proximal recordings at ∼5 μm distance from the soma (Fig. 1E) (∼1670 pS · μm−2) and decreased toward more distal parts to an asymptotic value of ∼400 pS · μm−2. Fitting a monoexponential function to the data points revealed that the channel density decayed with a length constant of ∼9 μm. A similar hot spot of Nav channels close to the axon hillock was found in hippocampal mossy fibers of juvenile rats using immunohistochemistry (Kress et al., 2008). Thus, peak Nav channel density at a hot spot close to the soma is 9 times the somatic density. However, with increasing distance from the soma, density decreases to an asymptotic value of about twice the somatic density, leading to an average density ratio of ∼5 throughout the proximal 35 μm of the axon. This is substantially less than the density ratio that has been estimated for the axon initial segment of neocortical and CA1 pyramidal cells (∼40–50 times, Kole et al., 2008; Hu et al., 2009; Lorincz and Nusser, 2010).
Axonal Nav channels activate and inactivate at more negative potentials than somatic channels
It has previously been reported that axonal Nav channels of neocortical pyramidal cells activate and inactivate at more negative potentials than somatic channels (Colbert and Pan, 2002; Kole et al., 2008; Hu et al., 2009). To test whether a similar difference can be found in granule cell axons as well, we measured peak activation and steady-state inactivation of axonal and somatic Nav channels (Fig. 2). To increase the current amplitudes recorded from somatic patches, larger pipettes were used for somatic recordings during these experiments. Midpoint potentials (Vmid) of the peak activation curve were −46 mV in the axon and −40 mV in the soma, corresponding to a voltage shift of about −6 mV for axonal Nav channels (Fig. 2A–C). Similarly, the midpoint potential of the steady-state inactivation curve in the axon (red) was shifted by −8 mV relative to the soma (black; Fig. 2D). In contrast to recent findings in layer V pyramidal cells (Hu et al., 2009), we could not detect a gradual decrease of Vmid of activation along the axon (supplemental Fig. 2, available at www.jneurosci.org as supplemental material). The low activation threshold of axonal Nav channels starting at the very beginning of the axon is consistent with recent immunohistochemical data revealing that the Nav 1.6 isoform is the dominant subunit in granule cell axon initial segments, while Nav 1.2 cannot be detected (Kress et al., 2010). In summary, both activation and inactivation occur at more negative potentials in proximal granule cell axons.
Axonal Nav channels activate and inactivate faster than somatic channels
In addition to the voltage dependence, Nav channel function critically depends on gating kinetics. To examine differences in gating kinetics between the proximal axon and the soma of granule cells, we measured activation, deactivation, and inactivation time constants (Fig. 3). Activation time constants τa were determined by fitting a monoexponential function with delayed onset to the rising phase of Na+ currents (Fig. 3A). The delay was necessary because currents activated with a sigmoidal onset that could not be accounted for by the recording circuit (see Materials and Methods). Deactivation time constants τd were measured by fitting monoexponential functions to the decay of tail currents evoked by a brief 300 μs voltage pulse to 0 mV, followed by test pulses to voltages between −110 and −50 mV (Fig. 3B). To account for the temporal overlap of deactivation and inactivation at voltages ≥−70 mV, a monoexponential function was fitted to the late part of the current and subtracted before fitting the monoexponential function to the tail current (Oxford, 1981; Engel and Jonas, 2005). These measurements revealed that in the voltage range above the AP threshold of about −50 mV (Kress et al., 2008), activation was substantially faster in the axon (Fig. 3E). At a membrane potential of −40 mV, the activation time constant was 93 ± 7 μs in the axon and 183 ± 22 μs in the soma (n = 12; p < 0.005) (supplemental Table 1, available at www.jneurosci.org as supplemental material). Additionally, activation delay was significantly shorter in the axon at −20 and 0 mV (Fig. 3D). By contrast, deactivation time constants did not significantly differ for voltages >−100 mV (Fig. 3B,E). Inactivation time constants τi were obtained by fitting a monoexponential function to the decay of Na+ currents at command voltages between −50 and 20 mV (Fig. 3C). As with activation kinetics, inactivation was significantly faster in the axon than in the soma in a voltage range between −50 and −10 mV (Fig. 3F). In contrast to more depolarized membrane potentials, inactivation kinetics at subthreshold potentials was very similar (supplemental Fig. 4, available at www.jneurosci.org as supplemental material). Together, our data indicate that both activation and inactivation kinetics of Nav channel gating differs substantially in a voltage range of −50 to −10 mV corresponding to the early phase of action potentials, including onset and early upstroke. This suggests that functional properties of axonal Nav channels are specialized for efficient action potential initiation.
Nav channel gating and AP initiation can be faithfully reproduced using an eight-state Nav channel gating model
To test the functional consequences of the differences in Nav channel gating between the axon and the soma, we developed kinetic gating models of Nav channel gating (Fig. 4). We first tried to fit a standard Hodgkin–Huxley type gating model (HH model) to the experimental data using three identical gates for the activation and a single one for the inactivation process (Hodgkin and Huxley, 1952). However, the HH model failed to fit τa and τd at −50 mV, which were very close to each other in the experimental data (supplemental Fig. 5, available at www.jneurosci.org as supplemental material), as observed previously for somatic and axonal Nav channels (Oxford, 1981; Baranauskas and Martina, 2006). To represent the Na+ current data more precisely, we developed an eight-state gating model with different transition rates between multiple closed states, inactivated states, and a single open state (Fig. 4A; see supplemental material, available at www.jneurosci.org) (Vandenberg and Bezanilla, 1991; Baranauskas and Martina, 2006). In contrast to HH-type gating schemes, this model may account for the cooperative gating of different voltage-sensing domains of the Nav channel protein (Chanda et al., 2004). The model could well describe all experimental observations, including τa and τd at −50 mV (Fig. 4B–F).
We then inserted these Na+ conductances into a compartmental cable model of a granule cell (Schmidt-Hieber et al., 2007), using an average channel density in the proximal axon of 940 pS · μm−2 according to our measured data (Fig. 1D,E). Similarly, the spatial distribution of Nav channels was identical to the experimentally observed densities with an asymptotic value of 400 pS · μm−2 for the distal axon (Fig. 5A). This compartmental model faithfully reproduced the granule cell AP phenotype, including maximal rate of rise and decay (Fig. 5B,C; supplemental Fig. 6, available at www.jneurosci.org as supplemental material). The AP was initiated at a distance of ∼30 μm from the soma and propagated along the axon with a speed of ∼0.25 m s−1, similar to experimentally determined values (Fig. 5D; supplemental Fig. 6, available at www.jneurosci.org as supplemental material) (Kress et al., 2008; Schmidt-Hieber et al., 2008). Thus, the compartmental model faithfully reproduces experimental observations using low-intensity stimulation (Schmidt-Hieber et al., 2008).
Fast activation kinetics provides a robust initiation site upon weak and strong synaptic stimuli
Next, we tested how the compartmental model would perform when stronger excitatory synaptic stimuli were used to evoke APs. To assess the effects of activation kinetics under these conditions, we used two further axonal gating models with slower activation time constants (Fig. 6A): a model of Nav channel gating in MFBs in the distal granule cell axon (E&J model) (Engel and Jonas, 2005) and a gating model that had previously been fitted to Na+ current recordings from the soma and the axon of layer V pyramidal cells (M&S model) (Mainen et al., 1995; Kole et al., 2008). We then distributed 500 synapses over the dendritic tree of a model cell and measured the AP initiation site while increasing the conductance of the synapses between simulation runs (Fig. 6B). When the experimentally observed Nav channel kinetics and densities (940 pS · μm−2 in the proximal axon) were used, APs were initiated in the proximal axon throughout the tested range of synaptic conductances from 0.1 to 1.5 nS per synapse. By contrast, when the somatic channel kinetics with slower activation time constants was used uniformly throughout the cell, the initiation site swapped into the dendrites at ∼1 nS per synapse (Fig. 6B, continuous green curve).
Similarly, when the M&S model was used to simulate Nav channels of layer 5 pyramidal cells, APs were initiated in the dendrites from ∼1 nS per synapse on (Fig. 6B, continuous orange curve). This was not due to a difference in voltage dependence of gating models because similar results were observed when the voltage dependence of the M&S model was shifted by 10 mV to the left to match the peak activation curve of our eight-state gating model (supplemental Fig. 7, available at www.jneurosci.org as supplemental material). By contrast, when Nav channel density was increased to 2250 pS · μm−2, similar to the values reported for the AIS in layer V pyramidal cells, the M&S model was able to faithfully reproduce the robust axonal AP initiation in mossy fiber axons (Fig. 6B, dashed orange curve). Replacing the original eight-state model with the E&J model throughout the axon led to dendritic AP initiation from ∼0.5 nS per synapse on (Fig. 6B, continuous blue curve), suggesting that Nav channels in the proximal and the distal axon are specialized to serve distinct functions.
We next examined how the low density of Nav channels in the AIS would affect the impact of proximal axonal inhibition on AP output (Fig. 6C). Thirty inhibitory synaptic conductance changes were distributed along the proximal 30 μm of the axon to simulate synaptic contacts between axoaxonic cells and granule cell axons (Buhl et al., 1994a). The conductance of the inhibitory synapses was increased between simulation runs while APs were evoked by near-threshold excitatory synapses with a constant conductance. When the eight-state model with fast axonal gating was used, AP generation was abolished at ∼1 nS per inhibitory synapse (continuous black curve in Fig. 6C). By contrast, the inhibitory conductance had to be increased 1.7-fold when the M&S model with a high axonal Nav channel density was used (2250 pS · μm−2) (dashed orange curve in Fig. 6C).
Together, these results show that fast axonal activation kinetics is required to reproduce experimental findings obtained with low- and high-intensity synaptic stimulation (Schmidt-Hieber et al., 2008) with a relatively low density of Nav channels in the proximal mossy fiber axon. At the same time, the low density of Nav channels allows axoaxonic interneurons to efficiently modulate AP output.
Nav channel gating promotes energy efficiency of AP initiation
It has recently been suggested that AP generation in principal neurons is relatively energy efficient because of a small overlap of Na+ and K+ currents during the time course of the AP (Alle et al., 2009; Carter and Bean, 2009). In mossy fiber boutons, the Na+ charge transfer during an AP has been reported to be 1.3 times the theoretical minimum that is required to depolarize the membrane (Alle et al., 2009). To assess the energy efficiency of AP initiation in the proximal mossy fiber axon, we simulated APs at room temperature in a long cylinder (Fig. 7; see supplemental material, available at www.jneurosci.org). Energy efficiency was quantified as the ratio of total Na+ charge transfer during an AP and the theoretical minimal charge that is required to depolarize the cell membrane to the AP peak amplitude (Carter and Bean, 2009). When the experimentally determined inactivation kinetics was used, energy efficiency was 1.24 (Fig. 7A), similar to values that have recently been reported for the soma and axon of principal neurons (Alle et al., 2009; Carter and Bean, 2009). Simulations performed at physiological temperatures yielded similar results (supplemental Fig. 8, available at www.jneurosci.org as supplemental material). To determine the influence of Nav channel inactivation kinetics on energy efficiency, we varied inactivation rates by multiplying the original rates with a factor fh (Fig. 7B–D). These simulations revealed that the energy efficiency obtained with the experimentally determined inactivation kinetics (Fig. 7D, blue line) was very close to the minimal value of 1.21 (Fig. 7D, red line). Interestingly, energy efficiency worsened when inactivation kinetics was accelerated beyond fh > 1.5 because of fast recovery from inactivation (Fig. 7C), suggesting that Nav channels are operating in an optimized kinetic range. Since we adjusted ḡK in each simulation individually so that the decay and half duration of the simulated AP closely matched experimentally determined values, Kv channel activation kinetics should only have little impact on our energy efficiency estimates (see supplemental Material for details). We tested this prediction by multiplying Kv channel activation rates with a factor of 0.5 (dotted-dashed curve in Fig. 7D) or 1.5 (dashed curve in Fig. 7D). Indeed, energy efficiency varied little for the different Kv channel kinetics. Most importantly, the energy efficiency obtained with the experimentally determined inactivation kinetics (fh = 1) remained very close to the minimal value independent of the Kv channel activation kinetics.
Fast APs may decrease energy efficiency because of incomplete Nav channel inactivation (Carter and Bean, 2009). In MFBs in the distal granule cell axon, APs have a half-duration of <1 ms (Bischofberger et al., 2002) (recorded at room temperature), but energy efficiency is close to the optimum (Alle et al., 2009). When we used our eight-state model from the proximal axon to replicate a fast AP in a compartmental model of an MFB, we failed to reproduce this finding (Fig. 7E). By contrast, when we used the E&J model with ∼50% faster inactivation kinetics, energy efficiency was close to the values that have been found experimentally (Fig. 7F). Hence, inactivation kinetics in the proximal and distal axon appears to be adapted to the slow and rapid AP repolarization time courses, respectively. Our energy efficiency simulations show that the specialized kinetics of Nav channels in the initial segment of nonmyelinated axons allows initiating APs with little ion flow, thereby promoting energy efficiency. Together, low density and optimized gating kinetics of axonal Nav channels in hippocampal mossy fibers support robust AP initiation with minimal current flow.
Discussion
Our experimental data and simulations indicate that a moderately high Nav channel density and fast gating kinetics in the axon support robust and energy-efficient AP initiation in hippocampal mossy fibers. Nonmyelinated CNS axons such as hippocampal mossy fibers do not possess an anatomically defined axon initial segment or a first node of Ranvier. Nevertheless, our data show that a localized spike initiation zone is generated in the proximal part of the axon (∼15–30 μm from soma) by an enhanced Nav channel density (∼2 to 9 times). In addition, we found that gating kinetics is up to twice as fast in the axon as in the granule cell soma, contributing to robust axonal AP initiation.
Comparing our data with Nav channel gating kinetics in MFBs (Engel and Jonas, 2005) suggests that Nav channel function is compartmentalized throughout the length of the mossy fiber axon. While fast activation kinetics in the proximal axon supports robust AP initiation, fast inactivation kinetics in more distal MFBs contributes to efficient Na+ entry during fast APs (Fig. 8). This heterogeneous distribution of functional Nav channel properties raises the question of the underlying subunit composition. It has recently been shown that Nav 1.6 is the predominant α-subunit in the proximal granule cell axon, while Nav 1.2 cannot be detected (Kress et al., 2010). In good agreement with this finding, channels composed of the Nav 1.6 subunit have been reported to exhibit more hyperpolarized peak activation and steady-state inactivation curves than Nav 1.2 channels (Rush et al., 2005). However, Nav 1.6 expression alone cannot account for the faster activation and inactivation kinetics in the axon (Rush et al., 2005), suggesting that auxiliary β-subunits might contribute to the kinetic differences. β1–3 subunits have been reported to accelerate the inactivation time course when coexpressed with α-subunits (Isom et al., 1992; Smith et al., 1998; Morgan et al., 2000). Of these, β2 and β3 are strongly expressed in the dentate gyrus (Morgan et al., 2000; Yu et al., 2003). Therefore, it seems likely that heteromers of Nav 1.6 α-subunits with a varying proportion of β2/3 auxiliary subunits could explain the compartmentalized function of Nav channels along the mossy fiber axon. Since these axons have very thin diameters of <500 nm, electron microscopy along with immunohistochemistry may be necessary to address this question.
Our results suggest that the mechanisms governing action potential initiation in thin nonmyelinated mossy fiber axons are substantially different from those in myelinated axons of cortical pyramidal cells. The peak Nav channel density in the axon appears to be smaller in granule cells (∼1670 pS · μm−2 vs 2500 pS · μm−2 in pyramidal neurons) (Kole et al., 2008; Hu et al., 2009), and the spatial extent of this enhanced density differs markedly. In pyramidal cells, the Nav channel density appears to be high throughout the AIS until ∼50 μm from the axon hillock. By contrast, in mossy fibers, the density decays within 10–35 μm from soma to values of ∼400 pS · μm−2, leading to an average density of ∼1000 pS · μm−2 within the proximal 35 μm of the axon. Furthermore, the Nav channel densities measured in outside-out patches excised from the proximal axon are perfectly consistent with the properties of the AP waveform measured in whole-cell configuration (supplemental Fig. 6, available at www.jneurosci.org as supplemental material), indicating that axonal Na+ currents can be reliably measured in outside-out patches (Hu et al., 2009).
On the other hand, somatic Nav channel density in granule cells appears to be higher than in layer V pyramidal cells (∼200 pS · μm−2 vs ∼80 pS · μm−2) (Hu et al., 2009). One reason for this might be that the axon initial segment of layer V pyramidal cells has a substantially larger diameter (∼1 μm beyond 40 μm distance from the soma) (Kole et al., 2008) than the proximal mossy fiber axon (0.3–0.4 μm beyond ∼30 μm distance from the soma). Comparing compartmental models of a layer V pyramidal cell (Kole et al., 2008) and a granule cell (Schmidt-Hieber et al., 2007), we estimate that the surface area of the pyramidal cell axon initial segment is ∼6 times larger than that of the proximal granule cell axon within 40 μm distance from the soma. In combination with the higher Nav channel density, axial current flow from the axon into the soma during the AP onset is therefore expected to be substantially larger in pyramidal cells than in granule cells. It seems plausible that the higher Nav channel density in the granule cell soma compensates for the smaller axial current flow provided by the axon.
The higher Nav channel density in pyramidal cell axons might be related to the fact that the APs are not only initiated in the proximal part of the axon, but also have to be reliably propagated toward the first node of Ranvier via axial current flow. By contrast, in nonmyelinated axons of granule cells, axonal Nav channels in neighboring compartments actively contribute to continuous axonal AP propagation. Interestingly, the asymptotic value for axonal Nav channel density in our study (∼400 pS · μm−2) (Fig. 1E) is very close to what was reported for Nav channels in mossy fiber boutons (490 pS · μm−2) (Engel and Jonas, 2005), indicating that the channel density might be homogeneously distributed along the distal mossy fiber axon. The relatively low Nav channel density in the proximal granule cell axon is also in good agreement with a lower PanNav immunostaining density as compared to CA3 pyramidal cell axons (Kress et al., 2010). When compared to a model of Nav channel gating in pyramidal cell axons, activation kinetics of axonal Nav channels is substantially faster in granule cells than in pyramidal cells (Fig. 6A). The faster gating appears to compensate for the lower density of Nav channels at the AP initiation site. As a consequence, AP threshold is always reached first in the axon, even during high background synaptic activity (Fig. 6B).
It has recently been suggested that Nav and Kv channel gating in MFBs is optimized for energy-efficient charging of the axonal membrane (Alle et al., 2009). Here we show that inactivation kinetics of Nav channels in the proximal mossy fiber axon is tuned for efficient depolarization of the membrane during AP initiation, resulting in efficient current flow close to the theoretical minimum. It remains to be determined whether these findings generalize to unmyelinated proximal axons of other cell types. The proximal part of cortical myelinated axons and nonmyelinated mossy fibers is preferentially targeted by axoaxonic cells (Han et al., 1993; Buhl et al., 1994a,b). If APs are initiated with minimal charge, this population of GABAergic interneurons can powerfully modulate the timing and rate of spike output by relatively small conductance changes.
The results show that a moderately higher density of Nav channels with specialized gating properties provide a reliable and unique spike initiation zone in the proximal mossy fiber axon in the absence of a morphologically defined axon initial segment. As a consequence, synaptic inputs from thousands of different dendritic locations are propagated toward the soma and evaluated in a remarkably energy-efficient manner at a single point in the axon close to the axon hillock.
Note added in proof.
Recently, another study investigated Na+ influx in axons of neocortical pyramidal cells using Na+-sensitive fluorescent imaging (Fleidervish et al., 2010). The authors find that Nav channel density in initial segments of these myelinated axons might be ∼3 times the somatic density, similar to our estimates for unmyelinated axons.
Footnotes
C.S.-H. is a Feodor Lynen scholar of the Alexander von Humboldt Foundation. This work was supported by grants from the Deutsche Forschungsgemeinschaft (SFB 780/B3 and Bi 642/2) and the Gatsby Charitable Foundation. We thank Peter Jonas for his continuous and generous support, helpful suggestions, and comments on this manuscript. We thank Henrik Alle, Boris Gutkin, Marco Martina, and Greg Stuart for their helpful comments on previous versions of this manuscript. We thank S. Becherer, M. Northemann, and K. Winterhalter for technical assistance.
- Correspondence should be addressed to Christoph Schmidt-Hieber, Wolfson Institute for Biomedical Research, University College London, London, WC1E 6BT, UK. c.schmidt-hieber{at}ucl.ac.uk