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.     

Nucleation and the transition state of the SH3 domain

We present a verified computational model of the SH3 domain transition state (TS) ensemble. This model was built for three separate SH3 domains using experimental phi-values as structural constraints in all-atom protein folding simulations. While averaging over all conformations incorrectly considers non-TS conformations as transition states, quantifying structures as pre-TS, TS, and post-TS by measurement of their transmission coefficient ("probability to fold", or p(fold)) allows for rigorous conclusions regarding the structure of the folding nucleus and a full mechanistic analysis of the folding process. Through analysis of the TS, we observe a highly polarized nucleus in which many residues are solvent-exposed. Mechanistic analysis suggests the hydrophobic core forms largely after an early nucleation step. SH3 presents an ideal system for studying the nucleation-condensation mechanism and highlights the synergistic relationship between experiment and simulation in the study of protein folding.     

Experimental investigation of the frequency and substitution dependence of negative phi-values in two-state proteins

Negative phi-values, which arise, for example, when a mutation stabilizes the folding transition state while destabilizing the native state, have been the focus of significant theoretical interest. Here we survey the experimental folding kinetics literature to ascertain the frequency with which negative phi-values occur in two-state proteins and describe the detailed experimental characterization of a negative phi-value previously reported to be among the most statistically significant. We find that, while almost 9% of more than 500 reported phi-values (from a set of 16, well-characterized two-state proteins) fall below zero, many of these do not represent statistically significant observations. For example, only 6% of the phi-values for which estimates of precision are available fall even one reported "error bar" below zero, and only 4% are simultaneously negative, significant at this level and associated with free energy changes at or above 2.5 kJ/mol (below which phi-value analysis is widely considered unreliable). Moreover, given the asymmetric distribution of phi-values around zero and given that reported error bars may significantly underestimate true confidence intervals, the actual number of negative phi-values may be much smaller still. We have also performed detailed characterization of one of the most statistically significant negative phi-values reported in the literature to date, the V55F mutant of FynSH3. We find that substitution of the wild-type valine to other hydrophobic residues often increases folding rates without significantly altering folding free energy. This in turn leads to poorly defined phi-values, some of which are formally negative but only one or two of which fall statistically significantly below zero. In contrast, substitution to polar residues significantly destabilizes both the transition and native states, generally producing small but statistically significant positive phi-values of approximately 0.1. Thus, unlike other previously characterized phi-values, the negative phi-value associated with position 55 of the FynSH3 domain appears to be strongly dependent on the substitution employed to measure it, suggesting that subtlety will be required in order to develop a theoretical model of such behavior.     

Investigation of folding/unfolding kinetics of apomyoglobin

Apomyoglobin kinetic and equilibrium unfolding and folding processes were studied at pH 6.2, 11 degrees C by stopped-flow tryptophan fluorescence. There are two distinct consecutive processes in apomyoglobin folding process, namely, the protein fast transition between the unfolded (U) and an intermediate (I) states (U <----> I) and slow transition between the intermediate and the native (N) states (I <----> N). Accumulation of the intermediate state was observed in the wide range of urea concentrations. The presence of the intermediate state was shown even beyond the middle transition on the unfolding limb. The dependence of observed folding/unfolding rates on urea concentration (chevron plot) was obtained. The shape of this dependence was compared with that of two-state proteins, folding from the U to N state.     

Conformational plasticity in folding of the split beta-alpha-beta protein S6: evidence for burst-phase disruption of the native state

An increasing number of folding studies of two-state proteins shows that point mutations sometimes change the kinetic m-values, leading to kinks and curves in the chevron plots. The molecular origin of these changes is yet unclear although it is speculated that they are linked to structural rearrangement of the transition state or to accumulation of meta-stable intermediates. To shed more light on this issue, we present here a combined m and phi-value analysis of the split beta-alpha-beta protein S6. Wild-type S6 displays classical two-state kinetics with v-shaped chevron plot, but a majority of its mutants display distinct m-value changes or curved chevrons. We observe that this kinetic aberration of S6 is linked to mutations that are clustered in distinct regions of the native structure. The most pronounced changes, i.e. decrease in the m-value for the unfolding rate constant, are seen upon truncation of interactions between the N and C termini, whereas mutations in the centre of the hydrophobic core show smaller or even opposed effects. As a consequence, the calculated phi-values display a systematic increase upon addition of denaturant. In the case of S6, the phenomenon seems to arise from a general plasticity of the different species on the folding pathway. That is, the structure of the denatured ensemble, the transition state, and the native ground-state for unfolding seem to change upon mutation. From these changes, it is concluded that interactions spanning the centre of the hydrophobic core form early in folding, whereas the entropically disfavoured interactions linking the N and C termini consolidate very late, mainly on the down-hill-side of the folding barrier.     

