Robust and convenient analysis of protein thermal and chemical stability

We present the software CDpal that is used to analyze thermal and chemical denaturation data to obtain information on protein stability. The software uses standard assumptions and equations applied to two‐state and various types of three‐state denaturation models in order to determine thermodynamic parameters. It can analyze denaturation monitored by both circular dichroism and fluorescence spectroscopy and is extremely flexible in terms of input format. Furthermore, it is intuitive and easy to use because of the graphical user interface and extensive documentation. As illustrated by the examples herein, CDpal should be a valuable tool for analysis of protein stability.     

Fourier-transform infrared spectroscopic investigation of the thermal denaturation of hen egg-white lysozyme dissolved in aqueous buffer and glycerol

The thermal denaturation of lysozyme dissolved in aqueous phosphate buffer (pH 5.1) and glycerol was studied by Fourier-transform infrared (FTIR) spectroscopy. In both solvents, a single temperature-induced conformational transition was observed but at the distinctly different temperatures of 73 °C in aqueous buffer and 94 ± 2 °C in glycerol. No changes in the secondary structure were observed in glycerol up to 90 °C. Thus, FTIR data were consistent with the formation of a highly ordered molten globule state at temperatures below 90 °C followed by lysozyme unfolding at higher temperatures in glycerol.     

Quantitative analysis of the kinetics of denaturation and renaturation of barstar in the folding transition zone.

The fluorescence-monitored kinetics of folding and unfolding of barstar by guanidine hydrochloride (GdnHCl) in the folding transition zone, at pH 7, 25 degrees C, have been quantitatively analyzed using a 3-state mechanism: U(S)<-->UF<-->N. U(S) and UF are slow-refolding and fast-refolding unfolded forms of barstar, and N is the native protein. U(S) and UF probably differ in possessing trans and cis conformations, respectively, of the Tyr 47-Pro 48 bond. The 3-state model could be used because the kinetics of folding and unfolding of barstar show 2 phases, a fast phase and a slow phase, and because the relative amplitudes of the 2 phases depend only on the final refolding conditions and not on the initial conditions. Analysis of the observed kinetics according to the 3-state model yields the values of the 4 microscopic rate constants that describe the transitions between the 3 states at different concentrations of GdnHCl. The value of the equilibrium unfolded ratio U(S):UF (K21) and the values of the rate constants of the U(S)-->UF and UF-->U(S) reactions, k12 and k21, respectively, are shown to be independent of the concentration of GdnHCl. K21 has a value of 2.1 +/- 0.1, and k12 and k21 have values of 5.3 x 10(-3) s-1 and 11.2 x 10(-3) s-1, respectively. Double-jump experiments that monitor reactions that are silent to fluorescence monitoring were used to confirm the values of K21, k12, and k21 obtained from the 3-state analysis and thereby the validity of the 3-state model.(ABSTRACT TRUNCATED AT 250 WORDS)     

Alpha-crystallin binds to the aggregation-prone molten-globule state of alkaline protease: implications for preventing irreversible thermal denaturation

Alpha-crystallin, the major eye-lens protein with sequence homology with heat-shock proteins (HSPs), acts like a molecular chaperone by suppressing the aggregation of damaged crystallins and proteins. To gain more insight into its chaperoning ability, we used a protease as the model system that is known to require a propeptide (intramolecular chaperone) for its proper folding. The protease ("N" state) from Conidiobolus macrosporus (NCIM 1298) unfolds at pH 2.0 ("U" state) through a partially unfolded "I" state at pH 3.5 that undergoes transition to a molten globule-(MG) like "I(A)" state in the presence of 0.5 M sodium sulfate. The thermally-stressed I(A) state showed complete loss of structure and was prone to aggregation. Alpha-crystallin was able to bind to this state and suppress its aggregation, thereby preventing irreversible denaturation of the enzyme. The alpha-crystallin-bound I(A) state exhibited native-like secondary and tertiary structure showing the interaction of alpha-crystallin with the MG state of the protease. 8-Anilinonaphthalene sulphonate (ANS) binding studies revealed the involvement of hydrophobic interactions in the formation of the complex of alpha-crystallin and protease. Refolding of acid-denatured protease by dilution to pH 7.5 resulted in aggregation of the protein. Unfolding of the protease in the presence of alpha-crystallin and its subsequent refolding resulted in the generation of a near-native intermediate with partial secondary and tertiary structure. Our studies represent the first report of involvement of a molecular chaperone-like alpha-crystallin in the unfolding and refolding of a protease. Alpha-crystallin blocks the unfavorable pathways that lead to irreversible denaturation of the alkaline protease and keeps it in a near-native, folding-competent intermediate state.     

