The folding transition state of the cold shock protein is strongly polarized

It has been shown recently that an 11-residue peptide fragment of transthyretin, TTR(105-115), can form amyloid fibrils in vitro by adopting an extended beta-strand conformation. We used molecular dynamics simulations on systems of TTR(105-115) peptides, for a total length of about 5 micros, to explore the process of self-assembly and the structures of the resulting aggregates. Our results suggest that an antiparallel association of the beta-strands is more probable than a parallel one and that the central residues (T106-L111) in a beta-strand have a high propensity to form inter-peptide hydrogen bonds. The study of the dynamics of self-association indicated that, for this peptide, trajectories leading to conformations with high alpha-helical content are off-pathway from those leading to aggregates with high beta-structure content. We also show that the diverse oligomeric structures that form spontaneously in the molecular dynamics simulations are, to a large extent, compatible with solid-state NMR experimental measurements, including chemical shifts, on fully formed fibrils. The strategy that we present may therefore be used in the design of new experiments to determine the structure of amyloid fibrils, such as those involving site-specific isotope labelling schemes to measure key inter-atomic distances.     

The role of Trp-82 in the folding of intestinal fatty acid binding protein

Multiple phases have been observed during the folding and unfolding of intestinal fatty acid binding protein (WT-IFABP) by stopped-flow fluorescence. Site-directed mutagenesis has been used to examine the role of each of the two tryptophans of this protein in these processes. The unfolding and refolding kinetics of the mutant protein containing only tryptophan 82 (W6Y-IFABP) showed that the tryptophan at this location was critical to the fluorescence signal changes observed throughout the unfolding reaction and early in the refolding reaction. However, the kinetic patterns of the mutant protein containing only tryptophan 6 (W82Y-IFABP) indicated that the tryptophan at this location participated in the fluorescence signal changes observed early in the unfolding reaction and late in the refolding reaction. Together, these data suggest that native-like structure was formed first in the vicinity of tryptophan 82, near the center of the hydrophobic core of this beta-sheet protein, prior to formation of native-like structure in the periphery of the protein.     

Destabilization of the Escherichia coli RNase H kinetic intermediate: switching between a two-state and three-state folding mechanism

Escherichia coli RNase H folds through a partially folded kinetic intermediate that mirrors a rarely populated, partially unfolded form detectable by native-state hydrogen exchange under equilibrium conditions. Residue 53 is at the interface of two helices known to be structured in this intermediate. Kinetic refolding studies on mutant proteins varying in size and hydrophobicity at residue 53 support a contribution of hydrophobicity to the stabilities of the kinetic intermediate and the transition state. Packing interactions also play a significant role in the stability of these two states, though they play a much larger role in the native-state stability. One dramatic mutation, I53D, results in the conversion from a three-state to a two-state folding mechanism, which is explained most easily through a simple destabilization of the kinetic intermediate such that it is no longer stable with respect to the unfolded state. These results demonstrate that interactions that stabilize an intermediate can accelerate folding if these same interactions are present in the transition state. Our results are consistent with a hierarchical model of folding, where the intermediate consists of native-like interactions, is on-pathway, and is productive for folding.     

Cooperativity in two-state protein folding kinetics

We present a solvable model that predicts the folding kinetics of two-state proteins from their native structures. The model is based on conditional chain entropies. It assumes that folding processes are dominated by small-loop closure events that can be inferred from native structures. For CI2, the src SH3 domain, TNfn3, and protein L, the model reproduces two-state kinetics, and it predicts well the average Phi-values for secondary structures. The barrier to folding is the formation of predominantly local structures such as helices and hairpins, which are needed to bring nonlocal pairs of amino acids into contact.     

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.     

Role of electrostatic interactions for the stability and folding behavior of cold shock protein

The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Predicting repeat protein folding kinetics from an experimentally determined folding energy landscape

The Notch ankyrin domain is a repeat protein whose folding has been characterized through equilibrium and kinetic measurements. In previous work, equilibrium folding free energies of truncated constructs were used to generate an experimentally determined folding energy landscape (Mello and Barrick, Proc Natl Acad Sci USA 2004;101:14102–14107). Here, this folding energy landscape is used to parameterize a kinetic model in which local transition probabilities between partly folded states are based on energy values from the landscape. The landscape‐based model correctly predicts highly diverse experimentally determined folding kinetics of the Notch ankyrin domain and sequence variants. These predictions include monophasic folding and biphasic unfolding, curvature in the unfolding limb of the chevron plot, population of a transient unfolding intermediate, relative folding rates of 19 variants spanning three orders of magnitude, and a change in the folding pathway that results from C‐terminal stabilization. These findings indicate that the folding pathway(s) of the Notch ankyrin domain are thermodynamically selected: the primary determinants of kinetic behavior can be simply deduced from the local stability of individual repeats.     

Local interactions and the optimization of protein folding

The role of local interactions in protein folding has recently been the subject of some controversy. Here we investigate an extension of Zwanzig's simple and general model of folding in which local and nonlocal interactions are represented by functions of single and multiple conformational degrees of freedom, respectively. The kinetics and thermodynamics of folding are studied for a series of energy functions in which the energy of the native structure is fixed, but the relative contributions of local and nonlocal interactions to this energy are varied over a broad range. For funnel shaped energy landscapes, we find that 1) the rate of folding increases, but the stability of the folded state decreases, as the contribution of local interactions to the energy of the native structure increases, and 2) the amount of native structure in the unfolded state and the transition state vary considerably with the local interaction strength. Simple exponential kinetics and a well-defined free energy barrier separating folded and unfolded states are observed when nonlocal interactions make an appreciable contribution to the energy of the native structure; in such cases a transition state theory type approximation yields reasonably accurate estimates of the folding rate. Bumps in the folding funnel near the native state, which could result from desolvation effects, side chain freezing, or the breaking of nonnative contacts, significantly alter the dependence of the folding rate on the local interaction strength: the rate of folding decreases when the local interaction strength is increased beyond a certain point. A survey of the distribution of strong contacts in the protein structure database suggests that evolutionary optimization has involved both kinetics and thermodynamics: strong contacts are enriched at both very short and very long sequence separations. Proteins 29:282–291, 1997. © 1997 Wiley-Liss, Inc.     

