GABAA receptor-mediated seizure liabilities: a mixed-methods screening approach.

GABAA receptors, members of the pentameric ligand-gated ion channel superfamily, are widely expressed in the central nervous system and mediate a broad range of pharmaco-toxicological effects including bidirectional changes to seizure threshold. Thus, detection of GABAA receptor-mediated seizure liabilities is a big, partly unmet need in early preclinical drug development. This is in part due to the plethora of allosteric binding sites that are present on different subtypes of GABAA receptors and the critical lack of screening methods that detect interactions with any of these sites. To improve in silico screening methods, we assembled an inventory of allosteric binding sites based on structural data. Pharmacophore models representing several of the binding sites were constructed. These models from the NeuroDeRisk IL Profiler were used for in silico screening of a compiled collection of drugs with known GABAA receptor interactions to generate testable hypotheses. Amoxapine was one of the hits identified and subjected to an array of in vitro assays to examine molecular and cellular effects on neuronal excitability and in vivo locomotor pattern changes in zebrafish larvae. An additional level of analysis for our compound collection is provided by pharmacovigilance alerts using FAERS data. Inspired by the Adverse Outcome Pathway framework, we postulate several candidate pathways leading from specific binding sites to acute seizure induction. The whole workflow can be utilized for any compound collection and should inform about GABAA receptor-mediated seizure risks more comprehensively compared to standard displacement screens, as it rests chiefly on functional data.

A de novo missense variant in GABRA4 alters receptor function in an epileptic and neurodevelopmental phenotype.

Variants in γ-aminobutyric acid A (GABAA ) receptor genes cause different forms of epilepsy and neurodevelopmental disorders. To date, GABRA4, encoding the α4-subunit, has not been associated with a monogenic condition. However, preclinical evidence points toward seizure susceptibility. Here, we report a de novo missense variant in GABRA4 (c.899C>T, p.Thr300Ile) in an individual with early-onset drug-resistant epilepsy and neurodevelopmental abnormalities. An electrophysiological characterization of the variant, which is located in the pore-forming domain, shows accelerated desensitization and a lack of seizure-protective neurosteroid function. In conclusion, our findings strongly suggest an association between de novo variation in GABRA4 and a neurodevelopmental disorder with epilepsy.

A Benzodiazepine Ligand with Improved GABAA Receptor α5-Subunit Selectivity Driven by Interactions with Loop C.

The family of GABAA receptors is an important drug target group in the treatment of sleep disorders, anxiety, epileptic seizures, and many others. The most frequent GABAA receptor subtype is composed of two α-, two ß-, and one γ2-subunit, whereas the nature of the α-subunit critically determines the properties of the benzodiazepine binding site of those receptors. Nearly all of the clinically relevant drugs target all GABAA receptor subtypes equally. In the past years, however, drug development research has focused on studying α5-containing GABAA receptors. Beyond the central nervous system, α5-containing GABAA receptors in airway smooth muscles are considered as an emerging target for bronchial asthma. Here, we investigated a novel compound derived from the previously described imidazobenzodiazepine SH-053-2'F-R-CH3 (SH53d-ester). Although SH53d-ester is only moderately selective for α5-subunit-containing GABAA receptors, the derivative SH53d-acid shows superior (>40-fold) affinity selectivity and is a positive modulator. Using two-electrode voltage clamp electrophysiology in Xenopus laevis oocytes and radioligand displacement assays with human embryonic kidney 293 cells, we demonstrated that an acid group as substituent on the imidazobenzodiazepine scaffold leads to large improvements of functional and binding selectivity for α5ß3γ2 over other αxß3γ2 GABAA receptors. Atom level structural studies provide hypotheses for the improved affinity to this receptor subtype. Mutational analysis confirmed the hypotheses, indicating that loop C of the GABAA receptor α-subunit is the dominant molecular determinant of drug selectivity. Thus, we characterize a promising novel α5-subunit-selective drug candidate. SIGNIFICANCE STATEMENT: In the current study we present the detailed pharmacological characterization of a novel compound derived from the previously described imidazobenzodiazepine SH-053-2'F-R-CH3. We describe its superior (>40-fold) affinity selectivity for α5-containing GABAA receptors and show atom-level structure predictions to provide hypotheses for the improved affinity to this receptor subtype. Mutational analysis confirmed the hypotheses, indicating that loop C of the GABAA receptor α-subunit is the dominant molecular determinant of drug selectivity.

GABAA Receptor Ligands Often Interact with Binding Sites in the Transmembrane Domain and in the Extracellular Domain-Can the Promiscuity Code Be Cracked?

