第8章 · 8.11

8.11 从谱图反推结构

Inferring structures from spectra
本页范围与学习目标
比较谱检索、约束枚举、候选重排和学习生成,用C₄H₁₀O逐步筛选展示分子式与实验约束的作用,并分别评价生成覆盖与排序。
先读图解和核心解释;需要细节时再展开。
先看图 · 阅读入口

谱到结构的反问题先规定答案范围

模型提出的结构是需要核对的候选,最终报告要解释完整证据。

谱到结构的反问题先规定答案范围 步骤 1 准备观测约束 分子式、位移与关系 保留不完整状态 步骤 2 检索或生成 建立具体候选分子 记录范围与输入先验 步骤 3 排序并检验 反向预测全部主要证据 设计候选区分实验 谱到结构的反问题先规定答案范围 步骤 1 准备观测约束 分子式、位移与关系 保留不完整状态 步骤 2 检索或生成 建立具体候选分子 记录范围与输入先验 步骤 3 排序并检验 反向预测全部主要证据 设计候选区分实验
图 8.11-1 本页概念关系图。箭头表示本页采用的检查顺序。

读图时保留这项条件:排名高低只能比较已提出的候选,连接与立体正确性也需分开评价。

8.11.1 反问题先规定允许的答案

正向预测给定结构而计算谱;反向解析给定谱而寻找结构。后者常有多个解:有限分辨率、峰重叠、对称与条件平均使不同结构可能具有接近的观测。因此,程序输出一个结构字符串只完成了候选提出,能否作为指认结论还要看它解释了哪些数据、排除了哪些替代方案。

展开说明与推导

实际任务可以分为库内检索、给定候选排序、约束下从头生成,以及学习模型的结构生成。库内检索受数据库覆盖限制;排序受输入候选集合限制;约束枚举受价态与约束设定限制;学习生成受训练范围、解码与采样限制。四类任务都可提供候选,但它们对“还有没有其他可能”的回答不同。

先说明目标是连接结构、相对构型还是绝对构型。通常非手性溶液中的对映体不能仅由普通位移和 J 区分。模型若输出带绝对立体标记的字符串,还要检查这些标记来自独立手性信息、输入先验或训练偏好。连接正确而立体标记未经支持,应分别报告。

输入也必须保留证据层次。原始 FID、处理谱、峰表、积分和人工指认并不相同。由峰表输入的模型通常没有收到完整线形;由谱图像输入的模型仍受谱宽、像素、参考和显示归一化影响。名称中写“原始谱”时,还需分清频域谱信号与仪器采集的时域 FID。

定义 · 搜索空间

搜索空间是允许提出的结构集合。要评价候选是否充分,必须同时说明这个集合怎样建立,以及程序实际搜索了其中多大范围。

8.11.2 分子式缩小什么范围

可靠分子式规定元素与原子数量,并通过不饱和度限制环和多重键的总量。它没有给出连接次序,也没有把所有同分异构体压成一个答案。分子式若来自质谱,需要先确认离子类型、加合物、荷电状态与同位素信息,避免将测到的离子组成直接当成中性待测物组成。

展开说明与推导

以 C₄H₁₀O 为教学例,不饱和度为零。假定单一中性闭壳层、通常碳四价氧二价、全部单键且分子连通,就只需考察无环连接图。四个碳与一个氧共有十种可能原子对,从中选四条边,先检查连通和价态,再按碳编号置换去重。

本书用独立脚本枚举得到一百一十二个符合条件的带编号图,去除编号差异后是七种连接:四种醇和三种醚。四种醇为丁醇的两个直链位置异构体、异丁醇与叔丁醇;三种醚为乙氧基乙烷、甲氧基丙烷与甲氧基异丙烷。枚举没有展开立体异构体,所以七是连接数。

现在给出经独立编辑实验与指认确认的碳类型:两个甲基、一个亚甲基、一个次甲基,没有无氢碳。这保留丁-2-醇与异丁醇两种连接。叔丁醇因有无氢碳被排除,丁-1-醇因有三个亚甲基被排除,三种醚的碳类型也均不符合。即使暂时没有观察到 OH,类型约束仍可完成这一步。

再用 HSQC 将亚甲基氢对应到其碳,并观察到该亚甲基与一个甲基的 COSY 相关。这里还通过分裂、选择性实验与完整相关链确认它属于相邻碳之间的主要三键耦合。丁-2-醇具有 CH(OH)—CH₂—CH₃ 局部链,异丁醇的 CH₂OH 直接连接次甲基,没有相邻甲基。由此支持丁-2-醇;若只有一个未经核实的弱交叉峰,则可能涉及长程耦合,不能直接完成排除。

图 8.11-2 分子式与实验约束逐步筛选C₄H₁₀O · 限定价态 · 7 种连接2 CH₃ + 1 CH₂ + 1 CH类型计数,非仅数可见峰丁-2-醇异丁醇CH₂ ↔ CH₃:保留无该相邻链:排除
图 8.11-2 每次筛选依赖明确条件。最后得到的是连接结构;丁-2-醇的绝对构型尚未确定。
展开说明与推导

