Structure of the ion channel peptide antibiotic gramicidin A

The crystal structure of the uncomplexed orthorhombic form of gramicidin A has been determined at 0.86 A resolution. The polypeptide crystallizes from ethanol as a left-handed, double-stranded, antiparallel beta 5.6-helical dimer that is 31 A long and an average of 4.8 A in diameter. The uncomplexed channel does not contain ions or solvent molecules, and its diameter is not uniform but varies from a minimum of 3.85 A to a maximum of 5.47 A. There are three empty cavities in the channel that have a diameter exceeding 5.25 A and appear to be large enough to accommodate water molecules or potassium ions in a chemically reasonable coordination environment. The observed crystal structure does not offer any obvious clues as to why an antiparallel beta 5.6-helix cannot function as an ion channel in lipid bilayers.     

The kinetics of ion movements in the gramicidin channel.

The main features of the ion permeability of gramicidin channels are summarized. The significance of maximums in the single channel conductance-concentration curves, of concentration-dependent permeability ratios, and or current-voltage curves with concentration-dependent form, as well as of other features, is discussed in terms of the mechanism of the ion transfer processes. The observations are then shown to be accounted for by rate theory expressions derived for a model pore consisting of two sites in series and in which ions are not permitted to pass each other. The status of other models is briefly reviewed.     

Modeling negative ion defect migration through the gramicidin A channel

The results of potential of mean force (PMF) calculations for the distinct stages of proton conduction through the gramicidin A channel, including proton migration, reorientation of the water file and negative ion defect migration, are presented. The negative ion defect migration mechanism was hypothesized in experimental studies but was not considered previously in molecular dynamics simulations. The model system consisted of the peptide chains constructed on the base of the structure PDBID:1JNO, the inner file of nine water molecules and external clusters of water molecules placed at both ends of the channel. Potential energy functions were computed with the CHARMM/PM6/TIP3P parameters. The results obtained for proton migration and water file reorientation are basically consistent with those reported previously by Pómès and Roux (Biophys J 82:2304, 2002) within the similar approach. For the newly considered mechanism of negative ion defect migration from the channel center to the end of the water file we obtain the energy 3.8 kcal mol⁻¹ which is not considerably different from the activation energy of water reorientation, 5.4 kcal mol⁻¹. Therefore this mechanism may principally compete for the rate-limiting step in proton conduction in gramicidin.     

Thallium-205 NMR determination of the thermodynamics of the interaction between the thallous ion and gramicidin dimers incorporated into micelles

A study has been made of the thermodynamics of the interaction between the thallous ion and gramicidin dimers incorporated into micelles using thallium-205 NMR spectroscopy. The chemical shift data obtained are interpreted interms of a model in which the dimer has only one tight binding site. The variation of the binding constant over the temperature range 303-323 K is used to determine the changes in enthalpy and entropy of binding giving values of -11.3 kcal/mole and -16 e.u. at 303 K, respectively.     

Structure and dynamics of hydronium in the ion channel gramicidin A.

The effects of the hydronium ion, H(3)0+, on the structure of the ion channel gramicidin A and the hydrogen-bonded network of waters within the channel were studied to help elucidate a possible mechanism for proton transport through the channel. Several classical molecular dynamics studies were carried out with the hydronium in either the center of a gramicidin monomer or in the dimer junction. Structural reorganization of the channel backbone was observed for different hydronium positions, which were most apparent when the hydronium was within the monomer. In both cases the average O-O distance between the hydronium ion and its nearest neighbor water molecule was found to be approximately 2.55 A, indicating a rather strong hydrogen bond. Importantly, a subsequent break in the hydrogen-bonded network between the nearest neighbor and the next-nearest neighbor(approximately 2.7 -3.0 A) was repeatedly observed. Moreover, the carbonyl groups of gramicidin A were found to interact with the charge on the hydronium ion, helping in its stabilization. These facts may have significant implications for the proton hopping mechanism. The presence of the hydronium ion in the channel also inhibits to some degree the reorientational motions of the channel water molecules.     

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.     

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.     

