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

Study of the Myosin Relay Helix Peptide by Molecular Dynamics Simulations, Pump-Probe and 2D Infrared Spectroscopy

The Davydov model was conjectured to describe how an amide I excitation created during ATP hydrolysis in myosin might be significant in providing energy to drive myosin's chemomechanical cycle. The free energy surfaces of the myosin relay helix peptide dissolved in 2,2,2-trifluoroethanol (TFE), determined by metadynamics simulations, demonstrate local minima differing in free energy by only ~2 kT, corresponding to broken and stabilized hydrogen bonds, respectively. Experimental pump-probe and 2D infrared spectroscopy were performed on the peptide dissolved in TFE. The relative heights of two peaks seen in the pump-probe data and the corresponding relative volumes of diagonal peaks seen in the 2D-IR spectra at time delays between 0.5 ps and 1 ps differ noticeably from what is seen at earlier or later time delays or in the linear spectrum, indicating that a vibrational excitation may influence the conformational state of this helix. Thus, it is possible that the presence of an amide I excitation may be a direct factor in the conformational state taken on by the myosin relay helix following ATP hydrolysis in myosin.

Theoretical basis for interpreting heterodyne chirality-selective sum frequency generation spectra of water

Chirality-selective vibrational sum frequency generation (chiral SFG) spectroscopy has emerged as a powerful technique for the study of biomolecular hydration water due to its sensitivity to the induced chirality of the first hydration shell. Thus far, water O-H vibrational bands in phase-resolved heterodyne chiral SFG spectra have been fit using one Lorentzian function per vibrational band, and the resulting fit has been used to infer the underlying frequency distribution. Here, we show that this approach may not correctly reveal the structure and dynamics of hydration water. Our analysis illustrates that the chiral SFG responses of symmetric and asymmetric O-H stretch modes of water have opposite phase and equal magnitude and are separated in energy by intramolecular vibrational coupling and a heterogeneous environment. The sum of the symmetric and asymmetric responses implies that an O-H stretch in a heterodyne chiral SFG spectrum should appear as two peaks with opposite phase and equal amplitude. Using pairs of Lorentzian functions to fit water O-H stretch vibrational bands, we improve spectral fitting of previously acquired experimental spectra of model β-sheet proteins and reduce the number of free parameters. The fitting allows us to estimate the vibrational frequency distribution and thus reveals the molecular interactions of water in hydration shells of biomolecules directly from chiral SFG spectra.

Refining protein amide I spectrum simulations with simple yet effective electrostatic models for local wavenumbers and dipole derivative magnitudes

Analysis of the amide I band of proteins is probably the most wide-spread application of bioanalytical infrared spectroscopy. Although highly desirable for a more detailed structural interpretation, a quantitative description of this absorption band is still difficult. This work optimized several electrostatic models with the aim to reproduce the effect of the protein environment on the intrinsic wavenumber of a local amide I oscillator. We considered the main secondary structures - α-helices, parallel and antiparallel β-sheets - with a maximum of 21 amide groups. The models were based on the electric potential and/or the electric field component along the CO bond at up to four atoms in an amide group. They were bench-marked by comparison to Hessian matrices reconstructed from density functional theory calculations at the BPW91, 6-31G** level. The performance of the electrostatic models depended on the charge set used to calculate the electric field and potential. Gromos and DSSP charge sets, used in common force fields, were not optimal for the better performing models. A good compromise between performance and the stability of model parameters was achieved by a model that considered the electric field at the positions of the oxygen, nitrogen, and hydrogen atoms of the considered amide group. The model describes also some aspects of the local conformation effect and performs similar on its own as in combination with an explicit implementation of the local conformation effect. It is better than a combination of a local hydrogen bonding model with the local conformation effect. Even though the short-range hydrogen bonding model performs worse, it captures important aspects of the local wavenumber sensitivity to the molecular surroundings. We improved also the description of the coupling between local amide I oscillators by developing an electrostatic model for the dependency of the dipole derivative magnitude on the protein environment.

Unpolarized laser method for infrared spectrum calculation of amide I CO bonds in proteins using molecular dynamics simulation

