Spectroscopic studies of protein folding: linear and nonlinear methods

How quickly can a β-hairpin fold from its transition state?

Understanding the structural nature of the free energy bottleneck(s) encountered in protein folding is essential to elucidating the underlying dynamics and mechanism. For this reason, several techniques, including Φ-value analysis, have previously been developed to infer the structural characteristics of such high free-energy or transition states. Herein we propose that one (or few) appropriately placed backbone and/or side chain cross-linkers, such as disulfides, could be used to populate a thermodynamically accessible conformational state that mimics the folding transition state. Specifically, we test this hypothesis on a model β-hairpin, Trpzip4, as its folding mechanism has been extensively studied and is well understood. Our results show that cross-linking the two β-strands near the turn region increases the folding rate by an order of magnitude, to about (500 ns)⁻¹, whereas cross-linking the termini results in a hyperstable β-hairpin that has essentially the same folding rate as the uncross-linked peptide. Taken together, these findings suggest that cross-linking is not only a useful strategy to manipulate folding free energy barriers, as shown in other studies, but also, in some cases, it can be used to stabilize a folding transition state analogue and allow for direct assessment of the folding process on the downhill side of the free energy barrier. The calculated free energy landscape of the cross-linked Trpzip4 also supports this picture. An empirical analysis further suggests, when folding of β-hairpins does not involve a significant free energy barrier, the folding time (τ) follows a power law dependence on the number of hydrogen bonds to be formed (n_H), namely, τ = τ₀n_H^α, with τ₀ = 20 ns and α = 2.3.

Discriminating trpzip2 and trpzip4 peptides' folding landscape using the two-dimensional infrared spectroscopy: a simulation study

We analyzed, based on the theoretical spectroscopic modeling, how the differences in the folding landscapes of two β-hairpin peptides trpzip2 and trpzip4 are reflected in their thermal unfolding infrared measurements. The isotope-edited equilibrium FTIR and two dimensional infrared spectra of the two peptides were calculated, using the nonlinear exciton propagation method, at a series of temperatures. The spectra calculations were based on the configuration distributions generated using the GB(OBC) implicit solvent MD simulation and the integrated tempering sampling technique. Conformational analysis revealed the different local thermal stabilities for these two peptides, which suggested the different folding landscapes. Our study further suggested that the ellipticities of the isotope peaks in the coherent IR signals are more sensitive to these local stability differences compared with other spectral features such as the peak intensities. Our technique can thus be combined with the relevant experimental measurements to achieve a better understanding of the peptide folding behaviors.

Ultrafast hydrogen exchange reveals specific structural events during the initial stages of folding of cytochrome c

Many proteins undergo a sharp decrease in chain dimensions during early stages of folding, prior to the rate-limiting step in folding. However, it remains unclear whether compact states are the result of specific folding events or a general hydrophobic collapse of the poly peptide chain driven by the change in solvent conditions. To address this fundamental question, we extended the temporal resolution of NMR-detected H/D exchange labeling experiments into the microsecond regime by adopting a microfluidics approach. By observing the competition between H/D exchange and folding as a function of labeling pH, coupled with direct measurement of exchange rates in the unfolded state, we were able to monitor hydrogen-bond formation for over 50 individual backbone NH groups within the initial 140 microseconds of folding of horse cytochrome c. Clusters of solvent-shielded amide protons were observed in two α-helical segments in the C-terminal half of the protein, while the N-terminal helix remained largely unstructured, suggesting that proximity in the primary structure is a major factor in promoting helix formation and association at early stages of folding, while the entropically more costly long-range contacts between the N- and C-terminal helices are established only during later stages. Our findings clearly indicate that the initial chain condensation in cytochrome c is driven by specific interactions among a subset of α-helical segments rather than a general hydrophobic collapse.

Using D-Amino Acids to Delineate the Mechanism of Protein Folding: Application to Trp-cage

