アルファ線治療の精度向上へ：アスタチン-211甲状腺線量評価の新手法
解説対象論文: Monte Carlo microdosimetry of 211at in multi-species thyroid follicle models: AlfaMC versus MCNPX and Geant4
EJNMMI Physics｜被引用数: 0（2026年8月18日時点）
がん治療の新たな希望として注目されるアルファ線治療において、アスタチン-211（211At）は強力な放射性核種です。しかし、その治療効果を最大化し、副作用を最小限に抑えるためには、生体内での線量評価、特にリスク臓器である甲状腺への正確な被ばく線量推定が不可欠です。本記事では、この重要な課題に取り組む最新の研究を紹介します。新たなモンテカルロシミュレーションコード「AlfaMC」が、既存のコードであるMCNPXやGeant4と比較され、その計算精度と効率性が検証されました。この研究は、より安全で効果的なアルファ線治療の実現に向けた大きな一歩となるでしょう。
はじめに：アルファ線治療と甲状腺被ばくの課題
どうも、Beyond the Pixelです。がんの転移性疾患に対する治療法として、放射性核種を用いた治療（放射性核種治療）が注目されています。特に、アルファ粒子を放出するアスタチン-211（211At）で標識された薬剤は、小さな腫瘍に対して高い吸収線量を集中させることができるため、非常に有望です。211Atから放出されるアルファ粒子は、ごく短い範囲（水中で最大70µm）で全エネルギーを放出し、隣接する細胞を含む局所的な照射を引き起こします。これにより、高い細胞死をもたらすことが期待されています。しかし、ヨウ素と同じハロゲン族である211Atは、生体内、特に遊離型の場合に甲状腺に蓄積する性質があります。甲状腺は、治療における主要なリスク臓器となるため、211Atによる甲状腺の放射線吸収線量を正確に評価することは極めて重要です。本研究は、甲状腺モデルにおける211Atのマイクロドシメトリー量を、モンテカルロ（MC）コードであるAlfaMCとGeant4を用いて評価し、先行研究でMCNPXを用いたデータと比較することで、AlfaMCの有効性を検証することを目的としています。特に、AlfaMCはアルファ粒子のマイクロドシメトリー計算に特化して開発されており、計算時間の短縮が期待されています。
研究の背景：アスタチン-211と甲状腺被ばくの課題
放射性核種治療は、転移性腫瘍疾患の対象者に対し、腫瘍組織への吸収線量を正常組織への線量に比べて十分に高めることを目指す治療法です。211Atは半減期7.2時間のアルファ粒子放出核種で、2つの崩壊経路を経て安定な鉛-207になります。放出されるアルファ粒子の主なエネルギーは5.867 MeV（41.8%）と7.450 MeV（58.2%）です。アルファ粒子は、その短い飛程と高い線エネルギー付与（LET）により、DNA二重鎖切断（DSB）やより複雑なDNA損傷を高頻度で引き起こし、高い細胞死率につながると考えられています。さらに、短い飛程のため、隣接する正常組織への吸収線量は非常に低いという利点もあります。これらの特性から、211Atはがん治療、特に小さな腫瘍に対する治療において注目されています。しかし、遊離型の211Atは、ヨウ素と同様に甲状腺に大きく蓄積することがマウスやラットでの研究で示されており、ヒトにおいても甲状腺疾患の対象者で高い取り込みが報告されています。そのため、甲状腺は211At標識薬剤を用いた放射性核種治療における副作用のリスクが最も高い正常臓器の一つと考えられています。アルファ粒子の短い飛程と高いLETの組み合わせは、隣接する細胞におけるエネルギー付与の大きなばらつきを生じさせます。そのため、放射線による生物学的効果を推定する際に、腫瘍やリスク臓器の平均吸収線量を用いることは不正確である可能性があります。マイクロドシメトリーは、付与されたエネルギーの確率的性質を考慮するため、アルファ粒子照射にはより適切です。
マイクロドシメトリー計算手法の比較
マイクロドシメトリー量を計算する方法はいくつかありますが、本研究では特にモンテカルロ（MC）シミュレーション技術に焦点を当てています。MC技術は、乱数を用いて放射線輸送の確率的プロセスをシミュレートし、各相互作用プロセスの結果を確率密度関数で記述します。アルファ粒子輸送をシミュレートできるMCコードはいくつか存在します。ロスアラモス国立研究所で開発された汎用コードMCNPXや、CERNで開発されたオブジェクト指向MCツールキットGeant4は、広く利用されてきました。Geant4-DNAプロジェクトは、Geant4ツールキットを分子レベルの計算に拡張し、放射線生物学や放射線治療に役立つものとなっています。これらの汎用MCコードは、複雑な形状で30種類以上の粒子輸送をシミュレートできますが、コードのインストール、物理プロセスの選択、入力ファイルのデザインには専門知識とプログラミングスキルが必要です。例えば、MCNPXは入力カードでタスクを定義し、Geant4はC++オブジェクト指向プログラミング言語に関する深い知識を必要とします。また、MCNPの配布は制限されています。さらに、汎用コードの速度は、AlfaMCのような専門コードが対象とする特定の問題に対しては最適ではない場合があります。AlfaMCは、リスボン大学のLuis PeraltaとAlina Louroによって開発された、アルファ粒子に特化したMCコードです。元々は細胞照射用に開発され、マイクロドシメトリー量の計算にカスタマイズされています。連続減速近似（CSDA）を使用し、エネルギーおよび角度の揺らぎ（ストラグリング）を考慮しています。また、組み合わせ幾何学に基づいた強力な幾何学ソフトウェアパッケージが含まれています。他のより複雑なMCコードと比較して、計算時間が短いことが大きな利点です。コード内で最も頻繁な操作である阻止能計算は、事前に計算されたルックアップテーブルを適用することで最適化されています。典型的なシミュレーションでは、10万個のアルファ粒子を追跡するのに約2分しかかかりません。AlfaMCは、運動エネルギー1 MeVから12 MeVのアルファ粒子の薄膜へのエネルギー付与の推定において、SRIMコードと良好な一致を示しています。これまで、ヒト、ラット、マウスの甲状腺モデルは、211Atや他の放射性ヨウ素の線量評価にMCNPXを用いて開発・適用されてきました。実験的な検証が不可能なため、独立した理論的データの相互比較は重要です。本研究の目的は、AlfaMCとGeant4を用いて、異なる甲状腺濾胞モデルにおける均一および不均一な211At分布のマイクロドシメトリー量を決定し、先行研究のMCNPXデータと比較することでした。
研究方法：多種甲状腺濾胞モデルとアスタチン-211の分布
本研究では、先行研究で発表されたヒト、ラット、マウスの単一および複数甲状腺濾胞モデルを、AlfaMCおよびGeant4コードにそれぞれ対応する幾何学ツールを用いて定義しました。
単一濾胞モデル
ヒト、ラット、マウスの単一甲状腺濾胞モデルは、球状の甲状腺濾胞、球状の濾胞内腔、濾胞内腔を囲む濾胞細胞層、および中心から等距離に左右対称に分布する6つの球状濾胞細胞核から構成されています（
）。各構成要素は液体水の密度（1.0 g/cm3）を持ちます。濾胞内腔の直径、濾胞細胞層の厚さ、濾胞細胞核の直径は、それぞれヒトで150、10、8 μm、ラットで70、8、6 μm、マウスで50、6、4 μmでした。
複数濾胞モデル
各生物種に合わせた先行研究の複数甲状腺濾胞モデルは、甲状腺濾胞が隣接する濾胞に囲まれている場合の単純化された表現です（
左）。このモデルでは、隣接する細胞の濾胞細胞が、中心濾胞の濾胞細胞の周囲と、それらが合体した内腔の周囲の両方に球状のシェルを形成しています（（図2）右）。周囲の濾胞内腔と周囲の細胞の厚さは、ヒトで150 μmと10 μm、ラットで70 μmと8 μm、マウスで50 μmと6 μmでした。
線源の分布
本研究では、放出されるアルファ粒子の主要な2つのエネルギーである5.867 MeV（41.8%）と7.450 MeV（58.2%）のみを考慮しました。アルファ粒子の等方性放出を仮定し、以下の4種類の211At分布でシミュレーションを行いました。
1. 濾胞内腔内に均一に分布する211At
2. 濾胞細胞内に均一に分布する211At
3. 濾胞細胞核内に均一に分布する211At
4. 濾胞内腔の中心に位置する同心球の表面に均一に分布する211At（球の半径は0から150 μmまで5 μm刻みで変化）。これは不均一な211At分布を表します。
これらの4つの線源局在は、遊離型アスタチンに対して考えられるもので、特に1、2、4が最も可能性が高いとされています。アスタチンが薬剤に結合している場合は、濾胞細胞の基底膜を表す同心球の表面への分布（4番）が最も適切です。すべてのシミュレーションにおいて、ターゲットは6つの濾胞細胞核でした。
モンテカルロシミュレーション
平均比エネルギー、単一ヒット平均比エネルギー、単一ヒット比エネルギー分布など、いくつかのマイクロドシメトリーパラメータを決定しました。MCシミュレーションはAlfaMCとGeant4を使用して行い、MCNPXを用いた先行発表データと比較しました。
AlfaMCでは、10^7個のアルファ粒子を追跡し、カットオフエネルギーは1 keVに設定されました。各計算から6つのターゲットに対して6つの独立したマイクロドシメトリー推定値が得られ、統計的精度を高めるためにシミュレーションを2回繰り返すことで12の独立した値が決定されました。
Geant4では、単一甲状腺濾胞モデルの場合に、線源が濾胞内腔、細胞、6つの核に均一に分布し、ターゲットが6つの核となるように実装されました。アルファ粒子輸送は、低エネルギー輸送に推奨される標準的な物理コンストラクタを使用してシミュレートされました。
結果：AlfaMCの計算精度
本研究で計算されたAlfaMC、Geant4の結果と、MCNPXを用いた先行研究の結果を比較しました。
単一甲状腺濾胞モデル
単一甲状腺モデルにおいて、211Atが細胞核、濾胞細胞、または内腔に均一に分布していると仮定した場合の、平均比エネルギーと単一ヒット平均比エネルギーの値が示されています（
）。AlfaMCとMCNPXの結果の差は、最大で4.6%でした。また、AlfaMCとGeant4の結果の差は、最大で1.7%でした。標準偏差は、すべてのMC結果において0.4%以内でした。AlfaMCの結果は、MCNPXよりもGeant4の結果とより良い一致を示しました。特に、211Atが細胞核内に均一に分布している場合、線源とターゲットが同一であるため、平均比エネルギーと単一ヒット平均比エネルギーの値は等しくなります。
3つの線源分布と各生物種モデルに対する単一ヒット比エネルギー分布は、AlfaMC、MCNPX、およびGeant4で決定された結果と一致しました（
）。確率分布は曲線下の面積が1になるように正規化されています。
同心球の表面に211Atが均一に分布している場合のAlfaMCの結果（
a）と、MCNPXによる対応する結果（（図4）b）も比較されました。AlfaMCとMCNPXのデータ間の比率は、ほとんどのデータポイントで1に近く、ヒトモデルの一部のデータポイント（半径10μmと15μmの表面線源）でそれぞれ2.2と1.2の比率が得られました。標準偏差は各測定で計算され、誤差棒として示されていますが、その値が小さいため、図には見えない箇所もあります。
複数甲状腺濾胞モデル
複数甲状腺濾胞モデルにおいて、211Atが濾胞内腔に位置している場合と、同心球の表面に均一に分布している場合の単一ヒット平均比エネルギーと平均比エネルギーが示されています（
）。AlfaMCとMCNPXのデータ間で強い一致が観察され、特に平均比エネルギーの比較では一貫して高い類似性が見られました。
考察：シミュレーションコード間の差異と臨床的意義
本研究の全体的な目的は、MCコードAlfaMCによって決定されたマイクロドシメトリーパラメータをMCNPXおよびGeant4の結果と比較することでした。平均比エネルギー、単一ヒット平均比エネルギー、および単一ヒット比エネルギー分布が調査されました。平均比エネルギーはゼロヒットイベントを含む質量あたりの平均比エネルギーを表し、平均吸収線量のS値計算に利用できます。単一ヒット平均比エネルギーは、単一イベントによってのみ付与された平均比エネルギーに関する情報を含んでいます。平均吸収線量は、ヒットしたアルファ粒子の数、あるいは確率分布と乱数発生器を用いて、各アルファ粒子の寄与をサンプリングすることで計算できます。
この研究は、ヒト、ラット、マウスの甲状腺濾胞モデルにおける211At曝露に適用されました。使用したモデルはすべて多細胞性で、単一または複数の甲状腺濾胞をシミュレートしていました。多濾胞モデルは甲状腺全体の状態を反映しており、生体内または臨床状況を最もよく表す興味深いモデルでした。また、濾胞内の異なる線源分布からの線量寄与を提示することは、異なる生物種間で得られた違いをよりよく理解するためにも重要であると考えられました。これは、対象者が211At標識薬剤による治療を受ける際の副作用を避けるために、甲状腺への吸収線量を最適化する際に役立つ可能性があります。最近のいくつかの論文では、類似の目的でin vitroの細胞に関する実験的および/またはマイクロドシメトリー研究が報告されています。さらに、他の放射性核種に対する異なる放射性核種分布に関する線量評価研究も、単一甲状腺濾胞モデルで行われています。したがって、単一濾胞モデルの結果も提示することで、研究間や放射性核種間の比較を容易にすることにしました。アスタチンは多くの価数状態をとり得るハロゲンであり、その生理学的挙動はヒトでは非常に限られており、動物でも限定的です。アスタチンはヨウ素と類似の挙動を示すと考えられていることもありますが、動物では多くの違いも示されています。本研究で考慮された4つの211At分布は、濾胞細胞および濾胞内腔における遊離型211Atの潜在的な局在をシミュレートしており、これらは可能性の高い分布です。ただし、細胞核内の局在は可能性が低いとされています。さらに、薬剤に結合した211Atは、濾胞細胞の基底膜を表す同心球の表面への分布によってシミュレートできます。
AlfaMC、MCNPX、Geant4の結果間の違いは、アルファ粒子の輸送とエネルギー付与の記述方法の違いに起因すると考えられます。AlfaMCが基づいているASTARデータベースは、高エネルギー（>2 MeV）アルファ粒子の衝突阻止能における不確実性が1～4%であるとされています。ブラッグピークテールのエネルギー領域である100 keV以下のエネルギーでは、不確実性は10%を超えます。したがって、ブラッグピークテールにおけるわずかな違いがエネルギー付与の変動につながる可能性があります。MCNPXおよびGeant4は、表形式または解析的な断面積を使用する物理モデルを採用しています。相互作用モデリングにおける近似と不確実性が、マイクロドシメトリー量の評価に違いを生じさせる可能性があります。MCNPXは、多くの衝突プロセスを単一のステップにグループ化する「凝縮履歴アプローチ」に基づいており、プロトンからの阻止能データをアルファ粒子にスケーリングして使用していました。Geant4で計算された阻止能は、AlfaMCで適用されたNIST/ASTAR阻止能データベースと同じものと比較されました。二次電子の輸送は、1マイクロメートルに設定された生成カットオフとステップ制限によって制御されました。カットオフ制限以下の飛程を持つ二次電子は、そのエネルギーを局所的に付与します。Geant4 DNAプロセスでは、すべての相互作用が明示的にシミュレートされるため、生成カットは使用されませんでした。核ターゲット領域のみでGeant4 DNA Physicsを有効にしても、Geant4の結果は実質的に変化しませんでした。MCNPXとの差は減少しましたが、統計的不確実性の範囲内でした。Geant4 DNA物理モデル全体を実装すると計算時間が長くなりすぎるため、本研究では実施されませんでした。1000個の粒子をシミュレートするのに約1時間かかり、統計的不確実性を1%未満にするには少なくとも5×10^7個の粒子が必要でした。ナノメートルスケールでの飛跡評価とエネルギー付与には、アルファ粒子のイベントごとの相互作用をシミュレートするDNA Physicsの使用が一般的に認識されていますが、マイクロメートルスケールではまだ未解決の疑問があります。許容可能な計算時間を達成するために、異なる物理モデルを持つ領域を組み合わせる効果は、さらなる調査が必要です。特に、DNA Physics領域をターゲットを超えて拡張する効果は、将来の研究で取り組むことができます。シミュレーションジオメトリも結果に影響を与える可能性のある要因でした。ファントムのボクセルサイズがMCNPXとGeant4のマイクロドシメトリー結果に影響を与えることが示されていますが、本研究で得られたヒト、ラット、マウスの結果間の系統的な違いは、サイズ効果が無視できることを示唆しています。異なるモデルとMCコードを用いてアルファ粒子輸送をシミュレートして得られたマイクロドシメトリーパラメータの相互比較は、マイクロドシメトリーモデルの検証に対する効果的な方法として用いられました。公表された結果と5%以内の合意は良好と見なされます。分析モデルとGeant4-DNAの結果の間で、直径1〜10 μmの球状ジオメトリで最大9%の偏差は許容範囲内でした。オーガー電子放出核種（99mTc）のサブセルラーレベルでの線量評価研究では、MCNP6コードが甲状腺濾胞モデルを用いたMCNPXの結果と比較して3%以内の一致を示しました。マイクロドシメトリーが提供する、特定のエネルギー分布に関する情報は、治療応用やアルファ粒子に対する細胞応答の理解において重要です。計算集約的な問題に関連するマイクロドシメトリーの欠点は、ターゲットアルファ治療のための高度な線量評価システムの可能性を開くことで減少しつつあります。巨視的レベルでの粒子輸送は、対象者固有の線量送達と生物学的効果の評価のために、マイクロドシメトリー運動モデルと結合されています。
甲状腺の被ばくに関するマイクロドシメトリー推定は、211At標識薬剤が非甲状腺がんの治療に使用される場合に重要です。甲状腺は、骨髄および潜在的に腎臓とともに、この治療法におけるリスク臓器です。甲状腺に焦点を当てる理由は、「遊離型211At」が放出され、NIS（ヨウ化ナトリウム共輸送体）を介して甲状腺に大きく蓄積するためです。臨床試験は多くの場合、マウスやラットを用いた前臨床研究に先行します。本研究で示されているように、甲状腺濾胞内腔のサイズ、細胞の厚さ、および細胞核の直径の違いは、アルファ粒子の飛程が限られているため、特に濾胞内腔の211Atからの寄与に関して、齧歯類からヒトへの甲状腺への吸収線量の外挿が実質的に異なることを明確に示しています。ヒトにおけるより現実的な線量評価のためには、ヒトの大きな濾胞サイズを含む現実的な濾胞モデルが重要です。なぜなら、211Atが甲状腺全体に均一に分布すると仮定した場合よりも、濾胞細胞とその核への吸収線量は低くなるからです。それでも、対象者の放射線安全のためには、甲状腺における211Atの取り込みと保持を減らすべきであり、これは安定ヨウ素や過塩素酸塩のような「甲状腺ブロッカー」を単回または複数回投与することで大いに達成できます。異なる幾何学的および生物動態学的/機能的甲状腺モデルを評価して、さらなる甲状腺保護戦略を検討する際には、迅速で一般的なユーザーフレンドリーなMCシミュレーションが望ましいです。本研究で示されたデータは、そのような研究にはMCNPXやGeant4よりもAlfaMCが好ましいことを示しています。
結論：AlfaMCの可能性
本研究により、AlfaMCの線量評価結果は、MCNPXおよびGeant4の結果と5%以内で一致することが確認されました。これは、AlfaMCコードが甲状腺の細胞およびサブセルラーレベルでのマイクロドシメトリー計算において正確であることを示しています。MCNPXおよびGeant4と比較して、AlfaMCの計算時間がはるかに短いこと、およびアルファ粒子輸送に特化してカスタマイズされていることから、AlfaMCは、生物種間のマイクロドシメトリーパラメータを比較する場合や、211At標識薬剤で治療する対象者の甲状腺への吸収線量を減らすための異なる戦略を評価する場合に、より好ましく使用できるでしょう。
用語集
- 211At (アスタチン-211): アルファ粒子を放出する放射性核種で、がんの標的治療に利用される。
- AlfaMC: アルファ粒子輸送とマイクロドシメトリー計算に特化し、高速計算が可能なモンテカルロシミュレーションコード。
- MCNPX (Monte Carlo N-Particle eXtended): 様々な粒子輸送をシミュレートできる汎用的なモンテカルロコード。
- Geant4 (GEometry ANd Tracking): 粒子輸送をシミュレートするためのオブジェクト指向モンテカルロツールキット。
- マイクロドシメトリー: 細胞や細胞内小器官などミクロなレベルでの放射線エネルギー付与の確率的性質を扱う線量測定法。
- アルファ粒子: ヘリウム原子核と同じ構造を持つ高エネルギー粒子で、短い飛程で高い線エネルギー付与能を持つ。
- 線エネルギー付与 (LET): 放射線が物質中を通過する際に単位飛程あたりに失うエネルギー量。LETが高いほど生物学的効果も高い傾向にある。
- 甲状腺濾胞: 甲状腺の基本構造単位で、ホルモンを貯蔵する内腔とその周囲の濾胞細胞からなる。
- 平均比エネルギー: ターゲットに付与された全エネルギーをそのターゲットの質量で割った平均値。
- 単一ヒット平均比エネルギー: ターゲットにエネルギー付与があったイベント（ヒット）のみを考慮した平均比エネルギー。
- 確率密度関数: 連続的な確率変数の相対的な発生確率を示す関数。
- 連続減速近似 (CSDA): 粒子が連続的にエネルギーを失うと仮定して、その飛程を計算する近似法。
- ストラグリング: 放射線が物質中を通過する際に、エネルギー損失や飛程にばらつきが生じる現象。
Leave a Reply