JN Watch the video to see how APS reaches out to developing nations.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 91: 346-357, 2004. First published September 17, 2003; doi:10.1152/jn.00062.2003
0022-3077/04 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
91/1/346    most recent
00062.2003v1
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 (19)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Komendantov, A. O.
Right arrow Articles by Canavier, C. C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Komendantov, A. O.
Right arrow Articles by Canavier, C. C.

A Modeling Study Suggests Complementary Roles for GABAA and NMDA Receptors and the SK Channel in Regulating the Firing Pattern in Midbrain Dopamine Neurons

Alexander O. Komendantov1, Olena G. Komendantova2, Steven W. Johnson3 and Carmen C. Canavier2

1Center for Computational Science, Tulane University, New Orleans 70118; 2Department of Psychology, University of New Orleans, New Orleans, Louisiana 70148; and 3Department of Neurology, Oregon Health Sciences University, Portland, Oregon 97201

Submitted 22 January 2003; accepted in final form 8 September 2003


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 REFERENCES
 
Midbrain dopaminergic (DA) neurons in vivo exhibit two major firing patterns: single-spike firing and burst firing. The firing pattern expressed is dependent on both the intrinsic properties of the neurons and their excitatory and inhibitory synaptic inputs. Experimental data suggest that the activation of N-methyl-D-aspartate (NMDA) and GABAA receptors is a crucial contributor to the initiation and suppression of burst firing, respectively, and that blocking Ca2+-activated potassium SK channels can facilitate burst firing. A multi-compartmental model of a DA neuron with a branching structure was developed and calibrated based on in vitro experimental data to explore the effects of different levels of activation of NMDA and GABAA receptors as well as the modulation of the SK current on the firing activity. The simulated tonic activation of GABAA receptors was calibrated by taking into account the difference in the electrotonic properties in vivo versus in vitro. Although NMDA-evoked currents are required for burst generation in the model, currents evoked by GABAA-receptor activation can also regulate the firing pattern. For example, the model predicts that increasing the level of NMDA receptor activation can produce excessive depolarization that prevents burst firing, but a concurrent increase in the activation of GABAA receptors can restore burst firing. Another prediction of the model is that blocking the SK channel current in vivo will facilitate bursting, but not as robustly as blocking the GABAA receptors.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 REFERENCES
 
Midbrain dopaminergic (DA) neurons play an important role in voluntary movement, reward mediated learning, working memory, cognition, and motivation (Kitai et al. 1999Go; Schultz 1998Go; Spanagel and Weiss 1999Go). They are also implicated in neuropsychiatric disorders including schizophrenia (Weinberger 1987Go), drug abuse (Koob et al. 1987Go), and Parkinson's disease (Obeso et al. 2000Go). This relatively homogenous population of neurons (but see Liss et al. 2001Go; Neuhoff et al. 2002Go) is located in the ventral tegmental area (VTA or A10), the substantia nigra pars compacta (SNC or A9), and the retrorubral area (A8).

Midbrain DA neurons in vivo exhibit at least two major firing patterns: single-spike firing, which is often irregular, and burst firing (Freeman et al. 1985Go; Grace and Bunney 1984aGo,bGo; Tepper et al. 1995Go; Wilson et al. 1977Go). Activation of N-methyl-D-aspartate (NMDA) receptors has been implicated in burst firing in vivo (Chergui et al. 1993Go). The mechanisms by which the firing pattern is regulated are of interest because the firing pattern is likely to encode information that is transmitted to target neurons. Transitions between burst firing and single-spike firing have been observed in single-neuron recordings from behaving rats (Freeman et al. 1985Go). The transition to bursting appeared to be correlated with an orienting response to a sensory stimulus (Freeman et al. 1985Go). When spikes are clustered into bursts, the increases in levels of extracellular dopamine in the projection area are much larger per spike than those observed for regularly spaced trains of action potentials at the same average frequency due to nonlinear summation as release outpaces uptake (Wightman and Zimmerman 1990Go). Furthermore, a bursting pattern of stimulation of the medial forebrain bundle increased the expression of one immediate early gene (cfos) (Chergui et al. 1996Go) and increased the mRNA expression of several immediate early genes (Chergui et al. 1997Go) in the dopamine-innervated brain areas of the rat, whereas a regular single-spike stimulation pattern did not.

In contrast to the activity exhibited in vivo, dopamine neurons in vitro generally fire in a regular single spike-firing pattern. Bath application of NMDA can induce burst firing in vitro under some circumstances (Johnson and North 1992Go), and increasing the activation of GABAA receptors in vitro by applying the appropriate agonist can convert NMDA-induced burst firing to single spike firing (Johnson and Canavier 1998Go; Paladini et al. 1999bGo). Blocking the calcium-activated small conductance (SK) potassium channel facilitates bursting in vitro (Seutin et al. 1993Go) and can putatively induce burst firing in the absence of NMDA (Nedergaard et al. 1993Go; Ping and Shepard 1996Go) as well. Under other conditions, an SK channel blocker such as apamin can induce irregular single-spike firing (Ping and Shepard 1996Go).

Spontaneous synaptic potentials from excitatory amino acid (EAA) and GABAergic inputs are not usually observed in the slice preparation, although spontaneous events that appear to be GABAA-mediated postsynaptic potentials (PSPs) can be observed if extracellular potassium concentration is increased (Johnson and North 1992Go). Thus the major difference between dopamine neurons in vivo versus in vitro is that the synaptic afferents are intact in vivo, leading to a much lower input resistance as measured with sharp electrodes [192 ± 12 (SE) M{Omega} (Johnson and North 1992Go), 183 ± 5 (SE) M{Omega} (Ping and Shepard 1996Go), 152 ± 14 (SE) M{Omega} (Scroggs et al. 2001Go), 292 ± 11 (SE) M{Omega} (Seutin et al. 1997Go), and 117.5 ± 4.8 (SD) M{Omega} (Yung et al. 1991Go) versus 31.2 ± 7.4 (SD) M{Omega} (Grace and Bunney 1983Go)], and a shorter time constant [33 ± 2 (SE) ms (Johnson and North 1992Go) and 37.5 ± 1.9 (SD) ms (Yung et al. 1991Go) versus 12.1 ± 3.2 (SD) ms (Grace and Bunney 1983Go)]. We have assumed that the difference in input resistance is primarily due to the tonic activation of the GABAA receptors in vivo because 70-90% of all the inputs to the DA neurons are GABAergic (Bolam and Smith 1990Go; Smith and Bolam 1990Go). The main sources of GABAergic inputs to DA neurons originate from the substantia nigra pars reticulata, striatum, and pallidum, which all exert their effects via GABAA receptors (Celada et al. 1999Go).

The application of the GABAA receptor antagonists such as picrotoxin and bicuculline strongly promotes burst firing in vivo (Paladini and Tepper 1999Go). However, bicuculline also blocks the apamin-sensitive SK channel (Johnson and Seutin 1997Go), and blockade of this channel has also been recently shown to promote burst firing in vivo (Seutin et al. 2002Go). We used a multi-compartmental model of a DA neuron with a realistic morphological structure and electrotonic properties to explore the effects of varying the levels of activation of NMDA and GABAA receptors as well as the level of modulation of the SK current on the firing pattern because experimental data suggest that all three play a major role. Burst firing mechanisms (and oscillatory mechanisms in general) require a fast positive feedback process and a slow negative feedback process (Rinzel and Ermentrout 1998Go). As in previous models (Canavier 1999Go; Komendantov and Canavier 2002Go; Li et al. 1996Go), we assumed that the fast positive feedback process is the voltage-dependent opening and closing of NMDA receptor channels in the presence of tonic activation of those receptors, whereas the slow negative feedback is provided by the sodium activation of the sodium pump resulting from changes in sodium concentration in the dendrites due to sodium entry via NMDA receptors and sodium removal by the pump itself (Johnson and North 1992Go). This mechanism is controversial, but we believe it is the best available explanation of NMDA-induced burst firing in vitro and have incorporated it into our in vivo model as well, although we acknowledge that it may not be the sole bursting mechanism exhibited by these neurons (Amini et al. 1999Go). A preliminary report of our findings has been published in abstract form (Komendantov and Canavier 2001Go).


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 REFERENCES
 
Experimental procedures

TISSUE PREPARATION. Male Sprague-Dawley rats (150-300 g; Bantin and Kingman, Seattle, WA) were anesthetized with halothane and killed by severing major thoracic vessels. The whole brain was quickly removed and submerged in ice-cold artificial cerebral fluid (ACSF). A vibrating microtome (Lancer) was used to cut slices of midbrain (300 µm) in the horizontal plane. A slice was then placed on a nylon mesh in a recording chamber (volume: 0.5 µl) and secured by two titanium electron microscopy grids that were held in place by small pieces of platinum wire (0.5 mm diam). The slice was submerged in continuously flowing (2 ml/min) ACSF of the following composition (in mM): 126 NaCl, 2.5 KCl, 1.2 NaH2PO4, 1.2 MgCl2, 2.4 CaCl2, 10 glucose, and 18 NaHCO3; gassed with 95% O2-5% CO2, pH 7.4, at 36°C.

