JN Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 90: 3888-3901, 2003; doi:10.1152/jn.00477.2003
0022-3077/03 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (13)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lee, S. C.
Right arrow Articles by Ishida, A. T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lee, S. C.
Right arrow Articles by Ishida, A. T.

Availability of Low-Threshold Ca2+ Current in Retinal Ganglion Cells

Sherwin C. Lee, Yuki Hayashida and Andrew T. Ishida

Section of Neurobiology, Physiology, and Behavior, University of California, Davis, California 95616-8519

Submitted 19 May 2003; accepted in final form 31 July 2003


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
Spiking in central neurons depends on the availability of inward and outward currents activated by depolarization and on the activation and priming of currents by hyperpolarization. Of these processes, priming by hyperpolarization is the least described. In the case of T-type Ca2+ current availability, the interplay of hyperpolarization and depolarization has been studied most completely in expression systems, in part because of the difficulty of pharmacologically separating the Ca2+ currents of native neurons. To facilitate understanding of this current under physiological conditions, we measured T-type current of isolated goldfish retinal ganglion cells with perforated-patch voltage-clamp methods in solutions containing a normal extracellular Ca2+ concentration. The voltage sensitivities and rates of current activation, inactivation, deactivation, and recovery from inactivation were similar to those of expressed {alpha}1G (CaV3.1) Ca2+ channel clones, except that the rate of deactivation was significantly faster. We reproduced the amplitude and kinetics of measured T currents with a numerical simulation based on a kinetic model developed for an {alpha}1G Ca2+ channel. Finally, we show that this model predicts the increase of T-type current made available between resting potential and spike threshold by repetitive hyperpolarizations presented at rates that are within the bandwidth of signals processed in situ by these neurons.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
Low-threshold, T-type Ca2+ current (IT) displays transient activation by depolarization when initiated from relatively hyperpolarized membrane potentials and inactivation during sustained depolarization (for review, Huguenard 1996Go; Perez-Reyes 2003Go). This combination of properties enables IT to influence the shape and timing of action potentials in several ways. The most well known of these is low-threshold spiking in which activation of IT triggers a brief burst of action potentials by transiently depolarizing cells. Second, the activation and steady-state inactivation profiles of IT often overlap near the resting potential of cells, forming a "window current" that can augment cell excitability (Williams et al. 1997Go). Third, IT activated during normal action potentials can both shape spikes (Lambert et al. 1998Go) and influence spike frequency (Wolfart and Roeper 2002Go).

Although a T-type current underlies all of these effects, variations in time course and pharmacology in different systems suggested that there were multiple isoforms of IT (Huguenard 1996Go). This was confirmed by the recent isolation of three T-type Ca2+ channel clones, termed CaV3.1, CaV3.2, and CaV3.3 (or {alpha}1G, {alpha}1H, and {alpha}1I, respectively, in the older nomenclature) (for review, Perez-Reyes 2003Go). Whole cell and single-channel recordings of these clones in expression systems have shown that they differ substantially in their rates of activation, inactivation, and recovery from inactivation, while continuing to demonstrate the voltage dependence that characterizes this channel family. These studies of expressed channels have provided an explanation for the variation in kinetic properties and pharmacology found in studies of natural IT (Huguenard 1996Go), and it has been possible, in a few instances, to match a specific cloned channel with the native channel in a particular cell (Chemin et al. 2002Go; Satin and Cribbs 2000Go). Still, accounting for the availability and influence of IT during spiking in particular cell types requires measurements of these properties from native channel populations under meticulous physiological conditions.

The best studied of the native systems are the various classes of thalamocortical cells that exhibit low-threshold spiking. IT was shown to support the low-threshold spike, both directly and in numerical simulations (Coulter et al. 1989Go; Crunelli et al. 1989Go; Huguenard and McCormick 1992Go; Huguenard and Prince 1992Go; McCormick and Huguenard 1992Go; Suzuki and Rogawski 1989Go; Wang et al. 1991Go). This and other work (Bal et al. 1995Go; Thomson 1988Go) have demonstrated that the ability of IT to support this behavior depends critically on the interaction of cell membrane potential and the timing and frequency of hyperpolarizing inhibitory postsynaptic potentials (IPSPs).

To determine the properties that underlie this function and to examine the role of IT in another type of central neuron, we report here on the voltage-dependent control and kinetics of IT in retinal ganglion cells, a class of central neuron that has displayed this current in every species examined to date (rat: Guenther et al. 1994Go; Karschin and Lipton 1989Go; turtle: Liu and Lasater 1994Go; goldfish: Bindokas and Ishida 1996Go; cat: Huang and Robinson 1998Go; salamander: Henderson and Miller 2003Go). Activation threshold and inactivation kinetics, as well as some pharmacology, have been described for the IT in retinal ganglion cells of postnatal rat (Guenther et al. 1994Go; Karschin and Lipton 1989Go), turtle (Liu and Lasater 1994Go) and goldfish (Bindokas and Ishida 1996Go), but the other biophysical properties of these currents have yet to be studied. To ensure that we obtained the basal properties of the current, we used isolated cells in short-term culture to avoid the influence of surrounding cells; we used perforated-patch methods to maintain cytoplasmic integrity; and we measured currents in normal physiological Ca2+ concentration to avoid surface charge effects and augmented current amplitudes. Under these conditions, we have found in goldfish retinal ganglion cells that the voltage sensitivities and rates of current activation, inactivation, and recovery from inactivation resemble those of T-type currents from a number of other tissues, and that deactivation is two to three times faster. We could reproduce the amplitude and kinetics of IT activated by typical depolarizing voltage jumps with a kinetic model developed for an {alpha}1G (CaV3.1) Ca2+ channel clone. Moreover, we show that this model predicts the increase of IT that can be recorded between resting potential and spike threshold after volleys of hyperpolarizations, as might occur during recurrent inhibitory input.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
The voltage-clamp currents described here were measured in retinal ganglion cell somata isolated from adult common goldfish (Carassius auratus). The methods used in this study have been described in detail elsewhere (Bindokas and Ishida 1996Go; Hidaka and Ishida 1998Go; Tabata and Ishida 1996Go, 1999Go). All animal care and experimental protocols were approved by the Animal Use and Care Administrative Advisory Committee of the University of California, Davis. Cells were dissociated by brief incubation of freshly dissected retinas in either protease-containing saline or, in a few instances, an enzyme-free, low-Ca2+ solution (Hayashida, Partida, and Ishida; unpublished), followed by thorough rinsing in saline supplemented with bovine serum albumin and trituration with a large-bore glass pipette. Cells were identified either by retrograde labeling, via the optic nerve, of rhodamine B isothiocyanate-coupled dextran (Tabata and Ishida 1996Go) or by the nucleolar expansion observed exclusively in ganglion cells subsequent to crush of the optic nerve (Ishida and Cohen 1988Go). We found no statistically significant difference between these cell groups in terms of the current properties reported here (activation threshold, kinetics, voltage sensitivity of steady-state inactivation, rate of recovery from inactivation). Therefore all data from these cells have been pooled. Voltage-gated Ca2+ current was recorded in perforated-patch mode (Horn and Marty 1988Go), using an Axopatch-1D patch-clamp amplifier (Axon Instruments, Union City, CA) and pCLAMP software (v. 8.1.01, Axon Instruments) for voltage-jump protocol generation and data acquisition.

