Study of the compact denatured state of a protein by molecular dynamics simulation]

The native state can be considered as a unique conformation of the protein molecule with the lowest free energy of residue contacts. In this case, all other conformations correspond to the denatured state. The degree of their compactness varies significantly. Under folding conditions, the compact denatured state rather than the random coil is in equilibrium with native protein. The balance between the main forces of protein folding, the solvophobic interactions and conformational entropy, suggests that some properties of the compact denatured state are close to those of native protein, whereas other properties are close to those of the random coil. To investigate the molecular structure of the compact denatured state, the method of molecular dynamics simulation seems to be very useful.     

Tryptophan fluorescence reveals the presence of long-range interactions in the denatured state of ribonuclease Sa

Characterizing the denatured state ensemble is crucial to understanding protein stability and the mechanism of protein folding. The aim of this research was to see if fluorescence could be used to gain new information on the denatured state ensemble. Ribonuclease Sa (RNase Sa) contains no Trp residues. We made five variants of RNase Sa by adding Trp residues at locations where they are found in other members of the microbial ribonuclease family. To better understand the protein denatured state, we also studied the fluorescence properties of the following peptides: N-acetyl-Trp-amide (NATA), N-acetyl-Ala-Trp-Ala-amide (AWA), N-acetyl-Ala-Ala-Trp-Ala-Ala-amide (AAWAA), and the five pentapeptides with the same sequence as the Trp substitution sites in RNase Sa. The major conclusions are: 1), the wavelength of maximum fluorescence intensity, λ_max, does not differ significantly for the peptides and the denatured proteins; 2), the fluorescence intensity at λ_max, I_F, differs significantly for the five Trp containing variants of RNase Sa; 3), the I_F differences for the denatured proteins are mirrored in the peptides, showing that the short-range effects giving rise to the I_F differences in the peptides are also present in the proteins; 4) the I_F values for the denatured proteins are more than 30% greater than for the peptides, showing the presence of long-range effects in the proteins; 5), fluorescence quenching of Trp by acrylamide and iodide is more than 50% greater in the peptides than in the denatured proteins, showing that long-range effects limit the accessibility of the quenchers to the Trp side chains in the proteins; and 6), these results show that nonlocal effects in the denatured states of proteins influence Trp fluorescence and accessibility significantly.     

The denatured state of Engrailed Homeodomain under denaturing and native conditions

Protein folding starts from the elusive form of the denatured state that is present under conditions that favour the native state. We have studied the denatured state of Engrailed Homeodomain (En-HD) under mildly and strongly denaturing conditions at the level of individual residues by NMR and more globally by conventional spectroscopy and solution X-ray scattering. We have compared these states with a destabilized mutant, L16A, which is predominantly denatured under conditions where the wild-type is native. This engineered denatured state, which could be directly studied under native conditions, was in genuine equilibrium with the native state, which could be observably populated by changing the conditions or introducing a stabilizing mutation. The denatured state had extensive native secondary structure and was significantly compact and globular. But, the side-chains and backbone were highly mobile. Non-cooperative melting of the residual structure on the denatured state of En-HD was observed, both at the residue and the molecular level, with increasingly denaturing conditions. The absence of a co-operative transition could result from the denatured state ensemble progressing through a series of intermediates or from a more general slide (second-order transition) from the compact form under native conditions to the more extended at highly denaturing conditions. In either case, the starting point for folding under native conditions is highly structured and already poised to adopt the native structure.     

Uncovering Specific Electrostatic Interactions in the Denatured States of Proteins

The stability and folding of proteins are modulated by energetically significant interactions in the denatured state that is in equilibrium with the native state. These interactions remain largely invisible to current experimental techniques, however, due to the sparse population and conformational heterogeneity of the denatured-state ensemble under folding conditions. Molecular dynamics simulations using physics-based force fields can in principle offer atomistic details of the denatured state. However, practical applications are plagued with the lack of rigorous means to validate microscopic information and deficiencies in force fields and solvent models. This study presents a method based on coupled titration and molecular dynamics sampling of the denatured state starting from the extended sequence under native conditions. The resulting denatured-state pK_as allow for the prediction of experimental observables such as pH- and mutation-induced stability changes. I show the capability and use of the method by investigating the electrostatic interactions in the denatured states of wild-type and K12M mutant of NTL9 protein. This study shows that the major errors in electrostatics can be identified by validating the titration properties of the fragment peptides derived from the sequence of the intact protein. Consistent with experimental evidence, our simulations show a significantly depressed pK_a for Asp⁸ in the denatured state of wild-type, which is due to a nonnative interaction between Asp⁸ and Lys¹². Interestingly, the simulation also shows a nonnative interaction between Asp⁸ and Glu⁴⁸ in the denatured state of the mutant. I believe the presented method is general and can be applied to extract and validate microscopic electrostatics of the entire folding energy landscape.     

pKa values and the pH dependent stability of the N-terminal domain of L9 as probes of electrostatic interactions in the denatured state. Differentiation between local and nonlocal interactions

pKa values were measured for the 6 carboxylates in the N-terminal domain of L9 (NTL9) by following NMR chemical shifts as a function of pH. The contribution of each carboxylate to the pH dependent stability of NTL9 was estimated by comparing the pKa values for the native and denatured state of the protein. A set of peptides with sequences derived from NTL9 were used to model the denatured state. In the protein fragments, the pKa values measured for the aspartates varied between 3.8 and 4.1 and the pKa values measured for the glutamates varied between 4.1 and 4.6. These results indicate that the local sequence can significantly influence pKa values in the denatured state and highlight the difficulties in using standard pKa values derived from small compounds. Calculations based on the measured pKa values suggest that the free energy of unfolding of NTL9 should decrease by 4.4 kcal mol-1 when the pH is lowered from 6 to 2. In contrast, urea and thermal denaturation experiments indicate that the stability of the protein decreases by only 2.6 kcal mol-1 when the carboxylates are protonated. This discrepancy indicates that the protein fragments are not a complete representation of the denatured state and that nonlocal sequence effects perturb the pKa's in the denatured state. Increasing the salt concentration from 100 to 750 mM NaCl removes the discrepancy between the stabilities derived from denaturation experiments and the stability changes calculated from the pKa values. At high concentrations of salt there is also less variation of the pKa values measured in the protein fragments. Our results argue that in the denatured state of NTL9 there are electrostatic interactions between groups both local and nonlocal in primary sequence.     

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.     

Electrostatic screening and backbone preferences of amino acid residues in urea-denatured ubiquitin

Local structures in denatured proteins may be important in guiding a polypeptide chain during the folding and misfolding processes. Existence of local structures in chemically denatured proteins is a highly controversial issue. NMR parameters [coupling constants (3) J(H(alpha),H(N)) and chemical shifts] of chemically denatured proteins in general deviate little from their values in small peptides. These peptides were presumed to be completely unstructured; therefore, it was considered that chemically denatured proteins are random coils. But recent experimental studies show that small peptides adopt relatively stable structures in aqueous solutions. Small deviations of the NMR parameters from their values in small peptides may thus actually indicate the existence of local structures in chemically denatured proteins. Using NMR data and theoretical predictions we show here that fluctuating beta-strands exist in urea-denatured ubiquitin (8 M urea at pH 2). Residues in such beta-strands populate more frequently the left side of the broad beta region of -psi space. Urea-denatured ubiquitin contains no detectable beta-sheet secondary structures; nevertheless, the fluctuating beta-strands in urea-denatured ubiquitin coincide to the beta-strands in the native state. Formation of beta-strands is in accord with the electrostatic screening model of unfolded proteins. The free energy of a residue in an unfolded protein is in this model determined by the local backbone electrostatics and its screening by backbone solvation. These energy terms introduce strong electrostatic coupling between neighboring residues, which causes cooperative formation of beta-strands in denatured proteins. We propose that fluctuating beta-strands in denatured proteins may serve as initiation sites to form fibrils.     

