Ion solvation and structural stability in a sodium channel investigated by molecular dynamics calculations

The stability and ion binding properties of the homo-tetrameric pore domain of a prokaryotic, voltage-gated sodium channel are studied by extensive all-atom molecular dynamics simulations, with the channel protein being embedded in a fully hydrated lipid bilayer. It is found that Na(+) ion presents in a mostly hydrated state inside the wide pore of the selectivity filter of the sodium channel, in sharp contrast to the nearly fully dehydrated state for K(+) ions in potassium channels. Our results also indicate that Na(+) ions make contact with only one or two out of the four polypeptide chains forming the selectivity filter, and surprisingly, the selectivity filter exhibits robust stability for various initial ion configurations even in the absence of ions. These findings are quite different from those in potassium channels. Furthermore, an electric field above 0.5V/nm is suggested to be able to induce Na(+) permeation through the selectivity filter.     

The influence of geometry, surface character, and flexibility on the permeation of ions and water through biological pores

A hydrophobic constriction site can act as an efficient barrier to ion and water permeation if its diameter is less than the diameter of an ion's first hydration shell. This hydrophobic gating mechanism is thought to operate in a number of ion channels, e.g. the nicotinic receptor, bacterial mechanosensitive channels (MscL and MscS) and perhaps in some potassium channels (e.g. KcsA, MthK and KvAP). Simplified pore models allow one to investigate the primary characteristics of a conduction pathway, namely its geometry (shape, pore length, and radius), the chemical character of the pore wall surface, and its local flexibility and surface roughness. Our extended (about 0.1 micros) molecular dynamic simulations show that a short hydrophobic pore is closed to water for radii smaller than 0.45 nm. By increasing the polarity of the pore wall (and thus reducing its hydrophobicity) the transition radius can be decreased until for hydrophilic pores liquid water is stable down to a radius comparable to a water molecule's radius. Ions behave similarly but the transition from conducting to non-conducting pores is even steeper and occurs at a radius of 0.65 nm for hydrophobic pores. The presence of water vapour in a constriction zone indicates a barrier for ion permeation. A thermodynamic model can explain the behaviour of water in nanopores in terms of the surface tensions, which leads to a simple measure of 'hydrophobicity' in this context. Furthermore, increased local flexibility decreases the permeability of polar species. An increase in temperature has the same effect, and we hypothesize that both effects can be explained by a decrease in the effective solvent-surface attraction which in turn leads to an increase in the solvent-wall surface free energy.     

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.     

Molecular dynamics study of the KcsA potassium channel

The structural, dynamical, and thermodynamic properties of a model potassium channel are studied using molecular dynamics simulations. We use the recently unveiled protein structure for the KcsA potassium channel from Streptomyces lividans. Total and free energy profiles of potassium and sodium ions reveal a considerable preference for the larger potassium ions. The selectivity of the channel arises from its ability to completely solvate the potassium ions, but not the smaller sodium ions. Self-diffusion of water within the narrow selectivity filter is found to be reduced by an order of magnitude from bulk levels, whereas the wider hydrophobic section of the pore maintains near-bulk self-diffusion. Simulations examining multiple ion configurations suggest a two-ion channel. Ion diffusion is found to be reduced to approximately (1)/(3) of bulk diffusion within the selectivity filter. The reduced ion mobility does not hinder the passage of ions, as permeation appears to be driven by Coulomb repulsion within this multiple ion channel.     

Ion leakage through transient water pores in protein-free lipid membranes driven by transmembrane ionic charge imbalance

