计算化学公社

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

[Amber] 求助,在Amber模拟中控制两体系水盒子数目相同

[复制链接 Copy URL]

138

帖子

0

威望

445

eV
积分
583

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
背景介绍:
目前在分子对接后想用Amber做一下能量最小化来粗略比较对接结果哪个更好
但比较结果后发现能量的差异大多来自水盒子中水分子数目的不同导致

因为对接的小分子过于巨大(> 170 atom)所以在autodock中无法使用全柔性对接,半柔性对接的结果我又怕是局部最优解想做一下想做一下minimization

目前问题:
(1)将tleap之后的结果保存为pdb
(2)按照(比如我要10000个水分子,实际TIP3PBOX中有10400个,就需要删除400个)多出的部分,在最尾端删除对应数目的水分子数
(3)重新tleap,此时不再添加水盒子(不使用solvateoct)可以获得top和crd输出
(4)进行amber能量最小化计算,报错:
        |  ERROR:   Box parameters not found in inpcrd file!

PS:对于对接结果的进一步的优化和判断是否有更好的方法?
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

1

帖子

0

威望

64

eV
积分
65

Level 2 能力者

8#
发表于 Post on 2022-12-19 08:56:16 | 只看该作者 Only view this author
洁洁lili 发表于 2021-11-8 10:22
**** 作者被禁止或删除 内容自动屏蔽 ****

可以用cpptraj的parmstrip删除Amber topology file里面的residue number的xxx的水,选择residue的语法和Amber input file里面的一样:
parm protein_solv.prmtop
parmstrip :xxx
parmwrite out protein_solv.truncated.parm7
再用strip删除coordinate file里面对应的水:
parm protein_solv.prmtop
trajin protein_solv.inpcrd
strip :xxx,xxx,xxx-xxx
trajout protein_solv.truncated.ncrst
最后用protein_solv.truncated.parm7和protein_solv.truncated.ncrst作为输入文件跑MD就行了。
pdb文件格式里没法存储周期性溶剂盒子的相关信息,所以就算可以手动删除pdb里面的水,pdb文件也没法用作Amber输入文件。

4

帖子

0

威望

29

eV
积分
33

Level 2 能力者

7#
发表于 Post on 2021-11-8 10:22:03 | 只看该作者 Only view this author
sobereva 发表于 2021-11-4 22:16
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620(http: ...

老师不好意思。描述的过于草率。
我再描述一下:
最近我在准备用amber跑靶向分子动力学,但是要求始态与终态的分子量一致
但是amber是直接用以半径加水盒子的,所以准备好的文件显示水分子的数量是不一致的
请问有没有办法可以加入同等数量的水分子呢?
(我导出了pdb文件,直接删除了一定数量的水分子,然后重新准备不加水盒子,但是会报error)
有网友给我推荐Solvate.sh这个处理,我准备去试试,非常感谢老师的回复!

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

6#
发表于 Post on 2021-11-4 22:16:34 | 只看该作者 Only view this author
洁洁lili 发表于 2021-11-4 17:54
老师,请问您知道如何控制两个体系的水分子量一样嘛?

在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620http://bbs.keinsci.com/thread-25787-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

4

帖子

0

威望

29

eV
积分
33

Level 2 能力者

5#
发表于 Post on 2021-11-4 17:54:05 | 只看该作者 Only view this author
sobereva 发表于 2020-7-28 06:00
靠能量极小化判断这个意义不大。光从能量极小化之后的能量来对比,还不如多数对接程序给的打分函数判断来得 ...

老师,请问您知道如何控制两个体系的水分子量一样嘛?

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

4#
发表于 Post on 2020-7-31 01:28:27 | 只看该作者 Only view this author
DoubeeTwT 发表于 2020-7-30 10:35
好的,谢谢Sob老师原本的实验计划是用Minimization粗略筛选一遍再用MMPDBSA计算
那现在就不需要初筛了, ...

没必要分别做几次取平均,直接跑一条较长的轨迹取平均就够了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

138

帖子

0

威望

445

eV
积分
583

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2020-7-30 10:35:56 | 只看该作者 Only view this author
本帖最后由 DoubeeTwT 于 2020-7-30 10:39 编辑
sobereva 发表于 2020-7-28 06:00
靠能量极小化判断这个意义不大。光从能量极小化之后的能量来对比,还不如多数对接程序给的打分函数判断来得 ...

好的,谢谢Sob老师原本的实验计划是用Minimization粗略筛选一遍再用MMPDBSA计算
那现在就不需要初筛了,直接上MMPBSA

顺便再问一下老师:
之前我做MMPBSA时发现,并列实验时计算得到的结合自由能会有上下 10Kcal/mol左右的波动
那我一次用多长时间的模拟做MMPBSA,并行几次取均值比较可信?(目前是 50ns x 3)
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

2#
发表于 Post on 2020-7-28 06:00:08 | 只看该作者 Only view this author
靠能量极小化判断这个意义不大。光从能量极小化之后的能量来对比,还不如多数对接程序给的打分函数判断来得靠谱。
应该跑MD,然后做MMPBSA
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

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

GMT+8, 2026-2-25 00:52 , Processed in 0.180387 second(s), 27 queries , Gzip On.

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