Structural basis for alcohol modulation of a pentameric ligand-gated ion channel

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.

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

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

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

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

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.

Conformational selection and adaptation to ligand binding in T4 lysozyme cavity mutants

The studies presented here explore the relationship between protein packing and molecular flexibility using ligand-binding cavity mutants of T4 lysozyme. Although previously reported crystal structures of the mutants investigated show single conformations that are similar to the WT protein, site-directed spin labeling in solution reveals additional conformational substates in equilibrium exchange with a WT-like population. Remarkably, binding of ligands, including the general anesthetic halothane shifts the population to the WT-like state, consistent with a conformational selection model of ligand binding, but structural adaptation to the ligand is also apparent in one mutant. Distance mapping with double electron-electron resonance spectroscopy and the absence of ligand binding suggest that the new substates induced by the cavity-creating mutations represent alternate packing modes in which the protein fills or partially fills the cavity with side chains, including the spin label in one case; external ligands compete with the side chains for the cavity space, stabilizing the WT conformation. The results have implications for mechanisms of anesthesia, the response of proteins to hydrostatic pressure, and protein engineering.

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.

Structural models of ligand-gated ion channels: sites of action for anesthetics and ethanol

The molecular mechanism(s) of action of anesthetic, and especially, intoxicating doses of alcohol (ethanol [EtOH]) have been of interest even before the advent of the Research Society on Alcoholism. Recent physiological, genetic, and biochemical studies have pin-pointed molecular targets for anesthetics and EtOH in the brain as ligand-gated ion channel (LGIC) membrane proteins, especially the pentameric (5 subunit) Cys-loop superfamily of neurotransmitter receptors including nicotinic acetylcholine (nAChRs), GABAA (GABAA Rs), and glycine receptors (GlyRs). The ability to demonstrate molecular and structural elements of these proteins critical for the behavioral effects of these drugs on animals and humans provides convincing evidence for their role in the drugs' actions. Amino acid residues necessary for pharmacologically relevant allosteric modulation of LGIC function by anesthetics and EtOH have been identified in these channel proteins. Site-directed mutagenesis revealed potential allosteric modulatory sites in both the trans-membrane domain (TMD) and extracellular domain (ECD). Potential sites of action and binding have been deduced from homology modeling of other LGICs with structures known from crystallography and cryo-electron microscopy studies. Direct information about ligand binding in the TMD has been obtained by photoaffinity labeling, especially in GABAA Rs. Recent structural information from crystallized procaryotic (ELIC and GLIC) and eukaryotic (GluCl) LGICs allows refinement of the structural models including evaluation of possible sites of EtOH action.

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

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

Molecular mechanism underlying ethanol activation of G-protein-gated inwardly rectifying potassium channels

Alcohol (ethanol) produces a wide range of pharmacological effects on the nervous system through its actions on ion channels. The molecular mechanism underlying ethanol modulation of ion channels is poorly understood. Here we used a unique method of alcohol-tagging to demonstrate that alcohol activation of a G-protein-gated inwardly rectifying potassium (GIRK or Kir3) channel is mediated by a defined alcohol pocket through changes in affinity for the membrane phospholipid signaling molecule phosphatidylinositol 4,5-bisphosphate. Surprisingly, hydrophobicity and size, but not the canonical hydroxyl, were important determinants of alcohol-dependent activation. Altering levels of G protein Gβγ subunits, conversely, did not affect alcohol-dependent activation, suggesting a fundamental distinction between receptor and alcohol gating of GIRK channels. The chemical properties of the alcohol pocket revealed here might extend to other alcohol-sensitive proteins, revealing a unique protein microdomain for targeting alcohol-selective therapeutics in the treatment of alcoholism and addiction.

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.

Selection for drinking in the dark alters brain gene coexpression networks

BACKGROUND

Heterogeneous stock (HS/NPT) mice have been used to create lines selectively bred in replicate for elevated drinking in the dark (DID). Both selected lines routinely reach a blood ethanol (EtOH) concentration (BEC) of 1.00 mg/ml or greater at the end of the 4-hour period of access in Day 2. The mechanisms through which genetic differences influence DID are currently unclear. Therefore, the current study examines the transcriptome, the first stage at which genetic variability affects neurobiology. Rather than focusing solely on differential expression (DE), we also examine changes in the ways that gene transcripts collectively interact with each other, as revealed by changes in coexpression patterns.

METHODS

Naïve mice (N = 48/group) were genotyped using the Mouse Universal Genotyping Array, which provided 3,683 informative markers. Quantitative trait locus (QTL) analysis used a marker-by-marker strategy with the threshold for a significant logarithm of odds (LOD) set at 10.6. Gene expression in the ventral striatum was measured using the Illumina Mouse 8.2 array. Differential gene expression and the weighted gene coexpression network analysis (WGCNA) were implemented largely as described elsewhere.

RESULTS

Significant QTLs for elevated BECs after DID were detected on chromosomes 4, 14, and 16; the latter 2 were associated with gene-poor regions. None of the QTLs overlapped with known QTLs for EtOH preference drinking. Ninety-four transcripts were detected as being differentially expressed in both selected lines versus HS controls; there was no overlap with known preference genes. The WGCNA revealed 2 modules as showing significant effects of both selections on intramodular connectivity. A number of genes known to be associated with EtOH phenotypes (e.g., Gabrg1, Glra2, Grik1, Npy2r, and Nts) showed significant changes in connectivity.

CONCLUSIONS

We found marked and consistent effects of selection on coexpression patterns; DE changes were more modest and less concordant. The QTLs and differentially expressed genes detected here are distinct from the preference phenotype. This is consistent with behavioral data and suggests that the DID and preference phenotypes are markedly different genetically.

The relative orientation of the TM3 and TM4 domains varies between α1 and α3 glycine receptors

Glycine receptors (GlyRs) are anion-conducting members of the pentameric ligand-gated ion channel family. We previously showed that the dramatic difference in glycine efficacies of α1 and α3 GlyRs is largely attributable to their nonconserved TM4 domains. Because mutation of individual nonconserved TM4 residues had little effect, we concluded that the efficacy difference was a distributed effect of all nonconserved TM4 residues. We therefore hypothesized that the TM4 domains of α1 and α3 GlyRs differ in structure, membrane orientation, and/or molecular dynamic properties. Here we employed voltage-clamp fluorometry to test whether their TM4 domains interact differently with their respective TM3 domains. We found a rhodamine fluorophore covalently attached to a homologous TM4 residue in each receptor interacts differentially with a conserved TM3 residue. We conclude that the α1 and α3 GlyR TM4 domains are orientated differently relative to their TM3 domains. This may underlie their differential ability to influence glycine efficacy.

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

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

Multisite binding of a general anesthetic to the prokaryotic pentameric Erwinia chrysanthemi ligand-gated ion channel (ELIC)

Pentameric ligand-gated ion channels (pLGICs), such as nicotinic acetylcholine, glycine, γ-aminobutyric acid GABA_A/C receptors, and the Gloeobacter violaceus ligand-gated ion channel (GLIC), are receptors that contain multiple allosteric binding sites for a variety of therapeutics, including general anesthetics. Here, we report the x-ray crystal structure of the Erwinia chrysanthemi ligand-gated ion channel (ELIC) in complex with a derivative of chloroform, which reveals important features of anesthetic recognition, involving multiple binding at three different sites. One site is located in the channel pore and equates with a noncompetitive inhibitor site found in many pLGICs. A second transmembrane site is novel and is located in the lower part of the transmembrane domain, at an interface formed between adjacent subunits. A third site is also novel and is located in the extracellular domain in a hydrophobic pocket between the β7–β10 strands. Together, these results extend our understanding of pLGIC modulation and reveal several specific binding interactions that may contribute to modulator recognition, further substantiating a multisite model of allosteric modulation in this family of ion channels.

