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.     

Osmolyte-induced protein folding free energy changes

Upon addition of protecting osmolyte to an aqueous solution of an intrinsically unstructured protein, spectral observables are often seen to change in a sigmoid fashion as a function of increasing osmolyte concentration. Commonly, such data are analyzed using the linear extrapolation model (LEM), a method that defines a scale from 0%-100% folded species at each osmolyte concentration by means of extending pre- and post-folding baselines into the transition region. Defining the 0%-100% folding scale correctly for each osmolyte is an important part of the analysis, leading to evaluation of the fraction of folded protein existing in the absence of osmolytes. In this study, we used reduced and carboxyamidated RNase T1 (RCAM-T1) as an intrinsically unstructured protein, and determined the thermodynamic stability of RCAM-T1 induced by naturally occurring osmolytes. Because the folded fraction of the protein population determined by experiments of thermal and urea-induced denaturation is nonzero in the absence of osmolytes at 15 degrees C, the commonly used LEM can lead to false values of DeltaG[stackD-->N0] for protein folding due to the arbitrary assumption that the protein is 100% unfolded in the presence of buffer alone. To correct this problem, titration of the protein solution with urea and extrapolating back to zero urea concentration gives the spectral value for 100% denatured protein. With fluorescence as the observable we redefine F/F0 to F/F0extrap = 1.0 and require that the denatured-state baseline have this value as its intercept. By so doing, the 0%-100% scale-corrected DeltaG[D-->N0] values of RCAM-T1 folding in the presence of various osmolytes are then found to be identical, with small error, demonstrating that DeltaG[D-->N0] is independent of the osmolytes used. Such a finding is an important step in validating this quantity derived from the LEM as having the properties expected of an authentic thermodynamic parameter. The rank order of osmolyte efficacies in stabilizing RCAM-T1 is sarcosine > sucrose > sorbitol > proline > betaine > glycerol.     

Exploring folding free energy landscapes using computational protein design

Recent advances in computational protein design have allowed exciting new insights into the sequence dependence of protein folding free energy landscapes. Whereas most previous studies have examined the sequence dependence of protein stability and folding kinetics by characterizing naturally occurring proteins and variants of these proteins that contain a small number of mutations, it is now possible to generate and characterize computationally designed proteins that differ significantly from naturally occurring proteins in sequence and/or structure. These computer-generated proteins provide insights into the determinants of protein structure, stability and folding, and make it possible to disentangle the properties of proteins that are the consequence of natural selection from those that reflect the fundamental physical chemistry of polypeptide chains.     

Reliable protein folding on complex energy landscapes: the free energy reaction path

A theoretical framework is developed to study the dynamics of protein folding. The key insight is that the search for the native protein conformation is influenced by the rate r at which external parameters, such as temperature, chemical denaturant, or pH, are adjusted to induce folding. A theory based on this insight predicts that 1), proteins with complex energy landscapes can fold reliably to their native state; 2), reliable folding can occur as an equilibrium or out-of-equilibrium process; and 3), reliable folding only occurs when the rate r is below a limiting value, which can be calculated from measurements of the free energy. We test these predictions against numerical simulations of model proteins with a single energy scale.     

Non-linear rate-equilibrium free energy relationships and Hammond behavior in protein folding

Non-linear rate-equilibrium relationships upon mutation or changes in solvent conditions are frequently observed in protein folding reactions and are usually interpreted in terms of Hammond behavior. Here we first give a general overview over the concept of transition state movements in chemical reactions and discuss its application to protein folding. We then show examples for genuine Hammond behavior and for apparent transition state movements caused by other effects like changes in the rate-limiting step of the folding reaction or ground state effects, i.e. structural changes in either the native state or the unfolded state. These examples show that apparent transition state movements can easily be mistaken for Hammond behavior. We describe experimental tests using self- and cross-interaction parameters to distinguish between structural changes in a single transition state following Hammond behavior and apparent transition state movements caused by other effects.     

Protein folding: the free energy surface

Quantitative models and experiments are revealing how the folding free energy surface of a protein is sculpted by sequence and environment. The sometimes conflicting demands of folding, structure and function determine which folding pathways, if any, dominate. Recent advances include experimental estimates of diffusive barrier-crossing times, the observation of ultrafast folders amenable to full-atom simulation, the use of thermodynamic tuning and nonconservative mutations to probe 'hidden' parts of the free energy surface, and a complete microscopic theory of folding.     

