γ-aminobutyric acid-type A (GABAA) 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 α1β2γ2 subtype GABAA 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 β3 homopentamer GABAA 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 α1β2γ2 GABAA 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.
Keywords: AChBP – Acetylcholine Binding Protein; BZDs – Benzodiazepines; CNS – Central Nervous System; ECD – Extracellular domain; ELIC – Erwinia ligand-gated ion channel; GABA – Gamma-Aminobutyric acid; GABAAR; GABAARs – Gamma-Aminobutyric Acid type A receptors; GLIC – Gloeobacter ligand-gated ion channel; MD – Molecular Dynamics; POPC – 1-palmitoyl-2-oleoyl-phophatidylcholine; SCAM – Substituted cysteine accessibility method; TMD – Transmembrane domain; benzodiazepines; docking; homology modelling; i-BZDs – Imidazo-Benzodiazepines; molecular dynamics.