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

Ligand-Dependent Intramolecular Motion of Native Nicotinic Acetylcholine Receptors Determined in Living Myotube Cells via Diffracted X-ray Tracking

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that play an important role in signal transduction at the neuromuscular junction (NMJ). Movement of the nAChR extracellular domain following agonist binding induces conformational changes in the extracellular domain, which in turn affects the transmembrane domain and opens the ion channel. It is known that the surrounding environment, such as the presence of specific lipids and proteins, affects nAChR function. Diffracted X-ray tracking (DXT) facilitates measurement of the intermolecular motions of receptors on the cell membranes of living cells, including all the components involved in receptor function. In this study, the intramolecular motion of the extracellular domain of native nAChR proteins in living myotube cells was analyzed using DXT for the first time. We revealed that the motion of the extracellular domain in the presence of an agonist (e.g., carbamylcholine, CCh) was restricted by an antagonist (i.e., alpha-bungarotoxin, BGT).

Structural Elucidation of Ivermectin Binding to α7nAChR and the Induced Channel Desensitization

The α7 nicotinic acetylcholine receptor (α7nAChR) mediates signaling in the central nervous system and cholinergic anti-inflammatory pathways. Ivermectin is a positive allosteric modulator of a full-length α7nAChR and an agonist of the α7nAChR construct containing transmembrane (TMD) and intracellular (ICD) domains, but structural insights of the binding have not previously been determined. Here, combining nuclear magnetic resonance as a primary experimental tool with Rosetta comparative modeling and molecular dynamics simulations, we have revealed details of ivermectin binding to the α7nAChR TMD + ICD and corresponding structural changes in an ivermectin-induced desensitized state. Ivermectin binding was stabilized predominantly by hydrophobic interactions from interfacial residues between adjacent subunits near the extracellular end of the TMD, where the inter-subunit gap was substantially expanded in comparison to the apo structure. The ion-permeation pathway showed a profile distinctly different from the resting-state profile but similar to profiles of desensitized α7nAChR. The ICD also exhibited structural changes, including reorientation of the MX and h3 helices relative to the channel axis. The resulting structures of the α7nAChR TMD + ICD in complex with ivermectin provide opportunities for discovering new modulators of therapeutic potential and exploring the structural basis of cytoplasmic signaling under different α7nAChR functional states.

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

Conformational decoupling in acid-sensing ion channels uncovers mechanism and stoichiometry of PcTx1-mediated inhibition

Acid-sensing ion channels (ASICs) are trimeric proton-gated cation channels involved in fast synaptic transmission. Pharmacological inhibition of ASIC1a reduces neurotoxicity and stroke infarct volumes, with the cysteine knot toxin psalmotoxin-1 (PcTx1) being one of the most potent and selective inhibitors. PcTx1 binds at the subunit interface in the extracellular domain (ECD), but the mechanism and conformational consequences of the interaction, as well as the number of toxin molecules required for inhibition, remain unknown. Here, we use voltage-clamp fluorometry and subunit concatenation to decipher the mechanism and stoichiometry of PcTx1 inhibition of ASIC1a. Besides the known inhibitory binding mode, we propose PcTx1 to have at least two additional binding modes that are decoupled from the pore. One of these modes induces a long-lived ECD conformation that reduces the activity of an endogenous neuropeptide. This long-lived conformational state is proton-dependent and can be destabilized by a mutation that decreases PcTx1 sensitivity. Lastly, the use of concatemeric channel constructs reveals that disruption of a single PcTx1 binding site is sufficient to destabilize the toxin-induced conformation, while functional inhibition is not impaired until two or more binding sites are mutated. Together, our work provides insight into the mechanism of PcTx1 inhibition of ASICs and uncovers a prolonged conformational change with possible pharmacological implications.

Markov state models of proton- and pore-dependent activation in a pentameric ligand-gated ion channel

Ligand-gated ion channels conduct currents in response to chemical stimuli, mediating electrochemical signaling in neurons and other excitable cells. For many channels, the details of gating remain unclear, partly due to limited structural data and simulation timescales. Here, we used enhanced sampling to simulate the pH-gated channel GLIC, and construct Markov state models (MSMs) of gating. Consistent with new functional recordings, we report in oocytes, our analysis revealed differential effects of protonation and mutation on free-energy wells. Clustering of closed- versus open-like states enabled estimation of open probabilities and transition rates, while higher-order clustering affirmed conformational trends in gating. Furthermore, our models uncovered state- and protonation-dependent symmetrization. This demonstrates the applicability of MSMs to map energetic and conformational transitions between ion-channel functional states, and how they reproduce shifts upon activation or mutation, with implications for modeling neuronal function and developing state-selective drugs.

Asymmetric opening of the homopentameric 5-HT3A serotonin receptor in lipid bilayers

