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.     

Destabilization of the Escherichia coli RNase H kinetic intermediate: switching between a two-state and three-state folding mechanism

Escherichia coli RNase H folds through a partially folded kinetic intermediate that mirrors a rarely populated, partially unfolded form detectable by native-state hydrogen exchange under equilibrium conditions. Residue 53 is at the interface of two helices known to be structured in this intermediate. Kinetic refolding studies on mutant proteins varying in size and hydrophobicity at residue 53 support a contribution of hydrophobicity to the stabilities of the kinetic intermediate and the transition state. Packing interactions also play a significant role in the stability of these two states, though they play a much larger role in the native-state stability. One dramatic mutation, I53D, results in the conversion from a three-state to a two-state folding mechanism, which is explained most easily through a simple destabilization of the kinetic intermediate such that it is no longer stable with respect to the unfolded state. These results demonstrate that interactions that stabilize an intermediate can accelerate folding if these same interactions are present in the transition state. Our results are consistent with a hierarchical model of folding, where the intermediate consists of native-like interactions, is on-pathway, and is productive for folding.     

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

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

The apomyoglobin folding pathway revisited: structural heterogeneity in the kinetic burst phase intermediate

Extensive analysis of accurate quench-flow hydrogen exchange results indicates that the burst phase kinetic intermediate in the folding of apomyoglobin (apoMb) from urea is structurally heterogeneous. The structural variability is associated with the partial folding of the E helix during the burst phase (<6.4ms) of the folding process. Analysis of the effects of exchange-out of amide proton labels during the labeling pulse ( approximately pH 10) of the quench-flow process indicates that three of the amide protons in the E helix are in fact largely protected in the burst phase of folding, while the remainder of the E helix has a substantial complement of amide protons that show biphasic kinetics, i.e. are protected partly during the burst phase and partly during the slow phase of folding. The locations of these amide protons can be used to map the sites of structural heterogeneity in the kinetic molten globule. These sites include, besides the E helix, the ends of the A and B helices and part of the C helix. Our results give significant support to the hypothesis that the kinetic molten globule intermediate of apoMb is native-like.     

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.     

Time-resolved small-angle X-ray scattering study of the folding dynamics of barnase

Structural changes of barnase during folding were investigated using time-resolved small-angle X-ray scattering (SAXS). The folding of barnase involves a burst-phase intermediate, sometimes designated as the denatured state under physiological conditions, D_phys, and a second hidden intermediate. Equilibrium SAXS measurements showed that the radius of gyration (R_g) of the guanidine unfolded state (U) is 26.9 ± 0.7 Å, which remains largely constant over a wide denaturant concentration range. Time-resolved SAXS measurements showed that the R_g value extrapolated from kinetic R_g data to time zero, R_g,0, is 24.3 ± 0.1 Å, which is smaller than that of U but which is expanded from that of folding intermediates of other proteins with similar chain lengths (19 Å). After the burst-phase change, a single-exponential reduction in R_g² was observed, which corresponds to the formation of the native state for the major component containing the native trans proline isomer. We estimated R_g of the minor component of D_phys containing the non-native cis proline isomer (D_phys,cis) to be 25.7 ± 0.6 Å. Moreover, R_g of the major component of D_phys containing the native proline isomer (D_phys,tra) was estimated as 23.9 ± 0.2 Å based on R_g,0. Consequently, both components of the burst-phase intermediate of barnase (D_phys,tra and D_phys,cis) are still largely expanded. It was inferred that D_phys possesses the N-terminal helix and the center of the β-sheet formed independently and that the formation of the remainder of the protein occurs in the slower phase.     

Structure and heterogeneity of the one- and two-disulfide folding intermediates of tick anticoagulant peptide

