Simulation of amide I′ band profiles of trans polyproline based on an excitonic coupling model

Triglycine (GGG) Adopts a Polyproline II (pPII) Conformation in Its Hydrated Crystal Form: Revealing the Role of Water in Peptide Crystallization

Polyproline II (pPII) is a left-handed 31-helix conformation, which has been observed to be the most abundant secondary structure in unfolded peptides and proteins compared to α-helix and β-sheet. Although pPII has been reported as the most stable conformation for several unfolded short chain peptides in aqueous solution, it is rarely observed in their solid state. Here, we show for the first time a glycine homopeptide (gly-gly-gly) adopting the pPII conformation in its crystalline dihydrate structure. The single crystal X-ray structure with molecular dynamic simulation suggests that a network of water and the charged carboxylate group is critical in stabilizing the pPII conformation in solid state, offering an insight into the structures of unfolded regions of proteins and the role of water in peptide crystallization.

Bottom-Up Approach to Assess the Molecular Structure of Aqueous Poly(N-Isopropylacrylamide) at Room Temperature via Infrared Spectroscopy

The structure of poly(N-isopropylacrylamide) (PNIPAM) in solution is still an unresolved topic. Here, the PNIPAM structure in water was investigated using a bottom-up approach, involving the monomer, dimer, and trimer, and a combination of infrared (IR) spectroscopies as well as molecular dynamics simulations. The experiments show that the monomer and oligomers exhibit a broad and asymmetric amide I band with two underlying transitions, while PNIPAM presents the same major transitions and a minor one. Analysis of the 2D IR spectra and theoretical modeling of the amide I band indicates that the two transitions of the monomer do not have the same molecular origin as the oligomers and the polymer. In the monomer, the two bands originate from the ultrafast rotation of its ethyl group, which leads to different solvation structures for the various rotational conformers. In the case of the oligomers, the asymmetry and splitting of the amide I band is caused by the vibrational coupling among adjacent amide side chains. Moreover, it is deduced from the simulations that the oligomers have three distinct backbone conformations for neighboring amides. In particular, two of the backbone conformations have a closed and compact structure, while in the third, the backbone is open and elongated. The bottom-up approach allowed us to infer that such backbone conformations exist in PNIPAM as well. Consequently, the two major amide I transitions of the polymer are also assigned to split amide I transitions resulting from the vibrationally coupled nearest-neighboring amides. In contrast, the additional minor transition observed in PNIPAM is assigned to unsolvated amide units of the polymer. The proposed molecular model successfully describes that PNIPAM amide I band changes with temperature in terms of its molecular structure. This new model strongly suggests that PNIPAM does not have a completely random backbone structure, but has distinct backbone conformers between neighboring amides.

A crystal structure of an oligoproline PPII-helix, at last

The first crystal structure of an oligoproline adopting an all-trans polyproline II (PPII) helix is presented. The high-resolution structure provides detailed insight into the dimensions and conformational properties of oligoprolines that are important for, e.g., their use as "molecular rulers" and "molecular scaffolds". The structure also showed that the amides interact with each other within a PPII helix and that water is not necessary for PPII helicity.

Vibrational circular dichroism as a probe of fibrillogenesis: the origin of the anomalous intensity enhancement of amyloid-like fibrils