Distinct sensitivity of slo1 channel proteins to ethanol

Ethanol levels reached in circulation during moderate-to-heavy alcohol intoxication (50-100 mM) modify Ca(2+)- and voltage-gated K(+) (BK) channel steady-state activity, eventually altering both physiology and behavior. Ethanol action on BK steady-state activity solely requires the channel-forming subunit slo1 within a bare lipid environment. To identify the protein regions that confer ethanol sensitivity to slo1, we tested the ethanol sensitivity of heterologously expressed slo1 and structurally related channels. Ethanol (50 mM) increased the steady-state activities of mslo1 and Ca(2+)-gated MthK, the latter after channel reconstitution into phospholipid bilayers. In contrast, 50-100 mM ethanol failed to alter the steady-state activities of Na(+)/Cl(-)-gated rslo2, H(+)-gated mslo3, and an mslo1/3 chimera engineered by joining the mslo1 region encompassing the N terminus to S6 with the mslo3 cytosolic tail domain (CTD). Collectively, data indicate that the slo family canonical design, which combines a transmembrane 6 (TM6) voltage-gated K(+) channel (K(V)) core with CTDs that empower the channel with ion-sensing, does not necessarily render ethanol sensitivity. In addition, the region encompassing the N terminus to the S0-S1 cytosolic loop (missing in MthK) is not necessary for ethanol action. Moreover, incorporation of both this region and an ion-sensing CTD to TM6 K(V) cores (a design common to mslo1, mslo3, and the mslo1/mslo3 chimera) is not sufficient for ethanol sensitivity. Rather, a CTD containing Ca(2+)-sensing regulator of conductance for K(+) domains seems to be critical to bestow K(V) structures, whether of TM2 (MthK) or TM6 (slo1), with sensitivity to intoxicating ethanol levels.

Modeling complex biological systems: From solution chemistry to membranes and channels

Complex biological systems are intimately linked to their environment, a very crowded and equally complex solution compartmentalized by fluid membranes. Modeling such systems remains challenging and requires a suitable representation of these solutions and their interfaces. Here, we focus on particle-based modeling at an atomistic level using molecular dynamics (MD) simulations. As an example, we discuss important steps in modeling the solution chemistry of an ion channel of the ligand-gated ion channel receptor family, a major target of many drugs including anesthetics and addiction treatments. The bacterial pentameric ligand-gated ion channel (pLGIC) called GLIC provides clues about the functional importance of solvation, in particular for mechanisms such as permeation and gating. We present some current challenges along with promising novel modeling approaches.

Structure and pharmacology of pentameric receptor channels: from bacteria to brain

Orthologs of the pentameric receptor channels that mediate fast synaptic transmission in the central and peripheral nervous systems have been found in several bacterial species and in a single archaea genus. Recent X-ray structures of bacterial and invertebrate pentameric receptors point to a striking conservation of the structural features within the whole family, even between distant prokaryotic and eukaryotic members. These structural data reveal general principles of molecular organization that allow allosteric membrane proteins to mediate chemoelectric transduction. Notably, several conformations have been solved, including open and closed channels with distinct global tertiary and quaternary structure. The data reveal features of the ion channel architecture and of diverse categories of binding sites, such as those that bind orthosteric ligands, including neurotransmitters, and those that bind allosteric modulators, such as general anesthetics, ivermectin, or lipids. In this review, we summarize the most recent data, discuss insights into the mechanism of action in these systems, and elaborate on newly opened avenues for drug design.

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

