计算化学公社

 找回密码 Forget password
 注册 Register

求助:分子在.gro文件和.itp文件中原子顺序不同导致的原子名称不匹配

查看数: 6115 | 评论数: 6 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-3-3 17:14

正文摘要:

我做的是石墨基底上有机小分子自组装方向的模拟,我在MS上搭建了整个体系的模型,然后导出来放到vmd中生成的.gro坐标文件;有机小分子的.itp文件是用sob老师分享的网站 LigParGen Server (yale.edu)生成的oplsaa力场 ...

回复 Reply

Rosefinch 发表于 Post on 2022-3-3 22:52:12
sobereva 发表于 2022-3-3 22:36
基于ligpargen给的gro文件去建模(毫无必要非得用M$。packmol、insert-molecules都可以插入小分子,你自己 ...

好的,谢谢sob老师,我去试一下
sobereva 发表于 Post on 2022-3-3 22:36:32
基于ligpargen给的gro文件去建模(毫无必要非得用M$。packmol、insert-molecules都可以插入小分子,你自己在VMD手动挪动小分子亦可)

另外,用sobtop(http://sobereva.com/soft/Sobtop)产生小分子拓扑文件可以直接避免原子序号被重排的问题,sobtop产生拓扑文件的过程中不会去重排序号
Rosefinch 发表于 Post on 2022-3-3 21:48:57
牧生 发表于 2022-3-3 21:33
68个也不算很多。如果用chemdraw画分子,并存为pdb,确实有可能和ligpargen给出的顺序不一样,但是也没关 ...

ligpargen给出的itp和gro里的原子顺序确实是对应的,那请问如何用它给的这个单分子的gro文件去组装成我一开始搭建好的自组装模型呀,我这个模型只需要4个这个分子。
牧生 发表于 Post on 2022-3-3 21:33:59
Rosefinch 发表于 2022-3-3 21:26
感谢您的回复,但是我这个有机分子的原子数比较多,有68个,我按照您说的在Multiwfn中理过顺序之后,在li ...

68个也不算很多。如果用chemdraw画分子,并存为pdb,确实有可能和ligpargen给出的顺序不一样,但是也没关系,ligpargen会同时给出itp和gro,这二者一定是能对应起的。

如果实在强迫症要用pdb,那就gmx editconf -f XXX.gro -o XXX2.pdb转为pdb就可以了
Rosefinch 发表于 Post on 2022-3-3 21:26:24
牧生 发表于 2022-3-3 17:29
我以前也有类似错误,但现在不会犯了
看第九楼
http://bbs.keinsci.com/thread-24586-1-1.html

感谢您的回复,但是我这个有机分子的原子数比较多,有68个,我按照您说的在Multiwfn中理过顺序之后,在ligpargen网站上生成的.itp文件里的原子顺序依旧是它自己定义的一个顺序,和提交的.pdb文件中的原子顺序没有关系。我也不知道.itp里面的H原子什么的对应着哪一个H,请问您知道怎么解决嘛?
牧生 发表于 Post on 2022-3-3 17:29:48
本帖最后由 牧生 于 2022-3-3 18:36 编辑

我以前也有类似错误,但现在不会犯了
看第九楼
http://bbs.keinsci.com/thread-24586-1-1.html

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 19:28 , Processed in 0.235607 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list