Modeling analysis of depolarization-assisted afterdischarge in hippocampal mossy fibers

A strong repetitive stimulus can occasionally enhance axonal excitability, leading to the generation of afterdischarge. This afterdischarge outlasts the stimulus period and originates either from the physiological spike initiation site, typically the axon initial segment, or from ectopic sites for spike generation. One of the possible mechanisms underlying the stimulus-induced ectopic afterdischarge is the local depolarization due to accumulated potassium ions surrounding the axonal membranes of the distal portion. In this study, the mechanisms were explored by computational approaches using a simple model of hippocampal mossy fibers implemented with the structure of en passant axons and experimentally obtained properties of ionic conductances. When slight depolarization of distal axons was given in conjunction with the high-frequency stimulus, robust afterdischarges were triggered after cessation of the repetitive stimulus and lasted for a prolonged period after the stimulus. Each spike during the afterdischarge recorded from distal axons precedes that recorded from the soma, suggesting that the afterdischarge was ectopically generated from distal axons and propagated antidromically toward the soma. Notably, when potassium channels in the model are replaced with non-inactivating ones, repetitive stimuli fail to induce afterdischarge. These results suggested that the inactivating property of axonal potassium channels plays a crucial role in generating the afterdischarge. Accumulated inactivation of potassium channels during strong repetitive stimulation may alter mossy fiber excitability, leading to ectopic afterdischarges from sites distinct from the physiological spike initiation region.

Introduction

Neurons are thought to be an analog integrator of excitation and inhibition from multiple inputs converging onto them. Once the membrane potential exceeds the threshold membrane potentials, action potentials are triggered from the spike initiation sites located at the proximal portions of axons (Schmidt-Hieber et al., 2008) or the axon initial segments (Meeks and Mennerick, 2007a; Ohura and Kamiya, 2016; Alpizar et al., 2019). In addition to the canonical mode of spike signaling along axons, it has been suggested that a non-canonical mode called retroaxonal action potentials is triggered from axons distal from the soma in both physiological and pathological conditions (Gutnick and Prince, 1972; Stasheff and Wilson, 1990; Pinault, 1995; Kapoor et al., 1997; Amir et al., 2005; Bucher and Goaillard, 2011; Bähner et al., 2011; Sheffield et al., 2011; Dugladze et al., 2012; Thome et al., 2018; Liu et al., 2022; Kramer et al., 2022). Massive neuronal activity often has an impact to drastically raise neuronal excitability, which causes the non-canonical mode of signaling in the brain, namely ectopic firings from different sites of physiological spike initiation, in an activity-dependent manner (Stasheff and Wilson, 1990; Stasheff et al., 1993; Sheffield et al., 2011). High-frequency activity has been shown to lead to retroaxonal barrage firings of some central neurons such as cortical inhibitory neurons (Sheffield et al., 2011; Suzuki et al., 2014). Induction mechanisms involve astrocytes (Meeks et al., 2005; Deemyad et al., 2018), hyperpolarization-activated cyclic-nucleotide-gated (HCN) channels (Elgueta et al., 2015; Rózsa et al., 2023; see also Roth and Hu, 2020), and K⁺ accumulation by the repetitive stimulus (Smith, 1983; Meeks and Mennerick, 2004, 2007b; Foley et al., 2010; Beswick-Jones et al., 2023; Rózsa et al., 2023). It was found that the hippocampal mossy fiber also displays stimulus-induced prolonged hyperexcitability accompanying the enhancement of spontaneous discharges (Kamiya, 2017), although the mechanisms remain obscure since the hippocampal mossy fiber axons do not significantly express HCN channel subtypes HCN1-4 (Notomi and Shigemoto, 2004). To test for the involvement of membrane depolarization caused by K⁺ accumulation surrounding the axons by the repetitive stimulus in the stimulus-induced afterdischarge, the computational approaches using a NEURON simulation environment were adopted in this study (Hines and Carnevale, 1997). Using a simple and realistic model of hippocampal mossy fiber with the structure of the en passant axons implemented with conductance characterized by the direct electrophysiological recordings from the axonal boutons (Engel and Jonas, 2005; Ohura and Kamiya, 2018), the conditions indispensable to generating stimulus-induced bursts were explored. It was found that a repetitive stimulus with depolarization of distal axons to −70 mV elicited persistent afterdischarge lasting after cessation of the stimulus. The afterdischarge was supposed to be generated ectopically from the distal axons, different from the physiological spike initiation site at the proximal axons or the axon initial segment, suggested by the reversal of the order of the somatic and the axonal spikes after the induction of the afterdischarge. Intriguingly, the stimulus-induced afterdischarge was abolished when the K⁺ channel model was replaced with a non-inactivating one. Taken together, our simulation data suggested that inactivating properties of K⁺ channels are involved in the depolarization-assisted generation of afterdischarge originating from the ectopic site from the physiological generating site of action potentials at the hippocampal mossy fiber synapse.

Materials and methods

Simulation

Computer simulations of axonal membrane potentials (V_m) were performed using NEURON 8.2 for Windows (Hines and Carnevale, 1997). In brief, the structure of the mossy fiber (Acsády et al., 1998; Henze et al., 2000) was approximated by a soma (diameter, 10 μm), 11 axonal cylinders (diameter, 0.2 μm; length, 100 μm), and 10 en passant boutons (diameter, 4 μm). The number of segments was 1 μm⁻¹ in all simulations. For the majority of simulations in this study, the time step was set at 0.1 ms to describe the fast action potential kinetics and the underlying currents. In some simulations for calculating the time derivative of membrane potential (dV_m/dt), the time step was set at 0.01 ms to accurately trace dV_m/dt. The passive electrical properties of the axon were assumed to be uniform, with a specific membrane capacitance C_m of 1 μF cm⁻², a specific membrane resistance R_m of 10,000 Ω cm², and an intracellular resistivity R_i of 110 Ω cm (Engel and Jonas, 2005; Hallermann et al., 2003). The resting membrane potential was set at −80 mV unless otherwise stated.

