Molecular mechanism for the dual alcohol modulation of Cys-loop receptors

Molecular Imaging Studies of Alcohol Use Disorder

Alcohol use disorder (AUD) is a serious public health problem in many countries, bringing a gamut of health risks and impairments to individuals and a great burden to society. Despite the prevalence of a disease model of AUD, the current pharmacopeia does not present reliable treatments for AUD; approved treatments are confined to a narrow spectrum of medications engaging inhibitory γ-aminobutyric acid (GABA) neurotransmission and possibly excitatory N-methyl-D-aspartate (NMDA) receptors, and opioid receptor antagonists. Molecular imaging with positron emission tomography (PET) and single-photon emission computed tomography (SPECT) can open a window into the living brain and has provided diverse insights into the pathology of AUD. In this narrative review, we summarize the state of molecular imaging findings on the pharmacological action of ethanol and the neuropathological changes associated with AUD. Laboratory and preclinical imaging results highlight the interactions between ethanol and GABA A-type receptors (GABAAR), but the interpretation of such results is complicated by subtype specificity. An abundance of studies with the glucose metabolism tracer fluorodeoxyglucose (FDG) concur in showing cerebral hypometabolism after ethanol challenge, but there is relatively little data on long-term changes in AUD. Alcohol toxicity evokes neuroinflammation, which can be tracked using PET with ligands for the microglial marker translocator protein (TSPO). Several PET studies show reversible increases in TSPO binding in AUD individuals, and preclinical results suggest that opioid-antagonists can rescue from these inflammatory responses. There are numerous PET/SPECT studies showing changes in dopaminergic markers, generally consistent with an impairment in dopamine synthesis and release among AUD patients, as seen in a number of other addictions; this may reflect the composite of an underlying deficiency in reward mechanisms that predisposes to AUD, in conjunction with acquired alterations in dopamine signaling. There is little evidence for altered serotonin markers in AUD, but studies with opioid receptor ligands suggest a specific up-regulation of the μ-opioid receptor subtype. Considerable heterogeneity in drinking patterns, gender differences, and the variable contributions of genetics and pre-existing vulnerability traits present great challenges for charting the landscape of molecular imaging in AUD.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

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.

Effect of ethanol on Munc13-1 C1 in Membrane: A Molecular Dynamics Simulation Study

BACKGROUND

EtOH has a significant effect on synaptic plasticity. Munc13-1 is an essential presynaptic active zone protein involved in priming the synaptic vesicle and releasing neurotransmitter in the brain. It is a peripheral membrane protein and binds to the activator, diacylglycerol (DAG)/phorbol ester at its membrane-targeting C1 domain. Our previous studies identified Glu-582 of C1 domain as the alcohol-binding residue (Das, J. et al, J. Neurochem., 126, 715-726, 2013).

METHODS

Here, we describe a 250 ns molecular dynamics (MD) simulation study on the interaction of EtOH and the activator-bound Munc13-1 C1 in the presence of varying concentrations of phosphatidylserine (PS).

RESULTS

In this study, Munc13-1 C1 shows higher conformational stability in EtOH than in water. It forms fewer hydrogen bonds with phorbol 13-acetate in the presence of EtOH than in water. EtOH also affected the interaction between the protein and the membrane and between the activator and the membrane. Similar studies in a E582A mutant suggest that these effects of EtOH are mostly mediated through Glu-582.

CONCLUSIONS

EtOH forms hydrogen bonds with Glu-582. While occupancy of the EtOH molecules at the vicinity (4Å) of Glu-582 is 34.4%, the occupancy in the E582A mutant is 26.5% of the simulation time. In addition, the amount of PS in the membrane influences the conformational stability of the C1 domain and interactions in the ternary complex. This study is important in providing the structural basis of EtOH's effects on synaptic plasticity.

The intriguing effect of ethanol and nicotine on acetylcholine-sensitive potassium current IKAch: Insight from a quantitative model

