Cunning simplicity of a hierarchical folding

A hierarchic scheme of protein folding does not solve the Levinthal paradox since it cannot provide a simultaneous explanation for major features observed for protein folding: (i) folding within non-astronomical time, (ii) independence of the native structure on large variations in the folding rates of given protein under different conditions, and (iii) co-existence, in a visible quantity, of only the native and the unfolded molecules during folding of moderate size (single-domain) proteins. On the contrary, a nucleation mechanism can account for all these major features simultaneously and resolves the Levinthal paradox.     

Spelling Protein Structure

Abstract

We would be tempted to state that there has never been a Levinthal paradox. Indeed, Levinthal raised an interesting problem about protein folding, as he realized that proteins have no time to explore exhaustively their conformational space on the way to their native structure. He did not seem to find this paradoxical and immediately proposed a straightforward solution, which has essentially never been refuted. In other words, Levinthal solved his own paradox.     

Loop Fold Structure of Proteins: Resolution of Levinthal's Paradox

Abstract

According to Levinthal a protein chain of ordinary size would require enormous time to sort its conformational states before the final fold is reached. Experimentally observed time of folding suggests an estimate of the chain length for which the time would be sufficient. This estimate by order of magnitude fits to experimentally observed universal closed loop elements of globular proteins—25–30 residues.     

Learning about protein folding via potential functions

Over the last few years we have developed an empirical potential function that solves the protein structure recognition problem: given the sequence for an n-residue globular protein and a collection of plausible protein conformations, including the native conformation for that sequence, identify the correct, native conformation. Having determined this potential on the basis of only some 6500 native/nonnative pairs of structures for 58 proteins, we find it recognizes the native conformation for essentially all compact, soluble, globular proteins having known native conformations in comparisons with 10⁴ to 10⁶ reasonable alternative conformations apiece. In this sense, the potential encodes nearly all the essential features of globular protein conformational preference. In addition it “knows” about many additional factors in protein folding, such as the stabilization of multimeric proteins, quaternary structure, the role of disulfide bridges and ligands, proproteins vs. processed proteins, and minimal strand lengths in globular proteins. Comparisons are made with other sorts of protein folding problems, and applications in protein conformational determination and prediction are discussed. © 1994 Wiley-Liss, Inc.     

Structure and heterogeneity of the one- and two-disulfide folding intermediates of tick anticoagulant peptide

Tick anticoagulant peptide (TAP) is a factor Xa-specific inhibitor and is structurally homologous to bovine pancreatic trypsin inhibitor (BPTI). The fully reduced TAP refolds spontaneously to form the native structure under a wide variation of redox buffers. The folding intermediates of TAP consist of at least 22 fractions of one-disulfide, two-disulfide, and three-disulfide scrambled isomers. Three species of well-populated one- and two-disulfide intermediates were isolated and structurally characterized. The predominant one-disulfide species contains TAP-(Cys33—Cys55). Two major two-disulfide isomers were TAP-(Cys33—Cys55, Cys15—Cys39) and TAP-(Cys33—Cys55, Cys5—Cys39). Both Cys33—Cys55 and Cys15—Cys39 are native disulfides of TAP. These three species are structural counterparts of BPTI-(Cys30—Cys51), BPTI-(Cys30—Cys51, Cys14—Cys38), and BPTI-(Cys30—Cys51,Cys5—Cys38), which have been shown to be the major intermediates of BPTI folding. In addition, time-course-trapped folding intermediates of TAP, consisting of about 47% one-disulfide species and 30% two-disulfide species, were collectively digested with thermolysin, and fragmented peptides were analyzed by Edman sequencing and mass spectrometry in order to characterize the disulfide-containing peptides. Among the 15 possible single-disulfide pairings of TAP, 10 (2 native and 8 nonnative) were found as structural components of its one- and two-disulfide folding intermediates. The results demonstrate that the major folding intermediates of TAP bear structural homology to those of BPTI. However, the folding pathway of TAP differs from that of BPTI by (a) a higher degree of heterogeneity of one- and two-disulfide intermediates and (b) the presence of three-disulfide scrambled isomers as folding intermediates. Mechanism(s) that may account for these diversities are proposed and discussed.     

