RESUMO
Molecular Dynamics (MD) simulations were performed on a gramicidin A dimer model representing a transmembrane channel. Different from previous simulations the peptide was in contact with bulk water at both ends of the dimer to guarantee a realistic description of the hydration of the biomolecule. The flexible BJH model for water was employed in the simulations and the gramicidin-water, gramicidin-ion and ion-water potentials used are based on molecular orbital calculations. The water structure near the gramicidin was investigated first by a simulation without ions, while for the energy profiles of the ion transport through the channel a potassium or a sodium ion was added. These investigations provide a detailed and conclusive picture on a molecular level of the role of water in the ion transport through a gramicidin A channel and can explain the experimental results on the selectivity between alkali ions, their double or even triple occupancy, the exclusion or permeability of anions depending upon cation concentration and the consequences of differences in the ionic charge. The investigation demonstrate that the water molecules around the gramicidin behave as an integral part of the peptide and the functionality is the result of the whole complex biomolecule-water.
