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.     

Modeling noncompetitive antagonism of a nicotinic acetylcholine receptor

Models of closed and open channel pores of a muscle-type nicotinic acetylcholine receptor (nAChR) channel comprising M1 and M2 segments are presented. A model of the closed channel is proposed in which hydrophobic residues of the Equatorial Leucine ring screen the oxygen domain formed by the Serine ring, thereby preventing ion flux without completely occluding the pore. This model demonstrates a high similarity with the structure derived from a recent electron microscopy study. We propose that hydrophobic residues of the Equatorial Leucine ring are retracted when the pore is open. Our models provide a possible resolution of the nAChR gate controversy. We have also obtained explanations for the complex mechanisms underlying inhibition of nAChR by philanthotoxins (PhTXs). PhTX-343, containing a spermine moiety with a charge of +3, binds deep in the pore near the Serine ring where classical open channel blockers of nAChR bind. In contrast, PhTX-(12), which has a single charged amino group is unable to reach deeply located rings because of steric restrictions. Both philanthotoxins may bind to a hydrophobic site located close to the external entrance of the pore in a region that includes residues associated with the regulation of desensitization.     

Ion channel gates: comparative analysis of energy barriers

The energetic profile of an ion translated along the axis of an ion channel should reveal whether the structure corresponds to a functionally open or closed state of the channel. In this study, we explore the combined use of Poisson–Boltzmann electrostatic calculations and evaluation of van der Waals interactions between ion and pore to provide an initial appraisal of the gating state of a channel. This approach is exemplified by its application to the bacterial inward rectifier potassium channel KirBac3.1, where it reveals the closed gate to be formed by a ring of leucine (L124) side chains. We have extended this analysis to a comparative survey of gating profiles, including model hydrophobic nanopores, the nicotinic acetylcholine receptor, and a number of potassium channel structures and models. This enables us to identify three gating regimes, and to show the limitation of this computationally inexpensive method. For a (closed) gate radius of 0.4 nm < R < 0.8 nm, a hydrophobic gate may be present. For a gate radius of 0.2 nm < R < 0.4 nm, both electrostatic and van der Waals interactions will contribute to the barrier height. Below R = 0.2 nm, repulsive van der Waals interactions are likely to dominate, resulting in a sterically occluded gate. In general, the method is more useful when the channel is wider; for narrower channels, the flexibility of the protein may allow otherwise-unsurmountable energetic barriers to be overcome.     

Computational study of the transmembrane domain of the acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein whose transmembrane domain (TM-domain) is believed to be responsible for channel gating via a hydrophobic effect. In this work, we perform molecular dynamics and Brownian dynamics simulations to investigate the effect of transmembrane potential on the conformation and water occupancy of TM-domain, and the resulting ion permeation events. The results show that the behavior of the hydrophobic gate is voltage-dependent. Large hyperpolarized membrane potential can change the conformation of TM-domain and water occupancy in this region, which may enable ion conduction. An electrostatic gating mechanism is also proposed from our simulations, which seems to play a role in addition to the well-known hydrophobic effect.     

Molecular dynamics study of water and Na+ ions in models of the pore region of the nicotinic acetylcholine receptor.

The nicotinic acetylcholine receptor (nAChR) is an integral membrane protein that forms ligand-gated and cation-selective channels. The central pore is lined by a bundle of five approximately parallel M2 helices, one from each subunit. Candidate model structures of the solvated pore region of a homopentameric (alpha7)5 nAChR channel in the open state, and in two possible forms of the closed state, have been studied using molecular dynamics simulations with restraining potentials. It is found that the mobility of the water is substantially lower within the pore than in bulk, and the water molecules become aligned with the M2 helix dipoles. Hydrogen-bonding patterns in the pore, especially around pore-lining charged and hydrophilic residues, and around exposed regions of the helix backbone, have been determined. Initial studies of systems containing both water and sodium ions together within the pore region have also been conducted. A sodium ion has been introduced into the solvated models at various points along the pore axis and its energy profile evaluated. It is found that the ion causes only a local perturbation of the water structure. The results of these calculations have been used to examine the effectiveness of the central ring of leucines as a component of a gate in the closed-channel model.     

An energy-efficient gating mechanism in the acetylcholine receptor channel suggested by molecular and Brownian dynamics

