Integrative network pharmacology and experimental validation reveal emodin derivatives as potential therapeutics for hepatocellular carcinoma.

Background: Hepatocellular carcinoma (HCC) is a major global health concern due to its high prevalence and mortality rate. Although emodin, a natural anthraquinone derivative, has demonstrated in vitro anticancer activity against HCC cells, its specific molecular targets in HCC remain unclear.

Method: This study used an integrated approach combining in silico network pharmacology, molecular docking, molecular dynamics simulations (MDS), and in vitro cytotoxicity assays to evaluate three emodin derivatives: emodin, 3-acetyl emodin (ACE), and 1,3,8-triacetyl emodin (TAEM). Target predictions were performed using the SwissTargetPrediction database, and HCC-related genes were retrieved from cBioPortal. Functional annotations (Gene Ontology and Reactome) identified EGFR and KIT as key targets. Docking simulations were conducted to assess binding affinities, followed by 100 ns MDS to evaluate stability. Cytotoxic effects on HepG2 cells were also assessed.

Result: TAEM showed the strongest binding affinity to both EGFR and KIT and demonstrated the highest cytotoxicity against HepG2 cells (IC50 = 0.021 mM). MDS results indicated that the KIT-TAEM complex was the most stable among all tested combinations, supported by RMSD, RMSF, Rg, protein-ligand distance, and MM-GBSA binding energy analyses.

Conclusion: These findings highlight TAEM as a promising therapeutic candidate for HCC. The study demonstrates the value of integrating computational predictions with experimental validation in early-stage drug discovery.