Using chimeric immunity proteins to explore the energy landscape for alpha-helical protein folding

Direct examination of the relevance for folding, binding and electron transfer of a conserved protein folding intermediate

Near the minimum free energy basin of proteins where the native ensemble resides, partly unfolded conformations of slightly higher energy can be significantly populated under native conditions. It has been speculated that they play roles in molecular recognition and catalysis, but they might represent contemporary features of the evolutionary process without functional relevance. Obtaining conclusive evidence on these alternatives is difficult because it requires comparing the performance of a given protein when populating and when not populating one such intermediate, in otherwise identical conditions. Wild type apoflavodoxin populates under native conditions a partly unfolded conformation (10% of molecules) whose unstructured region includes the binding sites for the FMN cofactor and for redox partner proteins. We recently engineered a thermostable variant where the intermediate is no longer detectable. Using the wild type and variant, we assess the relevance of the intermediate comparing folding kinetics, cofactor binding kinetics, cofactor affinity, X-ray structure, intrinsic dynamics, redox potential of the apoflavodoxin-cofactor complex (Fld), its affinity for partner protein FNR, and electron transfer rate within the Fld/FNR physiological complex. Our data strongly suggest the intermediate state, conserved in long-chain apoflavodoxins, is not required for the correct assembly of flavodoxin nor does it contribute to shape its electron transfer properties. This analysis can be applied to evaluate other native basin intermediates.

Protein folding, misfolding and aggregation: The importance of two-electron stabilizing interactions

Proteins associated with neurodegenerative diseases are highly pleiomorphic and may adopt an all-α-helical fold in one environment, assemble into all-β-sheet or collapse into a coil in another, and rapidly polymerize in yet another one via divergent aggregation pathways that yield broad diversity of aggregates’ morphology. A thorough understanding of this behaviour may be necessary to develop a treatment for Alzheimer’s and related disorders. Unfortunately, our present comprehension of folding and misfolding is limited for want of a physicochemical theory of protein secondary and tertiary structure. Here we demonstrate that electronic configuration and hyperconjugation of the peptide amide bonds ought to be taken into account to advance such a theory. To capture the effect of polarization of peptide linkages on conformational and H-bonding propensity of the polypeptide backbone, we introduce a function of shielding tensors of the C^α atoms. Carrying no information about side chain-side chain interactions, this function nonetheless identifies basic features of the secondary and tertiary structure, establishes sequence correlates of the metamorphic and pH-driven equilibria, relates binding affinities and folding rate constants to secondary structure preferences, and manifests common patterns of backbone density distribution in amyloidogenic regions of Alzheimer’s amyloid β and tau, Parkinson’s α-synuclein and prions. Based on those findings, a split-intein like mechanism of molecular recognition is proposed to underlie dimerization of Aβ, tau, αS and PrP^C, and divergent pathways for subsequent association of dimers are outlined; a related mechanism is proposed to underlie formation of PrP^Sc fibrils. The model does account for: (i) structural features of paranuclei, off-pathway oligomers, non-fibrillar aggregates and fibrils; (ii) effects of incubation conditions, point mutations, isoform lengths, small-molecule assembly modulators and chirality of solid-liquid interface on the rate and morphology of aggregation; (iii) fibril-surface catalysis of secondary nucleation; and (iv) self-propagation of infectious strains of mammalian prions.

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

Protein destabilisation in ionic liquids: the role of preferential interactions in denaturation

The preferential binding of anions and cations in aqueous solutions of the ionic liquids (ILs) 1-butyl-3-methylimidazolium ([C4mim](+)) and 1-ethyl-3-methylimidazolium ([C2mim](+)) chloride and dicyanamide (dca(-)) with the small alpha-helical protein Im7 was investigated using a combination of differential scanning calorimetry, NMR spectroscopy and molecular dynamics (MD) simulations. Our results show that direct ion interactions are crucial to understand the effects of ILs on the stability of proteins and that an anion effect is dominant. We show that the binding of weakly hydrated anions to positively charged or polar residues leads to the partial dehydration of the backbone groups, and is critical to control stability, explaining why dca(-) is more denaturing than Cl(-). Direct cation-protein interactions also mediate stability; cation size and hydrophobicity are relevant to account for destabilisation as shown by the effect of [C4mim](+) compared to [C2mim](+). The specificity in the interaction of IL ions with protein residues established by weak favourable interactions is confirmed by NMR chemical shift perturbation, amide hydrogen exchange data and MD simulations. Differences in specificity are due to the balance of interaction established between ion pairs and ion-solvent that determine the type of residues affected. When the interaction of both cation and anion with the protein is strong the net result is similar to a non-specific interaction, leading ultimately to unfolding. Since the nature of the ions is a determinant of the level of interaction with the protein towards denaturation or stabilisation, ILs offer a unique possibility to modulate protein stabilisation or even folding events.

Understanding how small helical proteins fold: conformational dynamics of Im proteins relevant to their folding landscapes

Understanding the mechanism of folding of small proteins requires characterization of their starting unfolded states and any partially unfolded states populated during folding. Here, we review what is known from NMR about these states of Im7, a 4-helix bundle protein that folds via an on-pathway intermediate, and show that there is an alignment of non-native structure in urea-unfolded Im7 with the helices of native Im7 that is a consequence of hydrophobic helix-promoting residues also promoting cluster-formation in the unfolded protein. We suggest that this kind of alignment is present in other proteins and is relevant to how native state topology determines folding rates.

Nuclease colicins and their immunity proteins

It is more than 80 years since Gratia first described 'a remarkable antagonism between two strains of Escherichia coli'. Shown subsequently to be due to the action of proteins (or peptides) produced by one bacterium to kill closely related species with which it might be cohabiting, such bacteriocins have since been shown to be commonplace in the internecine warfare between bacteria. Bacteriocins have been studied primarily from the twin perspectives of how they shape microbial communities and how they penetrate bacteria to kill them. Here, we review the modes of action of a family of bacteriocins that cleave nucleic acid substrates in E. coli, known collectively as nuclease colicins, and the specific immunity (inhibitor) proteins that colicin-producing organisms make in order to avoid committing suicide. In a process akin to targeting in mitochondria, nuclease colicins engage in a variety of cellular associations in order to translocate their cytotoxic domains through the cell envelope to the cytoplasm. As well as informing on the process itself, the study of nuclease colicin import has also illuminated functional aspects of the host proteins they parasitize. We also review recent studies where nuclease colicins and their immunity proteins have been used as model systems for addressing fundamental problems in protein folding and protein-protein interactions, areas of biophysics that are intimately linked to the role of colicins in bacterial competition and to the import process itself.

NMR characterisation of the relationship between frustration and the excited state of Im7

Previous work shows that Im9 folds in a two-state transition while its homologue Im7 folds in a three-state transition via an on-pathway kinetic intermediate state (KIS), with this difference being related to frustration in the structure of Im7. We have used NMR spectroscopy to study conformational dynamics connected to the frustration. A combination of equilibrium peptide N(1)H/N(2)H exchange, model-free analyses of backbone NH relaxation data and relaxation dispersion (RD)-NMR shows that the native state of Im7 is in equilibrium with an intermediate state that is lowly populated [equilibrium intermediate state (EIS)]. Comparison of kinetic and thermodynamic parameters describing the EIS native-state equilibrium obtained by RD-NMR with previously reported parameters describing the KIS native-state equilibrium obtained from stopped-flow fluorescence studies of refolding His-tagged Im7 shows that the KIS and the EIS are the same species. (15)N chemical shifts of the EIS obtained from the RD-NMR analysis show that residues forming helix III in the native state are unstructured in the EIS while other residues experiencing frustration in the native state are in structured regions of the EIS. We show that binding of Im7 and its L53A/I54A variant (which resembles the EIS as shown in previous work) to the cognate partner for Im7, the DNase domain of colicin E7, causes the dynamic processes associated with the frustration to be dampened.

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.

Computer simulations of the refolding of sperm whale apomyoglobin from high-temperature denaturated state

The refolding mechanism of apomyoglobin (apoMb) subsequent to high-temperature unfolding has been examined using computer simulations with atomic level detail. The folding of this protein has been extensively studied experimentally, providing a large database of folding parameters which can be probed using simulations. In the present study, 4-folding trajectories of apoMb were computed starting from coiled structures. A crystal structure of sperm whale myoglobin taken from the Protein Data Bank was used to construct the final native conformation by removal of the heme group followed by energy optimization. The initial unfolded conformations were obtained from high-temperature molecular dynamics simulations. Room-temperature refolding trajectories at neutral pH were obtained using the stochastic difference equation in length algorithm. The folding trajectories were compared with experimental results and two previous molecular dynamics studies at low pH. In contrast to the previous simulations, an extended intermediate with large helical content was not observed. In the present study, a structural collapse occurs without formation of helices or native contacts. Once the protein structure is more compact (radius of gyration<18 A) secondary and tertiary structures appear. These results suggest that apoMb follows a different folding pathway after high-temperature denaturation.

Identification and characterization of protein folding intermediates

In order to understand the mechanism by which a polypeptide chain folds into its functionally active native state it is necessary to characterize in detail all the species accumulated along the pathway. The elusive nature of protein folding intermediates poses their identification and characterization as an extremely difficult task in the protein folding field. In the case of small single domain proteins, the direct measurement of the thermodynamics and structural parameters of protein folding intermediates has provided new insights on the nature of the forces involved in the stabilization of nascent protein structures. Here we summarize some of the experimental approaches aimed at the detection and characterization of folding intermediates along with a discussion of some general structural features emerging from these studies.

Engineering of the Escherichia coli Im7 immunity protein as a loop display scaffold

Protein scaffolds derived from non-immunoglobulin sources are increasingly being adapted and engineered to provide unique binding molecules with a diverse range of targeting specificities. The ColE7 immunity protein (Im7) from Escherichia coli is potentially one such molecule, as it combines the advantages of (i) small size, (ii) stability conferred by a conserved four anti-parallel alpha-helical framework and (iii) availability of variable surface loops evolved to inactivate members of the DNase family of bacterial toxins, forming one of the tightest known protein-protein interactions. Here we describe initial cloning and protein expression of Im7 and its cognate partner the 15 kDa DNase domain of the colicin E7. Both proteins were produced efficiently in E.coli, and their in vitro binding interactions were validated using ELISA and biosensor. In order to assess the capacity of the Im7 protein to accommodate extensive loop region modifications, we performed extensive molecular modelling and constructed a series of loop graft variants, based on transfer of the extended CDR3 loop from the IgG1b12 antibody, which targets the gp120 antigen from HIV-1. Loop grafting in various configurations resulted in chimeric proteins exhibiting retention of the underlying framework conformation, as measured using far-UV circular dichroism spectroscopy. Importantly, there was low but measurable transfer of antigen-specific affinity. Finally, to validate Im7 as a selectable scaffold for the generation of molecular libraries, we displayed Im7 as a gene 3 fusion protein on the surface of fd bacteriophages, the most common library display format. The fusion was successfully detected using an anti-Im7 rabbit polyclonal antibody, and the recombinant phage specifically recognized the immobilized DNase. Thus, Im7 scaffold is an ideal protein display scaffold for the future generation and for the selection of libraries of novel binding proteins.

Sulfate-Induced Effects in the On-Pathway Intermediate of the Bacterial Immunity Protein Im7*

Intermediates have now been identified in the folding of a number of small, single-domain proteins. Here we describe experiments to determine the effect of Na(2)SO(4) on the properties of the on-pathway intermediate formed early during the folding of the four-helical protein, Im7. This intermediate, studied previously in 0.4 M Na(2)SO(4), contains three of the four native helices and is fascinating in that several residues in helices I, II, and IV make non-native interactions that stabilize this state. Whether these contacts form as a consequence of the presence of Na(2)SO(4), however, remained unresolved. Using kinetic analysis of the effect of Na(2)SO(4) on the unfolding and refolding kinetics of Im7*, combined with detailed analysis of the resulting chevron plots, we show that decreasing the concentration of Na(2)SO(4) from 0.4 to 0 M destabilizes the intermediate and rate-limiting transition (TS2) states by 7 and 10 kJ mol(-)(1), respectively, and has little effect on the relative compactness of these states compared with that of the unfolded ensemble (beta(I) approximately 0.8, beta(TS2) approximately 0.9 in 0 to 0.4 M Na(2)SO(4)). Analysis of 10 variants of the protein in 0.2 M Na(2)SO(4) using Phi-values showed that the structural properties of the intermediate and TS2 are not altered significantly by the concentration of the kosmotrope. The data demonstrate that the rapid formation of a compact intermediate stabilized by non-native interactions during Im7* folding is not induced by high concentrations of the stabilizing salt, but is a generic feature of the folding of this protein.

The transition state for folding of a peripheral subunit-binding domain contains robust and ionic-strength dependent characteristics

The denaturant dependencies of the folding and unfolding kinetics were used to characterize the structure of the transition state for folding of E3BD, a peripheral subunit-binding domain. For the majority of E3BD mutants, the Phi-values calculated at 298 K from the analysis of chevron plots were in good agreement with those previously determined at 325 K using Arrhenius analysis. This agreement further demonstrates the general robustness of Phi-value analyses, since different experiments, methods of denaturation and thermodynamic assumptions were used to determine each set of Phi(F) values. The structure of the transition state for folding was grossly conserved at 298 K and 325 K, with residues in Helix I playing a lesser role in folding than those located in the 3(10) helix, disordered loop and Helix II. However, the energetic contributions of a cluster of basic residues close to the N-terminus and Helix I, which are an integral part of the ligand-binding site, were susceptible to ionic strength effects because of electrostatic strain in native and transition states of E3BD at low ionic strength. We found no evidence of the downhill folding previously proposed for E3BD, even though the conditions employed in this study significantly increased the energetic bias towards the native state.

Spectrin R16: broad energy barrier or sequential transition states?

A number of models have been proposed to account for nonlinearity in the relation between observed rate constants for folding and/or unfolding and denaturant concentration. Where curvature is seen principally in the arm of a chevron plot, three explanations are proposed: a change in the ground state at increasing concentration of urea, movement of the transition state along a broad energy barrier, and a switch between two sequential transition states separated by an on-pathway high-energy intermediate. Here we demonstrate that the latter two models in particular can be used to describe the data for the all-alpha protein spectrin R16. Further, whatever the method of analysis, the pattern of Phi-values seen is robust; thus we would draw the same conclusions from our data set independently of the method used for analysis. While this is not a novel observation, this is the first systematic study where a comparison has been made between Phi-values calculated using the broad and sequential transition state models.

Helix stability and hydrophobicity in the folding mechanism of the bacterial immunity protein Im9

Recent models suggest that the mechanism of protein folding is determined by the balance between the stability of secondary structural elements and the hydrophobicity of the sequence. Here we determine the role of these factors in the folding kinetics of Im9* by altering the secondary structure propensity or hydrophobicity of helices I, II or IV by the substitution of residues at solvent exposed sites. The folding kinetics of each variant were measured at pH 7.0 and 10 degrees C, under which conditions wild-type Im9* folds with two-state kinetics. We show that increasing the helicity of these sequences in regions known to be structured in the folding intermediate of Im7*, switches the folding of Im9* from a two- to three-state mechanism. By contrast, increasing the hydrophobicity of helices I or IV has no effect on the kinetic folding mechanism. Interestingly, however, increasing the hydrophobicity of solvent-exposed residues in helix II stabilizes the folding intermediate and the rate-limiting transition state, consistent with the view that this helix makes significant non-native interactions during folding. The results highlight the generic importance of intermediates in folding and show that such species can be populated by increasing helical propensity or by stabilizing inter-helix contacts through non-native interactions.

The folding of spectrin domains II: phi-value analysis of R16

Studies on the folding of helical proteins have shown a wide range of different mechanisms and highlighted the importance of helical propensity as a factor in determining folding mechanism. Here, we contribute to this interesting field with the protein engineering phi-value analysis of the 16th domain of chicken brain alpha-spectrin, R16. The fortuitous curvature seen in the unfolding arm of the chevron plot allows us to investigate both early and late events in folding. R16 is the first two-state helical protein for which this has been possible.

Equilibrium hydrogen exchange reveals extensive hydrogen bonded secondary structure in the on-pathway intermediate of Im7

The four-helical immunity protein Im7 folds through an on-pathway intermediate that has a specific, but partially misfolded, hydrophobic core. In order to gain further insight into the structure of this species, we have identified the backbone hydrogen bonds formed in the ensemble by measuring the amide exchange rates (under EX2 conditions) of the wild-type protein and a variant, I72V. In this mutant the intermediate is significantly destabilised relative to the unfolded state (deltadeltaG(ui) = 4.4 kJ/mol) but the native state is only slightly destabilised (deltadeltaG(nu) = 1.8 kJ/mol) at 10 degrees C in 2H2O, pH* 7.0 containing 0.4 M Na2SO4, consistent with the view that this residue forms significant non-native stabilising interactions in the intermediate state. Comparison of the hydrogen exchange rates of the two proteins, therefore, enables the state from which hydrogen exchange occurs to be identified. The data show that amides in helices I, II and IV in both proteins exchange slowly with a free energy similar to that associated with global unfolding, suggesting that these helices form highly protected hydrogen-bonded helical structure in the intermediate. By contrast, amides in helix III exchange rapidly in both proteins. Importantly, the rate of exchange of amides in helix III are slowed substantially in the Im7* variant, I72V, compared with the wild-type protein, whilst other amides exchange more rapidly in the mutant protein, in accord with the kinetics of folding/unfolding measured using chevron analysis. These data demonstrate, therefore, that local fluctuations do not dominate the exchange mechanism and confirm that helix III does not form stable secondary structure in the intermediate. By combining these results with previously obtained Phi-values, we show that the on-pathway folding intermediate of Im7 contains extensive, stable hydrogen-bonded structure in helices I, II and IV, and that this structure is stabilised by both native and non-native interactions involving amino acid side-chains in these helices.

Improved Gō-like Models Demonstrate the Robustness of Protein Folding Mechanisms Towards Non-native Interactions

The use of simple theoretical models has provided a considerable contribution to our present understanding of the means by which proteins adopt their native fold from the plethora of available unfolded states. A common assumption in building computationally tractable models has been the neglect of stabilizing non-native interactions in the class of models described as "Gō-like." The focus of this study is the characterization of the folding of a number of proteins via a Gō-like model, which aims to map a maximal amount of information reflecting the protein sequence onto a "minimalist" skeleton. This model is shown to contain sufficient information to reproduce the folding transition states of a number of proteins, including topologically analogous proteins that fold via different transition states. Remarkably, these models also demonstrate consistency with the general features of folding transition states thought to be stabilized by non-native interactions. This suggests that native interactions are the primary determinant of most protein folding transition states, and that non-native interactions lead only to local structural perturbations. A prediction is also included for an asymmetrical folding transition state of bacteriophage lambda protein W, which has yet to be subjected to experimental characterization.

Structural determinants of the rate of protein folding

To understand the mechanism of protein folding and to assist rational design of fast-folding, non-aggregating and stable artificial enzymes, it is essential to determine the structural parameters which govern the rate constants of folding, kf. It has been found that -logkf is a linear function of the so-called chain topology parameter (CTP) within the range of 10(-1)s(-1)< or = kf < or =10(8)s(-1). The correlation between -logkf and CTP is much improved than using previously published contact order (CO) method. It has been further suggested that short sequence separations may be preferred for the establishment of stable interactions for the design of novel artificial enzymes and the modification of slow-folding proteins with aggregating intermediates.

