| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
1 Sloan-Swartz Center for Theoretical Neurobiology, Salk Institute, La Jolla, California 92037 2 Howard Hughes Medical Institute, and Computational Neurobiology Lab, Salk Institute, La Jolla, California 92037 3 Department of Physics and Astronomy, University of North Carolina, Chapel Hill, North Carolina 27599 4 Department of Biology, University of California San Diego, La Jolla, California 92037 5 Institute for Theoretical Biology, Humboldt-University Berlin, D-10115 Berlin, Germany
Submitted 9 June 2003; accepted in final form 15 September 2003
 |
ABSTRACT |
|---|
|
|
|
INTRODUCTION |
|---|
|
; Magee 2002; Mainen and Sejnowski 1996; Marder et al. 1996; Turrigiano et al. 1994). In this study we explore the influence of ionic conductances on the reliability of the timing of spikes of cortical cells. Robustness of spike timing to physiological noise is the prerequisite for a spike timing–based code, and has recently been investigated (Beierholm et al. 2001; Brette and Guigon 2003; Fellous et al. 2001; Fricker and Miles 2000; Gutkin et al. 2003; Mainen and Sejnowski 1995; Reinagel and Reid 2002; Tiesinga et al. 2002). It has been found experimentally that different types of neurons are tuned to different stimuli with respect to spike timing reliability. For example, cortical interneurons show maximum reliability in response to higher-frequency sinusoidal stimuli, whereas pyramidal cells respond more reliably to lower-frequency sinusoidal inputs (Fellous et al. 2001). An important difference between those types of neurons is the composition of their ion channels. Taking into account that effective numbers of ion channels can be adjusted on short time scales through neuromodulation, changes in ion channels may also provide a useful way for a neuron to dynamically maximize spike timing reliability according to the properties of the input.
Spike timing reliability is enhanced with increasing stimulus amplitude (Mainen and Sejnowski 1995). In the intermediate amplitude regime, the frequency content of the stimulus is an important factor determining reliability (Fellous et al. 2001; Haas and White 2002; Hunter and Milton 2003; Hunter et al. 1998; Jensen 1998; Nowak et al. 1997; Tiesinga 2002). Spike timing reliability of a neuron is maximal for those stimuli that contain frequencies matching the intrinsic frequency of a neuron (Hunter et al. 1998). The intrinsic (or preferred) frequency is given by the firing rate of a neuron in response to the DC component of the stimulus. Because of the relation to the DC firing rate of a neuron, both the DC value (whether the stimulus mean or additional synaptic input) and the conductances of a cell can be expected to influence the spike timing frequency preference. The former was recently shown by Hunter and Milton (2003). The influence of conductances (rather than injected current) on spike timing reliability through changes in the neuronal activity according to the resonance effect is the focus of the first part of this paper, see RESULTS (Influence of conductances on the frequency preference).
In this part we specifically seek to understand which ionic conductances of cortical neurons can mediate changes of the preferred frequency (with respect to spike timing reliability) over a broad range of frequencies. Reliability is assessed on the basis of the robustness of spike timing to noise (of amplitude smaller than the stimulus amplitude). Injecting sinusoidal currents on top of a DC current into conductance-based model neurons, we confirm that spike timing reliability is frequency-dependent as predicted by the resonance effect. We show that reliability can be regulated at the level of ion channel populations, and identify the slow potassium channels as powerful to influence the preferred frequency. Our simulations support that the influence of ion channels on spike timing reliability also holds for more realistic rhythmic stimulus waveforms. Dynamic-clamp experiments in slices of rat prefrontal cortex confirm the theoretical prediction that slow potassium channels can mediate a change in spike timing reliability, dependent on the frequency of the input.
In the second part of the RESULTS section (Influence of conductances on spike timing reliability at the preferred frequency) we explore the influence of ion channels beyond changes in preferred frequency attributed to the resonance effect. Different neurons may have the same preferred frequency (i.e., the same DC firing rate) but different composition of ion channels. We analyze the influence of slow potassium channels and persistent sodium channels on spike timing reliability for neurons with the same preferred frequency. We find that both channel types significantly influence spike timing reliability. Slow potassium channels increase reliability, whereas persistent sodium channels lower it.
|
|
METHODS |
|---|
|
The single-compartment conductance-based model neurons were implemented in NEURON. In the basic implementation, the neurons contained fast sodium channels (Na), delayed-rectifier potassium channels (Kdr), leak channels (leak), slow potassium channels (Ks), and persistent sodium channels (NaP). The time resolution of the numerical simulation was 0.1 ms. The kinetic parameters of the 5 basic channel types and reversal potentials were taken from a model of a cortical pyramidal cell Golomb and Amitai (1997), apart from the reversal potential of the leak channels, which was set to –80 mV (to avoid spikes in the absence of input and noise). The conductances of the cell we will refer to as the reference cell were (in mS/cm2): gNa = 24, gKdr = 3, gleak = 0.02, gKs = 1, and gNaP = 0.07 (Golomb and Amitai 1997). Its input resistance was 186 M
. The slow potassium conductance represented potassium channels with an activation time on the order of several tens to hundreds of milliseconds (here 75 ms). In the model it is responsible for a spike frequency adaptation to a current step, which is experimentally observed in cortical pyramidal neurons (Connors and Gutnick 1990; McCormick et al. 1985). We also investigated cells where the Ks channels were replaced by muscarinic potassium channels (KM) and by calcium-dependent potassium channels (KCa). The muscarinic channel KM was a slow noninactivating potassium channel with Hodgkin–Huxley style kinetics (Barkai et al. 1994; Storm 1990). The calcium-dependent conductance KCa was based on first-order kinetics and was responsible for a slow afterhyperpolarization (Tanabe et al. 1998). This channel was activated by intracellular calcium and did not depend on voltage. Because of the dependency of KCa on calcium, we also inserted an L-type calcium channel as well as a simple Ca-ATPase pump and internal buffering of calcium. For the parameters of these additional currents see APPENDIX. Kinetic parameters of all channels used were set to 36°C.
Stimulus waveforms
The stimuli used to characterize spike timing reliability of individual cells consisted of 2 components. The first component was a constant depolarizing current IDC, which was the same for all model cells (apart from the simulations designed to study of the influence of the DC), and which also remained fixed throughout experimental recording of a cell. The second component was a sine wave with frequency f
 |
