A Unified Model of the GABAA Receptor Comprising Agonist and Benzodiazepine Binding Sites

We present a full-length α₁β₂γ₂ 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 α₁R66, β₂T202, α₁T129, β₂E155, β₂Y205 and the backbone of β₂S156. Muscimol is predicted to bind similarly, however, with minor differences rationalized with quantum mechanical energy calculations. Muscimol key interactions are predicted to be α₁R66, β₂T202, α₁T129, β₂E155, β₂Y205 and β₂F200. Furthermore, we argue that a water molecule could mediate further interactions between muscimol and the backbone of β₂S156 and β₂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 α₁T206 and γ₂T142, similar to the acidic moiety of GABA. The chlorine atom of DZP is placed near the important α₁H101 and the N-methyl group near α₁Y159, α₁T206, and α₁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.