A hydrophobic gate in an ion channel: the closed state of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the 'Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid and serotonin. Cryo-electron microscopy has yielded a three-dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height about 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united-atom force field for the protein instead of an all-atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.     

A data-driven model of a modal gated ion channel: The inositol 1,4,5-trisphosphate receptor in insect Sf9 cells

The inositol 1,4,5-trisphosphate (IP(3)) receptor (IP(3)R) channel is crucial for the generation and modulation of intracellular Ca(2+) signals in animal cells. To gain insight into the complicated ligand regulation of this ubiquitous channel, we constructed a simple quantitative continuous-time Markov-chain model from the data. Our model accounts for most experimentally observed gating behaviors of single native IP(3)R channels from insect Sf9 cells. Ligand (Ca(2+) and IP(3)) dependencies of channel activity established six main ligand-bound channel complexes, where a complex consists of one or more states with the same ligand stoichiometry and open or closed conformation. Channel gating in three distinct modes added one complex and indicated that three complexes gate in multiple modes. This also restricted the connectivity between channel complexes. Finally, latencies of channel responses to abrupt ligand concentration changes defined a model with specific network topology between 9 closed and 3 open states. The model with 28 parameters can closely reproduce the equilibrium gating statistics for all three gating modes over a broad range of ligand concentrations. It also captures the major features of channel response latency distributions. The model can generate falsifiable predictions of IP(3)R channel gating behaviors and provide insights to both guide future experiment development and improve IP(3)R channel gating analysis. Maximum likelihood estimates of the model parameters and of the parameters in the De Young-Keizer model yield strong statistical evidence in favor of our model. Our method is simple and easily applicable to the dynamics of other ion channels and molecules.     

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.     

The molecular mechanism by which PIP(2) opens the intracellular G-loop gate of a Kir3.1 channel

Inwardly rectifying potassium (Kir) channels are characterized by a long pore comprised of continuous transmembrane and cytosolic portions. A high-resolution structure of a Kir3.1 chimera revealed the presence of the cytosolic (G-loop) gate captured in the closed or open conformations. Here, we conducted molecular-dynamics simulations of these two channel states in the presence and absence of phosphatidylinositol bisphosphate (PIP(2)), a phospholipid that is known to gate Kir channels. Simulations of the closed state with PIP(2) revealed an intermediate state between the closed and open conformations involving direct transient interactions with PIP(2), as well as a network of transitional inter- and intrasubunit interactions. Key elements in the G-loop gating transition involved a PIP(2)-driven movement of the N-terminus and C-linker that removed constraining intermolecular interactions and led to CD-loop stabilization of the G-loop gate in the open state. To our knowledge, this is the first dynamic molecular view of PIP(2)-induced channel gating that is consistent with existing experimental data.     

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.     

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.     

Voltage-dependent hydration and conduction properties of the hydrophobic pore of the mechanosensitive channel of small conductance

A detailed picture of water and ion properties in small pores is important for understanding the behavior of biological ion channels. Several recent modeling studies have shown that small, hydrophobic pores exclude water and ions even if they are physically large enough to accommodate them, a mechanism called hydrophobic gating. This mechanism has been implicated in the gating of several channels, including the mechanosensitive channel of small conductance (MscS). Although the pore in the crystal structure of MscS is wide and was initially hypothesized to be open, it is lined by hydrophobic residues and may represent a nonconducting state. Molecular dynamics simulations were performed on MscS to determine whether or not the structure can conduct ions. Unlike previous simulations of hydrophobic nanopores, electric fields were applied to this system to model the transmembrane potential, which proved to be important. Although simulations without a potential resulted in a dehydrated, occluded pore, the application of a potential increased the hydration of the pore and resulted in current flow through the channel. The calculated channel conductance was in good agreement with experiment. Therefore, it is likely that the MscS crystal structure is closer to a conducting than a nonconducting state.     

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.     

