Identifying Residual Structure in Intrinsically Disordered Systems: A 2D IR Spectroscopic Study of the GVGXPGVG Peptide

Structural and Dynamical Properties of Elastin-Like Peptides near Their Lower Critical Solution Temperature

Elastin-like peptides (ELPs) are artificially derived intrinsically disordered proteins (IDPs) mimicking the hydrophobic repeat unit in the protein elastin. ELPs are characterized by a lower critical solution temperature (LCST) in aqueous media. Here, we investigate the sequence GVG(VPGVG)₃ over a wide range of temperatures (below, around, and above the LCST) and peptide concentrations employing all-atom molecular dynamics simulations, where we focus on the role of intra- and interpeptide interactions. We begin by investigating the structural properties of a single peptide that demonstrates a hydrophobic collapse with temperature, albeit moderate, because the sequence length is short. We observe a change in the interaction between two peptides from repulsive to attractive with temperature by evaluating the potential of mean force, indicating an LCST-like behavior. Next, we explore dynamical and structural properties of peptides in multichain systems. We report the formation of dynamical aggregates with coil-like conformation, in which valine central residues play an important role. Moreover, the lifetime of contacts between chains strongly depends on the temperature and can be described by a power-law decay that is consistent with the LCST-like behavior. Finally, the peptide translational and internal motion are slowed by an increase in the peptide concentration and temperature.

Evaluating and improving the reliability of the UV-persulfate method for the determination of TOC/DOC in surface waters

The UV-persulfate oxidation method is widely used for determining the total organic carbon concentration of aqueous samples (denoted for convenience as UVP-TOC). However, for some surface water samples, the measurement of TOC by this method can be unreliable, deviating significantly from the true carbon content. In this study, the performance of the UVP-TOC method has been investigated by comparing the results from the analysis of a variety of aqueous samples that included two kinds of surface water samples and related surface water model substances: bovine serum albumin (BSA), sodium alginate (SA), humic acid (HA), tannic acid (TA), benzoic acid (BA) and citric acid (CA), with those from a high-temperature combustion method (elemental analysis); the latter providing the true carbon content value. By comparing the above data, it was found that the UVP-TOC method significantly underestimated the TOC value of the surface water samples, and it was also found that the model components BSA (protein) and HA (humic substances, HS) had a substantial influence on the TOC underestimation, while the SA (polysaccharide), TA (complex organic molecule) and CA/BA (small molecules) had little effect. The results showed that the agglomeration within and between BSA and HA molecules was an important reason for the inaccurate UVP-TOC values of BSA and HA. A further limitation was that for BSA, surfactants (e.g. sodium dodecylbenzene sulfonate, SDBS) and other surfactant-like substances, foam was formed during the CO₂ removal purging process by N₂ that seriously interfered with the determination of TOC. The study provides new information and insight into the causes of inaccuracies in the UVP-TOC analysis of surface waters and possible approaches to improve the accuracy.

Lessons from combined experimental and theoretical examination of the FTIR and 2D-IR spectroelectrochemistry of the amide I region of cytochrome c

Amide I difference spectroscopy is widely used to investigate protein function and structure changes. In this article, we show that the common approach of assigning features in amide I difference signals to distinct secondary structure elements in many cases may not be justified. Evidence comes from Fourier transform infrared (FTIR) and 2D-IR spectroelectrochemistry of the protein cytochrome c in the amide I range, in combination with computational spectroscopy based on molecular dynamics (MD) simulations. This combination reveals that each secondary structure unit, such as an alpha-helix or a beta-sheet, exhibits broad overlapping contributions, usually spanning a large part of the amide I region, which in the case of difference absorption experiments (such as in FTIR spectroelectrochemistry) may lead to intensity-compensating and even sign-changing contributions. We use cytochrome c as the test case, as this small electron-transferring redox-active protein contains different kinds of secondary structure units. Upon switching its redox-state, the protein exhibits a different charge distribution while largely retaining its structural scaffold. Our theoretical analysis suggests that the change in charge distribution contributes to the spectral changes and that structural changes are small. However, in order to confidently interpret FTIR amide I difference signals in cytochrome c and proteins in general, MD simulations in combination with additional experimental approaches such as isotope labeling, the insertion of infrared labels to selectively probe local structural elements will be required. In case these data are not available, a critical assessment of previous interpretations of protein amide I 1D- and 2D-IR difference spectroscopy data is warranted.

