The chicken-egg scenario of protein folding revisited

What is the first step in protein folding - hydrophobic collapse (compaction) or secondary structure formation? It is still not clear if the major driving force in protein folding is hydrogen bonding or hydrophobic interactions or both. We analyzed data on the conformational characteristics of 41 globular proteins in native and partially folded conformational states. Our analysis shows that a good correlation exists between relative decrease in hydrodynamic volume and increase in secondary structure content. No compact equilibrium intermediates lacking secondary structure, or highly ordered non-compact species, were found. This correlation provides experimental support for the hypothesis that hydrophobic collapse occurs simultaneously with formation of secondary structure in the early stages of the protein folding.     

Specific non-native hydrophobic interactions in a hidden folding intermediate: implications for protein folding

Structures of intermediates and transition states in protein folding are usually characterized by amide hydrogen exchange and protein engineering methods and interpreted on the basis of the assumption that they have native-like conformations. We were able to stabilize and determine the high-resolution structure of a partially unfolded intermediate that exists after the rate-limiting step of a four-helix bundle protein, Rd-apocyt b(562), by multidimensional NMR methods. The intermediate has partial native-like secondary structure and backbone topology, consistent with our earlier native state hydrogen exchange results. However, non-native hydrophobic interactions exist throughout the structure. These and other results in the literature suggest that non-native hydrophobic interactions may occur generally in partially folded states. This can alter the interpretation of mutational protein engineering results in terms of native-like side chain interactions. In addition, since the intermediate exists after the rate-limiting step and Rd-apocyt b(562) folds very rapidly (k(f) approximately 10(4) s(-1)), these results suggest that non-native hydrophobic interactions, in the absence of topological misfolding, are repaired too rapidly to slow folding and cause the accumulation of folding intermediates. More generally, these results illustrate an approach for determining the high-resolution structure of folding intermediates.     

Comparing the folding and misfolding energy landscapes of phosphoglycerate kinase

Partitioning of polypeptides between protein folding and amyloid formation is of outstanding pathophysiological importance. Using yeast phosphoglycerate kinase as model, here we identify the features of the energy landscape that decide the fate of the protein: folding or amyloidogenesis. Structure formation was initiated from the acid-unfolded state, and monitored by fluorescence from 10 ms to 20 days. Solvent conditions were gradually shifted between folding and amyloidogenesis, and the properties of the energy landscape governing structure formation were reconstructed. A gradual transition of the energy landscape between folding and amyloid formation was observed. In the early steps of both folding and misfolding, the protein searches through a hierarchically structured energy landscape to form a molten globule in a few seconds. Depending on the conditions, this intermediate either folds to the native state in a few minutes, or forms amyloid fibers in several days. As conditions are changed from folding to misfolding, the barrier separating the molten globule and native states increases, although the barrier to the amyloid does not change. In the meantime, the native state also becomes more unstable and the amyloid more stable. We conclude that the lower region of the energy landscape determines the final protein structure.     

Protein stability and surface electrostatics: a charged relationship

Engineering proteins to withstand a broad range of conditions continues to be a coveted objective, holding the potential to advance biomedicine, industry, and our understanding of disease. One way of achieving this goal lies in elucidating the underlying interactions that define protein stability. It has been shown that the hydrophobic effect, hydrogen bonding, and packing interactions between residues in the protein interior are dominant factors that define protein stability. The role of surface residues in protein stability has received much less attention. It has been believed that surface residues are not important for protein stability particularly because their interactions with the solvent should be similar in the native and unfolded states. In the case of surface charged residues, it was sometimes argued that solvent exposure meant that the high dielectric of the solvent will further decrease the strength of the charge-charge interactions. In this paper, we challenge the notion that the surface charged residues are not important for protein stability. We computationally redesigned sequences of five different proteins to optimize the surface charge-charge interactions. All redesigned proteins exhibited a significant increase in stability relative to their parent proteins, as experimentally determined by circular dichroism spectroscopy and differential scanning calorimetry. These results suggest that surface charge-charge interactions are important for protein stability and that rational optimization of charge-charge interactions on the protein surface can be a viable strategy for enhancing protein stability.     

Im7 folding mechanism: misfolding on a path to the native state

Many proteins populate collapsed intermediate states during folding. In order to elucidate the nature and importance of these species, we have mapped the structure of the on-pathway intermediate of the four-helix protein, Im7, together with the conformational changes it undergoes as it folds to the native state. Kinetic data for 29 Im7 point mutants show that the intermediate contains three of the four helices found in the native structure, packed around a specific hydrophobic core. However, the intermediate contains many non-native interactions; as a result, hydrophobic interactions become disrupted in the rate-limiting transition state before the final helix docks onto the developing structure. The results of this study support a hierarchical mechanism of protein folding and explain why the misfolding of Im7 occurs. The data also demonstrate that non-native interactions can play a significant role in folding, even for small proteins with simple topologies.     

