Toward a mechanism for GroEL.GroES chaperone activity: an ATPase-gated and -pulsed folding and annealing cage.

Free GroEL binds denatured proteins very tightly: it retards the folding of barnase 400-fold and catalyzes unfolding fluctuations in native barnase and its folding intermediate. GroEL undergoes an allosteric transition from its tight-binding T-state to a weaker binding R-state on the cooperative binding of nucleotides (ATP/ADP) and GroES. The preformed GroEL.GroES.nucleotide complex retards the folding of barnase by only a factor of 4, and the folding rate is much higher than the ATPase activity that releases GroES from the complex. Binding of GroES and nucleotides to a preformed GroEL.denatured-barnase complex forms an intermediately fast-folding complex. We propose the following mechanism for the molecular chaperone. Denatured proteins bind to the resting GroEL.GroES.nucleotide complex. Fast-folding proteins are ejected as native structures before ATP hydrolysis. Slow-folding proteins enter chaperoning cycles of annealing and folding after the initial ATP hydrolysis. This step causes transient release of GroES and formation of the GroEL.denatured-protein complexes with higher annealing potential. The intermediately fast-folding complex is formed on subsequent rebinding of GroES. The ATPase activity of GroEL.GroES is thus the gatekeeper that selects for initial entry of slow-folding proteins to the chaperone action and then pumps successive transitions from the faster-folding R-states to the tighter-binding/stronger annealing T-states. The molecular chaperone acts as a combination of folding cage and an annealing machine.

Barriers to protein folding: formation of buried polar interactions is a slow step in acquisition of structure.

In the MYL mutant of the Arc repressor dimer, sets of partially buried salt-bridge and hydrogen-bond interactions mediated by Arg-31, Glu-36, and Arg-40 in each subunit are replaced by hydrophobic interactions between Met-31, Tyr-36, and Leu-40. The MYL refolding/dimerization reaction differs from that of wild type in being 10- to 1250-fold faster, having an earlier transition state, and depending upon viscosity but not ionic strength. Formation of the wild-type salt bridges in a hydrophobic environment clearly imposes a kinetic barrier to folding, which can be lowered by high salt concentrations. The changes in the position of the transition state and viscosity dependence can be explained if denatured monomers interact to form a partially folded dimeric intermediate, which then continues folding to form the native dimer. The second step is postulated to be rate limiting for wild type. Replacing the salt bridge with hydrophobic interactions lowers this barrier for MYL. This makes the first kinetic barrier rate limiting for MYL refolding and creates a downhill free-energy landscape in which most molecules which reach the intermediate state continue to form native dimers.

Conformations and folding of lysozyme ions in vacuo.

Proton transfer reactivity of isolated charge states of the protein hen egg-white lysozyme shows that multiple distinct conformations of this protein are stable in the gas phase. The reactivities of the 9+ and 10+ charge state ions, formed by electrospray ionization of "native" (disulfide-intact) and "denatured" (disulfide-reduced) solutions, are consistent with values calculated for ions in their crystal structure and fully denatured conformations, respectively. Charge states below 8+ of both forms, formed by proton stripping, have similar or indistinguishable reactivities, indicating that the disulfide-reduced ions fold in the gas phase to a more compact conformation.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

The protein-folding activity of chaperonins correlates with the symmetric GroEL14(GroES7)2 heterooligomer.

Chaperonins GroEL and GroES form, in the presence of ATP, two types of heterooligomers in solution: an asymmetric GroEL14GroES7 "bullet"-shaped particle and a symmetric GroEL14(GroES7)2 "football"-shaped particle. Under limiting concentrations of ATP or GroES, excess ADP, or in the presence of 5'-adenylyl imidodiphosphate, a correlation is seen between protein folding and the amount of symmetric GroEL14(GroES7)2 particles in a chaperonin solution, as detected by electron microscopy or by chemical crosslinking. Kinetic analysis suggests that protein folding is more efficient when carried out by a chaperonin solution populated with a majority of symmetric GroEL14(GroES7)2 particles than by a majority of asymmetric GroEL14GroES7 particles. The symmetric heterooligomer behaves as a highly efficient intermediate of the chaperonin protein folding cycle in vitro.

