Rhipicephalus (Boophilus) microplus is one of the most widespread ticks causing a massive loss to livestock production. The long-term use of acaracides rapidly develops acaracide resistance. In R. microplus, enhancing the metabolic activity of glutathione S-transferase (RmGST) is one of the mechanisms underlying acaracide resistance. RmGST catalyzes the conjugation of glutathione (GSH) to insecticides causing an easy-to-excrete conjugate. The active RmGST dimer contains two active sites (hydrophobic co-substrate binding site (H-site) and GSH binding site (G-site)) in each monomer. To preserve the insecticide efficacy, s-hexyl glutathione (GTX), a GST inhibitor, has been used as a synergist. To date, no molecular information on the RmGST-GSH/GTX complex is available. The insight is important for developing a novel RmGST inhibitor. Therefore, in this work, molecular dynamics simulations (MD) were performed to explore the binding of GTX and GSH to RmGST. GSH binds tighter and sits rigidly inside the G-site, while flexible GTX occupies both active sites. In GSH, the backbone mainly interacts with W8, R43, W46, K50, N59, L60, Q72, and S73, while its thiol group directs to Y7. In contrast, the aliphatic hexyl of GTX protrudes into the H-site and allows a flexible peptide core to form various interactions. Such high GTX flexibility and the protrusion of its hexyl moiety to the H-site suggest the dual role of GTX in preventing the conjugation reaction and the binding of acaracide. This insight can provide a better understanding of an important insecticide-resistance mechanism, which may in turn facilitate the development of novel approaches to tick control.
Keywords: MD simulation; R. microplus; glutathione; glutathione s-transferase; s-hexyl glutathione.