Back to Units of Protein Folding

Abstract

In response to the criticism by A. Finkelstein (J. Biomol. Struct. Dyn. 20, 311–314, 2002) of our Communication (J. Biomol. Struct. Dyn. 20, 5–6, 2002) several issues are dealt with. Importance of the notion of elementary folding unit, its size and structure, and the necessity of further characterization of the units for the elucidation of the protein folding in vivo are discussed. The criticism (J. Biomol. Struct. Dyn. 20, 311–314, 2002) on the hierarchical protein folding is also briefly addressed.     

Pathway of Oxidative Folding of a 3-Disulfide α-Lactalbumin May Resemble Either BPTI Model or Hirudin Model

Pathways of oxidative folding of disulfide proteins display a high degree of diversity and vary among two extreme models. The BPTI model is defined by limited species of folding intermediates adopting mainly native disulfide bonds. The hirudin model is characterized by highly heterogeneous folding intermediates containing mostly non-native disulfide bonds. αLA-IIIA is a 3-disulfide variant of α-lactalbumin (αLA) with a 3-D conformation essentially identical to that of intact αLA. αLA-IIIA contains 3 native disulfide bonds of αLA, two of them are located at the calcium binding β-subdomain (Cys⁶¹–Cys⁷⁷ and Cys⁷³–Cys⁹¹) and the third bridge is located within the α-helical domain of the molecule (Cys²⁸–Cys¹¹¹). We investigate here the pathway of oxidative folding of fully reduced αLA-IIIA with and without stabilization of its β-subdomain by calcium binding. In the absence of calcium, the folding pathway of αLA-IIIA was shown to resemble that of hirudin model. Upon stabilization of β-sheet domain by calcium binding, the folding pathway of αLA-IIIA exhibits a striking similarity to that of BPTI model. Three predominant folding intermediates of αLA-IIIA containing exclusively native disulfide bonds were isolated and structurally characterized. Our results further demonstrate that stabilization of subdomains in a protein may dictate its folding pathway and represent a major cause for the existing diversity in the folding pathways of the disulfide-containing proteins.     

Recent Advances in GFP Folding Reporter and Split-GFP Solubility Reporter Technologies. Application to Improving the Folding and Solubility of Recalcitrant Proteins from Mycobacterium tuberculosis

We have improved our green fluorescent protein (GFP) folding reporter technology [Waldo et al., (1999) Nat. Biotechnol. 17, 691–695] to evolve recalcitrant proteins from Mycobacterium tuberculosis. The target protein is inserted into the scaffolding of the GFP, eliminating false-positive artifacts caused by expression of truncated protein variants from internal cryptic ribosome binding sites in the target RNA. In parallel, we have developed a new quantitative fluorescent protein tagging and detection system based on micro-domains of GFP. This split-GFP system, which works both in vivo and in vitro, is amenable to high-throughput assays of protein expression and solubility [Cabantous et al., (2005) Nat. Biotechnol. 23, 102–107]. Together, the GFP folding reporter and split-GFP technologies offer a comprehensive system for manipulating and improving protein folding and solubility.     

Structure of a partially unfolded form of Escherichia coli dihydrofolate reductase provides insight into its folding pathway

Proteins frequently fold via folding intermediates that correspond to local minima on the conformational energy landscape. Probing the structure of the partially unfolded forms in equilibrium under native conditions can provide insight into the properties of folding intermediates. To elucidate the structures of folding intermediates of Escherichia coli dihydrofolate reductase (DHFR), we investigated transient partial unfolding of DHFR under native conditions. We probed the structure of a high‐energy conformation susceptible to proteolysis (cleavable form) using native‐state proteolysis. The free energy for unfolding to the cleavable form is clearly less than that for global unfolding. The dependence of the free energy on urea concentration (m‐value) also confirmed that the cleavable form is a partially unfolded form. By assessing the effect of mutations on the stability of the partially unfolded form, we found that native contacts in a hydrophobic cluster formed by the F‐G and Met‐20 loops on one face of the central β‐sheet are mostly lost in the partially unfolded form. Also, the folded region of the partially unfolded form is likely to have some degree of structural heterogeneity. The structure of the partially unfolded form is fully consistent with spectroscopic properties of the near‐native kinetic intermediate observed in previous folding studies of DHFR. The findings suggest that the last step of the folding of DHFR involves organization in the structure of two large loops, the F‐G and Met‐20 loops, which is coupled with compaction of the rest of the protein.     

