Loop length, intramolecular diffusion and protein folding

Exploration of the transition state for tertiary structure formation between an RNA helix and a large structured RNA

Docking of the P1 duplex into the pre-folded core of the Tetrahymena group I ribozyme exemplifies the formation of tertiary interactions in the context of a complex, structured RNA. We have applied Phi-analysis to P1 docking, which compares the effects of modifications on the rate constant for docking (k(dock)) with the effects on the docking equilibrium (K(dock)). To accomplish this we used a single molecule fluorescence resonance energy transfer assay that allows direct determination of the rate constants for formation of thermodynamically favorable, as well as unfavorable, states. Modification of the eight groups of the P1 duplex that make tertiary interactions with the core and changes in solution conditions decrease K(dock) up to 500-fold, whereas k(dock) changes by </=2-fold. The absence of effects on k(dock), both from atomic modifications and global perturbations, strongly suggests that the transition state for docking is early and does not closely resemble the docked state. These results, the slow rate of docking of 3s(-1), and the observation that a modification that is expected to increase the degrees of freedom between the P1 duplex and the ribozyme core accelerates docking, suggest a model in which a kinetic trap(s) slows docking substantially. Nonetheless, urea does not increase k(dock), suggesting that there is little change in the exposed surface area between the trapped, undocked state and the transition state. The findings highlight that urea and temperature dependencies can be inadequate to diagnose the presence of kinetic traps in a folding process. The results described here, combined with previous work, provide an in-depth view of an RNA tertiary structure formation event and suggest that large, highly structured RNAs may have local regions that are misordered.

The unfolding story of three-dimensional domain swapping

Three-dimensional domain swapping is the event by which a monomer exchanges part of its structure with identical monomers to form an oligomer where each subunit has a similar structure to the monomer. The accumulating number of observations of this phenomenon in crystal structures has prompted speculation as to its biological relevance. Domain swapping was originally proposed to be a mechanism for the emergence of oligomeric proteins and as a means for functional regulation, but also to be a potentially harmful process leading to misfolding and aggregation. We highlight experimental studies carried out within the last few years that have led to a much greater understanding of the mechanism of domain swapping and of the residue- and structure-specific features that facilitate the process. We discuss the potential biological implications of domain swapping in light of these findings.

Low free energy cost of very long loop insertions in proteins

Long insertions into a loop of a folded host protein are expected to have destabilizing effects because of the entropic cost associated with loop closure unless the inserted sequence adopts a folded structure with amino- and carboxy-termini in close proximity. A loop entropy reduction screen based on this concept was used in an attempt to retrieve folded sequences from random sequence libraries. A library of long random sequences was inserted into a loop of the SH2 domain, displayed on the surface of M13 phage, and the inserted sequences that did not disrupt SH2 function were retrieved by panning using beads coated with a phosphotyrosine containing SH2 peptide ligand. Two sequences of a library of 2 x 10(8) sequences were isolated after multiple rounds of panning, and were found to have recovery levels similar to the wild-type SH2 domain and to be relatively intolerant to further mutation in PCR mutagenesis experiments. Surprisingly, although these inserted sequences exhibited little nonrandom structure, they do not significantly destabilize the host SH2 domain. Additional insertion variants recovered at lower levels in the panning experiments were also found to have a minimal effect on the stability and peptide-binding function of the SH2 domain. The additional level of selection present in the panning experiments is likely to involve in vivo folding and assembly, as there was a rough correlation between recovery levels in the phage-panning experiments and protein solubility. The finding that loop insertions of 60-80 amino acids have minimal effects on SH2 domain stability suggests that the free energy cost of inserting long loops may be considerably less than polymer theory estimates based on the entropic cost of loop closure, and, hence, that loop insertion may have provided an evolutionary route to multidomain protein structures.

The topomer search model: A simple, quantitative theory of two-state protein folding kinetics