Using the miniprotein Trp-cage as a model, we show that D-amino acids can be used to facilitate the delineation of protein folding mechanism. Specifically, we study the folding-unfolding kinetics of three Trp-cage mutants where the native glycine residue near the C-terminus of the α-helix is replaced by a D-amino acid. A previous study showed that these mutations increase the Trp-cage stability, due to a terminal capping effect. Our results show that the stabilizing effect of D-asparagine and D-glutamine originates almost exclusively from a decrease in the unfolding rate, while the D-alanine mutation results in a similar decrease in the unfolding rate, but it also increases the folding rate. Together, these results support a folding mechanism wherein the α-helix formation in the transition state is nucleated at the N-terminus, whereas those long-range native interactions stabilizing this helix are developed at the downhill side of the folding free energy barrier.

Using VIPT-jump to distinguish between different folding mechanisms: application to BBL and a Trpzip

Protein folding involves a large number of sequential molecular steps or conformational substates. Thus, experimental characterization of the underlying folding energy landscape for any given protein is difficult. Herein, we present a new method that can be used to determine the major characteristics of the folding energy landscape in question, e.g., to distinguish between activated and barrierless downhill folding scenarios. This method is based on the idea that the conformational relaxation kinetics of different folding mechanisms at a given final condition will show different dependences on the initial condition. We show, using both simulation and experiment, that it is possible to differentiate between disparate kinetic folding models by comparing temperature jump (T-jump) relaxation traces obtained with a fixed final temperature and varied initial temperatures, which effectively varies the initial potential (VIP) of the system of interest. We apply this method (hereafter refer to as VIPT-jump) to two model systems, tryptophan zipper (Trpzip)-2c and BBL, and our results show that BBL exhibits characteristics of barrierless downhill folding, whereas Trpzip-2c folding encounters a free energy barrier. In addition, using the T-jump data of BBL we are able to provide, via Langevin dynamics simulations, a realistic estimate of its conformational diffusion coefficient.

α-helix to β-hairpin transition of human amylin monomer

The human islet amylin polypeptide is produced along with insulin by pancreatic islets. Under some circumstances, amylin can aggregate to form amyloid fibrils, whose presence in pancreatic cells is a common pathological feature of Type II diabetes. A growing body of evidence indicates that small, early stage aggregates of amylin are cytotoxic. A better understanding of the early stages of the amylin aggregation process and, in particular, of the nucleation events leading to fibril growth could help identify therapeutic strategies. Recent studies have shown that, in dilute solution, human amylin can adopt an α-helical conformation, a β-hairpin conformation, or an unstructured coil conformation. While such states have comparable free energies, the β-hairpin state exhibits a large propensity towards aggregation. In this work, we present a detailed computational analysis of the folding pathways that arise between the various conformational states of human amylin in water. A free energy surface for amylin in explicit water is first constructed by resorting to advanced sampling techniques. Extensive transition path sampling simulations are then employed to identify the preferred folding mechanisms between distinct minima on that surface. Our results reveal that the α-helical conformer of amylin undergoes a transformation into the β-hairpin monomer through one of two mechanisms. In the first, misfolding begins through formation of specific contacts near the turn region, and proceeds via a zipping mechanism. In the second, misfolding occurs through an unstructured coil intermediate. The transition states for these processes are identified. Taken together, the findings presented in this work suggest that the inter-conversion of amylin between an α-helix and a β-hairpin is an activated process and could constitute the nucleation event for fibril growth.

Monitoring guanidinium-induced structural changes in ribonuclease proteins using Raman spectroscopy and 2D correlation analysis

Assessing the stability of proteins by comparing their unfolding profiles is a very important characterization and quality control step for any biopharmaceutical, and this is usually measured by fluorescence spectroscopy. In this paper we propose Raman spectroscopy as a rapid, noninvasive alternative analytical method and we shall show this has enhanced sensitivity and can therefore reveal very subtle protein conformational changes that are not observed with fluorescence measurements. Raman spectroscopy is a powerful nondestructive method that has a strong history of applications in protein characterization. In this work we describe how Raman microscopy can be used as a fast and reliable method of tracking protein unfolding in the presence of a chemical denaturant. We have compared Raman spectroscopic data to the equivalent samples analyzed using fluorescence spectroscopy in order to validate the Raman approach. Calculations from both Raman and fluorescence unfolding curves of [D]50 values and Gibbs free energy correlate well with each other and more importantly agree with the values found in the literature for these proteins. In addition, 2D correlation analysis has been performed on both Raman and fluorescence data sets in order to allow further comparisons of the unfolding behavior indicated by each method. As many biopharmaceuticals are glycosylated in order to be functional, we compare the unfolding profiles of a protein (RNase A) and a glycoprotein (RNase B) as measured by Raman spectroscopy and discuss the implications that glycosylation has on the stability of the protein.

