Folding of chymotrypsin inhibitor 2. 2. Influence of proline isomerization on the folding kinetics and thermodynamic characterization of the transition state of folding

Folding kinetics of the all-beta-sheet protein human basic fibroblast growth factor, a structural homolog of interleukin-1beta

The refolding and unfolding kinetics of the all-beta-sheet protein human basic fibroblast growth factor (hFGF-2) were studied by fluorescence spectroscopy. The kinetics of the unfolding transition are monophasic. The refolding reaction at high and low guanidinium chloride (GdmCl) concentrations is best described by mono- and biphasic folding, respectively. Refolding and unfolding of hFGF-2 (155 amino acids) is very slow compared with other non-disulfide-bonded monomeric proteins of similar size. For example, the rate constant for unfolding at 4.5 mol.liter(-1) GdmCl is 0.006 s(-1), and the refolding rate constants at 0.4 mol.liter(-1) GdmCl are 0.01 s(-1) and 0.0009 s(-1) (15 degrees C, pH 7.0). A characterization of the thermodynamic nature of the folding process using transition state theory revealed that the slow refolding is almost exclusively controlled by entropic factors, namely the strong loss of conformational freedom during refolding. The rate of the slow unfolding kinetics is mainly (and at low denaturant concentrations exclusively) controlled by the large positive change in enthalpy. hFGF-2 shows similar slow folding kinetics to that of its structural homolog interleukin-1beta. Since both proteins show very little sequence identity, it is suggested that their slow folding kinetics are determined by the complex beta-sheet arrangement of the native molecules.

Folding pathway of FKBP12 and characterisation of the transition state

The folding pathway of human FKBP12, a 12 kDa FK506-binding protein (immunophilin), has been characterised. Unfolding and refolding rate constants have been determined over a wide range of denaturant concentrations and data are shown to fit to a two-state model of folding in which only the denatured and native states are significantly populated, even in the absence of denaturant. This simple model for folding, in which no intermediate states are significantly populated, is further supported from stopped-flow circular dichroism experiments in which no fast "burst" phases are observed. FKBP12, with 107 residues, is the largest protein to date which folds with simple two-state kinetics in water (kF=4 s(-1)at 25 degrees C). The topological crossing of two loops in FKBP12, a structural element suggested to cause kinetic traps during folding, seems to have little effect on the folding pathway. The transition state for folding has been characterised by a series of experiments on wild-type FKBP12. Information on the thermodynamic nature of, the solvent accessibility of, and secondary structure in, the transition state was obtained from experiments measuring the unfolding and refolding rate constants as a function of temperature, denaturant concentration and trifluoroethanol concentration. In addition, unfolding and refolding studies in the presence of ligand provided information on the structure of the ligand-binding pocket in the transition state. The data suggest a compact transition state relative to the unfolded state with some 70 % of the surface area buried. The ligand-binding site, which is formed mainly by two loops, is largely unstructured in the transition state. The trifluoroethanol experiments suggest that the alpha-helix may be formed in the transition state. These results are compared with results from protein engineering studies and molecular dynamics simulations (see the accompanying paper).

Development of enzymatic activity during protein folding. Detection of a spectroscopically silent native-like intermediate of muscle acylphosphatase

The recovery of enzymatic activity during the folding of muscle acylphosphatase and two single residue mutants (proline 54 to alanine and proline 71 to alanine) from 7 M urea has been monitored and compared with the development of intrinsic fluorescence emission. Fluorescence measurements reveal the presence in the wild-type protein of a major rapid refolding phase followed by a second low amplitude slow phase. The slow phase is absent in the fluorescence trace acquired with the proline 54 to alanine mutant, suggesting the involvement of this proline residue in the fluorescence-detected slow phase of the wild-type protein. The major kinetic phase is associated with a considerable recovery of enzymatic activity, indicating that a large fraction of molecules refolds with effective two-state behavior. The use of time-resolved enzymatic activity as a probe to follow the folding process reveals, however, the presence of another exponential slow phase arising from proline 71. This slow phase is not observable by utilizing optical probes, indicating that, unlike proline 54, the cis to trans isomerization of proline 71 can take place in an intermediate possessing a native-like fold. We suggest that, although spectroscopically silent and structurally insignificant, the cis-trans interconversion of proline residues in native-like intermediates may be crucial for the generation of enzymatic activity of functional enzymes.

Pseudo-prolines: induction of cis/trans-conformational interconversion by decreased transition state barriers

Molecular events such as cis/trans isomerization of Xaa-Pro tertiary amide bonds in peptides and proteins are slow on the overall time scale of the formation of a final biostructure and are, therefore, rate limiting. In order to pursue a better understanding of the molecular events underlying such slow interconversions, we applied the recently introduced pseudo-proline (PsiPro) concept as a tool to study the dynamics of Xaa-Pro bonds by determining the kinetics and thermodynamics of cis/trans isomerism. We show that enhanced isomerization rates of tertiary amide bonds prior to a PsiPro unit in short model peptides is due to lowered transition state barriers. In addition, pronounced effects upon the dynamics of the reversible transition between helix I and II of oligoprolines containing one or several PsiPro units were observed.

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.

Formation of fibrous aggregates from a non-native intermediate: the isolated P22 tailspike beta-helix domain

In the assembly pathway of the trimeric P22 tailspike protein, the protein conformation critical for the partitioning between productive folding and off-pathway aggregation is a monomeric folding intermediate. The central domain of tailspike, a large right-handed parallel beta-helix, is essentially structured in this species. We used the isolated beta-helix domain (Bhx), expressed with a hexahistidine tag, to investigate the mechanism of aggregation without the two terminal domains present in the complete protein. Although Bhx has been shown to fold reversibly at low ionic strength conditions, increased ionic strength induced aggregation with a maximum at urea concentrations corresponding to the midpoint of urea-induced folding transitions. According to size exclusion chromatography, aggregation appeared to proceed via a linear polymerization mechanism. Circular dichroism indicated a secondary structure content of the aggregates similar to that of the native state, but at the same time their tryptophan fluorescence was largely quenched. Microscopic analysis of the aggregates revealed a variety of morphologies; among others, fibrils with fine structure were observed that exhibited bright green birefringence if viewed under cross-polarized light after staining with Congo red. These observations, together with the effects of folding mutations on the aggregation process, indicate the involvement of a partially structured intermediate distinct from both unfolded and native Bhx.

