Computer simulation of the folding-unfolding transition of island-model proteins--folding pathway, transition process, and fluctuations

Monte Carlo computer simulations were performed to elucidate the dynamic aspects of the folding and unfolding transitions of island-model protein. Five different types of model proteins were designed, according to characteristics of backbone structure. The computer simulations clearly show that the unfolding and folding transitions are all-or-none processes between the N-and U-states. They are typical Poisson processes. From the Arrhenius plots of rate constants, the activation enthalpies of folding and unfolding were determined. In addition, the folding pathways were determined along the reaction coordinate. Formations of several local structures along a polypeptide chain are almost simultaneous, but the most probable time sequence of events exists at the moment of transition. That is the most probable folding pathway. The unfolding pathway was found to be just the reverse process of the most probable folding pathway. The relationship between the fluctuations in each equilibrium state and the transition process was considered. In contrast to the theory of absolute reaction rate, the transient states are widely distributed along the reaction coordinate. From analysis of the “transient process,” we tried to determine the critical states from which the transient process starts. As a result, we found that the unfolding transition occurs at the stage near the N-state. During the U-state, large joined blocks rarely appear, but they appear in the transient process towards the N-state. However, the “branch point” between the N- and U-states lies near the N-state, and joined blocks tend to unfold prior to passing over the branch point. We concluded that the stability of later folding intermediates is important for selection of the folding pathway, while preferential selection of an early folding intermediate is important in acceleration of the folding rate. The effects of intrachain cross-linking and peptide fragment binding on the rate constants were examined by using computer simulations of model proteins. In general, a small-sized loop formed by cross-linking accelerates the folding rate and a large-sized loop contributes much to the stabilization of the native conformation. We also found that peptide fragment binding contributes little to the acceleration of the folding rate of the residual protein.     

(Un)Folding Mechanisms of the FBP28 WW Domain in Explicit Solvent Revealed by Multiple Rare Event Simulation Methods

We report a numerical study of the (un)folding routes of the truncated FBP28 WW domain at ambient conditions using a combination of four advanced rare event molecular simulation techniques. We explore the free energy landscape of the native state, the unfolded state, and possible intermediates, with replica exchange molecular dynamics. Subsequent application of bias-exchange metadynamics yields three tentative unfolding pathways at room temperature. Using these paths to initiate a transition path sampling simulation reveals the existence of two major folding routes, differing in the formation order of the two main hairpins, and in hydrophobic side-chain interactions. Having established that the hairpin strand separation distances can act as reasonable reaction coordinates, we employ metadynamics to compute the unfolding barriers and find that the barrier with the lowest free energy corresponds with the most likely pathway found by transition path sampling. The unfolding barrier at 300 K is ∼17 k_BT ≈ 42 kJ/mol, in agreement with the experimental unfolding rate constant. This work shows that combining several powerful simulation techniques provides a more complete understanding of the kinetic mechanism of protein folding.     

Kinetics and Reaction Coordinates of the Reassembly of Protein Fragments Via Forward Flux Sampling

We studied the mechanism of the reassembly and folding process of two fragments of a split lattice protein by using forward flux sampling (FFS). Our results confirmed previous thermodynamics and kinetics analyses that suggested that the disruption of the critical core (of an unsplit protein that folds by a nucleation mechanism) plays a key role in the reassembly mechanism of the split system. For several split systems derived from a parent 48-mer model, we estimated the reaction coordinates in terms of collective variables by using the FFS least-square estimation method and found that the reassembly transition is best described by a combination of the total number of native contacts, the number of interchain native contacts, and the total conformational energy of the split system. We also analyzed the transition path ensemble obtained from FFS simulations using the estimated reaction coordinates as order parameters to identify the microscopic features that differentiate the reassembly of the different split systems studied. We found that in the fastest folding split system, a balanced distribution of the original-core amino acids (of the unsplit system) between protein fragments propitiates interchain interactions at early stages of the folding process. Only this system exhibits a different reassembly mechanism from that of the unsplit protein, involving the formation of a different folding nucleus. In the slowest folding system, the concentration of the folding nucleus in one fragment causes its early prefolding, whereas the second fragment tends to remain as a detached random coil. We also show that the reassembly rate can be either increased or decreased by tuning interchain cooperativeness via the introduction of a single point mutation that either strengthens or weakens one of the native interchain contacts (prevalent in the transition state ensemble).     