Equilibrium and kinetic folding of an alpha-helical Greek key protein domain: caspase recruitment domain (CARD) of RICK

We have characterized the equilibrium and kinetic folding of a unique protein domain, caspase recruitment domain (CARD), of the RIP-like interacting CLARP kinase (RICK) (RICK-CARD), which adopts a alpha-helical Greek key fold. At equilibrium, the folding of RICK-CARD is well described by a two-state mechanism representing the native and unfolded ensembles. The protein is marginally stable, with a DeltaG(H)()2(O) of 3.0 +/- 0.15 kcal/mol and an m-value of 1.27 +/- 0.06 kcal mol(-1) M(-1) (30 mM Tris-HCl, pH 8, 1 mM DTT, 25 degrees C). While the m-value is constant, the protein stability decreases in the presence of moderate salt concentrations (below 200 mM) and then increases at higher salt concentrations. The results suggest that electrostatic interactions are stabilizing in the native protein, and the favorable Coulombic interactions are reduced at low ionic strength. Above 200 mM salt, the results are consistent with Hofmeister effects. The unfolding pathway of RICK-CARD is complex and contains at least three non-native conformations. The refolding pathway of RICK-CARD also is complex, and the data suggest that the unfolded protein folds via two intermediate conformations prior to reaching the native state. Overall, the data suggest the presence of kinetically trapped, or misfolded, species that are on-pathway both in refolding and in unfolding.

Construction and characterization of a chimeric myoglobin

In order to investigate the functional and structural role of modular structure in globins, we have engineered a chimeric myoglobin (ChimMb) in which the first and third exon come from the gene coding for the sperm whale Mb and the second exon from the gene coding for Aplysia limacina Mb. This ChimMb, fused to the Maltose Binding Protein (MBP) and expressed in Escherichia coli as an apoprotein, binds protoheme in a 1:1 stoichiometric ratio. Based on some functional and spectroscopic properties, we conclude that the central core of the ChimMb (which derives from A. limacina) is native-like. On the other hand, the ChimMb deprived (by proteolytic digestion) of the fused MBP displays a considerably reduced stability. These results suggest that the sperm whale A-G-H nucleus does not contribute significantly to the overall stability of the ChimMb.

Structural analysis of the rate-limiting transition states in the folding of Im7 and Im9: similarities and differences in the folding of homologous proteins

The bacterial immunity proteins Im7 and Im9 fold with mechanisms of different kinetic complexity. Whilst Im9 folds in a two-state transition at pH 7.0 and 10 degrees C, Im7 populates an on-pathway intermediate under these conditions. In order to assess the role of sequence versus topology in the folding of these proteins, and to analyse the effect of populating an intermediate on the landscape for folding, we have determined the conformational properties of the rate-limiting transition state for Im9 folding/unfolding using Phi(F)-value analysis and have compared the results with similar data obtained previously for Im7. The data show that the rate-limiting transition states for Im9 and Im7 folding/unfolding are similar: both are compact (beta(T)=0.94 and 0.89, respectively) and contain three of the four native helices docked around a specific hydrophobic core. Significant differences are observed, however, in the magnitude of the Phi(F)-values obtained for the two proteins. Of the 20 residues studied in both proteins, ten have Phi(F)-values in Im7 that exceed those in Im9 by more than 0.2, and of these five differ by more than 0.4. The data suggest that the population of an intermediate in Im7 results in folding via a transition state ensemble that is conformationally restricted relative to that of Im9. The data are consistent with the view that topology is an important determinant of folding. Importantly, however, they also demonstrate that while the folding transition state may be conserved in homologous proteins that fold with two and three-state kinetics, the population of an intermediate can have a significant effect on the breadth of the transition state ensemble.

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.

Fast folding of a four-helical bundle protein

