Mutational analysis of the BPTI folding pathway: I. Effects of aromatic-->leucine substitutions on the distribution of folding intermediates.

The roles of aromatic residues in determining the folding pathway of bovine pancreatic trypsin inhibitor (BPTI) were analyzed mutationally by examining the distribution of disulfide-bonded intermediates that accumulated during the refolding of protein variants in which tyrosine or phenylalanine residues were individually replaced with leucine. The eight substitutions examined all caused significant changes in the intermediate distribution. In some cases, the major effect was to decrease the accumulation of intermediates containing two of the three disulfides found in the native protein, without affecting the distribution of earlier intermediates. Other substitutions, however, led to much more random distributions of the intermediates containing only one disulfide. These results indicate that the individual residues making up the hydrophobic core of the native protein make clearly distinguishable contributions to conformation and stability early in folding: The early distribution of intermediates does not appear to be determined by a general hydrophobic collapse. The effects of the substitutions were generally consistent with the structures of the major intermediates determined by NMR studies of analogs, confirming that the distribution of disulfide-bonded species is determined by stabilizing interactions within the ordered regions of the intermediates. The plasticity of the BPTI folding pathway implied by these results can be described using conformational funnels to illustrate the degree to which conformational entropy is lost at different stages in the folding of the wild-type and mutant proteins.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Predicting repeat protein folding kinetics from an experimentally determined folding energy landscape

The Notch ankyrin domain is a repeat protein whose folding has been characterized through equilibrium and kinetic measurements. In previous work, equilibrium folding free energies of truncated constructs were used to generate an experimentally determined folding energy landscape (Mello and Barrick, Proc Natl Acad Sci USA 2004;101:14102–14107). Here, this folding energy landscape is used to parameterize a kinetic model in which local transition probabilities between partly folded states are based on energy values from the landscape. The landscape‐based model correctly predicts highly diverse experimentally determined folding kinetics of the Notch ankyrin domain and sequence variants. These predictions include monophasic folding and biphasic unfolding, curvature in the unfolding limb of the chevron plot, population of a transient unfolding intermediate, relative folding rates of 19 variants spanning three orders of magnitude, and a change in the folding pathway that results from C‐terminal stabilization. These findings indicate that the folding pathway(s) of the Notch ankyrin domain are thermodynamically selected: the primary determinants of kinetic behavior can be simply deduced from the local stability of individual repeats.     

Pathway and stability of protein folding.

We describe an experimental approach to the problem of protein folding and stability which measures interaction energies and maps structures of intermediates and transition states during the folding pathway. The strategy is based on two steps. First, protein engineering is used to remove interactions that stabilize defined positions in barnase, the RNAse from Bacillus amyloliquefaciens. The consequent changes in stability are measured from the changes in free energy of unfolding of the protein. Second, each mutation is used as a probe of the structure around the wild-type side chain during the folding process. Kinetic measurements are made on the folding and unfolding of wild-type and mutant proteins. The kinetic and thermodynamic data are combined and analysed to show the role of individual side chains in the stabilization of the folded, transition and intermediate states of the protein. The protein engineering experiments are corroborated by nuclear magnetic resonance studies of hydrogen exchange during the folding process. Folding is a multiphasic process in which alpha-helices and beta-sheet are formed relatively early. Formation of the hydrophobic core by docking helix and sheet is (partly) rate determining. The final steps involve the forming of loops and the capping of the N-termini of helices.     

Conformation of P22 tailspike folding and aggregation intermediates probed by monoclonal antibodies.

The partitioning of partially folded polypeptide chains between correctly folded native states and off-pathway inclusion bodies is a critical reaction in biotechnology. Multimeric partially folded intermediates, representing early stages of the aggregation pathway for the P22 tailspike protein, have been trapped in the cold and isolated by nondenaturing polyacrylamide gel electrophoresis (PAGE) (speed MA, Wang DIC, King J. 1995. Protein Sci 4:900-908). Monoclonal antibodies against tailspike chains discriminate between folding intermediates and native states (Friguet B, Djavadi-Ohaniance L, King J, Goldberg ME. 1994. J Biol Chem 269:15945-15949). Here we describe a nondenaturing Western blot procedure to probe the conformation of productive folding intermediates and off-pathway aggregation intermediates. The aggregation intermediates displayed epitopes in common with productive folding intermediates but were not recognized by antibodies against native epitopes. The nonnative epitope on the folding and aggregation intermediates was located on the partially folded N-terminus, indicating that the N-terminus remained accessible and nonnative in the aggregated state. Antibodies against native epitopes blocked folding, but the monoclonal directed against the N-terminal epitope did not, indicating that the conformation of the N-terminus is not a key determinant of the productive folding and chain association pathway.     

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.     

