Modified Donnan potentials for ion transport through biological ion channels

Cellular communication among smooth muscle cells: The role of membrane potential via connexins

Communication via action potentials among neurons has been extensively studied. However, effective communication without action potentials is ubiquitous in biological systems, yet it has received much less attention in comparison. Multi-cellular communication among smooth muscles is crucial for regulating blood flow, for example. Understanding the mechanism of this non-action potential communication is critical in many cases, like synchronization of cellular activity, under normal and pathological conditions. In this paper, we employ a multi-scale asymptotic method to derive a macroscopic homogenized bidomain model from the microscopic electro-neutral (EN) model. This is achieved by considering different diffusion coefficients and incorporating nonlinear interface conditions. Subsequently, the homogenized macroscopic model is used to investigate communication in multi-cellular tissues. Our computational simulations reveal that the membrane potential of syncytia, formed by interconnected cells via connexins, plays a crucial role in propagating oscillations from one region to another, providing an effective means for fast cellular communication. Statement of Significance: In this study, we investigated cellular communication and ion transport in vascular smooth muscle cells, shedding light on their mechanisms under normal and abnormal conditions. Our research highlights the potential of mathematical models in understanding complex biological systems. We developed effective macroscale electro-neutral bi-domain ion transport models and examined their behavior in response to different stimuli. Our findings revealed the crucial role of connexinmediated membrane potential changes and demonstrated the effectiveness of cellular communication through syncytium membranes. Despite some limitations, our study provides valuable insights into these processes and emphasizes the importance of mathematical modeling in unraveling the complexities of cellular communication and ion transport.

Salt Effect on Donnan Equilibrium in Montmorillonite Demonstrated with Molecular Dynamics Simulations

Donnan equilibrium governs the distribution of ions in many systems such as ion exchange membranes and biological cells in contact with an external electrolyte. Herein, Donnan equilibrium between bulk salt solution and bihydrated montmorillonite was investigated because such a system is of great importance for many nuclear waste disposal concepts. Specifically, we used molecular dynamics simulations to determine the partition coefficient of chloride, which was achieved by calculating the free-energy difference of chloride in the interlayer and the bulk using enhanced sampling methodology. Montmorillonite in equilibrium with either NaCl or CaCl₂ was examined to elucidate the general difference between 1:1 and 2:1 salts. The concentration dependence of the partition coefficient for each salt was determined using three and four concentrations for NaCl and CaCl₂, respectively. In the case of NaCl, we found that the partition coefficient increased linearly with the concentration, while for CaCl₂, the increase was proportional to the square root of the concentration. A derivation of the partition coefficient using general Donnan theory that includes excess free energy contributions beyond the electrostatic Donnan potential is also presented. For both salts, the agreement between the partition coefficient from the simulations and Donnan theory was excellent. Although Donnan theory is a continuum theory derived without any reference to atomistic details, the present results justify its application to systems with nanoscale pores.

Bounds for systems of coupled advection-diffusion equations with application to the Poisson-Nernst-Plank equations

We consider coupled systems of advection-diffusion equations with initial and boundary conditions and determine conditions on the advection terms that allow us to obtain solutions that can be explicitly bounded above and below using the initial and boundary conditions. Given the advection terms, using our methodology one can easily check if such bounds can be obtained and then one can construct the necessary nonlinear transformation to allow the bounds to be determined. We apply this technique to determine bounding quantities for a number of examples. In particular, we show that the three-ion electroneutral Poisson-Nernst-Planck system of equations can be transformed into a system, which allows for the use of our techniques and we determine the bounding quantities. In addition, we determine the general form of advection terms that allow these techniques to be applied and show that our method can be applied to a very wide class of advection-diffusion equations.

Analyzing the effect of ion binding to the membrane-surface on regulating the light-induced transthylakoid electric potential (ΔΨm)

