| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Facultad de Ciencias, Departamento de Quimica Fisica, Universidad de Granada, 18071 Granada, Spain
Correspondence: Address reprint requests to Dr. Beatriz Ibarra-Molero or to Dr. Jose M. Sanchez-Ruiz, Facultad de Ciencias, Departamento de Quimica Fisica, Universidad de Granada, 18071 Granada, Spain. Tel.: 34-958-243-189; Fax: 34-958-272-879; E-mail: beatriz{at}ugr.es; sanchezr{at}ugr.es.
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONCONCLUDING REMARKSAPPENDIX 1: ON THE...APPENDIX 2: ON THE...ACKNOWLEDGEMENTSREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONCONCLUDING REMARKSAPPENDIX 1: ON THE...APPENDIX 2: ON THE...ACKNOWLEDGEMENTSREFERENCES |
|---|
; Sundd et al., 2002; Laurents et al., 2003) and that the surface-charge distribution may be rationally designed for enhanced protein stability and for optimized intermolecular interactions (Grimsley et al., 1999; Ibarra-Molero et al., 1999a; Loladze et al., 1999; Pace, 2000; Perl et al., 2000; Spector et al., 2000; Lee and Tidor, 2001; Nohaile et al., 2001; Perl and Schmid, 2001; Sanchez-Ruiz and Makhatadze, 2001; Ibarra-Molero and Sanchez-Ruiz, 2002; Marshall et al., 2002; Martin et al., 2002; Makhatadze et al., 2003). The experimental estimation of the contribution from charge-charge interactions to processes involving proteins is, therefore, an issue of considerable importance. This estimation is often carried out on the basis of the sodium chloride dependence of the energetic parameters for the process under study, under the assumption that charge-charge interactions are effectively screened out at moderate NaCl concentrations (below 1 M). However, evidence has accumulated over the years indicating that some surface salt bridges can be insensitive to NaCl, even when there is a clear charge-charge, coulombic interaction (Perutz et al., 1985; Yu et al., 1996; Kao et al., 2000; Luisi et al., 2003). Also, the recent work of Dominy et al. (2002) supports that NaCl is likely to screen more efficiently long-range interactions over short-range ones. It appears, therefore, that complete screening of all charge-charge interactions in proteins by moderate concentrations of a given salt (such as NaCl) cannot, in general, be taken for granted.
Here we address a different, but related, issue: the extent to which the charge screening efficiency depends on the nature of the salt. It must be noted that several analyses of experimental data support that denaturants tend to interact preferentially with protein surfaces, whereas stabilizers tend to be preferentially excluded (Arakawa and Timasheff, 1984, 1985; Makhatadze and Privalov, 1992; Bolen and Baskarov, 2001; Courtenay et al., 2001). It appears likely then that denaturing salts (such as guanidinium chloride) accumulate near the surface of proteins and, as a result, that they are very efficient at screening charge-charge interactions in proteins. Indeed, for several proteins, the guanidinium chloride dependence of the denaturation free energy shows abrupt deviations from linearity at low denaturant concentrations (<
1 M), which have been attributed to the screening of interactions involving charged groups (Santoro and Bolen, 1992; Monera et al., 1994; Ibarra-Molero and Sanchez-Ruiz, 1996; Ibarra-Molero et al., 1999a; Bolen and Yang, 2000; Garcia-Mira and Sanchez-Ruiz, 2001) (in some cases, deviations have been interpreted in terms of ion binding (Greene and Pace, 1974; Santoro and Bolen, 1988; Pace et al., 1990; Hagihara et al., 1993; Mayr and Schmid, 1993; Makhatadze et al., 1998)).
We report here a detailed experimental characterization of the sodium chloride and guanidinium chloride effects on the thermodynamic stability of Escherichia coli thioredoxin and a T14K variant (designed for improved charge-charge interactions on the surface and slightly more stable than the wild-type (WT) form). Our results provide evidence for significant screening by the denaturing salt but not by NaCl and suggest that the efficiency of different salts to screen charge-charge interactions in proteins does correlate with their denaturing strength and, ultimately, with the position of the constituent ions in the Hofmeister series (Baldwin, 1996). We derive additional support for this proposal from experiments on the effect of guanidinium thiocyanate (a stronger denaturant than guanidinium chloride) on the stability of E. coli thioredoxin.
We use in this work a simple electrostatic model (of the Tanford-Kirkwood type) to derive rough theoretical estimates of the contribution from charge-charge interactions to thioredoxin stability and, also, to design the T14K variant. We wish to emphasize, however, that, to estimate the screening efficiency of the different salts, we use an experimental approach based on the premise that screening of charge interactions is reflected in a clearly nonlinear (exponential-like) dependence of denaturation free energy in the 0–1 M salt concentration range. Experimental detection of such "exponential-like" dependence is straightforward in the case of nondenaturing salts, such as NaCl. In the case of guanidinium chloride and guanidinium thiocyanate, on the other hand, the screening contribution to denaturation 
G is superimposed on a large decrease of G with salt concentration associated with the denaturing character of these salts; detection of the screening effect in these cases can be more conveniently carried out on the basis of m values (the derivatives -
G/[salt]) derived from experimental differential scanning calorimetry (DSC) data, as we showed a few years ago (Ibarra-Molero and Sanchez-Ruiz, 1996).
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONCONCLUDING REMARKSAPPENDIX 1: ON THE...APPENDIX 2: ON THE...ACKNOWLEDGEMENTSREFERENCES |
|---|
). Calculations on variants of thioredoxin, in which neutral-residue
{Glu or Lys} mutations (see lower panel in Fig. 1) were introduced, were based on structures modeled (starting from the WT-thioredoxin structure) with the SwissViewer v3.7b2 program. No structure optimization was performed, but calculations were carried out for all the sterically allowed rotamers of the new side chains.
|
), which randomly places 2000 points in the expanded van der Waals sphere representing each atom. A radius of 1.4 Å for the solvent probe and the Chothia set (Chothia, 1976) for the protein atoms were used. Residue accessibilities were calculated as the ratio between the side-chain ASA in the native structure and that in Gly-X-Gly tripeptide.
Site-directed mutagenesis
Oligonucleotides used for mutagenesis were obtained from Genotek (Ottawa, Canada). Mutation in the codon corresponding to position 14 in the amino acid sequence of thioredoxin was introduced by the QuikChange site-directed mutagenesis method developed by Stratagene (La Jolla, CA). Briefly, the QuikChange method is based on polymerase chain reaction amplification using two complementary oligonucleotide primers containing the desired mutation. The parental nonmutated DNA is finally digested by an endonuclease. Mutation was verified by DNA sequence analysis.
Protein expression and purification
Plasmid pTK100 encoding wild-type thioredoxin (a gift from Dr. Maria Luisa Tasayco) was transformed into E. coli JF521 strain for protein overexpression. Cells were grown, starting from single colonies, at 37°C in Luria broth with 40 µg/mL of kanamycin to select for the plasmid-bearing cells. The final 750-fold dilution of the cell broth was allowed to grow during 12 h after stationary phase was reached. After centrifugation, cell pellets were frozen at -20°C until purification. Protein purification protocol was as follows. Briefly, cells were thawed, resuspended in 1 mM EDTA, 30 mM TRIZMA buffer, pH 8.3, and lyzed using a French press. The cell debris was centrifuged and the supernatant was collected and stirred with streptomycin sulfate (10% w/v) at 4°C overnight to precipitate nucleic acids. The filtered supernatant was then loaded onto a 2 L Sephacryl S-100 high resolution (Amersham Pharmacia Biotech, Upsala, Sweden) gel filtration column equilibrated in 1 mM EDTA, 30 mM TRIZMA buffer, pH 8.3. Thioredoxin fractions were identified by SDS-PAGE, pooled, and applied to a 250 mL Fractogel EMD DEAE (M) (Merck, Darmstadt, Germany) ion exchange column equilibrated in 1 mM EDTA, 30 mM TRIZMA buffer, pH 8.3. The protein was eluted by a linear gradient between 0 and 0.5 M NaCl. The proteins were pure as measured by SDS-PAGE gel densitometry. The molecular weight of pure proteins was confirmed by mass spectrometry. Thioredoxin concentration was determined spectrophotometrically at 280 nm using a published value of the extinction coefficient (Holmgren and Reichard, 1967).
Reagents and experimental conditions
Guanidinium chloride was ultrapure grade from Pierce (Rockford, IL). NaCl was analysis grade from Merck. Deionized water was used throughout. Aqueous stock solutions of WT and T14K thioredoxin were prepared by exhaustive dialysis against 5 mM HEPES, pH 7.0. Stock solutions of 6 M guanidinium chloride in HEPES buffer and 6 M NaCl in HEPES buffer were prepared as described previously (Ibarra-Molero et al., 1999b). Guanidinium chloride concentrations were determined from refraction index measurements (Pace et al., 1989) using an Atago (Tokyo, Japan) R 5000 hand refractometer. Guanidinium thiocyanate was ultrapure grade from Sigma (St. Louis, MO) and its solutions were prepared by weight.
Measurements of pH values for guanidinium salts solutions were carried out after calibration of the glass-electrode cell with aqueous standard buffers. No pH corrections (Garcia-Mira and Sanchez-Ruiz, 2001; Acevedo et al., 2002) were applied and, therefore, the pH value of 7 given for the guanidinium salts solutions is actually a pH-meter reading value or "apparent" pH value (see Garcia-Mira and Sanchez-Ruiz, 2001, for further discussion). We note that, in any case, the pH-dependence of thioredoxin denaturation energetics appears to be rather small, in the neighborhood of pH 7 (see Georgescu et al., 2001).
Circular dichroism
The experiments were carried out using a Jasco (Tokyo, Japan) J-715 spectropolarimeter equipped with a PTC-348WI temperature control unit.
Two sets of thermal unfolding experiments were collected to study guanidinium chloride concentration effects in the far and near ultraviolet (UV) regions, respectively. The change in circular dichroism (CD) signal within the temperature range of 20–85°C was monitored at both 222 and 280 nm, using an equilibration time of 60 s, signal-averaging time of 4 s, a bandwidth of 2 nm, and 1.0 nm step size. Thermal scans were done at 0, 0.5, 1, and 1.6 M guanidinium chloride. In the far-UV CD experiments, the protein concentration was 0.2 mg/mL and 1 mm cell path length was used. In the near-UV CD thermal melts, the protein concentration was 0.9 mg/mL and 10 mm cell path length was used. All transitions were highly reversible as was shown by the recovery of the CD signal after cooling the protein solution from 85°C to 20°C. Fittings of the two-state model to the CD signal versus temperature profiles were carried out assuming linear pre- and posttransition baselines and that the denaturation enthalpy does not significantly change within the narrow temperature range of the transition (see Ibarra-Molero and Sanchez-Ruiz, 1997, for further details).
Additionally, far-UV CD spectra (from 260 to 210 nm) and near-UV CD spectra (from 325 to 260 nm) of wild-type and T14K thioredoxin were monitored in the absence and in the presence of guanidinium chloride at 20°C and 85°C, using a bandwidth of 1 nm, an average of 4 scans, and 1.0 nm step size.
Differential scanning calorimetry
DSC experiments were carried out with a VP-DSC calorimeter from MicroCal (Northampton, MA) at a scan rate of 1.5 K/min. Protein solutions for the calorimetric experiments were prepared by exhaustive dialysis against the buffer (5 mM HEPES, pH 7.0). The samples were degassed at room temperature before the calorimetric experiments. Calorimetric cells (operating volume 0.5 ml) were kept under an excess pressure of 30 psi to prevent degassing during the scan. In all measurements, the buffer from the last dialysis step was used in the reference cell of the calorimeter. Several buffer-buffer baselines were obtained before each run with a protein solution to ascertain proper equilibration of the instrument. In most experiments, a reheating run was carried out to determine the reversibility of the denaturation process. Finally, an additional buffer-buffer baseline was obtained immediately after the protein runs to check that no significant change in instrumental baseline had occurred. When working with aqueous solutions, the level of instrumental baseline reproducibility attained was excellent and similar to that we have recently described (see Fig. 2 in Irun et al., 2001). However, as we have pointed out (Plaza del Pino and Sanchez-Ruiz, 1995; Ibarra-Molero et al., 1999b) baseline reproducibility is significantly poorer in the presence of cosolvents. This prevents us from obtaining absolute heat capacity values in NaCl and guanidinium salt solutions, although it does not compromise the calculation of denaturation enthalpies and denaturation temperatures from the analysis of the transitions.
|
). Therefore, we carried out all the DSC experiments at comparatively low protein concentrations: 0.5 mg/mL or below in some cases. We found no protein concentration effects on denaturation energetics within the 0.1–0.5 mg/mL range. Fittings of theoretical models to the heat capacity profiles were performed using programs written by us in the MLAB environment (Civilized Software, Silver Spring, MD). The general approach used in the two-state fittings was as described previously (Ibarra-Molero et al., 1999b).
Gibbs energy calculations
Protein stability curves (plots of denaturation G versus temperature) for aqueous solutions (i.e., in the absence of NaCl or guanidinium salts) were calculated from DSC data under two different assumptions: I), assuming that the denaturation heat capacity change is temperature independent and using the constant-CP integrated Gibbs-Helmholtz equation with the CP value obtained by extrapolating to the transition Tm the pre- and posttransition baselines; and II), taking into account the potential temperature dependence of CP and assuming, for the purposes of CP(T) calculation, that the heat capacity of the denatured state is equal to that expected for a fully solvated unfolded state as estimated from the group contributions given by Makhatadze and Privalov (1990) (the native-state and unfolded-state heat capacities we used are those shown in Fig. 5 of Georgescu et al., 2001). Further details about the calculation of protein stability under different assumptions for the temperature dependence of CP can be found elsewhere (Ibarra-Molero and Sanchez-Ruiz, 1996; Ibarra-Molero et al., 1999b). Here, we simply point out that the procedures used (I and II) involve different assumptions and that the comparison of the two stability curves obtained provides an estimate of the uncertainty involved in the G calculation from DSC data in this case (see Results for further details).
|
K mutation and of NaCl on denaturation G were very small and, therefore, we could use the Schellman equation (Schellman, 1987) to calculate them without introducing significant errors. In the terminology of this work, the Schellman equation can be written as,
 | (1) |