The investigation of the strong infrared (IR)-active amide I modes of peptides and proteins has received considerable attention because a wealth of detailed information on hydrogen bonding, dipole-dipole interactions, and the conformations of the peptide backbone can be derived from the amide I bands. The interpretation of experimental spectra typically requires substantial theoretical support, such as direct ab-initio molecular dynamics simulation or mixed quantum-classical description. However, considering the difficulties associated with these theoretical methods and their applications are limited in small peptides, it is highly desirable to develop a simple yet efficient approach for simulating the amide I modes of any large proteins in solution. In this work, we proposed a comprehensive computational method that extends the well-established molecular dynamics (MD) simulation method to include an unpolarized IR laser for exciting the CO bonds of proteins. We showed the amide I frequency corresponding to the frequency of the laser pulse which resonated with the CO bond vibration. At this frequency, the protein energy and the CO bond length fluctuation were maximized. Overall, the amide I bands of various single proteins and amyloids agreed well with experimental data. The method has been implemented into the AMBER simulation package, making it widely available to the scientific community. Additionally, the application of the method to simulate the transient amide I bands of amyloid fibrils during the IR laser-induced disassembly process was discussed in details.

Ultrafast Two-Dimensional Infrared Spectroscopy Resolved a Structured Lysine 159 on the Cytoplasmic Surface of the Microbial Photoreceptor Bacteriorhodopsin

Bacteriorhodopsin (bR) is a light-driven microbial receptor, and lysine 159 (K159) is a charged residue on the cytoplasmic (CP) side of its E-F loop. However, its conformation and function remain unknown due to fast surface dynamics. By utilizing a ¹³C, ¹⁵N-labeled lysine (K) as an isotope probe, we created a network of site-specific amide-I vibrational signatures (backbone carbonyl stretch) to identify the frequency contribution of the labeled residues to the amide-I excitonic band structure. Thus, the red-shifted amide-I frequency in the ¹³C, ¹⁵N-lysine-labeled bR (uK-bR) to the unlabeled bR (WT-bR) could be differentiated and examined by ultrafast two-dimensional vibrational echo infrared (2D IR) spectroscopy. Our results showed that the backbone carbonyl of K159 is located at a high frequency of ca. 1693 cm^-1 and has a vibrational excited-state relaxation time shorter than the bulk helical amide-I mode at the same frequency, suggesting that K159 may possess a hydrogen-bonded γ-turn structure with E161, one of the carboxylate residues on the CP surface of bR. The 2D solid-state NMR study of uK-bR also revealed conformational dependent lysine residues, from which K159 was found to involve the turn motif. This γ-turn structure maintained by K159 may help to stabilize the E-F loop and support E161 in attracting protons from the bulk during the late stage of the bR photocycle. The combined spectroscopic approach illustrated in this work may be applied to map residue-specific local structures and dynamics of other receptors and large proteins.

Does liquid-liquid phase separation drive peptide folding?

Proline-arginine (PR) dipeptide repeats have been shown to undergo liquid-liquid phase separation and are an example of a growing number of intrinsically disordered proteins that can assemble into membraneless organelles. These structures have been posited as nucleation sites for pathogenic protein aggregation. As such, a better understanding of the effects that the increased local concentration and volumetric crowding within droplets have on peptide secondary structure is necessary. Herein we use Fourier transform infrared (FTIR) and two-dimensional infrared (2DIR) spectroscopy to show that formation of droplets by PR20 accompanies changes in the amide-I spectra consistent with folding into poly-proline helical structures.

The Amide I Spectrum of Proteins-Optimization of Transition Dipole Coupling Parameters Using Density Functional Theory Calculations

The amide I region of the infrared spectrum is related to the protein backbone conformation and can provide important structural information. However, the interpretation of the experimental results is hampered because the theoretical description of the amide I spectrum is still under development. Quantum mechanical calculations, for example, using density functional theory (DFT), can be used to study the amide I spectrum of small systems, but the high computational cost makes them inapplicable to proteins. Other approaches that solve the eigenvalues of the coupled amide I oscillator system are used instead. An important interaction to be considered is transition dipole coupling (TDC). Its calculation depends on the parameters of the transition dipole moment. This work aims to find the optimal parameters for TDC in three major secondary structures: α-helices, antiparallel β-sheets, and parallel β-sheets. The parameters were suggested through a comparison between DFT and TDC calculations. The comparison showed a good agreement for the spectral shape and for the wavenumbers of the normal modes for all secondary structures. The matching between the two methods improved when hydrogen bonding to the amide oxygen was considered. Optimal parameters for individual secondary structures were also suggested.

Insight into the internal structure of amyloid-β oligomers by isotope-edited Fourier transform infrared spectroscopy

