pH-jump-induced folding and unfolding studies of barstar: evidence for multiple folding and unfolding pathways.

Equilibrium and kinetic characterization of the high pH-induced unfolding transition of the small protein barstar have been carried out in the pH range 7-12. A mutant form of barstar, containing a single tryptophan, Trp 53, completely buried in the core of the native protein, has been used. It is shown that the protein undergoes reversible unfolding above pH 10. The pH 12 form (the D form) appears to be as unfolded as the form unfolded by 6 M guanidine hydrochloride (GdnHCl) at pH 7 (the U form): both forms have similar fluorescence and far-UV circular dichroism (CD) signals and have similar sizes, as determined by dynamic light scattering and size-exclusion chromatography. No residual structure is detected in the D form: addition of GdnHCl does not alter its fluorescence and far-UV CD properties. The fluorescence signal of Trp 53 has been used to monitor folding and unfolding kinetics. The kinetics of folding of the D form in the pH range 7-11 are complex and are described by four exponential processes, as are the kinetics of unfolding of the native state (N state) in the pH range 10.5-12. Each kinetic phase of folding decreases in rate with increase in pH from 7 to 10.85, and each kinetic phase of unfolding decreases in rate with decrease in pH from 12 to 10.85. At pH 10.85, the folding and unfolding rates for any particular kinetic phase are identical and minimal. The two slowest phases of folding and unfolding have identical kinetics whether measured by Trp 53 fluorescence or by mean residue ellipticity at 222 nm. Direct determination of the increase in the N state with time of folding at pH 7 and of the D form with time of unfolding at pH 12, by means of double-jump assays, show that between 85 and 95% of protein molecules fold or unfold via fast pathways between the two forms. The remaining 5-15% of protein molecules appear to fold or unfold via slower pathways, on which at least two intermediates accumulate. The mechanism of folding from the high pH-denatured D form is remarkably similar to the mechanism of folding from the urea or GdnHCl-denatured U form.

GroEL channels the folding of thioredoxin along one kinetic route.

Many proteins display complex folding kinetics, which represent multiple parallel folding pathways emanating from multiple unfolded forms and converging to the unique native form. The small protein thioredoxin from Escherichia coli is one such protein. The effect of the chaperonin GroEL on modulating the complex energy landscape that separates the unfolded ensemble from the native state of thioredoxin has been studied. It is shown that while the fluorescence change accompanying folding occurs in five kinetic phases in the absence of GroEL, only the two slowest kinetic phases are discernible in the presence of saturating concentrations of GroEL. This result is shown to be consistent with only one out of several available folding routes being operational in the presence of GroEL. It is shown that native protein, which forms via fast as well as slow routes in the absence of GroEL, forms only via a slow route in its presence. The effect of GroEL on the folding of thioredoxin is shown to be the consequence of it binding differentially to the many folding-competent forms. While some of these forms can continue folding when bound to GroEL, others cannot. All molecules are then drawn into the operational folding route by the law of mass action. This observation indicates a new role for GroEL, which is to bias the energy landscape of a folding polypeptide towards fewer available pathways. It is suggested that such channeling might be a mechanism to avoid possible aggregation-prone routes available to a refolding polypeptide in vivo.

Folding of horse cytochrome c in the reduced state.