Amyloid fibrils are affiliated with various human pathologies. Knowledge of their molecular architecture is necessary for a detailed understanding of the mechanism of fibril formation. Vibrational circular dichroism (VCD) spectroscopy has recently shown sensitivity to amyloid fibrils [Ma et al. J. Am. Chem. Soc. 2007, 129, 12364 and Measey et al. J. Am. Chem. Soc. 2009, 131, 18218]. In particular, amyloid fibrils give rise to an intensity enhanced signal in the amide I band region of the corresponding VCD spectrum, offering promise of utilizing such a method for probing fibrillogenesis and the chiral structure of fibrils. Herein, we further investigate this phenomenon and demonstrate the use of VCD to probe the fibril formation kinetics of a short alanine-rich peptide. To elucidate the origin of the anomalous VCD intensity enhancement, we use an excitonic coupling model to simulate the VCD spectrum of stacked β-sheets containing one (Ising-like model) and two amide I oscillators per strand, as models for the underlying amyloid-fibril secondary structure. With this simple model, we show that the VCD intensity enhancement of amyloid-like fibrils results from intrasheet and, to a more limited extent, also from intersheet vibrational coupling between stacked β-sheets. The enhancement requires helically twisted sheets and is most pronounced for arrangements with parallel-oriented strands. Both the intersheet distance and the orientation of the amide I transition dipole moments of neighboring sheets are found to modulate the intensity enhancement of the amide I VCD signal. Moreover, our simulations suggest that, depending on the three-dimensional arrangement of the β-strands, the sign of the VCD signal of amyloid-like fibrils can be used to distinguish between right- and left-handed helical twists of parallel-oriented β-sheets. We compare the results of our simulation to experimental spectra of two short peptides, GNNQQNY, the N-terminal peptide fragment of the yeast prion protein Sup35, and an amyloidogenic alanine-rich peptide, AKY8. Our results demonstrate the advantages of using VCD spectroscopy to probe the kinetics of peptide and protein aggregation as well as the chirality of the resulting supramolecular structure.

Conformations of phenylalanine in the tripeptides AFA and GFG probed by combining MD simulations with NMR, FTIR, polarized Raman, and VCD spectroscopy

Conformational properties of small, flexible peptides are a matter of ongoing interest since they can be considered as models for unfolded proteins. However, the investigation of the conformations of small peptides is challenging as they are ensembles of rapidly interconverting conformers; moreover, the different methods used are prone to different approximations and errors. In order to obtain more reliable results, it is prudent to combine different techniques; here, molecular dynamics (MD) simulations together with nuclear magnetic resonance (NMR), Fourier transform IR (FTIR), polarized Raman, and vibrational circular dichroism (VCD) measurements were used to study the conformational propensity of phenylalanine in the tripeptides AFA and GFG, motivated by the relevance of phenylalanine for the self-aggregation of peptides. The results of this analysis indicate that the F residue predominantly populates the beta-strand (beta) and polyproline II (PPII) conformations in both AFA and GFG. However, while phenylalanine exhibits a propensity for beta-strand conformations in GFG (0.40 < or = beta population < or = 0.69 and 0.29 < or = PPII population < or = 0.42), the substitution of terminal glycines with alanine residues induces a higher population of PPII (0.31 < or = beta population < or = 0.50 and 0.37 < or = PPII population < or = 0.57).

Intrinsic propensities of amino acid residues in GxG peptides inferred from amide I' band profiles and NMR scalar coupling constants

A reliable intrinsic propensity scale of amino acid residues is indispensable for an assessment of how local conformational distributions in the unfolded state can affect the folding of peptides and proteins. Short host-guest peptides, such as GxG tripeptides, are suitable tools for probing such propensities. To explore the conformational distributions sampled by the central amino acid residue in these motifs, we combined vibrational (IR, Raman, and VCD) with NMR spectroscopy. The data were analyzed in terms of a superposition of two-dimensional Gaussian distribution functions in the Ramachandran space pertaining to subensembles of polyproline II, beta-strand, right- and left-handed helical, and gamma-turn-like conformations. The intrinsic propensities of eight amino acid residues (x = A, V, F, L, S, E, K, and M) in GxG peptides were determined as mole fractions of these subensembles. Our results show that alanine adopts primarily (approximately 80%) a PPII-like conformation, while valine and phenylalanine were found to sample PPII and beta-strand-like conformations equally. The centers of the respective beta-strand distributions generally do not coincide with canonical values of dihedral angles of residues in parallel or antiparallel beta-strands. In fact, the distributions for most residues found in the beta-region significantly overlap the PPII-region. A comparison with earlier reported results for trivaline reveals that the terminal valines increase the beta-strand propensity of the central valine residue even further. Of the remaining investigated amino acids, methionine preferred PPII the most (0.64), and E, S, L, and K exhibit moderate (0.56-0.45) PPII propensities. Residues V, F, S, E, and L sample, to a significant extent, a region between the canonical PPII and (antiparallel) beta-strand conformations. This region coincides with the sampling reported for L and V using theoretical predictions (Tran et al. Biochemistry 2005, 44, 11369). The distributions of all investigated residues differ from coil library and computationally predicted distributions in that they do not exhibit a substantial sampling of helical conformations. We conclude that this sampling of helical conformations arises from the context dependence, for example, neighboring residues, in proteins and longer peptides, some of which is long-range.

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.

