使用DBSCAN算法预测质子与碳离子造成的DNA损伤产额并评估束流质量

更锐利的癌症束流

现代癌症治疗越来越依赖带电粒子束（如质子和碳离子）来攻击肿瘤，同时尽量保护周围的健康组织。然而相同的物理剂量并不总会造成相同程度的生物学损伤。本文探讨一个实用问题：能否用比当前复杂计算更简单、更快速的方法预测给定粒子束对DNA的“破坏力”？

为何DNA断裂至关重要

当辐射穿过细胞时，会在水和DNA中留下微小的能量沉积。这些事件可能会切断DNA双螺旋的一条或两条链。单链断裂通常可被修复，而双链断裂——尤其是聚集在一起时——更可能导致细胞死亡或产生突变。临床上主要使用物理量（如剂量和线性能量传递率 LET）来制定治疗计划，但这些量无法完全解释严重DNA损伤发生的频率。若能建立束流性质与DNA断裂之间更直接的联系，就能帮助设计既更有效杀伤肿瘤又对患者更安全的粒子治疗方案。

将微小打击聚为有意义的损伤

作者借鉴数据科学中的一个想法：聚类分析。他们不再模拟辐射击中水后每一步的化学过程，而只模拟质子和碳离子在液态水中留下的初始“轨迹结构”能量沉积。然后采用广泛使用的聚类算法DBSCAN来识别损伤点群。任何沉积至少17.5电子伏特的相互作用都被计为潜在的链断裂；若有至少两个此类点相距约2.1纳米——与DNA直径相近——则被归为同一簇，解释为双链断裂。孤立的点则视为单链断裂。通过调整该距离使模型重现详细的基准模拟，研究团队将原始轨迹转化为对简单与复杂DNA损伤产额的估计。

为束流打分的新途径

基于聚类结果，作者引入一项称为束流质量（Quality of Beam，QoB）的新指标：每个粒子每微米路径产生的簇数。将其按粒子沿路径沉积的能量归一化，得到以“每单位能量簇数”为单位的量。对于能量范围在0.5到200兆电子伏特的治疗质子，这一归一化QoB与一个可信且更为复杂模型预测的双链断裂数量呈显著的线性关系。也就是说，一个简单的换算因子可以将归一化QoB直接转换为双链和单链断裂产额，从而绕过完整的水放射解离模拟，同时与既有研究保持一致。

比较质子与碳离子

相同框架也用于碳离子——其轨迹更致密，在一些专门的癌症中心被采用。使用针对质子优化的参数，模型在一定LET范围内（约160–200千电子伏特/微米）仍发现归一化QoB与碳离子双链断裂之间的紧密线性关系。超过该范围后，趋势出现拐弯：额外的能量不再持续增加新簇的数量，这就是所谓的“过杀”效应。在这种情况下，大量能量集中到已损伤区域，额外的电离对生物学影响贡献甚微。值得注意的是，质子与碳离子的归一化QoB对LET的曲线在细胞相对生物效应（RBE）的已发表测量中表现出相似特征——上升、一个较宽的峰值，以及在极高LET处的下降——而传统仅用LET无法充分捕捉这些变化。

这对未来治疗的意义

对非专业读者而言，核心信息是：相同物理强度的辐射并不总以相同方式伤害细胞。关键在于能量在纳米尺度上怎样分布于DNA周围。本研究表明，通过将辐射轨迹视为数据点并用聚类算法分组，可以快速估算严重DNA断裂发生的频率，并定义一个更能反映生物学影响的束流“质量”指标。对质子而言，该方法可用单一因子直接预测单链与双链断裂产额；对更重的离子仍需一定调整，但同样的方法能突出诸如过杀效应等重要现象。长远来看，这类具生物学信息的束流度量有望帮助细化粒子治疗计划，将致瘤杀伤力精确集中到需要的位置，同时减少对健康组织的非预期伤害。

引用: Chaibura, S., Liamsuwan, T. Predicting DNA damage yields and assessing beam quality for protons and carbon ions using a DBSCAN algorithm. Sci Rep 16, 10327 (2026). https://doi.org/10.1038/s41598-026-40571-x

关键词: 质子治疗, 碳离子治疗, DNA损伤, 辐射质量, 聚类分析