请选择 进入手机版 | 继续访问电脑版
第13届北京科音初级量子化学培训班将于10月5~8日于北京举办,请点击此链接查看详情。这是新人一次性正确、完整学习量子化学计算的最好、最快机会,能少走无数弯路,欢迎参加并相互转告!(已报满)

计算化学公社

 找回密码
 现在注册!
查看: 118|回复: 13

[GROMACS] 金属蛋白生成top文件时,CA离子被当作了C的问题求助

[复制链接]

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

发表于 6 天前 | 显示全部楼层 |阅读模式
请问在使用amberFF99sb力场生成蛋白top文件,分别生成了蛋白top和金属top,但是金属中的CA离子被当作了C处理,导致我后面进行添加溶剂时会出现warning,对于这种情况,我该如何处理才对?求助
M(MAG428UUYC$(LHMJQ~R2Y.png

topol_Ion_chain_A2.itp

1.21 KB, 下载次数: 0

120

帖子

0

威望

365

eV
积分
485

Level 3 能力者

发表于 5 天前 | 显示全部楼层
1 这个金属离子是抗衡离子还是配位离子,抗衡的话(就是不确定坐标,在哪里都可以)在制作top时用命令自动添加
2 如果你是配位的,必须在某一特定位置
    1)检查你的pdb文件里的CA离子那行的最后一列是不是C和上下对齐而不是A,你缩进一个空格就行(这个如果可以你能不能发一下pdb的文件,CA这里截个图)
    2)实在搞不定就用第一步加一个随机位置的抗衡CA离子然后把原来配位的坐标复制过来(下下策)
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

 楼主| 发表于 5 天前 | 显示全部楼层
DoubeeTwT 发表于 2020-9-15 08:46
1 这个金属离子是抗衡离子还是配位离子,抗衡的话(就是不确定坐标,在哪里都可以)在制作top时用命令自动 ...

你好,我这个是配位离子,我上传了我的pdb文件,你看看,我在生产top时,他就变成c。

targetx.pdb

110.59 KB, 下载次数: 2

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

 楼主| 发表于 5 天前 | 显示全部楼层
DoubeeTwT 发表于 2020-9-15 08:46
1 这个金属离子是抗衡离子还是配位离子,抗衡的话(就是不确定坐标,在哪里都可以)在制作top时用命令自动 ...

这是pdb的截图,后面是四个金属离子

3N%Z%V01XRVTIX43}O8K@CJ.png

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

 楼主| 发表于 5 天前 | 显示全部楼层
DoubeeTwT 发表于 2020-9-15 08:46
1 这个金属离子是抗衡离子还是配位离子,抗衡的话(就是不确定坐标,在哪里都可以)在制作top时用命令自动 ...

原子里面有种写法C0就代表CA,是没错的,是我自己理解错了,谢谢你

120

帖子

0

威望

365

eV
积分
485

Level 3 能力者

发表于 5 天前 | 显示全部楼层
本帖最后由 DoubeeTwT 于 2020-9-15 14:45 编辑
琦锅锅 发表于 2020-9-15 10:46
原子里面有种写法C0就代表CA,是没错的,是我自己理解错了,谢谢你

所以现在解决了吗?看pdb似乎没什么问题
试试第四列(第二个CA往前缩进一格)

上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

2万

帖子

25

威望

3万

eV
积分
58547

管理员

公社社长

发表于 5 天前 | 显示全部楼层
不要24小时内同时在公社QQ群和公社论坛里同一问题,注意看群文件里群规,以及论坛置顶的计算化学公社论坛新人必读。
我都已经在群里回复过了。下次再这样就记违规了
Clipboard01.png
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

 楼主| 发表于 5 天前 | 显示全部楼层
sobereva 发表于 2020-9-15 14:54
不要24小时内同时在公社QQ群和公社论坛里同一问题,注意看群文件里群规,以及论坛置顶的计算化学公社论坛新 ...

好的,以后会注意的,谢谢sob老师

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

 楼主| 发表于 5 天前 | 显示全部楼层
DoubeeTwT 发表于 2020-9-15 14:38
所以现在解决了吗?看pdb似乎没什么问题
试试第四列(第二个CA往前缩进一格)

嗯嗯,解决了。我想请教一下,我的金属在活性位点处,为防止金属在MD过程中跑飞,要进行距离限制,一般是不是要在top文件中加入[distance_constraints]在最后面么?请问一般是什么格式

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

 楼主| 发表于 5 天前 | 显示全部楼层
DoubeeTwT 发表于 2020-9-15 14:38
所以现在解决了吗?看pdb似乎没什么问题
试试第四列(第二个CA往前缩进一格)

这是我的top文件,我在网上看了挺多关于金属限制的,但是,对于它的处理还是很懵懂,所以,想请教一下你,该如何建立距离限制它,非常感谢

topol.top

680.17 KB, 下载次数: 1

2万

帖子

25

威望

3万

eV
积分
58547

管理员

公社社长

发表于 4 天前 | 显示全部楼层
琦锅锅 发表于 2020-9-15 20:35
嗯嗯,解决了。我想请教一下,我的金属在活性位点处,为防止金属在MD过程中跑飞,要进行距离限制,一般是 ...

格式在手册里的拓扑文件部分都有说明
先尝试直接跑,如果呆不住,再考虑加位置限制
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

 楼主| 发表于 4 天前 | 显示全部楼层
sobereva 发表于 2020-9-16 09:20
格式在手册里的拓扑文件部分都有说明
先尝试直接跑,如果呆不住,再考虑加位置限制

好的,谢谢sob老师

120

帖子

0

威望

365

eV
积分
485

Level 3 能力者

发表于 前天 10:21 | 显示全部楼层
琦锅锅 发表于 2020-9-16 14:59
好的,谢谢sob老师

对于限制性MD,其实不是自然的情况,一般有两种使用到的情形
1)拉伸动力学:比如你知道一段loop区存在开合两种构象, 你像研究两种构象之间的转化时长,受力,能量变化。这时候外加限制去移动它
2)构建体系:有时候复合物体系没有可信的晶体结构,同源建模的时候也是空载蛋白为模板,但我们通过实验上的突变已经可以确定口袋和活性位点,此时对接后通过限制性MD在对接结构下让受体和配体稳定一下有利于之后放开限制后的模拟体系更加稳定

因此对照上述,考察你的体系是否需要限制,毕竟是外力,加入要慎重,即使情况2也只是通过限制合理化初始结构,最后还是需要在无约束条件下模拟的
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

30

帖子

0

威望

43

eV
积分
73

Level 2 能力者

 楼主| 发表于 前天 21:15 | 显示全部楼层
DoubeeTwT 发表于 2020-9-18 10:21
对于限制性MD,其实不是自然的情况,一般有两种使用到的情形
1)拉伸动力学:比如你知道一段loop区存在开 ...

好的,我在研究研究,非常感谢你的解答
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-9-20 22:22 , Processed in 0.177597 second(s), 27 queries .

快速回复 返回顶部 返回列表