Characterisation of the transition states for protein folding: towards a new level of mechanistic detail in protein engineering analysis

The Mechanism of Folding of Human Frataxin in Comparison to the Yeast Homologue - Broad Energy Barriers and the General Properties of the Transition State

The funneled energy landscape theory suggests that the folding pathway of homologous proteins should converge at the late stages of folding. In this respect, proteins displaying a broad energy landscape for folding are particularly instructive, allowing inferring both the early, intermediate and late stages of folding. In this paper we explore the folding mechanisms of human frataxin, an essential mitochondrial protein linked to the neurodegenerative disorder Friedreich's ataxia. Building upon previous studies on the yeast homologue, the folding pathway of human frataxin is thoroughly examined, revealing a mechanism implying the presence of a broad energy barrier, reminiscent of the yeast counterpart. Through an extensive site-directed mutagenesis, we employed a Φ -value analysis to map native-like contacts in the folding transition state. The presence of a broad energy barrier facilitated the exploration of such contacts in both early and late folding events. We compared results from yeast and human frataxin providing insights into the impact of native topology on the folding mechanism and elucidating the properties of the underlying free energy landscape. The findings are discussed in the context of the funneled energy landscape theory of protein folding.

Multistep molecular mechanism of amyloid-like aggregation of nucleic acid-binding domain of TDP-43

TDP-43 protein is associated with many neurodegenerative diseases and has been shown to adopt various oligomeric and fibrillar states. However, a detailed kinetic understanding of the structural transformation of the native form of the protein to the fibrillar state is missing. In this study, we delineate the temporal sequence of structural events during the amyloid-like assembly of the functional nucleic acid-binding domain of TDP-43. We kinetically mapped the aggregation process using multiple probes such as tryptophan and thioflavin T (ThT) fluorescence, circular dichroism (CD), and dynamic light scattering (DLS) targeting different structural events. Our data reveal that aggregation occurs in four distinct steps-very fast, fast, slow, and very slow. The "very fast" change results in partially unfolded forms that undergo conformational conversion, oligomerization and bind to ThT in the "fast step" to form higher order intermediates (HOI). The temporal sequence of the formation of ThT binding sites and conformational conversion depends upon the protein concentration. The HOI further undergoes structural rearrangement to form protofibrils in the "slow" step, which, consequently, assembles in the "very slow" step to form an amyloid-like assembly. The spectroscopic properties of the amyloid-like assembly across the protein concentration remain similar. Additionally, we observe no lag phase across protein concentration for all the probes studied, suggesting that the aggregation process follows a linear polymerization reaction. Overall, our study demonstrates that the amyloid-like assembly forms in multiple steps, which is also supported by the temperature dependence of the kinetics.

A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip

Gupta et al. use single-channel electrophysiology to investigate the gating mechanism of acetylcholine receptor ion channels. They propose that channel opening starts at the M2–M3 linker and ligand-binding sites and proceeds through four brief intermediate conformations before ending with the collapse of a gate bubble.

Nicotinic acetylcholine receptors are allosteric proteins that generate membrane currents by isomerizing (“gating”) between resting and active conformations under the influence of neurotransmitters. Here, to explore the mechanisms that link the transmitter-binding sites (TBSs) with the distant gate, we use mutant cycle analyses to measure coupling between residue pairs, phi value analyses to sequence domain rearrangements, and current simulations to reproduce a microsecond shut component (“flip”) apparent in single-channel recordings. Significant interactions between amino acids separated by >15 Å are rare; an exception is between the αM2–M3 linkers and the TBSs that are ∼30 Å apart. Linker residues also make significant, local interactions within and between subunits. Phi value analyses indicate that without agonists, the linker is the first region in the protein to reach the gating transition state. Together, the phi pattern and flip component suggest that a complete, resting↔active allosteric transition involves passage through four brief intermediate states, with brief shut events arising from sojourns in all or a subset. We derive energy landscapes for gating with and without agonists, and propose a structure-based model in which resting→active starts with spontaneous rearrangements of the M2–M3 linkers and TBSs. These conformational changes stabilize a twisted extracellular domain to promote transmembrane helix tilting, gate dilation, and the formation of a “bubble” that collapses to initiate ion conduction. The energy landscapes suggest that twisting is the most energetically unfavorable step in the resting→active conformational change and that the rate-limiting step in the reverse process is bubble formation.

Understanding the frustration arising from the competition between function, misfolding, and aggregation in a globular protein

Folding and function may impose different requirements on the amino acid sequences of proteins, thus potentially giving rise to conflict. Such a conflict, or frustration, can result in the formation of partially misfolded intermediates that can compromise folding and promote aggregation. We investigate this phenomenon by studying frataxin, a protein whose normal function is to facilitate the formation of iron-sulfur clusters but whose mutations are associated with Friedreich's ataxia. To characterize the folding pathway of this protein we carry out a Φ-value analysis and use the resulting structural information to determine the structure of the folding transition state, which we then validate by a second round of rationally designed mutagenesis. The analysis of the transition-state structure reveals that the regions involved in the folding process are highly aggregation-prone. By contrast, the regions that are functionally important are partially misfolded in the transition state but highly resistant to aggregation. Taken together, these results indicate that in frataxin the competition between folding and function creates the possibility of misfolding, and that to prevent aggregation the amino acid sequence of this protein is optimized to be highly resistant to aggregation in the regions involved in misfolding.

Unfolding simulations of holomyoglobin from four mammals: identification of intermediates and β-sheet formation from partially unfolded states

