ION-TRANSPORT IN THE GRAMICIDIN CHANNEL - FREE-ENERGY OF THE SOLVATEDRIGHT-HANDED DIMER IN A MODEL MEMBRANE

The free energy profile of a Na+ ion along the axis of the gramicidin A dimer channel is calculated by molecular dynamics simulations. A cyl indrical system is used that is based on a detailed atomic model of th e gramicidin dimer, which is solvated at both ends of the channel and surrounded by a simplified model membrane. The entire system consists of 967 atoms and is treated with cylindrical periodic boundary conditi ons. The gramicidin channel was modeled from the right-handed head-to- head beta-helix dimer determined experimentally by two-dimensional NMR . Simulations were made to estimate the free energy profile in three r egions: the transition from bulk-like solvation to single file, the si ngle-file region at the entrance of the channel, and the intermonomer junction of the dimer channel. In each region, approximate corrections are introduced to account for hysteresis effects in the simulations. Binding sites are found at the extremities of the channel in agreement with experiment. The principal energy barrier is at the entrance of t he channel in the single-rile region. No large activation energy barri er is introduced by the dehydration process, because the transformatio n from bulk to single-file solvation is progressive and takes place ov er a distance of 5-6 angstrom. The free energy simulation results obta ined in this paper are joined with those calculated from a periodic mo del of the interior of the pore to construct an approximate free energ y profile of the entire channel system.