Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.     

Dimensions,energetics, and denaturant effects of the protein unstructured state

Determining the energetics of the unfolded state of a protein is essential for understanding the folding mechanics of ordered proteins and the structure–function relation of intrinsically disordered proteins. Here, we adopt a coil‐globule transition theory to develop a general scheme to extract interaction and free energy information from single‐molecule fluorescence resonance energy transfer spectroscopy. By combining protein stability data, we have determined the free energy difference between the native state and the maximally collapsed denatured state in a number of systems, providing insight on the specific/nonspecific interactions in protein folding. Both the transfer and binding models of the denaturant effects are demonstrated to account for the revealed linear dependence of inter‐residue interactions on the denaturant concentration, and are thus compatible under the coil‐globule transition theory to further determine the dimension and free energy of the conformational ensemble of the unfolded state. The scaling behaviors and the effective θ‐state are also discussed.     

Hammond behavior versus ground state effects in protein folding: evidence for narrow free energy barriers and residual structure in unfolded states

Apparent transition state movement upon mutation or changes in solvent conditions is frequently observed in protein folding and is often interpreted in terms of Hammond behavior. This led to the conclusion that barrier regions in protein folding are broad maxima on the free energy landscape. Here, we use the concept of self-interaction and cross-interaction parameters to test experimental data of 21 well-characterized proteins for Hammond behavior. This allows us to characterize the origin of transition state movements along different reaction coordinates. Only one of the 21 proteins shows a small but coherent transition state movement in agreement with the Hammond postulate. In most proteins the structure of the transition state is insensitive to changes in protein stability. The apparent change in the position of the transition state upon mutation, which is frequently observed in phi-value analysis, is in most cases due to ground-state effects caused by structural changes in the unfolded state. This argues for significant residual structure in unfolded polypeptide chains of many proteins. Disruption of these residual interactions by mutation often leads to decreased folding rates, which implies that these interactions are still present in the transition state. The failure to detect Hammond behavior shows that the free energy barriers encountered by a folding polypeptide chain are generally rather narrow and robust maxima for all experimentally explorable reaction coordinates.     

Structural dynamics, stability and folding of proteins

The present concepts of protein folding in vitro are reviewed. According to these concepts, amino acid sequence of protein, which has appeared a result of evolutionary selection, determines the native structure of protein, the pathway of protein folding, and the existence of free energy barrier between native and denatured states of protein. The latter means that protein macromolecule can exist in either native or denatured state. And all macromolecules in the native state are identical but for structural fluctuations due to Brownian motion of their atoms. Identity of all molecules in native state is of primary importance for their correct functioning. The dependence of protein stability, which is measured as the difference between free energy of protein in native and denatured states, on temperature and denaturant concentration is discussed. The modern approaches characterizing transition state and nucleation are regarded. The role of intermediate and misfolded states in amorphous aggregate and amyloid fibril formation is discussed.     

Unfolding of the cold shock protein studied with biased molecular dynamics

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.     

Snapshots of a dynamic folding nucleus in zinc-substituted Pseudomonas aeruginosa azurin

Zinc-substituted Pseudomonas aeruginosa azurin folds in two-state equilibrium and kinetic reactions. In the unfolded state, the zinc ion remains bound to the unfolded polypeptide via two native-state ligands (His117 and Cys112). The significantly curved Chevron plot for zinc-substituted azurin was earlier ascribed to movement of the folding-transition state. At low concentrations of denaturant, the transition state occurs early in the folding reaction (low Tanford beta-value), whereas at high-denaturant concentration, it moves closer to the native structure (high Tanford beta-value). Here, we use this movement to track the formation and growth of zinc-substituted azurin's folding nucleus with atomic resolution using protein engineering. The average phi (phi) value for 17 positions (covering all secondary-structure elements) goes from 0.25 in 0 M GuHCl (beta approximately 0.46) to 0.76 in 4 M GuHCl (beta approximately 0.86); a phi-value of 1 or 0 indicates native-like or unfolded-like interactions, respectively. Analysis of individual phi-values reveals a delocalized nucleus where structure condenses around a leading density centered on Leu50 in the core. The diffuse moving transition state for zinc-substituted azurin is in sharp contrast to the fixed polarized folding nucleus observed for apo-azurin. The dramatic difference in apparent kinetic behavior for the two forms of azurin can be rationalized as a minor alteration on a common free-energy profile that exhibits a broad activation barrier.     

