What is the role of thermodynamics on protein stability?

The most challenging and emerging field of biotechnology is the tailoring of proteins to attain the desired characteristic properties. In order to increase the stability of proteins and to study the function of proteins, the mechanism by which proteins fold and unfold should be known. It has been debated for a long time how exactly the linear form of a protein is converted into a stable 3-dimensional structure. The literature showed that many theories support the fact that protein folding is a thermodynamically controlled process. It is also possible to predict the mechanism of protein deactivation and stability to an extent from thermodynamic studies. This article reviewed various theories that have been proposed to explain the process of protein folding after its biosynthesis in ribosomes. The theories of the determination of the thermodynamic properties and the interpretation of thermodynamic data of protein stability are also discussed in this article.     

A model study of protein nascent chain and cotranslational folding using hydrophobic-polar residues

To study protein nascent chain folding during biosynthesis, we investigate the folding behavior of models of hydrophobic and polar (HP) chains at growing length using both two-dimensional square lattice model and an optimized three-dimensional 4-state discrete off-lattice model. After enumerating all possible sequences and conformations of HP heteropolymers up to length N = 18 and N = 15 in two and three-dimensional space, respectively, we examine changes in adopted structure, stability, and tolerance to single point mutation as the nascent chain grows. In both models, we find that stable model proteins have fewer folded nascent chains during growth, and often will only fold after reaching full length. For the few occasions where partial chains of stable proteins fold, these partial conformations on average are very similar to the corresponding parts of the final conformations at full length. Conversely, we find that sequences with fewer stable nascent chains and sequences with native-like folded nascent chains are more stable. In addition, these stable sequences in general can have many more point mutations and still fold into the same conformation as the wild type sequence. Our results suggest that stable proteins are less likely to be trapped in metastable conformations during biosynthesis, and are more resistant to point-mutations. Our results also imply that less stable proteins will require the assistance of chaperone and other factors during nascent chain folding. Taken together with other reported studies, it seems that cotranslational folding may not be a general mechanism of in vivo protein folding for small proteins, and in vitro folding studies are still relevant for understanding how proteins fold biologically.     

A framework for describing topological frustration in models of protein folding

In a natively folded protein of moderate or larger size, the protein backbone may weave through itself in complex ways, raising questions about what sequence of events might have to occur in order for the protein to reach its native configuration from the unfolded state. A mathematical framework is presented here for describing the notion of a topological folding barrier, which occurs when a protein chain must pass through a hole or opening, formed by other regions of the protein structure. Different folding pathways encounter different numbers of such barriers and therefore different degrees of frustration. A dynamic programming algorithm finds the optimal theoretical folding path and minimal degree of frustration for a protein based on its natively folded configuration. Calculations over a database of protein structures provide insights into questions such as whether the path of minimal frustration might tend to favor folding from one or from many sites of folding nucleation, or whether proteins favor folding around the N terminus, thereby providing support for the hypothesis that proteins fold co-translationally. The computational methods are applied to a multi-disulfide bonded protein, with computational findings that are consistent with the experimentally observed folding pathway. Attention is drawn to certain complex protein folds for which the computational method suggests there may be a preferred site of nucleation or where folding is likely to proceed through a relatively well-defined pathway or intermediate. The computational analyses lead to testable models for protein folding.     

Folding studies of immunoglobulin-like beta-sandwich proteins suggest that they share a common folding pathway.

BACKGROUND: Are folding pathways conserved in protein families? To test this explicitly and ask to what extent structure specifies folding pathways requires comparison of proteins with a common fold. Our strategy is to choose members of a highly diverse protein family with no conservation of function and little or no sequence identity, but with structures that are essentially the same. The immunoglobulin-like fold is one of the most common structural families, and is subdivided into superfamilies with no detectable evolutionary or functional relationship. RESULTS: We compared the folding of a number of immunoglobulin-like proteins that have a common structural core and found a strong correlation between folding rate and stability. The results suggest that the folding pathways of these immunoglobulin-like proteins share common features. CONCLUSIONS: This study is the first to compare the folding of structurally related proteins that are members of different superfamilies. The most likely explanation for the results is that interactions that are important in defining the structure of immunoglobulin-like proteins are also used to guide folding.     

