Simulations of oligopeptide vibrational CD: effects of isotopic labeling

Increased accuracy of vibrational circular dichroism calculations for isotopically labeled helical peptides

Vibrational circular dichroism (VCD) spectra have been computed with qualitatively correct sign patterns for α-helical peptides using various methods, ranging from empirical models to ab initio quantum mechanical computations. However, some details, such as deuteration effects and isotope substitution shifts and sign patterns for the resultant amide I' band shape, have remained a predictive challenge. Fully optimized computations for a 25-residue Ala-rich peptide, including implicit solvent corrections and explicit side chains that experimentally stabilize these model helical peptides in water, have been carried out using density functional theory (DFT). These fully minimized structures show minor changes in the (ϕ,ψ) torsions at the termini and yield an extra negative band to the low energy side of the characteristic amide I' couplet VCD, in agreement with experiments. Additionally, these calculations give the right sign and relative intensity patterns, as compared to experimental results, for several 13C=O substituted variants. The differences from previously reported computations that used ideal helical structures and vacuum conditions imply that inclusion of distorted termini and solvent effects can have an impact on the final detailed spectral patterns. Inclusion of side chains in these calculations had very little effect on the computed amide I' IR and VCD. Tests of constrained geometries, varying dielectric, and different functionals indicate that each can affect the band shapes, particularly for the 12C=O components, but these aspects do not fully explain the difference from previous spectral simulations. Inclusion of long-range amide coupling, as obtained from DFT computation of the full structure, or transfer of parameters from a somewhat longer peptide model, rather than shorter model, seems to be more important for the final detailed band shape under isotopic substitution. However, these corrections can also induce other changes, suggesting that previously reported, limited calculations may have been qualitatively useful due to a balance of errors. This may also explain the success of simple empirical IR models.

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.

Advances of Vibrational Circular Dichroism (VCD) in bioanalytical chemistry. A review

Vibrational Circular Dichroism (VCD) is a unique and relatively new spectroscopic technique that is capable of determining an absolute configuration of chiral molecules. VCD can be also used to determine structure of large macromolecules. This review highlights the most recent advances of VCD in bioanalytical chemistry. It shows that VCD is capable of unraveling supramolecular organization of peptides, proteins, saccharides, glycerophospholipids, polypeptide microcrystals, as well as amyloid fibrils and DNA. This review also demonstrates how VCD can be utilized to explore molecule-molecule interactions that determine mechanisms of chiral separations in chromatography. It aims to attract attention of scientists from all different research areas demonstrating the strength and capability of this very powerful spectroscopic technique.

PICVib: an accurate, fast and simple procedure to investigate selected vibrational modes and evaluate infrared intensities

The generalization of the PICVib approach for calculating selected infrared intensities is shown to be successful and to preserves its easiness of implementation and parallelization, flexibility and treatment of large systems and/or at high theoretical levels.

Infrared, vibrational circular dichroism, and Raman spectral simulations for β-sheet structures with various isotopic labels, interstrand, and stacking arrangements using density functional theory

Infrared (IR), Raman, and vibrational circular dichroism (VCD) spectral variations for different β-sheet structures were studied using simulations based on density functional theory (DFT) force field and intensity computations. The DFT vibrational parameters were obtained for β-sheet fragments containing nine-amides and constrained to a variety of conformations and strand arrangements. These were subsequently transferred onto corresponding larger β-sheet models, normally consisting of five strands with ten amides each, for spectral simulations. Further extension to fibril models composed of multiple stacked β-sheets was achieved by combining the transfer of DFT parameters for each sheet with dipole coupling methods for interactions between sheets. IR spectra of the amide I show different splitting patterns for parallel and antiparallel β-sheets, and their VCD, in the absence of intersheet stacking, have distinct sign variations. Isotopic labeling by (13)C of selected residues yields spectral shifts and intensity changes uniquely sensitive to relative alignment of strands (registry) for antiparallel sheets. Stacking of multiple planar sheets maintains the qualitative spectral character of the single sheet but evidences some reduction in the exciton splitting of the amide I mode. Rotating sheets with respect to each other leads to a significant VCD enhancement, whose sign pattern and intensity is dependent on the handedness and degree of rotation. For twisted β-sheets, a significant VCD enhancement is computed even for sheets stacked with either the same or opposite alignments and the inter-sheet rotation, depending on the sense, can either further increase or weaken the enhanced VCD intensity. In twisted, stacked structures (without rotation), similar VCD amide I patterns (positive couplets) are predicted for both parallel and antiparallel sheets, but different IR intensity distributions still enable their differentiation. Our simulation results prove useful for interpreting experimental vibrational spectra in terms of β-sheet and fibril structure, as illustrated in the accompanying paper.

