Role of the C-terminal helix in the folding and stability of yeast phosphoglycerate kinase

Unraveling the Mechanical Unfolding Pathways of a Multidomain Protein: Phosphoglycerate Kinase

Phosphoglycerate kinase (PGK) is a highly conserved enzyme that is crucial for glycolysis. PGK is a monomeric protein composed of two similar domains and has been the focus of many studies for investigating interdomain interactions within the native state and during folding. Previous studies used traditional biophysical methods (such as circular dichroism, tryptophan fluorescence, and NMR) to measure signals over a large ensemble of molecules, which made it difficult to observe transient changes in stability or structure during unfolding and refolding of single molecules. Here, we unfold single molecules of PGK using atomic force spectroscopy and steered molecular dynamic computer simulations to examine the conformational dynamics of PGK during its unfolding process. Our results show that after the initial forced separation of its domains, yeast PGK (yPGK) does not follow a single mechanical unfolding pathway; instead, it stochastically follows two distinct pathways: unfolding from the N-terminal domain or unfolding from the C-terminal domain. The truncated yPGK N-terminal domain unfolds via a transient intermediate, whereas the structurally similar isolated C-terminal domain has no detectable intermediates throughout its mechanical unfolding process. The N-terminal domain in the full-length yPGK displays a strong unfolding intermediate 13% of the time, whereas the truncated domain (yPGKNT) transitions through the intermediate 81% of the time. This effect indicates that the mechanical properties of yPGK cannot be simply deduced from the mechanical properties of its constituents. We also find that Escherichia coli PGK is significantly less mechanically stable as compared to yPGK, contrary to bulk unfolding measurements. Our results support the growing body of observations that the folding behavior of multidomain proteins is difficult to predict based solely on the studies of isolated domains.

Protein Stability, Folding and Misfolding in Human PGK1 Deficiency

Conformational diseases are often caused by mutations, altering protein folding and stability in vivo. We review here our recent work on the effects of mutations on the human phosphoglycerate kinase 1 (hPGK1), with a particular focus on thermodynamics and kinetics of protein folding and misfolding. Expression analyses and in vitro biophysical studies indicate that disease-causing mutations enhance protein aggregation propensity. We found a strong correlation among protein aggregation propensity, thermodynamic stability, cooperativity and dynamics. Comparison of folding and unfolding properties with previous reports in PGKs from other species suggests that hPGK1 is very sensitive to mutations leading to enhance protein aggregation through changes in protein folding cooperativity and the structure of the relevant denaturation transition state for aggregation. Overall, we provide a mechanistic framework for protein misfolding of hPGK1, which is insightful to develop new therapeutic strategies aimed to target native state stability and foldability in hPGK1 deficient patients.

Role of domain interactions in the collective motion of phosphoglycerate kinase

Protein function is governed by the underlying conformational dynamics of the molecule. The experimental and theoretical work leading to contemporary understanding of enzyme dynamics was mostly restricted to the large-scale movements of single-domain proteins. Collective movements resulting from a regulatory interplay between protein domains is often crucial for enzymatic activity. It is not clear, however, how our knowledge could be extended to describe collective near-equilibrium motions of multidomain enzymes. We examined the effect of domain interactions on the low temperature near equilibrium dynamics of the native state, using phosphoglycerate kinase as model protein. We measured thermal activation of tryptophan phosphorescence quenching to explore millisecond-range protein motions. The two protein domains of phosphoglycerate kinase correspond to two dynamic units, but interdomain interactions link the motion of the two domains. The effect of the interdomain interactions on the activation of motions in the individual domains is asymmetric. As the temperature of the frozen protein is increased from the cryogenic, motions of the N domain are activated first. This is a partial activation, however, and the full dynamics of the domain becomes activated only after the activation of the C domain.

Shotgun redox proteomics identifies specifically modified cysteines in key metabolic enzymes under oxidative stress in Saccharomyces cerevisiae