Equilibrium and kinetic folding studies of horse cytochrome c in the reduced state have been carried out under strictly anaerobic conditions at neutral pH, 10 degrees C, in the entire range of aqueous solubility of guanidinium hydrochloride (GdnHCl). Equilibrium unfolding transitions observed by Soret heme absorbance, excitation energy transfer from the lone tryptophan residue to the ferrous heme, and far-UV circular dichroism (CD) are all biphasic and superimposable, implying no accumulation of structural intermediates. The thermodynamic parameters obtained by two-state analysis of these transitions yielded DeltaG(H2O)=18.8(+/-1.45) kcal mol(-1), and C(m)=5.1(+/-0.15) M GdnHCl, indicating unusual stability of reduced cytochrome c. These results have been used in conjunction with the redox potential of native cytochrome c and the known stability of oxidized cytochrome c to estimate a value of -164 mV as the redox potential of the unfolded protein. Stopped-flow kinetics of folding and unfolding have been recorded by Soret heme absorbance, and tryptophan fluorescence as observables. The refolding kinetics are monophasic in the transition region, but become biphasic as moderate to strongly native-like conditions are approached. There also is a burst folding reaction unobservable in the stopped-flow time window. Analyses of the two observable rates and their amplitudes indicate that the faster of the two rates corresponds to apparent two-state folding (U<-->N) of 80-90 % of unfolded molecules with a time constant in the range 190-550 micros estimated by linear extrapolation and model calculations. The remaining 10-20 % of the population folds to an off-pathway intermediate, I, which is required to unfold first to the initial unfolded state, U, in order to refold correctly to the native state, N (I<-->U<-->N). The slower of the two observable rates, which has a positive slope in the linear functional dependence on the denaturant concentration indicating that an unfolding process under native-like conditions indeed exists, originates from the unfolding of I to U, which rate-limits the overall folding of these 10-20 % of molecules. Both fast and slow rates are independent of protein concentration and pH of the refolding milieu, suggesting that the off-pathway intermediate is not a protein aggregate or trapped by heme misligation. The nature or type of unfolded-state heme ligation does not interfere with refolding. Equilibrium pH titration of the unfolded state yielded coupled ionization of the two non-native histidine ligands, H26 and H33, with a pK(a) value of 5.85. A substantial fraction of the unfolded population persists as the six-coordinate form even at low pH, suggesting ligation of the two methionine residues, M65 and M80. These results have been used along with the known ligand-binding properties of unfolded cytochrome c to propose a model for heme ligation dynamics. In contrast to refolding kinetics, the unfolding kinetics of reduced cytochrome c recorded by observation of Soret absorbance and tryptophan fluorescence are all slow, simple, and single-exponential. In the presence of 6.8 M GdnHCl, the unfolding time constant is approximately 300(+/-125) ms. There is no burst unfolding reaction. Simulations of the observed folding-unfolding kinetics by numerical solutions of the rate equations corresponding to the three-state I<-->U<-->N scheme have yielded the microscopic rate constants.

Structure is lost incrementally during the unfolding of barstar.

Coincidental equilibrium unfolding transitions observed by multiple structural probes are taken to justify the modeling of protein unfolding as a two-state, N <==> U, cooperative process. However, for many of the large number of proteins that undergo apparently two-state equilibrium unfolding reactions, folding intermediates are detected in kinetic experiments. The small protein barstar is one such protein. Here the two-state model for equilibrium unfolding has been critically evaluated in barstar by estimating the intramolecular distance distribution by time-resolved fluorescence resonance energy transfer (TR-FRET) methods, in which fluorescence decay kinetics are analyzed by the maximum entropy method (MEM). Using a mutant form of barstar containing only Trp 53 as the fluorescence donor and a thionitrobenzoic acid moiety attached to Cys 82 as the fluorescence acceptor, the distance between the donor and acceptor has been shown to increase incrementally with increasing denaturant concentration. Although other probes, such as circular dichroism and fluorescence intensity, suggest that the labeled protein undergoes two-state equilibrium unfolding, the TR-FRET probe clearly indicates multistate equilibrium unfolding. Native protein expands progressively through a continuum of native-like forms that achieve the dimensions of a molten globule, whose heterogeneity increases with increasing denaturant concentration and which appears to be separated from the unfolded ensemble by a free energy barrier.

Reversible formation of on-pathway macroscopic aggregates during the folding of maltose binding protein.

