Protein folding in the cell: molecular chaperones pave the way

In vitro, many unfolded polypeptides are able to fold to the native state spontaneously, indicating that the amino acid sequence of a protein contains all the information necessary to specify its three-dimensional conformation. It had been assumed that protein folding in vivo also generally occurs in a spontaneous process. This view has changed only recently due to the discovery of a number of proteins, now commonly called 'molecular chaperones', which are essential for cellular protein folding and occur ubiquitously in eubacteria, archaebacteria and in eukaryotic cells.     

A wheat germ cell-free system is a novel way to screen protein folding and function

For high-throughput protein structural analysis, it is indispensable to develop a reliable protein overexpression system. Although many protein overexpression systems, such as that involving Escherichia coli cells, have been developed, the number of overexpressed proteins showing the same biological activities as those of the native proteins is limited. A novel wheat germ cell-free protein synthesis system was developed recently, and most of the proteins functioning in solution were synthesized as soluble forms. This suggests the applicability of this protein synthesis method to determination of the solution structures of functional proteins. To examine this possibility, we have synthesized two (15)N-labeled proteins and obtained (1)H-(15)N HSQC spectra for them. The structural analysis of these proteins has already progressed with an E. coli overexpression system, and (1)H-(15)N HSQC spectra for biologically active proteins have already been obtained. Comparing the spectra, we have shown that proteins synthesized with a wheat germ cell-free system have the proper protein folding and enough biological activity. This is the first experimental evidence of the applicability of the wheat germ cell-free protein synthesis system to high-throughput protein structural analysis.     

A simple way to measure protein refolding rates in water.

Refolding of proteins is traditionally carried out either by diluting the denaturant-unfolded protein into buffer (GdmCl-jump) or by mixing the acid-denatured protein with strong buffer (pH-jump). The first method does not allow direct measurement of folding rates in water since the GdmCl cannot be infinitely diluted, and the second method suffers from the limitation that many proteins cannot be pH-denatured. Further, some proteins do not refold reversibly from low pH where they get trapped as aggregation prone intermediates. Here, we present an alternative approach for direct measurement of refolding rates in water, which does not rely on extrapolation. The protein is denatured in SDS, and is then mixed with alpha-cyclodextrin, which rapidly strips SDS molecules from the protein, leaving the naked unfolded protein to refold.     

Diffusional barrier crossing in a two-state protein folding reaction.

There has been some debate as to whether protein folding involves diffusive chain motions and thus depends on solvent viscosity. The interpretation of folding kinetics in viscous solvents has remained difficult and controversial, in that viscogenic agents affect folding rates not only by increasing solvent viscosity but also by increasing protein stability. By carefully choosing experimental conditions, we can now eliminate the effect on stability and show that the folding dynamics of the cold shock protein CspB are viscosity dependent. Thus Kramers' theory of reaction rates rather than transition state theory should be used to describe this folding reaction.     

Increasing stability reduces conformational heterogeneity in a protein folding intermediate ensemble

A multi-site, time-resolved fluorescence resonance energy transfer methodology has been used to study structural heterogeneity in a late folding intermediate ensemble, IL, of the small protein barstar. Four different intra-molecular distances have been measured within the structural components of IL. The IL ensemble is shown to consist of different sub-populations of molecules, in each of which one or more of the four distances are native-like and the remaining distances are unfolded-like. In very stable conditions that favor formation of IL, all four distances are native-like in most molecules. In less stable conditions, one or more distances are unfolded-like. As stability is decreased, the proportion of molecules with unfolded-like distances increases. Thus, the results show that protein folding intermediates are ensembles of different structural forms, and they demonstrate that conformational entropy increases as structures become less stable. These observations provide direct experimental evidence in support of a basic tenet of energy landscape theory for protein folding, that available conformational space, as represented by structural heterogeneity in IL, becomes restricted as the stability is increased. The results also vindicate an important prediction of energy landscape theory, that different folding pathways may become dominant under different folding conditions. In more stable folding conditions, uniformly native-like compactness is achieved during folding to IL, whereas in less stable conditions, uniformly native-like compactness is achieved only later during the folding of IL to N.     

Structural basis for altering the stability of homologous RNAs from a mesophilic and a thermophilic bacterium

