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.     

On the Prediction of Folding Nuclei in Globular Proteins

An approach to predicting folding nuclei in globular proteins with known three-dimensional structures is proposed. This approach is based on the pinpointing of the lowest saddle points on the barrier between the unfolded state and native structure on the free-energy landscape of a protein chain; the proposed technique uses the dynamic programming method. A comparison of calculation results with experimental data on the folding nuclei of 21 proteins shows that the model provides good Φ value predictions for protein structures determined by X-ray analysis and, less successfully, in structures determined by nuclear magnetic resonance. Consideration of the whole ensemble of transition states provides a better prediction of folding nuclei than consideration of only transition states with lowest free energies. In addition, we predict the location of folding nuclei in three-dimensional structures of some proteins whose folding kinetics is being studied, but there is no experimental evidence concerning their folding nuclei.     

Sensitivity of the folding pathway to the details of amino acid sequence

Experimental data on the structure of transition state demonstrate that the proteins with the same topology as a rule have similar folding nuclei (the structured formed part of the transition state). In this review we discuss the experimental works which show that the position of folding nuclei is different among proteins with the same topology. These facts underline that the folding pathway is sensitive to the details of amino-acid sequence.     

Structural analysis of the rate-limiting transition states in the folding of Im7 and Im9: similarities and differences in the folding of homologous proteins

The bacterial immunity proteins Im7 and Im9 fold with mechanisms of different kinetic complexity. Whilst Im9 folds in a two-state transition at pH 7.0 and 10 degrees C, Im7 populates an on-pathway intermediate under these conditions. In order to assess the role of sequence versus topology in the folding of these proteins, and to analyse the effect of populating an intermediate on the landscape for folding, we have determined the conformational properties of the rate-limiting transition state for Im9 folding/unfolding using Phi(F)-value analysis and have compared the results with similar data obtained previously for Im7. The data show that the rate-limiting transition states for Im9 and Im7 folding/unfolding are similar: both are compact (beta(T)=0.94 and 0.89, respectively) and contain three of the four native helices docked around a specific hydrophobic core. Significant differences are observed, however, in the magnitude of the Phi(F)-values obtained for the two proteins. Of the 20 residues studied in both proteins, ten have Phi(F)-values in Im7 that exceed those in Im9 by more than 0.2, and of these five differ by more than 0.4. The data suggest that the population of an intermediate in Im7 results in folding via a transition state ensemble that is conformationally restricted relative to that of Im9. The data are consistent with the view that topology is an important determinant of folding. Importantly, however, they also demonstrate that while the folding transition state may be conserved in homologous proteins that fold with two and three-state kinetics, the population of an intermediate can have a significant effect on the breadth of the transition state ensemble.     

Identifying folding nucleus based on residue contact networks of proteins

In the native structure of a protein, all the residues are tightly parked together in a specific order following its folding and every residue contacts with some spatially neighbor residues. A residue contact network can be constructed by defining the residues as nodes and the native contacts as edges. During the folding of small single-domain proteins, there is a set of contacts (or bonds), defined as the folding nucleus (FN), which is formed around the transition state, i.e., a rate-limiting barrier located at about the middle between the unfolded states and the native state on the free energy landscape. Such a FN plays an essential role in the folding dynamics and the residues, which form the related contacts called as folding nucleus residues (FNRs). In this work, the FNRs in proteins are identified by using quantities which characterize the topology of residue contact networks of proteins. By comparing the specificities of residues with the network quantities K(R), L(R), and D(R), up to 90% FNRs of six typical proteins found experimentally are identified. It is found that the FNRs behave the full-closeness centrals rather than degree or closeness centers in the residue contact network, implying that they are important to the folding cooperativity of proteins. Our study shows that the FNRs can be identified solely from the native structures of proteins based on the analysis of residue contact network without any knowledge of the transition state ensemble.     

Outlining folding nuclei in globular proteins

Our theoretical approach for prediction of folding/unfolding nuclei in three-dimensional protein structures is based on a search for free energy saddle points on networks of protein unfolding pathways. Under some approximations, this search is performed rapidly by dynamic programming and results in prediction of Phi values, which can be compared with those found experimentally. In this study, we optimize some details of the model (specifically, hydrogen atoms are taken into account in addition to heavy atoms), and compare the theoretically obtained and experimental Phi values (which characterize involvement of residues in folding nuclei) for all 17 proteins, where Phi values are now known for many residues. We show that the model provides good Phi value predictions for proteins whose structures have been determined by X-ray analysis (the average correlation coefficient is 0.65), with a more limited success for proteins whose structures have been determined by NMR techniques only (the average correlation coefficient is 0.34), and that the transition state free energies computed from the same model are in a good anticorrelation with logarithms of experimentally measured folding rates at mid-transition (the correlation coefficient is -0.73).     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

