Folding of a de novo designed native-like four-helix bundle protein

The folding of a model native-like dimeric four-helix bundle protein, (alpha(2))(2), was investigated using guanidine hydrochloride, hydrostatic pressure, and low temperature. Unfolding by guanidine hydrochloride followed by circular dichroism and intrinsic fluorescence spectroscopy revealed a highly cooperative transition between the native-like and unfolded states, with free energy of unfolding determined from CD data, DeltaG(unf) = 14.3 +/- 0.8 kcal/mol. However, CD and intrinsic fluorescence data were not superimposable, indicating the presence of an intermediate state during the folding transition. To stabilize the folding intermediate, we used hydrostatic pressure and low temperature. In both cases, dissociation of the dimeric native-like (alpha(2))(2) into folded monomers (alpha(2)) was observed. van't Hoff analysis of the low temperature experiments, assuming a two-state dimer 171-monomer transition, yielded a free energy of dissociation of (alpha(2))(2) of DeltaG(diss) = 11.4 +/- 0.4 kcal/mol, in good agreement with the free energy determined from pressure dissociation experiments (DeltaG(diss) = 10.5 +/- 0.1 kcal/mol). Binding of the hydrophobic fluorescent probe 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid (bis-ANS) to the pressure- and cold-dissociated states of (alpha(2))(2) indicated the existence of molten-globule monomers. In conclusion, we demonstrate that the folding pathway of (alpha(2))(2) can be described by a three-state transition including a monomeric molten globule-like state.     

Shape of the free energy barriers for protein folding probed by multiple perturbation analysis

The characterization of the free energy barriers has been a major goal in studies on the mechanism of protein folding. Testing the effect of mutations or denaturants on protein folding reactions revealed that transition state movement is rare, suggesting that folding barriers are robust and narrow maxima on the free energy landscape. Here we demonstrate that the application of multiple perturbations allows the observation of small transition state movements that escape detection in single perturbation experiments. We used tendamistat as a model protein to test the broadness of the free energy barriers. Tendamistat folds over two consecutive transition states and through a high-energy intermediate. Measuring the combined effect of temperature and denaturant on the position of the transition state in the wild-type protein and in several mutants revealed that the early transition state shows significant transition state movement. Its accessible surface area state becomes more native-like with destabilization of the native state by temperature. To the same extent, the entropy of the early transition state becomes more native-like with increasing denaturant concentration, in accordance with Hammond behavior. The position of the late transition state, in contrast, is much less sensitive to the applied perturbations. These results suggest that the barriers in protein folding become increasingly narrow as the folding polypeptide chain approaches the native state.     

Temperature dependence of the free energy landscape of the src-SH3 protein domain

The temperature dependence of the free energy landscape of the src-SH3 protein domain is investigated through fully atomic simulations in explicit solvent. Simulations are performed above and below the folding transition temperature, enabling an analysis of both protein folding and unfolding. The transition state for folding and unfolding, identified from the free energy surfaces, is found to be very similar, with structure in the central hydrophobic sheet and little structure throughout the rest of the protein. This is a result of a polarized folding (unfolding) mechanism involving early formation (late loss) of the central hydrophobic sheet at the transition state. Unfolding simulations map qualitatively well onto low-temperature free energy surfaces but appear, however, to miss important features observed in folding simulations. In particular, details of the folding mechanism involving the opening and closing of the hydrophobic core are not captured by unfolding simulations performed under strongly denaturing conditions. In addition, free energy surfaces at high temperatures do not display a desolvation barrier found at lower temperatures, involving the expulsion of water molecules from the hydrophobic core.     

Exploring structures in protein folding funnels with free energy functionals: the transition state ensemble

We use free energy functionals that account for the partial ordering of residues in the transition state ensemble to characterize the free energy surfaces for fast folding proteins. We concentrate on chymotrypsin inhibitor and lambda-repressor. We show how the explicit cooperativity that can arise from many body forces, such as side-chain ordering or hydrophobic surface burial, determines the crossover from folding with a large delocalized nucleus and the specific small classical nucleus of the type envisioned in nucleation growth scenarios. We compare the structural correlations present in the transition state ensemble obtained from free energy functionals with those inferred from experiment using extrathermodynamic free energy relations for folding time obtained via protein engineering kinetics experiments. We also use the free energy functionals to examine both the size of barriers and multidimensional representations of the free energy profiles in order to address the question of appropriate reaction coordinates for folding.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Effects of a helix substitution on the folding mechanism of bovine alpha-lactalbumin