Chemically crosslinked protein dimers: stability and denaturation effects.

Nine single substitution cysteine mutants of staphylococcal nuclease (nuclease) were preferentially crosslinked at the introduced cysteine residues using three different bifunctional crosslinking reagents; 1,6-bismaleimidohexane (BMH), 1,3-dibromo-2-propanol (DBP), and the chemical warfare agent, mustard gas (bis(2-chloroethyl)sulfide; mustard). BMH and mustard gas are highly specific reagents for cysteine residues, whereas DBP is not as specific. Guanidine hydrochloride (GuHCl) denaturations of the resulting dimeric proteins exhibited biphasic unfolding behavior that did not fit the two-state model of unfolding. The monofunctional reagent, epsilon-maleimidocaproic acid (MCA), was used as a control for the effects of alkylation. Proteins modified with MCA unfolded normally, showing that this unusual unfolding behavior is due to crosslinking. The data obtained from these crosslinked dimers was fitted to a three-state thermodynamic model of two successive transitions in which the individual subunits cooperatively unfold. These two unfolding transitions were very different from the unfolding of the monomeric protein. These differences in unfolding behavior can be attributed in large part to changes in the denatured state. In addition to GuHCl titrations, the crosslinked dimers were also thermally unfolded. In contrast to the GuHCl denaturations, analysis of this data fit a two-state model well, but with greatly elevated van't Hoff enthalpies in many cases. However, clear correlations between the thermal and GuHCl denaturations exist, and the differences in thermal unfolding can be rationalized by postulating interactions of the denatured crosslinked proteins.     

Molecular mechanism underlying the thermal stability and pH-induced unfolding of CHABII

The 37-residue alpha/beta protein CHABII was previously demonstrated to undergo a gradual pH-induced unfolding. It has been shown that even at pH 4.0 CHABII still retained a highly native-like secondary structure and tertiary topology although its tight side-chain packing was severely disrupted, typical of the molten globule state. Here, we have expressed and refolded the recombinant proteins of CHABII and its mutant [Phe21]-CHABII, and subsequently conducted extensive CD and NMR characterizations. The results indicated: (1) replacement of His21 by Phe in [Phe21]-CHABII eliminated the pH-induced unfolding from pH 6.5 to 4.0, indicating that His21 was responsible for the observed pH-induced unfolding of CHABII. Further examinations revealed that although the pH-induced unfolding of CHABII was also triggered by the protonation of the His residue as previously uncovered for apomyoglobin, their molecular mechanisms are different. (2) Monitoring the pH-induced unfolding by 1H-15N HSQC spectroscopy allowed us to visualize the gradual development of the CHABII molten globule. At pH 4.0, the HSQC spectrum of CHABII was poorly dispersed with dispersions of approximately 1 ppm over proton dimension and 10 ppm over 15N dimension, characteristic of severely or even "completely unfolded" proteins. One the other hand, unambiguous assignments of the NOESY spectra of CHABII led to the identification of the persistent medium and long-range NOEs at pH 4.0, which define a highly native-like secondary structure and tertiary packing. This implies that the degree of the native-like topology might be underestimated in the previous characterization of partially folded and even completely unfolded proteins. (3) Replacement of His21 by Phe with higher side-chain hydrophobicity only caused a minor structural rearrangement but considerably enhanced the packing interaction of the hydrophobic core, as evident from a dramatic increase in NOE contacts in [Phe21]-CHABII. The enhancement led to an increase of the thermal stability of [Phe21]-CHABII by approximately 17 deg. C.     

Conformational stability of HPr: the histidine-containing phosphocarrier protein from Bacillus subtilis.

