Molecular dynamics-derived pharmacophores of Schistosoma glutathione transferase in complex with bromosulfophthalein: Screening and analysis of potential inhibitors.

Schistosoma glutathione transferases (GSTs) have been identified as attractive drug targets for the design of novel antischistosomals. Here, we used in silico methods to validate the discriminative inhibitory properties of bromosulfophthalein (BSP) against the 26-kDa GST from S. japonicum (Sj26GST), and the 28-kDa GST from S. haematobium (Sh28GST), versus human GST (hGST) isoforms alpha (hGSTA), mu (hGSTM) and pi (hGSTP). The use of BSP as an archetypal selective inhibitor was harnessed to produce molecular dynamics-derived pharmacophores of the two targets. Pharmacophore-based screening using a large dataset of experimental and approved drug compounds was performed to produce a shortlist of candidates. The top candidate for each target was prioritised via molecular docking, yielding guanosine-3'-monophosphate-5'-diphosphate (G3D) for Sj26GST, and quercetin-3'-O-phosphate (Q3P) for Sh28GST. Comparative molecular dynamics studies of both candidates compared to BSP showed similar characteristics of binding stability and strength, suggesting their potential to emulate the inhibitory effects of BSP.

Comparative structural analysis of the human and Schistosoma glutathione transferase dimer interface using selective binding of bromosulfophthalein.

The binding channel of Schistosoma glutathione transferase (SGST) has been identified to possess a non-substrate site implicated in enzyme inhibition. This binding channel is formed by the interface of the GST dimer. We produced a comparative characterization of the SGST dimer interface with respect to that of human GST (hGST) analogues using the selective binding of bromosulfophthalein (BSP). First, two SGST and three hGST structures were used as search queries to assemble a data set of 48 empirical GST structures. Sequence alignment to generate a universal residue indexing scheme was then performed, followed by local superposition of the dimer interface. Principal component analysis revealed appreciable variation of the dimer interface, suggesting the potential for selective inhibition of SGST. BSP was found to dock invariably in the dimer interface core pocket, placing it in proximity to the GST catalytic domains, through which it may exert its inhibitory behavior. Binding poses across the GST forms were distinguished with ligand interaction profiling, where SGST complexes showed stabilization of ligand aromatic- and sulfonate moieties, which altogether anchor the ligand and produce a tight association. In comparison, missing aromatic stabilization in the hGST complexes impart large bonding distances, causing mobile poses likely to dissociate. Altogether, this study illustrates the potential for selective inhibition of SGST, rationalizes the selective behavior of the BSP inhibitor, and produces a reliable metric for construction and validation of pharmacophore models of the SGST binding channel.