Ca2+ entry through NaV channels generates submillisecond axonal Ca2+ signaling
Abstract
Calcium ions (Ca2+) are essential for many cellular signaling mechanisms and enter the cytosol mostly through voltage-gated calcium channels. Here, using high-speed Ca2+ imaging up to 20 kHz in the rat layer five pyramidal neuron axon we found that activity-dependent intracellular calcium concentration ([Ca2+]i) in the axonal initial segment was only partially dependent on voltage-gated calcium channels. Instead, [Ca2+]i changes were sensitive to the specific voltage-gated sodium (NaV) channel blocker tetrodotoxin. Consistent with the conjecture that Ca2+ enters through the NaV channel pore, the optically resolved ICa in the axon initial segment overlapped with the activation kinetics of NaV channels and heterologous expression of NaV1.2 in HEK-293 cells revealed a tetrodotoxin-sensitive [Ca2+]i rise. Finally, computational simulations predicted that axonal [Ca2+]i transients reflect a 0.4% Ca2+ conductivity of NaV channels. The findings indicate that Ca2+ permeation through NaV channels provides a submillisecond rapid entry route in NaV-enriched domains of mammalian axons.
eLife digest
Nerve cells communicate using tiny electrical impulses called action potentials. Special proteins termed ion channels produce these electric signals by allowing specific charged particles, or ions, to pass in or out of cells across its membrane. When a nerve cell ‘fires’ an action potential, specific ion channels briefly open to let in a surge of positively charged ions which electrify the cell. Action potentials begin in the same place in each nerve cell, at an area called the axon initial segment. The large number of sodium channels at this site kick-start the influx of positively charged sodium ions ensuring that every action potential starts from the same place.
Previous research has shown that, when action potentials begin, the concentration of calcium ions at the axon initial segment also increases, but it was not clear which ion channels were responsible for this entry of calcium. Channels that are selective for calcium ions are the prime candidates for this process. However, research in squid nerve cells gave rise to an unexpected idea by suggesting that sodium channels may not exclusively let in sodium but also allow some calcium ions to pass through. Hanemaaijer, Popovic et al. therefore wanted to test the routes that calcium ions take and see whether the sodium channels in mammalian nerve cells are also permeable to calcium.
Experiments using fluorescent dyes to track the concentration of calcium in rat and human nerve cells showed that calcium ions accumulated at the axon initial segment when action potentials fired. Most of this increase in calcium could be stopped by treating the neurons with a toxin that prevents sodium channels from opening. Electrical manipulations of the cells revealed that, in this context, the calcium ions were effectively behaving like sodium ions. Human kidney cells were then engineered to produce the sodium channel protein. This confirmed that calcium and sodium ions were indeed both passing through the same channel.
These results shed new light on the relationship between calcium ions and sodium channels within the mammalian nervous system and that this interplay occurs at the axon initial segment of the cell. Genetic mutations that ‘nudge’ sodium channels towards favoring calcium entry are also found in patients with autism spectrum disorders, and so this new finding may contribute to our understanding of these conditions.
Introduction
Ca2+ ions crossing the neuronal plasma membrane are critically involved in depolarization and distribute in the cytosol in spatial microdomains and organelles to regulate a wide range of processes ranging from gene expression to fast transmitter release (Berridge, 2006; Neher and Sakaba, 2008). In axons, voltage-gated Ca2+ (CaV) channels at presynaptic terminals open in response to a single action potential (AP), raising intracellular Ca2+ concentrations ([Ca2+]i) in nanodomains from ~50 nM up to ~10 µM to increase transmitter vesicle release rates by the power of ~4 (Helmchen et al., 1997; Schneggenburger and Neher, 2000). In response to APs, large and local [Ca2+]i transients are typically also observed in the axon initial segment (AIS) and nodes of Ranvier (Callewaert et al., 1996; Bender and Trussell, 2009; Yu et al., 2010; Gründemann and Clark, 2015; Zhang and David, 2016; Clarkson et al., 2017). At these sites, Ca2+ currents have been implicated in AP initiation and propagation by a local depolarizing action of the inward current or by activating the large conductance, Ca2+- and voltage-dependent K+ (BKCa) channels modulating burst firing probability and limiting frequency-dependent AP failure rates (Bender and Trussell, 2009; Yu et al., 2010; Hirono et al., 2015). The CaV channel subtypes identified in axons are both cell type- and species-dependent and include the T-, P/Q- or N-type CaV channels (Callewaert et al., 1996; Bender and Trussell, 2009; Yu et al., 2010; Gründemann and Clark, 2015; Zhang and David, 2016). At the AIS in particular the T-type Ca2+ channel mediates AP-dependent Ca2+ influx (Bender and Trussell, 2009; Martinello et al., 2015; Fukaya et al., 2018; Jin et al., 2019). However, in the prefrontal cortical pyramidal neuron AIS about 70% of the AP-evoked [Ca2+]i remains following pharmacological block of T-type CaV channels (Clarkson et al., 2017). Furthermore, evidence for a clustering of T-type CaV channels at the AIS is ambiguous and immunofluorescence or immuno-gold labeling studies show a density which is comparable to somatodendritic or spine compartments (McKay et al., 2006; Martinello et al., 2015).
Several other mechanisms may contribute to axoplasmic [Ca2+]i elevations in the AIS. Firstly, Ca2+ levels could rise due to Ca2+-induced Ca2+ release from intracellular sources such as the endoplasmic reticulum (ER). Most AISs contain ER cisternae organelles consisting of stacks of membranes expressing the store-operated ryanodine receptor (RyR), inositol 1,4,5-triphosphate receptor 1(IP3R1) and sarcoplasmic ER Ca2+ ATPase (SERCA) pumps (Benedeczky et al., 1994; King et al., 2014; Antón-Fernández et al., 2015). The coupling of transmembrane Ca2+ entry with intracellular store release may generate a local activity-dependent rise of [Ca2+]i. However, a contribution of ER stores to AIS [Ca2+]i remains to be directly demonstrated. Secondly, near the peak of the AP the electrogenic Na+-Ca2+ exchanger (NCX) reverses direction and imports Ca2+. A reverse mode of operation has not only been implicated in pathological [Ca2+]i elevations in axons during hypoxia and injury (Stys et al., 1991; Iwata et al., 2004), but also occurs during trains of APs in nodes and neighboring internodes (Zhang and David, 2016). Finally, one alternative pathway that has yet to be directly examined in mammalian cortical axons involves the voltage-gated Na+ (NaV) channels. Studies in the squid giant axon combining electrophysiological recordings with Ca2+ imaging showed that an early component of depolarization-induced Ca2+ entry is tetrodotoxin (TTX)-sensitive (Baker et al., 1971; Meves and Vogel, 1973; Brown et al., 1975). Voltage-clamp recordings from axons and perfusing distinct ionic solutions provided a quantitative estimate that NaV channels may pass divalent Ca2+ ions with permeability ratios (PCa/PNa) up to 0.10 (Hille, 1972; Meves and Vogel, 1973). Ca2+ permeability of NaV channels has also been shown in cardiac cells and hippocampal neurons (Akaike and Takahashi, 1992; Aggarwal et al., 1997; Santana et al., 1998) but whether this extends to the cortical axons remains to be examined.
Here, using wide-field Ca2+ imaging with a high-speed CCD camera enabling detection of [Ca2+]i changes at high sensitivity and high temporal resolution (Jaafari et al., 2014; Ait Ouares et al., 2016), we explored the various pathways of Ca2+ entry in axons of rat thick-tufted neocortical layer 5 (L5) pyramidal neurons. We found that during subthreshold depolarizations [Ca2+]i transients were highly compartmentalized to the AIS and nodes of Ranvier. While these transients were amplified by ER store release, the trigger was only modestly accounted for by CaV channels. The largest fraction of activity-dependent [Ca2+]i was TTX-sensitive and overlapped with the rapid gating of NaV channels. Experiments in HEK-293 cells transfected with the human NaV1.2 channel confirmed that TTX-sensitive Na+ currents were sufficient to generate [Ca2+]i elevations. Together, the data suggest that [Ca2+]i dynamics in the mammalian AIS are predominantly mediated by a rapid Ca2+ entry through NaV channels.
Results
Activity-dependent compartmentalized Ca2+ entry in layer five axons
Thick-tufted L5 pyramidal neurons, also called L5B or pyramidal tract neurons, are the largest pyramidal neurons in the cortex and integrate synaptic inputs from all cortical layers, playing a central role in cognitive tasks including perception (Groh et al., 2010; Ramaswamy and Markram, 2015; Takahashi et al., 2016). Their large axons (~1.5 µm in diameter) send long-range output projections to the thalamus, striatum and spinal cord, but within the cortex branch sparsely and have a trajectory perpendicular to the pia providing an excellent anatomical arrangement to image and record from. To optically record the spatial profile of axonal [Ca2+]i we made somatic whole-cell patch-clamp recordings from neurons filled with the high-affinity Ca2+ indicator Oregon Green BAPTA 1 (OGB-1, 100 µM) and imaged epifluorescence signals along the proximal region of the main axon (Figure 1). We first used subthreshold depolarizations evoked by artificial excitatory postsynaptic potentials (a-EPSPs, 100 Hz, peak depolarization 17.0 ± 0.6 mV, n = 15; Figure 1a). Examination of the spatial profile revealed that Ca2+ signals were observed in the AIS and hot spots separated with regular distances along the axon (locations 2, 4 and 6; Figure 1a). In order to examine whether the [Ca2+]i hot spots corresponded to nodes of Ranvier, we post-hoc stained for βIV-spectrin and biocytin, and found indeed overlap between subthreshold [Ca2+]i rise and spectrin-enriched sites (Figure 1b). In the same cells we examined the spatial profile of Δ[Ca2+]i in response to a single AP evoked with a brief square current injection (Figure 1c). As expected from back- and forward-propagating APs with much higher depolarizations (~100 mV), large [Ca2+]i transients were observed widespread throughout all axonal and somatodendritic domains. Population analysis showed that AP-induced [Ca2+]i transients were similar between AIS and nodes (one-way ANOVA followed by Tukey’s multiple comparison test, p<0.0001, differences between all groups were significant (p<0.05) except between AIS and node (p=0.13) and between internode (IN) and dendrites (p=0.85); Figure 1c,d). Interestingly, also during a-EPSPs the [Ca2+]i transients in the AIS and the first nodes were highly comparable, while [Ca2+]i signals in the internodal and dendritic domains were an order of magnitude smaller (one-way ANOVA with Tukey’s multiple comparison test, p<0.0001, differences between all groups were significant (p<0.0001), except between AIS and node (p=0.38) and IN and Dend (p=0.97); Figure 1a,d). Similar experiments in L5 neocortical pyramidal neurons in slices from human temporal cortex also revealed a-EPSP evoked Δ[Ca2+]i in the AIS, but not in the dendrite, suggesting that subthreshold sensitive [Ca2+]i transients are conserved across mammalian species (Figure 1—figure supplement 1). Together, these results show that activity-dependent [Ca2+]i transients are spatiotemporally compartmentalized and Ca2+ entry dynamics are similar in the axoplasm of the AIS and nodes.
Giant saccular organelle amplifies activity-dependent [Ca2+]i transients in the AIS
The thick-tufted L5 pyramidal neuron AIS contains a unique variant of cisternal organelle characterized by a continuous tubular organization of smooth ER, called the giant saccular organelle (Antón-Fernández et al., 2015). Cisternal organelles with smooth ER express synaptopodin (synpo), RyR, the IP3 receptor 1, and SERCA that are thought to contribute to Ca2+ release, buffering and storage (Bas Orth et al., 2007; King et al., 2014). We hypothesized that these organelles could generate Ca2+-induced Ca2+-release, thereby contributing to domain-selective activity-dependent [Ca2+]i transients (Figure 1). Triple immunostaining for synpo, Ankyrin G and biocytin confirmed that the cisternal organelle was present along the entire axis of the AIS and spatially overlapped with the subthreshold-evoked [Ca2+]i transients (n = 19; Figure 2a,b). However, while subthreshold depolarization-induced Ca2+ transients were present in the nodes, synaptopodin expression was not detected (n = 10 nodes; Figures 1 and 2a). To experimentally test whether AIS Ca2+-store release contributes to activity-dependent [Ca2+]i transients we performed experiments with standard intracellular solution and subsequently re-patched the same cell with a solution containing ryanodine (200 µM, blocking RyR-mediated Ca2+ release) and heparin (5 mg/ml, competitively inhibiting IP3-evoked Ca2+ release; Figure 2c). Blocking Ca2+ release significantly lowered ΔF/F Ca2+ peak transients in the AIS, both for the subthreshold- and AP-evoked [Ca2+]i changes (a-EPSP, 53.2% reduction, p=0.006; AP, 34.3% reduction, p=0.02, one-tailed ratio paired t-tests, n = 5; Figure 2d). Consistent with the AIS-specific location of the giant saccular organelle, store blockers had no effect on AP-evoked Δ[Ca2+]i in the basal dendrite (p=0.48, n = 3; Figure 2d). Furthermore, since the stores contribute significantly to AIS Ca2+ levels, blocking store release could act as a low-pass filter for Ca2+ level kinetics, reducing rise and decay times. However, blocking Ca2+-store release did not alter the rise- or decay time kinetics in the AIS (τact, p=0.52; τde-act, p=0.18, two-tailed paired t-tests, n = 5; Figure 2e). These data suggest that the giant saccular organelle amplifies activity-dependent [Ca2+]i changes selectively in the AIS.
CaV channels and NCX have a limited role in activity-dependent Ca2+ entry at the AIS
Ca2+ release from internal stores is likely triggered by Ca2+ entry via neuronal voltage-dependent plasmalemmal routes. To test whether AIS [Ca2+]i changes require Ca2+ from the extracellular space, we bath applied 2.5 mM of the Ca2+ chelator EGTA which effectively lowered the extracellular Ca2+ concentration ([Ca2+]o) from 2 mM to ~437 nM, thereby reducing the driving force for Ca2+ (see Materials and methods). Ca2+ imaging at the AIS (OGB-1, 100 µM) showed that EGTA almost fully abolished the subthreshold-evoked Δ[Ca2+]i (90.7% reduction, one-tailed ratio paired t-test, p=0.0031, n = 4; Figure 3a,b). Similarly, the AP-generated Δ[Ca2+]i was almost extinguished after bath application of EGTA (92.8% reduction, p=0.0011, n = 4; Figure 3a,b).
Next, we hypothesized that the transmembrane pathway for Ca2+ entry in the AIS during subthreshold stimuli could be explained by the low-voltage gated CaV channels (T- and R-type). However, bath application of the highly selective T-type (CaV3.1–3) blocker TTA-P2 (1 µM, Choe et al., 2011) or nickel (Ni2+, 100 µM) did not significantly reduce Ca2+ signals (one-tailed ratio paired t-tests; TTA-P2, p=0.17, n = 4; Ni2+, p=0.063, n = 5; Figure 3b,c). We next blocked R-type CaV channels, by puffing SNX-482 (500 nM) locally to the AIS, but this did not lead to a reduction in subthreshold [Ca2+]i rise either (SNX-482, p=0.11, n = 3). Furthermore, consistent with their more depolarized voltage range of activation, the L-type channels did not affect subthreshold Δ[Ca2+]i (20 µM isradipine, p=0.14; 10 µM nimodipine, p=0.41, both n = 4; Figure 3b) and the block of N-type and P/Q-type channels, by local application of ω-conotoxin MVIIC (2 µM) to the AIS, also failed to reduce Ca2+ signals (p=0.42, n = 5; Figure 3b). Furthermore, a combined block of T- and L-type channels did not affect the peak ΔF/F in the AIS (TTA-P2 and isradipine, p=0.12, n = 3; Figure 3b).
Although application of the T-type blockers TTA-P2 and Ni2+ was ineffective to block subthreshold [Ca2+]i rise, in the same neurons it did reduce the peak ΔF/F evoked by a single AP by almost 20% (TTA-P2, 18.7%, p=0.021, n = 4; Ni2+, 19.7% block, p=0.013, n = 5; Figure 3b,c). In addition, isradipine reduced the AP-evoked Δ[Ca2+]i in the AIS by 14.9% (isradipine, p=0.0070, n = 4) and the alternative L-type blocker nimodipine showed a non-significant blocking trend (nimodipine, p=0.060, n = 4; Figure 3b). A combined application of T- and L-type channel blockers (1 µM TTA-P2 and 20 µM isradipine) caused a 27.2% reduction of peak ΔF/F, showing a sublinear summation of two blocking agents (TTA-P2 and isradipine, p=0.0071, n = 5; Figure 3b). In contrast, the R-type CaV channel blocker SNX-482 (500 nM) did not reduce the AP-evoked Δ[Ca2+]i (SNX-482, p=0.29, n = 3; Figure 3b). Local application of the ω-conotoxin MVIIC (2 µM) showed a non-significant trend to block the peak ΔF/F (6.8%, p=0.064, n = 5; Figure 3b). As a positive control experiment, we imaged a collateral bouton of the same neuron and used local application of ω-conotoxin MVIIC which almost completely abolished the peak ΔF/F by 91.5%, consistent with the presence of N- and P/Q-type CaV channel subtypes in presynaptic terminals (p=0.021, n = 3; Figure 3b). Finally, given the unexpected remaining [Ca2+]i transients in the AIS in the presence of CaV channel blockers, we hypothesized that NCX may contribute to activity-dependent [Ca2+]i increase in the AIS. At the resting membrane potential NCX exports Ca2+ from the cytoplasm to maintain [Ca2+]i near ~50 nM, however Na+ entry may promote instantaneous Ca2+ influx by a reverse mode of operation (Yu and Choi, 1997; Stys and LoPachin, 1997; Figure 3b). To examine its contribution, we pharmacologically blocked NCX by combined bath application of KB-R7943 (20 µM) and SN-6 (10 µM). The results showed, however, no change in the subthreshold nor AP-evoked Δ[Ca2+]i (two-tailed ratio paired t-tests, a-EPSP, p=0.16; AP, p=0.13, n = 5, respectively; Figure 3b).
In summary, these data show that while transmembrane Ca2+ influx from the extracellular space is required for activity-evoked Δ[Ca2+]i, none of the CaV channels played a role in the subthreshold depolarization, whereas T- and L-type CaV channels partially contributed to the AP-evoked influx.
Subthreshold- and AP-evoked Ca2+ entry at the AIS requires TTX-sensitive channels
What could be the source of the remaining component of Ca2+ influx at the AIS? Both in hippocampal neurons and heart muscle cells, Ca2+ currents have been described which are not blocked by Ni2+ nor by other known CaV channel antagonists, but instead are sensitive to the highly selective NaV channel blocker (TTX), and therefore called ICa(TTX) (Akaike and Takahashi, 1992; Aggarwal et al., 1997). In Na+-free extracellular solution ICa(TTX) resembles the Na+ current and activates at potentials as negative as –70 mV while peaking at –30 mV (Akaike and Takahashi, 1992). To examine the presence of ICa(TTX) in L5 axons we took advantage of the low-affinity indicator Oregon Green BAPTA 5N (OGB-5N, 1 mM; Figure 4a), which gives smaller fluorescent signals but is linear over a wider range of [Ca2+]i compared to OGB-1 (Kd 20 µM vs. 170 nM, respectively). We used the voltage-clamp configuration and injected depolarizing ramps of 50 ms with increasing slopes (from 0.0 to 0.55 mV ms–1) with a maximum peak at ~95% of the AP threshold (Figure 4b), thereby studying the same depolarization range and duration as the a-EPSPs used in Figures 1–3. The results showed that Ca2+ influx was strongly compartmentalized to the AIS and nodal axolemma (Figure 4a) consistent with the a-EPSP evoked OGB-1 transients (Figure 1). Remarkably, bath application of TTX almost completely abolished [Ca2+]i elevations, even at depolarizations above the AP threshold (at –54.5 mV, 92.8% block, Cohen’s d: 1.88, two-way ANOVA, p<0.0001, n = 4; Figure 4b,c). As an alternative to TTX we used the quaternary amine NaV channel inhibitor QX-314, which plugs the open state of the NaV channel from the internal side. Similar to TTX, with 6 mM QX-314 added to the pipette solution voltage ramps did not evoke Ca2+ transients (at –54.5 mV, 94.8% block, Cohen’s d: 1.92, control vs. QX-314, p<0.0001, n = 4, TTX vs. QX-314, p=0.97, n = 4; Figure 4c). Although QX-314 at this concentration has been reported to also block voltage-gated Ca2+ currents (Talbot and Sayer, 1996), subthreshold-evoked [Ca2+]i was not mediated by CaV channels (Figure 3b). The near complete block by two distinct NaV channel blockers therefore indicates an important role of NaV channels in mediating subthreshold axonal Ca2+ influx.
We next investigated whether NaV channels also contribute to AP-evoked Δ[Ca2+]i (Figure 5a). To dissociate a putative role of NaV channels to pass Ca2+ ions from generating the AP depolarization of ~100 mV we first imaged Ca2+ at the AIS in current-clamp, subsequently applied 1 µM TTX and imaged Ca2+ transients evoked by the recorded AP waveform injected as a voltage command (‘AP-clamp’). The results showed a near complete abolishment of Δ[Ca2+]i in the presence of TTX (one-way ANOVA with Tukey’s multiple comparisons test, CC vs. VC, 89.5% reduction, p<0.0001, n = 7; Figure 5a,b). However, this [Ca2+]i peak amplitude reduction could also be due to an incomplete voltage- and space-clamp of the AIS for fast voltage transients. The small diameter of the axon has a high axial resistance, acting as a low-pass filter for the antidromic AIS action potentials (Hamada et al., 2016) which also may attenuate the orthodromic voltage spread into the axon. To examine the possibility that axonal APs are attenuated in the somatic AP clamp configuration, we optically recorded JPW3028, a fast fluorescent voltage indicator that remains stable over long recording periods and is highly linear over a large voltage range (~250 mV, Figure 5—figure supplement 1). Consistent with the voltage loss, when we injected the AP-clamp in the soma in the presence of TTX and optically recorded Vm in the AIS, we observed a significant ~2 fold reduction in the AP peak amplitude (one-way ANOVA with Tukey’s multiple comparisons test, VC vs. CC, p=0.014; Figure 5a,c). To restore peak depolarization in the presence of TTX and reliably compare the Ca2+ transients evoked by equal depolarization, we doubled the amplitude of the somatic AP-clamp (VC ×2). With this protocol both the peak depolarization and AP half-width in the AIS were indistinguishable from the control APs (peak JPW, VC ×2 vs. CC, one-way ANOVA with Tukey’s multiple comparisons test, p=0.75, n = 4, Figure 5c, half-width in JPW, VC ×2 vs. CC, one-way ANOVA with Tukey’s multiple comparisons test, p=0.36, n = 4, not shown). TTX blocked 65.5% of the AP-evoked Δ[Ca2+]i (peak OGB-5N, Cohen’s d: 5.49, VC ×2 vs. CC, p<0.0001, n = 8; Figure 5b). Taken together, these data suggest that a large fraction of both subthreshold-depolarization and AP-evoked Ca2+ ions enter the axoplasm through TTX-sensitive channels at the AIS.
Axonal [Ca2+]i follows gating kinetics of NaV channels
Whether ICa(TTX) is carried by a specific TTX-sensitive CaV channel or reflects Ca2+ permeating directly through the NaV channel remains debated (Santana et al., 1998; Cruz et al., 1999; Heubach et al., 2000; Chen-Izu et al., 2001; Guatimosim et al., 2001). We hypothesized that if Ca2+ ions enter the AIS cytoplasm by flowing through NaV channels, ICa(TTX) should reflect the time course of INa. To measure submillisecond rapid events with fluorescence in the small axon (diameter ~1.5 µm), we optimized multiple imaging parameters enabling the acquisition of fluorescence at 20 kHz (see Materials and methods). Using a Na+-sensitive indicator (sodium-binding benzofuran isophthalate, SBFI, 1 mM) in combination with OGB-5N (1 mM) showed that the two indicators were indistinguishable in their rising phase during an AP, suggesting that Ca2+ entry at the AIS may be as rapid as Na+ entry (Figure 6—figure supplement 1).
In order to quantify the kinetics of INa and ICa more directly, we next used a voltage-clamp approach. Near ~20°C and –40 mV, NaV channels open at least one order of magnitude faster compared to the T-type CaV channels (~200 µs [Schmidt-Hieber and Bischofberger, 2010] vs. ~5 ms [Perez-Reyes et al., 1998], respectively) which may be sufficiently different to compare against the kinetics of optically recorded [Ca2+]i at the AIS. To determine the specific activation kinetics of INa and T-type ICa (ICaT) in L5 pyramidal neurons we measured total inward current (INa + ICa,) by depolarizing the soma with a step to –35 mV and pharmacologically isolated Na+ and Ca2+ current components by bath application of 1 µM TTX or 100 µM Ni2+, respectively (Figure 6a). The activation time constant of the total inward current was identical to INa (single exponential fit τtotal = 438.2 µs vs. τNa = 440.3 µs, one-way ANOVA with Tukey’s multiple comparison test, p>0.999, n = 6), whereas the total current was substantially faster in comparison to ICaT (τtotal = 438.2 µs vs. τCaT = 4.8 ms, p<0.0001, n = 5; Figure 6a). The initial fraction of the inward current was thus primarily generated by INa. The large difference in gating could provide a temporal window to distinguish Ca2+ entry via NaV channels or T-type CaV channels. Theoretical and experimental work show that low-affinity Ca2+ indicators, like OGB-5N, are capable of tracking rapidly activating Ca2+ currents when imaged at high speed: the first time derivative of ΔF/F (dΔF/F dt–1) overlaps with electrically recorded ICa, providing a mean to optically resolve the time course of ICa (Sabatini and Regehr, 1998; Jaafari et al., 2014; Ait Ouares et al., 2016).
Imaging OGB-5N (1 mM) at 20 kHz in the AIS we observed that ΔF/F comprised of two separate time courses, a fast initial rise followed by a slower rising phase (Figure 6b). Both components were almost completely abolished by TTX, leaving only a small transient reflecting putatively ICaT (n = 6). We quantitatively compared the activation time constants of INa, ICa_opt (dΔF/F dt–1) and ICaT by resampling the electrically recorded INa and ICaT to 20 kHz and filtering both electrical and optical traces identically (see Materials and methods). Multiple hallmarks of INa matched with ICa_opt: both traces showed a rapid inward component, followed by inactivation and a persistent component (Figure 6b,c). In comparison, ICaT lacked both the rapid activation and inactivation time constants (Figure 6a,c). Given the lower signal-to-noise ratio in the optical traces we fitted Boltzmann sigmoid functions to the rising phase to compare the slopes of the optically and electrically recorded currents (Figure 6c). The average slope of ICa_opt was significantly faster compared to the activation of ICaT (~500 µs vs. ~2.5 ms, respectively, one-way ANOVA with Tukey’s multiple comparison test, p<0.0001) and slower compared to INa (~500 µs vs. ~200 µs, ICa_opt vs. INap<0.0001; Figure 6c,d). The small difference between INa and ICa_opt may be explained by the equilibration time of OGB-5N (~200 µs) (Ait Ouares et al., 2016), local differences between NaV channels in the soma and AIS or the presence of Ca2+-store release in the AIS. Together, the findings indicate that the current mediating [Ca2+]i at the AIS resembles NaV channel kinetics.
Calcium influx through NaV1.2 channels
The results suggest that Ca2+ ions could enter the cytoplasm by permeation through the NaV channel pore. Previous studies showed Ca2+ influx through the cardiac NaV1.5 channel (Cruz et al., 1999; Guatimosim et al., 2001). To examine whether NaV channel isoforms of the axon initial segment also enable Ca2+ influx we performed experiments in HEK-293 cells which were transfected with the human gene SCN2A encoding NaV1.2 channel with auxiliary β1 and β2 subunits and EGFP tag (Materials and methods; Figure 7a). Whole-cell recording revealed Na+ currents in EGFP+ cells but not in non-transfected cells (average peak current density –115.7 ± 28.4 pA/pF, n = 10 vs. –3.7 ± 1.9 pA/pF at –20 mV, n = 5, respectively; Figure 7b). The inward currents were completely abolished by 1 µM TTX (96.1 ± 1.3%, n = 7, one-tailed Wilcoxon matched-pairs signed rank test, p=0.0078, Figure 7c) and the voltage-dependence of activation and inactivation revealed midpoints at –25.4 ± 2.1 mV and –74.1 ± 3.9 mV, respectively (n = 10, Figure 7d,e), consistent with previous work (Ben-Shalom et al., 2017), indicating a highly selective expression of NaV1.2 channels. Next, we filled the transfected cells with 100 µM OGB-1 and imaged the fluorescence changes in response to a train of depolarizing pulses (200 Hz for 1 s, –120 to –30 mV steps; Figure 7f). We observed an increase in ΔF/F in every EGFP+ cell, indicating an influx of Ca2+ (average peak 0.46 ± 0.18% ΔF/F, range: 0.06–1.4% ΔF/F, n = 7; Figure 7f–h). To test whether the [Ca2+]i increase required NaV channel opening we bath applied TTX (1 µM), revealing a significant decrease in the peak ΔF/F (92.1 ± 3.8% reduction, one-tailed Mann Whitney test, p=0.012, n = 4; Figure 7h). The results indicate that molecular expression and opening of NaV1.2 channels suffices to mediate transmembrane Ca2+ influx.
Estimating Ca2+ conductivity of NaV channels with computational modeling
Our findings are in agreement with the depolarization-induced Ca2+ entry in the squid axon which is tetrodotoxin (TTX)-sensitive and reflected a 1% conductivity of NaV channels for Ca2+ ions (Baker et al., 1971). To estimate the conductivity ratios (gCa/gNa) in L5 axons we performed computational simulations. Ca2+ entry through NaV channels was implemented by adding an ohmic Ca2+ ion mechanism into a mathematical 8-state NaV channel model that calculated the current carried by Ca2+ (ICa(Na)) and Na+ (INa) (see Materials and methods). A single compartment containing INa and ICa(Na) together with high voltage-gated and T type-gated CaV channel models (ICaH and ICaT, respectively) showed that with an axonal AP waveform ICa(Na) is activated during the first microseconds of AP onset, rapidly inactivates and is temporally separated from ICaH and ICaT (Figure 8—figure supplement 1). Next, to estimate the gCa/gNa we made a multicompartmental model of a L5 pyramidal neuron (Figure 8a, including detailed reconstructions of the AIS and nodal domains (see Figure 2a within Hamada et al., 2016). Based on multiple experimentally recorded parameters we constrained the model AP and found that a peak conductance density of NaV channels of 16,000 and 850 pS µm–2 in the AIS and soma, respectively, reproduced the recorded AP and matched with AP-evoked Δ[Na+]i imaged in the AIS (see Materials and methods and Figure 8a). Subsequently, [Ca2+]i was simulated based on mathematical equations representing Ca2+ diffusion and extrusion, endogenous stationary Ca2+ buffers (taken together as κs) and was supplemented with the buffering capacities of the specific Ca2+ indicators (Fink et al., 2000) (see Materials and methods). The Ca2+ extrusion threshold and rates were adjusted to approximate the experimentally imaged peak and decay time course of measured OGB-5N ΔF/F in the AIS (Figure 8—figure supplement 2). To determine the absolute rise in [Ca2+]i produced exclusively by NaV channels, we performed additional experiments in which we imaged [Ca2+]i while blocking CaV channels that contributed to AP-evoked Δ[Ca2+]i: T- and L-type calcium channels (TTA-P2 and isradipine, respectively, see Figure 3b). Using calibrated ratiometric bis-Fura-2 (200 μM) imaging, we found that during 1 AP, Ca2+ entry though NaV channels induces a peak Δ[Ca2+]i of 55.6 nM (n = 4; Figure 8b, Figure 8—figure supplement 2). Since ~35% of AP-evoked Δ[Ca2+]i is caused by internal store amplification (Figure 2d) ~36 nM is mediated by transmembrane Ca2+ entry via NaV channels (Figure 8b). We subsequently simulated these experiments in the multicompartmental model by removing CaV channels and including the buffering properties of 200 µM bis-Fura-2. Varying endogenous buffering (κs) between 1 and 100 we updated gCa/gNa to obtain a 36 nM rise of free [Ca2+]i at the AIS. A κs of ~100 corresponds to dendritic buffering capacities (Cornelisse et al., 2007), while axonal buffering capacities are reported to be lower (10–40, Klingauf and Neher, 1997; Jackson and Redman, 2003; Delvendahl et al., 2015). When changing κs between 10–40 the gCa/gNa ratio was 0.38% (κs = 10: 0.37%, κs = 40: 0.39%; Figure 8c).
Spatiotemporal distribution of Ca2+ entry routes under physiological conditions
Using the 0.38% conductivity ratio we next evaluated how Ca2+ currents through NaV and CaV channels spatiotemporally varied across the neuronal compartments without the buffering capacities of externally applied Ca2+ dyes (Figure 9a,b and Figure 8—figure supplement 2). The simulations showed that the Δ[Ca2+]i from one AP reached a peak concentration of ~800 nM in the AIS (Figure 9a,b). Due to the high density of NaV channels in the AIS they contribute to the majority of Δ[Ca2+]i and cause a rise of [Ca2+]i within submillisecond from the start of the AP (450 nM within <150 µs from AP threshold, red arrow in Figure 9b). These results are likely to provide an underestimate of the total Δ[Ca2+]i since in the model Ca2+ release from giant saccular organelle was not simulated, which would result in a total AP-evoked Δ[Ca2+]i of ~1.2 µM. In the basal dendritic branches the AP has a slower rise time and broader half-width, causing dendritic [Ca2+]i to accumulate slower and to higher concentrations, consistent with our experimental findings (Figure 1). Because the dendritic NaV channel density is substantially lower, their contribution to the total [Ca2+]i is negligible. The distinct contribution of NaV and CaV channels to [Ca2+]i is clearly visible when comparing the different Ca2+ currents in the AIS, showing that the majority of the total ICa during an AP is carried by ICa(Na) (Figure 9c). Simulations predict that the ICa(Na) rapidly inactivates during the AP while ICa activates more slowly and has an incomplete inactivation during the AP repolarization, likely becoming the dominant contribution to [Ca2+]i during the afterdepolarization and high-frequency spike generation (Figure 9c and Figure 8—figure supplement 1).
Our experiments and computational simulations show that [Ca2+]i changes mediated by CaV or NaV channels act at distinct spatiotemporal scales. To experimentally test the differential impact of Ca2+ influx on the AP waveform, we analyzed the somatically recorded APs when using distinct blockers most effective in modulating [Ca2+]i in the AIS (Figure 3). Blocking both T- and L-type CaV channels, contributing to ~27% to AP-evoked Ca2+ at the AIS (Figure 3), significantly reduced the afterdepolarization and showed a trend to reduce the afterhyperpolarization (ADP, p=0.039 and AHP, p=0.055, respectively, two-tailed paired t-tests, n = 5), without affecting other AP properties (p>0.30; Figure 9d and Table 1). In contrast, when lowering [Ca2+]o with EGTA, which abolished all AIS Ca2+ influx (Figure 3), the AP half-width significantly increased and the AHP was reduced (AP half-width, p=0.035 and AHP, p=0.043, two-tailed paired t-tests, n = 4; Figure 9d, Table 1). These results are consistent with the temporal differences in AP-evoked AIS [Ca2+]i and suggest that NaV-mediated Ca2+ entry may act to open BKCa channels (Figure 9e), thereby driving K+ efflux and facilitating axonal AP repolarization.
Discussion
In the present study we identified Ca2+ permeation through NaV channels as a source for activity-dependent Ca2+ entry in mammalian axons. The findings were supported by independent and converging lines of evidence, ranging from anatomical compartmentalized [Ca2+]i transients at sites with high NaV channel densities, a pharmacological block by TTX, an overlap of optically recorded ICa with NaV channel gating and molecular evidence for Ca2+ influx mediated by the NaV1.2 channel. In axonal domains NaV channels are thus not only involved in electrically generating the upstroke of the action potential, but also contribute to cytoplasmic Ca2+ signaling.
Ca2+ entry in the L5 pyramidal neuron AIS was in part mediated by T- and L-type CaV channels (Figure 3 and Figure 5) in keeping with previous studies showing that CaV channel subtypes mediate activity-dependent Ca2+ changes in both central- and peripheral nervous system axons (Callewaert et al., 1996; Bender and Trussell, 2009; Yu et al., 2010; Gründemann and Clark, 2015; Zhang and David, 2016; Clarkson et al., 2017). However, when quantifying the specific fraction of block CaV channels explained only ~35% of the total Ca2+ entry during a single AP (Figure 3). These results are consistent with recent 2-photon Ca2+ imaging from the prefrontal cortical pyramidal neuron AIS, showing that ~70% of the [Ca2+]i transients evoked by a train of three APs remains in the presence of T-type CaV channel block (Clarkson et al., 2017). Here, we found that for APs the remaining ~65% of [Ca2+]i increase is actually TTX-sensitive and this accounted even for >90% of the subthreshold-induced [Ca2+]i changes (Figure 4 and Figure 5). That Ca2+ enters through NaV channels builds on landmark studies showing that the initial component of depolarization-induced Ca2+ entry in the squid giant axon is tetrodotoxin (TTX)-sensitive (Baker et al., 1971; Meves and Vogel, 1973; Brown et al., 1975). Squid axons even generate rapid spikes in the sole presence of Ca2+ ions (Watanabe et al., 1967). A TTX-sensitive Ca2+ current has also been identified in hippocampal neurons (Akaike and Takahashi, 1992) and more extensively investigated in cardiac myocytes (Aggarwal et al., 1997; Santana et al., 1998; Cruz et al., 1999; Heubach et al., 2000; Chen-Izu et al., 2001). The incomplete selectivity of toxins and blockers for Ca2+ channels continued, however, to cast doubt about the precise identity of the TTX-sensitive Ca2+ current (Cruz et al., 1999; Chen-Izu et al., 2001; Sun et al., 2008). Indeed, an alternative explanation for some of the present results is that TTX blocks a CaV channel subtype. This has been reported at very high concentrations (30 µM) for CaV3.3 (Sun et al., 2008), which is higher than what we used (1 µM). Although we cannot exclude the presence of a TTX-sensitive CaV channel at the AIS that is not blocked by one of the compounds used in the pharmacological screening (Figure 3), our optically recorded ICa provides biophysical evidence that the TTX-sensitive Ca2+ current follows the same rapid activation time course as the Na+ channel pore, incompatible with T-type channel kinetics (Figure 6). Importantly, in further support that NaV channels give rise to cytoplasmic Ca2+ changes, heterologous expression of α- and β-subunits of NaV1.2, known to be expressed in the rodent and human AIS (Garrido et al., 2003; Hu et al., 2009; Tian et al., 2014), showed that the channel proteins expressed in isolation were sufficient to mediate Ca2+ influx (Figure 7). These results are in support of the findings by Lederer and colleagues showing that heterologous expression of hNaV1.5 channel produces depolarization-evoked [Ca2+]i, if expressed with its β subunits (Cruz et al., 1999; Guatimosim et al., 2001).
From an evolutionary point of view, some Ca2+ permeation of NaV channels is not surprising. NaV channels evolved from the CaV channel superfamily and share molecular structure both in their pore sequence and intracellular regulatory domains (Zakon, 2012; Ben-Johny et al., 2014). Furthermore, the Born radii for Na+ and Ca2+ are comparable (1.68 and 1.73 Å, respectively) and Ca2+ ions are known to enter the NaV channel pore to block Na+ permeation in a concentration-dependent manner (Lewis, 1979; Armstrong and Cota, 1999). Interestingly, single residue mutations in the selectivity filter of NaV channels suffices to increase Ca2+ ion permeation (Heinemann et al., 1992; Naylor et al., 2016). Our calculations indicate that the conductivity of the NaV channel for Ca2+ is ~0.4% (Figure 7). If we assume there exists proportionality between permeability and conduction we can apply the equation of permeability ~ conductance/concentration × valency2 (Baker et al., 1971; Meves and Vogel, 1973). With our extracellular solutions [Ca2+]o/[Na+]o being 0.0148 and the valence ratio (Ca2+/Na+)2 being four we can calculate that NaV channels in mammalian axons have a PCa/PNa ratio of 0.06. Notably, the value is in range of direct recordings for PCa/PNa in the sciatic nerve and squid axons (0.10 and 0.14, Hille, 1972; Meves and Vogel, 1973) as well as recordings from NaV1.5 channels revealing PCa/PNa ratios of ~0.04 (Cruz et al., 1999). Such permeability is orders of magnitude lower than PCa/PNa ratios of acetylcholine receptors or NMDA receptors (1.0 and 17, respectively) (Lewis, 1979; Iino et al., 1997). An independence of Na+ and Ca2+ ions, modeled as ohmic conductances, will be a major simplification of the NaV channel under multi-ion conditions. Molecular dynamic studies of NaV channels showed that ionic interactions between Na+ and Ca2+ at the channel pore are complex (Corry, 2013; Boiteux et al., 2014; Naylor et al., 2016). The energy barrier of the selectivity filter strongly favors Na+ ions but can be flexible, changing in conformational states and consistent with modest Ca2+ permeation (Corry, 2013; Boiteux et al., 2014; Naylor et al., 2016). In future experiments it will be interesting to obtain more detailed permeability ratios (PCa/PNa) by recording changes in Erev with varying intra- and extracellular concentrations, fitting the data to mathematical solutions such as the electro-diffusion theory of Goldman-Hodgkin-Katz (GHK) extended with surface charge potentials (Lewis, 1979; Campbell et al., 1988), or using Eyring–Läuger theory based on individual ionic rate constants (Läuger, 1973). Such experiments would in particular be interesting for NaV1.6, the main isoform expressed in axonal domains (Lorincz and Nusser, 2010; Kole and Stuart, 2012).
Although the Ca2+ conductivity of the channels is small it achieves near–micromolar Ca2+ changes in axons as NaV channels are clustered to very high densities (~1000 channels µm–2) at the AIS and nodes of Ranvier (Neumcke and Stämpfli, 1982; Lorincz and Nusser, 2010; Kole and Stuart, 2012). Consistent with this idea, our imaging experiments showed that also nodes of Ranvier produced subthreshold-activated Ca2+ entry (Figure 1 and Figure 4), suggesting that NaV channel mediated Ca2+ entry could play similar roles in these domains. At the AIS, the rapid opening of high densities of NaV channels may further act as a trigger to amplify [Ca2+]i via the activation of ryanodine receptors, mediating ER store release of Ca2+ from the giant saccular organelle, which extends continuously along the AIS of thick-tufted L5 pyramidal neurons (Figure 2).
The rapid inactivation of NaV channels compared to the slow inactivation of CaV channels will lead to voltage- and time-dependent changes in the relative contribution of NaV and CaV channels to [Ca2+]i. In axons, a single AP will mostly lead to NaV-mediated Ca2+ entry while an increasing number of APs, or longer sustained depolarization, will lead to an accumulation of Ca2+ mediated byCaV channel activation. Indeed, Ca2+ entry via CaV channels has been identified as a major contributor to Ca2+ entry in sciatic nerve or Purkinje axons during trains of APs or prolonged depolarization for hundreds of milliseconds (Callewaert et al., 1996; Zhang and David, 2016). However, in vivo recordings from L5 pyramidal neurons show that they typically fire sparsely and on average ~1–4 Hz (de Kock et al., 2007) and the half-width of axonal APs is about ~300 µs (Kole et al., 2007). In this view, NaV-mediated Ca2+ entry may be the main source for activity-dependent [Ca2+]i in the excitable domains of axons under physiological conditions.
What could be the functional role of NaV-mediated Ca2+ entry in axon initial segments and nodes of Ranvier? One downstream target of submembranous axoplasmic [Ca2+]i may be regulation of NaV inactivation kinetics via their Ca2+/calmodulin domain at their C-terminus as has been demonstrated for multiple NaV subtypes, including NaV1.2 and NaV1.6 (Sarhan et al., 2012; Reddy Chichili et al., 2013; Ben-Johny et al., 2014; Wang et al., 2014). Another target of NaV-mediated Ca2+ could be to open axonal large-conductance BKCa channels. The BKCa channel opens with the cooperative action of membrane voltage and [Ca2+]i ≥10 µM to repolarize APs and shorten their duration (Berkefeld et al., 2010). Across cell types there are considerable variations in the magnitude and time course of BKCa currents due to differences in nanodomain coupling with CaV channel isoforms. BKCa channels are exclusively activated by P/Q-type CaV channels in cerebellar Purkinje neurons with short AP durations (Womack et al., 2004), but in rat chromaffin cells with wider APs, BKCa channels are coupled with the Q- and slower activating L-type CaV channels (Prakriya and Lingle, 1999). Also in L5 pyramidal neurons BKCa activation shortens the duration of somatic APs from ~1 ms to ~600 µs (Yu et al., 2010; Bock and Stuart, 2016; Roshchin et al., 2018). Considering the brief duration of axonal APs in the L5 pyramidal neurons (~300 µs, Kole et al., 2007) NaV channels may provide both a precisely-timed voltage-dependent activation, via Na+ current, as well as a [Ca2+]i rise within 150 µs (Figure 9), to rapidly open BKCa channels and shape axonal AP repolarization. In agreement with this conjecture our data show that T- and L-type CaV channels are too slow to mediate somatic AP repolarization, leaving open the possibility that a NaV-BKCa channel nanodomain coupling provides the required Ca2+ signal. Firm evidence for such interaction would require mutating the selectivity filter of NaV channels to abolish Ca2+ but not Na+ permeation. Recent two-photon Ca2+ uncaging experiments already showed that the [Ca2+]i rise at the first node of Ranvier in L5 axons opens nodal BKCa channels to shorten the AP duration and facilitate the generation of high firing rates in the proximal axon (Roshchin et al., 2018). Further downstream from the initiation site nodal BKCa channels in Purkinje axons play a role in augmenting the hyperpolarization following APs and facilitate recovery from NaV inactivation to prevent propagation failures (Hirono et al., 2015). The dual permeation of axonal NaV channels for Na+ and Ca2+ ions may thus serve a common function; mediating the rapid electrical upstroke of the AP and via Ca2+ signaling activating K+ efflux to recover from inactivation and accelerating NaV channel availability for the next AP, representing a fine-tuning specifically to the needs of axonal AP generation and conduction fidelity.
Materials and methods
Ethical approval
Request a detailed protocolAll animal experiments were performed in compliance with the European Communities Council Directive 2010/63/EU effective from 1 January 2013. They were evaluated and approved by the national CCD authority (license AVD8010020172426) and by the KNAW animal welfare and ethical guidelines and protocols (DEC NIN 14.49, DEC NIN 12.13, IvD NIN 17.21.01 and 17.21.03). Written informed consent was obtained from patients and all procedures on human tissue were performed with the approval of the Medical Ethical Committee of the Amsterdam UMC, location VuMC and in accordance with Dutch license procedures and the Declaration of Helsinki. All data were anonymized.
Tissue collection
Request a detailed protocolYoung-adult male Wistar rats (RjHan:WI) were used at an age between P21 and P35 (Charles River Laboratories and Janvier labs). Animals were deeply anaesthetized by 3% isoflurane inhalation, decapitated and 300 µm parasagittal slices containing the primary somatosensory cortex were cut with a Vibratome (1200S, Leica Microsystems B.V.) within ice-cold artificial cerebrospinal fluid (ACSF) of the following composition (in mM): 125 NaCl, 3 KCl, 25 glucose, 25 NaHCO3, 1.25 Na2H2PO4, 1 CaCl2, 6 MgCl2, saturated with 95% O2 and 5% CO2 (pH 7.4). Following a recovery period at 35°C for 35–45 min slices were stored at room temperature in the ACSF. Human slices were obtained from non-pathological cortex removed for the surgical treatment of deeper brain structures for mesial temporal lobe epilepsy. After resection, a block of the temporal lobe was placed within 30 s in ice-cold artificial cerebrospinal fluid (ACSF) slicing solution which contained in (mM): 110 choline chloride, 26 NaHCO3, 10 D-glucose, 11.6 sodium ascorbate, 7 MgCl2, 3.1 sodium pyruvate, 2.5 KCl, 1.25 NaH2PO4, and 0.5 CaCl2 (300 mOsm) and transported to the laboratory, as described in detail previously (Testa-Silva et al., 2014). Transition time between resection of the tissue and preparation of the slices was <15 min. Neocortical slices (~350 µm thickness) were cut in an ice-cold slicing solution, stored for 30 min at 34°C, and afterwards switched to room temperature in standard ACSF. Slices were subsequently transported (<15 min) towards the NIN (KNAW) in continuously carbogenated ACSF.
Cell lines
Request a detailed protocolHuman embryonic kidney 293 cells (HEK 293T/17 cell line, CRL-11268 obtained from ATCC) were cultured in growth medium consisting of equal parts of Dulbecco’s modified Eagle’s medium (DMEM) (DMEM Glutamax, Gibco, Thermo Fisher Scientific) and Ham’s F10 nutrient mix (Gibco, Thermo Fisher Scientific), supplemented with 10% fetal calf serum (FCS) and 1% penicillin–streptomycin. Cells were split twice a week by trypsinization and grown at 37°C with a humidified atmosphere containing 5% CO2. STR profiling confirmed a 100% match with the HEK 293T cellline (ATCC). Human SCN2A (D-splice variant), encoding for the alpha subunit of the NaV1.2 channel was cloned in pcDNA3.1-IRES-GFP, and SCN1B/SCN2B, encoding for beta subunits 1 and 2, was cloned into pcDNA3.1. These vectors were described previously (Ben-Shalom et al., 2017) and obtained from Genscript (Genscript, USA). The constructs were amplified in Stbl3 bacteria (Genscript, USA) and were purified using the GeneJET Plasmid Maxiprep kit (ThermoFisher, USA) according to the manufacturer’s protocols. The plasmids were transiently transfected into 70% confluent HEK-293 cells plated in 12-well plates. Per well, the transfection cocktail contained 500 ng pcDNA3.1-SCN2A-IRES-GFP, 290 ng pcDNA3.1-SCN1B- IRES-SCN2B and 5 μL of polyethylenimine (PEI) diluted in 100 μL 1% saline, incubated for 20 min at room temperature before addition to the culture medium. Cells were incubated with 100 μL of transfection cocktail in 1 mL of culture medium for 24 hr at 37°C in a humidified atmosphere containing 5% CO2. Cells were trypsinised and used for electrophysiological recording typically 48 hr after transfection.
Electrophysiological recording from neurons
Request a detailed protocolFor patch-clamp recording, slices were transferred to a customized upright microscope (BX51WI, Olympus Nederland BV, or LNscope, Luigs and Neumann, Ratingen, Germany). The transmitted light path consisted of a custom made 850 nm Light Emitting Diode (LED) light source (LZ1-10R602, LED Engin, CA), collimated using an aspheric condenser lens (ACL50832U-B, Thorlabs, Germany) and adapted to the microscope lamp port using a custom 3D printed adapter and passed through an oblique or Dodt illumination condenser (WI-OBCD, Olympus). The top 50 µm of the slice surface was visualized using an optical pathway consisting of a 60× water immersion objective (N.A. 1.0, LUMPLFLN60XW, Olympus or N.A. 1.1, LUMPLFLN60XW, Olympus), 2× intermediate zoom attachment (U-ECA, Olympus), camera splitter (U-TRU, Olympus) with inbuilt 180 mm tube lens on the back port and a 0.63× demagnifier (U-TV0.63XC, Olympus) projected the final image onto a high resolution CCD camera (CoolSNAP-EZ, Photometrics), which was operated using μManager (Edelstein et al., 2014). Based on the bright-field image large L5 neurons with an intact axon parallel and close to the surface were targeted for recording. Current-clamp recordings were made with Dagan BVC-700A amplifiers (Dagan Corporation, MN, USA) or AxoClamp 900A (Molecular Devices Limited, UK). An Axopatch 200B (Molecular Devices) was used for voltage-clamp and AP-clamp experiments. The microscope bath was perfused with oxygenated (95% O2, 5% CO2) ACSF consisting of (in mM): 125 NaCl, 3 KCl, 25 glucose, 25 NaHCO3, 1.25 Na2H2PO4, 2 CaCl2, and 1 MgCl2. Patch pipettes were pulled from borosilicate glass (Harvard Apparatus, Edenbridge, Kent, UK) pulled to an open tip of 3–6 MΩ resistance. For all current-clamp, subthreshold voltage-clamp ramp and AP-clamp recordings the intracellular solution contained (in mM): 130 K-Gluconate, 10 KCl, 4 Mg-ATP, 0.3 Na2-GTP, 10 HEPES and 10 Na2-phosphocreatine (pH 7.25 adjusted with KOH, 280 mOsmol kg−1). The liquid junction potential difference of –13.5 mV was corrected in all recordings. For morphological reconstruction, 5 mg ml−1 biocytin was routinely added. Voltage recordings were analogue low-pass filtered at 10 kHz (Bessel) and digitally sampled at 100 kHz using A-D converter (ITC-18, HEKA Elektronik Dr. Schulze GmbH, Germany) and data acquisition software Axograph X (v.1.5.4, Axograph Scientific, NSW, Australia). Bridge-balance and capacitances were fully compensated based on small current injections leading to minimal voltage errors. The recording temperature was 33 ± 1°C. Only cells with a stable bridge-balance (<25 MΩ), resting membrane potential and AP shape throughout the recording session were included in the analysis.
For voltage-clamp recordings of INa and ICa (Figure 5) the bath was perfused with oxygenated (95% O2, 5% CO2) extracellular recording solution consisting of (in mM): 100 NaCl, 3 KCl, 25 glucose, 25 NaHCO3, 1.25 Na2H2PO4, 2 CaCl2, 1 MgCl2, 5 4-AP, 20 TEA-Cl, 0.02 CNQX, 0.05 D-AP5, 0.02 ZD-7288, 0.01 XE991 and 0.003 Gabazine (SR-95531). The intracellular solution contained (in mM): 130 CsCl, 10 TEA-Cl, 10 HEPES, 4 Mg-ATP, 5 Na2-phosphocreatine and 0.3 Na2-GTP (pH 7.25 adjusted with CsOH, 280 mOsmol kg−1). A liquid junction potential of –5.6 mV was applied to the recordings. Series resistance was routinely compensated to >75% and the linear leak and capacitance off-line subtracted using a P/9 protocol with 10-fold scaled pulses. Current recordings were analogue low-pass filtered at 10 kHz (Bessel) and digitally sampled at 100 kHz using A-D converter (ITC-18, HEKA Elektronik Dr. Schulze GmbH, Germany) and data acquisition software Axograph X (v.1.5.4, Axograph Scientific, NSW, Australia). To improve voltage-clamp of the large and rapid Na+ currents the recordings were made at room temperature (~20°C).
Electrophysiological recordings from HEK-293 cells
Request a detailed protocolFor recordings from HEK-293 cells they were transferred to a recording chamber which was continuously perfused with extracellular solution, containing (in mM): 135 NaCl, 4.5 KCl, 2 CaCl2, 1 MgCl, 10 HEPES and 11 Glucose. The intracellular solution contained (in mM): 110 CsF, 10 NaCl, 20 EGTA and 10 HEPES. In the OGB-1 experiments, we added 100 µM OGB-1, EGTA was omitted and CsF raised to 120 mM instead of 110 mM. The liquid junction potential difference of –10 mV was corrected for. Whole-cell patch-clamp recordings were made ~48 hr after transfection. Cells were recorded at room temperature (~20°C) and continuously perfused with extracellular solution at a flow rate of 1.5 mL∙min–1. Patch pipettes were pulled to a resistance of 2–3 MΩ. Round, isolated cells with a diameter >10 μm, a smooth cell surface and a moderate EGFP fluorescent signal were selected for recordings (Figure 7a). HEK-293 cells had an average capacitance of 9.19 ± 0.70 pF (n = 16). The holding potential was –70 mV and voltage dependence of activation of NaV1.2 was determined by an activating protocol consisting of a hyperpolarizing pulse to –130 mV (20 ms) followed by step pulses from –80 mV to +50 mV with increments of 10 mV for 20 ms. Voltage dependence of inactivation was assessed with voltage pulses from –130 mV to –30 mV with increments of 10 mV for 100 ms duration, followed by a depolarizing pulse to –20 mV for 20 ms. A P/5 leak subtraction protocol (10-fold scaling) was used to subtract remaining capacitive and leak currents. Series resistance was not compensated.
Blockers and toxins
Request a detailed protocolEGTA and blockers were added to the appropriate concentration to the ACSF and perfused. The extracellular Ca2+ ([Ca2+]o) was lowered by bath application of 2.5 mM EGTA and using the online maxchelator tool (https://somapp.ucdmc.ucdavis.edu/pharmacology/bers/maxchelator/CaMgATPEGTA-TS.htm; Bers et al., 2010) we calculated the remaining [Ca2+]o to be 437 nM, based on a recording temperature of 35°C, a pH of 7.4 and an ionic strength of 0.15 M of the free ions in our solution. To limit hyperexcitability in the presence of EGTA, we added synaptic blockers to the ACSF (20 µM CNXQ and 50 µM D-AP5) and kept the a-EPSP voltage peak amplitude constant by reducing the amplitude of the current injections in both control and EGTA measurements (Figure 2). To prevent precipitation of Ni2+ we used phosphate-free extracellular solutions containing (in mM): 126.25 NaCl, 3 KCl, 25 glucose, 25 NaHCO3, 2 CaCl2, 1 MgCl2 and 0.1 Ni2+. A > 10% increase in fluorescence baseline was observed in some experiments (3 out of 8 recordings), which were subsequently excluded. Two blockers (SNX-478 and ω-conotoxin MVIIC) were not perfused but were locally puffed using a Picospritzer III (Intracel) for 3 s ending 0.5 s before imaging to avoid vibration artifacts. Bovine serum albumin (0.1 mg/ml) was added to the rACSF before ω-conotoxin MVIIC was introduced to minimize non-specific binding of the drug.
Ca2+ and Na+ imaging
Request a detailed protocolTo optically record [Ca2+]i and [Na+]i in axons, membrane impermeable Ca2+ and/or Na+ indicators were added to intracellular solutions. For Ca2+ we used OGB-1 (100 µM), OGB-5N (1 mM) or bis-Fura-2 (200 µM) and for Na+ imaging we used sodium-binding benzofuran isophthalate (SBFI, 1–1.5 mM). Patch pipettes were first filled with dye-free solution for half of the tapered part of the pipette tip, then backfilled with the dye-containing solution. Fluorescence intensity at the AIS was monitored during dye loading and imaging started only when the indicators were fully equilibrated (typically after 0.5–1 hr). Optical recordings of Ca2+ or Na+ dye fluorescence changes were obtained with wide-field epifluorescence microscopy. Fluorescence was collected by the same 60× water immersion objective, passed through the microscope tube lens (U-TR30IR, Olympus) and projected onto a rapid data-acquisition camera with relatively low spatial resolution (80 × 80 pixels) but high dynamic range (14 bits) and low read noise (NeuroCCD-SM, RedShirtImaging LLC, Decatur, GA) via a 0.1× or 0.35× demagnifier. The CCD frame corresponded to an area of approximately 320 or 91 µm2 in the object plane with each individual pixel receiving light from an area of ∼4 × 4 or 1.1 × 1.1 µm2, respectively. High-speed recordings (20 kHz) were always performed with the 0.35× demagnifier and with 3 × 3 binning of pixels on the chip, the ultimate pixel sizes in these recordings were ∼3.4 × 3.4 µm2. The two recordings from human cells were performed under a 100× NA 1.1 Nikon objective (MRL07920) in combination with a 0.1× demagnifier, resulting in a pixel size of ∼2.4 × 2.4 µm2 (see Figure 1—figure supplement 1).
The epifluorescence light path consisted of an excitation LED light source, which was collimated using an aspheric lens (ACL5040U-A, Thorlabs, Germany) and the appropriate excitation filter, dichroic mirror and emission filter. For OGB-1 and OGB-5N excitation a 470 nm LED was used (SP-01-B4, Luxeon Star LEDs, Canada), the excitation light was filtered with 475/30 nm (475 nm center wavelength, 30 nm wide) excitation filter, reflected to the preparation by a dichroic mirror with a central wavelength of 500 nm and the fluorescence light was passed through a 520 nm barrier filter (U-MWB2 cube, Olympus). For SBFI excitation, LED light (365 nm LED, LZ1-10UV00, Ledengin, USA) was filtered by 357/44 nm filter (FF01-357/44-25, Semrock), a dichroic mirror with edge at 415 nm (Di03-R405-t1−25 × 36, Semrock) reflected excitation light to the sample and the emission light was then passed through a long-pass colored glass filter with the edge at 400 nm (FGL400, Thorlabs). For combined Na+ and Ca2+ imaging (Figure 6—figure supplement 1), the light from the 365 and 470 LEDs was combined by a dichroic mirror with edge at 458 nm (FF458-Di02, Semrock) and the filter set switched between trials. Light was directed through a fluorescence illuminator (BX-RFA, Olympus) equipped with a rectangular field stop providing an open area of 150 × 250 µm to reduce phototoxicity (U-RFSS, Olympus). The cell body was positioned just outside the field stop and the axon in the middle parallel to the long side (see e.g. Figure 1).
The critical benefit of epifluorescence measurements over two-photon imaging is increased light collection (~90% quantum efficiency, low-read noise of the CCD camera) enabling a high sensitivity and temporal fidelity. We optimized all imaging parameters to obtain maximal signal to noise ratio, which allowed us to image at the maximum acquisition rate of 20 kHz. In addition to the light collection optimization and selective targeting of superficial neurons, multiple trials were averaged to improve signal-to-noise ratio (typically 20–40). Fluorescence signals were temporally aligned to the electrophysiological voltage or current signals. For optical recording of ICa, which requires the transformation of ΔF/F into the first time derivative, Ca2+ binding to the indicator must be proportional to [Ca2+]i and the endogenous buffering capacity to be low. Based on the submillisecond equilibration time of OGB-5N and imaging at the maximally possible frame rate of 20 kHz it was recently shown that in CA1 hippocampal neurons these conditions are met and optically recorded ICa tracks electrically recorded ICa, enabling the identification of CaV channel subtypes in dendrites (Jaafari et al., 2014). Considering the low buffering capacities of endogenous buffers in the axon (κs ≈ 20) (Jackson and Redman, 2003; Delvendahl et al., 2015), we employed this technique in the AIS. Ca2+ imaging in HEK-293 cells was performed with 0.1 mM OGB1 added to the HEK-293 cell intracellular solution (from which EGTA was omitted). The fluorescence was recorded during a 200 Hz 1 s train of depolarizing pulses from –120 to –30 or –20 mV (corrected for liquid junction potential).
Voltage imaging
Request a detailed protocolVoltage imaging in neurons was performed as reported previously (Hamada et al., 2017). Neurons were filled with intracellular solution containing JPW3028 (0.8 mM) for typically 1 hr at room temperature, after which the patch pipette was retracted and the dye was left to diffuse into the lipid membranes for 1–4 hr. Subsequently the bath temperature was increased to 35°C and the cell was re-patched with normal intracellular solution. A 530 nm LED (SP-05-G4, Luxeon Star LEDs, Canada) was used for excitation of the dye. The excitation light was filtered with a 530/20 nm filter (BP510-550, Olympus), reflected to the sample by a dicroic mirror with a center wavelength of 570 nm (DM570, Olympus) and the emission light filtered by a 590 longpass filter (BA590, Olumpys). Data were collected at 20 kHz and low-pass filtered by a binomial filter (one pass) and averaged over 20–30 trials. Voltage imaging in HEK-293 cells was performed identically, with the only exception being that the experiments were performed at 20°C and the dye diffused equally in the small round cells, so imaging experiments were initiated 20 min after obtaining whole-cell configuration. The cells were maintained at –75 mV holding potential and the fluorescence of JPW3028 recorded at 1 kHz. The voltage command consisted of 100 ms steps of 50 mV increasing steps with a maximal step of +250 mV relative to holding potential. The average ΔF/F per voltage step was defined as the first 20 frames of each bleach corrected and normalized voltage step.
Imaging data analysis
Request a detailed protocolImaging data analysis was performed using Neuroplex (Redshirt imaging), Axograph and Excel. Fluorescence signals were always background-subtracted. To correct for bleach effects, every 5th trial was recorded without current injection. A first order exponential was fitted to the average of the bleach trials and normalized to the peak. The average of the signal trials was divided by this trace to correct for bleach decay. Values for each ROI were defined as a fractional fluorescence change (ΔF/Fbaseline), where Fbaseline is the raw intensity average of 10 frames before the signal (subthreshold or AP) was initiated. Pixels were color coded with ‘physics’ color scheme from FIJI image processing software (NIH, USA) (Schindelin et al., 2012). For both OGB-1 and OGB-5N, we recorded the fluorescence in response to subthreshold stimuli, single APs and multiple APs. The ΔF/F response to subthreshold stimuli and a single AP was always below dye saturation. For OGB-1, we recorded trains of APs and observed that a single AP was 25 ± 2.7% of dye saturation (n = 8) and for OGN-5N, we observed a linear increase from 1 to 3 APs (n = 3), indicating that the fluorescence of a single AP was far from dye saturation.
Optical current measurements
Request a detailed protocolTo compare kinetics between electrically and optically recorded currents, electrical currents were first downsampled to 20 kHz (optical acquisition rate). The optical trace was differentiated and then inverted (to mimic INa, which is conventionally depicted as a negative, inward current). All current traces were then filtered with a 3-window binomial filter of 50–150 iterations (generally 100). The traces were baselined to the current before the onset of the fast current and normalized to the peak of the current. Because the ICa_opt traces were too noisy to be fitted with a single exponential fit, conventionally used to obtain activation rise time, we used a Boltzmann sigmoid function to obtain the slopes of the traces.
We fitted all traces using Axograph and the slope values (k) were used to compare activation kinetics. Although the INa peak amplitude recorded at the soma highly varied between neurons (–36.61 ± 6.53 nA), the slope was nearly constant (207 ± 0.007 µs, n = 17).
Ratiometric imaging
Request a detailed protocolTo estimate the absolute [Na+]i and [Ca2+]i in response to a single AP, we used a ratiometric imaging approach. Patch pipettes were front filled with clear intracellular solution and back-filled with intracellular solution containing 1.5 mM SBFI (Invitrogen) or 200 µM bis-Fura-2 hexapotassium salt (bF2, Biotium). Fluorescent emission of ratiometric indicators depends on the ionic concentration and the excitation wavelength where an increase in [Ca2+]i produces an increase in bF2 fluorescence with the wavelength of 340 nm but a decrease with 385 nm. On the other hand, with SBFI an increase in [Na+]i decreases SBFI fluorescence at 340 nm but does not alter at 385 nm excitation wavelength. Using the ratio (R) corrects for differences in cytosolic volume or dye concentration differences along imaged compartments (Langer and Rose, 2009). The sources of excitation light were two LEDs (Thorlabs) with peaks at 340 nm and 385 nm, fitted with band pass excitation filters at 340/22 and 387/11 nm (FF01-340/22-25 and FF01-387/11-25, Semrock) and combined by a dichroic mirror with a central wavelength of 376 nm (FF376-Di01−25 × 36, Semrock). The excitation light was reflected to the sample by a dichroic mirror with a central wavelength at 405 nm (Di01-R405−25 × 36, Semrock) and passed through the objective to the sample. The fluorescent emission signals were passed through a 420 long pass filter (Thorlabs). Ratiometric imaging was performed by alternatingly triggering each LED at the frame rate of the camera, as described previously (Miyazaki and Ross, 2015). This was achieved by combining custom designed Arduino/Parallax machines with Cyclops LED drivers. These hardware solutions allowed us to digitally control the voltage driving the LED, thus having maximum control over excitation light intensity. Fluorescence emission signals originating from each LED were separated with custom written software (FrameSplitter.txt, Battefeld et al., 2018). The camera operated at 0.5–1.0 kHz, resulting in a ratiometric frame rate of 0.25–0.5 kHz. Per experiment the fluorescent signals were averaged for 40 to 120 trials. The ratio R was defined as F1/F2, where F1 and F2 are the background-subtracted fluorescence intensities at excitation with 340 nm and 385 nm, respectively.
Calibration of ratiometric imaging
Request a detailed protocolIn order to scale ratiometric bF2 signals to absolute changes in Ca2+ concentration, we used the standard equation for ratiometric measurements (equation 1 in Figure 8—figure supplement 2; Grynkiewicz et al., 1985), which depends on KD, the dissociation constant, Rmin and Rmax, the ratio in zero and dye-saturating Ca2+, respectively and the scaling factor (Sf2/Sb2), defined as the fluorescence intensity at excitation with 385 nm of zero Ca2+ divided by saturating Ca2+. These values were determined in an ex situ calibration, by measuring the ratiometric signal of solutions containing 0 Ca2+ and a high [Ca2+] (Figure 8—figure supplement 2). The solutions closely mimicked intracellular solutions and contained (in mM): 110 K-gluconate, 4.4 or 21 KCl, 0 or 10 CaCl2, 3.8 or 5.36 MgCl2, 10 HEPES, 4 Mg-ATP, 0.3 Na2-GTP, 10 Na2-phosphocreatine, 10 EGTA, and 0.2 bF2. The final free [Ca2+] depends on interaction between Ca2+, Mg2+ and EGTA and was calculated using the maxchelator tool (Bers et al., 2010).
[Ca2+]min was 0 and [Ca2+]max was 4.39 µM. We repeated the calibration experiment three times. In our experimental setting, the KD of bF2 was 507.3 ± 5.7 nM, which matches with the reported value of 525 nM in the presence of Mg2+ (Molecular Probes Handbook, Thermofischer). We then performed ratiometric imaging in the AIS in response to a single AP. ∆R/R was calculated by dividing every ratio by the average of the baseline ratio before the onset of the AP (similar to the conventional ∆F/F). We used the KD as determined from our calibration experiments. The Rmin was scaled to be ~95% of Rbaseline to result in resting [Ca2+]i of 50 nM and Rmax as established in the calibration experiments. If the LED intensity of the cellular recording was different from the calibration intensity used during the calibration experiments, the Rmax was corrected linearly, assuming that Rmin/Rmax was constant. These experiments showed that after a single AP, [Ca2+]free in the AIS rises with 55.6 ± 12.6 nM (see Figure 8—figure supplement 2).
Because Rmin and Rmax were not measured in situ we verified the Δ[Ca2+]free with an alternative analysis that is independent of the exact values for KD, Rmin and Rmax (equation 2 in Figure 8—figure supplement 2; Langer and Rose, 2009). In this approach, changes in fluorescence ratio ΔR/R0 (%) are plotted versus [Ca2+]free, showing a nearly linear increase in ΔR/R0 (%) for low [Ca2+]free (see Figure 8—figure supplement 2). A linear fit to the region of [Ca2+]i between 0 and 193 nM indicated that a 1% increase in ΔR/RbF2 corresponded to a Δ[Ca2+]i of ~10.4 nM (R2 = 0.99, six concentrations, n = 3 repetitions). We measured an AP-evoked Δ[Ca2+]free of 52.5 ± 12.2 nM (Figure 8—figure supplement 2), in good support of the standard Ca2+ measurement approach. We analyzed the ratiometric SBFI data only using this second approach, which is standard for SBFI measurements (Langer and Rose, 2009). The two base calibration solutions contained (in mM): 130 K-Gluconate or Na-Gluconate, 10 KCl of NaCl, 0.3 Tris2-GTP or Na2-GTP, 10 HEPES, 4 Mg2+ATP, 10 Tris2-phospocreatine or Na2-phospocreatine, and 1.5 SBFI, pH 7.25 adjusted with Tris base. These two base solutions provided a range of 0–160.6 mM [Na+]. When normalized to the ratio obtained in Na+-free solution (R0), a 1% increase in ΔR/RSBFI corresponded to a Δ[Na+]i of 0.35 mM for Δ[Na+]i between 0 and 48 mM (R2 = 0.98; eight concentrations, n = 3 repetitions, Figure 8—figure supplement 2). An AP evoked a Δ[Na+]i of 1.49 ± 0.2 mM (Figure 8—figure supplement 2).
Immunofluorescence staining
Request a detailed protocolFollowing imaging experiments, the slices were fixed using 4% PFA in 0.1 M phosphate-buffered saline (PBS), pH 7.4 for 20 min and stored in 0.1 M PBS, pH 7.4 at 4°C. For triple immunohistological labeling the slices were washed three times in PBS and then incubated in a blocking solution (10% normal goat serum, 0.5% Triton X-100 in PBS) at room temperature for two hours, followed by 24 hr incubation at room temperature in the blocking solution containing primary antibodies: streptavidin Alexa-488 conjugate (1:500; Invitrogen), primary antibody for giant saccular organelle Synaptopodin (rabbit; 1:500; Sigma-Aldrich Chemie) and antibody for AIS marker: Ankyrin G (mouse; 1:100; Neuromab) or ßIV-spectrin (mouse; 1:250; Neuromab, see also Key Resources Table). The slices were 3x washed in 0.1 M PBS and then incubated with secondary antibodies: Alexa-555 goat anti rabbit IgG (1:500; Invitrogen) and Alexa 633 goat anti mouse IgG (1:500; Invitrogen). Subsequently, the slices were 3x washed in 0.1 M PBS and mounted with Vectashield mounting medium with 4,6-diamidino-2-phenylindole (DAPI; Vector Laboratories). Images (bit depth, 8) were collected as described previously (Hamada et al., 2016). To align confocal images and the Ca2+ fluorescence images of the RedShirt CCD camera (Figures 1, 2 and 4) we used the original calibrated images of the two systems. The maximum Ca2+ fluorescence image was calibrated within the original optical path. We overlaid the maximum Ca2+ fluorescence image of the RedShirt camera and the maximum projection of the streptavidin image of the neuron morphology from confocal microscopy within ImageJ and applied only a rotation translation to visually match the two images based on the AP-evoked Ca2+ signals spreading into dendrites and axons.
Model simulations with single compartment
Request a detailed protocolAll model simulations were performed with NEURON (v.7.5) (Hines and Carnevale, 2001). A single compartment was created with length and diameter dimensions of 10 µm and nseg = 10, with specific membrane capacitance of 1.0 µF cm–2, specific membrane resistance of 25 kΩ cm2 and specific axial resistivity of 150 Ω cm. The resting membrane potential set to –77 mV using e_pass. Conductance models for Ca2+ were based on the high-voltage activated (CaH) and a T-type CaV channel model (CaT) obtained from ModelDB (https://senselab.med.yale.edu/ModelDB/) (Mainen and Sejnowski, 1996). Ca2+ conductivity of NaV channels was modeled by including a standard ohmic Ca2+ ion mechanism with a reversal potential (eca) of +140 mV into a mathematical 8-state Na+ conductance model, computing simultaneously voltage- and time-dependence of the Ca2+ current ICa(Na) and INa based on experimentally constrained rate constants of somatodendritic and axonal INa (Schmidt-Hieber and Bischofberger, 2010). The kinetics of the voltage-gated conductance models was examined by fitting the current rise times with an exponential function for a –35 mV command potential, resampling the simulated traces to 20 kHz. The results showed that the INa in the model activated with 240 µs and ICa(Na) (0.5% conductivity ratio) with 280 µs. In comparison, ICaT activated with 4.88 ms and ICaH with 6.51 ms. These time constants are well in range of the experimentally determined values for the TTX- and Ni2+-sensitive components recorded at the soma (Figure 5). For Figure 8—figure supplement 1 we used an AP recorded from the L5 pyramidal neuron AIS at 100 kHz (threshold-to-peak, 94 mV; half-width duration of 285 µs [Hallermann et al., 2012] as the command potential in VectorPlay linked to the SEClamp function in NEURON (with Rs being infinitely small). Single compartmental models were run at dt of 10 µs at a nominal temperature of 33°C.
Model simulations with a multicompartmental model
Request a detailed protocolConductance-based multi-compartmental simulations were performed with an anatomically realistic reconstructed rat L5 pyramidal neuron (NeuroMorpho.Org ID: NMO_75667, Neuron Name 2014-04-01_1). The morphology was acquired with a confocal microscope at 2048 × 2048 pixels (1.0 µm z-steps, Leica SP8) using a 40× oil immersion objective (NA 1.3) scanning both the biocytin-streptavidin fluorescence and the ßIV-spectrin fluorescence. Uncompressed image stacks (~20 GB) were imported and reconstructed into Neurolucida (v.10, MBF Bioscience Inc, Germany), compartmentalized for the AIS and nodes as described previously (Hamada et al., 2016) and imported with the 3D import tool in NEURON (Carnevale and Hines, 2006). Multicompartmental simulations were performed to estimate the detailed ionic accumulation, concentration and diffusion in the proximal sites of the axon and match our experimental recordings as close as possible. Ca2+ diffusion, buffering and pump (cdp) mechanisms were implemented based on the algorithms described in the NEURON book (Chapter 9, example 9.8 in Carnevale and Hines, 2006) and on a previously published Ca2+ model (Fink et al., 2000) (available at ModelDB, accession number 125745, https://senselab.med.yale.edu/ModelDB/). We implemented cdp.mod (Fink et al., 2000) with the following alterations: we removed all SERCA related parameters, updated some starting values to our experimental conditions and extended the models to report not only [Ca2+]i, but also to simulate the Ca2+ indicator response ΔF/F, using the equation:
With c being a constant to scale simulated ΔF/F. Because the equation was used to match the simulation to experimental data with regard to the temporal dynamics of Ca2+ extrusion, the absolute amplitude of ΔF/F was not used and c was set to a nominal value of 6. The different Ca2+ indicators used experimentally were implemented by adjusting the concentration of the exogenous buffer, and its known or measured KD. Static Ca2+ buffering properties of endogenous organelles (κs) were simulated with a TBufs of 100–400 µM and KDs of 10 µM, to mimic a κs of 10–40 (Jackson and Redman, 2003; Delvendahl et al., 2015).
To constrain the peak Na+ conductance densities (Na) we injected a 3 ms square current pulse in the somatic compartment and iteratively adjusted Na and K. We varied both their peak conductance densities as well as the voltage-dependence of activation of NaV and KV channels by constraining the model AP to the experimentally recorded AP of the same neuron, with the aim to match the AP both in the V-t as well as the phase-plane dimensions recorded and simulated at 100 kHz (Figure 7c). To further constrain Na we compared the AP-evoked [Na+] with the experimental recordings using the ratiometric indicator SBFI, yielding a Δ[Na+]i of on average ~1.5 ± 0.2 mM (n = 5, imaged at 0.5 kHz; Figure 8—figure supplement 2). CaV channels were incorporated based on previously published models (Mainen and Sejnowski, 1996) and CaV channel conductance was separated in high- and low-voltage (T-type) activated channels and was varied between 2 and 4 pS µm–2 in the AIS, 8 and 4 pS µm–2 in the soma and ranged between 0.5 and 4 pS µm–2 in the dendrites.
Statistics and data availability
Request a detailed protocolAll statistical tests were done in GraphPad Prism 8 (version 8.1.2, GraphPad Software, Inc). Sample sizes for the pharmacological experiments were estimated based on the following assumptions: to observe a 50% block (based on Bender and Trussell, 2009) with 25% standard deviation (relative to mean) with a power of 0.8 and a type I error probability of 0.05, we would need a minimum of 4 paired recordings per treatment (PS Software version 3.1.6).
The cutoff significance level (P) was 0.05. Control peak ΔF/F values at the AIS in response to both subthreshold and AP signals were tested for normality. Since both data sets passed the D’Agostino and Pearson normality test, parametric tests were used to test all differences between peak OGB-1 ΔF/F. To compare the spatial differences in signal amplitude we used one-way ANOVAs with multiple comparisons with Tukey correction for false positives. A linear regression was used to assess the synaptopodin and AIS marker length (Ankyrin G or ß4-spectrin) relationship. We used one-tailed ratio (compared log differences in the data set) paired t-test when analyzing all our toxin data. Differences between toxin and control give a measure of absolute reduction; differences between logarithms give a measure of relative reduction log toxin – log control = log (toxin/control). One tailed test was used on the premise that toxins reduce Ca2+ signals. The exception was the effect of NCX for which we used a two-tailed ratio paired t-test. The OGB-5N peaks in response to subthreshold depolarizations passed the D’Agostino and Pearson normality test and to compare the effects of TTX and QX-314 on subthreshold Ca2+ responses we used two-way-ANOVA with Sidak’s correction for false positives. The following data sets passed the Shapiro-Wilk test for normality, so we compared the means using parametric tests: rise and decay times of OGB1 before and after Ca2+ store release block, peak ΔF/F JPW3028, peak ΔF/F OGB-5N and slopes of INa, ICa_opt and ICa, the ratio of subthreshold peak to AP peak between sodium (SBFI) and calcium (OGB-5N) fluorescence. The following data did not pass the Shapiro-Wilk test for normality, so we compared the means using nonparametric tests: the slope AP peak between sodium (SBFI) and calcium (OGB-5N) fluorescence, the peak INa and peak OGB-1 ΔF/F in HEK-293 cells. All data generated or analyzed are in the manuscript or supporting files. The source data files are provided for Figures 1–7 and Table 1. The NEURON model morphology is available at NeuroMorpho.Org ID: NMO_75667, Neuron Name 2014-04-01_1 and the mod file used to model Ca2+ diffusion and buffering is available at ModelDB, accession 125745 (https://senselab.med.yale.edu/ModelDB/), with adjustments described in ‘Model simulations with a multicompartmental model’.
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 1 to 7 and Table 1.
References
-
A generalised method to estimate the kinetics of fast ca(2+) currents from ca(2+) imaging experimentsJournal of Neuroscience Methods 268:66–77.https://doi.org/10.1016/j.jneumeth.2016.05.005
-
Tetrodotoxin-sensitive calcium-conducting channels in the rat hippocampal CA1 regionThe Journal of Physiology 450:529–546.https://doi.org/10.1113/jphysiol.1992.sp019141
-
Depolarization and calcium entry in squid giant axonsThe Journal of Physiology 218:709–755.https://doi.org/10.1113/jphysiol.1971.sp009641
-
Loss of the cisternal organelle in the axon initial segment of cortical neurons in synaptopodin-deficient miceThe Journal of Comparative Neurology 504:441–449.https://doi.org/10.1002/cne.21445
-
A versatile and Open-Source rapid LED switching system for One-Photon imaging and Photo-ActivationFrontiers in Cellular Neuroscience 12:530.https://doi.org/10.3389/fncel.2018.00530
-
Ca2+-activated K+ channels: from protein complexes to functionPhysiological Reviews 90:1437–1459.https://doi.org/10.1152/physrev.00049.2009
-
Calcium microdomains: organization and functionCell Calcium 40:405–412.https://doi.org/10.1016/j.ceca.2006.09.002
-
A practical guide to the preparation of ca(2+) buffersMethods in Cell Biology 99:1–26.https://doi.org/10.1016/B978-0-12-374841-6.00001-3
-
The impact of BK channels on cellular excitability depends on their subcellular locationFrontiers in Cellular Neuroscience 10:206.https://doi.org/10.3389/fncel.2016.00206
-
Axonal calcium entry during fast 'sodium' action potentials in rat cerebellar Purkinje neuronesThe Journal of Physiology 495:641–647.https://doi.org/10.1113/jphysiol.1996.sp021622
-
Reversal potential of the calcium current in bull-frog atrial myocytesThe Journal of Physiology 403:267–286.https://doi.org/10.1113/jphysiol.1988.sp017249
-
D3 receptors regulate excitability in a unique class of prefrontal pyramidal cellsThe Journal of Neuroscience 37:5846–5860.https://doi.org/10.1523/JNEUROSCI.0310-17.2017
-
Layer- and cell-type-specific suprathreshold stimulus representation in rat primary somatosensory cortexThe Journal of Physiology 581:139–154.https://doi.org/10.1113/jphysiol.2006.124321
-
Advanced methods of microscope control using μmanager softwareJournal of Biological Methods 1:e10.https://doi.org/10.14440/jbm.2014.36
-
An image-based model of calcium waves in differentiated neuroblastoma cellsBiophysical Journal 79:163–183.https://doi.org/10.1016/S0006-3495(00)76281-3
-
Tonotopic variation of the T-Type Ca2+ Current in Avian Auditory Coincidence Detector NeuronsThe Journal of Neuroscience 38:335–346.https://doi.org/10.1523/JNEUROSCI.2237-17.2017
-
A new generation of Ca2+ indicators with greatly improved fluorescence propertiesThe Journal of Biological Chemistry 260:3440–3450.
-
Molecular identification of a TTX-sensitive ca(2+) currentAmerican Journal of Physiology. Cell Physiology 280:C1327–C1339.https://doi.org/10.1152/ajpcell.2001.280.5.C1327
-
State and location dependence of action potential metabolic cost in cortical pyramidal neuronsNature Neuroscience 15:1007–1014.https://doi.org/10.1038/nn.3132
-
Loss of saltation and presynaptic action potential failure in demyelinated axonsFrontiers in Cellular Neuroscience 11:45.https://doi.org/10.3389/fncel.2017.00045
-
The permeability of the sodium channel to metal cations in myelinated nerveThe Journal of General Physiology 59:637–658.https://doi.org/10.1085/jgp.59.6.637
-
NEURON: a tool for neuroscientistsThe Neuroscientist 7:123–135.https://doi.org/10.1177/107385840100700207
-
Imaging fast calcium currents beyond the limitations of electrode techniquesBiophysical Journal 107:1280–1288.https://doi.org/10.1016/j.bpj.2014.07.059
-
Calcium dynamics, buffering, and buffer saturation in the boutons of dentate granule-cell axons in the hilusThe Journal of Neuroscience 23:1612–1621.https://doi.org/10.1523/JNEUROSCI.23-05-01612.2003
-
A unique ion channel clustering domain on the axon initial segment of mammalian neuronsJournal of Comparative Neurology 522:2594–2608.https://doi.org/10.1002/cne.23551
-
Synaptically induced sodium signals in hippocampal astrocytes in situThe Journal of Physiology 587:5859–5877.https://doi.org/10.1113/jphysiol.2009.182279
-
Ion transport through pores: a rate-theory analysisBiochimica Et Biophysica Acta (BBA) - Biomembranes 311:423–441.https://doi.org/10.1016/0005-2736(73)90323-4
-
Ca(V)3 T-type calcium channel isoforms differentially distribute to somatic and dendritic compartments in rat central neuronsEuropean Journal of Neuroscience 24:2581–2594.https://doi.org/10.1111/j.1460-9568.2006.05136.x
-
Calcium inward currents in internally perfused giant axonsThe Journal of Physiology 235:225–265.https://doi.org/10.1113/jphysiol.1973.sp010386
-
Simultaneous sodium and calcium imaging from dendrites and axonsEneuro 2:ENEURO.0092-15.2015.https://doi.org/10.1523/ENEURO.0092-15.2015
-
Molecular basis of ion permeability in a voltage-gated sodium channelThe EMBO Journal 35:820–830.https://doi.org/10.15252/embj.201593285
-
Sodium currents and sodium-current fluctuations in rat myelinated nerve fibresThe Journal of Physiology 329:163–184.https://doi.org/10.1113/jphysiol.1982.sp014296
-
Anatomy and physiology of the thick-tufted layer 5 pyramidal neuronFrontiers in Cellular Neuroscience 9:233.https://doi.org/10.3389/fncel.2015.00233
-
Optical measurement of presynaptic calcium currentsBiophysical Journal 74:1549–1563.https://doi.org/10.1016/S0006-3495(98)77867-1
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Fast sodium channel gating supports localized and efficient axonal action potential initiationJournal of Neuroscience 30:10233–10242.https://doi.org/10.1523/JNEUROSCI.6335-09.2010
-
Differential interactions of na+ channel toxins with T-type Ca2+ channelsJournal of General Physiology 132:101–113.https://doi.org/10.1085/jgp.200709883
-
Active cortical dendrites modulate perceptionScience 354:1587–1590.https://doi.org/10.1126/science.aah6066
-
Intracellular QX-314 inhibits calcium currents in hippocampal CA1 pyramidal neuronsJournal of Neurophysiology 76:2120–2124.https://doi.org/10.1152/jn.1996.76.3.2120
-
Molecular identity of axonal sodium channels in human cortical pyramidal cellsFrontiers in Cellular Neuroscience 8:297.https://doi.org/10.3389/fncel.2014.00297
-
Na(+)-Ca2+ exchange currents in cortical neurons: concomitant forward and reverse operation and effect of glutamateThe European Journal of Neuroscience 9:1273–1281.https://doi.org/10.1111/j.1460-9568.1997.tb01482.x
-
Stimulation-induced ca(2+) influx at nodes of ranvier in mouse peripheral motor axonsThe Journal of Physiology 594:39–57.https://doi.org/10.1113/JP271207
Article and author information
Author details
Funding
Nederlandse Organisatie voor Wetenschappelijk Onderzoek (865.17.003)
- Maarten HP Kole
Stichting voor Fundamenteel Onderzoek der Materie (16NEPH02)
- Maarten HP Kole
National Multiple Sclerosis Society (RG 4924A1/1)
- Maarten HP Kole
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors are indebted to Stefan Hallermann (Univ. Leipzig), Christian Lohmann (NIN–KNAW) and Christiaan de Kock (CNCR, VU) for critically reading earlier versions of the manuscript. We thank Huib Mansvelder (CNCR, VU) for sharing the human tissue. Part of the work was supported by the National Multiple Sclerosis Society (RG 4924A1/1) and the Netherlands Organization for Scientific Research (NWO), including an NWO Vici grant 865.17.003 and a Program Grant 16NEPH02 of the Foundation for Fundamental Research on Matter (FOM) to MK OPA was supported by the Erasmus + Traineeship Program and the Graduate School of Systemic Neurosciences.
Ethics
Human subjects: Written informed consent was obtained from patients and all procedures on human tissue were performed with the approval of the Medical Ethical Committee of the Amsterdam UMC, location VuMC and in accordance with Dutch license procedures and the Declaration of Helsinki. All data were anonymized.
Animal experimentation: All animal experiments were performed in compliance with the European Communities Council Directive 2010/63/EU effective from 1 January 2013. They were evaluated and approved by the national CCD authority (license AVD8010020172426) and by the KNAW animal welfare and ethical guidelines and protocols (DEC NIN 14.49, DEC NIN 12.13, IvD NIN 17.21.01 and 17.21.03).
Copyright
© 2020, Hanemaaijer et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,480
- views
-
- 722
- downloads
-
- 25
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Neuroscience
Williams syndrome (WS; OMIM#194050) is a rare disorder, which is caused by the microdeletion of one copy of 25–27 genes, and WS patients display diverse neuronal deficits. Although remarkable progresses have been achieved, the mechanisms for these distinct deficits are still largely unknown. Here, we have shown that neural progenitor cells (NPCs) in WS forebrain organoids display abnormal proliferation and differentiation capabilities, and synapse formation. Genes with altered expression are related to neuronal development and neurogenesis. Single cell RNA-seq (scRNA-seq) data analysis revealed 13 clusters in healthy control and WS organoids. WS organoids show an aberrant generation of excitatory neurons. Mechanistically, the expression of transthyretin (TTR) are remarkably decreased in WS forebrain organoids. We have found that GTF2IRD1 encoded by one WS associated gene GTF2IRD1 binds to TTR promoter regions and regulates the expression of TTR. In addition, exogenous TTR can activate ERK signaling and rescue neurogenic deficits of WS forebrain organoids. Gtf2ird1-deficient mice display similar neurodevelopmental deficits as observed in WS organoids. Collectively, our study reveals critical function of GTF2IRD1 in regulating neurodevelopment of WS forebrain organoids and mice through regulating TTR-ERK pathway.
-
- Neuroscience
Specialized chemosensory signals elicit innate social behaviors in individuals of several vertebrate species, a process that is mediated via the accessory olfactory system (AOS). The AOS comprising the peripheral sensory vomeronasal organ has evolved elaborate molecular and cellular mechanisms to detect chemo signals. To gain insight into the cell types, developmental gene expression patterns, and functional differences amongst neurons, we performed single-cell transcriptomics of the mouse vomeronasal sensory epithelium. Our analysis reveals diverse cell types with gene expression patterns specific to each, which we made available as a searchable web resource accessed from https://www.scvnoexplorer.com. Pseudo-time developmental analysis indicates that neurons originating from common progenitors diverge in their gene expression during maturation with transient and persistent transcription factor expression at critical branch points. Comparative analysis across two of the major neuronal subtypes that express divergent GPCR families and the G-protein subunits Gnai2 or Gnao1, reveals significantly higher expression of endoplasmic reticulum (ER) associated genes within Gnao1 neurons. In addition, differences in ER content and prevalence of cubic membrane ER ultrastructure revealed by electron microscopy, indicate fundamental differences in ER function.