We have employed atomic-scale molecular dynamics simulations to address ion leakage through transient water pores in protein-free phospholipid membranes. Our results for phospholipid membranes in aqueous solution with NaCl and KCl salts show that the formation of transient water pores and the consequent ion leakage can be induced and be driven by a transmembrane ionic charge imbalance, an inherent feature in living cells. These processes take place if the gradient is large enough to develop a sufficiently significant potential difference across the membrane. The transport of cations and anions through the water pores is then seen; it discharges the transmembrane potential, considerably reduces the size of a water pore, and makes the water pore metastable, leading eventually to its sealing. The ion transport is found to be sensitive to the type of ions. It turns out that Na(+) and Cl(-) ions leak through a membrane at approximately the same ratio despite the fact that Na(+) ions are expected to experience a lower potential barrier for the permeation through the pore. This is because of strong interactions of sodium ions with the carbonyl region of a phospholipid membrane as well as with lipid headgroups forming pore "walls," considerably slowing down the permeation of sodium ions. In contrast, we observed a pronounced selectivity of a phospholipid membrane to the permeation of potassium ions as compared to chloride ions: Potassium ions, being larger than sodium ions, interact only weakly with phospholipid headgroups, so that these interactions are not able to compensate for a large difference in free-energy barriers for permeation of K(+) and Cl(-) ions. These findings are found to be robust to a choice of force-field parameters for ions (tested by Gromacs and Charmm force-fields for ions). What is more, a potassium ion is found to be able to permeate a membrane along an alternate, "water-defect-mediated" pathway without actual formation of a pore. The "water-defect-mediated" leakage involves formation of a single water defect only and is found to be at least one order of magnitude faster than the pore-mediated ion leakage.     

Energy barriers for passage of ions through channels. Exact solution of two electrostatic problems

Exact solutions are given to two electrostatic problems relevant to ion permeation through pores in membranes. The first assesses the importance of the pore forming molecule as a dielectric shield. It is shown on the basis of structural and dielectric considerations alone (neglecting effects attributable to possible charge distribution at the interior surface of the pre-former) that the minimum electrostatic barrier for monovalent ion passage through a gramicidin-like channel is 11 kT. It is further shown that given favorable circumstances, dielectric shielding might dramatically reduce the barrier to ion passage through potassium channels. The second problem considers the error introduced by treating ions as point charges. It is shown that for structureless pores the point charge approximation introduces no meaningful error, even if the ratio of ion radius to pore radius is as great as 0.95.     

Ion selectivity in potassium channels

Potassium channels are tetrameric membrane-spanning proteins that provide a selective pore for the conduction of K(+) across the cell membranes. One of the main physiological functions of potassium channels is efficient and very selective transport of K(+) ions through the membrane to the cell. Classical views of ion selectivity are summarized within a historical perspective, and contrasted with the molecular dynamics (MD) simulations free energy perturbation (FEP) performed on the basis of the crystallographic structure of the KcsA phospholipid membrane. The results show that the KcsA channel does not select for K(+) ions by providing a binding site of an appropriate (fixed) cavity size. Rather, selectivity for K(+) arises directly from the intrinsic local physical properties of the ligands coordinating the cation in the binding site, and is a robust feature of a pore symmetrically lined by backbone carbonyl groups. Further analysis reveals that it is the interplay between the attractive ion-ligand (favoring smaller cation) and repulsive ligand-ligand interactions (favoring larger cations) that is the basic element governing Na(+)/K(+) selectivity in flexible protein binding sites. Because the number and the type of ligands coordinating an ion directly modulate such local interactions, this provides a potent molecular mechanism to achieve and maintain a high selectivity in protein binding sites despite a significant conformational flexibility.     

Dynamic hydration numbers for biologically important ions

The role of ionized groups in biological systems is determined by their affinity for water [Biophys. J. 72 (1997) 65-76]. The tightly bound water associated with biologically important ions increases their apparent size. We define the apparent dynamic hydration number of an ion here as the number of tightly bound water molecules that must be assigned to the ion to explain its apparent molecular weight on a Sephadex G-10 size exclusion column, and report the first accurate determination of tightly bound water for 23 ions of biological significance, including H(+) and HO(-). We also calculate the radius of the equivalent hydrated sphere (r(h)) for each ion. We find that the ratio of the hydrated volumes of two ions approximates the ratio of the square of the charges of the same two ions. Since the 'ionic strength' of the solution also depends upon the square of the charges on the ions, our results suggest that ionic strength effects may largely arise from local effects related to the hydrated volume of the ion--that is, from space filling, osmotic, water activity, surface tension and hydration shell overlap effects rather than from long-range electric field effects.     

