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     

Ab initio modeling of amide I coupling in antiparallel beta-sheets and the effect of 13C isotopic labeling on infrared spectra

Isotopic substitution with 13C on the amide C=O has become an important means of determining localized structural information about peptide conformations with vibrational spectroscopy. Various approaches to the modeling of the interactions between labeled amide sites, specifically for antiparallel two-stranded, beta-forming peptides, were investigated, including different force fields [dipole-dipole interaction vs density functional theory (DFT) treatments], basis sets, and sizes of model peptides used for ab initio calculations, as well as employing models of solvation. For these beta-sheet systems the effect of the relative positions of the 13C isotopic labels in each strand on their infrared spectra was investigated. The results suggest that the interaction between labeled amide groups in different strands can be used as an indicator of local beta-structure formation, because coupling between close-lying C=O groups on opposing chains leads to the largest frequency shifts, yet some alternate placements can lead to intensity enhancements. The basic character of the coupling interaction between labeled modes on opposing strands is independent of changes in peptide length, water solvent environment, twisting of the sheet structure, and basis set used in the calculations, although the absolute frequencies and detailed coupling magnitudes change under each of these perturbations. In particular, two strands of three amides each contain the basic interactions needed to simulate larger sheets, with the only exception that the C=O groups forming H-bonded rings at the termini can yield different coupling values than central ones of the same structure. Spectral frequencies and intensities were modeled ab initio by DFT primarily at the BPW91/6-31G** level for pairs of three, four, and six amide strands. Comparison to predictions of a classical coupled oscillator model show qualitative but not quantitative agreement with these DFT results.     

The anharmonic vibrational potential and relaxation pathways of the amide I and II modes of N-methylacetamide

We investigate the influence of isotopic substitution and solvation of N-methylacetamide (NMA) on anharmonic vibrational coupling and vibrational relaxation of the amide I and amide II modes. Differences in the anharmonic potential of isotopic derivatives of NMA in D2O and DMSO-d6 are quantified by extraction of the anharmonic parameters and the transition dipole moment angles from cross-peaks in the two-dimensional infrared (2D-IR) spectra. To interpret the effects of isotopic substitution and solvent interaction on the anharmonic potential, density functional theory and potential energy distribution calculations are performed. It is shown that the origin of anharmonic variation arises from differing local mode contributions to the normal modes of the NMA isotopologues, particularly in amide II. The time domain manifestation of the coupling is the coherent exchange of excitation between amide modes seen as the quantum beats in femtosecond pump-probes. The biphasic behavior of population relaxation of the pump-probe and 2D-IR experiments can be understood by the rapid exchange of strongly coupled modes within the peptide backbone, followed by picosecond dissipation into weakly coupled modes of the bath.     

A density functional theory study of vibrational coupling in the amide I band of beta-sheet models

We report the first molecular orbital/density functional theory (DFT) calculations on the vibrational frequencies involved in the amide I band of completely geometrically optimized models for beta-sheet peptides based upon (up to 16) glycine residues. These calculations use the B3LYP/D95** level of DFT. The primary means of vibrational coupling occurs through H bond, rather than through space, interactions, which is consistent with a previous report on alpha-helical polyalanines and H-bonding chains of both formamides and 4-pyridones. We decoupled the C=O stretching vibrations using selected 14C substitutions to probe the coupling mechanism and to determine "natural" frequencies for individual 14C=Os. The intermolecular H-bonding interactions affect the geometries of the amide groups. Those near the center of H-bonding chains have long C=O bonds. The C=O bond lengths correlate with these "natural" frequencies, The frequencies obtained from the DFT calculations are generally more coupled, and the most intense are more red shifted than those calculated by transition dipole coupling (TDC). TDC inverts the order of the shifted frequencies compared to DFT in several cases.     

Simulation of infrared spectra for beta-hairpin peptides stabilized by an Aib-Gly turn sequence: correlation between conformational fluctuation and vibrational coupling

