Sensitivity of the Folding Pathway to the Details of Amino Acid Sequence

Experimental data on the structure of the transition state demonstrate that proteins with the same topology as a rule have similar folding nuclei (the structured formed part of the transition state). In this review discussed are the experimental works showing that the position of folding nuclei is different among proteins with the same topology. These facts emphasize that the folding pathway is sensitive to the details of amino acid sequence.     

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.     

Role of electrostatic interactions for the stability and folding behavior of cold shock protein

The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

GroEL Can Unfold Late Intermediates Populated on the Folding Pathways of Monellin

The modulation of the folding mechanism of the small protein single-chain monellin (MNEI) by the Escherichia coli chaperone GroEL has been studied. In the absence of the chaperone, the folding of monellin occurs via three parallel routes. When folding is initiated in the presence of a saturating concentration of GroEL, only 50-60% of monellin molecules fold completely. The remaining 40-50% of the monellin molecules remain bound to the GroEL and are released only upon addition of ATP. It is shown that the basic folding mechanism of monellin is not altered by the presence of GroEL, but that it occurs via only one of the three available routes when folding is initiated in the presence of saturating concentrations of GroEL. Two pathways become nonoperational because GroEL binds very tightly to early intermediates that populate these pathways in a manner that makes the GroEL-bound intermediates incompetent to fold. This accounts for the monellin molecules that remain GroEL-bound at the end of the folding reaction. The third pathway remains operational because the GroEL-bound early intermediate on this pathway is folding-competent, suggesting that this early intermediate binds to GroEL in a manner that is different from that of the binding of the early intermediates on the other two pathways. It appears, therefore, that the same protein can bind GroEL in more than one way. The modulation of the folding energy landscape of monellin by GroEL occurs because GroEL binds folding intermediates on parallel folding pathways, in different ways, and with different affinities. Moreover, when GroEL is added to refolding monellin at different times after commencement of refolding, the unfolding of two late kinetic intermediates on two of the three folding pathways can be observed. It appears that the unfolding of late folding intermediates is enabled by a thermodynamic coupling mechanism, wherein GroEL binds more tightly to an early intermediate than to a late intermediate on a folding pathway, with preferential binding energy being larger than the stability of the late intermediate. Hence, it is shown that GroEL can inadvertently and passively cause, through its ability to bind different folding intermediates differentially, the unfolding of late productive intermediates on folding pathways, and that its unfolding action is not restricted solely to misfolded or kinetically trapped intermediates.     

Pathway of Oxidative Folding of a 3-Disulfide α-Lactalbumin May Resemble Either BPTI Model or Hirudin Model

Pathways of oxidative folding of disulfide proteins display a high degree of diversity and vary among two extreme models. The BPTI model is defined by limited species of folding intermediates adopting mainly native disulfide bonds. The hirudin model is characterized by highly heterogeneous folding intermediates containing mostly non-native disulfide bonds. αLA-IIIA is a 3-disulfide variant of α-lactalbumin (αLA) with a 3-D conformation essentially identical to that of intact αLA. αLA-IIIA contains 3 native disulfide bonds of αLA, two of them are located at the calcium binding β-subdomain (Cys⁶¹–Cys⁷⁷ and Cys⁷³–Cys⁹¹) and the third bridge is located within the α-helical domain of the molecule (Cys²⁸–Cys¹¹¹). We investigate here the pathway of oxidative folding of fully reduced αLA-IIIA with and without stabilization of its β-subdomain by calcium binding. In the absence of calcium, the folding pathway of αLA-IIIA was shown to resemble that of hirudin model. Upon stabilization of β-sheet domain by calcium binding, the folding pathway of αLA-IIIA exhibits a striking similarity to that of BPTI model. Three predominant folding intermediates of αLA-IIIA containing exclusively native disulfide bonds were isolated and structurally characterized. Our results further demonstrate that stabilization of subdomains in a protein may dictate its folding pathway and represent a major cause for the existing diversity in the folding pathways of the disulfide-containing proteins.     

The nature of protein folding pathways: The classical versus the new view

