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.     

The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H

Elucidation of the high-resolution structures of folding intermediates is a necessary but difficult step toward the ultimate understanding of the mechanism of protein folding. Here, using hydrogen-exchange-directed protein engineering, we populated the folding intermediate of the Thermus thermophilus ribonuclease H, which forms before the rate-limiting transition state, by removing the unfolded regions of the intermediate, including an α-helix and two β-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 Å). It has a native-like backbone topology and shows some local deviations from the native structure, revealing that the structure of the folded region of an early folding intermediate can be as well defined as the native structure. The topological parameters calculated from the structures of the intermediate mimic and the native state predict that the intermediate should fold on a millisecond time scale or less and form much faster than the native state. Other factors that may lead to the slow folding of the native state and the accumulation of the intermediate before the rate-limiting transition state are also discussed.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

Early Closure of a Long Loop in the Refolding of Adenylate Kinase: A Possible Key Role of Non-Local Interactions in the Initial Folding Steps

Most globular protein chains, when transferred from high to low denaturant concentrations, collapse instantly before they refold to their native state. The initial compaction of the protein molecule is assumed to have a key effect on the folding pathway, but it is not known whether the earliest structures formed during or instantly after collapse are defined by local or by non-local interactions—that is, by secondary structural elements or by loop closure of long segments of the protein chain. Stable closure of one or several long loops can reduce the chain entropy at a very early stage and can prevent the protein from following non-productive pathways whose number grows exponentially with the length of the protein chain. In Escherichia coli adenylate kinase (AK), about seven long loops define the topology of the native structure. We selected four loop-forming sections of the chain and probed the time course of loop formation during refolding of AK. We labeled the termini of the loop segments with tryptophan and cysteine-5-amidosalicylic acid. This donor-acceptor pair of probes used with fluorescence resonance excitation energy transfer spectroscopy (FRET) is suitable for detecting very short distances and thus is able to distinguish between random and specific compactions. Refolding of AK was initiated by stopped-flow mixing, followed simultaneously by donor and acceptor fluorescence, and analyzed in terms of energy transfer efficiency and distance. In the collapsed state of AK, observed after the 5-ms dead time of the instrument, one of the selected segments shows a native-like separation of its termini; it forms a loop already in the collapsed state. A second segment that includes the first but is longer by 15 residues shows an almost native-like separation of its termini. In contrast, a segment that is shorter but part of the second segment shows a distance separation of its termini as high as a segment that spans almost the whole protein chain. We conclude that a specific network of non-local interactions, the closure of one or several loops, can play an important role in determining the protein folding pathway at its early phases.     

The cold denatured state is compact but expands at low temperatures: hydrodynamic properties of the cold denatured state of the C-terminal domain of L9

A point mutation of a small globular protein, the C-terminal domain of L9 destabilizes the protein and leads to observable cold-denaturation at temperatures above zero. The cold denatured state is in slow exchange with the native state on the NMR time scale, and this allows the hydrodynamic properties of the cold unfolded state and the native state to be measured under identical conditions using pulsed-field gradient NMR diffusion measurements. This provides the first experimental measurement of the hydrodynamic properties of a cold unfolded protein and its folded form under identical conditions. Hydrodynamic radii of the cold-induced unfolded states were measured for a set of temperatures ranging from 2 °C to 25 °C at pD 6.6 in the absence of denaturant. The cold unfolded state is compact compared to the urea or acid unfolded state and a trend of increasing radii of hydration is observed as the temperature is lowered. These observations are confirmed by experiments on the same protein at pD 8.0, where it is more stable, in the presence of a modest concentration of urea. The expansion of the cold-denatured state at lower temperatures is consistent with the temperature dependence of hydrophobic interactions.     

Population of on-pathway intermediates in the folding of ubiquitin