G is the perturbation Gibbs energy, 
and 
are the denaturation temperature and denaturation enthalpy change (at the denaturation temperature) for the unperturbed protein, and Tm is the perturbation effect on denaturation temperature. Note that Eq. 1 provides the G value at the temperature . If the perturbation is the T14K mutation, the unperturbed protein is the WT form, and Tm is the difference Tm(T14K) - Tm(WT). If the perturbation is the presence of NaCl, the unperturbed protein is the WT form in the absence of salt and Tm is the difference Tm(WT in the presence of salt) - Tm(WT in the absence of salt); then application of Eq. 1 allows us to obtain G = G(WT in the presence of NaCl) - G(WT in the absence of salt); the calculation actually yields G(WT in the presence of NaCl), since at the denaturation temperature of the unperturbed protein G(WT in the absence of salt) = 0. The NaCl concentration dependence of G for WT thioredoxin denaturation was obtained in this way.
Monte Carlo estimates of the errors associated to the reported energetic parameters
All error intervals given in this work have been obtained using the Monte Carlo method. That is, several replicas of each given original data set were randomly generated using suitable distribution functions for the errors associated to the original data; subsequently, the replica data sets were processed in the same manner as the original set and the statistical analysis of the results obtained led to the errors associated to the derived energetic parameters. As an illustrative example, we explain below the Monte Carlo calculation of the errors for the guanidinium chloride m1/2 values.
The original data set for the m value calculation (see Discussion for details) consists of Tm and Hm values for different guanidinium chloride concentrations (the experimental C values). The experimental Tm versus C and Hm versus C dependencies could be adequately described by first- and third-order polynomials, respectively. We took those polynomials as the starting point for replica generation. That is, we used the polynomials to calculate, for the experimental C values, "error-free" Tm and Hm values. To these we added errors randomly generated according to Gaussian distributions of zero mean and standard deviation of 0.36° (for Tm) and 10 kJ/mol (for Hm) (these representative standard deviation values were obtained from the analysis of several DSC experiments carried out in the absence of denaturant). In this way, we generated 20 replicas of the original data, which were subjected to the same type of data processing: fitting of third-order polynomial to the Tm versus C dependence, from which we obtained the derivatives dTm/dC, which were used, together with the generated values of Tm and Hm, to obtain the m1/2 values. This procedure yielded 20 m1/2 values for each experimental denaturant concentration so that standard deviations could be calculated; these are given as the errors associated to the original m1/2 values.
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONCONCLUDING REMARKSAPPENDIX 1: ON THE...APPENDIX 2: ON THE...ACKNOWLEDGEMENTSREFERENCES |
|---|
) of the Tanford-Kirkwood model (Tanford and Kirkwood, 1957) (see Materials and Methods for details). The plot shown is similar to those we have previously reported for ubiquitin and other proteins (Ibarra-Molero et al., 1999a; Sanchez-Ruiz and Makhatadze, 2001): a positive value of the interaction energy for a given group means that the group is involved in predominantly destabilizing interactions with groups of alike charge; conversely, negative values for the interaction energy indicate stabilizing interactions with groups of opposite charge. In the case of the ubiquitin molecule, we found (Ibarra-Molero et al., 1999a) several groups with positive interaction energy and, in fact, charge-deletion and charge-reversal mutations on the corresponding positions led to the expected stability enhancements (Loladze et al., 1999). According to the calculations we report here (Fig. 1, upper panel), most ionizable groups in the thioredoxin molecule are involved in clearly stabilizing interactions at pH 7. In fact, it appears that, in this case, the charge distribution is already optimized for stability to a significant extent and the total charge-charge interaction energy (sum of all pairwise interaction energies) for native thioredoxin is estimated to be -22 kJ/mol, that is, significant when compared with typical denaturation Gibbs energy values and stabilizing. Of course, it is conceivable that electrostatic interactions may also occur in the denatured state (Pace et al., 2000; Guzman-Casado et al., 2003) and contribute to the denaturation Gibbs energy change. Nevertheless, it seems at least reasonable to assume in this case that denatured-state charge-charge interactions do not fully cancel the strongly stabilizing native-state interactions. Accordingly, we may expect a significant (and stabilizing) contribution from charge-charge interactions to the thermodynamic stability of thioredoxin.
The design of the T14K variant
Our electrostatic calculations on native thioredoxin (Fig. 1, upper panel) do not suggest any positions in which charge-deletion or charge-reversal mutations could likely lead to very large stability enhancement. We thus turned to consider the creation of additional favorable charge-charge interactions via the introduction of new charged groups (see Sanchez-Ruiz and Makhatadze, 2001, for a general discussion). We carried out electrostatic calculations on modeled variants of thioredoxin (see Materials and Methods for details) in which Glu or Lys residues had been substituted for surface polar residues (taken here as the polar residues with side-chain solvent accessibility >0.5). Calculations were performed for all sterically allowed rotamers of the new side chains, and the results are summarized in the lower panel of Fig. 1 as the difference between the total charge-charge interaction energy calculated for the variant and that corresponding to the WT. Again, all values for this difference are moderate, although the more promising variant appears to be T14K, which shows a calculated charge-charge interaction energy in the native state 2.5 kJ/mol lower than that for the WT form. This variant has been obtained and characterized in this work and it is, in fact, slightly more stable than the WT form (see further below). It is to be noted that the newly introduced Lys residue in this variant is fully exposed to the solvent: modeling based on the WT structure shows that all its rotamers are sterically allowed and our ASA calculations indicate an accessibility to the solvent of 0.81 ± 0.04 (average of the values obtained for all rotamers) when taking Gly-X-Gly tripeptides as reference (see Materials and Methods for details). Thus, comparison of the stability of the variant T14K with that of WT thioredoxin provides a suitable model system to test the efficiency of salts to screen a well-exposed charge.
Circular dichroism studies
Far-UV and near-UV CD spectra of the WT form of thioredoxin and the T14K variant under conditions in which both proteins are in the native state (pH 7, 20°C, guanidinium chloride concentrations within the 0–1.6 M range) are shown in Fig. 2. It appears that both the T14K mutation and the guanidinium chloride concentration have little effect on these native-state CD spectra.
Due to technical limitations, we could not carry out CD spectra determinations at temperatures >90°C. The spectra for the denatured states of WT and T14K shown in Fig. 2 were obtained at 85°C and in the presence of 1 M and 1.6 M guanidinium chloride (so that the denaturation temperature is clearly <70°C; see further below). The denatured-state far-UV CD spectra are similar for WT and T14K and suggestive of some kind of residual structure in the denatured state. Actually, there appears to be only a small difference between the native-state and the denatured-state far-UV CD spectra in the 210–260 nm region. On the other hand, the near-UV signal is essentially absent in the denatured state for both proteins (Fig. 2) and the ellipticity at 280 nm provides a sensitive probe to follow thermal denaturation, as is shown in Fig. 3.
|
), reproducibility of the instrumental baseline recorded with guanidinium chloride solutions in the calorimetric cells is poor, a fact that prevents us from calculating protein absolute heat capacities values from the DSC thermograms obtained in the presence of this denaturant. On the other hand, baseline reproducibility is excellent for aqueous buffers (see Materials and Methods) and we found no significant effect of the T14K mutation on the absolute heat capacities of the native and thermally denatured states (results not shown), although the variant is slightly more stable than the WT form as shown by a somewhat higher value of the denaturation temperature (Fig. 4, lower panel). It must be noted that the difference in denaturation temperature between T14K and WT thioredoxin is small, but significant and reproducible: we carried out three sets of DSC experiments with the two protein forms and found the three following values for the denaturation temperature difference (Tm = Tm(T14K) - Tm(WT)): 1.0 K, 1.0 K, and 1.1 K.
|
).
The effect of guanidinium thiocyanate on the thermal denaturation of WT thioredoxin as followed by DSC
We have carried out DSC experiments for WT thioredoxin at pH 7 and in the presence of several guanidinium thiocyanate concentrations within the 0–0.8 M range. It must be noted that, for guanidinium thiocyanate concentrations >1 M, the denaturation temperature values get close to room temperature, the denaturation enthalpies become small (which implies very broad transitions) and no useful information can be derived from the DSC thermograms—hence, the comparatively narrow range of guanidinium thiocyanate concentration studied.
The statements we have made in the previous section regarding baseline reproducibility, reversibility, and the two-state fittings also apply here. The values for the denaturation temperature and the denaturation enthalpy of WT thioredoxin in the presence of guanidinium thiocyanate are given in Fig. 5.
The effect of sodium chloride on the thermal denaturation of WT thioredoxin and the T14K variant as followed by DSC
We have carried out DSC experiments for WT thioredoxin and the T14K variant at pH 7 and in the presence of several sodium chloride concentrations within the 0–2 M range (see Fig. 4 for representative examples), and fittings of the DSC profiles were carried out on the basis of the two-state model. Only a very small effect of NaCl on the DSC profiles was found. The denaturation temperature values show a rather small (and essentially linear) increase with NaCl concentration (see upper panel of Fig. 5) and the denaturation enthalpy values do not change significantly within the 0–1 M NaCl concentration range, although they appear to decrease slightly with NaCl concentration >1 M (results not shown).
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONCONCLUDING REMARKSAPPENDIX 1: ON THE...APPENDIX 2: ON THE...ACKNOWLEDGEMENTSREFERENCES |
|---|
G versus denaturant concentration, known as the m value (Greene and Pace, 1974; Myers et al., 1995), has found widespread application in protein folding studies. Although m values are routinely determined from chemical denaturation experiments, we showed a few years ago (Ibarra-Molero and Sanchez-Ruiz, 1996) that they can also be obtained from DSC experiments performed at different denaturant concentrations using a procedure that is straightforward and, to a large extent, model-independent. Since these "calorimetrically determined" m values play a key role in the discussion of the results reported in this work, we briefly summarize in this section some essential features about their calculation and interpretation.
We take all thermodynamic changes for an equilibrium denaturation process as functions of both temperature (T) and denaturant concentration (C); thus, the denaturation change in Gibbs energy is expressed as
 | (2) |