INTRACELLULAR RECORDINGS. Microelectrodes were made from borosilicate capillary tubing (1.0 mm OD, 0.5 mm ID; Dagan) using a P-97 Flaming-Brown micropipette puller (Sutter Instrument, Novato, CA). Microelectrodes were filled with either 2 M potassium chloride or potassium acetate with resistances of 45-150 M{Omega}. Membrane voltage was amplified with an Axoclamp-2B amplifier (Axon Instruments, Foster, CA) and recorded on an IBM-compatible personal computer running Axotape (Axon Instruments) software. To ensure accurate measurement of voltage, recordings were made with an active bridge circuit that was frequently checked for proper balance by passing small (50 pA) current steps while voltage output was monitored on an oscilloscope. Membrane input resistance was calculated by measuring the change in membrane potential in response to small (20-100 pA) hyperpolarizing current pulses.

IDENTIFICATION OF DOPAMINE NEURONS. Using a dissection microscope, the VTA was identified as the region lateral to the fasciculus retroflexus and medial to the medial terminal nucleus of the accessory optic tract, whereas the SNC was located immediately rostral and caudal to the medial terminal nucleus of the accessory optic tract (Paxinos and Watson 1986Go). Neurons were identified as dopaminergic using well-established electrophysiological and pharmacological criteria (Grace and Onn 1989Go; Lacey et al. 1989Go; Yung et al. 1991Go). Briefly, dopamine-containing neurons were identified as such by their broad action potentials (>2 ms), spontaneous pacemaker-like firing pattern (1-5 Hz), by the "sag" in membrane potentials recorded during hyperpolarizing current pulses (signifying H current), and by 5- to 15-mV hyperpolarization evoked by superfusing the slice with dopamine (30 µM). Because we found no differences in results from neurons recorded in the SNC and VTA, data from these regions were pooled.

PHARMACOLOGICAL MATERIALS. All drugs were added to the superfusate. Drug solutions required ~2 min for equilibration due to the time required by pass though a heat exchanger. Bicuculline-free base (BFB) was dissolved in dimethyl sulfoxide as a stock solution, which was subsequently diluted 1:1,000 in ACSF prior to use. We have shown previously that this concentration of dimethyl sulfoxide has no effect on burst firing (Johnson and Seutin 1997Go). Stock solutions of all other drugs were mad in water and diluted 1:1,000 in ACSF prior to use. Bicuculline methiodine (BMI), BFB, NMDA, and isoguvacine were obtained from Research Biochemicals International (Natick, MA), whereas apamin, dopamine HCl, and tetrodotoxin (TTX) were obtained from Sigma Chemical (St. Louis, MO).

BURST FIRING. Dopamine neurons were induced to fire in bursts by perfusing the slice with NMDA (20 µM) and apamin (100 nM) as described previously (Johnson et al. 1992Go). Burst firing was defined as a repeating pattern in which each "burst" of action potentials, consisting of at least three spikes firing at a frequency of >=10 Hz, is separated by an interburst hyperpolarization of >=5 mV (Seutin et al. 1993Go). A constant hyperpolarizing current (50-250 pA) was passed though the electrode as needed to prevent excessive NMDA-induced depolarization.

Modeling procedures

MODEL DEVELOPMENT. We used a modified version of compartmental model of DA neuron as described by Komendantov and Canavier (2002Go) but with a more realistic morphological structure, which consists conceptually of a soma and four identical branched dendrites with a single proximal and two distal branches (Fig. 1) and allowed us to use realistic values for the passive properties in intensive units, while still matching the macroscopic properties such as the input resistance and the time constant. The distal portions of the dendrites can sustain an oscillation in current even when the soma is voltage clamped (Johnson and North 1992Go), so the need for the distal dendritic compartment is clear. The inclusion of a proximal dendritic compartment allows the electrotonic effects of coupling between the soma and distal dendrites to be simulated more realistically with minimum computational overhead. It is certainly possible to reduce the model to two compartments, or even one, but in the case of a single compartment, direct comparison of the model with physiological data recorded from the soma would no longer be possible.



View larger version (19K):
[in this window]
[in a new window]
 
FIG. 1. Compartmental model. In this model, the soma has its own cylindrical compartment (S) with 4 symmetrical dendrites subdivided into a proximal dendrite (PD) and 2 distal dendrites (DD). The soma was modeled as a cylinder 15 µm in diameter and 25 µm long to which were connected the 4 dendrites. Each dendrite with branching structure was subdivided into 1 proximal compartment and 2 connected to it distal compartments. Four proximal dendrites and 8 distal dendrites were also modeled as cylinders 3 and 1.5 µm in diameter and 150 and 350 µm in length, respectively.

 
The revised model maintains a better correspondence than our previous models between its geometric and morphological parameters (Juraska et al. 1977Go; Preston et al. 1981Go; Tepper et al. 1987Go) and those of real DA neurons in vitro and in vivo. The electrical activity (Johnson et al. 1992Go; Paladini et al. 1999bGo; Wang et al. 1994Go) is also comparable to that of physiological DA neurons. The membrane potential of each compartment of the model neuron obeys the current balance equation



where INa,i is the fast sodium current, INaP,i is the sodium pump current, IA,i is the outward potassium current, IKDR,i is the potassium delay rectifier current, IL,i is the linear leak current, Isp, Ips, Ipd, and Idp are the intercompartmental coupling currents, INMDA,i is the NMDA-induced current, IGABAA,i is the GABAA-induced current, Cm is the specific membrane capacitance. The subscripts "d", "p", and "s" indicate distal dendritic compartment, proximal dendritic compartment and somatic compartment, respectively. The subscript "i" indicates the compartment. The change in morphological structure required that the passive parameters of the model and the parameters associated with INMDA be recalibrated, and a better fit was produced. In particular, the larger surface area afforded by including the secondary branching structure allowed the specific capacitance (Cm) in each compartment to be set to the more realistic value of 1 µF/cm2, while still matching the experimentally observed time constant and input resistance. The NMDA-gated currents were placed on the dendrites because morphological studies have shown that synaptic inputs that release glutamate are more frequently located on dendrites than on the somata of dopamine neurons (Sesack and Pickel 1992Go; Smith et al. 1996Go).

Calcium dynamics were added to the soma, which is described by a calcium balance equation

where fCa is the constant fraction of free calcium, and F is Faraday's constant. The calcium dynamics include a low-threshold calcium current ICa,T = gCa,TdTfT(Vs - ECa), an N-type calcium current ICa,N = gCa,NdN fCa,N(Vs - ECa), a nifedipine-sensitive calcium current ICa,L = gCa,LdL fCa,L(Vs - ECa), a calcium pump , a calcium leakage current IL,Ca,s = gL,Ca,s(Vs - ECa), and an apamin-sensitive current IK,Ca,SK = gK,Ca,SK {1/[1+ (KM,K,Ca/[Ca2+]in)4]}(Vs - EK). For all gating variables X = dT, dN, dL, fT, the dynamics of X was described by dX/dt = (X{infty}(Vs) - X)/{tau}x, where {tau}x is the time constant. The steady-state voltage dependence was determined using X{infty}(Vs) = 1.0/{1.0 + exp[-(Vs - Vhalf,x)/{theta}x]}, where Vhalf,x is the half activation voltage for the gating variable X, {theta}x is the slope factor for that variable. The time constants were described by {tau}x = Ax x exp(-[(Vs - Vhalf,{tau}x)/{sigma}x]2) + Bx. For Z = fCa,L, fCa,N, the steady-state calcium dependent inactivation function was modeled as Z{infty}(Ca2+) = KM,Z/(KM,Z + [Ca2+]in). Although voltage-activated calcium channels are present on the dendrites in these neurons (Wilson and Callaway 2000Go), we omitted calcium dynamics in the dendrites as the simplest way of implementing the segregation of calcium entry through NMDA receptors from SK channels that is observed experimentally (Canavier and Johnson 1999Go; Morikawa et al. 2003Go). The descriptions of the calcium currents were taken from Amini et al. (1999Go) with a few modifications as noted in the following text. Several currents included in Amini et al. (1999Go) were omitted for simplicity: a high-voltage threshold calcium current, the hyperpolarization-activated cation current, and the sodium-calcium exchanger. In addition, the slow component of ICa,T described by Amini et al. (1999Go) was ignored, the form of the calcium-dependency of IK,Ca,SK was changed to accommodate the data presented in this paper, and some parameters of IK,DR and ICa,L were adjusted. The NMDA-induced current (INMDA) was simulated using a Goldman-Hodgkin-Katz equation (Canavier 1999Go; Komendantov and Canavier 2002Go) and an instantaneous magnesium block (p{infty}) described by p{infty} = 0.0225 + 0.9775/[1 + ([Mg]out/Km,Mg) exp(-V/q)].

