Two-state expansion and collapse of a polypeptide

The initial phase of folding for many proteins is presumed to be the collapse of the polypeptide chain from expanded to compact, but still denatured, conformations. Theory and simulations suggest that this collapse may be a two-state transition, characterized by barrier-crossing kinetics, while the collapse of homopolymers and random heteropolymers is continuous and multi-phasic. A new rapid-mixing flow technique has been used to resolve the late stages of polypeptide collapse, at time scales >/=45 microseconds. We have used a laser temperature-jump with fluorescence spectroscopy to resolve the complete time-course of the collapse of denatured cytochrome c with nanosecond time resolution. We find the process to be exponential in time and thermally activated, with an apparent activation energy approximately 9 k(B)T (after correction for solvent viscosity). These results indicate that polypeptide collapse is kinetically a two-state transition. Because of the observed free energy barrier, the time scale of polypeptide collapse is dramatically slower than is predicted by Langevin models for homopolymer collapse.     

Folding under inequilibrium conditions as a possible reason for partial irreversibility of heat-denatured proteins: computer simulation study

Using computer simulations we have studied possible effects of heating and cooling at different scan rates on unfolding and refolding of macromolecules. We have shown that even the simplest two-state reversible transition can behave irreversibly when an unfavorable combination of cooling rate, relaxation time and activation energy of refolding occurs. On the basis of this finding we suppose that apparent irreversibility of some proteins denatured by heat may result from slow relaxation on cooling rather than thermodynamic instability and/or irreversible alterations of the polypeptide chain. Using this kinetic reversible two-state model, we estimated the effects of the scan rate and kinetic parameters of the macromolecule on its unfolding-refolding process. A few recommendations are suggested on how to reach maximal possible recovery after denaturation if refolding appears to be under kinetic control.     

The dynamical contact order: Protein folding rate parameters based on quantum conformational transitions

Protein folding is regarded as a quantum transition between the torsion states of a polypeptide chain. According to the quantum theory of conformational dynamics, we propose the dynamical contact order (DCO) defined as a characteristic of the contact described by the moment of inertia and the torsion potential energy of the polypeptide chain between contact residues. Consequently, the protein folding rate can be quantitatively studied from the point of view of dynamics. By comparing theoretical calculations and experimental data on the folding rate of 80 proteins, we successfully validate the view that protein folding is a quantum conformational transition. We conclude that (i) a correlation between the protein folding rate and the contact inertial moment exists; (ii) multi-state protein folding can be regarded as a quantum conformational transition similar to that of two-state proteins but with an intermediate delay. We have estimated the order of magnitude of the time delay; (iii) folding can be classified into two types, exergonic and endergonic. Most of the two-state proteins with higher folding rate are exergonic and most of the multi-state proteins with low folding rate are endergonic. The folding speed limit is determined by exergonic folding.     

Calorimetric study of the heat and cold denaturation of beta-lactoglobulin.

Temperature-induced changes of the states of beta-lactoglobulin have been studied calorimetrically. In the presence of a high concentration of urea this protein shows not only heat but also cold denaturation. Its heat denaturation is approximated very closely by a two-state transition, while the cold denaturation deviates considerably from the two-state transition and this deviation increases as the temperature decreases. The heat effect of cold denaturation is opposite in sign to that of heat denaturation and is noticeably larger in magnitude. This difference in magnitude is caused by the temperature-dependent negative heat effect of additional binding of urea to the polypeptide chain of the protein upon its unfolding, which decreases the positive enthalpy of heat denaturation and increases the negative enthalpy of cold denaturation. The binding of urea considerably increases the partial heat capacity of the protein, especially in the denatured state. However, when corrected for the heat capacity effect of urea binding, the partial heat capacity of the denatured protein is close in magnitude to that expected for the unfolded polypeptide chain in aqueous solution without urea but only for temperatures below 10 degrees C. At higher temperatures, the heat capacity of the denatured protein is lower than that expected for the unfolded polypeptide chain. It appears that at temperatures above 10 degrees C not all the surface of the beta-lactoglobulin polypeptide chain is exposed to the solvent, even in the presence of 6 M urea; i.e., the denatured protein is not completely unfolded and unfolds only at temperatures lower than 10 degrees C.(ABSTRACT TRUNCATED AT 250 WORDS)     

Computing the transition state populations in simple protein models

