咨询方案

从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

核心看点:

剑桥大学 Goodman 团队最新发表在 Chemical Science 上的工作提出 DP5q——用图卷积神经网络(GCNN)+ 分位数回归替代 DFT 计算,显著缩短了NMR结构确证的计算时间,实现数个数量级的加速。

这项工作与 ACD/Labs 的 CASE(计算机辅助结构解析)专家系统形成天然的互补关系:DP5q 擅长"快速确证候选结构",而 CASE 擅长"从零生成并验证所有可能的结构"。两者结合,构成了 AI 时代 NMR 结构解析的完整解决方案。

1|NMR结构确证的困境

核磁共振波谱(NMR)是小分子表征的核心技术。当一个新化合物被合成或分离后,如果没有参考谱图可供比对,结构确认就变得异常困难。这导致了许多后来被纠正的错误指认案例——据不完全统计,天然产物领域已报道超过1000例结构修正。

传统的人工 NMR 解析流程面临三大挑战:

从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

2|DP5q:用GCNN把DFT"踢出局"

2.1 从DP4到DP5再到DP5q

Goodman 团队在 NMR 辅助结构解析领域深耕多年:

  • DP4 方法(2010):基于 DFT 预测位移与实验位移的误差分布,从候选列表中选择最佳匹配结构。局限性在于必须假设正确结构已在列表中。

  • DP5 方法(2024):更进一步,构建了包含数万个原子环境表示及其 DFT 误差的数据集,对每个候选结构单独计算绝对概率置信度。但仍然依赖昂贵的 DFT 计算。

  • DP5q 方法(本文):使用图卷积神经网络(GCNN)+ 分位数回归直接预测化学位移的概率分布,完全消除对DFT计算的依赖。

2.2 核心创新:分位数回归的不确定性校准

与以往方法不同,DP5q 并非只预测化学位移的单一值,而是预测完整的概率分布(第1至99百分位共99个分位数):

损失函数 Ltotal = Lquantiles(分位数预测损失)+ Lorder(单调性正则化)

关键设计要点:

  • 采用修改后的 Huber 损失函数进行分位数回归

  • 通过单调性正则化确保预测分位数严格递增(符合累积分布函数性质)

  • 训练数据来自 NMR Shift DB 整理后的 22,248个唯一分子

  • 模型架构为 CASCADE 图卷积网络

为什么选分位数回归而不是正态分布?

以往方法(如 IMPRESSION)假设预测服从正态分布,需要显式拟合分布参数。分位数回归直接构建概率分布,且推理速度不随参考数据集规模增加而变慢,可任意扩展训练数据。

2.3 速度飞跃:数量级加速

这是最引人注目的结果——DP5q(基于GCNN)在保持准确度接近 DFT 级别的同时,实现了数个数量级的速度提升:

原文核心结论(Chemical Science, 2026):

  • "accelerating structural elucidation by several orders of magnitude by eliminating costly DFT calculations" —— 通过消除昂贵的DFT计算,将结构解析加速了几个数量级

  • "previously required CPU months, now can be done in seconds" —— 以前需要CPU运行数月的工作,现在可以在几秒钟内完成

  • DP5q 分析耗时约35分钟(含构象搜索),而传统 DP5 流程中仅 DFT 相关步骤就需要数周至数月不等

需要注意的是,上述加速幅度针对的是 DP5 分析流程本身(即从候选结构到概率评分这一步)。完整的结构解析流程还包含前置的构象搜索等步骤,这部分时间两种方法相同。

2.4 实测效果

研究团队在两类极具挑战性的测试上验证了 DP5q 的有效性:

结构修订案例(24对分子)

从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

图 8.  若干分子的实例展示,这些分子的首次公开报道结构均经过后续实验与计算方法重新鉴定。其中,化合物命名中的字母 a(例如 S1a)代表最初发表的推测结构,字母 b(例如 S1b)则代表经结构确证后修正的最终结构。

选取已发表后被修正的天然产物/合成分子案例(原始提案→后续通过重新合成或X-ray确认被推翻):

  • 使用单一构象(力场优化)时:22/24 案例的修正结构 DP5q 分数高于原始错误提案

  • 使用 ETKDGv3 多构象优化后:24/24 全部案例中正确结构的分数均高于原始提案

  • 唯一有挑战的是 S13(对称双萘结构 vs 非对称修正结构,NMR环境数目不同),多构象处理后也成功区分

相对立体化学解析(42个分子)

从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

对比基准:DP4-AI 仅 15/42,DFT-DP5 仅 16/42。DP5q (ETKDG) 不仅超越了所有基于DFT的方法,还避免了所有DFT计算成本。

推荐做法:日常使用推荐 ETKDGv3 + DP5q 分数阈值 0.2,可在可靠排除错误结构提案的同时实现极速分析。