Myoglobin (Mb) is a centrally important, widely studied mammalian protein. While much work has investigated multi-step unfolding of apoMb using acid or denaturant, holomyoglobin unfolding is poorly understood despite its biological relevance. We present here the first systematic unfolding simulations of holoMb and the first comparative study of unfolding of protein orthologs from different species (sperm whale, pig, horse, and harbor seal). We also provide new interpretations of experimental mean molecular ellipticities of myoglobin intermediates, notably correcting for random coil and number of helices in intermediates. The simulated holoproteins at 310 K displayed structures and dynamics in agreement with crystal structures (R g ~1.48-1.51 nm, helicity ~75%). At 400 K, heme was not lost, but some helix loss was observed in pig and horse, suggesting that these helices are less stable in terrestrial species. At 500 K, heme was lost within 1.0-3.7 ns. All four proteins displayed exponentially decaying helix structure within 20 ns. The C- and F-helices were lost quickly in all cases. Heme delayed helix loss, and sperm whale myoglobin exhibited highest retention of heme and D/E helices. Persistence of conformation (RMSD), secondary structure, and ellipticity between 2-11 ns was interpreted as intermediates of holoMb unfolding in all four species. The intermediates resemble those of apoMb notably in A and H helices, but differ substantially in the D-, E- and F-helices, which interact with heme. The identified mechanisms cast light on the role of metal/cofactor in poorly understood holoMb unfolding. We also observed β-sheet formation of several myoglobins at 500 K as seen experimentally, occurring after disruption of helices to a partially unfolded, globally disordered state; heme reduced this tendency and sperm-whale did not display any sheet propensity during the simulations.

Analyses of the folding properties of ferredoxin-like fold proteins by means of a coarse-grained Gō model: relationship between the free energy profiles and folding cores

The folding mechanisms of proteins with multi-state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on-going research issues. In this study, we investigate ferredoxin-like fold proteins which have a simple topology and multi-state transitions. We analyze the folding processes by means of a coarse-grained Gō model. We are able to reproduce the differences in the folding mechanisms between U1A, which has a high-free-energy intermediate state, and ADA2h and S6, which fold into the native structure through two-state transitions. The folding pathways of U1A, ADA2h, S6, and the S6 circular permutant, S6_p54-55, are reproduced and compared with experimental observations. We show that the ferredoxin-like fold contains two common regions consisting folding cores as predicted in other studies and that U1A produces an intermediate state due to the distinct cooperative folding of each core. However, because one of the cores of S6 loses its cooperativity and the two cores of ADA2h are tightly coupled, these proteins fold into the native structure through a two-state mechanism.

Protein folding: adding a nucleus to guide helix docking reduces landscape roughness

The elongated three-helix‐bundle spectrin domains R16 and R17 fold and unfold unusually slowly over a rough energy landscape, in contrast to the homologue R15, which folds fast over a much smoother, more typical landscape. R15 folds via a nucleation–condensation mechanism that guides the docking of the A and C-helices. However, in R16 and R17, the secondary structure forms first and the two helices must then dock in the correct register. Here, we use variants of R16 and R17 to demonstrate that substitution of just five key residues is sufficient to alter the folding mechanism and reduce the landscape roughness. We suggest that, by providing access to an alternative, faster, folding route over their landscape, R16 and R17 can circumvent their slow, frustrated wild-type folding mechanism.

The H2A-H2B dimeric kinetic intermediate is stabilized by widespread hydrophobic burial with few fully native interactions

The H2A-H2B histone heterodimer folds via monomeric and dimeric kinetic intermediates. Within ∼5 ms, the H2A and H2B polypeptides associate in a nearly diffusion limited reaction to form a dimeric ensemble, denoted I₂ and I₂*, the latter being a subpopulation characterized by a higher content of nonnative structure (NNS). The I₂ ensemble folds to the native heterodimer, N₂, through an observable, first-order kinetic phase. To determine the regions of structure in the I₂ ensemble, we characterized 26 Ala mutants of buried hydrophobic residues, spanning the three helices of the canonical histone folds of H2A and H2B and the H2B C-terminal helix. All but one targeted residue contributed significantly to the stability of I₂, the transition state and N₂; however, only residues in the hydrophobic core of the dimer interface perturbed the I₂* population. Destabilization of I₂* correlated with slower folding rates, implying that NNS is not a kinetic trap but rather accelerates folding. The pattern of Φ values indicated that residues forming intramolecular interactions in the peripheral helices contributed similar stability to I₂ and N₂, but residues involved in intermolecular interactions in the hydrophobic core are only partially folded in I₂. These findings suggest a dimerize-then-rearrange model. Residues throughout the histone fold contribute to the stability of I₂, but after the rapid dimerization reaction, the hydrophobic core of the dimer interface has few fully native interactions. In the transition state leading to N₂, more native-like interactions are developed and nonnative interactions are rearranged.

Stable folding core in the folding transition state of an alpha-helical integral membrane protein

Defining the structural features of a transition state is important in understanding a folding reaction. Here, we use Φ-value and double mutant analyses to probe the folding transition state of the membrane protein bacteriorhodopsin. We focus on the final C-terminal helix, helix G, of this seven transmembrane helical protein. Φ-values could be derived for 12 amino acid residues in helix G, most of which have low or intermediate values, suggesting that native structure is disrupted at these amino acid positions in the transition state. Notably, a cluster of residues between E204 and M209 all have Φ-values close to zero. Disruption of helix G is further confirmed by a low Φ-value of 0.2 between residues T170 on helix F and S226 on helix G, suggesting the absence of a native hydrogen bond between helices F and G. Φ-values for paired mutations involved in four interhelical hydrogen bonds revealed that all but one of these bonds is absent in the transition state. The unstructured helix G contrasts with Φ-values along helix B that are generally high, implying native structure in helix B in the transition state. Thus helix B seems to constitute part of a stable folding nucleus while the consolidation of helix G is a relatively late folding event. Polarization of secondary structure correlates with sequence position, with a structured helix B near the N terminus contrasting with an unstructured C-terminal helix G.

Cutting off functional loops from homodimeric enzyme superoxide dismutase 1 (SOD1) leaves monomeric β-barrels

Demetallation of the homodimeric enzyme Cu/Zn-superoxide dismutase (SOD1) is known to unleash pronounced dynamic motions in the long active-site loops that comprise almost a third of the folded structure. The resulting apo species, which shows increased propensity to aggregate, stands out as the prime disease precursor in amyotrophic lateral sclerosis (ALS). Even so, the detailed structural properties of the apoSOD1 framework have remained elusive and controversial. In this study, we examine the structural interplay between the central apoSOD1 barrel and the active-site loops by simply cutting them off; loops IV and VII were substituted with short Gly-Ala-Gly linkers. The results show that loop removal breaks the dimer interface and leads to soluble, monomeric β-barrels with high structural integrity. NMR-detected nuclear Overhauser effects are found between all of the constituent β-strands, confirming ordered interactions across the whole barrel. Moreover, the breathing motions of the SOD1 barrel are overall insensitive to loop removal and yield hydrogen/deuterium protection factors typical for cooperatively folded proteins (i.e. the active-site loops act as a "bolt-on" domain with little dynamic influence on its structural foundation). The sole exceptions are the relatively low protection factors in β-strand 5 and the turn around Gly-93, a hot spot for ALS-provoking mutations, which decrease even further upon loop removal. Taken together, these data suggest that the cytotoxic function of apoSOD1 does not emerge from its folded ground state but from a high energy intermediate or even from the denatured ensemble.

