Folding of S6 structures with divergent amino acid composition: pathway flexibility within partly overlapping foldons

Studies of circular permutants have demonstrated that the folding reaction of S6 from Thermus thermophilus (S6_T) is malleable and responds in an ordered manner to changes of the sequence separation between interacting residues: the S6_T permutants retain a common nucleation pattern in the form of a two-strand-helix motif that can be recruited from different parts of the structure. To further test the robustness of the two-strand-helix nucleus we have here determined the crystal structure and folding reaction of an evolutionary divergent S6 protein from the hyperthermophilic bacterium Aquifex aeolicus (S6_A). Although the overall topology of S6_A is very similar to that of S6_T the architecture of the hydrophobic core is radically different by containing a large proportion of stacked Phe side-chains. Despite this disparate core composition, the folding rate constant and the kinetic m values of S6_A are identical to those of S6_T. The folding nucleus of S6_A is also found to retain the characteristic two-strand-helix motif of the S6_T permutants, but with a new structural emphasis. The results suggest that the protein folding reaction is linked to topology only in the sense that the native-state topology determines the repertoire of accessible nucleation motifs. If the native structure allows several equivalent ways of recruiting a productive nucleus the folding reaction is free to redistribute within these topological constraints.     

Are turns required for the folding of ribonuclease T1?

Ribonuclease T1 (RNase T1) is a small, globular protein of 104 amino acids for which extensive thermodynamic and structural information is known. To assess the specific influence of variations in amino acid sequence on the mechanism for protein folding, circularly permuted variants of RNase T1 were constructed and characterized in terms of catalytic activity and thermodynamic stability. The disulfide bond connecting Cys-2 and Cys-10 was removed by mutation of these residues to alanine (C2, 10A) to avoid potential steric problems imposed by the circular permutations. The original amino-terminus and carboxyl-terminus of the mutant (C2, 10A) were subsequently joined with a tripeptide linker to accommodate a reverse turn and new termini were introduced throughout the primary sequence in regions of solvent-exposed loops at Ser-35 (cp35S1), Asp-49 (cp49D1), Gly-70 (cp70G1), and Ser-96 (cp96S1). These circularly permuted RNase T1 mutants retained 35-100% of the original catalytic activity for the hydrolysis of guanylyl(3'-->5')cytidine, suggesting that the overall tertiary fold of these mutants is very similar to that of wild-type protein. Chemical denaturation curves indicated thermodynamic stabilities at pH 5.0 of 5.7, 2.9, 2.6, and 4.6 kcal/mol for cp35S1, cp49D1, cp70G1, and cp96S1, respectively, compared to a value of 10.1 kcal/mol for wild-type RNase T1 and 6.4 kcal/mol for (C2, 10A) T1. A fifth set of circularly permuted variants was attempted with new termini positioned in a tight beta-turn between Glu-82 and Gln-85. New termini were inserted at Asn-83 (cp83N1), Asn-84 (cp84N1), and Gln-85 (cp85Q1). No detectable amount of protein was ever produced for any of the mutations in this region, suggesting that this turn may be critical for the proper folding and/or thermodynamic stability of RNase T1.     

Transition states for folding of circular-permuted proteins

To explore the role of entropy and chain connectivity in protein folding, a particularly interesting scheme, namely, the circular permutation, has been used. Recently, experimental observations showed that there are large differences in the folding mechanisms between the wild-type proteins and their circular permutants. These differences are strongly related to the change in the intrachain connectivity. Some results obtained by molecular dynamics simulations also showed a good agreement with the experimental findings. Here, we use a topology-based free-energy functional method to study the role of the chain connectivity in folding by comparing features of transition states of the wild-type proteins with those of their circular permutants. We concentrate our study on 3 small globular proteins, namely, the alpha-spectrin SH3 domain (SH3), the chymotrypsin inhibitor 2 (CI2), and the ribosomal protein S6, and obtain exciting results that are consistent with the available experimental and simulation results. A heterogeneity of the interaction energies between contacts for protein CI2 and for protein S6 is also introduced, which characterizes the strong interactions between contacts with long loops, as speculated from experiments for protein S6. The comparison between the folding nucleus of the wild-type proteins and those of their circular permutants indicates that chain connectivity affects remarkably the shapes of the energy profiles and thus the folding mechanism. Further comparisons between our theoretical calculated phi(th) values and the experimental observed phi(exp) values for the 3 proteins and their permutants show that our results are in good agreement with experimental ones and that correlations between them are high. These indicate that the free-energy functional method really provides a way to analyze the folding behavior of the circular-permuted proteins and therefore the folding mechanism of the wild-type proteins.     