Perturbations of the denatured state ensemble: modeling their effects on protein stability and folding kinetics.

By considering the denatured state of a protein as an ensemble of conformations with varying numbers of sequence-specific interactions, the effects on stability, folding kinetics, and aggregation of perturbing these interactions can be predicted from changes in the molecular partition function. From general considerations, the following conclusions are drawn: (1) A perturbation that enhances a native interaction in denatured state conformations always increases the stability of the native state. (2) A perturbation that promotes a non-native interaction in the denatured state always decreases the stability of the native state. (3) A change in the denatured state ensemble can alter the kinetics of aggregation and folding. (4) The loss (or increase) in stability accompanying two mutations, each of which lowers (or raises) the free energy of the denatured state, will be less than the sum of the effects of the single mutations, except in cases where both mutations affect the same set of partially folded conformations. By modeling the denatured state as the ensemble of all non-native conformations of hydrophobic-polar (HP) chains configured on a square lattice, it can be shown that the stabilization obtained from enhancement of native interactions derives in large measure from the avoidance of non-native interactions in the D state. In addition, the kinetic effects of fixing single native contacts in the denatured state or imposing linear gradients in the HH contact probabilities are found, for some sequences, to significantly enhance the efficiency of folding by a simple hydrophobic zippering algorithm. Again, the dominant mechanism appears to be avoidance of non-native interactions. These results suggest stabilization of native interactions and imposition of gradients in the stability of local structure are two plausible mechanisms involving the denatured state that could play a role in the evolution of protein folding and stability.     

Effect of a concentrated "inert" macromolecular cosolute on the stability of a globular protein with respect to denaturation by heat and by chaotropes: a statistical-thermodynamic model

An equilibrium statistical-thermodynamic model for the effect of volume exclusion arising from high concentrations of stable macromolecules upon the stability of a trace globular protein with respect to denaturation by heat and by chaotropes is presented. The stable cosolute and the native form of the trace protein are modeled by effective hard spherical particles. The denatured state of the trace protein is represented as an ensemble of substates modeled by random coils having the same contour length but different rms end-to-end distances (i.e., different degrees of compaction). The excess or nonideal chemical potential of the native state and of each denatured substate is calculated as a function of the concentration of stable cosolute, leading to an estimate of the relative abundance of each state and substate, and the ensemble average free energy of the transition between native and denatured protein. The effect of the addition of stable cosolute upon the temperature of half-denaturation and upon the concentration of chaotrope required to half-denature the tracer at constant temperature is then estimated. At high cosolute concentration (>100 g/l) these effects are predicted to be large and readily measurable experimentally, provided that an experimental system exhibiting a fully reversible unfolding equilibrium at high total macromolecular concentration can be developed.     

Hydrogen exchange in native and denatured states of hen egg-white lysozyme.

The hydrogen exchange kinetics of 68 individual amide protons in the native state of hen lysozyme have been measured at pH 7.5 and 30 degrees C by 2D NMR methods. These constitute the most protected subset of amides, with exchange half lives some 10(5)-10(7) times longer than anticipated from studies of small model peptides. The observed distribution of rates under these conditions can be rationalized to a large extent in terms of the hydrogen bonding of individual amides and their burial from bulk solvent. Exchange rates have also been measured in a reversibly denatured state of lysozyme; this was made possible under very mild conditions, pH 2.0 35 degrees C, by lowering the stability of the native state through selective cleavage of the Cys-6-Cys-127 disulfide cross-link (CM6-127 lysozyme). In this state the exchange rates for the majority of amides approach, within a factor of 5, the values anticipated from small model peptides. For a few amides, however, there is evidence for significant retardation (up to nearly 20-fold) relative to the predicted rates. The pattern of protection observed under these conditions does not reflect the behavior of the protein under strongly native conditions, suggesting that regions of native-like structure do not persist significantly in the denatured state of CM6-127 lysozyme. The pattern of exchange rates from the native protein at high temperature, pH 3.8 69 degrees C, resembles that of the acid-denatured state, suggesting that under these conditions the exchange kinetics are dominated by transient global unfolding. The rates of folding and unfolding under these conditions were determined independently by magnetization transfer NMR methods, enabling the intrinsic exchange rates from the denatured state to be deduced on the basis of this model, under conditions where the predominant equilibrium species is the native state. Again, in the case of most amides these rates showed only limited deviation from those predicted by a simple random coil model. This reinforces the view that these denatured states of lysozyme have little persistent residual order and contrasts with the behavior found for compact partially folded states of proteins, including an intermediate detected transiently during the refolding of hen lysozyme.     