Voltage-gated Na⁺ channels, K⁺ channels, and leakage channels were inserted into the soma, axon, and boutons, as described by the previous study (Kamiya, 2019). The Na⁺ conductance density was set at 50 mS cm⁻² for the axon and boutons and 10 mS cm⁻² for the soma. The K⁺ conductance density was set at 36 mS cm⁻² throughout all parts of the neurons. Action potentials were evoked by the injection of depolarizing current into the soma (2 ms, 0.2 nA). The equilibrium potentials for Na⁺ and K⁺ ions were assumed to be +50 and −85 mV, respectively. In some simulations, the non-inactivating K⁺ channel model, which lacks inactivation properties, was used instead of the inactivating K⁺ channel model. The hhmfb.mod, a Hodgkin–Huxley type model for the set of sodium, potassium, and leakage channels adapted to channels in mossy fiber terminals, was downloaded from ModelDB (accession no. 128079) for a non-inactivating I_K model for sodium and potassium conductances. For an inactivating I_K model, the potassium conductance in hhmfb.mod was omitted and replaced with KIn.mod, a Hodgkin–Huxley model of inactivating K channels representing the kinetics of K_V1.4 (Wissmann et al., 2003), which was downloaded from the same ModelDB (accession no. 128079).

Action potentials at the mossy fiber boutons, as well as the granule cell somata evoked by the repetitive action potentials of 50 times at 50 Hz, 60 times at 20 Hz, and 50 times at 100 Hz, were calculated for testing the use-dependent initiation of ectopic afterdischarge of the mossy fiber axon. To check for the threshold of stimulation number required for burst generation, stimuli of 40 times at 50 Hz, 50 times at 20 Hz, and 40 times at 100 Hz were also tested, as shown in Figure 1. For local depolarization, the resting membrane potential was changed by altering the el value (equilibrium potential of leak conductance). All model and simulation files will be available from the ModelDB database as accession no. 2018003.¹

Figure 1

Figure 1. Depolarization-assisted induction of afterdischarge by repetitive stimulation of the hippocampal mossy fiber model. (A) Effect of repetitive stimulation of 50 times at 50 Hz at a resting potential of −80 mV. Membrane potential (V_m) at the 10th bouton is shown. (B) When the distal portion of the axon (the 10th bouton and the axons on both sides, shown in green in the schematic drawing) was depolarized to −70 mV, the same repetitive stimulation evoked afterdischarge persistent after the termination of the repetitive stimulus, as shown in the blue trace and an asterisk (*). (C) Depolarization-assisted afterdischarge was not induced with the stimulus of 40 times at 50 Hz with the same local depolarization. (D) The stimulus at 20 Hz applied 60 times, with slight depolarization of the distal axon to −70 mV also induced persistent afterdischarge (blue trace, **). (E) The stimuli at 20 Hz applied 50 times did not induce burst firings. (F) Stimulus-induced afterdischrge was also induced by 100 Hz stimuli 50 times (blue trace, ***). (G) The stimulus at 100 Hz applied 40 times is insufficient for generating afterdischarge. (H) Effect of 100 Hz stimuli applied 50 times to the model whose g_Na value was increased from the original value of 50 mS cm⁻² to 60 mS cm⁻². All the results suggest that local small depolarization is sufficient for assisting in boosting regenerative afterdischarge.

Results

Depolarization-assisted induction of afterdischarge by repetitive stimulation of the mossy fiber model

To test for the possible induction of stimulus-induced hyperexcitability of the axon, a repetitive stimulus of 50 times at 50 Hz was applied to the model of mossy fiber axons. In the normal condition at the resting membrane potentials at −80 mV (Engel and Jonas, 2005), repetitive stimulation of 50 times at 50 Hz did not induce persistent discharges after the termination of the repetitive stimulus (Figure 1A). In contrast, with slight depolarization of the distal part of the axon (the 10th bouton to −70 mV), the same repetitive stimulus elicited persistent firings after cessation of the repetitive stimulus (Figure 1B). The frequency of the afterdischarge induced by 50 Hz 50 times stimuli was 15.8 Hz. This stimulus-induced burst shows a regenerative and paroxysmal nature in its generation (Coggan et al., 2010) since a shorter repetitive stimulus of 40 times at 50 Hz did not elicit the burst discharge (Figure 1C). Even when the stimulus frequencies are changed to 20 Hz (Figures 1D,E) and 100 Hz (Figures 1F,G), similar results were obtained, i.e., repetitive stimuli induced persistent firings after the termination of the stimulus train above the threshold stimulus number. For example, 60 times at 20 Hz stimuli with depolarization to −70 mV induced the stimulus-induced afterdischarge of 14.7 Hz (Figure 1D), and 50 times at 100 Hz with the same depolarization elicited afterdischarge of 16.1 Hz (Figure 1F). The frequency of afterdischarge did not vary much among the stimulus frequencies of 20, 50, and 100 Hz. Shorter repetitive stimuli of 50 times at 20 Hz (Figure 1E) or 40 times at 100 Hz (Figure 1G) did not induce persistent firings after the termination of the repetitive stimulus. Notably, the action potentials during the 100 Hz stimulus were suppressed slightly at the initial phase and showed irregularity of the amplitudes at the later phase (green bars in Figures 1F,G). The initial decrease may be caused by incomplete recovery from the inactivation of Na⁺ channels since lower frequency stimulus (20–50 Hz) did not cause the depression of action potentials. The irregular responses at the late phase are caused by propagation failures that appear only at the late phase of the repetitive stimuli. One possible explanation is that accumulated fractional inactivation of sodium channels or insufficient recovery from inactivation of voltage-dependent sodium channels at short intervals (10 ms) reduced the available sodium channels and caused propagation failures. In support of this notion, it was found that the propagation failure was abolished if the g_Na value was increased from the original value of 50 mS cm⁻² to 60 mS cm⁻², as shown in Figure 1H. These findings suggest that a small local depolarization is sufficient to boost regenerative afterdischarge originating from an ectopic spike initiation site in the axon.

Antidromic propagation of the stimulus-induced afterdischarge of the mossy fiber model

