Theoretical study of internal field effects on peptide amide I modes

Charged terminal groups or polar side chains of amino acids create spatially nonuniform electrostatic potential around intramolecular peptide bonds and induce amide I mode frequency shifts in polypeptides. By carrying out a series of quantum chemistry calculation studies of various ionic di- and tripeptides as well as dipeptides of 20 different amino acids, these internal field effects on vibrational properties are theoretically investigated. The amide I local and normal mode frequencies and dipole and rotational strengths determining IR and vibrational circular dichroism intensities, respectively, are found to depend on the polar nature of side chains, whereas the vibrational coupling strength weakly does so. The empirical correction and fragment analysis methods were used to theoretically calculate the amide I local mode frequencies and dipole and rotational strengths. These values were directly compared with ab initio and density functional theory calculation results, and the agreements were found to be quantitative.     

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     

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.     

Amide I vibrational circular dichroism of dipeptide: Conformation dependence and fragment analysis

The amide I vibrational circular dichroic response of alanine dipeptide analog (ADA) was theoretically investigated and the density functional theory calculation and fragment analysis results are presented. A variety of vibrational spectroscopic properties, local and normal mode frequencies, coupling constant, dipole, and rotational strengths, are calculated by varying two dihedral angles determining the three-dimensional ADA conformation. Considering two monopeptide fragments separately, we show that the amide I vibrational circular dichroism of the ADA can be quantitatively predicted. For several representative conformations of the model ADA, vibrational circular dichroism spectra are calculated by using both the density functional theory calculation and fragment analysis methods.     

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.     

振动圆二色谱: 一种确定手性分子绝对构型的新方法

手性分子绝对构型的确定是一个极其重要且长期存在的问题. 振动圆二色谱是在红外波长区域测定分子圆二色性的一种新方法, 极大地扩展了圆二色谱的应用范围. 振动圆二色谱法通过构象搜索、量子化学计算等手段准确预测手性分子的振动圆二色谱图, 进而与实测谱图进行比较确定其绝对构型. 该方法已经得到了越来越广泛的应用, 必将成为一种有效测定手性分子绝对构型的常规方法.     

Discrimination between peptide 3(10)- and alpha-helices. Theoretical analysis of the impact of alpha-methyl substitution on experimental spectra

Detailed spectral simulations based on ab initio density functional theory computations of the amide I and II infrared (IR) and vibrational circular dichroism (VCD) spectra for Ac-(Ala)(4)-NH(2), Ac-(Aib-Ala)(2)-NH(2), and Ac-(Aib)(4)-NH(2) constrained to 3(10)- and alpha-helical conformations are presented. Parameters from these ab initio calculations are transferred onto corresponding larger oligopeptides to simulate the spectra for dodecamers. The differences between conformations and for different Aib substitution patterns within a conformation are reflected in observable spectral patterns where data are available. Simulated IR spectra show small frequency shifts in the amide I maxima between 3(10)- and alpha-helices, but the same magnitude shifts occur within one conformation upon Aib substitution. Thus, from a computational basis, the frequency of the amide I maximum does not discriminate between the 3(10)- or alpha-helical conformations. Calculated VCD band shapes for 3(10)-helices showed more significant changes in amplitude, with change in the fraction of Aib, than those for alpha-helices. Generally, with increasing Aib content, the overall amide I VCD intensity becomes weaker and the amide I couplet becomes more conservative, while the amide II VCD is less affected. Although the detailed band shape is shown to be sensitive to alpha-Me substitution, the basic pattern of amide I and II relative VCD intensities still differs between alpha- and 3(10)-helices and, as a consequence, successfully discriminates between them. These predictions are all borne out in experimental spectra of Aib, mixed Aib-Ala, and Ala-based helical peptides, where available.     

Optical Spectra of Unfolded Proteins: A Partially Ordered Polymer Problem

Proteins are hetero-sequence polypeptides with complex folded structures whose topology and structural details are vital to their biological function. In this paper, uses of electronic circular dichroism (ECD) in the uv, infrared (IR) absorption and its chiroptical variant, vibrational circular dichroism (VCD) for study of residual structure in peptide models and unfolded proteins under various conditions are addressed. The complementary information gained from analysis of the short range vibrational coupling (with IR and VCD) vs. long range dipole coupling (from ECD) allow detection of partial ordering. The vibrational techniques show an ability to detect local order often missed by ECD (or fluorescence). Perturbation with thermal and solvent variation can then probe stabilities. Furthermore the rapid timescales allow such spectra to follow dynamic changes. Examples from protein folding of cytochrome c and various beta containing proteins are given.     

Amide I vibrational circular dichroism of polypeptides: generalized fragmentation approximation method

