综合网络药理学与分子动力学揭示美洲鳄梨中的木犀草素作为多癌SRC/GSK3β抑制剂

《International Journal of Molecular Sciences》：Integrated Network Pharmacology and Molecular Dynamics Reveal Luteolin from Persea americana as a Multi-Cancer SRC/GSK3β Inhibitor

　　

癌症进展由多种恶性肿瘤中失调的激酶信号传导和凋亡逃避驱动。尽管靶向激酶抑制剂改善了预后，但耐药性和毒性仍然是主要挑战。天然植物化学物质提供了有希望的多靶点治疗潜力。美洲鳄梨（Persea americana）含有多种生物活性化合物；然而，它在多癌激酶靶向中的作用仍未充分探索。本研究旨在使用综合网络药理学和分子模拟方法，识别并验证美洲鳄梨植物成分在五种癌症（肺癌、乳腺癌、宫颈癌、结直肠癌和前列腺癌）中的抗癌激酶靶点。癌症相关基因从Open Targets Platform获取，并通过基因本体论（Gene Ontology，GO）分析进行优先排序。识别了与208个预测的美洲鳄梨人类靶点的重叠靶点。蛋白质-蛋白质相互作用网络揭示了枢纽基因，随后进行了基于癌症基因组图谱（The Cancer Genome Atlas，TCGA）的验证。25种植物成分对SRC（一种非受体酪氨酸激酶）和糖原合成酶激酶3β（GSK3β）进行了分子对接，顶级复合物进行了100纳秒分子动力学模拟。富集分析突出了激酶活性和凋亡。SRC成为泛癌枢纽，而GSK3β在乳腺癌中突出。木犀草素（Luteolin）对SRC显示出最强的结合（?11.9 kcal/mol），优于共晶抑制剂，而瓦伦烯（Valencene）对GSK3β显示出亲和力（?8.8 kcal/mol）。模拟证实了稳定的相互作用。木犀草素表现出强大的多靶点激酶抑制，特别是针对SRC，支持其作为泛癌治疗候选药物的潜力。

**论文解读：基于网络药理学与分子动力学揭示美洲鳄梨中木犀草素作为多癌SRC/GSK3β抑制剂**

**研究背景与问题**

癌症是全球发病率和死亡率的主要原因之一，2022年全球新发病例近2000万，死亡970万。肺癌、乳腺癌、结直肠癌、前列腺癌和宫颈癌这五种主要癌症占全球癌症发病率的近40%，且在不同收入国家中均构成重大负担。尽管靶向激酶抑制剂在临床治疗中取得了一定进展，但耐药性和毒副作用仍是主要挑战。天然产物因其多靶点调节特性而备受关注，其中美洲鳄梨（Persea americana）含有多种生物活性化合物，如黄酮类、萜类等，具有抗炎、抗氧化和抗癌活性。然而，其植物成分在多癌激酶靶向中的系统研究尚不充分。为此，研究人员开展本研究，旨在通过整合网络药理学和分子模拟方法，系统评估美洲鳄梨植物成分对五种主要癌症（肺癌、乳腺癌、宫颈癌、结直肠癌、前列腺癌）的抗癌潜力，并识别关键激酶靶点。

**研究方法概述**

本研究采用的主要关键技术方法包括：从Open Targets Platform获取五种癌症的癌症相关基因数据集，并通过基因本体论（GO）富集分析进行功能优先排序；将优先排序后的基因与美洲鳄梨的208个预测人类靶点进行交集分析；利用STRING数据库构建蛋白质-蛋白质相互作用（PPI）网络，并通过Cytoscape的CytoNCA插件识别枢纽基因；基于TCGA数据库进行基因表达验证，以SRC和GSK3β作为最终靶点；使用AutoDock Vina进行分子对接，并对最优复合物进行100纳秒分子动力学模拟（Desmond模块，OPLS3e力场）。样本队列来源于TCGA数据库。

**研究结果**

**2.1 癌症相关基因数据集**：通过检索Open Targets Platform，获得各癌症相关基因数量：肺癌15702个，乳腺癌17229个，宫颈癌8452个，结直肠癌15057个，前列腺癌11229个。乳腺癌基因数量最多，反映了其分子异质性，而宫颈癌基因数量最少，提示其病因相对明确。

**2.2 功能富集与靶点优先排序**：GO富集分析显示，在所有五种癌症中，蛋白激酶活性、凋亡过程调控和细胞表面组分均显著富集（FDR值极低，范围为10^?36至10^?164），表明这些癌症共享由失调的激酶信号、凋亡逃逸和膜相关靶点构成的致癌架构。

**2.3 美洲鳄梨与癌症靶点的共同基因识别**：通过维恩图分析，美洲鳄梨靶点与各癌症优先基因的交集数量为：肺癌100个，乳腺癌98个，宫颈癌86个，结直肠癌98个，前列腺癌96个。约3-4%的重叠率支持选择性富集，提示美洲鳄梨可能通过共同调控机制发挥多癌治疗作用。

**2.4 蛋白质-相互作用网络分析与枢纽基因识别**：PPI网络分析显示，各癌症网络均高度互联。EGFR、SRC、ESR1、HSP90AB1和GSK3β在所有五种癌症中均位列前五枢纽基因，其中MMP9在宫颈癌中特异性出现。表明这些枢纽基因是关键的致癌信号调节因子。

