Anxiety over GABA(A) receptor structure relieved by AChBP

Rationalizing the binding and α subtype selectivity of synthesized imidazodiazepines and benzodiazepines at GABAA receptors by using molecular docking studies

The pharmacological actions exerted by benzodiazepines are dependent on the discrete α protein subunits of the γ-aminobutyric acid type A receptor (GABA_A R). Recent developments via a cryo-EM structure of the α1β3γ2L GABA_A R ion channel provide crucial insights into ligand efficacy and binding affinity at this subtype. We investigated the molecular interactions of diazepam and alprazolam bound GABA_A R structures (6HUP and 6HUO) to determine key binding interaction domains. A halogen bond between the chlorine atoms of diazepam and alprazolam with the group on the backbone of the α1 histidine amino acid 102 is important to the positive allosteric modulatory actions of diazepam and alprazolam in the α1β3γ2L GABA_A R ion channel. In order to gain insight into α subtype selectivity we designed and synthesized close structural analogs of diazepam and alprazolam. These compounds were then docked into the recently publish cryo-EM structures of GABA_A Rs (6HUP and 6HUO). This modeling along with radio-ligand binding data resulted in the conclusion that the non-classical bioisosteric replacement of the chlorine atom at C7 with an ethinyl group (compound 5) resulted in an 11-fold gain in α5 binding selectivity over the α1 subtype. Moreover, the potency of compound 5 resulted in a ligand with less sedation than diazepam, while still maintaining the same anxiolytic potency. These modeling data extend our understanding of the structural requirements for α-subtype-selective compounds that can be utilized to achieve improved medical treatments. It is clear that the ethinyl group in place of a halogen atom decreases the affinity and efficacy of benzodiazepines and imidazodiazepines at α1 subtypes, which results in less sedation and ataxia.

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

Traditional drug discovery is a lengthy process which involves a huge amount of resources. Modern-day drug discovers various multidisciplinary approaches amongst which, computational ligand and structure-based drug designing methods contribute significantly. Structure-based drug designing techniques require the knowledge of structural information of drug target and drug-target complexes. Proper understanding of drug-target binding requires the flexibility of both ligand and receptor to be incorporated. Molecular docking refers to the static picture of the drug-target complex(es). Molecular dynamics, on the other hand, introduces flexibility to understand the drug binding process.

OBJECTIVE

The aim of the present study is to provide a systematic review on the usage of molecular dynamics simulations to aid the process of structure-based drug design.

METHOD

This review discussed findings from various research articles and review papers on the use of molecular dynamics in drug discovery. All efforts highlight the practical grounds for which molecular dynamics simulations are used in drug designing program. In summary, various aspects of the use of molecular dynamics simulations that underline the basis of studying drug-target complexes were thoroughly explained.

RESULTS

This review is the result of reviewing more than a hundred papers. It summarizes various problems that use molecular dynamics simulations.

CONCLUSION

The findings of this review highlight how molecular dynamics simulations have been successfully implemented to study the structure-function details of specific drug-target complexes. It also identifies the key areas such as stability of drug-target complexes, ligand binding kinetics and identification of allosteric sites which have been elucidated using molecular dynamics simulations.

Exploiting the Indole Scaffold to Design Compounds Binding to Different Pharmacological Targets

Several indole derivatives have been disclosed by our research groups that have been collaborating for nearly 25 years. The results of our investigations led to a variety of molecules binding selectively to different pharmacological targets, specifically the type A γ-aminobutyric acid (GABA_A) chloride channel, the translocator protein (TSPO), the murine double minute 2 (MDM2) protein, the A_2B adenosine receptor (A_2B AR) and the Kelch-like ECH-associated protein 1 (Keap1). Herein, we describe how these works were conceived and carried out thanks to the versatility of indole nucleus to be exploited in the design and synthesis of drug-like molecules.

GABAA Receptor β2E155 Residue Located at the Agonist-Binding Site Is Involved in the Receptor Gating

GABA_A receptors (GABA_ARs) play a crucial role in mediating inhibition in the adult brain. In spite of progress in describing (mainly) the static structures of this receptor, the molecular mechanisms underlying its activation remain unclear. It is known that in the α₁β₂γ_2L receptors, the mutation of the β₂E155 residue, at the orthosteric binding site, strongly impairs the receptor activation, but the molecular and kinetic mechanisms of this effect remain elusive. Herein, we investigated the impact of the β₂E155C mutation on binding and gating of the α₁β₂γ_2L receptor. To this end, we combined the macroscopic and single-channel analysis, the use of different agonists [GABA and muscimol (MSC)] and flurazepam (FLU) as a modulator. As expected, the β₂E155C mutation caused a vast right shift of the dose–response (for GABA and MSC) and, additionally, dramatic changes in the time course of current responses, indicative of alterations in gating. Mutated receptors showed reduced maximum open probability and enhanced receptor spontaneous activity. Model simulations for macroscopic currents revealed that the primary effect of the mutation was the downregulation of the preactivation (flipping) rate. Experiments with MSC and FLU further confirmed a reduction in the preactivation rate. Our single-channel analysis revealed the mutation impact mainly on the second component in the shut times distributions. Based on model simulations, this finding further confirms that this mutation affects mostly the preactivation transition, supporting thus the macroscopic data. Altogether, we provide new evidence that the β₂E155 residue is involved in both binding and gating (primarily preactivation).

GABAA receptor family: overview on structural characterization

The pentameric γ-aminobutyric acid type A receptors are ion channels activated by ligands, which intervene in the rapid inhibitory transmission in the mammalian CNS. Due to their rich pharmacology and therapeutic potential, it is essential to understand their structure and function thoroughly. This deep characterization was hampered by the lack of experimental structural information for many years. Thus, computational techniques have been extensively combined with experimental data, in order to undertake the study of γ-aminobutyric acid type A receptors and their interaction with drugs. Here, we review the exciting journey made to assess the structures of these receptors and outline major outcomes. Finally, we discuss the brand new structure of the α1β2γ2 subtype and the amazing advances it brings to the field.

The role of loops B and C in determining the potentiation of GABAA receptors by midazolam

Many benzodiazepines are positive allosteric modulators (PAMs) of GABAA receptors that cause sedation, hypnosis, and anxiolysis. Benzodiazepines bind GABAA receptors at the extracellular interface of the α and γ subunits. Within the α subunit, the benzodiazepine binding site is defined by three highly conserved structural loops, loops A-C. Although previous mutagenesis studies have identified His102 in Loop A as important for benzodiazepine modulation of GABAA receptors, the functional roles of many of the other conserved residues in loops A-C remain incompletely understood. In this study, we made single mutations in loops A-C of the benzodiazepine binding-site across all six α subunits. We used whole-cell patch clamp recording to measure the functional effects of these mutations on midazolam potentiation. The results showed that mutating the threonine in loop B and serine in loop C (Thr163 and S206 in human α1) did not abolish the receptors' responsiveness to midazolam, as the α1(H102R) mutation did. The loop C mutations exhibited a novel array of α-isoform specific effects on midazolam potentiation. The α3(S230I) and α5(S209I) mutations had the largest effect on midazolam potentiation, increasing the efficacy of midazolam. Novel benzodiazepines targeting loop C may represent a future direction for designing new drugs that specifically alter the activity of α3- and α5-containing GABAA receptors.

GABAA receptor: Positive and negative allosteric modulators

