|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1Department of Epileptology, University of Bonn Medical Center, Bonn, Germany; 2Helmholtz-Institute for Radiation and Nuclear Physics and 3Interdisciplinary Center for Complex Systems, University of Bonn, Bonn, Germany; and 4Department of Physiology, Hebrew University-Hadassah School of Medicine, Jerusalem, Israel
Submitted 5 March 2008; accepted in final form 10 July 2008
|
|
ABSTRACT |
|---|
|
|
|
INTRODUCTION |
|---|
|
What factors cause the axon initial segment (AIS) to have the lowest spike threshold? One factor may be a relatively high density of Na+ channels in this region as evidenced in different types of neurons by immunolabelings of Na+ channel proteins (Boiko et al. 2001
, 2003
; Catterall 1981
; Hossain et al. 2005
; Pan et al. 2006
). Indeed the AIS contains a machinery to concentrate certain types of ion channels. Ankyrin G is a key player in this process as it was shown to be both necessary and sufficient to direct different types of Na+ channels (Garrido et al. 2003
; Zhou et al. 1998
) as well as KV7 (KCNQ) K+ channels (Pan et al. 2006
), to the AIS. Although previous electrophysiological studies using cell-attached recordings have proclaimed a uniform transient Na+ current (INaT) density at AIS and soma (Colbert and Johnston 1996
; Colbert and Pan 2002
), a more recent study employing also Na+ imaging has argued that INaT density is in fact higher at the AIS than in the soma (Kole et al. 2008
). The specific biophysical properties of the Na+ channels expressed at the AIS also may play a role in localizing the spike trigger zone to this region (Colbert and Pan 2002
; Naundorf et al. 2006
). In particular, it was found that Na+ channels at the AIS of cortical neurons exhibit a voltage dependence of activation that is shifted by
8 mV in a hyperpolarized direction compared with somatic Na+ channels (Colbert and Pan 2002
). However, the molecular basis for this functional specialization remains unresolved.
At the AIS, NaV1.1, NaV1.2, and NaV1.6 channels have been detected on the protein level (Boiko et al. 2001
, 2003
; Garrido et al. 2003
; Hossain et al. 2005
; Ogiwara et al. 2007
; Van Wart and Matthews 2006
; Van Wart et al. 2007
). The functional role of NaV1.6 subunits in particular have been assessed in number of investigations in mutant mice lacking NaV1.6 channels, for instance in cerebellar and globus pallidus neurons, as well as dorsal root and trigeminal ganglion cells (Levin et al. 2006
; Mercer et al. 2007
; Raman et al. 1997
). The results argue for a role of NaV1.6 subunits in mediating resurgent and persistent Na+ currents in these cells with a resulting effect on repetitive firing behavior.
A striking biophysical peculiarity of NaV1.6 subunits is its hyperpolarized voltage of activation compared with other Na+ channel isoforms. This finding has been obtained in mouse dorsal root ganglion neurons overexpressing a TTX-insensitive variant of NaV1.6, and thus allowing assessment of the properties of these channel isoforms in isolation in a neuronal cell (Rush et al. 2005
; but see Smith et al. 1998
). We therefore hypothesized that a preponderance of NaV1.6 expression at the AIS may contribute to its low spike threshold in addition to affecting repetitive discharge behavior. We explored the role of this channel subunit in firing behavior of CA1 pyramidal neurons using mice lacking functional NaV1.6 subunits (Scn8amed mice) as well as with computational modeling approaches. Our results indicate a critical role for NaV1.6 in setting the low spike threshold at the AIS of CA1 pyramidal neurons.
|
|
METHODS |
|---|
|
Experiments were performed on mice deficient in functional NaV1.6
-subunits bearing the recessive muscle endplate disease (med) mutation in the Scn8a gene. This mutation causes the expression of a truncated nonfunctional form of the protein by altering mRNA splicing due to insertion of a LINE element in exon 2 (Kohrman et al. 1996
). Heterozygous breeding pairs of Scn8amed/wt mice (C3HeB/FeJ-Scn8amed/J; Stock No. 003798) were acquired from Jackson Laboratories (Bar Harbor, ME). Wild-type (Scn8awt) or mutant (Scn8amed) homozygous littermate offspring animals aged 17–21 days were used in all experiments. All animal experiments were conducted in accordance with the guidelines of the Animal Care and Use Committee of the University of Bonn. For all experiments, animals were heart-perfused with 1–3°C cold sucrose-based artificial cerebrospinal fluid (ACSF) containing (in mM) 56 NaCl, 100 sucrose, 2.5 KCl, 1.25 NaH2PO4, 30 NaHCO3, 1 CaCl2, 5 MgCl2, 1 kynurenic acid, and 20 glucose (95% O2-5% CO2) under deep anesthesia with ketamine (100 mg/kg, Pfizer) and xylazine (15 mg/kg, Bayer). After perfusion mice were decapitated, the brain was quickly removed, and 300-, 400-, or 600-µm-thick transverse hippocampal slices were cut with a vibratome (MICROM) for electrophysiological or immunohistochemical studies.
Immunohistochemistry
Freshly cut 600-µm hippocampal slices were placed in a tissue boat, submerged under Tissue-Tec (Sakura) and carefully frozen over liquid nitrogen before being stored at –80°C. From the frozen tissue 12-µm-thick sections were cut with a cryostat (MICROM) and mounted to either DAKO-slides (DAKO) or Superfrost-plus-slides (Menzel) on which they were allowed to rest for 15 min at 20°C. Then the slides were fixed by submerging them for 2 min into a 1:1 mixture of ethanol and acetone (Merck) and left to dry overnight at 20°C. Finally the slides were stored in a –20°C freezer until the staining experiments were conducted.
Slides were thawed for 30 min at 20°C and afterward briefly washed in PBS (Biochrom AG,). To avoid unspecific antibody binding, the slices were incubated for 2 h at 20°C in blocking solution consisting of PBS, Triton-X100 (0.1%), fetal calf serum (10%; PAA Laboratories), and normal goat serum (5%; Vector, Burlingame, CA). All primary antibodies were diluted 1:200 in blocking solution, and the binding reaction was allowed to take place at 4°C for 12–16 h. For double immunolabelings, primary antibodies were applied together. The primary antibodies used were a monoclonal mouse anti-Ankyrin G antibody directed against the spectrin binding domain of Ankyrin G (Zymed, San Francisco, CA), a polyclonal rabbit anti-NaV1.6 directed against amino acids 1042-1061 of the rat NaV1.6 protein (Alomone Labs), a monoclonal mouse anti-PanNaV antibody and a polyclonal rabbit anti-PanNaV antibody, both raised against amino acids 1491–1508 of the rat NaV1.1 protein with the antigen for the polyclonal antibody containing an additional cystein (Noda et al. 1986
), a sequence identical in all mammalian NaV
-subunits (Sigma-Aldrich). It should be noted that the polyclonal antibody also produced a robust immunolabeling of neuronal somata in the hippocampus, which was absent with the monoclonal antibody (cf. Fig. 1, Ab and B). Labeling of AIS, however, was similar with both antibodies. Excessive unbound primary antibodies were washed away three times at 20°C for 5 min with PBS. Subsequently, slices were incubated for 2 h at 20°C in the dark with FITC- and CY3-conjugated secondary antibodies (Dianova). Secondary antibodies were also diluted 1:200 in blocking solution and applied synchronously. Finally the slides were washed again 3 times in PBS for 5 min at 20°C and furnished with cover slips using a 1:1 mixture of Vectashield-Harding and Vectashield-Harding with DAPI cover media (Vector). The slides were then stored light protected at 4°C.
|
2 h prior to the imaging session. All images were taken in one continuous imaging session, where apart from focal plane all laser and microscope settings remained untouched. The pinhole was set to 0.83 Airy units. Detector gain was set to
60%. To determine mean Na+ channel density at AIS, we first defined regions of interest (ROI) corresponding to AIS based on the Ankyrin G staining. The mean staining intensity for both Ankyrin G and PanNaV was measured. From each section values for ten AIS were determined. We calculated the intensity of PanNaV staining as a ratio of the average intensity in the PanNaV channel divided by the corresponding average intensity in the Ankyrin G channel. Storage of slices and preparation of dissociated neurons
For electrophysiological experiments, freshly cut slices were first placed into a storage chamber with room temperature (20°C) sucrose-based ACSF containing (in mM) 60 NaCl, 100 sucrose, 2.5 KCl, 1.25 NaH2PO4, 26 NaHCO3, 1 CaCl2, 5 MgCl2, 1 kynurenic acid, and 20 glucose (95% O2-5% CO2) and gradually warmed to 36°C during 30 min. Subsequently, slices were equilibrated in a chamber with sucrose-free ACSF containing (in mM) 125 NaCl, 3.5 KCl, 1.25 NaH2PO4, 26 NaHCO3, 2 CaCl2, 2 MgCl2, and 15 glucose (95% O2-5% CO2) for
30 min at 20°C. For recordings of identified CA1 neurons in the slice preparation, 300 µm slices were used.
For preparation of dissociated neurons, 400 µm slices were placed in 5 ml of trituration solution containing (in mM) 145 Na-methanesulfonate, 3 KCl, 10 N-2-hydroxy-ethylpiperazine-N'-2-ethane sulfonic acid (HEPES), 0.5 CaCl2, 1 MgCl2, and 15 glucose. Solution pH was adjusted to 7.4 with NaOH. Pronase (protease type XIV; 2 mg/ml; Sigma, St. Louis, MO) was added to the oxygenated buffer (100% O2). After two incubation periods, 10 min at 35°C and followed by 10 min at 20°C, slices were washed with pronase-free buffer saline of identical composition and transferred to a Petri dish containing 5 poly-L-lysine-coated cover slips. The CA1 region was microdissected under a binocular and triturated with fire-polished glass pipettes of decreasing aperture. Cells were allowed to settle for
10 min before removing cover slips and placing them into a submerged chamber mounted on the headstage of an upright microscope (Axioskop F-2, Zeiss). Cells were equilibrated for further 10 min before recording was attempted. Whole cell recordings of dissociated neurons were performed only on pyramidal-shaped neurons with a smooth surface and a three-dimensional contour. All cells recorded possessed a clearly identifiable apical dendrite and remnants of basal dendrites and the axon.
Electrophysiology
Patch pipettes with a resistance of 3–5 M
were pulled from borosilicate glass capillaries (1.5 mm OD, 1 mm ID; Science Products) on a Narishige PP-830 puller (Narishige, Tokyo, Japan) and filled with the appropriate intracellular (IC) solution. Voltage- and current-clamp recordings were conducted at 20 and 30°C, respectively. Data were recorded and stored by a personal computer using a data-acquisition system (Digidata 1322A) and the pClamp9.0 software (Molecular Devices). Unless otherwise indicated data were filtered at 10 kHz and digitized at 100 kHz. Passive membrane properties were quantified as follows. The input resistance was determined in voltage clamp mode according to Ohm's law from the steady-state current response to 5- or 10-mV voltage steps (200 ms) from a –85-mV holding potential and was not significantly different between the mice from both genotypes (Scn8amed 342.52 ± 79.00 M
, Scn8awt 300.60 ± 25.28 M
). Cell capacitance was determined by quantifying the charge (Qc) required to fully charge the membrane. Qc was measured as the total area under the current response to the abovementioned voltage steps, minus the charge flowing across the membrane resistance. Cell capacitance was then calculated as Qc/V, where V is the size of the voltage step (Scn8amed 111.55 ± 15.22 pF, Scn8awt 100.99 ± 8.23 pF; n = 12 and n = 22, respectively).
Electrophysiology
CURRENT-CLAMP RECORDINGS. Current-clamp recordings were performed in intact CA1 neurons in the slice preparation, using a Multiclamp 700B amplifier (Molecular Devices). Whole cell configuration was obtained in voltage-clamp mode before switching to current-clamp mode, where pipette capacitance and bridge balance were monitored and carefully compensated. Cells with native membrane potential more positive than –60 mV were excluded. Subsequently, the slow current-clamp circuit of the amplifier set to 5 s was used to set the initial membrane potential prior to current injection steps to defined values. The intracellular solution used was (in mM) 130 K-gluconate, 20 KCl, 10 HEPES, 0.16 ethylene glycol-bis (2-aminoethylether)-N,N,N',N'-tetraacetic acid (EGTA), 2 Mg-adenosine 5'-triphosphate (ATP), and 2 Na2-ATP; pH was titrated to 7.25 with KOH; osmolality was adjusted to 295 mosM using sucrose. For bath solution, a modified ACSF was used containing (in mM) 124 NaCl, 3.5 KCl, 26 NaHCO3, 1.6 CaCl2, 2 MgCl2, and 10 glucose (95% O2-5% CO2). Temperature for current-clamp recordings was maintained at 30 ± 1°C. The liquid junction potential determined for these solutions was –15 mV, and all values and figures were corrected accordingly.
Analysis of current-clamp recordings
The measured resting membrane potential was not different between Scn8amed (–72.80 ± 1.25 mV) and Scn8awt mice (–73.98 ± 0.66 mV). Spike thresholds were determined by measuring the voltage at which the increase in slope of the voltage trace is maximal. This time point corresponds to the maximum of the second derivation of the voltage step (d2V/dt2) and was determined as the time at which the third derivation of the voltage trace became zero. Spike amplitude was measured as the difference between resting membrane potential and the peak of the spike. The maximal rates of rise and decay were determined as the peak and antipeak of the second derivation of the voltage trace. Spikes during prolonged (600 ms) current injections vary systematically, depending on the time of occurrence during the current injection and the number of prior spikes. We analyzed the first, second and subsequent spikes in an action potential train separately. Analysis of these spike parameters for spikes elicited by 4-ms current injection was done using Clampfit 9.0. Repetitive firing was analyzed using an automated Igor routine that detected spikes and measured their properties.
In addition to these parameters, we determined the axo-somatic delay by assessing the delay between the two peaks observed in the second derivation of the voltage-trace. This assessment was carried out with an automated IGOR detection routine. Each automatically analyzed spike was subsequently inspected. In some cases, the automated IGOR detection routine failed to detect two peaks because of overlap between the two peaks. In these cases, an estimate of the axo-somatic delay had to be obtained by a manual determination.
The size of the spike afterdepolarization (ADP) was determined by measuring the area under the ADP starting from the beginning of the fast afterhyperpolarization to the time when membrane voltage returned to the holding potential. This delivers a value that incorporates both active and passive portions of the ADP. To evaluate the magnitude of the active portion of the ADP, we first estimated the contribution of passive components by obtaining voltage responses to subthreshold current injections of identical duration. These passive voltage responses were scaled so that the peak of the passive response was superimposed to the action potential threshold. The corresponding area approximates the passive response of the neuron, and was subtracted from the total ADP area, yielding the active component of the ADP.
VOLTAGE-CLAMP RECORDINGS. Voltage-clamp recordings of transient Na+ current (INaT) were carried out in dissociated CA1 neurons to obtain a reliable voltage control and to minimize space-clamp problems. Even in dissociated neurons, the large amplitude of INaT necessitated a reduction of the Na+ gradient between bath and intracellular solutions. The following intracellular solution was used (in mM): 110 CsF, 10 HEPES-Na, 11 EGTA, 20 tetraethylammonium-Cl, 2 MgCl2, 0.5 guanosine 5'-triphosphate-tris(hydroxyl-methyl)-aminomethane (GTP-Tris), and 5 ATP-Na2. Osmolality was adjusted with sucrose to 295 mosM; pH to 7.25 with CsOH. The oxygenated bath consisted of (in mM) 30 Na-methanesulphonate, 120 tetraethylammonium-Cl, 10 HEPES, 1.6 CaCl2, 2 MgCl2, 0.2 CdCl2, 5 4-aminopyridine (Acros Organics), and 15 glucose. The pH was adjusted to 7.4 with HCl, osmolality was adjusted to 310 mosM with sucrose, and temperature was maintained at 20 ± 1°C. The liquid junction potential between intra- and extracellular solution was +10 mV.
Recordings of the persistent Na+ current (INaP) were carried out in intact neurons in the slice preparation with intracellular solution containing (in mM) 110 CsF, 10 HEPES-Na, 11 EGTA, 2 MgCl2, 0.5 GTP-Tris, and 2 ATP-Na2. Osmolality was adjusted with mannitol to 295 mosM; pH was adjusted to 7.25 (CsOH). The bath solution consisted of (in mM) 100 Na-methanesulfonate, 40 tetraethylammonium-Cl, 10 HEPES, 2 CaCl2, 3 MgCl2, 0.2 CdCl2, 5 4-aminopyridine, and 15 glucose. pH 7.4, NaOH; osmolality was adjusted to 305 mosM with sucrose. Liquid junction potential was +10.0 mV.
Recordings of the resurgent Na+ current (INaR) were carried out in dissociated neurons with the intracellular solution containing (in mM) 110 CsF, 10 HEPES-Na, 11 EGTA, 2 MgCl2, 0.5 GTP-Tris, and 2 ATP-Na2. Osmolality was adjusted with mannitol to 295 mosM; pH was adjusted to 7.25 using CsOH. The bath solution consisted of (mM) 100 NaCl, 40 tetraethylammonium-Cl, 10 HEPES, 2 CaCl2, 3 MgCl2, 0.2 CdCl2, 5 4-aminopyridine (Acros Organics), and 15 glucose (pH 7.4, NaOH; osmolality was adjusted to 305 mosM with sucrose). Liquid junction potential was –9.99 mV.
Recordings of T-type Ca2+ currents (ICaT) was carried out in slices that had been preincubated for 1 h in 5 ml oxygenated bath containing: omega-CgTx GVIA (2 µM), omega-CgTx MVIIC (3 µM), omega-AgaTx IVA (0.2 µM; Biotrend), and cytochrome C (2 mg/ml) to block N- and P/Q-type Ca2+ channels. Following transfer of the slices to the recording chamber, recordings were carried out with intracellular solution containing (in mM) 105 Cs-methanesulfonate, 25 tetraethylammonium-Cl, 10 HEPES, 5 EGTA, 2 MgCl2, 2 CaCl2, 25 sucrose, 4 ATP-Na2, and 0.3 GTP-Tris; pH was adjusted to 7.2 with CsOH; osmolality with sucrose to 295 mosM. The bath solution contained (in mM) 115 Na-methanesulfonate, 25 tetraethylammonium-Cl, 3.5 KCl, 2 MgCl2, 2 CaCl2, 4 4-aminopyridine, 10 HEPES, 25 glucose, 0.005 tetrodotoxin (Biotrend), and 0.01 nifedipine (pH 7.4, NaOH; osmolality was adjusted to 310 mosM with sucrose). Liquid junction potential was –5.0 mV.
Tight seal whole cell recordings were obtained with a seal resistance >1 G
in all recordings using an Axopatch 200B amplifier (Molecular Devices). Series resistance was 6 ± 2 M
. To improve voltage control, the prediction and compensation dials of the amplifier's series resistance compensation were set between 70 and 90% to achieve a maximal residual voltage error <2 mV (<0.5 mV for recordings of INaP, INaR, and ICaT). All other recordings were excluded. Currents were recorded with the pClamp acquisition and analysis program, sampled at 100 kHz and filtered at 10 kHz (20 and 1 kHz for INaP). All potentials shown were corrected for liquid junction potentials. Recording temperature was 20°C for all voltage-clamp recordings. Unless otherwise indicated, all chemicals or drugs were obtained from Sigma.
Analysis of voltage-clamp recordings
The voltage dependence activation of INaT was determined using standard protocols (see Fig. 3A, inset). Peak currents were fitted using the following Boltzmann function
![]() | (1) |
Peak currents were then converted to conductance G(v) using
![]() | (2) |
The voltage dependence of steady-state inactivation was determined using standard procedures with prepulses (500 ms) to various voltages, followed by a 10-ms test pulse to 0 mV (see Fig. 3C, inset). The peak currents were fitted using
![]() | (3) |
where Imax is the maximal Na+ current, V1/2 is membrane potential at which I(V) is half of Imax and km is the slope at V1/2.
To determine the voltage dependent activation of INaP, the TTX-subtracted current responses to the voltage ramp (Fig. 4A) were converted to conductance using Eq. 2 and subsequently fitted using Eq. 3 (Fig. 4C). In all cases, fitting was done using a Levenberg-Marquardt algorithm.
The magnitude of INaR was determined by analyzing the current responses to different 100-ms test pulses (–100 to –10 mV) following a 15-ms prepulse to 20 mV from a holding potential of 100 mV (Fig. 5A, inset). The amplitude of INaR was determined as the peak current during the test pulse minus the steady-state current at the end of the test pulse (see Fig. 5A).
The amplitude of ICaT was determined by fitting the tail current following a 20-ms depolarization with a biexponential function (see Fig. 6A, inset) using a Levenberg-Marquardt algorithm. Under our recording conditions, the faster deactivating current component represents R-type Ca2+ currents, while the slower component is due to deactivation of T-type Ca2+ currents (Sochivko et al. 2002
). The amplitude corresponding to the slower deactivating component was derived by extrapolation of the fitted curve to the end of the depolarizing voltage step.
All data are presented as averages ± SE. For comparison of means, a two-tailed Student's t-test was performed as appropriate. Differences between axo-somatic spike delay and input-output relations between Scn8amed and Scn8awt mice, were analyzed by MANOVA. For all tests, the significance level was set at P < 0.05. All data analyses were done with the Clampfit 9.0 software (Molecular Devices), Origin 7 (OriginLab, Northampton, MA), IGOR (Wavemetrics, Lake Oswego, OR), SPSS 14.0 (SPSS) and Excel 2003 at a Windows based PC-system (Microsoft, Redmond, WA).
Modeling of a CA1 pyramidal neuron
We have created a model of a CA1 neuron with a realistic morphology and different voltage- and Ca2+-dependent currents with differential subcellular distribution. The modeling environment was Microsoft Windows XP, running on a dual core processor, each Intel Processor with 2.39 GHz, 1.97 GB Ram. The simulation was implemented within the simulation software NEURON (Carnevale and Hines 2006
). The integration time steps were fixed at 0.01 ms. The general approach to model the properties of different ionic currents is based on a Hodgkin-Huxley-type formalism (Hodgkin and Huxley 1952
), where the voltage and time dependence of currents flowing through ion channels is governed by gating particles that determine the opening and closing of the channel pore. The time- and place-dependent total current density im(x, t) through a cell membrane is given by
![]() |
The dynamics of gating particles is governed by the differential equation
![]() |
denotes the equilibrium state, and
p the time constant of the dynamics. In general p
and
p can be dependent on the membrane voltage and ionic concentrations. The functional dependencies are given in descriptions of the individual currents. Abbreviations for variables and constants are explained in Table 1. The maximum conductances
with which the currents occur in the different parts of the model neuron are given in Table 2. The current through an ion channel is then given by Ohm's law. For the Ca2+ currents Ohm's law was replaced by the Goldman-Hodgkin-Katz-equation.
|
|
The morphology of the CA1 model neuron is adapted from Varona et al. (2000)
and comprises 265 sections (829 segments) with branched basal and apical dendrite, soma, and an axon. It is based on a detailed morphometric study of average compartment dimensions, branching pattern, and tapering (Bannister and Larkman 1995
).
Passive electrophysiological properties
Passive parameters were also adapted from Varona et al. (2000)
and include values for the specific membrane capacitance, the membrane resistivity, and the resistivity of the cytoplasm. Leak currents were assumed to have a reversal potential of –70 mV.
Temperature dependence
The dependence of ion channel dynamics on the environmental temperature T can be expressed by Q(T) =
varies for different ion channels and can be different for activation (Q10,activation), inactivation (Q10,inactivation), and current amplitude (Q10,amplitude). The values for T0 are given in the description of the individual currents. The dependence of the ion channel dynamics on Q(T) was applied according to published data (see citations in the description of the individual currents). Simulations were performed for a temperature T = 30°C.
Na+ currents
The equilibrium potential for Na+ was ENa = 55 mV.
Transient Na+ current
The somatic iNaT was modeled according to Migliore et al. (1999)
![]() |
The equations describing activation were as follows
![]() |
![]() |
![]() |
m < 0.02 ms then
m=0.02 ms
![]() |
The parameter
V1/2 was used to introduce a shift in the midpoint of the activation curve. This parameter was zero for the somatic iNaT.
The equations describing fast inactivation were as follows
![]() |
![]() |
![]() |
h < 0.5 ms then
h=0.5
![]() |
The equations describing slow inactivation were as follows
![]() |
![]() |
![]() |
![]() |
s < 10 ms then
s=10 ms
We assumed T0 = 24°C, the Q10 values were derived from Migliore et al. (1999)
.
The Na+ current at the AIS was identical to the somatic Na+ current but lacked the slow inactivation process. The parameter
V1/2, which produces a shift of the activation behavior, was systematically varied as described in RESULTS.
Persistent Na+ current
The persistent Na+ current (iNaP) is a fast activating and noninactivating current
![]() |
The dynamics of the activating gating particle m are
![]() |
![]() |
iNaT In corresponds to a Na+ current with intermediate inactivation kinetics, which is observed in CA1 neurons (Yue et al. 2005
)
![]() |
The dynamics of the activation gate particle are
![]() |
![]() |
The inactivation dynamics were derived from Magistretti and Alonso (1999)
![]() |
![]() |
![]() |
![]() |
K+ currents
The equilibrium potential for K+ was EK = –95 mV.
Delayed rectifier K+ current
The delayed rectifier K+ current (iKDR) was modeled according to Golomb et al. (2006)
![]() |
with the following activation dynamics
![]() |
![]() |
A-type K+ current
The A-type K+ current (iKA) was modeled according to Golomb et al. (2006)
![]() |
The activation dynamics were as follows
![]() |
![]() |
The inactivation dynamics were as follows
![]() |
![]() |
M-type K+ current
The M current (iKM) was modeled according to Warman et al. (1994)
![]() |
We assumed T0 = 23°C. Q10 values for iKM were derived from Halliwell and Adams (1982)
.
Activation dynamics
![]() |
![]() |
![]() |
![]() |
Voltage- and Ca2+-dependent K+ current
This K+ current (iKCT) adapted from Stacey and Durand (2000)
is dependent both on the intracellular Ca2+ concentration Ca]I,1 and on the membrane potential E. For the dynamics of the Ca2+ ions see following text
![]() |
The Ca2+ dependence was implemented as follows
![]() |
![]() |
![]() |
![]() |
![]() |
The voltage dependence of gating was defined as follows
![]() |
![]() |
![]() |
![]() |
Ca2+-dependent K+ current
The gating properties of this K+ current (iKAHP) are only dependent on the intracellular Ca2+ concentration [Ca2+]1,2 and is therefore in our model restricted to the somatic compartment
![]() |
Activation dynamics
![]() |
![]() |
In both rate functions, [Ca2+]I,2 is given in mM
![]() |
![]() |
These dynamics were implemented according to Stacey and Durand (2000)
and Warman et al. (1994)
.
Ca2+ currents
The maximal permeabilities
of the various Ca2+ channels were chosen from investigations reported in Takahashi and Akaike (1991)
and Su et al. (2002)
.
T-type Ca2+ current
The T-type Ca2+ current (iCaT) is mainly based on findings reported in Lee et al. (1999)
and Klöckner et al. (1999)
![]() |
For this and the other Ca2+ currents, E < 10-4 mV 1/1 – exp(E) was approximated by the first terms of a Taylor expansion –1 + E/2 because the term exp(2FE/RT) is present at the denominator of the preceding equation, and so the denominator would become 0 when E = 0.
Activation dynamics
![]() |
![]() |
Inactivation dynamics
![]() |
![]() |
For the temperature dependence, we assumed T0 = 23°C. Q10 values for iCaT were derived from Coulter et al. (1989)
.
R-type Ca2+ current
The R-type Ca2+ current (iCaR) was modeled with current parameters taken from Sochivko et al. (2003)
and Randall and Tsien (1997)
![]() |
Activation dynamics
![]() |
![]() |
![]() |
![]() |
Inactivation dynamics
![]() |
![]() |
![]() |
![]() |
For the temperature dependence, we assumed T0. Q10 values were derived from McAllister-Williams and Kelly (1995)
.
L-type Ca2+ current
The L-type Ca2+ current (iCaL) was modeled as follows:
Activation dynamics
![]() |
![]() |
![]() |
![]() |
![]() |
For the temperature dependence, we assumed T0 = 21°C. X10 values were derived from McAllister-Williams and Kelly (1995)
.
N- and P/Q–type Ca2+ current
The high-threshold Ca2+ currents (iCanpq) mediated by the N- and P/Q-type were summarized into a single current with the following properties
![]() |
Activation dynamics
![]() |
![]() |
![]() |
![]() |
Temperature dependence: T0 = 21°C for qampl and T0 = 22°C for qm. Q10 values were derived from McAllister-Williams and Kelly (1995)
.
Ca2+ dynamics
As in Warman et al. (1994)
, the intracellular Ca2+ dynamics were modeled assuming two distinct intracellular Ca2+ pools with appropriate dynamics, given by
![]() |
The particular pool is indexed by n. Apart from the diffusion contribution, [Ca2+]n.i is changed due to the total Ca2+ current density iCa, which is the sum of iCanpq, iCaR, iCaL, and iCaT. [Ca2+]
,n denotes the intracellular Ca2+ concentration for large times and closed Ca2+ channels.
Ca,n is the associated time constant of diffusion. fn denotes the fraction of the Ca2+ current density that is active in pool n and iCa/
n is the rate of Ca2+ removal per volume. We assume an inner shell thickness
n, which is filled with Ca2+. The parameters of the two pools are as follows
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Pool 1 is present in the soma and in the dendrites; pool 2 is only present in the soma. The extracellular Ca2+ concentration was set to [Ca]o = 2 mM.
Hyperpolarization activated h-current
This unspecific cationic current ih is activated by hyperpolarization and modeled according to Gasparini et al. (2004)
![]() |
Activation dynamics
![]() |
![]() |
For the temperature dependence, we assumed T0 = 33°C. Q10 values were derived from Gasparini et al. (2004)
.
Induction of spiking
Current injections were introduced into the soma at t = 100 ms for 4 ms. Stimulus intensity was increased in steps of 0.01 nA. For analysis, we chose the lowest stimulus amplitude to which the model neuron responded with a spike to the current injection.
|
|
RESULTS |
|---|
|
NaV1.6 channels are strongly concentrated at AIS of different types of neurons in the CNS (Boiko et al. 2003
; Hossain et al. 2005
; Van Wart and Matthews 2006
; Van Wart et al. 2007
). We examined whether NaV1.6 is similarly expressed in CA1 pyramidal cells using double immunolabeling for Ankyrin G (a marker for AIS) (see for instance Garrido et al. 2003
) and for NaV1.6 in hippocampal sections. In the CA1 region, NaV1.6 subunits were clearly aggregated at the AIS (Fig. 1Aa, see insets for larger magnification of individual AIS, stratum pyramidale, oriens, and alveus indicated by SP, SO, and AL, respectively). Additionally, double immunolabeling with a PanNaV antibody and the NaV1.6 antibody revealed a concentration of both immunolabels at AIS (Fig. 1Ab). Mice lacking functional NaV1.6 channels due to a truncation mutation in exon 2 of the Scn8a gene (Scn8amed) were devoid of NaV1.6 immunoreactivity, but PanNaV immunoreactive AIS were still present (Fig. 1B). These experiments also revealed that NaV1.6 channel aggregation at AIS constitutes a general feature of cortical neurons as it was also found in dentate granule and CA3 pyramidal cells, and in subicular and neocortical neurons (Fig. 1C).
In Scn8amed mice, AIS were present in undiminished numbers and did not appear altered in immunolabelings for Ankyrin G (Fig. 2A, compare leftmost micrographs). Moreover, immunolabeling with the PanNaV antibody produced a robust signal at AIS of Scn8amed mice (Fig. 2A, compare rightmost micrographs, see also Fig. 1, Ab and B). We therefore analyzed the density of Na+ channel proteins at Scn8awt and Scn8amed mice AIS in a semi-quantitative manner. Double immunolabelings for Ankyrin G and PanNaV allowed us to demarcate individual AIS in the Ankyrin G channel (Fig. 2A, leftmost panels, see insets for larger magnification) and to determine the intensity of the PanNaV immunolabeling within this region of interest (see METHODS, Scn8awt: 8 slices from 5 mice, Scn8amed: 10 slices from 5 mice; 10 AIS were analyzed in each slice). Indeed the ratios of PanNaV to Ankyrin G immunolabeling intensities at Scn8awt and Scn8amed mice AIS were the same (1.5 ± 0.1 and 1.3 ± 0.04, respectively, P > 0.05; Fig. 2B). The lack of a significant reduction in PanNaV staining at AIS of Scn8amed mice suggests that the absence of NaV1.6 subunits allows other NaV subunits to accumulate at the AIS. This is in good agreement with Van Wart and Mathews. (2006
) and supported by the finding that NaV1.6 channels share with NaV1.1 and NaV1.2 channels Ankyrin G binding motifs that confer targeting to the axon (Garrido et al. 2003
; Pan et al. 2006
).
|
It has been previously hypothesized that INaT at the AIS activates at more negative voltages than INaT at the soma, causing spikes to commence at or close to the AIS (Colbert and Pan 2002
). This peculiarity may be due to selective accumulation of NaV1.6 channels at the AIS because these channels were shown to activate at more negative voltages than other Na+ channels when expressed in cultured dorsal root ganglion neurons (Rush et al. 2005
). If this is the case, loss of NaV1.6 channels in native CA1 neurons should lead to a depolarizing shift in INaT activation curve. To test this, we performed whole cell recordings of INaT in dissociated CA1 pyramidal neurons of Scn8awt and Scn8amed mice (n = 6 and n = 7, respectively). A representative family of INaT traces evoked by increasing voltage steps in Scn8awt (topmost traces) and Scn8amed (bottom traces) neurons are shown in Fig. 3A (voltage protocols shown in the inset). From this data, we constructed the INaT activation curve for each of the tested neurons by fitting it with a Boltzmann function (see METHODS). The peak conductance of INaT was not significantly different between the groups of neurons (Scn8awt: 69.2 ± 10.5 nS, n = 8; Scn8amed: 59.8 ± 6.5 nS, P > 0.05, n = 11). The averaged normalized data for each group of neurons are provided in Fig. 3B. We found that the INaT activation curve was
5 mV more positive in mutant neurons (V1/2 = –25.00 ± 1.18 mV) than in wild-type neurons (V1/2 = –29.77 ± 1.00 mV, P = 0.008; Fig. 3B). The steepness of the activation curve was not significantly different between the two groups (slope factor km = 5.87 ± 0.48 mV for Scn8awt and km = 5.37 ± 0.39 mV for Scn8amed neurons; P > 0.05).
|
Absence of NaV1.6 reduces the persistent Na+ current INaP
Recombinant NaV1.6 channels generate a conspicuous persistent Na+ current (INaP) component (Rush et al. 2005
), and published data suggest that these subunits underlie a significant proportion of INaP in different neuronal cell types (Do and Bean 2004
; Maurice et al. 2001
). Ramp commands (50 mV/s) applied to CA1 pyramidal neurons recorded in hippocampal slices (Fig. 4Aa) revealed a prominent inward current corresponding to INaP that was blocked by application of 1 µM TTX (Fig. 4Ab). INaP was isolated by subtracting recordings in the presence of TTX from control recordings (Fig. 4A, c and d, for Scn8awt and Scn8amed mice, respectively). The maximal INaP conductance was 1.9 ± 0.1 nS in Scn8awt neurons (n = 11) and 1.1 ± 0.2 nS in Scn8amed neurons (n = 16), corresponding to a reduction of INaP in the latter group to 58.1% of wild-type levels (P = 0.01, Fig. 4B). At the same time, the voltage- dependence of INaP was similar in the two groups (Fig. 4C; Scn8awt neurons: V1/2 = –38.6 ± 2.4 mV and km = 4.1 ± 0.3 mV, n = 11; Scn8amed neurons: V1/2 = –39.8 ± 1.3 mV and km = 3.4 ± 0.3 mV, n = 16, P > 0.05).
|
NaV1.6 subunits have been shown to contribute to resurgent Na+ (INaR) in expression systems (Smith et al. 1998
) and cerebellar neurons (Raman et al. 1997
). We first tested whether INaR is present in CA1 neurons of Scn8awt mice. Following inactivation of Na+ currents during a 15-ms prepulse to 20 mV, repolarization with a test pulse to various potentials from –100 to –10 mV gave rise to a resurgent current component within the voltage range of –50 to –10 mV (Fig. 5A, current trace at test pulse of –30 mV, in Scn8awt mouse). The magnitude of the resurgent current INaR was assessed by subtracting the steady-state current component at the end of the test pulse (Iss) from the peak of the test pulse current (Fig. 5A). Representative current families from Scn8awt and Scn8amed mice are shown in Fig. 5B, a and b, respectively. The magnitude of the resurgent current INaR proved to be significantly smaller in Scn8amed (–20.77 ± 4.1 pA, n = 11) compared with Scn8awt mice (–71.29 ± 17.04 pA, n = 10, P < 0.01, see Fig. 5C for cumulative probability plot of INaR amplitudes at –30-mV test pulses, and Fig. 5D for mean values). The voltage dependence of INaR does not appear different when comparing both genotypes (Fig. 5E). These experiments indicate that NaV1.6 subunits localized at the AIS generate resurgent currents in CA1 pyramidal cells.
|
The loss of NaV1.6 has been shown to lead to compensatory regulation of other subthreshold inward currents, notably T-type Ca2+ currents (ICaT) in Purkinje cells (Swensen and Bean 2005
). We isolated ICaT current pharmacologically in intact CA1 neurons in the slice preparation using a cocktail of Ca2+ channel blockers and TTX (see METHODS). T-type currents were discriminated on the basis of their slow deactivation kinetics in Ca2+ tail current recordings (Fig. 6A) (Sochivko et al. 2002
). T-type current mediated tail current amplitudes in CA1 were not different at all tested command voltages (Fig. 6B). For instance, average maximal current amplitudes were –388.51 ± 55.9 pA in Scn8amed (n = 8) and –373.30 ± 106.7 pA in Scn8awt neurons (n = 7).
|
The pronounced depolarizing shift in the voltage dependence of the transient Na+ current INaT predicts a depolarizing shift in spike threshold. To test this prediction, we performed whole cell current-clamp recordings in CA1 pyramidal cells in the slice preparation. Spikes were evoked by injecting brief (4 ms) depolarizing current pulses from a membrane potential of –80 mV imposed with slow current clamp (see METHODS, Fig. 7, A and B). Spike thresholds were significantly more depolarized in Scn8amed compared with Scn8awt neurons (–56.7 ± 1.0 mV, n = 14 compared with –60.4 ± 0.9 mV; n = 22), respectively. This corresponds to a statistically significant 3.7 mV shift (Fig. 7C; P = 0.011). Changes of similar magnitude were also observed when spikes were elicited from other holding potentials within the range of –65 to –80 mV (Fig. 7D), and, for instance, amounted to 4.9 mV for spikes elicited from –70 mV. We also measured other parameters of single spikes. When spikes were elicited by brief current injections, spike amplitude and the maximal rate of depolarization during spike upstroke were the same in the two groups of neurons (118.5 ± 0.4 mV and 419.6 ± 5.1 mV/ms in Scn8awt and 117.7 ± 0.7 mV and 405.9 ± 5.8 mV/ms in Scn8amed neurons, P > 0.05), as expected from the lack of difference in maximal Na+ conductance. We did find a statistically significant, albeit small, increase in the maximal rate of spike repolarization in Scn8amed versus Scn8awt neurons (–86.2 ± 1.9 vs. –78.2 ± 2.0 mV/ms, respectively, for spikes evoked from a holding potential of –80 mV; P = 0.007). The active spike afterdepolarization (spike ADP, see METHODS) was not different when comparing Scn8awt (168.3 ± 6.8 mV * ms, n = 22) and Scn8amed neurons (149.3 ± 11.3 mV * ms, n = 14).
|
|
Both the depolarizing shift in the spike threshold, as well as potentially the diminished resurgent Na+ current (Raman and Bean 1997
; Raman et al. 1997
) would be expected to reduce the spike gain of CA1 neurons in Scn8amed mice. We therefore tested whether spike gain is affected by applying prolonged (600 ms) depolarizing current pulses of increasing magnitude (from 20 to 120 pA) and examining the number of spikes evoked by equivalent current injection steps in seven Scn8awt and 5 Scn8amed neurons (Fig. 8, A and B, respectively). Indeed the relation of current injection to the corresponding spike frequency was significantly steeper in Scn8awt compared with Scn8amed neurons (P < 0.01; Fig. 8C).
NaV1.6 contributes to axonal spike initiation
Spike initiation occurs within the axon in most types of cortical neurons (Colbert and Pan 2002
; Khaliq and Raman 2006
; Palmer and Stuart 2006
; Stuart and Hausser 1994
; Stuart and Sakmann 1994
; Stuart et al. 1997
), and more precise attempts at localization have revealed an initiation site at the most distal portion of the AIS in layer 5 cortical pyramidal neurons (Palmer and Stuart 2006
) and CA3 pyramidal neurons (Meeks and Mennerick 2007
). In action potentials elicited by prolonged current injection, phase plots (dV/dt vs. V) allowed to distinguish a first phase of spike upstroke due to spike propagation from the AIS into the soma (McCormick et al. 2007
; Shu et al. 2007
), and a second phase, caused by the somatic spike (Fig. 9, A and B). This phenomenon was observed both for the first spike as well as for spikes occurring later during prolonged (600 ms) current injections (Fig. 9, A: 1st spike in train; B: 5th spike in train, multiple spikes from individual cells are shown). The initiation of spikes in Scn8awt neurons (n = 7; Fig. 9, A and B, top traces) appeared more abrupt than in Scn8amed mice (n = 5; Fig. 9, A and B, bottom traces). This abrupt initiation was previously described in neocortical neurons as "kink" and is a consequence of the invasion of the soma by an axonal spike (McCormick et al. 2007
; Shu et al. 2007
). The "abruptness" of the voltage change at the onset of a spike can be quantified as a maximum of the second derivation of the voltage trace. We calculated the second derivation of the voltage traces (see Coombs et al. 1957
); voltage recordings in Fig. 10, Aa and Ba, first and second derivation in Ab and Bb, second derivation depicted in gray, corresponding to the rate of change of dV/dt), in which a first (axonal) component and the second (somatic) component could be discriminated (Fig. 10, Ab and Bb). When this analysis was performed, the amplitude of the first peak in the second derivation of the voltage traces was significantly smaller in Scn8amed neurons, regardless of which spike in a train was evaluated (Fig. 10C, comparable results obtained when action potentials were binned into 100-ms bins according to the time of occurrence after onset of the current injection, data not shown), reflecting the less abrupt rise of the voltage trace at the initiation of spikes seen in the phase plots (see Fig. 9, Ab and Bb).
|
|
Computer simulations of spike initiation at the AIS
Our electrophysiological results described above strongly suggest that in CA1 pyramidal cells, the high density of NaV1.6 channels imposes a low spike threshold at the AIS so that spikes are initiated in this region before they appear in the soma. Another factor that may influence spike threshold and spike trigger zone is the overall density of Na+ channels at the AIS compared with that at the soma. Studies using cell-attached patch-clamp recordings to compare INaT densities at AIS versus soma membranes have reported either equal densities (Colbert and Johnston 1996
; Colbert and Pan 2002
) or a much higher densities at the AIS (Kole et al. 2008
). To explore the consequences of systematically altering INaT density and/or its voltage dependence on spike threshold and trigger zone, we performed simulations in a realistic computer model of a CA1 neuron (see Fig. 11C for morphology; see METHODS for detailed description of conductances). This approach also allowed us to directly compare voltage traces at axonal and somatic sites. The incorporated in this model is shown in Fig. 11A (see METHODS for parameters). This current was incorporated in the axonal and somatic compartments. We then varied the voltage of half-maximal activation (V1/2) systematically at the AIS, such that it was
7 mV more hyperpolarized than at the soma (
V1/2: shift of V1/2 of activation relative to somatic iNaT, activation curves are depicted for
V1/2 of 0, –4 and –7 mV shown in Fig. 11B). As a second parameter, we varied iNaT density at the AIS. Figure 11D shows exemplary somatic spikes elicited by brief current injection at the soma of the model neuron (iNaT densities at the AIS and at the soma were equal;
V1/2 was 0 and –7 mV, as indicated, detailed description of spike properties for different iNaT densities and
V1/2 in Supplementary Fig. S2).1
|
V1/2 (0 to –7 mV) and iNaT density at the AIS (from 0.02 to 1 S/cm2, corresponding to a 0.2- to 10-fold difference in iNaT density relative to the somatic iNaT density of 0.1 S/cm2). The axo-somatic delay was then calculated as the delay between the time points at which the slope of rise in both compartments was maximal. A delay could also be derived from somatic voltage traces alone in our model, similar to the in vitro recordings. Derivations of simulated somatic voltage traces also revealed two distinct peaks under most conditions. The values of the axo-somatic delay obtained in this manner from the somatic recording alone showed a strong linear correlation to the values derived as a delay between AIS and somatic spikes (R2 = 0.9388).
|
V1/2. A pronounced delay from axonal to somatic spike initiation was observed at values of
V1/2 from –7 to –4 mV. When
V1/2 was reduced further, the axo-somatic delay showed a steep reduction (examples for
V1/2 of 0 and –7 mV in Fig. 12Ab, results for all values of
V1/2 in Fig. 12Ba, gray data points). A higher iNaT density at the AIS as suggested by Kole et al. (2008)
10-fold increase relative to the soma implemented in our model) always led to a spike initiation at the AIS, and a stereotypical axo-somatic delay of
0.15 ms, irrespective of
V1/2 (Figs. 12Ac and 9Ba, black symbols). Conversely, a reduced iNaT density at the AIS (0.2-fold of somatic iNaT density) caused the spike to arise almost simultaneously in both compartments for all values of
V1/2 (Fig. 12, Aa and Ba, open symbols). Thus a
V1/2 of more than –4 mV strongly promotes spike initiation at the AIS, even when the iNaT densities at the AIS and soma were uniform. This phenomenon was also clear when we plotted the axo-somatic delay versus the relative iNaT density at the AIS (Fig. 12Ca). This analysis revealed that for
V1/2 of 0 mV, the axo-somatic delay increased gradually with an increasing density of axonal iNaT. When
V1/2 was increased, this relation began to show a steeper increase. As a consequence, a
V1/2 of –4 to –7 mV strongly affected spike initiation site over a wide range of AIS Na+ channel density ratios (from
0.5-fold to 3-fold somatic density, Fig. 12Ca).
The voltage dependence of activation of iNaT at the AIS also influenced spike threshold as observed experimentally. When iNaT densities at the AIS and soma were equal, the firing threshold was dependent on
V1/2, such that an increase in
V1/2 led to a more hyperpolarized spike threshold (examples for
V1/2 of 0 and –7 mV in Fig. 12Ab, results for all values of
V1/2 in Fig. 12Bb, gray data points). At a very high iNaT density at the AIS, spike threshold was always hyperpolarized, irrespective of
V1/2 (Fig. 12, Ac and Bb, black symbols). Conversely, very low iNaT density at the AIS led to a depolarized spike threshold without dependence on
V1/2 (Fig. 12, Aa and Bb, open symbols).
In Scn8amed mice, we observed a significant reduction of INaP and INaR current. Of these two current components, INaP might conceivably contribute to action potential initiation. We have therefore repeated the modeling experiment with iNaP reduced to 60% in all compartments in which it was present (soma: reduction to 0.6 mS/cm2, AIS: 0.3 mS/cm2, Fig. 13). In additional experiments, we reduced iNaP only at the AIS (Supplementary Fig. S1). Under both conditions, the impact of varying iNaT was similar to those depicted in Fig. 12. In both cases, changing the voltage dependence of activation of iNaT at the AIS still influenced the axo-somatic delay (Fig. 13Aa and Supplementary Fig. S1Aa) and spike threshold (Fig. 13Ab and Supplementary Fig. S1Ab). Varying the density of iNaT at the AIS also caused changes in axo-somatic delay and spike threshold that were well comparable to the data obtained without reduction in iNaP (Fig. 13B and Supplementary Fig. S1B, cf. Fig. 12C).
|
|
|
DISCUSSION |
|---|
|
Regarding spike initiation, two major changes were observed in Scn8amed mice. First we observed a significant depolarizing shift in spike threshold in mice lacking NaV1.6 channels. In addition, deletion of NaV1.6 channels from the AIS significantly reduced the temporal separation between axonal and somatic components of spike initiation in repetitive firing. Previous studies have shown that spike initiation occurs within the distal portion of the AIS in cortical neurons (Meeks and Mennerick 2007
; Palmer and Stuart 2006
) or the first node of Ranvier in Purkinje neurons (Clark et al. 2005
). Interplay between several factors likely endows these subcellular compartments with a particularly low spike threshold. First the passive electrical properties of axon versus soma may play an important role. Modeling and physiological studies suggest that charging of the AIS capacitance by inward current is rapid with the much larger somatic capacitance being charged with a significant delay (McCormick et al. 2007
; Meeks and Mennerick 2007
; Shu et al. 2007
). Second, a high density of AIS Na+ channels was suggested to subserve AIS spike initiation in modeling and electrophysiological studies. Several studies have shown a high density of Na+ channel proteins at the AIS (Boiko et al. 2001
, 2003
; Catterall 1981
; Hossain et al. 2005
; Pan et al. 2006
; Van Wart and Matthews 2006
), but how far this correlates with AIS Na+ current density is a matter of current debate (Colbert and Pan 2002
; Kole et al. 2008
; Palmer and Stuart 2006
). Finally, the more negative activation voltages of AIS Na+ channels are thought to lower spike threshold (Colbert and Pan 2002
). Clearly these factors are not mutually exclusive; rather, it is likely that these three factors in combination localize the spike trigger zone to the AIS. The most likely interpretation of the reduced axo-somatic delay in our view is that the site of spike initiation is located closer to the soma. This is also suggested by the modeling data, where removing the voltage shift of INaT caused a simultaneous spike initiation in soma and AIS (cf. Fig. 12Ab, equal density of INaT at AIS and soma).
In mice lacking the AIS Na+ channel subunit Nav1.6, we found a pronounced depolarizing shift in the half-maximal activation of INaT in CA1 neurons. This finding is consistent with studies that have examined the properties of NaV1.2 or NaV1.6 channels by overexpressing them in mammalian cells. These experiments have indicated that the activation curve of NaV1.6 channels is shifted in a hyperpolarized direction compared with NaV1.2 (Rush et al. 2005
). It should be noted that such a shift was not observed when NaV subunits were expressed in oocytes, for unknown reasons (Smith et al. 1998
). A shift in the voltage dependence of activation was also not observed in globus pallidus neurons (Mercer et al. 2007
), cerebellar neurons (Raman et al. 1997
) or mesencephalic trigeminal neurons (Enomoto et al. 2007
) from Scn8amed mice. Regarding the voltage-dependence of inactivation, a more hyperpolarized voltage dependence of INaT was observed for NaV1.6 channels compared with NaV1.2 channels (Rush et al. 2005
), but no changes in this biophysical parameter were observed in different cell types in mice lacking functional NaV1.6 channels (Enomoto et al. 2007
; Mercer et al. 2007
; and this study, but see Raman et al. 1997
). The reasons for these disparate findings are currently unknown but may indicate both cell-specific regulation of NaV1.6 channels as well as potential compensatory changes following loss of NaV1.6 channels. Regardless of these discrepancies, our results indicate that in CA1 neurons, NaV1.6 subunits contribute a Na+ channel component that activates at more hyperpolarized voltages than the remainder of the cellular Na+ currents. Our and published immunolabeling experiments (Boiko et al. 2003
; Garrido et al. 2003
; Van Wart and Matthews 2006
; Van Wart et al. 2007
) indicate that these channels are located at the AIS of different types of principal neurons, suggesting that they may underlie biophysical specialization of AIS Na+ channels (Colbert and Pan 2002
). It should be noted, however, that our recordings of the biophysical properties of INaT in Scn8amed and Scn8awt mice were performed in dissociated CA1 neurons, which may contain variable portions of axonal membrane. We cannot therefore exclude that NaV1.6 channels at the AIS might have properties distinct from somatic NaV1.6 channels, perhaps via specific interactions with AIS proteins (Shirahata et al. 2006
). Nevertheless, the most parsimonious explanation for our results is that NaV1.6 channels with a hyperpolarized threshold of activation aggregate at the AIS.
We did not quantitatively assess if the density of Na+ channels at the AIS is altered in CA1 neurons and therefore cannot exclude a reduction in the overall density of AIS Na+ channels in SCN8amed neurons. However, our immunohistochemical data suggest that there is no dramatic loss of AIS Na+ channels in these neurons. Relative to Ankyrin G as an AIS marker, we did not observe a reduction in PanNaV immunolabeling in Scn8amed neurons. This is similar to the results reported by van Wart et al. (2006)
, indicating a compensation of the loss of NaV1.6 subunits at the AIS by other subunits, in particular NaV1.2. A mild reduction in Na+ channel density might not be detected using immunolabeling, but would be unlikely to exclusively account for the observed changes in spike initiation.
Our modeling data allowed us to further address the interplay of the density and the voltage dependence of AIS Na+ channels in spike initiation. We show that a hyperpolarized voltage dependence of AIS Na+ currents influences spike initiation over a wide range of AIS Na+ channel densities (from
0.5- to 3-fold of somatic density), If the density of Na+ channels at the AIS becomes even higher, the initiation site is less affected by the biophysical properties of these channels. The threshold for generation of a spike was differently affected by altering AIS Na+ channels. In this case, even at very high AIS Na+ channel densities (
10x somatic density), a shift in voltage-dependent Na+ channel activation still influenced spike threshold (see Fig. 11Cb). At the same time, increasing the density of AIS channels always led to a more hyperpolarized somatic spike threshold. Thus the effects of varying the voltage dependence of AIS Na+ channels on spike threshold and spike trigger zone were robust over a large range of AIS Na+ current densities. These data indicate that the biophysical properties of AIS INaT are an important determinant of spike threshold and are consistent with the view that the voltage dependence of AIS NaV1.6 is an important factor in spike initiation of CA1 pyramidal neurons. In addition to the changes in INaT, we also found a reduction of INaP in Scn8amed mice. It is conceivable that NaV1.6-mediated INaP could, by virtue of its hyperpolarized threshold of activation, contribute to spike initiation. However, modeling experiments showed that the influence of this current component on spike threshold and axo-somatic delay is likely to be much smaller than the influence of INaT.
The changes in INaP (by 41%) and INaR (by 69.2%) we observed in Scn8amed mice are similar to the results reported for mesencephalic trigeminal neurons in NaV1.6 null mice (39% reduction in INaP, 76% reduction in INaR) (Enomoto et al. 2007
), DRG neuron cultures (complete ablation of INaR) (Cummins et al. 2005
), subthalamic nucleus neurons (63% reduction in INaR, 55% reduction in INaP) (Do and Bean 2004
), or cerebellar neurons (Raman and Bean 1997
). Globus pallidus neurons in mice lacking NaV1.6, surprisingly, show no reduction in INaP, but INaR is reduced (Mercer et al. 2007
). Taken together, these results suggest that a significant portion of INaP and INaR is mediated by axonal NaV1.6 channels. In addition to these neuron types in the cerebellum, diencephalon and brain stem, the presence of INaR was also reported in cortical pyramidal neurons of the perirhinal and entorhinal cortex, as well as in dentate granule cells and CA1 pyramidal neurons of ventral hippocampus (Castelli et al. 2007a,b). Both INaR and INaP mediated by NaV1.6 have been shown to affect repetitive firing and spike output gain (Levin et al. 2006
; Mercer et al. 2007
; Raman et al. 1997
). In addition, the changes in spike threshold would also be expected to have a similar effect. Indeed we also found a large reduction in spike output gain in Scn8amed compared with Scn8awt mice. It is likely that the changes in INaR, INaP, and INaT conspire in CA1 neurons to produce changes in output gain. These results are also interesting because they imply that a substantial portion of INaP and INaR may be generated at the AIS of different types of central neurons, as shown with physiological techniques (Astman et al. 2006
; Castelli et al. 2007a
).
INaP has also been shown to contribute strongly to spike afterdepolarizations in CA1 pyramidal neurons from adult animals (Yue et al. 2005
). In young animals comparable to the age range employed in this study, not only INaP but also dendritic voltage-gated Ca2+ currents strongly amplify spike afterdepolarizations and cause the generation of spike bursts (Chen et al. 2005
). In this age range, blocking either voltage-gated Ca2+ currents at the dendrites or INaP in the perisomatic region pharmacologically reduces spike afterdepolarizations and associated burst discharges. Surprisingly, spike afterdepolarizations were not reduced in Scn8amed mice despite a reduction of INaP by 41.9%. One explanation for this unexpected finding might be that a partial reduction of INaP in young animals is not sufficient to affect the magnitude of the spike afterdepolarization, given the important contribution of voltage-gated Ca2+ currents at this age (Chen et al. 2005
). An alternative explanation would be compensatory regulation of other voltage-gated ion channels occurring as a consequence of the constitutive lack of function of NaV1.6. Indeed, functional deletion of NaV1.6 in Scn8amed mice causes compensatory upregulation of T-type Ca2+ channels in Purkinje neurons (Swensen and Bean 2005
). In contrast, changes in K+ channels were subtle, with only small changes in the voltage dependence of K+ currents highly sensitive to TEA in Scn8amed mice (Khaliq et al. 2003
). We did not find a compensatory upregulation of T-type Ca2+ channels, indicating that different compensatory changes may be invoked in different neuron types.
Taken together, our results indicate that the presence of NaV1.6 endows AIS Na+ channels with a hyperpolarized voltage dependence of activation that is important for the low threshold for spike initiation at the AIS. Furthermore, axonal NaV1.6 channels contribute to INaP and INaR. The contribution of NaV1.6 to these three current components plays a significant role in regulating neuronal repetitive discharge behavior. Our findings may be pertinent to many other types of brain neurons because NaV1.6 subunit aggregation at the AIS has been demonstrated in neocortical, subicular, and hippocampal pyramidal neurons (Van Wart and Matthews 2006
and this study), as well as in cochlear (Hossain et al. 2005
), retinal ganglion (Boiko et al. 2003
), and Purkinje cells (Van Wart and Matthews 2006
). The role of NaV1.6 in controlling neuronal firing behavior is consistent with the elevated seizure thresholds observed in heterozygous Scn8amed/wt mice (Martin et al. 2007
). This study also suggests that reduced function of Scn8a limits hyperexcitability in a mouse model of severe myoclonic epilepsy of infancy, suggesting a role for this gene as a disease modifier in epilepsy. This study further underscores the important role of NaV1.6 channels in controlling neuronal excitability on a systems level.
|
|
GRANTS |
|---|
|
|
|
FOOTNOTES |
|---|
1 The online version of this article contains supplemental data. ![]()
Address for reprint requests and other correspondence: H. Beck, Dept. of Epileptology, University of Bonn, Sigmund-Freud-Str. 25, 53105 Bonn, Germany (E-mail: heinz.beck{at}ukb.uni-bonn.de)
|
|
REFERENCES |
|---|
|
Bannister NJ, Larkman AU. Dendritic morphology of CA1 pyramidal neurones from the rat hippocampus. I. Branching patterns. J Comp Neurol 360: 150–160, 1995.[CrossRef][Web of Science][Medline]
Boiko T, Rasband MN, Levinson SR, Caldwell JH, Mandel G, Trimmer JS, Matthews G. Compact myelin dictates the differential targeting of two sodium channel isoforms in the same axon. Neuron 30: 91–104, 2001.[CrossRef][Web of Science][Medline]
Boiko T, Van Wart A, Caldwell JH, Levinson SR, Trimmer JS, Matthews G. Functional specialization of the axon initial segment by isoform-specific sodium channel targeting. J Neurosci 23: 2306–2313, 2003.
Burgess DL, Kohrman DC, Galt J, Plummer NW, Jones JM, Spear B, Meisler MH. Mutation of a new sodium channel gene, Scn8a, in the mouse mutant "motor endplate disease." Nat Genet 10: 461–465, 1995.[CrossRef][Web of Science][Medline]
Carnevale NT, Hines M. In: The NEURON Book. Cambridge, UK: Cambridge Univ. Press, 2006.
Castelli L, Biella G, Toselli M, Magistretti J. Resurgent Na+ current in pyramidal neurons of rat perirhinal cortex: axonal location of channels and contribution to depolarizing drive during repetitive firing. J Physiol 582: 1179–1193, 2007s.
Castelli L, Nigro MJ, Magistretti J. Analysis of resurgent sodium-current expression in rat parahippocampal cortices and hippocampal formation. Brain Res 1163: 44–55, 2007b.[CrossRef][Web of Science][Medline]
Catterall WA. Localization of sodium channels in cultured neural cells. J Neurosci 1: 777–783, 1981.[Abstract]
Chen S, Yue C, Yaari Y. A transitional period of Ca2+-dependent spike afterdepolarization and bursting in developing rat CA1 pyramidal cells. J Physiol 567: 79–93, 2005.
Clark BA, Monsivais P, Branco T, London M, Hausser M. The site of action potential initiation in cerebellar Purkinje neurons. Nat Neurosci 8: 137–139, 2005.[CrossRef][Web of Science][Medline]
Colbert CM, Johnston D. Axonal action-potential initiation and Na+ channel densities in the soma and axon initial segment of subicular pyramidal neurons. J Neurosci 16: 6676–6686, 1996.
Colbert CM, Pan E. Ion channel properties underlying axonal action potential initiation in pyramidal neurons. Nat Neurosci 5: 533–538, 2002.[CrossRef][Web of Science][Medline]
Coombs JS, Curtis DR, Eccles JC. The generation of impulses in motoneurons. J Physiol 139: 232–249, 1957.
Coulter DA, Huguenard JR, Prince DA. Calcium currents in rat thalamocortical relay neurons: kinetic properties of the transient, low-threshold current. J Physiol 414: 587–604, 1989.
Cummins TR, Dib-Hajj SD, Herzog RI, Waxman SG. Nav1.6 channels generate resurgent sodium currents in spinal sensory neurons. FEBS Lett 579: 2166–2170, 2005.[CrossRef][Web of Science][Medline]
Do MT, Bean BP. Sodium currents in subthalamic nucleus neurons from Nav1.6-null mice. J Neurophysiol 92: 726–733, 2004.
Enomoto A, Han JM, Hsiao CF, Chandler SH. Sodium currents in mesencephalic trigeminal neurons from Nav1.6 null mice. J Neurophysiol 98: 710–719, 2007.
Garrido JJ, Giraud P, Carlier E, Fernandes F, Moussif A, Fache MP, Debanne D, Dargent B. A targeting motif involved in sodium channel clustering at the axonal initial segment. Science 300: 2091–2094, 2003.
Gasparini S, Migliore M, Magee JC. On the initiation and propagation of dendritic spikes in CA1 pyramidal neurons. J Neurosci 24: 11046–11056, 2004.
Golomb D, Yue C, Yaari Y. Contribution of persistent Na+ current and M-type K+ current to somatic bursting in CA1 pyramidal cells: combined experimental and modeling study. J Neurophysiol 96: 1912–1926, 2006.
Halliwell JV, Adams PR. Voltage-clamp analysis of muscarinic excitation in hippocampal neurons. Brain Res 250: 71–92, 1982.[CrossRef][Web of Science][Medline]
Hamann M, Meisler MH, Richter A. Motor disturbances in mice with deficiency of the sodium channel gene Scn8a show features of human dystonia. Exp Neurol 184: 830–838, 2003.[CrossRef][Web of Science][Medline]
Hodgkin AL, Huxley AF. A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol 117: 500–544, 1952.
Hossain WA, Antic SD, Yang Y, Rasband MN, Morest DK. Where is the spike generator of the cochlear nerve? Voltage-gated sodium channels in the mouse cochlea. J Neurosci 25: 6857–68, 2005.
Khaliq ZM, Gouwens NW, Raman IM. The contribution of resurgent sodium current to high-frequency firing in Purkinje neurons: an experimental and modeling study. J Neurosci 23: 4899–4912, 2003.
Khaliq ZM, Raman IM. Relative contributions of axonal and somatic Na channels to action potential initiation in cerebellar Purkinje neurons. J Neurosci 26: 1935–1944, 2006.
Klöckner U, Lee JH, Cribbs LL, Daud A, Hescheler J, Pereverzev A, Perez-Reyes E, Schneider T. Comparison of the Ca2+ currents induced by expression of three cloned
1 subunits,
1G,
1H and
1I, of low-voltage-activated T-type Ca2+ channels. Eur J Neurosci 11: 4171–4178, 1999.[CrossRef][Web of Science][Medline]
Kohrman DC, Harris JB, Meisler MH. Mutation detection in the med and medJ alleles of the sodium channel Scn8a. Unusual splicing due to a minor class AT-AC intron. J Biol Chem 271: 17576–81, 1996.
Kole MH, Ilschner SU, Kampa BM, Williams SR, Ruben PC, Stuart GJ. Action potential generation requires a high sodium channel density in the axon initial segment. Nat Neurosci 11: 178–186, 2008.[CrossRef][Web of Science][Medline]
Lee JH, Daud AN, Cribbs LL, Lacerda AE, Pereverzev A, Klockner U, Schneider T, Perez-Reyes E. Cloning and expression of a novel member of the low voltage-activated T- type calcium channel family. J Neurosci 19: 1912–1921, 1999.
Levin SI, Khaliq ZM, Aman TK, Grieco TM, Kearney JA, Raman IM, Meisler MH. Impaired motor function in mice with cell-specific knockout of sodium channel Scn8a (NaV1.6) in cerebellar purkinje neurons and granule cells. J Neurophysiol 96: 785–793, 2006.
Magistretti J, Alonso A. Biophysical properties and slow voltage-dependent inactivation of a sustained sodium current in entorhinal cortex layer-II principal neurons: a whole-cell and single-channel study. J Gen Physiol 114: 491–509, 1999.
Martin MS, Tang B, Papale LA, Yu FH, Catterall WA, Escayg A. The voltage-gated sodium channel Scn8a is a genetic modifier of severe myoclonic epilepsy of infancy. Hum Mol Genet 16: 2892–2899, 2007.
Maurice N, Tkatch T, Meisler M, Sprunger LK, Surmeier DJ. D1/D5 dopamine receptor activation differentially modulates rapidly inactivating and persistent sodium currents in prefrontal cortex pyramidal neurons. J Neurosci 21: 2268–2277, 2001.
McAllister-Williams RH, Kelly JS. The temperature dependence of high-threshold calcium channel currents recorded from adult rat dorsal raphe neurons. Neuropharmacology 34: 1479–1490, 1995.[CrossRef][Web of Science][Medline]
McCormick DA, Shu Y, Yu Y. Neurophysiology: Hodgkin and Huxley model–still standing? Nature 445: E1–E2, 2007.[CrossRef][Web of Science][Medline]
Meeks JP, Mennerick S. Action potential initiation and propagation in CA3 pyramidal axons. J Neurophysiol 97: 3460–3472, 2007.
Mercer JN, Chan CS, Tkatch T, Held J, Surmeier DJ. Nav1.6 sodium channels are critical to pacemaking and fast spiking in globus pallidus neurons. J Neurosci 27: 13552–13566, 2007.
Migliore M, Hoffman DA, Magee JC, Johnston D. Role of an A-type K+ conductance in the back-propagation of action potentials in the dendrites of hippocampal pyramidal neurons. J Comput Neurosci 7: 5–15, 1999.[CrossRef][Web of Science][Medline]
Naundorf B, Wolf F, Volgushev M. Unique features of action potential initiation in cortical neurons. Nature 440: 1060–1063, 2006.[CrossRef][Web of Science][Medline]
Noda M, Ikeda T, Suzuki H, Takeshima H, Takahashi T, Kuno M, Numa S. Expression of functional sodium channels from cloned cDNA. Nature 322: 826–828, 1986.[CrossRef][Web of Science][Medline]
Ogiwara I, Miyamoto H, Morita N, Atapour N, Mazaki E, Inoue I, Takeuchi T, Itohara S, Yanagawa Y, Obata K, Furuichi T, Hensch TK, Yamakawa K. Na(v)1.1 localizes to axons of parvalbumin-positive inhibitory interneurons: a circuit basis for epileptic seizures in mice carrying an scn1a gene mutation. J Neurosci 27: 5903–5914, 2007.
Palmer LM, Stuart GJ. Site of action potential initiation in layer 5 pyramidal neurons. J Neurosci 26: 1854–1863, 2006.
Pan Z, Kao T, Horvath Z, Lemos J, Sul JY, Cranstoun SD, Bennett V, Scherer SS, Cooper EC. A common ankyrin-G-based mechanism retains KCNQ and NaV channels at electrically active domains of the axon. J Neurosci 26: 2599–2613, 2006.
Raman IM, Bean BP. Resurgent sodium current and action potential formation in dissociated cerebellar Purkinje neurons. J Neurosci 17: 4517–4526, 1997.
Raman IM, Sprunger LK, Meisler MH, Bean BP. Altered subthreshold sodium currents and disrupted firing patterns in Purkinje neurons of Scn8a mutant mice. Neuron 19: 881–891, 1997.[CrossRef][Web of Science][Medline]
Randall AD, Tsien RW. Contrasting biophysical and pharmacological properties of T-type and R- type calcium channels. Neuropharmacology 36: 879–893, 1997.[CrossRef][Web of Science][Medline]
Rush AM, Dib-Hajj SD, Waxman SG. Electrophysiological properties of two axonal sodium channels, Nav1.2 and Nav1.6, expressed in mouse spinal sensory neurones. J Physiol 564: 803–815, 2005.
Shirahata E, Iwasaki H, Takagi M, Lin C, Bennett V, Okamura Y, Hayasaka K. Ankyrin-G regulates inactivation gating of the neuronal sodium channel, Nav1.6. J Neurophysiol 96: 1347–1357, 2006.
Shu Y, Duque A, Yu Y, Haider B, McCormick DA. Properties of action-potential initiation in neocortical pyramidal cells: evidence from whole cell axon recordings. J Neurophysiol 97: 746–760, 2007.
Smith MR, Smith RD, Plummer NW, Meisler MH, Goldin AL. Functional analysis of the mouse Scn8a sodium channel. J Neurosci 18: 6093–6102, 1998.
Sochivko D, Chen J, Becker A, Beck H. Blocker-resistant Ca2+ currents in rat CA1 hippocampal pyramidal neurons. Neuroscience 116: 629–638, 2003.[CrossRef][Web of Science][Medline]
Sochivko D, Pereverzev A, Smyth N, Gissel C, Schneider T, Beck H. The Cav2.3 calcium channel subunit contributes to R-type calcium currents in murine hippocampal and neocortical neurones. J Physiol 542.3: 600–710, 2002.
Stacey WC, Durand DM. Stochastic resonance improves signal detection in hippocampal CA1 neurons. J Neurophysiol 83: 1394–1402, 2000.
Stuart G, Hausser M. Initiation and spread of sodium action potentials in cerebellar Purkinje cells. Neuron 13: 703–712, 1994.[CrossRef][Web of Science][Medline]
Stuart G, Schiller J, Sakmann B. Action potential initiation and propagation in rat neocortical pyramidal neurons. J Physiol 505: 617–632, 1997a.
Stuart GJ, Sakmann B. Active propagation of somatic action potentials into neocortical pyramidal cell dendrites. Nature 367: 69–72, 1994.[CrossRef][Web of Science][Medline]
Su H, Sochivko D, Becker A, Chen J, Jiang Y, Yaari Y, Beck H. Upregulation of a T-Type Ca2+ channel causes a long-lasting modification of neuronal firing mode after status epilepticus. J Neurosci 22: 3645–3655, 2002.
Swensen AM, Bean BP. Robustness of burst firing in dissociated purkinje neurons with acute or long-term reductions in sodium conductance. J Neurosci 25: 3509–3520, 2005.
Takahashi K, Akaike N. Calcium antagonist effects on low-threshold (T-type) calcium current in rat isolated hippocampal CA1 pyramidal neurons. J Pharmacol Exp Ther 256: 169–175, 1991.
Van Wart A, Matthews G. Impaired firing and cell-specific compensation in neurons lacking nav1.6 sodium channels. J Neurosci 26: 7172–7180, 2006.
Van Wart A, Trimmer JS, Matthews G. Polarized distribution of ion channels within microdomains of the axon initial segment. J Comp Neurol 500: 339–352, 2007.[CrossRef][Web of Science][Medline]
Varona P, Ibarz JM, Lopez-Aguado L, Herreras O. Macroscopic and subcellular factors shaping population spikes. J Neurophysiol 83: 2192–2208, 2000.
Warman EN, Durand DM, Yuen GLF. Reconstruction of hippocampal CA1 pyramidal cell electrophysiology by computer simulation. Neurophysiology 2033–2045, 1994.
Yue C, Remy S, Su H, Beck H, Yaari Y. Proximal persistent Na+ channels drive spike afterdepolarizations and associated bursting in adult CA1 pyramidal cells. J Neurosci 25: 9704–9720, 2005.
Zhou D, Lambert S, Malen PL, Carpenter S, Boland LM, Bennett V. AnkyrinG is required for clustering of voltage-gated Na channels at axon initial segments and for normal action potential firing. J Cell Biol 143: 1295–1304, 1998.
This article has been cited by other articles:
![]() |
J. Golowasch, G. Thomas, A. L. Taylor, A. Patel, A. Pineda, C. Khalil, and F. Nadim Membrane Capacitance Measurements Revisited: Dependence of Capacitance Value on Measurement Method in Nonisopotential Neurons J Neurophysiol, October 1, 2009; 102(4): 2161 - 2175. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |