计算化学公社

 找回密码 Forget password
 注册 Register

求助:蛋白配体复合物模拟过程中用acpype处理小分子mol2文件报错

查看数: 2494 | 评论数: 10 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-11-22 23:48

正文摘要:

请问各位老师: 我在使用acpype生成小分子的拓扑文件,出现下面的报错,请问该如何解决(附带有我的配体分子结构和mol2文件) 【跟着sob老师在北京的课件,跑过3ATL,里面BEN.mol2这一步是完全可以跑出来的】 === ...

回复 Reply

wanwan5339 发表于 Post on 2024-5-11 01:01:39
sobereva 发表于 2024-5-11 00:25
当然会变。难道MD过程中构象都不变?为什么不能变?
普通有机配体根本就不需要基于Hessian去计算力常数 ...

好的,十分感谢老师解答我的疑惑
sobereva 发表于 Post on 2024-5-11 00:25:37
wanwan5339 发表于 2024-5-11 00:23
很抱歉,老师。学生笨拙,我还想问您一下,您在这里提到“以对接得到的结构为初猜结构做个常规的配体分子 ...

当然会变。难道MD过程中构象都不变?为什么不能变?
普通有机配体根本就不需要基于Hessian去计算力常数,认真把sobtop主页的例子和FAQ部分看完整
wanwan5339 发表于 Post on 2024-5-11 00:23:37
sobereva 发表于 2024-5-10 16:56
不改变坐标还优化什么
不想改变坐标就不优化
用什么结构算RESP电荷,441博文末尾该说的都说了

很抱歉,老师。学生笨拙,我还想问您一下,您在这里提到“以对接得到的结构为初猜结构做个常规的配体分子的几何优化,然后算RESP电荷”,这样的话不还是需要Gaussian进行几何优化吗,如果进行几何优化我的配体坐标和构象还是会改变呀?
我现在还是不明白这个问题。而且我在sobtop里bonded项会基于Hessian计算,这个时候还需要fchk文件,如果不进行结构优化,我也没办法得到这个文件。

优化.png (67.72 KB, 下载次数 Times of downloads: 11)

优化.png
sobereva 发表于 Post on 2024-5-10 16:56:25
wanwan5339 发表于 2024-5-10 11:51
老师,我看完后还是不太明白。
我对接得到的结构会因为分子优化而改变它的坐标以及构象。这个时候我生成 ...

不改变坐标还优化什么
不想改变坐标就不优化
用什么结构算RESP电荷,441博文末尾该说的都说了
wanwan5339 发表于 Post on 2024-5-10 11:51:53
本帖最后由 wanwan5339 于 2024-5-10 11:52 编辑
sobereva 发表于 2024-5-10 07:31
如果是让sobtop指认GAFF原子类型和GAFF的力场参数,完全不依赖于几何结构,仅看连接关系
RESP电荷会受到 ...

老师,我看完后还是不太明白。
我对接得到的结构会因为分子优化而改变它的坐标以及构象。这个时候我生成的RESP电荷和拓扑文件都是基于优化后的分子构型坐标得到的。这样我的对接结果就没有意义了。
我对接的结果该拿什么软件进行分子优化才能不改变它的坐标呢?还是说我生成拓扑文件后,在文件里改成对接的时得到的坐标?
sobereva 发表于 Post on 2024-5-10 07:31:38
wanwan5339 发表于 2024-5-9 00:21
老师,使用sobtop之前,需要对我的结果进行优化,但是我的结构需要有特定的坐标(分子对接得到的坐标), ...

如果是让sobtop指认GAFF原子类型和GAFF的力场参数,完全不依赖于几何结构,仅看连接关系
RESP电荷会受到构象影响,用什么结构看此文最后一节
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441http://bbs.keinsci.com/thread-10880-1-1.html

wanwan5339 发表于 Post on 2024-5-9 00:21:28
sobereva 发表于 2022-11-23 00:08
acpype已经过时了,用sobtop强得多,看http://sobereva.com/soft/Sobtop

早期的北京科音分子动力学与GRO ...

老师,使用sobtop之前,需要对我的结果进行优化,但是我的结构需要有特定的坐标(分子对接得到的坐标),经过分子优化,坐标必定发生改变,这个时候我继续生成拓扑文件,能不能在生成.gro等文件里把坐标替换成我分子对接里得到的坐标呢?
喵星大佬 发表于 Post on 2022-11-23 10:50:10
acpype对于未明确定义的成键项参数怎么生成的,和sobtop猜的似乎差不少
松果seven 发表于 Post on 2022-11-23 09:40:56
sobereva 发表于 2022-11-23 00:08
acpype已经过时了,用sobtop强得多,看http://sobereva.com/soft/Sobtop

早期的北京科音分子动力学与GRO ...

好的,谢谢老师
sobereva 发表于 Post on 2022-11-23 00:08:21
acpype已经过时了,用sobtop强得多,看http://sobereva.com/soft/Sobtop

早期的北京科音分子动力学与GROMACS培训班里用acpype的地方,在最新一期里大多都已经替换成了sobtop了

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

GMT+8, 2026-2-17 23:20 , Processed in 0.175294 second(s), 25 queries , Gzip On.

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