Equilibrium and kinetic measurements of the conformational transition of thioredoxin in urea

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

Natural selection for kinetic stability is a likely origin of correlations between mutational effects on protein energetics and frequencies of amino acid occurrences in sequence alignments

It appears plausible that natural selection constrains, to some extent at least, the stability in many natural proteins. If, during protein evolution, stability fluctuates within a comparatively narrow range, then mutations are expected to be fixed with frequencies that reflect mutational effects on stability. Indeed, we recently reported a robust correlation between the effect of 27 conservative mutations on the thermodynamic stability (unfolding free energy) of Escherichia coli thioredoxin and the frequencies of residues occurrences in sequence alignments. We show here that this correlation likely implies a lower limit to thermodynamic stability of only a few kJ/mol below the unfolding free energy of the wild-type (WT) protein. We suggest, therefore, that the correlation does not reflect natural selection of thermodynamic stability by itself, but of some other factor which is linked to thermodynamic stability for the mutations under study. We propose that this other factor is the kinetic stability of thioredoxin in vivo, since( i) kinetic stability relates to irreversible denaturation, (ii) the rate of irreversible denaturation in a crowded cellular environment (or in a harsh extracellular environment) is probably determined by the rate of unfolding, and (iii) the half-life for unfolding changes in an exponential manner with activation free energy and, consequently, comparatively small free energy effects can have deleterious consequences for kinetic stability. This proposal is supported by the results of a kinetic study of the WT form and the 27 single-mutant variants of E. coli thioredoxin based on the global analyses of chevron plots and equilibrium unfolding profiles determined from double-jump unfolding assays. This kinetic study suggests, furthermore, one of the factors that may contribute to the high activation free energy for unfolding in thioredoxin (required for kinetic stability), namely the energetic optimization of native-state residue environments in regions, which become disrupted in the transition state for unfolding.

Thermodynamics of replacing an alpha-helical Pro residue in the P40S mutant of Escherichia coli thioredoxin

Escherichia coli thioredoxin is a 108 amino acid oxidoreductase and contains a single Met residue at position 37. The protein contains a long alpha-helical stretch between residues 32 and 49. The central residue of this helix, Pro40, has been replaced by Ser. The stabilities of the oxidized states of two proteins, the single mutant M37L and the double mutant M37L,P40S, have been characterized by differential scanning calorimetry (DSC) and also by a series of isothermal guanidine hydrochloride (GuHCl) melts in the temperature range of 277 to 333 K. The P40S mutation was found to stabilize the protein at all temperatures upto 340 K though both proteins had similar Tm values of about 356 K. At 298 K, the M37L,P40S mutant was found to be more stable than M37L by 1.5 kcal/mol. A combined analysis of GuHCl and calorimetric data was carried out to determine the enthalpy, entropy, and heat capacity change upon unfolding. At 298 K there was a large, stabilizing enthalpic effect in P40S though significant enthalpy-entropy compensation was observed and the two proteins had similar values of deltaCp. Thus, replacement of a Pro in the interior of an alpha helix can have substantial effects on protein stability.

A protein folding intermediate of ribonuclease T1 characterized at high resolution by 1D and 2D real-time NMR spectroscopy

The rate-limiting step during the refolding of S54G/P55N ribonuclease T1 is determined by the slow trans-->cis prolyl isomerisation of Pro39. We investigated the refolding of this variant by one-dimensional (1D) and two-dimensional (2D) real-time NMR spectroscopy, initiated by a tenfold dilution from 6 M guanidine hydrochloride at 10 degreesC. Two intermediates could be resolved with the 1D approach. The minor intermediate, which is only present early during refolding, is largely unfolded. The major intermediate, with an incorrect trans Pro39 peptide bond, is highly structured with 33 amide protons showing native chemical shifts and native NOE patterns. They could be assigned in a real-time 2D-NOESY (nuclear Overhauser enhancement spectroscopy) by using a new assignment strategy to generate positive and negative signal intensities for native and non-native NOE cross-peaks, respectively. Surprisingly, amide protons with non-native environments are located not only close to Tyr38-Pro39, but are spread throughout the entire protein, including the C-terminal part of the alpha-helix, beta-strands 3 and 4 and several loop regions. Native secondary and tertiary structure was found for the major intermediate in the N-terminal beta-strands 1 and 2 and the C terminus (connected by the disulfide bonds), the N-terminal part of the alpha-helix, and the loops between beta-strands 4/5 and 5/6. Implications of these native and non-native structure elements of the intermediate for the refolding of S54G/P55N ribonuclease T1 and for cis/trans isomerizations are discussed.