Tick anticoagulant peptide (TAP) is a factor Xa-specific inhibitor and is structurally homologous to bovine pancreatic trypsin inhibitor (BPTI). The fully reduced TAP refolds spontaneously to form the native structure under a wide variation of redox buffers. The folding intermediates of TAP consist of at least 22 fractions of one-disulfide, two-disulfide, and three-disulfide scrambled isomers. Three species of well-populated one- and two-disulfide intermediates were isolated and structurally characterized. The predominant one-disulfide species contains TAP-(Cys33—Cys55). Two major two-disulfide isomers were TAP-(Cys33—Cys55, Cys15—Cys39) and TAP-(Cys33—Cys55, Cys5—Cys39). Both Cys33—Cys55 and Cys15—Cys39 are native disulfides of TAP. These three species are structural counterparts of BPTI-(Cys30—Cys51), BPTI-(Cys30—Cys51, Cys14—Cys38), and BPTI-(Cys30—Cys51,Cys5—Cys38), which have been shown to be the major intermediates of BPTI folding. In addition, time-course-trapped folding intermediates of TAP, consisting of about 47% one-disulfide species and 30% two-disulfide species, were collectively digested with thermolysin, and fragmented peptides were analyzed by Edman sequencing and mass spectrometry in order to characterize the disulfide-containing peptides. Among the 15 possible single-disulfide pairings of TAP, 10 (2 native and 8 nonnative) were found as structural components of its one- and two-disulfide folding intermediates. The results demonstrate that the major folding intermediates of TAP bear structural homology to those of BPTI. However, the folding pathway of TAP differs from that of BPTI by (a) a higher degree of heterogeneity of one- and two-disulfide intermediates and (b) the presence of three-disulfide scrambled isomers as folding intermediates. Mechanism(s) that may account for these diversities are proposed and discussed.     

Fast folding of a prototypic polypeptide: the immunoglobulin binding domain of streptococcal protein G.

The folding of the small (56 residues) highly stable B1 immunoglobulin binding domain (GB1) of streptococcal protein G has been investigated by quenched-flow deuterium-hydrogen exchange. This system represents a paradigm for the study of protein folding because it exhibits no complicating features superimposed upon the intrinsic properties of the polypeptide chain. Collapse to a semicompact state exhibiting partial order, reflected in protection factors for ND-NH exchange up to 10-fold higher than that expected for a random coil, occurs within the dead time (< or = 1 ms) of the quenched flow apparatus. This is followed by the formation of the fully native state, as monitored by the fractional proton occupancy of 26 backbone amide groups spread throughout the protein, in a single rapid concerted step with a half-life of 5.2 ms at 5 degrees C.     

ALASKA: A Mathematica package for two-state kinetic analysis of protein folding reactions

A Mathematica package (ALASKA) has been developed to simplify the measurement of protein folding kinetics by analysis of ¹H NMR lineshape analysis. This package reads NMR data in ASCII format and can simulate an aromatic ¹ NMR spectrum with or without lineshape broadening from chemical exchange. We describe the analysis of a urea denaturation series of a fast-folding protein, the G46A/G48A variant of monomeric repressor.     

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.     

Rapid collapse precedes the fast two-state folding of the cold shock protein

The cold shock protein Bc-Csp folds very rapidly in a reaction that is well described by a kinetic two-state mechanism without intermediates. We measured the shortening of six intra-protein distances during folding by F?rster resonance energy transfer (FRET) in combination with stopped-flow experiments. Single tryptophan residues were engineered into the protein as the donors, and single 5-(((acetylamino)ethyl)amino)naphthalene-1-sulfonate (AEDANS) residues were placed as the acceptors at solvent-exposed sites of Bc-Csp. Their R0 value of about 22 A was well suited for following distance changes during the folding of this protein with a high sensitivity. The mutagenesis and the labeling did not alter the refolding kinetics. The changes in energy transfer during folding were monitored by both donor and acceptor emission and reciprocal effects were found. In two cases the donor-acceptor distances were similar in the unfolded and the folded state and, as a consequence, the kinetic changes in energy transfer upon folding were very small. For four donor/acceptor pairs we found that > or =50% of the increase in energy transfer upon folding occurred prior to the rate-limiting step of folding. This reveals that about half of the shortening of the intra-molecular distances upon folding has occurred already before the rate-limiting step and suggests that the fast two-state folding reaction of Bc-Csp is preceded by a very rapid collapse.     

Fast and slow intermediate accumulation and the initial barrier mechanism in protein folding