Since the depolarization was applied to the distal portion of the axon in the previous simulation, it was hypothesized that the afterdischarge was triggered from the depolarized part of the distal axon and propagated antidromically to the soma. To test this hypothesis, the membrane potentials at the soma (black traces) and the 10th bouton (red traces) were calculated (Figure 2A). Both at the first (*) and the last action potentials (**) during the repetitive stimulus of 50 times at 50 Hz, somatic action potentials (black traces) were followed by those at the 10th bouton (red trace), indicating that action potentials were generated near the soma and propagated orthodontically toward the distal axons (Figure 2B, left and middle). In contrast, during the stimulus-induced persistent afterdischarge (***) following the repetitive stimulation, however, the order was reversed as red to black sequence (Figure 2B, right), as expected for the ectopic site for the burst firings from the depolarized portions of the axon and subsequent antidromic propagation toward the soma. Antidromic propagation of the stimulus-induced burst was also supported by the different action potential kinetics during and after the repetitive stimulus (Figures 2C,D). For traces in Figures 2C,D, simulations were performed at a faster sampling rate of 100 kHz to accurately monitor the kinetics of membrane potential (V_m) and its time derivative (dV_m/dt). The action potential kinetics was slower during the stimulus-induced afterdischarge. The half-duration of the somatic action potentials at the first (*) and the last (**) of the repetitive stimulus was 1.11 ms and 2.69 ms, respectively, while it was 2.40 ms after the repetitive stimulus for 50 times at 50 Hz. A slower action potential kinetics was also observed for the 10th bouton (Figure 2D). The half-duration of the action potentials at the 10th bouton during the same repetitive stimulus was 0.86 ms and 1.69 ms, respectively, while it was 1.53 ms after the repetitive stimulus. Taken together, all the evidence so far supported the generation of stimulus-induced afterdischarge from the ectopic site at the distal axons whose membrane potentials were slightly depolarized. The ectopic afterdischarge was propagated antidromically back to the soma. It should be noted that it is evident that the inflection of the waveform of the first derivative of the somatic membrane potential (dV_m/dt) exists in the rising phase. This inflection (asterisk in Figure 2C) may represent IS-SD inflection (Eccles et al., 1966; Stasheff et al., 1993), which may be caused by the impedance mismatch of the small initial segment (IS) and a large volume of somatodendritic compartment (SD) and thus supposed to indicate the antidromic invasion of axonal spikes to the soma after the generation of the ectopic afterdischarge from the distal axon. All the findings are in support of the antidromic propagation of the stimulus-induced afterdischarge of the mossy fiber model.

Figure 2

Figure 2. Ectopic spiking of the distal axon during depolarization-assisted stimulus-induced afterdischarge. (A) Stimulus-induced afterdischarge induced by repetitive stimulation of 50 times at 50 Hz with slight depolarization of distal axon to −70 mV, as shown in Figure 1B. Membrane potentials (V_m) at the soma and the 10th bouton are shown in black and red traces, respectively. (B) The time-expanded traces of somatic (black trace) and axonal action potentials (red trace) at the first (*; left), the last (**; middle), and after (***; right) the repetitive stimulus. (C) Comparison of the action potential kinetics (V_m, upper traces) and the time derivative of V_m (dV_m/dt, lower traces) at the soma (black traces). The asterisk (*) shows IS-SD inflections characteristic of antidromic spikes. (D) The similar data set of V_m at the 10th bouton (red traces) as in panel (C). The data for (C) and (D) are collected by the simulations with a faster sampling rate at 100 kHz to trace dV_m/dt with sufficient temporal resolution. These data suggest ectopic origin and subsequent antidromic propagation of stimulus-induced afterdischarge.

Inactivating the property of K⁺ channels assists the stimulus-induced afterdischarge of the mossy fiber model

Action potentials propagating along the hippocampal mossy fibers display use-dependent modification of the duration by slowing the falling phase (Geiger and Jonas, 2000). This unique property of axonal spike signaling imparts additional computing ability for fine-tuning the network function. As for the mechanisms, it was demonstrated that accumulated inactivation of K⁺ channels during the repetitive action potentials was critically involved (Geiger and Jonas, 2000; see also Zheng and Kamiya, 2023). It was inferred by analogy that the inactivation of K⁺ channels may underlie the generation of the stimulus-induced afterdischarge of the hippocampal mossy fibers. To test for the contribution of the inactivating property of axonal K⁺ channels, the effects of removal of inactivation from the axonal K⁺ channel model were explored (Zheng and Kamiya, 2023). A repetitive stimulation applied 50 times at 50 Hz induced the stimulus-induced afterdischarge, which outlasted the repetitive stimulus (Figure 3A). During the repetitive stimulation, action potentials both at the soma (black traces) and the 10th bouton (red traces) prolonged the kinetics by slowing the falling phase (Figure 3B). The half-duration of the first somatic action potential during the repetitive stimulus was 1.10 ms (black thin trace), while the half-duration of the last (50th) action potential was 2.68 ms (black thick trace). At the 10th bouton, the half-duration of the first action potentials during the repetitive stimulus was 0.85 ms (red thin trace), while the half-duration of the 50th action potential was 1.68 ms (red thick trace). These findings confirmed the results of our previous exploration (Zheng and Kamiya, 2023).

Figure 3

Figure 3. Inactivation of K⁺ channels assists induction of the stimulus-induced afterdischarge. (A) Simulation in the same conditions as in Figure 2A. (B) The time-expanded traces of action potentials (V_m) at the soma (black trace) and distal axon (red trace). The first (thin traces) and the last action potentials (thick traces) during the repetitive stimulus are superimposed. (C) Similar data as in panel (A) except that the model of K⁺ channels was replaced with those removed inactivation (non-inactivating K⁺ channel model, blue trace). (D) Similar data as in panel (B) with the non-inactivating type K⁺ channel model. (E,F) The model of K⁺ channels was replaced with a mixture of 90% inactivating type and 10% non-inactivating type and displayed as (A) and (B). (G,H) Similar data as in panels (E) and (F) when the distal portion of the axon was depolarized to −60 mV. (I,J) Simulation with longer repetitive stimulus (50 Hz 100 times) with strong depolarization of the distal portion to −60 mV with the non-inactivating K⁺ channel model. All the results suggest that the accumulated inactivation of K⁺ channels during the repetitive stimulus assists the generation of the afterdischarge in conjunction with local axonal depolarization of the distal axon.