Summary Pulsed hydrogen exchange and other studies of the kinetic refolding pathways of several small proteins have established that folding intermediates with native-like secondary structures are well populated, but these studies have also shown that the folding kinetics are not well synchronized. Older studies of the kinetics of formation of the native protein, monitored by optical probes, indicate that the folding kinetics should be synchronized. The model commonly used in these studies is the simple sequential model, which postulates a unique folding pathway with defined and sequential intermediates. Theories of the folding process and Monte Carlo simulations of folding suggest that neither the folding pathway nor the set of folding intermediates is unique, and that folding intermediates accumulate because of kinetic traps caused by partial misfolding. Recent experiments with cytochrome c lend support to this new view of folding pathways. These different views of the folding process are discussed. Misfolding and consequent slowing down of the folding process as a result of cis-trans isomerization about prolyl peptide bonds in the unfolded protein are well known; isomerization occurs before refolding is initiated. The occurrence of equilibrium intermediates on the kinetic folding pathways of some proteins, such as -lactalbumin and apomyoglobin, argues that these intermediates are not caused by kinetic traps but rather are stable intermediates under certain conditions, and this conclusion is consistent with a sequential model of folding. Folding reactions with successive kinetic intermediates, in which late intermediates are more highly folded than early intermediates, indicate that folding is hierarchical. New experiments that test the predictions of the classical and the new views are needed.     

Mutational analysis of the BPTI folding pathway: II. Effects of aromatic-->leucine substitutions on folding kinetics and thermodynamics.

The rates of the individual steps in the disulfide-coupled folding and unfolding of eight BPTI variants, each containing a single aromatic to leucine amino acid replacement, were measured. From this analysis, the contributions of the four phenylalanine and four tyrosine residues to the stabilities of the native protein and the disulfide-bonded folding intermediates were determined. While the substitutions were found to destabilize the native protein by 2 to 7 kcal/mol, they had significantly smaller effects on the intermediates that represent the earlier stages of folding, even when the site of the substitution was located within the ordered regions of the intermediates. These results suggest that stabilizing interactions contribute less to conformational stability in the context of a partially folded intermediate than in a fully folded native protein, perhaps because of decreased cooperativity among the individual interactions. The kinetic analysis also provides new information about the transition states associated with the slowest steps in folding and unfolding, supporting previous suggestions that these transition states are extensively unfolded. Although the substitutions caused large changes in the distribution of folding intermediates and in the rates of some steps in the folding pathway, the kinetically-preferred pathway for all of the variants involved intramolecular disulfide rearrangements, as observed previously for the wild-type protein. These results suggest that the predominance of the rearrangement mechanism reflects conformational constraints present relatively early in the folding pathway.     

Models for protein folding and nature's choice of protein as catalyst

The study of protein folding and unfolding pathways lends a fascinating dimension to protein biochemistry. Several models for protein folding have been postulated. Two powerful probes used in protein folding study are far UV-CD monitored stopped flow kinetics and pulse hydrogen exchange in conjunction with NMR. The formation of molten globule, which is an intermediate possessing secondary structure but not a well packed tertiary structure, is now emerging as a common feature on the folding pathway of many proteins. The molten globule is recognized by a class of molecules called chaperones which act as accelerators of protein folding. This article ends by elucidating why proteins are Nature's choice as catalysts.     

Substitutions of prolines examine their role in kinetic trap formation of the caspase recruitment domain (CARD) of RICK