Conformational plasticity in folding of the split beta-alpha-beta protein S6: evidence for burst-phase disruption of the native state

An increasing number of folding studies of two-state proteins shows that point mutations sometimes change the kinetic m-values, leading to kinks and curves in the chevron plots. The molecular origin of these changes is yet unclear although it is speculated that they are linked to structural rearrangement of the transition state or to accumulation of meta-stable intermediates. To shed more light on this issue, we present here a combined m and phi-value analysis of the split beta-alpha-beta protein S6. Wild-type S6 displays classical two-state kinetics with v-shaped chevron plot, but a majority of its mutants display distinct m-value changes or curved chevrons. We observe that this kinetic aberration of S6 is linked to mutations that are clustered in distinct regions of the native structure. The most pronounced changes, i.e. decrease in the m-value for the unfolding rate constant, are seen upon truncation of interactions between the N and C termini, whereas mutations in the centre of the hydrophobic core show smaller or even opposed effects. As a consequence, the calculated phi-values display a systematic increase upon addition of denaturant. In the case of S6, the phenomenon seems to arise from a general plasticity of the different species on the folding pathway. That is, the structure of the denatured ensemble, the transition state, and the native ground-state for unfolding seem to change upon mutation. From these changes, it is concluded that interactions spanning the centre of the hydrophobic core form early in folding, whereas the entropically disfavoured interactions linking the N and C termini consolidate very late, mainly on the down-hill-side of the folding barrier.     

Mapping the folding pathway of an immunoglobulin domain: structural detail from Phi value analysis and movement of the transition state

BACKGROUND: Do proteins that have the same structure fold by the same pathway even when they are unrelated in sequence? To address this question, we are comparing the folding of a number of different immunoglobulin-like proteins. Here, we present a detailed protein engineering phi value analysis of the folding pathway of TI I27, an immunoglobulin domain from human cardiac titin. RESULTS: TI I27 folds rapidly via a kinetic intermediate that is destabilized by most mutations. The transition state for folding is remarkably native-like in terms of solvent accessibility. We use phi value analysis to map this transition state and show that it is highly structured; only a few residues close to the N-terminal region of the protein remain completely unfolded. Interestingly, most mutations cause the transition state to become less native-like. This anti-Hammond behavior can be used as a novel means of obtaining additional structural information about the transition state. CONCLUSIONS: The residues that are involved in nucleating the folding of TI I27 are structurally equivalent to the residues that form the folding nucleus in an evolutionary unrelated fibronectin type III protein. These residues form part of the common structural core of Ig-like domains. The data support the hypothesis that interactions essential for defining the structure of these beta sandwich proteins are also important in nucleation of folding.     

Alteration of the disulfide-coupled folding pathway of BPTI by circular permutation

The kinetics of disulfide-coupled folding and unfolding of four circularly permuted forms of bovine pancreatic trypsin inhibitor (BPTI) were studied and compared with previously published results for both wild-type BPTI and a cyclized form. Each of the permuted proteins was found to be less stable than either the wild-type or circular proteins, by 3-8 kcal/mole. These stability differences were used to estimate effective concentrations of the chain termini in the native proteins, which were 1 mM for the wild-type protein and 2.5 to 4000 M for the permuted forms. The circular permutations increased the rates of unfolding and caused a variety of effects on the kinetics of refolding. For two of the proteins, the rates of a direct disulfide-formation pathway were dramatically increased, making this process as fast or faster than the competing disulfide rearrangement mechanism that predominates in the folding of the wild-type protein. These two permutations break the covalent connectivity among the beta-strands of the native protein, and removal of these constraints appears to facilitate direct formation and reduction of nearby disulfides that are buried in the folded structure. The effects on folding kinetics and mechanism do not appear to be correlated with relative contact order, a measure of overall topological complexity. These observations are consistent with the results of other recent experimental and computational studies suggesting that circular permutation may generally influence folding mechanisms by favoring or disfavoring specific interactions that promote alternative pathways, rather than through effects on the overall topology of the native protein.     