这套计数有前提。若样品是混合物、类型实验灵敏度不足或相关峰是伪影,不能照搬筛选结论。原子数与环境数也不能混淆:对称等价可让两个甲基只出现一个碳环境,重叠则会让不等价环境暂时无法分开。没有可靠类型信息时,应保留更宽的候选集合。

8.11.3 检索与特征预测的区别

直接谱检索将待测谱与库中参考谱比较,再返回近邻。分子式可以先过滤库,溶剂、温度和参考条件可以帮助解释差异。若真实化合物不在库中,最近邻只能提示类似环境或骨架,不能因为总有一个第一名就把它认作待测物。应检查库外样本的处理方式。

展开说明与推导

SMART 2.0 的研究以 HSQC 数据支持天然产物候选检索与快速注释,用学习表示寻找相关结构。[1] 这里的 HSQC 输入反映直接相连氢碳的局部环境;未带氢碳、完整长程连接和立体关系还需其他证据。该论文中的应用可以帮助缩小探索范围,但不等于任何一张 HSQC 都足以唯一恢复完整结构。

DeepSAT 从 HSQC 学习化学指纹、分子量与结构类别,再利用这些特征搜索结构库。[2] 结构库条目不必逐一具有实测参考谱,这缓解了谱库覆盖限制;输出仍受结构库覆盖与特征预测误差影响。指纹相近表示某些局部特征相近,不保证连接完全相同。

例如甲氧基片段在多个候选中均存在,一个特征模型识别到它后可以排除明显不含此片段的方案,却不能单凭这一特征决定甲氧基接在哪个位置。应把预测子结构作为待验证约束,检查是否能与分子式、所有核数和相关链同时共存,避免将分别高分的片段拼成价态不合法的结构。

约束驱动的 Sherlock 则将结构生成与去重复发现等功能接入工作流,结合明确的实验约束、指认及数据库信息。[3] 使用时要区分已知数据检索和约束生成分别贡献了什么。它的发表能力与所用核种有具体范围,不能泛化为只输入任意氢谱即可自动完成所有结构问题。

8.11.4 从头生成与多模态学习

从头生成可以逐条建立连接图,也可以生成表示结构的 SMILES 序列。约束枚举逐步剪除不满足条件的连接;学习生成则根据训练得到的分布提出较可能的序列。束搜索或随机采样保留有限候选,因此其返回集合通常没有覆盖所有化学上可行连接的证明。

展开说明与推导

Hu 等提出从常规一维氢、碳谱到结构的多任务学习方法,并训练子结构相关任务帮助生成。[4] 研究中的谱到结构训练与测试主要使用预测谱,并在特定分子大小与元素范围内评价。本节只介绍已核实的任务设置,不把其模拟测试结果作为未知实验样品的成功率。

该工作的输入不依赖预先给定分子式,说明模型可以同时提出组成与连接。但在实际实验中,可靠的分子式仍是有价值的独立约束。若生成结构与已确认分子式不同,应显示冲突并检查输入或模型;不能因为方法允许无分子式运行,就忽略已经取得的组成证据。

NMRViT 将一维氢、碳谱和二维 HSQC 的频域谱信号分块编码,生成候选结构,并研究分子式提示与基于位移的后续重排。[5] 该研究分别考察模拟数据训练与实验数据适配,体现二者之间存在分布差异。它处理的频域输入也不能被误称为已经直接读取所有厂商 FID。

多模态方法利用各实验提供的互补信息:碳谱约束碳环境,HSQC 建立直接氢碳配对,氢谱包含积分与耦合线形。如果再结合准确质量、IR 或其他数据,应记录每种输入的处理与来源。模态增加不保证错误减少;某一谱属于不同组分或参考错位时,融合会引入互相矛盾的信息。

方法 · 生成后的检查

依次检查结构字符串可解析、价态与组成合法、候选去重、核数符合、相关约束可满足,以及模拟谱与实验残差。前两项通过只说明化学表示有效。

8.11.5 排名怎样变成证据

候选排序可依据谱相似度、预测残差、约束满足程度或学习评分。不同分数并不天然可相加:一个是余弦相似度,另一个是 ppm 残差,再一个是生成对数权重时,需要规定尺度与组合依据。正式报告应保存原始分项与权重,不能只保留一个不明来历的总分。

展开说明与推导

硬约束与软约束也要区分。可靠组成冲突通常要求排除或回查组成;预测位移略有偏差常可容许。HMBC 强度与路径不是严格一一对应,漏掉一个峰未必就是连接不存在。把不确定观测全部设为硬约束,可能让真实结构在排序前已被删掉。

