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

技術要約：肝細胞癌におけるイリシンの分子メカニズムの解明

問題提起
肝細胞癌（HCC）は、依然として世界的に癌関連死亡の主要な原因であり、多くの場合、根治的な治療オプションが限られている進行した段階で診断される。運動が疾患リスクを低減し、生活の質を向上させることは知られているが、運動によって誘導されるマイオカイン、特にイリシンが、どのような分子メカニズムを通じてHCCの病態形成を調節しているのかについては、その詳細なメカニズムが十分に解明されていない。既存の文献では、イリシンのHCCに対する役割に関して相反する証拠が示されている。イリシンがPI3K/Akt経路を介して増殖や移動を促進すると示唆する研究もあれば、腫瘍組織において影響がない、あるいはむしろ発現が低下していると報告する研究もある。さらに、正常対照と比較したHCC組織におけるイリシンの発現レベルについても合意が得られていない。本研究は、システム生物学的なアプローチを用いて、HCCにおけるイリシンの潜在的な分子標的とシグナル伝経路を包括的に特定することにより、これらの不一致を解消することを目的としている。

方法論
本研究では、ネットワーク毒性学と分子ドッキングを統合した多段階のバイオインフォマティクス・パイプラインを用いた：

1. 標的の収集:

   • イリシンの標的: SwissTargetPrediction、ChEMBL、およびSTITCHデータベースから、イリシンによって調節される可能性のあるヒト標的遺伝子を回収した。
   • HCCの標的: GeneCardsおよびOMOMから、肝臓癌に関連する遺伝子を精査した。
   • 発現変動: TCGA-LIHCコホート（腫瘍サンプル374、正常サンプル50）のRNA-seqデータを用い、DESeq2を使用して、発現変動遺伝子（DEGs）を特定した（カットオフ値：adj. P < 0.05, |log₂FC| > 1）。
   • 交差解析: イリシンの標的、HCC関連遺伝子、およびTCGA由来のDEGsの交差部分から、93個の潜在的な標的を導出した。

2. ネットワーク構築とハブ遺伝子の特定:

   • タンパク質相互作用（PPI）ネットワークをSTRINGデータベース（信頼度スコア > 0.4）を用いて構築し、Cytoscapeで可視化した。
   • モジュール解析: MCODEプラグインを用いて、最も重要なサブネットワークを特定した。
   • ハブの選択: cytoHubbaプラグイン内の5つのアルゴリズム（Degree, Betweenness, Closeness, MNC, Radiality）を用いて遺伝子をランク付けした。これらのアルゴリズムによる上位遺伝子の交差部分を、最終的な「ハブ遺伝子」と定義した。

3. 予後モデリング:

   • 特定されたハブ遺伝子の発現レベルに基づき、多変量Cox回帰を用いてリスクスコアモデルを構築した。
   • 患者を高リスク群と低リスク群に層別化した。カプラン・マイヤー生存分析およびROC曲線を用いて、予後予測精度を評価した。
   • 多変量Cox回帰を用いて、年齢、ステージ、性別といった臨床変数とともに、リスクスコアが独立した予後因子であるかを評価した。

4. 機能および経路の濃縮解析:

   • clusterProfilerおよびReactomePAパッケージを使用し、生物学的プロセスおよび経路を解明するために、Gene Ontology (GO)、KEGG、Reactome、およびDisease Ontology (DO) 解析を実施した。

5. 分子ドッキングによる検証:

   • ハブ遺伝子のタンパク質構造をAlphaFold 3を用いて予測した。
   • イリシンとハブタンパク質との間の結合親和性を評価するため、分子ドッキングシミュレーション（AutoDock Vina）を実施した。相互作用の安定性を評価するために、結合エネルギーを算出した。

主な結果

• 標的の特定: 本研究では、イリシンとHCCを結びつける93個の潜在的な標的を特定した。
• ハブ遺伝子: 8つのコアハブ遺伝子（RACGAP1, BIRC5, NEK2, UBE2C, TYMS, EZH2, AURKB, CDK1）を特定した。これら8つの遺伝子はすべて、HCC組織において有意にアップレギュレートされており、不良な予後（ハザード比 > 1）と関連していた。なお、AURKAは当初の交差解析で強調されていたが、最終的なコア遺伝子セットからは除外され、残りの8つの遺伝子が保持された。
• 予後価値: イリシンに基づくリスクシグネチャーは、HCC患者における独立した予後因子として機能し、これは年齢、ステージ、性別とは異なるものである。高リスク患者は、有意に低い全生存率を示した（p < 0.001）。モデルの予測能は中程度であり、AUC値は0.76（1年）、0.66（3年）、0.582（5年）であった。
• 機能的メカニズム: 濃縮解析により、イリシンに関連する標的が、イリシン誘発性HCCにおける4つのメカニズムに関与している可能性が明らかになった：
  • 有糸分裂の安定性: 紡錘体微小管のキネトコアへの結合および姉妹染色分体の分離の調節。
  • 細胞周期の調節不全: G2/M期への移行およびAPC/Cによる細胞周期タンパク質の分解への関与。
  • 腫瘍形成活性: タンパク質セリン/スレオニンキナーゼ活性およびシャペロン結合の増強。
  • 細胞生存能: 腫瘍増殖および化学療法剤（例：抗葉酸剤）に対する耐性の促進。
• 分子ドッキング: すべての8つのハブタンパク質は、イリシンに対して強い結合親和性を示した。7つのタンパク質は−7.2 kcal/molを下回る結合エネルギーを示し、RACGAP1は−6.8 kcal/molの結合エネルギーを示した。これは、すべての標的に対して強固な結合安定性と高い親和性を持つことを示唆している。

意義と主張
著者らは、本研究が単一遺伝子の報告による矛盾を超え、ネットワークベースの理解を通じて、HCCにおけるイシンの役割に関するシステムレベルの視点を提供するものであると主張している。主な貢献は以下の通りである：

1. コア標的の特定: 本研究は、イリシンによるHCCの進行において極めて重要な8つの特定のハブ遺伝子（RACGAP1, BIRC5, NEK2, UBE2C, TYMS, EZH2, AURKB, CDK1）を明らかにしている。
2. メカニズムの洞察: 4つの潜在的なメカニズム（有糸分裂中の染色体不安定性/異数性、細胞周期の調節不全、腫瘍形成タンパク質活性の増強、およびがん細胞生存能の向上）を特定した。
3. 予後ツール: イリシンに関連する8遺伝子のシグネチャーは、患者の生存に対する独立した予後指標として検証されている。
4. 治療の可能性: 分子ドッキングによって確認された強い結合親和性は、これらのハブタンパク質が将来の治療戦略における実行可能な標的であることを示唆している。

限界と謙虚な姿勢
著者らは、本研究が計算科学的研究であることを明示的に認めている。彼らは以下の点に留意している：

• 様々なデータベースからの標的データセットは、予測アルゴリズムや基準の違いによる不均一性を伴う可能性がある。
• 結合相互作用はin silico（コンピュータ上）で予測されたものであり、実際の親和性を確認するには実験的な検証が必要である。
• 特定されたメカニズムのすべてが、イリシン誘発性HCCの文脈において実験的に検証されているわけではない。
• 本研究は、イリシンによるHCCの病態形成における正確な役割を完全に解明するために、厳格な実験的検証を求めている。

結論として、本論文は、イリシンが8つのハブ遺伝子を介した特定の有糸分裂および細胞周期経路を通じて、HCCの進行を調節している可能性があると提唱しており、実験的な確認を待たずとも、肝癌発生の理解および予防・治療のための標的特定に向けた新しい枠組みを提供している。