The cooperativity of burst phase reactions explored.

The denaturant-dependence of the major, observable relaxation rates for folding (kobs) of ribonuclease HI from Escherichia coli (RNase H) and phage T4 lysozyme (T4L) reveal that, for both proteins, folding begins with the rapid and transient accumulation of intermediate species in a "burst phase" which precedes the rate-limiting formation of the native state; this is evidenced by a "rollover" in the folding limb of the rate profiles (kobs versus denaturant, or chevron plot). These rate profiles are most simply described by a three-state mechanism (unfolded-to-intermediate-to-native), which implies that the burst phase represents a transition between two distinct thermodynamic states. It is shown here that the equilibrium properties of these burst phase reactions can be equally well modeled by a mechanism involving a continuum of states where the free energy of each state is linearly related to its m-value (the parameter describing the linear relationship between free energy and denaturant). A numerical model is also developed to describe the time evolution of such a system, which exhibits nearly perfect exponential behavior. Both models emphasize how a continuum of states operating under a linear free energy relationship may behave like a two state system. Such a scheme finds experimental justification from an interpretation of recent native state hydrogen exchange data. The analytical model described for a continuum can account for the observed kinetic profiles of several other model proteins. The results, however, appear context specific, suggesting that burst phase reactions are not entirely random and non-specific. The results reported in this study have important implications for the concept of cooperativity in protein folding reactions.     

Hydration and packing along the folding pathway of SH3 domains by pressure-dependent NMR

The volumetric properties associated with protein folding transitions reflect changes in protein packing and hydration of the states that participate in the folding reaction. Here, NMR spin relaxation techniques are employed to probe the folding-unfolding kinetics of two SH3 domains as a function of pressure so that the changes in partial molar volumes along the folding pathway can be measured. The two domains fold with rates that differ by approximately 3 orders of magnitude, so their folding dynamics must be probed using different NMR relaxation experiments. In the case of the drkN SH3 domain that folds via a two-state mechanism on a time scale of seconds, nitrogen magnetization exchange spectroscopy is employed, while for the G48M mutant of the Fyn SH3 domain where the folding occurs on the millisecond time scale (three-step reaction), relaxation dispersion experiments are utilized. The NMR methodology is extremely sensitive to even small changes in equilibrium and rate constants, so reliable estimates of partial molar volumes can be obtained using low pressures (1-120 bar), thus minimizing perturbations to any of the states along the folding reaction coordinate. The volumetric data that were obtained are consistent with a similar folding mechanism for both SH3 domains, involving early chain compaction to states that are at least partially hydrated. This work emphasizes the role of NMR spin relaxation in studying dynamic processes over a wide range of time scales.     

Picosecond dynamical changes on denaturation of yeast phosphoglycerate kinase revealed by quasielastic neutron scattering

Quasielastic neutron scattering experiments performed on yeast phosphoglycerate kinase in the native form and denatured in 1.5 M guanidinium chloride reveal a change in the fast (picosecond time scale) diffusive internal dynamics of the protein. The momentum and energy transfer dependences of the scattering for both states are fitted by an analytical model in which, on the experimentally accessible picosecond time scale and angstrom length scale, the dynamics of a fraction of the nonexchangeable hydrogens in the protein is described as a superposition of vibrations with uniform diffusion in a sphere, the rest of the hydrogens undergoing only vibrational motion. The fraction diffusing changes, from ≈60% in the native protein to ≈82% in the denatured protein. The radius of the sphere also changes slightly, from ≈1.8 Å in the native protein to ≈2.2 Å in the denatured protein. Possible implications of these results for the general protein folding problem are discussed. Proteins 28:380–387, 1997 © 1997 Wiley-Liss, Inc.     

Cunning simplicity of protein folding landscapes

Funnel-like landscapes are widely used to visualize protein folding. It might seem that any funnel-like energy landscape helps to avoid the 'Levinthal paradox', i.e. to avoid sampling the impossibly large number of conformations for a folding protein. This cunning suggestion, reinforced by beautiful drawings of the energy funnels, stimulated some simple models of protein folding; one of them [D.J. Bicout and A. Szabo (2000) Protein Sci., 9, 452-465] is especially straightforward and instructive. A thorough analysis of this strict funnel model (which does not consider a nucleation of phase separation in the course of folding) shows that it cannot provide a simultaneous explanation for both major features observed for protein folding: (i) folding within non-astronomical time, and (ii) co-existence of the native and the unfolded states during the folding process. On the contrary, the nucleation mechanism of protein folding can account for both these major features simultaneously.     