The conformational stability of the histidine-containing phosphocarrier protein (HPr) from Bacillus subtilis has been determined using a combination of thermal unfolding and solvent denaturation experiments. The urea-induced denaturation of HPr was monitored spectroscopically at fixed temperatures and thermal unfolding was performed in the presence of fixed concentrations of urea. These data were analyzed in several different ways to afford a measure of the cardinal parameters (delta Hg, Tg, delta Sg, and delta Cp) that describe the thermodynamics of folding for HPr. The method of Pace and Laurents (Pace CN, Laurents DV, 1989, Biochemistry 28:2520-2525) was used to estimate delta Cp as was a global analysis of the thermal- and urea-induced unfolding data. Each method used to analyze the data gives a similar value for delta Cp (1,170 +/- 50 cal mol-1K-1). Despite the high melting temperature for HPr (Tg = 73.5 degrees C), the maximum stability of the protein, which occurs at 26 degrees C, is quite modest (delta Gs = 4.2 kcal mol-1). In the presence of moderate concentrations of urea, HPr exhibits cold denaturation, and thus a complete stability curve for HPr, including a measure of delta Cp, can be achieved using the method of Chen and Schellman (Chen B, Schellman JA, 1989, Biochemistry 28:685-691). A comparison of the different methods for the analysis of solvent denaturation curves is provided and the effects of urea on the thermal stability of this small globular protein are discussed. The methods presented will be of general utility in the characterization of the stability curve for many small proteins.     

Increasing and decreasing protein stability: effects of revertant substitutions on the thermal denaturation of phage lambda repressor

The thermal denaturations of five revertant lambda repressors containing single amino acid substitutions in their N-terminal domains have been studied by differential scanning calorimetry. Two substitutions slightly decrease stability, and the remaining three render the protein more stable than wild type. The Gly48----Asn and Gly48----Ser proteins are 4 degrees C more stable than wild type. These two substitutions replace an alpha helical residue, and in each case a poor helix forming residue, glycine, is replaced by a residue with a higher helical propensity. We also present data showing that one revertant, Tyr22----Phe, has reduced operator DNA binding affinity despite its enhanced stability.     

pH Dependence of structural stability of interleukin-2 and granulocyte colony-stimulating factor

After a cytokine binds to its receptor on the cell surface (pH approximately 7), the complex is internalized into acidic endosomal compartments (pH approximately 5-6), where partially unfolded intermediates can form. The nature of these structural transitions was studied for wild-type interleukin-2 (IL-2) and wild-type granulocyte colony-stimulating factor (G-CSF). A noncoincidence of denaturation transitions in the secondary and tertiary structure of IL-2 and tertiary structural perturbations in G-CSF suggest the presence of an intermediate state for each, a common feature of this structural family of four-helical bundle proteins. Unexpectedly, both IL-2 and G-CSF display monotonic increases in stability as the pH is decreased from 7 to 4. We hypothesize that such cytokines with cell-based clearance mechanisms in vivo may have evolved to help stabilize endosomal complexes for sorting to lysosomal degradation. We show that mutants of both IL-2 and G-CSF have differential stabilities to their wild-type counterparts as a function of pH, and that these differences may explain the differences in ligand trafficking and depletion. Further understanding of the structural changes accompanying unfolding may help guide cytokine design with respect to ligand binding, endocytic trafficking, and, consequently, therapeutic efficacy.     

Molecular basis of the structural stability of hemochromatosis factor E: A combined molecular dynamic simulation and GdmCl‐induced denaturation study

Hemochromatosis factor E (HFE) is a member of class I MHC family and plays a significant role in the iron homeostasis. Denaturation of HFE induced by guanidinium chloride (GdmCl) was measured by monitoring changes in [θ]₂₂₂ (mean residue ellipticity at 222 nm), intrinsic fluorescence emission intensity at 346 nm (F₃₄₆) and the difference absorption coefficient at 287 nm (Δε₂₈₇) at pH 8.0 and 25°C. Coincidence of denaturation curves of these optical properties suggests that GdmCl‐induced denaturation (native (N) state ? denatured (D) state) is a two‐state process. The GdmCl‐induced denaturation was found reversible in the entire concentration range of the denaturant. All denaturation curves were analyzed for , Gibbs free energy change associated with the denaturation equilibrium (N state ? D state) in the absence of GdmCl, which is a measure of HFE stability. We further performed molecular dynamics simulation for 40 ns to see the effect of GdmCl on the structural stability of HFE. A well defined correlation was established between in vitro and in silico studies. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 133–142, 2016.     