Currents evoked by GABAA receptors activation were modeled as a linear leak IGABAA,i = gGABAA,i(Vi - ECl). The simulated distribution of the GABAA receptors in different compartments (soma, proximal dendrites, and distal dendrites) was based on the assumption that synapses on cell body are mainly inhibitory (Kandel and Siegelbaum 2000Go) and evidence of a preferential effect of GABA on the somata of dopamine neurons compared with the dendrites (Kalivas 1993Go). We used the following ratio between the conductances for GABAA-induced currents in all simulations: gGABAA,s:gGABAA,p:gGABAA,d = 1:0.1:0.1. The leakage current (IL) composed of calcium (IL,Ca in soma only), sodium (IL,Na), and potassium (IL,K) was adjusted along with other passive parameters to obtain better correspondence to the input resistance (215 M{Omega}) (Johnson and North 1992Go; Ping and Shepard 1996Go), and time constant ({tau} = 31 ms) (Johnson and North 1992Go) experimentally observable in vitro. These values were measured by applying 94 pA of hyperpolarizing current via a simulated electrode in the soma to hold the membrane potential at -66 mV to suppress the spontaneous SOP, then applying an additional 71 pA of hyperpolarizing current for 600 ms to measure the input resistance and the time constant. The dynamics of the model neuron were described by 14 first-order differential equations in the soma (V, [Na]in, [Ca]in, m, h, n, q, s, dT, fT, dN, dL, fCa,L, fCa,N) and by eight first-order differential equations each in the proximal and distal dendrites (V, [Na]in, m, h, n, p, q, and s). As in Canavier (1999Go), symmetry was invoked so that only a single instance of the somatic, proximal dendritic, and distal dendritic compartments was integrated. The additional currents from Amini et al. (1999Go) were inserted into the soma only of the single neuron model described in the appendix of Komendantov and Canavier (2002Go). The values of the parameters that are not the same as those given in Komendantov and Canavier (2002Go) or Amini et al. (1999Go) are presented in Table 1.


View this table:
[in this window]
[in a new window]
 
TABLE 1. Model parameters

 
CALIBRATION OF THE APAMIN-SENSITIVE AND NMDA CURRENTS IN THE MODEL. The apamin-sensitive current was modeled with a Michaelis-Menten dependence on cytosolic Ca2+ concentration with a half-activation at 190 nM and a Hill coefficient of 4. The somatic conductance of this current (900 µS/cm2) was set to produce agreement with experimental values. A simulated voltage-clamp experiment that corresponds to experimental data (Canavier and Johnson 1999Go) with a 350-ms voltage step to -30 mV from a holding potential -60 mV is shown in Fig. 2, A1. Simulated block of the SK channels eliminates the peak current of ~1,200 pA with a latency of ~200 ms, and the peak tail current of 450 pA.



View larger version (12K):
[in this window]
[in a new window]
 
FIG. 2. Model calibration. A: simulated voltage-clamp data. A1: simulated current recorded before, during, and after a simulated 350-ms step to -30 mV from holding potential of -60 mV in the absence (—) and simulated blockage of IK,Ca,SK (- · -). · · ·, 0 current. A2: comparison of simulated N-methyl-D-aspartate (NMDA) difference currents recorded at the end of a simulated 800-ms step to the test potential from a holding potential of -60 mV. - · -, experimental data (Wu and Johnson 1996Go); —, simulated data. All model concentrations were equilibrated at -60 mV before each step. B: simulation of pacemaker activity in the absence of NMDA and GABAA receptor activation (PNMDA = 0; gGABAA,i = 0). - - -, 0 mV. B1: spiking activity in control. B2: spiking activity under simulated SK channel block (gK,Ca,SK = 0). B3: slow oscillations of membrane potential under simulated TTX block of spike generation (gNA,i = 0).

 
The parameters of model NMDA-induced current were fit to the experimental data of Wu and Johnson (1996Go) using a simulated voltage clamp in the soma. The conductances for the L-type calcium current (gCa,L) and the fast sodium current (gNa,i) were set to zero to simulate the application of nifedipine and tetrodotoxin, respectively. The holding potential was -60 mV, the pulse duration was 800 ms, and the incremental voltage step was 10 mV. The control current was recorded at the end of the test potential. The series of experiments in the presence of 10 µM NMDA was reproduced with PNMDA = 0.4 · 10-6 cm/s, and the recorded current was subtracted from the control current to obtain the difference current. The ratio of PNMDA,d: PNMDA,p = 1:1 on the distal and proximal dendrites produced the best fit to the experimental data. Figure 2A2 shows that the simulated NMDA difference currents (—) closely resemble those (- · -) obtained by Wu and Johnson (1996Go) in 10 µM NMDA. The region of negative slope in the current-voltage curve from -90 to -50 mV causes positive feedback, which is essential to the generation of the oscillation underlying burst firing (Mereu et al. 1991Go; Overton and Clark 1997Go; Wu and Johnson 1996Go).

The model produces spontaneous spiking in a simulation of an in vitro preparation in normal Ringer solution (PNMDA,i = 0; gGABAA,i = 0; Fig. 2B1). The simulated application of apamin by setting gK,Ca,SK = 0 (Fig. 2B2) decreased the afterhyperpolarization (AHP) after the action potential by 15 mV and increased spike frequency, which is consistent with experimental data (Johnson and Seutin 1997Go; Ping and Shepard 1996Go; Shepard and Bunney 1991). If application of TTX is simulated instead by setting gNA,i = 0, slow oscillations in membrane potential with an amplitude of 10 mV (Fig. 2B3) result from the interaction of Ca2+ and Ca2+-dependent conductances (Kang and Kitai 1993Go).

MODEL IMPLEMENTATION. As in previous research (Canavier 1999Go; Komendantov and Canavier 2002Go), the simulations and programs for analysis were coded in the C programming language. Numerical integrations of simulations were performed using the Fortran implementation of an implicit Runge-Kutta method of order five with variable step size (Hairer and Wanner 1996Go).

ANALYSIS OF MODEL ACTIVITY. The analysis of model responses used the time courses of V and interspike intervals (ISIs); the first 30 s of a simulation was considered as a transient period and excluded. The next 20 s of simulations were recorded and analyzed. The control parameters were the permeability for the NMDA-induced current in the dendrites (PNMDA), corresponding to the level of excitation, and the conductance for the GABAA receptor current in soma (gGABAA,s), corresponding to the level of inhibition. The automated determination of the mode of activity (regular and irregular spiking, regular and irregular bursting) was made according to algorithms used for the previous model (Komendantov and Canavier 2002Go) with minor changes to the evaluation criteria for burst determination. Bursts were defined as groups of at least two spikes terminated by an interburst interval of >=160 ms (Grace and Bunney 1984bGo) with hyperpolarization of >=5 mV (Johnson and Seutin 1997Go). The selection criteria for detection of irregularity in the model activity were determined empirically and confirmed by visual observation of the membrane potential waveform in many instances.

The dependence of input resistance on gGABAA was calculated with the NMDA-induced permeability (PNMDA) set to zero using the methods described in the preceding text. The values of the conductance for the GABAA receptor currents were varied from 0 to 2,750 µS/cm2 for the soma and from 0 to 275 µS/cm2 for the dendritic compartments. We selected the range of gGABAA that produced a physiologically acceptable range of input resistance in vivo.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 REFERENCES
 
Effect of a GABAA agonist on the firing pattern in in vitro experiments

We used a concentration of 100 µM isoguvacine because preliminary studies showed that this concentration of this GABAA receptor agonist is near the EC50 value for increasing conductance in DA neurons. In all dopamine neurons tested (n = 3), isoguvacine markedly inhibited NMDA-dependent burst firing and membrane oscillations, as seen in Fig. 3A.



View larger version (27K):
[in this window]
[in a new window]
 
FIG. 3. Effects of GABAA receptor stimulation on NMDA-induced burst firing in vitro. A: experimental activation of GABAA receptors interferes with burst firing and the membrane oscillation induced by NMDA (20 µM) and apamin (100 nM). Isoguvacine (100 µM) converts burst firing to single-spike firing in a ventral tegmental area (VTA) dopamine neuron. The inhibitory effect of isoguvacine is prevented by pretreatment with the GABAA receptor antagonist bicuculline-free base (BFB; 30 µM). Voltage recordings were made with a pipette that contained potassium acetate; a constant hyperpolarizing current (180 pA) was passed through the pipette. - - -, -50 mV. A1: control. A2: isoguvacine. A3: recovery in BFB. A4: isoguvacine in BFB. B: computer simulations. The simulations were done with injection of 180 pA of constant hyperpolarizing current and under PNMDA = 0.75 · 10-6 cm/s and gK,Ca,SK = 0 µS/cm2, that corresponds to application of NMDA and apamin. - - -, 0 mV. B1: burst firing, gGABAA,s = 0. B2: simulated application of isoguvacine converts burst firing to single-spike firing, gGABAA,s = 350 µS/cm 2.

 
When recording with electrodes filled with potassium acetate, isoguvacine converted burst firing to a single-spike pattern (n = 3) with little or no change in average membrane potential (see Fig. 3A). In contrast, this GABAA receptor agonist caused membrane depolarization (10-20 mV) and blocked all action potentials when recordings were made with microelectrodes filled with potassium chloride (n = 4). In the absence of NMDA and apamin, isoguvacine reduced the cellular input resistance by 58 ± 6% (n = 3). Isoguvacine acted within 3 min of starting superfusion. Recovery was observed 10-15 min after washout of isoguvacine.

Simulation of effects of GABAA receptor stimulation on bursting in vitro