Characterization of the hydrodynamic properties of the folding transition state of an SH3 domain by magnetization transfer NMR spectroscopy

Protein folding kinetic data have been obtained for the marginally stable N-terminal Src homology 3 domain of the Drosophila protein drk (drkN SH3) in an investigation of the hydrodynamic properties of its folding transition state. Due to the presence of NMR resonances of both folded and unfolded states at equilibrium, kinetic data can be derived from NMR magnetization transfer techniques under equilibrium conditions. Kinetic analysis as a function of urea (less than approximately 1 M) and glycerol enables determination of alpha values, measures of the energetic sensitivity of the transition state to the perturbation relative to the end states of the protein folding reaction (the folded and unfolded states). Both end states have previously been studied experimentally by NMR spectroscopic and other biophysical methods in great detail and under nondenaturing conditions. Combining these results with the kinetic folding data obtained here, we can characterize the folding transition state without requiring empirical models for the unfolded state structure. We are thus able to give a reliable measure of the solvent-accessible surface area of the transition state of the drkN SH3 domain (4730 +/- 360 A(2)) based on urea titration data. Glycerol titration data give similar results and additionally demonstrate that folding of this SH3 domain is dependent on solvent viscosity, which is indicative of at least partial hydration of the transition state. Because SH3 domains appear to fold by a common folding mechanism, the data presented here provide valuable insight into the transition states of the drkN and other SH3 domains.     

Folding/Unfolding Kinetics of Apomyoglobin

The equilibrium and kinetic folding/unfolding of apomyoglobin (ApoMb) were studied at pH 6.2, 11 °C by recording tryptophan fluorescence. The equilibrium unfolding of ApoMb in the presence of urea was shown to involve accumulation of an intermediate state, which had a higher fluorescence intensity as compared with the native and unfolded states. The folding proceeded through two kinetic phases, a rapid transition from the unfolded to the intermediate state and a slow transition from the intermediate to the native state. The accumulation of the kinetic intermediate state was observed in a wide range of urea concentrations. The intermediate was detected even in the region corresponding to the unfolding limb of the chevron plot. Urea concentration dependence was obtained for the observed folding/unfolding rate. The shape of the dependence was compared with that of two-state proteins characterized by a direct transition from the unfolded to the native state.     

Comparison of the transition states for folding of two Ig-like proteins from different superfamilies

In the "fold approach" proteins with a similar fold but different sequences are compared in order to investigate the relationship between native state structure and folding behaviour. Here we compare the properties of the transition states for folding of TI I27, the 27th immunoglobulin domain from human cardiac titin, and that of TNfn3, the third fibronectin type III domain from human tenascin. Experimental phi-values were used as restraints in molecular dynamics simulations to determine the structures that make up the transition state ensembles (TSEs) for folding of the two proteins. The restrained simulations that we present allow a detailed structural comparison of the two TSEs to be made. Further calculations show explicitly that for both proteins the formation of the interactions involving the residues in the folding nucleus is sufficient for the establishment of the topology of the Ig-like fold. We found that, although the folding nuclei of the two proteins are similar, the packing of the folding nucleus of TI I27 is much tighter than that of TNfn3, reflecting the higher experimental phi-values and beta(T) (Tanford Beta) of TI I27. These results suggest that the folding nucleus can be significantly deformed to accommodate extensive sequence variation while conserving the same folding mechanism.     

Probing the influence on folding behavior of structurally conserved core residues in P. aeruginosa apo-azurin

The effects on folding kinetics and equilibrium stability of core mutations in the apo-mutant C112S of azurin from Pseudomonas aeruginosa were studied. A number of conserved residues within the cupredoxin family were recognized by sequential alignment as constituting a common hydrophobic core: I7, F15, L33, W48, F110, L50, V95, and V31. Of these, I7, V31, L33, and L50 were mutated for the purpose of obtaining information on the transition state and a potential folding nucleus. In addition, residue V5 in the immediate vicinity of the common core, as well as T52, separate from the core, were mutated as controls. All mutants exhibited a nonlinear dependence of activation free energy of folding on denaturant concentration, although the refolding kinetics of the V31A/C112S mutant indicated that the V31A mutation destabilizes the transition state enough to allow folding via a parallel transition state ensemble. Phi-values could be calculated for three of the six mutants, V31A/C112S, L33A/C112S, and L50A/C112S, and the fractional values of 0.63, 0.33, and 0.50 (respectively) obtained at 0.5 M GdmCl suggest that these residues are important for stabilizing the transition state. Furthermore, a linear dependence of ln k(obs)(H2O) on DeltaG(U-N)(H2O) of the core mutations and the putative involvement of ground-state effects suggest the presence of native-like residual interactions in the denatured state that bias this ensemble toward a folding-competent state.     