**2.5 基因表达验证与最终靶点优先排序**：基于TCGA的表达验证表明，SRC在肺癌、宫颈癌、结直肠癌和前列腺癌中表现出最广泛的表达变化和上四分位数，而在乳腺癌中，GSK3β显示出更强的表达优势。因此，SRC被确定为泛癌靶点，GSK3β作为乳腺癌特异性靶点。

**2.6 活性位点识别**：通过PLIP分析，SRC（PDB ID: 3F3V）的ATP结合口袋包含Leu273、Glu310、Met341、Asp404等关键残基；GSK3β（PDB ID: 7SXJ）的活性位点包括Ile62、Val70、Ala83、Tyr134、Val135等。

**2.7 分子对接结果**：25种植物成分对接显示，木犀草素（Luteolin）对SRC结合亲和力最强（?11.9 kcal/mol），优于共晶抑制剂（?9.0 kcal/mol）；瓦伦烯（Valencene）对GSK3β结合最强（?8.8 kcal/mol），优于共晶配体（?7.7 kcal/mol）。这些结果支持木犀草素的高亲和力多靶点抑制潜力。

**2.8 分子动力学模拟分析**：100 ns模拟显示，SRC-木犀草素复合物的蛋白骨架RMSD稳定在2.2-3.0 ?，配体RMSD低于2.1 ?，远优于SRC-瓦伦烯（配体RMSD高达10 ?）和SRC-共晶复合物。GSK3β-木犀草素复合物同样表现出最低的RMSD和RMSF，配体RMSF在0.7-0.9 ?，活性位点残基波动小。蛋白-配体接触分析进一步证实，木犀草素在两种激酶中均维持了关键的铰链区氢键和催化残基相互作用，稳定性和持久性接近或优于共晶对照。

**讨论与结论**

讨论部分总结了本研究的系统工作流程：从大规模癌症靶点挖掘、功能优先排序、交集分析、网络拓扑评估、转录组验证，到分子对接和长时间尺度分子动力学模拟，最终将数千个癌症相关基因缩小至两个机制上关键的激酶靶点SRC和GSK3β。功能富集一致显示五种癌症在蛋白激酶活性、凋亡控制和膜信号传导上的汇聚，表明共享的致癌架构。交集分析显示选择性且可重复的重叠，支持美洲鳄梨的多癌相关性。网络中心性分析识别了SRC和GSK3β作为主导调控枢纽，TCGA表达验证确认SRC为泛癌失调靶点，GSK3β在乳腺癌中突出。结构对接和分子动力学模拟验证了木犀草素和瓦伦烯的高亲和力结合，其中木犀草素在两种激酶系统中均表现出优越的结构稳定性、持久的铰链区氢键和稳定的催化位点占据。

研究结论部分的翻译如下：在本研究中，研究人员采用整合的系统药理学和分子模拟框架，阐明了美洲鳄梨植物成分的多癌治疗潜力。通过结合大规模癌症靶点挖掘、功能基因优先排序、交集分析、网络拓扑评估、转录组验证、分子对接和长时间尺度分子动力学模拟，研究人员系统地将数千个癌症相关基因缩小至两个机制上中心的激酶靶点——SRC和GSK3β。五种主要恶性肿瘤的功能富集一致显示在失调的蛋白激酶活性、凋亡控制通路和膜相关信号上的汇聚，强调了共享的致癌架构。交集分析显示癌症相关基因与预测的美洲鳄梨靶点之间选择性和可重复的重叠，支持有生物学意义的多癌相关性。网络中心性分析识别出SRC和GSK3β作为主导调控枢纽，基于TCGA的基因表达验证进一步确认SRC是泛癌失调靶点，而GSK3β在乳腺癌中表现出癌症特异性突出。基于结构的对接表明，木犀草素和瓦伦烯在SRC和GSK3β的ATP结合催化口袋中表现出与共晶抑制剂相当或更优的结合亲和力。重要的是，100纳秒的分子动力学模拟提供了机制上的验证，揭示木犀草素在两种激酶系统中均维持了优越的结构稳定性、持久的铰链区氢键、最小的骨架扰动和持续的催化位点结合。蛋白-配体接触分析进一步确认了与关键催化残基的保守相互作用，增强了对接预测的可靠性。总体而言，这些发现表明，美洲鳄梨中的木犀草素是一种计算预测的、动态稳定的多靶点激酶抑制剂，具有潜在的泛癌适用性。网络中心性、转录组失调、高结合亲和力和长期构象稳定性的汇聚为支持其潜在治疗相关性提供了强有力的计算证据。然而，这些发现需要通过体外和体内研究进行进一步的实验验证。这项工作推进了对鳄梨来源植物化学物质在肿瘤学中的理解，并展示了一个稳健的整合计算框架，用于合理的多癌药物发现。所识别的植物成分，特别是木犀草素，代表了基于计算预测的激酶靶向癌症治疗中的潜在先导化合物。尽管如此，这些发现应被视为假设生成性，需要全面的体外和体内实验验证，才能确认其治疗潜力。