Solvent-tuning the collapse and helix formation time scales of lambda(6-85)*

The lambda(6-85)(*) pseudo-wild type of lambda repressor fragment is a fast two-state folder (k(f) approximately 35 microsec(-1) at 58 degrees C). Previously, highly stable lambda(6-85)(*) mutants with k(f) > 30 microsec(-1) have been engineered to fold nearly or fully downhill. Stabilization of the native state by solvent tuning might also tune lambda(6-85)(*) away from two-state folding. We test this prediction by examining the folding thermodynamics and kinetics of lambda(6-85)(*) in a stabilizing solvent, 45% by weight aqueous ethylene glycol at -28 degrees C. Detection of kinetics by circular dichroism at 222 nm (sensitive to helix content) and small angle X-ray scattering (measuring the radius of gyration) shows that refolding from guanidine hydrochloride denatured conditions exhibits very different time scales for collapse and secondary structure formation: the two processes become decoupled. Collapse remains a low-barrier activated process, while the fastest of several secondary structure formation time scales approaches the downhill folding limit. Two-state folding of lambda(6-85)(*) is not a robust process.     

Tuning lambda6-85 towards downhill folding at its melting temperature

The five-helix bundle lambda6-85* is a fast two-state folder. Several stabilized mutants have been reported to fold kinetically near-downhill or downhill. These mutants undergo a transition to two-state folding kinetics when heated. It has been suggested that this transition is caused by increased hydrophobicity at higher temperature. Here we investigate two histidine-containing mutants of lambda6-85* to see if a weaker hydrophobic core can extend the temperature range of downhill folding. The very stable lambdaHA is the fastest-folding lambda repressor to date (k(f)(-1) approximately k(obs)(-1)=2.3 micros at 44 degrees C). It folds downhill at low temperature, but transits back to two-state folding at its unfolding midpoint. lambdaHG has a weakened hydrophobic core. It is less stable than some slower folding mutants of lambda6-85*, and it has more exposed hydrophobic surface area in the folded state. This mutant nonetheless folds very rapidly, and has the non-exponential folding kinetics of an incipient downhill folder even at the unfolding midpoint (k(m)(-1) approximately 2 micros, k(a)(-1)=15 micros at 56 degrees C). We also compare the thermodynamic melting transition of lambdaHG with the nominal two-state folding mutant lambdaQG, which has a similar melting temperature. Unlike lambdaQG, lambdaHG yields fluorescence wavelength-dependent cooperativities and probe-dependent melting temperatures. This result combined with previous work shows that the energy landscapes of lambda repressor mutants support all standard folding mechanisms.     

Rate-temperature relationships in lambda-repressor fragment lambda 6-85 folding

Two classes of lambda(6-85) mutants (those richer in alanine, and those richer in glycine) have very similar slopes in an Arrhenius plot of the unfolding rates but very different temperature dependencies of the folding rates. Temperature-dependent interactions (e.g., hydrophobicity) play a large role in the initial stages of folding but not in the initial stages of unfolding of lambda(6-85). Placement of the transition state in terms of its surface exposure and entropy shows that at least two reaction coordinates are required to describe folding of all mutants over the full temperature range. The unusual Arrhenius plots of the very fastest mutant provide an additional kinetic signature for downhill folding.     

The transition state for folding of a peripheral subunit-binding domain contains robust and ionic-strength dependent characteristics

The denaturant dependencies of the folding and unfolding kinetics were used to characterize the structure of the transition state for folding of E3BD, a peripheral subunit-binding domain. For the majority of E3BD mutants, the Phi-values calculated at 298 K from the analysis of chevron plots were in good agreement with those previously determined at 325 K using Arrhenius analysis. This agreement further demonstrates the general robustness of Phi-value analyses, since different experiments, methods of denaturation and thermodynamic assumptions were used to determine each set of Phi(F) values. The structure of the transition state for folding was grossly conserved at 298 K and 325 K, with residues in Helix I playing a lesser role in folding than those located in the 3(10) helix, disordered loop and Helix II. However, the energetic contributions of a cluster of basic residues close to the N-terminus and Helix I, which are an integral part of the ligand-binding site, were susceptible to ionic strength effects because of electrostatic strain in native and transition states of E3BD at low ionic strength. We found no evidence of the downhill folding previously proposed for E3BD, even though the conditions employed in this study significantly increased the energetic bias towards the native state.     