|
= 20 pA] on each individual trial. For the reference cell the noise resulted in voltage fluctuations of about 1.3 mV SD at rest. The noise was generated from a Gaussian distribution and filtered with an alpha function with a time constant of 
= 3 ms. Although overall reliability systematically decreased with the size of the noise, neither the frequency content of the noise nor the absolute size of the noise (in that range) significantly changed the results. Spike times were determined as the time when the voltage crossed –20 mV from below. The input resistance for the model cells was estimated by application of a depolarizing DC current step sufficiently large to depolarize the cell by 
10 mV.
Model neurons were also tested with a stimulus where power was distributed around one dominant frequency. These more realistic stimuli were constructed to have a peak in the power spectrum in either the theta- or the gamma-frequency range. These waveforms mimic theta- and gamma-type inputs and were created by inverse Fourier transform of the power spectrum (with random phases). For the theta-rich wave, the power spectrum consisted of a large peak at 8 Hz (Gaussian, = 1 Hz) and a small peak at 50 Hz (Gaussian, = 6 Hz). For gamma-rich waves, the power spectrum had a large peak at 30 or 50 Hz (with = 3 and 6 Hz, respectively), and a small peak at 8 Hz ( = 1 Hz). These waveforms were first normalized to have a root-mean-square (RMS) value of 1 and were then used with different scaling factors (yielding different RMS values). The DC component was added after scaling. These stimuli were presented for 10 s and when evaluating reliability, the first 500 ms after stimulus onset were discarded.
The reliability measure
Spike timing reliability was calculated from the neuronal responses to repeated presentations of the same stimulus. For the model studies this implied the same initial conditions, but different noise for each trial. Reliability was quantified by a correlation-based measure, which relies on the structure of individual trials and does not require the definition of a priori events. For a more detailed discussion of the method see Schreiber et al. (2003). The spike trains obtained from N repeated presentations of the same stimulus were smoothed with a Gaussian filter of width 2t, and then pairwise correlated. The normalized value of the correlation was averaged over all pairs. The correlation measure Rcorr, based on the smoothed spike trains, 
), is
 |