An expanded view of the protein folding landscape of PDZ domains

Most protein domains fold in an apparently co-operative and two-state manner with only the native and denatured states significantly populated at any experimental condition. However, the protein folding energy landscape is often rugged and different transition states may be rate limiting for the folding reaction under different conditions, as seen for the PDZ protein domain family. We have here analyzed the folding kinetics of two PDZ domains and found that a previously undetected third transition state is rate limiting under conditions that stabilize the native state relative to the denatured state. In light of these results, we have re-analyzed previous folding data on PDZ domains and present a unified folding mechanism with three distinct transition states separated by two high-energy intermediates. Our data show that sequence composition tunes the relative stabilities of folding transition states within the PDZ family, while the overall mechanism is determined by topology. This model captures the kinetic folding mechanism of all PDZ domains studied to date.     

Confirmation of the hierarchical folding of RNase H: a protein engineering study.

The kinetic intermediate of RNase H is structured in a core region of the protein. To probe the role of this intermediate in the folding of RNase H, the folding kinetics of mutant proteins with altered native state stabilities were investigated. Mutations within the folding core destabilize the kinetic intermediate and slow refolding in a manner consistent with an obligatory intermediate model. Mutations outside of the folding core, however, do not affect the stability of the kinetic intermediate but do perturb the native state and transition state. These results indicate that interactions formed in the intermediate persist in the transition and native states and that RNase H folds through a hierarchical mechanism.     

Parallel protein-unfolding pathways revealed and mapped

Theoretical studies of protein folding suggest that multiple folding pathways should exist, but there is little experimental evidence to support this. Here we demonstrate changes in the flux between different transition states on parallel folding pathways, resulting in unprecedented upward curvature in the denaturant-dependent unfolding kinetics of a beta-sandwich protein. As denaturant concentration increases, the highly compact transition state of one pathway becomes destabilized and the dominant flux of protein molecules shifts toward another pathway with a less structured transition state. Furthermore, point mutations alter the relative accessibility of the pathways, allowing the structure of two transition states on separate, direct folding pathways to be mapped by systematic Phi-value analysis. It has been suggested that pathways with diffuse rather than localized transition states are evolutionarily selected to prevent misfolding, and indeed we find that the transition state favored at high concentrations of denaturant is more polarized than the physiologically relevant one.     

Intramolecular cross-linking evaluated as a structural probe of the protein folding transition state

We examine the utility of intramolecular covalent cross-linking to identify the structure present in the folding transition state. In mammalian ubiquitin, cysteine residues located across two beta-strands are cross-linked with dichloroacetone. The kinetic effects of these covalent cross-links in ubiquitin, and engineered disulfide bonds in src SH3 (Grantcharova, V. P., Riddle, D. S., and Baker, D. (2000) Proc. Natl. Acad. Sci. U.S.A. 97, 7084-7089), are compared to the results of psi-analysis where strand association is stabilized by metal ion binding to engineered bihistidine sites (Krantz, B. A., Dothager, R. S., and Sosnick, T. R. (2004) J. Mol. Biol. 337, 463-75) at the same positions. The results for the two methods agree at some of the sites. The cross-linking phi crosslink-values agree with their corresponding psi-values when they have both have values of zero or one, which represent the absence and presence of native structure, respectively. When phi crosslink > psi, the apparent inconsistency is rationalized by the difference between each method's mode of stabilization; cross-linking reduces the configurational entropy of the unfolded state whereas metal binding directly stabilizes the native state. However, when the cross-linking phi-values are smaller than their corresponding psi-values, the apparent underestimation of structure formation is difficult to rationalize while retaining the assumption that the cross-link exclusively affects the entropy of the unfolded state. The interpretation also is problematic for data on cross-links located across strands which are not hairpins, and hence, these sites are likely to be of limited utility in folding studies. We conclude that cross-linking data for sites on hairpins generally report on the amount of structure formed within the enclosed loop while the metal binding data report on the amount structure formed at the site itself.     

Commitment and nucleation in the protein G transition state

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

Mutational analysis of the folding transition state of the C-terminal domain of ribosomal protein L9: a protein with an unusual beta-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.     

