Theoretical UV circular dichroism of cyclo(L-Proline-L-Proline).

MP2, DFT, and molecular mechanics (AMBER, CVFF, and CFF91) geometry optimizations were performed on the cyclic dipeptide cyclo(L-Pro-L-Pro) starting from crystal structure data. Three stable conformations were identified as energy minima by all methods, but assignment of relative energy varied between the methods. The pi-pi transition feature of the UV circular dichroic (CD) spectrum was predicted for each minimized structure using the classical physics method of the dipole interaction model. The model was sensitive to the different conformations. The UV-CD predictions were compared individually and as a Boltzmann-weighted composite with published experimental CD spectra [Bowman, R. L.; Kellerman, M.; Johnson, W. C., Jr. Biopolymers 1983, 22, 1045]. For all structures, the original parameters of Applequist [Applequist, J. J. Chem. Phys. 1979, 71, 4324] with a bandwidth of 3000 cm(-1) most accurately replicated experiment, except for the CFF91 structures, which matched experiment best with a bandwidth of 4000 cm(-1). The inclusion of solvent by a continuum model did not significantly alter the minimized geometries obtained by molecular or quantum mechanics, but it did have an effect on the relative predicted energies of CFF91 and B3LYP conformations. The overall effect of solvent inclusion was negligible when Boltzmann-weighted spectra were considered. Gas-phase CFF91 structures were also reasonably good for prediction of CD spectra, and when water was included via a continuum model for energy calculations, the weighting scheme resembled that of the higher-level weightings. The CD calculated using the MP2/6-311G structures and energies for weighting were most descriptive of the 180 nm negative band in the experimental CD, but red-shifted the location of the 205 nm band. DFT structures were comparably, though not identically, as descriptive of the first pi-pi band, and did a better job with placement of the second (positive) pi-pi band. DFT calculations were less sensitive to basis set effect than the MP2 calculations, with 6-31G results in close agreement with 6-311G. The results suggest that it is possible to use geometries obtained from a variety of different methods (molecular mechanical or quantum mechanical) with the classical physics dipole interaction model to qualitatively reproduce the UV CD of model amides.