IR Study of Cross-Strand Coupling in a β-Hairpin Peptide Using Isotopic Labels

Chasing weakly-bound biological water in aqueous environment near the peptide backbone by ultrafast 2D infrared spectroscopy

There has been a long-standing debate as to how many hydrogen bonds a peptide backbone amide can form in aqueous solution. Hydrogen-bonding structural dynamics of N-ethylpropionamide (a β-peptide model) in water was examined using infrared (IR) spectroscopy. Two amide-I sub bands arise mainly from amide C=O group that forms strong H-bonds with solvent water molecules (SHB state), and minorly from that involving one weak H-bond with water (WHB state). This picture is supported by molecular dynamics simulations and ab-initio calculations. Further, thermodynamics and kinetics of the SHB and WHB species were examined mainly by chemical-exchange two-dimensional IR spectroscopy, yielding an activation energy for the SHB-to-WHB exchange of 13.25 ± 0.52 kJ mol‒1, which occurs in half picosecond at room temperature. Our results provided experimental evidence of an unstable water molecule near peptide backbone, allowing us to gain more insights into the dynamics of the protein backbone hydration.

Frequency Changes in Terminal Alkynes Provide Strong, Sensitive, and Solvatochromic Raman Probes of Biochemical Environments

The C≡C stretching frequencies of terminal alkynes appear in the "clear" window of vibrational spectra, so they are attractive and increasingly popular as site-specific probes in complicated biological systems like proteins, cells, and tissues. In this work, we collected infrared (IR) absorption and Raman scattering spectra of model compounds, artificial amino acids, and model proteins that contain terminal alkyne groups, and we used our results to draw conclusions about the signal strength and sensitivity to the local environment of both aliphatic and aromatic terminal alkyne C≡C stretching bands. While the IR bands of alkynyl model compounds displayed surprisingly broad solvatochromism, their absorptions were weak enough that alkynes can be ruled out as effective IR probes. The same solvatochromism was observed in model compounds' Raman spectra, and comparisons to published empirical solvent scales (including a linear regression against four meta-aggregated solvent parameters) suggested that the alkyne C≡C stretching frequency mainly reports on local electronic interactions (i.e., short-range electron donor-acceptor interactions) with solvent molecules and neighboring functional groups. The strong solvatochromism observed here for alkyne stretching bands introduces an important consideration for Raman imaging studies based on these signals. Raman signals for alkynes (especially those that are π-conjugated) can be exceptionally strong and should permit alkynyl Raman signals to function as probes at very low concentrations, as compared to other widely used vibrational probe groups like azides and nitriles. We incorporated homopropargyl glycine into a transmembrane helical peptide via peptide synthesis, and we installed p-ethynylphenylalanine into the interior of the Escherichia coli fatty acid acyl carrier protein using a genetic code expansion technique. The Raman spectra from each of these test systems indicate that alkynyl C≡C bands can act as effective and unique probes of their local biomolecular environments. We provide guidance for the best possible future uses of alkynes as solvatochromic Raman probes, and while empirical explanations of the alkyne solvatochromism are offered, open questions about its physical basis are enunciated.

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

Computing the frequency fluctuation dynamics of highly coupled vibrational transitions using neural networks

The description of frequency fluctuations for highly coupled vibrational transitions has been a challenging problem in physical chemistry. In particular, the complexity of their vibrational Hamiltonian does not allow us to directly derive the time evolution of vibrational frequencies for these systems. In this paper, we present a new approach to this problem by exploiting the artificial neural network to describe the vibrational frequencies without relying on the deconstruction of the vibrational Hamiltonian. To this end, we first explored the use of the methodology to predict the frequency fluctuations of the amide I mode of N-methylacetamide in water. The results show good performance compared with the previous experimental and theoretical results. In the second part, the neural network approach is used to investigate the frequency fluctuations of the highly coupled carbonyl stretch modes for the organic carbonates in the solvation shell of the lithium ion. In this case, the frequency fluctuation predicted by the neural networks shows a good agreement with the experimental results, which suggests that this model can be used to describe the dynamics of the frequency in highly coupled transitions.

Template-assisted design of monomeric polyQ models to unravel the unique role of glutamine side chains in disease-related aggregation

PolyQ model peptides reveal the effect of individual glutamine side chains on fibril formation.

