Reversible unfolding and refolding behavior of a monomeric aldolase from Staphylococcus aureus.

Thermal and GdmCl-induced unfolding transitions of aldolase from Staphylococcus aureus are reversible under a variety of solvent conditions. Analysis of the transitions reveals that no partially folded intermediates can be detected under equilibrium conditions. The stability of the enzyme is very low with a delta G0 value of -9 +/- 2 kJ/mol at 20 degrees C. The kinetics of unfolding and refolding of aldolase are complex and comprise at least one fast and two slow reactions. This complexity arises from prolyl isomerization reactions in the unfolded chain, which are kinetically coupled to the actual folding reaction. Comparison with model calculations shows that at least two prolyl peptide bonds give rise to the observed slow folding reactions of aldolase and that all of the involved bonds are presumably in the trans conformation in the native state. The rate constant of the actual folding reaction is fast with a relaxation time of about 15 s at the midpoint of the folding transition at 15 degrees C. The data presented on the folding and stability of aldolase are comparable to the properties of much smaller proteins. This might be connected with the simple and highly repetitive tertiary structure pattern of the enzyme, which belongs to the group of alpha/beta barrel proteins.     

Equilibrium and kinetic folding pathways of a TIM barrel with a funneled energy landscape

The role of native contact topology in the folding of a TIM barrel model based on the alpha-subunit of tryptophan synthase (alphaTS) from Salmonella typhimurium (Protein Data Bank structure 1BKS) was studied using both equilibrium and kinetic simulations. Equilibrium simulations of alphaTS reveal the population of two intermediate ensembles, I1 and I2, during unfolding/refolding at the folding temperature, Tf = 335 K. Equilibrium intermediate I1 demonstrates discrete structure in regions alpha0-beta6 whereas intermediate I2 is a loose ensemble of states with N-terminal structure varying from at least beta1-beta3 (denoted I2A) to alpha0-beta4 at most (denoted I2B). The structures of I1 and I2 match well with the two intermediate states detected in equilibrium folding experiments of Escherichia coli alphaTS. Kinetic folding simulations of alphaTS reveal the sequential population of four intermediate ensembles, I120Q, I200Q, I300Q, and I360Q, during refolding. Kinetic intermediates I120Q, I200Q, and I300Q are highly similar to equilibrium alphaTS intermediates I2A, I2B, and I1, respectively, consistent with kinetic experiments on alphaTS from E. coli. A small population (approximately 10%) of kinetic trajectories are trapped in the I120Q intermediate ensemble and require a slow and complete unfolding step to properly refold. Both the on-pathway and off-pathway I120Q intermediates show structure in beta1-beta3, which is also strikingly consistent with kinetic folding experiments of alphaTS. In the off-pathway intermediate I(120Q), helix alpha2 is wrapped in a nonnative chiral arrangement around strand beta3, sterically preventing the subsequent folding step between beta3 and beta4. These results demonstrate the success of combining kinetic and equilibrium simulations of minimalist protein models to explore TIM barrel folding and the folding of other large proteins.     

Folding of ribonuclease T1. 2. Kinetic models for the folding and unfolding reactions

The slow refolding of ribonuclease T1 was investigated by different probes. Structural intermediates with secondary structure are formed early during refolding, as indicated by the rapid regain of a native-like circular dichroism spectrum in the amide region. This extensive structure formation is much faster than the slow steps of refolding, which are limited in rate by the reisomerization of incorrect proline isomers. The transient folding intermediates were also detected by unfolding assays, which make use of the reduced stability of folding intermediates relative to that of the native protein. The results of this and the preceding paper [Kiefhaber et al. (1990) Biochemistry (preceding paper in this issue)] were used to propose kinetic models for the unfolding and refolding of ribonuclease T1. The unfolding mechanism is based on the assumption that, after the structural unfolding step, the slow isomerizations of two X-Pro peptide bonds occur independently of each other in the denatured protein. At equilibrium a small amount of fast-folding species coexists with three slow-folding species: two with one incorrect proline isomer each and another, dominant species with both these prolines in the incorrect isomeric state. In the mechanism for refolding we assume that all slow-folding molecules can rapidly regain most of the secondary and part of the tertiary structure early in folding. Reisomerizations of incorrect proline peptide bonds constitute the slow, rate-limiting steps of refolding. A peculiar feature of the kinetic model for refolding is that the major unfolded species with two incorrect proline isomers can enter two alternative folding pathways, depending on which of the two reisomerizes first. The relative rates of reisomerization of the respective proline peptide bonds at the stage of the rapidly formed intermediate determine the choice of pathway. It is changed in the presence of prolyl isomerase, because this enzyme catalyzes these two isomerizations with different efficiency and consequently leads to a shift from the very slow to the intermediate refolding pathway.     