Site-specific structure of Aβ(25-35) peptide: isotope-assisted vibrational circular dichroism study

We investigated the site-specific local structure of an amyloid peptide, NH(2)-GSNKGAIIGLM-COOH [Aβ(25-35)], one of the active fragments of amyloid β peptide that is known to be responsible for Alzheimer's disease, in the fibrillar aggregated state. Isotope-assisted infrared vibrational circular dichroism (VCD) and absorption (VA) spectroscopy were used for the parent Aβ(25-35) peptide, along with doubly (13)C labeled peptides at the carbonyl groups of residues 29 (Gly) and 30 (Ala) [Aβ(25-35:(13)C-29/30)] and at the carbonyl groups of residues 33 (Gly) and 34 (Leu) [Aβ(25-35:(13)C-33/34)]. The present results confirm that Aβ(25-35) peptide fibrils adopt a β-sheet structure and isotopic dilution experiments suggest a parallel β-sheet structure. The isotopic shifts suggest that the microenvironment of residues 29 (Gly) and 30 (Ala) could be different from that of residues 33 (Gly) and 34 (Leu). An unusual enhancement for the amide II' VCD intensities of Aβ(25-35:(13)C-29/30) and Aβ(25-35:(13)C-33/34) peptide fibrils, considered to originate from inter-strand coupling, was found for the first time. The structural information reported in this manuscript has important implications in understanding the role of this peptide in the development of Alzheimer's disease.

Impact ofβ-Turn Sequence onβ-Hairpin Dynamics Studied with Infrared-Detected Temperature Jump

Folding dynamics forβ-structure loss and disordered structure gain were studied in a modelβ-hairpin peptide based on Cochran’s tryptophan zipper peptide Trpzip2, but with an altered Thr-Gly (TG) turn sequence, that is, SWTWETGKWTWK, using laser-induced temperature-jump (T-jump) kinetics with IR detection. As has been shown previously, the TG turn sequence reduces the thermodynamicβ-hairpin stability as compared to the Asn-Gly sequence used in Trpzip2 (TZ2-NG). In this study, we found that the TG-turn slows down the overall relaxation dynamics as compared to TZ2-NG, which were studied at higher temperatures where the time constants show little difference between relaxation of theβ-strand and the disordered conformation. These time constants become equivalent at lower temperatures for TZ2-TG than was seen for TZ2-NG. The correlation of thermodynamic stability and rates of relaxation suggests that the change from NG to TG turn results in a slowing of folding, lowerkf, with less change of the unfolding rate,ku, assuming two state behavior at higher temperatures.

Isotope-assisted vibrational circular dichroism investigations of amyloid β peptide fragment, Aβ(16-22)

Isotope-assisted vibrational circular dichroism (VCD) investigations have been used to probe the site specific local structure of an amyloid peptide for the first time. A seven residue peptide, NH(2)-KLVFFAE-COOH, which represents the Aβ(16-22) fragment of the Alzheimer's amyloid β peptide, was used for these investigations. (13)C labels were introduced separately at the carbonyl group of leucine (residue 17), alanine (residue 21) and also at both sites together. Since VCD spectra provide structure dependent signs, band shapes and frequencies, the isotope-assisted VCD spectroscopy revealed information on site specific secondary structure of the polypeptide. Isotope dilution VCD experiments provided a means to distinguish between parallel and anti-parallel nature of the β-sheet structure formed by the Aβ(16-22) fragment. The current results establish the usefulness of isotope-assisted VCD analysis in determining the site specific secondary structure of amyloid peptides.

