The folding of an enzyme. VI. The folding pathway of barnase: comparison with theoretical models.

The sequence of events in the refolding pathway of barnase has been analysed to search for general principles in protein folding. There appears to be a correlation between burying hydrophobic surface area and early folding events. All the regions that fold early interact extensively with the beta-sheet. These interactions involve predominantly hydrophobic interactions and the burial of very extensive hydrophobic areas in which multiple, close, hydrophobic-hydrophobic contacts are established around a central group of aliphatic residues. There is no burial of hydrophilic residues in these regions; those that are partly screened from the solvent make hydrogen bonds. All the regions or interactions that are made late in the folding pathway do not make extensive contacts with the beta-sheet. Their buried hydrophobic regions lack a central hydrophobic residue or residues around which other hydrophobic residues pack. Further, in some of these regions there is an extensive burial of hydrophilic residues. The results are consistent with one of the earlier events in protein folding being the local formation of native-like secondary structure elements driven by local hydrophobic surface burial. A possible candidate for an initiation site is a beta-hairpin between beta-strands 3 and 4 that is conserved in the microbial ribonuclease family. A comparison of structures in this family shows that those regions that can be superimposed, or have sequence homology, correspond to elements of structure that are formed and interact with each other early in the folding pathway, suggesting that some of these residues could be involved in directing the folding process. The data on barnase combined with results from other laboratories suggest the following tentative conclusions for the refolding of small monomeric proteins. (1) The refolding pathway is, at least in part, sequential and of compulsory order. (2) Secondary structure formation is driven by local hydrophobic surface burial and precedes the formation of most tertiary interactions. These elements are then stabilized and sometimes elongated by tertiary interactions. It is plausible that there are stop signals encoded in the linear sequence that prevent the elongation of isolated secondary structure elements in solution to a larger extent than is found in the folded protein. (3) Many tertiary interactions are not very constrained in the intermediate but become more and more defined as the hydrophobic cores consolidate, loop structures form and the configuration of surface residues takes place. The interactions between different elements of secondary structure are the last ones to be consolidated while the interactions within the secondary structure elements are consolidated earlier.(ABSTRACT TRUNCATED AT 400 WORDS)     

The folding of an enzyme. V. H/2H exchange-nuclear magnetic resonance studies on the folding pathway of barnase: complementarity to and agreement with protein engineering studies.

Two major methods are currently being used to characterize transient intermediates during protein folding at the level of individual residues. Nuclear magnetic resonance (n.m.r.) measurements on the protection of peptide NH hydrogens against exchange with solvent during refolding can provide information about secondary structure formation. Protein engineering and kinetics can provide direct information about intramolecular interactions of protein side-chains and indirect evidence on secondary structure. These procedures have provided the most complete pictures so far about protein folding intermediates. Both methods have been applied to the characterization of an intermediate in the refolding of barnase. Although the two methods give complementary information, there are some regions of the protein where the methods overlap well. We show that, with one possible exception that is obscure, n.m.r. and protein engineering give identical results for those interactions that can be analysed by both methods. This suggests that these are valid approaches for the study of protein folding intermediates in the case of barnase and that the combination of the methods is a powerful analytical procedure. Information provided by n.m.r. data that is complementary to the protein engineering experiments is: (1) early formation of the C terminus of helix2; (2) early formation of helix3; (3) early formation of several beta-turns (46-49, 101-104 in loop5); and (5) partial formation of loop5. Confirmatory evidence of protein engineering data on the intermediate is: (1) helix1 is complete from residues 10 to 18; (2) the interactions between all beta-strands are present; (3) part of loop2 is not formed; (4) part of loop3 is formed; and (5) some specific tertiary interactions are not made. For some interactions the protein engineering and H/2H exchange methods overlap directly. The information obtained for direct overlap is self consistent.     

The folding of an enzyme. II. Substructure of barnase and the contribution of different interactions to protein stability.