Patch electrodes were pulled from borosilicate glass capillaries to tip resistances of ca. 5 M{Omega}. The tips of these pipettes were filled with either of two "pipette solutions." One formulation (ca. 60% of experiments) contained (in mM) 15 Na methanesulfonate, 120 CsOH, 30 tetraethylammonium chloride (TEA), 1 ethyleneglycoltetraacetic acid (EGTA), 0.34 CaCl2, 2.6 MgCl2, and 5 HEPES; pH adjusted to 7.40 with methanesulfonic acid. The second formulation (ca. 40% of experiments) was the same except that TEA was replaced with 30 CsOH, 15 sucrose, and additional methanesulfonic acid to adjust pH. We observed no differences between the Ca2+ currents recorded with these pipette solutions. The osmolality of both solutions was 285–290 mosmol/kg. The pipette shanks were filled with the same solution supplemented with a 1:200 dilution of a mixture of 2 mg amphotericin B and 3 mg Pluronic F-127 (Molecular Probes; Eugene, OR) in 60 µL DMSO.

The extracellular ("bath") solution contained (in mM) 120 NaCl, 3 CsCl, 2.5 CaCl2, 1 MgCl2, 30 TEA-Cl, 10 glucose, and 5 HEPES. The pH was adjusted to 7.40 with NaOH, and the osmolality was 290–300 mosmol/kg. Just prior to use, the bath solution was supplemented with 1 µM tetrodotoxin.

Reagents were obtained from the following sources: Na methanesulfonate, methanesulfonic acid (Aldrich Chemical, Milwaukee, WI); CsOH, CsCl (ICN, Aurora, OH); NaCl, TEA, MgCl2, EGTA, HEPES, amphotericin B, and DMSO (Sigma Chemical, St. Louis, MO); CaCl2 (BDH; Dorset, UK); and tetrodotoxin (Sankyo, Tokyo).

Experiments were performed at 20–24°C within 36 h of cell isolation. Recordings were made from cells adhering to a glass coverslip attached to the bottom of a 35-mm tissue culture dish. Cells recorded the day of isolation or after overnight culture gave equivalent results; there was no drift in T-current density or properties with time in culture. The current monitor output of the patch-clamp amplifier was analog-filtered by the built-in 4-pole Bessel filter (Fc = 2 or 5 kHz) and digitally sampled (usually at 10 kHz). An agar bridge was used to ground the bath solution, and the membrane potentials reported here have been corrected for liquid junction potentials that arose from differences between the bath and pipette solutions. Linear leak subtraction was performed with the P/N protocol in pCLAMP, using holding and command potentials that were never more negative than –100 mV, and never more positive than –60 mV. Series resistance compensation was not used.

Data were analyzed in pCLAMP and Excel (Excel 2000; Microsoft, Redmond, WA), and data curves were fitted in SigmaPlot (version 5.0.5; SPSS, Chicago, IL). Means are presented as ±SE.

Computer simulations were carried out using the simulation software NEURON (ver. 4.3.1 and 5.2 by J. W. Moore, M. Hines, and T. Carnevale) (see Hines and Carnevale 1997Go, 2000Go). Because the dendrites and axons of our retinal ganglion cells were typically shorn off by the dissociation procedures, leaving just the isolated somata, we modeled the ganglion cell as a single spherical compartment, 20 µm in diameter.

Current flow through the T-channel in the open state (IT, units of A/cm2) was calculated with an equation

as used in other studies (e.g., De Schutter and Smolen 1998Go; Hille 1984Go; Huguenard and McCormick 1992Go). Here, PCa is the maximum permeability (units of cm/s); O is the fraction of channels in the open state (see following text); z is the valence (+2); [Ca2+]in and [Ca2+]out are Ca2+ concentrations inside and outside the cell; Vm is the membrane voltage; F, R, and T are the Faraday constant, gas constant, and absolute temperature, respectively. In this equation, the concentration-dependent rectification of current is described by the Goldman-Hodgkin-Katz constant field equation (Lewis 1979Go). Because our experiments utilized physiological concentrations of internal Na+ (15 mM) and extracellular Ca2+ (2.5 mM) and were confined to a limited range of negative membrane potentials, we could ignore monovalent cation permeability through the T channel (Fukushima and Hagiwara 1985Go; Lux et al. 1990Go). Similarly, as we used the same physiological concentration of external Mg2+ (1 mM) throughout, we did not explicitly incorporate the relatively minor effect of Mg2+ block on the T channel (Fukushima and Hagiwara 1985Go; Serrano et al. 2000Go). In the present simulations, an estimate of 7 x 10-6 cm/s was used for PCa to produce 5 pA/pF of current density (see RESULTS); and [Ca2+]in and [Ca2+]out were assumed to be constant at 100 nM and 2.5 mM, respectively. Within the time scale of our simulations, the amplitude and kinetics of calculated currents were not markedly altered by changes in [Ca2+]in between 10 nM and 1 µM (data not shown).

The fraction of channels in the open state was calculated using a Markovian kinetic model shown in Fig. 8A. The kinetic scheme was adapted from one developed for a cloned T-channel (Serrano et al. 1999Go). In brief, the channel occupies any of 12 states (5 closed states, 1 open state, and 6 inactivated states). As membrane potential depolarizes, the probability of channel occupancy is higher for the states in the right-hand side of figure. The transitions through closed states (C1 to C5) and most of the inactivated-inactivated transitions (I1 to I5) are voltage dependent. The closed-open transition (C5 to O) and one inactivated-inactivated transition (I5 to IO) are voltage-independent. The closed/open-inactivated transitions (C1–5/O to I1-O; parallel vertical transitions in figure) are also voltage independent. However, for the first three closed-closed transitions (C1 to C4), the corresponding inactivation process is coupled allosterically to activation (cf. Kuo and Bean 1994Go). Thus inactivation develops with membrane depolarization, and saturates as the channel reaches the last two closed states and open state.



View larger version (38K):
[in this window]
[in a new window]
 
