TY - JOUR
T1 - Kinetic lattice grand canonical Monte Carlo simulation for ion current calculations in a model ion channel system
AU - Hwang, Hyonseok
AU - Schatz, George C.
AU - Ratner, Mark A.
N1 - Funding Information:
This work was supported by the Network for Computational Nanotechnology (NCN) through a grant from the National Science Foundation, and by the Northwestern Materials Research Center (NSF Grant No. DMR-0520513).
PY - 2007
Y1 - 2007
N2 - 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.
AB - 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.
UR - http://www.scopus.com/inward/record.url?scp=34547161294&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=34547161294&partnerID=8YFLogxK
U2 - 10.1063/1.2748373
DO - 10.1063/1.2748373
M3 - Article
C2 - 17640144
AN - SCOPUS:34547161294
VL - 127
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
SN - 0021-9606
IS - 2
M1 - 024706
ER -