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.

Interaction of GroEL with a highly structured folding intermediate: iterative binding cycles do not involve unfolding.

The GroE proteins are molecular chaperones involved in protein folding. The general mechanism by which they facilitate folding is still enigmatic. One of the central open questions is the conformation of the GroEL-bound nonnative protein. Several suggestions have been made concerning the folding stage at which a protein can interact with GroEL. Furthermore, the possibility exists that binding of the nonnative protein to GroEL results in its unfolding. We have addressed these issues that are basic for understanding the GroE-mediated folding cycle by using folding intermediates of an Fab antibody fragment as molecular probes to define the binding properties of GroEL. We show that, in addition to binding to an early folding intermediate, GroEL is able to recognize and interact with a late quaternary-structured folding intermediate (Dc) without measurably unfolding it. Thus, the prerequisite for binding is not a certain folding stage of a nonnative protein. In contrast, general surface properties of nonnative proteins seem to be crucial for binding. Furthermore, unfolding of a highly structured intermediate does not necessarily occur upon binding to GroEL. Folding of Dc in the presence of GroEL and ATP involves cycles of binding and release. Because in this system no off-pathway reactions or kinetic traps are involved, a quantitative analysis of the reactivation kinetics observed is possible. Our results indicate that the association reaction of Dc and GroEL in the presence of ATP is rather slow, whereas in the absence of ATP association is several orders of magnitude more efficient. Therefore, it seems that ATP functions by inhibiting reassociation rather than promoting release of the bound substrate.

Submillisecond folding of monomeric lambda repressor.

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.

An optimal Mg2+ concentration for kinetic folding of the Tetrahymena ribozyme

Divalent metal ions, such as Mg2+, are generally required for tertiary structure formation in RNA. Although the role of Mg2+ binding in RNA-folding equilibria has been studied extensively, little is known about the role of Mg2+ in RNA-folding kinetics. In this paper, we explore the effect of Mg2+ on the rate-limiting step in the kinetic folding pathway of the Tetrahymena ribozyme. Analysis of these data reveals the presence of a Mg2+-stabilized kinetic trap that slows folding at higher Mg2+ concentrations. Thus, the Tetrahymena ribozyme folds with an optimal rate at 2 mM Mg2+, just above the concentration required for stable structure formation. These results suggest that thermodynamic and kinetic folding of RNA are cooptimized at a Mg2+ concentration that is sufficient to stabilize the folded form but low enough to avoid kinetic traps and misfolding.

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.

Delayed Activation of the Store-operated Calcium Current Induced by Calreticulin Overexpression in RBL-1 Cells

Calreticulin (CRT) is a high-capacity, low-affinity Ca2+-binding protein located in the lumen of the endoplasmic reticulum (ER) of all eukaryotic cells investigated so far. Its high level of conservation among different species suggests that it serves functions fundamental to cell survival. The role originally proposed for CRT, i.e., the main Ca2+ buffer of the ER, has been obscured or even casted by its implication in processes as diverse as gene expression, protein folding, and cell adhesion. In this work we seek the role of CRT in Ca2+ storing and signaling by evaluating its effects on the kinetics and amplitude of the store-operated Ca2+ current (ICRAC). We show that, in the rat basophilic leukemia cell line RBL-1, overexpression of CRT, but not of its mutant lacking the high-capacity Ca2+-binding domain, markedly retards the ICRAC development, however, only when store depletion is slower than the rate of current activation. On the contrary, when store depletion is rapid and complete, overexpression of CRT has no effect. The present results are compatible with a major Ca2+-buffering role of CRT within the ER but exclude a direct, or indirect, role of this protein on the mechanism of ICRAC activation.

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.

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.

Functional expression of the Wilson disease protein reveals mislocalization and impaired copper-dependent trafficking of the common H1069Q mutation

Wilson disease is an autosomal recessive disorder of hepatic copper metabolism caused by mutations in a gene encoding a copper-transporting P-type ATPase. To elucidate the function of the Wilson protein, wild-type and mutant Wilson cDNAs were expressed in a Menkes copper transporter-deficient mottled fibroblast cell line defective in copper export. Expression of the wild-type cDNA demonstrated trans-Golgi network localization and copper-dependent trafficking of the Wilson protein identical to previous observations for the endogenously expressed protein in hepatocytes. Furthermore, expression of the Wilson cDNA rescued the mottled phenotype as evidenced by a reduction in copper accumulation and restoration of cell viability. In contrast, expression of an H1069Q mutant Wilson cDNA did not rescue the mottled phenotype, and immunofluorescence studies showed that this mutant Wilson protein was localized in the endoplasmic reticulum. Consistent with these findings, pulse–chase analysis demonstrated a 5-fold decrease in the half-life of the H1069Q mutant as compared with the wild-type protein. Maintenance of these transfected cell lines at 28°C resulted in localization of the H1069Q protein in the trans-Golgi network, suggesting that a temperature-sensitive defect in protein folding followed by degradation constitutes the molecular basis of Wilson disease in patients harboring the H1069Q mutation. Taken together, these studies describe a tractable expression system for elucidating the function and localization of the copper-transporting ATPases in mammalian cells and provide compelling evidence that the Wilson protein can functionally substitute for the Menkes protein, supporting the concept that these proteins use common biochemical mechanisms to effect cellular copper homeostasis.

Synechocystis HSP17 is an amphitropic protein that stabilizes heat-stressed membranes and binds denatured proteins for subsequent chaperone-mediated refolding

The small heat shock proteins (sHSPs) are ubiquitous stress proteins proposed to act as molecular chaperones to prevent irreversible protein denaturation. We characterized the chaperone activity of Synechocystis HSP17 and found that it has not only protein-protective activity, but also a previously unrecognized ability to stabilize lipid membranes. Like other sHSPs, recombinant Synechocystis HSP17 formed stable complexes with denatured malate dehydrogenase and served as a reservoir for the unfolded substrate, transferring it to the DnaK/DnaJ/GrpE and GroEL/ES chaperone network for subsequent refolding. Large unilamellar vesicles made of synthetic and cyanobacterial lipids were found to modulate this refolding process. Investigation of HSP17-lipid interactions revealed a preference for the liquid crystalline phase and resulted in an elevated physical order in model lipid membranes. Direct evidence for the participation of HSP17 in the control of thylakoid membrane physical state in vivo was gained by examining an hsp17− deletion mutant compared with the isogenic wild-type hsp17+ revertant Synechocystis cells. We suggest that, together with GroEL, HSP17 behaves as an amphitropic protein and plays a dual role. Depending on its membrane or cytosolic location, it may function as a “membrane stabilizing factor” as well as a member of a multichaperone protein-folding network. Membrane association of sHSPs could antagonize the heat-induced hyperfluidization of specific membrane domains and thereby serve to preserve structural and functional integrity of biomembranes.

Hydrophobic moments of protein structures: Spatially profiling the distribution

It is generally accepted that globular proteins fold with a hydrophobic core and a hydrophilic exterior. Might the spatial distribution of amino acid hydrophobicity exhibit common features? The hydrophobic profile detailing this distribution from the protein interior to exterior has been examined for 30 relatively diverse structures obtained from the Protein Data Bank, for 3 proteins of the 30S ribosomal subunit, and for a simple set of 14 decoys. A second-order hydrophobic moment has provided a simple measure of the spatial variation. Shapes of the calculated spatial profiles of all native structures have been found to be comparable. Consequently, profile shapes as well as particular profile features should assist in validating predicted protein structures and in discriminating between different protein-folding pathways. The spatial profiles of the 14 decoys are clearly distinguished from the profiles of their native structures.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.