Towards the complete structural characterization of a protein folding pathway: the structures of the denatured, transition and native states for the association/folding of two complementary fragments of cleaved chymotrypsin inhibitor 2. Direct evidence for a nucleation-condensation mechanism

Autonomous Sequences in Myoglobin Emerging from X-ray Structure of Holomyoglobin

The proposed continuous folding structure units are fundamental to analyze protein structure. Here, we could elucidate for the first time two types of hydrophobic core networks in apomyoglobin using continuous folding structure units. In myoglobin, two autonomous sequences emerged clearly. We could thus characterize the autonomous sequences using well-defined hydrophobic core networks within respective semifolds. A hydrophobic core is defined as a pair of topology-local hydrophobic amino acids in different folding structures. Hydrophobic core formation is indispensable to stabilize the different folding structures via an efficient hydrophobic interaction. Autonomous sequences in myoglobin encode tertiary structure information for semifolds. These sequences fold autonomously into small sets of continuous folding structure units to grow separate semifolds on each separate framework. The autonomous sequence can be defined as the local sequence assigned to the small set of continuous folding structure units. They create the discrete hydrophobic region in a semifold by assembly of their hydrophobic regions. Semifolds were characterized by discrete hydrophobic regions stabilized by respective type I hydrophobic core networks, which were present within each semifold. The discrete hydrophobic region of a semifold propagated itself with that of a different semifold by hydrophobic interactions in type II hydrophobic core network, which was present between different semifolds, as observed by the X-ray structures of semifolds. The most significant feature of semifolds in apomyoglobin was that they could be verified by the X-ray structure of holomyoglobin regardless of the instability of folds characteristic to autonomous sequence fragments. This work presents the first description of autonomous sequences.

The folding of a family of three-helix bundle proteins: spectrin R15 has a robust folding nucleus, unlike its homologous neighbours

Three homologous spectrin domains have remarkably different folding characteristics. We have previously shown that the slow-folding R16 and R17 spectrin domains can be altered to resemble the fast folding R15, in terms of speed of folding (and unfolding), landscape roughness and folding mechanism, simply by substituting five residues in the core. Here we show that, by contrast, R15 cannot be engineered to resemble R16 and R17. It is possible to engineer a slow-folding version of R15, but our analysis shows that this protein neither has a rougher energy landscape nor does change its folding mechanism. Quite remarkably, R15 appears to be a rare example of a protein with a folding nucleus that does not change in position or in size when its folding nucleus is disrupted. Thus, while two members of this protein family are remarkably plastic, the third has apparently a restricted folding landscape.

Take home lessons from studies of related proteins

The 'Fold Approach' involves a detailed analysis of the folding of several topologically, structurally and/or evolutionarily related proteins. Such studies can reveal determinants of the folding mechanism beyond the gross topology, and can dissect the residues required for folding from those required for stability or function. While this approach has not yet matured to the point where we can predict the native conformation of any polypeptide chain in silico, it has been able to highlight, amongst others, the specific residues that are responsible for nucleation, pathway malleability, kinetic intermediates, chain knotting, internal friction and Paracelsus switches. Some of the most interesting discoveries have resulted from the attempt to explain differences between homologues.

Protein stability, flexibility and function

Proteins rely on flexibility to respond to environmental changes, ligand binding and chemical modifications. Potentially, a perturbation that changes the flexibility of a protein may interfere with its function. Millions of mutations have been performed on thousands of proteins in quests for a delineation of the molecular details of their function. Several of these mutations interfered with the binding of a specific ligand with a concomitant effect on the stability of the protein scaffold. It has been ambiguous and not straightforward to recognize if any relationships exist between the stability of a protein and the affinity for its ligand. In this review, we present examples of proteins where changes in stability results in changes in affinity and of proteins where stability and affinity are uncorrelated. We discuss the possibility for a relationship between stability and binding. From the data presented is it clear that there are specific sites (flexibility hotspots) in proteins that are important for both binding and stability. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Deciphering the folding transition state structure and denatured state properties of nucleophosmin C-terminal domain