gamma-Aminobutyric acid (GABA)-mediated inhibitory neurotransmission and the gene products involved were discovered during the mid-twentieth century. Historically, myriad existing nervous system drugs act as positive and negative allosteric modulators of these proteins, making GABA a major component of modern neuropharmacology, and suggesting that many potential drugs will be found that share these targets. Although some of these drugs act on proteins involved in synthesis, degradation, and membrane transport of GABA, the GABA receptors Type A (GABAAR) and Type B (GABABR) are the targets of the great majority of GABAergic drugs. This discovery is due in no small part to Professor Norman Bowery. Whereas the topic of GABABR is appropriately emphasized in this special issue, Norman Bowery also made many insights into GABAAR pharmacology, the topic of this article. GABAAR are members of the ligand-gated ion channel receptor superfamily, a chloride channel family of a dozen or more heteropentameric subtypes containing 19 possible different subunits. These subtypes show different brain regional and subcellular localization, age-dependent expression, and potential for plastic changes with experience including drug exposure. Not only are GABAAR the targets of agonist depressants and antagonist convulsants, but most GABAAR drugs act at other (allosteric) binding sites on the GABAAR proteins. Some anxiolytic and sedative drugs, like benzodiazepine and related drugs, act on GABAAR subtype-dependent extracellular domain sites. General anesthetics including alcohols and neurosteroids act at GABAAR subunit-interface trans-membrane sites. Ethanol at high anesthetic doses acts on GABAAR subtype-dependent trans-membrane domain sites. Ethanol at low intoxicating doses acts at GABAAR subtype-dependent extracellular domain sites. Thus GABAAR subtypes possess pharmacologically specific receptor binding sites for a large group of different chemical classes of clinically important neuropharmacological agents. This article is part of the "Special Issue Dedicated to Norman G. Bowery".

Orthosteric and benzodiazepine cavities of the α1β2γ2 GABAA receptor: insights from experimentally validated in silico methods

γ-aminobutyric acid-type A (GABA_A) receptors mediate fast synaptic inhibition in the central nervous system of mammals. They are modulated via several sites by numerous compounds, which include GABA, benzodiazepines, ethanol, neurosteroids and anaesthetics among others. Due to their potential as targets of novel drugs, a detailed knowledge of their structure-function relationships is needed. Here, we present the model of the α₁β₂γ₂ subtype GABA_A receptor in the APO state and in complex with selected ligands, including agonists, antagonists and allosteric modulators. The model is based on the crystallographic structure of the human β₃ homopentamer GABA_A receptor. The complexes were refined using atomistic molecular dynamics simulations. This allowed a broad description of the binding modes and the detection of important interactions in agreement with experimental information. From the best of our knowledge, this is the only model of the α₁β₂γ₂ GABA_A receptor that represents altogether the desensitized state of the channel and comprehensively describes the interactions of ligands of the orthosteric and benzodiazepines binding sites in agreement with the available experimental data. Furthermore, it is able to explain small differences regarding the binding of a variety of chemically divergent ligands. Finally, this new model may pave the way for the design of focused experimental studies that will allow a deeper description of the receptor.

γ-Aminobutyric acid (GABA) signalling in plants

The role of γ-aminobutyric acid (GABA) as a signal in animals has been documented for over 60 years. In contrast, evidence that GABA is a signal in plants has only emerged in the last 15 years, and it was not until last year that a mechanism by which this could occur was identified-a plant 'GABA receptor' that inhibits anion passage through the aluminium-activated malate transporter family of proteins (ALMTs). ALMTs are multigenic, expressed in different organs and present on different membranes. We propose GABA regulation of ALMT activity could function as a signal that modulates plant growth, development, and stress response. In this review, we compare and contrast the plant 'GABA receptor' with mammalian GABAA receptors in terms of their molecular identity, predicted topology, mode of action, and signalling roles. We also explore the implications of the discovery that GABA modulates anion flux in plants, its role in signal transduction for the regulation of plant physiology, and predict the possibility that there are other GABA interaction sites in the N termini of ALMT proteins through in silico evolutionary coupling analysis; we also explore the potential interactions between GABA and other signalling molecules.

The dual modulatory effects of efavirenz on GABAA receptors are mediated via two distinct sites

Efavirenz is a widely prescribed medicine used to treat type 1 human immunodeficiency virus (HIV-1), the most prevalent pathogenic strain of the virus responsible for the acquired immune deficiency syndrome (AIDS) pandemic. Under prescribed dosing conditions, either alone or in combination therapy, efavirenz-induced CNS disturbances are frequently reported. Efavirenz was recently reported to interact in a similar concentration range with a number of receptors, transporters and ion channels including recombinant rat α1β2γ2 GABA_A receptors whose actions were potentiated (Gatch et al., 2013; Dalwadi et al., 2016). Now we report on the molecular mechanism of efavirenz on GABA_A receptors as a function of concentration and subunit composition via whole-cell recordings of GABA-activated currents from HEK293 cells expressing varying subunit configurations of GABA_A receptors. Efavirenz elicited dual effects on the GABA response; it allosterically potentiated currents at low concentrations, whereas it inhibited currents at higher concentrations. The allosteric potentiating action on GABA_A receptors was pronounced in the α1β2γ2, α2β2γ2 and α4β2γ2 configurations, greatly diminished in the α6β2γ2 configuration, and completely absent in the α3β2γ2 or α5β2γ2 configuration. In stark contrast, the inhibitory modulation of efavirenz at higher concentrations was evident in all subunit configurations examined. Moreover, efavirenz-induced modulatory effects were dependent on GABA concentration ([GABA]), with a pronounced impact on currents activated by low [GABA] but little effect at saturating [GABA]. Mutation of a highly-conserved threonine to phenylalanine in transmembrane domain 2 of the α1 subunit abolished the inhibitory effect of efavirenz in α1β2 receptors. Finally, mutations of any of the three conserved extracellular residues in α1/2/4 subunits to the conserved residues at the corresponding positions in α3/5 subunits (i.e., R84P, M89L or I120L) completely eliminated the potentiating effect of efavirenz in α1β2γ2 configuration. These findings demonstrate that efavirenz's positive allosteric modulation of the GABA_A receptor is mediated via a novel allosteric site associated with the extracellular domain of the receptor.

Loop G in the GABAA receptor α1 subunit influences gating efficacy

KEY POINTS

The functional importance of residues in loop G of the GABA_A receptor has not been investigated. D43 and T47 in the α1 subunit are of particular significance as their structural modification inhibits activation by GABA. While the T47C substitution had no significant effect, non-conservative substitution of either residue (D43C or T47R) reduced the apparent potency of GABA. Propofol potentiated maximal GABA-evoked currents mediated by α1(D43C)β2γ2 and α1(T47R)β2γ2 receptors. Non-stationary variance analysis revealed a reduction in maximal GABA-evoked P_open , suggesting impaired agonist efficacy. Further analysis of α1(T47R)β2γ2 receptors revealed that the efficacy of the partial agonist THIP (4,5,6,7-tetrahydroisoxazolo[5,4-c]pyridine-3-ol) relative to GABA was impaired. GABA-, THIP- and propofol-evoked currents mediated by α1(T47R)β2γ2 receptors deactivated faster than those mediated by α1β2γ2 receptors, indicating that the mutation impairs agonist-evoked gating. Spontaneous gating caused by the β2(L285R) mutation was also reduced in α1(T47R)β2(L285R)γ2 compared to α1β2(L285R)γ2 receptors, confirming that α1(T47R) impairs gating independently of agonist activation.

ABSTRACT

The modification of cysteine residues (substituted for D43 and T47) by 2-aminoethyl methanethiosulfonate in the GABA_A α1 subunit loop G is known to impair activation of α1β2γ2 receptors by GABA and propofol. While the T47C substitution had no significant effect, non-conservative substitution of either residue (D43C or T47R) reduced the apparent potency of GABA. Propofol (1 μm), which potentiates sub-maximal but not maximal GABA-evoked currents mediated by α1β2γ2 receptors, also potentiated maximal currents mediated by α1(D43C)β2γ2 and α1(T47R)β2γ2 receptors. Furthermore, the peak open probabilities of α1(D43C)β2γ2 and α1(T47R)β2γ2 receptors were reduced. The kinetics of macroscopic currents mediated by α1(D43C)β2γ2 and α1(T47R)β2γ2 receptors were characterised by slower desensitisation and faster deactivation. Similar changes in macroscopic current kinetics, together with a slower activation rate, were observed with the loop D α1(F64C) substitution, known to impair both efficacy and agonist binding, and when the partial agonist THIP (4,5,6,7-tetrahydroisoxazolo[5,4-c]pyridine-3-ol) was used to activate WT or α1(T47R)β2γ2 receptors. Propofol-evoked currents mediated by α1(T47R)β2γ2 and α1(F64C)β2γ2 receptors also exhibited faster deactivation than their WT counterparts, revealing that these substitutions impair gating through a mechanism independent of orthosteric binding. Spontaneous gating caused by the introduction of the β2(L285R) mutation was also reduced in α1(T47R)β2(L285R)γ2 compared to α1β2(L285R)γ2 receptors, confirming that α1(T47R) impairs gating independently of activation by any agonist. These findings implicate movement of the GABA_A receptor α1 subunit's β1 strand during agonist-dependent and spontaneous gating. Immobilisation of the β1 strand may provide a mechanism for the inhibition of gating by inverse agonists such as bicuculline.

