Computational spectroscopy of ubiquitin: comparison between theory and experiments

Using the constrained molecular dynamics simulation method in combination with quantum chemistry calculation, Hessian matrix reconstruction, and fragmentation approximation methods, the authors have established computational schemes for numerical simulations of amide I IR absorption, vibrational circular dichroism (VCD), and two-dimensional (2D) IR photon echo spectra of the protein ubiquitin in water. Vibrational characteristic features of these spectra in the amide I vibration region are discussed. From the semiempirical quantum chemistry calculation results on an isolated ubiquitin, amide I local mode frequencies and vibrational coupling constants were fully determined. It turns out that the amide I local mode frequencies of ubiquitin in both gas phase and aqueous solution are highly heterogeneous and site dependent. To directly test the quantitative validity of thus obtained spectroscopic properties, they compared the experimentally measured amide I IR, 2D IR, and electronic circular dichroism spectra with experiments, and found good agreements between theory and experiments. However, the simulated VCD spectrum is just qualitatively similar to the experimentally measured one. This indicates that, due to delicate cancellations between the positive and negative VCD contributions, the prediction of protein VCD spectrum is critically relied on quantitative accuracy of the theoretical model for predicting amide I local mode frequencies. On the basis of the present comparative investigations, they found that the site dependency of amide I local mode frequency, i.e., diagonal heterogeneity of the vibrational Hamiltonian matrix in the amide I local mode basis, is important. It is believed that the present computational methods for simulating various vibrational and electronic spectra of proteins will be of use in further refining classical force fields and in addressing the structure-spectra relationships of proteins in solution.     

Simulation of vibrational energy transfer in two-dimensional infrared spectroscopy of amide I and amide II modes in solution

Population transfer between vibrational eigenstates is important for many phenomena in chemistry. In solution, this transfer is induced by fluctuations in molecular conformation as well as in the surrounding solvent. We develop a joint electrostatic density functional theory map that allows us to connect the mixing of and thereby the relaxation between the amide I and amide II modes of the peptide building block N-methyl acetamide. This map enables us to extract a fluctuating vibrational Hamiltonian from molecular dynamics trajectories. The linear absorption spectrum, population transfer, and two-dimensional infrared spectra are then obtained from this Hamiltonian by numerical integration of the Schrodinger equation. We show that the amide I/amide II cross peaks in two-dimensional infrared spectra in principle allow one to follow the vibrational population transfer between these two modes. Our simulations of N-methyl acetamide in heavy water predict an efficient relaxation between the two modes with a time scale of 790 fs. This accounts for most of the relaxation of the amide I band in peptides, which has been observed to take place on a time scale of 450 fs in N-methyl acetamide. We therefore conclude that in polypeptides, energy transfer to the amide II mode offers the main relaxation channel for the amide I vibration.     

Vibrational relaxation pathways of amide I and amide II modes in N-methylacetamide

We studied the vibrational energy relaxation mechanisms of the amide I and amide II modes of N-methylacetamide (NMA) monomers dissolved in bromoform using polarization-resolved femtosecond two-color vibrational spectroscopy. The results show that the excited amide I vibration transfers its excitation energy to the amide II vibration with a time constant of 8.3 ± 1 ps. In addition to this energy exchange process, we observe that the excited amide I and amide II vibrations both relax to a final thermal state. For the amide I mode this latter process dominates the vibrational relaxation of this mode. We find that the vibrational relaxation of the amide I mode depends on frequency which can be well explained from the presence of two subbands with different vibrational lifetimes (~1.1 ps on the low frequency side and ~2.7 ps on the high frequency side) in the amide I absorption spectrum.     

Residual native structure in a thermally denatured beta-hairpin

We investigate the thermal denaturation of trpzip2 between 15 and 82 degrees C using two-dimensional infrared (2D IR) vibrational spectroscopy, dispersed vibrational echo (DVE) spectroscopy, and Fourier transform infrared (FTIR) spectroscopy. The FTIR and DVE spectra of trpzip2 show in the amide I region of the spectrum two resonances, which arise primarily from the interstrand coupling between local amide I oscillators along the peptide backbone. The coupling is seen directly in the 2D IR spectra as the formation of cross-peak ridges. Although small shifts of these frequencies occur on heating the sample, the existence of cross-peak ridges at all temperatures indicates that stable hydrogen bond interactions persist between the two beta-strands. These observations indicate a significant amount of native structure in the thermally denatured state of trpzip2.     

