A comparison of experimental and computational methods for mapping the interactions present in the transition state for folding of FKBP12.

The folding pathway of FKBP12, a 107 residue α/ß protein, has been characterised in detail using a combination of experimental and computational techniques. FKBP12 follows a two-state model of folding in which only the denatured and native states are significantly populated; no intermediate states are detected. The refolding rate constant in water is 4 s(-1) at 25 °C. Two different experimental strategies were employed for studying the transition state for folding. In the first case, a non-mutagenic approach was used and the unfolding and refolding of the wild-type protein measured as a function of experimental conditions such as temperature, denaturant, ligand and trifluoroethanol (TFE) concentration. These data suggest a compact transition state relative to the unfolded state with some 70% of the surface area buried. The ligand-binding site, whichis mainly formed by two long loops, is largely unstructured in the transition state. TFE experiments suggest that the α-helix may be formed in the transition state. The second experimental approach involved using protein engineering techniques with φ-value analysis. Residue-specific information on the structure and energetics of the transition state can be obtained by this method. 34 mutations were made at sites throughout the protein to probe the extent of secondary and tertiary structure in the transition state. In contrast to some other proteins of this size, no element of structure is fully formed in the transition state, instead, the transition state is similar to that found for smaller, single-domain proteins, such as chymotrypsin inhibitor 2 and the SH3 domainfrom α-spectrin. For FKBP12, the central three strands of the ß-sheet (2, 4 and 5), comprise the most structured region of the transition state. In particular Val 101, which is one of the most highly buried residues and located in the middle of the central ß-strand,makes approximately 60% of its native interactions. The outer ß-strands, and the ends of the central ß-strands are formed to a lesser degree. The short α-helix is largely unstructured in the transition state as are the loops. The data are consistent with a nucleation-condensation model of folding, the nucleus of which is formed by side chains within ß-strands 2, 4 and 5 and the C-terminus of the α-helix. These residues are distant in the primary sequence, demonstrating the importance of tertiary interactions in the transition state. High-temperature molecular dynamic simulations on the unfoldingpathway of FKBP12 are in good agreement with the experimental results.

Folding pathway of FKBP12 and characterisation of the transition state.

The folding pathway of human FKBP12, a 12 kDa FK506-binding protein (immunophilin), has been characterised. Unfolding and refolding rate constants have been determined over a wide range of denaturant concentrations and data are shown to fit to a two-state model of folding in which only the denatured and native states are significantly populated, even in the absence of denaturant. This simple model for folding, in which no intermediate states are significantly populated, is further supported from stopped-flow circular dichroism experiments in which no fast "burst" phases are observed. FKBP12, with 107 residues, is the largest protein to date which folds with simple two-state kinetics in water (kF=4 s(-1)at 25 degrees C). The topological crossing of two loops in FKBP12, a structural element suggested to cause kinetic traps during folding, seems to have little effect on the folding pathway. The transition state for folding has been characterised by a series of experiments on wild-type FKBP12. Information on the thermodynamic nature of, the solvent accessibility of, and secondary structure in, the transition state was obtained from experiments measuring the unfolding and refolding rate constants as a function of temperature, denaturant concentration and trifluoroethanol concentration. In addition, unfolding and refolding studies in the presence of ligand provided information on the structure of the ligand-binding pocket in the transition state. The data suggest a compact transition state relative to the unfolded state with some 70 % of the surface area buried. The ligand-binding site, which is formed mainly by two loops, is largely unstructured in the transition state. The trifluoroethanol experiments suggest that the alpha-helix may be formed in the transition state. These results are compared with results from protein engineering studies and molecular dynamics simulations (see the accompanying paper).

Mapping the interactions present in the transition state for unfolding/folding of FKBP12.

The structure of the transition state for folding/unfolding of the immunophilin FKBP12 has been characterised using a combination of protein engineering techniques, unfolding kinetics, and molecular dynamics simulations. A total of 34 mutations were made at sites throughout the protein to probe the extent of secondary and tertiary structure in the transition state. The transition state for folding is compact compared with the unfolded state, with an approximately 30 % increase in the native solvent-accessible surface area. All of the interactions are substantially weaker in the transition state, as probed by both experiment and molecular dynamics simulations. In contrast to some other proteins of this size, no element of structure is fully formed in the transition state; instead, the transition state is similar to that found for smaller, single-domain proteins, such as chymotrypsin inhibitor 2 and the SH3 domain from alpha-spectrin. For FKBP12, the central three strands of the beta-sheet, beta-strand 2, beta-strand 4 and beta-strand 5, comprise the most structured region of the transition state. In particular Val101, which is one of the most highly buried residues and located in the middle of the central beta-strand, makes approximately 60 % of its native interactions. The outer beta-strands and the ends of the central beta-strands are formed to a lesser degree. The short alpha-helix is largely unstructured in the transition state, as are the loops. The data are consistent with a nucleation-condensation model of folding, the nucleus of which is formed by side-chains within beta-strands 2, 4 and 5, and the C terminus of the alpha-helix. The precise residues involved in the nucleus differ in the two simulated transition state ensembles, but the interacting regions of the protein are conserved. These residues are distant in the primary sequence, demonstrating the importance of tertiary interactions in the transition state. The two independently derived transition state ensembles are structurally similar, which is consistent with a Bronsted analysis confirming that the transition state is an ensemble of states close in structure.

Context-dependent nature of destabilizing mutations on the stability of FKBP12.

The context-dependent nature in which mutations affect protein stability was investigated using the FK506-binding protein, FKBP12. Thirty-four mutations were made at sites throughout the protein, including residues located in the hydrophobic core, the beta-sheet, and the solvent-exposed face of the alpha-helix. Urea-induced denaturation experiments were used to measure the change in stability of the mutants relative to that of the wild type (Delta DeltaGU-F). The results clearly show that the extent of destabilization, or stabilization, is highly context-dependent. Correlations were sought in order to link Delta DeltaGU-F to various structural parameters. The strongest correlation found was between Delta DeltaGU-F and N, the number of methyl(ene) groups within a 6 A radius of the group(s) deleted. For mutations of buried hydrophobic residues, a correlation coefficient of 0.73 (n = 16,where n is the number of points) was obtained. This increased to 0.81 (n = 24) on inclusion of mutations of partially buried hydrophobic residues. These data could be superimposed on data obtained for other proteins for which similarly detailed studies have been performed. Thus, the contribution to stability from hydrophobic side chains, independent of the extent to which a side chain is buried, can be estimated quantitatively using N. This correlation appears to be a general feature of all globular proteins. The effect on stability of mutating polar and charged residues in the alpha-helix and beta-sheet was also found to be highly context-dependent. Previous experimental and statistical studies have shown that specific side chains can stabilize the N-caps of alpha-helices in proteins. Substitutions of Ile56 to Thr and Asp at the N-cap of the alpha-helix of FKBP12, however, were found to be highly destabilizing. Thus, the intrinsic propensities of an amino acid for a particular element of secondary structure can easily be outweighed by tertiary packing factors. This study highlights the importance of packing density in determining the contribution of a residue to protein stability. This is the most important factor that should be taken into consideration in protein design.