The transthylakoid membrane potential (ΔΨm) is essential because it can drive the ATP synthesis through the CF0-CF1 type of ATP-synthase in chloroplasts as an energetic equivalent similar to ΔpH. In addition, a high fraction of proton motive force (PMF) stored as the ΔΨm component is physiologically important in the acclimation of photosynthesis to environmental stresses. It has been shown that ΔΨm is the sum of the Donnan potential difference (ΔΨdn) and the diffusion potential difference (ΔΨd). Specifically, ΔΨdn, ΔΨd, and ΔΨm are strongly associated with the ionic activities near the membrane surface, particularly, the extent of ion binding to the charged/neutral sites adjacent to the membrane surface. However, an in-depth analysis of the effect of altered cationic binding to the membrane surface on adjusting the transthylakoid electric potentials (ΔΨdn, ΔΨd, and ΔΨm) is still missing. This lack of a mechanistic understanding is due to the experimental difficulty of closely observing cations binding to the membrane surface in vivo. In this work, a computer model was proposed to investigate the transthylakoid electric phenomena in the chloroplast focusing on the interaction between cations and the negative charges close to the membrane surface. By employing the model, we simulated the membrane potential and consequently, the measured ECS traces, proxing the ΔΨm, were well described by the computing results on continuous illumination followed by a dark-adapted period. Moreover, the computing data clarified the components of transthylakoid membrane potential, unraveled the functional consequences of altered cationic attachment to the membrane surface on adjusting the transthylakoid electric potential, and further revealed the key role played by Donnan potential in regulating the energization of the thylakoid membrane. The current model for calculating electric potentials can function as a preliminary network for the further development into a more detailed theoretical model by which multiple important variables involved in photosynthesis can be explored.

Selectivity of the KcsA potassium channel: Analysis and computation

Ion channels regulate the flux of ions through cell membranes and play significant roles in many physiological functions. Most of the existing literature focuses on computational approaches based on molecular dynamics simulation or numerical solution of the modified Poisson-Nernst-Planck (PNP) system. In this paper, we present an analytical and computational study of a mathematical model of the KcsA potassium channel, including the effects of ion size (Bikerman model) and solvation energy (Born model). Under equilibrium conditions, we obtain an analytical solution of our modified PNP system, which is used to explain selectivity of KcsA of various ions (K^{+}, Na^{+}, Cl^{-}, Ca^{2+}, and Ba^{2+}) due to negative permanent charges inside the filter region and the effect of ion sizes. Our results show that K^{+} is always selected over Na^{+}, as smaller Na^{+} ions have larger solvation energy. As the amount of negative charges in the filter exceeds a critical value, divalent ions (Ca^{2+} and Ba^{2+}) can enter the filter region and block the KcsA channel. For the nonequilibrium cases, due to difficulties associated with a pure analytical or numerical approach, we use a hybrid analytical-numerical method to solve the modified PNP system. Our predictions of selectivity of KcsA channels and saturation phenomenon of the current-voltage (I-V) curve agree with experimental observations.

Electroneutral models for dynamic Poisson-Nernst-Planck systems

The Poisson-Nernst-Planck (PNP) system is a standard model for describing ion transport. In many applications, e.g., ions in biological tissues, the presence of thin boundary layers poses both modeling and computational challenges. In this paper, we derive simplified electroneutral (EN) models where the thin boundary layers are replaced by effective boundary conditions. There are two major advantages of EN models. First, it is much cheaper to solve them numerically. Second, EN models are easier to deal with compared to the original PNP system; therefore, it would also be easier to derive macroscopic models for cellular structures using EN models. Even though the approach used here is applicable to higher-dimensional cases, this paper mainly focuses on the one-dimensional system, including the general multi-ion case. Using systematic asymptotic analysis, we derive a variety of effective boundary conditions directly applicable to the EN system for the bulk region. This EN system can be solved directly and efficiently without computing the solution in the boundary layer. The derivation is based on matched asymptotics, and the key idea is to bring back higher-order contributions into the effective boundary conditions. For Dirichlet boundary conditions, the higher-order terms can be neglected and the classical results (continuity of electrochemical potential) are recovered. For flux boundary conditions, higher-order terms account for the accumulation of ions in boundary layer and neglecting them leads to physically incorrect solutions. To validate the EN model, numerical computations are carried out for several examples. Our results show that solving the EN model is much more efficient than the original PNP system. Implemented with the Hodgkin-Huxley model, the computational time for solving the EN model is significantly reduced without sacrificing the accuracy of the solution due to the fact that it allows for relatively large mesh and time-step sizes.

