Location of ion-binding sites in the gramicidin channel by X-ray diffraction.

We report the first X-ray diffraction on gramicidin in its membrane-active form by using uniformly aligned multilayer samples of membranes containing gramicidin and ions (T1+, K+, Ba2+, Mg2+ or without ions). From the difference electron density profiles, we found a pair of symmetrically located ion-binding sites for T1- at 9.6 (+/- 0.3) A and for Ba2+ at 13.0 (+/- 0.2) A from the midpoint of the gramicidin channel. The location of Ba(2+)-binding sites is near the ends of the channel, consistent with the experimental observation that divalent cations do not permeate but block the channel. The location of T1(+)-binding sites is somewhat of a surprise. It was generally thought that monovalent cations bind to the first turn of the helix from the mouth of the channel. (It is now generally accepted that the gramicidin channel is a cylindrical pore formed by two monomers, each a single-stranded beta 6.3 helix and hydrogen-bonded head-to-head at their N termini.) But our experiment shows that the T1(+)-binding site is either near the bottom of or below the first helix turn.     

Gramicidin channel selectivity. Molecular mechanics calculations for formamidinium, guanidinium, and acetamidinium.

Empirical energy function calculations were used to evaluate the effects of minimization on the structure of a gramicidin A channel and to analyze the energies of interaction between three cations (guanidinium, acetamidinium, formamidinium) and the channel as a function of position along the channel axis. The energy minimized model of the gramicidin channel, which was based on the results of Venkatachalam and Urry (1983), has a constriction at the channel entrance. If the channel is not allowed to relax in the presence of the ions (rigid model), there is a large potential energy barrier for all three cations. The barrier varies with cation size and is due to high van der Waals and ion deformation energies. If the channel is minimized in the presence of the ions, the potential energy barrier to formamidinium entry is almost eliminated, but a residual barrier remains for guanidinium and acetamidinium. The residual barrier is primarily due, not to the expansion of the helix, but, to the disruption of hydrogen bonds between the terminal ethanoloamine and the next turn of the helix which occurs when the carbonyls of the outer turn of the helix librate inward toward the ion as it enters the channel. The residual potential energy barriers could be a possible explanation for the measured selectivity of gramicidin for formamidinium over guanidinium. The results of this full-atomic model address the applicability of the size-exclusion concept for the selectivity of the gramicidin channel.     

Binding of alkaline cations to the double-helical form of gramicidin.

Gramicidin is a polypeptide antibiotic that forms monovalent cation-specific channels in membrane environments. In organic solvents and in lipids containing unsaturated fatty acid chains, it forms a double-helical "pore" structure, in which two monomers are intertwined. This form of gramicidin can bind two cations inside its lumen, and the crystal structures of both an ion complex and an ion-free form have been determined. In this study, we have used circular dichroism (CD) spectroscopy to examine the binding mechanism and the binding constants (K1 and K2) of cations to gramicidin in the double helical form in methanol solution. The dramatic change in optical rotation in the far-ultraviolet CD spectrum of gramicidin provides a useful tool for monitoring the binding. The binding mechanism appears to involve a large conformation change associated with the binding of ions to the first of the two sites. The calculated values for the K1 binding constants for alkaline cations are considerably smaller than the K2 binding constants. The order of binding affinity for alkaline cations is similar to that for the helical dimer "channel" form of gramicidin, i.e., Cs+ approximately Rb+ > > K+ > Li+, but in comparison to the helical dimer form, the binding to double-helical dimers is dominated by a cation size-dependent conformational change in the gramicidin structure.     

Gramicidin A directly inhibits mammalian Na+/K+-ATPase

The linear pentadecapeptide gramicidin A forms an ion channel in the lipid bilayer to selectively transport monovalent cations. Nevertheless, we have surprisingly found that gramicidin A directly inhibits mammalian Na(+)/K(+)-ATPase. Gramicidin A inhibited ATP hydrolysis by Na(+)/K(+)-ATPase from porcine cerebral cortex at the IC(50) value of 8.1 microM, while gramicidin S was approximately fivefold less active. The synthetic gramicidin A analog lacking N-terminal formylation and C-terminal ethanolamine exhibited a weaker inhibitory effect on the ATP-hydrolyzing activity of Na(+)/K(+)-ATPase than gramicidin A, indicating that these end modifications are necessary for gramicidin A to inhibit Na(+)/K(+)-ATPase activity. Moreover, Lineweaver-Burk analysis showed that gramicidin A exhibits a mixed type of inhibition. In addition to the most well-studied ionophore activity, our present study has disclosed a novel biological function of gramicidin A as a direct inhibitor of mammalian Na(+)/K(+)-ATPase activity.     

