Equilibrium unfolding thermodynamics of beta2-microglobulin analyzed through native-state H/D exchange

The exchange rates for the amide hydrogens of β₂-microglobulin, the protein responsible for dialysis-related amyloidosis, were measured under native conditions at different temperatures ranging from 301 to 315 K. The pattern of protection factors within different regions of the protein correlates well with the hydrogen-bonding pattern of the deposited structures. Analysis of the exchange rates indicates the presence of mixed EX1- and EX2-limit mechanisms. The measured parameters are consistent with a two-process model in which two competing pathways, i.e., global unfolding in the core region and partial openings of the native state, determine the observed exchange rates. These findings are analyzed with respect to the amyloidogenic properties of the protein.     

Characterization of the unfolding of ribonuclease a by a pulsed hydrogen exchange study: evidence for competing pathways for unfolding

The unfolding of ribonuclease A was studied in 5.2 M guanidine hydrochloride at pH 8 and 10 degrees C using multiple optical probes, native-state hydrogen exchange (HX), and pulse labeling by hydrogen exchange. First, native-state HX studies were used to demonstrate that the protein exists in two slowly interconverting forms under equilibrium native conditions: a predominant exchange-incompetent N form and an alternative ensemble of conformations, N(I), in which some amide hydrogens are fully exposed to exchange. Pulsed HX studies indicated that, during unfolding, the rates of exposure to exchange with solvent protons were similar for all backbone NH probe protons. It is shown that two parallel routes of unfolding are available to the predominant N conformation as soon as it encounters strong unfolding conditions. A fraction of molecules appears to rapidly form N(I) on one route. On the other route an exchange-incompetent intermediate state ensemble, I(U)(2), is formed. The kinetics of unfolding measured by far-UV circular dichroism (CD) were faster than those measured by near-UV CD and intrinsic tyrosine fluorescence of the protein. The logarithms of the rate constants of the unfolding reaction measured by all three optical probes also showed a nonlinear dependence on GdnHCl concentration. All of the data suggest that N(I) and I(U)(2) are nativelike in their secondary and tertiary structures. While N(I) unfolds directly to the fully exchange-competent unfolded state (U), I(U)(2) forms another intermediate I(U)(3) which then unfolds to U. I(U)(3) is devoid of all native alpha-helical secondary structure and has only 30% of the tertiary interactions still intact. Since the rates of global unfolding measured by near-UV CD and fluorescence agree well with the rates of exposure determined for all of the backbone NH probe protons, it appears that the rate-limiting step for the unfolding of RNase A is the dissolution of the entire native tertiary structure and penetration of water into the hydrophobic core.     

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.     

Unfolding rates of barstar determined in native and low denaturant conditions indicate the presence of intermediates

The folding and unfolding rates of the small protein, barstar, have been monitored using stopped-flow measurements of intrinsic tryptophan fluorescence at 25 degrees C, pH 8.5, and have been compared over a wide range of urea and guanidine hydrochloride (GdnHCl) concentrations. When the logarithms of the rates of folding from urea and from GdnHCl unfolded forms are extrapolated linearly with denaturant concentration, the same rate is obtained for folding in zero denaturant. Similar linear extrapolations of rates of unfolding in urea and GdnHCl yield, however, different unfolding rates in zero denaturant, indicating that such linear extrapolations are not valid. It has been difficult, for any protein, to determine unfolding rates under nativelike conditions in direct kinetic experiments. Using a novel strategy of coupling the reactivity of a buried cysteine residue with 5,5'-dithiobis(2-nitrobenzoic acid) (DTNB) to the unfolding reaction of barstar, the global unfolding and refolding rates have now been determined in low denaturant concentrations. The logarithms of unfolding rates obtained at low urea and GdnHCl concentrations show a markedly nonlinear dependence on denaturant concentration and converge to the same unfolding rate in the absence of denaturant. It is shown that the native protein can sample the fully unfolded conformation even in the absence of denaturant. The observed nonlinear dependences of the logarithms of the refolding and unfolding rates observed for both denaturants are shown to be due to the presence of (un)folding intermediates and not due to movements in the position of the transition state with a change in denaturant concentration.     

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. Proteins 32:241–247, 1998. © 1998 Wiley-Liss, Inc.     

Folding subdomains of thioredoxin characterized by native-state hydrogen exchange

