关键词GROMACS分子动力学MM-PBSA黄芪多糖分子对接一、文章简要介绍黄芪的根是常用中药茎秆却大量废弃。这篇论文把黄芪茎里的多糖提出来喂给泌乳中期的奶山羊看能不能改善抗氧化能力和产奶性能。作者先做提取条件优化再做60天饲喂试验然后上血清代谢组学找差异代谢物最后用分子对接和GROMACS分子动力学验证关键代谢物与靶蛋白的结合。论文发表在Scientific Reports 2026年第16卷12762篇DOI10.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分子对接⑤代谢组学与计算模拟联合分析。需要做药物-靶点或代谢物-酶结合验证的欢迎私信聊聊。
)
