General treatment of vibrations of helical molecules and application to transition dipole coupling in amide I and amide II modes of α-helical poly(l-alanine)

Broadband Multidimensional Spectroscopy Identifies the Amide II Vibrations in Silkworm Films

We used two-dimensional infrared spectroscopy to disentangle the broad infrared band in the amide II vibrational regions of Bombyx mori native silk films, identifying the single amide II modes and correlating them to specific secondary structure. Amide I and amide II modes have a strong vibrational coupling, which manifests as cross-peaks in 2D infrared spectra with frequencies determined by both the amide I and amide II frequencies of the same secondary structure. By cross referencing with well-known amide I assignments, we determined that the amide II (N-H) absorbs at around 1552 and at 1530 cm-1 for helical and β-sheet structures, respectively. We also observed a peak at 1517 cm-1 that could not be easily assigned to an amide II mode, and instead we tentatively assigned it to a Tyrosine sidechain. These results stand in contrast with previous findings from linear infrared spectroscopy, highlighting the ability of multidimensional spectroscopy for untangling convoluted spectra, and suggesting the need for caution when assigning silk amide II spectra.

Cross-α/β polymorphism of PSMα3 fibrils

Phenol soluble modulins (PSMs) are an important class of peptides secreted by Staphylococcus aureus bacteria. The toxicity to human cells and unique ability of one such peptide, PSMα3, to aggregate into an α-helical amyloid-like structure may hold a key to a better understanding of the virulence of dangerous pathogens such as methicillin resistant S. aureus. In reporting a detailed two-dimensional infrared (2DIR) analysis of PSMα3, we found direct evidence of multiple aggregate architectures existing in equilibrium with one another. We also discovered a unique and characteristic 2DIR spectroscopic signature that unambiguously reports on the presence of the unusual and highly cytotoxic cross-α amyloid structure.

The formation of ordered cross-β amyloid protein aggregates is associated with a variety of human disorders. While conventional infrared methods serve as sensitive reporters of the presence of these amyloids, the recently discovered amyloid secondary structure of cross-α fibrils presents new questions and challenges. Herein, we report results using Fourier transform infrared spectroscopy and two-dimensional infrared spectroscopy to monitor the aggregation of one such cross-α–forming peptide, phenol soluble modulin alpha 3 (PSMα3). Phenol soluble modulins (PSMs) are involved in the formation and stabilization of Staphylococcus aureus biofilms, making sensitive methods of detecting and characterizing these fibrils a pressing need. Our experimental data coupled with spectroscopic simulations reveals the simultaneous presence of cross-α and cross-β polymorphs within samples of PSMα3 fibrils. We also report a new spectroscopic feature indicative of cross-α fibrils.

Quantification of transition dipole strengths using 1D and 2D spectroscopy for the identification of molecular structures via exciton delocalization: application to α-helices