Kinetics and thermodynamics of thermal denaturation in acyl carrier protein.

The denaturation of Escherichia coli acyl carrier protein (ACP) in buffers containing both monovalent and divalent cations was followed by variable-temperature NMR and differential scanning calorimetry. Both high concentrations of monovalent salts (Na+) and moderate concentrations of divalent salts (Ca2+) raise the denaturation temperature, but calorimetry indicates that a significant increase in the enthalpy of denaturation is obtained only with the addition of a divalent salt. NMR experiments in both low ionic strength monovalent buffers and low ionic strength monovalent buffers containing calcium ions show exchange between native and denatured forms to be slow on the NMR time scale. However, in high ionic strength monovalent buffers, where the temperature of denaturation is elevated as it is in the presence of Ca2+, the transition is fast on the NMR time scale. These results suggest that monovalent and divalent cations may act to stabilize ACP in different ways. Monovalent ions may nonspecifically balance the intrinsic negative charge of this protein in a way that is similar for native, denatured, and intermediate forms. Divalent cations provide stability by binding to specific sites present only in the native state.     

Cation binding effects on the pH, thermal and urea denaturation transitions in alpha-lactalbumin

The binding of monovalent (Na+, K+) and divalent (Ca2+, Mg2+) cations to bovine alpha-lactalbumin at 20 and 37 degrees C has been studied by means of intrinsic protein fluorescence. The values of apparent binding constants for these ions obtained at 37 degrees C are about one order of magnitude lower than those measured at 20 degrees C. Urea and alkali (pH greater than 10) induce unfolding transitions which involve stable partially unfolded intermediates for all metal ion-bound forms of alpha-lactalbumin. Heating induces similar partially unfolded states. Nevertheless, the partially unfolded states induced by heating, urea, alkaline or acidic treatments are somewhat different in their tryptophan residue environment properties. The results have been interpreted in terms of a simple scheme of equilibria between metal-free and metal-bound forms in their native, partially unfolded and unfolded states. The scheme provides an approach to the quantitative interpretation of any transition equilibrium shift induced by a low molecular mass species able to be bound by a protein.     

Irreversible thermal denaturation of Torpedo californica acetylcholinesterase.

Thermal denaturation of Torpedo californica acetylcholinesterase, a disulfide-linked homodimer with 537 amino acids in each subunit, was studied by differential scanning calorimetry. It displays a single calorimetric peak that is completely irreversible, the shape and temperature maximum depending on the scan rate. Thus, thermal denaturation of acetylcholinesterase is an irreversible process, under kinetic control, which is described well by the two-state kinetic scheme N-->D, with activation energy 131 +/- 8 kcal/mol. Analysis of the kinetics of denaturation in the thermal transition temperature range, by monitoring loss of enzymic activity, yields activation energy of 121 +/- 20 kcal/mol, similar to the value obtained by differential scanning calorimetry. Thermally denatured acetylcholinesterase displays spectroscopic characteristics typical of a molten globule state, similar to those of partially unfolded enzyme obtained by modification with thiol-specific reagents. Evidence is presented that the partially unfolded states produced by the two different treatments are thermodynamically favored relative to the native state.     

Thermodynamics of staphylococcal nuclease denaturation. II. The A-state.

Staphylococcal nuclease, at low pH and in the presence of high salt concentrations, has previously been proposed to exist in a partially folded or molten globule form called the "A-state" (Fink et al., 1993, Protein Sci 2:1155-1160). We have found that the A-state of nuclease at pH 2.1 in the presence of moderate to high salt concentrations and at low temperature exists in a substantially folded form structurally more similar to a native state. The A-state has the far-UV circular dichroism spectra characteristic of the native protein, which indicates that it has a large degree of secondary structure. Upon heating, the A-state denatures with a sigmoidal change in far-UV ellipticity and an observable peak in a differential scanning calorimeter trace, indicating that it is thermodynamically distinct from the denatured state. Three different mutations in a residue normally buried in the protein's core stabilize or destabilize the A-state in the same way as they affect the denaturation of the native state. The A-state must, therefore, contain at least some tertiary packing of side chains. Unlike the native state, which shows cold denaturation at low temperatures, the A-state is most stable at temperatures below 0 degrees C.     