The structure, stability, and unfolding-refolding kinetics of a chimeric protein, in which the amino acid sequence of the flexible loop region (residues 105-110) comes from equine lysozyme and the remainder of the sequence comes from bovine alpha-lactalbumin were studied by circular dichroism spectroscopy and stopped-flow measurements, and the results were compared with those of bovine alpha-lactalbumin. The substitution of the flexible loop in bovine alpha-lactalbumin with the helix D of equine lysozyme destabilizes the molten globule state, although the native state is significantly stabilized by substitution of the flexible loop region. The kinetic refolding and unfolding experiments showed that the chimeric protein refolds significantly faster and unfolds substantially slower than bovine alpha-lactalbumin. To characterize the transition state between the molten globule and the native states, we investigated the guanidine hydrochloride concentration dependence of the rate constants of refolding and unfolding. Despite the significant differences in the stabilities of both the molten globule and native states between the chimeric protein and bovine alpha-lactalbumin, the free energy level of the transition state is not affected by the amino acid substitution in the flexible loop region. Our results suggest that the destabilization in the molten globule state of the chimeric protein is caused by the disruption of the non-native interaction in the flexible loop region and that the disruption of the non-native interaction reduces the free energy barrier of refolding. We conclude that the non-native interaction in the molten globule state may act as a kinetic trap for the folding of alpha-lactalbumin.     

Protein folding kinetics beyond the phi value: using multiple amino acid substitutions to investigate the structure of the SH3 domain folding transition state

The SH3 domain folding transition state structure contains two well-ordered turn regions, known as the diverging turn and the distal loop. In the Src SH3 domain transition state, these regions are stabilized by a hydrogen bond between Glu30 in the diverging turn and Ser47 in the distal loop. We have examined the effects on folding kinetics of amino acid substitutions at the homologous positions (Glu24 and Ser41) in the Fyn SH3 domain. In contrast to most other folding kinetics studies which have focused primarily on non-disruptive substitutions with Ala or Gly, here we have examined the effects of substitutions with diverse amino acid residues. Using this approach, we demonstrate that the transition state structure is generally tolerant to amino acid substitutions. We also uncover a unique role for Ser at position 41 in facilitating folding of the distal loop, which can only be replicated by Asp at the same position. Both these residues appear to accelerate folding through the formation of short-range side-chain to backbone hydrogen bonds. The folding of the diverging turn region is shown to be driven primarily by local interactions. The diverging turn and distal loop regions are found to interact in the transition state structure, but only in the context of particular mutant backgrounds. This work demonstrates that studying the effects of a variety of amino acid substitutions on protein folding kinetics can provide unique insights into folding mechanisms which cannot be obtained by standard Phi value analysis.     

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.     

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.     

Interpretation of protein folding psi values

The structural characterization of transition states is essential for understanding the mechanism of protein folding. Analyzing the effect of mutations on protein stability and folding kinetics in phi-value analysis is commonly used to gain information about the presence of side-chain interactions in transition states. Recently, specific binding of ligands to engineered binding sites was applied to monitor the formation of local structures in transition states (psi analysis). A surprising result from psi analysis was the presence of parallel folding pathways in all reported studies and a major discrepancy between phi and psi values measured in the same protein. Here, we show that psi values cannot be analyzed in the same way as other rate-equilibrium free energy relationships due to the involvement of bimolecular reactions that may have different dissociation constants for the native, unfolded and transition state. As a consequence, psi values reflect the relative binding energy (kappa) of the transition state only for the extreme values of kappa=0 or kappa=1. In all other cases, non-linear rate-equilibrium free-energy relationships (Leffler plots) are observed. This apparently indicates the presence of parallel folding pathways even if folding occurs over a single homogeneous transition state. Consequently, the results from Leffler plots do not yield information about the structural properties of the transition state. This explains the lack of agreement between results from psi analysis and other methods used to characterize protein folding transition states. We further show that the same considerations apply for the analysis of the effect of pH on protein folding.     

Does trifluoroethanol affect folding pathways and can it be used as a probe of structure in transition states?

Nonaqueous co-solvents, particularly 2,2,2-trifluoroethanol (TFE), have been used as tools to study protein folding. By analyzing FKBP12, an alpha/beta-protein that folds with two-state kinetics, we have been able to address three key questions concerning the use of TFE. First, does TFE perturb the folding pathway? Second, can the observed changes in the rate of folding and unfolding in TFE be attributed to a change in free energy of a single state? Finally, can TFE be used to infer information on secondary structure formation in the transition state? Protein engineering experiments on FKBP12, coupled with folding and unfolding experiments in 0% and 9.6% TFE, conclusively show that TFE does not perturb the folding pathway of this protein. Our results also suggest that the changes in folding and unfolding rates observed in 9.6% TFE are due to a global effect of TFE on the protein, rather than the stabilization of any elements of secondary structure in the transition state. Thus, studies with TFE and other co-solvents can be accurately interpreted only when combined with other techniques.     

