Ion transport in the gramicidin channel: molecular dynamics study of single and double occupancy.

The structural and thermodynamic factors responsible for the singly and doubly occupied saturation states of the gramicidin channel are investigated with molecular dynamics simulations and free energy perturbation methods. The relative free energy of binding of all of the five common cations Li+, Na+, K+, Rb+, and Cs+ is calculated in the singly and doubly occupied channel and in bulk water. The atomic system, which includes the gramicidin channel, a model membrane made of neutral Lennard-Jones particles and 190 explicit water molecules to form the bulk region, is similar to the one used in previous work to calculate the free energy profile of a Na+ ion along the axis of the channel. In all of the calculations, the ions are positioned in the main binding sites located near the entrances of the channel. The calculations reveal that the doubly occupied state is relatively more favorable for the larger ions. Thermodynamic decomposition is used to show that the origin of the trend observed in the calculations is due to the loss of favorable interactions between the ion and the single file water molecules inside the channel. Small ions are better solvated by the internal water molecules in the singly occupied state than in the doubly occupied state; bigger ions are solvated almost as well in both occupation states. Water-channel interactions play a role in the channel response. The observed trends are related to general thermodynamical properties of electrolyte solutions.     

Water transport and ion-water interaction in the gramicidin channel

The diffuse permeability and the diffusion coefficient of water (Dw) in the gramicidin channel is determined from the osmotic water permeability of the channel and "single file" pore theory. Dw is about 7% of the self-diffusion coefficient of bulk water. The diffusion coefficient of a single water molecule alone in the channel is also determined and is about equal to the value in bulk water. This provides an estimate of the mobility of water on the channel walls in the absence of water-water interaction. Since the gramicidin channel walls should be representative of uncharged polar protein surfaces, this result provides direct evidence that the presence of a cation in the channel reduces the hydraulic water permeability by a factor ranging from 60 for Tl+ to 5 for Na+. The diffusion coefficient of a cation (Dc) in the channel is estimated and compared with Dw. For Na+ it is found that Dc approximately equal to Dw, which implies that the movement of the row of water molecules through the channel determines the local mobility of Na+. Thus, it seems that short range ion-wall interactions are not important in determining the channel conductance for Na+. In contrast, for Li+, local ion-wall interactions probably do limit the conductance.     

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.     

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.     

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.     

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.     

Ion passage pathways and thermodynamics of the amphotericin B membrane channel

Amphotericin B is a polyene macrolide antibiotic used to treat systemic fungal infections. Amphotericin B's chemotherapeutic action requires the formation of transmembrane channels, which are known to transmit monovalent ions. We have investigated the ion passage pathways through the pore of a realistic model structure of the channel and computed the associated thermodynamic properties. Our calculations combined the free energy computations using the Poisson equation with a continuum solvent model and the molecular simulations in which solvent molecules were present explicitly. It was found that there are no substantial structural barriers to a single sodium or chloride ion passage. Thermodynamic free energy calculations showed that the path along which the ions prefer to move is off center from the channel's central axis. In accordance with experiments, Monte Carlo molecular simulations established that sodium ions can pass through the pore. When it encounters a chloride anion in the channel, the sodium cation prefers to form a solvent-bridged pair configuration with the anion.     

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.     

Structure of gramicidin A.

Gramicidin A, a hydrophobic linear polypeptide, forms channels in phospholipid membranes that are specific for monovalent cations. Nuclear Magnetic Resonance (NMR) spectroscopy provided the first direct physical evidence that the channel conformation in membranes is an amino terminal-to-amino terminal helical dimer, and circular dichroism (CD) spectroscopy has shown the sensitivity of its conformation to different environments and the structural consequences of ion binding. The three-dimensional structure of a gramicidin/cesium complex has been determined by x-ray diffraction of single crystals using single wavelength anomalous scattering for phasing. The left-handed double helix in this crystal form corresponds to one of the intermediates in the process of folding and insertion into membranes. Co-crystals of gramicidin and lipid that appear to have gramicidin in their membrane channel conformation have also been formed and are presently under investigation. Hence, we have used a combination of spectroscopic and diffraction techniques to examine the conformation and functionally-related structural features of gramicidin A.     

Voltage-dependent behavior of a "ball-and-chain" gramicidin channel.

The channel-forming properties of two analogs of gramicidin, gramicidin-ethylenediamine (gram-EDA), and gramicidin-N,N-dimethylethylenediamine (gram-DMEDA) were studied in planar lipid bilayers, using protons as the permeant ion. These peptides have positively charged amino groups tethered to their C-terminal ends via a linker containing a carbamate group. Gram-DMEDA has two extra methyl groups attached to the terminal amino group, making it a bulkier derivative. The carbamate groups undergo thermal cis-trans isomerization on the 10-100-ms time scale. The conductance behavior of gram-EDA is found to be markedly voltage dependent, whereas the behavior of gram-DMEDA is not. In addition, voltage affects the cis-trans ratios of the carbamate groups of gram-EDA, but not those of gram-DMEDA. A model is proposed to account for these observations, in which voltage can promote the binding of the terminal amino group of gram-EDA to the pore in a "ball-and-chain" fashion. The bulkiness of the gram-DMEDA derivative prevents this binding.     

The total electrostatic potential in a gramicidin channel

Summary This paper describes a parameter free model of the electrostatic structure of gramicidin channels incorporated into uncharged lipid bilayer membranes. The electrical potential due to all sources is calculated for singly and doubly occupied channels. The model is consistent with all channel properties that are elearly dependent on coulombic interactions. The calculated value of the translocation rate constant and of the binding constant ratio for single and double occupancy are in excellent accord with experiment.     

