| 查看: 1080 | 回复: 3 | ||
[求助]
急！急！急！用Gromacs时在赋力场这步出现了错误，不知道该怎么解决？ 已有1人参与
|
|
我在用Gromacs时，使用program pdb2gmx准备给酶赋上charmm27的力场，却这出现了错误，错误的信息截图附在后面。 由于我的酶是人工带上一个铜离子，可能Gromacs不能给附加离子的酶赋力场，那么请问在加入离子那步该怎么输入命令，才能把额外的铜离子给添加进去？  错误信息.png  查看UNK位置.png |
回复此楼» 猜你喜欢
286求调剂
已经有5人回复
-
085601求调剂
已经有3人回复
-
302求调剂
已经有5人回复
-
考研化学学硕调剂，一志愿985
已经有5人回复
-
328求调剂，英语六级551，有科研经历
已经有4人回复
-
机械专硕325，寻找调剂院校
已经有5人回复
-
材料专硕306英一数二
已经有6人回复
-
伙伴们，祝我生日快乐吧
已经有26人回复
-
333求调剂
已经有7人回复
-
0854控制工程 359求调剂 可跨专业
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用ANSYS--DYNA做锚杆
已经有2人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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
