Conformational plasticity in folding of the split beta-alpha-beta protein S6: evidence for burst-phase disruption of the native state

An increasing number of folding studies of two-state proteins shows that point mutations sometimes change the kinetic m-values, leading to kinks and curves in the chevron plots. The molecular origin of these changes is yet unclear although it is speculated that they are linked to structural rearrangement of the transition state or to accumulation of meta-stable intermediates. To shed more light on this issue, we present here a combined m and phi-value analysis of the split beta-alpha-beta protein S6. Wild-type S6 displays classical two-state kinetics with v-shaped chevron plot, but a majority of its mutants display distinct m-value changes or curved chevrons. We observe that this kinetic aberration of S6 is linked to mutations that are clustered in distinct regions of the native structure. The most pronounced changes, i.e. decrease in the m-value for the unfolding rate constant, are seen upon truncation of interactions between the N and C termini, whereas mutations in the centre of the hydrophobic core show smaller or even opposed effects. As a consequence, the calculated phi-values display a systematic increase upon addition of denaturant. In the case of S6, the phenomenon seems to arise from a general plasticity of the different species on the folding pathway. That is, the structure of the denatured ensemble, the transition state, and the native ground-state for unfolding seem to change upon mutation. From these changes, it is concluded that interactions spanning the centre of the hydrophobic core form early in folding, whereas the entropically disfavoured interactions linking the N and C termini consolidate very late, mainly on the down-hill-side of the folding barrier.     

Refolding rate of stability-enhanced cytochrome c is independent of thermodynamic driving force.

N52I iso-2 cytochrome c is a variant of yeast iso-2 cytochrome c in which asparagine substitutes for isoleucine 52 in an alpha helical segment composed of residues 49-56. The N52I substitution results in a significant increase in both stability and cooperativity of equilibrium unfolding, and acts as a "global suppressor" of destabilizing mutations. The equilibrium m-value for denaturant-induced unfolding of N52I iso-2 increases by 30%, a surprisingly large amount for a single residue substitution. The folding/unfolding kinetics for N52I iso-2 have been measured by stopped-flow mixing and by manual mixing, and are compared to the kinetics of folding/unfolding of wild-type protein, iso-2 cytochrome c. The results show that the observable folding rate and the guanidine hydrochloride dependence of the folding rate are the same for iso-2 and N52I iso-2, despite the greater thermodynamic stability of N52I iso-2. Thus, there is no linear free-energy relationship between mutation-induced changes in stability and observable refolding rates. However, for N52I iso-2 the unfolding rate is slower and the guanidine hydrochloride dependence of the unfolding rate is smaller than for iso-2. The differences in the denaturant dependence of the unfolding rates suggest that the N52I substitution decreases the change in the solvent accessible hydrophobic surface between the native state and the transition state. Two aspects of the results are inconsistent with a two-state folding/unfolding mechanism and imply the presence of folding intermediates: (1) observable refolding rate constants calculated from the two-state mechanism by combining equilibrium data and unfolding rate measurements deviate from the observed refolding rate constants; (2) kinetically unresolved signal changes ("burst phase") are observed for both N52I iso-2 and iso-2 refolding. The "burst phase" amplitude is larger for N52I iso-2 than for iso-2, suggesting that the intermediates formed during the "burst phase" are stabilized by the N52I substitution.     

BPPred: a Web-based computational tool for predicting biophysical parameters of proteins

We exploit the availability of recent experimental data on a variety of proteins to develop a Web-based prediction algorithm (BPPred) to calculate several biophysical parameters commonly used to describe the folding process. These parameters include the equilibrium m-values, the length of proteins, and the changes upon unfolding in the solvent-accessible surface area, in the heat capacity, and in the radius of gyration. We also show that the knowledge of any one of these quantities allows an estimate of the others to be obtained, and describe the confidence limits with which these estimations can be made. Furthermore, we discuss how the kinetic m-values, or the Beta Tanford values, may provide an estimate of the solvent-accessible surface area and the radius of gyration of the transition state for protein folding. Taken together, these results suggest that BPPred should represent a valuable tool for interpreting experimental measurements, as well as the results of molecular dynamics simulations.     

