In silico simulations of STAT1 and STAT3 inhibitors predict SH2 domain cross-binding specificity

Signal transducers and activators of transcription (STATs) comprise a family of transcription factors that are structurally related and which participate in signaling pathways activated by cytokines, growth factors and pathogens. Activation of STAT proteins is mediated by the highly conserved Src homology 2 (SH2) domain, which interacts with phosphotyrosine motifs for specific contacts between STATs and receptors and for STAT dimerization. By generating new models for human (h)STAT1, hSTAT2 and hSTAT3 we applied comparative in silico docking to determine SH2-binding specificity of the STAT3 inhibitor stattic, and of fludarabine (STAT1 inhibitor). Thus, we provide evidence that by primarily targeting the highly conserved phosphotyrosine (pY+0) SH2 binding pocket stattic is not a specific hSTAT3 inhibitor, but is equally effective towards hSTAT1 and hSTAT2. This was confirmed in Human Micro-vascular Endothelial Cells (HMECs) in vitro, in which stattic inhibited interferon-α-induced phosphorylation of all three STATs. Likewise, fludarabine inhibits both hSTAT1 and hSTAT3 phosphorylation, but not hSTAT2, by competing with the highly conserved pY+0 and pY-X binding sites, which are less well-preserved in hSTAT2. Moreover we observed that in HMECs in vitro fludarabine inhibits cytokine and lipopolysaccharide-induced phosphorylation of hSTAT1 and hSTAT3 but does not affect hSTAT2. Finally, multiple sequence alignment of STAT-SH2 domain sequences confirmed high conservation between hSTAT1 and hSTAT3, but not hSTAT2, with respect to stattic and fludarabine binding sites. Together our data offer a molecular basis that explains STAT cross-binding specificity of stattic and fludarabine, thereby questioning the present selection strategies of SH2 domain-based competitive small inhibitors.