Abstract

Patch clamp studies of the potassium-transport proteins TRK1,2 in Saccharomyces cerevisiae have revealed large chloride efflux currents: at clamp voltages negative to -100 mV, and intracellular chloride concentrations >10 mM (J. Membr. Biol. 198:177, 2004). Stationary-state current-voltage analysis led to an in-series two-barrier model for chloride activation: the lower barrier (α) being 10–13 kcal/mol located ~30% into the membrane from the cytoplasmic surface; and the higher one (β) being 12–16 kcal/mol located at the outer surface. Measurements carried out with lyotrophic anions and osmoprotective solutes have now demonstrated the following new properties: (1) selectivity for highly permeant anions changes with extracellular pH; at pHo=5.5: I−≈Br−>Cl−>SCN−>NO3−, and at pHo 7.5: I−≈Br−>SCN−>NO3−>Cl−. (2) NO2− acts like “superchoride”, possibly enhancing the channel's intrinsic permeability to Cl−. (3) SCN− and NO3− block chloride permeability. (4) The order of selectivity for several slightly permeant anions (at pHo=5.5 only) is formate > gluconate > acetate phosphate−1. (5) All anion conductances are modulated (choked) by osmoprotective solutes. (6) The data and descriptive two-barrier model evoke a hypothetical structure (Biophys. J. 77:789, 1999) consisting of an intramembrane homotetramer of fungal TRK molecules, arrayed radially around a central cluster of four single helices (TM7) from each monomer. (7) That tetrameric cluster would resemble the hydrophobic core of (pentameric) ligand-gated ion channels, and would suggest voltage-modulated hydrophobic gating to underlie anion permeation.

Introduction

Uncoupled, channel-like, ion fluxes through certain ion-coupled substrate transporters began to be recognized about 20 years ago [84], and subsequently were widely observed in amine transporters and excitatory amino acid transporters of vertebrate neural tissues [27, 54, 70], as well as in a variety of invertebrates (reviewed in Ref. [76]). Some such movements were shown to be triggered by application of the principal substrate, but found to occur in non-stoichiometric quantities (transport-associated currents); others arose under special conditions in the absence of substrate (constitutive leakage currents). In different transporters, these uncoupled fluxes comprise mainly chloride ions or sodium ions, which augment the normal transport fluxes.

A discovery of patently outsized currents in the yeast Saccharomyces cerevisiae [9], associated with the potassium-uptake proteins Trk1p and Trk2p, led to a demonstration of uncoupled chloride fluxes through these transporters, as well [50]. That finding proved particularly interesting because TRK proteins are sequence similar to microbial potassium channels [24, 44], and have been shown to fold like K+ channels [23, 33, 48, 59, 84, 91]. That is, each single molecule of yeast Trk1p or Trk2p is folded into an internal tetramer with four pairs of transmembrane (TM) helices bracketing four recurrent strands to form a selectivity filter for K+. In Fig. 1a, these K+ filters are depicted as magenta clusters. The most distal pair of TM helices in each molecule, designated M1D and M2D (colored bars in Fig. 1a), is far better conserved among fungal species than are the other three pairs [61], and was postulated to coordinate a supertetramer of TRK monomers [23]. The axial pore of this structure would be lined by four copies of M1D, and could provide a pathway for chloride currents [68].

Models for TRK structure and for chloride conduction. a Tetramer co-organized by proximation of M1D helices in four separate TRK monomers: cylinder and wire diagrams in green, gray, blue, and red (counterclockwise) depict the backbone structures as originally...

The steady-state voltage-dependence of these outsized TRK currents resembles that of the mammalian glutamate transporter EAAT5 [1, 88] expressed in Xenopus oocytes. However, the implied Cl− flux through the yeast TRK proteins swamps those proteins' normal K+ transport function and is certainly not initiated by K+ transport [50]. The findings on Saccharomyces have proven consistent with the observation that KCl in conventional penetrating microelectrodes (1−3 M) strongly depolarizes the membranes of yeast-like spherocytes of Neurospora, a normally large-celled mycelial organism [10]. Given expected yeast membrane voltages near −200 mV, the underlying chloride permeability via TRK proteins probably serves to keep cytosolic Cl− below ~1 mM, even with extracellular concentrations >1 M. However, the actual functional reason for conspicuous cross-species sequence conservation in the M1D and M2D helices remains to be determined (see “Discussion” section).

In principle, chloride could be transported by mechanisms other than conventional ion channels or carriers since it behaves both as a weak chaotropic ion—capable of disrupting protein structure [17, 86]—and as a lipid-soluble ion [64, 85] which might diffuse through the hydrophobic regions in membrane proteins. Indeed, chaotropic anions are known to inhibit P-type transport ATPases [34, 49, 67, 83], to bias the Na+, K+ ATPase toward its E1P conformation [49, 80], and to modulate nerve and muscle behavior by means of altered membrane surface potential [20, 37, 40]. At least part of their lipophilicity is thought to lie in an ability to combine with protons and to diffuse through phospholipids as neutral molecules (viz., HCl, HNO3, HSCN, etc. [31, 64, 85]). Such action could dissipate transmembrane pH gradients, uncouple proton-coupled ion transport reactions, and induce pseudo-anion currents associated with membrane depolarization.

