Symmetry and asymmetry of permeation through toxin-modified Na+ channels.

Single Na+ channels from rat skeletal muscle were inserted into planar lipid bilayers in the presence of either 200 nM batrachotoxin (BTX) or 50 microM veratridine (VT). These toxins, in addition to their ability to shift inactivation of voltage-gated Na+ channels, may be used as probes of ion conduction in these channels. Channels modified by either of the toxins have qualitatively similar selectivity for the alkali cations (Na+ approximately Li+ greater than K+ greater than Rb+ greater than Cs+). Biionic reversal potentials, for example, were concentration independent for all ions studied. Na+/K+ and Na+/Rb+ reversal potentials, however, were dependent on the orientation of the ionic species with respect to the intra- or extracellular face of the channel, whereas Na+/Li+ biionic reversal potentials were not orientation dependent. A simple, four-barrier, three-well, single-ion occupancy model was used to generate current-voltage relationships similar to those observed in symmetrical solutions of Na, K, or Li ions. The barrier profiles for Na and Li ions were symmetric, whereas that for K ions was asymmetric. This suggests the barrier to ion permeation for K ions may be different than that for Na and Li ions. With this model, these hypothetical energy barrier profiles could predict the orientation-dependent reversal potentials observed for Na+/K+ and Na+/Rb+. The energy barrier profiles, however, were not capable of describing biionic Na/Li ion permeation. Together these results support the hypothesis that Na ions have a different rate determining step for ion permeation than that of K and Rb ions.     

Modeling ion permeation through batrachotoxin-modified Na+ channels from rat skeletal muscle with a multi-ion pore.

The mechanism of ion permeation through Na+ channels that have been modified by batrachotoxin (BTX) and inserted into planar bilayers has been generally described by models based on single-ion occupancy, with or without an influence of negative surface charge, depending on the tissue source. For native Na+ channels there is evidence suggestive of a multi-ion conduction mechanism. To explore the question of ion occupancy, we have reexamined permeation of Na+, Li+, and K+ through BTX-modified Na+ channels from rat skeletal muscle. Single-channel current-voltage (I-V) behavior was studied in neutral lipid bilayers in the presence of symmetrical Na+ concentrations ranging from 0.5 to 3,000 mM. The dependence of unitary current on the mole fraction of Na+ was also examined in symmetrical mixtures of Na(+)-Li+ and Na(+)-K+ at a constant total ionic strength of 206 and 2,006 mM. The dependence of unitary conductance on symmetrical Na+ concentration does not exhibit Michaelis-Menten behavior characteristic of single-ion occupancy but can be simulated by an Eyring-type model with three barriers and two sites (3B2S) that includes double occupancy and ion-ion repulsion. Best-fit energy barrier profiles for Na+, Li+, and K+ were obtained by nonlinear curve fitting of I-V data using the 3B2S model. The Na(+)-Li+ and Na(+)-K+ mole-fraction experiments do not exhibit an anomalous mole-fraction effect. However, the 3B2S model is able to account for the biphasic dependence of unitary conductance on symmetrical [Na+] that is suggestive of multiple occupancy and the monotonic dependence of unitary current on the mole fraction of Na+ that is compatible with single or multiple occupancy. The best-fit 3B2S barrier profiles also successfully predict bi-ionic reversal potentials for Na(+)-Li+ and Na(+)-K+ in both orientations across the channel. Our experimental and modeling results reconcile the dual personality of ion permeation through Na+ channels, which can display features of single or multiple occupancy under various conditions. To a first approximation, the 3B2S model developed for this channel does not require corrections for vestibule surface charge. However, if negative surface charges of the protein do influence conduction, the conductance behavior in the limit of low [Na+] does not correspond to a Gouy-Chapman model of planar surface charge.     

Stochastic theory of ion movement in channels with single-ion occupancy. Application to sodium permeation of gramicidin channels.

The electrodiffusion equations were solved for the one-ion channel both by the analytical method due to Levitt and also by Brownian dynamic simulations. For both types of calculations equilibration of ion distribution between the bath and the ends of the channel was assumed. Potential profiles were found that give good fits to published data on Na+ permeation of gramicidin channels. The data were best fit by profiles that have no relative energy maximum at the mouth of the channel. This finding suggests that alignment of waters or channel charged groups inside the channel in response to an ion's approach may provide an energetically favorable situation for entry sufficient to overcome the energy required for removing bulk waters of hydration. An alternative possibility is that the barrier to ion entry is situated outside the region restricted to single-ion occupancy. Replacement of valine with more polar amino acids at the No. 1 location was found to correspond to a deepening of the potential minima near the channel mouths, an increase in height of the central barrier to ion translocation across the channel, and possibly a reduction in the mobility of the ion-water complex in the channel. The Levitt theory was extended to calculate passage times for ions to cross the channel and the blocking effects of ions that entered the channel but didn't cross. These quantities were also calculated by the Brownian dynamics method.     