Amyloid fibrils formed by the programmed cell death regulator Bcl-xL

The accumulation of amyloid fibers due to protein misfolding is associated with numerous human diseases. For example, the formation of amyloid deposits in neurodegenerative pathologies is correlated with abnormal apoptosis. We report here the in vitro formation of various types of aggregates by Bcl-xL, a protein of the Bcl-2 family involved in the regulation of apoptosis. Bcl-xL forms aggregates in three states, micelles, native-like fibrils, and amyloid fibers, and their biophysical characterization has been performed in detail. Bcl-xL remains in its native state within micelles and native-like fibrils, and our results suggest that native-like fibrils are formed by the association of micelles. Formation of amyloid structures, that is, nonnative intermolecular β-sheets, is favored by the proximity of proteins within fibrils at the expense of the Bcl-xL native structure. Finally, we provide evidence of a direct relationship between the amyloid character of the fibers and the tertiary-structure stability of the native Bcl-xL. The potential causality between the accumulation of Bcl-xL into amyloid deposits and abnormal apoptosis during neurodegenerative diseases is discussed.     

Monte Carlo Simulations of Proteins in Cages: Influence of Confinement on the Stability of Intermediate States

We study the folding of small proteins inside confining potentials. Proteins are described using an effective potential model that contains the Ramachandran angles as degrees of freedom and does not need any a priori information about the native state. Hydrogen bonds, dipole-dipole-, and hydrophobic interactions are taken explicitly into account. An interesting feature displayed by this potential is the presence of metastable intermediates between the unfolded and native states. We consider different types of confining potentials to describe proteins folding inside cages with repulsive or attractive walls. Using the Wang-Landau algorithm, we determine the density of states and analyze in detail the thermodynamical properties of the confined proteins for different sizes of the cages. We show that confinement dramatically reduces the phase space available to the protein and that the presence of intermediate states can be controlled by varying the properties of the confining potential. Cages with strongly attractive walls destabilize the intermediate states and lead to a two-state folding into a configuration that is less stable than the native structure. However, cages with slightly attractive walls enhance the stability of native structure and induce a folding process, which occurs through intermediate configurations.     

Folding and stability of helical bundle proteins from coarse‐grained models

We develop a coarse‐grained model where solvent is considered implicitly, electrostatics are included as short‐range interactions, and side‐chains are coarse‐grained to a single bead. The model depends on three main parameters: hydrophobic, electrostatic, and side‐chain hydrogen bond strength. The parameters are determined by considering three level of approximations and characterizing the folding for three selected proteins (training set). Nine additional proteins (containing up to 126 residues) as well as mutated versions (test set) are folded with the given parameters. In all folding simulations, the initial state is a random coil configuration. Besides the native state, some proteins fold into an additional state differing in the topology (structure of the helical bundle). We discuss the stability of the native states, and compare the dynamics of our model to all atom molecular dynamics simulations as well as some general properties on the interactions governing folding dynamics. Proteins 2013; 81:1200–1211. © 2013 Wiley Periodicals, Inc.     

Coarse-grained strategy for modeling protein stability in concentrated solutions

We present a coarse-grained approach for modeling the thermodynamic stability of single-domain globular proteins in concentrated aqueous solutions. Our treatment derives effective protein-protein interactions from basic structural and energetic characteristics of the native and denatured states. These characteristics, along with the intrinsic (i.e., infinite dilution) thermodynamics of folding, are calculated from elementary sequence information using a heteropolymer collapse theory. We integrate this information into Reactive Canonical Monte Carlo simulations to investigate the connections between protein sequence hydrophobicity, protein-protein interactions, protein concentration, and the thermodynamic stability of the native state. The model predicts that sequence hydrophobicity can affect how protein concentration impacts native-state stability in solution. In particular, low hydrophobicity proteins are primarily stabilized by increases in protein concentration, whereas high hydrophobicity proteins exhibit richer nonmonotonic behavior. These trends appear qualitatively consistent with the available experimental data. Although factors such as pH, salt concentration, and protein charge are also important for protein stability, our analysis suggests that some of the nontrivial experimental trends may be driven by a competition between destabilizing hydrophobic protein-protein attractions and entropic crowding effects.     

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

A Simple Lattice Model That Captures Protein Folding,Aggregation and Amyloid Formation