Native-state hydrogen exchange (HX) studies, used in conjunction with NMR spectroscopy, have been carried out on Escherichia coli thioredoxin (Trx) for characterizing two folding subdomains of the protein. The backbone amide protons of only the slowest-exchanging 24 amino acid residues, of a total of 108 amino acid residues, could be followed at pH 7. The free energy of the opening event that results in an amide hydrogen exchanging with solvent (DeltaG(op)) was determined at each of the 24 amide hydrogen sites. The values of DeltaG(op) for the amide hydrogens belonging to residues in the helices alpha(1), alpha(2), and alpha(4) are consistent with them exchanging with the solvent only when the fully unfolded state is sampled transiently under native conditions. The denaturant-dependences of the values of DeltaG(op) provide very little evidence that the protein samples partially unfolded forms, lower in energy than the unfolded state. The amide hydrogens belonging to the residues in the beta strands, which form the core of the protein, appear to have higher values of DeltaG(op) than amide hydrogens belonging to residues in the helices, suggesting that they might be more stable to exchange. This apparently higher stability to HX of the beta strands might be either because they exchange out their amide hydrogens in a high energy intermediate preceding the globally unfolded state, or, more likely, because they form residual structure in the globally unfolded state. In either case, the central beta strands-beta(3,) beta(2), and beta(4)-would appear to form a cooperatively folding subunit of the protein. The native-state HX methodology has made it possible to characterize the free energy landscape that Trx can sample under equilibrium native conditions.     

Role of native-state structure in rubredoxin native-state hydrogen exchange

Amide exchange rates were measured for Pyrococcus furiosus (Pf) rubredoxin substituted with either Zn(II), Ga(III), or Ge(IV). Base-catalyzed exchange rate constants increase up to 3000-fold per unit charge for the highly protected amides surrounding the active site metal, yielding apparent residue-specific conformational energy decreases of more than 8 kcal/mol in a comparison of the Zn(II)- and Ge(IV)-substituted proteins. However, the exchange kinetics for many of the other amides of the protein are insensitive to these metal substitutions. These differential rates are inversely correlated with the distance between the amide nitrogen and the metal in the X-ray structure, out to a distance of at least 12 A, consistent with an electrostatic potential-dependent shifting of the amide nitrogen pK. This strongly correlated distance dependence is consistent with a nativelike structure for the exchange-competent conformations. The electric field potential within the interior of the rubredoxin structure gives rise to a change of as much as a million-fold in the rate for the exchange-competent state of the individual amide hydrogens. Nevertheless, the strength of these electrostatic interactions in Pf rubredoxin appears to be comparable to those previously reported within other proteins. As a result, contrary to the conventional analysis of hydrogen exchange data, for exchange processes that occur via nonglobal transitions, the residual conformational structure will often modulate the observed rates. Although this necessarily complicates the estimation of the conformational equilibria of these exchange-competent states, this dependence on residual structure can provide insight into the conformation of these transient states.     

A two-process model describes the hydrogen exchange behavior of cytochrome c in the molten globule state with various extents of acetylation

Acetylation of Lys residues of horse cytochrome c steadily stabilizes the molten globule state in 18 mM HCl as more Lys residues are acetylated [Goto and Nishikiori (1991) J. Mol. Biol. 222, 679-686]. The dynamic features of the molten globule state were characterized by hydrogen/deuterium exchange of amide protons, monitored by mass spectrometry as each deuteration increased the protein mass by 1 Da. Electrospray mass spectrometry enabled us to monitor simultaneously the exchange kinetics of more than seven species with a different number of acetyl groups. One to four Lys residue-acetylated cytochrome c showed almost no protection of the amide protons from rapid exchange. The transition from the unprotected to the protected state occurred between five and eight Lys residue-acetylated species. For species with more than nine acetylated Lys residues, the exchange kinetics were independent of the extent of acetylation, and 26 amide protons were protected at 60 min of exchange, indicating the formation of a rigid hydrophobic core with hydrogen-bonded secondary structures. The apparent transition to the protected state required a higher degree of acetylation than the conformational transition measured by circular dichroism, which had a midpoint at about four acetylated residues. This difference in the transitions suggested a two-process model in which the exchange occurs either from the protected folded state or from the unprotected unfolded state through global unfolding. On the basis of a two-process model and with the reported values of the exchange and stability parameters, we simulated the exchange kinetics of a series of acetylated cytochrome c species. The simulated kinetics reproduced the observed kinetics well, indicating validity of this model for hydrogen exchange of the molten globule state.     