Nucleophosmin (NPM1), one of the most abundant nucleolar proteins, is a frequent target of oncogenic mutations in acute myeloid leukaemia (AML). Mutation-induced changes at the C-terminal domain of NPM1 (Cter-NPM1) compromise its stability and cause the aberrant translocation of NPM1 to the cytosol. Hence, this protein represents a suitable candidate to investigate the relations between folding and disease. Since Cter-NPM1 folds via a compact denatured state, stabilization of the folded state of the mutated variants demands detailed structural information on both the native and denatured states. Here, we present the characterization of the complete folding pathway of Cter-NPM1 and provide molecular details for both the transition and the denatured states. The structure of the transition state was assessed by Phi-value analysis, whereas residual structure in the denatured state was mapped by evaluating the effect of mutations as modulated by conditions promoting denatured state compaction. Data reveal that folding of Cter-NPM1 proceeds via an extended nucleus and that the denatured state retains significant malleable structure at the interface between the second and third helices. Our observations constitute the essential prerequisite for structure-based drug-design studies, aimed at identifying molecules that may rescue pathological NPM1 mutants by stabilizing the native-like state.

In silico protein fragmentation reveals the importance of critical nuclei on domain reassembly

Protein complementation assays (PCAs) based on split protein fragments have become powerful tools that facilitate the study and engineering of intracellular protein-protein interactions. These assays are based on the observation that a given protein can be split into two inactive fragments and these fragments can reassemble into the original properly folded and functional structure. However, one experimentally observed limitation of PCA systems is that the folding of a protein from its fragments is dramatically slower relative to that of the unsplit parent protein. This is due in part to a poor understanding of how PCA design parameters such as split site position in the primary sequence and size of the resulting fragments contribute to the efficiency of protein reassembly. We used a minimalist on-lattice model to analyze how the dynamics of the reassembly process for two model proteins was affected by the location of the split site. Our results demonstrate that the balanced distribution of the "folding nucleus," a subset of residues that are critical to the formation of the transition state leading to productive folding, between protein fragments is key to their reassembly.

Protein reconstitution and three-dimensional domain swapping: benefits and constraints of covalency

The phenomena of protein reconstitution and three-dimensional domain swapping reveal that highly similar structures can be obtained whether a protein is comprised of one or more polypeptide chains. In this review, we use protein reconstitution as a lens through which to examine the range of protein tolerance to chain interruptions and the roles of the primary structure in related features of protein structure and folding, including circular permutation, natively unfolded proteins, allostery, and amyloid fibril formation. The results imply that noncovalent interactions in a protein are sufficient to specify its structure under the constraints imposed by the covalent backbone.

High-affinity fragment complementation of a fibronectin type III domain and its application to stability enhancement

The tenth fibronectin type III (FN3) domain of human fibronectin (FNfn10), a prototype of the ubiquitous FN3 domain, is a small, monomeric beta-sandwich protein. In this study, we have bisected FNfn10 in each loop to generate a total of six fragment pairs. We found that fragment pairs bisected at multiple loops of FNfn10 show complementation in vivo as tested with a yeast two-hybrid system. The dissociation constant of these fragment pairs determined in vitro were as low as 3 nM, resulting in one of the tightest fragment complementation systems reported so far. Furthermore, we show that the affinity of fragment complementation is correlated with the stability of the uncut parent protein. Exploring this correlation, we screened a yeast two-hybrid library of one fragment and identified mutations that suppress the effect of a destabilizing mutation in the other fragment. One of the identified mutations significantly increased the stability of the uncut wild-type protein, proving that fragment complementation can be used as a novel strategy for the selection of proteins with enhanced stability.

Are residues in a protein folding nucleus evolutionarily conserved?