Right‐ and left‐handed three‐helix proteins. I. Experimental and simulation analysis of differences in folding and structure

Despite the large number of publications on three‐helix protein folding, there is no study devoted to the influence of handedness on the rate of three‐helix protein folding. From the experimental studies, we make a conclusion that the left‐handed three‐helix proteins fold faster than the right‐handed ones. What may explain this difference? An important question arising in this paper is whether the modeling of protein folding can catch the difference between the protein folding rates of proteins with similar structures but with different folding mechanisms. To answer this question, the folding of eight three‐helix proteins (four right‐handed and four left‐handed), which are similar in size, was modeled using the Monte Carlo and dynamic programming methods. The studies allowed us to determine the orders of folding of the secondary‐structure elements in these domains and amino acid residues which are important for the folding. The obtained data are in good correlation with each other and with the experimental data. Structural analysis of these proteins demonstrated that the left‐handed domains have a lesser number of contacts per residue and a smaller radius of cross section than the right‐handed domains. This may be one of the explanations of the observed fact. The same tendency is observed for the large dataset consisting of 332 three‐helix proteins (238 right‐ and 94 left‐handed). From our analysis, we found that the left‐handed three‐helix proteins have some less‐dense packing that should result in faster folding for some proteins as compared to the case of right‐handed proteins.Proteins 2013; © 2013 Wiley Periodicals, Inc.     

Distinguishing between sequential and nonsequentially folded proteins: implications for folding and misfolding.

We describe here an algorithm for distinguishing sequential from nonsequentially folding proteins. Several experiments have recently suggested that most of the proteins that are synthesized in the eukaryotic cell may fold sequentially. This proposed folding mechanism in vivo is particularly advantageous to the organism. In the absence of chaperones, the probability that a sequentially folding protein will misfold is reduced significantly. The problem we address here is devising a procedure that would differentiate between the two types of folding patterns. Footprints of sequential folding may be found in structures where consecutive fragments of the chain interact with each other. In such cases, the folding complexity may be viewed as being lower. On the other hand, higher folding complexity suggests that at least a portion of the polypeptide backbone folds back upon itself to form three-dimensional (3D) interactions with noncontiguous portion(s) of the chain. Hence, we look at the mechanism of folding of the molecule via analysis of its complexity, that is, through the 3D interactions formed by contiguous segments on the polypeptide chain. To computationally splice the structure into consecutively interacting fragments, we either cut it into compact hydrophobic folding units or into a set of hypothetical, transient, highly populated, contiguous fragments ("building blocks" of the structure). In sequential folding, successive building blocks interact with each other from the amino to the carboxy terminus of the polypeptide chain. Consequently, the results of the parsing differentiate between sequentially vs. nonsequentially folded chains. The automated assessment of the folding complexity provides insight into both the likelihood of misfolding and the kinetic folding rate of the given protein. In terms of the funnel free energy landscape theory, a protein that truly follows the mechanism of sequential folding, in principle, encounters smoother free energy barriers. A simple sequentially folded protein should, therefore, be less error prone and fold faster than a protein with a complex folding pattern.     

Matching simulation and experiment: a new simplified model for simulating protein folding.

Simulations of simplified protein folding models have provided much insight into solving the protein folding problem. We propose here a new off-lattice bead model, capable of simulating several different fold classes of small proteins. We present the sequence for an alpha/beta protein resembling the IgG-binding proteins L and G. The thermodynamics of the folding process for this model are characterized using the multiple multihistogram method combined with constant-temperature Langevin simulations. The folding is shown to be highly cooperative, with chain collapse nearly accompanying folding. Two parallel folding pathways are shown to exist on the folding free energy landscape. One pathway contains an intermediate--similar to experiments on protein G, and one pathway contains no intermediates-similar to experiments on protein L. The folding kinetics are characterized by tabulating mean-first passage times, and we show that the onset of glasslike kinetics occurs at much lower temperatures than the folding temperature. This model is expected to be useful in many future contexts: investigating questions of the role of local versus nonlocal interactions in various fold classes, addressing the effect of sequence mutations affecting secondary structure propensities, and providing a computationally feasible model for studying the role of solvation forces in protein folding.     

