Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Folding of a designed simple ankyrin repeat protein

Ankyrin repeats (AR) are 33-residue motifs containing a beta-turn, followed by two alpha-helices connected by a loop. AR occur in tandem arrangements and stack side-by-side to form elongated domains involved in very different cellular tasks. Recently, consensus libraries of AR repeats were constructed. Protein E1_5 represents a member of the shortest library, and consists of only a single consensus repeat flanked by designed N- and C-terminal capping repeats. Here we present a biophysical characterization of this AR domain. The protein is compactly folded, as judged from the heat capacity of the native state and from the specific unfolding enthalpy and entropy. From spectroscopic data, thermal and urea-induced unfolding can be modeled by a two-state transition. However, scanning calorimetry experiments reveal a deviation from the two-state behavior at elevated temperatures. Folding and unfolding at 5 degrees C both follow monoexponential kinetics with k(folding) = 28 sec(-1) and k(unfolding) = 0.9 sec(-1). Kinetic and equilibrium unfolding parameters at 5 degrees C agree very well. We conclude that E1_5 folds in a simple two-state manner at low temperatures while equilibrium intermediates become populated at higher temperatures. A chevron-plot analysis indicates that the protein traverses a very compact transition state along the folding/unfolding pathway. This work demonstrates that a designed minimal ankyrin repeat protein has the thermodynamic and kinetic properties of a compactly folded protein, and explains the favorable properties of the consensus framework.     

Folding and association of an extremely stable dimeric protein from Sulfolobus islandicus

ORF56 is a plasmid-encoded protein from Sulfolobus islandicus, which probably controls the copy number of the pRN1 plasmid by binding to its own promotor. The protein showed an extremely high stability in denaturant, heat, and pH-induced unfolding transitions, which can be well described by a two-state reaction between native dimers and unfolded monomers. The homodimeric character of native ORF56 was confirmed by analytical ultracentrifugation. Far-UV circular dichroism and fluorescence spectroscopy gave superimposable denaturant-induced unfolding transitions and the midpoints of both heat as well as denaturant-induced unfolding depend on the protein concentration supporting the two-state model. This model was confirmed by GdmSCN-induced unfolding monitored by heteronuclear 2D NMR spectroscopy. Chemical denaturation was accomplished by GdmCl and GdmSCN, revealing a Gibbs free energy of stabilization of -85.1 kJ/mol at 25 degrees C. Thermal unfolding was possible only above 1 M GdmCl, which shifted the melting temperature (t(m)) below the boiling point of water. Linear extrapolation of t(m) to 0 M GdmCl yielded a t(m) of 107.5 degrees C (5 microM monomer concentration). Additionally, ORF56 remains natively structured over a remarkable pH range from pH 2 to pH 12. Folding kinetics were followed by far-UV CD and fluorescence after either stopped-flow or manual mixing. All kinetic traces showed only a single phase and the two probes revealed coincident folding rates (k(f), k(u)), indicating the absence of intermediates. Apparent first-order refolding rates depend linearly on the protein concentration, whereas the unfolding rates do not. Both lnk(f) and lnk(u) depend linearly on the GdmCl concentration. Together, folding and association of homodimeric ORF56 are concurrent events. In the absence of denaturant ORF56 refolds fast (7.0 x 10(7)M(-1)s(-1)) and unfolds extremely slowly (5.7 year(-1)). Therefore, high stability is coupled to a slow unfolding rate, which is often observed for proteins of extremophilic organisms.     

Kinetic analysis of folding and unfolding the 56 amino acid IgG-binding domain of streptococcal protein G.

The 56 amino acid B domain of protein G (GB) is a stable globular folding unit with no disulfide cross-links. The physical properties of GB offer extraordinary flexibility for evaluating the energetics of the folding reaction. The protein is monomeric and very soluble in both folded and unfolded forms. The folding reaction has been previously examined by differential scanning calorimetry (Alexander et al., 1992) and found to exhibit two-state unfolding behavior over a wide pH range with an unfolding transition near 90 degrees C (GB1) at neutral pH. Here, the kinetics of folding and unfolding two naturally occurring versions of GB have been measured using stopped-flow mixing methods and analyzed according to transition-state theory. GB contains no prolines, and the kinetics of folding and unfolding can be fit to a single, first-order rate constant over the temperature range of 5-35 degrees C. The major thermodynamic changes going from the unfolded state to the transition state are (1) a large decrease in heat capacity (delta Cp), indicating that the transition state is compact and solvent inaccessible relative to the unfolded state; (2) a large loss of entropy; and (3) a small increase in enthalpy. The most surprising feature of the folding of GB compared to that of previously studied proteins is that its folding approximates a rapid diffusion controlled process with little increase in enthalpy going from the unfolded to the transition state.     