Vibrational spectra of a 12-residue beta-hairpin peptide, RYVEVBGKKILQ (HBG), stabilized by an Aib-Gly turn sequence (B = Aib) were investigated theoretically using a combination of molecular dynamics (MD) and density functional theory (DFT) calculations. Selected conformations of HBG were extracted from a classical MD trajectory and used for spectral simulations. DFT calculations, based on the Cartesian coordinate spectral property transfer protocol, were carried out for peptide structures in which all residues are replaced with Ala, except for the Aib and Gly residues, but the backbone (phi, psi, omega) structure of the original configuration is retained. The simulations provide a basis for interpretation of the HBG amide I infrared spectra in terms of structural variables such as detailed secondary structure and thermal conformational fluctuation as well as vibrational coupling as indicated by spectra of 13C isotope-labeled variants. The characteristic amide I band shape of such small beta-hairpin peptides appears to arise from the structure of the short antiparallel beta-sheet strands. The role of structural parameter fluctuation in vibrational coupling is evaluated by comparison of DFT-derived amide coupling constants for selected configurations and from transition dipole coupling calculations of coupling parameters between (13)C isotopically labeled residues for a MD-derived ensemble of configurations. Calculated results were compared with the experimentally obtained spectra for several (13)C isotope-labeled peptides of this sequence.     

Competition between amide stacking and intramolecular H bonds in γ-peptide derivatives: controlling nearest-neighbor preferences

Resonant two-photon ionization (R2PI), IR-UV holeburning (IR-UV), and resonant ion-dip infrared spectroscopy (RIDIRS) have been used to record mass-selected, single-conformation ultraviolet and infrared spectra of three simple diamide derivatives of γ-amino acids as isolated molecules cooled in a supersonic expansion. This work builds on an earlier study of Ac-γ(2)-hPhe-NHMe (James, W. H., III, et al. J. Am. Chem. Soc. 2009, 131, 14243), which showed that this methyl-capped γ-peptide forms amide-stacked conformations that are similar in stability to H-bonded conformations containing a C9 ring and more stable than C7 H-bonded ring structures. Among the γ-peptides discussed here, Ac-γ(2)-hPhe-N(Me)(2) contains an additional methyl group relative to the previously studied Ac-γ(2)-hPhe-NHMe and therefore lacks the amide NH group responsible for C9 ring formation. Three conformations of Ac-γ(2)-hPhe-N(Me)(2) are observed, all of which are amide-stacked structures. In a second new molecule, Ac-γ(2)-hPhe-NH(iPr), the C-terminal NHMe group of Ac-γ(2)-hPhe-NHMe is replaced with an NH(iPr) group. Three conformations of Ac-γ(2)-hPhe-NH(iPr) are observed, all of which are C9 H-bonded structures. The dramatic difference between C-terminal NHMe and NH(iPr) reveals the delicate balance of noncovalent forces within these γ-peptides. The third molecule we examined is a gabapentin-derived diamide (designated 1), which contains a phenylacyl group at the N-terminus and an N(Me)(2) group at the C-terminus; the latter precludes C9 H bonding. Comparison of 1 with Ac-γ(2)-hPhe-N(Me)(2) allows us to examine the impact of the backbone substitution pattern (monosubstitution at carbon-2 vs disubstitution at carbon-3) on the competition between the C7 H-bonded and the amide-stacked conformation. In this case, only C7 rings are observed. The different gas-phase behaviors observed among the molecules analyzed here offer insight on the intrinsic conformational propensities of the γ-peptide backbone, information that provides a foundation for future foldamer design efforts.     

Modeling the amide I bands of small peptides

