哇卡卡卡卡卡 发表于 2022-11-24 11:28
sobereva老师，您好！不好意思，基于前面的ONIOM（QM:MM）与QM/MM方法，我最近学习遇到一些疑惑。1、您前 ...

要善用簇模型，不要盲目用ONIOM算蛋白质-小分子相互作用问题
http://sobereva.com/597

sobereva 发表于 2022-11-4 08:55
1 ONIOM(QM:MM)和QM/MM都可以用，前者一大麻烦之处是你还得给底物分子定义力场参数

2 需要根据提示把缺 ...

sobereva老师，您好！不好意思，基于前面的ONIOM（QM:MM）与QM/MM方法，我最近学习遇到一些疑惑。1、您前面提到oniom方法计算需要补充参数，但是因为我需要去找过渡态，要在Gaussian中进行，然后要用oniom方法，现在遇到的问题就是缺失的参数应该从哪里去找呢，目前com文件中有些原子是缺少参数的，这个com文件我是从amber的rst文件转换而来的，是用sobtop软件中的参数还是什么呢？我查看了sobtop中有amber.itp文件，里面有一些atomtype，但是有很多，不知道该怎么选。（com文件如图所示）
2、我自己填了一下原子类型进去，然后运行的时候出现错误： MM function not complete，还是缺少参数，我看您前面回复其他人类似问题，需要补充参数，具体怎么补充去看exploring第三版的例子，老师，在网上并没有找到具体在gif怎么补充参数
3、因为我这个体系有2000多个原子，out文件报错缺失好几百个参数，这难道要一个一个补充嘛？4、老师，您说的用簇模型，这个具体是指什么呀，具体应该怎么操作呢，我已经用amber平衡好了底物+溶剂盒子（底物是含有金属的），然后需要去找一下过渡态，这里我除了想到用oniom分层的方法去计算在高斯中找过渡态，这种情况也可以用簇模型去计算吗，簇模型具体去哪里看呢，请老师指导！

sobereva 发表于 2022-11-6 09:57
1 如果只source了leaprc.gaff没source AMBER的，当前就是用的GAFF

2 取决于AMBER程序里怎么操作的

好的好的，谢谢sobereva老师

哇卡卡卡卡卡 发表于 2022-11-4 10:55
谢谢sobereva老师的指导！老师，我还有三个小问题：
①在文章中作者用Oniom(M06-2X/6-31G(d):amber)在高 ...

1 如果只source了leaprc.gaff没source AMBER的，当前就是用的GAFF

2 取决于AMBER程序里怎么操作的
Gaussian不自带GAFF力场，需要自己补充参数

3 对于有机结构，UFF质量整体不如GAFF，而且也没有标配的像样的算原子电荷的方法

sobereva 发表于 2022-11-4 08:55
1 ONIOM(QM:MM)和QM/MM都可以用，前者一大麻烦之处是你还得给底物分子定义力场参数

2 需要根据提示把缺 ...

谢谢sobereva老师的指导！老师，我还有三个小问题：
①在文章中作者用Oniom(M06-2X/6-31G(d):amber)在高斯上去计算过渡态，我了解到在amber程序中跑MD的时候，往往涉及到有机小分子配体时，会source leaprc.gaff，这里的GAFF力场，其函数形式和AMBER力场相同，与AMBER力场完全兼容（也就是说，GAFF力场和amber力场是并列关系）。可不可以理解为这里只是在amber程序中运行，但是真正用到的力场还是GAFF力场。
②如果在高斯中用Oniom(M06-2X/6-31G(d):gaff)去计算过渡态，然后在amber程序中跑MD，此时用的力场就是GAFF力场嘛，这里的Oniom(QM:MM)与QM/MM中的MM指得就是GAFF力场是吗？但是好像高斯内置的力场并不包含GAFF力场
③我们组用Oniom(QM:MM)分层方法处理大体系时(没有显性溶剂)，一般用到的力场都是uff力场，GAFF力场和uff力场处理溶剂二氯甲烷的话，这两者是不是都可以，没有什么特别大的好坏之分？

1 ONIOM(QM:MM)和QM/MM都可以用，前者一大麻烦之处是你还得给底物分子定义力场参数

2 需要根据提示把缺失的参数在gjf里都补充，让Gaussian读入。这就是我说的ONIOM用起来麻烦之处

3 如果MM区是蛋白质，在Gaussian自带的力场中AMBER是最佳选择。如果是其它分子，MM力场就不需要是AMBER了（很多情况也根本不适合AMBER），可以比如是GAFF，如果是水的话还可以用比如OPC3的专门的水模型。在同一篇文章里，AMBER程序里做QM/MM和Gaussian里做ONIOM(QM:MM)用的力场最好一样，省得别人comment。

4 AMBER不支持ONIOM。Gaussian不支持QM/MM。