The role that intermediate states play in protein folding is the subject of intense investigation and in the case of ubiquitin has been controversial. We present fluorescence-detected kinetic data derived from single and double mixing stopped-flow experiments to show that the F45W mutant of ubiquitin (WT*), a well-studied single-domain protein and most recently regarded as a simple two-state system, folds via on-pathway intermediates. To account for the discrepancy we observe between equilibrium and kinetic stabilities and m-values, we show that the polypeptide chain undergoes rapid collapse to an intermediate whose presence we infer from a fast lag phase in interrupted refolding experiments. Double-jump kinetic experiments identify two direct folding phases that are not associated with slow isomerisation reactions in the unfolded state. These two phases are explained by kinetic partitioning which allows molecules to reach the native state from the collapsed state via two possible competing routes, which we further examine using two destabilised ubiquitin mutants. Interrupted refolding experiments allow us to observe the formation and decay of an intermediate along one of these pathways. A plausible model for the folding pathway of ubiquitin is presented that demonstrates that obligatory intermediates and/or chain collapse are important events in restricting the conformational search for the native state of ubiquitin.     

Coulomb forces control the density of the collapsed unfolded state of barstar

Although it has been recently shown that unfolded polypeptide chains undergo a collapse on transfer from denaturing to native conditions, the forces determining the dynamics and the size of the collapsed form have not yet been understood. Here, we use single-molecule fluorescence resonance energy transfer experiments on the small protein barstar to characterize the unfolded chain in guanidinium chloride (GdmCl) and urea. The unfolded protein collapses on decreasing the concentration of denaturants. Below the critical concentration of 3.5 M denaturant, the collapse in GdmCl leads to a more dense state than in urea. Since it is known that GdmCl suppresses electrostatic interactions, we infer that Coulomb forces are the dominant forces acting in the unfolded barstar under native conditions. This hypothesis is clearly buttressed by the finding of a compaction of the unfolded barstar by addition of KCl at low urea concentrations.     

Characterization of transient intermediates in lysozyme folding with time-resolved small-angle X-ray scattering.

We have used synchrotron radiation, together with stopped-flow and continuous-flow mixing techniques to monitor refolding of lysozyme at pH 5.2. From data measured at times which range from 14 ms to two seconds, we can monitor changes in the size, the shape and the pair distribution function of the polypeptide chain during the folding process. Comparison of the results with the properties of native and GdmCl-unfolded lysozyme shows that a major chain collapse occurs in the dead-time of mixing. During this process about 50 % of the change in radius of gyration between the unfolded protein and the native state occurs and the polypeptide chain adopts a globular shape. Time-resolved fluorescence spectra of this collapsed state suggest that the hydrophobic side-chains are still highly solvent accessible. A subsequently formed intermediate with helical structure in the alpha-domain is nearly identical in size and shape with native lysozyme and has a solvent-inaccessible hydrophobic core. Despite its native-like properties, this intermediate is only slightly more stable (DeltaG0=-4 kJ/mol) than the collapsed state and still much less stable than native lysozyme (DeltaDeltaG0=36 kJ/mol) at 20 degrees C.     

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.     

Molecular Determinants of a Native-State Prolyl Isomerization

Prolyl cis/trans isomerizations determine the rates of many protein-folding reactions, and they can serve as molecular switches and timers. The energy required to shift the prolyl cis/trans equilibrium during these processes originates from conformational reactions that are linked structurally and energetically with prolyl isomerization. We used the N2 domain of the gene-3-protein of phage fd to elucidate how such an energetic linkage develops in the course of folding. The Asp160-Pro161 bond at the tip of a β hairpin of N2 is cis in the crystal structure, but in fact, it exists as a mixture of conformers in folded N2. During refolding, about 10 kJ mol^− 1 of conformational energy becomes available for a 75-fold shift of the cis/trans equilibrium constant at Pro161, from 7/93 in the unfolded to 90/10 in the folded form. We combined single- and double-mixing kinetic experiments with a mutational analysis to identify the structural origin of this proline shift energy and to elucidate the molecular path for the transfer of this energy to Pro161. It originates largely, if not entirely, from the two-stranded β sheet at the base of the Pro161 hairpin. The two strands improve their stabilizing interactions when Pro161 is cis, and this stabilization is propagated to Pro161, because the connector peptides between the β strands and Pro161 are native-like folded when Pro161 is cis. In the presence of a trans-Pro161, the connector peptides are locally unfolded, and thus, Pro161 is structurally and energetically uncoupled from the β sheet. Such interrelations between local folding and prolyl isomerization and the potential modulation by prolyl isomerases might also be used to break and reestablish slow communication pathways in proteins.     

