Peptide conformational heterogeneity revealed from nonlinear vibrational spectroscopy and molecular-dynamics simulations

Simulation Studies of Amide I IR Absorption and Two-Dimensional IR Spectra of β Hairpins in Liquid Water

Amide I IR absorption and two-dimensional (2D) IR photon echo spectra of a model beta hairpin in aqueous solution are theoretically studied and simulated by combining semiempirical quantum chemistry calculations and molecular dynamics simulation methods. The instantaneous normal-mode analysis of the beta hairpin in solution is performed to obtain the density of states and the inverse participation ratios of the one-exciton states. The motional and exchange narrowing processes are taken into account by employing the time-correlation function theory for the linear and nonlinear response functions. Numerically simulated IR absorption and 2D spectra are then found to be determined largely by the amide I normal modes delocalized on the peptides in the two strands. The site-specific isotope-labeling effects on the IR and 2D IR spectra are discussed. The simulation results for the ideal (A17) beta hairpin are directly compared with those of the realistic 16-residue (GB1) beta hairpin from an immunoglobulin G-binding protein. It was found that the characteristic features in IR and 2D spectra of both the ideal (A17) beta hairpin and the GB1 beta hairpin are the same. The simulated IR spectrum of the GB1 beta hairpin is found to be in good agreement with experiment, which demonstrates that the present computational method is quantitatively reliable.

Dissecting coherent vibrational spectra of small proteins into secondary structural elements by sensitivity analysis

The response of proteins to sequences of femtosecond infrared pulses provides a multidimensional view into their equilibrium distribution of structures and snapshot pictures of fast-triggered dynamical events. Analyzing these experiments requires advanced computational tools for assigning regions in the resulting multi-dimensional correlation plots to specific secondary-structure elements and their couplings. A differential sensitivity analysis technique based on a perturbation of the local (real space) Hamiltonian is developed to achieve that goal. Application to the amide I region of a small globular protein reveals regions associated with the alpha-helix, beta-sheet, and their coupling. Comparison of signals generated in different directions shows that the double-quantum-coherence signal has a higher sensitivity to the couplings compared with the single-quantum-coherence (photon echo) technique.

Quantum corrections in vibrational and electronic condensed phase spectroscopy: line shapes and echoes

Various linear and nonlinear vibrational and electronic spectroscopy experiments in liquids are usually analyzed within the second-cumulant approximation, and therefore the fundamental quantity of interest is the equilibrium time-correlation function of the fluctuating transition frequency. In the usual approach the "bath" variables responsible for the fluctuating frequency are treated classically, leading to a classical time-correlation function. Alternatively, sometimes a quantum correction appropriate for relatively high temperatures is included, which adds an imaginary part to the classical time-correlation function. This approach, although appealing, does not satisfy detailed balance. One can consider a similar correction, but where detailed balance is satisfied, by using the harmonic quantum correction factor. In this article, we compare these approaches for a model system and two realistic examples. Our conclusion is that for linear spectroscopy the classical result is usually adequate, whereas for nonlinear spectroscopy it can be more important to include quantum corrections.

Collective Solvent Coordinates for the Infrared Spectrum of HOD in D2O Based on an ab Initio Electrostatic Map

An ab initio MP2 vibrational Hamiltonian of HOD in an external electrostatic potential parametrized by the electric field and its gradient-tensor is constructed. By combining it with the fluctuating electric field induced by the D(2)O solvent obtained from molecular dynamics simulations, we calculate the infrared absorption of the O-H stretch. The resulting solvent shift and infrared line shape for three force fields (TIP4P, SPC/E, and SW) are in good agreement with the experiment. A collective coordinate response for the solvent effect is constructed by identifying the main electrostatic field and gradient components contributing to the line shape. This allows a realistic stochastic Liouville equation simulation of the line shapes which is not restricted to Gaussian frequency fluctuations.

Stochastic Liouville equation simulation of multidimensional vibrational line shapes of trialanine

The line shapes detected in coherent femtosecond vibrational spectroscopies contain direct signatures of peptide conformational fluctuations through their effect on vibrational frequencies and intermode couplings. These effects are simulated in trialanine using a Green's function solution of a stochastic Liouville equation constructed for four collective bath coordinates (two Ramachandran angles affecting the mode couplings and two diagonal energies). We find that fluctuations of the Ramachandran angles which hardly affect the linear absorption can be effectively probed by two-dimensional spectra. The signal generated at k(1)+k(2)-k(3) is particularly sensitive to such fluctuations.

