天然产物 Zerumol A 的计算机辅助结构修正

天然产物的分离与结构鉴定是全球天然产物化学研究的核心内容,其驱动力主要源于对新结构生物活性分子的持续探索1。在新型有机物结构解析中,研究人员通常综合运用一维/二维核磁共振(NMR)、高分辨率质谱(HRMS)、红外光谱(IR)及紫外-可见光谱(UV–vis)等技术。然而,基于波谱数据的结构解析至今仍是一项复杂的逻辑组合问题,需处理大量谱图与结构之间的关联信息,导致实践者并不总能成功解析纯化单体的结构,有时甚至造成错误化学结构的发表。研究表明,约有10%的已报道天然产物结构存在误判2,这对后续研究带来严重困扰——例如基于错误结构开展生物学研究时,目标分子可能实为已知化合物或具有完全不同的骨架。
从方法论角度看,基于波谱数据的结构解析属于"逆问题"3,即从观测数据反推结构模型,理论上可能存在多个解。通过引入额外约束条件(如空间构象、化学环境等),可排除与实验证据矛盾的候选结构,从而逼近唯一可靠解。该思路已实现在计算机辅助结构解析(Computer Assisted Structure Elucidation,CASE)系统中,其能够生成所有符合现有谱图信息的化学结构,并筛选出最合理的候选分子3。进一步借助密度泛函理论(DFT)方法4,5对候选结构进行化学位移预测与验证,可提供关键的理论佐证。大量研究案例表明,CASE与DFT相结合的策略是发现和修正结构误判的有效工具。需要强调的是,在此过程中必须充分融入化学直觉与结构规则6–10,才能提升解析结果的可靠性。
近期,Zhang等11报道了从草果(Alpinia zerumbet)果实中分离得到的新型天然产物zerumol A (1)的结构(图1)。

图1. 化合物1的结构通过¹H NMR、¹³C NMR、HMBC、COSY和ROESY谱图数据,以及质谱分析(ESI-HRMS测得分子式为C20H22O5)共同解析确定。相对构型依据ROESY谱中H-3与H-5的相关信号确定。C-1与C-2间双键的构型根据H-1较大的耦合常数判定为Z式。绝对构型通过实验与计算的ECD谱图对比确定。
作者用于解析化合物1结构的光谱数据列于表1。
表1. Zerumol A结构解析所用光谱数据11。表中同时列出基于HOSE 算法及DFT DU8
ML方法12,13计算的¹³C NMR化学位移预测结果。

值得注意的是,含有氧杂环丁烷(oxetane)环的结构在天然产物中较为罕见14,且在二芳基庚烷类天然产物中尚未见报道15。根据《天然产物词典》检索结果,此类结构在已知天然化合物中的占比不足0.5%16。此外,氧杂环丁烷中氧-碳键的¹³C NMR化学位移值与环状及无环醚类化合物的范围高度重叠,这显著增加了氧杂环丁烷片段的识别难度16。
基于上述分析,我们首先采用ACD/Structure Elucidator专家系统17,18对结构1的正确性进行验证。该系统通过三种经验方法预测了1的¹³C NMR化学位移(图2)。

图2. Zhang等11提出的zerumol A (1)结构,其¹³C化学位移预测分别采用基于HOSE算法、神经网络及增量法三种方法完成。图中标注的化学位移值为实验测定结果。三种计算方法得到的¹³C NMR化学位移平均偏差依次表示为dA、dN和dI。各原子按实验值与计算值的差异进行颜色标记:绿色表示偏差在0–3 ppm范围内,黄色表示偏差大于3 ppm但不超过15 ppm。红色箭头指示非标准HMBC相关信号(NSC),即相关路径跨越三个以上化学键(nJCH,n > 3)。
经验预测值与实验值的平均偏差(1.0–2.5 ppm)处于正确结构可接受的范围,这一结果通常可视为合理。然而,HMBC谱中由红色箭头标示的两处相关信号跨越了五个化学键,而此类长程相关在HMBC谱中通常信号极弱甚至无法观测。值得注意的是,在HMBC谱图(图S1)中,对应的两组交叉峰(6.63/77.8 ppm和6.52/77.8 ppm,见表1)强度相当,且与常规2–3JCH相关信号的强度相近,这一矛盾表明结构1很可能存在误判。
进一步采用DFT DU8ML方法12,13进行¹³C NMR化学位移预测(表1),所得均方根偏差(rmsd)高达2.92 ppm,明确指示结构1不正确。
为深入验证该结构或寻找可能的替代结构,我们再次使用ACD/Structure Elucidator进行分析。通过对图2的直观分析发现,三个绿色标记的片段在结构1中清晰可辨(图3)。在这些片段中,¹³C NMR化学位移的预测值与实验值偏差均在3 ppm以内。需特别指出的是,偏差最大的碳原子均位于氧杂环丁烷环内(另见表S1)。

图3. Zerumol A的候选结构示意图,其中三个独立结构片段以红色圆圈标出,表示这些片段的¹³C NMR化学位移预测值与实验值较为接近。
我们进一步探究了如何基于所选片段及假设构成氧杂环丁烷环的原子生成可能结构。将绿色片段分别输入Structure Elucidator后,结合分子式构建了分子连接图(MCD)(图4)。为全面生成所有可能的结构式,我们暂时解除了HMBC和COSY相关信号所施加的约束条件。