The posttranslational concept of protein folding: how valid is it?

Two concepts of protein folding are known. One of them, the cotranslational concept, states that a protein folds during the synthesis of the polypeptide chain on the ribosome. According to the other, the posttranslational concept, the protein starts to fold just after the synthesis of its polypeptide chain. This article attempts to show that the posttranslational concept is hardly suited to solve the problem of protein folding. In our opinion, polypeptide chains cannot be represented as random coils. They are stiff chain-like macromolecules rather than flexible ones: the single bond rotational barriers of a polypeptide substantially exceed the accepted standard values; even in strong denaturing conditions, a protein possesses a considerable amount of residual folded structures. We believe that the popular "hierarchical" models for the protein folding mechanism are not realistic because the formation of secondary and tertiary structures of proteins occurs simultaneously and cooperatively. The time for the elongation of a polypeptide chain by one amino acid residue during biosynthesis exceeds considerably the time of the formation of alpha-helices and beta-sheets in proteins as well as the time supposed for the spatial structure formation of a native protein during renaturation. Thus, we believe that the mechanism of protein folding in vivo cannot be clarified by denaturation-renaturation experiments. In our opinion, the phenomenon of protein renaturation is no more than the restoration of native protein conformation (which initially forms cotranslationally) disrupted during denaturation, and thus denaturation-renaturation experiments cannot serve as a model to clarify the mechanism of protein folding.     

Analyses of the folding properties of ferredoxin‐like fold proteins by means of a coarse‐grained Gō model: Relationship between the free energy profiles and folding cores

The folding mechanisms of proteins with multi‐state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on‐going research issues. In this study, we investigate ferredoxin‐like fold proteins which have a simple topology and multi‐state transitions. We analyze the folding processes by means of a coarse‐grained Gō model. We are able to reproduce the differences in the folding mechanisms between U1A, which has a high‐free‐energy intermediate state, and ADA2h and S6, which fold into the native structure through two‐state transitions. The folding pathways of U1A, ADA2h, S6, and the S6 circular permutant, S6_p54‐55, are reproduced and compared with experimental observations. We show that the ferredoxin‐like fold contains two common regions consisting folding cores as predicted in other studies and that U1A produces an intermediate state due to the distinct cooperative folding of each core. However, because one of the cores of S6 loses its cooperativity and the two cores of ADA2h are tightly coupled, these proteins fold into the native structure through a two‐state mechanism. Proteins 2014; 82:954–965. © 2013 Wiley Periodicals, Inc.     

How the folding rates of two‐ and multistate proteins depend on the amino acid properties

Proteins fold by either two‐state or multistate kinetic mechanism. We observe that amino acids play different roles in different mechanism. Many residues that are easy to form regular secondary structures (α helices, β sheets and turns) can promote the two‐state folding reactions of small proteins. Most of hydrophilic residues can speed up the multistate folding reactions of large proteins. Folding rates of large proteins are equally responsive to the flexibility of partial amino acids. Other properties of amino acids (including volume, polarity, accessible surface, exposure degree, isoelectric point, and phase transfer energy) have contributed little to folding kinetics of the proteins. Cysteine is a special residue, it triggers two‐state folding reaction and but inhibits multistate folding reaction. These findings not only provide a new insight into protein structure prediction, but also could be used to direct the point mutations that can change folding rate. Proteins 2014; 82:2375–2382. © 2014 Wiley Periodicals, Inc.     

First principles prediction of protein folding rates