Recent experimental work has revealed unusual features of the effect of certain drugs on cardiac inwardly rectifying potassium currents, including the constitutively active and acetylcholine-induced components of acetylcholine-sensitive current (I_KAch). These unusual features have included alternating susceptibility of the current components to activation and inhibition induced by ethanol or nicotine applied at various concentrations, and significant correlation between the drug effect and the current magnitude measured under drug-free conditions. To explain these complex drug effects, we have developed a new type of quantitative model to offer a possible interpretation of the effect of ethanol and nicotine on the I_KAch channels. The model is based on a description of I_KAch as a sum of particular currents related to the populations of channels formed by identical assemblies of different α-subunits. Assuming two different channel populations in agreement with the two reported functional I_KAch-channels (GIRK1/4 and GIRK4), the model was able to simulate all the above-mentioned characteristic features of drug-channel interactions and also the dispersion of the current measured in different cells. The formulation of our model equations allows the model to be incorporated easily into the existing integrative models of electrical activity of cardiac cells involving quantitative description of I_KAch. We suppose that the model could also help make sense of certain observations related to the channels that do not show inward rectification. This new ionic channel model, based on a concept we call population type, may allow for the interpretation of complex interactions of drugs with ionic channels of various types, which cannot be done using the ionic channel models available so far.

Mutagenesis computer experiments in pentameric ligand-gated ion channels: the role of simulation tools with different resolution

Pentameric ligand-gated ion channels (pLGICs) are an important class of widely expressed membrane neuroreceptors, which play a crucial role in fast synaptic communications and are involved in several neurological conditions. They are activated by the binding of neurotransmitters, which trigger the transmission of an electrical signal via facilitated ion flux. They can also be activated, inhibited or modulated by a number of drugs. Mutagenesis electrophysiology experiments, with natural or unnatural amino acids, have provided a large body of functional data that, together with emerging structural information from X-ray spectroscopy and cryo-electron microscopy, are helping unravel the complex working mechanisms of these neuroreceptors. Computer simulations are complementing these mutagenesis experiments, with insights at various levels of accuracy and resolution. Here, we review how a selection of computational tools, including first principles methods, classical molecular dynamics and enhanced sampling techniques, are contributing to construct a picture of how pLGICs function and can be pharmacologically targeted to treat the disorders they are responsible for.

Alcohol reduces muscle fatigue through atomistic interactions with nicotinic receptors

Alcohol consumption affects many organs and tissues, including skeletal muscle. However, the molecular mechanism of ethanol action on skeletal muscle remains unclear. Here, using molecular dynamics simulations and single channel recordings, we show that ethanol interacts with a negatively charged amino acid within an extracellular region of the neuromuscular nicotinic acetylcholine receptor (nAChR), thereby altering its global conformation and reducing the single channel current amplitude. Charge reversal of the negatively charged amino acid abolishes the nAChR-ethanol interaction. Moreover, using transgenic animals harboring the charge-reversal mutation, ex vivo measurements of muscle force production show that ethanol counters fatigue in wild type but not homozygous αE83K mutant animals. In accord, in vivo studies of motor coordination following ethanol administration reveal an approximately twofold improvement for wild type compared to homozygous mutant animals. Together, the converging results from molecular to animal studies suggest that ethanol counters muscle fatigue through its interaction with neuromuscular nAChRs.

Molecular basis for potentiation of Cx36 gap junction channel conductance by n-alcohols and general anesthetics

In our recent study, we have demonstrated that short carbon chain n-alcohols (up to octanol) stimulated while long carbon chain n-alcohols inhibited the conductance of connexin (Cx) 36 (Cx36) gap junction (GJ) channels. In contrast, GJ channels composed of other types of Cxs all were inhibited by n-alcohols independent of their carbon chain length. To identify the putative structural domains of Cx36, responsible for the dual effect of n-alcohols, we performed structural modeling of Cx36 protein docking with hexanol and isoflurane that stimulated as well as nonanol and carbenoxolone that inhibited the conductance of Cx36 GJs and revealed their multiple common docking sites and a single pocket accessible only to hexanol and isoflurane. The pocket is located in the vicinity of three unique cysteine residues, namely C264 in the fourth, and C92 and C87 in the second transmembrane domain of the neighboring Cx36 subunits. To examine the hypothesis that disulphide bonding might be involved in the stimulatory effect of hexanol and isoflurane, we generated cysteine substitutions in Cx36 and demonstrated by a dual whole-cell patch-clamp technique that in HeLa (human cervix carcinoma cell line) and N2A (mouse neuroblastoma cell line) cells these mutations reversed the stimulatory effect of hexanol and isoflurane to inhibitory one, typical of other Cxs that lack respective cysteines and a specific docking pocket for these compounds. Our findings suggest that the stimulatory effect of hexanol and isoflurane on Cx36 GJ conductance could be achieved by re-shuffling of the inter-subunit disulphide bond between C264 and C92 to the intra-subunit one between C264 and C87.

