计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2066|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[Amber] 求助Amber如何消除蛋白的平动和转动

[复制链接 Copy URL]

13

帖子

0

威望

185

eV
积分
198

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
大佬们好,纯新人,Amber22进行蛋白和小分子配体之间的NVT模拟时发现蛋白会跑出水盒,请教应该如何改善这个问题,最好是能详细一些,感谢大家

附上正式模拟的文件




05_md.in

279 Bytes, 下载次数 Times of downloads: 4

4

帖子

0

威望

207

eV
积分
211

Level 3 能力者

6#
发表于 Post on 前天 11:16 | 只看该作者 Only view this author
casea 发表于 2022-12-3 09:04
cpptraj 导入轨迹之后 autoimage

请问autoimage怎么操作呢

13

帖子

0

威望

185

eV
积分
198

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2022-12-3 12:45:28 | 只看该作者 Only view this author
casea 发表于 2022-12-3 09:04
cpptraj 导入轨迹之后 autoimage

十分感谢您的解答,这样确实有所改善

67

帖子

3

威望

1545

eV
积分
1672

Level 5 (御坂)

4#
发表于 Post on 2022-12-3 09:04:07 | 只看该作者 Only view this author
亚切ss 发表于 2022-12-3 04:07
十分感谢您的耐心解答,可能是我的描述不够准确,我现在遇到的问题是模拟中蛋白和配体会有跨越周期性边界 ...

cpptraj 导入轨迹之后 autoimage

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
亚切ss + 4 谢谢

查看全部评分 View all ratings

13

帖子

0

威望

185

eV
积分
198

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-12-3 04:07:55 | 只看该作者 Only view this author
Frozen-Penguin 发表于 2022-12-2 20:15
准备好轨迹文件(下面记为trj)和参数文件(下面记为prmtop),用cpptraj读取

然后使用align消除平动转 ...

十分感谢您的耐心解答,可能是我的描述不够准确,我现在遇到的问题是模拟中蛋白和配体会有跨越周期性边界的问题,因此我想在模拟中消除蛋白的整体运动,请问这种情况应该怎么解决呢?

310

帖子

0

威望

1493

eV
积分
1803

Level 5 (御坂)

2#
发表于 Post on 2022-12-2 20:15:13 | 只看该作者 Only view this author
本帖最后由 Frozen-Penguin 于 2022-12-2 20:16 编辑

准备好轨迹文件(下面记为trj)和参数文件(下面记为prmtop),用cpptraj读取
  1. cppraj -p prmtop -y trj
复制代码

然后使用align消除平动转动(如果有跨越周期性边界的情况要在消除平动转动之前处理)
  1. align :1-200@CA
复制代码

这样操作就可以消除1-200号残基的C-alpha原子的平动转动
最后把结果保存
  1. trajout traj.mdcrd
复制代码

cpptraj中的指令输入之后再输入go或者run才会执行
  1. go
复制代码


更多用法可以用help+命令的方式查看,例如
  1. help align
复制代码


评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
亚切ss + 4 谢谢

查看全部评分 View all ratings

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

GMT+8, 2025-8-18 04:23 , Processed in 0.180652 second(s), 25 queries , Gzip On.

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