Experimental studies have demonstrated that many small, single-domain proteins fold via simple two-state kinetics. We present a first principles approach for predicting these experimentally determined folding rates. Our approach is based on a nucleation-condensation folding mechanism, where the rate-limiting step is a random, diffusive search for the native tertiary topology. To estimate the rates of folding for various proteins via this mechanism, we first determine the probability of randomly sampling a conformation with the native fold topology. Next, we convert these probabilities into folding rates by estimating the rate that a protein samples different topologies during diffusive folding. This topology-sampling rate is calculated using the Einstein diffusion equation in conjunction with an experimentally determined intra-protein diffusion constant. We have applied our prediction method to the 21 topologically distinct small proteins for which two-state rate data is available. For the 18 beta-sheet and mixed alpha-beta native proteins, we predict folding rates within an average factor of 4, even though the experimental rates vary by a factor of approximately 4 x 10(4). Interestingly, the experimental folding rates for the three four-helix bundle proteins are significantly underestimated by this approach, suggesting that proteins with significant helical content may fold by a faster, alternative mechanism. This method can be applied to any protein for which the structure is known and hence can be used to predict the folding rates of many proteins prior to experiment.     

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.     

The protein folding problem: when will it be solved?

The protein folding problem can be viewed as three different problems: defining the thermodynamic folding code; devising a good computational structure prediction algorithm; and answering Levinthal's question regarding the kinetic mechanism of how proteins can fold so quickly. Once regarded as a grand challenge, protein folding has seen much progress in recent years. Folding codes are now being used to successfully design proteins and non-biological foldable polymers; aided by the Critical Assessment of Techniques for Structure Prediction (CASP) competition, protein structure prediction has now become quite good. Even the once-challenging Levinthal puzzle now seems to have an answer--a protein can avoid searching irrelevant conformations and fold quickly by making local independent decisions first, followed by non-local global decisions later.     

Methods to identify the substrates of thiol‐disulfide oxidoreductases

The formation of a disulfide bond is a critical step in the folding of numerous secretory and membrane proteins and catalyzed in vivo. A variety of mechanisms and protein structures have evolved to catalyze oxidative protein folding. Those enzymes that directly interact with a folding protein to accelerate its oxidative folding are mostly thiol‐disulfide oxidoreductases that belong to the thioredoxin superfamily. The enzymes of this class often use a CXXC active‐site motif embedded in their thioredoxin‐like fold to promote formation, isomerization, and reduction of a disulfide bond in their target proteins. Over the past decade or so, an increasing number of substrates of the thiol‐disulfide oxidoreductases that are present in the ER of mammalian cells have been discovered, revealing that the enzymes play unexpectedly diverse physiological functions. However, functions of some of these enzymes still remain unclear due to the lack of information on their substrates. Here, we review the methods used by researchers to identify the substrates of these enzymes and provide data that show the importance of using trichloroacetic acid in sample preparation for the substrate identification, hoping to aid future studies. We particularly focus on successful studies that have uncovered physiological substrates and functions of the enzymes in the periplasm of Gram‐negative bacteria and the endoplasmic reticulum of mammalian cells. Similar approaches should be applicable to enzymes in other cellular compartments or in other organisms.     

Barnase and barstar: two small proteins to fold and fit together

Barnase and barstar are the extracellular ribonuclease and its intracellular inhibitor produced by Bacillus amyloliquefaciens. Both are small single-chain proteins and thus are suitable for application to the study of how a protein's sequence directs its fold. Barnase has neither disulfide bonds nor non-peptide components and unfolds reversibly in what closely approximates a two-state reaction. The genes for both these proteins have been cloned in E. coli. Expression of barstar is necessary to counter the lethal effect of expressed active barnase. Site-directed mutagenesis is being used to answer specific and general questions relating to protein folding and protein-protein interaction.     

Exposing the kinetic traps in RNA folding.

Large ribozymes fold on a 'glacial' timescale compared to the folding of their protein counterparts. The sluggish folding exhibited by large RNAs results from the formation of kinetically trapped, misfolded intermediates, which are nonessential features of the folding mechanism. Newly developed mutant ribozymes that avoid kinetic traps should facilitate the study of the RNA folding problem.     

The In Vivo Association of BiP with Newly Synthesized Proteins Is Dependent on the Rate and Stability of Folding and Not Simply on the Presence of Sequences That Can Bind to BiP