Tertiary RNA structures from thermophilic bacteria generally are more stable than their mesophilic homologs. To understand the structural basis of the increase in stability, we investigated equilibrium folding of the specificity domain (S-domain) of RNase P RNA from a mesophilic (Escherichia coli) and a thermophilic (Thermus thermophilus) bacterium. Equilibrium folding of both S-domains is described by a minimal, three-state folding scheme, U-to-I-to-N. In the I-to-N transition of the thermophilic S-domain, more structure forms and protections are stronger against T1 nuclease and hydroxyl radical reactions. Phylogenetic comparison in the context of the native structure reveals that among 39 nucleotide differences between these S-domains, 12 likely contribute to higher stability. These residues participate in extensive networks of hydrogen bonding, stacking, and metal ion coordination throughout the molecule. The thermophilic S-domain achieves higher stability by mutating strategic base pairs to G-C, decreasing surface accessibility of the native state, and increasing the amount of structure formation in the native folding transition. An E. coli S-domain mutant containing these 12 nt has the same stability and folding cooperativity as the T. thermophilus S-domain. E. coli S-domain mutants containing a subset of 4 or 6 nt have the same stability as the T. thermophilus S-domain but the same folding cooperativity as the E. coli S-domain. These results show that increasing stability can be accomplished by mutations within a local structure, but increasing folding cooperativity needs concerted changes among multiple structural units.     

Practicality and time complexity of a sparsified RNA folding algorithm

Commonly used RNA folding programs compute the minimum free energy structure of a sequence under the pseudoknot exclusion constraint. They are based on Zuker's algorithm which runs in time O(n(3)). Recently, it has been claimed that RNA folding can be achieved in average time O(n(2)) using a sparsification technique. A proof of quadratic time complexity was based on the assumption that computational RNA folding obeys the "polymer-zeta property". Several variants of sparse RNA folding algorithms were later developed. Here, we present our own version, which is readily applicable to existing RNA folding programs, as it is extremely simple and does not require any new data structure. We applied it to the widely used Vienna RNAfold program, to create sibRNAfold, the first public sparsified version of a standard RNA folding program. To gain a better understanding of the time complexity of sparsified RNA folding in general, we carried out a thorough run time analysis with synthetic random sequences, both in the context of energy minimization and base pairing maximization. Contrary to previous claims, the asymptotic time complexity of a sparsified RNA folding algorithm using standard energy parameters remains O(n(3)) under a wide variety of conditions. Consistent with our run-time analysis, we found that RNA folding does not obey the "polymer-zeta property" as claimed previously. Yet, a basic version of a sparsified RNA folding algorithm provides 15- to 50-fold speed gain. Surprisingly, the same sparsification technique has a different effect when applied to base pairing optimization. There, its asymptotic running time complexity appears to be either quadratic or cubic depending on the base composition. The code used in this work is available at: .     

Growing up in a dangerous environment: a network of multiple targeting and folding pathways for nascent polypeptides in the cytosol

The first events in the lives of proteins are the most hazardous. Starting at the ribosome, nascent polypeptides undergo complex folding processes endangered by aggregation reactions. Proteins with organellar destinations require correct targeting to the translocation machineries and prevention from premature folding. The high precision and speed of these processes is ensured by a cystosolic system consisting of molecular chaperones, folding catalysts and targeting factors. This review focuses on the interactions of this system with nascent polypeptides and discusses new concepts for protein folding in the cytosol. It is proposed that folding and targeting are promoted by a flexible network of multiple unassisted and assisted pathways.     

An on-pathway intermediate in the folding of a PDZ domain

The folding pathways of some proteins include the population of partially structured species en route to the native state. Identification and characterization of these folding intermediates are particularly difficult as they are often only transiently populated and play different mechanistic roles, being either on-pathway productive species or off-pathway kinetic traps. To define the role of folding intermediates, a quantitative analysis of the folding and unfolding rate constants over a wide range of denaturant concentration is often required. Such a task is further complicated by the reversible nature of the folding reaction, which implies the observed kinetics to be governed by a complex combination of different microscopic rate constants. Here, we tackled this problem by measuring directly the folding rate constant under highly denaturing conditions, namely by inducing the folding of a PDZ domain through a quasi-irreversible binding reaction with a specific peptide. In analogy with previous works based on hydrogen exchange experiments, we present evidence that the folding pathway of the PDZ domain involves the formation of an obligatory on-pathway intermediate. The results presented exemplify a novel type of kinetic test to detect an on-pathway folding intermediate.     

Squeezed exponential kinetics to describe a nonglassy downhill folding as observed in a lattice protein model

