A kinetic role for ionizable sites in membrane channel proteins

Electrically charged residues in a membrane channel protein will certainly have a direct effect upon its gating and selectivity if they are near the channel pore. It is customary to regard the charged state of such residues as a fixed feature of the channel. In this paper it is argued that far from being fixed, the charged state of ionizable residues near the pore will very probably change rapidly in response to the channel opening and to ions passing through it. Calculations are presented using simple models which demonstrate that changes in the dielectric environment and changes in the distances to other charged groups resulting from channel opening can shift the effective pK values of the sites by 3 or 4 units leading to switching of its charged state. Examples are given of how this time dependent charge state of ionizable residues may play an important role in the functioning of channels. Also, by considering the influence of the electric field due to the mobile ion upon the charge state of a residue in the channel wall, it is shown that a channel lined with acid residues may very effectively block the passage of cations while allowing the passage of anions.     

Effect of pore structure on energy barriers and applied voltage profiles. I. Symmetrical channels.

This paper presents calculations of the image potential for an ion in an aqueous pore spanning a lipid membrane and for the electric field produced in such a pore when a transmembrane potential is applied. The pore diameter may be variable. As long as the length-to-radius ratio in the narrow portion of a channel is large enough, the image potential for an ion in or near the mouth of a channel is determined by the geometry of the mouth. Within the constriction, the image potential of the ion-pore system may be reasonably approximated by constructing an "equivalent pore" of uniform diameter spanning a somewhat thinner membrane. When a transmembrane potential is applied the electric field within a constricted, constant radius, section of the model pore is constant. If the length-to-radius ratio of the narrow part of the channel is not too large or the channel ensemble has wide mouths, the field extends a significant distance into the aqueous region. The method is used to model features of the gramicidin A channel. The energy barrier for hydration (for exiting the channel) is identified with the activation energy for gramicidin conductance (Bamberg and Läuger, 1974, Biochim. Biophys. Acta. 367:127).     

Correlated ion flux through parallel pores: Application to channel subconductance states

Summary Many ion channels that normally gate fully open or shut have recently been observed occasionally to display well-defined subconductance states with conductances much less than those of the fully open channel. One model of this behavior is a channel consisting of several parallel pores with a strong correlation between the flux in each pore such that, normally, they all conduct together but, under special circumstances, the pores may transfer to a state in which only some of them conduct. This paper introduces a general technique for modeling correlated pores, and explores in detail by computer simulation a particular model based upon electric interaction between the pores. Correlation is obtained when the transient electric field of ions passing through the pores acts upon a common set of ionizable residues of the channel protein, causing transient changes in their effective pK and hence in their charged state. The computed properties of such a correlated parallel pore channel with single occupation of each pore are derived and compared to those predicted for a single pore that can contain more than one ion at a time and also to those predicted for a model pore with fluctuating barriers. Experiments that could distinguish between the present and previous models are listed.R.M.B. is grateful to the S.E.R.C. for the award of a graduate studentship.     

A simple method for the determination of the pore radius of ion channels in planar lipid bilayer membranes

Abstract A new method of pore size determination is presented. The results of applying this simple method to ion channels formed by staphylococcal α-toxin and its N-terminal fragment as well as to cholera toxin channels are shown. The advantages and the difficulties of this method are discussed. It was found that (i) the mobility of ions in solutions depends only on the percentage of concentration of added non-electrolytes and practically not on their chemical nature (sugars or polyglycols) and molecular size; (ii) the proportional change of both ion channel conductance and bulk solution conductivity by low M . non-electrolytes may be used as an indication of a diffusion mechanism of ion transport through channels; (iii) the slope of the dependence of the ion channel conductance on the bulk conductivity of solutions containing different concentrations of non-electrolyte is a good measure of channel permeability for non-electrolytes.     

Molecular dynamics simulations of the gramicidin A-dimyristoylphosphatidylcholine system with an ion in the channel pore region