Acetylcholine receptors mediate electrical signaling between nerve and muscle by opening and closing a transmembrane ion conductive pore. Molecular and Brownian dynamics simulations are used to shed light on the location and mechanism of the channel gate. Four separate 5 ns molecular dynamics simulations are carried out on the imaged structure of the channel, a hypothetical open structure with a slightly wider pore and a mutant structure in which a central ring of hydrophobic residues is replaced by polar groups. Water is found to partially evacuate the pore during molecular simulations of the imaged structure, whereas ions face a large energy barrier and do not conduct through the channel in Brownian dynamics simulations. The pore appears to be in a closed configuration despite containing an unobstructed pathway across the membrane as a series of hydrophobic residues in the center of the channel provide an unfavorable home to water and ions. When the channel is widened slightly, water floods into the channel and ions conduct at a rate comparable to the currents measured experimentally in open channels. The pore remains permeable to ions provided the extracellular end of the pore-lining helix is restrained near the putative open configuration to mimic the presence of the ligand binding domain. Replacing some of the hydrophobic residues with polar ones decreases the barrier for ion permeation but does not result in significant currents. The channel is posited to utilize an energy efficient gating mechanism in which only minor conformational changes of the hydrophobic region of the pore are required to create macroscopic changes in conductance.     

Water dynamics and dewetting transitions in the small mechanosensitive channel MscS

The dynamics of confined water in capillaries and nanotubes suggests that gating of ion channels may involve not only changes of the pore geometry, but also transitions between water-filled and empty states in certain locations. The recently solved heptameric structure of the small mechanosensitive channel of Escherichia coli, MscS, has revealed a relatively wide (7-15 A) yet highly hydrophobic transmembrane pore. Continuum estimations based on the properties of pore surface suggest low conductance and a thermodynamic possibility of dewetting. To test the predictions we performed molecular dynamics simulations of MscS filled with flexible TIP3P water. Irrespective to the initial conditions, several independent 6-ns simulations converged to the same stable state with the pore water-filled in the wider part, but predominantly empty in the narrow hydrophobic part, displaying intermittent vapor-liquid transitions. The polar gain-of-function substitution L109S in the constriction resulted in a stable hydration of the entire pore. Steered passages of Cl(-) ions through the narrow part of the pore consistently produced partial ion dehydration and required a force of 200-400 pN to overcome an estimated barrier of 10-20 kcal/mole, implying negligibly low conductance. We conclude that the crystal structure of MscS does not represent an open state. We infer that MscS gate, which is similar to that of the nicotinic ACh receptor, involves a vapor-lock mechanism where limited changes of geometry or surface polarity can locally switch the regime between water-filled (conducting) and empty (nonconducting) states.     

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.     

Asp433 in the closing gate of ASIC1 determines stability of the open state without changing properties of the selectivity filter or Ca2+ block

A constriction formed by the crossing of the second transmembrane domains of ASIC1, residues G432 to G436, forms the narrowest segment of the pore in the crystal structure of chicken ASIC1, presumably in the desensitized state, suggesting that it constitutes the "desensitization gate" and the "selectivity filter." Residues Gly-432 and Asp-433 occlude the pore, preventing the passage of ions from the extracellular side. Here, we examined the role of Asp-433 and Gly-432 in channel kinetics, ion selectivity, conductance, and Ca(2+) block in lamprey ASIC1 that is a channel with little intrinsic desensitization in the pH range of maximal activity, pH 7.0. The results show that the duration of open times depends on residue 433, with Asp supporting the longest openings followed by Glu, Gln, or Asn, whereas other residues keep the channel closed. This is consistent with residue Asp-433 forming the pore's closing gate and the properties of the side chain either stabilizing (hydrophobic amino acids) or destabilizing (Asp) the gate. The data also show residue 432 influencing the duration of openings, but here only Gly and Ala support long openings, whereas all other residues keep channels closed. The negative charge of Asp-433 was not required for block of the open pore by Ca(2+) or for determining ion selectivity and unitary conductance. We conclude that the conserved residue Asp-433 forms the closing gate of the pore and thereby determines the duration of individual openings while desensitization, defined as the permanent closure of all or a fraction of channels by the continual presence of H(+), modulates the on or off position of the closing gate. The latter effect depends on less conserved regions of the channel, such as TM1 and the extracellular domain. The constriction made by Asp-433 and Gly-432 does not select for ions in the open conformation, implying that the closing gate and selectivity filter are separate structural elements in the ion pathway of ASIC1. The results also predict a significantly different conformation of TM2 in the open state that relieves the constriction made by TM2, allowing the passage of ions unimpeded by the side chain of Asp-433.     

Homology modeling and molecular dynamics simulations of transmembrane domain structure of human neuronal nicotinic acetylcholine receptor