Post-translational redox modification of thiol groups can form the molecular basis of antioxidative protection and redox control. We have implemented a shotgun redox proteomic technique to identify the precise cysteines reversibly oxidised in key proteins. The method was applied to Saccharomyces cerevisiae subjected to peroxide treatment. Enrichment by covalent redox affinity chromatography allowed the isolation of a "redox subpeptidome" that was analysed by LC-MS/MS. Unique peptides containing specific reversibly oxidised cysteines were used to identify over 70 proteins in control and treated samples of which 27 were consistently present in all replicates. In most cases, the redox modification negatively affects their function and slows down their metabolic pathways. Integration of the data provides a snapshot consistent with a metabolic defensive strategy, regulating key enzymes by redox modification, redirecting energy toward ribulose-5-phosphate recycling for NADPH production and antioxidative defence.This generally applicable method has allowed us to discover new redox regulated proteins (DAHP and carbamoylphosphate synthases, Doa1p) and to precisely identify target cysteines in a number of known ones.

The influence of interdomain interactions on the intradomain motions in yeast phosphoglycerate kinase: a molecular dynamics study

A 3-ns molecular dynamics simulation in explicit solvent was performed to examine the inter- and intradomain motions of the two-domain enzyme yeast phosphoglycerate kinase without the presence of substrates. To elucidate contributions from individual domains, simulations were carried out on the complete enzyme as well as on each isolated domain. The enzyme is known to undergo a hinge-bending type of motion as it cycles from an open to a closed conformation to allow the phosphoryl transfer occur. Analysis of the correlation of atomic movements during the simulations confirms hinge bending in the nanosecond timescale: the two domains of the complete enzyme exhibit rigid body motions anticorrelated with respect to each other. The correlation of the intradomain motions of both domains converges, yielding a distinct correlation map in the enzyme. In the isolated domain simulations-in which interdomain interactions cannot occur-the correlation of domain motions no longer converges and shows a very small correlation during the same simulation time. This result points to the importance of interdomain contacts in the overall dynamics of the protein. The secondary structure elements responsible for interdomain contacts are also discussed.

Domain interactions direct misfolding and amyloid formation of yeast phosphoglycerate kinase

There are proteins that are built of two structural domains and are deposited full-length in amyloid plaques formed in various diseases. In spite of the known differences in the mechanisms of folding of single- and multidomain proteins, no published studies can be found that address the role of the domain-domain interactions during misfolding and amyloid formation. By the discovery of the role of domain-domain interactions, here we provide important insight in the submolecular mechanism of amyloid formation. A model system based on yeast phosphoglycerate kinase was designed. This system includes the wild-type yeast phosphoglycerate kinase and single-tryptophan mutants of the individual N and C terminal domains and the complete protein. Electron microscopic measurements proved that amyloid fibrils grow from all mutants under identical conditions as for the wild-type protein. Misfolding and amyloid formation was followed in stopped-flow and manual mixing experiments on the 1 ms to 4 days timescale. Tryptophan fluorescence was used for selective detection of conformational changes accompanying the formation of the amyloidogenic intermediates and the growth of amyloid fibrils. The interactions between the polypeptide chains of the two domains direct the misfolding process from the early steps to the amyloid formation, and influence the final structure. The kinetics of misfolding is different for the individual domains, pointing to the significance of the amino acid sequence. Misfolding of the domains within the complete protein is synchronized indicating that domain-domain interactions direct the misfolding and amyloid formation mechanism.

Random gene dissection: a tool for the investigation of protein structural organization

To investigate the domain structure of proteins and the function of individual domains, proteins are usually subjected to limited proteolysis, followed by isolation of protein fragments and determination of their functions. We have developed an approach we call random gene dissection (RGD) for the identification of functional protein domains and their interdomain regions as well as their in vivo complementing fragments. The approach was tested on a two-domain protein, the type IIS restriction endonuclease BfiI. The collection of BfiI insertional mutants was screened for those that are endonucleolytically active and thus induce the SOS DNA repair response. Sixteen isolated mutants of the wild-type specificity contained insertions that were dispersed in a relatively large region of the target recognition domain. They split the gene into two complementing parts that separately were unable to induce the SOS DNA repair response. In contrast, all 19 mutants of relaxed specificity contained the cassette inserted into a very narrow interdomain region that connects BfiI domains responsible for DNA recognition and for cleavage. As expected, only the N-terminal fragment of BfiI was required to induce SOS response. Our results demonstrate that RGD can be used as a general method to identify complementing fragments and functional domains in enzymes.

Asymmetric effect of domain interactions on the kinetics of folding in yeast phosphoglycerate kinase