Mechanical unfolding of human telomere G-quadruplex DNA probed by integrated fluorescence and magnetic tweezers spectroscopy

Single-molecule techniques facilitate analysis of mechanical transitions within nucleic acids and proteins. Here, we describe an integrated fluorescence and magnetic tweezers instrument that permits detection of nanometer-scale DNA structural rearrangements together with the application of a wide range of stretching forces to individual DNA molecules. We have analyzed the force-dependent equilibrium and rate constants for telomere DNA G-quadruplex (GQ) folding and unfolding, and have determined the location of the transition state barrier along the well-defined DNA-stretching reaction coordinate. Our results reveal the mechanical unfolding pathway of the telomere DNA GQ is characterized by a short distance (<1 nm) to the transition state for the unfolding reaction. This mechanical unfolding response reflects a critical contribution of long-range interactions to the global stability of the GQ fold, and suggests that telomere-associated proteins need only disrupt a few base pairs to destabilize GQ structures. Comparison of the GQ unfolded state with a single-stranded polyT DNA revealed the unfolded GQ exhibits a compacted non-native conformation reminiscent of the protein molten globule. We expect the capacity to interrogate macromolecular structural transitions with high spatial resolution under conditions of low forces will have broad application in analyses of nucleic acid and protein folding.     

Kinetic pathways of beta-hairpin (un)folding in explicit solvent

We examine the dynamical (un)folding pathways of the C-terminal beta-hairpin of protein G-B1 at room temperature in explicit solvent, by employing transition path sampling algorithms. The path ensembles contain information on the folding kinetics, including solvent motion. We determine the transition state ensembles for the two main transitions: 1), the hydrophobic collapse; and 2), the backbone hydrogen bond formation. In both cases the transition state ensembles are characterized by a layer (1) or a strip (2) of water molecules in between the two hairpin strands, supporting the hypothesis of the solvent as lubricant in the folding process. The transition state ensembles do not correspond with saddle points in the equilibrium free-energy landscapes. The kinetic pathways are thus not completely determined by the free-energy landscape. This phenomenon can occur if the order parameters obey different timescales. Using the transition interface sampling technique, we calculate the rate constants for (un)folding and find them in reasonable agreement with experiments, thus supporting the validation of using all-atom force fields to study protein folding.     

Characterization of a slow folding reaction for the alpha subunit of tryptophan synthase

The equilibria and kinetics of urea-induced unfolding and refolding of the alpha subunit of tryptophan synthase of E. coli have been examined for their dependences on viscosity, pH, and temperature in order to investigate the properties of one of the rate-limiting steps, domain association. A viscosity enhancer, 0.58 M sucrose, was found to slow unfolding and accelerate refolding. This apparently anomalous result was shown to be due to the stabilizing effect of sucrose on the folding reaction. After accounting for this stabilization effect by using linear free-energy plots, the unfolding and refolding kinetics were found to have a viscosity dependence. A decrease in pH was found to stabilize the domain association reaction by increasing the refolding rate and decreasing the unfolding rate. This effect was accounted for by protonation of a single residue with a pK value of 8.8 in the native state and 7.1 in the intermediate, in which the two domains are not yet associated. The activation energy of unfolding is 4.8 kcal/mol, close to the diffusion limit. The negative activation entropy of unfolding, -47 cal/deg-mol, which controls this reaction, may result from ordering of solvent about the newly exposed domain interface of the transition state. These results may provide information on the types of noncovalent interactions involved in domain association and improve the ability to interpret the folding of mutants with single amino-acid substitutions at the interface.     

Kinetic analysis of folding and unfolding the 56 amino acid IgG-binding domain of streptococcal protein G.