Use of weak acids to determine the bulk diffusion limitation of H+ ion conductance through the gramicidin channel.

The addition of 2 M formic acid at pH 3.75 increased the single channel H+ ion conductance of gramicidin channels 12-fold at 200 mV. Other weak acids (acetic, lactic, oxalic) produce a similar, but smaller increase. Formic acid (and other weak acids) also blocks the K+ conductance at pH 3.75, but not at pH 6.0 when the anion form predominates. This increased H+ conductance and K+ block can be explained by formic acid (HF) binding to the mouth of the gramicidin channel (Km = 1 M) and providing a source of H+ ions. A kinetic model is derived, based on the equilibrium binding of formic acid to the channel mouth, that quantitatively predicts the conductance for different mixtures of H+, K+, and formic acid. The binding of the neutral formic acid to the mouth of the gramicidin channel is directly supported by the observation that a neutral molecule with a similar structure, formamide (and malonamide and acrylamide), blocks the K+ conductance at pH 6.0. The H+ conductance in the presence of formic acid provides a lower bound for the intrinsic conductance of the gramicidin channel when there is no diffusion limitation at the channel mouth. The 12-fold increase in conductance produced by formic acid suggests that greater than 90% of the total resistance to H+ results from diffusion limitation in the bulk solution.     

Stimulation of cation transport in mitochondria by gramicidin and truncated derivatives

Gramicidin and the truncated derivatives desformylgramicidin (desfor) and des(formylvalyl)gramicidin (desval) stimulate monovalent cation transport in rat liver mitochondria. Cation fluxes were compared indirectly from the effect of cations on the membrane potential at steady state (state 4) or from the associated stimulation of electron transport. Rb+ transport was measured directly from the uptake of 86Rb. The truncated gramicidins show enhanced selectivity for K+ and Rb+ when compared to gramicidin. Moreover, the pattern of selectivity within the alkali cation series is altered, i.e., Rb+ greater than K+ greater than Cs+ greater than Na+ greater than Li+ for desfor and desval as compared to Cs+ greater than Rb+ greater than K+ = Na+ greater than Li+ for gramicidin. The cation fluxes through the truncated derivatives are more strongly dependent on the cation concentration. The presence of high concentrations of permeating cation enhances the transport of other cations through the truncated derivative channels, suggesting that cations are required for stabilizing the channel structure. In high concentrations of KCl, desfor and desval are nearly as effective as gramicidin in collapsing the mitochondrial membrane potential, and, consequently, in the uncoupling of oxidative phosphorylation and enhancement of ATP hydrolysis. Preliminary experiments with liposomes show that 86Rb exchange is stimulated by desfor and desval almost to the same extent as gramicidin. These results strongly suggest that the truncated gramicidins form a novel conducting channel which differs from the gramicidin head-to-head, single-stranded beta 6.3-helical dimer ("channel") in its conductance characteristic and its structure.(ABSTRACT TRUNCATED AT 250 WORDS)     

Open channel noise. V. Fluctuating barriers to ion entry in gramicidin A channels.

We have measured the fluctuations in the current through gramicidin A (GA) channels in symmetrical solutions of monovalent cations of various concentrations, and compared the spectral density values with those computed using E. Frehland's theory for noise in discrete transport systems (Frehland, E. 1978. Biophys. Chem. 8:255-265). The noise for the transport of NH4+ and Na+ ions in glycerol-monooleate/squalene membranes could be accounted for entirely by "shot noise" in the process of transport through a single-filing pore with two ion binding sites. However, in confirmation of results in a previous paper (Sigworth, F. J., D. W. Urry, and K. U. Prasad. 1987. Biophys. J. 52:1055-1064) currents of Cs+ showed a substantial excess noise at low ion concentrations, as did currents of K+ and Rb+. The excess noise was increased in thicker membranes. The observations are accounted for by a theory that postulates fluctuations of the entry rates of ions into the channel on a time scale of approximately 1 microsecond. These fluctuations occur preferentially when the channel is empty; the presence of bound ions stabilizes the "high conductance" conformation of the channel. The fluctuations are sensed to different degrees by the various ion species, and their kinetics depend on membrane thickness.     

Energy profile of a channel formed by alpha-latrotoxin in a bilayer phospholipid membrane