Characterization of the structure and dynamics of a near-native equilibrium intermediate in the unfolding pathway of an all beta-barrel protein

The structure and dynamics of equilibrium intermediate in the unfolding pathway of the human acidic fibroblast growth factor (hFGF-1) are investigated using a variety of biophysical techniques including multidimensional NMR spectroscopy. Guanidinium hydrochloride (GdnHCl)-induced unfolding of hFGF-1 proceeds with the accumulation of a stable intermediate state. The transition from the intermediate state to the unfolded state(s) is cooperative without the accumulation of additional intermediate(s). The intermediate state induced maximally in 0.96 m GdnHCl is found to be obligatory in the folding/unfolding pathway of hFGF-1. Most of the native tertiary structure interactions are preserved in the intermediate state. (1)H-(15)N chemical shift perturbation data suggest that the residues in the C-terminal segment including those located in the beta-strands IX, X, and XI undergo the most discernible structural change(s) in the intermediate state in 0.96 m GdnHCl. hFGF-1 in the intermediate state (0.96 m GdnHCl) does not bind to its ligand, sucrose octasulfate. Limited proteolytic digestion experiments and hydrogen-deuterium exchange monitored by (15)N heteronuclear single quantum coherence (HSQC) spectra show that the conformational flexibility of the protein in the intermediate state is significantly higher than in the native conformation. (15)N spin relaxation experiments show that many residues located in beta-strands IX, X, and XI exhibit conformational motions in the micro- to millisecond time scale. Analysis of (15)N relaxation data in conjunction with the amide proton exchange kinetics suggests that residues in the beta-strands II, VIII, and XII possibly constitute the stability core of the protein in the near-native intermediate state.     

Hydrogen exchange in a large 29 kD protein and characterization of molten globule aggregation by NMR

The nature of denatured ensembles of the enzyme human carbonic anhydrase (HCA) has been extensively studied by various methods in the past. The protein constitutes an interesting model for folding studies that does not unfold by a simple two-state transition, instead a molten globule intermediate is highly populated at 1.5 M GuHCl. In this work, NMR and H/D exchange studies have been conducted on one of the isozymes, HCA I. The H/D exchange studies, which were enabled by the previously obtained resonance assignment of HCA I, have been used to identify unfolded forms that are accessible from the native state. In addition, the GuHCl-induced unfolded states of HCA I have also been characterized by NMR at GuHCl concentrations in the 0-5 M range. The most important findings in this work are as follows: (1) Amide protons located in the center of the beta-sheet require global unfolding events for efficient H/D exchange. (2) The molten globule and the native state give similar protection against H/D exchange for all of the observable amide protons (i.e., water seems not to efficiently penetrate the interior of the molten globule). (3) At high protein concentrations, the molten globule can form large aggregates, which are not detectable by solution-state NMR methods. (4) The unfolded state (U), present at GuHCl concentrations above 2 M, is composed of an ensemble of conformations having residual structures with different stabilities.     

Hydrogen exchange study of canine milk lysozyme: stabilization mechanism of the molten globule

The native state (1)H, (15)N resonance assignment of 123 of the 128 nonproline residues of canine milk lysozyme has enabled measurements of the amide hydrogen exchange of over 70 amide hydrogens in the molten globule state. To elucidate the mechanism of protein folding, the molten globule state has been studied as a model of the folding intermediate state. Lysozyme and alpha-lactalbumin are homologous to each other, but their equilibrium unfolding mechanisms differ. Generally, the folding mechanism of lysozyme obeys a two-state model, whereas that of alpha-lactalbumin follows a three-state model. Exceptions to this rule are equine and canine milk lysozymes, which exhibit a partially unfolded state during the equilibrium unfolding; this state resembles the molten globule state of alpha-lactalbumin but with extreme stability. Study of the molten globules of alpha-lactalbumin and equine milk lysozyme showed that the stabilities of their alpha-helices are similar, despite the differences in the thermodynamic stability of their molten globule states. On the other hand, our hydrogen exchange study of the molten globule of canine milk lysozyme showed that the alpha-helices are more stabilized than in alpha-lactalbumin or equine milk lysozyme and that this enhanced stability is caused by the strengthened cooperative interaction between secondary structure elements. Thus, our results underscore the importance of the cooperative interaction in the stability of the molten globule state.     