The topomer-sampling model of protein folding

Clearly, a protein cannot sample all of its conformations (e.g., approximately 3(100) approximately 10(48) for a 100 residue protein) on an in vivo folding timescale (<1 s). To investigate how the conformational dynamics of a protein can accommodate subsecond folding time scales, we introduce the concept of the native topomer, which is the set of all structures similar to the native structure (obtainable from the native structure through local backbone coordinate transformations that do not disrupt the covalent bonding of the peptide backbone). We have developed a computational procedure for estimating the number of distinct topomers required to span all conformations (compact and semicompact) for a polypeptide of a given length. For 100 residues, we find approximately 3 x 10(7) distinct topomers. Based on the distance calculated between different topomers, we estimate that a 100-residue polypeptide diffusively samples one topomer every approximately 3 ns. Hence, a 100-residue protein can find its native topomer by random sampling in just approximately 100 ms. These results suggest that subsecond folding of modest-sized, single-domain proteins can be accomplished by a two-stage process of (i) topomer diffusion: random, diffusive sampling of the 3 x 10(7) distinct topomers to find the native topomer ( approximately 0.1 s), followed by (ii) intratopomer ordering: nonrandom, local conformational rearrangements within the native topomer to settle into the precise native state.

Rapid folding with and without populated intermediates in the homologous four-helix proteins Im7 and Im9

The kinetics and thermodynamics of the folding of the homologous four-helix proteins Im7 and Im9 have been characterised at pH 7.0 and 10 degrees C. These proteins are 60 % identical in sequence and have the same three-dimensional structure, yet appear to fold by different kinetic mechanisms. The logarithm of the folding and unfolding rates of Im9 change linearly as a function of urea concentration and fit well to an equation describing a two-state mechanism (with a folding rate of 1500 s-1, an unfolding rate of 0. 01 s-1, and a highly compact transition state that has approximately 95 % of the native surface area buried). By contrast, there is clear evidence for the population of an intermediate during the refolding of Im7, as indicated by a change in the urea dependence of the folding rate and the presence of a significant burst phase amplitude in the refolding kinetics. Under stabilising conditions (0.25 M Na2SO4, pH 7.0 and 10 degrees C) the folding of Im9 remains two-state, whilst under similar conditions (0.4 M Na2SO4, pH 7.0 and 10 degrees C) the intermediate populated during Im7 refolding is significantly stabilised (KUI=125). Equilibrium denaturation experiments, under the conditions used in the kinetic measurements, show that Im7 is significantly less stable than Im9 (DeltaDeltaG 9.3 kJ/mol) and the DeltaG and m values determined accord with those obtained from the fit to the kinetic data. The results show, therefore, that the population of an intermediate in the refolding of the immunity protein structure is defined by the precise amino acid sequence rather than the global stability of the protein. We discuss the possibility that the intermediate of Im7 is populated due to differences in helix propensity in Im7 and Im9 and the relevance of these data to the folding of helical proteins in general.

Stability and folding of the tumour suppressor protein p16

The tumour suppressor p16 is a member of the INK4 family of inhibi tors of the cyclin D-dependent kinases, CDK4 and CDK6, that are involved in the key growth control pathway of the eukaryotic cell cycle. The 156 amino acid residue protein is composed of four ankyrin repeats (a helix-turn-helix motif) that stack linearly as two four-helix bundles resulting in a non-globular, elongated molecule. The thermodynamic and kinetic properties of the folding of p16 are unusual. The protein has a very low free energy of unfolding, Delta GH-2O/D-N, of 3.1 kcal mol-1 at 25 degreesC. The rate-determining transition state of folding/unfolding is very compact (89% as compact as the native state). The other unusual feature is the very rapid rate of unfolding in the absence of denaturant of 0.8 s-1 at 25 degreesC. Thus, p16 has both thermodynamic and kinetic instability. These features may be essential for the regulatory function of the INK4 proteins and of other ankyrin-repeat-containing proteins that mediate a wide range of protein-protein interactions. The mechanisms of inactivation of p16 by eight cancer-associated mutations were dissected using a systematic method designed to probe the integrity of the secondary structure and the global fold. The structure and folding of p16 appear to be highly vulnerable to single point mutations, probably as a result of the protein's low stability. This vulnerability provides one explanation for the striking frequency of p16 mutations in tumours and in immortalised cell lines.

Directed evolution converts subtilisin E into a functional equivalent of thermitase

We used directed evolution to convert Bacillus subtilis subtilisin E into an enzyme functionally equivalent to its thermophilic homolog thermitase from Thermoactinomyces vulgaris. Five generations of random mutagenesis, recombination and screening created subtilisin E 5-3H5, whose half-life at 83 degrees C (3.5 min) and temperature optimum for activity (Topt, 76 degrees C) are identical with those of thermitase. The Topt of the evolved enzyme is 17 degrees C higher and its half-life at 65 degrees C is >200 times that of wild-type subtilisin E. In addition, 5-3H5 is more active towards the hydrolysis of succinyl-Ala-Ala-Pro-Phe-p-nitroanilide than wild-type at all temperatures from 10 to 90 degrees C. Thermitase differs from subtilisin E at 157 amino acid positions. However, only eight amino acid substitutions were sufficient to convert subtilisin E into an enzyme equally thermostable. The eight substitutions, which include known stabilizing mutations (N218S, N76D) and also several not previously reported, are distributed over the surface of the enzyme. Only two (N218S, N181D) are found in thermitase. Directed evolution provides a powerful tool to unveil mechanisms of thermal adaptation and is an effective and efficient approach to increasing thermostability without compromising enzyme activity.

Stability and folding of the cell cycle regulatory protein, p13(suc1)

