Not only enthalpy: large entropy contribution to ion permeation barriers in single-file channels

The effect of channel length on the barrier for potassium ion permeation through single-file channels has been studied by means of all-atom molecular dynamics simulations. Using series of peptidic gramicidin-like and simplified ring-structured channels, both embedded in model membranes, we obtained two distinct types of behavior: saturation of the central free energy barriers for peptidic channels and a linear increase in simplified ring-structured channels with increasing channel length. The saturation of the central free energy barrier for the peptidic channels occurs at relatively short lengths, and it is correlated with the desolvation from the bulk water. Remarkably, decomposition of free energy barriers into enthalpic and entropic terms reveals an entropic cost for ion permeation. Furthermore, this entropic cost dominates the ion permeation free energy barrier, since the corresponding free energy contribution is higher than the enthalpic barrier. We conclude that the length dependence of the free energy is enthalpy-dominated, but the entropy is the major contribution to the permeation barrier. The decrease in rotational water motion and the reduction of channel mobility are putative origins for the overall entropic penalty.     

Properties of amphotericin B channels in a lipid bilayer

Properties of individual ionic channels formed by polyene antibiotic Amphotericin B were studied on brain phospholipid membranes containing cholesterol. The ionic channels have a closed state and an open one (with conductance of about 6.5 pS in 2 M KCl). The conductance value of an open channel is independent of cholesterol concentration in the membrane and of pH in the range from 3.5 to 8.0. The voltage-current characteristics of a single channel are superlinear. Zero current potential value in the case of different KCl concentrations in the two solutions indicates preferential but not ideal anionic selectivity of a single channel. Channel conductivity grows as the electrolyte concentration is increased and tends to a limiting value at high concentrations. A simple model having only one site for an ion was shown to represent satisfactorily an open channel behaviour under different conditions. An individual ionic channel performs a large number of transitions between the open and closed states during its life-time of several minutes. Rate constants of these transitions depend on the kind and concentration of salt in aqueous solutions. The switching system functioning is not influenced by an ion situated inside the pore.     

Kinetic characteristics of the excitability-inducing material channel in oxidized cholesterol and brain lipid bilayer membranes

The kinetic characteristics of the opening and closing of the excitability-inducing material (EIM) channel in oxidized cholesterol and in brain lipid bilayers are compared. The kinetics of the opening and closing of individual ion-conducting channels in bilayers doped with small amounts of EIM are determined from discrete fluctuations in ionic current. The kinetics for approach to steady-state conductance are determined for lipid bilayers containing many channels. Steady-state and kinetic characteristics for the EIM channel incorporated in brain lipid bilayers can be accounted for by the model developed for the EIM channel incorporated in oxidized cholesterol membranes. Relaxation time, calculated from rate constants of single-channel membranes or directly measured in many-channel membranes is strongly temperature dependent, and is always shorter in brain lipid membranes. Changes in temperature do not affect the interaction of the electric field and the open channel, but the open configuration of the EIM channel in brain lipid bilayers is stablized with increasing temperature. The configurational energy difference between the open and closed channel, calculated from temperature studies, is larger in brain lipid bilayers. The energy barrier which separates the two configurations of the channel is larger in oxidized cholesterol bilayers.     

Statistical properties of ion channel records. Part II: estimation from the macroscopic current

Macroscopic ion channel current is the summation of the stochastic records of individual channel currents and therefore relates to their statistical properties. As a consequence of this relationship, it may be possible to derive certain statistical properties of single channel records or even generate some estimates of the records themselves from the macroscopic current when the direct measurement of single channel currents is not applicable. We present a procedure for generating the single channel records of an ion channel from its macroscopic current when the stochastic process of channel gating has the following two properties: (I) the open duration is independent of the time of opening event and has a single exponential probability density function (pdf), (II) all the channels have the same probability to open at time t. The application of this procedure is considered for cases where direct measurement of single channel records is difficult or impossible. First, the probability density function (pdf) of opening events, a statistical property of single channel records, is derived from the normalized macroscopic current and mean channel open duration. Second, it is shown that under the conditions (I) and (II), a non-stationary Markov model can represent the stochastic process of channel gating. Third, the non-stationary Markov model is calibrated using the results of the first step. The non-stationary formulation increases the model ability to generate a variety of different single channel records compared to common stationary Markov models. The model is then used to generate single channel records and to obtain other statistical properties of the records. Experimental single channel records of inactivating BK potassium channels are used to evaluate how accurately this procedure reconstructs measured single channel sweeps.     

