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.

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.

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

Variants in γ‐aminobutyric acid A (GABA_A) 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.

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 GABA_A 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 GABA_A 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 GABA_A 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 GABA_A 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 GABA_A 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 GABA_A 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.

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

The family of GABA_A receptors is an important drug target group in the treatment of sleep disorders, anxiety, epileptic seizures, and many others. The most frequent GABA_A 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 GABA_A receptor subtypes equally. In the past years, however, drug development research has focused on studying α5-containing GABA_A receptors. Beyond the central nervous system, α5-containing GABA_A 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 GABA_A 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 GABA_A 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 GABA_A 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 GABA_A 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 GABA_A receptor α-subunit is the dominant molecular determinant of drug selectivity.