Simple model of protein folding kinetics.

A simple model of the kinetics of protein folding is presented. The reaction coordinate is the "correctness" of a configuration compared with the native state. The model has a gap in the energy spectrum, a large configurational entropy, a free energy barrier between folded and partially folded states, and a good thermodynamic folding transition. Folding kinetics is described by a master equation. The folding time is estimated by means of a local thermodynamic equilibrium assumption and then is calculated both numerically and analytically by solving the master equation. The folding time has a maximum near the folding transition temperature and can have a minimum at a lower temperature.

Prediction of protein folding class using global description of amino acid sequence.

We present a method for predicting protein folding class based on global protein chain description and a voting process. Selection of the best descriptors was achieved by a computer-simulated neural network trained on a data base consisting of 83 folding classes. Protein-chain descriptors include overall composition, transition, and distribution of amino acid attributes, such as relative hydrophobicity, predicted secondary structure, and predicted solvent exposure. Cross-validation testing was performed on 15 of the largest classes. The test shows that proteins were assigned to the correct class (correct positive prediction) with an average accuracy of 71.7%, whereas the inverse prediction of proteins as not belonging to a particular class (correct negative prediction) was 90-95% accurate. When tested on 254 structures used in this study, the top two predictions contained the correct class in 91% of the cases.

Negative activation enthalpies in the kinetics of protein folding.

Although the rates of chemical reactions become faster with increasing temperature, the converse may be observed with protein-folding reactions. The rate constant for folding initially increases with temperature, goes through a maximum, and then decreases. The activation enthalpy is thus highly temperature dependent because of a large change in specific heat (delta Cp). Such a delta Cp term is usually presumed to be a consequence of a large decrease in exposure of hydrophobic surfaces to water as the reaction proceeds from the denatured state to the transition state for folding: the hydrophobic side chains are surrounded by "icebergs" of water that melt with increasing temperature, thus making a large contribution to the Cp of the denatured state and a smaller one to the more compact transition state. The rate could also be affected by temperature-induced changes in the conformational population of the ground state: the heat required for the progressive melting of residual structure in the denatured state will contribute to delta Cp. By examining two proteins with different refolding mechanisms, we are able to find both of these two processes; barley chymotrypsin inhibitor 2, which refolds from a highly unfolded state, fits well to a hydrophobic interaction model with a constant delta Cp of activation, whereas barnase, which refolds from a more structured denatured state, deviates from this ideal behavior.

Cyclophilin catalyzes protein folding in yeast mitochondria.

Cyclophilins are a family of ubiquitous proteins that are the intracellular target of the immunosuppressant drug cyclosporin A. Although cyclophilins catalyze peptidylprolyl cis-trans isomerization in vitro, it has remained open whether they also perform this function in vivo. Here we show that Cpr3p, a cyclophilin in the matrix of yeast mitochondria, accelerates the refolding of a fusion protein that was synthesized in a reticulocyte lysate and imported into the matrix of isolated yeast mitochondria. The fusion protein consisted of the matrix-targeting sequence of subunit 9 of F1F0-ATPase fused to mouse dihydrofolate reductase. Refolding of the dihydrofolate reductase moiety in the matrix was monitored by acquisition of resistance to proteinase K. The rate of refolding was reduced by a factor of 2-6 by 2.5 microM cyclosporin A. This reduced rate of folding was also observed with mitochondria lacking Cpr3p. In these mitochondria, protein folding was insensitive to cyclosporin A. The rate of protein import was not affected by cyclosporin A or by deletion of Cpr3p.

The folding of GroEL-bound barnase as a model for chaperonin-mediated protein folding.

We have analyzed the pathway of folding of barnase bound to GroEL to resolve the controversy of whether proteins can fold while bound to chaperonins (GroEL or Cpn60) or fold only after their release into solution. Four phases in the folding were detected by rapid-reaction kinetic measurements of the intrinsic fluorescence of both wild type and barnase mutants. The phases were assigned from their rate laws, sensitivity to mutations, and correspondence to regain of catalytic activity. At high ratios of denatured barnase to GroEL, 4 mol of barnase rapidly bind per 14-mer of GroEL. At high ratios of GroEL to barnase, 1 mol of barnase binds with a rate constant of 3.5 x 10(7) s-1.M-1. This molecule then refolds with a low rate constant that changes on mutation in parallel with the rate constant for the folding in solution. This rate constant corresponds to the regain of the overall catalytic activity of barnase and increases 15-fold on the addition of ATP to a physiologically relevant value of approximately 0.4 s-1. The multiply bound molecules of barnase that are present at high ratios of GroEL to barnase fold with a rate constant that is also sensitive to mutation but is 10 times higher. If the 110-residue barnase can fold when bound to GroEL and many moles can bind simultaneously, then smaller parts of large proteins should be able to fold while bound.

