Energy landscape theory is a powerful tool for understanding the structure and dynamics of complex molecular systems, in particular biological macromolecules1. The primary sequence of a protein defines its free energy landscape, and thus determines the folding pathway and the rate constants of folding and unfolding, as well as its native structure. Theory has shown that roughness in the energy landscape will lead to slower folding1, but derivation of detailed experimental descriptions of this landscape is challenging. Simple folding models2,3 show that folding is significantly influenced by chain entropy; proteins where the contacts are local fold fast, due to the low entropy cost of forming stabilising, native contacts during folding4,5. For some protein families, stability is also a determinant of folding rate constants6. Where these simple metrics fail to predict folding behaviour it is probable that there are features in the energy landscape that are unusual. Such general observations cannot explain the folding behaviour of the R15, R16 and R17 domains of α-spectrin. R15 folds ~3000 times faster than its homologues, although they have similar structures, stabilities and, as far as can be determined, transition state stabilities7-10. Here we show that landscape roughness (internal friction) is responsible for the slower folding and unfolding of R16 and R17. We use chimeric domains to demonstrate that this internal friction is a property of the core, and suggest that frustration in the landscape of the slow folding spectrin domains may be due to mis-docking of the long helices during folding. Although theoretical studies have suggested that rugged landscapes will result in slower folding, this is the first time that such a phenomenon has been shown experimentally to directly influence the folding kinetics of a “normal” protein with a significant energy barrier – one which folds on a relatively slow ms-s timescale.

The elongated three-helix bundle domains spectrin R16 and R17 fold some two to three orders of magnitude more slowly than their homologue R15. We have shown that this slow folding is due, at least in part, to roughness in the free-energy landscape of R16 and R17. We have proposed that this roughness is due to a frustrated search for the correct docking of partly preformed helices. However, this accounts for only a small part of the slowing of folding and unfolding. Five residues on the A helix of R15, when inserted together into R16 or R17, increase the folding rate constants, reduce landscape roughness, and alter the folding mechanism to one resembling R15. The effect of each of these mutations individually is investigated here. No one mutation causes the behavior seen for the five in combination. However, two mutations, E18F and K25V, significantly increase the folding and unfolding rates of both R16 and R17 but without a concomitant loss in landscape roughness. E18F has the greatest effect on the kinetics, and a Φ-value analysis of the C helix reveals that the folding mechanism is unchanged. For both E18F and K25V the removal of the charge and resultant transition state stabilization is the main origin of the faster folding. Consequently, the major cause of the unusually slow folding of R16 and R17 is the non-native burial of the two charged residues in the transition state. The slowing due to landscape roughness is only about fivefold.free energy landscape | frustration | phi value | protein folding T he 15th, 16th, and 17th domains of chicken brain α-spectrin (R15, R16, and R17) have very similar structures, stabilities, and β-Tanford (β T ) values (which reflect the compactness of the transition state for folding and unfolding) (1-6). However, the folding of R15 differs from that of R16 and R17 in a number of respects. R15 folds and unfolds two orders of magnitude faster than R16 and three orders of magnitude faster than R17. R16 and R17 have two sequential transition states, and for both domains the first of these (TS1) shows significant landscape roughness (or "internal friction") (7). This has not been seen for any other domain of comparable size and folding kinetics, although theory has long predicted the possibility of such landscape roughness (8-15). R15 has a broad transition state (characterized by "rollover" in the unfolding limb for some mutants and for wild type under some conditions), but due to the speed of folding and unfolding it is not known whether it has two sequential transition states (16). However, the early transition state of R15 (which corresponds to TS1 of R16 and R17 and will be referred to as such) has a smoother, less frustrated landscape (7).Φ-value analysis shows that for all three domains the A and C helices are partially structured at TS1, whereas the B helix is relatively unstructured (16-18). R16 and R17 fold via a framework-like mechanism, with some tertiary contacts formed but more extensive secondary structure that extends throughout the A and C helices. In contrast, R1...

Interdomain interactions of spectrin are critical for maintenance of the erythrocyte cytoskeleton. In particular, "head-to-head" dimerization occurs when the intrinsically disordered C-terminal tail of β-spectrin binds the N-terminal tail of α-spectrin, folding to form the "spectrin tetramer domain". This non-covalent three-helix bundle domain is homologous in structure and sequence to previously studied spectrin domains. We find that this tetramer domain is surprisingly kinetically stable. Using a protein engineering Φ-value analysis to probe the mechanism of formation of this tetramer domain, we infer that the domain folds by the docking of the intrinsically disordered β-spectrin tail onto the more structured α-spectrin tail.

In this study, the contribution of intramembrane hydrogen bonding at the interface between polypeptide and cofactor is explored in the native lipid environment by use of model bacteriochlorophyll proteins. In the peripheral antenna complex, LH2, large portions of the transmembrane helices, which make up the dimeric bacteriochlorophyll-binding site, are replaced by simplified, alternating alanine-leucine stretches. Replacement of either one of the two helices with the helices containing the model sequence at a time results in the assembly of complexes with nearly native light harvesting properties. In contrast, replacement of both helices results in the loss of antenna complexes from the membrane. The assembly of such doubly modified complexes is restored by a single intramembrane serine residue at position ؊4 relative to the liganding histidine of the ␣-subunit. In situ analysis of the spectral properties in a series of site-directed mutants reveals a critical dependence of the model complex assembly on the side chain of the residue at this position in the helix. A hydrogen bond between the hydroxy group of the serine and the 13 1 keto group of one of the central bacteriochlorophylls of the complexes is identified by Raman spectroscopy in the model antenna complex containing one of the alanine-leucine helices. The additional OH group of the serine residue, which participates in hydrogen bonding, increases the thermal stability of the model complexes in the native membrane. Intramembrane hydrogen bonding is thus shown to be a key factor for the binding of bacteriochlorophyll and assembly of this model cofactor-polypeptide site.

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.