藤黄根酚类化合物抗乳腺癌潜力揭秘：网络药理学与分子模拟研究

《Journal of Pharmaceutical Innovation》：Unveiling the Potential of Garcinia Atroviridis Root-Derived Phenolic Compounds Against Breast Cancer: Isolation, Network Pharmacology and In-Silico Approaches

本研究针对乳腺癌治疗耐药性难题，通过系统分离鉴定藤黄根中5种酚类化合物，结合网络药理学与分子模拟技术，发现morelloflavone可通过靶向HIF1A信号通路发挥显著抗乳腺癌活性。研究创新性地揭示了藤黄来源酚类化合物的多靶点作用机制，为开发新型乳腺癌治疗药物提供了重要理论依据。

　　

在全球女性癌症发病谱中，乳腺癌始终占据着令人担忧的重要位置。据世界卫生组织数据显示，2022年全球新发乳腺癌病例约230万例，导致约67万人死亡。尽管他莫昔芬等药物已成为雌激素受体阳性乳腺癌的标准治疗方案，但肿瘤转移仍然是乳腺癌患者死亡的主要原因。近年来，癌症干细胞理论的提出为解释治疗耐药性提供了新视角——这群处于静息状态的细胞能够逃避常规治疗，成为肿瘤复发和转移的"罪魁祸首"。

面对日益严峻的耐药性挑战，天然产物研究为抗肿瘤药物开发提供了新的希望。藤黄属植物作为藤黄科中最大的属，包含300多种物种，在传统医学中具有广泛应用历史。其中，藤黄（Garcinia atroviridis）在马来西亚被称为"asam gelugur"，其根茎传统上用于治疗咳嗽、促进血液循环，并在产前产后护理中发挥作用。前期研究表明，该植物叶片精油对MCF-7人乳腺癌细胞具有显著细胞毒性，而根部乙酸乙酯提取物也显示出抗癌潜力。然而，其活性成分的系统分离及其抗乳腺癌分子机制仍有待深入探索。

在这项发表于《Journal of Pharmaceutical Innovation》的研究中，研究人员开展了一项综合性研究，旨在从藤黄根中分离活性成分并评估其抗乳腺癌潜力。研究团队采用色谱技术从藤黄根乙酸乙酯提取物中成功分离得到五种酚类化合物，包括atrovirisidone、1,3,5-三羟基-2-甲氧基口山酮、GB1a、volkensiflavone和morelloflavone，其中volkensiflavone为该物种中首次发现的新化合物。

研究采用网络药理学方法，通过SuperPred和SwissTargetPrediction数据库预测化合物靶点，发现58个共同作用靶点。蛋白互作网络分析识别出HSP90AA1、HIF1A、MAPK1、PIK3R1等14个核心靶点。基因本体论分析显示这些靶点显著富集于癌症相关通路，而京都基因与基因组百科全书通路分析进一步证实其与前列腺癌、癌症通路、病毒致癌等16条信号通路密切相关。

分子对接结果显示，所有五种化合物与低氧诱导因子1α的结合亲和力均优于对照配体，结合能范围为-33.02至-51.19 kcal/mol。其中，morelloflavone表现出最强的结合能力，结合能为-51.19 kcal/mol。该化合物与HIF1A活性位点形成氢键、π-烷基、π-π堆积和静电相互作用，共同稳定了复合物结构。

为验证结合稳定性，研究进行了150纳秒分子动力学模拟。结果显示，morelloflavone-HIF1A复合物具有稳定的均方根偏差值（3.29±0.47 ?），低于3.0 ?的阈值标准，表明复合物结构保持稳定。分子力学/广义波恩表面积分析计算的结合自由能为-53.94±5.14 kcal/mol，显著低于对照复合物（-19.51±4.97 kcal/mol），进一步证实了morelloflavone与HIF1A的强效结合。