Two-dimensional infrared spectroscopy of antiparallel beta-sheet secondary structure

We investigate the sensitivity of femtosecond Fourier transform two-dimensional infrared spectroscopy to protein secondary structure with a study of antiparallel beta-sheets. The results show that 2D IR spectroscopy is more sensitive to structural differences between proteins than traditional infrared spectroscopy, providing an observable that allows comparison to quantitative models of protein vibrational spectroscopy. 2D IR correlation spectra of the amide I region of poly-l-lysine, concanavalin A, ribonuclease A, and lysozyme show cross-peaks between the IR-active transitions that are characteristic of amide I couplings for polypeptides in antiparallel hydrogen-bonding registry. For poly-l-lysine, the 2D IR spectrum contains the eight-peak structure expected for two dominant vibrations of an extended, ordered antiparallel beta-sheet. In the proteins with antiparallel beta-sheets, interference effects between the diagonal and cross-peaks arising from the sheets, combined with diagonally elongated resonances from additional amide transitions, lead to a characteristic "Z"-shaped pattern for the amide I region in the 2D IR spectrum. We discuss in detail how the number of strands in the sheet, the local configurational disorder in the sheet, the delocalization of the vibrational excitation, and the angle between transition dipole moments affect the position, splitting, amplitude, and line shape of the cross-peaks and diagonal peaks.     

Effects of intermolecular vibrational coupling and liquid dynamics on the polarized Raman and two-dimensional infrared spectral profiles of liquid N,N-dimethylformamide analyzed with a time-domain computational method

A time-domain method for calculating polarized Raman and two-dimensional infrared (2D-IR) spectra that includes the effects of both the diagonal frequency modulations (of individual molecules in the system) and the off-diagonal (intermolecular) vibrational coupling is presented and applied to the case of the amide I band of liquid N,N-dimethylformamide. It is shown that the effect of the resonant off-diagonal vibrational coupling and the resulting delocalization of vibrational modes is clearly seen as the noncoincidence effect in the polarized Raman spectrum and some spectral features (especially as asymmetric intensity patterns) in the 2D-IR spectra. The type of 2D-IR spectra (concerning the polarization condition) most appropriate for observing this effect is discussed. On the basis of the agreement between the observed and calculated band profiles of the polarized Raman spectrum, the time dependence of the transient IR absorption anisotropy is also calculated. The method of evaluating the extent of delocalization of vibrational modes that is relevant to the features of these optical signals in the time and frequency domains is discussed. The nature of the molecular motions (concerning the liquid dynamics) that are effective on the diagonal frequency modulations is also examined.     

Calculations of intermode coupling constants and simulations of amide I,II, and III vibrational spectra of dipeptides

Amide I, II, and III vibrations of polypeptides are important marker modes whose vibrational spectra can provide critical information on structure and dynamics of proteins in solution. The extent of delocalization and vibrational properties of amide normal mode can be described by the amide local mode frequencies and intermode coupling constants between a pair of amide local modes. To determine these fundamental quantities, the previous Hessian matrix reconstruction method has been generalized here and applied to the density functional theory results for various dipeptide conformers. The calculation results are then used to simulate IR absorption, vibrational circular dichroism, and 2D IR spectra of dipeptides. The relationships between dipeptide backbone conformations and these vibrational spectra are discussed. It is believed that the present computational method and results will be of use to quantitatively simulate vibrational spectra of complicated polypeptides beyond simple dipeptides     

Site-specific vibrational dynamics of the CD3zeta membrane peptide using heterodyned two-dimensional infrared photon echo spectroscopy

Heterodyned two-dimensional infrared (2D IR) spectroscopy has been used to study the amide I vibrational dynamics of a 27-residue peptide in lipid vesicles that encompasses the transmembrane domain of the T-cell receptor CD3zeta. Using 1-(13)C[Double Bond](18)O isotope labeling, the amide I mode of the 49-Leucine residue was spectroscopically isolated and the homogeneous and inhomogeneous linewidths of this mode were measured by fitting the 2D IR spectrum collected with a photon echo pulse sequence. The pure dephasing and inhomogeneous linewidths are 2 and 32 cm(-1), respectively. The population relaxation time of the amide I band was measured with a transient grating, and it contributes 9 cm(-1) to the linewidth. Comparison of the 49-Leucine amide I mode and the amide I band of the entire CD3zeta peptide reveals that the vibrational dynamics are not uniform along the length of the peptide. Possible origins for the large amount of inhomogeneity present at the 49-Leucine site are discussed.     