Structure of electric double layers in capacitive systems and to what extent (classical) density functional theory describes it

Ongoing scientific interest is aimed at the properties and structure of electric double layers (EDLs), which are crucial for capacitive energy storage, water treatment, and energy harvesting technologies like supercapacitors, desalination devices, blue engines, and thermocapacitive heat-to-current converters. A promising tool to describe their physics on a microscopic level is (classical) density functional theory (DFT), which can be applied in order to analyze pair correlations and charge ordering in the primitive model of charged hard spheres. This simple model captures the main properties of ionic liquids and solutions and it predicts many of the phenomena that occur in EDLs. The latter often lead to anomalous response in the differential capacitance of EDLs. This work constructively reviews the powerful theoretical framework of DFT and its recent developments regarding the description of EDLs. It explains to what extent current approaches in DFT describe structural ordering and in-plane transitions in EDLs, which occur when the corresponding electrodes are charged. Further, the review briefly summarizes the history of modeling EDLs, presents applications, and points out limitations and strengths in present theoretical approaches. It concludes that DFT as a sophisticated microscopic theory for ionic systems is expecting a challenging but promising future in both fundamental research and applications in supercapacitive technologies.

Singular perturbation solutions of steady-state Poisson-Nernst-Planck systems

We study the Poisson-Nernst-Planck (PNP) system with an arbitrary number of ion species with arbitrary valences in the absence of fixed charges. Assuming point charges and that the Debye length is small relative to the domain size, we derive an asymptotic formula for the steady-state solution by matching outer and boundary layer solutions. The case of two ionic species has been extensively studied, the uniqueness of the solution has been proved, and an explicit expression for the solution has been obtained. However, the case of three or more ions has received significantly less attention. Previous work has indicated that the solution may be nonunique and that even obtaining numerical solutions is a difficult task since one must solve complicated systems of nonlinear equations. By adopting a methodology that preserves the symmetries of the PNP system, we show that determining the outer solution effectively reduces to solving a single scalar transcendental equation. Due to the simple form of the transcendental equation, it can be solved numerically in a straightforward manner. Our methodology thus provides a standard procedure for solving the PNP system and we illustrate this by solving some practical examples. Despite the fact that for three ions, previous studies have indicated that multiple solutions may exist, we show that all except for one of these solutions are unphysical and thereby prove the existence and uniqueness for the three-ion case.

Critical assessment of OmpF channel selectivity: merging information from different experimental protocols

The ion selectivity of a channel can be quantified in several ways by using different experimental protocols. A wide, mesoscopic channel, the OmpF porin of the outer membrane of E. coli, serves as a case study for comparing and analysing several measures of the channel cation-anion permeability in chlorides of alkali metals (LiCl, NaCl, KCl, CsCl). We show how different insights can be gained and integrated to rationalize the global image of channel selectivity. To this end, reversal potential, channel conductance and bi-ionic potential (two different salts with a common anion on each side of the channel but with the same concentration) experiments are discussed in light of an electrodiffusion model based on the Poisson-Nernst-Planck formalism. Measurements and calculations based on the atomic crystal structure of the channel show that each protocol displays a particular balance between the different sources of selectivity.

Energy variational analysis of ions in water and channels: Field theory for primitive models of complex ionic fluids