Very fast folding and association of a trimerization domain from bacteriophage T4 fibritin

The foldon domain constitutes the C-terminal 30 amino acid residues of the trimeric protein fibritin from bacteriophage T4. Its function is to promote folding and trimerization of fibritin. We investigated structure, stability and folding mechanism of the isolated foldon domain. The domain folds into the same trimeric beta-propeller structure as in fibritin and undergoes a two-state unfolding transition from folded trimer to unfolded monomers. The folding kinetics involve several consecutive reactions. Structure formation in the region of the single beta-hairpin of each monomer occurs on the submillisecond timescale. This reaction is followed by two consecutive association steps with rate constants of 1.9(+/-0.5)x10(6)M(-1)s(-1) and 5.4(+/-0.3)x10(6)M(-1)s(-1) at 0.58 M GdmCl, respectively. This is similar to the fastest reported bimolecular association reactions for folding of dimeric proteins. At low concentrations of protein, folding shows apparent third-order kinetics. At high concentrations of protein, the reaction becomes almost independent of protein concentrations with a half-time of about 3 ms, indicating that a first-order folding step from a partially folded trimer to the native protein (k=210 +/- 20 s(-1)) becomes rate-limiting. Our results suggest that all steps on the folding/trimerization pathway of the foldon domain are evolutionarily optimized for rapid and specific initiation of trimer formation during fibritin assembly. The results further show that beta-hairpins allow efficient and rapid protein-protein interactions during folding.     

pH-dependent stability and folding kinetics of a protein with an unusual alpha-beta topology: the C-terminal domain of the ribosomal protein L9

The folding kinetics and thermodynamics of the isolated C-terminal domain of the ribosomal protein L9 (CTL9) have been studied as a function of pH. CTL9 is an alpha-beta protein that contains a single beta-sheet with an unusual mixed parallel, anti-parallel topology. The folding is fully reversible and two-state over the entire pH range. Stopped-flow fluorescence and CD experiments yield the same folding rate, and the chevron plots have the characteristic V-shape expected for two-state folding. The values of DeltaG*(H2O) and the m value calculated from the kinetic experiments are in excellent agreement with the equilibrium measurements. The extrapolated initial amplitudes of both the stopped-flow fluorescence and CD measurements show that there is no detectable burst phase intermediate. The domain contains three histidine residues, two of which are largely buried in the native state. They do not participate in salt-bridges or take part in a hydrogen bonded network. NMR measurements reveal that the buried histidine residues have significantly perturbed pK(a) values in the native state. The equilibrium stability and the folding rate are found to be strongly dependent upon their ionization state. There is a linear relationship between the log of the folding rate and DeltaG* (H2O) . The protein is much more stable and folds noticeably faster at pH values above the native state pK(a) values. DeltaG*(H2O) of unfolding increases from 2.90 kcal mol(-1) at pH 5.0 to 6.40 kcal mol(-1) at pH 8.0 while the folding rate increases from 0.60 to 18.7 s(-1). Tanford linkage analysis revealed that the interactions involving the two histidine residues are largely developed in the transition state. The results are compared to other studies of the pH-dependence of folding.     

Thermodynamics, kinetics, and salt dependence of folding of YopM, a large leucine-rich repeat protein

Small globular proteins have many contacts between residues that are distant in primary sequence. These contacts create a complex network between sequence-distant segments of secondary structure, which may be expected to promote the cooperative folding of globular proteins. Although repeat proteins, which are composed of tandem modular units, lack sequence-distant contacts, several of considerable length have been shown to undergo cooperative two-state folding. To explore the limits of cooperativity in repeat proteins, we have studied the unfolding of YopM, a leucine-rich repeat (LRR) protein of over 400 residues. Despite its large size and modular architecture (15 repeats), YopM equilibrium unfolding is highly cooperative, and shows a very strong dependence on the concentration of urea. In contrast, kinetic studies of YopM folding indicate a mechanism that includes one or more transient intermediates. The urea dependence of the folding and unfolding rates suggests a relatively small transition state ensemble. As with the urea dependence, we have found an extreme dependence of the free energy of unfolding on the concentration of salt. This salt dependence likely results from general screening of a large number of unfavorable columbic interactions in the folded state, rather than from specific cation binding.     