The origin of pH-dependent changes in m-values for the denaturant-induced unfolding of proteins

Denaturant-induced unfolding is one of the most prevalent means of evaluating the structural stability of proteins and of determining the energetic consequences of mutations or changes in solution conditions. In spite of the widespread use of this approach, controversies and inconsistencies still persist with regard to the interpretation of the results of such studies. For example, most proteins show either a significant increase or a decrease (as much as 100 %) in the denaturant-dependence of the free energy of unfolding (i.e. the m-value) under increasingly acidic conditions. The pH dependence of the m-value is given different interpretations depending on whether the m-values increase or decrease with decreasing pH. In cases where m-values decrease, the decrease is attributed to the presence of an intermediate that becomes transiently stabilized during the unfolding transition at low pH. Cases where m-values increase as pH is lowered are usually interpreted in terms of an increase in the amount of surface area exposed by the denatured state at low pH. We have developed a general thermodynamic model that accounts for both types of behavior in terms of an intermediate that is populated throughout the unfolding transition. The model provides a unified framework for explaining both types of observed behavior, and the validity of the model was tested through the analysis of the pH dependence of m-values of staphylococcal nuclease. According to the model, the observed increase in m-values with decreasing pH is consistent with the existence of an intermediate that is populated during urea and guanidine unfolding. The intermediate becomes less populated during the unfolding transition at lower pH values giving rise to the apparent increase in m-values. These results argue that the prevailing interpretation need not apply to all proteins.     

Chemical denaturation: potential impact of undetected intermediates in the free energy of unfolding and m-values obtained from a two-state assumption.

The chemical unfolding transition of a protein was simulated, including the presence of an intermediate (I) in equilibrium with the native (N) and unfolded (U) states. The calculations included free energies of unfolding, DeltaGuw, in the range of 1.4 kcal/mol to 10 kcal/mol and three different global m-values. The simulations included a broad range of equilibrium constants for the N left arrow over right arrow I process. The dependence of the N <--> I equilibrium on the concentration of denaturant was also included in the simulations. Apparent DeltaGuw and m-values were obtained from the simulated unfolding transitions by fitting the data to a two-state unfolding process. The potential errors were calculated for two typical experimental situations: 1) the unfolding is monitored by a physical property that does not distinguish between native and intermediate states (case I), and 2) the physical property does not distinguish between intermediate and unfolded states (case II). The results obtained indicated that in the presence of an intermediate, and in both experimental situations, the free energy of unfolding and the m-values could be largely underestimated. The errors in DeltaGuw and m-values do not depend on the m-values that characterize the global N <--> U transition. They are dependent on the equilibrium constant for the N <--> I transition and its characteristic m1-value. The extent of the underestimation increases for higher energies of unfolding. Including no random error in the simulations, it was estimated that the underestimation in DeltaGuw could range between 25% and 35% for unfolding transitions of 3-10 kcal/mol (case I). In case II, the underestimation in DeltaGuw could be even larger than in case I. In the same energy range, a 50% error in the m-value could also take place. The fact that most of the mutant proteins are characterized by both a lower m-value and a lower stability than the wild-type protein suggests that in some cases the results could have been underestimated due to the application of the two-state assumption.     

Formation of the molten globule-like state of cytochrome c induced by n-alkyl sulfates at low concentrations

The molten globule state of cytochrome c is the major intermediate of protein folding. Elucidation of the thermodynamic mechanism of conformational stability of the molten globule state would enhance our understanding of protein folding. The formation of the molten globule state of cytochrome c was induced by n-alkyl sulfates including sodium octyl sulfate, SOS; sodium decyl sulfate, SDeS; sodium dodecyl sulfate, SDS; and sodium tetradecyl sulfate, STS, at low concentrations. The refolding states of the protein were monitored by spectroscopic techniques including circular dichroism (CD), visible absorbance and fluorescence. The effect of n-alkyl sulfates on the structure of acid-unfolded horse cytochrome c at pH 2 was utilized to investigate the contribution of hydrophobic interactions to the stability of the molten globule state. The addition of n-alkyl sulfates to the unfolded state of cytochrome c appears to support the stabilized form of the molten globule. The m-values of the refolded state of cytochrome c by SOS, SDeS, SDS, and STS showed substantial variation. The enhancement of m-values as the stability criterion of the molten globule state corresponded with increasing chain length of the cited n-alkyl sulfates. The compaction of the molten globule state induced by SDS, as a prototype for other n-alkyl sulfates, relative to the unfolded state of cytochrome c was confirmed by Stokes radius and thermal transition point (T(m)) measured by microviscometry and differential scanning calorimetry (DSC), respectively. Thus, hydrophobic interactions play an important role in stabilizing the molten globule state.     

