进行蛋白质-配体对接时出现 no such moleculetype ligand

admin

整理一个实际报错场景，看看大家一般会怎么排查。

操作流程如下：1. （进行蛋白质拓扑）选择的力场为6，即amber99sb-ildn2.获得3nsn.gro posre.itp topol.top文件3.使用top处理配体文件得到 ligand.gro ligand.top ligand.itp文件4.将3nsn.gro 复制并重命名为complex.gro，然后将ligand.gro坐标部分复制粘贴到complex.gro……

关键命令：

1. gmx pdb2gmx -f 3nsn.pdb -o 3nsn.gro -water spc -ignh
2. gmx editconf -f complex.gro -o newbox.gro -bt dodecahedron -d 1.0

复制代码

核心报错/现象：

1. no such moleculetype ligand

复制代码

主要想确认：这个现象通常应该优先检查哪些文件或参数，怎样定位根因？

admin

这个报错通常不是坐标文件本身的问题，而是 topol.top 里没有正确读到配体的 [ moleculetype ]。建议先检查两处：一是 ligand.itp 是否在 [ system ] 和 [ molecules ] 之前被 #include；二是 [ molecules ] 中写的名称是否与 ligand.itp 里 [ moleculetype ] 第一列完全一致。

可以先让 grompp 展开拓扑看看：

1. gmx grompp -f em.mdp -c newbox.gro -p topol.top -pp processed.top

复制代码

再在 processed.top 里搜索 ligand。如果搜不到，基本就是 include 路径或 include 顺序的问题。