All-atom molecular dynamics simulations have been applied in the recent past to explore the free energetics underlying ion transport processes in biological ion channels. Roux and co-workers, Kuyucak and coworkers, Busath and coworkers, and others have performed rather elegant and extended timescale molecular dynamics simulations using current state-of-the-art fixed-charge (non-polarizable) force fields in order to calculate the potential of mean force defining the equilibrium flux of ions through prototypical channels such as Gramicidin A. An inescapable conclusion of such studies has been the gross overestimation of the equilibrium free energy barrier, generally predicted to be from 10 – 20 kcal/mole depending on the force field and simulation protocol used in the calculation; this translates to an underestimation of experimentally measurable single channel conductances by several orders of magnitude. Next-generation polarizable force fields have been suggested as possible alternatives for more quantitative predictions of the underlying free energy surface in such systems. Presently, we consider ion permeation energetics in the gramicidin A channel using a novel polarizable force field. Our results predict a peak barrier height of 6 kcal/mole relative to the channel entrance; this is significantly lower than the uncorrected value of 12 kcal/mol for non-polarizable force fields such as GROMOS and CHARMM27 which do not account for electronic polarization. These results provide promising initial indications substantiating the long-conjectured importance of polarization effects in describing ion-protein interactions in narrow biological channels.
展开▼