A three-dimensional model of the transmembrane domain of a neuronal-type nicotinic acetylcholine receptor (nAChR), (alpha4)2(beta2)3, was constructed from a homology structure of the muscle-type nAChR recently determined by cryo-electron microscopy. The neuronal channel model was embedded in a fully hydrated DMPC lipid bilayer, and molecular-dynamics simulations were performed for 5 ns. A comparative analysis of the neuronal- versus muscle-type nAChR models revealed many conserved pore-lining residues, but an important difference was found near the periplasmic mouth of the pore. A flickering salt-bridge of alpha4-E266 with its adjacent beta2-K260 was observed in the neuronal-type channel during the course of the molecular-dynamics simulations. The narrowest region, with a pore radius of approximately 2 A formed by the salt-bridges, does not seem to be the restriction site for a continuous water passage. Instead, two hydrophobic rings, formed by alpha4-V259, alpha4-L263, and the homologous residues in the beta2-subunits, act as the gates for water flow, even though the region has a slightly larger pore radius. The model offers new insight into the water transport across the (alpha4)2(beta2)3 nAChR channel, and may lead to a better understanding of the structures, dynamics, and functions of this family of ion channels.     

Molecular restraints in the permeation pathway of ion channels

Ion channels assist and control the diffusion of ions through biological membranes. The conduction process depends on the structural characteristics of these nanopores, among which are the hydrophobicity and the afforded diameter of the conduction pathway. In this contribution, we use full atomistic free-energy molecular dynamics simulations to estimate the effect of such characteristics on the energetics of ion conduction through the activation gate of voltage-gated potassium (Kv) channels. We consider specifically the ionic translocation through three different permeation pathways, corresponding to the activation gate of an atomistic model of Shaker channels in closed and partially opened conformations, and that of the open conformation of the Kv1.2 channel. In agreement with experiments, we find that the region of Val(478) constitutes the main gate. The conduction is unfavorable through this gate when the constriction is smaller than an estimated threshold of 4.5-5.0 A, mainly due to incomplete coordination-hydration of the ion. Above this critical size, e.g., for the Kv1.2, the valine gate is wide enough to allow fully coordination of the ion and therefore its diffusion on a flat energy surface. Similar to other ion channels, Kv channels appear therefore to regulate diffusion by constricting hydrophobic regions of the permeation pathway.     

Probing structure-function relationships and gating mechanisms in the CorA Mg2+ transport system

Recent crystal structures of the CorA Mg(2+) transport protein from Thermotoga maritima (TmCorA) revealed an unusually long ion pore putatively gated by hydrophobic residues near the intracellular end and by universally conserved asparagine residues at the periplasmic entrance. A conformational change observed in an isolated funnel domain structure also led to a proposal for the structural basis of gating. Because understanding the molecular mechanisms underlying ion channel and transporter gating remains an important challenge, we have undertaken a structure-guided engineering approach to probe structure-function relationships in TmCorA. The intracellular funnel domain is shown to constitute an allosteric regulatory module that can be engineered to promote an activated or closed state. A periplasmic gate centered about a proline-induced kink of the pore-lining helix is described where "helix-straightening" mutations produce a dramatic gain-of-function. Mutation to the narrowest constriction along the pore demonstrates that a hydrophobic gate is operational within this Mg(2+)-selective transport protein and likely forms an energetic barrier to ion flux. We also provide evidence that highly conserved acidic residues found in the short periplasmic loop are not essential for TmCorA function or Mg(2+) selectivity but may be required for proper protein folding and stability. This work extends our gating model for the CorA-Alr1-Mrs2 superfamily and reveals features that are characteristic of an ion channel. Aspects of these results that have broader implications for a range of channel and transporter families are highlighted.     

Structural effects of divalent calcium cations on the α7 nicotinic acetylcholine receptor: A molecular dynamics simulation study

The α7 subtype of neuronal nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein that is vital to various neurological functions, including modulation of neurotransmitter release. A relatively high concentration of extracellular Ca²⁺ in the neuronal environment is likely to exert substantial structural and functional influence on nAChRs, which may affect their interactions with agonists and antagonists. In this work, we employed atomistic molecular dynamics (MD) simulations to examine the effects of elevated Ca²⁺ on the structure and dynamics of α7 nAChR embedded in a model phospholipid bilayer. Our results suggest that the presence of Ca²⁺ in the α7 nAChR environment results in closure of loop C-in the extracellular ligand-binding domain, a motion normally associated with agonist binding and receptor activation. Elevated Ca²⁺ also alters the conformation of key regions of the receptor, including the inter-helical loops, pore-lining helices and the “gate” residues, and causes partial channel opening in the absence of an agonist, leading to an attendant reduction in the free energy of Ca²⁺ permeation through the pore as elucidated by umbrella sampling simulations. Overall, the structural and permeability changes in α7 nAChR suggest that elevated Ca²⁺ induces a partially activated receptor state that is distinct from both the resting and the agonist-activated states. These results are consistent with the notion that divalent ions can serve as a potentiator of nAChRs, resulting in a higher rate of receptor activation (and subsequent desensitization) in the presence of agonists, with possible implications for diseases involving calcium dysregulation.     

Stabilizing the closed S6 gate in the Shaker Kv channel through modification of a hydrophobic seal