Isotope-edited infrared spectroscopy reveals the structural unit of amyloid-β oligomers.

Not All β-Sheets Are the Same: Amyloid Infrared Spectra, Transition Dipole Strengths, and Couplings Investigated by 2D IR Spectroscopy

We report the transition dipole strengths and frequencies of the amyloid β-sheet amide I mode for the aggregated proteins amyloid-β_1-40, calcitonin, α-synuclein, and glucagon. According to standard vibrational coupling models for proteins, the frequencies of canonical β-sheets are set by their size and structural and environmental disorder, which determines the delocalization length of the vibrational excitons. The larger the delocalization the lower the frequency of the main infrared-allowed transition, A_⊥. The models also predict an accompanying increase in transition dipole strength. For the proteins measured here, we find no correlation between transition dipole strengths and amyloid β-sheet transition frequency. To understand this observation, we have extracted from the protein data bank crystal structures of amyloid peptides from which we calculate the amide I vibrational couplings, and we use these in a model β-sheet Hamiltonian to simulate amyloid vibrational spectra. We find that the variations in amyloid β-sheet structures (e.g., dihedral angles, interstrand distances, and orientations) create significant differences in the average values for interstrand and nearest neighbor couplings, and that those variations encompass the variation in measured A_⊥ frequencies. We also find that off-diagonal disorder about the average values explains the range of transition dipole strengths observed experimentally. Thus, we conclude that the lack of correlation between transition dipole-strength and frequency is caused by variations in amyloid β-sheet structure. Taken together, these results indicate that the amide I frequency is very sensitive to amyloid β-sheet structure, the β-sheets of these 4 proteins are not identical, and the assumption that frequency of amyloids scales with β-sheet size cannot be adopted without an accompanying measurement of transition dipole strengths.

Effective representation of amide III, II, I, and A modes on local vibrational modes: Analysis of ab initio quantum calculation results

The Hamiltonian matrix for the first excited vibrational states of a protein can be effectively represented by local vibrational modes constituting amide III, II, I, and A modes to simulate various vibrational spectra. Methods for obtaining the Hamiltonian matrix from ab initio quantum calculation results are discussed, where the methods consist of three steps: selection of local vibrational mode coordinates, calculation of a reduced Hessian matrix, and extraction of the Hamiltonian matrix from the Hessian matrix. We introduce several methods for each step. The methods were assessed based on the density functional theory calculation results of 24 oligopeptides with four different peptide lengths and six different secondary structures. The completeness of a Hamiltonian matrix represented in the reduced local mode space is improved by adopting a specific atom group for each amide mode and reducing the effect of ignored local modes. The calculation results are also compared to previous models using C=O stretching vibration and transition dipole couplings. We found that local electric transition dipole moments of the amide modes are mainly bound on the local peptide planes. Their direction and magnitude are well conserved except amide A modes, which show large variation. Contrary to amide I modes, the vibrational coupling constants of amide III, II, and A modes obtained by analysis of a dipeptide are not transferable to oligopeptides with the same secondary conformation because coupling constants are affected by the surrounding atomic environment.

Transition Dipoles from 1D and 2D Infrared Spectroscopy Help Reveal the Secondary Structures of Proteins: Application to Amyloids

Transition dipoles are an underutilized quantity for probing molecular structures. The transition dipole strengths in an extended system like a protein are modulated by the couplings and thus probe the structures. Here we measure the absolute transition dipole strengths of human and rat amylin in their solution, aggregated, membrane, and micelleular bound forms, using a combination of 1D and 2D infrared spectroscopy. We find that the vibrational modes of amyloid fibers made of human amylin can extend across as many as 12 amino acids, reflecting very ordered β-sheets in the most carefully prepared samples. Rat amylin has FTIR spectra that are nearly identical in solution, micelles, and membranes. We show that the transition dipoles of rat amylin are much larger when bound to micelles and membranes than when in solution, consistent with rat amylin adopting an α-helical structure. We interpret the transition dipole strengths as experimental measurements of the inverse participation ratio often calculated in theoretical studies. The structure of aggregating and membrane-bound proteins can be difficult to identify with existing techniques, especially during kinetics. These results demonstrate how absolute transition dipoles measured via our 1D/2D spectroscopy method can provide important structural information.