Analyzing the effect of homogeneous frustration in protein folding

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a C_α structure‐based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well‐separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non‐native contact interactions in different folding scenarios. These findings strongly correlate with the protein free‐energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.Proteins 2013; 81:1727–1737. © 2013 Wiley Periodicals, Inc.     

Folding kinetics of phage 434 Cro protein

Folding kinetics for phage 434 Cro protein are examined and compared with those reported for lambda(6-85), the N-terminal domain of the repressor of phage lambda. The two proteins have similar all-helical structures consisting of five helices but different stabilities. In contrast to lambda(6-85), sharp and distinct aromatic (1)H NMR signals without exchange broadening characterize the native and urea-denatured 434 Cro forms at equilibrium at 20 degrees C, indicating slow interconversion on the NMR time scale. Stopped-flow fluorescence data using the single 434 Cro tryptophan indicate strongly urea-dependent refolding rates and smaller urea dependencies of the unfolding rates, suggesting a native-like transition state ensemble. Refolding rates are slower and unfolding rates considerably faster at pH 4 than at pH 6. This accounts for the lower stability of 434 Cro at pH 4 and suggests the existence of pH-dependent, possibly salt bridge interactions that are more stabilizing at pH 6. At <2 M urea, decreased folding amplitudes and nonlinear urea dependencies that are apparent at pH 6 indicate deviation from two-state behavior and suggest the formation of an early folding intermediate. The folding behavior of 434 Cro and why it folds 2 orders of magnitude slower than lambda(6-85) are rationalized in terms of the lower intrinsic helix stabilities and putative charge interactions in 434 Cro.     

Folding pathway dependence on energetic frustration and interaction heterogeneity for a three-dimensional hydrophobic protein model

Monte Carlo simulations of a hydrophobic protein model of 40 monomers in the cubic lattice are used to explore the effect of energetic frustration and interaction heterogeneity on its folding pathway. The folding pathway is described by the dependence of relevant conformational averages on an appropriate reaction coordinate, pfold, defined as the probability for a given conformation to reach the native structure before unfolding. We compare the energetically frustrated and heterogeneous hydrophobic potential, according to which individual monomers have a higher or lower tendency to form contacts unspecifically depending on their hydrophobicities, to an unfrustrated homogeneous Go-type potential with uniformly attractive native interactions and neutral non-native interactions (called Go1 in this study), and to an unfrustrated heterogeneous potential with neutral non-native interactions and native interactions having the same energy as the hydrophobic potential (called Go2 in this study). Folding kinetics are slowed down dramatically when energetic frustration increases, as expected and previously observed in a two-dimensional model. Contrary to our previous results in two dimensions, however, it appears that the folding pathway and transition state ensemble can be significantly dependent on the energy function used to stabilize the native structure. The sequence of events along the reaction coordinate, or the order along this coordinate in which different regions of the native conformation become structured, turns out to be similar for the hydrophobic and Go2 potentials, but with analogous events tending to occur at lower pfold values in the first case. In particular, the transition state obtained from the ensemble around pfold = 0.5 is more structured for the hydrophobic potential. For Go1, not only the transition state ensemble but the order of events itself is modified, suggesting that interaction heterogeneity, in addition to energetic frustration, can have significant effects on the folding mechanism, most likely by modifying the probability of different contacts in the unfolded state, the starting point for the folding reaction. Although based on a simple model, these results provide interesting insight into how sequence-dependent switching between folding pathways might occur in real proteins.     

