- 著者: Dai-Yan Zhou, Pei Xie, Si-Yuan Wang, Wei-Lu Liu, Meng-Di Hao, Zi-Yi Wang, Xin-Xin Zeng, Xu-Hui Zhang
- Corresponding author: Xin-Xin Zeng & Xu-Hui Zhang (Guangdong Second Provincial General Hospital / Southern Medical University)
- 雑誌: Cancer Cell International
- 発行年: 2026
- Epub日: N/A
- Article種別: Original Article
- DOI: 10.1186/s12935-026-04400-5
背景
Epidermal growth factor receptor (EGFR) 変異は NSCLC (non-small cell lung cancer) の約 10〜15% に認められ、アジア人では最大 50% に達する。osimertinib (OSI、オシメルチニブ) は EGFR 変異 NSCLC の標準一次治療として確立されているが、約 50% の患者が 20 ヶ月以内に疾患進行を経験し、acquired resistance (後天性耐性) が主要な臨床的課題となっている。既知の耐性機構として MET 増幅 (30〜50%) と C797S 変異 (10〜26%) が知られているが、約半数の症例では耐性ドライバーが未解明のまま残されており (Leonetti et al. BrJCancer 2019)、新規耐性機構の探索が急務である。
extracellular matrix (ECM、細胞外マトリックス) はコラーゲンやインテグリンを介して EGFR チロシンキナーゼ阻害薬 (TKI) への耐性を dose-dependent かつ reversible に変調しうることが示されており、腫瘍微小環境が耐性形成に積極的に関与することが注目されている。しかし従来の 2D 細胞培養系では、腫瘍の三次元構造・細胞多様性・微小環境の相互作用を十分に再現できないため、2D モデルでは見逃されがちな ECM 依存性耐性メカニズムを同定するための新たなプラットフォームが不足していた (Ricordel et al. AnnOncol 2018)。
Angiopoietin-like protein 4 (ANGPTL4) はヘパリン/ヘパラン硫酸プロテオグリカン依存的機構で ECM に結合する分泌型糖タンパク質であり、肺がんを含む複数のがんで増殖・浸潤・転移を促進することが報告されている。しかしオシメルチニブ耐性における ANGPTL4 の具体的役割は未解明であり、3D オルガノイドモデルを活用した探索的研究が必要とされていた (Lee et al. Cell 2018)。
目的
EGFR 変異 NSCLC の 3D オルガノイド耐性モデルを用いたマルチトランスクリプトーム解析により ANGPTL4 をオシメルチニブ耐性関連候補として同定し、その機能的寄与と N-myc Downregulated Gene 1 (NDRG1)/ECM 軸を介した分子機構を解明するとともに、ANGPTL4 を標的とする候補低分子を計算科学的にスクリーニングすること。
結果
3D オルガノイドモデルの構築と ANGPTL4 の候補遺伝子同定：PC9 細胞 (EGFR 19-del) を Matrigel (70%) 内に 1×10⁶ cells/mL で包埋培養し、7〜10 日後に直径 100〜200 μm の三次元腺様構造を形成率 >90%・7 日生存率 >85% で確立した (Fig. 1A)。3D オルガノイドのオシメルチニブ IC50 は 2D 培養細胞の 3 倍超に達し、より臨床実態に近い耐性表現型を示した (Fig. 1B)。単純比較では捉えられない 3D 微小環境由来の耐性ドライバーを特定するため、① 2D 耐性 vs 2D 感受性、② 3D 耐性 vs 3D 感受性、③ 3D 耐性 vs 2D 耐性の 3 方向トランスクリプトーム比較を実施し、Venn 図解析で共通上方制御遺伝子群を抽出した (Fig. 1G-H)。続いて STRING タンパク質ネットワーク + OncoPredict IC50 相関 (GSE130160) + TCGA 予後解析の 3 段スクリーニングにより ANGPTL4 が最有力候補として選出された (Fig. 1I-K)。RT-qPCR による検証では ANGPTL4 mRNA 発現が osimertinib-resistant PC9 cells (PC9-OSI) で感受性細胞に比べ約 15 倍上昇しており (n=3 replicates)、Western Blot でもタンパク質レベルの強い濃縮が確認された (Fig. 2A-B)。
ANGPTL4 の機能的オシメルチニブ耐性促進と臨床相関：PC9 オルガノイドへの ANGPTL4 stable overexpression cells (OE-ANGPTL4) は IC50 を有意に増大させてオシメルチニブ感受性を低下させた (n=3 replicates、Fig. 2G-I)。逆に PC9-OSI 耐性オルガノイドへの shRNA 媒介 ANGPTL4 ノックダウン (shANGPTL4) はオシメルチニブ感受性を部分的に回復し、IC50 値が低下した (Fig. 2J-L)。この ANGPTL4 依存性耐性表現型は H1975 モデル (EGFR L858R/T790M) でも再現され、osimertinib-resistant H1975 cells (H1975-OSI) では ANGPTL4 発現が親株より増大、安定過剰発現 overexpression (OE) で IC50 増大・shRNA ノックダウンで IC50 低下が確認された (Fig. 2M-P)。臨床的相関として、オシメルチニブ治療前後の対応 NSCLC 患者組織ペア IHC 解析では、耐性相/進行相の組織で ANGPTL4 タンパク質発現が顕著に亢進していた (Fig. 2C)。6 つの独立単細胞 NSCLC データセットの統合解析では、ANGPTL4 が特定の悪性細胞サブポピュレーションに濃縮しており、drug sensitivity prediction algorithm (PERCEPTION) による予測耐性スコアも同集団で高値を示した (Fig. 2D-F)。
ANGPTL4 による ECM 経路活性化と機能的救済実験：PC9-OSI オルガノイドの shANGPTL4 後の全転写解析では、ANGPTL4 ノックダウンにより有意に発現変動した遺伝子群 (|log2FC| > 1, adj. p < 0.05) が GO で「extracellular matrix organization」「collagen fiber organization」、KEGG で「ECM-receptor interaction」「adhesion spots」に顕著に濃縮された (Fig. 3A-C)。単細胞 scRNA-seq データ (AUCell アルゴリズム) でも ANGPTL4 高発現悪性細胞の ECM 経路活性スコアが低発現群より有意に高く (Fig. 3D)、scTenifoldKnk による仮想 ANGPTL4 ノックアウトは ECM 関連経路活性の有意な抑制を予測した (Fig. 3E)。ウェスタンブロット・qPCR による分子検証では、ECM コア成分の Collagen I と laminin γ1 (LAMC1) が PC9-OSI 耐性オルガノイドで感受性細胞より有意に高値を示し (Fig. 3F-H)、ANGPTL4 KD で両タンパク質が低下・OE で増大した (Fig. 3J)。Pyrintegrin (ECM 活性化剤) を shANGPTL4 耐性オルガノイドに投与すると ANGPTL4 ノックダウンによる感作効果が逆転して IC50 が再上昇し (Fig. 3K)、一方 Cilengitide (ECM 経路阻害薬) で OE-ANGPTL4 オルガノイドの ECM を阻害するとオシメルチニブ耐性が回復された (Fig. 3L)。これらの rescue 実験は ANGPTL4-ECM 軸の機能的因果関係を直接証明した。
ANGPTL4-NDRG1 相互作用の同定と ECM 調節機能：ANGPTL4 ノックダウンにより同時に有意に低下した遺伝子として STRING タンパク質ネットワーク + 公開 ANGPTL4 免疫沈降質量分析 (IP-MS) データのスクリーニングで N-myc Downregulated Gene 1 (NDRG1) が浮上した (Fig. 4A-C)。ANGPTL4 結晶構造を受容体・NDRG1 予測構造をリガンドとした分子ドッキング計算では両分子が低結合自由エネルギーの安定複合体を形成することが示された (Fig. 4D)。Co-免疫沈降 (co-IP) 実験では ANGPTL4 抗体で沈降した複合体に NDRG1 が検出され、逆方向の NDRG1 抗体でも ANGPTL4 が引き下げられ、細胞内での直接物理的相互作用が確認された (Fig. 4F)。臨床肺がん組織切片での多重免疫蛍光染色・共焦点顕微鏡解析では ANGPTL4 と NDRG1 が細胞質内で広範な colocalization を示し (Fig. 4E)、6 データセット統合単細胞解析で両分子 mRNA に有意な正相関が認められた (Fig. 4G)。ANGPTL4・NDRG1 「ダブルポジティブ」細胞は「シングルポジティブ」や「ダブルネガティブ」細胞より ECM 経路スコアが有意に高値であった (Fig. 4H)。OE-ANGPTL4 細胞への NDRG1 ノックダウンは Collagen I・LAMC1 発現を有意に低下させ、NDRG1 が ANGPTL4 依存的 ECM 活性化の必須メディエーターであることが示された (Fig. 4J)。
バーチャルスクリーニングによる ANGPTL4 阻害候補化合物の同定：ZINC20 データベースから 208 万個の「drug-like」低分子ライブラリを構築し、AutoDock Vina を用いた高スループット分子ドッキングで結合エネルギーを主スクリーニング指標として top 50 分子を選定した (Fig. 5B-C)。続いてOpenMM (AMBER14 力場、TIP3P 水モデル) で 10 ns 分子動力学 (MD) シミュレーションを実施し、MDAnalysis により RMSD・結合界面揺動・水素結合占有率・疎水性接触面積を評価後、MM-GBSA 法で結合自由エネルギーを算出した。最終的に予測結合親和性・MD 安定性ともに良好な 5 候補化合物が選出された (Fig. 5D-F)。これらは計算科学的ノミネートであり、生化学・細胞レベルの実験的検証は今後の課題である。
考察/結論
① 先行研究との違い：従来のオシメルチニブ耐性研究は MET 増幅・C797S 変異・EGFR 非依存性バイパス経路に焦点を当てていたが、これらは全症例の約 50% しか説明できず、残り半数の耐性機構の解明が課題であった。本研究は 3D オルガノイドという臨床 in vivo 環境により近いプラットフォームを活用することで、2D 系では検出困難な ECM 依存性耐性プログラムとこれを制御する ANGPTL4/NDRG1 軸をこれまでにない形で明示した点で、既存の耐性機構研究と異なる (Ricordel et al. AnnOncol 2018)。
② 新規性：本研究で新規に示したのは、① ANGPTL4 が NDRG1 と直接物理的に結合して ECM 遺伝子活性化を促進するという新規なシグナル軸 (Co-IP・分子ドッキング・単細胞解析による三層証拠) と、② 3D オルガノイド特異的な ECM 活性化（3D 耐性 vs 2D 耐性比較で初めて検出）がオシメルチニブ耐性の機能的基盤となること、③ PC9 と H1975 という 2 異なる EGFR 変異背景で同一の耐性パターンが再現されることである。特に分泌型タンパク質 ANGPTL4 の細胞間相互作用依存性の機能が 2D では見えにくかった理由を 3D モデルで初めて解明した。
③ 臨床応用：臨床的意義として、ANGPTL4 は耐性患者組織 IHC で有意に高発現しており、治療前後の液性バイオマーカーや組織マーカーとしての応用が期待される。また ANGPTL4/NDRG1/ECM 軸は既存の EGFR TKI と ECM 経路阻害薬 (integrin アンタゴニスト等) の併用療法の合理的標的となり得る。臨床現場でのオシメルチニブ耐性後の次治療戦略として、ECM リモデリングを標的とした介入が新たな選択肢となる可能性がある。計算科学でノミネートした 5 候補化合物はさらなる実験的検証を経て ANGPTL4 阻害薬開発の出発点となり得る (Leonetti et al. BrJCancer 2019)。
④ 残された課題：今後の研究として、PC9・H1975 以外の EGFR 変異モデル・患者由来オルガノイド・in vivo PDX モデルでの検証が必要である。また ANGPTL4 が NDRG1 に作用する具体的な分子レベル (タンパク質安定性・ユビキチン化・転写因子活性化等) の解明と、ECM リモデリング下流の effector 同定が今後の課題として残された。仮想スクリーニングで選出した 5 化合物の ANGPTL4 直接結合・シグナル阻害・ECM 効果・in vivo 薬効・安全性の実験検証も不可欠である。さらに ANGPTL4 と PERCEPTION 予測耐性スコアの相関が有意ながらも中等度であることは、ANGPTL4/NDRG1/ECM 軸が広域耐性ネットワークの一要素に過ぎないことを示唆し、EMT・低酸素・薬物透過性変化等との協調解析が求められる。
方法
PC9 (EGFR exon 19 欠失変異) および H1975 (EGFR L858R/T790M) のヒト NSCLC 細胞株を使用した。オシメルチニブ耐性株 (PC9-OSI: 1 μM 維持、H1975-OSI: 0.5 μM 維持) を 6 ヶ月超の段階的濃度漸増法 (10 nM→1 μM) で樹立した。3D オルガノイドは Matrigel (70%, 1×10⁶ cells/mL, 50 μL/well) を 37°C で 15〜20 分固化後、organoid 専用培地 (Advanced DMEM/F12 + B27 + 20 ng/mL EGF + 10 ng/mL FGF) で培養し 7〜10 日で構造形成を確認した。安定発現株は慢性ウイルスベクターシステム (shANGPTL4、OE-ANGPTL4、shNC、VectOSI) + ピューロマイシン選択で作製した。トランスクリプトーム解析は BGI・Illumina NovaSeq で実施し、GRCh38 アラインメント後 featureCounts でカウント、DESeq2 で差次発現解析 (|log2FC| > 1、adj. p < 0.05 を有意とした)。GO 機能注釈 + KEGG 経路濃縮解析 (clusterProfiler)、単細胞 scRNA-seq は Seurat で前処理 + AUCell で ECM 経路活性スコアリングを実施した。細胞生存率は CCK-8 法 (3000〜5000 cells/well、96 ウェルプレート、72 h、450 nm 吸光度測定) で評価し、IC50 は GraphPad Prism 9 による非線形回帰で算出した。RT-qPCR は QuantStudio 6 Flex、SYBR Green、2^(-ΔΔCt) 法。タンパク質解析は RIPA 溶解液 + BCA 定量 + SDS-PAGE 後 PVDF 転写 (抗体: anti-ANGPTL4 1:1000/Abcam、anti-NDRG1 1:1500/CST、anti-Collagen I 1:1000/Proteintech、anti-LAMC1 1:1000/Santa Cruz、anti-GAPDH 1:5000/Proteintech)。バーチャルスクリーニングは ZINC20 (208 万分子) を AutoDock Vina でドッキング後 top 50 を選別し、OpenMM (AMBER14 + TIP3P、10 ns MD) + MDAnalysis (RMSD・MM-GBSA) で安定性評価した。統計は Student’s t-test・one-way ANOVA、Kaplan-Meier + log-rank 法を用い、p < 0.05 を有意水準とした (全実験 n ≥ 3 回独立実施、データは mean ± SD で表示)。