The values of T and C for which G = 0 (and, consequently, the equilibrium constant for the process is unity) define an equilibrium line in C versus T or T versus C plots. Note that this equilibrium line can be viewed in two entirely equivalent ways: 1), as the effect of denaturant concentration on the denaturation temperature (effect of C on Tm); and 2), as the effect of temperature on the mid-point denaturant concentration (effect of T on C1/2).
The partial derivative of G with respect to C gives, by definition, the denaturant m value,
 | (3) |
G versus C for that temperature will be linear down to C = 0 and the linear extrapolation method will be valid. Note, however, that we do not make here any assumptions regarding the denaturant concentration dependence of the m values (actually, we aim at determining such dependence from the experimental DSC data).
We will refer to the m values corresponding to the equilibrium line (i.e., to C and T conditions for which G = 0) as m1/2 values. These m1/2 values can be calculated from experimental DSC profiles using the following equation (for a derivation, see Ibarra-Molero and Sanchez-Ruiz, 1996),
 | (4) |
Equation 4 is rigorous (for two-state equilibrium denaturation) and, we note again, is not based on the linear extrapolation approximation. The calculation of m1/2 from Eq. 4 only requires values of the denaturation enthalpy change at the Tm (equilibrium-line Hm values) and the effect of denaturant concentration on denaturation temperature (so that the derivative dTm/dC can be computed). In particular, the value of the denaturation heat capacity is not required.
Equation 4 can be used to calculate m1/2 values for different denaturant concentrations, provided that DSC experiments at those denaturant concentrations have been performed (of course, the m1/2 values belong to the equilibrium line and thus they correspond to different denaturant concentrations and to different temperatures; we expect, however, the effect of denaturant concentration to dominate the change of m1/2 along the equilibrium line: see further below in this Discussion for an illustration). For hen egg-white lysozyme at pH 4.5 (see Fig. 6 in Ibarra-Molero and Sanchez-Ruiz, 1996) we found that the m1/2 values were constant within the experimental scatter for guanidinium chloride concentrations >1 M; <1 M, however, the m1/2 values increased sharply as the denaturant concentration approached zero. These results indicate that the guanidinium chloride dependence of G for lysozyme denaturation at pH 4.5 is actually linear over an extended denaturant concentration range, but also that a strong deviation from linearity occurs at low guanidinium chloride concentrations, deviation which is likely due to the screening of charge-charge interactions (for a more detailed discussion, see Ibarra-Molero and Sanchez-Ruiz, 1996, and Ibarra-Molero et al., 1999a). As we discuss below, the same general kind of behavior is found for thioredoxin.
|
Hm and Tm data of Fig. 5 (the calculation of the derivative dTm/dC was based on a polynomial fitting to the Tm versus C dependence: see legend to Fig. 5 for details). Clearly, for both proteins, m1/2 changes along the equilibrium line, in particular for low denaturant concentrations (<1 M). It is to be noted (see Fig. 6) that the m1/2 values at high C1/2 agree with the value reported in the literature from chemical denaturation studies (Kelley et al., 1987; Santoro and Bolen, 1992). The m1/2 values for both proteins appear to be the same within the experimental scatter and can be adequately described (see Fig. 6, upper panel) by the following empirical equation:
 | (5) |