Preformed secondary structure drives the association reaction of GCN4-p1, a model coiled-coil system

The structure of the transition state for the rate-limiting step in the folding and association of the homodimeric coiled-coil peptide GCN4-p1, was probed by mutational analysis. A series of quadruple amino acid replacements that spanned the helix propensity scale were made at the four external f positions in the heptad repeat. Equilibrium and kinetic circular dichroism studies demonstrate that both the stability and the unfolding and refolding rate constants vary with helix propensity but also reflect interactions of the altered side-chains with their local environments. Pairwise replacements and fragment studies show that the two C-terminal heptads are the likely source of the nucleating helices. Helix-helix recognition between preformed elements of secondary structure plays an important role in this fundamental folding reaction.     

Effects of the phenylalanine-22----leucine, glutamic acid-49----methionine, glycine-234----aspartic acid, and glycine-234----lysine mutations on the folding and stability of the alpha subunit of tryptophan synthase from Escherichia coli

The effects of four single amino acid replacements on the stability and folding of the alpha subunit of tryptophan synthase from Escherichia coli have been investigated by ultraviolet differences spectroscopy. In previous studies [Miles, E. W., Yutani, K., & Ogasahara, K. (1982) Biochemistry 21, 2586], it had been shown that the urea-induced unfolding at pH 7.8, 25 degrees C, proceeds by the initial unfolding of the less stable carboxyl domain (residues 189-268) followed by the unfolding of the more stable amino domain (residues 1-188). The effects of the Phe-22----Leu, Glu-49----Met, Gly-234----Asp, and Gly-234----Lys mutants on the equilibrium unfolding process can all be understood in terms of the domain unfolding model. With the exception of the Glu-49----Met replacement, the effects on stability are small. In contrast, the effects of three of the four mutations on the kinetics of interconversion of the native form and one of the stable partially folded intermediates are dramatic. The results for the Phe-22----Leu and Gly-234----Asp mutations indicate that these residues play a key role in the rate-limiting step. The Glu-49----Met mutation increases the stability of the native form with respect to that of the intermediate but does not affect the rate-limiting step. The Gly-234----Lys mutation does not affect either the stability or the kinetics of folding for the transition between native and intermediate forms. The changes in stability calculated from the unfolding and refolding rate constants agree quantitatively with those obtained from the equilibrium data. When considered with the results from a previous study on the Gly-211----Glu replacement [Matthews, C. R., Crisanti, M. M., Manz, J. T., & Gepner G. L. (1983) Biochemistry 22, 1445], it can be concluded that the rate-limiting step in the conversion of the intermediate to the native conformation involves either domain association or some other type of molecule-wide phenomenon.     

Pathways and kinetic barriers in mechanical unfolding and refolding of RNA and proteins

Using self-organized polymer models, we predict mechanical unfolding and refolding pathways of ribozymes, and the green fluorescent protein. In agreement with experiments, there are between six and eight unfolding transitions in the Tetrahymena ribozyme. Depending on the loading rate, the number of rips in the force-ramp unfolding of the Azoarcus ribozymes is between two and four. Force-quench refolding of the P4-P6 subdomain of the Tetrahymena ribozyme occurs through a compact intermediate. Subsequent formation of tertiary contacts between helices P5b-P6a and P5a/P5c-P4 leads to the native state. The force-quench refolding pathways agree with ensemble experiments. In the dominant unfolding route, the N-terminal alpha helix of GFP unravels first, followed by disruption of the N terminus beta strand. There is a third intermediate that involves disruption of three other strands. In accord with experiments, the force-quench refolding pathway of GFP is hierarchic, with the rate-limiting step being the closure of the barrel.     

