Exploring the toxicological impact of DEHP exposure on colorectal cancer through network toxicology, machine learning and bioinformatics analysis

OBJECTIVE: Colorectal cancer (CRC) ranks among the most prevalent malignant tumors, yet its underlying mechanisms remain not fully understood. The role of environmental factors is critical in its development and progression. Di-(2-ethylhexyl) phthalate (DEHP), a widespread environmental contaminant, poses significant hazards to human health. Prior research indicates that exposure to DEHP can disrupt cellular processes like proliferation, differentiation, and apoptosis, which may contribute to cancer development. Nevertheless, the specific role and molecular pathways involved with DEHP in CRC are still poorly defined. METHODS: We integrated multiple public databases to identify overlapping targets of DEHP and CRC. Machine learning methods were applied to prioritize core targets, whose expression levels and diagnostic performance were then validated in the Gene Expression Omnibus (GEO) datasets. The association between core targets and immune cell infiltration was evaluated using CIBERSORT-based deconvolution. Finally, molecular docking was performed to predict the interactions between DEHP and the hub proteins, and, based on the docking results, an explicit-solvent molecular dynamics (MD) simulation of the top-ranked DEHP–target complex was conducted to further assess binding stability. RESULTS: Analysis of publicly available databases revealed 57 potential targets that were significantly enriched in signaling pathways related to “microRNAs in cancer” and “chemical carcinogenesis–receptor activation.” Machine learning methods identified five primary targets: CASP3, BCL6, BRD4, PPARA, and PRKCD. The expression levels of these genes in CRC tissues were significantly different from those in control tissues and showed reasonable diagnostic performance, as well as correlations with immune cell infiltration. Molecular docking suggested that DEHP can bind all five hub proteins, with PPARA exhibiting the most favorable binding affinity (− 6.7 kcal/mol). The 100-ns MD simulation further supported this binding mode by demonstrating that the DEHP–PPARA complex maintains a dynamically stable conformation with persistent key interactions. CONCLUSIONS: Our findings generate preliminary mechanistic hypotheses on how DEHP exposure may contribute to CRC by perturbing known cancer-related pathways and reshaping the intestinal immune microenvironment, with PPARA emerging as a key candidate mediator. The study highlights the utility of combining network toxicology, machine learning, immune deconvolution, molecular docking, and MD simulations to evaluate the potential carcinogenic risks of environmental pollutants.