[0; 1]. Rcorr = 1 indicates the highest reliability and Rcorr = 0 the lowest. For all model cell studies, t = 1.8 ms and for the experimental data t = 3 ms. The value of t for model cells was chosen such that, given the noise level, reliability values Rcorr exploited the possible range of its values [0; 1], allowing for better discrimination between reliable and unreliable spike timing. The experimental data proved more noisy and therefore a larger t was chosen to yield a good distinction between reliable and unreliable states. All evaluation of model and experimental data (beyond obtaining spike times) was performed in Matlab. Firing rate analysis
For the firing rate analysis, the full parameter space of Na, NaP, Kdr, Ks, and leak conductances was analyzed. DC firing rates were obtained for all possible parameter combinations within the parameter space of the 5 conductances considered (see APPENDIX). The maximum change in firing rate achievable by one ion channel type was characterized (for each combination of the other 4 conductances) as the difference between the maximum and minimum (nonzero) firing rates achievable by variation of the ion channel conductance of interest, keeping the other 4 conductances fixed. If a cell never fired despite variation in one conductance, it was excluded from the parameter space (<5% of the total 4-dimensional conductance space for any channel type tested). The distribution of maximum changes in firing rate achievable by variation of the density of one ion channel type over all combinations of the other 4 densities is presented in the paper.
Experimental protocols
Coronal slices of rat prelimbic and infra limbic areas of prefrontal cortex were obtained from 2- to 4-wk-old Sprague-Dawley rats. Rats were anesthetized with Isoflurane (Abbott Laboratories, North Chicago, IL) and decapitated. Brains were removed and cut into 350-µm-thick slices using standard techniques. Patch-clamp was performed under visual control at 30–32°C. In most experiments Lucifer yellow (RBI, 0.4%) or Biocytin (Sigma, 0.5%) was added to the internal solution. In all experiments, synaptic transmission was blocked by D-2-amino-5-phosphonovaleric acid (D-APV; 50 µM), 6,7-dinitroquinoxaline-2,3, dione (DNQX; 10 µM), and biccuculine methiodide (Bicc; 20 µM). All drugs were obtained from RBI or Sigma, freshly prepared in artificial cerebrospinal fluid, and bath applied. Whole cell patch-clamp recordings were achieved using glass electrodes (4–10 M) containing (in mM): KMeSO4, 140; Hepes, 10; NaCl, 4; EGTA, 0.1; MgATP, 4; MgGTP, 0.3; phosphocreatine, 14. Data were acquired in current-clamp mode using an Axoclamp 2A amplifier (Axon Instruments, Foster City, CA).
Data were acquired using 2 computers. The first computer was used for standard data acquisition and current injection. Programs were written using Labview 6.1 (National Instrument, Austin, TX) and data were acquired with a PCI16E1 data acquisition board (National Instrument). Data acquisition rate was either 10 or 20 kHz. The second computer was dedicated to dynamic clamp. Programs were written using either a Labview RT 5.1 (National Instrument) or a Dapview (Microstar Laboratory, Bellevue, WA) frontend and a C language backend. Dynamic clamp (Hughes et al. 1998; Jaeger and Bower 1999; Sharp et al. 1993) was implemented using a DAP5216a board (Microstar Laboratory) at a rate of 10 kHz. A dynamic clamp was achieved by implementing a rapid (0.1-ms) acquisition/injection loop in current-clamp mode. All experiments were carried in accordance with animal protocols approved by the N.I.H.
Stimuli consisted of sine waves of 30 different frequencies (1–30 Hz) presented for 2 s. Only one amplitude was tested. No additional noise was injected. The first 500 ms were discarded for analysis of reliability.
|
|
RESULTS |
|---|
|
Model cells were stimulated with a set of sine waves on top of a fixed DC. Reliability values for each individual stimulus and cell, based on correlation of responses to repeated presentation of a stimulus each with an independent realization of the noise, were derived as a function of the frequency f and the amplitude of the sine component C. Figure 1, B–E show examples of 2 stimuli used and responses to those stimuli obtained from the reference cell. Figure 1F shows the complete set of reliability values as a function of frequency and amplitude of the sine component of the input.
Distinct, tongue-shaped regions of high reliability, so-called Arnold tongues (Beierholm et al. 2001), arising from the resonance effect of spike timing reliability, are visible. Figure 1F also shows that the degree of reliability depended on the power of the input at the resonant frequency of a neuron. The higher the amplitude at the resonant frequency, the more pronounced was the reliability. At high amplitudes, frequencies close to the resonant frequency also showed enhanced reliability. The Arnold tongues were approximately vertical, so that the frequency of maximum reliability showed only a weak dependency on the amplitude of the sine component. The difference in input frequency for maximal reliability, as the amplitude C varied from 0.05 to 0.15 nA, was usually <2 Hz. In most examples presented in this study, the strongest resonance was found at a 1:1 locking to the stimulus, where one spike per cycle of the sine wave was elicited. Additional regions of enhanced reliability could be observed at harmonics of the main resonant frequency (1:2, 1:3, and 1:4 phase locking, in order of decreasing strength), and at the 1st subharmonic (2:1 phase locking).
The location of the strongest Arnold tongue in frequency space revealed the preferred frequency of a neuron, which was well approximated by the firing rate of the neuron in response to the DC component alone. Figure 1G shows a strong correlation between the preferred frequency (i.e., position of the strongest Arnold tongue on the frequency axis determined by the frequency of highest reliability for a given amplitude, C) and the DC firing rate of a cell for a wide range of conductance values in the model (see APPENDIX). In all (but 2) cases the resonant frequency was close to the DC firing rate. Usually, the resonant frequency at the lowest amplitude of the sine component was closest to the DC firing frequency. For the 2 outliers the highest value of reliability was achieved at the subharmonic, or the 1st harmonic of the DC firing frequency. The importance of the DC firing rate in generating phase-locked firing patterns was previously emphasized (see e.g. Coombes and Bressloff 1999; Hunter et al. 1998; Keener et al. 1981; Knight 1972; Rescigno et al. 1970). The resonant frequency is referred to as preferred frequency throughout the paper.
Influence of conductances on the frequency preference
Because ionic conductances are known to influence neuronal activity levels, we investigated the ability of ion channels to modulate the preferred frequency in the first part of this study.
SIMULATION RESULTS FOR A CORTICAL SINGLE-COMPARTMENT MODEL CELL. We started from the model of a cortical neuron (the reference cell). First, we varied one channel density at a time, keeping the densities of the other channels fixed. The Arnold plots of cells whose leak density and slow potassium density were varied respectively are shown in Fig. 2. Example spike shapes (at DC stimulation) are shown next to the Arnold plots.
|
|
In all cases studied, for a given channel density the reliability at the preferred frequency was also higher than it would have been at this stimulus frequency for most other values of channel densities. Comparably high values were achieved only for channel densities where the frequency at the 1st harmonic or the subharmonic Arnold tongue coincided with the stimulus frequency.
We also analyzed the influence of 2 other potassium channels with slower kinetics on frequency preference of the reference cell—a muscarinic potassium channel KM and a calcium-dependent potassium channel KCa (for details see APPENDIX. For both cases, we substituted Ks by the new potassium conductance, KM or KCa, respectively. The results of the Arnold plot analysis are shown in Fig. 3B. For both channel types, an increase of their conductance shifted the preferred frequency over a broad range of frequencies. The lowest achievable frequency at a given DC depended on the time constant of the slow potassium conductance. If 2 or more slow potassium conductances were present at high densities, the broad tuning effect was diminished and eventually suppressed at high conductance levels (data not shown). Figure 3C presents the preferred frequency as a function of Ks conductance for different Ks. The slower the kinetics of the Ks channel, the lower the minimum achievable frequency and the broader the frequency range accessible through variation of the slow potassium conductance.
For completeness we analyzed all combinations of Na, NaP, Kdr, Ks, and leak conductances. In this case, we relied on the DC firing rate as an estimate of the preferred frequency. The distribution of maximum changes in firing rate (i.e., preferred frequency) achievable by variation of the density of one ion channel type over all combinations of the other 4 densities is presented in Fig. 4, which shows one curve for each ion channel type. For a more detailed description of this analysis see METHODS. Variation of Ks had a significant effect on the firing frequency in almost all parameter regimes. Its influence was weakest when another potassium channel, Kdr in this case, was present at high density. The mean change achieved with Ks was around 20 Hz. The mean change achieved by the other ion channels was <10 Hz.
|
20 Hz. Within the parameter space investigated, this was true for only a minority of values of the other 4 conductances. Figure 4B shows 4 examples of parameter regimes where these channels significantly changed the preferred frequency. For example, this occurred for Kdr when Ks was not present or present only in small amounts. NaP could cause a large frequency shift when both potassium conductances, Kdr and Ks, were low. Na was potent in changing the frequency when both potassium conductances and NaP were low. Its influence in these cases weakened further with a higher density of leak channels. Leak channel variation also gave rise to higher frequency shifts when both potassium conductances were low and the sodium conductances were not too large. In general, higher densities of leak channels tended to lower the minimum achievable frequency. To illustrate that regulation of ionic conductances on spike timing reliability frequency preference would allow a cell to dynamically adjust its spike timing reliability, the effect of a temporary increase in Ks conductance on spike timing reliability is presented in Fig. 5. The conductance step was chosen such that the preferred frequency of the cell after the conductance increase matched the stimulus frequency. Spike timing reliability during elevation of the Ks conductance was significantly enhanced.
|
|
EXPERIMENTAL RESULTS. To test the effects of slow potassium channels on preferred frequency physiologically, we performed patch-clamp recordings in slices of rat prefrontal cortex. We used the dynamic-clamp technique, which allows time-dependent currents to be injected that experimentally simulate conductances through on-line feedback. Thus we were able to artificially introduce Ks currents (with the same dynamics as the Ks reference channel used in the model simulations above). We first characterized its spike timing reliability as a function of frequency of a sine wave input in the "native" conductance state of the cell (the control state). Then we recorded responses to the same set of stimuli, in the presence of an additional Ks conductance in the cell, which was introduced with dynamic clamp (the Ks state). The DC level was kept constant and was the same for the control and the Ks experiment. Figure 7 shows spike timing reliability as a function of the frequency of the input for the control and Ks states. In the control experiment the cell had a preferred frequency at about 17 Hz. When Ks was introduced, the preferred frequency was lowered to 7 Hz and the spike timing reliability was increased beyond the maximum reliability of the control state. We observed a negative shift in the preferred frequency for all cells recorded (n = 4). For most cells, the shift in the preferred frequency resulted in an increase of spike timing reliability at that value (n = 3).
|
Our analysis supports that changes in ionic conductance can mediate changes in the preferred frequency through variation of the firing rate. However, we also found an effect of individual conductances on spike timing reliability not involving changes of the preferred frequency. We show that the overall level of reliability at the preferred frequency also depends on individual conductances, in addition to the general frequency preference discussed in the first part.
We created a set of cells with the same preferred frequency (at IDC = 0.3 nA), but different time constants of the Ks channel. The peak conductance gKs of each cell was adjusted such that the DC firing rate (i.e., preferred frequency) amounted to 20 Hz. All other conductance parameters were those of the reference cell. Sine wave stimulation of these cells confirmed that cells within one set had the same preferred frequency (Fig. 8A).
|
Ks [i.e., faster kinetics and higher values of gKs, Fig. 8, C and D]. This effect was more pronounced at the 1st harmonic of the preferred frequency, despite a small deviation from the general trend around Ks = 25 ms. The width also increased with faster kinetics and larger values of gKs. Figure 8, F and G show the voltage traces at the 1st harmonic for Ks values of 10 and 150 ms. The increase in reliability with larger and faster Ks conductance was accompanied by a deepening of the afterhyperpolarization following a spike (see Fig. 8) and a more efficient shutdown of sodium currents after a spike (not shown).
To distinguish the contribution of the peak conductance and the time constant of the slow potassium channel in this observation, 2 separate sets of neurons with a preferred frequency of 20 Hz were created. Cells in the 1st set varied in the Ks peak conductance. This time, the DC was adjusted to keep the firing rate at 20 Hz. Cells in the 2nd set differed in the time constant Ks. Again, adjustment of the DC kept the firing rate constant. The amplitude of the sine wave component of the input was the same for all cells.
Figure 9A shows a significant influence of Ks peak conductance on reliability at the preferred frequency for all 3 amplitude levels tested. Also, the width of the reliable region around the preferred frequency increased with gKs (Fig. 9D). For low conductance values, the curves were nonmonotonic. For these values, the difference between the peak reliability and that at surrounding frequencies was very small, which as a consequence broadened the estimate of width at the preferred frequency. The time constants also influenced spike timing reliability at the preferred frequency. Apart from a shallow maximum around Ks = 30 ms, the general trend confirmed that faster Ks channels increase the peak reliability and broaden the width. Overall, the effect of the Ks peak conductance was higher than that of the time constant.
|
|
|
DISCUSSION |
|---|
|
We first showed that neurons exhibit a preferred frequency for spike timing reliability, which correlates with the DC firing rate and can be regulated at the level of ion channels. Reliability is increased around the preferred frequency, but also around the harmonics and subharmonics of the preferred frequency. These findings confirm the resonance effect described by Hunter et al. (1998) in models of leaky integrate-and-fire neurons and in recordings from Aplysia motoneurons. We analyzed the influence of conductances in cortical cells and provided numerical and experimental support that slow potassium channels can regulate the preferred frequency. We showed that conductance-based regulation of reliability also extends to more realistic rhythm-like inputs that contain more than one frequency.
In a second step we studied the influence of ionic conductances and their kinetics on spike timing reliability beyond the location of the resonance effect of spike timing reliability in frequency space. We compared cells of the same DC firing rate and found that more and faster slow potassium channels enhance absolute reliability levels at the preferred frequency, whereas persistent sodium channels impair reliability.
Adjustment of preferred frequency through potassium channels
Analysis of pyramidal model neurons has revealed that the slow potassium conductance has a significant potential to tune the preferred frequency for spike timing reliability. Potassium currents counteract depolarization and set a limit to overall excitability. Controlling the interspike interval during repetitive firing (Hille 2001), they set the firing rate and thus preferred frequency. This observation was confirmed for a general model of a Ks channel, as well as a muscarinic potassium channel (KM) and a calcium-dependent potassium channel (KCa). Furthermore, we saw that the range of achievable firing rates through variation in Ks depended on the time constant of the channel. The faster the channel, the smaller the range of possible preferred frequencies.
Given the relation between preferred frequency and firing rate, it is not surprising that the other conductances are important too. First, they determine the DC firing rate if no slow potassium channels are present and thus set the upper limit of preferred frequencies that can be achieved by adjustment of Ks. Second, sodium leak and faster potassium channels also mediate larger shifts in frequency. Generally, the other channels extended their influence on preferred frequency with decreasing Ks conductance. The presence of (at least one) large potassium conductance, on the contrary, tended to dominate excitability, slower channels being more effective. Along those lines, the influence of Ks itself was weakest when the delayed rectifier conductance Kdr was strong (data not shown).
Pyramidal neurons in cortex show a strong adaptation, presumably mediated by a slow potassium conductance. Therefore we conclude that in these cells, the slow potassium conductance can efficiently regulate the spike timing reliability frequency preference. The firing rate analysis, however, depends on the parameter range of the pyramidal cell model. The experiments in slices of rat prefrontal cortex support the potential of the slow potassium conductance to shift frequency preference.
Potassium channels improve spike timing beyond a shift in preferred frequency
Analyzing neurons of the same preferred frequency but different Ks densities and kinetics revealed that ionic conductances can increase or decrease spike timing reliability without involving changes in overall excitability and preferred frequency. The effect was most pronounced at the preferred frequency and its harmonics. However, we also analyzed reliability in response to the DC alone and found a small effect. For a pure DC input, phase locking is not possible because the only time-varying signal is the noise, which is different in each trial and therefore randomizes spike timing for large times after stimulus onset. Nevertheless, spike timing reliability during the 1st second after stimulus onset (not counting the initial 500 ms) was significantly more precise for neurons with the higher and faster Ks conductance (data not shown). This suggests that the improvement of reliability with more and faster Ks channels results from a reduced sensitivity to noise after a spike. The deepening of the afterhyperpolarization and the more efficient shut down of sodium channels observed after a spike support this hypothesis. The separation of the effects of Ks time constant and Ks peak conductance showed that both parameters influence spike timing reliability, although the effect of peak conductance seemed somewhat more dominant. Variation of the time constant seemed to exhibit a small nonmonotonic influence on reliability for lower values, which will be subject to further investigation.
Potassium channels are promising candidates for regulating both the preferred frequency and the reliability at the preferred frequency. The trends for tunable frequency range (slower Ks channels are better) and the value of reliability at the preferred frequency (faster Ks channels are better) are opposed. Indeed, potassium channels are the ion channels with the greatest diversity, with gating kinetics that span several orders of magnitude (Hille 2001). They are regulated by a multitude of cellular signals and neuromodulators, examples of which are cAMP-dependent phosphorylation of slow potassium channels in the heart or regulation of muscarinic potassium channels through acetylcholine.
Functional importance of reliability modulation
In principle, the frequency tuning of spike timing reliability may be important for cells that are stimulated by inputs with dominant frequencies that do not change with time, such as frequency-tuned cells at the periphery of the auditory system. Differential expression of channel densities and different types of ion channels to match the prevalent input frequencies is not unknown, although it has not yet been shown to regulate spike timing reliability. One example for a different mechanism is the gradient of calcium and potassium channels in cochlear hair cells of lower vertebrates, which support electrical tuning to sound pressure waves of particular frequencies (Hille 2001). In addition, cells could in principle also regulate spike timing reliability dynamically. Important central signals are, for example, hippocampal, theta and gamma rhythms. There, cells may change their spike timing frequency preference through cellular messengers or neuromodulators on relatively short time scales. As our examples of more rhythm-like inputs suggest (Fig. 6), adjustment of spike timing frequency preference can also increase reliability for more natural types of inputs, such as theta- or gamma-type oscillations.
Future directions
Because the preferred frequency correlates with the DC firing rate of a cell, the preferred frequency does not only depend on intrinsic properties of a cell, but also on the DC level of the input, or in other words, the mean level of depolarization, as a recent study of Hunter and Milton (2003) has shown. Thus intrinsic spike timing preferences depend on the signal mean as well as synaptic input, reflecting the interplay between intrinsic properties and properties of the network. However, extrinsic changes of the DC such as modification of synaptic inputs can influence the value of the preferred frequency only according to the status of the intrinsic parameters of a cell. Changes of reliability at the preferred frequency cannot be achieved without tuning of intrinsic properties of a cell.
A recent investigation has shown that the intrinsic currents determine the state of synchrony in neural networks (Pfeuty et al. 2003). When the network is in a synchronous state and each neuron receives periodic input, an individual neuron in the network maintains synchrony by phase locking, which produces reliable spike timing. Interestingly, the results of Pfeuty et al. (2003) indicate that slow potassium currents promote synchrony in a network of conductance-based neurons, whereas persistent sodium channels impede it, which is in accordance with our findings of the influence of these channels on spike timing reliability. They also point out that synchrony varies with the firing rate depending on the intrinsic currents, which may be attributed to a resonance effect of synchrony, similar to that of spike timing reliability.
An aspect that we did not explore in the analysis is the possibility that channels may be regulated together, as a recent study of MacLean and colleagues (2003) has shown. If more than one channel type is regulated together (by gene expression or neuromodulation), a more efficient combination of conductances could be obtained.
The signals that we analyzed are suprathreshold mean signals. For lower levels of depolarization (signals with subthreshold mean) the prediction of the preferred frequency through the DC firing fails. It is probable that in this regime other intrinsic properties of cells, for instance determined by conductances that are activated at lower levels of depolarization, may shape spike timing frequency preferences of neurons. Channels like IH or IKA, which do not significantly influence suprathreshold signals (data not shown), are likely to gain importance.
In conclusion, regulation of conductances of a cell can efficiently affect spike timing reliability for inputs with dominant frequencies and suprathreshold mean. Changes in conductances can therefore be exploited to adjust spike timing reliability to the dominant input frequency. Moreover, ion channels also can influence spike timing reliability without changes to the mean firing rate. Spike timing reliability is the result of several mechanisms (of intrinsic and network origin), including the mean level of depolarization (Hunter and Milton 2003), which may be regulated through the background synaptic input, or well timed inhibitory or excitatory input (see e.g. Mehta et al. 2002).
|
|
APPENDIX |
|---|
|
Table A1 summarizes the conductance values underlying the normalized conductances (Fig. 3). The second column states the parameter values used for the firing rate analysis (Fig. 4A).
|
|
For most channels (Na, NaP, Kleak, Ks, and Kdr) the details are given in the work of Golomb and Amitai (1997). The details of the KM and channels KCa are given here.
The slow, noninactivating muscarinic potassium current KM was defined as
 |
n = –35 mV, n = 10 mV, n = –35 mV, 1n = 40 mV, 2n = 20 mV, and T = 36°C (cf. Gutfreund et al. 1995; Yamada et al. 1989). Values of n < 0.001 were set to 0.001. Exponentials with arguments larger than 50 were set to exp(50).
The calcium-dependent potassium current Kca was a nonvoltage-dependent potassium channel, activated by intracellular calcium. It was based on a first-order kinetic scheme
 |
M = 0.03 Hz (backward rate constant), Kca = 5.0e-04 mM (midactivation concentration), T = 36°C [for reference also see Destexhe et al. (1994)]. If M < 0.1 ms, then M was set to 0.1 ms.
The functionality of Kca requires intracellular calcium, whose dynamics include pumping of calcium ions and internal buffering, as well as a (voltage-dependent) calcium channel, CaL. We chose an L-type calcium channel with the following dynamics
 |
 |
ca = 0.1 mS/cm2 (maximum conductance of CaL), z = 5,040 ms, x = 1 mV, u = 11.3 mV, z = 18.9 mV. If abs (Vu) < 1e-04 mV then cu was set to cu = 0.00001 mV; otherwise cu = 0 mV. If abs(Vz) < 1e-04 mV then cz was set to cz = 0.00001 mV; otherwise cz = 0 mV.
The removal of intracellular calcium dynamics was accomplished by an ATPase pump
 |
= 2.4e-04 mM (equilibrium intracellular calcium concentration), and r = 1e10
[PDB]
ms (time constant of calcium removal by buffering). The time constant r was designed to model a first-order removal of intracellular calcium, which we neglected by setting r to a very large value. Kpump = 4e-04 mM/ms (time constant of the ATPase), Keq = 0.02 mM (equilibrium value for the ATPase), and F = 96485.309C mol–1 (Faraday constant). If wch was smaller than 0 then wch was set to 0 (no inward pumping) [for reference see Destexhe et al. (1993)]. |
|
ACKNOWLEDGMENTS |
|---|
|
GRANTS
This work was supported by the Daimler-Benz Foundation, the Sloan-Swartz Center for Theoretical Neurobiology, the German National Merit Foundation, and the Howard Hughes Medical Institute.
|
|
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: T. J. Sejnowski, Computational Neurobiology Laboratory, The Salk Institute, 10010 N. Torrey Pines Road, La Jolla, CA 92037 (E-mail: terry{at}salk.edu).
|
|
REFERENCES |
|---|
|
Beierholm U, Nielsen CD, Ryge J, Alstrom P, and Kiehn O. Characterization of reliability of spike timing in spinal interneurons during oscillating inputs. J Neurophysiol 86: 1858–1868, 2001.
Brette R and Guigon E. Reliability of spike timing is a general property of spiking model neurons. Neural Comput 15: 279–308, 2003.[CrossRef][Web of Science][Medline]
Connors B and Gutnick M. Intrinsic firing patterns of diverse neocortical neurons. Trends Neurosci 13: 365–366, 1990.[Web of Science][Medline]
Coombes S and Bressloff PC. Mode locking and Arnold tongues in integrate-and-fire neural oscillators. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Top 60: 2086–2096, 1999.[Medline]
Destexhe A, Babloyantz A, and Sejnowski TJ. Ionic mechanisms for intrinsic slow: oscillations in thalamic relay neurons. Biophys J 65: 1538–1552, 1993.[Web of Science][Medline]
Destexhe A, Contreras D, and Sejnowski TJ. A model of spindle rhythmicity in the isolated thalamic reticular nucleus. J Neurophys 72: 803–818, 1994.
Fellous JM, Houweling AR, Modi RH, Rao RP, Tiesinga PH, and Sejnowski TJ. Frequency dependence of spike timing reliability in cortical pyramidal cells and interneurons. J Neurophysiol 85: 1782–1787, 2001.
Fricker D and Miles R. EPSP amplification and the precision of spike timing in hippocampal neurons. Neuron 28: 559–569, 2000.[CrossRef][Web of Science][Medline]
Goldman MS, Golowasch J, Marder E, and Abbott LF. Global structure, robustness, and modulation of neuronal models. J Neurosci 21: 5229–5238, 2001.
Golomb D and Amitai Y. Propagating neuronal discharges in neocortical slices: computational and experimental study. J Neurophysiol 78: 1199–1211, 1997.
Gutfreund Y, Yarom Y, and Segev I. Subthreshold ocsillations and resonant frequency in guinea pig cortical neurons - physiology and modeling. J Physiol 483: 621–640, 1995.
Gutkin B, Ermentrout G, and Rudolph M. Spike generating dynamics and the conditions for spike-time precision in cortical neurons. J Comput Neurosci 15: 91–103, 2003.[CrossRef][Web of Science][Medline]
Haas JS and White JA. Frequency selectivity of layer II stellate cells in the medial entorhinal cortex. J Neurophys 88: 2422–2429, 2002.
Hille B. Ion Channels of Excitable Membranes (3rd ed.). Sunderland, MA: Sinauer Associates, 2001.
Hughes SW, Cope DW, and Crunelli V. Dynamic clamp study of Ih modulation of burst firing and delta oscillations in thalamocortical neurons in vitro. Neuroscience 87: 541–550, 1998.[CrossRef][Web of Science][Medline]
Hunter JD and Milton G. Amplitude and frequency dependence of spike timing: implications for dynamic regulation. J Neurophysiol 90: 000–000, 2003.
Hunter JD, Milton JG, Thomas PJ, and Cowan JD. Resonance effect for neural spike time reliability. J Neurophysiol 80: 1427–1438, 1998.
Jaeger D and Bower JM. Synaptic control of spiking in cerebellar purkinje cells: dynamic current clamp based on model conductances. J Neurosci 19: 6090–6101, 1999.
Jensen RV. Synchronization of randomly driven nonlinear oscillators. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Top 58: 6907–6910, 1998.
Keener J, Hoppensteadt F, and Rinzel J. Integrate-and-fire models of nerve membrane response to oscillatory input. SIAM J Appl Math 41: 503–517, 1981.[CrossRef]
Knight BW. The relationship between the firing rate of a single neuron and the level of activity in a population of neurons. Experimental evidence for resonant enhancement in the population response. J Gen Physiol 59: 767–778, 1972.
MacLean JN, Zhang Y, Johnson BR, and Harris-Warrick RM. Activity-independent homeostasis in rhythmically active neurons. Neuron 37: 109–120, 2003.[CrossRef][Web of Science][Medline]
Magee JC. A prominent role for intrinsic neuronal properties in temporal coding. Trends Neurosci 26: 14–16, 2002.[CrossRef][Web of Science]
Mainen Z and Sejnowski T. Reliability of spike timing in neocortical neurons. Science 168: 1503–1506, 1995.
Mainen Z and Sejnowski T. Influence of dendritic structure on firing pattern in model neocortical neurons. Nature 382: 363–366, 1996.[CrossRef][Medline]
Marder E, Abbott LF, Turrigiano GG, Liu Z, and Golowasch J. Memory from the dynamics of intrinsic membrane currents. Proc Natl Acad Sci USA 93: 13481–13486, 1996.
McCormick D, Connors B, Lighthall J, and Prince DA. Comparative electrophysiology of pyramidal and sparsely stellate neurons of the neocortex. J Neurophysiol 54: 782–805, 1985.
Mehta MR, Lee AK, and Wilson MA. Role of experience and oscillations in transforming a rate code into a temporal code. Nature 417: 741–746, 2002.[CrossRef][Medline]
Nowak LG, Sanchez-Vives MV, and McCormick DA. Influence of low and high frequency inputs on spike timing in visual cortical neurons. Cereb Cortex 7: 487–501, 1997.
Pfeuty B, Mato G, Golomb D, and Hansel D. Electrical synapses and synchrony: The role of intrinsic currents. J Neurosci 23: 6280–6294, 2003.
Reinagel P and Reid RC. Precise firing events are conserved across neurons. J Neurosci 22: 6837–6841, 2002.
Rescigno A, Stein R, Purple R, and Poppele R. A neuronal model for the discharge patterns produced by cyclic inputs. Bull Math Biophys 32: 337–353, 1970.[CrossRef][Web of Science][Medline]
Schreiber S, Fellous J-M, Whitmer D, Tiesinga PHE, and Sejnowski T. A new correlation-based measure of spike timing reliability. Neurocomputing 2003, in press.
Sharp AA, O'Neil MB, Abbott LF, and Marder E. Dynamic clamp: computer-generated conductances in real neurons. J Neurophysiol 69: 992–995, 1993.
Storm JF. Potassium currents in hippocampal pyramidal cells. Prog Brain Res 83: 161–187, 1990.[Web of Science][Medline]
Tanabe M, Gahwiler BH, and Gerber U. L-type Ca2+ channels mediate the slow Ca2+-dependent afterhyperpolarization current in rat CA3 pyramidal cells in vitro. J Neurophysiol 80: 2268–2273, 1998.
Tiesinga PHE. Precision and reliability of periodically and quasiperiodically driven integrate-and-fire neurons. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Top 65: 041913, 2002.
Tiesinga PHE, Fellous J-M, and Sejnowski TJ. Attractor reliability reveals deterministic structure in neuronal spike trains. Neural Comput 14: 1629–1650, 2002.[CrossRef][Web of Science][Medline]
Turrigiano G, Abbott LF, and Marder E. Activity-dependent changes in the intrinsic properties of cultured neurons. Science 264: 974–977, 1994.
Yamada Y, Nakazato Y, and Ohga A. Ouabain distinguishes between nicotinic and muscarinic receptor-mediated catecholamine secretions in perfused adrenal-glands of cat. Br J Pharmacol 96: 470–479, 1989.[Web of Science][Medline]
This article has been cited by other articles:
|
|
 |
|
  S. Schreiber, I. Samengo, and A. V.M. Herz Two Distinct Mechanisms Shape the Reliability of Neural Responses J Neurophysiol, May 1, 2009; 101(5): 2239 - 2251. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. A. Prescott and T. J. Sejnowski Spike-Rate Coding and Spike-Time Coding Are Affected Oppositely by Different Adaptation Mechanisms J. Neurosci., December 10, 2008; 28(50): 13649 - 13661. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
T. A. Engel, L. Schimansky-Geier, A.V.M. Herz, S. Schreiber, and I. Erchova Subthreshold Membrane-Potential Resonances Shape Spike-Train Patterns in the Entorhinal Cortex J Neurophysiol, September 1, 2008; 100(3): 1576 - 1589. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 D. A. Markowitz, F. Collman, C. D. Brody, J. J. Hopfield, and D. W. Tank Rate-specific synchrony: Using noisy oscillations to detect equally active neurons PNAS, June 17, 2008; 105(24): 8422 - 8427. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. F. M. van Brederode and A. J. Berger Spike-Firing Resonance in Hypoglossal Motoneurons J Neurophysiol, June 1, 2008; 99(6): 2916 - 2928. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
R. F. Galan, G. B. Ermentrout, and N. N. Urban Optimal Time Scale for Spike-Time Reliability: Theory, Simulations, and Experiments J Neurophysiol, January 1, 2008; 99(1): 277 - 283. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
S. E. Street and P. B. Manis Action Potential Timing Precision in Dorsal Cochlear Nucleus Pyramidal Cells J Neurophysiol, June 1, 2007; 97(6): 4162 - 4172. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 A. D. Dorval The Rhythmic Consequences of Ion Channel Stochasticity Neuroscientist, October 1, 2006; 12(5): 442 - 448. [Abstract] [PDF]
|
|
|
|
 |
|
 J.-M. Goaillard and E. Marder Dynamic Clamp Analyses of Cardiac, Endocrine, and Neural Function Physiology, June 1, 2006; 21(3): 197 - 207. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C. P. Billimoria, R. A. DiCaprio, J. T. Birmingham, L. F. Abbott, and E. Marder Neuromodulation of spike-timing precision in sensory neurons. J. Neurosci., May 31, 2006; 26(22): 5910 - 5919. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. J. Lawrence, Z. M. Grinspan, J. M. Statland, and C. J. McBain Muscarinic receptor activation tunes mouse stratum oriens interneurones to amplify spike reliability J. Physiol., March 15, 2006; 571(3): 555 - 562. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. Bauer, R. Oostenveld, M. Peeters, and P. Fries Tactile Spatial Attention Enhances Gamma-Band Activity in Somatosensory Cortex and Reduces Low-Frequency Activity in Parieto-Occipital Areas J. Neurosci., January 11, 2006; 26(2): 490 - 501. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
V. Zsiros and S. Hestrin Background Synaptic Conductance and Precision of EPSP-Spike Coupling at Pyramidal Cells J Neurophysiol, June 1, 2005; 93(6): 3248 - 3256. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
R. K Powers, Y Dai, B. M Bell, D. B Percival, and M. D Binder Contributions of the input signal and prior activation history to the discharge behaviour of rat motoneurones J. Physiol., February 1, 2005; 562(3): 707 - 724. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
K. Diba, H. A. Lester, and C. Koch Intrinsic Noise in Cultured Hippocampal Neurons: Experiment and Modeling J. Neurosci., October 27, 2004; 24(43): 9723 - 9733. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
R. Corlew, M. M Bosma, and W. J Moody Spontaneous, synchronous electrical activity in neonatal mouse cortical neurones J. Physiol., October 15, 2004; 560(2): 377 - 390. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J.-M. Fellous, P. H. E. Tiesinga, P. J. Thomas, and T. J. Sejnowski Discovering Spike Patterns in Neuronal Responses J. Neurosci., March 24, 2004; 24(12): 2989 - 3001. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |
TOP
ABSTRACT
INTRODUCTION
