Conduction through the Inward Rectifier Potassium Channel, Kir2.1, Is Increased by Negatively Charged Extracellular Residues
http://www.100md.com
《普通生理学杂志》
1 Department of Physiology and Medicine and Department of Heart and Stroke, Richard Lewar Centre, University of Toronto, Toronto, Ontario M5S 3E2, Canada
2 Institute of Molecular Cardiobiology, Johns Hopkins University, Baltimore, MD 21218
3 Department of Chemistry and Department of Physics, University of Guelph, Guelph, Ontario N1G 2W1, Canada
ABSTRACT
Ion channel conductance can be influenced by electrostatic effects originating from fixed "surface" charges that are remote from the selectivity filter. To explore whether surface charges contribute to the conductance properties of Kir2.1 channels, unitary conductance was measured in cell-attached recordings of Chinese hamster ovary (CHO) cells transfected with Kir2.1 channels over a range of K+ activities (4.6–293.5 mM) using single-channel measurements as well as nonstationary fluctuation analysis for low K+ activities. K+ ion concentrations were shown to equilibrate across the cell membrane in our studies using the voltage-sensitive dye DiBAC4(5). The dependence of on the K+ activity (aK) was fit well by a modified Langmuir binding isotherm, with a nonzero intercept as aK approaches 0 mM, suggesting electrostatic surface charge effects. Following the addition of 100 mM N-methyl-D-glucamine (NMG+), a nonpermeant, nonblocking cation or following pretreatment with 50 mM trimethyloxonium (TMO), a carboxylic acid esterifying agent, the –aK relationship did not show nonzero intercepts, suggesting the presence of surface charges formed by glutamate or aspartate residues. Consistent with surface charges in Kir2.1 channels, the rates of current decay induced by Ba2+ block were slowed with the addition of NMG or TMO. Using a molecular model of Kir2.1 channels, three candidate negatively charged residues were identified near the extracellular mouth of the pore and mutated to cysteine (E125C, D152C, and E153C). E153C channels, but not E125C or D152C channels, showed hyperbolic –aK relationships going through the origin. Moreover, the addition of MTSES to restore the negative charges in E53C channels reestablished wild-type conductance properties. Our results demonstrate that E153 contributes to the conductance properties of Kir2.1 channels by acting as a surface charge.
Key Words: surface charges Kir2.1 conduction NMG TMO
N. D'Avanzo and H.C. Cho contributed equally to this project.
Abbreviations used in this paper: CHO, Chinese hamster ovarian; DiBAC4(5), bis-(1,3-dibutylbarbituric acid)pentamethine oxonol; GCS, Gouy-Chapman-Stern; NAChR, nicotinic acetylcholine receptor; NMG, N-methyl-D-glucamine; NSFA, nonstationary fluctuation analysis; TMO, trimethyloxonium; WT, wild-type.
INTRODUCTION
Electrostatic forces originating from charged groups on protein surfaces influence a variety of macromolecular functions. For example, charges on the surface of superoxide dismutases provides electrostatic guidance for substrates entering the active site (Getzoff et al., 1983, 1992; Sharp et al., 1987). Surface charges have also been shown to affect channel conductance in a variety of ion channels, such as neuronal Na+ channels (Sigworth and Spalding, 1980; Worley et al., 1986), Ca2+-activated K+ (BK) channels (MacKinnon and Miller, 1989; MacKinnon et al., 1989), and nicotinic acetylcholine receptors (NAChR) (Imoto et al., 1988; Konno et al., 1991), presumably by influencing the permeant ion concentration of ions at the entrance of the selectivity filter (Apell et al., 1977). In BK channels, surface charges also act to increase the affinity of cationic pore blockers by altering the local electrostatic potential at the pore entrance (MacKinnon and Miller, 1989; MacKinnon et al., 1989), while the "ring" of negative surface charges in NAChRs modifies the affinity of the pore blocker, Mg2+ (Imoto et al., 1988; Konno et al., 1991).
IK1 currents are found in many tissues and are generated by Kir2.x channels. In heart, a major molecular component of IK1 is the Kir2.1 channel (Zobel et al., 2003) and mutations in Kir2.1 channels have been linked to a form of long QT syndrome (LQT7) known as Andersen's syndrome (Plaster et al., 2001). Levels of IK1 in cardiomyocytes provide background K+ currents that are involved in setting the resting membrane potential, preventing membrane hyperpolarization due to Na+ pump activity, influencing propagation velocity, altering the electrical space constant, and promoting late phase repolarization. Consistent with an important role for IK1 in the heart, reductions in IK1 have been implicated in heart failure, and triggered arrhythmias resulting from early or delayed after-depolarizations observed in heart disease (Beuckelmann et al., 1993; Kb et al., 1996; Puglisi et al., 2000; Pogwizd et al., 2001). On the other hand, blockade of IK1 is a potent inhibitor of reentry type arrhythmias and ventricular fibrillation by increasing electrical space constants (Qi et al., 2002). IK1 channels also play a critical role in conditions of ischemia and hypoxia, where [K+]o is often elevated from 4 mM to as high as 20 mM (Weiss and Lamp, 1989), resulting in membrane potential depolarization as a consequence of the anomalous crossover properties of IK1 channels (Matsuda, 1991; Lopatin and Nichols, 1996).
The recent crystallography of the bacterial inward rectifier channel KcsA (Doyle et al., 1998) has provided important new insights into the molecular and biophysical mechanisms underlying the function of K+ channels. Indeed, crystallography and molecular dynamic models have revealed important channel properties such as ion occupancy and unitary conductance as well mechanisms for ion permeation (Berneche and Roux, 2001, 2003; Morais-Cabral et al., 2001). While previous reports of the versus K+ activity relationship in Kir2.1 channels were dedicated to studying the validity of various permeation models (Stampe et al., 1998), no studies have specifically examined the role of surface charges in Kir2.1 channel function.
In this study, we examined the conduction properties of the Kir2.1 channel over a range of K+ activities using single-channel recordings and nonstationary fluctuation analysis. We employed three experimental strategies to alter surface charges in Kir2.1 channels: variation of the bulk solution ionic strength, surface charge neutralization, and site-directed mutagenesis. Our results establish a critical role for surface charges near the extracellular entrance of the pore in Kir2.1 channels originating from E153 residues to the conduction properties of K+ ions. It will be necessary to consider the effect of E153 residues on ion occupancy and permeation in future studies.
MATERIALS AND METHODS
Expression of Mouse Kir2.1 Channels in Chinese Hamster Ovary (CHO) K1 Cells
Mouse Kir2.1 cDNA (a gift from L. Jan, Howard Hughes Medical Institute, San Francisco, CA) was subcloned into a bidirectional paBi-G vector (CLONTECH Laboratories, Inc.) modified by us to contain GFP (paBi-EGFP) (Zobel et al., 2003) in a TetOFF system requiring activation of the promoter by tetracycline transactivator (tTA) (CLONTECH Laboratories, Inc.). CHO-K1 cells were incubated at 37°C, 5% CO2 in F12 media (Sigma-Aldrich), supplemented with 5% FBS and 1% penicillin/streptomycin (GIBCO BRL) and transfected with 2 μg of Kir2.1 cDNA plus 1 μg of tTA cDNA using 4 μl of Lipofectamine 2000 (GIBCO BRL) following the supplied protocol. Transfected cells were incubated in supplemented culture media (see below) for 24 h before the electrophysiological recordings. Single mutants were created on a four cysteine knockout background (4CKO = C89S/C101S/C149Y/C169S) via the Kunkel method.
Cell-attached Macropatch Recordings for Unitary Conductance
Electrophysiological recordings were done using the cell-attached patch clamp technique. Bath solutions contained (in mM) 2 MgCl2, 10 HEPES, and 5–400 K+ (Cl– and OH–) (pH 7.4, with KOH). The osmolarity was maintained when necessary by the addition of varying amounts of glucose. Pipette solutions had the same composition as the bath solutions, except MgCl2 was removed and 1–30 μM Ba2+ was added in order to induce time-dependent current decay, allowing for nonstationary fluctuation analysis to be used for unitary current estimates (see below).
For some experiments, recording solutions also contained 100 mM N-methyl-D-glucamine chloride (NMG+Cl–) to "screen" charges. In some studies, transfected cells were pretreated in our culture medium (see below) for 5 min with 50 mM trimethyloxonium (TMO), an agent used to esterify carboxylic acid residues of negatively charged amino acid side chains (Sigworth and Spalding, 1980; Dudley and Baumgarten, 1993). To replace negative charges in some studies using mutant channels, 1 mM MTSES (in DMSO) was added for 10–15 min and then thoroughly washed out with fresh bath solution before recording. TW140F-4 patch pipettes (World Precision Instruments) were pulled and fire polished to give resistances of 6–10 M for single channel recordings and resistances 2 M for fluctuation analysis recordings. Pipettes were coated with Sylgard (Dow Corning Corporation) to reduce their capacitance. Data were acquired using Axopatch 200A amplifier, collected at 20 kHz by Clampex 8 (pCLAMP; Axon Instruments) with a 2 kHz low-pass Bessel filter and analyzed using Clampfit 8 (pCLAMP; Axon Instruments).
To obtain estimates for the unitary conductance () at varied levels of symmetric K+ in our cell attached recordings, it is desirable for the K+ concentrations to equilibrate across the cell membrane and for the resting membrane potential to be 0 mV. For thermodynamic reasons, this cannot be readily established directly from our cell-attached recordings. Therefore, to demonstrate K+ ion equilibration, we used the voltage-sensitive dye bis-(1,3-dibutylbarbituric acid)pentamethine oxonol (DiBAC4(5)) (Molecular Probes), which increases its fluorescence as the membrane potential depolarizes. To perform these studies, DiBAC4(5) was dissolved in 0.2% DMSO and added to bath solutions at a final concentration of 200 nM. Initially transfected cells were incubated for 20 min with a solution containing 140 mM KCl, 10 mM HEPES, 2 mM MgCl2, 200 nM DiBAC4(5), pH 7.4. Next, the fluorescence (at 616 nm) generated by DiBAC4(5) illuminated at 590 nm was recorded every 3 s from cells expressing Kir2.1 channels, identified by GFP fluorescence. Fig. 1 shows typical fluorescence recordings from DiBAC4(5) in response to changes in K+ levels from 140 mM to 20 mM K+ (20 mM K+, 10 mM HEPES, 220 mM glucose, 2 mM MgCl2, 200 nM DiBAC4(5), pH 7.4). As expected from the properties of DiBAC4(5), the application of 20 mM K+ results in rapid fluorescence decreases, indicative of membrane hyperpolarization, which was followed by a slow return of the fluorescence to baseline, recorded when 140 mM K+ was present. Furthermore, reintroduction of a 140 mM K+ containing solution, caused a rapid increase in fluorescence, indicating membrane depolarization followed again by decay back to baseline. Identical patterns of DiBAC4(5) fluorescence changes were observed when other solution changes were applied. In general, these observations can only be obtained if K+ ions equilibrate and have equal concentrations across the cell membrane. This conclusion is bolstered by the reasonable assumption that the initial membrane potential in our 140 mM K+ containing solution was 0 mV (since intracellular environment of our cells should be typical of other mammalian cells that contain 140 mM K+), thereby indicating that the baseline fluorescence in our recordings corresponds to a membrane potential of 0 mV.
The observation that the membrane potential equilibrates to 0 mV and the K+ levels are equal across the cell membrane (regardless of the applied external K+ level) is entirely expected for the following reasons. First, the conductance of these cells following transfection is very high (i.e., >40 nS) and K+ selective, which promotes K+ equilibration across the membrane. Second, since Na+ was excluded from the media, thereby inhibiting the Na+ pump, there is no biochemical mechanism (that we are aware of) to maintain K+ gradients across the membrane.
Nonstationary Fluctuation Analysis (NSFA) on Cell-attached Multichannel Recordings
Unitary currents were estimated by the nonstationary fluctuation analysis (NSFA) method using the equation (Sigworth, 1980)
(1)
where 2 is the variance of the measured macroscopic current, I, is the mean of I, N is the number of channels in the patch, and L is the contribution of leak currents. Since Kir2.1 channels do not vary greatly in response to voltage steps as is required for NSFA, Ba2+ was included in the extracellular solutions to induce time- and voltage-dependent block of Kir2.1 currents (Neyton and Miller, 1988; Jiang and MacKinnon, 2000). In the presence of Ba2+, 200–300 successive K+ current sweeps were recorded in response to voltage steps from a holding potential of +5 mV to –60, –80, –100, –120, and –140 mV. Variance and mean currents were estimated using consecutive pairs of traces (to minimize the complicating effects of current drift) and these estimates for 2 and were averaged. Estimates for 2 versus were fit with Eq. 1 using Maple 7.00 (Waterloo Maple Inc.) to estimate i. The single-channel conductance () was calculated from the dependence of unitary current estimates on voltage using Ohm's law (Kubo et al., 1993; Fig. 2 B).
Unitary Conductance–K+ Activity Relationship
To properly account for the effect of counter-ions in the solution on the driving force for K+ ion entry into the pore (i.e., the chemical potential), the activity of K+ (aK) was calculated using the relationship
(2)
To accurately estimate the contribution of long-range electrostatic interactions to the activity coefficients, fKact, in our studies we used the following equations. When the ionic strength (I) was <10 mM (i.e., I < 10 mM), fKact was calculated using the relationship (Debye and Hückel, 1923)
(3)
where z+ is the valence of the major cation, z– the major anion, and I is the ionic strength of the solution, which is defined as (Debye and Hückel, 1923)
(4)
For 10 mM < I < 100 mM, the Extended Debye-Hückel relationship was used (Kielland, 1937; MacInnes, 1961):
(5)
where B 0.33 in water at 25°C and a is an adjustable parameter corresponding to ionic size (=3 x 10–8 cm for K+). For I > 100 mM, we used Davies' approximation (MacInnes, 1961),
(6)
The dependence of on aK (i.e., the –aK relationship) was initially fit with a modified Langmuir binding isotherm (Fig. 3 B):
(7)
where aK is defined by Eq. 3, 5, or 6 as required, max is the maximal conductance measured at saturating aK, and EC50 is a measure of aK when is max /2. The parameter S is required to account for the nonzero intercept of the –aK relationship which is predicted to originate from the presence of surface charges (see below).
Surface Charge Calculations
To quantify the contribution of surface charges to Kir2.1 channel properties, we assumed that the conductance is a function of the activity of K+ ions at the channel surface (aKs) and not aK in the bulk solution. aKs was estimated by the relationship
(8)
where z, e, kB, and T have their usual meanings, and s is the electrostatic potential at the surface. In our studies, we estimated s using the relationship (McLaughlin et al., 1970)
(9)
where control is the channel conductance in the presence of surface charges, and TMO is the unitary conductance in following pretreatment with TMO that was assumed to neutralize all negative surface charges. Evidence for the validity of this assumption is provided in the results section.
To estimate the number of charged residues that might contribute to the observed surface charge effects, the surface charge density () was calculated using the following equation (Grahame, 1947; McLaughlin, 1977):
(10)
where aa = 0.15 M–1 (Eisenberg et al., 1979), aK is the activity in the bulk solution, r is assumed to be the dielectric of water (80), and o is the permittivity of free space. Values of for every s and aKs were averaged (±SEM) to determine the final estimate.
Kinetics of Barium Block
Since Ba2+ blocks Kir2.1 channels by binding within the selectivity filter, time constants () for Kir2.1 current decay in the presence of Ba2+ are predicted to be influenced by surface charges in a manner similar to K+ conductance (see above). By assuming that Ba2+ block is adequately described by a single binding site model (Neyton and Miller, 1988; Alagem et al., 2001) where only the on-rate is influenced by the surface potential combined with the assumptions that the rate of Ba2+ entry into the channel is given by the Ba2+ activity at the channel surface (aBas) and that Ba2+ does not alter the surface potential at the low concentrations used (see Eqs. 3 and 4), we can write
(11)
where kon(V) and koff(V) are the rates of Ba2+ binding and unbinding, respectively, at a given voltage, and aBas is the external Ba2+ activity at the surface of the channel, which is related to Ba2+ activity in the bulk (aBa) as described above for K+ ions in Eq. 6. Assuming that kon(V) depends exponentially on voltage (Reuveny et al., 1996), it can be shown that
(12)
where (1/) is the difference in 1/ for a given [Ba2+] in the presence and absence of NMG, and on is an estimate of the electrical distance for the position of the peak of the free energy barrier for Ba2+ binding to the channel pore. kon(0) and on were obtained from the intercept and slope of the relationship between 1/ and voltage (–80, –100, and –120 mV) in the presence of NMG (i.e., where the surface charges are screened) (Reuveny et al., 1996). Estimates of s using this method were made at each voltage and averaged in order to determine Ba.
Data Analysis
With the exception of NSFA, all other data were analyzed using Origin 6.0 (Microcal Software Inc.). This includes mathematical calculations, graph development, curve fitting with the appropriate functions, as well as statistical analysis.
RESULTS
Kir2.1 Single-channel Conductances at Varying K+ Concentrations
To investigate the existence of surface charges in Kir2.1 channels, we examined the conductance as a function of K+ activity. Fig. 2 A illustrates typical single channel currents (i) recorded in CHO cells transfected with Kir2.1 channels at several K+ concentrations. As shown in Fig. 2 B for 20 and 140 mM K+, i depends linearly on voltage (i.e., the channels are ohmic) over the range of voltages studied as previously reported (Kubo et al., 1993). Thus, unitary conductances () were estimated using Ohm's Law. Since the unitary current could not be reliably measured using single-channel recordings with <20 mM K+ for Kir2.1 channels, NSFA of currents recorded in cell-attached patches containing Ba2+ were used to estimate i. Typical cell-attached recordings are shown in Fig. 3 A, along with a typical plot of the mean variance (<2>) versus mean current () assessed from analysis of 200 sweeps, which allowed estimates of i from the parabolic function in Eq. 1. Fig. 3 B plots the estimated unitary conductance from single channel recordings (closed squares) and NSFA (open squares) as a function of the bulk K+ activity. Note that the estimates using NSFA were identical (P < 0.01) to the values estimated using single-channel recordings (at 20 and 140 mM K+). These data were fit (R2 = 0.98) using the modified Langmuir binding isotherm (see MATERIALS AND METHODS) to estimate max (26.0 ± 1.5 pS, 10 < n < 20) and EC50 (50 ± 15 mM aK) as well as the intercept, S (4.4 ± 1.4 pS). These results establish a clear saturation of unitary conductance at high K+ concentrations as reported previously (Stampe et al., 1998) and suggest the existence of a surface charge.
To ensure that our conductance estimates did not arise from experimental limitations of the NSFA method, we performed Monte Carlo simulations of channel blocking events (Hanna et al., 1996) using experimentally derived mean blocking and unblocking times in the presence of 30 μM Ba2+ (Table I) assuming a single site blocking model, as established previously (Neyton and Miller, 1988; Alagem et al., 2001). Typical simulated single-channel events are shown in Fig. 4 A (bottom two traces), while Fig. 4 B shows the predicted macroscopic current generated by superimposing 100 sweeps, as well as the predicted dependence of <2> on . Clearly, the simulated results in Fig. 4 B (left) closely resemble the experimental results in Fig. 3 A. More important, the estimated unitary current estimated from the <2> versus plot is identical to the assumed unitary current (Makhina et al., 1994, Krapivinsky et al., 1998).
Charge Screening Decreases Kir2.1 Single Channel Conductance at a Low K+
The presence of the nonzero intercept, S, in the versus aK relationship suggests the presence of a surface charge. To investigate this possibility, the versus aK relationship was measured in the presence of 100 mM NMG to screen surface charges. As shown in Fig. 5 B, with 200 mM symmetrical K+ (see MATERIALS AND METHODS), the inclusion of 100 mM NMG did not affect the Kir2.1 single channel conductance (solid line is without NMG), confirming that NMG does not block Kir2.1 channels. By contrast, at 5 mM K+, the single channel conductance was reduced (P < 0.01) to 2.0 ± 0.3 pS in the presence of 100 mM NMG compared with 6.6 ± 1.5 pS in the absence of NMG (Fig. 5 A). Reductions in were also observed at 10, 15, 20, and 40 mM K+, when 100 mM NMG was added (P < 0.05). More important, fitting the –K+ activity relationship with Eq. 7, when 100 mM NMG was present, reduced the conductance intercept to a value not different from zero (i.e., 0.8 ± 1.3, P = 0.36) (Fig. 5 C).
As an alternate method for eliminating surface charges, CHO cells were treated with TMO, a carboxylate esterifying agent. The unitary conductance with TMO was 2.4 ± 1.0 pS with 5 mM K+ which was below (P < 0.01) the value of estimated without TMO treatment (Fig. 5 A). Furthermore, the unitary conductance in 5 mM K+ after treatment with TMO was not reduced further by the inclusion of 100 mM NMG (2.9 ± 1.0 pS) (Fig. 5 A), and the intercept for the versus aK relationship (0.8 ± 1.9 pS) was not different (P = 0.32) from zero. Taken together, these results suggest that negatively charged amino acid residues generate surface charges in Kir2.1 channels, thereby increasing unitary currents at low permeant ion concentrations.
Electrostatic Theory Predictions of Surface Charge Effects
Using our results above, it is possible to calculate the surface charge potential (s) for Kir2.1 channels (Eq. 9) and thereby the activity of K+ ions at the mouth of the channel pore (i.e., aKs) (using Eqs. 8 and 9). Fig. 6 shows the relationship between and aKS (closed squares), which can be well fit (R2 = 0.99) with Eq. 7 going through the origin (S = 0.7 ± 1.5 pS, P > 0.05). Notice that the relationship of versus aKs in control conditions is identical to the versus aK (i.e., bulk K+ activity) relationship for channels in the presence of NMG (open triangles), consistent with a role for surface charges in concentrating permeant K+ ions at the channel surface.
Our estimates of s were also used to calculate the apparent surface charge density () of Kir2.1 channels (Eq. 10). Values of charge density in Kir2.1 channels were calculated for each aKs examined assuming a dielectric (r) of 80 (that of water), and then averaged to determine a of –0.49 ± 0.08 e nm–2. When NMG was present, the estimated surface charge was reduced to –0.06 ± 0.02 e nm–2, which represents less than one net charge per channel subunit. This is consistent with the screening of surface charges in Kir2.1 channels by NMG. Using current molecular models of Kir2.1 channels (Cho et al., 2000), we estimated the channel's surface area to be 11 nm2, which predicted the presence of approximately –5.4 e net charges per channel. Since Kir2.1 channels are tetramers, these calculations suggest that the surface charge arises from one free net charge per monomer.
Kinetics of Extracellular Barium Block Are Altered by Charge Screening
Ba2+ is a potent voltage- and time-dependent blocker of K+ current in Kir2.1, by binding within the selectivity filter (Jiang and MacKinnon, 2000). Thus, the presence of surface charges is predicted to alter the blocking time constant () of Ba2+. Indeed, as shown in Fig. 7 A, treatment with NMG or TMO significantly increased at 5 mM K+ compared with control. The same was observed for 20 mM (Fig. 7, compare B and C) and 140 mM [K+]. Using the voltage dependence of kon in the absence of surface charges (Fig. 7 D) and Eqs. 11 and 12, s was calculated using Ba2+ blocking kinetics for 20 mM K+, which again allows an estimate of the charge density (Ba) in Kir2.1. Ba was estimated to be –0.49 ± 0.12 e nm–2, which is indistinguishable (P < 0.01) from the value obtained using the conductance data, confirming the presence of four negative charges on the surface of Kir2.1 channels.
Identification of Extracellular Charges Acting as Surface Charges
Using the molecular model for Kir2.1 channels, which is based on the KcsA crystal structure (Cho et al., 2000), we identified three candidate charged residues (E125, D152, and E153) that might act as surface charges (Fig. 8 A). We replaced these charged residues with cysteine in order to examine the potential role of these charges in the electrostatic effects of Kir2.1 channels. Cysteine was chosen to allow charge replacement using MTS compounds. The dependence of on aK was unaltered for E125C and D152C mutant channels, compared with wild-type (WT) channels at 20 mM K+ (Fig. 8 B). On the other hand, E153C channels showed a smaller unitary conductance compared with WT, which was not different from the unitary conductance observed in WT channels pretreated with TMO (Fig. 8 B) (P < 0.01). Fig. 8 C shows that the same results were observed for 5 mM K+, whereas, at 140 mM K+, WT channels, with and without TMO treatment, and all three mutant channels had similar values. Treatment of E153C channels with MTSES, in order to replace the negative charge, increased the to levels similar to WT at 20 mM K+ (Fig. 8 B). Taken together, these results support the conclusion that E153 acts as the primary source of surface charges in Kir2.1 channels.
DISCUSSION
Our results establish that Kir2.1 channels display classical surface charge properties (Sigworth and Spalding, 1980; Worley et al., 1986; MacKinnon and Miller, 1989; MacKinnon et al., 1989) that depend on E153 residues predicted from a previous molecular model to be located on the outer face of the channel (Cho et al., 2000). Investigations into the possible role of surface charges in our studies necessitated the conversion of K+ ion concentrations into chemical activity (i.e., chemical potential) in order to account for the effects of ionic strength on the thermodynamic driving force of K+ ions into Kir2.1 channels. K+ activity was calculated using a combination of the linear Debye-Huckel, the Extended Debye-Huckel, and the Davies' theories to account for long-range electrostatic interactions between ions in solution under our experimental conditions (Debye and Hückel, 1923; MacInnes, 1961). The use of multiple equations was necessary because no single analytical expression can comprehensively estimate ionic activities over the range of experimental conditions used in our studies (MacInnes, 1961). The use of K+ activity was essential in our studies since a wide range of ionic strengths were used to alter surface potentials in Kir2.1 channels.
In our studies, it was necessary to determine the unitary current and unitary conductance () with varied symmetric levels of K+. This could not be achieved using whole-cell recordings since the current magnitudes were too large to allow NSFA to be effectively used (see below). Therefore, to circumvent this problem, we resorted to cell-attached recordings. To ensure the levels of K+ were symmetric in our cell-attached recordings, we monitored the membrane potential using a voltage-sensitive dye. These studies showed that K+ levels equilibrate within 2–3 min following changes of the external K+ levels under our experimental conditions, probably due to the presence of very high levels of Kir2.1 channel expression (i.e., >40 nS in whole-cell recordings with 140 mM K+) combined with the omission of Na+, thereby inhibiting the Na+ pump.
Clear evidence for the presence of surface charges was obtained in our studies using two methods that yielded comparable estimates for the electrostatic surface potential in Kir2.1 channels. In the first method, the dependence of unitary current () on the K+ activity was used. These studies required the use of NSFA at low K+ activities, due to small unitary currents, as done previously in HRK1 (Makhina et al., 1994) and Kir7.1 (Krapivinsky et al., 1998) channels. We validated the NSFA method in our studies using Monte Carlo simulations to model Ba2+ blocking events. In the versus aK plots, we found that does not approach zero as K+ activities are reduced. While it is possible that the finite limiting intercept originates from finite permeability of Kir2.1 channels to other cations in our solutions (e.g., Ba2+), Kir2.1 channels are impermeant to these other ions (Alagem et al., 2001). More important, the limiting was reduced to zero at low aK following the addition of the impermanent cations NMG+ to screen surface charges or following pretreatment of the channels with TMO to specifically neutralize carboxylate residues (Parsons et al., 1969; Nakayama et al., 1970). The nonzero conductance intercept at low K+ activity, combined with a zero conductance intercept at high ionic strength or following TMO pretreatment are signature features of local surface charges (Manning, 1978; Fixman, 1979; Klein et al., 1981; Cai and Jordan, 1990).
To analyze and interpret our versus aK results we used the Gouy-Chapman-Stern (GCS) theory for surface charges, which assumes the surface charge density is uniformly distributed in a flat plane. While this assumption is generally expected to be invalid for ion channel proteins, the GCS theory accurately accounts for surface charge effects in many proteins, possibly due to cancellation of errors in the functional expansions of rigorously correct canonical partition functions (Fixman, 1979; Klein et al., 1981). We anticipate that the applicability of the GCS theory is assisted by the symmetrical arrangement of four negative point charges around the pore opening of tetrameric Kir2.1 channels that should mimic more closely a uniform charge distribution than a single point charge located near the pore opening. Regardless, the application of the GCS theory to our results predicted that depends uniquely on the local K+ ion activity at the entrance to the Kir2.1 channel pore. Indeed, the dependence of on aK at the channel surface (i.e., aKs) could be accurately described by a Langmuir binding isotherm going through the origin. The GCS theory further predicted a surface charge density of –0.49 ± 0.08 e nm–2. Using recent molecular models of Kir2.1 channels (Cho et al., 2000) to estimate channel area, we calculated a net of approximately four charges on the surface of Kir2.1 channels in the absence of charge screening agents, which equated to nearly one charge per monomer.
Thus, our results strongly support the conclusion that surface charges influence the Kir2.1 channel conductance. Similar investigations and analyses have been applied to other channels, including Na+ channels (Sigworth and Spalding, 1980; Worley et al., 1986), Ca2+-activated K+ channel (Moczydlowski et al., 1985; MacKinnon and Miller, 1989; MacKinnon et al., 1989; Park et al., 2003), NAChRs (Imoto et al., 1988; Konno et al., 1991), to establish the presence of surface charges. On the other hand, versus concentration in KcsA channels show two asymptotic conductance values at low and high K+ levels (Morais-Cabral et al., 2001; Nimigean et al., 2003). The limiting value of KcsA at low K+ levels appears remarkably similar to that observed in Kir2.1 channels, suggesting that surface charges may also play an important role in the conductance properties of these channels as well.
Since Ba2+ blocks Kir2.1 channels by binding within the channel pore, Ba2+ should be affected by the presence of surface charges in a comparable manner to K+ ions. Specifically, the rate of Ba2+ entry into the channel should be proportional to the activity of Ba2+ at the channel entry. This provided a second method for estimating surface charges in Kir2.1. Indeed, when charges are screened by NMG, or neutralized by TMO, the kinetics of Ba2+ block is slowed significantly. Furthermore, after accounting for ionic strength (i.e., activity) of Ba2+ in our recording conditions, our estimates of the surface charge density from our Ba2+ block studies (Ba), using the GCS theory, were identical to those estimated from our conductance studies, once again supporting the existence of four negative charges on the surface of Kir2.1 channel proteins. This contrasts with a previous report examining the effects of lipid charges on BK channels (Park et al., 2003), where alterations in lipid charge (i.e., lipid surface charges) affected BK channel conductance, but not Ba2+ block properties. This led the authors to conclude, consistent with molecular models of channel proteins, that lipid charges are too remote from the channel pore entrance to produce classical surface charge effects (Park et al., 2003). While we did not examine the effects of altered lipid composition on the Kir2.1 channel properties, our results are consistent with the conclusions of Park et al. (2003) since surface charge effects in Kir2.1 channels originate from charged E153 (see below), predicted to be more adjacent to the channel pore entrance than charged lipids. Further studies will be required to examine whether lipid charge also affects Kir2.1 channel properties.
More definitive evidence for the existence of surface charges in Kir2.1 channels was provided by our mutagenesis studies. Somewhat surprisingly, E125C and D152C mutant channels showed conductance properties identical to WT channels. Previous studies have suggested that E125 contributes a superficial site for binding of extracellular divalent cations such as Mg2+, Ba2+, and Ca2+ (Murata et al., 2002). Indeed, in the presence of divalents, E125Q mutant channels had decreased single-channel conductance compared with WT channels. Furthermore, the E125 was also postulated to play a key role in dehydration of Ba2+ before entering the pore (Alagem et al., 2001). Consistent with these results, Ba2+ blocking kinetics were increased in E125C mutant channels compared with WT channels. The significance of these observations is unclear but could arise from altered binding affinities for these divalents at this site.
Our observations in D152C channels were unexpected, based on our model of Kir2.1 channels (Cho et al., 2000). This may indicate that the D152 residue binds K+ ions with high affinity, thereby not contributing to the surface charge effect. Alternatively, our molecular model (Cho et al., 2000) may require revision with the D152 residue either buried or oriented in such a way that its side chain does not contribute to the surface potential. Although less likely, it is conceivable that the sulfhydryl side chains of the C125 and C152 residues may be ionized in the E125C and D152C mutant channels. Consequently, the residues would not have led to charge neutralization, thereby having maintained their surface potential. However, this possibility would be inconsistent with the expectation that polar residues are less likely (not more likely) to be ionized (i.e., have a higher pKa) when localized in regions with reduced dielectric constants, such as in proteins compared with the aqueous phase.
Unlike E125C and D152C mutant channels, E153C channels have reduced single-channel conductance compared with WT channels at low ionic activities and are identical in conductance properties to WT channels pretreated with 50 mM TMO (Fig. 8 B). Moreover, addition of 1 mM MTSES to restore negative charges restored WT behavior to E153C channels (Fig. 8 B). These results support the conclusion that E153 is responsible for the surface charge effects observed in Kir2.1 channels by contributing approximately one negative charge per monomer at the channel surface. In accordance with electrostatic theory, our results support the conclusion that ionized carboxyl groups increase the conduction at low permeant ion activities, by attracting the permeant ions near the conduction pathway.
Our results also establish that for Kir2.1 channels, increases with elevations in aK, exhibiting a clear maximum conductance >160 mM (K+ concentrations >300 mM) with an apparent aK EC50 of 50 ± 15 mM. These properties of Kir2.1 channels resemble those of ROMK1 (Kir1.1) channels that have a similar EC50, although ROMK1 channels exhibit a maximum conductance as aK is elevated above 300 mM K+ (Heginbotham and MacKinnon, 1993; Lu and MacKinnon, 1994). It is conceivable that of Kir2.1 channels might also reach a maximum at higher K+ concentrations. In stark contrast to Kir1.1 and Kir2.1 channels, the –[K+] relationship in KcsA channels is nearly linear from 50 to 500 mM K+ (Morais-Cabral et al., 2001), with evidence of saturation at K+ concentrations >1600 mM K+ (Nimigean et al., 2003). The basis for the differences in the properties of various inward rectifier channels will clearly require further investigations.
While our studies focused on the contribution of surface charges from E153 residues to Kir2.1 channel conductance, it is conceivable that surface charges may affect other channel properties. Indeed, since polyamine and Mg2+ block Kir2.1 channels by binding within the channel pore, alterations in surface charges could conceivably alter binding stability of these blockers. Thus, surface charge effects may contribute to the crossover effect observed in Kir2.1 channels with increased extracellular [K+].
Our findings support the presence of approximately four negative charges on the surface of Kir2.1 channels, originating from E153 (but not E125 or D152) residues. The results of our studies should help guide the development of new molecular models of Kir2.1 channels in the absence of crystal structures. It is apparent that new models of Kir2.1 channels must reexamine the location, orientation, and ionization of these three residues to account for the surface charge effect we attribute to the E153 residue. As well, future studies of ion occupancy, permeation, and perhaps rectification in Kir2.1 channels need to consider surface charge effects as described in our studies.
ACKNOWLEDGMENTS
We thank L.Y. Jan for providing us with the Kir2.1 clone and A. Prendergast for his technical assistance.
This study was supported by Heart and Stroke Foundation of Ontario (grant T-4179) to P.H. Backx. P.H. Backx is a Career Investigator of the Heart and Stroke Foundation of Ontario. H.C. Cho was supported by a Canadian Institutes of Health Research studentship. N. D'Avanzo holds a studentship from Heart & Stroke Richard Lewar Center of Excellence at the University of Toronto.
David C. Gadsby served as editor.
Submitted: 25 August 2004
Accepted: 18 March 2005
REFERENCES
Alagem, N., M. Dvir, and E. Reuveny. 2001. Mechanism of Ba2+ block of a mouse inwardly rectifying K+ channel: differential contribution by two discrete residues. J. Physiol. 534:381–393.
Apell, H.J., E. Bamberg, H. Alpes, and P. Lauger. 1977. Formation of ion channels by a negatively charged analog of gramicidin A. J. Membr. Biol. 31:171–188.
Berneche, S., and B. Roux. 2001. Energetics of ion conduction through the K+ channel. Nature. 414:73–77.
Berneche, S., and B. Roux. 2003. A microscopic view of ion conduction through the K+ channel. Proc. Natl. Acad. Sci. USA. 100:8644–8648.
Beuckelmann, D.J., M. Nabauer, and E. Erdmann. 1993. Alterations of K+ currents in isolated human ventricular myocytes from patients with terminal heart failure. Circ. Res. 73:379–385.
Cai, M., and P.C. Jordan. 1990. How does vestibule surface charge affect ion conduction and toxin binding in a sodium channel. Biophys. J. 57:883–891.
Cho, H.C., R.G. Tsushima, T. Nguyen, H.R. Guy, and P.H. Backx. 2000. Two critical cysteine residues implicated in disulfide bond formation and proper folding of Kir2.1. Biochemistry. 39:4649–4657.
Debye, P., and E. Hückel. 1923. Zur Theorie der Elektrolyte. II. Das Grenzgesetz fur die elektrische Leitfahigkeit. Physiol. Z. 24:305–325.
Doyle, D.A., J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, and R. MacKinnon. 1998. The structure of the potassium channel: molecular basis of K+ conduction and selectivity. Science. 280:69–77.
Dudley, S.C., Jr., and C.M. Baumgarten. 1993. Modification of cardiac sodium channels by carboxyl reagents. Trimethyloxonium and water-soluble carbodiimide. J. Gen. Physiol. 101:651–671.
Eisenberg, M., T. Gresalfi, T. Riccio, and S. McLaughlin. 1979. Adsorption of monovalent cations to bilayer membranes containing negative phospholipids. Biochemistry. 18:5213–5223.
Fixman, M. 1979. The Poisson-Boltzmann equation and its application to polyelectrolytes. J. Chem. Physiol. 70:4995–5005.
Getzoff, E.D., D.E. Cabelli, C.L. Fisher, H.E. Parge, M.S. Viezzoli, L. Banci, and R.A. Hallewell. 1992. Faster superoxide dismutase mutants designed by enhancing electrostatic guidance. Nature. 358:347–351.
Getzoff, E.D., J.A. Tainer, P.K. Weiner, P.A. Kollman, J.S. Richardson, and D.C. Richardson. 1983. Electrostatic recognition between superoxide and copper, zinc superoxide dismutase. Nature. 306:287–290.
Grahame, D.C. 1947. The electrical double layer and the theory of electrocapiliarity. Chem. Rev. 41:441–501.
Hanna, W.J., R.G. Tsushima, R. Sah, L.J. McCutcheon, E. Marban, and P.H. Backx. 1996. The equine periodic paralysis Na+ channel mutation alters molecular transitions between the open and inactivated states. J. Physiol. 497(Pt 2):349–364.
Heginbotham, L., and R. MacKinnon. 1993. Conduction properties of the cloned Shaker K+ channel. Biophys. J. 65:2089–2096.
Imoto, K., C. Busch, B. Sakmann, M. Mishina, T. Konno, J. Nakai, H. Bujo, Y. Mori, K. Fukuda, and S. Numa. 1988. Rings of negatively charged amino acids determine the acetylcholine receptor channel conductance. Nature. 335:645–648.
Jiang, Y., and R. MacKinnon. 2000. The barium site in a potassium channel by x-ray crystallography. J. Gen. Physiol. 115:269–272.
Kb, S., H.B. Nuss, N. Chiamvimonvat, B. O'Rourke, P.H. Pak, D.A. Kass, E. Marban, and G.F. Tomaselli. 1996. Ionic mechanism of action potential prolongation in ventricular myocytes from dogs with pacing-induced heart failure. Circ. Res. 78:262–273.
Kielland, J. 1937. Individual activity coefficients of ions in aqueous solutions. J. Am. Chem. Soc. 59:1675–1678.
Klein, B.K., C.F. Anderson, and M.T. Record. 1981. Condensation model expressions for the colligative properties of cylindrical polyions. Biopolymers. 20:2263–2280.
Konno, T., C. Busch, E. Von Kitzing, K. Imoto, F. Wang, J. Nakai, M. Mishina, S. Numa, and B. Sakmann. 1991. Rings of anionic amino acids as structural determinants of ion selectivity in the acetylcholine receptor channel. Proc. R. Soc. Lond. B. Biol. Sci. 244:69–79.
Krapivinsky, G., I. Medina, L. Eng, L. Krapivinsky, Y. Yang, and D.E. Clapham. 1998. A novel inward rectifier K+ channel with unique pore properties. Neuron. 20:995–1005.
Kubo, Y., T.J. Baldwin, Y.N. Jan, and L.Y. Jan. 1993. Prmary structure and functional expression of a mouse inward rectifier potassium channel. Nature. 362:127–133.
Lopatin, A.N., and C.G. Nichols. 1996. [K+] dependence of polyamine-induced rectification in inward rectifier potassium channels (IRK1, Kir2.1). J. Gen. Physiol. 108:105–113.
Lu, Z., and R. MacKinnon. 1994. A conductance maximum observed in an inward-rectifier potassium channel. J. Gen. Physiol. 104:477–486.
MacInnes, D.A. 1961. The Principles of Electrochemistry, New edition. Dover, New York. 478 pp.
MacKinnon, R., and C. Miller. 1989. Functional modification of a Ca2+-activated K+ channel by trimethyloxonium. Biochemistry. 28:8087–8092.
MacKinnon, R., R. Latorre, and C. Miller. 1989. Role of surface electrostatics in the operation of a high-conductance Ca2+-activated K+ channel. Biochemistry. 28:8092–8099.
Makhina, E.N., A.J. Kelly, A.N. Lopatin, R.W. Mercer, and C.G. Nichols. 1994. Cloning and expression of a novel human brain inward rectifier potassium channel. J. Biol. Chem. 269:20468–20474.
Manning, G.S. 1978. The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides. Q. Rev. Biophys. 11:179–246.
Matsuda, H. 1991. Effects of external and internal K+ ions on magnesium block of inwardly rectifying K+ channels in guinea-pig heart cells. J. Physiol. 435:83–99.
McLaughlin, S. 1977. Electrostatic potentials at membrane-solution interfaces. Curr. Top. Membr. Transp. 9:71–144.
McLaughlin, S.G., G. Szabo, G. Eisenman, and S.M. Ciani. 1970. Surface charge and the conductance of phospholipid membranes. Proc. Natl. Acad. Sci. USA. 67:1268–1275.
Moczydlowski, E., O. Alvarez, C. Vergara, and R. Latorre. 1985. Effect of phospholipid surface charge on the conductance and gating of a Ca2+-activated K+ channel in planar lipid bilayers. J. Membr. Biol. 83:273–282.
Morais-Cabral, J.H., Y. Zhou, and R. MacKinnon. 2001. Energetic optimization of ion conduction rate by the K+ selectivity filter. Nature. 414:37–42.
Murata, Y., Y. Fujiwara, and Y. Kubo. 2002. Identification of a site involved in the block by extracellular Mg2+ and Ba2+ as well as permeation of K+ in the Kir2.1 K+ channel. J. Physiol. 544:665–677.
Nakayama, H., K. Tanizawa, and Y. Kanaoka. 1970. Modification of carboxyl groups in the binding sites of trypsin with the Meerwein reagent. Biochem. Biophys. Res. Commun. 40:537–541.
Neyton, J., and C. Miller. 1988. Potassium blocks barium permeation through a calcium-activated potassium channel. J. Gen. Physiol. 92:549–567.
Nimigean, C.M., J.S. Chappie, and C. Miller. 2003. Electrostatic tuning of ion conductance in potassium channels. Biochemistry. 42:9263–9268.
Park, J.B., H.J. Kim, P.D. Ryu, and E. Moczydlowski. 2003. Effect of phosphatidylserine on unitary conductance and Ba2+ block of the BK Ca2+-activated K+ channel: re-examination of the surface charge hypothesis. J. Gen. Physiol. 121:375–397.
Parsons, S.M., L. Jao, F.W. Dahlquist, C.L. Borders Jr., J. Racs, T. Groff, and M.A. Raftery. 1969. The nature of amino acid side chains which are critical for the activity of lysozyme. Biochemistry. 8:700–712.
Plaster, N.M., R. Tawil, M. Tristani-Firouzi, S. Canun, S. Bendahhou, A. Tsunoda, M.R. Donaldson, S.T. Iannaccone, E. Brunt, R. Barohn, et al. 2001. Mutations in Kir2.1 cause the developmental and episodic electrical phenotypes of Andersen's syndrome. Cell. 105:511–519.
Pogwizd, S.M., K. Schlotthauer, L. Li, W. Yuan, and D.M. Bers. 2001. Arrhythmogenesis and contractile dysfunction in heart failure: roles of sodium-calcium exchange, inward rectifier potassium current, and residual -adrenergic responsiveness. Circ. Res. 88:1159–1167.
Puglisi, J.L., S.M. Pogwizd, W. Yuan, and D.M. Bers. 2000. Increased Na/Ca exchange and reduced IK1 facilitate triggered action potentials in a rabbit model of heart failure. Biophys. J. 78:55A.
Qi, X., P. Varma, D. Newman, N. Mamalias, and P. Dorian. 2002. Terikalant and barium decrease the area of vulnerability to ventricular fibrillation induction by T-wave shocks. J. Cardiovasc. Pharmacol. 39:242–250.
Reuveny, E., Y.N. Jan, and L.Y. Jan. 1996. Contributions of a negatively charged residue in the hydrophobic domain of the IRK1 inwardly rectifying K+ channel to K+-selective permeation. Biophys. J. 70:754–761.
Sharp, K., R. Fine, and B. Honig. 1987. Computer simulations of the diffusion of a substrate to an active site of an enzyme. Science. 236:1460–1463.
Sigworth, F.J. 1980. The variance of sodium current fluctuations at the node of Ranvier. J. Physiol. 307:97–129.
Sigworth, F.J., and B.C. Spalding. 1980. Chemical modification reduces the conductance of sodium channels in nerve. Nature. 283:293–295.
Stampe, P., J. Arreola, P. Perez-Cornejo, and T. Begenisich. 1998. Nonindependent K+ movement through the pore in IRK1 potassium channels. J. Gen. Physiol. 112:475–484.
Weiss, J.N., and S.T. Lamp. 1989. Cardiac ATP-sensitive K+ channels. Evidence for preferential regulation by glycolysis. J. Gen. Physiol. 94:911–935.
Worley, J.F., III, R.J. French, and B.K. Krueger. 1986. Trimethyloxonium modification of single batrachotoxin-activated sodium channels in planar bilayers. Changes in unit conductance and in block by saxitoxin and calcium. J. Gen. Physiol. 87:327–349.
Zobel, C., H.C. Cho, T.T. Nguyen, R. Pekhletski, R.J. Diaz, G.J. Wilson, and P.H. Backx. 2003. Molecular dissection of the inward rectifier potassium current (IK1) in rabbit cardiomyocytes: evidence for heteromeric co-assembly of Kir2.1 and Kir2.2. J. Physiol. 550:365–372.(Nazzareno D'Avanzo, Hee C)
2 Institute of Molecular Cardiobiology, Johns Hopkins University, Baltimore, MD 21218
3 Department of Chemistry and Department of Physics, University of Guelph, Guelph, Ontario N1G 2W1, Canada
ABSTRACT
Ion channel conductance can be influenced by electrostatic effects originating from fixed "surface" charges that are remote from the selectivity filter. To explore whether surface charges contribute to the conductance properties of Kir2.1 channels, unitary conductance was measured in cell-attached recordings of Chinese hamster ovary (CHO) cells transfected with Kir2.1 channels over a range of K+ activities (4.6–293.5 mM) using single-channel measurements as well as nonstationary fluctuation analysis for low K+ activities. K+ ion concentrations were shown to equilibrate across the cell membrane in our studies using the voltage-sensitive dye DiBAC4(5). The dependence of on the K+ activity (aK) was fit well by a modified Langmuir binding isotherm, with a nonzero intercept as aK approaches 0 mM, suggesting electrostatic surface charge effects. Following the addition of 100 mM N-methyl-D-glucamine (NMG+), a nonpermeant, nonblocking cation or following pretreatment with 50 mM trimethyloxonium (TMO), a carboxylic acid esterifying agent, the –aK relationship did not show nonzero intercepts, suggesting the presence of surface charges formed by glutamate or aspartate residues. Consistent with surface charges in Kir2.1 channels, the rates of current decay induced by Ba2+ block were slowed with the addition of NMG or TMO. Using a molecular model of Kir2.1 channels, three candidate negatively charged residues were identified near the extracellular mouth of the pore and mutated to cysteine (E125C, D152C, and E153C). E153C channels, but not E125C or D152C channels, showed hyperbolic –aK relationships going through the origin. Moreover, the addition of MTSES to restore the negative charges in E53C channels reestablished wild-type conductance properties. Our results demonstrate that E153 contributes to the conductance properties of Kir2.1 channels by acting as a surface charge.
Key Words: surface charges Kir2.1 conduction NMG TMO
N. D'Avanzo and H.C. Cho contributed equally to this project.
Abbreviations used in this paper: CHO, Chinese hamster ovarian; DiBAC4(5), bis-(1,3-dibutylbarbituric acid)pentamethine oxonol; GCS, Gouy-Chapman-Stern; NAChR, nicotinic acetylcholine receptor; NMG, N-methyl-D-glucamine; NSFA, nonstationary fluctuation analysis; TMO, trimethyloxonium; WT, wild-type.
INTRODUCTION
Electrostatic forces originating from charged groups on protein surfaces influence a variety of macromolecular functions. For example, charges on the surface of superoxide dismutases provides electrostatic guidance for substrates entering the active site (Getzoff et al., 1983, 1992; Sharp et al., 1987). Surface charges have also been shown to affect channel conductance in a variety of ion channels, such as neuronal Na+ channels (Sigworth and Spalding, 1980; Worley et al., 1986), Ca2+-activated K+ (BK) channels (MacKinnon and Miller, 1989; MacKinnon et al., 1989), and nicotinic acetylcholine receptors (NAChR) (Imoto et al., 1988; Konno et al., 1991), presumably by influencing the permeant ion concentration of ions at the entrance of the selectivity filter (Apell et al., 1977). In BK channels, surface charges also act to increase the affinity of cationic pore blockers by altering the local electrostatic potential at the pore entrance (MacKinnon and Miller, 1989; MacKinnon et al., 1989), while the "ring" of negative surface charges in NAChRs modifies the affinity of the pore blocker, Mg2+ (Imoto et al., 1988; Konno et al., 1991).
IK1 currents are found in many tissues and are generated by Kir2.x channels. In heart, a major molecular component of IK1 is the Kir2.1 channel (Zobel et al., 2003) and mutations in Kir2.1 channels have been linked to a form of long QT syndrome (LQT7) known as Andersen's syndrome (Plaster et al., 2001). Levels of IK1 in cardiomyocytes provide background K+ currents that are involved in setting the resting membrane potential, preventing membrane hyperpolarization due to Na+ pump activity, influencing propagation velocity, altering the electrical space constant, and promoting late phase repolarization. Consistent with an important role for IK1 in the heart, reductions in IK1 have been implicated in heart failure, and triggered arrhythmias resulting from early or delayed after-depolarizations observed in heart disease (Beuckelmann et al., 1993; Kb et al., 1996; Puglisi et al., 2000; Pogwizd et al., 2001). On the other hand, blockade of IK1 is a potent inhibitor of reentry type arrhythmias and ventricular fibrillation by increasing electrical space constants (Qi et al., 2002). IK1 channels also play a critical role in conditions of ischemia and hypoxia, where [K+]o is often elevated from 4 mM to as high as 20 mM (Weiss and Lamp, 1989), resulting in membrane potential depolarization as a consequence of the anomalous crossover properties of IK1 channels (Matsuda, 1991; Lopatin and Nichols, 1996).
The recent crystallography of the bacterial inward rectifier channel KcsA (Doyle et al., 1998) has provided important new insights into the molecular and biophysical mechanisms underlying the function of K+ channels. Indeed, crystallography and molecular dynamic models have revealed important channel properties such as ion occupancy and unitary conductance as well mechanisms for ion permeation (Berneche and Roux, 2001, 2003; Morais-Cabral et al., 2001). While previous reports of the versus K+ activity relationship in Kir2.1 channels were dedicated to studying the validity of various permeation models (Stampe et al., 1998), no studies have specifically examined the role of surface charges in Kir2.1 channel function.
In this study, we examined the conduction properties of the Kir2.1 channel over a range of K+ activities using single-channel recordings and nonstationary fluctuation analysis. We employed three experimental strategies to alter surface charges in Kir2.1 channels: variation of the bulk solution ionic strength, surface charge neutralization, and site-directed mutagenesis. Our results establish a critical role for surface charges near the extracellular entrance of the pore in Kir2.1 channels originating from E153 residues to the conduction properties of K+ ions. It will be necessary to consider the effect of E153 residues on ion occupancy and permeation in future studies.
MATERIALS AND METHODS
Expression of Mouse Kir2.1 Channels in Chinese Hamster Ovary (CHO) K1 Cells
Mouse Kir2.1 cDNA (a gift from L. Jan, Howard Hughes Medical Institute, San Francisco, CA) was subcloned into a bidirectional paBi-G vector (CLONTECH Laboratories, Inc.) modified by us to contain GFP (paBi-EGFP) (Zobel et al., 2003) in a TetOFF system requiring activation of the promoter by tetracycline transactivator (tTA) (CLONTECH Laboratories, Inc.). CHO-K1 cells were incubated at 37°C, 5% CO2 in F12 media (Sigma-Aldrich), supplemented with 5% FBS and 1% penicillin/streptomycin (GIBCO BRL) and transfected with 2 μg of Kir2.1 cDNA plus 1 μg of tTA cDNA using 4 μl of Lipofectamine 2000 (GIBCO BRL) following the supplied protocol. Transfected cells were incubated in supplemented culture media (see below) for 24 h before the electrophysiological recordings. Single mutants were created on a four cysteine knockout background (4CKO = C89S/C101S/C149Y/C169S) via the Kunkel method.
Cell-attached Macropatch Recordings for Unitary Conductance
Electrophysiological recordings were done using the cell-attached patch clamp technique. Bath solutions contained (in mM) 2 MgCl2, 10 HEPES, and 5–400 K+ (Cl– and OH–) (pH 7.4, with KOH). The osmolarity was maintained when necessary by the addition of varying amounts of glucose. Pipette solutions had the same composition as the bath solutions, except MgCl2 was removed and 1–30 μM Ba2+ was added in order to induce time-dependent current decay, allowing for nonstationary fluctuation analysis to be used for unitary current estimates (see below).
For some experiments, recording solutions also contained 100 mM N-methyl-D-glucamine chloride (NMG+Cl–) to "screen" charges. In some studies, transfected cells were pretreated in our culture medium (see below) for 5 min with 50 mM trimethyloxonium (TMO), an agent used to esterify carboxylic acid residues of negatively charged amino acid side chains (Sigworth and Spalding, 1980; Dudley and Baumgarten, 1993). To replace negative charges in some studies using mutant channels, 1 mM MTSES (in DMSO) was added for 10–15 min and then thoroughly washed out with fresh bath solution before recording. TW140F-4 patch pipettes (World Precision Instruments) were pulled and fire polished to give resistances of 6–10 M for single channel recordings and resistances 2 M for fluctuation analysis recordings. Pipettes were coated with Sylgard (Dow Corning Corporation) to reduce their capacitance. Data were acquired using Axopatch 200A amplifier, collected at 20 kHz by Clampex 8 (pCLAMP; Axon Instruments) with a 2 kHz low-pass Bessel filter and analyzed using Clampfit 8 (pCLAMP; Axon Instruments).
To obtain estimates for the unitary conductance () at varied levels of symmetric K+ in our cell attached recordings, it is desirable for the K+ concentrations to equilibrate across the cell membrane and for the resting membrane potential to be 0 mV. For thermodynamic reasons, this cannot be readily established directly from our cell-attached recordings. Therefore, to demonstrate K+ ion equilibration, we used the voltage-sensitive dye bis-(1,3-dibutylbarbituric acid)pentamethine oxonol (DiBAC4(5)) (Molecular Probes), which increases its fluorescence as the membrane potential depolarizes. To perform these studies, DiBAC4(5) was dissolved in 0.2% DMSO and added to bath solutions at a final concentration of 200 nM. Initially transfected cells were incubated for 20 min with a solution containing 140 mM KCl, 10 mM HEPES, 2 mM MgCl2, 200 nM DiBAC4(5), pH 7.4. Next, the fluorescence (at 616 nm) generated by DiBAC4(5) illuminated at 590 nm was recorded every 3 s from cells expressing Kir2.1 channels, identified by GFP fluorescence. Fig. 1 shows typical fluorescence recordings from DiBAC4(5) in response to changes in K+ levels from 140 mM to 20 mM K+ (20 mM K+, 10 mM HEPES, 220 mM glucose, 2 mM MgCl2, 200 nM DiBAC4(5), pH 7.4). As expected from the properties of DiBAC4(5), the application of 20 mM K+ results in rapid fluorescence decreases, indicative of membrane hyperpolarization, which was followed by a slow return of the fluorescence to baseline, recorded when 140 mM K+ was present. Furthermore, reintroduction of a 140 mM K+ containing solution, caused a rapid increase in fluorescence, indicating membrane depolarization followed again by decay back to baseline. Identical patterns of DiBAC4(5) fluorescence changes were observed when other solution changes were applied. In general, these observations can only be obtained if K+ ions equilibrate and have equal concentrations across the cell membrane. This conclusion is bolstered by the reasonable assumption that the initial membrane potential in our 140 mM K+ containing solution was 0 mV (since intracellular environment of our cells should be typical of other mammalian cells that contain 140 mM K+), thereby indicating that the baseline fluorescence in our recordings corresponds to a membrane potential of 0 mV.
The observation that the membrane potential equilibrates to 0 mV and the K+ levels are equal across the cell membrane (regardless of the applied external K+ level) is entirely expected for the following reasons. First, the conductance of these cells following transfection is very high (i.e., >40 nS) and K+ selective, which promotes K+ equilibration across the membrane. Second, since Na+ was excluded from the media, thereby inhibiting the Na+ pump, there is no biochemical mechanism (that we are aware of) to maintain K+ gradients across the membrane.
Nonstationary Fluctuation Analysis (NSFA) on Cell-attached Multichannel Recordings
Unitary currents were estimated by the nonstationary fluctuation analysis (NSFA) method using the equation (Sigworth, 1980)
(1)
where 2 is the variance of the measured macroscopic current, I, is the mean of I, N is the number of channels in the patch, and L is the contribution of leak currents. Since Kir2.1 channels do not vary greatly in response to voltage steps as is required for NSFA, Ba2+ was included in the extracellular solutions to induce time- and voltage-dependent block of Kir2.1 currents (Neyton and Miller, 1988; Jiang and MacKinnon, 2000). In the presence of Ba2+, 200–300 successive K+ current sweeps were recorded in response to voltage steps from a holding potential of +5 mV to –60, –80, –100, –120, and –140 mV. Variance and mean currents were estimated using consecutive pairs of traces (to minimize the complicating effects of current drift) and these estimates for 2 and were averaged. Estimates for 2 versus were fit with Eq. 1 using Maple 7.00 (Waterloo Maple Inc.) to estimate i. The single-channel conductance () was calculated from the dependence of unitary current estimates on voltage using Ohm's law (Kubo et al., 1993; Fig. 2 B).
Unitary Conductance–K+ Activity Relationship
To properly account for the effect of counter-ions in the solution on the driving force for K+ ion entry into the pore (i.e., the chemical potential), the activity of K+ (aK) was calculated using the relationship
(2)
To accurately estimate the contribution of long-range electrostatic interactions to the activity coefficients, fKact, in our studies we used the following equations. When the ionic strength (I) was <10 mM (i.e., I < 10 mM), fKact was calculated using the relationship (Debye and Hückel, 1923)
(3)
where z+ is the valence of the major cation, z– the major anion, and I is the ionic strength of the solution, which is defined as (Debye and Hückel, 1923)
(4)
For 10 mM < I < 100 mM, the Extended Debye-Hückel relationship was used (Kielland, 1937; MacInnes, 1961):
(5)
where B 0.33 in water at 25°C and a is an adjustable parameter corresponding to ionic size (=3 x 10–8 cm for K+). For I > 100 mM, we used Davies' approximation (MacInnes, 1961),
(6)
The dependence of on aK (i.e., the –aK relationship) was initially fit with a modified Langmuir binding isotherm (Fig. 3 B):
(7)
where aK is defined by Eq. 3, 5, or 6 as required, max is the maximal conductance measured at saturating aK, and EC50 is a measure of aK when is max /2. The parameter S is required to account for the nonzero intercept of the –aK relationship which is predicted to originate from the presence of surface charges (see below).
Surface Charge Calculations
To quantify the contribution of surface charges to Kir2.1 channel properties, we assumed that the conductance is a function of the activity of K+ ions at the channel surface (aKs) and not aK in the bulk solution. aKs was estimated by the relationship
(8)
where z, e, kB, and T have their usual meanings, and s is the electrostatic potential at the surface. In our studies, we estimated s using the relationship (McLaughlin et al., 1970)
(9)
where control is the channel conductance in the presence of surface charges, and TMO is the unitary conductance in following pretreatment with TMO that was assumed to neutralize all negative surface charges. Evidence for the validity of this assumption is provided in the results section.
To estimate the number of charged residues that might contribute to the observed surface charge effects, the surface charge density () was calculated using the following equation (Grahame, 1947; McLaughlin, 1977):
(10)
where aa = 0.15 M–1 (Eisenberg et al., 1979), aK is the activity in the bulk solution, r is assumed to be the dielectric of water (80), and o is the permittivity of free space. Values of for every s and aKs were averaged (±SEM) to determine the final estimate.
Kinetics of Barium Block
Since Ba2+ blocks Kir2.1 channels by binding within the selectivity filter, time constants () for Kir2.1 current decay in the presence of Ba2+ are predicted to be influenced by surface charges in a manner similar to K+ conductance (see above). By assuming that Ba2+ block is adequately described by a single binding site model (Neyton and Miller, 1988; Alagem et al., 2001) where only the on-rate is influenced by the surface potential combined with the assumptions that the rate of Ba2+ entry into the channel is given by the Ba2+ activity at the channel surface (aBas) and that Ba2+ does not alter the surface potential at the low concentrations used (see Eqs. 3 and 4), we can write
(11)
where kon(V) and koff(V) are the rates of Ba2+ binding and unbinding, respectively, at a given voltage, and aBas is the external Ba2+ activity at the surface of the channel, which is related to Ba2+ activity in the bulk (aBa) as described above for K+ ions in Eq. 6. Assuming that kon(V) depends exponentially on voltage (Reuveny et al., 1996), it can be shown that
(12)
where (1/) is the difference in 1/ for a given [Ba2+] in the presence and absence of NMG, and on is an estimate of the electrical distance for the position of the peak of the free energy barrier for Ba2+ binding to the channel pore. kon(0) and on were obtained from the intercept and slope of the relationship between 1/ and voltage (–80, –100, and –120 mV) in the presence of NMG (i.e., where the surface charges are screened) (Reuveny et al., 1996). Estimates of s using this method were made at each voltage and averaged in order to determine Ba.
Data Analysis
With the exception of NSFA, all other data were analyzed using Origin 6.0 (Microcal Software Inc.). This includes mathematical calculations, graph development, curve fitting with the appropriate functions, as well as statistical analysis.
RESULTS
Kir2.1 Single-channel Conductances at Varying K+ Concentrations
To investigate the existence of surface charges in Kir2.1 channels, we examined the conductance as a function of K+ activity. Fig. 2 A illustrates typical single channel currents (i) recorded in CHO cells transfected with Kir2.1 channels at several K+ concentrations. As shown in Fig. 2 B for 20 and 140 mM K+, i depends linearly on voltage (i.e., the channels are ohmic) over the range of voltages studied as previously reported (Kubo et al., 1993). Thus, unitary conductances () were estimated using Ohm's Law. Since the unitary current could not be reliably measured using single-channel recordings with <20 mM K+ for Kir2.1 channels, NSFA of currents recorded in cell-attached patches containing Ba2+ were used to estimate i. Typical cell-attached recordings are shown in Fig. 3 A, along with a typical plot of the mean variance (<2>) versus mean current () assessed from analysis of 200 sweeps, which allowed estimates of i from the parabolic function in Eq. 1. Fig. 3 B plots the estimated unitary conductance from single channel recordings (closed squares) and NSFA (open squares) as a function of the bulk K+ activity. Note that the estimates using NSFA were identical (P < 0.01) to the values estimated using single-channel recordings (at 20 and 140 mM K+). These data were fit (R2 = 0.98) using the modified Langmuir binding isotherm (see MATERIALS AND METHODS) to estimate max (26.0 ± 1.5 pS, 10 < n < 20) and EC50 (50 ± 15 mM aK) as well as the intercept, S (4.4 ± 1.4 pS). These results establish a clear saturation of unitary conductance at high K+ concentrations as reported previously (Stampe et al., 1998) and suggest the existence of a surface charge.
To ensure that our conductance estimates did not arise from experimental limitations of the NSFA method, we performed Monte Carlo simulations of channel blocking events (Hanna et al., 1996) using experimentally derived mean blocking and unblocking times in the presence of 30 μM Ba2+ (Table I) assuming a single site blocking model, as established previously (Neyton and Miller, 1988; Alagem et al., 2001). Typical simulated single-channel events are shown in Fig. 4 A (bottom two traces), while Fig. 4 B shows the predicted macroscopic current generated by superimposing 100 sweeps, as well as the predicted dependence of <2> on . Clearly, the simulated results in Fig. 4 B (left) closely resemble the experimental results in Fig. 3 A. More important, the estimated unitary current estimated from the <2> versus plot is identical to the assumed unitary current (Makhina et al., 1994, Krapivinsky et al., 1998).
Charge Screening Decreases Kir2.1 Single Channel Conductance at a Low K+
The presence of the nonzero intercept, S, in the versus aK relationship suggests the presence of a surface charge. To investigate this possibility, the versus aK relationship was measured in the presence of 100 mM NMG to screen surface charges. As shown in Fig. 5 B, with 200 mM symmetrical K+ (see MATERIALS AND METHODS), the inclusion of 100 mM NMG did not affect the Kir2.1 single channel conductance (solid line is without NMG), confirming that NMG does not block Kir2.1 channels. By contrast, at 5 mM K+, the single channel conductance was reduced (P < 0.01) to 2.0 ± 0.3 pS in the presence of 100 mM NMG compared with 6.6 ± 1.5 pS in the absence of NMG (Fig. 5 A). Reductions in were also observed at 10, 15, 20, and 40 mM K+, when 100 mM NMG was added (P < 0.05). More important, fitting the –K+ activity relationship with Eq. 7, when 100 mM NMG was present, reduced the conductance intercept to a value not different from zero (i.e., 0.8 ± 1.3, P = 0.36) (Fig. 5 C).
As an alternate method for eliminating surface charges, CHO cells were treated with TMO, a carboxylate esterifying agent. The unitary conductance with TMO was 2.4 ± 1.0 pS with 5 mM K+ which was below (P < 0.01) the value of estimated without TMO treatment (Fig. 5 A). Furthermore, the unitary conductance in 5 mM K+ after treatment with TMO was not reduced further by the inclusion of 100 mM NMG (2.9 ± 1.0 pS) (Fig. 5 A), and the intercept for the versus aK relationship (0.8 ± 1.9 pS) was not different (P = 0.32) from zero. Taken together, these results suggest that negatively charged amino acid residues generate surface charges in Kir2.1 channels, thereby increasing unitary currents at low permeant ion concentrations.
Electrostatic Theory Predictions of Surface Charge Effects
Using our results above, it is possible to calculate the surface charge potential (s) for Kir2.1 channels (Eq. 9) and thereby the activity of K+ ions at the mouth of the channel pore (i.e., aKs) (using Eqs. 8 and 9). Fig. 6 shows the relationship between and aKS (closed squares), which can be well fit (R2 = 0.99) with Eq. 7 going through the origin (S = 0.7 ± 1.5 pS, P > 0.05). Notice that the relationship of versus aKs in control conditions is identical to the versus aK (i.e., bulk K+ activity) relationship for channels in the presence of NMG (open triangles), consistent with a role for surface charges in concentrating permeant K+ ions at the channel surface.
Our estimates of s were also used to calculate the apparent surface charge density () of Kir2.1 channels (Eq. 10). Values of charge density in Kir2.1 channels were calculated for each aKs examined assuming a dielectric (r) of 80 (that of water), and then averaged to determine a of –0.49 ± 0.08 e nm–2. When NMG was present, the estimated surface charge was reduced to –0.06 ± 0.02 e nm–2, which represents less than one net charge per channel subunit. This is consistent with the screening of surface charges in Kir2.1 channels by NMG. Using current molecular models of Kir2.1 channels (Cho et al., 2000), we estimated the channel's surface area to be 11 nm2, which predicted the presence of approximately –5.4 e net charges per channel. Since Kir2.1 channels are tetramers, these calculations suggest that the surface charge arises from one free net charge per monomer.
Kinetics of Extracellular Barium Block Are Altered by Charge Screening
Ba2+ is a potent voltage- and time-dependent blocker of K+ current in Kir2.1, by binding within the selectivity filter (Jiang and MacKinnon, 2000). Thus, the presence of surface charges is predicted to alter the blocking time constant () of Ba2+. Indeed, as shown in Fig. 7 A, treatment with NMG or TMO significantly increased at 5 mM K+ compared with control. The same was observed for 20 mM (Fig. 7, compare B and C) and 140 mM [K+]. Using the voltage dependence of kon in the absence of surface charges (Fig. 7 D) and Eqs. 11 and 12, s was calculated using Ba2+ blocking kinetics for 20 mM K+, which again allows an estimate of the charge density (Ba) in Kir2.1. Ba was estimated to be –0.49 ± 0.12 e nm–2, which is indistinguishable (P < 0.01) from the value obtained using the conductance data, confirming the presence of four negative charges on the surface of Kir2.1 channels.
Identification of Extracellular Charges Acting as Surface Charges
Using the molecular model for Kir2.1 channels, which is based on the KcsA crystal structure (Cho et al., 2000), we identified three candidate charged residues (E125, D152, and E153) that might act as surface charges (Fig. 8 A). We replaced these charged residues with cysteine in order to examine the potential role of these charges in the electrostatic effects of Kir2.1 channels. Cysteine was chosen to allow charge replacement using MTS compounds. The dependence of on aK was unaltered for E125C and D152C mutant channels, compared with wild-type (WT) channels at 20 mM K+ (Fig. 8 B). On the other hand, E153C channels showed a smaller unitary conductance compared with WT, which was not different from the unitary conductance observed in WT channels pretreated with TMO (Fig. 8 B) (P < 0.01). Fig. 8 C shows that the same results were observed for 5 mM K+, whereas, at 140 mM K+, WT channels, with and without TMO treatment, and all three mutant channels had similar values. Treatment of E153C channels with MTSES, in order to replace the negative charge, increased the to levels similar to WT at 20 mM K+ (Fig. 8 B). Taken together, these results support the conclusion that E153 acts as the primary source of surface charges in Kir2.1 channels.
DISCUSSION
Our results establish that Kir2.1 channels display classical surface charge properties (Sigworth and Spalding, 1980; Worley et al., 1986; MacKinnon and Miller, 1989; MacKinnon et al., 1989) that depend on E153 residues predicted from a previous molecular model to be located on the outer face of the channel (Cho et al., 2000). Investigations into the possible role of surface charges in our studies necessitated the conversion of K+ ion concentrations into chemical activity (i.e., chemical potential) in order to account for the effects of ionic strength on the thermodynamic driving force of K+ ions into Kir2.1 channels. K+ activity was calculated using a combination of the linear Debye-Huckel, the Extended Debye-Huckel, and the Davies' theories to account for long-range electrostatic interactions between ions in solution under our experimental conditions (Debye and Hückel, 1923; MacInnes, 1961). The use of multiple equations was necessary because no single analytical expression can comprehensively estimate ionic activities over the range of experimental conditions used in our studies (MacInnes, 1961). The use of K+ activity was essential in our studies since a wide range of ionic strengths were used to alter surface potentials in Kir2.1 channels.
In our studies, it was necessary to determine the unitary current and unitary conductance () with varied symmetric levels of K+. This could not be achieved using whole-cell recordings since the current magnitudes were too large to allow NSFA to be effectively used (see below). Therefore, to circumvent this problem, we resorted to cell-attached recordings. To ensure the levels of K+ were symmetric in our cell-attached recordings, we monitored the membrane potential using a voltage-sensitive dye. These studies showed that K+ levels equilibrate within 2–3 min following changes of the external K+ levels under our experimental conditions, probably due to the presence of very high levels of Kir2.1 channel expression (i.e., >40 nS in whole-cell recordings with 140 mM K+) combined with the omission of Na+, thereby inhibiting the Na+ pump.
Clear evidence for the presence of surface charges was obtained in our studies using two methods that yielded comparable estimates for the electrostatic surface potential in Kir2.1 channels. In the first method, the dependence of unitary current () on the K+ activity was used. These studies required the use of NSFA at low K+ activities, due to small unitary currents, as done previously in HRK1 (Makhina et al., 1994) and Kir7.1 (Krapivinsky et al., 1998) channels. We validated the NSFA method in our studies using Monte Carlo simulations to model Ba2+ blocking events. In the versus aK plots, we found that does not approach zero as K+ activities are reduced. While it is possible that the finite limiting intercept originates from finite permeability of Kir2.1 channels to other cations in our solutions (e.g., Ba2+), Kir2.1 channels are impermeant to these other ions (Alagem et al., 2001). More important, the limiting was reduced to zero at low aK following the addition of the impermanent cations NMG+ to screen surface charges or following pretreatment of the channels with TMO to specifically neutralize carboxylate residues (Parsons et al., 1969; Nakayama et al., 1970). The nonzero conductance intercept at low K+ activity, combined with a zero conductance intercept at high ionic strength or following TMO pretreatment are signature features of local surface charges (Manning, 1978; Fixman, 1979; Klein et al., 1981; Cai and Jordan, 1990).
To analyze and interpret our versus aK results we used the Gouy-Chapman-Stern (GCS) theory for surface charges, which assumes the surface charge density is uniformly distributed in a flat plane. While this assumption is generally expected to be invalid for ion channel proteins, the GCS theory accurately accounts for surface charge effects in many proteins, possibly due to cancellation of errors in the functional expansions of rigorously correct canonical partition functions (Fixman, 1979; Klein et al., 1981). We anticipate that the applicability of the GCS theory is assisted by the symmetrical arrangement of four negative point charges around the pore opening of tetrameric Kir2.1 channels that should mimic more closely a uniform charge distribution than a single point charge located near the pore opening. Regardless, the application of the GCS theory to our results predicted that depends uniquely on the local K+ ion activity at the entrance to the Kir2.1 channel pore. Indeed, the dependence of on aK at the channel surface (i.e., aKs) could be accurately described by a Langmuir binding isotherm going through the origin. The GCS theory further predicted a surface charge density of –0.49 ± 0.08 e nm–2. Using recent molecular models of Kir2.1 channels (Cho et al., 2000) to estimate channel area, we calculated a net of approximately four charges on the surface of Kir2.1 channels in the absence of charge screening agents, which equated to nearly one charge per monomer.
Thus, our results strongly support the conclusion that surface charges influence the Kir2.1 channel conductance. Similar investigations and analyses have been applied to other channels, including Na+ channels (Sigworth and Spalding, 1980; Worley et al., 1986), Ca2+-activated K+ channel (Moczydlowski et al., 1985; MacKinnon and Miller, 1989; MacKinnon et al., 1989; Park et al., 2003), NAChRs (Imoto et al., 1988; Konno et al., 1991), to establish the presence of surface charges. On the other hand, versus concentration in KcsA channels show two asymptotic conductance values at low and high K+ levels (Morais-Cabral et al., 2001; Nimigean et al., 2003). The limiting value of KcsA at low K+ levels appears remarkably similar to that observed in Kir2.1 channels, suggesting that surface charges may also play an important role in the conductance properties of these channels as well.
Since Ba2+ blocks Kir2.1 channels by binding within the channel pore, Ba2+ should be affected by the presence of surface charges in a comparable manner to K+ ions. Specifically, the rate of Ba2+ entry into the channel should be proportional to the activity of Ba2+ at the channel entry. This provided a second method for estimating surface charges in Kir2.1. Indeed, when charges are screened by NMG, or neutralized by TMO, the kinetics of Ba2+ block is slowed significantly. Furthermore, after accounting for ionic strength (i.e., activity) of Ba2+ in our recording conditions, our estimates of the surface charge density from our Ba2+ block studies (Ba), using the GCS theory, were identical to those estimated from our conductance studies, once again supporting the existence of four negative charges on the surface of Kir2.1 channel proteins. This contrasts with a previous report examining the effects of lipid charges on BK channels (Park et al., 2003), where alterations in lipid charge (i.e., lipid surface charges) affected BK channel conductance, but not Ba2+ block properties. This led the authors to conclude, consistent with molecular models of channel proteins, that lipid charges are too remote from the channel pore entrance to produce classical surface charge effects (Park et al., 2003). While we did not examine the effects of altered lipid composition on the Kir2.1 channel properties, our results are consistent with the conclusions of Park et al. (2003) since surface charge effects in Kir2.1 channels originate from charged E153 (see below), predicted to be more adjacent to the channel pore entrance than charged lipids. Further studies will be required to examine whether lipid charge also affects Kir2.1 channel properties.
More definitive evidence for the existence of surface charges in Kir2.1 channels was provided by our mutagenesis studies. Somewhat surprisingly, E125C and D152C mutant channels showed conductance properties identical to WT channels. Previous studies have suggested that E125 contributes a superficial site for binding of extracellular divalent cations such as Mg2+, Ba2+, and Ca2+ (Murata et al., 2002). Indeed, in the presence of divalents, E125Q mutant channels had decreased single-channel conductance compared with WT channels. Furthermore, the E125 was also postulated to play a key role in dehydration of Ba2+ before entering the pore (Alagem et al., 2001). Consistent with these results, Ba2+ blocking kinetics were increased in E125C mutant channels compared with WT channels. The significance of these observations is unclear but could arise from altered binding affinities for these divalents at this site.
Our observations in D152C channels were unexpected, based on our model of Kir2.1 channels (Cho et al., 2000). This may indicate that the D152 residue binds K+ ions with high affinity, thereby not contributing to the surface charge effect. Alternatively, our molecular model (Cho et al., 2000) may require revision with the D152 residue either buried or oriented in such a way that its side chain does not contribute to the surface potential. Although less likely, it is conceivable that the sulfhydryl side chains of the C125 and C152 residues may be ionized in the E125C and D152C mutant channels. Consequently, the residues would not have led to charge neutralization, thereby having maintained their surface potential. However, this possibility would be inconsistent with the expectation that polar residues are less likely (not more likely) to be ionized (i.e., have a higher pKa) when localized in regions with reduced dielectric constants, such as in proteins compared with the aqueous phase.
Unlike E125C and D152C mutant channels, E153C channels have reduced single-channel conductance compared with WT channels at low ionic activities and are identical in conductance properties to WT channels pretreated with 50 mM TMO (Fig. 8 B). Moreover, addition of 1 mM MTSES to restore negative charges restored WT behavior to E153C channels (Fig. 8 B). These results support the conclusion that E153 is responsible for the surface charge effects observed in Kir2.1 channels by contributing approximately one negative charge per monomer at the channel surface. In accordance with electrostatic theory, our results support the conclusion that ionized carboxyl groups increase the conduction at low permeant ion activities, by attracting the permeant ions near the conduction pathway.
Our results also establish that for Kir2.1 channels, increases with elevations in aK, exhibiting a clear maximum conductance >160 mM (K+ concentrations >300 mM) with an apparent aK EC50 of 50 ± 15 mM. These properties of Kir2.1 channels resemble those of ROMK1 (Kir1.1) channels that have a similar EC50, although ROMK1 channels exhibit a maximum conductance as aK is elevated above 300 mM K+ (Heginbotham and MacKinnon, 1993; Lu and MacKinnon, 1994). It is conceivable that of Kir2.1 channels might also reach a maximum at higher K+ concentrations. In stark contrast to Kir1.1 and Kir2.1 channels, the –[K+] relationship in KcsA channels is nearly linear from 50 to 500 mM K+ (Morais-Cabral et al., 2001), with evidence of saturation at K+ concentrations >1600 mM K+ (Nimigean et al., 2003). The basis for the differences in the properties of various inward rectifier channels will clearly require further investigations.
While our studies focused on the contribution of surface charges from E153 residues to Kir2.1 channel conductance, it is conceivable that surface charges may affect other channel properties. Indeed, since polyamine and Mg2+ block Kir2.1 channels by binding within the channel pore, alterations in surface charges could conceivably alter binding stability of these blockers. Thus, surface charge effects may contribute to the crossover effect observed in Kir2.1 channels with increased extracellular [K+].
Our findings support the presence of approximately four negative charges on the surface of Kir2.1 channels, originating from E153 (but not E125 or D152) residues. The results of our studies should help guide the development of new molecular models of Kir2.1 channels in the absence of crystal structures. It is apparent that new models of Kir2.1 channels must reexamine the location, orientation, and ionization of these three residues to account for the surface charge effect we attribute to the E153 residue. As well, future studies of ion occupancy, permeation, and perhaps rectification in Kir2.1 channels need to consider surface charge effects as described in our studies.
ACKNOWLEDGMENTS
We thank L.Y. Jan for providing us with the Kir2.1 clone and A. Prendergast for his technical assistance.
This study was supported by Heart and Stroke Foundation of Ontario (grant T-4179) to P.H. Backx. P.H. Backx is a Career Investigator of the Heart and Stroke Foundation of Ontario. H.C. Cho was supported by a Canadian Institutes of Health Research studentship. N. D'Avanzo holds a studentship from Heart & Stroke Richard Lewar Center of Excellence at the University of Toronto.
David C. Gadsby served as editor.
Submitted: 25 August 2004
Accepted: 18 March 2005
REFERENCES
Alagem, N., M. Dvir, and E. Reuveny. 2001. Mechanism of Ba2+ block of a mouse inwardly rectifying K+ channel: differential contribution by two discrete residues. J. Physiol. 534:381–393.
Apell, H.J., E. Bamberg, H. Alpes, and P. Lauger. 1977. Formation of ion channels by a negatively charged analog of gramicidin A. J. Membr. Biol. 31:171–188.
Berneche, S., and B. Roux. 2001. Energetics of ion conduction through the K+ channel. Nature. 414:73–77.
Berneche, S., and B. Roux. 2003. A microscopic view of ion conduction through the K+ channel. Proc. Natl. Acad. Sci. USA. 100:8644–8648.
Beuckelmann, D.J., M. Nabauer, and E. Erdmann. 1993. Alterations of K+ currents in isolated human ventricular myocytes from patients with terminal heart failure. Circ. Res. 73:379–385.
Cai, M., and P.C. Jordan. 1990. How does vestibule surface charge affect ion conduction and toxin binding in a sodium channel. Biophys. J. 57:883–891.
Cho, H.C., R.G. Tsushima, T. Nguyen, H.R. Guy, and P.H. Backx. 2000. Two critical cysteine residues implicated in disulfide bond formation and proper folding of Kir2.1. Biochemistry. 39:4649–4657.
Debye, P., and E. Hückel. 1923. Zur Theorie der Elektrolyte. II. Das Grenzgesetz fur die elektrische Leitfahigkeit. Physiol. Z. 24:305–325.
Doyle, D.A., J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, and R. MacKinnon. 1998. The structure of the potassium channel: molecular basis of K+ conduction and selectivity. Science. 280:69–77.
Dudley, S.C., Jr., and C.M. Baumgarten. 1993. Modification of cardiac sodium channels by carboxyl reagents. Trimethyloxonium and water-soluble carbodiimide. J. Gen. Physiol. 101:651–671.
Eisenberg, M., T. Gresalfi, T. Riccio, and S. McLaughlin. 1979. Adsorption of monovalent cations to bilayer membranes containing negative phospholipids. Biochemistry. 18:5213–5223.
Fixman, M. 1979. The Poisson-Boltzmann equation and its application to polyelectrolytes. J. Chem. Physiol. 70:4995–5005.
Getzoff, E.D., D.E. Cabelli, C.L. Fisher, H.E. Parge, M.S. Viezzoli, L. Banci, and R.A. Hallewell. 1992. Faster superoxide dismutase mutants designed by enhancing electrostatic guidance. Nature. 358:347–351.
Getzoff, E.D., J.A. Tainer, P.K. Weiner, P.A. Kollman, J.S. Richardson, and D.C. Richardson. 1983. Electrostatic recognition between superoxide and copper, zinc superoxide dismutase. Nature. 306:287–290.
Grahame, D.C. 1947. The electrical double layer and the theory of electrocapiliarity. Chem. Rev. 41:441–501.
Hanna, W.J., R.G. Tsushima, R. Sah, L.J. McCutcheon, E. Marban, and P.H. Backx. 1996. The equine periodic paralysis Na+ channel mutation alters molecular transitions between the open and inactivated states. J. Physiol. 497(Pt 2):349–364.
Heginbotham, L., and R. MacKinnon. 1993. Conduction properties of the cloned Shaker K+ channel. Biophys. J. 65:2089–2096.
Imoto, K., C. Busch, B. Sakmann, M. Mishina, T. Konno, J. Nakai, H. Bujo, Y. Mori, K. Fukuda, and S. Numa. 1988. Rings of negatively charged amino acids determine the acetylcholine receptor channel conductance. Nature. 335:645–648.
Jiang, Y., and R. MacKinnon. 2000. The barium site in a potassium channel by x-ray crystallography. J. Gen. Physiol. 115:269–272.
Kb, S., H.B. Nuss, N. Chiamvimonvat, B. O'Rourke, P.H. Pak, D.A. Kass, E. Marban, and G.F. Tomaselli. 1996. Ionic mechanism of action potential prolongation in ventricular myocytes from dogs with pacing-induced heart failure. Circ. Res. 78:262–273.
Kielland, J. 1937. Individual activity coefficients of ions in aqueous solutions. J. Am. Chem. Soc. 59:1675–1678.
Klein, B.K., C.F. Anderson, and M.T. Record. 1981. Condensation model expressions for the colligative properties of cylindrical polyions. Biopolymers. 20:2263–2280.
Konno, T., C. Busch, E. Von Kitzing, K. Imoto, F. Wang, J. Nakai, M. Mishina, S. Numa, and B. Sakmann. 1991. Rings of anionic amino acids as structural determinants of ion selectivity in the acetylcholine receptor channel. Proc. R. Soc. Lond. B. Biol. Sci. 244:69–79.
Krapivinsky, G., I. Medina, L. Eng, L. Krapivinsky, Y. Yang, and D.E. Clapham. 1998. A novel inward rectifier K+ channel with unique pore properties. Neuron. 20:995–1005.
Kubo, Y., T.J. Baldwin, Y.N. Jan, and L.Y. Jan. 1993. Prmary structure and functional expression of a mouse inward rectifier potassium channel. Nature. 362:127–133.
Lopatin, A.N., and C.G. Nichols. 1996. [K+] dependence of polyamine-induced rectification in inward rectifier potassium channels (IRK1, Kir2.1). J. Gen. Physiol. 108:105–113.
Lu, Z., and R. MacKinnon. 1994. A conductance maximum observed in an inward-rectifier potassium channel. J. Gen. Physiol. 104:477–486.
MacInnes, D.A. 1961. The Principles of Electrochemistry, New edition. Dover, New York. 478 pp.
MacKinnon, R., and C. Miller. 1989. Functional modification of a Ca2+-activated K+ channel by trimethyloxonium. Biochemistry. 28:8087–8092.
MacKinnon, R., R. Latorre, and C. Miller. 1989. Role of surface electrostatics in the operation of a high-conductance Ca2+-activated K+ channel. Biochemistry. 28:8092–8099.
Makhina, E.N., A.J. Kelly, A.N. Lopatin, R.W. Mercer, and C.G. Nichols. 1994. Cloning and expression of a novel human brain inward rectifier potassium channel. J. Biol. Chem. 269:20468–20474.
Manning, G.S. 1978. The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides. Q. Rev. Biophys. 11:179–246.
Matsuda, H. 1991. Effects of external and internal K+ ions on magnesium block of inwardly rectifying K+ channels in guinea-pig heart cells. J. Physiol. 435:83–99.
McLaughlin, S. 1977. Electrostatic potentials at membrane-solution interfaces. Curr. Top. Membr. Transp. 9:71–144.
McLaughlin, S.G., G. Szabo, G. Eisenman, and S.M. Ciani. 1970. Surface charge and the conductance of phospholipid membranes. Proc. Natl. Acad. Sci. USA. 67:1268–1275.
Moczydlowski, E., O. Alvarez, C. Vergara, and R. Latorre. 1985. Effect of phospholipid surface charge on the conductance and gating of a Ca2+-activated K+ channel in planar lipid bilayers. J. Membr. Biol. 83:273–282.
Morais-Cabral, J.H., Y. Zhou, and R. MacKinnon. 2001. Energetic optimization of ion conduction rate by the K+ selectivity filter. Nature. 414:37–42.
Murata, Y., Y. Fujiwara, and Y. Kubo. 2002. Identification of a site involved in the block by extracellular Mg2+ and Ba2+ as well as permeation of K+ in the Kir2.1 K+ channel. J. Physiol. 544:665–677.
Nakayama, H., K. Tanizawa, and Y. Kanaoka. 1970. Modification of carboxyl groups in the binding sites of trypsin with the Meerwein reagent. Biochem. Biophys. Res. Commun. 40:537–541.
Neyton, J., and C. Miller. 1988. Potassium blocks barium permeation through a calcium-activated potassium channel. J. Gen. Physiol. 92:549–567.
Nimigean, C.M., J.S. Chappie, and C. Miller. 2003. Electrostatic tuning of ion conductance in potassium channels. Biochemistry. 42:9263–9268.
Park, J.B., H.J. Kim, P.D. Ryu, and E. Moczydlowski. 2003. Effect of phosphatidylserine on unitary conductance and Ba2+ block of the BK Ca2+-activated K+ channel: re-examination of the surface charge hypothesis. J. Gen. Physiol. 121:375–397.
Parsons, S.M., L. Jao, F.W. Dahlquist, C.L. Borders Jr., J. Racs, T. Groff, and M.A. Raftery. 1969. The nature of amino acid side chains which are critical for the activity of lysozyme. Biochemistry. 8:700–712.
Plaster, N.M., R. Tawil, M. Tristani-Firouzi, S. Canun, S. Bendahhou, A. Tsunoda, M.R. Donaldson, S.T. Iannaccone, E. Brunt, R. Barohn, et al. 2001. Mutations in Kir2.1 cause the developmental and episodic electrical phenotypes of Andersen's syndrome. Cell. 105:511–519.
Pogwizd, S.M., K. Schlotthauer, L. Li, W. Yuan, and D.M. Bers. 2001. Arrhythmogenesis and contractile dysfunction in heart failure: roles of sodium-calcium exchange, inward rectifier potassium current, and residual -adrenergic responsiveness. Circ. Res. 88:1159–1167.
Puglisi, J.L., S.M. Pogwizd, W. Yuan, and D.M. Bers. 2000. Increased Na/Ca exchange and reduced IK1 facilitate triggered action potentials in a rabbit model of heart failure. Biophys. J. 78:55A.
Qi, X., P. Varma, D. Newman, N. Mamalias, and P. Dorian. 2002. Terikalant and barium decrease the area of vulnerability to ventricular fibrillation induction by T-wave shocks. J. Cardiovasc. Pharmacol. 39:242–250.
Reuveny, E., Y.N. Jan, and L.Y. Jan. 1996. Contributions of a negatively charged residue in the hydrophobic domain of the IRK1 inwardly rectifying K+ channel to K+-selective permeation. Biophys. J. 70:754–761.
Sharp, K., R. Fine, and B. Honig. 1987. Computer simulations of the diffusion of a substrate to an active site of an enzyme. Science. 236:1460–1463.
Sigworth, F.J. 1980. The variance of sodium current fluctuations at the node of Ranvier. J. Physiol. 307:97–129.
Sigworth, F.J., and B.C. Spalding. 1980. Chemical modification reduces the conductance of sodium channels in nerve. Nature. 283:293–295.
Stampe, P., J. Arreola, P. Perez-Cornejo, and T. Begenisich. 1998. Nonindependent K+ movement through the pore in IRK1 potassium channels. J. Gen. Physiol. 112:475–484.
Weiss, J.N., and S.T. Lamp. 1989. Cardiac ATP-sensitive K+ channels. Evidence for preferential regulation by glycolysis. J. Gen. Physiol. 94:911–935.
Worley, J.F., III, R.J. French, and B.K. Krueger. 1986. Trimethyloxonium modification of single batrachotoxin-activated sodium channels in planar bilayers. Changes in unit conductance and in block by saxitoxin and calcium. J. Gen. Physiol. 87:327–349.
Zobel, C., H.C. Cho, T.T. Nguyen, R. Pekhletski, R.J. Diaz, G.J. Wilson, and P.H. Backx. 2003. Molecular dissection of the inward rectifier potassium current (IK1) in rabbit cardiomyocytes: evidence for heteromeric co-assembly of Kir2.1 and Kir2.2. J. Physiol. 550:365–372.(Nazzareno D'Avanzo, Hee C)