p13(suc1) (suc1) is a member of the CDC28 kinase specific family of cell cycle regulatory proteins that bind to the cyclin-dependent kinase CDK2 and regulate its activity. suc1 has two distinct conformational and assembly states, a compact globular monomer and a beta strand-exchanged dimer. The dimerisation is an example of domain-swapping, and is mediated by a molecular hinge mechanism that is conserved across the entire CKS family. It has been proposed that the function of suc1 may be modulated by the dimerisation process with monomer-dimer switching occurring in response to a change in the cell environment. We have investigated the stability and folding of suc1 as a first step in determining the mechanism and functional role of the strand exchange. Suc1 unfolds reversibly at equilibrium in a two-state manner with a free energy of unfolding of 7.2 kcal mol-1. The kinetics of folding and unfolding are complex, and double-jump stopped-flow methods revealed that there are at least three parallel folding pathways arising from distinct unfolded and partly folded, intermediate states. The major population of unfolded species fold rapidly according to a three-state mechanism, D1->I1->N, with a rate constant for the formation of native species, N, from the intermediate, I1, of 65 s-1 in water. Two minor populations of unfolded molecules fold more slowly. Folding of one population is limited by proline isomerisation in a partly folded state, and some expansion of the protein is required for isomerisation to occur. The other population could be assigned to rate-limiting isomerisation of the peptidyl-proline bond of residue 90, which is located in the molecular hinge. A minor, fast phase was detected in the unfolding kinetics that corresponds to unfolding of a small population of a distinct native-like form. Heterogeneity was removed upon mutation of Pro90 to Ala. The unfolding kinetics of the strand-exchanged dimer were also investigated and showed that the dimer unfolds at the same rate as the monomer.

pH dependence of the activation parameters for chymopapain unfolding: influence of ion pairs on the kinetic stability of proteins

We studied the irreversible thermal denaturation of chymopapain, a papain-related cysteine proteinase. It was found that this process follows simple first-order kinetics under all conditions tested. Rate constants determined by monitoring ellipticity changes at 220 or 279 nm are essentially identical, indicating that denaturation involves global unfolding of the protein. Enthalpies (DeltaH(double dagger)) and entropies (DeltaS(double dagger)) of activation for unfolding were determined at various pH values from the temperature dependence of the rate constant. In the pH range 1.1-3.0, a large variation of both DeltaH(double dagger) and DeltaS(double dagger) was observed. For the few proteins studied so far (lysozyme, trypsin, barnase) it is known that activation parameters for unfolding vary little with pH. It is proposed that this contrasting behavior of chymopapain originates from the numerous ion pairs - especially those with low solvent accessibilities - present in its molecular structure. In contrast, fewer, more exposed ion pairs are present in the other proteins mentioned above. Our results were analyzed in terms of differences in the protonation behavior of carboxylic groups between the transition (TS) and native (N) states of the protein. For this purpose, a model of independently titrating sites was assumed, which explained reasonably well the pH dependence of activation parameters, as well as the protonation properties of native chymopapain. According to these calculations, pK values of carboxyls in TS are shifted 0.6-0.9 units upwards with respect to those in N. In addition, some groups in TS appear to be protonated with unusually large enthalpy changes.

A pulse-chase-competition experiment to determine if a folding intermediate is on or off-pathway: application to ribonuclease A

A modified pulse-chase experiment is applied to determine if the native-like intermediate IN of ribonuclease A is on or off-pathway. The 1H label retained in the native protein is compared when separate samples of 1H-labeled IN and unfolded protein are allowed to fold to native in identical conditions. The solvent is 2H2O and the pH* is such that the unfolded protein rapidly exchanges its peptide NH protons with solvent, and IN does not. If IN is on-pathway, more 1H-label will be retained in the test sample starting with IN than in the control sample starting with unfolded protein. The results show that IN is a productive (on-pathway) intermediate. Application of the modified pulse-chase experiment to the study of rapidly formed folding intermediates may be possible when a rapid mixing device is used.

Folding under inequilibrium conditions as a possible reason for partial irreversibility of heat-denatured proteins: computer simulation study

Using computer simulations we have studied possible effects of heating and cooling at different scan rates on unfolding and refolding of macromolecules. We have shown that even the simplest two-state reversible transition can behave irreversibly when an unfavorable combination of cooling rate, relaxation time and activation energy of refolding occurs. On the basis of this finding we suppose that apparent irreversibility of some proteins denatured by heat may result from slow relaxation on cooling rather than thermodynamic instability and/or irreversible alterations of the polypeptide chain. Using this kinetic reversible two-state model, we estimated the effects of the scan rate and kinetic parameters of the macromolecule on its unfolding-refolding process. A few recommendations are suggested on how to reach maximal possible recovery after denaturation if refolding appears to be under kinetic control.

Side-chain effects on peptidyl-prolyl cis/trans isomerisation

Peptidyl-prolyl cis/trans isomerisation has been frequently found as a rate limiting step in the folding of proteins. In order to determine whether the nature of the amino acid preceding proline controls the probability of cis prolyl bonds in native proteins, systematic studies on the thermodynamics and kinetics of the prolyl isomerisation in the pentapeptide series Ac-Ala-Xaa-Pro-Ala-Lys-NH2 were performed. All proteinogenic amino acids were substituted in the position preceding proline. When measured by 1H-NMR and CD spectroscopy both isomers proved to be devoid of ordered structure in the whole series of the oligopeptides in aqueous solution. Thus, isomerization rates and cis/trans ratios calculated from solvent jump and 1H-NMR magnetisation transfer experiments exclusively reflect the side-chain effects of the Xaa position in the peptide series. There is a rough correlation between the cis content in the oligopeptides and the propensity of Xaa-Pro cis prolyl bonds in proteins. This correlation suggests that the prolyl bond conformation is mainly determined by local effects in proteins. The rate constants kc-->t of pentapeptides containing unionised amino acids preceding proline range from 3.2 x 10(-3) s-1 (Xaa = Ala) to 0.5 x 10(-3) s-1 (Xaa = Trp) at 4 degrees C. Proline clustering led to an isomerisation cycle indicating considerable influence on the isomerisation rates of the peptide bond conformations flanking the rotating bond. Both tyrosine and histidine specifically reduce isomerisation rates severalfold by deprotonation of their respective side-chains.

Kinetic role of electrostatic interactions in the unfolding of hyperthermophilic and mesophilic rubredoxins