Some recommendations for the practitioner to improve the precision of experimentally determined protein folding rates and phi values

The mechanism by which proteins fold from an initially random conformation into a functional, native structure remains a major unsolved question in molecular biology. Of particular interest to the protein folding community is the structure that the protein adopts in the folding transition state (the highest free energy state on the pathway from unfolded to folded), as that state forms the barrier that defines the folding pathway. Unfortunately, however, unlike those of the initial, unfolded state and the final, folded state of the protein, the structure in the transition state cannot be directly assessed via experiment. Instead, experimentalists infer the structure of the transition state, often by estimating changes in its free energy by measuring the effects of amino acid substitutions on folding and unfolding rates (Phi-value analysis). In this article we show how to obtain more efficient estimates of these important quantities via improved experimental designs, and how to avoid common pitfalls in the analysis of kinetic data during the extraction of these parameters.     

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.     

A study of intermediates involved in the folding pathway for recombinant human macrophage colony-stimulating factor (M-CSF): evidence for two distinct folding pathways.

The folding pathway for a 150-amino acid recombinant form of the dimeric cytokine human macrophage colony-stimulating factor (M-CSF) has been studied. All 14 cysteine residues in the biologically active homodimer are involved in disulfide linkages. The structural characteristics of folding intermediates blocked with iodoacetamide reveal a rapid formation of a small amount of a non-native dimeric intermediate species followed by a slow progression via both monomeric and dimeric intermediates to the native dimer. The transition from monomer to fully folded dimer is complete within 25 h at room temperature at pH 9.0. The blocked intermediates are stable under conditions of sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and thus represent various dimeric and folded monomeric species of the protein with different numbers of disulfide bridges. Peptide mapping and electrospray ionization mass spectrometry revealed that a folded monomeric species of M-CSF contained three of the four native disulfide bridges, and this folded monomer also showed some biological activity in a cell-based assay. The results presented here strongly suggest that M-CSF can fold via two different pathways, one involving monomeric intermediates and another involving only dimeric intermediates.     

On the role of some conserved and nonconserved amino acid residues in the transitional state and intermediate of apomyoglobin folding

The contributions of some amino acid residues in the A, B, G, and H helices to the formation of the folding nucleus and folding intermediate of apomyoglobin were estimated. The effects of point substitutions of Ala for hydrophobic amino acid residues on the structural stability of the native (N) protein and its folding intermediate (I), as well as on the folding/unfolding rates for four mutant apomyoglobin forms, were studied. The equilibrium and kinetic studies of the folding/unfolding rates of these mutant proteins in a wide range of urea concentrations demonstrated that their native state was considerably destabilized as compared with the wild-type protein, whereas the stability of the intermediate state changed moderately. It was shown that the amino acid residues in the A, G, and H helices contributed insignificantly to the stabilization of the apomyoglobin folding nucleus in the rate-limiting I ? N transition, taking place after the formation of the intermediate, whereas the residue of the B helix was of great importance in the formation of the folding nucleus in this transition.     

Local interactions and the optimization of protein folding

The role of local interactions in protein folding has recently been the subject of some controversy. Here we investigate an extension of Zwanzig's simple and general model of folding in which local and nonlocal interactions are represented by functions of single and multiple conformational degrees of freedom, respectively. The kinetics and thermodynamics of folding are studied for a series of energy functions in which the energy of the native structure is fixed, but the relative contributions of local and nonlocal interactions to this energy are varied over a broad range. For funnel shaped energy landscapes, we find that 1) the rate of folding increases, but the stability of the folded state decreases, as the contribution of local interactions to the energy of the native structure increases, and 2) the amount of native structure in the unfolded state and the transition state vary considerably with the local interaction strength. Simple exponential kinetics and a well-defined free energy barrier separating folded and unfolded states are observed when nonlocal interactions make an appreciable contribution to the energy of the native structure; in such cases a transition state theory type approximation yields reasonably accurate estimates of the folding rate. Bumps in the folding funnel near the native state, which could result from desolvation effects, side chain freezing, or the breaking of nonnative contacts, significantly alter the dependence of the folding rate on the local interaction strength: the rate of folding decreases when the local interaction strength is increased beyond a certain point. A survey of the distribution of strong contacts in the protein structure database suggests that evolutionary optimization has involved both kinetics and thermodynamics: strong contacts are enriched at both very short and very long sequence separations. Proteins 29:282–291, 1997. © 1997 Wiley-Liss, Inc.     