FIG. 8. Kinetic model and simulated T current. A: Markovian kinetic scheme developed originally by Serrano et al. (1999Go). C1 to C5 are the closed states; O is the open state; I1 to IO represent the inactivated states. kV and k–V are functions of membrane voltage, and the other transition rates (kO, k–O, kI, and k–I) as well as the allosteric coupling factors (f and h) are constants. B: voltage-dependent activation and decay of inward current elicited by voltage steps from –92 mV to –62 through –37 mV, as simulated by the model (black). The corresponding superimposed experimental data were taken from Fig. 1A (gray). Time constant for the command voltage was set to 0.7 ms in the simulation. Model parameters were the same as Table 1 except k-V0 was 0.015. C: activation and deactivation of T current calculated by the model (black). The corresponding superimposed data were taken from Fig. 3A (gray). The time constant for the command voltage was set to 0.25 ms in the simulation. The effect of this delay is shown by the thin black trace (marked, *) that shows simulated current for instantaneous command steps. The model parameters were the same as Table 1 except k-V0 was 0.039.

 



View larger version (31K):
[in this window]
[in a new window]
 
FIG. 1. IT currents. Ca2+ currents elicited by 100-ms depolarization from a holding potential (Ehold) of –92 mV. The voltage steps were increased in 5-mV increments starting at –62 mV, as shown in the upper part of the figure. P/N subtraction was used to remove linear leak and capacitive currents. Inset: peak current vs. test voltage for these data ({bullet}). —, Boltzmann fit, described in the text. The peak current was normalized to the value obtained at –37 mV.

 

View this table:
[in this window]
[in a new window]
 
TABLE 1. Backward and forward rates of activation

 


View larger version (25K):
[in this window]
[in a new window]
 
FIG. 3. IT deactivation kinetics. A: IT decay currents activated by the voltage protocol shown (top). Each trace is the average of 3 trials with P/N leak subtraction. A period approximately equal to the charging time of this cell (0.25 ms for 30 M{Omega} access x 8.4 pF capacitance, perforated patch) was blanked at the beginning of the repolarization. The traces were staggered to avoid overlap, and a sample exponential fit is shown in the lower right (–92 mV). Fitting was done in Clampfit. B: average deactivation time constants as a function of voltage. —, a simple exponential fit.

 
In this model, there are a total of 10 free parameters to be estimated. The constraints on parameter estimation were as follows: 1) kV0, the rate of voltage sensor movement for activation, was selected to match the time-to-peak at a depolarized potential (–37 mV, Fig. 2A); 2)k–V0, the rate of voltage sensor movement for deactivation, was selected to reproduce V1/2 for the current activation curve (Fig. 5B); 3) VkV, Vk–V, the gating charge parameters associated with voltage sensor movement, were selected to generate the slope factor for the I-V relation (Fig. 5B); 4) kO, the channel opening rate, was set to match the time-to-peak at an intermediate potential (–52 mV, Fig. 2A); 5) k–O, the channel closing rate, was determined directly from the time constant of deactivation at a hyperpolarized potential (Fig. 3B); 6) kI, the inactivation rate, was taken directly from the time constant of the current decay at a depolarized potential, extrapolated from Fig. 2B; 7) k–I, the reverse rate of inactivation, was set to reproduce V1/2 for the steady-state inactivation curve (Fig. 5B); 8) f, the allosteric coupling factor for inactivation, was selected to yield the time constant of entry to inactivation measured near resting potential (Fig. 4C); 9) h, the allosteric coupling factor for recovery from inactivation, was selected to reproduce the time constant of recovery from inactivation measured at a hyperpolarized potential (Fig. 7B). When simulations were repeated to optimize the parameter values, the values of k–O and kI were fixed and kV0, k–V0, VkV, Vk–V, and kO were adjusted first, and then k–I, f and h were tuned. The iteration and annealing were carried out manually until all constraints were reasonably satisfied.



View larger version (12K):
[in this window]
[in a new window]
 
FIG. 2. IT activation and inactivation kinetics. A: time to peak current as a function of test voltage. B: time constant of decay of the current as a function of test voltage as determined from a simple exponential fit. {bullet}, mean; error bars represent ±1 SE; n = 19 cells. —, simple exponential fits.

 


View larger version (26K):
[in this window]
[in a new window]
 
FIG. 5. IT steady-state inactivation and window current. A: steady-state inactivation: inward currents elicited by 100-ms voltage steps to –42 mV from holding potentials that varied from –92 to –62 mV. Ehold was changed in 5-mV increments starting at –92 mV (top). Each new membrane potential was held for 15 s before the test step to –42 mV. P/N subtraction was used to remove linear leak and capacitive currents. Inset: peak current vs. holding potential for this cell ({circ}). —, Boltzmann fit to these data as described in the text. Peak current at each holding potential is normalized to the value obtained by the test step from –92 mV. B: IT window current. Activated current (peak I-V): {bullet}, the average normalized peak current vs. test voltage in 9 cells with the voltage protocol shown in Fig. 1. The peak current elicited by the step to –37 mV was used for normalization for each cell. Error bars represent ±1 SE but are mostly covered by the data markers in this figure. —, Boltzmann fit as described in the text. Residual current (steady-state inactivation): {circ} average normalized peak current at one fixed test potential (–42 mV) vs. holding potential for the same 9 cells. The peak current produced from Ehold –92 mV was used for normalization of the data from each cell. Window current is the overlap region of the extrapolated Boltzmann fits and represents a voltage region in which IT could be persistently active.

 


View larger version (15K):
[in this window]
[in a new window]
 
FIG. 4. Inactivation of IT at –62 mV. A: overlaid raw current traces produced by the voltage protocol presented in the upper part of this figure. From Ehold (–92 mV), a control current was elicited by a 50-ms test step to –42 mV; after 2.5 s at –92 mV, voltage was stepped to –62 mV for varying lengths of time (/ /) before another 50-ms test step to –42 mV to check for residual current. Capacitive currents were clipped from these traces manually rather than by P/N subtraction. B: the currents elicited by the test step to –42 mV, at higher resolution, after the indicated times at –62 mV. C: plot of the peak current of each test step divided by its matched control vs. time at –62 mV ({bullet}) for this cell. —, a simple exponential fit; the time constant for onset of inactivation in this cell was 103 ms.

 


View larger version (19K):
[in this window]
[in a new window]
 
FIG. 7. Recovery from depolarized potential IT inactivation. A: overlaid currents produced by the voltage protocol (top). From Ehold (–92 mV), IT was measured and subsequently inactivated by a 100-ms test step to –42 mV; cell voltage was then stepped back to –92 mV for varying lengths of time (/ /) before testing the recovery of IT with another test step to –42 mV for 30 ms. With P/N subtraction. B: plot of the peak current activated by each test step, divided by its matched control vs. time at –92 mV ({bullet}); this is the average of 3 identical protocols, 1 of which is shown in A. —, a simple exponential fit ({tau} 121 ms) described in the text.

 