Applications of the Cartesian coordinate tensor transfer technique in the simulations of vibrational circular dichroism spectra of oligonucleotides

The application of the Cartesian coordinate tensor transfer (CCT) technique for simulations of the IR absorption and vibrational circular dichroism (VCD) spectra of relatively large nucleic acid fragments is demonstrated on several case studies. The approach is based on direct ab initio calculations of atomic tensors, determining molecular properties, for relatively small fragments, and subsequent transfer of these tensors to the larger systems in Cartesian coordinates. This procedure enables precise computations of vibrational spectra for large biomolecular systems, currently with up to several thousands of atoms. The versatile ability of the CCT methods is emphasized on the examples of VCD and IR absorption spectra calculations for B- and Z-forms of DNA, single-, double-, and triple-stranded RNA helices and DNA structures with different base content and sequences. The development and recent improvements of the methodology are followed, including utilization of the constrained normal mode optimization (NMO) strategy and combined quantum mechanics and molecular dynamics simulations. Advantages, drawbacks, and recommendations for future improvements of the CCT method as applied to nucleic acid spectra calculations are discussed.

Bioinspired interactions studied by vibrational circular dichroism

Vibrational circular dichroism (VCD) spectra are reliable indicators of the spatial structure of chiral molecules. The specific and characteristic feature of vibrational spectroscopy, and therefore also of VCD, where the energy of some vibrational modes is predominantly focused to a specific part of the molecule, enables monitoring both the structure of the molecule dissolved in different solvents and under different physicochemical conditions and molecular interactions. This minireview deals with recent contributions covering structural information on the bioinspired interactions obtained by means of VCD, especially in the following areas: interaction of DNA with biomolecules and biogenic metals, guanine tetramers and quadruplexes, biointeractions of bile pigments, and polypeptide and protein interactions with other biomolecules.

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.

Simulations of structure and vibrational spectra of deoxyoctanucleotides

Combined molecular dynamics and ab initio computations were applied for analysis of infrared absorption and vibrational circular dichroism spectra of deoxyoctanucleotides. Unlike for previous idealized models, molecular geometries of these shorter DNA fragments in solutions were obtained as dynamic averages from simulations in a periodic water box. Vibrational spectra for the whole octamers including hydrogen-bonded solvent molecules were simulated on the basis of density-functional calculations on small fragments and subsequent transfer of molecular property tensors. Explicit and continuum solvent models were compared. Apparently, the DNA segments retain an approximate B-conformation in the aqueous solutions, but the terminal base pairs significantly deviate from the planar arrangement and the vibrational circular dichroism spectrum for (CG)(2) nucleotide indicates a larger average helical twist. Sodium counterions moved freely around the molecule during the simulation and do not influence spectral intensities. Simulated absorption spectra faithfully reproduced the experimental signal of principal functional groups, while only qualitative agreement was obtained for the dependence on the basis sequence.

DNA oligonucleotide-cis-platin binding: Ab initio interpretation of the vibrational spectra