The nature of the free energy barriers to two-state folding

We introduce a simple procedure to analyze the temperature dependence of the folding and unfolding rates of two-state proteins. We start from the simple transition-state-like rate expression: k = D(eff)exp(-DeltaG(TS)/RT), in which upper and lower bounds for the intra-chain effective diffusion coefficient (D(eff)) are obtained empirically using the timescales of elementary processes in protein folding. From the changes in DeltaG(TS) as a function of temperature, we calculate enthalpies and heat capacities of activation, together with the more elusive entropies of activation. We then estimate the conformational entropy of the transition state by extrapolation to the temperature at which the solvation entropy vanishes by cancellation between polar and apolar terms. This approach is based on the convergence temperatures for the entropy of solvating apolar (approximately 385 K) and polar groups (approximately 335 K), the assumption that the structural properties of the transition state are somewhere in between the unfolded and folded states, and the established relationship between observed heat capacity and solvent accessibility.1 To circumvent the lack of structural information about transition states, we use the empirically determined heat capacities of activation as constraints to identify the extreme values of the transition state conformational entropy that are consistent with experiment. The application of this simple approach to six two-state folding proteins for which there is temperature-dependent data available in the literature provides important clues about protein folding. For these six proteins, we obtain an average equilibrium cost in conformational entropy of -4.3 cal x mol(-1)K(-1)per residue, which is in close agreement to previous empirical and computational estimates of the same quantity. Furthermore, we find that all these proteins have a conformationally diverse transition state, with more than half of the conformational entropy of the unfolded state. In agreement with predictions from theory and computer simulations, the transition state signals the change from a regime dominated by loss in conformational entropy to one driven by the gain in stabilization free energy (i.e., including protein interactions and solvation effects). Moreover, the height of the barrier is determined by how much stabilization free energy is realized at that point, which is related to the relative contribution of local versus non-local interactions. A remarkable observation is that the fraction of conformational entropy per residue that is present in the transition state is very similar for the six proteins in this study. Based on this commonality, we propose that the observed change in thermodynamic regime is connected to a change in the pattern of structure formation: from one driven by formation of pairwise interactions to one dominated by coupling of the networks of interactions involved in forming the protein core. In this framework, the barrier to two-state folding is crossed when the folding protein reaches a "critical native density" that allows expulsion of remaining interstitial water and consolidation of the core. The principle of critical native density should be general for all two-state proteins, but can accommodate different folding mechanisms depending on the particularities of the structure and sequence.     

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.     

Mutational analysis of acylphosphatase suggests the importance of topology and contact order in protein folding.

Muscle acylphosphatase (AcP) is a small protein that folds very slowly with two-state behavior. The conformational stability and the rates of folding and unfolding have been determined for a number of mutants of AcP in order to characterize the structure of the folding transition state. The results show that the transition state is an expanded version of the native protein, where most of the native interactions are partially established. The transition state of AcP turns out to be remarkably similar in structure to that of the activation domain of procarboxypeptidase A2 (ADA2h), a protein having the same overall topology but sharing only 13% sequence identity with AcP. This suggests that transition states are conserved between proteins with the same native fold. Comparison of the rates of folding of AcP and four other proteins with the same topology, including ADA2h, supports the concept that the average distance in sequence between interacting residues (that is, the contact order) is an important determinant of the rate of protein folding.     

Absence of a stable intermediate on the folding pathway of protein A.