We describe the master equation method for computing the kinetics of protein folding. We illustrate the method using a simple Go model. Presently most models of two-state fast-folding protein folding kinetics invoke the classical idea of a transition state to explain why there is a single exponential decay in time. However, if proteins fold via funnel-shaped energy landscapes, as predicted by many theoretical studies, then it raises the question of what is the transition state. Is it a specific structure, or a small ensemble of structures, as is expected from classical transition state theory? Or is it more like the denatured states of proteins, a very broad ensemble? The answer that is usually obtained depends on the assumptions made about the transition state. The present method is a rigorous way to find transition states, without assumptions or approximations, even for very nonclassical shapes of energy landscapes. We illustrate the method here, showing how the transition states in two-state protein folding can be very broad ensembles. © 2002 Wiley Periodicals, Inc. Biopolymers 68: 35–46, 2003     

Water and urea interactions with the native and unfolded forms of a beta-barrel protein

A fundamental understanding of protein stability and the mechanism of denaturant action must ultimately rest on detailed knowledge about the structure, solvation, and energetics of the denatured state. Here, we use (17)O and (2)H magnetic relaxation dispersion (MRD) to study urea-induced denaturation of intestinal fatty acid-binding protein (I-FABP). MRD is among the few methods that can provide molecular-level information about protein solvation in native as well as denatured states, and it is used here to simultaneously monitor the interactions of urea and water with the unfolding protein. Whereas CD shows an apparently two-state transition, MRD reveals a more complex process involving at least two intermediates. At least one water molecule binds persistently (with residence time >10 nsec) to the protein even in 7.5 M urea, where the large internal binding cavity is disrupted and CD indicates a fully denatured protein. This may be the water molecule buried near the small hydrophobic folding core at the D-E turn in the native protein. The MRD data also provide insights about transient (residence time <1 nsec) interactions of urea and water with the native and denatured protein. In the denatured state, both water and urea rotation is much more retarded than for a fully solvated polypeptide. The MRD results support a picture of the denatured state where solvent penetrates relatively compact clusters of polypeptide segments.     

An expanded view of the protein folding landscape of PDZ domains

Most protein domains fold in an apparently co-operative and two-state manner with only the native and denatured states significantly populated at any experimental condition. However, the protein folding energy landscape is often rugged and different transition states may be rate limiting for the folding reaction under different conditions, as seen for the PDZ protein domain family. We have here analyzed the folding kinetics of two PDZ domains and found that a previously undetected third transition state is rate limiting under conditions that stabilize the native state relative to the denatured state. In light of these results, we have re-analyzed previous folding data on PDZ domains and present a unified folding mechanism with three distinct transition states separated by two high-energy intermediates. Our data show that sequence composition tunes the relative stabilities of folding transition states within the PDZ family, while the overall mechanism is determined by topology. This model captures the kinetic folding mechanism of all PDZ domains studied to date.     

Polymer principles of protein calorimetric two-state cooperativity

The experimental calorimetric two-state criterion requires the van't Hoff enthalpy DeltaH(vH) around the folding/unfolding transition midpoint to be equal or very close to the calorimetric enthalpy DeltaH(cal) of the entire transition. We use an analytical model with experimental parameters from chymotrypsin inhibitor 2 to elucidate the relationship among several different van't Hoff enthalpies used in calorimetric analyses. Under reasonable assumptions, the implications of these DeltaH(vH)'s being approximately equal to DeltaH(cal) are equivalent: Enthalpic variations among denatured conformations in real proteins are much narrower than some previous lattice-model estimates, suggesting that the energy landscape theory "folding to glass transition temperature ratio" T(f) /T(g) may exceed 6.0 for real calorimetrically two-state proteins. Several popular three-dimensional lattice protein models, with different numbers of residue types in their alphabets, are found to fall short of the high experimental standard for being calorimetrically two-state. Some models postulate a multiple-conformation native state with substantial pre-denaturational energetic fluctuations well below the unfolding transition temperature, or predict a significant post-denaturational continuous conformational expansion of the denatured ensemble at temperatures well above the transition point, or both. These scenarios either disagree with experiments on protein size and dynamics, or are inconsistent with conventional interpretation of calorimetric data. However, when empirical linear baseline subtractions are employed, the resulting DeltaH(vH)/DeltaH(cal)'s for some models can be increased to values closer to unity, and baseline subtractions are found to correspond roughly to an operational definition of native-state conformational diversity. These results necessitate a re-assessment of theoretical models and experimental interpretations.     

Is it always possible to distinguish two- and three-state systems by evaluating the van't Hoff enthalpy?

