This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.105.074633v1 90/11/3941    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Bastug, T. Articles by Kuyucak, S. Search for Related Content PubMed PubMed Citation Articles by Bastug, T. Articles by Kuyucak, S.

Energetics of Ion Permeation, Rejection, Binding, and Block in Gramicidin A from Free Energy Simulations

Turgut Batu and Serdar Kuyucak

School of Physics, University of Sydney, Sydney, Australia

Correspondence: Address reprint requests to Dr. Serdar Kuyucak, Tel.: 61-2-93-51-3162; E-mail: serdar{at}physics.usyd.edu.au.

The rigid force fields currently used in molecular dynamics (MD) simulations of biomolecules are optimized for globular proteins. Whether they can also be used in MD simulations of membrane proteins is an important issue that needs to be resolved. Here we address this issue using the gramicidin A channel, which provides an ideal test case because of the simplicity of its structure and the availability of a wealth of functional data. Permeation properties of gramicidin A can be summarized as "it conducts monovalent cations, rejects anions, and binds divalent cations." Hence, a comprehensive test should consider the energetics of permeation for all three types of ions. To that end, we construct the potential of mean force for K⁺, Cl^–, and Ca²⁺ ions along the channel axis. For an independent check of the potential-of-mean-force results, we also calculate the free energy differences for these ions at the channel center and binding sites relative to bulk. We find that "rejection of anions" is satisfied but there are difficulties in accommodating the other two properties using the current MD force fields.