Folding and unfolding of gammaTIM monomers and dimers

Kinetic simulations of the folding and unfolding of triosephosphate isomerase (TIM) from yeast were conducted using a single monomer gammaTIM polypeptide chain that folds as a monomer and two gammaTIM chains that fold to the native dimer structure. The basic protein model used was a minimalist Gō model using the native structure to determine attractive energies in the protein chain. For each simulation type--monomer unfolding, monomer refolding, dimer unfolding, and dimer refolding--thirty simulations were conducted, successfully capturing each reaction in full. Analysis of the simulations demonstrates four main conclusions. First, all four simulation types have a similar "folding order", i.e., they have similar structures in intermediate stages of folding between the unfolded and folded state. Second, despite this similarity, different intermediate stages are more or less populated in the four different simulations, with 1), no intermediates populated in monomer unfolding; 2), two intermediates populated with beta(2)-beta(4) and beta(1)-beta(5) regions folded in monomer refolding; 3), two intermediates populated with beta(2)-beta(3) and beta(2)-beta(4) regions folded in dimer unfolding; and 4), two intermediates populated with beta(1)-beta(5) and beta(1)-beta(5) + beta(6) + beta(7) + beta(8) regions folded in dimer refolding. Third, simulations demonstrate that dimer binding and unbinding can occur early in the folding process before complete monomer-chain folding. Fourth, excellent agreement is found between the simulations and MPAX (misincorporation proton alkyl exchange) experiments. In total, this agreement demonstrates that the computational Gō model is accurate for gammaTIM and that the energy landscape of gammaTIM appears funneled to the native state.     

Stabilisation of alpha-helices by site-directed mutagenesis reveals the importance of secondary structure in the transition state for acylphosphatase folding

The effects of stabilising mutations on the folding process of common-type acylphosphatase have been investigated. The mutations were designed to increase the helical propensity of the regions of the polypeptide chain corresponding to the two alpha-helices of the native protein. Various synthetic peptides incorporating the designed mutations were produced and their helical content estimated by circular dichroism. The most substantial increase in helical content is found for the peptide carrying five mutations in the second alpha-helix. Acylphosphatase variants containing the corresponding mutations display, to different extents, enhanced conformational stabilities as indicated by equilibrium urea denaturation experiments monitored by changes of intrinsic fluorescence. All the protein variants studied here refold with apparent two-state kinetics. Mutations in the first alpha-helix are responsible for a small increase in the refolding rate, accompanied by a marked decrease in the unfolding rate. On the other hand, multiple mutations in the second helix result in a considerable increase in the refolding rate without any significant effect on the unfolding rate. Addition of trifluoroethanol was found to accelerate the folding of the acylphosphatase variants, the extent of the acceleration being inversely proportional to the intrinsic rate of folding of the corresponding mutant. The trifluoroethanol-induced acceleration is far less marked for those variants whose alpha-helical structure is efficiently stabilised by amino acid replacements. This observation suggests that trifluoroethanol acts in a similar manner to the stabilising mutations in promoting native-like secondary structure. Analysis of the kinetic data indicates that the second helix is fully consolidated in the transition state for folding of acylphosphatase, whereas the first helix is only partially formed. These data suggest that the second helix is an important element in the folding process of the protein.     

Unfolding/refolding studies of smooth muscle tropomyosin. Evidence for a chain exchange mechanism in the preferential assembly of the native heterodimer