Most small, single-domain proteins fold with the uncomplicated, single-exponential kinetics expected for diffusion on a smooth energy landscape. Despite this energetic smoothness, the folding rates of these two-state proteins span a remarkable million-fold range. Here, we review the evidence in favor of a simple, mechanistic description, the topomer search model, which quantitatively accounts for the broad scope of observed two-state folding rates. The model, which stipulates that the search for those unfolded conformations with a grossly correct topology is the rate-limiting step in folding, fits observed rates with a correlation coefficient of approximately 0.9 using just two free parameters. The fitted values of these parameters, the pre-exponential attempt frequency and a measure of the difficulty of ordering an unfolded chain, are consistent with previously reported experimental constraints. These results suggest that the topomer search process may dominate the relative barrier heights of two-state protein-folding reactions.

Mechanism of fast protein folding

An explosion of in vitro experimental data on the folding of proteins has revealed many examples of folding in the millisecond or faster timescale, often occurring in the absence of stable intermediate states. We review experimental methods for measuring fast protein folding kinetics, and then discuss various analytical models used to interpret these data. Finally, we classify general mechanisms that have been proposed to explain fast protein folding into two catagories, heterogeneous and homogeneous, reflecting the nature of the transition state. One heterogeneous mechanism, the diffusion-collision mechanism, can be used to interpret experimental data for a number of proteins.

Experimental evaluation of topological parameters determining protein-folding rates

Recent work suggests that structural topology plays a key role in determining protein-folding rates and pathways. The refolding rates of small proteins that fold without intermediates are found to correlate with simple structural parameters such as relative contact order, long-range order, or the fraction of short-range contacts. To test and evaluate the role of structural topology experimentally, a set of circular permutants of the ribosomal protein S6 from Thermus thermophilus was analyzed. Despite a wide range of relative contact order, the permuted proteins all fold with similar rates. These results suggest that alternative topological parameters may better describe the role of topology in protein-folding rates.

Protein folding kinetics beyond the phi value: using multiple amino acid substitutions to investigate the structure of the SH3 domain folding transition state

The SH3 domain folding transition state structure contains two well-ordered turn regions, known as the diverging turn and the distal loop. In the Src SH3 domain transition state, these regions are stabilized by a hydrogen bond between Glu30 in the diverging turn and Ser47 in the distal loop. We have examined the effects on folding kinetics of amino acid substitutions at the homologous positions (Glu24 and Ser41) in the Fyn SH3 domain. In contrast to most other folding kinetics studies which have focused primarily on non-disruptive substitutions with Ala or Gly, here we have examined the effects of substitutions with diverse amino acid residues. Using this approach, we demonstrate that the transition state structure is generally tolerant to amino acid substitutions. We also uncover a unique role for Ser at position 41 in facilitating folding of the distal loop, which can only be replicated by Asp at the same position. Both these residues appear to accelerate folding through the formation of short-range side-chain to backbone hydrogen bonds. The folding of the diverging turn region is shown to be driven primarily by local interactions. The diverging turn and distal loop regions are found to interact in the transition state structure, but only in the context of particular mutant backgrounds. This work demonstrates that studying the effects of a variety of amino acid substitutions on protein folding kinetics can provide unique insights into folding mechanisms which cannot be obtained by standard Phi value analysis.

Application of the diffusion-collision model to the folding of three-helix bundle proteins

The diffusion-collision model has been successful in explaining many features of protein folding kinetics, particularly for helical proteins. In the model the folding reaction is described in terms of coupled chemical kinetic (Master) equations of coarse grained entities, called microdomains. Here, the diffusion-collision model is applied to compute the folding kinetics of four three-helix bundle proteins, all of which fold on a time scale of tens of microseconds and appear to have two-state folding. The native structure and the stability of the helical microdomains are used to determine the parameters of the model. The formulation allows computation of the overall rate and determination of the importance of kinetic intermediates. The proteins considered are the B domain of protein A (1BDC), the Engrailed Homeodomain (1ENH), the peripheral sub-unit-binding domain (1EBD C-chain) and the villin headpiece subdomain (1VII). The results for the folding time of protein A, the Engrailed Homeodomain, and 1EBD C-chain are in agreement with experiment, while 1VII is not stable in the present model. In the three proteins that are stable, two-state folding is predicted by the diffusion-collision model. This disagrees with published assertions that multistate kinetics would be obtained from the model. The contact order prediction agrees with experiment for protein A, but yields values that are a factor of 40, 30 and 15 too slow for 1ENH, 1EBD C-chain and 1VII. The effect of mutants on folding is described for protein A and it is demonstrated that significant intermediate concentrations (i.e. deviation from two-state folding) can occur if the stability of some of the helical microdomains is increased. A linear relationship between folding time and the length of the loop between helices B and C in protein A is demonstrated; this is not evident in the contact order description.