Fragment analyses of vibrational circular dichroic response of dipeptides were carried out recently [Choi and Cho, J. Chem. Phys. 120, 4383 (2004)]. In the present paper, by using a minimal size unit peptide containing two chiral carbons covalently bonded to the peptide group, a generalized fragmentation approximation method is discussed and applied to the calculations of infrared-absorption and vibrational circular dichroism (VCD) intensities of amide I vibrations in various secondary structure polypeptides. Unlike the dipole strength determining IR-absorption intensity, the rotational strength is largely determined by the cross terms that are given by the inner product between the transition electric dipole and the transition magnetic dipole of two different peptides. This explains why the signs and magnitudes of VCD peaks are far more sensitive to the relative orientation and distance between different peptide bonds in a given protein. In order to test the validity of fragmentation approximation, three different segments in a globular protein ubiquitin, i.e., right-handed alpha-helix, beta-sheet, and beta-turn regions, were chosen for density-functional theory (DFT) calculations of amide I vibrational properties and the numerically simulated IR-absorption and VCD spectra by using the fragmentation method are directly compared with DFT results. It is believed that the fragmentation approximation method will be of use in numerically simulating vibrational spectra of proteins in solutions.     

Phosphorylation effect on the GSSS peptide conformation in water: infrared, vibrational circular dichroism, and circular dichroism experiments and comparisons with molecular dynamics simulations

The phosphorylation effect on the small peptide conformation in water has not been clearly understood yet, despite the widely acknowledged notion that control of protein activity by phosphorylation works mainly by inducing conformational change. To elucidate the detailed mechanism, we performed infrared (IR) absorption and vibrational and electronic circular dichroism studies of both unphosphorylated and phosphorylated tetrapeptides, GSSS 1 and GSSpS 2. The solution structure of the tetrapeptide is found to be little dependent on the presence of the neutral or negatively charged phosphoryl group, and to be a mixture of extended structures including polyproline II (PII) and beta-sheet conformations. The additional band at 1598 cm(-1) in the amide I IR spectrum of the phosphorylated peptide GSSpS at neutral pD appears to be clear spectroscopic evidence for direct intramolecular hydrogen-bonding interaction between the side chain dianionic phosphoryl group and the backbone amide proton. On the basis of amide I IR band analyses, the authors found that the probability of finding the phosphoryl group strongly H bonded to the backbone proton in GSSpS is about 43% at pD 7.0 and 37 degrees C. Such a H-bonding interaction in GSSpS has the biological standard enthalpy and entropy of -15.1 kJ/mol and -51.2 J/K mol, respectively. Comparisons between the experimentally measured IR and VCD spectra and the numerically simulated ones suggested that the currently available force field parameters need to be properly modified. The results in this paper may shed light on an unknown mechanism of controlling the peptide conformation by phosphorylation.     

Analysis of vibrational circular dichroism spectra of (s)-(+)-2-butanol by rotational strengths expressed in local symmetry coordinates

The vibrational circular dichroism (VCD) spectra of (S)-(+)-2-butanol have been observed in dilute CS(2) solutions. Two strong VCD bands are assigned mainly to the OH bending modes with the aid of quantum chemical calculations. The calculated VCD spectra corresponding to these bands are shown to depend on the conformation of the OH group. To understand this feature, we have calculated the contribution of each local vibrational mode to the rotational strengths and concluded that the coupling of the group vibrations between the in-plane and out-of-plane modes about the locally assumed symmetry planes play a significant role in VCD. This finding has provided a new scope of VCD in relation to molecular vibrations.     

Infrared and vibrational CD spectra of partially solvated alpha-helices: DFT-based simulations with explicit solvent

Theoretical simulations are used to investigate the effects of aqueous solvent on the vibrational spectra of model alpha-helices, which are only partly exposed to solvent to mimic alpha-helices in proteins. Infrared absorption (IR) and vibrational circular dichroism (VCD) amide I' spectra for 15-amide alanine alpha-helices are simulated using density functional theory (DFT) calculations combined with the property transfer method. The solvent is modeled by explicit water molecules hydrogen bonded to the solvated amide groups. Simulated spectra for two partially solvated model alpha-helices, one corresponding to a more exposed and the other to a more buried structure, are compared to the fully solvated and unsolvated (gas phase) simulations. The dependence of the amide I spectra on the orientation of the partially solvated helix with respect to the solvent and effects of solvation on the amide I' of 13C isotopically substituted alpha-helices are also investigated. The partial exposure to solvent causes significant broadening of the amide I' bands due to differences in the vibrational frequencies of the explicitly solvated and unsolvated amide groups. The different degree of partial solvation is reflected primarily in the frequency shifts of the unsolvated (buried) amide group vibrations. Depending on which side of the alpha-helix is exposed to solvent, the simulated IR band-shapes exhibit significant changes, from broad and relatively featureless to distinctly split into two maxima. The simulated amide I' VCD band-shapes for the partially solvated alpha-helices parallel the broadening of the IR and exhibit more sign variation, but generally preserve the sign pattern characteristic of the alpha-helical structures and are much less dependent on the alpha-helix orientation with respect to the solvent. The simulated amide I' IR spectra for the model peptides with explicitly hydrogen-bonded water are consistent with the experimental data for small alpha-helical proteins at very low temperatures, but overestimate the effects of solvent on the protein spectra at ambient temperatures, where the peptide-water hydrogen bonds are weakened by thermal motion.     

