＜報	  告＞	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 Report	 	 日本放射線影響学会第 59 回大会ワークショップ	 「計算シミュレーションによる放射線生物研究」開催報告	 日本原子力研究開発機構 1、電力中央研究所 2	 佐藤	  達彦 1*、浜田	  信行 2	 Report on the Workshop “Radiation Research Based on Computer Simulation”Held at the 59th Annual Meeting of the Japanese Radiation Research Society1Japan Atomic Energy Agency, 2Central Research Institute of Electric Power IndustryTatsuhiko Sato1*, Nobuyuki Hamada2	 	  放射線生物研究に関するシミュレーションに焦点を当てたワークショップを日本放射線影響学会の年会で開催した。ワークショップには約 50 名の参加があり、国内でシミュレーションによる放射線生物研究を実施している研究者 6 名がそれぞれの研究テーマを紹介し、今後の連携可能性について議論した。本稿では、ワークショップの概要をまとめるとともに、各発表内容を簡単に紹介する。	 
キーワード：シミュレーション、数理モデル、分野間連携	 	 A workshop entitled “Radiation Research Based on Computer Simulation” was held at the 59th annualmeeting of the Japanese Radiation Research Society. Approximately 50 participants attended the workshop, andsix speakers gave a talk about their computational radiation biology studies. In the workshop, possible futurecollaborations were intensively discussed. This report briefly summarizes the outline of the workshop.	 
Key words: Computer simulation, mathematical model, interdisciplinary interaction	 	 *	 〒319-1195	  茨城県那珂郡東海村白方 2-4	 2-4, Shirakata, Tokai, Naka, Ibaraki 319-1195 JapanTEL: +81-29-282-5803, FAX: +81-29-282-6768, E-mail: sato.tatsuhiko@jaea.go.jp-	 397	 -	 ワークショップの概要	 近年の計算科学技術の進展により、放射線生物研究においても計算シミュレーションの重要性が世界的に高まっている。日本でも、これまで DNA 損傷・修復、細胞致死、治療計画、発がんなど様々な研究テーマに対応したシミュレーション技術が開発されてきた。しかし、それらの開発は個別に進められることが多く、異なるテーマ間でシミュレーション研究に関して統合的に議論する機会はなかった。そこで、2016 年 10 月 26 日～28 日にかけて広島市で開催された日本放射線影響学会第 59 回大会において、著者らが座長となり「計算シミュレーションによる放射線生物研究」と題したワークショップを開催した。ワークショップには約 50 名の参加があり、国内でシミュレーションによる放射線生物研究を実施している研究者 6 名が講演した。以下、そのプログラムを示す（敬称略）。	 
1.

「DNA 損傷・飛跡構造解析に関する研究」渡辺立子（量子科学技術研究開発機構: 以下、量研機構）

2.

「照射・非照射細胞混在環境を模擬した確率的モデルによる細胞応答に関する研究」佐々木恒平（北海道科学大学）	 

3.

「低線量放射線生体影響研究における数理モデルの有用性」真鍋勇一郎（大阪大学）	 

4.

「細胞生存率モデルを用いた粒子線治療に関する研究」稲庭拓（量研機構）	 

5.

「白内障の自然発症モデルの構築」坂下哲哉（量研機構）	 

6.

「時空間的異質性を考慮した発がん数理モデル解析」川口勇生（量研機構）	 全体として、ミクロからマクロな影響に対するシミュレーション研究となるよう演題順番を決定した。また、ワークショップでは、各テーマで必要となる入力情報やそこから得られる出力情報を整理することにより、テーマ間の連携可能性について重点的に議論した。その議論をふまえて作成した各研究テーマの相関を図１に示す。	 
	 