Gramicidin channel selectivity. Molecular mechanics calculations for formamidinium, guanidinium, and acetamidinium.

Empirical energy function calculations were used to evaluate the effects of minimization on the structure of a gramicidin A channel and to analyze the energies of interaction between three cations (guanidinium, acetamidinium, formamidinium) and the channel as a function of position along the channel axis. The energy minimized model of the gramicidin channel, which was based on the results of Venkatachalam and Urry (1983), has a constriction at the channel entrance. If the channel is not allowed to relax in the presence of the ions (rigid model), there is a large potential energy barrier for all three cations. The barrier varies with cation size and is due to high van der Waals and ion deformation energies. If the channel is minimized in the presence of the ions, the potential energy barrier to formamidinium entry is almost eliminated, but a residual barrier remains for guanidinium and acetamidinium. The residual barrier is primarily due, not to the expansion of the helix, but, to the disruption of hydrogen bonds between the terminal ethanoloamine and the next turn of the helix which occurs when the carbonyls of the outer turn of the helix librate inward toward the ion as it enters the channel. The residual potential energy barriers could be a possible explanation for the measured selectivity of gramicidin for formamidinium over guanidinium. The results of this full-atomic model address the applicability of the size-exclusion concept for the selectivity of the gramicidin channel.     

Is arginine charged in a membrane?

“Charged” amino acids play countless important roles in protein structure and function. Yet when these side chains come into contact with membranes we do not fully understand their behavior. This is highlighted by a recent model of voltage-gated ion channel activity and translocon-based experiments that suggest small penalties to expose these side chains to lipids, opposing the prevailing view in membrane biophysics. Here we employ a side chain analog as well as a transmembrane helix model to determine the free energy as a function of protonation state and position for a lipid-exposed arginine (Arg) residue across a membrane. We observe high free energy barriers for both the charged and neutral states. Due to the stabilizing influence of membrane deformations for the protonated form, the Arg side chain experiences a pK_a shift of ≤4.5 units and remains mostly protonated. The cost for exposing Arg to lipid hydrocarbon is prohibitively high with implications for many membrane translocating processes and the activation mechanisms of voltage-gated ion channels.     

Conversion of a porin-like peptide channel into a gramicidin-like channel by glycine to D-alanine substitutions

The beta-barrel and beta-helix formation, as in porins and gramicidin, respectively, represent two distinct mechanisms for ion channel formation by beta-sheet proteins in membranes. The design of beta-barrel proteins is difficult due to incomplete understanding of the basic principles of folding. The design of gramicidin-like beta-helix relies on an alternating pattern of L- and D-amino acid sequences. Recently we noticed that a short beta-sheet peptide (xSxG)(6), can form porin-like channels via self-association in membranes. Here, we proposed that glycine to D-alanine substitutions of the N-formyl-(xSxG)(6) would transform the porin-like channel into a gramicidin-like beta(12)-helical channel. The requirement of an N-formyl group for channel activity, impermeability to cations with a diameter >4 A, high monovalent cation selectivity, and the absence of either voltage gating or subconductance states upon D-alanine substitution support the idea of a gramicidin-like channel. Moreover, the circular dichroism spectrum in membranes is different, indicating a change in regular beta-sheet backbone structure. The conversion of a complex porin-like channel into a gramicidin-like channel provides a link between two different mechanisms of beta-sheet channel formation in membranes and emphasizes the importance of glycine and D-amino acid residues in protein folding and function and in the engineering of ion channels.     