Expanded polyglutamine (polyQ) sequences cause numerous neurodegenerative diseases which are accompanied by the formation of polyQ fibrils. The unique role of glutamines in the aggregation onset is undoubtedly accepted and a lot structural data of the fibrils have been acquired, however side-chain specific structural dynamics inducing oligomerization are not well understood yet. To analyze spectroscopically the nucleation process, we designed various template-assisted glutamine-rich β-hairpin monomers mimicking the structural motif of a polyQ fibril. In a top-down strategy, we use a template which forms a well-defined stable hairpin in solution, insert polyQ-rich sequences into each strand and monitor the effects of individual glutamines by NMR, CD and IR spectroscopic approaches. The design was further advanced by alternating glutamines with other amino acids (T, W, E, K), thereby enhancing the solubility and increasing the number of cross-strand interacting glutamine side chains. Our spectroscopic studies reveal a decreasing hairpin stability with increased glutamine content and demonstrate the enormous impact of only a few glutamines – far below the disease threshold – to destabilize structure. Furthermore, we could access sub-ms conformational dynamics of monomeric polyQ-rich peptides by laser-excited temperature-jump IR spectroscopy. Both, the increased number of interacting glutamines and higher concentrations are key parameters to induce oligomerization. Concentration-dependent time-resolved IR measurements indicate an additional slower kinetic phase upon oligomer formation. The here presented peptide models enable spectroscopic molecular analyses to distinguish between monomer and oligomer dynamics in the early steps of polyQ fibril formation and in a side-chain specific manner.

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.

Enhanced Sensitivity to Local Dynamics in Peptides by Use of Temperature-Jump IR Spectroscopy and Isotope Labeling

Site-specific isotopic labeling of molecules is a widely used approach in IR spectroscopy to resolve local contributions to vibrational modes. The induced frequency shift of the corresponding IR band depends on the substituted masses, as well as on hydrogen bonding and vibrational coupling. The impact of these different factors was analyzed with a designed three-stranded β-sheet peptide and by use of selected ¹³ C isotope substitutions at multiple positions in the peptide backbone. Single-strand labels give rise to isotopically shifted bands at different frequencies, depending on the specific sites; this demonstrates sensitivity to the local environment. Cross-strand double- and triple-labeled peptides exhibited two resolved bands that could be uniquely assigned to specific residues, the equilibrium IR spectra of which indicated only weak local-mode coupling. Temperature-jump IR laser spectroscopy was applied to monitor structural dynamics and revealed an impressive enhancement of the isotope sensitivity to both local positions and coupling between them, relative to that of equilibrium FTIR spectroscopy. Site-specific relaxation rates were altered upon the introduction of additional cross-strand isotopes. Likewise, the rates for the global β-sheet dynamics were affected in a manner dependent on the distinct relaxation behavior of the labeled oscillator. This study reveals that isotope labels provide not only local structural probes, but rather sense the dynamic complexity of the molecular environment.

Implications of aromatic-aromatic interactions: From protein structures to peptide models

With increasing structural information on proteins, the opportunity to understand physical forces governing protein folding is also expanding. One of the significant non-covalent forces between the protein side chains is aromatic-aromatic interactions. Aromatic interactions have been widely exploited and thoroughly investigated in the context of folding, stability, molecular recognition, and self-assembly processes. Through this review, we discuss the contribution of aromatic interactions to the activity and stability of thermophilic, mesophilic, and psychrophilic proteins. Being hydrophobic, aromatic amino acids tend to reside in the protein hydrophobic interior or transmembrane segments of proteins. In such positions, it can play a diverse role in soluble and membrane proteins, and in α-helix and β-sheet stabilization. We also highlight here some excellent investigations made using peptide models and several approaches involving aryl-aryl interactions, as an increasingly popular strategy in protein and peptide engineering. A recent survey described the existence of aromatic clusters (trimer, tetramer, pentamer, and higher order assemblies), revealing the self-associating property of aryl groups, even in folded protein structures. The application of this self-assembly of aromatics in the generation of modern bionanomaterials is also discussed.

Temperature dependence of C-terminal carboxylic group IR absorptions in the amide I' region