In the present case, however, such essentially physical-chemical transfer processes can be largely discounted [50], because (a) deletion of both of the yeast TRK proteins eliminates more than 90% of the chloride-dependent currents; (b) the current–voltage relationship in the presence of chloride is strongly rectifying, with exponential (Eyring) kinetics, rather than Nernst–Planck kinetics; and (c) elevated intracellular chloride evokes inward currents— implying net chloride efflux; whereas elevated extracellular chloride evokes no corresponding outward currents.

The present experiments were undertaken to determine the selectivity of the TRK-Cl− pathway, particularly with regard to other halide and chaotropic ions; and to determine the extent of validity of the two-barrier kinetic model previously developed for Cl- currents. The resultant data demonstrate that other halides permeate even better than chloride itself; that more chaotropic ions can supplant/block halide ions; and that osmoprotective solutes can choke the currents due either to halide ions (Cl−) or to strongly chaotropic ions (SCN−). All of these findings are compatible with anion conduction through a central (fifth) pore in the K+ transport proteins, which would be formed by spontaneous membrane assembly of TRK tetramers [23], and which bears a strong structural resemblance to the well-known ligand-gated anion channels.

Methods

Strains

Four strains of S. cerevisiae were used in these experiments: our principle “wild-type” strain EBC202; its parents BS202/203; the derived double knockout strain for genes TRK1 and TRK2: EBK202; and a separate “wild-type,” K837. The latter strain was provided by Dr. Kyle Cunningham (Department of Biology, Johns Hopkins University, Baltimore); and BS202/203, by Dr. Albert Smith (Yale Department of Molecular, Cellular, and Developmental Biology; now of the Icelandic University, Reykjavik). The two EB strains were constructed in this laboratory by Dr. Esther Bashi using standard yeast genetics techniques [47]. All of these strains are equipped with an array of useful markers, and their full genotypes are given at the bottom of this paragraph. The TRK1 and TRK2 genes, deleted from EBK202, encode the two proteins which mediate active transport of potassium in Saccharomyces, and—under special conditions—channel the conspicuous anion effluxes described in detail below [50, 68]. The PMR1 gene, deleted from K837 and from the EB strains, encodes a Golgi-resident calcium pump, whose deletion results in useful enlargement of both intact cells and spheroplasts of yeast, so that cells ~8 μm in diameter can routinely be selected for patch clamping. Detailed genotypes are as follows.

Growth media, spheroplasting

Cells were normally grown in shaking liquid YPD-Ca medium (1% yeast extract, 2% bactopeptone, 2% glucose, plus 10 mM CaCl2), at 30°C, to a final OD600 of 0.8–1.2; the added calcium was required for normal growth of pmr 1 Δ strains. Also, for the trk1Δ and trk1Δ −trk2Δ strains, 90 mM KCl was added to YPD-Ca. In experiments with molar concentrations of intracellular osmoprotective solutes (Figs. 5 and ​and8),8), 1 M glycerol was added to growth medium and was also added to the extracellular media for protoplasting and recording. Glycerol, at concentrations up to 1 M, has no readily observable deleterious effect on yeast, other than slight retardation of growth (doubling times of 1.9–2.1 h, instead of the normal 1.75 h, at 30°C), and has been found useful as a chemical chaperone [26, 91].

Osmoprotective solutes choke currents of chaotropic ions. I–V plots with buffer G-SCN supplemented by molar concentrations of glycerol or glycine betaine (pHo 5.5 only). All three plots fit by Eq. 1, with all parameters found in common except...

Patch recording

Borosilicate pipettes, prepared as previously described [7], were connected to an EPC9 amplifier [68], which was managed by “Pulse” software (HEKA Elektronik, Lambrecht, Germany) on a PowerMacG4 computer (Apple Computer Inc., Cupertino, CA, USA). Yeast membrane currents were measured in the “whole-cell” recording mode, using a staircase of 1.5 or 2.5-s voltage pulses (from a holding voltage of −40 mV) to +100 mV, +80…−180, −200 mV, in 20 mV decrements. Steady-state current-voltage curves (I–V curves) were generated by averaging the values sampled between 50% and 90% of each trace and plotting the averages against the clamped voltage for each pulse. Data were usually collected between +100 and −180 mV, but the actual data displays were truncated at +20 mV, since positive (outward) currents were unrelated to the TRK proteins [7, 50]. Resulting I–V curves were corrected for current leakage by subtracting an ohmic line corresponding to the minimal slope near the origin [50]. Small residual scatter was observed in the apparent short-circuit currents due to the very low slope conductance near-zero voltage, but this could be minimized by adding an error term (Ie) to the descriptive transport model fit to each I–V plot (see Eq. 1 below). The common sealing adjuvant, potassium fluoride, was tested (at mM concentrations) in pipette solutions for a few experiments, but had no effect on leakage currents; higher concentrations were avoided because of likely interference with chloride permeability. Data were routinely collected at 2 kHz and filtered at 250 Hz.

The model

Previous quantitative analysis of TRK-chloride currents showed that most of the steady-state data can be described equally well by either of two models. The first and more familiar of the two is a conventional voltage-gated channel whose characteristic gating voltage is strongly pH-sensitive: shifting ~37 mV positive per pH unit, from −270 mVat pH 7.5 to −160 mVat pH 4.5.

