Experiment and theory highlight role of native state topology in SH3 folding.

We use a combination of experiments, computer simulations and simple model calculations to characterize, first, the folding transition state ensemble of the src SH3 domain, and second, the features of the protein that determine its folding mechanism. Kinetic analysis of mutations at 52 of the 57 residues in the src SH3 domain revealed that the transition state ensemble is even more polarized than suspected earlier: no single alanine substitution in the N-terminal 15 residues or the C-terminal 9 residues has more than a two-fold effect on the folding rate, while such substitutions at 15 sites in the central three-stranded beta-sheet cause significant decreases in the folding rate. Molecular dynamics (MD) unfolding simulations and ab initio folding simulations on the src SH3 domain exhibit a hierarchy of folding similar to that observed in the experiments. The similarity in folding mechanism of different SH3 domains and the similar hierarchy of structure formation observed in the experiments and the simulations can be largely accounted for by a simple native state topology-based model of protein folding energy landscapes.     

Mapping the folding pathway of an immunoglobulin domain: structural detail from Phi value analysis and movement of the transition state

BACKGROUND: Do proteins that have the same structure fold by the same pathway even when they are unrelated in sequence? To address this question, we are comparing the folding of a number of different immunoglobulin-like proteins. Here, we present a detailed protein engineering phi value analysis of the folding pathway of TI I27, an immunoglobulin domain from human cardiac titin. RESULTS: TI I27 folds rapidly via a kinetic intermediate that is destabilized by most mutations. The transition state for folding is remarkably native-like in terms of solvent accessibility. We use phi value analysis to map this transition state and show that it is highly structured; only a few residues close to the N-terminal region of the protein remain completely unfolded. Interestingly, most mutations cause the transition state to become less native-like. This anti-Hammond behavior can be used as a novel means of obtaining additional structural information about the transition state. CONCLUSIONS: The residues that are involved in nucleating the folding of TI I27 are structurally equivalent to the residues that form the folding nucleus in an evolutionary unrelated fibronectin type III protein. These residues form part of the common structural core of Ig-like domains. The data support the hypothesis that interactions essential for defining the structure of these beta sandwich proteins are also important in nucleation of folding.     

COSMIC analysis of the major alpha-helix of barnase during folding

The structures of transition states and intermediates in protein folding may be analysed by protein engineering methods that remove simple interactions that stabilize the folded state. We have now extended the range and reliability of the procedure by using the COSMIC (Combination of Sequential Mutant Interaction Cycles) technique, in which a series of double-mutant cycles is constructed. In each cycle, the side-chains of two amino acid residues that interact in the folded state are mutated separately and together. Kinetic and equilibrium measurements on folding for each cycle show unambiguously whether or not two residues interact during protein folding. A series of such cycles has been constructed to leapfrog along the major alpha-helix of barnase, comprising residues 6 to 18. The helix is found to be intact from its C terminus to residue 12 but begins to unwind towards the N terminus in both the transition state for unfolding and in a folding intermediate.     

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3?kcal?mol^?1. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.     

Autonomous folding of the excised coenzyme-binding domain of D-glyceraldehyde 3-phosphate dehydrogenase from Thermotoga maritima.

An important question in protein folding is whether compact substructures or domains are autonomous units of folding and assembly. The protomer of the tetrameric D-glyceraldehyde-3-phosphate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima has a complex coenzyme-binding domain, in which residues 1-146 form a compact substructure with the last 31 residues (313-333). Here it is shown that the gene of a single-chain protein can be expressed in Escherichia coli after deleting the 163 codons corresponding to the interspersed catalytic domain (150-312). The purified gene product is a soluble, monomeric protein that binds both NAD+ and NADH strongly and possesses the same unfolding transition induced by guanidinium chloride as the native tetramer. The autonomous folding of the coenzyme-binding domain has interesting implications for the folding, assembly, function, and evolution of the native enzyme.     

Generation of a Highly Active Folding Enzyme by Combining a Parvulin-Type Prolyl Isomerase from SurA with an Unrelated Chaperone Domain

Parvulins are small prolyl isomerases and serve as catalytic domains of folding enzymes. SurA (survival protein A) from the periplasm of Escherichia coli consists of an inactive (Par1) and an active (Par2) parvulin domain as well as a chaperone domain. In the absence of the chaperone domain, the folding activity of Par2 is virtually abolished. We created a chimeric protein by inserting the chaperone domain of SlyD, an unrelated folding enzyme from the FKBP family, into a loop of the isolated Par2 domain of SurA. This increased its folding activity 450-fold to a value higher than the activity of SurA, in which Par2 is linked with its natural chaperone domain. In the presence of both the natural and the foreign chaperone domain, the folding activity of Par2 was 1500-fold increased. Related and unrelated chaperone domains thus are similarly efficient in enhancing the folding activity of the prolyl isomerase Par2. A sequence analysis of various chaperone domains suggests that clusters of exposed methionine residues in mobile chain regions might be important for a generic interaction with unfolded protein chains. This binding is highly dynamic to allow frequent transfer of folding protein chains between chaperone and catalytic domains.     