Cys-loop receptors constitute a superfamily of pentameric ligand-gated ion channels (pLGICs), including receptors for acetylcholine, serotonin, glycine and γ-aminobutyric acid. Several bacterial homologues have been identified that are excellent models for understanding allosteric binding of alcohols and anesthetics in human Cys-loop receptors. Recently, we showed that a single point mutation on a prokaryotic homologue (GLIC) could transform it from a channel weakly potentiated by ethanol into a highly ethanol-sensitive channel. Here, we have employed molecular simulations to study ethanol binding to GLIC, and to elucidate the role of the ethanol-enhancing mutation in GLIC modulation. By performing 1-µs simulations with and without ethanol on wild-type and mutated GLIC, we observed spontaneous binding in both intra-subunit and inter-subunit transmembrane cavities. In contrast to the glycine receptor GlyR, in which we previously observed ethanol binding primarily in an inter-subunit cavity, ethanol primarily occupied an intra-subunit cavity in wild-type GLIC. However, the highly ethanol-sensitive GLIC mutation significantly enhanced ethanol binding in the inter-subunit cavity. These results demonstrate dramatic effects of the F(14′)A mutation on the distribution of ligands, and are consistent with a two-site model of pLGIC inhibition and potentiation.

Communication from one nerve cell to the next is an essential process for brain and muscle function. Nerve impulses result in release of transmitter molecules from one cell that bind to receptors on the next cell. Transmitter binding opens a pore in each receptor and ions flow across the membrane, leading to either enhancement or inhibition of new nerve impulses. These receptors are modulated by numerous drugs, including alcohols and anesthetics; identifying the precise location of modulator binding is critical for drug development. We have used computer simulation methods to model alcohol diffusion and binding to a receptor. By modifying a single residue in the receptor, we were able to move the location of the binding site and dramatically alter alcohol modulation, which supports a model with two separate binding sites for enhancement and inhibition in this family of receptors.

Conformational transitions underlying pore opening and desensitization in membrane-embedded Gloeobacter violaceus ligand-gated ion channel (GLIC)

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

Structure of the pentameric ligand-gated ion channel GLIC bound with anesthetic ketamine

Pentameric ligand-gated ion channels (pLGICs) are targets of general anesthetics, but a structural understanding of anesthetic action on pLGICs remains elusive. GLIC, a prokaryotic pLGIC, can be inhibited by anesthetics, including ketamine. The ketamine concentration leading to half-maximal inhibition of GLIC (58 μM) is comparable to that on neuronal nicotinic acetylcholine receptors. A 2.99 Å resolution X-ray structure of GLIC bound with ketamine revealed ketamine binding to an intersubunit cavity that partially overlaps with the homologous antagonist-binding site in pLGICs. The functional relevance of the identified ketamine site was highlighted by profound changes in GLIC activation upon cysteine substitution of the cavity-lining residue N152. The relevance is also evidenced by changes in ketamine inhibition upon the subsequent chemical labeling of N152C. The results provide structural insight into the molecular recognition of ketamine and are valuable for understanding the actions of anesthetics and other allosteric modulators on pLGICs.

General anesthetics predicted to block the GLIC pore with micromolar affinity

Although general anesthetics are known to modulate the activity of ligand-gated ion channels in the Cys-loop superfamily, there is at present neither consensus on the underlying mechanisms, nor predictive models of this modulation. Viable models need to offer quantitative assessment of the relative importance of several identified anesthetic binding sites. However, to date, precise affinity data for individual sites has been challenging to obtain by biophysical means. Here, the likely role of pore block inhibition by the general anesthetics isoflurane and propofol of the prokaryotic pentameric channel GLIC is investigated by molecular simulations. Microscopic affinities are calculated for both single and double occupancy binding of isoflurane and propofol to the GLIC pore. Computations are carried out for an open-pore conformation in which the pore is restrained to crystallographic radius, and a closed-pore conformation that results from unrestrained molecular dynamics equilibration of the structure. The GLIC pore is predicted to be blocked at the micromolar concentrations for which inhibition by isofluorane and propofol is observed experimentally. Calculated affinities suggest that pore block by propofol occurs at signifcantly lower concentrations than those for which inhibition is observed: we argue that this discrepancy may result from binding of propofol to an allosteric site recently identified by X-ray crystallography, which may cause a competing gain-of-function effect. Affinities of isoflurane and propofol to the allosteric site are also calculated, and shown to be 3 mM for isoflurane and for propofol; both anesthetics have a lower affinity for the allosteric site than for the unoccupied pore.