Current-voltage characteristics are obtained for a channel formed by alpha-latrotoxin when inserting into a bilayer lipid membrane separating solutions of different ionic composition. They are used for determining parameters of a two-barrier model of this channel energy profile. It is shown that selectivity of these channels is based on the same principles that in the calcium channels of biological membranes and is mainly determined by the ion binding inside the channel. Affinity of the channel for the penetrating ions of alkali-earth metals decreases in the sequence: Mg2+ greater than Ca2+ greater than Sr2+ approximately equal to Ba2+ and the blocking ability of the cations of transition metals increases in the series: Mn2+ less than Zn2+ approximately less than Ni2+ approximately less than Co2+ less than Cd2+ much less than La3+. The channel gives the monovalent ions to pass through as well, its permeability being dependent on the concentration of divalent ions from the cis-, but not from the trans-side of the membrane.     

Investigations of the metal ion-binding sites of yeast inorganic pyrophosphatase

Yeast inorganic pyrophosphatase was found to bind two Mn2+ per subunit in the absence of phosphate and three Mn2+ per subunit in the presence of phosphate. Kinetic studies of the pyrophosphatase-catalyzed hydrolysis of Cr(NH3)4PP and Cr(H2O)4PP were carried out with Mn2+ and with Mg2+ as activators. The results from these studies suggest that three divalent cations per pyrophosphatase active site are required for catalysis. NMR and EPR studies were conducted to evaluate the relative location of the metal ion binding sites on the enzyme. The two Mn2+ ions bound to the free enzyme are in close enough proximity to magnetically interact. Analysis of the NMR and EPR data in terms of a dipolar relaxation mechanism between Mn2+ ions provides an estimate of the distance between them of 10-14 A. When the diamagnetic substrate analog [Co(NH3)4PNP]- or intermediate analog [Co(NH3)4 (P)2]- are bound to pyrophosphatase, two Mn2+ ions still bind to the enzyme and their magnetic interaction increases. In the presence of these Co3+ complexes, the Mn2+--Mn2+ separation decreases to 7-9 A. Several NMR and EPR experiments were conducted at low Mn2+ to pyrophosphatase ratios (approximately 0.3), where only one Mn2+ ion binds per subunit, in the presence of Cr3+ or Co3+ complexes of PNP or PP. Analysis of the Mn2+--Cr3+ dipolar relaxation evident in proton NMR and EPR data provided for the calculation of Mn2+--Cr3+ distances. When the substrate analog CrPNP was present, the Mn2+--Cr3+ distance was congruent to 7 A whereas, when Cr(P)2 was bound to pyrophosphatase, the Mn2+--Cr3+ distance was congruent to 5 A. These results strongly support a model for the catalytic site of pyrophosphatase that involves three metal ion cofactors.     

Ion-specific diffusion rates through transmembrane protein channels. A molecular dynamics study

The migration of different alkali metal cations through a transmembrane model channel is simulated by means of the molecular dynamics technique. The parameters of the model are chosen in close relation to the gramicidin A channel. Coulomb- and van der Waals-type potentials between the ions and flexible carbonyl groups of the pore-forming molecule are used to describe the ion channel interaction. The diffusion properties of the ions are obtained from three-dimensional trajectory calculations. The diffusion rates for the different ions Li+, Na+, K+ and Rb+ are affected not only by the mass of the particles but also very strongly by their size. The latter effect is more pronounced for rigid channels, i.e., for binding vibrational frequencies of the CO groups with v greater than 400 cm-1. In this range the selectivity sequence for the diffusion rates is the inverse of that expected from normal rate theory but agrees with that found in experiments for gramicidin A.     

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.     

Thermodynamic parameters for the binding of divalent cations to gramicidin A incorporated into a lipid environment by Tl-205 nuclear magnetic resonance.

Thermodynamic parameters, enthalpy and entropy, for the binding of the divalent cations, Mg+2, Ca+2, Sr+2, Ba+2, and Cd+2, to gramicidin A, incorporated into lysophosphatidylcholine, have been determined using a combination of Tl-205 nuclear magnetic resonance spectroscopy and competition binding. The binding process is thermodynamically driven by the enthalpy and not the entropy. The enthalpy values are related to the process involving the transfer of cations from an aqueous environment to an amide environment. A comparison is made between the thermodynamic parameters for the binding of monovalent and divalent cations to gramicidin A to illustrate the channel blocking ability of the divalent cations with respect to monovalent cation transport.     

Effects of calcium on the gramicidin A single channel in phosphatidylserine membranes

In phosphatidylserine membranes the decrease in the conductance of the gramicidin A single channel caused by calcium is attributed to a reduction of surface potential and to a direct blocking of the pore (Apell et al. 1979).The aim of this paper is to make a, quantitative evaluation of these two effects. We recorded the conductance of gramicidin single channels in 100 mM KCl in the presence of different amounts of CaCl₂, MgCl₂ or TEACl.The ionic activities at the channel mouth were calculated using the Gouy-Chapman-Stern theory. Our experiments showed that even when the K⁺ activity at the channel mouth was estimated to be the same, the single channel conductance was lower if divalent cations were present. This effect is attributed to a blocking action of these ions.Abbreviations PS phosphatidylserine - TEA tetraethylammonium     