Squeezed exponential kinetics to describe a nonglassy downhill folding as observed in a lattice protein model

We previously studied the so-called strange kinetics in the two-dimensional lattice HP model. To further study the strange kinetics, folding processes of a 27-mer cubic lattice protein model with Gō potential were investigated by simulating how the bundle of folding trajectories, consisting of a number of independent Monte Carlo simulations, evolves as the folding reaction proceeds, covering a wide range of temperature. Three realms of folding kinetics were observed depending on temperature. Although at temperatures where folding was two-state-like, the kinetics was conventional single exponential, we found that the time course data were well represented by a squeezed (or "shrunken") exponential function, exp [-(t/tau)beta] with beta > 1, at temperatures lower than the folding temperature, where folding was fastest and of a nonglassy downhill type. The squeezed exponential kinetics was found to pertain to the subdiffusion on the nonglassy downhill free energy surface and presents a marked contrast both to the single exponential kinetics and to the stretched exponential kinetics that was observed at lower temperatures where folding was also downhill but topological frustration came into effect. The observed temperature dependence of the folding kinetics suggests that some small single-domain proteins may follow the squeezed exponential kinetics at about the room temperature.     

Multiple routes and milestones in the folding of HIV-1 protease monomer

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV-1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV-1 protease.     

Topological and energetic factors: what determines the structural details of the transition state ensemble and "en-route" intermediates for protein folding? An investigation for small globular proteins

Recent experimental results suggest that the native fold, or topology, plays a primary role in determining the structure of the transition state ensemble, at least for small, fast-folding proteins. To investigate the extent of the topological control of the folding process, we studied the folding of simplified models of five small globular proteins constructed using a Go-like potential to retain the information about the native structures but drastically reduce the energetic frustration and energetic heterogeneity among residue-residue native interactions. By comparing the structure of the transition state ensemble (experimentally determined by Phi-values) and of the intermediates with those obtained using our models, we show that these energetically unfrustrated models can reproduce the global experimentally known features of the transition state ensembles and "en-route" intermediates, at least for the analyzed proteins. This result clearly indicates that, as long as the protein sequence is sufficiently minimally frustrated, topology plays a central role in determining the folding mechanism.     

Temperature-dependent downhill unfolding of ubiquitin. I. Nanosecond-to-millisecond resolved nonlinear infrared spectroscopy

Transient thermal unfolding of ubiquitin is investigated using nonlinear infrared spectroscopy after a nanosecond laser temperature jump (T-jump). The abrupt change in the unfolding free energy surface and the ns time resolution allow us to observe a fast response on ns to micros time-scales, which we attribute to downhill unfolding, before a cross-over to ms kinetics. The downhill unfolding by a sub-population of folded proteins is induced through a shift of the barrier toward the native state. By adjusting the T-jump width, the effect of the initial (T(i)) and final (T(f)) temperature on the unfolding dynamics can be separated. From the amplitude of the fast downhill unfolding, the fractional population prepared at the unfolding transition state is obtained. This population increases with both T(i) and with T(f). A two-state kinetic analysis of the ms refolding provides thermodynamic information about the barrier height. By a combination of the fast and slow unfolding and folding parameters, a quasi-two-state kinetic analysis is performed to calculate the time-dependent population changes of the folded state. This calculation coincides with the experimentally obtained population changes at low temperature but deviations are found in the T-jump from 67 to 78 degrees C. Using temperature-dependent barrier height changes, a temperature Phi value analysis is performed. The result shows a decreasing trend of Phi(T) with temperature, which indicates an increase of the heterogeneity of the transition state. We conclude that ubiquitin unfolds along a well-defined pathway at low temperature which expands with increasing temperature to include multiple routes.     

Probe-dependent and nonexponential relaxation kinetics: unreliable signatures of downhill protein folding

