Noninteracting local-structure model of folding and unfolding transition in globular proteins. I. Formulation

A statistical-mechanical model (a noninteracting local structure model) of folding and unfolding transition in globular proteins is described and a formulation is given to calculate the partition function. The process of transition is discussed in this model within the framework of equilibrium statistical mechanics. In order to clarify the range of applicability of such an approach, the characteristics of the folding and unfolding transition in globular proteins are analyzed from the statistical-physical point of view. A theoretical advantage is pointed out in studying folding and unfolding processes taking place as conformational fluctuations in individual protein molecules under macroscopic equilibrium at the melting temperature. In this case, paths of folding and unfolding are shown to be identical in the statistical sense. A key to the noninteracting local structure model lies in the concept of local structures and the assumption of the absence of interactions between local structures. A local structure is defined as a continuous section of the chain which takes the same or similar local conformation as in the native conformation. The assumption of the absence of inter-actions between local structures endows the model with the remarkable character that its partition function can be calculated exactly; thereby the equilibrium population of various conformations along the folding and unfolding paths can be discussed only by a knowledge of the folded native conformation.     

The effect of amino acid substitution on protein-folding and -unfolding transition studied by computer simulation

Protein-folding and -unfolding transitions were studied by the method of computer simulation. The protein was modeled as a two-dimensional lattice polymer. Various energy terms were assumed to be operative between units composing the polymer. But hydrophobic interactions were neglected explicitly. Both thermodynamic and kinetic quantities were obtained from the simulation, and from their temperature dependence in the transition zone characteristics of the conformational transition of proteins were discussed. Two amino acid substituted models, differing in the location of substitution, were studied and compared with the original in order to clarify the effect of substitution on conformational transition of proteins. The following conclusions were reached in this study: (1) The relaxation time of the slow mode, which reflects the overall folding and unfolding processes, shows a peak near the transition temperature, while that of the fast mode is almost independent of temperature. The peak of the slow mode occurs at a slightly lower temperature than the transition temperature. (2) The dependence of the logarithm of the rate constants on the inverse of temperature (Arrhenius plot) is linear. Therefore, the plot of the free energy of activation vs temperature is linear. (3) The values of kinetic parameters obtained suggest that in the activated state the intramolecular interactions are half broken, while the state is close to the native state on the entropy axis. (4) The amino acid substitution, which is modeled as having slightly unfavorable short-range interactions, causes the substituted ones to be slightly unstable. Moreover, it causes the folding transition to slow. From the analysis of the way slowing down is observed in the two substituted models, we conclude that a structure, designed to model a β-sheet, is formed before it gets assembled with other structures, which are designed to model α-helices. The process of assembly occurs nearly at the activated state of the folding and unfolding transition. (5) It is suggested from this study that the maximum of folding rate constant in the Arrhenius plot that has been observed experimentally in real proteins is likely due to hydrophobic interactions.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.     

A measure of conformational entropy change during thermal protein unfolding using neutron spectroscopy

Thermal unfolding of proteins at high temperatures is caused by a strong increase of the entropy change which lowers Gibbs free energy change of the unfolding transition (DeltaG(unf) = DeltaH - TDeltaS). The main contributions to entropy are the conformational entropy of the polypeptide chain itself and ordering of water molecules around hydrophobic side chains of the protein. To elucidate the role of conformational entropy upon thermal unfolding in more detail, conformational dynamics in the time regime of picoseconds was investigated with neutron spectroscopy. Confined internal structural fluctuations were analyzed for alpha-amylase in the folded and the unfolded state as a function of temperature. A strong difference in structural fluctuations between the folded and the unfolded state was observed at 30 degrees C, which increased even more with rising temperatures. A simple analytical model was used to quantify the differences of the conformational space explored by the observed protein dynamics for the folded and unfolded state. Conformational entropy changes, calculated on the basis of the applied model, show a significant increase upon heating. In contrast to indirect estimates, which proposed a temperature independent conformational entropy change, the measurements presented here, demonstrated that the conformational entropy change increases with rising temperature and therefore contributes to thermal unfolding.     

Essential dynamics of reversible peptide folding: memory-free conformational dynamics governed by internal hydrogen bonds

A principal component analysis has been applied on equilibrium simulations of a beta-heptapeptide that shows reversible folding in a methanol solution. The analysis shows that the configurational space contains only three dense sub-states. These states of relatively low free energy correspond to the "native" left-handed helix, a partly helical intermediate, and a hairpin-like structure. The collection of unfolded conformations form a relatively diffuse cloud with little substructure. Internal hydrogen-bonding energies were found to correlate well with the degree of folding. The native helical structure folds from the N terminus; the transition from the major folding intermediate to the native helical structure involves the formation of the two most C-terminal backbone hydrogen bonds. A four-state Markov model was found to describe transition frequencies between the conformational states within error limits, indicating that memory-effects are negligible beyond the nanosecond time-scale. The dominant native state fluctuations were found to be very similar to unfolding motions, suggesting that unfolding pathways can be inferred from fluctuations in the native state. The low-dimensional essential subspace, describing 69% of the collective atomic fluctuations, was found to converge at time-scales of the order of one nanosecond at all temperatures investigated, whereas folding/unfolding takes place at significantly longer time-scales, even above the melting temperature.     