Two-Dimensional Vibrational Spectroscopy of a Dissipative System with the Optimized Mean-Trajectory Approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions representing radiation-matter interactions. Here we apply this method to an anharmonic chromophore coupled to a harmonic bath. A forward-backward trajectory implementation of the OMT method is described that addresses the numerical challenges of applying the OMT to large systems with disparate frequency scales. The OMT is shown to well reproduce line shapes and waiting time dynamics in the pure dephasing limit of weak coupling to an off-resonant bath. The OMT is also shown to describe a case where energy transfer is the predominant source of line broadening.

Vibrational coherence and energy transfer in two-dimensional spectra with the optimized mean-trajectory approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions that represent radiation-matter interactions. Here, we extend the OMT to include additional vibrational coherence and energy transfer processes. This generalized approximation is applied to a pair of anharmonic chromophores coupled to a bath. The resulting 2D spectra are shown to reflect coherence transfer between normal modes.

The non-uniform early structural response of globular proteins to cold denaturing conditions: a case study with Yfh1

The mechanism of cold denaturation in proteins is often incompletely understood due to limitations in accessing the denatured states at extremely low temperatures. Using atomistic molecular dynamics simulations, we have compared early (nanosecond timescale) structural and solvation properties of yeast frataxin (Yfh1) at its temperature of maximum stability, 292 K (Ts), and the experimentally observed temperature of complete unfolding, 268 K (Tc). Within the simulated timescales, discernible "global" level structural loss at Tc is correlated with a distinct increase in surface hydration. However, the hydration and the unfolding events do not occur uniformly over the entire protein surface, but are sensitive to local structural propensity and hydrophobicity. Calculated infrared absorption spectra in the amide-I region of the whole protein show a distinct red shift at Tc in comparison to Ts. Domain specific calculations of IR spectra indicate that the red shift primarily arises from the beta strands. This is commensurate with a marked increase in solvent accessible surface area per residue for the beta-sheets at Tc. Detailed analyses of structure and dynamics of hydration water around the hydrophobic residues of the beta-sheets show a more bulk water like behavior at Tc due to preferential disruption of the hydrophobic effects around these domains. Our results indicate that in this protein, the surface exposed beta-sheet domains are more susceptible to cold denaturing conditions, in qualitative agreement with solution NMR experimental results.

Discriminating trpzip2 and trpzip4 peptides' folding landscape using the two-dimensional infrared spectroscopy: a simulation study

We analyzed, based on the theoretical spectroscopic modeling, how the differences in the folding landscapes of two β-hairpin peptides trpzip2 and trpzip4 are reflected in their thermal unfolding infrared measurements. The isotope-edited equilibrium FTIR and two dimensional infrared spectra of the two peptides were calculated, using the nonlinear exciton propagation method, at a series of temperatures. The spectra calculations were based on the configuration distributions generated using the GB(OBC) implicit solvent MD simulation and the integrated tempering sampling technique. Conformational analysis revealed the different local thermal stabilities for these two peptides, which suggested the different folding landscapes. Our study further suggested that the ellipticities of the isotope peaks in the coherent IR signals are more sensitive to these local stability differences compared with other spectral features such as the peak intensities. Our technique can thus be combined with the relevant experimental measurements to achieve a better understanding of the peptide folding behaviors.

Amide-I characteristics of helical β-peptides by linear infrared measurement and computations

In this work, we have examined the amide-I characteristics of three β-peptide oligomers in typical helical conformations (two in 14-helix and one in 12/10-helix), solvated in water, methanol, and chloroform, respectively. Local-mode frequencies and their distributions were computed using a molecular-mechanics force field based frequency map that was constructed on the basis of molecular dynamics simulations. The local-mode frequencies were found to be determined primarily by peptide backbone and side chain, rather by solvent, suggesting their local structural sensitivities. Intermode vibrational couplings computed using a transition dipole scheme were found to be very sensitive to peptide conformation, with their signs and magnitudes varying periodically along the peptide chain. Linear infrared absorption spectra of the three peptides, simulated using a frequency-frequency time-correlation function method, were found to be in fair agreement with experimental results. Normalized potential energy distribution analysis indicated that the amide-I mode can delocalize over a few amide units. However, the IR band structure appears to be more sophisticated in helical β-peptides than in helical α-peptides.

Development and validation of transferable amide I vibrational frequency maps for peptides

