GSFM: A genome-scale functional module transformation to represent drug efficacy for in silico drug discovery

Pharmacotranscriptomic profiles, which capture drug-induced changes in gene expression, offer vast potential for computational drug discovery and are widely used in modern medicine. However, current computational approaches neglected the associations within gene-gene functional networks and unrevealed the systematic relationship between drug efficacy and the reversal effect. Here, we developed a new genome-scale functional module (GSFM) transformation framework to quantitatively evaluate drug efficacy for in silico drug discovery. GSFM employs four biologically interpretable quantifiers: GSFM_Up, GSFM_Down, GSFM_ssGSEA, and GSFM_TF to comprehensively evaluate the multi-dimension activities of each functional module (FM) at gene-level, pathway-level, and transcriptional regulatory network-level. Through a data transformation strategy, GSFM effectively converts noisy and potentially unreliable gene expression data into a more dependable FM active matrix, significantly outperforming other methods in terms of both robustness and accuracy. Besides, we found a positive correlation between RS_GSFM and drug efficacy, suggesting that RS_GSFM could serve as representative measure of drug efficacy. Furthermore, we identified WYE-354, perhexiline, and NTNCB as candidate therapeutic agents for the treatment of breast-invasive carcinoma, lung adenocarcinoma, and castration-resistant prostate cancer, respectively. The results from in vitro and in vivo experiments have validated that all identified compounds exhibit potent anti-tumor effects, providing proof-of-concept for our computational approach.