Crystal structures and increased stabilization of the protein G variants with switched folding pathways NuG1 and NuG2

We recently described two protein G variants (NuG1 and NuG2) with redesigned first hairpins that were almost twice as stable, folded 100-fold faster, and had a switched folding mechanism relative to the wild-type protein. To test the structural accuracy of our design algorithm and to provide insights to the dramatic changes in the kinetics and thermodynamics of folding, we have now determined the crystal structures of NuG1 and NuG2 to 1.8 A and 1.85 A, respectively. We find that they adopt hairpin structures that are closer to the computational models than to wild-type protein G; the RMSD of the NuG1 hairpin to the design model and the wild-type structure are 1.7 A and 5.1 A, respectively. The crystallographic B factor in the redesigned first hairpin of NuG1 is systematically higher than the second hairpin, suggesting that the redesigned region is somewhat less rigid. A second round of structure-based design yielded new variants of NuG1 and NuG2, which are further stabilized by 0.5 kcal/mole and 0.9 kcal/mole.     

Intermediates and the folding of proteins L and G

We use a minimalist protein model, in combination with a sequence design strategy, to determine differences in primary structure for proteins L and G, which are responsible for the two proteins folding through distinctly different folding mechanisms. We find that the folding of proteins L and G are consistent with a nucleation-condensation mechanism, each of which is described as helix-assisted beta-1 and beta-2 hairpin formation, respectively. We determine that the model for protein G exhibits an early intermediate that precedes the rate-limiting barrier of folding, and which draws together misaligned secondary structure elements that are stabilized by hydrophobic core contacts involving the third beta-strand, and presages the later transition state in which the correct strand alignment of these same secondary structure elements is restored. Finally, the validity of the targeted intermediate ensemble for protein G was analyzed by fitting the kinetic data to a two-step first-order reversible reaction, proving that protein G folding involves an on-pathway early intermediate, and should be populated and therefore observable by experiment.     

The role of sequence and structure in protein folding kinetics; the diffusion-collision model applied to proteins L and G

The diffusion-collision model (DCM) is applied to the folding kinetics of protein L and protein G. In the DCM, the two proteins are treated as consisting of two beta-hairpins and one alpha-helix, so that they are isomorphous with the three-helix bundle DCM model. In the absence of sequence dependent factors, both proteins would fold in the same way in the DCM, with the coalescence of the N-terminal hairpin and the helix slightly favored over the C-terminal hairpin and the helix because the former are closer together than the latter. However, sequence dependent factors make the N-terminal hairpin of protein L and the C-terminal hairpin of protein G more stable in the ensemble of unfolded conformations. This difference in the stabilities gives rise to the difference in the calculated folding behavior, in agreement with experiment.     

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

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

Commitment and nucleation in the protein G transition state

An accurate characterization of the transition state ensemble (TSE) is central to furthering our understanding of the protein folding reaction. We have extensively tested a recently reported method for studying a protein's TSE, utilizing phi-value data from protein engineering experiments and computational studies as restraints in all-atom Monte Carlo (MC) simulations. The validity of interpreting experimental phi-values as the fraction of native contacts made by a residue in the TSE was explored, revealing that this definition is unable to uniquely specify a TSE. The identification of protein G's second hairpin, in both pre and post-transition conformations demonstrates that high experimental phi-values do not guarantee a residue's importance in the TSE. An analysis of simulations based on structures restrained by experimental phi-values is necessary to yield this result, which is not obvious from a simplistic interpretation of individual phi-values. The TSE that we obtain corresponds to a single, specific nucleation event, characterized by six residues common to all three observed, convergent folding pathways. The same specific nucleus was independently identified from computational and experimental data, and "Conservation of Conservation" analysis in the protein G fold. When associated strictly with complete nucleus formation and concomitant chain collapse, folding is a well-defined two state event. Once the nucleus has formed, the folding reaction enters a slow relaxation process associated with side-chain packing and small, local backbone rearrangements. A detailed analysis of phi-values and their relationship to the transition state ensemble allows us to construct a unified theoretical model of protein G folding.     

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.     

