Toward Accurate Conformational Energies of Smaller Peptides and Medium-Sized Macrocycles: MPCONF196 Benchmark Energy Data Set

0489837 - ÚOCHB 2019 RIV US eng J - Journal Article
Řezáč, Jan - Bím, Daniel - Gutten, Ondrej - Rulíšek, Lubomír
Toward Accurate Conformational Energies of Smaller Peptides and Medium-Sized Macrocycles: MPCONF196 Benchmark Energy Data Set.
Journal of Chemical Theory and Computation. Roč. 14, č. 3 (2018), s. 1254-1266. ISSN 1549-9618. E-ISSN 1549-9626
R&D Projects: GA ČR(CZ) GA17-24155S; GA ČR(CZ) GJ16-11321Y
Institutional support: RVO:61388963
Keywords : density functional theory * main group thermochemistry * generalized gradient approximation
OECD category: Physical chemistry
Impact factor: 5.313, year: 2018

A carefully selected set of acyclic and cyclic model peptides and several other macrocycles, comprising 13 compounds in total, has been used to calibrate the accuracy of the DFT(-D3) method for conformational energies, employing BP86, PBEO, PBE, B3LYP, BLYP, TPSS, TPSSh, M06-2X, B97-D, OLYP, revPBE, M06-L, SCAN, revTPSS, BH-LYP, and omega B97X-D3 functionals. Both high-and low-energy conformers, 15 or 16 for each compound adding to 196 in total, denoted as the MPCONF196 data set, were included, and the reference values were obtained by the composite protocol, yielding the CCSD(T)/CBS extrapolated energies or their DLPNO-CCSD(T)/CBS equivalents in the case of larger systems. The latter was shown to be in near-quantitative (similar to 0.10-0.15 kcal.mol(-1)) agreement with the canonical CCSD(T), provided the TightPNO setting is used, and, therefore, can be used as the reference for larger systems (likely up to 150-200 atoms) for the problem studied here. At the same time, it was found that many D3-corrected DFT functionals provide results of similar to 1 kcal.mol(-1) accuracy, which we consider as quite encouraging. This result implies that DFT-D3 methods can be, for example, safely used in efficient conformational sampling algorithms. Specifically, the DFT-D3/DZVP-DFT level of calculation seems to be the best trade-off between computational cost and accuracy. Based on the calculated data, we have not found any cheaper variant for the treatment of conformational energies, since the semiempirical methods (including DFTB) provide results of inferior accuracy (errors of 3-5 kcal.mol(-1)).
Permanent Link: http://hdl.handle.net/11104/0284170