Ionic solutions are mixtures of interacting anions and cations. They hardly resemble dilute gases of uncharged noninteracting point particles described in elementary textbooks. Biological and electrochemical solutions have many components that interact strongly as they flow in concentrated environments near electrodes, ion channels, or active sites of enzymes. Interactions in concentrated environments help determine the characteristic properties of electrodes, enzymes, and ion channels. Flows are driven by a combination of electrical and chemical potentials that depend on the charges, concentrations, and sizes of all ions, not just the same type of ion. We use a variational method EnVarA (energy variational analysis) that combines Hamilton's least action and Rayleigh's dissipation principles to create a variational field theory that includes flow, friction, and complex structure with physical boundary conditions. EnVarA optimizes both the action integral functional of classical mechanics and the dissipation functional. These functionals can include entropy and dissipation as well as potential energy. The stationary point of the action is determined with respect to the trajectory of particles. The stationary point of the dissipation is determined with respect to rate functions (such as velocity). Both variations are written in one Eulerian (laboratory) framework. In variational analysis, an "extra layer" of mathematics is used to derive partial differential equations. Energies and dissipations of different components are combined in EnVarA and Euler-Lagrange equations are then derived. These partial differential equations are the unique consequence of the contributions of individual components. The form and parameters of the partial differential equations are determined by algebra without additional physical content or assumptions. The partial differential equations of mixtures automatically combine physical properties of individual (unmixed) components. If a new component is added to the energy or dissipation, the Euler-Lagrange equations change form and interaction terms appear without additional adjustable parameters. EnVarA has previously been used to compute properties of liquid crystals, polymer fluids, and electrorheological fluids containing solid balls and charged oil droplets that fission and fuse. Here we apply EnVarA to the primitive model of electrolytes in which ions are spheres in a frictional dielectric. The resulting Euler-Lagrange equations include electrostatics and diffusion and friction. They are a time dependent generalization of the Poisson-Nernst-Planck equations of semiconductors, electrochemistry, and molecular biophysics. They include the finite diameter of ions. The EnVarA treatment is applied to ions next to a charged wall, where layering is observed. Applied to an ion channel, EnVarA calculates a quick transient pile-up of electric charge, transient and steady flow through the channel, stationary "binding" in the channel, and the eventual accumulation of salts in "unstirred layers" near channels. EnVarA treats electrolytes in a unified way as complex rather than simple fluids. Ad hoc descriptions of interactions and flow have been used in many areas of science to deal with the nonideal properties of electrolytes. It seems likely that the variational treatment can simplify, unify, and perhaps derive and improve those descriptions.

Self-consistent analytic solution for the current and the access resistance in open ion channels

A self-consistent analytic approach is introduced for the estimation of the access resistance and the current through an open ion channel for an arbitrary number of species. For an ion current flowing radially inward from infinity to the channel mouth, the Poisson-Boltzmann-Nernst-Planck equations are solved analytically in the bulk with spherical symmetry in three dimensions, by linearization. Within the channel, the Poisson-Nernst-Planck equation is solved analytically in a one-dimensional approximation. An iterative procedure is used to match the two solutions together at the channel mouth in a self-consistent way. It is shown that the current-voltage characteristics obtained are in good quantitative agreement with experimental measurements.

Rectification in synthetic conical nanopores: a one-dimensional Poisson-Nernst-Planck model

Ion transport in biological and synthetic nanochannels is characterized by phenomena such as ion current fluctuations and rectification. Recently, it has been demonstrated that nanofabricated synthetic pores can mimic transport properties of biological ion channels [P. Yu. Apel, Nucl. Instrum Methods Phys. Res. B 184, 337 (2001); Z. Siwy, Europhys. Lett. 60, 349 (2002)]. Here, the ion current rectification is studied within a reduced one-dimensional (1D) Poisson-Nernst-Planck (PNP) model of synthetic nanopores. A conical channel of a few nm to a few hundred nm in diameter, and of a few mum long is considered in the limit where the channel length considerably exceeds the Debye screening length. The rigid channel wall is assumed to be weakly charged. A one-dimensional reduction of the three-dimensional problem in terms of corresponding entropic effects is put forward. The ion transport is described by the nonequilibrium steady-state solution of the 1D Poisson-Nernst-Planck system within a singular perturbation treatment. An analytic formula for the approximate rectification current in the lowest order perturbation theory is derived. A detailed comparison between numerical results and the singular perturbation theory is presented. The crucial importance of the asymmetry in the potential jumps at the pore ends on the rectification effect is demonstrated. This so constructed 1D theory is shown to describe well the experimental data in the regime of small-to-moderate electric currents.

Poisson-Nernst-Planck model of ion current rectification through a nanofluidic diode