Folding kinetics of the cooperatively folded subdomain of the IκBα ankyrin repeat domain

The ankyrin repeat (AR) domain of IκBα consists of a cooperative folding unit of roughly four ARs (AR1-AR4) and of two weakly folded repeats (AR5 and AR6). The kinetic folding mechanism of the cooperative subdomain, IκBα(67-206), was analyzed using rapid mixing techniques. Despite its apparent architectural simplicity, IκBα(67-206) displays complex folding kinetics, with two sequential on-pathway high-energy intermediates. The effect of mutations to or away from the consensus sequences of ARs on folding behavior was analyzed, particularly the GXTPLHLA motif, which have not been examined in detail previously. Mutations toward the consensus generally resulted in an increase in folding stability, whereas mutations away from the consensus resulted in decreased overall stability. We determined the free energy change upon mutation for three sequential transition state ensembles along the folding route for 16 mutants. We show that folding initiates with the formation of the interface of the outer helices of AR3 and AR4, and then proceeds to consolidate structure in these repeats. Subsequently, AR1 and AR2 fold in a concerted way in a single kinetic step. We show that this mechanism is robust to the presence of AR5 and AR6 as they do not strongly affect the folding kinetics. Overall, the protein appears to fold on a rather smooth energy landscape, where the folding mechanism conforms a one-dimensional approximation. However, we note that the AR does not necessarily act as a single folding element.

Mutational studies uncover non-native structure in the dimeric kinetic intermediate of the H2A-H2B heterodimer

The folding pathway of the histone H2A-H2B heterodimer minimally includes an on-pathway, dimeric, burst-phase intermediate, I(2). The partially folded H2A and H2B monomers populated at equilibrium were characterized as potential monomeric kinetic intermediates. Folding kinetics were compared for initiation from isolated, folded monomers and the heterodimer unfolded in 4 M urea. The observed rates were virtually identical above 0.4 M urea, exhibiting a log-linear relationship on the final denaturant concentration. Below approximately 0.4 M urea (concentrations inaccessible from the 4-M urea unfolded state), a rollover in the rates was observed; this suggests that a component of the I(2) ensemble contains non-native structure that rearranges/isomerizes to a more native-like species. The contribution of helix propensity to the stability of the I(2) ensemble was assessed with a set of H2A-H2B mutants containing Ala and Gly replacements at nine sites, focusing mainly on the long, central alpha2 helix. Equilibrium and kinetic folding/unfolding data were collected to determine the effects of the mutations on the stability of I(2) and the transition state between I(2) and N(2). This limited mutational study indicated that residues in the alpha2 helices of H2A and H2B as well as alpha1 of H2B and both the C-terminus of alpha3 and the short alphaC helix of H2A contribute to the stability of the I(2) burst-phase species. Interestingly, at least eight of the nine targeted residues stabilize I(2) by interactions that are non-native to some extent. Given that destabilizing I(2) and these non-native interactions does not accelerate folding, it is concluded that the native and non-native structures present in the I(2) ensemble enable efficient folding of H2A-H2B.

Is protein folding rate dependent on number of folding stages? Modeling of protein folding with ferredoxin-like fold

Statistical analysis of protein folding rates has been done for 84 proteins with available experimental data. A surprising result is that the proteins with multi-state kinetics from the size range of 50-100 amino acid residues (a.a.) fold as fast as proteins with two-state kinetics from the same size range. At the same time, the proteins with two-state kinetics from the size range 101-151 a.a. fold faster than those from the size range 50-100 a.a. Moreover, it turns out unexpectedly that usually in the group of structural homologs from the size range 50-100 a.a., proteins with multi-state kinetics fold faster than those with two-state kinetics. The protein folding for six proteins with a ferredoxin-like fold and with a similar size has been modeled using Monte Carlo simulations and dynamic programming. Good correlation between experimental folding rates, some structural parameters, and the number of Monte Carlo steps has been obtained. It is shown that a protein with multi-state kinetics actually folds three times faster than its structural homologs.

Experimental evidence for a frustrated energy landscape in a three-helix-bundle protein family

Energy landscape theory is a powerful tool for understanding the structure and dynamics of complex molecular systems, in particular biological macromolecules. 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 the protein's native structure. Theory has shown that roughness in the energy landscape will lead to slower folding, but derivation of detailed experimental descriptions of this landscape is challenging. Simple folding models show that folding is significantly influenced by chain entropy; proteins in which the contacts are local fold quickly, owing to the low entropy cost of forming stabilizing, native contacts during folding. For some protein families, stability is also a determinant of folding rate constants. 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 alpha-spectrin. R15 folds approximately 3,000 times faster than its homologues, although they have similar structures, stabilities and, as far as can be determined, transition-state stabilities. Here we show that landscape roughness (internal friction) is responsible for the slower folding and unfolding of R16 and R17. We use chimaeric 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 misdocking of the long helices during folding. Theoretical studies have suggested that rugged landscapes will result in slower folding; here we show experimentally that such a phenomenon directly influences the folding kinetics of a 'normal' protein, that is, one with a significant energy barrier that folds on a relatively slow, millisecond-second, timescale.

Distinguishing between smooth and rough free energy barriers in protein folding