Does liquid-liquid phase separation drive peptide folding?

Proline–arginine (PR) dipeptide repeats have been shown to undergo liquid–liquid phase separation and are an example of a growing number of intrinsically disordered proteins that can assemble into membraneless organelles. These structures have been posited as nucleation sites for pathogenic protein aggregation. As such, a better understanding of the effects that the increased local concentration and volumetric crowding within droplets have on peptide secondary structure is necessary. Herein we use Fourier transform infrared (FTIR) and two-dimensional infrared (2DIR) spectroscopy to show that formation of droplets by PR₂₀ accompanies changes in the amide-I spectra consistent with folding into poly-proline helical structures.

Two-dimensional infrared spectroscopy reveals folding of an intrinsically disordered peptide when sequestered into a model “membrane-less” organelle.

Structural Analysis of Nonapeptides Derived from Elastin

Elastin-derived peptides are released from the extracellular matrix remodeling by numerous proteases and seem to regulate many biological processes, notably cancer progression. The canonical elastin peptide is VGVAPG, which harbors the XGXXPG consensus pattern, allowing interaction with the elastin receptor complex located at the surface of cells. Besides these elastokines, another class of peptides has been identified. This group of bioactive elastin peptides presents the XGXPGXGXG consensus sequence, but the reason for their bioactivity remains unexplained. To better understand their nature and structure-function relationships, herein we searched the current databases for this nonapeptide motif and observed that the XGXPGXGXG elastin peptides define a specific group of tandemly repeated patterns. Further, we focused on four tandemly repeated human elastin nonapeptides, i.e., AGIPGLGVG, VGVPGLGVG, AGVPGLGVG, and AGVPGFGAG. These peptides were analyzed by means of optical spectroscopies and molecular dynamics. Ultraviolet-circular dichroism and Raman spectra are consistent with a mixture of β-turn, β-strand, and random-chain secondary elements in aqueous media. Quantitative analysis of their conformations suggested that turns corresponded to half of the total population of structural elements, whereas the remaining half were equally distributed between β-strand and unordered chains. These distributions were confirmed by molecular dynamics simulations. Altogether, our data suggest that these highly dynamic peptides harbor a type II β-turn located in their central part. We hypothesize that this structural element could explain their specific bioactivity.

Selenium-chelating corn oligopeptide as a potential antioxidant supplement: investigation of the protein conformational changes and identification of the antioxidant fragment composition

A selenium-chelating corn oligopeptide (Se-COP) with high protein and low molecular weight was prepared as a selenium supplement. We utilized infrared (IR), ultraviolet (UV), and circular dichroism (CD) spectroscopy, 1-anilinonaphthalene-8-sulfonate (ANS)-binding fluorescence spectra, and isothermal titration calorimetry (ITC) to analyze and describe Se-COP and its reactions. It was concluded that the chelation reaction was a spontaneous process driven by enthalpy and entropy, with ΔH=3.79 × 104 ± 4075 cal/mol, ΔS = 146 cal/mol, ΔG = –23356.30 ± 126.94 cal/mol, binding constant Ka = 1.18 × 104 ± 855 M–1, and binding site number n = 0.13 ± 0.0126, and described as coordination bonds forming and hydrophobic interaction, as well as protein conformational changes including secondary and tertiary hydrophobic structure. Se-COP had strong antioxidant capacity, and mass spectrometry (MS) was used to identify the antioxidant peptide fragment, which was characterized as LLPPY and quantified at 428.95 ng/mg. This study indicated that Se-COP prepared by chelation may be a Se supplement with antioxidant capacity that can be applied in functional foods or ingredients.

Temperature-Induced Collapse of Elastin-like Peptides Studied by 2DIR Spectroscopy

Elastin-like peptides are hydrophobic biopolymers that exhibit a reversible coacervation transition when the temperature is raised above a critical point. Here, we use a combination of linear infrared spectroscopy, two-dimensional infrared spectroscopy, and molecular dynamics simulations to study the structural dynamics of two elastin-like peptides. Specifically, we investigate the effect of the solvent environment and temperature on the structural dynamics of a short (5-residue) elastin-like peptide and of a long (450-residue) elastin-like peptide. We identify two vibrational energy transfer processes that take place within the amide I' band of both peptides. We observe that the rate constant of one of the exchange processes is strongly dependent on the solvent environment and argue that the coacervation transition is accompanied by a desolvation of the peptide backbone where up to 75% of the water molecules are displaced. We also study the spectral diffusion dynamics of the valine(1) residue that is present in both peptides. We find that these dynamics are relatively slow and indicative of an amide group that is shielded from the solvent. We conclude that the coacervation transition of elastin-like peptides is probably not associated with a conformational change involving this residue.