Pentameric ligand-gated ion channels (pLGICs) of the Cys-loop receptor family are key players in fast signal transduction throughout the nervous system. They have been shown to be modulated by the lipid environment, however the underlying mechanism is not well understood. We report three structures of the Cys-loop 5-HT_3A serotonin receptor (5HT₃R) reconstituted into saposin-based lipid bilayer discs: a symmetric and an asymmetric apo state, and an asymmetric agonist-bound state. In comparison to previously published 5HT₃R conformations in detergent, the lipid bilayer stabilises the receptor in a more tightly packed, ‘coupled’ state, involving a cluster of highly conserved residues. In consequence, the agonist-bound receptor conformation adopts a wide-open pore capable of conducting sodium ions in unbiased molecular dynamics (MD) simulations. Taken together, we provide a structural basis for the modulation of 5HT₃R by the membrane environment, and a model for asymmetric activation of the receptor.

Pentameric ligand-gated ion channels (pLGICs) are key players in neurotransmission and have been shown to be modulated by the lipid environment, however the underlying mechanism is not well understood. Here, the authors report structures of the pLGIC 5-HT_3A serotonin receptor reconstituted into lipid bilayer discs and reveal lipid–protein interactions as well as asymmetric activation of the homopentameric receptor.

In silico studies of the interactions between propofol and fentanyl using Gaussian accelerated molecular dynamics

Fentanyl is a potent opioid analgesic, which for decades has been used routinely in surgical and therapeutic applications. In addition to its analgesic properties, fentanyl also possesses anesthetic properties, which are not well understood. Fentanyl is used in the general anesthesia process to induce and maintain anesthesia in combination with the general anesthetic propofol, which fentanyl is known to potentiate. As the atomic-level mechanism behind the potentiation of propofol is unclear, we have used classical molecular dynamics simulations to study the interactions of these drugs with the Gloeobacter violaceus ion channel (GLIC). This ion channel has been identified as a target for many anesthetic drugs. We identified multiple binding sites using flooding style and Gaussian accelerated molecular dynamics (GaMD) simulations, showing fentanyl acting as a stabiliser that holds propofol within binding sites. Our extensive GaMD simulations were also able to show the pathway by which propofol blocks the channel pore, which has previously been suggested as a mechanism for ion channel modulation. General anesthesia is a multi-drug process and this study provides the first insight into the interactions between two different drugs in the anesthesia process in a relevant biological environment.Communicated by Ramaswamy H. Sarma.

Nicotinic receptor pharmacology in silico: Insights and challenges

Nicotinic acetylcholine receptors (nAChR) are homo- or hetero-pentameric ligand-gated ion channels of the Cys-loop superfamily and play important roles in the nervous system and muscles. Studies on nAChR benefit from in silico modeling due to the lack of high-resolution structures for most receptor subtypes and challenges in experiments addressing the complex mechanism of activation involving allosteric sites. Although there is myriad of computational modeling studies on nAChR, the multitude of the methods and parameters used in these studies makes modeling nAChR a daunting task, particularly for the non-experts in the field. To address this problem, the modeling literature on Torpedo nAChR and α7 nAChR were focused on as examples of heteromeric and homomeric nAChR, and the key in silico modeling studies between the years 1995-2019 were concisely reviewed. This was followed by a critical analysis of these studies by comparing the findings with each other and with the emerging experimental and computational data on nAChR. Based on these critical analyses, suggestions were made to guide the future researchers in the field of in silico modeling of nAChR. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

Modulation of the Gloeobacter violaceus Ion Channel by Fentanyl: A Molecular Dynamics Study

Fentanyl is an opioid analgesic, which is routinely used in general surgery to suppress the sensation of pain and as the analgesic component in the induction and maintenance of anesthesia. Fentanyl is also used as the main component to induce anesthesia and as a potentiator to the general anesthetic propofol. The mechanism by which fentanyl induces its anesthetic action is still unclear, and we have therefore employed fully atomistic molecular dynamics simulations to probe this process by simulating the interactions of fentanyl with the Gloeobacter violaceus ligand-gated ion channel (GLIC). In this paper, we identify multiple extracellular fentanyl binding sites, which are different from the transmembrane general anesthetic binding sites observed for propofol and other general anesthetics. Our simulations identify a novel fentanyl binding site within the GLIC that results in conformational changes that inhibit conduction through the channel.

Closed-Locked and Apo-Resting State Structures of the Human α7 Nicotinic Receptor: A Computational Study