Conformational stability and mechanism of folding of ribonuclease T1

Urea and thermal unfolding curves for ribonuclease T1 (RNase T1) were determined by measuring several different physical properties. In all cases, steep, single-step unfolding curves were observed. When these results were analyzed by assuming a two-state folding mechanism, the plots of fraction unfolded protein versus denaturant were coincident. The dependence of the free energy of unfolding, delta G (in kcal/mol), on urea concentration is given by delta G = 5.6 - 1.21 (urea). The parameters characterizing the thermodynamics of unfolding are: midpoint of the thermal unfolding curve, Tm = 48.1 degrees C, enthalpy change at Tm, delta Hm = 97 kcal/mol, and heat capacity change, delta Cp = 1650 cal/mol deg. A single kinetic phase was observed for both the folding and unfolding of RNase T1 in the transition and post-transition regions. However, two slow kinetic phases were observed during folding in the pre-transition region. These two slow phases account for about 90% of the observed amplitude, indicating that a faster kinetic phase is also present. The slow phases probably result from cis-trans isomerization at the 2 proline residues that have a cis configuration in folded RNase T1. These results suggest that RNase T1 folds by a highly cooperative mechanism with no structural intermediates once the proline residues have assumed their correct isomeric configuration. At 25 degrees C, the folded conformation is more stable than the unfolded conformations by 5.6 kcal/mol at pH 7 and by 8.9 kcal/mol at pH 5, which is the pH of maximum stability. At pH 7, the thermodynamic data indicate that the maximum conformational stability of 8.3 kcal/mol will occur at -6 degrees C.     

Modulation of the multistate folding of designed TPR proteins through intrinsic and extrinsic factors

Tetratricopeptide repeats (TPRs) are a class of all alpha-helical repeat proteins that are comprised of 34-aa helix-turn-helix motifs. These stack together to form nonglobular structures that are stabilized by short-range interactions from residues close in primary sequence. Unlike globular proteins, they have few, if any, long-range nonlocal stabilizing interactions. Several studies on designed TPR proteins have shown that this modular structure is reflected in their folding, that is, modular multistate folding is observed as opposed to two-state folding. Here we show that TPR multistate folding can be suppressed to approximate two-state folding through modulation of intrinsic stability or extrinsic environmental variables. This modulation was investigated by comparing the thermodynamic unfolding under differing buffer regimes of two distinct series of consensus-designed TPR proteins, which possess different intrinsic stabilities. A total of nine proteins of differing sizes and differing consensus TPR motifs were each thermally and chemically denatured and their unfolding monitored using differential scanning calorimetry (DSC) and CD/fluorescence, respectively. Analyses of both the DSC and chemical denaturation data show that reducing the total stability of each protein and repeat units leads to observable two-state unfolding. These data highlight the intimate link between global and intrinsic repeat stability that governs whether folding proceeds by an observably two-state mechanism, or whether partial unfolding yields stable intermediate structures which retain sufficient stability to be populated at equilibrium.     

Measuring the stability of partly folded proteins using TMAO

Standard methods for measuring free energy of protein unfolding by chemical denaturation require complete folding at low concentrations of denaturant so that a native baseline can be observed. Alternatively, proteins that are completely unfolded in the absence of denaturant can be folded by addition of the osmolyte trimethylamine N-oxide (TMAO), and the unfolding free energy can then be calculated through analysis of the refolding transition. However, neither chemical denaturation nor osmolyte-induced refolding alone is sufficient to yield accurate thermodynamic unfolding parameters for partly folded proteins, because neither method produces both native and denatured baselines in a single transition. Here we combine urea denaturation and TMAO stabilization as a means to bring about baseline-resolved structural transitions in partly folded proteins. For Barnase and the Notch ankyrin domain, which both show two-state equilibrium unfolding, we found that DeltaG degrees for unfolding depends linearly on TMAO concentration, and that the sensitivity of DeltaG degrees to urea (the m-value) is TMAO independent. This second observation confirms that urea and TMAO exert independent effects on stability over the range of cosolvent concentrations required to bring about baseline-resolved structural transitions. Thermodynamic parameters calculated using a global fit that assumes additive, linear dependence of DeltaG degrees on each cosolvent are similar to those obtained by standard urea-induced unfolding in the absence of TMAO. Finally, we demonstrate the applicability of this method to measurement of the free energy of unfolding of a partly folded protein, a fragment of the full-length Notch ankyrin domain.     