Do stable intermediates form very early in the protein folding process? New results and a quantity of literature that bear on this issue are examined here. Results available provide little support for early intermediate accumulation before an initial search-dependent nucleation barrier.     

On the origin of the cooperativity of protein folding: Implications from model simulations

There is considerable experimental evidence that the cooperativity of protein folding resides in the transition from the molten globule to the native state. The objective of this study is to examine whether simplified models can reproduce this cooperativity and if so, to identify its origin. In particular, the thermodynamics of the conformational transition of a previously designed sequence (A. Kolinski, W. Galazka, and J. Skolnick, J. Chem. Phys. 103: 10286–10297, 1995), which adopts a very stable Greek-key β-barrel fold has been investigated using the entropy Monte Carlo sampling (ESMC) technique of Hao and Scheraga (M.-H. Hao and H.A. Scheraga, J. Phys. Chem. 98: 9882–9883, 1994). Here, in addition to the original potential, which includes one body and pair interactions between side chains, the force field has been supplemented by two types of multi-body potentials describing side chain interactions. These potentials facilitate the proteinlike pattern of side chain packing and consequently increase the cooperativity of the folding process. Those models that include an explicit cooperative side chain packing term exhibit a well-defined all-or-none transition from a denatured, random coil state to a high-density, well-defined, nativelike low-energy state. By contrast, models lacking such a term exhibit a conformational transition that is essentially continuous. Finally, an examination of the conformations at the free-energy barrier between the native and denatured states reveals that they contain a substantial amount of native-state secondary structure, about 50% of the native contacts, and have an average root mean square radius of gyration that is about 15% larger than native. © 1996 Wiley-Liss, Inc.     

The folding transition state of the cold shock protein is strongly polarized

It has been shown recently that an 11-residue peptide fragment of transthyretin, TTR(105-115), can form amyloid fibrils in vitro by adopting an extended beta-strand conformation. We used molecular dynamics simulations on systems of TTR(105-115) peptides, for a total length of about 5 micros, to explore the process of self-assembly and the structures of the resulting aggregates. Our results suggest that an antiparallel association of the beta-strands is more probable than a parallel one and that the central residues (T106-L111) in a beta-strand have a high propensity to form inter-peptide hydrogen bonds. The study of the dynamics of self-association indicated that, for this peptide, trajectories leading to conformations with high alpha-helical content are off-pathway from those leading to aggregates with high beta-structure content. We also show that the diverse oligomeric structures that form spontaneously in the molecular dynamics simulations are, to a large extent, compatible with solid-state NMR experimental measurements, including chemical shifts, on fully formed fibrils. The strategy that we present may therefore be used in the design of new experiments to determine the structure of amyloid fibrils, such as those involving site-specific isotope labelling schemes to measure key inter-atomic distances.     

Tryptophan replacements in the trp aporepressor from Escherichia coli: probing the equilibrium and kinetic folding models.

Mutants of the dimeric Escherichia coli trp aporepressor are constructed by replacement of the two tryptophan residues in each subunit in order to assess the effects on equilibrium and kinetic fluorescence properties of the folding reaction. The three kinetic phases detected by intrinsic tryptophan fluorescence in refolding of the wild-type aporepressor are also observed in folding of both Trp 19 to Phe and Trp 99 to Phe single mutants, demonstrating that these phases correspond to global rather than local conformational changes. Comparison of equilibrium fluorescence (Royer, C.A., Mann, C.J., & Matthews, C.R., 1993, Protein Sci. 2, 1844-1852) and circular dichroism transition curves induced by urea shows that replacement of either Trp 19 or Trp 99 results in noncoincident behavior. Unlike the wild-type protein (Gittelman, M.S. & Matthews, C.R., 1990, Biochemistry 29, 7011-7020), tertiary and/or quaternary structures are disrupted at lower denaturant concentration than is secondary structure. The equilibrium results can be interpreted in terms of enhancement in the population of a monomeric folding intermediate in which the lone tryptophan residue is highly exposed to solvent, but in which substantial secondary structure is retained. The location of both mutations at the interface between the two subunits (Zhang, R.G., et al., 1987, Nature 327, 591-597) provides a simple explanation for this phenomenon.     

Protein folding in the secretory pathway of animal cells