= 13.6 ± 1.2 kJ·mol-1·M-1, ß = 13.7 ± 3.6 kJ·mol-1·M-1, and 
= 2.3 ± 0.6 M-1. It is important to note again that C1/2 changes with temperature (see Fig. 5, upper panel) and, therefore, there could be some doubt as to whether Eq. 5 actually reflects the guanidinium chloride dependence of the m1/2 values or, rather, an effect of temperature on them. This later interpretation, however, is disfavored by the following illustrative calculations (see Appendix 1 for a more formal and rigorous analysis):
Let us assume for the sake of the argument that Eq. 5 reflects exclusively a temperature effect and, consequently, that the m values are denaturant-concentration independent for any given temperature. Then, the linear extrapolation method will be valid and denaturation Gibbs energy at zero denaturant concentration can be calculated (Ibarra-Molero and Sanchez-Ruiz, 1996) as m1/2 x C1/2. That is:
 | (6) |
GLEM(T, C = 0) is the linear-extrapolation estimate of G at the temperature T and zero denaturant concentration, and C1/2(T) is the C1/2 value at the temperature T.
If, on the other hand, we assume that Eq. 5 represents the actual denaturant-concentration effect, then the temperature effect is not significant, m1/2 and C1/2 in Eq. 5 can be taken simply as m and C, and integration yields the denaturant-concentration dependence of G at any temperature:
 | (7) |
