刘建川・GROMACS 实战 | 生成CHARMM力场
一、Charmm力场基础介绍
Charmm力场由哈佛大学Martin Karplus教授研发，这位教授也是2013年诺贝尔化学奖得主。该力场适用范围极广，可用于有机分子、溶液体系、高分子聚合物、生物生化分子等研究场景。
除有机金属分子以外，Charmm力场能够给出高精度计算结果，可精准模拟分子结构、相互作用能、构型能、转动与振动能量、自由能等各类物理化学参数。
Charmm力场独有Urey-Bradley作用项，也是它区别于其他力场的核心特色，主要用于表征三重原子之间的键-角耦合振动效应。
二、Charmm高阶力场下载与安装
GROMACS自身已内置Charmm基础版本，若需要更高版本完整版Charmm力场，可进入官网下载： http://mackerell.umaryland.edu/charmm_ff.shtml#gromacs
下载对应版本压缩包解压后，将力场文件夹放入GROMACS力场目录： ~/share/gromacs/top/ 也可直接放置在模拟工作目录下，程序可自动识别调用。
三、小分子Charmm力场ITP/拓扑文件生成方法
生物大分子如蛋白、核酸，可直接使用GROMACS自带pdb2gmx模块一键生成拓扑；有机小分子则需要借助CGenFF网页工具生成，详细步骤如下：
1. 注册登录CGenFF官网
访问网址完成注册登录： https://cgenff.umaryland.edu/
2. 准备小分子mol2文件
使用GaussView搭建分子结构并保存为mol2格式； 注意：若文件中存在Ar原子标识，需手动替换为小写ar，无此项则无需修改。
3. 上传文件获取参数文件
登录CGenFF网站后，上传准备好的mol2文件； 无需勾选：Guess bond orders from connectivity 和 Include parameters that are already in CGenFF。
上传后页面出现注释行警告属于正常现象，无需处理。点击页面蓝色链接，复 制全部内容，新建文本文件（如CH2.str）粘贴保存。
4. 看懂CGenFF惩罚分数
CGenFF会为每一组力场参数给出惩罚分数，用来评估参数可靠性：
分数＜10：参数可靠，可直接用于科研模拟 10～50：需要人工验证拓扑结构与电荷分配 ＞50：参数偏差大，必须手动重新参数化
建议养成习惯，一定要核对原子类型、电荷分布，避免劣质拓扑影响模拟结果、浪费科研时间。
5. 下载转换脚本与力场包
网页顶部点击 More Info & Tools → 选择 Utilities； 进入后下载 GROMACS conversion program，得到脚本：cgenff_charmm2gmx.py。
同时从Charmm官网下载对应版本力场文件夹，例如 charmm36-jul2020.ff。 将mol2文件、str参数文件、python转换脚本、Charmm力场文件夹统一放在同一个工作目录。
6. 命令行一键生成GROMACS拓扑
先打开CH2.str文件，记下RESI后方的残基名称，在Linux或配置Python环境的Windows终端执行命令：
./cgenff_charmm2gmx.py Molecule CH2.mol2 CH2.str charmm36-jul2020.ff运行完成后，自动生成GROMACS可直接使用的top、itp、gro等全套结构与拓扑文件。
四、常见报错及完整解决办法
1. 缺少依赖模块报错
脚本运行依赖 numpy 和 networkx 库，未安装会直接报错，不同系统安装命令如下：
Ubuntu 系统：
pip install numpypip install networkxCentos 系统：
yum install numpyyum install python-networkx2. networkx版本兼容报错
出现报错：TypeError: add_node() takes exactly 2 arguments (3 given) 原因是高版本networkx不兼容脚本，需降级安装 1.x 旧版本。
下载地址：https://github.com/networkx/networkx/releases
以 networkx-1.10 为例，安装命令：
pip uninstall networkxtar -zxvf networkx-networkx-1.10.tar.gzcd networkx-networkx-1.10sudo pip install .按照以上整套流程操作，就能从零完成Charmm力场安装、小分子参数化、拓扑文件生成，产出的文件可直接投入GROMACS分子动力学模拟、自由能计算、构象分析等科研场景使用，完全满足日常科研建模需求。