The exit of newly-synthesized proteins from the lumen of the endoplasmic reticulum (ER) is the rate-determining step in protein secretion. Only correctly-folded and fully-assembled proteins exit the ER and progress along the secretory pathway. Folding and assembly in the ER are mediated by a variety of factors including folding catalysts and molecular chaperones. The properties of these factors, and the nature of their interactions with folding substrates, are beginning to be clarified. Little work has been done to characterize these processes and these factors in cell lines employed for large-scale cell culture. Manipulation of these process may permit improvement in yield or productivity of recombinant proteins by cultured animal cells.     

Study of a major intermediate in the oxidative folding of leech carboxypeptidase inhibitor: contribution of the fourth disulfide bond

The oxidative folding pathway of leech carboxypeptidase inhibitor (LCI; four disulfide bonds) proceeds through the formation of two major intermediates (III-A and III-B) that contain three native disulfide bonds and act as strong kinetic traps in the folding process. The III-B intermediate lacks the Cys19-Cys43 disulfide bond that links the beta-sheet core with the alpha-helix in wild-type LCI. Here, an analog of this intermediate was constructed by replacing Cys19 and Cys43 with alanine residues. Its oxidative folding follows a rapid sequential flow through one, two, and three disulfide species to reach the native form; the low accumulation of two disulfide intermediates and three disulfide (scrambled) isomers accounts for a highly efficient reaction. The three-dimensional structure of this analog, alone and in complex with carboxypeptidase A (CPA), was determined by X-ray crystallography at 2.2A resolution. Its overall structure is very similar to that of wild-type LCI, although the residues in the region adjacent to the mutation sites show an increased flexibility, which is strongly reduced upon binding to CPA. The structure of the complex also demonstrates that the analog and the wild-type LCI bind to the enzyme in the same manner, as expected by their inhibitory capabilities, which were similar for all enzymes tested. Equilibrium unfolding experiments showed that this mutant is destabilized by approximately 1.5 kcal mol(-1) (40%) relative to the wild-type protein. Together, the data indicate that the fourth disulfide bond provides LCI with both high stability and structural specificity.     

Protein folding

The importance of protein folding in the biosynthesis of proteins is reviewed.     

Oxidative folding of bovine pancreatic ribonuclease A: insight into the overall catalysis of the refolding pathway by phosphate

The effects of the strong stabilizing anion, phosphate, on the oxidative folding of bovine pancreatic ribonuclease A were examined. Phosphate was found to catalyze several steps involved in the oxidative folding process at pH 8.0 and 25°C, resulting in an increase in the rate of pre-equilibration of unstructured species on the folding pathway. In the presence of 400 mM phosphate, the overall increase in the rate of regeneration of native protein was caused primarily by the increased formation and stabilization of tertiary structure in the nativelike intermediates, des-[40-95] and des-[65-72], involved in the rate-determining step. Based on the regeneration of native protein and the stability of Cys Ala substituted mutant analogs of the des-species, (C40A, C95A) and (C65A, C72A), it is suggested that the primary role of phosphate is to catalyze the overall regeneration of native protein through nonspecific electrostatic and hydrogen-bonding effects on the protein and solvent.     

Outlining folding nuclei in globular proteins

Our theoretical approach for prediction of folding/unfolding nuclei in three-dimensional protein structures is based on a search for free energy saddle points on networks of protein unfolding pathways. Under some approximations, this search is performed rapidly by dynamic programming and results in prediction of Phi values, which can be compared with those found experimentally. In this study, we optimize some details of the model (specifically, hydrogen atoms are taken into account in addition to heavy atoms), and compare the theoretically obtained and experimental Phi values (which characterize involvement of residues in folding nuclei) for all 17 proteins, where Phi values are now known for many residues. We show that the model provides good Phi value predictions for proteins whose structures have been determined by X-ray analysis (the average correlation coefficient is 0.65), with a more limited success for proteins whose structures have been determined by NMR techniques only (the average correlation coefficient is 0.34), and that the transition state free energies computed from the same model are in a good anticorrelation with logarithms of experimentally measured folding rates at mid-transition (the correlation coefficient is -0.73).