Figure 3B shows that the simulated addition of isoguvacine in the bath in the simulated presence of both NMDA and apamin converts burst firing to single-spike firing as in the in vitro experiments (Fig. 3A, 1 and 2). The application of isoguvacine (100 µM) was simulated by adding an ohmic GABAA receptor current that reversed at -70 mV with gGABAA,s = 350 µS/cm2, gGABAA,p = 35 µS/cm2, gGABAA,d = 35 µS/cm2. The value of the reversal potential for this current was chosen based on its correspondence to the chloride reversal potential determined experimentally in DA neuron using micropipettes filled with potassium acetate (Johnson and North 1992Go). GABAA receptor-mediated depolarizations in dopamine neurons have been shown to be due to the shift in the Nernst potential for chloride, which is -72 ± 4 mV with potassium acetate-filled microelectrodes but shifts to -36 ± 6 mV when potassium chloride-filled microelectrodes raise the intracellular concentration of chloride (Johnson and North 1992Go).

The simulated application of isoguvacine in the absence of NMDA and apamin decreases the input resistance by 63% (from its control value of 215 to 80 M{Omega}), which is within the range of the decrease in input resistance (58 ± 6%) that we observed experimentally. In both experiment and simulation, a constant hyperpolarizing current of 180 pA was applied, and the average membrane potential did not change significantly.

Joint modulation of the firing pattern by the activity of NMDA and GABAA receptors

The firing pattern of DA neurons depends on the activation of both excitatory inputs and inhibitory inputs. Komendantov and Canavier (2002Go) investigated the dependence of the firing pattern of an in vitro model of DA neuron on the simulated level of activation of the excitatory NMDA inputs. As PNMDA was increased, the firing pattern changed from low-frequency single spiking to bursting and then to high-frequency spiking. We observed similar transitions here. An increase in PNMDA leads to an increase in the firing frequency as well as transitions from low-frequency spiking through regular and irregular bursting to high-frequency spiking. These sequences correspond to those observed in in vitro experiments as the concentration of NMDA in the bath is increased (Johnson et al. 1992Go; Paladini et al. 1999bGo; Wang et al. 1994Go; Wu and Johnson 2001Go) except that those studies show depolarization block rather than high-frequency spiking.

The conductance for IGABAA,s was also varied at a fixed value of PNMDA = 2.3 · 10-6 cm/s, which results in burst firing in the absence of GABAA receptor activation. At this fixed value of PNMDA, the different regimes may be distinguished as gGABAA,s is increased: regular bursting (Fig. 4A, 1 and 3), irregular bursting (Fig. 4A2), single spiking (Fig. 4A, 4 and 5) and a silent regime (gGABAA,s > 2,500 µS/cm2). The amplitude of the slow wave decreased as gGABAA,s was increased (Fig. 4B). On the diagram, the minima corresponding to the AHPs after action potentials were extracted using the first and second time derivatives of the membrane potentials (dVs/dt = 0, d2Vs/dt2 > 0). Of these minima, only the most and least hyperpolarized minima at each value of gGABAA,s were plotted on Fig. 4B, thus the upper branch represents the shallowest hyperpolarization during any burst and the lower branch represents the deepest hyperpolarization during any interburst interval.



View larger version (23K):
[in this window]
[in a new window]
 
FIG. 4. Dynamic responses of a model neuron as a function of the simulated tonic level of GABAA receptor activation (gGABAA,s) at a constant level of NMDA receptor activation (PNMDA = 2.3 · 10-6 cm/s). A: simulated firing pattern at various levels of GABAA receptor activation. - - -, 0 mV. A1: bursting at gGABAA,s = 500 µS/cm2. A2: bursting at gGABAA,s = 950 µS/cm2. A3: bursting at gGABAA,s = 1,100 µS/cm2. A4: regular spike firing at gGABAA,s = 1,800 µS/cm2. A5: low-frequency spike firing at gGABAA,s = 2,500 µS/cm2. B: bifurcation diagram for membrane afterhyperpolarizations (AHP). According to the quantitative criteria for analysis of firing activity (see METHODS), the region below 1,600 µS/cm2 corresponds to bursting, between 1,600 and 1,780 µS/cm2 to irregular spiking, between 1,700 and 2,500 µS/cm2 to regular single spiking, and above to silent hyperpolarization.

 
We investigated numerically the dependence of the firing pattern on the activation levels of NMDA and GABAA receptors simultaneously using a two-parameter bifurcation diagram (Fig. 5A). Each square on the diagram represents the type of activity that can be obtained at given values of gGABAA,s and PNMDA. The following modes were observed: silence (hyperpolarized, nonfiring), regular spiking (RS), irregular spiking (IS), regular bursting (RB), and irregular bursting (IB). Similar modes are also observed experimentally in vivo (Dai and Tepper 1998Go; Grace and Bunney 1984aGo,bGo).



View larger version (41K):
[in this window]
[in a new window]
 
FIG. 5. Simulated effect of simultaneously varying the activation of GABAA and NMDA receptors on the dynamics of model DA neuron. Each diagram (A, C, and D) is the result of 2,016 integrations of model equations with different values of parameters gGABAA,i and PNMDA (56 and 36) and detection of different modes. Incremental steps for gGABAA,s and PNMDA are 50 µS/cm2 and 0.1 · 10-6 cm/s, respectively. Blue bars near the scales show ranges of the values (input resistance, average and maximal frequency, number spikes per burst) taken from in vivo experimental data (see text). A: gGABAA,s - PNMDA phase diagram of model DA neuron. Model investigation reveals the following regimes: silent (S), light gray squares; irregular spiking (IS), light blue squares; regular spiking (RS), blue squares; irregular bursting (IB), yellow squares; regular bursting (RB), red squares. At very high values of PNMDA, depolarization block is achieved (not shown). B: input resistance of model neuron as a function of gGABAA,s. Green bar indicates the range of gGABAA,s that corresponds to the range of input resistance measured in vivo (blue bar). C: effect of gGABAA,s and PNMDA on average frequency of the model neuron. Here and in D the range of shaded values is indicated by the associated column of color bars on the right. Thick black lines border the region of bursting activity. D: effect of gGABAA,s and PNMDA on number of spikes per burst in the model neuron. Gray areas correspond to nonbursting modes of activity.

 
Along the y axis where gGABAA,s is equal to zero, the firing pattern corresponds to the in vitro preparation in the absence of any tonic activation of GABAA receptors.

Calibration of the GABAergic inputs

Grace and Bunney (1983Go) showed that nigral dopaminergic neurons have an input resistance in vivo in the range of 18-45 M{Omega}. Figure 5B shows the dependence of input resistance of model DA neuron on the value of conductance of the linear GABAA receptor current, gGABAA,s. The range of gGABAA,s from 850 to 2,750 µS/cm2 (green bar) corresponds to the range of input resistance experimentally observable in vivo (blue bar). Thus we can use Fig. 5A to explain the effects of the application of GABAA antagonists in vivo (Paladini and Tepper 1999Go; Paladini et al. 1999aGo). Every single instance of single-spike firing in the in vivo range (to the right of the dashed line) can be converted to burst firing when gGABAA,s is decreased by 41-65% of its original in vivo value, simulating the in vivo application of a GABAA channel blocker such as picrotoxin. This is consistent with the results of Paladini and Tepper (1999Go) that revealed that application of GABAA antagonists in vivo (picrotoxin, 200-1,000 µM, and bicuculline methiodide, 200-400 µM) "caused a robust change to a bursty pattern regardless of the baseline firing pattern," although the exact fraction of GABAA receptors blocked and the change in input resistance after the application of picrotoxin the in vivo experiments is not known precisely.

Quantitative analysis of model activity

We computed the average (Fig. 5C) spike frequency for each pair of values of the control parameters, gGABAA,s and PNMDA, for which the model demonstrated spiking or bursting spontaneous activity. Also, we counted the number spikes per burst (Fig. 5D) for each pair of the control parameters when the activity was determined as bursting. In the simulated physiological range of gGABAA,s from 850 to 2,750 µS/cm2, average spiking frequencies are always less 6.5 Hz for low-frequency single-spike firing, whereas they range from 5.0 to 16 Hz for burst firing. Thus there is some correlation between firing pattern and average frequency (Fig. 5, A and C). On each diagram C and D, the vertical dashed line bounds the area which corresponds to the putative in vivo of gGABAA,s. In in vivo experiments, the average frequency has been reported as 0-7.1 Hz (Richards et al. 1997Go), 0 to 7.0-7.2 Hz (Grace and Bunney 1984bGo), and <10 Hz (Dai and Tepper 1998Go); the number spikes per burst lies in the range from 2 to 10, and the maximal number spikes per burst is 23 (Grace and Bunney 1984).

Figures 6, A1 and A2, shows two possible effects of blocking the IK,Ca,SK current on the in vivo activity of DA neuron: a transition from single spiking (A1) to bursting (A2) and a transition from short bursts that contain few spikes to longer bursts containing more spikes (A, 3 and 4).



View larger version (40K):
[in this window]
[in a new window]
 
