Monte Carlo Simulation of Electrolyte Solutions in Biology

Past and present of electrochemical science in Hungary

The electrochemistry-related scientific activities in Hungary over the past 3 decades are reviewed. In the first section, we summarize those research areas that are already ceased; in the next section, the ongoing research is discussed; finally, the trends and outlook are highlighted. A special emphasis is put on new experimental methods elaborated in the country.

Individual Ion Activity Coefficients in Aqueous Electrolytes from Explicit-Water Molecular Dynamics Simulations

We compute individual ion activity coefficients (IIACs) in aqueous NaCl, KCl, NaF, and KF solutions from explicit-water molecular dynamics simulations. Free energy changes are obtained from insertion of single ions-accompanied by uniform neutralizing backgrounds-into solution by gradually turning on first Lennard-Jones interactions, followed by Coulombic interactions using Ewald electrostatics. Simulations are performed at multiple system sizes, and all results are extrapolated to the thermodynamic limit, thus eliminating any possible artifacts from the neutralizing backgrounds. Because of controversies associated with measurements of IIACs from electrochemical cells with ion-selective electrodes, the reported experimental data are not widely accepted; thus there remains a knowledge gap with respect to the contributions of individual ions to solution nonidealities. Our results are in good qualitative agreement with these reported measurements, though significantly larger in magnitude. In particular, the relative positioning for the activity coefficients of anions and cations matches the experimental ordering for all four systems. This work establishes a robust thermodynamic framework, without a need to invoke extra hypotheses, that sheds light on the behavior of individual ions and their contributions to nonidealities of aqueous electrolyte solutions.

Modeling the Device Behavior of Biological and Synthetic Nanopores with Reduced Models

Biological ion channels and synthetic nanopores are responsible for passive transport of ions through a membrane between two compartments. Modeling these ionic currents is especially amenable to reduced models because the device functions of these pores, the relation of input parameters (e.g., applied voltage, bath concentrations) and output parameters (e.g., current, rectification, selectivity), are well defined. Reduced models focus on the physics that produces the device functions (i.e., the physics of how inputs become outputs) rather than the atomic/molecular-scale physics inside the pore. Here, we propose four rules of thumb for constructing good reduced models of ion channels and nanopores. They are about (1) the importance of the axial concentration profiles, (2) the importance of the pore charges, (3) choosing the right explicit degrees of freedom, and (4) creating the proper response functions. We provide examples for how each rule of thumb helps in creating a reduced model of device behavior.

Multiscale analysis of the effect of surface charge pattern on a nanopore's rectification and selectivity properties: From all-atom model to Poisson-Nernst-Planck

We report a multiscale modeling study for charged cylindrical nanopores using three modeling levels that include (1) an all-atom explicit-water model studied with molecular dynamics, and reduced models with implicit water containing (2) hard-sphere ions studied with the Local Equilibrium Monte Carlo simulation method (computing ionic correlations accurately), and (3) point ions studied with Poisson-Nernst-Planck theory (mean-field approximation). We show that reduced models are able to reproduce device functions (rectification and selectivity) for a wide variety of charge patterns, that is, reduced models are useful in understanding the mesoscale physics of the device (i.e., how the current is produced). We also analyze the relationship of the reduced implicit-water models with the explicit-water model and show that diffusion coefficients in the reduced models can be used as adjustable parameters with which the results of the explicit- and implicit-water models can be related. We find that the values of the diffusion coefficients are sensitive to the net charge of the pore but are relatively transferable to different voltages and charge patterns with the same total charge.

Do Bistable Steric Poisson-Nernst-Planck Models Describe Single-Channel Gating?