Domain organization and function of Salmonella FliK, a flagellar hook-length control protein

Salmonella hook-length control protein FliK, which consists of 405 amino acid residues, switches substrate specificity of the type III flagellar protein export apparatus from rod/ hook-type to filament-type by causing a conformational change in the cytoplasmic domain of FlhB (FlhB(C)) upon completion of the hook assembly. An N-terminal region of FliK contains an export signal, and a highly conserved C-terminal region consisting of amino acid residues 265-405 (FliK((265-405))) is directly involved in the switching of FlhB(C). Here, we have investigated the structural properties of FliK. Gel filtration chromatography, multi-angle light scattering and analytical ultracentrifugation showed that FliK is monomeric in solution and has an elongated shape. Limited proteolysis showed that FliK consists of two domains, the N-terminal (FliK(N)) and C-terminal domains (FliK(C)), and that the first 203 and the last 35 amino acid residues are partially unfolded and subjected to proteolysis. Both FliK(N) and FliK(C) are more globular than full-length FliK, suggesting that these domains are connected in tandem. Overproduced His-FliK((199-405)) failed to switch export specificity of the export apparatus. Affinity blotting revealed that FlhB(C) binds to FliK and FliK((1-147)), but not to FliK((265-405)). Based on these results, we propose that FliK(N) within the central channel of the hook-basal body during the export of FliK is the sensor and transmitter of hook completion information and that the binding interaction of FliK(C) to FlhB(C) is structurally regulated by FliK(N) so as to occur only when the hook has reached a preset length. The conformational flexibility of FliK(C) may play a role in interfering with switching at an inappropriate point of flagellar assembly.     

The Folding of a Family of Three-Helix Bundle Proteins: Spectrin R15 Has a Robust Folding Nucleus,Unlike Its Homologous Neighbours

Three homologous spectrin domains have remarkably different folding characteristics. We have previously shown that the slow-folding R16 and R17 spectrin domains can be altered to resemble the fast folding R15, in terms of speed of folding (and unfolding), landscape roughness and folding mechanism, simply by substituting five residues in the core. Here we show that, by contrast, R15 cannot be engineered to resemble R16 and R17. It is possible to engineer a slow-folding version of R15, but our analysis shows that this protein neither has a rougher energy landscape nor does change its folding mechanism. Quite remarkably, R15 appears to be a rare example of a protein with a folding nucleus that does not change in position or in size when its folding nucleus is disrupted. Thus, while two members of this protein family are remarkably plastic, the third has apparently a restricted folding landscape.     

Buried hydrophobic side-chains essential for the folding of the parallel beta-helix domains of the P22 tailspike

The processive beta-strands and turns of a polypeptide parallel beta-helix represent one of the topologically simplest beta-sheet folds. The three subunits of the tailspike adhesin of phage P22 each contain 13 rungs of a parallel beta-helix followed by an interdigitated section of triple-stranded beta-helix. Long stacks of hydrophobic residues dominate the elongated buried core of these two beta-helix domains and extend into the core of the contiguous triple beta-prism domain. To test whether these side-chain stacks represent essential residues for driving the chain into the correct fold, each of three stacked phenylalanine residues within the buried core were substituted with less bulky amino acids. The mutant chains with alanine in place of phenylalanine were defective in intracellular folding. The chains accumulated exclusively in the aggregated inclusion body state regardless of temperature of folding. These severe folding defects indicate that the stacked phenylalanine residues are essential for correct parallel beta-helix folding. Replacement of the same phenylalanine residues with valine or leucine also impaired folding in vivo, but with less severity. Mutants were also constructed in a second buried stack that extends into the intertwined triple-stranded beta-helix and contiguous beta-prism regions of the protein. These mutants exhibited severe defects in later stages of chain folding or assembly, accumulating as misfolded but soluble multimeric species. The results indicate that the formation of the buried hydrophobic stacks is critical for the correct folding of the parallel beta-helix, triple-stranded beta-helix, and beta-prism domains in the tailspike protein.     

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.     

Key role of coulombic interactions for the folding transition state of the cold shock protein

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.     

Subdomain interactions as a determinant in the folding and stability of T4 lysozyme.