Microsecond hydrophobic collapse in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein

Using small-angle X-ray scattering combined with a continuous-flow mixing device, we monitored the microsecond compaction dynamics in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein. A significant collapse of the radius of gyration from 30 A to 23.2 A occurs within 300 micros after the initiation of refolding by a urea dilution jump. The subsequent folding after the major chain collapse occurs on a considerably longer time-scale. The protein folding trajectories constructed by comparing the development of the compactness and the secondary structure suggest that the specific hydrophobic collapse model rather than the framework model better explains the experimental observations. The folding trajectory of this alpha/beta-type protein is located between those of alpha-helical and beta-sheet proteins, suggesting that native structure determines the folding landscape.     

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.     

Evidence for Initial Non-specific Polypeptide Chain Collapse During the Refolding of the SH3 Domain of PI3 Kinase

Refolding of the SH3 domain of PI3 kinase from the guanidine hydrochloride (GdnHCl)-unfolded state has been probed with millisecond (stopped flow) and sub-millisecond (continuous flow) measurements of the change in fluorescence, circular dichroism, ANS fluorescence and three-site fluorescence resonance energy transfer (FRET) efficiency. Fluorescence measurements are unable to detect structural changes preceding the rate-limiting step of folding, whereas measurements of changes in ANS fluorescence and FRET efficiency indicate that polypeptide chain collapse precedes the major structural transition. The initial chain collapse reaction is complete within 150 μs. The collapsed form at this time possesses hydrophobic clusters to which ANS binds. Each of the three measured intra-molecular distances has contracted to an extent predicted by the dependence of the FRET signal in completely unfolded protein on denaturant concentration, indicating that contraction is non-specific. The extent of contraction of each intra-molecular distance in the collapsed product of sub-millisecond folding increases continuously with a decrease in [GdnHCl]. The gradual contraction is continuous with the gradual contraction seen in completely unfolded protein, and its dependence on [GdnHCl] is not indicative of an all-or-none collapse reaction. The dependence of the extent of contraction on [GdnHCl] was similar for the three distances, indicating that chain collapse occurs in a synchronous manner across different segments of the polypeptide chain. The sub-millisecond measurements of folding in GdnHCl were unable to determine whether hydrophobic cluster formation, probed by ANS fluorescence measurement, precedes chain contraction probed by FRET. To determine whether hydrogen bonding plays a role in initial chain collapse, folding was initiated by dilution of the urea-unfolded state. The extent of contraction of at least one intra-molecular distance in the collapsed product of sub-millisecond folding in urea is similar to that seen in GdnHCl, and the initial contraction in urea too appears to be gradual.     

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.     

The unfolded state of the villin headpiece helical subdomain: computational studies of the role of locally stabilized structure

The 36 residue villin headpiece helical subdomain (HP36) is one of the fastest cooperatively folding proteins, folding on the microsecond timescale. HP36's simple three helix topology, fast folding and small size have made it an attractive model system for computational and experimental studies of protein folding. Recent experimental studies have explored the denatured state of HP36 using fragment analysis coupled with relatively low-resolution spectroscopic techniques. These studies have shown that there is apparently only a small tendency to form locally stabilized secondary structure. Here, we complement the experimental studies by using replica exchange molecular dynamics with explicit solvent to investigate the structural features of these peptide models of unfolded HP36. To ensure convergence, two sets of simulations for each fragment were performed with different initial structures, and simulations were continued until these generated very similar final ensembles. These simulations reveal low populations of native-like structure and early folding events that cannot be resolved by experiment. For each fragment, calculated J-coupling constants and helical propensities are in good agreement with experimental trends. HP-1, corresponding to residues 41 to 53 and including the first alpha-helix, contains the highest helical population. HP-3, corresponding to residues 62 through 75 and including the third alpha-helix, contains a small population of helical turn residing at the N terminus while HP-2, corresponding to residues 52 through 61 and including the second alpha-helix, formed little to no structure in isolation. Overall, HP-1 was the only fragment to adopt a native-like conformation, but the low population suggests that formation of significant structure only occurs after formation of specific tertiary interactions.     

Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

Influence of the internal disulfide bridge on the folding pathway of the CL antibody domain