Refinement of the spatial structure of the gramicidin A ion channel]

The spatial structure of the gramicidin A (GA) transmembrane ion-channel was refined on the base of cross-peak volumes measured in NOESY spectra (mixing time tau m = 100 and 200 ms). The refinement methods included the comparison of experimental cross-peak volumes with those calculated for low-energy GA conformations, dynamic averaging of the low-energy conformation set and restrained energy minimization. Accuracy of the spatial structure determination was estimated by the penalty function Fr defined as a root mean square deviation of interproton distances corresponding to the calculated and experimental cross-peak volumes. As the initial conformation we used the right-handed pi 6,3 LD pi 6,3 LD helix established on the base of NMR data regardless of the cross-peak volumes. The conformation is in a good agreement with NOE cross-peak volumes (Fr 0.2 to 0.5 A depending on NOESY spectrum). For a number of NOEs formed by the side chain protons, distances errors were found as much as 0.5-2.0 A. Restrained energy minimization procedure had little further success. However some of these errors were eliminated by the change in torsional angle chi 2 of D-Leu12 and dynamic averaging of the Val7 side chain conformations. Apparently, majority of deviations of the calculated and experimental cross-peak volumes are due to the intramolecular mobility of GA and cannot be eliminated within the framework of rigid globule model. In summary the spatial structure of GA ion-channel can be thought as a set of low-energy conformations, differing by the side chain torsion angles chi 1 Val7 and chi 2 D-Leu4 and D-Leu10 and the orientation of the C-terminal ethanolamine group. Root mean square differences between the atomic coordinates of conformations are in the range of 0.3-0.8 A.     

Redox-regulated ion channel activity of a cysteine-containing gramicidin A analogue

According to recent data, gramicidin A analogues having positively charged amino acid sequences at the C-termini exhibit two types of channel activity in lipid membranes: classical cation-selective channels and large unselective pores. The induction of unselective pores was shown here to strongly depend on the redox state of the membrane-bathing solution, if the gramicidin analogue contained a cysteine residue in the sequence GSGPKKKRKVC attached to the C-terminus. In particular, the addition of H2O2 led to an increase in the transmembrane current and the loss of cationic selectivity on planar bilayer lipid membranes and an increase in the carboxyfluorescein leakage of liposomes. The effect was observed at high concentration of the peptide while was absent at the single-channel level. It was concluded that oxidation led to possible formation of dimers of the peptide, which promoted the formation of large unselective pores.     

Redox-regulated ion channel activity of a cysteine-containing gramicidin A analogue

According to recent data, gramicidin A analogues having positively charged amino acid sequences at the C-termini exhibit two types of channel activity in lipid membranes: classical cation-selective channels and large unselective pores. The induction of unselective pores was shown here to strongly depend on the redox state of the membrane-bathing solution, if the gramicidin analogue contained a cysteine residue in the sequence GSGPKKKRKVC attached to the C-terminus. In particular, the addition of H₂O₂ led to an increase in the transmembrane current and the loss of cationic selectivity on planar bilayer lipid membranes and an increase in the carboxyfluorescein leakage of liposomes. The effect was observed at high concentration of the peptide while was absent at the single-channel level. It was concluded that oxidation led to possible formation of dimers of the peptide, which promoted the formation of large unselective pores.     

Computed energy profiles for ion transport in the gramicidin A channel

Energy profiles have been established for the transport of Na⁺, K⁺ and Cs⁺ through the gramicidin A channel. In Urry’s head-to-headβ _3.3 ^6.3 left-handed helical dimer structure, using a refined methodology for the calculation of intra and intermolecular interactions. The computations show the important role, for the energy profile and the position of the possible binding site, of the flexible ethanolamine chain, whose significance was till now overlooked. The calculations indicate that the barriers at the entrance and at the center of the channel should be in the order Na⁺ XXX K⁺ XXX Cs⁺ but predict also that the energies of the binding sites should be the greatest for Na⁺ and, then, comparable for K⁺ and Cs⁺. The indications concerning the barriers are confirmed by experiment.     