The temperature dependence of the unfolding kinetics of rubredoxins from the hyperthermophile Pyrococcus furiosus (RdPf) and the mesophile Clostridium pasteurianum (RdCp) has been studied. Results show that RdPf unfolds much more slowly, under all experimentally accessible temperature regimes, than RdCp and other typical mesophilic proteins. Rates of RdCp and RdPf unfolding decrease upon increasing the pH above 2 and diverge dramatically at pH 7. As shown by detailed electrostatic energy calculations, this is the result of a differential degree of protonation of the negatively charged amino acids, which causes distinct electrostatic configurations as a function of pH. We propose that ion pairs, particularly those that are placed in key surface positions, may play a kinetic role by mildly clamping the protein and thereby influencing the nature and the number of the vibrational normal modes that are thermally accessible upon unfolding. More generally, these modes are also likely to be affected by the favorable electrostatic configurations, which we have shown to be directly linked to the extremely slow unfolding rates of RdPf at neutral pH. Even at pH 2, in the absence of any salt bridges, the unfolding rates of RdPf are much smaller than those of RdCp. This is ascribed to presently unidentified structural elements of nonelectrostatic nature. Since electrostatic effects influence the unfolding kinetics of both mesophilic and thermophilic rubredoxins, these findings may be of general significance for proteins.

Folding kinetics of the SH3 domain of PI3 kinase by real-time NMR combined with optical spectroscopy

The refolding kinetics of the chemically denatured SH3 domain of phosphatidylinositol 3'-kinase (PI3-SH3) have been monitored by real-time one-dimensional 1H NMR coupled with a variety of other biophysical techniques. These experiments indicate that the refolding kinetics of PI3-SH3 are biphasic. The slow phase (27 (+/- 8)% amplitude) is due to a population of substantially unfolded molecules with an incorrectly configured cis proline residue. The fast phase (73 (+/- 8)% amplitude) corresponds to the folding of protein molecules with proline residues in a trans configuration in the unfolded state. NMR experiments indicate that the first species populated after the initiation of folding exhibit poor chemical shift dispersion and have spectra very similar to that of the denatured protein in 8 M guanidine hydrochloride. Linear combinations of the first spectrum and of the spectrum of the native protein accurately reconstruct all of the spectra acquired during refolding. Consistent with this, native side-chain and backbone H alpha atom packing (NMR), secondary structure (far-UV circular dichroism), burial of aromatic residues (near-UV circular dichroism), intrinsic fluorescence and peptide binding activity are all recovered with effectively identical kinetics. Equilibrium unfolding and folding/unfolding kinetics yield, within experimental error, identical values for the free energy of unfolding (delta Gu-H2O = 3.38 kcal mol-1) and for the slope of the free energy of unfolding versus denaturant concentration (meq = 2.33 kcal mol-1 M-1). Together, these data provide strong evidence that PI3-SH3 folds without significant population of kinetic well-structured intermediates. That PI3-SH3 folds slowly (time constant 2.8 seconds in H2O at 20 degrees C) indicates that slow refolding is not always a consequence of kinetic traps but may be observed even when a protein appears to fold via a simple, two-state mechanism.

A folding pathway for betapep-4 peptide 33mer: from unfolded monomers and beta-sheet sandwich dimers to well-structured tetramers

It was recently reported that a de novo designed peptide 33mer, betapep-4, can form well-structured beta-sheet sandwich tetramers (Ilyina E, Roongta V, Mayo KH, 1997b, Biochemistry 36:5245-5250). For insight into the pathway of betapep-4 folding, the present study investigates the concentration dependence of betapep-4 self-association by using 1H-NMR pulsed-field gradient (PFG)-NMR diffusion measurements, and circular dichroism. Downfield chemically shifted alphaH resonances, found to arise only from the well-structured betapep-4 tetramer state, yield the fraction of tetramer within the oligomer equilibrium distribution. PFG-NMR-derived diffusion coefficients, D, provide a means for deriving the contribution of monomer and other oligomer states to this distribution. These data indicate that tetramer is the highest oligomer state formed, and that inclusion of monomer and dimer states in the oligomer distribution is sufficient to explain the concentration dependence of D values for betapep-4. Equilibrium constants calculated from these distributions [2.5 x 10(5) M(-1) for M-D and 1.2 x 10(4) M(-1) for D-T at 313 K] decrease only slightly, if at all, with decreasing temperature indicating a hydrophobically mediated, entropy-driven association/folding process. Conformational analyses using NMR and CD provide a picture where "random coil" monomers associate to form molten globule-like beta-sheet sandwich dimers that further associate and fold as well-structured tetramers. Betapep-4 folding is thermodynamically linked to self-association. As with folding of single-chain polypeptides, the final folding step to well-structured tetramer betapep-4 is rate limiting.

Protein folding in the landscape perspective: chevron plots and non-Arrhenius kinetics

We use two simple models and the energy landscape perspective to study protein folding kinetics. A major challenge has been to use the landscape perspective to interpret experimental data, which requires ensemble averaging over the microscopic trajectories usually observed in such models. Here, because of the simplicity of the model, this can be achieved. The kinetics of protein folding falls into two classes: multiple-exponential and two-state (single-exponential) kinetics. Experiments show that two-state relaxation times have "chevron plot" dependences on denaturant and non-Arrhenius dependences on temperature. We find that HP and HP+ models can account for these behaviors. The HP model often gives bumpy landscapes with many kinetic traps and multiple-exponential behavior, whereas the HP+ model gives more smooth funnels and two-state behavior. Multiple-exponential kinetics often involves fast collapse into kinetic traps and slower barrier climbing out of the traps. Two-state kinetics often involves entropic barriers where conformational searching limits the folding speed. Transition states and activation barriers need not define a single conformation; they can involve a broad ensemble of the conformations searched on the way to the native state. We find that unfolding is not always a direct reversal of the folding process.

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.

Molecular recognition in the HIV-1 capsid/cyclophilin A complex

