- 著者: Roeder I, Horn M, Glauche I, Hochhaus A, Mueller MC, Loeffler M
- Corresponding author: Ingo Roeder (Institute for Medical Informatics, Statistics and Epidemiology, University of Leipzig, Germany)
- 雑誌: Nature Medicine
- 発行年: 2006
- Epub日: 2006-10-01
- Article種別: Original Article (Letter / mathematical modeling)
- PMID: 17013383
背景
慢性骨髄性白血病 (CML; chronic myeloid leukemia) は、Philadelphia 染色体 (Ph; Philadelphia chromosome) を生む染色体転座を起点とする造血器腫瘍である。この転座により BCR-ABL1 融合遺伝子 (Breakpoint Cluster Region-Abelson fusion gene) が形成され、その産物である構成的活性型チロシンキナーゼが悪性クローンを拡大させて正常造血を駆逐する。標準治療である TKI (tyrosine kinase inhibitor; チロシンキナーゼ阻害薬) のイマチニブ、すなわち当初化合物識別番号 STI571 (Signal Transduction Inhibitor 571) は、BCR-ABL1 陽性細胞に選択的な増殖抑制とアポトーシス誘導を介して、多くの症例で速やかな血液学的・細胞遺伝学的奏効をもたらす (Tanaka et al. JExpMed 2020 が示すように、イマチニブの作用は腫瘍免疫を含む多面的なものである)。しかし悪性クローンの完全根絶は稀で、治療中止後には急速再発が高頻度に観察される。
イマチニブ治療中の BCR-ABL1 転写産物は特徴的な二相性低下 (初期の急峻な低下 → 後期の緩徐かつ持続的な低下) を示す。先行研究 (Michor et al. Nature 2005) は、この二相性を「イマチニブが幹細胞層には実質的に作用しない」とする幹細胞非感受性モデル (悪性・正常細胞の 2 つの非相互作用集団の動態) で解釈した。しかし、この解釈が唯一の説明とは限らない。また、休眠 (quiescent) BCR-ABL1 陽性幹細胞がイマチニブの除去効果を逃れるという臨床観察 (in vitro で Ph 陽性 quiescent 幹細胞が STI571 に非感受性であるとの報告) もある。これらを単一の定量的・力学的枠組みへ統合する試みが手薄であった点が、本研究のギャップ (gap in knowledge) である。悪性クローン残存の機序と、それを踏まえた治療最適化への含意は未確立で、幹細胞層への作用様式の理解が不足していた。
目的
マウス造血で検証済みの自己組織化幹細胞モデルをヒト CML に適応し、(1) イマチニブ治療下の BCR-ABL1 二相性低下、(2) 個体間の応答不均一性、(3) 治療中止後の急速再発、という臨床観察を単一の定量的枠組みで一貫して説明できるかを検証する。さらに、悪性クローン残存を「増殖型 BCR-ABL1 陽性幹細胞への選択的作用 + 休眠型幹細胞の逃避」として捉えた場合に、イマチニブと増殖刺激の併用が悪性クローン根絶を加速し得るかを予測することを目的とした。
結果
2 独立コホートでの二相性動態の定量的再現： モデルは、性質の異なる 2 つの独立コホート — 既報の選択コホート (n=68 patients、再増加例を除外、BCR-ABL1/BCR%) と、ドイツ IRIS (International Randomized study of Interferon and STI571) コホート (n=69 patients、未選択、BCR-ABL1/ABL1%) — の双方で観察された二相性 BCR-ABL1 低下を定量的に再現した (Fig 1a)。初期の急峻な低下は、除去パラメータ r_deg による増殖中 BCR-ABL1 陽性細胞の急速な減少として説明され、後期の緩徐な低下は幹細胞数の縮小に伴う系の調節応答変化として解釈された。2 コホート間の定量差は、イマチニブ効果の大きさの小さな変動だけで説明可能で、選択コホートでは r_deg=0.033、IRIS コホートでは r_deg=0.028 という近接した最適値でデータに適合した (Fig 1a)。臨床データは合計 894 末梢血サンプル、median follow-up 60 months (range 6-66 months) に基づく。
個体間応答不均一性のパラメータ表現： 6 例の定性的に異なる BCR-ABL1 動態パターン (各コホートから 3 例ずつ) が、モデル仮定を変えることなく、阻害強度 r_inh と除去強度 r_deg の個体間差だけで再現された (Fig 1b)。急速奏効例は高い r_deg を、緩徐奏効例は低い r_deg を要した。フィッティングに用いた r_inh は 0.005-0.05 と約 10-fold、r_deg は 0.012-0.0375 と約 3.1-fold の幅をとり、これが臨床応答の広いスペクトルに対応した。各パラメータは「1 時間の時間ステップ内に個々の細胞が影響を受ける確率」として定義される。モデルはまた、CML 発症前に必要な 5-7 年の潜在期や、治療開始時に 99% 超の細胞が BCR-ABL1 陽性である状態も整合的に再現した。
休眠幹細胞の逃避と治療中止後の急速再発機序： モデルの最重要な定性的予測は、G0 (休眠) 状態の BCR-ABL1 陽性幹細胞が増殖依存性の除去効果 (r_deg) を本質的に逃れて生存し続ける、という点である。これらの休眠 Ph 陽性幹細胞は再活性化の潜在能を保持したまま治療中も残存するため、治療中止後にこれらが細胞周期へ復帰し急速再発を駆動すると予測された (Fig 1c)。シミュレーション曲線とは対照的に、実患者は中止前後で 0.01%-5% という幅広い BCR-ABL1/BCR% を示し、この個体差も残存幹細胞数の違いとして枠組み内で説明された。重要なのは、この機序が幹細胞に「一般的非感受性」を仮定せず、増殖状態に依存した「選択的感受性」だけで成立する点であり (Michor et al. Nature 2005 の幹細胞非感受性説と概念的に相違する)、休眠を逃避ニッチとする幹細胞動態の理解と整合する (Zeng et al. Nature 2026)。
長期動態と耐性クローンの寄与： 耐性が生じない条件での長期シミュレーションでは、BCR-ABL1 転写産物の持続的低下と悪性クローンの最終的根絶が予測された (Fig 2c, 黒線; r_deg=0.028, r_inh=0.050)。一方、治療開始 3 年後に 10 個の耐性幹細胞からなるクローンを導入すると、完全耐性 (r_deg=r_inh=0) では急速再発が、部分耐性 (r_deg=0.0021, r_inh=0.001) では中程度の再発が予測され、耐性強度に依存して再発速度が連続的に変化した (Fig 2c)。この多様性は、イマチニブ治療 CML で臨床的に観察される長期転帰の幅と一致する。
増殖刺激併用による根絶加速の予測： 休眠 Ph 陽性・Ph 陰性幹細胞を 0.1%/hour の速度で細胞周期へ動員する非特異的増殖刺激をイマチニブに併用したシミュレーションでは、単剤に比べて BCR-ABL1 消失が早期化し、治療中止後の腫瘍再増殖も大幅に減速すると予測された (Fig 2a; r_deg=0.033, r_inh=0.050)。この利益は完全耐性クローンが存在する場合には耐性度に応じて減弱するが (Fig 2b)、加速的な腫瘍縮小はそれ自体が新規耐性変異の発生確率を低下させると考察された。
考察/結論
本研究は、イマチニブ治療下 CML の BCR-ABL1 二相性動態・個体間不均一性・中止後再発を、幹細胞ニッチ競合の自己組織化モデルで一貫して再現した数理モデル研究である。先行の Michor et al. Nature 2005 (既報) は「イマチニブが幹細胞に無効」とする幹細胞非感受性モデルで同一観察を説明したが、本研究はこれとは概念的に異なり (in contrast)、対照的に「増殖型 BCR-ABL1 陽性幹細胞への選択的感受性 + 休眠型幹細胞の残存・逃避」という代替機序を提示した点が novel (新規) である。同一の臨床観察データが、生物学的に大きく異なる 2 つの解釈で説明され得ることを明示したこと自体が、本研究で初めて定量的に示された重要な知見である。
両モデルの相違は治療戦略に直結する。幹細胞非感受性説では既存 TKI 以外の新規薬剤が原理的に必要となるが、本モデルは「休眠幹細胞を増殖へ動員すれば既存イマチニブの除去効果に曝露できる」と予測し、増殖刺激併用による根絶加速という臨床応用 (clinical implication) の理論的根拠を提供する。これは bench-to-bedside の橋渡しとして、その後の TKI 中止 (treatment-free remission; TFR) 研究や、イマチニブと G-CSF (granulocyte colony-stimulating factor; 顆粒球コロニー刺激因子) ・インターフェロン等の増殖系刺激の併用戦略の概念的基盤となった。休眠幹細胞を再発リスクの本体と位置づける本枠組みは、がん幹細胞の休眠・治療抵抗性に関する現在の理解 (Batlle et al. NatMed 2017) とも高度に整合し、CML を超えた一般性を持つ。
残された課題 (future research) として、著者らは (1) 耐性変異の出現タイミングと確率の定量的予測、(2) 耐性クローン出現後のクローン進化の予測、(3) イマチニブの薬物動態を組み込んだモデル精緻化、を挙げる。さらに、本モデルが依拠する「ニッチ競合」が唯一の発症機序とは限らず (悪性細胞が正常細胞のアクセスできない微小環境へ侵入し競合を回避する可能性)、相対値 (BCR-ABL1/BCR または BCR-ABL1/ABL1) のみに基づく現行データではこれら仮説を判別できないという limitation も明示された。両者の区別には正常クローンの絶対サイズの経時測定が必要である。本研究は TKI 時代の CML 数理モデリングの礎石であり、BCR-ABL1 動態解析・TFR 達成条件の予測・併用治療設計に関する後続研究群の概念的枠組みを形成した。
方法
モデル設計: マウス造血幹細胞キメラ (DBA/2 ↔ C57BL/6 系統) のクローン競合動態データで実装・検証済みの自己組織化幹細胞モデル (Michor et al. Nature 2005 とは独立の系譜) をヒト CML へ適応した。モデルは 2 種の増殖微小環境を定義する — 増殖が起こらない休眠ニッチ環境 A (G0 相に相当) と、増殖が可能なニッチ外環境 O。各細胞は細胞固有の親和性 a を持ち、環境 O では a を漸減、環境 A では a を再獲得する。a が閾値 a_min を下回ると 0 に固定され、その細胞は分化細胞として一過性増幅後に固定寿命で死滅する。環境間の移行は確率過程 (stochastic process) としてモデル化され、移行確率は a と、環境 A・O の細胞数に依存する移行特性 f_a・f_o によって決まる。CML は、正常 (Ph 陰性) と悪性 (Ph 陽性) 幹細胞の f_a・f_o の量的差として表現され、これが悪性細胞のニッチ競合 (niche competition) 上の優位性を生む。イマチニブ効果は Ph 陽性細胞特異的な 2 パラメータ — 細胞周期進入抑制の阻害強度 r_inh と増殖中幹細胞の除去強度 r_deg — として組み込んだ。
検証データセットと統計: 2 つの独立コホートの median と interquartile range (IQR; 四分位範囲) を観察値とし、シミュレーション曲線との定量的適合を比較した。コホート 1 は既報の選択コホート (n=68 patients、BCR-ABL1/BCR%)、コホート 2 はドイツ IRIS (International Randomized study of Interferon and STI571) 試験コホート (n=69 patients、未選択、BCR-ABL1/ABL1%) で、新規診断 CML が定量的逆転写 PCR (RT-PCR; reverse transcriptase polymerase chain reaction) で 3 か月間隔モニタリングされた。臨床測定は計 894 末梢血サンプル、median observation time 60 months (range 6-66 months)。BCR-ABL1% は非増殖分化細胞集団の細胞数から算出した。本研究は Declaration of Helsinki に従い、各施設倫理委員会の承認と全患者の書面同意を得て実施された。
シミュレーションシナリオ: (i) 併用治療 — 休眠 Ph 陽性・Ph 陰性幹細胞を 0.1%/hour で細胞周期へ動員する非特異的増殖刺激をイマチニブに併用。(ii) 耐性 — 治療開始 3 年後に 10 個の耐性幹細胞を導入し、完全耐性 (r_deg=r_inh=0) と部分耐性 (r_deg=0.0021, r_inh=0.001) の長期動態を比較。シミュレーションプログラムのソースコードは著者から入手可能とされた。