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

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.