Kinetic lattice grand canonical Monte Carlo simulation for ion current calculations in a model ion channel system

Research output: Contribution to journalArticle

11 Citations (Scopus)

Abstract

An algorithm in which kinetic lattice grand canonical Monte Carlo simulations are combined with mean field theory (KLGCMC/MF) is presented to calculate ion currents in a model ion channel system. In this simulation, the relevant region of the system is treated by KLGCMC simulations, while the rest of the system is described by modified Poisson-Boltzmann mean field theory. Calculation of reaction field due to induced charges on the channel/water and membrane/water boundaries is carried out using a basis-set expansion method [Im and Roux, J. Chem. Phys. 115, 4850 (2001)]. Calculation of ion currents, electrostatic potentials, and ion concentrations, as obtained from the KLGCMC/MF simulations, shows good agreement with Poisson-Nernst-Planck (PNP) theory predictions when the channel and membrane have the same dielectric constant as water. If the channel and membrane have a lower dielectric constant than water, however, there is a considerable difference between the KLGCMC/MF and PNP predictions. This difference is attributed to the reaction field, which is missing in PNP theory. It is demonstrated that the reaction field as well as fixed charges in the channel play key roles in selective ion transport. Limitations and further development of the current KLGCMC/MF approach are also discussed.

Original languageEnglish
Article number024706
JournalJournal of Chemical Physics
Volume127
Issue number2
DOIs
Publication statusPublished - Aug 1 2007

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Fingerprint Dive into the research topics of 'Kinetic lattice grand canonical Monte Carlo simulation for ion current calculations in a model ion channel system'. Together they form a unique fingerprint.

  • Cite this