Analysis of curved chevron plots is a powerful tool in investigating protein folding pathways, as the curvatures can be used to gain information about both early and late folding events. When and if accumulation of low-energy intermediates can be ruled out, two different models have classically been applied to describe curved chevron plots, namely , (i) Hammond effects along smooth barrier profiles and (ii) changes in the rate-limiting step between two discrete transition states. The two models lead to very similar numerical solutions, which are generally indistinguishable. This is not surprising, since the smooth barrier assumption approximates barrier profiles with a more complex topology involving multiple local maxima that are too close, or too broad, to yield clear-cut kinks in the chevron data. In this work, we have reconstructed the transition state shifts as a function of protein stability over a wide stability range for three small globular proteins, to screen for fingerprints more sensitive for different barrier profiles. We show that such an analysis represents a valuable test for the discrimination between the two different scenarios.

Topology is the principal determinant in the folding of a complex all-alpha Greek key death domain from human FADD

In order to elucidate the relative importance of secondary structure and topology in determining folding mechanism, we have carried out a phi-value analysis of the death domain (DD) from human FADD. FADD DD is a 100 amino acid domain consisting of six anti-parallel alpha helices arranged in a Greek key structure. We asked how does the folding of this domain compare with that of (a) other all-alpha-helical proteins and (b) other Greek key proteins? Is the folding pathway determined mainly by secondary structure or is topology the principal determinant? Our Φ-value analysis reveals a striking resemblance to the all-beta Greek key immunoglobulin-like domains. Both fold via diffuse transition states and, importantly, long-range interactions between the four central elements of secondary structure are established in the transition state. The elements of secondary structure that are less tightly associated with the central core are less well packed in both cases. Topology appears to be the dominant factor in determining the pathway of folding in all Greek key domains.

The transition state for integral membrane protein folding

Biology relies on the precise self-assembly of its molecular components. Generic principles of protein folding have emerged from extensive studies on small, water-soluble proteins, but it is unclear how these ideas are translated into more complex situations. In particular, the one-third of cellular proteins that reside in biological membranes will not fold like water-soluble proteins because membrane proteins need to expose, not hide, their hydrophobic surfaces. Here, we apply the powerful protein engineering method of Phi-value analysis to investigate the folding transition state of the alpha-helical membrane protein, bacteriorhodopsin, from a partially unfolded state. Our results imply that much of helix B of the seven-transmembrane helical protein is structured in the transition state with single-point alanine mutations in helix B giving Phi values >0.8. However, residues Y43 and T46 give lower Phi values of 0.3 and 0.5, respectively, suggesting a possible reduction in native structure in this region of the helix. Destabilizing mutations also increase the activation energy of folding, which is accompanied by an apparent movement of the transition state toward the partially unfolded state. This apparent transition state movement is most likely due to destabilization of the structured, unfolded state. These results contrast with the Hammond effect seen for several water-soluble proteins in which destabilizing mutations cause the transition state to move toward, and become closer in energy to, the folded state. We thus introduce a classic folding analysis method to membrane proteins, providing critical insight into the folding transition state.

Proline 54 trans-cis isomerization is responsible for the kinetic partitioning at the last-step photocycle of photoactive yellow protein

Photoactive yellow protein (PYP), a blue-light photoreceptor for Ectothiorhodospira halophila, has provided a unique system for studying protein folding that is coupled with a photocycle. Upon receptor activation by blue light, PYP proceeds through a photocycle that includes a partially folded signaling state. The last-step photocycle is a thermal recovery reaction from the signaling state to the native state. Bi-exponential kinetics had been observed for the last-step photocycle; however, the slow phase of the bi-exponential kinetics has not been extensively studied. Here we analyzed both fast and slow phases of the last-step photocycle in PYP. From the analysis of the denaturant dependence of the fast and slow phases, we found that the last-step photocycle proceeds through parallel channels of the folding pathway. The burial of the solvent-accessible area was responsible for the transition state of the fast phase, while structural rearrangement from the compact state to the native state was responsible for the transition state of the slow phase. The photocycle of PYP was linked to the thermodynamic cycle that includes both unfolding and refolding of the fast- and slow-phase intermediates. In order to test the hypothesis of proline-limited folding for the slow phase, we constructed two proline mutants: P54A and P68A. We found that only a single phase of the last-step photocycle was observed in P54A. This suggests that there is a low energy barrier between trans to cis conformation in P54 in the light-induced state of PYP, and the resulting cis conformation of P54 generates a slow-phase kinetic trap during the photocycle-coupled folding pathway of PYP.

The folding process of acylphosphatase from Escherichia coli is remarkably accelerated by the presence of a disulfide bond

The acylphosphatase from Escherichia coli (EcoAcP) is the first AcP so far studied with a disulfide bond. A mutational variant of the enzyme lacking the disulfide bond has been produced by substituting the two cysteine residues with alanine (EcoAcP mutational variant C5A/C49A, mutEcoAcP). The native states of the two protein variants are similar, as shown by far-UV and near-UV circular dichroism and dynamic light-scattering measurements. From unfolding experiments at equilibrium using intrinsic fluorescence and far-UV circular dichroism as probes, EcoAcP shows an increased conformational stability as compared with mutEcoAcP. The wild-type protein folds according to a two-state model with a very fast rate constant (k(F)(H2O)=72,600 s(-1)), while mutEcoAcP folds ca 1500-fold slower, via the accumulation of a partially folded species. The correlation between the hydrophobicity of the polypeptide chain and the folding rate, found previously in the AcP-like structural family, is maintained only when considering the mutant but not the wild-type protein, which folds much faster than expected from this correlation. Similarly, the correlation between the relative contact order and the folding rate holds only for mutEcoAcP. The correlation also holds for EcoAcP, provided the relative contact order value is recalculated by considering the disulfide bridge as an alternate path for the backbone to determine the shortest sequence separation between contacting residues. These results indicate that the presence of a disulfide bond in a protein is an important determinant of the folding rate and allows its contribution to be determined in quantitative terms.

Temperature-dependent Hammond behavior in a protein-folding reaction: analysis of transition-state movement and ground-state effects