Synthesis and Characterization of a Novel γ-Aminobutyric Acid Type A (GABAA) Receptor Ligand That Combines Outstanding Metabolic Stability, Pharmacokinetics, and Anxiolytic Efficacy

1,4-Benzodiazepines are used in the treatment of anxiety disorders but have limited long-term use due to adverse effects. HZ-166 (2) has been shown to have anxiolytic-like effects with reduced sedative/ataxic liabilities. A 1,3-oxazole KRM-II-81 (9) was discovered from a series of six bioisosteres with significantly improved pharmacokinetic and pharmacodynamic properties as compared to 2. Oxazole 9 was further characterized and exhibited improved anxiolytic-like effects in a mouse marble burying assay and a rat Vogel conflict test.

A role for loop G in the β1 strand in GABAA receptor activation

KEY POINTS

The role of the β1 strand in GABAA receptor function is unclear. It lies anti-parallel to the β2 strand, which is known to participate in receptor activation. Molecular dynamics simulation revealed solvent accessible residues within the β1 strand of the GABAA β3 homopentamer that might be amenable to analysis using the substituted Cys accessibility method. Cys substitutions from Asp43 to Thr47 in the GABAA α1 subunit showed that D43C and T47C reduced the apparent potency of GABA. F45C caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Cys43 and Cys47 were accessible to 2-aminoethyl methanethiosulphonate (MTSEA) modification, whereas Cys45 was not. Both GABA and the allosteric agonist propofol reduced MTSEA modification of Cys43 and Cys47. By contrast, modification of Cys64 in the β2 strand loop D was impeded by GABA but unaffected by propofol. These data reveal movement of β1 strand loop G residues during agonist activation of the GABAA receptor.

ABSTRACT

The GABAA receptor α subunit β1 strand runs anti-parallel to the β2 strand, which contains loop D, known to participate in receptor activation and agonist binding. However, a role for the β1 strand has yet to be established. We used molecular dynamics simulation to quantify the solvent accessible surface area (SASA) of β1 strand residues in the GABAA β3 homopentamer structure. Residues in the complementary interface equivalent to those between Asp43 and Thr47 in the α1 subunit have an alternating pattern of high and low SASA consistent with a β strand structure. We investigated the functional role of these β1 strand residues in the α1 subunit by individually replacing them with Cys residues. D43C and T47C substitutions reduced the apparent potency of GABA at α1β2γ2 receptors by 50-fold and eight-fold, respectively, whereas the F45C substitution caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Receptors with D43C or T47C substitutions were sensitive to 2-aminoethyl methanethiosulphonate (MTSEA) modification. However, GABA-evoked currents mediated by α1(F45C)β2γ2 receptors were unaffected by MTSEA, suggesting that this residue is inaccessible. Both GABA and the allosteric agonist propofol reduced MTSEA modification of α1(D43C)β2γ2 and α1(T47C)β2γ2 receptors, indicating movement of the β1 strand even during allosteric activation. This is in contrast to α1(F64C)β2γ2 receptors, where only GABA, but not propofol, reduced MTSEA modification. These findings provide the first functional evidence for movement of the β1 strand during gating of the receptor and identify residues that are critical for maintaining GABAA receptor function.

Assembly, trafficking and function of α1β2γ2 GABAA receptors are regulated by N-terminal regions, in a subunit-specific manner

GABAA receptors are pentameric ligand-gated ion channels that mediate inhibitory fast synaptic transmission in the central nervous system. Consistent with recent pentameric ligand-gated ion channels structures, sequence analysis predicts an α-helix near the N-terminus of each GABAA receptor subunit. Preceding each α-helix are 8-36 additional residues, which we term the N-terminal extension. In homomeric GABAC receptors and nicotinic acetylcholine receptors, the N-terminal α-helix is functionally essential. Here, we determined the role of the N-terminal extension and putative α-helix in heteromeric α1β2γ2 GABAA receptors. This role was most prominent in the α1 subunit, with deletion of the N-terminal extension or further deletion of the putative α-helix both dramatically reduced the number of functional receptors at the cell surface. Conversely, deletion of the β2 or γ2 N-terminal extension had little effect on the number of functional cell surface receptors. Additional deletion of the putative α-helix in the β2 or γ2 subunits did, however, decrease both functional cell surface receptors and incorporation of the γ2 subunit into mature receptors. In the β2 subunit only, α-helix deletions affected GABA sensitivity and desensitization. Our findings demonstrate that N-terminal extensions and α-helices make key subunit-specific contributions to assembly, consistent with both regions being involved in inter-subunit interactions. N-terminal α-helices and preceding sequences of eukaryotic pentameric ligand-gated ion channels are absent in prokaryotic homologues, suggesting they may not be functionally essential. Here, we show that in heteropentameric α1β2γ2 GABAA receptors, the role of these segments is highly subunit dependent. The extension preceding the α-helix in the α subunit is crucial for assembly and trafficking, but is of little importance in β and γ subunits. Indeed, robust receptor levels remain when the extension and α-helix are removed in β or γ subunits.

Molecular basis for convergent evolution of glutamate recognition by pentameric ligand-gated ion channels

Glutamate is an indispensable neurotransmitter, triggering postsynaptic signals upon recognition by postsynaptic receptors. We questioned the phylogenetic position and the molecular details of when and where glutamate recognition arose in the glutamate-gated chloride channels. Experiments revealed that glutamate recognition requires an arginine residue in the base of the binding site, which originated at least three distinct times according to phylogenetic analysis. Most remarkably, the arginine emerged on the principal face of the binding site in the Lophotrochozoan lineage, but 65 amino acids upstream, on the complementary face, in the Ecdysozoan lineage. This combined experimental and computational approach throws new light on the evolution of synaptic signalling.

Role of the ρ1 GABA(C) receptor N-terminus in assembly, trafficking and function

The GABAC receptor and closely related GABAA receptor are members of the pentameric ligand-gated ion channels (pLGICs) superfamily and mediate inhibitory fast synaptic transmission in the nervous system. Each pLGIC subunit comprises an N-terminal extracellular agonist-binding domain followed by a channel domain and a variable intracellular domain. Available structural information shows that the core of the agonist-binding domain is a β sandwich of ten β-strands, which form the agonist-binding pocket at the subunit interface. This β-sandwich is preceded by an N-terminal α-helix in eukaryotic structures but not in prokaryotic structures. The N-terminal α-helix has been shown to be functionally essential in α7 nicotinic acetylcholine receptors. Sequence analysis of GABAC and GABAA receptors predicts an α-helix in a similar position but preceded by 8 to 46 additional residues, of unknown function, which we term the N-terminal extension. To test the functional role of both the N-terminal extension and the putative N-terminal α-helix in the ρ1 GABAC receptor, we created a series of deletions from the N-terminus. The N-terminal extension was not functionally essential, but its removal did reduce both cell surface expression and cooperativity of agonist-gated channel function. Further deletion of the putative N-terminal α-helix abolished receptor function by preventing cell-surface expression. Our results further demonstrate the essential role of the N-terminal α-helix in the assembly and trafficking of eukaryotic pLGICs. They also provide evidence that the N-terminal extension, although not essential, contributes to receptor assembly, trafficking and conformational changes associated with ligand gating.

Allosteric modulation of GABAA receptors via multiple drug-binding sites