The Energy Landscape,Folding Pathways and the Kinetics of a Knotted Protein

The folding pathway and rate coefficients of the folding of a knotted protein are calculated for a potential energy function with minimal energetic frustration. A kinetic transition network is constructed using the discrete path sampling approach, and the resulting potential energy surface is visualized by constructing disconnectivity graphs. Owing to topological constraints, the low-lying portion of the landscape consists of three distinct regions, corresponding to the native knotted state and to configurations where either the N or C terminus is not yet folded into the knot. The fastest folding pathways from denatured states exhibit early formation of the N terminus portion of the knot and a rate-determining step where the C terminus is incorporated. The low-lying minima with the N terminus knotted and the C terminus free therefore constitute an off-pathway intermediate for this model. The insertion of both the N and C termini into the knot occurs late in the folding process, creating large energy barriers that are the rate limiting steps in the folding process. When compared to other protein folding proteins of a similar length, this system folds over six orders of magnitude more slowly.     

The Role of High-Dimensional Diffusive Search,Stabilization, and Frustration in Protein Folding

Proteins are polymeric molecules with many degrees of conformational freedom whose internal energetic interactions are typically screened to small distances. Therefore, in the high-dimensional conformation space of a protein, the energy landscape is locally relatively flat, in contrast to low-dimensional representations, where, because of the induced entropic contribution to the full free energy, it appears funnel-like. Proteins explore the conformation space by searching these flat subspaces to find a narrow energetic alley that we call a hypergutter and then explore the next, lower-dimensional, subspace. Such a framework provides an effective representation of the energy landscape and folding kinetics that does justice to the essential characteristic of high-dimensionality of the search-space. It also illuminates the important role of nonnative interactions in defining folding pathways. This principle is here illustrated using a coarse-grained model of a family of three-helix bundle proteins whose conformations, once secondary structure has formed, can be defined by six rotational degrees of freedom. Two folding mechanisms are possible, one of which involves an intermediate. The stabilization of intermediate subspaces (or states in low-dimensional projection) in protein folding can either speed up or slow down the folding rate depending on the amount of native and nonnative contacts made in those subspaces. The folding rate increases due to reduced-dimension pathways arising from the mere presence of intermediate states, but decreases if the contacts in the intermediate are very stable and introduce sizeable topological or energetic frustration that needs to be overcome. Remarkably, the hypergutter framework, although depending on just a few physically meaningful parameters, can reproduce all the types of experimentally observed curvature in chevron plots for realizations of this fold.     

An expanded view of the protein folding landscape of PDZ domains

Most protein domains fold in an apparently co-operative and two-state manner with only the native and denatured states significantly populated at any experimental condition. However, the protein folding energy landscape is often rugged and different transition states may be rate limiting for the folding reaction under different conditions, as seen for the PDZ protein domain family. We have here analyzed the folding kinetics of two PDZ domains and found that a previously undetected third transition state is rate limiting under conditions that stabilize the native state relative to the denatured state. In light of these results, we have re-analyzed previous folding data on PDZ domains and present a unified folding mechanism with three distinct transition states separated by two high-energy intermediates. Our data show that sequence composition tunes the relative stabilities of folding transition states within the PDZ family, while the overall mechanism is determined by topology. This model captures the kinetic folding mechanism of all PDZ domains studied to date.     

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.     

Probing the folding and unfolding dynamics of secondary and tertiary structures in a three-helix bundle protein

Fast relaxation kinetics studies of the B-domain of staphylococcal protein A were performed to characterize the folding and unfolding of this small three-helix bundle protein. The relaxation kinetics were initiated using a laser-induced temperature jump and probed using time-resolved infrared spectroscopy. The kinetics monitored within the amide I' absorbance of the polypeptide backbone exhibit two distinct kinetics phases with nanosecond and microsecond relaxation times. The fast kinetics relaxation time is close to the diffusion limits placed on protein folding reactions. The fast kinetics phase is dominated by the relaxation of the solvated helix (nu = 1632 cm(-1)), which reports on the fast relaxation of the individual helices. The slow kinetics phase follows the cooperative relaxation of the native helical bundle core that is monitored by both solvated (nu = 1632 cm(-1)) and buried helical IR bands (nu = 1652 cm(-1)). The folding rates of the slow kinetics phase calculated over an extended temperature range indicate that the core formation of this protein follows a pathway that is energetically downhill. The unfolding rates are much more strongly temperature-dependent indicating an activated process with a large energy barrier. These results provide significant insight into the primary process of protein folding and suggest that fast formation of helices can drive the folding of helical proteins.     