Folding rate prediction using total contact distance

Linear regression analysis found that either contact order (CO) or long-range order (LRO) parameter has a significant correlation with the logarithms of folding rates. This suggests that sequence separation per contact and total number of contacts are both important in determining the rate of folding. Here, the two factors are incorporated into a new parameter, total contact distance (TCD). Using a database of 28 two-state or weakly three-state folding proteins, TCD is found to be the most accurate among the three parameters (CO, LRO, and TCD) in terms of correlation and prediction. It provides even more accurate prediction than the best neural network results with two descriptors (contact order and stability per residue). The improvement is achieved in all three-structural classes (all alpha, beta, and mixed). The accuracy of total contact distance in predicting folding rates is essentially unchanged if "short"-ranged contacts (absolute value of i - j < or = 14) are not included in calculation. Thus, only long-range contacts with a sequence separation of more than 14 residues are important in determining the rate of folding. This is consistent with the results from the long-range order parameter. One of the significant outliers in prediction is found to be associated with the only protein in the database that involves nonlocal disulfide bonds. Removing the protein leads to a correlation coefficient of 0.89 between experimental observed and predicted folding rates in jackknife cross validation. The corresponding values for CO and LRO are 0.71 and 0.80, respectively.

Folding of circular permutants with decreased contact order: general trend balanced by protein stability

To examine the influence of contact order and stability on the refolding rate constant for two-state proteins, we have analysed the folding kinetics of the small beta-alpha-beta protein S6 and two of its circular permutants with relative contact orders of 0.19, 0.15 and 0.12. Data reveal a small but significant increase of the refolding rate constant (log k(f)) with decreasing contact order. At the same time, the decreased contact order is correlated to losses in global stability and alterations of the folding nucleus. When the differences in stability are accounted for by addition of Na2SO4 or by comparison of the folding kinetics at the transition mid-point, the dependence between log k(f) and contact order becomes stronger and follows the general correlation for two-state proteins. The observation emphasizes the combined action of topology and stability in controlling the rate constant of protein folding.

Role of loops in the folding and stability of yeast phosphoglycerate kinase

Yeast phosphoglycerate kinase (yPGK) is a monomeric two domain protein used as folding model representative of large proteins. We inserted short unstructured sequences (four Gly or four Thr) into the connections between secondary structure elements and studied the consequences of these insertions on the folding process and stability of yPGK. All the mutated proteins can refold efficiently. The effect per residue on stability is larger for the first inserted residue. Insertion in two long betaalpha loops (at residue positions 71 and 129) is more destabilizing than an insertion in a short alphabeta loop (at residue position 89) located on the opposite side of the N-terminal domain. The effect on stability is mainly due to a large increase of the unfolding rate rather than a decrease of the folding rate. This suggests that these connections between secondary structure elements do not play an active role in directing the folding process. Insertion into the short alphabeta loop (position 89) has limited effects on stability and results in the detection of a kinetic phase not previously seen with the wild-type protein, suggesting that insertions in this particular loop do qualitatively affect the folding process without a large effect on folding efficiency. For the two long betaalpha loops (positions 71 and 129) located in the inner surface of the N-terminal domain, the effects on stability are possibly associated with decoupling of the two domains as observed by differential scanning calorimetry during thermal unfolding.

Protein folding theory: from lattice to all-atom models

This review focuses on recent advances in understanding protein folding kinetics in the context of nucleation theory. We present basic concepts such as nucleation, folding nucleus, and transition state ensemble and then discuss recent advances and challenges in theoretical understanding of several key aspects of protein folding kinetics. We cover recent topology-based approaches as well as evolutionary studies and molecular dynamics approaches to determine protein folding nucleus and analyze other aspects of folding kinetics. Finally, we briefly discuss successful all-atom Monte-Carlo simulations of protein folding and conclude with a brief outlook for the future.