Vibrational and electronic transition dipole strengths are often good probes of molecular structures, especially in excitonically coupled systems of chromophores. One cannot determine transition dipole strengths using linear spectroscopy unless the concentration is known, which in many cases it is not. In this paper, we report a simple method for measuring transition dipole moments from linear absorption and 2D IR spectra that does not require knowledge of concentrations. Our method is tested on several model compounds and applied to the amide I(') band of a polypeptide in its random coil and α-helical conformation as modulated by the solution temperature. It is often difficult to confidently assign polypeptide and protein secondary structures to random coil or α-helix by linear spectroscopy alone, because they absorb in the same frequency range. We find that the transition dipole strength of the random coil state is 0.12 ± 0.013 D(2), which is similar to a single peptide unit, indicating that the vibrational mode of random coil is localized on a single peptide unit. In an α-helix, the lower bound of transition dipole strength is 0.26 ± 0.03 D(2). When taking into account the angle of the amide I(') transition dipole vector with respect to the helix axis, our measurements indicate that the amide I(') vibrational mode is delocalized across a minimum of 3.5 residues in an α-helix. Thus, one can confidently assign secondary structure based on exciton delocalization through its effect on the transition dipole strength. Our method will be especially useful for kinetically evolving systems, systems with overlapping molecular conformations, and other situations in which concentrations are difficult to determine.

Toward detecting the formation of a single helical turn by 2D IR cross peaks between the amide-I and -II modes

We have combined two-dimensional infrared (2D IR) spectroscopy and isotope substitutions to reveal the vibrational couplings between a pair of amide-I and -II modes that are several residues away but directly connected through a hydrogen bond in a helical peptide. This strategy is demonstrated on a 3(10)-helical hexapeptide, Z-Aib-L-Leu-(Aib)2-Gly-Aib-OtBu, and its 13C=18O-Leu monolabeled and 13C=18O-Leu/15N-Gly bis-labeled isotopomers in CDCl3. The isotope-dependent amide-I/II cross peaks clearly show that the second and fourth peptide linkages are vibrationally coupled as they are in proximity, forming a 3(10)-helical turn. The experimental spectra are compared to simulations based on a vibrational exciton Hamiltonian model that fully takes into account the amide-I and -II modes. The amide-II local mode frequency is evaluated by a new model based on the effects of hydrogen-bond geometry and sites. Ab initio nearest-neighbor coupling maps of the amide-I/I, -I/II, -II/I and -II/II modes are generated by isotopically isolating the local modes of N-acetyl-glycine N'-methylamide (AcGlyNHMe). Longer range couplings are modeled by transition charge interactions. The effects of the capping groups are incorporated and isotope effects are analyzed based on ab initio calculations of six model compounds. The main features of the 2D IR spectra are reproduced by this modeling. The conformational sensitivity of the isotope-dependent amide-I/II cross peaks is discussed in comparison with the calculated spectra for a semiextended structure. Our experimental and theoretical study demonstrates that the combination of 2D IR and 13C=18O/15N labeling is a useful structural method for detecting helical turn formation with residue-level specificity.

Conformation dependence of the C(alpha)D(alpha) stretch mode in peptides. II. explicitly hydrated alanine peptide structures

Our previous studies of the potential utility of the C(alpha)D(alpha) stretch frequency, nu(CD), as a tool for determining conformation in peptide systems (Mirkin and Krimm, J Phys Chem A 2004, 108, 10923-10924; 2007, 111, 5300-5303) dealt with the spectroscopic characteristics of isolated alanine peptides with alpha(R), beta, and polyproline II structures. We have now extended these ab initio calculations to include various explicit-water environments interacting with such conformers. We find that the structure-discriminating feature of this technique is in fact enhanced as a result of the conformation-specific interactions of the bonding waters, in part due to our finding (Mirkin and Krimm, J Phys Chem B 2008, 112, 15268) that C(alpha)--D(alpha)...O(water) hydrogen bonds can be present in addition to those expected between water and the CO and NH of the peptide groups. In fact, nu(CD) is hardly affected by the latter bonding but can be shifted by up to 70 cm(-1) by the former hydrogen bonds. We also discuss the factors that will have to be considered in developing the molecular dynamics (MD) treatment needed to satisfactorily take account of the influence of outer water layers on the structure of the first-layer water molecules that hydrogen bond to the peptide backbone.

Amide-I relaxation-induced hydrogen bond distortion: An intermediate in electron capture dissociation mass spectrometry of alpha-helical peptides?

Electron capture dissociation (ECD) of peptides and proteins in the gas phase is a powerful tool in tandem mass spectrometry whose current description is not sufficient to explain many experimental observations. Here, we attempt to bridge the current understanding of the vibrational dynamics in alpha-helices with the recent experimental results on ECD of alpha-helical peptides through consideration of amide-I relaxation-induced hydrogen bond distortion. Based on a single spine of H-bonded peptide units, we assume that charge neutralization upon electron capture by a charged alpha-helix excites a nearby amide-I mode, which relaxes over a few picoseconds due to Fermi resonances with intramolecular normal modes. The amide-I population plays the role of an external force, which drives the displacements of each peptide unit. It induces a large immobile contraction of the H bonds surrounding the excited site whose lifetime is about the amide-I lifetime. In addition, it creates two lattice deformations describing H bond stretchings, which propagate from the excited region toward both termini of the alpha-helix, get reflected at the termini and yield H bond contractions which move back to the excited region. Consequently, we show that H bonds experience rather large contractions whose amplitude depends on general features such as the position of the amide-I mode, the peptide length and the H bond force constants. When an H bond contraction is sufficiently large, it may promote a hydrogen atom transfer between two neighboring peptide units leading to the formation of a radical at charge site remote carbonyl carbon which is known to be a precursor to the rupture of the corresponding N[Single Bond]C(alpha) bond. The introduced here way of excitation energy generation and transfer may significantly advance ECD understanding and complement existing ECD mechanisms.

Amide-I lifetime-limited vibrational energy flow in a one-dimensional lattice of hydrogen-bonded peptide units

A time-convolutionless master equation is established for describing the amide-I vibrational energy flow in a lattice of H-bonded peptide units. The dynamics is addressed within the small polaron formalism to account for the strong coupling between the amide-I vibron and the phonons describing the H-bond vibrations. Therefore, special attention is paid to characterize the influence of the amide-I relaxation on the polaron transport properties. This relaxation is modeled by assuming that each amide-I mode interacts with a bath of intramolecular normal modes whose displacements are strongly localized on the C=O groups. It has been shown that the energy relaxation occurs over a very short time scale which prevents any significant delocalization of the polaron. At biological temperature, the polaron explores a finite region around the excited site whose size is about one or two lattice parameters. However, two regimes occur depending on whether the vibron-phonon coupling is weak or strong. For a weak coupling, the energy propagates coherently along the lattice until the polaron disappears. By contrast, for a strong coupling, a diffusive regime occurs so that the polaron explores a finite size region incoherently. In both cases, the finite polaron lifetime favors the localization of the vibron density whose amplitude decreases exponentially.

Infrared spectroscopy of proteins

This review discusses the application of infrared spectroscopy to the study of proteins. The focus is on the mid-infrared spectral region and the study of protein reactions by reaction-induced infrared difference spectroscopy.

Detecting protein-protein interactions by isotope-edited infrared spectroscopy: a numerical approach

We present a theoretical and numerical analysis of the vibrational coupling between isotope-edited amino acids in protein dimers. Depending on the presence and magnitude of coupling between 13Calpha=O peptide bond oscillators, characteristic level splittings of vibrational eigenstates are predicted. For the example of the Gramicidin A ion channel polypeptide, we observe typical IR fingerprints for the head-to-head and the antiparallel double-helical conformation of the dimer. We suggest that these findings can be used to clearly identify the structure of polypeptide aggregates using a particularly simple isotope substitution pattern.

UV−Resonance Raman Thermal Unfolding Study of Trp-Cage Shows That It Is Not a Simple Two-State Miniprotein

Trp-cage, a synthetic 20 residue polypeptide, is proposed to be an ultrafast folding synthetic miniprotein which utilizes tertiary contacts to define its native conformation. We utilized UV resonance Raman spectroscopy (UVRS) with 204 and 229 nm excitation to follow its thermal melting. Our results indicate that Trp-cage melting is complex, and it is not a simple two-state process. Using 204 nm excitation we probe the peptide secondary structure and find the Trp-cage's alpha-helix shows a broad melting curve where on average four alpha-helical amide bonds melt upon a temperature increase from 4 to 70 degrees C. Using 229 nm excitation we probe the environment of the Trp side chain and find that its immediate environment becomes more compact as the temperature is increased from 4 to 20 degrees C; however, further temperature increases lead to exposure of the Trp to water. The chi(2) angle of the Trp side chain remains invariant throughout the entire temperature range. Previous kinetic results indicated a single-exponential decay in the 4-70 degrees C temperature range, suggesting that Trp-cage behaves as a two-state folder. However, this miniprotein does not show clear two-state behavior in our steady-state studies. Rather it shows a continuous distribution of steady-state spectral parameters. Only the alpha-helix melting curve even hints of a cooperative transition. Possibly, the previous kinetic results monitor only a small region of the Trp-cage which locally appears two-state. This would then argue for spatially decoupled folding even for this small peptide.

Phonon spectra and thermodynamic properties of the infinite polyalanine alpha helix: a density-functional-theory-based harmonic vibrational analysis

We have performed a density-functional theory harmonic vibrational analysis of the infinite polyalanine alpha helix. The calculated phonon dispersion spectrum shows excellent agreement to available experimental data, except for the high frequency hydrogen stretching modes which show characteristic shifts due to anharmonic effects. A major advantage compared to previously performed empirical force field studies is that long range effects such as electrostatic interaction and polarization are intrinsically taken into account for characterizing hydrogen bond formation in the helix. Our results indicate that these effects are crucial to accurately describe the low frequency acoustical branches and lead to a significantly better agreement with experiment for the specific heat in the low temperature range.

Uncoupled Peptide Bond Vibrations in α-Helical and Polyproline II Conformations of Polyalanine Peptides

We examined the 204-nm UV resonance Raman (UVR) spectra of the polyproline II (PPII) and alpha-helical states of a 21-residue mainly alanine peptide (AP) in different H2O/D2O mixtures. Our hypothesis is that if the amide backbone vibrations are coupled, then partial deuteration of the amide N will perturb the amide frequencies and Raman cross sections since the coupling will be interrupted; the spectra of the partially deuterated derivatives will not simply be the sum of the fully protonated and deuterated peptides. We find that the UVR spectra of the AmIII and AmII' bands of both the PPII conformation and the alpha-helical conformation (and also the PPII AmI, AmI', and AmII bands) can be exactly modeled as the linear sum of the fully N-H protonated and N-D deuterated peptides. Negligible coupling occurs for these vibrations between adjacent peptide bonds. Thus, we conclude that these peptide bond Raman bands can be considered as being independently Raman scattered by the individual peptide bonds. This dramatically simplifies the use of these vibrational bands in IR and Raman studies of peptide and protein structure. In contrast, the AmI and AmI' bands of the alpha-helical conformation cannot be well modeled as a linear sum of the fully N-H protonated and N-D deuterated derivatives. These bands show evidence of coupling between adjacent peptide bond vibrations. Care must be taken in utilizing the AmI and AmI' bands for monitoring alpha-helical conformations since these bands are likely to change as the alpha-helical length changes and the backbone conformation is perturbed.

