DeepDTF: Dual-Branch Transformer Fusion for Multi-Omics Anticancer Drug Response Prediction
Authors: Yuhan Zhao, Jacob Tennant, James Yang, Zhishan Guo, Young Whang, Ning Sui
Organizations: Department of Computer Science, North Carolina State University, Raleigh, NC, USA · UNC School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA · Department of Molecular and Structural Biochemistry, North Carolina State University, Raleigh, NC, USA
Cancer drug response varies widely across tumors due to multi-layer molecular heterogeneity, motivating computational decision support for precision oncology. Despite recent progress in deep CDR models, robust alignment between high-dimensional multi-omics and chemically structured drugs remains challenging due to cross-modal misalignment and limited inductive bias. We present DeepDTF, an end-to-end dual-branch Transformer fusion framework for joint log(IC50) regression and drug sensitivity classification. The cell-line branch uses modality-specific encoders for multi-omics profiles with Transformer blocks to capture long-range dependencies, while the drug branch represents compounds as molecular graphs and encodes them with a GNN-Transformer to integrate local topology with global context. Omics and drug representations are fused by a Transformer-based module that models cross-modal interactions and mitigates feature misalignment. On public pharmacogenomic benchmarks under 5-fold cold-start cell-line evaluation, DeepDTF consistently outperforms strong baselines across omics settings, achieving up to RMSE=1.248, R^2=0.875, and AUC=0.987 with full multi-omics inputs, while reducing classification error (1-ACC) by 9.5%. Beyond accuracy, DeepDTF provides biologically grounded explanations via SHAP-based gene attributions and pathway enrichment with pre-ranked GSEA.
Figures & tables
Fig. 1: Overall workflow of the DeepDTF, including (i) multi-omics and drug-structure inputs, (ii) transformer model (detailed in Fig. 2 ), (iii) dual task prediction (IC50 regression and sensitivity classification), and (iv) an interpretability module based on SHAP → GSEA. SHAP produces a signed gene ranking, which is used for pre-ranked GSEA: The top panel illustrates a GSEA enrichment plot (running ES with pathway “hit” positions along the ranked list), and the bottom panel shows the corresponding null/random ES distribution for assessing enrichment direction and significance.
Fig. 2: Detailed architecture of DeepDTF. The left block is a dual-branch feature extractor for multi-omics and drug structure. The right block is a feature fusion joint prediction module with a Fusion-Transformer and two task-specific heads for IC50 regression and sensitivity classification.
Setting
Method
RMSE
R2
PCC
ACC
SEN
SPEC
AUC
MUT+CNV
tCNN [ 13 ]
1.459
0.725
0.851
0.952
0.706
0.977
0.963
DeepCDR [ 14 ]
1.401
0.770
0.878
0.953
0.732
0.975
0.969
DeepTTA [ 17 ]
1.420
0.775
0.880
0.953
0.732
0.978
0.970
DeepDTF
1.403
0.780
0.883
0.954
0.738
0.978
0.975
MUT+CNV+GE
CDRscan [ 12 ]
1.483
0.715
0.846
0.950
0.667
0.979
0.960
tCNN [ 13 ]
1.459
0.725
0.851
0.952
0.706
0.977
0.963
TABLE I: Evaluation of the Proposed Method and Baselines
Fig. 3: Ablation studies on omics integration and model components, where (a) and (c) show the regression performance under different combinations of omics inputs and model variants; while (b) and (d) report the classification metrics obtained from various omics integration strategies and model variants.
Fig. 4: Cartesian-regime biological explanations via SHAP–GSEA. (a) Top genes with positive SHAP (sensitivity-associated). (b) Bottom genes with negative SHAP (resistance-associated). (c) Top enriched Hallmark pathways from signed-SHAP rankings (positive ES; sensitivity-associated programs). (d) Top enriched Hallmark pathways from negative signed-SHAP rankings (negative ES; resistance-associated programs)
Predicting how a cell's transcriptome responds to a drug it has never seen is a core, hard problem in computational cell biology: recent benchmarks show complex models often fail to beat trivial baselines once test compounds are held out by chemistry. We study one cell line and assay, THP-1 cells profiled by DRUG-seq, scored by the active-compound weighted MSE(wMSE) of the VCPI prediction contest. We propose a staged approach: dumb baselines (untreated control and mean training-compound response) that the field keeps failing to beat; non-parametric retrieval (a Tanimoto-weighted average of a held-out compound's nearest training compounds); and a fusion stage combining a frozen chemistry embedding with retrieval-support features to predict the residual over the mean, with an uncertainty head and gene programs. On the released VCPI THP-1 drug-seq data (14,026 training compounds), under a Bemis-Murcko scaffold split, the model ranking inverts depending on the metric. Under an inverse-variance per-gene proxy, a regularized linear regression on Morgan fingerprints appears to win over the deep models, retrieval, and ChemBERTa -- the textbook "simple baselines win" result. But under the contest's true active-set metric (per-(gene, compound) Mejia weights, validated against the official scorer; mean baseline 0.535 vs the organizers' 0.507 reference), that reverses: the deep models win, our fusion decoder significantly beats the linear fingerprint baseline (-0.012 wMSE, paired bootstrap p < 10^-4), and the proxy's winner becomes the worst chemistry-aware predictor. Picking the metric picks the winner -- to our knowledge the first demonstration on real held-out drug chemistry of the metric-calibration effect established largely on genetic perturbation. We release a reproducible pipeline wired to the official scorer that emits a valid submission over the real 1064 x 12,995 grid.
Predicting cancer drug response from transcriptomic profiles is a cornerstone of precision oncology, yet the scientific value of machine learning models hinges not solely on predictive accuracy, but also on their capacity to generate reliable biological insights. Current explainability approaches in this setting are computationally costly, lack robustness, and reduce complex drug response to univariate gene importance scores, overlooking the coordinated gene activity that drives sensitivity and resistance. In this work, we present ILLUME+, a scalable post-hoc explainability framework that moves beyond single-gene assessments to capture multiple, complementary forms of explanation. Integrated into our end-to-end pipeline, ILLUME+ produces more stable gene importance scores than existing baselines, recovers established drug-gene associations and mechanisms of action, and enables AI-assisted hypothesis generation to uncover novel interaction-driven molecular signals in cancer biology.
Accurate prediction of patient-specific therapeutic response from pre-treatment transcriptomes is hindered by the scarcity of matched clinical response labels and post-treatment molecular profiles. Preclinical transfer-learning models can simulate drug-induced expression changes but are often hard to interpret and unstable, whereas knowledge-graph methods provide mechanistic context yet remain static and fail to capture drug-induced transcriptomic perturbation dynamics. We propose PREDIKTOR, a patient-centered multi-view framework that aligns a personalized network view with a transferable transcriptomic perturbation view to predict clinical drug response. For each patient, we construct an individualized gene regulatory network from tumor expression using DysRegNet and augment it with drug-target links from DrugBank; a graph neural encoder yields a drug-centric, mechanistically grounded embedding. In parallel, a frozen condition-specific gene-gene attention model pretrained on LINCS L1000 generates a simulated post-perturbation transcriptomic profile for the same patient-drug pair. We align the two views in a shared latent space via a CLIP-style contrastive objective with drug-context hard negatives, then concatenate the representations for end-to-end response classification. On TCGA, PREDIKTOR consistently outperforms state-of-the-art baselines under patient-, drug-, and tissue-split evaluations, and transfers zero-shot to the I-SPY2 trial, improving AUROC by 5.6% over competing methods. The aligned embeddings yield stable gene and pathway attributions that recover known mechanisms, supporting actionable and interpretable precision oncology.
Dongmin Bang, Sugyun An, Inyoung Sung +3
Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, Republic of Korea · AIGENDRUG Co., Ltd., Seoul, Republic of Korea · BK21 FOUR Intelligence Computing, Seoul National University, Seoul, Republic of Korea +2