The folding of an enzyme. III. Structure of the transition state for unfolding of barnase analysed by a protein engineering procedure.

The structure of the first significant transition state on the unfolding pathway of barnase has been analysed in detail by protein engineering methods. Over 50 mutations placed strategically over the whole protein have been used as probes to report on the local structure in the transition state. Several different probes for many regions of the protein give consistent results as do multiple probes at the same site. The overall consistency of phi values indicates that the mutations have not produced changes in the protein that significantly alter the transition state for unfolding. A fine-structure analysis of interactions has also been conducted by removing different parts of the same side-chains. Many of the results of simple mutations fall nicely into the two clear-cut cases of phi = 1 or 0, indicating that the local noncovalent bonds are either fully broken or fully made in the transition state. Much of the structure of barnase in the transition state for unfolding is very similar to that in the folded protein. Both major alpha-helices fray at the N terminus. The last two turns in helix1 are certainly intact, as is the C terminus of helix2. The general picture of the beta-sheet is that the three central beta-strands are completely intact while the two edge beta-strands are mainly present but certainly weakened. The first five residues of the protein unwind but the C terminus remains folded. Three of the five loops are unfolded. The edges of the main hydrophobic core (core1) are significantly weakened, however, and their breaking appears partly rate determining. The centre of the small hydrophobic core3 remains intact. Core2 is completely disrupted. The first events in unfolding are thus: the unfolding of several loops, the unwinding of the helices from the N termini, and the weakening and disruption of the hydrophobic cores. The values of phi are found to be substantially the same under conditions that favour folding as under conditions that are highly denaturing, and so the structure of the unfolding transition state is substantially the same in water as in the presence of denaturant. The structure of the final kinetically significant transition state for refolding is identical to that for unfolding. The final events in refolding are, accordingly, the consolidation of the hydrophobic cores, the closing of many loops and the capping of the N termini of the helices.     

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.     

Sequence evolution and the mechanism of protein folding

The impact on protein evolution of the physical laws that govern folding remains obscure. Here, by analyzing in silico-evolved sequences subjected to evolutionary pressure for fast folding, it is shown that: First, a subset of residues in the thermodynamic folding nucleus is mainly responsible for modulating the protein folding rate. Second and most important, the protein topology itself is of paramount importance in determining the location of these residues in the structure. Further stabilization of the interactions in this nucleus leads to fast folding sequences. Third, these nucleation points restrict the sequence space available to the protein during evolution. Correlated mutations between positions around these hot spots arise in a statistically significant manner, and most involve contacting residues. When a similar analysis is carried out on real proteins, qualitatively similar results are obtained.     

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3?kcal?mol^?1. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.     

Protein β-Sheet Nucleation Is Driven by Local Modular Formation

Despite its central role in the protein folding process, the specific mechanism(s) behind β-sheet formation has yet to be determined. For example, whether the nucleation of β-sheets, often containing strands separated in sequence by many residues, is local or not remains hotly debated. Here, we investigate the initial nucleation step of β-sheet formation by performing an analysis of the smallest β-sheets in a non-redundant dataset on the grounds that the smallest sheets, having undergone little growth after nucleation, will be enriched for nucleating characteristics. We find that the residue propensities are similar for small and large β-sheets as are their interstrand pairing preferences, suggesting that nucleation is not primarily driven by specific residues or interacting pairs. Instead, an examination of the structural environments of the two-stranded sheets shows that virtually all of them are contained in single, compact structural modules, or when multiple modules are present, one or both of the chain termini are involved. We, therefore, find that β-nucleation is a local phenomenon resulting either from sequential or topological proximity. We propose that β-nucleation is a result of two opposite factors; that is, the relative rigidity of an associated folding module that holds two stretches of coil close together in topology coupled with sufficient chain flexibility that enables the stretches of coil to bring their backbones in close proximity. Our findings lend support to the hydrophobic zipper model of protein folding (Dill, K. A., Fiebig, K. M., and Chan, H. S. (1993) Proc. Natl. Acad. Sci. U.S.A. 90, 1942–1946). Implications for protein folding are discussed.     