Many small globular proteins are traditionally classified as thermodynamical two-state systems, i.e., the protein is either in the native, active state (folded) or in the denatured state (unfolded). We challenge this view and show that there may exist (protein) systems for which a van't Hoff analysis of experimental data cannot determine whether the system corresponds to two or three thermodynamical states when only temperatures in a narrow temperature region around the transition are considered. We generalize a widely employed two-state protein folding model to include a third, transition state. For this three-state system we systematically study the deviation of the calorimetric enthalpy (heat of transition) from the van't Hoff enthalpy, a measure of the two-stateness of a transition. We show that under certain conditions the heat capacity of the three-state system can be almost indistinguishable from the heat capacity for the two-state system over a broad temperature interval. The consequence may be that some three-state (or even more than three-states) systems have been misinterpreted as two-state systems when the conclusion is drawn solely upon the van't Hoff enthalpy. These findings are important not only for proteins, but also for the interpretation of thermodynamical systems in general.     

Perchlorate-induced conformational transition of Staphylococcal nuclease: evidence for an equilibrium unfolding intermediate

The sodium perchlorate-induced conformational transition of Staphylococcal nuclease has been monitored by both circular dichroism (CD) and fluorescence spectroscopy. The perchlorate-induced transition is cooperative as observed by both spectroscopic signals. However, the protein loses only about one-third of its native far-UV CD signal at high perchlorate concentrations, indicating that a significant amount of secondary structure remains in the post-transition state. The remaining CD signal can be further diminished in a cooperative manner by the addition of the strong denaturant, urea. Near-UV CD spectra clearly show that the protein loses its tertiary structure in the perchlorate-induced denatured state. The perchlorate-induced transition curves were fit to the standard two-state model and the standard free energy change and m value of the transition are 2.3kcal/mol and 1.8kcal/(molM), respectively. By comparison, the urea-induced unfolding of Staphylococcal nuclease (in the absence of perchlorate) yields an unfolding free energy change, DeltaG(0,un), of 5.6kcal/mol and an m value of 2.3kcal/(molM). Thus, the thermodynamic state obtained in the post-transition region of perchlorate-induced conformation transition has a significantly lower free energy change, a high content of secondary structure, and diminished tertiary structure. These results suggest that the perchlorate-induced denatured state is a partially folded equilibrium state. Whether this intermediate is relevant to the folding/unfolding path under standard conditions is unknown at this time.     

Engineering stabilising beta-sheet interactions into a conformationally flexible region of the folding transition state of ubiquitin

Protein engineering studies suggest that the transition state for the folding of ubiquitin is highly polarised towards the N-terminal part of the sequence and involves a nucleus of residues within the beta-hairpin (residues 1-17) and main alpha-helix (residues 23-34). In contrast, the observation of small phi-values for residues in the C-terminal portion of the sequence (residues 35-76), coupled with a folding topology that results in a much higher contact order, suggests that fast folding of ubiquitin is dependent upon configurational flexibility in the C-terminal part of the polypeptide chain to ensure passage down a relatively smooth folding funnel to the native state. We show that the introduction of a small mini-hairpin motif as an extension of the native 43-50 hairpin stabilises local interactions in the C-terminal part of the sequence, resulting largely in a deceleration of the unfolding kinetics without perturbing the apparent two-state folding mechanism. However, a single-point Leu-->Phe substitution within the engineered hairpin sequence leads to the premature collapse of the denatured ensemble through the stabilisation of non-native interactions and the population of a compact intermediate. Non-linear effects in the kinetic data at low concentrations of denaturant suggest that the collapsed state, which is further stabilised in the presence of cosmotropic salts, may subsequently fold directly to the native state through a "triangular" reaction scheme involving internal rearrangement rather than unfolding and refolding.     

RNA tertiary interactions mediate native collapse of a bacterial group I ribozyme

Large RNAs collapse into compact intermediates in the presence of counterions before folding to the native state. We previously found that collapse of a bacterial group I ribozyme correlates with the formation of helices within the ribozyme core, but occurs at Mg2+ concentrations too low to support stable tertiary structure and catalytic activity. Here, using small-angle X-ray scattering, we show that Mg2+-induced collapse is a cooperative folding transition that can be fit by a two-state model. The Mg2+ dependence of collapse is similar to the Mg2+ dependence of helix assembly measured by partial ribonuclease T1 digestion and of an unfolding transition measured by UV hypochromicity. The correspondence between multiple probes of RNA structure further supports a two-state model. A mutation that disrupts tertiary contacts between the L9 tetraloop and its helical receptor destabilized the compact state by 0.8 kcal/mol, while mutations in the central triplex were less destabilizing. These results show that native tertiary interactions stabilize the compact folding intermediates under conditions in which the RNA backbone remains accessible to solvent.     

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.     