Quenching of p-Cyanophenylalanine Fluorescence by Various Anions

To expand the spectroscopic utility of the non-natural amino acid p-cyanophenylalanine (Phe_CN), we examine the quenching efficiencies of a series of commonly encountered anions toward its fluorescence. We find that iodide exhibits an unusually large Stern-Volmer quenching constant, making it a convenient choice in Phe_CN fluorescence quenching studies. Indeed, using the villin headpiece subdomain as a testbed we demonstrate that iodide quenching of Phe_CN fluorescence offers a convenient means to reveal protein conformational heterogeneity. Furthermore, we show that the amino group of Phe_CN strongly quenches its fluorescence, suggesting that Phe_CN could be used as a local pH sensor.

A fresh look at the male-specific region of the human Y chromosome

The Chromosome-centric Human Proteome Project (C-HPP) aims to systematically map the entire human proteome with the intent to enhance our understanding of human biology at the cellular level. This project attempts simultaneously to establish a sound basis for the development of diagnostic, prognostic, therapeutic, and preventive medical applications. In Iran, current efforts focus on mapping the proteome of the human Y chromosome. The male-specific region of the Y chromosome (MSY) is unique in many aspects and comprises 95% of the chromosome's length. The MSY continually retains its haploid state and is full of repeated sequences. It is responsible for important biological roles such as sex determination and male fertility. Here, we present the most recent update of MSY protein-encoding genes and their association with various traits and diseases including sex determination and reversal, spermatogenesis and male infertility, cancers such as prostate cancers, sex-specific effects on the brain and behavior, and graft-versus-host disease. We also present information available from RNA sequencing, protein-protein interaction, post-translational modification of MSY protein-coding genes and their implications in biological systems. An overview of Human Y chromosome Proteome Project is presented and a systematic approach is suggested to ensure that at least one of each predicted protein-coding gene's major representative proteins will be characterized in the context of its major anatomical sites of expression, its abundance, and its functional relevance in a biological and/or medical context. There are many technical and biological issues that will need to be overcome in order to accomplish the full scale mapping.

Solute's perspective on how trimethylamine oxide, urea, and guanidine hydrochloride affect water's hydrogen bonding ability

While the thermodynamic effects of trimethylamine oxide (TMAO), urea, and guanidine hydrochloride (GdnHCl) on protein stability are well understood, the underlying mechanisms of action are less well characterized and, in some cases, even under debate. Herein, we employ the stretching vibration of two infrared (IR) reporters, i.e., nitrile (C≡N) and carbonyl (C═O), to directly probe how these cosolvents mediate the ability of water to form hydrogen bonds with the solute of interest, e.g., a peptide. Our results show that these three agents, despite having different effects on protein stability, all act to decrease the strength of the hydrogen bonds formed between water and the infrared probe. While the behavior of TMAO appears to be consistent with its protein-protecting ability, those of urea and GdnHCl are inconsistent with their role as protein denaturants. The latter is of particular interest as it provides strong evidence indicating that although urea and GdnHCl can perturb the hydrogen-bonding property of water their protein-denaturing ability does not arise from a simple indirect mechanism.

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.

Using thioamides to site-specifically interrogate the dynamics of hydrogen bond formation in β-sheet folding

Thioamides are sterically almost identical to their oxoamide counterparts, but they are weaker hydrogen bond acceptors. Therefore, thioamide amino acids are excellent candidates for perturbing the energetics of backbone-backbone H-bonds in proteins and hence should be useful in elucidating protein folding mechanisms in a site-specific manner. Herein, we validate this approach by applying it to probe the dynamic role of interstrand H-bond formation in the folding kinetics of a well-studied β-hairpin, tryptophan zipper. Our results show that reducing the strength of the peptide's backbone-backbone H-bonds, except the one directly next to the β-turn, does not change the folding rate, suggesting that most native interstrand H-bonds in β-hairpins are formed only after the folding transition state.