Ultrafast vibrational spectroscopy of water and aqueous N-methylacetamide: Comparison of different electronic structure/molecular dynamics approaches

Kwac and Cho [J. Chem. Phys. 119, 2247 (2003)] have recently developed a combined electronic structure/molecular dynamics approach to vibrational spectroscopy in liquids. The method involves fitting ab initio vibrational frequencies for a solute in a cluster of solvent molecules to a linear combination of the electrostatic potentials on the solute atoms due to the charges on the solvent molecules. These authors applied their method to the N-methylacetamide-D/D(2)O system. We (S. A. Corcelli, C. P. Lawrence, and J. L. Skinner, [J. Chem. Phys. 120, 8107 (2004)]) have recently explored a closely related method, where instead of the electrostatic potential, the solute vibrational frequencies are fit to the components of the electric fields on the solute atoms due to the solvent molecules. We applied our method to the HOD/D(2)O and HOD/H(2)O systems. In order to make a direct comparison of these two approaches, in this paper we apply their method to the water system, and our method to the N-methylacetamide system. For the water system we find that the electric field method is superior to the potential approach, as judged by comparison with experiments for the absorption line shape. For the N-methylacetamide system the two methods are comparable.

Transient two-dimensional infrared spectroscopy: Exploring the polarization dependence

We present a general expression for the polarization dependence of transient two-dimensional IR spectroscopy (T2D-IR), a technique designed to measure 2D-IR spectra of transient species. T2D-IR is a UV pump narrowband-IR-pump broadband-IR-probe experiment of fifth order in the laser field which involves up to three different transition dipole moments. The UV pulse adds an additional degree of freedom in polarization as compared to 2D-IR spectroscopy and increases the versatility of signal manipulation and the potential structural information content of the signals. The polarization conditions leading to a maximum of structural information are discussed. Important special cases of polarization conditions are formulated. The application of polarization selectivity is demonstrated for different types of T2D-IR experiments on photo triggered metal-to-ligand charge transfer in the model system [Re(CO)(3)(dmbpy)Cl].

Double-resonance versus pulsed Fourier transform two-dimensional infrared spectroscopy: An experimental and theoretical comparison

In this study we focus on the differences and analogies of two experimental implementations of two-dimensional infrared (2D-IR) spectroscopy: double-resonance or dynamic hole burning 2D-IR spectroscopy and pulsed Fourier transform or heterodyne detected photon echo spectroscopy. A comparison is done theoretically as well as experimentally by contrasting data obtained from both methods. As an example we have studied the strongly coupled asymmetric and symmetric carbonyl stretching vibrations of dicarbonylacetylacetonato rhodium dissolved in hexane. Both methods yield the same peaks in a 2D-IR spectrum. Within certain approximations we derive an analytic expression which shows that the 2D-IR spectra are broadened in one frequency dimension in the double-resonance experiment by convolution with the pump pulse spectral width, while the spectral resolution in the other frequency direction is the same in both cases.

A two-dimensional infrared study of localization, structure, and dynamics of a dipeptide in membrane environment

Using methods of time-resolved two-dimensional infrared (2D-IR) spectroscopy, we approach the problem of the structural characterization of small polypeptide systems in a membrane environment. The 2D-IR spectra recorded for a model dipeptide in different environments demonstrated a significant change in the homogeneous and inhomogeneous broadenings of the amide I resonances when the molecule inserts either into a surfactant or a phospholipid membrane. Besides the change in the diagonal features in the 2D-IR response, we observe both intramolecular and intermolecular crosspeaks between the carbonyls of the dipeptide and the phospholipid. Considering the character of the diagonal peaks and the presence of the crosspeaks, we discuss the localization of the dipeptide moieties in the membrane. Using both the anisotropy and relative intensity of the observed intramolecular crosspeaks between the two amide I modes, we provide observables that help to determine the phi/psi-dihedral angles for the backbone of the dipeptide. Time dependent studies revealed slower conformational fluctuations of the dipeptide backbone in a membrane as compared to that in an aqueous environment.

Combined electronic structure/molecular dynamics approach for ultrafast infrared spectroscopy of dilute HOD in liquid H2O and D2O

