| 查看: 1240 | 回复: 3 | |||
[求助]
急！急！急！用Gromacs时在赋力场这步出现了错误，不知道该怎么解决？ 已有1人参与
|
|
我在用Gromacs时，使用program pdb2gmx准备给酶赋上charmm27的力场，却这出现了错误，错误的信息截图附在后面。 由于我的酶是人工带上一个铜离子，可能Gromacs不能给附加离子的酶赋力场，那么请问在加入离子那步该怎么输入命令，才能把额外的铜离子给添加进去？  错误信息.png  查看UNK位置.png |
回复此楼» 猜你喜欢
能否退出参与的面上项目解除限项
已经有6人回复
-
什么时候开奖？
已经有11人回复
-
范进中举一文的中心思想
已经有6人回复
-
建议基金发布提前给出明确的时间点
已经有14人回复
-
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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
40