Many allosteric binding sites that modulate gamma aminobutyric acid (GABA) effects have been described in heteropentameric GABA type A (GABAA) receptors, among them sites for benzodiazepines, pyrazoloquinolinones and etomidate. Diazepam not only binds at the high affinity extracellular "canonical" site, but also at sites in the transmembrane domain. Many ligands of the benzodiazepine binding site interact also with homologous sites in the extracellular domain, among them the pyrazoloquinolinones that exert modulation at extracellular α+/ß- sites. Additional interaction of this chemotype with the sites for etomidate has also been described. We have recently described a new indole-based scaffold with pharmacophore features highly similar to pyrazoloquinolinones as a novel class of GABAA receptor modulators. Contrary to what the pharmacophore overlap suggests, the ligand presented here behaves very differently from the identically substituted pyrazoloquinolinone. Structural evidence demonstrates that small changes in pharmacophore features can induce radical changes in ligand binding properties. Analysis of published data reveals that many chemotypes display a strong tendency to interact promiscuously with binding sites in the transmembrane domain and others in the extracellular domain of the same receptor. Further structural investigations of this phenomenon should enable a more targeted path to less promiscuous ligands, potentially reducing side effect liabilities.

Probes for the heterogeneity of muscimol binding sites in rat brain.

Introduction: The plant-based alkaloid muscimol is a potent agonist of inhibitory GABAA-neurotransmitter receptors. GABAA receptors are a heterogeneous family of pentameric complexes, with 5 out of 19 subunits assembling around the central anion pore. Muscimol is considered to bind to all receptor subtypes at the orthosteric drug binding site at the ß+/α- interface. Recently, we observed that the antipsychotic drugs clozapine (CLZ), loxapine (LOX) and chlorpromazine (CPZ) although exerting functional inhibition on multiple GABAA receptor subtypes showed diverging results in displacing 3H-muscimol. While a complete displacement could be observed in hippocampal membranes by bicuculline (BIC), and no displacement with CPZ, the compounds CLZ and LOX competed partially. Non-sigmoidal, complex dose response curves were indicative of multiple sites. In the current study we now aimed to investigate more extensively this heterogeneity of bicuculline sensitive muscimol sites in rat brain. Methods: We tested membranes from four different brain regions (hippocampus, cerebellum, thalamus and striatum) and selected recombinantly expressed subunit combinations with displacement assays. 3H-muscimol displacement was tested with BIC, LOX, CLZ and CPZ. In silico ligand structural analysis and computational docking was performed. Results: We observed a unique pharmacology of each tested compound in the studied brain regions. Combining two of the tested ligands suggests that in striatum all CLZ sites are contained in the pool of LOX sites, while the CPZ sites may in part be non-overlapping with LOX sites. Experiments on recombinantly expressed receptors indicate, that BIC can displace 3H-muscimol from all tested receptors, while LOX and CLZ display different and variable competition indicative of multiple sites. Molecular docking produced structural correlates of the observed diversity of muscimol sites on the basis of bicuculline bound experimental structures. Discussion: These findings indicate that 3H-muscimol binding sites in rat brain are heterogeneous, with different populations of receptors, which are CPZ, LOX or CLZ sensitive or insensitive. These binding sites show a varying distribution in different rat brain regions. Molecular docking suggests that the so-called loop F region of α subunits drives the observed differences.

Navigating the complex landscape of benzodiazepine- and Z-drug diversity: insights from comprehensive FDA adverse event reporting system analysis and beyond.

Introduction: Medications which target benzodiazepine (BZD) binding sites of GABAA receptors (GABAARs) have been in widespread use since the nineteen-sixties. They carry labels as anxiolytics, hypnotics or antiepileptics. All benzodiazepines and several nonbenzodiazepine Z-drugs share high affinity binding sites on certain subtypes of GABAA receptors, from which they can be displaced by the clinically used antagonist flumazenil. Additional binding sites exist and overlap in part with sites used by some general anaesthetics and barbiturates. Despite substantial preclinical efforts, it remains unclear which receptor subtypes and ligand features mediate individual drug effects. There is a paucity of literature comparing clinically observed adverse effect liabilities across substances in methodologically coherent ways. Methods: In order to examine heterogeneity in clinical outcome, we screened the publicly available U.S. FDA adverse event reporting system (FAERS) database for reports of individual compounds and analyzed them for each sex individually with the use of disproportionality analysis. The complementary use of physico-chemical descriptors provides a molecular basis for the analysis of clinical observations of wanted and unwanted drug effects. Results and Discussion: We found a multifaceted FAERS picture, and suggest that more thorough clinical and pharmacoepidemiologic investigations of the heterogenous side effect profiles for benzodiazepines and Z-drugs are needed. This may lead to more differentiated safety profiles and prescription practice for particular compounds, which in turn could potentially ease side effect burden in everyday clinical practice considerably. From both preclinical literature and pharmacovigilance data, there is converging evidence that this very large class of psychoactive molecules displays a broad range of distinctive unwanted effect profiles - too broad to be explained by the four canonical, so-called "diazepam-sensitive high-affinity interaction sites". The substance-specific signatures of compound effects may partly be mediated by phenomena such as occupancy of additional binding sites, and/or synergistic interactions with endogenous substances like steroids and endocannabinoids. These in turn drive the wanted and unwanted effects and sex differences of individual compounds.