Importance of hydration and dynamics on the selectivity of the KcsA and NaK channels

Fundamental concepts governing ion selectivity in narrow pores are reviewed and the microscopic factors responsible for the lack of selectivity of the NaK channel, which is structurally similar to the K+-selective KcsA channel, are elucidated on the basis of all-atom molecular dynamics free energy simulations. The results on NaK are contrasted and compared with previous studies of the KcsA channel. Analysis indicates that differences in hydration of the cation in the pore of NaK is at the origin of the lack of selectivity of NaK.     

Microscopic model for selective permeation in ion channels.

Ionic permeation in the selectivity filter of ion channels is analyzed by a microscopic model based on molecular kinetic theory. The energy and flux equations are derived by assuming that: (a) the selectivity filter is formed by a symmetrical array of carbonyl groups; (b) ion movement is near the axis of the channel; (c) a fraction of water molecules is separated from the ion while it moves across the selectivity filter; (d) the applied voltage drops linearly across the selectivity filter; (e) ions move independently. Energy profiles, single channel conductances, and the degree of hydration of K+ in a hypothetical K+ channel are examined by varying the following microscopic parameters: ion radius and mass, channel radius, number of effective water dipoles, and number of carbonyl groups. The i-V curve is linear up to +/- 170 mV. If the positions of energy maxima and minima are fixed, this linear range is reduced to +/- 50 mV. Channel radius and ion-water interactions are found to be two major channel structural determinants for selectivity sequences. Both radius and mass of an ion are important in selectivity mediated by these interactions. The theory predicts a total of 15 possible kinetic selectivity sequences for alkali cations in ion channels with a single selectivity filter.     

Exploring the open pore of the potassium channel from Streptomyces lividans

The tetrameric potassium channel from Streptomyces lividans (KcsA) embedded in planar bilayers exhibits the following electrophysiological characteristics: (i) K+ ions can cross the pore in a highly hydrated state (nH2O > or = 6), (ii) the selectivity for K+ exceeds that for Na+ ions by 11 times, and both Ca2+ and Mg2+ are permeant, (iii) the internal side is blocked by Ba2+ ions in a voltage-dependent manner, (iv) intrinsic rectification is due to gating, depending on the direction of the electric field, (v) the internal side is pH-sensitive, and (vi) the open pore has a diameter of approximately 5.8 A. In conclusion, our results show that ion conduction and selectivity of KcsA cannot easily be reconciled with the properties deduced from the rigid crystal structure [Doyle et al., Science 280 (1998) 69-77], which must be concluded to have the pore trapped in its closed state.     

Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations

We have performed simulations of both a single potassium ion and a single sodium ion within the pore of the bacterial potassium channel KcsA. For both ions there is a dehydration energy barrier at the cytoplasmic mouth suggesting that the crystal structure is a closed conformation of the channel. There is a potential energy barrier for a sodium ion in the selectivity filter that is not seen for potassium. Radial distribution functions for both ions with the carbonyl oxygens of the selectivity filter indicate that sodium may interact more tightly with the filter than does potassium. This suggests that the key to the ion selectivity of KcsA is the greater dehydration energy of Na⁺ ions, and helps to explain the block of KcsA by internal Na⁺ ions.     

Hydrophobic Gating of Ion Permeation in Magnesium Channel CorA

Ion channels catalyze ionic permeation across membranes via water-filled pores. To understand how changes in intracellular magnesium concentration regulate the influx of Mg²⁺ into cells, we examine early events in the relaxation of Mg²⁺ channel CorA toward its open state using massively-repeated molecular dynamics simulations conducted either with or without regulatory ions. The pore of CorA contains a 2-nm-long hydrophobic bottleneck which remained dehydrated in most simulations. However, rapid hydration or “wetting” events concurrent with small-amplitude fluctuations in pore diameter occurred spontaneously and reversibly. In the absence of regulatory ions, wetting transitions are more likely and include a wet state that is significantly more stable and more hydrated. The free energy profile for Mg²⁺ permeation presents a barrier whose magnitude is anticorrelated to pore diameter and the extent of hydrophobic hydration. These findings support an allosteric mechanism whereby wetting of a hydrophobic gate couples changes in intracellular magnesium concentration to the onset of ionic conduction.     