DFT, ab initio and FT-IR studies of the structure of sulfonamide triazenes

We present an in-depth investigation of the structural, infrared spectra and visible spectra of various triazenes with a sulfonamide moiety. The preparation of 4-sulfonyl amide benzenediazonium chloride with cyclic amines of various ring sizes (pyrrolidine, piperidine, 4-methylpiperidine, N-methylpiperazine, morpholine and hexamethyleneimine) are theoretically investigated using the density functional theory (DFT) and Hartree-Fock (HF) levels of theory with the standard 6-31G* basis set. The calculated vibrational frequencies are evaluated via comparison with experimental values. The vibrational spectral data obtained from solid-phase FT-IR spectra are assigned modes based on the results of the theoretical calculations. The UV-Vis spectrum of each compound is measured in various solvents with a wide range of polarity to examine the role of solvent chemistry on observed spectral changes. The calculated spectra are in good agreement with the experimental spectra.     

Two-dimensional infrared spectral signatures of 3(10)- and alpha-helical peptides

Two-dimensional infrared (2D IR) spectra of Calpha-alkylated model octapeptides Z-(Aib)8-OtBu, Z-(Aib)5-L-Leu-(Aib)2-OMe, and Z-[L-(alphaMeVal)]8-OtBu have been measured in the amide I region to acquire 2D spectral signatures characteristic of 3(10)- and alpha-helical conformations. Phase-adjusted 2D absorptive spectra recorded with parallel polarizations are dominated by intense diagonal peaks, whereas 2D rephasing spectra obtained at the double-crossed polarization configuration reveal cross-peak patterns that are essential for structure determination. In CDCl3, all three peptides are of the 3(10)-helix conformation and exhibit a doublet cross-peak pattern. In 1,1,1,3,3,3-hexafluoroisopropanol, Z-[L-(alphaMeVal)]8-OtBu undergoes slow acidolysis and 3(10)-to-alpha-helix transition. In the course of this conformational change, its 2D rephasing spectrum evolves from an elongated doublet, characteristic of a distorted 3(10)-helix, to a multiple-peak pattern, after becoming an alpha-helix. The linear IR and 2D absorptive spectra are much less informative in discerning the structural changes. The experimental spectra are compared to simulations based on a vibrational exciton Hamiltonian model. The through-bond and through-space vibrational couplings are modeled by ab initio coupling maps and transition dipole interactions. The local amide I frequency is evaluated by a new approach that takes into account the effects of hydrogen-bond geometry and sites. The static diagonal and off-diagonal disorders are introduced into the Hamiltonian through statistical models to account for conformational fluctuations and inhomogeneous broadening. The sensitivity of cross-peak patterns to different helical conformations and the chain length dependence of the spectral features for short 3(10)- and alpha-helices are discussed.     

Time-resolved vibrational spectroscopy of a molecular shuttle

Time-resolved vibrational spectroscopy is used to investigate the inter-component motion of an ultraviolet-triggered two-station molecular shuttle. The operation cycle of this molecular shuttle involves several intermediate species, which are observable in the amide I and amide II regions of the mid-IR spectrum. Using ab initio calculations on specific parts of the rotaxane, and by comparing the transient spectra of the normal rotaxane with that of the N-deuterated version, we can assign the observed vibrational modes of each species occurring during the shuttling cycle in an unambiguous way. The complete time- and frequency-dependent data set is analyzed using singular value decomposition (SVD). Using a kinetic model to describe the time-dependent concentrations of the transient species, we derive the absorption spectra associated with each stage in the operation cycle of the molecular shuttle, including the recombination of the charged species.     

Signatures of beta-sheet secondary structures in linear and two-dimensional infrared spectroscopy

Using idealized models for parallel and antiparallel beta sheets, we calculate the linear and two-dimensional infrared spectra of the amide I vibration as a function of size and secondary structure. The model assumes transition-dipole coupling between the amide I oscillators in the sheet and accounts for the anharmonic nature of these oscillators. Using analytical and numerical methods, we show that the nature of the one-quantum vibrational eigenstates, which govern the linear spectrum, is, to a large extent, determined by the symmetry of the system and the relative magnitude of interstrand interactions. We also find that the eigenstates, in particular their trends with system size, depend sensitively on the secondary structure of the sheet. While in practice these differences may be difficult to distinguish in congested linear spectra, we demonstrate that they give rise to promising markers for secondary structure in the two-dimensional spectra. In particular, distinct differences occur between the spectra of parallel and antiparallel beta sheets and between beta hairpins and extended beta sheets.     