Cys-loop receptors on cannabinoids: All high?

Endocannabinoids (eCBS) are endogenously derived lipid signaling molecules that serve as tissue hormones and interact with multiple targets, mostly within the endocannabinoid system (ECS). The ECS is a highly conserved regulatory system involved in homeostatic regulation, organ formation, and immunomodulation of chordates. The term "cannabinoid" evolved from the distinctive class of plant compounds found in Cannabis sativa, an ancient herb, due to their action on CB1 and CB2 receptors. CB1/2 receptors are the primary targets for eCBs, but their effects are not limited to the ECS. Due to the high interest and extensive research on the ECS, knowledge on its constituents and physiological role is substantial and still growing. Crosstalk and multiple targeting of molecules are common features of endogenous and plant compounds. Cannabimimetic molecules can be divided according to their origin, natural or synthetic, including phytocannabinoids (pCB's) or synthetic cannabinoids (sCB's). The endocannabinoid system (ECS) consists of receptors, transporters, enzymes, and signaling molecules. In this review, we focus on the effects of cannabinoids on Cys-loop receptors. Cys-loop receptors belong to the class of membrane-bound pentameric ligand gated ion channels, each family comprising multiple subunits. Mammalians possess GABA type A receptors (GABAAR), glycine receptors (GlyR), serotonin receptors type 3 (5-HT3R), and nicotinic acetylcholine receptors (nAChR). Several studies have shown different modulatory effects of CBs on multiple members of the Cys-loop receptor family. We highlight the existing knowledge, especially on subunits and protein domains with conserved binding sites for CBs and their possible pharmacological and physiological role in epilepsy and in chronic pain. We further discuss the potential for cannabinoids as first line treatments in epilepsy, chronic pain and other neuropsychiatric conditions, indicated by their polypharmacology and therapeutic profile.

Molecular Mingling: Multimodal Predictions of Ligand Promiscuity in Pentameric Ligand-Gated Ion Channels.

Background: Human pentameric ligand-gated ion channels (pLGICs) comprise nicotinic acetylcholine receptors (nAChRs), 5-hydroxytryptamine type 3 receptors (5-HT3Rs), zinc-activated channels (ZAC), γ-aminobutyric acid type A receptors (GABAARs) and glycine receptors (GlyRs). They are recognized therapeutic targets of some of the most prescribed drugs like general anesthetics, anxiolytics, smoking cessation aids, antiemetics and many more. Currently, approximately 100 experimental structures of pLGICs with ligands bound exist in the protein data bank (PDB). These atomic-level 3D structures enable the generation of a comprehensive binding site inventory for the superfamily and the in silico prediction of binding site properties. Methods: A panel of high throughput in silico methods including pharmacophore screening, conformation analysis and descriptor calculation was applied to a selection of allosteric binding sites for which in vitro screens are lacking. Variant abundance near binding site forming regions and computational docking complement the approach. Results: The structural data reflects known and novel binding sites, some of which may be unique to individual receptors, while others are broadly conserved. The membrane spanning domain, comprising four highly conserved segments, contains ligand interaction sites for which in vitro assays suitable for high throughput screenings are critically lacking. This is also the case for structurally more variable novel sites in the extracellular domain. Our computational results suggest that the phytocannabinoid Δ9-tetrahydrocannabinol (Δ9-THC) can utilize multiple pockets which are likely to exist on most superfamily members. Conclusion: With this study, we explore the potential for polypharmacology among pLGICs. Our data suggest that ligands can display two forms of promiscuity to an extent greater than what has been realized: 1) Ligands can interact with homologous sites in many members of the superfamily, which bears toxicological relevance. 2) Multiple pockets in distinct localizations of individual receptor subtypes share common ligands, which counteracts efforts to develop selective agents. Moreover, conformational states need to be considered for in silico drug screening, as certain binding sites display considerable flexibility. In total, this work contributes to a better understanding of polypharmacology across pLGICs and provides a basis for improved structure guided in silico drug development and drug derisking.