Infrared (IR) spectroscopy of the amide I band has been widely utilized for the analysis of peptides and proteins. Theoretical modeling of IR spectra of proteins requires an accurate and efficient description of the amide I frequencies. In this paper, amide I frequency maps for protein backbone and side chain groups are developed from experimental spectra and vibrational lifetimes of N-methylacetamide and acetamide in different solvents. The frequency maps, along with established nearest-neighbor frequency shift and coupling schemes, are then applied to a variety of peptides in aqueous solution and reproduce experimental spectra well. The frequency maps are designed to be transferable to different environments; therefore, they can be used for heterogeneous systems, such as membrane proteins.

Simulation of the amide I absorption of stacked β-sheets

Aggregated β-sheet structures are associated with amyloid and prion diseases. Techniques capable of revealing detailed structural and dynamical information on β-sheet structure are thus of great biomedical and biophysical interest. In this work, the infrared (IR) amide I spectral characteristics of stacked β-sheets were modeled using the transition dipole coupling model. For a test set of β-sheet stacks, the simulated amide I spectrum was analyzed with respect to the following parameters; intersheet distance, relative rotation of the sheets with respect to each other and the effect of number of sheets stacked. The amide I maximum shifts about 5 cm(-1) to higher wavenumbers when the intersheet distance between two identical β-sheets decreases from 20 to 5 Å. Rotation around the normal of one of the sheets relative to the other results in maximum intersheet coupling near 0° and 180°. Upon of rotation from 0° to 90° at an intersheet distance of 9 Å, the amide I maximum shifts about 3 cm(-1). Tilting of one of the sheets by 30° from the normal results in a shift of the amide I maximum by less than 1 cm(-1). When stacking several β-sheets along the normal, the amide I maximum shifts to higher wavenumbers with increasing stack size. The amide I maximum shifts about 6 cm(-1) when stacking four sheets with an intersheet distance of 9 Å. The study provides an aid in the interpretation of the IR amide I region for experiments involving β-sheets and creates awareness of the many effects that determine the spectrum of β-sheet structures.

Probing the folding of mini-protein Beta3s by two-dimensional infrared spectroscopy; simulation study

We propose to use infrared coherent two-dimensional correlation spectroscopy (2DCS) to characterize the folding mechanism of the mini-protein Beta3s. In this study Beta3s was folded by molecular dynamics (MD) simulation and intermediate conformational ensembles were identified. The one and two-dimensional correlation spectrum was calculated for the intermediate and native states of the mini-protein. A direct structure-spectra relationship was determined by analysis of conformational properties and specific residue contributions. We identified the structural origin of diagonal and off-diagonal peaks in the 2DCS spectra for the native and intermediate conformational ensembles in the folding mechanism. This work supports the implementation of computational techniques in conjunction with experimental 2DCS to study the folding mechanism of proteins. In addition to exploring the folding mechanism the work presented here can be applied in combination with experiment to refine and validate current molecular dynamics force fields.

PACS Codes: 87.15.Cc, 87.15.hm, 87.15.hp

Integrated and dispersed photon echo studies of nitrile stretching vibration of 4-cyanophenol in methanol

By means of integrated and dispersed IR photon echo measurement methods, the vibrational dynamics of C-N stretch modes in 4-cyanophenol and 4-cyanophenoxide in methanol is investigated. The vibrational frequency-frequency correlation function (FFCF) is retrieved from the integrated photon echo signals by assuming that the FFCF is described by two exponential functions with about 400 fs and a few picosecond components. The excited state lifetimes of the C-N stretch modes of neutral and anionic 4-cyanophenols are 1.45 and 0.91 ps, respectively, and the overtone anharmonic frequency shifts are 25 and 28 cm(-1). At short waiting times, a notable underdamped oscillation, which is attributed to a low-frequency intramolecular vibration coupled to the CN stretch, in the integrated and dispersed vibrational echo as well as transient grating signals was observed. The spectral bandwidths of IR absorption and dispersed vibrational echo spectra of the 4-cyanophenoxide are significantly larger than those of its neutral form, indicating that the strong interaction between phenoxide and methanol causes large frequency fluctuation and rapid population relaxation. The resonance effects in a paradisubstituted aromatic compound would be of interest in understanding the conjugation effects and their influences on chemical reactivity of various aromatic compounds in organic solvents.

Deconstructing the native state: energy landscapes, function, and dynamics of globular proteins