Mapping two neurosteroid-modulatory sites in the prototypic pentameric ligand-gated ion channel GLIC

Neurosteroids are endogenous sterols that potentiate or inhibit pentameric ligand-gated ion channels (pLGICs) and can be effective anesthetics, analgesics, or anti-epileptic drugs. The complex effects of neurosteroids on pLGICs suggest the presence of multiple binding sites in these receptors. Here, using a series of novel neurosteroid-photolabeling reagents combined with top-down and middle-down mass spectrometry, we have determined the stoichiometry, sites, and orientation of binding for 3α,5α-pregnane neurosteroids in the Gloeobacter ligand-gated ion channel (GLIC), a prototypic pLGIC. The neurosteroid-based reagents photolabeled two sites per GLIC subunit, both within the transmembrane domain; one site was an intrasubunit site, and the other was located in the interface between subunits. By using reagents with photoreactive groups positioned throughout the neurosteroid backbone, we precisely map the orientation of neurosteroid binding within each site. Amino acid substitutions introduced at either site altered neurosteroid modulation of GLIC channel activity, demonstrating the functional role of both sites. These results provide a detailed molecular model of multisite neurosteroid modulation of GLIC, which may be applicable to other mammalian pLGICs.

Full mutational mapping of titratable residues helps to identify proton-sensors involved in the control of channel gating in the Gloeobacter violaceus pentameric ligand-gated ion channel

The Gloeobacter violaceus ligand-gated ion channel (GLIC) has been extensively studied by X-ray crystallography and other biophysical techniques. This provided key insights into the general gating mechanism of pentameric ligand-gated ion channel (pLGIC) signal transduction. However, the GLIC is activated by lowering the pH and the location of its putative proton activation site(s) still remain(s) unknown. To this end, every Asp, Glu, and His residue was mutated individually or in combination and investigated by electrophysiology. In addition to the mutational analysis, key mutations were structurally resolved to address whether particular residues contribute to proton sensing, or alternatively to GLIC-gating, independently of the side chain protonation. The data show that multiple residues located below the orthosteric site, notably E26, D32, E35, and D122 in the lower part of the extracellular domain (ECD), along with E222, H235, E243, and H277 in the transmembrane domain (TMD), alter GLIC activation. D122 and H235 were found to also alter GLIC expression. E35 is identified as a key proton-sensing residue, whereby neutralization of its side chain carboxylate stabilizes the active state. Thus, proton activation occurs allosterically to the orthosteric site, at the level of multiple loci with a key contribution of the coupling interface between the ECD and TMD.

Pentameric ligand-gated ion channels are an important class of receptors that are involved in many neurological diseases. They have been extensively studied but a full understanding of their mechanism of action has yet to be achieved. In an effort to bypass obstacles in the research of human receptors, bacterial versions have been used to characterize the family’s structure-function relationship. One key bacterial receptor, known as GLIC, has lead the way in structural resolution of various mechanistic states along the gating pathway, yet its activation by protons is significantly less understood than its human counterparts. To define the site(s) involved in proton gating, we systematically mutated all titratable residues near the pH₅₀ of activation: Asp, Glu, and His. We determined that a previously established His residue in the transmembrane domain is structurally important but likely plays little or no role in proton gating. We instead found that proton activation is a complex multiple loci mechanism, with the key contribution stemming from the coupling interface between the extracellular and transmembrane domain, with E35 acting as a key proton-sensing residue.

Permeating disciplines: Overcoming barriers between molecular simulations and classical structure-function approaches in biological ion transport

Ion translocation across biological barriers is a fundamental requirement for life. In many cases, controlling this process-for example with neuroactive drugs-demands an understanding of rapid and reversible structural changes in membrane-embedded proteins, including ion channels and transporters. Classical approaches to electrophysiology and structural biology have provided valuable insights into several such proteins over macroscopic, often discontinuous scales of space and time. Integrating these observations into meaningful mechanistic models now relies increasingly on computational methods, particularly molecular dynamics simulations, while surfacing important challenges in data management and conceptual alignment. Here, we seek to provide contemporary context, concrete examples, and a look to the future for bridging disciplinary gaps in biological ion transport. This article is part of a Special Issue entitled: Beyond the Structure-Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.

