This Article Full Text Full Text (PDF) 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 HighWire Citing Articles via Google Scholar Google Scholar Articles by Mamonov, A. B. Articles by Kurnikova, M. G. Search for Related Content PubMed PubMed Citation Articles by Mamonov, A. B. Articles by Kurnikova, M. G.

The Role of the Dielectric Barrier in Narrow Biological Channels: A Novel Composite Approach to Modeling Single-Channel Currents

Artem B. Mamonov ^*, Rob D. Coalson ^*, Abraham Nitzan  and Maria G. Kurnikova 

^* Chemistry Department, University of Pittsburgh, Pittsburgh, Pennsylvania 15260; School of Chemistry, University of Tel Aviv, Israel; and Chemistry Department, Marquette University, Milwaukee, Wisconsin 53202

Correspondence: Address reprint requests to Maria G. Kurnikova, Todd Wehr Chemistry Building, Marquette University, Milwaukee, WI 53233. Tel.: 414-288-3538; E-mail: maria.kurnikova{at}marquette.edu.

A composite continuum theory for calculating ion current through a protein channel of known structure is proposed, which incorporates information about the channel dynamics. The approach is utilized to predict current through the Gramicidin A ion channel, a narrow pore in which the applicability of conventional continuum theories is questionable. The proposed approach utilizes a modified version of Poisson-Nernst-Planck (PNP) theory, termed Potential-of-Mean-Force-Poisson-Nernst-Planck theory (PMFPNP), to compute ion currents. As in standard PNP, ion permeation is modeled as a continuum drift-diffusion process in a self-consistent electrostatic potential. In PMFPNP, however, information about the dynamic relaxation of the protein and the surrounding medium is incorporated into the model of ion permeation by including the free energy of inserting a single ion into the channel, i.e., the potential of mean force along the permeation pathway. In this way the dynamic flexibility of the channel environment is approximately accounted for. The PMF profile of the ion along the Gramicidin A channel is obtained by combining an equilibrium molecular dynamics (MD) simulation that samples dynamic protein configurations when an ion resides at a particular location in the channel with a continuum electrostatics calculation of the free energy. The diffusion coefficient of a potassium ion within the channel is also calculated using the MD trajectory. Therefore, except for a reasonable choice of dielectric constants, no direct fitting parameters enter into this model. The results of our study reveal that the channel response to the permeating ion produces significant electrostatic stabilization of the ion inside the channel. The dielectric self-energy of the ion remains essentially unchanged in the course of the MD simulation, indicating that no substantial changes in the protein geometry occur as the ion passes through it. Also, the model accounts for the experimentally observed saturation of ion current with increase of the electrolyte concentration, in contrast to the predictions of standard PNP theory.

This article has been cited by other articles:

				D. Boda, W. Nonner, M. Valisko, D. Henderson, B. Eisenberg, and D. Gillespie Steric Selectivity in Na Channels Arising from Protein Polarization and Mobile Side Chains Biophys. J., September 15, 2007; 93(6): 1960 - 1980. [Abstract] [Full Text] [PDF]

				H. Chen, B. Ilan, Y. Wu, F. Zhu, K. Schulten, and G. A. Voth Charge Delocalization in Proton Channels, I: The Aquaporin Channels and Proton Blockage Biophys. J., January 1, 2007; 92(1): 46 - 60. [Abstract] [Full Text] [PDF]

				H. Miedema, M. Vrouenraets, J. Wierenga, D. Gillespie, B. Eisenberg, W. Meijberg, and W. Nonner Ca2+ Selectivity of a Chemically Modified OmpF with Reduced Pore Volume Biophys. J., December 15, 2006; 91(12): 4392 - 4400. [Abstract] [Full Text] [PDF]

				T. W. Allen, O. S. Andersen, and B. Roux Ion Permeation through a Narrow Channel: Using Gramicidin to Ascertain All-Atom Molecular Dynamics Potential of Mean Force Methodology and Biomolecular Force Fields Biophys. J., May 15, 2006; 90(10): 3447 - 3468. [Abstract] [Full Text] [PDF]

				T. Bastug, A. Gray-Weale, S. M. Patra, and S. Kuyucak Role of Protein Flexibility in Ion Permeation: A Case Study in Gramicidin A Biophys. J., April 1, 2006; 90(7): 2285 - 2296. [Abstract] [Full Text] [PDF]

				M. H. Cheng, M. Cascio, and R. D. Coalson Theoretical Studies of the M2 Transmembrane Segment of the Glycine Receptor: Models of the Open Pore Structure and Current-Voltage Characteristics Biophys. J., September 1, 2005; 89(3): 1669 - 1680. [Abstract] [Full Text] [PDF]

				T. W. Allen, O.S. Andersen, and B. Roux On the Importance of Atomic Fluctuations, Protein Flexibility, and Solvent in Ion Permeation J. Gen. Physiol., November 29, 2004; 124(6): 679 - 690. [Abstract] [Full Text] [PDF]

				V. L. Dorman and P. C. Jordan Ionic Permeation Free Energy in Gramicidin: A Semimicroscopic Perspective Biophys. J., June 1, 2004; 86(6): 3529 - 3541. [Abstract] [Full Text] [PDF]

				G. V. Miloshevsky and P. C. Jordan Anion Pathway and Potential Energy Profiles along Curvilinear Bacterial ClC Cl- Pores: Electrostatic Effects of Charged Residues Biophys. J., February 1, 2004; 86(2): 825 - 835. [Abstract] [Full Text] [PDF]