A mass spectrometry-based probe of equilibrium intermediates in protein-folding reactions

Described here is a mass spectrometry- and H/D exchange-based approach for the detection of equilibrium intermediate state(s) in protein-folding reactions. The approach utilizes the stability of unpurified proteins from rates of H/D exchange (SUPREX) technique to measure the m value (i.e., delta DeltaG/delta [denaturant] value) associated with the folding reaction of a protein. Such SUPREX m-value analyses can be made over a wide range of denaturant concentrations. Thus, the described approach is well-suited for the detection of high-energy intermediates that might be populated at low denaturant concentrations and hard to detect in conventional chemical denaturation experiments using spectroscopic probes. The approach is demonstrated on four known non-two-state folding proteins, including alpha-lactalbumin, cytochrome c, intestinal fatty acid binding protein (IFABP), and myoglobin. The non-two-state folding behavior of each model protein system was detected by the described method. The cytochrome c, myoglobin, and IFABP systems each had high-energy intermediate states that were undetected in conventional optical spectroscopy-based studies and previously required other more specialized biophysical approaches (e.g., nuclear magnetic resonance spectroscopy-based methods and protease protection assays) for their detection. The SUPREX-based approach outlined here offers an attractive alternative to these other approaches, because it has the advantage of speed and the ability to analyze both purified and unpurified protein samples in either concentrated or dilute solution.     

Kinetics and mechanism of the folding of cytochrome c.

The reversible folding of cytochrome c in urea at pH 4.0 was investigated by repetitive pressure perturbation kinetics and by equilibrium spectroscopic methods. Two folding reactions were observed in the 1 ms to 10 s time range. The rates and amplitudes of these reactions depend on urea concentration in a complex manner, which is different for each process. The absorbance spectra of the kinetic amplitudes of the two reactions also differ from each other. A model with a three-state mechanism can quantitatively account for all of the kinetic and equilibrium data, and it enables us to determine the rate constants and volume changes of the two steps. If a rapid protonation step is added to the mechanism, the analysis can be extended to calculate the pH dependence of the rate and amplitude of the faster folding step. This pH dependence is in excellent agreement with previously published data [Tsong, T. Y. (1977) J. Biol. Chem. 252, 8778-8780]. Kinetic experiments in the 695-nm band show clearly that the axial ligand methionine-80 is involved in the slow folding process and the other axial ligand, histidine-18, is involved in the fast process. Additional experiments with a cyanogen bromide fragment of the protein, and fluorescence detection of the folding kinetics of the intact protein, support an interpretation of the model in terms of known structural elements of cytochrome c. This work provides new information about the mechanism of the folding of cytochrome c, resolves conflicts in earlier interpretations, and demonstrates the applicability of the repetitive pressure perturbation kinetics method to protein folding.     

Thermal denaturation of iso-1-cytochrome c variants: comparison with solvent denaturation.