Toward an outline of the topography of a realistic protein-folding funnel.

Experimental information on the structure and dynamics of molten globules gives estimates for the energy landscape's characteristics for folding highly helical proteins, when supplemented by a theory of the helix-coil transition in collapsed heteropolymers. A law of corresponding states relating simulations on small lattice models to real proteins possessing many more degrees of freedom results. This correspondence reveals parallels between "minimalist" lattice results and recent experimental results for the degree of native character of the folding transition state and molten globule and also pinpoints the needs of further experiments.

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.

Lattice model for rapidly folding protein-like heteropolymers.

Protein folding is a relatively fast process considering the astronomical number of conformations in which a protein could find itself. Within the framework of a lattice model, we show that one can design rapidly folding sequences by assigning the strongest attractive couplings to the contacts present in a target native state. Our protein design can be extended to situations with both attractive and repulsive contacts. Frustration is minimized by ensuring that all the native contacts are again strongly attractive. Strikingly, this ensures the inevitability of folding and accelerates the folding process by an order of magnitude. The evolutionary implications of our findings are discussed.

Regulation of the dolichol pathway in human fibroblasts by the endoplasmic reticulum unfolded protein response

Accumulation of unfolded proteins within the endoplasmic reticulum (ER) of eukaryotic cells triggers the unfolded protein response (UPR), which activates transcription of several genes encoding ER chaperones and folding enzymes. This study reports that conversion of dolichol-linked Man2–5GlcNAc2 intermediates into mature Glc3Man9GlcNAc2 oligosaccharides in primary human adult dermal fibroblasts is also stimulated by the UPR. This stimulation was not evident in several immortal cell lines and did not require a cytoplasmic stress response. Inhibition of dolichol-linked Glc3Man9GlcNAc2 synthesis by glucose deprivation could be counteracted by the UPR, improving the transfer of Glc3Man9GlcNAc2 to asparagine residues on nascent polypeptides. Glycosidic processing of asparagine-linked Glc3Man9GlcNAc2 in the ER leads to the production of monoglucosylated oligosaccharides that promote interaction with the lectin chaperones calreticulin and calnexin. Thus, control of the dolichol-linked Glc3Man9GlcNAc2 supply gives the UPR the potential to maintain efficient protein folding in the ER without new synthesis of chaperones or folding enzymes.

Two-dimensional IR correlation spectroscopy: Sequential events in the unfolding process of the λ Cro-V55C repressor protein

A question often posed in protein folding/unfolding studies is whether the process is fully cooperative or whether it contains sequential elements. To address this question, one needs tools capable of resolving different events. It seems that, at least in certain cases, two-dimensional (2D) IR correlation spectroscopy can provide answers to this question. To illustrate this point, we have turned to the Cro-V55C dimer of the λ Cro repressor, a protein known to undergo thermal unfolding in two discrete steps through a stable equilibrium intermediate. The secondary structure of this intermediate is compatible with that of a partially unfolded protein and involves a reorganization of the N terminus, whereas the antiparallel β-ribbon formed by the C-terminal part of each subunit remains largely intact. To establish whether the unfolding process involves sequential events, we have performed a 2D correlation analysis of IR spectra recorded over the temperature range of 20–95°C. The 2D IR correlation analysis indeed provides evidence for a sequential formation of the stable intermediate, which is created in three (closely related) steps. A first step entails the unfolding of the short N-terminal β-strand, followed by the unfolding of the α-helices in a second step, and the third step comprises the reorganization of the remaining β-sheet and of some unordered segments in the protein. The complete unfolding of the stable intermediate at higher temperatures also undergoes sequential events that ultimately end with the breaking of the H bonds between the two β-strands at the dimer interface.

An autonomous folding unit mediates the assembly of two-stranded coiled coils

Subunit oligomerization of many proteins is mediated by coiled-coil domains. Although the basic features contributing to the thermodynamic stability of coiled coils are well understood, the mechanistic details of their assembly have not yet been dissected. Here we report a 13-residue sequence pattern that occurs with limited sequence variations in many two-stranded coiled coils and that is absolutely required for the assembly of the Dictyostelium discoideum actin-bundling protein cortexillin I and the yeast transcriptional activator GCN4. The functional relationship between coiled-coil “trigger” sequences was manifested by replacing the intrinsic trigger motif of GCN4 with the related sequence from cortexillin I. We demonstrate that these trigger sequences represent autonomous helical folding units that, in contrast to arbitrarily chosen heptad repeats, can mediate coiled-coil formation. Aside from being of general interest for protein folding, trigger motifs should be of particular importance in the protein de novo design.

Mammalian Transcription Factor ATF6 Is Synthesized as a Transmembrane Protein and Activated by Proteolysis in Response to Endoplasmic Reticulum Stress