Calculation of the free energy of association for protein complexes.

We have developed a method for calculating the association energy of quaternary complexes starting from their atomic coordinates. The association energy is described as the sum of two solvation terms and an energy term to account for the loss of translational and rotational entropy. The calculated solvation energy, using atomic solvation parameters and the solvent accessible surface areas, has a correlation of 96% with experimentally determined values. We have applied this methodology to examine intermediates in viral assembly and to assess the contribution isomerization makes to the association energy of molecular complexes. In addition, we have shown that the calculated association can be used as a predictive tool for analyzing modeled molecular complexes.     

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.     

Dynamics and cooperativity of Trp-cage folding

In this study, multiple independent molecular dynamics (MD) simulations on Trp-cage folding were performed at 300, 325 and 375 K using generalized Born (GB) implicit solvent model. The orientational movement of the side-chain of Trp6 to form a hydrophobic core with 3₁₀-helix was observed. The breaking/formation of a salt bridge between Asp9 and Arg16 was proposed to be the prerequisite for Trp-cage folding/refolding. Our results demonstrate that the cooperation between the salt bridge and the Trp6 orientation leads to a stable tertiary structure of Trp-cage. Analyses on backbone concerted motions at different temperatures indicate that interactions between Trp6 and 3₁₀-helix & Pro18 and between Pro12 and Pro17 & Pro18 are weakened at 375 K but strengthened at lower temperatures, suggesting that they could be the potential driving force of hydrophobic collapse.     

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

We discuss the formulation of free energy functionals that describe the formation of structure in partially folded proteins. These free energy functionals take into account the inhomogeneous nature of contact energies, chain entropy and cooperative contributions reflecting the many body character of some folding forces like hydrophobicity, but do not directly account for non-native contacts because they assume the validity of the minimal frustration principle. We show how the free energy functionals can be used to interpret experiments on partially folded proteins that probe the fractional occupancy of specific local structures. In particular, we study the hydrogen protection factors in lysozyme studied in transient experiments by Gladwin and Evans and by Nash and Jonas using equilibrium pressure denaturation and the NMR order parameters measured by Dobson and Kim for the homologous protein alpha-lactalbumin.     

Molecular simulations of cotranslational protein folding: fragment stabilities, folding cooperativity, and trapping in the ribosome

Although molecular simulation methods have yielded valuable insights into mechanistic aspects of protein refolding in vitro, they have up to now not been used to model the folding of proteins as they are actually synthesized by the ribosome. To address this issue, we report here simulation studies of three model proteins: chymotrypsin inhibitor 2 (CI2), barnase, and Semliki forest virus protein (SFVP), and directly compare their folding during ribosome-mediated synthesis with their refolding from random, denatured conformations. To calibrate the methodology, simulations are first compared with in vitro data on the folding stabilities of N-terminal fragments of CI2 and barnase; the simulations reproduce the fact that both the stability and thermal folding cooperativity increase as fragments increase in length. Coupled simulations of synthesis and folding for the same two proteins are then described, showing that both fold essentially post-translationally, with mechanisms effectively identical to those for refolding. In both cases, confinement of the nascent polypeptide chain within the ribosome tunnel does not appear to promote significant formation of native structure during synthesis; there are however clear indications that the formation of structure within the nascent chain is sensitive to location within the ribosome tunnel, being subject to both gain and loss as the chain lengthens. Interestingly, simulations in which CI2 is artificially stabilized show a pronounced tendency to become trapped within the tunnel in partially folded conformations: non-cooperative folding, therefore, appears in the simulations to exert a detrimental effect on the rate at which fully folded conformations are formed. Finally, simulations of the two-domain protease module of SFVP, which experimentally folds cotranslationally, indicate that for multi-domain proteins, ribosome-mediated folding may follow different pathways from those taken during refolding. Taken together, these studies provide a first step toward developing more realistic methods for simulating protein folding as it occurs in vivo.     

Network and graph analyses of folding free energy surfaces

Protein folding is governed by a complex free energy surface whose entropic contributions are relevant because of the large number of degrees of freedom involved. Such complexity, in particular the conformational heterogeneity of the denatured state, is hidden in projections onto one or two order parameters (e.g. fraction of native contacts and/or radius of gyration), which usually results in relatively smooth surfaces. Recent approaches borrowed from network and graph theory have yielded quantitative unprojected representations of the free energy surfaces of a beta-hairpin and a three-stranded beta-sheet peptide using equilibrium folding-unfolding molecular dynamics simulations. Interestingly, the network and graph analyses of these structured peptides have revealed a very heterogeneous denatured state ensemble. It includes high-enthalpy, high-entropy conformations with fluctuating non-native secondary structure, as well as low-enthalpy, low-entropy traps.     

