Topological determinants of protein folding

The folding of many small proteins is kinetically a two-state process that represents overcoming the major free-energy barrier. A kinetic characteristic of a conformation, its probability to descend to the native state domain in the amount of time that represents a small fraction of total folding time, has been introduced to determine to which side of the free-energy barrier a conformation belongs. However, which features make a protein conformation on the folding pathway become committed to rapidly descending to the native state has been a mystery. Using two small, well characterized proteins, CI2 and C-Src SH3, we show how topological properties of protein conformations determine their kinetic ability to fold. We use a macroscopic measure of the protein contact network topology, the average graph connectivity, by constructing graphs that are based on the geometry of protein conformations. We find that the average connectivity is higher for conformations with a high folding probability than for those with a high probability to unfold. Other macroscopic measures of protein structural and energetic properties such as radius of gyration, rms distance, solvent-accessible surface area, contact order, and potential energy fail to serve as predictors of the probability of a given conformation to fold.     

Excluded volume, local structural cooperativity, and the polymer physics of protein folding rates

A coarse-grained variational model is used to investigate the polymer dynamics of barrier crossing for a diverse set of two-state folding proteins. The model gives reliable folding rate predictions provided excluded volume terms that induce minor structural cooperativity are included in the interaction potential. In general, the cooperative folding routes have sharper interfaces between folded and unfolded regions of the folding nucleus and higher free energy barriers. The calculated free energy barriers are strongly correlated with native topology as characterized by contact order. Increasing the rigidity of the folding nucleus changes the local structure of the transition state ensemble nonuniformly across the set of proteins studied. Nevertheless, the calculated prefactors k(0) are found to be relatively uniform across the protein set, with variation in 1/k(0) less than a factor of 5. This direct calculation justifies the common assumption that the prefactor is roughly the same for all small two-state folding proteins. Using the barrier heights obtained from the model and the best-fit monomer relaxation time 30 ns, we find that 1/k(0) approximately 1-5 mus (with average 1/k(0) approximately 4 micros). This model can be extended to study subtle aspects of folding such as the variation of the folding rate with stability or solvent viscosity and the onset of downhill folding.     

Balancing energy and entropy: a minimalist model for the characterization of protein folding landscapes

Coarse-grained models have been extremely valuable in promoting our understanding of protein folding. However, the quantitative accuracy of existing simplified models is strongly hindered either from the complete removal of frustration (as in the widely used Gō-like models) or from the compromise with the minimal frustration principle and/or realistic protein geometry (as in the simple on-lattice models). We present a coarse-grained model that "naturally" incorporates sequence details and energetic frustration into an overall minimally frustrated folding landscape. The model is coupled with an optimization procedure to design the parameters of the protein Hamiltonian to fold into a desired native structure. The application to the study of src-Src homology 3 domain shows that this coarse-grained model contains the main physical-chemical ingredients that are responsible for shaping the folding landscape of this protein. The results illustrate the importance of nonnative interactions and energetic heterogeneity for a quantitative characterization of folding mechanisms.     

The folding mechanism of larger model proteins: role of native structure.

The folding mechanism of a 125-bead heteropolymer model for proteins is investigated with Monte Carlo simulations on a cubic lattice. Sequences that do and do not fold in a reasonable time are compared. The overall folding behavior is found to be more complex than that of models for smaller proteins. Folding begins with a rapid collapse followed by a slow search through the semi-compact globule for a sequence-dependent stable core with about 30 out of 176 native contacts which serves as the transition state for folding to a near-native structure. Efficient search for the core is dependent on structural features of the native state. Sequences that fold have large amounts of stable, cooperative structure that is accessible through short-range initiation sites, such as those in anti-parallel sheets connected by turns. Before folding is completed, the system can encounter a second bottleneck, involving the condensation and rearrangement of surface residues. Overly stable local structure of the surface residues slows this stage of the folding process. The relation of the results from the 125-mer model studies to the folding of real proteins is discussed.     

Prediction of protein-folding mechanisms from free-energy landscapes derived from native structures

