# Systematic Investigation of Tumor Immune Microenvironment Modulation by Cynomorium songaricum Against Breast Cancer Through Integrated Chemomics, Network Pharmacology and Molecular Docking

**Authors:** Ze-An Mao, Mei-Ling Zhang, Zi-Yi An, Wei-Lin Jin

PMC · DOI: 10.3390/ph19020314 · 2026-02-13

## TL;DR

This study explores how Cynomorium songaricum, a traditional Chinese medicine, modulates the tumor immune environment in breast cancer using chemical profiling and computational methods.

## Contribution

The study systematically identifies bioactive compounds in Cynomorium songaricum and their immunomodulatory mechanisms against breast cancer.

## Key findings

- 1100 compounds were identified, with 84 meeting drug-likeness criteria, including flavonoids and terpenes.
- Network pharmacology revealed 776 overlapping targets linked to breast cancer and immune regulation.
- Molecular docking and MD simulations confirmed stable interactions between compounds and immune-related targets.

## Abstract

Background/Objectives: Breast cancer remains a leading cause of cancer-related mortality in women, with therapeutic resistance frequently arising from tumor heterogeneity and an immunosuppressive tumor immune microenvironment (TIME). While Cynomorium songaricum Rupr. (CS) has been used traditionally in Chinese medicine and exhibits preliminary anti-tumor activity, its bioactive constituents and precise mechanisms against breast cancer remain to be elucidated. Methods: The chemical constituents of CS were systematically profiled using ultra-high-performance liquid chromatography coupled with Q Exactive Orbitrap mass spectrometry (UHPLC-Q-Exactive Orbitrap MS/MS). Network pharmacology and functional enrichment analyses were performed to identify immuno-related targets and pathways, followed by molecular docking to prioritize component–target pairs. Molecular dynamics (MD) simulations were conducted to validate the stability of a representative docked complex and to characterize binding stability, interaction persistence, molecular mechanics/(Poisson–Boltzmann) surface area (MM/(P)BSA) energetics, and principal component analysis (PCA)-based conformational landscapes. Results: We identified 1100 compounds, of which 84 satisfied the in silico drug-likeness criteria, including 12 phenylpropanoids, 4 terpenes, 35 flavonoids, 2 quinones, 1 phenol, 3 alkaloids, and other phytochemicals. Network pharmacology analysis revealed 776 overlapping targets associated with both breast cancer and immune regulation. Functional enrichment analysis underscored significant involvement in immune-related pathways, and molecular docking studies supported high-affinity interactions between the components and their targets. MD analyses further supported a stable bound ensemble for the representative SRC–Tomentogenin complex during the equilibrated window, with persistent pocket occupancy, consistent interaction signatures, favorable MM/(P)BSA binding energetics, and a concentrated low-energy basin on the PCA-based free energy landscape. Conclusions: These findings elucidate the chemical basis of CS and uncover its immunomodulatory mechanism against breast cancer, offering a foundation for developing CS-based immunotherapeutic strategies and supporting multi-target drug discovery from traditional medicines.

## Linked entities

- **Proteins:** SRC (SRC proto-oncogene, non-receptor tyrosine kinase)
- **Chemicals:** phenol (PubChem CID 996)
- **Diseases:** breast cancer (MONDO:0004989)

## Full-text entities

- **Genes:** PIK3CB (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta) [NCBI Gene 5291] {aka P110BETA, PI3K, PI3KBETA, PIK3C1}, PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha) [NCBI Gene 5290] {aka CCM4, CLAPO, CLOVE, CWS5, HMH, MCAP}, PCSK1 (proprotein convertase subtilisin/kexin type 1) [NCBI Gene 5122] {aka BMIQ12, NEC1, PC1, PC1/3, PC3, SPC3}, MTOR (mechanistic target of rapamycin kinase) [NCBI Gene 2475] {aka FRAP, FRAP1, FRAP2, RAFT1, RAPT1, SKS}, TGFB1 (transforming growth factor beta 1) [NCBI Gene 7040] {aka CAEND1, CED, DPD1, IBDIMDE, LAP, TGF-beta1}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, ERBB2 (erb-b2 receptor tyrosine kinase 2) [NCBI Gene 2064] {aka CD340, HER-2, HER-2/neu, HER2, MLN 19, MLN-19}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, CS (citrate synthase) [NCBI Gene 1431], KRT6B (keratin 6B) [NCBI Gene 3854] {aka CK-6B, CK6B, K6B, KRTL1, PC2, PC4}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, AKT1 (AKT serine/threonine kinase 1) [NCBI Gene 207] {aka AKT, PKB, PKB-ALPHA, PRKBA, RAC, RAC-ALPHA}, SRC (SRC proto-oncogene, non-receptor tyrosine kinase) [NCBI Gene 6714] {aka ASV, SRC1, THC6, c-SRC, p60-Src}, STAT3 (signal transducer and activator of transcription 3) [NCBI Gene 6774] {aka ADMIO, ADMIO1, APRF, HIES}, IL10 (interleukin 10) [NCBI Gene 3586] {aka CSIF, GVHDS, IL-10, IL10A, TGIF}, NR4A1 (nuclear receptor subfamily 4 group A member 1) [NCBI Gene 3164] {aka GFRP1, HMR, N10, NAK-1, NGFIB, NP10}
- **Diseases:** Breast Cancer (MESH:D001943), cytotoxicity (MESH:D064420), oncogenesis (MESH:D063646), injury to (MESH:D014947), inflammatory (MESH:D007249), MD (MESH:D000092242), Cancer (MESH:D009369)
- **Chemicals:** steroid (MESH:D013256), Andropanolide (MESH:C510040), Hydrogen (MESH:D006859), Lys (MESH:D008239), ice (MESH:D007053), FA (MESH:D005492), Glycosides (MESH:D006027), Flavonoids (MESH:D005419), CS (-), Aurantiamide acetate (MESH:C011670), Amino Acids (MESH:D000596), Asp (MESH:D001224), Carbohydrates (MESH:D002241), Arg (MESH:D001120), Phenols (MESH:D010636), Phe (MESH:D010649), Sesamolin (MESH:C054124), phenol (MESH:D019800), water (MESH:D014867), Terpenoids (MESH:D013729), Peptides (MESH:D010455), paeoniflorin (MESH:C015423), phenolic acids (MESH:C017616), liquiritigenin (MESH:C083152), Alkaloids (MESH:D000470), Asn (MESH:D001216), methanol (MESH:D000432), NaCl (MESH:D012965), quinones (MESH:D011809), oxygen (MESH:D010100), formic acid (MESH:C030544), nitrogen (MESH:D009584), polysaccharides (MESH:D011134), acetonitrile (MESH:C032159)
- **Species:** Bacillus sp. SA (species) [taxon 1168094], Cynomorium songaricum (species) [taxon 627609], Homo sapiens (human, species) [taxon 9606]

## Figures

8 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12944387/full.md

---
Source: https://tomesphere.com/paper/PMC12944387