On the origin of the cooperativity of protein folding: Implications from model simulations

There is considerable experimental evidence that the cooperativity of protein folding resides in the transition from the molten globule to the native state. The objective of this study is to examine whether simplified models can reproduce this cooperativity and if so, to identify its origin. In particular, the thermodynamics of the conformational transition of a previously designed sequence (A. Kolinski, W. Galazka, and J. Skolnick, J. Chem. Phys. 103: 10286–10297, 1995), which adopts a very stable Greek-key β-barrel fold has been investigated using the entropy Monte Carlo sampling (ESMC) technique of Hao and Scheraga (M.-H. Hao and H.A. Scheraga, J. Phys. Chem. 98: 9882–9883, 1994). Here, in addition to the original potential, which includes one body and pair interactions between side chains, the force field has been supplemented by two types of multi-body potentials describing side chain interactions. These potentials facilitate the proteinlike pattern of side chain packing and consequently increase the cooperativity of the folding process. Those models that include an explicit cooperative side chain packing term exhibit a well-defined all-or-none transition from a denatured, random coil state to a high-density, well-defined, nativelike low-energy state. By contrast, models lacking such a term exhibit a conformational transition that is essentially continuous. Finally, an examination of the conformations at the free-energy barrier between the native and denatured states reveals that they contain a substantial amount of native-state secondary structure, about 50% of the native contacts, and have an average root mean square radius of gyration that is about 15% larger than native. © 1996 Wiley-Liss, Inc.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

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.     

The nature of the free energy barriers to two-state folding

We introduce a simple procedure to analyze the temperature dependence of the folding and unfolding rates of two-state proteins. We start from the simple transition-state-like rate expression: k = D(eff)exp(-DeltaG(TS)/RT), in which upper and lower bounds for the intra-chain effective diffusion coefficient (D(eff)) are obtained empirically using the timescales of elementary processes in protein folding. From the changes in DeltaG(TS) as a function of temperature, we calculate enthalpies and heat capacities of activation, together with the more elusive entropies of activation. We then estimate the conformational entropy of the transition state by extrapolation to the temperature at which the solvation entropy vanishes by cancellation between polar and apolar terms. This approach is based on the convergence temperatures for the entropy of solvating apolar (approximately 385 K) and polar groups (approximately 335 K), the assumption that the structural properties of the transition state are somewhere in between the unfolded and folded states, and the established relationship between observed heat capacity and solvent accessibility.1 To circumvent the lack of structural information about transition states, we use the empirically determined heat capacities of activation as constraints to identify the extreme values of the transition state conformational entropy that are consistent with experiment. The application of this simple approach to six two-state folding proteins for which there is temperature-dependent data available in the literature provides important clues about protein folding. For these six proteins, we obtain an average equilibrium cost in conformational entropy of -4.3 cal x mol(-1)K(-1)per residue, which is in close agreement to previous empirical and computational estimates of the same quantity. Furthermore, we find that all these proteins have a conformationally diverse transition state, with more than half of the conformational entropy of the unfolded state. In agreement with predictions from theory and computer simulations, the transition state signals the change from a regime dominated by loss in conformational entropy to one driven by the gain in stabilization free energy (i.e., including protein interactions and solvation effects). Moreover, the height of the barrier is determined by how much stabilization free energy is realized at that point, which is related to the relative contribution of local versus non-local interactions. A remarkable observation is that the fraction of conformational entropy per residue that is present in the transition state is very similar for the six proteins in this study. Based on this commonality, we propose that the observed change in thermodynamic regime is connected to a change in the pattern of structure formation: from one driven by formation of pairwise interactions to one dominated by coupling of the networks of interactions involved in forming the protein core. In this framework, the barrier to two-state folding is crossed when the folding protein reaches a "critical native density" that allows expulsion of remaining interstitial water and consolidation of the core. The principle of critical native density should be general for all two-state proteins, but can accommodate different folding mechanisms depending on the particularities of the structure and sequence.     

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq-UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq-UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq-UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq-UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq-UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq-UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq-UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq-UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq-UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq-UIM, appears to behave very close to an equilibrium two-state.