In contrast, if the model of axonal K⁺ channels was replaced with those without inactivation (Zheng and Kamiya, 2023), the stimulus-induced burst was not induced (Figure 3C). In this condition, action potentials do not display use-dependent broadening during the repetitive stimulation of 50 times at 50 Hz (Figure 3D). The half-duration of the first somatic action potential during the repetitive stimulus was 1.10 ms (blue thin trace), while the half-duration of the 50th action potential was 1.27 ms (blue thick trace). At the 10th bouton, the half-duration of the first action potentials was 0.88 ms (red thin trace), while the half-duration for the 50th action potential was 0.90 ms during the repetitive stimulus (red thick trace). All the results are consistent with the notion that the accumulated inactivation of K⁺ channels during the repetitive stimulus assists the generation of the stimulus-induced afterdischarge of the mossy fiber model in conjunction with local small depolarization of the distal axon.

Although mossy fiber boutons predominantly express inactivating K⁺ channels (Geiger and Jonas, 2000), multiple types of K⁺ channels have been identified that shape action potentials, including some that do not inactivate (Alle et al., 2011). To investigate whether stimulus-induced afterdischarge could occur with a mixture of inactivating and non-inactivating channels, a model combining 90% inactivating and 10% non-inactivating channels was tested. The same depolarization of the distal portions to the −70 mV did not induce the afterdischarge (Figure 3E), although action potentials displayed a partial broadening of duration (Figure 3F). The half-duration of the first somatic action potential during the repetitive stimulus was 1.14 ms (black thin trace), while the half-duration of the last (50th) action potential was 1.91 ms (black thick trace). At the 10th bouton, the half-duration of the first action potential during the repetitive stimulus was 0.88 ms (red thin trace), while the half-duration of the 50th action potential was 1.14 ms (red thick trace). With this combination of inactivating and non-inactivating K⁺ channels, a stimulus-induced afterdischarge was elicited when the depolarization of the distal portions was increased to −60 mV, as shown in Figure 3G. The afterdischarge frequency was 31.7 Hz. The half-duration of the first somatic action potential during the repetitive stimulus was 1.14 ms (black thin trace), while the half-duration of the last (50th) action potential was 1.91 ms (black thick trace). At the 10th bouton, the half-duration of the first action potential during the repetitive stimulus was 0.88 ms (red thin trace), while the half-duration of the 50th action potential was 1.29 ms (red thick trace) (Figure 3H). These results confirmed the causal relationships between the inactivation properties of K⁺ channels and the generation of stimulus-induced afterdischarge.

One may argue that non-inactivating K⁺ channels may simply raise the threshold, but not be indispensable, for the generation of afterdischarge. To address this point, an additional simulation was performed to apply longer repetitive stimuli of 100 times at 50 Hz with the depolarization of distal portions to −60 mV. Even though this condition was expected to enhance the possibility of generating afterdischarge, the repetitive stimuli did not induce afterdischarge (Figure 3I) nor broadening of action potentials (Figure 3J). Taken together, the inactivation of axonal K⁺ channels is critically important for the stimulus-induced afterdischarge at the hippocampal mossy fibers.

Stimulus-induced ectopic afterdischarge is triggered along the course of mossy fibers

To date, the role of local depolarization of distal portions of the mossy fiber axons has been explored, as shown in Figure 4A. To further confirm the robustness of the notion that local depolarization assisted the generation of ectopic afterdischarge, the effects of the depolarization of the distal, middle, and proximal portions of mossy fiber axons and the uniform depolarization of the course of mossy fiber axons were explored. When the middle portion (the sixth boutons and both side axons) was depolarized to −70 mV, the stimulus-induced afterdischarge was not elicited (Figure 4B), while the same repetitive stimulation of 50 times at 50 Hz caused afterdischarge when depolarization was increased to −60 mV (Figure 4C). The frequency of the afterdischarge was 20.0 Hz. Similar results were obtained when the proximal portion (the second boutons and both side axons) was depolarized to −70 mV (Figure 4D) or to −60 mV (Figure 4E). The frequency of the afterdischarge was 15.5 Hz. It should be noted that the uniform depolarization of the axons (the 2nd–10th boutons and the axons on both sides) also elicited stimulus-induced afterdischarge with weaker depolarization to −75 mV (Figure 4F). The frequency of the afterdischarge was 19.2 Hz. These results suggest that depolarization-assisted generation of the stimulus-induced afterdischarge can be initiated from the whole axon.

Figure 4

Figure 4. Effects of the location of depolarization on stimulus-induced afterdischarge. (A) Simulation in the same conditions as in Figure 1B. (B,C) Responses to repetitive stimulation of 50 times at 50 Hz with slight depolarization of the middle portion of the axon (the 6th bouton and the axons on both sides) to −70 mV (B) and −60 mV (C). (D,E) Simulation with depolarization of the proximal portion of the axon (the 2nd bouton and the axons on both sides) to −70 mV (D) and −60 mV (E). (F) Simulation with uniform depolarization of the axon (the 2nd–10th boutons and the axons on both sides) to −75 mV. The results suggest that depolarization-assisted generation of the stimulus-induced afterdischarges initiated from the whole axon.

In this series of simulations, local depolarization was given to the portion of the axon by changing the equilibrium potential of the leak conductance. It should be noted that the steady-state level of depolarization did not reach the value since the depolarized compartment was connected to the non-depolarized compartment. To explore the threshold membrane potentials for the generation of afterdischarge, the steady-state membrane potentials at the depolarized compartment were calculated as −74.5 mV, −75.7 mV, −72.9 mV, −76.0 mV, −73.1 mV, and −75.3 mV for the simulation conditions of Figures 4A–F, respectively. It might be intriguing to speculate that the threshold for generating afterdischarge lies near −75 mV, although this is a too simplified notion not considering complex spatiotemporal dynamics of the axonal excitability. More importantly, it implicates the physiological significance that local small depolarization may cause a drastic impact to shift to the non-canonical mode of axonal signaling.