The unfolded protein response (UPR) controls the levels of molecular chaperones and enzymes involved in protein folding in the endoplasmic reticulum (ER). We recently isolated ATF6 as a candidate for mammalian UPR-specific transcription factor. We report here that ATF6 constitutively expressed as a 90-kDa protein (p90ATF6) is directly converted to a 50-kDa protein (p50ATF6) in ER-stressed cells. Furthermore, we showed that the most important consequence of this conversion was altered subcellular localization; p90ATF6 is embedded in the ER, whereas p50ATF6 is a nuclear protein. p90ATF6 is a type II transmembrane glycoprotein with a hydrophobic stretch in the middle of the molecule. Thus, the N-terminal half containing a basic leucine zipper motif is oriented facing the cytoplasm. Full-length ATF6 as well as its C-terminal deletion mutant carrying the transmembrane domain is localized in the ER when transfected. In contrast, mutant ATF6 representing the cytoplasmic region translocates into the nucleus and activates transcription of the endogenous GRP78/BiP gene. We propose that ER stress-induced proteolysis of membrane-bound p90ATF6 releases soluble p50ATF6, leading to induced transcription in the nucleus. Unlike yeast UPR, mammalian UPR appears to use a system similar to that reported for cholesterol homeostasis.

A specific transition state for S-peptide combining with folded S-protein and then refolding

We measured the folding and unfolding kinetics of mutants for a simple protein folding reaction to characterize the structure of the transition state. Fluorescently labeled S-peptide analogues combine with S-protein to form ribonuclease S analogues: initially, S-peptide is disordered whereas S-protein is folded. The fluorescent probe provides a convenient spectroscopic probe for the reaction. The association rate constant, kon, and the dissociation rate constant, koff, were both determined for two sets of mutants. The dissociation rate constant is measured by adding an excess of unlabeled S-peptide analogue to a labeled complex (RNaseS*). This strategy allows kon and koff to be measured under identical conditions so that microscopic reversibility applies and the transition state is the same for unfolding and refolding. The first set of mutants tests the role of the α-helix in the transition state. Solvent-exposed residues Ala-6 and Gln-11 in the α-helix of native RNaseS were replaced by the helix destabilizing residues glycine or proline. A plot of log kon vs. log Kd for this series of mutants is linear over a very wide range, with a slope of −0.3, indicating that almost all of the molecules fold via a transition state involving the helix. A second set of mutants tests the role of side chains in the transition state. Three side chains were investigated: Phe-8, His-12, and Met-13, which are known to be important for binding S-peptide to S-protein and which also contribute strongly to the stability of RNaseS*. Only the side chain of Phe-8 contributes significantly, however, to the stability of the transition state. The results provide a remarkably clear description of a folding transition state.

Folding simulations of a three-stranded antiparallel β-sheet peptide

Protein folding is a grand challenge of the postgenomic era. In this paper, 58 folding events sampled during 47 molecular dynamics trajectories for a total simulation time of more than 4 μs provide an atomic detail picture of the folding of a 20-residue synthetic peptide with a stable three-stranded antiparallel β-sheet fold. The simulations successfully reproduce the NMR solution conformation, irrespective of the starting structure. The sampling of the conformational space is sufficient to determine the free energy surface and localize the minima and transition states. The statistically predominant folding pathway involves the formation of contacts between strands 2 and 3, starting with the side chains close to the turn, followed by association of the N-terminal strand onto the preformed 2–3 β-hairpin. The folding mechanism presented here, formation of a β-hairpin followed by consolidation, is in agreement with a computational study of the free energy surface of another synthetic three-stranded antiparallel β-sheet by Bursulaya and Brooks [(1999) J. Am. Chem. Soc. 121, 9947–9951]. Hence, it might hold in general for antiparallel β-sheets with short turns.

Absence of stable intermediates on the folding pathway of barnase

Barnase is one of the few protein models that has been studied extensively for protein folding. Previous studies led to the conclusion that barnase folds through a very stable submillisecond intermediate (≈3 kcal/mol). The structure of this intermediate was characterized intensively by using a protein engineering approach. This intermediate has now been reexamined with three direct and independent methods. (i) Hydrogen exchange experiments show very small protection factors (≈2) for the putative intermediate, indicating a stability of ≈0.0 kcal/mol. (ii) Denaturant-dependent unfolding of the putative intermediate is noncooperative and indicates a stability less than 0.0 kcal/mol. (iii) The logarithm of the unfolding rate constant of native barnase vs. denaturant concentrations is not linear. Together with the measured rate (“I” to N), this nonlinear behavior accounts for almost all of the protein stability, leaving only about 0.3 kcal/mol that could be attributed to the rapidly formed intermediate. Other observations previously interpreted to support the presence of an intermediate are now known to have alternative explanations. These results cast doubts on the previous conclusions on the nature of the early folding state in barnase and therefore should have important implications in understanding the early folding events of barnase and other proteins in general.