Quantitative analysis of the kinetics of denaturation and renaturation of barstar in the folding transition zone.

The fluorescence-monitored kinetics of folding and unfolding of barstar by guanidine hydrochloride (GdnHCl) in the folding transition zone, at pH 7, 25 degrees C, have been quantitatively analyzed using a 3-state mechanism: U(S)<-->UF<-->N. U(S) and UF are slow-refolding and fast-refolding unfolded forms of barstar, and N is the native protein. U(S) and UF probably differ in possessing trans and cis conformations, respectively, of the Tyr 47-Pro 48 bond. The 3-state model could be used because the kinetics of folding and unfolding of barstar show 2 phases, a fast phase and a slow phase, and because the relative amplitudes of the 2 phases depend only on the final refolding conditions and not on the initial conditions. Analysis of the observed kinetics according to the 3-state model yields the values of the 4 microscopic rate constants that describe the transitions between the 3 states at different concentrations of GdnHCl. The value of the equilibrium unfolded ratio U(S):UF (K21) and the values of the rate constants of the U(S)-->UF and UF-->U(S) reactions, k12 and k21, respectively, are shown to be independent of the concentration of GdnHCl. K21 has a value of 2.1 +/- 0.1, and k12 and k21 have values of 5.3 x 10(-3) s-1 and 11.2 x 10(-3) s-1, respectively. Double-jump experiments that monitor reactions that are silent to fluorescence monitoring were used to confirm the values of K21, k12, and k21 obtained from the 3-state analysis and thereby the validity of the 3-state model.(ABSTRACT TRUNCATED AT 250 WORDS)     

Insights into conformation and dynamics of protein GB1 during folding and unfolding by NMR

Understanding protein stability requires characterization of structural determinants of the folded and unfolded states. Many proteins are capable of populating partially folded states under specific solution conditions. Occasionally, coexistence of the folded and an unfolded state under non- or mildly denaturing conditions can be observed by NMR, allowing us to structurally probe these states under identical conditions. Here we report on a destabilized mutant of the B1 domain of protein G (GB1) whose equilibrium unfolding was systematically investigated. Backbone amide residual dipolar couplings (RDCs), the tryptophan Nepsilon-H resonance and the amide nitrogen transverse relaxation rates (R2s) for varying pH values and different temperatures were measured. The backbone amide RDCs indicate that prior to complete unfolding, two melting hot spots are formed at the turn around T11, L12 and K13 and the N terminus of the helix at A24 and T25. The RDCs for the low pH, thermally unfolded state of GB1 are very small and do not indicate the presence of any native-like structure. Amide nitrogen transverse relaxation rates for GB1 in the folded state at different temperatures exhibit large contributions from exchange processes and the associated dynamics display considerable heterogeneity. Our data provide clear evidence for intermediate conformations and multi-state equilibrium un/folding for this GB1 variant.     

pH dependent thermodynamic and amide exchange studies of the C-terminal domain of the ribosomal protein L9: implications for unfolded state structure

It is now recognized that unfolded states of globular proteins are not random coils but instead can contain significant amounts of residual structure. Here, we combine amide H/D exchange studies and thermodynamic measurements to probe pH dependent structure in the unfolded state of the small, mixed alpha-beta protein CTL9. The m value measured by urea denaturation is strongly dependent upon pD, increasing by 40% from pD 7.5 to 4.85. Likewise, the change in heat capacity upon unfolding, deltaCp(o), increases significantly from pD 7.5 to 5.5. These studies argue that the unfolded state contains interactions, presumably hydrophobic in nature, that lead to a more compact state at high pH. The expansion at lower pH correlates with the estimated unfolded state pKa values of the three histidines in CTL9 with additional contributions from acid side chains at the lower pH. Amide H/D exchange studies were conducted at pD 5.0, 6.0, and 7.0. At pD 5.0, the exchange rates could be measured for 44 residues, 29 of which exchanged by global unfolding. No evidence was found for any super protected sites, that is, sites that exchange at rates slower than those expected for global exchange. The estimated precision for the experiments limits detection to residues that are protected 2.3-fold above the intrinsic exchange rate. Thirty-seven residues could be followed at pD 6 and 27 residues at pD 7. Again no evidence for a significant super protected structure was observed. The properties of CTL9(11) are compared to other structured denatured states.     

Test for cooperativity in the early kinetic intermediate in lysozyme folding