Bi-directional propagation of the stimulus-induced burst of the mossy fiber model

To date, the effects of local small depolarization in conjunction with the repetitive stimulus of the granule cells were examined, and it was found that the stimulus-induced burst discharges with ectopic origin were triggered after the repetitive stimulus and propagated antidromically to the soma. To explore the possible orthodromic propagation of the ectopic burst, repetitive stimulus for 50 times at 50 Hz in conjunction with depolarization to −70 mV of the sixth bouton and the axons on both sides did not induce the afterdischarge by the repetitive stimulus, as shown in Figure 5A. In contrast, the same repetitive stimulus with depolarization of the same compartments to −60 mV triggered the afterdischarge by the repetitive stimulus, as shown in Figure 5B. To determine the order of action potentials at the soma, the 6th bouton, and the 10th bouton, time-expanded traces were analyzed, as shown in Figure 5C. During the stimulus train (left and middle traces, * and **), the action potential at the soma (black trace) was followed by the action potential at the 6th (blue trace) and the 10th boutons (red trace), consistent with the orthodromic propagation from the soma to the distal axon. On the other hand, during the stimulus-induced afterdischarge, the order of these action potentials changed. The action potential at the 6th bouton where local depolarization was applied (blue trace) was followed by the action potentials at the 10th bouton (red trace) and the soma (black trace), as shown in the right traces of Figure 5C (***). It should be noted that the action potential kinetics at the soma were different during and after the repetitive stimulus (Figure 5D). The half-duration of the somatic action potentials at the first (*) and the last (**) of the repetitive stimulus was 1.11 ms (*) and 2.68 ms (**), respectively. After 50 repetitions at 50 Hz (***), the half-duration was 2.44 ms after the repetitive stimulus. A slower action potential kinetics was also observed for the 6th (Figure 5E) and the 10th boutons (Figure 5F). The half-duration of the action potentials at the 6th bouton at the first (*) and the last (**) of the repetitive stimulus was 0.85 ms and 1.85 ms, respectively, while it was 1.84 ms after the repetitive stimulus (**). The half-duration of the action potentials at the 10th bouton at the first (*) and the last (**) of the repetitive stimulus was 0.88 ms and 1.73 ms, respectively, while it was 1.62 ms after the repetitive stimulus (**). The different action potential kinetics during and after the repetitive stimulus were consistent with the notion that the stimulus-induced afterdischarge by the repetitive stimuli was triggered ectopically at different sites from the physiological spike initiation and propagated bi-directionally following the ectopic afterdischarge generation at the depolarized portions.

Figure 5

Figure 5. Bi-directional propagation of the afterdischarge from the site of local depolarization. (A) Effect of repetitive stimulation of 50 times at 50 Hz with depolarization of the 6th bouton and the axons on both sides to −70 mV. Membrane potentials (V_m) at the soma, the 6th bouton, and the 10th bouton are shown in black, blue, and red traces, respectively. (B) When the 6th bouton and the axons on both sides, shown in green in the schematic drawing, were depolarized to −60 mV, stimulus-induced afterdischarge was induced after cessation of the repeated stimulus. (C) The time-expanded traces of action potentials at the soma (black trace), the 6th bouton (blue trace), and the 10th bouton (red trace) at the first (*; left), the last (**; middle), and after (***; right) the repetitive stimulus. (D) Comparison of the action potential kinetics (V_m, upper traces) and the time derivative of V_m (dVm/dt, lower traces) at the soma (black traces). (E,F) The similar data set of V_m at the 6th bouton (blue traces) and 10th bouton (red traces) as in panel (D). The data for (D–F) are collected by the simulations with a faster sampling rate at 100 kHz to trace dV_m/dt with sufficient temporal resolution. The results suggest bi-directional propagation of the afterdischarge from the site of initiation.

Discussion

In this study, a numerical simulation approach was adopted to examine the mechanisms underlying the stimulus-induced afterdischarge at the hippocampal mossy fiber, which does not substantially express HCN channels, key molecules for this form of plasticity of excitability of axons. A repetitive stimulus, in conjunction with a small depolarization of distal portions of the mossy fiber, triggered a prolonged afterdischarge in the simple biophysical model of mossy fiber. The stimulus-induced burst firings were supposed to be generated at ectopic sites of spike initiation different from the physiological spike initiation sites, due to the reversed order of action potentials at the soma and the distal axon and the different waveform of action potentials during and after the repetitive stimulus. Depolarization-assisted stimulus-induced bursts would be an additional mechanism for the retroaxonal barrage firings observed experimentally.

Stimulus-induced burst firings induced in the mossy fiber model

Recent studies have shown that strong repetitive stimulation occasionally triggers afterdischarge persistent even after the cessation of the repetitive stimulus in some neurons. Some cortical inhibitory neurons have been shown to exhibit robust plasticity in the intrinsic excitability of their axons (Sheffield et al., 2011). This form of plasticity of the axonal excitability called retroaxonal barrage firings is facilitated by the HCN channels on the axons (Elgueta et al., 2015; Rózsa et al., 2023; see also Roth and Hu, 2020) and K⁺ accumulation by the repetitive stimulus (Smith, 1983; Meeks and Mennerick, 2004, 2007b; Foley et al., 2010; Beswick-Jones et al., 2023; Rózsa et al., 2023) or surrounding glia (Meeks et al., 2005; Deemyad et al., 2018). Repeated stimulation of the mossy fiber has been demonstrated to induce the stimulus-induced prolonged repeated firings (Kamiya, 2017), although the mechanisms remain uncertain since the hippocampal mossy fiber axons do not substantially express HCN channel subtypes HCN1-4, which has been shown to mediate the induction of burst discharges (Notomi and Shigemoto, 2004). A recent study has indicated that a local increase in potassium concentration also promoted the induction of stimulus-induced afterdischarge in the cortical inhibitory interneurons (Rózsa et al., 2023). In analogy with this observation, I wondered if local subthreshold depolarization, assuming a local increase in potassium concentration, can assist the induction of the bursts. As expected from this notion, it was confirmed that local depolarization of the axonal portion to −70 mV from the resting membrane potential of −80 mV can generate stimulus-induced afterdischarge, as shown in Figure 1. The paroxysmal nature was demonstrated by the finding that the stimulus numbers above the threshold are needed for the induction of the stimulus-induced burst, as also shown in Figure 1.

