Incorporation of inhomogeneous ion diffusion coefficients into kinetic lattice grand canonical monte carlo simulations and application to ion current calculations in a simple model ion channel

To deal with inhomogeneous diffusion coefficients of ions without altering the lattice spacing in the kinetic lattice grand canonical Monte Carlo (KLGCMC) simulation, an algorithm that incorporates diffusion coefficient variation into move probabilities is proposed and implemented into KLGCMC calculations. Using this algorithm, the KLGCMC simulation method is applied to the calculation of ion currents in a simple model ion channel system. Comparisons of ion currents and ion concentrations from these simulations with Poisson-Nernst-Planck (PNP) results show good agreement between the two methods for parameters where the latter method is expected to be accurate.     

Role of the intracellular cavity in potassium channel conductivity

The role of several fragments of the potassium channel KcsA has been examined by the Poisson-Nernst-Planck (PNP) theory. The efficiency of the computational method allowed comparing a large number of channel models, with different intracellular gate openings, partial atomic charges, and amino acid sequences. Perhaps counter-intuitively, the calculated ion current decreases when the mean radius of the entrance cavity increases. Widening of the vestibule, in fact, increases the volume accessible to water, which is the volume with a high dielectric constant. In turn, water screens the attractive charges of the P-loop backbone. The backbone charges of the M2 helixes instead oppose the entrance of potassium ions through a complicated mechanism that can be separated in the activity of two interfering dipoles. The conductance of the KcsA models increased when two neutral residues in M2 were mutated to glutamic acid, in agreement with experimental results (Brelidze, T. I.; Niu, X.; Magleby, K. L. PNAS 2003, 100, 9017-9022). As a general conclusion, a relation between channel conductance and potassium concentration in the intracellular cavity emerged. Although the ion transport is the result of the fine balance of a number of different effects, the experimental results can be reproduced quantitatively only on the basis of electrostatic forces, which are the only driving forces modeled by the PNP theory.     

Ion current calculations based on three dimensional Poisson-Nernst-Planck theory for a cyclic peptide nanotube

Ion current calculations based on Poisson-Nernst-Planck (PNP) theory are performed for a synthetic cyclic peptide nanotube that consists of eight or ten cyclo[(-L-Trp-D-Leu-)4] embedded in a lipid bilayer membrane to investigate the ion transport properties of the nanotube. To explore systems with arbitrary geometries, three-dimensional PNP theory is implemented using a finite difference method. The influence of dipolar lipid molecules on the ion currents is also examined by turning on or off the charges of the lipid dipoles in dipalmitoylphosphatidylcholine (DPPC). Comparisons between the calculated and experimentally measured ion currents show that the PNP approach agrees well with the measurements at low ion concentrations but overestimates the currents at higher concentrations. Concentration profiles reveal the selectivity of the peptide nanotube to cations, which is attributed to the negatively charged carbonyl oxygens inside the nanotube. The dominant cation and the minimum anion concentrations inside the cyclic peptide nanotube suggest that these cyclic peptide nanotubes can be employed as ion sensors. In the case of the polar DPPC bilayer, smaller currents are obtained in the calculation. The variation of current with polarity of the lipids implies that both polar and nonpolar lipid bilayer membranes can be utilized to regulate ion currents in the peptide nanotube and other ion channels. Strengths and limitations of the PNP theory are also discussed.     

A parallel finite element simulator for ion transport through three‐dimensional ion channel systems

A parallel finite element simulator, ichannel, is developed for ion transport through three‐dimensional ion channel systems that consist of protein and membrane. The coordinates of heavy atoms of the protein are taken from the Protein Data Bank and the membrane is represented as a slab. The simulator contains two components: a parallel adaptive finite element solver for a set of Poisson–Nernst–Planck (PNP) equations that describe the electrodiffusion process of ion transport, and a mesh generation tool chain for ion channel systems, which is an essential component for the finite element computations. The finite element method has advantages in modeling irregular geometries and complex boundary conditions. We have built a tool chain to get the surface and volume mesh for ion channel systems, which consists of a set of mesh generation tools. The adaptive finite element solver in our simulator is implemented using the parallel adaptive finite element package Parallel Hierarchical Grid (PHG) developed by one of the authors, which provides the capability of doing large scale parallel computations with high parallel efficiency and the flexibility of choosing high order elements to achieve high order accuracy. The simulator is applied to a real transmembrane protein, the gramicidin A (gA) channel protein, to calculate the electrostatic potential, ion concentrations and I – V curve, with which both primitive and transformed PNP equations are studied and their numerical performances are compared. To further validate the method, we also apply the simulator to two other ion channel systems, the voltage dependent anion channel (VDAC) and α‐Hemolysin (α‐HL). The simulation results agree well with Brownian dynamics (BD) simulation results and experimental results. Moreover, because ionic finite size effects can be included in PNP model now, we also perform simulations using a size‐modified PNP (SMPNP) model on VDAC and α‐HL. It is shown that the size effects in SMPNP can effectively lead to reduced current in the channel, and the results are closer to BD simulation results. © 2013 Wiley Periodicals, Inc.     