The HIV-1 capsid protein (CA) makes an essential interaction with the human peptidyl prolyl isomerase, cyclophilin A (CypA), that results in packaging of CypA into the virion at a CA to CypA stoichiometry of approximately 10:1. The 231 amino acid residue capsid protein is composed of an amino-terminal CypA binding domain (1 to approximately 151; CA151) and a carboxyl-terminal dimerization domain (approximately 151 to 231). We find that CypA binds dimeric CA and monomeric CA151 with identical intrinsic affinities (K[d] = 16(+/-4) microM). This result demonstrates that capsid dimerization and cyclophilin A binding are not thermodynamically coupled and suggests that the substoichiometric ratio of CypA in the HIV-1 virion results from the intrinsic stability of the CA/CypA complex. In the known co-crystal structure of the CA151/CypA complex, CypA binding is mediated exclusively by an exposed capsid loop that spans residues Pro85 to Pro93. The energetic contributions to CypA binding were quantified for each residue in this loop, and the results demonstrate that the Gly89-Pro90 dipeptide is the primary cyclophilin A recognition motif, with Pro85, Val86, His87, Ala88, and Pro93 also making energetically favorable contacts. These studies reveal that the active site of CypA, which can catalyze the isomerization of proline residues in vitro, also functions as a sequence-specific, protein-binding motif in HIV-1 replication.

Diffusion control in an elementary protein folding reaction

The cold-shock protein CspB (from Bacillus subtilis), a very small protein of 67 residues, folds extremely fast in a reversible N &lrharr; U two-state reaction. Both unfolding and refolding are strongly decelerated when the viscosity of the solvent is increased by adding ethylene glycol or sucrose. The folding of CspB thus seems to follow Kramers' model for reactions in which the reactants must diffuse together. It indicates that the compaction of the protein chain occurs in the rate-limiting step of folding. Chain diffusion to a productively collapsed form and the crossing of a high energy barrier are thus tightly coupled in this folding reaction, and the measured reaction rate depends on both the diffusion of the protein chain in the solvent and the magnitude of the activation energy. We suggest that in protein folding an energetic barrier is essential to separate the native from the unfolded conformations of a protein. This barrier protects the ordered structure of a native protein against continuous unfolding by diffusive chain motions and leads to apparent two-state behavior.

Folding of the disulfide-bonded beta-sheet protein tendamistat: rapid two-state folding without hydrophobic collapse

We investigated the reversible folding and unfolding reactions of the small 74 amino acid residue protein tendamistat. The secondary structure of tendamistat contains only beta-sheets and loop regions and the protein contains two disulfide bonds. Fluorescence-detected refolding kinetics of tendamistat (disulfide bonds intact) comprise of a major rapid fast reaction (tau = 10 ms in water) and two minor slow reactions. In the fast reaction 80% of the unfolded molecules are converted to native protein. The two slow reactions are part of a parallel slow folding pathway. On this pathway the rate-limiting step in the formation of native molecules is cis to trans isomerization of at least one of the three trans Xaa-Pro peptide bonds. This reaction is catalyzed efficiently by the enzyme peptidyl-prolyl cis-trans isomerase. Comparison of kinetic data with equilibrium unfolding transitions shows that the fast folding pathway follows a two-state process without populated intermediate states. Additionally, various sensitive tests did not detect any rapid chain collapse during tendamistat folding prior to the acquisition of the native three-dimensional structure. These results show that pre-formed disulfide bonds do not prevent efficient and rapid protein folding.

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.

Extreme stabilization of a thermolysin-like protease by an engineered disulfide bond

The thermal inactivation of broad specificity proteases such as thermolysin and subtilisin is initiated by partial unfolding processes that render the enzyme susceptible to autolysis. Previous studies have revealed that a surface-located region in the N-terminal domain of the thermolysin-like protease produced by Bacillus stearothermophilus is crucial for thermal stability. In this region a disulfide bridge between residues 8 and 60 was designed by molecular modelling, and the corresponding single and double cysteine mutants were constructed. The disulfide bridge was spontaneously formed in vivo and resulted in a drastic stabilization of the enzyme. This stabilization presents one of the very few examples of successful stabilization of a broad specificity protease by a designed disulfide bond. We propose that the success of the present stabilization strategy is the result of the localization and mutation of an area of the molecule involved in the partial unfolding processes that determine thermal stability.

The folding pathway of a protein at high resolution from microseconds to seconds

We have documented the folding pathway of the 10-kDa protein barstar from the first few microseconds at the resolution of individual residues from its well characterized denatured state. The denatured state had been shown from NMR to have flickering native-like structure in the first two of its four alpha-helices. phi-value analysis shows that the first helix becomes substantially consolidated as the intermediate is formed in a few hundred microseconds, as does the second to a lesser extent. A native-like structure then is formed in a few hundred milliseconds as the whole structure consolidates. Peptide fragments corresponding to sequences containing the first two helices separately and together as a helix-loop-helix motif have little helical structure under conditions that favor folding. The early stages of folding fit the nucleation-condensation model that was proposed for the smaller chymotrypsin inhibitor 2, which is a single module of structure and folds by two-state kinetics. The early stages of the multistate folding of the larger, multimodular, barnase have proved experimentally inaccessible. The folding pathway of barstar links those of CI2 and barnase to give a unified scheme for folding.

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.

Cyclophilin A complexed with a fragment of HIV-1 gag protein: insights into HIV-1 infectious activity

BACKGROUND

Cyclophilin A (CyPA), a receptor of the immunosuppressive drug cyclosporin A, catalyzes the cis-trans isomerization of peptidyl-prolyl bonds and is required for the infectious activity of human immunodeficiency virus type 1 (HIV-1). The crystal structure of CyPA complexed with a fragment of the HIV-1 gag protein should provide insights into the nature of CyPA-gag interactions and may suggest a role for CyPA in HIV-1 infectious activity.

RESULTS

The crystal structure of CyPA complexed with a 25 amino acid peptide of HIV-1 gag capsid protein (25-mer) was determined and refined to an R factor of 0.195 at 1.8 A resolution. The sequence Ala88-Gly89-Pro90-Ile91 of the gag fragment is the major portion to bind to the active site of CyPA. Two residues of the 25-mer (Pro90-Ile91) bind to CyPA in a similar manner to two residues (Pro-Phe) of the CyPA substrate, succinyl-Ala-Ala-Pro-Phe-p-nitroanilide (AAPF). However, the N-terminus of the 25-mer (Ala88-Gly89) exhibits a different hydrogen-bonding pattern and molecular conformation than AAPF. The peptidyl-prolyl bond between Gly89 and Pro90 of the 25-mer has a trans conformation, in contrast to the cis conformation observed in other known CyPA-peptide complexes. The residue preceding proline, Gly89, has an unfavorable backbone conformation usually only adopted by glycine.

