A breakdown of symmetry in the folding transition state of protein L

The 62 residue IgG binding domain of protein L consists of a central alpha-helix packed on a four-stranded beta-sheet formed by N and C-terminal beta-hairpins. The overall topology of the protein is quite symmetric: the beta-hairpins have similar lengths and make very similar interactions with the central helix. Characterization of the effects of 70 point mutations distributed throughout the protein on the kinetics of folding and unfolding reveals that this symmetry is completely broken during folding; the first beta-hairpin is largely structured while the second beta-hairpin and helix are largely disrupted in the folding transition state ensemble. The results are not consistent with a "hydrophobic core first" picture of protein folding; the first beta-hairpin appears to be at least as ordered at the rate limiting step in folding as the hydrophobic core.     

Diffusion-collision of foldons elucidates the kinetic effects of point mutations and suggests control strategies of the folding process of helical proteins

In this article we use mutation studies as a benchmark for a minimal model of the folding process of helical proteins. The model ascribes a pivotal role to the collisional dynamics of a few crucial residues (foldons) and predicts the folding rates by exploiting information drawn from the protein sequence. We show that our model rationalizes the effects of point mutations on the kinetics of folding. The folding times of two proteins and their mutants are predicted. Stability and location of foldons have a critical role as the determinants of protein folding. This allows us to elucidate two main mechanisms for the kinetic effects of mutations. First, it turns out that the mutations eliciting the most notable effects alter protein stability through stabilization or destabilization of the foldons. Secondly, the folding rate is affected via a modification of the foldon topology by those mutations that lead to the birth or death of foldons. The few mispredicted folding rates of some mutants hint at the limits of the current version of the folding model proposed in the present article. The performance of our folding model declines in case the mutated residues are subject to strong long-range forces. That foldons are the critical targets of mutation studies has notable implications for design strategies and is of particular interest to address the issue of the kinetic regulation of single proteins in the general context of the overall dynamics of the interactome.     

Understanding the mechanism of beta-hairpin folding via phi-value analysis

The folding kinetics of a 16-residue beta-hairpin (trpzip4) and five mutants were studied by a laser-induced temperature-jump infrared method. Our results indicate that mutations which affect the strength of the hydrophobic cluster lead to a decrease in the thermal stability of the beta-hairpin, as a result of increased unfolding rates. For example, the W45Y mutant has a phi-value of approximately zero, implying a folding transition state in which the native contacts involving Trp45 are not yet formed. On the other hand, mutations in the turn or loop region mostly affect the folding rate. In particular, replacing Asp46 with Ala leads to a decrease in the folding rate by roughly 9 times. Accordingly, the phi-value for D46A is determined to be approximately 0.77, suggesting that this residue plays a key role in stabilizing the folding transition state. This is most likely due to the fact that the main chain and side chain of Asp46 form a characteristic hydrogen bond network with other residues in the turn region. Taken together, these results support the folding mechanism we proposed before, which suggests that the turn formation is the rate-limiting step in beta-hairpin folding and, consequently, a stronger turn-promoting sequence increases the stability of a beta-hairpin primarily by increasing its folding rate, whereas a stronger hydrophobic cluster increases the stability of a beta-hairpin primarily by decreasing its unfolding rate. In addition, we have examined the compactness of the thermally denatured and urea-denatured states of another 16-residue beta-hairpin, using the method of fluorescence resonance energy transfer. Our results show that the thermally denatured state of this beta-hairpin is significantly more compact than the urea-denatured state, suggesting that the very first step in beta-hairpin folding, when initiated from an extended conformation, probably corresponds to a process of hydrophobic collapse.     

Thermodynamics of a beta-hairpin structure: evidence for cooperative formation of folding nucleus

To elucidate early nucleation stages in protein folding, multi-probed thermodynamic characterization was applied to the beta-hairpin structural formation of G-peptide, which is a C-terminal fragment of the B1 domain of streptococcal protein G. The segment corresponding to the sequence of G-peptide is believed to act as a nucleus during the folding process of the B1 domain. In spite of the broad thermal transition of G-peptide, nuclear magnetic resonance (NMR) melting measurements combined with our original analytical theory enabled us to obtain the thermodynamic properties of the beta-hairpin formation with considerable accuracy. Additionally, all the thermodynamic properties determined by every NMR probe on both the main-chain and the side-chains were quite similar, and also comparable to the values that were independently determined by calorimetric analysis of G-peptide. These results demonstrate that G-peptide folds cooperatively throughout the molecule. In other words, the formation of the beta-hairpin is interpreted as the fashion of a first-order phase transition between two states without any distinguishable intermediates. This cooperative formation of the short linear peptide consisting of only 16 residues provides insight into not only the first folding events of the B1 domain, but also the general principles of proteins in terms of structural hierarchy, stability and folding mechanism.     

Deciphering the Hidden Informational Content of Protein Sequences: FOLDABILITY OF PROINSULIN HINGES ON A FLEXIBLE ARM THAT IS DISPENSABLE IN THE MATURE HORMONE*

Protein sequences encode both structure and foldability. Whereas the interrelationship of sequence and structure has been extensively investigated, the origins of folding efficiency are enigmatic. We demonstrate that the folding of proinsulin requires a flexible N-terminal hydrophobic residue that is dispensable for the structure, activity, and stability of the mature hormone. This residue (Phe^B1 in placental mammals) is variably positioned within crystal structures and exhibits ¹H NMR motional narrowing in solution. Despite such flexibility, its deletion impaired insulin chain combination and led in cell culture to formation of non-native disulfide isomers with impaired secretion of the variant proinsulin. Cellular folding and secretion were maintained by hydrophobic substitutions at B1 but markedly perturbed by polar or charged side chains. We propose that, during folding, a hydrophobic side chain at B1 anchors transient long-range interactions by a flexible N-terminal arm (residues B1–B8) to mediate kinetic or thermodynamic partitioning among disulfide intermediates. Evidence for the overall contribution of the arm to folding was obtained by alanine scanning mutagenesis. Together, our findings demonstrate that efficient folding of proinsulin requires N-terminal sequences that are dispensable in the native state. Such arm-dependent folding can be abrogated by mutations associated with β-cell dysfunction and neonatal diabetes mellitus.     

Beta-hairpin folding simulations in atomistic detail using an implicit solvent model

We have used distributed computing techniques and a supercluster of thousands of computer processors to study folding of the C-terminal beta-hairpin from protein G in atomistic detail using the GB/SA implicit solvent model at 300 K. We have simulated a total of nearly 38 micros of folding time and obtained eight complete and independent folding trajectories. Starting from an extended state, we observe relaxation to an unfolded state characterized by non-specific, temporary hydrogen bonding. This is followed by the appearance of interactions between hydrophobic residues that stabilize a bent intermediate. Final formation of the complete hydrophobic core occurs cooperatively at the same time that the final hydrogen bonding pattern appears. The folded hairpin structures we observe all contain a closely packed hydrophobic core and proper beta-sheet backbone dihedral angles, but they differ in backbone hydrogen bonding pattern. We show that this is consistent with the existing experimental data on the hairpin alone in solution. Our analysis also reveals short-lived semi-helical intermediates which define a thermodynamic trap. Our results are consistent with a three-state mechanism with a single rate-limiting step in which a varying final hydrogen bond pattern is apparent, and semi-helical off-pathway intermediates may appear early in the folding process. We include details of the ensemble dynamics methodology and a discussion of our achievements using this new computational device for studying dynamics at the atomic level.     

Role of side-chains in the cooperative beta-hairpin folding of the short C-terminal fragment derived from streptococcal protein G

A short C-terminal fragment of immunoglobulin-binding domain of streptococcal protein G is known to form nativelike beta-hairpin at physiological conditions. To understand the cooperative folding of the short peptide, eight Ala-substituted mutants of the fragment were investigated with respect to their structural stabilities by analyzing temperature dependence of NMR signals. On comparison of the obtained thermodynamic parameters, we found that the nonpolar residues Tyr45 and Phe52 and the polar residues Asp46 and Thr49 are crucial for the beta-hairpin folding. The results suggest a strong interaction between the nonpolar side chains that participates in a putative hydrophobic cluster and that the polar side chains form a fairly rigid conformation around the loop (46-51). We also investigated the complex formation of the mutants with N-terminal fragment at the variety of temperature to get their thermal unfolding profiles and found that the mutations on the residues Asp46 and Thr49 largely destabilized the complexes, while substitution of Asp47 slightly stabilized the complex. From these results, we deduced that both the hydrophobic cluster formation and the rigidity of the loop (46-51) cooperatively stabilize the beta-hairpin structure of the fragment. These interactions which form a stable beta-hairpin may be the initial structural scaffold which is important in the early folding events of the whole domain.     

