基于第一性原理计算的225Ac-DOTA螯合放射性药物在衰变途径中的结构与性质研究
一、 研究作者、机构与发表信息
本研究由 Artem R. Khabibullin, Aleksandra Karolak, Mikalai M. Budzevich, Mark L. McLaughlin, David L. Morse 和 Lilia M. Woods 共同完成。研究团队主要来自美国南佛罗里达大学（University of South Florida）等机构。该研究于2018年在线发表在皇家化学会（Royal Society of Chemistry）旗下的期刊 *MedChemComm*（*Medicinal Chemistry Communications*）上，作为一期关于“抗生素耐药性”的主题特刊的一部分。文章最终版本在线发布于2018年6月6日。
二、 学术背景与研究目的
本研究属于核医学、放射性药物化学与计算材料科学的交叉领域，具体聚焦于靶向α粒子治疗（Targeted Alpha-particle Therapy, TAT）中核心药物载体的稳定性问题。
研究背景： 靶向α粒子疗法是一种新兴的癌症治疗策略，其核心是将发射高能α粒子的放射性核素通过靶向分子（如抗体或肽段）精准递送至肿瘤细胞。α粒子具有高线性能量转移（Linear Energy Transfer, LET）和短射程（约40-60微米）的特点，能高效引发DNA双链断裂，同时对周围健康组织损伤较小。放射性核素锕-225（225Ac）因其约10天的半衰期、易于被DOTA（1,4,7,10-四氮杂环十二烷-1,4,7,10-四乙酸）螯合剂稳定结合，并能衰变产生连续4个α粒子的“纳米发生器”效应，被视为TAT的理想候选核素之一。
关键问题与挑战： 然而，225Ac的复杂衰变链（225Ac → 221Fr → 217At → 213Bi → …）带来了一个核心科学问题：在每次α衰变事件中，产生的巨大反冲能（约100-130 keV）是否会破坏DOTA螯合剂的结构，导致子代核素（如221Fr, 217At, 213Bi）从螯合物中释放？这种释放可能导致子代核素非靶向分布，引起肾脏等器官的毒性，并降低肿瘤局部累积的α粒子剂量，从而影响疗效和安全性。尽管有初步实验观察（如Schwartz等人的研究）提示213Bi可能在体内从225Ac-DOTA-抗体中释放，但由于实验方法的局限性和复杂性，这一现象尚未得到明确证实，其背后的物理化学机制更是一片空白。传统的化学结合能（约几十电子伏特）远低于衰变反冲能（约十万电子伏特），但从量子力学角度看，反冲能是否必然导致化学键断裂和核素位移尚不确定。
研究目的： 为了从根本上理解α衰变对225Ac-DOTA靶向放射性药物稳定性的影响，本研究旨在从第一性原理出发，系统研究DOTA与225Ac及其衰变链初始子代核素（221Fr, 217At, 213Bi）形成的螯合物的结构、稳定性及电子性质。同时，研究也包含了作为对照的、已广泛用于磁共振成像（MRI）对比剂的68Ga-DOTA螯合物。研究目标是通过计算模拟，揭示这些螯合物的结合能、键长、电子局域函数、与水分子相互作用等关键性质，评估其在衰变过程中的潜在稳定性，并为未来的实验研究提供理论指导和解释框架。
三、 研究详细工作流程
本研究完全采用计算模拟方法，基于密度泛函理论（Density Functional Theory, DFT）进行。整个工作流程可分为以下几个核心步骤：
1. 模型构建与初始优化： * 研究对象： 构建了五个DOTA金属螯合物模型：68Ga-DOTA、213Bi-DOTA、217At-DOTA、221Fr-DOTA和225Ac-DOTA。其中，213Bi-DOTA的初始晶体结构参考自文献（Csajbók et al., 2003）。对于其他螯合物，通过在Avogadro分子建模与可视化软件中替换中心金属原子，并随后进行初步的几何优化来获得初始结构。 * 研究扩展： 为了模拟生理水环境，对每个螯合物模型额外添加了一个水分子（H2O），构建了X-DOTA-H2O（X代表金属核素）体系，以研究水分子配位对结构稳定性的影响。
2. 第一性原理计算设置： * 计算软件： 所有DFT计算均在维也纳从头算模拟软件包（Vienna Ab-initio Simulation Package, VASP）中完成。 * 计算参数： 采用投影缀加平面波（Projector Augmented Wave, PAW）方法，交换关联泛函使用PBE（Perdew-Burke-Ernzerhof）。每个分子被放置在一个足够大的立方晶胞中心进行单k点计算。能量和力的收敛标准分别设置为10^-4 eV和10^-5 eV/Å。 * 关键物理效应考量： 鉴于所研究金属元素（尤其是Ac、Fr、At、Bi）质量大，研究特别考察了两种关键效应对结果的影响： * 自旋轨道耦合（Spin-Orbit Coupling, SOC）： 通过VASP中的非共线磁性计算来实现，以考虑重元素的相对论效应。 * 范德华（van der Waals, vdW）色散相互作用： 采用vdW-DF2泛函（Langreth-Lundqvist方案）进行修正，以更准确地描述分子间的弱相互作用。 * 电子性质分析： 计算了电子局域化函数（Electron Localization Function, ELF）和电荷密度差（Charge Density Difference），以直观分析化学键类型和电荷转移情况。
3. 计算与分析方法： * 结合能计算： 计算了金属核素与DOTA配体的结合能（Eb），公式为 Eb = E_tot - (E_DOTA + E_X)，其中E_tot为螯合物总能量，E_DOTA为自由DOTA配体能量，E_X为自由金属原子的能量。结合能负值越大，表示螯合物越稳定。 * 结构参数分析： 测量并分析了金属原子与最近邻氧原子（X-O）和氮原子（X-N）的平均键长。同时，对于X-DOTA-H2O体系，定义了多个关键几何参数（如图3所示）：金属到O4/N4平面的距离（dX-O， dX-N）、O4与N4平面间距（dN-O）、金属与水分子距离（dX-H2O）、开口角（φ， 金属与O4平面夹角）和扭转角（ω， O4与N4平面的相对扭转角），用以量化配位结构。 * 水分子结合能： 计算了水分子与螯合物的结合能（E_bind, H2O），公式为 E_bind,H2O = E_tot+H2O - (E_X-DOTA + E_H2O)。 * 反冲能估算： 基于经典非相对论运动学，利用美国国家核数据中心（National Nuclear Data Center）提供的α粒子动能（Qα）数据，估算了子代核素（221Fr, 217At, 213Bi）在α衰变中获得的动能（反冲能，E_rec），公式为 E_rec = Qα * m_α / (m_α + m_daughter)。
四、 主要研究结果
1. DOTA螯合物的结构与稳定性： * 结合能与键长： 计算结果表明（见表1），不考虑SOC和vdW时，225Ac-DOTA和68Ga-DOTA的结合能最大（最负），最为稳定；而221Fr-DOTA的结合能最小，最不稳定。引入SOC和vdW修正后，稳定性排序不变，但具体数值发生变化。例如，SOC使221Fr-DOTA的结合能增加了0.90 eV（更稳定），但使其他几个螯合物的结合能绝对值略有减小。vdW相互作用显著增强了68Ga-DOTA的稳定性，但轻微削弱了213Bi-DOTA、217At-DOTA和225Ac-DOTA的稳定性。 * 键长分析： 平均键长显示，221Fr-DOTA的Fr-O和Fr-N键长均超过3 Å，显著长于其他螯合物，表明结合非常松散。217At-DOTA的键长则相对较短。与已知离子半径和（RIS）比较发现，221Fr-DOTA和217At-DOTA的键长与RIS预期值偏差较大，进一步提示其结构脆弱。 * 化学键性质： 通过电子局域化函数（ELF）和电荷密度差分析（图1，图2），揭示了化学键的本质。68Ga-DOTA、213Bi-DOTA和225Ac-DOTA主要呈现离子键特征，金属与配体原子间存在明显的电荷转移（电子从金属向N、O原子转移）。217At-DOTA显示出一定的共价-离子混合键特征。而221Fr-DOTA的ELF和电荷密度差图显示，Fr原子与周围O、N原子之间几乎没有电子共享或明显的电荷转移区域，表明其与DOTA之间主要是极弱的分子间作用，而非真正的化学键。这从电子结构层面证实了221Fr-DOTA极不稳定的预测。
2. 水分子配位的影响： * 结构稳定性： 在X-DOTA-H2O体系的弛豫计算中，221Fr-DOTA-H2O发生了结构崩解（图4c,d），无法保持完整构型，这直接证明了221Fr-DOTA在含水环境中极不稳定。 * 结构参数与配位模式： 对于其他稳定的体系（表2），计算的结构参数（如键长、角度）与已有68Ga-DOTA的实验数据大体吻合，验证了计算方法的可靠性。特别重要的是开口角（φ）和扭转角（ω）。研究发现，68Ga-DOTA-H2O和225Ac-DOTA-H2O在考虑vdW作用后，其φ角分别达到143.7°和148.3°，均大于文献中提出的镧系元素DOTA螯合物发生水分子配位的阈值（~135°）。同时，它们与水分子的结合能（E_bind, H2O）也相对较强（负值更大）。这表明它们处于扭曲四方反棱柱（Twisted Square-Antiprismatic, TSA）构型，有利于水分子配位和快速水交换，这对于MRI对比剂功能至关重要。 * 213Bi-DOTA-H2O的特殊性： 213Bi-DOTA-H2O的φ角仅为116.1°，远低于135°阈值，且Bi原子位于DOTA空腔深处，与水分子的距离（dX-H2O）较远（~3.79 Å）。其水分子结合能几乎为零（0.03 eV）。电荷密度差分析（图6b）也显示Bi-DOTA与水分子之间几乎没有电荷转移。这些结果一致表明，213Bi-DOTA不与内层水分子配位，水溶性差，这与之前的实验观察相符。 * 217At-DOTA-H2O： 其水分子结合能尚可，但φ角（123.3°）仍低于阈值，且水分子以氢原子朝向DOTA的O4平面，主要通过偶极-偶极相互作用结合，配位作用较弱。
3. 反冲能估算： 基于经典运动学，计算得出221Fr、217At和213Bi在α衰变中获得的动能（反冲能）分别为103.64 keV、114.78 keV和130.30 keV。这个能量量级（~100 keV）远高于DOTA与金属核素之间的化学结合能（~几十eV）。然而，文章特别强调，α衰变是一个量子力学过程，不能简单套用经典力学模型。反冲能量很可能主要以子核激发能的形式耗散，而非全部转化为破坏化学键的动能。因此，子核是否一定会从螯合剂中逃逸，以及螯合剂本身是否会被反冲破坏，仍是一个开放性问题，需要更深入的研究。
五、 研究结论与价值
结论： 本研究通过系统的第一性原理计算，首次在原子尺度上深入揭示了225Ac衰变链中关键DOTA螯合物的结构、稳定性及与水相互作用的性质。 1. 稳定性排序： 在研究的螯合物中，225Ac-DOTA本身具有最强的结合能，是最稳定的复合物。其子代产物中，221Fr-DOTA（无论是否含水）是最不稳定的，其结合能最低，键长过长，且电子结构分析表明其与DOTA之间几乎不存在有效的化学键，在含水环境中结构会崩溃。 2. 水合行为差异： 225Ac-DOTA和68Ga-DOTA倾向于形成TSA构型，开口角大，能与水分子有效配位。而213Bi-DOTA不与内层水配位，水溶性差。217At-DOTA与水的相互作用较弱。 3. 对衰变链稳定性的启示： 研究结果表明，在225Ac的衰变链中，221Fr的释放风险最高。一旦225Ac发生衰变生成221Fr，由于221Fr-DOTA极不稳定，Fr核素很可能从螯合剂中解离。而后续的217At和213Bi，其DOTA螯合物的稳定性相对较高，如果它们在衰变事件后仍能停留在DOTA配位场内（尽管存在反冲能的不确定性），则可能被部分保留。
科学价值与应用意义： * 理论指导价值： 该研究为理解α发射体放射性药物在衰变过程中的化学行为提供了重要的量子力学视角和定量数据。它超越了简单的能量比较，从电子结构层面解释了不同核素-螯合剂稳定性的差异。 * 指导实验设计： 研究结果明确指出221Fr是衰变链中的“薄弱环节”，这为后续的实验研究（如体外稳定性实验、体内生物分布研究）提供了明确的关注点和解释依据。例如，在评估225Ac-DOTA药物体内分布时，应重点关注221Fr和213Bi的子体逃逸情况。 * 药物设计优化： 研究提示，若要提高225Ac靶向药物的体内稳定性，可能需要开发对Fr有更强螯合能力的新型配体，或采用能更好地束缚衰变子体的递送策略（如纳米载体包封）。 * 跨领域方法学应用： 本研究展示了将计算材料科学方法（DFT， 包括SOC和vdW修正）应用于核医学和放射性药物化学问题的可行性，为相关领域的研究提供了方法学范例。
六、 研究亮点
七、 其他有价值内容
研究还对作为对照的68Ga-DOTA进行了计算，其结果与大量实验数据吻合良好，这间接验证了本研究采用的计算方法和参数的可靠性。此外，文章在讨论部分对先前Schwartz等人的实验发现进行了批判性评述，指出了其数据解读中可能存在的局限性（如未考虑半衰期差异对测量比例的影响、未能区分循环中与肿瘤内药物的差异等），体现了严谨的学术态度。最后，作者团队声明无利益冲突，并感谢了美国能源部的经费支持（DE FG02-06ER46297）。