The aim of this work is to shed more light on the effect of domain-domain interactions on the kinetics and the pathway of protein folding. A model protein system consisting of several single-tryptophan variants of the two-domain yeast phosphoglycerate kinase (PGK) and its individual domains was studied. Refolding was initiated from the guanidine-unfolded state by stopped-flow or manual mixing and monitored by tryptophan fluorescence from 1 msec to 1000 sec. Denaturant titrations of both individual domains showed apparent two-state unfolding transitions. Refolding kinetics of the individual domains from different denaturant concentrations, however, revealed the presence of intermediate structures during titration for both domains. Refolding of the same domains within the complete protein showed that domain-domain interactions direct the folding of both domains, but in an asymmetric way. Folding of the N domain was already altered within 1 msec, while detectable changes in the folding of the C domain occurred only 60-100 msec after initiating refolding. All mutants showed a hyperfluorescent kinetic intermediate. Both the disappearance of this intermediate and the completion of the folding were significantly faster in the individual N domain than in the complete protein. On the contrary, folding of the individual C domain was slower than in the complete protein. The presence of the C domain directs the refolding of the N domain along a completely different pathway than that of the individual N domain, while folding of the individual C domain follows the same path as within the complete protein.

Tuning the Heterogeneous Early Folding Dynamics of Phosphoglycerate Kinase

We recently reported stretched kinetics during the formation of a collapsed, long-lived intermediate state of the large two-domain enzyme phosphoglycerate kinase (PGK). It was postulated that intrinsic roughness of the energy landscape on the way downhill to the intermediate causes the lack of a single time-scale. Here, we investigate several alternative explanations for stretched refolding dynamics in more detail: tyrosine fluorescence, multiple tryptophan probes, and rate differences between independently folding domains. To this end, we systematically simplify PGK in several steps from the full protein with two tryptophan residues and all tyrosine residues probed, to a single domain with only one tryptophan residue and no tyrosine residue probed. The kinetics in the 10 micros to 10 ms range are revealed by laser-induced temperature-jump relaxation experiments. The isolated N-terminal domain forms an intermediate by nearly single-exponential kinetics, but the isolated C-terminal domain shows strongly non-exponential kinetics. Thus, domain interaction and a cis-proline residue between the two domains are ruled out as the sole contributors to heterogeneity during the earliest folding dynamics of the C-terminal domain. We apply two limiting models for the roughness of the energy landscape. A sequential three-state model lumps all the roughness into a single trap. The "strange kinetics" model with logarithmic oscillations developed by Klafter and co-workers distributes the roughness over a larger number of states. Both models explain our data about equally well, but the coincidental values of rate constants in all of our double-exponential fits, and the absence of a spectroscopic signature distinct from the endpoints of the folding process favors more roughness than can be explained by just a single trap.

A 1.8 A resolution structure of pig muscle 3-phosphoglycerate kinase with bound MgADP and 3-phosphoglycerate in open conformation: new insight into the role of the nucleotide in domain closure

3-phosphoglycerate kinase (PGK) is a typical kinase with two structural domains. The domains each bind one of the two substrates, 3-phosphoglycerate (3-PG) and MgATP. For the phospho-transfer reaction to take place the substrates must be brought closer by a hinge-bending domain closure. Open and closed structures of the enzyme with different relative domain positions have been determined from different species, but a comprehensive description of this conformational transition is yet to be attained. Crystals of pig muscle PGK in complex with MgADP and 3-phosphoglycerate were grown under the conditions which have previously resulted in crystals of the closed, catalytically competent conformation of Trypanosoma brucei PGK. The X-ray structure of the pig muscle ternary complex was determined at 1.8 A and the model was refined to R=20.8% and Rfree=24.1%. Contrary to expectation, however, it represents an essentially open conformation compared to that of T. brucei PGK. In addition, the beta-phosphate group of ADP is mobile in the new structure, in contrast to its well-defined position in T. brucei PGK. An extensive comparison of the ternary complexes from these remote species has been carried out in order to establish general differences between the two conformations and is reported here. A second pair of the open and closed structures was also compared. These analyses have made it possible to define several characteristic changes which accompany the structural transition, in addition to those identified previously: (1) the operation of a hinge at beta-strand L in the inter-domain region which greatly affects the relative domain positions; (2) the rearrangement and movement of helix 8, regulated through the interactions with the nucleotide phosphate; and (3) the existence of another hinge between helix 14 and the rest of the C-terminal part of the chain, which allows fine adjustment of the N-domain position. The main hinge at beta-strand L acts in concert with the C-terminal hinge at helix 7 described previously. Simultaneous interactions of the nucleotide phosphate groups with the loop that precedes helix 8, beta-strand J and the N terminus of helix 13 are required for propagation of the nucleotide effect towards the beta-strand L molecular hinge. A detailed description of the role of nucleotide binding in the hinge operation is presented.

