Interacting ions in biophysics: real is not ideal

Setting Boundaries for Statistical Mechanics

Statistical mechanics has grown without bounds in space. Statistical mechanics of noninteracting point particles in an unbounded perfect gas is widely used to describe liquids like concentrated salt solutions of life and electrochemical technology, including batteries. Liquids are filled with interacting molecules. A perfect gas is a poor model of a liquid. Statistical mechanics without spatial bounds is impossible as well as imperfect, if molecules interact as charged particles, as nearly all atoms do. The behavior of charged particles is not defined until boundary structures and values are defined because charges are governed by Maxwell's partial differential equations. Partial differential equations require boundary structures and conditions. Boundary conditions cannot be defined uniquely 'at infinity' because the limiting process that defines 'infinity' includes such a wide variety of structures and behaviors, from elongated ellipses to circles, from light waves that never decay, to dipolar fields that decay steeply, to Coulomb fields that hardly decay at all. Boundaries and boundary conditions needed to describe matter are not prominent in classical statistical mechanics. Statistical mechanics of bounded systems is described in the EnVarA system of variational mechanics developed by Chun Liu, more than anyone else. EnVarA treatment does not yet include Maxwell equations.

The Role of Energy Scales for the Structure of Ionic Liquids at Electrified Interfaces: A Theory-Based Approach

Ionic liquids offer unique bulk and interfacial characteristics as battery electrolytes. Our continuum approach naturally describes the electrolyte on a macroscale. An integral formulation for the molecular repulsion, which can be quantitatively determined by both experimental and theoretical methods, models the electrolyte on the nanoscale. In this article, we perform a systematic series expansion of this integral formulation, derive a description of chemical potentials in terms of higher-order concentration gradients, and rationalize the appearance of fourth-order derivative operators in modified Poisson equations, as recently proposed in this context. In this way, we formulate a rigorous multiscale methodology from atomistic quantum chemistry calculations to phenomenological continuum models. We apply our generalized framework to ionic liquids near electrified interfaces and perform analytical asymptotic analysis. Three energy scales describing electrostatic forces between ions, molecular repulsion, and thermal motion determine the shape and width of the long-ranging charged double layer. We classify the charge screening mechanisms dependent on the system parameters as dielectricity, ion size, interaction strength, and temperature. We find that the charge density of electrochemical double layers in ionic liquids either decays exponentially, for negligible molecular repulsion, or oscillates continuously. Charge ordering across several ion diameters occurs if the repulsion between molecules is comparable with thermal energy and Coulomb interactions. Eventually, phase separation of the bulk electrolyte into ionic layers emerges once the molecular repulsion becomes dominant. Our framework predicts the exact phase boundaries among these three phases as a function of temperature, dielectricity, and ion size.

Dielectric Properties of Water in Charged Nanopores

In this study, we examine the spectral dielectric properties of liquid water in charged nanopores over a wide range of frequencies (0.3 GHz to 30 THz) and pore widths (0.3 to 5 nm). This has been achieved using classical molecular dynamics simulations of hydrated Na-smectite, the prototypical swelling clay mineral. We observe a drastic (20-fold) and anisotropic decrease in the static relative permittivity of the system as the pore width decreases. This large decrement in static permittivity reflects a strong attenuation of the main Debye relaxation mode of liquid water. Remarkably, this strong attenuation entails very little change in the time scale of the collective relaxation. Our results indicate that water confined in charged nanopores is a distinct solvent with a much weaker collective nature than bulk liquid water, in agreement with recent observations of water in uncharged nanopores. Finally, we observe remarkable agreement between the dielectric properties of the simulated clay system against a compiled set of soil samples at various volumetric water contents. This implies that saturation may not be the sole property dictating the dielectric properties of soil samples, rather that the pore-size distribution of fully saturated nanopores may also play a critically important role.

Direct Measurement of Interhelical DNA Repulsion and Attraction by Quantitative Cross-Linking