Optimization of rates of protein folding: the nucleation-condensation mechanism and its implications.

Small, single-module proteins that fold in a single cooperative step may be paradigms for understanding early events in protein-folding pathways generally. Recent experimental studies of the 64-residue chymotrypsin inhibitor 2 (CI2) support a nucleation mechanism for folding, as do some computer stimulations. CI2 has a nucleation site that develops only in the transition state for folding. The nucleus is composed of a set of adjacent residues (an alpha-helix), stabilized by long-range interactions that are formed as the rest of the protein collapses around it. A simple analysis of the optimization of the rate of protein folding predicts that rates are highest when the denatured state has little residual structure under physiological conditions and no intermediates accumulate. This implies that any potential nucleation site that is composed mainly of adjacent residues should be just weakly populated in the denatured state and become structured only in a high-energy intermediate or transition state when it is stabilized by interactions elsewhere in the protein. Hierarchical mechanisms of folding in which stable elements of structure accrete are unfavorable. The nucleation-condensation mechanism of CI2 fulfills the criteria for fast folding. On the other hand, stable intermediates do form in the folding of more complex proteins, and this may be an unavoidable consequence of increasing size and nucleation at more than one site.

Submillisecond events in protein folding.

The pathway of protein folding is now being analyzed at the resolution of individual residues by kinetic measurements on suitably engineered mutants. The kinetic methods generally employed for studying folding are typically limited to the time range of > or = 1 ms because the folding of denatured proteins is usually initiated by mixing them with buffers that favor folding, and the dead time of rapid mixing experiments is about a millisecond. We now show that the study of protein folding may be extended to the microsecond time region by using temperature-jump measurements on the cold-unfolded state of a suitable protein. We are able to detect early events in the folding of mutants of barstar, the polypeptide inhibitor of barnase. A preliminary characterization of the fast phase from spectroscopic and phi-value analysis indicates that it is a transition between two relatively solvent-exposed states with little consolidation of structure.

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.

Kinetic traps in lysozyme folding.

Folding of lysozyme from hen egg white was investigated by using interrupted refolding experiments. This method makes use of a high energy barrier between the native state and transient folding intermediates, and, in contrast to conventional optical techniques, it enables one to specifically monitor the amount of native molecules during protein folding. The results show that under strongly native conditions lysozyme can refold on parallel pathways. The major part of the lysozyme molecules (86%) refold on a slow kinetic pathway with well-populated partially folded states. Additionally, 14% of the molecules fold faster. The rate constant of formation of native molecules on the fast pathway corresponds well to the rate constant expected for folding to occur by a two-state process without any detectable intermediates. The results suggest that formation of the native state for the major fraction of lysozyme molecules is retarded compared with the direct folding process. Partially structured intermediates that transiently populate seem to be kinetically trapped in a conformation that can only slowly reach the native structure.

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.

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.

On constructing folding heteropolymers.

Simplified models of the protein-folding process have led to valuable insights into the generic properties of the folding of heteropolymers. On the basis of theoretical arguments, Shakhnovich and Gutin [(1993) Proc. Natl. Acad. Sci. USA 90, 7195-7199] have proposed a specific method to generate folding sequences for one of these. Here we present a model of folding in heteropolymers that is comparable in simplicity but different in spirit to the one studied by Shakhnovich and Gutin. In our model, the proposed recipe for constructing folding sequence fails. We find that, as a rule, the construction of folding sequences is impossible to achieve by looking at the native conformation only. Rather, competing conformations have to be taken into account too. An evolutionary algorithm that generates folding sequences by optimizing both stability of the native state and folding time is described. Remarkably, this algorithm produces, among others, sequences that fold reproducibly to metastable states.

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.

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.