Unraveling the molecular mechanisms of irisin in liver cancer through network toxicology and molecular docking strategies

技术摘要：揭示鸢尾素（Irisin）在肝癌中的分子机制

问题陈述
肝细胞癌（HCC）仍是全球范围内导致癌症相关死亡的主要原因，且往往在诊断时已处于晚期，导致治愈选择有限。虽然运动已知能降低疾病风险并改善生活质量，但运动诱导的肌因子（特别是鸢尾素）调节 HCC 发病机制的具体分子机制仍缺乏充分表征。现有文献对鸢尾素在 HCC 中的作用存在矛盾证据：一些研究表明它通过 PI3K/Akt 通路增强增殖和迁移，而另一些研究则报告其无影响甚至在肿瘤组织中表达下调。此外，关于 HCC 组织与正常对照组中鸢尾素表达水平是否存在共识也尚未达成。本研究旨在通过系统生物学方法，全面识别鸢尾素在 HCC 中的潜在分子靶点及信号通路，以解决这些不一致性。

研究方法
本研究采用了整合网络毒理学和分子对接的多步生物信息学流程：

1. 靶点收集：

   • 鸢尾素靶点： 从 SwissTargetPrediction、ChEMBL 和 STITCH 数据库中检索受鸢尾素调控的潜在人类靶基因。
   • HCC 靶点： 从 GeneCards 和 OMIM 中提取与肝癌相关的基因。
   • 差异表达： 使用 DESeq2 对来自 TCGA-LIHC 队列（374 例肿瘤样本和 50 例正常样本）的 RNA-seq 数据进行分析，以识别差异表达基因（DEGs）（截断值：adj. P < 0.05, |log₂FC| > 1）。
   • 交集： 鸢尾素靶点、HCC 相关基因以及 TCGA 衍生 DEGs 的交集产生了 93 个潜在靶点。

2. 网络构建与核心基因识别：

   • 利用 STRING 数据库构建蛋白质-蛋白质相互作用（PPI）网络（置信度分数 > 0.4），并使用 Cytoscape 进行可视化。
   • 模块分析： 使用 MCODE 插件识别最关键的子网络。
   • 核心基因选择： 使用 cytoHubba 插件中的五种算法（度、介数、紧密性、MNC、径向度）对基因进行排名。这五种算法中排名靠前的基因的交集定义为最终的“核心基因（hub genes）”。

3. 预后建模：

   • 基于所识别核心基因的表达水平，通过多变量 Cox 回归构建风险评分模型。
   • 将患者分为高风险组和低风险组。使用 Kaplan-Meier 生存分析和 ROC 曲线评估预后准确性。
   • 通过多变量 Cox 回归评估风险评分作为独立于临床变量（年龄、分期、性别）的独立预后因素的效力。

4. 功能与通路富集：

   • 通过 clusterProfiler 和 ReactomePA 包进行基因本体（GO）、KEGG、Reactome 和疾病本体（DO）分析，以阐明生物学过程和通路。

5. 分子对接验证：

   • 使用 AlphaFold 3 预测核心基因的蛋白质结构。
   • 进行分子对接模拟（AutoDock Vina），以评估鸢尾素与核心蛋白之间的结合亲和力。计算结合能以评估相互作用的稳定性。

主要结果

• 靶点识别： 研究确定了 93 个将鸢尾素与 HCC 联系起来的潜在靶点。
• 核心基因： 识别出八个核心基因：RACGAP1, BIRC5, NEK2, UBE2C, TYMS, EZH2, AURKB, 以及 CDK1。所有八个基因在 HCC 组织中均显著上调，并与不良预后相关（风险比 HR > 1）。尽管 AURKA 最初在交集分析中被提及，但它被排除在最终的核心基因集之外，而其余八个基因被保留了下来。
• 预后价值： 基于鸢尾素的风险特征被证实是 HCC 患者的独立预后因素，且独立于年龄、分期和性别。高风险患者的总体生存率显著较低（p < 0.001）。该模型显示出中等的预测能力，AUC 值分别为 0.76（1 年）、0.66（3 年）和 0.582（5 年）。
• 功能机制： 富集分析表明，鸢尾素相关的靶点可能参与了鸢尾素诱导 HCC 的四种机制：
  • 有丝分裂稳定性： 调节纺锤体微管向动粒的附着以及姐妹染色单体的分离。
  • 细胞周期失调： 涉及 G2/M 期转换以及 APC/C 介导的细胞周期蛋白降解。
  • 致癌活性： 增强蛋白质丝氨酸/苏氨酸激酶活性和伴侣分子结合。
  • 细胞活力： 促进肿瘤生长并产生对化学治疗药物（如抗叶酸药物）的耐药性。
• 分子对接： 所有八个核心蛋白均表现出与鸢尾素的强结合亲和力。其中七个蛋白的结合能低于 −7.2 kcal/mol，而 RACGAP1 的结合能为 −6.8 kcal/mol，表明所有靶点均具有稳健的结合稳定性和高亲和力。

意义与主张
作者声称，本研究从系统层面提供了关于鸢尾素在 HCC 中作用的视角，超越了不一致的单基因报道，实现了基于网络的理解。其主要贡献包括：

1. 识别核心靶点： 本研究发现了八个特定的核心基因（RACGAP1, BIRC5, NEK2, UBE2C, TYMS, EZH2, AURKB, CDK1），它们在鸢尾素介导的 HCC 进展中起关键作用。
2. 机制见解： 确定了四种潜在机制：有丝分裂期间的染色体不稳定/非整倍体、细胞周期失调、增强致癌蛋白活性以及增加癌细胞活力。
3. 预后工具： 验证了基于鸢尾素的 8 基因特征是患者生存情况的独立预后指标。
4. 治疗潜力： 分子对接证实的强结合亲和力表明，这些核心蛋白是未来治疗策略的有效靶点。

局限性与谦逊态度
作者明确承认这是一项计算研究。他们指出：

• 来自不同数据库的靶点数据集可能由于预测算法和标准的差异而存在异质性。
• 结合相互作用是基于计算机模拟（in silico）预测的，需要实验验证以确认实际亲和力。
• 并非所有识别出的机制都在鸢尾素诱导 HCC 的背景下得到了实验验证。
• 研究呼吁进行严格的实验验证，以充分阐明鸢尾素在 HCC 发病机制中的精确作用。

总之，本文认为，鸢尾素可能通过由八个核心基因介导的特定有丝分裂和细胞周期通路来调节 HCC 的进展，为理解肝细胞癌发生提供了新的框架，并为预防和治疗提供了潜在靶点，尚需实验验证。