各発表について	 渡辺らは、放射線の飛跡構造を解析するモンテカルロシミュレーションコードを開発し、ミクロな領域の線量分布やそれに伴うラジカル生成・化学反応・DNA 損傷の計算を可能とした(1,2)。そして、放射線の線量や線質が異なる条件下で生成する鎖切断、塩基損傷などの DNA 損傷の量,	 そして,	 二本鎖切断や複雑なクラスター損傷を形成する空間分布を推定した。これらの情報は、以下に示す全てのシミュレーション研究の入力情報となり、それをどのような形で他のシミュレーションに渡すかが、今後の連携研究を進展させる重要な鍵となると思われる。	 
佐々木らは、非標的効果による DNA 損傷がどのようなメカニズムで引き起こされるか解明するため、照射・非照射細胞混在環境でのバイスタンダーシグナル伝達シミュレーションを実施した(3)。
その際、ギャップジャンクション経由と培地経由で伝達する 2 種類のシグナルがあると仮定するこ-	 398	 -	 とにより、照射・非照射細胞毎の DNA 損傷生成率の時間変化を再現することに成功した。この結果と渡辺らのシミュレーションによって計算される標的効果による DNA 損傷を加えることにより、様々なシミュレーションで非標的効果まで考慮した計算が可能になると期待できる。	 
	  真鍋らは、低線量率被ばくにおける直線しきい値なし（LNT）仮説の真偽を解明するため、生体機能（細胞増殖・修復機能・細胞死など）を考慮して外部刺激と生体との相互作用（突然変異発生頻度）を表現する数理モデルを開発した(4,5)。その結果、様々な生物種の突然変異発生頻度を汎用的な数理モデルで再現できることが確認でき、低線量被ばくでは総線量ではなく線量率が重要であることが示唆された。今後、外部刺激として DNA 損傷数などを利用したり、得られた突然変異発生頻度を発がん数理モデルに入力したりすることにより、さらに定量的な議論が行えるようになると考えられる。	 
	  稲庭らは、細胞核中の直径 1	 µm 程度のターゲットに付与されるエネルギーから細胞致死に繋がる DNA 損傷数を推定し、その時間変化を表現する微分方程式を解くことにより細胞致死率を評価する Microdosimetric	 Kinetic モデル（通称 MK モデル）の開発・改良を行っている(6)。MK モデルを炭素線治療計画ソフトウェアに組み込むことにより、照射時間や照射の時間構造が治療効果にもたらす影響の評価が可能となる。本研究は、放射線生物に関する基礎研究が実社会に直接的に役立てられた貴重な成功例と言えるであろう。	 
	  坂下らは、眼の水晶体の細胞群動態と白内障の発症の関係を数理モデルの観点から明らかにすることを目的とした研究を実施している。これまでに、細胞群が徐々に分化していく過程において外部刺激によりダメージが蓄積して白内障が発生する、という仮説に基づく確率論的数理モデルを構築し、疫学研究から得られた白内障自然発生データ(7)を再現することに成功した。今後、外部刺激によるダメージと DNA 損傷などを関連づけることにより、放射線被ばくによりヒトに誘発される白内障の線量や線量率応答が評価可能となると期待できる。	 
	  川口らは、線量率や線質の違いによる線量分布の時空間的な異質性に着目し、2 次元格子シミュレーションに基づく発がん数理モデルを開発している(8,9)。本モデルでは、線量分布の異質性とフリーラジカル生成やそれに伴う DNA 損傷の関係や、DNA 損傷と細胞の致死率や突然変異率の関係が入力情報として必要となり、現在は、簡易的な仮説に基づいてそれらの関係性を推定している。
今後、それらの関係性をこれまでに解説した様々なモデルで得られた知見に基づいて更新することにより、本数理モデルは、放射線防護分野でも利用可能な発がんリスク評価モデルに発展可能であると思われる。	 
	 まとめ	 
	  シミュレーション研究に焦点を当てたワークショップを日本放射線影響学会において初めて開催した。これまでの放射線影響研究は、DNA 損傷・修復、放射線治療効果、発がんなど生物学的なエ-	 399	 -	 ンドポイントを主題として縦割り的に実施される傾向にあり、分野間の連携が十分とは言えない。
