Structure of a potentially open state of a proton-activated pentameric ligand-gated ion channel

Conformational changes in acetylcholine binding protein investigated by temperature accelerated molecular dynamics

Despite the large number of studies available on nicotinic acetylcholine receptors, a complete account of the mechanistic aspects of their gating transition in response to ligand binding still remains elusive. As a first step toward dissecting the transition mechanism by accelerated sampling techniques, we study the ligand-induced conformational changes of the acetylcholine binding protein (AChBP), a widely accepted model for the full receptor extracellular domain. Using unbiased Molecular Dynamics (MD) and Temperature Accelerated Molecular Dynamics (TAMD) simulations we investigate the AChBP transition between the apo and the agonist-bound state. In long standard MD simulations, both conformations of the native protein are stable, while the agonist-bound structure evolves toward the apo one if the orientation of few key sidechains in the orthosteric cavity is modified. Conversely, TAMD simulations initiated from the native conformations are able to produce the spontaneous transition. With respect to the modified conformations, TAMD accelerates the transition by at least a factor 10. The analysis of some specific residue-residue interactions points out that the transition mechanism is based on the disruption/formation of few key hydrogen bonds. Finally, while early events of ligand dissociation are observed already in standard MD, TAMD accelerates the ligand detachment and, at the highest TAMD effective temperature, it is able to produce a complete dissociation path in one AChBP subunit.

Seeking structural specificity: direct modulation of pentameric ligand-gated ion channels by alcohols and general anesthetics

Alcohols and other anesthetic agents dramatically alter neurologic function in a wide range of organisms, yet their molecular sites of action remain poorly characterized. Pentameric ligand-gated ion channels, long implicated in important direct effects of alcohol and anesthetic binding, have recently been illuminated in renewed detail thanks to the determination of atomic-resolution structures of several family members from lower organisms. These structures provide valuable models for understanding and developing anesthetic agents and for allosteric modulation in general. This review surveys progress in this field from function to structure and back again, outlining early evidence for relevant modulation of pentameric ligand-gated ion channels and the development of early structural models for ion channel function and modulation. We highlight insights and challenges provided by recent crystal structures and resulting simulations, as well as opportunities for translation of these newly detailed models back to behavior and therapy.

Insights into the binding of GABA to the insect RDL receptor from atomistic simulations: a comparison of models

The resistance to dieldrin (RDL) receptor is an insect pentameric ligand-gated ion channel (pLGIC). It is activated by the neurotransmitter γ-aminobutyric acid (GABA) binding to its extracellular domain; hence elucidating the atomistic details of this interaction is important for understanding how the RDL receptor functions. As no high resolution structures are currently available, we built homology models of the extracellular domain of the RDL receptor using different templates, including the widely used acetylcholine binding protein and two pLGICs, the Erwinia Chrysanthemi ligand-gated ion channel (ELIC) and the more recently resolved GluCl. We then docked GABA into the selected three dimensional structures, which we used as starting points for classical molecular dynamics simulations. This allowed us to analyze in detail the behavior of GABA in the binding sites, including the hydrogen bond and cation-π interaction networks it formed, the conformers it visited and the possible role of water molecules in mediating the interactions; we also estimated the binding free energies. The models were all stable and showed common features, including interactions consistent with experimental data and similar to other pLGICs; differences could be attributed to the quality of the models, which increases with increasing sequence identity, and the use of a pLGIC template. We supplemented the molecular dynamics information with metadynamics, a rare event method, by exploring the free energy landscape of GABA binding to the RDL receptor. Overall, we show that the GluCl template provided the best models. GABA forming direct salt-bridges with Arg211 and Glu204, and cation-π interactions with an aromatic cage including Tyr109, Phe206 and Tyr254, represents a favorable binding arrangement, and the interaction with Glu204 can also be mediated by a water molecule.

The impact of human hyperekplexia mutations on glycine receptor structure and function

Hyperekplexia is a rare neurological disorder characterized by neonatal hypertonia, exaggerated startle responses to unexpected stimuli and a variable incidence of apnoea, intellectual disability and delays in speech acquisition. The majority of motor defects are successfully treated by clonazepam. Hyperekplexia is caused by hereditary mutations that disrupt the functioning of inhibitory glycinergic synapses in neuromotor pathways of the spinal cord and brainstem. The human glycine receptor α1 and β subunits, which predominate at these synapses, are the major targets of mutations. International genetic screening programs, that together have analysed several hundred probands, have recently generated a clear picture of genotype-phenotype correlations and the prevalence of different categories of hyperekplexia mutations. Focusing largely on this new information, this review seeks to summarise the effects of mutations on glycine receptor structure and function and how these functional alterations lead to hyperekplexia.

Allosteric regulation of pentameric ligand-gated ion channels: an emerging mechanistic perspective

Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communications in the nervous system by converting the binding of a chemical messenger—a neurotransmitter—into an ion flux through the postsynaptic membrane. They are oligomeric assemblies that provide prototypical examples of allosterically regulated integral membrane proteins. Here, we present an overview of the most recent advances on the signal transduction mechanism based on the X-ray structures of both prokaryotic and invertebrate eukaryotic pLGICs and atomistic Molecular Dynamics simulations. The present results suggest that ion gating involves a large structural reorganization of the molecule mediated by two distinct quaternary transitions, a global twisting and the blooming of the extracellular domain, which can be modulated by ligand binding at the topographically distinct orthosteric and allosteric sites. The emerging model of gating is consistent with a wealth of functional studies and will boost the development of novel pharmacological strategies.

NMR structures of the human α7 nAChR transmembrane domain and associated anesthetic binding sites

The α7 nicotinic acetylcholine receptor (nAChR), assembled as homomeric pentameric ligand-gated ion channels, is one of the most abundant nAChR subtypes in the brain. Despite its importance in memory, learning and cognition, no structure has been determined for the α7 nAChR TM domain, a target for allosteric modulators. Using solution state NMR, we determined the structure of the human α7 nAChR TM domain (PDB ID: 2MAW) and demonstrated that the α7 TM domain formed functional channels in Xenopus oocytes. We identified the associated binding sites for the anesthetics halothane and ketamine; the former cannot sensitively inhibit α7 function, but the latter can. The α7 TM domain folds into the expected four-helical bundle motif, but the intra-subunit cavity at the extracellular end of the α7 TM domain is smaller than the equivalent cavity in the α4β2 nAChRs (PDB IDs: 2LLY; 2LM2). Neither drug binds to the extracellular end of the α7 TM domain, but two halothane molecules or one ketamine molecule binds to the intracellular end of the α7 TM domain. Halothane and ketamine binding sites are partially overlapped. Ketamine, but not halothane, perturbed the α7 channel-gate residue L9'. Furthermore, halothane did not induce profound dynamics changes in the α7 channel as observed in α4β2. The study offers a novel high-resolution structure for the human α7 nAChR TM domain that is invaluable for developing α7-specific therapeutics. It also provides evidence to support the hypothesis: only when anesthetic binding perturbs the channel pore or alters the channel motion, can binding generate functional consequences.

Conduits of life's spark: a perspective on ion channel research since the birth of neuron