Immunoglobulin heavy chain-binding protein (BiP) is a member of the hsp70 family of chaperones and one of the most abundant proteins in the ER lumen. It is known to interact transiently with many nascent proteins as they enter the ER and more stably with protein subunits produced in stoichiometric excess or with mutant proteins. However, there also exists a large number of secretory pathway proteins that do not apparently interact with BiP. To begin to understand what controls the likelihood that a nascent protein entering the ER will associate with BiP, we have examined the in vivo folding of a murine λI immunoglobulin (Ig) light chain (LC). This LC is composed of two Ig domains that can fold independent of the other and that each possess multiple potential BiP-binding sequences. To detect BiP binding to the LC during folding, we used BiP ATPase mutants, which bind irreversibly to proteins, as “kinetic traps.” Although both the wild-type and mutant BiP clearly associated with the unoxidized variable region domain, we were unable to detect binding of either BiP protein to the constant region domain. A combination of in vivo and in vitro folding studies revealed that the constant domain folds rapidly and stably even in the absence of an intradomain disulfide bond. Thus, the simple presence of a BiP-binding site on a nascent chain does not ensure that BiP will bind and play a role in its folding. Instead, it appears that the rate and stability of protein folding determines whether or not a particular site is recognized, with BiP preferentially binding to proteins that fold slowly or somewhat unstably.     

A comparison of the folding of two knotted proteins: YbeA and YibK

The extraordinary topology of proteins belonging to the alpha/beta-knot superfamily of proteins is unexpected, due to the apparent complexities involved in the formation of a deep trefoil knot in a polypeptide backbone. Despite this, an increasing number of knotted structures are being identified; how such proteins fold remains a mystery. Studies on the dimeric protein YibK from Haemophilus influenzae have led to the characterisation of its folding pathway in some detail. To complement research into the folding of YibK, and to address whether folding pathways are conserved for members of the alpha/beta-knot superfamily, the structurally similar knotted protein YbeA from Escherichia coli has been studied. A comprehensive thermodynamic and kinetic analysis of the folding of YbeA is presented here, and compared to that of YibK. Both fold via an intermediate state populated under equilibrium conditions that is monomeric and considerably structured. The unfolding/refolding kinetics of YbeA are simpler than those found for YibK and involve two phases attributed to the formation of a monomeric intermediate state and a dimerisation step. In contrast to YibK, a change in the rate-determining step on the unfolding pathway for YbeA is observed with a changing concentration of urea. Despite this difference, both proteins fold by a mechanism involving at least one sequential monomeric intermediate that has properties similar to that observed during the equilibrium unfolding. The rate of dimerisation observed for YbeA and YibK is very similar, as is the rate constant for formation of the kinetic monomeric intermediate that precedes dimerisation. The findings suggest that relatively slow folding and dimerisation may be common attributes of knotted proteins.     

Chaperonin 60: a paradoxical,evolutionarily conserved protein family with multiple moonlighting functions

Chaperonin 60 is the prototypic molecular chaperone, an essential protein in eukaryotes and prokaryotes, whose sequence conservation provides an excellent basis for phylogenetic analysis. Escherichia coli chaperonin 60 (GroEL), the prototype of this family of proteins, has an established oligomeric‐structure‐based folding mechanism and a defined population of folding partners. However, there is a growing number of examples of chaperonin 60 proteins whose crystal structures and oligomeric composition are at variance with GroEL, suggesting that additional complexities in the protein‐folding function of this protein should be expected. In addition, many organisms have multiple chaperonin 60 proteins, some of which have lost their protein‐folding ability. It is emerging that this highly conserved protein has evolved a bewildering variety of additional biological functions – known as moonlighting functions – both within the cell and in the extracellular milieu. Indeed, in some organisms, it is these moonlighting functions that have been left after the loss of the protein‐folding activity. This highlights the major paradox in the biology of chaperonin 60. This article reviews the relationship between the folding and non‐folding (moonlighting) activities of the chaperonin 60 family and discusses current knowledge on their molecular evolution focusing on protein domains involved in the non‐folding chaperonin functions in an attempt to understand the emerging biology of this evolutionarily ancient protein family.