The thermal and the urea-induced unfolding profiles of the coiled-coil alpha-helix of native and refolded tropomyosin from chicken gizzard were studied by circular dichroism. Refolding of tropomyosin at low temperature from alpha + beta subunits, dissociated by guanidinium chloride, urea, or high temperature, predominantly produced alpha alpha + beta beta homodimers in agreement with earlier studies of refolding from guanidinium chloride (Graceffa, P. (1989) Biochemistry 28, 1282-1287). The presence of two unfolding transitions in low salt solutions with about equal helix loss verified the composition with the first unfolding transition of the homodimer mixture originating from alpha alpha. In contrast, refolding by equilibrating at temperatures close to physiological, however, produced the native alpha beta heterodimer, which unfolded in a single transition. The refolding kinetics of dissociated alpha + beta subunits indicated that beta beta homodimers form first, leading to alpha alpha homodimers both of which are relatively stable against chain exchange below approximately 25 degrees C. Equilibrating the homodimer mixture at 37-40 degrees C for long times, however, produced the native alpha beta molecule via chain exchange. The equilibria involved indicate that the free energy of formation from subunits of alpha beta is much less than that of (alpha alpha + beta beta)/2. In vivo folding of alpha beta from the two separate alpha and beta gene products is, therefore, thermodynamically favored over the formation of homodimers and biological factors need not be considered to explain the native preferred alpha beta composition.     

Characterization of the folding and unfolding reactions of single-chain monellin: evidence for multiple intermediates and competing pathways

The mechanisms of folding and unfolding of the small plant protein monellin have been delineated in detail. For this study, a single-chain variant of the natively two-chain monellin, MNEI, was used, in which the C terminus of chain B was connected to the N terminus of chain A by a Gly-Phe linker. Equilibrium guanidine hydrochloride (GdnHCl)-induced unfolding experiments failed to detect any partially folded intermediate that is stable enough to be populated at equilibrium to a significant extent. Kinetic experiments in which the refolding of GdnHCl-unfolded protein was monitored by measurement of the change in the intrinsic tryptophan fluorescence of the protein indicated the accumulation of three transient partially structured folding intermediates. The fluorescence change occurred in three kinetic phases: very fast, fast, and slow. It appears that the fast and slow changes in fluorescence occur on competing folding pathways originating from one unfolded form and that the very fast change in fluorescence occurs on a third parallel pathway originating from a second unfolded form of the protein. Kinetic experiments in which the refolding of alkali-unfolded protein was monitored by the change in the fluorescence of the hydrophobic dye 8-anilino-1-naphthalenesulfonic acid (ANS), consequent to the dye binding to the refolding protein, as well as by the change in intrinsic tryptophan fluorescence, not only confirmed the presence of the three kinetic intermediates but also indicated the accumulation of one or more early intermediates at a few milliseconds of refolding. These experiments also exposed a very slow kinetic phase of refolding, which was silent to any change in the intrinsic tryptophan fluorescence of the protein. Hence, the spectroscopic studies indicated that refolding of single-chain monellin occurs in five distinct kinetic phases. Double-jump, interrupted-folding experiments, in which the accumulation of folding intermediates and native protein during the folding process could be determined quantitatively by an unfolding assay, indicated that the fast phase of fluorescence change corresponds to the accumulation of two intermediates of differing stabilities on competing folding pathways. They also indicated that the very slow kinetic phase of refolding, identified by ANS binding, corresponds to the formation of native protein. Kinetic experiments in which the unfolding of native protein in GdnHCl was monitored by the change in intrinsic tryptophan fluorescence indicated that this change occurs in two kinetic phases. Double-jump, interrupted-unfolding experiments, in which the accumulation of unfolding intermediates and native protein during the unfolding process could be determined quantitatively by a refolding assay, indicated that the fast unfolding phase corresponds to the formation of fully unfolded protein via one unfolding pathway and that the slow unfolding phase corresponds to a separate unfolding pathway populated by partially unfolded intermediates. It is shown that the unfolded form produced by the fast unfolding pathway is the one which gives rise to the very fast folding pathway and that the unfolded form produced by the slower unfolding pathway is the one which gives rise to the slow and fast folding pathways.     

Helix mutations stabilize a late productive intermediate on the folding pathway of ubiquitin