The other model, illustrated in Fig. 1b, postulates that anions transiting the TRK proteins must cross two activation–energy barriers, twice acquiring excess energy in a manner similar to ordinary chemical reactions [12, 38, 46, 52, 53]. This two-barrier model has proven the more cogent, for two main reasons: first, occasional yeast spheroplasts which tolerate very large clamping voltages (−220 to −300 mV) all display exponentially increasing conductances at the extreme voltages (Fig. 13 in Ref. 68). This feature is predicted by the two-barrier model, whereas a limiting (ohmic) conductance is predicted by the gated-channel model.

Second, transit through two barriers accords qualitatively with the association of four TRK molecules, described in Fig. 1a, as originally proposed by Durell and Guy [23]. The resulting central pore would be lined by four copies of the sequence 1075-LRKQLSFDLWFLFLGLFIICIC-1096 (M1D, in Saccharomyces Trk1p). The side chains of Trp1084s and Phe1091s would point directly into the channel and form tight rings, viz., energy barriers against the passage of chloride ions. Furthermore, Arg1076 and Lys1077 would lie in a wide vestibule at the inner end of the central pore, where their eight positive charges should act as an electrostatic sieve for anions.

A simple mathematical formulation of the two-barrier model expresses voltage-driven transit of single chloride ions, under stationary-state conditions, as Equation 1:

where, Im is the measured current, and Ie represents scattered error in the short-circuit current for each I–V curve. Cli and Clo are the intracellular (pipette) and extracellular concentrations of chloride ions, with their valence zc=−1; zcF converts chemical flux to current, and FVm/RT is the “reduced” (dimensionless) membrane voltage. a1 represents the fractional distance of the first barrier (α) from the inner face of the membrane dielectric, b2 represents the fractional distance of the second barrier (β) from the outer face of the membrane dielectric, A=exp(Eα/RT) is the reciprocal of the Boltzmann expression relating reaction velocity to activation energy (Eα) at barrier α, and B=exp(Eβ/RT) is the corresponding term for barrier β. The membrane's electric field is assumed to be constant, under steady-state conditions, so that the voltage offsets at barriers α and β are specified by the products a1Vm and b2Vm, respectively.

These equations contain four intrinsic parameters determining the detailed shapes of current–voltage curves under different experimental conditions: a1, b2, A, and B. However, only three of these actually vary; in trial optimizations, b2 usually assumes its minimal value (zero) and has therefore been fixed at zero, for all of the analysis described below (see Fig. 5 in Rivetta et al. [68] and Fig. 1b, here). Finally, ϕ is an experimental scaling factor used mainly to rationalize data from spheroplasts of different size. It cannot be determined independently of A and B, and must be fixed (usually=1.0) for at least one plot in each I–V dataset. Additional details are given at the end of this section, and in the legend to Fig. 1.

[Caveat: a perhaps unconventional feature of this model is the required variation a1 of barrier α, upon large changes in subtrate ion concentration or upon the addition of neutral osmoprotective solutes (see, e.g., Figs. 2 and ​and6).6). That requirement could be suppressed by introducing a third activation barrier, with its own position and amplitude: viz., two more parameters. But since the environment of biological membranes: including pH, ionic strength, electric field, boundary solutes, etc., can affect detailed membrane structure, the usual preference for modeling with fixed-position barriers seems arbitrary. For the present measurements and description, it is also unnecessary.]

Osmoprotective solutes choke halide currents through TRK. a I–V plots with buffer G-Cl supplemented by molar concentrations of five different neutral solutes (pHo 5.5 only). All six plots fit by Eq. 1, with all parameters found in common, except...

Equation 1 extrapolates along an exponential asymptote at large negative membrane voltages [68], and— because b2=0—along a horizontal asymptote at large positive voltages (Equation 2):

limVm→−∞Im=ϕ[zcFCliA]exp(zcFa1Vm/RT),limVm→+∞Im=ϕ[−zcFCloB]⋅

(2a, b)