Titration Properties and Thermodynamics of the Transition State for Folding: Comparison of Two-state and Multi-state Folding Pathways

CI2 folds and unfolds as a single cooperative unit by simple two-state kinetics, which enables the properties of the transition state to be measured from both the forward and backward rate constants. We have examined how the free energy of the transition state for the folding of chymotrypsin inhibitor 2 (CI2) changes with pH and temperature. In addition to the standard thermodynamic quantities, we have measured the overall acid-titration properties of the transition state and its heat capacity relative to both the denatured and native states. We were able to determine the latter by a method analogous to a well-established procedure for measuring the change in heat capacity for equilibrium unfolding: the enthalpy of activation of unfolding at different values of acid pH were plotted against the average temperature of each determination. Our results show that the transition state of CI2 has lost most of the electrostatic and van der Waals' interactions that are found in the native state, but it remains compact and this prevents water molecules from entering some parts of the hydrophobic core. The properties of the transition state of CI2 are then compared with the major folding transition state of the larger protein barnase, which folds by a multi-state mechanism, with the accumulation of a partly structured intermediate (D^physor I). CI2 folds from a largely unstructured denatured state under physiological conditionsviaa transition state which is compact but relatively uniformly unstructured, with tertiary and secondary structure being formed in parallel. We term this an expanded pathway. Conversely, barnase folds from a largely structured denatured state in which elements of structure are well formed through a transition state that has islands of folded elements of structure. We term this a compact pathway. These two pathways may correspond to the two extreme ends of a continuous spectrum of protein folding mechanisms. Although the properties of the two transition states are very different, the activation barrier for folding (D^phys→3 ) is very similar for both proteins.     

Modeling protein density of states: additive hydrophobic effects are insufficient for calorimetric two-state cooperativity

A well-established experimental criterion for two-state thermodynamic cooperativity in protein folding is that the van't Hoff enthalpy DeltaH(vH) around the transition midpoint is equal, or very nearly so, to the calorimetric enthalpy DeltaH(cal) of the entire transition. This condition is satisfied by many small proteins. We use simple lattice models to provide a statistical mechanical framework to elucidate how this calorimetric two-state picture may be reconciled with the hierarchical multistate scenario emerging from recent hydrogen exchange experiments. We investigate the feasibility of using inverse Laplace transforms to recover the underlying density of states (i.e., enthalpy distribution) from calorimetric data. We find that the constraint imposed by DeltaH(vH)/DeltaH(cal) approximately 1 on densities of states of proteins is often more stringent than other "two-state" criteria proposed in recent theoretical studies. In conjunction with reasonable assumptions, the calorimetric two-state condition implies a narrow distribution of denatured-state enthalpies relative to the overall enthalpy difference between the native and the denatured conformations. This requirement does not always correlate with simple definitions of "sharpness" of a transition and has important ramifications for theoretical modeling. We find that protein models that assume capillarity cooperativity can exhibit overall calorimetric two-state-like behaviors. However, common heteropolymer models based on additive hydrophobic-like interactions, including highly specific two-dimensional Gō models, fail to produce proteinlike DeltaH(vH)/DeltaH(cal) approximately 1. A simple model is constructed to illustrate a proposed scenario in which physically plausible local and nonlocal cooperative terms, which mimic helical cooperativity and environment-dependent hydrogen bonding strength, can lead to thermodynamic behaviors closer to experiment. Our results suggest that proteinlike thermodynamic cooperativity may require a cooperative interplay between local and nonlocal interactions. The prospect of using calorimetric data to constrain Z-scores of knowledge-based potentials is discussed.     

Solvent-tuning the collapse and helix formation time scales of lambda(6-85)*

The lambda(6-85)(*) pseudo-wild type of lambda repressor fragment is a fast two-state folder (k(f) approximately 35 microsec(-1) at 58 degrees C). Previously, highly stable lambda(6-85)(*) mutants with k(f) > 30 microsec(-1) have been engineered to fold nearly or fully downhill. Stabilization of the native state by solvent tuning might also tune lambda(6-85)(*) away from two-state folding. We test this prediction by examining the folding thermodynamics and kinetics of lambda(6-85)(*) in a stabilizing solvent, 45% by weight aqueous ethylene glycol at -28 degrees C. Detection of kinetics by circular dichroism at 222 nm (sensitive to helix content) and small angle X-ray scattering (measuring the radius of gyration) shows that refolding from guanidine hydrochloride denatured conditions exhibits very different time scales for collapse and secondary structure formation: the two processes become decoupled. Collapse remains a low-barrier activated process, while the fastest of several secondary structure formation time scales approaches the downhill folding limit. Two-state folding of lambda(6-85)(*) is not a robust process.     