Although general anesthesia is performed every day on thousands of people, its detailed microscopic mechanisms are not known. What is known is that general anesthetic drugs modulate the activity of ion channels in the central nervous system. These channels are proteins that open in response to binding of neurotransmitter molecules, creating an electric current through the cell membrane and thus propagating nerve impulses between cells. One possible mechanism for ion channel inhibition by anesthetics is that the drugs bind inside the pore of the channels, blocking ion current. Here we investigate such a pore block mechanism by computing the strength of the drugs' interaction with the pore – and hence the likelihood of binding, in the case of GLIC, a bacterial channel protein. The results, obtained from numerical simulations of atomic models of GLIC, indicate that the anesthetics isoflurane and propofol have a tendency to bind in the pore that is strong enough to explain blocking of the channel, even at low concentration of the drugs.

Discovery of a novel allosteric modulator of 5-HT3 receptors: inhibition and potentiation of Cys-loop receptor signaling through a conserved transmembrane intersubunit site

The ligand-gated ion channels in the Cys-loop receptor superfamily mediate the effects of neurotransmitters acetylcholine, serotonin, GABA, and glycine. Cys-loop receptor signaling is susceptible to modulation by ligands acting through numerous allosteric sites. Here we report the discovery of a novel class of negative allosteric modulators of the 5-HT(3) receptors (5-HT(3)Rs). PU02 (6-[(1-naphthylmethyl)thio]-9H-purine) is a potent and selective antagonist displaying IC(50) values of ~1 μM at 5-HT(3)Rs and substantially lower activities at other Cys-loop receptors. In an elaborate mutagenesis study of the 5-HT(3)A receptor guided by a homology model, PU02 is demonstrated to act through a transmembrane intersubunit site situated in the upper three helical turns of TM2 and TM3 in the (+)-subunit and TM1 and TM2 in the (-)-subunit. The Ser(248), Leu(288), Ile(290), Thr(294), and Gly(306) residues are identified as important molecular determinants of PU02 activity with minor contributions from Ser(292) and Val(310), and we propose that the naphthalene group of PU02 docks into the hydrophobic cavity formed by these. Interestingly, specific mutations of Ser(248), Thr(294), and Gly(306) convert PU02 into a complex modulator, potentiating and inhibiting 5-HT-evoked signaling through these mutants at low and high concentrations, respectively. The PU02 binding site in the 5-HT(3)R corresponds to allosteric sites in anionic Cys-loop receptors, which emphasizes the uniform nature of the molecular events underlying signaling through the receptors. Moreover, the dramatic changes in the functional properties of PU02 induced by subtle changes in its binding site bear witness to the delicate structural discrimination between allosteric inhibition and potentiation of Cys-loop receptors.

Desensitization mechanism in prokaryotic ligand-gated ion channel

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

Characterization of two mutations, M287L and Q266I, in the α1 glycine receptor subunit that modify sensitivity to alcohols