To investigate the process of ion permeation in an ion channel systematically, we performed molecular dynamics (MD) simulations on a gramicidin A (GA)-phospholipid model system with an ion in the channel pore region. Each of the three types of ions (Ca2+, Na+ Cl-) was placed at five different positions along the channel axis by replacing a water molecule. MD simulations were performed on each system at constant pressure and constant temperature. The MD trajectories showed that the Ca2+ and Na+ ions could stably fluctuate in the pore region, but the Cl- ion was pushed out because of the unfavorable interaction with the channel. This result is consistent with experimental data. It was also found that the conformation of the GA channel underwent a significant change due to the presence of the ion, and the two ends of the GA monomer were more flexible than its middle region. In particular, the dramatic change of local pore radius near the ion indicated this kind of deformation. The strong interaction between the ion and carbonyl oxygen atoms of GA was the major contributor to this change. Furthermore, it was found that the ethanolamine group of the GA molecule was the most flexible group in the GA channel and often observed to block the entrance of GA. These results imply that the deformation of channel structure plays a very important factor in ion permeation, and the ethanolamine group may play a key role in regulating ion entry into the pore. In conclusion, our results indicate that the ion has a dominant influence on the structure of the GA channel and that the flexibility of the ion channel is a crucial factor in the ion permeation process.     

A microscopic electrostatic model for the amphotericin B channel.

A microscopic model of an amphotericin B channel is proposed. The structure of the pores is generated using the atomic coordinates of the molecule in the structure determined experimentally by X-ray diffraction. The net charges of the atoms are determined by Mulliken analysis. With these charges the electrostatic energy profiles are calculated for a monovalent ion passing through the channels formed by different number of antibiotic molecules having different radii. The water inside the channel was considered through a continuum medium using the dielectric constant of the bulk, and the membrane contribution was included using the virtual images of the pore in a dielectric slab of epsilon = 3. The model satisfactorily explains the permeability and selectivity characteristics as well as other observations yet unexplained. The electrostatic profiles obtained reinforce the hypothesis of the existence of channels formed by a variable number of units.     

The Behaviour of Ions in Narrow Water-Filled Pores

Today, the equilibrium behavior of ions in solution may be predicted with some confidence, essentially because rapid ionic diffusion over small distances ensures homogeneity throughout the solution. Equilibrium concepts such as ionic strength and pH apply. However, when attempting to understand the behavior of ions passing rapidly through narrow pores such as ion channels, no such equilibrium state may be assumed. The passing solution may have been in equilibrium with conditions at the mouth of the pore but will not be in equilibrium with charged molecules on the pore wall. In addition, the water in narrow pores will be partially ordered by contact with the pore walls and will not behave like bulk water.To illustrate this difference, a simple equilibrium calculation of the ion concentrations near a plastic sheet penetrated by narrow pores and containing in its surface partially ionized carboxyl groups is shown to be in good agreement with experiment. However, to predict the non-equilibrium behavior within the narrow pores is much more difficult. To illustrate the difficulty, a Monte Carlo computer model is described which attempts to predict the rapid switching of ion current observed experimentally with these narrow pores.     

Theoretical evaluation of cell membrane ion channel activation by applied magnetic fields

This letter re-examines a recently published calculation of the forces exerted on a membrane ion channel by a cation passing through in the presence of an externally applied magnetic field. We show here, in contradiction to the originally published calculation, that the forces generated due to the Lorentz force of the magnetic field on the cation are negligible compared with the forces required to activate an ion channel protein conformation change associated with the gating of the channel. Received: 11 August 1998 / Revised version: 25 October 1998 / Accepted: 11 November 1998     

Cryo‐electron microscopy structure of CLHM1 ion channel from Caenorhabditis elegans

Calcium homeostasis modulators (CALHMs/CLHMs) comprise a family of pore‐forming protein complexes assembling into voltage‐gated, Ca²⁺‐sensitive, nonselective channels. These complexes contain an ion‐conduction pore sufficiently wide to permit the passing of ATP molecules serving as neurotransmitters. While their function and structure information is accumulating, the precise mechanisms of these channel complexes remain to be full understood. Here, we present the structure of the Caenorhabditis elegans CLHM1 channel in its open state solved through single‐particle cryo‐electron microscopy at 3.7‐Å resolution. The transmembrane region of the channel structure of the dominant class shows an assembly of 10‐fold rotational symmetry in one layer, and its cytoplasmic region is involved in additional twofold symmetrical packing in a tail‐to‐tail manner. Furthermore, we identified a series of amino acid residues critical for the regulation of CeCLHM1 channel using functional assays, electrophysiological analyses as well as structural‐based analysis. Our structure and function analyses provide new insights into the mechanisms of CALHM channels.     

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.     

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.     

Theoretical studies of the M2 transmembrane segment of the glycine receptor: models of the open pore structure and current-voltage characteristics