Key role of coulombic interactions for the folding transition state of the cold shock protein

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.     

Ionic-strength-dependent effects in protein folding: analysis of rate equilibrium free-energy relationships and their interpretation

The traditional approach to studying protein folding involves applying a perturbation, usually denaturant or mutation, and determining the effect upon the free energy of folding, DeltaG0, and the activation free energy, DeltaG(not equal). Data collected as a function of the perturbation can be used to construct rate equilibrium free-energy relationships, which report on the development of interactions in the transition state for folding. We examine the use of the ionic-strength-dependent rate equilibrium free-energy relationship in protein folding using the N-terminal domain of L9, a small alpha-beta protein, as a model system. Folding is two-state for the range of ionic strength examined, 0.045-1.52 M. The plot of DeltaG(not equal) versus DeltaG0 is linear (r2= 0.918), with a slope equal to 0.45. The relatively low value of the slope indicates that the ionic-strength-dependent interactions are modestly developed in the transition state. The slope is, however, greater than that of a plot of DeltaG(not equal) versus DeltaG0 constructed by varying pH, thus demonstrating directly that ionic-strength-dependent studies probe more than simple electrostatic interactions. Potential transition movement was probed by analysis of the denaturant, ionic strength cross-interaction parameters. The values are small but nonzero and positive, suggesting a small shift of the transition state toward the native state as the protein is destabilized, i.e., Hammond behavior. The complications that arise in the interpretation of ionic-strength-dependent rate equilibrium free-energy relationships are discussed, and it is concluded that the ionic-strength-dependent studies do not provide a reliable indicator of the role of electrostatic interactions. Complications include incomplete screening of electrostatic interactions, specific ion binding, Hofmeister effects, and the potential presence of electrostatic interactions in the denatured state ensemble.     

Complexity of aromatic ring-flip motions in proteins: Y97 ring dynamics in cytochrome c observed by cross-relaxation suppressed exchange NMR spectroscopy

Dynamics of large-amplitude conformational motions in proteins are complex and less understood, although these processes are intimately associated with structure, folding, stability, and function of proteins. Here, we use a large set of spectra obtained by cross-relaxation suppressed exchange NMR spectroscopy (EXSY) to study the 180° flipping motion of the Y97 ring of horse ferricytochrome c as a function of near-physiological temperature in the 288–308 K range. With rising temperature, the ring-flip rate constant makes a continuous transition from Arrhenius to anti-Arrhenius behavior through a narrow Arrhenius-like zone. This behavior is seen not only for the native state of the protein, but also for native-like states generated by adding subdenaturing amounts of guanidine deuterochloride (GdnDCl). Moderately destabilizing concentrations of the denaturant (1.5 M GdnDCl) completely removes the Arrhenius-like feature from the temperature window employed. The Arrhenius to anti-Arrhenius transition can be explained by the heat capacity model where temperature strengthens ground state interactions, perhaps hydrophobic in nature. The effect of the denaturant may appear to arise from direct protein-denaturant interactions that are structure-stabilizing under subdenaturing conditions. The temperature distribution of rate constants under different stability conditions also suggests that the prefactor in Arrhenius-like relations is temperature dependent. Although the use of the transition state theory (TST) offers several challenges associated with data interpretation, the present results and a consideration of others published earlier provide evidence for complexity of ring-flip dynamics in proteins.     

Folding of a de novo designed native-like four-helix bundle protein

