Pt(II)取代速率何时可信？基于DFT的四个可靠性诊断框架及其在顺铂取代动力学中的应用

《International Journal of Molecular Sciences》：When Can a Computed Pt(II) Substitution Rate Be Trusted? A DFT Framework with Four Reliability Diagnostics Applied to Cisplatin

　　

本研究提出了一种用于Pt(II)配合物取代动力学（substitution kinetics）的密度泛函理论（DFT）计算方案，并以顺铂（cisplatin）作为案例研究进行应用。研究人员旨在证明，一个与实验相符的计算势垒本身并不能验证动力学模型，且需要进一步诊断以评估结果的可靠性。研究采用M06-2X/x2c-TZVPall计算，结合隐式溶剂化（implicit solvation）、显式处理一个配位水分子（coordinating water molecule）以及一系列物理修正（physical corrections），研究了水合（aquation）、硫醇化（thiolation）（采用封端的半胱氨酸硫醇盐模型capped cysteine-thiolate model）和鸟嘌呤取代（guanine substitution）（通过9-甲基鸟嘌呤9-methylguanine）路径。第一次水合的预复合体到过渡态（precomplex-to-transition state）势垒比从实验推断的Eyring势垒低2.0 kcal?mol^?1；应用对称因子σ = 2后，计算得到的观测速率常数比实验值高56倍。硫醇化的形式低浓度二级系数（formal low-concentration second-order coefficient）比已报道的谷胱甘肽速率系数（glutathione rate coefficient）高1.5–1.8个数量级，而省略物理修正会使计算速率预测偏移约5–8个数量级。研究引入了一个直接选择性描述符（direct selectivity descriptor, S_direct），表明其由预复合项（precomplex term）而非活化步骤控制，使其易受隐式溶剂对预复合体过度稳定化（precomplex overstabilization by the implicit solvent）的影响。研究提出了四个应伴随任何计算速率常数或选择性的诊断标准：预复合自由能（precomplex free energy）的符号、K_pre的大小、接近扩散极限（diffusion limit）的程度，以及将选择性分解为预复合和活化贡献。这些为评估计算预测的可靠性提供了系统基础。

顺铂（cisplatin）于1844年由Michele Peyrone首次合成，但其生物活性直至1965年才由Rosenberg及其同事在E. coli中偶然发现。研究人员阐明，顺铂的抗癌机制主要依赖水合形式与DNA的相互作用，特别是与鸟嘌呤N7位结合形成1,2-链内交联（cross-links），这些DNA损伤通过p53信号通路和caspase级联激活诱导细胞周期阻滞与细胞凋亡。后续研究认识到，水合铂物种还可与细胞内强亲核试剂——尤其是含硫生物分子如谷胱甘肽（glutathione）、金属硫蛋白（metallothioneins）及其他半胱氨酸富集蛋白发生作用。此外，硫醚给体如甲硫氨酸（methionine）能够可逆地配位至铂中心，铜转运蛋白CTR1的甲硫氨酸富集基序也被证实可与顺铂直接相互作用并形成混合配体铂-鸟嘌呤加合物。值得注意的是，Pt(II)配合物可直接与硫醇亲核试剂反应，无需预先水合激活，从而提供一条不依赖激活的铂捕获途径，降低可用于DNA铂化的铂比例。这些铂-硫醇加合物不仅参与后续鸟嘌呤导向转化，更反复与铂类药物解毒和耐药机制相关联，成为临床肿瘤学的主要障碍。

在分子尺度上，计算化学家长期致力于Pt(II)配合物研究，建立了显式溶剂化（explicit solvation）和隐式溶剂化（implicit solvation）对水解势垒解释的重要性等理论框架。然而，既有研究多聚焦于单一步骤的反应能与活化势垒，而对相遇复合体（encounter complex）形成、标准态效应（standard-state effects）及扩散对实验观测速率的影响关注不足。尽管水合与后续取代均在动力学控制下进行，但计算势垒与实验相符是否等同于模型物理正确，这一核心问题尚未得到系统回答。研究人员指出，这种相符可能源于误差抵消（cancellation of errors），特别是在省略解释缔合物理修正的情况下。因此，以顺铂为概念验证，研究人员开展本研究旨在证明：计算势垒与实验一致本身并不能验证动力学模型，必须引入进一步诊断以评估计算结果的可靠性。该研究为计算Pt(II)取代动力学提供了标准化协议，并建立了系统性的可靠性评估框架，对避免虚假模型验证具有重要意义。