We previously studied the so-called strange kinetics in the two-dimensional lattice HP model. To further study the strange kinetics, folding processes of a 27-mer cubic lattice protein model with Gō potential were investigated by simulating how the bundle of folding trajectories, consisting of a number of independent Monte Carlo simulations, evolves as the folding reaction proceeds, covering a wide range of temperature. Three realms of folding kinetics were observed depending on temperature. Although at temperatures where folding was two-state-like, the kinetics was conventional single exponential, we found that the time course data were well represented by a squeezed (or "shrunken") exponential function, exp [-(t/tau)beta] with beta > 1, at temperatures lower than the folding temperature, where folding was fastest and of a nonglassy downhill type. The squeezed exponential kinetics was found to pertain to the subdiffusion on the nonglassy downhill free energy surface and presents a marked contrast both to the single exponential kinetics and to the stretched exponential kinetics that was observed at lower temperatures where folding was also downhill but topological frustration came into effect. The observed temperature dependence of the folding kinetics suggests that some small single-domain proteins may follow the squeezed exponential kinetics at about the room temperature.     

Evidence for identity between the equilibrium unfolding intermediate and a transient folding intermediate: a comparative study of the folding reactions of alpha-lactalbumin and lysozyme

The refolding kinetics of alpha-lactalbumin at different concentrations of guanidine hydrochloride have been investigated by means of kinetic circular dichroism and stopped-flow absorption measurements. The refolding reaction consists of at least two stages, the instantaneous accumulation of the transient intermediate that has peptide secondary structure and the subsequent slow process associated with formation of tertiary structure. The transient intermediate is compared with the well-characterized equilibrium intermediate observed during the denaturant-induced unfolding. Stabilities of the secondary structures against the denaturant, affinities for Ca2+, and tryptophan absorption properties of the transient and equilibrium intermediates were investigated. In all of these respects, the transient intermediate is identical with the equilibrium one, demonstrating the validity of the use of the equilibrium intermediate as a model of the folding intermediate. Essentially the same transient intermediate was also detected in the folding of lysozyme, the protein known to be homologous to alpha-lactalbumin but whose equilibrium unfolding is represented as a two-state reaction. The stability and cooperativity of the secondary structure of the intermediate of lysozyme are compared with those of alpha-lactalbumin. The results show that the protein folding occurring via the intermediate is not limited to the proteins that show equilibrium intermediates. Although the unfolding equilibria of most proteins are well approximated as a two-state reaction, the two-state hypothesis may not be applicable to the folding reaction under the native condition. Two models of protein folding, intermediate-controlled folding model and multiple-pathway folding model, which are different in view of the role of the intermediate in determining the pathway of folding, are also discussed.     

Desolvation Barrier Effects Are a Likely Contributor to the Remarkable Diversity in the Folding Rates of Small Proteins

The variation in folding rate among single-domain natural proteins is tremendous, but common models with explicit representations of the protein chain are either demonstrably insufficient or unclear as to their capability for rationalizing the experimental diversity in folding rates. In view of the critical role of water exclusion in cooperative folding, we apply native-centric, coarse-grained chain modeling with elementary desolvation barriers to investigate solvation effects on folding rates. For a set of 13 proteins, folding rates simulated with desolvation barriers cover ∼ 4.6 orders of magnitude, spanning a range essentially identical to that observed experimentally. In contrast, folding rates simulated without desolvation barriers cover only ∼ 2.2 orders of magnitude. Following a Hammond-like trend, the folding transition-state ensemble (TSE) of a protein model with desolvation barriers generally has a higher average number of native contacts and is structurally more specific, that is, less diffused, than the TSE of the corresponding model without desolvation barriers. Folding is generally significantly slower in models with desolvation barriers because of their higher overall macroscopic folding barriers as well as slower conformational diffusion speeds in the TSE that are ≈ 1/50 times those in models without desolvation barriers. Nonetheless, the average root-mean-square deviation between the TSE and the native conformation is often similar in the two modeling approaches, a finding suggestive of a more robust structural requirement for the folding rate-limiting step. The increased folding rate diversity in models with desolvation barriers originates from the tendency of these microscopic barriers to cause more heightening of the overall macroscopic folding free-energy barriers for proteins with more nonlocal native contacts than those with fewer such contacts. Thus, the enhancement of folding cooperativity by solvation effects is seen as positively correlated with a protein's native topological complexity.     