Maltose binding protein (MBP) is widely used as a model for protein folding and export studies. We show here that macroscopic aggregates form transiently during the refolding of MBP at micromolar protein concentrations. Disaggregation occurs spontaneously without any aid, and the refolded material has structure and activity identical to those of the native, nondenatured protein. A considerable fraction of protein undergoing folding partitions into the aggregate phase and can be manually separated from the soluble phase by centrifugation. The separated MBP precipitate can be resolubilized and yields active, refolded protein. This demonstrates that both the soluble and aggregate phases contribute to the final yield of refolded protein. SecB, the cognate Escherichia coli cytosolic chaperone in vivo for MBP, reduces but does not entirely prevent aggregation, whereas GroEL and a variety of other control proteins have no effect. Kinetic studies using a variety of spectroscopic probes show that aggregation occurs through a collapsed intermediate with some secondary structure. The aggregate formed during refolding can convert directly to a near native state without going through the unfolded state. Further, optical and electron microscopic studies indicate that the MBP precipitate is not an amyloid.

Backbone dynamics of free barnase and its complex with barstar determined by 15N NMR relaxation study.

Backbone dynamics of uniformly 15N-labeled free barnase and its complex with unlabelled barstar have been studied at 40 degrees C, pH 6.6, using 15N relaxation data obtained from proton-detected 2D [1H]-15N NMR spectroscopy. 15N spin-lattice relaxation rate constants (R1), spin-spin relaxation rate constants (R2), and steady-state heteronuclear [1H]-15N NOEs have been measured at a magnetic field strength of 14.1 Tesla for 91 residues of free barnase and for 90 residues out of a total of 106 in the complex (excluding three prolines and the N-terminal residue) backbone amide 15N sites of barnase. The primary relaxation data for both the cases have been analyzed in the framework of the model-free formalism using both isotropic and axially symmetric models of the rotational diffusion tensor. As per the latter, the overall rotational correlation times (tau(m)) are 5.0 and 9.5 ns for the free and complexed barnase, respectively. The average order parameter is found to be 0.80 for free barnase and 0.86 for the complex. However, the changes are not uniform along the backbone and for about 5 residues near the binding interface there is actually a significant decrease in the order parameters on complex formation. These residues are not involved in the actual binding. For the residues where the order parameter increases, the magnitudes vary significantly. It is observed that the complex has much less internal mobility, compared to free barnase. From the changes in the order parameters, the entropic contribution of NH bond vector motion to the free energy of complex formation has been calculated. It is apparent that these motion's cause significant unfavorable contributions and therefore must be compensated by many other favorable contributions to effect tight complex formation. The observed variations in the motion and their different locations with regard to the binding interface may have important implications for remote effects and regulation of the enzyme action.

Backbone dynamics of barstar: a (15)N NMR relaxation study.

Backbone dynamics of uniformly (15)N-labeled barstar have been studied at 32 degrees C, pH 6.7, by using (15)N relaxation data obtained from proton-detected 2D (1)H-(15)N NMR spectroscopy. (15)N spin-lattice relaxation rate constants (R(1)), spin-spin relaxation rate constants (R(2)), and steady-state heteronuclear (1)H-(15)N NOEs have been determined for 69 of the 86 (excluding two prolines and the N-terminal residue) backbone amide (15)N at a magnetic field strength of 14.1 Tesla. The primary relaxation data have been analyzed by using the model-free formalism of molecular dynamics, using both isotropic and axially symmetric diffusion of the molecule, to determine the overall rotational correlation time (tau(m)), the generalized order parameter (S(2)), the effective correlation time for internal motions (tau(e)), and NH exchange broadening contributions (R(ex)) for each residue. As per the axially symmetric diffusion, the ratio of diffusion rates about the unique and perpendicular axes (D( parallel)/D( perpendicular)) is 0.82 +/- 0.03. The two results have only marginal differences. The relaxation data have also been used to map reduced spectral densities for the NH vectors of these residues at three frequencies: 0, omega(H), and omega(N), where omega(H),(N) are proton and nitrogen Larmor frequencies. The value of tau(m) obtained from model-free analysis of the relaxation data is 5.2 ns. The reduced spectral density analysis, however, yields a value of 5.7 ns. The tau(m) determined here is different from that calculated previously from time-resolved fluorescence data (4.1 ns). The order parameter ranges from 0.68 to 0.98, with an average value of 0.85 +/- 0.02. A comparison of the order parameters with the X-ray B-factors for the backbone nitrogens of wild-type barstar does not show any considerable correlation. Model-free analysis of the relaxation data for seven residues required the inclusion of an exchange broadening term, the magnitude of which ranges from 2 to 9.1 s(-1), indicating the presence of conformational averaging motions only for a small subset of residues.