欧州では、MELODI（Multidisciplinary	 European	 LOw	 Dose	 Initiative）と呼ばれる低線量影響に関する学際的イニシアチブがあり、シミュレーションや数理モデルを介した分野間の連携研究が活発化する傾向にある。一方、日本では、まだそのような組織的な活動には至っていない。本ワークショップや本稿がシミュレーションによる放射線生物研究やその連携研究の活性化に繋がれば幸いである。	 

謝辞	 
	  本ワークショップの開催の機会をご恵与いただいた日本放射線影響学会第 59 回大会委員の先生方、演者の方々、聴衆として参加・ご質問していただいた方々,	 そして,	 本報告の機会をいただいた編集委員長の松本義久先生に感謝いたします。	 
	 
参考文献	 
1.

Watanabe	  R,	  Rahmanian	  S,	  Nikjoo	  H.	  Spectrum	  of	  radiation-induced	  clustered	  non-DSB	 damage	  –	  A	  Monte	  Carlo	  track	  structure	  modelling	  and	  calculation.	  Radiat	  Res	  2015;	 183:	 525-40.	 

2.

渡辺立子.	 重粒子線トラックと DNA 損傷（解説）.	 放射線化学	 2015;89:13-9.	 

3.

Sasaki	 K,	 Wakui	 K,	 Tsutsumi	 K,	 Itoh	 A,	 Date	 H.	
 A	 simulation	 study	 of	 the	 radiationinduced	  bystander	  effect:	  modeling	  with	  stochastically	  defined	  signal	  reemission.	 
Comput	 Math	 Methods	 Med	 2012;	 389095.	 

4.

Manabe	  Y,	  Bando	  M.	  Comparison	  of	  data	  on	  mutation	  frequencies	  of	  mice	  caused	  by	 radiation	 with	 low	 dose	 model.	 J	 Phys	 Soc	 Jpn	 2013;	 82(9):	 094004.	 

5.

Manabe	 Y,	 Wada	 T,	 Tsunoyama	 Y,	 Nakajima	 H,	 Nakamura	 I,	 Bando	 M.	 Whack-A-Mole	 Model:	 Towards	  unified	  description	  of	  biological	  effect	  caused	  by	  radiation-exposure.	  J	 Phys	 Soc	 Jpn	 2015;	 84(4):	 044002.	 

6.

Inaniwa	 T,	 Suzuki	 M,	 Furukawa	 T,	 Kase	 Y,	 Kanematsu	 N,	 Shirai	 T,	 Hawkins	 RB.	 Effects	 
of	 dose-delivery	 time	 structure	 on	 relative	 biological	 effectiveness	 for	 therapeutic	 carbon-ion	 beams	 evaluated	 with	 microdosimetric	 kinetic	 model.	 Radiat	 Res	 2013;	 180:	 44-59.	 

7.

Klein	  BE,	  Klein	  R,	  Lee	  KE.	  Incidence	  of	  age-related	  cataract:	  the	  Beaver	  Dam	  Eye	 Study.	 Arch	 Ophthalmol	 1998;	 116(2):	 219-25.	 

8.

Kawaguchi	 I,	 Sasaki,	 A.	 The	 wave	 speed	 of	 intergradation	 zone	 in	 two-species	 lattice	 Muellerian	 mimicry	 model,	 J	 Theor	 Biol	 2006;	 243(4):	 594-603.	 

-	 400	 -	 9.

Doi	  M,	  Kawaguchi	  I.	  Ecological	  impacts	  of	  umbrella	  effects	  of	  radiation	  on	  the	 individual	 members.	
 J	 Environ	 Radioact	 2007;	 96(1/3):	 32-8.	 	 

	 
図１	  本ワークショップで紹介した各研究の相関。現在は、このようなネットワークは形成されていないが、連携可能性のあるテーマや入出力情報は矢印で相関づけた。	 

-	 401	 -