Overview of mutations affecting the function of the GABA_AR or the binding of orthosteric ligands.

<p>All data in the table are in agreement with our GABA_AR model.</p

DZP binding mode.

<p>A) The assumed biologically active binding mode of DZP (gray) at the interface between the α₁ (smudge) and γ₂ (firebrick) subunits. In this conformation the C-3 points upwards and the pending phenyl substituent is directed inwards. B) Covalently attached DZP-NCS (cyan) overlaid with DZP (gray). Only moderate differences between the docked and the covalently attached ligands exist.</p

Conserved template residues.

<p>The figure shows residues that are conserved or homologous to GABA_AR binding site residues from the GluCl X-ray structure (PBD ID 3RIF) as grey sticks and the bacterial Cys-Loop receptor homolog, ELIC (PDB ID 2VL0) as purple sticks. Glutamate as co-crystallized with GluCl is shown in yellow, where the structure corresponding to GABA is shown as sticks and the α-carboxylic acid removed prior to homology modeling is shown as lines. A) GABA binding site and B) BZD binding site.</p

Some classical GABA_A receptor ligands.

<p>GABA is the endogenous GAB_AR agonist, muscimol a classical high-affinity agonist and THIP a muscimol analogue. Although not a GAB_AR ligand, glutamate is included to illustrate the resemblance to GABA. Diazepam (DZP) belongs to the benzodiazepine class of compounds, which are allosteric GABA_A modulators. The DZP-NCS analogue attaches covalently to GABA_ARs and is included for validation purposes.</p

Overview of mutations affecting the function of the GABA_AR or the ability of BZD binding site ligands to bind to the BZD binding site.

<p>All data in the table is in agreement with our GABA_AR model.</p

Conformational energy profile for dihedral drive of the amino-methyl side chain of muscimol.

<p>B3LYP/6-31G** energies.</p

Region suited for a tightly bound water molecule identified in agonist site.

<p>A GRID calculation at the agonist binding site, using the water probe, identified two regions of strong binding interaction energy (−11 kcal/mol). One region is overlapping with the acidic moiety of agonists and the other region is situated next to the backbone of β₂S156 and β₂Y157 (grey mesh). The calculation was performed in absence of agonist in the binding site. In the picture, the site has been optimized for muscimol as described in the methods section. A) When a water molecule is placed between muscimol and the B-loop backbone, perfect hydrogen bonding distances are obtained, resulting in optimal interactions between the high affinity ligand muscimol and the GABA receptor. B) When also GABA is included in the site, it is obvious that the water molecule would make up for the backbone interaction that GABA is predicted to make.</p

Alignment of protein sequences from GluCl, GLIC, ELIC and the human α₁, β₂, and γ₂ GABA_AR subunits.

<p>The GluCl sequence was used as template for homology modeling throughout the GABA_AR subunits, and ELIC was included as a template in the regions marked with blue boxes. The secondary structure deduced from the X-ray structure of GluCl is shown above the alignment (red shapes denote α-helices, and green arrows denote β-strands), whereas historically assigned loop regions are indicated below the alignment. As in the GluCl structure the M3–M4 intracellular loop of GABA_AR sequences was replaced by an AGT tri-peptide linker. Residues comprising the binding sites (within 8 Å of Glu in the GluCl structure and pointing towards the binding site) are colored pink (Glu and GABA binding site) and cyan (BZD binding site). Binding site residues conserved with respect to the templates are indicated as follows: <b>+</b> conserved in both GABA and BZD binding sites; <b>*</b> conserved in the GABA binding sites; <b>.</b> conserved in the BZD binding site. For details on calculation of binding site sequence identities see <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0052323#pone.0052323.s002" target="_blank">Figure S1</a>.</p

Illustration of the GABA_AR structural composition.

<p><b>A</b>) Top view showing the pentameric assembly of α₁, β₂ and γ₂ subunits and the location of binding sites for GABA and BZDs; <b>B</b>) Side view illustrating the extracellular domain (ECD) where agonists and benzodiazepines bind and the transmembrane domain (TMD); <b>C</b>) Zooming in on a GABA binding site at the subunit interface between β₂ and α₁ subunits, loop regions A–F mentioned in the text are shown (A: yellow, B: orange, C: red, D: purple, E: blue and F: pink).</p

Agonist binding modes determined by induced fit docking.

<p>A) GABA (green), B) THIP (pink) and C) muscimol (cyan) are shown in the orthosteric binding site at the interface between the β₂ subunit (teal) and the α₁ subunit (smudge).</p