We have investigated ion current rectification properties of a recently prepared bipolar nanofluidic diode. This device is based on a single conically shaped nanopore in a polymer film whose pore walls contain a sharp boundary between positively and negatively charged regions. A semiquantitative model that employs Poisson and Nernst-Planck equations predicts current-voltage curves as well as ionic concentrations and electric potential distributions in this system. We show that under certain conditions the rectification degree, defined as a ratio of currents recorded at the same voltage but opposite polarities, can reach values of over 1000 at a voltage range <-2 V,+2 V>. The role of thickness and position of the transition zone on the ion current rectification is discussed as well. We also show that the rectification degree scales with the applied voltage.

Electrostatic properties and macroscopic electrodiffusion in OmpF porin and mutants

The bacterial porin OmpF found in the outer membrane of E. coli is a wide channel, characterized by its poor selectivity and almost no ion specificity. It has an asymmetric structure, with relatively large entrances and a narrow constriction. By applying continuum electrostatic methods we determine the ionization states of titratable amino acid residues in the protein and calculate self-consistently the electric potential 3-D distribution within the channel. The average electrostatic properties are then represented by an effective fixed charge distribution along the pore which is the input for a macroscopic electrodiffusion model. The theoretical predictions agree with measurements performed under different salt gradients and pH. The sensitivity of reversal potential and conductance to the direction of the salt gradient and the solution pH is captured by the model. The theory is also able to explain the influence of the lipid membrane charge. The same methodology is satisfactorily applied to some OmpF mutants involving slight structural changes but a large number of net charges. The correlation found between atomic structure and ionic selectivity shows that the transport characteristics of wide channels like OmpF and its mutants are mainly regulated by the collective action of a large number of residues, rather than by the specific interactions of residues at particular locations.

Computing numerically the access resistance of a pore

The access resistance (AR) of a channel is an important component of the conductance of ion channels, particularly in wide and short channels, where it accounts for a substantial fraction of the total resistance to the movement of ions. The AR is usually calculated by using a classical and simple expression derived by Hall from electrostatics (J.E. Hall 1975 J. Gen. Phys. 66:531-532), though other expressions, both analytical and numerical, have been proposed. Here we report some numerical results for the AR of a channel obtained by solving the Poisson-Nernst-Planck equations at the entrance of a circular pore. Agreement is found between numerical calculations and analytical results from Hall's equation for uncharged pores in neutral membranes. However, for channels embedded in charged membranes, Hall's expression overestimates the AR, which is much lower and can even be neglected in some cases. The weak dependence of AR on the pore radius for charged membranes at low salt concentration can be exploited to separate the channel and the access contributions to the measured conductance.

Improved 3D continuum calculations of ion flux through membrane channels

A continuum model, based on the Poisson-Nernst-Planck (PNP) theory, is applied to simulate steady-state ion flux through protein channels. The PNP equations are modified to explicitly account (1) for the desolvation of mobile ions in the membrane pore and (2) for effects related to ion sizes. The proposed algorithm for a three-dimensional self-consistent solution of PNP equations, in which final results are refined by a focusing technique, is shown to be suitable for arbitrary channel geometry and arbitrary protein charge distribution. The role of the pore shape and protein charge distribution in formation of basic electrodiffusion properties, such as channel conductivity and selectivity, as well as concentration distributions of mobile ions in the pore region, are illustrated by simulations on model channels. The influence of the ionic strength in the bulk solution and of the externally applied electric field on channel properties are also discussed.

Synthetic nanopores with fixed charges: an electrodiffusion model for ionic transport

Synthetic nanopores with fixed charges exhibit ionic equilibrium and transport properties that resemble those displayed by biological ion channels. We present an electrodiffusion model based on the Nernst-Planck flux equations, which allows for a qualitative description of the steady state ionic transport through a nanopore when the membrane fixed charges and all mobile carriers (including the water ions) are properly taken into account. In particular, we study the current-voltage curve, the electrical conductance, the reversal potential (a measure of the nanopore ionic selectivity), as well as the flux inhibition by protons and divalent cations in the nanopore. The model clearly shows how the changes in the ionization state of the fixed charges with pH and salt concentration dictate the electrical properties of the nanopore. The agreement between the model predictions and previous experimental data allows us to identify which are the main characteristics that permit a simple description of this complex system.