Abstract
The recently obtained cryo-electron microscopy structure (PDB code: 7SK2) of the human γ-aminobutyric acid transporter type 1 (hGAT-1) in complex with the antiepileptic drug, tiagabine, revealed a rather unexpected binding mode for this inhibitor in an inward-open state of the transporter. The simultaneously released crystal structures of the modified dopamine transporter with mutations mimicking hGAT-1 indicated an alternative binding mode for the tiagabine analogues that were found to block the transporter in an outward-open state, which is more consistent with the results of previous biological and molecular modeling studies. In view of the above discrepancies, our study compares different hypothetical tiagabine binding modes using classical and accelerated molecular dynamics simulations, as well as MM-GBSA free binding energy (dG) calculations. The results indicate that the most stable and energetically favorable binding mode of tiagabine is the one where the nipecotic acid fragment is located in the main binding site (S1) and the aromatic rings are arranged within the S2 site of the hGAT-1 transporter in an outward-open state, confirming the previous molecular modelling findings. The position of tiagabine bound to hGAT-1 in an inward-open state, partially within the intracellular release pathway, was significantly less stable and the dG values calculated for this complex were higher. Furthermore, analysis of the cryo-electron map for the 7SK2 structure shows that the model does not appear to fit into the map optimally at the ligand binding site. These findings suggest that the position of tiagabine found in the 7SK2 structure is rather ambiguous and requires further experimental verification. The identification of the main, high-affinity binding site for tiagabine and its analogues is crucial for the future rational design of the GABA transporter inhibitors.