Dihedral angles of tripeptides in solution directly determined by polarized Raman and FTIR spectroscopy

The amide I mode of the peptide linkage is highly delocalized in peptides and protein segments due to through-bond and through-space vibrationally coupling between adjacent peptide groups. J. Phys. Chem. B. 104:11316-11320) used coherent femtosecond infrared (IR) spectroscopy to determine the excitonic coupling energy and the orientational angle between the transition dipole moments of the interacting amide I modes of cationic tri-alanine in D(2)O. Recently, the same parameters were determined for all protonation states of tri-alanine by analyzing the amide I bands in the respective IR and isotropic Raman spectra (. J. Am. Chem. Soc. 119:1720-1726.). In both studies, the dihedral angles phi and psi were then obtained by utilizing the orientational dependence of the coupling energy obtained from ab initio calculations on tri-glycine in vacuo (. J. Raman Spectrosc. 29:81-86) to obtain an extended 3(1) helix-like structure for the tripeptide. In the present paper, a novel algorithm for the analysis of excitonic coupling between amide I modes is presented, which is based on the approach by Schweitzer-Stenner et al. but avoids the problematic use of results from ab initio calculations. Instead, the dihedral angles are directly determined from infrared and visible polarized Raman spectra. First, the interaction energy and the corresponding degree of wave-function mixing were obtained from the amide I profile in the isotropic Raman spectrum. Second, the depolarization ratios and the amide I profiles in the anisotropic Raman and IR-absorption spectra were used to determine the orientational angle between the peptide planes and the transition dipole moments, respectively. Finally, these two geometric parameters were utilized to determine the dihedral angles phi and psi between the interacting peptide groups. Stable extended conformations with dihedral angles in the beta-sheet region were obtained for all protonation states of tri-alanine, namely phi(+) = -126 degrees, psi(+) = 178 degrees; phi(+/-) = -110 degrees, psi(+/-) = 155 degrees; and phi(-) = -127 degrees, psi(-) = 165 degrees for the cationic, zwitterionic, and anionic state, respectively. These values reflect an extended beta-helix structure. Tri-glycine was found to be much more heterogeneous in that different extended conformers coexist in the cationic and zwitterionic state, which yield a noncoincidence between isotropic and anisotropic Raman scattering. Our study introduces vibrational spectroscopy as a suitable tool for the structure analysis of peptides in solution and tripeptides as suitable model systems for investigating the role of local interactions in determining the propensity of peptide segments for distinct secondary structure motifs.

