[论文解读] Mathematical Model of DSB formation by Ionizing Radiation

J. Barilla, M. Lokajíček+1 more|ArXiv.org|Jan 31, 2008

Radiation Effects and Dosimetry参考文献 2被引用 5

一句话总结

本论文提出了一种关于电离辐射诱导DNA双链断裂（DSB）形成的数学模型，重点研究水中的自由基簇引发的间接损伤。通过使用优化速率常数和氧气溶解度效应的速率常数模拟，该模型成功再现了不同氧气浓度下的实验DSB产率，表明OH•自由基在DSB形成中占主导地位（α_OH : α_H ≈ 93.53 : 5.95），而直接作用可忽略不计（p₁ ≈ 0.0003）。

ABSTRACT

The understanding of inactivation radiobiological mechanism in individual cells is important when from one side the application of ionizing radiation to tumour therapy and from the other side the protection against radiation are to be effectively improved. Important part of this mechanism is double-strand break (DSB) formation; these DSB may be formed directly by impacting ionizing particles or indirectly by different secondary radicals. The latter kind of formation is much more frequent when cells contain normal water content. Mathematical model of the corresponding chemical stage will be presented with the aim to demonstrate how the individual radicals, but also other chemical agents present in the water during irradiation may influence DSB formation.

研究动机与目标

• 理解放射生物学损伤的化学阶段，特别是通过水中次级自由基形成的DSB。
• 量化放射敏感剂或放射防护剂在低LET辐射照射下对DSB产率的影响。
• 建立动力学模型，准确再现不同氧气浓度下实验测得的DSB形成数据。
• 确定不同自由基（OH•、H•、e⁻aq、O₂）对DNA中DSB形成的相对贡献。
• 实现对不同化学环境下DSB产率的预测，支持放射防护与放疗优化。

提出的方法

• 建立常微分方程组，模拟电离辐射形成的球形自由基簇内的自由基反应。
• 通过假设自由基密度随时间因空间扩散而降低，引入扩散效应。
• 利用文献中已知的相对反应速率比约束参数，绝对值通过MINUIT优化。
• 采用非线性函数[O₂] = A₀{1 − exp[−(cσ₀)τ₀]}模拟水中的氧气溶解度，以反映溶解度饱和效应。
• 校准初始自由基数量（N_H• = 18.11，N_OH• = 31.9，N_e⁻ = 35.0）和反应速率比（α_H : α_OH : α_e : α_O = 5.95 : 93.53 : 0.01 : 0.58）。
• 将模型输出与实验DSB产率曲线进行验证，尤其关注低氧浓度下的非线性响应。

实验结果

研究问题

• RQ1 在低LET辐射下，不同水衍生自由基（OH•、H•、e⁻aq、O₂）对DNA中DSB形成的相对贡献是什么？
• RQ2 氧气浓度如何影响DSB产率？是否可考虑溶解度饱和效应进行定量建模？
• RQ3 导致DSB形成的自由基簇的平均尺寸和组成是什么？
• RQ4 与直接电离相比，间接机制在DSB形成中占主导地位的程度如何？
• RQ5 自由基簇内扩散与反应动力学如何影响最终的DSB产率？

主要发现

• 通过电离粒子直接撞击DNA导致DSB形成的概率可忽略不计（p₁ ≈ 0.0003）。
• OH•自由基对DSB形成贡献最大，其相对速率常数为93.53，远高于H•（5.95）。
• 初始自由基簇在t₀时刻包含约18.11个H•、31.9个OH•和35.0个e⁻aq自由基。
• 该模型准确再现了实验观察到的高氧浓度下DSB产率曲线平坦的现象，归因于溶解度饱和。
• 氧气溶解度函数[O₂] = 0.0004{1 − exp[−(c×43.68)^1.7]}与实验数据吻合良好，其中A₀ = 0.0004 nm⁻³与纯水的溶解度一致。
• 反应速率常数k₁至k₁₇已优化，其中k₁₄ = 325.54和k₁₅ = 0.098在涉及O₂和自由基的二级反应中起关键作用。