Roles of physical interactions in determining protein-folding mechanisms: molecular simulation of protein G and alpha spectrin SH3

Protein-folding mechanisms of two small globular proteins, IgG binding domain of protein G and alpha spectrin SH3 domain are investigated via Brownian dynamics simulations with a model made of coarse-grained physical energy functions responsible for sequence-specific interactions and weak Gō-like energies. The folding pathways of alpha spectrin SH3 are known to be mainly controlled by the native topology, while protein G folding is anticipated to be more sensitive to the sequence-specific effects than native topology. We found in the folding of protein G that the C terminal beta hairpin is formed earlier and is rigid, once ordered, in the presence of an intact C terminal turn. The alpha helix is found to exhibit repeated partial formations/deformations during folding and to be stabilized via the tertiary contact with preformed beta sheets. This predicted scenario is fully consistent with experimental phi value data. Moreover, we found that the folding route is critically affected when the hydrophobic interaction is excluded from physical energy terms, suggesting that the hydrophobicity critically contributes to the folding propensity of protein G. For the folding of alpha spectrin SH3, we found that the distal beta hairpin and diverging turn are parts formed early, fully in harmony with previous results of simple Gō-like and experimental analysis, supporting that the folding route of SH3 domain is robust and coded by the native topology. The hybrid method provides useful tools for analyzing roles of physical interactions in determining folding mechanisms.     

A molecular dynamics simulation study of segment B1 of protein G

The immunoglobulin binding protein, segment B1 of protein G, has been studied experimentally as a paradigm for protein folding. This protein consists of 56 residues, includes both β sheet and α helix and contains neither disulfide bonds nor proline residues. We report an all-atom molecular dynamics study of the native manifold of the protein in explicit solvent. A 2-ns simulation starting from the nuclear magnetic resonance (NMR) structure and a 1-ns control simulation starting from the x-ray structure were performed. The difference between average structures calculated over the equilibrium portion of trajectories is smaller than the difference between their starting conformations. These simulation averages are structurally similar to the x-ray structure and differ in systematic ways from the NMR-determined structure. Partitioning of the fluctuations into fast (<20 ps) and slow (<20 ps) components indicates that the β sheet displays greater long-time mobility than does the α helix. Clore and Gronenborn [J. Mol. Biol. 223:853–856, 1992] detected two long-residence water molecules by NMR in a solution structure of segment B1 of protein G. Both molecules were found in the fully exposed regions and were proposed to be stabilized by bifurcated hydrogen bonds to the protein backbone. One of these long-residence water molecules, found near an exposed loop region, is identified in both of our simulations, and is seen to be involved in the formation of a stable water-mediated hydrogen bond bridge. The second water molecule, located near the middle of the α helix, is not seen with an exceptional residence time in either as a result of the conformation being closer to the x-ray structure in this region of the protein. Proteins 29:193–202, 1997. © 1997 Wiley-Liss, Inc.     

Importance of sequence specificity for predicting protein folding pathways: perturbed Gaussian chain model

Recent experimental and theoretical studies suggest that rates and pathways of protein folding are largely decided by topology of the native structures, at least for small proteins. However, some exceptions are known; for example, protein L and protein G have the same topology, but exhibit different characteristics of the TSE. Thus, folding pathways of some proteins are critically affected by detailed information on amino acid sequences. To investigate the sequence specificity, we calculate folding pathways of 20 small proteins using the perturbed Gaussian chain model developed by Portman et al. (Phys Rev Lett 1998;81:5237-5240; J Chem Phys 2001;114:5069-5081). Characteristics of the TSE predicted by the model are in good agreement with experimental phi-value data for many proteins at coarse-grained level. Especially, estimation of folding TSE for protein G and protein L based on both topology and additional sequence information are consistent with experimental phi-value data. With only topology information, however, the model predicts the TSE of protein G incorrectly. Moreover, the model that uses topology and sequence information describes free energy profiles of two-state and three-state folders consistently with experiment, whereas the topology only model predicts free energy profiles of some proteins incorrectly. This indicates that sequence specificity also has critical roles in determining the folding pathways for some proteins.     