CONCLUSIONS

The unfavorable backbone conformation of Gly89 of the gag 25-mer fragment suggests that binding between HIV-1 gag protein and CyPA requires a special sequence, Gly-Pro. Thus, in HIV-1 infectivity, CyPA is likely to function as a chaperone, rather than as a cis-trans isomerase. However, the observation of similarities between the C termini of the 25-mer and the substrate AAPF means that the involvement of the cis-trans isomerase activity of CyPA cannot be completely ruled out.

Protein folding kinetics: timescales, pathways and energy landscapes in terms of sequence-dependent properties

BACKGROUND

Recent experimental and theoretical studies have revealed that protein folding kinetics can be quite complex and diverse depending on various factors such as size of the protein sequence and external conditions. For example, some proteins fold apparently in a kinetically two-state manner, whereas others follow complex routes to the native state. We have set out to provide a theoretical basis for understanding the diverse behavior seen in the refolding kinetics of proteins in terms of properties that are intrinsic to the sequence.

RESULTS

The folding kinetics of a number of sequences for off-lattice continuum models of proteins is studied using Langevin simulations at two different values of the friction coefficient. We show for these models that there is a remarkable correlation between folding time, tau F, and sigma = (T theta - TF)/T theta, where T theta and TF are the equilibrium collapse and folding transition temperatures, respectively. The microscopic dynamics reveals that several scenarios for the kinetics of refolding arise depending on the range of values of sigma. For relatively small sigma, the chain reaches the native conformation by a direct native conformation nucleation collapse (NCNC) mechanism without being trapped in any detectable intermediates. For moderate and large values of sigma, the kinetics is described by the kinetic partitioning mechanism, according to which a fraction of molecules phi (kinetic partition factor) reach the native conformation via the NCNC mechanism. The remaining fraction attains the native state by off-pathway processes that involve trapping in several misfolded structures. The rate-determining step in the off-pathway processes is the transition from the misfolded structures to the native state. The partition factor phi is also determined by sigma: the smaller the value of sigma, the larger is phi. The qualitative aspects of our results are found to be independent of the friction coefficient. The simulation results and theoretical arguments are used to obtain estimates for timescales for folding via the NCNC mechanism in small proteins, those with less than about 70 amino acid residues.

CONCLUSIONS

We have shown that the various scenarios for folding of proteins, and possibly other biomolecules, can be classified solely in terms of sigma. Proteins with small values of sigma reach the native conformation via a nucleation collapse mechanism and their energy landscape is characterized by having one dominant native basin of attraction (NBA). On the other hand, proteins with large sigma get trapped in competing basins of attraction (CBAs) in which they adopt misfolded structures. Only a small fraction of molecules access the native state rapidly when sigma is large. For these sequences, the majority of the molecules approach the native state by a three-stage multipathway mechanism in which the rate-determining step involves a transition from one of the CBAs to the NBA.

Unfolding kinetics of bovine trypsinogen

The unfolding kinetics of bovine trypsinogen were studied by a fluorescence-detected stopped-flow technique at pH 5.8. Trypsinogen unfolding appeared to be a rather complex reaction. Two phases, fast (with a time constant in the millisecond range) and slow, were detected in the range 2-7 M guanidium chloride (GdmCl). The natural logarithm of the rate constant of the slow phase exhibited strong dependence on [GdmCl], changing from hundreds of seconds at low denaturant concentration to about 20 ms at 7 M GdmCl. The curvature of this dependence further suggests a complex mechanism of unfolding. Generally, similar kinetics were observed for the trypsinogen.Ca complex. Small differences could be noticed, however, for the fast phase. In agreement, Ca2+ influenced only this stage of the reaction. Analysis of the dependence of the time constant of the fast phase on [CaCl2] indicates that at 4 M GdmCl, trypsinogen.Ca unfolds about sixfold slower than free zymogen, and that native trypsinogen at 4 M GdmCl still exhibits high affinity for Ca2+. Limited data on trypsin unfolding show virtually an identical dependence of the slow phase on [GdmCl]; the fast phase, however was not observed. Moreover, in the 3-4.5 M GdmCl range, a separate phase was detected. It is postulated that this phase is a manifestation of the activation-domain unfolding. The Eyring plots for the fast phase of . trypsinogen and trypsinogen.Ca unfolding are linear, indicating little change in heat capacity for this stage of reaction. The slow step of unfolding, however, shows significant curvature which indicates a substantial increase in heat capacity.

The dimeric DNA binding domain of the human papillomavirus E2 protein folds through a monomeric intermediate which cannot be native-like

The dimeric DNA binding domain of the human papillomavirus E2 protein displays a two-state concerted unfolding and dissociation, with no detectable monomeric intermediate species accumulated at equilibrium. We investigated the kinetic folding mechanism of the dimeric domain using stopped-flow spectroscopic techniques and observed a fast forming monomeric intermediate, followed by a slower bimolecular reaction. Both phases involve secondary structure rearrangements of similar magnitude. Our results support a folding pathway in which the formation of an early monomeric intermediate, with characteristics of hydrophobic collapse, is followed by a bimolecular step encompassing association and folding. The interwoven folding topology of this particular type of dimeric beta-barrel found in the E2 DNA binding domain strongly suggests that any monomeric species formed could not be native-like.

Protein folding thermodynamics applied to the photocycle of the photoactive yellow protein