The 56 amino acid B domain of protein G (GB) is a stable globular folding unit with no disulfide cross-links. The physical properties of GB offer extraordinary flexibility for evaluating the energetics of the folding reaction. The protein is monomeric and very soluble in both folded and unfolded forms. The folding reaction has been previously examined by differential scanning calorimetry (Alexander et al., 1992) and found to exhibit two-state unfolding behavior over a wide pH range with an unfolding transition near 90 degrees C (GB1) at neutral pH. Here, the kinetics of folding and unfolding two naturally occurring versions of GB have been measured using stopped-flow mixing methods and analyzed according to transition-state theory. GB contains no prolines, and the kinetics of folding and unfolding can be fit to a single, first-order rate constant over the temperature range of 5-35 degrees C. The major thermodynamic changes going from the unfolded state to the transition state are (1) a large decrease in heat capacity (delta Cp), indicating that the transition state is compact and solvent inaccessible relative to the unfolded state; (2) a large loss of entropy; and (3) a small increase in enthalpy. The most surprising feature of the folding of GB compared to that of previously studied proteins is that its folding approximates a rapid diffusion controlled process with little increase in enthalpy going from the unfolded to the transition state.     

Pressure-jump studies of the folding/unfolding of trp repressor.

The dimeric protein, trp apo-repressor of Escherichia coli has been subjected to high hydrostatic pressure under a variety of conditions, and the effects have been monitored by fluorescence spectroscopic and infra-red absorption techniques. Under conditions of micromolar protein concentration and low, non-denaturing concentrations of guanidinium hydrochloride (GuHCl), tryptophan and 8-anilino-1-naphthalene sulfonate (ANS) fluorescence detected high pressure profiles demonstrate that pressures below 3 kbar result in dissociation of the dimer to a monomeric species that presents no hydrophobic binding sites for ANS. The FTIR-detected high pressure profile obtained under significantly different solution conditions (30 mM trp repressor in absence of denaturant) exhibits a much smaller pressure dependence than the fluorescence detected profiles. The pressure-denatured form obtained under the FTIR conditions retains about 50 % alpha-helical structure. From this we conclude that the secondary structure present in the high pressure state achieved under the conditions of the fluorescence experiments is at least as disrupted as that achieved under FTIR conditions. Fluorescence-detected pressure-jump relaxation studies in the presence of non-denaturing concentrations of GuHCl reveal a positive activation volume for the association/folding reaction and a negative activation volume for dissociation/unfolding reaction, implicating dehydration as the rate-limiting step for association/folding and hydration as the rate-limiting step for unfolding. The GuHCl concentration dependence of the kinetic parameters place the transition state at least half-way along the reaction coordinate between the unfolded and folded states. The temperature dependence of the pressure-jump fluorescence-detected dissociation/unfolding reaction in the presence of non-denaturing GuHCl suggests that the curvature in the temperature dependence of the stability arises from non-Arrhenius behavior of the folding rate constant, consistent with a large decrease in heat capacity upon formation of the transition state from the unfolded state. The decrease in the equilibrium volume change for folding with increasing temperature (due to differences in thermal expansivity of the folded and unfolded states) arises from a decrease in the absolute value for the activation volume for unfolding, thus indicating that the thermal expansivity of the transition state is similar to that of the unfolded state.     

Kinetics and thermodynamics of the unfolding and refolding of the three-stranded alpha-helical coiled coil, Lpp-56

Temperature-induced reversible unfolding and refolding of the three-stranded alpha-helical coiled coil, Lpp-56, were studied by kinetic and thermodynamic methods, using CD spectroscopy, dynamic light scattering, and scanning calorimetry. It was found that both unfolding and refolding reactions of this protein in neutral solution in the presence of 100 mM NaCl are characterized by unusually slow kinetics, which permits detailed investigation of the mechanism of these reactions. Kinetic analyses show that the unfolding of this coiled coil represents a single-stage first-order reaction, while the refolding represents a single-stage third-order reaction. The activation enthalpy and entropy for unfolding do not depend noticeably on temperature and are both significantly greater than those for the folding reaction, which show a significant dependence on temperature. The activation heat capacity change for the unfolding reaction is close to zero, while it is quite significant for the folding reaction. The correlation between the activation and structural parameters obtained for the Lpp-56 coiled coil suggests that interhelical van der Waals interactions are disrupted in the transition state, which is nevertheless still compact, and water has not yet penetrated into the interface; the transition from the transient state to the unfolded state results in hydration of exposed apolar groups of the interface and the disruption of helices. The low propensity for the Lpp-56 strands to fold and associate is caused by the high number of charged groups at neutral pH. On one hand, these charges give rise to considerable repulsive forces destabilizing the helical conformation of the strands. On the other hand, they align the folded helices in parallel and in register so that the apolar sides face each other, and the oppositely charged groups may form salt links, which are important for the formation of the trimeric coiled coil. A decrease in pH, which eliminates the salt links, dramatically decreases the stability of Lpp-56; its structure becomes less rigid and unfolds much faster.     