G at zero denaturant concentration:
 | (8) |
Both Eqs. 6 and 8 provide denaturation G values in the absence of denaturant; they, however, are based upon different assumptions. The validity of these assumptions may be assessed by comparing the G values calculated on the basis of Eqs. 6 and 8 with those obtained from the DSC profile in the absence of denaturant by using standard thermodynamic procedures (see Materials and Methods for details). Such a comparison is shown in the lower panel of Fig. 6 and supports clearly the validity of Eq. 8 and its underlying assumptions.
The guanidinium chloride concentration dependence of the denaturation Gibbs energy for thioredoxin
From the above calculations, we conclude that, to an acceptable degree of approximation, the dependence of m1/2 with C1/2 shown in the upper panel of Fig. 6 (and described by Eq. 5) reflects the actual denaturant-concentration dependence of the m values. Accordingly, we are justified in using Eq. 7 to calculate the denaturant-concentration dependence of G at constant temperature. The results of such calculation are given in Fig. 7 (upper panel) for the denaturation temperature of WT thioredoxin in the absence of denaturant (note that, for that temperature, C1/2 = 0). As was to be expected from the m1/2 data (Fig. 4, upper panel), the dependence of G with C is linear over an extended denaturant-concentration range, but that there is a clear deviation from linearity below C 1 M, a behavior similar to that we have previously found for other protein systems (Ibarra-Molero and Sanchez-Ruiz, 1996; Ibarra-Molero et al., 1999a).
|
 | (9) |