Barnase is described anatomically in terms of its substructures and their mode of packing. The surface area of hydrophobic residues buried on formation and packing of the structural elements has been calculated. Changes in stability have been measured for 64 mutations, 41 constructed in this study, strategically located over the protein. The purpose is to provide: (1) information on the magnitudes of changes in stabilization energy for mutations of residues that are important in maintaining the structure; and (2) probes for the folding pathway to be used in subsequent studies. The majority of mutations delete functional moieties of side-chains or make isosteric changes. The energetics of the interactions are variable and context-dependent. The following general conclusions may be drawn, however, from this study about the classes of interactions that stabilize the protein. (1) Truncation of buried hydrophobic side-chains has, in general, the greatest effect on stability. For fully buried residues, this averages at 1.5 kcal mol-1 per methylene group with a standard deviation of +/- 0.6 kcal mol-1. Truncation of partly exposed leucine, isoleucine or valine residues that are in the range of 50 to 80 A2 of solvent-accessible area (30 to 50% of the total solvent-accessible area on a Gly-X-Gly tripeptide, i.e. those packed against the surface) has a smaller, but relatively constant effect on stability, at 0.81 kcal mol-1 per methylene group with a statistical standard deviation of +/- 0.18 kcal mol-1. (2) There is a very poor correlation between hydrophobic surface area buried and the free energy change for an extensive data set of hydrophobic mutants. The best correlation is found to be between the free energy change and the number of methylene groups within a 6 A radius of the hydrophobic groups deleted. (3) Burial of the hydroxyl group of threonine in a pocket that is intended for a gamma-methyl group of valine costs 2.5 kcal mol-1, in the range expected for the loss of two hydrogen bonds.(ABSTRACT TRUNCATED AT 400 WORDS)     

The analysis of protein folding kinetic data produced in protein engineering experiments

Over the past decade, the "protein engineering method" has been used to investigate the folding pathways of more than 20 different proteins. This method involves measuring the folding and unfolding rates of mutant proteins with single amino acid substitutions spread across the sequence. Comparison of folding rates of the mutant proteins to that of the wild-type protein allows the calculation of the phi value, which can be used to evaluate the stabilizing contribution of an amino acid side chain to the structure of the folding transition state. Here, we review the methodology for analysing data collected in protein engineering folding kinetics studies. We discuss the calculation of folding rates and kinetic m values, the estimation of errors in folding kinetics experiments, phi value calculation including potential pitfalls of the analysis, Br?nsted plots, detecting Hammond behaviour, and the analysis of curved chevron plots.     

Pathway and stability of protein folding.

We describe an experimental approach to the problem of protein folding and stability which measures interaction energies and maps structures of intermediates and transition states during the folding pathway. The strategy is based on two steps. First, protein engineering is used to remove interactions that stabilize defined positions in barnase, the RNAse from Bacillus amyloliquefaciens. The consequent changes in stability are measured from the changes in free energy of unfolding of the protein. Second, each mutation is used as a probe of the structure around the wild-type side chain during the folding process. Kinetic measurements are made on the folding and unfolding of wild-type and mutant proteins. The kinetic and thermodynamic data are combined and analysed to show the role of individual side chains in the stabilization of the folded, transition and intermediate states of the protein. The protein engineering experiments are corroborated by nuclear magnetic resonance studies of hydrogen exchange during the folding process. Folding is a multiphasic process in which alpha-helices and beta-sheet are formed relatively early. Formation of the hydrophobic core by docking helix and sheet is (partly) rate determining. The final steps involve the forming of loops and the capping of the N-termini of helices.     

The folding of an enzyme. III. Structure of the transition state for unfolding of barnase analysed by a protein engineering procedure.

The structure of the first significant transition state on the unfolding pathway of barnase has been analysed in detail by protein engineering methods. Over 50 mutations placed strategically over the whole protein have been used as probes to report on the local structure in the transition state. Several different probes for many regions of the protein give consistent results as do multiple probes at the same site. The overall consistency of phi values indicates that the mutations have not produced changes in the protein that significantly alter the transition state for unfolding. A fine-structure analysis of interactions has also been conducted by removing different parts of the same side-chains. Many of the results of simple mutations fall nicely into the two clear-cut cases of phi = 1 or 0, indicating that the local noncovalent bonds are either fully broken or fully made in the transition state. Much of the structure of barnase in the transition state for unfolding is very similar to that in the folded protein. Both major alpha-helices fray at the N terminus. The last two turns in helix1 are certainly intact, as is the C terminus of helix2. The general picture of the beta-sheet is that the three central beta-strands are completely intact while the two edge beta-strands are mainly present but certainly weakened. The first five residues of the protein unwind but the C terminus remains folded. Three of the five loops are unfolded. The edges of the main hydrophobic core (core1) are significantly weakened, however, and their breaking appears partly rate determining. The centre of the small hydrophobic core3 remains intact. Core2 is completely disrupted. The first events in unfolding are thus: the unfolding of several loops, the unwinding of the helices from the N termini, and the weakening and disruption of the hydrophobic cores. The values of phi are found to be substantially the same under conditions that favour folding as under conditions that are highly denaturing, and so the structure of the unfolding transition state is substantially the same in water as in the presence of denaturant. The structure of the final kinetically significant transition state for refolding is identical to that for unfolding. The final events in refolding are, accordingly, the consolidation of the hydrophobic cores, the closing of many loops and the capping of the N termini of the helices.     