Caspase recruitment domains (CARDs) are small helical protein domains that adopt the Greek key fold. For the two CARDs studied to date, RICK-CARD and caspase-1-CARD (CP1-CARD), the proteins unfold by an apparent two-state process at equilibrium. However, the folding kinetics are complex for both proteins and may contain kinetically trapped species on the folding pathway. In the case of RICK-CARD, the time constants of the slow refolding phases are consistent with proline isomerism. RICK-CARD contains three prolines, P47 in turn 3, and P85 and P87. The latter two prolines constitute a nonconserved PxP motif in helix 6. To examine the role of the prolines in the complex folding kinetics of RICK-CARD, we generated seven proline-to-alanine mutants, including three single mutants, three double mutants, and one triple mutant. We examined the spectroscopic properties, equilibrium folding, binding to CP1-CARD, and folding kinetics. The results show that P85 is critical for maintaining the function of the protein and that all mutations decrease the stability. Results from single mixing and sequential mixing stopped-flow studies strongly suggest the presence of parallel folding pathways consisting of at least two unfolded populations. The mutations affect the distribution of the two unfolded species, thereby affecting the population that folds through each channel. The two conformations also are present in the triple mutant, demonstrating that interconversion between them is not due to prolyl isomerism. Overall, the data show that the complex folding pathway, especially formation of kinetically trapped species, is not due to prolyl isomerism.     

Incorporating into a Cα Go model the effects of geometrical restriction on Cα atoms caused by side chain orientations

Coarse‐grained Go models have been widely used for studying protein‐folding mechanisms. Despite the simplicity of the model, these can reproduce the essential features of the folding process of a protein. However, it is also known that side chains significantly contribute to the folding mechanism. Hence, it is desirable to incorporate the side chain effects into a coarse‐grained Go model. In this study, to distinguish the effects of side chain orientation and to understand how these effects contribute to folding mechanisms, we incorporate into a Cα Go model not only heterogeneous contact energies but also geometrical restraints around two Cα atoms in contact with each other. We confirm that the heterogeneity of contact energies governs the folding pathway of a protein and that the geometric constraints attributed to side chains reproduce cooperative transitions in folding. Proteins 2013; 81:1434–1445. © 2013 Wiley Periodicals, Inc.     

Simulation and experiment at high temperatures: ultrafast folding of a thermophilic protein by nucleation-condensation

We used Phi-value analysis to characterise the transition state for folding of a thermophilic protein at the relatively high temperature of 325 K. PhiF values for the folding of the three-helix bundle, peripheral subunit binding domain from Bacillus stearothermophilus (E3BD) were determined by temperature-jump experiments in the absence of chemical denaturants. E3BD folded in microseconds through a highly diffuse transition state. Excellent agreement was observed between experiment and the results from eight (independent) molecular dynamics simulations of unfolding at 373 K. We used a combination of heteronuclear NMR experiments and molecular dynamics simulations to characterise the denatured ensemble, and found that it contained very little persistent, residual structure. However, those regions that adopt helical structure in the native state were found by simulation to be poised for helix formation in the denatured state. These regions also had significant structure in the transition state for folding. The overall folding pathway appears to be nucleation-condensation.     

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 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.     

In vitro insulin refolding: Characterization of the intermediates and the putative folding pathway

The in vitro refolding process of the double-chain insulin was studied based on the investigation of in vitro single-chain insulin refolding. Six major folding intermediates, named P1A, P2B, P3A, P4B, P5B, and P6B, were captured during the folding process. The refolding experiments indicate that all of these intermediates are on-pathway. Based on these intermediates and the formation of hypothetic transients, we propose a two-stage folding pathway of insulin. (1) At the early stage of the folding process, the reduced A chain and B chain individually formed the intermediates two A chain intermediates (P1A and P3A), and four B chain intermediates (P2B, P4B, P5B, and P6B). (2) In the subsequent folding process, transient Ⅰ was formed from P3A through thiol/disulfide exchange reaction; then, transients Ⅱ and Ⅲ, each containing two native disulfides, were formed through the recognition and interaction of transient Ⅰ with P4B or P6B and the thiol group's oxidation reaction mainly using GSSG as oxidative reagent; finally, transients Ⅱ and Ⅲ, through thiol/mixture disulfide exchange reaction, formed the third native disulfide of insulin to complete the folding.     

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.     

Differentiation between two-state and multi-state folding proteins based on sequence