In this paper different floating oscillator models for describing the amide I band of peptides and proteins are compared with density functional theory (DFT) calculations. Models for the variation of the frequency shifts of the oscillators and the nearest-neighbor coupling between them with respect to conformation are constructed from DFT normal mode calculations on N-acetyl-glycine-N(')-methylamide. The calculated frequencies are compared with those obtained from existing electrostatic models. Furthermore, a new transition charge coupling model is presented. We suggest a model which combines the nearest-neighbor maps with long-range interactions accounted for using the new transition charge model and an existing electrostatic map for long-range interaction frequency shifts. This model and others, which account for the frequency shifts by electrostatic maps exclusively, are tested by comparing the predicted IR spectra with those from DFT calculations on the pentapeptide [Leu]-enkephalin. The new model described above gives the best agreement and, after a systematic blueshift is accounted for, reproduces the DFT frequencies to within 3.5 cm(-1). The correlation of the intensities for this model with intensities from DFT calculations is 0.94.     

Computational spectroscopy of ubiquitin: comparison between theory and experiments

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

Peptide bond vibrational coupling

Neutral trialanine (Ala3), which is geometrically constrained to have its peptide bond at Phi and Psi angles of alpha-helix and PPII-like conformers, are studied at the B3LYP/6-31+G(d,p) level of theory to examine vibrational interactions between adjacent peptide units. Delocalization of the amide I, amide II, and amide III3 vibrations are analyzed by calculating their potential energy distributions (PED). The vibrational coupling strengths are estimated from the frequency shifts between the amide vibrations of Ala3 and the local amide bond vibrations of isotopically substituted Ala3 derivatives. Our calculations show the absence of vibrational coupling of the amide I and amide II bands in the PPII conformations. In contrast, the alpha-helical conformation shows strong coupling between the amide I vibrations due to the favorable orientation of the C=O bonds and the strong transitional dipole coupling. The amide III3 vibration shows weak coupling in both the alpha-helix and PPII conformations; this band can be treated as a local independent vibration. Our calculated results in general agree with our previous experimental UV Raman studies of a 21-residue mainly alanine-based peptide (AP).     

Electrostatic DFT map for the complete vibrational amide band of NMA

An anharmonic vibrational Hamiltonian for the amide I, II, III, and A modes of N-methyl acetamide (NMA), recast in terms of the 19 components of an external electric field and its first and second derivative tensors (electrostatic DFT map), is calculated at the DFT(BPW91/6-31G(d,p)) level. Strong correlations are found between NMA geometry and the amide frequency fluctuations calculated using this Hamiltonian together with the fluctuating solvent electric field obtained from the MD simulations in TIP3 water. The amide I and A frequencies are strongly positively correlated with the C=O and N-H bond lengths. The C=O and C-N amide bond lengths are negatively correlated, suggesting the solvent-induced fluctuations of the contribution of zwitterionic resonance form. Sampling the global electric field in the entire region of the transition charge densities (TCDs) is required for accurate infrared line shape simulations. Collective electrostatic solvent coordinates which represent the fluctuations of the 10 lowest amide fundamental and overtone states are reported. Normal-mode analysis of an NMA-3H(2)O cluster shows that the 660 cm(-1) to 1100 cm(-1) oscillation found in the frequency autocorrelation functions of the amide modes may be ascribed to the two bending vibrations of intermolecular hydrogen bonds with the amide oxygen of NMA.     

Vibrational spectral signatures of peptide secondary structures: N-methylation and side chain hydrogen bond in cyclosporin A

Molecular dynamics simulations are performed to explore important conformations of cyclosporin A, an immunosuppressive cyclic undecapeptide drug, in different media including gas-phase, chloroform, and acetonitrile. Density functional theory calculations are used to refine the low-lying conformers and to predict their infrared and vibrational circular dichroism spectra. Vibrational spectral signatures in the important amide II, I, and A regions are identified for typical peptide secondary structures including β-turn (type II' or I), antiparallel β-sheet (flat or twisted), inverse γ-turn, N-methylated peptide bond, and side chain H-bond. New insights into the spectral signatures of secondary structures especially with N-methylation and side chain hydrogen bond are provided, which can be very useful for further peptide conformation analysis in general.     

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.     

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.     

Computation of the amide I band of polypeptides and proteins using a partial Hessian approach

A partial Hessian approximation for the computation of the amide I band of polypeptides and proteins is introduced. This approximation exploits the nature of the amide I band, which is largely localized on the carbonyl groups of the backbone amide residues. For a set of model peptides, harmonic frequencies computed from the Hessian comprising only derivatives of the energy with respect to the displacement of the carbon, oxygen, and nitrogen atoms of the backbone amide groups introduce mean absolute errors of 15 and 10 cm(-1) from the full Hessian values at the Hartree-Fock/STO-3G and density functional theory EDF16-31G(*) levels of theory, respectively. Limiting the partial Hessian to include only derivatives with respect to the displacement of the backbone carbon and oxygen atoms yields corresponding errors of 24 and 22 cm(-1). Both approximations reproduce the full Hessian band profiles well with only a small shift to lower wave number. Computationally, the partial Hessian approximation is used in the solution of the coupled perturbed Hartree-Fock/Kohn-Sham equations and the evaluation of the second derivatives of the electron repulsion integrals. The resulting computational savings are substantial and grow with the size of the polypeptide. At the HF/STO-3G level, the partial Hessian calculation for a polypeptide comprising five tryptophan residues takes approximately 10%-15% of the time for the full Hessian calculation. Using the partial Hessian method, the amide I bands of the constituent secondary structure elements of the protein agitoxin 2 (PDB code 1AGT) are calculated, and the amide I band of the full protein estimated.     

Determination of conformational preferences of dipeptides using vibrational spectroscopy

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

Single-conformation ultraviolet and infrared spectroscopy of model synthetic foldamers: beta-peptides Ac-beta3-hPhe-NHMe and Ac-beta3-hTyr-NHMe

The conformational preferences and infrared and ultraviolet spectral signatures of two model beta-peptides, Ac-beta3-hPhe-NHMe (1) and Ac-beta3-hTyr-NHMe (2), have been explored under jet-cooled, isolated molecule conditions. The mass-resolved, resonant two-photon ionization spectra of the two molecules were recorded in the region of the S0-S1 origin of the phenyl or phenol ring substituents, respectively. UV-UV hole-burning spectroscopy was used to determine that two conformations of 1 are present, with the transitions due to conformer A, with S0-S1 origin at 34431 cm(-1), being almost 20 times larger than those due to conformer B, with S0-S1 origin at 34404 cm(-1). Only one conformation of 2 was observed. Resonant ion-dip infrared spectroscopy provided single-conformation infrared spectra in the 3300-3700 cm(-1) region. The spectra of conformer A of both molecules have H-bonded and free amide NH stretch infrared transitions at 3400 and 3488 cm(-1), respectively, while conformer B of 1 possesses bands at 3417 and 3454 cm(-1). For comparison with experiment, full optimizations of all low-lying minima of 1 were carried out at the DFT B3LYP/6-31+G* and RIMP2/aug-cc-pVDZ levels of theory, and single point MP2/6-31+G* calculations at the DFT geometries. On the basis of the comparison with previous studies in solution and the calculated results, conformer A of 1 and 2 were assigned to a C6 conformer, while conformer B of 1 was assigned to a unique C8 structure with a weak intramolecular H-bond. The reasons for the preference for C6 over C8 structures and the presence of only two conformations in the jet-cooled spectrum are discussed in light of the predictions from calculations.     

Entropy-driven population distributions in a prototypical molecule with two flexible side chains: O-(2-acetamidoethyl)-N-acetyltyramine

Resonant two-photon ionization (R2PI), resonant ion-dip infrared (RIDIR), and UV-UV hole-burning spectroscopies have been employed to obtain conformation-specific infrared and ultraviolet spectra under supersonic expansion conditions for O-(2-acetamidoethyl)-N-acetyltyramine (OANAT), a doubly substituted aromatic in which amide-containing alkyl and alkoxy side chains are located in para positions on a phenyl ring. For comparison, three single-chain analogs were also studied: (i) N-phenethyl-acetamide (NPEA), (ii) N-(p-methoxyphenethyl-acetamide) (NMPEA), and (iii) N-(2-phenoxyethyl)-acetamide (NPOEA). Six conformations of OANAT have been resolved, with S(0)-S(1) origins ranging from 34,536 to 35,711 cm(-1), denoted A-F, respectively. RIDIR spectra show that conformers A-C each possess an intense, broadened amide NH stretch fundamental shifted below 3400 cm(-1), indicative of the presence of an interchain H bond, while conformers D-F have both amide NH stretch fundamentals in the 3480-3495 cm(-1) region, consistent with independent-chain structures with two free NH groups. NPEA has a single conformer with S(0)-S(1) origin at 37,618 cm(-1). NMPEA has three conformers, two that dominate the R2P1 spectrum, with origin transitions between 35,580 and 35,632 cm(-1). Four conformations, one dominate and three minor, of NPOEA have been resolved with origins between 35,654 and 36,423 cm(-1). To aid the making of conformational assignments, the geometries of low-lying structures of all four molecules have been optimized and the associated harmonic vibrational frequencies calculated using density functional theory (DFT) and RIMP2 methods. The S(0)-S(1) adiabatic excitation energies have been calculated using the RICC2 method and vertical excitation energies using single-point time-dependent DFT. The sensitivity of the S(0)-S(1) energy separation in OANAT and NPOEA primarily arises from different orientations of the chain attached to the phenoxy group. Using the results of the single-chain analogs, tentative assignments have been made for the observed conformers of OANAT. The RIMP2 calculations predict that interchain H-bonded conformers of OANAT are 25-30 kJ/mol more stable than the extended-chain structures. However, the free energies of the interchain H-bonded and extended structures calculated at the preexpansion temperature (450 K) differ by less than 10 kJ/mol, and the number of extended structures far outweighs the number of H-bonded conformers. This entropy-driven effect explains the presence of the independent-chain conformers in the expansion, and cautions future studies that rely solely on relative energies of conformers in considering possible assignments.     

A polarizable force field for computing the infrared spectra of the polypeptide backbone

The shapes of the amide bands in the infrared (IR) spectra of proteins and peptides are caused by electrostatically coupled vibrations within the polypeptide backbone and code the structures of these biopolymers. A structural decoding of the amide bands has to resort to simplified models because the huge size of these macromolecules prevents the application of accurate quantum mechanical methods such as density functional theory (DFT). Previous models employed transition-dipole coupling methods that are of limited accuracy. Here we propose a concept for the computation of protein IR spectra, which describes the molecular mechanics (MM) of polypeptide backbones by a polarizable force field of "type II". By extending the concepts of conventional polarizable MM force fields, such a PMM/II approach employs field-dependent parameters not only for the electrostatic signatures of the molecular components but also for the local potentials modeling the stiffness of chemical bonds with respect to elongations, angle deformations, and torsions. Using a PMM/II force field, the IR spectra of the polypeptide backbone can be efficiently calculated from the time dependence of the backbone's dipole moment during a short (e.g., 100 ps) MD simulation by Fourier transformation. PMM/II parameters are derived for harmonic bonding potentials of amide groups in polypeptides from a series of DFT calculations on the model molecule N-methylacetamide (NMA) exposed to homogeneous external electric fields. The amide force constants are shown to vary by as much as 20% for relevant field strengths. As a proof of principle, it is shown that the large solvatochromic effects observed in the IR spectra of NMA upon transfer from the gas phase into aqueous solution are not only excellently reproduced by DFT/MM simulations but are also nicely modeled by the PMM/II approach. The tasks remaining for a proof of practice are specified.     

Single-conformation ultraviolet and infrared spectroscopy of model synthetic foldamers: beta-peptides Ac-beta3-hPhe-beta3-hAla-NHMe and Ac-beta3-hAla-beta3-hPhe-NHMe

The conformational preferences and infrared and ultraviolet spectral signatures of two model beta-peptides, Ac-beta3-hPhe-beta3-hAla-NHMe (1) and Ac-beta3-hAla-beta3-hPhe-NHMe (2), have been explored under jet-cooled, isolated-molecule conditions. The mass-resolved, resonant two-photon ionization spectra of the two molecules were recorded in the region of the S0-S1 origin of the phenyl substituents (37,200-37,800 cm(-1)). UV-UV hole-burning spectroscopy was used to determine the ultraviolet spectral signatures of five conformational isomers of both 1 and 2. Transitions due to two conformers (labeled A and B) dominate the R2PI spectra of each molecule, while the other three are minor conformers (C-E) with transitions a factor of 3-5 smaller. Resonant ion-dip infrared spectroscopy was used to obtain single-conformation infrared spectra in the 3300-3700 cm(-1) region. The infrared spectra showed patterns of NH stretch transitions characteristic of the number and type of intramolecular H-bonds present in the beta-peptide backbone. For comparison with experiment, full optimizations of low-lying minima of both molecules were carried out at DFT B3LYP/6-31+G*, followed by single point MP2/6-31+G* and selected MP2/aug-cc-pVDZ calculations at the DFT optimized geometries. Calculated harmonic vibrational frequencies and infrared intensities for the amide NH stretch vibrations were used to determine the beta-peptide backbone structures for nine of the ten observed conformers. Conformers 1B, 1D, and 2A were assigned to double ring structures containing two C6 H-bonded rings (C6a/C6a), conformers 1A and 2B are C10 single H-bonded rings, conformers 1C and 2D are double ring structures composed of two C8 H-bonded rings (C8/C8), and conformers 1E and 2E are double ring/double acceptor structures in which two NH groups H-bond to the same C=O group, thereby weakening both H-bonds. Both 1E and 2E are tentatively assigned to C6/C8 double ring/double acceptor structures, although C8/C12 structures cannot be ruled out unequivocally. Finally, no firm conformational assignment has been made for conformer 2C whose unusual infrared spectrum contains one very strong H-bond with NH stretch frequency at 3309 cm(-1), a second H-bonded NH stretch fundamental of more typical value (3399 cm(-1)), and a third fundamental at 3440 cm(-1), below that typical of a branched-chain free NH. The single conformation spectra provide characteristic wavenumber ranges for the amide NH stretch fundamentals ascribed to C6 (3378-3415 cm(-1)), C8 (3339-3369 cm(-1)), and C10 (3381-3390 cm(-1)) H-bonded rings.     

Amide I vibrational frequencies of alpha-helical peptides based upon ONIOM and density functional theory (DFT) studies

We present ONIOM and pure DFT calculations on infrared spectra of alpha-helical-capped polyalanines. The calculations used two-layer ONIOM (B3LYP/D95**:AM1) calculations of the amide I vibrational frequencies for acetyl(ala)NNH2 (N=8, 10, 12-18) whose structures have been previously completely optimized by the same method. These are the first such calculations based upon structures of alpha-helical peptides that are completely optimized using DFT or molecular orbital methods. As the peptide becomes longer, the amide I band becomes both more intense and more red shifted. However, the individual absorptions that contribute most to the band vary between three patterns: one very intense absorption, two absorptions of similar intensity, and two strong absorptions where one is roughly twice as intense as the other. This pattern appears to be related to the relative number of H bonds in the individual H-bonding chains; however, there is one exception. Using 14C=O's to selectively decouple specific C=O's, we found that the couplings between the C=O's within each of the three individual H-bonding chains within the helices follow the same pattern previously reported for planar H-bonding chains of formamides. The coupling between the H-bonding chains appears to involve through-space coupling between the H-bonding chains. While decoupling individual C=O's always decreases the intensity of the amide I band, it leads to complex changes in the individual amide I absorptions that contribute to the band. Depending upon the position of the 14C=O, the amide I band can either red or blue shift. Moreover, the individual absorptions that contribute to the band can increase or decrease in intensity as well as shift. The patterns of the individual absorptions (mentioned above) also change. Using the C=O stretch of acetamide as a reference, we calculate the red shifts for the most intense absorptions to be much greater than predicted by the transition dipole method.