The B-domain of protein A has one of the simplest protein topologies, a three-helix bundle. Its folding has been studied as a model for elementary steps in the folding of larger proteins. Earlier studies suggested that folding might occur by way of a helical hairpin intermediate. Equilibrium hydrogen exchange measurements indicate that the C-terminal helical hairpin could be a potential folding intermediate. Kinetic refolding experiments were performed using stopped-flow circular dichroism and NMR hydrogen-deuterium exchange pulse labeling. Folding of the entire molecule is essentially complete within the 6 ms dead time of the quench-flow apparatus, indicating that the intermediate, if formed, progresses rapidly to the final folded state. Site-directed mutagenesis of the isoleucine residue at position 16 was used to generate a variant protein containing tryptophan (the 116 W mutant). The formation of the putative folding intermediate was expected to be favored in this mutant at the expense of the native folded form, due to predicted unfavorable steric interactions of the bulky tryptophan side chain in the folded state. The 116 W mutant refolds completely within the dead time of a stopped-flow fluorescence experiment. No partly folded intermediate could be detected by either kinetic or equilibrium measurements. Studies of peptide fragments suggest that the protein A sequence has an intrinsic propensity to form a helix II/helix III hairpin. However, its stability appears to be marginal (of the order of 1/2 kT) and it could not be an obligatory intermediate on a defined folding pathway. These results explicitly demonstrate that the protein A B domain folds extremely rapidly by an apparent two-state mechanism without formation of stable partly folded intermediates. Similar mechanisms may also be involved in the rapid folding of subdomains of larger proteins to form the compact molten globule intermediates that often accumulate during the folding process.     

Very fast folding and association of a trimerization domain from bacteriophage T4 fibritin

The foldon domain constitutes the C-terminal 30 amino acid residues of the trimeric protein fibritin from bacteriophage T4. Its function is to promote folding and trimerization of fibritin. We investigated structure, stability and folding mechanism of the isolated foldon domain. The domain folds into the same trimeric beta-propeller structure as in fibritin and undergoes a two-state unfolding transition from folded trimer to unfolded monomers. The folding kinetics involve several consecutive reactions. Structure formation in the region of the single beta-hairpin of each monomer occurs on the submillisecond timescale. This reaction is followed by two consecutive association steps with rate constants of 1.9(+/-0.5)x10(6)M(-1)s(-1) and 5.4(+/-0.3)x10(6)M(-1)s(-1) at 0.58 M GdmCl, respectively. This is similar to the fastest reported bimolecular association reactions for folding of dimeric proteins. At low concentrations of protein, folding shows apparent third-order kinetics. At high concentrations of protein, the reaction becomes almost independent of protein concentrations with a half-time of about 3 ms, indicating that a first-order folding step from a partially folded trimer to the native protein (k=210 +/- 20 s(-1)) becomes rate-limiting. Our results suggest that all steps on the folding/trimerization pathway of the foldon domain are evolutionarily optimized for rapid and specific initiation of trimer formation during fibritin assembly. The results further show that beta-hairpins allow efficient and rapid protein-protein interactions during folding.     

Structure of a partially unfolded form of Escherichia coli dihydrofolate reductase provides insight into its folding pathway

Proteins frequently fold via folding intermediates that correspond to local minima on the conformational energy landscape. Probing the structure of the partially unfolded forms in equilibrium under native conditions can provide insight into the properties of folding intermediates. To elucidate the structures of folding intermediates of Escherichia coli dihydrofolate reductase (DHFR), we investigated transient partial unfolding of DHFR under native conditions. We probed the structure of a high‐energy conformation susceptible to proteolysis (cleavable form) using native‐state proteolysis. The free energy for unfolding to the cleavable form is clearly less than that for global unfolding. The dependence of the free energy on urea concentration (m‐value) also confirmed that the cleavable form is a partially unfolded form. By assessing the effect of mutations on the stability of the partially unfolded form, we found that native contacts in a hydrophobic cluster formed by the F‐G and Met‐20 loops on one face of the central β‐sheet are mostly lost in the partially unfolded form. Also, the folded region of the partially unfolded form is likely to have some degree of structural heterogeneity. The structure of the partially unfolded form is fully consistent with spectroscopic properties of the near‐native kinetic intermediate observed in previous folding studies of DHFR. The findings suggest that the last step of the folding of DHFR involves organization in the structure of two large loops, the F‐G and Met‐20 loops, which is coupled with compaction of the rest of the protein.     

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.     

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.