24小时热门版块排行榜    

查看: 787  |  回复: 3

maggie凤

木虫 (正式写手)

[求助] Gromacs把蛋白的二级结构跑散了怎么办 已有1人参与

RT
用的1UBI蛋白 按照下面的步骤跑了10ns模拟
pdb2gmx -f 1UBI.pdb -o 1UBI.gro
editconf -f 1UBI.gro -o conf.gro -c -box 5 5 5
genbox -cp conf.gro -cs -p topol.top -o box.gro

grompp -c box.gro -p topol.top -f em.mdp -o em.tpr -maxwarn 10
mdrun -deffnm em

grompp -f nvt.mdp -c em.gro -p topol.top -o nvt.tpr -maxwarn 10
mdrun -deffnm nvt

grompp -f npt.mdp -t nvt.cpt -c nvt.gro -p topol.top -o npt.tpr -maxwarn 10
mdrun -deffnm npt

grompp -f md.mdp -c npt.gro -p topol.top -o md.tpr -maxwarn 10
mdrun -deffnm md

最后把.xtc文件做成pdb 在pymol里显示的蛋白里的四个beta折叠全都少了一大截 有的时候还会彻底消失
那些消失的beta折叠的位置的RMSF都偏大 但是不如loop区那么大
求解这是怎么回事 万谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maggie凤

木虫 (正式写手)

试了好几个力场GROMOS OPLSAA AMBER都用过了 效果有一些区别 但是也全都跑散了
这蛋白本身是中性的也不用加离子平衡啊
到底是怎么回事 求大神救急
2楼2014-12-11 13:20:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maggie凤

木虫 (正式写手)

求帮忙啊嘤嘤婴
3楼2014-12-12 18:34:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awaken2013

铜虫 (小有名气)

【答案】应助回帖

就算蛋白是中性的,也要看处理日志里面的提示,有无带电。
跑散了也有可能是本身蛋白就不够稳定。或则你可以使用束缚性条件跑跑看
4楼2014-12-16 11:51:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 maggie凤 的主题更新
信息提示
请填处理意见