To better understand the forces that mediate nucleic acid compaction in biology, we developed the disulfide cross-linking approach xHEED (X-linking of Helices to measure Electrostatic Effects at Distance) to measure the distance-dependent encounter frequency of two DNA helices in solution. Using xHEED, we determined the distance that the electrostatic potential extends from DNA helices, the dependence of this distance on ionic conditions, and the magnitude of repulsion when two helices approach one another. Across all conditions tested, the potential falls to that of the bulk solution within 15 Å of the major groove surface. For separations of ∼30 Å, we measured a repulsion of 1.8 kcal/mol in low monovalent ion concentration (30 mM Na+), with higher Na+ concentrations ameliorating this repulsion, and 2 M Na+ or 100 mM Mg2+ eliminating it. Strikingly, we found full screening at very low Co3+ concentrations and net attraction at higher concentrations, without the higher-order DNA condensation that typically complicates studies of helical attraction. Our measurements define the relevant distances for electrostatic interactions of nucleic-acid helices in biology and introduce a new method to propel further understanding of how these forces impact biological processes.

Underscreening and hidden ion structures in large scale simulations of concentrated electrolytes

The electrostatic screening length predicted by Debye-Hückel theory decreases with increasing ionic strength, but recent experiments have found that the screening length can instead increase in concentrated electrolytes. This phenomenon, referred to as underscreening, is believed to result from ion-ion correlations and short-range forces such as excluded volume interactions among ions. We use Brownian Dynamics to simulate a version of the Restrictive Primitive Model for electrolytes over a wide range of ion concentrations, ionic strengths, and ion excluded volume radii for binary electrolytes. We measure the decay of the charge-charge correlation among ions in the bulk and compare it against scaling trends found experimentally and determined in certain weak coupling theories of ion-ion correlation. Moreover, we find that additional large scale ion structures emerge at high concentrations. In this regime, the frequency of oscillations computed from the charge-charge correlation function is not dominated by electrostatic interactions but rather by excluded volume interactions and with oscillation periods on the order of the ion diameter. We also find the nearest neighbor correlation of ions sharing the same charge transitions from negative at small concentrations to positive at high concentrations, representing the formation of small, like-charge ion clusters. We conclude that the increase in local charge density due to the formation of these clusters and the topological constraints of macroscopic charged surfaces can help explain the degree of underscreening observed experimentally.

Long-Range Electrostatic Colloidal Interactions and Specific Ion Effects in Deep Eutectic Solvents

While the traditional consensus dictates that high ion concentrations lead to negligible long-range electrostatic interactions, we demonstrate that electrostatic correlations prevail in deep eutectic solvents where intrinsic ion concentrations often surpass 2.5 M. Here we present an investigation of intermicellar interactions in 1:2 choline chloride:glycerol and 1:2 choline bromide:glycerol using small-angle neutron scattering. Our results show that long-range electrostatic repulsions between charged colloidal particles occur in these solvents. Interestingly, micelle morphology and electrostatic interactions are modulated by specific counterion condensation at the micelle interface despite the exceedingly high concentration of the native halide from the solvent. This modulation follows the trends described by the Hofmeister series for specific ion effects. The results are rationalized in terms of predominant ion-ion correlations, which explain the reduction in the effective ionic strength of the continuum and the observed specific ion effects.

Electrodiffusion Phenomena in Neuroscience and the Nernst–Planck–Poisson Equations

This work is aimed to give an electrochemical insight into the ionic transport phenomena in the cellular environment of organized brain tissue. The Nernst–Planck–Poisson (NPP) model is presented, and its applications in the description of electrodiffusion phenomena relevant in nanoscale neurophysiology are reviewed. These phenomena include: the signal propagation in neurons, the liquid junction potential in extracellular space, electrochemical transport in ion channels, the electrical potential distortions invisible to patch-clamp technique, and calcium transport through mitochondrial membrane. The limitations, as well as the extensions of the NPP model that allow us to overcome these limitations, are also discussed.

Surface Forces and Structure in a Water-in-Salt Electrolyte

Water-in-salt electrolytes are a fascinating new class of highly concentrated aqueous solutions with wide electrochemical stability windows that make them viable as aqueous battery electrolytes. However, the high ion concentration of water-in-salt electrolytes means that these systems are poorly understood when compared to more dilute electrolyte solutions. Here, we present direct surface force measurements across thin films of a water-in-salt electrolyte at several concentrations. We find that the electrolyte adopts a layered structure at charged interfaces composed of a nanostructure of a hydrated cation and nonaqueous anion-rich domains. These observations will aid in the interpretation of capacitance and double-layer behavior of water-in-salt electrolytes with consequences for their use in energy storage devices.