The folding of a model native-like dimeric four-helix bundle protein, (alpha(2))(2), was investigated using guanidine hydrochloride, hydrostatic pressure, and low temperature. Unfolding by guanidine hydrochloride followed by circular dichroism and intrinsic fluorescence spectroscopy revealed a highly cooperative transition between the native-like and unfolded states, with free energy of unfolding determined from CD data, DeltaG(unf) = 14.3 +/- 0.8 kcal/mol. However, CD and intrinsic fluorescence data were not superimposable, indicating the presence of an intermediate state during the folding transition. To stabilize the folding intermediate, we used hydrostatic pressure and low temperature. In both cases, dissociation of the dimeric native-like (alpha(2))(2) into folded monomers (alpha(2)) was observed. van't Hoff analysis of the low temperature experiments, assuming a two-state dimer 171-monomer transition, yielded a free energy of dissociation of (alpha(2))(2) of DeltaG(diss) = 11.4 +/- 0.4 kcal/mol, in good agreement with the free energy determined from pressure dissociation experiments (DeltaG(diss) = 10.5 +/- 0.1 kcal/mol). Binding of the hydrophobic fluorescent probe 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid (bis-ANS) to the pressure- and cold-dissociated states of (alpha(2))(2) indicated the existence of molten-globule monomers. In conclusion, we demonstrate that the folding pathway of (alpha(2))(2) can be described by a three-state transition including a monomeric molten globule-like state.     

Temperature-dependent Hammond behavior in a protein-folding reaction: analysis of transition-state movement and ground-state effects

Characterization of the transition-state ensemble and the nature of the free-energy barrier for protein folding are areas of intense activity and some controversy. A key issue that has emerged in recent years is the width of the free-energy barrier and the susceptibility of the transition state to movement. Here we report denaturant-induced and temperature-dependent folding studies of a small mixed α-β protein, the N-terminal domain of L9 (NTL9). The folding of NTL9 was determined using fluorescence-detected stopped-flow fluorescence measurements conducted at seven different temperatures between 11 and 40 °C. Plots of the log of the observed first-order rate constant versus denaturant concentration, “chevron plots,” displayed the characteristic V shape expected for two-state folding. There was no hint of deviation from linearity even at the lowest denaturant concentrations. The relative position of the transition state, as judged by the Tanford β parameter, β_T, shifts towards the native state as the temperature is increased. Analysis of the temperature dependence of the kinetic and equilibrium m values indicates that the effect is due to significant movement of the transition state and also includes a contribution from temperature-dependent ground-state effects. Analysis of the Leffler plots, plots of ΔG^‡versus ΔG°, and their cross-interaction parameters confirms the transition-state movement. Since the protein is destabilized at high temperature, the shift represents a temperature-dependent Hammond effect. This provides independent confirmation of a recent theoretical prediction. The magnitude of the temperature-denaturant cross-interaction parameter is larger for NTL9 than has been reported for the few other cases studied. The implications for temperature-dependent studies of protein folding are discussed.     

Sequential barriers and an obligatory metastable intermediate define the apparent two-state folding pathway of the ubiquitin-like PB1 domain of NBR1