Elastin-like Peptide in Confinement: FT-IR and NMR T 1 Relaxation Data

We employed FT-IR and NMR experiments to investigate the influence of a cell-mimicking crowding environment on the structure and dynamics of an elastin-like peptide (ELP) with the sequence GVG(VPGVG)3, which – due to a high number of hydrophobic amino acid side chains – exhibits an inverse temperature transition (ITT). As simplified crowding agent, we used 30 wt% Ficoll. The FT-IR data revealed the well-known broad ITT above ~25°C, as observed by the decrease of the relative population of random coil structures and the concomitant increase of type II β-turns. Interestingly, the addition of Ficoll leads to a destabilizing effect of type II β-turn structures. This is in contrast to the expected excluded-volume effect of the macromolecular crowder, but can be explained by weak interactions of the peptide with the polysaccharide chains of the crowding agent. Further, the crowding agent leads to the onset of a reversal of the folding transition at high temperatures. The full assignment of the ELP allowed for a residue-specific investigation of the dynamic behavior of ELP by NMR. Due to a strong change of microscopic viscosity between native/buffered conditions and crowded conditions, relaxation data remain inconclusive with respect to the observation of an ITT. Hence, no quantitative details in terms of internal conformational changes can be obtained. However, temperature dependent differences in the 13C relaxation behavior between core and terminal parts of the peptide indicate temperature induced changes in the internal dynamics with generally higher internal mobility at chain ends: This is in full agreement with FT-IR data. In harmony with the FT-IR analysis, macromolecular crowding does not lead to significant changes in the relaxation behavior.

Hydrophobic Collapse in N-Methylacetamide-Water Mixtures

Aqueous N-methylacetamide solutions were investigated by polarization-resolved pump-probe and 2D infrared spectroscopy (2D IR), using the amide I mode as a reporter. The 2D IR results are compared with molecular dynamics simulations and spectral calculations to gain insight into the molecular structures in the mixture. N-Methylacetamide and water molecules tend to form clusters with "frozen" amide I dynamics. This is driven by a hydrophobic collapse as the methyl groups of the N-methylacetamide molecules cluster in the presence of water. Since the studied system can be considered as a simplified model for the backbone of proteins, the present study forms a convenient basis for understanding the structural and vibrational dynamics in proteins. It is particularly interesting to find out that a hydrophobic collapse as the one driving protein folding is observed in such a simple system.

Revealing the Solvation Structure and Dynamics of Carbonate Electrolytes in Lithium-Ion Batteries by Two-Dimensional Infrared Spectrum Modeling

Carbonate electrolytes in lithium-ion batteries play a crucial role in conducting lithium ions between two electrodes. Mixed solvent electrolytes consisting of linear and cyclic carbonates are commonly used in commercial lithium-ion batteries. To understand how the linear and cyclic carbonates introduce different solvation structures and dynamics, we performed molecular dynamics simulations of two representative electrolyte systems containing either linear or cyclic carbonate solvents. We then modeled their two-dimensional infrared (2DIR) spectra of the carbonyl stretching mode of these carbonate molecules. We found that the chemical exchange process involving formation and dissociation of lithium-ion/carbonate complexes is responsible for the growth of 2DIR cross peaks with increasing waiting time. In addition, we also found that cyclic carbonates introduce faster dynamics of dissociation and formation of lithium-ion/carbonate complexes than linear carbonates. These findings provide new insights into understanding the lithium-ion mobility and its interplay with solvation structure and ultrafast dynamics in carbonate electrolytes used in lithium-ion batteries.

Interplay between Hydrogen Bonding and Vibrational Coupling in Liquid N-Methylacetamide

