Dataset Information

Poisson-Boltzmann-Nernst-Planck model.

ABSTRACT: The Poisson-Nernst-Planck (PNP) model is based on a mean-field approximation of ion interactions and continuum descriptions of concentration and electrostatic potential. It provides qualitative explanation and increasingly quantitative predictions of experimental measurements for the ion transport problems in many areas such as semiconductor devices, nanofluidic systems, and biological systems, despite many limitations. While the PNP model gives a good prediction of the ion transport phenomenon for chemical, physical, and biological systems, the number of equations to be solved and the number of diffusion coefficient profiles to be determined for the calculation directly depend on the number of ion species in the system, since each ion species corresponds to one Nernst-Planck equation and one position-dependent diffusion coefficient profile. In a complex system with multiple ion species, the PNP can be computationally expensive and parameter demanding, as experimental measurements of diffusion coefficient profiles are generally quite limited for most confined regions such as ion channels, nanostructures and nanopores. We propose an alternative model to reduce number of Nernst-Planck equations to be solved in complex chemical and biological systems with multiple ion species by substituting Nernst-Planck equations with Boltzmann distributions of ion concentrations. As such, we solve the coupled Poisson-Boltzmann and Nernst-Planck (PBNP) equations, instead of the PNP equations. The proposed PBNP equations are derived from a total energy functional by using the variational principle. We design a number of computational techniques, including the Dirichlet to Neumann mapping, the matched interface and boundary, and relaxation based iterative procedure, to ensure efficient solution of the proposed PBNP equations. Two protein molecules, cytochrome c551 and Gramicidin A, are employed to validate the proposed model under a wide range of bulk ion concentrations and external voltages. Extensive numerical experiments show that there is an excellent consistency between the results predicted from the present PBNP model and those obtained from the PNP model in terms of the electrostatic potentials, ion concentration profiles, and current-voltage (I-V) curves. The present PBNP model is further validated by a comparison with experimental measurements of I-V curves under various ion bulk concentrations. Numerical experiments indicate that the proposed PBNP model is more efficient than the original PNP model in terms of simulation time.

SUBMITTER: Zheng Q

PROVIDER: S-EPMC3122111 | biostudies-literature | 2011 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Poisson-Boltzmann-Nernst-Planck model.

Zheng Qiong Q Wei Guo-Wei GW

The Journal of chemical physics 20110501 19

The Poisson-Nernst-Planck (PNP) model is based on a mean-field approximation of ion interactions and continuum descriptions of concentration and electrostatic potential. It provides qualitative explanation and increasingly quantitative predictions of experimental measurements for the ion transport problems in many areas such as semiconductor devices, nanofluidic systems, and biological systems, despite many limitations. While the PNP model gives a good prediction of the ion transport phenomenon ...[more]

PMID: 21599038

Similar Datasets

Project description:The Poisson Nernst-Planck (PNP) theory is a simplified continuum model for a wide variety of chemical, physical and biological applications. Its ability of providing quantitative explanation and increasingly qualitative predictions of experimental measurements has earned itself much recognition in the research community. Numerous computational algorithms have been constructed for the solution of the PNP equations. However, in the realistic ion-channel context, no second order convergent PNP algorithm has ever been reported in the literature, due to many numerical obstacles, including discontinuous coefficients, singular charges, geometric singularities, and nonlinear couplings. The present work introduces a number of numerical algorithms to overcome the abovementioned numerical challenges and constructs the first second-order convergent PNP solver in the ion-channel context. First, a Dirichlet to Neumann mapping (DNM) algorithm is designed to alleviate the charge singularity due to the protein structure. Additionally, the matched interface and boundary (MIB) method is reformulated for solving the PNP equations. The MIB method systematically enforces the interface jump conditions and achieves the second order accuracy in the presence of complex geometry and geometric singularities of molecular surfaces. Moreover, two iterative schemes are utilized to deal with the coupled nonlinear equations. Furthermore, extensive and rigorous numerical validations are carried out over a number of geometries, including a sphere, two proteins and an ion channel, to examine the numerical accuracy and convergence order of the present numerical algorithms. Finally, application is considered to a real transmembrane protein, the Gramicidin A channel protein. The performance of the proposed numerical techniques is tested against a number of factors, including mesh sizes, diffusion coefficient profiles, iterative schemes, ion concentrations, and applied voltages. Numerical predictions are compared with experimental measurements.

Project description:Protegrin peptides are potent antimicrobial agents believed to act against a variety of pathogens by forming nonselective transmembrane pores in the bacterial cell membrane. We have employed 3D Poisson-Nernst-Planck (PNP) calculations to determine the steady-state ion conduction characteristics of such pores at applied voltages in the range of -100 to +100 mV in 0.1 M KCl bath solutions. We have tested a variety of pore structures extracted from molecular dynamics (MD) simulations based on an experimentally proposed octomeric pore structure. The computed single-channel conductance values were in the range of 290-680 pS. Better agreement with the experimental range of 40-360 pS was obtained using structures from the last 40 ns of the MD simulation, where conductance values range from 280 to 430 pS. We observed no significant variation of the conductance with applied voltage in any of the structures that we tested, suggesting that the voltage dependence observed experimentally is a result of voltage-dependent channel formation rather than an inherent feature of the open pore structure. We have found the pore to be highly selective for anions, with anionic to cationic current ratios (I(Cl-)/I(K+)) on the order of 10(3). This is consistent with the highly cationic nature of the pore but surprisingly in disagreement with the experimental finding of only slight anionic selectivity. We have additionally tested the sensitivity of our PNP model to several parameters and found the ion diffusion coefficients to have a significant influence on conductance characteristics. The best agreement with experimental data was obtained using a diffusion coefficient for each ion set to 10% of the bulk literature value everywhere inside the channel, a scaling used by several other studies employing PNP calculations. Overall, this work presents a useful link between previous work focused on the structure of protegrin pores and experimental efforts aimed at investigating their conductance characteristics.

Dataset Information

Poisson-Boltzmann-Nernst-Planck model.

Publications

Poisson-Boltzmann-Nernst-Planck model.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets