Searching for multiple folding pathways of a nearly symmetrical protein: temperature dependent phi-value analysis of the B domain of protein A

The B domain of protein A (BdpA) is a popular paradigm for simulating protein folding pathways. The discrepancies between so many simulations and subsequent experimental testing may be attributable to the protein being highly symmetrical: changing experimental conditions could perturb the subtle interplay between the effects of symmetry in the native structure and the effects of asymmetry from specific interactions in a given sequence. If the protein folds via multiple pathways, perturbations, such as temperature, denaturant concentration, and mutation, should change the flux of micro pathways, leading to changes in the bulk properties of the transition state. We tested this hypothesis by conducting a Phi-analysis of BdpA as a function of temperature from 25.0 degrees C to 60.0 degrees C. The Phi-values had no significant dependence on temperature and the values at 55.0 degrees C (denaturing conditions) are very similar to those at 25.0 degrees C (folding conditions), indicating the structure of the transition state does not significantly change although the experimental conditions are considerably altered. The results suggest that BdpA folds via a single dominant folding pathway.     

Phi-analysis at the experimental limits: mechanism of beta-hairpin formation

The 37-residue Formin-binding protein, FBP28, is a canonical three-stranded beta-sheet WW domain. Because of its small size, it is so insensitive to chemical denaturation that it is barely possible to determine accurately a denaturation curve, as the transition spans 0-7 M guanidinium hydrochloride (GdmCl). It is also only marginally stable, with a free energy of denaturation of just 2.3 kcal/mol at 10 degrees Celsius so only small changes in energy upon mutation can be tolerated. But these properties and relaxation times for folding of 25 micros-400 micros conspire to allow the rapid acquisition of accurate and reproducible kinetic data for Phi-analysis using classical temperature-jump methods. The transition state for folding is highly polarized with some regions having Phi-values of 0 and others 1, as readily seen in chevron plots, with Phi-values of 0 having the refolding arms overlaying and those of 1 the unfolding arms superimposable. Good agreement is seen with transition state structures identified from independent molecular dynamics (MD) simulations at 60, 75, and 100 degrees Celsius, which allows us to explore further the details of the folding and unfolding pathway of FBP28. The first beta-turn is near native-like in the transition state for folding (experimental) and unfolding (MD and experiment). The simulations show that there are transient contacts between the aromatic side-chains of the beta-strands in the denatured state and that these interactions provide the driving force for folding of the first beta-hairpin of this three-stranded sheet. Only after the backbone hydrogen bonds are formed between beta1 and beta2 does a hydrogen bond form to stabilize the intervening turn, or the first beta-turn.     

The structure of the major transition state for folding of an FF domain from experiment and simulation

We have analysed the transition state of folding of the four-helix FF domain from HYPA/FBP11 by high-resolution experiment and simulation as part of a continuing effort to understand the principles of folding and the refinement of predictive methods. The major transition state for folding was subjected to a Phi-value analysis utilising 50 mutants. The transition state contained a nucleus for folding centred around the end of helix 1 (H1) and the beginning of helix 2 (H2). Secondary structure in this region was fully formed (PhiF=0.9-1) and tertiary interactions were well developed. Interactions in the distal part of the native structure were weak (PhiF=0-0.2). The hydrophobic core and other parts of the protein displayed intermediate Phi-values, suggesting that interactions coalesce as the end of H1 and beginning of H2 are in the process of being formed. The distribution of Phi-values resembled that of barnase, which folds via an intermediate, rather than that of CI2 which folds by a concerted nucleation-condensation mechanism. The overall picture of the transition state structure identified in molecular dynamics simulations is in qualitative agreement, with the turn connecting H1 and H2 being formed, a loosened core, and H4 partially unfolded and detached from the core. There are some differences in the details and interpretation of specific Phi-values.     

Amino Acid Insertion Reveals a Necessary Three-Helical Intermediate in the Folding Pathway of the Colicin E7 Immunity Protein Im7

