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.     

Thermodynamics of denaturation of hisactophilin, a beta-trefoil protein

Hisactophilin is a histidine-rich pH-dependent actin-binding protein from Dictyostelium discoideum. The structure of hisactophilin is typical of the beta-trefoil fold, a common structure adopted by diverse proteins with unrelated primary sequences and functions. The thermodynamics of denaturation of hisactophilin have been measured using fluorescence- and CD-monitored equilibrium urea denaturation curves, pH-denaturation, and thermal denaturation curves, as well as differential scanning calorimetry. Urea denaturation is reversible from pH 5.7 to pH 9.7; however, thermal denaturation is highly reversible only below pH approximately 6.2. Reversible denaturation by urea and heat is well fit using a two-state transition between the native and the denatured states. Urea denaturation curves are best fit using a quadratic dependence of the Gibbs free energy of unfolding upon urea concentration. Hisactophilin has moderate, roughly constant stability from pH 7.7 to pH 9.7; however, below pH 7.7, stability decreases markedly, most likely due to protonation of histidine residues. Enthalpic effects of histidine ionization upon unfolding also appear to be involved in the occurrence of cold unfolding of hisactophilin under relatively mild solution conditions. The stability data for hisactophilin are compared with data on hisactophilin function, and with data for two other beta-trefoil proteins, human interleukin-1beta, and basic fibroblast growth factor.     

pH and urea dependence of amide hydrogen-deuterium exchange rates in the beta-trefoil protein hisactophilin.

Amide hydrogen/deuterium exchange rates were measured as a function of pH and urea for 37 slowly exchanging amides in the beta-trefoil protein hisactophilin. The rank order of exchange rates is generally maintained under different solution conditions, and trends in the pH and urea dependence of exchange rates are correlated with the rank order of exchange rates. The observed trends are consistent with the expected behavior for exchange of different amides via global and/or local unfolding. Analysis of the pH dependence of exchange in terms of rate constants for structural opening and closing reveals a wide range of rates in different parts of the hisactophilin structure. The slowest exchanging amides have the slowest opening and closing rates. Many of the slowest exchanging amides are located in trefoil 2, but there are also some slow exchanging amides in trefoils 1 and 3. Slow exchangers tend to be near the interface between the beta-barrel and the beta-hairpin triplet portions of this single-domain structure. The pattern of exchange behaviour in hisactophilin is similar to that observed previously in interleukin-1 beta, indicating that exchange properties may be conserved among beta-trefoil proteins. Comparisons of opening and closing rates in hisactophilin with rates obtained for other proteins reveal clear trends for opening rates; however, trends in closing rates are less apparent, perhaps due to inaccuracies in the values used for intrinsic exchange rates in the data fitting. On the basis of the pH and urea dependence of exchange rates and optical measurements of stability and folding, EX2 is the main exchange mechanism in hisactophilin, but there is also evidence for varying levels of EX1 exchange at low and high pH and high urea concentrations. Equilibrium intermediates in which subglobal portions of structure are cooperatively disrupted are not apparent from analysis of the urea dependence of exchange rates. There is, however, a strong correlation between the Gibbs free energy of opening and the denaturant dependence of opening for all amides, which suggests exchange from a continuum of states with different levels of structure. Intermediates are not very prominent either in equilibrium exchange experiments or in quenched-flow kinetic studies; hence, hisactophilin may not form partially folded states as readily as IL-1 beta and other beta-trefoil proteins.     

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 stability of trp aporepressor from Escherichia coli

Equilibrium and kinetic studies of the urea-induced unfolding of trp aporepressor from Escherichia coli were performed to probe the folding mechanism of this intertwined, dimeric protein. The equilibrium unfolding transitions at pH 7.6 and 25 degrees C monitored by difference absorbance, fluorescence, and circular dichroism spectroscopy are coincident within experimental error. All three transitions are well described by a two-state model involving the native dimer and the unfolded monomer; the free energy of folding in the absence of denaturant and under standard-state conditions is estimated to be 23.3 +/- 0.9 kcal/mol of dimer. The midpoint of the equilibrium unfolding transition increases with increasing protein concentration in the manner expected from the law of mass action for the two-state model. We find no evidence for stable folding intermediates. Kinetic studies reveal that unfolding is governed by a single first-order reaction whose relaxation time decreases exponentially with increasing urea concentration and also decreases with increasing protein concentration in the transition zone. Refolding involves at least three phases that depend on both the protein concentration and the final urea concentration in a complex manner. The relaxation time of the slowest of these refolding phases is identical with that for the single phase in unfolding in the transition zone, consistent with the results expected for a reaction that is kinetically reversible. The two faster refolding phases are presumed to arise from slow isomerization reactions in the unfolded form and reflect parallel folding channels.     