Microsecond Barrier-Limited Chain Collapse Observed by Time-Resolved FRET and SAXS

It is generally held that random-coil polypeptide chains undergo a barrier-less continuous collapse when the solvent conditions are changed to favor the fully folded native conformation. We test this hypothesis by probing intramolecular distance distributions during folding in one of the paradigms of folding reactions, that of cytochrome c. The Trp59-to-heme distance was probed by time-resolved Förster resonance energy transfer in the microsecond time range of refolding. Contrary to expectation, a state with a Trp59–heme distance close to that of the guanidinium hydrochloride (GdnHCl) denatured state is present after ~ 27 μs of folding. A concomitant decrease in the population of this state and an increase in the population of a compact high-FRET (Förster resonance energy transfer) state (efficiency > 90%) show that the collapse is barrier limited. Small-angle X-ray scattering (SAXS) measurements over a similar time range show that the radius of gyration under native favoring conditions is comparable to that of the GdnHCl denatured unfolded state. An independent comprehensive global thermodynamic analysis reveals that marginally stable partially folded structures are also present in the nominally unfolded GdnHCl denatured state. These observations suggest that specifically collapsed intermediate structures with low stability in rapid equilibrium with the unfolded state may contribute to the apparent chain contraction observed in previous fluorescence studies using steady-state detection. In the absence of significant dynamic averaging of marginally stable partially folded states and with the use of probes sensitive to distance distributions, barrier-limited chain contraction is observed upon transfer of the GdnHCl denatured state ensemble to native-like conditions.     

Kinetic analysis of local structure formations in protein folding

The protein folding process is described by a cluster model based on the assumption that local structures or clusters are formed at an early stage in different regions of the polypeptide chain. Possible local structural elements in a globular protein are helices, bends, and hydrophobic cores whose formation is presumably determined by the interaction with the environment. Thus the tendency of local structure formation is expressed by a surface free energy of the cluster, which is assigned to the interface between the cluster and its environment. The probability of finding the chain of N residues with k clusters and m residues in the cluster is represented by a cluster distribution map. The cluster model exhibits a distinct two-state-like equilibrium transition, which can be seen on this map as well-separated native and denatured populations at the midpoint of the transition. The native population is localized at k ≈ 1 and m ≈ N, while the position of the denatured population can vary significantly depending on the surface free energy of the cluster. If the surface free energy is strong, the denatured population is localized near k = 0 and m = 0. On the other hand, if the surface free energy is weak, the denatured population is localized at high k and m values. The dynamics of the cluster model are treated as a stochastic process involving the transition from a state (k,m) to one of its six neighbors. The transition probability for each transition is determined by the free energy difference between two states; thus no activation process is assumed. However, the conversion of the two macrostates, native and denatured populations, involves the free energy activation due to the cooperative interaction of the macrosystem. The dynamics are analyzed by following the time evolution of the population profile on the cluster distribution map. Kinetic schemes are proposed to describe the multistep mechanism of protein folding and unfolding.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Effects of excluded volume upon protein stability in covalently cross-linked proteins with variable linker lengths

To explore the effects of molecular crowding and excluded volume upon protein stability, we used a series of cross-linking reagents with nine different single-cysteine mutants of staphylococcal nuclease to make covalently linked dimers. These cross-linkers ranged in length from 10.5 to 21.3 A, compelling separations which would normally be found only in the most concentrated protein solutions. The stabilities of the dimeric proteins and monomeric controls were determined by guanidine hydrochloride and thermal denaturation. Dimers with short linkers tend to exhibit pronounced three-state denaturation behavior, as opposed to the two-state behavior of the monomeric controls. Increasing linker length leads to less pronounced three-state behavior. The three-state behavior is interpreted in a three-state model where cross-linked native protein dimer, N-N, interconverts in a two-state transition with a dimer where one protein subunit is denatured, N-D. The remaining native protein in turn can denature in another two-state transition to a state, D-D, in which both tethered proteins are denatured. Three-state behavior is best explained by excluded volume effects in the denatured state. For many dimers, linkers longer than 17 A removed most three-state character. This sets a limit on the flexibility and size of the denatured state. Notably, in contradiction to theoretical predictions, these cross-linked dimers were not stabilized. The failure of these predictions is possibly due to neglect of the alteration in hydrophobic exposure that accompanies any significant reduction in the conformational space of the denatured state.