The Glycine Receptor-A Functionally Important Primary Brain Target of Ethanol

Identification of ethanol's (EtOH) primary molecular brain targets and determination of their functional role is an ongoing, important quest. Pentameric ligand-gated ion channels, that is, the nicotinic acetylcholine receptor, the γ-aminobutyric acid type A receptor, the 5-hydroxytryptamine₃ , and the glycine receptor (GlyR), are such targets. Here, aspects of the structure and function of these receptors and EtOH's interaction with them are briefly reviewed, with special emphasis on the GlyR and the importance of this receptor and its ligands for EtOH pharmacology. It is suggested that GlyRs are involved in (i) the dopamine-activating effect of EtOH, (ii) regulating EtOH intake, and (iii) the relapse preventing effect of acamprosate. Exploration of the GlyR subtypes involved and efforts to develop subtype specific agonists or antagonists may offer new pharmacotherapies for alcohol use disorders.

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.

Simulation of Ligand Binding to Membrane Proteins

Membrane proteins are involved in a large variety of functions. Most of these protein functions are regulated by ligand binding with diverse modes of action: agonists, partial agonists, antagonists, and allosteric modulators, potentiators and inhibitors. From the pharmacological point of view, membrane proteins are one if not the major target for drug development. However, experimental structure determination of membrane proteins in complex or in free form still represents a great challenge. Molecular dynamics (MD) simulations commonly reach the microsecond scale on membrane systems. This numerical tool is mature enough to predict and add molecular details on the different ligand-binding modes. In the present chapter, I will present the different steps to design, simulate, and analyze a MD simulation system containing a protein embedded in a membrane and surrounded by water and ligand. As an illustration, the simulation of the ligand-gated ion channel γ-aminobutyric acid type A receptor (GABA_AR) surrounded by one of its allosteric potentiators, bromoform, will be presented and discussed.

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.

Effect of ethanol at clinically relevant concentrations on atrial inward rectifier potassium current sensitive to acetylcholine

Alcohol intoxication tends to induce arrhythmias, most often the atrial fibrillation. To elucidate arrhythmogenic mechanisms related to alcohol consumption, the effect of ethanol on main components of the ionic membrane current is investigated step by step. Considering limited knowledge, we aimed to examine the effect of clinically relevant concentrations of ethanol (0.8-80 mM) on acetylcholine-sensitive inward rectifier potassium current I K(Ach). Experiments were performed by the whole-cell patch clamp technique at 23 ± 1 °C on isolated rat and guinea-pig atrial myocytes, and on expressed human Kir3.1/3.4 channels. Ethanol induced changes of I K(Ach) in the whole range of concentrations applied; the effect was not voltage dependent. The constitutively active component of I K(Ach) was significantly increased by ethanol with the maximum effect (an increase by ∼100 %) between 8 and 20 mM. The changes were comparable in rat and guinea-pig atrial myocytes and also in expressed human Kir3.1/3.4 channels (i.e., structural correlate of I K(Ach)). In the case of the acetylcholine-induced component of I K(Ach), a dual ethanol effect was apparent with a striking heterogeneity of changes in individual cells. The effect correlated with the current magnitude in control: the current was increased by eth-anol in the cells showing small current in control and vice versa. The average effect peaked at 20 mM ethanol (an increase of the current by ∼20 %). Observed changes of action potential duration agreed well with the voltage clamp data. Ethanol significantly affected both components of I K(Ach) even in concentrations corresponding to light alcohol consumption.

Inter- and Intra-Subunit Butanol/Isoflurane Sites of Action in the Human Glycine Receptor