Nonadiabatic effects on peptide vibrational dynamics induced by conformational changes

Quantum dynamical simulations of vibrational spectroscopy have been carried out for glycine dipeptide (CH(3)-CO-NH-CH(2)-CO-NH-CH(3)). Conformational structure and dynamics are modeled in terms of the two Ramachandran dihedral angles of the molecular backbone. Potential energy surfaces and harmonic frequencies are obtained from electronic structure calculations at the density functional theory (DFT) [B3LYP/6-31+G(d)] level. The ordering of the energetically most stable isomers (C(7) and C(5)) is reversed upon inclusion of the quantum mechanical zero point vibrational energy. Vibrational spectra of various isomers show distinct differences, mainly in the region of the amide modes, thereby relating conformational structures and vibrational spectra. Conformational dynamics is modeled by propagation of quantum mechanical wave packets. Assuming a directed energy transfer to the torsional degrees of freedom, transitions between the C(7) and C(5) minimum energy structures occur on a sub-picosecond time scale (700...800 fs). Vibrationally nonadiabatic effects are investigated for the case of the coupled, fundamentally excited amide I states. Using a two state-two mode model, the resulting wave packet dynamics is found to be strongly nonadiabatic due to the presence of a seam of the two potential energy surfaces. Initially prepared adiabatic vibrational states decay upon conformational change on a time scale of 200...500 fs with population transfer of more than 50% between the coupled amide I states. Also the vibrational energy transport between localized (excitonic) amide I vibrational states is strongly influenced by torsional dynamics of the molecular backbone where both enhanced and reduced decay rates are found. All these observations should allow the detection of conformational changes by means of time-dependent vibrational spectroscopy.     

Theoretical calculations of infrared absorption, vibrational circular dichroism, and two-dimensional vibrational spectra of acetylproline in liquids water and chloroform

Infrared absorption, vibrational circular dichroism, and two-dimensional infrared pump-probe and photon echo spectra of acetylproline solutions are theoretically calculated and directly compared with experiments. In order to quantitatively determine interpeptide interaction-induced amide I mode frequency shifts, high-level quantum chemistry calculations were performed. The solvatochromic amide I mode frequency shift and fluctuation were taken into account by carrying out molecular dynamics simulations of acetylproline dissolved in liquids water and chloroform and by using the extrapolation method developed recently. We first studied correlation time scales of the two amide I vibrational frequency fluctuations, cross correlation between the two fluctuating local mode frequencies, ensemble averaged conformations of the acetylproline molecule in liquids water and chloroform. The corresponding conformations of the acetylproline in liquids water and chloroform are close to the ideal 3(10) helix and the C(7) structure, respectively. A few methods proposed to determine the angle between the two transition dipoles associated with the amide I vibrations were tested and their limitations are discussed.     

Determination of conformational preferences of dipeptides using vibrational spectroscopy

The NMR coupling constants ((3)J(H(N), H(alpha))) of dipeptides indicate that the backbone conformational preferences vary strikingly among dipeptides. These preferences are similar to those of residues in small peptides, denatured proteins, and the coil regions of native proteins. Detailed characterization of the conformational preferences of dipeptides is therefore of fundamental importance for understanding protein structure and folding. Here, we studied the conformational preferences of 13 dipeptides using infrared and Raman spectroscopy. The main advantage of vibrational spectroscopy over NMR spectroscopy is in its much shorter time scale, which enables the determination of the conformational preferences of short-lived states. Accuracy of structure determination using vibrational spectroscopy depends critically on identification of the vibrational parameters that are sensitive to changes in conformation. We show that the frequencies of the amide I band and the A12 ratio of the amide I components of dipeptides correlate with the (3)J(H(N), H(alpha)). These two infrared vibrational parameters are thus analogous to (3)J(H(N), H(alpha)), indicators for the preference for the dihedral angle phi. We also show that the intensities of the components of the amide III bands in infrared spectra and the intensities of the skeletal vibrations in Raman spectra are indicators of populations of the P(II), beta, and alpha(R) conformations. The results show that alanine dipeptide adopts predominantly a PII conformation. The population of the beta conformation increases in valine dipeptides. The populations of the alpha(R) conformation are generally small. These data are in accord with the electrostatic screening model of conformational preferences.