Interaction and dimerization energies in methyl-blocked alpha,gamma-peptide nanotube segments

The building blocks of a promising class of peptide nanotubes composed of alternating D-alpha-amino acids and (1R,3S)-3-aminocyclohexane (or cyclopentane) carboxylic acid (D-gamma-Ach or D-gamma-Acp) were explored by computational methods. Specifically, density functional theory (DFT) calculations on monomers and dimers of gamma-Ach-based and gamma-Acp-based alpha,gamma-cyclo-hexapeptides and cyclo-octapeptides were carried out to investigate the experimentally observed preference for alpha-alpha over gamma-gamma dimerization, associated with the two types of stacking patterns present in these peptide nanotubes, as well as the preference for heterodimerization versus homodimerization. Full geometry optimizations were performed at the B3LYP/6-31G(d) level, and single point calculations were subsequently carried out with the B3LYP and M05-2X functionals and the 6-31+G(d,p) basis set. The calculations predict that the interaction energies in the alpha-alpha species are quite similar to those in the gamma-gamma dimers. However, a comparison of dimerization energies (i.e., interaction energies plus deformation energies of monomers) shows that alpha-alpha dimerization is energetically favored over gamma-gamma dimerization. The calculations strongly suggest that the preference for alpha-alpha binding is governed by differences between the deformation energies in the alpha and gamma monomers, rather than by differences between the relative strengths of the alpha-alpha and gamma-gamma hydrogen-bonding patterns. Calculations based on local properties of the electron density support the previous suggestion that the H-N bonds of the alpha-amino acids are more polarized than those of the gamma-amino acids.