Characterization of the transition-state ensemble and the nature of the free-energy barrier for protein folding are areas of intense activity and some controversy. A key issue that has emerged in recent years is the width of the free-energy barrier and the susceptibility of the transition state to movement. Here we report denaturant-induced and temperature-dependent folding studies of a small mixed alpha-beta protein, the N-terminal domain of L9 (NTL9). The folding of NTL9 was determined using fluorescence-detected stopped-flow fluorescence measurements conducted at seven different temperatures between 11 and 40 degrees C. Plots of the log of the observed first-order rate constant versus denaturant concentration, "chevron plots," displayed the characteristic V shape expected for two-state folding. There was no hint of deviation from linearity even at the lowest denaturant concentrations. The relative position of the transition state, as judged by the Tanford beta parameter, beta(T), shifts towards the native state as the temperature is increased. Analysis of the temperature dependence of the kinetic and equilibrium m values indicates that the effect is due to significant movement of the transition state and also includes a contribution from temperature-dependent ground-state effects. Analysis of the Leffler plots, plots of Delta G versus Delta G degrees, and their cross-interaction parameters confirms the transition-state movement. Since the protein is destabilized at high temperature, the shift represents a temperature-dependent Hammond effect. This provides independent confirmation of a recent theoretical prediction. The magnitude of the temperature-denaturant cross-interaction parameter is larger for NTL9 than has been reported for the few other cases studied. The implications for temperature-dependent studies of protein folding are discussed.

Simulating temperature jumps for protein folding

We present a new computational methodology aimed to calculate the thermodynamics and kinetics of peptide folding. We focus in particular on temperature jump experiments of folding rates and show how a combination of replica exchange molecular dynamics (REMD) followed by multiplexed molecular dynamics starting from structures taken from the REMD runs can be used to extract properties in line with experiments. A model system, alanine20, is studied in this article as a proof of principle and used to describe the methodology.

Crosstalk between the protein surface and hydrophobic core in a core-swapped fibronectin type III domain

Two homologous fibronectin type III (fnIII) domains, FNfn10 (the 10th fnIII domain of human fibronectin) and TNfn3 (the third fnIII domain of human tenascin), have essentially the same backbone structure, although they share only ∼ 24% sequence identity. While they share a similar folding mechanism with a common core of key residues in the folding transition state, they differ in many other physical properties. We use a chimeric protein, FNoTNc, to investigate the molecular basis for these differences. FNoTNc is a core-swapped protein, containing the “outside” (surface and loops) of FNfn10 and the hydrophobic core of TNfn3. Remarkably, FNoTNc retains the structure of the parent proteins despite the extent of redesign, allowing us to gain insight into which components of each parent protein are responsible for different aspects of its behaviour. Naively, one would expect properties that appear to depend principally on the core to be similar to TNfn3, for example, the response to mutations, folding kinetics and side-chain dynamics, while properties apparently determined by differences in the surface and loops, such as backbone dynamics, would be more like FNfn10. While this is broadly true, it is clear that there are also unexpected crosstalk effects between the core and the surface. For example, the anomalous response of FNfn10 to mutation is not solely a property of the core as we had previously suggested.

Folding versus aggregation: polypeptide conformations on competing pathways

Protein aggregation has now become recognised as an important and generic aspect of protein energy landscapes. Since the discovery that numerous human diseases are caused by protein aggregation, the biophysical characterisation of misfolded states and their aggregation mechanisms has received increased attention. Utilising experimental techniques and computational approaches established for the analysis of protein folding reactions has ensured rapid advances in the study of pathways leading to amyloid fibrils and amyloid-related aggregates. Here we describe recent experimental and theoretical advances in the elucidation of the conformational properties of dynamic, heterogeneous and/or insoluble protein ensembles populated on complex, multidimensional protein energy landscapes. We discuss current understanding of aggregation mechanisms in this context and describe how the synergy between biochemical, biophysical and cell-biological experiments are beginning to provide detailed insights into the partitioning of non-native species between protein folding and aggregation pathways.

Stable intermediate states and high energy barriers in the unfolding of GFP

We present a study of the denaturation of a truncated, cycle3 variant of green fluorescent protein (GFP). Chemical denaturation is used to unfold the protein, with changes in structure being monitored by the green fluorescence, tyrosine fluorescence and far-UV circular dichroism. The results show that the denaturation behaviour of GFP is complex compared to many small proteins: equilibrium is established only very slowly, over the time course of weeks, suggesting that there are high folding/unfolding energy barriers. Unfolding kinetics confirm that the rates of unfolding at low concentrations of denaturant are very low, consistent with the slow establishment of the equilibrium. In addition, we find that GFP significantly populates an intermediate state under equilibrium conditions, which is compact and stable with respect to the unfolded state (m(IU)=4.6 kcal mol(-1) M(-1) and Delta G(IU)=12.5 kcal mol(-1)). The global and local stability of GFP was probed further by measuring the hydrogen/deuterium (H/D) NMR exchange rates of more than 157 assigned amide protons. Analysis at two different values of pH showed that amide protons within the beta-barrel structure exchange at the EX2 limit, consequently, free energies of exchange could be calculated and compared to those obtained from the denaturation-curve studies providing further support for the three-state model and the existence of a stable intermediate state. Analysis reveals that amide protons in beta-strands 7, 8, 9 and 10 have, on average, higher exchange rates than others in the beta-barrel, suggesting that there is greater flexibility in this region of the protein. Forty or so amide protons were found which do not undergo significant exchange even after several months and these are clustered into a core region encompassing most of the beta-strands, at least at one end of the barrel structure. It is likely that these residues play an important role in stabilizing the structure of the intermediate state. The intermediate state observed in the chemical denaturation studies described here, is similar to that observed at pH 4 in other studies.

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.

Malleability of protein folding pathways: a simple reason for complex behaviour

Although the structures of native proteins are generally unique, the pathways by which they form are often free to vary. Some proteins fold by a multitude of different pathways, whereas others seem restricted to only one choice. An explanation for this variation in folding behaviour has recently emerged from studies of transition state changes: the number of accessible pathways is linked to the number of nucleation motifs contained within the native topology. We refer to these nucleation motifs as 'foldons', as they approach the size of an independent cooperative unit. Thus, with respect to pathway malleability and the composition of the folding funnel, proteins can be seen as modular assemblies of competing foldons. For the split beta-alpha-beta fold, these foldons are two-strand-helix motifs coupled by spatial overlap.

Mutational Analysis of the Folding Transition State of the C-Terminal Domain of Ribosomal Protein L9: A Protein with an Unusual β-Sheet Topology

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