Folding pathway dependence on energetic frustration and interaction heterogeneity for a three-dimensional hydrophobic protein model

Monte Carlo simulations of a hydrophobic protein model of 40 monomers in the cubic lattice are used to explore the effect of energetic frustration and interaction heterogeneity on its folding pathway. The folding pathway is described by the dependence of relevant conformational averages on an appropriate reaction coordinate, pfold, defined as the probability for a given conformation to reach the native structure before unfolding. We compare the energetically frustrated and heterogeneous hydrophobic potential, according to which individual monomers have a higher or lower tendency to form contacts unspecifically depending on their hydrophobicities, to an unfrustrated homogeneous Go-type potential with uniformly attractive native interactions and neutral non-native interactions (called Go1 in this study), and to an unfrustrated heterogeneous potential with neutral non-native interactions and native interactions having the same energy as the hydrophobic potential (called Go2 in this study). Folding kinetics are slowed down dramatically when energetic frustration increases, as expected and previously observed in a two-dimensional model. Contrary to our previous results in two dimensions, however, it appears that the folding pathway and transition state ensemble can be significantly dependent on the energy function used to stabilize the native structure. The sequence of events along the reaction coordinate, or the order along this coordinate in which different regions of the native conformation become structured, turns out to be similar for the hydrophobic and Go2 potentials, but with analogous events tending to occur at lower pfold values in the first case. In particular, the transition state obtained from the ensemble around pfold = 0.5 is more structured for the hydrophobic potential. For Go1, not only the transition state ensemble but the order of events itself is modified, suggesting that interaction heterogeneity, in addition to energetic frustration, can have significant effects on the folding mechanism, most likely by modifying the probability of different contacts in the unfolded state, the starting point for the folding reaction. Although based on a simple model, these results provide interesting insight into how sequence-dependent switching between folding pathways might occur in real proteins.     

Is protein unfolding the reverse of protein folding? A lattice simulation analysis.

Simulations and experiments that monitor protein unfolding under denaturing conditions are commonly employed to study the mechanism by which a protein folds to its native state in a physiological environment. Due to the differences in conditions and the complexity of the reaction, unfolding is not necessarily the reverse of folding. To assess the relevance of temperature initiated unfolding studies to the folding problem, we compare the folding and unfolding of a 125-residue protein model by Monte Carlo dynamics at two temperatures; the lower one corresponds to the range used in T -jump experiments and the higher one to the range used in unfolding simulations of all-atom models. The trajectories that lead from the native state to the denatured state at these elevated temperatures are less diverse than those observed in the folding simulations. At the lower temperature, the system unfolds through a mandatory intermediate that corresponds to a local free energy minimum. At the higher temperature, no such intermediate is observed, but a similar pathway is followed. The structures contributing to the unfolding pathways resemble most closely those that make up the "fast track" of folding. The transition state for unfolding at the lower temperature (above Tm) is determined and is found to be more structured than the transition state for folding below the melting temperature. This shift towards the native state is consistent with the Hammond postulate. The implications for unfolding simulations of higher resolution models and for unfolding experiments of proteins are discussed.     

Folding rates and low-entropy-loss routes of two-state proteins

We develop a simple model for computing the rates and routes of folding of two-state proteins from the contact maps of their native structures. The model is based on the graph-theoretical concept of effective contact order (ECO). The model predicts that proteins fold by "zipping up" in a sequence of small-loop-closure events, depending on the native chain fold. Using a simple equation, with a few physical rate parameters, we obtain a good correlation with the folding rates of 24 two-state folding proteins. The model rationalizes data from Phi-value analysis that have been interpreted in terms of delocalized or polarized transition states. This model indicates how much of protein folding may take place in parallel, not along a single reaction coordinate or with a single transition state.     

