|
||
Original Article |
Physical Basis for Lyotropic Anion Selectivity Patterns
b From the Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109
Department of Physiology & Pharmacology, Oregon Health Sciences University, 3181 SW Sam Jackson Park Rd., L-334, Portland, OR 97201.503-494-4352
dawsonda{at}ohsu.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: hydration energy anion binding pseudohalides ion-selective electrodes anion channels
| INTRODUCTION |
|---|
|
|
|---|
0.03) (Tabcharani et al. 1997
To gain insight into the nature of the physical interactions that are reflected in the peak and well energies that characterize anion permeation, we compared the selectivity of CFTR to that of a well-characterized, synthetic anion-selective membrane composed of plasticized poly(vinyl chloride) (PVC) doped with tridodecylmethylammonium chloride (TDMAC) for which the physical basis for the response to anion substitution is more readily apparent (Morf 1981
). The selectivity patterns for CFTR and the synthetic membrane differed by a multiplicative constant and both could be predicted by a continuum electrostatic model, based on the Born charging energy, that modeled anion–channel interactions as the energy of a charged sphere in a polarizable medium.
The continuum electrostatic approach provides a unified, quantitative interpretation of the observed energetics of permeation and block and offers a plausible explanation for the differential effects of mutations on these two processes that may be useful in understanding the physical nature of the conduction path, and in evaluating proposed structural models for the pore domain.
| MATERIALS AND METHODS |
|---|
|
|
|---|
r2). PC Spartan lacks basis values for gold so we could not perform the full set of calculations for Au(CN)2; however, we were able to model Ag(CN)2, which should have approximately the same dimensions. One of the advantages of the molecular orbital calculation method is that it permits an assessment of the charge distribution within the molecule. A hallmark of the halides and pseudohalides is that the negative charge is uniformly distributed (Douglas et al. 1994
|
(Bockris and Reddy 1970
= 80). There is reasonable agreement between the hydration energy calculated this way and the measured values for the halides and pseudohalides (Marcus 1997
Ghyd| plotted versus reciprocal ionic radius. It is apparent that the polyatomic anions (often used for sizing anion-selective pores; Bormann et al. 1987
|
The electrode polymer membrane was composed of 1 wt% TDMAC, 33 wt% PVC, and 66 wt% ortho-nitrophenyloctyether (o-NPOE), a plasticizer. The ions were tested at a 10-mM concentration (as a Na or K salt) in distilled/deionized water, as well as in a 10-mM HEPES solution buffered at pH 7.4. The values were calculated from the average of six electrodes. The selectivity coefficient was calculated using the separate solution method (Morf 1981
; Ammann 1986
) as follows (
):
![]() | (1) |
The measured potentiometric selectivity coefficients (logK pot) for the ions used in this study are presented in Table . For ions like dicyanamide and tricyanomethanide, for which values for the free energy of hydration (
Ghyd) were not available, we used the measured logK pot to estimate
Ghyd values by fitting a line through the data for ClO4, SCN, I, NO3, Br, and Cl using reported values of
Ghyd (Marcus 1997
) and measured Kpot values. The extrapolated values are given in Table .
RNA Synthesis and Xenopus Oocyte Expression
Human wild-type CFTR (wtCFTR) in a pBluescript vector (Stratagene Inc.) (Drumm et al. 1991
) was linearized with XhoI and used as template for the generation of cRNA using the mMessage Machine protocol (Ambion, Inc.). The cRNA was resuspended in DEPC-treated water and maintained at –70°C before injection into oocytes. Female Xenopus (Xenopus-1) were anesthetized by immersion in ice water containing 3-aminobenzoic acid ethyl ester (Tricaine, 3 mg/ml; Sigma Chemical Co.) and oocytes were removed through an abdominal incision. The follicular membrane was removed by mechanical agitation (2–3 h) in a nominally Ca2+-free collagenase solution containing (mM): 82.5 NaCl, 2 KCl, 1 MgCl, 10 HEPES, pH 7.5, and 2.5 mg/ml collagenase (GIBCO BRL). Defolliculated oocytes were maintained in a modified Barth's solution (MBSH) containing (mM): 88 NaCl, 1 KCl, 2.4 NaHCO3, 0.82 MgSO4, 0.33 Ca(NO3)2, 0.41 CaCl2, 10 HEPES, pH 7.5, and 150 mg/liter gentamicin sulfate. Oocytes were maintained at 18°C in a humidified incubator. 1 d after isolation, oocytes were injected with cRNA (diluted to give 50–250 µS of stimulated conductance:
0.15 ng/oocyte in a 50-nl volume) using a microinjector (Drummond Scientific Co.) and beveled injection needles (
10-µm tip diameter). Injected oocytes were maintained in MBSH and used for electrophysiological analysis 2–6 d after injection.
Electrophysiology
Individual oocytes were perfused with an amphibian Ringer's solution containing (mM): 100.5 NaCl, 2 KCl, 1.8 CaCl, 1 MgCl, and 5 HEPES, pH 7.5. The oocyte was impaled with two microelectrodes with tips pulled (P-97; Sutter Instruments Co.) to give 0.5–1.5 M
of resistance when filled with 3 M KCl. The open circuit membrane potential was continuously monitored on a strip chart recorder (Kipp & Zonen), and periodically the membrane was clamped (TEV-200; Dagan Corp.) and using a computer-driven protocol (Clampex; Axon Instruments), ramped from –120 to +60 mV at a rate of 100 mV/s for most analyses, although a step protocol (from –120 to +40 mV in 10-mV steps, 200 ms/step) was also used to check for time-dependent currents. The membrane conductance was calculated using the slope conductance over a 20-mV range centered on the reversal potential, and using chord conductances at various voltages. For ramp data, a correction for the capacitive transient was estimated by comparing the current measured at the holding potential to that determined at the same potential within the ramp, and the difference was subtracted from the entire event. The data was analyzed using an Excel (Microsoft Co.) spreadsheet, and secondary analyses were performed using Sigmaplot (SPSS Inc.).
After the oocyte recovered from impalement, CFTR was stimulated by adding a cocktail containing 10 mM forskolin and 1 mM 3-isobutyl-methylxanthine (IBMX) (Research Biochemicals, Inc.) to the perfusate. For ion substitution protocols, the basic amphibian Ringer's was modified to reduce interference from the endogenous Ca2+-activated Cl channel (Ackerman et al. 1994
; Kowdley et al. 1994
; Machaca and Hartzell 1998
) and contained (mM): 98 Na-anion, 2 K-aspartate, 1 Mg-aspartate, 1.8 Ba-acetate, and 10 HEPES, pH 7.5. Some of the anions used are only available as K salts, in which case the solution was modified to contain 98 K-anion and 2 Na-aspartate; otherwise identical (all salts were from Aldrich Chemical Co., Sigma Chemical Co., or Strem Chemicals, Inc.). Typically, an uninjected, mock-injected, or CFTR-expressing, but unstimulated, oocyte exhibited a background conductance of 0.5–1.5 µS, with a slightly higher conductance at depolarized potentials due to the endogenous Ca2+-activated Cl channel in parallel with an endogenous voltage-gated Ca2+ channel. Three of the substitute anions [Au(CN)2, C(CN)3, and N(CN)2] produced a current in uninjected, mock-injected or CFTR-expressing, but unstimulated, oocytes. If the background current was >2% of the current in the stimulated condition, then the background currents were subtracted. The CFTR-independent currents seen in the presence of Au(CN)2, C(CN)3, and N(CN)2 were completely reversible and the magnitude could be faithfully reproduced with repeated exposure. The CFTR-independent current for Au(CN)2, C(CN)3, and N(CN)2 were moderately sensitive to flufenamic acid (
50% inhibition with 250 µM flufenamic acid), suggesting that at least a portion of the conductance is through an endogenous, anion-selective pathway. A portion of the residual current may be due to partitioning into the bilayer, as suggested by the lipophilicity of these compounds (Finkelstein and Cass 1968
). There were no apparent "toxic" effects as judged by the ability of the CFTR conductance to completely recover after exposure to Au(CN)2. Ag(CN)2 appeared to be somewhat toxic to the oocytes and, therefore, its permeation was not characterized.
Calculation of Permeability Ratios, Relative Barrier Height, Relative Binding, and Relative Well Depth
Permeability ratios were calculated using the Goldman-Hodgkin-Katz equation (
) as follows:
![]() | (2) |
Erev is the measured shift in zero current potential after Clo is replaced with the substitute ion, subo, [Cl]'o is the bath concentration of Cl, [Cl]o is the residual Cl in the substituted solution, [sub]o is the concentration of the substitute ion, and R, T, z, and F have their usual meaning. The central goal of the analysis presented here was to use measurements of relative anion permeability and relative anion blockade to estimate the energies associated with transferring an anion from water into the channel. This required that we adopt a model (or models) for the anion translocation process that would permit us to estimate the energetic significance of differences in anion permeability or binding. It is important to note that the primary aim was not to arrive at absolute values for these energies, but rather to determine the trend in the change in the energies from one anion to another so that this trend could be compared with the change in anion size.
The relation used to interpret permeability ratios in terms of energy differences can be obtained using either of two complimentary approaches. Rate theory models for permeation predict that permeability ratios are determined by the difference in peak height for the two ions. For example, for Cl and a substitute anion the permeability ratio is given by
:
![]() | (3) |
(
G)peak is the difference in free energy between Cl and the substitute ion. This calculation applies to any number of barriers as long as the difference in barrier height,
(
G)peak, is the same for all barriers (peak energy offset condition; Hille 1975a
A second route to the relation described by
is to use a lumped Nernst-Planck model in which the permeability is expressed (
):
![]() | (4) |
x is the length of the channel. If it is assumed that the long-range, anion–channel interaction is reflected in the apparent value of β; and that D is approximately equal for different anions, then (
![]() | (5) |
(
G)eq is the difference in the equilibrium transfer free energy between Cl and the substitute ion. In the discussion, a simple model is used to show how these two approaches, one focusing on peak energy and the other on an equilibrium transfer energy, can yield the same result; i.e.,
(
G)peak =
(
G)eq.
Blocking efficacy was determined by exposing the oocyte to a 5-mM concentration of the substitute anion (some of the salts are only available as K salts and, for those experiments, a 5-mM KCl control was added to the protocol). Percent block was determined by measuring the decrease in the slope conductance (+/– 10 mV) at the reversal potential. Each anion was characterized by a half-maximal inhibition constant, K i1/2, calculated by assuming blockade to be a unimolecular binding event that can described by Michaelis-Menten kinetics (Mansoura et al. 1998
) (
):
![]() | (6) |
![]() |
![]() | (7) |
(
G)well is the difference in well depth between the two ions,
GAwell –
GBwell.
We note two qualifications as regards this approach to determining relative well depth. First, the blocking anions used were all permeant to varying extents and were expected, therefore, to contribute to the measured current. Errors due to permeation of the blocking anion were minimized by using a blocker concentration (5 mM) that was
5% that of Cl (105 mM). Second, multiple anion occupancy of CFTR (Tabcharani et al. 1993
) could in principle affect the estimation of relative peak heights and well depths, but under the conditions used in the present experiments (external substitution), only monotonically decreasing conductance was seen with increasing concentration of the blocking anion (see also Mansoura et al. 1998
), so we assumed that channels were singly occupied at all times.
| RESULTS |
|---|
|
|
|---|
(
G)peak so that each number this column reflects the increase in peak height seen by each anion over that seen by C(CN)3. The table also contains the K1/2 for blockade of Cl currents for each anion and the ratio of each K1/2 to that measured for C(CN)3. The latter values were used to calculate
(
G)well, which represents the change in well depth for each anion relative to C(CN)3, so that positive values indicate a shallower well and negative values a deeper well. Also tabulated in Table are values for logKpot, obtained as described in MATERIALS AND METHODS from anion substitution protocols using the PVC-TDMAC membrane and expressed relative to Cl and to C(CN)3. As with permeability ratios and inhibitory constants, the values for logKpot were used to calculate
(
G)trans (see ), which represents the increase in transfer free energy from water to the synthetic membrane for each of the anions with respect to C(CN)3.
|
0.4–5 kJ/mol. As expected (Eisenman and Horn 1983
In Fig. 2, the energy differences associated with anion permeation (A) and block (B) are plotted versus the corresponding energies derived from the response of the PVC-TDMAC membrane. The high correlation of these values (with the exception of iodide, see below) indicates that the selectivity pattern exhibited by CFTR, as judged by either relative permeability or relative binding, is qualitatively identical to that of the synthetic membrane, differing in each case only by a multiplicative constant. Anions that see a barrier height that is increased relative to that of C(CN)3 also experience a more positive (less favorable) transfer free energy between water and the synthetic membrane. Similarly, anions that bind less tightly than C(CN)3 are those for which the water-synthetic membrane transfer free energy is less favorable. It is apparent from Fig. 2 that the peak and well energies change in a parallel fashion. SCN, for example, sees an energy barrier to entering the CFTR channel that is lower than that of Cl, and also sees an equilibrium free energy associated with partitioning into the synthetic membrane that is more favorable than that of Cl. Similarly, the tighter binding of SCN (relative to Cl) is correlated with ease of partitioning into the synthetic membrane. These results are, perhaps, not surprising in that the selectivity patterns for both CFTR and the PVC-TDMAC membrane have both been previously identified as being consistent with the "lyotropic" or Hofmeister series, which is ordered according to relative free energy of hydration (Anderson et al. 1991
; Linsdell et al. 1997a
,Linsdell et al. 1997b
; Tabcharani et al., 1997b). Anions that are more readily dehydrated than Cl exhibit higher permeability ratios and bind more tightly within the CFTR pore, and also partition more readily into the PVC-TDMAC membrane.
|
(
G)trans (see ) are plotted versus the reciprocal of ionic radius (Table ). In Fig. 3 B, the relative heights of the energy barriers associated with entering the CFTR channel obtained from permeability ratios (Table ) are plotted versus the reciprocal of the anionic radius. Because C(CN)3, the largest and most permeant ion, was chosen as the reference anion for both plots, for each anion, either the increase in equilibrium transfer energy (synthetic membrane) or the increase in barrier height (CFTR) relative to that seen by C(CN)3 is plotted versus 1/r. In both cases, the energy difference increases linearly with 1/r. Lewis and Stevens 1983
|
(
Gtrans), associated with either differential anion partitioning into the synthetic membrane (see ) or the barrier to entry into the channel, can be written as in
![]() | (8) |
(
Ghyd) is the relative hydration energy and
(
Gsolv) the relative solvation energy in the membrane, both calculated using C(CN)3 as a reference (i.e.,
![]() | (9) |
![]() | (10) |
GhydA,
GhydCCN3,
GsolvA, and
GsolvCCN3are the energies of hydration and solvation, respectively for anion, A, and C(CN)3 expressed with respect to a vacuum phase.
The solid lines in Fig. 3A and Fig. B, represent the values of |
(
Ghyd)| relative to C(CN)3 calculated as the Born charging energy (see ) using corrected radii as described in MATERIALS AND METHODS (Table ). Note that, although the values of
(
Ghyd) are negative, they are plotted here as absolute values for convenience of comparison with
(
Gtrans). The hydration component of the energy difference will be identical for any channel or membrane, as the necessity to remove an ion from bulk water is a universal constant for any permeation process (Krasne and Eisenman 1973
; Eisenman and Horn 1983
).
The linear, monotonic behavior of
(
Gtrans) and
(
Ghyd) suggested a straightforward calculation of
(
Gsolv) as the difference between these two functions, and the predicted behavior of the solvation energy is indicated by the dotted lines in Fig. 2A and Fig. B. The selectivity patterns exhibited by the PVC-TDMAC membrane and CFTR can be readily understood in terms of the differences in the relative interaction energy of the anions with water,
(
Ghyd), and with the membrane,
(
Gsolv). In the synthetic membrane, the work required to transfer an anion from water to the plasticized PVC decreases with increasing anion radius because
(
Ghyd) is a steeper function of 1/r than is
(
Gsolv). In other words, anions larger than Cl experience weaker interactions with water and with the synthetic membrane, but they partition into the membrane more readily because they see the smallest difference between these two energies.
The linear relation between the apparent solvation energy for the PVC-TDMAC membrane and reciprocal anionic radius (Fig. 3 A) suggests that
Gsolv, the anion-membrane interaction energy for the PVC-TDMAC membrane, behaves precisely as predicted by the Born energy (see
) for a spherical anion contained within a polarizable medium having a dielectric constant somewhat less than that of water. The slope of the plot predicts an effective dielectric constant,
eff, for the synthetic membrane of 4.1. The measured dielectric constant of the PVC-TDMAC membrane was
= 11.5 ± 0.6, as compared with a published value of
= 14 for o-NPOE plasticized PVC membranes constructed without TDMAC (Armstrong and Horvai 1990
). In other words, the synthetic membrane is somewhat more selective than would be predicted by the bulk dielectric constant determined from the apparent capacitance. Fig. 3 A (inset) provides some appreciation of the impact on
Gsolv of changes in the dielectric constant.
Gsolv for a 1-Å sphere increases rapidly as
increases from 1 to
20, but changes are minimal for
> 20. From the perspective of the continuum analysis of solvation energy, a medium having a dielectric constant >20 is effectively as polarizable as water (
= 80).
Fig. 3 B shows the behavior of
(
Gpeak),
(
Ghyd), and
(
Gsolv) for CFTR. It is immediately apparent from Fig. 3 B that the modest permeability selectivity of CFTR can be attributed to the fact that the energies of hydration and solvation differ very little over the range of anion sizes examined. In other words, a visiting anion is solvated within the CFTR pore nearly as well as it is in bulk water. Accordingly, the solvation energy predicts an effective dielectric constant within the pore of
19. The near identity of the value of
Ghyd and
Gsolv justifies treating the energies associated with anion entry as a near equilibrium process. The point is made more explicitly in Fig. 4, in which are shown the predicted values expressed with respect to a vacuum reference phase for
Gpeak,
Ghyd, and
Gsolv, calculated using a value of 19 for the effective dielectric constant within the channel. This plot predicts a peak energy for Cl of 14.5 kJ/mol (5.86 RT), which agrees well with the values derived by Linsdell et al. 1997a
for a multi-site model, and is about half the value of 27.8 kJ/mol (11.2 RT) predicted from a symmetric two-barrier, one-site model using a well depth of 8.2 kJ/mol (3.3 RT) (based on the dissociation constant of 38 mM for chloride as determined by Tabcharani et al. 1997
) and constraining the single channel conductance to 10 pS (Dawson et al. 1999
). Although these values are likely to be significantly affected by the ambiguity as to the appropriate value for the prefactor in the Eyring rate equations (Andersen 1999
; Levitt 1999
; Nonner et al. 1999
), it is apparent that the absolute barrier heights predicted from continuum analysis fall in a range consistent with observed transport rates. The error in the prefactor incurred by using kT/h has been estimated, using a discrete approximation to a continuum model, to be of the order of 102, which would translate into an error in the calculated value of
G of
4.6 RT (Andersen and Koeppe 1992
).
|
12.4, because the selectivity is somewhat higher. The absolute barrier for Cl is predicted to be 19.6 kJ/mol (7.9 RT), which agrees fairly well with the estimates of Bormann et al. 1987
|
6.5 kJ/mol, which agrees well with the value of 8.2 kJ/mol calculated from the dissociation constant of 38 mM reported by Tabcharani et al. 1997
|
The dotted line in Fig. 6 B represents the anion–channel interaction energy derived from relative permeability measurements taken from Fig. 4, plotted for comparison and to emphasize an important point, namely that the radius-dependent portion of the anion-channel interaction energy is similar, regardless of whether it is defined by relative permeability (peak height) or relative blocking affinity (well depth). This plot shows why larger anions experience not only a reduced barrier to entering the channel, but also a deeper energy well within the channel.
Iodide Permeation and Block
In terms of permeability selectivity, iodide stands out as an anomaly. In Fig. 2 A, it can be seen that the value of PI/PCl determined for human wtCFTR expressed in Xenopus oocytes is well below that predicted for an ion that is easier to dehydrate than Cl. Tabcharani et al. 1997
provided evidence that CFTR can exhibit a higher value of PI/PCl (2.1; Fig. 2 A, asterisk), but that the ion causes a rapid modification of CFTR that leads to the lower value of 0.2–0.4 that is most often reported (Anderson et al. 1991
; Sheppard et al. 1993
; Mansoura et al. 1998
). The higher value of PI/PCl would place iodide in its predicted position on the 1/r plot (Fig. 2 A). Using human wtCFTR expressed in Xenopus oocytes, we have not detected any evidence of this higher PI/PCl using a panel of various mole fractions of I:Cl, as well as various voltage clamping protocols (data not shown). It may be, however, that the conversion is simply too fast to resolve in this experimental setting. In accord with the hypothesis of Tabcharani et al. 1997
that I is inducing a modification of CFTR, we have observed a block of CFTR by I that is qualitatively distinct from the block seen with the other permeant ions (see Mansoura et al. 1998
). Block by a small amount of I in the external bath is weakly voltage dependent in the negative quadrant, as shown in Fig. 7, whereas the block seen with most permeant ions is largely voltage independent (Mansoura et al. 1998
). The efficacy of block is enhanced if the external concentration of Cl is lowered (by substitution with aspartate), as shown in Fig. 7 B, suggesting that Cl and I ions are competing for the same site(s). Iodide has a tendency to form interhalogens, particularly triiodide, a reactive species known to attack cysteine, tyrosine, and histidine residues within proteins, but the addition of sodium thiosulfate, which reduces I3, did not alter the observed effects. Furthermore, increasing concentrations of I3 (made by adding I2 to NaI solution) appeared to be quite toxic above 500 µM I3, leading to uncontrollable increases in oocyte conductance. While in solutions near or below 100 µM, the I3 appeared to be unstable in our standard amphibian Ringer's (as judged colorimetrically by clarification of the yellowish-red color of the I3), it was possible to make a fresh I3 solution and feel relatively confident that the oocyte was seeing 10–100 µM I3, in which case the effect on CFTR was an irreversible, voltage-independent block. However, we cannot rule out the possibility that I3 may be forming within the lumen, leading to a mild chemical modification of CFTR. Data from mutagenesis experiments support the view that iodide permeation is most properly thought of as a special case in that mutations that affect PI/PCl do not appear to dramatically alter the permeability ratios for other permeant ions (Anderson et al. 1991
; Tabcharani et al. 1997
; Mansoura et al. 1998
).
|
| DISCUSSION |
|---|
|
|
|---|
The analysis presented here provides a physical basis for "weak field strength" selectivity. In the continuum electrostatic model, the field strength is that due to the polarization of dipolar entities that are subject to the field of the visiting anion. The magnitude of the ion–channel interaction energy depends on the size of the anion and the effective dielectric constant experienced by the anion when it resides within the channel. It is important to point out that the value of the effective dielectric constant would not be expected to be governed solely by the properties of amino acid side chains that might line the pore. Although the polarizability of such entities would contribute to
eff, other contributions would be expected from the remainder of the protein, including side chains and the peptide backbone, water molecules that may reside within the channel, the surrounding lipid, and even water bathing the membrane. The term "effective dielectric constant" embraces this notion (Sham et al. 1998
). It is clear from this analysis that the implication of the term "weak field strength" is not that the anion–channel energies are small. They are, in fact, predicted to be approximately equal to the anion–water energies, but to increase slightly less steeply with decreasing anion size.
The Effective Dielectric Constant and the Solvation of Permeant Ions
The Born-type model employed here is, from both a conceptual and computational perspective, the simplest approach to accounting for the apparent solvation energies of anions that traverse the pore. The model represents the heterogeneous ensemble of components that comprise the environment of an anion as an equivalent continuum of infinite extent characterized by an effective dielectric constant,
eff. The most important result of the analysis is not the value of
eff, however. Rather, it is the fact that the lyotropic selectivity pattern can be predicted by presuming that anions interact with the channel much as they do with water, such that the stabilization energy is linearly related to reciprocal anion radius.
Attempts to measure or predict the value of the apparent dielectric constant seen by an anion inside a channel have produced widely varying results. In a large bore (minimum diameter
10 Å), anion-selective porin (phoE, for example), Gutman et al. 1992
used time-resolved fluorometry to determine a value of
eff for the pore of 24, whereas Karshikoff et al. 1994
predicted a value for the same channel of 30 using a macroscopic, multidielectric model. The analysis conducted by Sham et al. 1998
and Sansom et al. 1997
, however, suggests that the definition of a unique
eff may be elusive because the value is likely to depend on the type of electrostatic interaction used to define it. As indicated in MATERIALS AND METHODS, the value of
eff could also depend on the approach taken to define the radius of the ion (Grunwald 1997
). It is probably most prudent to consider
eff, determined here, as an empirical parameter that provides a measure of the relative ability of the channel to solvate a permeating ion that offers a useful first approach to conceptualizing the electrostatic origins of the selectivity pattern.
The behavior of ions in physical or biological systems has been analyzed previously using an approach similar to the one adopted here. An analysis of the swelling and shrinking of gels by Büchner et al. 1932
led to a relationship between the action of salts on colloidal systems and reciprocal ion radius not unlike that used here. In addition, Lewis and Stevens 1983
presented essentially similar findings for cation permeation in the acetylcholine receptor. The Born analysis, based on their permeability ratios, would predict
eff = 55. More recently, Aqvist and Warshel 1989
used two methods, free energy perturbation and protein dipoles Langevin dipoles simulations, to analyze the polarizable interior of the gramicidin channel to predict the energy barrier to permeation for Na, deriving values of 9.6 (3.9 RT) and 16.7 (6.7 RT) kJ/mol, respectively. Using the data of Urban et al. 1980
, the analysis used here would predict an
eff of 21.5–27 for gramicidin A, and an absolute barrier for Na entry of 10–14 kJ/mol. Partenskii et al. 1994
reported that, within the framework of a three-dielectric, continuum model in which the dielectric constant of the pore region was restricted to between 2 and 5, it was difficult to account for the apparent energy barrier that characterizes cation flow through the gramicidin channel.
Structural Implications for Cl-selective Channels
As regards the structure of the CFTR pore domain, the most important implication of the analysis presented here is the prediction that the permeability selectivity pattern characteristic of CFTR, as well as several other nonhomologous anion channels, does not depend on the interaction of the permeating anion with some specific component of the channel. The basic selectivity pattern, common to anion permeability and anion binding, can be viewed as the result of the interaction of the anion with a volume that exhibits the generalized property of dielectric polarizability. This type of environment could presumably be provided by a configuration of the membrane spanning segments of CFTR which, along with the resident water molecules, forms a polarizable "tunnel" through which the anions can pass (Sansom et al. 1996
, Sansom et al. 1997
). The present results do not speak to the location of any "selectivity filter" and, in fact, suggest that the very idea that the anion permeability selectivity of CFTR is associated with any specific structure, as in the K channel visualized by Doyle et al. 1998
, may be moot. The barrier to anion entry (or the energy associated with equilibrium partitioning) is predicted to be a generalized feature of this dielectric tunnel that is not highly dependent on the details of the structure of the transmembrane segments. This prediction has, in fact, been borne out by the results of mutational analysis, which have shown that permeability ratios are largely insensitive to mutations (Anderson et al. 1991
; Mansoura et al. 1998
). In addition, the observation that the same analysis predicts the relative pattern for anion channels that are likely to be divergent in their structure lends further support to the notion that the anion–channel interaction energy is the result of a generalized interaction of the anion with the pore.
The selectivity of anion binding by CFTR differs from relative permeability in at least two ways. First, in contrast to permeability, binding is highly sensitive to point mutations, particularly in transmembrane segments 5 and 6 (Anderson et al. 1991
; Tabcharani et al. 1993
; Linsdell et al. 1997b
; Smith et al. 1997
; Mansoura et al. 1998
). Second, the analysis of energetics presented here suggests that, in order to predict the selectivity pattern associated with anion binding, the underlying anion–channel interaction energy must be envisioned as comprising two components: one due to the same radius-dependent, dielectric stabilization that appears to determine barrier height, and a second that is similar for all anions, roughly independent of size. These two attributes, sensitivity to structural perturbations and two-component energetics, are easily reconciled by assuming that anion binding, represented by the energy well in the rate theory description, is dependent on the existence within the conduction path of a narrow region in which the residues that line the channel are able to come into a more intimate, inner sphere and contact the visiting ion. If this region were relatively rigid so that it could not change conformation in the presence of different permeant ions, then the energy associated with residing there could be roughly independent of anion size. In contrast, the contribution to the anion–channel energy due to the region of the protein surrounding the cavity would exhibit the dependence on anion size predicted for a polarizable medium. Mutations could alter the nature of this cavity by changing its size or, equally likely, by changing the charge or orientation of amino acid side chains that are required for a favorable anion–channel interaction. Reduction or elimination of the anion–channel contact would leave the underlying, radius-dependent anion–channel energy intact so that permeability ratios would not be greatly affected even though anion binding was nearly lost.
A Working Model for Anion Permeation and Binding
The analysis of permeation and binding energetics undertaken here provides the basis for a working model of the anion conduction process in CFTR and perhaps the GABAR, GlyR, and T84-ORCC as well. In the case of CFTR, it is possible to envision two sorts of CFTR pores: those that bind anions, exemplified by the wild-type channel, and those that do not, exemplified by mutant CFTRs like G314E or Q (Mansoura et al. 1998
) and R347D (Tabcharani et al. 1993
). This structural dichotomy suggests that it is reasonable to treat these two facets of the conduction process as representing two distinct types of physical interaction of the anion with the channel.
Consider first a channel that does not bind anions (Fig. 8 A). It is useful to envision permeant anions in the bulk solution as coordinated by an inner sphere of water molecules and surrounded by an outer sphere or shell that is the remainder of the bulk solution (Sharp and Honig 1990
; Andersen and Koeppe 1992
; Marcus 1994
). Upon colliding with the mouth of the channel, the anion and channel form a transition state complex that leads to the anion, along with most of its inner shell water, residing within the channel. In this state, the channel, to a first approximation, has replaced the outer shell waters (i.e., the bulk solution), but the inner sphere waters remain associated with the anion. The energy profile expected for the process is depicted as the trapezoidal line in Fig. 8 D and represents an equilibrium transfer free energy (reflected in the partition coefficient, β) that, for Cl, would be unfavorable by
6 RT (see RESULTS) due to the fact that the effective dielectric constant of the channel protein, the lipid membrane, and surrounding water is somewhat less than that of the outer sphere of water molecules in bulk solution. This represents the energy barrier seen by an anion entering the channel.
|
11.4 kJ/mol (4.6 RT) that is of the same order as that seen in the weakly selective channels considered here. If these two profiles are summed (Fig. 8 C), they give rise to the familiar two-barrier, one-well profile that is depicted as the dotted line in D. Although the shape of the profile as depicted in Fig. 8 D is largely arbitrary, the diagram makes the point that it is possible, in principle, to account for the energetics of anion permeation through CFTR in a relatively straightforward way and illustrates how the summing of an equilibrium transfer energy with a single, localized energy well could give rise to the familiar two-barrier, one-site channel model.
Fig. 8 E illustrates how radius-dependent variation in the outer sphere anion energy could vary the peak heights and well depths in a parallel fashion if the inner-sphere contribution in the binding region of the channel was roughly size independent. Large anions that experience a reduced barrier height, and enter the channel more readily, also see a deeper energy well and bind more tightly because of the reduced, radius-dependent equilibrium transfer energy.
The Significance of Ion Selectivity Patterns: from Bias to Recognition
The permeability selectivity exhibited by CFTR, and shared by GABAR, GlyR, and T84-ORCC, represents the most primitive type of ion discrimination, characteristic of a permeation path that functions as a polarizable tunnel that can stabilize a partially dehydrated ion as it passes through. This situation is a striking contrast to that envisioned for the bacterial K channel that selects for K by means of a clearly identifiable structure, the selectivity filter, consisting of a tetrahedral array of oxygen ligands (Armstrong 1975
, Armstrong 1989
; Doyle et al. 1998
). In a recent review of the physical factors that govern anion separations, Moyer and Bonnesen 1997
suggested the utility of viewing selectivity within the framework of two limiting types. In one, exemplified in biology by the K channel, the host (channel) is structurally specialized to recognize the guest (ion). In the other, here exemplified by CFTR, selectivity is based simply on a physical bias that is imposed by the physics of hydration and a primitive form of solvation of the guest by the host molecule. This approach to categorizing selectivity types provides a useful framework for interpreting selectivity patterns characteristic of different channel families (Halm 1998
). The physical forces that produce a bias-type selectivity must impinge on any ion–channel interaction, but channels may be expected to vary in the extent to which this bias is overshadowed by a specific recognition component. The muscle Cl channel, ClC-1, may be a case in point. It has been reported to exhibit a distinctly nonlyotropic pattern of permeability selectivity (Cl > SCN > Br > NO3 > I) (Fahlke et al. 1997
) that may be an indication of some structural specialization that has evolved to recognize the Cl anion. It is of interest in this regard that mutation of several residues in ClC-1 resulted in a reversion of the selectivity sequence to the more primitive, lyotropic order (Fahlke et al. 1997
). What is striking in the case of CFTR (and, likely, GlyR, GABAR, and T84-ORCC) is that the compound effects of bias-based permeability selectivity and anion binding result in channels that are "optimized" Cl filters.
Permeability and Anion Size
The analysis employed here emphasizes one aspect of the relation of anion size to permeation, namely that halides and pseudohalides having an equivalent sphere radius larger than that of Cl enter the pore more readily due to lower anion-water interaction energies. There is clearly a limit, however, to any "larger is better" theory of permeation. As anion size increases, the physical dimension of the pore, its effective diameter, must become limiting. Several of the "larger" molecules that were the focus of this study are roughly cylindrical in shape and the actual physical diameter of the cylinder is less than that of the diameter of the equivalent sphere (dicyanoaurate, for example, has a cylindrical diameter of
3.4 Å at its widest point, while the equivalent sphere diameter of the molecule is
6 Å; Table ). On the other hand, the trimmest right cylinder into which tricyanomethanide, which has flattened pyramidal geometry, could fit would be
7.4 Å in diameter, its widest dimension, due to the fact that there is no way to "twist" the pyramid to fit it into a smaller cylinder. This may seem to be inconsistent with the effective pore diameter of
5.5 Å determined by Linsdell et al. 1997b
, Linsdell et al. 1998
; however, if one imagines that the shape of the pore is elliptical, such that the widest part of the pore is on the order of 7.5 Å, then the narrow portion of the pore could be on the order of 5.5 Å, and thereby accommodate tricyanomethanide.
It may be necessary, however, to exercise some caution in imputing effective pore size from the behavior of poorly permeant ions. As indicated in Fig. 1, polyatomic molecules like gluconate are characterized by a hydration energy that is much larger than their physical size would predict, an effect that is presumably due in part to a nonuniform charge distribution (Gilson and Honig 1988
; Marcus 1997
). On the basis of Fig. 1, it is possible to assign to gluconate an apparent "thermochemical radius" of
1.65 Å. Using the value of
eff = 19 for the CFTR pore predicts a peak barrier height of 15.9 kJ/mol (6.4 RT), and the difference in this value and that for Cl [14.5 kJ/mol (5.86 RT), see RESULTS] predicts a permeability ratio (Pgluconate/PCl) of
0.58. This may be compared with the experimentally determined values of 0.071 and 0.013 reported by Linsdell and Hanrahan 1998
, determined by substitution of the solution on the cytoplasmic and extracellular side of the patch, respectively. The larger of these two values (0.071) predicts a peak energy of 21.1 kJ/mol, 5.2 kJ/mol greater than the simple electrostatic model. This increased barrier height could point to a size exclusion effect. Marcus 1997
, however, points out that asymmetric charge distribution could have other consequences that are not directly related to size, because the polar and nonpolar portion of the molecule can experience very different interactions with the immediate environment.
| An Electrostatic Analysis of the Selectivity of the Plasticized PVC-TDMAC Membrane Anion-selective membranes were constructed of a mixture of 66% (by weight) o-NPOE, 33% PVC, and 1% TDMAC. Within the 100-µm-thick membrane, TDMAC functions as a completely dissociated anion exchanger that favors the extraction of mobile anions and excludes cations by a Donnan mechanism (Teorell 1953
We treat the potential across this "thick" membrane as comprised of three components, two phase boundary potentials,
Anion equilibrium between the aqueous and membrane phases is defined by the equality of the anion chemical potential in the two phases (
1, is given by
The first term on the right hand side represents the contribution to
The second term in
Outside of the space charge region, electroneutrality dictates that the concentration of the mobile anions be equal to that of the dissociated ionophore, as shown in
1, is expressed as the sum of two components: one the Donnan potential and the other due to the differential partitioning of the anion between the aqueous and membrane phases. Now, we consider the change in 1, ![]() 1, that occurs when a test anion, A, is substituted for the reference anion, Cl. For a complete equimolar replacement of Cl by S, the first term in
![]() 1 is given by (
In the analytical chemistry literature (Morf 1981
Recalling
µ*A and µ*Cl represent the difference in standard chemical potential between the aqueous and membrane phases for the test anion, A and Cl, respectively.
Each of these differences in standard chemical potential can be expressed as the difference in the free energy associated with bringing the anion from a reference phase (defined as a vacuum) to the aqueous phase or to the membrane phase, respectively, so that for the substitute ion and Cl,
For any anion, the difference in standard chemical potential, and hence the partition coefficient, may be regarded as a reflection of the difference between the energy of interaction of the anion with water,
An Electrostatic Model for
Each of these phases is represented as being an infinite, structureless continuum, characterized by a dielectric constant,
is the dielectric constant, and K is a constant given by
0 is the permittivity of free space.
|
Dr. Smith's and Dr. Dawson's present address is Department of Physiology & Pharmacology, Oregon Health Sciences University, Portland, OR 97201
Abbreviations used in this paper: CFTR, cystic fibrosis transmembrane conductance regulator; GABAR, gamma amino butyric acid receptor; GlyR, glycine receptor; o-NPOE, ortho-nitrophenyloctyether; ORCC, outwardly rectifying chloride channel; PVC, poly(vinyl chloride); TDMAC, tridodecylmethylammonium chloride; wtCFTR, wild-type CFTR.
Here we assume that we may treat anion entry into CFTR as a quasi equilibrium process, as implied in the Nernst-Planck interpretation of the permeability ratio (see MATERIALS AND METHODS and DISCUSSION).
| ACKNOWLEDGMENTS |
|---|
This work was supported by grants from the National Institutes of Health (DK45880 to D.C. Dawson and GM-28882 to M.E. Meyerhoff), the University of Michigan G.I. Peptide Center, and The Center for Membrane Toxicity Studies at the Mount Desert Island Biological Laboratory.
Submitted: 1 April 1999
Revised: 29 October 1999
Accepted: 1 November 1999
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Fatehi and P. Linsdell State-dependent Access of Anions to the Cystic Fibrosis Transmembrane Conductance Regulator Chloride Channel Pore J. Biol. Chem., March 7, 2008; 283(10): 6102 - 6109. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. D. Parker, M. T. Young, C. M. Daly, R. W. Meech, W. F. Boron, and M. J. A. Tanner A conductive pathway generated from fragments of the human red cell anion exchanger AE1 J. Physiol., May 15, 2007; 581(1): 33 - 50. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. G. Palmer and G. Frindt Cl- channels of the distal nephron Am J Physiol Renal Physiol, December 1, 2006; 291(6): F1157 - F1168. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Liu, C. Alexander, J. Serrano, E. Borg, and D. C. Dawson Variable Reactivity of an Engineered Cysteine at Position 338 in Cystic Fibrosis Transmembrane Conductance Regulator Reflects Different Chemical States of the Thiol J. Biol. Chem., March 24, 2006; 281(12): 8275 - 8285. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Linsdell Mechanism of chloride permeation in the cystic fibrosis transmembrane conductance regulator chloride channel Exp Physiol, January 1, 2006; 91(1): 123 - 129. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Liu, D. Kozono, Y. Kato, P. Agre, A. Hazama, and M. Yasui From The Cover: Conversion of aquaporin 6 from an anion channel to a water-selective channel by a single amino acid substitution PNAS, February 8, 2005; 102(6): 2192 - 2197. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Qu and C. Hartzell Determinants of Anion Permeation in the Second Transmembrane Domain of the Mouse Bestrophin-2 Chloride Channel J. Gen. Physiol., September 27, 2004; 124(4): 371 - 382. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Ryan, A. D. Mitrovic, and R. J. Vandenberg The Chloride Permeation Pathway of a Glutamate Transporter and Its Proximity to the Glutamate Translocation Pathway J. Biol. Chem., May 14, 2004; 279(20): 20742 - 20751. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Qu, R. Fischmeister, and C. Hartzell Mouse Bestrophin-2 Is a Bona fide Cl- Channel: Identification of a Residue Important in Anion Binding and Conduction J. Gen. Physiol., March 29, 2004; 123(4): 327 - 340. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Gong and P. Linsdell Mutation-induced Blocker Permeability and Multiion Block of the CFTR Chloride Channel Pore J. Gen. Physiol., November 24, 2003; 122(6): 673 - 687. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Gong and P. Linsdell Coupled movement of permeant and blocking ions in the CFTR chloride channel pore J. Physiol., June 1, 2003; 549(2): 375 - 385. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Gong and P. Linsdell Molecular determinants and role of an anion binding site in the external mouth of the CFTR chloride channel pore J. Physiol., June 1, 2003; 549(2): 387 - 397. [Abstract] [Full Text] [PDF] |
||||
![]() |
V E Wotring, T S Miller, and D S Weiss Mutations at the GABA receptor selectivity filter: a possible role for effective charges J. Physiol., April 15, 2003; 548(2): 527 - 540. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Zhou, S. Hu, and T.-C. Hwang Probing an Open CFTR Pore with Organic Anion Blockers J. Gen. Physiol., October 29, 2002; 120(5): 647 - 662. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Van der Kloot, J. Molgo, R. Cameron, and C. Colasante Vesicle size and transmitter release at the frog neuromuscular junction when quantal acetylcholine content is increased or decreased J. Physiol., June 1, 2002; 541(2): 385 - 393. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Gong, S. M Burbridge, E. A Cowley, and P. Linsdell Molecular determinants of Au(CN)2- binding and permeability within the cystic fibrosis transmembrane conductance regulator Cl- channel pore J. Physiol., April 1, 2002; 540(1): 39 - 47. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Linsdell and X. Gong Multiple inhibitory effects of Au(CN)2- ions on cystic fibrosis transmembrane conductance regulator Cl- channel currents J. Physiol., April 1, 2002; 540(1): 29 - 38. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. S. Smith, X. Liu, Z.-R. Zhang, F. Sun, T. E. Kriewall, N. A. McCarty, and D. C. Dawson Cftr: Covalent and Noncovalent Modification Suggests a Role for Fixed Charges in Anion Conduction J. Gen. Physiol., October 1, 2001; 118(4): 407 - 432. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. A. McCarty and Z.-R. Zhang Identification of a region of strong discrimination in the pore of CFTR Am J Physiol Lung Cell Mol Physiol, October 1, 2001; 281(4): L852 - L867. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Zhou, S. Hu, and T.-C. Hwang Voltage-dependent flickery block of an open cystic fibrosis transmembrane conductance regulator (CFTR) channel pore J. Physiol., April 15, 2001; 532(2): 435 - 448. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Linsdell Relationship between anion binding and anion permeability revealed by mutagenesis within the cystic fibrosis transmembrane conductance regulator chloride channel pore J. Physiol., February 15, 2001; 531(1): 51 - 66. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Y Rychkov, M. Pusch, M. L Roberts, and A. H Bretag Interaction of hydrophobic anions with the rat skeletal muscle chloride channel ClC-1: effects on permeation and gating J. Physiol., February 1, 2001; 530(3): 379 - 393. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. O'Connor and K. Beyenbach Chloride channels in apical membrane patches of stellate cells of Malpighian tubules of Aedes aegypti J. Exp. Biol., January 1, 2001; 204(2): 367 - 378. [Abstract] [PDF] |
||||
![]() |
Z. Qu and H. C. Hartzell Anion Permeation in Ca2+-Activated Cl- Channels J. Gen. Physiol., December 1, 2000; 116(6): 825 - 844. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. McCarty Permeation through the CFTR chloride channel J. Exp. Biol., January 7, 2000; 203(13): 1947 - 1962. [Abstract] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|