Role of hydrophobic interactions in yeast phosphoglycerate kinase stability

Cold denaturation of yeast phosphoglycerate kinase (yPGK) was investigated by a combination of far UV circular dichroism (CD), steady-state and time-resolved fluorescence, and small angle X-ray scattering. It was shown that cold denaturation of yPGK cannot be accounted for by a simple two-state process and that an intermediate state can be stabilized under mild denaturing conditions. Comparison between far UV CD and fluorescence shows that in this state the protein displays a fluorescence signal corresponding mainly to exposed tryptophans, whereas its CD signal is only partially modified. Comparison with spectroscopic data obtained from a mutant missing the last 12 amino-acids (yPGK delta404) suggests that lowering the temperature mainly results in a destabilization of hydrophobic interactions between the two domains. Small angle X-ray scattering measurements give further information about this stabilized intermediate. At 4 degrees C and in the presence of 0.45 M Gdn-HCl, the main species corresponds to a protein as compact as native yPGK, whereas a significant proportion of ellipticity has been lost. Although various techniques have shown the existence of residual structures in denatured proteins, this is one example of a compact denatured state devoid of its main content in alpha helices.

Folding funnels and conformational transitions via hinge-bending motions

In this article we focus on presenting a broad range of examples illustrating low-energy transitions via hinge-bending motions. The examples are divided according to the type of hinge-bending involved; namely, motions involving fragments of the protein chains, hinge-bending motions involving protein domains, and hinge-bending motions between the covalently unconnected subunits. We further make a distinction between allosterically and nonallosterically regulated proteins. These transitions are discussed within the general framework of folding and binding funnels. We propose that the conformers manifesting such swiveling motions are not the outcome of "induced fit" binding mechanism; instead, molecules exist in an ensemble of conformations that are in equilibrium in solution. These ensembles, which populate the bottoms of the funnels, a priori contain both the "open" and the "closed" conformational isomers. Furthermore, we argue that there are no fundamental differences among the physical principles behind the folding and binding funnels. Hence, there is no basic difference between funnels depicting ensembles of conformers of single molecules with fragment, or domain motions, as compared to subunits in multimeric quaternary structures, also showing such conformational transitions. The difference relates only to the size and complexity of the system. The larger the system, the more complex its corresponding fused funnel(s). In particular, funnels associated with allosterically regulated proteins are expected to be more complicated, because allostery is frequently involved with movements between subunits, and consequently is often observed in multichain and multimolecular complexes. This review centers on the critical role played by flexibility and conformational fluctuations in enzyme activity. Internal motions that extend over different time scales and with different amplitudes are known to be essential for the catalytic cycle. The conformational change observed in enzyme-substrate complexes as compared to the unbound enzyme state, and in particular the hinge-bending motions observed in enzymes with two domains, have a substantial effect on the enzymatic catalytic activity. The examples we review span the lipolytic enzymes that are particularly interesting, owing to their activation at the water-oil interface; an allosterically controlled dehydrogenase (lactate dehydrogenase); a DNA methyltransferase, with a covalently-bound intermediate; large-scale flexible loop motions in a glycolytic enzyme (TIM); domain motion in PGK, an enzyme which is essential in most cells, both for ATP generation in aerobes and for fermentation in anaerobes; adenylate kinase, showing large conformational changes, owing to their need to shield their catalytic centers from water; a calcium-binding protein (calmodulin), involved in a wide range of cellular calcium-dependent signaling; diphtheria toxin, whose large domain motion has been shown to yield "domain swapping;" the hexameric glutamate dehydrogenase, which has been studied both in a thermophile and in a mesophile; an allosteric enzyme, showing subunit motion between the R and the T states (aspartate transcarbamoylase), and the historically well-studied lac repressor. Nonallosteric subunit transitions are also addressed, with some examples (aspartate receptor and BamHI endonuclease). Hence, using this enzyme-catalysis-centered discussion, we address energy funnel landscapes of large-scale conformational transitions, rather than the faster, quasi-harmonic, thermal fluctuations.