Correspondence between anomalous m- and DeltaCp-values in protein folding

Proteins folding according to a classical two-state system characteristically show V-shaped chevron plots. We have previously interpreted the symmetrically curved chevron plot of the protein U1A as denaturant-dependent movements in the position of the transition state ensemble (TSE). S6, a structural analog of U1A, shows a classical V-shaped chevron plot indicative of straightforward two-state kinetics, but the mutant LA30 has a curved unfolding limb, which is most consistent with TSE mobility. The kinetic m-values (derivatives of the rate constants with respect to denaturant concentration) in themselves depend on denaturant concentration. To obtain complementary information about putative mobile TSEs, we have carried out a thermodynamic analysis of the three proteins, based on data for refolding and unfolding over the range 10 degrees C to 70 degrees C. The data at all temperatures can be fitted to two-state model systems. Importantly, for all three proteins the activation heat capacities are, within error, identical to the heat capacities measured in independent experiments under equilibrium conditions. Although the equilibrium heat capacities are essentially invariant with regard to denaturant concentration, the activation heat capacities, similar to the structurally equivalent kinetic m-values, show marked denaturant dependence. Furthermore, the values of beta++ at different denaturant concentrations measured by m-values and by heat capacity values are very similar. These observations are consistent with significant transition state movements within the framework of two-state folding. The basis for TSE movement appears to be enthalpic rather than entropic, suggesting that the binding energy of denaturant-protein interactions is a major determinant of the response of energy landscape contours to changing environments.     

Perchlorate-induced conformational transition of Staphylococcal nuclease: evidence for an equilibrium unfolding intermediate

The sodium perchlorate-induced conformational transition of Staphylococcal nuclease has been monitored by both circular dichroism (CD) and fluorescence spectroscopy. The perchlorate-induced transition is cooperative as observed by both spectroscopic signals. However, the protein loses only about one-third of its native far-UV CD signal at high perchlorate concentrations, indicating that a significant amount of secondary structure remains in the post-transition state. The remaining CD signal can be further diminished in a cooperative manner by the addition of the strong denaturant, urea. Near-UV CD spectra clearly show that the protein loses its tertiary structure in the perchlorate-induced denatured state. The perchlorate-induced transition curves were fit to the standard two-state model and the standard free energy change and m value of the transition are 2.3kcal/mol and 1.8kcal/(molM), respectively. By comparison, the urea-induced unfolding of Staphylococcal nuclease (in the absence of perchlorate) yields an unfolding free energy change, DeltaG(0,un), of 5.6kcal/mol and an m value of 2.3kcal/(molM). Thus, the thermodynamic state obtained in the post-transition region of perchlorate-induced conformation transition has a significantly lower free energy change, a high content of secondary structure, and diminished tertiary structure. These results suggest that the perchlorate-induced denatured state is a partially folded equilibrium state. Whether this intermediate is relevant to the folding/unfolding path under standard conditions is unknown at this time.     

Apparent pKa shifts of titratable residues at high denaturant concentration and the impact on protein stability

Urea and guanidine-hydrochloride (GdnHCl) are frequently used for protein denaturation in order to determine the Gibbs free energy of folding and kinetic folding/unfolding parameters. Constant pH value is applied in the folding/unfolding experiments at different denaturant concentrations and steady protonation state of titratable groups is assumed in the folded and unfolded protein, respectively. The apparent side-chain pKa values of Asp, Glu, His and Lys in the absence and presence of 6 M urea and GdnHCl, respectively, have been determined by 1H-NMR. pKa values of all four residues are up-shifted by 0.3-0.5 pH units in presence of 6 M urea by comparison with pKa values of the residues dissolved in water. In the presence of 6 M GdnHCl, pKa values are down-shifted by 0.2-0.3 pH units in the case of acidic and up-shifted by 0.3-0.5 pH units in the case of basic residues. Shifted pKa values in the presence of denaturant may have a pronounced effect on the outcome of the protein stability obtained from denaturant unfolding experiments.     

Thermodynamic and kinetic stability of a large multi-domain enzyme from the hyperthermophile Aeropyrum pernix

