计算化学公社

标题: 求助关于蛋白配体结合MD的问题 [打印本页]

作者
Author: 溪临昭煌 时间: 2021-11-3 00:35
标题: 求助关于蛋白配体结合MD的问题
各位老师好，我使用gromacs对我的酶——底物复合体系进行MD，体系来源于共晶体，使用的力场是CHARMM36和CGENFF，配体的电荷使用了RESP电荷，在300K，1bar放开限制跑了15ns以后发现小分子跑出了活性中心，开始在溶剂里面乱窜，目前排除了电荷设置错误的问题和结合位点的问题，请问各位老师这是什么原因造成的，该如何解决呢？
(, 下载次数 Times of downloads: 6)

作者
Author: sobereva 时间: 2021-11-3 00:42
mdp没明显问题，从这上看不出什么来。没有更多信息没法判断

作者
Author: k64_cc 时间: 2021-11-3 10:13
CGENFF和RESP组合本身就不合适。

CHARMM系列的charge是拟合分子和水的相互作用得到的，不是重现ESP。看原文献Figure 2。
https://onlinelibrary.wiley.com/doi/10.1002/jcc.21367

作者
Author: 溪临昭煌 时间: 2021-11-3 10:49

sobereva 发表于 2021-11-3 00:42
mdp没明显问题，从这上看不出什么来。没有更多信息没法判断

谢谢老师，这是我的配体的str文件， (, 下载次数 Times of downloads: 4) ，还想请问一下老师如果模拟正确，是否蛋白和配体会始终保持结合状态？

作者
Author: 溪临昭煌 时间: 2021-11-3 10:50

k64_cc 发表于 2021-11-3 10:13
CGENFF和RESP组合本身就不合适。

CHARMM系列的charge是拟合分子和水的相互作用得到的，不是重现ESP。看 ...

谢谢老师回复，我先使用了cgenff自带的电荷进行模拟，但是出现了配体飞出，因此更换了RESP，但是也飞出了，不过RESP电荷飞出去时候模拟的时间确实比用自带的电荷进行模拟要短

作者
Author: k64_cc 时间: 2021-11-3 13:28

溪临昭煌发表于 2021-11-3 10:50
谢谢老师回复，我先使用了cgenff自带的电荷进行模拟，但是出现了配体飞出，因此更换了RESP，但是也飞出了 ...

有没有那么一种可能，这俩事本身没关系呢？

要不你换Amber和GAFF吧。

作者
Author: 溪临昭煌 时间: 2021-11-3 15:10

k64_cc 发表于 2021-11-3 13:28
有没有那么一种可能，这俩事本身没关系呢？

要不你换Amber和GAFF吧。

确实有可能...不过这样对CHARMM就有点儿心理阴影了...我试试去，谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)