研究人员采用ORCA 6.1.0软件包开展量子化学计算。全局最小值结构采样通过GFN2-xTB方法和ALPB隐式水模型结合GOAT构象搜索模块完成；孤立结构、预复合体（precomplex）、后复合体（postcomplex）及过渡态（transition state）的优化在M06-2X泛函、精确两分量相对论哈密顿量、x2c-TZVPall基组及SMD隐式水溶液模型下进行，并采用DEFGRID3数值积分网格与TightSCF收敛标准。热化学评估在同一理论水平下于298.15 K完成，所有过渡态均通过内禀反应坐标（IRC）计算验证连接预期的预复合体与后复合体。预复合体、过渡态和后复合体结构均以单个配位水分子氢键连接至离去基团构建，以捕获特异性第一层相互作用。关联熵通过解析方式处理，采用标准态变化（S_SS）、溶剂浓度变化（W_SCA）和Benson-Okuno修正（Benson–Okuno correction）三项物理修正。半胱氨酸建模为乙酰半胱氨酰-N-甲基酰胺硫醇盐（Ac-Cys-NHMe^?），鸟嘌呤采用中性9-甲基鸟嘌呤（9-methylguanine）。速率常数通过常规过渡态理论（TST）计算，并经由Collins-Kimball方程处理扩散极限效应；接触离子对缔合常数以Fuoss方程（Fuoss equation）独立估算。研究未涉及生物样本队列来源。

2.1. 反应路径：Cisplatin与Water、Thiolate和Guanine的反应性

研究人员计算了顺铂第一次水合、硫醇盐（thiolate）和9-甲基鸟嘌呤对非水合及单水合顺铂的进攻路径。第一次水合的预复合体到过渡态势垒为22.0 kcal?mol^?1，较由实验速率常数经Eyring方程推断的活化自由能24.0 kcal?mol^?1低2.0 kcal?mol^?1。引入对称因子σ = 2后，计算观测速率常数为9.06 × 10^?4 s^?1，较实验值高约56倍。非水合态下，硫醇盐预复合体稳定化达7.3 kcal?mol^?1（K_pre = 2.4 × 10⁵ M^?1），而鸟嘌呤预复合体仅稳定化1.3 kcal?mol^?1（K_pre = 9.2 M^?1），直接选择性描述符S_direct为3.7，且该偏好由相遇复合体而非本征活化步骤主导。水合后，硫醇盐取代偏好进一步增强：单水合态预复合体到过渡态势垒降至12.6 kcal?mol^?1，整体势垒降至4.3 kcal?mol^?1，计算表观二级速率常数k_app接近扩散极限（3.13 × 10⁹ M^?1?s^?1），而鸟嘌呤因预复合体去稳定化（ΔG = 1.0 kcal?mol^?1），整体势垒为13.8 kcal?mol^?1，S_direct升至6.8。硫醇化的形式低浓度二级系数k_app × f_active ≈ 0.87 M^?1 s^?1，较已报道的谷胱甘肽速率系数高1.5–1.8个数量级。温度校正至310 K后，该高估仍不少于约0.2个数量级，表明温度差异不能解释此偏差。此外，文献报道单水合顺铂的酸解离常数pK_a = 6.56（37°C），在pH 7.4时仅约13%以单水合形式存在，但计算仅纳入单水合形态。

2.2. 协议设置敏感性分析