We present a new approach that combines electronic structure methods and molecular dynamics simulations to investigate the infrared spectroscopy of condensed phase systems. This approach is applied to the OH stretch band of dilute HOD in liquid D2O and the OD stretch band of dilute HOD in liquid H2O for two commonly employed models of water, TIP4P and SPC/E. Ab initio OH and OD anharmonic transition frequencies are calculated for 100 HOD x (D2O)n and HOD x(H2O)n (n = 4-9) clusters randomly selected from liquid water simulations. A linear empirical relationship between the ab initio frequencies and the component of the electric field from the solvent along the bond of interest is developed. This relationship is used in a molecular dynamics simulation to compute frequency fluctuation time-correlation functions and infrared absorption line shapes. The normalized frequency fluctuation time-correlation functions are in good agreement with the results of previous theoretical approaches. Their long-time decay times are 0.5 ps for the TIP4P model and 0.9 ps for the SPC/E model, both of which appear to be somewhat too fast compared to recent experiments. The calculated line shapes are in good agreement with experiment, and improve upon the results of previous theoretical approaches. The methods presented are simple, and transferable to more complicated systems.

Signatures of β-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.

The origin of vibrational mode couplings in various secondary structural motifs of polypeptides

Electrostatic (through-space) and covalent (through-bond) contributions to couplings involving the C[double bond]O and C[bond]N vibrational stretching modes of the amide group in the alpha-helix and the parallel and antiparallel beta-sheet structures of alanine polypeptides are analyzed. Coupling constants computed at the density functional theory level are compared with the transition dipole coupling model and the complete electrostatic interaction between transition densities. We find that the transition densities of C[double bond]O modes are localized, and the electrostatic mechanism then holds. In contrast, the C[bond]N mode transition densities are delocalized, and covalent contributions to the coupling are significant.

Two-dimensional Raman and infrared vibrational spectroscopy for a harmonic oscillator system nonlinearly coupled with a colored noise bath

Multidimensional vibrational response functions of a harmonic oscillator are reconsidered by assuming nonlinear system-bath couplings. In addition to a standard linear-linear (LL) system-bath interaction, we consider a square-linear (SL) interaction. The LL interaction causes the vibrational energy relaxation, while the SL interaction is mainly responsible for the vibrational phase relaxation. The dynamics of the relevant system are investigated by the numerical integration of the Gaussian-Markovian Fokker-Planck equation under the condition of strong couplings with a colored noise bath, where the conventional perturbative approach cannot be applied. The response functions for the fifth-order nonresonant Raman and the third-order infrared (or equivalently the second-order infrared and the seventh-order nonresonant Raman) spectra are calculated under the various combinations of the LL and the SL coupling strengths. Calculated two-dimensional response functions demonstrate that those spectroscopic techniques are very sensitive to the mechanism of the system-bath couplings and the correlation time of the bath fluctuation. We discuss the primary optical transition pathways involved to elucidate the corresponding spectroscopic features and to relate them to the microscopic sources of the vibrational nonlinearity induced by the system-bath interactions. Optical pathways for the fifth-order Raman spectroscopies from an "anisotropic" medium were newly found in this study, which were not predicted by the weak system-bath coupling theory or the standard Brownian harmonic oscillator model.

Investigations into sequence and conformational dependence of backbone entropy, inter-basin dynamics and the Flory isolated-pair hypothesis for peptides

The populations and transitions between Ramachandran basins are studied for combinations of the standard 20 amino acids in monomers, dimers and trimers using an implicit solvent Langevin dynamics algorithm and employing seven commonly used force-fields. Both the basin populations and inter-conversion rates are influenced by the nearest neighbor's conformation and identity, contrary to the Flory isolated-pair hypothesis. This conclusion is robust to the choice of force-field, even though the use of different force-fields produces large variations in the populations and inter-conversion rates between the dominant helical, extended beta, and polyproline II basins. The computed variation of conformational and dynamical properties with different force-fields exceeds the difference between explicit and implicit solvent calculations using the same force-field. For all force-fields, the inter-basin transitions exhibit a directional dependence, with most transitions going through extended beta conformation, even when it is the least populated basin. The implications of these results are discussed in the context of estimates for the backbone entropy of single residues, and for the ability of all-atom simulations to reproduce experimental protein folding data.