Activation by divalent cations of a Ca2+-activated K+ channel from skeletal muscle membrane

Several divalent cations were studied as agonists of a Ca2+-activated K+ channel obtained from rat muscle membranes and incorporated into planar lipid bilayers. The effect of these agonists on single-channel currents was tested in the absence and in the presence of Ca2+. Among the divalent cations that activate the channel, Ca2+ is the most effective, followed by Cd2+, Sr2+, Mn2+, Fe2+, and Co2+. Mg2+, Ni2+, Ba2+, Cu2+, Zn2+, Hg2+, and Sn2+ are ineffective. The voltage dependence of channel activation is the same for all the divalent cations. The time-averaged probability of the open state is a sigmoidal function of the divalent cation concentration. The sigmoidal curves are described by a dissociation constant K and a Hill coefficient N. The values of these parameters, measured at 80 mV are: N = 2.1, K = 4 X 10(-7) mMN for Ca2+; N = 3.0, K = 0.02 mMN for Cd2+; N = 1.45, K = 0.63 mMN for Sr2+; N = 1.7, K = 0.94 mMN for Mn2+; N = 1.1, K = 3.0 mMN for Fe2+; and N = 1.1 K = 4.35 mMN for Co2+. In the presence of Ca2+, the divalent cations Cd2+, Co2+, Mn2+, Ni2+, and Mg2+ are able to increase the apparent affinity of the channel for Ca2+ and they increase the Hill coefficient in a concentration-dependent fashion. These divalent cations are only effective when added to the cytoplasmic side of the channel. We suggest that these divalent cations can bind to the channel, unmasking new Ca2+ sites.     

The effects of calcium ions on double helical forms of gramicidin

The effects of binding calcium ions to the double helical forms of gramicidin present in methanol solution were examined using circular dichroism (CD) and nuclear magnetic resonance (NMR) spectroscopies. It was found that calcium ions principally alter the relative composition of the equilibrium mixture of gramicidin conformers present in the solvent. In the absence of calcium, both parallel and antiparallel double helices are present. However, the addition of small amounts of Ca²⁺ shifts the equilibrium towards the left-handed parallel double helical form. This conformational change prevents monovalent cations (caesiums) from binding to the gramicidin double helix, and even converts the shorter, wider anti-parallel double helical form normally produced in the presence of caesium into the longer, narrower parallel double helical form. Furthermore, a temperature study showed that calcium ions tend to stabilize this form relative to the ion-free forms. The conformation of gramicidin is further changed, becoming a disordered structure, when the concentration of Ca²⁺ is raised. Thus, the binding of divalent calcium ions has a number of dramatic effects on the conformations of gramicidin present in solution.     

Streaming potentials in gramicidin channels measured with ion-selective microelectrodes.

Streaming potentials have been measured for gramicidin channels with a new method employing ion-selective microelectrodes. It is shown that ideally ion-selective electrodes placed at the membrane surface record the true streaming potential. Using this method for ion concentrations below 100 mM, approximately seven water molecules are transported whenever a sodium, potassium, or cesium ion, passes through the channel. This new method confirms earlier measurements (Rosenberg, P.A., and A. Finkelstein. 1978. Interaction of ions and water in gramicidin A channels. J. Gen. Physiol. 72:327-340) in which the streaming potentials were calculated as the difference between electrical potentials measured in the presence of gramicidin and in the presence of the ion carriers valinomycin and nonactin.     

Equilibrium binding constants for Tl+ with gramicidins A, B and C in a lysophosphatidylcholine environment determined by 205Tl nuclear magnetic resonance spectroscopy.

Nuclear Magnetic Resonance (NMR) 205Tl spectroscopy has been used to monitor the binding of Tl+ to gramicidins A, B, and C packaged in aqueous dispersions of lysophosphatidylcholine. For 5 mM gramicidin dimer in the presence of 100 mM lysophosphatidylcholine, only approximately 50% or less of the gramicidin appears to be accessible to Tl+. Analysis of the 205Tl chemical shift as a function of Tl+ concentration over the 0.65-50 mM range indicates that only one Tl+ ion can be bound by gramicidin A, B, or C under these experimental conditions. In this system, the Tl+ equilibrium binding constant is 582 +/- 20 M-1 for gramicidin 1949 +/- 100 M-1 for gramicidin B, and 390 +/- 20 M-1 for gramicidin C. Gramicidin B not only binds Tl+ more strongly but it is also in a different conformational state than that of A and C, as shown by Circular Dichroism spectroscopy. The 205Tl NMR technique can now be extended to determinations of binding constants of other cations to gramicidin by competition studies using a 205Tl probe.     