Dipeptide structure determination by vibrational circular dichroism combined with quantum chemistry calculations.

The solution structure and the local solvation environments of alanine dipeptide (AD, 1 a) and its isotopomer (AD*, 1 b, 13C on the acetyl end C==O) are studied by using infrared (IR) spectroscopy and vibrational circular dichroism (VCD). From the amide I IR spectra of AD* in various protic solvents, it is found that each of the two carbonyl groups is fully H-bonded to two water molecules. However, the number of alcohol molecules H-bonded to each C==O varies from one to two, and the local solvation environments are asymmetric around the two peptides of AD* in alcohol solutions. The amide I VCD spectra of AD and AD* in D2O are also measured, and a series of density functional theory (DFT, B3LYP/6-311++G**) calculations are performed to obtain the amide I normal-mode rotational strengths of AD and the intrinsic rotational strengths of its two peptide fragments. By combining the VCD-measurement and DFT-calculation results and employing a coupled oscillator theory, we show that the aqueous-solution structure of the dipeptide can be determined. We believe that the present method will be of use in building up a library of dipeptide solution structures in water.     

Beta Sheets with a Twist: The Conformation of Helical Polyisocyanopeptides Determined by Using Vibrational Circular Dichroism

Detailed information on the architecture of polyisocyanopeptides based on vibrational circular dichroism (VCD) spectroscopy in combination with DFT calculations is presented. It is demonstrated that the screw sense of the helical polyisocyanides can be determined directly from the C?N‐stretch vibrational region of the VCD spectrum. Analysis of the VCD signals associated with the amide I and amide II modes provides detailed information on the peptide side‐chain arrangement in the polymer and indicates the presence of a helical β‐sheet architecture, in which the dihedral angles are slightly different to those of natural β‐sheet helices.     

Conformational manifold of alpha-aminoisobutyric acid (Aib) containing alanine-based tripeptides in aqueous solution explored by vibrational spectroscopy, electronic circular dichroism spectroscopy, and molecular dynamics simulations

Replacement of the alpha-proton of an alanine residue to generate alpha-aminoisobutyric acid (Aib) in alanine-based oligopeptides favors the formation of a 3(10) helix when the length of the oligopeptide is about four to six residues. This research was aimed at experimentally identifying the structural impact of an individual Aib residue in an alanine context of short peptides in water and Aib's influence on the conformation of nearest-neighbor residues. The amide I band profile of the IR, isotropic and anisotropic Raman, and vibrational circular dichroism (VCD) spectra of Ac-Ala-Ala-Aib-OMe, Ac-Ala-Aib-Ala-OMe, and Ac-Aib-Ala-Ala-OMe were measured and analyzed in terms of different structural models by utilizing an algorithm that exploits the excitonic coupling between amide I' modes. The conformational search was guided by the respective 1H NMR and electronic circular dichroism spectra of the respective peptides, which were also recorded. From these analyses, all peptides adopted multiple conformations. Aib predominantly sampled the right-handed and left-handed 3(10)-helix region and to a minor extent the bridge region between the polyproline (PPII) and the helical regions of the Ramachandran plot. Generally, alanine showed the anticipated PPII propensity, but its conformational equilibrium was shifted towards helical conformations in Ac-Aib-Ala-Ala-OMe, indicating that Aib can induce helical conformations of neighboring residues positioned towards the C-terminal direction of the peptide. An energy landscape exploration by molecular dynamics simulations corroborated the results of the spectroscopic studies. They also revealed the dynamics and pathways of potential conformational transitions of the corresponding Aib residues.     

A matrix isolation study on Ac-Gly-NHMe and Ac-l-Ala-NHMe, the simplest chiral and achiral building blocks of peptides and proteins