Three-state protein folding: experimental determination of free-energy profile

When considering protein folding with a transient intermediate, a difficulty arises as to determination of the rates of separate transitions. Here we overcome this problem, using the kinetic studies of the unfolding/refolding reactions of the three-state protein apomyoglobin as a model. Amplitudes of the protein refolding kinetic burst phase corresponding to the transition from the unfolded (U) to intermediate (I) state, that occurs prior to the native state (N) formation, allow us to estimate relative populations of the rapidly converting states at various final urea concentrations. On the basis of these proportions, a complicated experimental chevron plot has been deconvolved into the urea-dependent rates of the I<-->N and U<-->N transitions to give the dependence of free energies of the main transition state and of all three (N, I, and U) stable states on urea concentration.     

ALASKA: A Mathematica package for two-state kinetic analysis of protein folding reactions

A Mathematica package (ALASKA) has been developed to simplify the measurement of protein folding kinetics by analysis of ¹H NMR lineshape analysis. This package reads NMR data in ASCII format and can simulate an aromatic ¹ NMR spectrum with or without lineshape broadening from chemical exchange. We describe the analysis of a urea denaturation series of a fast-folding protein, the G46A/G48A variant of monomeric repressor.     

Sub-microsecond protein folding

We have investigated the structure, equilibria, and folding kinetics of an engineered 35-residue subdomain of the chicken villin headpiece, an ultrafast-folding protein. Substitution of two buried lysine residues by norleucine residues stabilizes the protein by 1 kcal/mol and increases the folding rate sixfold, as measured by nanosecond laser T-jump. The folding rate at 300 K is (0.7 micros)(-1) with little or no temperature dependence, making this protein the first sub-microsecond folder, with a rate only twofold slower than the theoretically predicted speed limit. Using the 70 ns process to obtain the effective diffusion coefficient, the free energy barrier height is estimated from Kramers theory to be less than approximately 1 kcal/mol. X-ray crystallographic determination at 1A resolution shows no significant change in structure compared to the single-norleucine-substituted molecule and suggests that the increased stability is electrostatic in origin. The ultrafast folding rate, very accurate X-ray structure, and small size make this engineered villin subdomain an ideal system for simulation by atomistic molecular dynamics with explicit solvent.     

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.     

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.     

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.     

Early collapse is not an obligate step in protein folding

The dimensions and secondary structure content of two proteins which fold in a two-state manner are measured within milliseconds of denaturant dilution using synchrotron-based, stopped-flow small-angle X-ray scattering and far-UV circular dichroism spectroscopy. Even upon a jump to strongly native conditions, neither ubiquitin nor common-type acylphosphatase contract prior to the major folding event. Circular dichroism and fluorescence indicate that negligible amounts of secondary and tertiary structures form in the burst phase. Thus, for these two denatured states, collapse and secondary structure formation are not energetically downhill processes even under aqueous, low-denaturant conditions. In addition, water appears to be as good a solvent as that with high concentrations of denaturant, when considering the over-all dimensions of the denatured state. However, the removal of denaturant does subtly alter the distribution of backbone dihedral phi,psi angles, most likely resulting in a shift from the polyproline II region to the helical region of the Ramachandran map. We consider the thermodynamic origins of these behaviors along with implications for folding mechanisms and computer simulations thereof.     

Reassessing the folding of the KIX domain: Evidence for a two‐state mechanism

The debate about the presence and role of intermediates in the folding of proteins has been a critical issue, especially for fast folders. One of the classical methodologies to identify such metastable species is the “burst-phase analysis,” whereby the observed signal amplitude from stopped-flow traces is determined as a function of denaturant concentration. However, a complication may arise when folding is sufficiently fast to jeopardize the reliability of the stopped-flow technique. In this study, we reassessed the folding of the KIX domain from cAMP Response Element-Binding (CREB)-binding protein, which has been proposed to involve the formation of an intermediate that accumulates in the dead time of the stopped flow. By using an in-house-built capillary continuous flow with a 50-μs dead time, we demonstrate that this intermediate is not present; the problem arose because of the instrumental limitation of the standard stopped flow to assess very fast refolding rate constants (e.g., ≥500 s⁻¹).     

Fast folding of a prototypic polypeptide: the immunoglobulin binding domain of streptococcal protein G.

The folding of the small (56 residues) highly stable B1 immunoglobulin binding domain (GB1) of streptococcal protein G has been investigated by quenched-flow deuterium-hydrogen exchange. This system represents a paradigm for the study of protein folding because it exhibits no complicating features superimposed upon the intrinsic properties of the polypeptide chain. Collapse to a semicompact state exhibiting partial order, reflected in protection factors for ND-NH exchange up to 10-fold higher than that expected for a random coil, occurs within the dead time (< or = 1 ms) of the quenched flow apparatus. This is followed by the formation of the fully native state, as monitored by the fractional proton occupancy of 26 backbone amide groups spread throughout the protein, in a single rapid concerted step with a half-life of 5.2 ms at 5 degrees C.