The mechanism of beta-hairpin formation

Beta-hairpins constitute an important class of connecting protein secondary structures. Several groups have postulated that such structures form early in the folding process and serve to nucleate the formation of extended beta-sheet structures. Despite the importance of beta-hairpins in protein folding, little is known about the mechanism of formation of these structures. While it is well established that there is a complex interplay between the stability of a beta-hairpin and loop conformational propensity, loop length, and the formation of stabilizing cross-strand interactions (H-bonds and hydrophobic interactions), the influence of these factors on the folding rate is poorly understood. Peptide models provide a simple framework for exploring the molecular details of the formation of beta-hairpin structures. We have explored the fundamental processes of folding in two linear peptides that form beta-hairpin structures, having a stabilizing hydrophobic cluster connected by loops of differing lengths. This approach allows us to evaluate existing models of the mechanism of beta-hairpin formation. We find a substantial acceleration of the folding rate when the connecting loop is made shorter (i.e., the hydrophobic cluster is moved closer to the turn). Analysis of the folding kinetics of these two peptides reveals that this acceleration is a direct consequence of the reduced entropic cost of the smaller loop search.     

Are residues in a protein folding nucleus evolutionarily conserved?

Protein is the working molecule of the cell, and evolution is the hallmark of life. It is important to understand how protein folding and evolution influence each other. Several studies correlating experimental measurement of residue participation in folding nucleus and sequence conservation have reached different conclusions. These studies are based on assessment of sequence conservation at folding nucleus sites using entropy or relative entropy measurement derived from multiple sequence alignment. Here we report analysis of conservation of folding nucleus using an evolutionary model alternative to entropy-based approaches. We employ a continuous time Markov model of codon substitution to distinguish mutation fixed by evolution and mutation fixed by chance. This model takes into account bias in codon frequency, bias-favoring transition over transversion, as well as explicit phylogenetic information. We measure selection pressure using the ratio omega of synonymous versus non-synonymous substitution at individual residue site. The omega-values are estimated using the PAML method, a maximum-likelihood estimator. Our results show that there is little correlation between the extent of kinetic participation in protein folding nucleus as measured by experimental phi-value and selection pressure as measured by omega-value. In addition, two randomization tests failed to show that folding nucleus residues are significantly more conserved than the whole protein, or the median omega value of all residues in the protein. These results suggest that at the level of codon substitution, there is no indication that folding nucleus residues are significantly more conserved than other residues. We further reconstruct candidate ancestral residues of the folding nucleus and suggest possible test tube mutation studies for testing folding behavior of ancient folding nucleus.     

Neutral evolution of model proteins: diffusion in sequence space and overdispersion.

We stimulate the evolution of model protein sequences subject to mutations. A mutation is considered neutral if it conserves (1) the structure of the ground state, (2) its thermodynamic stability and (3) its kinetic accessibility. All other mutations are considered lethal and are rejected. We adopt a lattice model, amenable to a reliable solution of the protein folding problem. We prove the existence of extended neutral networks in sequence space-sequences can evolve until their similarity with the starting point is almost the same as for random sequences. Furthermore, we find that the rate of neutral mutations has a broad distribution in sequence space. Due to this fact, the substitution process is overdispersed (the ratio between variance and mean is larger than 1). This result is in contrast with the simplest model of neutral evolution, which assumes a Poisson process for substitutions, and in qualitative agreement with the biological data.     

Understanding the folding and stability of a zinc finger-based full sequence design protein with replica exchange molecular dynamics simulations

Full sequence design protein FSD-1 is a designed protein based on the motif of zinc finger protein. In this work, its folding mechanism and thermal stability are investigated using the replica exchange molecular dynamics model with the water molecules being treated explicitly. The results show that the folding of the FSD-1 is initiated by the hydrophobic collapse, which is accompanied with the formation of the C-terminal alpha-helix. Then the folding proceeds with the formation of the beta-hairpin and the further package of the hydrophobic core. Compared with the beta-hairpin, the alpha-helix has much higher stability. It is also found that the N-capping motif adopted by the FSD-1 contributes to the stability of the alpha-helix dramatically. The hydrophobic contacts made by the side chain of Tyr3 in the native state are essential for the stabilization of the beta-hairpin. It is also found that the folding of the N-terminal beta-hairpin and the C-terminal alpha-helix exhibits weak cooperativity, which is consistent with the experimental data. Meanwhile, the folding pathway is compared between the FSD-1 and the target zinc finger peptide, and the possible role of the zinc ion on the folding pathway of zinc finger is proposed. Proteins 2007. (c) 2007 Wiley-Liss, Inc.     