Thermal denaturation studies as a function of pH were carried out on wild-type iso-1-cytochrome c and three variants of this protein at the solvent-exposed position 73 of the sequence. By examining the enthalpy and Tm at various pH values, the heat capacity increment (delta Cp), which is dominated by the degree of change in nonpolar hydration upon protein unfolding, was found for the wild type where lysine 73 is normally present and for three variants. For the Trp 73 variant, the delta Cp value (1.15 +/- 0.17 kcal/mol K) decreased slightly relative to wild-type iso-1-cytochrome c (1.40 +/- 0.06 kcal/mol K), while for the Ile 73 (1.65 +/- 0.07 kcal/mol K) and the Val 73 (1.50 +/- 0.06 kcal/mol K) variants, delta Cp increased slightly. In previous studies, the Trp 73, Ile 73, and Val 73 variants have been shown to have decreased m-values in guanidine hydrochloride denaturations relative to the wild-type protein (Hermann L, Bowler BE, Dong A, Caughey WS. 1995. The effects of hydrophilic to hydrophobic surface mutations on the denatured state of iso-1-cytochrome c: Investigation of aliphatic residues. Biochemistry 34:3040-3047). Both the m-value and delta Cp are related to the change in solvent exposure upon unfolding and other investigators have shown a correlation exists between these two parameters. However, for this subset of variants of iso-1-cytochrome c, a lack of correlation exists which implies that there may be basic differences between the guanidine hydrochloride and thermal denaturations of this protein. Spectroscopic data are consistent with different denatured states for thermal and guanidine hydrochloride unfolding. The different response of m-values and delta Cp for these variants will be discussed in this context.     

Origin of unusual phi-values in protein folding: evidence against specific nucleation sites

phi(f)-value analysis is one of the most common methods to characterize the structure of protein folding transition states. It compares the effects of mutations on the folding kinetics with the respective effects on equilibrium stability. The interpretation of the results usually focuses on a few unusual phi(f)-values, which are either particularly high or which are larger than 1 or smaller than 0. These mutations are believed to affect the most important regions for the folding process. A major uncertainty in experimental phi(f)-values is introduced by the commonly used analysis of only a single mutant at various positions in a protein (two-point analysis). To test the reliability of two-point phi(f)-values we used reference data from three positions in two different proteins at which multiple mutations have been introduced. The results show that two-point phi(f)-values are highly inaccurate if the difference in stability between two variants is less than 7 kJ/mol, corresponding to a 20-fold difference in equilibrium constant. Comparison with reported phi(f)-values for 11 proteins shows that most unusual phi(f)-values are observed in mutants which show changes in protein stability that are too small to allow a reliable analysis. The results argue against specific nucleation sites in protein folding and give a picture of transition states as distorted native states for the major part of a protein or for large substructures.     

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.     

Analysis of the stability of multimeric proteins by effective DeltaG and effective m-values

Analyzing the stability of a multimeric protein is challenging because of the intrinsic difficulty in handling the mathematical model for the folded multimer-unfolded monomer equilibrium. To circumvent this problem, we introduce the concept of effective stability, DeltaGeff (= -RTlnKeff), where Keff is the equilibrium constant expressed in monomer units. Analysis of the denaturant effect on DeltaGeff gives new insight into the stability of multimeric proteins. When a multimeric protein is mostly folded, the dependence of effective stability on denaturant concentration (effective m-value) is simply the m-value of its monomeric unit. However, when the protein is mostly unfolded, its stability depends on denaturant concentration with the m-value of its multimeric form. We also find that the effective m-value at the Cm is a good approximation of the apparent m-value determined by fitting the equilibrium unfolding data from multimeric proteins with a two-state monomer model. Moreover, when the m-value of a monomeric unit is estimated from its size, the effective stability of a multimeric protein can be determined simply from Cm and this estimated m-value. These simple and intuitive approaches will allow a facile analysis of the stability of multimeric proteins. These analyses are also applicable for high-throughput analysis of protein stability on a proteomic scale.     

Metrics that differentiate the origins of osmolyte effects on protein stability: a test of the surface tension proposal

Osmolytes that are naturally selected to protect organisms against environmental stresses are known to confer stability to proteins via preferential exclusion from protein surfaces. Solvophobicity, surface tension, excluded volume, water structure changes and electrostatic repulsion are all examples of forces proposed to account for preferential exclusion and the ramifications exclusion has on protein properties. What has been lacking is a systematic way of determining which force(s) is(are) responsible for osmolyte effects. Here, we propose the use of two experimental metrics for assessing the abilities of various proposed forces to account for osmolyte-mediated effects on protein properties. Metric 1 requires prediction of the experimentally determined ability of the osmolyte to bring about folding/unfolding resulting from the application of the force in question (i.e. prediction of the m-value of the protein in osmolyte). Metric 2 requires prediction of the experimentally determined ability of the osmolyte to contract or expand the Stokes radius of the denatured state resulting from the application of the force. These metrics are applied to test separate claims that solvophobicity/solvophilicity and surface tension are driving forces for osmolyte-induced effects on protein stability. The results show clearly that solvophobic/solvophilic forces readily account for protein stability and denatured state dimensional effects, while surface tension alone fails to do so. The agreement between experimental and predicted m-values involves both positive and negative m-values for three different proteins, and as many as six different osmolytes, illustrating that the tests are robust and discriminating. The ability of the two metrics to distinguish which forces account for the effects of osmolytes on protein properties and which do not, provides a powerful means of investigating the origins of osmolyte-protein effects.     