Engineering and characterization of a superfolder green fluorescent protein

Existing variants of green fluorescent protein (GFP) often misfold when expressed as fusions with other proteins. We have generated a robustly folded version of GFP, called 'superfolder' GFP, that folds well even when fused to poorly folded polypeptides. Compared to 'folding reporter' GFP, a folding-enhanced GFP containing the 'cycle-3' mutations and the 'enhanced GFP' mutations F64L and S65T, superfolder GFP shows improved tolerance of circular permutation, greater resistance to chemical denaturants and improved folding kinetics. The fluorescence of Escherichia coli cells expressing each of eighteen proteins from Pyrobaculum aerophilum as fusions with superfolder GFP was proportional to total protein expression. In contrast, fluorescence of folding reporter GFP fusion proteins was strongly correlated with the productive folding yield of the passenger protein. X-ray crystallographic structural analyses helped explain the enhanced folding of superfolder GFP relative to folding reporter GFP.     

There is more than one way to turn a spherical cellular monolayer inside out: type B embryo inversion in Volvox globator

Background

Epithelial folding is a common morphogenetic process during the development of multicellular organisms. In metazoans, the biological and biomechanical processes that underlie such three-dimensional (3D) developmental events are usually complex and difficult to investigate. Spheroidal green algae of the genus Volvox are uniquely suited as model systems for studying the basic principles of epithelial folding. Volvox embryos begin life inside out and then must turn their spherical cell monolayer outside in to achieve their adult configuration; this process is called 'inversion.' There are two fundamentally different sequences of inversion processes in Volvocaceae: type A and type B. Type A inversion is well studied, but not much is known about type B inversion. How does the embryo of a typical type B inverter, V. globator, turn itself inside out?

Results

In this study, we investigated the type B inversion of V. globator embryos and focused on the major movement patterns of the cellular monolayer, cell shape changes and changes in the localization of cytoplasmic bridges (CBs) connecting the cells. Isolated intact, sectioned and fragmented embryos were analyzed throughout the inversion process using light microscopy, confocal laser scanning microscopy, scanning electron microscopy and transmission electron microscopy techniques. We generated 3D models of the identified cell shapes, including the localizations of CBs. We show how concerted cell-shape changes and concerted changes in the position of cells relative to the CB system cause cell layer movements and turn the spherical cell monolayer inside out. The type B inversion of V. globator is compared to the type A inversion in V. carteri.

Conclusions

Concerted, spatially and temporally coordinated changes in cellular shapes in conjunction with concerted migration of cells relative to the CB system are the causes of type B inversion in V. globator. Despite significant similarities between type A and type B inverters, differences exist in almost all details of the inversion process, suggesting analogous inversion processes that arose through parallel evolution. Based on our results and due to the cellular biomechanical implications of the involved tensile and compressive forces, we developed a global mechanistic scenario that predicts epithelial folding during embryonic inversion in V. globator.     

Trehalose favors a cutinase compact intermediate off-folding pathway

The folding of cutinase, an enzyme displaying lipolytic activity, has been studied in the presence of trehalose. Equilibrium unfolding data show that trehalose increases the free energy change between folded and unfolded states. Unfolding kinetics reveal the presence of an intermediate which is ca. 60% folded in terms of solvent exposure. Trehalose stabilizes this intermediate relative to the folded state. In contrast, the intermediate revealed by folding kinetics is more compact than the transition state, as shown by the positive slope observed at low denaturant concentration in the chevron plot, as well as the decrease in the observable rate constant for folding with the increase in trehalose concentration. This intermediate displays more than 50% of area buried from the solvent (relative to the native state) compared to around 40% for the transition state for folding and therefore appears to be off the folding pathway. Trehalose stabilizes and guanidine hydrochloride destabilizes this compact intermediate. Both unfolding and folding kinetics show that compact conformational states are stabilized by trehalose, in agreement with current models on the effect of compatible solutes. This effect occurs even for compact states that decelerate the folding as in the case of the intermediate revealed by folding kinetics.     

Observation of Continuous Contraction and a Metastable Misfolded State during the Collapse and Folding of a Small Protein