Folding thermodynamics and kinetics of the leucine-rich repeat domain of the virulence factor Internalin B

Although the folding of alpha-helical repeat proteins has been well characterized, much less is known about the folding of repeat proteins containing beta-sheets. Here we investigate the folding thermodynamics and kinetics of the leucine-rich repeat (LRR) domain of Internalin B (InlB), an extracellular virulence factor from the bacterium Lysteria monocytogenes. This domain contains seven tandem leucine-rich repeats, of which each contribute a single beta-strand that forms a continuous beta-sheet with neighboring repeats, and an N-terminal alpha-helical capping motif. Despite its modular structure, InlB folds in an equilibrium two-state manner, as reflected by the identical thermodynamic parameters obtained by monitoring its sigmoidal urea-induced unfolding transition by different spectroscopic probes. Although equilibrium two-state folding is common in alpha-helical repeat proteins, to date, InlB is the only beta-sheet-containing repeat protein for which this behavior is observed. Surprisingly, unlike other repeat proteins exhibiting equilibrium two-state folding, InlB also folds by a simple two-state kinetic mechanism lacking intermediates, aside from the effects of prolyl isomerization on the denatured state. However, like other repeat proteins, InlB also folds significantly more slowly than expected from contact order. When plotted against urea, the rate constants for the fast refolding and single unfolding phases constitute a linear chevron that, when fitted with a kinetic two-state model, yields thermodynamic parameters matching those observed for equilibrium folding. Based on these kinetic parameters, the transition state is estimated to comprise 40% of the total surface area buried upon folding, indicating that a large fraction of the native contacts are formed in the rate-limiting step to folding.     

Folding and unfolding mechanism of highly stable full-consensus ankyrin repeat proteins

Full-consensus designed ankyrin repeat proteins were designed with one to six identical repeats flanked by capping repeats. These proteins express well in Escherichia coli as soluble monomers. Compared to our previously described designed ankyrin repeat protein library, randomized positions have now been fixed according to sequence statistics and structural considerations. Their stability increases with length and is even higher than that of library members, and those with more than three internal repeats are resistant to denaturation by boiling or guanidine hydrochloride. Full denaturation requires their heating in 5 M guanidine hydrochloride. The folding and unfolding kinetics of the proteins with up to three internal repeats were analyzed, as the other proteins could not be denatured. Folding is monophasic, with a rate that is nearly identical for all proteins (∼ 400-800 s^− 1), indicating that essentially the same transition state must be crossed, possibly the folding of a single repeat. In contrast, the unfolding rate decreases by a factor of about 10⁴ with increasing repeat number, directly reflecting thermodynamic stability in these extraordinarily slow denaturation rates. The number of unfolding phases also increases with repeat number. We analyzed the folding thermodynamics and kinetics both by classical two-state and three-state cooperative models and by an Ising-like model, where repeats are considered as two-state folding units that can be stabilized by interacting with their folded nearest neighbors. This Ising model globally describes both equilibrium and kinetic data very well and allows for a detailed explanation of the ankyrin repeat protein folding mechanism.     

Conformational stability of apoflavodoxin.

Flavodoxins are alpha/beta proteins that mediate electron transfer reactions. The conformational stability of apoflavodoxin from Anaboena PCC 7119 has been studied by calorimetry and urea denaturation as a function of pH and ionic strength. At pH > 12, the protein is unfolded. Between pH 11 and pH 6, the apoprotein is folded properly as judged from near-ultraviolet (UV) circular dichroism (CD) and high-field 1H NMR spectra. In this pH interval, apoflavodoxin is a monomer and its unfolding by urea or temperature follows a simple two-state mechanism. The specific heat capacity of unfolding for this native conformation is unusually low. Near its isoelectric point (3.9), the protein is highly insoluble. At lower pH values (pH 3.5-2.0), apoflavodoxin adopts a conformation with the properties of a molten globule. Although apoflavodoxin at pH 2 unfolds cooperatively with urea in a reversible fashion and the fluorescence and far-UV CD unfolding curves coincide, the transition midpoint depends on the concentration of protein, ruling out a simple two-state process at acidic pH. Apoflavodoxin constitutes a promising system for the analysis of the stability and folding of alpha/beta proteins and for the study of the interaction between apoflavoproteins and their corresponding redox cofactors.     

