计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 求助粗粒化的膜蛋白能量Min的问题

求助粗粒化的膜蛋白能量Min的问题

查看数: 10268 | 评论数: 5 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

路过的一只

发布时间: 2019-5-17 18:58

正文摘要:

拜托各位大神看过来……求助问题，我用Martini和insane的脚本做了一个粗粒化的膜蛋白（脂质成分较复杂，已溶剂化），在进行minimization时的结果不太理解……我的mdp文件如图，是教程的文件略改

回复 Reply

路过的一只 发表于 Post on 2019-5-19 09:45:53

sobereva 发表于 2019-5-19 00:24
本来蛋白质、磷脂在水环境下就并不都是中性的，整体带净电荷是司空见惯的，不加抗衡离子显然不行。我建议 ...

谢谢老师回复！好的老师~

sobereva 发表于 Post on 2019-5-19 00:24:22

路过的一只发表于 2019-5-18 15:11
老师，我搜到了这个命令行：$gmx/genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -np 225 ...

本来蛋白质、磷脂在水环境下就并不都是中性的，整体带净电荷是司空见惯的，不加抗衡离子显然不行。我建议你先把普通力场的蛋白质动力学基本流程搞明白再去做粗粒化的。盒子整体电荷必须通过加抗衡离子使之为0，否则化学意义上说不通，程序计算上也会遇到问题

路过的一只 发表于 Post on 2019-5-18 15:11:42

sobereva 发表于 2019-5-17 22:38
盒子里没加抗衡离子？

老师，我搜到了这个命令行：$gmx/genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -np 225 -nname CL 加入了225个NA之后，不会出现那个note了！~但是我还不是很理解为什么会出现这个问题？现在加的na+很明显比cl-多出好多……

路过的一只 发表于 Post on 2019-5-18 09:58:57

sobereva 发表于 2019-5-17 22:38
盒子里没加抗衡离子？

谢谢老师回复！我在用insane生成膜&盒子时有加入0.15的nacl，之后直接进行的这步min，我看到note说non-zero total charge，但我不是很清楚为什么和怎么加这个抗衡离子呢？

sobereva 发表于 Post on 2019-5-17 22:38:16

盒子里没加抗衡离子？

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register