FIG. 6. Changes in the firing pattern resulting from the simulated block of IK,Ca,SK. A: time courses of membrane potentials at different levels of PNMDA and gGABAA,s. · · ·, 0 mV. A, 1 and 2: PNMDA = 1.5 · 10-6 cm/s, gGABAA,s = 1,000 µS/cm2. A1: gK,Ca,SK = 900 µS/cm2, average frequency is 4.93 Hz. A2: gK,Ca,SK = 0, average frequency is 6.36 Hz. A, 3 and 4: PNMDA = 2.5 · 10-6 cm/s, gGABAA,s = 1,450 µS/cm2. A3: gK,Ca,SK = 900 µS/cm2, average frequency is 8.54 Hz. A4: gK,Ca,SK = 0, average frequency is 11.41 Hz. B: effect of SK channel blockers on the phase diagram. The letters mark zones with different regimes of activity according to the abbreviations presented in the legend of Fig. 5A. B1: the phase diagram in gGABAA,s and PNMDA coordinates corresponding to the in vivo part of the diagram in Fig. 5A. The diagram is the result of 1,040 integrations of model equations with different pairs of values of parameters gGABAA,s and PNMDA. B2: the same phase diagram obtained with IK,Ca,SK set to 0 to simulate blocking the SK current: about a quarter of the single spiking neurons observed in B1 are converted to bursters by this manipulation.

 
Figure 6B1 reproduces Fig. 5A but only for gGABAA,s from 850 to 2,750 µS/cm2 and PNMDA from 1.0 · 10-6 to 3.5 · 10-6cm/s. Figure 6B2 displays the transformation of this phase diagram when gK,Ca,SK is set to zero. This simulated application of an SK channel blocker converts 59 of the 270 instances of single-spike firing shown in B1, or ~22%, to burst firing and prolongs burst duration of the instances of burst firing. The number of instances of irregular spiking in the in vivo range of parameters increased by 60%, from 45 to 72.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 REFERENCES
 
Relevance to behavior

Dopamine neurotransmission probably operates on at least three time scales (Schultz 1998Go): reward responses on the order of tens and hundreds of milliseconds; the processing of a wide array of rewards, feeding, drinking, punishments, stress, and social behavior of a time scale of seconds to minutes; and an enabling function of a variety of motor, cognitive, and motivational processes on a slow, and possibly tonic, time scale. Phasic firing in DA neurons increases within 100 ms after the onset of a salient sensory event (Horvitz 2002Go). Our results suggest that, due to the relatively slow rate of change of sodium concentration, a DA neuron must already be biased near the bursting range in order for a transient input to cause bursting in such a short interval. The changes in firing pattern we present here are representative of persistent activity after all transients dissipate, under conditions of inputs that remain nearly constant for seconds to minutes.

Role of different receptors in the activation of glutamatergic inputs

The main sources of EAA inputs (for a review, see Kitai et al. 1999Go; Overton and Clark 1997Go) to the midbrain DA neurons are the prefrontal cortex (PFC) (Sesack and Pickel 1992Go), the subthalamic nucleus (STN) (Kita and Kitai 1987Go), and the pedunculopontine nucleus (PPN) (Tokuno et al. 1988Go). Minor EAA projections include the amygdala (McDonald 1996Go; Phillipson 1979Go), the laterodorsal tegmental nucleus (Clements and Grant 1990Go; Gould et al. 1989Go), and the habenula nucleus (Matsuda and Fujimara 1992Go). Recently a novel EAA projection from the bed nucleus of stria terminalis has been identified (Georges and Aston-Jones 2002Go). Although glutamate also activates non-NMDA ionotropic receptors (AMPA and kainate) and metabotropic receptors, we have chosen to focus initially on the NMDA receptor, based on the following evidence. The glutamate-induced excitation of DA neurons is more sensitive to the selective blockade of NMDA versus non-NMDA receptors because iontophoresis of glutamate and the selective NMDA receptor antagonist (±)-3-(-2-carboxypiperazin-4-yl)-propyl-1-phosphonic acid (CPP) produced an 80% reduction of the glutamate-induced increase in firing (Christoffersen and Meltzer 1995Go), whereas glutamate and the AMPA receptor antagonist 2,3-dihydroxy-6-nitro-7-sulfamoyl-benzo-(f)quinoxaline (NBQX) produced only a 26% inhibition of the glutamate-induced excitation. Similarly, Chergui et al. (1993Go) showed that bursting was blocked by the local application (iontophoresis and pressure ejection) of the NMDA receptor antagonist 2-amino-5-phosphonovaleric acid (AP-5) but not by the AMPA/kainate receptor antagonist 6-cyano-7-nitroquinoxalene-2,3-dione (CNQX) (see also Overton and Clark 1997Go). The duration of the excitatory postsynaptic current (EPSC) evoked by AMPA is ~20 ms (Bonci and Malenka 1999Go), whereas that evoked by NMDA is much longer lasting (Mercuri et al. 1996Go; Overton and Clark 1997Go), on the order of 600 ms (Destexhe et al. 1998Go). Glutamate affinity for AMPA receptors is lower (EC50 = 19 µM) than for NMDA receptors (EC50 = 2.3 µM) (Patneau and Mayer 1990Go), thus lower concentrations of glutamate in the cleft are required to activate NMDA receptors, which accounts at least in part for the longer duration of the activation of these receptors. Recent simulations (Canavier and Landry 2002Go) suggest that only EAA afferent firing frequencies at which NMDA but not AMPA inputs summate temporally are effective in evoking burst firing. Thus modeling the activation of NMDA receptors only is a reasonable first approximation of the effects of glutamatergic inputs.

In the model, we observed high-frequency bursts with a large number of spikes per burst at high values of PNMDA that have not been observed experimentally in vivo. We hypothesize that in vivo the presence of AMPA receptors prevents this pattern of activity. Similar activity with a large number of spikes per burst was obtained in vitro using bath application of NMDA (Wu and Johnson 2001Go). The inclusion of AMPA receptor activation in the model as a sum of linear leaks for Na+ and K+ (see Canavier 1999Go) converts this high-frequency spiking into depolarization block (for example, for PNMDA = 3.4 · 10-6 cm/s, gGABAA,i = 0 µS/cm2, gAMPA,K = 30 µS/cm2, gAMPA,Na = 60 µS/cm2). We conclude that AMPA receptors play a more significant role at high levels of glutamate due to the lower affinity of AMPA receptors (Patneau and Mayer 1990Go) and suggest that AMPA receptors can cause depolarization block at high levels of glutamate rather than the very high-frequency spiking or bursting (~28 Hz) observed in the model.

The model predicts some correlation (Fig. 5, A and C) between firing pattern and firing frequency; however, the presence of such a correlation is controversial (see Freeman et al. 1985Go; Grace and Bunney 1984bGo; Paladini and Tepper 1999Go). Paladini and Tepper (1999Go) did not observe any correlation in vivo when the application of bicucculine or picrotoxin converted single-spike firing to burst firing. Our model, with its assumption of constant parameters that correspond to activation of NMDA and GABAA receptors by the average concentration of transmitter in the cleft, cannot account for this lack of correlation. Therefore additional mechanisms must be postulated to account for it, and we hypothesize that the transient nature of the synaptic dynamics observed in vivo may provide these additional mechanisms.

Possible in vivo analogues of blockade of apamin-sensitive current

Although no endogenous neuromodulator is known to act directly on the receptor on the SK channel that binds apamin, the modulation of this current is nonetheless relevant to the modulation of the firing pattern of dopamine neurons by synaptic afferents. It is not necessary for a neuromodulator to act directly on this current to block it because reducing its access to calcium activation will effectively attenuate it. For example, muscarine (50 µM) reduces the amplitude of the slow AHP that follow action potentials in most substantia nigra DA neurons (Scroggs et al. 2001Go). The SK current is responsible for the AHP (Shepard and Bunney 1991), and muscarine is believed to exert its effects by reducing calcium entry. Another possible modulatory pathway involves inhibitory PSPs evoked by activation of metabotropic glutamate receptors (mGluR). This results in production of inositol triphosphate, IP3, which causes release of calcium from intracellular stores and, in turn, the activation of the SK channels (Fiorillo and Williams 1998Go). There is evidence that activation of {alpha}1 adrenergic receptors or M1 muscarinic receptors interferes with the release of calcium from intracellar stores (Fiorillo and Williams 2000Go; Paladini et al. 2001Go), decreasing the amount of SK current evoked by mGluR activation. Thus both the noradrenergic afferents from the locus cerulus and the cholinergic afferents from the PPN may attenuate the apamin-sensitive SK current, and there have been suggestions (Brodie et al. 1999Go) that serotonin, possibly acting via IP3-coupled 5HT2 receptors, may attenuate the SK current as well, which would tend to promote burst firing. Also, dopamine, which is released somatodendritically by DA neurons, has some affinity for the {alpha}1 receptor. Thus the results of our simulations showing that blockage of Ca2+-sensitive potassium current may promote or facilitate bursting in DA neuron in vivo (Fig. 6) are probably physiologically relevant.