During the folding of many proteins, collapsed globular states are formed prior to the native structure. The role of these states for the folding process has been widely discussed. Comparison with properties of synthetic homo and heteropolymers had suggested that the initial collapse represented a shift of the ensemble of unfolded conformations to more compact states without major energy barriers. We investigated the folding/unfolding transition of a collapsed state, which transiently populates early in lysozyme folding. This state forms within the dead-time of stopped-flow mixing and it has been shown to be significantly more compact and globular than the denaturant-induced unfolded state. We used the GdmCl-dependence of the dead-time signal change to characterize the unfolding transition of the burst phase intermediate. Fluorescence and far-UV CD give identical unfolding curves, arguing for a cooperative two-state folding/unfolding transition between unfolded and collapsed lysozyme. These results show that collapse leads to a distinct state in the folding process, which is separated from the ensemble of unfolded molecules by a significant energy barrier. NMR, fluorescence and small angle X-ray scattering data further show that some local interactions in unfolded lysozyme exist at denaturant concentrations above the coil-collapse transition. These interactions might play a crucial role in the kinetic partitioning between fast and slow folding pathways.     

Revealing a Concealed Intermediate that Forms after the Rate-limiting Step of Refolding of the SH3 Domain of PI3 Kinase

Kinetic and equilibrium studies of the folding and unfolding of the SH3 domain of the PI3 kinase, have been used to identify a folding intermediate that forms after the rate-limiting step on the folding pathway. Folding and unfolding, in urea as well as in guanidine hydrochloride (GdnHCl), were studied by monitoring changes in the intrinsic fluorescence or in the far-UV circular dichroism (CD) of the protein. The two probes yield non-coincident equilibrium transitions for unfolding in urea, indicating that an intermediate, I, exists in equilibrium with native (N) and unfolded (U) protein, during unfolding. Hence, the equilibrium unfolding data were analyzed according to a three-state N ↔ I ↔ U mechanism. An intermediate is observed also in kinetic unfolding studies, and its presence leads to the unfolding reaction in urea as well as in GdnHCl, occurring in two steps. The fast step is complete within the initial 11 ms of unfolding and manifests itself in a burst phase change in fluorescence. At high concentrations of GdnHCl, the entire change in fluorescence during unfolding occurs during the 11 ms burst phase. CD measurements indicate, however, that I retains N-like secondary structure. An analysis of the kinetic and thermodynamic data, according to a minimal three-state N ↔ I ↔ U mechanism, positions I after the rate-limiting transition state, TS1, of folding, on the reaction coordinate of folding in GdnHCl. Hence, I is not revealed when folding is commenced from U, regardless of the nature of the probe used to follow the folding reaction. Interrupted unfolding experiments, in which the protein is unfolded transiently in GdnHCl for various lengths of time before being refolded, showed that I refolds to N much faster than does U, confirms the analysis of the direct folding and unfolding experiments, that I is formed after the rate-limiting step of refolding in GdnHCl.     

pH dependence of the stability of barstar to chemical and thermal denaturation.

Equilibrium unfolding of barstar with guanidine hydrochloride (GdnHCl) and urea as denaturants as well as thermal unfolding have been carried out as a function of pH using fluorescence, far-UV and near-UV CD, and absorbance as probes. Both GdnHCl-induced and urea-induced denaturation studies at pH 7 show that barstar unfolds through a two-state F<->U mechanism and yields identical values for delta GU, the free energy difference between the fully folded (F) and unfolded (U) forms, of 5.0 +/- 0.5 kcal.mol-1 at 25 degrees C. Thermal denaturation of barstar also follows a two-state F<->U unfolding transition at pH 7, and the value of delta GU at 25 degrees C is similar to that obtained from chemical denaturation. The pH dependence of denaturation by GdnHCl is complex. The Cm value (midpoint of the unfolding transition) has been used as an index for stability in the pH range 2-10, because barstar does not unfold through a two-state transition on denaturation by GdnHCl at all pH values studied. Stability is maximum at pH 2-3, where barstar exists in a molten globule-like form that forms a large soluble oligomer. The stability decreases with an increase in pH to 5, the isoelectric pH of the protein. Above pH 5, the stability increases as the pH is raised to 7. Above pH 8, it again decreases as the pH is raised to 10. The decrease in stability from pH 7 to 5 in wild-type (wt) barstar, which is shown to be characterized by an apparent pKa of 6.2 +/- 0.2, is not observed in H17Q, a His 17-->Gln 17 mutant form of barstar. This decrease in stability has therefore been correlated with the protonation of His 17 in barstar. The decrease in stability beyond pH 8 in wt barstar, which is characterized by an apparent pKa of 9.2 +/- 0.2, is not detected in BSCCAA, the Cys 40 Cys 82-->Ala 40 Ala 82 double mutant form of barstar. Thus, this decrease in stability has been correlated with the deprotonation of at least one of the two cysteines present in wt barstar. The increase in stability from pH 5 to 3 is characterized by an apparent pKa of 4.6 +/- 0.2 for wt barstar and BSCCAA, which is similar to the apparent pKa that characterizes the structural transition leading to the formation of the A form. The use of Cm as an index of stability has been supported by thermal denaturation studies.(ABSTRACT TRUNCATED AT 400 WORDS)     