Temperature-dependent conformational transitions and hydrogen-bond dynamics of the elastin-like octapeptide GVG(VPGVG): a molecular-dynamics study

A joint experimental/theoretical investigation of the elastin-like octapeptide GVG(VPGVG) was carried out. In this article a comprehensive molecular-dynamics study of the temperature-dependent folding and unfolding of the octapeptide is presented. The current study, as well as its experimental counterpart (see companion article in this issue) find that this peptide undergoes an inverse temperature transition (ITT), leading to a folding at approximately 40-60 degrees C. In addition, an unfolding transition is identified at unusually high temperatures approaching the normal boiling point of water. Due to the small size of the system, two broad temperature regimes are found: the ITT regime at approximately 10-60 degrees C and the unfolding regime at approximately T > 60 degrees C, where the peptide has a maximum probability of being folded at T approximately 60 degrees C. A detailed molecular picture involving a thermodynamic order parameter, or reaction coordinate, for this process is presented along with a time-correlation function analysis of the hydrogen-bond dynamics within the peptide as well as between the peptide and solvating water molecules. Correlation with experimental evidence and ramifications on the properties of elastin are discussed.     

Folding kinetics of the all-beta-sheet protein human basic fibroblast growth factor, a structural homolog of interleukin-1beta

The refolding and unfolding kinetics of the all-beta-sheet protein human basic fibroblast growth factor (hFGF-2) were studied by fluorescence spectroscopy. The kinetics of the unfolding transition are monophasic. The refolding reaction at high and low guanidinium chloride (GdmCl) concentrations is best described by mono- and biphasic folding, respectively. Refolding and unfolding of hFGF-2 (155 amino acids) is very slow compared with other non-disulfide-bonded monomeric proteins of similar size. For example, the rate constant for unfolding at 4.5 mol.liter(-1) GdmCl is 0.006 s(-1), and the refolding rate constants at 0.4 mol.liter(-1) GdmCl are 0.01 s(-1) and 0.0009 s(-1) (15 degrees C, pH 7.0). A characterization of the thermodynamic nature of the folding process using transition state theory revealed that the slow refolding is almost exclusively controlled by entropic factors, namely the strong loss of conformational freedom during refolding. The rate of the slow unfolding kinetics is mainly (and at low denaturant concentrations exclusively) controlled by the large positive change in enthalpy. hFGF-2 shows similar slow folding kinetics to that of its structural homolog interleukin-1beta. Since both proteins show very little sequence identity, it is suggested that their slow folding kinetics are determined by the complex beta-sheet arrangement of the native molecules.     

pH dependent unfolding characteristics of DLC8 dimer: Residue level details from NMR

Environment dependence of folding and unfolding of a protein is central to its function. In the same vein, knowledge of pH dependence of stability and folding/unfolding is crucial for many biophysical equilibrium and kinetic studies designed to understand protein folding mechanisms. In the present study we investigated the guanidine induced unfolding transition of dynein light chain protein (DLC8), a cargo adaptor of the dynein complex in the pH range 7-10. It is observed that while the protein remains a dimer in the entire pH range, its stability is somewhat reduced at alkaline pH. Global unfolding features monitored using fluorescence spectroscopy revealed that the unfolding transition of DLC8 at pH 7 is best described by a three-state model, whereas, that at pH 10 is best described by a two-state model. Chemical shift perturbations due to pH change provided insights into the corresponding residue level structural perturbations in the DLC8 dimer. Likewise, backbone (15)N relaxation measurements threw light on the corresponding motional changes in the dimeric protein. These observations have been rationalized on the basis of expected changes with increasing pH in the protonation states of the titratable residues on the structure of the protein. These, in turn provide an explanation for the change from three-state to two-state guanidine induced unfolding transition as the pH is increased from 7 to 10. All these results exemplify and highlight the role of environment vis-à-vis the sequence and structure of a given protein in dictating its folding/unfolding characteristics.     

Unfolding of the cold shock protein studied with biased molecular dynamics

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.     