Nicotinic acetylcholine receptors, belonging to the Cys-loop superfamily of ligand-gated ion channels (LGICs), are membrane proteins present in neurons and at neuromuscular junctions. They are responsible for signal transmission, and their function is regulated by neurotransmitters, agonists, and antagonists drugs. A detailed knowledge of their conformational transition in response to ligand binding is critical to understanding the basis of ligand-receptor interaction, in view of new pharmacological approaches to control receptor activity. However, the scarcity of experimentally derived structures of human channels makes this perspective extremely challenging. To contribute overcoming this issue, we have recently reported structural models for the open and the desensitized states of the human α7 nicotinic receptor. Here, we provide all-atom structural models of the same receptor in two different nonconductive states. The first structure, built via homology modeling and relaxed with extensive Molecular Dynamics simulations, represents the receptor bound to the natural antagonist α-conotoxin ImI. After comparison with available experimental data and computational models of other eukaryotic LGICs, we deem it consistent with the "closed-locked" state. The second model, obtained with simulations from the spontaneous relaxation of the open, agonist-bound α7 structure after ligand removal, recapitulates the characteristics of the apo-resting state of the receptor. These results add to our previous work on the active and desensitized state conformations, contributing to the structural characterization of the conformational landscape of the human α7 receptor and suggesting benchmarks to discriminate among conformations found in experiments or in simulations of LGICs. In particular key interactions at the interface between the extracellular domain and the transmembrane domain are identified, that could be critical to the α7 receptor function.

Allosteric potentiation of a ligand-gated ion channel is mediated by access to a deep membrane-facing cavity

Theories of general anesthesia have shifted in focus from bulk lipid effects to specific interactions with membrane proteins. Target receptors include several subtypes of pentameric ligand-gated ion channels; however, structures of physiologically relevant proteins in this family have yet to define anesthetic binding at high resolution. Recent cocrystal structures of the bacterial protein GLIC provide snapshots of state-dependent binding sites for the common surgical agent propofol (PFL), offering a detailed model system for anesthetic modulation. Here, we combine molecular dynamics and oocyte electrophysiology to reveal differential motion and modulation upon modification of a transmembrane binding site within each GLIC subunit. WT channels exhibited net inhibition by PFL, and a contraction of the cavity away from the pore-lining M2 helix in the absence of drug. Conversely, in GLIC variants exhibiting net PFL potentiation, the cavity was persistently expanded and proximal to M2. Mutations designed to favor this deepened site enabled sensitivity even to subclinical concentrations of PFL, and a uniquely prolonged mode of potentiation evident up to ∼30 min after washout. Dependence of these prolonged effects on exposure time implicated the membrane as a reservoir for a lipid-accessible binding site. However, at the highest measured concentrations, potentiation appeared to be masked by an acute inhibitory effect, consistent with the presence of a discrete, water-accessible site of inhibition. These results support a multisite model of transmembrane allosteric modulation, including a possible link between lipid- and receptor-based theories that could inform the development of new anesthetics.

Identification of binding sites contributing to volatile anesthetic effects on GABA type A receptors

Most general anesthetics enhance GABA type A (GABA_A) receptor activity at clinically relevant concentrations. Sites of action of volatile anesthetics on the GABA_A receptor remain unknown, whereas sites of action of many intravenous anesthetics have been identified in GABA_A receptors by using photolabeling. Here, we used photoactivatable analogs of isoflurane (AziISO) and sevoflurane (AziSEVO) to locate their sites on α₁β₃γ_2L and α₁β₃ GABA_A receptors. As with isoflurane and sevoflurane, AziISO and AziSEVO enhanced the currents elicited by GABA. AziISO and AziSEVO each labeled 10 residues in α₁β₃ receptors and 9 and 8 residues, respectively, in α₁β₃γ_2L receptors. Photolabeled residues were concentrated in transmembrane domains and located in either subunit interfaces or in the interface between the extracellular domain and the transmembrane domain. The majority of these transmembrane residues were protected from photolabeling with the addition of excess parent anesthetic, which indicated specificity. Binding sites were primarily located within α+/β- and β+/α- subunit interfaces, but residues in the α+/γ- interface were also identified, which provided a basis for differential receptor subtype sensitivity. Isoflurane and sevoflurane did not always share binding sites, which suggests an unexpected degree of selectivity.-Woll, K. A., Zhou, X., Bhanu, N. V., Garcia, B. A., Covarrubias, M., Miller, K. W., Eckenhoff, R. G. Identification of binding sites contributing to volatile anesthetic effects on GABA type A receptors.

Recent progress on the molecular pharmacology of propofol

The precise mechanism by which propofol enhances GABAergic transmission remains unclear, but much progress has been made regarding the underlying structural and dynamic mechanisms. Furthermore, it is now clear that propofol has additional molecular targets, many of which are functionally influenced at concentrations achieved clinically. Focusing primarily on molecular targets, this brief review attempts to summarize some of this recent progress while pointing out knowledge gaps and controversies. It is not intended to be comprehensive but rather to stimulate further thought, discussion, and study on the mechanisms by which propofol produces its pleiotropic effects.

