计算化学公社

 找回密码 Forget password
 注册 Register
Views: 7362|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[Amber] amber MD骨架RMSD求助

[复制链接 Copy URL]

18

帖子

0

威望

111

eV
积分
129

Level 2 能力者

本帖最后由 Freshfish 于 2021-8-2 22:18 编辑

各位老师好,

根据amber教程A3 MM-PBSA计算中(中文版本AMBER高级教程A3:MM-PBSA|Jerkwin),单次运行时长为500ps,运行4次获得总时长2ns的运算(抱歉,图片说明应该是500ps,我错误标成了ns)。

我有以下疑问:
(1)运行时长可以直接叠加吗,如果我要获得10ns的运行时长,运行20次即可?
(2)我自己的MD计划运行100ns,于是我运行了2个50ns的,可是我用cpptraj计算骨架RMSD时的情况如下,在第一个50ns和第二个50ns中间有一个明显的落差,而教程中的骨架RMSD则没有明显的落差。请问是我理解有问题,还是我利用measure_equil_rmsd.ptraj脚本提取有问题呢(脚本内容与教程中几乎一致)?

202108022212362246..png (31.99 KB, 下载次数 Times of downloads: 45)

教程500ns,分4次运行

教程500ns,分4次运行

202108022210529255..png (16.59 KB, 下载次数 Times of downloads: 38)

自己100ns分2次运行

自己100ns分2次运行

206

帖子

0

威望

2251

eV
积分
2457

Level 5 (御坂)

2#
发表于 Post on 2021-8-3 00:43:50 | 只看该作者 Only view this author
Are you sure you selected the right coord file, it looks like for the second part of the simulation you take the initial coordinate intead of the previous coordinate

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
Freshfish + 4 正解

查看全部评分 View all ratings

18

帖子

0

威望

111

eV
积分
129

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2021-8-3 09:25:54 | 只看该作者 Only view this author
rpestana94 发表于 2021-8-3 00:43
Are you sure you selected the right coord file, it looks like for the second part of the simulation  ...

噢,我突然明白了。我第二次运行的确是以初始的坐标参数运行的,万分感谢!
又有一个小问题,如果这样的话,分多次运行的优势是什么呢,直接设置总时长1次运行不就好了吗?

206

帖子

0

威望

2251

eV
积分
2457

Level 5 (御坂)

4#
发表于 Post on 2021-8-3 09:44:25 | 只看该作者 Only view this author
Freshfish 发表于 2021-8-2 20:25
噢,我突然明白了。我第二次运行的确是以初始的坐标参数运行的,万分感谢!
又有一个小问题,如果这样的 ...

You can do both and is fine, as far I know, but depending in what you want, you can use one or another, check this two links, are discussions about that:

https://www.researchgate.net/pos ... -multiple-short-MDs

https://www.researchgate.net/pos ... ltiple_run_dynamics

88

帖子

0

威望

2390

eV
积分
2478

Level 5 (御坂)

5#
发表于 Post on 2021-8-3 11:34:05 | 只看该作者 Only view this author
restart file 可以設定小一點...不過你的RMSD變動真的太大了..要檢查一下結構

1171

帖子

7

威望

6862

eV
积分
8173

Level 6 (一方通行)

6#
发表于 Post on 2021-8-3 11:41:41 | 只看该作者 Only view this author
Freshfish 发表于 2021-8-3 09:25
噢,我突然明白了。我第二次运行的确是以初始的坐标参数运行的,万分感谢!
又有一个小问题,如果这样的 ...

你这种情况多次运行没有什么特别的优势,并且会放大initial bias的影响

18

帖子

0

威望

111

eV
积分
129

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2021-8-4 15:58:21 | 只看该作者 Only view this author
c00jsw00 发表于 2021-8-3 11:34
restart file 可以設定小一點...不過你的RMSD變動真的太大了..要檢查一下結構

你好,进一步请教一下, restart file设定小一点是什么意思?
RMSD变动大是因为这是一个蛋白-蛋白复合体,可能有一些构象变化,我将restrain_wt设定值提高了一些,RMSD变动小了很多。但是收敛时间依然是不太短。我尝试将restrain_wt设为500或者1000,结果会报错:“Coordinate resetting (SHAKE) cannot be accomplished, deviation is too large.”

18

帖子

0

威望

111

eV
积分
129

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2021-8-4 16:00:23 | 只看该作者 Only view this author
fhh2626 发表于 2021-8-3 11:41
你这种情况多次运行没有什么特别的优势,并且会放大initial bias的影响

我根据rpestana94推荐的网站,发现部分情况分多次运行更有好处。根据您的建议,单次运行更好,请教一下,是因为MD的运行对象有区别或者是有其他判断标准吗?
十分感谢。

1171

帖子

7

威望

6862

eV
积分
8173

Level 6 (一方通行)

9#
发表于 Post on 2021-8-8 14:07:41 | 只看该作者 Only view this author
Freshfish 发表于 2021-8-4 16:00
我根据rpestana94推荐的网站,发现部分情况分多次运行更有好处。根据您的建议,单次运行更好,请教一下, ...

具有随机性的非平衡模拟才需要多次运行,平衡模拟单次、长时间运行更好

16

帖子

0

威望

91

eV
积分
107

Level 2 能力者

10#
发表于 Post on 2023-12-28 15:42:54 | 只看该作者 Only view this author
Freshfish 发表于 2021-8-3 09:25
噢,我突然明白了。我第二次运行的确是以初始的坐标参数运行的,万分感谢!
又有一个小问题,如果这样的 ...

请问第二次运行怎么选择以第一次的初始坐标计算RMSD呀?

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 22:50 , Processed in 1.677865 second(s), 24 queries , Gzip On.

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