FTIR studies of collagen model peptides: complementary experimental and simulation approaches to conformation and unfolding

X-ray crystallography of collagen model peptides has provided high-resolution structures of the basic triple-helical conformation and its water-mediated hydration network. Vibrational spectroscopy provides a useful bridge for transferring the structural information from X-ray diffraction to collagen in its native environment. The vibrational mode most useful for this purpose is the amide I mode (mostly peptide bond C=O stretch) near 1650 cm-1. The current study refines and extends the range of utility of a novel simulation method that accurately predicts the infrared (IR) amide I spectral contour from the three-dimensional structure of a protein or peptide. The approach is demonstrated through accurate simulation of the experimental amide I contour in solution for both a standard triple helix, (Pro-Pro-Gly)10, and a second peptide with a Gly --> Ala substitution in the middle of the chain that models the effect of a mutation in the native collagen sequence. Monitoring the major amide I peak as a function of temperature gives sharp thermal transitions for both peptides, similar to those obtained by circular dichroism spectroscopy, and the Fourier transform infrared (FTIR) spectra of the unfolded states were compared with polyproline II. The simulation studies were extended to model early stages of thermal denaturation of (Pro-Pro-Gly)10. Dihedral angle changes suggested by molecular dynamics simulations were made in a stepwise fashion to generate peptide unwinding from each end, which emulates the effect of increasing temperature. Simulated bands from these new structures were then compared to the experimental bands obtained as temperature was increased. The similarity between the simulated and experimental IR spectra lends credence to the simulation method and paves the way for a variety of applications.

Time-Domain Calculations of the Infrared and Polarized Raman Spectra of Tetraalanine in Aqueous Solution

The IR and polarized (isotropic and anisotropic) Raman spectra are calculated for the amide I band of tetraalanine ((Ala)4) in aqueous solution by using a time-domain computational method, which includes the effects of the diagonal frequency modulations (of individual peptide groups), the off-diagonal (interpeptide) vibrational couplings, and structural dynamics. It is shown that the calculated band profiles, especially the existence of a large negative noncoincidence effect (i.e., large frequency separations between the IR, isotropic Raman, and anisotropic Raman bands, with the isotropic Raman being higher in frequency), are in reasonable agreement with the experimental results. This negative noncoincidence effect derives from two conditions: the positive coupling between the amide I vibrations of peptide groups and the angle larger than 90 degrees between the transition dipoles of the coupled vibrations. This result means that the dynamically changing structures mainly in the polyproline II and beta-type conformations containing some repeated interconversions obtained from the molecular dynamics calculation are consistent with the existence of a large negative noncoincidence effect, as far as the structures satisfy the above two conditions. It is also shown that the electric fields from solvent water molecules induce larger frequency shifts than those of intrachain interactions, with rapid underdamped oscillatory modulations ( approximately 43 fs) due to the librational motions of water molecules that give rise to motional narrowing effect on the spectra. The reason for the difference from the behavior seen for the O-H stretching mode of liquid water is discussed. The time-domain analysis of the mode identity shows that the system proceeds halfway to complete mode mixing with a similar time scale ( approximately 60 fs), suggesting the importance of the nonadiabatic effect, which is included in a natural way in the present computational method.