The ability of many proteins to convert from their functional soluble state to amyloid fibrils can be attributed to inter-molecular beta strand formation. Such amyloid formation is associated with neurodegenerative disorders like Alzheimer''s and Parkinson''s. Molecular modelling can play a key role in providing insight into the factors that make proteins prone to fibril formation. However, fully atomistic models are computationally too expensive to capture the length and time scales associated with fibril formation. As the ability to form fibrils is the rule rather than the exception, much insight can be gained from the study of coarse-grained models that capture the key generic features associated with amyloid formation. Here we present a simple lattice model that can capture both protein folding and beta strand formation. Unlike standard lattice models, this model explicitly incorporates the formation of hydrogen bonds and the directionality of side chains. The simplicity of our model makes it computationally feasible to investigate the interplay between folding, amorphous aggregation and fibril formation, and maintains the capability of classic lattice models to simulate protein folding with high specificity. In our model, the folded proteins contain structures that resemble naturally occurring beta-sheets, with alternating polar and hydrophobic amino acids. Moreover, fibrils with intermolecular cross-beta strand conformations can be formed spontaneously out of multiple short hydrophobic peptide sequences. Both the formation of hydrogen bonds in folded structures and in fibrils is strongly dependent on the amino acid sequence, indicating that hydrogen-bonding interactions alone are not strong enough to initiate the formation of beta sheets. This result agrees with experimental observations that beta sheet and amyloid formation is strongly sequence dependent, with hydrophobic sequences being more prone to form such structures. Our model should open the way to a systematic study of the interplay between the factors that lead to amyloid formation.     

Thermodynamics of protein denatured states

Recent work on the thermodynamics of protein denatured states is providing insight into the stability of residual structure and the conformational constraints that affect the disordered states of proteins. Current data from native state hydrogen exchange and the pH dependence of protein stability indicate that residual structure can modulate the stability of the denatured state by up to 4 kcal mol(-1). NMR structural data have emphasized the role of hydrophobic clusters in stabilizing denatured state residual structures, however recent results indicate that electrostatic interactions, both favorable and unfavorable, are also important modulators of the stability of the denatured state. Thermodynamics methods that take advantage of histidine-heme ligation chemistry have also been developed to probe the conformational constraints that act on denatured states. These methods have provided insights into the role of excluded volume, chain stiffness, and loop persistence in modulating the conformational preferences of highly disordered proteins. New insights into protein folding and novel methods to manipulate protein stability are emerging from this work.     

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.     

Diffusion-collision model study of misfolding in a four-helix bundle protein.

Proteins with complex folding kinetics will be susceptible to misfolding at some stage in the folding process. We simulate this problem by using the diffusion-collision model to study non-native kinetic intermediate misfolding in a four-helix bundle protein. We find a limit on the size of the pairwise hydrophobic area loss in non-native intermediates, such that burying above this limit creates long-lasting non-native kinetic intermediates that would disrupt folding and prevent formation of the native state. Our study of misfolding suggests a method for limiting the production of misfolded kinetic intermediates for helical proteins and could, perhaps, lead to more efficient production of proteins in bulk.     

Direct analysis of backbone-backbone hydrogen bond formation in protein folding transition states