Mutational analysis demonstrates that specific electrostatic interactions can play a key role in the denatured state ensemble of proteins

The nature of the denatured state ensemble has been controversial for decades owing, in large part, to the difficulty in characterizing the structure and energetics of denatured state interactions. There is increasing evidence for relatively non-specific hydrophobic clustering in the denatured states of some proteins but other types of interactions are much less well characterized. Here, we report the characterization of highly specific electrostatic interactions in the denatured state of a small alpha-beta protein, the N-terminal domain of the ribosomal protein L9 (NTL9). Mutation of Lys12 to Met has been shown to increase the stability of NTL9 significantly through the disruption of denatured state interactions. Here, we describe the analysis of the pH-dependent stability of 13 mutants designed to probe the nature of the Lys12 denatured state interaction. Lys12 is located in a lysine-rich region of the protein but analysis of a set of Lys to Met mutants shows that it plays a unique role in the denatured state. Analysis of mutants of all of the acidic residues in NTL9 shows that Lys12 forms a specific non-native electrostatic interaction with Asp8 in the denatured state ensemble. Thus the distribution of charge-charge interactions in the denatured state ensemble of NTL9 appears to be biased by few key interactions and is very different from that expected in a random coil. We propose that these interactions are not encoded by local sequence effects but rather reflect interactions among residues more distant in sequence. These results demonstrate that electrostatic as well as hydrophobic interactions can play an important role in the denatured state ensemble.     

Denatured state aggregation parameters derived from concentration dependence of protein stability

Protein aggregation is a major issue affecting the long-term stability of protein preparations. Proteins exist in equilibrium between the native and denatured or partially denatured conformations. Often denatured or partially denatured conformations are prone to aggregate because they expose to solvent the hydrophobic core of the protein. The aggregation of denatured protein gradually shifts the protein equilibrium toward increasing amounts of denatured and ultimately aggregated protein. Recognizing and quantitating the presence of denatured protein and its aggregation at the earliest possible time will bring enormous benefits to the identification and selection of optimal solvent conditions or the engineering of proteins with the best stability/aggregation profile. In this article, a new approach that allows simultaneous determination of structural stability and the amount of denatured and aggregated protein is presented. This approach is based on the analysis of the concentration dependence of the Gibbs energy (ΔG) of protein stability. It is shown that three important quantities can be evaluated simultaneously: (i) the population of denatured protein, (ii) the population of aggregated protein, and (iii) the fraction of denatured protein that is aggregated.     

Network and graph analyses of folding free energy surfaces

Protein folding is governed by a complex free energy surface whose entropic contributions are relevant because of the large number of degrees of freedom involved. Such complexity, in particular the conformational heterogeneity of the denatured state, is hidden in projections onto one or two order parameters (e.g. fraction of native contacts and/or radius of gyration), which usually results in relatively smooth surfaces. Recent approaches borrowed from network and graph theory have yielded quantitative unprojected representations of the free energy surfaces of a beta-hairpin and a three-stranded beta-sheet peptide using equilibrium folding-unfolding molecular dynamics simulations. Interestingly, the network and graph analyses of these structured peptides have revealed a very heterogeneous denatured state ensemble. It includes high-enthalpy, high-entropy conformations with fluctuating non-native secondary structure, as well as low-enthalpy, low-entropy traps.     