Conformational changes in serpins: I. The native and cleaved conformations of alpha(1)-antitrypsin

The serpins (SERine Proteinase INhibitors) are a family of proteins with important physiological roles, including but not limited to the inhibition of chymotrypsin-like serine proteinases. The inhibitory mechan- ism involves a large conformational change known as the S-->R (stressed-->relaxed) transition. The largest structural differences occur in a region around the scissile bond called the reactive centre loop: In the native (S) state, the reactive centre is exposed, and is free to interact with proteinases. In inhibitory serpins, in the cleaved (R) state the reactive centre loop forms an additional strand within the beta-sheet. The latent state is an uncleaved state in which the intact reactive centre loop is integrated into the A sheet as in the cleaved form, to give an alternative R state.The serpin structures illustrate detailed control of conformation within a single protein. Serpins are also an unusual family of proteins in which homologues have native states with different folding topologies. Determination of the structures of inhibitory serpins in multiple conformational states permits a detailed analysis of the mechanism of the S-->R transition, and of the way in which a single sequence can form two stabilised states of different topology.Here we compare the conformations of alpha(1)-antitrypsin in native and cleaved states. Many protein conformational changes involve relative motions of large rigid subunits. We determine the rigid subunits of alpha(1)-antitrypsin and analyse the changes in their relative position and orientation. Knowing that the conformational change is initiated by cleavage at the reactive centre, we describe a mechanism of the S-->R transition as a logical sequence of mechanical effects, even though the transition likely proceeds in a concerted manner.     

Conformational changes in serpins: I. The native and cleaved conformations of alpha(1)-antitrypsin

The serpins (SERine Proteinase INhibitors) are a family of proteins with important physiological roles, including but not limited to the inhibition of chymotrypsin-like serine proteinases. The inhibitory mechan- ism involves a large conformational change known as the S-->R (stressed-->relaxed) transition. The largest structural differences occur in a region around the scissile bond called the reactive centre loop: In the native (S) state, the reactive centre is exposed, and is free to interact with proteinases. In inhibitory serpins, in the cleaved (R) state the reactive centre loop forms an additional strand within the beta-sheet. The latent state is an uncleaved state in which the intact reactive centre loop is integrated into the A sheet as in the cleaved form, to give an alternative R state.The serpin structures illustrate detailed control of conformation within a single protein. Serpins are also an unusual family of proteins in which homologues have native states with different folding topologies. Determination of the structures of inhibitory serpins in multiple conformational states permits a detailed analysis of the mechanism of the S-->R transition, and of the way in which a single sequence can form two stabilised states of different topology.Here we compare the conformations of alpha(1)-antitrypsin in native and cleaved states. Many protein conformational changes involve relative motions of large rigid subunits. We determine the rigid subunits of alpha(1)-antitrypsin and analyse the changes in their relative position and orientation. Knowing that the conformational change is initiated by cleavage at the reactive centre, we describe a mechanism of the S-->R transition as a logical sequence of mechanical effects, even though the transition likely proceeds in a concerted manner.     

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.     

Some recommendations for the practitioner to improve the precision of experimentally determined protein folding rates and phi values

The mechanism by which proteins fold from an initially random conformation into a functional, native structure remains a major unsolved question in molecular biology. Of particular interest to the protein folding community is the structure that the protein adopts in the folding transition state (the highest free energy state on the pathway from unfolded to folded), as that state forms the barrier that defines the folding pathway. Unfortunately, however, unlike those of the initial, unfolded state and the final, folded state of the protein, the structure in the transition state cannot be directly assessed via experiment. Instead, experimentalists infer the structure of the transition state, often by estimating changes in its free energy by measuring the effects of amino acid substitutions on folding and unfolding rates (Phi-value analysis). In this article we show how to obtain more efficient estimates of these important quantities via improved experimental designs, and how to avoid common pitfalls in the analysis of kinetic data during the extraction of these parameters.     

A catalytic loop within Pseudomonas aeruginosa exotoxin A modulates its transferase activity.

