“Designing Out” Disulfide Bonds: Thermodynamic Properties of 30−51 Cystine Substitution Mutants of Bovine Pancreatic Trypsin Inhibitor

Abstract
We have used a combination of spectroscopic and calorimetric techniques to assess the thermodynamic and extrathermodynamic consequences of paired amino acid substitutions at positions 30 and 51 in bovine pancreatic trypsin inhibitor (BPTI). Correctly folded, wild type BPTI contains a disulfide at the 30−51 positions, with the nonbackbone atoms of this cystine being relatively solvent inaccessible. Mutants missing this buried 30−51 disulfide adopt a conformation very similar to that of the native state of wild type BPTI (Eigenbrot et al., 1990, 1992), although they are severely destabilized relative to the wild type molecule (Hurle et al., 1990). We have conducted a systematic effort to find the energetically most favorable substitution for this buried 30−51 disulfide in BPTI. To this end, we have studied and characterized the thermally induced and guanidine hydrochloride-induced denaturation transitions for a family of mutants in which the amino acid residue(s) at positions 30 and/or 51 have been systematically altered. Specifically, we studied the unfolding transitions of the following series of residue 30/residue 51 paired substitution mutants: C30A/C51A, C30V/C51A, C30G/C51A, C30S/C51A, C30T/C51A, C30A/C51S, C30S/C51S, and C30G/C51M. For this series of mutants, comparisons between the relative stabilization free energies, derived from analysis of the denaturation profiles, allow us to reach the following conclusions: (a) side chains containing polar moieties (Ser and Thr) are destabilizing, with this effect being position dependent (i.e., a serine substitution is more destabilizing at position 51 than at position 30); (b) the destabilizing effects of two serine substitutions are approximately additive, suggesting that side chain−side chain hydrogen bonds between the two serine hydroxyl groups probably are weak or nonexistent; (c) the thermodynamic impact of a Gly30 substitution is consistent with a glycine-induced increase in the configurational entropy of the unfolded state; (d) the C30G/C51M mutant is highly destabilized relative to the C30A/C51A mutant despite the fact that, based on considerations of hydrophobicity and steric fit, substitution of a buried disulfide by Gly30 and Met51 would be expected to be optimal. Methionine may be particularly ill-suited as a buried disulfide substitute due to the large loss of side chain conformational entropy it undergoes during the transition from the unfolded to the native state. In the aggregate, our data provide insight into the residue-, position-, and context-dependent consequences on protein stability of “designing out” the buried 30−51 disulfide bond in the BPTI molecule. These data also suggest that a previously unrecognized component of disulfide bridge stabilization of proteins is the relatively minor penalty in side chain conformational entropy incurred by cystine residues during folding due to their severely restricted rotation even in the unfolded state.