The Structure of the Major Transition State for Folding of an FF Domain from Experiment and Simulation

Fractional 13C enrichment of isolated carbons using [1-13C]- or [2- 13C]-glucose facilitates the accurate measurement of dynamics at backbone Calpha and side-chain methyl positions in proteins

A simple labeling approach is presented based on protein expression in [1-(13)C]- or [2-(13)C]-glucose containing media that produces molecules enriched at methyl carbon positions or backbone C(alpha) sites, respectively. All of the methyl groups, with the exception of Thr and Ile(delta1) are produced with isolated (13)C spins (i.e., no (13)C-(13)C one bond couplings), facilitating studies of dynamics through the use of spin-spin relaxation experiments without artifacts introduced by evolution due to large homonuclear scalar couplings. Carbon-alpha sites are labeled without concomitant labeling at C(beta) positions for 17 of the common 20 amino acids and there are no cases for which (13)C(alpha)-(13)CO spin pairs are observed. A large number of probes are thus available for the study of protein dynamics with the results obtained complimenting those from more traditional backbone (15)N studies. The utility of the labeling is established by recording (13)C R (1rho) and CPMG-based experiments on a number of different protein systems.

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

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

Probing a moving target with a plastic unfolding intermediate of an ankyrin-repeat protein

Repeat proteins are composed of tandem arrays of 30- to 40-residue structural motifs and are characterized by short-range interactions between residues close in sequence. Here we have investigated the equilibrium unfolding of D34, a 426-residue fragment of ankyrinR that comprises 12 ankyrin repeats. We show that D34 unfolds via an intermediate in which the C-terminal half of the protein is structured and the N-terminal half is unstructured. Surprisingly, however, we find that we change the unfolding process when we attempt to probe it. Single-site, moderately destabilizing mutations at the C terminus result in different intermediates dominating. The closer to the C terminus the mutation, the fewer repeats are structured in the intermediate; thus, structure in the intermediate frays from the site of the mutation. This behavior contrasts with the robust unfolding of globular proteins in which mutations can destabilize an intermediate but do not cause a different intermediate to be populated. We suggest that, for large repeat arrays, the energy landscape is very rough, with many different low-energy species containing varying numbers of folded modules so the species that dominates can be altered easily by single, conservative mutations. The multiplicity of partly folded states populated in the equilibrium unfolding of D34 is also mirrored by the kinetic folding mechanism of ankyrin-repeat proteins in which we have observed that parallel pathways are accessible from different initiation sites in the structure.

A single-quantum methyl 13C-relaxation dispersion experiment with improved sensitivity

A pulse sequence is described for recording single-quantum (13)C-methyl relaxation dispersion profiles of (13)C-selectively labeled methyl groups in proteins that offers significant improvements in sensitivity relative to existing approaches where initial magnetization derives from (13)C polarization. Sensitivity gains in the new experiment are achieved by making use of polarization from (1)H spins and (1)H --> (13)C --> (1)H type magnetization transfers. Its utility has been established by applications involving three different protein systems ranging in molecular weight from 8 to 28 kDa, produced using a number of different selective labeling approaches. In all cases exchange parameters from both (13)C-->(1)H and (1)H --> (13)C --> (1)H classes of experiment are in good agreement, with gains in sensitivity of between 1.7 and 4-fold realized using the new scheme.

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.

Folding transition-state and denatured-state ensembles of FSD-1 from folding and unfolding simulations

Characterization of the folding transition-state ensemble and the denatured-state ensemble is an important step toward a full elucidation of protein folding mechanisms. We report herein an investigation of the free-energy landscape of FSD-1 protein by a total of four sets of folding and unfolding molecular dynamics simulations with explicit solvent. The transition-state ensemble was initially identified from unfolding simulations at 500 K and was verified by simulations at 300 K starting from the ensemble structures. The denatured-state ensemble and the early-stage folding were studied by a combination of unfolding simulations at 500 K and folding simulations at 300 K starting from the extended conformation. A common feature of the transition-state ensemble was the substantial formation of the native secondary structures, including both the alpha-helix and beta-sheet, with partial exposure of the hydrophobic core in the solvent. Both the native and non-native secondary structures were observed in the denatured-state ensemble and early-stage folding, consistent with the smooth experimental melting curve. Interestingly, the contact orders of the transition-state ensemble structures were similar to that of the native structure and were notably lower than those of the compact structures found in early-stage folding, implying that chain and topological entropy might play significant roles in protein folding. Implications for FSD-1 folding mechanisms and the rate-limiting step are discussed. Analyses further revealed interesting non-native interactions in the denatured-state ensemble and early-stage folding and the possibility that destabilization of these interactions could help to enhance the stability and folding rate of the protein.

