Validation challenge of density-functional theory for peptides-example of Ac-Phe-Ala5-LysH(+).

Journal Article

We assess the performance of a group of exchange-correlation functionals for predicting the secondary structure of peptide chains, up to a new many-body dispersion corrected hybrid density functional, dubbed PBE0+MBD* by its original authors. For the purpose of validation, we first compare to published, high-level benchmark conformational energy hierarchies (coupled cluster at the singles, doubles, and perturbative triples level, CCSD(T)) for 73 conformers of small three-residue peptides, establishing that the van der Waals corrected PBE0 functional yields an average error of only ∼20 meV (∼0.5 kcal/mol). This compares to ∼40-50 meV for nondispersion corrected PBE0 and 40-100 meV for different empirical force fields (estimated for the alanine tetrapeptide). For longer peptide chains that form a secondary structure, CCSD(T) level benchmark data are currently unaffordable. We thus turn to the experimentally well studied Ac-Phe-Ala5-LysH(+) peptide, for which four closely competing conformers were established by infrared spectroscopy. For comparison, an exhaustive theoretical conformational space exploration yields at least 11 competing low energy minima. We show that (i) the many-body dispersion correction, (ii) the hybrid functional nature of PBE0+MBD*, and (iii) zero-point corrections are needed to reveal the four experimentally observed structures as the minima that would be populated at low temperature.

Full Text

Duke Authors

Cited Authors

  • Rossi, M; Chutia, S; Scheffler, M; Blum, V

Published Date

  • September 4, 2014

Published In

Volume / Issue

  • 118 / 35

Start / End Page

  • 7349 - 7359

PubMed ID

  • 24405171

Electronic International Standard Serial Number (EISSN)

  • 1520-5215

Digital Object Identifier (DOI)

  • 10.1021/jp412055r

Language

  • eng

Conference Location

  • United States