Heartbeats, muscle twitches, and lightning-fast thoughts are all manifestations of bioelectricity and rely on the activity of a class of membrane proteins known as ion channels. The basic function of an ion channel can be distilled into, "The hole opens. Ions go through. The hole closes." Studies of the fundamental mechanisms by which this process happens and the consequences of such activity in the setting of excitable cells remains the central focus of much of the field. One might wonder after so many years of detailed poking at such a seemingly simple process, is there anything left to learn?

Intra-subunit flexibility underlies activation and allosteric modulation of neuronal nicotinic acetylcholine receptors

Allosteric modulation is a general feature of nicotinic acetylcholine receptors, yet the structural components and movements important for conversions among functional states are not well understood. In this study, we examine the communication between the binding sites for agonist and the modulator morantel (Mor) of neuronal α3β2 receptors, measuring evoked currents of receptors expressed in Xenopus oocytes with the two-electrode voltage-clamp method. We hypothesized that movement along an interface of β sheets connecting the agonist and modulator sites is necessary for allosteric modulation. To address this, we created pairs of substituted cysteines that span the cleft formed where the outer β sheet meets the β sheet constituting the (-)-face of the α3 subunit; the three pairs were L158C-A179C, L158C-G181C and L158C-K183C. Employing a disulfide trapping approach in which bonds are formed between neighboring cysteines under oxidation conditions, we found that oxidation treatments decreased the amplitude of currents evoked by either the agonist (ACh) or co-applied agonist and modulator (ACh + Mor), by as much as 51%, consistent with the introduced bond decreasing channel efficacy. Reduction treatment increased evoked currents up to 89%. The magnitude of the oxidation effects depended on whether agonists were present during oxidation and on the cysteine pair. Additionally, the cysteine mutations themselves decreased Mor potentiation, implicating these residues in modulation. Our findings suggest that these β sheets in the α3 subunit move with respect to each other during activation and modulation, and the residues studied highlight the contribution of this intramolecular allosteric pathway to receptor function.

Crystal structures of a pentameric ligand-gated ion channel provide a mechanism for activation

Pentameric ligand-gated ion channels mediate fast chemical transmission of nerve signals. The structure of a bacterial proton-gated homolog has been established in its open and locally closed conformations at acidic pH. Here we report its crystal structure at neutral pH, thereby providing the X-ray structures of the two end-points of the gating mechanism in the same pentameric ligand-gated ion channel. The large structural variability in the neutral pH structure observed in the four copies of the pentamer present in the asymmetric unit has been used to analyze the intrinsic fluctuations in this state, which are found to prefigure the transition to the open state. In the extracellular domain (ECD), a marked quaternary change is observed, involving both a twist and a blooming motion, and the pore in the transmembrane domain (TMD) is closed by an upper bend of helix M2 (as in locally closed form) and a kink of helix M1, both helices no longer interacting across adjacent subunits. On the tertiary level, detachment of inner and outer β sheets in the ECD reshapes two essential cavities at the ECD-ECD and ECD-TMD interfaces. The first one is the ligand-binding cavity; the other is close to a known divalent cation binding site in other pentameric ligand-gated ion channels. In addition, a different crystal form reveals that the locally closed and open conformations coexist as discrete ones at acidic pH. These structural results, together with site-directed mutagenesis, physiological recordings, and coarse-grained modeling, have been integrated to propose a model of the gating transition pathway.

Structural basis for allosteric coupling at the membrane-protein interface in Gloeobacter violaceus ligand-gated ion channel (GLIC)

Ligand binding at the extracellular domain of pentameric ligand-gated ion channels initiates a relay of conformational changes that culminates at the gate within the transmembrane domain. The interface between the two domains is a key structural entity that governs gating. Molecular events in signal transduction at the interface are poorly defined because of its intrinsically dynamic nature combined with functional modulation by membrane lipid and water vestibules. Here we used electron paramagnetic resonance spectroscopy to delineate protein motions underlying Gloeobacter violaceus ligand-gated ion channel gating in a membrane environment and report the interface conformation in the closed and the desensitized states. Extensive intrasubunit interactions were observed in the closed state that are weakened upon desensitization and replaced by newer intersubunit contacts. Gating involves major rearrangements of the interfacial loops, accompanied by reorganization of the protein-lipid-water interface. These structural changes may serve as targets for modulation of gating by lipids, alcohols, and amphipathic drug molecules.

Insights into distinct modulation of α7 and α7β2 nicotinic acetylcholine receptors by the volatile anesthetic isoflurane