First-Principles Protein Folding Simulations

Abstract

It is widely believed that the prediction of the three-dimensional structures of proteins from the first principles is impossible. This view is based on the fact that the number of possible structures for each protein is astronomically large. The question is then why a protein folds into its native structure with the proper biological functions in the time scale of milliseconds to minutes, and this is called Levinthal's paradox. In this article I will discuss our strategy for attacking the protein folding problem. Our approach consists of two elements: the inclusion of accurate solvent effects and the development of powerful simulation algorithms that can avoid getting trapped in states of energy local minima. For the former, we discuss several models varying in nature from crude (distance-dependent dielectric function) to rigorous (reference interaction site model). For the latter, we show the effectiveness of Monte Carlo simulated annealing and generalized-ensemble algorithms.     

Loop Folds in Proteins and Evolutionary Conservation of Folding Nuclei

Abstract

We show that loops of close contacts involving hydrophobic residues are important in protein folding. Contrary to Berezovsky and Trifonov (J. Biomol. Struct. Dyn. 20, 5–6, 2002) the loops important in protein folding usually are much larger in size than 23–31 residues, being instead comparable to the size of the protein for single domain proteins. Additionally what is important are not single loop contacts, but a highly interconnected network of such loop contacts, which provides extra stability to a protein fold and which leads to their conservation in evolution.     

Peptide fragment studies on the folding elements of dihydrofolate reductase from Escherichia coli

One of the necessary conditions for a protein to be foldable is the presence of a complete set of “folding elements” (FEs) that are short, contiguous peptide segments distributed over an amino acid sequence. The FE‐assembly model of protein folding has been proposed, in which the FEs play a role in guiding structure formation through FE–FE interactions early in folding. However, two major issues remain to be clarified regarding the roles of the FEs in determining protein foldability. Are the FEs AFUs that can form nativelike structures in isolation? Is the presence of only the FEs without mutual connections a sufficient condition for a protein to be foldable? Here, we address these questions using peptide fragments corresponding to the FEs of DHFR from Escherichia coli. We show by CD measurement that the FE peptides are unfolded under the native conditions, and some of them have the propensities toward non‐native helices. MD simulations also show the non‐native helical propensities of the peptides, and the helix contents estimated from the simulations are well correlated with those estimated from the CD in TFE. Thus, the FEs of DHFR are not AFUs, suggesting the importance of the FEs in nonlocal interactions. We also show that equimolar mixtures of the FE peptides do not induce any structural formation. Therefore, mutual connections between the FEs, which should strengthen the nonlocal FE–FE interactions, are also one of the necessary conditions for a protein to be foldable. Proteins 2006. © 2005 Wiley‐Liss, Inc.     

In vitro insulin refolding: Characterization of the intermediates and the putative folding pathway

The in vitro refolding process of the double-chain insulin was studied based on the investigation of in vitro single-chain insulin refolding. Six major folding intermediates, named P1A, P2B, P3A, P4B, P5B, and P6B, were captured during the folding process. The refolding experiments indicate that all of these intermediates are on-pathway. Based on these intermediates and the formation of hypothetic transients, we propose a two-stage folding pathway of insulin. (1) At the early stage of the folding process, the reduced A chain and B chain individually formed the intermediates two A chain intermediates (P1A and P3A), and four B chain intermediates (P2B, P4B, P5B, and P6B). (2) In the subsequent folding process, transient Ⅰ was formed from P3A through thiol/disulfide exchange reaction; then, transients Ⅱ and Ⅲ, each containing two native disulfides, were formed through the recognition and interaction of transient Ⅰ with P4B or P6B and the thiol group's oxidation reaction mainly using GSSG as oxidative reagent; finally, transients Ⅱ and Ⅲ, through thiol/mixture disulfide exchange reaction, formed the third native disulfide of insulin to complete the folding.     