Glycine receptors (GlyRs) mediate inhibitory neurotransmission and are targets for alcohols and anesthetics in brain. GlyR transmembrane (TM) domains contain critical residues for alcohol/anesthetic action: amino acid A288 in TM3 forms crosslinks with TM1 (I229) in the adjacent subunit as well as TM2 (S267) and TM4 (Y406, W407, I409, Y410) in the same subunit. We hypothesized that these residues may participate in intra-subunit and inter-subunit sites of alcohol/anesthetic action. The following double and triple mutants of GLRA1 cDNA (encoding human glycine receptor alpha 1 subunit) were injected into Xenopus laevis oocytes: I229C/A288C, I229C/A288C/C290S, A288C/Y406C, A288C/W407C, A288C/I409C, and A288C/Y410C along with the corresponding single mutants and wild-type GLRA1. Butanol (22 mM) or isoflurane (0.6 mM) potentiation of GlyR-mediated currents before and after application of the cysteine crosslinking agent HgCl₂ (10 μM) was measured using two-electrode voltage clamp electrophysiology. Crosslinking nearly abolished butanol and isoflurane potentiation in the I229C/A288C and I229C/A288C/C290S mutants but had no effect in single mutants or wild-type. Crosslinking also inhibited butanol and isoflurane potentiation in the TM3-4 mutants (A288C/Y406C, A288C/W407C, A288C/I409C, A288C/Y410C) with no effect in single mutants or wild-type. We extracted proteins from oocytes expressing I229C/288C, A288C/Y410C, or wild-type GlyRs, used mass spectrometry to verify their expression and possible inter-subunit dimerization, plus immunoblotting to investigate the biochemical features of proposed crosslinks. Wild-type GlyR subunits measured about 50 kDa; after crosslinking, the dimeric/monomeric 100:50 kDa band ratio was significantly increased in I229C/288C but not A288C/Y410C mutants or wild-type, providing support for TM1-3 inter-subunit and TM3-4 intra-subunit crosslinking. A GlyR homology model based on the GluCl template provides further evidence for a multi-site model for alcohol/anesthetic interaction with human GLRA1.

Identification of an Inhibitory Alcohol Binding Site in GABAA ρ1 Receptors