Recently, Wolfart and Roeper (2002Go) showed that in response to 20-ms hybrid-clamp depolarizations that evoked action potentials, only calcium entering via the T-type calcium channel activates the apamin-sensitive current. Thus neuromodulators that reduce the T-type calcium current would also indirectly reduce the SK current under similar circumstances. It is possible that the experimental conditions did not mimic the activation pattern of calcium channels exhibited during spontaneous, pacemaker-like firing, which has a contribution from additional calcium currents, such as the L-type calcium current (Kang and Kitai 1993Go; Mercuri et al. 1994Go; Nedergaard et al. 1993Go; Takada et al. 2001Go). On the other hand, Wolfart and Roeper have suggested that SK channels, T-type calcium channels, and intracellular calcium stores might be selectively colocalized, forming a specialized calcium-signaling complex in which T-type calcium channel influx triggers secondary calcium release that specifically activates SK channels, whereas other types of calcium influx would not. To address these issues, more detailed models that incorporate compartmentalization of calcium are required. Adding voltage-gated calcium channels and calcium dynamics to the dendritic compartments may enhance the ability of the SK channel current to suppress burst firing, as well as the ability of the model to exhibit irregular firing (Wilson and Callaway 2000Go) due to the desynchronization of calcium oscillations in the various compartments.

How the various manipulations induce or suppress burst firing

This is the first model in which we have combined sodium and calcium dynamics and can simulate regular, calcium-dependent spike firing as well as NMDA-induced burst firing and, as such, is our most complete model of the regulation of the firing pattern in dopamine neurons to date. This model is more focused on the interaction of the various synaptic and intrinsic currents in producing single-spike firing versus bursting than other models of dopamine neurons (Kotter and Feizelmeier 1998Go; Li et al. 1996Go; Penney and Britton 2002Go; Wilson and Callaway 2000Go). In previous papers (Canavier 1999Go; Komendantov and Canavier 2002Go), we used nullcline analyses to gain insight into the mechanisms underlying NMDA-induced burst firing. Although it is not practical to use nullcline analysis on the complex model in this paper, many of the principles from our earlier work still apply. A regenerative, nonlinear current such as the NMDA-induced current is required for a bursting oscillation, whereas linear GABAA receptor currents oppose any oscillation. Therefore the many of the modeling results presented herein are likely to be robust and not highly dependent on the precise mechanism of burst repolarization.

It has been suggested (Paladini et al. 1999bGo) that GABAA receptor activation prevents the observation of bursting at the soma by a shunt mechanism. If the GABAA receptor current was merely acting as a shunt, that would imply that the distal dendrites continue to oscillate as they did before the GABAA agonist was applied but that the current leaks out of the GABAA channels and so does not reach the soma. In the model, the application of the excess linear current actually suppresses the oscillation in the dendrites, so there is no burst to spread passively (or actively) to the soma. By analogy to classical mechanics, we call this alternative mechanism linear damping. There are other considerations, such as bias. A higher level of NMDA receptor activation does not always lead to more bursting. Excessive depolarization causes the neuron to leave the voltage range in which the nonlinearity is active due to positive feedback is active because all available NMDA channels are already open. Sometimes hyperpolarization that is not excessive can facilitate bursting by bringing the neuron back into the range in which the positive feedback promotes an oscillation. Finally, the spikes themselves may promote or abolish burst firing via complex nonlinear mechanisms. Therefore increases in gGABAA evoke complex nonlinear activity (Fig. 4) that eventually leads to the suppression of bursting. Whereas in general increasing glutamatergic input may promote burst firing and in general increasing GABAergic input reduces burst firing, Fig. 5 shows that the ability of a change in one input to effect a change in pattern often depends on the value of the other input. An interesting possibility is that somatic hyperpolarization promotes or actually enables burst firing in the presence of dendritic glutamatergic input by countering excessive depolarization that prevents the dendrites from bursting.

The SK channel current is known to contribute to the regularity of single-spike firing in dopamine neurons because it mediates the repolarization of the slow oscillatory potential that underlies pacemaker-like firing, thus blocking this current disrupts regular pacemaker-like firing. The nonlinearity associated with the apamin-sensitive current is not regenerative and therefore works to suppress any burst generating nonlinearity (Amini et al. 1999Go).

The blockade by picrotoxin of GABAA receptors is maximally effective in inducing bursting. Bicuculline blocks not only GABAA receptors, but also the SK channels, the blockage of which also promotes burst firing both in vitro (Seutin et al. 1993Go) and in vivo (Seutin et al. 2002Go), so it is not surprising that bicuculline is also maximally effective in inducing burst firing in vivo (Paladini and Tepper 1999Go). A salient prediction of the model is that blocking the SK channel current in vivo will facilitate bursting but not as robustly as blocking GABAA receptors. If the model is shown to underestimate this effect, the contribution of the SK channel on the dendrites must be reconsidered.


    DISCLOSURES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 REFERENCES
 
This work was supported by National Institutes of Health Grants NS-37963 (to C. C. Canavier) and MH-40416 (to S. W. Johnson). A. Komendantov received support from U. S. Department of Energy Grant DE-FG02-01ER63119.


    FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: C. C. Canavier, Dept. of Psychology, University of New Orleans, New Orleans, LA 70148 (E-mail: ccanavie{at}uno.edu).


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 REFERENCES
 
Amini B, Clark JW, and Canavier CC. Calcium dynamics underlying pacemaker-like and burst firing in midbrain dopaminergic neurons: a computational study. J Neurophysiol 83: 2249-2261, 1999.

Bolam JP and Smith Y. The GABA and substance P input to dopaminergic neurones in the substantia nigra of the rat. Brain Res 529: 57-78, 1990.[CrossRef][Web of Science][Medline]

Bonci A and Malenka RC. Properties and plasticity of excitatory synapses on dopaminergic and GABAergic cells in the ventral tegmental area. J Neurosci 19: 3723-3730, 1999.[Abstract/Free Full Text]

Brodie MS, McElvain MA, Bunney EB, and Appel SB. Pharmacological reduction of small conductance calcium-activated potassium current (SK) potentiates the excitatory effect of ethanol on ventral tegmental area dopamine neurons. J Pharmacol Exp Ther 290: 325-333, 1999.[Abstract/Free Full Text]

Canavier CC. Sodium dynamics underlying burst firing and putative mechanisms for the regulation of the firing pattern in midbrain dopamine neurons: a computational approach. J Comput Neurosci 6: 49-69, 1999.[CrossRef][Web of Science][Medline]

Canavier CC and Johnson SW. Computer simulation suggest a mechanism by which apamin promotes NMDA-induced burst firing in rat midbrain dopamine neurons. Soc Neurosci Abstr 25: 765.5, 1999.

Canavier CC and Landry RS. Firing pattern of model dopamine neuron depends upon the frequency of unpatterned excitatory amino acid inputs (CD-ROM). Soc Neurosci Abstr 764.14, 2002.

Celada P, Paladini CA, and Tepper JM. GABAergic control of rat substantia nigra dopaminergic neurons: role of globus pallidus and substantia nigra pars reticulata. Neuroscience 89: 813-825, 1999.[CrossRef][Web of Science][Medline]

Chergui K, Charlety PJ, Akaoka H, Saunier CF, Brunet J-L, Buda M, Svensson TH, and Chouvet G. Tonic activation of NMDA receptors causes spontaneous burst discharge of rat midbrain dopamine neurons in vivo. Eur J Neurosci 5: 137-144, 1993.[CrossRef][Web of Science][Medline]

Chergui K, Nomikos GG, Mathe JM, Gonon F, and Svensson TH. Burst stimulation of the medial forebrain bundle selectively increases Fos-like immunoreactivity in the limbic forebrain of the rat. Neuroscience 72: 141-156, 1996.[CrossRef][Web of Science][Medline]

Chergui K, Svenningsson P, Nomikos GG, Gonon F, Fredholm BB, and Svensson TH. Increased expression of NGFI-A mRNA in the rat striatum following burst stimulation of the medial forebrain bundle. Eur J Neurosci 9: 2370-2382, 1997.[CrossRef][Web of Science][Medline]

Christoffersen CL and Meltzer LT. Evidence for N-methyl-D-aspartate and AMPA subtypes of the glutamate receptor on substantia nigra dopamine neurons: possible preferential role for N-methyl-D-aspartate receptors. Neuroscience 67: 373-381, 1995.[CrossRef][Web of Science][Medline]

Clements JR and Grant S. Glutamate-like immunoreactivity in neurons of the laterodorsal tegmental and pedunculopontine nucleus in the rat. Neurosci Lett 120: 70-73, 1990.[CrossRef][Web of Science][Medline]

Dai M and Tepper JM. Do silent dopaminergic neurons exist in rat substantia nigra in vivo? Neuroscience 85: 1089-1099, 1998.[CrossRef][Web of Science][Medline]

Destexhe A, Mainen ZF, and Sejnowski TJ. Kinetic models of synaptic transmission. In: Methods in Neuronal Modeling—From Ions to Networks (2nd ed.), edited by Koch C and Segev I. Cambridge, MA: The MIT Press, 1998, p. 1-25.

Fiorillo CD and Williams JT. Glutamate mediates an inhibitory postsynaptic potential in dopamine neurons. Nature 394: 19-21, 1998.[CrossRef][Medline]

Fiorillo CD and Williams JT. Cholinergic inhibition of ventral midbrain dopamine neurons. J Neurosci 20: 7855-7860, 2000.[Abstract/Free Full Text]

Freeman AS, Meltzer LT, and Bunney BS. Firing properties of substantia nigra dopaminergic neurons in freely moving rats. Life Sci 36: 1983-1994, 1985.[CrossRef][Web of Science][Medline]

