Proteins, channels and crowded ions

New insights into the effects of small permanent charge on ionic flows: A higher order analysis

This study investigated how permanent charges influence the dynamics of ionic channels. Using a quasi-one-dimensional classical Poisson-Nernst-Planck (PNP) model, we investigated the behavior of two distinct ion species-one positively charged and the other negatively charged. The spatial distribution of permanent charges was characterized by zero values at the channel ends and a constant charge $ Q_0 $ within the central region. By treating the classical PNP model as a boundary value problem (BVP) for a singularly perturbed system, the singular orbit of the BVP depended on $ Q_0 $ in a regular way. We therefore explored the solution space in the presence of a small permanent charge, uncovering a systematic dependence on this parameter. Our analysis employed a rigorous perturbation approach to reveal higher-order effects originating from the permanent charges. Through this investigation, we shed light on the intricate interplay among boundary conditions and permanent charges, providing insights into their impact on the behavior of ionic current, fluxes, and flux ratios. We derived the quadratic solutions in terms of permanent charge, which were notably more intricate compared to the linear solutions. Through computational tools, we investigated the impact of these quadratic solutions on fluxes, current-voltage relations, and flux ratios, conducting a thorough analysis of the results. These novel findings contributed to a deeper comprehension of ionic flow dynamics and hold potential implications for enhancing the design and optimization of ion channel-based technologies.

Spontaneous pattern formation in monolayers of binary mixtures with competing interactions

A model for a monolayer of two types of particles spontaneously forming ordered patterns is studied using a mesoscopic theory and MC simulations. We assume hard-cores of the same size a for both components. For r > a, like particles attract and repel each other at short and large distances, respectively, with the same potential u(r) for both species, and the cross-interaction is -u(r). The model is inspired by oppositely charged particles or macromolecules with preferential solubility in different components of a solvent that is close to a miscibility critical point, in particular by inclusions in biological membranes. We obtain the phase diagram in the chemical potentials and temperature variables as well as in the concentration, density and temperature variables, using the mean-field one-shell approximation. We find that the presence of the second component significantly extends the temperature range of stability of the ordered phases. We obtain three stable phases with periodic concentration: the lamellar L phase with alternating stripes of the two components for similar chemical potentials, and a hexagonal arrangement of the clusters of the minority component in the liquid of the majority component. The latter two phases, however, are stable only at relatively high temperatures. At lower temperatures, the L phase coexists with a disordered one-component fluid or with very dilute gas with mixed components. At still lower temperatures, the one-component phase coexisting with the L phase can be disordered or ordered, depending on the chemical potentials. The theoretical results are confirmed by MC simulations for selected thermodynamic states.

Finite ion size effects on I-V relations via Poisson-Nernst-Planck systems with two cations: A case study

We consider a quasi-one-dimensional Poisson-Nernst-Planck model with two cations having the same valances and one anion. Bikerman's local hard-sphere potential is included to account for ion size effects. Under some further restrictions on the boundary conditions of the two cations, we obtain approximations of the I-V (current-voltage) relations by treating the ion sizes as small parameters. Critical potentials are identified, which play critical roles in characterizing finite ion size effects on ionic flows. Nonlinear interplays between system parameters, such as boundary concentrations and diffusion coefficients, are analyzed. To provide more intuitive illustrations of our analytical results and better understanding of the dynamics of ionic flows through membrane channels, numerical simulations are performed.

Mathematical Analysis on Current-Voltage Relations via Classical Poisson-Nernst-Planck Systems with Nonzero Permanent Charges under Relaxed Electroneutrality Boundary Conditions

We focus on a quasi-one-dimensional Poisson-Nernst-Planck model with small permanent charges for ionic flows of two oppositely charged ion species through an ion channel. Of particular interest is to examine the dynamics of ionic flows in terms of I-V (current-voltage) relations with boundary layers due to the relaxation of neutral conditions on boundary concentrations. This is achieved by employing the regular perturbation analysis on the solutions established through geometric singular perturbation analysis. Rich dynamics are observed, particularly, the nonlinear interplays among different physical parameters are characterized. Critical potentials are identified, which play critical roles in the study of ionic flows and can be estimated experimentally. Numerical simulations are performed to further illustrate and provide more intuitive understandings of our analytical results.

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.

Competition between Cations via Classical Poisson-Nernst-Planck Models with Nonzero but Small Permanent Charges

