J Physiol Physiology in Press
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Physiol Volume 567, Number 3, 869-890, September 15, 2005 DOI: 10.1113/jphysiol.2005.092056
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
567/3/869    most recent
jphysiol.2005.092056v1
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
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 Google Scholar
Google Scholar
Right arrow Articles by Newell, E. W
Right arrow Articles by Schlichter, L. C
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Newell, E. W
Right arrow Articles by Schlichter, L. C

Integration of K+ and Cl currents regulate steady-state and dynamic membrane potentials in cultured rat microglia

Evan W Newell1 and Lyanne C Schlichter1

1 Division of Cellular and Molecular Biology, Toronto Western Research Institute, and Department of Physiology, University of Toronto, Toronto, Ontario, Canada


    Abstract
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
The role of ion channels and membrane potential (Vm) in non-excitable cells has recently come under increased scrutiny. Microglia, the brain's resident immune cells, express voltage-gated Kv1.3 channels, a Kir2.1-like inward rectifier, a swelling-activated Cl current and several other channels. We previously showed that Kv1.3 and Cl currents are needed for microglial cell proliferation and that Kv1.3 is important for the respiratory burst. Although their mechanisms of action are unknown, one general role for these channels is to maintain a negative Vm. An impediment to measuring Vm in non-excitable cells is that many have a very high electrical resistance, which makes them extremely susceptible to leak-induced depolarization. Using non-invasive Vm-sensitive dyes, we show for the first time that the membrane resistance of microglial cells is several gigaohms; much higher than the seal resistance during patch-clamp recordings. Surprisingly, we observed that small current injections can evoke large Vm oscillations in some microglial cells, and that injection of sinusoidal currents of varying frequency exposes a strong intrinsic electrical resonance in the 5- to 20-Hz frequency range in all microglial cells tested. Using a dynamic current clamp that we developed to actively compensate for the damage done by the patch-clamp electrode, we found that the Vm oscillations and resonance were more prevalent and larger. Both types of electrical behaviour required Kv1.3 channels, as they were eliminated by the Kv1.3 blocker, agitoxin-2. To further determine how the ion currents integrate in these cells, voltage-clamp recordings from microglial cells displaying these behaviours were used to analyse the biophysical properties of the Kv1.3, Kir and Cl currents. A mathematical model that incorporated only these three currents reproduced the observed Vm oscillations and electrical resonance. Thus, the electrical behaviour of this ‘non-excitable’ cell type is much more complex than previously suspected, and might reflect a more common oversight in high resistance cells.

(Received 3 June 2005; accepted after revision 8 July 2005; first published online 14 July 2005)
Corresponding author L. C. Schlichter: Toronto Western Research Institute, 399 Bathurst Street MC9-417, Toronto, Ontario, Canada, M5T 2S8. Email: schlicht{at}uhnres.utoronto.ca


    Introduction
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Microglia are immune cells that reside in the central nervous system (CNS). They respond to essentially all insults (damage or disease) by undergoing complex changes, collectively called ‘activation’. As several functions of activated microglia can be harmful to bystander cells, understanding their mechanisms of activation will be important for developing therapies to regulate CNS inflammation. In this study, we test a strategy that uses potent and specific blockers of the voltage-gated K+ channel, Kv1.3 (Kotecha & Schlichter, 1999; Khanna et al. 2001; Schlichter & Khanna, 2002). Kv1.3, which is widely expressed in immune cells, plays a clear role in activation of human natural killer cells and T lymphocytes (Schlichter et al. 1986; Cahalan et al. 2001). Indeed, Kv1.3 appears to be a promising target for autoimmune diseases that involve T-cell activation (Beeton et al. 2001; Chandy et al. 2004). In T cells, one mechanistic role for Kv1.3 channels is to maintain a negative membrane potential (Vm), which promotes calcium influx (Krasznai et al. 1995; Cahalan et al. 2001; Chandy et al. 2004; Robbins et al. 2005). The same might be true for other immune cells. Numerous ion currents have been identified in microglia, including several K+ channels, a Cl current, ionotropic purinergic receptors, an H+ current and the cation channel, TRPM7 (see Eder, 1998; Eder & Decoursey, 2001; Schlichter & Khanna, 2002 for reviews; Cayabyab & Schlichter, 2002; Jiang et al. 2003). Little is known about how ion channels affect microglial cell functions and how Vm is controlled. In whole-cell recordings from cultured microglial cells, the predominant currents are Kv1.3 and a Kir2.1-like inward rectifier (Norenberg et al. 1994a; Visentin et al. 1995; Schlichter et al. 1996; Schilling et al. 2000). However, in comparing freshly isolated with cultured microglia, we showed that activity of Kv1.3, Kv1.5 and Cl channels promotes microglia proliferation (Schlichter et al. 1996; Kotecha & Schlichter, 1999) and have recently found that microglial Kv1.3 channel activity increases their ability to kill neurones (Fordyce et al. 2005). In addition, Kv1.3 and Ca2+-activated K+ channels regulate the respiratory burst of microglia (Colton et al. 1994; Spranger et al. 1998; Khanna et al. 2001), and Cl channels are involved in nitric oxide production (Pyo et al. 1997) and shape changes (Eder et al. 1997).

Why are so many channels expressed in ‘non-excitable’ microglia, and how do the channels integrate to contribute to the cells' functions? One proposed role for the inward rectifier current is to maintain a negative membrane potential (Vm) (Chung et al. 1999), which drives calcium entry when purinergic receptors are activated (Franchini et al. 2004). Clearly, elucidating the integrated roles of ion channels in microglial functions will require accurate measurements of Vm. Previous Vm measurements, which employed whole-cell recordings in the current-clamp mode have yielded diverse values: (i) bimodal distributions at –35 and –70 mV (Norenberg et al. 1994a), –61 and –29 mV, and –52 and –23 mV in brain slices (Boucsein et al. 2003); (ii) a relatively depolarized potential at –32 mV in unstimulated cells but bimodal values at about –70 and –44 mV after activation by lipopolysaccharide (Chung et al. 1999); or (iii) spontaneous jumps between about –35 and –74 mV (Norenberg et al. 1994a). Whole-cell recording itself can affect Vm by artificially altering intracellular ion concentrations and diluting channel regulatory molecules, which can affect the current prevalence, amplitudes and properties. For patch-clamp recordings to accurately measure Vm, the membrane resistance must be substantially lower than the seal resistance. Obvious signs of a leak artifact include more negative Vm values in recordings that have higher resistances, and a relatively linear current–voltage relation under current clamp. Conversely, if the leak allows Ca2+ entry, activation of Ca2+-activated K+ (KCa) channels (Eder, 1998; Khanna et al. 2001) could hyperpolarize Vm. Attempts to overcome leak artifacts include the historical use of sharp electrodes to report the Vm value at the instant of impalement, before the cells depolarize (Schlichter, 1983; Ince et al. 1983). More recently, the cell-attached configuration in voltage-clamp (Verheugen & Vijverberg, 1995) or current-clamp mode (Mason et al. 2005) has been used to non-invasively measure Vm; however, the latter measurements require relatively low membrane resistances. Vm values for microglial cells have not been verified using non-invasive methods.

Using a combination of fluorescent dyes, current-clamp and voltage-clamp recordings and mathematical modelling, we show: (i) several biophysical features of the Kv1.3, Kir2.1 and swelling-sensitive Cl currents in cultured rat microglia; (ii) that microglia have very high resistances in the physiological Vm range; (iii) that Vm measurements and current–voltage relations are easily compromised by leak around the electrode; (iv) that the Cl current contributes to Vm in intact cells; and (v) that these three currents cause strongly non-linear current–voltage relations and are capable of promoting Vm oscillations and sharp frequency selection.


    Methods
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Cells