Biophysical Characterisation of the Small Ankyrin Repeat Protein Myotrophin

The 118 residue protein myotrophin is composed of four ankyrin repeats that stack linearly to form an elongated, predominantly alpha-helical structure. The protein folds via a two-state mechanism at equilibrium. The free energy change of unfolding in water (DeltaG(U-N)(H(2)O)) is 5.8 kcal.mol(-1). The chevron plot reveals that the folding reaction has a broad energy barrier and that it conforms to a two-state mechanism. The rate of folding in water (k(f)(H(2)O)) of 95 s(-1) is several orders of magnitude slower than the value predicted by topological calculations. Proline mutants were used to show that the minor kinetic phases observed for myotrophin arise from heterogeneity of the ground states due to cis-trans isomerisation of prolyl as well as non-prolyl peptide bonds. Myotrophin is the first example of a naturally occurring ankyrin repeat protein that conforms to an apparent two-state mechanism at equilibrium and under kinetic conditions, making it highly suitable for high resolution protein folding studies.

Folding of S6 structures with divergent amino acid composition: pathway flexibility within partly overlapping foldons

Studies of circular permutants have demonstrated that the folding reaction of S6 from Thermus thermophilus (S6(T)) is malleable and responds in an ordered manner to changes of the sequence separation between interacting residues: the S6(T) permutants retain a common nucleation pattern in the form of a two-strand-helix motif that can be recruited from different parts of the structure. To further test the robustness of the two-strand-helix nucleus we have here determined the crystal structure and folding reaction of an evolutionary divergent S6 protein from the hyperthermophilic bacterium Aquifex aeolicus (S6(A)). Although the overall topology of S6(A) is very similar to that of S6(T) the architecture of the hydrophobic core is radically different by containing a large proportion of stacked Phe side-chains. Despite this disparate core composition, the folding rate constant and the kinetic m values of S6(A) are identical to those of S6(T). The folding nucleus of S6(A) is also found to retain the characteristic two-strand-helix motif of the S6(T) permutants, but with a new structural emphasis. The results suggest that the protein folding reaction is linked to topology only in the sense that the native-state topology determines the repertoire of accessible nucleation motifs. If the native structure allows several equivalent ways of recruiting a productive nucleus the folding reaction is free to redistribute within these topological constraints.

Phi-value analysis of apo-azurin folding: comparison between experiment and theory

Pseudomonas aeruginosa azurin is a 128-residue beta-sandwich metalloprotein; in vitro kinetic experiments have shown that it folds in a two-state reaction. Here, we used a variational free energy functional to calculate the characteristics of the transition state ensemble (TSE) for folding of the apo-form of P. aeruginosa azurin and investigate how it responds to thermal and mutational changes. The variational method directly yields predicted chevron plots for wild-type and mutant apo-forms of azurin. In parallel, we performed in vitro kinetic-folding experiments on the same set of azurin variants using chemical perturbation. Like the wild-type protein, all apo-variants fold in apparent two-state reactions both in calculations and in stopped-flow mixing experiments. Comparisons of phi (phi) values determined from the experimental and theoretical chevron parameters reveal an excellent agreement for most positions, indicating a polarized, highly structured TSE for folding of P. aeruginosa apo-azurin. We also demonstrate that careful analysis of side-chain interactions is necessary for appropriate theoretical description of core mutants.

Folding of Cu/Zn superoxide dismutase suggests structural hotspots for gain of neurotoxic function in ALS: parallels to precursors in amyloid disease

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease linked to misfolding of the ubiquitous enzyme Cu/Zn superoxide dismutase (SOD). In contrast to other protein-misfolding disorders with similar neuropathogenesis, ALS is not always associated with the in vivo deposition of protein aggregates. Thus, under the assumption that all protein-misfolding disorders share at primary level a similar disease mechanism, ALS constitutes an interesting disease model for identifying the yet-mysterious precursor states from which the cytotoxic pathway emerges. In this study, we have mapped out the conformational repertoire of the apoSOD monomer through analysis of its folding behavior. The results allow us to target the regions of the SOD structure that are most susceptible to unfolding locally under physiological conditions, leading to the exposure of structurally promiscuous interfaces that are normally hidden in the protein's interior. The structure of this putative ALS precursor is strikingly similar to those implicated in amyloid disease.

Effects of surface tethering on protein folding mechanisms

The folding mechanisms of proteins are increasingly being probed through single-molecule experiments in which the protein is immobilized on a surface. Nevertheless, a clear understanding of how the surface might affect folding, and whether or not it changes folding from its bulk behavior, is lacking. In this work, we use molecular dynamics simulations of a model beta-barrel protein tethered to a surface to systematically investigate how the surface impacts folding. In the bulk, this protein folds in a three-state manner through a compact intermediate state, and its transition state (TS) has a well formed hydrophobic core. Upon tethering, we find that folding rates and stability are impacted differently by the surface, with dependencies on both the length and location of the tether. Significant changes in folding times are observed for tether points that do not alter the folding temperature. Tethering also locally enhances the formation of structure for residues proximal to the tether point. We find that neither the folding mechanism nor the TS of this protein are altered if the tether is in a fully structured or completely unstructured region of the TS. By contrast, tethering in a partially structured region of the TS leads to dramatic changes. For one such tether point, the intermediate present in bulk folding is eliminated, leading to a two-state folding process with a heterogeneous, highly unstructured TS ensemble. These results have implications for both the design of single-molecule experiments and biotechnological applications of tethered proteins.

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.

Identification of the minimal protein-folding nucleus through loop-entropy perturbations

To explore the plasticity and structural constraints of the protein-folding nucleus we have constructed through circular permutation four topological variants of the ribosomal protein S6. In effect, these topological variants represent entropy mutants with maintained spatial contacts. The proteins were characterized at two complementary levels of detail: by phi-value analysis estimating the extent of contact formation in the transition-state ensemble and by Hammond analysis measuring the site-specific growth of the folding nucleus. The results show that, although the loop-entropy alterations markedly influence the appearance and structural location of the folding nucleus, it retains a common motif of one helix docking against two strands. This nucleation motif is built around a shared subset of side chains in the center of the hydrophobic core but extends in different directions of the S6 structure following the permutant-specific differences in local loop entropies. The adjustment of the critical folding nucleus to alterations in loop entropies is reflected by a direct correlation between the phi-value change and the accompanying change in local sequence separation.