Ketamine Inhibition of the Pentameric Ligand-Gated Ion Channel GLIC

Ketamine inhibits pentameric ligand-gated ion channels (pLGICs), including the bacterial pLGIC from Gloeobacter violaceus (GLIC). The crystal structure of GLIC shows R-ketamine bound to an extracellular intersubunit cavity. Here, we performed molecular dynamics simulations of GLIC in the absence and presence of R- or S-ketamine. No stable binding of S-ketamine in the original cavity was observed in the simulations, largely due to its unfavorable access to residue D154, which provides important electrostatic interactions to stabilize R-ketamine binding. Contrary to the symmetric binding shown in the crystal structure, R-ketamine moved away from some of the binding sites and was bound to GLIC asymmetrically at the end of simulations. The asymmetric binding is consistent with the experimentally measured negative cooperativity of ketamine binding to GLIC. In the presence of R-ketamine, all subunits showed changes in structure and dynamics, irrespective of binding stability; the extracellular intersubunit cavity expanded and intersubunit electrostatic interactions involved in channel activation were altered. R-ketamine binding promoted a conformational shift toward closed GLIC. Conformational changes near the ketamine-binding site were propagated to the interface between the extracellular and transmembrane domains, and further to the pore-lining TM2 through two pathways: pre-TM1 and the β1-β2 loop. Both signaling pathways have been predicted previously using the perturbation-based Markovian transmission model. The study provides a structural and dynamics basis for the inhibitory modulation of ketamine on pLGICs.

Signalling assemblies: the odds of symmetry

The assembly of proteins into complexes is fundamental to nearly all biological signalling processes. Symmetry is a dominant feature of the structures of experimentally determined protein complexes, observed in the vast majority of homomers and many heteromers. However, some asymmetric structures exist, and asymmetry also often forms transiently, intractable to traditional structure determination methods. Here, we explore the role of protein complex symmetry and asymmetry in cellular signalling, focusing on receptors, transcription factors and transmembrane channels, among other signalling assemblies. We highlight a recurrent tendency for asymmetry to be crucial for signalling function, often being associated with activated states. We conclude with a discussion of how consideration of protein complex symmetry and asymmetry has significant potential implications and applications for pharmacology and human disease.

A membrane-embedded pathway delivers general anesthetics to two interacting binding sites in the Gloeobacter violaceus ion channel

General anesthetics exert their effects on the central nervous system by acting on ion channels, most notably pentameric ligand-gated ion channels. Although numerous studies have focused on pentameric ligand-gated ion channels, the details of anesthetic binding and channel modulation are still debated. A better understanding of the anesthetic mechanism of action is necessary for the development of safer and more efficacious drugs. Herein, we present a computational study identifying two anesthetic binding sites in the transmembrane domain of the Gloeobacter violaceus ligand-gated ion channel (GLIC) channel, characterize the putative binding pathway, and observe structural changes associated with channel function. Molecular simulations of desflurane reveal a binding pathway to GLIC via a membrane-embedded tunnel using an intrasubunit protein lumen as the conduit, an observation that explains the Meyer-Overton hypothesis, or why the lipophilicity of an anesthetic and its potency are generally proportional. Moreover, employing high concentrations of ligand led to the identification of a second transmembrane site (TM2) that inhibits dissociation of anesthetic from the TM1 site and is consistent with the high concentrations of anesthetics required to achieve clinical effects. Finally, asymmetric binding patterns of anesthetic to the channel were found to promote an iris-like conformational change that constricts and dehydrates the ion pore, creating a 13.5 kcal/mol barrier to ion translocation. Together with previous studies, the simulations presented herein demonstrate a novel anesthetic binding site in GLIC that is accessed through a membrane-embedded tunnel and interacts with a previously known site, resulting in conformational changes that produce a non-conductive state of the channel.

The cellular membrane as a mediator for small molecule interaction with membrane proteins

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

A Novel Bifunctional Alkylphenol Anesthetic Allows Characterization of γ-Aminobutyric Acid, Type A (GABAA), Receptor Subunit Binding Selectivity in Synaptosomes