Guided by recent experimental results suggesting that protein-folding rates and mechanisms are determined largely by native-state topology, we develop a simple model for protein folding free-energy landscapes based on native-state structures. The configurations considered by the model contain one or two contiguous stretches of residues ordered as in the native structure with all other residues completely disordered; the free energy of each configuration is the difference between the entropic cost of ordering the residues, which depends on the total number of residues ordered and the length of the loop between the two ordered segments, and the favorable attractive interactions, which are taken to be proportional to the total surface area buried by the ordered residues in the native structure. Folding kinetics are modeled by allowing only one residue to become ordered/disordered at a time, and a rigorous and exact method is used to identify free-energy maxima on the lowest free-energy paths connecting the fully disordered and fully ordered configurations. The distribution of structure in these free-energy maxima, which comprise the transition-state ensemble in the model, are reasonably consistent with experimental data on the folding transition state for five of seven proteins studied. Thus, the model appears to capture, at least in part, the basic physics underlying protein folding and the aspects of native-state topology that determine protein-folding mechanisms.     

Symmetry and frustration in protein energy landscapes: a near degeneracy resolves the Rop dimer-folding mystery

Protein folding has become one of the best understood biochemical reactions from a kinetic viewpoint. The funneled energy landscape, a consequence of the minimal frustration achieved by evolution in sequences, explains how most proteins fold efficiently and robustly to their functional structure and allows robust prediction of folding kinetics. The folding of Rop (repressor of primer) dimer is exceptional because some of its mutants with a redesigned hydrophobic core both fold and unfold much faster than the WT protein, which seems to conflict with a simple funneled energy landscape for which topology mainly determines the kinetics. We propose that the mystery of Rop folding can be unraveled by assuming a double-funneled energy landscape on which there are two basins that correspond to distinct but related topological structures. Because of the near symmetry of the molecule, mutations can cause a conformational switch to a nearly degenerate yet distinct topology or lead to a mixture of both topologies. The topology predicted to have the lower free-energy barrier height for folding was further found by all-atom modeling to give a better structural fit for those mutants with the extreme folding and unfolding rates. Thus, the non-Hammond effects can be understood within energy-landscape theory if there are in fact two different but nearly degenerate structures for Rop. Mutations in symmetric and regular structures may give rise to frustration and thus result in degeneracy.     

Structure-function-folding relationship in a WW domain

Protein folding barriers result from a combination of factors including unavoidable energetic frustration from nonnative interactions, natural variation and selection of the amino acid sequence for function, and/or selection pressure against aggregation. The rate-limiting step for human Pin1 WW domain folding is the formation of the loop 1 substructure. The native conformation of this six-residue loop positions side chains that are important for mediating protein-protein interactions through the binding of Pro-rich sequences. Replacement of the wild-type loop 1 primary structure by shorter sequences with a high propensity to fold into a type-I' beta-turn conformation or the statistically preferred type-I G1 bulge conformation accelerates WW domain folding by almost an order of magnitude and increases thermodynamic stability. However, loop engineering to optimize folding energetics has a significant downside: it effectively eliminates WW domain function according to ligand-binding studies. The energetic contribution of loop 1 to ligand binding appears to have evolved at the expense of fast folding and additional protein stability. Thus, the two-state barrier exhibited by the wild-type human Pin1 WW domain principally results from functional requirements, rather than from physical constraints inherent to even the most efficient loop formation process.     

Kinetic traps in lysozyme folding.

Folding of lysozyme from hen egg white was investigated by using interrupted refolding experiments. This method makes use of a high energy barrier between the native state and transient folding intermediates, and, in contrast to conventional optical techniques, it enables one to specifically monitor the amount of native molecules during protein folding. The results show that under strongly native conditions lysozyme can refold on parallel pathways. The major part of the lysozyme molecules (86%) refold on a slow kinetic pathway with well-populated partially folded states. Additionally, 14% of the molecules fold faster. The rate constant of formation of native molecules on the fast pathway corresponds well to the rate constant expected for folding to occur by a two-state process without any detectable intermediates. The results suggest that formation of the native state for the major fraction of lysozyme molecules is retarded compared with the direct folding process. Partially structured intermediates that transiently populate seem to be kinetically trapped in a conformation that can only slowly reach the native structure.     

Conformation, energy, and folding ability of selected amino acid sequences.