Mechanism of Tet repressor induction by tetracyclines: length compensates for sequence in the alpha8-alpha9 loop

Natural Tet repressor (TetR) variants are alpha-helical proteins bearing a large loop between helices 8 and 9, which is variable in sequence and length. We have deleted this loop consisting of 14 amino acid residues in TetR(D) and rebuilt it stepwise with up to 42 alanine residues. All except the mutant with the longest alanine loop show wild-type repression, but none is inducible with tetracycline. This demonstrates the importance of the alpha8-alpha9 loop and its amino acid sequence for induction. The induction efficiencies increase with loop length, when the more tightly binding inducer anhydrotetracycline is used. The largest increase of inducibility was observed for TetR mutants with loop lengths between eight and 17 alanine residues. Since loop residues Asp/Glu157 and Arg158 are conserved in the natural TetR sequence variants, we constructed a mutant in which all other residues of the loop were replaced by alanine. This mutant exhibits increased anhydrotetracycline induction compared to the corresponding alanine variant. Thus, these residues are important for induction. Binding constants for the anhydrotetracycline-TetR interaction are below the detection level of 10(5) M(-1) for the mutant with a loop of two alanine residues and increase sharply until a loop size of ten residues is reached. TetR variants with longer loops have similar anhydrotetracycline-binding constants, ranging between 2.6 x 10(9) M(-1) and 8.0 x 10(9) M(-1), about 500-fold lower than wild-type TetR. The increase of the affinity occurs at shorter loop lengths than that of inducibility. We conclude that the induction defect of the polyalanine variants arises from two increments: (i) the loop must have a minimal length-to allow efficient inducer binding; (ii) the loop must structurally participate in the conformational change associated with induction.

The cost of exposing a hydrophobic loop and implications for the functional role of 4.5 S RNA in the Escherichia coli signal recognition particle

The signal recognition particle (SRP) is an RNA-protein complex that directs ribosomes to the rough endoplasmic reticulum membrane by binding to targeting signals found on the nascent chain of proteins destined for export to the endoplasmic reticulum. We found evidence from studies with fragments of the protein component of the Escherichia coli SRP that a long hydrophobic loop (the so-called "finger loop") is detrimental to the stability of its signal peptide-binding domain, the M domain. This hydrophobic loop is highly conserved and thus may have a critical role in the function of the SRP. Given our previously reported evidence that 4.5 S RNA stabilizes the tertiary fold of the M domain (Zheng, N., and Gierasch, L. M. (1997) Mol. Cell 1, 79-87), we now propose that the functional requirement for 4.5 S RNA resides in its ability to counteract the destabilizing influence of the finger loop.

Keeping it in the family: folding studies of related proteins

Investigators have recently turned to studies of protein families to shed light on the mechanism of protein folding. In small proteins for which detailed analysis has been performed, recent studies show that transition-state structure is generally conserved. The number and structures of populated folding intermediates have been found to vary in homologous families of larger (greater than 100-residue) proteins, reflecting a balance of local and global interactions.

Intracellular folding pathway of the cystine knot-containing glycoprotein hormone alpha-subunit

Three of the five disulfide bonds in the glycoprotein hormone alpha-subunit (GPH-alpha) form a cystine knot motif that stabilizes a three-loop antiparallel structure. Previously, we described a mutant (alpha(k)) that contained only the three knot disulfide bonds and demonstrated that the cystine knot was necessary and sufficient for efficient GPH-alpha folding and secretion. In this study, we used alpha(k) as a model to study the intracellular GPH-alpha folding pathway. Cystine knot formation proceeded through a 1-disulfide intermediate that contained the 28-82 disulfide bond. Formation of disulfide bond 10-60, then disulfide bond 32-84, followed the formation of 28-82. Whether the two non-cystine knot bonds 7-31 and 59-87 could form independent of the knot was also tested. Disulfide bond 7-31 formed rapidly, whereas 59-87 did not form when all cysteine residues of the cystine knot were converted to alanine, suggesting that 7-31 forms early in the folding pathway and that 59-87 forms during or after cystine knot formation. Finally, loop 2 of GPH-alpha has been shown to be very flexible, suggesting that loop 2 does not actively drive GPH-alpha folding. To test this, we replaced residues 36-55 in the flexible loop 2 with an artificially flexible glycine chain. Consistent with our hypothesis, folding and secretion were unaffected when loop 2 was replaced with the glycine chain. Based on these findings, we describe a model for the intracellular folding pathway of GPH-alpha and discuss how these findings may provide insight into the folding mechanisms of other cystine knot-containing proteins.