A Bidomain Model for Lens Microcirculation

There exists a large body of research on the lens of the mammalian eye over the past several decades. The objective of this work is to provide a link between the most recent computational models and some of the pioneering work in the 1970s and 80s. We introduce a general nonelectroneutral model to study the microcirculation in the lens of the eye. It describes the steady-state relationships among ion fluxes, between water flow and electric field inside cells, and in the narrow extracellular spaces between cells in the lens. Using asymptotic analysis, we derive a simplified model based on physiological data and compare our results with those in the literature. We show that our simplified model can be reduced further to the first-generation models, whereas our full model is consistent with the most recent computational models. In addition, our simplified model captures in its equations the main features of the full computational models. Our results serve as a useful link intermediate between the computational models and the first-generation analytical models. Simplified models of this sort may be particularly helpful as the roles of similar osmotic pumps of microcirculation are examined in other tissues with narrow extracellular spaces, such as cardiac and skeletal muscle, liver, kidney, epithelia in general, and the narrow extracellular spaces of the central nervous system, the "brain." Simplified models may reveal the general functional plan of these systems before full computational models become feasible and specific.

Poisson-Fermi modeling of ion activities in aqueous single and mixed electrolyte solutions at variable temperature

The combinatorial explosion of empirical parameters in tens of thousands presents a tremendous challenge for extended Debye-Hückel models to calculate activity coefficients of aqueous mixtures of the most important salts in chemistry. The explosion of parameters originates from the phenomenological extension of the Debye-Hückel theory that does not take steric and correlation effects of ions and water into account. By contrast, the Poisson-Fermi theory developed in recent years treats ions and water molecules as nonuniform hard spheres of any size with interstitial voids and includes ion-water and ion-ion correlations. We present a Poisson-Fermi model and numerical methods for calculating the individual or mean activity coefficient of electrolyte solutions with any arbitrary number of ionic species in a large range of salt concentrations and temperatures. For each activity-concentration curve, we show that the Poisson-Fermi model requires only three unchanging parameters at most to well fit the corresponding experimental data. The three parameters are associated with the Born radius of the solvation energy of an ion in electrolyte solution that changes with salt concentrations in a highly nonlinear manner.

Relative dielectric constants and selectivity ratios in open ionic channels

We investigate the effects of the relative dielectric coefficient on ionic flows in open ion channels, using mathematical analysis of reasonably general Poisson-Nernst-Planck type models that can include the finite sizes of ions. The value of the relative dielectric coefficient is of course a crucial parameter for ionic behavior in general. Using the powerful theory of singularly perturbed problems in applied mathematics, we show that some properties of open channels are quite insensitive to variation in the relative dielectric coefficient, thereby explaining such effects seen unexpectedly in simulations. The ratio between the total number of one ion species and that of another ion species, and the ratio between the flux of one ion species and that of another ion species do not depend significantly on the relative dielectric coefficient.

Dynamics of Current, Charge and Mass

