Ab initio quantum mechanical models of peptide helices and their vibrational spectra

Structural parameters for standard peptide helices (alpha, 3(10), 3(1) left-handed) were fully ab initio optimized for Ac-(L-Ala)(9)-NHMe and for Ac-(L-Pro)(9)-NHMe (poly-L-proline-PLP I and PLP II-forms), in order to better understand the relative stability and minimum energy geometries of these conformers and the dependence of the ir absorption and vibrational CD (VCD) spectra on detailed variation in these conformations. Only the 3(10)-helical Ala-based conformation was stable in vacuum for this decaamide structure, but both Pro-based conformers minimized successfully. Inclusion of solvent effects, by use of the conductor-like screening solvent model (COSMO), enabled ab initio optimizations [at the DFT/B3LYP/SV(P) level] without any constraints for the alpha- and 3(10)-helical Ala-based peptides as well as the two Pro-based peptides. The geometries obtained compare well with peptide chain torsion angles and hydrogen-bond distances found for these secondary structure types in x-ray structures of peptides and proteins. For the simulation of VCD spectra, force field and intensity response tensors were obtained ab initio for the complete Ala-based peptides in vacuum, but constrained to the COSMO optimized torsional angles, due to limitations of the solvent model. Resultant spectral patterns reproduce well many aspects of the experimental spectra and capture the differences observed for these various helical types.