Sequential domain refolding of pig muscle 3-phosphoglycerate kinase: kinetic analysis of reactivation

BACKGROUND

Slow refolding of 3-phosphoglycerate kinase is supposed to be caused mainly by its domain structure: folding of the C-terminal domain and/or domain pairing has been suggested to be the rate-limiting step. A slow isomerization has been observed during refolding of the isolated C-terminal proteolytic fragment (larger than the C-domain of about 22 kDa by 5 kDa) of the pig muscle enzyme. Here, the role of this step in the reformation of the active enzyme species is investigated.

RESULTS

The time course of reactivation during refolding of 3-phosphoglycerate kinase or its complementary proteolytic fragments (residues 1-155 and 156-416) exhibits a pronounced lag-phase indicating the formation of an inactive folding intermediate. The whole process, which leads to a high (60-85%) recovery of the enzyme activity, can be described by two consecutive first-order steps (with rate constants 0.012+/-0.0035 and 0.007+/-0.0020 s(-1)). A prior renaturation of the C-fragment restores MgATP binding by the C-domain and abolishes the faster step, allowing the separate observation of the slower step. In accordance with this, refolding of the C-domain as monitored by a change in Trp fluorescence occurs at a rate similar to that of the faster step.

CONCLUSIONS

In addition to the previously observed slow refolding step (0.012 s(-1)) within the C-domain, the occurrence of another slow step (0.007 s(-1)), probably within the N-domain, is detected. The independence of the folding of the C-domain is demonstrated whereas, from the comparative kinetic analysis, independent folding of the N-domain looks less probable. Our data are more compatible with a sequential, rather than random, mechanism and suggest that folding of the C-domain, leading to an inactive intermediate, occurs first, followed by folding of the N-domain.

Purification and characterisation of the phosphoglycerate kinase isoenzymes of Trypanosoma brucei expressed in Escherichia coli

The Trypanosoma brucei phosphoglycerate kinase (PGK) glycosomal and cytosolic isoenzymes have been overexpressed in Escherichia coli and purified to near-homogeneity. Both enzymes were similar to the corresponding natural proteins with respect to their physicochemical and kinetic properties. In addition, a mutant of the glycosomal PGK lacking the 20 amino acid long C-terminal extension was overexpressed and purified. Various properties of this truncated glycosomal PGK were examined and it was found that in some aspects the protein behaved quite differently when compared with its natural counterpart. This was notably the case for the apparent Km for 3-phosphoglyceric acid, its sensitivity to inhibitors and its response to salts and guanidine HCl. However, its Vmax was found to be similar to that of the natural glycosomal PGK. These results suggest that the changes in the C-terminus caused a conformational change effecting the 3-phosphoglyceric acid binding site located at the N-terminal domain of the protein.

Repulsive interparticle interactions in a denatured protein solution revealed by small angle neutron scattering

In order to investigate the effect of concentration in biological processes such as protein folding, small angle neutron scattering measurements were used to determine the second virial coefficient of solutions of both native and strongly denatured phosphoglycerate kinase and the radius of gyration of the protein at zero concentration. The value of the second virial coefficient is a good probe of the non-ideality of a solution. The present results show that the unfolding of the protein leads to a drastic change in the repulsive intermolecular interactions. We conclude that these interactions are due mainly to the behaviour of the denatured polypeptide chain as an excluded volume polymer.

An engineered amino-terminal domain of yeast phosphoglycerate kinase with native-like structure

Previous studies have suggested that the carboxy-terminal peptide (residues 401-415) and interdomain helix (residues 185-199) of yeast phosphoglycerate kinase, a two-domain enzyme, play a role in the folding and stability of the amino-terminal domain (residues 1-184). A deletion mutant has been created in which the carboxy-terminal peptide is attached to the amino-terminal domain (residues 1-184) plus interdomain helix (residues 185-199) through a flexible peptide linker, thus eliminating the carboxy-terminal domain entirely. CD, fluorescence, gel filtration, and NMR experiments indicated that, unlike versions described previously, this isolated N-domain is soluble, monomeric, compactly folded, native-like in structure, and capable of binding the substrate 3-phosphoglycerate with high affinity in a saturable manner. The midpoint of the guanidine-induced unfolding transition was the same as that of the native two-domain protein (Cm approximately 0.8 M). The free energy change associated with guanidine-induced unfolding was one-third that of the native enzyme, in agreement with previous studies that evaluated the intrinsic stability of the N-domain and the contribution of domain-domain interactions to the stability of PGK. These observations suggest that the C-terminal peptide and interdomain helix are sufficient for maintaining a native-like fold of the N-domain in the absence of the C-domain.