Nicotinic acetylcholine receptors (nAChRs) are targets of general anesthetics, but functional sensitivity to anesthetic inhibition varies dramatically among different subtypes of nAChRs. Potential causes underlying different functional responses to anesthetics remain elusive. Here we show that in contrast to the α7 nAChR, the α7β2 nAChR is highly susceptible to inhibition by the volatile anesthetic isoflurane in electrophysiology measurements. Isoflurane-binding sites in β2 and α7 were found at the extracellular and intracellular end of their respective transmembrane domains using NMR. Functional relevance of the identified β2 site was validated via point mutations and subsequent functional measurements. Consistent with their functional responses to isoflurane, β2 but not α7 showed pronounced dynamics changes, particularly for the channel gate residue Leu-249(9'). These results suggest that anesthetic binding alone is not sufficient to generate functional impact; only those sites that can modulate channel dynamics upon anesthetic binding will produce functional effects.

Macroscopic kinetics of pentameric ligand gated ion channels: comparisons between two prokaryotic channels and one eukaryotic channel

Electrochemical signaling in the brain depends on pentameric ligand-gated ion channels (pLGICs). Recently, crystal structures of prokaryotic pLGIC homologues from Erwinia chrysanthemi (ELIC) and Gloeobacter violaceus (GLIC) in presumed closed and open channel states have been solved, which provide insight into the structural mechanisms underlying channel activation. Although structural studies involving both ELIC and GLIC have become numerous, thorough functional characterizations of these channels are still needed to establish a reliable foundation for comparing kinetic properties. Here, we examined the kinetics of ELIC and GLIC current activation, desensitization, and deactivation and compared them to the GABA_A receptor, a prototypic eukaryotic pLGIC. Outside-out patch-clamp recordings were performed with HEK-293T cells expressing ELIC, GLIC, or α₁β₂γ_2L GABA_A receptors, and ultra-fast ligand application was used. In response to saturating agonist concentrations, we found both ELIC and GLIC current activation were two to three orders of magnitude slower than GABA_A receptor current activation. The prokaryotic channels also had slower current desensitization on a timescale of seconds. ELIC and GLIC current deactivation following 25 s pulses of agonist (cysteamine and pH 4.0 buffer, respectively) were relatively fast with time constants of 24.9±5.1 ms and 1.2±0.2 ms, respectively. Surprisingly, ELIC currents evoked by GABA activated very slowly with a time constant of 1.3±0.3 s and deactivated even slower with a time constant of 4.6±1.2 s. We conclude that the prokaryotic pLGICs undergo similar agonist-mediated gating transitions to open and desensitized states as eukaryotic pLGICs, supporting their use as experimental models. Their uncharacteristic slow activation, slow desensitization and rapid deactivation time courses are likely due to differences in specific structural elements, whose future identification may help uncover mechanisms underlying pLGIC gating transitions.

Site-directed spin labeling reveals pentameric ligand-gated ion channel gating motions

Pentameric ligand-gated ion channels (pLGICs) are neurotransmitter-activated receptors that mediate fast synaptic transmission. In pLGICs, binding of agonist to the extracellular domain triggers a structural rearrangement that leads to the opening of an ion-conducting pore in the transmembrane domain and, in the continued presence of neurotransmitter, the channels desensitize (close). The flexible loops in each subunit that connect the extracellular binding domain (loops 2, 7, and 9) to the transmembrane channel domain (M2–M3 loop) are essential for coupling ligand binding to channel gating. Comparing the crystal structures of two bacterial pLGIC homologues, ELIC and the proton-activated GLIC, suggests channel gating is associated with rearrangements in these loops, but whether these motions accurately predict the motions in functional lipid-embedded pLGICs is unknown. Here, using site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy and functional GLIC channels reconstituted into liposomes, we examined if, and how far, the loops at the ECD/TMD gating interface move during proton-dependent gating transitions from the resting to desensitized state. Loop 9 moves ∼9 Å inward toward the channel lumen in response to proton-induced desensitization. Loop 9 motions were not observed when GLIC was in detergent micelles, suggesting detergent solubilization traps the protein in a nonactivatable state and lipids are required for functional gating transitions. Proton-induced desensitization immobilizes loop 2 with little change in position. Proton-induced motion of the M2–M3 loop was not observed, suggesting its conformation is nearly identical in closed and desensitized states. Our experimentally derived distance measurements of spin-labeled GLIC suggest ELIC is not a good model for the functional resting state of GLIC, and that the crystal structure of GLIC does not correspond to a desensitized state. These findings advance our understanding of the molecular mechanisms underlying pLGIC gating.

Ligand-gated ion channels reside in the membranes of nerve and muscle cells. These proteins form channels that span the membrane, where they transduce chemical signals into changes in electrical excitability. Neurotransmitters bind to the extracellular surface of these proteins to trigger global structural rearrangements that open the channel, allowing ions to flow across the cell membrane. In the continued presence of neurotransmitters, the channels desensitize and close. Channel opening and closing regulate muscle contraction and signaling in the brain, and defects in these channels lead to a variety of diseases. While crystal structures have provided frozen snapshots of these proteins in presumed closed and open channel states, little is known about how the channels desensitize and move during actual signaling events. Here, we applied a technique to investigate the structure and local dynamics of proteins known as site-directed spin labeling to a prototypical ligand-gated channel, GLIC. We directly quantified ligand-induced motions in regions at the boundary between the binding domain (loops 2 and 9) and the channel domain (M2–M3 loop). We show that a large movement of loop 9 and an immobilization of loop 2, which rearranges the interface between the binding and channel domains, accompanies GLIC channel gating transitions into a desensitized state. These data provide new insights into the protein movements that underlie electrochemical transmission of signals between cells.

Qualitative computational bioanalytics: assembly of viral channel-forming peptides around mono and divalent ions

A fine-grained docking protocol was used to generate a bundle-like structure of the bitopic membrane protein Vpu from HIV-1. Vpu is a type I membrane protein with 81 amino acids. It is proposed that Vpu forms ion- and substrate-conducting bundles, which are located at the plasma membrane in the infected cell. The Vpu1-32 peptide that includes the transmembrane domain (TMD) is assembled into homo-pentameric bundles around prepositioned Na, K, Ca or Cl ions. For bundles with the lowest energy, the TMDs generate a hydrophobic pore. Bundles in which Ser-24 faces the pore have higher energy. The tilt of the helices in the lowest energy bundles is larger than bundles with serines facing the pore. Left-handed bundles are lowest in energy where the ions are located at the serines.

Gating of the proton-gated ion channel from Gloeobacter violaceus at pH 4 as revealed by X-ray crystallography

Cryoelectron microscopy and X-ray crystallography have recently been used to generate structural models that likely represent the unliganded closed-channel conformation and the fully liganded open-channel conformation of different members of the nicotinic-receptor superfamily. To characterize the structure of the closed-channel conformation in its liganded state, we identified a number of positions in the loop between transmembrane segments 2 (M2) and 3 (M3) of a proton-gated ortholog from the bacterium Gloeobacter violaceus (GLIC) where mutations to alanine reduce the liganded-gating equilibrium constant, and solved the crystal structures of two such mutants (T25'A and Y27'A) at pH ~4.0. At the level of backbone atoms, the liganded closed-channel model presented here differs from the liganded open-channel structure of GLIC in the pre-M1 linker, the M3-M4 loop, and much more prominently, in the extracellular half of the pore lining, where the more pronounced tilt of the closed-channel M2 α-helices toward the pore's long axis narrows the permeation pathway. On the other hand, no differences between the liganded closed-channel and open-channel models could be detected at the level of the extracellular domain, where conformational changes are expected to underlie the low-to-high proton-affinity switch that drives gating of proton-bound channels. Thus, the liganded closed-channel model is nearly indistinguishable from the recently described "locally closed" structure. However, because cross-linking strategies (which could have stabilized unstable conformations) and mutations involving ionizable side chains (which could have affected proton-gated channel activation) were purposely avoided, we favor the notion that this structure represents one of the end states of liganded gating rather than an unstable intermediate.

Structural basis for potentiation by alcohols and anaesthetics in a ligand-gated ion channel

Ethanol alters nerve signalling by interacting with proteins in the central nervous system, particularly pentameric ligand-gated ion channels. A recent series of mutagenesis experiments on Gloeobacter violaceus ligand-gated ion channel, a prokaryotic member of this family, identified a single-site variant that is potentiated by pharmacologically relevant concentrations of ethanol. Here we determine crystal structures of the ethanol-sensitized variant in the absence and presence of ethanol and related modulators, which bind in a transmembrane cavity between channel subunits and may stabilize the open form of the channel. Structural and mutagenesis studies defined overlapping mechanisms of potentiation by alcohols and anaesthetics via the inter-subunit cavity. Furthermore, homology modelling show this cavity to be conserved in human ethanol-sensitive glycine and GABA(A) receptors, and to involve residues previously shown to influence alcohol and anaesthetic action on these proteins. These results suggest a common structural basis for ethanol potentiation of an important class of targets for neurological actions of ethanol.

Analysis of hyperekplexia mutations identifies transmembrane domain rearrangements that mediate glycine receptor activation

Pentameric ligand-gated ion channels (pLGICs) mediate numerous physiological processes and are therapeutic targets for a wide range of clinical indications. Elucidating the structural differences between their closed and open states may help in designing improved drugs that bias receptors toward the desired conformational state. We recently showed that two new hyperekplexia mutations, Q226E and V280M, induced spontaneous activity in α1 glycine receptors. Gln-226, located near the top of transmembrane (TM) 1, is closely apposed to Arg-271 at the top of TM2 in the neighboring subunit. Using mutant cycle analysis, we inferred that Q226E induces activation via an enhanced electrostatic attraction to Arg-271. This would tilt the top of TM2 toward TM1 and hence away from the pore axis to open the channel. We also concluded that the increased side chain volume of V280M, in the TM2-TM3 loop, exerts a steric repulsion against Ile-225 at the top of TM1 in the neighboring subunit. We infer that this steric repulsion would tilt the top of TM3 radially outwards against the stationary TM1 and thus provide space for TM2 to relax away from the pore axis to create an open channel. Because the transmembrane domain movements inferred from this functional analysis are consistent with the structural differences evident in the x-ray atomic structures of closed and open state bacterial pLGICs, we propyose that the model of pLGIC activation as outlined here may be broadly applicable across the eukaryotic pLGIC receptor family.

Identification of a possible secondary picrotoxin-binding site on the GABA(A) receptor

The type A GABA receptors (GABARs) are ligand-gated ion channels (LGICs) found in the brain and are the major inhibitory neurotransmitter receptors. Upon binding of an agonist, the GABAR opens and increases the intraneuronal concentration of chloride ions, thus hyperpolarizing the cell and inhibiting the transmission of the nerve action potential. GABARs also contain many other modulatory binding pockets that differ from the agonist-binding site. The composition of the GABAR subunits can alter the properties of these modulatory sites. Picrotoxin is a noncompetitive antagonist for LGICs, and by inhibiting GABAR, picrotoxin can cause overstimulation and induce convulsions. We use addition of picrotoxin to probe the characteristics and possible mechanism of an additional modulatory pocket located at the interface between the ligand-binding domain and the transmembrane domain of the GABAR. Picrotoxin is widely regarded as a pore-blocking agent that acts at the cytoplasmic end of the channel. However, there are also data to suggest that there may be an additional, secondary binding site for picrotoxin. Through homology modeling, molecular docking, and molecular dynamics simulations, we show that binding of picrotoxin to this interface pocket correlates with these data, and negative modulation occurs at the pocket via a kinking of the pore-lining helices into a more closed orientation.

Identification of a probable pore-forming domain in the multimeric vacuolar anion channel AtALMT9

Aluminum-activated malate transporters (ALMTs) form an important family of anion channels involved in fundamental physiological processes in plants. Because of their importance, the role of ALMTs in plant physiology is studied extensively. In contrast, the structural basis of their functional properties is largely unknown. This lack of information limits the understanding of the functional and physiological differences between ALMTs and their impact on anion transport in plants. This study aimed at investigating the structural organization of the transmembrane domain of the Arabidopsis (Arabidopsis thaliana) vacuolar channel AtALMT9. For that purpose, we performed a large-scale mutagenesis analysis and found two residues that form a salt bridge between the first and second putative transmembrane α-helices (TMα1 and TMα2). Furthermore, using a combination of pharmacological and mutagenesis approaches, we identified citrate as an "open channel blocker" of AtALMT9 and used this tool to examine the inhibition sensitivity of different point mutants of highly conserved amino acid residues. By this means, we found a stretch within the cytosolic moiety of the TMα5 that is a probable pore-forming domain. Moreover, using a citrate-insensitive AtALMT9 mutant and biochemical approaches, we could demonstrate that AtALMT9 forms a multimeric complex that is supposedly composed of four subunits. In summary, our data provide, to our knowledge, the first evidence about the structural organization of an ion channel of the ALMT family. We suggest that AtALMT9 is a tetramer and that the TMα5 domains of the subunits contribute to form the pore of this anion channel.

Intermediate closed state for glycine receptor function revealed by cysteine cross-linking

Pentameric ligand-gated ion channels (pLGICs) mediate signal transmission by coupling the binding of extracellular ligands to the opening of their ion channel. Agonist binding elicits activation and desensitization of pLGICs, through several conformational states, that are, thus far, incompletely characterized at the structural level. We previously reported for GLIC, a prokaryotic pLGIC, that cross-linking of a pair of cysteines at both sides of the extracellular and transmembrane domain interface stabilizes a locally closed (LC) X-ray structure. Here, we introduced the homologous pair of cysteines on the human α1 glycine receptor. We show by electrophysiology that cysteine cross-linking produces a gain-of-function phenotype characterized by concomitant constitutive openings, increased agonist potency, and equalization of efficacies of full and partial agonists. However, it also produces a reduction of maximal currents at saturating agonist concentrations without change of the unitary channel conductance, an effect reversed by the positive allosteric modulator propofol. The cross-linking thus favors a unique closed state distinct from the resting and longest-lived desensitized states. Fitting the data according to a three-state allosteric model suggests that it could correspond to a LC conformation. Its plausible assignment to a gating intermediate or a fast-desensitized state is discussed. Overall, our data show that relative movement of two loops at the extracellular-transmembrane interface accompanies orthosteric agonist-mediated gating.

Nicotinic acetylcholine receptor and the structural basis of neuromuscular transmission: insights from Torpedo postsynaptic membranes

The nicotinic acetylcholine (ACh) receptor, at the neuromuscular junction, is a neurotransmitter-gated ion channel that has been fine-tuned through evolution to transduce a chemical signal into an electrical signal with maximum efficiency and speed. It is composed from three similar and two identical polypeptide chains, arranged in a ring around a narrow membrane pore. Central to the design of this assembly is a hydrophobic gate in the pore, more than 50 Å away from sites in the extracellular domain where ACh binds. Although the molecular properties of the receptor have been explored intensively over the last few decades, only recently have structures emerged revealing its complex architecture and illuminating how ACh entering the binding sites opens the distant gate. Postsynaptic membranes isolated from the (muscle-derived) electric organ of the Torpedo ray have underpinned most of the structural studies: the membranes form tubular vesicles having receptors arranged on a regular surface lattice, which can be imaged directly in frozen physiological solutions. Advances in electron crystallographic techniques have also been important, enabling analysis of the closed- and open-channel forms of the receptor in unreacted tubes or tubes reacted briefly with ACh. The structural differences between these two forms show that all five subunits participate in a concerted conformational change communicating the effect of ACh binding to the gate, but that three of them (α_γ, β and δ) play a dominant role. Flexing of oppositely facing pore-lining α-helices is the principal motion determining the closed/open state of the gate. These results together with the findings of biochemical, biophysical and other structural studies allow an integrated description of the receptor and of its mode of action at the synapse.

A gating mechanism of pentameric ligand-gated ion channels

Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communication in the nervous system and are involved in fundamental processes such as attention, learning, and memory. They are oligomeric protein assemblies that convert a chemical signal into an ion flux through the postsynaptic membrane, but the molecular mechanism of gating ions has remained elusive. Here, we present atomistic molecular dynamics simulations of the prokaryotic channels from Gloeobacter violaceus (GLIC) and Erwinia chrysanthemi (ELIC), whose crystal structures are thought to represent the active and the resting states of pLGICs, respectively, and of the eukaryotic glutamate-gated chloride channel from Caenorhabditis elegans (GluCl), whose open-channel structure was determined complexed with the positive allosteric modulator ivermectin. Structural observables extracted from the trajectories of GLIC and ELIC are used as progress variables to analyze the time evolution of GluCl, which was simulated in the absence of ivermectin starting from the structure with bound ivermectin. The trajectory of GluCl with ivermectin removed shows a sequence of structural events that couple agonist unbinding from the extracellular domain to ion-pore closing in the transmembrane domain. Based on these results, we propose a structural mechanism for the allosteric communication leading to deactivation/activation of the GluCl channel. This model of gating emphasizes the coupling between the quaternary twisting and the opening/closing of the ion pore and is likely to apply to other members of the pLGIC family.

Cysteine substitutions define etomidate binding and gating linkages in the α-M1 domain of γ-aminobutyric acid type A (GABAA) receptors

Etomidate is a potent general anesthetic that acts as an allosteric co-agonist at GABAA receptors. Photoreactive etomidate derivatives labeled αMet-236 in transmembrane domain M1, which structural models locate in the β+/α- subunit interface. Other nearby residues may also contribute to etomidate binding and/or transduction through rearrangement of the site. In human α1β2γ2L GABAA receptors, we applied the substituted cysteine accessibility method to α1-M1 domain residues extending from α1Gln-229 to α1Gln-242. We used electrophysiology to characterize each mutant's sensitivity to GABA and etomidate. We also measured rates of sulfhydryl modification by p-chloromercuribenzenesulfonate (pCMBS) with and without GABA and tested if etomidate blocks modification of pCMBS-accessible cysteines. Cys substitutions in the outer α1-M1 domain impaired GABA activation and variably affected etomidate sensitivity. In seven of eight residues where pCMBS modification was evident, rates of modification were accelerated by GABA co-application, indicating that channel activation increases water and/or pCMBS access. Etomidate reduced the rate of modification for cysteine substitutions at α1Met-236, α1Leu-232 and α1Thr-237. We infer that these residues, predicted to face β2-M3 or M2 domains, contribute to etomidate binding. Thus, etomidate interacts with a short segment of the outer α1-M1 helix within a subdomain that undergoes significant structural rearrangement during channel gating. Our results are consistent with in silico docking calculations in a homology model that orient the long axis of etomidate approximately orthogonal to the transmembrane axis.

Open-channel structures of the human glycine receptor α1 full-length transmembrane domain

Glycine receptors play a major role in mediating fast inhibitory neurotransmission in the spinal cord and brain stem, yet their high-resolution structures remain unsolved. We determined open-channel structures of the full-length transmembrane domain (TMD) of the human glycine receptor α1-subunit (hGlyR-α1) using nuclear magnetic resonance (NMR) spectroscopy and electron micrographs. hGlyR-α1 TMD spontaneously forms pentameric Cl(-)-conducting channels, with structures sharing overall topology observed in crystal structures of homologous bacterial and nematode pentameric ligand-gated ion channels (pLGICs). However, the mammalian hGlyR-α1 structures present several distinctive features, including a shorter, pore-lining TM2 helix with helical unwinding near the C-terminal end, a TM3 helical kink at A288 that partially overlaps with the homologous ivermectin-binding site in GluCl, and a highly dynamic segment between S267(15') of TM2 and A288 that likely affects allosteric modulations of channel function. Our structures provide additional templates for identifying potential drug targets in GlyRs and other mammalian pLGICs.

Functional validation of virtual screening for novel agents with general anesthetic action at ligand-gated ion channels

GABA(A) receptors play a crucial role in the actions of general anesthetics. The recently published crystal structure of the general anesthetic propofol bound to Gloeobacter violaceus ligand-gated ion channel (GLIC), a bacterial homolog of GABA(A) receptors, provided an opportunity to explore structure-based ligand discovery for pentameric ligand-gated ion channels (pLGICs). We used molecular docking of 153,000 commercially available compounds to identify molecules that interact with the propofol binding site in GLIC. In total, 29 compounds were selected for functional testing on recombinant GLIC, and 16 of these compounds modulated GLIC function. Active compounds were also tested on recombinant GABA(A) receptors, and point mutations around the presumed binding pocket were introduced into GLIC and GABA(A) receptors to test for binding specificity. The potency of active compounds was only weakly correlated with properties such as lipophilicity or molecular weight. One compound was found to mimic the actions of propofol on GLIC and GABA(A), and to be sensitive to mutations that reduce the action of propofol in both receptors. Mutant receptors also provided insight about the position of the binding sites and the relevance of the receptor's conformation for anesthetic actions. Overall, the findings support the feasibility of the use of virtual screening to discover allosteric modulators of pLGICs, and suggest that GLIC is a valid model system to identify novel GABA(A) receptor ligands.

Nicotinic receptor transduction zone: invariant arginine couples to multiple electron-rich residues

Gating of the muscle-type acetylcholine receptor (AChR) channel depends on communication between the ACh-binding site and the remote ion channel. A key region for this communication is located within the structural transition zone between the ligand-binding and pore domains. Here, stemming from β-strand 10 of the binding domain, the invariant αArg209 lodges within the hydrophobic interior of the subunit and is essential for rapid and efficient channel gating. Previous charge-reversal experiments showed that the contribution of αArg209 to channel gating depends strongly on αGlu45, also within this region. Here we determine whether the contribution of αArg209 to channel gating depends on additional anionic or electron-rich residues in this region. Also, to reconcile diverging findings in the literature, we compare the dependence of αArg209 on αGlu45 in AChRs from different species, and compare the full agonist ACh with the weak agonist choline. Our findings reveal that the contribution of αArg209 to channel gating depends on additional nearby electron-rich residues, consistent with both electrostatic and steric contributions. Furthermore, αArg209 and αGlu45 show a strong interdependence in both human and mouse AChRs, whereas the functional consequences of the mutation αE45R depend on the agonist. The emerging picture shows a multifaceted network of interdependent residues that are required for communication between the ligand-binding and pore domains.

Covalent trapping of methyllycaconitine at the α4-α4 interface of the α4β2 nicotinic acetylcholine receptor: antagonist binding site and mode of receptor inhibition revealed

The α4β2 nicotinic acetylcholine receptors (nAChRs) are widely expressed in the brain and are implicated in a variety of physiological processes. There are two stoichiometries of the α4β2 nAChR, (α4)2(β2)3 and (α4)3(β2)2, with different sensitivities to acetylcholine (ACh), but their pharmacological profiles are not fully understood. Methyllycaconitine (MLA) is known to be an antagonist of nAChRs. Using the two-electrode voltage clamp technique and α4β2 nAChRs in the Xenopus oocyte expression system, we demonstrate that inhibition by MLA occurs via two different mechanisms; that is, a direct competitive antagonism and an apparently insurmountable mechanism that only occurs after preincubation with MLA. We hypothesized an additional MLA binding site in the α4-α4 interface that is unique to this stoichiometry. To prove this, we covalently trapped a cysteine-reactive MLA analog at an α4β2 receptor containing an α4(D204C) mutation predicted by homology modeling to be within reach of the reactive probe. We demonstrate that covalent trapping results in irreversible reduction of ACh-elicited currents in the (α4)3(β2)2 stoichiometry, indicating that MLA binds to the α4-α4 interface of the (α4)3(β2)2 and providing direct evidence of ligand binding to the α4-α4 interface. Consistent with other studies, we propose that the α4-α4 interface is a structural target for potential therapeutics that modulate (α4)3(β2)2 nAChRs.

Specificity of intersubunit general anesthetic-binding sites in the transmembrane domain of the human α1β3γ2 γ-aminobutyric acid type A (GABAA) receptor

GABA type A receptors (GABAAR), the brain's major inhibitory neurotransmitter receptors, are the targets for many general anesthetics, including volatile anesthetics, etomidate, propofol, and barbiturates. How such structurally diverse agents can act similarly as positive allosteric modulators of GABAARs remains unclear. Previously, photoreactive etomidate analogs identified two equivalent anesthetic-binding sites in the transmembrane domain at the β(+)-α(-) subunit interfaces, which also contain the GABA-binding sites in the extracellular domain. Here, we used R-[(3)H]5-allyl-1-methyl-5-(m-trifluoromethyl-diazirynylphenyl) barbituric acid (R-mTFD-MPAB), a potent stereospecific barbiturate anesthetic, to photolabel expressed human α1β3γ2 GABAARs. Protein microsequencing revealed that R-[(3)H]mTFD-MPAB did not photolabel the etomidate sites at the β(+)-α(-) subunit interfaces. Instead, it photolabeled sites at the α(+)-β(-) and γ(+)-β(-) subunit interfaces in the transmembrane domain. On the (+)-side, α1M3 was labeled at Ala-291 and Tyr-294 and γ2M3 at Ser-301, and on the (-)-side, β3M1 was labeled at Met-227. These residues, like those in the etomidate site, are located at subunit interfaces near the synaptic side of the transmembrane domain. The selectivity of R-etomidate for the β(+)-α(-) interface relative to the α(+)-β(-)/γ(+)-β(-) interfaces was >100-fold, whereas that of R-mTFD-MPAB for its sites was >50-fold. Each ligand could enhance photoincorporation of the other, demonstrating allosteric interactions between the sites. The structural heterogeneity of barbiturate, etomidate, and propofol derivatives is accommodated by varying selectivities for these two classes of sites. We hypothesize that binding at any of these homologous intersubunit sites is sufficient for anesthetic action and that this explains to some degree the puzzling structural heterogeneity of anesthetics.

Structural and functional characterization of an anesthetic binding site in the second cysteine-rich domain of protein kinase Cδ*

Elucidating the principles governing anesthetic-protein interactions requires structural determinations at high resolutions not yet achieved with ion channels. Protein kinase C (PKC) activity is modulated by general anesthetics. We solved the structure of the phorbol-binding domain (C1B) of PKCδ complexed with an ether (methoxymethylcycloprane) and with an alcohol (cyclopropylmethanol) at 1.36-Å resolution. The cyclopropane rings of both agents displace a single water molecule in a surface pocket adjacent to the phorbol-binding site, making van der Waals contacts with the backbone and/or side chains of residues Asn-237 to Ser-240. Surprisingly, two water molecules anchored in a hydrogen-bonded chain between Thr-242 and Lys-260 impart elasticity to one side of the binding pocket. The cyclopropane ring takes part in π-acceptor hydrogen bonds with the amide of Met-239. There is a crucial hydrogen bond between the oxygen atoms of the anesthetics and the hydroxyl of Tyr-236. A Tyr-236-Phe mutation results in loss of binding. Thus, both van der Waals interactions and hydrogen-bonding are essential for binding to occur. Ethanol failed to bind because it is too short to benefit from both interactions. Cyclopropylmethanol inhibited phorbol-ester-induced PKCδ activity, but failed to do so in PKCδ containing the Tyr-236-Phe mutation.

Propofol binding to the resting state of the gloeobacter violaceus ligand-gated ion channel (GLIC) induces structural changes in the inter- and intrasubunit transmembrane domain (TMD) cavities

General anesthetics exert many of their CNS actions by binding to and modulating membrane-embedded pentameric ligand-gated ion channels (pLGICs). The structural mechanisms underlying how anesthetics modulate pLGIC function remain largely unknown. GLIC, a prokaryotic pLGIC homologue, is inhibited by general anesthetics, suggesting anesthetics stabilize a closed channel state, but in anesthetic-bound GLIC crystal structures the channel appears open. Here, using functional GLIC channels expressed in oocytes, we examined whether propofol induces structural rearrangements in the GLIC transmembrane domain (TMD). Residues in the GLIC TMD that frame intrasubunit and intersubunit water-accessible cavities were individually mutated to cysteine. We measured and compared the rates of modification of the introduced cysteines by sulfhydryl-reactive reagents in the absence and presence of propofol. Propofol slowed the rate of modification of L240C (intersubunit) and increased the rate of modification of T254C (intrasubunit), indicating that propofol binding induces structural rearrangements in these cavities that alter the local environment near these residues. Propofol acceleration of T254C modification suggests that in the resting state propofol does not bind in the TMD intrasubunit cavity as observed in the crystal structure of GLIC with bound propofol (Nury, H., Van Renterghem, C., Weng, Y., Tran, A., Baaden, M., Dufresne, V., Changeux, J. P., Sonner, J. M., Delarue, M., and Corringer, P. J. (2011) Nature 469, 428-431). In silico docking using a GLIC closed channel homology model suggests propofol binds to intersubunit sites in the TMD in the resting state. Propofol-induced motions in the intersubunit cavity were distinct from motions associated with channel activation, indicating propofol stabilizes a novel closed state.

Exploring volatile general anesthetic binding to a closed membrane-bound bacterial voltage-gated sodium channel via computation

Despite the clinical ubiquity of anesthesia, the molecular basis of anesthetic action is poorly understood. Amongst the many molecular targets proposed to contribute to anesthetic effects, the voltage gated sodium channels (VGSCs) should also be considered relevant, as they have been shown to be sensitive to all general anesthetics tested thus far. However, binding sites for VGSCs have not been identified. Moreover, the mechanism of inhibition is still largely unknown. The recently reported atomic structures of several members of the bacterial VGSC family offer the opportunity to shed light on the mechanism of action of anesthetics on these important ion channels. To this end, we have performed a molecular dynamics "flooding" simulation on a membrane-bound structural model of the archetypal bacterial VGSC, NaChBac in a closed pore conformation. This computation allowed us to identify binding sites and access pathways for the commonly used volatile general anesthetic, isoflurane. Three sites have been characterized with binding affinities in a physiologically relevant range. Interestingly, one of the most favorable sites is in the pore of the channel, suggesting that the binding sites of local and general anesthetics may overlap. Surprisingly, even though the activation gate of the channel is closed, and therefore the pore and the aqueous compartment at the intracellular side are disconnected, we observe binding of isoflurane in the central cavity. Several sampled association and dissociation events in the central cavity provide consistent support to the hypothesis that the "fenestrations" present in the membrane-embedded region of the channel act as the long-hypothesized hydrophobic drug access pathway.

Signal transduction pathways in the pentameric ligand-gated ion channels

The mechanisms of allosteric action within pentameric ligand-gated ion channels (pLGICs) remain to be determined. Using crystallography, site-directed mutagenesis, and two-electrode voltage clamp measurements, we identified two functionally relevant sites in the extracellular (EC) domain of the bacterial pLGIC from Gloeobacter violaceus (GLIC). One site is at the C-loop region, where the NQN mutation (D91N, E177Q, and D178N) eliminated inter-subunit salt bridges in the open-channel GLIC structure and thereby shifted the channel activation to a higher agonist concentration. The other site is below the C-loop, where binding of the anesthetic ketamine inhibited GLIC currents in a concentration dependent manner. To understand how a perturbation signal in the EC domain, either resulting from the NQN mutation or ketamine binding, is transduced to the channel gate, we have used the Perturbation-based Markovian Transmission (PMT) model to determine dynamic responses of the GLIC channel and signaling pathways upon initial perturbations in the EC domain of GLIC. Despite the existence of many possible routes for the initial perturbation signal to reach the channel gate, the PMT model in combination with Yen's algorithm revealed that perturbation signals with the highest probability flow travel either via the β1–β2 loop or through pre-TM1. The β1–β2 loop occurs in either intra- or inter-subunit pathways, while pre-TM1 occurs exclusively in inter-subunit pathways. Residues involved in both types of pathways are well supported by previous experimental data on nAChR. The direct coupling between pre-TM1 and TM2 of the adjacent subunit adds new insight into the allosteric signaling mechanism in pLGICs.

Controlling potassium selectivity and proton blocking in a hybrid biological/solid-state polymer nanoporous membrane

Specific separations of protons and cations are usually performed by electromembrane processes, which require external electric energy. An easier process would be using a membrane able to separate both entities by passive diffusion. Presently, such synthetic nanoporous membranes do not exist. Here, we report the production of a robust hybrid biological/artificial solid-state membrane, which allows selective permeation of alkali metal cations without competing or concurrent permeation of protons. This membrane is simple to prepare and is based on the hydrophobic nature of the polymeric pore walls, and the confined gramicidin A molecules within. This work opens a new route for separation in the domain of nanobiofiltration, especially for tunable nanodevices based on differential ion conduction, with a fundamental understanding of the confinement mechanism.

Experimental determination of the vertical alignment between the second and third transmembrane segments of muscle nicotinic acetylcholine receptors

Nicotinic acetylcholine receptors (nAChR) are members of the Cys-loop ligand-gated ion channel superfamily. Muscle nAChR are heteropentamers that assemble from two α, and one each of β, γ, and δ subunits. Each subunit is composed of three domains, extracellular, transmembrane and intracellular. The transmembrane domain consists of four α-helical segments (M1-M4). Pioneering structural information was obtained using electronmicroscopy of Torpedo nAChR. The recently solved X-ray structure of the first eukaryotic Cys-loop receptor, a truncated (intracellular domain missing) glutamate-gated chloride channel α (GluClα) showed the same overall architecture. However, a significant difference with regard to the vertical alignment between the channel-lining segment M2 and segment M3 was observed. Here, we used functional studies utilizing disulfide trapping experiments in muscle nAChR to determine the spatial orientation between M2 and M3. Our results are in agreement with the vertical alignment as obtained when using the GluClα structure as a template to homology model muscle nAChR, however, they cannot be reconciled with the current Torpedo nAChR model. The vertical M2-M3 alignments as observed in X-ray structures of prokaryotic Gloeobacter violaceus ligand-gated ion channel and GluClα are in agreement. Our results further confirm that this alignment in Cys-loop receptors is conserved between prokaryotes and eukaryotes.

An outline of desensitization in pentameric ligand-gated ion channel receptors

Pentameric ligand-gated ion channel (pLGIC) receptors exhibit desensitization, the progressive reduction in ionic flux in the prolonged presence of agonist. Despite its pathophysiological importance and the fact that it was first described over half a century ago, surprisingly little is known about the structural basis of desensitization in this receptor family. Here, we explain how desensitization is defined using functional criteria. We then review recent progress into reconciling the structural and functional basis of this phenomenon. The extracellular-transmembrane domain interface is a key locus. Activation is well known to involve conformational changes at this interface, and several lines of evidence suggest that desensitization involves a distinct conformational change here that is incompatible with activation. However, major questions remain unresolved, including the structural basis of the desensitization-induced agonist affinity increase and the mechanism of pore closure during desensitization.

Structural sensitivity of a prokaryotic pentameric ligand-gated ion channel to its membrane environment

Although the activity of the nicotinic acetylcholine receptor (nAChR) is exquisitely sensitive to its membrane environment, the underlying mechanisms remain poorly defined. The homologous prokaryotic pentameric ligand-gated ion channel, Gloebacter ligand-gated ion channel (GLIC), represents an excellent model for probing the molecular basis of nAChR sensitivity because of its high structural homology, relative ease of expression, and amenability to crystallographic analysis. We show here that membrane-reconstituted GLIC exhibits structural and biophysical properties similar to those of the membrane-reconstituted nAChR, although GLIC is substantially more thermally stable. GLIC, however, does not possess the same exquisite lipid sensitivity. In particular, GLIC does not exhibit the same propensity to adopt an uncoupled conformation where agonist binding is uncoupled from channel gating. Structural comparisons provide insight into the chemical features that may predispose the nAChR to the formation of an uncoupled state.

Influences of membrane mimetic environments on membrane protein structures

The number of membrane protein structures in the Protein Data Bank is becoming significant and growing. Here, the transmembrane domain structures of the helical membrane proteins are evaluated to assess the influences of the membrane mimetic environments. Toward this goal, many of the biophysical properties of membranes are discussed and contrasted with those of the membrane mimetics commonly used for structure determination. Although the mimetic environments can perturb the protein structures to an extent that potentially gives rise to misinterpretation of functional mechanisms, there are also many structures that have a native-like appearance. From this assessment, an initial set of guidelines is proposed for distinguishing native-like from nonnative-like membrane protein structures. With experimental techniques for validation and computational methods for refinement and quality assessment and enhancement, there are good prospects for achieving native-like structures for these very important proteins.

Structural basis for ion permeation mechanism in pentameric ligand-gated ion channels

To understand the molecular mechanism of ion permeation in pentameric ligand-gated ion channels (pLGIC), we solved the structure of an open form of GLIC, a prokaryotic pLGIC, at 2.4 Å. Anomalous diffraction data were used to place bound anions and cations. This reveals ordered water molecules at the level of two rings of hydroxylated residues (named Ser6' and Thr2') that contribute to the ion selectivity filter. Two water pentagons are observed, a self-stabilized ice-like water pentagon and a second wider water pentagon, with one sodium ion between them. Single-channel electrophysiology shows that the side-chain hydroxyl of Ser6' is crucial for ion translocation. Simulations and electrostatics calculations complemented the description of hydration in the pore and suggest that the water pentagons observed in the crystal are important for the ion to cross hydrophobic constriction barriers. Simulations that pull a cation through the pore reveal that residue Ser6' actively contributes to ion translocation by reorienting its side chain when the ion is going through the pore. Generalization of these findings to the pLGIC family is proposed.

The structural mechanism of the Cys-loop receptor desensitization

The cys-loop receptors are neurotransmitter-operated ion channels, which mediate fast synaptic transmission for communication between neurons. However, prolonged exposure to the neurotransmitter drives the receptor to a desensitization state, which plays an important role in shaping synaptic transmission. Much progress has been made through more than half a century's research since Katz and Thesleff first descried desensitization for muscle nicotinic acetylcholine receptor. In this review, we summarized recent research developments of receptor desensitization. Now, it has been identified that many parts of the receptor, such as the pore domain (including the hinge in the M2-M3 linker), the binding domain, the coupling region, and the intracellular domain, are all involved in the cys-loop receptor desensitization and that uncoupling between the amino-terminal domain and channel lining domain seems to play a central role in desensitization. This uncoupling is mainly governed by the balance between coupling strength and relative tightness of gating machinery and influenced by other parts of the receptor. Agonist binding induces conformational change to overcome the gating barrier to open the channel through the stressed coupling region, which is subsequently broken, causing receptor desensitization. With rapid advancement in structural biology of membrane receptors, final validation of this mechanism is expected to occur in the near future when the high-resolution structure of the desensitized state is available.

Asymmetric ligand binding facilitates conformational transitions in pentameric ligand-gated ion channels

The anesthetic propofol inhibits the currents of the homopentameric ligand-gated ion channel GLIC, yet the crystal structure of GLIC with five propofol molecules bound symmetrically shows an open-channel conformation. To address this dilemma and determine if the symmetry of propofol binding sites affects the channel conformational transition, we performed a total of 1.5 μs of molecular dynamics simulations for different GLIC systems with propofol occupancies of 0, 1, 2, 3, and 5. GLIC without propofol binding or with five propofol molecules bound symmetrically, showed similar channel conformation and hydration status over multiple replicates of 100-ns simulations. In contrast, asymmetric binding to one, two or three equivalent sites in different subunits accelerated the channel dehydration, increased the conformational heterogeneity of the pore-lining TM2 helices, and shifted the lateral and radial tilting angles of TM2 toward a closed-channel conformation. The results differentiate two groups of systems based on the propofol binding symmetry. The difference between symmetric and asymmetric groups is correlated with the variance in the propofol-binding cavity adjacent to the hydrophobic gate and the force imposed by the bound propofol. Asymmetrically bound propofol produced greater variance in the cavity size that could further elevate the conformation heterogeneity. The force trajectory generated by propofol in each subunit over the course of a simulation exhibits an ellipsoidal shape, which has the larger component tangential to the pore. Asymmetric propofol binding creates an unbalanced force that expedites the channel conformation transitions. The findings from this study not only suggest that asymmetric binding underlies the propofol functional inhibition of GLIC, but also advocate for the role of symmetry breaking in facilitating channel conformational transitions.

The energy and work of a ligand-gated ion channel

Ligand-gated ion channels are allosteric membrane proteins that isomerize between C(losed) and O(pen) conformations. A difference in affinity for ligands in the two states influences the C↔O "gating" equilibrium constant. The energies associated with adult-type mouse neuromuscular nicotinic acetylcholine receptor (AChR) channel gating have been measured by using single-channel electrophysiology. Without ligands, the free energy, enthalpy and entropy of gating are ΔG0=+8.4, ΔH0=+10.9 and TΔS0=+2.5kcal/mol (-100mV, 23°C). Many mutations throughout the protein change ΔG0, including natural ones that cause disease. Agonists and most mutations change approximately independently the ground-state energy difference; thus, it is possible to forecast and engineer AChR responses simply by combining perturbations. The free energy of the low↔high affinity change for the neurotransmitter at each of two functionally equivalent binding sites is ΔGB(ACh)=-5.1kcal/mol. ΔGB(ACh) is set mainly by interactions of ACh with just three binding site aromatic groups. For a series of structurally related agonists, there is a correlation between the energies of low- and high-affinity binding, which implies that gating commences with the formation of the low-affinity complex. Brief, intermediate states in binding and gating have been detected. Several proposals for the nature of the gating transition-state energy landscape and the isomerization mechanism are discussed.

Probing microscopic material properties inside simulated membranes through spatially resolved three-dimensional local pressure fields and surface tensions

Cellular lipid membranes are spatially inhomogeneous soft materials. Materials properties such as pressure and surface tension thus show important microscopic-scale variation that is critical to many biological functions. We present a means to calculate pressure and surface tension in a 3D-resolved manner within molecular-dynamics simulations and show how such measurements can yield important insight. We also present the first corrections to local virial and pressure fields to account for the constraints typically used in lipid simulations that otherwise cause problems in highly oriented systems such as bilayers. Based on simulations of an asymmetric bacterial ion channel in a POPC bilayer, we demonstrate how 3D-resolved pressure can probe for both short-range and long-range effects from the protein on the membrane environment. We also show how surface tension is a sensitive metric for inter-leaflet equilibrium and can be used to detect even subtle imbalances between bilayer leaflets in a membrane-protein simulation. Since surface tension is known to modulate the function of many proteins, this effect is an important consideration for predictions of ion channel function. We outline a strategy by which our local pressure measurements, which we make available within a version of the GROMACS simulation package, may be used to design optimally equilibrated membrane-protein simulations.

Identification of propofol binding sites in a nicotinic acetylcholine receptor with a photoreactive propofol analog

Propofol, a widely used intravenous general anesthetic, acts at anesthetic concentrations as a positive allosteric modulator of γ-aminobutyric acid type A receptors and at higher concentration as an inhibitor of nicotinic acetylcholine receptors (nAChRs). Here, we characterize propofol binding sites in a muscle-type nAChR by use of a photoreactive analog of propofol, 2-isopropyl-5-[3-(trifluoromethyl)-3H-diazirin-3-yl]phenol (AziPm). Based upon radioligand binding assays, AziPm stabilized the Torpedo nAChR in the resting state, whereas propofol stabilized the desensitized state. nAChR-rich membranes were photolabeled with [(3)H]AziPm, and labeled amino acids were identified by Edman degradation. [(3)H]AziPm binds at three sites within the nAChR transmembrane domain: (i) an intrasubunit site in the δ subunit helix bundle, photolabeling in the nAChR desensitized state (+agonist) δM2-18' and two residues in δM1 (δPhe-232 and δCys-236); (ii) in the ion channel, photolabeling in the nAChR resting, closed channel state (-agonist) amino acids in the M2 helices (αM2-6', βM2-6' and -13', and δM2-13') that line the channel lumen (with photolabeling reduced by >90% in the desensitized state); and (iii) at the γ-α interface, photolabeling αM2-10'. Propofol enhanced [(3)H]AziPm photolabeling at αM2-10'. Propofol inhibited [(3)H]AziPm photolabeling within the δ subunit helix bundle at lower concentrations (IC50 = 40 μm) than it inhibited ion channel photolabeling (IC50 = 125 μm). These results identify for the first time a single intrasubunit propofol binding site in the nAChR transmembrane domain and suggest that this is the functionally relevant inhibitory binding site.