Ion transport across the chick ileum: a good model for transport studies.

1. The young chick (5-8 days) has been found to be an excellent preparation for the study of transepithelial intestinal ion transport. Due to the thinness of the intestinal tissue, it is not necessary to remove the serosal layers (serosal membranes, circular, and longitudinal muscles), thus circumventing the problems inherent in "stripping" the tissue. 2. The intact chick ileum had a significantly greater short-circuit current (Isc) and lower resistance than did intact adult ileum and transport parameters remained stable over the 6 hr experimental period. 3. Compared to the adult tissue, unidirectional fluxes of Na and Cl were greater in the chick ileum. Net flux of Na (absorption) was about 3 times greater in the chick ileum and the flux was equivalent to the Isc, thus this preparation appears to be characterized by electrogenic Na absorption. 4. Several ileal preparations from day old chicks were studied over an 18 hr period and these preparations were found to remain viable for this period of time with the Isc at the end of 18 hr being nearly identical to that at 2 hr. 5. Besides the advantage of not having to strip the intestinal tissue, and the long-term viability of the tissue, the chick is very inexpensive and easy to obtain and maintain.     

Orientation of gramicidin A transmembrane channel. Infrared dichroism study of gramicidin in vesicles.

Polarized infrared spectroscopy has been used to investigate the orientation of gramicidin A incorporated in dimyristoylphosphatidylcholine liposomes. Dichroism measurements of the major lipid (C = O ester, PO2-, CH2) and peptide (amide A, I, II) bands were performed on liposomes (with or without gramicidin) oriented by air-drying. The mean orientation of the lipid groups and of the pi LD helix chain in the gramicidin has been determined. It can be inferred from infrared frequencies of gramicidin that the dominant conformation of the peptide in liposomes cannot be identified to the antiparallel double-helical dimer found in organic solution. No shift in lipid frequencies was observed upon incorporation of gramicidin in the liposomes. However, a slight reorganization of the lipid hydrocarbon chains which become oriented more closely to the normal to the bilayer is evidenced by a change in the dichroism of the CH2 vibrations. The infrared dichroism results of gramicidin imply a perpendicular orientation of the gramicidin transmembrane channel with the pi LD helix axis at less than 15 degrees with respect to the normal to the bilayer.     

Spin-labeled gramicidin a: channel formation and dissociation

Gramicidin A was studied by continuous wave electron spin resonance (CW-ESR) and by double-quantum coherence electron spin resonance (DQC-ESR) in several lipid membranes (using samples that were macroscopically aligned by isopotential spin-dry ultracentrifugation) and vesicles. As a reporter group, the nitroxide spin-label was attached at the C-terminus yielding the spin-labeled product (GAsl). ESR spectra of aligned membranes containing GAsl show strong orientation dependence. In DPPC and DSPC membranes at room temperature the spectral shape is consistent with high ordering, which, in conjunction with the observed high polarity of the environment of the nitroxide, is interpreted in terms of the nitroxide moiety being close to the membrane surface. In contrast, spectra of GAsl in DMPC membranes indicate deeper embedding and tilt of the NO group. The GAsl spectrum in the DPPC membrane at 35 degrees C (the gel to Pbeta phase transition) exhibits sharp changes, and above this temperature becomes similar to that of DMPC. The dipolar spectrum from DQC-ESR clearly indicates the presence of pairs in DMPC membranes. This is not the case for DPPC, rapidly frozen from the gel phase; however, there are hints of aggregation. The interspin distance in the pairs is 30.9 A, in good agreement with estimates for the head-to-head GAsl dimer (the channel-forming conformation), which matches the hydrophobic thickness of the DMPC bilayer. Both DPPC and DSPC, apparently as a result of hydrophobic mismatch between the dimer length and bilayer thickness, do not favor the channel formation in the gel phase. In the Pbeta and Lalpha phases of DPPC (above 35 degrees C) the channel dimer forms, as evidenced by the DQC-ESR dipolar spectrum after rapid freezing. It is associated with a lateral expansion of lipid molecules and a concomitant decrease in bilayer thickness, which reduces the hydrophobic mismatch. A comparison with studies of dimer formation by other physical techniques indicates the desirability of using low concentrations of GA (approximately 0.4-1 mol %) accessible to the ESR methods employed in the study, since this yields non-interacting dimer channels.     

Effective pore radius of the gramicidin channel. Electrostatic energies of ions calculated by a three-dielectric model.

Electrostatic calculation of the gramicidin channel is performed on the basis of a three-dielectric model in which the peptide backbone of the channel is added as a third dielectric region to the conventional two-dielectric channel model (whose pore radius is often referred to as the effective pore radius reff). A basic principle for calculating electrostatic fields in three-dielectric models is introduced. It is shown that the gramicidin channel has no unique value of reff. The reff with respect to the "self-image energy" (i.e., the image energy in the presence of a single ion) is 2.6-2.7 A, slightly depending upon the position of the ion (the least-square value over the whole length of the pore is 2.6 A). In contrast, the reff with respect to the electric potential due to an ion (and hence the reff with respect to the interaction energy between two ions) is dependent upon the distance s of separation; it ranges from 2.6 to greater than 5 A, increasing with an increase in s. However, for the purpose of rough estimation, the reff with respect to the self-image energy can also be used in calculating the electric potential and the interaction energy, because the error introduced by this approximation is an overestimation of the order of 30% at most. It is also shown that the apparent dielectric constant for the interaction between two charges depends markedly upon the positions of the charges. In the course of this study, the dielectric constant and polarizability of the peptide backbone in the beta-sheet structure is estimated to be 10 and 8.22 A3.     

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.