Dual-inhibitors of STAT5 and STAT3: studies from molecular docking and molecular dynamics simulations

Although molecularly targeted therapy with imatinib has improved treatments of chronic myeloid leukemia (CML), clinical resistance gradually develops in patients with accelerated or blast phase CML. The inability of imatinib to cure CML suggests that inactivation of BCR-ABL kinase activity alone is not sufficient to control the disease. Aberrant STAT signaling and constitutive STAT5 or STAT3 activation are frequently found in both acute and chronic leukemia. Constitutive activation of STAT5 and STAT3 are associated with imatinib resistance on leukemia cells. Development of drugs targeting SH2 domains of STAT5 and STAT3 provides a novel strategy for the treatment of the imatinib-resistant CML. Here, molecular docking and molecular dynamics simulations were used to investigate the interactions of the drugs targeting STAT3 and STAT5 receptors at molecular level. The calculated binding free energies are consistent with the ranking of the experimental affinities and our simulations also explained their differences in binding energy. Then virtual screening based on molecular docking and molecular dynamics was applied to screen a set of ~1500 compounds for dual inhibitors of the SH2 domains of STAT5 and STAT3. Three top score compounds obtained in virtual screening were compound 660, 304, and 561. Results show that the three predicted dual-inhibitors are well fitted within the two binding domains and are predicted to present improved STAT5 and STAT3 SH2 inhibitory activity.