Protein is the working molecule of the cell, and evolution is the hallmark of life. It is important to understand how protein folding and evolution influence each other. Several studies correlating experimental measurement of residue participation in folding nucleus and sequence conservation have reached different conclusions. These studies are based on assessment of sequence conservation at folding nucleus sites using entropy or relative entropy measurement derived from multiple sequence alignment. Here we report analysis of conservation of folding nucleus using an evolutionary model alternative to entropy-based approaches. We employ a continuous time Markov model of codon substitution to distinguish mutation fixed by evolution and mutation fixed by chance. This model takes into account bias in codon frequency, bias-favoring transition over transversion, as well as explicit phylogenetic information. We measure selection pressure using the ratio omega of synonymous versus non-synonymous substitution at individual residue site. The omega-values are estimated using the PAML method, a maximum-likelihood estimator. Our results show that there is little correlation between the extent of kinetic participation in protein folding nucleus as measured by experimental phi-value and selection pressure as measured by omega-value. In addition, two randomization tests failed to show that folding nucleus residues are significantly more conserved than the whole protein, or the median omega value of all residues in the protein. These results suggest that at the level of codon substitution, there is no indication that folding nucleus residues are significantly more conserved than other residues. We further reconstruct candidate ancestral residues of the folding nucleus and suggest possible test tube mutation studies for testing folding behavior of ancient folding nucleus.

Topological determinants of protein folding

The folding of many small proteins is kinetically a two-state process that represents overcoming the major free-energy barrier. A kinetic characteristic of a conformation, its probability to descend to the native state domain in the amount of time that represents a small fraction of total folding time, has been introduced to determine to which side of the free-energy barrier a conformation belongs. However, which features make a protein conformation on the folding pathway become committed to rapidly descending to the native state has been a mystery. Using two small, well characterized proteins, CI2 and C-Src SH3, we show how topological properties of protein conformations determine their kinetic ability to fold. We use a macroscopic measure of the protein contact network topology, the average graph connectivity, by constructing graphs that are based on the geometry of protein conformations. We find that the average connectivity is higher for conformations with a high folding probability than for those with a high probability to unfold. Other macroscopic measures of protein structural and energetic properties such as radius of gyration, rms distance, solvent-accessible surface area, contact order, and potential energy fail to serve as predictors of the probability of a given conformation to fold.

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.

Reconstitution of a native-like SH2 domain from disordered peptide fragments examined by multidimensional heteronuclear NMR

The N-terminal SH2 domain from the p85alpha subunit of phosphatidylinositol 3' kinase is cleaved specifically into 9- and 5-kD fragments by limited proteolytic digestion with trypsin. The noncovalent SH2 domain complex and its constituent tryptic peptides have been investigated using high-resolution heteronuclear magnetic resonance (NMR). These studies have established the viability of the SH2 domain as a fragment complementation system. The individual peptide fragments are predominantly unstructured in solution. In contrast, the noncovalent 9-kD + 5-kD complex shows a native-like (1)H-(15)N HSQC spectrum, demonstrating that the two fragments fold into a native-like structure on binding. Chemical shift analysis of the noncovalent complex compared to the native SH2 domain reveals that the highest degree of perturbation in the structure occurs at the cleavage site within a flexible loop and along the hydrophobic interface between the two peptide fragments. Mapping of these chemical shift changes on the structure of the domain reveals changes consistent with the reduction in affinity for the target peptide ligand observed in the noncovalent complex relative to the intact protein. The 5-kD fragment of the homologous Src protein is incapable of structurally complementing the p85 9-kD fragment, either in complex formation or in the context of the full-length protein. These high-resolution structural studies of the SH2 domain fragment complementation features establish the suitability of the system for further protein-folding and design studies.

Conversion of two-state to multi-state folding kinetics on fusion of two protein foldons

Chymotrypsin inhibitor 2 (CI2) is the archetypal single-foldon protein that folds in simple two-state kinetics without the accumulation of a folding intermediate. To model the effects of fusion of single foldons to give a multi-foldon protein, we engineered a "double-CI2" protein, in which another CI2 polypeptide was inserted into the loop region of the parent CI2. CD and HSQC spectra demonstrated that while the double-CI2 protein adopted two kinds of native conformations, CI2-like structure was almost preserved in both the domains of double-CI2. In the folding kinetic studies, double-CI2 exhibited a remarkable rollover of the observed folding rates at low denaturant concentrations, indicating that double-CI2 accumulated a kinetic folding intermediate. The different folding mechanisms between WT-CI2 and double-CI2 support the present view that protein size or number of domains is an important determinant for formation of folding intermediates.