<Emphasis Type="Italic">In vitro</Emphasis> insulin refolding: Characterization of the intermediates and the putative folding pathway

The in vitro refolding process of the double-chain insulin was studied based on the investigation of in vitro single-chain insulin refolding. Six major folding intermediates, named P1A, P2B, P3A, P4B, P5B, and P6B, were captured during the folding process. The refolding experiments indicate that all of these intermediates are on-pathway. Based on these intermediates and the formation of hypothetic transients, we propose a two-stage folding pathway of insulin. (1) At the early stage of the folding process, the reduced A chain and B chain individually formed the intermediates: two A chain intermediates (P1A and P3A), and four B chain intermediates (P2B, P4B, P5B, and P6B). (2) In the subsequent folding process, transient I was formed from P3A through thiol/disulfide exchange reaction; then, transients II and III, each containing two native disulfides, were formed through the recognition and interaction of transient I with P4B or P6B and the thiol group’s oxidation reaction mainly using GSSG as oxidative reagent; finally, transients II and III, through thiol/mixture disulfide exchange reaction, formed the third native disulfide of insulin to complete the folding.     

On the mechanism of protein fold‐switching by a molecular sensor

Alternate frame folding (AFF) is a mechanism by which conformational change can be engineered into a protein. The protein structure switches from the wild‐type fold (N) to a circularly‐permuted fold (N′), or vice versa, in response to a signaling event such as ligand binding. Despite the fact that the two native states have similar structures, their interconversion involves folding and unfolding of large parts of the molecule. This rearrangement is reported by fluorescent groups whose relative proximities change as a result of the order–disorder transition. The nature of the conformational change is expected to be similar from protein to protein; thus, it may be possible to employ AFF as a general method to create optical biosensors. Toward that goal, we test basic aspects of the AFF mechanism using the AFF variant of calbindin D_9k. A simple three‐state model for fold switching holds that N and N′ interconvert through the unfolded state. This model predicts that the fundamental properties of the switch—calcium binding affinity, signal response (i.e., fluorescence change upon binding), and switching rate—can be controlled by altering the relative stabilities of N and N′. We find that selectively destabilizing N or N′ changes the equilibrium properties of the switch (binding affinity and signal response) in accordance with the model. However, kinetic data indicate that the switching pathway does not require whole‐molecule unfolding. The rate is instead limited by unfolding of a portion of the protein, possibly in concert with folding of a corresponding region. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Folding Nuclei in Proteins

When a protein folds or unfolds, it passes through many half-folded microstates. Only a few of them can accumulate and be seen experimentally, and this happens only when the folding (or unfolding) occurs far from the point of thermodynamic equilibrium between the native and denatured states. The universal features of folding, though, are observed in the vicinity of the equilibrium point. Here the two-state transition proceeds without any accumulation of metastable intermediates, and only the transition state (folding nucleus) is outlined by its key influence on the folding/unfolding kinetics. This review covers recent experimental and theoretical studies of folding nuclei.     

Acceleration of oxidative protein folding by curcumin through novel non-redox chemistry

Curcumin, the major constituent of turmeric is a known antioxidant. We have examined the oxidative folding of the model four-disulfide-bond-containing protein bovine pancreatic ribonuclease A (RNase A) in its presence; results indicate that RNase A regeneration rate increases in a curcumin-dependent manner. Examination of the native tendency of the fully-reduced polypeptide and the stability of key folding intermediates suggests that the increased oxidative folding rate can be attributed to native-like elements induced within the fully-reduced polypeptide and the stabilization of native-like species by this non-redox-active natural product. Our results provide a template for the design of curcuminoid-based synthetic small-molecule fold catalysts that accelerate the folding of ER-processed proteins; this assumes significance given that nitrosative stress and dysfunction of the ER-resident oxidoreductase protein disulfide isomerise due to S-nitrosylation are factors associated with the pathogenesis of Alzheimer’s and Parkinson’s diseases.     