Mapping the folding pathway of an immunoglobulin domain: structural detail from Phi value analysis and movement of the transition state

BACKGROUND: Do proteins that have the same structure fold by the same pathway even when they are unrelated in sequence? To address this question, we are comparing the folding of a number of different immunoglobulin-like proteins. Here, we present a detailed protein engineering phi value analysis of the folding pathway of TI I27, an immunoglobulin domain from human cardiac titin. RESULTS: TI I27 folds rapidly via a kinetic intermediate that is destabilized by most mutations. The transition state for folding is remarkably native-like in terms of solvent accessibility. We use phi value analysis to map this transition state and show that it is highly structured; only a few residues close to the N-terminal region of the protein remain completely unfolded. Interestingly, most mutations cause the transition state to become less native-like. This anti-Hammond behavior can be used as a novel means of obtaining additional structural information about the transition state. CONCLUSIONS: The residues that are involved in nucleating the folding of TI I27 are structurally equivalent to the residues that form the folding nucleus in an evolutionary unrelated fibronectin type III protein. These residues form part of the common structural core of Ig-like domains. The data support the hypothesis that interactions essential for defining the structure of these beta sandwich proteins are also important in nucleation of folding.     

Hydrophobic core packing in the SH3 domain folding transition state.

How tightly packed is the hydrophobic core of a folding transition state structure? We have addressed this question by characterizing the effects on folding kinetics of > 40 substitutions of both large and small amino acids in the hydrophobic core of the Fyn SH3 domain. Our results show that residues at three positions, which we designate as the 'core folding nucleus', are tightly packed in the transition state, and substitutions at these positions cause the largest changes in the folding rate. The other six positions examined appear to be loosely packed; thus, substitutions at these positions with larger hydrophobic residues generally accelerate folding, presumably by increasing the rate of nonspecific hydrophobic collapse. Surprisingly, the folding rate can be greatly accelerated by residues that also significantly destabilize the native state structure. Furthermore, mutants with identical thermodynamic stability can differ by up to 55-fold in their folding rates. These results highlight the importance of hydrophobic core composition, as opposed to only topology, in determining the folding rate of a protein. They also provide a new explanation for the 'abnormal' phi-values observed in many protein folding kinetics studies.     

Posttransition state desolvation of the hydrophobic core of the src-SH3 protein domain

The folding thermodynamics of the src-SH3 protein domain were characterized under refolding conditions through biased fully atomic molecular dynamics simulations with explicit solvent. The calculated free energy surfaces along several reaction coordinates revealed two barriers. The first, larger barrier was identified as the transition state barrier for folding, associated with the formation of the first hydrophobic sheet of the protein. phi values calculated from structures residing at the transition state barrier agree well with experimental phi values. The microscopic information obtained from our simulations allowed us to unambiguously assign intermediate phi values as the result of multiple folding pathways. The second, smaller barrier occurs later in the folding process and is associated with the cooperative expulsion of water molecules between the hydrophobic sheets of the protein. This posttransition state desolvation barrier cannot be observed through traditional folding experiments, but is found to be critical to the correct packing of the hydrophobic core in the final stages of folding. Hydrogen exchange and NMR experiments are suggested to probe this barrier.     

Influence of Trp mutation on native, intermediate, and transition states of goat alpha-lactalbumin: an equilibrium and kinetic study

Equilibrium circular dichroism and kinetic stopped-flow fluorescence studies on the stability and the folding kinetics of a set of Trp to Phe mutants of goat alpha-lactalbumin (GLA) were used to characterize the native, intermediate, and transition states of these constructs. GLA contains four tryptophan residues, three of which, Trp26, Trp104, and Trp118, are located in the alpha-domain, while the fourth, Trp60, is located in the beta-domain. Trp26, Trp60, and Trp104 are part of a hydrophobic cluster, whereas Trp118 is situated in a more flexible region near the C-terminal end of the protein. In each case, the mutation leads to a reduction in the overall stability, but only for W26F and W60F is an equilibrium intermediate observed in guanidine hydrochloride-induced unfolding experiments. In kinetic refolding experiments, however, for all samples a burst phase is observed, the amplitude of which depends on the specific mutation. Refolding and unfolding kinetics can adequately be described by a sequential three-state mechanism. phi value analysis showed that the local structure around Trp26, Trp60, and Trp104 is formed in the intermediate and in the transition state of the folding reaction, while around Trp118 no persistent native contacts are observed. From these findings, we conclude that, although hydrophobicity is a major driving force for folding, minor steric changes induced by point mutation can considerably influence the overall stability and the folding process of the protein.