计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 求助：模拟双链短RNA，跑完md后周期性边界条件始终处理 ...

求助：模拟双链短RNA，跑完md后周期性边界条件始终处理不好

查看数: 1085 | 评论数: 8 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

Roosa

发布时间: 2023-9-15 14:21

正文摘要:

我对一个序列为5'-G6-3'和3'-CU-5'配对的RNA（6表示R6A，N-6甲基修饰的腺嘌呤），使用gromacs带R6A修饰的力场进行500ns的模拟，盒子大小设置的1nm。跑完后使用-pbc mol -center处理轨迹，得到的RNA跑飞了，便-pbc ...

回复 Reply

Christ 发表于 Post on 2023-12-2 05:24:23

我用Amber跑一个dimer的时候也出现这种现象，目前我不知道怎么解决

beyond 发表于 Post on 2023-9-18 21:22:56

做PBC correction，在用center的时候，最好选择的原子在一条链上的，如果选择2条链，很容易两条莲分离

看你的分子也不大，应该比较好调，你可以试一下不用center，或者center的group选择一条链上中间的1-2个氨基酸
-pbc mol 应该就可以了，要实在还有1-2帧弄不好，就vmd手动调一下也可以的。。。

sobereva 发表于 Post on 2023-9-18 10:38:58

Roosa 发表于 2023-9-18 10:33
我后来处理轨迹时先用了-pbc cluster，导出pdb后查看发现大多数两条链完整，偶尔有几帧一条链跑飞了，rms ...

-pbc atom作甚，只会导致分子不完整
-pbc cluster只要用对了（诸如不是一开始就是分家的），不可能两条链在某些帧分家。只要两条链完整挨在一起，算RMSD什么毛病都不可能有

Roosa 发表于 Post on 2023-9-18 10:33:51

sobereva 发表于 2023-9-18 10:29
先用-pbc cluster，把两条链的二聚体整体保持完整再说

我后来处理轨迹时先用了-pbc cluster，导出pdb后查看发现大多数两条链完整，偶尔有几帧一条链跑飞了，rmsd是帖子主楼那张图，我再-pbc atom -center，得到的结果和只-pbc cluster差不多，rmsd高的都往1.5nm走

sobereva 发表于 Post on 2023-9-18 10:29:05

先用-pbc cluster，把两条链的二聚体整体保持完整再说

Roosa 发表于 Post on 2023-9-18 10:23:26

sobereva 发表于 2023-9-16 01:33
缺乏经验的话老老实实用矩形盒子，省得实际计算考虑的周期性和显示的结构的关系弄不清楚

我用的就是矩形的盒子，有关命令是gmx editconf -f name.gro -o name_newbox.gro -c -d 1.0 -bt cubic，还是说您指的矩形是triclinic这种？

sobereva 发表于 Post on 2023-9-16 01:33:10

缺乏经验的话老老实实用矩形盒子，省得实际计算考虑的周期性和显示的结构的关系弄不清楚

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register