Structure and function in rhodopsin: correct folding and misfolding in point mutants at and in proximity to the site of the retinitis pigmentosa mutation Leu-125-->Arg in the transmembrane helix C.

L125R is a mutation in the transmembrane helix C of rhodopsin that is associated with autosomal dominant retinitis pigmentosa. To probe the orientation of the helix and its packing in the transmembrane domain, we have prepared and studied the mutations E122R, I123R, A124R, S127R, L125F, and L125A at, and in proximity to, the above mutation site. Like L125R, the opsin expressed in COS-1 cells from E122R did not bind 11-cis-retinal, whereas those from I123R and S127R formed the rhodopsin chromophore partially. A124R opsin formed the rhodopsin chromophore (lambda max 495 nm) in the dark, but the metarhodopsin II formed on illumination decayed about 6.5 times faster than that of the wild type and was defective in transducin activation. The mutant opsins from L125F and L125A bound 11-cis-retinal only partially, and in both cases, the mixtures of the proteins produced were separated into retinal-binding and non-retinal-binding (misfolded) fractions. The purified mutant rhodopsin from L125F showed lambda max at 500 nm, whereas that from L125A showed lambda max at 503 nm. The mutant rhodopsin L125F showed abnormal bleaching behavior and both mutants on illumination showed destabilized metarhodopsin II species and reduced transducin activation. Because previous results have indicated that misfolding in rhodopsin is due to the formation of a disulfide bond other than the normal disulfide bond between Cys-110 and Cys-187 in the intradiscal domain, we conclude from the misfolding in mutants L125F and L125A that the folding in vivo in the transmembrane domain is coupled to that in the intradiscal domain.

How optimization of potential functions affects protein folding.

The relationship between the optimization of the potential function and the foldability of theoretical protein models is studied based on investigations of a 27-mer cubic-lattice protein model and a more realistic lattice model for the protein crambin. In both the simple and the more complicated systems, optimization of the energy parameters achieves significant improvements in the statistical-mechanical characteristics of the systems and leads to foldable protein models in simulation experiments. The foldability of the protein models is characterized by their statistical-mechanical properties--e.g., by the density of states and by Monte Carlo folding simulations of the models. With optimized energy parameters, a high level of consistency exists among different interactions in the native structures of the protein models, as revealed by a correlation function between the optimized energy parameters and the native structure of the model proteins. The results of this work are relevant to the design of a general potential function for folding proteins by theoretical simulations.

Chaperonin-facilitated protein folding: optimization of rate and yield by an iterative annealing mechanism.

We develop a heuristic model for chaperonin-facilitated protein folding, the iterative annealing mechanism, based on theoretical descriptions of "rugged" conformational free energy landscapes for protein folding, and on experimental evidence that (i) folding proceeds by a nucleation mechanism whereby correct and incorrect nucleation lead to fast and slow folding kinetics, respectively, and (ii) chaperonins optimize the rate and yield of protein folding by an active ATP-dependent process. The chaperonins GroEL and GroES catalyze the folding of ribulose bisphosphate carboxylase at a rate proportional to the GroEL concentration. Kinetically trapped folding-incompetent conformers of ribulose bisphosphate carboxylase are converted to the native state in a reaction involving multiple rounds of quantized ATP hydrolysis by GroEL. We propose that chaperonins optimize protein folding by an iterative annealing mechanism; they repeatedly bind kinetically trapped conformers, randomly disrupt their structure, and release them in less folded states, allowing substrate proteins multiple opportunities to find pathways leading to the most thermodynamically stable state. By this mechanism, chaperonins greatly expand the range of environmental conditions in which folding to the native state is possible. We suggest that the development of this device for optimizing protein folding was an early and significant evolutionary event.

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.