Modulation of gramicidin A open channel lifetime by ion occupancy.

The hypothesis that the gramicidin A channel stability depends on the level of ion occupancy of the channel was used to derive a mathematical model relating channel lifetime to channel occupancy. Eyring barrier permeation models were examined for their ability to fit the zero-voltage conductance, current-voltage, as well as lifetime data. The simplest permeation model required to explain the major features of the experimental data consists of three barriers and four sites (3B4S) with a maximum of two ions occupying the channel. The average lifetime of the channel was calculated from the barrier model by assuming the closing rate constant to be proportional to the probability of the internal channel sites being empty. The link between permeation and lifetime has as its single parameter the experimentally determined averaged lifetime of gramicidin A channels in the limit of infinitely dilute solutions and has therefore no adjustable parameters. This simple assumption that one or more ions inside the channel completely stabilize the dimer conformation is successful in explaining the experimental data considering the fact that this model for stabilization is independent of ion species and configurational occupancy. The model is used to examine, by comparison with experimental data, the asymmetrical voltage dependence of the lifetime in asymmetrical solutions, the effects of blockers, and the effects of elevated osmotic pressure.     

Cadmium and thallous ion permeabilities through lipid bilayer membranes

Cadmium (Cd2+) and thallous ion (Tl+) permeabilities were measured in planar (Mueller-Rudin) lipid bilayer membranes made from diphytanoylphosphatidylcholine in decane. Permeabilities of the electroneutral Cl- complexes, measured with tracers (109Cd and 204Tl), were about 10(-8) cm X s-1 for CdCl2 and 10(-6) cm X s-1 for TlCl. Electrical conductance measurements showed that permeabilities to Cd2+ and Tl+ were approx. 10(-11) cm X s-1, similar to the Na+ permeability. The low permeabilities to both Cd2+ and CdCl2 are consistent with biological studies which suggest that Cd transport and toxicity are protein mediated and correlated with Cd2+, not CdCl2, concentration. However, the low bilayer permeability to Tl+ raises questions about recent reports that Tl+ is a lipid permeable cation in biological membranes and liposomes. An alternative explanation for the lipid permeable behavior of Tl+ is presented, based on the diffusion of TlCl and other complexes of Tl+ with inorganic and organic anions.     

Temperature dependence of single channel currents and the peptide libration mechanism for ion transport through the gramicidin A transmembrane channel

Summary A study of the temperature dependence of gramicidin A conductance of K⁺ in diphytanoyllecithin/n-decane membranes shows the plot of In (single channel conductance) as a function of reciprocal temperature to be nonlinear for the most probable set of conductance, states. These results are considered in terms of a series of barriers, of the dynamics of channel conformation,vis-a-vis the peptide libration mechanism, and of the effect of lipid viscosity on side chain motions again as affecting the energetics of peptide libration.     

Comparison of gramicidin A and gramicidin M channel conductance dispersities

To explore the possible role of Trp side chains in gramicidin channel conductance dispersity, we studied the dispersity of gramicidin M (gM), a gramicidin variant in which all four tryptophan residues are replaced with phenylalanine residues, and its enantiomer, gramicidin M(-) (gM(-)), and compared them to that of gramicidin A (gA). The conductances of highly purified gM and gM(-) were studied in alkali metal solutions at a variety of concentrations and voltages, in seven different types of lipid, and in the presence of detergent. Like gA channels, the most common gM channel conductance forms a narrow band. However, unlike gA channels, where the remaining 5-30% of channel conductances are broadly distributed below (and slightly above) the main band, in gM there is a narrow secondary band with <50% of the main peak conductance. This secondary peak was prominent in NaCl and KCl, but significantly diminished in CsCl and RbCl. Under some conditions, minor components can be observed with conductances yet lower than the secondary peak. Interconversions between the primary conductance state and these yet lower conductance states were observed. The current-voltage relations for both primary and secondary gM channel types have about the same curvature. The mean lifetime of the secondary channel type is below one third that of the primary type. The variants represent state deviations in the peptide or adjacent lipid structure.     

