24小时热门版块排行榜    

查看: 1202  |  回复: 1

41XY

铁虫 (初入文坛)

[求助] namd模拟中重原子的固定

新人求助:
初学分子模拟,老师教namd,给的任务是跑完2YHC这个蛋白,可是这个蛋白晶体结构有残缺,从102-115残基都缺,老师又让先学modeller补全,补好之后要用namd对补出的那段(102-115)进行优化,分别有两个配置文件NVT.conf 和NPT.conf,我截取了101-116进行优化,先固定重原子(occupancy列重原子设置成1,氢原子设置0),  依次跑了NVT和NPT,之后放开重原子(配置文件里的一个选项设置成off),又依次跑了NPT和NVT,四步过程中101和116的所有原子(包括氢原子)的occupancy列我都设置是1(我想的是这样模拟完成之后101和116的坐标都没变很好回拼),这样在最后的一步NVT模拟后我将得到的dcd文件转化成pdb文件希望能拼回去可是发现101和116残基的坐标已经变了,老师说要经过平移旋转等操作接回去,现在我技术还没掌握,不过我却对模拟操作有疑问。
模拟操作第三步和第四步中要放开重原子,然后要把conf文件中原子约束行设置成off,这样我就纠结了,那我把101和116的所有原子都固定有什么意义呢,原子约束行设置成off后就不识别pdb文件中occupancy列的数字了吧,namd不会理解我想固定101和116的想法,这样得到的最终结果中101和116的坐标肯定会变。然后老师说可以通过平移旋转等接回去,我看了下101和116的坐标,他们和原坐标不具备简单的统一平移关系,如果硬要把坐标变回去可能会破坏优化的结果,这样优化就没有意义了,所以我怀疑优化过程中约束重原子部分的操作有问题,另一方面就是可能关于回拼的技术平移旋转部分有深意我理解不够,希望得到大家指点。
本人namd新手,老师给个方向自己摸索,理解不对之处多多,前段时间也发过类似的帖子这次是想印证下自己想法,说不清楚所以写的略多。
本人新虫,金币不多,希望有想法的前辈不嫌弃给个回复。谢谢啦。
回复此楼

» 猜你喜欢

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

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

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

fhh2626

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
41XY: 金币+10, ★有帮助, 谢谢 2013-09-10 12:41:42
一般先固定重原子是把水跑平衡,让体系整体能量降低,然后在放开让体系自由运动
2楼2013-09-09 11:15:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 41XY 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 czAiolzR0IJC 2026-11-10 4/200 2026-11-12 02:12 by 9fOTqj074tIg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 4/200 2026-11-12 02:12 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 4/200 2026-11-12 02:12 by 9fOTqj074tIg
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 01:52 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 01:39 by 9fOTqj074tIg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 01:28 by 9fOTqj074tIg
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-12 01:06 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 01:03 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 00:55 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 cqv1zKnzd7hr 2026-11-09 10/500 2026-11-12 00:41 by 9fOTqj074tIg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 cqv1zKnzd7hr 2026-11-09 10/500 2026-11-12 00:35 by 9fOTqj074tIg
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 00:33 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 zSbvVX9JFl9J 2026-11-09 7/350 2026-11-12 00:33 by 9fOTqj074tIg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-11 17:58 by TgvkbZZl9tqB
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 C6CRn66RYYXS 2026-11-10 4/200 2026-11-11 17:53 by TgvkbZZl9tqB
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-11 11:41 by gRND4NVAZ8VG
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
信息提示
请填处理意见