Energetics of ion permeation, rejection, binding, and block in gramicidin A from free energy simulations

The rigid force fields currently used in molecular dynamics (MD) simulations of biomolecules are optimized for globular proteins. Whether they can also be used in MD simulations of membrane proteins is an important issue that needs to be resolved. Here we address this issue using the gramicidin A channel, which provides an ideal test case because of the simplicity of its structure and the availability of a wealth of functional data. Permeation properties of gramicidin A can be summarized as "it conducts monovalent cations, rejects anions, and binds divalent cations." Hence, a comprehensive test should consider the energetics of permeation for all three types of ions. To that end, we construct the potential of mean force for K(+), Cl(-), and Ca(2+) ions along the channel axis. For an independent check of the potential-of-mean-force results, we also calculate the free energy differences for these ions at the channel center and binding sites relative to bulk. We find that "rejection of anions" is satisfied but there are difficulties in accommodating the other two properties using the current MD force fields.     

Thallous ion movements through gramicidin channels incorporated in lipid monolayers supported by mercury.

The potential independent limiting flux of hydrated Tl(+) ions through gramicidin (GR) channels incorporated in phospholipid monolayers self assembled on a hanging mercury-drop electrode is shown to be controlled both by diffusion and by a dehydration step. Conversely, the potential independent limiting flux of dehydrated Tl(+) ions stemming from Tl amalgam electro-oxidation is exclusively controlled by diffusion of thallium atoms within the amalgam. Modulating the charge on the polar heads of dioleoylphosphatidylserine (DOPS) by changing pH affects the limiting flux of hydrated Tl(+) ions to a notable extent, primarily by electrostatic interactions. The dipole potential of DOPS and dioleoylphosphatidylcholine (DOPC), positive toward the hydrocarbon tails, does not hinder the translocation step of Tl(+) ions to such an extent as to make it rate limiting. Consequently, incorporation in the lipid monolayer of phloretin, which decreases such a positive dipole potential, does not affect the kinetics of Tl(+) flux through GR channels. In contrast, the increase in the positive dipole potential produced by the incorporation of ketocholestanol causes the translocation step to contribute to the rate of the overall process. A model providing a quantitative interpretation of the kinetics of diffusion, dehydration-hydration, translocation, and charge transfer of the Tl(+)/Tl(0)(Hg) couple through GC channels incorporated in mercury-supported phospholipid monolayers is provided. A cut-off disk model yielding the profile of the local electrostatic potential created by an array of oriented dipoles located in the lipid monolayer along the axis of a cylindrical ion channel is developed.     

The permeation properties of small organic cations in gramicidin A channels.

The conductance properties of organic cations in single gramicidin A channels were studied using planar lipid bilayers. From measurements at 10 mM and at 27 mV the overall selectivity sequence was found to be NH4+ > K+ > hydrazinium > formamidinium > Na+ > methylammonium, which corresponds to Eisenman polyatomic cation sequence X'. Methylammonium and formamidinium exhibit self block, suggesting multiple occupancy and single filing. Formamidinium has an apparent dissociation constant (which is similar to those of alkali metal cations) for the first ion being 22 mM from the Eadie-Hofstee plot (G0 vs. G0/C), 12 mM from the rate constants of a three-step kinetic model. The rate-limiting step for formamidinium is translocation judging from supralinear I-V relations at low concentrations. 1 M formamidinium solutions yields exceptionally long single channel lifetimes, 20-fold longer than methylammonium, which yields lifetimes similar to those found with alkali metal cations. The average lifetime in formamidinium solution significantly decreases with increasing voltage up to 100 mV but is relatively voltage independent between 100 and 200 mV. At lower voltages (< or = 100 mV), the temperature and concentration dependences of the average lifetime of formamidinium were steep. At very low salt concentrations (0.01 M, 100 mV), there was no significant difference in average lifetime from that formed with 0.01 M methylammonium or hydrazinium. We conclude that formamidinium very effectively stabilizes the dimeric channel while inside the channel and speculate that it does so by affecting tryptophan-reorientation or tryptophan-lipid interactions at binding sites.     

Na+ permeation and block of hERG potassium channels

