锌(II)配位模型的DFT初步筛选：胍、Galegine、二甲双胍与L-精氨酸的质子化、第一配位层水合及电子描述符

《Biophysica》：A Preliminary DFT Screening of Zn(II) Coordination Models for Guanidine, Galegine, Metformin, and L-Arginine: Protonation, First-Shell Hydration, and Electronic Descriptors

　　

锌稳态失调与神经退行性过程相关，但胍基和双胍基分子识别Zn(II)的分子层面模型必须考虑质子化、构象柔性以及与第一水合层的竞争。研究人员采用层次化密度泛函理论（DFT）策略研究胍、Galegine、二甲双胍和L-精氨酸。中性配体、相关质子化形式及简化[Zn(L)]²⁺加合物采用Becke三参数Lee–Yang–Parr（B3LYP）泛函、6-311+G(d,p)基组和积分方程形式极化连续介质模型（IEFPCM）水溶剂化（Set I）计算；[Zn(H₂O)₆]²⁺及显式水合的[Zn(L)(H₂O)₅]²⁺第一配位层取代模型采用B3LYP结合Grimme D3色散与Becke–Johnson阻尼[D3(BJ)]/def2-SVP/IEFPCM(水)（Set II）评价。谐波频率分析支持所报道的驻点结构在相应计算模型内为局域极小值。在Set I中，简化[Zn(L)]²⁺加合物的Zn-N距离为1.924–1.935 ?。在独立的Set II水合配位模型中，直接Zn-N配位得以保留，Zn-N距离为2.007–2.028 ?。质子化使胍基/双胍基片段周围的分子静电势（MEP）向正值移动，降低静电给体可达性。定向扭转检查进一步表明，柔性配体描述符依赖构象；对于中性L-精氨酸，最终分析采用的A6羧酸-OH旋转异构体的前沿轨道分离能为6.214 eV，而初始pre-A6构象异构体为5.984 eV。这些结果确定了一个初步的基于描述符的电子结构和第一配位层配位框架，而非对结合亲和力、主要生理形态或药理活性的预测。

神经退行性疾病如阿尔茨海默病（AD）和帕金森病（PD）为多因素疾病，遗传、代谢、炎症和环境过程交织。代谢功能障碍受到关注，AD中胰岛素信号受损、线粒体功能障碍、氧化应激和神经炎症与2型糖尿病和认知下降的联系相关；PD也有代谢改变报道，但因果关系未完全解决。Zn(II)是氧化还原惰性必需金属离子，参与胰岛素储存和信号、突触功能、酶催化、蛋白稳定和抗氧化防御；锌稳态紊乱与氧化应激、信号改变、线粒体功能障碍和异常蛋白聚集相关，并可调节淀粉样β、α-突触核蛋白等蛋白行为，但其生物效应依赖浓度、细胞区室、配体环境和内源给体竞争。胍、Galegine、二甲双胍和L-精氨酸提供结构相关系列，用于考察胍基或双胍基对质子化和金属配位的响应；胍-金属相互作用具有立体化学偏好，二甲双胍及相关双胍有实验记录的金属结合行为，L-精氨酸是生物上不同的含胍基氨基酸。问题在于：水合Zn(II)常见配位数为4–6，连续介质溶剂化不能显式描述Zn-O(水)接触或第一配位层配体取代；质子化会增强正静电特征并降低定域氮孤对电子对阳离子金属中心的可达性；柔性配体描述符还依赖构象。因此，研究人员用层次化DFT分离质子化效应、内在Zn-N配位、显式第一配位层水合和构象敏感性，并将几何、MEP/ESP和HOMO–LUMO分离能作为比较描述符，而非结合强度、热力学稳定性、主要生理形态或药理活性度量。该研究的结论是：简化[Zn(L)]²⁺保留短Zn-N接触，显式水合模型中直接Zn-N配位仍存在，但Zn-O接触更不均一；质子化使MEP正移；构象影响L-精氨酸前沿轨道分离能。其意义在于提供初步电子结构和第一配位层配位框架，用于生成后续合成和物理化学表征可检验的配位假设，而不是预测生理物种或药物活性。