A PDZ domain recapitulates a unifying mechanism for protein folding

A unifying view has been recently proposed according to which the classical diffusion-collision and nucleation-condensation models may represent two extreme manifestations of an underlying common mechanism for the folding of small globular proteins. We report here the characterization of the folding process of the PDZ domain, a protein that recapitulates the three canonical steps involved in this unifying mechanism, namely: (i) the early formation of a weak nucleus that determines the native-like topology of a large portion of the structure, (ii) a global collapse of the entire polypeptide chain, and (iii) the consolidation of the remaining partially structured regions to achieve the native state conformation. These steps, which are clearly detectable in the PDZ domain investigated here, may be difficult to distinguish experimentally in other proteins, which would thus appear to follow one of the two limiting mechanisms. The analysis of the (un)folding kinetics for other three-state proteins (when available) appears consistent with the predictions ensuing from this unifying mechanism, thus providing a powerful validation of its general nature.

Key Role of Coulombic Interactions for the Folding Transition State of the Cold Shock Protein

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.

Transition-states in protein folding kinetics: the structural interpretation of Phi values

Phi values are experimental measures of the effects of mutations on the folding kinetics of a protein. A central question is what structural information Phi values give about the transition-state of folding. Traditionally, a Phi value is interpreted as representing the "nativeness" of a mutated residue in the transition-state. However, this interpretation is often problematic. We present here a better structural interpretation of Phi values for mutations within a given helix. Our interpretation is based on a simple physical model that distinguishes between secondary and tertiary free energy contributions of helical residues. From a linear fit of the model to experimental data, we obtain two structural parameters: the extent of helix formation in the transition-state, and the nativeness of tertiary interactions in the transition-state. We apply the model to all proteins with well-characterized helices for which more than 10 Phi values are available: protein A, CI2, and protein L. The model is simple to apply to experimental data, captures nonclassical Phi values <0 or >1 in these helices, and explains how different mutations at a given site can lead to different Phi values.

The side chain interaction index as a tool for predicting fast-folding elements and the structure and stability of engineered peptides

The side chain interaction index (SCII) is a method of calculating the propensity for short-range interactions among side chains within a peptide sequence. Here, it is shown that the SCII values of secondary structure elements that have been shown to fold early and independently cluster separately from those of structures that fold later and/or are dependent on long-range interactions. In addition, the SCII values of engineered peptides that spontaneously adopt a particular desired fold in solution are significantly different from those of engineered peptides that fail to exhibit a stable conformation. Thus, the SCII, as a measure of local structural stability, constitutes a useful tool in folding prediction and in protein/peptide engineering. A program that allows rapid calculation of SCII values is presented.

Phi-analysis of the folding of the B domain of protein A using multiple optical probes

We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.

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.

The folding pathway of spectrin R17 from experiment and simulation: using experimentally validated MD simulations to characterize States hinted at by experiment

We present an experimental and computational analysis of the folding pathway of the 17th domain of chicken brain alpha-spectrin, R17. Wild-type R17 folds in a two-state manner and the chevron plot (plot of the logarithm of the observed rate constant against concentration of urea) shows essentially linear folding and unfolding arms. A number of mutant proteins, however, show a change in slope of the unfolding arm at high concentration of denaturant, hinting at complexity in the folding landscape. Through a combination of mutational studies and high temperature molecular dynamics simulations we show that the folding of R17 can be described by a model with two sequential transition states separated by an intermediate species. The rate limiting transition state for folding in water has been characterized both through experimental Phi-value analysis and by simulation. In contrast, a detailed analysis of the transition state predicted to dominate under highly denaturing conditions is only possible by simulation.