When comparing the model to a specific data set, time constants for the rise and fall of command voltage steps were set equal to the charging time constant measured from that cell. In all simulations, the temperature and the calculation time step were set to 23°C and 2–5 µs, respectively.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
We describe here the following properties of low-threshold Ca2+ current (IT) in retinal ganglion cells: voltage dependence and kinetics of activation and deactivation, voltage sensitivity of steady-state inactivation, rate of inactivation near resting potential, and rate of recovery from inactivation. To measure these properties, current was identified as IT on the basis of three features: 1) it was activated by step-wise depolarizations from holding potentials more negative than –65 mV, to test potentials more positive than –65 mV (Fig. 1); 2) it inactivated markedly, if not entirely, when the holding potential was shifted to values more positive than –60 mV (Fig. 5A); and 3) it was kinetically transient, i.e., it decayed during test depolarizations of 100 ms (Fig. 1). This T-type current is greatly reduced by lowering extracellular Ca2+ (data not shown), and it can be abolished by replacement of the extracellular Ca2+ with a mixture of 2.4 mM Co2+ and 0.1 mM Ca2+ (Bindokas and Ishida 1996Go). It is insensitive to the spider toxin {omega}-Aga-IIIA (Bindokas and Ishida 1996Go), so it is unlikely that any of it is the R-type Ca2+ current that sometimes shows overlapping voltage dependence and kinetic properties (Randall and Tsien 1997Go). Most (ca. 75%) of the cells recorded from in this study exhibited low-threshold Ca2+ current meeting these criteria, as we have reported previously (Bindokas and Ishida 1996Go).

All of the currents reported here were recorded under conditions designed to minimize contamination by other voltage-gated currents, optimize recording stability, and avoid space clamp artifacts due to neurites. Voltage-gated Na+ current was blocked by inclusion of 1 µM tetrodotoxin in the bath (Hidaka and Ishida 1998Go); hyperpolarization-activated cation current (Ih) was blocked by the presence of 3 mM Cs+, and the absence of K+, in the bath (Tabata and Ishida 1996Go); outward K+ currents were minimized by use of a Cs+-based, K+-free pipette solution and by inclusion of 30 mM TEA in the bath (Tabata and Ishida 1999Go); activation of high-threshold Ca2+ current and of TEA-resistant Cs+ efflux were avoided by use of sufficiently negative test potentials (Bindokas and Ishida 1996Go; Tabata and Ishida 1996Go, 1999Go); and leak-like Cl current was minimized by use of perforated-patch recording mode (Tabata and Ishida 1999Go). Perforated-patch mode was also used to avoid drifts in the voltage sensitivity of gating that occurred during ruptured-patch recording (Munckton, Pignatelli and Ishida, unpublished observations) and allay concerns about possible "rundown" (Wan et al. 1996Go). We restricted our investigation to voltages negative to –35 mV to avoid contamination by high-threshold Ca2+ current (see following text). Recordings showing signs of inadequate space clamp (in particular, delayed activation or distorted I-V relations) were discarded from the data set. For this reason, we did not study cells bearing neurites longer than 10 µm or so, and we did not attempt to compare the properties of currents recorded from neurite-free cell bodies versus cells with significant amounts of neurites. Instead, we limited our analyses to the amplitude and kinetics of currents where we felt our control of membrane potential was the best possible. As such, our study is concentrated on somatically expressed T current, and we cannot say anything about currents derived from other cell locations (Ahlijanian et al. 1990Go; Baldridge 1996Go; Henderson and Miller 2003Go). Even with all these precautions, about 1/2 of all cells with low-threshold current were discarded because of the presence of slower kinetic components in the voltage range we used to characterize IT. This left 103 cells that formed the basis for the following analysis.

Activation range

The range of membrane potentials that activate IT was measured by the voltage-jump protocol shown at the top of Fig. 1. Holding potential was set to the most negative value at which stable recordings were routinely obtained (-92 mV), and cells were depolarized to test potentials between –62 and –37 mV. Current versus voltage (I-V) curves were constructed by plotting the peak amplitude of the Ca2+ current versus the test potential at which the Ca2+ current was activated. To facilitate comparison of data collected from different cells, current amplitudes were normalized to the peak amplitude obtained in response to the step to –37 mV recorded in each cell. The normalized I-V relation for the data in Fig. 1 is plotted in the inset and fitted by Marquardt-Levenberg regression to the Boltzmann equation

where I is peak current at each test voltage, Imax is peak current at –37 mV, V is the test voltage, V1/2 is the midpoint of the Boltzmann fit, and Vc is the voltage for an e-fold change around V1/2 (i.e., slope factor). For this cell, V1/2 is –53 mV and Vc is 4.3 mV.

In most cells, we observed the onset of a slowly inactivating inward component at test potentials more positive than ca. –37 mV. This portion of the whole cell current was considered to be distinct from low-threshold Ca2+ current because it was high threshold and was only partially reduced in amplitude by depolarized holding potentials (see Bindokas and Ishida 1996Go). Because high-threshold Ca2+ currents were larger when elicited from a holding potential of –92 mV than from –62 mV, even in cells with no demonstrable IT, we could not subtract high-threshold current elicited from –62 mV to remove contamination, by high-threshold current, from IT elicited from –92 mV. Moreover, we found previously that there was not a satisfactory mix of blockers that could cleanly isolate IT from the other Ca2+ currents in goldfish retinal ganglion cells (Bindokas and Ishida 1996Go). Hence, we exploited the differential voltage dependence of IT activation to assess IT properties by using currents collected only at command potentials that did not elicit noticeable amounts of high-threshold Ca2+ current, as in other studies of native cells (e.g., Biagi and Enyeart 1991Go; Carbone and Lux 1987Go). In most cells, the current elicited by the step to –37 mV is probably not quite the maximum for IT. The Boltzmann fit in Figs. 1 and 5B suggests the true maximum current would occur between –35 and –30 mV. A more direct measure of channel activation would be to normalize peak current by the instantaneous current (e.g., Herrington and Lingle 1992Go; Serrano et al. 1999Go). However, our inability to isolate IT across a wider voltage range and the inherent limitations of perforated-patch recording prevented us from determining instantaneous currents. Compared with Figs. 1 and 5B, we would expect the channel activation curve to be broader and show attainment of maximum activation at a slightly more positive potential, as shown by the simulated activation curve corrected for the calculated Ca2+ driving force (dashed line) in Fig. 9C (see following text). These curves also indicate that the error introduced by our use of the peak I-V relation is small at the hyperpolarized potentials where both channel activation and total T current are small.



View larger version (19K):
[in this window]
[in a new window]
 