关键技术方法概括（不超过250字）：研究人员使用Gaussian 16 Revision C.01开展DFT计算。Set I在B3LYP/6-311+G(d,p)/IEFPCM(水)水平优化中性配体、相关质子化形式及简化[Zn(L)]²⁺加合物，并进行谐波频率验证；MEP/ESP映射于ρ=0.0004 a.u.电子密度等值面，HOMO/LUMO以0.02 a.u.等值面可视化，前沿轨道分离能按ΔE_HL=ε_LUMO?ε_HOMO计算。柔性配体在B3LYP/6-31+G(d,p)/IEFPCM(水)下以60°增量进行定向弛豫扭转扫描，A6再于最终Set I水平优化和频率验证。Set II在B3LYP-D3(BJ)/def2-SVP/IEFPCM(水)水平处理[Zn(H₂O)₆]²⁺和含五个显式水分子的[Zn(L)(H₂O)₅]²⁺模型，并进行几何优化、谐频验证、结构分析和前沿轨道分析。原文未涉及样本队列来源。

研究结果如下。3.1 分子模型的几何与结构响应：结构分析从孤立配体模型到简化Zn(II)加合物再到显式水合第一配位层模型。胍保持最刚性、近平面CN₃骨架，Galegine、二甲双胍和L-精氨酸外围构象自由度逐步增加；质子化不消除胍基/双胍基核心的共振稳定特征，但改变N-C-N键极化和外围N-H或取代基矢量取向。Galegine质子化位点经同级比较，亚氨基质子化异构体E_elec=?401.3129928 Eh、G₂₉₈=?401.1458141 Eh，氨基质子化替代结构E_elec=?401.2676897 Eh、G₂₉₈=?401.0986657 Eh，前者G₂₉₈低29.59 kcal mol^?1，且C11-N距离近等，故采用亚氨基质子化形式为参考。中性L-精氨酸A6经无约束最终水平优化和频率验证仍为独立局域极小值，G₂₉₈比初始中性构象低3.64 kcal mol^?1，被采用为工作结构。简化[Zn(L)]²⁺中四个中性配体均保留末端胍基/双胍基氮与Zn(II)配位，Zn-N距离分别为1.925、1.926、1.935和1.924 ?。显式水合模型中[Zn(H₂O)₆]²⁺的Zn-O距离为2.088–2.098 ?；配体取代后Zn-N距离为2.018、2.007、2.008和2.028 ?，Zn-O距离分别变为2.117–2.188、2.071–2.245、2.125–2.178和2.122–2.166 ?；Galegine的Zn-O范围最宽。由于Set I和Set II理论水平不同，这些距离仅在各自集内描述性解释。3.2 模型层级中的分子静电势：中性配体在N/O给体域存在局部负电势；质子化使全局尤其胍基/双胍基片段显著正移，降低质子化氮给体对阳离子Zn(II)的静电可达性；生理L-精氨酸模型仍静电异质，其含氧区域相对富电子而质子化氮位点偏正。中性L-精氨酸MEP采用A6参考构象，羧酸-OH取向相对初始构象改变，说明含氧末端局部静电特征可依赖构象，MEP被解释为A6结构描述符而非构象平均量。简化[Zn(L)]²⁺显示强正Zn中心表面及配体相关残余各向异性，显式水合模型则将静电场分布于配体和水第一配位层；简化加合物ESP数值限约±0.579、±0.576、±0.587和±0.580 a.u.，水合模型约±0.341、±0.318、±0.355和±0.324 a.u.，因协议不同不作跨方法定量比较。3.3 前沿分子轨道与HOMO-LUMO能隙：Set I中性配体ΔE_HL为6.589 eV（胍）、6.190 eV（Galegine）、6.220 eV（二甲双胍）和6.214 eV（L-精氨酸A6）；相应质子化/生理模型为8.248、6.493、6.161和6.630 eV。中性L-精氨酸pre-A6构象ΔE_HL=5.984 eV，采用A6后变为6.214 eV，变化约0.230 eV，说明前沿轨道分离能可依赖构象。Set I简化加合物ΔE_HL约5.277 eV（胍）、4.174 eV（Galegine）、4.003 eV（二甲双胍）和4.478 eV（L-精氨酸）。Set II中[Zn(H₂O)₆]²⁺的ΔE_HL约9.017 eV，第一配位层配体取代模型为7.016、6.617、5.939和6.361 eV，相对六水合参考的缩小主要来自占据前沿能级变化；这些值仅为优化模型描述符，不转化为配体排序或取代有利性代理。3.4 整合解释与模型局限：结果分离了三类效应：质子化控制给体可达性并强烈改变孤立配体静电特征；中性配体参考模型可考察胍基/双胍基给体对直接Zn-N配位的局部立体电子响应；显式第一配位层水合检验五个水O给体竞争下Zn-N接触是否仍为驻点结构基序。水合极小值中Zn-N持续存在是结构结果，不证明中性配体在生理pH下占优势或从[Zn(H₂O)₆]²⁺发生配体取代在热力学上有利。Set I提供孤立配体状态和简化加合物内部一致比较，但Pople型基组和无显式第一层水限制其描述溶剂化Zn(II)；Set II使用Zn兼容def2基组、D3(BJ)色散和显式第一层水，但仍为可处理簇-连续介质筛选模型，无抗衡离子、竞争内源配体、动态溶剂交换或pH依赖形态，def2-SVP为实用筛选水平而非基准质量基组。定向扭转敏感性分析表明柔性配体必须显式处理构象选择；一维扫描不等同穷尽多维构象集合，A6基L-精氨酸电子描述符应视为构象特异而非构象平均性质或唯一全局极小值证据。配体取代热力学测试需同级反应、配体质子化、构象态、标准态和微溶剂化一致处理，并需稳定常数、光谱、生物给体竞争和pH依赖形态连接溶液化学。