Bergerac-SH3: "frustation" induced by stabilizing the folding nucleus

The influence of an inserted exogenous independent folding element on the thermodynamics and folding properties of SH3 domain from alpha-spectrin has been investigated by creating a fused form between this small all-beta domain and a stable beta-hairpin (BH19). NMR analysis of synthetic peptides shows that insertion of BH19 nucleates formation of the original natural beta-hairpin (distal loop) that is part of the SH3 folding nucleus. The resulting protein (Bergerac-SHH) is more stable, folds faster and contains an elongated hairpin protruding from the globular domain as determined by 2D-NMR. "Protein engineering" analysis of the inserted region shows that it is folded in the transition state. Interestingly, stabilisation by insertion of the distal loop region results in the appearance of a compact intermediate revealed by a curved chevron plot at low denaturant concentration. This effect is eliminated at low salt concentrations by a single mutation of a hydrophobic residue within BH19 sequence, which is most probably involved in non-native interactions. Local stabilisation by enlargement and reinforcement of the folding nucleus, global compaction by the addition of salt and non-native interactions are shown to contribute to the observed deviation from the two-state behaviour.     

Weak cooperativity in the core causes a switch in folding mechanism between two proteins of the cks family

The human protein ckshs1 (cks1) is a 79 residue alpha/beta protein with low thermodynamic and kinetic stability. Its folding mechanism was probed by mutation at sites throughout the structure. Many of the mutations caused changes in the slope of the unfolding arm of the chevron plot. The effects can be rationalised in terms of either transition-state movement or native-state "breathing", and in either case, the magnitude of the effect enables the sequence of events in the folding reaction to be determined. Those sites that fold early exhibit a small perturbation, whilst those sites that fold late exhibit a large perturbation. The results show that cks1 folds sequential pairs of beta-strands first; beta1/beta2 and beta3/beta4. Subsequently, these pairs pack against each other and onto the alpha-helical region to form the core. The folding process of cks1 contrasts with that of the homologue, suc1. The 113 residue suc1 has the same beta-sheet core structure but, additionally, two large insertions that confer much greater thermodynamic and kinetic stability. The more extensive network of tertiary interactions in suc1 provides sufficient enthalpic gain to overcome the entropic cost of forming the core and thus tips the balance in favour of non-local interactions: the non-local, central beta-strand pair, beta2/beta4, forms first and the periphery strands pack on later. Moreover, the greater cooperativity of the core of suc1 protects its folding from perturbation and consequently the slope of the unfolding arm of the chevron plot is much less sensitive to mutation.     

Identification of a key structural element for protein folding within beta-hairpin turns

Specific residues in a polypeptide may be key contributors to the stability and foldability of the unique native structure. Identification and prediction of such residues is, therefore, an important area of investigation in solving the protein folding problem. Atypical main-chain conformations can help identify strains within a folded protein, and by inference, positions where unique amino acids may have a naturally high frequency of occurrence due to favorable contributions to stability and folding. Non-Gly residues located near the left-handed alpha-helical region (L-alpha) of the Ramachandran plot are a potential indicator of structural strain. Although many investigators have studied mutations at such positions, no consistent energetic or kinetic contributions to stability or folding have been elucidated. Here we report a study of the effects of Gly, Ala and Asn substitutions found within the L-alpha region at a characteristic position in defined beta-hairpin turns within human acidic fibroblast growth factor, and demonstrate consistent effects upon stability and folding kinetics. The thermodynamic and kinetic data are compared to available data for similar mutations in other proteins, with excellent agreement. The results have identified that Gly at the i+3 position within a subset of beta-hairpin turns is a key contributor towards increasing the rate of folding to the native state of the polypeptide while leaving the rate of unfolding largely unchanged.     

Interplay between hydrophobic cluster and loop propensity in beta-hairpin formation

