基于经典细胞动力学及细胞变化假设的放疗期间放射生物学参数与再增殖剂量当量动态估算模型

《Frontiers in Oncology》：Estimation of radiobiological parameters and repopulation dose equivalents using classical cell kinetics with hypotheses concerning cellular changes during radiotherapy

　　

目的：研究肿瘤反应及放射生物学参数变化，以更好指导不同组织学类型肿瘤放疗中总治疗时间问题和治疗延迟的决策。方法：研究人员将经典细胞动力学方程与放射生物学原则相结合，在包含非增殖性放射抵抗缺氧分数、再氧合以及由肿瘤体积倍增时间（tumour volume doubling time, TD）控制的细胞丢失因子（cell loss factor, φ）下降的模型中，研究放疗期间有效倍增时间（effective doubling time, T_eff）和再增殖剂量当量（repopulation dose equivalent, K_repop）的变化，使得较长TD者具有较长的加速再增殖（accelerated repopulation, AR）表观起始时间。结果：该模型发现，在组织学类型明确且具有典型细胞动力学特征的肿瘤中，T_eff缩短的时间与加速再增殖起始时间相一致，并伴随肿瘤缺氧程度随时间下降。K_repop值在治疗期间呈稳定增加，并趋近于各组织学类别治疗中断评估中通常使用的恒定K_repop值。初始φ较低的肿瘤类别显示缺氧更持久；初始TD很长的肿瘤类别在放疗期间或紧接放疗后不发生加速再增殖，其K_repop值保持低水平。结论：目前仅在关键AR起始时间之后使用恒定K_repop值，在治疗中断补偿中可能为次优，因为T_eff进行性缩短使K_repop可能随时间连续变化。对于较短疗程，使用高K值（如较长疗程所用）可能导致补偿性过量照射；应优先采用较低K_repop值获得的较低补偿剂量，且无需等待AR起始的滞后时间。非常长的中断可能需要与实施补偿时相对应的最高K_repop值。

研究背景方面，放疗总治疗时间延长和治疗中断可影响肿瘤控制。生物有效剂量（biological effective dose, BED）模型及每日补偿剂量K_repop已用于中断补偿，但K_repop常被假定为恒定值，并通常在AR起始后才设滞后时间，这可能过度简化AR过程。Fowler提出细胞丢失因子φ在放疗中可能下降，使T_eff缩短并解释AR；Jones and Dale进一步建模φ连续变化，但需对φ指数下降进行积分，方程应用较复杂。为此，研究人员提出一种替代性迭代肿瘤控制模型，纳入非增殖性缺氧细胞、增殖性含氧细胞、缺氧向含氧区的再氧合通量以及由TD控制的φ下降，以估算不同组织学肿瘤放疗期间T_eff与K_repop的动态变化。该研究得出结论：K_repop可能随时间连续变化，恒定K_repop及固定AR起始滞后用于中断补偿可能不理想，对总治疗时间和治疗延迟决策具有意义。

关键技术方法概括：研究人员将经典细胞动力学方程与放射生物学原则结合，采用线性二次模型（linear quadratic model, LQ）计算含氧与缺氧细胞存活分数，缺氧细胞按固定氧增强修正因子q_α、q_β调整α和β；以k_L通量模拟再氧合，φ按指数下降，速率参数z=1/TD，并在3×细胞周期时间（TC）后启动φ下降；T_eff按公式随φ缩短，肿瘤控制概率（tumour control probability, TCP）采用Poisson e^?N计算；使用Mathematica编程迭代，初始克隆原细胞数设为5×10⁸。该研究未使用临床样本队列，参数来自文献汇编并按组织学类型作示意性选择。

鳞状细胞癌模拟结果：在TD约38天、2 Gy×5次/周条件下，T_eff有效缩短，K_repop在20—30天上升至约0.3 Gy·day^-1并趋近0.6—0.9 Gy·day^-1，缺氧分数先升后降，TCP曲线陡峭；结论与Withers分析及SCC的AR预期相容。腺癌模拟结果：较快生长者TD=75天，约40—50天后K_repop显著增加并在最长疗程约0.3 Gy·day^-1；慢生长者TD>270天，70天内K_repop无显著增加且T_eff仍>50天。基底细胞癌模拟结果：φ₀=98%导致TD和T_eff很长，K_repop低，30天约0.03、55天约0.05 Gy·day^-1，70天T_eff约50天而初始>350天。胶质母细胞瘤模拟结果：T_pot=10天、φ₀约66%、TD约30天时，模型显示缺氧显著并持续。肉瘤模拟结果：快生长与慢生长肉瘤中缺氧随时间下降但仍保持较低显著水平；快生长者K_repop显著，慢生长者不显著。低分次鳞状细胞癌模拟结果：7 Gy×3次/周、49 Gy疗程末TCP明显增加，缺氧约第6天急剧下降，T_eff下降但14天K_repop仅0.16 Gy·day^-1，仍低。K_repop与TD关系模拟结果：TD=20天和130天者在70天预测K_repop分别为1.1和0.1 Gy·day^-1，且不依赖剂量分割方案。三维参数敏感性模拟结果：K_repop随φ₀和TD变化显示敏感性，部分极端值不现实但多数相关，结果不依赖剂量分割方案。

讨论部分总结：该模型与已知细胞动力学理论一致，TD被假设决定φ变化速率z及T_eff缩短，φ₀影响再氧合速率。模型局限包括依赖LQ、固定氧修正因子、无患者异质性、缺氧区完全不增殖但仍参与再氧合、固定克隆原细胞数、缺乏显式血管动力学及参数不确定性。T_pot与SF₂单独预后相关性不足，因其未纳入TD、缺氧分数和再氧合；BCC与胶质母细胞瘤T_pot可能相近，但TD可因φ₀差异约10倍。Tannock及Jones and Camplejohn的结果提示远离血管处T_pot和生长分数降低、缺氧区细胞不增殖；k_L对应细胞丢失半衰期约数天，SCC初始约4.2天、80天约11天，Warburg效应可解释较长半衰期。K_repop应随补偿时点变化，短疗程低分次中断宜考虑较低K_repop，长中断需与补偿时间相应的最高K_repop；α由0.3变为0.4 Gy^-1可使K_repop降低25%。研究结论部分翻译：该模型提供了关于肿瘤缺氧、再氧合、放射敏感性和最终增殖之间动态关系的重要见解，其时间序列由治疗前TD控制，而TD本身与φ及其变化率相关；AR的显现只能在φ充分降低并由再氧合促进后发生。过去过度强调孤立的增殖标志物，尽管其相关，但仍需了解φ及其时间相关变化，TD在此方面可能至关重要。由于该模型为确定性模型，未来版本可纳入随机参数异质性效应；将遗传突变数据与重要细胞动力学和放射敏感性参数连接可能带来进一步进展。总体结论是K_repop值可能随治疗持续时间变化，这对非计划治疗中断补偿尤其低分次方案具有一定启示。