We study a one-dimensional Poisson–Nernst–Planck system for ionic flow through a membrane channel. Nonzero but small permanent charge, the major structural quantity of an ion channel, is included in the model. Two cations with the same valences and one anion are included in the model, which provides more rich and complicated correlations/interactions between ions. The cross-section area of the channel is included in the system, and it provides certain information of the geometry of the three-dimensional channel, which is critical for our analysis. Geometric singular perturbation analysis is employed to establish the existence and local uniqueness of solutions to the system for small permanent charges. Treating the permanent charge as a small parameter, through regular perturbation analysis, we are able to derive approximations of the individual fluxes explicitly, and this allows us to study the competition between two cations, which is related to the selectivity phenomena of ion channels. Numerical simulations are performed to provide a more intuitive illustration of our analytical results, and they are consistent.

Dynamics of ionic flows via Poisson-Nernst-Planck systems with local hard-sphere potentials: Competition between cations

We study a quasi-one-dimensional steady-state Poisson-Nernst-Planck type model for ionic flows through a membrane channel with three ion species, two positively charged with the same valence and one negatively charged. Bikerman's local hard-sphere potential is included in the model to account for ion sizes. The problem is treated as a boundary value problem of a singularly perturbed differential system. Under the framework of a geometric singular perturbation theory, together with specific structures of this concrete model, the existence and uniqueness of solutions to the boundary value problem for small ion sizes is established. Furthermore, treating the ion sizes as small parameters, we derive an approximation of individual fluxes, from which one can further study the qualitative properties of ionic flows and extract concrete information directly related to biological measurements. Of particular interest is the competition between two cations due to the nonlinear interplay between finite ion sizes, diffusion coefficients and boundary conditions, which is closely related to selectivity phenomena of open ion channels with given protein structures. Furthermore, we are able to characterize the distinct effects of the nonlinear interplays between these physical parameters. Numerical simulations are performed to identify some critical potentials which play critical roles in examining properties of ionic flows in our analysis.

Molecular Mean-Field Theory of Ionic Solutions: A Poisson-Nernst-Planck-Bikerman Model

We have developed a molecular mean-field theory—fourth-order Poisson–Nernst–Planck–Bikerman theory—for modeling ionic and water flows in biological ion channels by treating ions and water molecules of any volume and shape with interstitial voids, polarization of water, and ion-ion and ion-water correlations. The theory can also be used to study thermodynamic and electrokinetic properties of electrolyte solutions in batteries, fuel cells, nanopores, porous media including cement, geothermal brines, the oceanic system, etc. The theory can compute electric and steric energies from all atoms in a protein and all ions and water molecules in a channel pore while keeping electrolyte solutions in the extra- and intracellular baths as a continuum dielectric medium with complex properties that mimic experimental data. The theory has been verified with experiments and molecular dynamics data from the gramicidin A channel, L-type calcium channel, potassium channel, and sodium/calcium exchanger with real structures from the Protein Data Bank. It was also verified with the experimental or Monte Carlo data of electric double-layer differential capacitance and ion activities in aqueous electrolyte solutions. We give an in-depth review of the literature about the most novel properties of the theory, namely Fermi distributions of water and ions as classical particles with excluded volumes and dynamic correlations that depend on salt concentration, composition, temperature, pressure, far-field boundary conditions etc. in a complex and complicated way as reported in a wide range of experiments. The dynamic correlations are self-consistent output functions from a fourth-order differential operator that describes ion-ion and ion-water correlations, the dielectric response (permittivity) of ionic solutions, and the polarization of water molecules with a single correlation length parameter.

Boundary Layer Effects on Ionic Flows Via Classical Poisson-Nernst-Planck Systems