We have investigated the relative placement of rate-limiting energy barriers and the role of productive or obstructive intermediates on the folding pathway of yeast wild-type ubiquitin ( wt-Ub) containing the F45W mutation. To manipulate the folding barriers, we have designed a family of mutants in which stabilizing substitutions have been introduced incrementally on the solvent-exposed surface of the main alpha-helix (residues 23-34), which has a low intrinsic helical propensity in the native sequence. Although the U --> I and I --> N transitions are not clearly delineated in the kinetics of wt-Ub, we show that an intermediate becomes highly populated and more clearly resolved as the predicted stability of the helix increases. The observed acceleration in the rate of folding correlates with helix stability and is consistent with the I-state representing a productive rather than misfolded state. A Leffler analysis of the effects on kinetics of changes in stability within the family of helix mutants results in a biphasic correlation in both the refolding and unfolding rates that suggest a shift from a nucleation-condensation mechanism (weakly stabilized helix) toward a diffusion-collision model (highly stabilized helix). Through the introduction of helix-stabilizing mutations, we are able to engineer a well-resolved I-state on the folding pathway of ubiquitin which is likely to be structurally distinct from that which is only weakly populated on the folding pathway of wild-type ubiquitin.     

Role of diffusion in the folding of the alpha subunit of tryptophan synthase from Escherichia coli

The rate-limiting step in the folding of the alpha subunit of tryptophan synthase has been proposed to be the association of two folding units. To probe the role of diffusion in this rate-limiting step, the urea-induced unfolding and refolding of the protein was examined in the presence of a number of viscosity-enhancing agents. The analysis was simplified by studying the effect of these agents on folding unit dissociation, the rate-limiting unfolding reaction, and the reverse of the rate-limiting step in refolding. In the presence of ethylene glycol, the relaxation times for unfolding to the same final conditions increased with increasing concentration of the cosolvent. When the effects of the cosolvent on protein stability were taken into account, the rates were found to show a unitary linear dependence on the viscosity of the solution. Similar results were obtained with glycerol and low concentrations of glucose, demonstrating that the effect is general and not specific to any viscogenic agent. These results clearly demonstrate that the rate-limiting folding unit association/dissociation reaction in the alpha subunit of tryptophan synthase involves a diffusional process.     

Formation of on- and off-pathway intermediates in the folding kinetics of Azotobacter vinelandii apoflavodoxin

The folding kinetics of the 179-residue Azotobacter vinelandii apoflavodoxin, which has an alpha-beta parallel topology, have been followed by stopped-flow experiments monitored by fluorescence intensity and anisotropy. Single-jump and interrupted refolding experiments show that the refolding kinetics involve four processes yielding native molecules. Interrupted unfolding experiments show that the two slowest folding processes are due to Xaa-Pro peptide bond isomerization in unfolded apoflavodoxin. The denaturant dependence of the folding kinetics is complex. Under strongly unfolding conditions (>2.5 M GuHCl), single exponential kinetics are observed. The slope of the chevron plot changes between 3 and 5 M denaturant, and no additional unfolding process is observed. This reveals the presence of two consecutive transition states on a linear pathway that surround a high-energy on-pathway intermediate. Under refolding conditions, two processes are observed for the folding of apoflavodoxin molecules with native Xaa-Pro peptide bond conformations, which implies the population of an intermediate. The slowest of these two processes becomes faster with increasing denaturant concentration, meaning that an unfolding step is rate-limiting for folding of the majority of apoflavodoxin molecules. It is shown that the intermediate that populates during refolding is off-pathway. The experimental data obtained on apoflavodoxin folding are consistent with the linear folding mechanism I(off) <==> U <==> I(on) <== > N, the off-pathway intermediate being the molten globule one that also populates during equilibrium denaturation of apoflavodoxin. The presence of such on-pathway and off-pathway intermediates in the folding kinetics of alpha-beta parallel proteins is apparently governed by protein topology.     

Slow formation of aggregation‐resistant β‐sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding‐related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large β‐helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all‐β‐sheet protein allows detailed analysis of the formation of β‐sheet structure in larger proteins. Using a combination of fluorescence and far‐UV circular dichroism spectroscopy, we show that the pertactin β‐helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and β‐sheet‐rich topology, pertactin refolding is reversible and not complicated by off‐pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate‐limiting step. Furthermore, site‐specific labeling experiments indicate that the β‐helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, β‐sheet‐rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, β‐sheet‐rich refolding intermediates. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Unfolding and refolding of bovine beta-lactoglobulin monitored by hydrogen exchange measurements.