The inactivation gating of hERG channels is important for the channel function and drug-channel interaction. Whereas hERG channels are highly selective for K+, we have found that inactivated hERG channels allow Na+ to permeate in the absence of K+. This provides a new way to directly monitor and investigate hERG inactivation. By using whole cell patch clamp method with an internal solution containing 135 mM Na+ and an external solution containing 135 mM NMG+, we recorded a robust Na+ current through hERG channels expressed in HEK 293 cells. Kinetic analyses of the hERG Na+ and K+ currents indicate that the channel experiences at least two states during the inactivation process, an initial fast, less stable state followed by a slow, more stable state. The Na+ current reflects Na+ ions permeating through the fast inactivated state but not through the slow inactivated state or open state. Thus the hERG Na+ current displayed a slow inactivation as the channels travel from the less stable, fast inactivated state into the more stable, slow inactivated state. Removal of fast inactivation by the S631A mutation abolished the Na+ current. Moreover, acceleration of fast inactivation by mutations T623A, F627Y, and S641A did not affect the hERG Na+ current, but greatly diminished the hERG K+ current. We also found that external Na+ potently blocked the hERG outward Na+ current with an IC50 of 3.5 mM. Mutations in the channel pore and S6 regions, such as S624A, F627Y, and S641A, abolished the inhibitory effects of external Na+ on the hERG Na+ current. Na+ permeation and blockade of hERG channels provide novel ways to extend our understanding of the hERG gating mechanisms.     

Formation of ion channels by a negatively charged analog of gramicidin A.

O-pyromellitylgramicidin is a derivative of gramicidin in which three carboxyl groups are introduced at the terminal hydroxyl end of the peptide. Experiments with artificial lipid membranes indicate that this negatively charged analog forms ion-permeable channels in a way similar to that of gramicidin. If O-pyromellitylgramicidin is added to only one aqueous solution, the membrane conductance remains small, but increases by several orders of magnitude if the same amount is also added to the other side. In accordance with the dimer model of the channel, the membrane conductance under symmetrical conditions is proportional to the square of the aqueous concentration of O-pyromellitylgramicidin over a wide range. The ratio lambdaPG/lambdaG of the single-channel conductance of O-pyromellitylgramicidin to that of gramicidin is close to unity at high ionic strength, but increases more than fivefold at smaller ionic strength (0.01 M). This observation is explained in terms of an electrostatic effect of the fixed negative charges localized near the mouth of the channel. In a mixture of O-pyromellitylgramicidin and gramicidin, unit conductance steps of intermediate size are observed in addition to the conductance steps corresponding to the pure compounds, indicating the formation of hybrid channels. Hybrid channels with preferred orientation may be formed if small amounts of gramicicin and O-pyromellitylgramicidin are added to opposite sides of the membrane. These hybrid channels show a distinct asymmetry in the current-voltage characteristic.     

Structure and dynamics of ion transport through gramicidin A.

Molecular dynamics calculations in which all atoms were allowed to move were performed on a water-filled ion channel of the polypeptide dimer gramicidin A (approximately 600 atoms total) in the head-to-head Urry model conformation. Comparisons were made among nine simulations in which four different ions (lithium, sodium, potassium, and cesium) were each placed at two different locations in the channel as well as a reference simulation with only water present. Each simulation lasted for 5 ps and was carried out at approximately 300 K. The structure and dynamics of the peptide and interior waters were found to depend strongly on the ion tested and upon its location along the pore. Speculations on the solution and diffusion of ions in gramicidin are offered based on the observations in our model that smaller ions tended to lie off axis and to distort the positions of the carbonyl oxygens more to achieve proper solvation and that the monomer-monomer junction was more distortable than the center of the monomer. With the potential energy surface used, the unique properties of the linear chain of interior water molecules were found to be important for optimal solvation of the various ions. Strongly correlated motions persisting over 25 A among the waters in the interior single-file column were observed.     

On the conductance heterogeneity in membrane channels formed by gramicidin A. A cooperative study.

The relative frequency of low-conductance variants of gramicidin A channels in lipid bilayers was determined in parallel experiments in two different laboratories. A common gramicidin stock solution was tested in both labs and, initially, gave rise to significantly different proportions (9% v. 23%) of "mini" channels in the two labs. The lipid and gramicidin solutions were exchanged to identify the source of the difference: When using solutions prepared in lab A (Andersen), lab B (Busath) observed 9% minis, consistent with the original findings in lab A; when using the gramicidin solution prepared in lab B, lab A observed 18% minis, consistent with the original findings in lab B. The experimental apparatus and analysis techniques are therefore not the source of the discrepancy; rather, the difference appears to stem from some factor(s) related to the gramicidin, lipid, and electrolyte solutions. It appears that the mini frequency cannot reflect intrinsic characteristics of the channel-forming peptide, but rather must, at least in part, reflect environmental modulations of channel properties. This has implications for the interpretation of multi-channel experiments on gramicidin A.     

