PNNL

Abstract

Conformations of the important model system, the alanine dipeptide, have been calculated by using high level, ab initio electronic structure theory. A Ramachandran plot, with the angle f in the range [-180o,90o] and the angle y in the range [-60o,180o], was generated by using density functional theory with the generalized-gradient BLYP functional and a polarized triple-z basis set (TZVP+). Six conformers, C7eq, C5, C7ax, b2, aL, and a', have been identified in this region of the Ramachandran plot. A second derivative (frequency) analysis showed that all conformers are stable at this level of theory. These structures were used as starting points for geometry optimizations at the MP2/aug-cc-pVDZ level. Single point energies were calculated at the MP2/aug-cc-pVTZ and MP2/aug-cc-pVQZ levels at the final MP2/aug-cc-pVDZ structures and together with the MP2/aug-cc-pVDZ results were used in extrapolations to the complete basis set limit. The N-HoooO, N-HoooN, and C-HoooO hydrogen bond interactions that are key to the energetics are discussed. In general, the results obtained at the BLYP/TZVP+, MP2/aug-cc-pVDZ, MP2/aug-cc-pVTZ//aug-cc-pVDZ, and MP2/aug-cc-pVQZ//aug-cc-pVDZ levels are in reasonable agreement with each other, except for the b2 conformation where there are significant differences in the structures. Although the same stability order is obtained at all levels of theory that were used, there are significant differences in the magnitude of the relative conformational energies.