As shown in the lower panel of Fig. 5, the change of the denaturation enthalpy along the equilibrium line can be adequately described as a linear dependence with the denaturation temperature:
 | (10) |
= 7.0 kJ·K-1·mol-1 and TH = 300 K. TH can be interpreted as the temperature at which Hm becomes zero for a given denaturant concentration (CH). On the other hand, we do not interpret as a denaturation heat capacity change, since both temperature and guanidinium chloride concentration change along the equilibrium line, and cosolvents (such as denaturants) may affect the denaturation enthalpy value (for a clear example of this, see Plaza del Pino and Sanchez-Ruiz, 1995). Substitution of Eq. 10 into Eq. 9, variable separation and integration (from {TH,CH} to any point in the equilibrium line: {Tm,C}) yields:
 | (11) |
 | (12) |
For the sake of convenience, we will refer to the left-hand-side term in Eq. 12 as F(Tm). Since we know the values of and TH, we can calculate F(Tm) for the several experimental Tm values and construct the plot of F(Tm) versus C. According to Eq. 12, if the linear G versus C dependence holds, this plot must be linear with a slope equal to minus the m value. The plot of F(Tm) versus C for WT thioredoxin is shown in the middle panel of Fig. 7. The data points corresponding to the three highest guanidinium chloride concentrations describe a straight line with a slope of -14.6 kJ·K-1·mol-1, in excellent agreement with the high-C m value we have obtained in this work (upper panel of Fig. 6) and with the m values derived from chemical denaturation experiments reported in the literature (Kelley et al., 1987; Santoro and Bolen, 1992). Note that the linear extrapolation value of F(Tm) is 35.3, consistent with a denaturation temperature of 84°C, whereas the experimental Tm value in the absence of denaturant is 88.8°C. That is, the linear extrapolation underestimates the Tm value.
There can be little doubt from all the above (upper panel in Fig. 6 and upper and middle panels in Fig. 7) that the dependence of denaturation G with guanidinium chloride concentration shows a clear deviation from linearity <1 M denaturant, and this deviation is such that linear extrapolation from high denaturant concentration underestimates the stability of thioredoxin in the absence of denaturant (see upper and middle panels in Fig. 7). It is intriguing then that Santoro and Bolen in their 1992 work reported the opposite for thioredoxin under similar solvent conditions (pH 7): a low-C deviation from linearity in the G versus C dependence, but in such a way that linear extrapolation overestimated the stability of the protein at zero denaturant concentration. The reason for this discrepancy is not clear to us; perhaps, it is related to the fact that Santoro and Bolen (1992) had to rely on a rather long constant CP Gibbs-Helmholtz extrapolation to obtain the G values at 25°C and low-C, as well as the fact that they had to use an estimate of CP based on the temperature-dependence of H values obtained by changing solvent composition. We emphasize again that the m1/2 calculation given here (based on our work in Ibarra-Molero and Sanchez-Ruiz, 1996) does not rely on Gibbs-Helmholtz extrapolations and does not require the use of a CP value.
The deviation from the linear G versus C relation that we find in this work cannot be explained in terms of specific binding of the denaturant ions to the native protein, since this would lead to increased values for denaturation G in the absence of denaturant upon linear extrapolation from high denaturant concentration (see Mayr and Schmid, 1993, and Appendix 2), which is the opposite effect of what we observe in the upper panel of Fig. 7. On the other hand, it would seem that a decreased linear-extrapolation G value (as shown in Fig. 7) could in principle be explained by the denaturant-binding model:
 | (13) |