Residue specific resolution of protein folding dynamics using isotope-edited infrared temperature jump spectroscopy

A major difficulty in experimental studies of protein folding is the lack of nonperturbing, residue specific probes of folding. Here, we demonstrate the ability to resolve protein folding dynamics at the level of a single residue using 13C=18O isotope-edited infrared spectroscopy. A single 13C=18O isotopic label was incorporated into the backbone of the 36 residue, three-helix bundle villin headpiece subdomain (HP36). The label was placed in a solvent protected region of the second alpha-helix of the protein. The 13C=18O isotopic label shifted the carbonyl stretching frequency to 1572.1 cm-1 in the folded state, well removed from the 12C=16O band of the unlabeled protein backbone. The unique IR signature of the 13C=18O label was exploited to probe the equilibrium thermal unfolding transition using temperature-dependent FTIR spectroscopy. The folding/unfolding dynamics were monitored using temperature-jump (T-jump) IR spectroscopy. The equilibrium unfolding studies showed conformational changes suggestive of a loss of helical structure in helix 2 prior to the global unfolding of the protein. T-jump relaxation kinetics probing both the labeled site and the 12C=16O band were found to be biphasic with similar relaxation rates. The slow relaxation phase (approximately 2 x 10(5) s-1) corresponds to the global folding transition. The location of the label, a buried position in helix 2, provides an important probe of the origin of the fast relaxation phase (approximately 10(7) s-1). This phase has significant amplitude for the labeled position even though it is well protected from solvent in the folded structure. The fast phase likely represents a rapid pre-equilibrium that involves solvent penetration around the label and possible partial unfolding of helix 2 prior to the global unfolding transition. This work represents the first experimental study of ultrafast folding dynamics with residue specific resolution.     

Molecular basis of co-operativity in protein folding.

The folding/unfolding transition of proteins is a highly co-operative process characterized by the presence of very few or no thermodynamically stable partially folded intermediate states. The purpose of this paper is to present a thermodynamic formalism aimed at describing quantitatively the co-operative folding behavior of proteins. In order to account for this behavior, a hierarchical algorithm aimed at evaluating the folding/unfolding partition function has been developed. This formalism defines the partition function in terms of multiple levels of interacting co-operative folding units. A co-operative folding unit is defined as a protein structural element that exhibits two-state folding/unfolding behavior. At the most fundamental level are those structural elements that behave co-operatively as a result of purely local interactions. Higher-order co-operative folding units are formed through interactions between different structural elements. The hierarchical formalism utilizes the crystallographic structure of the protein as a template to generate partially folded conformations defined in terms of co-operative folding units. The Gibbs free energy of those states and their corresponding statistical weights are then computed using experimental energetic parameters determined calorimetrically. This formalism has been applied to the case of myoglobin. It is shown that the hierarchical partition function correctly predicts the presence, energetics and co-operativity of the heat and cold denaturation transitions. The major contribution to the co-operative folding behavior arises from the solvent exposure of non-polar residues located in regions complementary to those that have undergone unfolding. This entropically uncompensated and energetically unfavorable solvent exposure characterizes all partially folded states but not the unfolded state, thus minimizing the population of partially folded intermediates throughout the folding/unfolding transition.     

Folding and stability of trp aporepressor from Escherichia coli

Equilibrium and kinetic studies of the urea-induced unfolding of trp aporepressor from Escherichia coli were performed to probe the folding mechanism of this intertwined, dimeric protein. The equilibrium unfolding transitions at pH 7.6 and 25 degrees C monitored by difference absorbance, fluorescence, and circular dichroism spectroscopy are coincident within experimental error. All three transitions are well described by a two-state model involving the native dimer and the unfolded monomer; the free energy of folding in the absence of denaturant and under standard-state conditions is estimated to be 23.3 +/- 0.9 kcal/mol of dimer. The midpoint of the equilibrium unfolding transition increases with increasing protein concentration in the manner expected from the law of mass action for the two-state model. We find no evidence for stable folding intermediates. Kinetic studies reveal that unfolding is governed by a single first-order reaction whose relaxation time decreases exponentially with increasing urea concentration and also decreases with increasing protein concentration in the transition zone. Refolding involves at least three phases that depend on both the protein concentration and the final urea concentration in a complex manner. The relaxation time of the slowest of these refolding phases is identical with that for the single phase in unfolding in the transition zone, consistent with the results expected for a reaction that is kinetically reversible. The two faster refolding phases are presumed to arise from slow isomerization reactions in the unfolded form and reflect parallel folding channels.     

