求助：蛋白-配体复合物在进行100ns的模拟后发现到60ns左右配体脱离结合口袋

linhan · 发表于 Post on 2024-10-25 15:55:01

我的体系是蛋白A-配体-蛋白B的模型，在进行初步的10ns模拟后较为稳定，于是用gromacs进行100ns的模拟，结果发现60ns左右有一端的配体从结合口袋脱离出来，导致最后蛋白体系就解离了，但是配体的这一段和结合口袋的互作是比较强（biotin和streptavidin），按理来说应该是不会脱离，我不清楚是怎么一回事，我用的力场是AMBER FF99SB，分别做了100ps的nvt和npt，放在cubic水盒子里面。

分子模拟刚刚入门，恳请各位大佬给予指点

sobereva · 发表于 Post on 2024-10-26 00:03:18

谈谈分子动力学模拟蛋白-配体复合物过程中配体发生脱离的原因
http://sobereva.com/632（http://bbs.keinsci.com/thread-27434-1-1.html）

linhan · 发表于 Post on 2024-10-26 14:49:28

sobereva 发表于 2024-10-26 00:03
谈谈分子动力学模拟蛋白-配体复合物过程中配体发生脱离的原因
http://sobereva.com/632（http://bbs.keins ...

谢谢sob老师，这个帖子在您发之前我已经学习过了，我准备换成amber 14sb力场去跑一下看看是不是力场的问题。因为这个体系别人用amber24中的amber14sb力场跑过，那个结果是0-100ns内配体和结合口袋的结合情况都是良好的，数据分析的rmsd和rmsf也是比较理想的。但我是用的gromacs的amberff99sb跑的。至于您说的“做生物体系的动力学我一般建议让参考温度从0 K开始缓慢线性上升，这样最为稳妥。”我不知道我的mdp文件有没有这样涉及到，但是在60ns后才发生配体跑开的话可能不是因为这个问题？还有您提到的可能是模型不合理甚至生物互作方面的不合理我觉得可能也是可以排除的，因为类似的复合物结构现实合成出来过。

HNUST · 发表于 Post on 2024-10-26 21:02:27

linhan 发表于 2024-10-26 14:49
谢谢sob老师，这个帖子在您发之前我已经学习过了，我准备换成amber 14sb力场去跑一下看看是不是力场的问 ...

99sb-ildn好像有点过时了。我有一个配体跑了三次结合很稳定，结果突发奇想换个14sb直接跑出去了，用desmond也显示不稳定，震惊

linhan · 发表于 Post on 2024-10-27 12:59:54

HNUST 发表于 2024-10-26 21:02
99sb-ildn好像有点过时了。我有一个配体跑了三次结合很稳定，结果突发奇想换个14sb直接跑出去了，用desm ...

是的我也觉得有点过时了，所以用14sb再跑一次，但是你说14sb也会跑出去

。因为用amber24的14sb力场跑是似乎是没有问题的，就是不清楚gromacs跑怎么就寄了

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

[GROMACS] 求助：蛋白-配体复合物在进行100ns的模拟后发现到60ns左右配体脱离结合口袋