In order to reconcile units on the right-hand side of Eqs. 1 and 2a, b with the units of current, it is convenient to assign units of permeability (cm/s) to the scaling factor ϕ. Then, because the data: measured currents, membrane voltage, Cl− concentration, and cell size are in pA/cell, mV, mM, and (μm2 surface, respectively, the practical (fitted) values of A and B must be rescaled by the factor 1.93×105. This is numerically equal to 2F, where 2 represents the standard yeast spheroplast surface area of 200 (μm2. Actual curve fitting has been carried out via a variant of the Marquardt algorithm [58], augmented with a parameter-coincidence matrix to accommodate multiple simultaneous plots.

Results

Confirming the chloride conductance

Previous experiments [68] showed that chloride currents through the yeast TRK proteins could be quantitatively described via the two-barrier model, allowing only a single parameter to vary, either with changes of [Cl-]i or with extracellular pH (pHo). The present repeat of those experiments has reinforced the general observation, and has demonstrated numerical stability of the actual parameters, as illustrated by Fig. 2. The current traces in Fig. 2a were obtained from two different spheroplasts: one patch perfused with 183 mM Cl−(buffer G-Cl), and the other with [Cl-]i lowered to 30 mM by potassium gluconate replacement of KCl. Two values of pHo were tested on the first spheroplast: 7.5 (buffer D) and 5.5 (buffer 5.5). The “inward” or negative currents, plotted downward, were clearly diminished by lowering [Cl-]i and by raising pHo. The latter effect mimicked an 80 mV negative offset of the I–V plot, with no obvious change of shape; this observation essentially coincides with an offset of 37 mV per pH unit reported previously [68].

Fits of Eq. 1 to the resulting I–V plots (Fig. 2b) were carried out with a1 set identically to the previously determined values at the two Cl− concentrations, and allowing only parameter B (height of barrier β) to vary between the two pHo values (documented in Fig. 2 legend). The purported barrier shapes and positions are diagramed in inset of Fig. 2b, to show that lowering pHo from 7.5 to 5.5 lowered the (β barrier by 1.9 Kcal/mol, and that reducing [Cl-]i from 183 to 30 mM moved the α peak deeper into the membrane dielectric (0.279→0.393). As noted by a caveat in “Methods” section, displacement of the α peak by changes of [Cl-]i could be suppressed by adding a third activation barrier (having two more parameters), but with little effect on the quality of fit. Thus, the two-barrier model describes these TRK-mediated chloride currents both efficiently and consistently with previous findings.

The same two intracellular ion concentrations, and extracellular pH values, were used as standard test points for further experiments.

Bromide, iodide

Currents of Br− and I− were assessed by replacing KCl in buffer G-Cl with either 175 mM KBr or 175 mM KI (buffers G-Br and G-I). Resulting patch clamp traces (Fig. 3) showed significantly larger currents than with chloride, at both pHo values tested and for both halides. Furthermore, the relative spacings of traces, per 20-mV step, differed for Br− and I− from that for Cl−; and quantitative analysis of this feature identified enlarged leakage currents relative to those for Cl−. Leakage–current amplitudes and conductances are indicated in Fig. 3 by the gray triangles, which also point to the null currents in each stack of traces. Finally, traces at the most negative voltages (−160, −180 mV) bore occasional bursts of noise with the character of conventional gated ion channels. Control measurements on yeast spheroplasts deleted of the TRK genes (trk1Δ−trk2Δ;Fig. 3b) also revealed enlarged non-TRK-related currents in bromide or iodide, relative to chloride.

Figure 4 demonstrates the leakage-corrected I–V plots for Br- (red diamonds) and I− (green squares) to be nearly identical with each other at both pHo values, but to be distinct from those for Cl− (black disks). Equation 1 could easily be fitted to these six datasets, when the intracellular mixture of 8 mM residual chloride plus 175 mM bromide or iodide was treated as if it represents a single ionic species at 183 mM (necessitated by the fact that Eq. 1 does not explicitly include terms for different ion species). All significant differences among the I–V plots could be generated by small shifts of activation energy at barrier β (outer face of the membrane; Fig. 4 legend). Clearly, there is no qualitative distinction among the voltage dependences of chloride, bromide, or iodide currents through the yeast TRK proteins. This peculiar anionic pathway, which was previously shown to be non-saturating for Cl− [68], is therefore a generalized halide channel, being about equally permeable to I− and Br−, and somewhat less permeable to Cl− in the range of experimentally accessible voltages.

The TRK inwardly rectifying pathway is permeable to all three common halides. Plots represent averaged I–V data from all halide records generated as in Fig. 3, for wild-type spheroplasts. Smooth curves calculated with a1=0.279 as previously determined...

Nitrite ions augment the chloride currents

Apart from very strongly chaotropic ions (see below), most other anions tested carried very much smaller currents through the TRK proteins than did chloride—typically less than 10%—when fully replacing KCl in the intracellular solution; also, they failed to block chloride currents when added at lower concentrations (~40 mM [50]). Tested ions included polyaspartate, divalent anions such as sulfate, selenate, and glutarate, and a variety of monovalents: azide, acetate, chloroacetate, isethionate, bromate, bicarbonate, and dihy-drogenphosphate. Formate and MES (zwitterionic) were minor exceptions—supplanting 15–20% of the chloride currents. However, nitrite proved a major exception: 30 mM nitrite, replacing equivalent KCl in the patch– pipette solution (buffer G-NO2), dramatically augmented the total TRK-mediated currents, as shown in Fig. 5.

How the recorded currents actually partition between chloride and nitrite was not determined, but measured current amplitudes roughly correspond to expectation for pure chloride concentrations of 500 mM at pH 7.5 and 250 mM at pH 5.5 [68]. Detailed analysis of the I–V plots via Eq. 1 yielded parameters fully consistent with the halide measurements in Fig. 4: in particular, the same position and value of barrier α sufficed for nitrite as for all the halides (a1=0.279, Eα=11.2 Kcal/mol). Only lowering of barrier β by less than 1 kcal/mol was required by the 30 mM NO2−, an effect which could most simply be related to increased permeability of the 153 mM Cl− still present.

Osmoprotectants modulate halide currents

The suggestion that nitrite ions might enhance chloride permeation through TRK proteins prompted systematic testing of reagents expected to have an opposite effect via membrane stabilization. Most familiar of these so-called kosmotropic agents are trehalose, glycine betaine, and glycerol [13, 19, 26, 60], to which we added also glucose and sorbitol, as controls for pure osmotic effects. Preliminary experiments, in which KCl in buffer G-Cl was partially replaced by sorbitol, had yielded unexpectedly small TRK-dependent currents [50], and molar glycerol had been shown to serve as a chemical chaperone [26, 91] without hampering growth. Data from a full experiment with molar solutes are shown in Fig. 6a.

Resulting I–V plots all had sensibly the same shape, without any offset along the voltage axis, and all could be described by the same parameter set allowing only the scaling factor ϕ to vary from curve-to-curve. Perhaps as a consequence of growth conditioning in glycerol, the α peak appeared slightly higher and deeper in the membrane (Eα= 11.2→11.5 Kcal/mol; a1=0.279→0.329) than for pure chloride experiments (cf. Fig. 5 legend with Fig. 3 legend); but Eβ was in the normal range. For the whole family of I–V plots, conductances declined exponentially toward zero voltage, and the conductance ratios (Gm/Gcontrol) were flat lines at the fitted values of ϕ, as shown in inset of Fig. 6a.

Apart from impossible changes of cell size, at least two distinct physical processes could simply scale down the TRK-dependent Cl− currents: (1) different osmoprotectants blocking different fractions of the TRK-dependent chloride pathways, each in all-or-nothing fashion; or (2) osmoprotectants reacting with all TRK molecules in the yeast membrane, in a manner to elevate Eα and Eβ by the same fractional amount (different fraction for each compound). This interpretation is represented in Fig. 6b and fits naturally into the general concept of the two-barrier model.

Inhibitory effects of strong chaotropic ions

Very chaotropic anions seem to destabilize Saccharomyces spheroplasts, but we were able to complete about two dozen current–voltage scans with 30–175 mM thiocyanate or nitrate. Current records for 30 mM thiocyanate displayed the usual spread due to large leakage currents (Fig. 7a, inset). After correction for this effect, SCN− actually reduced the total TRK-mediated current at pHo 5.5 and did so in manner to suggest that chloride conductance was blocked —despite the predominance of 153 mM Cl− in the pipette solution. Furthermore, identical inward currents were observed with 30 mM SCN−, whether the major anion was Cl− itself (diamonds) or was the barely conducted reference anion, gluconate (open circles at 145 mM, +8.3 mM Cl− from other salts in buffer G-Cl). Finally, transport appeared chemically saturated at 30 mM, showing no significant increase of plotted currents at 175 mM (triangles). [For reference, voltage-driven currents with chloride only (Fig. 7, control) not only did not saturate as a function of concentration, even out to 590 mM [68], but depended more steeply on membrane voltage—prominently so at pH 5.5.] These comparisons became more complicated at pH 7.5 (Fig. 7b), but again with no distinction between chloride itself and gluconate. And for 175 mM SCN−, currents were actually reduced, relatively to those for 30 mM SCN−.

I–V plots in the presence of 30 mM nitrate ions differed quantitatively from those with 30 mM SCN−, in being superposable for pHo 5.5 and 7.5 and lying roughly halfway between the two control plots for Cl− ions (data not shown). Although the general meaning of this result is unclear, decreased current relative to Cl− alone, at pH 5.5, implies that NO3− ions—like SCN− ions—also block chloride transport.

Equation 1 has not been fitted to the I–V plots for nitrate and thiocyanate because of difficulty in specifying concentrations, with more than one permeant species. It is nevertheless clear that these very chaotropic ions react quite differently with the TRK protein than halide ions do; most importantly, SCN− and NO3− ions block chloride currents. Easily demonstrated saturation in the latter effect implies a non-competitive mode of inhibition.

Osmoprotectants also choke thiocyanate currents

The results above imply that very chaotropic non-halide ions restructure the anion pathway through TRK proteins, likely by binding tightly within the channel. The question naturally arises, therefore, of whether or how osmoprotective solutes might affect conduction of the chaotropic ions. Tests were carried out with 1 M glycerol or 1 M glycine betaine included in the patch pipette, along with 30 mM thiocyanate. The mechanical stability of yeast spheroplasts was improved by both agents, and the SCN− -related currents (Fig. 8) were reduced, though not as effectively as Cl− currents had been (Fig. 6, above). Currents remaining in 1 M glycerol were 0.78×control, and in 1 M glycine betaine were 0.25 × control. The ratioed conductance lines in Fig. 8 inset show these two fractions clearly. Again, only small barrier elevations, Eα and Eβ, were required.

Computing ionic selectivity

Since the I–V relationship defined by Eq. 1 is strongly non-linear, conductance (dIm/dVm) and conventional permeability become convoluted measures of ionic selectivity. Therefore, practical comparison of the channel's selectivity for different ions requires another measure, one applicable over the testable range of membrane voltages. The I–V plots themselves—in Figs. 3, ​,4,4, ​,5,5, and ​and7—suggest7—suggest that chord conductances (Gm) at −180 mV, if corrected for ion concentrations and cell size, would reasonably mirror average ionic permeabilities between 0 and −180 mV. The appropriate calculations are summarized in Table 1.

With salts at high concentration, the highest permeabilities obtained for iodide, decreasing in the order I− ≈Br− >NO3− > SCN− >Cl− at pH 7.5; and in the order I− ≈Br− >Cl− >SCN− > NO3− at pHo 5.5. Decrease over the range is about twofold, and the increase of permeability from pH 7.5 to 5.5 is also approximately twofold. At the low concentration (30 mM), there are no data for bromide and iodide, but the other calculated permeabilities decline from nitrite, in the order NO2− ≈SCN− >NO3− at pH 7.5; and in the order NO2− > SCN− ≈NO3− >Cl− at pH 5.5. Again, decrease over the range is approximately twofold, but without a significant effect of pHo. These are not large differences in permeability, and bespeak a rather non-selective anion pathway through the TRK proteins. However, currents and calculated permeabilities were much smaller for a wide range of other anions tested at pH 5.5 [50] and ordered in Table 1; Cl−formate>gluconate>acetate phosphate−1.

The relatively large permeabilities calculated for 30 mM salts, compared with 175 mM salts (10–15×10−7 cm/s, cf. 2–5×10−7 cm/s), result from strong non-linearity between measured currents and concentrations, for both SCN− (Fig. 7) and NO3−.

Discussion

General features

The purpose of this investigation has been to complete a quantitative, stationary-state description of the surprising phenomenon of anion conduction through the TRK potassium transport proteins of S. cerevisiae [50, 68]. Stationary-state measurements have been enforced by the yeast, since all practicable conditions for whole-cell patch clamping have produced “always on” TRK-Cl− currents, or at least, on- and off-switching faster than the patch clamp amplifier.

Given the operational necessity of stationary-state analysis, the two-barrier model of Eq. 1 provides a simple, but quantitative description of the TRK-mediated anion currents, which can be shown to apply equally well to either protein, TRK1 or TRK2, and to their mixture in wild-type cells (Cl− experiments [50]). A remarkable additional fact is that the measured stationary-state I–V relationship for yeast TRK proteins strongly resembles the whole-cell I–V relationships for a wide variety of plant cells, viz., with conspicuous inward rectification at large negative membrane voltages [6, 18, 21, 30, 87].

Most surprising about the two-barrier formulation for TRK-mediated currents is the range of data which it can accommodate with almost all parameters held in common for different test conditions. Indeed, of the four parameters intrinsic to Eq. 1, the two positional parameters proved very stable: b2 being close to zero and completely insensitive to our experimental manipulations; and a1 responding only to large changes of permeant-ion concentration (intracellular; see esp. Fig. 2). The latter fact suggests that changes of ionic strength, in the cytoplasm (pipette solution), should contort the intracellular end of the channel.

Parameter B, representing the height of the outer barrier (β), has proven by far the most important parameter in the model. Modest changes in its value alone can accommodate all I–V differences related to changes of pHo, regardless of which permeant ions are involved (Figs. 2, ​,4,4, ​,5,5, ​,7);7); and as well, differences of selectivity between Cl− and Br− or I− (Fig. 4), or between chloride and its admixture with nitrite (Fig. 5). The observed practical significance of B owes both to its mathematical role as a determinant of shape in the midrange of I–V data (viz., −100 mV to −180 mV), and to the greater height of that barrier, relative to barrier α (Eβ>Eα by 2–4 Kcal/mol). The calculated sensitivity of B to changes of pHo is intuitively reasonable since barrier β would reside at the interface with extracellular solution; but adjustment of B in response to changes of the intracellular ionic species is more difficult to rationalize, and presumably must involve a cascade of structural adjustments, as different ions transit the channel.

Selectivity and structural homology

There are two dominant functional characteristics of TRK-mediated anion currents through the Saccharomyces plasma membrane, as demonstrated in Figs. 2, ​,4,4, ​,5,5, and ​and77 and Table 1. (1) Low selectivity among halide and non-halide chaotropic ions, but strong discrimination against all other anions which have been tested. (2) Very strong inward rectification: conduction is activated by large negative voltages and represents large anion effluxes, but is essentially blocked at positive membrane voltages. Another important characteristic is (3) insensitivity to conventional anion channel blockers such as disulfonic stilbenes [50].

The apparent order of selectivity, at the high salt concentration (175 mM) varies somewhat with pHo: I−≈ Br−>NO3−>SCN−>Cl−, at pH 7.5; and I−≈Br−>Cl−> SCN−>NO3−, at pHo 5.5; but the permeability change over the whole range is only approximately twofold, as is the change from pH 5.5 to 7.5. [Similar relationships hold at the low salt concentration (30 mM), even though the calculated permeabilities are several fold higher due to non-linearity of the I–V curves.] Thus, the TRK pathway for anions is essentially non-selective with respect to halides and chaotropic ions, even though it clearly discriminates against other small anions, viz., Cl−formate>gluconate>acetatephosphate−1. This sort of selectivity is unknown among ClC class channels and Cl− -conducting ABC proteins [16, 55, 81], but resembles that in glycine-receptors (GlyR) and related ligand-gated anion channels [12, 25, 56].

A very broad sequence search failed to identify any structures homologous to the fungal TRK-M1D helices, among bona fide chloride channels, transporters, or pumps. This negative information, plus the characteristic selectivity, noted above, of the TRK–anion pathway prompted us to review sequences and structures in ligand-gated ion channels (LGICs). The central pore of LGICs comprises a barrel stave array of five transmembrane helices (designated M2 or α2), one from each monomer within an assembled pentameric protein. And the M1D helix of yeast TRK proteins bears remarkable similarity to those LGIC barrel staves. Figure 9a compares two yeast M1D sequences (+flanks) with M2 sequences from two anion-selective LGICs (GlyR, GABAAR) and two cation-selective LGICs (nAchR, ELIC). Figure 9b compares the opposition of two pore-lining helices in the Durell and Guy model for SpTrk1 with two pore-lining helices from the recent crystal structure for ELIC (pentameric [11, 35]).

Charge selectivity (anion versus cation) of LGIC channels can be controlled by the residue at position −1′. When that residue is uncharged, like alanine in wild-type GlyR or GABAAR, the channel is anion selective; when it is negative, like glutamate in nAchR, ELIC, or certain mutants of GlyR, the channel is cation selective [63]. It seems likely that charge selectivity in the yeast structure would follow a similar rule. We are currently testing this hypothesis in ScTrk1 having glutamate substituted for arginine and lysine at positions −1′ and 0′, and also for the obstructing tryptophane and phenylalanine residues at positions 7′ and 14′ (mutants constructed by Dr. Jost Ludwig, University of South Bohemia, Czech Republic).

Both pores in Fig. 9b are depicted in presumed closed states, which would be enforced by hydrophobic exclusion of water, and therefore exclusion of hydrated ions. For LGICs, elaborate machinery has evolved in the surrounding protein to control gating; and there is much discussion of how that might work, focusing on Arg/Lys-0′ and adjacent residues, on the center of the membrane adjacent to Leu-9′, and on changes in tilt of the M2 helices [11, 36, 62, 65, 82]. For the yeast TRK-Cl− conductance, no organic ligands, TRK segments, or β subunits are presently known to effect gating. The only secure gating factor is strong membrane hyperpolarization, as modulated by extracellular pH: at pHo 4.5, the apparent gating voltage is near −160 mV; and at pHo 7.5, near −270 mV (Fig. 2 above; see also [68]). High-voltage gating, sometimes called punch through [29, 92], is qualitatively consistent with the predominance of bulky, hydrophobic side chains within the pore, as depicted in Fig. 9b (left panel). Indeed, that structure should be an excellent experimental model for hydrophobic gating [5, 3, 45, 90], wherein rapid conductance transitions and sharp thresholds are expected [4, 42].

Protein stabilization

Hydrophobic gating may also provide a plausible physical mechanism for osmoprotective agents to scale down Cl− and SCN− currents through the TRK proteins (Figs. 6, ​,8).8). Channels poised at the threshold for hydrophobic closure could be triggered by increasing the osmotic (diffusion) gradient between the channel interior and neutral but hydrophilic solutes in the surround. Thus, the number of channels conducting at any one voltage or time—not the intrinsic barrier amplitudes (see text for Fig. 6)—would be modulated by osmoprotective solutes.

Alternatively, suppression of TRK-chloride currents by organic osmolytes can be viewed as a specific case of such osmolytes' well-known ability to stabilize protein structures. In other kinds of experiments from the literature, polyalcohols and sugars—including the above-tested agents—have been shown to protect proteins from disruption by heat, salt, and freezing [2, 13, 43, 66, 71]. And some, especially glycerol and glycine betaine, have been widely adopted in evolution as “compatible” solutes, displacing cytoplasmic salt in salt-tolerant microorganisms and plants [51, 89], while trehalose has turned up acutely in most organisms which can survive complete dedhydration [41, 79].

In the specific case of the putative TRK oligomer, all four effective osmolytes could stabilize TRK either (1) in a monomeric state where the facultative chloride channel does not exist, or (2) in a water-excluded oligomer, viz., hydrophobically gated shut. Discrimination between these two possibilities should be revealed by elevating temperature. In the first case, the temperature coefficient of the effects in Figs. 6 and ​and88 should be unequivocally positive (suppression of currents with rising temperature), whereas in the second case an interval of negative temperature coefficient would be expected (enhancement of currents with rising temperature). Experiments into the influence of temperature on TRK-chloride currents are being undertaken. More generally, thermal denaturation may be a very useful way to characterize the underlying process of protein stabilization by kosmotropic agents.

[One trivial suggestion can easily be dismissed: these osmolyte effects upon TRK-Cl currents do not depend on altered ionic properties of the cytoplasmic (pipette) solutions. This was shown by a simple comparison of the observed current scaling (Figs. 6, ​,8)8) with effects of molar osmolytes on measured solution conductances and on actual salt concentrations, as affected by the partial volumes of osmolytes (A. Rivetta & C. L. Slayman, unpublished expts.). In particular, 1 M glucose and 1 M sorbitol both depressed buffer G-Cl conductance by ~35%; but glucose had no effect on the TRK-Cl currents, while sorbitol reduced those currents by ~63%. And, 1 M glucose and 1 M glycine betaine have the same concentrating effect on the KCl-to-water ratio (viz., the same partial molar volumes, ca. 100 ml/Mol), but only glycine betaine diminished the TRK-Cl currents, and that by ~89%.]