The 90-residue N-terminal Phox and Bem1p (PB1) domain of NBR1 forms an α/β ubiquitin-like fold. Kinetic analysis using stopped-flow fluorescence reveals two-state kinetics; however, nonlinear effects in the denaturant dependence of the unfolding data demonstrate changes in the position of the rate-limiting barrier along the folding coordinate as the folding conditions change. The kinetics of wt-PB1 and several mutants show that this curvature is consistent with a single-pathway mechanism involving sequential transition states (TS1 and TS2) separated by a transiently populated high-energy intermediate, rather than movement of the transition state on a broad energy plateau. We show that the two transition states within the sequential model represent structurally and thermodynamically distinct species. TS1 is a collapsed state (α_TS1 = 0.71) with a large enthalpic barrier to formation that is rate-limiting under conditions that strongly favour folding. TS2 is highly native-like (α_TS2 = 0.93) and represents a late entropic barrier to formation of the native state. In support of the sequential transition state mechanism, we show that the G62A helix 2 substitution stabilises TS1 and the intermediate to such an extent that the latter becomes significantly populated, leading to the observation of a fast kinetic phase representing the initial U → I transition, with TS2 (α_TS2 = 0.87) becoming rate-limiting. The folding rate is not retarded by populating an intermediate, which would be expected for a misfold state, but is accelerated, suggesting that the I state is productive and on-pathway. The results show that the apparent two-state folding of the wt-PB1 domain occurs along a well-defined pathway involving structurally and thermodynamically distinct sequential transition states and an obligatory metastable intermediate that represents a productive local minimum in the energy landscape that increases the efficiency of barrier crossing through favourable effects on the entropy of activation.     

Sensitivity of the folding/unfolding transition state ensemble of chymotrypsin inhibitor 2 to changes in temperature and solvent

To better characterize the transition state for folding/unfolding and its sensitivity to environmental changes, we have run multiple molecular dynamics simulations of chymotrypsin inhibitor 2 (CI2) under varying solvent conditions and temperature. The transition state structures agree well with experiment, and are similar under all of the conditions investigated here. Increasing the temperature leads to some movement in the position of the transition state along several reaction coordinates, as measured by changes in properties of the transition state structures. These structural changes are in the direction of a more native-like transition state as denaturation conditions become more severe, as expected for a Hammond effect. These structural changes are not, however, reflected in the global structure as measured by the total number of contacts or the average S-values. These results suggest that the small changes in average Phi-values with temperature seen by experiment may be due to an increase in the sensitivity of the transition state to mutation rather than a change in the average structure of the transition state. A simple analysis of the rates of unfolding indicates that the free energy barrier to unfolding decreases with increasing temperature, but even in our very high temperature simulations there is a small free energy barrier.     

Characterization of the critical state in protein folding. Effects of guanidine hydrochloride and specific Ca2+ binding on the folding kinetics of alpha-lactalbumin

The reversible unfolding and refolding kinetics of alpha-lactalbumin induced by concentration jump of guanidine hydrochloride were measured at pH 7.0 and 25 degrees C using tryptophan absorption at 292 nm, with varying concentrations of the denaturant and free Ca2+. The refolding reaction of alpha-lactalbumin from the fully unfolded (D) state occurs through the two stages: (1) instantaneous formation of a compact intermediate (the A state) that has a native-like secondary structure; (2) tight packing of the preformed secondary structure segments to lead finally to the native structure, this stage being the rate-determining step of the reaction and associated with acquisition of the specific structure necessary for strong Ca2+ binding. Under strongly native conditions, the observed kinetics of refolding is also complicated by the presence of a slow-folding species (10%) in the unfolded state. Considering these facts, the microscopic rate constants in folding and unfolding directions have been evaluated from the observed kinetics and from the equilibrium constants of the transitions among the native (N), A and D states. Close linear relationships have been found in the plots of the activation free energies, obtained from the microscopic rate constants, against the denaturant concentration. They are similar to the linear relationship between the free energy of unfolding and the denaturant concentration. It was demonstrated that the slope of the plots should be approximately proportional to a change in accessible surface area of the protein during the respective activation process, and that only a third of the difference in accessible surface area between A and N is buried in the critical activated state of folding. However, the selective effect of Ca2+ binding on the folding rate constant has been observed also, demonstrating that the specific Ca2+-binding substructure in the N state is already organized in the activated state. Thus, only a part of the protein molecule involving the Ca2+-binding region is organized in the activated state, with the other part of the molecule being left less organized, suggesting that the second stage of folding may be a sequential growing process of organized assemblage of the performed secondary structure segments.     

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.     