Critical nucleation size in the folding of small apparently two-state proteins

For apparently two-state proteins, we found that the size (number of folded residues) of a transition state is mostly encoded by the topology, defined by total contact distance (TCD) of the native state, and correlates with its folding rate. This is demonstrated by using a simple procedure to reduce the native structures of the 41 two-state proteins with native TCD as a constraint, and is further supported by analyzing the results of eight proteins from protein engineering studies. These results support the hypothesis that the major rate-limiting process in the folding of small apparently two-state proteins is the search for a critical number of residues with the topology close to that of the native state.     

Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

We perform a detailed analysis of the thermodynamics and folding kinetics of the SH3 domain fold with discrete molecular dynamic simulations. We propose a protein model that reproduces some of the experimentally observed thermodynamic and folding kinetic properties of proteins. Specifically, we use our model to study the transition state ensemble of the SH3 fold family of proteins, a set of unstable conformations that fold to the protein native state with probability 1/2. We analyze the participation of each secondary structure element formed at the transition state ensemble. We also identify the folding nucleus of the SH3 fold and test extensively its importance for folding kinetics. We predict that a set of amino acid contacts between the RT-loop and the distal hairpin are the critical folding nucleus of the SH3 fold and propose a hypothesis that explains this result.     

Transition state contact orders correlate with protein folding rates

We have used molecular dynamics simulations restrained by experimental phi values derived from protein engineering experiments to determine the structures of the transition state ensembles of ten proteins that fold with two-state kinetics. For each of these proteins we then calculated the average contact order in the transition state ensemble and compared it with the corresponding experimental folding rate. The resulting correlation coefficient is similar to that computed for the contact orders of the native structures, supporting the use of native state contact orders for predicting folding rates. The native contacts in the transition state also correlate with those of the native state but are found to be about 30% lower. These results show that, despite the high levels of heterogeneity in the transition state ensemble, the large majority of contributing structures have native-like topologies and that the native state contact order captures this phenomenon.     

On the role of some conserved and nonconserved amino acid residues in the transitional state and intermediate of apomyoglobin folding

The contributions of some amino acid residues in the A, B, G, and H helices to the formation of the folding nucleus and folding intermediate of apomyoglobin were estimated. The effects of point substitutions of Ala for hydrophobic amino acid residues on the structural stability of the native (N) protein and its folding intermediate (I), as well as on the folding/unfolding rates for four mutant apomyoglobin forms, were studied. The equilibrium and kinetic studies of the folding/unfolding rates of these mutant proteins in a wide range of urea concentrations demonstrated that their native state was considerably destabilized as compared with the wild-type protein, whereas the stability of the intermediate state changed moderately. It was shown that the amino acid residues in the A, G, and H helices contributed insignificantly to the stabilization of the apomyoglobin folding nucleus in the rate-limiting I ? N transition, taking place after the formation of the intermediate, whereas the residue of the B helix was of great importance in the formation of the folding nucleus in this transition.     

The origins of asymmetry in the folding transition states of protein L and protein G

Topology has been shown to be an important determinant of many features of protein folding; however, the delineation of sequence effects on folding remains obscure. Furthermore, differentiation between the two influences proves difficult due to their intimate relationship. To investigate the effect of sequence in the absence of significant topological differences, we examined the folding mechanisms of segment B1 peptostreptococcal protein L and segment B1 of streptococcal protein G. These proteins share the same highly symmetrical topology. Despite this symmetry, neither protein folds through a symmetrical transition state. We analyzed the origins of this difference using theoretical models. We found that the strength of the interactions present in the N-terminal hairpin of protein L causes this hairpin to form ahead of the C-terminal hairpin. The difference in chain entropy associated with the formation of the hairpins of protein G proves sufficient to beget initiation of folding at the shorter C-terminal hairpin. Our findings suggest that the mechanism of folding may be understood by examination of the free energy associated with the formation of partially folded microstates.     

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.     

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.     

Mutational analysis of acylphosphatase suggests the importance of topology and contact order in protein folding.