End notes

A major handicap to hard mechanistic analysis of this fungal TRK-anion channel is the lack of independent evidence for oligomerization within the membrane to form a pore that is not committed to K+ transport. While a tetramer is easiest to imagine, trimers and pentamers might also serve, or perhaps even a dimer. The latter has, in fact, been demonstrated in a crystal structure for the bacterial homologue VpTrkH (inVibrio parahemolyticus; 3.5Å resolution, M. Zhou, personal communication [14]). That structure confirms a K+ channel-like arrangement for the two potassium pathways [23] and shows the bacterial protein to be dimerized in the crystal at a complicated hydrophobic interface involving M1D and M2C. Thus, oligomerization of TRK-type proteins does occur. However, outside the actual K+ pathway, VpTrkH and the fungal TRK proteins share no sequence conservation; and even within the M1D helix itself, overlapping leucine (L/I) heptads found in the fungal proteins are not copied in the bacterial protein. Thus, identical oligomerization of TRK monomers in Saccharomyces compared with Vibrio seems unlikely.

Apart from the fact itself—of anion conduction through a K+-transport protein—two particular features of the stationary-state data will demand detailed mechanisms: the anomalous interactions of chloride with thiocyanate (Fig. 7), nitrate, or nitrite (Fig. 5); and the very strong inward rectification, which allows essentially zero outward current (Cl− influx) via the TRK-Cl− pathway. Anomalous (non-linear or paradoxical) permeabilities, frequently observed when two ionic species traverse the same transport pathway, have classically been attributed to interactions within multi-ion, single-file channels [39]. But other schemes also work, including charged wide-channel models [28], allosteric models, and simple enzyme kinetic models [22, 32]. The phenomenon of very strong rectification has also been frequently observed in stationary-state I–V measurements of ligand-gated ion channels, esp., during whole-cell recording [8, 69, 78]. As a feature of yeast TRK-Cl− currents, therefore, strong rectification provides a further argument in support of TRK oligomerization.