Unveiling a hidden folding intermediate in c-type cytochromes by protein engineering

Several investigators have highlighted a correlation between the basic features of the folding process of a protein and its topology, which dictates the folding pathway. Within this conceptual framework we proposed that different members of the cytochrome c (cyt c) family share the same folding mechanism, involving a consensus partially structured state. Pseudomonas aeruginosa cyt c(551) (Pa cyt c(551)) folds via an apparent two-state mechanism through a high energy intermediate. Here we present kinetic evidence demonstrating that it is possible to switch its folding mechanism from two to three state, stabilizing the high energy intermediate by rational mutagenesis. Characterization of the folding kinetics of one single-site mutant of the Pa cyt c(551) (Phe(7) to Ala) indeed reveals an additional refolding phase and a fast unfolding process which are explained by the accumulation of a partially folded species. Further kinetic analysis highlights the presence of two parallel processes both leading to the native state, suggesting that the above mentioned species is a non obligatory on-pathway intermediate. Determination of the crystallographic structure of F7A shows the presence of an extended internal cavity, which hosts three "bound" water molecules and a H-bond in the N-terminal helix, which is shorter than in the wild type protein. These two features allow us to propose a detailed structural interpretation for the stabilization of the native and especially the intermediate states induced by a single crucial mutation. These results show how protein engineering, x-ray crystallography and state-of-the-art kinetics concur to unveil a folding intermediate and the structural determinants of its stability.

An Experimental Investigation of Conformational Fluctuations in Proteins G and L

The B1 domains of streptococcal proteins G and L are structurally similar, but they have different sequences and they fold differently. We have measured their NMR spectra at variable temperature using a range of concentrations of denaturant. Many residues have curved amide proton temperature dependence, indicating that they significantly populate alternative, locally unfolded conformations. The results, therefore, provide a view of the locations of low-lying, locally unfolded conformations. They indicate approximately 4-6 local minima for each protein, all within ca. 2.5 kcal/mol of the native state, implying a locally rough energy landscape. Comparison with folding data for these proteins shows that folding involves most molecules traversing a similar path, once a transition state containing a beta hairpin has been formed, thereby defining a well-populated pathway down the folding funnel. The hairpin that directs the folding pathway differs for the two proteins and remains the most stable part of the folded protein.

Symmetric connectivity of secondary structure elements enhances the diversity of folding pathways

The influence of native connectivity of secondary structure elements (SSE) on folding is studied using coarse-grained models of proteins with mixed alpha and beta structure and the analysis of the structural database of wild-type proteins. We found that the distribution of SSE along a sequence determines the diversity of folding pathways. If alpha and beta SSE are localized in different parts of a sequence, the diversity of folding pathways is restricted. An even (symmetric) distribution of alpha and beta SSE with respect to sequence midpoint favors multiple folding routes. Simulations are supplemented by the database analysis of the distribution of SSE in wild-type protein sequences. On an average, two-thirds of wild-type proteins with mixed alpha and beta structure have symmetric distribution of alpha and beta SSE. The propensity for symmetric distribution of SSE is especially evident for large proteins with the number of SSE > or = 10. We suggest that symmetric SSE distribution in protein sequences may arise due to nearly random allocation of alpha and beta structure along wild-type sequences. The tendency of long sequences to misfold is perhaps compensated by the enhanced pathway diversity. In addition, folding pathways are shown to progress via hierarchic assembly of SSE in accordance with their proximity along a sequence. We demonstrate that under mild denaturation conditions folding and unfolding pathways are similar. However, the reversibility of folding/unfolding pathways is shown to depend on the distribution of SSE. If alpha and beta SSE are localized in different parts of a sequence, folding and unfolding pathways are likely to coincide.

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.

Stabilization of the ribosomal protein S6 by trehalose is counterbalanced by the formation of a putative off-pathway species

The effect of trehalose on folding and stability of the small ribosomal protein S6 was studied. Non-disruptive point mutations distributed along the protein structure were analyzed to characterize the stabilizing effect of trehalose and map the folding pathway of S6. On average, the stability of the wild-type and S6 mutants increases by 3 kcal/mol M trehalose. Despite the non-specific thermodynamic stabilization mechanism, trehalose particularly stabilizes the less destabilized mutants. Folding/unfolding kinetics shows clearly that trehalose induces the collapse of the unfolded state to an off-pathway intermediate with non-native diffuse contacts. This state is similar to the collapsed state induced by high concentrations of stabilizing salts, as previously reported. Although it leads to the accumulation of this off-pathway intermediate, trehalose does not change the compactness of the transition state ensemble. Furthermore, the productive folding pathway of S6 is not affected by trehalose as shown by a Phi-value analysis. The unfolded state ensemble of S6 should be more compact in the presence of trehalose and therefore destabilized due to decreased conformational entropy. Increased compaction of the unfolded state ensemble might also occur for more stable mutants of S6, thus explaining the synergistic effect of trehalose and point mutations on protein stabilization.

NMR Spectroscopic Characterization of Millisecond Protein Folding by Transverse Relaxation Dispersion Measurements

The cold shock protein CspB adopts its native and functional tertiary structure on the millisecond time scale. We employed transverse relaxation NMR methods, which allow a quantitative measurement of the cooperativity of this fast folding reaction on a residue basis. Thereby, chemical exchange contributions to the transverse relaxation rate (R(2)) were observed for every residue of CspB verifying the potential of this method to identify not only local dynamics but also to characterize global events. Toward this end, the homogeneity of the transition state of folding was probed by comparing Chevron plots (i.e., dependence of the apparent folding rate on the denaturant concentration) determined by stopped-flow fluorescence with Chevron plots of six residues acquired by R(2) dispersion experiments. The coinciding results obtained for probes at different locations in the three-dimensional structure of CspB indicate the ability and significance of transverse relaxation NMR to determine Chevron plots on a residue-by-residue basis providing detailed insights on the nature of the transition state of folding.

Snapshots of a Dynamic Folding Nucleus in Zinc-Substituted Pseudomonas aeruginosa Azurin