Ectopic origin of stimulus-induced burst firings from distal axon

Regarding the nature of the stimulus-induced afterdischarge, it was noted that the prolonged repeated discharges were generated from the distal axon where small depolarization was applied. This contrasts with the action potential initiation from the proximal site of the mossy fiber (Schmidt-Hieber et al., 2008) in the physiological condition. Emerging evidence suggests that the action potentials are generated not only from the proximal axon or the axon initial segments but also from the distal axon ectopically (Stasheff et al., 1993; Avoli et al., 1998; Sheffield et al., 2011; Dugladze et al., 2012). As in these experimental findings, the stimulus-induced afterdischarge was triggered with the assistance of depolarization of the distal portion of the axon. The ectopic origin of the afterdischarge in the hippocampal mossy fibers is supported by several lines of the findings of this study. First, the stimulus-induced bursts are triggered only when small depolarization is applied to distal portions of the axon which is far from the physiological spike initiation site at the proximal axon (Schmidt-Hieber et al., 2008). Second, the action potentials at the distal axon were followed by those at the soma during the bursts, indicating antidromic propagation of the burst firings, as shown in Figure 1. Third, the somatic action potential during the stimulus-induced bursts displayed the inflection on the rising phase, resembling the IS-SD inflection characteristic for invasion of antidromic propagation to the soma (Eccles et al., 1966; Stasheff et al., 1993), as shown in Figure 2C. These findings are all consistent with the ectopic origin of the stimulus-induced afterdischarge and the subsequent antidromic propagation toward the soma.

Since the simulation also indicates that the ectopically triggered afterdischarge propagates bi-directionally from the site of burst generation, as shown in Figure 5, the stimulus-induced ectopic bursts are supposed to cause the massive transmitter release and a barrage of postsynaptic currents on the postsynaptic CA3 pyramidal neurons and lead to raising the excitatory tones in the CA3 neuronal network (Suzuki et al., 2014). On the other hand, however, antidromic invasion of bursting action potentials may collide with and therefore suppress output signals from the soma. Synaptic plasticity caused by the stimulus-induced bursts (Bukalo et al., 2013) also needs to be considered for the net influence of the burst firings. Future investigations should determine the spatiotemporal consequences of stimulus-induced bursts on activity within the CA3 neuronal network.

Use-dependent broadening of action potentials assists induction of the stimulus-induced bursts

A characteristic feature of hippocampal mossy fibers is the activity-dependent broadening of action potentials, caused by the progressive slowing of the falling phase due to the cumulative inactivation of voltage-dependent potassium channels (Geiger and Jonas, 2000; see also Zheng and Kamiya, 2023). By removing the inactivation of the axonal K⁺ channels in models of computer simulations, the depolarization-assisted stimulus-induced burst firings were not induced in the same repetitive stimulation, as shown in Figure 4. It would be interesting to address this notion experimentally, although no pharmacological tools to selectively and completely suppress the inactivation of voltage-dependent K⁺ channels are available currently. Future studies are needed to test this idea in experiments.

In this study, a series of numerical simulations using a simple model of hippocampal mossy fiber was performed to test for the possible induction of the stimulus-induced burst firings persisting even after the termination of the repetitive stimulus. The simulation demonstrated that the repetitive stimulation with the aid of local small depolarization of the axon generates regenerative burst firings in the model. Ectopic generation of afterdischarge was implicated by the signs of antidromic propagation including reversed order of somatic and axonal action potentials and IS-SD inflection at the rising phase of the somatic action potentials. Use-dependent broadening of action potentials by accumulated inactivation of K⁺ channels has been suggested to mediate the generation of the stimulus-induced afterdischarge. Taking all pieces of evidence into account, these mechanisms are suggested to promote this form of plasticity of axonal excitability of hippocampal mossy fiber axons.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors without undue reservation.

Author contributions

HK: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by a Grant-in-Aid for Scientific Research (KAKENHI) from the Japan Society for the Promotion of Science (21K06434 and 24H02336 to HK).

Conflict of interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

1. ^https://modeldb.science/2018003

References

Acsády, L., Kamondi, A., Sík, A., Freund, T., and Buzsáki, G. (1998). GABAergic cells are the major postsynaptic targets of mossy fibers in the rat hippocampus. J. Neurosci. 18, 3386–3403. doi: 10.1523/JNEUROSCI.18-09-03386.1998

PubMed Abstract | Crossref Full Text | Google Scholar

Alle, H., Kubota, H., and Geiger, J. R. (2011). Sparse but highly efficient Kv3 outpace BKCa channels in action potential repolarization at hippocampal mossy fiber boutons. J. Neurosci. 31, 8001–8012. doi: 10.1523/JNEUROSCI.0972-11.2011

PubMed Abstract | Crossref Full Text | Google Scholar

Alpizar, S. A., Cho, I. H., and Hoppa, M. B. (2019). Subcellular control of membrane excitability in the axon. Curr. Opin. Neurobiol. 57, 117–125. doi: 10.1016/j.conb.2019.01.020

PubMed Abstract | Crossref Full Text | Google Scholar

Amir, R., Kocsis, J. D., and Devor, M. (2005). Multiple interacting sites of ectopic spike electrogenesis in primary sensory neurons. J. Neurosci. 25, 2576–2585. doi: 10.1523/JNEUROSCI.4118-04.2005

Crossref Full Text | Google Scholar

Avoli, M., Methot, M., and Kawasaki, H. (1998). GABA-dependent generation of ectopic action potentials in the rat hippocampus. Eur. J. Neurosci. 10, 2714–2722

PubMed Abstract | Google Scholar

Bähner, F., Weiss, E. K., Birke, G., Maier, N., Schmitz, D., Rudolph, U., et al. (2011). Cellular correlate of assembly formation in oscillating hippocampal networks in vitro. Proc. Natl. Acad. Sci. USA 108, E607–E616. doi: 10.1073/pnas.1103546108

