QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

HOME  |   SEARCH  |   ARCHIVE  |   SUBSCRIBE  |   CONTACT  |   HELP

This Article Abstract Full Text (PDF) Submit an eLetter Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (49) Citing Articles via Google Scholar Google Scholar Articles by D'Angelo, E. Articles by Naldi, G. Search for Related Content PubMed PubMed Citation Articles by D'Angelo, E. Articles by Naldi, G.

 Previous Article  |  Next Article 

The Journal of Neuroscience, February 1, 2001, 21(3):759-770

Theta-Frequency Bursting and Resonance in Cerebellar Granule Cells: Experimental Evidence and Modeling of a Slow K⁺-Dependent Mechanism

Egidio D'Angelo^{1, 2}, Thierry Nieus¹, Arianna Maffei¹, Simona Armano¹, Paola Rossi¹, Vanni Taglietti¹, Andrea Fontana³, and Giovanni Naldi⁴

¹  Department of Molecular/Cellular Physiology and Instituto Nazionale per la Fisica della Materia, University of Pavia, I-27100 Pavia, Italy, ²  Department of Evolutionary and Functional Biology, University of Parma, Parma, Italy, ³  Department of Nuclear and Theoretical Physics, University of Pavia, Pavia, Italy, and ⁴  Department of Mathematics and Applications, University of Milano Bicocca, Milan, Italy

	ABSTRACT

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Neurons process information in a highly nonlinear manner, generating oscillations, bursting, and resonance, enhancing responsiveness at preferential frequencies. It has been proposed that slow repolarizing currents could be responsible for both oscillation/burst termination and for high-pass filtering that causes resonance (Hutcheon and Yarom, 2000). However, different mechanisms, including electrotonic effects (Mainen and Sejinowski, 1996), the expression of resurgent currents (Raman and Bean, 1997), and network feedback, may also be important. In this study we report theta-frequency (3-12 Hz) bursting and resonance in rat cerebellar granule cells and show that these neurons express a previously unidentified slow repolarizing K⁺ current (I_K-slow). Our experimental and modeling results indicate that I_K-slow was necessary for both bursting and resonance. A persistent (and potentially a resurgent) Na⁺ current exerted complex amplifying actions on bursting and resonance, whereas electrotonic effects were excluded by the compact structure of the granule cell. Theta-frequency bursting and resonance in granule cells may play an important role in determining synchronization, rhythmicity, and learning in the cerebellum.

Key words: bursting; resonance; M-current; cerebellum; granule cell; modeling

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Neurons process information by generating action potentials organized either in regular discharges (fast repetitive firing) or in bursts (Connors and Gutnik, 1990), which can occur repetitively when they are sustained by slow membrane potential oscillations (Wang and Rinzel, 1999). Moreover, some neurons respond better to a preferential input frequency, a property called resonance (Hutcheon and Yarom, 2000). Oscillations, bursting, and resonance have been related to synchronization of neuronal activity and to the emergence of brain rhythms (Llinas, 1988).

Oscillations and bursting can arise from various mechanisms that involve slow depolarizing and repolarizing currents. Noteworthy examples are provided by the interaction of a persistent Na⁺ current (I_Na-p) with a slow K⁺ current (I_M) in hippocampal pyramidal neurons (Gutfreund et al., 1995; Pape and Driesang, 1998) or with an inward rectifier current (I_h) in entorhinal neurons (Alonso and Llinas, 1989; Dickson et al., 2000). In thalamic neurons, depolarization and delayed repolarization are determined by low-threshold Ca²⁺ current (I_T) activation and inactivation and are regulated by other currents, including I_h (McCormick and Huguenard, 1992). In invertebrate cells, Ca²⁺-dependent bursts are terminated by a Ca²⁺-dependent K⁺ current (I_AHP) (Wang and Rinzel, 1999). In addition, bursting emerges when a membrane potential difference is established between dendrites and soma, causing rebound depolarization after a spike (Mainen and Sejinowski, 1996). Rebound depolarization may also be caused by currents activated by spike repolarization, e.g., by a resurgent Na⁺ current (Raman and Bean, 1997). Despite the multitude of mechanisms that potentially generate oscillation and bursting, resonance essentially requires a slow repolarizing current that reduces neuronal excitability at low input frequency (Hutcheon et al., 1996a,b; Hutcheon and Yarom, 2000).

The most apparent discharge mode of cerebellar granule cells is fast repetitive firing (Gabbiani et al., 1994; D'Angelo et al., 1995). However, a more complex electrical behavior has been suggested by the observation of spike bursting unveiled by pharmacological manipulation (D'Angelo et al., 1998). In this study, we report that cerebellar granule cells also show theta-frequency resonance, and we suggest that both bursting and resonance are based on a slow K⁺ current. This conclusion is supported by a mathematical model that provides a realistic reconstruction of granule cell electroresponsiveness.