图4. Zerumol A候选结构的分子连接图(MCD)。碳原子杂化类型通过颜色标记:sp2为紫色,sp3为蓝色。标签"ob"(必需)和"fb"(禁止)由程序自动标注于与杂原子相邻的碳原子上。HMBC及COSY相关约束已关闭。
最终,系统在5秒内共生成14个候选结构。通过对这些结构进行¹³C NMR化学位移预测,并按dA偏差值升序排列输出结果(图S2),发现所有生成的可能结构均含有二氢吡喃、二氢呋喃、氧杂环丁烷或环氧环单元。排名前三的结构(2–4)如图5所示。

图5. 输出文件中排名前三的结构。需注意,在结构4中,HMBC相关信号13→3和9→3的路径长度为五键(红色箭头所示),而在结构2和3中,这些相关信号为三键(绿色箭头)。其他细节参见图2说明。
分析表明,结构2为最可能的正确结构,其在三种预测方法中的平均偏差均为最小值。结构3的偏差范围同样可接受,而文献中提出的zerumol A结构11(即结构4)仅位列第三。值得注意的是,结构4的平均偏差和最大偏差均显著高于结构2和3。
在结构2和3中,先前提及的5JCHHMBC相关路径为标准的三键长度,与其交叉峰强度相符(图S1)。进一步通过COSY谱验证结构2:谱图中不仅观察到作者已标注的相关信号,也完全符合结构2的预期相关模式(图6,另见图S3)。此外,将zerumol A (1)的¹H NMR谱与已知天然产物519(其结构与2高度相似,图6)进行对比,发现两者脂肪区信号几乎一致(图S4),进一步支持结构2的合理性。

图6. Zerumol A修正后的平面结构(结构2)及对应的COSY谱交叉峰。作者11提及的相关信号以蓝色箭头(标准相关)和红色箭头(非标准相关)标示。黄色虚线箭头表示修正后结构预期存在的相关信号。(关于颜色标识的说明,请参阅本文的网络版本。)
为进一步验证结构2的正确性,我们采用DFT DU8ML方法12,13对候选结构2和3(图5)进行了¹³C NMR化学位移预测。DU8ML计算结果显示,结构2的¹³C化学位移均方根偏差(rmsd)为1.38 ppm,而结构3的rmsd较高,为1.59 ppm。进一步通过HOSE 和DU8ML方法将实验数据与计算值进行对比(表S1,见SI),发现原结构4的化学位移预测偏差达7–9 ppm,而修正结构2的偏差不超过3 ppm。此外,DP4概率计算表明,结构2的正确性概率为100%。
针对Zhang等11报道的ROESY谱,结构4中观察到H-1与H-9、H-13之间的相关信号。已知当相互作用质子间距小于5 Å时,NOESY/ROESY相关信号可能被观测到20。基于DU8ML计算结果(见SI),我们计算了结构2低能构象中相关质子对的玻尔兹曼加权平均距离:H-1与H-13间距为3.46 Å,H-1与H-9间距为3.69 Å。这一结果进一步表明结构2与报道的ROESY谱数据相符。
此外,通过将计算ECD谱与zerumol A (1)的实验谱图11进行对比,确定了结构2的绝对构型。为此,我们采用含时密度泛函理论(TD-DFT)对修正结构2的电子圆二色谱(ECD)进行了模拟计算。计算在TD-B3LYP/6–311+G(d,p)理论级别下完成,并使用极化连续介质模型(PCM)以甲醇为溶剂模拟溶剂效应。几何优化和溶剂校正均在B3LYP/6–311+G(d,p)级别进行。最终ECD谱由七个低能构象的玻尔兹曼加权平均得到:对每个对映体构象,通过叠加20个激发态的计算跃迁并采用0.3 eV高斯展宽及12 nm波长校正获得其ECD贡献。TD-DFT模拟中应用波长校正在文献中已有成熟先例(参见21等)。图7显示,zerumol A修正结构的计算曲线与实验ECD谱高度吻合,从而证实原结构1应修正为顺式-(R,R)-异构体6。

图7. Zerumol A的实验ECD谱(蓝色实线)与修正结构6的TD-DFT计算ECD曲线(黑色虚线)及原误判结构1的计算ECD曲线(红色实线)对比图。原始ECD数据转载自Fitoterapia, 2025, 184, 106662, Z. Zhang, G. Li, X. Luo, F. Xiao, Y. Zhang, C. Tang, B. Hu, Z. Xie, S. Yao, Y. Ye. Diterpenoids and diarylheptanoids from the fruits of Alpinia zerumbet and their anti-inflammatory activity © 2025,经Elsevier授权使用。
综上,本研究将zerumol A的结构重新确认为结构6,其氧杂环丁烷部分实为含有顺式-(R,R)-立体构型的5,6-二氢-2H-吡喃环系。这一结论主要得益于计算化学方法的协同应用,特别是ACD/Structure Elucidator与DFT DU8ML的有机结合。