The folding of large, multidomain proteins involves the hierarchical assembly of individual domains. It remains unclear whether the stability and folding of small, single-domain proteins occurs through a comparable assembly of small, autonomous folding units. We have investigated the relationship between two subdomains of the protein T4 lysozyme. Thermodynamically, T4 lysozyme behaves as a cooperative unit and the unfolding transition fits a two-state model. The structure of the protein, however, resembles a dumbbell with two potential subdomains: an N-terminal subdomain (residues 13-75), and a C-terminal subdomain (residues 76-164 and 1-12). To investigate the effect of uncoupling these two subdomains within the context of the native protein, we created two circular permutations, both at the subdomain interface (residues 13 and 75). Both variants adopt an active wild-type T4 lysozyme fold. The protein starting with residue 13 is 3 kcal/mol less stable than wild type, whereas the protein beginning at residue 75 is 9 kcal/mol less stable, suggesting that the placement of the termini has a major effect on protein stability while minimally affecting the fold. When isolated as protein fragments, the C-terminal subdomain folds into a marginally stable helical structure, whereas the N-terminal subdomain is predominantly unfolded. ANS fluorescence studies indicate that, at low pH, the C-terminal subdomain adopts a loosely packed acid state. An acid state intermediate is also seen for all of the full-length variants. We propose that this acid state is comprised of an unfolded N-terminal subdomain and a loosely folded C-terminal subdomain.     

Molecular thermodynamic model to predict the alpha-helical secondary structure of polypeptide chains in solution.

The native state of a protein molecule in aqueous solutions represents one of the lowest states of Gibbs energy [Anfinsen, C.B. (1973) Science 181, 223-230]. Much progress has been made about the rules of protein folding [King, J. (1989) Chem. Eng. News 67, 32-54] and the dominant forces in protein folding [Dill, K.A. (1990) Biochemistry 29, 7133-7155]. However, the quantitative contributions of different Gibbs energy terms to protein stability remains a controversial issue [Moult, J., & Unger, R. (1991) Biochemistry 30, 3816-3824]. A molecular thermodynamic model has been proposed for the Gibbs energy of folding a residue in aqueous homopolypeptides from a random-coiled state to either the alpha-helix state or the beta-sheet state [Chen, C.-C., Zhu, Y., King, J.A., & Evans, L.B. (1992) Biopolymers 32, 1375-1392]. In this work, we present a generalization of the molecular thermodynamic model for the Gibbs energy of folding natural and synthetic heteropolypeptides from random-coiled conformations into alpha-helical conformations. The generalized model incorporates the intrinsic folding potential due to residue-solvent interactions, the cooperative folding effect due to residue-residue interactions, and the location and length of alpha-helices. The utility of the model was demonstrated by examining the stability of alpha-helical conformations of a number of natural polypeptides including C-peptide (residues 1-13) and S-peptide (residues 1-20) of RNase A (bovine pancreatic ribonuclease A), the P alpha fragment in BPTI (bovine pancreatic trypsin inhibitor), and synthetic polypeptides (the copolymers of different amino acid residues) including alanine-based peptides (16 or 17 residues long) in water. The computed Gibbs energies correspond well with the experimental data on helicity. The results also accounted for the effects of amino acid substitution and temperature on the stability of alpha-helical conformations of the test polypeptides.     

The SH3-fold family: experimental evidence and prediction of variations in the folding pathways

To investigate the relationships between protein topology, amino acid sequence and folding mechanisms, the folding transition state of the Sso7d protein has been characterised both experimentally and theoretically. Although Sso7d protein has a similar topology to that of the SH3 domains, the structure of its transition state is different from that of alpha-spectrin and src SH3 domains previously studied. The folding algorithm, Fold-X, including an energy function with specific sequence features, accounts for these differences and reproduces with a good agreement the set of experimental phi(double dagger-U) values obtained for the three proteins. Our analysis shows that taking into account sequence features underlying protein topology is critical for an accurate prediction of the folding process.     

Plasticity within the obligatory folding nucleus of an immunoglobulin-like domain

A number of β-sandwich immunoglobulin-like domains have been shown to fold using a set of structurally equivalent residues that form a folding nucleus deep within the core of the protein. Formation of this nucleus is sufficient to establish the complex Greek key topology of the native state. These nucleating residues are highly conserved within the immunoglobulin superfamily, but are less well conserved in the fibronectin type III (fnIII) superfamily, where the requirement is simply to have four interacting hydrophobic residues. However, there are rare examples where this nucleation pattern is absent. In this study, we have investigated the folding of a novel member of the fnIII superfamily whose nucleus appears to lack one of the four buried hydrophobic residues. We show that the folding mechanism is unaltered, but the folding nucleus has moved within the hydrophobic core.     