GABAA receptors are ligand-gated ion channels composed of five subunits that can be opened by GABA and be modulated by multiple pharmacologically and clinically important drugs. Over the time, hundreds of compounds from different structural classes have been demonstrated to modulate, directly activate, or inhibit GABAA receptors, and most of these compounds interact with more than one binding site at these receptors. Crystal structures of proteins and receptors homologous to GABAA receptors as well as homology modeling studies have provided insights into the possible location of ligand interaction sites. Some of these sites have been identified by mutagenesis, photolabeling, and docking studies. For most of these ligands, however, binding sites are not known. Due to the high flexibility of GABAA receptors and the existence of multiple drug-binding sites, the unequivocal identification of interaction sites for individual drugs is extremely difficult. The existence of multiple GABAA receptor subtypes with distinct subunit composition, the contribution of distinct subunit sequences to binding sites of different receptor subtypes, as well as the observation that even subunits not directly contributing to a binding site are able to influence affinity and efficacy of drugs, contribute to a unique pharmacology of each GABAA receptor subtype. Thus, each receptor subtype has to be investigated to identify a possible subtype selectivity of a compound. Although multiple binding sites make GABAA receptor pharmacology even more complicated, the exploitation of ligand interaction with novel-binding sites also offers additional possibilities for a subtype-selective modulation of GABAA receptors.

Principles of agonist recognition in Cys-loop receptors

Cys-loop receptors are ligand-gated ion channels that are activated by a structurally diverse array of neurotransmitters, including acetylcholine, serotonin, glycine, and GABA. After the term "chemoreceptor" emerged over 100 years ago, there was some wait until affinity labeling, molecular cloning, functional studies, and X-ray crystallography experiments identified the extracellular interface of adjacent subunits as the principal site of agonist binding. The question of how subtle differences at and around agonist-binding sites of different Cys-loop receptors can accommodate transmitters as chemically diverse as glycine and serotonin has been subject to intense research over the last three decades. This review outlines the functional diversity and current structural understanding of agonist-binding sites, including those of invertebrate Cys-loop receptors. Together, this provides a framework to understand the atomic determinants involved in how these valuable therapeutic targets recognize and bind their ligands.

α1F64 Residue at GABA(A) receptor binding site is involved in gating by influencing the receptor flipping transitions

GABA receptors (GABAARs) mediate inhibition in the adult brain. These channels are heteropentamers and their ligand binding sites are localized at the β+ / α- interfaces. As expected, mutations of binding-site residues affect binding kinetics but accumulating evidence indicates that gating is also altered, although the underlying mechanisms are unclear. We investigated the impact of the hydrophobic box residue localized at α1(-), F64 (α1F64), on the binding and gating of rat recombinant α1β1γ2 receptors. The analysis of current responses to rapid agonist applications confirmed a marked effect of α1F64 mutations on agonist binding and revealed surprisingly strong effects on gating, including the disappearance of rapid desensitization, the slowing of current onset, and accelerated deactivation. Moreover, nonstationary variance analysis revealed that the α1F64C mutation dramatically reduced the maximum open probability without altering channel conductance. Interestingly, for wild-type receptors, responses to saturating concentration of a partial agonist, P4S, showed no rapid desensitization, similar to GABA-evoked responses mediated by α1F64C mutants. For the α1F64L mutation, the application of the high-affinity agonist muscimol partially rescued rapid desensitization compared with responses evoked by GABA. These findings suggest that α1F64 mutations do not disrupt desensitization mechanisms but rather affect other gating features that obscure it. Model simulations indicated that all of our observations related to α1F64 mutations could be properly reproduced by altering the flipped state transitions that occurred after agonist binding but preceded opening. In conclusion, we propose that the α1F64 residue may participate in linking binding and gating by influencing flipping kinetics.

Gating of pentameric ligand-gated ion channels: structural insights and ambiguities

Pentameric ligand-gated ion channels (pLGICs) mediate fast synaptic communication by converting chemical signals into an electrical response. Recently solved agonist-bound and agonist-free structures greatly extend our understanding of these complex molecular machines. A key challenge to a full description of function, however, is the ability to unequivocally relate determined structures to the canonical resting, open, and desensitized states. Here, we review current understanding of pLGIC structure, with a focus on the conformational changes underlying channel gating. We compare available structural information and review the evidence supporting the assignment of each structure to a particular conformational state. We discuss multiple factors that may complicate the interpretation of crystal structures, highlighting the potential influence of lipids and detergents. We contend that further advances in the structural biology of pLGICs will require deeper insight into the nature of pLGIC-lipid interactions.

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.

Probing the orthosteric binding site of GABAA receptors with heterocyclic GABA carboxylic acid bioisosteres

The ionotropic GABAA receptors (GABAARs) are widely distributed in the central nervous system where they play essential roles in numerous physiological and pathological processes. A high degree of structural heterogeneity of the GABAAR has been revealed and extensive effort has been made to develop selective and potent GABAAR agonists. This review investigates the use of heterocyclic carboxylic acid bioisosteres within the GABAAR area. Several heterocycles including 3-hydroxyisoxazole, 3-hydroxyisoxazoline, 3-hydroxyisothiazole, and the 1- and 3-hydroxypyrazole rings have been employed in order to map the orthosteric binding site. The physicochemical properties of the heterocyclic moieties making them suitable for bioisosteric replacement of the carboxylic acid in the molecule of GABA are discussed. A variety of synthetic strategies for synthesis of the heterocyclic scaffolds are available. Likewise, methods for introduction of substituents into specific positions of the heterocyclic scaffolds facilitate the investigation of different regions in the orthosteric binding pocket in close vicinity of the core scaffolds of muscimol/GABA. The development of structural models, from pharmacophore models to receptor homology models, has provided more insight into the molecular basis for binding. Similar binding modes are proposed for the heterocyclic GABA analogues covered in this review by use of ligand-receptor docking into the receptor homology model and the presented structure-activity relationships. A network of interactions between the analogues and the binding pocket is leaving no room for substituents and underline the limited space in the GABAAR orthosteric binding site when in the agonist conformation.

Molecular cloning and characterization of novel glutamate-gated chloride channel subunits from Schistosoma mansoni

Cys-loop ligand-gated ion channels (LGICs) mediate fast ionotropic neurotransmission. They are proven drug targets in nematodes and arthropods, but are poorly characterized in flatworms. In this study, we characterized the anion-selective, non-acetylcholine-gated Cys-loop LGICs from Schistosoma mansoni. Full-length cDNAs were obtained for SmGluCl-1 (Smp_096480), SmGluCl-2 (Smp_015630) and SmGluCl-3 (Smp_104890). A partial cDNA was retrieved for SmGluCl-4 (Smp_099500/Smp_176730). Phylogenetic analyses suggest that SmGluCl-1, SmGluCl-2, SmGluCl-3 and SmGluCl-4 belong to a novel clade of flatworm glutamate-gated chloride channels (GluCl) that includes putative genes from trematodes and cestodes. The flatworm GluCl clade was distinct from the nematode-arthropod and mollusc GluCl clades, and from all GABA receptors. We found no evidence of GABA receptors in S. mansoni. SmGluCl-1, SmGluCl-2 and SmGluCl-3 subunits were characterized by two-electrode voltage clamp (TEVC) in Xenopus oocytes, and shown to encode Cl⁻-permeable channels gated by glutamate. SmGluCl-2 and SmGluCl-3 produced functional homomers, while SmGluCl-1 formed heteromers with SmGluCl-2. Concentration-response relationships revealed that the sensitivity of SmGluCl receptors to L-glutamate is among the highest reported for GluCl receptors, with EC₅₀ values of 7–26 µM. Chloride selectivity was confirmed by current-voltage (I/V) relationships. SmGluCl receptors are insensitive to 1 µM ivermectin (IVM), indicating that they do not belong to the highly IVM-sensitive GluClα subtype group. SmGluCl receptors are also insensitive to 10 µM meclonazepam, a schistosomicidal benzodiazepine. These results provide the first molecular evidence showing the contribution of GluCl receptors to L-glutamate signaling in S. mansoni, an unprecedented finding in parasitic flatworms. Further work is needed to elucidate the roles of GluCl receptors in schistosomes and to explore their potential as drug targets.

Schistosomiasis is a debilitating disease caused by blood flukes in the genus Schistosoma that afflicts over 200 million people worldwide. Treatment relies almost exclusively on a single drug, praziquantel. Reports of sub-optimal efficacy of praziquantel raise concerns about the prospect of drug resistance and highlight the need to develop new schistosomicidal drugs. Neuroactive receptors are recognized targets of insecticides and anthelmintics. Likewise, neuronal receptors of schistosomes are attractive targets for drug development. Lacking a coelom and a proper circulatory system, schistosomes are thought to lack the capacity for endocrine signaling, and therefore depend entirely on neuronal modulation to control functions vital to their survival and reproduction. We characterized a novel family of glutamate-gated chloride channel (GluCl) receptors from S. mansoni that are pharmacologically and evolutionarily distinct from GluCls in nematodes, insects and snails. Our phylogenetic analyses suggest that these receptors are also widely distributed in other flukes and tapeworms. This study provides the first molecular evidence for the contribution of an inhibitory component to glutamatergic signaling in S. mansoni. Our findings add to a growing body of evidence suggesting that glutamatergic signaling in schistosomes may be physiologically important, and could be targeted for chemotherapeutic intervention.