Multifrequency calorimetry of the folding/unfolding transition of cytochrome c.

The folding-unfolding transition of Fe(III) cytochrome c has been studied with the new technique of multifrequency calorimetry. Multifrequency calorimetry is aimed at measuring directly the dynamics of the energetic events that take place during a thermally induced transition by measuring the frequency dispersion of the heat capacity. This is done by modulating the folding/unfolding equilibrium using a variable frequency, small oscillatory temperature perturbation (approximately 0.05-0.1 degrees C) centered at the equilibrium temperature of the system. Fe(III) cytochrome c at pH 4 undergoes a fully reversible folding/unfolding transition centered at 67.7 degrees C and characterized by an enthalpy change of 81 kcal/mol and heat capacity difference between unfolded and folded states of 0.9 kcal/K*mol. By measuring the temperature dependence of the frequency dispersion of the heat capacity in the frequency range of 0.1-1 Hz it has been possible to examine the time regime of the enthalpic events associated with the transition. The multifrequency calorimetry results indicate that approximately 85% of the excess heat capacity associated with the folding/unfolding transition relaxes with a single relaxation time of 326 +/- 68 ms at the midpoint of the transition region. This is the first time that the time regime in which heat is absorbed and released during protein folding/unfolding has been measured.     

The protein import motor of mitochondria: unfolding and trapping of preproteins are distinct and separable functions of matrix Hsp70.

Mitochondrial heat shock protein 70 (mtHsp70) functions in unfolding, translocation, and folding of imported proteins. Controversial models of mtHsp70 action have been discussed: (1) physical trapping of preproteins is sufficient to explain the various mtHsp70 functions, and (2) unfolding of preproteins requires an active motor function of mtHsp70 ("pulling"). Intragenic suppressors of a mutant mtHsp70 separate two functions: a nonlethal folding defect caused by enhanced trapping of preproteins, and a conditionally lethal unfolding defect caused by an impaired interaction of mtHsp70 with the membrane anchor Tim44. Even enhanced trapping in wild-type mitochondria does not generate a pulling force. The motor function of mtHsp70 cannot be explained by passive trapping alone but includes an essential ATP-dependent interaction with Tim44 to generate a pulling force and unfold preproteins.     

Multiple interactions of a DNA-binding protein in vivo. II. Effects of host mutations on DNA replication of phage T4 gene 32 mutants.

A quantitative analysis has been made of the kinetics of disulphide bond formation, breakage, and rearrangement which occur during the folding and unfolding of the pancreatic trypsin inhibitor. The results have been used to infer the energetics of the protein conformational transitions which accompany each step.The folding transition is shown to be a co-operative process in which all intermediate states with one or two disulphide bonds are unstable relative to the unfolded, reduced protein and that in the fully folded conformation with three disulphide bonds. The approximate two-state nature of the transition at equilibrium was demonstrated experimentally. The folding transition of the inhibitor which has been determined kinetically is therefore analogous to that observed generally with other globular proteins.     

Folding of Escherichia coli DsbC: characterization of a monomeric folding intermediate

The homodimeric protein DsbC is a disulfide isomerase and a chaperone located in the periplasm of Escherichia coli. We have studied the guanidine hydrochloride (GdnHCl)-induced unfolding and refolding of DsbC using mutagenesis, intrinsic fluorescence, circular dichroism spectra, size-exclusion chromatography, and sedimentation velocity analysis. The equilibrium refolding and unfolding of DsbC was thermodynamically reversible. The equilibrium folding profile measured by fluorescence excited at 280 nm exhibited a three-state transition profile with a stable folding intermediate formed at 0-2.0 M GdnHCl followed by a second transition at higher GdnHCl concentrations. Sedimentation velocity data revealed dissociation of the dimer to the monomer over the concentration range of the first transition (0-2.0 M). In contrast, fluorescence emission data for DsbC excited at 295 nm showed a single two-state transition. Fluorescence emission data for the equilibrium unfolding of the monomeric G49R mutant, excited at either 295 or 280 nm, indicated a single two-state transition. Data obtained for the dimeric Y52W mutant indicated a strong protein concentration dependence of the first transition but no dependence of the second transition in equilibrium unfolding. This suggests that the fluorescence of Y52W sensitively reports conformational changes caused by dissociation of the dimer. Thus, the folding of DsbC follows a three-state transition model with a monomeric folding intermediate formed in 0-2.0 M GdnHCl. The folding of DsbC in the presence of DTT indicates an important role for the non-active site disulfide bond in stabilizing the conformation of the molecule. Dimerization ensures the performance of chaperone and isomerase functions of DsbC.     