Differentiation of beta-sheet-forming structures: ab initio-based simulations of IR absorption and vibrational CD for model peptide and protein beta-sheets

Ab initio quantum mechanical computations of force fields (FF) and atomic polar and axial tensors (APT and AAT) were carried out for triamide strands Ac-A-A-NH-CH(3) clustered into single-, double-, and triple-strand beta-sheet-like conformations. Models with phi, psi, and omega angles constrained to values appropriate for planar antiparallel and parallel as well as coiled antiparallel (two-stranded) and twisted antiparallel and parallel sheets were computed. The FF, APT, and AAT values were transferred to corresponding larger oligopeptide beta-sheet structures of up to five strands of eight residues each, and their respective IR and vibrational circular dichroism (VCD) spectra were simulated. The antiparallel planar models in a multiple-stranded assembly give a unique IR amide I spectrum with a high-intensity, low-frequency component, but they have very weak negative amide I VCD, both reflecting experimental patterns seen in aggregated structures. Parallel and twisted beta-sheet structures do not develop a highly split amide I, their IR spectra all being similar. A twist in the antiparallel beta-sheet structure leads to a significant increase in VCD intensity, while the parallel structure was not as dramatically affected by the twist. The overall predicted VCD intensity is quite weak but predominantly negative (amide I) for all conformations. This intrinsically weak VCD can explain the high variation seen experimentally in beta-forming peptides and proteins. An even larger variation was predicted in the amide II VCD, which had added complications due to non-hydrogen-bonded residues on the edges of the model sheets.