Changes of protein folding pathways by circular permutation. Overlapping nuclei promote global cooperativity

The evolved properties of proteins are not limited to structure and stability but also include their propensity to undergo local conformational changes. The latter, dynamic property is related to structural cooperativity and is controlled by the folding-energy landscape. Here we demonstrate that the structural cooperativity of the ribosomal protein S6 is optimized by geometric overlap of two competing folding nuclei: they both include the central beta-strand 1. In this way, folding of one nucleus catalyzes the formation of the other, contributing to make the folding transition more concerted overall. The experimental evidence is provided by an extended set of circular permutations of S6 that allows quantitative analysis of pathway plasticity at the level of individual side chains. Because similar overlap between competing nuclei also has been discerned in other proteins, we hypothesize that the coupling of several small nuclei into extended "supernuclei" represents a general principle for propagating folding cooperativity across large structural distances.     

Fragment reconstitution of a small protein: disulfide mutant of a short C-terminal fragment derived from streptococcal protein G

Hierarchical studies on the folding of protein G B1 domain have shown that the C-terminal fragment (C16) has a considerable amount of beta-hairpin structure that exchanges between the folded and unfolded states at room temperature, and that the C16 fragment binds noncovalently to an N-terminal fragment (N40) under physiological conditions. Those studies have led us to the hypothesis that the amphipathic beta-hairpin structure of C16 initiates folding of the domain. To obtain a more detailed understanding of the folding mechanism of the domain, we designed a mutant of C16 (SS16ox) with a disulfide bond between residues 41 and 56, and then examined the interaction of the mutant with N40 by surface plasmon resonance (SPR) and by thermal denaturation studies using circular dichroism. SS16ox strongly interacted with N40, with an equilibrium constant, KD, that was 7-fold higher than wild-type. The association rate constant, kon, of SS16ox was 8.7-fold higher than that of wild-type. This strong interaction can be explained by the entropic effect of the disulfide bond. The introduction of the disulfide bond into C16 stabilizes the beta-hairpin structure of C16, accelerates the association rate with N40, and then stabilizes the whole complex. These results support a hypothetical folding mechanism of protein G where the amphipathic beta-hairpin structure of C16 acts as a nucleus and accelerates folding of the whole molecule.     

Comparison of local and global stability of an analogue of a disulfide-folding intermediate with those of the wild-type protein in bovine pancreatic ribonuclease A: identification of specific regions of stable structure along the oxidative folding pathway

We have identified specific regions of the polypeptide chain of bovine pancreatic ribonuclease A (RNase A) that are critical for stabilizing the oxidative folding intermediate des-[40-95] (with three native disulfide bonds but lacking the fourth native Cys40-Cys95 disulfide bond) in an ensemble of largely disordered three-disulfide precursors (3S if des-[40-95]). A stable analogue of des-[40-95], viz., [C40A, C95A] RNase A, which contains three out of four native disulfide pairings, was previously found to have a three-dimensional structure very similar to that of the wild-type protein. However, it is determined here from GdnHCl denaturation experiments to have significantly reduced global stability, i.e., = 4.5 kcal /mol at 20 degrees C and pH 4.6. The local stability of [C40A, C95A] RNase A was also examined using site-specific amide (2)H/(1)H exchange measurements at pD 5.0 to determine the individual unfolding free energy of specific residues under both strongly native (12 degrees C) and more destabilizing (20 degrees C) conditions. Comparison of the relative stabilities at specific amide sites of [C40A, C95A] RNase A at both temperatures with the corresponding values for the wild-type protein at 35 degrees C corroborates previous experimental evidence that unidentified intramolecular contacts in the vicinity of the preferentially formed native one-disulfide (C65-C72) loop are crucial for stabilizing early folding intermediates, leading to des-[40-95]. Moreover, values of for residues at or near the third alpha-helix, and in part of the second beta-sheet of [C40A, C95A] RNase A, indicate that these two regions of regular backbone structure contribute to stabilizing the global chain fold of the des-[40-95] disulfide-folding intermediate in the wild-type protein. More significantly, we have identified numerous specific residues in the first alpha-helix and the first beta-sheet of the protein that are stabilized in the final step of the major oxidative regeneration pathway of RNase A (des-[40-95] --> N).     

