计算化学公社

 找回密码 Forget password
 注册 Register

如何正确筛除amber拓扑和轨迹文件中的水分子?

查看数: 1917 | 评论数: 3 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2024-2-12 19:32

正文摘要:

本帖最后由 yorksy 于 2024-2-12 21:12 编辑 各位老师好,我在用amber分析一次md的结果。因为这次模拟的水分子很多,我想用strip去除所有在我的配体7埃外的水分子。我在parmed和cpptraj都用了strip,reference是 ...

回复 Reply

Serious 发表于 Post on 2024-2-15 23:25:10
yorksy 发表于 2024-2-14 17:11
我也试过每一帧单独trajout,但这样就导致我不知道用什么拓扑结构。我的拓扑文件里面把所有水分子的编号 ...

1)如果要导出相应top文件,可以试试cpptraj里parmwrite 和 parmstrip;
2)不过这样子挺麻烦的。。如果只是看结构,导pdb挺好;VMD没用过,不知论坛里的同志有没有办法。。
yorksy 发表于 Post on 2024-2-14 17:11:27
Serious 发表于 2024-2-12 23:57
1)我觉得一帧一帧单独trajout比较好;
2)假定reference是3_prot_npt_0000.rst,strip是依照reference去 ...

我也试过每一帧单独trajout,但这样就导致我不知道用什么拓扑结构。我的拓扑文件里面把所有水分子的编号都写上了。如果不能导出对应的拓扑文件的话,就会导致我的拓扑和结构合不到一起去,VMD里面显示一团乱。
Serious 发表于 Post on 2024-2-12 23:57:27
1)我觉得一帧一帧单独trajout比较好;
2)假定reference是3_prot_npt_0000.rst,strip是依照reference去找保留的分子。比如第100位水分子在reference中是保留的,符合strip距离的设定,那么trajin的所有帧都保留第100位水分子,所以会有“发现在某一个时间点它们跟配体很近,大概三到四埃的样子。但是其他时间离配体都很远”的情况;总的来说,strip的分子是按照reference定的;

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-4-27 13:51 , Processed in 0.559561 second(s), 32 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list