Studies of structural changes in peptides and proteins using IR spectroscopy often rely on subtle changes in the amide I' band as a function of temperature. However, these changes can be obscured by the overlap with other absorptions, namely the side-chain and terminal carboxylic groups. The former were the subject of our previous report (Anderson et al., 2014). In this paper we investigate the IR spectra of the asymmetric stretch of α-carboxylic groups for amino acids representing all major types (Gly, Ala, Val, Leu, Ser, Thr, Asp, Glu, Lys, Asn, His, Trp, Pro) as well as the C-terminal groups of three dipeptides (Gly-Gly, Gly-Ala, Ala-Gly) in D₂O at neutral pH. Experimental temperature dependent IR spectra were analyzed by fitting of both symmetric and asymmetric pseudo-Voigt functions. Qualitatively the spectra exhibit shifts to higher frequency, loss in intensity and narrowing with increased temperature, similar to that observed previously for the side-chain carboxylic groups of Asp. The observed dependence of the band parameters (frequency, intensity, width and shape) on temperature is in all cases linear: simple linear regression is therefore used to describe the spectral changes. The spectral parameters vary between individual amino acids and show systematic differences between the free amino acids and dipeptides, particularly in the absolute peak frequencies, but the temperature variations are comparable. The relative variations between the dipeptide spectral parameters are most sensitive to the C-terminal amino acid, and follow the trends observed in the free amino acid spectra. General rules for modeling the α-carboxylic IR absorption bands in peptides and proteins as the function of temperature are proposed.

The effects of alpha-helical structure and cyanylated cysteine on each other

β-Thiocyanatoalanine, or cyanylated cysteine, is an artificial amino acid that can be introduced at solvent-exposed cysteine residues in proteins via chemical modification. Its facile post-translational synthesis means that it may find broad use in large protein systems as a probe of site-specific structure and dynamics. The C≡N stretching vibration of this artificial side chain provides an isolated infrared chromophore. To test both the perturbative effect of this side chain on local secondary structure and its sensitivity to structural changes, three variants of a model water-soluble alanine-repeat helix were synthesized containing cyanylated cysteine at different sites. The cyanylated cysteine side chain is shown to destabilize, but not completely disrupt, the helical structure of the folded peptide when substituted for alanine. In addition, the C≡N stretching bandwidth of the artificial side chain is sensitive to the helix−coil structural transition. These model system results indicate that cyanylated cysteine can be placed into protein sequences with a native helical propensity without destroying the helix, and further that the CN probe may be able to report local helix formation events even when it is water-exposed in both the ordered and disordered conformational states. These results indicate that cyanylated cysteine could be a widely useful probe of structure-forming events in proteins with large in vitro structural distributions.

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.

Carbon-deuterium vibrational probes of peptide conformation: alanine dipeptide and glycine dipeptide

The utility of alpha-carbon deuterium-labeled bonds (C(alpha)-D) as infrared reporters of local peptide conformation was investigated for two model dipeptide compounds: C(alpha)-D labeled alanine dipeptide (Adp-d(1)) and C(alpha)-D(2) labeled glycine dipeptide (Gdp-d(2)). These model compounds adopt structures that are analogous to the motifs found in larger peptides and proteins. For both Adp-d(1) and Gdp-d(2), we systematically mapped the entire conformational landscape in the gas phase by optimizing the geometry of the molecule with the values of phi and psi, the two dihedral angles that are typically used to characterize the backbone structure of peptides and proteins, held fixed on a uniform grid with 7.5 degrees spacing. Since the conformations were not generally stationary states in the gas phase, we then calculated anharmonic C(alpha)-D and C(alpha)-D(2) stretch transition frequencies for each structure. For Adp-d(1) the C(alpha)-D stretch frequency exhibited a maximum variability of 39.4 cm(-1) between the six stable structures identified in the gas phase. The C(alpha)-D(2) frequencies of Gdp-d(2) show an even more substantial difference between its three stable conformations: there is a 40.7 cm(-1) maximum difference in the symmetric C(alpha)-D(2) stretch frequencies and an 81.3 cm(-1) maximum difference in the asymmetric C(alpha)-D(2) stretch frequencies. Moreover, the splitting between the symmetric and asymmetric C(alpha)-D(2) stretch frequencies of Gdp-d(2) is remarkably sensitive to its conformation.

Cross-strand coupling and site-specific unfolding thermodynamics of a trpzip beta-hairpin peptide using 13C isotopic labeling and IR spectroscopy