Alcohols inhibit γ-aminobutyric acid type A ρ1 receptor function. After introducing mutations in several positions of the second transmembrane helix in ρ1, we studied the effects of ethanol and hexanol on GABA responses using two-electrode voltage clamp electrophysiology in Xenopus laevis oocytes. The 6' mutations produced the following effects on ethanol and hexanol responses: small increase or no change (T6'M), increased inhibition (T6'V), and small potentiation (T6'Y and T6'F). The 5' mutations produced mainly increases in hexanol inhibition. Other mutations produced small (3' and 9') or no changes (2' and L277 in the first transmembrane domain) in alcohol effects. These results suggest an inhibitory alcohol binding site near the 6' position. Homology models of ρ1 receptors based on the X-ray structure of GluCl showed that the 2', 5', 6', and 9' residues were easily accessible from the ion pore, with 5' and 6' residues from neighboring subunits facing each other; L3' and L277 also faced the neighboring subunit. We tested ethanol through octanol on single and double mutated ρ1 receptors [ρ1(I15'S), ρ1(T6'Y), and ρ1(T6'Y,I15'S)] to further characterize the inhibitory alcohol pocket in the wild-type ρ1 receptor. The pocket can only bind relatively short-chain alcohols and is eliminated by introducing Y in the 6' position. Replacing the bulky 15' residue with a smaller side chain introduced a potentiating binding site, more sensitive to long-chain than to short-chain alcohols. In conclusion, the net alcohol effect on the ρ1 receptor is determined by the sum of its actions on inhibitory and potentiating sites.

Beyond Membrane Protein Structure: Drug Discovery, Dynamics and Difficulties

Most of the previous content of this book has focused on obtaining the structures of membrane proteins. In this chapter we explore how those structures can be further used in two key ways. The first is their use in structure based drug design (SBDD) and the second is how they can be used to extend our understanding of their functional activity via the use of molecular dynamics. Both aspects now heavily rely on computations. This area is vast, and alas, too large to consider in depth in a single book chapter. Thus where appropriate we have referred the reader to recent reviews for deeper assessment of the field. We discuss progress via the use of examples from two main drug target areas; G-protein coupled receptors (GPCRs) and ion channels. We end with a discussion of some of the main challenges in the area.

Evidence for α-helices in the large intracellular domain mediating modulation of the α1-glycine receptor by ethanol and Gβγ

The α1-subunit containing glycine receptors (GlyRs) is potentiated by ethanol, in part, by intracellular Gβγ actions. Previous studies have suggested that molecular requirements in the large intracellular domain are involved; however, the lack of structural data about this region has made it difficult to describe a detailed mechanism. Using circular dichroism and molecular modeling, we generated a full model of the α1-GlyR, which includes the large intracellular domain and provides new information on structural requirements for allosteric modulation by ethanol and Gβγ. The data strongly suggest the existence of an α-helical conformation in the regions near transmembrane (TM)-3 and TM4 of the large intracellular domain. The secondary structure in the N-terminal region of the large intracellular domain near TM3 appeared critical for ethanol action, and this was tested using the homologous domain of the γ2-subunit of the GABAA receptor predicted to have little helical conformation. This region of γ2 was able to bind Gβγ and form a functional channel when combined with α1-GlyR, but it was not sensitive to ethanol. Mutations in the N- and C-terminal regions introduced to replace corresponding amino acids of the α1-GlyR sequence restored the ability to be modulated by ethanol and Gβγ. Recovery of the sensitivity to ethanol was associated with the existence of a helical conformation similar to α1-GlyR, thus being an essential secondary structural requirement for GlyR modulation by ethanol and G protein.

Molecular modeling of a tandem two pore domain potassium channel reveals a putative binding site for general anesthetics

Anesthetics are thought to mediate a portion of their activity via binding to and modulation of potassium channels. In particular, tandem pore potassium channels (K2P) are transmembrane ion channels whose current is modulated by the presence of general anesthetics and whose genetic absence has been shown to confer a level of anesthetic resistance. While the exact molecular structure of all K2P forms remains unknown, significant progress has been made toward understanding their structure and interactions with anesthetics via the methods of molecular modeling, coupled with the recently released higher resolution structures of homologous potassium channels to act as templates. Such models reveal the convergence of amino acid regions that are known to modulate anesthetic activity onto a common three- dimensional cavity that forms a putative anesthetic binding site. The model successfully predicts additional important residues that are also involved in the putative binding site as validated by the results of suggested experimental mutations. Such a model can now be used to further predict other amino acid residues that may be intimately involved in the target-based structure–activity relationships that are necessary for anesthetic binding.

Atomistic models of general anesthetics for use in in silico biological studies

While small molecules have been used to induce anesthesia in a clinical setting for well over a century, a detailed understanding of the molecular mechanism remains elusive. In this study, we utilize ab initio calculations to develop a novel set of CHARMM-compatible parameters for the ubiquitous modern anesthetics desflurane, isoflurane, sevoflurane, and propofol for use in molecular dynamics (MD) simulations. The parameters generated were rigorously tested against known experimental physicochemical properties including dipole moment, density, enthalpy of vaporization, and free energy of solvation. In all cases, the anesthetic parameters were able to reproduce experimental measurements, signifying the robustness and accuracy of the atomistic models developed. The models were then used to study the interaction of anesthetics with the membrane. Calculation of the potential of mean force for inserting the molecules into a POPC bilayer revealed a distinct energetic minimum of 4–5 kcal/mol relative to aqueous solution at the level of the glycerol backbone in the membrane. The location of this minimum within the membrane suggests that anesthetics partition to the membrane prior to binding their ion channel targets, giving context to the Meyer–Overton correlation. Moreover, MD simulations of these drugs in the membrane give rise to computed membrane structural parameters, including atomic distribution, deuterium order parameters, dipole potential, and lateral stress profile, that indicate partitioning of anesthetics into the membrane at the concentration range studied here, which does not appear to perturb the structural integrity of the lipid bilayer. These results signify that an indirect, membrane-mediated mechanism of channel modulation is unlikely.

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.

Stabilization of the GluCl ligand-gated ion channel in the presence and absence of ivermectin

Improving our understanding of the mechanisms and effects of anesthetics is a critically important part of neuroscience. The currently dominant theory is that anesthetics and similar molecules act by binding to Cys-loop receptors in the postsynaptic terminal of nerve cells and potentiate or inhibit their function. Although structures for some of the most important mammalian channels have still not been determined, a number of important results have been derived from work on homologous cationic channels in bacteria. However, partly due to the lack of a nervous system in bacteria, there are a number of questions about how these results relate to higher organisms. The recent determination of a structure of the eukaryotic chloride channel, GluCl, is an important step toward accurate modeling of mammalian channels, because it is more similar in function to human Cys-loop receptors such as GABAAR or GlyR. One potential issue with using GluCl to model other receptors is the presence of the large ligand ivermectin (IVM) positioned between all five subunits. Here, we have performed a series of microsecond molecular simulations to study how the dynamics and structure of GluCl change in the presence versus absence of IVM. When the ligand is removed, subunits move at least 2 Å closer to each other compared to simulations with IVM bound. In addition, the pore radius shrinks to 1.2 Å, all of which appears to support a model where IVM binding between subunits stabilizes an open state, and that the relaxed nonIVM conformations might be suitable for modeling other channels. Interestingly, the presence of IVM also has an effect on the structure of the important loop C located at the neurotransmitter-binding pocket, which might help shed light on its partial agonist behavior.

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.

Proteomic approaches and identification of novel therapeutic targets for alcoholism

Recent studies have shown that gene regulation is far more complex than previously believed and does not completely explain changes at the protein level. Therefore, the direct study of the proteome, considerably different in both complexity and dynamicity to the genome/transcriptome, has provided unique insights to an increasing number of researchers. During the past decade, extraordinary advances in proteomic techniques have changed the way we can analyze the composition, regulation, and function of protein complexes and pathways underlying altered neurobiological conditions. When combined with complementary approaches, these advances provide the contextual information for decoding large data sets into meaningful biologically adaptive processes. Neuroproteomics offers potential breakthroughs in the field of alcohol research by leading to a deeper understanding of how alcohol globally affects protein structure, function, interactions, and networks. The wealth of information gained from these advances can help pinpoint relevant biomarkers for early diagnosis and improved prognosis of alcoholism and identify future pharmacological targets for the treatment of this addiction.

Assessment of homology templates and an anesthetic binding site within the γ-aminobutyric acid receptor

BACKGROUND

Anesthetics mediate portions of their activity via modulation of the γ-aminobutyric acid receptor (GABAaR). Although its molecular structure remains unknown, significant progress has been made toward understanding its interactions with anesthetics via molecular modeling.

METHODS

The structure of the torpedo acetylcholine receptor (nAChRα), the structures of the α4 and β2 subunits of the human nAChR, the structures of the eukaryotic glutamate-gated chloride channel (GluCl), and the prokaryotic pH-sensing channels, from Gloeobacter violaceus and Erwinia chrysanthemi, were aligned with the SAlign and 3DMA algorithms. A multiple sequence alignment from these structures and those of the GABAaR was performed with ClustalW. The Modeler and Rosetta algorithms independently created three-dimensional constructs of the GABAaR from the GluCl template. The CDocker algorithm docked a congeneric series of propofol derivatives into the binding pocket and scored calculated binding affinities for correlation with known GABAaR potentiation EC50s.

RESULTS

Multiple structure alignments of templates revealed a clear consensus of residue locations relevant to anesthetic effects except for torpedo nAChR. Within the GABAaR models generated from GluCl, the residues notable for modulating anesthetic action within transmembrane segments 1, 2, and 3 converged on the intersubunit interface between α and β subunits. Docking scores of a propofol derivative series into this binding site showed strong linear correlation with GABAaR potentiation EC50.

CONCLUSION

Consensus structural alignment based on homologous templates revealed an intersubunit anesthetic binding cavity within the transmembrane domain of the GABAaR, which showed a correlation of ligand docking scores with experimentally measured GABAaR potentiation.

GABA(A) receptor transmembrane amino acids are critical for alcohol action: disulfide cross-linking and alkyl methanethiosulfonate labeling reveal relative location of binding sites

Alcohols and inhaled anesthetics modulate GABA(A) receptor (GABA(A)R) function via putative binding sites within the transmembrane regions. The relative position of the amino acids lining these sites could be either inter- or intra-subunit. We introduced cysteines in relevant TM locations and tested the proximity of cysteine pairs using oxidizing and reducing agents to induce or break disulfide bridges between cysteines, and thus change GABA-mediated currents in wild-type and mutant α1β2γ2 GABA(A)Rs expressed in Xenopus laevis oocytes. We tested for: (i) inter-subunit cross-linking: a cysteine located in α1TM1 [either α1(Q229C) or α1(L232C)] was paired with a cysteine in different positions of β2TM2 and TM3; (ii) intra-subunit cross-linking: a cysteine located either in β2TM1 [β2(T225C)] or in TM2 [β2(N265C)] was paired with a cysteine in different locations along β2TM3. Three inter-subunit cysteine pairs and four intra-subunits cross-linked. In three intra-subunit cysteine combinations, the alcohol effect was reduced by oxidizing agents, suggesting intra-subunit alcohol binding. We conclude that the structure of the alcohol binding site changes during activation and that potentiation or inhibition by binding at inter- or intra-subunit sites is determined by the specific receptor and ligand.

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.

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.

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.