| 查看: 1014 | 回复: 1 | ||
zhangmao511铁虫 (小有名气)
|
[求助]
膜蛋白模拟，能量最小化无法收敛的问题！！
|
|
这2天，在做膜蛋白模拟的时候碰到一个问题，查了半天也没解决。我用的POPE做膜，POPE bilayer with 340 lipids，然后把我的蛋白嵌进去。 在嵌蛋白的过程中，我是用inflategro.pl做的，先体积膨胀4个scaling factor，再能量最小化；再体积收缩0.95个scaling factor，不停的收缩直到area per lipid达到71 Å2。 但是在膨胀后的能量最小化的时候，无法收敛。。 如图1：  图1.jpg 按照上面的指示，9600和9606这2个原子的F max一下子升到了e+13次方。我查了一下这2个原子，是P上的氧原子和另外的H原子。然后，我查了一下，有教程这样说： 如图2：  图2.jpg 我猜测是Add [ exclusions ] 这种情况，所以在我的topol.top文件后面加上了： [ exclusions ] ;atom1 atom2 9606 9600 如图3：  图3.jpg 但是，发现还是能量最小化无法收敛，依然在atom=9600那一步Fmax陡然上升到e+13次方。 不知道，大侠们能不能给一点指点，感谢！！  |
回复此楼» 猜你喜欢
315求调剂
已经有8人回复
-
312求调剂
已经有5人回复
-
一志愿郑州大学材料与化工085600，求调剂
已经有4人回复
-
315分 085602 求调剂
已经有7人回复
-
298分 070300求调剂
已经有5人回复
-
349求调剂
已经有9人回复
-
302求调剂
已经有4人回复
-
化工调剂303分，过四级
已经有10人回复
-
321求调剂
已经有14人回复
-
296材料专硕求调剂
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
zhangmao511
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 33
- 红花: 2
- 帖子: 136
- 在线: 29.7小时
- 虫号: 2553312
- 注册: 2013-07-19
- 性别: GG
- 专业: 蛋白质组学
2楼2014-11-12 09:06:36