Conformational transitions of adenylate kinase: switching by cracking

Conformational heterogeneity in proteins is known to often be the key to their function. We present a coarse grained model to explore the interplay between protein structure, folding and function which is applicable to allosteric or non-allosteric proteins. We employ the model to study the detailed mechanism of the reversible conformational transition of Adenylate Kinase (AKE) between the open to the closed conformation, a reaction that is crucial to the protein's catalytic function. We directly observe high strain energy which appears to be correlated with localized unfolding during the functional transition. This work also demonstrates that competing native interactions from the open and closed form can account for the large conformational transitions in AKE. We further characterize the conformational transitions with a new measure Phi(Func), and demonstrate that local unfolding may be due, in part, to competing intra-protein interactions.     

Disulfide bonds and protein folding

The applications of disulfide-bond chemistry to studies of protein folding, structure, and stability are reviewed and illustrated with bovine pancreatic ribonuclease A (RNase A). After surveying the general properties and advantages of disulfide-bond studies, we illustrate the mechanism of reductive unfolding with RNase A, and discuss its application to probing structural fluctuations in folded proteins. The oxidative folding of RNase A is then described, focusing on the role of structure formation in the regeneration of the native disulfide bonds. The development of structure and conformational order in the disulfide intermediates during oxidative folding is characterized. Partially folded disulfide species are not observed, indicating that disulfide-coupled folding is highly cooperative. Contrary to the predictions of "rugged funnel" models of protein folding, misfolded disulfide species are also not observed despite the potentially stabilizing effect of many nonnative disulfide bonds. The mechanism of regenerating the native disulfide bonds suggests an analogous scenario for conformational folding. Finally, engineered covalent cross-links may be used to assay for the association of protein segments in the folding transition state, as illustrated with RNase A.     

Molecular mechanisms of pH-driven conformational transitions of proteins: insights from continuum electrostatics calculations of acid unfolding

The acid unfolding of staphylococcal nuclease (SNase) is very cooperative (Whitten and García-Moreno, Biochemistry 2000;39:14292-14304). As many as seven hydrogen ions (H+) are bound preferentially by the acid-unfolded state relative to the native (N) state in the pH range 3.2-3.9. To investigate the mechanism of acid unfolding, structure-based pKa calculations were performed with a variety of continuum electrostatic methods. The calculations reproduced successfully the H+ binding properties of the N state between pH 5 and 9, but they systematically overestimated the number of H+ bound upon acid unfolding. The calculated pKa values of all carboxylic residues in the N state were more depressed than they should be. The discrepancy between the observed and the calculated H+ uptake upon acid unfolding was not improved by using high protein dielectric constants, structures relaxed with molecular dynamics, or other empirical modifications implemented previously by others to maximize agreement between measured and calculated pKa values. This suggests an important role for conformational fluctuations of the backbone as important determinants of pKa values of carboxylic groups. Because no global or subglobal conformational changes have been observed previously for SNase under acidic conditions above the acid-unfolding region, these fluctuations must be local. The acid unfolding of SNase does not seem to involve the disruption of the N state by accruement of intramolecular repulsive interactions, nor the protonation of key ion paired carboxylic residues. It is more consistent with modest contributions from many H+ binding groups, with an important role for local conformational fluctuations in the coupling between H+ binding and the global structural transition.     

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.     

Pressure-jump small-angle x-ray scattering detected kinetics of staphylococcal nuclease folding

The kinetics of chain disruption and collapse of staphylococcal nuclease after positive or negative pressure jumps was monitored by real-time small-angle x-ray scattering under pressure. We used this method to probe the overall conformation of the protein by measuring its radius of gyration and pair-distance-distribution function p(r) which are sensitive to the spatial extent and shape of the particle. At all pressures and temperatures tested, the relaxation profiles were well described by a single exponential function. No fast collapse was observed, indicating that the rate limiting step for chain collapse is the same as that for secondary and tertiary structure formation. Whereas refolding at low pressures occurred in a few seconds, at high pressures the relaxation was quite slow, approximately 1 h, due to a large positive activation volume for the rate-limiting step for chain collapse. A large increase in the system volume upon folding implies significant dehydration of the transition state and a high degree of similarity in terms of the packing density between the native and transition states in this system. This study of the time-dependence of the tertiary structure in pressure-induced folding/unfolding reactions demonstrates that novel information about the nature of protein folding transitions and transition states can be obtained from a combination of small-angle x-ray scattering using high intensity synchrotron radiation with the high pressure perturbation technique.