Propofol, an intravenous anesthetic, is a positive modulator of the GABAA receptor, but the mechanistic details, including the relevant binding sites and alternative targets, remain disputed. Here we undertook an in-depth study of alkylphenol-based anesthetic binding to synaptic membranes. We designed, synthesized, and characterized a chemically active alkylphenol anesthetic (2-((prop-2-yn-1-yloxy)methyl)-5-(3-(trifluoromethyl)-3H-diazirin-3-yl)phenol, AziPm-click (1)), for affinity-based protein profiling (ABPP) of propofol-binding proteins in their native state within mouse synaptosomes. The ABPP strategy captured ∼4% of the synaptosomal proteome, including the unbiased capture of five α or β GABAA receptor subunits. Lack of γ2 subunit capture was not due to low abundance. Consistent with this, independent molecular dynamics simulations with alchemical free energy perturbation calculations predicted selective propofol binding to interfacial sites, with higher affinities for α/β than γ-containing interfaces. The simulations indicated hydrogen bonding is a key component leading to propofol-selective binding within GABAA receptor subunit interfaces, with stable hydrogen bonds observed between propofol and α/β cavity residues but not γ cavity residues. We confirmed this by introducing a hydrogen bond-null propofol analogue as a protecting ligand for targeted-ABPP and observed a lack of GABAA receptor subunit protection. This investigation demonstrates striking interfacial GABAA receptor subunit selectivity in the native milieu, suggesting that asymmetric occupancy of heteropentameric ion channels by alkylphenol-based anesthetics is sufficient to induce modulation of activity.

Common Internal Allosteric Network Links Anesthetic Binding Sites in a Pentameric Ligand-Gated Ion Channel

General anesthetics bind reversibly to ion channels, modifying their global conformational distributions, but the underlying atomic mechanisms are not completely known. We examine this issue by way of the model protein Gloeobacter violaceous ligand-gated ion channel (GLIC) using computational molecular dynamics, with a coarse-grained model to enhance sampling. We find that in flooding simulations, both propofol and a generic particle localize to the crystallographic transmembrane anesthetic binding region, and that propofol also localizes to an extracellular region shared with the crystallographic ketamine binding site. Subsequent simulations to probe these binding modes in greater detail demonstrate that ligand binding induces structural asymmetry in GLIC. Consequently, we employ residue interaction correlation analysis to describe the internal allosteric network underlying the coupling of ligand and distant effector sites necessary for conformational change. Overall, the results suggest that the same allosteric network may underlie the actions of various anesthetics, regardless of binding site.

Common Anesthetic-binding Site for Inhibition of Pentameric Ligand-gated Ion Channels

BACKGROUND

Identifying functionally relevant anesthetic-binding sites in pentameric ligand-gated ion channels (pLGICs) is an important step toward understanding the molecular mechanisms underlying anesthetic action. The anesthetic propofol is known to inhibit cation-conducting pLGICs, including a prokaryotic pLGIC from Erwinia chrysanthemi (ELIC), but the sites responsible for functional inhibition remain undetermined.

METHODS

We photolabeled ELIC with a light-activated derivative of propofol (AziPm) and performed fluorine-19 nuclear magnetic resonance experiments to support propofol binding to a transmembrane domain (TMD) intrasubunit pocket. To differentiate sites responsible for propofol inhibition from those that are functionally irrelevant, we made an ELIC-γ-aminobutyric acid receptor (GABAAR) chimera that replaced the ELIC-TMD with the α1β3GABAAR-TMD and compared functional responses of ELIC-GABAAR and ELIC with propofol modulations.

RESULTS

Photolabeling showed multiple AziPm-binding sites in the extracellular domain (ECD) but only one site in the TMD with labeled residues M265 and F308 in the resting state of ELIC. Notably, this TMD site is an intrasubunit pocket that overlaps with binding sites for anesthetics, including propofol, found previously in other pLGICs. Fluorine-19 nuclear magnetic resonance experiments supported propofol binding to this TMD intrasubunit pocket only in the absence of agonist. Functional measurements of ELIC-GABAAR showed propofol potentiation of the agonist-elicited current instead of inhibition observed on ELIC.

CONCLUSIONS

The distinctly different responses of ELIC and ELIC-GABAAR to propofol support the functional relevance of propofol binding to the TMD. Combining the newly identified TMD intrasubunit pocket in ELIC with equivalent TMD anesthetic sites found previously in other cationic pLGICs, we propose this TMD pocket as a common site for anesthetic inhibition of pLGICs.

Asymmetric dynamics of ion channel forming proteins - Hepatitis C virus (HCV) p7 bundles

Protein p7 of hepatitis C virus (HCV) is a short 63 amino acid membrane protein which homo-oligomerises in the lipid membrane to form ion and proton conducting bundles. Two different genotypes (GTs) of p7, 1a and 5a, are used to simulate hexameric bundles of the protein embedded in a fully hydrated lipid bilayer during 400 ns molecular dynamics (MD) simulations. Whilst the bundle of GT 1a is based on a fully computational derived structure, the bundle of GT 5a is based on NMR spectroscopic data. Results of a full correlation analysis (FCA) reveal that albeit structural differences both bundles screen local minima during the simulation. The collective motion of the protein domains is asymmetric. No 'breathing-mode'-like dynamics is observed. The presence of divalent ions, such as Ca-ions affects the dynamics of especially solvent exposed parts of the protein, but leaves the asymmetric domain motion unaffected.