Mixed osmolytes: the degree to which one osmolyte affects the protein stabilizing ability of another

Mixtures of organic osmolytes occur in cells of many organisms, raising the question of whether their actions on protein stability are independent or synergistic. To investigate this question it is desirable to develop a system that permits evaluation of the effect of one osmolyte on the efficacy of another to either force-fold or denature a protein. A means of evaluating the efficacy of an osmolyte is provided by its m-value, an experimental quantity that measures the ability of the osmolyte to force a protein to unfold or fold. An experimental system is presented that enables evaluations of the m-values of osmolytes in the presence and absence of a second osmolyte. The experimental system involves use of a marginally stable protein in 10 mM buffer (pH 7, 200 mM salt, and 34 degrees C) that is at the midpoint of its native to denatured transition. These conditions enable determination of m-values for protecting and denaturing osmolytes in the presence and absence of a second osmolyte, permitting assessment of the extent to which the two osmolytes affect each other's efficacy. The two osmolytes investigated in this work are the denaturing osmolyte, urea, and the protecting osmolyte, sarcosine. Results show unequivocally that neither osmolyte alters the efficacy of the other in forcing the protein to fold or unfold-the osmolytes act independently on the protein despite their combined concentrations being in the multi-molar range. These osmolytes avoid altering one another's efficacy at these high concentrations because the number of osmolyte interaction sites on the protein is large and the binding constants are quite small. Consequently, the site occupancies are low enough in number that the two osmolytes neither compete nor cooperate in interacting with the protein.     

Rapid unfolding of a domain populates an aggregation-prone intermediate that can be recognized by GroEL

Some amino acid substitutions in phage P22 coat protein cause a temperature-sensitive folding (tsf) phenotype. In vivo, these tsf amino acid substitutions cause coat protein to aggregate and form intracellular inclusion bodies when folded at high temperatures, but at low temperatures the proteins fold properly. Here the effects of tsf amino acid substitutions on folding and unfolding kinetics and the stability of coat protein in vitro have been investigated to determine how the substitutions change the ability of coat protein to fold properly. The equilibrium unfolding transitions of the tsf variants were best fit to a three-state model, N if I if U, where all species concerned were monomeric, a result confirmed by velocity sedimentation analytical ultracentrifugation. The primary effect of the tsf amino acid substitutions on the equilibrium unfolding pathway was to decrease the stability (DeltaG) and the solvent accessibility (m-value) of the N if I transition. The kinetics of folding and unfolding of the tsf coat proteins were investigated using tryptophan fluorescence and circular dichroism (CD) at 222 nm. The tsf amino acid substitutions increased the rate of unfolding by 8-14-fold, with little effect on the rate of folding, when monitored by tryptophan fluorescence. In contrast, when folding or unfolding reactions were monitored by CD, the reactions were too fast to be observed. The tsf coat proteins are natural substrates for the molecular chaperones, GroEL/S. When native tsf coat protein monomers were incubated with GroEL, they bound efficiently, indicating that a folding intermediate was significantly populated even without denaturant. Thus, the tsf coat proteins aggregate in vivo because of an increased propensity to populate this unfolding intermediate.     

Stability and cooperativity of individual tertiary contacts in RNA revealed through chemical denaturation