Protein dynamics: From the native to the unfolded state and back again

Simulations to study protein unfolding and folding were performed. The unfolding simulations make use of molecular dynamics and treat an atomic model of barnase in aqueous solvent. The cooperative nature of the unfolding transition and the important role of water are described. The folding simulations are based on a bead model of the protein on a cubic lattice. It is shown for the 27-mer model that a large energy gap between the lowest energy (native) state and the excited states is a necessary and sufficient condition for fast folding.     

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.     

Analysis on folding of misgurin using two-dimensional HP model

Misgurin is an antimicrobial peptide from the loach, while the hydrophobic-polar (HP) model is a way to study the folding conformations and native states in peptide and protein although several amino acids cannot be classified either hydrophobic or polar. Practically, the HP model requires extremely intensive computations, thus it has yet to be used widely. In this study, we use the two-dimensional HP model to analyze all possible folding conformations and native states of misgurin with conversion of natural amino acids according to the normalized amino acid hydrophobicity index as well as the shortest benchmark HP sequence. The results show that the conversion of misgurin into HP sequence with glycine as hydrophobic amino acid at pH 2 has 1212 folding conformations with the same native state of minimal energy -6; the conversion of glycine as polar amino acid at pH 2 has 13,386 folding conformations with three native states of minimal energy -5; the conversion of glycine as hydrophobic amino acid at pH 7 has 2538 folding conformations with three native states of minimal energy -5; and the conversion of glycine as polar amino acid at pH 7 has 12,852 folding conformations with three native states of minimal energy -4. Those native states can be ranked according to the normalized amino acid hydrophobicity index. The detailed discussions suggest two ways to modify misgurin.     

Side-chain interactions in the folding pathway of a Fyn SH3 domain mutant studied by relaxation dispersion NMR spectroscopy

A major challenge to the study of protein folding is the fact that intermediate states along the reaction pathway are generally unstable and thus difficult to observe. Recently developed NMR relaxation dispersion experiments present an avenue to accessing such states, providing kinetic, thermodynamic, and structural information for intermediates with small (greater than or equal to approximately 1%) populations at equilibrium. We have employed these techniques to study the three-state folding reaction of the G48M Fyn SH3 domain. Using (13)C-, (1)H-, and (15)N-based methods, we have characterized backbone and side-chain interactions in the folded, unfolded, intermediate, and transition states, thereby mapping the energy landscape of the protein. We find that the intermediate, populated to approximately 1%, contains nativelike structure in a central beta-sheet, and is disordered at the amino and carboxy termini. The intermediate is stabilized by side-chain van der Waals contacts, yet (13)C chemical shifts indicate that methyl-containing residues remain disordered. This state has a partially structured backbone and a collapsed yet mobile hydrophobic core and thus closely resembles a molten globule. Nonpolar side-chain contacts are formed in the unfolded-intermediate transition state; these interactions are disrupted in the intermediate-folded transition state, possibly allowing side chains to rearrange as they adopt the native packing configuration. This work illustrates the power of novel relaxation dispersion experiments in characterizing excited states that are "invisible" in even the most sensitive of NMR experiments.     

Pressure-jump-induced kinetics reveals a hydration dependent folding/unfolding mechanism of ribonuclease A

Pressure-jump (p-jump)-induced relaxation kinetics was used to explore the energy landscape of protein folding/unfolding of Y115W, a fluorescent variant of ribonuclease A. Pressure-jumps of 40 MPa amplitude (5 ms dead-time) were conducted both to higher (unfolding) and to lower (folding) pressure, in the range from 100 to 500 MPa, between 30 and 50 degrees C. Significant deviations from the expected symmetrical protein relaxation kinetics were observed. Whereas downward p-jumps resulted always in single exponential kinetics, the kinetics induced by upward p-jumps were biphasic in the low pressure range and monophasic at higher pressures. The relative amplitude of the slow phase decreased as a function of both pressure and temperature. At 50 degrees C, only the fast phase remained. These results can be interpreted within the framework of a two-dimensional energy surface containing a pressure- and temperature-dependent barrier between two unfolded states differing in the isomeric state of the Asn-113-Pro-114 bond. Analysis of the activation volume of the fast kinetic phase revealed a temperature-dependent shift of the unfolding transition state to a larger volume. The observed compensation of this effect by glycerol offers an explanation for its protein stabilizing effect.     