Georges F and Aston-Jones G. Activation of ventral tegmental area cells by the bed nucleus of the stria terminalis: a novel excitatory amino acid input to midbrain dopamine neurons. J Neurosci 22: 5173-5187, 2002.[Abstract/Free Full Text]

Gould E, Woolf N, and Butcher L. Cholinergic projections to the substantia nigra from the pedunculopontine and lateral dorsal tegmental nuclei. Neuroscience 28: 611-623, 1989.[CrossRef][Web of Science][Medline]

Grace AA and Bunney BS. Intracellular and extracellular electrophysiology of nigral dopaminergic neurons. I. Identification and characterization. Neuroscience 10: 301-315, 1983.[CrossRef][Web of Science][Medline]

Grace AA and Bunney BS. The control of firing pattern in nigral dopamine neurons: single spike firing. J Neurosci 11: 2866-2876, 1984a.

Grace AA and Bunney BS. The control of firing pattern in nigral dopamine neurons: burst firing. J Neurosci 11: 2877-2890, 1984b.

Grace AA and Onn SP. Morphology and electrophysiological properties of immunocytochemically identified rat dopamine neurons recorded in vitro. J Neurosci 14: 765-774, 1989.

Hairer E and Wanner E. Solving ordinary differential equations. II. Stiff and differential-algebraic problems. Springer Ser Comput Math 14: 118-130, 565-576, 1996.

Horvitz JC. Dopamine gating of glutamatergic sensorimotor and incentive motivational input signals to the striatum. Behav Brain Res 137: 65-74, 2002.[CrossRef][Web of Science][Medline]

Johnson SW and Canavier CC. GABAA receptor stimulation interferes with burst firing by increasing a linear conductance in rat midbrain dopamine neurons. Soc Neurosci Abstr 24: 1641, 1998.

Johnson SW and North RA. Two types of neurone in the rat ventral tegmental area and their synaptic inputs. J Physiol 450: 455-468, 1992.[Abstract/Free Full Text]

Johnson SW and Seutin V. Bicuculline methiodide potentiates NMDA-dependent burst firing in rat dopamine neurons by blocking apamin-sensitive Ca2+ -activated K+ currents. Neurosci Lett 231: 13-16, 1997.[CrossRef][Web of Science][Medline]

Johnson SW, Seutin V, and North RA. Burst firing in dopamine neurons induces by N-methyl-D-aspartate: role of electrogenic sodium pump. Science 258: 665-667, 1992.[Abstract/Free Full Text]

Juraska JM, Wilson CJ, and Groves PM. The substantia nigra of the rat: a Golgi study. J Comp Neurol 172: 585-600, 1977.[CrossRef][Web of Science][Medline]

Kalivas PW. Neurotransmitter regulation of dopamine neurons in the ventral tegmental area. Brain Res Rev 18: 75-113, 1993.[CrossRef][Medline]

Kandel ER and Siegelbaum SA. Synaptic integration. In: Principles of Neural Science (4th ed.). New York: McGraw-Hill/Appleton and Lange, 2000, p. 207-228.

Kang Y and Kitai ST. Calcium spike underlying rhythmic firing in dopaminergic neurons of the rat substantia nigra. Neurosci Res 18: 195-207, 1993.[CrossRef][Web of Science][Medline]

Kita H and Kitai ST. Efferent projections of the subthalamic nucleus in the rat: light and electron microscopic analysis with the PHA-L method. J Comp Neurol 260: 435-452, 1987.[CrossRef][Web of Science][Medline]

Kitai ST, Shepard PD, Callaway JC, and Scroggs R. Afferent modulation of dopamine neuron firing patterns. Curr Opin Neurobiol 9: 690-697, 1999.[CrossRef][Web of Science][Medline]

Komendantov AO and Canavier CC. A computational model of a midbrain dopamine neuron in vivo suggests mechanisms by which GABAA afferents regulate the firing pattern (CD-ROM). Soc Neurosci Abstr 290.5, 2001.

Komendantov AO and Canavier CC. Electrical coupling between model midbrain dopamine neurons: effect on firing pattern and synchrony. J Neurophysiol 87: 1526-1541, 2002.[Abstract/Free Full Text]

Koob GF, Le HT, and Creese I. The D1 dopamine receptor antagonist SCH 23390 increases cocaine self-administration in the rat. Neurosci Lett 79: 315-320, 1987.[CrossRef][Web of Science][Medline]

Kotter R and Feizelmeier M. Species-dependence and relationship of morphological and electrophysiological properties in nigral compacta neurons. Prog Neurobiol 54: 619-632, 1998.[CrossRef][Web of Science][Medline]

Lacey MG, Mercury NB, and North RA. Two cell types in rat substantia nigra zona compacta distinguished by membrane properties and the action of dopamine and opioids. J Neurosci 9: 1233-1241, 1989.[Abstract]

Li YX, Bertram R, and Rinzel J. Modeling N-methyl-D-aspartate-induced bursting in dopamine neurons. Neuroscience 71: 397-410, 1996.[CrossRef][Web of Science][Medline]

Liss B, Franz O, Sewing S, Bruns R, Neuhoff H, and Roeper J. Tuning pacemaker frequency of individual dopaminergic neurons by Kv4.3L and KChip3.1 transcription. EMBO J 20: 5715-5724, 2001.[CrossRef][Web of Science][Medline]

Matsuda Y and Fujimura K. Action of habenular efferents on ventral tegmental area neurons studies in vitro. Brain Res Bull 28: 743-749, 1992.[CrossRef][Web of Science][Medline]

McDonald AJ. Glutamate and aspartate immunoreactive neurons of the rat basolateral amygdala: colocalization of excitatory amino acids and projections to the limbic circuit. J Comp Neurol 365: 367-379, 1996.[CrossRef][Web of Science][Medline]

Mercuri NB, Bonci A, Calabresi P, Stratta F, Stefani A, and Bernardi G. Effects of dihydropyridine calcium antagonists on rat midbrain dopaminergic neurons. Br J Pharmacol 113: 831-838, 1994.[Web of Science][Medline]

Mercuri NB, Bonci A, Siniscalchi A, Stefani A, Calabresi P, and Bernardi G. Electrophysiological effects of monoamine oxidase inhibition on rat midbrain dopaminergic neurons: an in vitro study. Br J Pharmacol 117: 528-532, 1996.[Web of Science][Medline]

Mereu G, Costa E, Armstrong DM, and Vicini S. Glutamate receptor subtypes mediate excitatory synaptic currents of dopamine neurons in midbrain slices. J Neurosci 11: 1359-1366, 1991.[Abstract]

Morikawa H, Khodakhah K, and Williams JT. Two intracellular pathways mediate metabotropic glutamate receptor-induced Ca2+ mobilization in dopamine neurons. J Neurosci 23: 149-157, 2003.[Abstract/Free Full Text]

Nedergaard S, Flatman JA, and Engberg I. Nifedipine- and omega-conotoxin-sensitive Ca2+ conductances in guinea pig substantia nigra pars compacta neurons. J Physiol 466: 727-747, 1993.[Abstract/Free Full Text]

Neuhoff H, Neu A, Liss B, and Roeper J. Ih channels contribute to the different functional properties of identified dopaminergic subpopulations in the midbrain. J Neurosci 22: 1290-1302, 2002.[Abstract/Free Full Text]

Obeso JA, Rodriguez-Oroz MC, Rodriguez M, Lanciego JL, Artieda J, Gonzalo N, and Olanow CW. Pathophysiology of the basal ganglia in Parkinson's disease. Trends Neurosci 23: 8-19, 2000.[Medline]

Overton P and Clark D. Burst firing in midbrain dopaminergic neurons. Brain Res Rev 25: 312-334, 1997.[CrossRef][Medline]

Paladini CA, Celada P, and Tepper JM. Striatal, pallidal, and pars reticulata evoked inhibition of nigrostriatal dopaminergic neurons is mediated by GABAA receptors in vivo. Neuroscience 89: 799-812, 1999a.[CrossRef][Web of Science][Medline]

Paladini CA, Fiorillo CD, Morikawa H, and Williams JT. Amphetamine selectively blocks inhibitory glutamate transmission in dopamine neurons. Nat Neurosci 4: 275-281, 2001.[CrossRef][Web of Science][Medline]

Paladini CA, Iribe Y, and Tepper JM. GABAA receptor stimulation blocks NMDA-induced bursting of dopaminergic neurons in vitro by decreasing input resistance. Brain Res 832: 145-151, 1999b.[CrossRef][Web of Science][Medline]

Paladini CA and Tepper JM. GABAA and GABAB antagonists differentially affect the firing pattern of substantia nigra dopaminergic neurons in vivo. Synapse 32: 165-176, 1999.[CrossRef][Web of Science][Medline]

Patneau DK and Mayer ML. Structure-activity relations for amino acid transmitter candidates acting at N-methyl-D-aspartate and quisqualate receptors. J Neurosci 10: 2385-2399, 1990.[Abstract]

Paxinos G and Watson C. The Rat Brain in Stereotaxic Coordinates. Orlando, FL: Academic, 1986.

Penney MS and Britton NF. Modelling natural burst firing in nigral dopamine neurons. J Theor Biol 219: 207-223, 2002.[CrossRef][Web of Science][Medline]