Ion Concentration- and Voltage-Dependent Push and Pull Mechanisms of Potassium Channel Ion Conduction

The mechanism of ion conduction by potassium channels is one of the central issues in physiology. In particular, it is still unclear how the ion concentration and the membrane voltage drive ion conduction. We have investigated the dynamics of the ion conduction processes in the Kv1.2 pore domain, by molecular dynamics (MD) simulations with several different voltages and ion concentrations. By focusing on the detailed ion movements through the pore including selectivity filter (SF) and cavity, we found two major conduction mechanisms, called the III-IV-III and III-II-III mechanisms, and the balance between the ion concentration and the voltage determines the mechanism preference. In the III-IV-III mechanism, the outermost ion in the pore is pushed out by a new ion coming from the intracellular fluid, and four-ion states were transiently observed. In the III-II-III mechanism, the outermost ion is pulled out first, without pushing by incoming ions. Increases in the ion concentration and voltage accelerated ion conductions, but their mechanisms were different. The increase in the ion concentrations facilitated the III-IV-III conductions, while the higher voltages increased the III-II-III conductions, indicating that the pore domain of potassium channels permeates ions by using two different driving forces: a push by intracellular ions and a pull by voltage.     

Ion binding properties and structure stability of the NaK channel

Ion distribution in the selectivity filter and ion-water and ion-protein interactions of NaK channel are systematically investigated by all-atom molecular dynamics simulations, with the tetramer channel protein being embedded in a solvated phospholipid bilayer. Analysis of the simulation results indicates that K⁺ ions prefer to bind within the sites formed by two adjacent planes of oxygen atoms from the selectivity filter, while Na⁺ ions are inclined to bind to a single plane of four oxygen atoms. At the same time, both K⁺ and Na⁺ ions can diffuse in the vestibule, accompanying with movements of the water molecules confined in a complex formed by the vestibule together with four small grottos connecting to it. As a result, K⁺ ions show a wide range of coordination numbers (6-8), while Na⁺ ions display a constant coordination number of ∼ 6 in the selectivity filter, which may result in the loss of selectivity of NaK. It is also found that a Ca²⁺ can bind at the extracellular site as reported in the crystal structure in a partially hydrated state, or at a higher site in a full hydration state. Furthermore, the carbonyl group of Asp66 can reorient to point towards the center pore when an ion exists in the vestibule, while that of Gly65 always aligns tangentially to the channel axis, as in the crystallographic structures.     

Characterization of aqueous pores in plant cuticles and permeation of ionic solutes

Plant cuticles are lipid membranes with separate diffusion paths for lipophilic non-electrolytes and hydrated ionic compounds. Ions are lipid insoluble and require an aqueous pathway across cuticles. Based on experimental data, the aqueous pathway in cuticles has been characterized. Aqueous pores arise by hydration of permanent dipoles and ionic functional groups. They can be localized using ionic fluorescent dyes, silver nitrate, and mercuric chloride. Aqueous pores preferentially occur in cuticular ledges, at the base of trichomes, and in cuticles over anticlinal walls. Average pore radii ranged from 0.45 to 1.18 nm. Penetration of ions was a first order process as the fraction of the salt remaining on the cuticle surface decreased exponentially with time. Permeability of cuticles to ions depended on humidity and was highest at 100% humidity. Wetting agents increased rate constants by factors of up to 12, which indicates that the pore openings are surrounded by waxes. The pores in cuticular ledges of Helxine soleirolii allowed passage of berberine sulphate, which has a molecular weight of 769 g mol(-1). Increasing the molecular weight of solutes from 100 to 500 g mol(-1) decreased the rate constants of penetration by factors of 7 (Vicia faba) and 13 (Populus canescens), respectively. Half-times of penetration of inorganic salts and organic ions across Populus cuticles and Vicia leaf surfaces varied between 1 and 12 h. This shows that penetration of ionic compounds can be fairly rapid, and ions with molecular weights of up to 800 g mol(-1) can penetrate cuticles that possess aqueous pores.     