THEFASTPROTEINFOLDINGPROBLEM

During protein folding, many of the events leading to secondary and tertiary structure occur in milliseconds or faster. Modern nuclear magnetic resonance and laser detection techniques, coupled with fast initiation of the folding reaction, are probing these events in great detail. Theory, ranging from analytical models to molecular dynamics calculations, is beginning to match up with experiment. As a result, timescales, from such elementary steps as the addition of a residue to a helix to strange kinetics of collapsing protein backbones, can now be measured and interpreted.

Identification of side-chain clusters in protein structures by a graph spectral method

This paper presents a novel method to detect side-chain clusters in protein three-dimensional structures using a graph spectral approach. Protein side-chain interactions are represented by a labeled graph in which the nodes of the graph represent the Cbeta atoms and the edges represent the distance between the Cbeta atoms. The distance information and the non-bonded connectivity of the residues are represented in the form of a matrix called the Laplacian matrix. The constructed matrix is diagonalized and clustering information is obtained from the vector components associated with the second lowest eigenvalue and cluster centers are obtained from the vector components associated with the top eigenvalues. The method uses global information for clustering and a single numeric computation is required to detect clusters of interest. The approach has been adopted here to detect a variety of side-chain clusters and identify the residue which makes the largest number of interactions among the residues forming the cluster (cluster centers). Detecting such clusters and cluster centers are important from a protein structure and folding point of view. The crucial residues which are important in the folding pathway as determined by PhiF values (which is a measure of the effect of a mutation on the stability of the transition state of folding) as obtained from protein engineering methods, can be identified from the vector components corresponding to the top eigenvalues. Expanded clusters are detected near the active and binding site of the protein, supporting the nucleation condensation hypothesis for folding. The method is also shown to detect domains in protein structures and conserved side-chain clusters in topologically similar proteins.

Folding of a pressure-denatured model protein

The noncovalent complex formed by the association of two fragments of chymotrypsin inhibitor-2 is reversibly denatured by pressure in the absence of chemical denaturants. On pressure release, the complex returned to its original conformation through a biphasic reaction, with first-order rate constants of 0.012 and 0.002 s-1, respectively. The slowest phase arises from an interconversion of the pressure-denatured state, as revealed by double pressure-jump experiments. Below 5 microM, the process was concentration dependent with a second-order rate constant of 1,700 s-1 M-1. Fragment association at atmospheric pressure showed a similar break in the order of the reaction above 5 microM, but both first- and second-order folding/association rates are 2.5 times faster than those for the refolding of the pressure-denatured state. Although the folding rates of the intact protein and the association of the fragments displayed nonlinear Eyring behavior for the temperature dependence, refolding of the pressure-denatured complex showed a linear response. The negligible heat capacity of activation reflects a balance of minimal change in the burial of residues from the pressure-denatured state to the transition state. If we add the higher energy barrier in the refolding of the pressure-denatured state, the rate differences must lie in the structure of this state, which has to undergo a structural rearrangement. This clearly differs from the conformational flexibility of the isolated fragments or the largely unfolded denatured state of the intact protein in acid and provides insight into denatured states of proteins under folding conditions.

Protein folding in the absence of chemical denaturants. Reversible pressure denaturation of the noncovalent complex formed by the association of two protein fragments