3

DP5q很强大,但它解决不了"候选结构从哪来"的问题

细心的读者可能已经注意到一个关键细节——
DP5q /DP4/DP5的前提条件:用户必须先提供候选结构列表,DP5q 才能在这些候选之间做评分排序和概率判断。但问题来了:候选结构列表从哪来?

在 DP5q 的工作流程中,候选结构的来源有三种:

  1. 人工提出——研究者根据自己的化学直觉和经验,手动列出可能的异构体。这种方式高度依赖于个人经验,极易遗漏正确的结构(尤其是骨架新颖的天然产物)。

  2. 数据库检索——在 PubChem、SciFinder 等数据库中搜索已知化合物。但对全新结构无效。

  3. CASE 系统(计算机辅助结构解析)自动生成——这是唯一能做到无偏、系统地穷举所有可能结构的方式。

这正是 ACD/Labs 的 Structure Elucidator Suite (SES) 登场的时刻。

4

ACD/SES CASE专家系统:NMR结构解析的"全能选手"

从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

ACD/Structure Elucidator Suite (SES) 是 ACD/Labs 公司推出的计算机辅助结构解析(CASE)专家系统平台。它解决的问题正是 DP5q 的"前置环节"——如何从原始 NMR 数据出发,推导所有可能的结构候选,并进行系统性筛选。

4.1 CASE的核心工作流程

SE Core(SES的核心模块)的标准流程如下:

从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案
从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

关键优势——无偏向性(Unbiased):

CASE 系统不会像人类研究者那样"偏爱"某种骨架类型。它严格按照 NMR 实验数据的逻辑约束来生成和评估候选结构,这意味着即使正确结构与直觉相悖,也能被系统公平地考虑和排名。

4.2 SES平台的完整能力矩阵

从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

4.3 CASE的实际威力:几个经典案例

ACD/Labs 的 CASE 系统在大量实际应用中证明了其高效性及可靠性:

    5|DP5q 与 SES CASE:互补而非竞争

    理解了 DP5q 和 CASE 各自的能力边界后,二者的互补关系就变得非常清晰:

    从 DP5q 到 CASE:AI 加速 NMR 结构确证的完整解决方案

    5.1 CASE的"日常模式" vs "高难模式"

    需要特别说明的是,并非所有结构解析任务都需要引入 DP5q 或 DFT。根据 ACD/Labs 在超过 400 个已发表案例中的统计,CASE 的经验预测算法(HOSE码/神经网络/增量法)在日常工作中表现如下:

    日常模式——CASE 经验算法即可搞定(占绝大多数情况):

    • 骨架常见、结构相对简单的分子 → 数据库中已有大量相似片段参考

    • CASE 生成候选后,正确结构的平均偏差 d_A 通常 < 3 ppm

    • Top-1 排名即为正确结构的概率很高,无需额外验证

    • 典型场景:已知骨架的衍生物、常规合成中间体、大多数天然产物

    高难模式——需要 CASE + DP5q/DFT 协同(少数棘手案例):

    • 骨架极度新颖(如梯烷类 ladderane、新型杂环体系)→ 数据库中缺乏参考片段

    • 多个候选差异度非常接近,经验算法无法自信区分排名前几的结构

    • 缺氢/贫氢分子(H/C 重原子比 < 2),2D NMR 信息量有限

    • 立体化学争议,非对映体之间谱图差异极小

    简而言之:DFT 和 DP5q 不是 CASE 的"标配",而是面对疑难问题时的"火力支援"。

    5.2 理想的协同工作流

    CASE → DFT/GCNN 筛选 → 最终确认(三步走策略)

    Step 1:CASE 自动生成候选(SES SE Core)输入实验 2D NMR + MS 数据 → 系统穷举所有可能结构 → 用经验预测法快速初筛 → 输出 Top-N 候选(通常 N = 1-5)

    Step 2:高精度评分(DP5q 或 DFT)将 CASE 输出的 Top-N 候选送入 DP5q(秒级完成)或 DFT(小时级)进行精确化学位移预测和概率评分 → 排除假阳性

    Step 3:最终确认与报告综合 CASE 的差异度数据和 DP5q/DFT 的概率分数 → 得出最终结论 → 生成标准化报告

    这种协同模式已有成功先例。

    Buevich & Elyashberg (J. Nat. Prod. 2016) 首次展示了 CASE + DFT 协同策略:先用 CASE 从海量异构体中筛选出少量候选(3~15000个),再对这些有限的候选执行 DFT 高精度计算。这种"粗筛+精判"的组合比纯 DFT 方案效率高出数十到数千倍——而 DP5q 正好可以把 Step 2 中的 DFT 替换为 GCNN,使整个流程进一步加速。