The slow folding reaction of barstar: the core tryptophan region attains tight packing before substantial secondary and tertiary structure formation and final compaction of the polypeptide chain.

The slow folding of a single tryptophan-containing mutant of barstar has been studied in the presence of 2 M urea at 10 degrees C, using steady state and time-resolved fluorescence methods and far and near-UV CD measurements. The protein folds in two major phases: a fast phase, which is lost in the dead time of measurement during which the polypeptide collapses to a compact form, is followed by a slow observable phase. During the fast phase, the rotational correlation time of Trp53 increases from 2.2 ns to 7.2 ns, and its mean fluorescence lifetime increases from 2.3 ns to 3.4 ns. The fractional changes in steady-state fluorescence, far-UV CD, and near-UV CD signals, which are associated with the fast phase are, respectively, 36 %, 46 %, and 16 %. The product of the fast phase can bind the hydrophobic dye ANS. These observations together suggest that the folding intermediate accumulated at the end of the fast phase has: (a) about 20 % of the native-state secondary structure, (b) marginally formed or disordered tertiary structure, (c) a water-intruded and mobile protein interior; and (d) solvent-accessible patches of hydrophobic groups. Measurements of the anisotropy decay of Trp53 suggest that it undergoes two types of rotational motion in the intermediate: (i) fast (tau(r) approximately 1 ns) local motion of its indole side-chain, and (ii) a slower (tau(r) approximately 7.2 ns) motion corresponding to global tumbling of the entire protein molecule. The ability of the Trp53 side-chain to undergo fast local motion in the intermediate, but not in the fully folded protein where it is completely buried in the hydrophobic core, suggests that the core of the intermediate is still poorly packed. The global tumbling time of the fully folded protein is faster at 5.6 ns, suggesting that the volume of the intermediate is 25 % more than that of the fully folded protein. The rate of folding of this intermediate to the native state, measured by steady-state fluorescence, far-UV CD, and near-UV CD, is 0.07(+/-0.01) min(-1) This rate compares to a rate of folding of 0.03(+/-0.005) min(-1), determined by double-jump experiments which monitor directly formation of native protein; and to a rate of folding of 0.05 min(-1), when determined from time-resolved anisotropy measurements of the long rotational correlation time, which relaxes from an initial value of 7.2 ns to a final value of 5. 6 ns as the protein folds. On the other hand, the amplitude of the short correlation time decreases rapidly with a rate of 0.24(+/-0.06) min(-1). These results suggest that tight packing of residues in the hydrophobic core occurs relatively early during the observable slow folding reaction, before substantial secondary and tertiary structure formation and before final compaction of the protein.

A thermodynamic coupling mechanism can explain the GroEL-mediated acceleration of the folding of barstar.