How do chemical denaturants affect the mechanical folding and unfolding of proteins?

We present the first single-molecule atomic force microscopy study on the effect of chemical denaturants on the mechanical folding/unfolding kinetics of a small protein GB1 (the B1 immunoglobulin-binding domain of protein G from Streptococcus). Upon increasing the concentration of the chemical denaturant guanidinium chloride (GdmCl), we observed a systematic decrease in the mechanical stability of GB1, indicating the softening effect of the chemical denaturant on the mechanical stability of proteins. This mechanical softening effect originates from the reduced free-energy barrier between the folded state and the unfolding transition state, which decreases linearly as a function of the denaturant concentration. Chemical denaturants, however, do not alter the mechanical unfolding pathway or shift the position of the transition state for mechanical unfolding. We also found that the folding rate constant of GB1 is slowed down by GdmCl in mechanical folding experiments. By combining the mechanical folding/unfolding kinetics of GB1 in GdmCl solution, we developed the “mechanical chevron plot” as a general tool to understand how chemical denaturants influence the mechanical folding/unfolding kinetics and free-energy diagram in a quantitative fashion. This study demonstrates great potential in combining chemical denaturation with single-molecule atomic force microscopy techniques to reveal invaluable information on the energy landscape underlying protein folding/unfolding reactions.     

Detection and structure determination of an equilibrium unfolding intermediate of Rd-apocytochrome b562: native fold with non-native hydrophobic interactions

The absence of detectable kinetic and equilibrium folding intermediates by optical probes is commonly taken to indicate that protein folding is a two-state process. However, for some small proteins with apparent two-state behavior, unfolding intermediates have been identified in native-state hydrogen exchange or kinetic unfolding experiments monitored by nuclear magnetic resonance. Rd-apocytochrome b(562), a four-helix bundle, is one such protein. Here, we found another unfolding intermediate for Rd-apocytochrome b(562). It is based on a cooperative transition of (15)N chemical shifts of amide protons as a function of urea concentrations before the global unfolding. We have solved the high-resolution structure of the protein at 2.8 M urea, which is after this cooperative transition but before the global unfolding. All four helices remained intact, but a number of hydrophobic core residues repacked. This intermediate provides a possible structural interpretation for the kinetic unfolding intermediates observed using nuclear magnetic resonance methods for several proteins and has important implications for theoretical studies of protein folding.     

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.     

Circular permutation of 5-aminolevulinate synthase: effect on folding,conformational stability,and structure

The first and regulatory step of heme biosynthesis in mammals begins with the pyridoxal 5'-phosphate-dependent condensation reaction catalyzed by 5-aminolevulinate synthase. The enzyme functions as a homodimer with the two active sites at the dimer interface. Previous studies demonstrated that circular permutation of 5-aminolevulinate synthase does not prevent folding of the polypeptide chain into a structure amenable to binding of the pyridoxal 5'-phosphate cofactor and assembly of the two subunits into a functional enzyme. However, while maintaining a wild type-like three-dimensional structure, active, circularly permuted 5-aminolevulinate synthase variants possess different topologies. To assess whether the aminolevulinate synthase overall structure can be reached through alternative or multiple folding pathways, we investigated the guanidine hydrochloride-induced unfolding, conformational stability, and structure of active, circularly permuted variants in relation to those of the wild type enzyme using fluorescence, circular dichroism, activity, and size exclusion chromatography. Aminolevulinate synthase and circularly permuted variants folded reversibly; the equilibrium unfolding/refolding profiles were biphasic and, in all but one case, protein concentration-independent, indicating a unimolecular process with the presence of at least one stable intermediate. The formation of this intermediate was preceded by the disruption of the dimeric interface or dissociation of the dimer without significant change in the secondary structural content of the subunits. In contrast to the similar stabilities associated with the dimeric interface, the energy for the unfolding of the intermediate as well as the overall conformational stabilities varied among aminolevulinate synthase and variants. The unfolding of one functional permuted variant was protein concentration-dependent and had a potentially different folding mechanism. We propose that the order of the ALAS secondary structure elements does not determine the ability of the polypeptide chain to fold but does affect its folding mechanism.