Ion transport in a model gramicidin channel. Structure and thermodynamics.

The potential of mean force for Na+ and K+ ions as a function of position in the interior of a periodic poly(L,D)-alanine model for the gramicidin beta-helix is calculated with a detailed atomic model and realistic interactions. The calculated free energy barriers are 4.5 kcal/mol for Na+ and 1.0 kcal/mol for K+. A decomposition of the free energy demonstrates that the water molecules make a significant contribution to the free energy of activation. There is an increase in entropy at the transition state associated with greater fluctuations. Analysis reveals that the free energy profile of ions in the periodic channel is controlled not by the large interaction energy involving the ion but rather by the weaker water-water, water-peptide and peptide-peptide hydrogen bond interactions. The interior of the channel retains much of the solvation properties of a liquid in its interactions with the cations. Of particular importance is the flexibility of the helix, which permits it to respond to the presence of an ion in a fluidlike manner. The distortion of the helix is local (limited to a few carbonyls) because the structure is too flexible to transmit a perturbation to large distances. The plasticity of the structure (i.e., the property to deform without generating a large energy stress) appears to be an essential factor in the transport of ions, suggesting that a rigid helix model would be inappropriate.     

Formamidinium-induced dimer stabilization and flicker block behavior in homo- and heterodimer channels formed by gramicidin A and N-acetyl gramicidin A.

Compared to the N-formyl gramicidin A (GA), the N-acetyl gramicidin A (NAG) channel has unchanged conductance in 1 M NH4+ (gamma NN/gamma GG = 1, conductance ratio) but reduced conductance in 1 M K+ (gamma NN/gamma GG = 0.6) methylammonium (gamma NN/gamma GG = 0.3), and formamidinium (gamma NN/gamma GG = 0.1) solutions. Except with formamidinium, "flicker blocks" are evident even at low cutoff frequencies. For all cations studied, channel lifetimes of N-acetyl homodimers (NN) are approximately 50-fold shorter than those of the GA homodimer (GG). The novel properties of GA channels in formamidinium solution (supralinear current-voltage relations and dimer stabilization (Seoh and Busath, 1993)) also appear in NN channels. The average single channel lifetime in 1 M formamidinium solution at 100 mV is 6-7-fold longer than in K+ and methylammonium solutions and, like in the GA channel, significantly decreases with increasing membrane potential. Experiments with mixtures of the two peptides, GA and NAG, showed three main conductance peaks. Oriented hybrids were formed utilizing the principle that monomers remain in one leaflet of the bilayer (O'Connell et al., 1990). With GA at the polarized side and NAG at the grounded side, at positive potentials (in which case hybrids were designated GN) and at negative potentials (in which case hybrids were designated NG), channels had the same conductances and channel properties at all potentials studied. Flicker blocks were not evident in the hybrid channels, which suggests that both N-acetyl methyl groups at the junction of the dimer are required to cause flickers. Channel lifetimes in hybrids are only approximately threefold shorter than those of the GG channels, and channel conductances are similar to those of GG rather than NN channels. We suggest that acetyl-acetyl crowding at the dimeric junction in NN channels cause dimer destabilization, flickers, and increased selectivity in N-acetyl gramicidin channels.     

Effects of different cations on the hydrodynamic radius of DNA.

The effects of different cations on the hydrodynamic radius (RH) of a 48-bp synthetic DNA are measured by time-resolved fluorescence polarization anisotropy of intercalated ethidium. Relative statistical errors in RH are only approximately 1%. With increasing cation concentration, Na+ causes a small decrease in RH, Cs+ causes a somewhat larger decrease by up to 0.5 A at 100 mM, and (CH3CH2)4N+ causes an increase in RH by approximately 0.5 A at 100 mM. The qualitatively different effects of these monovalent cations indicates that the changes in RH with cation concentration do not arise primarily from electrolyte friction. Divalent cations cause much larger increases in RH with increasing cation concentration. Mg2+ causes an increase in RH by up to 1.0 A at 24.4 mM, and Mn2+ causes an increase in RH by up to 1.6 A at 24.4 mM. These effects are independent of DNA concentration. There is some positive correlation between the order of effects of the different cations on RH and the order of their effects on interhelical hydration forces. It is suggested that these different ions affect RH either by altering the hydration layer or possibly by some effect on DNA structure, such as stabilizing bends.