Despite extensive structural and kinetic studies, the mechanism by which the Escherichia coli chaperonin GroEL assists protein folding has remained somewhat elusive. It appears that GroEL might play an active role in facilitating folding, in addition to its role in restricting protein aggregation by secluding folding intermediates. We have investigated the kinetic mechanism of GroEL-mediated refolding of the small protein barstar. GroEL accelerates the observed fast (millisecond) refolding rate, but it does not affect the slow refolding kinetics. A thermodynamic coupling mechanism, in which the concentration of exchange-competent states is increased by the law of mass action, can explain the enhancement of the fast refolding rates. It is not necessary to invoke a catalytic role for GroEL, whereby either the intrinsic refolding rate of a productive folding transition or the unfolding rate of a kinetically trapped off-pathway intermediate is increased by the chaperonin.

Measurements of cysteine reactivity during protein unfolding suggest the presence of competing pathways.

Evidence that proteins may unfold utilizing complex competing pathways comes from a new pulse-labeling protocol in which the change in reactivity of a single cysteine residue in a protein during unfolding is measured, making use of its easily monitored reaction with the Ellman reagent, dithionitrobenzoic acid. The kinetics of unfolding of two single cysteine-containing mutant forms of the small protein barstar, C82A, which contains only Cys40, and C40A, which contains only Cys82, have been studied. The data suggest that unfolding occurs via two parallel pathways, each forming competing intermediates. In one of these early intermediates, Cys40 and Cys82 are already as reactive as they are in the fully unfolded protein, while in the other intermediate, the Cys thiol groups are unreactive. One more long-lived intermediate also needs to be included on the pathway defined by the early intermediate with unreactive Cys thiol groups to account for the difference in the rates of fluorescence change and of change in Cys40 reactivity. The demonstration of multiple intermediates and pathways for unfolding indicates that protein unfolding reactions can be as complex as protein folding reactions.

Observation of multistate kinetics during the slow folding and unfolding of barstar.

The kinetics of the slow folding and unfolding reactions of barstar, a bacterial ribonuclease inhibitor protein, have been studied at 23(+/-1) degrees C, pH 8, by the use of tryptophan fluorescence, far-UV circular dichroism (CD), near-UV CD, and transient mixing (1)H nuclear magnetic resonance (NMR) spectroscopic measurements in the 0-4 M range of guanidine hydrochloride (GdnHCl) concentration. The denaturant dependences of the rates of folding and unfolding processes, and of the initial and final values of optical signals associated with these kinetic processes, have been determined for each of the four probes of measurement. Values determined for rates as well as amplitudes are shown to be very much probe dependent. Significant differences in the intensities and rates of appearance and disappearance of several resolved resonances in the real-time one-dimensional NMR spectra have been noted. The NMR spectra also show increasing dispersion of chemical shifts during the slow phase of refolding. The denaturant dependences of rates display characteristic folding chevrons with distinct rollovers under strongly native as well as strongly unfolding conditions. Analyses of the data and comparison of the results obtained with different probes of measurement appear to indicate the accumulation of a myriad of intermediates on parallel folding and unfolding pathways, and suggest the existence of an ensemble of transition states. The energetic stabilities of the intermediates estimated from kinetic data suggest that they are approximately half as stable as the fully folded protein. The slowness of the folding and unfolding processes (tau = 10-333 s) and values of 20.5 (+/-1.4) and 18 (+/-0.5) kcal mol(-)(1) for the activation energies of the slow refolding and unfolding reactions suggest that proline isomerization is involved in these reactions, and that the intermediates accumulate and are therefore detectable because the slow proline isomerization reaction serves as a kinetic trap during folding.

SecB binds only to a late native-like intermediate in the folding pathway of barstar and not to the unfolded state.

SecB is a cytosolic, tetrameric chaperone of Escherichia coli which maintains precursor proteins in a translocation competent state. We have investigated the effect of SecB on the refolding kinetics of the small protein barstar in 1 M guanidine hydrochloride at pH 7.0 and 25 degreesC using fluorescence spectroscopy. We show that SecB does not bind either the native or the unfolded states of barstar but binds to a late near-native intermediate along the folding pathway. For barstar, polypeptide collapse and formation of a hydrophobic surface are required for binding to SecB. SecB does not change the apparent rate constant of barstar refolding. The kinetic data for SecB binding to barstar are not consistent with simple kinetic partitioning models.