Evolutionary selection of sequences is studied with a knowledge-based Hamiltonian to find the design principle for folding to a model protein structure. With sequences selected by naive energy minimization, the model structure tends to be unstable and the folding ability is low. Sequences with high folding ability have only the low-lying energy minimum but also an energy landscape which is similar to that found for the native sequence over a wide region of the conformation space. Though there is a large fluctuation in foldable sequences, the hydrophobicity pattern and the glycine locations are preserved among them. Implications of the design principle for the molecular mechanism of folding are discussed.     

Comparison of successive transition states for folding reveals alternative early folding pathways of two homologous proteins

The energy landscape theory provides a general framework for describing protein folding reactions. Because a large number of studies, however, have focused on two-state proteins with single well-defined folding pathways and without detectable intermediates, the extent to which free energy landscapes are shaped up by the native topology at the early stages of the folding process has not been fully characterized experimentally. To this end, we have investigated the folding mechanisms of two homologous three-state proteins, PTP-BL PDZ2 and PSD-95 PDZ3, and compared the early and late transition states on their folding pathways. Through a combination of Φ value analysis and molecular dynamics simulations we obtained atomic-level structures of the transition states of these homologous three-state proteins and found that the late transition states are much more structurally similar than the early ones. Our findings thus reveal that, while the native state topology defines essentially in a unique way the late stages of folding, it leaves significant freedom to the early events, a result that reflects the funneling of the free energy landscape toward the native state.     

Quantitative criteria for native energetic heterogeneity influences in the prediction of protein folding kinetics

Energy landscape theory requires that the protein-folding mechanism is generally globally directed or funneled toward the native state. The collective nature of transition state ensembles further suggests that sufficient averaging of the native interactions can occur so that the knowledge of the native topology may suffice for predicting the mechanism. Nevertheless, while simple homogeneously weighted native topology-based models predict the folding mechanisms for many proteins, for other proteins knowledge of the native topology, by itself, seems not to suffice in determining the folding mechanism. Simulations of proteins with differing topologies reveal that the failure of homogeneously weighted topology-based models can, however, be completely understood within the framework of a funneled energy landscape and can be quantified by comparing the fluctuation of entropy cost for forming contacts to the expected fluctuations in contact energy. To be precise, we find the transition state ensembles of proteins with all-α topologies, which are more uniform in the specific entropy cost of contact formation, have transition state ensembles that are more readily perturbed by differences in energetic weights than are the transition state ensembles of proteins with significant amounts of β-structure, where the specific entropy costs of contact formation are more widely distributed. This behavior is consistent with a random-field Ising model analogy that follows from the free energy functional approach to folding.     

Structures and relative free energies of partially folded states of proteins

The ability of proteins to fold to well defined compact structures is one of the most remarkable examples of the effect of natural selection on biological molecules. To understand their properties, including the stability, the mechanism of folding, and the possibilities of misfolding and association, it is necessary to know the protein free energy landscape. We use NMR data as restraints in a Monte Carlo sampling procedure to determine the ensemble of structures populated by human alpha-lactalbumin in the presence of increasing concentrations of urea. The ensembles of structures that represent the partially folded states of the protein show that two structural cores, corresponding to portions of the alpha and beta domains of the native protein, are preserved even when the native-like interactions that define their existence are substantially weakened. Analysis of the network of residual contacts reveals the presence of a complex interface region between the two structural cores and indicates that the development of specific interactions within this interface is the key step in achieving the native structure. The relative probabilities of the conformations determined from the NMR data are used to construct a coarse-grained free energy landscape for alpha-lactalbumin in the absence of urea. The form of the landscape, together with the existence of distinct cores, supports the concept that robustness and modularity are the properties that make possible the folding of complex proteins.     

Competition between native topology and nonnative interactions in simple and complex folding kinetics of natural and designed proteins