Electricity plays a special role in our lives and life. The dynamics of electrons allow light to flow through a vacuum. The equations of electron dynamics are nearly exact and apply from nuclear particles to stars. These Maxwell equations include a special term, the displacement current (of a vacuum). The displacement current allows electrical signals to propagate through space. Displacement current guarantees that current is exactly conserved from inside atoms to between stars, as long as current is defined as the entire source of the curl of the magnetic field, as Maxwell did.We show that the Bohm formulation of quantum mechanics allows the easy definition of the total current, and its conservation, without the dificulties implicit in the orthodox quantum theory. The orthodox theory neglects the reality of magnitudes, like the currents, during times that they are not being explicitly measured.We show how conservation of current can be derived without mention of the polarization or dielectric properties of matter. We point out that displacement current is handled correctly in electrical engineering by ‘stray capacitances’, although it is rarely discussed explicitly. Matter does not behave as physicists of the 1800’s thought it did. They could only measure on a time scale of seconds and tried to explain dielectric properties and polarization with a single dielectric constant, a real positive number independent of everything. Matter and thus charge moves in enormously complicated ways that cannot be described by a single dielectric constant,when studied on time scales important today for electronic technology and molecular biology. When classical theories could not explain complex charge movements, constants in equations were allowed to vary in solutions of those equations, in a way not justified by mathematics, with predictable consequences. Life occurs in ionic solutions where charge is moved by forces not mentioned or described in the Maxwell equations, like convection and diffusion. These movements and forces produce crucial currents that cannot be described as classical conduction or classical polarization. Derivations of conservation of current involve oversimplified treatments of dielectrics and polarization in nearly every textbook. Because real dielectrics do not behave in that simple way-not even approximately-classical derivations of conservation of current are often distrusted or even ignored. We show that current is conserved inside atoms. We show that current is conserved exactly in any material no matter how complex are the properties of dielectric, polarization, or conduction currents. Electricity has a special role because conservation of current is a universal law.Most models of chemical reactions do not conserve current and need to be changed to do so. On the macroscopic scale of life, conservation of current necessarily links far spread boundaries to each other, correlating inputs and outputs, and thereby creating devices.We suspect that correlations created by displacement current link all scales and allow atoms to control the machines and organisms of life. Conservation of current has a special role in our lives and life, as well as in physics. We believe models, simulations, and computations should conserve current on all scales, as accurately as possible, because physics conserves current that way. We believe models will be much more successful if they conserve current at every level of resolution, the way physics does.We surely need successful models as we try to control macroscopic functions by atomic interventions, in technology, life, and medicine. Maxwell’s displacement current lets us see stars. We hope it will help us see how atoms control life.

Electrodiffusion phenomena in neuroscience: a neglected companion

The emerging technological revolution in genetically encoded molecular sensors and super-resolution imaging provides neuroscientists with a pass to the real-time nano-world. On this small scale, however, classical principles of electrophysiology do not always apply. This is in large part because the nanoscopic heterogeneities in ionic concentrations and the local electric fields associated with individual ions and their movement can no longer be ignored. Here, we review basic principles of molecular electrodiffusion in the cellular environment of organized brain tissue. We argue that accurate interpretation of physiological observations on the nanoscale requires a better understanding of the underlying electrodiffusion phenomena.

Long range electrostatic forces in ionic liquids

Experimental evidence for long range surface forces in ionic liquids is collated and examined, key outstanding questions are identified, and possible mechanisms underpinning these long range forces are explored.

Pore size matters for potassium channel conductance

Ion channels are membrane proteins that mediate efficient ion transport across the hydrophobic core of cell membranes, an unlikely process in their absence. K⁺ channels discriminate K⁺ over cations with similar radii with extraordinary selectivity and display a wide diversity of ion transport rates, covering differences of two orders of magnitude in unitary conductance. The pore domains of large- and small-conductance K⁺ channels share a general architectural design comprising a conserved narrow selectivity filter, which forms intimate interactions with permeant ions, flanked by two wider vestibules toward the internal and external openings. In large-conductance K⁺ channels, the inner vestibule is wide, whereas in small-conductance channels it is narrow. Here we raise the idea that the physical dimensions of the hydrophobic internal vestibule limit ion transport in K⁺ channels, accounting for their diversity in unitary conductance.

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.

Introduction of effective dielectric constant to the Poisson-Nernst-Planck model

The Poisson-Nernst-Planck (PNP) model has been widely used for analyzing impedance or dielectric spectra observed for dilute electrolytic cells. In the analysis, the behavior of mobile ions in the cell under an external electric field has been explained by a conductive nature regardless of ionic concentrations. However, if the cell has parallel-plate blocking electrodes, the mobile ions may also play a role as a dielectric medium in the cell by the effect of space-charge polarization when the ionic concentration is sufficiently low. Thus the mobile ions confined between the blocking electrodes can have conductive and dielectric natures simultaneously, and their intensities are affected by the ionic concentration and the adsorption of solvent molecules on the electrodes. The balance of the conductive and dielectric natures is quantitatively determined by introducing an effective dielectric constant to the PNP model in the data analysis. The generalized PNP model with the effective dielectric constant successfully explains the anomalous frequency-dependent dielectric behaviors brought about by the mobile ions in dilute electrolytic cells, for which the conventional PNP model fails in interpretation.

Numerical methods for a Poisson-Nernst-Planck-Fermi model of biological ion channels