Stopped-flow NMR measurement of hydrogen exchange rates in reduced horse cytochrome c under strongly destabilizing conditions.

A procedure to measure exchange rates of fast exchanging protein amide hydrogens by time-resolved NMR spectroscopy following in situ initiation of the reaction by diluting a native protein solution into an exchanging deuterated buffer is described. The method has been used to measure exchange rates of a small set of amide hydrogens of reduced cytochrome c, maintained in a strictly anaerobic atmosphere, in the presence of an otherwise inaccessible range of guanidinium deuterochloride concentrations. The results for the measured protons indicate that hydrogen exchange in the unfolding transition region of cytochrome c reach the EX2 limit, but emphasize the difficulty in interpretation of the exchange mechanism in protein hydrogen exchange studies. Comparison of free energies of structure opening for the measured hydrogens with the global unfolding free energy monitored by far-UV CD measurements has indicated the presence of at least one partially unfolded equilibrium species of reduced cytochrome c. The results provide the first report of measurement of free energy of opening of structure to exchange in the 0-2-kcal/mol range.

Discrepancies between the NMR and X-ray structures of uncomplexed barstar: analysis suggests that packing densities of protein structures determined by NMR are unreliable.

The crystal structure of the C82A mutant of barstar, the intracellular inhibitor of the Bacillus amyloliquefaciens ribonuclease barnase, has been solved to a resolution of 2.8 A. The molecule crystallizes in the space group I41 with a dimer in the asymmetric unit. An identical barstar dimer is also found in the crystal structure of the barnase-barstar complex. This structure of uncomplexed barstar is compared to the structure of barstar bound to barnase and also to the structure of barstar solved using NMR. The free structure is similar to the bound state, and there are no significant main-chain differences in the 27-44 region involved in barstar binding to barnase. The C82A structure shows significant differences from the average NMR structure, both overall and in the binding region. In contrast to the crystal structure, the NMR structure shows an unusually high packing value based on the occluded surface algorithm, indicating errors in the packing of the structure. We show that the NMR structures of homologous proteins generally show large differences in packing value, while the crystal structures of such proteins have very similar packing values, suggesting that protein packing density is not well determined by NMR.

Multiple kinetic intermediates accumulate during the unfolding of horse cytochrome c in the oxidized state.

The unfolding kinetics of horse cytochrome c in the oxidized state has been studied at 10, 22, and 34 degreesC as a function of guanidine hydrochloride (GdnHCl) concentration. Rapid (millisecond) measurements of far-UV circular dichroism (CD) as well as fluorescence quenching due to tryptophan to heme excitation energy transfer have been used to monitor the unfolding process. At 10 degreesC, the decrease in far-UV CD signal that accompanies unfolding occurs in two phases. The unobservable burst phase is complete within 4 ms, while the slower phase occurs over tens to hundreds of milliseconds. The burst phase unfolding amplitude increases cooperatively with an increase in GdnHCl concentration, exhibiting a transition midpoint of 3.2 M at 10 degreesC. In contrast, no burst phase change in fluorescence occurs during unfolding at 10 degreesC. At 22 and 34 degreesC, both the fluorescence-monitored unfolding kinetics and the far-UV CD-monitored unfolding kinetics are biphasic. At both temperatures, the two probes yield burst phase unfolding transitions that are noncoincident with respect to the transition midpoints as well as the dependency of the burst phase amplitudes on GdnHCl concentration. The results suggest that at least two kinetic unfolding intermediates accumulate during unfolding. One burst phase intermediate, IU1, has lost virtually all the native-state secondary structure, while the other burst phase intermediate, IU2, has lost both secondary structure and native-like compactness. The presence of kinetic unfolding intermediates is also indicated by the nonlinear dependence of the logarithm of the apparent unfolding rate constant on GdnHCl concentration, which is particularly pronounced at 10 and 22 degreesC. Analysis of the burst phase unfolding transitions obtained using the two probes shows that the stabilities of IU1 and IU2 decrease steadily with an increase in temperature from 10 to 34 degreesC, suggesting that the structures present in them are stabilized principally by hydrogen bonding interactions.