On the molecular basis of ion permeation in the epithelial Na+ channel.

The epithelial Na+ channel (ENaC) is highly selective for Na+ and Li+ over K+ and is blocked by the diuretic amiloride. ENaC is a heterotetramer made of two alpha, one beta, and one gamma 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 beta and gamma subunits, which precedes the second putative transmembrane alpha helix (M2). A residue in the alpha subunit (alphaS589) 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 alpha, beta, and gamma subunits (alphaG587, betaG529, gammaS541) 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 betaG529 to Ser, Cys, or Asp made ENaC partially permeable to K+ and larger ions, similar to the previously reported alphaS589 mutations. We conclude that the residues alphaG587 to alphaS589 and homologous residues in the beta and gamma subunits form the selectivity filter, which tightly accommodates Na+ and Li+ ions and excludes larger ions like K+.     

Radiation inactivation of ion channels formed by gramicidin A. Protection by lipid double bonds and by alpha-tocopherol.

The conductance induced by the channel-forming peptide gramicidin A in lipid membranes is reduced by many orders of magnitude on exposure of the membrane and its aqueous environment to ionizing radiation. This results from an interaction of free radicals of water radiolysis with the tryptophan residues of gramicidin A. The sensitivity of the ion channels towards irradiation is strongly reduced in the presence of either vitamin E or of highly unsaturated lipids. An increase of the D37 dose up to a factor of 50 was found. The phenomena are interpreted via a reduction of the effective concentration of free radicals (such as OH.) in the membrane by reaction with unsaturated fatty acid residues or with vitamin E.     

Energetics of 22Na+ Absorption by Bean-leaf Slices

Sodium absorption by bean-leaf slices from o·1 mM ²²NaCldepends on metabolic energy, but is not enhanced by light. Lightalso does not effect the ATP content of the tissue. Inhibitorsof energy metabolism which decrease the ATP content also inhibit²²Na⁺ absorption. In darkness, anoxia severely affects the ATPcontent as well as ²²Na⁺ absorption; illumination restores bothto almost normal values.     

The gramicidin ion channel: A model membrane protein

The linear peptide gramicidin forms prototypical ion channels specific for monovalent cations and has been extensively used to study the organization, dynamics and function of membrane-spanning channels. In recent times, the availability of crystal structures of complex ion channels has challenged the role of gramicidin as a model membrane protein and ion channel. This review focuses on the suitability of gramicidin as a model membrane protein in general, and the information gained from gramicidin to understand lipid-protein interactions in particular. Special emphasis is given to the role and orientation of tryptophan residues in channel structure and function and recent spectroscopic approaches that have highlighted the organization and dynamics of the channel in membrane and membrane-mimetic media.     

Gating of cardiac Na+ channels in excised membrane patches after modification by alpha-chymotrypsin.

Single cardiac Na+ channels were investigated after intracellular proteolysis to remove the fast inactivation process in an attempt to elucidate the mechanisms of channel gating and the role of slow inactivation. Na+ channels were studied in inside-out patches excised from guinea-pig ventricular myocytes both before and after very brief exposure (2-4 min) to the endopeptidase, alpha-chymotrypsin. Enzyme exposure times were chosen to maximize removal of fast inactivation and to minimize potential nonspecific damage to the channel. After proteolysis, the single channel current-voltage relationship was approximately linear with a slope conductance of 18 +/- 2.5 pS. Na+ channel reversal potentials measured before and after proteolysis by alpha-chymotrypsin were not changed. The unitary current amplitude was not altered after channel modification suggesting little or no effect on channel conductance. Channel open times were increased after removal of fast inactivation and were voltage-dependent, ranging between 0.7 (-70 mV) and 3.2 (-10 mV) ms. Open times increased with membrane potential reaching a maximum at -10 mV; at more positive membrane potentials, open times decreased again. Fast inactivation appeared to be completely removed by alpha-chymotrypsin and slow inactivation became more apparent suggesting that fast and slow inactivation normally compete, and that fast inactivation dominates in unmodified channels. This finding is not consistent with a slow inactivated state that can only be entered through the fast inactivated state, since removal of fast inactivation does not eliminate slow inactivation. The data indicate that cardiac Na+ channels can enter the slow inactivated state by a pathway that bypasses the fast inactivated state and that the likelihood of entering the slow inactivated state increases after removal of fast inactivation.     