The primary activation gate in K+ channels is thought to reside near the intracellular entrance to the ion conduction pore. In a previous study of the S6 activation gate in Shaker (Hackos et al., 2002), we found that mutation of V478 to W results in a channel that cannot conduct ions even though the voltage sensors are competent to translocate gating charge in response to membrane depolarization. In the present study we explore the mechanism underlying the nonconducting phenotype in V478W and compare it to that of W434F, a mutation located in an extracellular region of the pore that is nonconducting because the channel is predominantly found in an inactivated state. We began by examining whether the intracellular gate moves using probes that interact with the intracellular pore and by studying the inactivation properties of heterodimeric channels that are competent to conduct ions. The results of these experiments support distinct mechanisms underlying nonconduction in W434F and V478W, suggesting that the gate in V478W either remains closed, or that the mutation has created a large barrier to ion permeation in the open state. Single channel recordings for heterodimeric and double mutant constructs in which ion conduction is rescued suggest that the V478W mutation does not dramatically alter unitary conductance. Taken together, our results suggest that the V478W mutation causes a profound shift of the closed to open equilibrium toward the closed state. This mechanism is discussed in the context of the structure of this critical region in K+ channels.     

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.     

Role of acetylcholine receptor domains in ion selectivity

The nicotinic acetylcholine receptor (nAChR) is a ligand gated ion channel protein, composed of three domains: a transmembrane domain (TM-domain), extracellular domain (EC-domain), and intracellular domain (IC-domain). Due to its biological importance, much experimental and theoretical research has been carried out to explore its mechanisms of gating and selectivity, but there are still many unresolved issues, especially on the ion selectivity. Moreover, most of the previous theoretical work has concentrated on the TM-domain or EC-domain of nAChR, which may be insufficient to understand the entire structure–function relation. In this work, we perform molecular dynamics, Brownian dynamics simulations and continuum electrostatic calculations to investigate the role of different nAChR domains in ion conduction and selectivity. The results show that although both the EC and IC domains contain strong negative charges that create large cation concentrations at either end of the pore, this alone is not sufficient to create the observed cation selectivity and may play a greater role in determining the channel conductance. The presence of cations in the wide regions of the pore can screen out the protein charge allowing anions to enter, meaning that local regions of the TM-domain are most likely responsible for discriminating between ions. These new results complement our understanding about the ion conduction and selectivity mechanism of nAChR.     

Electrostatics and the ion selectivity of ligand-gated channels.

The nicotinic acetylcholine receptor (nAChR) is a cation-selective ion channel that opens in response to acetylcholine binding. The related glycine receptor (GlyR) is anion selective. The pore-lining domain of each protein may be modeled as a bundle of five parallel M2 helices. Models of the pore-lining domains of homopentameric nAChR and GlyR have been used in continuum electrostatics calculations to probe the origins of ion selectivity. Calculated pKA values suggest that "rings" of acidic or basic side chains at the mouths of the nAChR or GlyR M2 helix bundles, respectively, may not be fully ionized. In particular, for the nAChR the ring of glutamate side chains at the extracellular mouth of the pore is predicted to be largely protonated at neutral pH, whereas those glutamate side chains in the intracellular and intermediate rings (at the opposite mouth of the pore) are predicted to be fully ionized. Inclusion of the other domains of each protein represented as an irregular cylindrical tube in which the M2 bundles are embedded suggests that both the M2 helices and the extramembrane domains play significant roles in determining ion selectivity.     

Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics.

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.     

Electrostatics of the intracellular vestibule of K+ channels

Previous calculations using continuum electrostatic calculations showed that a fully hydrated monovalent cation is electrostatically stabilized at the center of the cavity of the KcsA potassium channel. Further analysis demonstrated that this cavity stabilization was controlled by a balance between the unfavorable reaction field due to the finite size of the cavity and the favorable electrostatic field arising from the pore helices. In the present study, continuum electrostatic calculations are used to investigate how the stability of an ion in the intracellular vestibular cavity common to known potassium channels is affected as the inner channel gate opens and the cavity becomes larger and contiguous with the intracellular solution. The X-ray structure of the calcium-activated potassium channel MthK, which was crystallized in the open state, is used to construct models of the KcsA channel in the open state. It is found that, as the channel opens, the barrier at the helix bundle crossing decreases to approximately 0 kcal/mol, but that the ion in the cavity is also significantly destabilized. The results are compared and contrasted with additional calculations performed on the KvAP (voltage-activated) and KirBac1.1 (inward rectifier) channels, as well as models of the pore domain of Shaker in the open and closed state. In conclusion, electrostatic factors give rise to energetic constraints on ion permeation that have important functional consequences on the various K+ channels, and partly explain the presence or absence of charged residues near the inner vestibular entry.