Prediction of protein-folding rates follows different rules in two-state and multi-state kinetics. The prerequisite for the prediction is to recognize the folding kinetic pathway of proteins. Here, we use the logistic regression and support vector machine to discriminate between two-state and multi-state folding proteins. We find that chain length is sufficient to accurately recognize multi-state proteins. There is a transition boundary between two kinetic models. Protein folds with multi-state kinetics, if its length is larger than 112 residues. The logistic prediction from amino acid composition shows that the kinetic pathway of folding is closely related to amino acid volume. Small amino acids make two-state folding easier, and vice versa. However, cysteine, alanine, arginine, lysine, histidine, and methionine do not conform to this rule.     

Divalent metal cofactor binding in the kinetic folding trajectory of Escherichia coli ribonuclease HI

Proteins often require cofactors to perform their biological functions and must fold in the presence of their cognate ligands. Using circular dichroism spectroscopy. we investigated the effects of divalent metal binding upon the folding pathway of Escherichia coli RNase HI. This enzyme binds divalent metal in its active site, which is proximal to the folding core of RNase HI as defined by hydrogen/deuterium exchange studies. Metal binding increases the apparent stability of native RNase HI chiefly by reducing the unfolding rate. As with the apo-form of the protein, refolding from high denaturant concentrations in the presence of Mg2+ follows three-state kinetics: formation of a rapid burst phase followed by measurable single exponential kinetics. Therefore, the overall folding pathway of RNase HI is minimally perturbed by the presence of metal ions. Our results indicate that the metal cofactor enters the active site pocket only after the enzyme reaches its native fold, and therefore, divalent metal binding stabilizes the protein by decreasing its unfolding rate. Furthermore, the binding of the cofactor is dependent upon a carboxylate critical for activity (Asp10). A mutation in this residue (D10A) alters the folding kinetics in the absence of metal ions such that they are similar to those observed for the unaltered enzyme in the presence of metal.     

<Emphasis Type="Italic">In vitro</Emphasis> insulin refolding: Characterization of the intermediates and the putative folding pathway

The in vitro refolding process of the double-chain insulin was studied based on the investigation of in vitro single-chain insulin refolding. Six major folding intermediates, named P1A, P2B, P3A, P4B, P5B, and P6B, were captured during the folding process. The refolding experiments indicate that all of these intermediates are on-pathway. Based on these intermediates and the formation of hypothetic transients, we propose a two-stage folding pathway of insulin. (1) At the early stage of the folding process, the reduced A chain and B chain individually formed the intermediates: two A chain intermediates (P1A and P3A), and four B chain intermediates (P2B, P4B, P5B, and P6B). (2) In the subsequent folding process, transient I was formed from P3A through thiol/disulfide exchange reaction; then, transients II and III, each containing two native disulfides, were formed through the recognition and interaction of transient I with P4B or P6B and the thiol group’s oxidation reaction mainly using GSSG as oxidative reagent; finally, transients II and III, through thiol/mixture disulfide exchange reaction, formed the third native disulfide of insulin to complete the folding.     

Asymmetric effect of domain interactions on the kinetics of folding in yeast phosphoglycerate kinase

The aim of this work is to shed more light on the effect of domain-domain interactions on the kinetics and the pathway of protein folding. A model protein system consisting of several single-tryptophan variants of the two-domain yeast phosphoglycerate kinase (PGK) and its individual domains was studied. Refolding was initiated from the guanidine-unfolded state by stopped-flow or manual mixing and monitored by tryptophan fluorescence from 1 msec to 1000 sec. Denaturant titrations of both individual domains showed apparent two-state unfolding transitions. Refolding kinetics of the individual domains from different denaturant concentrations, however, revealed the presence of intermediate structures during titration for both domains. Refolding of the same domains within the complete protein showed that domain-domain interactions direct the folding of both domains, but in an asymmetric way. Folding of the N domain was already altered within 1 msec, while detectable changes in the folding of the C domain occurred only 60-100 msec after initiating refolding. All mutants showed a hyperfluorescent kinetic intermediate. Both the disappearance of this intermediate and the completion of the folding were significantly faster in the individual N domain than in the complete protein. On the contrary, folding of the individual C domain was slower than in the complete protein. The presence of the C domain directs the refolding of the N domain along a completely different pathway than that of the individual N domain, while folding of the individual C domain follows the same path as within the complete protein.