The alanine-rich XAO peptide adopts a heterogeneous population, including turn-like and polyproline II conformations

The solution structure of the hepta-alanine polypeptide Ac-X(2)A(7)O(2)-NH(2) (XAO) has been a matter of controversy in the current literature. On one side of the argument is a claim that the peptide adopts a mostly polyproline II (PPII) structure, with a <20% population of beta conformations at room temperature [Shi Z, Olson CA, Rose GA, Baldwin RL, Kallenbach NR (2002) Proc Natl Acad Sci USA 99:9190-9195], whereas the other side of the argument insists that the peptide exists as an ensemble of conformations, including multiple beta-turn structures [Makowska J, Rodziewicz-Motowidlo S, Baginska K, Vila JA, Liwo A, Chmurzynski L, Scheraga HA (2006) Proc Natl Acad Sci USA 103:1744-1749]. We have used an excitonic coupling model to simulate the amide I band of the FTIR, vibrational circular dichroism, and isotropic and anisotropic Raman spectra of XAO, where, for each residue, the backbone dihedral angle varphi was constrained by using the reported (3)J(CalphaHNH) values and a modified Karplus relation. The best reproduction of the experimental data could only be achieved by assuming an ensemble of conformations, which contains various beta-turn conformations ( approximately 26%), in addition to beta-strand ( approximately 23%) and PPII ( approximately 50%) conformations. PPII is the dominant conformation in segments not involved in turn formations. Most of the residues were found to sample the bridge region connecting the PPII and right-handed helix troughs in the Ramachandran plot, which is part of the very heterogeneous ensemble of conformations generally termed type IV beta-turn.

Conformations of Alanine-Based Peptides in Water Probed by FTIR, Raman, Vibrational Circular Dichroism, Electronic Circular Dichroism, and NMR Spectroscopy

We have used a combination of FTIR, VCD, ECD, Raman, and NMR spectroscopies to probe the solution conformations sampled by H-(AAKA)-OH by utilizing an excitonic coupling model and constraints imposed by the 3JCalphaHNH coupling constants of the central residues to simulate the amide I' profile of the IR, isotropic Raman, anisotropic Raman, and VCD spectra in terms of a mixture of three conformations, i.e., polyproline II, beta-strand and right-handed helical. The representative coordinates of the three conformations were obtained from published coil libraries. Alanine was found to exhibit PPII fractions of 0.60 or greater, mixed with smaller fractions of helices and beta-strand conformations. Lysine showed no clear conformational propensity in that it samples polyproline II, beta-strand, and helical conformations with comparable probability. This is at variance with results obtained earlier for ionized polylysine, which suggest a high polyproline II propensity. We reanalyzed previously investigated tetra- and trialanine by combining published vibrational spectroscopy data with 3JCalphaHNH coupling constants and obtained again blends dominated by PPII with smaller admixtures of beta-strand and right-handed helical conformations. The polyproline II propensity of alanine was found to be higher in tetraalanine than in trialanine. For all peptides investigated, our results rule out a substantial population of turn-like conformations. Our results are in excellent agreement with MD simulations on short alanine peptides by Gnanakaran and Garcia [(2003) J. Phys. Chem. B 107, 12555-12557] but at variance with multiple MD simulations particularly for the alanine dipeptide.

Aggregation of the Amphipathic Peptides (AAKA)n into Antiparallel β-Sheets

Helical wheel projections of peptides based on the repeating unit Ac-(AAKA)n-NH2 clearly illustrate an amphipathic nature. One should therefore expect these peptides to form helices if the number of residues exceeds a certain threshold value. Indeed, ECD measurements show that Ac-(AAKA)4-NH2 and, to a minor extent, also Ac-(AAKA)3-NH2 exhibit some helical content at millimolar concentrations in aqueous solution. Surprisingly, however, these peptides were found to form hydrogels with an antiparallel beta-sheet conformation at centimolar concentrations. This occurs despite the positively charged lysine side chain which would be expected to inhibit the formation of extended beta-sheet layers.