Another important yet unanswered question is the actual biological function of the TRK-anion pathway, viz., the reason for especial sequence conservation in helices M1D and M2D, across fungal species [24, 61]. The short answer is that we do not know. But a significant possibility is that it provides the normal escape route for metablic anions to balance charges during net proton ejection by the plasma-membrane proton pump ([72, 75]). In Neurospora crassa, for example, net proton extrusion rises approximately tenfold as pHo is raised over the range 5–9, with a high pH value >20, and a low value <2 mMol/kg cell water/min [74]. The high value corresponds to pump currents required to sustain membrane voltage; and it is roughly balanced by release of carboxylic acid anions, whereas charge balance at low pH is maintained predominantly by proton reentry coupled to nutrient uptake. The major carboxylic anions released under normal laboratory conditions for Saccharomyces are lactate, malate, and succinate [73]; and for Neurospora are malate, α-ketoglutarate, glutamate, fumarate, citrate, and succinate (C.L. Slayman & P. Kaminski, unpublished data; lactate not measured). Although none of these anions has been tested in the yeast TRK system, formate, acetate, and several inorganic divalent ions proved to be poor substitutes for Cl− [50]. Those trials were carried out at pHo=5.5, however; and because release of carboxylate anions rises steeply with pHo, the intriguing possibility exists that TRK permeability to metabolic divalent ions might actually increase with rising pHo. [Carboxylic-acid transporters, per se, have been only sparingly studied in fungi, and almost always in regard to uptake of those anions as alternative carbon sources [15, 57, 77]. There has been little interest in efflux, outside of industrial processing.]