The small (87-residue) α-helical protein Im7 (an inhibitor protein for colicin E7 that provides immunity to cells producing colicin E7) folds via a three-state mechanism involving an on-pathway intermediate. This kinetic intermediate contains three of four native helices that are oriented in a non-native manner so as to minimise exposed hydrophobic surface area at this point in folding. The short (6-residue) helix III has been shown to be unstructured in the intermediate ensemble and does not dock onto the developing hydrophobic core until after the rate-limiting transition state has been traversed. After helix III has docked, it adopts an α-helical secondary structure, and the side chains of residues within this region provide contacts that are crucial to native-state stability. In order to probe further the role of helix III in the folding mechanism of Im7, we created a variant that contains an eight-amino-acid polyalanine-like helix stabilised by a Glu-Arg salt bridge and an Asn-Pro-Gly capping motif, juxtaposed C-terminal to the natural 6-residue helix III. The effect of this insertion on the structure of the native protein and its folding mechanism were studied using NMR and ?-value analysis, respectively. The results reveal a robust native structure that is not perturbed by the presence of the extended helix III. Mutational analysis performed to probe the folding mechanism of the redesigned protein revealed a conserved mechanism involving the canonical three-helical intermediate. The results suggest that folding via a three-helical species stabilised by both native and non-native interactions is an essential feature of Im7 folding, independent of the helical propensity of helix III.     

Structural comparison of the two alternative transition states for folding of TI I27

TI I27, a beta-sandwich domain from the human muscle protein titin, has been shown to fold via two alternative pathways, which correspond to a change in the folding mechanism. Under physiological conditions, TI I27 folds by a classical nucleation-condensation mechanism (diffuse transition state), whereas at extreme conditions of temperature and denaturant it switches to having a polarized transition state. We have used experimental Phi-values as restraints in ensemble-averaged molecular dynamics simulations to determine the ensembles of structures representing the two transition states. The comparison of these ensembles indicates that when native interactions are substantially weakened, a protein may still be able to fold if it can access an alternative transition state characterized by a much larger entropic contribution. Analysis of the probability distribution of Phi-values derived from ensemble averaged simulations, enables us to identify residues that form contacts in some members of the ensemble but not in others illustrating that many interactions present in transition states are not strictly required for the successful completion of the folding process.     

The Folding Mechanism of BBL: Plasticity of Transition-State Structure Observed within an Ultrafast Folding Protein Family

Studies on members of protein families with similar structures but divergent sequences provide insights into the effects of sequence composition on the mechanism of folding. Members of the peripheral subunit-binding domain (PSBD) family fold ultrafast and approach the smallest size for cooperatively folding proteins. Φ-Value analysis of the PSBDs E3BD and POB reveals folding via nucleation-condensation through structurally very similar, polarized transition states. Here, we present a Φ-value analysis of the family member BBL and found that it also folds by a nucleation-condensation mechanism. The mean Φ values of BBL, E3BD, and POB were near identical, indicating similar fractions of non-covalent interactions being formed in the transition state. Despite the overall conservation of folding mechanism in this protein family, however, the pattern of Φ values determined for BBL revealed a larger dispersion of the folding nucleus across the entire structure, and the transition state was less polarized. The observed plasticity of transition-state structure can be rationalized by the different helix-forming propensities of PSBD sequences. The very strong helix propensity in the first helix of BBL, relative to E3BD and POB, appears to recruit more structure formation in that helix in the transition state at the expense of weaker interactions in the second helix. Differences in sequence composition can modulate transition-state structure of even the smallest natural protein domains.     

Different secondary structure elements as scaffolds for protein folding transition states of two homologous four-helix bundles

Comparison of the folding processes for homologue proteins can provide valuable information about details in the interactions leading to the formation of the folding transition state. Here the folding kinetics of 18 variants of yACBP and 3 variants of bACBP have been studied by Phi-value analysis. In combination with Phi-values from previous work, detailed insight into the transition states for folding of both yACBP and bACBP has been obtained. Of the 16 sequence positions that have been studied in both yACBP and bACBP, 5 (V12, I/L27, Y73, V77, and L80) have high Phi-values and appear to be important for the transition state formation in both homologues. Y31, A34, and A69 have high Phi-values only in yACBP, while F5, A9, and I74 have high Phi-values only in bACBP. Thus, additional interactions between helices A2 and A4 appear to be important for the transition state of yACBP, whereas additional interactions between helices A1 and A4 appear to be important for the transition state of bACBP. To examine whether these differences could be assigned to different packing of the residues in the native state, a solution structure of yACBP was determined by NMR. Small changes in the packing of the hydrophobic side-chains, which strengthen the interactions between helices A2 and A4, are observed in yACBP relative to bACBP. It is suggested that different structure elements serve as scaffolds for the folding of the 2 ACBP homologues.     

Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.     

Analysis of the folding pathway of chymotrypsin inhibitor by correlation of phi-values with inter-residue contacts

The transition state for folding of chymotrypsin inhibitor 2 (CI2) is investigated by correlating Phi-values with inter-residue contacts. In agreement with former work, the strongest consolidation of secondary structure is found in the alpha-helix. There are correlations for tertiary structure interactions between the residues Leu49, Ile57 and the helix which have been suggested to represent the main components of the nucleation site of CI2 folding. However, correlations for tertiary structure interactions of comparable magnitude are also found in the helix-strand2-strand1-motif and between strand3and strand4. Copyright 1999 Academic Press.     

The effect of increasing the stability of non-native interactions on the folding landscape of the bacterial immunity protein Im9

How stabilising non-native interactions influence protein folding energy landscapes is currently not well understood: such interactions could speed folding by reducing the conformational search to the native state, or could slow folding by increasing ruggedness. Here, we examine the influence of non-native interactions in the folding process of the bacterial immunity protein Im9, by exploiting our ability to manipulate the stability of the intermediate and rate-limiting transition state (TS) in the folding of this protein by minor alteration of its sequence or changes in solvent conditions. By analysing the properties of these species using Phi-value analysis, and exploration of the structural properties of the TS ensemble using molecular dynamics simulations, we demonstrate the importance of non-native interactions in immunity protein folding and demonstrate that the rate-limiting step involves partial reorganisation of these interactions as the TS ensemble is traversed. Moreover, we show that increasing the contribution to stability made by non-native interactions results in an increase in Phi-values of the TS ensemble without altering its structural properties or solvent-accessible surface area. The data suggest that the immunity proteins fold on multiple, but closely related, micropathways, resulting in a heterogeneous TS ensemble that responds subtly to mutation or changes in the solvent conditions. Thus, altering the relative strength of native and non-native interactions influences the search to the native state by restricting the pathways through the folding energy landscape.     

Folding of a LysM domain: entropy-enthalpy compensation in the transition state of an ideal two-state folder

Protein-engineering methods (Φ-values) were used to investigate the folding transition state of a lysin motif (LysM) domain from Escherichia coli membrane-bound lytic murein transglycosylase D. This domain consists of just 48 structured residues in a symmetrical βααβ arrangement and is the smallest αβ protein yet investigated using these methods. An extensive mutational analysis revealed a highly robust folding pathway with no detectable transition state plasticity, indicating that LysM is an example of an ideal two-state folder. The pattern of Φ-values denotes a highly polarised transition state, with significant formation of the helices but no structure within the β-sheet. Remarkably, this transition state remains polarised after circularisation of the domain, and exhibits an identical Φ-value pattern; however, the interactions within the transition state are uniformly weaker in the circular variant. This observation is supported by results from an Eyring analysis of the folding rates of the two proteins. We propose that the folding pathway of LysM is dominated by enthalpic rather than entropic considerations, and suggest that the lower entropy cost of formation of the circular transition state is balanced, to some extent, by the lower enthalpy of contacts within this structure.     

Proline residues in transmembrane alpha helices affect the folding of bacteriorhodopsin

Proline residues occur frequently in transmembrane alpha helices, which contrasts with their behaviour as helix-breakers in water-soluble proteins. The three membrane-embedded proline residues of bacteriorhodopsin have been replaced individually by alanine and glycine to give P50A, or P50G on helix B, P91A, or P91G on helix C, and P186A or P186G on helix F, and the effect on the protein folding kinetics has been investigated. The rate-limiting apoprotein folding step, which results in formation of a seven transmembrane, alpha helical state, was slower than wild-type protein for the Pro50 and Pro91 mutants, regardless of whether they were mutated to Ala or Gly. These proline residues give rise to several inter-helix contacts, which are therefore important in folding to the seven transmembrane helix state. No evidence for cis-trans isomerisations of the peptidyl prolyl bonds was found during this rate-limiting apoprotein folding step. Mutations of all three membrane-embedded proline residues affected the subsequent retinal binding and final folding to bacteriorhodopsin, suggesting that these proline residues contribute to formation of the retinal binding pocket within the helix bundle, again via helix/helix interactions. These results point to proline residues in transmembrane alpha helices being important in the folding of integral membrane proteins. The helix/helix interactions and hydrogen bonds that arise from the presence of proline residues in transmembrane alpha helices can affect the formation of transmembrane alpha helix bundles as well as cofactor binding pockets.     

Characterisation of the conformational properties of urea-unfolded Im7: implications for the early stages of protein folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.     

Core formation in apomyoglobin: probing the upper reaches of the folding energy landscape

An acid-destabilized form of apomyoglobin, the so-called E state, consists of a set of heterogeneous structures that are all characterized by a stable hydrophobic core composed of 30-40 residues at the intersection of the A, G, and H helices of the protein, with little other secondary structure and no other tertiary structure. Relaxation kinetics studies were carried out to characterize the dynamics of core melting and formation in this protein. The unfolding and/or refolding response is induced by a laser-induced temperature jump between the folded and unfolded forms of E, and structural changes are monitored using the infrared amide I' absorbance at 1648-1651 cm(-1) that reports on the formation of solvent-protected, native-like helix in the core and by fluorescence emission changes from apomyoglobin's Trp14, a measure of burial of the indole group of this residue. The fluorescence kinetics data are monoexponential with a relaxation time of 14 micros. However, infrared kinetics data are best fit to a biexponential function with relaxation times of 14 and 59 micros. These relaxation times are very fast, close to the limits placed on folding reactions by diffusion. The 14 micros relaxation time is weakly temperature dependent and thus represents a pathway that is energetically downhill. The appearance of this relaxation time in both the fluorescence and infrared measurements indicates that this folding event proceeds by a concomitant formation of compact secondary and tertiary structures. The 59 micros relaxation time is much more strongly temperature dependent and has no fluorescence counterpart, indicating an activated process with a large energy barrier wherein nonspecific hydrophobic interactions between helix A and the G and H helices cause some helix burial but Trp14 remains solvent exposed. These results are best fit by a multiple-pathway kinetic model when U collapses to form the various folded core structures of E. Thus, the results suggest very robust dynamics for core formation involving multiple folding pathways and provide significant insight into the primary processes of protein folding.     

Helix capping in the GCN4 leucine zipper.

Capping interactions associated with specific sequences at or near the ends of alpha-helices are important determinants of the stability of protein secondary and tertiary structure. We investigate here the role of the helix-capping motif Ser-X-X-Glu, a sequence that occurs frequently at the N termini of alpha helices in proteins, on the conformation and stability of the GCN4 leucine zipper. The 1.8 A resolution crystal structure of the capped molecule reveals distinct conformations, packing geometries and hydrogen-bonding networks at the amino terminus of the two helices in the leucine zipper dimer. The free energy of helix stabilization associated with the hydrogen-bonding and hydrophobic interactions in this capping structure is -1.2 kcal/mol, evaluated from thermal unfolding experiments. A single cap thus contributes appreciably to stabilizing the terminated helix and thereby the native state. These results suggest that helix capping plays a further role in protein folding, providing a sensitive connector linking alpha-helix formation to the developing tertiary structure of a protein.     

The folding pathway of an FF domain: characterization of an on-pathway intermediate state under folding conditions by (15)N, (13)C(alpha) and (13)C-methyl relaxation dispersion and (1)H/(2)H-exchange NMR spectroscopy