Microglia were isolated from brain explants of 2-to 3-day-old Wistar rats as previously described (Cayabyab et al. 2000; Khanna et al. 2001). Rat pups were killed by cervical dislocation in accordance with guidelines from the Canadian Institutes of Health Research and the University Health Network. Neopallial brain tissue was digested by agitation in minimal essential medium (MEM; Gibco Invitrogen, Burlington, Ontario, Canada) containing 2.5% trypsin and 100 U DNase I (Pharmacia, Baie d'Urfe, Quebec, Canada) for 30 min at room temperature (21–23°C), and then triturated and digested for a further 30 min until no tissue clumps were visible. This mixture was pelleted, resuspended and seeded into flasks with MEM containing 10% fetal bovine serum (FBS; Gibco BRL) and 100 µM gentamycin (Gibco BRL). After 2 days, cellular debris, non-adherent cells and supernatant were removed and fresh medium was added to the flask. The mixed cultures were allowed to grow for 10 days and then shaken overnight on an orbital shaker at 8–10 Hz in a standard tissue culture incubator. The non-adherent microglia were harvested, suspended in culture medium containing 2% FBS, plated on coverslips for 2 h, and then washed with culture medium. For experiments using flow cytometry and voltage-sensitive dyes, microglia were used immediately after shaking. After plating, some cell batches were treated with 5 ng ml–1 lipopolysaccharide (LPS; Sigma-Aldrich, Oakville, Ontario, Canada) for 12–16 h.

Patch-clamp electrophysiology

Microglia on coverslips were mounted in a perfusion chamber (Model RC-25, Warner Instruments, Hamden, CT, USA) and the tissue culture medium was replaced with an extracellular (bath) solution containing (mM): NaCl 125, KCl 5, MgCl2 1, CaCl2 1, glucose 5 and Hepes 10; pH adjusted to 7.4 with NaOH. Approximately 20 mM sucrose was added to adjust the osmolarity to 280–300 mosmol l–1, measured with a freezing-point depression osmometer (Advanced Instruments, Norwood, MA, USA). When desired, a large swelling-activated Cl current was elicited with a hypo-osmotic solution made by diluting the bath solution with distilled water to 60% of the normal osmolarity. The Cl current remained activated for several minutes after the bath solution was replaced by a solution containing isotonic N-methyl-D-glucamine (NMDG) chloride, which was useful for isolating the Cl current while maintaining the normal Cl gradient. The intracellular (pipette) solution contained (mM): potassium aspartate 100, KCl 40, MgCl2 1, CaCl2 1, EGTA 10, Hepes 10, MgATP 2; pH adjusted to 7.2 with KOH, osmolarity 280–300 mosmol l–1. Where indicated in the Results, 1 nM agitoxin-2 was added to the bath solution to block the Kv1.3 current.

Pipettes with resistances of 4–6 M{Omega} were pulled from borosilicate glass (WPI, Sarasota, FL, USA). A multiclamp 700A patch-clamp amplifier (Axon Instruments Molecular Devices, Sunnyvale, CA, USA) was used in either the current-clamp or voltage-clamp mode (filtered at 5 kHz), and recordings were compensated on-line for capacitance and series resistance. Data were digitized and acquired using a Digidata 1322 board with pCLAMP version 9.0 (Molecular Devices) and were analysed (see below) using Origin version 7.0 (OriginLabs, Northampton, MA, USA). Liquid–liquid junction potentials (e.g. between the NaCl bath and KCl/potassium aspartate pipette solution) were calculated with the utility in pCLAMP, confirmed by measuring the values using a 3-M KCl electrode (Barry & Lynch, 1991), and subtracted before data analysis.

Dynamic current clamp

Dynamic current clamp allows injection of a current whose amplitude is a function of the measured membrane potential (Vm) (Sharp et al. 1993). This method was used to compensate for leak through the electrode seal, as follows. We wrote a simple program in BASIC, running in DOS, and used a Tekmar Labmaster digitizer to rapidly (> 10 kHz) change the amplitude of the injected current as a linear (ohmic) function of Vm whenever it changed. Thus, the electrode leak was compensated for by injecting an opposing current of equal magnitude. The seal resistance was estimated by assuming that all of the current at the Nernst potential for K+ (EK) was due to an ohmic leak around the seal; a reasonable assumption when the swelling-activated Cl current has not been induced. Further validation of this assumption comes from non-invasive measurements wherein the mean membrane resistance of microglia was shown to be ~8 G{Omega} between the resting potential and ~0 mV (see Figs 4D and 5C).



View larger version (38K):
[in this window]
[in a new window]
 
Figure 4.  A voltage-sensitive dye, di-8-ANEPPS, reveals the high membrane resistance of microglia
A, verification of the response speed of di-8-ANEPPS using the whole-cell voltage clamp. A 25-ms long voltage step from –70 to 0 mV (top trace) was applied while measuring di-8-ANEPPS fluorescence of the same cell in normal bath solution (5 mM K+o; middle trace) or with 140 mM K+o (bottom trace). Dye traces are an average of 30 such steps. B, dependence of inward rectifier current (Kir) on extracellular K+ concentration. Currents in normal solution (5 mM K+o; smaller amplitude) or 140 mM K+o were recorded in response to voltage ramps from +120 to –140 mV from a holding potential of –10 mV. C, Kir currents were used to test the rate of change of the external solution. The lower panel indicates when the position of the pipette superfusing the cell was switched from 5 to 140 mM K+ solution for 250 ms. The top panel shows that elevating external K+ produced larger inward Kir currents during voltage-clamp steps (middle traces) between –90 and +70 (in 20-mV increments) from a holding potential of –10 mV. The time course of each change in current was fitted to a mono-exponential, and used to calculate the average time constant ({tau}) for the solution change. Di, a representative recording of di-8-ANEPPS fluorescence from a single microglial cell. Cell depolarization in response to a 250-ms long exposure to 140 mM K+ solution reduced the fluorescence intensity (measured as photodiode current; upper trace). Note the slower fluorescence response to repolarization when the solution was switched back to 5 mM K+. The dashed line indicates the drift in signal over the 25 averaged traces. Lower panel, the position of the superfusing pipette controlling the external K+ concentration, with an inset showing the time course of the resulting change in Kir current (at –50 mV) on the same time scale (lower panel). Dii, for each cell, the time constants ({tau}) of depolarization and repolarization, calculated from mono-exponential curve fits, are shown for depolarization (17.7 ± 2.1 ms) and repolarization (95.4 ± 9.0 ms). Values are mean ± S.E.M. for 22 cells (**P < 0.0001, Student's t test).

 


View larger version (17K):
[in this window]
[in a new window]
 
Figure 5.  Non-invasive dye measurements of membrane potential in microglia
A, initial evidence that a chloride conductance contributes to the membrane potential (Vm). The fast dye, Di-8-ANEPPS, was used in combination with fast solution changes as in Fig. 3. For each of five cells, the mean change in fluorescence (± S.E.M.) in response to 25 switches in the bath solution was divided by the total fluorescence (*P < 0.05; Student's t test). B, calibration of the more sensitive voltage-dependent dye, DiBAC4, using flow cytometry on microglia. First, the membrane potential of all cell batches was set to 0 mV using 1 µM gramicidin, and then varying concentrations of DiBAC4 were applied to the cells for ~20 min before monitoring cell fluorescence with a flow cytometer (see Methods). This representative experiment shows the linear relationship between intracellular dye concentration and mean fluorescence intensity, from which the actual Vm was calculated. C, effects of ion-channel blockers and K+ concentrations on the mean membrane potential (Vm) of microglia. Vm values were determined using DiBAC4 (see Methods and B). Values are mean ± S.E.M. for the number of separate cell batches indicated on each bar (*P < 0.05; ANOVA with Bonferroni correction for multiple comparisons). Bath solutions contained normal saline with 5 mM K+ and 134 mM Cl (control) and normal saline to which compounds were added (10 nM agitoxin-2, 150 µM flufenamic acid, 100 µM NPPB, 100 µM or 1 mM BaCl2) or saline in which the K+ concentration was raised to 30, 55 or 140 mM in exchange for Na+.

 
Electrical resonance measurements

A sinusoidal current waveform was applied to each microglial cell under current clamp using the ‘stimulus file’ feature of pCLAMP. To continually vary the frequency, we used a modified ‘zap’ function whose frequency decreased (forward) or increased (reverse) according to a hyperbolic function (Hutcheon et al. 1996; Hutcheon & Yarom, 2000). This produced a sufficient number of low frequency oscillations to allow analysis of slow responses. When desired, the sinusoidal current stimulus was combined with an offset current (see Results). Impedance plots were created by dividing the fast-Fourier transform (FFT) of the output (mV Hz–1) by the FFT of the input (pA Hz–1) (Hutcheon et al. 1996) using pCLAMP, and then fitting the curves to an equation describing the impedance of a resistance/inductance/capacitance (RLC) circuit where all three components are in parallel, as follows:


{tjp_1102_m1}

(A)
where:


{tjp_1102_m2}

(B)
In this equation, Z is impedance, Fn is the frequency where Z peaks and represents the natural (or resonant) frequency of the circuit and L is the inductance. Q, the quality factor, which is an indication of the strength of the resonance, is the proportion of the energy stored by inductance versus that dissipated by resistance. Bandwidth (BW = Fn/Q) indicates the resonance sharpness, and is the range of frequencies within which the impedance power is 50% or more of the maximum. Both strength and sharpness are indicators of a cell's ability to respond selectively to specific frequency ranges (i.e. cells with higher Q and smaller bandwidth will respond strongly to a smaller range of frequencies). As seen from these equations, in an ideal RLC circuit with all components in parallel, Q will decrease by the addition of a parallel resistance (seal resistance) and bandwidth will increase. In contrast, the resonant frequency should be unaffected, as Fn = ({surd}[1/(LC)])/2{pi} is independent of R.

Voltage-sensitive dye measurements of membrane potential

Two fluorescent voltage-sensitive dyes were used to monitor the membrane potential (Vm). To measure absolute Vm values, we used the slow, monovalent anionic voltage-sensitive dye, DiBAC4 (Molecular Probes, Invitrogen, Ontario, Canada) as previously described (Krasznai et al. 1995). DiBAC4 equilibrates across the cell membrane (i.e. reaches a Nernstian distribution) according to the membrane potential. A standard curve was created relating the intracellular DiBAC4 concentration to the measured fluorescence, as follows. Separate aliquots of freshly shaken, non-adherent microglia were treated with various concentrations of DiBAC4 (0–1000 nM) for 20 min together with 1 µM gramicidin to fully depolarize the cells and the fluorescence was measured using the FL1 emission filter on a FACScan flow cytometer (Becton Dickinson, Ontario, Canada). Then, the fluorescence of test microglial samples exposed for the same time (20 min) was measured using 500 nM DiBAC4 without gramicidin. Interpolation of the mean cell fluorescence onto the calibration curve yielded an apparent internal dye concentration, which was then used to calculate the Nernst potential for DiBAC4 (i.e. Vm).

For measuring fast membrane potential changes, the cells were labelled with 10 µM di-8-ANEPPS and 0.05% Pluronic F-127 at room temperature for 10 min (Bullen & Saggau, 1999) and then washed in Hepes-buffered bath solution (see ‘Patch-clamp electrophysiology’). Di-8-ANEPPS is an electrochromic dye that partitions into the membrane and undergoes fast shifts in excitation and emission spectra when Vm changes (Loew, 1982). Di-8-ANEPPS fluorescence was monitored with an inverted microscope (Diaphot TMD, Nikon, Canada) using a 530-nm (530BP10) excitation filter, a 550-nm dichroic mirror (550DCLP) and a 590-nm emission filter (590ALP), all from Omega Optical (Brattleboro, VT, USA). For each experiment, light emitted from a field containing a single microglial cell was detected with a cooled Avalanche photodiode connected to a Photomax 200 amplifier (Dagan, Minneapolis, MN, USA) filtered at 2 kHz. To rapidly change extracellular solutions, a piezoelectric actuator (EXFO Burleigh, Quebec, Canada) controlled the position of a theta-glass perfusion pipette (pulled to ~300 µm in diameter) positioned near the cell but outside the partially occluded field of the 100 x objective. To ensure that pipette movement did not affect the signal, some recordings from each cell were made near the isosbestic wavelength for di-8-ANEPPS using a 490-nm excitation filter (490DF20, Omega Optical).

Curve fitting, mathematical modelling and simulations

Curve fitting (see Results) was done using the non-linear least-squares reduction and Simplex method in Origin version 6.1 or 7.0. For simulations of membrane potential changes, the differential equations were entered into XPPAUT for Windows (WinPP) (Ermentrout, http://www.math.pitt.edu/~bard/xpp/xpp.html) and solved quantitatively. We used the RK4 method, which produced accurate simulations for a time step of 0.05 ms, and we confirmed that either longer (0.1 ms) or shorter (0.01 ms) time steps produced the same results. In order to add a sinusoidal current to the simulation, we used the same sine function that was used for the experimental stimulus file.

Reagents

DiBAC4 and Di-8-ANEPPS were obtained from Molecular Probes, through Invitrogen. Agitoxin-2 (AgTx-2) was obtained from Alomone Laboratories (Jerusalem, Israel). Flufenamic acid, 5-nitro-2-(3-phenylpropylamino)benzoic acid (NPPB), adenosine 5'-triphosphate, magnesium salt (MgATP) and other salts were obtained from Sigma-Aldrich (Oakville, Ontario, Canada).


    Results
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Measuring the currents to incorporate into a mathematical model

First, we used voltage-clamp analysis to determine the amplitude, voltage-dependence and kinetic properties of the currents in microglial cells, and incorporated these parameters into a mathematical model that describes their contribution to the membrane potential (Vm).

Kv1.3 current

The kinetics of native Kv1.3 currents (Fig. 1A) have previously been described in detail for microglia (Norenberg et al. 1994a) and T lymphocytes (Cahalan et al. 1985; Decoursey, 1990; Pahapill & Schlichter, 1992; Panyi et al. 1995). For the present analysis, we used a model for Kv1.3 (Marom & Levitan, 1994) in which gating is described by six states: four closed states that undergo sequential voltage-dependent transitions to an open state, and an inactivated state that is accessible only from the open state. This model was used to show that inactivation is state-dependent, not voltage-dependent, and then the state model was reproduced using the standard Hodgkin–Huxley formalism (Marom & Abbott, 1994). Before curve fitting the Kv1.3 currents in microglial cells, we subtracted a small linear leak current by assuming that all current at the K+ equilibrium potential was leak. Justification for omitting Cl current in this calculation is presented below.



View larger version (35K):
[in this window]
[in a new window]
 
Figure 1.  Measuring the current parameters used in the mathematical model
A, representative voltage-clamp recordings of the Kv1.3 current. Outward Kv1.3 currents elicited by voltage-clamp steps between –40 and +50 mV, from a holding potential of –100 mV. The smooth curves are n4j fits of eqns (6)–(9) (see Results) used to calculate time constants of opening and inactivation. The maximal conductance in this LPS-treated cell was 9.72 nS, calculated at +50 mV (i). Tail currents were elicited by voltage-clamp steps between –40 and –75 mV (in 5-mV increments) after a 100-ms pulse to +40 mV. The smooth curves are fits to eqns (6)–(9) and were used to calculate the time constant for channel closing (ii). Steady-state values of n (the activation gate) were obtained from the curve fits shown in i and fitted to a Boltzmann equation, yielding a slope of 8.2 mV and a voltage for half-maximal activation (V1/2) of –41.9 mV (iii). Time course of recovery from inactivation at –100 mV. Vm was held at +10 mV for 45 s to fully inactivate the Kv1.3 current, and then stepped for varying times to –100 mV to allow recovery from inactivation. The conductance was normalized to the maximal value after 18-s recovery time, and the time course was fitted with a mono-exponential function, yielding a time constant ({tau}) of 8.7 s (iv). Time constants for opening and closing as a function of voltage were obtained from the curve fits in i and iv. The smooth curve was fitted to the values calculated from the inverse of the sum of {alpha} and ß curves (see eqns (10) and (11)). The parameters obtained for the rate constants, {alpha} and ß, were: Aa, 0.00567; Ab, 41.5; Ac, 3.17; Ba, 0.0117; Bb, 68.6; Bc, 19.0 (v). B, voltage-clamp analysis of the inward rectifier (Kir) current. From a holding potential of 0 mV (chosen to inactivate the Kv1.3 current), 700-ms long voltage ramps were applied between 0 and 140 mV. To reduce the noise, 30 current–voltage traces were averaged. Then, the conductance was normalized to the maximal value at –140 mV, and the curve was fitted to a Boltzmann equation, yielding a slope of 12.4 mV and V1/2 of –81.7 mV. Inset, a representative example of the observed Kir current, overlaid with a trace of the calculated current (eqns (12) and (13)). C, swelling-activated Cl current. A voltage-clamp ramp was applied between –140 and +75 mV after the swelling-activated chloride current was induced with hypotonic (60%, 200 mosmol l–1) bath solution (see Methods). To eliminate Kir, the bath solution lacked potassium, and the holding potential was 0 mV to inactivate the Kv1.3 current. The current–voltage relation was very well fitted by a polynomial equation (smooth curve): I = (6.4)(10–7)(Vm4) + (3.3)(10–4)(Vm3) + (0.0553)(Vm2) + (6.061)(Vm) + 94.6 pA.

 
Kv1.3 currents were fitted from the first point after the end of the capacitance transient using the following equations.


{tjp_1102_m3}

(1)


{tjp_1102_m4}

(2)


{tjp_1102_m5}

(3)


{tjp_1102_m6}

(4)


{tjp_1102_m7}

(5)
n is the probability of each of the four gates being in the activated conformation; j is the probability that the state-dependent inactivation gate is open; and dj/dt depends on n4j (probability of each channel being open). For greater accuracy, activation and inactivation were fitted simultaneously at each voltage (Fig. 1Ai) using a constant-voltage solution to the previous equations, as follows:


{tjp_1102_m8}

(6)


{tjp_1102_m9}

(7)


{tjp_1102_m10}

(8)


{tjp_1102_m11}

(9)

The variable, xj, was used to describe state-dependent inactivation, which changes with time during the voltage pulse. The initial condition was achieved by holding Vm at –100 mV to remove all inactivation (i.e. set j ~1). Then, the current inactivation was fitted during a step to +50 mV (time constant {tau} = 1/{alpha}j). ßj, the inverse of the time constant of recovery from inactivation was calculated as follows (Fig. 1Aiv). Vm was held at +50 mV (to set j ~0) then stepped for variable durations to –100 mV to allow recovery from inactivation, and finally stepped back to +50 mV to determine the extent of recovery. For all curve fits, the maximal conductance (gmax) was determined by fitting the time course of the current at +50 mV (Fig. 1Ai), and then used to normalize the curve (Fig. 1Aiii). This simple method yielded the same results as calculating gmax from the measured peak current, because activation was much faster than inactivation.

From the curve fits to currents elicited by voltage steps between –40 and +50 mV, {tau}n, and n{infty} were estimated after fixing gmax, {alpha}i and ßi to the values obtained above (Fig. 1Ai). Similarly, tail currents at –50 mV and below (Fig. 1Aii) were fitted using an initial value of n = 1, because the maximal number of n gates is expected to accumulate in the open conformation during a 10-ms pulse to +50 mV. Initial values for j (j0) were calculated by solving eqns (6)–(9) at +40 mV and t = 100 ms, then the closing time constants were obtained from curve fits that included the continuing inactivation. The closing time constant for n gates is equivalent to a mono-exponential decay with a time constant of 4/{tau}n.

The n{infty} values obtained from fitting the activation curves (Fig. 1Ai) were fitted to a standard Boltzmann equation with a fixed minimum of 0 and maximum of 1 (Fig. 1Aiii). These values, combined with the {tau}n values, were used to construct curves describing {alpha}n and ßn as follows (see Fig. 1Av).


{tjp_1102_m12}

(10)


{tjp_1102_m13}

(11)

The curve fits for activation, inactivation and closing obtained from the n4j model accurately described the voltage-clamp traces (Fig. 1Ai and ii), except for small deviations near the threshold for activation of the Kv1.3 current, where the kinetics of the rising phase were closer to mono-exponential. However, we determined that this small deviation did not affect the simulations. That is, when opening was modelled as one fast transition instead of four (i.e. nj, instead of n4j), the results (including Vm oscillations) were very similar. Steady-state inactivation–voltage relations were examined on many cells (not shown, but see Kotecha & Schlichter, 1999) but because inactivation was modelled as a voltage-independent process, this information was not required for parameter estimation. Nevertheless, the voltage-independent parameters, {alpha}j and ßj, were estimated as described above and both inactivation and recovery from inactivation were included in the model and in all simulations. We found that inactivation played very little role in any of the simulations. This is not surprising, because the rate and extent of inactivation depend on the proportion of channels activated, which is small for the durations and voltages used in the simulations. Although one paper reported more complicated kinetics for recovery from inactivation (Levy & Deutsch, 1996), the differences were small and we found that the simulated Vm oscillations were not sensitive to this parameter.

Inward rectifier (Kir) current

The inward rectifier (Kir) current (Fig. 1B) is very similar to the cloned Kir2.1 channel, and Kir2.1 transcripts have been detected in microglia (Schilling et al. 2000). Although the Kir current is strongly inwardly rectifying, a small outward component is seen above EK. We approximated Kir ‘activation’ as time independent, which is justified because Kir2.1 is gated on a millisecond time scale by voltage-dependent relief of block by cations, including Mg2+, spermine and spermidine (Matsuda et al. 1987; Xie et al. 2002). Ionic block of Kir channels has been exploited in cell functional studies wherein Ba2+ (or less commonly, Cs+) is added to the bath solution (Schlichter et al. 1996; Franchini et al. 2004); however, high concentrations are needed to reduce outward current because the block is strongly voltage dependent. Our simulations did not extend more negative than EK (–85 mV), thus the slow voltage-dependent block by external Na+ at very negative potentials was not an issue. As no time-dependent changes were observed during voltage steps, when recording Kir current during voltage ramps its conductance–voltage relation was well described by a Boltzmann equation:


{tjp_1102_m14}

(12)


{tjp_1102_m15}

(13)
Since steady-state inactivation of the Kv1.3 current in microglia is complete at 0 mV (tested on same cells, data not shown; and Kotecha & Schlichter, 1999); a holding voltage of 0 mV was used to isolate the Kir current. Currents evoked by 20 consecutive voltage ramps between +40 and –160 mV were averaged before calculating the conductance in order to eliminate the exaggerated calculation errors near EK caused by current noise. An overlay of the recorded current (Fig. 1B, inset) and that calculated from eqns (12) and (13) shows excellent agreement.

Swelling-activated chloride current

Under the whole-cell recording conditions used in the present study, the Cl current (Fig. 1C) was very small and extremely difficult to distinguish from the leak. Thus, in order to determine its current–voltage relations and kinetics we used a hypotonic solution to exploit its swelling-induced activation (Schlichter et al. 1996). For these measurements, Kir current was eliminated at all physiologically relevant voltages by using a K+-free bath solution, and Kv1.3 was inactivated using a holding potential of 0 mV. The observed reversal potential of the Cl current was about –20 mV, as expected from the chloride gradient (40 mM Cl and 100 mM aspartate in the pipette; 134 mM Cl in the bath). As previously observed (Schlichter et al. 1996), the Cl current was neither voltage activated nor time dependent during voltage-clamp steps (not shown); thus, the current–voltage relationship was fitted to a polynomial (see legend, Fig. 1C).

Current-clamp recordings

First, we attempted to use whole-cell, current-clamp recordings to measure the resting membrane potential (Vm) and the Vm responses to current injections. Our measured Vm values (from –5 to –85 mV, n > 40) were highly variable, as were the responses to current injection; however, small current injections often evoked large Vm changes, indicating a high resistance. The membrane resistance, calculated by dividing the holding current at EK by the driving force (85 mV) for the leak, was ~2.5 G{Omega} for the cell in Fig. 2A, with a mean of 2.1 ± 0.5 G{Omega} recorded in six microglia. This value was not significantly different from the seal resistance (Rseal) measured in the cell-attached configuration in the same six cells (3.2 ± 0.8 G{Omega}). Consequently, the Vm of microglia, as in other high-resistance cells (Ince et al. 1983) including some neurones (Tyzio et al. 2003; Pun & Kleene, 2004; Wang et al. 2003), will be susceptible to electrode-induced leak even when the seal resistance is in the gigaohm range. For example, with an electrode-induced leak as small as 10 pA, the true membrane potential of the cell in Fig. 2A would have been nearly –80 mV (Fig. 2Aii), rather than the recorded value of –64 mV (Fig. 2Ai). There would be 10 pA of leak current in the –60 to –80 mV range if Rseal was 6–8 G{Omega}: a realistic value (see below). Using a –10-pA holding current to maintain a more negative Vm had no effect on the largely leak-dominated voltage–current relationship (Fig. 2Aii).



View larger version (35K):
[in this window]
[in a new window]
 
Figure 2.  Non-linear relationship between current and voltage in microglia reveals high-resistance regions
A, representative current-clamp recordings of membrane potential (Vm) from a microglial cell that had small Kv1.3 and inward rectifier (Kir) currents (see text). The uncorrected resting potential was –64 mV and Vm changes were evoked by steps of current (see inset) superimposed on a holding current of 0 pA (left) or –10 pA (right) (i). Voltage–current relations from the recordings in i, with a holding current of 0 pA ({blacksquare}) or –10 pA ({circ}). Vm values were measured at the 1000 ms time point and plotted against the injected current (ii). B, relationship between voltage- and current-clamp recordings in a representative untreated microglial cell that was also used for comparison in Fig. 3B. Under voltage-clamp, a voltage ramp from –120 to +80 mV was applied from a holding potential of –80 mV. This cell had relatively little Kir (inward current below –80 mV) and Kv1.3 (outward current above –40 mV). Note the extremely small current between –80 and –40 mV, indicating a total resistance of about 2.6 G{Omega}(i). (ii) Conventional current clamp. The traces are Vm responses to steps of injected current between –50 and +15 pA (in 5-pA increments), with a holding current of –40 pA, chosen to set the initial Vm to the Nernst potential for K+ (EK; about –85 mV).

 
In order to understand how the Kv1.3, Kir and Cl currents and inherently high membrane resistance of microglia affect current-clamp recordings of Vm, we made voltage-clamp and current-clamp recordings on some of the same cells. Figure 2B shows an untreated microglial cell which had small K+ currents. In response to a voltage ramp stimulus, the current–voltage relation was highly non-linear, with an extremely high membrane resistance between about –80 and –40 mV (~2.6 G{Omega}) where most Kv1.3 and Kir channels would be closed. Based on the observed high resistance in the region –80 to –40 mV, the Cl current must be very small, and this was verified by applying a Cl channel blocker. That is, at –85 mV, 150 µM flufenamic acid caused a non-significant change in holding current, from –40 ± 4.8 pA, which corresponds with a membrane resistance of 2.0 ± 0.4 G{Omega} (n = 6), to –34 ± 5.5 pA, which corresponds to 2.5 ± 0.8 G{Omega}. The flufenamic acid-sensitive current ranged from +18 to –1 pA, which was not significantly different from zero (i.e. when calculated for each cell and averaged; P = 0.1, one-sample t test). As the Cl current does not show voltage-dependent gating and its whole-cell current is nearly linear at negative potentials (Fig. 1C), it is easily mistaken for, or masked by, the electrode-induced leak. For the cell in Fig. 2Bii, small current injections produced large, nearly mono-exponential Vm excursions. Thus, opening of only a few channels would cause a substantial change in this Vm range. In contrast, at non-physiologically negative Vm values the membrane resistance was low and nearly constant (i.e. ~130 M{Omega} below –80 mV in this cell), as expected when the Kir conductance is high (see Fig. 1C). Above –40 mV, the resistance progressively decreased due to Kv1.3 channel opening, thus the evoked Vm deflections were reduced. In Fig. 2Ai, note the dramatic time-dependent depolarization when Vm reached about –35 mV, consistent with time-dependent Kv1.3 inactivation (Fig. 1A) and a consequent increase in resistance with time. As expected, the current–voltage relation (Fig. 2Aii) was not changed when the initial Vm was set to –80 mV by a –10 pA holding current.

Dynamic current clamp to compensate for electrode damage

To reduce the error caused by electrode damage, we employed a dynamic current clamp, which counteracts the leak by electronically adding an equal and opposite parallel resistance to that of the Rseal (see Methods). The simple assumptions were that the electrode-induced leak conductance was ohmic (i.e. constant with voltage), reversed at 0 mV and carried almost all of the current at EK. We recognize that dynamically compensating for leak through the patch-clamp seal may not be exact, as Rseal must be estimated and the assumption that the current at EK is due to leak is accurate only if the net ionic current is zero at EK. As shown in the previous section, these assumptions were well met under normal (non-swelling) recording conditions (i.e. there was very little current at EK).

To demonstrate the ability to subtract a parallel resistance (e.g. Rseal) we used a model cell (PATCH-1 U, Axon Instruments) with an effective resistance (Rm) of 500 M{Omega} and a membrane capacitance (Cm) of 33 pF. We chose to subtract a parallel Rseal value of 1 G{Omega} with the dynamic current clamp. The resulting Rm should be 1 G{Omega} (i.e. = 1/(1/500 M{Omega} – 1/1000 M{Omega})). As shown in Fig. 3A, an instantaneous step in command current elicited an instantaneous increase in applied current, followed by a time-dependent change in applied current as Vm changed. At the end of the command current step, the instantaneous decrease in applied current was followed by a time-dependent decrease that mirrored the change in Vm. Note that increasing the Rm increases the time constant for charging Vm because {tau} = RmCm (data not shown).



View larger version (35K):
[in this window]
[in a new window]
 
Figure 3.  Seal resistance compensation using dynamic current clamp
A, use of a model cell (Patch-1 U, Axon Instruments; Cm, 33 pF; Rm, 500 M{Omega}) to demonstrate parallel resistance compensation using dynamic current clamp (for details, see Methods). While compensating for a 1-G{Omega} seal resistance (Rseal), current steps from –100 to +100 pA were applied in 20-pA increments (upper traces). Middle traces show the actual applied current (IAPP) injected as a consequence of the resistance compensation (IAPP = ICMDVm/Rseal, where ICMD is the command current and Rseal is 1 G{Omega}). The bottom panel shows the Vm response of the model cell. B, recording from a microglial cell (same cell as in Fig. 2B) to show how dynamic current clamp can explain the existence of two quasi-stable resting potentials. The compensation current set the initial Vm to –85 mV (as in Fig. 2B), at which the seal resistance was calculated to be ~2.5 G{Omega}. Current steps from 0 to 25 pA were applied in 2.5-pA increments (top traces). The actual applied current (IAPP, middle traces) is a function of both the command current (ICMD) and the cell Vm (lower traces) as in A. C, current–voltage relations from the recordings using conventional current clamp (from Fig. 2B, {blacksquare}) and dynamic current clamp (from B, {circ}) on the same cell. Vm values were measured at 1000 ms.

 
A comparison of conventional (Fig. 2Bii) and dynamic current clamp (Fig. 3B bottom traces) in the same cell shows dramatic effects of compensating for the leak through the 2.2-G{Omega} seal resistance. The initial Vm was similar (about –87 mV, which is close to EK) and small current steps produced small Vm changes that were nearly exponential with time. However, when Vm reached about –62 mV, there was a dramatic depolarization to about –35 mV (approximately ECl) and a clear inflection point. It is important this threshold effect was never seen with a model circuit (Fig. 3A), but was sometimes observed in microglia under conventional current clamp. A sharp depolarization is expected as Kir channels close, and has been observed in other cells with Kir currents (Gallin & Livengood, 1981; Goodman & Art, 1996b). Above the all-or-nothing depolarization, Vm entered another low-resistance range, dominated by Kv1.3. This behaviour might account for the commonly observed bimodal resting Vm values at about –35 and –70 mV (Norenberg et al. 1994a; Chung et al. 1999; Boucsein et al. 2003) and occasional spontaneous jumps between these values (Norenberg et al. 1994a).

Again, these Vm responses and the resulting current–voltage relations (Fig. 3C) are consistent with the underlying currents. At very negative Vm values both current-clamp methods show a relatively low Rm due to open Kir channels. In contrast, at Vm values between –70 and –40 mV, there was a dramatically higher Rm (steeper slope) under dynamic current-clamp conditions. Above about –40 mV, both methods revealed a lower Rm region due to open Kv1.3 channels. Concomitantly, in response to current steps (Fig. 3C), the Vm changes above –40 mV reflect a decrease in Rm as the Kv1.3 channels open, then an increase in Rm as they inactivate with time. Of most importance, although a normal seal resistance was obtained in this cell (~2.2 G{Omega}), the current–voltage relation under standard current clamp was linearized by the seal resistance, and this masked the dynamic response seen in Fig. 3B.

Non-invasive measurements of membrane resistance

To more accurately measure the resistance of microglia, we devised a non-invasive method using a single photodiode and the extremely fast voltage-sensitive dye, di-8-ANEPPS (Loew, 1982). To demonstrate the properties of di-8-ANEPPS, Fig. 4A shows the fluorescence response of a single microglial cell to a voltage clamp step (whole-cell configuration) from –70 to 0 mV in normal bath solution (5 mM K+, middle trace) or 140 mM K+ (bottom trace). Clearly, in both solutions, the dye response is faster than the filtering frequency used (2 kHz). This approach exploited a rapid switch in external K+ concentration, which changed the driving force and caused a large (nearly instantaneous) increase in the inward Kir current at negative Vm values (Fig. 4B). Further validation in a microglial cell (Fig. 4C) is illustrated by measuring the di-8-ANEPPS voltage dependence and response time when the perfusion solution was changed with a fast piezoelectric device (see Methods) while monitoring the Kir current. Voltage-clamp steps of 800-ms duration were applied while rapidly switching the bath solution between 5 and 140 mM K+-containing saline. From the change in membrane current, we found that the time constant ({tau}) of solution exchange was 13.1 ± 2.3 ms when the bath was switched to 140 mM K+ and 19.8 ± 1.9 ms when it was switched to 5 mM K+ (n = 4 cells). As shown in Fig. 4A, the dye response is much faster than these values.

Finally, the di-8-ANEPPS responses were recorded from 22 individual microglial cells (e.g. Fig. 4Di), and the values of {tau} were calculated from mono-exponential curve fits (Fig. 4Dii). Notably, {tau} was 17.7 ± 2.1 ms during depolarization (with 140 mM K+), which is significantly faster than the 95.4 ± 9.0 ms value measured during repolarization (with 5 mM K+). For five of these cells, the fluorescence response of each intact cell was measured, then a whole-cell recording was immediately established and Cm was measured and used to calculate membrane resistance (Rm = {tau}/Cm). The value (8.8 ± 1.4 G{Omega}) during repolarization was almost identical to the 7.7 G{Omega} calculated from 25 additional whole-cell recordings of Cm (12.4 ± 1.5 pF) and the mean time constant obtained from the non-invasive di-8-ANEPPS measurements above ({tau} = 95.4 ms). It was not possible to calculate the membrane resistance during the depolarization phase (Fig. 4Di) because the observed response was nearly as fast as the kinetics of solution exchange (Fig. 4C). Nevertheless, the faster rate of change indicates a much lower resistance during depolarization. This difference is consistent with the position of each steady-state membrane potential on the current–voltage relation (e.g. Fig. 4B). That is, at the beginning of the depolarization caused by 140 mM K+, Rm is expected to be low because of the large Kir conductance, and hence the resulting time constant will be shorter. Conversely, after Vm is depolarized and the Kir conductance is very low (Rm high), the time constant for repolarization by low external K+ concentration will be longer.

Non-invasive measurements of membrane potential

Next, we used voltage-sensitive dyes to ascertain the contribution to the resting membrane potential (Vm) of the three ionic currents we had recorded: Kv1.3, Kir and Cl. To assess whether the Cl conductance contributes (Fig. 5A) we used di-8-ANEPPS (as in Fig. 4) and fast changes in the external Cl concentration (from 134 to 50 mM) which changes ECl to about 0 mV. The result was compared with switching external K+ from 5 to 140 mM, which changes EK to about 0 mV. Either condition is expected to depolarize the cells if the ion is significantly permeant. Switching to high external K+ concentration significantly decreased the fluorescence signal, indicating depolarization, and lowering external Cl concentration from 150 to 50 mM caused a smaller, yet significant depolarization. Although the fast dye, di-8-ANEPPS, was useful for avoiding intracellular ion changes that can occur with longer-term changes in extracellular ions, it was not useful for obtaining absolute Vm values, because it has a very low Vm sensitivity (< 5% change per 100 mV) and we observed a continuous drift in signal intensity (see dashed line in Fig. 4Di) which made calibration very difficult.

To obtain steady-state Vm values, we used DiBAC4, which trades off a high Vm sensitivity with slow response kinetics. An advantage is that it can be calibrated and used with flow cytometry (Krasznai et al. 1995) to report Vm values in a large number of cells individually. Figure 5B shows the dye calibration procedure for a representative batch of microglia exposed to normal external ion concentrations. Vm was first set to 0 mV using the ionophore, gramicidin, and then the cells were exposed to a range of external DiBAC4 concentrations, which were allowed to equilibrate across the membrane according to a Nernstian distribution. After using flow cytometry to measure the fluorescence intensity of 5000 cells at each DiBAC4 concentration, the average cellular fluorescence was plotted to create a standard curve. For experiments, the linear standard curve was used with the fluorescence measurements from microglia without gramicidin (Fig. 5C). In control bath solution, the mean Vm was –42 ± 2 mV (i.e. near the activation threshold of Kv1.3 current). Thus, we predicted that the Kv1.3 current would be very small, and concordantly, we observed that the potent Kv1.3 blocker, agitoxin-2 (Garcia et al. 1997) had almost no effect on Vm (P = 0.34). Previous patch-clamp studies have reported resting Vm values of –54 mV (Kettenmann et al. 1990) or two stable states at about –35 and –70 mV for untreated microglia (Norenberg et al. 1994a; Chung et al. 1999; Boucsein et al. 2003), and about –48 mV for lipopolysaccharide-stimulated microglia (Chung et al. 1998).

At first surprisingly, as it has been generally assumed that Kir is responsible for the negative resting Vm of microglia (Franchini et al. 2004), we found that even a very high concentration (1 mM) of the Kir blocker, Ba2+, did not significantly affect the Vm. However, while inward Kir current in microglia is very sensitive to block by external Ba2+ (Schlichter et al. 1996; Franchini et al. 2004), the voltage dependence of block means that much higher Ba2+ concentrations are needed to block outward Kir current. A further important result (Fig. 5C) refutes the common assumption that raising external K+ will depolarize Vm to the new EK value. This will not happen if other ion conductances are significant. We observed (Fig. 5C) that Vm did not reach EK when external K+ was elevated (i.e. to 30 mM (EK –39 versus Vm –19 mV), 55 mM (EK –23 versus Vm –15 mV) or 140 mM (EK 0 versus Vm –5.3 mV)). This finding is consistent with the lack of effect of agitoxin-2 and Ba2+. Instead, we found that both Cl channel blockers caused substantial depolarizations: to –14 ± 3 mV (NPPB) or –17 ± 8 mV (flufenamic acid), thus Cl channels contribute to the resting potential. Despite this, it is not surprising that isolating the Cl current was difficult in whole-cell recordings under non-swelling conditions. That is, a very small Cl current could account for this contribution to Vm when it is in the high-resistance region around –42 mV (for example, it is 6 pA when Rm is ~8 G{Omega}). Moreover, the Cl current could be larger in intact cells if it is regulated by second messenger pathways that are compromised by whole-cell recording.

Electrical resonance in microglia

The two K+ currents we recorded in microglia, Kv1.3 and Kir, are very similar to those in low-frequency cochlear hair cells, where they are responsible for electrical resonance (Crawford & Fettiplace, 1981; Goodman & Art, 1996a); that is, the ability to discriminate input frequency. Thus, we tested for resonance in microglia by applying sinusoidal currents of continuously varying frequency and recording the Vm response under whole-cell current clamp. For the cell shown in Fig. 6Ai, a holding current of (–15 pA) was applied to maintain the initial Vm near the Nernst potential for K+ (–85 mV). Then, with a 0-pA offset, the Vm returned to about –42 mV; the same value observed from non-invasive dye measurements (Fig. 5C). From this Vm, a sinusoidal stimulus (± 5 pA) was injected and the Vm response was used to quantify the electrical resonance. Note that immediately after the depolarizing current step the Vm oscillations were heavily damped. Then, as the frequency decreased to between 10 and 15 Hz, larger oscillations occurred in phase with the stimulus. At lower frequencies the in-phase Vm oscillations were smaller. To quantify this resonance effect, impedance plots were constructed (see Methods) by dividing the FFT of the output signal (Vm) by the FFT of the input signal (injected current) for each experiment. Impedance was plotted as a function of the frequency and fitted with an equation describing an RLC circuit with all components in parallel (see Methods, eqn (A)). Under conventional current clamp, the values from this curve fit were: resonant frequency, 11.2 Hz; peak impedance, 2.07 G{Omega}; quality factor, 0.76 (Fig. 6Aii). These values were the same whether the applied current frequency increased or decreased during the stimulus (not shown).



View larger version (32K):
[in this window]
[in a new window]
 
Figure 6.  Microglia display electrical resonance that requires Kv1.3 current
Left-hand panels, a sine wave current of decreasing frequency (upper traces in each panel) was injected into the same LPS-treated microglial cell under conventional current clamp (A) or dynamic current clamp (B), and the resulting membrane potential (Vm) responses (lower traces in each panel) were monitored. Ai, from a holding current of –15 pA, the offset was switched to 0 pA and a sine wave of amplitude ± 5 pA was injected. Aii, the right panel for each experiment shows an impedance plot, fitted to eqn (A) (see Methods) to obtain resonant frequency, peak impedance, bandwidth and quality factor. B, electrical resonance is enhanced by using the dynamic current clamp to compensate for the leak resistance. i, from a holding current of 0 pA (Vm –83 mV), the offset current was switched to +40 pA, and a sine wave of amplitude ± 5 pA was injected (as in A). ii, the corresponding impedance plot and fitted function (smooth curve; eqn (A); see Methods) display a narrower peak at a lower frequency and have higher peak impedance than with conventional current clamp (A). iii, hyperpolarization can eliminate the resonance. From the same holding current (0 pA; Vm, –84 mV), a smaller offset current was applied (+30 pA), along with the same sine wave (± 5 pA) as in A. iv, the lack of a peak in the corresponding impedance plot indicates there is no resonance.

 
Many factors can influence the electrical resonance properties of a cell. These include Cm and the time constants of channel opening and closing, which are often a function of voltage. As in previous reports on Kv1.3 (Cahalan et al. 1985; Marom & Levitan, 1994), we found that the kinetics of activation and closing of Kv1.3 channels in microglia have a very steep voltage dependence with relatively slow closing kinetics, and the slowest closing time constants occur between –40 and –60 mV. It is important to note that for every cell tested, the resonant frequency and even the occurrence of resonance was affected by varying the mean Vm (i.e. varying the offset current). Electrical resonance was only observed in a microglial cell if the Vm was between –50 and –35 mV, and within this range, the resonant frequency increased with depolarization (five cells tested; not shown).

The high Rm of microglia makes the leak conductance a particularly powerful variable that should alter the quality and maximal impedance. The measured resonant frequency (Fn) should not be affected by the leak conductance because the cell acts like a parallel RmLCm circuit, where Fn = ({surd}[1/LC])/2{pi}. We next used the dynamic current clamp (see Methods and Fig. 3) to assess the effects of compensating for the leak conductance (Fig. 6Bi). While applying the dynamic current clamp, the holding current was stepped to +40 pA to set the initial Vm to about –46 mV, and then a sinusoidal stimulus of ± 5 pA amplitude was applied, as before. The Vm response (Fig. 6Bi) was more complicated than under conventional current clamp (compare Fig. 6Ai). That is, immediately after the depolarizing current step, the endogenous oscillations were damped because they were out of phase with the command current oscillations. Then, for stimulation frequencies between about 6 and 8 Hz, there were large Vm responses in phase with the stimulation. At lower frequencies, the Vm responses were smaller and again, were out of phase. Thus, with dynamic current clamp, phase synchrony occurred over a much narrower frequency range (compare Fig. 6Aii with Fig. 6Bii). It is expected that compensating for Rseal will increase the total Rm and peak impedance, as well as the quality factor value (i.e. Q = R/(2{pi}FnL); see Methods, eqn (B)), while decreasing the bandwidth. Accordingly, in Fig. 6Bii, when the Vm responses under dynamic current clamp were fitted to the all-parallel RLC impedance function (Methods, eqn (A)), the maximal impedance increased (from 2.2 to ~3 G{Omega}), the frequency tuning sharpened (bandwidth decreased from 14.7 to 5.0 Hz), and the quality factor increased (from 0.76 to 1.4). The resonant frequency was lower in this example under dynamic current clamp (7.2 Hz), exactly as we expected from the hyperpolarized Vm and our observation that resonant frequency decreases with hyperpolarization. In principle, if Rseal is overcompensated under dynamic current clamp, the bandwidth and maximal impedance would decrease and the quality factor would increase; however, as illustrated below (see Fig. 8), the inaccuracies are much smaller than if leak compensation is omitted. As for conventional current clamp, resonance only occurred under dynamic current clamp when the mean Vm was in a permissive range. For instance, Fig. 6Biii shows that for the same sinusoidal stimulus as in Fig. 6A (or larger amplitude stimuli, not shown), no resonance peak was observed in the impedance plot (Fig. 6Biv) when a ‘subthreshold’ holding current (+ 30 pA) was applied to set the mean Vm to –77 mV.



View larger version (26K):
[in this window]
[in a new window]
 
Figure 8.  The quantitative model recapitulates the dynamic membrane potential responses
A, electrical resonance is recapitulated by the quantitative properties of Kv1.3 and Kir currents in microglia. i, simulated current injection of +45 pA and a superimposed sine wave (amplitude, ± 5 pA) of decreasing frequency (see Methods; upper trace). The calculated Vm response based on the actual conductances and kinetics of the Kv1.3 and Kir currents in the microglial cell shown in Fig. 6 (lower trace). ii, impedance plot of the data, calculated by dividing the FFToutput by the FFTinput and fitting to eqn (A) (see Methods). B, calculated membrane potential (Vm) response to a simulated current injection of +100 pA for a microglial cell with relatively large Kv1.3 and Kir currents, and a modest Cl current. These simulations used actual parameters from the cells shown in Figs 1 and 6A, determined as described in the text. The maximal conductances were 9.72 nS for Kv1.3, 22.3 nS for Kir and 0.5 nS for Cl, and the membrane capacitance was 22.5 pF. C, calculated Vm response to a simulated current injection (+10.5 pA) for a microglial cell with small Kv1.3 and Kir conductances. Values used for this simulation, obtained from the cell shown in Fig. 2B were maximal conductances of 1.5 nS for Kv1.3, 3.6 nS for Kir, 0.25 nS for ICl and a membrane capacitance of 10 pF. D, the time course of Kv1.3 current (thin trace), Kir current (dotted trace) and membrane potential during a single Vm oscillation (thick trace). The values were calculated from the simulation in B.

 
Membrane potential oscillations in microglia

To our surprise, when current pulses were injected to examine current–voltage relations and membrane resistance (as in Fig. 2), occasionally Vm oscillations were seen that were damped (Fig. 7A) or nearly stable (Fig. 7Bi and Cii). The oscillatory response appeared to require a very negative initial Vm and a permissive range of injected current as larger current stimuli caused faster damping (e.g. Fig. 7Aii), and some cells displayed only damped oscillations. In cells that displayed Vm oscillations, the recording mode was subsequently switched to voltage clamp and the ion currents were quantified. Under conventional current clamp, large, stable oscillations were elicited only in cells that had large Kv1.3 and Kir currents. Kv1.3 was necessary as the blocker, agitoxin-2, stopped the oscillations that could be repeatedly evoked by current injection before adding blocker. Figure 7B shows an example wherein stable Vm oscillations (Fig. 7Bi) were prevented after AgTx-2 was applied (Fig. 7Bii) and could no longer be evoked by a wide variety of depolarizing current injections. Instead, even a small depolarizing current evoked a dramatic depolarization, with an inflection point and all-or-nothing response at about –45 mV in this cell. At the end of the stimulus, the cell re-polarized only when a large hyperpolarizing current was injected, after which Vm remained very negative because the large Kir current was activated. Thus, Kv1.3 was necessary for the Vm oscillations and repolarization after depolarizing events, while Kir could maintain a very negative membrane potential. As the permissive voltage range falls within the high resistance range, oscillations and resonance will be easily shunted by the electrode-induced leak. This is a plausible explanation for the infrequency of observing spontaneous oscillations in patch-clamp recordings.



View larger version (30K):
[in this window]
[in a new window]
 
Figure 7.  Microglia can produce Vm oscillations that require Kv1.3 current
A, conventional current-clamp recordings from an LPS-treated microglial cell that had large Kv1.3 and Kir currents, and an apparent resting potential of –82 mV i, Vm oscillations were elicited when a +60-pA current was injected for the duration of the recording. ii, injecting a larger current (+80 pA) elicited oscillations that were rapidly damped during the current injection (shown by the bar). B, under conventional current clamp, Vm oscillations are prevented by the Kv1.3-channel blocker, agitoxin-2. i, quasi-stable Vm oscillations were elicited by a depolarizing current (indicated by bar). ii after adding 1 nM agitoxin-2 to the bath, no Vm oscillations were elicited by depolarizing currents (+20 pA shown). Injection of –100 pA produced a rapid hyperpolarization of Vm. C, the dynamic current clamp can expose Vm oscillations in response to current injections. i, a conventional current clamp recording with a holding current of –40 pA, chosen to set the initial Vm to –85 mV. Currents were injected between –5 pA and +20 pA, in 5-pA steps. Note the inflection point at about –60 mV. ii, a dynamic current-clamp recording, wherein compensation for the seal resistance set the initial Vm to –85 mV. Depolarizing current steps from +50 to +65 pA in 5-pA increments show quasi-stable Vm oscillations and an inflection point at about –58 mV.

 
Although spontaneous Vm oscillations were seen in only three cells under conventional current clamp, they were seen in three formerly non-oscillating cells when recording was switched to dynamic current clamp. Of importance, these Vm oscillations even occurred in cells with smaller Kv1.3 and Kir conductances, such as the cell shown in Fig. 7C, which nevertheless displayed a very negative initial Vm (about –78 mV). In this cell, under standard current clamp (Fig. 7Ci), the Vm responses to a series of depolarizing current steps showed an inflection at about –60 mV but, at most, only very small damped oscillations. In contrast, under dynamic current clamp (Fig. 7Cii), which did not change the initial Vm, this cell displayed stable Vm oscillations in response to sufficient depolarizing current. The Vm responses of cells with small K+ conductances are more likely to be short-circuited by the electrode-induced leak.

The quantitative model predicts electrical resonance

For all simulations, the cell was modelled as a single compartment and a uniform distribution of each channel type was assumed. Changes in membrane potential (Vm) with time were calculated from ionic and injected currents by numerically solving the conservation of current equation, as follows.


{tjp_1102_m16}

(14)
Cm is the capacitance of the cell and IApp simulates the injected current. Iseal was used in selected simulations to mimic damage caused by the patch electrode. It is a simple ohmic current of the form: Iseal = Vm/Rseal, where Rseal is the seal resistance. IKv and Iir were defined above (eqns (1)–(5) and (10)–(13)). First, the parameters obtained from the cell shown in Fig. 6 were used to test the hypothesis that the two prevalent K+ channels in microglia can cause the observed electrical resonance. The maximal conductance parameters for the simulation in Fig. 8A were Kv1.3: gKvmax, 9.4 nS; Kir: gKirmax, 7.3 nS. At the highest frequencies, the Vm response to the sinusoidal current stimulus (Fig. 8Ai) showed damped oscillations, and the maximal changes in Vm occurred over a narrow band of frequencies centred around 10 Hz. To quantify the electrical resonance, an impedance plot (Fig. 8Aii) was constructed by dividing the FFT of the