Characterization of equilibrium intermediates in denaturant-induced unfolding of ferrous and ferric cytochromes c using magnetic circular dichroism, circular dichroism, and optical absorption spectroscopies

Protein unfolding during guanidine HCl denaturant titration of the reduced and oxidized forms of cytochrome c is monitored with magnetic circular dichroism (MCD), natural CD, and absorption of the heme bands and far-UV CD of the amide bands. Direct MCD spectral evidence is presented for bis-histidinyl heme ligation in the unfolded states of both the reduced and oxidized protein. For both redox states, the unfolding midpoints measured with MCD, which is an indicator of tertiary structure, are significantly lower than those measured with far-UV CD, an indicator of secondary structure. The disparate titration curves are interpreted in terms of a compound mechanism for denaturant-induced folding and unfolding involving a molten globulelike intermediate state (MG) with near-native secondary structure and nonnative tertiary structure and heme ligation. A comparison of the dependence of the free energy of formation of the MG intermediate on the redox state with the known contributions from heme ligation and solvation suggests that the heme is significantly more accessible to solvent in the MG intermediate than it is in the native state.     

Beyond the EX1 limit: probing the structure of high-energy states in protein unfolding

Hydrogen exchange kinetics in native solvent conditions have been used to explore the conformational fluctuations of an immunoglobulin domain (CD2.domain1). The global folding/unfolding kinetics of the protein are unaltered between pH 4.5 and pH 9.5, allowing us to use the pH-dependence of amide hydrogen/deuterium exchange to characterise conformational states with energies up to 7.2kcal/mol higher than the folded ground state. The study was intended to search for discreet unfolding intermediates in this region of the energy spectrum, their presence being revealed by the concerted exchange behaviour of subsets of amide groups that become accessible at a given free energy, i.e. the spectrum would contain discreet groupings. Protection factors for 58 amide groups were measured across the pH range and the hydrogen-exchange energy profile is described.More interestingly, exchange behaviour could be grouped into three categories; the first two unremarkable, the third unexpected. (1) In 33 cases, amide exchange was dominated by rapid fluctuation, i.e. the free energy difference between the ground state and the rapidly accessed open state is sufficiently low that the contribution from crossing the unfolding barrier is negligible. (2) In 18 cases exchange is dominated by the global folding transition barrier across the whole pH range measured. The relationship between hydroxyl ion concentration and observed exchange rate is hyperbolic, with the limiting rate being that for global unfolding; the so-called EX1 limit. For these, the free energy difference between the folded ground state and any rapidly-accessed open state is too great for the proton to be exchanged through such fluctuations, even at the highest pH employed in this study. (3) For the third group, comprising five cases, we observe a behaviour that has not been described. In this group, as in category 2, the rate of exchange reaches a plateau; the EX1 limit. However, as the intrinsic exchange rate (k(int)) is increased, this limit is breached and the rate begins to rise again. This unintuitive behaviour does not result from pH instability, rather it is a consequence of amide groups experiencing two processes; rapid fluctuation of structure and crossing the global barrier for unfolding. The boundary at which the EX1 limit is overcome is determined by the equilibrium distribution of the fluctuating open and closed states (K(O/C)) and the rate constant for unfolding (k(u)). This critical boundary is reached when k(int)K(O/C)=k(u). Given that, in a simple transition state formalism: k(u)=K(#)k' (where K(#) describes the equilibrium distribution between the transition and ground state and k' describes the rate of a barrierless rearrangement), it follows that if the pH is raised to a level where k(int)=k', then the entire free energy spectrum from ground state to transition state could be sampled.