Electric-field-controlled water and ion permeation of a hydrophobic nanopore

The permeation of hydrophobic, cylindrical nanopores by water molecules and ions is investigated under equilibrium and out-of-equilibrium conditions by extensive molecular-dynamics simulations. Neglecting the chemical structure of the confining pore surface, we focus on the effects of pore radius and electric field on permeation. The simulations confirm the intermittent filling of the pore by water, reported earlier under equilibrium conditions for pore radii larger than a critical radius R(c). Below this radius, water can still permeate the pore under the action of a strong electric field generated by an ion concentration imbalance at both ends of the pore embedded in a structureless membrane. The water driven into the channel undergoes considerable electrostriction characterized by a mean density up to twice the bulk density and by a dramatic drop in dielectric permittivity which can be traced back to a considerable distortion of the hydrogen-bond network inside the pore. The free-energy barrier to ion permeation is estimated by a variant of umbrella sampling for Na(+), K(+), Ca(2+), and Cl(-) ions, and correlates well with known solvation free energies in bulk water. Starting from an initial imbalance in ion concentration, equilibrium is gradually restored by successive ion passages through the water-filled pore. At each passage the electric field across the pore drops, reducing the initial electrostriction, until the pore, of radius less than R(c), closes to water and hence to ion transport, thus providing a possible mechanism for voltage-dependent gating of hydrophobic pores.     

Electrokinetics in nanochannels: part II. Mobility dependence on ion density and ionic current measurements

In the first of this two-paper series, a new model was developed for calculating the electric potential field in a long, thin nanochannel with overlapped electric double layers. The model takes into account the dependence of ion mobility on local ion densities and pH. This model is used here to study and demonstrate the effect of ion density and pH on ionic current measurements. A comparison is shown of predictions based on each of three boundary conditions, as studied in Part I. The model developed in Part I is validated by comparing simulations with measurements of ionic current as a function of sodium borate concentration. Results show that predictions based on extended Debye-Hückel theory for ion mobility significantly improve the accuracy of simulations, but that these do not predict exact scaling behavior. A simple bulk conductivity measurement used as input parameter for the simulations, in place of the predicted bulk conductivity (K(0)), guarantees agreement with data in the thin EDL region. Results also indicate that the charge regulation boundary condition, complemented with an adequate bulk electrolyte model, provides better agreement with experimental trends than the specified zeta potential or specified surface net charge boundary conditions. Further, it is shown that currents due to advection (by electroosmotic flow) are in all cases studied less than 25% of the total current in the system.     

Electrostatic contribution to the ion solvation energy: cavity effects

Ion solvation process has been analysed for the spherically symmetrical system where an ion is located inside a cavity surrounded by an isotropic nonlocal dielectric medium. It has been proven that for any dielectric properties of the medium, the electric field outside the cavity as well as the ion solvation energy depend only on the total ion charge but not of the particular distribution of the ion charge density inside the cavity. These characteristics remain unchanged if the charge is displaced from the external boundary of the cavity into it. Analytical formulas for them have been derived for a particular model of the nonlocal dielectric function. Comparison of results for the solvation energy on the basis of this new theory and of the conventional approach (disregarding the existence of the cavity) shows a significant difference between their predictions if the ion charge is displaced inside the ion cavity.     

The influence of discrete surface charges on the force between charged surfaces