Reversible unfolding of sheep liver tetrameric serine hydroxymethyltransferase.

Equilibrium unfolding studies of the tetrameric serine hydroxymethyltransferase from sheep liver (SHMT, E.C.2.1.2.1) revealed that the holoenzyme, apoenzyme and the sodium borohydride-reduced holoenzyme had random coil structures in 8 M urea. In the presence of a non-ionic detergent, Brij-35, and polyethylene glycol, the 8 M urea unfolded protein could be completely (> 95%) refolded by a 20-fold dilution. The refolded enzyme was completely active and kinetically similar to the native enzyme. The midpoint of inactivation of the enzyme occurred at a urea concentration that was much below the urea concentration required to bring about a substantial loss of secondary structure. This observation suggested the occurrence of a 'predenaturation transition' in the unfolding pathway. The equilibrium urea-induced denaturation curve of holoSHMT showed two transitions. The midpoint of the first transition was 1.2 M, which was comparable to that required for 50% decrease in enzyme activity. Further, 50% release of the pyridoxal-5'-phosphate (PLP) from the active site, as monitored by decrease in absorbance at 425 nm, also occurred at about 1.2 M urea. Size exclusion chromatography showed that the tetrameric SHMT unfolds via the intermediate formation of dimers. This dissociation occurred at a much lower urea concentration (0.15 M) in the unfolding of the apoenzyme, and at a higher urea concentration (1.2 M) in the unfolding of holoenzyme, thereby demonstrating the involvement of PLP in stabilizing the quaternary structure of the enzyme. Size exclusion chromatography of the refolding intermediates demonstrated that the cofactor shifts the equilibrium towards the formation of the active tetramer. The reduced holoenzyme could also be refolded to its native structure, as observed by fluorescence and CD measurements, indicating that the presence of covalently linked PLP does not affect refolding. The results demonstrate clearly that the dimer is an intermediate in the urea-induced equilibrium unfolding/refolding of sheep liver SHMT; and PLP, in addition to its role in catalysis, is required for the stabilization of the tetrameric structure of the enzyme.

Two structural subdomains of barstar detected by rapid mixing NMR measurement of amide hydrogen exchange.

Equilibrium amide hydrogen exchange studies of barstar have been carried out at pH 6.7, 32 degrees C using one- and two-dimensional nuclear magnetic resonance. An unusually large fraction of the backbone amide hydrogens of barstar exchange too fast to be measured, and the exchange rates of only fifteen slow-exchanging amide sites including indole amides of two tryptophans could be measured in the presence of 0 to 1.8 M guanidine hydrochloride (GdnHCl). Measurement of exchange occurring in tens of seconds in the unfolding transition region was possible by the use of a fast stopped-flow mixing method. The observed exchange rates have been simulated in the EX2 limit according to a two-process model that incorporates two exchange-competent states: a transiently unfolded state (U*) in which many amide hydrogens are completely accessible to solvent-exchange, and a near-native locally unfolded state (N*), in which only one or a few amide hydrogens are completely accessible to solvent-exchange. The two-process model appears to account for the observed exchange behavior over the entire range of GdnHCl concentrations studied. For several measurable slow-exchanging amide hydrogens, the free energies of production of exchange-competent states from the exchange-incompetent native state are significantly higher than the free-energy of production of the equilibrium unfolded state from the native state, when the latter is determined from circular dichroism- or fluorescence-monitored equilibrium unfolding curves. The result implies that U*, which forms transiently in the strongly native-like conditions used for the hydrogen exchange studies, is higher in energy than the equilibrium-unfolded state. The higher energy of this transiently unfolded exchange-competent state can be attributed to either proline isomerization or to the presence of residual structure. On the basis of the free energies of production of exchange-competent states, the measured amide sites of barstar appear to define two structural subdomains--a three-helix unit and a two-beta-strand unit in the core of the protein.