Entropy and barrier-controlled fluctuations determine conformational viscoelasticity of single biomolecules

Biological macromolecules have complex and nontrivial energy landscapes, endowing them with a unique conformational adaptability and diversity in function. Hence, understanding the processes of elasticity and dissipation at the nanoscale is important to molecular biology and emerging fields such as nanotechnology. Here we analyze single molecule fluctuations in an atomic force microscope, using a generic model of biopolymer viscoelasticity that includes local "internal" conformational dissipation. Comparing two biopolymers, dextran and cellulose (polysaccharides with and without local bistable transitions), demonstrates that signatures of simple conformational change are minima in both the elastic and internal friction constants around a characteristic force. A novel analysis of dynamics on a bistable energy landscape provides a simple explanation: an elasticity driven by the entropy, and friction by a barrier-controlled hopping time of populations between states, which is surprisingly distinct to the well-known relaxation time. This nonequilibrium microscopic analysis thus provides a means of quantifying new dynamical features of the energy landscape of the glucopyranose ring, revealing an unexpected underlying roughness and information on the shape of the barrier of the chair-boat transition in dextran. The results presented herein provide a basis toward probing the viscoelasticity of macromolecular conformational transitions on more complex energy landscapes, such as during protein folding.     

A heteropolymer model study for the mechanism of protein folding

A heteropolymer model of randomly self-interacting chains in two dimensions is studied with numerical simulations in order to elucidate the folding mechanism of protein. We find that the model occasionally shows folding propensity depending on the sequence of random numbers given to the chain. We study the thermodynamic and kinematic roles in the folding mechanism by grouping the local energy minima found in the simulations into clusters according to the similarity of their conformations. It is suggested that the local minima to which some heteropolymers show a folding tendency are always the lowest energy states of the energy spectrum within a cluster, though which cluster is selected depends on the sequence. For the eight random sequences we study, we find that the energy gap between the ground state and excited states is little correlated with folding or nonfolding. We rather find that folding propensities are correlated with the global structure of the average energy surface, implying a dominant kinetic role in the folding mechanism, although thermal factors cannot be ignored as the mechanism of choosing the ground state within a cluster of states connected by small deformations. We suggest that a hierarchical cluster structure plays an important role in selecting a unique folded state out of the huge number of local minima of heteropolymers. © 1997 John Wiley & Sons, Inc.     

Toward minimalist models of larger proteins: a ubiquitin-like protein

Our recently developed off-lattice bead model capable of simulating protein structures with mixed alpha/beta content has been extended to model the folding of a ubiquitin-like protein and provides a means for examining the more complex kinetics involved in the folding of larger proteins. Using trajectories generated from constant-temperature Langevin dynamics simulations and sampling with the multiple multi-histogram method over five-order parameters, we are able to characterize the free energy landscape for folding and find evidence for folding through compact intermediates. Our model reproduces the observation that the C-terminus loop structure in ubiquitin is the last to fold in the folding process and most likely plays a spectator role in the folding kinetics. The possibility of a productive metastable intermediate along the folding pathway consisting of collapsed states with no secondary structure, and of intermediates or transition structures involving secondary structural elements occurring early in the sequence, is also supported by our model. The kinetics of folding remain multi-exponential below the folding temperature, with glass-like kinetics appearing at T/T(f) approximately 0.86. This new physicochemical model, designed to be predictive, helps validate the value of modeling protein folding at this level of detail for genomic-scale studies, and motivates further studies of other protein topologies and the impact of more complex energy functions, such as the addition of solvation forces.     

Untangling the complexity of membrane protein folding

Delineating the folding steps of helical-bundle membrane proteins has been a challenging task. Many questions remain unanswered, including the conformation and stability of the states populated during folding, the shape of the energy barriers between the states, and the role of lipids as a solvent in mediating the folding. Recently, theoretical frames have matured to a point that permits detailed dissection of the folding steps, and advances in experimental techniques at both single-molecule and ensemble levels enable selective modulation of specific steps for quantitative determination of the folding energy landscapes. We also discuss how lipid molecules would play an active role in shaping the folding energy landscape of membrane proteins, and how folding of multi-domain membrane proteins can be understood based on our current knowledge. We conclude this review by offering an outlook for emerging questions in the study of membrane protein folding.     

Interpretation of protein folding psi values

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