Vibrational spectral simulation for peptides of mixed secondary structure: method comparisons with the Trpzip model...

Infrared absorption and vibrational circular dichroism (IR and VCD) spectra of model fragments of TrpZip-style beta-hairpin structures are simulated using density functional theory (DFT) methods to estimate the influence of fragment size, end effects, conformational irregularities, peptide side chains, and solvent. Different fragmentation schemes, computing the strands and turn segments separately, were tested by varying the sizes of each and their respective overlaps. For suitably overlapping fragments, atomic property tensors were found to be reliably transferable, as tested by their ability to generate simulated spectra in good agreement with results from ab initio DFT computations for the entire peptide. This fragment approach significantly reduces computational times and opens up a wider range of systems that can be studied with a DFT-based approach as compared to previous methods based on uniform repeating sequences. However, vacuum calculations do not adequately represent the frequency dispersion of solvated molecules, and thus, some alternate strategies for solvation correction are explored for improving the simulation accuracy. Unlike for regular periodic secondary structure, the solvent significantly impacts the spectral shapes of hairpins, due to the different degrees of hydration of individual amide groups, which can be exposed to or shielded from water due to external vs internal hydrogen bonding. This is amplified by the shielding of selected amides from the solvent due to bulky side chains. The peptide plus solvent was structurally modeled with molecular dynamics methods, and then an electrostatic field-based parametrization correction was added to the force field and intensity tensors to compensate for the solvent dipolar field. The effect of the shielding and subsequent reordering of modes has a larger impact on VCD than IR band shapes.