The multi-domain enzyme isocitrate dehydrogenase from the hyperthermophile Aeropyrum pernix was studied by denaturant-induced unfolding. At pH 7.5, changes in circular dichroism ellipticity and intrinsic fluorescence showed a complex unfolding transition, whereas at pH 3.0, an apparently two-state and highly reversible unfolding occurred. Analytical ultracentrifugation revealed the dissociation from dimer to monomer at pH 3.0. The thermodynamic and kinetic stability were studied at pH 3.0 to explore the role of inter-domain interactions independently of inter-subunit interplay on the wild type and R211M, a mutant where a seven-membered inter-domain ionic network has been disrupted. The unfolding and folding transitions occurred at slightly different denaturant concentrations even after prolonged equilibration time. The difference between the folding and the unfolding profiles was decreased in the mutant R211M. The apparent Gibbs free energy decreased approximately 2 kcal/mol and the unfolding rate increased 4.3-fold in the mutant protein, corresponding to a decrease in activation free energy of unfolding of 0.86 kcal/mol. These results suggest that the inter-domain ionic network might be responsible for additional stabilization through a significant kinetic barrier in the unfolding pathway that could also explain the larger difference observed between the folding and unfolding transitions of the wild type.     

Stability and folding properties of a model beta-sheet protein, Escherichia coli CspA.

Although beta-sheets represent a sizable fraction of the secondary structure found in proteins, the forces guiding the formation of beta-sheets are still not well understood. Here we examine the folding of a small, all beta-sheet protein, the E. coli major cold shock protein CspA, using both equilibrium and kinetic methods. The equilibrium denaturation of CspA is reversible and displays a single transition between folded and unfolded states. The kinetic traces of the unfolding and refolding of CspA studied by stopped-flow fluorescence spectroscopy are monoexponential and thus also consistent with a two-state model. In the absence of denaturant, CspA refolds very fast with a time constant of 5 ms. The unfolding of CspA is also rapid, and at urea concentrations above the denaturation midpoint, the rate of unfolding is largely independent of urea concentration. This suggests that the transition state ensemble more closely resembles the native state in terms of solvent accessibility than the denatured state. Based on the model of a compact transition state and on an unusual structural feature of CspA, a solvent-exposed cluster of aromatic side chains, we propose a novel folding mechanism for CspA. We have also investigated the possible complications that may arise from attaching polyhistidine affinity tags to the carboxy and amino termini of CspA.     

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.     

Unfolding rates of barstar determined in native and low denaturant conditions indicate the presence of intermediates

The folding and unfolding rates of the small protein, barstar, have been monitored using stopped-flow measurements of intrinsic tryptophan fluorescence at 25 degrees C, pH 8.5, and have been compared over a wide range of urea and guanidine hydrochloride (GdnHCl) concentrations. When the logarithms of the rates of folding from urea and from GdnHCl unfolded forms are extrapolated linearly with denaturant concentration, the same rate is obtained for folding in zero denaturant. Similar linear extrapolations of rates of unfolding in urea and GdnHCl yield, however, different unfolding rates in zero denaturant, indicating that such linear extrapolations are not valid. It has been difficult, for any protein, to determine unfolding rates under nativelike conditions in direct kinetic experiments. Using a novel strategy of coupling the reactivity of a buried cysteine residue with 5,5'-dithiobis(2-nitrobenzoic acid) (DTNB) to the unfolding reaction of barstar, the global unfolding and refolding rates have now been determined in low denaturant concentrations. The logarithms of unfolding rates obtained at low urea and GdnHCl concentrations show a markedly nonlinear dependence on denaturant concentration and converge to the same unfolding rate in the absence of denaturant. It is shown that the native protein can sample the fully unfolded conformation even in the absence of denaturant. The observed nonlinear dependences of the logarithms of the refolding and unfolding rates observed for both denaturants are shown to be due to the presence of (un)folding intermediates and not due to movements in the position of the transition state with a change in denaturant concentration.     

Folding of staphylococcal nuclease A studied by equilibrium and kinetic circular dichroism spectra

The urea-induced unfolding of staphylococcal nuclease A has been studied by circular dichroism both at equilibrium and by the kinetics of unfolding and refolding (pH 7.0 and 4.5 degrees C), as a function of Ca2+ and thymidine 3',5'-diphosphate (pdTp) concentration. The results are as follows. (1) The unfolding transition is shifted to higher concentrations of urea by Ca2+ and pdTp, and the presence of both ligands further stabilizes the protein. (2) In the first stage of kinetic refolding, the peptide ellipticity changes rapidly within the dead time of stopped-flow measurement (15 ms), indicating accumulation of a transient intermediate. This intermediate is remarkably less stable than those of other globular proteins previously studied. (3) Dependence of the folding and unfolding rate constants on urea concentration indicates that the critical activated state of folding ("transition state") has considerable structural organization. The transition state does not, however, have the capacity to bind Ca2+ and pdTp, as indicated by the effects of these ligands on the unfolding rate constant. (4) There are at least four different phases in the refolding kinetics in native conditions below 1 M urea. In the absence of pdTp, there are two phases in unfolding, while in the presence of pdTp the unfolding kinetics show a single phase. Some characteristics of the transient intermediate and of the transition state for folding are discussed.     