Mutagenesis techniques were used to replace two loop regions within the catalytic domain of Pseudomonas aeruginosa exotoxin A (ETA) with functionally silent polyglycine loops. The loop mutant proteins, designated polyglycine Loops N and C, were both less active than the wild-type enzyme. However, the polyglycine Loop C mutant protein, replaced with the Gly(483)-Gly(490) loop, showed a much greater loss of enzymatic activity than the polyglycine Loop N protein. The former mutant enzyme exhibited an 18,000-fold decrease in catalytic turnover number (k(cat)), with only a marginal effect on the K(m) value for NAD(+) and the eukaryotic elongation factor-2 binding constant. Furthermore, alanine-scanning mutagenesis of this active-site loop region revealed the specific pattern of a critical region for enzymatic activity. Binding and kinetic data suggest that this loop modulates the transferase activity between ETA and eukaryotic elongation factor-2 and may be responsible for stabilization of the transition state for the reaction. Sequence alignment and molecular modeling also identified a similar loop within diphtheria toxin, a functionally and structurally related class A-B toxin. Based on these results and the similarities between ETA and diphtheria toxin, we propose that this catalytic subregion represents the first report of a diphthamide-specific ribosyltransferase structural motif. We expect these findings to further the development of pharmaceuticals designed to prevent ETA toxicity by disrupting the stabilization of the transition state during the ADP-ribose transfer event.     

Calculation of mutational free energy changes in transition states for protein folding

Recent advances in experimental and computational methods have made it possible to determine with considerable accuracy the structures whose formation is rate limiting for the folding of some small proteins-the transition state ensemble, or TSE. We present a method to analyze and validate all-atom models of such structures. The method is based on the comparison of experimental data with the computation of the change in free energy of the TSE resulting from specific mutations. Each mutation is modeled individually in all members of an ensemble of transition state structures using a method originally developed to predict mutational changes in the stability of native proteins. We first apply this method to six proteins for which we have determined the TSEs with a technique that uses experimental mutational data (Phi-values) as restraints in the structure determination and find a highly significant correlation between the calculated free energy changes and those derived from experimental kinetic data. We then use the procedure to analyze transition state structures determined by molecular dynamics simulations of unfolding, again finding a high correlation. Finally, we use the method to estimate changes in folding rates of several hydrophobic core mutants of Fyn SH3. Taken together, these results show that the procedure developed here is a tool of general validity for analyzing, assessing, and improving the quality of the structures of transition states for protein folding.     

The Dominant Folding Route Minimizes Backbone Distortion in SH3

Energetic frustration in protein folding is minimized by evolution to create a smooth and robust energy landscape. As a result the geometry of the native structure provides key constraints that shape protein folding mechanisms. Chain connectivity in particular has been identified as an essential component for realistic behavior of protein folding models. We study the quantitative balance of energetic and geometrical influences on the folding of SH3 in a structure-based model with minimal energetic frustration. A decomposition of the two-dimensional free energy landscape for the folding reaction into relevant energy and entropy contributions reveals that the entropy of the chain is not responsible for the folding mechanism. Instead the preferred folding route through the transition state arises from a cooperative energetic effect. Off-pathway structures are penalized by excess distortion in local backbone configurations and contact pair distances. This energy cost is a new ingredient in the malleable balance of interactions that controls the choice of routes during protein folding.     

Geometrical Frustration in Interleukin-33 Decouples the Dynamics of the Functional Element from the Folding Transition State Ensemble

Interleukin-33 (IL-33) is currently the focus of multiple investigations into targeting pernicious inflammatory disorders. This mediator of inflammation plays a prevalent role in chronic disorders such as asthma, rheumatoid arthritis, and progressive heart disease. In order to better understand the possible link between the folding free energy landscape and functional regions in IL-33, a combined experimental and theoretical approach was applied. IL-33 is a pseudo- symmetrical protein composed of three distinct structural elements that complicate the folding mechanism due to competition for nucleation on the dominant folding route. Trefoil 1 constitutes the majority of the binding interface with the receptor whereas Trefoils 2 and 3 provide the stable scaffold to anchor Trefoil 1. We identified that IL-33 folds with a three-state mechanism, leading to a rollover in the refolding arm of its chevron plots in strongly native conditions. In addition, there is a second slower refolding phase that exhibits the same rollover suggesting similar limitations in folding along parallel routes. Characterization of the intermediate state and the rate limiting steps required for folding suggests that the rollover is attributable to a moving transition state, shifting from a post- to pre-intermediate transition state as you move from strongly native conditions to the midpoint of the transition. On a structural level, we found that initially, all independent Trefoil units fold equally well until a Q_CA of 0.35 when Trefoil 1 will backtrack in order to allow Trefoils 2 and 3 to fold in the intermediate state, creating a stable scaffold for Trefoil 1 to fold onto during the final folding transition. The formation of this intermediate state and subsequent moving transition state is a result of balancing the difficulty in folding the functionally important Trefoil 1 onto the remainder of the protein. Taken together our results indicate that the functional element of the protein is geometrically frustrated, requiring the more stable elements to fold first, acting as a scaffold for docking of the functional element to allow productive folding to the native state.