Pore dilatation increases the bicarbonate permeability of CFTR, ANO1 and glycine receptor anion channels

KEY POINTS

Cellular stimuli can modulate the ion selectivity of some anion channels, such as CFTR, ANO1 and the glycine receptor (GlyR), by changing pore size. Ion selectivity of CFTR, ANO1 and GlyR is critically affected by the electric permittivity and diameter of the channel pore. Pore size change affects the energy barriers of ion dehydration as well as that of size-exclusion of anion permeation. Pore dilatation increases the bicarbonate permeability (P HC O3/ Cl ) of CFTR, ANO1 and GlyR. Dynamic change in P HC O3/ Cl may mediate many physiological and pathological processes.

ABSTRACT

Chloride (Cl(-) ) and bicarbonate (HCO3 (-) ) are two major anions and their permeation through anion channels plays essential roles in our body. However, the mechanism of ion selection by the anion channels is largely unknown. Here, we provide evidence that pore dilatation increases the bicarbonate permeability (P HC O3/ Cl ) of anion channels by reducing energy barriers of size-exclusion and ion dehydration of HCO3 (-) permeation. Molecular, physiological and computational analyses of major anion channels, such as cystic fibrosis transmembrane conductance regulator (CFTR), anoctamin-1(ANO1/TMEM16A) and the glycine receptor (GlyR), revealed that the ion selectivity of anion channels is basically determined by the electric permittivity and diameter of the pore. Importantly, cellular stimuli dynamically modulate the anion selectivity of CFTR and ANO1 by changing the pore size. In addition, pore dilatation by a mutation in the pore-lining region alters the anion selectivity of GlyR. Changes in pore size affected not only the energy barriers of size exclusion but that of ion dehydration by altering the electric permittivity of water-filled cavity in the pore. The dynamic increase in P HC O3/ Cl by pore dilatation may have many physiological and pathophysiological implications ranging from epithelial HCO3 (-) secretion to neuronal excitation.

Direct Pore Binding as a Mechanism for Isoflurane Inhibition of the Pentameric Ligand-gated Ion Channel ELIC

Pentameric ligand-gated ion channels (pLGICs) are targets of general anesthetics, but molecular mechanisms underlying anesthetic action remain debatable. We found that ELIC, a pLGIC from Erwinia chrysanthemi, can be functionally inhibited by isoflurane and other anesthetics. Structures of ELIC co-crystallized with isoflurane in the absence or presence of an agonist revealed double isoflurane occupancies inside the pore near T237(6′) and A244(13′). A pore-radius contraction near the extracellular entrance was observed upon isoflurane binding. Electrophysiology measurements with a single-point mutation at position 6′ or 13′ support the notion that binding at these sites renders isoflurane inhibition. Molecular dynamics simulations suggested that isoflurane binding was more stable in the resting than in a desensitized pore conformation. This study presents compelling evidence for a direct pore-binding mechanism of isoflurane inhibition, which has a general implication for inhibitory action of general anesthetics on pLGICs.

Conformational gating dynamics in the GluCl anion-selective chloride channel

Cys-loop receptors are central to propagation of signals in the nervous system. The gating of the membrane-spanning pore is triggered by structural rearrangements in the agonist-binding site, located some 50 Å away from the pore. A sequential conformational change, propagating from the ligand-binding site to the pore, has been proposed to govern gating in all Cys-loop receptors. Here, we identify structural and dynamic components of the conformational gating in the eukaryotic glutamate-gated chloride channel (GluCl) by means of molecular dynamics (MD) simulations with and without the l-glutamate agonist bound. A significant increase in pore opening and accompanying hydration is observed in the presence of glutamate. Potential of mean force calculations reveal that the barrier for ion passage drops from 15 kcal/mol to 5-10 kcal/mol with the agonist bound. This appears to be explained by agonist binding that leads to significant changes in the intersubunit hydrogen-bonding pattern, which induce a slight tilt of the extracellular domain relative to the transmembrane domain in the simulations. This rearrangement is subtle, but correspond to the direction of the quaternary twist observed as a key difference between open and closed X-ray structures. While the full reversible gating is still a much slower process, the observed structural dynamics sheds new light on the early stages of how the agonist influences the extracellular domain, how the extracellular domain interacts with the transmembrane domain, and how changes in the transmembrane domain alter the free energy of ion passage.