Muscle acylphosphatase (AcP) is a small protein that folds very slowly with two-state behavior. The conformational stability and the rates of folding and unfolding have been determined for a number of mutants of AcP in order to characterize the structure of the folding transition state. The results show that the transition state is an expanded version of the native protein, where most of the native interactions are partially established. The transition state of AcP turns out to be remarkably similar in structure to that of the activation domain of procarboxypeptidase A2 (ADA2h), a protein having the same overall topology but sharing only 13% sequence identity with AcP. This suggests that transition states are conserved between proteins with the same native fold. Comparison of the rates of folding of AcP and four other proteins with the same topology, including ADA2h, supports the concept that the average distance in sequence between interacting residues (that is, the contact order) is an important determinant of the rate of protein folding.     

Hammond behavior versus ground state effects in protein folding: evidence for narrow free energy barriers and residual structure in unfolded states

Apparent transition state movement upon mutation or changes in solvent conditions is frequently observed in protein folding and is often interpreted in terms of Hammond behavior. This led to the conclusion that barrier regions in protein folding are broad maxima on the free energy landscape. Here, we use the concept of self-interaction and cross-interaction parameters to test experimental data of 21 well-characterized proteins for Hammond behavior. This allows us to characterize the origin of transition state movements along different reaction coordinates. Only one of the 21 proteins shows a small but coherent transition state movement in agreement with the Hammond postulate. In most proteins the structure of the transition state is insensitive to changes in protein stability. The apparent change in the position of the transition state upon mutation, which is frequently observed in phi-value analysis, is in most cases due to ground-state effects caused by structural changes in the unfolded state. This argues for significant residual structure in unfolded polypeptide chains of many proteins. Disruption of these residual interactions by mutation often leads to decreased folding rates, which implies that these interactions are still present in the transition state. The failure to detect Hammond behavior shows that the free energy barriers encountered by a folding polypeptide chain are generally rather narrow and robust maxima for all experimentally explorable reaction coordinates.     

Prediction of folding transition-state position (betaT) of small, two-state proteins from local secondary structure content

Folding kinetics of proteins is governed by the free energy and position of transition states. But attempts to predict the position of folding transition state on reaction pathway from protein structure have been met with only limited success, unlike the folding-rate prediction. Here, we find that the folding transition-state position is related to the secondary structure content of native two-state proteins. We present a simple method for predicting the transition-state position from their alpha-helix, turn and polyproline secondary structures. The method achieves 81% correlation with experiment over 24 small, two-state proteins, suggesting that the local secondary structure content, especially for content of alpha-helix, is a determinant of the solvent accessibility of the transition state ensemble and size of folding nucleus.     

Prediction of folding rates and transition-state placement from native-state geometry

A variety of experimental and theoretical studies have established that the folding process of monomeric proteins is strongly influenced by the topology of the native state. In particular, folding times have been shown to correlate well with the contact order, a measure of contact locality. Our investigation focuses on identifying additional topologic properties that correlate with experimentally measurable quantities, such as folding rates and transition-state placement, for both two- and three-state folders. The validation against data from 40 experiments shows that a particular topological property that measures the interdependence of contacts, termed cliquishness or clustering coefficient, can account with statistically significant accuracy both for the transition state placement and especially for folding rates. The observed correlations can be further improved by optimally combining the distinct topological information captured by cliquishness and contact order.     

Nucleation‐based prediction of the protein folding rate and its correlation with the folding nucleus size

The problem of protein self‐organization is in the focus of current molecular biology studies. Although the general principles are understood, many details remain unclear. Specifically, protein folding rates are of interest because they dictate the rate of protein aggregation which underlies many human diseases. Here we offer predictions of protein folding rates and their correlation with folding nucleus sizes. We calculated free energies of the transition state and sizes of folding nuclei for 84 proteins and peptides whose other parameters were measured at the point of thermodynamic equilibrium between their unfolded and native states. We used the dynamic programming method where each residue was considered to be either as folded as in its native state or completely disordered. The calculated and measured folding rates showed a good correlation at the temperature mid‐transition point (the correlation coefficient was 0.75). Also, we pioneered in demonstrating a moderate (‐0.57) correlation coefficient between the calculated sizes of folding nuclei and the folding rates. Predictions made by different methods were compared. The established good correlation between the estimated free energy barrier and the experimentally found folding rate of each studied protein/peptide indicates that our model gives reliable results for the considered data set. Proteins 2012; © 2012 Wiley Periodicals, Inc.