评价模型要同时报告生成覆盖和排序能力。教学验证集有二十个分子,正确结构只进入十二个候选列表,其中八个排第一。全体首位正确率为八除二十,即 40%;候选召回率为十二除二十,即 60%;在已召回子集上的首位比例为八除十二,约 66.7%。只报告第三个数会掩盖另外八个分子根本没有生成正确候选。

图 8.11-3 生成与排序的分别评价20 个独立教学样本8 首位正确48 未召回4 个正确结构位于候选列表的其他位置候选召回 60% · 全体首位正确 40%
图 8.11-3 筛选和重排只能在已提出的候选上工作。图中数字是教学构造,不是任何研究模型的实测成绩。
展开说明与推导

不同字符串可能表示同一连接;互变异构、盐、立体标记和芳香表示也会影响“完全正确”的定义。评价前应规定标准化与去重方式,并区分精确结构命中、骨架命中和子结构相似。提高骨架相似度很有帮助,但不能将这种指标改写成完整结构的确定率。

新样品往往缺少真值,无法当场算正确率。此时要检查候选间关键差异,设计能排除替代方案的实验,并查看预测的不确定度与残差。若所有候选都与高可信数据冲突,应扩展搜索或回查数据,不能自动选择冲突最小者并写成“已确定”。

误判 · 首位候选

第一名表示在给定候选和评分规则下排名最高。候选覆盖不足、共同预测误差或错误输入,都可能使第一名仍然错误。

8.11.6 形成可复核的结论

完成结构分析时,将分子式依据、候选来源、搜索范围、分项评分、峰映射和排除理由一起保存。对每个关键连接列出支持观测,再说明未定的立体关系。本节的七种连接枚举已实际运行;所列学习模型与 Sherlock 仅核实论文及公开资料,没有在本书中运行。

展开说明与推导

一个可复核的结论应允许他人看到哪些观测改变了候选集合。对教学丁醇例,碳类型先将七种连接减到两种,可靠的亚甲基—甲基相关再选择连接;绝对构型仍保留未知。后续语言模型或工具代理可以帮助组织这张证据表,但不能用自然语言的流畅程度替代其中的实验判断。

8.11.7 习题与解答

习题8.11-1 分子式确定后还能有多少连接?

在本节规定的 C₄H₁₀O 范围内,已确认没有 CH 或无氢碳,是否能够立即认作乙醚?

解答不能。丁-1-醇有一个甲基与三个亚甲基;乙醚和甲氧基丙烷均有两个甲基与两个亚甲基,也无 CH 或无氢碳。需进一步可靠区分类型计数、氧相连环境、对称与相关链。分子式与这两项缺失条件不足以唯一指认。

习题8.11-2 评价的分母

三十个测试分子中,十八个被正确召回,九个首位正确。分别计算候选召回率、全体首位正确率与召回子集首位比例。

解答依次为 18/30=60%、9/30=30%、9/18=50%。报告必须注明分母;只写“首位正确 50%”而省略子集定义,会高估全体任务表现。

习题8.11-3 两个模型都给出同一结构

谱检索与生成模型给出同一第一候选,但该结构的碳数与可靠分子式不符,应如何处理?

解答先标出组成冲突并回查离子解释、结构标准化和输入;可靠分子式若仍成立,应排除该候选。两个模型可能共享数据或偏差,意见一致不能消除组成矛盾。需扩大候选来源或修正输入,再用完整实验约束检查。

8.11.8 术语

反问题
由观测数据推求产生观测的结构或参数。
库内检索
在给定数据库中寻找匹配或近邻。
从头生成
依据约束或学习分布构造候选结构。
候选重排
用额外评分调整已提出候选的次序。
化学指纹
编码结构局部特征的表示。
多模态输入
结合不同谱种或其他测量信息的输入。
束搜索
保留有限高评分序列分支的解码过程。
候选召回率
真实结构进入所返回集合的样本比例。
结构标准化
按规定统一组成、表示与比较形式。

8.11.9 参考文献

  1. Reher R, Kim H W, Zhang C, et al. A Convolutional Neural Network-Based Approach for the Rapid Annotation of Molecularly Diverse Natural Products. J Am Chem Soc, 2020, 142:4114–4120. SMART 2.0 原始研究。
  2. Kim H W, Zhang C, Reher R, et al. DeepSAT: Learning Molecular Structures from Nuclear Magnetic Resonance Data. J Cheminform, 2023, 15:71. 原始研究。
  3. Wenk M, Nuzillard J M, Steinbeck C. Sherlock: A Free and Open-Source System for Computer-Assisted Structure Elucidation. Molecules, 2023, 28:1448. 原始研究。
  4. Hu F, Chen M S, Rotskoff G M, Kanan M W, Markland T E. Accurate and Efficient Structure Elucidation from Routine One-Dimensional NMR Spectra Using Multitask Machine Learning. ACS Cent Sci, 2024, 10:2162–2170. 原始研究。
  5. Han C, Pan X, Zhang Y. End-to-end molecular structure elucidation from multimodal NMR spectra images using vision transformers. Chem Sci, 2026. NMRViT;作者代码。