The force between two parallel charged flat surfaces, with discrete surface charges, has been calculated with Monte Carlo simulations for different values of the electrostatic coupling. For low electrostatic coupling (small counterion valence, small surface charge, high dielectric constant, and high temperature) the total force is dominated by the entropic contribution and can be described by mean field theory, independent of the character of the surface charges. For moderate electrostatic coupling, counterion correlation effects lead to a smaller repulsion than predicted by mean field theory. This correlation effect is strengthened by discrete surface charges and the repulsive force is further reduced. For large electrostatic coupling the total force for smeared out surface charges is known to be attractive due to counterion correlations. If discrete surface charges are considered the attractive force is weakened and can even be turned into a repulsive force. This is due to the counterions being strongly correlated to the discrete surface charges forming effective, oppositely directed, dipoles on the two walls.     

Membrane potential across low-water-content charged membranes: Effect of ion pairing

In the present paper, we systematically examined the ion-pairing effect in low-water-content charged membranes. Cation- and anion-exchange membranes with various water contents and homogeneous fixed-charge distribution were prepared by radical copolymerization and then characterized by membrane potential measurements. The experimental results were analyzed by our recently developed theoretical model (Yamamoto, R.; Matsumoto, H.; Tanioka, A. J. Phys. Chem. B 2003, 107, 10615), which is based on the Donnan equilibrium, the Nernst-Planck equation for ion flux, and the Fuoss formalism for ion-pair formation between the fixed-charge group and the counterion in the membrane. The theoretical predictions agreed well with the experimental results for both cation- and anion-exchange membranes. This supported the belief that the ion-pairing effect was substantial in a low-water-content membrane system. Our theoretical analysis also showed the following results: (i) the dielectric constant in the membrane, epsilon(r), was smaller than the value in bulk water, (ii) the center-to-center distance of the ion pair, a, was independent of the water content of the membranes, and (iii) the charge effectiveness of all membranes, Q, was small (<0.35).     

Theory of space-charge polarization for determining ionic constants of electrolytic solutions

A theoretical expression of the complex dielectric constant attributed to space-charge polarization has been derived under an electric field calculated using Poisson's equation considering the effects of bound charges on ions. The frequency dependence of the complex dielectric constant of chlorobenzene solutions doped with tetrabutylammonium tetraphenylborate (TBATPB) has been analyzed using the theoretical expression, and the impact of the bound charges on the complex dielectric constant has been clarified quantitatively in comparison with a theory that does not consider the effect of the bound charges. The Stokes radius of TBA+(=TPB-) determined by the present theory shows a good agreement with that determined by conductometry in the past; hence, the present theory should be applicable to the direct determination of the mobility of ion species in an electrolytic solution without the need to measure ionic limiting equivalent conductance and transport number.     

反渗透膜UTC-70在水溶液中的介电谱及其解析

根据平面层状体系介电弛豫理论研究了反渗透膜UTC-70在各种浓度氯化钠和氯化钾溶液中的介电弛豫行为.利用计算机拟合的方法得到膜/溶液体系的介电参数,并由此计算得到了UTC-70膜相和水溶液相的相参数,获得了反映反渗透膜UTC-70荷电情况的信息及其与电解质溶液浓度的关系,介电解析的结果解释了介电弛豫的产生机制.     

Hierarchical multiscale simulation of electrokinetic transport in silica nanochannels at the point of zero charge

Effects of nanoscale confinement and partial charges that stem from quantum calculations are investigated in silica slit channels filled with 1 M KCl at the point of zero charge by using a hierarchical multiscale simulation methodology. Partial charges of both bulk and surface atoms from ab initio quantum calculations that take into account bond polarization and electronegativity are used in molecular dynamics (MD) simulations to obtain ion and water concentration profiles for channel widths of 1.1, 2.1, 2.75, and 4.1 nm. The interfacial electron density profiles of simulations matched well with that of recent X-ray reflectivity experiments. By simulating corresponding channels with no partial charges, it was observed that the partial charges affect the concentration profiles and transport properties such as diffusion coefficients and mobilities up to a distance of about 3 sigma(O)(-)(O) from the surface. Both in uncharged and partially charged cases, oscillations in concentration profiles of K(+) and Cl(-) ions give rise to an electro-osmotic flow in the presence of an external electric field, indicating the presence of an electric double layer at net zero surface charge, contrary to the expectations from classical continuum theory. I-V curves in a channel-bath system using ionic mobilities from MD simulations were significantly different for channels with and without partial charges for channel widths less than 4.1 nm.     

Incidence of partial charges on ion selectivity in potassium channels