The cis-platin binding to the d(CCTGGTCC)*d(GGACCAGG) model DNA octamer was monitored with infrared absorption (IR) and vibrational circular dichroism (VCD) spectroscopies. The spectra were modeled with the aid of density functional computations and a Cartesian coordinate-based transfer of molecular property tensors from smaller DNA fragments. Because of the fragmentation, the tensors could be calculated with a higher precision. Environmental effects, such as the presence of the solvent or the cis-platin ligand, could be included in the modeling. The solvent was modeled by an explicit inclusion of hydrogen-bound water molecules, positions of which were estimated from a molecular dynamics simulation, or by the polarized continuum COSMO model. The B3LYP and BPW91 functionals used for the calculations of the spectral parameters were combined with the relativistic LANL2DZ platinum pseudo-potentials. The simulations reproduced the main IR and VCD DNA spectral features and explained most of the changes observed experimentally upon metal binding. The results confirmed that the influence of the ligand on DNA vibrational properties is quite complex; it originates in the geometry deformation and normal mode coupling pattern changes of the platinated octamer, as well as in local perturbations of the electronic structure and force field of the GC base pairs to which the platinum is bound. Many of the local effects could be accounted for by a point charge used in place of the metal in the GC complex.

Through Hydrogen-Bond Vibrational Coupling in Hydrogen-Bonding Chains of 4-Pyridones with Implications for Peptide Amide I Absorptions: Density Functional Theory Compared with Transition Dipole Coupling

We present B3LYP/D95** calculations on the C=O and N-H couplings in H-bonded chains of 4-pyridones. 14C-substitutions are used to decouple various vibrations for purposes of illustration. The coupled C=O vibrations bear analogy to the amide I bands of proteins and peptides. The coupling of the C=O's occurs primarily via the cooperative H-bonds rather than transition dipole coupling (TDC), as demonstrated by the fact that (1) the couplings are greater than previously reported for similar studies on formamides despite the larger distance between the C=O's in the pyridone chains (TDC coupling decreases with distance) and (2) the red shifts (also greater than for formamides) can be attributed to the changes in the geometries (particularly the C=O bond lengths) of the individual 4-pyridones in the H-bonding chains induced by the H-bonds and resulting polarization of the monomers.

Optical spectroscopic elucidation of beta-turns in disulfide bridged cyclic tetrapeptides

Vibrational circular dichroism (VCD) spectroscopic features of type II beta-turns were characterized previously, but, criteria for differentiation between beta-turn types had not been established yet. Model tetrapeptides, cyclized through a disulfide bridge, were designed on the basis of previous experimental results and the observed incidence of amino acid residues in the i + 1 and i + 2 positions in beta-turns, to determine the features of VCD spectra of type I and II beta-turns. The results were correlated with electronic circular dichroism (ECD) spectra and VCD spectra calculated from conformational data obtained by molecular dynamics (MD) simulations. All cyclic tetrapeptides yielded VCD signals with a higher frequency negative and a lower frequency positive couplet with negative lobes overlapping. MD simulations confirmed the conformational homogeneity of these peptides in solution. Comparison with ECD spectroscopy, MD, and quantum chemical calculation results suggested that the low frequency component of VCD spectra originating from the tertiary amide vibrations could be used to distinguish between types of beta-turn structures. On the basis of this observation, VCD spectroscopic features of type II and VIII beta-turns and ECD spectroscopic properties of a type VIII beta-turn were suggested. The need for independent experimental as well as theoretical investigations to obtain decisive conformational information was recognized.

Elucidation of residue-level structure and dynamics of polypeptides via isotope-edited infrared spectroscopy

Infrared spectroscopy is a powerful tool for analyzing the structure of proteins and peptides. The amide I band is particularly sensitive to the strength and position of the hydrogen bonds that define secondary structure as well as dipole-dipole interactions that are affected by the geometry of the peptide backbone. The introduction of a single (13)C-labeled carbonyl into a peptide backbone results in a resolvable shoulder to the main amide I band, which can be analyzed as a separate peak. Thus, site-specific structural information can be obtained by sequential, systematic labeling of the backbone. This method of isotope-edited infrared spectroscopy is a tool for obtaining medium-resolution information about the backbone conformation and dynamics. This tool has been used to dissect the conformation and dynamics of alpha helices and amyloid aggregates, where the versatility of possible sampling with infrared spectroscopy is well-suited for studies of large-protein aggregates.

Vibrational Spectral Simulation for Peptides of Mixed Secondary Structure: Method Comparisons with the Trpzip Model Hairpin