Further analysis of counterion permeation through anion-selective glycine receptor channels

The functional role of ion channels, which allow counterion permeation, depends critically on their relative anion-cation relative selectivity. From whole-cell patch clamp reversal potential measurements under dilution potential conditions, we have already shown that anion-cation permeabilities of anion-selective wild-type (WT) and mutant (with larger pore diameter) glycine receptor (GlyR) channels in the presence of Li+, Na+ and Cs+ counterions, were inversely correlated with the equivalent hydration diameter of the counterion, with chloride-cation permeability increasing as counterion equivalent hydration diameter increased with respect to the channel minimum pore diameter. Corrected for liquid junction potentials (LJPs; using ion activities), the previous chloride-cation permeabilities for the alkali cations were 23.4 (Li+), 10.9 (Na+) and 5.0 (Cs+) for the smaller WT channel. Further analysis to incorporate an initial offset potential correction, to fully allow for slight differences between internal cell composition and external control salt solution, changed the above permeability ratios to 30.6 (Li+), 11.8 (Na+) and 5.0 (Cs+), adding enhanced support for the inverse correlation between anion-to-counterion permeability ratio and equivalent hydrated counterion diameter relative to channel pore diameter (erroneously ignoring LJPs reduces each permeability ratio to about 4). Also, new direct measurements of LJPs (for NaCl and LiCl salt dilutions) using a 3M KCl-agar reference salt bridge (with freshly-cut end for each solution composition change) have shown excellent agreement with calculated LJPs (using ion activities), validating calculated LJP values. We continue to suggest that counterion cations permeate with chloride ions as neutral pairs.     

A thermodynamic approach to alamethicin pore formation

The structure and energetics of alamethicin Rf30 monomer to nonamer in cylindrical pores of 5 to 11 Å radius are investigated using molecular dynamics simulations in an implicit membrane model that includes the free energy cost of acyl chain hydrophobic area exposure. Stable, low energy pores are obtained for certain combinations of radius and oligomeric number. The trimer and the tetramer formed 6 Å pores that appear closed while the larger oligomers formed open pores at their optimal radius. The hexamer in an 8 Å pore and the octamer in an 11 Å pore give the lowest effective energy per monomer. However, all oligomers beyond the pentamer have comparable energies, consistent with the observation of multiple conductance levels. The results are consistent with the widely accepted “barrel-stave” model. The N terminal portion of the molecule exhibits smaller tilt with respect to the membrane normal than the C terminal portion, resulting in a pore shape that is a hybrid between a funnel and an hourglass. Transmembrane voltage has little effect on the structure of the oligomers but enhances or decreases their stability depending on its orientation. Antiparallel bundles are lower in energy than the commonly accepted parallel ones and could be present under certain experimental conditions. Dry aggregates (without an aqueous pore) have lower average effective energy than the corresponding aggregates in a pore, suggesting that alamethicin pores may be excited states that are stabilized in part by voltage and in part by the ion flow itself.     

Ion behavior in the selectivity filter of HCN1 channels

Hyperpolarization-activated cyclic-nucleotide gated channels (HCNs) are responsible for the generation of pacemaker currents (I_f or I_h) in cardiac and neuronal cells. Despite the overall structural similarity to voltage-gated potassium (Kv) channels, HCNs show much lower selectivity for K⁺ over Na⁺ ions. This increased permeability to Na⁺ is critical to their role in membrane depolarization. HCNs can also select between Na⁺ and Li⁺ ions. Here, we investigate the unique ion selectivity properties of HCNs using molecular-dynamics simulations. Our simulations suggest that the HCN1 pore is flexible and dilated compared with Kv channels with only one stable ion binding site within the selectivity filter. We also observe that ion coordination and hydration differ within the HCN1 selectivity filter compared with those in Kv and cyclic-nucleotide gated channels. Additionally, the C358T mutation further stabilizes the symmetry of the binding site and provides a more fit space for ion coordination, particularly for Li⁺.