Intrinsically disordered proteins play an important role in biology, and unraveling their labile structure presents a vital challenge. However, the dynamical structure of such proteins thwarts their study by standard techniques such as X-ray diffraction and NMR spectroscopy. Here, we use a neat liquid composed of N-methylacetamide molecules as a model system to elucidate dynamical and structural properties similar to those one can expect to see in intrinsically disordered proteins. To examine the structural dynamics in the neat liquid, we combine molecular dynamics, response-function-based spectral simulations, and two-dimensional polarization-resolved infrared spectroscopy in the amide I (CO stretch) region. The two-dimensional spectra reveal a delicate interplay between hydrogen bonding and intermolecular vibrational coupling effects, observed through a fast anisotropy decay. The present study constitutes a general platform for understanding the structure and dynamics of highly disordered proteins.

Refining Disordered Peptide Ensembles with Computational Amide I Spectroscopy: Application to Elastin-Like Peptides

The characterization of intrinsically disordered protein (IDP) ensembles is complicated both by inherent heterogeneity and by the fact that many common experimental techniques function poorly when applied to IDPs. For this reason, the development of alternative structural tools for probing IDP ensembles has attracted considerable attention. Here we describe our recent work in developing experimental and computational tools for characterizing IDP ensembles using Amide I (backbone carbonyl stretch) vibrational spectroscopy. In this approach, the infrared (IR) absorption frequencies of isotope-labeled amide bonds probe their local electrostatic environments and structures. Empirical frequency maps allow us to use this spectroscopic data as a direct experimental test of atomistic structural models. We apply these methods to a family of short elastin-like peptides (ELPs), fragments of the elastin protein based around the Pro-Gly turn motif characteristic of the elastomeric segments of the full protein. Using a maximum entropy analysis that applies constraints from experimental spectra to weighting predicted spectra from molecular dynamics (MD) ensembles, we find that peptides with Ala or Val side chains preceding the Pro-Gly turn unit exhibit a stronger tendency toward extended structures than do Gly-Pro-Gly motifs, suggesting an important role for steric interactions in tuning the molecular properties of elastin.

Assessing Spectral Simulation Protocols for the Amide I Band of Proteins

We present a benchmark study of spectral simulation protocols for the amide I band of proteins. The amide I band is widely used in infrared spectroscopy of proteins due to the large signal intensity, high sensitivity to hydrogen bonding, and secondary structural motifs. This band has, thus, proven valuable in many studies of protein structure-function relationships. We benchmark spectral simulation protocols using two common force fields in combination with several electrostatic mappings and coupling models. The results are validated against experimental linear absorption and two-dimensional infrared spectroscopy for three well-studied proteins. We find two-dimensional infrared spectroscopy to be much more sensitive to the simulation protocol than linear absorption and report on the best simulation protocols. The findings demonstrate that there is still room for ideas to improve the existing models for the amide I band of proteins.

Application of two-dimensional infrared spectroscopy to benchmark models for the amide I band of proteins

In this paper, we present a novel benchmarking method for validating the modelling of vibrational spectra for the amide I region of proteins. We use the linear absorption spectra and two-dimensional infrared spectra of four experimentally well-studied proteins as a reference and test nine combinations of molecular dynamics force fields, vibrational frequency mappings, and coupling models. We find that two-dimensional infrared spectra provide a much stronger test of the models than linear absorption does. The best modelling approach in the present study still leaves significant room for future improvement. The presented benchmarking scheme, thus, provides a way of validating future protocols for modelling the amide I band in proteins.

Structural disorder of folded proteins: isotope-edited 2D IR spectroscopy and Markov state modeling

The conformational heterogeneity of the N-terminal domain of the ribosomal protein L9 (NTL91-39) in its folded state is investigated using isotope-edited two-dimensional infrared spectroscopy. Backbone carbonyls are isotope-labeled ((13)C=(18)O) at five selected positions (V3, V9, V9G13, G16, and G24) to provide a set of localized spectroscopic probes of the structure and solvent exposure at these positions. Structural interpretation of the amide I line shapes is enabled by spectral simulations carried out on structures extracted from a recent Markov state model. The V3 label spectrum indicates that the β-sheet contacts between strands I and II are well folded with minimal disorder. The V9 and V9G13 label spectra, which directly probe the hydrogen-bond contacts across the β-turn, show significant disorder, indicating that molecular dynamics simulations tend to overstabilize ideally folded β-turn structures in NTL91-39. In addition, G24-label spectra provide evidence for a partially disordered α-helix backbone that participates in hydrogen bonding with the surrounding water.

Coacervation of α-elastin studied by ultrafast nonlinear infrared spectroscopy