Zinc-substituted Pseudomonas aeruginosa azurin folds in two-state equilibrium and kinetic reactions. In the unfolded state, the zinc ion remains bound to the unfolded polypeptide via two native-state ligands (His117 and Cys112). The significantly curved Chevron plot for zinc-substituted azurin was earlier ascribed to movement of the folding-transition state. At low concentrations of denaturant, the transition state occurs early in the folding reaction (low Tanford beta-value), whereas at high-denaturant concentration, it moves closer to the native structure (high Tanford beta-value). Here, we use this movement to track the formation and growth of zinc-substituted azurin's folding nucleus with atomic resolution using protein engineering. The average phi (phi) value for 17 positions (covering all secondary-structure elements) goes from 0.25 in 0 M GuHCl (beta approximately 0.46) to 0.76 in 4 M GuHCl (beta approximately 0.86); a phi-value of 1 or 0 indicates native-like or unfolded-like interactions, respectively. Analysis of individual phi-values reveals a delocalized nucleus where structure condenses around a leading density centered on Leu50 in the core. The diffuse moving transition state for zinc-substituted azurin is in sharp contrast to the fixed polarized folding nucleus observed for apo-azurin. The dramatic difference in apparent kinetic behavior for the two forms of azurin can be rationalized as a minor alteration on a common free-energy profile that exhibits a broad activation barrier.

Reconstruction of the src-SH3 Protein Domain Transition State Ensemble using Multiscale Molecular Dynamics Simulations

We use an integrated computational approach to reconstruct accurately the transition state ensemble (TSE) for folding of the src-SH3 protein domain. We first identify putative TSE conformations from free energy surfaces generated by importance sampling molecular dynamics for a fully atomic, solvated model of the src-SH3 protein domain. These putative TSE conformations are then subjected to a folding analysis using a coarse-grained representation of the protein and rapid discrete molecular dynamics simulations. Those conformations that fold to the native conformation with a probability (P(fold)) of approximately 0.5, constitute the true transition state. Approximately 20% of the putative TSE structures were found to have a P(fold) near 0.5, indicating that, although correct TSE conformations are populated at the free energy barrier, there is a critical need to refine this ensemble. Our simulations indicate that the true TSE conformations are compact, with a well-defined central beta sheet, in good agreement with previous experimental and theoretical studies. A structured central beta sheet was found to be present in a number of pre-TSE conformations, however, indicating that this element, although required in the transition state, does not define it uniquely. An additional tight cluster of contacts between highly conserved residues belonging to the diverging turn and second beta-sheet of the protein emerged as being critical elements of the folding nucleus. A number of commonly used order parameters to identify the transition state for folding were investigated, with the number of native Cbeta contacts displaying the most satisfactory correlation with P(fold) values.

Scanning malleable transition state ensembles: comparing theory and experiment for folding protein U1A

Using a variational free energy functional, we calculate the characteristics of the transition state ensembles (TSE) for the folding of protein U1A and investigate how they respond to thermal and mutational changes. The functional directly yields predicted chevron plots both for the wild-type protein and for various mutants. The detailed variations of the TSE and changes in chevron plots predicted by the theory agree reasonably well with the results of the experiments. We also show how to visualize the folding nuclei using 3D isodensity plots.

Ubiquitin folds through a highly polarized transition state

The small alpha/beta protein ubiquitin has been used as a model system for experimental and computational studies on protein folding for many years. Here, we present a comprehensive phi-value analysis and characterize the structure and energetics of the transition state ensemble (TSE). Twenty-seven non-disruptive mutations are made throughout the structure and a range of phi-values from zero to one are observed. The values cluster such that medium and high values and found only in the N-terminal region of the protein, whilst the C-terminal region has consistently low phi-values. In the TSE, the main alpha-helix appears to be fully formed (two phi-values which specifically probe helical structure are one) and the helix is stabilized by packing against the first beta-turn, which is partially structured. In striking comparison, the phi-values in the C-terminal region are all very low, suggesting that this region of the protein is largely unstructured in the TSE. Data are consistent with a nucleation-condensation mechanism in which there is a highly polarized folding nucleus comprising the first beta-hairpin and the alpha-helix. Data presented from the protein engineering study and phi-value analysis are compared with results from other experimental studies and also computational studies.

Gating of acetylcholine receptor channels: brownian motion across a broad transition state

Acetylcholine receptor channels (AChRs) are proteins that switch between stable "closed" and "open" conformations. In patch clamp recordings, diliganded AChR gating appears to be a simple, two-state reaction. However, mutagenesis studies indicate that during gating dozens of residues across the protein move asynchronously and are organized into rigid body gating domains ("blocks"). Moreover, there is an upper limit to the apparent channel opening rate constant. These observations suggest that the gating reaction has a broad, corrugated transition state region, with the maximum opening rate reflecting, in part, the mean first-passage time across this ensemble. Simulations reveal that a flat, isotropic energy profile for the transition state can account for many of the essential features of AChR gating. With this mechanism, concerted, local structural transitions that occur on the broad transition state ensemble give rise to fractional measures of reaction progress (Phi values) determined by rate-equilibrium free energy relationship analysis. The results suggest that the coarse-grained AChR gating conformational change propagates through the protein with dynamics that are governed by the Brownian motion of individual gating blocks.

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.

The folding of spectrin domains I: wild-type domains have the same stability but very different kinetic properties

The study of proteins with the same architecture, but different sequence has proven to be a valuable tool in the protein folding field. As a prelude to studies on the folding mechanism of spectrin domains we present the kinetic characterisation of the wild-type forms of the 15th, 16th, and 17th domains of chicken brain alpha-spectrin (referred to as R15, R16 and R17, respectively). We show that the proteins all behave in a two-state manner, with different kinetic properties. The folding rate varies remarkably between different members, with a 5000-fold variation in folding rate and 3000-fold variation in unfolding rate seen for proteins differing only 1 kcal mol(-1) in stability. We show clear evidence for significant complexity in the energy landscape of R16, which shows a change in amplitude outside the stopped-flow timescale and curvature in the unfolding arm of the chevron plot. The accompanying paper describes the characterisation of the folding pathway of this domain.