Ionic selectivity, saturation and block in gramicidin A channels: I. Theory for the electrical properties of ion selective channels having two pairs of binding sites and multiple conductance states.

A model for the gramicidin A channel is proposed which extends existing models by adding a specific cationic binding site at each entrance to the channel. The binding of ions to these outer channel sites is assumed to shift the energy levels of the inner sites and barriers and thereby alter the channel conductance. The resulting properties are analyzed theoretically for the simplest case of two inner sites and a single energy barrier. This for-site model (two outer and two inner) predicts that the membrane potential at zero current (Uo) should be a Goldman-Hodgkin-Katz equation with concentration-dependent permeability ratios. The coefficients of the concentration-dependent terms are shown to be related to the peak energy shifts of the barrier and to the binding constants of the outer sites. The thory also predicts the channel conductance in symmetrical solutions to exhibit three limiting behaviors, from which the properties of the outer and inner sites can be characterized. In two-cation symmetrical mixtures the conductance as a function of mole fraction is shown to have a minimum, and the related phenomenon of inhibition and block exerted by one ion on the other is explained explicitly by the theory. These various phenomena, having ion interactions in a multiply occupied channel as a common physical basis, are all related (by the theory) through a set of measurable parameters describing the properties of the system.     

Ion selectivity mechanism in a bacterial pentameric ligand-gated ion channel

The proton-gated ion channel from Gloeobacter violaceus (GLIC) is a prokaryotic homolog of the eukaryotic nicotinic acetylcholine receptor that responds to the binding of neurotransmitter acetylcholine and mediates fast signal transmission. Recent emergence of a high-resolution crystal structure of GLIC captured in a potentially open state allowed detailed, atomic-level insight into ion conduction and selectivity mechanisms in these channels. Herein, we have examined the barriers to ion conduction and origins of ion selectivity in the GLIC channel by the construction of potential-of-mean-force profiles for sodium and chloride ions inside the transmembrane region. Our calculations reveal that the GLIC channel is open for a sodium ion to transport, but presents a ∼11 kcal/mol free energy barrier for a chloride ion. Our collective findings identify three distinct contributions to the observed preference for the permeant ions. First, there is a substantial contribution due to a ring of negatively charged glutamate residues (E-2′) at the narrow intracellular end of the channel. The negative electrostatics of this region and the ability of the glutamate side chains to directly bind cations would strongly favor the passage of sodium ions while hindering translocation of chloride ions. Second, our results imply a significant hydrophobic contribution to selectivity linked to differences in the desolvation penalty for the sodium versus chloride ions in the central hydrophobic region of the pore. This hydrophobic contribution is evidenced by the large free energy barriers experienced by Cl⁻ in the middle of the pore for both GLIC and the E-2′A mutant. Finally, there is a distinct contribution arising from the overall negative electrostatics of the channel.     

Free energy of a potassium ion in a model of the channel formed by an amphipathic leucine-serine peptide

We use molecular dynamics simulations to investigate the position-dependent free energy of a potassium ion in a model of an ion channel formed by the synthetic amphipathic leucine-serine peptide, LS3. The channel model is a parallel bundle of six LS3 helices around which are packed 146 methane-like spheres in order to mimic a membrane. At either end of and within the channel are 1051 water molecules, plus four ions (two potassium and two chloride). The free energy of a potassium ion in the channel was estimated using the weighted histogram analysis (WHAM) method. This is the first time to our knowledge that such a calculation has been carried out as a function of the position of an ion in three dimensions within a channel. The results indicate that for this channel, which is lined by hydrophilic serine sidechains, there is a relatively weak dependence of the free energy on the axial/off-axial position of the ion. There are some off-axis local minima, especially in the C-terminal half of the channel. Using the free energy results, a single channel current-voltage curve was estimated using a one-dimensional Nernst-Planck equation. Although reasonable agreement with experiment is achieved for K(+) ions flowing from the N-terminal to the C-terminal mouth, in the opposite direction the current is underestimated. This underestimation may be a consequence of under-sampling of the conformational dynamics of the channel. We suggest that our simulations may have captured, for example, a sub-conductance level (i.e. an incompletely open state) of the LS3 channel.     

