簇模型计算酶催化底物水解机制的opt和TS结构问题

xiaoche

各位前辈好。我有个问题想请教一下。
我在用簇模型计算蛋白酶催化底物的机制，起始结构中底物由羰基氧、酰胺氮原子和羧基氧原子分别与酶活性中心的2个Zn2+配位。用的理论方法是B3LYP-D3，基组是6-311g(d,p) ，对Zn用SDD赝势基组。用了隐式溶剂（水）模型，共288个原子。用gaussian16做完opt任务后，freq无虚频，但是底物与锌的距离变大，按经验判断已经没有配位情况了。我使用opt后的结构，初猜TS构象进一步做TS时，全部退回opt结构。
我主要有以下问题：
（1）opt后得到的结构中，为什么底物与双锌离子的距离变大以至配位作用减弱甚至消失？是因为用了隐式溶剂+簇模型截断，对底物的真实环境模拟不足吗？
（2）TS退回opt结构，是否是opt的结构不是真实的ES构象的原因？
（3）我是否可以在opt时，对Zn与配位原子间的距离添加约束，来保持合理的几何参数？这样做是否合理？
（4）用的理论方法和基组是否可以降低一个标准，比如用B3LYP-D3，6-31g(d,p)，或者def2-SVP基组？300个原子目前算起来有些费时间。我的体系中，主要考虑的是底物与蛋白间的金属配位及弱相互作用（如氢键、Pi-Pi、疏水）。
（5）是否有必要用QM/MM？

另外，对于我的ES优化，起始构象是基于晶体结构的MD模拟聚类结构。在选定的理论方法和基组下首次进行opt时，出现多次振荡或者其它报错。我每次都是选取该次opt的能量最小的一帧结构，然后调整或者不调整参数，重新提交作业，最后经过5-6次后，opt收敛且freq无虚频。随帖子附上我的opt输出文件（文件太大上传的百度网盘：optfreq.out
链接: https://pan.baidu.com/s/1WiKDb-amqcZy6EjdH0FVEA?pwd=6w3w 提取码: 6w3w

sobereva

我没时间具体看你的文件。可能初始结构不合理。如果你确信那些化学键本该有，初始结构摆合适后，把相应键长冻结后优化，然后再放开冻结进一步优化，没准能维持住。
建议改用PBE0-D3(BJ)
可以尝试先用GFN2-xTB做预优化

不重要的原子可以降到6-31G*，甚至碳用6-31G往往也可以接受
没必要折腾QM/MM，折腾成本太高

xiaoche

谢谢卢老师！我按您的建议走一遍试试

Huschein

对于cluster来说泛函有点问题 如果要做预优化的话 xtb在重金属的表现上未必有MLFF好 不过你都可以用 MLFF最大的优势是跑的比xtb还要快一两个数量级 可以参考MAPLE 来做cluster预优化 https://pubs.rsc.org/en/content/articlelanding/2026/sc/d6sc01279e

QMMM对于你的体系来说我觉得基本不用想了 cluster都算的费劲的体系 QMMM除非缩小QM区 不然不可能算动，如果非要做QMMM可以考虑做MLMM 参考 https://www.nature.com/articles/s41467-026-72904-9

xiaoche

Huschein 发表于 2026-6-28 01:11
对于cluster来说泛函有点问题 如果要做预优化的话 xtb在重金属的表现上未必有MLFF好 不过你都可以用 MLFF最 ...

谢谢老师建议，我了解一下MLFF。
另外对于泛函的选择，卢老师建议PBE0-D3(BJ)，我看相关文献还有用M06-2X的。您有什么建议呢？我的体系涉及二价锌离子、有机分子，主要考虑的是金属配位，疏水、Pi-Pi、氢键等弱相互作用，主要涉及亲核进攻反应。
谢谢！

xiaoche

sobereva 发表于 2026-6-26 21:57
我没时间具体看你的文件。可能初始结构不合理。如果你确信那些化学键本该有，初始结构摆合适后，把相应键长 ...

卢老师，我在
“#p opt freq PBE1PBE/genecp scrf=(iefpcm,solvent=water) nosymm em=gd3bj

。。。
。。。

222 240 F
222 236 F
235 263 F

Zn 0
SDD
****
2,5-8,10-19,27,30-33,35-42,65,68-74,76-83,98,101-103,105-270,282-292 0
6-311g(d,p)
****
1,3-4,9,20-26,28-29,34,43-64,66-67,75,84-97,99-100,104,271-281 0
6-31g(d,p)
****

Zn 0
SDD”
下先做了预优化，冻结了边界原子和3个化学键键长；opt后做freq有一个大小为 -14 的虚频。
随后去除3个键长的冻结，再次在相同的理论方法和基组下进行opt，这次仅冻结了边界原子；opt进行了31步后，这3个键的距离依然变大，分别由 2.9——3.55，2.5——3.36，2.01——2.18。

请问这样的距离变化是否合理？
是否可以带着约束做优化呢，审稿人是否认可呢？
谢谢！