Bovine beta-lactoglobulin (beta-LG) is a widely studied protein belonging to the lipocalin family, whose structural characterisation has been reported by X-ray crystallography and NMR studies at physiological and acidic pH, respectively. Bovine beta-LG consists of nine antiparallel beta-sheets and a terminal alpha-helix segment. The beta-sheets form a calyx structure with a hydrophobic buried cluster conferring stability to the protein while a hydrophobic surface patch provides stabilising interactions between the barrel and the flanking terminal helix. Here, the stability and the folding properties of bovine beta-LG in the presence of a chemical denaturant is probed. The analysis of the NMR spectra recorded in aqueous solution with increasing amounts of urea revealed that the intensities of the backbone cross-peaks decrease upon increasing urea concentration, while their secondary shifts do not change significantly on going from 0 to 5 M urea, thus suggesting the presence of slow exchange between native and unfolded protein. Hydrogen exchange measurements at different urea concentrations were performed in order to evaluate the exchange rates of individual backbone amide protons. The opening reactions that determine protein exchange can be computed for the most slowly exchanging hydrogen atoms, and the measured exchange rates and the corresponding free energies can be expressed in terms of the equilibrium energetic for the global transition and the local fluctuations. Most of the residues converge to define a common isotherm identifying a unique cooperative folding unit, encompassing all the strands, except strand betaI, and the terminal region of the helix. The amides that do not join the same global unfolding isotherm are characterised by low DeltaGH20op and especially by low m values, indicating that they are already substantially exposed in the native state. A two-state unfolding model N <==> U is therefore proposed for this rather big protein, in agreement with CD data. Renaturation studies show that the unfolding mechanism is reversible up to 6 M urea and suggest a similar unfolding and refolding pathway. Present results are discussed in light of the hypothesis of an alpha-->beta transition proposed for bovine beta-LG refolding.     

Mechanism of phage P22 tailspike protein folding mutations.

Temperature-sensitive folding (tsf) and global-tsf-suppressor (su) point mutations affect the folding yields of the trimeric, thermostable phage P22 tailspike endorhamnosidase at elevated temperature, both in vivo and in vitro, but they have little effect on function and stability of the native folded protein. To delineate the mechanism by which these mutations modify the partitioning between productive folding and off-pathway aggregation, the kinetics of refolding after dilution from acid-urea solutions and the thermal stability of folding intermediates were analyzed. The study included five tsf mutations of varying severity, the two known su mutations, and four tsf/su double mutants. At low temperature (10 degrees C), subunit-folding rates, measured as an increase in fluorescence, were similar for wild-type and mutants. At 25 degrees C, however, tsf mutations reduced the rate of subunit folding. The su mutations increased this rate, when present in the tsf-mutant background, but had no effect in the wild-type background. Conversely, tsf mutations accelerated, and su mutations retarded the irreversible off-pathway reaction, as revealed by temperature down-shifts after varied times during refolding at high temperature (40 degrees C). The kinetic results are consistent with tsf mutations destabilizing and su mutations stabilizing an essential subunit folding intermediate. In accordance with this interpretation, tsf mutations decreased, and su mutations increased the temperature resistance of folding intermediates, as disclosed by temperature up-shifts during refolding at 25 degrees C. The stabilizing and destabilizing effects were most pronounced early during refolding. However, they were not limited to subunit-folding intermediates and were also observable during thermal unfolding of the native protein.     

Unfolding of hen egg lysozyme by molecular dynamics simulations at 300K: insight into the role of the interdomain interface