Theta-frequency bursting and resonance in granule cells may play an important role in determining synchronization (Maex and DeShutter, 1998), rhythmicity (Pellerin and Lamarre, 1997; Hartmann and Bower, 1998), and learning (D'Angelo et al., 1999; Armano et al., 2000) at the major input stage of the cerebellum.

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Whole-cell patch-clamp recordings. Cerebellar granule cells were recorded in acute cerebellar slices obtained from 20 ±2-d-old rats. Slice preparation and patch-clamp recordings were performed as reported previously (Rossi et al., 1994, 1998; D'Angelo et al., 1995, 1997, 1998, 1999; Armano et al., 2000).

Current-clamp recordings were performed at 30°C. The extracellular solution contained (in mM): NaCl 120, KCl 2, MgSO₄ 1.2, NaHCO₃ 26, KH₂PO₄ 1.2, CaCl₂ 2, glucose 11, and bicuculline 0.01, and was equilibrated with 95% O₂ and 5% CO₂, pH 7.4. The pipette solution contained (in mM): K-gluconate 126 (or Cs₂SO₄ 78), KCl 4 (or CsCl 4), NaCl 4, MgSO₄ 1, CaCl₂ 0.02, BAPTA 0.1, glucose 15, ATP 3, GTP 0.1, HEPES 5; pH was adjusted to 7.2 with KOH (or CsOH).

Voltage-clamp recordings were performed at room temperature (25.5°C). The extracellular solution contained (in mM): NaCl 100, KCl 2, KH₂PO₄ 1.2, MgSO₄ 1.2, NaHCO₃ 26, glucose 11, Tetraethyl-ammonium Cl (TEA) 20, 4-amino-piridine (4-AP) 4, Ni²⁺ 2, tetrodotoxin (TTX) 0.001, and bicuculline 0.01, and was equilibrated with 95% O₂ and 5% CO₂, pH 7.4. The pipette solution contained (in mM): K-gluconate 126, NaCl 4, MgSO₄ 1, CaCl₂ 0.02, BAPTA 0.1, glucose 15, ATP 3, GTP 0.1, HEPES 5; pH was adjusted to 7.2 with KOH.

TEA, 4-AP, TTX, and bicuculline were obtained from Sigma (St. Louis, MO). The glutamate receptor antagonists D-2-amino-5-phosphonovaleric acid (APV), 7-chlorokinurenic acid (7-Cl-Kyn), and 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX) were obtained from Tocris Cookson (Bristol, UK).

Data were recorded with an Axopatch 200B amplifier, digitized with a Digidata 1200 interface (500 µsec/point), and analyzed with PClamp software (Axon Instruments, Foster City, CA). In voltage-clamp recordings, leak subtraction was performed by using a P/4 protocol. All data are reported as mean ± SD.

Mathematical modeling. A mathematical model of rat cerebellar granule cell electroresponsiveness (D'Angelo et al., 1995, 1998; Brickley et al., 1996) was constructed using the NEURON simulator (Hines and Carnevale, 1997). Because granule cells have a compact electrotonic structure (Silver et al., 1992; D'Angelo et al., 1993, 1995), a single-compartment model was used. The experimental value of membrane capacitance (3 pF; see refences above) was used to calculate the granule cell surface, assuming a spherical shape and a specific membrane capacitance of 1 µF/cm².

Mathematical methods. The mathematical problem in neuronal simulation is to solve the set of differential equations representing membrane voltage, intracellular Ca²⁺ concentration, and channel gating dynamics [see for example Yamada et al. (1998)]. Voltage was obtained as the time integral of the equation:

(1)

where V is membrane potential, C_m is membrane capacitance, g_i is ionic conductance, V_i is reversal potential (the subscript i indicates different channels), and i_inj is the injected current. Membrane conductances were represented using Hodgkin-Huxley-like models (Hodgkin and Huxley, 1952) of the type:

(2)

where G_maxi is the maximum ionic conductance, x_i and y_i are state variables (probabilities ranging from 0 to 1) for a gating particle, and z_i is the number of such gating particles in ionic channel i. x and y (with the suffix i omitted) were related to the first-order rate constants  and  by the equations:

(3)

(4)

where  and  are functions of voltage. The equations used to parameterize  and  and the state variables x, _x, y, and _y for different ionic channels are shown in Table 1. The state variable kinetics were:

(5)

(6)

The model included a leakage current and voltage-dependent Na⁺, Ca²⁺, and K⁺ conductances (see Table 1). Nernst equilibrium potentials were calculated from ionic concentrations used in current-clamp recordings. The Ca²⁺ equilibrium potential was updated after changes in the intracellular Ca²⁺ concentration.

All ionic currents used in the model have been identified in cerebellar granule cells in situ, when the excitable response has assumed its mature pattern (>P20) (D'Angelo et al., 1997). Gating kinetics were corrected using a Q₁₀ = 3 according to the relation Q₁₀ ^{(Tsim-Texp)/10} (Gutfreund et al., 1995) to account for differences between simulation temperature (Tsim = 30°C) and experimental temperature (T_exp). Maximum ionic conductances were corrected for ionic concentration differences between voltage- and current-clamp recordings. A further adjustment (usually <30%) of current densities allowed us to fine tune the excitable response [for further explanations see Traub and Llinas (1979); Traub et al. (1991); Vanier and Bower (1999)].

Leakage current. Mature rat cerebellar granule cells in situ have an aspecific and a GABA-A receptor-dependent leakage (Brickley et al., 1996). In the model, leakage consisted of a 5.68 × 10⁵ S/cm² conductance with reversal potential at 59 mV (g_L), and of a 2.17 × 10⁵ S/cm² Cl conductance with reversal potential at 65 mV (accounting for 28% of the total input conductance) (Armano et al., 2000). No qualitative difference was observed in model responses by setting GABA-A receptor leakage to zero (data not shown).

Na⁺ currents (I_Na-f, I_Na-p, I_Na-r). Mature rat cerebellar granule cells in situ express a fast and persistent Na⁺ current (I_Na-f and I_Na-p) (D'Angelo et al., 1998) and probably also a resurgent Na⁺ current (I_Na-) (E. D'Angelo and J. Magistretti, unpublished observations). The I_Na-f model was based on Gutfreund et al. (1995), and inactivation was slowed down around threshold to reproduce spike adaptation during bursts (Mainen et al., 1995). I_Na-f density was set to reproduce repetitive firing. I_Na-p was reproduced from Gutfreund et al. (1995) and shifted by 2 mV to match spike threshold. I_Na-p density (1/65 I_Na-f) was set to reproduce Na⁺-dependent plateaus (D'Angelo et al., 1998). I_Na-r was reconstructed from Raman and Bean (1997). I_Na-r density (initially 1/26 I_Na-f) was regulated in different simulations.

Ca²⁺ current (I_Ca). The I_Ca model was derived from high-threshold Ca²⁺ currents (mostly N-type) measured in mature rat cerebellar granule cells in situ (Rossi et al., 1994). I_Ca had fast second-order activation kinetics and slow voltage-dependent inactivation. I_Ca density was halved to account for different extracellular Ca²⁺ concentrations.

K⁺ currents (I_K-V, I_K-Ca, I_K-A, I_K-slow, I_K-IR). Mature rat cerebellar granule cells in situ express I_K-V, I_K-Ca, and I_K-A (Cull-Candy et al., 1989; Bardoni and Belluzzi, 1993), I_K-IR (Rossi et al., 1998), and I_K-slow (this study). I_K-V is a voltage-dependent K⁺ current resembling other neuronal delayed rectifiers, and its model has been adapted from Gutfreund et al. (1995). I_K-V was shifted by 5 mV to match I_Na-f. I_K-Ca is a voltage- and Ca²⁺-dependent K⁺ current corresponding to "big-K" channel recordings from granule cells in culture (Fagni et al., 1991), the kinetics of which are largely determined by those of the associated Ca²⁺ channels and intracellular Ca²⁺ fluctuations. The I_K-Ca model is the same as in Gabbiani et al. (1994), and I_K-Ca density was set close to that of I_K-V. I_A is a fast-activating, fast-inactivating voltage-dependent K⁺ current, which was reproduced from data reported by Bardoni and Belluzzi (1993). I_K-slow is a slow Ca²⁺-independent TEA-insensitive K⁺ current, which was reconstructed from data shown in Figure 5. I_K-slow voltage dependence was shifted by 10 mV (8 mV accounting for liquid-junction potential and 2 mV required to maintain a proper matching with I_Na-p). I_K-IR is a fast inward rectifier current that was reconstructed from data reported by Rossi et al. (1998).

Ca²⁺ dynamics. The intracellular Ca²⁺ concentration, [Ca²⁺], was calculated through the equation:

(7)

where d is the depth of a shell adjacent to the cell surface of area A, _Ca determines the loss of Ca²⁺ ions from the shell approximating the effect of fluxes, ionic pumps, diffusion, and buffers (Traub and Llinas, 1979; McCormick and Huguenard, 1992; DeSchutter and Smolen, 1998), and Ca0 is resting [Ca²⁺]. Once I_Ca and I_KCa had been set, Ca²⁺ dynamics were adapted to yield Ca²⁺ transients of ~1 µM, similar to those reported by Gabbiani et al. (1994), and to reproduce the excitable response (Traub and Llinas, 1979; Traub et al., 1991). Parameters used in Equation 7 were d = 200 nm, _Ca=1.5, and Ca0 = 100 nM. Ca0 was measured in rat cerebellar granule cells in culture (Irving et al., 1992; Marchetti et al., 1995) and experimentally maintained by appropriate BAPTA-Ca²⁺ buffers (see above).

Resting membrane potential. Resting membrane potential in the model settled at 80 mV, reflecting a prominent contribution of I_K-IR. Although resting membrane potential measured in rat cerebellar granule cells ranges from 60 to 85 mV (D'Angelo et al., 1995, 1998; Brickley et al., 1996; Watkins and Mathie, 1996; Rossi et al., 1998; Armano et al., 2000), at 80 mV the model was rather insensitive to manipulation of ionic conductances, which provided a useful reference potential for subsequent simulations.

	RESULTS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Bursting and resonance in granule cell

Intrinsic granule cell electroresponsiveness was investigated in current-clamp recordings. During step current injection (Fig. 1A), granule cells showed inward rectification in the hyperpolarizing direction. Just-threshold depolarizing currents generated spikes, which could be clustered in doublets-triplets or longer bursts occurring at a frequency of 3-10/sec (Fig. 1A,B) [also see D'Angelo et al. (1998), their Figs. 1, 2]. Spikes were followed by a fast afterhyperpolarization (AHP), an afterdepolarization, and a slow afterhyperpolarization (Fig. 1B). When stronger depolarizing currents were injected, firing became regular, and afterdepolarization and slow afterhyperpolarization were no longer observed (Fig. 1A, top trace). It should be noted that recordings were performed in the presence of 10 µM bicuculline preventing granule cell rhythmic inhibition by Golgi cells (Brickley et al., 1996) and that spontaneous EPSPs were too rare (~0.1/sec) to significantly affect spike generation (no difference was noted after application of the glutamate receptor blockers 10 µM CNQX, 100 µM APV, and 50 µM 7-Cl-kyn; n = 3; data not shown). Thus, intrinsic membrane mechanisms should generate spike bursts, slow afterhyperpolarization (Alonso and Llinas, 1989; Gutfreund et al., 1995; Pape and Driesang, 1998), and spike afterdepolarization (Azouz et al., 1996; Raman and Bean, 1997).

View larger version (29K): [in this window] [in a new window]   Figure 1.   Granule cell electroresponsiveness during step current injection. A, The injection of current steps (from 8 to 6 pA, resting potential = 62 mV) causes inward rectification in the hyperpolarizing direction. Spikes are activated around 40 mV. The tracing at 4 pA shows a single spike, the tracing at 6 pA shows spikes clustered in two bursts, and the tracing at 8 pA shows regular repetitive firing. B, Just-threshold response illustrating spike fast afterhyperpolarization (fAHP), slow afterhyperpolarization (sAHP), and afterdepolarization (ADP). The neuron in B is different from that in A (spikes are truncated). Recordings in this and the following figures were performed in the presence of 10 µM bicuculline.

Injection of sinusoidal currents of appropriate amplitude generated spike bursts in correspondence with the positive phase of the stimulus (Fig. 2A). Spike frequency within bursts increased and then decreased with injected current frequency, therefore showing resonance. The resonance frequency was 8.9 ± 3.2 Hz (n = 8; average data are shown in Fig. 10A). Resonance was also observed in the absence of spikes (1 µM TTX in the bath) (Fig. 2B). In this case, the maximum depolarization reached during the positive phase of the sinusoidal voltage response showed a resonance frequency of 8.1 ± 2.9 Hz (n = 8; average data are shown in Fig. 10B). Regarding both spike frequency and membrane potential measurements, the resonance frequency tended to increase slightly with the intensity of the injected current (Fig. 2A,B; see Fig. 10A,B).

View larger version (32K): [in this window] [in a new window]   Figure 2.   Resonance in a cerebellar granule cell (same cell as in Fig. 1A). A, Injection of sinusoidal currents at various frequencies (0.5-40 Hz) reveals resonance in burst spike frequency, which was measured by dividing the time period between the first and last spike in a burst by the number of interspike intervals. The plot shows that the resonance frequency was 6 Hz, with sinusoidal currents of ±6 pA () and 8 Hz with ±8 pA (). At frequencies higher than those shown in the plot, just one or no spikes were generated, and spike frequency fell to zero. B, After 1 µM TTX perfusion, injection of sinusoidal currents at various frequencies reveals resonance in the maximum depolarization reached during the positive phase of the sinusoidal voltage response. The plot shows that the resonance frequency was 6 Hz, with sinusoidal currents of ± 6 pA () and 8 Hz with ±8 pA (). Beyond the resonance peak, the sinusoidal voltage response decreased monotonically until 40 Hz (data not shown).

Evidence for a K⁺-dependent mechanism in oscillation and resonance

The results shown in Figures 1 and 2 suggest that granule cells combine membrane mechanisms that generate fast repetitive firing, with others responsible for slow oscillations and resonance in the theta-frequency range. Accordingly, previous observations unveiled oscillations that sustained spike bursting by reducing K⁺ conductances with TEA and showed that they depended on a persistent TTX-sensitive Na⁺ current (D'Angelo et al., 1998). Here we investigated the ionic dependence of the slow oscillatory mechanism and of resonance during bath application of 1 mM Ni²⁺, which fully blocks granule cell Ca²⁺ currents (Rossi et al., 1994; Tottene et al., 1996; D'Angelo et al., 1997, 1998), and of 4 mM 4-AP and 20 mM TEA, which block the granule cell K⁺ currents I_K-V, I_K-Ca, and I_K-A (Cull-Candy et al., 1989; Bardoni and Belluzzi, 1993).

When these ionic channel blockers were used, depolarizing current steps sustained large-size oscillations surmounted by a solitary spike when the patch pipette solution contained K⁺ (n = 5) but not when it contained Cs⁺ (n = 6) as the main intracellular cation (Fig. 3A). In these recordings, a marked adaptation prevented repetitive spike activation and bursting. Oscillations and solitary spikes were blocked by 1 µM TTX, unveiling their Na⁺ dependence. Likewise, resonance was observed when the patch pipette contained K⁺ (n = 3) but not when it contained Cs⁺ (n = 3) (Fig. 3B). Because Cs⁺ prevents K⁺ permeation through K⁺ channels, and because Ca²⁺ channels in these recordings are blocked, oscillation and resonance are likely to involve a TEA-insensitive Ca²⁺-independent K⁺ current.

View larger version (14K): [in this window] [in a new window]   Figure 3.   K+ dependence of slow oscillation and resonance. Current-clamp recordings were performed in the presence of 20 mM TEA, 4 mM 4-AP, and 1 mM Ni2+. A, A sustained slow oscillation is observed in granule cells recorded with K+-containing patch pipette during step current injection (10 pA from 80 mV). A solitary action potential is generated in a different cell recorded with a Cs+-containing patch pipette. In both cases, excitable responses were abolished by 1 µM TTX. B, Resonance curves in a cell recorded with K+ () and in another cell recorded with Cs+ () inside the patch pipette. Comparable voltage responses in neurons shown in A and B were obtained by properly adjusting the intensity of injected current (lower with Cs+- than with K+-containing pipettes).

Isolation of a slow K⁺ current in granule cells

A TEA-insensitive Ca²⁺-independent slow outward current (I_K-slow) was isolated by performing voltage-clamp recordings in the presence of 20 mM TEA, 4 mM 4-AP, 1 mM Ni²⁺, and 1 µM TTX (Fig. 4). I_K-slow could not be measured in cells internally perfused with Cs⁺ rather than K⁺ (data not shown), which revealed its K⁺ dependence, and was reversibly reduced by 48.7 ± 8.4% (n = 4) by extracellular perfusion of 1 mM Ba²⁺ (Fig. 4A).

View larger version (19K): [in this window] [in a new window]   Figure 4.   Isolation of a slow K+ current, IK-slow, in the presence of 20 mM TEA, 4 mM 4-AP, 1 mM Ni2+, and 1 µM TTX. A, IK-slow was generated by a voltage pulse from 80 to +30 mV. IK-slow was reversibly inhibited by application of 1 mM Ba2+. B, IK-slow activation was investigated by applying 1 sec, 10 mV depolarizing voltage steps from the holding potential of 80 mV (a short pre-step was applied to inactivate IK-A) (Bardoni and Belluzzi, 1993). The inset shows exponential fitting to the rising phase of a current recorded at 0 mV with the function I(t) = Iss * (1 - exp(t/act)), where Iss = 53 pA is the steady-state current and  act = 33.5 msec is the activation time constant. C, IK-slow deactivation was investigated by using voltage jumps to different potentials after a 300 msec conditioning pulse at +30 mV (holding potential = 80 mV). The inset shows exponential fitting to a tail current recorded at 10 mV with the function I(t) = Iss + I0 * exp(t/deact), where (Io + Iss) = 48.9 pA is the instantaneous current, Iss = 39.2 pA is the steady-state current, and deact = 46.9 msec is the deactivation time constant. D, Voltage dependence of steady-state amplitude of deactivation curves (Iss, ), and of time constants obtained by exponential fitting to activation (act, ) and deactivation (deact, ) curves. The inset shows intersection of the linear regression curve to instantaneous tail current amplitude with the voltage axis at 71.4 mV. Data in B-D were obtained from the same granule cell.

During application of depolarizing voltage steps from the holding potential of 80 mV (Fig. 4B), I_K-slow activated around 40 mV, and its amplitude increased by increasing the test potential. I_K-slow rising phase was well fitted by a single exponential function (Fig. 4B, inset), indicating first-order activation kinetics. Exponential time constants ranged from 10 to >100 msec for test potentials between 40 and 0 mV, reaching values two orders of magnitude higher than those of other granule cell outward currents (Cull-Candy et al., 1989; Bardoni and Belluzzi, 1993). I_K-slow persisted for >1 sec, but in six of nine cells it showed a slight tendency to inactivate at positive test potentials (<15% after 1 sec at +30 mV).

Voltage jumps to different potentials after a 300 msec conditioning pulse at +30 mV (when I_K-slow was almost fully activated) generated tail currents that relaxed with an exponential time course (Fig. 4C). Exponential fitting (Fig. 4B, inset) yielded the activation time constant, the instantaneous tail current amplitude (time = 0), and the steady-state current. According to the Nernstian equilibrium potential of K⁺ ions, the instantaneous tail current was zero at 74.8 ± 5.5 mV (n = 4) and shifted to 37.7 ± 6.1 mV (n = 4) when extracellular K⁺ was raised to 25 mM by equimolar Na⁺ substitution.

Voltage dependence of the steady-state current and the activation time constant obtained from data shown in Figure 4, A and B, are shown in Figure 4D. The I-V curve showed outward rectification, and the activation time constant showed a bell-shaped voltage dependence. It should be noted that time constants measured from activation and deactivation currents coincided, consistent with a first-order gating mechanism. However, deactivation currents were preferred to activation currents for gating reconstruction, because they covered a more extended membrane potential range and were free of potential contamination by I_K-A, which is not completely blocked even by 4 mM 4-AP in cerebellar granule cells (Bardoni and Belluzzi, 1993).

Figure 5 shows average data obtained from nine granule cell recordings. The I-V curve (Fig. 5A) was fitted with a Boltzman equation of the form [see Rossi et al. (1998)]:

(8)

where G_max = 0.8 nS is the maximum conductance, V_rev = 70.2 mV is the reversal potential, V_1/2 = 20 mV is the half-activation potential, and k = 6 mV¹ is the activation voltage dependence. These parameters were used to reconstruct steady-state I_K-slow activation (Fig. 5A) by using the equation:

(9)

The relationship between activation time constant and voltage (Fig. 5B) was fitted by Equation 4 using the kinetic constants shown in Figure 5C. As in other ionic current models (Gutfreund et al., 1995; Mainen et al., 1995), independent expressions for the activation time constant and steady-state activation improved data representation. These results allowed I_K-slow reconstruction according to a first-order Hodgkin-Huxley kinetic scheme (Yamada et al., 1998) (see Materials and Methods; Table 1).

View larger version (14K): [in this window] [in a new window]   Figure 5.   Gating properties of IK-slow (average data from 9 granule cells, mean ± SEM). A, Average I-V relationship () fitted with Equation 8 (solid line). The broken line is the normalized steady-state activation curve (x (K-slow)) obtained with Equation 9. B, Average activation time constant ((K-slow)) versus membrane potential (). The fitting line was obtained from Equation 4 and the kinetic functions shown in C. C, Voltage dependence of the kinetic constants  and  (see Eq. 3 and Table 1).

View this table: [in this window] [in a new window]   Table 1.   Voltage-dependent conductance parameters

Mathematical reconstruction of intrinsic excitability

To investigate the role of I_K-slow in oscillations, bursting, and resonance, we developed a mathematical model of granule cell excitability (see Materials and Methods; Table 1). With I_Na-f, I_Ca, I_K-V, I_K-Ca, I_K-A, and I_K-IR, the model generated inward rectification and fast repetitive firing (Gabbiani et al., 1994). The model was then endowed with a persistent Na⁺ current (I_Na-p) (D'Angelo et al., 1998) and I_K-slow, and was extended to test the potential contribution of a resurgent Na⁺ current (I_Na-r) (Raman and Bean, 1997). With I_Na-p, I_K-slow, and I_Na-r, just-threshold step current injection generated membrane potential oscillations and spike bursts, with the spikes showing fast- and slow-afterhyper-polarization and afterdepolarization (Fig. 6A,B). With stronger current injection, the model generated regular repetitive firing (Fig. 6A; see Fig. 8). Moreover, the model reproduced TEA-induced bursting (see Fig. 9) and intrinsic resonance (see Fig. 10).

View larger version (15K): [in this window] [in a new window]   Figure 6.   Mathematical modeling of granule cell excitability. A, Model responses to 2 pA step current injection from 80 mV. The model generates inward rectification in subthreshold responses, followed by regular repetitive firing with almost no adaptation. B, Slow oscillations, slow afterhyperpolarization, occasional uncoupling of spike prepotential from upstroke, and spike bursts can be generated by the model by using just-threshold stimulation (10.5 pA in the top and middle tracings, 12 pA in the bottom tracing). The shape of oscillations and bursting could be modified by changing the GNa-r or GK-slow intensities.

The ionic mechanism of slow oscillations

Model simulations showed that I_Na-p and I_K-slow were sufficient to generate regular theta-frequency oscillations, which were eliminated when either of these currents was turned off (Fig. 7A). During oscillations, I_Na-p and I_K-slow showed activation/deactivation cycles describing stable orbits in the phase plane (Fig. 7B). The histerisis observed in I_K-slow (and to a lesser extent in I_Na-p) trajectory reflected its delayed gating during membrane potential changes.

View larger version (16K): [in this window] [in a new window]   Figure 7.   The oscillatory mechanism in granule cells. A, Simulation of stable oscillations sustained by GK-slow and GNa-p during injection of an 11 pA current step. Oscillations are eliminated by turning off either GK-slow or GNa-p. All other active conductances were set to zero except GK-IR, which was used to keep input conductance close to its normal value. B, Time course and phase-plane trajectory of IK-slow and INa-p during membrane potential oscillations (dotted line). C, Voltage responses to an 11 pA current step simulating TEA (GK-V = 0, GK-Ca = 0), and Ni2+ (GCa = 0) application. The broken line simulates subsequent application of TTX (GNa = 0). D, Voltage responses to an 11 pA current step simulating TEA (GK-V = 0, GK-Ca = 0) and TTX (GNa = 0) application. Same calibration in A, C, and D.

The model generated repetitive oscillations when G_K-V, G_K-Ca, and G_Ca were set to zero, reproducing the pharmacological application of 20 mM TEA and 1 mM Ni²⁺ (Fig. 7C; compare Fig. 2). Moreover, the model generated a solitary Ca²⁺ action potential when G_K-V, G_K-Ca, and G_Na were set to zero, reproducing the pharmacological application of 20 mM TEA and 1 µM TTX (Fig. 7D) [compare Fig. 3 in D'Angelo et al. (1998) and Fig. 10 in D'Angelo et al. (1997)]. Modeling results therefore indicated that slow oscillations in granule cells depend on the interaction of I_Na-p and I_K-slow but are independent from Ca²⁺ currents.

Repetitive firing

Figure 8A shows repetitive firing in the model. As in experimentally recorded granule cell responses (D'Angelo et al., 1995, 1998; Brickley et al., 1996), spikes showed negligible adaptation, and first-spike latency decreased whereas spike frequency increased when the injected current intensity was raised. The frequency-intensity (f-I) plot was almost linear between 0 and 100 Hz, with a slope of 7.3 spikes per pA¹ · sec¹.

View larger version (31K): [in this window] [in a new window]   Figure 8.   Ionic mechanisms of repetitive firing. A, Model tracings show regular firing with negligible adaptation at >100 Hz. The top plot (f-I plot) reports firing frequency in control conditions () and after having turned off IK-A (broken line) or INa-r (dotted line). The bottom plot shows first-spike latency in control conditions () and after having turned off IK-A (broken line). B, Ionic currents and [Ca2+] changes during repetitive firing. The right set of tracings is an enlargement of currents associated with an action potential. Note that both INa-p and INa-r are activated after the spike.

Spikes were generated by a sudden I_Na-f raise followed by activation of I_K-V and the I_Ca/I_K-Ca system (Fig. 8B). The intracellular Ca²⁺ wave peaked after the spike and decayed to zero in ~3 msec. I_Na-f, I_K-V, I_Ca, and I_K-Ca accounted for most of the active current during the spike and the fast AHP, whereas other currents were two orders of magnitude smaller. I_Na-p increased during the spike prepotential and persisted for several milliseconds. I_Na-r increased just after the spike. I_K-A, I_K-IR, and I_K-slow were mostly active during the interspike trajectory. Thus, whereas I_Na-p and I_Na-r enhanced spike activation, I_K-A, I_K-IR, and I_K-slow delayed it. The regulatory action of I_Na-r and I_K-A on spike frequency is shown in the plots of Figure 8A.

TEA-induced bursting

Figure 9A illustrates theta-frequency bursting in the model after partial I_K-Ca blockage, simulating experimental TEA application (D'Angelo et al., 1998, their Figs. 6 and 8). A relatively modest I_K-Ca blockage caused spike doublets-triplets, whereas stronger blockage caused marked membrane potential oscillations surmounted by adapting spike bursts. Thus, I_K-Ca prevented activation of the I_Na-p/I_K-slow oscillatory mechanism. Bursting was enhanced by I_Na-r (Fig. 9A, compare left with right column), but no bursting was generated by I_Na-r alone when the I_Na-p/I_K-slow mechanism was turned off (Fig. 9A, bottom panel).

View larger version (34K): [in this window] [in a new window]   Figure 9.   Ionic mechanisms of bursting. A, Tracings show intensification of bursting and spike adaptation by progressively increasing GK-Ca inhibition during injection of 11 pA current steps. Stronger GK-Ca inhibition is needed to generate bursting when GNa-r is set to zero. No bursting is generated when GNa-p and GK-slow are turned off, but GNa-r is left active. B, Ionic currents and [Ca2+] changes during bursting elicited with IK-Ca reduction to 37% of its normal value. Note that IK-slow and INa-p are greatly enhanced during bursting compared with repetitive firing.

TEA-induced bursting was characterized by a remarkable rise in I_Na-p associated with a progressive I_K-slow activation (Fig. 9B, compare Fig. 7B). During the burst, I_Na-f inactivation caused spike amplitude adaptation.

Resonance

Injecting the granule cell model with sinusoidal currents of different frequency generated resonant responses. Resonance in burst spike frequency is shown in Figure 10A, and resonance in maximum membrane depolarization (with spikes blocked by setting I_Na-f, I_Na-p, and I_Na-r to zero mimicking TTX block) is shown in Figure 10B. In both cases a family of resonant curves is shown, one of which matches average experimental measurements from a set of five granule cells.

View larger version (27K): [in this window] [in a new window]   Figure 10.   Ionic mechanisms of resonance. A, Injection of sinusoidal currents causes oscillatory bursting in the model. Tracings are generated by a ±6 pA sinusoidal current superimposed on a 12 pA current step. Insets show higher spike frequency in bursts generated at 10 Hz than at 2 Hz. The plot shows model resonance with three different sinusoidal current intensities (±4, ±6, or ± 8 pA superimposed on a constant 12 pA current step). The curve generated with ±6 pA is a good match with the average experimental response (; mean ± SD; n = 5). B, Same as in A, except that maximum membrane depolarization during the positive phase of sinusoidal voltage responses is measured with INa-f, INa-p, and INa-r set to zero. This result is compared with experimental recordings in the presence of 1 µM TTX (; mean ± SD; n = 5, same cells as in A). As with real granule cells, the model shows resonance ~10 Hz.

The model was used to investigate the ionic mechanisms of resonance. Resonance was eliminated by blocking I_K-slow, which determined the ascending branch of resonance curves (Fig. 11A,B). The descending branch of resonance curves was determined by passive membrane filtering, as demonstrated by its persistence when the ionic channels involved in resonance were blocked. In addition, I_K-A markedly accelerated the descending branch. Enhanced activation of I_K-slow at low frequency and I_K-A at high frequency is shown in Figure 11C. It should be noted that resonance was not eliminated by blocking I_Na-p (Fig. 11A,B), which showed minor frequency-dependent changes (Fig. 11C). Finally, I_Na-r (in the case of burst spike frequency) decelerated the descending branch of resonance curves.

View larger version (29K): [in this window] [in a new window]   Figure 11.   Resonance regulation. This Figure shows resonance being regulated by injection of a ±6 pA sinusoidal current superimposed on a 12 pA current step. A, Resonance in burst spike frequency in different conditions: , control; , IK-slow = 0; , GK-A = 0; , GNa-p = 0; , GNa-r = 0. No resonance could be observed in the model when INa-p, IK-slow, and INa-r were turned off (thin dotted line). B, Resonance in maximum membrane depolarization during the positive phase of sinusoidal voltage responses in different conditions: control (, INa-p, INa-f, INa-r = 0; TTX condition); , IK-slow = 0, , IK-A = 0; , INa-p active. No resonance could be observed in the model when INa-p, IK-slow, and INa-r were turned off (thin dotted line). C, INa-p, IK-slow, and IK-A at three different frequencies (thick line, 8 Hz; thin line, 2 Hz; broken line, 14 Hz). Note the different frequency-dependent activation of these currents.

	DISCUSSION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

This study shows that in addition to generating fast repetitive firing, cerebellar granule cells generate oscillations, bursting, and resonance in the theta-frequency range. Experimental and modeling results indicated that these aspects of intrinsic excitability require a slow K⁺ current (I_K-slow) to be generated.

I_K-slow is a Ca²⁺-independent TEA-insensitive K⁺ current activating in the spike threshold region with slow kinetics (10-100 msec). A current like I_K-slow has not been reported previously in cerebellar granule cells, although a persistent TEA-resistant current component was noted in cell culture (Cull-Candy et al., 1989) [also see Bardoni and Belluzzi (1993), their Fig. 10B]. I_K-slow biophysical and pharmacological properties are similar to those of I_M of vertebrate neurons (Brown and Adams, 1980; Adams et al., 1982a,b) and are suitable to generate the delayed repolarizing feedback and high-pass filtering required for bursting and resonance. Indeed, I_M has been reported to sustain oscillations and resonance in amygdaloid neurons (Pape and Driesang, 1998) and in cortical pyramidal neurons (Gutfreund et al., 1995). Alternative mechanisms that might be invoked to explain bursting and resonance are unlikely to occur in cerebellar granule cells. (1) A slow inward rectifier current (I_h) (Dickson et al., 2000) and a slow Ca²⁺-dependent K⁺ current (I_AHP) (Wang and Rinzel, 1999) are apparently not expressed in granule cells. Granule cell inward rectification is fully explained by a fast K⁺-dependent inward rectifier (Fig. 6A) (Rossi et al., 1998), and granule cell slow afterhyperpolarization is fully explained by I_K-slow (Fig. 6B). It should also be noted that the I_AHP blocker apamin did not affect the excitable response, and oscillations and bursting persisted in the presence of Ca²⁺ channel blockers (Fig. 3) (D'Angelo et al., 1998). Moreover, contrary to what would be expected from I_AHP, granule cell slow AHP occurred with just-threshold stimulation disappearing rather than being enhanced with higher spike frequency and was not associated with any spike frequency adaptation. (2) A resurgent current, I_Na-r, facilitated but proved not sufficient to induce bursting and resonance (Fig. 8A). (3) Finally, return currents from the dendrites, which might generate somatic rebound depolarization and spike bursting, are unlikely to be effective because of the granule cell compact electrotonic structure. Actually, with a 10 M somatodendritic resistance (calculated assuming an axial resistance of 150 /cm in a 20 µm dendrite) and a somatodendritic ratio between 1 and 2 (Silver et al., 1992; D'Angelo et al., 1993; Gabbiani et al., 1994), no bursting is expected even if active conductances are expressed in the dendrites (Mainen and Sejinowsky, 1996).

Modeling reliability was based on the extensive characterization of membrane currents and the compact electrotonic structure of cerebellar granule cells (for details, see Materials and Methods). Indeed, patch-clamp recordings from mature cerebellar granule cells in situ have been used (1) to reconstruct I_Ca (Rossi et al., 1994), I_K-IR (Rossi et al., 1998), I_K-A (Bardoni and Belluzzi, 1993), and I_K-slow (this study), (2) to identify different Na⁺ current components (including I_Na-f, I_Na-p, and I_Na-r) (D'Angelo and Magistretti, unpublished results), and (3) to characterize the pharmacological properties of excitability (D'Angelo et al., 1998). Patch-clamp recordings from granule cells in culture were used to reconstruct I_K-Ca (Fagni et al., 1991; Gabbiani et al., 1994). Although further investigation is required to clarify Ca²⁺ dynamics and the biophysical properties of the Na⁺ current, mathematical modeling allowed us to investigate the role of I_Kslow in relationship to other excitable properties of the granule cell. In addition to modeling, I_Kslow involvement in oscillations, bursting, and resonance may be investigated by I_Kslow selective pharmacological blockage or by I_Kslow electronic antagonism/expression through a dynamic-clamp circuit (Hutcheon et al., 1996a,b).

Experimental and modeling observations allow the following reconstruction of the ionic mechanisms of rat cerebellar granule cell electroresponsiveness.

(1) I_Na-f and I_K-V were the core of a fast oscillatory mechanism sustaining fast repetitive firing, as in the classical Hodgkin-Huxley model (Hodgkin and Huxley, 1952). I_K-A increased spike latency (Connor and Stevens, 1971), and the I_Ca-I_K-Ca system stabilized repetitive firing by enhancing fast AHP and Na⁺ channel deinactivation [see also Gabbiani et al. (1994)].

(2) I_Na-p and I_K-slow were the core of a slow oscillatory mechanism sustaining theta-frequency oscillation and bursting. I_Na-p has been reported to sustain oscillations in association with I_M in neocortical pyramidal neurons (Gutfreund et al., 1995) and in association with I_h in enthorinal pyramidal neurons (Alonso and Llinas, 1989; Dickson et al., 2000). I_K-slow caused delayed repolarization, terminating the positive phase of the oscillation promoted by I_Na-p.

(3) Emergence of bursting was regulated by the I_Ca-I_K-Ca system (D'Angelo et al., 1998) through the fast AHP, which reduced the depolarizing action of I_Na-p (Azouz et al., 1996) and I_Na-r (Raman and Bean, 1997) after the spike.

(4) Resonance depended on high-pass filtering caused by I_K-slow (generating the ascending branch of the resonance curve) in association with low-pass filtering caused by passive membrane properties (generating the descending branch of the resonance curve) and was amplified by I_Na-p. Thus, resonance involved subthreshold changes in membrane excitability [corresponding to subthreshold impedance resonance reported by Hutcheon et al. (1996a,b)], consistent with the identification of I_K-slow as a "resonator" current and I_Na-p as an "amplifier" current" (Hutcheon and Yarom, 2000). In addition, the model suggested that I_K-A, which is inactivated at low but not high frequency, accelerated the descending branch of the resonance curve. Finally, the model predicted that resonance in burst spike frequency could be intensified by spike clustering promoted by I_Na-r. Thus, resonance and bursting may have in common their dependence on I_K-slow as well as that on I_Na-r.

Granule cells are excitatory interneurons relaying information conveyed by mossy fibers into the cerebellar cortex (Marr, 1969). Although repetitive firing implements a mechanism of linear frequency coding (Gabbiani et al., 1994), intrinsic bursting and resonance tune granule cells on theta-frequency, which is diffused in sensory-motor structures (Koch, 1999). Intrinsic bursting and resonance may favor granule cell phase-locking through the recurrent Golgi cell inhibitory circuit (Maex and DeShutter, 1998). Theta-frequency discharge has indeed been recorded from the granular layer during specific activity states in vivo (Pellerin and Lamarre, 1997; Hartmann and Bower, 1998). Finally, as observed in the hippocampus (Holsher et al., 1997), theta-frequency bursting may regulate the induction of synaptic plasticity in the mossy fiber-granule cell pathway (D'Angelo et al., 1999; Armano et al., 2000).

	FOOTNOTES

Received Aug. 10, 2000; revised Oct. 30, 2000; accepted Nov. 3, 2000.

This work was supported by European Community Grants PL97 0182 and PL97 6060, and by Instituto Nazionale per la Fisica della Materia. We acknowledge Marja-Leena Linne and Massimiliano Zaniboni for contributing to preliminary simulations, and Lia Forti and Elisabetta Sola for their helpful comments on this manuscript.

Correspondence should be addressed to Egidio D'Angelo, Department of Cellular/Molecular Physiology and Pharmacology, Via Forlanini 6, I-27100 Pavia, Italy. E-mail: dangelo{at}unipv.it.

	REFERENCES

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

• Adams PR, Brown DA, Constanti A (1982a) M-currents and other potassium currents in bullfrog sympathetic neurones. J Physiol (Lond) 330:537-572[Abstract/Free Full Text].
• Adams PR, Brown DA, Constanti A (1982b) Pharmacological inhibition of the M-current. J Physiol (Lond) 332:223-262[Abstract/Free Full Text].
• Alonso A, Llinas RR (1989) Subthreshold Na⁺-dependent theta-like rhythmicity in stellate cells of entorhinal cortex layer II. Nature 342:175-177[Medline].
• Armano S, Rossi P, Taglietti V, D'Angelo E (2000) Long-term potentiation of intrinsic excitability at the mossy fiber-granule cell synapse of rat cerebellum. J Neurosci 20:5208-5216[Abstract/Free Full Text].
• Azouz R, Jensen MS, Yaari Y (1996) Ionic basis of spike after-depolarization and burst generation in adult rat hippocampal CA1 pyramidal cells. J Physiol (Lond) 492:211-223[ISI][Medline].
• Bardoni R, Belluzzi O (1993) Kinetic study and numerical reconstruction of A-type current in granule cells of rat cerebellar slices. J Neurophysiol 69:2222-2231[Abstract/Free Full Text].
• Brickley SG, Cull-Candy SG, Farrant M (1996) Development of a tonic form of synaptic inhibition in rat cerebellar granule cells resulting from persistent activation of GABA_A receptors. J Physiol (Lond) 497:753-759[ISI][Medline].
• Brown DA, Adams PR (1980) Muscarinic suppression of a novel voltage-sensitive K⁺ current in a vertebrate neurone. Nature 283:673-676[Medline].
• Connor JA, Stevens CF (1971) Predictions of repetitive firing behaviour from voltage-clamp data on an isolated neuron soma. J Physiol (Lond) 213:31-53[Abstract/Free Full Text].
• Connors BW, Gutnick MJ (1990) Intrinsic firing patterns of diverse neocortical neurons. Trends Neurosci 13:99-104[ISI][Medline].
• Cull-Candy SG, Marshall CG, Ogden D (1989) Voltage-activated membrane currents in rat cerebellar granule neurones. J Physiol (Lond) 414:179-199[Abstract/Free Full Text].
• D'Angelo E, Rossi P, Taglietti V (1993) Different proportions of N-methyl-D-aspartate and non-N-methyl-D-aspartate receptor currents at the mossy fibre-granule cell synapse of developing rat cerebellum. Neuroscience 53:121-130[ISI][Medline].
• D'Angelo E, De Filippi G, Rossi P, Taglietti V (1995) Synaptic excitation of individual rat cerebellar granule cells in situ: evidence for the role of NMDA receptors. J Physiol (Lond) 484:397-413[ISI]