刘建川・Gromacs 实战｜如何获取并修改分子 ITP 文件中的原子电荷
对于蛋白体系，通过 pdb2gmx 生成的原子电荷基本合理，一般无需检查和修改；但小分子不一样，不同力场生成ITP的方式差异很大，经常出现原子电荷不合理、总电荷不为0、与分子实际带电量不匹配等问题，这种情况就需要我们手动计算并修改ITP电荷。
对于刚性分子模拟，MK、CHELPG 静电势拟合电荷完全够用；但对于存在多构象的柔性分子，MK/CHELPG电荷存在明显短板：构象依赖性强、单一构象拟合无法体现原子等价性、包埋原子电荷拟合偏差大。
而Kollman团队提出的RESP约束静电势电荷，是目前最适配柔性小分子分子动力学、构象分析、分子对接的原子电荷方案，完美解决了传统MK/CHELPG电荷的诸多缺陷。下面手把手教大家完整获取、拟合并替换RESP电荷。
一、高斯软件计算ESP静电势电荷
1. 分子文件预处理
用GaussView打开生成Top/ITP所用的PDB文件，先另存为 mol2格式，再用GaussView重新打开mol2文件，导出为 gjf 输入文件； 这样可以避免原PDB中的残基名称冗余信息写入gjf，保证量化计算纯净性。
2. 分子结构优化
用记事本打开分子 gjf 文件，修改基组与计算任务： 前三周期元素推荐基组：6-311G** 含第四周期及重原子：重原子采用SDD赝势+配套赝势基组
以DMSO分子为例，结构优化任务关键词设置：
#p opt b3lyp/6-311G**
删除gjf中多余的Windows路径与chk目录冗余配置，保存后提交高斯计算。
Linux系统下高斯运行命令：
g09 dmso.gjf
计算完成后生成 log 日志文件，可用GaussView直接打开优化后的分子结构，另存为 dmso_resp.gjf，用于下一步静电势电荷计算。
3. 静电势（ESP）电荷计算
拟合RESP电荷需在HF/6-31G(d) 级别下计算静电势，HF理论可适度高估偶极矩，等效隐含溶剂极化效应，是拟合标准范式。
电荷计算关键词设置：
#p HF/6-31G(d) Pop=MK Iop(6/33=2,6/41=10,6/42=17)
保存gjf文件后再次提交高斯计算：
g09 dmso_resp.gjf
运行结束得到 dmso_resp.log 量化输出文件，作为后续RESP电荷拟合的源文件。
二、AmberTools 拟合标准RESP电荷
调用AmberTools中 antechamber 工具，一键拟合RESP电荷，终端输入以下命令：
antechamber -i dmso_resp.log -fi gout -o dmso_resp.mol2 -fo mol2 -c resp -nc -1 -pf y
运行完成生成 dmso_resp.mol2 文件，原子电荷数据直接标注在文件中 MOL 字段后方，可直接提取使用。
三、手动替换Top/ITP文件原子电荷
用记事本打开小分子对应的 top或itp拓扑文件；从 dmso_resp.mol2中逐一复 制拟合好的RESP原子电荷；可利用记事本列复 制模式，批量批量替换原有错误电荷； 替换完成后核对分子总电荷，确保与理论带电量完全一致。
四、带电离子体系专属注意事项
若计算对象是带电离子体系，不能直接默认中性结构： 需要先在GaussView中手动设置分子体系带电量，再依次进行结构优化 → ESP静电势计算 → RESP电荷拟合，否则拟合出的电荷总量会严重出错，影响后续分子模拟精度。
整套流程从高斯量化结构优化、静电势计算，再到Amber拟合RESP电荷、手动替换ITP拓扑电荷，是小分子力场构建的标准实操方案，适配AMBER、GAFF、OPLS等主流力场。熟练掌握后，就能彻底解决小分子电荷错乱、总电荷不守恒、模拟相互作用失真等常见难题。