Conformational properties of a 12-residue tryptophan zipper (trpzip) beta-hairpin peptide (AWAWENGKWAWK-NH(2), a modification of the original trpzip2 sequence) are analyzed under equilibrium conditions using ECD and IR spectra of a series of variants, singly and doubly C(1)-labeled with (13)C on the amide CO. The characteristic features of the (13)CO component of the amide I' IR band and their sensitivity to the local structure of the peptide are used to differentiate stabilities for parts of the hairpin structure. Doubly labeled peptide spectra indicate that the ends of the beta-strands are frayed and that the center part is more stable as would be expected from formation of a stable hydrophobic core consisting of four tryptophan residues, and supported by MD simulations. NMR analyses were used to determine a best fit solution structure that is in close agreement with that of trpzip2, except for a small variation in the turn geometry. The distinct vibrational coupling patterns of the labeled sites based on this structure are also well matched by ab initio DFT-level calculations of their IR spectral patterns. Thermal unfolding of the peptides as studied with CD spectra could be fit with an apparent two-state transition model. ECD senses only the tryptophan interactions (tertiary-like) and their overall environment, as shown by TD-DFT modeling of the Trp-Trp pi-pi ECD. However, variation of the amide I IR spectra of (13)C-isotopomers showed that the thermal unfolding process is not cooperative in terms of the peptide backbone (secondary structure), since the transition temperatures sensed for labeled modes differ from those for the whole peptide. The thermal data also evidence dependence on concentration and pH but these cause little spectral variation. This study illustrates the consequences of multistate conformational change at the residue- or sequence-specific level in a system whose structure is dominated by hydrophobic collapse.

Very short peptides with stable folds: building on the interrelationship of Trp/Trp, Trp/cation, and Trp/backbone-amide interaction geometries

By combining a favorable turn sequence with a turn flanking Trp/Trp interaction and a C-terminal H-bonding interaction between a backbone amide and an i-2 Trp ring, a particularly stable (DeltaG(U) > 7 kJ/mol) truncated hairpin, Ac-WI-(D-Pro-D-Asn)-KWTG-NH(2), results. In this construct and others with a W-(4-residue turn)-W motif in severely truncated hairpins, the C-terminal Trp is the edge residue in a well-defined face-to-edge (FtE) aryl/aryl interaction. Longer hairpins and those with six-residue turns retain the reversed "edge-to-face" (EtF) Trp/Trp geometry first observed for the trpzip peptides. Mutational studies suggest that the W-(4-residue turn)-W interaction provides at least 3 kJ/mol of stabilization in excess of that due to the greater beta-propensity of Trp. The pi-cation, and Trp/Gly-H(N) interactions have been defined. The latter can give rise to >3 ppm upfield shifts for the Gly-H(N) in -WX(n)G- units both in turns (n = 2) and at the C-termini (n = 1) of hairpins. Terminal YTG units result in somewhat smaller shifts (extrapolated to 2 ppm for 100% folding). In peptides with both the EtF and FtE W/W interaction geometries, Trp to Tyr mutations indicate that Trp is the preferred "face" residue in aryl/aryl pairings, presumably because of its greater pi basicity.

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.

Diacid linkers that promote parallel beta-sheet secondary structure in water

We report the development of diacid units that promote formation of a two-stranded parallel beta-sheet secondary structure between peptide segments attached via their N-termini. These linker units are formed by attaching glycine to one carboxyl group of cis-1,2-cyclohexanedicarboxylic acid (CHDA). Parallel sheet formation in water is observed when l-residue strands are attached to the CHDA-Gly unit with either of the two absolute configurations.

Conformational dependence of anharmonic vibrations in peptides: amide-I modes in model dipeptide

The conformational dependence of a set of anharmonic vibrational parameters for the amide-I modes in peptides has been examined at the level of Hartree-Fock theory. By using glycine dipeptide as a model molecule, the phi-psi maps of diagonal and off-diagonal anharmonicities, local mode mixing angle, zero-order local mode frequencies, as well as intermode coupling, were calculated, and all were found to exhibit certain conformational sensitivities. The characteristics of the phi-psi maps of the diagonal and off-diagonal anharmonicities were found to be complementary to each other, and the latter was found to correlate well with that of the mixing angle, reflecting the fact that these anharmonic parameters are interconnected and all determined by the same set of underlying anharmonic force field. The mean values of the diagonal and off-diagonal anharmonicities were found to be 15.7 cm(-1) and 9.9 cm(-1) respectively. The mean value for the two local mode frequency differences was estimated to be 9.7 cm(-1), showing a nondegenerate local mode picture. For significant peptide conformations, the calculated anharmonic parameters were found to be in reasonable agreement with values obtained at the level of density functional theory as well as values obtained with recent two-dimensional infrared experiments.