Theory of protein folding

Protein folding should be complex. Proteins organize themselves into specific three-dimensional structures, through a myriad of conformational changes. The classical view of protein folding describes this process as a nearly sequential series of discrete intermediates. In contrast, the energy landscape theory of folding considers folding as the progressive organization of an ensemble of partially folded structures through which the protein passes on its way to the natively folded structure. As a result of evolution, proteins have a rugged funnel-like landscape biased toward the native structure. Connecting theory and simulations of minimalist models with experiments has completely revolutionized our understanding of the underlying mechanisms that control protein folding.     

Computer-based redesign of a protein folding pathway.

A fundamental test of our current understanding of protein folding is to rationally redesign protein folding pathways. We use a computer-based design strategy to switch the folding pathway of protein G, which normally involves formation of the second, but not the first, beta-turn at the rate limiting step in folding. Backbone conformations and amino acid sequences that maximize the interaction density in the first beta-hairpin were identified, and two variants containing 11 amino acid replacements were found to be approximately 4 kcal mol-1 more stable than wild type protein G. Kinetic studies show that the redesigned proteins fold approximately 100 x faster than wild type protein and that the first beta-turn is formed and the second disrupted at the rate limiting step in folding.     

Rational engineering of enzyme stability

During the past 15 years there has been a continuous flow of reports describing proteins stabilized by the introduction of mutations. These reports span a period from pioneering rational design work on small enzymes such as T4 lysozyme and barnase to protein design, and directed evolution. Concomitantly, the purification and characterization of naturally occurring hyperstable proteins has added to our understanding of protein stability. Along the way, many strategies for rational protein stabilization have been proposed, some of which (e.g. entropic stabilization by introduction of prolines or disulfide bridges) have reasonable success rates. On the other hand, comparative studies and efforts in directed evolution have revealed that there are many mutational strategies that lead to high stability, some of which are not easy to define and rationalize. Recent developments in the field include increasing awareness of the importance of the protein surface for stability, as well as the notion that normally a very limited number of mutations can yield a large increase in stability. Another development concerns the notion that there is a fundamental difference between the "laboratory stability" of small pure proteins that unfold reversibly and completely at high temperatures and "industrial stability", which is usually governed by partial unfolding processes followed by some kind of irreversible inactivation process (e.g. aggregation). Provided that one has sufficient knowledge of the mechanism of thermal inactivation, successful and efficient rational stabilization of enzymes can be achieved.     

Potential of genetic algorithms in protein folding and protein engineering simulations.

Genetic algorithms are very efficient search mechanisms which mutate, recombine and select amongst tentative solutions to a problem until a near optimal one is achieved. We introduce them as a new tool to study proteins. The identification and motivation for different fitness functions is discussed. The evolution of the zinc finger sequence motif from a random start is modelled. User specified changes of the lambda repressor structure were simulated and critical sites and exchanges for mutagenesis identified. Vast conformational spaces are efficiently searched as illustrated by the ab initio folding of a model protein of a four beta strand bundle. The genetic algorithm simulation which mimicked important folding constraints as overall hydrophobic packaging and a propensity of the betaphilic residues for trans positions achieved a unique fold. Cooperativity in the beta strand regions and a length of 3-5 for the interconnecting loops was critical. Specific interaction sites were considerably less effective in driving the fold.     

Energetics of membrane protein folding and stability

The critical role of membrane proteins in a myriad of biological and physiological functions has spawned numerous investigations over the past several decades with the long-term goal of identifying the molecular origins and energetic forces that stabilize these proteins within the membrane. Parallel structural and thermodynamics studies on several systems have provided significant insight regarding the driving forces governing folding, assembly, insertion, and translocation of membrane proteins. The present review surveys families of membrane-associated proteins including α-helical and β-barrel structures, viral surface receptors, and pore-forming toxins, citing representative proteins within each of these classes for further scrutiny in terms of structure-function relationships and global conformational stability. This overview presents seminal findings from pioneering studies on the energetics of membrane protein folding and stability to modern techniques that are exploiting the use of molecular genetics and single molecule studies. An overall consensus regarding the molecular origins of membrane protein stability is that a number of intrinsic properties resemble features of soluble proteins, yet there are distinct energetic differences arising from specific intra- and intermolecular interactions within the membrane. The combined efforts from structural, energetics, and dynamics approaches offer unique insights and improve our fundamental understanding of the driving forces dictating membrane protein folding and stability.     

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.     

