来，和我一起啃GROMACS在真实科研中的这些硬骨头

导读：大家好，我是刘建川，网名刘十三，做GROMACS分子动力学模拟的研究和教学十四年了。自2023年8月以来，我在仿真秀开通分子动力学仿真专栏这段时间，从GROMACS零基础入门课做起，慢慢发现一个很普遍的问题——很多朋友把基础的蛋白+水+离子体系跑通之后，一进到自己的科研课题就卡住了。膜蛋白怎么插膜、金属离子怎么处理、非标准残基力场怎么建、PMF自由能怎么算……这些都是真实科研里的硬骨头，基础教程覆盖不到。

刘建川主讲GROMACS分子动力学模拟线下培训

自今年开始，我会陆续推出GROMACS系列专题线上课程，来啃一个个硬骨头。就在上周，我一口气把四个专题，在仿真秀官网全部上架了，下面跟大家详细说说每门课讲什么、能解决什么问题。总而言之，这4 门专题课，就是围绕这些真实科研体系中的技术难点展开，欢迎关注在仿真秀官网搜索【刘建川】。

线下培训全年招生-可在仿真秀专栏报名

专题一：生物膜、膜蛋白的构建和模拟

这门课的重点不是"用网站生成一层膜"，而是处理膜体系真正容易出问题的环节。课程内容包括：

（1）混合脂质膜的构建，组分和比例怎么调，溶剂化后进入膜疏水区域的水分子怎么修正

（2）上下膜层分别建立索引、膜体系的semiisotropic压力耦合处理

（3）CHARMM-GUI膜蛋白建模完整流程

（4）packmol-memgen构建不对称膜

（5）Amber体系转GROMACS的处理方法

膜怎么建只是第一步，更重要的是膜蛋白放得对不对、疏水区合不合理、体系怎么稳定进入MD。这门课把膜体系从建模到平衡的整条路线讲透。

专题二：含有金属配体的蛋白建模和模拟

金属蛋白建模最重要的，不是先找一个Fe或Zn的参数，而是先判断金属属于哪一种状态。课程里把金属体系明确分成三类：

（1）非键模型：金属与蛋白不成键、可以自由运动

（2）配位限制模型：金属不成键，但必须保持在结合位点

（3）成键模型：金属与附近氨基酸存在明确的键连关系

对应三种模型分别怎么处理——非键模型怎么设、配位限制怎么加、成键模型的bond/angle/dihedral参数怎么写进ITP。不是所有金属都能套同一个模板，而是先判断"自由、限制还是成键"，再决定拓扑怎么写。

专题三：质心牵引、离子（分子）跨膜计算PMF

这门课从普通MD进入反应坐标与自由能计算，是很多做膜体系、做转运的朋友必须掌握的技术。完整技术路线：质心牵引→构型提取→伞形窗口→多窗口模拟→WHAM→PMF→ΔG

课程覆盖：

（1）distance、direction、direction-periodic三类Pull几何方式的区别和适用条件

（2）连续牵引之后怎么沿反应坐标提取构型

（3）伞形采样窗口怎么建立，怎么保证相邻窗口有合理重叠

（4）多窗口模拟的批量运行方法

（5）用gmx wham重建PMF、分析自由能变化

课程还扩展到离子/小分子跨膜PMF，把目标离子替换成研究分子，就能建立分子—膜反应坐标，研究分子穿过膜不同区域时的自由能变化。

本课程一套围绕 GROMACS Pull Code、伞形采样和自由能分析设计的专题实战课程。课程基于 GROMACS 2026.3 以上版本，从质心牵引的基本逻辑开始，逐步完成蛋白体系的牵引、反应坐标生成、伞形采样、多窗口计算、WHAM 分析，并进一步扩展到离子和小分子跨膜 PMF 计算。课程重点不是只给一套 Pull 参数，而是从“怎么拉”，一直做到“怎么得到 PMF 和自由能结果”。以下是课程大纲

第 1 讲｜质心牵引 Pull‑Code 原理与 mdp 参数详解

（1）Pull‑Code 质心牵引三种工作模式：拉伸、约束、恒力

（2）pull 相关 mdp 全部关键参数解析，牵引组、参考组概念

（3）反应坐标定义，跨膜体系索引分组构建

（4）膜体系预处理，为 PMF 计算做准备

第 2 讲｜Steered‑MD 质心拉伸模拟，生成伞形采样窗口构型