2-Guanidine-4-methylquinazoline acts as a novel competitive antagonist of A type γ-aminobutyric acid receptors

The pentameric A type γ-aminobutyric acid receptors (GABAARs) are the major inhibitory neurotransmitter receptors in the nervous system and have long been considered as important pharmaceutical targets for the treatment of multiple neurological or psychological disorders. Here, we show that 2-guanidine-4-methylquinazoline (GMQ), a recently identified acid-sensing ion channel (ASIC) modulator, strongly and preferentially inhibits GABAAR among the major neurotransmitter-gated ion channels in cultured rat hippocampal neurons. GMQ inhibited GABA (1 μM)-induced currents in a competitive manner, with an IC50 (0.39±0.05 μM) comparable to that of bicuculline. Schild analysis revealed a slope of 1.04±0.06 for GMQ on α1β2 GABAARs expressed in HEK293T cells. Single-channel analysis showed that GMQ decreased open probability of GABAARs without affecting conductance. Moreover, GMQ inhibited GABAergic neurotransmission in hippocampal neurons, while having no significant effect on the basal field excitatory postsynaptic potentials (fEPSPs) and the intrinsic excitability of neurons. Using site-directed mutagenesis, we further demonstrated that mutations at Glu155 of β2 subunit and Phe64 of α1 subunit, both located inside the GABA binding pocket, profoundly decreased the sensitivity of the receptor to both GABA and GMQ. Interestingly, these mutations did not significantly affect the inhibition by amiloride, a diuretic structurally similar to GMQ and a known GABAAR inhibitor. We conclude that GMQ represents a novel chemical structure that acts, possibly, by competing with GABA binding to GABAARs. It is anticipated that GMQ and its analogs will facilitate the development of new chemical probes for GABAARs.

SB-205384 is a positive allosteric modulator of recombinant GABAA receptors containing rat α3, α5, or α6 subunit subtypes coexpressed with β3 and γ2 subunits

Many drugs used to treat anxiety are positive modulators of GABAA receptors, which mediate fast inhibitory neurotransmission. The GABAA receptors can be assembled from a combination of at least 16 different subunits. The receptor's subunit composition determines its pharmacologic and functional properties, and subunit expression varies throughout the brain. A primary goal for new treatments targeting GABAA receptors is the production of subunit-selective modulators acting upon a discrete population of receptors. The anxiolytic 4-amino-7-hydroxy-2-methyl-5,6,7,8,-tetrahydrobenzo[b]thieno[2,3-b]pyridine-3-carboxylic acid, but-2-ynyl ester (SB-205384) is widely considered to be selective for α3-containing GABAA receptors. However, it has been tested only on α1-, α2-, and α3-containing receptors. We examined the activity of SB-205384 at recombinant receptors containing the six different α subunits and found that receptors containing the α3, α5, and α6 subunits were potentiated by SB-205384, with the α6 subunit conferring the greatest responsiveness. Properties associated with chimeric α1/α6 subunits suggested that multiple structural domains influence sensitivity to SB-205384. Point mutations of residues within the extracellular N-terminal domain identified a leucine residue located in loop E of the agonist binding site as an important determinant of high sensitivity to modulation. In the α6 subunit the identity of this residue is species-dependent, with the leucine found in rat subunits but not in human. Our results indicate that SB-205384 is not an α3-selective modulator, and instead acts at several GABAA receptor isoforms. These findings have implications for the side-effect profile of this anxiolytic as well as for its use in neuronal and animal studies as a marker for contribution from α3-containing receptors.

Physostigmine and galanthamine bind in the presence of agonist at the canonical and noncanonical subunit interfaces of a nicotinic acetylcholine receptor

Galanthamine and physostigmine are clinically used cholinomimetics that both inhibit acetylcholinesterase and also interact directly with and potentiate nAChRs. As with most nAChR-positive allosteric modulators, the location and number of their binding site(s) within nAChRs are unknown. In this study, we use the intrinsic photoreactivities of [(3)H]physostigmine and [(3)H]galanthamine upon irradiation at 312 nm to directly identify amino acids contributing to their binding sites in the Torpedo californica nAChR. Protein sequencing of fragments isolated from proteolytic digests of [(3)H]physostigmine- or [(3)H]galanthamine-photolabeled nAChR establish that, in the presence of agonist (carbamylcholine), both drugs photolabeled amino acids on the complementary (non-α) surface of the transmitter binding sites (γTyr-111/γTyr-117/δTyr172). They also photolabeled δTyr-212 at the δ-β subunit interface and γTyr-105 in the vestibule of the ion channel, with photolabeling of both residues enhanced in the presence of agonist. Furthermore, [(3)H]physostigmine photolabeling of γTyr-111, γTyr-117, δTyr-212, and γTyr-105 was inhibited in the presence of nonradioactive galanthamine. The locations of the photolabeled amino acids in the nAChR structure and the results of computational docking studies provide evidence that, in the presence of agonist, physostigmine and galanthamine bind to at least three distinct sites in the nAChR extracellular domain: at the α-γ interface (1) in the entry to the transmitter binding site and (2) in the vestibule of the ion channel near the level of the transmitter binding site, and at the δ-β interface (3) in a location equivalent to the benzodiazepine binding site in GABA(A) receptors.

A unified model of the GABA(A) receptor comprising agonist and benzodiazepine binding sites

We present a full-length α(1)β(2)γ(2) GABA receptor model optimized for agonists and benzodiazepine (BZD) allosteric modulators. We propose binding hypotheses for the agonists GABA, muscimol and THIP and for the allosteric modulator diazepam (DZP). The receptor model is primarily based on the glutamate-gated chloride channel (GluCl) from C. elegans and includes additional structural information from the prokaryotic ligand-gated ion channel ELIC in a few regions. Available mutational data of the binding sites are well explained by the model and the proposed ligand binding poses. We suggest a GABA binding mode similar to the binding mode of glutamate in the GluCl X-ray structure. Key interactions are predicted with residues α(1)R66, β(2)T202, α(1)T129, β(2)E155, β(2)Y205 and the backbone of β(2)S156. Muscimol is predicted to bind similarly, however, with minor differences rationalized with quantum mechanical energy calculations. Muscimol key interactions are predicted to be α(1)R66, β(2)T202, α(1)T129, β(2)E155, β(2)Y205 and β(2)F200. Furthermore, we argue that a water molecule could mediate further interactions between muscimol and the backbone of β(2)S156 and β(2)Y157. DZP is predicted to bind with interactions comparable to those of the agonists in the orthosteric site. The carbonyl group of DZP is predicted to interact with two threonines α(1)T206 and γ(2)T142, similar to the acidic moiety of GABA. The chlorine atom of DZP is placed near the important α(1)H101 and the N-methyl group near α(1)Y159, α(1)T206, and α(1)Y209. We present a binding mode of DZP in which the pending phenyl moiety of DZP is buried in the binding pocket and thus shielded from solvent exposure. Our full length GABA(A) receptor is made available as Model S1.

Ivermectin binding sites in human and invertebrate Cys-loop receptors

Ivermectin is a gold standard antiparasitic drug that has been used successfully to treat billions of humans, livestock and pets. Until recently, the binding site on its Cys-loop receptor target had been a mystery. Recent protein crystal structures, site-directed mutagenesis data and molecular modelling now explain how ivermectin binds to these receptors and reveal why it is selective for invertebrate members of the Cys-loop receptor family. Combining this with emerging genomic information, we are now in a position to predict species sensitivity to ivermectin and better understand the molecular basis of ivermectin resistance. An understanding of the molecular structure of the ivermectin binding site, which is formed at the interface of two adjacent subunits in the transmembrane domain of the receptor, should also aid the development of new lead compounds both as anthelmintics and as therapies for a wide variety of human neurological disorders.