An NMR view of the folding process of a CheY mutant at the residue level

The folding of CheY mutant F14N/V83T was studied at 75 residues by NMR. Fluorescence, NMR, and sedimentation equilibrium studies at different urea and protein concentrations reveal that the urea-induced unfolding of this CheY mutant includes an on-pathway molten globule-like intermediate that can associate off-pathway. The populations of native and denatured forms have been quantified from a series of 15N-1H HSQC spectra recorded under increasing concentrations of urea. A thermodynamic analysis of these data provides a detailed picture of the mutant's unfolding at the residue level: (1) the transition from the native state to the molten globule-like intermediate is highly cooperative, and (2) the unfolding of this state is sequential and yields another intermediate showing a collapsed N-terminal domain and an unfolded C-terminal tail. This state presents a striking similarity to the kinetic transition state of the CheY folding pathway.     

Kinetic stability plays a dominant role in the denaturant-induced unfolding of Erythrina indica lectin

The urea-induced denaturation of dimeric Erythrina indica lectin (EIL) has been studied at pH 7.2 under equilibrium and kinetic conditions in the temperature range of 40-55 degrees C. The structure of EIL is largely unaffected in this temperature range in absence of denaturant, and also in 8 M urea after incubation for 24 h at ambient temperature. The equilibrium denaturation of EIL exhibits a monophasic unfolding transition from the native dimer to the unfolded monomer as monitored by fluorescence, far-UV CD, and size-exclusion FPLC. The thermodynamic parameters determined for the two-state unfolding equilibrium show that the free energy of unfolding (DeltaGu, aq) remains practically same between 40 and 55 degrees C, with a value of 11.8 +/- 0.6 kcal mol(-1) (monomer units). The unfolding kinetics of EIL describes a single exponential decay pattern, and the apparent rate constants determined at different temperatures indicate that the rate of the unfolding reaction increases several fold with increase in temperature. The presence of probe like external metal ions (Mn2+, Ca2+) does not influence the unfolding reaction thermodynamically or kinetically; however, the presence of EDTA affects only kinetics. The present results suggest that the ability of EIL to preserve the structural integrity against the highly denaturing conditions is linked primarily to its kinetic stability, and the synergic action of heat and denaturant is involved in the unfolding of the protein.     

Unusual cold denaturation of a small protein domain

A thermal unfolding study of the 45-residue α-helical domain UBA(2) using circular dichroism is presented. The protein is highly thermostable and exhibits a clear cold unfolding transition with the onset near 290 K without denaturant. Cold denaturation in proteins is rarely observed in general and is quite unique among small helical protein domains. The cold unfolding was further investigated in urea solutions, and a simple thermodynamic model was used to fit all thermal and urea unfolding data. The resulting thermodynamic parameters are compared to those of other small protein domains. Possible origins of the unusual cold unfolding of UBA(2) are discussed.     

Correlation between mutational destabilization of phage T4 lysozyme and increased unfolding rates

The thermodynamics and kinetics of unfolding of 28 bacteriophage T4 lysozyme variants were compared by using urea gradient gel electrophoresis. The mutations studied cause a variety of sequence changes at different residues throughout the polypeptide chain and result in a wide range of thermodynamic stabilities. A striking relationship was observed between the thermodynamic and kinetic effects of the amino acid replacements: All the substitutions that destabilized the native protein by 2 kcal/mol or more also increased the rate of unfolding. The observed increases in unfolding rate corresponded to a decrease in the activation energy of unfolding (delta Gu) at least 35% as large as the decrease in thermodynamic stability (delta Gu). Thus, the destabilizing lesions bring the free energy of the native state closer to that of both the unfolded state and the transition state for folding and unfolding. Since a large fraction of the mutational destabilization is expressed between the transition state and the native conformation, the changes in folding energetics cannot be accounted for by effects on the unfolded state alone. The results also suggest that interactions throughout much of the folded structure are altered in the formation of the transition state during unfolding.     

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.