We compared folding properties of designed protein Top7 and natural protein S6 by using coarse-grained chain models with a mainly native-centric construct that accounted also for nonnative hydrophobic interactions and desolvation barriers. Top7 and S6 have similar secondary structure elements and are approximately equal in length and hydrophobic composition. Yet their experimental folding kinetics were drastically different. Consistent with experiment, our simulated folding chevron arm for Top7 exhibited a severe rollover, whereas that for S6 was essentially linear, and Top7 model kinetic relaxation was multiphasic under strongly folding conditions. The peculiar behavior of Top7 was associated with several classes of kinetic traps in our model. Significantly, the amino acid residues participating in nonnative interactions in trapped conformations in our Top7 model overlapped with those deduced experimentally. These affirmations suggest that the simple ingredients of native topology plus sequence-dependent nonnative interactions are sufficient to account for some key features of protein folding kinetics. Notably, when nonnative interactions were absent in the model, Top7 chevron rollover was not correctly predicted. In contrast, nonnative interactions had little effect on the quasi linearity of the model folding chevron arm for S6. This intriguing distinction indicates that folding cooperativity is governed by a subtle interplay between the sequence-dependent driving forces for native topology and the locations of favorable nonnative interactions entailed by the same sequence. Constructed with a capability to mimic this interplay, our simple modeling approach should be useful in general for assessing a designed sequence’s potential to fold cooperatively.     

Dynamics of one-state downhill protein folding

The small helical protein BBL has been shown to fold and unfold in the absence of a free energy barrier according to a battery of quantitative criteria in equilibrium experiments, including probe-dependent equilibrium unfolding, complex coupling between denaturing agents, characteristic DSC thermogram, gradual melting of secondary structure, and heterogeneous atom-by-atom unfolding behaviors spanning the entire unfolding process. Here, we present the results of nanosecond T-jump experiments probing backbone structure by IR and end-to-end distance by FRET. The folding dynamics observed with these two probes are both exponential with common relaxation times but have large differences in amplitude following their probe-dependent equilibrium unfolding. The quantitative analysis of amplitude and relaxation time data for both probes shows that BBL folding dynamics are fully consistent with the one-state folding scenario and incompatible with alternative models involving one or several barrier crossing events. At 333 K, the relaxation time for BBL is 1.3 μs, in agreement with previous folding speed limit estimates. However, late folding events at room temperature are an order of magnitude slower (20 μs), indicating a relatively rough underlying energy landscape. Our results in BBL expose the dynamic features of one-state folding and chart the intrinsic time-scales for conformational motions along the folding process. Interestingly, the simple self-averaging folding dynamics of BBL are the exact dynamic properties required in molecular rheostats, thus supporting a biological role for one-state folding.     

Folding funnels and frustration in off-lattice minimalist protein landscapes

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a β-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.     

Multiple routes lead to the native state in the energy landscape of the beta-trefoil family

In general, the energy landscapes of real proteins are sufficiently well designed that the depths of local energetic minima are small compared with the global bias of the native state. Because of the funneled nature of energy landscapes, models that lack energetic frustration have been able to capture the main structural features of the transition states and intermediates found in experimental studies of both small and large proteins. In this study we ask: Are the experimental differences in folding mechanisms among members of a particular structural family due to local topological constraints that deviate from the tertiary fold common to the family? The beta-trefoil structural family members IL-1beta, hisactophilin, and acidic/basic FGFs were chosen to address this question. It has been observed that the topological landscape of the beta-trefoils allows for the population of diverse, geometrically disconnected routes that provide energetically similar but structurally distinct ways for this family to fold. Small changes in topology or energetics can alter the preferred route. Taken together, these results indicate that the global fold of the beta-trefoil family determines the energy landscape but that the routes accessed on that landscape might differ as a result of functional requirements of the individual family members.     

Full distance-resolved folding energy landscape of one single protein molecule

Kinetic bulk and single molecule folding experiments characterize barrier properties but the shape of folding landscapes between barrier top and native state is difficult to access. Here, we directly extract the full free energy landscape of a single molecule of the GCN4 leucine zipper using dual beam optical tweezers. To this end, we use deconvolution force spectroscopy to follow an individual molecule’s trajectory with high temporal and spatial resolution. We find a heterogeneous energy landscape of the GCN4 leucine zipper domain. The energy profile is divided into two stable C-terminal heptad repeats and two less stable repeats at the N-terminus. Energies and transition barrier positions were confirmed by single molecule kinetic analysis. We anticipate that deconvolution sampling is a powerful tool for the model-free investigation of protein energy landscapes.