A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking. Homology modeling and Molecular Dynamics simulations are valuable tools to predict structures of unknown proteins, however, for the case of human nAchRs, they have been unsuccessful in providing a stable open structure so far. This is due to different problems with the template structures: on one side the homology with prokaryotic species is too low, while on the other the open eukaryotic GluCl proved itself unstable in several MD studies and collapsed to a dehydrated, non-conductive conformation, even when bound to an agonist. Aim of this work is to obtain, by a mixing of state-of-the-art homology and simulation techniques, a plausible prediction of the structure (still unknown) of the open state of human α7 nAChR complexed with epibatidine, from which it is possible to start structural and functional test studies. To prevent channel closure we employ a restraint that keeps the transmembrane pore open, and obtain in this way a stable, hydrated conformation. To further validate this conformation, we run four long, unbiased simulations starting from configurations chosen at random along the restrained trajectory. The channel remains stable and hydrated over the whole runs. This allows to assess the stability of the putative open conformation over a cumulative time of 1 μs, 800 ns of which are of unbiased simulation. Mostly based on the analysis of pore hydration and size, we suggest that the obtained structure has reasonable chances to be (at least one of the possible) structures of the channel in the open conformation.

Conformational Changes Underlying Desensitization of the Pentameric Ligand-Gated Ion Channel ELIC

Structural rearrangements underlying functional transitions of pentameric ligand-gated ion channels (pLGICs) are not fully understood. Using (19)F nuclear magnetic resonance and electron spin resonance spectroscopy, we found that ELIC, a pLGIC from Erwinia chrysanthemi, expanded the extracellular end and contracted the intracellular end of its pore during transition from the resting to an apparent desensitized state. Importantly, the contraction at the intracellular end of the pore likely forms a gate to restrict ion transport in the desensitized state. This gate differs from the hydrophobic gate present in the resting state. Conformational changes of the TM2-TM3 loop were limited to the N-terminal end. The TM4 helices and the TM3-TM4 loop appeared relatively insensitive to agonist-mediated structural rearrangement. These results indicate that conformational changes accompanying functional transitions are not uniform among different ELIC regions. This work also revealed the co-existence of multiple conformations for a given state and suggested asymmetric conformational arrangements in a homomeric pLGIC.

Investigating ion channel conformational changes using voltage clamp fluorometry

Ion channels are membrane proteins whose functions are governed by conformational changes. The widespread distribution of ion channels, coupled with their involvement in most physiological and pathological processes and their importance as therapeutic targets, renders the elucidation of these conformational mechanisms highly compelling from a drug discovery perspective. Thanks to recent advances in structural biology techniques, we now have high-resolution static molecular structures for members of the major ion channel families. However, major questions remain to be resolved about the conformational states that ion channels adopt during activation, drug modulation and desensitization. Patch-clamp electrophysiology has long been used to define ion channel conformational states based on functional criteria. It achieves this by monitoring conformational changes at the channel gate and cannot detect conformational changes occurring in regions distant from the gate. Voltage clamp fluorometry involves labelling cysteines introduced into domains of interest with environmentally sensitive fluorophores and inferring structural rearrangements from voltage or ligand-induced fluorescence changes. Ion channel currents are monitored simultaneously to verify the conformational status. By defining real time conformational changes in domains distant from the gate, this technique provides unexpected new insights into ion channel structure and function. This review aims to summarise the methodology and highlight recent innovative applications of this powerful technique. This article is part of the Special Issue entitled 'Fluorescent Tools in Neuropharmacology'.

The evolution of pentameric ligand-gated ion-channels and the changing family of anthelmintic drug targets

SUMMARY Pentameric ligand-gated ion-channels rapidly transduce synaptic neurotransmitter signals to an electrical response in post-synaptic neuronal or muscle cells and control the neuromusculature of a majority of multicellular animals. A wide range of pharmaceuticals target these receptors including ethanol, nicotine, anti-depressants and other mood regulating drugs, compounds that control pain and mobility and are targeted by a majority of anthelmintic drugs used to control parasitic infection of humans and livestock. Major advances have been made in recent years to our understanding of the structure, function, activity and the profile of compounds that can activate specific receptors. It is becoming clear that these anthelmintic drug targets are not fixed, but differ in significant details from one nematode species to another. Here we review what is known about the evolution of the pentameric ligand-gated ion-channels, paying particular attention to the nematodes, how we can infer the origins of such receptors and understand the factors that determine how they change both over time and from one species to another. Using this knowledge provides a biological framework in which to understand these important drug targets and avenues to identify new receptors and aid the search for new anthelmintic drugs.

POVME 2.0: An Enhanced Tool for Determining Pocket Shape and Volume Characteristics

Analysis of macromolecular/small-molecule binding pockets can provide important insights into molecular recognition and receptor dynamics. Since its release in 2011, the POVME (POcket Volume MEasurer) algorithm has been widely adopted as a simple-to-use tool for measuring and characterizing pocket volumes and shapes. We here present POVME 2.0, which is an order of magnitude faster, has improved accuracy, includes a graphical user interface, and can produce volumetric density maps for improved pocket analysis. To demonstrate the utility of the algorithm, we use it to analyze the binding pocket of RNA editing ligase 1 from the unicellular parasite Trypanosoma brucei, the etiological agent of African sleeping sickness. The POVME analysis characterizes the full dynamics of a potentially druggable transient binding pocket and so may guide future antitrypanosomal drug-discovery efforts. We are hopeful that this new version will be a useful tool for the computational- and medicinal-chemist community.