Engineering stabilising beta-sheet interactions into a conformationally flexible region of the folding transition state of ubiquitin

Protein engineering studies suggest that the transition state for the folding of ubiquitin is highly polarised towards the N-terminal part of the sequence and involves a nucleus of residues within the beta-hairpin (residues 1-17) and main alpha-helix (residues 23-34). In contrast, the observation of small phi-values for residues in the C-terminal portion of the sequence (residues 35-76), coupled with a folding topology that results in a much higher contact order, suggests that fast folding of ubiquitin is dependent upon configurational flexibility in the C-terminal part of the polypeptide chain to ensure passage down a relatively smooth folding funnel to the native state. We show that the introduction of a small mini-hairpin motif as an extension of the native 43-50 hairpin stabilises local interactions in the C-terminal part of the sequence, resulting largely in a deceleration of the unfolding kinetics without perturbing the apparent two-state folding mechanism. However, a single-point Leu-->Phe substitution within the engineered hairpin sequence leads to the premature collapse of the denatured ensemble through the stabilisation of non-native interactions and the population of a compact intermediate. Non-linear effects in the kinetic data at low concentrations of denaturant suggest that the collapsed state, which is further stabilised in the presence of cosmotropic salts, may subsequently fold directly to the native state through a "triangular" reaction scheme involving internal rearrangement rather than unfolding and refolding.     

The role of the turn in beta-hairpin formation during WW domain folding

The folding of WW domains is rate limited by formation of a beta-hairpin comprising residues from strands 1 and 2. Residues in the turn of this hairpin have reported Phi-values for folding close to 1 and have been proposed to nucleate folding. High Phi-values do not necessarily imply that the energetics of formation are a driving force for initiating folding. We demonstrate by NMR studies and molecular dynamics simulations that the first turn of the hYAP, FBP28, and PIN1 WW domains is structurally dynamic and solvent exposed in the native and folding transition states. It is, therefore, unlikely that the formation of the beta-turn per se provides the energetic driving force for hairpin folding. It is more likely that the turn acts as an easily formed hinge that facilitates the formation of the hairpin; it is a nucleus as defined by the nucleation-condensation mechanism whereby a diffuse nucleus is stabilized by associated interactions.     

Removal of kinetic traps and enhanced protein folding by strategic substitution of amino acids in a model alpha-helical hairpin peptide

The presence of non-native kinetic traps in the free energy landscape of a protein may significantly lengthen the overall folding time so that the folding process becomes unreliable. We use a computational model alpha-helical hairpin peptide to calculate structural free energy landscapes and relate them to the kinetics of folding. We show how protein engineering through strategic changes in only a few amino acid residues along the primary sequence can greatly increase the speed and reliability of the folding process, as seen experimentally. These strategic substitutions also prevent the formation of long-lived misfolded configurations that can cause unwanted aggregations of peptides. These results support arguments that removal of kinetic traps, obligatory or nonobligatory, is crucial for fast folding.     

Measuring cooperativity in the formation of a three-stranded beta sheet (double hairpin)