Teaching an old GABA receptor new tricks

The accompanying articles in this issue of the journal's special collection describe attempts to improve on the dynamics of distribution and reduce side effects of analogs of etomidate and benzodiazepines. Both classes of drugs have their principal sites of action on γ-aminobutyric acid type A receptors, although at very different binding sites and by different mechanisms of action. Herein, we review the structure of γ-aminobutyric acid type A receptors and describe the location of the 2 likely binding sites. In addition, we describe how these drugs can interact with the nervous system at a systems level. We leave it to other reviewers to discuss whether these new drugs offer true clinical improvements.

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.

New insights into the structural bases of activation of Cys-loop receptors

Neurotransmitter receptors of the Cys-loop superfamily mediate rapid synaptic transmission throughout the nervous system, and include receptors activated by ACh, GABA, glycine and serotonin. They are involved in physiological processes, including learning and memory, and in neurological disorders, and they are targets for clinically relevant drugs. Cys-loop receptors assemble either from five copies of one type of subunit, giving rise to homomeric receptors, or from several types of subunits, giving rise to heteromeric receptors. Homomeric receptors are invaluable models for probing fundamental relationships between structure and function. Receptors contain a large extracellular domain that carries the binding sites and a transmembrane region that forms the ion pore. How the structural changes elicited by agonist binding are propagated through a distance of 50Å to the ion channel gate is central to understanding receptor function. Depending on the receptor subtype, occupancy of either two, as in the prototype muscle nicotinic receptor, or three binding sites, as in homomeric receptors, is required for full activation. The conformational changes initiated at the binding sites are propagated to the gate through the interface between the extracellular and transmembrane domains. This region forms a network that relays structural changes from the binding site towards the pore, and also contributes to open channel lifetime and rate of desensitization. Thus, this coupling region controls the beginning and duration of a synaptic response. Here we review recent advances in the molecular mechanism by which Cys-loop receptors are activated with particular emphasis on homomeric receptors.

Three arginines in the GABAA receptor binding pocket have distinct roles in the formation and stability of agonist- versus antagonist-bound complexes

Binding of the agonist GABA to the GABA(A) receptor causes channel gating, whereas competitive antagonists that bind at the same site do not. The details of ligand binding are not well understood, including which residues interact directly with ligands, maintain the structure of the binding pocket, or transduce the action of binding into opening of the ion channel gate. Recent work suggests that the amine group of the GABA molecule may form a cation-π bond with residues in a highly conserved "aromatic box" within the binding pocket. Although interactions with the carboxyl group of GABA remain unknown, three positively charged arginines (α(1)Arg67, α(1)Arg132, and β(2)Arg207) just outside of the aromatic box are likely candidates. To explore their roles in ligand binding, we individually mutated these arginines to alanine and measured the effects on microscopic ligand binding/unbinding rates and channel gating. The mutations α(1)R67A or β(2)R207A slowed agonist binding and sped unbinding with little effect on gating, demonstrating that these arginines are critical for both formation and stability of the agonist-bound complex. In addition, α(1)R67A sped binding of the antagonist 2-(3-carboxypropyl)-3-amino-6-(4 methoxyphenyl)pyridazinium bromide (SR-95531), indicating that this arginine poses a barrier to formation of the antagonist-bound complex. In contrast, β(2)R207A and α(1)R132A sped antagonist unbinding, indicating that these arginines stabilize the antagonist-bound state. α(1)R132A also conferred a new long-lived open state, indicating that this arginine influences the channel gate. Thus, each of these arginines plays a unique role in determining interactions with agonists versus antagonists and with the channel gate.

The GABRA6 mutation, R46W, associated with childhood absence epilepsy, alters 6β22 and 6β2 GABA(A) receptor channel gating and expression

A GABA(A) receptor α6 subunit mutation, R46W, was identified as a susceptibility gene that may contribute to the pathogenesis of childhood absence epilepsy (CAE), but the molecular basis for alteration of GABA(A) receptor function is unclear. The R46W mutation is located in a region homologous to a GABA(A) receptor γ2 subunit missense mutation, R82Q, that is associated with CAE and febrile seizures in humans. To determine how this mutation reduces GABAergic inhibition, we expressed wild-type (α6β2γ2L and α6β2δ) and mutant (α6(R46W)β2γ2L and α6(R46W)β2δ) receptors in HEK 293T cells and characterize their whole-cell and single-channel currents, and surface and total levels. We demonstrated that gating and assembly of both α6(R46W)β2γ2L and α6(R46W)β2δ receptors were impaired. Compared to wild-type currents, α6(R46W)β2γ2L and α6(R46W)β2δ receptors had a reduced current density, α6(R46W)β2γ2L currents desensitized to a greater extent and deactivated at a slower rate, α6(R46W)β2δ receptors did not desensitize but deactivated faster and both α6(R46W)β2γ2L and α6(R46W)β2δ single-channel current mean open times and burst durations were reduced. Surface levels of coexpressed α6(R46W), β2 and δ, but not γ2L, subunits were decreased. 'Heterozygous' coexpression of α6(R46W) and α6 subunits with β2 and γ2L subunits produced intermediate macroscopic current amplitudes by increasing incorporation of wild-type and decreasing incorporation of mutant subunits into receptors trafficked to the surface. Finally, these findings suggest that similar to the γ2(R82Q) mutation, the CAE-associated α6(R46W) mutation could cause neuronal disinhibition and thus increase susceptibility to generalized seizures through a reduction of αβγ and αβδ receptor function and expression.

Discrete M3-M4 intracellular loop subdomains control specific aspects of γ-aminobutyric acid type A receptor function

The GABA type A receptor (GABA(A)R) is a member of the pentameric ligand gated ion channel (pLGIC) family that mediates ionotropic neurotransmission. Residues in the intracellular loop domain (ILD) have recently been shown to define part of the ion permeation pathway in several closely related members of the pentameric ligand gated ion channel family. In this study, we investigated the role the ILD of the GABA(A)R α1 subunit plays in channel function. Deletion of the α1 ILD resulted in a significant increase in GABA EC(50) and maximal current amplitude, suggesting that the ILD must be intact for proper receptor function. To test this hypothesis, we conducted a mutagenic screen of all amino acids harboring ionizable side chains within this domain to investigate the contribution of individual charged residues to ion permeation. Using macroscopic and single channel voltage-clamp recording techniques, we found that mutations within a subdomain of the α1 ILD near M3 altered GABA apparent affinity; interestingly, α1(K312E) exhibited reduced partial agonist efficacy. We introduced point mutations near M4, including α1(K383E) and α1(K384E), that enhanced receptor desensitization. Mutation of 5 charged residues within a 39-residue span contiguous with M4 reduced relative anion permeability of the channel and may represent a weak intracellular selectivity filter. Within this subdomain, the α1(K378E) mutation induced a significant reduction in single channel conductance, consistent with our hypothesis that the GABA(A)R α1 ILD contributes directly to the permeation pathway.

A tight coupling between β₂Y97 and β₂F200 of the GABA(A) receptor mediates GABA binding

The GABA(A) receptor is an oligopentameric chloride channel that is activated via conformation changes induced upon the binding of the endogenous ligand, GABA, to the extracellular inter-subunit interfaces. Although dozens of amino acid residues at the α/β interface have been implicated in ligand binding, the structural elements that mediate ligand binding and receptor activation are not yet fully described. In this study, double-mutant cycle analysis was employed to test for possible interactions between several arginines (α₁R67, α₁R120, α₁R132, and β₂R207) and two aromatic residues (β₂Y97 and β₂F200) that are present in the ligand-binding pocket and are known to influence GABA affinity. Our results show that neither α₁R67 nor α₁R120 is functionally coupled to either of the aromatics, whereas a moderate coupling exists between α₁R132 and both aromatic residues. Significant functional coupling between β₂R207 and both β₂Y97 and β₂F200 was found. Furthermore, we identified an even stronger coupling between the two aromatics, β₂Y97 and β₂F200, and for the first time provided direct evidence for the involvement of β₂Y97 and β₂F200 in GABA binding. As these residues are tightly linked, and mutation of either has similar, severe effects on GABA binding and receptor kinetics, we believe they form a single functional unit that may directly coordinate GABA.