Autonomously folding beta-hairpins have recently emerged as powerful tools for elucidating the origins of antiparallel beta-sheet folding preferences. Analysis of such model systems has suggested four potential sources of beta-sheet stability: (1) the conformational propensity of the loop segment that connects adjacent strands; (2) favorable contacts between side-chains on adjacent strands; (3) interstrand hydrogen bonds; and (4) the intrinsic beta-sheet propensities of the strand residues. We describe the design and analysis of a series of isomeric 20 residue peptides in which factors (1)-(4) are identical. Differences in beta-hairpin formation within this series demonstrate that these four factors, individually, are not sufficient to explain beta-sheet stability. In agreement with the prediction of a simple statistical mechanical model for beta-hairpin formation, our results show that the separation between the loop segment and an interstrand cluster of hydrophobic side-chains strongly influences beta-hairpin size and stability, with a smaller separation leading to greater stability.     

Mutational analysis of the folding transition state of the C-terminal domain of ribosomal protein L9: a protein with an unusual beta-sheet topology

The C-terminal domain of ribosomal protein L9 (CTL9) is a 92-residue alpha-beta protein which contains an unusual three-stranded mixed parallel and antiparallel beta-sheet. The protein folds in a two-state fashion, and the folding rate is slow. It is thought that the slow folding may be caused by the necessity of forming this unusual beta-sheet architecture in the transition state for folding. This hypothesis makes CTL9 an interesting target for folding studies. The transition state for the folding of CTL9 was characterized by phi-value analysis. The folding of a set of hydrophobic core mutants was analyzed together with a set of truncation mutants. The results revealed a few positions with high phi-values (> or = 0.5), notably, V131, L133, H134, V137, and L141. All of these residues were found in the beta-hairpin region, indicating that the formation of this structure is likely to be the rate-limiting step in the folding of CTL9. One face of the beta-hairpin docks against the N-terminal helix. Analysis of truncation mutants of this helix confirmed its importance in folding. Mutations at other sites in the protein gave small phi-values, despite the fact that some of them had major effects on stability. The analysis indicates that formation of the antiparallel hairpin is critical and its interactions with the first helix are also important. Thus, the slow folding is not a consequence of the need to fully form the unusual three-stranded beta-sheet in the transition state. Analysis of the urea dependence of the folding rates indicates that mutations modulate the unfolded state. The folding of CTL9 is broadly consistent with the nucleation-condensation model of protein folding.     

Optimal combination of theory and experiment for the characterization of the protein folding landscape of S6: how far can a minimalist model go?

The detailed characterization of the overall free energy landscape associated with the folding process of a protein is the ultimate goal in protein folding studies. Modern experimental techniques provide accurate thermodynamic and kinetic measurements on restricted regions of a protein landscape. Although simplified protein models can access larger regions of the landscape, they are oftentimes built on assumptions and approximations that affect the accuracy of the results. We present a new methodology that allows to combine the complementary strengths of theory and experiment for a more complete characterization of a protein folding landscape. We prove that this new procedure allows a simplified protein model to reproduce remarkably well (correlation coefficient > 0.9) all experimental data available on free energies differences upon single mutations for S6 ribosomal protein and two circular permutants. Our results confirm and quantify the hypothesis, recently formulated on the basis of experimental data, that the folding landscape of protein S6 is strongly affected by an atypical distribution of contact energies.     

Modeling the effects of mutations on the denatured states of proteins.

We develop a model for the reversible denaturation of proteins and for the effects of single-site mutations on the denatured states. The model is based on short chains of sequences of H (hydrophobic) and P (other) monomers configured as self-avoiding walks on the two-dimensional square lattice. The N (native) state is defined as the unique conformation of lowest contact energy, whereas the D (denatured) state is defined as the collection of all other conformations. With this model we are able to determine the exact partition function, and thus the exact native-denatured equilibrium for various solvent conditions, using the computer to exhaustively enumerate every possible configuration. Previous studies confirm that this model shows many aspects of protein-like behavior. The present study attempts to model how the denatured state (1) depends on the amino acid sequence, and (2) is changed by single-site mutations. The model accounts for two puzzling experimental results: (1) the replacement of a polar residue by a hydrophobic amino acid on the surface of a protein can destabilize a native protein, and (2) the "denaturant slope," m = partial delta G/partial c (where c is the concentration of denaturant--urea, guanidine hydrochloride), can sometimes change by as much as 30% due to a single mutation. The principal conclusion of the present study is that, under strong folding conditions, the denatured conformations that are in equilibrium with the native state are not open random configurations. Instead, they are an ensemble of highly compact conformations with a distribution that depends on the residue sequence and that can be substantially altered by single mutations. Most importantly, we conclude that mutations can exert their dominant effects on protein stability by changing the entropy of folding.