The folding pathway of reduced lysozyme.

Studies on the mechanism of the glutathione regeneration (Saxena, V.P., and Wetlaufer, D.B. (1970) Biochemistry 9, 5015-5023) of hen egg lysozyme have been carried out. The first two stoichiometric disulfides in lysozyme are formed about 8 times more rapidly than the second two. Almost no enzymic activity is regained until the first two disulfides are formed, thus ruling out an all-or-none mechanism. The disulfide peptides formed early in the regeneration have been isolated and identified. The results show a limited search of folding intermediates, and outline a folding pathway. The early disulfides involve cysteinyl residues III, IV, V, and VI. At the same time cysteinyl residues I, II, VII, and VIII are still reduced, as demonstrated by their isolation as S-alkylated derivatives. At slightly later times a peptide is found which contains the (native) disulfide between cysteinyl residues II and VII. It is likely, but as yet unproven, that formation of disulfide I-VIII completes the cross-linking of lysozyme.     

Sequence determinants of a protein folding pathway

Local hydrophobic collapse of the polypeptide chain and transient long-range interactions in unfolded states of apomyoglobin appear to occur in regions of the amino acid sequence which, upon folding, bury an above-average area of hydrophobic surface. To explore the role of these interactions in protein folding, we prepared and characterized apomyoglobins with compensating point mutations designed to change the average buried surface area in local regions of the sequence, while conserving as much as possible the constitution of the hydrophobic core. The behavior of the mutants in quench-flow experiments to determine the folding pathway was exactly as predicted by the changes in the buried surface area parameter calculated from the amino acid sequence. In addition, spin label experiments with acid-unfolded mutant apomyoglobin showed that the transient long-range contacts that occur in the wild-type protein are abolished in the mutant, while new contacts are observed between areas that now have above-average buried surface area. We conclude that specific groupings of amino acid side-chains, which can be predicted from the sequence, are responsible for early hydrophobic interactions in the first phase of folding in apomyoglobin, and that these early interactions determine the subsequent course of the folding process.     

Theory for the folding and stability of globular proteins

Using lattice statistical mechanics, we develop theory to account for the folding of a heteropolymer molecule such as a protein to the globular and soluble state. Folding is assumed to be driven by the association of solvophobic monomers to avoid solvent and opposed by the chain configurational entropy. Theory predicts a phase transition as a function of temperature or solvent character. Molecules that are too short or too long or that have too few solvophobic residues are predicted not to fold. Globular molecules should have a largely solvophobic core, but there is an entropic tendency for some residues to be "out of place", particularly in small molecules. For long chains, molecules comprised of globular domains are predicted to be thermodynamically more stable than spherical molecules. The number of accessible conformations in the globular state is calculated to be an exceedingly small fraction of the number available to the random coil. Previous estimates of this number, which have motivated kinetic theories of folding, err by many tens of orders of magnitude.     

The folding pathway of T4 lysozyme: an on-pathway hidden folding intermediate

T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. In earlier studies, an amide hydrogen/deuterium exchange pulse-labeling experiment detected a stable submillisecond intermediate that accumulates before the rate-limiting transition state. It involves the formation of structures in both the N and C-terminal regions. However, a native-state hydrogen exchange experiment subsequently detected an equilibrium intermediate that only involves the formation of the C-terminal domain. Here, using stopped-flow circular dichroism and fluorescence, amide hydrogen exchange-folding competition, and protein engineering methods, we re-examined the folding pathway of T4-lysozyme. We found no evidence for the existence of a stable folding intermediate before the rate-limiting transition state at neutral pH. In addition, using native-state hydrogen exchange-directed protein engineering, we created a mimic of the equilibrium intermediate. We found that the intermediate mimic folds with the same rate as the wild-type protein, suggesting that the equilibrium intermediate is an on-pathway intermediate that exists after the rate-limiting transition state.     

A critical residue in the folding pathway of an integral membrane protein

Although a number of common diseases are a direct consequence of membrane protein misfolding, studies of membrane protein folding and misfolding lag well behind those of soluble proteins. Here it is shown that an interfacial residue, Tyr16, of the integral membrane protein diacylglycerol kinase (DAGK) plays a critical role in the folding pathway of this protein. Properly folded Y16C exhibits kinetic parameters and stability similar to wild-type DAGK. However, when unfolded and then allowed to spontaneously fold in the presence of model membranes, Y16C exhibits dramatically lower rates and efficiencies of functional assembly compared to the wild-type protein. The Y16C mutant represents a class of mutations which may be commonly found in disease-related membrane proteins.