多肽自组装模拟：二苯丙氨酸纳米管成核自由能与临界核尺寸

本项目用 well-tempered metadynamics 模拟二苯丙氨酸（diphenylalanine, FF）在 12 mM 有效浓度下的自组装成核过程。体系 48 个肽，温度 300 K，序参量 Q 为 β-sheet 含量，偏置因子 12，模拟 22 ns 液相。分析层间堆垛、二级结构含量演化与经典成核自由能曲线。

案例分析

FF 肽形成反平行 β-sheet 堆叠的纳米管结构，层间距约 4.72 Å，与实验 XRD 一致。二级结构含量在 22 ns 内从随机卷曲迅速上升至 β-sheet 含量 0.81。经典成核理论拟合自由能曲线得到临界核尺寸 n*=5 个肽，成核势垒 ΔG* = 11.1 kJ·mol⁻¹。超过临界核后自由能单调下降，组装自发进行。层间堆垛以疏水苯环 π–π 相互作用为主，氢键稳定 β-sheet 骨架。

模型细节、验证与工程意义

模拟在 GROMACS 中运行，MARTINI 力场，well-tempered metadynamics 用 PLUMED 2.9 插件实现。序参量 Q 为 β-sheet 氢键数与相邻肽对数的综合指标。自由能曲线通过 reweighting 获得。验证时层间距 4.72 Å 与 FF 纳米管 XRD 峰值吻合；临界核 n*=5 与超速离心沉降实验测得的低聚体分布趋势一致。

总结教训

第一个教训是 metadynamics 偏置因子选择要平衡探索与精度。偏置因子 12 在本项目中使自由能在 22 ns 内收敛。第二个教训是临界核尺寸对浓度敏感，12 mM 下 n*=5；浓度升高会减小 n*。第三个教训是序参量 Q 必须同时反映聚集数与二级结构，仅用 RMSD 会误判无定形聚集体为核。

项目经验

1. 多肽自组装模拟二苯丙氨酸 FF 体系。
2. 48 个肽，12 mM 有效浓度，300 K。
3. 方法：well-tempered metadynamics，偏置因子 12。
4. 纳米管层间距约 4.72 Å。
5. 反平行 β-sheet 堆叠结构。
6. 22 ns 液相中 β-sheet 含量达 0.81。
7. 临界核尺寸 n* = 5 个肽。
8. 成核势垒 ΔG* = 11.1 kJ·mol⁻¹。
9. 自由能曲线符合经典成核理论。
10. 层间 π–π 相互作用驱动堆垛，氢键稳定骨架。
11. 偏置因子 12 可在 22 ns 内收敛自由能。
12. metadynamics 由 PLUMED 2.9 实现。
13. 自由能通过 reweighting 获得。
14. 本项目方法可迁移到淀粉样纤维、自组装水凝胶设计。

补充经验

1. 溶剂极性会显著影响 π–π 相互作用强度，从而改变层间距。
2. pH 影响 FF 肽的电荷状态，低 pH 有利于自组装。
3. 不同浓度下需重新计算临界核与势垒，不可直接外推。
4. 温度升高会降低成核势垒但增加临界核波动。

复现参数速查（案例图硬参数）

更多相关仿真经验，可参考本站的 MD 仿真案例栏目。