Thermodynamics of staphylococcal nuclease denaturation. I. The acid-denatured state.

Using high-sensitivity differential scanning calorimetry, we reexamined the thermodynamics of denaturation of staphylococcal nuclease. The denaturational changes in enthalpy and heat capacity were found to be functions of both temperature and pH. The denatured state of staphylococcal nuclease at pH 8.0 and high temperature has a heat capacity consistent with a fully unfolded protein completely exposed to solvent. At lower pH values, however, the heat capacity of the denatured state is lower, resulting in a lower delta Cp and delta H for the denaturation reaction. The acid-denatured protein can thus be distinguished from a completely unfolded protein by a defined difference in enthalpy and heat capacity. Comparison of circular dichroism spectra suggests that the low heat capacity of the acid-denatured protein does not result from residual helical secondary structure. The enthalpy and heat capacity changes of denaturation of a less stable mutant nuclease support the observed dependence of delta H on pH.     

pH-induced structural change of a multi-tryptophan protein MPT63 with immunoglobulin-like fold: identification of perturbed tryptophan residue/residues

The structural change of M. tuberculosis MPT63, which is predominantly a β-sheet protein having an immunoglobulin like fold, has been investigated in the pH range 7.5–1.5 using various biophysical techniques along with low-temperature phosphorescence (LTP) spectroscopy. MPT63 contains four Tryptophan (Trp) residues at 26, 48, 82, and 129. Although circular dichroism, steady-state and time-resolved fluorescence, time-resolved anisotropy, 1-aniline-8-naphthalene sulfonic (ANS) acid binding, and analytical ultracentrifuge depict more open largely unfolded structure of MPT63 at pH 1.5 and also more accessible nature of Trp residues to neutral quencher at pH 1.5, it is, however, not possible to assign the specific Trp residue/residues being perturbed. This problem has been resolved using LTP of MPT63, which shows optically resolved four distinct (0, 0) bands corresponding to four Trp residues in the pH range 7.5–3.0. LTP at pH 1.5 clearly reveals that the solvent-exposed Trp 82 and the almost buried Trp 129 are specifically affected compared with Trp 48 and Trp 26. Lys 8 and Lys 27 are predicted to affect Trp 129. Tyrosine residues are found to be silent even at pH 1.5. This type of specific perturbation in a multi-Trp protein has not been addressed before. LTP further indicates that partially exposed Trp 48 is preferentially quenched by acrylamide compared with other Trp residues at both pH 7.5 and 1.5. The solvent-exposed Trp 82 is surprisingly found to be not quenched by acrylamide at pH 7.5. All the results are obtained using micromolar concentration of protein and without using any Trp-substituted mutant.     

Thermodynamics of unfolding for turkey ovomucoid third domain: thermal and chemical denaturation.

We have used thermal and chemical denaturation to characterize the thermodynamics of unfolding for turkey ovomucoid third domain (OMTKY3). Thermal denaturation was monitored spectroscopically at a number of wave-lengths and data were subjected to van't Hoff analysis; at pH 2.0, the midpoint of denaturation (Tm) occurs at 58.6 +/- 0.4 degrees C and the enthalpy of unfolding at this temperature (delta Hm) is 40.8 +/- 0.3 kcal/mol. When Tm was perturbed by varying pH and denaturant concentration, the resulting plots of delta Hm versus Tm yield a mean value of 590 +/- 120 cal/(mol.K) for the change in heat capacity upon unfolding (delta Cp). A global fit of the same data to an equation that includes the temperature dependence for the enthalpy of unfolding yielded a value of 640 +/- 110 cal/(mol.K). We also performed a variation of the linear extrapolation method described by Pace and Laurents, which is an independent method for determining delta Cp (Pace, C.N. & Laurents, D., 1989, Biochemistry 28, 2520-2525). First, OMTKY3 was thermally denatured in the presence of a variety of denaturant concentrations. Linear extrapolations were then made from isothermal slices through the transition region of the denaturation curves. When extrapolated free energies of unfolding (delta Gu) were plotted versus temperature, the resulting curve appeared linear; therefore, delta Cp could not be determined. However, the data for delta Gu versus denaturant concentration are linear over an extraordinarily wide range of concentrations. Moreover, extrapolated values of delta Gu in urea are identical to values measured directly.     