The SH3-fold family: experimental evidence and prediction of variations in the folding pathways

To investigate the relationships between protein topology, amino acid sequence and folding mechanisms, the folding transition state of the Sso7d protein has been characterised both experimentally and theoretically. Although Sso7d protein has a similar topology to that of the SH3 domains, the structure of its transition state is different from that of alpha-spectrin and src SH3 domains previously studied. The folding algorithm, Fold-X, including an energy function with specific sequence features, accounts for these differences and reproduces with a good agreement the set of experimental phi(double dagger-U) values obtained for the three proteins. Our analysis shows that taking into account sequence features underlying protein topology is critical for an accurate prediction of the folding process.

Multiple roles of prolyl residues in structure and folding

To explore the ways that proline residues may influence the conformational options of a polypeptide backbone, we have characterized Pro-->Ala mutants of cellular retinoic acid-binding protein I (CRABP I). While all three Xaa-Pro bonds are in the trans conformation in the native protein and the equilibrium stability of each mutant is similar to that of the parent protein, each has distinct effects on folding and unfolding kinetics. The mutation of Pro105 does not alter the kinetics of folding of CRABP I, which indicates that the flexible loop containing this residue is passive in the folding process. By contrast, replacement of Pro85 by Ala abolishes the observable slow phase of folding, revealing that correct configuration of the 84-85 peptide bond is prerequisite to productive folding. Substitution of Pro39 by Ala yields a protein that folds and unfolds more slowly. Removal of the conformational constraint imposed by the proline ring likely raises the transition state barrier by increasing the entropic cost of narrowing the conformational ensemble. Additionally, the Pro-->Ala mutation removes a helix-termination signal that is important for efficient folding to the native state.

Long-range order in the src SH3 folding transition state

One of the outstanding questions in protein folding concerns the degree of heterogeneity in the folding transition state ensemble: does a protein fold via a large multitude of diverse "pathways," or are the elements of native structure assembled in a well defined order? Herein, we build on previous point mutagenesis studies of the src SH3 by directly investigating the association of structural elements and the loss of backbone conformational entropy during folding. Double-mutant analysis of polar residues in the distal beta-hairpin and the diverging turn indicates that the hydrogen bond network between these elements is largely formed in the folding transition state. A 10-glycine insertion in the n-src loop (which connects the distal hairpin and the diverging turn) and a disulfide crosslink at the base of the distal beta-hairpin exclusively affect the folding rate, showing that these structural elements are nearly as ordered in the folding transition state as in the native state. In contrast, crosslinking the base of the RT loop or the N and C termini dramatically slows down the unfolding rate, suggesting that dissociation of the termini and opening of the RT loop precede the rate-limiting step in unfolding. Taken together, these results suggest that essentially all conformations in the folding transition state ensemble have the central three-stranded beta-sheet formed, indicating that, for the src homology 3 domain, there is a discrete order to structure assembly during folding.

Transition-state structure as a unifying basis in protein-folding mechanisms: contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Phi-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and alpha-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Prediction of protein-folding mechanisms from free-energy landscapes derived from native structures

Guided by recent experimental results suggesting that protein-folding rates and mechanisms are determined largely by native-state topology, we develop a simple model for protein folding free-energy landscapes based on native-state structures. The configurations considered by the model contain one or two contiguous stretches of residues ordered as in the native structure with all other residues completely disordered; the free energy of each configuration is the difference between the entropic cost of ordering the residues, which depends on the total number of residues ordered and the length of the loop between the two ordered segments, and the favorable attractive interactions, which are taken to be proportional to the total surface area buried by the ordered residues in the native structure. Folding kinetics are modeled by allowing only one residue to become ordered/disordered at a time, and a rigorous and exact method is used to identify free-energy maxima on the lowest free-energy paths connecting the fully disordered and fully ordered configurations. The distribution of structure in these free-energy maxima, which comprise the transition-state ensemble in the model, are reasonably consistent with experimental data on the folding transition state for five of seven proteins studied. Thus, the model appears to capture, at least in part, the basic physics underlying protein folding and the aspects of native-state topology that determine protein-folding mechanisms.