The FF domain from the human protein HYPA/FBP11 folds via a low-energy on-pathway intermediate (I). Elucidation of the structure of such folding intermediates and denatured states under conditions that favour folding are difficult tasks. Here, we investigated the millisecond time-scale equilibrium folding transition of the 71-residue four-helix bundle wild-type protein by (15)N, (13)C(alpha) and methyl(13)C Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiments and by (1)H/(2)H-exchange measurements. The relaxation data for the wild-type protein fitted a simple two-site exchange process between the folded state (F) and I. Destabilization of F in mutants A17G and Q19G allowed the detection of the unfolded state U by (15)N CPMG relaxation dispersion. The dispersion data for these mutants fitted a three-site exchange scheme, U<-->I<-->F, with I populated higher than U. The kinetics and thermodynamics of the folding reaction were obtained via temperature and urea-dependent relaxation dispersion experiments, along with structural information on I from backbone (15)N, (13)C(alpha) and side-chain methyl (13)C chemical shifts, with further information from protection factors for the backbone amide groups from (1)H/(2)H-exchange. Notably, helices H1-H3 are at least partially formed in I, while helix H4 is largely disordered. Chemical shift differences for the methyl (13)C nuclei suggest a paucity of stable, native-like hydrophobic interactions in I. These data are consistent with Phi-analysis of the rate-limiting transition state between I and F. The combination of relaxation dispersion and Phi data can elucidate whole experimental folding pathways.     

Mechanism of protein folding

The high structural resolution of the main transition states for the formation of native structure for the six small proteins of which Phi-values for a large set of mutants have become available, barstar, barnase, chymotrypsin inhibitor 2, Arc repressor, the src SH3 domain, and a tetrameric p53 domain reveals that for the first 5 of these proteins: (1) Residues that belong to regular secondary structure have a significantly larger average fraction of native structural consolidation than residues in loops; (2) on the other hand, secondary and tertiary structures have built up to the same degree, or at least a high degree, but nonuniformly distributed over the molecule; (3) the most consolidated parts of each protein molecule in the transition state cluster together, and these clusters contain a significantly higher percentage of residues that belong to regular secondary structure than the rest of the molecule. These observations further reconcile the framework model with the nucleation-condensation mechanism for folding: The amazing speed of protein folding can be understood as caused by the catalytic effect of the formation of clusters of residues which have particularly high preferences for the early formation of regular secondary structure in the presence of significant amounts of tertiary structure interactions.     

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.     

Characterization of the folding kinetics of a three-helix bundle protein via a minimalist Langevin model.

We use a simple off-lattice Langevin model of protein folding to characterize the folding and unfolding of a fast-folding, 46 residue three-helix bundle. Under conditions at which the C-terminal helix is 30 % stable, we observe a clear three-state folding mechanism. In the on-pathway intermediate state, the middle and C-terminal helices are folded and in contact with each other, while the N-terminal region remains disordered. Nevertheless, under these conditions this intermediate is thermodynamically unstable relative to its unfolded state. The first and highest folding barrier corresponds to the organization of the hinge between the middle and C-terminal helices. A subsequent major barrier corresponds to the organization of the hinge between the middle and N-terminal helices. Hyperstabilizing the hinge regions leads to twice the folding rate that is obtained from hyperstabilizing the helices, even though much fewer contacts are involved in hinge hyperstabilization than in helix hyperstabilization. Unfolding follows single-exponential kinetics, even at temperatures only slightly above the folding transition temperature.     

Protein folding simulation with solvent-induced force field: folding pathway ensemble of three-helix-bundle proteins

We propose a coarse-grained model of proteins that take into account solvent effects and apply it for simulating folding of a three-helix-bundle protein. The energy functional form, refined from our previous work (Takada et al., J Chem Phys 1999;110:11616-11629), tries to closely imitate real physico-chemical interactions. In particular, the hydrogen bond that depends on local dielectric constant, the helix capping effect, and side-chain entropic effects are included. With use of the model, we simulate folding of the GA module of an albumin binding domain, 1prb(7-53), finding most trajectories reach at the native topology within 1 micros. In the simulation, helices 1 and 3 are mostly formed earlier accompanied by non-specific collapse, while second helix is intrinsically less stable and is formed with the help of tertiary contacts at later stage. We compute an analog of the transition state ensemble and compare it with those of other three-helix-bundle proteins. The transition state of 1prb(7-53) includes a few specific tertiary contacts of C terminus of helix 3 with the loop region between helices 1 and 2. This resembles, but is not equivalent to, an early formed region of fragment B of staphylococcal protein A, but is quite different from the folding transient structures of a de novo designed three-helix-bundle peptide.