Proteins are highly complex molecules with features exquisitely selected by nature to carry out essential biological functions. Physical chemistry and polymer physics provide us with the tools needed to make sense of this complexity. Upon translation, many proteins fold to a thermodynamically stable form known as the native state. The native state is not static, but consists of a hierarchy of conformations, that are continuously explored through dynamics. In this review we provide a brief introduction to some of the core concepts required in the discussion of the protein native dynamics using energy landscapes ideas. We first discuss recent works which have challenged the structure-function paradigm by demonstrating function in disordered proteins. Next we examine the hierarchical organization in the energy landscapes using atomistic molecular dynamics simulations and principal component analysis. In particular, the role of direct and water-mediated contacts in sculpting the landscape is elaborated. Another approach to studying the native state ensemble is based on choosing high-resolution order parameters for computing one- or two-dimensional free energy surfaces. We demonstrate that 2D free energy surfaces provide rich thermodynamic and kinetic information about the native state ensemble. Brownian dynamics simulations on such a surface indicate that protein conformational dynamics is weakly activated. Finally, we briefly discuss implicit and coarse-grained protein models and emphasize the solvent role in determining native state structure and dynamics.

Vibrational solvatochromism and electrochromism: coarse-grained models and their relationships

A theoretical description of vibrational solvatochromism and electrochromism is presented by using a coarse-grained model based on a distributed charge and multipole interaction theory. Solvatochromic frequency shift has been described by considering the interaction between distributed charges of a solute and electrostatic potential due to distributed charges of solvent molecules. Another approach was based on the expansion of the solvatochromic frequency shift in terms of solvent electric field and its gradient at distributed sites on solute. The relationship between these two approaches is elucidated and their validities are discussed. It is also shown that the distributed charge and multipole model for solvatochromism developed here can be used to describe vibrational Stark effects on frequency and transition dipole moment. The relationship between the vibrational Stark tuning rate and the parameters obtained from recent vibrational solvatochromism studies is clarified and used to determine the vibrational Stark tuning rates of a few stretching modes, which are then directly compared with experimentally measured values. We anticipate that the present theoretical model can be used to study a variety of vibrational solvatochromic and electrochromic phenomena and to extract critical information on local electrostatic environment around a small IR probe in solution or protein from linear and nonlinear IR spectroscopic studies.

Amide I Raman optical activity of polypeptides: fragment approximation

Vibrational optical activity (VOA) is an important property used to determine the absolute configuration of a chiral molecule in condensed phases. In particular, vibrational circular dichroism and Raman optical activity (ROA) are two representative VOA measurement techniques that have been extensively used to study structures and dynamics of biomolecules. Recently, the amide I vibrational circular dichroism of polypeptides was theoretically described by using fragment approximation methods, which are based on the assumption that amide I VOA can be described as a linear combination of those of constituent fragment peptide units. Here, we develop a fragment approximation theory applicable to numerical simulations of Raman and Raman optical activity spectra for the amide I vibrations in polypeptides. For an alanine dipeptide and pentapeptide analogs, we carried out density functional theory calculations of polarizability, magnetic dipole-, and electric quadrupole-ROA tensors. Numerically simulated spectra using the fragment approximation are directly compared to density functional theory results. Furthermore, the simulated ROA spectra of alanine-based right-handed alpha-helix and polyproline II polypeptides are directly compared to the previously reported experimental results. The agreements were found to be excellent, which suggests that the fragment approximation method developed for the numerical simulation of ROA spectrum of polypeptide in solution is valid and useful.

Empirical amide I vibrational frequency map: application to 2D-IR line shapes for isotope-edited membrane peptide bundles

The amide I vibrational mode, primarily associated with peptide-bond carbonyl stretches, has long been used to probe the structures and dynamics of peptides and proteins by infrared (IR) spectroscopy. A number of ab initio-based amide I vibrational frequency maps have been developed for calculating IR line shapes. In this paper, a new empirical amide I vibrational frequency map is developed. To evaluate its performance, we applied this map to a system of isotope-edited CD3-zeta membrane peptide bundles in aqueous solution. The calculated 2D-IR diagonal line widths vary from residue to residue and show an asymmetric pattern as a function of position in the membrane. The theoretical results are in fair agreement with experiments on the same system. Through analysis of the computed frequency time-correlation functions, it is found that the 2D-IR diagonal widths are dominated by contributions from the inhomogeneous frequency distributions, from which it follows that these widths are a good probe of the extent of local structural fluctuations. Thus, the asymmetric pattern of line widths follows from the asymmetric structure of the bundle in the membrane.