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.     

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.     

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.     

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 Conduction in Ligand-Gated Ion Channels: Brownian Dynamics Studies of Four Recent Crystal Structures

Four x-ray crystal structures of prokaryotic homologs of ligand-gated ion channels have recently been determined: ELIC from Erwinia chrysanthemi, two structures of a proton-activated channel from Gloebacter violaceus (GLIC1 and GLIC2) and that of the E221A mutant (GLIC1M). The availability of numerous structures of channels in this family allows for aspects of channel gating and ion conduction to be examined. Here, we determine the likely conduction states of the four structures as well as IV curves, ion selectivity, and steps involved in ion permeation by performing extensive Brownian dynamics simulations. Our results show that the ELIC structure is indeed nonconductive, but that GLIC1 and GLIC1M are both conductive of ions with properties different from those seen in experimental studies of the channel. GLIC2 appears to reflect an open state of the channel with a predicted conductance of 10.8-12.4 pS in 140 mM NaCl solution, which is comparable to the experimental value 8 ± 2 pS. The extracellular domain of the channel is shown to have an important influence on the channel current, but a less significant role in ion selectivity.     

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.     

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.     

Electrostatic couplings in OmpA ion-channel gating suggest a mechanism for pore opening

The molecular forces that drive structural transitions between the open and closed states of channels and transporters are not well understood. The gate of the OmpA channel is formed by the central Glu52-Arg138 salt bridge, which can open to form alternate ion pairs with Lys82 and Glu128. To gain deeper insight into the channel-opening mechanism, we measured interaction energies between the relevant side chains by double-mutant cycle analysis and correlated these with the channel activities of corresponding point mutants. The closed central salt bridge has a strong interaction energy of -5.6 kcal mol(-1), which can be broken by forming the open-state salt bridge Glu52-Lys82 (DeltaDeltaG(Inter) = -3.5 kcal mol(-1)) and a weak interaction between Arg138 and Glu128 (DeltaDeltaG(Inter) = -0.6 kcal mol(-1)). A covalent disulfide bond in place of the central salt bridge completely blocks the channel. Growth assays indicate that this gating mechanism could physiologically contribute to the osmoprotection of Escherichia coli cells from environmental stress.     

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.     

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.     

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 determinants of lateral gate opening in the protein translocon

The heterotrimeric SecY/Sec61 complex is a protein-conducting channel that provides a passage for proteins across the membrane as well as a means to integrate nascent proteins into the membrane. While the first function is common among membrane protein channels and transporters, the latter is unique. Insertion of nascent membrane proteins, one transmembrane segment at a time, by SecY likely occurs through a lateral gate in the channel. Molecular dynamics simulations have been used to investigate the mechanism of gate opening. Opening and closing the gate under different conditions allowed us to identify structural elements that resist opening as well as those that aid closure. SecE, considered to act as a clamp keeping the lateral gate closed, was found to play no such role. Loosening of the plug by lateral gate opening, a potential step in channel gating, was also observed. The simulations revealed that lipids on time scales of up to 1 micros do not flood channels with an open lateral gate.     

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.     

Chemically charging the pore constriction opens the mechanosensitive channel MscL

MscL is a bacterial mechanosensitive channel that protects the cell from osmotic downshock. We have previously shown that substitution of a residue that resides within the channel pore constriction, MscL's Gly-22, with all other 19 amino acids affects channel gating according to the hydrophobicity of the substitution (). Here, we first make a mild substitution, G22C, and then attach methanethiosulfonate (MTS) reagents to the cysteine under patch clamp. Binding MTS reagents that are positively charged ([2-(trimethylammonium)ethyl] methanethiosulfonate and 2-aminoethyl methanethiosulfonate) or negatively charged (sodium (2-sulfonatoethyl)methanethiosulfonate) causes MscL to gate spontaneously, even when no tension is applied. In contrast, the polar 2-hydroxyethyl methanethiosulfonate halves the threshold, and the hydrophobic methyl methanethiolsulfonate increases the threshold. These observations indicate that residue 22 is in a hydrophobic environment before gating and in a hydrophilic environment during opening to a substate, a finding consistent with our previous study. In addition, we have found that cysteine 22 is accessible to reagents from the cytoplasmic side only when the channel is opened whereas it is accessible from the periplasmic side even in the closed state. These results support the view that exposure of hydrophobic surfaces to a hydrophilic environment during channel opening serves as the barrier to gating.     

Desensitization mechanism in prokaryotic ligand-gated ion channel