Cross talk between activation and slow inactivation gates of Shaker potassium channels

This study addresses the energetic coupling between the activation and slow inactivation gates of Shaker potassium channels. To track the status of the activation gate in inactivated channels that are nonconducting, we used two functional assays: the accessibility of a cysteine residue engineered into the protein lining the pore cavity (V474C) and the liberation by depolarization of a Cs(+) ion trapped behind the closed activation gate. We determined that the rate of activation gate movement depends on the state of the inactivation gate. A closed inactivation gate favors faster opening and slower closing of the activation gate. We also show that hyperpolarization closes the activation gate long before a channel recovers from inactivation. Because activation and slow inactivation are ubiquitous gating processes in potassium channels, the cross talk between them is likely to be a fundamental factor in controlling ion flux across membranes.     

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.     

Some unexpected consequences of a simple physical mechanism for voltage-dependent gating in biological membranes

We consider a model for voltage-dependent gating of channels in which the gating charges are on the channel wall and move only a small distance. When this movement occurs across the closed gate, the charges move through the entire transmembrane potential, which is energetically equivalent to their moving across the entire membrane. The channel exists in two open states, O1 and O2, and two closed states, C1 and C2; each open and closed configuration is divided into two states because of the two possible positions of the gating charges. An unusual property of this model is that the electrical work in going from an open to a closed configuration (for example, in going from O1 to C2) is path dependent, and net work can result from going reversibly around a complete cycle. The model channel, like many biological channels, shows bursting activity. This flickering on and off of the channel enables the gate to sense the electric field and decide if it should be in the open or closed configuration. We prove here some general theorms concerning the electrical work associated with the movements of the walls of channels and the movements of charges on these walls.     

Voltage-dependent conformational changes in connexin channels

Channels formed by connexins display two distinct types of voltage-dependent gating, termed V(j)- or fast-gating and loop- or slow-gating. Recent studies, using metal bridge formation and chemical cross-linking have identified a region within the channel pore that contributes to the formation of the loop-gate permeability barrier. The conformational changes are remarkably large, reducing the channel pore diameter from 15 to 20? to less than 4?. Surprisingly, the largest conformational change occurs in the most stable region of the channel pore, the 3(10) or parahelix formed by amino acids in the 42-51 segment. The data provide a set of positional constraints that can be used to model the structure of the loop-gate closed state. Less is known about the conformation of the V(j)-gate closed state. There appear to be two different mechanisms; one in which conformational changes in channel structure are linked to a voltage sensor contained in the N-terminus of Cx26 and Cx32 and a second in which the C-terminus of Cx43 and Cx40 may act either as a gating particle to block the channel pore or alternatively to stabilize the closed state. The later mechanism utilizes the same domains as implicated in effecting pH gating of Cx43 channels. It is unclear if the two V(j)-gating mechanisms are related or if they represent different gating mechanisms that operate separately in different subsets of connexin channels. A model of the V(j)-closed state of Cx26 hemichannel that is based on the X-ray structure of Cx26 and electron crystallographic structures of a Cx26 mutation suggests that the permeability barrier for V(j)-gating is formed exclusively by the N-terminus, but recent information suggests that this conformation may not represent a voltage-closed state. Closed state models are considered from a thermodynamic perspective based on information from the 3.5? Cx26 crystal structure and molecular dynamics (MD) simulations. The applications of computational and experimental methods to define the path of allosteric molecular transitions that link the open and closed states are discussed. This article is part of a Special Issue entitled: The Communicating junctions, composition, structure and characteristics.     

Conformational Transitions Underlying Pore Opening and Desensitization in Membrane-embedded Gloeobacter violaceus Ligand-gated Ion Channel (GLIC)