Simulations of oligopeptide vibrational CD: effects of isotopic labeling

Simulated ir absorption and vibrational CD (VCD) spectra of four alanine-based octapeptides, each having its main chain constrained to a different secondary structure conformation, were analyzed and compared with experimental results for several different peptides. The octapeptide simulations were based on transfer of property tensors from a series of ab initio calculations for a short L-alanine based segment containing 3 peptide bonds with relative straight phi, psi angles fixed to those appropriate for alpha-helix, 3(10)-helix, ProII-like helix, and beta-sheet-like strand. The tripeptide force field (FF) and atomic polar tensors were obtained with density functional theory techniques at the BPW91/6-31G** level and the atomic axial tensor at the mixed BPW91/6-31G**/HF/6-31G level. Allowing for frequency correction due to the FF limitations, the octapeptide results obtained are qualitatively consistent with experimental observations for ir and VCD spectra of polypeptides and oligopeptides in established conformations. In all cases, the correct VCD sign patterns for the amide I and II bands were predicted, but the intensities did have some variation from the experimental patterns. Predicted VCD changes upon deuteration of either the peptide or side-chains as well as for (13)C isotopic labeling of the amide C=O at specific sites in the peptide chain were computed for analysis of experimental observations. A combination of theoretical modeling with experimental data for labeled compounds leads both to enhanced resolution of component transitions and added conformational applicability of the VCD spectra.