- 著者: Stein BD, Ferrarone JR (共同筆頭著者), Gardner EE, Chang JW, Wu D, Hollstein PE, Liang RJ, Yuan M, Chen Q, Coukos JS, Sindelar M, Ngo B, Gross SS, Shaw RJ, Zhang C, Asara JM, Moellering RE, Varmus H, Cantley LC
- Corresponding author: Lewis C. Cantley (lewis_cantley@dfci.harvard.edu), Dana-Farber Cancer Institute; Benjamin D. Stein (bds2005@med.cornell.edu); Harold Varmus (varmus@med.cornell.edu), Weill Cornell Medicine
- 雑誌: Cancer Discovery
- 発行年: 2023
- Epub日: 2023-05-01
- Article種別: Original Article
- PMID: 36715544
背景
KRAS は肺腺癌 (LUAD: lung adenocarcinoma) で最頻変異がん遺伝子であり (CancerGenomeAtlasResearchNetwork et al. Nature 2014)、TP53 または STK11/LKB1 (liver kinase B1) の機能喪失変異と頻繁に共起する。LKB1 は AMPKR (AMPK-related kinase) ファミリーを活性化してエネルギー恒常性・代謝ストレス応答を制御する腫瘍抑制キナーゼで、その喪失は代謝リプログラミング・免疫回避・免疫療法抵抗性と関連する。先行研究では以下のギャップが欠けていた・足りなかった: (i) KRAS/TP53/LKB1 三重変異がヒトコホートで期待値より有意に少ない現象の分子機序が不明であった (Zehir et al. Cancer Discov 2017)、(ii) マウス GEMM (genetically engineered mouse model) では同三遺伝子変異が高悪性度の肺腺癌を生むというヒトと逆の表現型矛盾の機序的説明が欠落していた、(iii) LKB1 直下の AMPKR (AMP-activated protein kinase related kinase) サブファミリーのうちどれが LUAD 関連代謝制御を担うか特定されていなかった、(iv) ヒト LKB1 変異肺癌に対する GEMM ベースの前臨床薬剤開発の妥当性が 批判的に検証されていなかったため、マウスで有効な治療がヒトで再現できない理由を説明する分子的枠組みが完全に不足していた。
TCGA 511 例の解析では KRAS・TP53・LKB1 三遺伝子の同時変異はわずか8例 (odds ratio = 0.35、p = 0.01) と期待値より有意に少ない。このヒトでの統計的稀少性 vs マウス GEMM での増悪促進という根本的矛盾を分子レベルで説明することが、ヒト治療開発における前臨床モデル妥当性の課題として急務であった。
目的
ヒト・マウス LUAD における KRAS/TP53/LKB1 三重変異の表現型差異の分子基盤を解明する。具体的には (a) ヒト LKB1 喪失が KRAS/TP53 変異LUADで selective lethality をもたらす機序、(b) その下流の AMPKR サブファミリー (SIK1/2/3 等) と基質の同定、(c) マウスで反対表現型を生む分子的差異 (進化的アミノ酸置換) の特定、(d) GEMM のヒト LKB1 変異 LUAD モデルとしての妥当性評価を目指す。
結果
ヒト KPL 肺腺癌は選択的増殖障害だがマウス KPL は逆に促進: TCGA + MSKCC で KRAS/TP53/LKB1 三重変異 odds ratio = 0.35、p = 0.01 (Fig. 1A) と稀少性を確認 (vs 期待値約2.86倍の差)。ヒト H358 KP 株は皮下ゼノグラフトで腫瘍を形成したが等質遺伝子型 KPL 株 (sgLkb1-2.1, sgLkb1-3.2) は 腫瘍形成不能 (n=10 mice/group, p < 0.001、Fig. 1C-E)。3D スフェロイドでも KPL が著明な増殖障害 (sphere area vs KP で-70%、Fig. 1F)。マウス KP・KPL は両方腫瘍形成可能で KPL ではむしろ増殖促進 (vs KP で+1.5-fold、Fig. 1G)。WT LKB1 再発現は KPL ゼノグラフト増殖を完全回復するが KI LKB1 では救済不可能で、LKB1 のキナーゼ活性必須性を証明。AMG-510 (KRAS G12C 阻害薬) はヒト KP スフェロイドを抑制 (IC50 ~10 nM) するが KPL では部分的救済 (vs vehicle で +40% sphere area)、LKB1 喪失が RAS 依存性を低下させる。
TPI1 Ser21 リン酸化が LKB1 依存的かつヒト特異的: 定量的リン酸化プロテオミクス (n=4 replicates per genotype) でヒト KP vs KPL でTPI1 Ser21 リン酸化が最も有意に低下 (fold change = -8.5、adj. p < 0.001、Fig. 2A)、TPI1 全タンパク量は不変。P-Ser21 特異的抗体で LKB1 依存性確認 (Fig. 2B、WT LKB1 再発現で完全回復)。マウス Tpi1 21 残基は Cys21 に置換されており KP・KPL 両方でリン酸化検出不能 (Fig. 2D)。哺乳類15種以上の進化的配列解析でほぼ全種で Ser21 保存だが、マウス・ラットのみ Cys21 に置換を確認 (Fig. 2E、種間 conservation score 14/15 vs 1/15)。
SIK1/2/3 が LKB1 下流の TPI1 Ser21 キナーゼ: AMPKR 各サブファミリー個別 KO で SIK TKO のみが LKB1 存在下での TPI1 Ser21 リン酸化を著明低下 (vs LKB1 alone で-92%、Fig. 3B)。SIK TKO は LKB1 喪失と同等の代謝欠陥 (OCR 最大容量 -45%、ECAR -38%、p < 0.01) とスフェロイド形成障害を再現 (Fig. 3D-F)。リン酸化プロテオミクスで SIK3 Ser551 (SIK 活性調節部位) が KPL 株で最も有意に低下したリン酸化サイトの一つ。CRTC3 の SIK 基質サイト (Ser62・Ser329・Ser370) も KPL で低下し SIK 軸活性化が中枢を担うことを多角的に確認。SIK1 mRNA は KPL でフィードバック増加。
TPI1 リン酸化が解糖/脂質産生フラックスの分岐を制御: ヒドロキシルアミントラッピングで KPL 細胞での DHAP (dihydroxyacetone phosphate) 蓄積 +3.2-fold (vs KP、p < 0.001) と G3P (glycerol-3-phosphate) +2.5-fold 増加 (Fig. 4A)。¹³C₂ グルコース MFA で KPL の M+2 G3P フラックスが KP と比較し全時点で有意に増加 (Spearman correlation r = 0.85 between LKB1 loss vs G3P flux, p < 0.05、Fig. 4C-E、cohort_n=3-5/timepoint)、TPI1 リン酸化が DHAP→GAP 変換を維持し解糖系炭素供給を担保していることを示す。LKB1 喪失時には DHAP が脂質合成経路へ迂回しグリセロール脂質合成 +2-fold 増加・解糖系効率低下を生む。ヒト KPL の正常グルコース ECAR は KP のグルコース制限条件と同等まで低下し著明な代謝障害を示す。BN-PAGE で LKB1 喪失により TPI1 二量体分率が増加 (vs control 1.8-fold、Fig. 4G)、thermal proteome profiling でリン酸化 TPI1 の Tm が ΔTm = 5.8°C 低下 (Fig. 4H、リン酸化が二量体不安定化と活性増進を生む構造モデルと整合)。
マウス Cys21 は酸化修飾で LKB1 非依存的代替調節: マウス KPL に過酸化水素処理 (200 μM, 30 min) で低グルコース下に TPI1 単量体シフトが LKB1 状態に依存せず独立して誘導 (BN-PAGE、Fig. 5C-D)。これはマウス Cys21 がスルフィン酸/スルホン酸形成による酸化修飾で TPI1 活性を調節し、LKB1/SIK に依存しない代替経路を持つことを示す。ヒト Ser21 はシステイン酸化調節経路を持たないため KRAS/TP53 変異下での代謝的柔軟性を LKB1/SIK 依存的リン酸化のみで担保せざるを得ず、LKB1 喪失が致死的となる選択的代謝的脆弱性が生まれる。
考察/結論
本研究は「LKB1 喪失がヒト KRAS/TP53 変異LUADでは致死的だがマウス GEMM では増殖を加速する」という長年の表現型矛盾に対し、TPI1 Ser21 vs Cys21 の進化的アミノ酸差異という一点で統合的説明を与えた点が新規 (novel) な貢献である。LKB1→SIK1/2/3→TPI1 Ser21 リン酸化→解糖/脂質産生分岐制御 という単一経路がヒト腫瘍細胞の代謝恒常性を担保し、これが喪失すると DHAP 蓄積と脂質産生不全が KRAS/TP53 変異の高エネルギー要求状態で致命的となる機序を分子・代謝・構造の3次元で解明した。
先行研究との対照・整合性: 先行の AMPK 中心の LKB1 下流解析と異なり、本研究は SIK ファミリー (SIK1/2/3) が肺腺癌特異的に LUAD 代謝制御の中枢を担うことを直接示した。先行の小細胞肺癌 (SCLC) における種間モデル比較研究 (Gay et al. CancerCell 2019) も human/mouse cross-comparison 重要性を示しているが、本研究は単一アミノ酸残基レベルでメカニズム的差異を pinpoint した点が異なる。
bench-to-bedside translation・臨床応用: KRAS/TP53/LKB1 変異ヒトLUADは現在も免疫療法 (PD-1/PD-L1 阻害薬) への反応が最不良サブタイプの一つで、LKB1-SIK-TPI1 軸の代謝脆弱性はこのサブタイプ特有の新規治療標的として高い臨床的意義を持つ。SIK 活性化薬または TPI1 機能補正は概念的治療介入方向となり、解糖中間体蓄積の細胞毒性を逆用する戦略 (解糖促進薬との組み合わせ) も検討に値する。KRAS G12C 阻害薬 (AMG-510 等) によるスフェロイド救済効果は、LKB1 喪失が RAS 依存性を変化させることを示し、KRAS 阻害薬感受性に対する LKB1 ステータスの影響評価を臨床試験設計に組み込む示唆を与える。さらにヒト Ser21 とマウス Cys21 の進化的分岐は GEMM がヒト LKB1 変異LUADの治療開発モデルとして根本的限界を持つことを示し、マウスで有効性が示された介入のヒト試験移行失敗リスクを説明する。
残された課題・limitations: (i) SIK アイソフォーム (SIK1/2/3) 個別寄与の組織特異性は未解明 (TKO のみで効果評価)、(ii) ヒト LUAD 臨床検体での TPI1 P-Ser21 発現と予後・治療応答との直接相関は未確立で、前向き観察研究での検証が残された課題、(iii) 過酸化水素処理は非生理的酸化ストレスで、in vivo マウス腫瘍でのCys21 酸化修飾レベルの定量データが不足、(iv) SIK 活性化薬・TPI1 機能補正薬の前臨床候補化合物がまだ開発初期段階、(v) NSCLC 以外のがん種 (KRAS変異膵癌・大腸癌) での同経路寄与は未検証、(vi) 免疫療法応答との関連性 (LKB1 変異 = ICI 抵抗性既知) における TPI1-DHAP 代謝産物による免疫微小環境への影響評価が残る。
方法
臨床コホート: TCGA Pan-Cancer Atlas (n=511 LUAD) と MSKCC 臨床コホート (n=1,357) で KRAS・TP53・LKB1 三遺伝子変異共存頻度を Fisher exact test で評価。
等質遺伝子型細胞株: ヒト KP 株 (H358、H2009、KRAS+TP53 変異) と GEMM 由来マウス KP 株を親株として、CRISPR-Cas9 で LKB1 をノックアウト (sgLkb1-2.1、sgLkb1-3.2 の 2 sgRNA) してヒト KPL・マウス KPL を樹立。WT LKB1 / KI (kinase-inactive K78I) LKB1 再発現 control 株も並列作製。
機能アッセイ: 3D スフェロイド形成 (n=8 wells/condition、3 biological replicates)、免疫不全マウス (NSG/Nude) 皮下ゼノグラフト・同所移植 (n=5-10 mice/group、合計 KP vs KPL で 30 匹超)、AMG-510 (KRAS G12C 阻害薬) 処理応答 (IC50 を sphere area で評価)。
定量的リン酸化プロテオミクス: 低グルコース条件 (0.5 mmol/L) で 4-6時間培養後、LC-MS/MS (Orbitrap Fusion Lumos) で TMT 10-plex 標識 + IMAC enrichment、n=3-4 biological replicates per genotype/species。TPI1 P-Ser21 特異的抗体を新規開発しウェスタンブロットで独立検証。
メタボロミクス・フラックス解析: ヒドロキシルアミン in situ 化学トラッピングメタボロミクス (DHAP・GAP・G3P 定量、Asara lab)。¹³C₂ 標識グルコース (1,2-¹³C₂) トレーサー MFA (metabolic flux analysis) で M+2 G3P フラックスを 0.5・1・2・5 分の経時採取で測定 (cohort_n=3-5 per timepoint)。Seahorse XF96 で OCR (oxygen consumption rate)・ECAR (extracellular acidification rate) を測定。
SIK ファミリー同定: AMPKR 各サブファミリー (MARK・NUAK・BRSK・AMPK・SNRK) の個別ノックアウト + LKB1 再発現で TPI1 Ser21 リン酸化を比較。SIK1/2/3 triple knockout (SIK TKO) で決定的評価。
TPI1 構造解析: BN-PAGE (blue-native PAGE) で TPI1 二量体/単量体比率、thermal proteome profiling で P-Ser21 TPI1 の Tm を測定 (n=3 replicates)。
進化解析: 哺乳類複数種 (>15 種) の TPI1 アミノ酸 21 残基を比較し Ser/Cys 種間分布を決定。マウス KPL に過酸化水素 (200 μM) 処理で Cys21 酸化修飾依存的 TPI1 活性調節を機能検証。