《Beni-Suef University Journal of Basic and Applied Sciences》：Network pharmacology and computational dissection of Solanum trilobatum bioactives reveals modulation of steroid receptor associated signalling in breast cancer
编辑推荐：
乳腺癌是女性最常见的恶性肿瘤，现有疗法副作用明显。本研究采用网络药理学、分子对接及分子动力学模拟等计算生物学技术，探究药用植物三裂茄（Solanum trilobatum）的抗癌潜力。研究筛选出16种具有抗乳腺肿瘤活性的先导化合物，明确了PI3K-Akt为关键调控通路，并揭示了其靶向KRAS、TP53、IGF1R等关键靶点的稳定结合模式，为开发新型植物源抗乳腺癌药物提供了重要的理论依据。
在女性健康领域，乳腺癌犹如一个挥之不去的阴霾，是全球女性发病率和死亡率最高的恶性肿瘤之一。尽管手术、放疗、化疗、靶向治疗和免疫治疗等现代医学手段不断发展，但现有的治疗方案，尤其是化疗，常常伴随着难以忽视的副作用。例如，常用的芳香化酶抑制剂可能导致潮热、抑郁和记忆障碍；选择性雌激素受体调节剂（SERM）有引发血栓的风险；而阿霉素等蒽环类药物则与心脏毒性相关。这些副作用不仅影响患者的生活质量，有时甚至迫使治疗中断。因此，从大自然中寻找高效、低毒的替代或辅助治疗药物，成为了一个极具吸引力的研究方向。
药用植物三裂茄（Solanum trilobatum）在传统医学中应用已久，现代研究也已证实其具有抗氧化、抗炎、抗菌和抗肿瘤等多种药理活性。特别是，利用其果提取物制备的银纳米颗粒对乳腺癌MCF-7细胞显示出细胞毒性。然而，一个关键的科学问题尚未得到充分解答：三裂茄中究竟是哪些具体的化学成分在发挥作用？它们又是通过怎样的分子“暗语”与人体内的癌症相关蛋白“对话”，从而抑制肿瘤生长的？其抗乳腺癌的分子作用机制仍然是一个“黑箱”。
为了破解这个“黑箱”，揭开三裂茄抗癌潜力的神秘面纱，一个研究团队在《Beni-Suef University Journal of Basic and Applied Sciences》上发表了一项综合性研究。他们没有直接从实验室的细胞培养皿或动物模型开始，而是首先走进了强大的“计算世界”，采用了一套整合的计算系统生物学策略。他们的目标很明确：通过虚拟筛选和模拟，系统性地阐明三裂茄生物活性成分对抗乳腺癌的分子机制，为后续的实验验证和药物开发铺平道路。
研究人员主要运用了网络药理学、分子对接和分子动力学模拟三大关键技术。网络药理学如同绘制一张“作战地图”，他们首先从多个数据库挖掘出与三裂茄及其同义词相关的425个植物化学成分。然后，通过药物相似性评分、类药五原则和ADMET（吸收、分布、代谢、排泄和毒性）性质预测，并结合预测的抗肿瘤（乳腺癌）活性，层层筛选，最终锁定了16个最有潜力的先导化合物。接着，他们“顺藤摸瓜”，找出这些化合物可能作用的蛋白质靶点，并与已知的2141个乳腺癌相关靶点取交集，构建了蛋白质-蛋白质相互作用网络，并通过KEGG通路富集和基因本体（Gene Ontology, GO）分析，识别出核心的信号通路和枢纽基因。分子对接技术则像是“锁钥匹配”实验的计算机版本，研究人员使用AutoDock Vina和Schr?dinger Glide软件，让筛选出的生物活性成分与预测的关键致癌靶点蛋白进行“虚拟结合”，计算结合能，评估结合模式的稳定性。最后，分子动力学模拟担任了“稳定性测试员”的角色，对结合能最优的蛋白-配体复合物进行了长达200-500纳秒的模拟，观察在实际溶液环境中复合物的结构动态、稳定性和结合强度，并通过MMPBSA（分子力学泊松-玻尔兹曼表面积）方法计算了精确的结合自由能。
3.1 网络药理学结果
3.1.1 植物化学成分的挖掘与筛选：从三裂茄及其8个同义词中，共鉴定出425个植物化学成分。经过基于类药性、ADMET性质和预测的抗肿瘤活性的多轮筛选，最终有16个化合物脱颖而出，进入后续研究。这些化合物包括瑞奥菊素（rheochrysin）、Hypoxiside A、细胞松弛素Z4（cytochalasin Z4）、4-对香豆酰奎宁酸（4-p-coumaroylquinic acid）等。
3.1.2 靶点识别与网络分析：研究发现，16个生物活性成分共同调节531个蛋白质靶点，其中与乳腺癌相关的有176条信号通路被显著富集。最重要的发现是，PI3K-Akt信号通路被确定为受影响最显著的核心通路。网络拓扑分析进一步识别出10个关键的枢纽基因（hub genes），包括AKT1、TP53、TNF、STAT3、MYC、IL6、HIF1A、EGFR、CTNNB1和CASP3。基因本体分析显示，这些靶点主要富集在“受体复合物”、“蛋白结合”和“对化学物质的反应”等生物学过程中，表明三裂茄成分可能通过干扰细胞表面受体信号和应激反应来发挥抗癌作用。约90.8%的KEGG通路基因与GO分析基因重叠，证明了靶点功能与通路关联的高度一致性。
3.2 分子对接结果
分子对接模拟揭示了多个高亲和力的结合对。使用AutoDock Vina模拟发现，瑞奥菊素与HRAS、Hypoxiside A与KRAS、细胞松弛素Z4与NCOA3、瑞奥菊素与PGR、以及4-对香豆酰奎宁酸与PIK3CB的结合能最优，表明它们能稳定地结合在这些致癌靶点的活性位点。进一步的Schr?dinger Glide精确对接显示，4-对香豆酰奎宁酸与KRAS、以及Hypoxiside A与TP53的结合最为稳固，对接得分最低。这些相互作用主要由氢键和π-π堆叠等分子间作用力维持。
3.3 分子动力学模拟与MMPBSA分析结果
为了验证对接结果的可靠性，研究人员对三个最有希望的复合物进行了长期的分子动力学模拟：KRAS–4-对香豆酰奎宁酸、IGF1R–Kwakhurin和TP53–Hypoxiside A。
- •
KRAS复合物：在200纳秒模拟中，复合物表现出良好的结构稳定性，根均方偏差（RMSD）和回旋半径（Rg）波动在合理范围内。能量分解显示，氨基酸PHE28、LYS147等对结合有主要贡献。
- •
IGF1R复合物：该复合物在模拟早期经历构象调整后趋于稳定。MMPBSA计算得到的结合自由能为+53.43 kcal/mol，提示Kwakhurin与IGF1R的结合在热力学上可能不利或较弱。
- •
TP53复合物：该复合物达到平衡所需时间较长，模拟延长至500纳秒。其RMSD波动略大，但整体结构保持完整。MMPBSA计算得到的结合自由能为+2.63 kcal/mol，也表明Hypoxiside A与TP53的结合亲和力一般。
值得注意的是，KRAS复合物的MMPBSA结合自由能为-3.82 kcal/mol，表明4-对香豆酰奎宁酸与KRAS的结合是自发且较稳定的。
4 讨论与结论
本研究通过整合计算生物学方法，系统揭示了药用植物三裂茄抗乳腺癌的多靶点、多通路作用机制。首先，网络药理学分析指出，其活性成分可能通过协同调控PI3K-Akt这一核心信号轴，以及MAPK、雌激素、Wnt、mTOR等多个与肿瘤增殖、存活、转移和耐药密切相关的通路来发挥作用。所识别的10个枢纽基因均在乳腺癌发生发展中扮演关键角色，如AKT1促进细胞存活，TP53调控基因组稳定和凋亡，EGFR驱动细胞增殖等。
其次，研究特别强调了三裂茄成分对类固醇激素受体相关信号网络的潜在调控作用。富集到的雌激素信号通路、内分泌抵抗通路等均与雌激素、孕激素受体活性密切相关。分子对接发现的与核受体共激活因子NCOA3和孕激素受体PGR的强结合，进一步支持了三裂茄可能通过调节激素受体介导的转录程序来影响乳腺癌进程。
最后，分子对接和动力学模拟为作用机制提供了原子层面的见解。尽管部分复合物的MMPBSA结合自由能计算值未显示出极强的结合力，但稳定的模拟轨迹和特定的残基相互作用表明，如4-对香豆酰奎宁酸与KRAS、Hypoxiside A与TP53等组合，确实可能形成具有一定稳定性的复合物，从而干扰这些关键癌蛋白或抑癌蛋白的功能。例如，靶向KRAS一直是药物研发的难点，本研究发现天然化合物能与之稳定结合，具有重要的启发性。
综上所述，该研究首次通过计算系统生物学方法，全面描绘了三裂茄抗乳腺癌的“作用地图”，筛选出多个有潜力的先导化合物，并明确了其通过干扰PI3K-Akt等信号通路、靶向KRAS、TP53、NCOA3等关键节点来发挥作用的分子机制。这些发现不仅为理解三裂茄的传统药用价值提供了现代科学依据，更重要的是，为后续针对乳腺癌，特别是激素受体相关型乳腺癌的新型植物源药物的设计与开发，奠定了坚实的理论和数据基础。
