Integrated Systems Toxicology Analysis and Molecular Dynamics Study Reveal the Potential Mechanism by which Benzo[a]pyrene Contributes to Depression in Lung Adenocarcinoma via RNASE1 Upregulation

Abstract

Benzo[a]pyrene (Bap), a prominent carcinogen found in tobacco smoke and air pollution, is closely associated with lung adenocarcinoma (LUAD) and subsequent depression, though the mechanism underlying this comorbidity remains unclear. LUAD patients were categorized based on depression-related gene expression. Prognostic genes identified through machine learning were cross-referenced with predicted Bap targets. Mendelian randomization (MR) evaluated the causal relationship between RNASE1 and depression. Single-cell RNA sequencing (scRNA-seq), molecular docking, and 100 ns molecular dynamics (MD) simulations (using the CHARMM36 force field and TIP3P water model) were conducted, followed by analyses of RMSD, RMSF, Rg, SASA, and Gibbs free energy landscape. A Bap-exposed mouse model was used to validate key gene expression and depression-like behaviors in vivo. Consensus clustering based on 52 depression-related genes stratified the TCGA-LUAD cohort into two gene signature-based molecular subtypes with distinct expression patterns and prognostic outcomes. Cluster 1 represented a favorable-prognosis subtype with better disease-free survival, whereas Cluster 2 represented a poor-prognosis subtype with stronger depression-related transcriptomic features; however, these clusters should not be interpreted as clinically confirmed depression subgroups because TCGA lacks psychiatric assessment data. Among 57 depression-related prognostic genes intersected with 375 Bap targets, five genes overlapped: CLEC7A, NEK2, RNASE1, ITGAL, and CPA3. MR confirmed that elevated lung RNASE1 increased depression risk (IVW: OR = 1.008, 95% CI: 1.002–1.014, P = 0.005), while CLEC7A and other genes showed no causal link. scRNA-seq localized RNASE1 upregulation to tumor epithelial cells, monocytes, and macrophages. Docking and MD analyses supported stable Bap–RNASE1 binding, involving hydrophobic interactions with key residues and a dominant low-energy conformational state. In vivo, Bap-exposed mice showed prolonged immobility time in the tail suspension test, fewer open-arm entries, and reduced open-arm time in the elevated plus-maze test, accompanied by elevated WBC, MPV, PDW, and increased lung RNASE1 protein expression. Bap may contribute to LUAD-related depression through RNASE1 upregulation. This provides molecular-level mechanistic insight into the cancer–depression comorbidity driven by environmental carcinogens.