Potassium channels are membrane proteins known to select potassium over sodium ions at a high diffusion rate. We conducted ab initio calculations on a filter model of KcsA of about 300 atoms at the Hartree-Fock level of theory. Partial charges were derived from the quantum mechanically determined electrostatic potential either with Merz-Kollman or Hinsen-Roux schemes. Large polarization and/or charge transfer occur on potassium ions located in the filter, while the charges on sodium ions remain closer to unity. As a result, a weaker binding is obtained for K(+) ions. Using a simplified version of a permeation model based on the concerted-motion mechanism for ion translocation within the single-file ion channel [P. H. Nelson, J. Chem. Phys. 117, 11396 (2002)], we discuss how differences in polarization effects in the adducts with K(+) and Na(+) can play a role as for ionic selectivity and conductance.     

Tangential streaming potential as a tool in modeling of ion transport through nanoporous membranes

Tangential streaming potential (TSP) measurements have been carried out so as to assess the electrokinetic properties of the active layer of organic nanofiltration (NF) membranes. Due to the porous structure of NF membranes, cares must be taken to convert the experimental data into zeta potential. Indeed, an assumption that is implicitly made in Smoluchowski's theory (or in related approaches accounting for the surface conduction phenomenon) is that both streaming and conduction currents involved in the streaming potential process flow through an identical path. Such an assumption does not hold with porous membranes since the conduction current is expected to flow wherever the electric conductivity differs from zero. Consequently, a non-negligible share of the conduction current is likely to flow through the membrane body filled with the electrolyte solution. This phenomenon has been taken into account by carrying out a series of TSP measurements at various channel heights. Experiments have been conducted with various electrolyte solutions. The inferred zeta potentials have been further converted into membrane volume charge densities which have been used to predict the membrane performances in terms of rejection rates. The conventional NF theory, i.e. based on a steric/Donnan exclusion mechanism, has been found to be unable to describe the experimental rejection rates. Using the volume charge density of the membrane as an adjustable parameter, it has been shown that the conventional theory even predicts the opposite sign for the membrane charge. On the other hand, the experimental rejection rates have been well described by including dielectric effects in the exclusion mechanism. In this case, a noticeable lowering of the effective dielectric constant of the electrolyte solution inside pores has been predicted (with respect to the bulk value).     

Charge density identification in ion channels

Biological channels permeate ions through cell membranes. Ion channels carry a permanent charge that has a significant role in determining channel's permeation properties such as selectivity to certain ions, current amplitude, etc. In this paper we deal with the determination of the permanent charge from current-voltage curves. The ion channel current behavior is modelled by Poisson-Nernst-Planck (PNP) equations system. Previous works on the fixed charge density identification problem contain several ill-posed steps and linearization of the nonlinear PNP system. We suggest here several methods to make these algorithms more stable and accurate.     

Effects of ion atmosphere relaxation on dipole isomerization reactions in electrolyte solutions

A theoretical study of the effects of ion atmosphere relaxation on the rate of a model dipole isomerization reaction in electrolyte solutions is presented. The time-dependent ion atmosphere friction is calculated by using a molecular hydrodynamic theory which properly includes the static and dynamic interionic correlations through ionic structure factors and van Hove functions. The rate constant is determined by employing the well-known Grote-Hynes theory. Numerical results are obtained for the time-dependent ion atmosphere friction and for the rate of isomerization reaction in electrolyte solutions of varying ion concentration and dielectric constant. It is found that the ion atmosphere friction can have significant effects in reducing the rate of isomerization below the prediction of equilibrium solvation transition state theory.     

Continuum treatment of electronic polarization effect

A continuum treatment of electronic polarization has been explored for in molecular mechanics simulations in implicit solvents. The dielectric constant for molecule interior is the only parameter in the continuum polarizable model. A value of 4 is found to yield optimal agreement with high-level ab initio quantum mechanical calculations for the tested molecular systems. Interestingly, its performance is not sensitive to the definition of molecular volume, in which the continuum electronic polarization is defined. In this model, quantum mechanical electrostatic field in different dielectric environments from vacuum, low-dielectric organic solvent, and water can be used simultaneously in atomic charge fitting to achieve consistent treatment of electrostatic interactions. The tests show that a single set of atomic charges can be used consistently in different dielectric environments and different molecular conformations, and the atomic charges transfer well from training monomers to tested dimers. The preliminary study gives us the hope of developing a continuum polarizable force field for more consistent simulations of proteins and nucleic acids in implicit solvents.