24小时热门版块排行榜    

查看: 3935  |  回复: 10

超人与小木虫

金虫 (小有名气)

[求助] Amber 分子动力学 求助!!!!

大家好!我用amber软件跑了100ns的动力学模拟,采用周期性边界条件,显性水溶剂模型,之前没有在输入文件中加入iwrap=1这个参数,跑到50ns的时候自动断掉!此时最后一帧构象中有的原子坐标为“****”,导致下一步不能读取.rst文件继续往下计算,(PS:原因是某些原子的坐标跑出了水盒子??)

加了iwrap=1这个命令后,可以跑完100ns,但是提取最后一帧构象发现,原来的蛋白构象(二聚体)跑散了!!!其中的配体也远离蛋白!!请问一下,这是怎么回事??这个问题该如何解决??急死人啊!!跑了很多天才跑完,无语了。。。。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

超人不会飞
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsingtaobeer

新虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
首先100ns你觉得有必要么?其次你的蛋白是两段吗?再则小分子跑出去也属于正常啊,结合力不牢

[ 发自手机版 http://muchong.com/3g ]
2楼2013-10-15 01:24:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超人与小木虫

金虫 (小有名气)

引用回帖:
2楼: Originally posted by tsingtaobeer at 2013-10-15 01:24:03
首先100ns你觉得有必要么?其次你的蛋白是两段吗?再则小分子跑出去也属于正常啊,结合力不牢

首先非常感谢你的回复!对于我这个体系,由于是探索二聚体构象变化,100ns已经算是比较少的了。蛋白是二聚体,就算分离也不会变得那么不靠谱。小分子是ATP,与其中一个单体结合非常牢固。是不是在加了iwrap=1这个参数后,需要对轨迹做什么特殊处理?恢复镜像??
超人不会飞
3楼2013-10-15 09:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaozhq

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
周期性边界也会跑出盒子? lz应该把 那50ns的轨迹好好看看再下定论
你做的是高温下的动力学么? 从来没加过iwrap=1这个选项
分开的原因还是很有可能电荷排斥吧
4楼2013-10-15 10:13:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超人与小木虫

金虫 (小有名气)

引用回帖:
4楼: Originally posted by yaozhq at 2013-10-15 10:13:00
周期性边界也会跑出盒子? lz应该把 那50ns的轨迹好好看看再下定论
你做的是高温下的动力学么? 从来没加过iwrap=1这个选项
分开的原因还是很有可能电荷排斥吧

的确采用的周期性边界条件,50ns后的rst文件中出现了*****,查了其他人遇到的问题,有人回复说是某些原子坐标超出盒子范围了。我模拟温度是343K,不加iwrap=1这个参数,二聚体不会分开。加了之后,可以继续跑,但我疯了。。。。。
超人不会飞
5楼2013-10-15 11:35:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lccbat

金虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
iwrap=1应该不会对模拟结果有很大影响。需要确定的是在提取构象的时候,有没有把分子先影射到原始的盒子中。考虑到模拟的温度很高,二聚体分开也算正常。
6楼2013-10-16 02:46:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超人与小木虫

金虫 (小有名气)

引用回帖:
6楼: Originally posted by lccbat at 2013-10-16 02:46:28
iwrap=1应该不会对模拟结果有很大影响。需要确定的是在提取构象的时候,有没有把分子先影射到原始的盒子中。考虑到模拟的温度很高,二聚体分开也算正常。

好的,非常感谢!!
超人不会飞
7楼2013-10-16 09:53:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1054950560

金虫 (正式写手)

楼主跑二聚体和小分子结合就要跑100ns,那我跑蛋白质蛋白质结合的那不是得好几百ns,我去,疯了
努力就会有回报,不管是好的还是坏的
8楼2013-10-21 15:59:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超人与小木虫

金虫 (小有名气)

引用回帖:
8楼: Originally posted by 1054950560 at 2013-10-21 15:59:09
楼主跑二聚体和小分子结合就要跑100ns,那我跑蛋白质蛋白质结合的那不是得好几百ns,我去,疯了

呵呵,这个不一定!我的小分子是ATP,主要想观察ATP的水解对蛋白构象的影响,所以时间必须长啊!如果你的研究目的不是这个,应该就不需要这么长时间的了。。。。根据自己的研究目的,加上参考文献来定吧。
超人不会飞
9楼2013-10-21 16:46:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1054950560

金虫 (正式写手)

引用回帖:
9楼: Originally posted by 超人与小木虫 at 2013-10-21 16:46:11
呵呵,这个不一定!我的小分子是ATP,主要想观察ATP的水解对蛋白构象的影响,所以时间必须长啊!如果你的研究目的不是这个,应该就不需要这么长时间的了。。。。根据自己的研究目的,加上参考文献来定吧。...

关键的是我这个没有文献报道,是头一次研究,哎命苦啊
努力就会有回报,不管是好的还是坏的
10楼2013-10-21 19:52:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 超人与小木虫 的主题更新
信息提示
请填处理意见