Synganglion transcriptome and developmental global gene expression in adult females of the American dog tick, Dermacentor variabilis (Acari: Ixodidae)

454 Pyrosequencing was used to characterize the expressed genes from the synganglion and associated neurosecretory organs of unfed and partially fed virgin and mated replete females of the American dog tick, Dermacentor variabilis. A total of 14,881 contiguous sequences (contigs) was assembled, with an average size of 229 bp. Gene ontology terms for Level 2 biological processes were assigned to 4366 contigs. Seven acetylcholinesterases, a muscarinic acetylcholine (ACh) receptor, two nicotinic ACh receptor β-subunits, two ACh unc-18 regulators, two dopamine receptors, two gamma aminobutyric acid (GABA) receptors, two GABA transporters, two norepinephrine transporters and an octopamine receptor are described. Microarrays were conducted to examine global gene expression and quantitative real-time polymerase chain reaction was used to verify expression of selected neuropeptides. Hierarchical clustering of all differentially expressed transcripts grouped part-fed and replete ticks as being more similar in terms of differentially expressed genes with unfed ticks as the outgroup. Nine putative neuropeptides (allatostatin, bursicon-β, preprocorazonin, glycoprotein hormone α, insulin-like peptide, three orcokinins, preprosulphakinin) and a gonadotropin releasing hormone receptor were differentially expressed, and their developmental expression and role in reproduction was investigated. The presence of eclosion hormone, corazonin and bursicon in the synganglion, which in insects regulate behaviour and cuticle development associated with moulting, suggest that this system may be used in ticks to regulate blood feeding, cuticle expansion and development related to female reproduction; adult ticks do not moult.

New insight into the central benzodiazepine receptor-ligand interactions: design, synthesis, biological evaluation, and molecular modeling of 3-substituted 6-phenyl-4H-imidazo[1,5-a][1,4]benzodiazepines and related compounds

3-Substituted 6-phenyl-4H-imidazo[1,5-a][1,4]benzodiazepines and related compounds were synthesized as central benzodiazepine receptor (CBR) ligands. Most of the compounds showed high affinity for bovine and human CBR, their K(i) values spanning from the low nanomolar to the submicromolar range. In particular, imidazoester 5f was able to promote a massive flow of (36)Cl(-) in rat cerebrocortical synaptoneurosomes overlapping its efficacy profile with that of a typical full agonist. Compound 5f was then examined in mice for its pharmacological effects where it proved to be a safe anxiolytic agent devoid of the unpleasant myorelaxant and amnesic effects of the classical 1,4-benzodiazepines. Moreover, the selectivity of some selected compounds has been assessed in recombinant α(1)β(2)γ(2)L, α(2)β(1)γ(2)L, and α(5)β(2)γ(2)L human GABA(A) receptors. Finally, some compounds were submitted to molecular docking calculations along with molecular dynamics simulations in the Cromer's GABA(A) homology model.

Ligand activation of the prokaryotic pentameric ligand-gated ion channel ELIC

While the pentameric ligand-gated ion channel ELIC has recently provided first insight into the architecture of the family at high resolution, its detailed investigation was so far prevented by the fact that activating ligands were unknown. Here we describe a study on the functional characterization of ELIC by electrophysiology and X-ray crystallography. ELIC is activated by a class of primary amines that include the neurotransmitter GABA at high micro- to millimolar concentrations. The ligands bind to a conserved site and evoke currents that slowly desensitize over time. The protein forms cation selective channels with properties that resemble the nicotinic acetylcholine receptor. The high single channel conductance and the comparably simple functional behavior make ELIC an attractive model system to study general mechanisms of ion conduction and gating in this important family of neurotransmitter receptors.

The Haemonchus contortus UNC-49B subunit possesses the residues required for GABA sensitivity in homomeric and heteromeric channels

Hco-UNC-49 is a GABA receptor from the parasitic nematode Haemonchus contortus that has a relatively low overall sequence similarity to vertebrate GABA receptors but is very similar to the UNC-49 receptor found in the free living nematode Caenorhabditis elegans. While the nematode receptors do share >80% sequence similarity they exhibit different sensitivities to GABA. In addition, the UNC-49C subunit appears to be a positive modulator of GABA sensitivity in the H. contortus heteromeric channel, but is a negative modulator in the C. elegans heteromeric channel. The cause(s) of these differences is currently unknown since the structural elements essential for GABA sensitivity in nematode receptors have been largely unexplored. Thus, the overall aim of this study was to investigate the residues that are important for UNC-49 receptor sensitivity through the use of homology modeling, site-directed mutagenesis, and two-electrode voltage clamp. This study revealed that Met(170) in Loop B of the GABA binding-site may partially account for the observed differences in GABA receptor sensitivity between the nematode species. Residues in Loops A-D that have been reported to form the GABA binding pocket in mammalian receptors, including those forming the conserved 'aromatic box', also appear to play analogous roles in Hco-UNC-49. In addition, the two mutations that produced the most significant reduction in GABA sensitivity were R66S and Y166S. Homology modeling indicates that these two residues share a hydrogen bond and are positioned close to the carboxyl end of the GABA molecule. However, of residues examined in this study, only those on the Hco-UNC-49B subunit and not its subunit partner, Hco-UNC-49C, appear important for GABA sensitivity. Overall, results from this study suggest that the binding site of the UNC-49 heteromeric GABA receptor exhibits some differences compared to classical vertebrate GABA(A) receptors.

A state-dependent salt-bridge interaction exists across the β/α intersubunit interface of the GABAA receptor

The GABA(A) receptor is a multisubunit protein that transduces the binding of a neurotransmitter at an intersubunit interface into the opening of a central ion channel. The structural components that mediate the steps involved in this action are poorly defined. A large amount of work has focused on clarifying the specific functions and interactions of residues believed to surround the GABA binding pocket. Here, we explored two charged residues (β(2)Asp163 and α(1)Arg120), which have been suggested by homology models to participate in a salt-bridge interaction. When mutated to alanine, both single mutants, as well as the double mutant, increase EC(50-GABA), decrease the GABA binding rate, and accelerate deactivation and GABA unbinding rates. Double-mutant cycle analysis demonstrates that the effects of each alanine mutation on the GABA binding rate were additive and independent. In contrast, a significant coupling energy was found during an analysis of deactivation time constants. Using kinetic modeling, we further demonstrated that the GABA unbinding rates, in particular, are strongly coupled. These data suggest that β(2)Asp163 and α(1)Arg120 form a state-dependent salt bridge, interacting when GABA is bound to the receptor but not when the receptor is in the unbound state.

Influence of GABAAR monoliganded states on GABAergic responses

To reach the open state, the GABA(A) receptor (GABA(A)R) is assumed to bind two agonist molecules. Although it is currently believed that GABA(A)R could also operate in the monoliganded state, the gating properties of singly bound GABA(A)R are poorly understood and their physiological role is still obscure. In the present study, we characterize for the first time the gating properties of singly bound GABA(A)Rs by using a mutagenesis approach and we propose that monoliganded GABA(A)R contribute in shaping synaptic responses. At saturating GABA concentrations, currents mediated by recombinant GABA(A)Rs with a single functional binding site display slow onset, fast deactivation kinetics, and slow rate of desensitization-resensitization. GABA(A)Rs with two binding sites activated by brief pulses of subsaturating GABA concentrations (in the range of the GABA concentration profile in the synaptic cleft) could also mediate fast deactivating currents, displaying deactivation kinetics similar to those mediated by GABA(A)Rs with a single functional binding site. Model simulations of receptors activated by realistic synaptic GABA waves revealed that a considerable proportion of GABA(A) receptors open in the monoliganded state during synaptic transmission, therefore contributing in shaping IPSCs.

New insights into the GABA(A) receptor structure and orthosteric ligand binding: receptor modeling guided by experimental data