Glycine receptors (GlyRs) are inhibitory ligand-gated ion channels. Ethanol potentiates glycine activation of the GlyR, and putative binding sites for alcohol are located in the transmembrane (TM) domains between and within subunits. To alter alcohol sensitivity of GlyR, we introduced two mutations in the GlyR α1 subunit, M287L (TM3) and Q266I (TM2). After expression in Xenopus laevis oocytes, both mutants showed a reduction in glycine sensitivity and glycine-induced maximal currents. Activation by taurine, another endogenous agonist, was almost abolished in the M287L GlyR. The ethanol potentiation of glycine currents was reduced in the M287L GlyR and eliminated in Q266I. Physiological levels of zinc (100 nM) potentiate glycine responses in wild-type GlyR and also enhance the ethanol potentiation of glycine responses. Although zinc potentiation of glycine responses was unchanged in both mutants, zinc enhancement of ethanol potentiation of glycine responses was absent in M287L GlyRs. The Q266I mutation decreased conductance but increased mean open time (effects not seen in M287L). Two lines of knockin mice bearing these mutations were developed. Survival of homozygous knockin mice was impaired, probably as a consequence of impaired glycinergic transmission. Glycine showed a decreased capacity for displacing strychnine binding in heterozygous knockin mice. Electrophysiology in isolated neurons of brain stem showed decreased glycine-mediated currents and decreased ethanol potentiation in homozygous knockin mice. Molecular models of the wild-type and mutant GlyRs show a smaller water-filled cavity within the TM domains of the Q266I α1 subunit. The behavioral characterization of these knockin mice is presented in a companion article (J Pharmacol Exp Ther 340:317-329, 2012).

The TM2 6' position of GABA(A) receptors mediates alcohol inhibition

Ionotropic GABA(A) receptors (GABA(A)Rs), which mediate inhibitory neurotransmission in the central nervous system, are implicated in the behavioral effects of alcohol and alcoholism. Site-directed mutagenesis studies support the presence of discrete molecular sites involved in alcohol enhancement and, more recently, inhibition of GABA(A)Rs. We used Xenopus laevis oocytes to investigate the 6' position in the second transmembrane region of GABA(A)Rs as a site influencing alcohol inhibition. We asked whether modification of the 6' position by substitution with larger residues or methanethiol labeling [using methyl methanethiosulfonate (MMTS)] of a substituted cysteine, reduced GABA action and/or blocked further inhibition by alcohols. Labeling of the 6' position in either α2 or β2 subunits reduced responses to GABA. In addition, methanol and ethanol potentiation increased after MMTS labeling or substitution with tryptophan or methionine, consistent with elimination of an inhibitory site for these alcohols. Specific alcohols, but not the anesthetic etomidate, competed with MMTS labeling at the 6' position. We verified a role for the 6' position in previously tested α2β2 as well as more physiologically relevant α2β2γ2s GABA(A)Rs. Finally, we built a novel molecular model based on the invertebrate glutamate-gated chloride channel receptor, a GABA(A)R homolog, revealing that the 6' position residue faces the channel pore, and modification of this residue alters volume and polarity of the pore-facing cavity in this region. These results indicate that the 6' positions in both α2 and β2 GABA(A)R subunits mediate inhibition by short-chain alcohols, which is consistent with the presence of multiple counteracting sites of action for alcohols on ligand-gated ion channels.

Allostery and the Monod-Wyman-Changeux model after 50 years

The Monod-Wyman-Changeux (MWC) model was conceived in 1965 to account for the signal transduction and cooperative properties of bacterial regulatory enzymes and hemoglobin. It was soon extended to pharmacological receptors for neurotransmitters and other macromolecular entities involved in intracellular and intercellular communications. Five decades later, the two main hypotheses of the model are reexamined on the basis of a variety of regulatory proteins with known X-ray structures: (a) Regulatory proteins possess an oligomeric structure with symmetry properties, and (b) the allosteric interactions between topographically distinct sites are mediated by a conformational transition established between a few preestablished states with conservation of symmetry and ligand-directed conformational selection. Several well-documented examples are adequately represented by the MWC model, yet a few possible exceptions are noted. New questions are raised concerning the dynamics of the allosteric transitions and more complex supramolecular ensembles.

Cys-loop receptor channel blockers also block GLIC