讨论部分总结：该层次化DFT研究区分质子化效应、内在Zn-N配位、显式第一配位层水合和构象敏感性，但轨道能隙、MEP分布和局域极小几何应仅用于形成机理假设，不能用于排序结合亲和力或宣称唯一溶液构象；未来实验应先经合成和结构表征确认候选配位基序，再评价水稳定性、pH依赖形态、与水和内源给体竞争及稳定常数，之后才可合理开展生理相关条件测试或配体设计应用。结论部分翻译/概括：该层次化DFT研究区分胍、Galegine、二甲双胍和L-精氨酸模型中的质子化效应、内在Zn-N配位、显式第一配位层水合和构象敏感性。定向扭转筛选确认Galegine和二甲双胍描述符结构靠近所探索坐标的低能区，完成生理质子化状态L-精氨酸模型的五个计划扫描，并识别出显著更低的中性L-精氨酸羧酸-OH旋转异构体A6；A6被采用为最终结构、MEP/ESP和前沿轨道分析的中性L-精氨酸参考构象，其ΔE_HL相对pre-A6由5.984变为6.214 eV，表明其前沿轨道描述符依赖构象。Galegine同级质子化位点比较中，亚氨基质子化形式G₂₉₈比氨基质子化替代结构低29.59 kcal mol^?1，故采用亚氨基质子化形式为质子化Galegine参考，氨基质子化异构体仅留作比较。质子化使胍基/双胍基MEP表面正移，降低质子化氮给体对阳离子Zn(II)的可达性；中性配体计算因而作为内在电子/配位参考，而非主要生理物种。所有简化[Zn(L)]²⁺模型在Set I内保留约1.92–1.94 ?的短Zn-N接触；独立Set II水合模型中直接Zn-N配位仍存在，Zn-N距离为2.007–2.028 ?，Zn-O接触较[Zn(H₂O)₆]²⁺更不均一；因Set I和Set II模型化学不同，这些值仅在各自集内描述性解释，而非同级水合诱导伸长度量。前沿轨道分析显示质子化、配位和水合引起进一步重组；Set II中配体取代相对六水合参考缩小HOMO-LUMO分离，主要通过占据前沿能级变化。这些结果仅支持初步电子结构和第一配位层配位框架，不建立结合亲和力、有利配体交换热力学、生理Zn(II)形态或药理活性；其作用是生成后续合成和物理化学表征可检验的配位假设，未来需结合同级取代自由能、构象和pH依赖形态、竞争内源配体、稳定常数测量、光谱和生物验证后再赋予更广泛生物或配体设计意义。