Rapid folding with and without populated intermediates in the homologous four-helix proteins Im7 and Im9

The kinetics and thermodynamics of the folding of the homologous four-helix proteins Im7 and Im9 have been characterised at pH 7.0 and 10 degrees C. These proteins are 60 % identical in sequence and have the same three-dimensional structure, yet appear to fold by different kinetic mechanisms. The logarithm of the folding and unfolding rates of Im9 change linearly as a function of urea concentration and fit well to an equation describing a two-state mechanism (with a folding rate of 1500 s-1, an unfolding rate of 0. 01 s-1, and a highly compact transition state that has approximately 95 % of the native surface area buried). By contrast, there is clear evidence for the population of an intermediate during the refolding of Im7, as indicated by a change in the urea dependence of the folding rate and the presence of a significant burst phase amplitude in the refolding kinetics. Under stabilising conditions (0.25 M Na2SO4, pH 7.0 and 10 degrees C) the folding of Im9 remains two-state, whilst under similar conditions (0.4 M Na2SO4, pH 7.0 and 10 degrees C) the intermediate populated during Im7 refolding is significantly stabilised (KUI=125). Equilibrium denaturation experiments, under the conditions used in the kinetic measurements, show that Im7 is significantly less stable than Im9 (DeltaDeltaG 9.3 kJ/mol) and the DeltaG and m values determined accord with those obtained from the fit to the kinetic data. The results show, therefore, that the population of an intermediate in the refolding of the immunity protein structure is defined by the precise amino acid sequence rather than the global stability of the protein. We discuss the possibility that the intermediate of Im7 is populated due to differences in helix propensity in Im7 and Im9 and the relevance of these data to the folding of helical proteins in general.     

Multiple interactions of a DNA-binding protein in vivo. II. Effects of host mutations on DNA replication of phage T4 gene 32 mutants.

A quantitative analysis has been made of the kinetics of disulphide bond formation, breakage, and rearrangement which occur during the folding and unfolding of the pancreatic trypsin inhibitor. The results have been used to infer the energetics of the protein conformational transitions which accompany each step.The folding transition is shown to be a co-operative process in which all intermediate states with one or two disulphide bonds are unstable relative to the unfolded, reduced protein and that in the fully folded conformation with three disulphide bonds. The approximate two-state nature of the transition at equilibrium was demonstrated experimentally. The folding transition of the inhibitor which has been determined kinetically is therefore analogous to that observed generally with other globular proteins.     

Kinetically robust monomeric protein from a hyperthermophile

Equilibrium and kinetic studies were carried out under denaturation conditions to clarify the energetic features of the high stability of a monomeric protein, ribonuclease HII, from a hyperthermophile, Thermococcus kodakaraensis (Tk-RNase HII). Guanidine hydrochloride (GdnHCl)-induced unfolding and refolding were measured with circular dichroism at 220 nm, and heat-induced denaturation was studied with differential scanning calorimetry. Both GdnHCl- and heat-induced denaturation are very reversible. It was difficult to obtain the equilibrated unfolding curve of Tk-RNase HII below 40 degrees C, because of the remarkably slow unfolding. The two-state unfolding and refolding reactions attained equilibrium at 50 degrees C after 2 weeks. The Gibbs energy change of GdnHCl-induced unfolding (DeltaG(H(2)O)) at 50 degrees C was 43.6 kJ mol(-1). The denaturation temperature in the DSC measurement shifted as a function of the scan rate; the denaturation temperature at a scan rate of 90 degrees C h(-1) was higher than at a scan rate of 5 degrees C h(-1). The unfolding and refolding kinetics of Tk-RNase HII were approximated as a first-order reaction. The ln k(u) and ln k(r) values depended linearly on the denaturant concentration between 10 and 50 degrees C. The DeltaG(H(2)O) value obtained from the rate constant in water using the two-state model at 50 degrees C, 44.5 kJ mol(-1), was coincident with that from the equilibrium study, 43.6 kJ mol(-1), suggesting the two-state folding of Tk-RNase HII. The values for the rate constant in water of the unfolding for Tk-RNase HII were much smaller than those of E. coli RNase HI and Thermus thermophilus RNase HI, which has a denaturation temperature similar to that of Tk-RNase HII. In contrast, little difference was observed in the refolding rates among these proteins. These results indicate that the stabilization mechanism of monomeric protein from a hyperthermophile, Tk-RNase HII, with reversible two-state folding is characterized by remarkably slow unfolding.     