Multiple intermediates and transition states during protein unfolding.

Rapid kinetic studies of the unfolding of the small protein barstar by urea have been used to demonstrate the presence of at least two unfolding intermediates on two competing unfolding pathways. One intermediate has native-like secondary structure but has a partially solvated hydrophobic core, while the other is devoid of considerable secondary structure but has an intact hydrophobic core. It is shown that the transition states on the two pathways are very dissimilar structurally, but very similar energetically.

DSC studies of the conformational stability of barstar wild-type.

The temperature induced unfolding of barstar wild-type of bacillus amyloliquefaciens (90 residues) has been characterized by differential scanning microcalorimetry. The process has been found to be reversible in the pH range from 6.4 to 8.3 in the absence of oxygen. It has been clearly shown by a ratio of delta HvH/delta Hcal near 1 that denaturation follows a two-state mechanism. For comparison, the C82A mutant was also studied. This mutant exhibits similar reversibility, but has a slightly lower transition temperature. The transition enthalpy of barstar wt (303 kJ mol-1) exceeds that of the C82A mutant (276 kJ mol-1) by approximately 10%. The heat capacity changes show a similar difference, delta Cp being 5.3 +/- 1 kJ mol-1 K-1 for the wild-type and 3.6 +/- 1 kJ mol-1 K-1 for the C82A mutant. The extrapolated stability parameters at 25 degrees C are delta G0 = 23.5 +/- 2 kJ mol-1 for barstar wt and delta G0 = 25.5 +/- 2 kJ mol-1 for the C82A mutant.

Thermodynamics of the complex protein unfolding reaction of barstar.

The complex unfolding reaction of barstar has been characterized by studying the apparent rate of unfolding, monitored by intrinsic Trp fluorescence, as a function of temperature and guanidine hydrochloride (GdnHCl) concentration. The kinetics of unfolding and folding of wild-type (wt) barstar at 5 degrees C were first studied in detail. It is shown that when unfolding is carried out using concentrations of GdnHCl in the posttransition zone of unfolding, the change in fluorescence that accompanies unfolding occurs in two phases: 30% of the change occurs in a burst phase that is complete within 4 ms, and 70% of the change occurs in a fast phase that is complete within 2 s. In contrast, when the protein is unfolded at 25 degrees C, no burst-phase change in fluorescence is observed. To confirm that a burst-phase change in fluorescence indeed accompanies unfolding at low temperature, unfolding studies were also carried out on a marginally destabilized mutant form of barstar for which the burst-phase change in fluorescence is shown to be as high as 70%. These results confirm a previous report [Nath et al., (1996), Nat. Struct. Biol. 3, 920-923], in which the detection of a burst-phase change in circular dichroism at 222 nm during unfolding at 25 degrees C led to the inclusion of a rapidly formed kinetic intermediate, IU, on the unfolding pathway. To characterize thermodynamically the unfolding pathway, apparent unfolding rates were then measured at six different concentrations of GdnHCl in the range 2.6 to 5.0 M, at five different temperatures from 5 to 46 degrees C. The subsequent analysis was done on the basis of the observation that a preequilibrium between the fully folded state (F) and IU gets established rapidly before further unfolding to the completely unfolded state (U). The results indicate that IU has a specific heat capacity similar to that of F and therefore suggest that IU is as compact as F, with practically no exposure of the hydrophobic core. On the other hand, the transition state of unfolding has a 45% greater heat capacity than F, indicating that significant hydration of the hydrophobic core occurs only after the rate-limiting step of unfolding.