|
||
Original Article |
laurent.schild{at}ipharm.unil.ch
| ABSTRACT |
|---|
|
|
|---|
, one β, and one
homologous subunits, each subunit comprising two transmembrane segments. Amino acid residues involved in binding of the pore blocker amiloride are located in the pre-M2 segment of β and
subunits, which precedes the second putative transmembrane
helix (M2). A residue in the
subunit (
S589) at the NH2 terminus of M2 is critical for the molecular sieving properties of ENaC. ENaC is more permeable to Li+ than Na+ ions. The concentration of half-maximal unitary conductance is 38 mM for Na+ and 118 mM for Li+, a kinetic property that can account for the differences in Li+ and Na+ permeability. We show here that mutation of amino acid residues at homologous positions in the pre-M2 segment of
, β, and
subunits (
G587, βG529,
S541) decreases the Li+/Na+ selectivity by changing the apparent channel affinity for Li+ and Na+. Fitting single-channel data of the Li+ permeation to a discrete-state model including three barriers and two binding sites revealed that these mutations increased the energy needed for the translocation of Li+ from an outer ion binding site through the selectivity filter. Mutation of βG529 to Ser, Cys, or Asp made ENaC partially permeable to K+ and larger ions, similar to the previously reported
S589 mutations. We conclude that the residues
G587 to
S589 and homologous residues in the β and
subunits form the selectivity filter, which tightly accommodates Na+ and Li+ ions and excludes larger ions like K+.
Key Words: epithelial Na+ channel Xenopus oocyte pore selectivity ion channel
| Introduction |
|---|
|
|
|---|
ENaC belongs to a new class of channel proteins called the ENaC/DEG superfamily, which includes a variety of proteins involved in mechanotransduction and neurotransmission, and is found in nematodes, flies, snails, and mammals (for review see Tavernarakis and Driscoll 1997
; Waldmann and Lazdunski 1998
). Subunits of this superfamily coassemble usually within subfamilies into Na+-preferring or -selective multimeric channels that are either constitutively active (e.g., ENaC), activated by mechanical stimuli (as postulated for Caenorhabditis elegans degenerins), activated by a peptide (FMRFamide peptide–gated Na+ channel, FANaC; Lingueglia et al. 1995
), or activated by protons (ASIC).
ENaC is a heterotetramer, and made of two
, one β, and one
homologous subunits arranged around the channel pore in an
β
configuration (Canessa et al. 1994b
; Firsov et al. 1998
). Each homologous subunit has two transmembrane spanning segments (M1 and M2) with intracellular NH2 and COOH termini leaving a large extracellular hydrophilic loop, as illustrated in Fig. 1 A (Canessa et al. 1994a
; Renard et al. 1994
). Based on sequence comparisons, current models predict that the second transmembrane spanning segment of ENaC forms an
helix starting with a conserved Ser residue (
S589) and extending 22 residues further downstream in the defined ENaC sequence (see Fig. 1 B). A pre-M2 segment can arbitrarily be defined as a sequence delineated by a conserved Gly (
G579) residue at the 5' end and
S589 at the 3' end that initiates the M2 segment. Fig. 1 C shows a model of the narrow pore region of ENaC, based on this and previous work. Previous mutagenesis experiments provided evidence that the pre-M2 forms the outer pore of ENaC (Schild et al. 1997
). Amino acid residues mutated in those experiments are presented in bold on a gray background in Fig. 1 B. These experiments showed that mutations of Gly residues in β (G525) and
(G537) subunits decrease the affinity for the pore blocker amiloride and change single-channel conductance. In addition, a Cys substitution at the homologous position of
ENaC (
S583C) generated a high affinity Zn2+ binding site that leads to channel block by Zn2+. Recently, we have shown that mutations of the conserved Ser residue
S589 allow larger ions such as K+, Rb+, and Cs+ as well as divalent cations to pass through the channel (Kellenberger et al. 1999
). Ion substitution experiments indicate that
S589 determines the molecular cutoff of the channel pore at its narrowest point, the selectivity filter. Thus, in the pre-M2 segment that lines the outer channel pore, the amiloride binding site precedes the selectivity filter in the sequence.
|
S583) and the selectivity filter (
S589) in channel permeation properties and blocking by amiloride. We have focused our interest on the conserved aromatic residues at the position corresponding to W585 in the
rENaC sequence, and on a cluster of Ser and Gly residues (
G587,
S588, and analogous amino acids in β and
subunits; see Fig. 1 B). | Methods |
|---|
|
|
|---|
, β,
subunit were synthesized in vitro. For binding experiments,
, β, and
subunits that had been tagged as described by Firsov et al. 1996
β
ENaC subunits at a total concentration of 100 ng/µl. For simplicity, ENaC mutants are named by the mutated subunit only, although always all three subunits (
, β, and
) were coexpressed.
Binding Experiments
The FLAG reporter octapeptide, which had been introduced in
, β, and
subunits, is recognized by the anti–FLAG M2 mouse monoclonal antibody (M2Ab) (Eastman Kodak Co.). M2Ab was iodinated as described by Firsov et al. 1996
. Iodinated M2Ab had a specific activity of 5–20 · 1017 cpm/mol and were used up to 2 mo after synthesis. On the day after mRNA injection, oocytes were transferred to a 2-ml Eppendorf tube containing modified Barth's saline (mM: 10 NaCl, 90 N-methyl-D-glutamine HCl, 0.8 MgSO4, 0.4 CaCl2, 5 HEPES, pH 7.2) supplemented with 10% heat-inactivated calf serum, and incubated for 30 min on ice. The binding was started upon addition of 12 nM 125I-M2Ab (final concentration) in a volume of 5–6 µl/oocyte. After 1 h of incubation on ice, the oocytes were washed eight times with 1 ml modified Barth's saline supplemented with 5% heat-inactivated calf serum, and then transferred individually into tubes for
counting containing 250 µl of the same solution. The samples were counted and the same oocytes were kept for subsequent measurement of the whole-cell current. Nonspecific binding was determined from parallel assays of noninjected oocytes. Theoretically, it might be possible that our mutations affect the accessibility of the FLAG epitope for the M2 antibody by changing the conformation of the extracellular loop. However, this possibility seems rather unlikely since most of the mutations did not affect 125I-M2Ab binding.
Electrophysiological Recording
Electrophysiological measurements were taken at 16–20 h after injection. Two-electrode voltage-clamp recordings were obtained using a TEV-200 amplifier (Dagan Corp.). The standard bath solution contained 110 mM NaCl, 1.8 mM CaCl2, 10 mM HEPES-NaOH, pH 7.35. For selectivity measurements, Na+ was replaced by Li+, K+, Rb+, or Cs+ at the same concentration. Macroscopic amiloride-sensitive currents (I) are defined as the difference between ionic currents obtained in the presence and absence of 5 µM (or higher concentrations for some mutants, as indicated) of amiloride (Sigma Chemical Co.) in the bath. All macroscopic currents shown are amiloride-sensitive currents as defined above. Pulses for current–voltage curves were applied, and data were acquired using a PC-based data acquisition system (Pulse; HEKA Electronik).
The cell-attached or outside-out configuration of the patch-clamp technique was used to obtain macropatch and single-channel data. Before recording, the vitelline layer of the oocyte was removed. For cell-attached patches, the oocytes were kept in a standard bath K+ solution to depolarize the membrane potential, and pipette solutions were Na+ or Li+, as described above. For outside-out patches, extracellular solutions were as described above. Changes of external solutions of outside-out patches were made using the Rapid Solution Changer RSC-200 (BioLogic International Ltd.). In this system, the perfusion solutions are driven by gravity to the rotating head containing in these experiments up to 11 glass tubes. The solution exchange is performed by a highly precise and fast rotation of the RSC head, which exposes the patch pipette to the flow of one of the tubes. Times of rotation from one tube to the adjacent one were a few milliseconds. The pipet solution contained 75 mM CsF, 17 mM N-methyl-D-glucamine, 10 mM EGTA, and 10 mM HEPES, pH 7.35. Pipettes were pulled from Borosilicate glass (World Precision Instruments, Inc.). In patch-clamp experiments, currents were recorded with a List EPC-9 patch clamp amplifier (HEKA Electronik) and filtered at 100 Hz for analysis. Data are shown as mean ± SEM, or as indicated.
Data Analysis
To analyze titration curves for inhibition of macroscopic Li+ or Na+ currents (I), the ratio I/I0 measured in the presence (I) of a particular blocker B to that in the absence of the blocker (I0) is described by the inhibition equation:
, where Ki is the inhibitory constant of the blocker, [B] is the concentration of the blocker, and n' is a pseudo–Hill coefficient. For the fit of Li+ block of Na+ current through
S541A ENaC (see Fig. 8), n' was set equal to 1, and a nonblockable fraction of the amiloride-sensitive current (I) was introduced. This nonblockable fraction of I was set equal to the amiloride-sensitive current carried by 140 mM Li+ alone, normalized to the amiloride-sensitive current carried by 20 mM Na+
.
|
The energy diagram in Fig. 9 A (below) summarizes the adjustable parameters of the model. The energies of the unoccupied channel at zero voltage are defined by three peak energies (G1, G2, and G3) and two wells or site energies (U1 and U2) for Na+ and Li+. The subscripts of the parameters refer to the position with respect to the inside solution as shown in Fig. 9 A. The program can deal with the simultaneous presence of three types of ions. The distances D1–D6 refer to the fraction of the electric field that separates peak and well positions with the requirement that
. In addition to the above energy and distance parameters, there is an interaction energy A(ion x, ion y), which models the effects of ion–ion interactions between the same or different types of ions. The reference energy state used in the model is 1 mol fraction, which corresponds to 55.5 M. Conversion to a 1-M reference state is readily accomplished by addition of 4.0 RT units to the reported peak and well energies. The translocation rate constant K for translocation from well U over peak G is related to the translocation barrier height
according to
, where k is Boltzmann's constant, T is the temperature in °K, and h is Planck's constant
.
|
, and this arrangement was used for fitting of all mutants. The quality of the fit was relatively insensitive to the ion–ion interaction parameter A. This can be readily explained by the fact that the data used for constructing the model was obtained at relatively low ion concentrations (typically <160 mM). In this concentration range, the probability of double occupancy of the channel was <0.01, even with A set to a low value (2.6). This is consistent with flux-ratio experiments indicating that ENaC forms a one-ion pore (Palmer 1982
The model obtained for Na+ permeation of wt channels is shown in Fig. 9 B (below). The model fits the conductance/ion concentration relationship well, as shown by the comparison of model prediction (solid line) and data points (symbols) in Fig. 7 (below). The outer well (U2) is deeper than the inner well (U1). This difference allowed fitting of mutant data, in which the apparent channel affinity for Na+ decreased (
S541A,
G587A) by just changing one well energy (one binding site) with regard to wt. This energy profile predicts an almost ohmic behavior of the channel, with outward currents <111% and >90% of inward currents over a voltage range of ±150 mV. Based on the energy profile obtained for the wt Na+ permeation, the profile best suited for Li+ permeation of wt ENaC was obtained by increasing U2 and G3 (see Fig. 7 and Fig. 9 A). Starting from these two models for Li+ and Na+ permeation, we attempted to fit the mutant data by changing only energy values of one well and one barrier at the extracellular side of the permeation pathway, either U2 and G2, or U2 and G3, assuming that single point mutations in the external channel pore would not affect binding sites on both sides of the highest barrier of the selectivity filter. Adjustment of at least one barrier was necessary to adjust the model to the experimental conductance values. By changing U2 and G2, we obtained better fits of our experimental data than by changing U2 and G3. For this reason, we report here models obtained by only the first method.
|
S541A mutant by increasing Li+ concentrations (see Fig. 8). In this case, the model takes into account the simultaneous presence of Na+ and Li+ in the external solution. The use of the 3B2S model requires that macroscopic currents are converted to unitary currents. The macroscopic current data of Li+ block in macropatches were normalized at each voltage to the condition with 20 mM Na+ alone and, for each voltage, these normalized currents were multiplied by the single-channel current for 20 mM Na+ predicted by the
S541A-Na+ model. These modified data were then used to obtain a model for Li+ permeation of the
S541A channel. | Results |
|---|
|
|
|---|
S583 and the corresponding residues βG525 and
G539 line the outer channel pore as shown in Fig. 1 C at a site where amiloride plugs the channel (Schild et al. 1997
S589 at the NH2 terminus of M2 is part of the narrowest region of the pore at the selectivity filter (Kellenberger et al. 1999
W585 and
V590 and their analogues in β and
subunits, are conserved and might be important for the ion permeation properties. We have mutated the
W585, βW527, and
W539 residues to various amino acids and did not observe any significant changes in ion selectivity, unitary conductance, or block by amiloride (Tables I–III). The level of cell surface expression of wt and mutant ENaC was determined on intact oocytes by specific binding of an iodinated monoclonal antibody to a FLAG epitope inserted in the extracellular domain of ENaC subunits (see METHODS). Some mutations resulted in a low expression of the channel at the cell surface (Table ). Similar to Trp mutations, Ala substitutions of the conserved
V590, βV532, and
V544 residues were without effects on the biophysical and pharmacological properties of the channel (Tables I–III). This lack of effects of mutations of the conserved Trp and Val residues suggests that these amino acids may not line the channel pore.
|
G587 and
S588 and the analogous amino acid residues in β and
subunits (Fig. 1 B). The
S588I mutation is known to increase single channel conductance for Na+ ions and to decrease channel affinity for blocking by amiloride (Waldmann et al. 1995b
S588 and the homologous residues βG530 and
G542 did not change ion conductance properties or channel affinity for amiloride (Tables I–III). The effect of
S588I on channel conductance corroborates the previous observation that
S588 is close to the selectivity filter (Kellenberger et al. 1999
Mutations of
G587, βG529, and
S541 Change Ion-permeation Properties and Channel Block by Amiloride
Conserved Gly residues in the
and β subunits at positions
587 and β529, and the unique Ser residue in the
subunit (
S541) were substituted with the polar Ser, the nonpolar residues Cys, Gly, or Ala or the negatively charged Asp. Analysis of the concentration-dependent block of macroscopic Na+ currents by amiloride (Table ) showed that specific mutations in
and β subunits decreased channel affinity for amiloride. The βG529A, βG529S, and βG529C mutations resulted in a 40- to 130-fold increase in Ki for amiloride, and a significant sevenfold increase in amiloride Ki was measured for
G587S. Some substitutions of these residues had no effect on block by amiloride (
G587A, βG529D, βG529R; Table ), suggesting that these amino acids are not part of the amiloride binding site, but rather are indirectly involved in channel block.
No INa was detected for
G587D and
S541N or
S541F mutants, and only a small INa (<0.3 µA) was detected for
S541R even using high (1 mM) amiloride or benzamil concentrations. In the case of the
G587D mutant, INa was not detected despite a 5.5-fold higher anti–FLAG binding signal compared with the background signal. The specific anti–FLAG binding to oocytes expressing the
G587D mutant was only slightly lower than binding to oocytes expressing the functional
G587A mutant (Table ). Similarly, specific anti–FLAG antibody binding was about the same for the nonconducting
S541N and
S541F mutants and for
S541R compared with the functional
S541A mutant (Table ). These observations indicate that
G587D,
S541N,
S541F, and
S541R mutants are expressed at the cell surface, but are nonconducting channels or channels with minuscule conductance (
S541R), consistent with the notion that mutations of
G587 and
S541 residues affect ion permeation through the channel. Thus it appears that ENaC can tolerate relatively conservative substitutions at these positions and that less conservative replacements can make the channel nonfunctional.
We tested whether
G587A, βG529A, or
S541A mutations affect channel selectivity for Li+ over Na+ ions. Representative current traces in the presence of extracellular Na+ solution at different voltages are shown in Fig. 2A and Fig. B, for wt and the
S541A mutant. When Li+ replaced Na+ in the external medium, the amiloride-sensitive inward current increased in oocytes expressing ENaC wt (Fig. 2 A), but clearly decreased in oocytes expressing
S541A (Fig. 2 B). The macroscopic current–voltage relationships of ENaC wt and the
G587A, βG529A, or
S541A mutants in the presence of external Na+, Li+, or K+ ions are plotted in Fig. 2 C. The currents of the I/V curves are normalized to INa of each channel type measured at –100 mV. In external Na+ solution, the I/V behavior of wt and mutant channels was identical. The ILi/INa ratio was lower for the three mutants compared with wt, and followed the sequence wt >
G587A > βG529A >
S541A. The current ratio ILi/INa at –100 mV of all the
G587, βG529, and
S541 mutants is shown in Table . In addition to mutants shown in Fig. 2, βG529C also exhibited a decreased ILi/INa ratio, whereas other mutants such as βG529D, βG529R, or
S541G maintained ILi/INa ratios > 1, similar to ENaC wt. Thus, the different effects by different βG529 substitutions cannot be correlated with properties of the amino acid side chain. The fact that substitution of βG529 by the much larger Arg did not change ion permeation indicates that the side chain may point away from the channel pore (see DISCUSSION).
|
|
G587A, βG529A, and
S541A mutants even at hyperpolarized membrane potentials (Fig. 2 C). However, βG529S generated a significant amiloride-sensitive inward K+ current as illustrated by the current–voltage behavior in Fig. 3. In the presence of external Na+ the I/V relationship of the βG529S mutant was identical to wt. When K+ replaced Na+ ions in the external medium, βG529S channels exhibited measurable inward currents at negative membrane potentials. At –100 mV, IK represented 22% of INa
. The βG529S mutant allows to a lesser extent Rb+ ions to pass through the channel
, but not larger ions like Cs+
. The permeability of the βG529S mutant to K+ and Rb+ suggests that this mutation increased the molecular cutoff of the channel. Similarly, the βG529C mutant was also slightly permeant to larger ions with an IK/INa ratio of 0.06 ± 0.0 (Table ) and an IRb/INa ratio of 0.03 ± 0.01
. The βG529D mutant, which exhibited a IK/INa ratio of 0.19 had a very low INa expression level and was not analyzed further. None of the
G587 or
S541 mutants showed significant K+ permeability (Table ). We conclude from the analysis of these macroscopic current data that (a) the residues
G587 and βG529 in the pre-M2 segment play a role in channel block by amiloride, (b)
G587, βG529, and
S541 are involved in defining channel Li+/Na+ permeability ratio, and (c) βG529 codetermines the molecular cutoff of the channel that normally prevents K+ ions from going through the channel.
|
G587, βG529, and
S541 Mutants
G587A and βG529A showed reduced unitary currents with Li+ ions as charge carrier, and no obvious changes in unitary Na+ current (iNa) or channel gating were observed. The corresponding
S541A mutation decreased iNa, and the unitary Li+ current (iLi) was considerably reduced so that transitions between open and closed states could not be precisely determined under our recording conditions (Fig. 4). Single-channel conductance values, which are summarized in Table , show that the unitary Li+ conductance (gLi) for the mutants follows the order wt >
G587A > βG529A >
S541A. The gLi/gNa ratios agree with the Li+/Na+ permeability ratios measured for macroscopic amiloride-sensitive currents (Table ). The higher Li+/Na+ permeability of the βG529R and the
S541G mutants compared with wt measured by the macroscopic ILi/INa ratio
, could not be confirmed by changes in either unitary currents for Na+ or Li+ ions, indicating that βG529R and
S541G mutations do not increase channel permeability to Li+ relative to Na+ ions.
|
|
G587 all known genes have a Gly residue, except for S541 in the ENaC
subunit. Substitution of the
S541 by Gly did not change the conductance or the ion selectivity of the channel. We also tested whether Gly and Ser residues are interchangeable at the homologous positions in
and β subunits. The
G587S mutant showed a decreased unitary Na+ and Li+ conductance, but the gLi/gNa ratio remained unchanged (Table ). The most dramatic effect was observed with βG529S. As shown before (Fig. 3), this mutation makes the channel permeable to K+ and preserves a macroscopic ILi/INa ratio >1. At the single-channel level, βG529S decreased unitary Na+ and Li+ current amplitudes to such an extent that transitions between open and closed states could not be detected. As shown in Fig. 5 in a representative experiment, channel openings or closures could not be detected under our conditions in excised outside-out patches containing channels with the βG529S mutation and exhibiting amiloride-sensitive Na+ and Li+ currents. This indicates that gLi and gNa of the βG529S mutant are more than likely <1 pS. Thus the analysis of the single-channel currents of the
G587, βG529, and
S541 mutants indicate that substitution of these residues by Ala, Ser, or Cys change the ion permeation properties by decreasing channel conductances for either Li+ ions alone or for both Li+ and Na+.
|
G587, βG529, and
S541 Residues
|
|
G587A mutant, which is equally permeant to Na+ and Li+, the apparent affinity for Na+ is slightly lower and the apparent affinity for Li+ ions is increased (Fig. 7). The KM of 114 mM for Na+ ions is close to that for Li+ ions (90 mM) (Table ). The
G587S mutant is characterized by a dramatic decrease in maximal conductances for Na+ and Li+ ions and a lower apparent KM for Li+. The unitary Na+ conductance of
G587S saturates at concentrations <100 mM. Due to the low gNa, the KM for Na+ permeation was not determined more precisely. The βG529A mutant shows a gLi that saturates between 2 and 3 pS with an apparent affinity for Li+ ions around 40 mM, thus the apparent affinity is increased threefold compared with wt. It is remarkable that the change in apparent affinity and maximal conductance of the βG529A mutant affects only Li+ ions, whereas the KM for Na+, and maximal gNa, remained basically unchanged (Table and Table ). Interestingly the βG529R mutation has no apparent effects on the permeation properties of Na+ and Li+ ions (Fig. 7). The barely detectable single-channel Li+ currents of the
S541A mutant (Fig. 4) made the determination of gLi impossible. The
S541A mutant showed a slight decrease in gNa, with a slightly lower apparent affinity for Na+ ions compared with ENaC wt. Thus mutations of
G587 and βG529 change the apparent affinity for Li+ and/or Na+ ions, resulting in alterations in the single channel conductance and the microscopic conductance ratio gLi/gNa.
What is the affinity of the
S541A mutant for Li+ ions? In a few excised outside-out macropatches, we measured macroscopic Li+ conductance of the
S541A mutant at five different Li+ concentrations, and obtained an apparent KM for Li+ as low as 11 ± 1 mM
. From this high apparent channel affinity for Li+ ions and low Li+ permeability of the
S541A mutant, we expected that Li+ ions should block Na+ currents. The experiment illustrated in Fig. 8 shows that this is indeed the case. Macroscopic, amiloride-sensitive inward currents were measured in 20 mM external Na+ and increasing Li+ concentrations from excised outside-out macropatches from oocytes expressing the
S541A mutant. 80 mM Li+ blocked 60% of the Na+ current in the presence of 20 mM external Na+, a Na+ concentration below the KM for Na+ of the
S541A mutant. The inward amiloride-sensitive current is not completely blocked by Li+ because Li+ still carries current in the
S541A channel. From the fit of the Li+ inhibition curve of the INa to a simple inhibition scheme that takes into account a nonblockable fraction (see METHODS), we obtained an apparent Li+ affinity for INa inhibition of 28 mM. In contrast to the situation in the
S541A mutant addition of 80 mM Li+ to the extracellular solution increased in the ENaC wt the amiloride-sensitive current by 3.1 ± 0.1-fold
at –150 mV. We conclude that the low unitary conductance of the
S541A mutant is related to a higher apparent channel affinity for Li+ ions.
Model of Na+ and Li+ Permeation in wt and Mutant ENaC
To interpret the changes in the maximal open state conductance and apparent affinity for Na+ and Li+ ions induced by
G587, βG529, or
S541 mutations, we have used a simple model of ion permeation through ENaC, based on reaction rate theory of ion permeation (Alvarez et al. 1992
; Hille 1992
; French et al. 1994
). According to this theory, the movement of ions through the single filing region of the pore is described as a series of discrete steps between energy minima (wells), separated by energy maxima (barriers). Our model incorporates a 3B2S kinetic scheme, outlined in Fig. 9 and described in METHODS. Single-channel current amplitudes, determined at various negative voltages and different concentrations of the permeant ion (either Na+ or Li+), were used to construct energy barrier profiles. The best fit of our I/V data and the conductance–ion concentration relationship to a 3B2S model provided energy profiles for wt and mutants shown in Fig. 9 B and listed in Table . The conductance–ion concentration relationship predicted by the model for the different channel types in the presence of Na+ or Li+ are shown as solid lines in Fig. 7. The extracellular part of the energy profile (G2, G3, U2, see Fig. 9 A) is much better defined by our data than the intracellular side (G1, U1) because we have obtained the energy barrier profiles by fitting exclusively inward currents. The best fit of our data predicts energy profiles for Na+ and Li+ permeation in wt ENaC that are not absolutely symmetrical, but still result in a near ohmic behavior of the I/V relationship over a voltage-range of ±150 mV. For the wt, the higher permeability of Li+ versus Na+ ions observed at concentrations >50 mM is mainly due to a higher energy of the outer well (U2) for Li+ permeation, corresponding to a lower affinity for Li+ binding.
The use of Eyring reaction rate theory has several limitations, among them the fact that it describes ion movements over well-to-well distances of fractions of Angstroms, whereas in our 3B2S model the ion translocation over the three energy barriers takes place over a distance of several Angstroms. In discrete-state permeation models, the calculation of the barrier energies depends on the use of a prefactor, kT/h according to the Eyring absolute reaction rate theory (Hille 1992
; see METHODS). In a different approach, an approximation to a continuum model, a prefactor including the ion's diffusion coefficient and information on the geometry of the barrier and well has been used to calculate barrier energies (Andersen 1989
; Andersen and Koeppe 1992
). The use of this prefactor yields barrier and well energies
4.6 RT units more negative than the values reported in Fig. 9 B and Table . Because of the uncertainty about absolute energies of barriers and wells, we illustrate in Fig. 9 C the changes in barrier and well energies (G2 and U2) obtained for the
G587A, βG529A, and
S541A mutants relative to wt. These relative changes are much less model dependent. The small changes in the Na+ permeation through these channel mutants are simply due to a parallel positive shift in G2 and U2 (Fig. 9 C), corresponding to a slight decrease in channel affinity for Na+ ions but a conserved maximal gNa. The gradual decrease in Li+ permeability and the changes in the gLi–[Li+] relationship observed in
G587A,
G587S, and βG529A mutants compared with wt are essentially due to a substantial increase in the peak energy of the middle barrier (G2) relative to wt and only small changes in the energy of the outer well (U2) (Fig. 9B and Fig. C, and Table ). This results in a higher energy (
G) required for Li+ ions to hop over G2 from an outer binding site U2 to the inner binding site U1, and a consequently slower translocation rate of Li+ ions from U2 to U1.
We have used the macroscopic current data of Na+ current block by external Li+ (Fig. 8) to obtain an energy profile for Li+ permeation through
S541A ENaC, as described in METHODS. The prediction of the Li+ block by the model is shown as a dotted line in Fig. 8, and the energy parameters for the fit are listed in Table and shown in Fig. 9 B. As for the
G587A and βG529A mutants, the corresponding energy profile predicted for Li+ permeation through
S541A shows a higher peak energy of the middle energy barrier that increases the energy required for translocation of Li+ ions from U2 to U1 (Fig. 9B and Fig. C). Thus, when Li+ occupies the pore of the
S541A in the single filing region, it slows the movement of Na+ ions.
Is the Pre-M2 Segment the Specific Target Region for Amiloride Block and Na+/Li+ Selectivity?
The question remains as to whether amiloride binds to a unique site in the external pore, and whether the ion selectivity filter is confined to the region we identified in the pre-M2 segment. Recently, Ismailov et al. 1997
identified a WYRFHY sequence rich in aromatic residues in the large extracellular loop of
ENaC as a putative amiloride-binding domain. Later, it was reported that conserved Lys residues at positions 550 and 561 in
ENaC are critical for channel sensitivity to amiloride and for ion selectivity (Fuller et al. 1997
). Since these studies were done on channels presumably made of
ENaC subunits reconstituted in lipid bilayers, we tested the relevance of these mutations for the pharmacological and functional properties of ENaC composed of the three types of subunits,
, β, and
.
Following expression of an
mutant with deletion of the WYRFHY sequence together with β and
subunits, no INa could be measured and no ENaC cell surface expression could be detected using anti–FLAG antibodies directed against β and
subunits. The point mutation within this sequence,
H282D, reported to be responsible for ENaC insensitivity to amiloride did not result in changes in amiloride sensitivity or macroscopic ionic selectivity of channels made of
H282Dβ
(Table and Table ). The
K550E and the
K561E mutants, when coexpressed with β and
subunits, exhibited the same sensitivity to block by amiloride as ENaC wt and showed normal macroscopic Na+/Li+ selectivity and unitary currents (Tables I–III). From these results, we conclude that the putative amiloride-binding sequence described by Ismailov et al. 1997
and the Lys residues in the large extracellular loop (Fuller et al. 1997
) are not relevant for the amiloride sensitivity and ionic selectivity of the native ENaC made of
, β, and
subunits.
| Discussion |
|---|
|
|
|---|
, β, and
ENaC subunits, located two residues upstream of the predicted NH2 terminus of the second transmembrane
helix, that play an important role in ion conduction. Mutations of these residues
G587, βG529, and
S541 resulted in the following changes in ion permeation and channel ionic selectivity. First, the mutations
G587D,
S541N,
S541F, and
S541R resulted in channels, which were expressed at the cell surface at normal densities, but did not conduct ionic current or showed only a minuscule ionic current (
S541R). Second,
G587A, βG529A, βG529C, and
S541A mutations decreased Li+ permeability with relatively small changes in Na+ permeability. Third, introduction of a Ser residue in the
(
G587S) and β subunits (βG529S) resulted in a decrease in both Na+ and Li+ unitary currents with no apparent changes in the Li+/Na+ permeability ratio. Fourth, the βG529S (and similarly βG529C and βG529D) mutant was permeant to K+ ions and to a lesser extend to Rb+ ions. Finally, the mutants
G587S, βG529A, βG529C, and βG529S with major changes in ion permeation or selectivity also showed resistance to amiloride block. These observations indicate that
G587, βG529, and
S541 interact closely with the permeant Na+ and Li+ ions, and participate at least in part in ion discrimination at the selectivity filter. The selectivity filter appears in close vicinity of the amiloride binding site.
The amiloride-sensitive ENaC, target of aldosterone action in the distal nephron and colon, is made of homologous
β
subunits. This channel is highly selective for Na+ over K+ ions and has the unique characteristic among the other members of the ENaC/DEG gene family of being more permeant to Li+ than to Na+ ions. The higher permeability of ENaC for Li+ compared with Na+ ions was evident from early macroscopic current measurements in flat epithelia like frog skin or toad urinary bladder. In these preparations, it was generally not possible to detect currents through ENaC carried by K+ ions. The only ion other than Na+ and Li+ allowed to pass through the channel was H+. The selectivity profile H+ > Li+ > Na+ of ENaC suggested that the pore discriminates among cations mainly on the basis of the size of the dehydrated ion allowing small cations to pass through the channel. Larger cations such as K+, Rb+, or Cs+ go only part-way along the ion conduction pathway, and block the channel pore at high concentrations, but are not able to pass the most constricted region of the pore, the selectivity filter (Palmer 1990
; Garty and Palmer 1997
).
In support of this hypothesis for ion selectivity in ENaC is our recent finding that mutation of the highly conserved
S589 residue alters the molecular sieving properties of ENaC (Kellenberger et al. 1999
). Specific
S589 mutants exhibit different IK/INa ratios. In particular, the
S589D mutant with the highest K+/Na+ permeability ratio showed a permeability profile for alkali metal cations following the order K+ > Rb+ > Cs+, indicating that
S589 mutations result in an increase in the molecular cutoff of ENaC. It was concluded that
S589 is part of the selectivity filter, which acts as a molecular sieve to discriminate between small Na+ or Li+ ions and larger K+ ions.
In the present study, the mutants that show a significant K+ permeability are βG529S, βG529C, and βG529D; the permeability relative to Na+ of these mutants was higher for K+ ions than for Rb+, and no current carried by Cs+ could be detected. The increase in the channel molecular cutoff by the βG529 mutations is consistent with the notion that, like the
S589 mutation, they enlarge the pore diameter at the selectivity filter and alter the molecular sieving properties of the channel (Kellenberger et al. 1999
). Thus the changes in the pore geometry at the selectivity filter that allow K+ or Rb+ ions to pass through the pore can be achieved by specific mutations of
S589 or βG529 residues. It suggests that these two residues are important for the steric selectivity of the channel by maintaining the proper pore geometry to tightly accommodate the permeating ion and exclude larger ions.
Beside the changes in ion selectivity of the βG529S, βG529C, and βG529D mutants, mutations of the
G587, βG529, and
S541 residues resulted in important changes in the permeability properties of the channel for Na+ and Li+ ions. These effects involve both the changes in the maximal open state conductance (Fig. 7) and/or in the apparent channel affinity for Na+ and Li+ ions as determined by KM values for channel conductance. Recordings of individual channels in the native tissue have shown that the open state conductance of ENaC saturates with increasing concentrations of Na+ and Li+ ions, indicating that Na+ and Li+ ions bind to specific sites in the ion permeation pathway. This saturation process arises when the binding–unbinding steps of the ion conduction at specific binding sites along the pore become rate limiting. At high ion concentrations, the channel pore is occupied most of the time and the rate of ion translocation is determined by the maximal rate at which the permeant ions dissociate from their binding site and the ion flux approaches saturation. For ENaC wt, the dependence of open state conductance on the concentration of the conducting ion gives an apparent affinity (KM) of 38 mM for Na+ and 118 mM for Li+. This difference in channel affinity for Na+ and Li+ ions can account for the higher Li+ over Na+ permeability usually observed at an ion concentration around 100 mM since the flux of Na+ ions already reaches saturation at these concentrations (see Fig. 7). Specific binding sites for Na+ have been postulated on the basis of competitive interactions between Na+ and amiloride or K+ ions that inhibit Na+ current through the channel in a voltage-dependent manner. In addition, the lack of voltage dependence of channel occupancy by Na+ ions was consistent with the presence of multiple ion binding sites along the channel pore. Reasonable models for Na+ permeation have been proposed that are consistent with these experimental observations. These models involve at least two common binding sites for Na+, Li+, K+, and amiloride in the outer mouth of the channel (Palmer 1990
).
How can the decrease in Li+ permeability relative to Na+ of the
G587A, βG529A, or
S541A mutants be explained in molecular terms? High affinity binding of ions in the channel pore tends to slow ion permeation provided that the channel is occupied by a single ion at a time. A low Li+ permeability was observed in
G587A, βG529A, or
S541A mutants characterized by a lower maximal gLi and an apparent higher affinity for Li+ compared with wt. Alternatively, the mutations might result in steric changes within the narrow region of the pore, making it more difficult for Li+ ions to pass through. For the interpretation of the changes in ion permeation in ENaC mutants, we need to consider ion conduction through ENaC consisting basically of three fundamental processes: (a) diffusion of an ion up to the outer entrance of the channel pore, (b) dehydration and solvation of the ion by polar chemical groups lining the pore, and (c) translocation through the selectivity filter. Models of ion permeation in terms of energy profile and Eyring rate theory can formulate these processes (Hille 1992
). We have empirically chosen an energy profile for Na+ and Li+ translocation through ENaC that consists of two energy wells representing two binding sites and three energy barriers for ion translocation. It has previously been shown that this type of model (3B2S) for ion conduction can account for the electrical properties of ENaC (Palmer and Andersen 1989
). In energy barrier models, the Kd of the permeant ion is defined by the value of the deepest energy well; thus, the deeper the well, the lower the concentration needed to reach saturation of the ion flux. By contrast, the higher the energy barriers, the lower the maximal conductance of the channel. While such models cannot provide an accurate representation of the structure of the ion conduction pathway, they can help us to determine whether alterations in ion conduction of
G587A, βG529A, or
S541A mutants can be related to particular permeation mechanisms. In the ENaC wt, the energy barrier profile obtained for Li+ permeation compared with Na+ is characterized by a higher outer energy well consistent with a lower affinity for Li+ at an external binding site (Fig. 9 and Table ), accounting for the differences in Li+ versus Na+ permeability. In our models for Li+ permeation through the
G587A, βG529A, or
S541A mutants, the predominant changes in Li+ versus Na+ conductance are due to an increase in the middle energy barrier, with little changes in the outer energy well. It suggests that translocation of Li+ ions in the single filing region at the central energy barrier requires more energy in mutants compared with wt. Thus, the
G587, βG529, or
S541 residues can be assigned to a region of the pore that represents a significant barrier for ion translocation. It is possible that βG529 and
S541 have preferential interactions with Li+ ions since mutations of these residues did not affect Na+ permeation much.
We observed that the βG529S, βG529C, and βG529D mutations change the channel molecular cutoff allowing K+ and Rb+ ions to pass through the channel. Thus, βG529 plays a role in defining the molecular cutoff of the channel and is at the same time important for the Li+/Na+ selectivity. The dual role of βG529 underlines that the discrimination between Li+ and Na+, and between Na+ and larger ions, occurs at overlying sites in the selectivity filter. Our data suggest that
G587A, βG529A, and
S541A mutations result in steric changes at the selectivity filter that impair translocation of Li+ ions through the narrowest part of the channel pore. In moving from the outer binding site U2 into the narrowest region of the pore, the permeant ion likely has to lose a few more molecules of its hydration shell. This loss of the hydration shell raises the energy of the permeant ion. Our data are consistent with the hypothesis that the energy cost for Li+ dehydration at the U2
G2 translocation is raised in the
G587A, βG529A, and
S541A mutants because it cannot be sufficiently compensated by the interactions between the permeant ion and the polar groups lining the pore. Mutations such as
G587S and βG529S decrease unitary conductance for both Na+ and Li+ ions, although the low unitary currents made it difficult for
G587S and impossible for βG529S to fit the I/V data to our 3B2S model. The energy profile of
G587S for Li+ conduction predicts an increased peak energy of the middle barrier, consistent with what we observed for the other mutants discussed above. The conductance data of βG529S are also compatible with a mutation resulting in a considerable increase in the peak energy of the middle barrier for Na+ and Li+ permeation, but a decrease in the energy barrier for K+ conduction. The binding site U2 cannot be assigned yet to a molecular structure, but is likely to be in the close vicinity of the selectivity filter.
It is somewhat surprising that the effects of substitutions of
G587, βG529, and
S541 did not correlate in a predictable way with general physical–chemical properties of the substituting amino acid residues. The fact that ion permeation is not altered by the βG529R mutation indicates that charges do not play a crucial role in ion coordination at the selectivity filter, and basically excludes the possibility that the βG529 side chain faces the lumen of the channel pore. As already suggested for
S589 mutations (Kellenberger et al. 1999
), the amino acid side chain of βG529 (and probably of its analogues) are likely to point away from the channel pore, in a similar way as has been shown for the selectivity filter of the KcsA channel (Doyle et al. 1998
). If this is the case, the side chains of these amino acid residues are likely to interact with other amino acid residues, to keep the backbone of
G587, βG529, and
S541 in place and maintaining the proper conformation of the pore at the selectivity filter. In ENaC wt, Gly, and Ser, amino acid residues with short side chains make these interactions, and the positive charge in the βG529R mutant (four carbon atoms between the
-carbon atom and NH3 group) might be too far away to disrupt this interaction. In mutant channels that are permeable to K+ (
S589 and some βG529 mutants), the mutated side chain might push parts of the pre-M2/M2 domains of the subunits away from each other or tilt them, resulting in a widening of the pore.
In summary, our mutagenesis experiments indicate that
G587, βG529,
S541, and
S589 residues are part of the narrow region of the pore that constitutes the selectivity filter of ENaC (Fig. 1 C). It seems that the important steps of ion conduction that involve ion dehydration, solvation, and translocation occur in a very restricted region of the ion permeation pathway that encompasses the
G587 and analogues and
S589 residues. For comparison, the crystal structure of the KcsA channel reveals that the narrow selectivity filter lined by the backbone of VGYG residues is only 12-Å long, whereas the overall length of the pore is 45 Å (Doyle et al. 1998
).
Finally, the mutations
G587S, βG529A, and βG529S that change ion selectivity and ion permeation properties also decrease channel affinity for amiloride, whereas the βG529 mutation without consequences on ion permeation properties did not change channel blocking by amiloride. It is unlikely that
G587 and βG529 are involved in specific binding interactions with amiloride since nonconserved substitutions have no effect on channel affinity for amiloride. The
G587A, βG529A, or
S541A mutations may indirectly impair the binding interaction of amiloride with the channel pore by steric changes that result in alterations in the pore geometry. According to its size, amiloride is not supposed to penetrate deep into the narrowest part of the selectivity filter. We propose that
G587 and βG529 are located at the external entrance of the narrow selectivity filter where steric alterations of the pore have consequences on amiloride binding in the channel outer vestibule.
Fig. 1 C, showing the ion pore of ENaC, illustrates some relevant experimental observations regarding the structures involved in ion conduction. The Gly 525 and 537 residues of β and
subunits and to a lesser extent the corresponding
S583 residue form the amiloride binding site in the pore (Schild et al. 1997
). Mutation of the two
S583 residues in the ENaC tetramer to Cys creates a high-affinity binding site for Zn2+ (Firsov et al. 1998
). Because the
subunits are most likely on opposite sides of the channel pore (Firsov et al. 1998
), and the optimal distance of two sulfhydryl groups for coordinated ligation of Zn2+ is
5 Å (Krovetz et al. 1997
), we estimate the pore diameter at this site to be around 5 Å. The pore then narrows down to its narrowest part that constitutes the selectivity filter. This narrow region is relatively short and may be lined by only three amino acid residues involving the residues Ser or Gly in the
β
subunits. We suggest that, in analogy to the KcsA channel, the cytoplasmic part of the pore is lined by residues from the second transmembrane
helix.
Portions of this work were previously published in abstract form (Kellenberger, S., N. Hoffmann-Pochon, E. Schneeberger, I. Gautschi, and L. Schild. 1998. J. Am. Soc. Nephrol. 9: A0193).
| ACKNOWLEDGMENTS |
|---|
This work was supported by grant 31-49654.96 from the Swiss National Foundation for Scientific Research (L. Schild).
Submitted: 5 March 1999
Revised: 23 April 1999
Accepted: 30 April 1999
| REFERENCES |
|---|
|
|
|---|
-bENaC regulate channel gating, ion selectivity, and sensitivity to amiloride, Biophys. J., 72, 1997, 1622–1632.[Medline]
, β,
subunits of the epithelial sodium channel (ENaC) involved in amiloride block and ion permeation, J. Gen. Physiol., 109, 1997, 15–26.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|