For proteins, understanding tertiary interactions involved in local versus global unfolding has become increasingly important for understanding the nature of the native state ensemble, the mechanisms of unfolding, and the stability of both the native and intermediate states in folding. In this work we have addressed related questions with respect to RNA structure by combining chemical denaturation and hydroxyl radical footprinting methods. We have determined unfolding isotherms for each of 26 discrete sites of protection located throughout the Tetrahymena thermophila group I ribozyme. The cooperativity of folding, m-value, and the free energy, DeltaG degrees N-U, associated with formation of each tertiary contact was determined by analysis of the isotherms. The DeltaG degrees N-U values measured in this study vary from 1.7 +/- 0.2 to 7. 6 +/- 1.2 kcal mol-1. Thus, the stability of these discrete tertiary contacts vary by almost 104. In addition, an intradomain contact and three interdomain contacts show high cooperativity (m-values of 1.1 +/- 0.2 to 1.7 +/- 0.3 kcal mol-1 M-1) indicating that these contacts exhibit global cooperatively in their folding behavior. This new approach to examining RNA stability provides an exciting comparison to our understanding of protein structure and folding mechanisms.     

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

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

Fast folding of cytochrome c.

Native iso-2 cytochrome c contains two residues (His 18, Met 80) coordinated to the covalently attached heme. On unfolding of iso-2, the His 18 ligand remains coordinated to the heme iron, whereas Met 80 is displaced by a non-native heme ligand, His 33 or His 39. To test whether non-native His-heme ligation slows folding, we have constructed a double mutant protein in which the non-native ligands are replaced by asparagine and lysine, respectively (H33N,H39K iso-2). The double mutant protein, which cannot form non-native histidine-heme coordinate bonds, folds significantly faster than normal iso-2 cytochrome c: gamma = 14-26 ms for H33N,H39K iso-2 versus gamma = 200-1,100 ms for iso-2. These results with iso-2 cytochrome c strongly support the hypothesis that non-native His-heme ligation results in a kinetic barrier to fast folding of cytochrome c. Assuming that the maximum rate of a conformational search is about 10(11) s-1, the results imply that the direct folding pathway of iso-2 involves passage through on the order of 10(9) or fewer partially folded conformers.     

Communication of stabilizing energy between substructures of a protein

Thermodynamic communication between protein substructures has been investigated by determining the stabilizing effect of mutations at position 52 in the least stable, N-yellow, substructure of cytochrome c on the second least stable, Red, and most stable, Blue, substructures of the protein. A Lys 73 --> His (H73) variant of iso-1-cytochrome c, containing these mutations was used to measure the stability of the Red substructure of cytochrome c through the pH and guanidine hydrochloride (gdnHCl) dependence of the His 73-mediated alkaline conformational transition. The stability of the Blue substructure was measured by global unfolding with gdnHCl and increased by 1 to 3.5 kcal/mol versus the H73 variant. The data demonstrate that the increase in stability of the Red substructure is similar to the increase in global stability, consistent with upward propagation of stabilizing energy from less (N-yellow) to more stable (Red and Blue) protein substructures. The result also supports sequential rather than independent unfolding of the N-yellow and Red substructures of cytochrome c. The data indicate that a leucine at position 52 alters the nature of partial unfolding of the Red substructure, a surprising effect for a single-site mutation. For all variants, the thermodynamics of formation of the Lys 79 alkaline state, which does not unfold the entire Red substructure, shows less stabilization of the portion of the protein unfolded relative to the stabilization of the Blue substructure, indicating that propagation of energy between substructures is somewhat disrupted when unfolding does not correspond to a natural substructure.     

Stop-flow kinetics studies of the interaction of surfactant, sodium dodecyl sulfate, with acid-denatured cytochrome c

Interactions of sodium dodecyl sulfate (SDS) at submicellar and micellar concentration, with the globular protein, horse heart cytochrome c, at low pH have been shown to stabilize two molten globule-like intermediates. These dynamic studies were performed using far-UV, near-UV, and Soret-band circular dichroism (CD) as well as fluorescence methods. Stopped-flow CD and fluorescence studies of acid-denatured cytochrome c refolding with SDS were performed at both submicellar and micellar concentrations. Distinctive refolding mechanisms (from analysis of both CD and fluorescence) were found under these two conditions, and an obvious refolding intermediate was evident in the fluorescence traces. In addition, stopped-flow CD in the Soret region showed multistep kinetics, suggesting that the spectral changes in this region are not only solvent effect related but also connected with the change of secondary structure. A possible folding mechanism is proposed to rationalize the kinetics results.