Infrared absorption and vibrational circular dichroism (IR and VCD) spectra of model fragments of TrpZip-style beta-hairpin structures are simulated using density functional theory (DFT) methods to estimate the influence of fragment size, end effects, conformational irregularities, peptide side chains, and solvent. Different fragmentation schemes, computing the strands and turn segments separately, were tested by varying the sizes of each and their respective overlaps. For suitably overlapping fragments, atomic property tensors were found to be reliably transferable, as tested by their ability to generate simulated spectra in good agreement with results from ab initio DFT computations for the entire peptide. This fragment approach significantly reduces computational times and opens up a wider range of systems that can be studied with a DFT-based approach as compared to previous methods based on uniform repeating sequences. However, vacuum calculations do not adequately represent the frequency dispersion of solvated molecules, and thus, some alternate strategies for solvation correction are explored for improving the simulation accuracy. Unlike for regular periodic secondary structure, the solvent significantly impacts the spectral shapes of hairpins, due to the different degrees of hydration of individual amide groups, which can be exposed to or shielded from water due to external vs internal hydrogen bonding. This is amplified by the shielding of selected amides from the solvent due to bulky side chains. The peptide plus solvent was structurally modeled with molecular dynamics methods, and then an electrostatic field-based parametrization correction was added to the force field and intensity tensors to compensate for the solvent dipolar field. The effect of the shielding and subsequent reordering of modes has a larger impact on VCD than IR band shapes.

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.

Coherent third-order spectroscopic probes of molecular chirality

The third-order optical response of a system of coupled localized anharmonic vibrations is studied using a Green's function solution of the nonlinear exciton equations for bosonized excitons, which are treated as interacting quasiparticles. The explicit calculation of two-exciton states is avoided and the scattering of quasiparticles provides the mechanism of optical nonlinearities. To first-order in the optical wave vector we find several rotationally invariant tensor components for isotropic ensembles which are induced by chirality. The nonlocal nonlinear susceptibility tensor is calculated for infinitely large periodic structures in momentum space, where the problem size reduces to the exciton interaction radius. Applications are made to alpha and 3(10) helical infinite peptides.

Solvent Effects on IR and VCD Spectra of Helical Peptides: DFT-Based Static Spectral Simulations with Explicit Water

Simulations of IR and VCD spectra are carried out for model alpha-helical, 3(10)-helical, and 3(1)-helical (polyProII-like) oligopeptides, with up to 21 amide groups, and including explicit consideration of effects of directly hydrogen-bonded solvent (water). Parameters used were obtained from ab initio density functional theory (DFT) computations of force field, atomic polar and axial tensors for oligopeptides of 5 to 7 amides, whose structures were constrained in (phi,psi) to target the secondary structure type but otherwise fully optimized. By comparison with experimental data as well as with calculations for identical but isolated (gas phase) peptides, the computed effects of an inner shell of aqueous solvent on the vibrational spectra of helical oligopeptides are illustrated. The interaction with solvent causes significant frequency shifts of the amide bands, but only minor changes in the characteristic IR intensity distributions and splittings and the VCD band shapes. Better agreement with experimental band shapes is achieved for the alpha-helical amide I' (N-deuterated) VCD by inclusion of explicit solvent in the calculations. Some improvements are also observed in theoretical VCD predictions for 13C labeled alpha-helical peptides when solvated models are used. However, the qualitative isotopic splitting patterns are preserved and just shifted in frequency due to consistent, solvent independent interamide coupling constants. The critical match of experiment and theory for relative positions of transitions in peptides with specifically separated 13C=O labels, including and neglecting solvent, confirms the stability of the coupling interactions. Despite these solvation effects, the calculated VCD band shape of the amide I mode is shown to be a reliable conformational probe, since it remains basically insensitive to frequency shifts caused by environment. Thus theoretical VCD simulations, even vacuum calculations, are shown to provide useful spectral predictions for solution-phase peptides.