The FK506-FKBP12 binding-domain of the kinase FRAP (FRB) forms a classic up-down four-helical bundle. The folding pathway of this protein has been investigated using a combination of equilibrium and kinetic studies. The native state of the protein is stable with respect to the unfolded state by some 7 kcal mol(-1) at pH 6.0, 10 degrees C. A kinetic analysis of unfolding and refolding rate constants as a function of chemical denaturant concentration suggests that an intermediate state may be populated during folding at low concentrations of denaturant. The presence of this intermediate state is confirmed by refolding experiments performed in the presence of the hydrophobic dye 8-anilinonaphthalene-1 sulfonate (ANS). ANS binds to the partially folded intermediate state populated during the folding of FRB and undergoes a large change in fluorescence that can be detected using stopped-flow techniques. Analysis of the kinetic data suggests that the intermediate state is compact and it may even be a misfolded species that has to partially unfold before it can reach the transition state. Folding and unfolding rate constants in water are approximately 150-200 s(-1) and 0.005-0.06 s(-1), respectively, at neutral pH and 10 degrees C. The folding of FRB is somewhat slower than for other all-helical proteins, probably as a consequence of the formation of a metastable intermediate state. The folding rate constant in the absence of any populated intermediate can be estimated to be 8800 s(-1). Despite the presence of an intermediate state, which effectively slows folding, the protein still folds rapidly with a half-life of 5 ms at 10 degrees C. The dependence of the rate constants on denaturant concentration indicates that the transition state for folding is compact with some 80% of the surface area exposed in the unfolded state buried in the transition state. Data presented for FRB is compared with kinetic data obtained for other all-helical proteins.

Ultrafast folding of WW domains without structured aromatic clusters in the denatured state

Ultrafast-folding proteins are important for combining experiment and simulation to give complete descriptions of folding pathways. The WW domain family comprises small proteins with a three-stranded antiparallel beta-sheet topology. Previous studies on the 57-residue YAP 65 WW domain indicate the presence of residual structure in the chemically denatured state. Here we analyze three minimal core WW domains of 38-44 residues. There was little spectroscopic or thermodynamic evidence for residual structure in either their chemically or thermally denatured states. Folding and unfolding kinetics, studied by using rapid temperature-jump and continuous-flow techniques, show that each domain folds and unfolds very rapidly in a two-state transition through a highly compact transition state. Folding half-times were as short as 17 micros at 25 degrees C, within an order of magnitude of the predicted maximal rate of loop formation. The small size and topological simplicity of these domains, in conjunction with their very rapid two-state folding, may allow us to reduce the difference in time scale between experiment and theoretical simulation.

Acidic conditions stabilise intermediates populated during the folding of Im7 and Im9

The helical bacterial immunity proteins Im7 and Im9 have been shown to fold via kinetic mechanisms of differing complexity, despite having 60 % sequence identity. At pH 7.0 and 10 degrees C, Im7 folds in a three-state mechanism involving an on-pathway intermediate, while Im9 folds in an apparent two-state transition. In order to examine the folding mechanisms of these proteins in more detail, the folding kinetics of both Im7 and Im9 (at 10 degrees C in 0.4 M sodium sulphate) have been examined as a function of pH. Kinetic modelling of the folding and unfolding data for Im7 between pH 5.0 and 8.0 shows that the on-pathway intermediate is stabilised by more acidic conditions, whilst the native state is destabilised. The opposing effect of pH on the stability of these states results in a significant population of the intermediate at equilibrium at pH 6.0 and below. At pH 7.0, the folding and unfolding kinetics for Im9 can be fitted adequately by a two-state model, in accord with previous results. However, under acidic conditions there is a clear change of slope in the plot of the logarithm of the folding rate constant versus denaturant concentration, consistent with the population of one or more intermediate(s) early during folding. The kinetic data for Im9 at these pH values can be fitted to a three-state model, where the intermediate ensemble is stabilised and the native state destabilised as the pH is reduced, rationalising previous results that showed that an intermediate is not observed experimentally at pH 7.0. The data suggest that intermediate formation is a general step in immunity protein folding and demonstrate that it is necessary to explore a wide range of refolding conditions in order to show that intermediates do not form in the folding of other small, single-domain proteins.