Direct structural insight into the mechanisms underlying activation and desensitization remain unavailable for the pentameric ligand-gated channel family. Here, we report the structural rearrangements underlying gating transitions in membrane-embedded GLIC, a prokaryotic homologue, using site-directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy. We particularly probed the conformation of pore-lining second transmembrane segment (M2) under conditions that favor the closed and the ligand-bound desensitized states. The spin label mobility, intersubunit spin-spin proximity, and the solvent-accessibility parameters in the two states clearly delineate the underlying protein motions within M2. Our results show that during activation the extracellular hydrophobic region undergoes major changes involving an outward translational movement, away from the pore axis, leading to an increase in the pore diameter, whereas the lower end of M2 remains relatively immobile. Most notably, during desensitization, the intervening polar residues in the middle of M2 move closer to form a solvent-occluded barrier and thereby reveal the location of a distinct desensitization gate. In comparison with the crystal structure of GLIC, the structural dynamics of the channel in a membrane environment suggest a more loosely packed conformation with water-accessible intrasubunit vestibules penetrating from the extracellular end all the way to the middle of M2 in the closed state. These regions have been implicated to play a major role in alcohol and drug modulation. Overall, these findings represent a key step toward understanding the fundamentals of gating mechanisms in this class of channels.     

Forced gating motions by a substituted titratable side chain at the bundle crossing of a potassium channel

Numerous inwardly rectifying potassium (Kir) channels possess an aromatic residue in the helix bundle crossing region, forming the narrowest pore constriction in crystal structures. However, the role of the Kir channel bundle crossing as a functional gate remains uncertain. We report a unique phenotype of Kir6.2 channels mutated to encode glutamate at this position (F168E). Despite a prediction of four glutamates in close proximity, Kir6.2(F168E) channels are predominantly closed at physiological pH, whereas alkalization causes rapid and reversible channel activation. These findings suggest that F168E glutamates are uncharged at physiological pH but become deprotonated at alkaline pH, forcing channel opening due to mutual repulsion of nearby negatively charged side chains. The potassium channel pore scaffold likely brings these glutamates close together, causing a significant pK(a) shift relative to the free side chain (as seen in the KcsA selectivity filter). Alkalization also shifts the apparent ATP sensitivity of the channel, indicating that forced motion of the bundle crossing is coupled to the ATP-binding site and may resemble conformational changes involved in wild-type Kir6.2 gating. The study demonstrates a novel mechanism for engineering extrinsic control of channel gating by pH and shows that conformational changes in the bundle crossing region are involved in ligand-dependent gating of Kir channels.     

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.     

Evidence for a deep pore activation gate in small conductance Ca2+-activated K+ channels

Small conductance calcium-gated potassium (SK) channels share an overall topology with voltage-gated potassium (K(v)) channels, but are distinct in that they are gated solely by calcium (Ca(2+)), not voltage. For K(v) channels there is strong evidence for an activation gate at the intracellular end of the pore, which was not revealed by substituted cysteine accessibility of the homologous region in SK2 channels. In this study, the divalent ions cadmium (Cd(2+)) and barium (Ba(2+)), and 2-aminoethyl methanethiosulfonate (MTSEA) were used to probe three sites in the SK2 channel pore, each intracellular to (on the selectivity filter side of) the region that forms the intracellular activation gate of voltage-gated ion channels. We report that Cd(2+) applied to the intracellular side of the membrane can modify a cysteine introduced to a site (V391C) just intracellular to the putative activation gate whether channels are open or closed. Similarly, MTSEA applied to the intracellular side of the membrane can access a cysteine residue (A384C) that, based on homology to potassium (K) channel crystal structures (i.e., the KcsA/MthK model), resides one amino acid intracellular to the glycine gating hinge. Cd(2+) and MTSEA modify with similar rates whether the channels are open or closed. In contrast, Ba(2+) applied to the intracellular side of the membrane, which is believed to block at the intracellular end of the selectivity filter, blocks open but not closed channels when applied to the cytoplasmic face of rSK2 channels. Moreover, Ba(2+) is trapped in SK2 channels when applied to open channels that are subsequently closed. Ba(2+) pre-block slows MTSEA modification of A384C in open but not in closed (Ba(2+)-trapped) channels. The findings suggest that the SK channel activation gate resides deep in the vestibule of the channel, perhaps in the selectivity filter itself.     