FIG. 9. Inactivation profiles of simulated T current. A: entry to inactivated state near the resting potential simulated by the model (black). The corresponding experimental data were taken from Fig. 4B (light red). Time constant for the command voltage was set to 0.4 ms in the simulation. B: recovery from open-state inactivation simulated by the model (black). The corresponding data were taken from Fig. 7A (light red). The time constant for the command voltage was set to 0.33 ms in the simulation. For both A and B, the model parameters were those in Table 1, except that kI was 0.082, f was 0.224, and h was 0.196. C: steady-state inactivation and normalized I-V curves, simulated by the model (black lines). The corresponding data were replotted with marks (light red) from Fig. 5B. Model parameters were as given in Table 1. A simulated channel activation curve, corrected for the calculated Ca2+ driving force, is shown by the dashed line.

 

At the test potentials used here (between –67 and –37 mV), IT ranged in amplitude up to ~200 pA. However, large currents were unusual, and the median peak current at –37 mV was 35 pA [44 ± 3.3 (mean ± SE) pA; n = 103]. Division of the maximum current amplitude recorded in each cell by the membrane capacitance, yielded current densities of 3–5 pA/pF [5 ± 0.7 pA/pF for cells identified by retrogradely transported dextran (n = 9); 3 ± 0.3 pA/pF (n = 94) for cells identified morphologically (see METHODS)].

Because IT was small in most cells, background noise could interfere with our ability to measure the voltage at which current first activated. Typically, we needed ca. 3 pA of inward current to detect this point of threshold. To ensure the sensitivity of this determination, we confined our analysis of activation threshold to cells with larger peak current (>60 pA at –37 mV), such that the current required to identify the threshold was not more than 5% of the peak. The activation threshold in such cells was –61 ± 0.5 mV (n = 21). This is consistent with the I-V relation for the cell in Fig. 1 and for the average I-V relation of nine other well-characterized cells shown in Fig. 5B (for these 9 cells: peak current –66 ± 9 pA at –37 mV; measured threshold –60 ± 0.8 mV). [The Boltzmann fit to the I-V plot for these cells meets the abscissa at around –70 mV (Fig. 5B), suggesting that small amounts of T current might activate at membrane potentials more negative than the measured threshold.] The test potential that elicited half-maximum current amplitudes (V1/2) in these nine cells was –51 ± 0.5 mV, and near this voltage the current amplitude grew e-fold in amplitude every 4 ± 0.2 mV (Vc; see Fig. 5B).

Activation, inactivation, and deactivation kinetics

The rate at which IT rose in amplitude was measured by the time-to-peak (time between the onset of the command depolarization and the current peak) as a function of voltage as shown in Fig. 2A. This yielded a smoothly increasing rate of rise with stronger depolarization. The rate at which IT fell in amplitude at these test voltages was estimated by the exponential time constant that fit the declining phase of currents (Fig. 2B). Stronger depolarizations produced faster decay of current until reaching a plateau at ca. –40 mV. The extrapolated limiting time constant for this inactivation was 15.4 ms.

To determine the rate of deactivation, IT was activated by a brief step to –42 mV, and the repolarization-induced decay of the current was recorded at potentials from –62 to –92 mV (Fig. 3A). Because deactivation proved to be a fairly rapid process, we restricted our analysis to cells in which the charging time of the cell (approximated by the product of access resistance and cell capacitance) was <0.5 ms. To reduce the access resistance and obtain the best quality measurements, some of these data were acquired in ruptured patch configuration immediately after patch rupture and before any significant change in the voltage dependence of IT properties. Figure 3B was generated from pooled perforated and ruptured patch data (7 ruptured and 5 perforated measurements as in A in 5 cells). Deactivation was faster as the final membrane voltage was made more negative. The average deactivation time constant at –92 mV was 0.82 ms.

Rate of inactivation near resting potential

Little or no T-type current is available for activation if the holding potential is more positive than –65 mV. Because this is near the resting potential of spiking neurons (including retinal ganglion cells; see DISCUSSION), we measured the rate at which IT inactivates (i.e., becomes unavailable for activation) at these voltages. To do this, we depolarized cells from a holding potential of –92 mV to a conditioning potential equal to a typical resting potential (-62 mV) using the voltage-jump protocol shown at the top of Fig. 4. To assess the rate of inactivation at this voltage, we measured the amount of Ca2+ current available for activation after conditioning depolarizations of various durations by depolarizations to a fixed test potential (e.g., Fig. 4, A and B). The amplitude of these currents was normalized to the amplitude of current activated from –92 mV and plotted against the duration of the conditioning depolarization to –62 mV. In all of three cells tested, plots of this type were fitted by an exponential decline in amplitude with an average time constant of 146 ± 25 ms (e.g., Fig. 4C).

Steady-state inactivation

The voltage sensitivity of inactivation was measured by the voltage-jump protocol shown in Fig. 5A, top. The holding potential was shifted to various values between –92 and –62 mV, and the amount of Ca2+ current available for activation was measured by depolarizing cells from each holding potential to a test potential of either –52 or –42 mV. Each holding potential was maintained for 15 s prior to initiating the test depolarizations, i.e., for a period well in excess of the inactivation time constant measured in Fig. 4. The current at each test potential from each holding potential was normalized to the current elicited by depolarizations from –92 mV to the same test potential, and fitted by least-squares regression to the Boltzmann equation (Fig. 5A, inset).

For the data in Fig. 5A, the conditioning potential that reduced the Ca2+ current to half-maximum amplitude (i.e., V1/2) was –74 mV. Near this voltage, the current amplitude fell e-fold every 3.9 mV. We found no dependence of the normalized currents on whether the holding potential was varied from –92 to –62 mV or from –62 to –92 mV. The average steady-state inactivation of nine cells (including this one) is plotted in Fig. 5B. Because we found no difference between the normalized currents collected at a test potential of –52 versus –42 mV, the average steady-state inactivation curve in Fig. 5B includes normalized current amplitudes that were measured at both test potentials. The average V1/2 for steady-state inactivation was –77 ± 0.9 mV and the average slope factor was 4 ± 0.2 mV. In three cells, IT was assessed from a holding potential of –102 mV. Currents activated by depolarizations from –102 mV were only ~3% larger than those activated in the same cells from –92 mV. The normalization error resulting from the use of –92 mV for the maximum current was therefore considered to be negligibly small.

Currents activating at membrane potentials more negative than the voltages that produce complete steady-state inactivation are termed "window" currents (Attwell et al. 1979Go). Superimposing the average peak I-V relation on the average steady-state inactivation curve provides an estimate of the voltage range for a window-type Ca2+ current in retinal ganglion cells (Fig. 5B). In the cells we recorded from, the overlap between these two curves was maximal at –64 mV. At this voltage, ~5% of IT is not inactivated (see DISCUSSION).

