24小时热门版块排行榜    

查看: 4033  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 超人与小木虫 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 东南大学364求调剂 +3 JasonYuiui 2026-03-15 3/150 2026-03-15 18:57 by 无际的草原
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
请填处理意见