The infrared absorption (IR) spectra of acetyl-N-methyl-glycine and acetyl-N-methyl-alanine have been recorded in dichloromethane and dimethyl sulfoxide-d(6) solution, as well as in Ar and Kr matrices. The spectra were assigned with the help of quantum chemical calculations. Based on the assignments of the matrix-isolation IR spectra, in line with theoretical predictions, two different hydrogen bonded conformers were identified, furthermore a third conformer is likely to be present, which cannot be unambiguously identified. In dichloromethane two conformers could be observed, while in dimethyl sulfoxide a single conformer could be identified. Vibrational circular dichroism (VCD) spectra of acetyl-N-methyl-l-alanine have also been recorded in solutions and matrices. These matrix-isolation VCD spectra not only support the assignments of the matrix-isolation IR spectra, but also demonstrate that these spectra can be interpreted much easier with the help of quantum chemical calculations than the VCD spectra recorded in solutions. It is also shown that the rotatory strength of some vibrational transitions changes rapidly as a function of the backbone torsional coordinates; hence the appearance of some regions in the VCD spectra is extremely dependent on any perturbations, e.g. weak intermolecular interactions.     

Nature of vibrational coupling in helical peptides: an isotopic labeling study

Infrared (IR) and vibrational circular dichroism (VCD) spectra were measured for a series of isotopically ((13)C on two or more amide Cdouble bond]O) labeled, 25 residue, alpha-helical peptides of the sequence Ac-(AAAAK)(4)AAAAY-NH(2) that were also studied in the previous paper. Theoretical IR and VCD simulations were performed for correspondingly isotopically labeled Ac-A(24)-NHCH(3) constrained to an alpha-helical conformation by use of property tensor transfer from density functional theory (DFT) calculations on Ac-A(10)-NHCH(3). The simulations predicted and experiments confirmed that the vibrational coupling constants between i, i + 1 and i, i + 2 residues differ in sign, thus leading to a reversal of the (13)C VCD pattern and explaining the large shift in the (13)C amide I frequency as reported in the previous paper. The sign of the coupling constant remained consistent for larger label separation (with the exception of i, i + 4) and for more labels with uniform separation. Such effects confirm that the isotopically labeled group vibrations are essentially only coupled to each other and are effectively uncoupled from those of the unlabeled groups. This development confirms the utility of isotopic labels for site-specific structural studies with vibrational spectra. Observed spectral effects cannot be explained by considering only transition dipole coupling (TDC) between amide oscillators, particularly for smaller label separations, but the TDC and ab initio predicted couplings roughly converge at large separation.     

A single chiroptical spectroscopic method may not be able to establish the absolute configurations of diastereomers: dimethylesters of hibiscus and garcinia acids

Electronic circular dichroism (ECD), optical rotatory dispersion (ORD), and vibrational circular dichroism (VCD) spectra of hibiscus acid dimethyl ester have been measured and analyzed in combination with quantum chemical calculations of corresponding spectra. These results, along with those reported previously for garcinia acid dimethyl ester, reveal that none of these three (ECD, ORD, or VCD) spectroscopic methods, in isolation, can unequivocally establish the absolute configurations of diastereomers. This deficiency is eliminated when a combined spectral analysis of either ECD and VCD or ORD and VCD methods is used. It is also found that the ambiguities in the assignment of absolute configurations of diastereomers may also be overcome when unpolarized vibrational absorption is included in the spectral analysis.     

Fenchone, camphor, 2-methylenefenchone and 2-methylenecamphor: a vibrational circular dichroism study

We report and discuss the infrared (IR) vibrational circular dichroism (VCD) spectra of the enantiomeric pairs of the olefin derivatives of fenchone (1,3,3-trimethyl-2-methylenebicyclo[2.2.1]heptane) and camphor (1,7,7-trimethyl-2-methylenebicyclo[2.2.1]heptane), respectively, together with those of the parent molecules. The VCD spectra were taken in three spectral regions: the mid-IR region, encompassing the fundamental deformation modes, the region of CH-stretching fundamental modes and the NIR-region between 1100 and 1300 nm, which corresponds to the second CH-stretching overtone. The VCD and absorption spectra in the first two regions are analyzed by use of current density functional theory (DFT) calculations. The NIR region is analyzed by a protocol that consists of the use of DFT-based calculations and in assuming local mode behavior: the local mode approach is found appropriate for interpreting the absorption spectra and, for the moment, acceptable for calculating NIR-VCD spectra. The analysis of the first region allows us to track the contribution of the C=O group in the vibrational optical activity of C-C stretching modes; notable differences are indeed found in olefins and ketones. On the contrary, in the other two regions the VCD spectra of olefins and ketones are more similar: in the normal mode region of CH stretching fundamentals the spectra are determined by the mutual orientation of the CH bonds; in the second overtone local mode region olefins and ketones signals show some differences.