Numerical methods are proposed for an advanced Poisson-Nernst-Planck-Fermi (PNPF) model for studying ion transport through biological ion channels. PNPF contains many more correlations than most models and simulations of channels, because it includes water and calculates dielectric properties consistently as outputs. This model accounts for the steric effect of ions and water molecules with different sizes and interstitial voids, the correlation effect of crowded ions with different valences, and the screening effect of polarized water molecules in an inhomogeneous aqueous electrolyte. The steric energy is shown to be comparable to the electrical energy under physiological conditions, demonstrating the crucial role of the excluded volume of particles and the voids in the natural function of channel proteins. Water is shown to play a critical role in both correlation and steric effects in the model. We extend the classical Scharfetter-Gummel (SG) method for semiconductor devices to include the steric potential for ion channels, which is a fundamental physical property not present in semiconductors. Together with a simplified matched interface and boundary (SMIB) method for treating molecular surfaces and singular charges of channel proteins, the extended SG method is shown to exhibit important features in flow simulations such as optimal convergence, efficient nonlinear iterations, and physical conservation. The generalized SG stability condition shows why the standard discretization (without SG exponential fitting) of NP equations may fail and that divalent Ca(2+) may cause more unstable discrete Ca(2+) fluxes than that of monovalent Na(+). Two different methods-called the SMIB and multiscale methods-are proposed for two different types of channels, namely, the gramicidin A channel and an L-type calcium channel, depending on whether water is allowed to pass through the channel. Numerical methods are first validated with constructed models whose exact solutions are known. The experimental data of both channels are then used to verify and explain novel features of PNPF as compared with previous PNP models. The PNPF currents are in accord with the experimental I-V (V for applied voltages) data of the gramicidin A channel and I-C (C for bath concentrations) data of the calcium channel with 10(-8)-fold bath concentrations that pose severe challenges in theoretical simulations.

Systematic interpretation of differential capacitance data

Differential capacitance (DC) data have been widely used to characterize the structure of electrolyte solutions near charged interfaces and as experimental validation of models for electrolyte structure. Fixing a large class of models of electrolyte free energy that incorporate finite-volume effects, a reduction is identified which permits the identification of all free energies within that class that return identical DC data. The result is an interpretation of DC data through the equivalence classes of nonideality terms, and associated boundary layer structures, that cannot be differentiated by DC data. Specifically, for binary salts, DC data, even if measured over a range of ionic concentrations, are unable to distinguish among models which exhibit charge asymmetry, charge reversal, and even ion crowding. The reduction applies to capacitors which are much wider than the associated Debye length and to finite-volume terms that are algebraic in charge density. However, within these restrictions the free energy is shown to be uniquely identified if the DC data are supplemented with measurements of the excess chemical potential of the system in the bulk state.

Long-range electrostatic screening in ionic liquids

Electrolyte solutions with high concentrations of ions are prevalent in biological systems and energy storage technologies. Nevertheless, the high interaction free energy and long-range nature of electrostatic interactions makes the development of a general conceptual picture of concentrated electrolytes a significant challenge. In this work, we study ionic liquids, single-component liquids composed solely of ions, in an attempt to provide a novel perspective on electrostatic screening in very high concentration (nonideal) electrolytes. We use temperature-dependent surface force measurements to demonstrate that the long-range, exponentially decaying diffuse double-layer forces observed across ionic liquids exhibit a pronounced temperature dependence: Increasing the temperature decreases the measured exponential (Debye) decay length, implying an increase in the thermally driven effective free-ion concentration in the bulk ionic liquids. We use our quantitative results to propose a general model of long-range electrostatic screening in ionic liquids, where thermally activated charge fluctuations, either free ions or correlated domains (quasiparticles), take on the role of ions in traditional dilute electrolyte solutions. This picture represents a crucial step toward resolving several inconsistencies surrounding electrostatic screening and charge transport in ionic liquids that have impeded progress within the interdisciplinary ionic liquids community. More broadly, our work provides a previously unidentified way of envisioning highly concentrated electrolytes, with implications for diverse areas of inquiry, ranging from designing electrochemical devices to rationalizing electrostatic interactions in biological systems.

Permeation redux: thermodynamics and kinetics of ion movement through potassium channels

The fundamental biophysics underlying the selective movement of ions through ion channels was launched by George Eisenman in the 1960s, using glass electrodes. This minireview examines the insights from these early studies and the explosive progress made since then.