ADMET（吸收、分布、代谢、排泄和毒性）特性预测表明，morelloflavone分子量为556.47 g/mol，具有11个氢键受体和7个氢键供体，拓扑极性表面积达198.12 ?²，这些特性可能导致口服生物利用度受限。然而，该化合物不被预测为P-糖蛋白底物，且仅抑制CYP2C9酶，这可能在减少药物相互作用方面具有优势。

主要技术方法

研究采用硅胶柱色谱和制备薄层色谱从藤黄根乙酸乙酯提取物中分离化合物，通过核磁共振技术进行结构鉴定。运用网络药理学方法（SuperPred和SwissTargetPrediction）预测化合物靶点，STRING数据库构建蛋白互作网络，DAVID进行基因本体论和京都基因与基因组百科全书富集分析。采用分子对接评估结合亲和力，进行150纳秒分子动力学模拟验证复合物稳定性，并通过MM/GBSA（分子力学/广义波恩表面积）计算结合自由能。

结构鉴定化合物

从藤黄根乙酸乙酯提取物中成功分离出五种酚类化合物，经鉴定分别为atrovirisidone（7.1 mg）、1,3,5-三羟基-2-甲氧基口山酮（15.4 mg）、GB1a（17.3 mg）、volkensiflavone（17.4 mg）和morelloflavone（15 mg）。核磁共振分析确认了化合物结构和纯度，其中volkensiflavone为首次从该物种中报道。

靶点基因鉴定

通过SuperPred和SwissTargetPrediction工具分别预测到191和322个靶基因，Venn图分析鉴定出58个共同靶基因。这些靶点包括PLAU、NOS2、CDK5、ITGB1、HIF1A等关键分子，其中15个靶基因在癌症相关通路中显著富集。

构建预测基因的蛋白互作网络

蛋白互作网络分析显示54个靶蛋白之间存在185个相互作用边，平均节点度为6.38。通过网络拓扑分析识别出14个核心枢纽基因，包括HSP90AA1、HIF1A、MAPK1、PIK3R1等，这些基因在乳腺癌发生发展中发挥关键作用。

靶点蛋白的基因本体论和京都基因与基因组百科全书通路富集分析

基因本体论分析显示靶基因显著富集于13个生物过程、15个细胞组分和19个分子功能类别。京都基因与基因组百科全书通路分析鉴定出16条显著富集通路，主要包括癌症通路、前列腺癌、HIF-1信号通路等，证实了靶基因在癌症生物学中的重要作用。

分子对接研究

分子对接结果显示，五种酚类化合物与HIF1A的结合能（-33.02至-51.19 kcal/mol）均优于对照配体（-24.55 kcal/mol）。其中morelloflavone结合能最低（-51.19 kcal/mol），与活性位点残基形成多种相互作用，表明其具有最强的结合亲和力。

分子动力学模拟和MM/GBSA

150纳秒分子动力学模拟表明，morelloflavone-HIF1A复合物结构稳定，均方根偏差值为3.29±0.47 ?，配体均方根偏差为1.14±0.36 ?。MM/GBSA计算结合自由能为-53.94±5.14 kcal/mol，显著低于对照复合物，证实了morelloflavone与HIF1A的强效结合。

ADMET研究

瑞士ADME预测显示，morelloflavone具有高极性表面积和氢键数量，可能限制其口服生物利用度，但非P-糖蛋白底物特性及选择性CYP450抑制特征可能减少药物相互作用风险。

本研究通过系统性的分离鉴定、网络药理学分析和分子模拟，揭示了藤黄根来源酚类化合物特别是morelloflavone的抗乳腺癌潜力。研究发现这些化合物通过多靶点作用机制影响癌症相关通路，其中对HIF1A的强效抑制可能是其发挥抗癌活性的关键机制。尽管morelloflavone在药代动力学特性方面存在一定局限性，但其独特的靶点选择性和作用机制为开发新型乳腺癌治疗药物提供了有前景的先导化合物。这些发现不仅深化了对藤黄植物药理活性的理解，也为天然产物抗肿瘤药物研发提供了新的思路和实验依据。