GABA(A) receptors (GABA(A)Rs) are ligand gated chloride ion channels that mediate overall inhibitory signaling in the CNS. A detailed understanding of their structure is important to gain insights in, e.g., ligand binding and functional properties of this pharmaceutically important target. Homology modeling is a necessary tool in this regard because experimentally determined structures are lacking. Here we present an exhaustive approach for creating a high quality model of the α(1)β(2)γ(2) subtype of the GABA(A)R ligand binding domain, and we demonstrate its usefulness in understanding details of orthosteric ligand binding. The model was constructed by using multiple templates and by incorporation of knowledge from biochemical/pharmacological experiments. It was validated on the basis of objective energy functions, its ability to account for available residue specific information, and its stability in molecular dynamics (MD) compared with that of the two homologous crystal structures. We then combined the model with extensive structure-activity relationships available from two homologous series of orthosteric GABA(A)R antagonists to create a detailed hypothesis for their binding modes. Excellent agreement with key experimental data was found, including the ability of the model to accommodate and explain a previously developed pharmacophore model. A coupling to agonist binding was thereby established and discussed in relation to activation mechanisms. Our results highlight the importance of critical evaluation and optimization of each step in the homology modeling process. The approach taken here can greatly aid in increasing the understanding of GABA(A)Rs and related receptors where structural insight is limited and reliable models are difficult to obtain.

Design, synthesis, and subtype selectivity of 3,6-disubstituted β-carbolines at Bz/GABA(A)ergic receptors. SAR and studies directed toward agents for treatment of alcohol abuse

A series of 3,6-disubstituted β-carbolines was synthesized and evaluated for their in vitro affinities at α(x)β(3)γ(2) GABA(A)/benzodiazepine receptor subtypes by radioligand binding assays in search of α(1) subtype selective ligands to treat alcohol abuse. Analogues of β-carboline-3-carboxylate-t-butyl ester (βCCt, 1) were synthesized via a CDI-mediated process and the related 6-substituted β-carboline-3-carboxylates 6 including WYS8 (7) were synthesized via a Sonogashira or Stille coupling processes from 6-iodo-βCCt (5). The bivalent ligands of βCCt (32 and 33) were also designed and prepared via a palladium-catalyzed homocoupling process to expand the structure-activity relationships (SAR) to larger ligands. Based on the pharmacophore/receptor model, a preliminary SAR study on 34 analogues illustrated that large substituents at position-6 of the β-carbolines were well tolerated. As expected, these groups are proposed to project into the extracellular domain (L(Di) region) of GABA(A)/Bz receptors (see 32 and 33). Moreover, substituents located at position-3 of the β-carboline nucleus exhibited a conserved stereo interaction in lipophilic pocket L(1), while N(2) presumably underwent a hydrogen bonding interaction with H(1). Three novel β-carboline ligands (βCCt, 3PBC and WYS8), which preferentially bound to α1 BzR subtypes permitted a comparison of the pharmacological efficacies with a range of classical BzR antagonists (flumazenil, ZK93426) from several different structural groups and indicated these β-carbolines were 'near GABA neutral antagonists'. Based on the SAR, the most potent (in vitro) α(1) selective ligand was the 6-substituted acetylenyl βCCt (WYS8, 7). Earlier both βCCt and 3PBC had been shown to reduce alcohol self-administration in alcohol preferring (P) and high alcohol drinking (HAD) rats but had little or no effect on sucrose self-administration.(1-3) Moreover, these two β-carbolines were orally active, and in addition, were anxiolytic in P rats but were only weakly anxiolytic in rodents. These data prompted the synthesis of the β-carbolines presented here.

Expression of GABAergic receptors in mouse taste receptor cells

Background

Multiple excitatory neurotransmitters have been identified in the mammalian taste transduction, with few studies focused on inhibitory neurotransmitters. Since the synthetic enzyme glutamate decarboxylase (GAD) for gamma-aminobutyric acid (GABA) is expressed in a subset of mouse taste cells, we hypothesized that other components of the GABA signaling pathway are likely expressed in this system. GABA signaling is initiated by the activation of either ionotropic receptors (GABA_A and GABA_C) or metabotropic receptors (GABA_B) while it is terminated by the re-uptake of GABA through transporters (GATs).

Methodology/Principal Findings

Using reverse transcriptase-PCR (RT-PCR) analysis, we investigated the expression of different GABA signaling molecules in the mouse taste system. Taste receptor cells (TRCs) in the circumvallate papillae express multiple subunits of the GABA_A and GABA_B receptors as well as multiple GATs. Immunocytochemical analyses examined the distribution of the GABA machinery in the circumvallate papillae. Both GABA_A-and GABA_B- immunoreactivity were detected in the peripheral taste receptor cells. We also used transgenic mice that express green fluorescent protein (GFP) in either the Type II taste cells, which can respond to bitter, sweet or umami taste stimuli, or in the Type III GAD67 expressing taste cells. Thus, we were able to identify that GABAergic receptors are expressed in some Type II and Type III taste cells. Mouse GAT4 labeling was concentrated in the cells surrounding the taste buds with a few positively labeled TRCs at the margins of the taste buds.

Conclusions/Significance

The presence of GABAergic receptors localized on Type II and Type III taste cells suggests that GABA is likely modulating evoked taste responses in the mouse taste bud.

Ligand- and subunit-specific conformational changes in the ligand-binding domain and the TM2-TM3 linker of {alpha}1 {beta}2 {gamma}2 GABAA receptors

Cys-loop receptor ligand binding sites are located at subunit interfaces where they are lined by loops A-C from one subunit and loops D-F from the adjacent subunit. Agonist binding induces large conformational changes in loops C and F. However, it is controversial as to whether these conformational changes are essential for gating. Here we used voltage clamp fluorometry to investigate the roles of loops C and F in gating the α1 β2 γ2 GABA(A) receptor. Voltage clamp fluorometry involves labeling introduced cysteines with environmentally sensitive fluorophores and inferring structural rearrangements from ligand-induced fluorescence changes. Previous attempts to define the roles of loops C and F using this technique have focused on homomeric Cys-loop receptors. However, the problem with studying homomeric receptors is that it is difficult to eliminate the possibility of bound ligands interacting directly with attached fluorophores at the same site. Here we show that ligands binding to the β2-α1 interface GABA binding site produce conformational changes at the adjacent subunit interface. This is most likely due to agonist-induced loop C closure directly altering loop F conformation at the adjacent α1-β2 subunit interface. However, as antagonists and agonists produce identical α1 subunit loop F conformational changes, these conformational changes appear unimportant for gating. Finally, we demonstrate that TM2-TM3 loops from adjacent β2 subunits in α1 β2 receptors can dimerize via K24'C disulfides in the closed state. This result implies unexpected conformational mobility in this crucial part of the gating machinery. Together, this information provides new insights into the activation mechanisms of Cys-loop receptors.

The structural basis of function in Cys-loop receptors

Cys-loop receptors are membrane-spanning neurotransmitter-gated ion channels that are responsible for fast excitatory and inhibitory transmission in the peripheral and central nervous systems. The best studied members of the Cys-loop family are nACh, 5-HT3, GABAA and glycine receptors. All these receptors share a common structure of five subunits, pseudo-symmetrically arranged to form a rosette with a central ion-conducting pore. Some are cation selective (e.g. nACh and 5-HT3) and some are anion selective (e.g. GABAA and glycine). Each receptor has an extracellular domain (ECD) that contains the ligand-binding sites, a transmembrane domain (TMD) that allows ions to pass across the membrane, and an intracellular domain (ICD) that plays a role in channel conductance and receptor modulation. Cys-loop receptors are the targets for many currently used clinically relevant drugs (e.g. benzodiazepines and anaesthetics). Understanding the molecular mechanisms of these receptors could therefore provide the catalyst for further development in this field, as well as promoting the development of experimental techniques for other areas of neuroscience.In this review, we present our current understanding of Cys-loop receptor structure and function. The ECD has been extensively studied. Research in this area has been stimulated in recent years by the publication of high-resolution structures of nACh receptors and related proteins, which have permitted the creation of many Cys loop receptor homology models of this region. Here, using the 5-HT3 receptor as a typical member of the family, we describe how homology modelling and ligand docking can provide useful but not definitive information about ligand interactions. We briefly consider some of the many Cys-loop receptors modulators. We discuss the current understanding of the structure of the TMD, and how this links to the ECD to allow channel gating, and consider the roles of the ICD, whose structure is poorly understood. We also describe some of the current methods that are beginning to reveal the differences between different receptor states, and may ultimately show structural details of transitions between them.