计算化学公社

标题: 分子动力学模拟锌离子的配位数限制问题 [打印本页]

作者
Author: poseur 时间: 2020-6-12 20:12
标题: 分子动力学模拟锌离子的配位数限制问题
我做的是锌酶体系，其晶体结构中金属中心原子配位数为四，我用Amber力场跑分子动力学总是得到六配位结构（多了两个水分子），用QM/MM MD方法验证四配位为稳定结构，请问用Amber力场有什么办法能限制住锌离子的配位数么？

作者
Author: mol 时间: 2020-6-12 20:13
分子力学没办法模拟配位键吧

作者
Author: poseur 时间: 2020-6-12 20:23

mol 发表于 2020-6-12 20:13
分子力学没办法模拟配位键吧

我模拟的是整个蛋白酶体系，不是模拟配位键，AMBER力场有金属锌的力场参数，但是其在水盒子中里弛豫后的得到的稳定结构好像总是倾向于六配位。

作者
Author: chenjinfeng850 时间: 2020-6-12 20:59

poseur 发表于 2020-6-12 20:23
我模拟的是整个蛋白酶体系，不是模拟配位键，AMBER力场有金属锌的力场参数，但是其在水盒子中里弛豫后的 ...

你是用的MCPB.py 建立的Zn的立场参数吗？
Zn上带的正电荷，水中的氧带的是负电荷，在跑MD的时候会吸在一起应该是正常现象。
简单粗暴的办法是你直接把Zn附近的水都删掉。
或者直接QM/MM MD

作者
Author: k64_cc 时间: 2020-6-12 21:06
Amber力场里离子的intermeoecular interaction都是各向同性的，而且也没特殊tune过，配位数不对也挺正常的。你都跑QM/MM了，就一直跑着吧……

作者
Author: sobereva 时间: 2020-6-12 23:06
适当调整键角力场参数之类，以维持MD过程中局部的四配位结构不变。没具体结构图不好说细节。

作者
Author: fhh2626 时间: 2020-6-13 09:44
可以用AMBER里面的MCPB固定配位数，也可以直接用固定配位数的Dummy-atom force field

作者
Author: poseur 时间: 2020-6-23 08:46
用MCPB.py解决了，谢谢各位

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