Universally conserved positions in protein folds: reading evolutionary signals about stability, folding kinetics and function

Here, we provide an analysis of molecular evolution of five of the most populated protein folds: immunoglobulin fold, oligonucleotide-binding fold, Rossman fold, alpha/beta plait, and TIM barrels. In order to distinguish between "historic", functional and structural reasons for amino acid conservations, we consider proteins that acquire the same fold and have no evident sequence homology. For each fold we identify positions that are conserved within each individual family and coincide when non-homologous proteins are structurally superimposed. As a baseline for statistical assessment we use the conservatism expected based on the solvent accessibility. The analysis is based on a new concept of "conservatism-of-conservatism". This approach allows us to identify the structural features that are stabilized in all proteins having a given fold, despite the fact that actual interactions that provide such stabilization may vary from protein to protein. Comparison with experimental data on thermodynamics, folding kinetics and function of the proteins reveals that such universally conserved clusters correspond to either: (i) super-sites (common location of active site in proteins having common tertiary structures but not function) or (ii) folding nuclei whose stability is an important determinant of folding rate, or both (in the case of Rossman fold). The analysis also helps to clarify the relation between folding and function that is apparent for some folds.

Protein redesign

It has been demonstrated that, for a number of proteins, it is possible to dramatically alter the connectivities between elements of secondary structure. Remarkably large loop insertions are tolerated and many redesigns have generated proteins that successfully fold to stable, active structures. Some redesigns have been entirely the choice of the investigators, whereas others have incorporated a randomization and selection step to identify optimal sequences. These studies have provided basic guidelines for the rational manipulation of protein structure and stability, they have allowed the dissection of folding pathways and they have generated proteins with the potential for practical therapeutic applications.

Matching theory and experiment in protein folding

There has been considerable progress made over the past year in linking experimental and theoretical approaches to protein folding. Recent results from several independent lines of investigation suggest that protein folding mechanisms and landscapes are largely determined by the topology of the native state and are relatively insensitive to details of the interatomic interactions. This dependence on low-resolution structural features, rather than high-resolution detail, suggests that it should be possible to describe the fundamental physics of the folding process using relatively low-resolution models. Recent experiments have set benchmarks for testing new models and progress has been made in developing theoretical models for interpreting and predicting experimental results.

Using loop length variants to dissect the folding pathway of a four-helix-bundle protein

Rop is a four-helix-bundle protein formed by the association of two helix-loop-helix monomers. The short helix-connecting loop was replaced with a series of polyglycine linkers of increasing length. These mutant proteins all appear to fold via the same general mechanism as that of the wild-type protein, even at the longest loop lengths. Replacement of the wild-type two-residue loop (Asp-Ala) with a (Gly-Gly) linker accelerates both unfolding and refolding rates. These changes in folding and unfolding kinetics likely indicate an alteration in the energy of the transition state. As the length of the glycine linker is further increased, the unfolding rate increases while the refolding rates decrease. The influence of loop length is not limited to these rates, but also impacts upon the stability of the folding intermediate. These dependences underscore the importance of loop closure and help refine the model for Rop's folding, implicating a dimeric intermediate involving hairpin formation. These observations show that loop alteration may be useful as a general technique for dissecting protein folding pathways.

High populations of non-native structures in the denatured state are compatible with the formation of the native folded state

The structures of the denatured states of the spectrin SH3 domain and a mutant designed to have a non-native helical tendency at the N terminus have been analyzed under mild acidic denaturing conditions by nuclear magnetic resonance methods with improved resolution. The wild-type denatured state has little residual structure. However, the denatured state of the mutant has an approximately 50% populated helical structure from residues 2 to 14, a region that forms part of the beta-sheet structure in the folded state. Comparison with a peptide corresponding to the same sequence shows that the helix is stabilized in the whole domain, likely by non-local interactions with other parts of the protein as suggested by changes in a region far from the mutated sequence. These results demonstrate that high populations of non-native secondary structure elements in the denatured state are compatible with the formation of the native folded structure.