Osmolyte effects on kinetics of FKBP12 C22A folding coupled with prolyl isomerization

Unfolding and refolding kinetics of human FKBP12 C22A were monitored by fluorescence emission over a wide range of urea concentration in the presence and absence of protecting osmolytes glycerol, proline, sarcosine and trimethylamine-N-oxide (TMAO). Unfolding is well described by a mono-exponential process, while refolding required a minimum of two exponentials for an adequate fit throughout the urea concentration range considered. The bi-exponential behavior resulted from complex coupling between protein folding, and prolyl isomerization in the denatured state in which the urea-dependent rate constant for folding was greater than, equal to, and less than the rate constants for prolyl isomerization within the urea concentration range of zero to five molar. Amplitudes and the observed folding and unfolding rate constants were fitted to a reversible three-state model composed of two sequential steps involving the native state and a folding-competent denatured species thermodynamically linked to a folding-incompetent denatured species. Excellent agreement between thermodynamic parameters for FKBP12 C22A folding calculated from the kinetic parameters and those obtained directly from equilibrium denaturation assays provides strong support for the applicability of the mechanism, and provides evidence that FKBP12 C22A folding/unfolding is two-state, with prolyl isomer heterogeneity in the denatured ensemble. Despite the chemical diversity of the protecting osmolytes, they all exhibit the same kinetic behavior of increasing the rate constant of folding and decreasing the rate constant for unfolding. Osmolyte effects on folding/unfolding kinetics are readily explained in terms of principles established in understanding osmolyte effects on protein stability. These principles involve the osmophobic effect, which raises the Gibbs energy of the denatured state due to exposure of peptide backbone, thereby increasing the folding rate. This effect also plays a key role in decreasing the unfolding rate when, as is often the case, the activated complex exposes more backbone than is exposed in the native state.     

The unfolding action of GroEL on a protein substrate

A molecular dynamics simulation of the active unfolding of denatured rhodanese by the chaperone GroEL is presented. The compact denatured protein is bound initially to the cis cavity and forms stable contacts with several of the subunits. As the cis ring apical domains of GroEL undergo the transition from the closed to the more open (ATP-bound) state, they exert a force on rhodanese that leads to the increased unfolding of certain loops. The contacts between GroEL and rhodanese are analyzed and their variation during the GroEL transition is shown. The major contacts, which give rise to the stretching force, are found to be similar to those observed in crystal structures of peptides bound to the apical domains. The results of the simulation show that multidomain interactions play an essential role, in accord with experiments. Implications of the results for mutation experiments and for the action of GroEL are discussed.     

Conformational analysis of peptide fragments derived from the peripheral subunit-binding domain from the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus: evidence for nonrandom structure in the unfolded state

There is currently a great deal of interest in the early events in protein folding. Two issues that have generated particular interest are the nature of the unfolded state under native conditions and the role of local interactions in folding. Here, we report the results of a study of a set of peptides derived from a small two-helix protein, the peripheral subunit-binding domain of the pyruvate dehydrogenase multienzyme complex. Five peptides of overlapping sequence were prepared, including sequences corresponding to each of the helices and to the region connecting them. The peptides were characterized by CD and, where possible, nmr. A peptide corresponding to the second helix is between 12 and 17% helical at neutral pH. CD also indicates a lower percentage of helical structure in the peptide corresponding to the first alpha-helix, although the values of the alpha-proton chemical shifts suggest some preference for nonrandom structure. Peptides corresponding to the interhelical loop, which in the full domain contains two overlapping beta-turns and a 5-residue 3(10)-helix, are less structured. There is no significant change in the helicity of any of these peptides with pH. To test for fragment complementation, CD spectra of the two peptides derived from each helix and the long connecting peptide were compared to the spectra of each possible pair, as well as to a mixture containing all three. No increase in structure was observed. We complement our peptide studies by characterizing a point mutant, D34V, which disrupts a critical hydrogen bonding network. This mutant is unable to fold and provides a useful model of the denatured state. The mutant is between 9 and 16% helical as judged by CD. The modest amount of helical structure formed in some of the peptide fragments and in the point mutant suggests that the denatured state of the peripheral subunit binding domain is not completely unstructured. This may contribute to the very rapid folding observed for the intact protein.     