The theoretical suggestion that some proteins may encounter no significant free energy barriers during their folding raises an important question: What experimental signature does this downhill folding produce? Several authors have argued that nonexponential (and especially stretched exponential) or probe‐dependent kinetics represent useful experimental signatures of a downhill free energy surface. Here we examine more closely the connection between unusual kinetics and downhill energy surfaces. Simulation of diffusive relaxation dynamics on a variety of generically downhill, two‐dimensional free energy surfaces shows that these surfaces do not necessarily produce either probe‐dependent or significantly nonexponential kinetics. Conversely, we find that two‐dimensional surfaces with significant (>3k_B T) energy barriers can readily give rise to nonexponential and probe‐dependent kinetics. These results show that downhill folding does not constitute a necessary or sufficient condition for nonexponential and/or probe‐dependent folding kinetics. One cannot easily prove or disprove that a protein folds downhill simply by studying its relaxation kinetics. Proteins 2007. © 2007 Wiley‐Liss, Inc.     

A Survey of λ Repressor Fragments from Two-State to Downhill Folding

We survey the two-state to downhill folding transition by examining 20 λ_6-85^? mutants that cover a wide range of stabilities and folding rates. We investigated four new λ_6-85^? mutants designed to fold especially rapidly. Two were engineered using the core remodeling of Lim and Sauer, and two were engineered using Ferreiro et al.'s frustratometer. These proteins have probe-dependent melting temperatures as high as 80 °C and exhibit a fast molecular phase with the characteristic temperature dependence of the amplitude expected for downhill folding. The survey reveals a correlation between melting temperature and downhill folding previously observed for the β-sheet protein WW domain. A simple model explains this correlation and predicts the melting temperature at which downhill folding becomes possible. An X-ray crystal structure with a 1.64-Å resolution of a fast-folding mutant fragment shows regions of enhanced rigidity compared to the full wild-type protein.     

Folding of the protein domain hbSBD

The folding of the alpha-helix domain hbSBD of the mammalian mitochondrial branched-chain alpha-ketoacid dehydrogenase complex is studied by the circular dichroism technique in absence of urea. Thermal denaturation is used to evaluate various thermodynamic parameters defining the equilibrium unfolding, which is well described by the two-state model with the folding temperature T(F) = 317.8 +/- 1.95 K and the enthalpy change DeltaH(G) = 19.67 +/- 2.67 kcal/mol. The folding is also studied numerically using the off-lattice coarse-grained Go model and the Langevin dynamics. The obtained results, including the population of the native basin, the free-energy landscape as a function of the number of native contacts, and the folding kinetics, also suggest that the hbSBD domain is a two-state folder. These results are consistent with the biological function of hbSBD in branched-chain alpha-ketoacid dehydrogenase.     

Energetic and structural consequences of desolvation/solvation barriers to protein folding/unfolding assessed from experimental unfolding rates

Theoretical work has suggested the existence of solvation/desolvation barriers in protein folding/unfolding processes. We propose that the energetic and structural consequences of such barriers for the folding transition state can be assessed from experimental unfolding rates using well-established structure-energetics relationships. For a set of proteins of size within the 60-130 number-of-residues range, we find energetic effects associated to solvation/desolvation on the order of 10(2) kJ/mol. This supports that the folding transition states may be characterized by large networks of water-unsatisfied, broken internal contacts. In terms of buried surface, we estimate the typical network size to be on the order of several thousands of A2, or approximately 50% of the total change in accessible surface area upon unfolding. The analyses reported here thus suggest a clear structural picture for the different energetic balance of native and folding transition states.     

Transition states in protein folding kinetics: modeling phi-values of small beta-sheet proteins