To obtain proper insight into how structure develops during a protein folding reaction, it is necessary to understand the nature and mechanism of the polypeptide chain collapse reaction, which marks the initiation of folding. Here, the time-resolved fluorescence resonance energy transfer technique, in which the decay of the fluorescence light intensity with time is used to determine the time evolution of the distribution of intra-molecular distances, has been utilized to study the folding of the small protein, monellin. It is seen that when folding begins, about one-third of the protein molecules collapse into a molten globule state (I_MG), from which they relax by continuous further contraction to transit to the native state. The larger fraction gets trapped into a metastable misfolded state. Exit from this metastable state occurs via collapse to the lower free energy I_MG state. This exit is slow, on a time scale of seconds, because of activation energy barriers. The trapped misfolded molecules as well as the I_MG molecules contract continuously and slowly as structure develops. A phenomenological model of Markovian evolution of the polymer chain undergoing folding, incorporating these features, has been developed, which fits well the experimentally observed time evolution of distance distributions. The observation that the “wrong turn” to the misfolded state has not been eliminated by evolution belies the common belief that the folding of functional protein sequences is very different from that of a random heteropolymer, and the former have been selected by evolution to fold quickly.     

Complex folding pathways in a simple beta-hairpin

The determination of the folding mechanisms of proteins is critical to understand the topological change that can propagate Alzheimer and Creutzfeld-Jakobs diseases, among others. The computational community has paid considerable attention to this problem; however, the associated time scale, typically on the order of milliseconds or more, represents a formidable challenge. Ab initio protein folding from long molecular dynamics simulations or ensemble dynamics is not feasible with ordinary computing facilities and new techniques must be introduced. Here we present a detailed study of the folding of a 16-residue beta-hairpin, described by a generic energy model and using the activation-relaxation technique. From a total of 90 trajectories at 300 K, three folding pathways emerge. All involve a simultaneous optimization of the complete hydrophobic and hydrogen bonding interactions. The first two pathways follow closely those observed by previous theoretical studies (folding starting at the turn or by interactions between the termini). The third pathway, never observed by previous all-atom folding, unfolding, and equilibrium simulations, can be described as a reptation move of one strand of the beta-sheet with respect to the other. This reptation move indicates that non-native interactions can play a dominant role in the folding of secondary structures. Furthermore, such a mechanism mediated by non-native hydrogen bonds is not available for study by unfolding and Gō model simulations. The exact folding path followed by a given beta-hairpin is likely to be influenced by its sequence and the solvent conditions. Taken together, these results point to a more complex folding picture than expected for a simple beta-hairpin.     

Using Caenorhabditis elegans as a Model System to Study Protein Homeostasis in a Multicellular Organism

The folding and assembly of proteins is essential for protein function, the long-term health of the cell, and longevity of the organism. Historically, the function and regulation of protein folding was studied in vitro, in isolated tissue culture cells and in unicellular organisms. Recent studies have uncovered links between protein homeostasis (proteostasis), metabolism, development, aging, and temperature-sensing. These findings have led to the development of new tools for monitoring protein folding in the model metazoan organism Caenorhabditis elegans. In our laboratory, we combine behavioral assays, imaging and biochemical approaches using temperature-sensitive or naturally occurring metastable proteins as sensors of the folding environment to monitor protein misfolding. Behavioral assays that are associated with the misfolding of a specific protein provide a simple and powerful readout for protein folding, allowing for the fast screening of genes and conditions that modulate folding. Likewise, such misfolding can be associated with protein mislocalization in the cell. Monitoring protein localization can, therefore, highlight changes in cellular folding capacity occurring in different tissues, at various stages of development and in the face of changing conditions. Finally, using biochemical tools ex vivo, we can directly monitor protein stability and conformation. Thus, by combining behavioral assays, imaging and biochemical techniques, we are able to monitor protein misfolding at the resolution of the organism, the cell, and the protein, respectively.     

Observing folding pathways and kinetics of a single sodium-proton antiporter from Escherichia coli

Mechanisms of folding and misfolding of membrane proteins are of interest in cell biology. Recently, we have established single-molecule force spectroscopy to observe directly the stepwise folding of the Na+/H+ antiporter NhaA from Escherichia coli in vitro. Here, we improved this approach significantly to track the folding intermediates of a single NhaA polypeptide forming structural segments such as the Na+-binding site, transmembrane alpha-helices, and helical pairs. The folding rates of structural segments ranged from 0.31 s(-1) to 47 s(-1), providing detailed insight into a distinct folding hierarchy of an unfolded polypeptide into the native membrane protein structure. In some cases, however, the folding chain formed stable and kinetically trapped non-native structures, which could be assigned to misfolding events of the antiporter.