Phillipson OT. Afferent projections to the ventral tegmental area of Tsai and interfascicular nucleus: a horseradish peroxidase study in the rat. J Comp Neurol 187: 117-143, 1979.[CrossRef][Web of Science][Medline]

Ping HX and Shepard PD. Apamin-sensitive Ca(2+)-activated K+ channels regulate pacemaker activity in nigral dopamine neurons. Neuroreport 7: 809-814, 1996.[Web of Science][Medline]

Preston RJ, McCrea RA, Chang HT, and Kitai ST. Anatomy and physiology of substantia nigra and retrorubral neurons studied by extra- and intracellular recording and by horseradish peroxidase labeling. Neuroscience 6: 331-344, 1981.[CrossRef][Web of Science][Medline]

Richards CD, Shiroyama T, and Kitai ST. Electrophysiological and immunocytochemical characterization of GABA and dopamine neurons in the substantia nigra of the rat. Neuroscience 80: 545-557, 1997.[CrossRef][Web of Science][Medline]

Rinzel J and Ermentrout B. Analysis of neuronal excitability and oscillations. In: Methods in Neuronal Modeling—From Ions to Networks (2nd ed.), edited by Koch C and Segev I. Cambridge, MA: MIT Press, 1998, p. 251-291.

Schultz W. Predictive reward signal of dopamine neurons. J Neurophysiol 80: 1-27, 1998.[Abstract/Free Full Text]

Scroggs RS, Cardenas CG, Whittaker JA, and Kitai ST. Muscarine reduces calcium-dependent electrical activity in substantia nigra dopaminergic neurons. J Neurophysiol 86: 2966-2972, 2001.[Abstract/Free Full Text]

Sesack SR and Pickel VM. Prefrontal cortical efferents in the rat synapse on unlabled neuronal target of catecholamine terminals in the nusleas accumbens septi and on dopamine neurons in the ventral tegmental area. J Comp Neural 320: 145-160, 1992.[CrossRef][Web of Science][Medline]

Seutin V, Johnson SW, and North RA. Apamin increases NMDA-induced burst-firing of rat mesencephalic dopamine neurons. Brain Res 10: 341-344, 1993.

Seutin VM, Massotte L, Liegeois J-F, and Scuvee-Moreau J. Effect of a new SK channel blocker, methyl-laudenosine, on the firing pattern of rat midbrain dopaminergic neurons in vivo (CD-ROM). Soc Neurosci Abstr 340.5, 2002.

Seutin V, Scuvee-Moreau J, and Dresse A. Evidence for a non-GABAergic action of quaternary salts of bicuculline on dopaminergic neurones. Neuropharmacology 36: 1653-1657, 1997.[CrossRef][Web of Science][Medline]

Smith Y and Bolam JP. The output neurones and the dopaminergic neurones of the substantia nigra receive a GABA-containing input from the globus pallidus in the rat. J Comp Neurol 296: 47-64, 1990.[CrossRef][Web of Science][Medline]

Smith Y, Charara A, and Parent A. Synaptic innervation of midbrain dopaminergic neurons by glutamate-enriched terminals in the squirrel monkey. J Comp Neurol 364: 231-253, 1996.[CrossRef][Web of Science][Medline]

Spanagel R and Weiss F. The dopamine hypothesis of reward: past and current status. Trends Neurosci 22: 521-527, 1999.[CrossRef][Web of Science][Medline]

Takada M, Kang Y, and Imanishi M. Immunohistochemical localization of voltage-gated calcium channels in substantia nigra dopamine neurons. Eur J Neurosci 13: 757-762, 2001.[CrossRef][Web of Science][Medline]

Tepper JM, Martin LP, and Anderson DR. GABAA receptor-mediated inhibition of rat substantia nigra dopaminergic neurons by pars reticulata projection neurons. J Neurosci 15: 3092-3103, 1995.[Abstract]

Tepper JM, Sawyer SF, and Groves PM. Electrophysiologically identified nigral dopaminergic neurons intracellularly labeled with HRP: light-microscopic analysis. J Neurosci 7: 2794-2806, 1987.[Abstract]

Tokuno H, Moriizumi T, Kudo M, and Nakamura Y. A morphological evidence for monosynaptic projections from the nucleus tegmenti pedunculopontinus pars compacta (TPC) to nigrostriatal projection neurons. Neurosci Lett 85: 1-4, 1988.[CrossRef][Web of Science][Medline]

Wang T, O'Connor WT, Ungerstedt U, and French ED. N-methyl-D-aspartic acid biphasically regulates the biochemical and electrophysiological response of A10 dopamine neurons in the ventral tegmental area: in vivo microdialysis and in vitro electrophysiological studies. Brain Res 666: 255-262, 1994.[CrossRef][Web of Science][Medline]

Weinberger DR. Implications of normal brain development for the pathogenesis of schizophrenia. Arch Gen Psychiatry 44: 660-669, 1987.[Abstract/Free Full Text]

Wightman RM and Zimmerman JB. Control of dopamine extracellular concentration in rat striatum by impulse flow and uptake. Brain Res Brain Res Rev 15: 135-144, 1990.[CrossRef][Medline]

Wilson CJ and Callaway JC. Coupled oscillator model of the dopaminergic neuron of the substantia nigra. J Neurophysiol 83: 3084-3100, 2000.[Abstract/Free Full Text]

Wilson CJ, Young SJ, and Groves PM. Statistical properties of neuronal spike trains in the substantia nigra: cell types and their interactions. Brain Res 136: 243-260, 1977.[CrossRef][Web of Science][Medline]

Wolfart J and Roeper J. Selective coupling of T-type calcium channels to SK potassium channels prevents intrinsic bursting in dopaminergic midbrain neurons. J Neurosci 22: 3404-13, 2002.[Abstract/Free Full Text]

Wu YN and Johnson SW. Pharmacological characterization of inward current evoked by N-methyl-D-aspartate in dopamine neurons in the rat brain slice. J Pharmacol Exp Ther 279: 457-463, 1996.[Abstract/Free Full Text]

Wu YN and Johnson SW. Comparison of burst firing induced by NMDA versus apamin in rat midbrain dopamine neuron (CD-ROM). Soc Neurosci Abstr, 290.19, 2001.

Yung WH, Häusser MA, and Jack JJ. Electrophysiology of dopaminergic and non-dopaminergic neurons of the guinea pig substabtia nigra pars compacta in vitro. J Physiol 436: 643-667, 1991.[Abstract/Free Full Text]




This article has been cited by other articles:


Home page
Proc. Natl. Acad. Sci. USAHome page
L. S. Zweifel, J. G. Parker, C. J. Lobb, A. Rainwater, V. Z. Wall, J. P. Fadok, M. Darvas, M. J. Kim, S. J. Y. Mizumori, C. A. Paladini, et al.
From the Cover: Feature Article: Disruption of NMDAR-dependent burst firing by dopamine neurons provides selective assessment of phasic dopamine-dependent behavior
PNAS, May 5, 2009; 106(18): 7281 - 7288.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
M. Migliore, C. Cannia, and C. C. Canavier
A Modeling Study Suggesting a Possible Pharmacological Target to Mitigate the Effects of Ethanol on Reward-Related Dopaminergic Signaling
J Neurophysiol, May 1, 2008; 99(5): 2703 - 2707.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
C. C. Canavier, S. A. Oprisan, J. C. Callaway, H. Ji, and P. D. Shepard
Computational Model Predicts a Role for ERG Current in Repolarizing Plateau Potentials in Dopamine Neurons: Implications for Modulation of Neuronal Activity
J Neurophysiol, November 1, 2007; 98(5): 3006 - 3022.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
F. W. Hopf, M. Martin, B. T. Chen, M. S. Bowers, M. M. Mohamedi, and A. Bonci
Withdrawal From Intermittent Ethanol Exposure Increases Probability of Burst Firing in VTA Neurons In Vitro
J Neurophysiol, October 1, 2007; 98(4): 2297 - 2310.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
C. C. Canavier and R. S. Landry
An Increase in AMPA and a Decrease in SK Conductance Increase Burst Firing by Different Mechanisms in a Model of a Dopamine Neuron In Vivo
J Neurophysiol, November 1, 2006; 96(5): 2549 - 2563.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
A. S. Kuznetsov, N. J. Kopell, and C. J. Wilson
Transient High-Frequency Firing in a Coupled-Oscillator Model of the Mesencephalic Dopaminergic Neuron
J Neurophysiol, February 1, 2006; 95(2): 932 - 947.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
C. Soto-Trevino, P. Rabbah, E. Marder, and F. Nadim
Computational Model of Electrically Coupled, Intrinsically Distinct Pacemaker Neurons
J Neurophysiol, July 1, 2005; 94(1): 590 - 604.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
C. A. McClung, K. Sidiropoulou, M. Vitaterna, J. S. Takahashi, F. J. White, D. C. Cooper, and E. J. Nestler
Regulation of dopaminergic transmission and cocaine reward by the Clock gene
PNAS, June 28, 2005; 102(26): 9377 - 9381.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
91/1/346    most recent
00062.2003v1
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 (19)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Komendantov, A. O.
Right arrow Articles by Canavier, C. C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Komendantov, A. O.
Right arrow Articles by Canavier, C. C.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2004 by the The American Physiological Society.