Unveiling Vibrational Couplings in Model Peptides in Solution by a Theoretical Approach

Vibrational analysis of peptides in solution and the theoretical determination of the effects of the microenvironment on infrared and Raman spectra are of key importance in many fields of chemical interest. In this work, we present a computational study combining static quantum mechanical calculations with ab initio molecular dynamics simulations to investigate the vibrational behavior of three peptide models in both the gas phase and in explicit water, under non-periodic boundary conditions. The vibrational spectra of the main amide bands, namely amide I-III and A, were analyzed using a time-frequency approach based on the wavelet transform, which allows the resolution of transient frequency shifts and mode couplings along the trajectories. This combined approach enabled us to perform a time-resolved vibrational analysis revealing how vibrational frequencies, especially of the C=O and N-H stretching modes, evolve over time due to dynamical microsolvation. These fluctuations modulate vibrational couplings and lead to spectral broadening and frequency shifts that correlate with the local structuring of the solvent. In conclusion, our results highlight how the proposed protocol allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, the peptide backbone, and its microenvironment.