计算化学公社

标题: Gromacs可以进行含有铜离子的蛋白质的分子动力学模拟吗? [打印本页]

作者
Author:
五支折断的箭    时间: 2020-5-2 18:36
标题: Gromacs可以进行含有铜离子的蛋白质的分子动力学模拟吗?
各位前辈好,我想请教一下Gromacs可以进行含有铜离子的蛋白质的分子动力学模拟吗?我在进行模拟时,在使用pdb2gmx指令后,会出现:如下错误
Fatal error:
Residue 'CU ‘ not found in residue topology database
有什么办法可以解决吗?是不是需要在atomtypes.atp中添加参数呢?请前辈指点。


作者
Author:
苏玖染    时间: 2020-5-2 20:49
金属蛋白需要考虑尝试Amber的做法,Amber有专门的教程
作者
Author:
五支折断的箭    时间: 2020-5-2 22:59
苏玖染 发表于 2020-5-2 20:49
金属蛋白需要考虑尝试Amber的做法,Amber有专门的教程

好的,谢谢了

作者
Author:
五支折断的箭    时间: 2020-5-2 23:16
苏玖染 发表于 2020-5-2 20:49
金属蛋白需要考虑尝试Amber的做法,Amber有专门的教程

你好,我还想问一下,我已经做了金属蛋白和有机小分子的对接了,结合的位置并不在铜的旁边,这种情况下如果把铜去掉做分子动力学模拟是否可行呢?感谢。
作者
Author:
sobereva    时间: 2020-5-3 08:18
gromacs可以很好地对含有金属的蛋白体系进行模拟。没有必要非得换成amber

当前报错是因为你用的力场目录下的rtp文件里没有叫CU的残基,让pdb里的残基名与rtp里对应上就能解决。诸如G54A7目录下的aminoacids.rtp里铜有CU+和CU2+两种,你需要根据实际情况将pdb里铜的残基名设为二者之一。
作者
Author:
znjnancy    时间: 2023-12-27 21:14
sobereva 发表于 2020-5-3 08:18
gromacs可以很好地对含有金属的蛋白体系进行模拟。没有必要非得换成amber

当前报错是因为你用的力场目录 ...

sob老师,我是非标准残基小肽和铜离子分子动力学,我已经构建好了amber99sb-ildn力场的非标准残基rtp文件,如果改成g54a7力场,虽然认识铜了,但是非标准残基又要报错了,请问这种情况如何选择力场呢?
作者
Author:
sobereva    时间: 2023-12-28 10:01
znjnancy 发表于 2023-12-27 21:14
sob老师,我是非标准残基小肽和铜离子分子动力学,我已经构建好了amber99sb-ildn力场的非标准残基rtp文件 ...

用G54A7做什么
Cu用Merz的离子参数,跟AMBER完全兼容
作者
Author:
HigherH    时间: 2024-4-19 10:41
苏玖染 发表于 2020-5-2 20:49
金属蛋白需要考虑尝试Amber的做法,Amber有专门的教程


请问能分享一下金属蛋白的amber教程吗,感谢




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3