Investigating the dielectric effects of channel pore water on the electrostatic barriers of the permeation ion by the finite difference Poisson-Boltzmann method

In this paper, the finite difference Poisson-Boltzmann (FDPB) method with four dielectric constants is developed to study the effect of dielectric saturation on the electrostatic barriers of the permeation ion. In this method, the inner shape of the channel pore is explicitly represented, and the fact that the dielectric constant inside the channel pore is different from that of bulk water is taken into account. A model channel system which is a right-handed twist bundle with four α-helical segments is provided for this study. From the FDPB calculations, it is found that the difference of the ionic electrostatic solvation energy for wider domains depends strongly on the pore radius in the vicinity of the ion when the pore dielectric constant is changed from 78 to 5. However, the electrostatic solvation energy of the permeation ion can not be significantly affected by the dielectric constant in regions with small pore radii. Our results indicate that the local electrostatic interactions inside the ion channel are of major importance for ion electrostatic solvation energies, and the effect of dielectric saturation on the electrostatic barriers is coupled to the interior channel dimensions. Received: 28 January 1997 / Accepted: 24 September 1997     

Collective Equilibrium Behaviour of Ion Channel Gating in Cell Membranes: An Ising Model Formulation

A statistical mechanical model for voltage-gated ion channels in cell membranes is proposed using the transfer matrix method. Equilibrium behavior of the system is studied. Representing the distribution of channels over the cellular membrane on a one-dimensional array with each channel having two states (open and closed) and incorporating channel–channel cooperative interactions, we calculate the fraction of channels in the open state at equilibrium. Experimental data obtained from batrachotoxin-modified sodium channels in the squid giant axon, using the cut-open axon technique, is best fit by the model when there is no interaction between the channels.     

Electrical properties of ionic channels formed by Helix pomatia hemocyanin in planar lipid bilayers

Helix pomatia hemocyanin forms ion-conducting channels in planar lipid bilayer membranes when added at mg/ml concentration. These channels have several original features. They fluctuate between one conducting and some poorly conducting states and fluctuations can be grouped in bursts. Different channels can have widely different conductance amplitudes. Both channel conductance and burst lifetime are dependent on the applied voltage. Fluctuations within a burst show a complex kinetic behaviour which has been explained developing a multistate model. The model calls for one single open state and six different closed states. Transitions are allowed only between one of the closed states and the open one and obey first order kinetics. This model is able to fit all our experimental curves obtained in single channel experiments.     

Analysis of fractal ion channel gating kinetics: kinetic rates, energy levels, and activation energies

Ion channels in the cell membranes of the corneal endothelium, hippocampal neurons, and fibroblasts, and gramicidin channels in lipid bilayers have open and closed times that can be fit, in whole or part, by power law distributions. That is, the gating is self-similar when viewed at different time scales. Hence, kinetic processes at slow and fast time scales are not independent but rather are interrelated. To study how such a relationship can arise we analyze a closed in equilibrium open channel with the fractal dimension for leaving the closed state DCO approximately 2 and the fractal dimension for leaving the open state DOC approximately 1. This special case can be analyzed because it can be represented by equivalent Markov processes. We show that it is equivalent to Markov chains with forward and backward kinetic rate constants approximately equal at each stage, and forming an approximate geometric progression along the different stages. These kinetic rates determine the energy levels and activation energy barriers separating those levels. We find that there are many conformational states (substates) separated by high activation energy barriers. This is similar to the energy structure found for globular proteins such as myoglobin. However, the novel feature reported here is that the activation energy barriers are not independent but are interrelated and form an arithmetic progression. Because of this relationship the fast processes across the low activation energy barriers are linked to slow processes across the high activation energy barriers.     

