The gramicidin A channel

Accelerated diffusion of Na+ in a hydrophobic region revealed by molecular dynamics simulations of a synthetic ion channel

To get insight into the significance of the hydrophobic lining on the ion permeation, we performed molecular dynamics simulations on a Na(+) permeation through a de novo synthetic hydrophobic channel. Electrophysiological study has suggested that the channel is formed from a tail-to-tail associated dimer of a cyclic octa-peptide coupled with hydrophobic acyl chains. The acyl chains line the channel pore while the cyclic peptide forms the channel entrance [Z. Qi, M. Sokabe, K. Donowaki, H. Ishida, Biophys. J. 76 (1999) 631]. Molecular dynamics simulation of water in the channel indicated that the inferred structure is physically reasonable [Z. Qi, M. Sokabe, Biophys. Chem. 71 (1998) 35]. In the present study, the potential energy profile of the Na(+) and the energy contributions from each component of the system at different positions along the channel axis were calculated. An energy well instead of a peak is located at the central hydrophobic cavity of the channel, due to its ability of accommodating at least five water molecules to hydrate the ion. Interestingly, the ion diffuses much faster in the hydrophobic acyl chain region, particularly in the central hydrophobic cavity, than it does in the peptide ring region and even surprisingly faster than that in the bulk phase. These results provide a physical basis for an idea that the hydrophobic lining of the K(+) channel [D.A. Doyle, J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, R. MacKinnon, Science 280 (1998) 69] plays an active role to facilitate the ion permeation through the channel pore.

Dynamic properties of individual water molecules in a hydrophobic pore lined with acyl chains: a molecular dynamics study

Recently, a certain class of synthetic molecules has been shown to form ion channels, the pore of which is lined with hydrophobic acyl chains [M. Sokabe, in: F. Oosawa, H. Hayashi, T. Yoshioka (Eds.), Transmembrane Signaling and Sensation, JSSP/VNU Science Press BV, Tokyo, 1984, p. 119; F. Hayashi, M. Sokabe, M. Takagi, K. Hayashi, U. Kishimoto, Biochim. Biophys. Acta, 510 (1978) 305; M.J. Pregel, L. Jullien, J. Canceill, L. Lacombe, J.M. Lehn, J. Chem. Soc. Perkin Trans., 2 (1995) 417; Y. Tanaka, Y. Kobuke, M. Sokabe, Angew. Chem. Int. Ed. Engl., 34 (1995) 693; M. Sokabe, Z. Qi, K. Donowaki, H. Ishida, K. Okubo, Biophys. J., 70 (1996) A201; H. Ishida, K. Donowaki, Y. Inoue, Z. Qi, M. Sokabe, Chem. Lett. (1997) p. 953]. As an initial step towards understanding the physical mechanisms of ion permeation across such a hydrophobic pore, systematic molecular dynamics simulations were performed to investigate dynamic and energetic properties of water molecules inside the pore using a dimer of alanine-N'-acylated cyclic peptide as a channel model. Dynamic energy profiles for water molecules indicated that the energy barrier at the middle region of the pore is approximately 2-3 kcal/mol higher than that in the cap water region which was defined as a vicinity region of the channel entrance. Energetics analyses demonstrated that the mutual interactions among intrapore water molecules are the major factor to give favorable interaction (negative energy contribution) for themselves. The pore, despite being lined with acyl chains, has a favorable van der Waals interaction with intrapore water molecules. These results may help to explain why water-filled channels can be formed by the hydrophobic helices in natural channels.

Solid-state 13C-NMR studies of the effects of sodium ions on the gramicidin A ion channel

End-to-end helical dimers of gramicidin A form transmembrane pores in lipid bilayers, through which monovalent ions may pass. The groups within the peptide that interact with these ions have been studied by application of solid-state spectroscopic methods to a series of gramicidin A analogues synthesized with 13C in selected peptide carbonyl groups. The resonances of D-Leu10, D-Leu12 and D-Leu14 analogues were perturbed in the presence of 0.16 M sodium ions, whereas the resonances of the carbonyls of Gly2, Ala3, D-Leu4 and Val7, which are closer to the formylated N-terminal end of the peptide, were unaffected. The observed changes in chemical shift anisotropy are indicative of a change in orientation of the abovementioned leucine carbonyls.

Energy profile of Cs+ in gramicidin A in the presence of water. Problem of the ion selectivity of the channel

The effect of water present at the mouth and inside the channel of Gramicidin A on the energy profile calculated for a caesium ion is determined. The total optimal interaction energy computed for the system GA-Cs+-(22 waters) leads to an energy profile characterized by a deep minimum at 11 A followed by an entrance energy barrier of 7 Kcal/mol expanding until 9 A from the center. After this point, a second minimum less deep than the previous one is observed, itself followed by a central barrier. The shape of the profile at the entrance is governed by the balance between the progressive desolvation process of the ion and the increase of favorable hydrogen bond interactions implying both the water molecules and GA. The comparison of this energy profile with that obtained in vacuo shows that the presence of water molecules does not modify the pathway of the ion which, owing to its size, is constrained essentially to remain on the channel axis. The comparison Na+ versus Cs+ indicates that although the phenomena involved are globally the same, differences between the two profiles appear due firstly to the difference in the affinity of the two ions for water and secondly to their respective size. This last difference implies that the number of water molecules present in the interior of the channel during the cation progression is reduced roughly by one in the case of caesium. The desolvation barrier computed for Cs+ is half the corresponding value for Na+, a result in agreement with the observed selectivity.

Energy profiles in the gramicidin A channel

Gramicidin A (GA) is a linear pentadecapeptide made of alternating D and L residues, in which the N-and C-terminals are blocked by a formyl group (head) and an ethanolamine end (tail), respectively (Sarges & Witkop, 1964):

Experimental and theoretical study of gramicidin P, an analog of gramicidin A with a methylamine C-terminal

Gramicidin P (a gramicidin A in which the ethanolamine C-terminus is replaced by methylamine) was synthesized and shown to have the same single-channel conductance behavior as gramicidin A. The results are discussed in connection with the energy profile computed in the presence of water in comparison with the corresponding profile for gramicidin A.