| 查看: 1276 | 回复: 3 | ||
[求助]
急！急！急！用Gromacs时在赋力场这步出现了错误，不知道该怎么解决？ 已有1人参与
|
|
我在用Gromacs时，使用program pdb2gmx准备给酶赋上charmm27的力场，却这出现了错误，错误的信息截图附在后面。 由于我的酶是人工带上一个铜离子，可能Gromacs不能给附加离子的酶赋力场，那么请问在加入离子那步该怎么输入命令，才能把额外的铜离子给添加进去？  错误信息.png  查看UNK位置.png |
回复此楼
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
屁桃粉: 金币+20, ★★★很有帮助 2015-11-04 10:35:02
感谢参与，应助指数 +1
屁桃粉: 金币+20, ★★★很有帮助 2015-11-04 10:35:02
|
你说得没错，是铜离子的缘故！因为缺少铜离子拓扑文件，所以报错！至于怎么加，就是要写一个铜离子的拓扑文件。你找到铜离子的相关参数，写一个铜离子的itp文件，然后引入到蛋白top文件中，当然你也可以直接在蛋白top文件中写入。参数你可以选用gaff立场参数。你也可以用量力化学计算得到，或者找文献。至于怎么写top文件，最简单的方法是，你按照gromacs教程下的protein-ligand systems 下的第二个教程，enzyme complex solvation 的方法。虽然他是钙离子，你是铜离子，你先按照钙离子的方法生存top文件，然后将钙离子的参数替换为铜就可以。或者你把铜离子看成一个配体，按照一个小分子的方法处理它即可。 发自小木虫IOS客户端 |
2楼2015-11-03 21:17:07
3楼2015-11-04 10:36:12
4楼2015-11-04 16:03:47