Structure ofN-Acetylproline Amide in Liquid Water: Experimentally Measured and Numerically Simulated Infrared and Vibrational Circular Dichroism Spectra†

A few experimental and theoretical studies on the molecular structure of N-acetylproline amide (AP) in D2O solution have been reported recently. However, there is no consensus of the precise structure of AP in D2O because spectroscopically determined structures and a theoretically simulated one have been found to be different from one another. To determine its aqueous solution structure, IR and vibrational circular dichroism spectra of both L- and D-form AP solutions were measured. Molecular dynamics simulations with two different force fields and density functional theory calculations for the trans and cis rotamers of AP were performed to numerically simulate those spectra. Comparisons between experimentally measured and computationally simulated spectra directly suggest that the AP in water adopts a polyproline II-like conformation and that the force field parameter ff03 in the AMBER 8 suite of programs is more realistic and reliable in predicting molecular structure of AP in water than the ff99 in AMBER 7.

Salmon calcitonin and amyloid beta: two peptides with amyloidogenic capacity adopt different conformational manifolds in their unfolded states

The molecular conformations of salmon calcitonin in aqueous solution have been investigated by exploiting the different influences of excitonic coupling on the amide I band profile in the isotropic and anisotropic Raman, FTIR, and vibrational circular dichroism spectra of a polypeptide. The N-terminal loop, caused by a disulfide bridge between cysteines at positions 1 and 7, was modeled by performing a conformational search by molecular mechanics calculations. The remaining part of the peptide chain was modeled as a mixture of three sequences containing different fractions of residues adopting poly-l-proline II (PPII), extended beta-strand, and alpha-helix-like conformations. This yielded an excellent reproduction of the experimentally observed amide I' band profiles. A comparison with recent data on the beta-amyloid fragment Abeta(1)(-)(28) revealed a lower PPII content and more conformational heterogeneity for calcitonin. Thus, our results underscore the notion that individual structural propensities of amino acid residues give rise to structural differences between the unfolded states of even long peptide chains, at variance with expectations based on a random or statistical coil model.

Intermediacy of poly(L-proline) II and beta-strand conformations in poly(L-lysine) beta-sheet formation probed by temperature-jump/UV resonance Raman spectroscopy

Ultraviolet resonance Raman spectroscopy (UVRR) in combination with a nanosecond temperature jump (T-jump) was used to investigate early steps in the temperature-induced alpha-helix to beta-sheet conformational transition of poly(L-lysine) [poly(K)]. Excitation at 197 nm from a tunable frequency-quadrupled Ti:sapphire laser provided high-quality UVRR spectra, containing multiple conformation-sensitive amide bands. Although un-ionized poly(K) (pH 11.6) is mainly alpha-helical below 30 degrees C, there is a detectable fraction (approximately 15%) of unfolded polypeptide, which is mainly in the poly(L-proline) II (PPII) conformation. However, deviations from the expected amide I and II signals indicate an additional conformation, suggested to be beta-strand. Above 30 degrees C un-ionized poly(K) forms a beta-sheet at a rate (minutes) which increases with increasing temperature. A 22-44 degrees C T-jump is accompanied by prompt amide I and II difference signals suggested to arise from a rapid shift in the PPII/beta-strand equilibrium. These signals are superimposed on a subsequently evolving difference spectrum which is characteristic of PPII, although the extent of conversion is low, approximately 2% at the 3 micros time limit of the experiment. The rise time of the PPII signals is approximately 250 ns, consistent with melting of short alpha-helical segments. A model is proposed in which the melted PPII segments interconvert with beta-strand conformation, whose association through interstrand H-bonding nucleates the formation of beta-sheet. The intrinsic propensity for beta-strand formation could be a determinant of beta-sheet induction time, with implications for the onset of amyloid diseases.

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.