Practical Aspects of Free-Energy Calculations: A Review

Free-energy calculations in the framework of classical molecular dynamics simulations are nowadays used in a wide range of research areas including solvation thermodynamics, molecular recognition, and protein folding. The basic components of a free-energy calculation, that is, a suitable model Hamiltonian, a sampling protocol, and an estimator for the free energy, are independent of the specific application. However, the attention that one has to pay to these components depends considerably on the specific application. Here, we review six different areas of application and discuss the relative importance of the three main components to provide the reader with an organigram and to make nonexperts aware of the many pitfalls present in free energy calculations.

Understanding nucleic acid-ion interactions

Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions.

Mechanisms of a human skeletal myotonia produced by mutation in the C-terminus of NaV1.4: is Ca2+ regulation defective?

Mutations in the cytoplasmic tail (CT) of voltage gated sodium channels cause a spectrum of inherited diseases of cellular excitability, yet to date only one mutation in the CT of the human skeletal muscle voltage gated sodium channel (hNa_V1.4_F1705I) has been linked to cold aggravated myotonia. The functional effects of altered regulation of hNa_V1.4_F1705I are incompletely understood. The location of the hNa_V1.4_F1705I in the CT prompted us to examine the role of Ca²⁺ and calmodulin (CaM) regulation in the manifestations of myotonia. To study Na channel related mechanisms of myotonia we exploited the differences in rat and human Na_V1.4 channel regulation by Ca²⁺ and CaM. hNa_V1.4_F1705I inactivation gating is Ca²⁺-sensitive compared to wild type hNa_V1.4 which is Ca²⁺ insensitive and the mutant channel exhibits a depolarizing shift of the V_1/2 of inactivation with CaM over expression. In contrast the same mutation in the rNa_V1.4 channel background (rNa_V1.4_F1698I) eliminates Ca²⁺ sensitivity of gating without affecting the CaM over expression induced hyperpolarizing shift in steady-state inactivation. The differences in the Ca²⁺ sensitivity of gating between wild type and mutant human and rat Na_V1.4 channels are in part mediated by a divergence in the amino acid sequence in the EF hand like (EFL) region of the CT. Thus the composition of the EFL region contributes to the species differences in Ca²⁺/CaM regulation of the mutant channels that produce myotonia. The myotonia mutation F1705I slows I_Na decay in a Ca²⁺-sensitive fashion. The combination of the altered voltage dependence and kinetics of I_Na decay contribute to the myotonic phenotype and may involve the Ca²⁺-sensing apparatus in the CT of Na_V1.4.

Energetics of discrete selectivity bands and mutation-induced transitions in the calcium-sodium ion channels family

We use Brownian dynamics (BD) simulations to study the ionic conduction and valence selectivity of a generic electrostatic model of a biological ion channel as functions of the fixed charge Q(f) at its selectivity filter. We are thus able to reconcile the discrete calcium conduction bands recently revealed in our BD simulations, M0 (Q(f)=1e), M1 (3e), M2 (5e), with a set of sodium conduction bands L0 (0.5e), L1 (1.5e), thereby obtaining a completed pattern of conduction and selectivity bands vs Q(f) for the sodium-calcium channels family. An increase of Q(f) leads to an increase of calcium selectivity: L0 (sodium-selective, nonblocking channel) → M0 (nonselective channel) → L1 (sodium-selective channel with divalent block) → M1 (calcium-selective channel exhibiting the anomalous mole fraction effect). We create a consistent identification scheme where the L0 band is putatively identified with the eukaryotic sodium channel The scheme created is able to account for the experimentally observed mutation-induced transformations between nonselective channels, sodium-selective channels, and calcium-selective channels, which we interpret as transitions between different rows of the identification table. By considering the potential energy changes during permeation, we show explicitly that the multi-ion conduction bands of calcium and sodium channels arise as the result of resonant barrierless conduction. The pattern of periodic conduction bands is explained on the basis of sequential neutralization taking account of self-energy, as Q(f)(z,i)=ze(1/2+i), where i is the order of the band and z is the valence of the ion. Our results confirm the crucial influence of electrostatic interactions on conduction and on the Ca(2+)/Na(+) valence selectivity of calcium and sodium ion channels. The model and results could be also applicable to biomimetic nanopores with charged walls.