多肽的分子动力学模拟：在溶剂、离子和膜环境中跑一条多肽链，水盒子里的每一颗钠离子都在改变构象分布

多肽分子动力学模拟是MD中最方便也最危险的一类体系。方便——一条20-40残基的肽链，加上水和离子，总共也就3-5万个原子，跑200ns在一张普通GPU上两天出结果。危险——多肽的构象景观是浅而宽的，自由能面上的极小值之间只隔了几个k_BT的高度，环境条件的细微偏移会让构象系综飘到完全不同的区域。

一个做抗菌肽膜作用机制的项目，用magainin 2（23个残基，GIGKFLHSAKKFGKAFVGEIMNS）在POPC脂双层中跑MD，第一轮结果看起来自洽——螺旋含量~65%，肽平行于膜表面。但对照实验的CD光谱和固态NMR数据：螺旋含量应该接近80%，而且肽的取向应该更倾斜地插入膜的疏水核心。

偏了15%的螺旋含量和十几度的取向角——不是实验误差，是仿真的初始条件选错了。

初始构象的偏置

magainin 2在纯水溶液里是无规卷曲为主，到了膜表面才折叠成两亲性α-螺旋。但初始构象应该用什么？做法一：从完全伸展的肽链开始——肽在膜表面自组装成螺旋，MD自己跑出二级结构。做法二：从NMR的螺旋模型开始——假设肽在膜上已经是螺旋构象，MD只对侧链和取向做采样。

做法一看似”更从头”，不需要实验输入——但200ns的时间尺度对于多肽折叠来说太短。magainin 2在水溶液中从无规卷曲折叠成螺旋的实验时间常数在微秒量级。200ns跑出来螺旋含量只有65%，是因为模拟还没走完折叠路径——不是平衡态螺旋含量就是65%。

切换到做法二：从NMR结构（PDB: 2MAG，脂胶束中两个螺旋段）出发，把肽嵌入预平衡的POPC双层中，螺旋轴初始平行于膜面。平衡200ns，螺旋含量稳定在78%——与CD实验一致。

初始构象偏置导致的15%螺旋含量差距需要强调：不是力场错了，是采样不足。如果做增强采样（如metadynamics或REMD），做法一也能收敛到接近80%——但要付出10-20倍的计算成本。

离子浓度在膜表面的分配

另一个容易忽视的变量是盐浓度。magainin 2带+4的净正电荷（4个Lys残基），在纯水里这些正电荷被Cl⁻反离子中和。但如果溶剂里加了150mM NaCl来模拟生理条件，Na⁺会在带负电的POPC膜表面积累形成一个离子双层。

这个离子双层的存在显著影响magainin 2的膜结合取向。Na⁺在膜表面的竞争性结合会部分屏蔽POPC头基的负电荷，减弱肽-头基的静电吸引力。在150mM NaCl下，magainin 2的肽-膜质心距离从纯水中的1.2nm增到了1.5nm——肽从”压”在膜上变成了”浮”在膜上。

这0.3nm的垂直位移对应着肽的疏水面和膜的疏水核心接触面积的变化——接触面积减小约25%，疏水驱动的螺旋稳定化效应随之减弱。无盐体系中magainin 2的螺旋含量是82%，150mM NaCl体系降到76%。

6%的差异——做抗菌肽的MD不指定盐浓度，不是在”简化模型”，是在改变物理。

力场的影响：CHARMM36 vs Amber ff19SB

最后也跑了一个力场交叉验证——同样的初始结构和模拟条件，用CHARMM36m力场跑出来螺旋含量78%，用Amber ff19SB跑出来71%。差别7%，属于肽模拟中两个主流力场的典型偏差范围。

差异的物理来源是力场对骨架Φ/Ψ二面角的偏好差异——CHARMM36m对α-螺旋区域的采样权重略高于Amber ff19SB。这不是哪个力场更”准”的问题——两个都在实验误差范围内（给定CD光谱的螺旋含量估算误差±8%），但如果你在螺旋含量70%和80%之间做机制推断（”这条肽的抗菌活性是孔形成还是地毯模型？”），这个差距就可能改变结论。

多肽的MD模拟要做对的不是某一个参数，而是三个互相牵制的选择：初始构象偏置决定你是否跑进了正确的自由能盆地；离子浓度决定肽-膜相互作用的强度；力场选择给了一个不可消除的±5-10%的系统不确定性。三个都对了，构象分布才能和实验在误差带内对齐。