A test of the relationship between sequence and structure in proteins: excision of the heme binding site in apocytochrome b5

The water-soluble domain of rat hepatic holocytochrome b5 is an alphabeta protein containing elements of secondary structure in the sequence beta1-alpha1-beta4-beta3-alpha2-alpha3-beta5- alpha4-alpha5-beta2-alpha6. The heme group is enclosed by four helices, a2, a3, a4, and a5. To test the hypothesis that a small b hemoprotein can be constructed in two parts, one forming the heme site, the other an organizing scaffold, a protein fragment corresponding to beta1-alpha1-beta4-beta3-lambda-beta2-alpha6 was prepared, where lambda is a seven-residue linker bypassing the heme binding site. The fragment ("abridged b5") was found to contain alpha and beta secondary structure by circular dichroism spectroscopy and tertiary structure by Trp fluorescence emission spectroscopy. NMR data revealed a species with spectral properties similar to those of the full-length apoprotein. This folded form is in slow equilibrium on the chemical shift time scale with other less folded species. Thermal denaturation, as monitored by circular dichroism, absorption, and fluorescence spectroscopy, as well as size-exclusion chromatography-fast protein liquid chromatography (SEC-FPLC), confirmed the coexistence of at least two distinct conformational ensembles. It was concluded that the protein fragment is capable of adopting a specific fold likely related to that of cytochrome b5, but does not achieve high thermodynamic stability and cooperativity. Abridged b5 demonstrates that the spliced sequence contains the information necessary to fold the protein. It suggests that the dominating influence to restrict the conformational space searched by the chain is structural propensities at a local level rather than internal packing. The sequence also holds the properties necessary to generate a barrier to unfolding.

Obligatory steps in protein folding and the conformational diversity of the transition state

We have analyzed the existence of obligatory steps in the folding reaction of the alpha-spectrin SH3 domain by mutating Asp 48 (D48G), which is at position i+3 of an isolated two-residue type II' beta-turn. Calorimetry and X-ray analysis show an entropic stabilizing effect resulting from local changes at the dihedral angles of the beta-turn. Kinetic analysis of D48G shows that this beta-turn is fully formed in the transition state, while there is no evidence of its formation in an isolated fragment. Introduction of several mutations in the D48G protein reveals that the local stabilization has not significantly altered the transition state ensemble. All these results, together with previous analysis of other alpha-spectrin and src SH3 mutants, indicate that: (i) in the folding reaction there could be obligatory steps which are not necessarily part of the folding nucleus; (ii) transition state ensembles in beta-sheet proteins could be quite defined and conformationally restricted ('mechanic folding nucleus'); and (iii) transition state ensembles in some proteins could be evolutionarily conserved.

A protein engineering analysis of the transition state for protein folding: simulation in the lattice model

BACKGROUND

Protein engineering has been used extensively to evaluate the properties of transition states in protein folding. Although the method has proved useful, its limitations and the details of interpretation of the obtained results remain largely unexplored.

RESULTS

Lattice model simulations are used to test and verify the protein engineering analysis of the transition state in protein folding. It is shown that in some cases - but not always - this method is able to determine the transition state with reasonable accuracy. Limitations of protein engineering are revealed and analyzed. In particular, the change in non-native interactions as a result of mutations is shown to influence the results of the protein engineering analysis. Furthermore, the temperature dependencies of phi values (which are a measure of the participation of a residue in the transition state) and the character of the transition state ensemble are studied. It is shown that as a general trend phi values decrease when the temperature decreases, a finding consistent with recent experimental results. Our analysis suggests that this trend results primarily from the formation of some contacts (native and non-native) in the unfolded state at a lower temperature, when the barrier for folding is energetic.

CONCLUSIONS

Our analysis helps to interpret the results of protein engineering and allows observed &phi; values to be directly related to structural features of the unfolded state, the transition state and the native state.