The pentameric glycine receptor (GlyR), a member of the nicotinicoid superfamily of ligand-gated ion channels, is an inhibitory Cl(-) channel that is gated by glycine. Using recently published NMR data of the second transmembrane segment (M2) of the human alpha1 GlyR, structural models of pentameric assemblies embedded in a lipid bilayer were constructed using a combination of experimentally determined constraints coupled with all-atom energy minimization. Based on this structure of the pentameric M2 "pore", Brownian dynamics simulations of ion permeation through this putative conducting open state of the channel were carried out. Simulated I-V curves were in good agreement with published experimental current-voltage curves and the anion/cation permeability ratio, suggesting that our open-state model may be representative of the conducting channel of the full-length receptor. These studies also predicted regions of chloride occupancy and suggested residues critical to anion permeation. Calculations of the conductance of the cation-selective mutant A251E channel are also consistent with experimental data. In addition, both rotation and untilting of the pore helices of our model were found to be broadly consistent with closing of the channel, albeit at distinct regions that may reflect alternate gates of the receptor.     

The proton ladder,a static mechanism for ion/proton coports and counterports

An ion/proton counterport is formed simply by locating a chain of ionizable residues connected by a proton conducting path near a passive ion pore which spans the membrane. The electric coupling between the ion in transit through the pore and the residues can ensure that for each ion passing through the pore in one direction a proton is driven along the chain of ionizable residues (the proton ladder) in the same or in the opposite direction. The mechanism is symmetrical in that a trans-membrane ion gradient may drive protons against their electrochemical potential gradient or a proton gradient may drive ions against theirs. The mechanism is applicable to cation or anion channels and to coports or counterports. No mechanical motion is required other than the motion of the ions and the protons. Monte Carlo computer simulations are performed on the model and its predicted properties are listed. The new type of counterport model is compared with currently used models. Offprint requests to: D. T Edmonds     

A new double-chamber model of ion channels. Beyond the Hodgkin and Huxley model

This paper proposes a new double-chamber model (DCM) of ion channels. The model ion channel consists of a series of three pores alternating with two chambers. The chambers are net negatively charged. The chamber's electric charge originates from dissociated amino acid side chains and is pH dependent. The chamber's net negative charge is compensated by cations present inside the chamber and in a diffuse electric layer outside the chamber. The pore's permeability is constant independent of time. One pore of the sodium channel and one of the potassium channel is a voltage-sensing pore. Due to the channel's structure, ions flow through the pores and chambers in a time-dependent manner. The model reproduces experimental voltage clamp and action potential data. The current flowing through a single sodium channel is less then one femtoampere. The DCM is considerably simpler then the Hodgkin and Huxley model (HHM) used to describe the electrophysiological properties of an axon. Unlike the HHM, the DCM can explain refractoriness, anode break excitation, accommodation and the effect of pH and temperature on the channels without additional parameters. In the DCM, the axon membrane shows repetitive activity depending on the channel density, sodium to potassium channel ratio and external potassium concentration. In the DCM, the action potential starts from 'hot spot areas' of higher channel densities and a higher sodium to potassium channel ratio, and then propagates through the whole axon.     

短杆菌肽A-DMPC通道内离子输运的分子动力学模拟

用最近提出的构建膜体系初始构象的有效方法 ,构建了在DMPC脂膜环境下短杆菌肽A通道模型 (GA -DMPC)。通过对Na 、Ca2 、Cl-三种不同离子在GA -DMPC通道内不同位置的分子动力学模拟 ,研究离子在通道内输运过程中与通道及通道内水分子的相互作用 ,从分子动力学的角度阐明离子在通道内的输运机制。主要计算结果表明 :(1)离子在通道内的输运使GA的构象发生变化 ,GA的柔性是离子在通道内通透的重要因素 ;(2)Cl- 离子可扩大通道半径 ,Na 离子和Ca2 离子则减小通道半径。Cl-离子不能在GA通道内通透 ;(3)离子的出现使通道内水分子的偶极方向发生变化。上述结果均与实验相符。     

Relationship between pore occupancy and gating in BK potassium channels

Permeant ions can have significant effects on ion channel conformational changes. To further understand the relationship between ion occupancy and gating conformational changes, we have studied macroscopic and single-channel gating of BK potassium channels with different permeant monovalent cations. While the slopes of the conductance-voltage curve were reduced with respect to potassium for all permeant ions, BK channels required stronger depolarization to open only when thallium was the permeant ion. Thallium also slowed the activation and deactivation kinetics. Both the change in kinetics and the shift in the GV curve were dependent on the thallium passing through the permeation pathway, as well as on the concentration of thallium. There was a decrease in the mean open time and an increase in the number of short flicker closing events with thallium as the permeating ion. Mean closed durations were unaffected. Application of previously established allosteric gating models indicated that thallium specifically alters the opening and closing transition of the channel and does not alter the calcium activation or voltage activation pathways. Addition of a closed flicker state into the allosteric model can account for the effect of thallium on gating. Consideration of the thallium concentration dependence of the gating effects suggests that the flicker state may correspond to the collapsed selectivity filter seen in crystal structures of the KcsA potassium channel under the condition of low permeant ion concentration.     