Simple model of ion transport through alamethicin channels in lipid membranes

The electrical characteristics of wide membrane channels such as those induced in lipid membranes by alamethicin have been analyzed using an electrodiffusion model. The channel is considered to be a water filled cylinder in which the potential energy barrier is a result of the difference in polarization energy of the ion environment when the ion is located inside as compared to outside of the channel. In addition, an electric field related to the channel structure is assumed. It is shown that without postulating any specific chemical ion-channel interaction one can reproduce experimental membrane potentials for NaCl, KCl, and CaCl₂ concentration gradients with a single set of channel parameters. The calculations also yield experimental J-V characteristics of discrete conduction states. In addition, a simple mechanism of interchannel coupling based on the above model is discussed. The model suggests a unifying approach to the problem of the origin of interionic selectivity of membrane channels induced by polyene antibiotics.     

Molecular dynamics simulation of water permeation through the alpha-hemolysin channel

The alpha-hemolysin (AHL) nanochannel is a non-selective channel that allows for uncontrolled transport of small molecules across membranes leading to cell death. Although it is a bacterial toxin, it has promising applications, ranging from drug delivery systems to nano-sensing devices. This study focuses on the transport of water molecules through an AHL nanochannel using molecular dynamics (MD) simulations. Our results show that AHL can quickly transport water across membranes. The first-passage time approach was used to estimate the diffusion coefficient and the mean exit time. To study the energetics of transport, the potential of mean force (PMF) of a water molecule along the AHL nanochannel was calculated. The results show that the energy barriers of water permeation across a nanopore are always positive along the channel and the values are close to thermal energy (k_BT). These findings suggest that the observed quick permeation of water is due to small energy barriers and a hydrophobic inner channel surface resulting in smaller friction. We speculate that these physical mechanisms are important in how AHL causes cell death.     

A semi-microscopic Monte Carlo study of permeation energetics in a gramicidin-like channel: the origin of cation selectivity.

The influence of a gramicidin-like channel former on ion free energy barriers is studied using Monte Carlo simulation. The model explicitly describes the ion, the water dipoles, and the peptide carbonyls; the remaining degrees of freedom, bulk electrolyte, non-polar lipid and peptide regions, and electronic (high frequency) permittivity, are treated in continuum terms. Contributions of the channel waters and peptide COs are studied both separately and collectively. We found that if constrained to their original orientations, the COs substantially increase the cationic permeation free energy; with or without water present, CO reorientation is crucial for ion-CO interaction to lower cation free energy barriers; the translocation free energy profiles for potassium-, rubidium-, and cesium-like cations exhibit no broad barriers; the lipid-bound peptide interacts more effectively with anions than cations; anionic translocation free energy profiles exhibit well defined maxima. Using experimental data to estimate transfer free energies of ions and water from bulk electrolyte to a non-polar dielectric (continuum lipid), we found reasonable ion permeation profiles; cations bind and permeate, whereas anions cannot enter the channel. Cation selectivity arises because, for ions of the same size and charge, anions bind hydration water more strongly.     

Model of multifractal gating of single ionic channels in biological membranes

A physico-mathematical model of the gating machinery of single ionic channels in biological membranes has been developed. In the paradigm of this model, gating particles are subjected to: (i) deterministic friction force responsible for interactions of gating particles with the surrounding solution; (ii) deterministic potential force depending on the structure and conformational state of the channel pore (the latter is controlled by the transmembrane voltage V and regulates the motion of particles overcoming potential barriers on going from the closed (open) to the open (closed) state of the channel); (iii) deterministic force responsible for interactions of water molecules with hydrophobic sites in the channel pore, and, finally, (iv) stochastic thermal fluctuation force. The model affords adequate approximation of experimental data.