n is the difference in the number of binding sites between the denatured and the native states, and Kb is the binding constant. However, the denaturant-binding model is not consistent with our experimental m values for thioredoxin denaturation (upper panel in Fig. 6). Thus, using the definition of m (Eq. 3) and Eq. 13, it is straightforward to arrive at the denaturant-concentration dependence of the m value predicted by the binding model:
 | (14) |
Note that this equation predicts that m will approach zero upon increasing denaturant concentration, in disagreement with the behavior of the experimental m values, which show clear evidence of approaching an m value of 15 kJ·mol-1·M-1 at high denaturant concentration. In fact, Eq. 14 is unable to yield an acceptable fit to the experimental m versus C dependence (results not shown).
Of course, the above analyses do no exhaust all possible models of denaturant action. However, the failure of the above models to account for our experimental data does suggest that the deviation from the linear G versus C dependence at low guanidinium chloride concentration is most likely associated to the screening of charge-charge interactions, a suggestion further supported by the two following facts: 1), for several proteins, sharp changes in calorimetrically determined m values at low-C are observed for guanidinium chloride-induced denaturation, but not for urea-induced denaturation (Ibarra-Molero et al., 2004); and 2), for ubiquitin, the deviations from the linear G versus C dependence at low guanidinium chloride concentration were found to change sign with pH (Ibarra-Molero et al., 1999a) in the way expected for a charge-charge contribution to protein stability.
It must be noted, nevertheless, that the size of the G deviation here is -7 kJ/mol, significantly smaller than the Tanford-Kirkwood estimate of the total energy of charge-charge interactions in the native structure of WT thioredoxin (-22 kJ/mol; see also Fig. 1). Several reasons may be adduced to account for this discrepancy: i), The Tanford-Kirkwood model is indeed a very simple one and, perhaps, we should only expect qualitative or semiquantitative predictions from it. ii), Electrostatic interactions in denatured states may be significant (Pace et al., 2000; Guzman-Casado et al., 2003); thus, even if we accepted as exact the Tanford-Kirkwood value for the energy of charge-charge interactions in native thioredoxin, this value would only provide an upper limit (in absolute value) to the contribution of charge-charge interactions to the denaturation Gibbs energy. iii), Guanidinium chloride may not be able to screen out all charge-charge interactions. Actually, this possibility is supported by the analysis of the guanidinium thiocyanate effects that we describe further below in this Discussion.
The sodium chloride concentration dependence of the denaturation Gibbs energy for thioredoxin
Sodium chloride concentrations within the 0–2 M range have a very small effect on the denaturation enthalpy and the denaturation temperature for WT thioredoxin and the T14K variant (see Fig. 5), and calculation of the denaturation G versus sodium chloride concentration for a temperature equal to the denaturation temperature in the absence of salt is straightforward (see Materials and Methods). Such a profile for WT thioredoxin is shown in Fig. 7. There is almost no effect of NaCl on G and, in particular, there is little evidence of a sharp change below 1 M salt that could be associated to screening of charge-charge interactions. It appears clear then that the experimental NaCl dependence of thioredoxin stability at pH 7 does not provide evidence for significant screening of charge-charge interactions (see, however, Concluding Remarks).
The effect of the T14K mutation on thioredoxin stability
The effect of NaCl and guanidinium chloride on the G value for the T14K mutation (G = G(T14K) - G(WT)) can be calculated from the experimental Tm and denaturation enthalpy values using a straightforward procedure described in Materials and Methods. The results of such calculation are shown in the lower panel of Fig. 7. It appears that both salts are able to screen out most of the stabilization afforded by the T14K mutation and, in addition, that they are roughly equally effective in this regard, possibly reflecting the fact that the lysine at position 14 in the T14K variant is well exposed to the solvent.
It is perhaps worth noting here that the stabilization afforded by the T14K mutation in the absence of salts, as measured by G, is 1 kJ/mol, that is, less than half the value predicted by the Tanford-Kirkwood calculation (2.5 kJ/mol; see Fig. 1). This discrepancy is approximately by the same factor than that we described above between estimate of the total charge-charge contribution to denaturation G and the corresponding Tanford-Kirkwood prediction, and, probably, the same reasons (except, of course, reason iii) can be adduced to explain it (see above).
Finally, it is interesting that a G of 1 kJ/mol in the absence of salt translates into an increment in denaturation temperature (Tm) of 1° only. To a good degree of approximation, mutation effects on G and Tm are related through the Schellman equation (Schellman, 1987) (see Materials and Methods),
 | (15) |