Small single-domain proteins often exhibit only a single free-energy barrier, or transition state, between the denatured and the native state. The folding kinetics of these proteins is usually explored via mutational analysis. A central question is which structural information on the transition state can be derived from the mutational data. In this article, we model and structurally interpret mutational Φ-values for two small β-sheet proteins, the PIN and the FBP WW domains. The native structure of these WW domains comprises two β-hairpins that form a three-stranded β-sheet. In our model, we assume that the transition state consists of two conformations in which either one of the hairpins is formed. Such a transition state has been recently observed in molecular dynamics folding-unfolding simulations of a small designed three-stranded β-sheet protein. We obtain good agreement with the experimental data 1), by splitting up the mutation-induced free-energy changes into terms for the two hairpins and for the small hydrophobic core of the proteins; and 2), by fitting a single parameter, the relative degree to which hairpins 1 and 2 are formed in the transition state. The model helps us to understand how mutations affect the folding kinetics of WW domains, and captures also negative Φ-values that have been difficult to interpret.     

An all-atom structure-based potential for proteins: bridging minimal models with all-atom empirical forcefields

Protein dynamics take place on many time and length scales. Coarse-grained structure-based (Go) models utilize the funneled energy landscape theory of protein folding to provide an understanding of both long time and long length scale dynamics. All-atom empirical forcefields with explicit solvent can elucidate our understanding of short time dynamics with high energetic and structural resolution. Thus, structure-based models with atomic details included can be used to bridge our understanding between these two approaches. We report on the robustness of folding mechanisms in one such all-atom model. Results for the B domain of Protein A, the SH3 domain of C-Src Kinase, and Chymotrypsin Inhibitor 2 are reported. The interplay between side chain packing and backbone folding is explored. We also compare this model to a C(alpha) structure-based model and an all-atom empirical forcefield. Key findings include: (1) backbone collapse is accompanied by partial side chain packing in a cooperative transition and residual side chain packing occurs gradually with decreasing temperature, (2) folding mechanisms are robust to variations of the energetic parameters, (3) protein folding free-energy barriers can be manipulated through parametric modifications, (4) the global folding mechanisms in a C(alpha) model and the all-atom model agree, although differences can be attributed to energetic heterogeneity in the all-atom model, and (5) proline residues have significant effects on folding mechanisms, independent of isomerization effects. Because this structure-based model has atomic resolution, this work lays the foundation for future studies to probe the contributions of specific energetic factors on protein folding and function.     

Folding propensities of synthetic peptide fragments covering the entire sequence of phage 434 Cro protein.

The phage 434 Cro protein, the N-terminal domain of its repressor (R1-69) and that of phage lambda (lambda6-85) constitute a group of small, monomeric, single-domain folding units consisting of five helices with striking structural similarity. The intrinsic helix stabilities in lambda6-85 have been correlated to its rapid folding behavior, and a residual hydrophobic cluster found in R1-69 in 7 M urea has been proposed as a folding initiation site. To understand the early events in the folding of 434 Cro, and for comparison with R1-69 and lambda6-85, we examined the conformational behavior of five peptides covering the entire 434 Cro sequence in water, 40% (by volume) TFE/water, and 7 M urea solutions using CD and NMR. Each peptide corresponds to a helix and adjacent residues as identified in the native 434 Cro NMR and crystal structures. All are soluble and monomeric in the solution conditions examined except for the peptide corresponding to the 434 Cro helix 4, which has low water solubility. Helix formation is observed for the 434 Cro helix 1 and helix 2 peptides in water, for all the peptides in 40% TFE and for none in 7 M urea. NMR data indicate that the helix limits in the peptides are similar to those in the native protein helices. The number of side-chain NOEs in water and TFE correlates with the helix content, and essentially none are observed in 7 M urea for any peptide, except that for helix 5, where a hydrophobic cluster may be present. The low intrinsic folding propensities of the five helices could account for the observed stability and folding behavior of 434 Cro and is, at least qualitatively, in accord with the results of the recently described diffusion-collision model incorporating intrinsic helix propensities.     

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.