关于分子氢作为抗氧化剂和辐射防护剂的蒙特卡洛模拟研究学术报告
一、 研究作者、机构及发表信息 本研究由加拿大舍布鲁克大学 (Université de Sherbrooke) 医学院医学影像与放射科学系的 Sumaiya Akhter Ria、Jintana Meesungnoen 和 Jean-Paul Jay-Gerin*（通讯作者）合作完成。研究论文《Molecular Hydrogen as an Antioxidant and Radioprotector: Mechanistic Insights from Monte Carlo Radiation-Chemical Simulations》于2025年8月27日发表在期刊 *Antioxidants*（2025年第14卷第1054期）。
二、 学术背景与研究目标 本研究属于辐射化学与放射生物学交叉领域，聚焦于利用计算模拟方法探究辐射防护剂的分子机制。其背景基于一个核心的放射生物学现象：生物体（约70-80%由水构成）在暴露于电离辐射时，细胞内的水分子会发生辐射分解（Radiolysis），产生一系列高活性反应物种，其中羟基自由基（•OH）因其极高的反应活性，被认为是造成DNA、脂质和蛋白质等关键生物分子氧化损伤的主要元凶。因此，清除•OH自由基是辐射防护和抗氧化策略的关键环节。
分子氢（H₂）近年来作为一种潜在的治疗性抗氧化剂受到关注，其被认为具有选择性中和•OH自由基的能力，且具备低毒性、高扩散性等优点。然而，与已知的高效辐射防护剂（如含硫化合物胱胺）相比，H₂在辐射化学层面的清除效率、作用机制及其在复杂生物化学环境中的最终产物的生物学意义，尚缺乏定量化的深入理解。
本研究旨在通过精确的蒙特卡洛径迹化学模拟，从辐射化学基础层面定量评估H₂作为抗氧化剂和辐射防护剂的效能。具体目标包括：1）模拟并比较在有无溶解H₂条件下，水辐解过程中•OH自由基随时间变化的产额；2）研究H₂浓度（从低于生理浓度到超饱和浓度）对其清除效率的影响；3）将H₂的清除效能与经典的辐射防护剂胱胺进行直接对比；4）探讨在含氧与无氧环境下，H₂清除•OH后产生的次级自由基（如H•原子）的归趋及其潜在的生物学意义。
三、 详细研究流程与方法 本研究采用理论计算与模拟方法，核心工具是作者实验室自主开发的“蒙特卡洛径迹化学模拟代码”——IonLYS-IRT。整个研究流程不涉及传统湿实验，而是通过构建物理化学模型并进行大规模随机模拟来获得数据。
1. 模拟体系与条件设定： 研究模拟了中性水溶液在受到300 MeV质子辐照下的辐射化学反应。选择300 MeV质子是因为其线性能量转移（Linear Energy Transfer, LET）值较低（约0.3 keV/µm），能够模拟⁶⁰Co γ射线或快电子等常规低LET辐射的辐射化学效应。模拟环境分为两种：脱气（无氧）条件和充气（含0.25 mM溶解氧，模拟生理条件）条件。在这两种基础环境中，分别添加不同浓度的溶质进行研究：分子氢（H₂）浓度范围为0至10 mM（涵盖并远超其在水中的饱和溶解度~0.78 mM）；作为对比，还模拟了相同浓度范围（0.01至10 mM）的胱胺（Cystamine）溶液。
2. 模拟代码与原理： 本研究使用的IonLYS-IRT代码是一个多层次、多阶段的模拟程序，它将辐射作用的早期物理过程与后续的化学反应动力学相结合。 * IonLYS模块（物理与物理化学阶段）： 该模块以事件驱动的方式，模拟从辐射能量沉积开始，直到约1皮秒（ps）内的早期过程。它追踪质子径迹中产生的所有初始物种（如水合电子e⁻ₐq、H₃O⁺、OH⁻、H•、H₂、•OH、H₂O₂等）的空间分布，这个分布是高度非均匀的，形成所谓的“刺团”（Spur）。 * IRT模块（化学阶段）： 此模块接收IonLYS模块输出的初始物种空间分布作为输入，模拟从1 ps到毫秒（ms）时间尺度的化学反应。它采用“独立反应时间”（Independent Reaction Times, IRT）方法，这是一种计算效率高的随机算法，避免了追踪每个粒子的随机行走轨迹，却能准确模拟物种的扩散及其相互之间、或与溶质（如O₂、H₂、胱胺）之间的化学反应。代码包含了详尽的化学反应网络及其对应的速率常数。对于含胱胺的体系，作者采用了经过先前研究验证的、包含17个额外反应的扩展反应方案（详见表1），该方案曾成功预测Fricke-胱胺剂量计的实验结果。 * 关键参数与假设： 模拟在25°C下进行。考虑了离子强度对离子-离子反应速率常数的影响（e⁻ₐq自复合反应除外）。采用了各物种（包括H₂和O₂）相应的扩散系数。由于溶质浓度远低于水浓度，其直接辐射效应被忽略。模拟通过计算大量（40-100次）质子径迹片段的化学产额并取平均，以获得具有统计可靠性的结果。辐射化学产额以每100 eV吸收能量产生的分子数（G值）报告。
3. 数据分析流程： 模拟的核心输出是不同辐射化学物种（特别是•OH、H•、H₂O₂、O₂•⁻、HO₂•等）的G值随时间（从1 ps到100 ms）演化的曲线。通过分析这些曲线，可以： * 量化清除效率： 观察添加H₂或胱胺后，•OH产额下降的幅度和速度。 * 比较作用机制： 通过对比不同溶质、不同浓度下的曲线，分析清除反应发生的主导时间尺度（是在非均匀的刺团内，还是在均匀的体相溶液中）。 * 评估次级效应： 追踪•OH被清除后，产生的H•原子在含氧/无氧条件下的进一步反应路径及最终产额。 * 进行对比研究： 将H₂与胱胺在相同浓度下对•OH的清除动力学进行直接对比。
四、 主要研究结果 模拟结果以一系列时间演化图（图1-4）清晰呈现，揭示了分子氢作为辐射防护剂的定量机制及其与胱胺的差异。
1. 脱气条件下H₂的作用（图1）： 在无氧水中添加1 mM H₂，导致•OH自由基的产额在1-100 µs时间范围内显著下降约2.5个G单位，并在约100 µs时趋近于零，这直接证实了反应（4）（•OH + H₂ → H• + H₂O）的发生。然而，一个关键发现是：•OH的清除伴随着H•原子产额的显著增加。在没有H₂时，10 ms时G(H•)约为2.8分子/100 eV；而添加1 mM H₂后，G(H•)上升至约5.2分子/100 eV。图1c和d进一步分解了H•原子的来源，明确显示在添加H₂后，反应（4）成为了H•的一个重要新来源。
2. 充气条件下H₂的作用（图2）： 在含氧（0.25 mM O₂）水中，溶解氧会迅速与水合电子和H•原子反应，分别生成超氧阴离子自由基（O₂•⁻）和氢过氧自由基（HO₂•）。添加1 mM H₂后，同样观察到了•OH被完全清除（约100 µs时）。此时，被清除的•OH通过反应（4）转化为H•，而H•又立即被O₂捕获生成HO₂•。因此，在长时间尺度上，体系中剩余的活性氧物种主要是HO₂•/O₂•⁻，其总产额在10 ms时达到约5.52分子/100 eV，与脱气条件下高H•产额（5.2分子/100 eV）相对应。H₂O₂的产额在不同条件下基本保持不变。
3. H₂浓度依赖性与作用时间尺度（图3）： 图3a显示，随着H₂浓度从0.01 mM增加到10 mM，•OH产额开始下降的时间点显著提前，从约60 µs（0.01 mM）提前到约30 ns（10 mM）。图3b进一步展示了在固定观察时间点（如0.1, 1, 10, 100 µs），•OH产额随H₂浓度增加而下降的曲线。这些结果定量证明了H₂对•OH的清除是浓度依赖性的。更重要的是，分析表明，在整个研究的浓度范围内，H₂对•OH的清除主要发生在辐射径迹消散后的“均匀化学阶段”。其清除的特征时间尺度可由伪一级速率常数倒数1/(k[H₂])估算，范围从0.01 mM时的~2.5 ms到10 mM时的~2.5 µs。
4. 与胱胺的效能对比（图4）： 图4展示了胱胺清除•OH的惊人效率。在10 mM浓度下，胱胺在约100 ps时就开始急剧降低•OH产额，并在约30 ns内将其降至接近零。相比之下，相同浓度的H₂在约30 ns时才开始降低•OH产额，直到约15 µs才完成清除。即使在低得多的浓度下（如0.01 mM），胱胺启动清除的时间（~100 ns）也远早于H₂。这表明胱胺的清除反应主要发生在辐射作用的“非均匀化学阶段”（即刺团内部），因其反应速率常数极高（与•OH的反应速率常数达1.7 × 10¹⁰ M⁻¹ s⁻¹），能更早、更快地拦截•OH。而H₂的反应速率常数（4 × 10⁷ M⁻¹ s⁻¹）较低，因此其作用主要体现在后期的均匀体相中。
5. 次级产物的生物学启示： 研究明确指出，H₂在清除•OH的同时，导致了H•原子产额的增加。在充气条件下，这些H•原子迅速被O₂转化为活性较低的HO₂•/O₂•⁻，后者可被内源性超氧化物歧化酶（SOD）进一步处理。然而，在缺氧环境中（如实体瘤内部），H•原子可能因缺乏O₂而积累，并参与其他化学反应（如还原生物分子或与金属离子相互作用），这可能带来复杂的、非单纯的抗氧化效应。
五、 研究结论与意义 本研究通过高保真的蒙特卡洛模拟，从辐射化学基本原理层面得出了以下核心结论： 1. 分子氢是一种有效的选择性•OH清除剂： 模拟定量证实，H₂能显著且浓度依赖性地降低水辐解产生的•OH自由基产额，即使在远低于饱和浓度的水平（0.01 mM）也有效。 2. H₂的清除效率低于经典辐射防护剂胱胺： 胱胺凭借其更高的化学反应活性，能在更早的时间尺度（刺团内）以更快的速度清除•OH，显示出更优的辐射化学防护效能。 3. H₂的作用伴随次级自由基生成： H₂清除•OH的化学产物是H•原子。在含氧环境中，H•被转化为HO₂•/O₂•⁻；在缺氧环境中，H•可能持续存在并引发其他反应。这意味着H₂的辐射防护作用并非简单地“消除”活性物种，而是将其转化为活性较低的形态，其最终生物学效应取决于细胞的氧化还原状态和抗氧化酶系。 4. H₂的突出优势在于其生物安全性： 尽管化学效能不及胱胺，但H₂具有无毒性、高生物膜渗透性、选择性（不干扰重要的信号分子如O₂•⁻、H₂O₂、•NO）以及已报道的抗炎作用。这些特性使其在需要长期或大剂量使用的临床场景（如放疗中的正常组织保护、慢性氧化应激疾病管理）中具有独特的应用潜力。
本研究的科学价值在于，它提供了一个定量的、机制性的框架来理解H₂作为辐射防护剂的作用。它超越了简单的“清除•OH”概念，揭示了伴随的化学反应链及其对最终生物效应的影响。应用价值在于为优化H₂在放射治疗（尤其是在含氧正常组织与缺氧肿瘤组织并存的复杂环境中）的应用策略提供了理论依据。
六、 研究亮点 1. 方法创新性： 成功应用自主开发的、经过广泛验证的蒙特卡洛径迹化学模拟代码（IonLYS-IRT），对H₂和胱胺的辐射防护机制进行了从皮秒到毫秒的全时间尺度、高空间分辨率的定量模拟。这是一种“计算实验”，能够揭示传统实验难以观测的早期瞬态过程和空间效应。 2. 对比研究的深度： 不仅研究了H₂，还将其与金标准辐射防护剂胱胺在相同模拟条件下进行直接、系统的对比，清晰揭示了二者在作用动力学和效率上的本质差异。 3. 对次级效应的关注： 研究没有停留在“清除•OH”这一初级结论，而是深入追踪并强调了清除反应产物（H•原子）的归趋，并讨论了其在缺氧条件下的潜在生物学意义，体现了对辐射化学复杂性的深刻认识。 4. 超越生理浓度的探索： 模拟了高达10 mM的H₂浓度，这虽然远超其物理溶解度，但从理论计算角度探索了极限情况下的化学趋势，有助于完整理解浓度-效应关系。
七、 其他有价值内容与未来展望 论文在讨论部分提出了两个明确的未来研究方向： 1. 高LET辐射的影响： 计划将模拟扩展到高LET辐射（如降低质子能量至150 keV，LET升至~72.2 keV/µm）。预计在高LET条件下，径迹内•OH浓度更高，其自身复合反应（•OH + •OH, •OH + H•）的竞争将增强，可能会降低H₂对•OH的净清除效率。 2. 高剂量率（FLASH放疗）效应： 计划研究剂量率的影响。在高剂量率下，体相中•OH浓度升高，可能促进•OH的相互复合，同样可能影响H₂的防护效果。
这些展望表明，本研究建立的模拟平台具有很好的可扩展性，可用于探索辐射质量（LET）和剂量率等关键参数如何影响抗氧化剂和辐射防护剂的效能，这对于发展新型放疗策略（如FLASH放疗）具有重要意义。此外，作者也指出，目前缺乏在含H₂条件下水辐解产物随时间变化的实验产额数据来直接验证模拟结果，这为未来的实验研究指明了方向。