Elastin is the main protein to confer elasticity to biological tissues, through the formation of a hierarchical network of fibres.

Characterizing Anharmonic Vibrational Modes of Quinones with Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared (2D IR) spectroscopy was used to study the vibrational modes of three quinones--benzoquinone, naphthoquinone, and anthraquinone. The vibrations of interest were in the spectral range of 1560-1710 cm(-1), corresponding to the in-plane carbonyl and ring stretching vibrations. Coupling between the vibrational modes is indicated by the cross peaks in the 2D IR spectra. The diagonal and off-diagonal anharmonicities range from 4.6 to 17.4 cm(-1) for the quinone series. In addition, there is significant vibrational coupling between the in-plane carbonyl and ring stretching vibrations. The diagonal anharmonicity, off-diagonal anharmonicity, and vibrational coupling constants are reported for benzoquinone, naphthoquinone, and anthraquinone.

Identification and relative quantification of tyrosine nitration in a model peptide using two-dimensional infrared spectroscopy

Nitration of tyrosine in proteins and peptides is a post-translational modification that occurs under conditions of oxidative stress. It is implicated in a variety of medical conditions, including neurodegenerative and cardiovascular diseases. However, monitoring tyrosine nitration and understanding its role in modifying biological function remains a major challenge. In this work, we investigate the use of electron-vibration-vibration (EVV) two-dimensional infrared (2DIR) spectroscopy for the study of tyrosine nitration in model peptides. We demonstrate the ability of EVV 2DIR spectroscopy to differentiate between the neutral and deprotonated states of 3-nitrotyrosine, and we characterize their spectral signatures using information obtained from quantum chemistry calculations and simulated EVV 2DIR spectra. To test the sensitivity of the technique, we use mixed-peptide samples containing various levels of tyrosine nitration, and we use mass spectrometry to independently verify the level of nitration. We conclude that EVV 2DIR spectroscopy is able to provide detailed spectroscopic information on peptide side-chain modifications and to detect nitration levels down to 1%. We further propose that lower nitration levels could be detected by introducing a resonant Raman probe step to increase the detection sensitivity of EVV 2DIR spectroscopy.

Observation of vibrational energy exchange in a type-III antifreeze protein

We performed time- and polarization-resolved pump-probe and two-dimensional infrared (2D-IR) experiments to study the dynamics of the amide I vibration of a 7 kDa type-III antifreeze protein. In the pump-probe experiments, we used femtosecond mid-infrared pulses to investigate the vibrational relaxation dynamics of the amide mode. The transient spectra show the presence of two spectral components that decay with different lifetimes, indicative of the presence of two distinct amide subbands. The 2D-IR experiments reveal the coupling between the two bands in the form of cross-peaks. On the basis of previous work by Demirdöven et al. ( J. Am. Chem. Soc. 2004 , 126 , 7981 - 7990 ), we assign the observed bands to the two infrared-active modes α(-) and α(+) found in protein β-sheets. The amplitudes of the cross-peak were found to increase with delay time, indicating that the cross-peaks originate from population transfer between the coupled amide oscillators. The time constant of the energy transfer was found to be 6-7 ps.

Linear absorption and two-dimensional infrared spectra of N-methylacetamide in chloroform revisited: polarizability and multipole effects

The effect of solvent polarizability and multipole effects on the amide I vibrational spectra of a peptide unit is investigated. Four molecular dynamics force fields of increasing complexity for the solvent are used to model both the linear absorption and two-dimensional infrared spectra. It is observed that, at least in chloroform solution, the predicted solvent shift is considerably improved when accounting for the polarizabiltiy and multipole effects. The latter are typically connected with halogen bonding. Significant deviations are still observed for more sensitive line shape parameters such as the spectral width and line skewness. However, the findings demonstrate that previously observed deviations have an origin in the force field treatment rather than in the electrostatic mapping procedure frequently employed to simulate linear absorption and two-dimensional infrared spectroscopy.

Electrostatic frequency shifts in amide I vibrational spectra: direct parameterization against experiment