The Gloeobacter ligand-gated ion channel (GLIC) is a bacterial homolog of vertebrate Cys-loop ligand-gated ion channels. Its pore-lining region in particular has a high sequence homology to these related proteins. Here we use electrophysiology to examine a range of compounds that block the channels of Cys-loop receptors to probe their pharmacological similarity with GLIC. The data reveal that a number of these compounds also block GLIC, although the pharmacological profile is distinct from these other proteins. The most potent compound was lindane, a GABA(A) receptor antagonist, with an IC₅₀ of 0.2 μM. Docking studies indicated two potential binding sites for this ligand in the pore, at the 9' or between the 0' and 2' residues. Similar experiments with picrotoxinin (IC₅₀ = 2.6 μM) and rimantadine (IC₅₀ = 2.6 μM) reveal interactions with 2'Thr residues in the GLIC pore. These locations are strongly supported by mutagenesis data for picrotoxinin and lindane, which are less potent in a T2'S version of GLIC. Overall, our data show that the inhibitory profile of the GLIC pore has considerable overlap with those of Cys-loop receptors, but the GLIC pore has a unique pharmacology.

Intramembrane proton binding site linked to activation of bacterial pentameric ion channel

Prokaryotic orthologs of eukaryotic Cys-loop receptor channels recently emerged as structural and mechanistic surrogates to investigate this superfamily of intercellular signaling proteins. Here, we examine proton activation of the prokaryotic ortholog GLIC using patch clamp electrophysiology, mutagenesis, and molecular dynamics (MD) simulations. Whole-cell current recordings from human embryonic kidney (HEK) 293 cells expressing GLIC show half-maximal activation at pH 6, close to the pK(a) of histidine, implicating the three native His residues in proton sensing linked to activation. The mutation H235F abolishes proton activation, H277Y is without effect, and all nine mutations of His-127 prevent expression on the cell surface. In the GLIC crystal structure, His-235 on transmembrane (TM) α-helix 2, hydrogen bonds to the main chain carbonyl oxygen of Ile-259 on TM α-helix 3. MD simulations show that when His-235 is protonated, the hydrogen bond persists, and the channel remains in the open conformation, whereas when His-235 is deprotonated, the hydrogen bond dissociates, and the channel closes. Mutations of the proximal Tyr-263, which also links TM α-helices 2 and 3 via a hydrogen bond, alter proton sensitivity over a 1.5 pH unit range. MD simulations show that mutations of Tyr-263 alter the hydrogen bonding capacity of His-235. The overall findings show that His-235 in the TM region of GLIC is a novel proton binding site linked to channel activation.

Functional prokaryotic-eukaryotic chimera from the pentameric ligand-gated ion channel family

Pentameric ligand-gated ion channels (pLGICs), which mediate chemo-electric signal transduction in animals, have been recently found in bacteria. Despite clear sequence and 3D structure homology, the phylogenetic distance between prokaryotic and eukaryotic homologs suggests significant structural divergences, especially at the interface between the extracellular (ECD) and the transmembrane (TMD) domains. To challenge this possibility, we constructed a chimera in which the ECD of the bacterial protein GLIC is fused to the TMD of the human α1 glycine receptor (α1GlyR). Electrophysiology in Xenopus oocytes shows that it functions as a proton-gated ion channel, thereby locating the proton activation site(s) of GLIC in its ECD. Patch-clamp experiments in BHK cells show that the ion channel displays an anionic selectivity with a unitary conductance identical to that of the α1GlyR. In addition, pharmacological investigations result in transmembrane allosteric modulation similar to the one observed on α1GlyR. Indeed, the clinically active drugs propofol, four volatile general anesthetics, alcohols, and ivermectin all potentiate the chimera while they inhibit GLIC. Collectively, this work shows the compatibility between GLIC and α1GlyR domains and points to conservation of the ion channel and transmembrane allosteric regulatory sites in the chimera. This provides evidence that GLIC and α1GlyR share a highly homologous 3D structure. GLIC is thus a relevant model of eukaryotic pLGICs, at least from the anionic type. In addition, the chimera is a good candidate for mass production in Escherichia coli, opening the way for investigations of "druggable" eukaryotic allosteric sites by X-ray crystallography.