（1）拉伸模拟 mdp 配置，牵引速率、力常数参数选择

（2）运行拉伸 MD，沿跨膜反应坐标拉动离子 / 小分子

（3）从拉伸轨迹批量提取各个伞形采样窗口构型

（4）窗口间隔设置原则，构型合理性检查

第 3 讲｜伞形采样 Umbrella Sampling 全流程实操

（1）每个采样窗口 mdp 文件配置要点

（2）批量生成各个窗口 tpr 输入文件

（3）伞形采样模拟运行，模拟过程监控与报错排查

（4）pull‑xvg 输出轨迹文件解读，评估每个窗口采样质量

第 4 讲｜gmx wham 后处理，重构 PMF 自由能曲线与结果校验

（1）wham 所需 tpr‑files.dat、pullf‑files.dat 输入文件准备

（2）gmx wham 命令实操，输出 PMF 平均力势曲线

（3）PMF 曲线解读、单位转换，结果可靠性判断

（4）常见异常曲线原因分析与调参思路；PMF 方法迁移拓展

课程覆盖知识点清单

✅GROMACS Pull‑Code 质心牵引原理与 mdp 参数 

✅牵引组、参考组、反应坐标索引构建 

✅Steered‑MD（拉伸分子动力学）实操

✅伞形采样窗口构型生成策略 

✅伞形采样模拟参数设置与批量处理

✅pull‑xvg 轨迹文件读取与采样质量评估 

✅gmx wham 工具完整使用 

✅PMF 平均力势曲线输出、解读与校验 

✅离子 / 小分子跨膜自由能计算完整工作流 

✅PMF 计算常见报错、曲线畸变问题排查

专题四：非标准残基力场构建

这门课解决的是蛋白模拟里非常典型的一类问题：pdb2gmx根本不认识这个残基怎么办？当然，GROMACS中的非标准残基并不是简单生成一个ITP就结束。如果力场中没有对应构建单元，就需要真正理解并扩展.rtp、.hdb、residuetypes.dat，以及pdb2gmx对力场数据库的读取逻辑。

课程完整流程：

（1）提取非标准残基、封端与质子化处理

（2）AmberTools/acpype参数化，生成ITP/MOL2

（3）构建RTP和HDB文件

（4）更新residuetypes.dat

（5）重新执行pdb2gmx，验证非标准残基能否被正确识别

不是得到一个孤立的小分子ITP，而是让非标准残基真正重新接回蛋白主链，成为GROMACS力场能够识别的残基。

四门课，四个核心能力：复杂膜体系构建与平衡、金属非键/限制与成键模型判断、从质心牵引做到自由能分析、自主扩展GROMACS蛋白力场数据库。此外，在大家学习过程中，我还提供课程VIP群答疑服务，欢迎加入我们。

踏踏实实深耕技术，认认真真做好教学。专研GROMACS，只做真实战。有需要的朋友，可以到我的仿真秀专栏看看，更多课程详情和目录都在里面，请点击文章末尾阅读原文关注。有问题也欢迎交流。

2026研发仿真开学季是仿真秀面向研发工程师推出的主题活动月，9月1日至10月31日以"从工具到智能——让AI接入仿真研发"为主线，覆盖结构、流体、电磁、多物理场和AI驱动仿真五个方向。仿真秀力邀企业资深研发工程师，围绕LS-DYNA、Workbench、Abaqus、Star-CCM+、HyperWorks、Maxwell、Fluent、CST、Flotherm、SolidWorks、ANSYS等主流软件，演示如何用LLM、Agent、Codex、Claude、DeepSeek等大模型能力打通建模、求解、后处理与优化全流程，让仿真从"手工操作"升级为"AI驱动"；

此外仿真秀精选100余门研发仿真精品视频课限时秒杀1元起售，并叠加专属券，秒杀课程覆盖从基础到前沿的完整学习路径；活动期间，仿真秀还发起案例征集，以"我的研发仿真故事——从工具到智能的进阶之路"为题，面向工程师征集攻坚实录、工具进化、AI接入、成长与团队四类真实故事，图文或视频，通过即获满意现金稿酬和百人千万"内容扶持计划，优秀作品还将在仿真秀APP、官网、公众 号、视频号四端专题展示并邀请连麦直播。整个主题月既是一场系统化的技术充电，也是工程师用AI重塑研发工作流的一次集中练兵。