Plants do it differently. A new basis for potassium/sodium selectivity in the pore of an ion channel

Understanding of the molecular architecture necessary for selective K(+) permeation through the pore of ion channels is based primarily on analysis of the crystal structure of the bacterial K(+) channel KcsA, and structure:function studies of cloned animal K(+) channels. Little is known about the conduction properties of a large family of plant proteins with structural similarities to cloned animal cyclic nucleotide-gated channels (CNGCs). Animal CNGCs are nonselective cation channels that do not discriminate between Na(+) and K(+) permeation. These channels all have the same triplet of amino acids in the channel pore ion selectivity filter, and this sequence is different from that of the selectivity filter found in K(+)-selective channels. Plant CNGCs have unique pore selectivity filters; unlike those found in any other family of channels. At present, the significance of the unique pore selectivity filters of plant CNGCs, with regard to discrimination between Na(+) and K(+) permeation is unresolved. Here, we present an electrophysiological analysis of several members of this protein family; identifying the first cloned plant channel (AtCNGC1) that conducts Na(+). Another member of this ion channel family (AtCNGC2) is shown to have a selectivity filter that provides a heretofore unknown molecular basis for discrimination between K(+) and Na(+) permeation. Specific amino acids within the AtCNGC2 pore selectivity filter (Asn-416, Asp-417) are demonstrated to facilitate K(+) over Na(+) conductance. The selectivity filter of AtCNGC2 represents an alternative mechanism to the well-known GYG amino acid triplet of K(+) channels that has been identified as the critical basis for K(+) over Na(+) permeation through the pore of ion channels.     

Brownian dynamics study of flux ratios in sodium channels

Measurements of unidirectional fluxes in ion channels provide one of the experimental methods for studying the steps involved in ion permeation in biological pores. Conventionally, the number of ions in the pore is inferred by fitting the ratio of inward and outward currents to an exponential function with an adjustable parameter known as the flux ratio exponent. Here we investigate the relationship between the number of ions in the pore and the flux ratio exponent in a model sodium channel under a range of conditions. Brownian dynamics simulations enable us to count the precise number of ions in the channel and at the same time measure the currents flowing across the pore in both directions. We show here that the values of the flux ratio exponent n′ ranges between 1 and 3 and is highly dependent on the ionic concentrations in which measurements are made. This is a consequence of the fact that both inward and outward currents are susceptible to saturation with increasing concentration. These results indicate that measurements of the flux ratio exponent cannot be directly related to the number of ions in the pore and that interpretation of such experimental measurements requires careful consideration of the conditions in which the study is made.     

Molecular dynamics investigation of the ω-current in the Kv1.2 voltage sensor domains

Voltage sensor domains (VSD) are transmembrane proteins that respond to changes in membrane voltage and modulate the activity of ion channels, enzymes, or in the case of proton channels allow permeation of protons across the cell membrane. VSDs consist of four transmembrane segments, S1-S4, forming an antiparallel helical bundle. The S4 segment contains several positively charged residues, mainly arginines, located at every third position along the helix. In the voltage-gated Shaker K(+) channel, the mutation of the first arginine of S4 to a smaller uncharged amino acid allows permeation of cations through the VSD. These currents, known as ω-currents, pass through the VSD and are distinct from K(+) currents passing through the main ion conduction pore. Here we report molecular dynamics simulations of the ω-current in the resting-state conformation for Kv1.2 and for four of its mutants. The four tested mutants exhibit various degrees of conductivity for K(+) and Cl(-) ions, with a slight selectivity for K(+) over Cl(-). Analysis of the ion permeation pathway, in the case of a highly conductive mutant, reveals a negatively charged constriction region near the center of the membrane that might act as a selectivity filter to prevent permeation of anions through the pore. The residues R1 in S4 and E1 in S2 are located at the narrowest region of the ω-pore for the resting state conformation of the VSD, in agreement with experiments showing that the largest increase in current is produced by the double mutation E1D and R1S.