A new metal-binding site for yeast phosphoglycerate kinase as determined by the use of a metal-ATP analog

Suicide substrate beta, gamma-bidentate Rh(III)ATP (RhATP) was used to map the metal ion-binding site in yeast phosphoglycerate kinase (PGK). Cleavage of the RhATP-inactivated enzyme with pepsin and subsequent separation of peptides by reverse-phase high-performance liquid chromatography gave two Rh-nucleotide bound peptides. One of the peptides corresponded to the C-terminal residues of PGK, and the other to a part of helix V. Of the four glutamates present in the C-terminal peptide, Glu 398 may be a likely metal coordination site. Therefore, importance of the C-terminal residues in PGK catalysis may be attributed, in part to the coordination of metal ion of the metal-ATP substrate. Metal coordination may then align the C-terminal peptide to extend toward the N-terminal domain and form the "closed" active site. Results presented in this paper suggest that one or more side chains of the enzyme may be coordinated to the metal ion in the PGK.3-phospho-D-glycerate-RhATP complex, and that exchange-inert metal-ATP analogs could be used to determine metal coordination sites on kinases and other metal-ATP-utilizing enzymes.

The carboxy-terminal region of human interleukin-5 is essential for maintenance of tertiary structure but not for dimerization

The C-terminal region of interleukin-5 has previously been suggested to be important for biological activity [Mackenzie et al., (1991), Mol. Immunol. 28, 155-158; Kodama et al. (1991), Biochem. Biophys. Res. Commun. 178, 514-519]. We have investigated this region by making a series of truncation mutants. The proteins were expressed in Escherichia coli, purified from inclusion bodies, and were able to refold with the disulfide homodimeric topology typical of interleukin-5. Analysis of the truncated carboxy-terminal proteins in an interleukin-5-dependent proliferation assay on TF-1 cells showed a rapid loss of activity as the C-terminal was shortened by more than two amino acids. This loss of biological activity correlated with a drop in binding affinity to both the alpha chain of the receptor and the high-affinity complex consisting of the alpha and beta subunits. Analysis of the proteins by 1H-NMR showed that the truncated mutants have higher exchange rates with solvent, indicating a less rigid structure. The carboxy-terminal region is therefore necessary to maintain the stability of the four-helix bundle and to orient correctly the important residues of the fourth helix. Inspection of the structure determined by X-ray crystallography shows that Trp-110 acts as the major residue in anchoring the fourth helix.

The role of the C-terminal lysine in the hinge bending mechanism of yeast phosphoglycerate kinase

Treatment of yeast phosphoglycerate kinase (PGK) with trypsin results in a fourfold increase in the Vmax of this enzyme, without affecting the Km. This activation is shown to be due to the removal of the C-terminal lysine residue. The C-terminal sequence folds back over the N-terminal domain and contacts the extreme N-terminal sequence which folds onto the C-terminal domain, thus making many of the inter-domain contacts in this two domain protein. Previous studies have shown that this C-terminal region is important in mediating the conformational changes required during catalysis by yeast PGK. Observation of the three-dimensional structure of this enzyme suggests that removal of the C-terminal lysine residue will strengthen the interaction between K5 and E413. This indicates that this salt bridge stabilises the enzyme in the higher activity form, while the presence of K415 reduces the strength of that interaction.

Denaturant m values and heat capacity changes: relation to changes in accessible surface areas of protein unfolding

Denaturant m values, the dependence of the free energy of unfolding on denaturant concentration, have been collected for a large set of proteins. The m value correlates very strongly with the amount of protein surface exposed to solvent upon unfolding, with linear correlation coefficients of R = 0.84 for urea and R = 0.87 for guanidine hydrochloride. These correlations improve to R = 0.90 when the effect of disulfide bonds on the accessible area of the unfolded protein is included. A similar dependence on accessible surface area has been found previously for the heat capacity change (delta Cp), which is confirmed here for our set of proteins. Denaturant m values and heat capacity changes also correlate well with each other. For proteins that undergo a simple two-state unfolding mechanism, the amount of surface exposed to solvent upon unfolding is a main structural determinant for both m values and delta Cp.