How strong are side chain interactions in the folding intermediate?

Influence of 12 nonpolar amino acids residues from the hydrophobic core of apomyoglobin on stability of its native state and folding intermediate was studied. Six of the selected residues are from the A, G and H helices; these are conserved in structure of the globin family, although nonfunctional, that is, not involved in heme binding. The rest are nonconserved hydrophobic residues that belong to the B, C, D, and E helices. Each residue was substituted by alanine, and equilibrium pH‐induced transitions in apomyoglobin and its mutants were studied by circular dichroism and fluorescent spectroscopy. The obtained results allowed estimating changes in their free energy during formation of the intermediate state. It was first shown that the strength of side chain interactions in the apomyoglobin intermediate state amounts to 15–50% of that in its native state for conserved residues, and practically to 0% for nonconserved residues. These results allow a better understanding of interactions occurring in the intermediate state and shed light on involvement of certain residues in protein folding at different stages.     

pH-dependent stability and folding kinetics of a protein with an unusual alpha-beta topology: the C-terminal domain of the ribosomal protein L9

The folding kinetics and thermodynamics of the isolated C-terminal domain of the ribosomal protein L9 (CTL9) have been studied as a function of pH. CTL9 is an alpha-beta protein that contains a single beta-sheet with an unusual mixed parallel, anti-parallel topology. The folding is fully reversible and two-state over the entire pH range. Stopped-flow fluorescence and CD experiments yield the same folding rate, and the chevron plots have the characteristic V-shape expected for two-state folding. The values of DeltaG*(H2O) and the m value calculated from the kinetic experiments are in excellent agreement with the equilibrium measurements. The extrapolated initial amplitudes of both the stopped-flow fluorescence and CD measurements show that there is no detectable burst phase intermediate. The domain contains three histidine residues, two of which are largely buried in the native state. They do not participate in salt-bridges or take part in a hydrogen bonded network. NMR measurements reveal that the buried histidine residues have significantly perturbed pK(a) values in the native state. The equilibrium stability and the folding rate are found to be strongly dependent upon their ionization state. There is a linear relationship between the log of the folding rate and DeltaG* (H2O) . The protein is much more stable and folds noticeably faster at pH values above the native state pK(a) values. DeltaG*(H2O) of unfolding increases from 2.90 kcal mol(-1) at pH 5.0 to 6.40 kcal mol(-1) at pH 8.0 while the folding rate increases from 0.60 to 18.7 s(-1). Tanford linkage analysis revealed that the interactions involving the two histidine residues are largely developed in the transition state. The results are compared to other studies of the pH-dependence of folding.     

An experimental investigation of conformational fluctuations in proteins G and L

The B1 domains of streptococcal proteins G and L are structurally similar, but they have different sequences and they fold differently. We have measured their NMR spectra at variable temperature using a range of concentrations of denaturant. Many residues have curved amide proton temperature dependence, indicating that they significantly populate alternative, locally unfolded conformations. The results, therefore, provide a view of the locations of low-lying, locally unfolded conformations. They indicate approximately 4-6 local minima for each protein, all within ca. 2.5 kcal/mol of the native state, implying a locally rough energy landscape. Comparison with folding data for these proteins shows that folding involves most molecules traversing a similar path, once a transition state containing a beta hairpin has been formed, thereby defining a well-populated pathway down the folding funnel. The hairpin that directs the folding pathway differs for the two proteins and remains the most stable part of the folded protein.     

The folding nucleus of a fibronectin type III domain is composed of core residues of the immunoglobulin-like fold

To identify the contacts that stabilise the rate-limiting transition state for folding of FNfn10 (the tenth fnIII domain of human fibronectin), 42 mutants have been analysed at 29 positions across this domain. An anomalous response to mutation means that structure formation in the A, B and G strands cannot be evaluated by this method. In all the residues analysed, phi-values are fractional and no completely structured region is observed. The analysis reveals that hydrophobic residues from the central strands of the beta-sandwich form a large core of interactions in the transition state. Br?nsted analysis shows that the stabilisation energy from the amino acid side-chains in the transition state is approximately 40 % of that in the native state. The protein folds by a nucleation-condensation mechanism, and tertiary interactions within the core make up the folding nucleus. Local interactions, in turns and loops, are apparently much less significant. Comparison with an homologous domain from human tenascin (TNfn3), shows that FNfn10 has a more extended, structured transition state spanning three different "layers" of the beta-sandwich. The results support the hypothesis that interactions in the common structural core guide the folding of these domains.