The interpretation of protein amide I infrared spectra has been greatly assisted by the observation that the vibrational frequency of a peptide unit reports on its local electrostatic environment. However, the interpretation of spectra remains largely qualitative due to a lack of direct quantitative connections between computational models and experimental data. Here, we present an empirical parameterization of an electrostatic amide I frequency map derived from the infrared absorption spectra of 28 dipeptides. The observed frequency shifts are analyzed in terms of the local electrostatic potential, field, and field gradient, evaluated at sites near the amide bond in molecular dynamics simulations. We find that the frequency shifts observed in experiment correlate very well with the electric field in the direction of the C=O bond evaluated at the position of the amide oxygen atom. A linear best-fit mapping between observed frequencies and electric field yield sample standard deviations of 2.8 and 3.7 cm(-1) for the CHARMM27 and OPLS-AA force fields, respectively, and maximum deviations (within our data set) of 9 cm(-1). These results are discussed in the broader context of amide I vibrational models and the effort to produce quantitative agreement between simulated and experimental absorption spectra.

Folding dynamics of the Trp-cage miniprotein: evidence for a native-like intermediate from combined time-resolved vibrational spectroscopy and molecular dynamics simulations

Trp-cage is a synthetic 20-residue miniprotein which folds rapidly and spontaneously to a well-defined globular structure more typical of larger proteins. Due to its small size and fast folding, it is an ideal model system for experimental and theoretical investigations of protein folding mechanisms. However, Trp-cage's exact folding mechanism is still a matter of debate. Here we investigate Trp-cage's relaxation dynamics in the amide I' spectral region (1530-1700 cm(-1)) using time-resolved infrared spectroscopy. Residue-specific information was obtained by incorporating an isotopic label ((13)C═(18)O) into the amide carbonyl group of residue Gly11, thereby spectrally isolating an individual 310-helical residue. The folding-unfolding equilibrium is perturbed using a nanosecond temperature-jump (T-jump), and the subsequent re-equilibration is probed by observing the time-dependent vibrational response in the amide I' region. We observe bimodal relaxation kinetics with time constants of 100 ± 10 and 770 ± 40 ns at 322 K, suggesting that the folding involves an intermediate state, the character of which can be determined from the time- and frequency-resolved data. We find that the relaxation dynamics close to the melting temperature involve fast fluctuations in the polyproline II region, whereas the slower process can be attributed to conformational rearrangements due to the global (un)folding transition of the protein. Combined analysis of our T-jump data and molecular dynamics simulations indicates that the formation of a well-defined α-helix precedes the rapid formation of the hydrophobic cage structure, implying a native-like folding intermediate, that mainly differs from the folded conformation in the orientation of the C-terminal polyproline II helix relative to the N-terminal part of the backbone. We find that the main free-energy barrier is positioned between the folding intermediate and the unfolded state ensemble, and that it involves the formation of the α-helix, the 310-helix, and the Asp9-Arg16 salt bridge. Our results suggest that at low temperature (T ≪ Tm) a folding path via formation of α-helical contacts followed by hydrophobic clustering becomes more important.

Modeling the thermal unfolding 2DIR spectra of a β-hairpin peptide based on the implicit solvent MD simulation

We simulated the equilibrium isotope-edited FTIR and 2DIR spectra of a β-hairpin peptide trpzip2 at a series of temperatures. The simulation was based on the configuration distributions generated using the GB(OBC) implicit solvent model and the integrated tempering sampling (ITS) technique. A soaking procedure was adapted to generate the peptide in explicit solvent configurations for the spectroscopy calculations. The nonlinear exciton propagation (NEP) method was then used to calculate the spectra. Agreeing with the experiments, the intensities and ellipticities of the isotope-shifted peaks in our simulated signals have the site-specific temperature dependences, which suggest the inhomogeneous local thermal stabilities along the peptide chain. Our simulation thus proposes a cost-effective means to understand a peptide's conformational change and related IR spectra across its thermal unfolding transition.

Thermally induced protein unfolding probed by isotope-edited IR spectroscopy

Infrared (IR) spectroscopy has been widely utilized for the study of protein folding, unfolding, and misfolding processes. We have previously developed a theoretical method for calculating IR spectra of proteins in the amide I region. In this work, we apply this method, in combination with replica-exchange molecular dynamics simulations, to study the equilibrium thermal unfolding transition of the villin headpiece subdomain (HP36). Temperature-dependent IR spectra and spectral densities are calculated. The spectral densities correctly reflect the unfolding conformational changes in the simulation. With the help of isotope labeling, we are able to capture the feature that helix 2 of HP36 loses its secondary structure before global unfolding occurs, in agreement with experiment.