芒柄花素通过上调GSTM2调控前列腺癌免疫微环境：基于ROS/STAT3介导的趋化因子重编程机制研究

《Frontiers in Immunology》：Integrated pharmacology and validation reveal Formononetin upregulates GSTM2 to modulate prostate cancer immune microenvironment via ROS/STAT3-driven chemokine reprogramming

　　

中文翻译： 背景： 前列腺癌（PCa）在全球男性癌症病例中占相当大比例。生化复发（BCR）是其疾病管理中的关键挑战。芒柄花素（Formononetin）是一种源自中药的天然异黄酮，已显示出潜在的抗癌特性。研究人员旨在确定芒柄花素靶向分子在前列腺癌中的转化潜力。 方法： 通过五个公共数据库鉴定了芒柄花素的推定靶点列表。重叠候选基因通过一系列生物信息学工具进一步筛选，包括差异表达谱分析、共表达网络构建和多种机器学习分类器。对接复合物和补充分析用于评估结合稳定性。为功能验证GSTM2，研究人员在体外评估了其对前列腺癌细胞生长、运动和侵袭的功能作用，以及CXCL10和CCL22趋化因子的表达。 结果： 整合分析确定了413个芒柄花素靶点，其中与335个前列腺癌相关基因重叠。通过加权基因共表达网络分析（WGCNA）和差异表达分析，确定了22个候选基因作为前列腺癌中芒柄花素相关基因。机器学习算法和生存分析进一步提示GSTM2是最重要的靶点。GSTM2与增强的免疫浸润和对比卡鲁胺（Bicalutamide）的敏感性呈正相关。单细胞分析显示，GSTM2阳性恶性细胞表现出上调的谷胱甘肽代谢和抑制的呼吸爆发通路。计算对接证明了芒柄花素与GSTM2之间形成稳定复合物。基于细胞的实验表明，芒柄花素以上调GSTM2并抑制前列腺癌恶性程度，且该作用依赖于GSTM2。在机制方面，芒柄花素通过GSTM2减弱ROS水平并抑制STAT3磷酸化，导致趋化因子表达重编程，其特征为CXCL10上调及CCL22下调。 结论： 芒柄花素通过GSTM2-ROS-STAT3轴重编程趋化因子表达以调节免疫微环境，提示其可能具有免疫调节作用，值得进一步的体内研究。

论文解读文章

一、研究背景与立题依据

前列腺癌是全球男性最常见的恶性肿瘤之一，根治性治疗后十年内生化复发（BCR）发生率可高达46%，但目前尚无专门针对术后这一阶段的获批药物。雄激素剥夺治疗（ADT）虽为高危患者的标准方案，却伴随代谢毒性且不可避免进展为去势抵抗性前列腺癌（CRPC）。中药活性成分因具有多靶点、低毒性和协同作用等特性，在肿瘤防治领域受到广泛关注。芒柄花素（Formononetin）是从黄芪中提取的一种异黄酮类化合物，已知具有抗氧化活性，但其在前列腺癌中的特异性分子靶点和调控网络尚缺乏系统性探索。本研究旨在通过整合药理学与实验验证策略，系统鉴定芒柄花素在前列腺癌中的关键靶点基因及其作用机制。

二、研究内容与主要结论

研究人员通过多层次整合分析，系统鉴定了GSTM2作为芒柄花素调控前列腺癌进展的核心靶点，并阐明了其通过GSTM2-ROS-STAT3轴重编程趋化因子表达、进而调节肿瘤免疫微环境的新机制。该研究不仅将GSTM2定位为芒柄花素的直接下游靶点，还揭示了GSTM2在前列腺癌早期进展中的抑癌功能及其阶段依赖性双重作用，为术后早期干预提供了新的理论依据和治疗策略思路。

三、关键技术方法

研究人员采用了多层次的整合研究策略。首先，通过TCMSP、CTD、HERB、SwissTarget和PharmMapper五个公共数据库预测芒柄花素靶点，结合GeneCards和OMIM数据库获取前列腺癌相关基因。利用TCGA_PRAD、DKFZ_2018和GSE54460三个转录组队列进行差异表达分析和WGCNA共表达网络构建。采用LASSO Cox回归、随机生存森林（RSF）和支持向量机-递归特征消除（SVM-RFE）三种机器学习算法筛选核心基因。通过分子对接和100 ns分子动力学（MD）模拟评估芒柄花素与GSTM2的结合稳定性。利用GSE157703单细胞RNA测序数据解析GSTM2在肿瘤微环境中的细胞类型分布和功能特征。最后通过体外细胞实验（包括LNCaP和VCaP细胞系）验证芒柄花素对GSTM2的调控作用及其功能后果。

四、研究结果

1. 1.
   芒柄花素靶点与前列腺癌进展相关基因的鉴定