Site-specific relaxation kinetics of a tryptophan zipper hairpin peptide using temperature-jump IR spectroscopy and isotopic labeling

Two antiparallel beta-strands connected by a turn make beta-hairpins an ideal model system to analyze the interactions and dynamics of beta-sheets. Site-specific conformational dynamics were studied by temperature-jump IR spectroscopy and isotopic labeling in a model based on the tryptophan zipper peptide, Trpzip2, developed by Cochran et al. (Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 5578). The modified Trpzip2C peptides have nearly identical equilibrium spectral behavior as Trpzip2 showing that they also form well-characterized beta-hairpin conformations in aqueous solution. Selective introduction of 13C=O groups on opposite strands lead to distinguishable cross-strand coupling of the labeled residues as monitored in the amide I' band. These frequency patterns reflect theoretical predictions, and the coupled 13C=O band loses intensity with increase in temperature and unfolding of the hairpin. Thermal relaxation kinetics were analyzed for unlabeled and cross-strand isotopically labeled variants. T-jumps of approximately 10 degrees C induce relaxation times of a few microseconds that decrease with increase of the peptide temperature. Differences in kinetic behavior for the loss of beta-strand and gain of disordered structure can be used to distinguish localized structure dynamics by comparison of nonlabeled and labeled amide I' components. Analysis of the data supports multistate dynamic and equilibrium behavior, but because of this process it is not possible to clearly define a folding and unfolding rate. Nonetheless, site-specific relaxation kinetics could be seen to be consistent with a hydrophobic collapse hypothesis for hairpin folding.

Nonstereogenic alpha-aminoisobutyryl-glycyl dipeptidyl unit nucleates type I' beta-turn in linear peptides in aqueous solution

The use of alpha,alpha-disubstituted amino acids represents a valuable strategy to exercise conformational control in peptides. Incorporation of the nonstereogenic alpha-aminoisobutyryl-glycyl (Aib-Gly) dipeptidyl sequence into i+1 and i+2 positions of an acyclic peptide sequence, originally designed and investigated by Gellman and coworkers, [H-Arg-Tyr-Val-Glu-Val-Yyy-Xxx-Orn-Lys-Ile-Leu-Gln-NH2] nucleates a stable [2:4] left-handed type I' beta-turn in water. NMR spectra show that this newly designed beta-hairpin does not aggregate in water up to a concentration of approximately 1 mM, and that its backbone conformation is superimposable on corresponding hairpins containing the DPro-Gly (literature) and Aib-DAla (this work) sequences. The Aib-Gly turn-inducer sequence eliminates complications because of cis-trans isomerization of Zzz-Pro bonds, and constitutes an attractive alternative to the proteogenic Asn-Gly and nonproteogenic DPro-Gly motifs previously suggested as turn-inducer sequences. These design principles could be exploited to prepare water-soluble beta-hairpin peptides with robust structures and novel function.

Spectrum–Structure Relationship in Carbohydrate Vibrational Circular Dichroism and Its Application to Glycoconjugates

Preliminary reports of the nature of the vibrational circular dichroism (VCD) peak at around 1145 cm(-1), which is characteristic of axial glycosidic sugars and is called the glycoside band (J. Am. Chem. Soc. 2004, 126, 9496), have been thoroughly examined. Through systematic carbohydrate measurements, it was found that the sign of the glycoside band reflects not only the anomeric configuration but also the pyranose conformation. Isotope and theoretical studies characterized its vibrational mode as C1-H1 deformation coupled with C1-O1 stretching, which indicates its applicability to more-complicated glycoconjugates. In this study, for the first time, carbohydrate VCD spectra were reliably predicted by means of density functional theory (DFT) calculations. The VCD technique was applied to glycopeptides, and simultaneous analysis of both the carbohydrate and aglycan parts was carried out.