Crystal structures of Gloeobacter violaceus ligand-gated ion channel (GLIC), a proton-gated prokaryotic homologue of pentameric ligand-gated ion channel (LGIC) from G. violaceus, have provided high-resolution models of the channel architecture and its role in selective ion conduction and drug binding. However, it is still unclear which functional states of the LGIC gating scheme these crystal structures represent. Much of this uncertainty arises from a lack of thorough understanding of the functional properties of these prokaryotic channels. To elucidate the molecular events that constitute gating, we have carried out an extensive characterization of GLIC function and dynamics in reconstituted proteoliposomes by patch clamp measurements and EPR spectroscopy. We find that GLIC channels show rapid activation upon jumps to acidic pH followed by a time-dependent loss of conductance because of desensitization. GLIC desensitization is strongly coupled to activation and is modulated by voltage, permeant ions, pore-blocking drugs, and membrane cholesterol. Many of these properties are parallel to functions observed in members of eukaryotic LGIC. Conformational changes in loop C, measured by site-directed spin labeling and EPR spectroscopy, reveal immobilization during desensitization analogous to changes in LGIC and acetylcholine binding protein. Together, our studies suggest conservation of mechanistic aspects of desensitization among LGICs of prokaryotic and eukaryotic origin.     

Mobility of Lower MA-Helices for Ion Conduction through Lateral Portals in 5-HT3A Receptors

The intracellular domain of the serotonin type 3A receptor, a pentameric ligand-gated ion channel, is crucial for regulating conductance. Ion permeation through the extracellular vestibule and the transmembrane channel is well understood, whereas the specific ion conduction pathway through the intracellular domain is less clear. The intracellular domain starts with a short loop after the third transmembrane segment, followed by a short α-helical segment, a large unstructured loop, and finally, the membrane-associated MA-helix that continues into the last transmembrane segment. The MA-helices from all five subunits form the extension of the transmembrane ion channel and shape what has been described as a “closed vestibule,” with their lateral portals obstructed by loops and their cytosolic ends forming a tight hydrophobic constriction. The question remains whether the lateral portals or cytosolic constriction conduct ions upon channel opening. In our study, we used disulfide bond formation between pairs of engineered cysteines to probe the proximity and mobility of segments of the MA-helices most distal to the membrane bilayer. Our results indicate that the proximity and orientation for cysteine pairs at I409C/R410C, in close proximity to the lateral windows, and L402C/L403C, at the cytosolic ends of the MA-helices, are conducive for disulfide bond formation. Although conformational changes associated with gating promote cross-linking for I409C/R410C, which in turn decreases channel currents, cross-linking of L402C/L403C is functionally silent in macroscopic currents. These results support the hypothesis that concerted conformational changes open the lateral portals for ion conduction, rendering ion conduction through the vertical portal unlikely.     

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.     

Locked on one side only: ground state dynamics of the outer membrane efflux duct TolC

Playing a major role in the expulsion of antibiotics and the secretion of cell toxins in conjunction with inner membrane transporters of three protein superfamilies, the outer membrane channel TolC occurs in at least two states blocking or permitting the passage of substrates. The details of the underlying gating mechanism are not fully understood. Addressing the questions of extracellular access control and periplasmic gating mechanism, we conducted a series of independent, unbiased 150-300 ns molecular dynamics simulations of wild-type TolC in a phospholipid membrane/150 mM NaCl water environment. We find that TolC opens and closes freely on the extracellular side, suggesting the absence of a gating mechanism on this side in the isolated protein. On the periplasmic side, we observe the outer periplasmic bottleneck region adopting in all simulations a conformation more open than the TolC wild-type crystal structures until in one run the successive binding of two sodium ions induces the transition to a conformation more closed than any of the available TolC X-ray structures. Concurrent with a heightened sodium residence probability near Asp374, the inner periplasmic bottleneck region at Asp374 remains closed throughout the simulations unless all NaCl is removed from the system, inducing a reopening of the outer and inner bottleneck. Our findings suggest that TolC is locked only on the periplasmic side in a sodium-dependent manner.     

Study,by use of coarse-grained models,of the functionally crucial residues and allosteric pathway of anesthetic regulation of the Gloeobacter violaceus ligand-gated ion channel

Although pentameric ligand-gated ion channels (pLGICs) have been found to be the targets of general anesthetics, the mechanism of the effects of anesthetics on pLGICs remains elusive. pLGICs from Gloeobacter violaceus (GLIC) can be inhibited by the anesthetic ketamine. X-ray crystallography has shown that the ketamine binding site is distant from the channel gate of the GLIC. It is still not clear how ketamine controls the function of the GLIC by long-range allosteric regulation. In this work, the functionally crucial residues and allosteric pathway of anesthetic regulation of the GLIC were identified by use of a coarse-grained thermodynamic method developed by our group. In our method, the functionally crucial sites were identified as the residues thermodynamically coupled with binding of ketamine. The results from calculation were highly consistent with experimental data. Our study aids understanding of the mechanism of the anesthetic action of ketamine on the GLIC by long-range allosteric modulation.