Na+ interacting with gramicidin D. A nuclear magnetic resonance study.

23Na nuclei in a milk-white emulsion composed to nonionic surfactant and higher alcohol in saline were characterized by single values of T1 and T2 and a single Larmor frequency. In the presence of small amounts of gramicidin D (Dubos), the relaxations of 23Na were greatly accelerated, and the transverse relaxation was a sum of two decaying exponentials. But only a single T1 was observed; it was roughly equal to the slow T2. The slow T2 accounted for about 40% of the total resonance intensity. The relaxation rates increased linearly with the increase of the gramicidin concentration. The absorption signal consisted of a narrow and a broad line, both centered at the same frequency. The present results suggest that nuclear magnetic resonance spectroscopy is a useful tool for studying the nature of ion-permeable channels of biological membranes, even when the channel has no ionizable groups.     

The relation between ion permeation and recovery from inactivation of ShakerB K+ channels.

We have studied the relation between permeation and recovery from N-type or ball-and-chain inactivation of ShakerB K channels. The channels were expressed in the insect cell line Sf9, by infection with a recombinant baculovirus, and studied under whole cell patch clamp. Recovery from inactivation occurs in two phases. The faster of the two lasts for approximately 200 ms and is followed by a slow phase that may require seconds for completion. The fast phase is enhanced by both permeant ions (K+, Rb+) and by the blocking ion Cs+, whereas the impermeant ions (Na+, Tris+, choline+) are ineffective. The relative potencies are K+ > Rb+ > Cs+ > NH4+ >> Na+ approximately choline+ approximately Tris+. Ion permeation through the channels is not essential for recovery. The results suggest that cations influence the fast phase of recovery by binding in a site with an electrical distance greater than 0.5. Recovery from fast inactivation is voltage-dependent. With Na+, choline+, or Tris+ outside, about 15% of the channels recover in the fast phase (-80 mV), and the other 85% apparently enter a second inactivated state from which recovery is very slow. Recovery in this phase is not influenced by external ions, but is speeded by hyperpolarization.     

Microscopic model for selective permeation in ion channels.

Ionic permeation in the selectivity filter of ion channels is analyzed by a microscopic model based on molecular kinetic theory. The energy and flux equations are derived by assuming that: (a) the selectivity filter is formed by a symmetrical array of carbonyl groups; (b) ion movement is near the axis of the channel; (c) a fraction of water molecules is separated from the ion while it moves across the selectivity filter; (d) the applied voltage drops linearly across the selectivity filter; (e) ions move independently. Energy profiles, single channel conductances, and the degree of hydration of K+ in a hypothetical K+ channel are examined by varying the following microscopic parameters: ion radius and mass, channel radius, number of effective water dipoles, and number of carbonyl groups. The i-V curve is linear up to +/- 170 mV. If the positions of energy maxima and minima are fixed, this linear range is reduced to +/- 50 mV. Channel radius and ion-water interactions are found to be two major channel structural determinants for selectivity sequences. Both radius and mass of an ion are important in selectivity mediated by these interactions. The theory predicts a total of 15 possible kinetic selectivity sequences for alkali cations in ion channels with a single selectivity filter.     

Characterization of micellar-packaged gramicidin A channels.

The effects of heating, on an aqueous gramicidin A lysolecithin system, were examined by carbon-13 nuclear magnetic resonance (¹³C-NMR), circular dichroism (CD), and sodium-23 nuclear magnetic resonance (²³Na-NMR), and the results are collectively interpreted to indicate micellar-packaging of gramicidin channels and cation occupancy in the channel. ¹³C-NMR of the gramicidin-lysolecithin system demonstrates a decrease in mobility of the micellar lipid on heating which is indicative of incorporation of gramicidin into the hydrophobic core of the micelle. A unique and reproducible CD spectrum is obtained for the heat incorporated state. Sodium-23 spin-lattice relaxation times (T₁) demonstrated sodium interaction to be dependent on heat incorporation. The T₁ identified interaction is blocked by silver ion which is known to block sodium transport through the channel in lipid bilayer studies. The temperature dependence of the sodium-23 line width defines an exchange process with an activation energy of 6.8 kcal/mole which is essentially the same as the activation energy reported for transport through the channel in lecithin bilayer studies, and the sodium exchange process is blocked by thallium ion which is also known to block sodium transport through the channel.