Two complementary aspects of the thermodynamics of the photoactive yellow protein (PYP), a new type of photoreceptor that has been isolated from Ectothiorhodospira halophila, have been investigated. First, the thermal denaturation of PYP at pH 3.4 has been examined by global analysis of the temperature-induced changes in the UV-VIS absorbance spectrum of this chromophoric protein. Subsequently, a thermodynamic model for protein (un)folding processes, incorporating heat capacity changes, has been applied to these data. The second aspect of PYP that has been studied is the temperature dependence of its photocycle kinetics, which have been reported to display an unexplained deviation from normal Arrhenius behavior. We have extended these measurements in two solvents with different hydrophobicities and have analyzed the number of rate constants needed to describe these data. Here we show that the resulting temperature dependence of the rate constants can be quantitatively explained by the application of a thermodynamic model which assumes that heat capacity changes are associated with the two transitions in the photocycle of PYP. This result is the first example of an enzyme catalytic cycle being described by a thermodynamic model including heat capacity changes. It is proposed that a strong link exists between the processes occurring during the photocycle of PYP and protein (un)folding processes. This permits a thermodynamic analysis of the light-induced, physiologically relevant, conformational changes occurring in this photoreceptor protein.

Alpha-helical, but not beta-sheet, propensity of proline is determined by peptide environment

Proline is established as a potent breaker of both alpha-helical and beta-sheet structures in soluble (globular) proteins. Thus, the frequent occurrence of the Pro residue in the putative transmembrane helices of integral membrane proteins, particularly transport proteins, presents a structural dilemma. We propose that this phenomenon results from the fact that the structural propensity of a given amino acid may be altered to conform to changes imposed by molecular environment. To test this hypothesis on proline, we synthesized model peptides of generic sequence H2N-(Ser-LyS)2-Ala- Leu-Z-Ala-Leu-Z-Trp-Ala-Leu-Z-(Lys-Ser)3-OH (Z = Ala and/or Pro). Peptide conformations were analyzed by circular dichroism spectroscopy in aqueous buffer, SDS, lysophosphatidylglycerol micelles, and organic solvents (methanol, trifluoroethanol, and 2-propanol). The helical propensity of Pro was found to be greatly enhanced in the membrane-mimetic environments of both lipid micelles and organic solvents. Proline was found to stabilize the alpha-helical conformation relative to Ala at elevated temperatures in 2-propanol, an observation that argues against the doctrine that Pro is the most potent alpha-helix breaker as established in aqueous media. Parallel studies in deoxycholate micelles of the temperature-induced conformational transitions of the single-spanning membrane bacteriophage IKe major coat protein, in which the Pro-containing wild type was compared with Pro30 --> Ala mutant, Pro was found to protect the helix, but disrupt the beta-sheet structure as effectively as it does to model peptides in water. The intrinsic capacity of Pro to disrupt beta-sheets was further reflected in a survey of porins where Pro was found to be selectively excluded from the core of membrane-spanning beta-sheet barrels. The overall data provide a rationale for predicting and understanding the structural consequences when Pro occurs in the context of a membrane.

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

BACKGROUND

Single-module proteins, such as chymotrypsin inhibitor 2 (CI2), fold as a single cooperative unit. To solve its folding pathway, we must characterize, under conditions that favour folding, its denatured state, its transition state, and its final folded structure. To obtain a "denatured state' that can readily be thus characterized, we have used a trick of cleaving CI2 into two complementary fragments that associate and fold in a similar way to intact protein.

RESULTS

Fragment CI2(1-40)-which contains the sequence of the single alpha-helix, spanning residues 12-24-and CI2(41-64), and mutants thereof, were analyzed by NMR spectroscopy, the transition state for association/folding was characterized by the protein engineering method, and the structure of the complex was solved by NMR and X-ray crystallography. Both isolated fragments are largely disordered. The transition state for association/folding is structured around a nucleus of a nearly fully formed alpha-helix, as is the transition state for the folding of intact CI2, from residues Ser12 to Leu21, Ala16, a residue from the helix whose sidechain is buried in the hydrophobic core, makes interactions with Leu49 and Ile57 in the other fragment. Ala16 makes its full interaction energy in the transition state for the association/folding reaction, just as found during the folding of the intact protein.

CONCLUSIONS

The specific contacts in the transition state from a nucleus that extends from one fragment to the next, but the nucleus is only "flickeringly' present in the denatured state. This is direct evidence for the nucleation-condensation mechanism in which the nucleus is only weakly formed in the ground state and develops in the transition state. The low conformational preferences in the denatured state are not enough to induce significant local secondary structure, but are reinforced by tertiary interactions during the rapid condensation around the nucleus.

On-pathway versus off-pathway folding intermediates

Rapidly formed molten globule intermediates accumulate at the start of the folding reactions of several small proteins. Opinion is sharply divided as to whether they are on-pathway or off-pathway intermediates. I discuss recent experiments aimed at resolving this issue. Specific points include whether a 'rollover' in the plot of folding rate versus denaturant concentration implies that a folding intermediate is or is not on-pathway; whether the failure to observe folding intermediates for some small proteins implies a different folding mechanism or only that the intermediates are less stable; possible interpretation of 'fast-track' folding of hen lysozyme; and the significance of recent results in the search for unfolding intermediates.

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

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

Submillisecond folding of monomeric lambda repressor

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.

Structural determinants of the stability of thermolysin-like proteinases

Thermolysin is a member of a family of homologous proteinases which differ in their resistance to thermally induced unfolding and subsequent autolytic degradation. Site-directed mutagenesis studies of the thermolysin-like proteinase (TLP) from Bacillus stearothermophilus (TLP-ste) show that its reduced resistance to thermally induced autolysis, as compared to thermolysin, is due to only some of the 44 naturally occurring amino-acid differences between them. In fact TLP-ste becomes more resistant than thermolysin by mutation of just a few of these amino-acids. The crucial differences are all localized to a solvent-exposed region in the N-terminal domain of TLP-ste.

The thermal denaturation of stem bromelain is consistent with an irreversible two-state model