Disulfide bridges are one of the most important factors stabilizing the native structure of a protein. Whereas the basis for their stabilizing effect is well understood, their role in a protein folding reaction still seems to require further attention. We used the constant domain of the antibody light chain (C(L)), a representative of the ubiquitous immunoglobulin (Ig)-superfamily, to delineate the kinetic role of its single buried disulfide bridge. Independent of its redox state, the monomeric C(L) domain adopts a typical Ig-fold under native conditions and does not retain significant structural elements when unfolded. Interestingly, its folding pathway is strongly influenced by the disulfide bridge. The more stable oxidized protein folds via a highly structured on-pathway intermediate, whereas the destabilized reduced protein populates a misfolded off-pathway species on its way to the native state. In both cases, the formation of the intermediate species is shown to be independent of the isomerization state of the Tyr(141)-Pro(142) bond. Our results demonstrate that the internal disulfide bridge in an antibody domain restricts the folding pathway by bringing residues of the folding nucleus into proximity thus facilitating the way to the native state.     

Topology and sequence in the folding of a TIM barrel protein: global analysis highlights partitioning between transient off-pathway and stable on-pathway folding intermediates in the complex folding mechanism of a (betaalpha)8 barrel of unknown function from B. subtilis

The relative contributions of chain topology and amino acid sequence in directing the folding of a (betaalpha)(8) TIM barrel protein of unknown function encoded by the Bacillus subtilis iolI gene (IOLI) were assessed by reversible urea denaturation and a combination of circular dichroism, fluorescence and time-resolved fluorescence anisotropy spectroscopy. The equilibrium reaction for IOLI involves, in addition to the native and unfolded species, a stable intermediate with significant secondary structure and stability and self-associated forms of both the native and intermediate states. Global kinetic analysis revealed that the unfolded state partitions between an off-pathway refolding intermediate and the on-pathway equilibrium intermediate early in folding. Comparisons with the folding mechanisms of two other TIM barrel proteins, indole-3-glycerol phosphate synthase from the thermophile Sulfolobus solfataricus (sIGPS) and the alpha subunit of Escherichia coli tryptophan synthase (alphaTS), reveal striking similarities that argue for a dominant role of the topology in both early and late events in folding. Sequence-specific effects are apparent in the magnitudes of the relaxation times and relative stabilities, in the presence of additional monomeric folding intermediates for alphaTS and sIGPS and in rate-limiting proline isomerization reactions for alphaTS.     

The denatured state of N-PGK is compact and predominantly disordered

The organisation of the structure present in the chemically denatured N-terminal domain of phosphoglycerate kinase (N-PGK) has been determined by paramagnetic relaxation enhancements (PREs) to define the conformational landscape accessible to the domain. Below 2.0 M guanidine hydrochloride (GuHCl), a species of N-PGK (denoted I_b) is detected, distinct from those previously characterised by kinetic experiments [folded (F), kinetic intermediate (I_k) and denatured (D)]. The transition to I_b is never completed at equilibrium, because F predominates below 1.0 M GuHCl. Therefore, the ability of PREs to report on transient or low population species has been exploited to characterise I_b. Five single cysteine variants of N-PGK were labelled with the nitroxide electron spin-label MTSL [(1-oxyl-2,2,5,5-tetramethyl-3-pyrroline-3-methyl)methanesulfonate] and the denaturant dependences of the relaxation properties of the amide NMR signals between 1.2 and 3.6 M GuHCl were determined. Significant PREs for I_b were obtained, but these were distributed almost uniformly throughout the sequence. Furthermore, the PREs indicate that no specific short tertiary contacts persist. The data indicate a collapsed state with no coherent three-dimensional structure, but with a restricted radius beyond which the protein chain rarely reaches. The NMR characteristics of I_b indicate that it forms from the fully denatured state within 100 μs, and therefore a rapid collapse is the initial stage of folding of N-PGK from its chemically denatured state. By extrapolation, I_b is the predominant form of the denatured state under native conditions, and the non-specifically collapsed structure implies that many non-native contacts and chain reversals form early in protein folding and must be broken prior to attaining the native state topology.     

Characterisation of the conformational properties of urea-unfolded Im7: implications for the early stages of protein folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.