A quasi-one-dimensional steady-state Poisson-Nernst-Planck model of two oppositely charged ion species through a membrane channel is analyzed. The model problem is treated as a boundary value problem of a singularly perturbed differential system. Our analysis is based on the geometric singular perturbation theory but, most importantly, on specific structures of this concrete model. The existence and (local ) uniqueness of solutions to the boundary value problem is established. In particular, an approximation of both the individual flux and the I-V (current-voltage) relation are derived explicitly from the zeroth order approximation (in ") solutions, from which the boundary layer effects on ionic flows are studied in great details.

Ion selectivity strategies of sodium channel selectivity filters

CONSPECTUS: Sodium ion channels selectively transport Na(+) cations across the cell membrane. These integral parts of the cell machinery are implicated in regulating the cardiac, skeletal and smooth muscle contraction, nerve impulses, salt and water homeostasis, as well as pain and taste perception. Their malfunction often results in various channelopathies of the heart, brain, skeletal muscles, and lung; thus, sodium channels are key drug targets for various disorders including cardiac arrhythmias, heart attack, stroke, migraine, epilepsy, pain, cancer, and autoimmune disorders. The ability of sodium channels to discriminate the native Na(+) among other competing ions in the surrounding fluids is crucial for proper cellular functions. The selectivity filter (SF), the narrowest part of the channel's open pore, lined with amino acid residues that specifically interact with the permeating ion, plays a major role in determining Na(+) selectivity. Different sodium channels have different SFs, which vary in the symmetry, number, charge, arrangement, and chemical type of the metal-ligating groups and pore size: epithelial/degenerin/acid-sensing ion channels have generally trimeric SFs lined with three conserved neutral serines and/or backbone carbonyls; eukaryotic sodium channels have EKEE, EEKE, DKEA, and DEKA SFs with an invariant positively charged lysine from the second or third domain; and bacterial voltage-gated sodium (Nav) channels exhibit symmetrical EEEE SFs, reminiscent of eukaryotic voltage-gated calcium channels. How do these different sodium channel SFs achieve high selectivity for Na(+) over its key rivals, K(+) and Ca(2+)? What factors govern the metal competition in these SFs and which of these factors are exploited to achieve Na(+) selectivity in the different sodium channel SFs? The free energies for replacing K(+) or Ca(2+) bound inside different model SFs with Na(+), evaluated by a combination of density functional theory and continuum dielectric calculations, have shed light on these questions. The SFs of epithelial and eukaryotic Nav channels select Na(+) by providing an optimal number and ligating strength of metal ligands as well as a rigid pore whose size fits the cognate Na(+) ideally. On the other hand, the SFs of bacterial Nav channels select Na(+), as the protein matrix attenuates ion-protein interactions relative to ion-solvent interactions by enlarging the pore and allowing water to enter, so the ion interacts indirectly with the conserved glutamates via bridging water molecules. This shows how these various SFs have adapted to the specific physicochemical properties of the native ion, using different strategies to select Na(+) among its contenders.

Energetics of Multi-Ion Conduction Pathways in Potassium Ion Channels

Potassium ion channels form pores in cell membranes, allowing potassium ions through while preventing the passage of sodium ions. Despite numerous high-resolution structures, it is not yet possible to relate their structure to their single molecule function other than at a qualitative level. Over the past decade, there has been a concerted effort using molecular dynamics to capture the thermodynamics and kinetics of conduction by calculating potentials of mean force (PMF). These can be used, in conjunction with the electro-diffusion theory, to predict the conductance of a specific ion channel. Here, we calculate seven independent PMFs, thereby studying the differences between two potassium ion channels, the effect of the CHARMM CMAP forcefield correction, and the sensitivity and reproducibility of the method. Thermodynamically stable ion-water configurations of the selectivity filter can be identified from all the free energy landscapes, but the heights of the kinetic barriers for potassium ions to move through the selectivity filter are, in nearly all cases, too high to predict conductances in line with experiment. This implies it is not currently feasible to predict the conductance of potassium ion channels, but other simpler channels may be more tractable.

Interacting ions in biophysics: real is not ideal

Ions in water are important throughout biology, from molecules to organs. Classically, ions in water were treated as ideal noninteracting particles in a perfect gas. Excess free energy of each ion was zero. Mathematics was not available to deal consistently with flows, or interactions with other ions or boundaries. Nonclassical approaches are needed because ions in biological conditions flow and interact. The concentration gradient of one ion can drive the flow of another, even in a bulk solution. A variational multiscale approach is needed to deal with interactions and flow. The recently developed energetic variational approach to dissipative systems allows mathematically consistent treatment of the bio-ions Na(+), K(+), Ca(2+), and Cl(-) as they interact and flow. Interactions produce large excess free energy that dominate the properties of the high concentration of ions in and near protein active sites, ion channels, and nucleic acids: the number density of ions is often >10 M. Ions in such crowded quarters interact strongly with each other as well as with the surrounding protein. Nonideal behavior found in many experiments has classically been ascribed to allosteric interactions mediated by the protein and its conformation changes. The ion-ion interactions present in crowded solutions-independent of conformation changes of the protein-are likely to change the interpretation of many allosteric phenomena. Computation of all atoms is a popular alternative to the multiscale approach. Such computations involve formidable challenges. Biological systems exist on very different scales from atomic motion. Biological systems exist in ionic mixtures (like extracellular and intracellular solutions), and usually involve flow and trace concentrations of messenger ions (e.g., 10(-7) M Ca(2+)). Energetic variational methods can deal with these characteristic properties of biological systems as we await the maturation and calibration of all-atom simulations of ionic mixtures and divalents.

Ionic Interactions Are Everywhere

Ionic solutions are dominated by interactions because they must be electrically neutral, but classical theory assumes no interactions. Biological solutions are rather like seawater, concentrated enough so that the diameter of ions also produces important interactions. In my view, the theory of complex fluids is needed to deal with the interacting reality of biological solutions.

Competition among Ca2+, Mg2+, and Na+ for model ion channel selectivity filters: determinants of ion selectivity

Because voltage-gated ion channels play critical biological roles, understanding how they can discriminate the native metal ion from rival cations in the milieu is of great interest. Although Ca(2+), Mg(2+), and Na(+) are present in comparable concentrations outside the cell, the factors governing the competition among these cations for the selectivity filter of voltage-gated Ca(2+) ion channel remain unclear. Using density functional theory combined with continuum dielectric methods, we evaluate the effect of (1) the number, chemical type, and charge of the ligands lining the pore, (2) the pore's rigidity, size, symmetry, and solvent accessibility, and (3) the Ca(2+) hydration number outside the selectivity filter on the competition among Ca(2+), Mg(2+), and Na(+) in model selectivity filters. The calculations show how the outcome of this competition depends on the interplay between electronic and solvation effects. Selectivity for monovalent Na(+) over divalent Ca(2+)/Mg(2+) is achieved when solvation effects outweigh electrostatic effects; thus filters comprising a few weak charge-donating groups such as Ser/Thr side chains, where electrostatic effects are relatively weak and are easily overcome by solvation effects, are Na(+)-selective. In contrast, selectivity for divalent Ca(2+)/Mg(2+) over monovalent Na(+) is achieved when metal-ligand electrostatic effects outweigh solvation effects. The key differences in selectivity between Mg(2+) and Ca(2+) lie in the pore size, oligomericity, and solvent accessibility. The results, which are consistent with available experimental data, reveal how the structure and composition of the ion channel selectivity pore had adapted to the specific physicochemical properties of the native metal ion to enhance the competitiveness of the native metal toward rival cations.

Multiple Scales in the Simulation of Ion Channels and Proteins

Computation of living processes creates great promise for the everyday life of mankind and great challenges for physical scientists. Simulations molecular dynamics have great appeal to biologists as a natural extension of structural biology. Once a biologist sees a structure, she/he wants to see it move. Molecular biology has shown that a small number of atoms, sometimes even one messenger ion, like Ca(2+), can control biological function on the scale of cells, organs, tissues, and organisms. Enormously concentrated ions-at number densities of ~20 M-in protein channels and enzymes are responsible for many of the characteristics of living systems, just as highly concentrated ions near electrodes are responsible for many of the characteristics of electrochemical systems. Here we confront the reality of the scale differences of ions. We show that the scale differences needed to simulate all the atoms of biological cells are 10(7) in linear dimension, 10(21) in three dimensions, 10(9) in resolution, 10(11) in time, and 10(13) in particle number (to deal with concentrations of Ca(2+)). These scales must be dealt with simultaneously if the simulation is to deal with most biological functions. Biological function extends across all of them, all at once in most cases. We suggest a computational approach using explicit multiscale analysis instead of implicit simulation of all scales. The approach is based on an energy variational principle EnVarA introduced by Chun Liu to deal with complex fluids. Variational methods deal automatically with multiple interacting components and scales. When an additional component is added to the system, the resulting Euler Lagrange field equations change form automatically-by algebra alone-without additional unknown parameters. Multifaceted interactions are solutions of the resulting equations. We suggest that ionic solutions should be viewed as complex fluids with simple components. Highly concentrated solutions-dominated by interactions of components-are easily computed by EnVarA. Successful computation of ions concentrated in special places may be a significant step to understanding the defining characteristics of biological and electrochemical systems. Indeed, computing ions near proteins and nucleic acids may prove as important to molecular biology and chemical technology as computing holes and electrons has been to our semiconductor and digital technology.

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.

Diffusion, exclusion, and specific binding in a large channel: a study of OmpF selectivity inversion

We find that moderate cationic selectivity of the general bacterial porin OmpF in sodium and potassium chloride solutions is inversed to anionic selectivity in concentrated solutions of barium, calcium, nickel, and magnesium chlorides. To understand the origin of this phenomenon, we consider several factors, which include the binding of divalent cations, electrostatic and steric exclusion of differently charged and differently sized ions, size-dependent hydrodynamic hindrance, electrokinetic effects, and significant "anionic" diffusion potential for bulk solutions of chlorides of divalent cations. Though all these factors contribute to the measured selectivity of this large channel, the observed selectivity inversion is mostly due to the following two. First, binding divalent cations compensates, or even slightly overcompensates, for the negative charge of the OmpF protein, which is known to be the main cause of cationic selectivity in sodium and potassium chloride solutions. Second, the higher anionic (versus cationic) transport rate expected for bulk solutions of chloride salts of divalent cations is the leading cause of the measured anionic selectivity of the channel. Interestingly, at high concentrations the binding of cations does not show any pronounced specificity within the divalent series because the reversal potentials measured in the series correlate well with the corresponding bulk diffusion potentials. Thus our study shows that, in contrast to the highly selective channels of neurophysiology that employ mostly the exclusion mechanism, quite different factors account for the selectivity of large channels. The elucidation of these factors is essential for understanding large channel selectivity and its regulation in vivo.

The anomalous mole fraction effect in calcium channels: a measure of preferential selectivity

The cause of the anomalous mole fraction effect (AMFE) in calcium-selective ion channels is studied. An AMFE occurs when the conductance through a channel is lower in a mixture of salts than in the pure salts at the same concentration. The textbook interpretation of the AMFE is that multiple ions move through the pore in coordinated, single-file motion. Instead of this, we find that at its most basic level an AMFE reflects a channel's preferential binding selectivity for one ion species over another. The AMFE is explained by considering the charged and uncharged regions of the pore as electrical resistors in series: the AMFE is produced by these regions of high and low ion concentration changing differently with mole fraction due to the preferential ion selectivity. This is demonstrated with simulations of a model L-type calcium channel and a mathematical analysis of a simplistic point-charge model. The particle simulations reproduce the experimental data of two L-type channel AMFEs. Conditions under which an AMFE may be found experimentally are discussed. The resistors-in-series model provides a fundamentally different explanation of the AMFE than the traditional theory and does not require single filing, multiple occupancy, or momentum-correlated ion motion.

Bubbles, gating, and anesthetics in ion channels

We suggest that bubbles are the bistable hydrophobic gates responsible for the on-off transitions of single channel currents. In this view, many types of channels gate by the same physical mechanism-dewetting by capillary evaporation-but different types of channels use different sensors to modulate hydrophobic properties of the channel wall and thereby trigger and control bubbles and gating. Spontaneous emptying of channels has been seen in many simulations. Because of the physics involved, such phase transitions are inherently sensitive, unstable threshold phenomena that are difficult to simulate reproducibly and thus convincingly. We present a thermodynamic analysis of a bubble gate using morphometric density functional theory of classical (not quantum) mechanics. Thermodynamic analysis of phase transitions is generally more reproducible and less sensitive to details than simulations. Anesthetic actions of inert gases-and their interactions with hydrostatic pressure (e.g., nitrogen narcosis)-can be easily understood by actions on bubbles. A general theory of gas anesthesia may involve bubbles in channels. Only experiments can show whether, or when, or which channels actually use bubbles as hydrophobic gates: direct observation of bubbles in channels is needed. Existing experiments show thin gas layers on hydrophobic surfaces in water and suggest that bubbles nearly exist in bulk water.

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.

Steric selectivity in Na channels arising from protein polarization and mobile side chains

Monte Carlo simulations of equilibrium selectivity of Na channels with a DEKA locus are performed over a range of radius R and protein dielectric coefficient epsilon(p). Selectivity arises from the balance of electrostatic forces and steric repulsion by excluded volume of ions and side chains of the channel protein in the highly concentrated and charged (approximately 30 M) selectivity filter resembling an ionic liquid. Ions and structural side chains are described as mobile charged hard spheres that assume positions of minimal free energy. Water is a dielectric continuum. Size selectivity (ratio of Na+ occupancy to K+ occupancy) and charge selectivity (Na+ to Ca2+) are computed in concentrations as low as 10(-5) M Ca2+. In general, small R reduces ion occupancy and favors Na+ over K+ because of steric repulsion. Small epsilon(p) increases occupancy and favors Na+ over Ca2+ because protein polarization amplifies the pore's net charge. Size selectivity depends on R and is independent of epsilon(p); charge selectivity depends on both R and epsilon(p). Thus, small R and epsilon(p) make an efficient Na channel that excludes K+ and Ca2+ while maximizing Na+ occupancy. Selectivity properties depend on interactions that cannot be described by qualitative or verbal models or by quantitative models with a fixed free energy landscape.

Kinetic lattice grand canonical Monte Carlo simulation for ion current calculations in a model ion channel system

An algorithm in which kinetic lattice grand canonical Monte Carlo simulations are combined with mean field theory (KLGCMC/MF) is presented to calculate ion currents in a model ion channel system. In this simulation, the relevant region of the system is treated by KLGCMC simulations, while the rest of the system is described by modified Poisson-Boltzmann mean field theory. Calculation of reaction field due to induced charges on the channel/water and membrane/water boundaries is carried out using a basis-set expansion method [Im and Roux, J. Chem. Phys. 115, 4850 (2001)]. Calculation of ion currents, electrostatic potentials, and ion concentrations, as obtained from the KLGCMC/MF simulations, shows good agreement with Poisson-Nernst-Planck (PNP) theory predictions when the channel and membrane have the same dielectric constant as water. If the channel and membrane have a lower dielectric constant than water, however, there is a considerable difference between the KLGCMC/MF and PNP predictions. This difference is attributed to the reaction field, which is missing in PNP theory. It is demonstrated that the reaction field as well as fixed charges in the channel play key roles in selective ion transport. Limitations and further development of the current KLGCMC/MF approach are also discussed.

Dissecting ICRAC, a store-operated calcium current

The use of Ca(2+) for intracellular signalling necessitates tight local and global control of cytoplasmic Ca(2+) concentration, and mechanisms for maintaining the net Ca(2+) balance. It has long been recognized that intracellular Ca(2+) stores exert control over Ca(2+) influx at the plasma membrane through a process of store-operated Ca(2+) entry (SOCE). The Ca(2+) current I(CRAC) is the best characterized instance of SOCE, but the elements of the pathway leading to I(CRAC) have eluded biochemical definition for more than a decade. However, the recent identification of key proteins underlying I(CRAC)--STIM1 and Orai1--has led to several insights into this ER-to-plasma membrane signalling system and to the recognition that it is an ancient and conserved mechanism in multicellular organisms.

Ca2+ selectivity of a chemically modified OmpF with reduced pore volume

We studied an E. coli OmpF mutant (LECE) containing both an EEEE-like locus, typical of Ca(2+) channels, and an accessible and reactive cysteine. After chemical modification with the cysteine-specific, negatively charged (-1e) reagents MTSES or glutathione, this LECE mutant was tested for Ca(2+) versus alkali metal selectivity. Selectivity was measured by conductance and zero-current potential. Conductance measurements showed that glutathione-modified LECE had reduced conductance at Ca(2+) mole fractions <10(-3). MTSES-modified LECE did not. Apparently, the LECE protein is (somehow) a better Ca(2+) chelator after modification with the larger glutathione. Zero-current potential measurements revealed a Ca(2+) versus monovalent cation selectivity that was highest in the presence of Li(+) and lowest in the presence of Cs(+). Our data clearly show that after the binding of Ca(2+) the LECE pore (even with the bulky glutathione present) is spacious enough to allow monovalent cations to pass. Theoretical computations based on density functional theory combined with Poisson-Nernst-Planck theory and a reduced pore model suggest a functional separation of ionic pathways in the pore, one that is specific for small and highly charged ions, and one that accepts preferentially large ions, such as Cs(+).

Orai1 is an essential pore subunit of the CRAC channel

Stimulation of immune cells causes depletion of Ca2+ from endoplasmic reticulum (ER) stores, thereby triggering sustained Ca2+ entry through store-operated Ca2+ release-activated Ca2+ (CRAC) channels, an essential signal for lymphocyte activation and proliferation. Recent evidence indicates that activation of CRAC current is initiated by STIM proteins, which sense ER Ca2+ levels through an EF-hand located in the ER lumen and relocalize upon store depletion into puncta closely associated with the plasma membrane. We and others recently identified Drosophila Orai and human Orai1 (also called TMEM142A) as critical components of store-operated Ca2+ entry downstream of STIM. Combined overexpression of Orai and Stim in Drosophila cells, or Orai1 and STIM1 in mammalian cells, leads to a marked increase in CRAC current. However, these experiments did not establish whether Orai is an essential intracellular link between STIM and the CRAC channel, an accessory protein in the plasma membrane, or an actual pore subunit. Here we show that Orai1 is a plasma membrane protein, and that CRAC channel function is sensitive to mutation of two conserved acidic residues in the transmembrane segments. E106D and E190Q substitutions in transmembrane helices 1 and 3, respectively, diminish Ca2+ influx, increase current carried by monovalent cations, and render the channel permeable to Cs+. These changes in ion selectivity provide strong evidence that Orai1 is a pore subunit of the CRAC channel.

Ion transport through membrane-spanning nanopores studied by molecular dynamics simulations and continuum electrostatics calculations

Narrow hydrophobic regions are a common feature of biological channels, with possible roles in ion-channel gating. We study the principles that govern ion transport through narrow hydrophobic membrane pores by molecular dynamics simulation of model membranes formed of hexagonally packed carbon nanotubes. We focus on the factors that determine the energetics of ion translocation through such nonpolar nanopores and compare the resulting free-energy barriers for pores with different diameters corresponding to the gating regions in closed and open forms of potassium channels. Our model system also allows us to compare the results from molecular dynamics simulations directly to continuum electrostatics calculations. Both simulations and continuum calculations show that subnanometer wide pores pose a huge free-energy barrier for ions, but a small increase in the pore diameter to approximately 1 nm nearly eliminates that barrier. We also find that in those wider channels the ion mobility is comparable to that in the bulk phase. By calculating local electrostatic potentials, we show that the long range Coulomb interactions of ions are strongly screened in the wide water-filled channels. Whereas continuum calculations capture the overall energetics reasonably well, the local water structure, which is not accounted for in this model, leads to interesting effects such as the preference of hydrated ions to move along the pore wall rather than through the center of the pore.

Salting out the ionic selectivity of a wide channel: the asymmetry of OmpF

Although the crystallographic structure of the bacterial porin OmpF has been known for a decade, the physical mechanisms of its ionic selectivity are still under investigation. We address this issue in a series of experiments with varied pH, salt concentrations, inverted salt gradient, and charged and uncharged lipids. Measuring reversal potential, we show that OmpF selectivity (traditionally regarded as slightly cationic) depends strongly on pH and salt concentration and is conditionally asymmetric, that is, the calculated selectivity is sensitive to the direction of salt concentration gradient. At neutral pH and subdecimolar salt concentrations the channel exhibits nearly ideal cation selectivity (t(G)(+)=0.98+/-0.01). Substituting neutral DPhPC with DPhPS, we demonstrate that the fixed charge of the host lipid has a small but measurable effect on the channel reversal potential. The available structural information allows for a qualitative explanation of our experimental findings. These findings now lead us to re-examine the ionization state of 102 titratable sites of the OmpF channel. Using standard methods of continuum electrostatics tailored to our particular purpose, we find the charge distribution in the channel as a function of solution acidity and relate the pH-dependent asymmetry in channel selectivity to the pH-dependent asymmetry in charge distribution. In an attempt to find a simple phenomenological description of our results, we also discuss different macroscopic models of electrodiffusion through large channels.

Saturation of conductance in single ion channels: the blocking effect of the near reaction field

The ionic current flowing through a protein channel in the membrane of a biological cell depends on the concentration of the permeant ion, as well as on many other variables. As the concentration increases, the rate of arrival of bath ions to the channel's entrance increases, and typically so does the net current. This concentration dependence is part of traditional diffusion and rate models that predict Michaelis-Menten current-concentration relations for a single ion channel. Such models, however, neglect other effects of bath concentrations on the net current. The net current depends not only on the entrance rate of ions into the channel, but also on forces acting on ions inside the channel. These forces, in turn, depend not only on the applied potential and charge distribution of the channel, but also on the long-range Coulombic interactions with the surrounding bath ions. In this paper, we study the effects of bath concentrations on the average force on an ion in a single ion channel. We show that the force of the reaction field on a discrete ion inside a channel embedded in an uncharged lipid membrane contains a blocking (shielding) term that is proportional to the square root of the ionic bath concentration. We then show that different blocking strengths yield different behavior of the current-concentration and conductance-concentration curves. Our theory shows that at low concentrations, when the blocking force is weak, conductance grows linearly with concentration, as in traditional models, e.g., Michaelis-Menten formulations. As the concentration increases to a range of moderate shielding, conductance grows as the square root of concentration, whereas at high concentrations, with high shielding, conductance may actually decrease with increasing concentrations: the conductance-concentration curve can invert. Therefore, electrostatic interactions between bath ions and the single ion inside the channel can explain the different regimes of conductance-concentration relations observed in experiments.

Permeation properties of an engineered bacterial OmpF porin containing the EEEE-locus of Ca2+ channels

The selectivity filter of the bacterial porin OmpF carries a small net charge close to -1 e and is therefore only slightly cation-selective. Calcium channels, on the other hand, contain four negatively charged glutamates, the EEEE-locus, and are among the most selective cation channels known. We aimed to turn the essentially nonselective OmpF into a Ca2+-selective channel. To that end, two additional glutamates (R42E and R132E) were introduced in the OmpF constriction zone that already contains D113 and E117. Mutant OmpF containing this DEEE-locus has a high Ca2+ over Cl- selectivity and a Na+ current with a strongly increased sensitivity to 1 mM Ca2+. The charge/space competition model, initially applied to the L-type Ca2+ channel, identifies the fixed charge and filter volume as key determinants of ion selectivity, with the precise atomic arrangement having only second-order effects. By implication, the reproduction of fixed charge and filter volume should transform two channels into channels of similar selectivity, even if the two belong to entirely different ion channel families, as is the case for OmpF and the L-type Ca2+ channel. The results presented here fit quite well in the framework of charge/space competition theory.

Maximum entropy formulation of the Kirkwood superposition approximation

Using a variational formulation, we derive the Kirkwood superposition approximation for systems at equilibrium in the thermodynamic limit. We define the entropy of the triplet correlation function and show that the Kirkwood closure brings the entropy to its maximal value. This approach leads to a different interpretation for the Kirkwood closure relation, usually explained by probabilistic considerations of dependence and independence of particles. The Kirkwood closure is generalized to finite volume systems at equilibrium by computing the pair correlation function in finite domains. Closure relations for high order correlation functions are also found using a variational approach. In particular, maximizing the entropy of quadruplets leads to the high order closure g(1234)=g(123)g(124)g(134)g(234)/[g(12)g(13)g(14)g(23)g(24)g(34)] used in the Born-Green-Yvon 2 equations which are a pair of integral equations for the triplet and pair correlation functions.

Electrodiffusion near an ion channel and the effect of mobile buffer

A numerical method was set up to calculate the dynamic concentration behavior of charged particles in the vicinity of an ion channel. It takes into account the electric potential due to the charge of the transported ions. Additionally, the finite on- and off-kinetics of mobile ion-buffers such as EGTA can be added to the simulations. The calculations were carried out using a modified Crank-Nicolson algorithm to solve the partial differential equations describing the problem. It was found that the electrostatic effect on the concentration of permeating ions is negligible in the presence of physiological salt concentrations. Nevertheless, there are electrostatic effects on other ion species near the channel mouth. Studies on the effect of a Ba2+ -current through a Ca2+ -channel onto the Ca2+ -concentration in the bath, and on the amplification of the Ca2+ -effect on the BK-channel due to the K+ -flux are presented. Additionally, the effect of mobile buffers was simulated and the numerical results are compared with some common analytical approximations.

Density functional theory of charged, hard-sphere fluids

An approximate electrostatic (ES) excess free energy functional for charged, hard sphere fluids is presented. This functional is designed for systems with large density variations, but may also be applied to systems without such variations. Based on the Rosenfeld method of perturbation about a bulk (homogeneous) reference fluid [Y. Rosenfeld, J. Chem. Phys. 98, 8126 (1993)], the new ES functional replaces the reference fluid densities with a functional of the particle densities, called the RFD functional. The first-order direct correlation function (DCF) in the particle densities is computed using as input the first- and second-order DCFs in [rho(i)(x)], the inhomogeneous densities defined by the RFD functional. Because this formulation imposes no a priori constraints on the form of the RFD functional-it is valid for any choice of [rho(i)(x)]-the RFD functional may be chosen (1) so that the input DCFs (that is, DCFs in [rho(i)(x)]) may be approximated and (2) so the combination of [rho(i)(x)] and input DCFs yields a good estimate of the first-order DCF in the particle densities. In this way, the general problem of finding the excess free energy functional has been replaced by the specific problem of choosing a RFD functional. We present a particular RFD functional that, together with bulk formulations for the input DCFs, accurately reproduces the results of Monte Carlo simulations.

Dielectric boundary force and its crucial role in gramicidin

In an electrostatic problem with nonuniform geometry, a charge Q in one region induces surface charges [called dielectric boundary charges (DBC)] at boundaries between different dielectrics. These induced surface charges, in return, exert a force [called dielectric boundary force (DBF)] on the charge Q that induced them. The DBF is often overlooked. It is not present in standard continuum theories of (point) ions in or near membranes and proteins, such as Gouy-Chapman, Debye-Huckel, Poisson-Boltzmann or Poisson-Nernst- Planck. The DBF is important when a charge Q is near dielectric interfaces, for example, when ions permeate through protein channels embedded in biological membranes. In this paper, we define the DBF and calculate it explicitly for a planar dielectric wall and for a tunnel geometry resembling the ionic channel gramicidin. In general, we formulate the DBF in a form useful for continuum theories, namely, as a solution of a partial differential equation with boundary conditions. The DBF plays a crucial role in the permeation of ions through the gramicidin channel. A positive ion in the channel produces a DBF of opposite sign to that of the fixed charge force (FCF) produced by the permanent charge of the gramicidin polypeptide, and so the net force on the positive ion is reduced. A negative ion creates a DBF of the same sign as the FCF and so the net (repulsive) force on the negative ion is increased. Thus, a positive ion can permeate the channel, while a negative ion is excluded from it. In gramicidin, it is this balance between the FCF and DBF that allows only singly charged positive ions to move into and through the channel. The DBF is not directly responsible, however, for selectivity between the alkali metal ions (e.g., Li+, Na+, K+): we prove that the DBF on a mobile spherical ion is independent of the ion's radius.