The thermal denaturation of bromelain, a cysteine proteinase from the papain family, was studied by means of circular dichroism (CD) and differential scanning calorimetry (DSC). It was found that this process is completely irreversible and apparently follows a simple two-state mechanism of the type N-->D. The activation energy, E, that characterizes this reaction was calculated by the use of different approaches: (i) the effect of heating rate on the temperature at which the transition is half completed; (ii) analysis of individual transition curves; (iii) kinetic studies at fixed temperatures; and (iv) single DSC tracings. The obtained values for E were rather similar to one another, varying from 164 to 226 kJ/mol. In comparison, the total calorimetric enthalpy change was 334 kJ/mol. When a more complex mechanism is considered (N<-->U-->D), which takes into account the presence of a reversibly unfolded state (U), our results suggest that the rate-limiting step is precisely the formation of U. Calculation of the corresponding activation enthalpy and entropy also seems to support this proposal.

Folding of a nascent polypeptide chain in vitro: cooperative formation of structure in a protein module

We have prepared a family of peptide fragments of the 64-residue chymotrypsin inhibitor 2, corresponding to its progressive elongation from the N terminus. The growing polypeptide chain has little tendency to form stable structure until it is largely synthesized, and what structures are formed are nonnative and lack, in particular, the native secondary structural elements of alpha-helix and beta-sheet. These elements then develop as sufficient tertiary interactions are made in the nearly full-length chain. The growth of structure in the small module is highly cooperative and does not result from the hierarchical accretion of substructures.

Funnels, pathways, and the energy landscape of protein folding: a synthesis

The understanding, and even the description of protein folding is impeded by the complexity of the process. Much of this complexity can be described and understood by taking a statistical approach to the energetics of protein conformation, that is, to the energy landscape. The statistical energy landscape approach explains when and why unique behaviors, such as specific folding pathways, occur in some proteins and more generally explains the distinction between folding processes common to all sequences and those peculiar to individual sequences. This approach also gives new, quantitative insights into the interpretation of experiments and simulations of protein folding thermodynamics and kinetics. Specifically, the picture provides simple explanations for folding as a two-state first-order phase transition, for the origin of metastable collapsed unfolded states and for the curved Arrhenius plots observed in both laboratory experiments and discrete lattice simulations. The relation of these quantitative ideas to folding pathways, to uniexponential vs. multiexponential behavior in protein folding experiments and to the effect of mutations on folding is also discussed. The success of energy landscape ideas in protein structure prediction is also described. The use of the energy landscape approach for analyzing data is illustrated with a quantitative analysis of some recent simulations, and a qualitative analysis of experiments on the folding of three proteins. The work unifies several previously proposed ideas concerning the mechanism protein folding and delimits the regions of validity of these ideas under different thermodynamic conditions.

Theoretical studies of protein folding and unfolding

The mechanism of protein folding is being investigated theoretically by the use of both simplified and all-atom models of the polypeptide chain. Lattice heteropolymer simulations of the folding process have led to proposals for the folding mechanism and for the resolution of the Levinthal paradox. Both stability and rapid folding have been shown in model studies to result from the presence of a pronounced global energy minimum corresponding to the native state. Concomitantly, molecular dynamics simulations with detailed atomic models have been used to analyze the initial stages of protein unfolding. Results concerning possible folding intermediates and the role of water in the unfolding process have been obtained. The two types of theoretical approaches are providing information essential for an understanding of the mechanism of protein folding and are useful for the design of experiments to study the mechanism in different proteins.

The structure of the transition state for the association of two fragments of the barley chymotrypsin inhibitor 2 to generate native-like protein: implications for mechanisms of protein folding

Possible early events in protein folding may be studied by dissecting proteins into complementary fragments. Two fragments of chymotrypsin inhibitor 2 [CI2-(20-59) and CI2-(60-83)] associate to form a native-like structure in a second-order reaction that combines collision and rearrangement. The transition state of the reaction, analyzed by the protein engineering method on 17 mutants, is remarkably similar to that for the folding of the intact protein--a structure that resembles an expanded version of the folded structure with most interactions significantly weakened. The exception is that the N-terminal region of the single alpha-helix (the N-capping box) is completely formed in the transition state for association of the fragments, whereas it is reasonably well formed for the intact protein. Preliminary evidence on the structures of the individual fragments indicates that both are mainly nonnative, lacking native secondary structure and having regions of nonnative buried hydrophobic clusters. The association reaction does not result from the collision of a subpopulation of two fully native-like fragments but involves a considerable rearrangement of structure.

Characterization of the transition state of protein unfolding by use of molecular dynamics: chymotrypsin inhibitor 2

Temperature-induced unfolding of chymotrypsin inhibitor 2 in water was investigated by molecular dynamics simulations. The major transition state of unfolding was identified on the basis of structural and conformational changes in the protein during the unfolding reaction. The native tertiary contacts in the hydrophobic core were considerably disrupted in the transition state, whereas the secondary structure was partially intact. The extent of structural change of the protein around a particular residue was represented quantitatively by the ratio of the number of contacts the residue makes in the transition state relative to the native state, phi MD, which allows quantitative comparison with the experimentally determined phi F values. For the region of the unfolding trajectory that is identified as the transition state, the phi MD and phi F values are in good agreement, suggesting that the transition state identified in the unfolding simulation corresponds to that probed with protein engineering methods. Although speculative, the transition state identified in the simulation is consistent with available experimental data and provides an in-depth view of what the transition state of unfolding may look like.

Structure of the transition state for the folding/unfolding of the barley chymotrypsin inhibitor 2 and its implications for mechanisms of protein folding

The equilibrium and kinetics of folding of the single-domain protein chymotrypsin inhibitor 2 conform to the simple two-state model. The structure of the rate-determining transition state has been mapped out at the resolution of individual side chains by using the protein engineering method on 74 mutants that have been constructed at 37 of the 64 residues. The structure contains no elements of secondary structure that are fully formed. The majority of interactions are weakened by > 50% in the transition state, although most regions do have some very weak structure. The structure of the transition state appears to be an expanded form of the native state in which secondary and tertiary elements have been partly formed concurrently. This is consistent with a "global collapse" model of folding rather than a framework model in which folding is initiated from fully preformed local secondary structural elements. This may be a general feature for the folding of proteins lacking a folding intermediate and is perhaps representative of the early stages of folding for multidomain or multimodule proteins. The major transition state for the folding of barnase, for example, has some fully formed secondary and tertiary structural elements in the major transition state, and barnase appears to form by a framework process. However, the fully formed framework may be preceded by a global collapse, and a unified folding scheme is presented.