Small monomeric proteins are the best models for studying protein folding, but they are often too stable for denaturation using pressure as the sole perturbant. In the present work we subject [CI-2(1-40).(41-64)], a noncovalent complex formed by the association of two complementary fragments of the chymotrypsin inhibitor-2, to high pressure to investigate the folding mechanism of a model protein. Pressures up to 3.5 kilobar do not affect the intact protein, but it can be unfolded reversibly by pressure in the presence of subdenaturing concentrations of guanidine chloride, with free energy and molar volume changes of 2.5 kcal mol-1 and 42.5 ml mol-1, respectively. In contrast, the complex can be reversibly denatured by high pressure without the addition of chemical denaturants. However, the process is clearly independent of the protein concentration, indicating lack of dissociation. We determined a change in the free energy of 1.4 kcal mol-1 and a molar volume change of 35 ml mol-1 for the pressure denaturation of the complex. A persistent quenching of the tryptophan adds further evidence for the presence of residual structure in the high pressure-denatured state. This state also appears to be compact as the small volume change indicates, compared with pressure denaturation of naturally occurring dimers. Based on observations of a number of pressure-denatured states and on characteristics of large CI-2 fragments with a solvent accessible core but maintaining tertiary interactions, the structure of the pressure-denatured state of the CI-2 complex could be explained by an ordered molten globule-like conformation.

Crystal structure of a dimeric chymotrypsin inhibitor 2 mutant containing an inserted glutamine repeat

We have constructed mutants of chymotrypsin inhibitor 2 with short glutamine repeats inserted into its inhibitory loop. These mutants oligomerize when expressed in Escherichia coli. The dimer of a mutant with four glutamines now has been crystallized, and its structure has been solved by molecular replacement by using the wild-type monomer as a search model. The structure of each half of the dimer is found to be the same as that of the wild-type monomer, except around the glutamine insertion. It was proposed that the components of the oligomers are held together by hydrogen bonds between the main-chain and side-chain amides of the glutamine repeats. Instead, they appear to form by swapping domains on folding in E. coli, and the glutamine repeats connecting the components of the dimers are disordered.

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.

Synergy between simulation and experiment in describing the energy landscape of protein folding

Experimental data from protein engineering studies and NMR spectroscopy have been used by theoreticians to develop algorithms for helix propensity and to benchmark computer simulations of folding pathways and energy landscapes. Molecular dynamic simulations of the unfolding of chymotrypsin inhibitor 2 (CI2) have provided detailed structural models of the transition state ensemble for unfolding/folding of the protein. We now have used the simulated transition state structures to design faster folding mutants of CI2. The models pinpoint a number of unfavorable local interactions at the carboxyl terminus of the single alpha-helix and in the protease-binding loop region of CI2. By removing these interactions or replacing them with stabilizing ones, we have increased the rate of folding of the protein up to 40-fold (tau = 0.4 ms). This correspondence, and other examples of agreement between experiment and theory in general, Phi-values and molecular dynamics simulations, in particular, suggest that significant progress has been made toward describing complete folding pathways at atomic resolution by combining experiment and simulation.

Calculations on folding of segment B1 of streptococcal protein G

We present an investigation of the folding thermodynamics and mechanism of segment B1 of streptococcal protein G. Molecular dynamics simulations of the fully solvated protein are used to probe thermodynamically significant states at different stages of folding. We performed several unfolding simulations to generate a database of initial conditions. The database is analyzed and clustered. The cluster centers extracted from this database were then used as starting points for umbrella sampling of the folding free energy landscape under folding conditions. The resulting sampling was combined with the weighted histogram analysis method. One and two-dimensional free energy surfaces were constructed along several order parameters and used to analyze the folding process. Our findings indicate that an initial collapse precedes the formation of significant native structure. Elements of local structure originate in the regions of the protein shown to have higher H/2H exchange protection factors in early stages of folding. A non-native contact, observed experimentally at the N terminus of the alpha-helix in a peptide excised from the protein, is seen to pre-organize the chain in early stages of folding. Collapse and early structure formation yields a compact globule with a significant number of water molecules present. Desolvation of the protein core is coincident with the final stages of folding from the compact state.

Molecular picture of folding of a small alpha/beta protein

We characterize, at the atomic level, the mechanism and thermodynamics of folding of a small alpha/beta protein. The thermodynamically significant states of segment B1 of streptococcal protein G (GB1) are probed by using the statistical mechanical methods of importance sampling and molecular dynamics. From a thermodynamic standpoint, folding commences with overall collapse, accompanied by formation of approximately 35% of the native structure. Specific contacts form at the loci experimentally inferred to be structured early in folding kinetics studies. Our study reveals that these initially structured regions are not spatially adjacent. As folding progresses, fluid-like nonlocal native contacts form, with many contacts forming and breaking as the structure searches for the native conformation. Although the alpha-helix forms early, the beta-sheet forms concomitantly with the overall topology. Water is present in the protein core up to a late stage of folding, lubricating conformational transitions during the search process. Once 80% of the native contacts have formed, water is squeezed from the protein interior and the structure descends into the native manifold. Examination of the onset of side-chain mobility within our model indicates side-chain motion is most closely linked to the overall volume of the protein and no sharp order-disorder transition appears to occur. Exploration of models for hydrogen deuterium exchange show qualitative agreement with equilibrium measurement of hydrogen/deuterium protection factors.

"New view" of protein folding reconciled with the old through multiple unfolding simulations

Twenty-four molecular dynamics trajectories of chymotrypsin inhibitor 2 provide a direct demonstration of the diversity of unfolding pathways. Comparison with experiments suggests that the transition state region for folding and unfolding occurs early with only 25 percent of the native contacts and that the root-mean-square deviations between contributing structures can be as large as 15 angstroms. Nevertheless, a statistically preferred unfolding pathway emerges from the simulations; disruption of tertiary interactions between the helix and a two-stranded portion of the beta sheet is the primary unfolding event. The results suggest a synthesis of the "new" and the classical view of protein folding with a preferred pathway on a funnel-like average energy surface.

Complementation of peptide fragments of the single domain protein chymotrypsin inhibitor 2

Chymotrypsin inhibitor 2 (CI2) folds kinetically as a single domain protein. It has been shown that elements of native secondary structure do not significantly form in fragments as the 64 residue protein is progressively increased in length from its N terminus, until at least 60 residues are present. Here, we analyse peptides of increasing length from the C terminus and find that native-like structure is not present even in the largest, fragment (7-64). We have examined sets of peptides of the form (1 - x) and ((x + 1)-64) to detect complementation. The only pair that readily complements and gives native-like structure is (1-40) and (41-64), where cleavage occurs in the protease-binding loop of CI2. But, all the pairs of peptides (1 - x) + (41-64) complement for x > 40, as do all pairs of (1-40) + (x-64), where x < 40. The resultant complexes appear to be equivalent to (1-40). (41-64) with the overlapping sequence being unstructured. Thus, the folding of CI2 is extremely co-operative, and interactions have to be made between subdomains (1-40) and (41-64). This is consistent with the mechanism proposed for the folding pathway of intact CI2 in which a diffuse nucleus is formed in the transition state between the alpha-helix in the N-terminal region of the protein and conserved hydrophobic contacts in the C-terminal region of the polypeptide. It is with these protein design features that CI2 can be an effective protease inhibitor.

Glutamine, alanine or glycine repeats inserted into the loop of a protein have minimal effects on stability and folding rates

Natural proteins can contain flexible regions in their polypeptide chain. We have investigated the effects of glycine, alanine and glutamine repeats on the stability and folding of a protein by inserting stretches of 7 to 13 residues into a suitable position in a model system, the chymotrypsin inhibitor-2 (CI2). This folds by residues (1-40) docking with residues (41-64) to form a folding nucleus. The peptides GQ4GM, GQ6GM, GQ8GM, GQ10GM, GA2SA4SA2GM and G3SG4SG3M were inserted after residue 40. The stability of the mutant proteins changes only weakly with chain length and nature of insertion, suggesting that the presence of unstructured polypeptide chains in a protein does not have a great energetic penalty. This has implications in catalysis, for example, where floppy regions have been noted in active sites, and in DNA transcription where activators, transcription factors and intermediary proteins all show long repeats of glycine/serine and/or glutamine, which are thought to be important for function. We find that the rate of folding is very insensitive to the length of the linker. The changes in rate are close to those predicted from polymer theory for the loss of configuration entropy on closing a loop. This implies that all the diffusion steps are relatively rapid.