Pentameric Ligand-gated Ion Channels : Insights from Computation

Pentameric ligand-gated ion channels (pLGICs) conduct upon the binding of an agonist and are fundamental to neurotransmission. New insights into the complex mechanisms underlying pLGIC gating, ion selectivity, and modulation have recently been gained via a series of crystal structures in prokaryotes and C .elegans, as well as computational studies relying on these structures. Here we review contributions from a variety of computational approaches, including normal mode analysis, automated docking, and fully atomistic molecular dynamics simulation. Examples from our own research, particularly concerning interactions with general anesthetics and lipids, are used to illustrate predictive results complementary to crystallographic studies.

Asymmetric perturbations of signalling oligomers

This review focuses on rapid and reversible noncovalent interactions for symmetric oligomers of signalling proteins. Symmetry mismatch, transient symmetry breaking and asymmetric perturbations via chemical (ligand binding) and physical (electric or mechanic) effects can initiate the signalling events. Advanced biophysical methods can reveal not only structural symmetries of stable membrane-bound signalling proteins but also asymmetric functional transition states. Relevant techniques amenable to distinguish between symmetric and asymmetric architectures are discussed including those with the capability of capturing low-populated transient conformational states. Typical examples of signalling proteins are overviewed for symmetry breaking in dimers (GPCRs, growth factor receptors, transcription factors); trimers (acid-sensing ion channels); tetramers (voltage-gated cation channels, ionotropic glutamate receptor, CNG and CHN channels); pentameric ligand-gated and mechanosensitive channels; higher order oligomers (gap junction channel, chaperonins, proteasome, virus capsid); as well as primary and secondary transporters. In conclusion, asymmetric perturbations seem to play important functional roles in a broad range of communicating networks.

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.

Photoaffinity labeling the propofol binding site in GLIC

Propofol, an intravenous general anesthetic, produces many of its anesthetic effects in vivo by potentiating the responses of GABA type A receptors (GABAAR), members of the superfamily of pentameric ligand-gated ion channels (pLGICs) that contain anion-selective channels. Propofol also inhibits pLGICs containing cation-selective channels, including nicotinic acetylcholine receptors and GLIC, a prokaryotic proton-gated homologue from Gloeobacter violaceus . In the structure of GLIC cocrystallized with propofol at pH 4 (presumed open/desensitized states), propofol was localized to an intrasubunit pocket at the extracellular end of the transmembrane domain within the bundle of transmembrane α-helices (Nury, H, et al. (2011) Nature 469, 428-431). To identify propofol binding sites in GLIC in solution, we used a recently developed photoreactive propofol analogue (2-isopropyl-5-[3-(trifluoromethyl)-3H-diazirin-3-yl]phenol or AziPm) that acts as an anesthetic in vivo and potentiates GABAAR in vitro. For GLIC expressed in Xenopus oocytes, propofol and AziPm inhibited current responses at pH 5.5 (EC20) with IC50 values of 20 and 50 μM, respectively. When [(3)H]AziPm (7 μM) was used to photolabel detergent-solubilized, affinity-purified GLIC at pH 4.4, protein microsequencing identified propofol-inhibitable photolabeling of three residues in the GLIC transmembrane domain: Met-205, Tyr-254, and Asn-307 in the M1, M3, and M4 transmembrane helices, respectively. Thus, for GLIC in solution, propofol and AziPm bind competitively to a site in proximity to these residues, which, in the GLIC crystal structure, are in contact with the propofol bound in the intrasubunit pocket.

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.

Inhibition versus potentiation of ligand-gated ion channels can be altered by a single mutation that moves ligands between intra- and intersubunit sites

Pentameric ligand-gated ion channels (pLGICs) are similar in structure but either inhibited or potentiated by alcohols and anesthetics. This dual modulation has previously not been understood, but the determination of X-ray structures of prokaryotic GLIC provides an ideal model system. Here, we show that a single-site mutation at the F14' site in the GLIC transmembrane domain turns desflurane and chloroform from inhibitors to potentiators, and that this is explained by competing allosteric sites. The F14'A mutation opens an intersubunit site lined by N239 (15'), I240 (16'), and Y263. Free energy calculations confirm this site is the preferred binding location for desflurane and chloroform in GLIC F14'A. In contrast, both anesthetics prefer an intrasubunit site in wild-type GLIC. Modulation is therefore the net effect of competitive binding between the intersubunit potentiating site and an intrasubunit inhibitory site. This provides direct evidence for a dual-site model of allosteric regulation of pLGICs.

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.

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.