H and S indicates that they correspond to the denaturation temperature. According to Eq. 15, the low Tm value is associated to comparatively high value for the denaturation enthalpy at the denaturation temperature (450 kJ/mol for thioredoxin). Indeed, it is well known that small proteins with low denaturation enthalpy values are more sensitive (in terms of denaturation temperature) to stabilizing effects of mutations and the environment, and often display high Tm values (Alexander et al., 1992; Ibarra-Molero et al., 2000). It is clear that the achievement of significant increments in Tm for not-so-small proteins via optimization of charge-charge interactions must rely on the cumulative effect of several mutations (for a discussion, see Sanchez-Ruiz and Makhatadze, 2001), that is, in the design of the surface charge distribution (see Ibarra-Molero and Sanchez-Ruiz, 2002).
Guanidinium thiocyanate m1/2 values for WT thioredoxin
As originally planned, this work was meant to consist in a detailed energetic characterization of the effect of sodium chloride and guanidinium chloride on thioredoxin denaturation addressed at determining the relative charge-screening efficiencies of these salts. However, since our results point to a relation between the screening efficiency and the denaturing strength, we deemed convenient to include some experimental data on the effect of guanidinium thiocyanate (a stronger denaturant than guanidinium chloride) on thioredoxin stability.
Guanidinium thiocyanate m1/2 values for WT thioredoxin denaturation were calculated in the same manner as the guanidinium chloride m1/2 values (although, in this case, a linear dependence sufficed to describe the denaturant-concentration dependence of Tm: see Fig. 5) and are shown in the upper panel of Fig. 6. Due to the small range of guanidinium thiocyanate concentration studied (see Results), we could not fully characterize the m versus C dependence. It is clear, however, that the m values for guanidinium thiocyanate are much larger than those for guanidinium chloride and that they decrease with concentration in a sharper manner, suggesting a higher charge-charge screening efficiency.
|
CONCLUDING REMARKS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONCONCLUDING REMARKSAPPENDIX 1: ON THE...APPENDIX 2: ON THE...ACKNOWLEDGEMENTSREFERENCES |
|---|
, and Dominy et al., 2002) often show an exponential-like dependence in the 0–1 M concentration range (attributed to screening of charge-charge interactions by the salt) and a gradual and almost linear dependence at higher concentrations. In the case of E. coli thioredoxin at pH 7, only a moderate linear increase in G is observed upon increasing NaCl concentration; that is, our experimental data do not provide evidence for significant NaCl screening of charge-charge interactions. The simplest explanation for this is, of course, that screening does not take place in this case (since we do not find evidence for it). However, there is an alternative explanation that is, at least, plausible. Dominy et al. (2002) have pointed out that NaCl is expected to screen more efficiently long-range interactions over short-range ones, and that the former mainly correlate with the total charge of the protein (Dominy et al., 2002; Zhou and Dong, 2003), which, for thioredoxin at pH 7, is -5 units. We might conceive then that unequal screening of short-range (predominantly stabilizing) interactions and long-range (predominantly destabilizing) interactions may perhaps yield a close-to-zero Gibbs energy balance in the low-NaCl concentration range, in such a way that screening would not be apparent in the NaCl dependence of the thermodynamic stability. Our results support a relation between screening and the accessibility to solvent of the charged groups. Thus, 1 M NaCl (as well as 1 M guanidinium chloride) does appear to screen significantly the interactions of a well-exposed Lys group in a T14K variant of thioredoxin. It must be noted that Lys-14 in this variant has an accessibility to solvent (0.81, taking a Gly-Lys-Gly peptide as reference) higher than that for the charged groups involved in strongly stabilizing interactions in WT thioredoxin (the accessibilities of K3, K18, K36, K57, K69, K82, K90, K96, and K100 in the native structure of WT thioredoxin are 0.34,
