黄芪茎多糖代谢物与DAO结合GROMACS模拟

关键词：GROMACS；分子动力学；MM-PBSA；黄芪多糖；分子对接

一、文章简要介绍

黄芪的根是常用中药，茎秆却大量废弃。这篇论文把黄芪茎里的多糖提出来，喂给泌乳中期的奶山羊，看能不能改善抗氧化能力和产奶性能。作者先做提取条件优化，再做60天饲喂试验，然后上血清代谢组学找差异代谢物，最后用分子对接和GROMACS分子动力学验证关键代谢物与靶蛋白的结合。论文发表在Scientific Reports 2026年第16卷12762篇，DOI：10.1038/s41598-026-39922-5。

图1 研究流程图：提取优化、饲喂试验、代谢组学、对接与MD验证

图1是整个研究的流程。从黄芪茎预处理、热水提取多糖，到奶山羊饲喂、血液和乳样分析，再到代谢组学筛选、对接和动力学验证，一条线走下来。

二、仿真步骤

步骤一：锁定代谢物和靶蛋白。血清代谢组学（UPLC-MS/MS）筛出307个差异代谢物，134个上调、173个下调，富集在甘油磷脂代谢、D-氨基酸代谢等20条通路。挑出6个关键代谢物：LPC(22:3)、L-鸟氨酸、L-色氨酸、咖啡酸、莽草酸和苯乙酰-L-谷氨酰胺，分别对应PLA2G2A、OTC、IDO1、SHMT1、SHMT1和DAO六个靶蛋白。

步骤二：AutoDock Vina对接。靶蛋白结构从PDB取人源同源结构（羊没有晶体结构，用人源做代理）：3U8I、1OTH、8ABX、1BJ4、1BJ4、7U9U。配体从PubChem拿SDF，用AutoDock Tools转成PDBQT。对接结果显示六个复合物结合能-4.687到-8.339千卡每摩尔，多数小于等于-5.0。最强的是苯乙酰-L-谷氨酰胺对DAO，-8.339千卡每摩尔，三个氢键撑着（Gly312、Leu316、Thr317）。这个复合物被选去做分子动力学。

图2 六个关键代谢物与通路对应蛋白的分子对接：LPC(22:3)-PLA2G2A、鸟氨酸-OTC、色氨酸-IDO1、咖啡酸-SHMT1、莽草酸-SHMT1、苯乙酰谷氨酰胺-DAO

图2是六个对接pose。每个代谢物都待在各自靶蛋白的口袋里。苯乙酰-L-谷氨酰胺与DAO的结合最牢，后面所有动力学分析都围绕它展开。

步骤三：GROMACS搭建体系。蛋白用AMBER14SB力场，配体用Antechamber算AM1-BCC电荷、GAFF2原子类型，再经ACPYPE转成GROMACS拓扑，离子参数走Joung-Cheatham。复合物放进截断八面体盒子，蛋白表面到盒边至少1.2纳米，TIP3P水溶剂，补0.15摩尔每升氯化钠。先最陡下降能量最小化，力小于1000千焦每摩尔纳米收敛，再分两段各200皮秒做NVT和NPT平衡，温度298K。

步骤四：100纳秒生产模拟。NPT系综，2飞秒步长，Verlet积分，PME处理长程静电，范德华和库仑截断1.2纳米，氢键用LINCS约束，Nosé-Hoover控温、Parrinello-Rahman控压，298K、1巴。轨迹每10皮秒存一帧，分析用GROMACS自带工具加VMD和PyMOL，结合自由能用gmx_MMPBSA算。

三、关键结果与发现

结果一：复合物很稳。骨架RMSD大约20纳秒后稳定在0.38加减0.05纳米，配体RMSD始终低于0.10纳米，空蛋白稳定在0.25纳米（图3A）。回转半径Rg在2.75纳米附近小幅波动，整体紧凑度没变（图3B）。RMSF显示308到322号残基的环区活动大，最大波动约0.60纳米（图3C）。配体和结合口袋的质心距离恒定在1.45加减0.03纳米，SASA均值6加减0.5平方纳米，界面稳定（图3D-E）。

图3 GROMACS分子动力学模拟总览：A-O共15个子图，涵盖稳定性、能量与构象分析

图3是MD全程的总览，A到O十五个子图。RMSD收敛、Rg平稳、氢键维持在3个左右、PCA只采样两个主导构象簇，这些信号指向同一个结论：苯乙酰-L-谷氨酰胺在DAO口袋里待得住。

图4 复合物稳定性分析：(A)RMSD，(B)Rg，(C)RMSF，(D)配体-催化位距离，(E)SASA

图4把稳定性指标单独放大看。A到E分别是RMSD、回转半径、残基波动、配体距离和溶剂可及面积，从头到尾没有漂移。

结果二：能量上有利。MM-PBSA算出的总结合自由能是负92.051加减0.366千焦每摩尔。拆开看，范德华项贡献最大，负158.225千焦每摩尔，静电项负37.614，非极性溶剂化负17.688，极性溶剂化正121.475抵消掉一部分，熵项正8.288，修正后结合自由能负83.764。范德华主导，疏水作用跟上来，静电帮忙但不占大头（图3G-H）。

结果三：关键残基和氢键网络。残基分解指出Tyr224和Tyr314是两个热点，π-π堆积、T形π-π和π-烷基相互作用挂在Tyr224和Tyr228上，范德华接触涉及Ile230、Tyr314、Thr44。全程氢键数平均3个，具体接触有六个：Gly312、Gly313、Gly315、Leu316、Thr317和Leu51（图3I-J）。静电表面显示结合沟带正电，正好配配体的极性基团（图3L）。自由能景观只有一个深盆地，体系就一个稳定构象（图3M）。

图5 结合能量与构象分析：(G)真空能量，(H)MM-PBSA残基分解，(I)氢键数，(J)3D视图，(K)PCA，(L)静电表面，(M)自由能图

图5是能量和构象部分。H图的MM-PBSA残基分解里，Tyr224和Tyr314的柱子明显高过其他残基，这两个位点是将来做定点突变的候选。

四、我们提供的服务

我们可提供以下仿真服务：①GROMACS全流程分子动力学（建模、平衡、生产、分析）；②AMBER力场配体参数化（Antechamber、ACPYPE）；③MM-PBSA结合自由能计算与残基分解；④AutoDock Vina分子对接；⑤代谢组学与计算模拟联合分析。需要做药物-靶点或代谢物-酶结合验证的，欢迎私信聊聊。

如需相关服务，欢迎通过公众号“320科技工作室”联系我们。