The rate of isomerisation of peptidyl-proline bonds as a probe for interactions in the physiological denatured state of chymotrypsin inhibitor 2

There are four peptidyl-proline bonds in the 64-residue protein chymotrypsin inhibitor 2 (CI2), all of which are in the trans conformation in the native structure. The isomerisation of one or more of these peptidyl-proline bonds to the cis conformation in the denatured state gives rise to heterogeneity, leading to both fast and slow-folding species. The refolding of the fast-folding species, which has all trans peptidyl-proline bonds, is much faster than that of the slow-folding species, which have one or more cis peptidyl-proline bonds. In CI2, the slow-folding species can be classified into two groups by their rates of refolding, temperature-dependence, pH-dependence and [GdmCl]-dependence of the rate constants and the effect of peptidyl-prolyl isomerase on the rate constants. The replacement of Pro6 by Ala removes one of the slow refolding phases, suggesting that the cis peptidyl-Pro6 conformation is solely responsible for one of the slow-folding species. Pro6 is located in a region of the protein where non-random interactions have been found in a series of N-terminal fragments of CI2 (residues 1 to 13, 1 to 25, 1 to 28 and 1 to 40). In addition, NMR studies on a mutant fragment, (1-40)T3A, have confirmed that this non-native interaction is associated with the bulky side-chain of Trp5. The atypical rate of cis to trans isomerisation of the peptidyl-Pro bond is indicative of the presence of a similar hydrophobic cluster in the physiological denatured state of intact CI2.

Following co-operative formation of secondary and tertiary structure in a single protein module

We have prepared a family of peptide fragments of the 64 amino acid protein chymotrypsin inhibitor (CI2), corresponding to progressive elongation from the N terminus, in order to elucidate the basis of conformational preferences in single-domain proteins and to obtain insights into their conformational pathway. Structural analysis of the fragment comprising the first 50 residues, CI2(1-50), indicates that it is mainly disordered, with patches of hydrophobic residues exposed to the solvent. Structural characterisation of the fragment CI2(1-63) which lacks only the C-terminal glycine, Gly64, shows native-like structure in all regions of the fragment. The study provides insights into the contribution of specific residues to the stability and co-operativity of the intact protein. We define a phiNMR value, derived from chemical shift analysis, which describes the build-up of structure at the level of individual residues (protons). All the macroscopic probes used to study the growth of structure in CI2 on elongation of the chain (circular dichroism, fluorescence and gel filtration) are in agreement with the residue-by-residue description by NMR. It is seen that secondary and tertiary structure build up in parallel in the fragments and show similar structures to those developed in the transition state for folding of the intact protein.

Nucleation mechanisms in protein folding

Experiment and theory are converging on the importance of nucleation mechanisms in protein folding. These mechanisms do not use classic nuclei, which are well formed elements of structure present in ground states, but they use diffuse, extended regions, which are observed in transition states.

From Levinthal to pathways to funnels

While the classical view of protein folding kinetics relies on phenomenological models, and regards folding intermediates in a structural way, the new view emphasizes the ensemble nature of protein conformations. Although folding has sometimes been regarded as a linear sequence of events, the new view sees folding as parallel microscopic multi-pathway diffusion-like processes. While the classical view invoked pathways to solve the problem of searching for the needle in the haystack, the pathway idea was then seen as conflicting with Anfinsen's experiments showing that folding is pathway-independent (Levinthal's paradox). In contrast, the new view sees no inherent paradox because it eliminates the pathway idea: folding can funnel to a single stable state by multiple routes in conformational space. The general energy landscape picture provides a conceptual framework for understanding both two-state and multi-state folding kinetics. Better tests of these ideas will come when new experiments become available for measuring not just averages of structural observables, but also correlations among their fluctuations. At that point we hope to learn much more about the real shapes of protein folding landscapes.