Minimization and optimization of designed beta-hairpin folds

Minimized beta hairpins have provided additional data on the geometric preferences of Trp interactions in TW-loop-WT motifs. This motif imparts significant fold stability to peptides as short as 8 residues. High-resolution NMR structures of a 16- (KKWTWNPATGKWTWQE, DeltaG(U)(298) >or= +7 kJ/mol) and 12-residue (KTWNPATGKWTE, DeltaG(U)(298) = +5.05 kJ/mol) hairpin reveal a common turn geometry and edge-to-face (EtF) packing motif and a cation-pi interaction between Lys(1) and the Trp residue nearest the C-terminus. The magnitude of a CD exciton couplet (due to the two Trp residues) and the chemical shifts of a Trp Hepsilon3 site (shifted upfield by 2.4 ppm due to the EtF stacking geometry) provided near-identical measures of folding. CD melts of representative peptides with the -TW-loop-WT- motif provided the thermodynamic parameters for folding, which reflect enthalpically driven folding at laboratory temperatures with a small DeltaC(p) for unfolding (+420 J K(-)(1)/mol). In the case of Asx-Pro-Xaa-Thr-Gly-Xaa loops, mutations established that the two most important residues in this class of direction-reversing loops are Asx and Gly: mutation to alanine is destabilizing by about 6 and 2 kJ/mol, respectively. All indicators of structuring are retained in a minimized 8-residue construct (Ac-WNPATGKW-NH(2)) with the fold stability reduced to DeltaG(U)(278) = -0.7 kJ/mol. NMR and CD comparisons indicate that -TWXNGKWT- (X = S, I) sequences also form the same hairpin-stabilizing W/W interaction.

An infrared spectroscopic study of the conformational transition of elastin-like polypeptides

The infrared spectroscopy of elastin-like polypeptides and the relation to the inverse thermal transition are discussed. To correlate the spectroscopic observations with structure a density function theory model was created that captures the essential hydrogen bonding and packing of the beta-spiral structure proposed for elastin and elastin-like polypeptides. The infrared spectrum was calculated using periodic boundary conditions and a method for estimating the difference dipole moment permits both frequencies and intensities to be obtained for the modeling of spectra. The two observed amide I bands at 1615 cm(-1) and 1656 cm(-1) are shown to arise from the beta-spiral structure. The increase in intensity of these bands with increasing salt concentration and temperature is assigned to the closer association of strands of the beta-spiral. The sharp inverse temperature transition is observed within 1 degrees C and involves a change in secondary structure that involves formation of interstrand beta-sheets for approximately 25% of the amino acids. This conclusion is consistent with available data and simulations that have been reported to date.

Amide I two-dimensional infrared spectroscopy of beta-hairpin peptides

In this report, spectral simulations and isotope labeling are used to describe the two-dimensional IR spectroscopy of beta-hairpin peptides in the amide I spectral region. 2D IR spectra of Gramicidin S, PG12, Trpzip2 (TZ2), and TZ2-T3(*)T10(*), a dual (13)C(') isotope label, are qualitatively described by a model based on the widely used local mode amide I Hamiltonian. The authors' model includes methods for calculating site energies for individual amide oscillators on the basis of hydrogen bonding, nearest neighbor and long-range coupling between sites, and disorder in the site energy. The dependence of the spectral features on the peptide backbone structure is described using disorder-averaged eigenstates, which are visualized by mapping back onto the local amide I sites. beta-hairpin IR spectra are dominated by delocalized vibrations that vary by the phase of adjacent oscillators parallel and perpendicular to the strands. The dominant nu(perpendicular) band is sensitive to the length of the hairpin and the amount of twisting in the backbone structure, while the nu(parallel) band is composed of several low symmetry modes that delocalize along the strands. The spectra of TZ2-T3(*)T10(*) are used to compare coupling models, from which we conclude that transition charge coupling is superior to transition dipole coupling for amide groups directly hydrogen bound across the beta strands. The 2D IR spectra of TZ2-T3(*)T10(*) are used to resolve the redshifted amide I band and extract the site energy of the labeled groups. This allows the authors to compare several methods for calculating the site energies used in excitonic treatments of the amide I band. Gramicidin S is studied in dimethyl sulfoxide to test the role of solvent on the spectral simulations.

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.