Here we investigate the role of backbone-backbone hydrogen bonding interactions in stabilizing the protein folding transition states of two model protein systems, the B1 domain of protein L (ProtL) and the P22 Arc repressor. A backbone modified analogue of ProtL containing an amide-to-ester bond substitution between residues 105 and 106 was prepared by total chemical synthesis, and the thermodynamic and kinetic parameters associated with its folding reaction were evaluated. Ultimately, these parameters were used in a Phi-value analysis to determine if the native backbone-backbone hydrogen bonding interaction perturbed in this analogue (i.e. a hydrogen bond in the first beta-turn of ProtL's beta-beta-alpha-beta-beta fold) was formed in the transition state of ProtL's folding reaction. Also determined were the kinetic parameters associated with the folding reactions of two Arc repressor analogues, each containing an amide-to-ester bond substitution in the backbone of their polypeptide chains. These parameters were used together with previously established thermodynamic parameters for the folding of these analogues in Phi-value analyses to determine if the native backbone-backbone hydrogen bonding interactions perturbed in these analogues (i.e. a hydrogen bond at the end of the intersubunit beta-sheet interface and hydrogen bonds at the beginning of the second alpha-helix in Arc repressor's beta-alpha-alpha structure) were formed in the transition state of Arc repressor's folding reaction. Our results reveal that backbone-backbone hydrogen bonding interactions are formed in the beta-turn and alpha-helical transition state structures of ProtL and Arc repressor, respectively; and they were not formed in the intersubunit beta-sheet interface of Arc repressor, a region of Arc repressor's polypeptide chain previously shown to have other non-native-like conformations in Arc's protein folding transition state.     

The stability and folding process of amyloidogenic mutant human lysozymes.

Amyloid deposits are frequently formed by mutant proteins that have a lower stability than the wild-type proteins. Some reports, however, have shown that mutant-induced thermodynamic destabilization is not always a general mechanism of amyloid formation. To obtain a better understanding of the mechanism of amyloid fibril formation, we show in this study that equilibrium and kinetic refolding-unfolding reaction experiments with two amyloidogenic mutant human lysozymes (I56T and D67H) yield folding pathways that can be drawn as Gibbs energy diagrams. The equilibrium stabilities between the native and denatured states of both mutant proteins were decreased, but the degrees of instability were different. The Gibbs energy diagrams of the folding process reveal that the Gibbs energy change between the native and folding intermediate states was similar for both proteins, and also that the activation Gibbs energy change from the native state to the transition state decreased. Our results confirm that the tendency to favor the intermediate of denaturation facilitates amyloid formation by the mutant human lysozymes more than equilibrium destabilization between the native and completely denatured states does.     

Effective interactions between chaotropic agents and proteins

Chaotropic agents are cosolutes that can disrupt the hydrogen bonding network between water molecules and reduce the stability of the native state of proteins by weakening the hydrophobic effect. In this work, we represent the chaotropic agent as a factor that reduces the amount of order in the structures formed by water molecules, both in the bulk and the hydration shells around hydrophobic amino acids. In this framework we show that low chaotrope concentrations lead to a destabilization of the native state of proteins, and that high concentrations induce complete denaturation. We also find that the reduction of the number of bulk ordered states of water molecules can give origin to an effective interaction between chaotropic molecules and proteins.     

不同类型氨基酸网络参量与蛋白质折叠的关系

理论和实验研究表明,蛋白质天然拓扑结构对其折叠过程具有重要的影响．采用复杂网络的方法分析蛋白质天然结构的拓扑特征,并探索蛋白质结构特征与折叠速率之间的内在联系．分别构建了蛋白质氨基酸网络、疏水网、亲水网、亲水-疏水网以及相应的长程网络,研究了这些网络的匹配系数(assortativity coefficient)和聚集系数(clustering coefficient)的统计特性．结果表明,除了亲水-疏水网,上述各网络的匹配系数均为正值,并且氨基酸网和疏水网的匹配系数与折叠速率表现出明显的线性正相关,揭示了疏水残基间相互作用的协同性有助于蛋白质的快速折叠．同时,研究发现疏水网的聚集系数与折叠速率有明显的线性负相关关系,这表明疏水残基间三角结构(triangle construction)的形成不利于蛋白质快速折叠．还进一步构建了相应的长程网络,发现序列上间距较远的残基接触对的形成将使蛋白质折叠进程变慢．     

Forces stabilizing proteins

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site-directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a –CH₂− group on folding contributes 1.1 ± 0.5 kcal/mol to protein stability. (2) The burial of non-polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1 ± 0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.     

Domain interactions direct misfolding and amyloid formation of yeast phosphoglycerate kinase

There are proteins that are built of two structural domains and are deposited full-length in amyloid plaques formed in various diseases. In spite of the known differences in the mechanisms of folding of single- and multidomain proteins, no published studies can be found that address the role of the domain-domain interactions during misfolding and amyloid formation. By the discovery of the role of domain-domain interactions, here we provide important insight in the submolecular mechanism of amyloid formation. A model system based on yeast phosphoglycerate kinase was designed. This system includes the wild-type yeast phosphoglycerate kinase and single-tryptophan mutants of the individual N and C terminal domains and the complete protein. Electron microscopic measurements proved that amyloid fibrils grow from all mutants under identical conditions as for the wild-type protein. Misfolding and amyloid formation was followed in stopped-flow and manual mixing experiments on the 1 ms to 4 days timescale. Tryptophan fluorescence was used for selective detection of conformational changes accompanying the formation of the amyloidogenic intermediates and the growth of amyloid fibrils. The interactions between the polypeptide chains of the two domains direct the misfolding process from the early steps to the amyloid formation, and influence the final structure. The kinetics of misfolding is different for the individual domains, pointing to the significance of the amino acid sequence. Misfolding of the domains within the complete protein is synchronized indicating that domain-domain interactions direct the misfolding and amyloid formation mechanism.