Recently validated chemical shift measures of hairpin structuring have been applied to a series of turn mutants of the Schenck-Gellman three-strand beta-sheet model with the aim of measuring the entropic advantage associated with aligning an additional strand onto an existing hairpin versus aligning the same two strands in an initial hairpin formation. In a four-state analysis (unfolded, 2 single hairpins, and the double hairpin fold in equilibrium) a cooperativity index can be defined as the factor by which the equilibrium constant for hairpin formation is improved when one strand is prestructured. This cooperativity index is 2.7 +/- 0.7 for hairpin formation about a stable D-Pro-Gly turn locus and increases to 7.6 +/- 1.2 for an Asn-Gly turn locus. The latter corresponds to a cooperativity induced DeltaDeltaG increment of 4.9 kJ/mol for the folding of a hairpin. Although larger than previous experimental measures of folding cooperativity in three-stranded sheets, the magnitude of this effect (which is considerably less than the TDeltaDeltaS expectation for prestructuring three or more beta-strand residue sites) likely reflects the intrinsic preference of these designed sequences for extended conformations. If similar or larger effects apply to protein beta-sheet folding, it is not surprising that particularly favorable hairpin alignments serve as nucleation sites in protein folding pathways.     

Critical role of beta-hairpin formation in protein G folding

Comparison of the folding mechanisms of proteins with similar structures but very different sequences can provide fundamental insights into the determinants of protein folding mechanisms. Despite very little sequence similarity, the approximately 60 residue IgG binding domains of protein G and protein L both consist of a single helix packed against a four-stranded sheet formed by two symmetrically disposed beta-hairpins. We demonstrate that, as in the case of protein L, one of the two beta-turns of protein G is formed and the other disrupted in the folding transition state. Unlike protein L, however, in protein G it is the second beta-turn that is formed in the folding transition state ensemble. Substitution of an Asp residue by Ala in protein G that eliminates an i,i+2 side chain-main chain hydrogen bond in the second beta-turn slows the folding rate approximately 20-fold but has virtually no effect on the unfolding rate. Taken together with previous results, these findings suggest that the presence of an intact beta-turn in the folding transition state is a consequence of the overall topology of protein L and protein G, but the particular hairpin that is formed is determined by the detailed interatomic interactions that determine the free energies of formation of the isolated beta-hairpins.     

Cooperativity in two-state protein folding kinetics

We present a solvable model that predicts the folding kinetics of two-state proteins from their native structures. The model is based on conditional chain entropies. It assumes that folding processes are dominated by small-loop closure events that can be inferred from native structures. For CI2, the src SH3 domain, TNfn3, and protein L, the model reproduces two-state kinetics, and it predicts well the average Phi-values for secondary structures. The barrier to folding is the formation of predominantly local structures such as helices and hairpins, which are needed to bring nonlocal pairs of amino acids into contact.     

Protein folding forces

We investigate the average inter-residue folding forces derived from mutational data of the 15 proteins: barstar, barnase, chymotrypsin inhibitor 2 (CI2), Src SH3 domain, spectrin R16 domain, Arc repressor, apo-azurin, cold shock protein B (cspB), C-terminal domain of ribosomal protein L9 (CTL9), FKBP12, α-lactalbumin, colicin E7 immunity protein 7 (IM7), colicin E9 immunity protein 9 (IM9), spectrin R17 domain, and ubiquitin. The residue-specific contributions to folding in most of the 15 protein molecules are highly non-uniformly distributed and are typically about 1 piconewton (pN) per interaction. The strongest folding forces often occur in some of the helices and strands of folding nuclei which suggests that folding nucleation−condensation is partially directed by formation of some secondary structure interactions. The correlation of the energy changes of mutants with inter-residue contact maps of the protein molecules provides a higher resolution than assigning the mutant data to certain positions in the polypeptide strand alone. In contrast to previous Φ-value analysis, we now can partially resolve folding motions. Compaction of at least one α-helix along its axis mediated by internal hydrogen bonds and stabilized by diffuse tertiary structure interactions appears to be one important molecular event during early folding in barstar, CI2, spectrin R16 domain, Arc repressor, α-lactalbumin, IM7, IM9, and spectrin R17 domain. A lateral movement of at least two strands neighbored in sequence towards each other appears to be involved in early folding of the SH3 domain, cspB, CTL9, and FKBP12.     

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.     

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.