关于《抗氧化剂》期刊2025年发表论文《FLASH放疗中的氧耗竭与细胞抗氧化剂作用：来自蒙特卡洛辐射化学建模的机制见解》的学术研究报告
一、 研究作者、机构与发表信息
本研究的主要作者为Israth Rabeya, Jintana Meesungnoen和Jean-Paul Jay-Gerin*（通讯作者）。三位作者均来自加拿大舍布鲁克大学（Université de Sherbrooke）医学与健康科学学院医学影像与放射科学系。该研究以题为《Oxygen depletion and the role of cellular antioxidants in flash radiotherapy: mechanistic insights from monte carlo radiation-chemical modeling》的论文形式，发表于学术期刊《Antioxidants》2025年第14卷第406期。论文于2025年2月25日收稿，3月17日修订，3月27日接受，并于3月28日正式在线发表。
二、 学术背景与研究目的
本研究隶属于放射生物学与放射化学的交叉领域，聚焦于一种新兴的肿瘤放疗技术——FLASH放疗（FLASH radiotherapy, FLASH-RT）。FLASH放疗的核心特征在于使用远超常规放疗的超高平均剂量率（通常超过40-150 Gy/s），在毫秒级时间内递送高剂量辐射。大量临床前研究表明，FLASH放疗能在有效杀伤肿瘤的同时，显著减轻对周围正常组织的损伤，此现象被称为“FLASH效应”。然而，导致FLASH效应的精确分子机制尚不明确，这严重制约了其临床转化与应用。
在众多解释FLASH效应的假说中，放射分解氧耗竭（radiolytic oxygen depletion, ROD）假说是最受关注的机制之一。该假说认为，超高剂量率的辐射脉冲会瞬时消耗细胞内的氧气，造成急性缺氧状态，从而暂时性地增强了细胞（尤其是正常细胞）对辐射的抗性。然而，ROD假说存在争议，焦点在于脉冲照射期间消耗的氧气量是否足以在初始氧合良好的组织中引发显著的缺氧，从而完全解释观察到的保护效应。
除了氧耗竭，细胞内的抗氧化防御系统也可能在FLASH效应中扮演关键角色。细胞内的抗氧化剂，如谷胱甘肽、抗坏血酸、一氧化氮和α-生育酚，能够清除辐射诱导产生的活性氧和自由基，从而减轻氧化损伤。但它们在FLASH放疗特定条件下的具体作用机制，尤其是在与ROD假说的关联中，尚未被充分阐明。
因此，本研究旨在通过计算建模的方法，从一个辐射化学的视角，系统评估ROD假说，并深入探究细胞内源性抗氧化剂在FLASH放疗中的作用。研究团队开发了一个计算模型，模拟在类似于受限细胞空间的含水环境中，受到瞬时高能质子脉冲照射时的辐射化学过程。具体目标包括：1）量化超高剂量率下细胞内氧气的瞬态消耗过程及其主要反应路径；2）评估过氧自由基（ROO•）的命运，特别是其通过自由基-自由基重组反应被中和的可能性；3）阐明谷胱甘肽、抗坏血酸、一氧化氮和α-生育酚等关键抗氧化剂在清除有害自由基、中断过氧化链式反应中的作用；4）确定FLASH效应（按ROD假说预测）得以完全显现所需的临界剂量率阈值；5）综合评估ROD假说单独解释FLASH效应的充分性。
三、 详细研究流程与方法
本研究是一项纯粹的计算模拟研究，未涉及湿实验。其核心流程基于一个名为 IONLYS-IRT 的蒙特卡洛多径迹化学模拟代码。该代码由研究团队长期开发并优化，专门用于模拟水溶液（在此扩展为模拟细胞环境）在辐射作用下的物理化学过程。研究流程可概括为以下几个关键步骤：
1. 模型构建与参数设定： * 模拟环境： 将细胞简化为一个均质的、充氧的含水介质。尽管承认真实细胞的异质性（如细胞质、细胞核、膜结构），但本研究采用均相动力学模型作为初步探索工具，以识别关键反应物和反应。 * 化学成分： 模型介质主要包含水，并添加了以下溶质以模拟细胞内环境： * 碳基生物分子（RH）： 代表DNA、蛋白质、脂质等，浓度为1 M（为考察敏感性，也测试了0.5 M）。 * 辐射诱导的生物自由基（R•）： 初始浓度设定为2.5 µM，模拟由RH直接电离产生的碳中心自由基。 * 溶解氧（O₂）： 浓度为30 µM，模拟正常哺乳动物组织中的典型氧浓度。 * 关键抗氧化剂： 谷胱甘肽（GSH， 6.5 mM）、抗坏血酸（Ascorbate, AH⁻， 1 mM）、一氧化氮（•NO， 1 µM）和α-生育酚（TOH， 0.2 mM）。这些浓度均基于文献报道的生理或应激状态下的水平设定。 * 辐射条件： 模拟使用300 MeV能量的质子进行照射，其线性能量传递（linear energy transfer, LET）约为0.3 keV/µm，模拟类似钴-60 γ射线或快电子的低LET辐射。采用“瞬时脉冲”（Dirac）模型，即假设脉冲持续时间为零，所有化学物种瞬间形成。通过改变每个脉冲中入射质子数（N = 20, 30, 40, 50）来模拟不同的剂量率。其中，N=20对应约10⁶–10⁷ Gy/s的瞬时剂量率，是典型的FLASH照射条件。
2. 化学反应方案集成： 在IONLYS-IRT代码原有的纯水辐解反应方案基础上，本研究集成了一个包含48个额外反应的扩展方案，以描述细胞环境中的关键化学反应。这些反应涵盖了： * 水辐解初级产物（如水合电子e⁻ₐq、羟基自由基•OH）与生物分子RH、抗氧化剂的反应。 * 生物自由基R•的形成（主要通过•OH从RH上抽氢）及其与O₂反应生成过氧自由基ROO•。 * 过氧自由基ROO•与各种抗氧化剂（GSH, AH⁻, •NO, TOH）的清除反应。 * 抗氧化剂之间的相互作用及次级反应（如谷胱甘肽硫自由基GS•的形成与反应）。 * 氧消耗的关键反应路径（如R• + O₂ → ROO• 和 GSSG•⁻ + O₂ → GSSG + O₂•⁻）。 所有反应的速率常数均从已发表的辐射化学和生物化学文献中选取。
3. 蒙特卡洛模拟执行： * 物理与物理化学阶段（IONLYS程序）： 模拟辐射能量沉积的初始过程（~1皮秒内），生成水辐解初级产物（e⁻ₐq, H•, •OH, H₂, H₂O₂等）在径迹内的非均匀空间分布。 * 化学阶段（IRT程序）： 采用“独立反应时间”（Independent Reaction Times, IRT）随机模拟方法，计算从1皮秒到1秒时间尺度内，所有反应物种的扩散及它们之间、以及与溶液中均匀分布的溶质（RH, O₂, 抗氧化剂等）之间的化学反应。该方法高效且准确，能够处理从径迹非均匀动力学到溶液均相动力学的完整过程。
4. 数据分析与输出： 模拟的核心输出是随时间变化的化学产额（G值），定义为每吸收100 eV能量所形成或消耗的分子数。研究重点关注以下关键物种的G值时间演化曲线： * 反应性物种：e⁻ₐq, •OH, R•, ROO• * 底物消耗：RH, O₂ * 抗氧化剂相关中间体 通过分析这些曲线，以及计算特定反应对总G值的贡献（ΔG），可以定量揭示不同反应路径的重要性。此外，利用公式 [–O₂] ≈ 5.3 × 10⁻⁴ × N × LET × G(–O₂) 将氧消耗的G值转换为实际消耗的氧浓度（µM），从而直接评估缺氧程度。
5. 敏感性分析： 为了评估模型假设的稳健性，研究进行了两项关键的敏感性分析： * 生物分子浓度影响： 将RH的浓度从1 M降低至0.5 M，重复模拟以观察对氧消耗结果的影响。 * 扩散系数影响： 考虑到细胞内分子扩散可能比在纯水中慢，将所有反应物种的扩散系数降低为纯水中的1/100，重复模拟以评估扩散环境变化对结果的影响。
四、 主要研究结果
1. 过氧自由基的命运与抗氧化剂的核心作用： 模拟结果显示，在FLASH照射后的微秒至毫秒时间尺度内，形成了显著的过氧自由基ROO•（见图1，ROO•曲线呈钟形）。然而，与Labarbe等人先前假设不同，这些ROO•并非主要通过自由基-自由基重组反应（如ROO• + ROO• 或 ROO• + R•）被中和。 图2b清晰地表明，ROO•的消耗主要由其与抗氧化剂的反应驱动。其中，抗坏血酸（AH⁻）和一氧化氮（•NO）是清除ROO•最有效的抗氧化剂，而谷胱甘肽（GSH）和α-生育酚（TOH）的作用相对较小。这一发现挑战了先前认为FLASH放射保护主要源于高浓度自由基间重组的观点，并突出了快速反应的细胞抗氧化剂，尤其是AH⁻和•NO，在有效清除辐射诱导的过氧自由基、从而防止有害的过氧化链式反应和减轻氧化损伤方面的关键作用。
2. 氧气消耗的机制与定量分析： 模拟揭示了细胞内氧气消耗的动态过程（图3a）。氧气消耗主要始于微秒级时间，通过两个主要反应路径进行（图3b）： * 路径一（已知）： R• + O₂ → ROO•。这与先前研究一致，R•自由基是消耗O₂的主要贡献者。 * 路径二（新发现）： GSSG•⁻ + O₂ → GSSG + O₂•⁻。本研究发现，谷胱甘肽二硫化物自由基阴离子（GSSG•⁻）对氧消耗的贡献与R•路径大致相当（各贡献约2个G值单位）。这一路径此前在FLASH背景下未被充分认识。 在典型的FLASH条件（N=20，对应~30 Gy剂量，~10⁶–10⁷ Gy/s剂量率）下，模拟计算出的氧消耗浓度[–O₂]随时间上升，在约1秒时达到一个平台，约为13.6 µM（图4）。这意味着，在初始氧浓度为30 µM的条件下，仅有约45.3%的细胞内氧气被消耗，未能实现完全耗竭。
3. 临界剂量率阈值的识别： 模型预测存在一个临界剂量率阈值，低于此阈值，ROD假说所预测的FLASH效应无法完全显现。如图4所示，只有当每个脉冲的质子数N增加到约50时，模拟的氧消耗浓度才能达到约30 µM，即实现细胞内氧气的完全（100%）瞬时耗竭。而在N=20（代表典型FLASH参数）时，仅能实现部分耗竭。这一结果表明，要达到ROD假说所设想的、足以引起显著急性缺氧的氧耗竭水平，所需的剂量率可能比当前许多FLASH实验中所用的更高。
4. 敏感性分析结果： * 降低RH浓度（从1 M到0.5 M）： 导致氧消耗量减少约10%（例如，N=20时从13.6 µM降至12.3 µM，图5）。这是因为更少的RH导致•OH生成更少的R•，进而减少了通过R•路径的氧消耗。这支持了主要结论，并表明模型结果对RH浓度变化不极度敏感。 * 降低扩散系数（降至纯水的1/100）： 对结果产生了显著影响（图6）。例如，在N=30时，1秒时的氧消耗浓度从20.4 µM大幅降低至10.3 µM（降幅超50%）。这突显了细胞内扩散环境对化学反应动力学的潜在重要影响，并进一步表明，在更真实的扩散限制环境下，实现完全氧耗竭将更加困难。
5. 对ROD假说的综合评估： 综合上述结果，本研究得出结论：在模拟的、接近生理条件的细胞环境中，使用典型的FLASH放疗参数（~30 Gy， ~10⁶–10⁷ Gy/s），仅能导致细胞内氧气部分消耗（~45%），而非完全耗竭。 同时，模型识别出的临界剂量率阈值（N≈50）高于常用FLASH条件。这些计算结果表明，仅凭瞬时的、放射分解引起的氧耗竭（即ROD假说本身）可能不足以完全解释在临床相关剂量和剂量率下观察到的FLASH组织保护效应。 研究认为，ROD可能部分参与了FLASH效应，但必须有其他机制或因素共同作用。
五、 研究结论与价值
本研究的核心结论是：虽然瞬时的放射分解氧耗竭（ROD）可能在FLASH效应中扮演一定角色，但它本身并不是一个充分的解释机制。研究通过精细的蒙特卡洛辐射化学建模，首次在量化层面上揭示了在生理相关的细胞化学成分和FLASH照射条件下，实现完全氧耗竭的挑战性。同时，研究首次明确地将细胞内源性抗氧化剂网络，特别是抗坏血酸和一氧化氮，置于FLASH放射保护机制的中心位置，指出它们通过高效清除过氧自由基，在减轻辐射损伤中发挥了至关重要的作用。
本研究的科学价值在于： 1. 提供了新的机制见解： 将研究焦点从单一的氧耗竭假说，扩展到了氧化应激与抗氧化防御的动态平衡，为理解FLASH效应提供了更全面的化学视角。 2. 提出了可验证的预测： 明确了临界剂量率阈值的存在，以及GSSG•⁻作为重要氧消耗途径的新角色，这些都可以为后续的实验研究提供具体的验证方向。 3. 建立了强大的计算工具： 开发的IONLYS-IRT模型及其扩展的细胞化学反应方案，为未来研究不同辐射品质、剂量率、细胞化学成分（如不同抗氧化剂水平、缺氧条件）下的FLASH效应提供了一个可扩展的、定量的理论框架。 4. 促进了跨学科对话： 连接了辐射化学、计算物理学和放射生物学，强调了在分子反应动力学层面理解复杂生物效应的重要性。
六、 研究亮点
七、 其他有价值内容
论文在引言部分对FLASH放疗的现状、ROD假说的历史与争议、细胞水辐解的特点以及四种关键抗氧化剂（•NO, GSH, AH⁻, TOH）的生物学功能和反应动力学进行了详尽而深入的综述，为不熟悉该交叉领域的读者提供了扎实的背景知识。此外，研究详细列出了所采用的48个化学反应及其速率常数（表1），并讨论了模型简化（如均质假设）的局限性，体现了研究的透明性和可重复性，为其他研究者进行类似或更复杂的模拟奠定了坚实基础。作者也指出，他们正在开发一个更详细的区室化细胞模型，以更好地模拟细胞异质性和不同细胞组分（如DNA与脂质）对辐射损伤响应的差异，这指明了未来研究的方向。