Conformational stability of dimeric proteins: quantitative studies by equilibrium denaturation.

The conformational stability of dimeric globular proteins can be measured by equilibrium denaturation studies in solvents such as guanidine hydrochloride or urea. Many dimeric proteins denature with a 2-state equilibrium transition, whereas others have stable intermediates in the process. For those proteins showing a single transition of native dimer to denatured monomer, the conformational stabilities, delta Gu (H2O), range from 10 to 27 kcal/mol, which is significantly greater than the conformational stability found for monomeric proteins. The relative contribution of quaternary interactions to the overall stability of the dimer can be estimated by comparing delta Gu (H2O) from equilibrium denaturation studies to the free energy associated with simple dissociation in the absence of denaturant. In many cases the large stabilization energy of dimers is primarily due to the intersubunit interactions and thus gives a rationale for the formation of oligomers. The magnitude of the conformational stability is related to the size of the polypeptide in the subunit and depends upon the type of structure in the subunit interface. The practical use, interpretation, and utility of estimation of conformational stability of dimers by equilibrium denaturation methods are discussed.     

Protein denaturation with guanidine hydrochloride or urea provides a different estimate of stability depending on the contributions of electrostatic interactions.

The objective of this study was to address the question of whether or not urea and guanidine hydrochloride (GdnHCl) give the same estimates of the stability of a particular protein. We previously suspected that the estimates of protein stability from GdnHCl and urea denaturation data might differ depending on the electrostatic interactions stabilizing the proteins. Therefore, 4 coiled-coil analogs were designed, where the number of intrachain and interchain electrostatic attractions (A) were systematically changed to repulsions (R): 20A, 15A5R, 10A10R, and 20R. The GdnHCl denaturation data showed that the 4 coiled-coil analogs, which had electrostatic interactions ranging from 20 attractions to 20 repulsions, had very similar [GdnHCl]1/2 values (average of congruent to 3.5 M) and, as well, their delta delta Gu values were very close to 0 (0.2 kcal/mol). In contrast, urea denaturation showed that the [urea]1/2 values proportionately decreased with the stepwise change from 20 electrostatic attractions to 20 repulsions (20A, 7.4 M; 15A5R, 5.4 M; 10A10R, 3.2 M; and 20R, 1.4 M), and the delta delta Gu values correspondingly increased with the increasing differences in electrostatic interactions (20A-15A5R, 1.5 kcal/mol; 20A-10A10R, 3.7 kcal/mol; and 20A-20R, 5.8 kcal/mol). These results indicate that the ionic nature of GdnHCl masks electrostatic interactions in these model proteins, a phenomenon that was absent when the unchanged urea was used. Thus, GdnHCl and urea denaturations may give vastly different estimates of protein stability, depending on how important electrostatic interactions are to the protein.     

Kinetic folding and unfolding of staphylococcal nuclease and its six mutants studied by stopped-flow circular dichroism

Kinetics of refolding and unfolding of staphylococcal nuclease and its six mutants, each carrying single or double amino acid substitutions, are studied by stopped-flow circular dichroism measurements. A transient kinetic intermediate formed within 10 ms after refolding starts possesses a substantial part of the N-domain core β-structure, whereas helices are formed at the later stages. The structure of the kinetic intermediate is less organized than the structure that is known to be formed by a nuclease 1-136 fragment. Only the refolding kinetics are affected by the mutations in all the mutants except two in which the mutations have changed the native structure. From this result and also from the locations of the mutation sites, the major N-terminal domain of the nuclease in the transition state of folding has a structure nearly identical to the native one. On the other hand, the minor C-terminal domain has previously been shown to be still disorganized in the transition state. The effects of the amino acid substitutions on the stability of the native and the transition states are in good agreement with the changes in the hydration free energy, expected for the corresponding amino acid replacements in the unfolded polypeptide. Since side chains of all the mutated residues are not accessible to solvent in the native structure, the result suggests that it is the unfolded state that is mainly affected by the mutations. © 1995 Wiley-Liss, Inc.