Comparative Network Pharmacology Analysis of Ethyl Acetate, n-Butanol and Residual Fractions of Stachytarpheta jamaicensis Against Lung Cancer

技術要約：肺がんに対する Stachytarpheta jamaicensis 分画の比較ネットワーク薬理学解析

問題提起
肺がんは依然として世界の癌死亡率における主要な原因であり、転移性非小細胞肺癌（NSCLC）は、獲得耐性と遺伝的多様性によって、5年生存率が低迷している。Stachytarpheta jamaicensis (L.) Vahl は、抗酸化作用、抗関節炎作用、抗菌作用が報告されている熱帯産の薬用植物であるが、その抗腫瘍能に関する評価は乳がんおよび子宮頸がん細胞株への検証に限定されていた。本種に対して肺がんへの活性を評価した研究はなく、ネットワーク薬理学の枠組みも適用されていない。さらに、既存の植物化学研究は通常、単一の粗抽出物の特性を記述するにとどまっており、どの溶媒極性範囲に特定の疾患ターゲットに関連する化合物が濃縮されているかという理解には空白がある。本研究は、エチルアセテート、n-ブタノール、および残渣分画を比較評価し、それらの薬物らしさ（drug-likeness）、ターゲットの収束性、および肺がんに対する結合能を決定することで、これらの空白を埋めるものである。

方法論
本研究では、分子ドッキングを組み合わせた系統的なネットワーク薬理学アプローチを採用した：

• 抽出およびプロファイリング： ナイジェリアで採取された新鮮な葉を認証し、順次 n-ヘキサン、エチルアセテート、n-ブタノール、および残渣分画へと分画した。完全な LC-MS 特性評価は付随する論文に記載されているが、本研究ではそのプロファイルを出発点として利用した。
• 薬物らしさのスクリーニング： 植物化学成分を、5つのルールベースの基準（Lipinski, Ghose, Veber, Egan, および Muegge）に基づき SwissADME を用いてフィルタリングした。違反が1つ以下の化合物のみを保持した。Silychristin や Marchantin J を含む11個の化合物は、ルール違反があったものの、報告されている抗腫瘍関連性を考慮して優先的に保持された。
• ターゲット予測とネットワーク構築： 推定ターゲットは SwissTargetPrediction を通じて予測され、5種類の癌種に関連する疾患遺伝子セット（GeneCards および Open Targets）と交差させた。重複は、超幾何分布検定を用いて評価された。タンパク質相互作用（PPI）ネットワークは STRING（信頼度 ≥ 0.7）を用いて構築され、ハブ遺伝子は CytoHubba の Maximal Clique Centrality (MCC) を通じて特定された。
• パスウェイ濃縮解析： 遺伝子オントロジー（GO）および KEGG パスウェイ濃縮解析を ShinyGO v0.85.1 で行い、偽発見率（FDR）に基づき上位50個のパスウェイを保持した。
• 分子ドッキング： 肺がんにおける役割と、定義されたリガンド結合部位を持つ実験的3D構造の可用性に基づき、5つのハブ遺伝子（HIF1A, AKT1, SRC, EGFR, BCL2）をドッキング対象として優先した。ドッキングは、Standard Precision (SP) および Extra Precision (XP) スコアリングを用いる Schrödinger Suite (Glide) を使用して行われた。受容体グリッドは、共結晶リガンドの再ドッキングによって検証された。参照阻害剤（dasatinib, erlotinib, ipatasertib, venetoclax, acriflavine）はベンチマークとして使用された。

主な結果

• 薬物らしさとターゲットの重複： 同定された27個の薬物様化合物の中で、11個が優先された。薬物様化合物の割合は、エチルアセテート分画（37.5%）と残渣分画（36.4%）で最も高く、n-ブタノール分画（26.1%）で最も低かった。極めて重要なことに、3つの分画すべてにおいて、肺がん関連遺伝子との重複率はほぼ同一であった：エチルアセテート 83.06%、n-ブタノール 82.54%、残渣 82.23%。これは、抗腫瘍に関連するターゲットが、単一の極性範囲に集中しているのではなく、抽出物全体に分散していることを示している。
• ハブ遺伝子の収束： HIF1A は、3つの分画すべてにおいて MCC のトップ10にランクされた唯一のハブ遺伝子であった。AKT1, SRC, EGFR, BCL2 は、それぞれ3つの分画のうち2つの分画に現れた。
• パスウェイ濃縮： すべての分画が「Pathways in Cancer」、アポトーシス、および EGFR チロシンキナーゼ阻害剤耐性に収束した。n-ブタノール分画は、アポトーシスに対して最も強い濃縮（17.0倍）を示し、エチルアセテート分画は小細胞肺がんのパスウェイを特異的に濃縮した。
• 分子ドッキング：
  • 検証： キナーゼ標的（SRC, AKT1, EGFR）における共結晶リガンドの再ドッキングは、高い精度（RMSD 0.21–0.69 Å）を示した。しかし、タンパク質間相互作用表面（BCL2, HIF1A）における再ドッキングは信頼性が低く、参照阻害剤が一部の植物化学成分よりも高いスコアを得られなかった。また、HIF1A にはドッキング可能なポケットが存在しなかった（Dscore 0.442）。したがって、HIF1A と BCL2 は分画間の比較検討から除外された。
  • リード化合物： Silychristin（n-ブタノール分画由来）が、検証されたキナーゼ標的において臨床用阻害剤に匹敵するドッキングスコアを達成し、リーディング化合物として浮上した（SRC において −9.43 kcal/mol [dasatinib −9.51]、EGFR において −9.30 kcal/mol [erlotinib −9.25]）。これに続くのは、Secoisolariciresinol（残渣分画）であり、SRC において −8.91 kcal/mol のスコアを示した。
  • 分画比較： n-ブタノールは単一の最も良好なスコアを生み出したものの、分画間の差異は、分画内の化合物間の差異と比較して僅かであった。どの植物化学成分も、対応する臨床用参照阻害剤を上回ることはなかった。

意義と主張
本論文は、S. jamaicensis に対してネットワーク薬理学の枠組みを適用した最初の事例であると主張している。その主な意義は、化学的に異なる溶媒分画間での一貫したターゲット重複とパスウェイ濃縮によって示されるように、植物の抗腫瘍能が抽出物全体の分散的特性であることを実証した点にある。

本研究は、Silychristin および Secoisolariciresinol を、構造的に検証された部位における臨床用キナーゼ阻害剤に近いドッキングスコアによって支持される、実験的検証のための優先候補として特定した。さらに、本論文は重要な方法論的知見を提示している。すなわち、ネットワーク中心性（トポロジー）は構造的な扱いやすさを保証しないということである。HIF1A は唯一の共通ハブ遺伝子であったにもかかわらず、ドッキング可能なポケットを持たず、計算による優先順位付けは構造的な実現可能性を考慮しなければならないことを示した。

著者らは、これらの知見が S. jamaicensis の抗腫瘍能に対するメカニズム的かつターゲットレベルの根拠を確立するものであり、肺がん細胞株における今後の in vitro 検証に向けた特定の化合物を特定するものであると結論付けている。同時に、ドッキングスコアは生物学的な確認を必要とする計算上の推定値であることにも注意を促している。