Evaluation of surface tension and ion occupancy effects on gramicidin A channel lifetime.

The surface tension of glycerylmonooleate-hexadecane lipid bilayer membranes and the lifetime of gramicidin A channels were measured at various concentrations of the surrounding solutions. For HCl the surface tension is essentially constant at approximately 5 mN/m up to approximately 1 M, whereas the average lifetime increases approximately 40-fold. At higher concentrations the surface tension decreases markedly. For CsCl the surface tension is constant up to about 1 M then increases with salt level. The average lifetime in this case increases about sixfold. In both cases the lifetime levels off and even decreases at higher salt levels. The increase in lifetime observed with ion activity is therefore qualitatively different from, and not explained by, the established dependence of lifetime on membrane properties (Elliot, J.R., D. Needham, J.P. Dilger, and D.A. Haydon. 1983. Biochim. Biophys. Acta. 735:95-103). We have previously proposed that ion occupancy is a determinant of channel stability, and to test this hypothesis the voltage dependence of channel lifetime was measured in asymmetrical solutions. For the case of a potassium chloride solution on one side of the membrane and a hydrogen chloride solution, on the other, the voltage dependence of the lifetime is asymmetrical. The asymmetry is such that when the electrical field is applied in the direction of the chemical gradient for each of the ions, the channel lifetime approaches, at increasing field strengths, that of a symmetrical solution of the respective ion. The voltage dependence of the surface tension, on the other hand, is negligible for the range of voltages used.(ABSTRACT TRUNCATED AT 250 WORDS)     

Movement of thallous ion across the ascites cell membrane

Summary The movement of thallous ion (Tl⁺ across the ascites cell membrane has been characterized. Analogous to previous findings for⁸⁶Rb⁺ (used as a tracer for K⁺),²⁰⁴Tl⁺-influx could be resolved into three components: a ouabain-inhibitable pump flux, a passive flux, and a furosemide- or NO ₃ ^– -sensitive exchange flux. Although Tl⁺ moved approximately nine times faster across the membrane than K⁺, the pump/leak ratio was equal for the two ions. This suggests that the pump- and leak-pathways share a common rate-limiting step. The exchange mechanism was shown to provide close coupling between the Tl⁺- and K⁺-gradients.     

Molecular dynamics simulations of the gramicidin A-dimyristoylphosphatidylcholine system with an ion in the channel pore region

To investigate the process of ion permeation in an ion channel systematically, we performed molecular dynamics (MD) simulations on a gramicidin A (GA)-phospholipid model system with an ion in the channel pore region. Each of the three types of ions (Ca2+, Na+ Cl-) was placed at five different positions along the channel axis by replacing a water molecule. MD simulations were performed on each system at constant pressure and constant temperature. The MD trajectories showed that the Ca2+ and Na+ ions could stably fluctuate in the pore region, but the Cl- ion was pushed out because of the unfavorable interaction with the channel. This result is consistent with experimental data. It was also found that the conformation of the GA channel underwent a significant change due to the presence of the ion, and the two ends of the GA monomer were more flexible than its middle region. In particular, the dramatic change of local pore radius near the ion indicated this kind of deformation. The strong interaction between the ion and carbonyl oxygen atoms of GA was the major contributor to this change. Furthermore, it was found that the ethanolamine group of the GA molecule was the most flexible group in the GA channel and often observed to block the entrance of GA. These results imply that the deformation of channel structure plays a very important factor in ion permeation, and the ethanolamine group may play a key role in regulating ion entry into the pore. In conclusion, our results indicate that the ion has a dominant influence on the structure of the GA channel and that the flexibility of the ion channel is a crucial factor in the ion permeation process.