Solution structure of Phrixotoxin 1, a specific peptide inhibitor of Kv4 potassium channels from the venom of the theraphosid spider Phrixotrichus auratus

Animal toxins block voltage-dependent potassium channels (Kv) either by occluding the conduction pore (pore blockers) or by modifying the channel gating properties (gating modifiers). Gating modifiers of Kv channels bind to four equivalent extracellular sites near the S3 and S4 segments, close to the voltage sensor. Phrixotoxins are gating modifiers that bind preferentially to the closed state of the channel and fold into the Inhibitory Cystine Knot structural motif. We have solved the solution structure of Phrixotoxin 1, a gating modifier of Kv4 potassium channels. Analysis of the molecular surface and the electrostatic anisotropy of Phrixotoxin 1 and of other toxins acting on voltage-dependent potassium channels allowed us to propose a toxin interacting surface that encompasses both the surface from which the dipole moment emerges and a neighboring hydrophobic surface rich in aromatic residues.     

Molecular dynamics simulation of TMEM16A channel: Linking structure with gating

TMEM16A, the calcium-activated chloride channel, is broadly expressed and plays pivotal roles in diverse physiological processes. To understand the structural and functional relationships of TMEM16A, it is necessary to fully clarify the structural basis of the gating of the TMEM16A channel. Herein, we performed the protein electrostatic analysis and molecular dynamics simulation on the TMEM16A in the presence and absence of Ca²⁺. Data showed that the separation of TM4 and TM6 causes pore expansion, and Q646 may be a key residue for the formation of π-helix in the middle segment of TM6. Moreover, E705 was found to form a group of H-bond interactions with D554/K588/K645 below the hydrophobic gate to stabilize the closed conformation of the pore in the Ca²⁺-free state. Interestingly, in the Ca²⁺ bound state, the E705 side chain swings 100^o to serve as Ca²⁺-binding coordination and released K645. K645 is closer to the hydrophobic gate in the calcium-bound state, which facilitates the provision of electrostatic forces for chloride ions as the ions pass through the hydrophobic gate. Our findings provide the structural-based insights to understanding the mechanisms of gating of TMEM16A.     

Gating of the ATP-sensitive K+ channel by a pore-lining phenylalanine residue

ATP-sensitive K(+) (K(ATP)) channels are gated by intracellular ATP, proton and phospholipids. The pore-forming Kir6.2 subunit has all essential machineries for channel gating by these ligands. It is known that channel gating involves the inner helix bundle of crossing in which a phenylalanine residue (Phe168) is found in the TM2 at the narrowest region of the ion-conduction pathway in the Kir6.2. Here we present evidence that Phe168-Kir6.2 functions as an ATP- and proton-activated gate via steric hindrance and hydrophobic interactions. Site-specific mutations of Phe168 to a small amino acid resulted in losses of the ATP- and proton-dependent gating, whereas the channel gating was well maintained after mutation to a bulky tryptophan, supporting the steric hindrance effect. The steric hindrance effect, though necessary, was insufficient for the gating, as mutating Phe168 to a bulky hydrophilic residue severely compromised the channel gating. Single-channel kinetics of the F168W mutant resembled the wild-type channel. Small residues increased P(open), and displayed long-lasting closures and long-lasting openings. Kinetic modeling showed that these resulted from stabilization of the channel to open and long-lived closed states, suggesting that a bulky and hydrophobic residue may lower the energy barrier for the switch between channel openings and closures. Thus, it is likely that the Phe168 acts as not only a steric hindrance gate but also potentially a facilitator of gating transitions in the Kir6.2 channel.