Unfolding of 4-aminobutyrate aminotransferase equilibrium and kinetic studies

The unfolding of pig liver 4-aminobutyrate aminotransferase by urea has been investigated at equilibrium. The overall process was reversible as judged from the recovery of catalytic activity after dilution of urea-treated samples. Unfolding of the enzyme was monitored by circular dichroism and fluorescence spectroscopy. The steepness of the fluorescence and CD changes between 2 and 8 M urea, and the lack of any discernible plateau suggests that unfolding of the protein is a cooperative process. The unfolding of 4-aminobutyrate aminotransferase as a function of urea concentration was monitored by fluorescence measurements of the tryptophanyl residues. The kinetic results indicate that the aminotransferase unfolds in a single kinetic phase. Unfolded 4-aminobutyrate aminotransferase recovers immediately its catalytic activity upon dilution with buffers of neutral pH. Based on equilibrium and kinetic results, a two-state model for the unfolding of the aminotransferase is proposed.

Improved synthetic methods for the selective deuteration of aromatic amino acids: applications of selective protonation towards the identification of protein folding intermediates through nuclear magnetic resonance

In this report we describe several novel methods for the preparation of selectively deuterated aromatic amino acids. New syntheses for [2,3,5,6-2H4]phenylalanine and [2,4,6,7-2H4]tryptophan, as well as improved catalytic exchange methods for [2,3,5,6-2H4]tyrosine and [2,3,4,5,6-2H5]phenylalanine are presented. Isotopic substitution levels for all compounds are generally found to be greater than 95%. Biosynthetic incorporation of these amino acids is also shown to be possible with little or no evidence of isotopic scrambling. The products from these new syntheses, in combination with other selectively deuterated aromatic amino acids, are found to permit group-specific 'single-proton' labelling of proteins. This highly-efficient and very cost-effective method of selective protonation is shown to produce greatly simplified 1H-NMR spectra of the aromatic region of proteins. The utility of this approach to isotopic editing is demonstrated with the identification of a transient folding intermediate of Escherichia coli thioredoxin which is undetectable by standard 2-D NMR techniques.

Use of analytical gel chromatography to analyze tertiary and quaternary structural changes in E. coli aspartate transcarbamylase

E. coli aspartate transcarbamylase (ATCase) is a large (310 kDa) protein that undergoes major changes in quaternary structure when substrates and regulatory nucleotides bind. We have used analytical gel chromatography to detect quaternary structure changes in both the holoenzyme and its catalytic subunit (c3), to characterize the quaternary structure of single site mutant proteins and to monitor urea-induced dissociation and unfolding of c3. Binding of the bisubstrate analog PALA (N-(phosphonacetyl)-L-aspartate) to ATCase and c3 has been shown to alter s20.w by -3.3% and + 1.4%, respectively [Howlett, G.J. and Schachman, H.K. (1977), Biochemistry 23, 5077-5083]. The corresponding changes in the chromatographic partition coefficient (sigma) are -2.6 +/- 0.3% and 5.5 +/- 1.9% on Sephacryl S400HR and S200, respectively. Partition coefficients of mutant ATCases with single site mutations in the c chain differ from those of the wild-type protein by +/- 0.5% in small zone experiments; for example, mutations Arg 269----Gly and Glu 239----Gln alter the partition coefficient by 0.4% and -0.5%, respectively. The partition coefficient of mutant Glu 50----Gln is identical to the wild type enzyme. In the presence of saturating PALA, partition coefficients of Glu 50----Gln and Arg 269----Gly, but not Glu 239----Gln are identical to those of the wild type. Results for Glu 239----Gln are consistent with measurements of activity, small angle X-ray scattering and sedimentation coefficient that indicate that mutations at this site shift the quaternary structure towards the R state [Ladjimi and Kantrowitz (1988), Biochemistry 27, 276-83; Vachette and Hervé, cited by Kantrowitz and Lipscomb (1988), Science 241, 669-674; Newell and Schachman (1988), FASEB J. 2, A551]. Results for Glu 50----Gln are also consistent with measurements of activity (Ladjimi et al. (1988), Biochemistry 27, 268-276). The changes in tertiary and quaternary structure that result from urea-induced denaturation of c3 result in larger changes in the partition coefficient. Dissociation into folded monomers in 1-1.75 M urea is accompanied by a 4.6% increase in partition coefficient, while denaturation at greater than 5 M urea gives rise to a 43% decrease on S-300 Sephacryl. The bisubstrate analog PALA suppresses dissociation and increases the cooperativity of the unfolding reaction.