Ab Initio Modeling of Amide I Coupling in Antiparallel β-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.

Circular dichroism spectroscopic study of non-covalent interactions of poly-L-glutamic acid with a porphyrin derivative in aqueous solutions

The interactions of poly-L-glutamic acid and a cationic porphyrin derivative in aqueous solutions were studied by the combination of vibrational circular dichroism (VCD) and electronic circular dichroism (ECD) spectroscopies. It was found that non-covalent interactions between both agents influence the structure of the polymeric matrix and the guest porphyrins and vice versa, but the physico-chemical properties of the solutions, especially the pH and the relative permittivity of the solvent, play a key role in the structure of the polypeptide part of the formed complexes. It was shown that the interaction with porphyrins prevents the precipitation of poly-L-glutamic acid in aqueous solution at acidic pH. In special conditions, the porphyrins attached to the polypeptide probably possess face-to-face interaction as demonstrated by the enhancement of the characteristic ECD signal and the appearance of sidebands on its short and long wavelength sides.

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.

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.

Optical spectroscopic investigations of model beta-sheet hairpins in aqueous solution

In this contribution we report optical spectroscopic data on a series of designed beta hairpins previously shown by NMR to contain a substantial population of beta-sheet structure. These models contain a designed hydrophobic cluster and a (D)Pro-Gly sequence to promote formation of a turn geometry. FTIR, electronic and vibrational CD (ECD and VCD) spectra for these small peptides are comparable to expected bandshapes for peptides of high beta-sheet content. The (D)Pro-Gly sequence provides a better turn motif than Asn-Gly as measured by its beta-sheet spectral characteristics. IR and VCD spectra are in qualitative agreement with theoretical simulations based on transfer of parameters from ab initio quantum mechanical force field and intensity computations for the turn and strands. These calculations provide assignments for some distinguishing modes in both IR and VCD spectra. Increased sheet structure can be induced in these hairpins by use of mixed solvent conditions. Thermal denaturation studies reveal that these hairpins undergo very broad unfolding transitions. Guanidine hydrochloride unfolding transitions for the selected hairpin models are similarly broad. However, the "end-states" of temperature and chaotropic denaturation are spectroscopically differentiable.

Discriminating 3(10)- from alpha-helices: vibrational and electronic CD and IR absorption study of related Aib-containing oligopeptides

Model peptides based on -(Aib-Ala)(n)-, and (Aib)(n)-Leu-(Aib)(2) sequences, which have varying amounts of 3(10)-helical character, were studied by use of vibrational and electronic circular dichroism (VCD and ECD) and Fourier transform infrared (FTIR) absorption spectroscopies to test the correlation of spectral response and conformation. The data indicate that these peptides, starting from a length of about four to six residues, predominantly adopt a 3(10)-helical conformation at room temperature. The longest model peptides, depending on the series, may evidence some alpha-helical contribution to the spectra, while the shorter ones, with less than six residues, have much less order. The IR absorption spectra (as supported by theory) showed only small frequency changes between 3(10)- and alpha-helices. By contrast, solvent effects are a source of much bigger perturbations. The ECD results show that the intensity ratio for the approximately 222-nm to approximately 208-nm bands, while useful for distinguishing between these two helical types in some sequences, may have a narrower range of application than VCD. However, the VCD data presented here continue to support the proposed discrimination between alpha- and 3(10)-helices based on qualitative amide I and II bandshape differences. The present study shows the intensities of the 3(10)-helical amide I (peak-to-peak) to its amide II VCD to be of the same order and useful for discriminating them from alpha-helices, whose amide I dominates the amide II in intensity. This qualitative result is experimentally independent of the amount of alphaMe-substituted residues in the sequence. These experimental VCD results are consistent in detail with theoretical spectral simulations for Ac-(Ala)(8)-NH(2), Ac-(Aib-Ala)(4)-NH(2), and Ac-(Aib)(8)-NH(2) in 3(10)- and alpha-helical conformations.