24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4068  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 超人与小木虫 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 309求调剂 +7 谁不是少年 2026-03-29 7/350 2026-03-29 13:17 by mumin1990
[考研] 297求调剂 +11 田洪有 2026-03-26 11/550 2026-03-29 13:14 by peike
[考研] 南京大学化学调剂 +5 景随风 2026-03-29 8/400 2026-03-29 13:00 by 无际的草原
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
[考研] 一志愿华东师范大学有机化学专业,初试351分,复试被刷求调剂! +3 真名有冰 2026-03-29 4/200 2026-03-29 08:47 by qingfeng258
[考研] 0856求调剂 +13 zhn03 2026-03-25 14/700 2026-03-29 08:13 by fmesaito
[考研] 压国家一区线,求导师收留,有恩必谢! +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
[考研] 材料与化工272求调剂 +9 阿斯蒂芬2004 2026-03-28 9/450 2026-03-28 15:21 by 1018329917
[考研] 调剂 +3 好好读书。 2026-03-28 3/150 2026-03-28 12:04 by 王保杰33
[考研] 295求调剂 +5 1428151015 2026-03-27 6/300 2026-03-28 04:04 by fmesaito
[考研] 265求调剂 +8 小木虫085600 2026-03-27 8/400 2026-03-27 22:16 by 无际的草原
[有机交流] 高温高压反应求助 10+4 chibby 2026-03-25 4/200 2026-03-27 21:08 by BT20230424
[考研] 07化学280分求调剂 +10 722865 2026-03-23 10/500 2026-03-27 15:51 by Plutoqq
[考研] 287求调剂 +10 land xuxu 2026-03-26 10/500 2026-03-27 15:33 by 帕尔马拉特
[考研] 安徽大学专硕生物与医药专业(086000)324分,英语已过四六级,六级521,求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 0703化学求调剂 +3 丹青奶盖 2026-03-26 5/250 2026-03-26 20:11 by macy2011
[考研] 271求调剂 +6 生如夏花… 2026-03-22 6/300 2026-03-26 16:48 by 张凯十八号
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +4 我爱学电池 2026-03-23 4/200 2026-03-25 00:59 by 1027_324
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
信息提示
请填处理意见