The stability of Taq DNA polymerase results from a reduced entropic folding penalty; identification of other thermophilic proteins with similar folding thermodynamics

The thermal stability of Taq DNA polymerase is well known, and is the basis for its use in PCR. A comparative thermodynamic characterization of the large fragment domains of Taq (Klentaq) and E. coli (Klenow) DNA polymerases has been performed by obtaining full Gibbs‐Helmholtz stability curves of the free energy of folding (ΔG) versus temperature. This analysis provides the temperature dependencies of the folding enthalpy and entropy (ΔH and ΔS), and the heat capacity (ΔC_p) of folding. If increased or enhanced non‐covalent bonding in the native state is responsible for enhanced thermal stabilization of a protein, as is often proposed, then an enhanced favourable folding enthalpy should, in general, be observed for thermophilic proteins. However, for the Klenow–Klentaq homologous pair, the folding enthalpy (ΔH_fold) of Klentaq is considerably less favorable than that of Klenow at all temperatures. In contrast, it is found that Klentaq's extreme free energy of folding (ΔG_fold) originates from a significantly reduced entropic penalty of folding (ΔS_fold). Furthermore, the heat capacity changes upon folding are similar for Klenow and Klentaq. Along with this new data, comparable extended analysis of available thermodynamic data for 17 other mesophilic–thermophilic protein pairs (where enough applicable thermodynamic data exists) shows a similar pattern in seven of the 18 total systems. When analyzed with this approach, the more familiar “reduced ΔC_p mechanism” for protein thermal stabilization (observed in a different six of the 18 systems) frequently manifests as a temperature dependent shift from enthalpy driven stabilization to a reduced‐entropic‐penalty model. Proteins 2014; 82:785–793. © 2013 Wiley Periodicals, Inc.     

Protein folding ‐ seeing is deceiving

This Perspective is intended to raise questions about the conventional interpretation of protein folding. According to the conventional interpretation, developed over many decades, a protein population can visit a vast number of conformations under unfolding conditions, but a single dominant native population emerges under folding conditions. Accordingly, folding comes with a substantial loss of conformational entropy. How is this price paid? The conventional answer is that favorable interactions between and among the side chains can compensate for entropy loss, and moreover, these interactions are responsible for the structural particulars of the native conformation. Challenging this interpretation, the Perspective introduces a proposal that high energy (i.e., unfavorable) excluding interactions winnow the accessible population substantially under physical–chemical conditions that favor folding. Both steric clash and unsatisfied hydrogen bond donors and acceptors are classified as excluding interactions, so called because conformers with such disfavored interactions will be largely excluded from the thermodynamic population. Both excluding interactions and solvent factors that induce compactness are somewhat nonspecific, yet together they promote substantial chain organization. Moreover, proteins are built on a backbone scaffold consisting of α‐helices and strands of β‐sheet, where the number of hydrogen bond donors and acceptors is exactly balanced. These repetitive secondary structural elements are the only two conformers that can be both completely hydrogen‐bond satisfied and extended indefinitely without encountering a steric clash. Consequently, the number of fundamental folds is limited to no more than ~10,000 for a protein domain. Once excluding interactions are taken into account, the issue of “frustration” is largely eliminated and the Levinthal paradox is resolved. Putting the “bottom line” at the top: it is likely that hydrogen‐bond satisfaction represents a largely under‐appreciated parameter in protein folding models.     

162 Dominant forces in protein folding

In the first part of this talk, I will discuss the need for a paradigm shift from hydrophobic (HφO) to a hydrophilic ((HφI) based theory of protein folding. Next, I will discuss the various types of solvent-induced forces that are exerted on various groups on the protein. It is argued, both theoretically and by simulations, that the HφI–HφI solvent-induced forces are likely to be the strongest. Therefore, it is suggested that these forces are also the forces that force the protein to fold, in a short time, along a narrow range of pathways. This paradigm shift also answers Levinthal’s question about the factors that “speed” and “guide” the folding of proteins.