Folding Intermediate and Folding Nucleus for I→N and U→I→N Transitions in Apomyoglobin: Contributions by Conserved and Nonconserved Residues

Kinetic investigation on the wild-type apomyoglobin and its 12 mutants with substitutions of hydrophobic residues by Ala was performed using stopped-flow fluorescence. Characteristics of the kinetic intermediate I and the folding nucleus were derived solely from kinetic data, namely, the slow-phase folding rate constants and the burst-phase amplitudes of Trp fluorescence intensity. This allowed us to pioneer the ?-analysis for apomyoglobin. As shown, these mutations drastically destabilized the native state N and produced minor (for conserved residues of G, H helices) or even negligible (for nonconserved residues of B, C, D, E helices) destabilizing effect on the state I. On the other hand, conserved residues of A, G, H helices made a smaller contribution to stability of the folding nucleus at the rate-limiting I→N transition than nonconserved residues of B, D, E helices. Thus, conserved side chains of the A-, G-, H-residues become involved in the folding nucleus before crossing the main barrier, whereas nonconserved side chains of the B-, D-, E-residues join the nucleus in the course of the I→N transition.     

Different circular permutations produced different folding nuclei in proteins: a computational study

There have been many studies about the effect of circular permutation on the transition state/folding nucleus of proteins, with sometimes conflicting conclusions from different proteins and permutations. To clarify this important issue, we have studied two circular permutations of a lattice protein model with side-chains. Both permuted sequences have essentially the same native state as the original (wild-type) sequence. Circular permutant 1 cuts at the folding nucleus of the wild-type sequence. As a result, the permutant has a drastically different nucleus and folds more slowly than wild-type. In contrast, circular permutant 2 involves an incision at a site unstructured in the wild-type transition state, and the wild-type nucleus is largely retained in the permutant. In addition, permutant 2 displays both two-state and multi-state folding, with a native-like intermediate state occasionally populated. Neither the wild-type nor permutant 1 has a similar intermediate, and both fold in an apparently two-state manner. Surprisingly, permutant 2 folds at a rate identical with that of the wild-type. The intermediate in permutant 2 is stabilised by native and non-native interactions, and cannot be classified simply as on or off-pathway. So we advise caution in attributing experimental data to on or off-pathway intermediates. Finally, our work illuminates the results on alpha-spectrin SH3, chymotrypsin inhibitor 2 and beta-lactoglobulin, and supports a key assumption in the experimental efforts to locate potential nucleation sites of real proteins via circular permutations.     

Random circular permutation of DsbA reveals segments that are essential for protein folding and stability

One of the key questions in protein folding is whether polypeptide chains require unique nucleation sites to fold to the native state. In order to identify possible essential polypeptide segments for folding, we have performed a complete circular permutation analysis of a protein in which the natural termini are in close proximity. As a model system, we used the disulfide oxidoreductase DsbA from Escherichia coli, a monomeric protein of 189 amino acid residues. To introduce new termini at all possible positions in its polypeptide chain, we generated a library of randomly circularly permuted dsbA genes and screened for active circularly permuted variants in vivo. A total of 51 different active variants were identified. The new termini were distributed over about 70 % of the polypeptide chain, with the majority of them occurring within regular secondary structures. New termini were not found in approximately 30 % of the DsbA sequence which essentially correspond to four alpha-helices of DsbA. Introduction of new termini into these "forbidden segments" by directed mutagenesis yielded proteins with altered overall folds and strongly reduced catalytic activities. In contrast, all active variants analysed so far show structural and catalytic properties comparable with those of DsbA wild-type. We suggest that random circular permutation allows identification of contiguous structural elements in a protein that are essential for folding and stability.     