PubMed Abstract | Crossref Full Text | Google Scholar

Beswick-Jones, H., Hopper, A. J., and Brown, A. M. (2023). Activity-dependent fluctuations in interstitial [K⁺]: investigations using ion-sensitive microelectrodes. Molecules 28:523. doi: 10.3390/molecules28020523

PubMed Abstract | Crossref Full Text | Google Scholar

Bucher, D., and Goaillard, J. M. (2011). Beyond faithful conduction: short-term dynamics, neuromodulation, and long-term regulation of spike propagation in the axon. Prog. Neurobiol. 94, 307–346. doi: 10.1016/j.pneurobio.2011.06.001

PubMed Abstract | Crossref Full Text | Google Scholar

Bukalo, O., Campanac, E., Hoffman, D. A., and Fields, R. D. (2013). Synaptic plasticity by antidromic firing during hippocampal network oscillations. Proc. Natl. Acad. Sci. USA 110, 5175–5180. doi: 10.1073/pnas.1210735110

PubMed Abstract | Crossref Full Text | Google Scholar

Coggan, J. S., Prescott, S. A., Bartol, T. M., and Sejnowski, T. J. (2010). Imbalance of ionic conductances contributes to diverse symptoms of demyelination. Proc. Natl. Acad. Sci. USA 107, 20602–20609. doi: 10.1073/pnas.1013798107

PubMed Abstract | Crossref Full Text | Google Scholar

Deemyad, T., Lüthi, J., and Spruston, N. (2018). Astrocytes integrate and drive action potential firing in inhibitory subnetworks. Nat. Commun. 9:4336. doi: 10.1038/s41467-018-06338-3

PubMed Abstract | Crossref Full Text | Google Scholar

Dugladze, T., Schmitz, D., Whittington, M. A., Vida, I., and Gloveli, T. (2012). Segregation of axonal and somatic activity during fast network oscillations. Science 336, 1458–1461. doi: 10.1126/science.1222017

PubMed Abstract | Crossref Full Text | Google Scholar

Eccles, J. C., Llinás, R., and Sasaki, K. (1966). The action of antidromic impulses on the cerebellar Purkinje cells. J. Physiol. 182, 316–345. doi: 10.1113/jphysiol.1966.sp007826

PubMed Abstract | Crossref Full Text | Google Scholar

Elgueta, C., Köhler, J., and Bartos, M. (2015). Persistent discharges in dentate gyrus perisoma-inhibiting interneurons require hyperpolarization-activated cyclic nucleotide-gated channel activation. J. Neurosci. 35, 4131–4139. doi: 10.1523/JNEUROSCI.3671-14.2015

PubMed Abstract | Crossref Full Text | Google Scholar

Engel, D., and Jonas, P. (2005). Presynaptic action potential amplification by voltage-gated Na⁺ channels in hippocampal mossy fiber boutons. Neuron 45, 405–417. doi: 10.1016/j.neuron.2004.12.048

PubMed Abstract | Crossref Full Text | Google Scholar

Foley, J., Nguyen, H., Bennett, C. B., and Muschol, M. (2010). Potassium accumulation as dynamic modulator of neurohypophysial excitability. Neuroscience 169, 65–73. doi: 10.1016/j.neuroscience.2010.04.049

PubMed Abstract | Crossref Full Text | Google Scholar

Geiger, J. R., and Jonas, P. (2000). Dynamic control of presynaptic Ca²⁺ inflow by fast-inactivating K⁺ channels in hippocampal mossy fiber boutons. Neuron 28, 927–939. doi: 10.1016/S0896-6273(00)00164-1

PubMed Abstract | Crossref Full Text | Google Scholar

Gutnick, M. J., and Prince, D. A. (1972). Thalamocortical relay neurons: antidromic invasion of spikes from a cortical epileptogenic focus. Science 176, 424–426. doi: 10.1126/science.176.4033.424

PubMed Abstract | Crossref Full Text | Google Scholar

Hallermann, S., Pawlu, C., Jonas, P., and Heckmann, M. (2003). A large pool of releasable vesicles in a cortical glutamatergic synapse. Proc. Natl. Acad. Sci. USA 100, 8975–8980. doi: 10.1073/pnas.1432836100

PubMed Abstract | Crossref Full Text | Google Scholar

Henze, D. A., Urban, N. N., and Barrionuevo, G. (2000). The multifarious hippocampal mossy fiber pathway: a review. Neuroscience 98, 407–427. doi: 10.1016/S0306-4522(00)00146-9

PubMed Abstract | Crossref Full Text | Google Scholar

Hines, M. L., and Carnevale, N. T. (1997). The NEURON simulation environment. Neural Comput. 9, 1179–1209. doi: 10.1162/neco.1997.9.6.1179

Crossref Full Text | Google Scholar

Kamiya, H. (2017). Prolonged hyperexcitability of hippocampal mossy fibers after repetitive high frequency stimulation. Soc. Neurosci. Abstract 474:05.

Google Scholar

Kamiya, H. (2019). Modeling analysis of axonal after potential at hippocampal mossy fibers. Front. Cell. Neurosci. 13:210. doi: 10.3389/fncel.2019.00210

PubMed Abstract | Crossref Full Text | Google Scholar

Kapoor, R., Li, Y. G., and Smith, K. J. (1997). Slow sodium-dependent potential oscillations contribute to ectopic firing in mammalian demyelinated axons. Brain 120, 647–652. doi: 10.1093/brain/120.4.647

PubMed Abstract | Crossref Full Text | Google Scholar

Kramer, P. F., Brill-Weil, S. G., Cummins, A. C., Zhang, R., Camacho-Hernandez, G. A., Newman, A. H., et al. (2022). Synaptic-like axo-axonal transmission from striatal cholinergic interneurons onto dopaminergic fibers. Neuron 110, 2949–2960.e4. doi: 10.1016/j.neuron.2022.07.011

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, C., Cai, X., Ritzau-Jost, A., Kramer, P. F., Li, Y., Khaliq, Z. M., et al. (2022). An action potential initiation mechanism in distal axons for the control of dopamine release. Science 375, 1378–1385. doi: 10.1126/science.abn0532