通过五个数据库的整合靶点预测共获得413个芒柄花素靶点，其中335个与前列腺癌相关基因重叠。TCGA_PRAD队列的差异表达分析鉴定出3070个差异表达基因（DEGs）。DKFZ_2018队列的WGCNA分析将11647个基因划分为22个共表达模块，其中与BCR相关性|r|>0.3的模块包含1324个基因。三者重叠后最终获得22个候选基因。

1. 2.
   候选靶点的功能图谱

GO富集分析显示22个候选基因主要涉及异生物质代谢、激素调节和细胞周期进程。KEGG通路主要富集于细胞色素P450药物代谢、类固醇激素生物合成和细胞周期通路。蛋白质-蛋白质相互作用（PPI）网络分析表明BCL2、PPARGC1A、KIT和PGR可能作为枢纽基因。

1. 3.
   核心靶点GSTM2的筛选与预后验证

三种机器学习算法（LASSO Cox回归、RSF和SVM-RFE）筛选出两个重叠基因：TK1和GSTM2。单因素Cox回归分析显示GSTM2与BCR风险显著负相关（HR=0.340，95%CI 0.206–0.563，p<0.001）。在DKFZ_2018和GSE54460两个队列中，GSTM2高表达组均表现出显著更优的无BCR生存期，而TK1的预后分层能力较弱且不一致。BCR患者中GSTM2表达显著低于非BCR患者。

1. 4.
   GSTM2的免疫浸润特征与药物敏感性

ESTIMATE分析显示GSTM2高表达与更高的基质评分、免疫评分和ESTIMATE评分相关。ssGSEA分析在三个队列中一致显示GSTM2高表达肿瘤具有更高的固有免疫细胞（肥大细胞、中性粒细胞、NK、γδT）和Th1细胞浸润，同时Tregs减少。药物敏感性分析显示GSTM2表达与比卡鲁胺敏感性呈正相关。

1. 5.
   GSTM2上调减少ROS蓄积

单细胞分析鉴定出七种主要细胞类型，GSTM2表达主要限于上皮（恶性）、内皮和平滑肌细胞。GSEA分析显示GSTM2阳性恶性细胞显著上调谷胱甘肽代谢和异生物质解毒通路，而下调呼吸爆发通路，表明其维持低氧化应激状态。

1. 6.
   芒柄花素与GSTM2结合的计算机分析

分子对接显示芒柄花素与GSTM2的结合自由能为?7.5 kcal/mol，主要通过疏水接触稳定。100 ns MD模拟显示蛋白和配体的RMSD分别在60 ns和95 ns后趋于稳定，Rg和SASA变化有限，RMSF值大多低于0.23 nm，自由能景观（FEL）识别出稳定的低能构象状态。

1. 7.
   芒柄花素以GSTM2依赖方式抑制前列腺癌细胞恶性程度

芒柄花素以剂量依赖方式上调LNCaP和VCaP细胞中GSTM2蛋白水平。CCK-8和集落形成实验表明芒柄花素显著抑制细胞增殖，GSTM2敲低可逆转此抑制作用。Transwell实验证实芒柄花素损害细胞迁移和侵袭能力，且同样依赖GSTM2。

1. 8.
   GSTM2介导芒柄花素诱导的ROS/STAT3驱动的趋化因子重编程

在低剂量H₂O₂处理条件下，芒柄花素显著抑制ROS应答，而GSTM2敲低可逆转此效应。H₂O₂处理选择性升高pSTAT3(Tyr705)而不影响总STAT1、STAT3或pSTAT1(Tyr701)。芒柄花素显著降低pSTAT3水平，GSTM2敲低可逆转此抑制。qPCR和ELISA证实芒柄花素通过GSTM2在转录和分泌水平调控CXCL10上调和CCL22下调，形成有利于Th1/NK而非Treg募集的趋化因子转换。

五、讨论与结论

讨论部分指出，芒柄花素通过直接杀伤肿瘤细胞和重塑免疫微环境两种互补策略发挥抗肿瘤作用。值得注意的是，既往在多发性骨髓瘤中的研究显示芒柄花素通过升高ROS抑制STAT3组成性激活，而本研究发现芒柄花素通过上调GSTM2清除ROS来抑制STAT3磷酸化，这种看似矛盾的调控模式提示芒柄花素对ROS/STAT3轴的调节具有细胞类型依赖性。GSTM2在前列腺癌中表现出阶段依赖性双重功能：在未经治疗的早期前列腺癌中，GSTM2作为看门基因维持抗氧化防御，其表观沉默促进恶性转化；而在CRPC中，肿瘤细胞重新激活GSTM2以抵抗治疗相关氧化应激并获得耐药性。研究局限性包括缺乏siRNA单独处理组的独立表型效应评估、缺乏体内实验验证、免疫微环境重塑需要更充分的验证等。

结论部分：研究人员鉴定GSTM2为芒柄花素在前列腺癌中的直接靶点。GSTM2将芒柄花素介导的氧化还原调节与CXCL10/CCL22重编程联系起来。这种趋化因子信号的重塑有助于创造活跃的免疫微环境。