Identification of human CYP2C19 residues that confer S-mephenytoin 4'-hydroxylation activity to CYP2C9

CYP2C19 is selective for the 4'-hydroxylation of S-mephenytoin while the highly similar CYP2C9 has little activity toward this substrate. To identify critical amino acids determining the specificity of human CYP2C19 for S-mephenytoin 4'-hydroxylation, we constructed chimeras by replacing portions of CYP2C9 containing various proposed substrate recognition sites (SRSs) with those of CYP2C19 and mutating individual residues by site-directed mutagenesis. Only a chimera containing regions encompassing SRSs 1--4 was active (30% of wild-type CYP2C19), indicating that multiple regions are necessary to confer specificity for S-mephenytoin. Mutagenesis studies identified six residues in three topological components of the proteins required to convert CYP2C9 to an S-mephenytoin 4'-hydroxylase (6% of the activity of wild-type CYP2C19). Of these, only the I99H difference located in SRS 1 between helices B and C reflects a change in a side chain that is predicted to be in the substrate-binding cavity formed above the heme prosthetic group. Two additional substitutions, S220P and P221T residing between helices F and G but not in close proximity to the substrate binding site together with five differences in the N-terminal portion of helix I conferred S-mephenytoin 4'-hydroxylation activity with a K(M) similar to that of CYP2C19 but a 3-fold lower K(cat). Three residues in helix I, S286N, V292A, and F295L, were essential for S-mephenytoin 4'-hydroxylation activity. On the basis of the structure of the closely related enzyme CYP2C5, these residues are unlikely to directly contact the substrate during catalysis but are positioned to influence the packing of substrate binding site residues and likely substrate access channels in the enzyme.     

Are the same or different amino acid residues responsible for correct and incorrect protein folding?

It has been shown for 20 proteins that amino acid residues included into the protein folding nucleus, determined experimentally, are often involved in the theoretically determined amyloidogenic fragments. For 18 proteins, Φ-values indicative of the extent of residue involvement into the folding nucleus are on average higher for amino acid residues within amyloidogenic regions. Amyloidogenic fragments were predicted for 20 proteins by two methods chosen from four on the basis of comparison of prediction of amyloidogenic regions known from experimental data. Since theoretical folding nuclei are detected by the protein three-dimensional structure and amyloidogenic regions by the protein chain primary structure, the detected regularity makes possible predictions of folding nucleation sites on the basis of amino acid sequence.     

Time-resolved backbone desolvation and mutational hot spots in folding proteins

A method is presented to identify hot mutational spots and predict the extent of surface burial at the transition state relative to the native fold in two-state folding proteins. The method is based on ab initio simulations of folding histories in which transitions between coarsely defined conformations and pairwise interactions are dependent on the solvent environments created by the chain. The highly conserved mammalian ubiquitin is adopted as a study case to make predictions. The evolution in time of the chain topology suggests a nucleation process with a critical point signaled by a sudden quenching of structural fluctuations. The occurrence of this nucleus is shown to be concurrent with a sudden escalation in the number of three-body correlations whereby hydrophobic units approach residue pairs engaged in amide-carbonyl hydrogen bonding. These correlations determine a pattern designed to structure the surrounding solvent, protecting intramolecular hydrogen bonds from water attack. Such correlations are shown to be required to stabilize the nucleus, with kinetic consequences for the folding process. Those nuclear residues that adopt the dual role of protecting and being protected while engaged in hydrogen bonds are predicted to be the hottest mutational spots. Some such residues are shown not to retain the same protecting role in the native fold. This kinetic treatment of folding nucleation is independently validated vis-a-vis a Phi-value analysis on chymotrypsin inhibitor 2, a protein for which extensive mutational data exists.