Recovery from inactivation

The normal resting potential of these retinal ganglion cells is in the range of –60 to –75 mV (Vaquero et al. 2001Go; C. F. Vaquero and A. T. Ishida, unpublished results); similar values have been reported for carp, salamander, turtle, guinea pig, and cat retinal ganglion cells of a variety of functional subtypes (Baylor and Fettiplace 1979Go; Coleman and Miller 1989Go; Demb et al. 2001Go; Mittman et al. 1990Go; Murakami and Shimoda 1977Go; O'Brien et al. 2002Go; Wiesel 1959Go). In this voltage region, only a few percent of IT is not inactivated. To estimate how quickly the availability of IT can increase, we measured the rate at which it recovers from inactivation at a hyperpolarized potential. Because this required us to compare test current amplitudes against control amplitudes, we used two different protocols to guard against effects of the control depolarization on the rate of recovery itself. In both, the holding potential was set to –92 mV, and IT was elicited by a test depolarization that started at various times after the end of a conditioning depolarization. In the first method, shown schematically at the top of Fig. 6A and with representative data in Fig. 6, A and B, the conditioning potential was –62 mV (to mimic resting potential). The amplitude of the test IT was compared with that of IT elicited by a control depolarization presented before the conditioning depolarization. The second method, shown schematically in Fig. 7A, top, and with representative data in Fig. 7A, was the more widely used "two-pulse protocol," in which the conditioning potential was set equal to the test potential, and the conditioning depolarization was used as the control depolarization. The rate of recovery from inactivation was estimated by normalizing the amplitude of current elicited by the test depolarization to that elicited by the control depolarization, and plotting this amplitude against the duration of the time elapsed between the end of the conditioning depolarization and the beginning of the test depolarization. Recovery, measured with either of the methods, followed an exponential time course. In the example data shown here, the time constants of recovery were 76 ms (Fig. 6C) and 121 ms (Fig. 7B). Overall, there was no difference between the recovery time course after inactivation at –62 mV (119 ± 11 ms, n = 8) or more depolarized potentials (112 ± 9 ms, n = 9). In three cells in which both protocols were performed, the time constant was 116 ± 9 ms when inactivation was induced by depolarizations to –62 mV versus 104 ± 11 ms when inactivation was induced by depolarizations to –52 or –42 mV. Combining all data obtained by both protocols gives an average time constant for recovery of 115 ± 7 ms (n = 17).



View larger version (18K):
[in this window]
[in a new window]
 
FIG. 6. Recovery from resting potential IT inactivation. A: overlaid currents elicited by the voltage protocol (top). From Ehold (–92 mV), a control current was generated by a 50-ms test step to –47 mV; after 2.5 s at –92 mV, the potential was set at –62 mV for 3.0 s and then stepped back to –92 mV for varying lengths of time (/ /) before testing the recovery of IT with another 50-ms test step to –47 mV. These traces are the average of 2 identical protocols performed in this cell, with P/N subtraction. B: at higher resolution, the currents elicited by the test step to –47 mV after the indicated times at –92 mV. C: plot of the peak current activated by each test step, divided by its matched control vs. time at –92 mV ({bullet}) for this cell. —, simple exponential fit({tau} 76 ms), described in the text.

 

We generally observed that the recovery from inactivation was well-described by a single exponential time constant. In some other preparations, the recovery from inactivation of IT is best described by two exponentials, with the faster component being most comparable to that observed here (Huguenard 1996Go). The amount of the slower component is increased by extending the length of the inactivating voltage step. However, we did not see evidence for a second recovery component in experiments comparing inactivation due to a 0.1 s step to –42 mV versus 3 s at –62 mV. We did not test the effect of holding the potential for many seconds at depolarized potentials, and we explicitly avoided any such time dependence in the determination of steady-state inactivation by holding at each potential for a long time.

Kinetic model for T current in isolated retinal ganglion cells

The preceding experiments characterized the fundamental kinetic properties of the goldfish retinal ganglion cell IT. To assess whether these properties are consistent with those established for other T currents and to provide a tool for predicting the behavior of IT in different circumstances, we developed a numerical model for IT using the simulation program NEURON (see METHODS). Previous studies have shown that Hodgkin-Huxley models are valid as practical descriptions of whole cell IT (Huguenard and McCormick 1992Go; Wang et al. 1991Go; see also Destexhe and Huguenard 2000Go), although they are not mechanistically accurate, especially in their description of inactivation processes. The recent cloning of T channels allowed detailed studies of their biophysical properties (Cribbs et al. 1998Go; Lee et al. 1999aGo; Perez-Reyes et al. 1998Go) and provided sufficient experimental data to develop plausible kinetic models of their behavior (Burgess et al. 2002Go; Frazier et al. 2001Go; Serrano et al. 1999Go). A similar kinetic model has also been applied to the IT of a native neuron (Kuo and Yang 2001Go). In the following sections, we show that, with one specific modification, the Markovian kinetic model previously developed by Serrano et al. (1999Go) for a cloned T channel can generate satisfactory fits to our experimental data and that these fits extend over a variety of conditions.

Figure 8A shows the kinetic scheme for gating of IT used in the present simulation. The behavior of this model is determined by 10 parameters (Table 1). The numerical values of some of these parameters were estimated directly from the present experimental data, whereas others were derived from iterative optimization of the overall model within the constraints imposed by our data (see METHODS for the fitting procedure). The qualitative characteristics of these parameters are the same as those described by Serrano et al. (1999Go; especially their Fig. 14) except that the transition rate from O to C5 (k–O) was assumed to be voltage independent. This assumption reduced the number of free parameters in the model, while yielding satisfactory fits to our experimental data (see following text). Fits between the simulations and experimental data did not improve when voltage-dependent rates of k–O [e.g., k–O = k–O0 exp(Vm/Vk–O)] were used in the simulations (results not shown). Hence, only the forward and backward rates of activation (kV and k–V) are voltage-dependent variables, and the other transition rates (kO, k–O, kI and k–I), as well as the allosteric coupling factors (f and h), are constants.

Voltage-clamp simulations

To evaluate the model, calculated T currents were compared with actual data from the present experiments. For each data set, the time course of the current elicited by the appropriate voltage protocol was simulated, with minor adjustment of the values in Table 1 for each individual cell.

Figure 8 shows simulations of the activation and decay of IT during sustained depolarization (Fig. 8B) and the deactivation of IT after repolarization (Fig. 8C). The simulation results (black) superimpose well on the corresponding experimental data from Figs. 1 and 3, respectively (reproduced in gray). For Fig. 8, B and C, the value of k–V0 was adjusted slightly, whereas the other parameters were fixed at the values in Table 1. Such adjustments could affect other properties of the current, but the changes were self-consistent. For example, although the adjustment applied in the simulation in Fig. 8C caused a shift in V1/2 for the I-V relation from –51 to ca. –46 mV for the simulated data, this value of V1/2 is close to the value measured experimentally in this particular cell (ca. –47 mV, data not shown).

The time course of entry into (Fig. 9A) and recovery from inactivation (Fig. 9B) were also reproduced well by the model. Figure 9, A and B, shows the simulation results (black) corresponding to the experimental data from Figs. 4 and 7, respectively (light red). For these, the values of kI, f, and h were adjusted, and the same set of parameter values was used in both simulations. A small discrepancy between the simulation and experiment was seen in the amplitude of current available after 32 ms of recovery in B. This might be due to an initial delay in the recovery phase, noted by others (Burgess et al. 2002Go; Kuo and Yang 2001Go; Satin and Cribbs 2000Go). Simulation of the other voltage protocol we employed (see Fig. 6A) reproduced almost the same time constant for recovery from inactivation (results not shown). The model also predicts that the apparent rate of recovery from inactivation will be slower at less hyper-polarized potentials. In particular, the model says the time constant for recovery will be 1.8 times longer at –77 mV than –92 mV. In two cells where such a comparison was made, the rate of recovery at –77 mV was 115 ± 7 versus 78 ± 2 ms at –92 mV, 1.5 times longer.

Figure 9C shows steady-state inactivation and I-V curves as generated by the model (black lines), corresponding to the average experimental data in Fig. 5B (light red marks). The values used were those in Table 1 without modification. As shown in Fig. 9C, simulated curves for both inactivation and I-V gave reasonable fit to the experimental data and reproduced the "window" current region. A simulated channel activation curve is shown by the dashed line.

IT availability

Because low-threshold current is almost completely inactivated at resting potential, the availability of IT under physiological conditions will depend on hyperpolarizing inputs that bring about relief from inactivation (see DISCUSSION below). Because ganglion cells are buffeted by a constant stream of excitatory and inhibitory inputs, and real hyperpolarizing influences are almost certainly transient (e.g., Sakai and Naka 1987Go), the availability of IT will depend not only on the magnitude of the net voltage shift but also on the rates of deactivation and inactivation near the resting potential and the rate of recovery from inactivation at hyperpolarized potentials. We therefore asked if the model utilized in this study, which reproduces the amplitude and time course of currents elicited by stepwise depolarizations, can be extended to describe the currents activated after a sequence of brief hyperpolarizations.

Figure 10 presents a comparison of simulated and experimental data examining the availability of IT after cycling the membrane potential from –62 to –82 mV a variable number of times. The membrane potential was held at –62 and –82 mV for the same amount of time during each conditioning cycle (10, 30, or 50 ms in this example), and then IT was activated by a step from –82 to –42 mV as shown schematically in Fig. 10, top. The frequencies of these conditioning cycles (generally 5–50 Hz) were chosen to be consistent with the time course of psychophysically measured temporal contrast sensitivity (Woodhouse and Barlow 1982Go), and the kinetics of GABA inhibitory postsynaptic currents (Protti et al. 1997Go). The simulated data (middle) predict that the amount of IT activated by the step to –42 mV increases with additional cycles until a maximum value is reached and that this value increases as the steps in the conditioning cycle are lengthened. When the step length is shorter, many more conditioning steps are needed to reach 90% of this limiting value, but the cumulative conditioning time to 90% is similar for each step length.



View larger version (33K):
[in this window]
[in a new window]
 
FIG. 10. Priming and relief from inactivation. The availability of IT was assessed after a variable number of repetitive hyperpolarizing steps from –62 to –82 mV, as shown schematically (top). The membrane potential was set alternately at –82 and –62 mV for equal times (10, 30, or 50 ms) during each cycle; and the final test step was from –82 to –42 mV. Middle: simulations of the normalized Prediction for the test current as a function of the number of conditioning cycles (number shown next to trace). The traces labeled "Experiment" are currents recorded in response to the same protocol. The time axis represents the cumulative length of the protocols. The conditioning cycles were produced with the conditioning train option in Clampex, which does not record the currents elicited by the conditioning protocol. However, other experiments in which "conditioning" currents were saved showed little or no activation of IT by this procedure. The time constant for the command voltage was set to 1.48 ms in the simulation. The model parameters were those in Table 1, except that k-V0 was 0.018, k-I was 0.13 x 10–3, f was 0.224, and h was 0.125. The prediction was made after adjusting the model parameters to reproduce the voltage dependence of activation and inactivation for this particular cell, based on its measured activation and inactivation kinetics, rate of recovery from inactivation, and activation and steady-state inactivation curves. The experimental priming data were not used for this parameter adjustment.

 

Representative experimental data obtained by this protocol are shown in the bottom row of Fig. 10. There is good correspondence between the prediction and experiment traces in each column, suggesting that the model satisfactorily implements the "priming" of IT by recurrent membrane potential oscillation. We recorded such data from three cells with equivalent results. As a fraction of the peak current elicited by a step to –42 mV from a holding potential of –82 mV, the average limiting availability of IT for these cells increased from 0.29 ± 0.01 (n = 2) to 0.59 ± 0.03 (n = 3) for conditioning steps of 10 and 100 ms, respectively.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
We have characterized and modeled the biophysical behavior of IT in goldfish retinal ganglion cells, with particular attention to the voltage dependence and kinetics of activation, inactivation and recovery from inactivation. In the following text, we compare the kinetics of ganglion cell IT to that of the T-channel clones and discuss the role and interaction of membrane potential and IT in ganglion cells at rest and during light responses.

Comparison with cloned T-channel currents

We found that the rate of activation was strongly voltage dependent, that the rate of inactivation reached a plateau with strong depolarization, and that deactivation and recovery from inactivation were fast. These are general characteristics of IT in many different native systems (Huguenard 1996Go) as well as the three primary cloned channels that produce T-type currents in expression systems (Perez-Reyes 2003Go). The kinetic parameters of IT of the goldfish retinal ganglion cell are plotted together with those of cloned rat and human T channels in Fig. 11. Goldfish IT kinetic parameters are comparable to these other values, although goldfish values tend to be among the fastest, especially for deactivation. This could reflect some unique property or accessory subunit association of the goldfish T channel (cf. Hobom et al. 2000Go), or it might simply reflect conditions used to make the different studies. Although the measurements shown in Fig. 11 were made under reasonably consistent ionic conditions, almost all the studies were done at room temperature, which is cold for the mammalian channels, but normal or a bit warm for the goldfish. Because the Q10's of T-type current properties are generally >2 (Coulter et al. 1989Go; Takahashi et al. 1991Go), mammalian T channels at appropriate physiological temperature would gate substantially faster, and their deactivation kinetics would be closer to that of the goldfish IT. The kinetics of both the goldfish and mammalian IT would then be comparable.



View larger version (22K):
[in this window]
[in a new window]
 
FIG. 11. Kinetic parameters summary. Compares kinetic parameters obtained for IT from goldfish retinal ganglion cells with those obtained for rat and human T-channel clones in expression systems. {circ}, values determined for {alpha}1G clones; {square}, for {alpha}1H; {triangleup}, for {alpha}1I. , the values for goldfish IT. {tau} (h) is the time constant for inactivation, {tau} (recov) the time constant for recovery from inactivation (at ca. –90 mV), t-t-p is time-to-peak (at ca. –40 mV), and {tau} (deact) is the time constant for deactivation (at ca. –90 mV). Note log scale for ordinate. Clone data were taken from Chemin et al. (2001Go, 2002Go), Cribbs et al. (2000Go), Klöckner et al. (1999Go), Monteil et al. (2000Go), Satin and Cribbs (2000Go), and Serrano et al. (1999Go).

 

Figure 11 shows that the kinetic parameters of the IT subtypes tend to cluster in a characteristic manner. The parameters obtained for goldfish IT show significant overlap with the {alpha}1G and {alpha}1H subtypes but not with {alpha}1I. Between the overlapping subtypes, two observations suggest that the goldfish retinal ganglion cell IT is of the {alpha}1G subtype: the kinetics of recovery from inactivation are markedly faster for IT in this work and for {alpha}1G than for {alpha}1H (Fig. 11) and the {alpha}1H subtype is blocked by Ni2+ at low concentrations (IC50 13 µM) (Lee et al. 1999bGo), but both cloned {alpha}1G (IC50 250 µM) (Lee et al. 1999bGo) and this IT (38% block at 100 µM) (Bindokas and Ishida 1996Go) are relatively insensitive to Ni2+. In the only other retinal ganglion cells where Ni2+ sensitivity has been assessed, T-type current was reduced ~85% by 20 µM Ni2+ in Xenopus (Akopian and Witkovsky 1996Go). This higher sensitivity to Ni2+ suggests that this T current is produced by a different channel subtype, but the recording conditions of this study were sufficiently different that a direct comparison might not be valid.

Consistent with this similarity of T-current properties between the retinal ganglion cell and the cloned {alpha}1G subtype (Fig. 11), most of the parameter values of the kinetic model used in our simulations were roughly near the values for cloned {alpha}1G (compare Table 1 vs. the legend for Fig. 14 of Serrano et al. 1999Go). The rate constants for channel opening and closing were significantly different from those for {alpha}1G; these parameters needed to be larger to account for the relatively rapid deactivation of T current in the retinal ganglion cell. The same kinetic model has been applied successfully to the {alpha}1I T-type channel (Frazier et al. 2001Go). Although this subtype displays faster deactivation kinetics than {alpha}1G (Fig. 11), its other rate constants are substantially slower (Frazier et al. 2001Go). A simulation using rate constants similar to {alpha}1I exhibits an order of magnitude slower increase in IT availability during priming and yields a much poorer fit to our data (not shown). Overall, the kinetics and pharmacology of ganglion cell IT are in agreement with the properties of the {alpha}1G subtype, although we cannot exclude that other subtypes may be present.

Physiological role

IT may contribute to setting the resting [Ca2+] level. In almost all cell types expressing IT, including the goldfish retinal ganglion cell, the extrapolated window current coincides with the approximate resting potential (Huguenard 1996Go). This suggests that a small residual Ca2+ influx due to an IT window could help set resting intracellular [Ca2+]. There are two different approaches we can use to estimate the Ca2+ influx due to persistently active IT in retinal ganglion cells. The "window" current arising from the overlap of the I-V and steady-state inactivation curves would be on the order of 0.1 pA. [i.e., 0.05 (fraction of noninactivated IT at –64 mV) x 0.02–0.05 (approximate open probability of IT at –64 mV) x 51 pA (mean peak IT adjusted for driving force at –64 mV)] in goldfish retinal ganglion cells. Although a very small current, this would still bring in the equivalent of ~150 nM/s homogeneously distributed [Ca2+] in a 20 µm diameter cell (~12.5 pF, modeled as a sphere). This "window" current calculation formally assumes the independence of activation and inactivation. On the other hand, in the kinetic model we used to simulate our results (Serrano et al. 1999Go), activation and inactivation are not independent. In this scheme, inactivation is incomplete and 0.1–0.2% (simulated with our parameter values) of T channels remain open at all voltages where IT activates. Although different in mechanism and voltage dependence from the "window" current, the magnitude of persistent Ca2+ current calculated by the kinetic model is very similar to the window at –64 mV (data not shown). In either case, with measured resting intracellular [Ca2+]of ~120 nM (Ishida et al. 1991Go), a continuous Ca2+ influx of this magnitude could contribute in setting [Ca2+] levels in goldfish retinal ganglion cells.

Another possibility is that T window (or persistent) current could contribute directly to a form of input signal amplification as described in thalamocortical neurons (Williams et al. 1997Go). However, this phenomenon depends on IT being relatively large compared with the leak, and satisfying the relation that the maximum slope of dIT/dV in the vicinity of the window is greater than the leak conductance (Williams et al. 1997Go). Estimated leak conductance in goldfish retinal ganglion cells is 0.9 nS (Lee and Ishida, unpublished). A generous estimate of window dIT/dV is 0.1 pA/5 mV = 0.02 nS, well short of the threshold for this effect. It therefore seems unlikely that the amounts of T current we have recorded could contribute to excitability by this mechanism in retinal ganglion cells. Cells with more T current might be different.

Voltage- and time-dependent changes in the availability of IT might also play a role in spike generation. Two issues must be considered when assessing the possible contribution that IT could make. First, do ganglion cells hyperpolarize enough to prime IT? Second, after priming, how does IT compare in amplitude with other subthreshold currents in ganglion cells?

Because the resting membrane potential of isolated goldfish retinal ganglion cells is in the range of –60 to –75 mV (Vaquero et al. 2001Go; C. F. Vaquero and A. T. Ishida, unpublished data), and because most IT is inactivated in this voltage range (Fig. 5), a mechanism to hyperpolarize ganglion cells, and thereby relieve the inactivation, must come into play for IT to contribute to excitability. Intracellular recordings have shown that suitably arranged stimuli hyperpolarize ganglion cells in various species by as much as 15 mV beyond their resting potential (Ammermüller and Kolb 1995Go; Slaughter and Bai 1989Go), and recent recordings have indicated that the reversal potential for these voltage changes is around –95 mV (e.g., Demb et al. 2001Go; Zaghloul et al. 2003Go). Assuming