Identification of equilibrium and kinetic intermediates involved in folding of urea-denatured creatine kinase

The unfolding transition and kinetic refolding of dimeric creatine kinase after urea denaturation were monitored by intrinsic fluorescence and far ultraviolet circular dichroism. An equilibrium intermediate and a kinetic folding intermediate were identified and characterized. The fluorescence intensity of the equilibrium intermediate is close to that of the unfolded state, whereas its ellipticity at 222 nm is about 50% of the native state. The transition curves measured by these two methods are therefore non-coincident. The kinetic folding intermediate, formed during the burst phase of refolding under native-like conditions, possesses 75% of the native secondary structure, but is mostly lacking in native tertiary structure. In moderate concentrations of urea, only the initial, rapid change in fluorescence intensity or negative ellipticity is observed, and the final state values do not reach the equivalent unfolding values. The unfolding and refolding transition curves measured under identical conditions are non-coincident within the transition from intermediate to fully unfolded state. It is observed by SDS-PAGE that disulfide bond-linked dimeric or oligomeric intermediates are formed in moderate urea concentrations, especially in the refolding reaction. These rapidly formed, soluble intermediates represent an off-pathway event that leads to the hysteresis in the refolding transition curves.     

Microsecond folding of the cold shock protein measured by a pressure-jump technique

A pressure-jump apparatus was employed in investigating the kinetics of protein unfolding and refolding. In the reaction cell, the pressure can be increased or decreased by 100-160 bar within 50-100 microseconds and then held constant. Thus, unfolding and refolding reactions in the time range from 70 microseconds to 70 s can be followed with this technique. Measurements are possible in the transition regions of thermally or denaturant-induced folding in a wide range of temperatures and solvent conditions. We used this pressure-jump method to determine the temperature dependence of the rate constants of unfolding and refolding of the cold shock protein of Bacillus subtilis and of three variants thereof with Phe --> Ala substitutions in the central beta-sheet region. For all variants, the change in heat capacity occurred in refolding between the unfolded and activated states, suggesting that the overall native-like character of the activated state of folding was not changed by the deletion of individual Phe side chains. The Phe27Ala mutation affected the rate of unfolding only; the Phe15Ala and Phe17Ala mutations changed the kinetics of both unfolding and refolding. Although the activated state of folding of the cold shock protein is overall native-like, individual side chains are still in a non-native environment.     

Distinct unfolding and refolding pathways of ribonuclease a revealed by heating and cooling temperature jumps