We present the results of two 1.2 ns molecular dynamics (MD) unfolding simulations on hen egg lysozyme in water at 300K, performed using a new procedure called PEDC (Path Exploration With Distance Constraints). This procedure allows exploration of low energy structures as a function of increasing RMSD from the native structure, and offers especially the possibility of extensive exploration of the conformational space during the initial unfolding stages. The two independent MD simulations gave similar chronology of unfolding events: disruption of the active site, kinking of helix C, partial unfolding of the three-stranded beta-sheet to a two-stranded sheet (during which the helices A, B, and D remain to a great extent native), and finally unfolding of the beta-domain and partial unfolding of the alpha-domain in which hydrophobic clusters persist. We show particularly that the loss of hydrophobic contacts between the beta-sheet turn residues Leu55 and Ile56 and the hydrobic patch of the alpha-domain destabilizes the beta-domain and leads to its unfolding, suggesting that the correct embedding of these residues in the alpha-beta interface may constitute the rate limiting step in folding. These results are in accord with experimental observations on the folding/unfolding behavior of hen egg lysozyme at room temperature. They would also explain the loss of stability and the tendency to aggregation observed for the mutant Leu55Thr, and the slow refolding kinetics observed in the analogous amyloidogenic variant of human lysozyme.     

Mechanistic studies of the folding of human lysozyme and the origin of amyloidogenic behavior in its disease-related variants.

The unfolding and refolding properties of human lysozyme and two amyloidogenic variants (Ile56Thr and Asp67His) have been studied by stopped-flow fluorescence and hydrogen exchange pulse labeling coupled with mass spectrometry. The unfolding of each protein in 5.4 M guanidine hydrochloride (GuHCl) is well described as a two-state process, but the rates of unfolding of the Ile56Thr variant and the Asp67His variant in 5.4 M GuHCl are ca. 30 and 160 times greater, respectively, than that of the wild type. The refolding of all three proteins in 0.54 M GuHCl at pH 5.0 proceeds through persistent intermediates, revealed by multistep kinetics in fluorescence experiments and by the detection of well-defined populations in quenched-flow hydrogen exchange experiments. These findings are consistent with a predominant mechanism for refolding of human lysozyme in which one of the structural domains (the alpha-domain) is formed in two distinct steps and is followed by the folding of the other domain (the beta-domain) prior to the assembly of the two domains to form the native structure. The refolding kinetics of the Asp67His variant are closely similar to those of the wild-type protein, consistent with the location of this mutation in an outer loop of the beta-domain which gains native structure only toward the end of the refolding process. By contrast, the Ile56Thr mutation is located at the base of the beta-domain and is involved in the domain interface. The refolding of the alpha-domain is unaffected by this substitution, but the latter has the effect of dramatically slowing the folding of the beta-domain and the final assembly of the native structure. These studies suggest that the amyloidogenic nature of the lysozyme variants arises from a decrease in the stability of the native fold relative to partially folded intermediates. The origin of this instability is different in the two variants, being caused in one case primarily by a reduction in the folding rate and in the other by an increase in the unfolding rate. In both cases this results in a low population of soluble partially folded species that can aggregate in a slow and controlled manner to form amyloid fibrils.     

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.     

Highly divergent dihydrofolate reductases conserve complex folding mechanisms.

To test the hypothesis that protein folding mechanisms are better conserved than amino acid sequences, the mechanisms for dihydrofolate reductases (DHFR) from human (hs), Escherichia coli (ec) and Lactobacillus casei (lc) were elucidated and compared using intrinsic Trp fluorescence and fluorescence-detected 8-anilino-1-naphthalenesulfonate (ANS) binding. The development of the native state was monitored using either methotrexate (absorbance at 380 nm) or NADPH (extrinsic fluorescence) binding. All three homologs displayed complex unfolding and refolding kinetic mechanisms that involved partially folded states and multiple energy barriers. Although the pairwise sequence identities are less than 30 %, folding to the native state occurs via parallel folding channels and involves two types of on-pathway kinetic intermediates for all three homologs. The first ensemble of kinetic intermediates, detected within a few milliseconds, has significant secondary structure and exposed hydrophobic cores. The second ensemble is obligatory and has native-like side-chain packing in a hydrophobic core; however, these intermediates are unable to bind active-site ligands. The formation of the ensemble of native states occurs via three channels for hsDHFR, and four channels for lcDHFR and ecDHFR. The binding of active-site ligands (methotrexate and NADPH) accompanies the rate-limiting formation of the native ensemble. The conservation of the fast, intermediate and slow-folding events for this complex alpha/beta motif provides convincing evidence for the hypothesis that evolutionarily related proteins achieve the same fold via similar pathways.