PubMed Abstract | Crossref Full Text | Google Scholar

Meeks, J. P., Jiang, X., and Mennerick, S. (2005). Action potential fidelity during normal and epileptiform activity in paired soma-axon recordings from rat hippocampus. J. Physiol. 566, 425–441. doi: 10.1113/jphysiol.2005.089086

PubMed Abstract | Crossref Full Text | Google Scholar

Meeks, J. P., and Mennerick, S. (2004). Selective effects of potassium elevations on glutamate signaling and action potential conduction in hippocampus. J. Neurosci. 24, 197–206. doi: 10.1523/JNEUROSCI.4845-03.2004

PubMed Abstract | Crossref Full Text | Google Scholar

Meeks, J. P., and Mennerick, S. (2007a). Action potential initiation and propagation in CA3 pyramidal axons. J. Neurophysiol. 97, 3460–3472. doi: 10.1152/jn.01288.2006

PubMed Abstract | Crossref Full Text | Google Scholar

Meeks, J. P., and Mennerick, S. (2007b). Astrocyte membrane responses and potassium accumulation during neuronal activity. Hippocampus 17, 1100–1108. doi: 10.1002/hipo.20344

PubMed Abstract | Crossref Full Text | Google Scholar

Notomi, T., and Shigemoto, R. (2004). Immunohistochemical localization of Ih channel subunits, HCN1-4, in the rat brain. J. Comp. Neurol. 471, 241–276. doi: 10.1002/cne.11039

PubMed Abstract | Crossref Full Text | Google Scholar

Ohura, S., and Kamiya, H. (2016). Excitability tuning of axons in the central nervous system. J. Physiol. Sci. 66, 189–196. doi: 10.1007/s12576-015-0415-2

PubMed Abstract | Crossref Full Text | Google Scholar

Ohura, S., and Kamiya, H. (2018). Short-term depression of axonal spikes at the mouse hippocampal mossy fibers and sodium channel-dependent modulation. eNeuro 5:ENEURO.0415-17.2018. doi: 10.1523/ENEURO.0415-17.2018

Crossref Full Text | Google Scholar

Pinault, D. (1995). Backpropagation of action potentials generated at ectopic axonal loci: hypothesis that axon terminals integrate local environmental signals. Brain Res. Brain Res. Rev. 21, 42–92. doi: 10.1016/0165-0173(95)00004-M

PubMed Abstract | Crossref Full Text | Google Scholar

Roth, F. C., and Hu, H. (2020). An axon-specific expression of HCN channels catalyzes fast action potential signaling in GABAergic interneurons. Nat. Commun. 11:2248. doi: 10.1038/s41467-020-15791-y

PubMed Abstract | Crossref Full Text | Google Scholar

Rózsa, M., Tóth, M., Oláh, G., Baka, J., Lákovics, R., Barzó, P., et al. (2023). Temporal disparity of action potentials triggered in axon initial segments and distal axons in the neocortex. Sci. Adv. 9:eade4511. doi: 10.1126/sciadv.ade4511

PubMed Abstract | Crossref Full Text | Google Scholar

Schmidt-Hieber, C., Jonas, P., and Bischofberger, J. (2008). Action potential initiation and propagation in hippocampal mossy fibre axons. J. Physiol. 586, 1849–1857. doi: 10.1113/jphysiol.2007.150151

PubMed Abstract | Crossref Full Text | Google Scholar

Sheffield, M. E., Best, T. K., Mensh, B. D., Kath, W. L., and Spruston, N. (2011). Slow integration leads to persistent action potential firing in distal axons of coupled interneurons. Nat. Neurosci. 14, 200–207. doi: 10.1038/nn.2728

PubMed Abstract | Crossref Full Text | Google Scholar

Smith, D. O. (1983). Extracellular potassium levels and axon excitability during repetitive action potentials in crayfish. J. Physiol. 336, 143–157. doi: 10.1113/jphysiol.1983.sp014573

PubMed Abstract | Crossref Full Text | Google Scholar

Stasheff, S. F., Hines, M., and Wilson, W. A. (1993). Axon terminal hyperexcitability associated with epileptogenesis in vitro. I. Origin of ectopic spikes. J. Neurophysiol. 70, 961–975. doi: 10.1152/jn.1993.70.3.961

PubMed Abstract | Crossref Full Text | Google Scholar

Stasheff, S. F., and Wilson, W. A. (1990). Increased ectopic action potential generation accompanies epileptogenesis in vitro. Neurosci. Lett. 111, 144–150. doi: 10.1016/0304-3940(90)90359-H

PubMed Abstract | Crossref Full Text | Google Scholar

Suzuki, N., Tang, C. S., and Bekkers, J. M. (2014). Persistent barrage firing in cortical interneurons can be induced in vivo and may be important for the suppression of epileptiform activity. Front. Cell. Neurosci. 8:76. doi: 10.3389/fncel.2014.00076

PubMed Abstract | Crossref Full Text | Google Scholar

Thome, C., Roth, F. C., Obermayer, J., Yanez, A., Draguhn, A., and Egorov, A. V. (2018). Synaptic entrainment of ectopic action potential generation in hippocampal pyramidal neurons. J. Physiol. 596, 5237–5249. doi: 10.1113/JP276720

PubMed Abstract | Crossref Full Text | Google Scholar

Wissmann, R., Bildl, W., Oliver, D., Beyermann, M., Kalbitzer, H. R., Bentrop, D., et al. (2003). Solution structure and function of the "tandem inactivation domain" of the neuronal A-type potassium channel Kv1.4. J. Biol. Chem. 278, 16142–16150. doi: 10.1074/jbc.M210191200

PubMed Abstract | Crossref Full Text | Google Scholar

Zheng, F., and Kamiya, H. (2023). Simulation test for impartment of use-dependent plasticity by inactivation of axonal potassium channels on hippocampal mossy fibers. Front. Cell. Neurosci. 17:1154910. doi: 10.3389/fncel.2023.1154910

PubMed Abstract | Crossref Full Text | Google Scholar