分子内相互作用的局域能量分解：CovaLED方法及在生物分子组装体分子识别中的应用

《ACS Central Science》：Local Energy Decomposition of Intramolecular Interactions: The CovaLED Approach and Its Application to Molecular Recognition in Biomolecular Assemblies

　　

实现对复杂组装体中调控分子功能与生物活性的非共价相互作用的定量解析，是量子化学分析面临的重大挑战，尤其在核酸等标准力场表现不佳的体系中。局域能量分解（LED）可为分子间相互作用提供金标准耦合簇（coupled-cluster）描述，但当相互作用片段通过共价键连接时（如嵌入核酸骨架的配体或通过键合网络连接的官能团），其化学可解释性会显著下降——片段边界共享的电子密度会模糊能量贡献的物理意义。研究人员在此提出CovaLED，这是LED方案的扩展，可在LED框架内严格处理共价连接性。将该方案应用于核酸基识别体系验证了方法的能力：在核糖开关RNA-配体复合物中，CovaLED揭示了鸟嘌呤甲基化如何导致稳定氢键相互作用的损失，与实验结合亲和力趋势一致；在人类DNA片段中，该方法可精确量化主链内共价连接的核苷酸间相互作用。因此，CovaLED能够在共价结构与非共价识别本质交织的真实生物分子体系中实现耦合簇水平的能量分解。

《ACS Central Science》发表的这项研究针对生物分子识别能量解析的长期瓶颈，开发了共价适配的局域能量分解（CovaLED）方法，填补了高精度波函数理论与化学可解释性之间的空白。现有方法存在三重局限：经典力场对RNA等体系精度不足，传统量子化学能量分解无法处理共价连接片段的“离子化偏差”，人工智能结构预测缺乏物理机制解释。研究人员基于域基对自然轨道（DLPNO）框架扩展了LED方案，通过共享轨道均衡分配消除人工电荷效应，实现了完整生物分子中共价与非共价相互作用的耦合簇（CCSD(T)）水平分解。该方法在鸟嘌呤核糖开关中成功量化了O6-甲基化导致的氢键损失，解释了结合自由能下降三个数量级的实验现象；在DNA双链分析中验证了与简化模型的兼容性，为生物分子相互作用研究提供了无偏工具。

关键技术方法包括：1. 基于DLPNO-CCSD(T)的多层计算策略，对配体及结合口袋采用高精度基组（def2-TZVP(-f)），其余RNA区域采用HF/MINIX简化；2. CovaLED共价处理模块，通过指定共价键自动识别跨片段共享轨道，采用均裂分配系数（n_kx=n_ky=0.5）避免极化依赖；3. 片段对局域能量分解（fp-LED）映射技术，将电子制备能重分布为片段间相互作用项；4. 完整体系计算采用GFN2-xTB优化几何结构，ONIOM分层模型平衡精度与成本，所有计算通过ORCA软件包实现。

研究结果分为四个部分：

The CovaLED Approach：传统LED将共价共享轨道完全分配给单一片段，导致人为电荷分离。CovaLED通过轨道布居分析识别跨片段共享轨道|k?，将其分解为|k_x?和|k_y?组分，引入分配系数n_kx和n_ky（和为1）。默认采用均裂分配（n=0.5）避免极化敏感性，重构了Hartree-Fock能量（E^HF）和关联能（E^C）表达式，使共价连接的片段无需切断或添加连接原子即可获得无偏能量分解。

An Illustrative Case Study: The Ethane-Na⁺Case：乙烷-钠离子体系测试显示，标准LED因甲基片段带电产生符号相反的Na⁺···CH₃相互作用（-28.1和+13.3 kcal mol^–1），而CovaLED均等分配C-C键电子密度后，两个Na⁺···CH₃作用均为-13.3 kcal mol^–1，共价C-C键作用反映其与钠离子结合后的弱化程度（19.3 kcal mol^–1），符合化学直觉。

The Interaction of DNA Filaments in Human DNA：人类DNA双链（序列5′-CTGAGGA-3′/3′-GACTCCT-5′）分析表明，去除主链的简化模型（14个碱基片段）用标准LED得到总相互作用能-175.3 kcal mol^–1，氢键与堆叠作用模式与已知一致；保留主链的完整模型（16个片段）经CovaLED处理后总能为-181.0 kcal mol^–1，碱基间相互作用（-206.1 kcal mol^–1）与简化模型（-209.5 kcal mol^–1）几乎一致，差异仅来自新增的主链相互作用，证明方法无伪影引入。

RNA-Ligand Interaction: Guanine Riboswitches：鸟嘌呤核糖开关复合物计算显示，天然配体鸟嘌呤（6UBU）结合能为-114.1 kcal mol^–1，O6-甲基鸟嘌呤（3F06）为-103.8 kcal mol^–1，与实验K_d值（0.004 μM vs 23 μM）趋势吻合。fp-CovaLED图谱揭示：6UBU中C74残基与配体羰基（片段2）作用达-18.5 kcal mol^–1，U51形成两个氢键（-15.5 kcal mol^–1）；甲基化后C74氢键从3个减至2个，与片段2作用降至-0.4 kcal mol^–1，氢键网络破坏是结合减弱的主因。

讨论与结论部分指出，CovaLED消除了传统片段划分的“离子化偏差”，首次实现了完整生物分子体系的耦合簇水平能量分解。DNA体系验证表明，即使不切断主链共价键，也能获得与连接原子模型一致的碱基相互作用模式。核糖开关应用证明该方法能将微小的化学修饰（O6-甲基化）转化为可量化的物理相互作用变化，明确氢键损失是结合能下降的核心因素。这一方法架起了高精度电子结构分析与化学可解释性之间的桥梁，可推广至酶催化、材料化学等共价与非共价作用交织的领域，为生物分子识别机制研究提供了通用工具。