Experiments measuring currents through single protein channels show unstable currents, a phenomena called the gating of a single channel. Channels switch between an "open" state with a well-defined single amplitude of current and "closed" states with nearly zero current. The existing mean-field theory of ion channels focuses almost solely on the open state. The physical modeling of the dynamical features of ion channels is still in its infancy and does not describe the transitions between open and closed states nor the distribution of the duration times of open states. One hypothesis is that gating corresponds to noise-induced fast transitions between multiple steady (equilibrium) states of the underlying system. In this work, we aim to test this hypothesis. Particularly, our study focuses on the (high-order) steric Poisson-Nernst-Planck (PNP)-Cahn-Hilliard model since it has been successful in predicting permeability and selectivity of ionic channels in their open state and since it gives rise to multiple steady states. We show that this system gives rise to a gatinglike behavior, but that important features of this switching behavior are different from the defining features of gating in biological systems. Furthermore, we show that noise prohibits switching in the system of study. The above phenomena are expected to occur in other PNP-type models, strongly suggesting that one has to go beyond overdamped (gradient flow) Nernst-Planck type dynamics to explain the spontaneous gating of single channels.

Simulation of a model nanopore sensor: Ion competition underlies device behavior

We study a model nanopore sensor with which a very low concentration of analyte molecules can be detected on the basis of the selective binding of the analyte molecules to the binding sites on the pore wall. The bound analyte ions partially replace the current-carrier cations in a thermodynamic competition. This competition depends both on the properties of the nanopore and the concentrations of the competing ions (through their chemical potentials). The output signal given by the device is the current reduction caused by the presence of the analyte ions. The concentration of the analyte ions can be determined through calibration curves. We model the binding site with the square-well potential and the electrolyte as charged hard spheres in an implicit background solvent. We study the system with a hybrid method in which we compute the ion flux with the Nernst-Planck (NP) equation coupled with the Local Equilibrium Monte Carlo (LEMC) simulation technique. The resulting NP+LEMC method is able to handle both strong ionic correlations inside the pore (including finite size of ions) and bulk concentrations as low as micromolar. We analyze the effect of bulk ion concentrations, pore parameters, binding site parameters, electrolyte properties, and voltage on the behavior of the device.

Incorporating Born solvation energy into the three-dimensional Poisson-Nernst-Planck model to study ion selectivity in KcsA K^{+} channels

Potassium channels are much more permeable to potassium than sodium ions, although potassium ions are larger and both carry the same positive charge. This puzzle cannot be solved based on the traditional Poisson-Nernst-Planck (PNP) theory of electrodiffusion because the PNP model treats all ions as point charges, does not incorporate ion size information, and therefore cannot discriminate potassium from sodium ions. The PNP model can qualitatively capture some macroscopic properties of certain channel systems such as current-voltage characteristics, conductance rectification, and inverse membrane potential. However, the traditional PNP model is a continuum mean-field model and has no or underestimates the discrete ion effects, in particular the ion solvation or self-energy (which can be described by Born model). It is known that the dehydration effect (closely related to ion size) is crucial to selective permeation in potassium channels. Therefore, we incorporated Born solvation energy into the PNP model to account for ion hydration and dehydration effects when passing through inhomogeneous dielectric channel environments. A variational approach was adopted to derive a Born-energy-modified PNP (BPNP) model. The model was applied to study a cylindrical nanopore and a realistic KcsA channel, and three-dimensional finite element simulations were performed. The BPNP model can distinguish different ion species by ion radius and predict selectivity for K^{+} over Na^{+} in KcsA channels. Furthermore, ion current rectification in the KcsA channel was observed by both the PNP and BPNP models. The I-V curve of the BPNP model for the KcsA channel indicated an inward rectifier effect for K^{+} (rectification ratio of ∼3/2) but indicated an outward rectifier effect for Na^{+} (rectification ratio of ∼1/6).

Nonlocal Poisson-Fermi model for ionic solvent

We propose a nonlocal Poisson-Fermi model for ionic solvent that includes ion size effects and polarization correlations among water molecules in the calculation of electrostatic potential. It includes the previous Poisson-Fermi models as special cases, and its solution is the convolution of a solution of the corresponding nonlocal Poisson dielectric model with a Yukawa-like kernel function. The Fermi distribution is shown to be a set of optimal ionic concentration functions in the sense of minimizing an electrostatic potential free energy. Numerical results are reported to show the difference between a Poisson-Fermi solution and a corresponding Poisson solution.