DFT-Based Simulations of Amide I′ IR Spectra of a Small Protein in Solution Using Empirical Electrostatic Map with a Continuum Solvent Model

Computing infrared spectra of proteins using the exciton model

The ability to compute from first principles the infrared spectrum of a protein in solution phase representing a biological system would provide a useful connection to atomistic models of protein structure and dynamics. Indeed, such calculations are a vital complement to 2DIR experimental measurements, allowing the observed signals to be interpreted in terms of detailed structural and dynamical information. In this article, we have studied nine structurally and spectroscopically well-characterized proteins, representing a range of structural types. We have simulated the equilibrium conformational dynamics in an explicit point charge water model. Using the resulting trajectories based on MD simulations, we have computed the one and two dimensional infrared spectra in the Amide I region, using an exciton approach, in which a local mode basis of carbonyl stretches is considered. The role of solvent in shifting the Amide I band (by 30 to 50 cm-1 ) is clearly evident. Similarly, the conformational dynamics contribute to the broadening of peaks in the spectrum. The inhomogeneous broadening in both the 1D and 2D spectra reflects the significant conformational diversity observed in the simulations. Through the computed 2D cross-peak spectra, we show how different pulse schemes can provide additional information on the coupled vibrations. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Temperature dependence of C-terminal carboxylic group IR absorptions in the amide I' region

Studies of structural changes in peptides and proteins using IR spectroscopy often rely on subtle changes in the amide I' band as a function of temperature. However, these changes can be obscured by the overlap with other absorptions, namely the side-chain and terminal carboxylic groups. The former were the subject of our previous report (Anderson et al., 2014). In this paper we investigate the IR spectra of the asymmetric stretch of α-carboxylic groups for amino acids representing all major types (Gly, Ala, Val, Leu, Ser, Thr, Asp, Glu, Lys, Asn, His, Trp, Pro) as well as the C-terminal groups of three dipeptides (Gly-Gly, Gly-Ala, Ala-Gly) in D₂O at neutral pH. Experimental temperature dependent IR spectra were analyzed by fitting of both symmetric and asymmetric pseudo-Voigt functions. Qualitatively the spectra exhibit shifts to higher frequency, loss in intensity and narrowing with increased temperature, similar to that observed previously for the side-chain carboxylic groups of Asp. The observed dependence of the band parameters (frequency, intensity, width and shape) on temperature is in all cases linear: simple linear regression is therefore used to describe the spectral changes. The spectral parameters vary between individual amino acids and show systematic differences between the free amino acids and dipeptides, particularly in the absolute peak frequencies, but the temperature variations are comparable. The relative variations between the dipeptide spectral parameters are most sensitive to the C-terminal amino acid, and follow the trends observed in the free amino acid spectra. General rules for modeling the α-carboxylic IR absorption bands in peptides and proteins as the function of temperature are proposed.

Infrared, vibrational circular dichroism, and Raman spectral simulations for β-sheet structures with various isotopic labels, interstrand, and stacking arrangements using density functional theory

Infrared (IR), Raman, and vibrational circular dichroism (VCD) spectral variations for different β-sheet structures were studied using simulations based on density functional theory (DFT) force field and intensity computations. The DFT vibrational parameters were obtained for β-sheet fragments containing nine-amides and constrained to a variety of conformations and strand arrangements. These were subsequently transferred onto corresponding larger β-sheet models, normally consisting of five strands with ten amides each, for spectral simulations. Further extension to fibril models composed of multiple stacked β-sheets was achieved by combining the transfer of DFT parameters for each sheet with dipole coupling methods for interactions between sheets. IR spectra of the amide I show different splitting patterns for parallel and antiparallel β-sheets, and their VCD, in the absence of intersheet stacking, have distinct sign variations. Isotopic labeling by (13)C of selected residues yields spectral shifts and intensity changes uniquely sensitive to relative alignment of strands (registry) for antiparallel sheets. Stacking of multiple planar sheets maintains the qualitative spectral character of the single sheet but evidences some reduction in the exciton splitting of the amide I mode. Rotating sheets with respect to each other leads to a significant VCD enhancement, whose sign pattern and intensity is dependent on the handedness and degree of rotation. For twisted β-sheets, a significant VCD enhancement is computed even for sheets stacked with either the same or opposite alignments and the inter-sheet rotation, depending on the sense, can either further increase or weaken the enhanced VCD intensity. In twisted, stacked structures (without rotation), similar VCD amide I patterns (positive couplets) are predicted for both parallel and antiparallel sheets, but different IR intensity distributions still enable their differentiation. Our simulation results prove useful for interpreting experimental vibrational spectra in terms of β-sheet and fibril structure, as illustrated in the accompanying paper.

Structural analyses of experimental 13C edited amide I' IR and VCD for peptide β-sheet aggregates and fibrils using DFT-based spectral simulations

In the preceding paper, computational models based on density functional theory (DFT) were presented to characterize the sensitivity of vibrational spectroscopic methods (IR, VCD, and Raman) to structural features of β-sheets. Isotopically edited amide I' IR for peptides labeled with (13)C in multiple different sites provides the most structurally distinct signatures of strand alignment, while VCD is sensitive to the sheet twist and intersheet stacking. In this report, we simulate the IR and VCD spectra for models approximating structures of four β-sheet forming peptides previously experimentally studied using these methods with (13)C isotopic editing. Various register alignments are tested. Agreement with experiment is evaluated based on frequency shifts of both the (12)C and (13)C IR amide I' signals, relative intensity patterns, and VCD spectra where available. While for the simulation of IR spectra canonical planar sheets provide a sufficient model system, for VCD simulation twisted, stacked sheets are required in order to reproduce strong couplet-like amide I' VCD. Effects of the solvent (water) and amino acid side chains are also tested by using a simplified, electrostatic solvent model and atomic partial charges for the side chains. Very good agreement with experimental spectra is obtained, particularly for the relative (12)C and (13)C band frequencies. All four peptide models are shown to be antiparallel as had previously been assumed. However, in some cases our simulations are consistent with different register alignment of strands than originally proposed.