Escherichia coli thioredoxin folds into two compact forms of different stability to urea denaturation

The urea-induced denaturation of Escherichia coli thioredoxin and thioredoxin variants has been examined by electrophoresis on urea gradient slab gels by the method of Creighton [Creighton, T. (1986) Methods Enzymol. 131, 156-172]. Thioredoxin has only two cysteine residues, and these form a redox-active disulfide at the active site. Oxidized thioredoxin-S2 and reduced thioredoxin-(SH)2 each show two folded isomers with a large difference in stability to urea denaturation. The difference in stability is greater for the isomers of oxidized than for the isomers of reduced thioredoxin. At 2 degrees C, the urea concentrations at the denaturation midpoint are approximately 8 and 4.3 M for the oxidized isomers and 4.8 and 3.7 M for the reduced isomers. The difference between the gel patterns of samples applied in native versus denaturing buffer, and at 2 and 25 degrees C, is characteristic for the involvement of a cis-proline-trans-proline isomerization. The data very strongly suggest that the two folded forms of different stabilities correspond to the cis and trans isomers of the highly conserved Pro 76 peptide bond, which is cis in the crystal structure of oxidized thioredoxin. Urea gel experiments with the mutant thioredoxin P76A, with alanine substituted for proline at position 76, corroborate this interpretation. The electrophoretic banding pattern diagnostic for an involvement of proline isomerization in urea denaturation is not observed for oxidized P76A. In broad estimates of delta G degree for the native-denatured transition, the difference in delta G degree (no urea) between the putative cis and trans isomers of the Ile 75-Pro 76 peptide bond is approximately 3 kcal/mol for oxidized thioredoxin and approximately 1.5 kcal/mol for reduced thioredoxin. Since cis oxidized thioredoxin is much more stable than trans, folded oxidized thioredoxin is essentially all cis. In folded reduced thioredoxin, cis and trans interconvert slowly, on the minute time scale at 2 and 25 degrees C. In the absence of urea, the folded reduced thioredoxin is less than a few percent trans. Three additional mutants with additions or substitutions at the active site also show electrophoresis banding patterns consistent with a difference in stability between cis and trans isomers.

Equilibrium and kinetic measurements of the conformational transition of reduced thioredoxin

The single disulfide bond in Escherichia coli thioredoxin was reduced by reaction with a 20-fold excess of reduced dithiothreitol at neutral pH and 25 degrees C. For some measurements, reduced thioredoxin was further reacted with iodoacetamide to alkylate the cysteinyl residues. The denaturation transitions of oxidized, reduced, and reduced alkylated thioredoxin were observed by using far-ultraviolet circular dichroic and exclusion chromatographic measurements. Cleavage of the disulfide bond lowers the stability of the native thioredoxin to denaturation by about 2.4 kcal/mol, and subsequent alkylation lowers the stability by a further 1.6 kcal/mol. The kinetics of the conformational change of reduced thioredoxin in guanidine hydrochloride were observed by using exclusion chromatography at moderate pressure and 2 degrees C. Analyses of single and multimixing protocols are consistent with a predominant nonnative configuration in the denatured state and the transient accumulation of a compact nativelike intermediate during refolding. The intermediate can incorporate the nonnative configuration and can accommodate its isomerization. No compelling chromatographic evidence was found for a conformation having an elution time different from that characteristic for either the native or the denatured protein.