企业官网建设流程全解析

关键词：GROMACS；分子动力学；直链淀粉；双螺旋；淀粉糊化

一、文章简要介绍

米饭煮久了会糊，淀粉在热水里慢慢失去结晶结构，这个过程叫糊化。糊化的分子机制一直没完全讲清楚。这篇论文用GROMACS分子动力学模拟，专门看A型直链淀粉双螺旋在水里加热时稳不稳。他们搭了六个体系，链长有6、12、24个葡萄糖单元，还分单根双螺旋和六根双螺旋堆在一起的十二聚体两种组合。每个体系跑5次独立的1微秒模拟。温度主要设在100度，另有一组40度做对比。论文发表在Scientific Reports，DOI：10.1038/s41598-026-66415-2。

图1 六个淀粉寡聚体模拟体系：6、12、24个葡萄糖单元，单根双螺旋或六根双螺旋堆积的十二聚体

图1就是六个体系的构成。葡萄糖单元越短，体系越小，水分子数从2176到38077不等。链长、根数、盒子形状这些差别，就是后面稳定性差异的源头。

二、仿真步骤

步骤一：搭六个淀粉体系。模型从A型淀粉晶体结构出发，用doGlycans工具把GLYCAM06力场参数分配给糖分子，水模型用TIP3P。短链体系用截断八面体盒子，24个葡萄糖的长链因为形状细长，改用三斜盒子。

图2 373K下RMSD随时间演化：每个体系5次独立1微秒重复，O4原子拟合到晶体结构

图2是5次独立1微秒模拟的RMSD曲线。判断结构崩没崩，就看这条线涨不涨。

步骤二：GROMACS跑微秒级模拟。软件用GROMACS 2021，静电用PME方法，截断半径10埃，范德华截断也是10埃，所有键长用P-LINCS约束。流程是先能量最小化，再固定淀粉原子分别跑100皮秒的NVT和NPT平衡，让周围水分子和盒子密度先松弛好，然后放开约束做生产模拟。生产阶段用NPT系综，温度用随机速度重标度控制，时间常数0.1皮秒，压力用Parrinello-Rahman压浴，时间常数2皮秒，步长2飞秒。每个体系用不同初始速度跑5次独立1微秒。

图3 6Glc_dodecamer稳定性：(A)六根双螺旋各自的RMSD，(B)与两根以上双螺旋接触的水分子数，(C)1微秒终态结构

图3是短链十二聚体的崩塌过程。看A图，六根双螺旋的RMSD在0.1微秒前就开始往上走，B图的水分子数在0.4到0.8微秒才明显增加。结构先乱，水后进来，这个顺序很有意思。

步骤三：怎么判断双螺旋散没散。主要看RMSD，论文只取O4原子，也就是形成糖苷键的氧原子，先拟合到初始结构再算偏差。RMSD涨到0.5纳米以上就算明显偏离。另外统计与两根以上双螺旋同时接触的水分子数，衡量水往寡聚体内部渗了多少。补充材料里验证了吡喃环构象和糖苷键扭转角没有出现异常状态，螺旋内的静电和范德华能变化也跟RMSD走势一致，说明RMSD这个指标靠谱。

三、关键结果与发现

第一个结果：六个葡萄糖的短链扛不住100度。6个葡萄糖的体系，不管是单根双螺旋还是十二聚体，在100度水里都塌了。dodecamer的RMSD在0.1微秒前就开始明显偏离初始结构，到1微秒时双螺旋基本散开。水分子数的激增主要发生在模拟中后期，说明短链是先失去有序结构，水再大量涌入。这个时序和实验上观察到的糊化过程对得上。

图4 12Glc_dodecamer稳定性：(A)六根双螺旋各自的RMSD，(B)内部水分子数，(C)1微秒终态结构

图4是12个葡萄糖的十二聚体，六根双螺旋的RMSD基本稳在低位，只有末端能看到局部松散，内部水分子数也没有明显上涨。

第二个结果：十二、二十四个葡萄糖的十二聚体能撑住。12个和24个葡萄糖的十二聚体，在1微秒里基本保持堆积结构，只是末端出现fraying，也就是双螺旋端头散开几圈。但同样12个葡萄糖的孤立双螺旋dimer在100度下却塌了，24个葡萄糖的dimer偏差也明显比dodecamer大。这说明六根双螺旋并排堆在一起，彼此之间的横向接触本身就是一种稳定因素，不能拿孤立双螺旋的行为去推断堆积体系。

图5 十二聚体1微秒终态MD结构：左侧棒状模型、右侧表面模型，水分子以球体显示

图5是十二聚体在1微秒时的终态结构，棒状和表面两种画法。长链体系整体还保持着双螺旋堆积，水分子主要在表面和缝隙里。

第三个结果：链长和堆积决定谁先崩。40度的对比模拟里，孤立12Glc_dimer保持有序，但6Glc_dodecamer照样塌。合起来看，链长和横向堆积都影响有序结构在水里的存活时间。支链淀粉结晶片层里的无分支段一般认为有12到16个葡萄糖单元，12Glc和24Glc十二聚体相对能撑，正好跟这个区间对应，算是给较长的无分支段为什么更抗糊化提供了一个分子层面的解释。

图6 糖苷键扭转角分布（补充图S2）：模拟中没有出现新的扭转角状态

图6是糖苷键扭转角分布，模拟轨迹里没有出现晶体结构之外的新扭转角状态，排除了构象采样异常的干扰。

图7 6Glc_dodecamer双螺旋内静电与范德华相互作用能（补充图S3）：变化趋势支持RMSD作为解旋指标

图7是双螺旋内部的静电和范德华相互作用能随时间的变化，和RMSD的涨落能对上，说明RMSD反映的就是双螺旋本身的解旋，不是别的因素。

四、我们提供的服务

我们可提供以下仿真服务：①GROMACS分子动力学模拟（含糖类、蛋白、核酸等生物大分子体系）；②力场参数分配与体系搭建（GLYCAM06、doGlycans等工具链）；③RMSD、氢键、水渗透、相互作用能等轨迹分析；④淀粉、多糖等材料的微观机制研究。如果你也在研究这类问题，欢迎私信聊聊。