Acknowledgments

The authors are indebted to Dr. Esther Bashi for technical assistance throughout these experiments and especially for the construction of the EB series of strains, and to Mr. Kenneth Allen for always ready advice and assistance in the management of yeast cultures. We are also grateful to Dr. Albert Smith and to Dr. Kyle Cunningham (Johns Hopkins University) for provision of starting strains; and to Dr. Ming Zhou (Columbia University) for prepublication review of his manuscript on the structure of VpTrkH. Finally, Drs. Fred Sigworth, Edward Moczdlowski (Clarkson University), and Bertl Hille (University of Washington) provided much creative and helpful criticism of the manuscript. The work was supported by Research Grant R01-GM60696 (to CLS) from the US National Institute of General Medical Sciences and by an Overseas Research Scholarship from the Japanese Ministry of Education, Culture, Sports, Science, and Technology (to TK).

15. Cássio F, Leão C. A comparative study on the transort of L (−)malic acid and other short-chain carboxylic acids in the yeast Candida utilis: evidence for a general organic acid permease. Yeast. 1993;9:743–752.[PubMed]

84. Tholema N, Bakker EP, Suzuki A, Nakamura T. Change to alanine of one out of four selectivity filter glycines in KtrB causes a two orders of magnitude decrease in the affinities for both K+ and Na+ of the Na+-dependent K+-uptake system KtrAB from Vibrio alginolyticus. FEBS Lett. 1999;450:217–220.[PubMed]