研究表明，移除三项物理修正使水合计算速率降低5.3个数量级、硫醇化速率降低8.1个数量级。水合反应中预复合体已饱和（K_pre[H₂O] ≈ 670），观测速率常数独立于K_pre；而硫醇化中形式二级常数k_Cys = k_app × f_active正比于K_pre，修正全额传递。扩散相关设置（Collins-Kimball、Debye因子、对称因子）对速率和S_direct影响有限（分别不超过0.3和0.14个数量级）。用ΔE_el+ZPE替代ΔG则使水合速率增加约1.0个数量级。显式水分子数量改变1个可导致标准态与Benson-Okuno修正改变约4.4 kcal?mol^?1。与Fuoss接触离子对估计值K_F相比，硫醇盐预复合体超出7.4（非水合）和7.8（单水合）kcal?mol^?1。当K_pre上限设为10² M^?1时，S_direct分别从6.84和3.67降至2.88和0.30；精确匹配实验谷胱甘肽速率需硫醇盐预复合体去稳定2.1–2.5 kcal?mol^?1，即使如此仍高于Fuoss估计4.9–5.4 kcal?mol^?1，暗示静电相互作用不足以解释其稳定化。消除硫醇盐偏好则需硫醇盐轮廓较鸟嘌呤去稳定5.0（非水合）或9.3（单水合）kcal?mol^?1。pH 7.0–7.8和硫醇pK_a 8.5–9.2的变动对k_Cys影响不超过0.4个数量级，且S_direct不变。

2.3. 四个可靠性诊断

基于上述结果，研究人员提出任何计算速率常数或选择性应伴随的四个诊断标准。其一，预复合自由能的符号：正值ΔG^pre?isol表明预复合体形成热力学不利，K_pre低于1 M^?1不应解释为稳定相遇复合体的证据。其二，K_pre的大小：应与Fuoss非特异性接触离子对估计K_F比较，超出数kcal?mol^?1的稳定化可能提示隐式溶剂过度稳定化或需论证的特异性相互作用。其三，接近扩散极限：当k_equiv接近k_D时，Collins-Kimball修正确压缩速率常数及由此衍生的选择性描述符，任何进一步预复合体稳定化将被掩盖。其四，选择性的分解：应将选择性分解为预复合项（S_pre = log₁₀(K_pre,Cys/K_pre,Gua)）与活化项，若预复合项主导，结果应谨慎报告。

讨论与结论部分

研究人员强调，计算协议明确规定了每项修正，使其可独立测试。第一次水合的预复合体到过渡态势垒虽仅比Eyring推断值低2.0 kcal?mol^?1，但经对称因子校正后，计算观测速率常数仍为实验值的56倍。硫醇化的形式低浓度二级系数较实验高1.5–1.8个数量级，但这不构成预复合体饱和条件下对观测速率的浓度匹配预测。由于水合速率对相遇复合体自由能不敏感，该基准仅验证势垒，而未验证主导硫醇化速率及硫醇/鸟嘌呤比率的项。匹配谷胱甘肽二级速率系数要求硫醇盐预复合体去稳定2.1–2.5 kcal?mol^?1，且即使如此仍较Fuoss接触对估计高4.9–5.4 kcal?mol^?1，表明静电相互作用无法单独解释该稳定化；其可能源于真实特异性相互作用或隐式溶剂过度稳定化，但现有数据无法区分。硫醇/鸟嘌呤偏好的符号对此量级误差稳健，而其量级并非如此；因此，此类计算仅能自信报告偏好的方向。

据此，研究人员提出计算Pt(II)取代动力学的五步协议：（1）采用GFN2-xTB/ALPB（GOAT）进行构象采样；（2）在M06-2X/x2c-TZVPall与SMD溶剂化水平下优化稳定点；（3）应用标准态、溶剂浓度及Benson-Okuno修正；（4）经由过渡态理论结合Collins-Kimball扩散极限计算速率常数；（5）报告四项诊断以限定结果。若无这些修正，计算速率可能偏离实验达数个数量级。计算势垒与实验相符，若计算省略修正或采用调整的显式水分子数，不应被视为计算正确的证据。