Novel alpha6 preferring GABA-A receptor ligands based on loreclezole.

The family of GABA-A receptors contains nineteen mammalian subunits from which pentameric, GABA gated anion channels are assembled. The subunit encoded by the GABRA6 gene is highly expressed in the cerebellum and the receptors to which it contributes have recently been demonstrated to be a promising candidate as a novel drug target. Here we examined a series of loreclezole derivatives for potentially selective action at α6ß3γ2 receptors with the help of computational methods and functional testing with the two-electrode voltage clamp technique. The synthetic routes to some previously published ligands were improved, and a new derivative was synthesized based on computational docking results. This new loreclezole derivative, [3-(2-chloro-4-methylphenyl)-3-methylbutanenitrile (40)], was shown to display stronger modulatory action in concatenated α6ß3γ2 receptors compared to their α1ß3γ2 counterpart. The hypothetical bound state structure provides valuable guidance for future design of selective therapeutics.

Tricyclic antipsychotics and antidepressants can inhibit α5-containing GABAA receptors by two distinct mechanisms.

BACKGROUND AND PURPOSE: Many psychotherapeutic drugs, including clozapine, display polypharmacology and act on GABAA receptors. Patients with schizophrenia show alterations in function, structure and molecular composition of the hippocampus, and a recent study demonstrated aberrant levels of hippocampal α5 subunit-containing GABAA receptors. The purpose of this study is to investigate the effects of tricyclic compounds on α5 subunit-containing receptor subtypes. EXPERIMENTAL APPROACH: Functional studies of effects by seven antipsychotic and antidepressant medications were performed in several GABAA receptor subtypes by two-electrode voltage-clamp electrophysiology using Xenopus laevis oocytes. Computational structural analysis was employed to design mutated constructs of the α5 subunit, probing a novel binding site. Radioligand displacement data complemented the functional and mutational findings. KEY RESULTS: The antipsychotic drugs clozapine and chlorpromazine exerted functional inhibition on multiple GABAA receptor subtypes, including those containing α5-subunits. Based on a chlorpromazine binding site observed in a GABA-gated bacterial homologue, we identified a novel site in α5 GABAA receptor subunits and demonstrate differential usage of this and the orthosteric sites by these ligands. CONCLUSION AND IMPLICATIONS: Despite high molecular and functional similarities among the tested ligands, they reduce GABA currents by differential usage of allosteric and orthosteric sites. The chlorpromazine site we describe here is a new potential target for optimizing antipsychotic medications with beneficial polypharmacology. Further studies in defined subtypes are needed to substantiate mechanistic links between the therapeutic effects of clozapine and its action on certain GABAA receptor subtypes.

Structure-Guided Computational Methods Predict Multiple Distinct Binding Modes for Pyrazoloquinolinones in GABAA Receptors.

Pyrazoloquinolinones (PQs) are a versatile class of GABAA receptor ligands. It has been demonstrated that high functional selectivity for certain receptor subtypes can be obtained by specific substitution patterns, but so far, no clear SAR rules emerge from the studies. As is the case for many GABAA receptor targeting chemotypes, PQs can interact with distinct binding sites on a given receptor pentamer. In pentamers of αßγ composition, such as the most abundant α1ß2γ2 subtype, many PQs are high affinity binders of the benzodiazepine binding site at the extracellular α+/γ2- interfaces. There they display a functionally near silent, flumazenil-like allosteric activity. More recently, interactions with extracellular α+/ß- interfaces have been investigated, where strong positive modulation can be steered toward interesting subtype preferences. The most prominent examples are functionally α6-selective PQs. Similar to benzodiazepines, PQs also seem to interact with sites in the transmembrane domain, mainly the sites used by etomidate and barbiturates. This promiscuity leads to potential contributions from multiple sites to net modulation. Developing ligands that interact exclusively with the extracellular α+/ß- interfaces would be desired. Correlating functional profiles with binding sites usage is hampered by scarce and heterogeneous experimental data, as shown in our meta-analysis of aggregated published data. In the absence of experimental structures, bound states can be predicted with pharmacophore matching methods and with computational docking. We thus performed pharmacophore matching studies for the unwanted sites, and computational docking for the extracellular α1,6+/ß3- interfaces. The results suggest that PQs interact with their binding sites with diverse binding modes. As such, rational design of improved ligands needs to take a complex structure-activity landscape with branches between sub-series of derivatives into account. We present a workflow, which is suitable to identify and explore potential branching points on the structure-activity landscape of any small molecule chemotype.