| 查看: 1105 | 回复: 3 | ||
[求助]
急！急！急！用Gromacs时在赋力场这步出现了错误，不知道该怎么解决？ 已有1人参与
|
|
我在用Gromacs时，使用program pdb2gmx准备给酶赋上charmm27的力场，却这出现了错误，错误的信息截图附在后面。 由于我的酶是人工带上一个铜离子，可能Gromacs不能给附加离子的酶赋力场，那么请问在加入离子那步该怎么输入命令，才能把额外的铜离子给添加进去？  错误信息.png  查看UNK位置.png |
回复此楼» 猜你喜欢
281求调剂
已经有6人回复
-
华南理工0703化学，总分336求调剂
已经有6人回复
-
285求调剂
已经有12人回复
-
085600材料与化工301分求调剂院校
已经有13人回复
-
295求调剂
已经有11人回复
-
材料与化工专硕306分找合适调剂
已经有3人回复
-
0703调剂，一志愿天津大学319分
已经有11人回复
-
求调剂
已经有21人回复
-
材料专硕(0856) 339分求调剂
已经有10人回复
-
319求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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
