Machine Learning-based Virtual Screening for STAT3 Anticancer Drug Target

Signal transducers and activators of the transcription (STAT) family is composed of seven structurally similar and highly conserved members, including STAT1, STAT2, STAT3, STAT4, STAT5a, STAT5b, and STAT6. The STAT3 signaling cascade is activated by upstream kinase signals and undergoes phosphorylation, homo-dimerization, nuclear translocation, and DNA binding, resulting in the expression of target genes involved in tumor cell proliferation, metastasis, angiogenesis, and immune editing. STAT3 hyperactivation has been documented in a number of tumors, including head and neck, breast, lung, liver, kidney, prostate, pancreas cancer, multiple myeloma, and acute myeloid leukemia. Drug discovery is a timeconsuming and costly process; it may take ten to fifteen years to bring a single drug to the market. Machine learning algorithms are very fast and effective and commonly used in the field, such as drug discovery. These algorithms are ideal for the virtual screening of large compound libraries to classify molecules as active or inactive.

The present work aims to perform machine learning-based virtual screening for the STAT3 drug target.

Machine learning models, such as k-nearest neighbor, support vector machine, Gaussian naïve Bayes, and random forest for classifying the active and inactive inhibitors against a STAT3 drug target, were developed. Ten-fold cross-validation was used for model validation. Then the test dataset prepared from the zinc database was screened using the random forest model. A total of 20 compounds with 88% accuracy was predicted as active against STAT3. Furthermore, these twenty compounds were docked into the active site of STAT3. The two complexes with good docking scores as well as the reference compound were subjected to MD simulation. A total of 100ns MD simulation was performed.

Compared to all other models, the random forest model revealed better results. Compared to the standard reference compound, the top two hits revealed greater stability and compactness.

In conclusion, our predicted hits have the ability to inhibit STAT3 overexpression to combat STAT3-associated diseases.