Heating and cooling temperature jumps (T-jumps) were performed using a newly developed technique to trigger unfolding and refolding of wild-type ribonuclease A and a tryptophan-containing variant (Y115W). From the linear Arrhenius plots of the microscopic folding and unfolding rate constants, activation enthalpy (ΔH^#), and activation entropy (ΔS^#) were determined to characterize the kinetic transition states (TS) for the unfolding and refolding reactions. The single TS of the wild-type protein was split into three for the Y115W variant. Two of these transition states, TS1 and TS2, characterize a slow kinetic phase, and one, TS3, a fast phase. Heating T-jumps induced protein unfolding via TS2 and TS3; cooling T-jumps induced refolding via TS1 and TS3. The observed speed of the fast phase increased at lower temperature, due to a strongly negative ΔH^# of the folding-rate constant. The results are consistent with a path-dependent protein folding/unfolding mechanism. TS1 and TS2 are likely to reflect X-Pro¹¹⁴ isomerization in the folded and unfolded protein, respectively, and TS3 the local conformational change of the β-hairpin comprising Trp¹¹⁵. A very fast protein folding/unfolding phase appears to precede both processes. The path dependence of the observed kinetics is suggestive of a rugged energy protein folding funnel.     

Folding pathway of FKBP12 and characterisation of the transition state.

The folding pathway of human FKBP12, a 12 kDa FK506-binding protein (immunophilin), has been characterised. Unfolding and refolding rate constants have been determined over a wide range of denaturant concentrations and data are shown to fit to a two-state model of folding in which only the denatured and native states are significantly populated, even in the absence of denaturant. This simple model for folding, in which no intermediate states are significantly populated, is further supported from stopped-flow circular dichroism experiments in which no fast "burst" phases are observed. FKBP12, with 107 residues, is the largest protein to date which folds with simple two-state kinetics in water (kF=4 s(-1)at 25 degrees C). The topological crossing of two loops in FKBP12, a structural element suggested to cause kinetic traps during folding, seems to have little effect on the folding pathway.The transition state for folding has been characterised by a series of experiments on wild-type FKBP12. Information on the thermodynamic nature of, the solvent accessibility of, and secondary structure in, the transition state was obtained from experiments measuring the unfolding and refolding rate constants as a function of temperature, denaturant concentration and trifluoroethanol concentration. In addition, unfolding and refolding studies in the presence of ligand provided information on the structure of the ligand-binding pocket in the transition state. The data suggest a compact transition state relative to the unfolded state with some 70 % of the surface area buried. The ligand-binding site, which is formed mainly by two loops, is largely unstructured in the transition state. The trifluoroethanol experiments suggest that the alpha-helix may be formed in the transition state. These results are compared with results from protein engineering studies and molecular dynamics simulations (see the accompanying paper).     

Three-state protein folding: experimental determination of free-energy profile

When considering protein folding with a transient intermediate, a difficulty arises as to determination of the rates of separate transitions. Here we overcome this problem, using the kinetic studies of the unfolding/refolding reactions of the three-state protein apomyoglobin as a model. Amplitudes of the protein refolding kinetic burst phase corresponding to the transition from the unfolded (U) to intermediate (I) state, that occurs prior to the native state (N) formation, allow us to estimate relative populations of the rapidly converting states at various final urea concentrations. On the basis of these proportions, a complicated experimental chevron plot has been deconvolved into the urea-dependent rates of the I<-->N and U<-->N transitions to give the dependence of free energies of the main transition state and of all three (N, I, and U) stable states on urea concentration.     

A folding model of α-lactalbumin deduced from the three-state denaturation mechanism

It has been shown that α-lactalbumin undergoes a three-state denaturation, involving a helical intermediate state, on treatment with guanidine hydrochloride. The unfolding of the protein and the characteristics of the intermediate state are examined by means of circular dichroism, difference spectra and pH-jump measurements to investigate the temperature dependence and kinetic properties of the unfolding and refolding, the pH dependence of the transition between the intermediate and the fully unfolded states, and the effect of disulphide bond reduction on the stabilization of the intermediate.The results show that the long-range specific interactions such as specific electrostatic interactions and disulphide linkages are not important for stabilizing the intermediate, and that the transition between the intermediate and the fully unfolded states is extremely rapid (a relaxation time of less than one millisecond) and may correspond to the helix-coil transition of a polypeptide backbone. On the other hand, the activation parameters of the transition between the native and the intermediate states have suggested that the final stabilization by charge-pair interactions is preceded by hydrophobic interactions in the process of going from the intermediate to the native state.The mechanism of folding of the protein is discussed, and the folding process from the fully unfolded to the native state is apparently divided into at least three main steps: (1) the formation of incipient helical structures dictated by local interactions; (2) the packing of the helical segments accompanied with hydrophobic interactions; (3) the final stabilization by the electrostatic interactions. The relevance to the current theoretical results on protein folding is also discussed.     

Fast collapse but slow formation of secondary structure elements in the refolding transition of E. coli adenylate kinase

The various models proposed for protein folding transition differ in their order of appearance of the basic steps during this process. In this study, steady state and time-resolved dynamic non-radiative excitation energy transfer (FRET and trFRET) combined with site specific labeling experiments were applied in order to characterize the initial transient ensemble of Escherichia coli adenylate kinase (AK) molecules upon shifting conditions from those favoring denaturation to refolding and from folding to denaturing. Three sets of labeled AK mutants were prepared, which were designed to probe the equilibrium and transient distributions of intramolecular segmental end-to-end distances. A 176 residue section (residues 28-203), which spans most of the 214 residue molecule, and two short secondary structure chain segments including an alpha-helix (residues 169-188) and a predominantly beta-strand region (residues 188-203), were labeled. Upon fast change of conditions from denaturing to folding, the end-to-end distance of the 176 residue chain section showed an immediate collapse to a mean value of 26 A. Under the same conditions, the two short secondary structure elements did not respond to this shift within the first ten milliseconds, and retained the characteristics of a fully unfolded state. Within the first 10 ms after changes of the solvent from folding to denaturing, only minor changes were observed at the local environments of residues 203 and 169. The response of these same local environments to the shift of conditions from denaturing to folding occurred within the dead time of the mixing device. Thus, the response of the CORE domain of AK to fast transfer from folding to unfolding conditions is slow at all three conformational levels that were probed, and for at least a few milliseconds the ensemble of folded molecules is maintained under unfolding conditions. A different order of the changes was observed upon initiation of refolding. The AK molecules undergo fast collapse to an ensemble of compact structures where the local environment of surface probes seems to be native-like but the two labeled secondary structure elements remain unfolded.     

The mechanism of folding of globular proteins. Equilibria and kinetics of conformational transitions of penicillinase from Staphylococcus aureus involving a state of intermediate conformation.

1. The thermodynamically reversible unfolding and refolding of penicillinase between the native and fully unfolded states were followed by using guanidinium chloride as denaturant. 2. The equilibria, studied by optical rotation, u.v. absorption, viscosity and enzyme activity, show the presence of a state of intermediate conformation, termed state H, which is stable at 20 degrees C in 0.8 M-guanidinium chloride. 3. The physical properties of this state show that it is slightly expanded with an intrinsic viscosity of 8 ml-g-1, that the 13 tyrosine residues, which are distributed through the primary sequence, are maximally exposed to the solvent and that the helix content is the same as that of the native state. 4. The kinetics of the transition between the native state, state H and the fully unfolded state were followed by u.v. absorption and by optical rotation. They are interpreted as showing that state H lies on the folding pathway between the native and fully unfolded states. 5. The transition between the native state and state H exhibits monophasic unfolding kinetics and biphasic refolding kinetics. This indicates that there must be at least two intermediate states in this process, at least one of which lies on the folding pathway which may also involve cul-de-sac paths. 6. The results are discussed in terms of a mechanism involving rapid stabilization of nucleation regions in a moderately compact but internally solvated structure, with 'native format' [Anfinsen (1973) Science 181, 233-230] secondary structure stabilized by tertiary interaction. The final and rate-limiting step in refolding involves shuffling of these structural elements into the native state. 7. This model is discussed in relation to folding in vivo.     

Direct observation of microscopic reversibility in single-molecule protein folding

Both folded and unfolded conformations should be observed for a protein at its melting temperature (T(m)), where DeltaG between these states is zero. In an all-atom molecular dynamics simulation of chymotrypsin inhibitor 2 (CI2) at its experimental T(m), the protein rapidly loses its low-temperature native structure; it then unfolds before refolding to a stable, native-like conformation. The initial unfolding follows the unfolding pathway described previously for higher-temperature simulations: the hydrophobic core is disrupted, the beta-sheet pulls apart and the alpha-helix unravels. The unfolded state reached under these conditions maintains a kernel of structure in the form of a non-native hydrophobic cluster. Refolding simply reverses this path, the side-chain interactions shift, the helix refolds, and the native packing and hydrogen bonds are recovered. The end result of this refolding is not the initial crystal structure; it contains the proper topology and the majority of the native contacts, but the structure is expanded and the contacts are long. We believe this to be the native state at elevated temperature, and the change in volume and contact lengths is consistent with experimental studies of other native proteins at elevated temperature and the chemical denaturant equivalent of T(m).     

Equilibrium and kinetic studies on folding of the authentic and recombinant forms of human alpha-lactalbumin by circular dichroism spectroscopy

The equilibrium and kinetics of the unfolding and refolding of authentic and recombinant human alpha-lactalbumin, the latter of which had an extra methionine residue at the N-terminus, were studied by circular dichroism spectroscopy, and the results were compared with the results for bovine and goat alpha-lactalbumins obtained in our previous studies. As observed in the bovine and goat proteins, the presence of the extra methionine residue in the recombinant protein remarkably destabilized the native state, and the destabilization was entirely ascribed to an increase in the rate of unfolding. The thermodynamic stability of the native state against the unfolded state was lower, and the thermodynamic stability of the molten globule state against the unfolded state was higher for the human protein than for the other alpha-lactalbumins previously studied. Thus, the population of the molten globule intermediate was higher during the equilibrium unfolding of human alpha-lactalbumin by guanidine hydrochloride. Unlike the molten globule states of the bovine and goat proteins, the human alpha-lactalbumin molten globule showed remarkably more intense circular dichroism ellipticity than the native state in the far-ultraviolet region below 225 nm. During refolding from the unfolded state, human alpha-lactalbumin thus exhibited overshoot kinetics, in which the alpha-helical peptide ellipticity exceeded the native value when the molten globule folding intermediate was formed in the burst phase. The subsequent folding involved reorganization of nonnative secondary structures. It should be noted that the rate constant of the major refolding phase was approximately the same among the three types of alpha-lactalbumin and that the rate constant of unfolding was accelerated 18-600 times in the human protein, and these results interpreted the lower thermodynamic stability of this protein.     

Protein Folding Mechanism of the Dimeric AmphiphysinII/Bin1 N-BAR Domain

The human AmphyphisinII/Bin1 N-BAR domain belongs to the BAR domain superfamily, whose members sense and generate membrane curvatures. The N-BAR domain is a 57 kDa homodimeric protein comprising a six helix bundle. Here we report the protein folding mechanism of this protein as a representative of this protein superfamily. The concentration dependent thermodynamic stability was studied by urea equilibrium transition curves followed by fluorescence and far-UV CD spectroscopy. Kinetic unfolding and refolding experiments, including rapid double and triple mixing techniques, allowed to unravel the complex folding behavior of N-BAR. The equilibrium unfolding transition curve can be described by a two-state process, while the folding kinetics show four refolding phases, an additional burst reaction and two unfolding phases. All fast refolding phases show a rollover in the chevron plot but only one of these phases depends on the protein concentration reporting the dimerization step. Secondary structure formation occurs during the three fast refolding phases. The slowest phase can be assigned to a proline isomerization. All kinetic experiments were also followed by fluorescence anisotropy detection to verify the assignment of the dimerization step to the respective folding phase. Based on these experiments we propose for N-BAR two parallel folding pathways towards the homodimeric native state depending on the proline conformation in the unfolded state.     

Unfolding rates of barstar determined in native and low denaturant conditions indicate the presence of intermediates

The folding and unfolding rates of the small protein, barstar, have been monitored using stopped-flow measurements of intrinsic tryptophan fluorescence at 25 degrees C, pH 8.5, and have been compared over a wide range of urea and guanidine hydrochloride (GdnHCl) concentrations. When the logarithms of the rates of folding from urea and from GdnHCl unfolded forms are extrapolated linearly with denaturant concentration, the same rate is obtained for folding in zero denaturant. Similar linear extrapolations of rates of unfolding in urea and GdnHCl yield, however, different unfolding rates in zero denaturant, indicating that such linear extrapolations are not valid. It has been difficult, for any protein, to determine unfolding rates under nativelike conditions in direct kinetic experiments. Using a novel strategy of coupling the reactivity of a buried cysteine residue with 5,5'-dithiobis(2-nitrobenzoic acid) (DTNB) to the unfolding reaction of barstar, the global unfolding and refolding rates have now been determined in low denaturant concentrations. The logarithms of unfolding rates obtained at low urea and GdnHCl concentrations show a markedly nonlinear dependence on denaturant concentration and converge to the same unfolding rate in the absence of denaturant. It is shown that the native protein can sample the fully unfolded conformation even in the absence of denaturant. The observed nonlinear dependences of the logarithms of the refolding and unfolding rates observed for both denaturants are shown to be due to the presence of (un)folding intermediates and not due to movements in the position of the transition state with a change in denaturant concentration.     

High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

Energetics of protein structure and folding

The available experimental date on the kinetics of unfolding and refolding of small proteins are reviewed. Excluding slow transitions in the unfolded protein due to cis–trans isomerization of peptide bonds, the rate-limiting transition state in both unfolding and refolding is concluded to be a high-energy distortion of the fully folded state. Partially folded intermediates are undoubtedly important for folding, but their formation is normally not rate limiting. A simple model is used to illustrate some of the aspects of protein-folding energetics.     

Charge-charge interactions in the denatured state influence the folding kinetics of ribonuclease Sa

Gaining a better understanding of the denatured state ensemble of proteins is important for understanding protein stability and the mechanism of protein folding. We studied the folding kinetics of ribonuclease Sa (RNase Sa) and a charge-reversal variant (D17R). The refolding kinetics are similar, but the unfolding rate constant is 10-fold greater for the variant. This suggests that charge-charge interactions in the denatured state and the transition state ensembles are more favorable in the variant than in RNase Sa, and shows that charge-charge interactions can influence the kinetics and mechanism of protein folding.     

Is protein unfolding the reverse of protein folding? A lattice simulation analysis.

Simulations and experiments that monitor protein unfolding under denaturing conditions are commonly employed to study the mechanism by which a protein folds to its native state in a physiological environment. Due to the differences in conditions and the complexity of the reaction, unfolding is not necessarily the reverse of folding. To assess the relevance of temperature initiated unfolding studies to the folding problem, we compare the folding and unfolding of a 125-residue protein model by Monte Carlo dynamics at two temperatures; the lower one corresponds to the range used in T -jump experiments and the higher one to the range used in unfolding simulations of all-atom models. The trajectories that lead from the native state to the denatured state at these elevated temperatures are less diverse than those observed in the folding simulations. At the lower temperature, the system unfolds through a mandatory intermediate that corresponds to a local free energy minimum. At the higher temperature, no such intermediate is observed, but a similar pathway is followed. The structures contributing to the unfolding pathways resemble most closely those that make up the "fast track" of folding. The transition state for unfolding at the lower temperature (above Tm) is determined and is found to be more structured than the transition state for folding below the melting temperature. This shift towards the native state is consistent with the Hammond postulate. The implications for unfolding simulations of higher resolution models and for unfolding experiments of proteins are discussed.     

Reversible denaturation of the gene V protein of bacteriophage f1

The guanidine hydrochloride (GuHCl)-induced denaturation of the gene V protein of bacteriophage f1 has been studied, using the chemical reactivity of a cysteine residue that is buried in the folded protein and the circular dichroism (CD) at 211 and 229 nm as measures of the fraction of polypeptide chains in the folded form. It is found that this dimeric protein unfolds in a single cooperative transition from a folded dimer to two unfolded monomers. A folded, monomeric form of the gene V protein was not detected at equilibrium. The kinetics of unfolding of the gene V protein in 3 M GuHCl and the refolding in 2 M GuHCl are also consistent with a transition between a folded dimer and two unfolded monomers. The GuHCl concentration dependence of the rates of folding and unfolding suggests that the transition state for folding is near the folded conformation.     

Characterization of the folding and unfolding reactions of single-chain monellin: evidence for multiple intermediates and competing pathways

The mechanisms of folding and unfolding of the small plant protein monellin have been delineated in detail. For this study, a single-chain variant of the natively two-chain monellin, MNEI, was used, in which the C terminus of chain B was connected to the N terminus of chain A by a Gly-Phe linker. Equilibrium guanidine hydrochloride (GdnHCl)-induced unfolding experiments failed to detect any partially folded intermediate that is stable enough to be populated at equilibrium to a significant extent. Kinetic experiments in which the refolding of GdnHCl-unfolded protein was monitored by measurement of the change in the intrinsic tryptophan fluorescence of the protein indicated the accumulation of three transient partially structured folding intermediates. The fluorescence change occurred in three kinetic phases: very fast, fast, and slow. It appears that the fast and slow changes in fluorescence occur on competing folding pathways originating from one unfolded form and that the very fast change in fluorescence occurs on a third parallel pathway originating from a second unfolded form of the protein. Kinetic experiments in which the refolding of alkali-unfolded protein was monitored by the change in the fluorescence of the hydrophobic dye 8-anilino-1-naphthalenesulfonic acid (ANS), consequent to the dye binding to the refolding protein, as well as by the change in intrinsic tryptophan fluorescence, not only confirmed the presence of the three kinetic intermediates but also indicated the accumulation of one or more early intermediates at a few milliseconds of refolding. These experiments also exposed a very slow kinetic phase of refolding, which was silent to any change in the intrinsic tryptophan fluorescence of the protein. Hence, the spectroscopic studies indicated that refolding of single-chain monellin occurs in five distinct kinetic phases. Double-jump, interrupted-folding experiments, in which the accumulation of folding intermediates and native protein during the folding process could be determined quantitatively by an unfolding assay, indicated that the fast phase of fluorescence change corresponds to the accumulation of two intermediates of differing stabilities on competing folding pathways. They also indicated that the very slow kinetic phase of refolding, identified by ANS binding, corresponds to the formation of native protein. Kinetic experiments in which the unfolding of native protein in GdnHCl was monitored by the change in intrinsic tryptophan fluorescence indicated that this change occurs in two kinetic phases. Double-jump, interrupted-unfolding experiments, in which the accumulation of unfolding intermediates and native protein during the unfolding process could be determined quantitatively by a refolding assay, indicated that the fast unfolding phase corresponds to the formation of fully unfolded protein via one unfolding pathway and that the slow unfolding phase corresponds to a separate unfolding pathway populated by partially unfolded intermediates. It is shown that the unfolded form produced by the fast unfolding pathway is the one which gives rise to the very fast folding pathway and that the unfolded form produced by the slower unfolding pathway is the one which gives rise to the slow and fast folding pathways.     

Test for cooperativity in the early kinetic intermediate in lysozyme folding

During the folding of many proteins, collapsed globular states are formed prior to the native structure. The role of these states for the folding process has been widely discussed. Comparison with properties of synthetic homo and heteropolymers had suggested that the initial collapse represented a shift of the ensemble of unfolded conformations to more compact states without major energy barriers. We investigated the folding/unfolding transition of a collapsed state, which transiently populates early in lysozyme folding. This state forms within the dead-time of stopped-flow mixing and it has been shown to be significantly more compact and globular than the denaturant-induced unfolded state. We used the GdmCl-dependence of the dead-time signal change to characterize the unfolding transition of the burst phase intermediate. Fluorescence and far-UV CD give identical unfolding curves, arguing for a cooperative two-state folding/unfolding transition between unfolded and collapsed lysozyme. These results show that collapse leads to a distinct state in the folding process, which is separated from the ensemble of unfolded molecules by a significant energy barrier. NMR, fluorescence and small angle X-ray scattering data further show that some local interactions in unfolded lysozyme exist at denaturant concentrations above the coil-collapse transition. These interactions might play a crucial role in the kinetic partitioning between fast and slow folding pathways.