Correlation between protein stability cores and protein folding kinetics: a case study on Pseudomonas aeruginosa apo-azurin

This paper reports a combined computational and experimental study of the correlation between protein stability cores and folding kinetics. An empirical potential function was developed, and it was used for analyzing interaction energies among secondary structure elements. Studies on a beta sandwich protein, Pseudomonas aeruginosa azurin, showed that the computationally identified substructure with the strongest interactions in the native state is identical to the "interlocked pair" of beta strands, an invariant motif found in most sandwich-like proteins. Moreover, previous and new in vitro folding results revealed that the identified substructure harbors most residues that form native-like interactions in the folding transition state. These observations demonstrate that the potential function is effective in revealing the relative strength of interactions among various protein parts; they also strengthen the suggestion that the most stable regions in native proteins favor stable interactions early during folding.     

Reliable protein folding on complex energy landscapes: the free energy reaction path

A theoretical framework is developed to study the dynamics of protein folding. The key insight is that the search for the native protein conformation is influenced by the rate r at which external parameters, such as temperature, chemical denaturant, or pH, are adjusted to induce folding. A theory based on this insight predicts that 1), proteins with complex energy landscapes can fold reliably to their native state; 2), reliable folding can occur as an equilibrium or out-of-equilibrium process; and 3), reliable folding only occurs when the rate r is below a limiting value, which can be calculated from measurements of the free energy. We test these predictions against numerical simulations of model proteins with a single energy scale.     

A polar, solvent-exposed residue can be essential for native protein structure

BACKGROUND: A large energy gap between the native state and the non-native folded states is required for folding into a unique three-dimensional structure. The features that define this energy gap are not well understood, but can be addressed using de novo protein design. Previously, alpha(2)D, a dimeric four-helix bundle, was designed and shown to adopt a native-like conformation. The high-resolution solution structure revealed that this protein adopted a bisecting U motif. Glu7, a solvent-exposed residue that adopts many conformations in solution, might be involved in defining the unique three-dimensional structure of alpha(2)D. RESULTS: A variety of hydrophobic and polar residues were substituted for Glu7 and the dynamic and thermodynamic properties of the resulting proteins were characterized by analytical ultracentrifugation, circular dichroism spectroscopy, and nuclear magnetic resonance spectroscopy. The majority of substitutions at this solvent-exposed position had little affect on the ability to fold into a dimeric four-helix bundle. The ability to adopt a unique conformation, however, was profoundly modulated by the residue at this position despite the similar free energies of folding of each variant. CONCLUSIONS: Although Glu7 is not involved directly in stabilizing the native state of alpha(2)D, it is involved indirectly in specifying the observed fold by modulating the energy gap between the native state and the non-native folded states. These results provide experimental support for hypothetical models arising from lattice simulations of protein folding, and underscore the importance of polar interfacial residues in defining the native conformations of proteins.     

Searching for multiple folding pathways of a nearly symmetrical protein: temperature dependent phi-value analysis of the B domain of protein A

The B domain of protein A (BdpA) is a popular paradigm for simulating protein folding pathways. The discrepancies between so many simulations and subsequent experimental testing may be attributable to the protein being highly symmetrical: changing experimental conditions could perturb the subtle interplay between the effects of symmetry in the native structure and the effects of asymmetry from specific interactions in a given sequence. If the protein folds via multiple pathways, perturbations, such as temperature, denaturant concentration, and mutation, should change the flux of micro pathways, leading to changes in the bulk properties of the transition state. We tested this hypothesis by conducting a Phi-analysis of BdpA as a function of temperature from 25.0 degrees C to 60.0 degrees C. The Phi-values had no significant dependence on temperature and the values at 55.0 degrees C (denaturing conditions) are very similar to those at 25.0 degrees C (folding conditions), indicating the structure of the transition state does not significantly change although the experimental conditions are considerably altered. The results suggest that BdpA folds via a single dominant folding pathway.