Disulfide bond effects on protein stability: designed variants of Cucurbita maxima trypsin inhibitor-V

Attempts to increase protein stability by insertion of novel disulfide bonds have not always been successful. According to the two current models, cross-links enhance stability mainly through denatured state effects. We have investigated the effects of removal and addition of disulfide cross-links, protein flexibility in the vicinity of a cross-link, and disulfide loop size on the stability of Cucurbita maxima trypsin inhibitor-V (CMTI-V; 7 kD) by differential scanning calorimetry. CMTI-V offers the advantage of a large, flexible, and solvent-exposed loop not involved in extensive intra-molecular interactions. We have uncovered a negative correlation between retention time in hydrophobic column chromatography, a measure of protein hydrophobicity, and melting temperature (T(m)), an indicator of native state stabilization, for CMTI-V and its variants. In conjunction with the complete set of thermodynamic parameters of denaturation, this has led to the following deductions: (1) In the less stable, disulfide-removed C3S/C48S (Delta Delta G(d)(50 degrees C) = -4 kcal/mole; Delta T(m) = -22 degrees C), the native state is destabilized more than the denatured state; this also applies to the less-stable CMTI-V* (Delta Delta G(d)(50 degrees C) = -3 kcal/mole; Delta T(m) = -11 degrees C), in which the disulfide-containing loop is opened by specific hydrolysis of the Lys(44)-Asp(45) peptide bond; (2) In the less stable, disulfide-inserted E38C/W54C (Delta Delta G(d)(50 degrees C) = -1 kcal/mole; Delta T(m) = +2 degrees C), the denatured state is more stabilized than the native state; and (3) In the more stable, disulfide-engineered V42C/R52C (Delta Delta G(d)(50 degrees C) = +1 kcal/mole; Delta T(m) = +17 degrees C), the native state is more stabilized than the denatured state. These results show that a cross-link stabilizes both native and denatured states, and differential stabilization of the two states causes either loss or gain in protein stability. Removal of hydrogen bonds in the same flexible region of CMTI-V resulted in less destabilization despite larger changes in the enthalpy and entropy of denaturation. The effect of a cross-link on the denatured state of CMTI-V was estimated directly by means of a four-state thermodynamic cycle consisting of native and denatured states of CMTI-V and CMTI-V*. Overall, the results show that an enthalpy-entropy compensation accompanies disulfide bond effects and protein stabilization is profoundly modulated by altered hydrophobicity of both native and denatured states, altered flexibility near the cross-link, and residual structure in the denatured state.     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

Denatured-state energy landscapes of a protein structural database reveal the energetic determinants of a framework model for folding

Position-specific denatured-state thermodynamics were determined for a database of human proteins by use of an ensemble-based model of protein structure. The results of modeling denatured protein in this manner reveal important sequence-dependent thermodynamic properties in the denatured ensembles as well as fundamental differences between the denatured and native ensembles in overall thermodynamic character. The generality and robustness of these results were validated by performing fold-recognition experiments, whereby sequences were matched with their respective folds based on amino acid propensities for the different energetic environments in the protein, as determined through cluster analysis. Correlation analysis between structure and energetic information revealed that sequence segments destined for β-sheet in the final native fold are energetically more predisposed to a broader repertoire of states than are sequence segments destined for α-helix. These results suggest that within the subensemble of mostly unstructured states, the energy landscapes are dominated by states in which parts of helices adopt structure, whereas structure formation for sequences destined for β-strand is far less probable. These results support a framework model of folding, which suggests that, in general, the denatured state has evolutionarily evolved to avoid low-energy conformations in sequences that ultimately adopt β-strand. Instead, the denatured state evolved so that sequence segments that ultimately adopt α-helix and coil will have a high intrinsic structure formation capability, thus serving as potential nucleation sites.