Multivariate Analysis of Spectral Data with Frequency Shifts: Application to Temperature Dependent Infrared Spectra of Peptides and Proteins

Multi-Probe Equilibrium Analysis of Gradual (Un)Folding Processes

Studies of small proteins that exhibit noncooperative, gradual (un)folding can offer unique insights into the rarely accessible intermediate stages of the protein folding processes. Detailed experimental characterization of these intermediate states requires approaches that utilize multiple site-specific probes of the local structure. Isotopically edited infrared (IR) spectroscopy has emerged as a powerful methodology capable of providing such high-resolution structural information. Labeling of selected amide carbonyls with ¹³C results in detectable side-bands of amide I' vibrations, which are sensitive to local conformation and/or solvent exposure without introducing any significant structural perturbation to the protein. Incorporation of isotopically labeled amino acids at specific positions can be achieved by the chemical synthesis of the studied proteins. We describe the basic procedures for synthesis of ¹³C isotopically edited protein samples, experimental IR spectroscopic measurements and analysis of the site-specific equilibrium thermal unfolding of a small protein from the temperature-dependent IR data.

Compensation of temperature effects on spectra through evolutionary rank analysis

Spectra measured in various ranges of temperature are usually slightly different from each other in shape and position of the bands. Although the displayed inconsistencies are rather small, yet may lead to incorrect analysis and interpretation of the collected spectrothermal data. Thus the unspecific spectral effects induced by temperature, in particular the thermal shifts and broadening of the bands, have to be compensated. In the paper, a simple two-step method of thermospectral dataset uniformisation is presented. Thermally induced 'movement' of the bands is approximated as a linear function of the difference of temperatures, so the co-shifting of the spectra is done linearly. Thermal broadening is mimicked by convoluting the low-temperature signal (spectrum) with a Gaussian or Lorentzian spreading filter. Proper widths (values of FWHM) of these filters, used to uniform the whole dataset, are assumed to depend on the difference of temperatures, in a form of one-parameter functions. This assumption, which has been empirically confirmed, is a fundamental premise of the method of Partial Compensation for Thermal Broadening (PCTB). Optimal values of the parameters of all the functions, used to compensate both thermal shifting and broadening, are found by the Evolutionary Rank Analysis (ERA) applied on an evolving data matrix. Efficiency of the proposed approach was verified on the UV-Vis thermospectral dataset of one-component model systems. In addition, since the method is aimed at making uniformed the thermospectral datasets of multi-component systems with similar spectral properties of individual components, the two-component conformer system of t-APE (trans-1-(2'-anthryl)-2-phenylethene) has also been analyzed.

Effect of Mutations on the Global and Site-Specific Stability and Folding of an Elementary Protein Structural Motif

Understanding the folding mechanism of proteins requires detailed knowledge of the roles of individual amino acid residues in stabilization of specific elements and local segments of the native structure. Recently, we have utilized the combination of circular dichroism (CD) and site-specific ¹³C isotopically edited infrared spectroscopy (IR) coupled with the Ising-like model for protein folding to map the thermal unfolding at the residue level of a de novo designed helix-turn-helix motif αtα. Here we use the same methodology to study how the sequence of local thermal unfolding is affected by selected mutations introduced into the most and least stable parts of the motif. Seven different mutants of αtα are screened to find substitutions with the most pronounced effects on the overall stability. Subsequently, thermal unfolding of two mutated αtα sequences is studied with site-specific resolution, using four distinct ¹³C isotopologues of each. The data are analyzed with the Ising-like model, which builds on a previous parametrization for the original αtα sequence and tests different ways of incorporating the amino acid substitution. We show that for both more and less stable mutants only the adjustment of all interaction parameters of the model can yield a satisfactory fit to the experimental data. The stabilizing and destabilizing mutations result, respectively, in a similar increase and decrease of the stability of all probed local segments, irrespective of their position with respect to the mutation site. Consequently, the relative order of their unfolding remains essentially unchanged. These results underline the importance of the interconnectivity of the stabilizing interaction network and cooperativity of the protein structure, which is evident even in a small motif with apparently noncooperative, heterogeneous unfolding. Overall, our findings are consistent with the native structure being the dominant factor in determining the folding mechanism, regardless of the details of its overall or local thermodynamic stabilization.

Interplay between Lipid Interaction and Homo-coiling of Membrane-Tethered Coiled-Coil Peptides

The designed coiled-coil-forming peptides E [(EIAALEK)3] and K [(KIAALKE)3] are known to trigger efficient membrane fusion when they are tethered to lipid vesicles in the form of lipopeptides. Knowledge of their secondary structure is a key element in understanding their role in membrane fusion. Special conditions can be found at the interface of the membrane, where the peptides are confined in close proximity to other peptide molecules as well as to the lipid interface. Consequently, different structural states were proposed for the peptides when tethered to this interface. Due to the multitude of possible states, determining the structure solely on the basis of circular dichroism (CD) spectra at a single temperature can be misleading. In addition, it has not yet been possible to unambiguously distinguish between the membrane-bound and the coiled-coil states of these peptides by means of infrared (IR) spectroscopy due to their very similar amide I' bands. Here, the molecular basis of this similarity is investigated by means of site-specific (13)C-labeled FTIR spectroscopy. Structural similarities between the membrane-interacting helix of K and the homo-coiled-coil-forming helix of E are shown to cause the similar spectroscopic properties. Furthermore, the peptide structure is investigated using temperature-dependent CD and IR spectroscopy, and it is shown that the different states can be distinguished on the basis of their thermal behavior. It is shown that the two peptides behave fundamentaly differently when tethered to the lipid membrane, which implies that their role during membrane fusion is different and the mechanism of this process is asymmetric.

Sequence, structure, and cooperativity in folding of elementary protein structural motifs

Residue-level unfolding of two helix-turn-helix proteins--one naturally occurring and one de novo designed--is reconstructed from multiple sets of site-specific (13)C isotopically edited infrared (IR) and circular dichroism (CD) data using Ising-like statistical-mechanical models. Several model variants are parameterized to test the importance of sequence-specific interactions (approximated by Miyazawa-Jernigan statistical potentials), local structural flexibility (derived from the ensemble of NMR structures), interhelical hydrogen bonds, and native contacts separated by intervening disordered regions (through the Wako-Saitô-Muñoz-Eaton scheme, which disallows such configurations). The models are optimized by directly simulating experimental observables: CD ellipticity at 222 nm for model proteins and their fragments and (13)C-amide I' bands for multiple isotopologues of each protein. We find that data can be quantitatively reproduced by the model that allows two interacting segments flanking a disordered loop (double sequence approximation) and incorporates flexibility in the native contact maps, but neither sequence-specific interactions nor hydrogen bonds are required. The near-identical free energy profiles as a function of the global order parameter are consistent with expected similar folding kinetics for nearly identical structures. However, the predicted folding mechanism for the two motifs is different, reflecting the order of local stability. We introduce free energy profiles for "experimental" reaction coordinates--namely, the degree of local folding as sensed by site-specific (13)C-edited IR, which highlight folding heterogeneity and contrast its overall, average description with the detailed, local picture.

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

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