计算化学公社

 找回密码 Forget password
 注册 Register

求助:使用Gromacs模拟双链DNA后轨迹处理问题

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

正文摘要:

本帖最后由 HZW 于 2022-5-19 12:23 编辑 我在使用gromacs模拟一个36bp的双链DNA在150mM的NaCl溶液中的运动后,模拟的时间为300ns,由于在MD时未设置对核酸结构进行消除质心平动和转动的参数;对得到的轨迹使用卢 ...

回复 Reply

tzypla 发表于 Post on 2025-2-15 20:44:42
HZW 发表于 2022-5-24 10:57
前天试过了,用-pbc mol处理后的轨迹虽然没有完整的双链结构,但是单条链是完整的,然后用-pbc nojump处 ...

请问你解决了DNA链分开的情况吗?能不能分享一下是怎么解决的
xuhuiting 发表于 Post on 2023-8-6 09:52:28
HZW 发表于 2022-5-24 10:57
前天试过了,用-pbc mol处理后的轨迹虽然没有完整的双链结构,但是单条链是完整的,然后用-pbc nojump处 ...

我模拟DNA遇到一模一样的问题,请问你是怎么解决的?
HZW 发表于 Post on 2022-5-24 10:57:36
sobereva 发表于 2022-5-19 21:54
先-pbc mol搞一次,再-pbc nojump搞一次。如果此时DNA运动已经连续了而且保持完整了,VMD里对DNA相对于第一 ...

前天试过了,用-pbc mol处理后的轨迹虽然没有完整的双链结构,但是单条链是完整的,然后用-pbc nojump处理后的轨迹则直接由于PBC切割而断成了几截了。算了,我去gromacs forum上问下,谢谢卢老师的建议。
HZW 发表于 Post on 2022-5-24 10:52:41
本帖最后由 HZW 于 2022-5-24 12:47 编辑

唉,在MD运行中加了对核酸组的移除质心的平动和转动,在模拟的100ns以内没出现两条链分开的情况,但是100ns以后的轨迹就又是上述情况。
comm-mode               = Angular
comm-grps                 = DNA
看来只能通过增大水盒子的方式去解决问题了
喵星大佬 发表于 Post on 2022-5-21 02:24:44
波波波 发表于 2022-5-20 11:27
我最近发现-center选项没用,选择的组没有被移动到盒子中心,用的是2020.7版本,是不是bug呢?

双链DNA的链质量差不多,如果一个在盒子这边一个在那边质心也在中间,所以要先移动到一起再center
波波波 发表于 Post on 2022-5-20 11:27:43
我最近发现-center选项没用,选择的组没有被移动到盒子中心,用的是2020.7版本,是不是bug呢?
sobereva 发表于 Post on 2022-5-19 22:51:21
HZW 发表于 2022-5-19 22:02
好的,谢谢卢老师,是不是我后面跑双链核酸,在MD参数中使用comm-mode = angular对核酸分子的group消除平 ...

通常是
HZW 发表于 Post on 2022-5-19 22:02:33
Puying 发表于 2022-5-19 18:47
实在不行就用vmd里面的pbc 显示多个就好了~~~

好的,谢谢
HZW 发表于 Post on 2022-5-19 22:02:06
sobereva 发表于 2022-5-19 21:54
先-pbc mol搞一次,再-pbc nojump搞一次。如果此时DNA运动已经连续了而且保持完整了,VMD里对DNA相对于第一 ...

好的,谢谢卢老师,是不是我后面跑双链核酸,在MD参数中使用comm-mode = angular对核酸分子的group消除平动和转动,后续轨迹就不会出现这样呢。
sobereva 发表于 Post on 2022-5-19 21:54:51
先-pbc mol搞一次,再-pbc nojump搞一次。如果此时DNA运动已经连续了而且保持完整了,VMD里对DNA相对于第一帧做个align,就一直在中央了
Puying 发表于 Post on 2022-5-19 18:47:12
HZW 发表于 2022-5-19 18:20
我昨晚就把trjconv -pbc mol ; -pbc atom ;  -pbc whole ; -pbc nojump等所有都试了,没用,依旧不能完全 ...

实在不行就用vmd里面的pbc 显示多个就好了~~~
HZW 发表于 Post on 2022-5-19 18:20:37
Frozen-Penguin 发表于 2022-5-19 17:13
可以试试用trjconv -pbc nojump

我昨晚就把trjconv -pbc mol ; -pbc atom ;  -pbc whole ; -pbc nojump等所有都试了,没用,依旧不能完全消除PBC;然后今天试了论坛搜到的卢老师的回帖,然后又试了网上的VMD处理:pbc wrap -compound res -all,
pbc box 都没用。
对了还有楼上说的用-trans, 这个试错法
Frozen-Penguin 发表于 Post on 2022-5-19 17:13:38
可以试试用trjconv -pbc nojump
Puying 发表于 Post on 2022-5-19 13:19:47
HZW 发表于 2022-5-19 13:00
使用-trans 指定的x y z 的vector怎么确定呢

我没有用过这个命令,不过你可以试试。比如你给的这个图,向上平移是x方向,平移2nm,那就是2 0 0,以此类推
我不确定,你试一下~~

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

GMT+8, 2025-8-16 11:41 , Processed in 0.671125 second(s), 31 queries , Gzip On.

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