Tuning lambda6-85 towards downhill folding at its melting temperature

The five-helix bundle lambda6-85* is a fast two-state folder. Several stabilized mutants have been reported to fold kinetically near-downhill or downhill. These mutants undergo a transition to two-state folding kinetics when heated. It has been suggested that this transition is caused by increased hydrophobicity at higher temperature. Here we investigate two histidine-containing mutants of lambda6-85* to see if a weaker hydrophobic core can extend the temperature range of downhill folding. The very stable lambdaHA is the fastest-folding lambda repressor to date (k(f)(-1) approximately k(obs)(-1)=2.3 micros at 44 degrees C). It folds downhill at low temperature, but transits back to two-state folding at its unfolding midpoint. lambdaHG has a weakened hydrophobic core. It is less stable than some slower folding mutants of lambda6-85*, and it has more exposed hydrophobic surface area in the folded state. This mutant nonetheless folds very rapidly, and has the non-exponential folding kinetics of an incipient downhill folder even at the unfolding midpoint (k(m)(-1) approximately 2 micros, k(a)(-1)=15 micros at 56 degrees C). We also compare the thermodynamic melting transition of lambdaHG with the nominal two-state folding mutant lambdaQG, which has a similar melting temperature. Unlike lambdaQG, lambdaHG yields fluorescence wavelength-dependent cooperativities and probe-dependent melting temperatures. This result combined with previous work shows that the energy landscapes of lambda repressor mutants support all standard folding mechanisms.     

Unfolding thermodynamics of Trp-cage, a 20 residue miniprotein, studied by differential scanning calorimetry and circular dichroism spectroscopy

Small proteins provide convenient models for computational studies of protein folding and stability, which are usually compared with experimental data. Until recently, the unfolding of Trp-cage was considered to be a two-state process. However, no direct experimental evidence for this has been presented, and in some cases, the contrary has been suggested. To elucidate a detailed unfolding mechanism, we studied the thermodynamics of unfolding of Trp-cage by differential scanning calorimetry (DSC) and circular dichroism (CD) spectroscopy. The observation that at low temperatures only approximately 90-95% of Trp-cage exists in the native conformation presented an analytical challenge. Nevertheless, it was found that the DSC and CD data can be fitted simultaneously to the same set of thermodynamic parameters. The major uncertainty in such a global fit is the heat capacity change upon unfolding, DeltaCp. This can be circumvented by obtaining DeltaCp directly from the difference between heat capacity functions of the native and unfolded states. Using such an analysis it is shown that Trp-cage unfolding can be represented by a two-state model with the following thermodynamic parameters: Tm = 43.9 +/- 0.8 degrees C, DeltaH(Tm) = 56 +/- 2 kJ/mol, DeltaCp = 0.3 +/- 0.1 kJ/(mol.K). Using these thermodynamic parameters it is estimated that Trp-cage is marginally stable at 25 degrees C, DeltaG(25 degrees C) = 3.2 +/- 0.2 kJ/mol, which is only 30% more than the thermal fluctuation energy at this temperature.     

Folding and stability of a primitive protein

We have previously attempted to simulate domain creation in early protein evolution by recombining polypeptide segments from non-homologous proteins, and we have described the structure of one such de novo protein, 1b11, a segment-swapped tetramer with novel architecture. Here, we have analyzed the thermodynamic stability and folding kinetics of the 1b11 tetramer and its monomeric and dimeric intermediates, and of 1b11 mutants with changes at the domain interface. Denatured 1b11 polypeptides fold into transient, folded monomers with marginal stability (DeltaG<1kcalmol(-1)) which convert rapidly ( approximately 6x10(4)M(-1)s(-1)) into dimers (DeltaG=9.8kcal/mol) and then more slowly ( approximately 3M(-1)s(-1)) into tetramers (DeltaG=28kcalmol(-1)). Segment swapping takes place during dimerization, as suggested by mass spectroscopic analysis of covalently linked peptides derived from proteolysis of a disulfide-linked dimer. Our results confirm that segment swapping and associated oligomerization are both powerful ways of stabilizing proteins, and we suggest that this may have been a feature of early protein evolution.