Glycine transporters are interesting therapeutic targets as they play significant roles in glycinergic and glutamatergic systems. The search for new selective inhibitors of particular types of glycine transporters (GlyT-1 and GlyT-2) with beneficial kinetics is hampered by limited knowledge about the spatial structure of these proteins. In this study, a pool of homology models of GlyT-1 and GlyT-2 in different conformational states was constructed using the crystal structures of related transporters from the SLC6 family and the recently revealed structure of GlyT-1 in the inward-open state, in order to investigate their binding sites. The binding mode of the known GlyT-1 and GlyT-2 inhibitors was determined using molecular docking studies, molecular dynamics simulations, and MM-GBSA free energy calculations. The results of this study indicate that two amino acids, Gly373 and Leu476 in GlyT-1 and the corresponding Ser479 and Thr582 in GlyT-2, are mainly responsible for the selective binding of ligands within the S1 site. Apart from these, one pocket of the S2 site, which lies between TM3 and TM10, may also be important. Moreover, selective binding of noncompetitive GlyT-1 inhibitors in the intracellular release pathway is affected by hydrophobic interactions with Ile399, Met382, and Leu158. These results can be useful in the rational design of new glycine transporter inhibitors with desired selectivity and properties in the future.
Keywords: GlyT-1 inhibitors; GlyT-2 inhibitors; glycine transporters modeling; molecular docking; molecular dynamics.