Relating the infrared absorption of proteins to their structure - dependence of the transition dipole moment on the protein structure
||Relating the infrared absorption of proteins to their structure - dependence of the transition dipole moment on the protein structure|
||SNIC Medium Compute|
||Andreas Barth <Andreas.Barth@dbb.su.se>|
||2021-07-27 – 2021-12-01|
||10603 10402 10302|
The leading application of infrared spectroscopy in the life sciences is the analysis of protein structure. It exploits the conformational sensitivity of the amide I vibrations, which are collective vibrations of the amide groups of the protein backbone. However, their theoretical description is currently limited and therefore the interpretation of experimental results is guided by (sometimes conflicting) empirical correlations. Therefore, a reliable theoretical description is desirable and can lead to additional structural insight (e.g. Baronio, Baldassarre, & Barth, PCCP 2019).
The large number of atoms in proteins prevents calculating their vibrational spectrum from density functional theory (DFT). Therefore, present calculations of the amide I spectrum use a different approach. They mostly consider each amide group as a local amide I oscillator, which has an intrinsic frequency and which couples with other amide groups. The intrinsic frequencies are affected by the local backbone conformation and electrostatic effects including hydrogen bonding. The former influences also the coupling between nearest neighbors. Long-range coupling is described by transition dipole coupling (TDC) between the transition dipole moments (TDMs) of close-by amide I oscillators. Once the intrinsic frequencies and the coupling constants are known, the spectrum of the delocalized amide I normal modes can be calculated. This approach describes general trends well, but needs to be improved for detailed structural analysis. These improvements rely on DFT calculations of small peptide systems (up to 21 amide groups). We used the last round of our calculations to improve the description of the intrinsic frequencies and to assess the magnitude of mechanical coupling (i.e. coupling via the oscillating atoms) between adjacent amide groups. In the latter study, we obtained experience in isolating the vibrations of specific amide groups by increasing the mass of adjacent atoms. We will use this approach to complete our original work plan with the following studies that are fundamental for the calculation of the amide I spectrum:
a) An effect of conformation on the magnitude and direction of the TDM was found in our recent DFT calculations (Baronio & Barth, J. Phys. Chem. B 2020). Implementing secondary structure dependent TDM parameters into our amide I calculations might cause abrupt changes of these parameters upon relatively small changes in backbone structure. To avoid that, we will explore whether the different dihedral angles are the underlying reason for the different TDM parameters using DFT calculations of tetrapeptides (three complete amide groups) in different conformations. We will start with model structures of helices, parallel and antiparallel β-sheets for which geometry optimized structures are already available from our previous work.
b) An effect of hydrogen bonding on the TDM parameters is also indicated by our previous work and we will use the same approach to study it, but using larger structures (10–21 amide groups), which are already geometry optimized.