Plett Christoph, Grimme Stefan, Hansen Andreas
Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstraße 4, 53115 Bonn, Germany.
J Chem Theory Comput. 2024 Sep 11. doi: 10.1021/acs.jctc.4c00801.
Simulating peptides and proteins is becoming increasingly important, leading to a growing need for efficient computational methods. These are typically semiempirical quantum mechanical (SQM) methods, force fields (FFs), or machine-learned interatomic potentials (MLIPs), all of which require a large amount of accurate data for robust training and evaluation. To assess potential reference methods and complement the available data, we introduce two sets, DipCONFL and DipCONFS, which cover large parts of the conformational space of 17 amino acids and their 289 possible dipeptides in aqueous solution. The conformers were selected from the exhaustive PeptideCS dataset by Andris et al. [ 2022, 126, 5949-5958]. The structures, originally generated with GFN2-xTB, were reoptimized using the accurate rSCAN-3c density functional theory (DFT) composite method including the implicit CPCM water solvation model. The DipCONFS benchmark set contains 918 conformers and is one of the largest sets with highly accurate coupled cluster conformational energies so far. It is employed to evaluate various DFT and wave function theory (WFT) methods, especially regarding whether they are accurate enough to be used as reliable reference methods for larger datasets intended for training and testing more approximated SQM, FF, and MLIP methods. The results reveal that the originally provided BP86-D3(BJ)/DGauss-DZVP conformational energies are not sufficiently accurate. Among the DFT methods tested as an alternative reference level, the revDSD-PBEP86-D4 double hybrid performs best with a mean absolute error (MAD) of 0.2 kcal mol compared with the PNO-LCCSD(T)-F12b reference. The very efficient rSCAN-3c composite method also shows excellent results, with an MAD of 0.3 kcal mol, similar to the best-tested hybrid ωB97M-D4. With these findings, we compiled the large DipCONFL set, which includes over 29,000 realistic conformers in solution with reasonably accurate rSCAN-3c reference conformational energies, gradients, and further properties potentially relevant for training MLIP methods. This set, also in comparison to DipCONFS, is used to assess the performance of various SQM, FF, and MLIP methods robustly and can complement training sets for those.
模拟肽和蛋白质变得越来越重要,这导致对高效计算方法的需求不断增长。这些方法通常是半经验量子力学(SQM)方法、力场(FF)或机器学习原子间势(MLIP),所有这些方法都需要大量准确的数据进行稳健的训练和评估。为了评估潜在的参考方法并补充可用数据,我们引入了两组数据,DipCONFL和DipCONFS,它们涵盖了17种氨基酸及其在水溶液中的289种可能二肽的大部分构象空间。这些构象是从Andris等人[2022, 126, 5949 - 5958]详尽的PeptideCS数据集中选择的。最初使用GFN2 - xTB生成的结构,使用包括隐式CPCM水溶剂化模型的精确rSCAN - 3c密度泛函理论(DFT)复合方法进行了重新优化。DipCONFS基准集包含918个构象,是迄今为止具有高精度耦合簇构象能量的最大数据集之一。它用于评估各种DFT和波函数理论(WFT)方法,特别是关于它们是否准确到足以用作更大数据集的可靠参考方法,这些数据集旨在训练和测试更近似的SQM、FF和MLIP方法。结果表明,最初提供的BP86 - D3(BJ)/DGauss - DZVP构象能量不够准确。在作为替代参考水平测试的DFT方法中,revDSD - PBEP86 - D4双杂化方法表现最佳,与PNO - LCCSD(T) - F12b参考相比,平均绝对误差(MAD)为0.2 kcal/mol。非常高效的rSCAN - 3c复合方法也显示出优异的结果,MAD为0.3 kcal/mol,与测试最好的杂化ωB97M - D4相似。基于这些发现,我们编制了大型DipCONFL集,其中包括超过29,000个溶液中的实际构象,具有合理准确的rSCAN - 3c参考构象能量、梯度以及可能与训练MLIP方法相关的其他性质。与DipCONFS相比,这个数据集也用于稳健地评估各种SQM、FF和MLIP方法的性能,并可以补充这些方法的数据训练集。