Probing the structure of folding intermediates

Simulating Protein Motions with Rigidity Analysis

Protein motions, ranging from molecular flexibility to large-scale conformational change, play an essential role in many biochemical processes. Despite the explosion in our knowledge of structural and functional data, our understanding of protein movement is still very limited. In previous work, we developed and validated a motion planning based method for mapping protein folding pathways from unstructured conformations to the native state. In this paper, we propose a novel method based on rigidity theory to sample conformation space more effectively, and we describe extensions of our framework to automate the process and to map transitions between specified conformations. Our results show that these additions both improve the accuracy of our maps and enable us to study a broader range of motions for larger proteins. For example, we show that rigidity-based sampling results in maps that capture subtle folding differences between protein G and its mutants, NuG1 and NuG2, and we illustrate how our technique can be used to study large-scale conformational changes in calmodulin, a 148 residue signaling protein known to undergo conformational changes when binding to Ca(2+). Finally, we announce our web-based protein folding server which includes a publicly available archive of protein motions: (http://parasol.tamu.edu/foldingserver/).

Role of unfolded state heterogeneity and en-route ruggedness in protein folding kinetics

In order to improve our understanding of the physical bases of protein folding, there is a compelling need for better connections between experimental and computational approaches. This work addresses the role of unfolded state conformational heterogeneity and en-route intermediates, as an aid for planning and interpreting protein folding experiments. The expected kinetics were modeled for different types of energy landscapes, including multiple parallel folding routes, preferential paths dominated by one primary folding route, and distributed paths with a wide spectrum of microscopic folding rate constants. In the presence of one or more preferential routes, conformational exchange among unfolded state populations slows down the observed rates for native protein formation. We find this to be a general phenomenon, taking place even when unfolded conformations interconvert much faster than the "escape" rate constants to folding. Dramatic kinetic deceleration is expected in the presence of an increasing number of folding-incompetent unfolded conformations. This argues for the existence of parallel folding paths involving several folding-competent unfolded conformations, during the early stages of protein folding. Deviations from single-exponential behavior are observed for unfolded conformations exchanging at comparable rates or more slowly than folding events. Analysis of the effect of en-route (on-path) intermediate formation and landscape ruggedness on folding kinetics leads to the following unexpected conclusions: (1) intermediates, which often retard native state formation, may in some cases accelerate folding, and (2) rugged landscapes, usually associated with stretched exponentials, display single-exponential behavior in the presence of late high-friction paths.

Conformational chemistry of surface-attached calmodulin detected by acoustic shear wave propagation

A thickness shear-mode acoustic wave device, operated in a flow-through format, was used to detect the binding of ions or peptides to surface-attached calmodulin. On-line surface attachment of the protein was achieved by immobilisation of the biotinylated molecule via a neutravidin-biotin linkage onto the surface of the gold electrode of the detector. The interaction between calmodulin, and calcium and magnesium ions induced an increase in resonant frequency and a decrease in motional resistance, which were reversible on washing with buffer. Interestingly, the changes in resonant frequency and motional resistance induced by the binding were opposite to the normal operation of the detector. The response was interpreted as a decrease in surface coupling (partial slip at the liquid/solid interface) instigated by exposure of hydrophobic domains on the protein, and an increase in the thickness, and hence effective wavelength, of the acoustic device, corresponding to an increase in the length of calmodulin by 1.5 A. This result is consistent with the literature value of 4 A. In addition, the interaction of the protein with peptide together with calcium ions was detected successfully, despite the relatively low molecular mass of the 2-kDa peptide. These results confirm the potential of acoustic wave physics for the detection of changes in the conformational chemistry of monolayer of biochemical macromolecules at the solid/liquid interface.

Influence of Fluoro, Chloro and Alkyl Alcohols on the Folding Pathway of Human Serum Albumin

Urea-induced equilibrium unfolding of human serum albumin (HSA) when studied by mean residue ellipticity at 222 nm (MRE(222)) or intrinsic fluorescence measurements showed a two-step, three-state transition with a stable intermediate around 4.6-5.2 M urea. The presence of 2,2,2-trifluoroethanol (TFE) resulted in a single-step, two-state transition with a significant shift towards higher urea concentration, suggesting the stabilizing effect of TFE. The free energy of stabilization (DeltaDeltaG(D)(H(2)O)) in the presence of 3.0 M TFE was determined to be 2.68 and 2.72 kcal/mol by MRE(222) and fluorescence measurements, respectively. The stabilizing potential of other alcohols on the refolding behavior of HSA at 5.0 M urea (where the intermediate exists) as studied by MRE(222) and intrinsic fluorescence measurements showed the following order: 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) > TFE > 2-chloroethanol > tert-butanol > iso-propanol > ethanol > methanol. Further, the extent of refolding at the highest concentration of alcohol was similar in all cases. The stabilizing effect of TFE on guanidine hydrochloride (GdnHCl)-induced unfolding of HSA was nearly equal to that found for urea denaturation, as reflected in the DeltaDeltaG(D)(H(2)O) value (2.38 kcal/mol). Taken together, these results suggest that the stabilizing effect of TFE and other alcohols on urea/GdnHCl-induced unfolding of HSA is higher for alcohols that contain bulky groups or fluorine atoms.

Intimate View of a Kinetic Protein Folding Intermediate: Residue-resolved Structure, Interactions, Stability, Folding and Unfolding Rates, Homogeneity

A cytochrome c kinetic folding intermediate was studied by hydrogen exchange (HX) pulse labeling. Advances in the technique and analysis made it possible to define the structured and unstructured regions, equilibrium stability, and kinetic opening and closing rates, all at an amino acid-resolved level. The entire N-terminal and C-terminal helices are formed and docked together at their normal native positions. They fray in both directions from the interaction region, due to a progression in both unfolding and refolding rates, leading to the surprising suggestion that helix propagation may proceed very slowly in the condensed milieu. Several native-like beta turns are formed. Some residues in the segment that will form the native 60s helix are protected but others are not, suggesting energy minimization to some locally non-native conformation in the transient intermediate. All other regions are unprotected, presumably dynamically disordered. The intermediate resembles a partially constructed native state. It is early, on-pathway, and all of the refolding molecules pass through it. These and related results consistently point to distinct, homogeneous, native-like intermediates in a stepwise sequential pathway, guided by the same factors that determine the native structure. Previous pulse labeling efforts have always assumed EX2 exchange during the labeling pulse, often leading to the suggestion of heterogeneous intermediates in alternative parallel pathways. The present work reveals a dominant role for EX1 exchange in the high pH labeling pulse, which will mimic heterogeneous behavior when EX2 exchange is assumed. The general problem of homogeneous versus heterogeneous intermediates and pathways is discussed.

Reversible denaturation of the soybean Kunitz trypsin inhibitor

The soybean Kunitz trypsin inhibitor (SKTI) is a beta-sheet protein with unusual stability to chemical and thermal denaturation. Different spectroscopic criteria were used to follow the thermal denaturation and renaturation of SKTI. Upon heating to 70 degrees C, changes in UV difference spectra showed increased absorbance at 292 and 297 nm, attributable to perturbation of aromatic residues. Cooling the protein resulted in restoration of the native spectrum unless reduced with dithiothreitol. Far- and near-UV CD spectra also indicate thermal unfolding involving the core tryptophan and tyrosine residues. Both CD and UV-absorbance data suggest a two-state transition with the midpoint at approximately 65 degrees C. CD data along with the increased fluorescence intensity of the reporter fluorophore, 1-anilino-8-naphthalenesulfonate with SKTI, between 60 and 70 degrees C, are consistent with a transition of the native inhibitor to an alternate conformation with a more molten state. Even after heating to 90 degrees C, subsequent cooling of SKTI resulted in >90% of native trypsin inhibition potential. These results indicate that thermal denaturation of SKTI is readily reversible to the native form upon cooling and may provide a useful system for future protein folding studies in the class of disordered beta-sheet proteins.

Partly folded states of members of the lysozyme/lactalbumin superfamily: a comparative study by circular dichroism spectroscopy and limited proteolysis

The partly folded states of protein members of the lysozyme (LYS)/alpha-lactalbumin (LA) superfamily have been analyzed by circular dichroism (CD) measurements and limited proteolysis experiments. Hen, horse, dog, and pigeon LYSs and bovine LA were used in the present study. These are related proteins of 123- to 129-amino-acid residues with similar three-dimensional structures but low similarity in amino acid sequences. Moreover, notable differences among them reside in their calcium-binding properties and capability to adopt partly folded states or molten globules in acid solution (A-state) or on depletion of calcium at neutral pH (apo-state). Far- and near-UV CD measurements revealed that although the structures of hen and dog LYS are rather stable in acid at pH 2.0 or at neutral pH in the absence of calcium, conformational transitions to various extents occur with all other LYS/LA proteins herewith investigated. The most significant perturbation of tertiary structure in acid was observed with bovine LA and LYS from horse milk and pigeon egg-white. Pepsin and proteinase K were used as proteolytic probes, because these proteases show broad substrate specificity, and therefore, their sites of proteolysis are dictated not by the specific amino acid sequence of the protein substrate but by its overall structure and dynamics. Although hen LYS at pH 2.0 was fully resistant to proteolysis by pepsin, the other members of the LYS/LA superfamily were cleaved at different rates at few sites of the polypeptide chain and thus producing rather large protein fragments. The apo-form of bovine LA, horse LYS, and pigeon LYS were attacked by proteinase K at pH 8.3, whereas dog and hen LYSs were resistant to proteolysis when reacted under identical experimental conditions. Briefly, it has been found that the proteolysis data correlate well with the extent of conformational transitions inferred from CD spectra and with existing structural informations regarding the proteins herewith investigated, mainly derived from NMR and hydrogen exchange measurements. The sites of initial proteolytic cleavages in the LYS variants occur at the level of the beta-subdomain (approximately chain region 34-57), in analogy to those observed with bovine LA. Proteolysis data are in agreement with the current view that the molten globule of the LYS/LA proteins is characterized by a structured alpha-domain and a largely disrupted beta-subdomain. Our results underscore the utility of the limited proteolysis approach for analyzing structure and dynamics of proteins, even if adopting an ensemble of dynamic states as in the molten globule.

Molten globule of bovine alpha-lactalbumin at neutral pH induced by heat, trifluoroethanol, and oleic acid: a comparative analysis by circular dichroism spectroscopy and limited proteolysis

The calcium-depleted form of alpha-lactalbumin (alpha-LA) at neutral pH can be induced to adopt a partly folded state or molten globule upon moderate heating, by dissolving the protein in aqueous TFE or by adding oleic acid. This last folding variant of the protein, named HAMLET, can induce apoptosis in tumor cells. The aim of the present work was to unravel from circular dichroism (CD) measurements and proteolysis experiments structural features of the molten globule of apo-alpha-LA at neutral pH. CD spectra revealed that the molten globule of apo-alpha-LA can be obtained upon mild heating at 45 degrees C, as well as at room temperature in the presence of 15% TFE or by adding to the protein solution 7.5 equivalents of oleic acid. Under these various conditions the far- and near-UV CD spectra of apo-alpha-LA are essentially identical to those of the most studied molten globule of alpha-LA at pH 2.0 (A-state). Proteolysis of the 123-residue chain of apo-alpha-LA by proteinase K at 4 degrees C occurs slowly as an all-or-none process leading to small peptides only. At 37 degrees C, proteinase K preferentially cleaves apo-alpha-LA at peptide bonds Ser34-Gly35, Gln39-Ala40, Gln43-Asn44, Phe53-Gln54, and Asn56-Asn57. All these peptide bonds are located at level of the beta-subdomain of the protein (chain region 34-57). Similar sites of preferential cleavage have been observed with the TFE- and oleic acid-induced molten globule of apo-alpha-LA. A protein species given by the N-terminal fragment 1-34 linked via the four disulfide bridges to the C-terminal fragment 54-123 or 57-123 can be isolated from the proteolytic mixture. The results of this study indicate that the same molten globule state of apo-alpha-LA can be obtained at neutral pH under mildly denaturing conditions, as indicated by using a classical spectroscopic technique such as CD and a simple biochemical approach as limited proteolysis. We conclude that the molten globule of alpha-LA maintains a native-like tertiary fold characterized by a rather well-structured alpha-domain and a disordered chain region encompassing the beta-subdomain 34-57 of the protein.

Comparison of protein fragments identified by limited proteolysis and by computational cutting of proteins

Here we present a comparison between protein fragments produced by limited proteolysis and those identified by computational cutting based on the building block folding model. The principles upon which the two methods are based are different. Limited proteolysis of natively folded proteins occurs at flexible sites and never at the level of chain segments of regular secondary structure such as alpha-helices. Therefore, the targets for limited proteolysis are locally unfolded regions. In contrast, the computational cutting algorithm considers the compactness of the fragments, their nonpolar buried surface area, and their isolatedness, that is, the surface area which was buried prior to the cutting and becomes exposed subsequently. Despite the different criteria, there is an overall correspondence between sites or regions of limited proteolysis with those identified by computational cutting. The computational cutting method has been applied to several model proteins for which detailed limited proteolysis data are available, namely apomyoglobin, cytochrome c, ribonuclease A, alpha-lactalbumin, and thermolysin. As expected, more cuts are obtained computationally than experimentally and the agreement is better when a number of proteolytic enzymes are used. For example, cytochrome c is cleaved by thermolysin at 56-57, 45-46, and at 80-81, and by proteinase K at 48-49 and 50-51. Incubation of the noncovalent and native-like complex of cytochrome c fragments 1-56 and 57-104 with proteinase K yielded the gapped protein species 1-48/57-104 and finally 1-40/57-104. Computational cutting of cytochrome c reproduced the major experimental observations, with cuts at 47, 64-65 or 65-66 and 80-81 and an unstable 32-47 region not assigned to any building block. The next step, not addressed in this work, is to probe the ability of the generated fragments to fold independently. Since both the computational algorithm and limited proteolysis attempt to dissect the protein folding problem, the general agreement between the two procedures is gratifying. This consistency allows us to propose the use of limited proteolysis to produce protein fragments that can adopt an independent folding and, therefore, to study folding intermediates. The results of the present study appear to validate the building block folding model and are in line with the proposal that protein folding is a hierarchical process, where parts constituting local minima of energy fold first, with their subsequent association and mutual stabilization to finally yield the global fold.

New stochastic strategy to analyze helix folding

We propose an alternative stochastic strategy to search secondary structures based on the generalized simulated annealing (GSA) algorithm, by using conformational preferences based on the Ramachandran map. We optimize the search for polypeptide conformational space and apply to peptides considered to be good alpha-helix promoters above a critical number of residues. Our strategy to obtain conformational energies consist in coupling a classical force field (THOR package) with the GSA procedure, biasing the Phi x Psi backbone angles to the allowed regions in the Ramachandran map. For polyalanines we obtained stable alpha-helix structures when the number of residues were equal or exceeded 13 amino acids residues. We also observed that the energy gap between the global minimum and the first local minimum tends to increase with the polypeptide size. These conformations were generated by performing 2880 stochastic molecular optimizations with a continuum medium approach. When compared with molecular dynamics or Monte Carlo methods, GSA can be considered the fastest.

Monitoring folding/unfolding transitions of proteins by capillary zone electrophoresis: measurement of deltaG and its variation along the pH scale

Free-solution capillary zone electrophoresis (CZE) can be used to monitor folding/unfolding transitions of proteins and to construct the classical sigmoidal transition curve describing this isomerization process. By performing a series of CZE experiments along the pH scale (here between pH 2.5 and 6.0) it is possible to measure the parameter [urea]1/2, which represents the concentration of urea at the midpoint of each transition curve, and its dependence from the local pH value. The [urea]1/2 parameter provides an idea of the stability of the protein at a given pH; in the case of cytochrome c, for example, it shows that at and below pH 2 the protein will spontaneously unfold even in the absence of a denaturant. The equation describing the sigmoidal folding/unfolding transition can be used for deriving the term deltaG degrees, which refers to the intrinsic difference in the Gibb's free energy between the (total or partial) denatured state and the reference state, taken usually as the native configuration of a protein. The variation of deltaG degrees between the two extremes of our measurements (pH 2.5 and 6.0) along the stated pH interval has been measured (and theoretically calculated) to be of the order of 7-10 kcal/mol and is here interpreted by assuming that at pH 2.5 and below there is an additionally stretching of the polypeptide coil due to coulombic repulsion, as the unfolded chain looses its zwitterionic character and assumes a pure (or very nearly so) cationic surface. Given the minute amounts of sample required, the fully automated state of the analysis, the rapidity and ease of operation, it is hoped that the CZE technique will become more and more popular in the years to come for monitoring folding/unfolding transitions of proteins.

Folding/unfolding/refolding of proteins: present methodologies in comparison with capillary zone electrophoresis

A series of techniques for monitoring protein folding/unfolding/misfolding equilibria are here assessed and compared with capillary zone electrophoresis (CZE). They include spectroscopic techniques, such as circular dichroism, intrinsic fluorescence, nuclear magnetic resonance, Fourier transform infrared and Raman spectroscopy, small-angle X-ray scattering, as well as techniques based on biological assays, such as limited proteolysis and immunochemical analysis of different conformational states. Some unusual probes, such as mass spectrometry for probing unfolding transitions, are also discussed. Size-exclusion chromatography is also evaluated in view of the fact that this technique, like all electrophoretic techniques, and unlike spectroscopic probes, which can only see an average signal in mixed populations, can indeed physically separate folded vs. unfolded macromolecules, especially in the case of slow equilibria. Particular emphasis is devoted to electrophoretic techniques, such as gel-slab electrophoresis in transverse urea or thermal gradients, and CZE. In the latter case, a number of applications are shown, demonstrating the excellent correlation of CZE with more traditional probes, such as intrinsic fluorescence monitoring. It is additionally shown that CZE can be used for measuring the deltaG degrees of unfolding over the pH scale, in good agreement with theoretical calculations on the electrostatic free energy of folding vs. pH, as calculated with a linearized Poisson-Boltzmann equation. Finally, it is demonstrated that CZE can probe also aggregate formation in the presence of helix-inducing agents, such as trifluorethanol.

Stepwise proteolytic removal of the beta subdomain in alpha-lactalbumin. The protein remains folded and can form the molten globule in acid solution

Bovine alpha-lactalbumin (alpha-LA) is an alpha/beta protein which adopts partly folded states when dissolved at low pH (A-state), by removal of the protein-bound calcium at neutral pH and low salt concentration (apo-state), as well as in aqueous trifluoroethanol. Previous spectroscopic studies have indicated that the A-state of alpha-LA at pH 2.0, considered a prototype molten globule, has a native-like fold in which the helical core is mostly retained, while the beta subdomain is less structured. Here, we investigate the conformational features of three derivatives of alpha-LA characterized by a single peptide bond fission or a deletion of 12 or 19/22 amino-acid residues of the beta subdomain of the native protein (approximately from residue 34 to 57). These alpha-LA derivatives were obtained by limited proteolysis of the protein in its partly folded state(s). A nicked alpha-LA species consisting of fragments 1-,3-40 and 41-123 (nicked-LA) was prepared by thermolytic digestion of the 123-residue chain of alpha-LA in 50% (v/v) aqueous trifluoroethanol. Two truncated or gapped protein species given by fragments 1-40 and 53-123 (desbeta1-LA) or fragments 1-34 and 54-,57-123 (desbeta2-LA) were obtained by digestion of alpha-LA with pepsin in acid or with proteinase K at neutral pH in its apo-state, respectively. The two protein fragments of nicked or gapped alpha-LA are covalently linked by the four disulfide bridges of the native protein. CD measurements revealed that, in aqueous solution at neutral pH and in the presence of calcium, the three protein species maintain the helical secondary structure of intact alpha-LA, while the tertiary structure is strongly affected by the proteolytic cleavages of the chain. Temperature effects of CD signals in the far- and near-UV region reveal a much more labile tertiary structure in the alpha-LA derivatives, while the secondary structure is mostly retained even upon heating. In acid solution at pH 2.0, the three alpha-LA variants adopt a conformational state essentially identical to the molten globule displayed by intact alpha-LA, as demonstrated by CD measurements. Moreover, they bind strongly the fluorescent dye 8-anilinonaphthalene-1-sulfonate, which is considered a diagnostic feature of the molten globule of proteins. Therefore, the beta subdomain can be removed from the alpha-LA molecule without impairing the capability of the rest of the chain to adopt a molten globule state. The results of this protein dissection study provide direct experimental evidence that in the alpha-LA molten globule only the alpha domain is structured.

The structural basis of protein folding and its links with human disease

The ability of proteins to fold to their functional states following synthesis in the intracellular environment is one of the most remarkable features of biology. Substantial progress has recently been made towards understanding the fundamental nature of the mechanism of the folding process. This understanding has been achieved through the development and concerted application of a variety of novel experimental and theoretical approaches to this complex problem. The emerging view of folding is that it is a stochastic process, but one biased by the fact that native-like interactions between residues are on average more stable than non-native ones. The sequences of natural proteins have emerged through evolutionary processes such that their unique native states can be found very efficiently even in the complex environment inside a living cell. But under some conditions proteins fail to fold correctly, or to remain correctly folded, in living systems, and this failure can result in a wide range of diseases. One group of diseases, known as amyloidoses, which includes Alzheimer's and the transmissible spongiform encephalopathies, involves deposition of aggregated proteins in a variety of tissues. These diseases are particularly intriguing because evidence is accumulating that the formation of the highly organized amyloid aggregates is a generic property of polypeptides, and not simply a feature of the few proteins associated with recognized pathological conditions. That such aggregates are not normally found in properly functional biological systems is again a testament to evolution, in this case of a variety of mechanisms inhibiting their formation. Understanding the nature of such protective mechanisms is a crucial step in the development of strategies to prevent and treat these debilitating diseases.

Limited proteolysis of bovine alpha-lactalbumin: isolation and characterization of protein domains

The partly folded states of alpha-lactalbumin (alpha-LA) exposed to acid solution at pH 2.0 (A-state) or at neutral pH upon EDTA-mediated removal of the single protein-bound calcium ion (apo form) have been probed by limited proteolysis experiments. These states are nowadays commonly considered to be molten globules and thus protein-folding intermediates. Pepsin was used for proteolysis at acid pH, while proteinase K and chymotrypsin at neutral pH. The expectations were that these proteolytic probes would detect sites and/or chain regions in the partly folded states of alpha-LA sufficiently dynamic, or even unfolded, capable of binding and adaptation to the specific stereochemistry of the protease's active site. A time-course analysis of the proteolytic events revealed that the fast, initial proteolytic cuts of the 123-residue chain of alpha-LA in its A-state or apo form by the three proteases occur at the same chain region 39-54, the actual site(s) of cleavage depending upon the protease employed. This region in native alpha-LA encompasses the beta-sheets of the protein. Subsequent cleavages occur mostly at chain regions 31-35 and 95-105. Four fragment species of alpha-LA have been isolated by reverse-phase high-performance liquid chromatography, and their conformational properties examined by circular dichroism and fluorescence emission spectroscopy. The single chain fragment 53-103, containing all the binding sites for calcium in native alpha-LA and cross-linked by two disulfide bridges, maintains in aqueous buffer and in the presence of calcium ions a folded structure characterized by the same content of alpha-helix of the corresponding chain segment in native alpha-LA. Evidence for some structure was also obtained for the two-chain species 1-40 and 104-123, as well as 1-31 and 105-123, both systems being covalently linked by two disulfide bonds. In contrast, the protein species given by fragment 1-34 connected to fragment 54-123 or 57-123 via four disulfide bridges adopts in solution a folded structure with the helical content expected for a native-like conformation. Of interest, the proteolytic fragment species herewith isolated correspond to the structural domains and subdomains of alpha-LA that can be identified by computational analysis of the three-dimensional structure of native alpha-LA (Siddiqui AS, Barton GI, 1995, Protein Sci 4:872-884). The fast, initial cleavages at the level of the beta-sheet region of native alpha-LA indicate that this region is highly mobile or even unfolded in the alpha-LA molten globule(s), while the rest of the protein chain maintains sufficient structure and rigidity to prevent extensive proteolysis. The subsequent cleavages at chain segment 95-105 indicate that also this region is somewhat mobile in the A-state or apo form of the protein. It is concluded that the overall domain topology of native alpha-LA is maintained in acid or at neutral pH upon calcium depletion. Moreover, the molecular properties of the partly folded states of alpha-LA deduced here from proteolysis experiments do correlate with those derived from previous NMR and other physicochemical measurements.

Mapping the interactions present in the transition state for unfolding/folding of FKBP12

The structure of the transition state for folding/unfolding of the immunophilin FKBP12 has been characterised using a combination of protein engineering techniques, unfolding kinetics, and molecular dynamics simulations. A total of 34 mutations were made at sites throughout the protein to probe the extent of secondary and tertiary structure in the transition state. The transition state for folding is compact compared with the unfolded state, with an approximately 30 % increase in the native solvent-accessible surface area. All of the interactions are substantially weaker in the transition state, as probed by both experiment and molecular dynamics simulations. In contrast to some other proteins of this size, no element of structure is fully formed in the transition state; instead, the transition state is similar to that found for smaller, single-domain proteins, such as chymotrypsin inhibitor 2 and the SH3 domain from alpha-spectrin. For FKBP12, the central three strands of the beta-sheet, beta-strand 2, beta-strand 4 and beta-strand 5, comprise the most structured region of the transition state. In particular Val101, which is one of the most highly buried residues and located in the middle of the central beta-strand, makes approximately 60 % of its native interactions. The outer beta-strands and the ends of the central beta-strands are formed to a lesser degree. The short alpha-helix is largely unstructured in the transition state, as are the loops. The data are consistent with a nucleation-condensation model of folding, the nucleus of which is formed by side-chains within beta-strands 2, 4 and 5, and the C terminus of the alpha-helix. The precise residues involved in the nucleus differ in the two simulated transition state ensembles, but the interacting regions of the protein are conserved. These residues are distant in the primary sequence, demonstrating the importance of tertiary interactions in the transition state. The two independently derived transition state ensembles are structurally similar, which is consistent with a Bronsted analysis confirming that the transition state is an ensemble of states close in structure.

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).

Fast events in protein folding: the time evolution of primary processes

Most experimental studies on the dynamics of protein folding have been confined to timescales of 1 ms and longer. Yet it is obvious that many phenomena that are obligatory elements of the folding process occur on much faster timescales. For example, it is also now clear that the formation of secondary and tertiary structures can occur on nanosecond and microsecond times, respectively. Although fast events are essential to, and sometimes dominate, the overall folding process, with a few exceptions their experimental study has become possible only recently with the development of appropriate techniques. This review discusses new approaches that are capable of initiating and monitoring the fast events in protein folding with temporal resolution down to picoseconds. The first important results from those techniques, which have been obtained for the folding of some globular proteins and polypeptide models, are also discussed.

Unfolding and refolding of dimeric creatine kinase equilibrium and kinetic studies

Equilibrium and kinetic studies of the guanidine hydrochloride induced unfolding-refolding of dimeric cytoplasmic creatine kinase have been monitored by intrinsic fluorescence, far ultraviolet circular dichroism, and 1-anilinonaphthalene-8-sulfonate binding. The GuHCl induced equilibrium-unfolding curve shows two transitions, indicating the presence of at least one stable equilibrium intermediate in GuHCl solutions of moderate concentrations. This intermediate is an inactive monomer with all of the thiol groups exposed. The thermodynamic parameters obtained by analysis using a three-state model indicate that this intermediate is similar in energy to the fully unfolded state. There is a burst phase in the refolding kinetics due to formation of an intermediate within the dead time of mixing (15 ms) in the stopped-flow apparatus. Further refolding to the native state after the burst phase follows biphasic kinetics. The properties of the burst phase and equilibrium intermediates were studied and compared. The results indicate that these intermediates are similar in some respects, but different in others. Both are characterized by pronounced secondary structure, compact globularity, exposed hydrophobic surface area, and the absence of rigid side-chain packing, resembling the "molten globule" state. However, the burst phase intermediate shows more secondary structure, more exposed hydrophobic surface area, and more flexible side-chain packing than the equilibrium intermediate. Following the burst phase, there is a fast phase corresponding to folding of the monomer to a compact conformation. This is followed by rapid assembly to form the dimer. Neither of the equilibrium unfolding transitions are protein concentration dependent. The refolding kinetics are also not concentration dependent. This suggests that association of the subunits is not rate limiting for refolding, and that under equilibrium conditions, dissociation occurs in the region between the two unfolding transitions. Based upon the above results, schemes of unfolding and refolding of creatine kinase are proposed.

Amide protection in an early folding intermediate of cytochrome c

BACKGROUND

For many proteins, compact states appear long before the rate-limiting step in the formation of the native structure. A key issue is whether the initial collapse of the chain is driven by random or more specific hydrophobic interactions.

RESULTS

Hydrogen-exchange labeling coupled with NMR was used to monitor the formation of stable hydrogen-bonded and solvent-excluded structure in horse cytochrome c (cyt c). Protection was measured using a hydrogen exchange/folding competition protocol at variable pH and short competition time (2 ms). Protection factors of threefold to eightfold were observed in all three alpha helices of cyt c, whereas other regions showed no significant protection. This suggests that the compact states that are present contain segments of marginally stable hydrogen-bonded structure. When the intermediate(s) are destabilized, only amide protons from Cys14, Ala15 and His18 show significant protection, indicating a region of persistent residual structure near the covalently bound heme group in the unfolded protein. Fluorescence-detected stopped-flow studies showed that the maximum protection factor in the early intermediate is consistent with its unfolding equilibrium constant.

CONCLUSIONS

Together with previous fluorescence and CD results, the observed pattern of amide protection is consistent with the early formation of an alpha-helical core domain in an ensemble of compact states, indicating that efficient folding is facilitated by stepwise acquisition of native structural elements. These specific early interactions are established on the sub-millisecond time scale, prior to the rate-limiting step for folding.

Detection of conformational changes in an immobilized protein using surface plasmon resonance

Utilizing surface plasmon resonance (SPR), we have developed novel methodology for the detection of conformational change(s) in immobilized proteins. A genetically altered E. coli dihydrofolate reductase (DHFR-ASC) was attached to a carboxymethyldextran matrix layer covering the sensor surface of an SPR biosensor through a disulfide linkage at the engineered protein's C-terminus. The DHFR-ASC-immobilized surface exhibited a larger response to acid treatment than reference surfaces lacking immobilized proteins. The SPR signal of the tethered protein and the molar ellipticity of DHFR-ASC in solution responded similarly to pH changes, consistent with the interpretation that changes in the SPR signal reflect conformational changes occurring during acid denaturation. A pH shift observed between the SPR signal and ellipticity changes may reflect a difference between surface and bulk pH. The tethered protein sensor surface was stable to repeated acid treatment using solutions in the pH range of 0.12-7.80 and yielded reproducible measurements. This is the first demonstration of detection of conformational changes in an immobilized protein using an SPR biosensor. This technique has potential for developing novel sensors and/or switching devices in response to protein conformational changes.

Evidence for barrier-limited protein folding kinetics on the microsecond time scale

Although important structural events in protein folding are known to occur on the submillisecond time scale, the limited time resolution of conventional kinetic methods has precluded direct observation of the initial collapse of the polypeptide chain. A continuous-flow capillary mixing method recently developed by us made it possible to account for the entire fluorescence change associated with refolding of cytochrome c from approximately 5-10(-5)-10(2) s, including the previously unresolved quenching of Trp 59 fluorescence (burst phase) indicative of the formation of compact states. The kinetics of folding exhibits a major exponential process with a time constant of approximately 50 micros, independent of initial conditions and heme ligation state, indicating that a common free energy barrier is encountered during the initial collapse of the polypeptide chain. The resulting loosely packed intermediate accumulates prior to the rate-limiting formation of specific tertiary interactions, confirming previous indications that folding involves at least two distinct stages.

The origin of the alpha-domain intermediate in the folding of hen lysozyme

Stopped-flow fluorescence and circular dichroism spectroscopy have been used in conjunction with quenched-flow hydrogen exchange labelling, monitored by electrospray ionization mass spectrometry, to compare the refolding kinetics of hen egg-white lysozyme at 20 degrees C and 50 degrees C. At 50 degrees C there is clear evidence for distinct fast and slow refolding populations, as observed at 20 degrees C, although folding occurs significantly more rapidly. The folding process is, however, substantially more cooperative at the higher temperature. In particular, the transient intermediate on the major refolding pathway at 20 degrees C, having persistent native-like structure in the alpha-helical domain of the protein, is not detected by hydrogen exchange labelling at 50 degrees C. In addition, the characteristic maximum in negative ellipticity and the minimum in fluorescence intensity observed in far UV CD and intrinsic fluorescence experiments at 20 degrees C, respectively, are not seen at 50 degrees C. Addition of 2 M NaCl to the refolding buffer at 50 degrees C, however, regenerates both the hydrogen exchange and optical properties associated with the alpha-domain intermediate but has no significant effect on the overall refolding kinetics. Together with previous findings, these results indicate that non-native interactions within the alpha-domain intermediate are directly responsible for the unusual optical properties observed during refolding, and that this intermediate accumulates as a consequence of its intrinsic stability in a folding process where the formation of stable structure in the beta-domain constitutes the rate-limiting step for the majority of molecules.

Definition of amide protection factors for early kinetic intermediates in protein folding

Hydrogen-deuterium exchange experiments have been used previously to investigate the structures of well defined states of a given protein. These include the native state, the unfolded state, and any intermediates that can be stably populated at equilibrium. More recently, the hydrogen-deuterium exchange technique has been applied in kinetic labeling experiments to probe the structures of transiently formed intermediates on the kinetic folding pathway of a given protein. From these equilibrium and nonequilibrium studies, protection factors are usually obtained. These protection factors are defined as the ratio of the rate of exchange of a given backbone amide when it is in a fully solvent-exposed state (usually obtained from model peptides) to the rate of exchange of that amide in some state of the protein or in some intermediate on the folding pathway of the protein. This definition is straightforward for the case of equilibrium studies; however, it is less clear-cut for the case of transient kinetic intermediates. To clarify the concept for the case of burst-phase intermediates, we have introduced and mathematically defined two different types of protection factors: one is P struc, which is more related to the structure of the intermediate, and the other is P app, which is more related to the stability of the intermediate. Kinetic hydrogen-deuterium exchange data from disulfide-intact ribonuclease A and from cytochrome c are discussed to explain the use and implications of these two definitions.

Molecular dynamics simulations of hydrophobic collapse of ubiquitin

Nine nonnative conformations of ubiquitin, generated during two different thermal denaturation trajectories, were simulated under nearly native conditions (62 degrees C). The simulations included all protein and solvent atoms explicitly, and simulation times ranged from 1-2.4 ns. The starting structures had alpha-carbon root-mean-square deviations (RMSDs) from the crystal structure of 4-12 A and radii of gyration as high as 1.3 times that of the native state. In all but one case, the protein collapsed when the temperature was lowered and sampled conformations as compact as those reached in a control simulation beginning from the crystal structure. In contrast, the protein did not collapse when simulated in a 60% methanol:water mixture. The behavior of the protein depended on the starting structure: during simulation of the most native-like starting structures (<5 A RMSD to the crystal structure) the RMSD decreased, the number of native hydrogen bonds increased, and the secondary and tertiary structure increased. Intermediate starting structures (5-10 A RMSD) collapsed to the radius of gyration of the control simulation, hydrophobic residues were preferentially buried, and the protein acquired some native contacts. However, the protein did not refold. The least native starting structures (10-12 A RMSD) did not collapse as completely as the more native-like structures; instead, they experienced large fluctuations in radius of gyration and went through cycles of expansion and collapse, with improved burial of hydrophobic residues in successive collapsed states.

Folding intermediates in cytochrome c

Folding of cytochrome c from its low pH guanidine hydrochloride (Gdn-HCl) denatured state revealed a new intermediate, a five-coordinate high spin species with a water molecule coordinated to the heme. Incorporation of this five-coordinated intermediate into the previously reported ligand exchange model can quantitatively account for the observed folding kinetics. In this new model, unfolded cytochrome c is converted to its native structure through an obligatory folding intermediate, the histidine-water coordination state, whereas the five-coordinate state and a bis-histidine state are off-pathway intermediates. When the concentration of Gdn-HCl in the refolding solution was increased, an acceleration of the conversion from the bis-histidine coordinated state to the histidine-water coordinated state was observed, demonstrating that the reaction requires unfolding of the mis-organized polypeptide structure associated with the bis-histidine state.

Snake venom cardiotoxins-structure, dynamics, function and folding

Snake cardiotoxins are highly basic (pI > 10) small molecular weight (approximately 6.5 kDa), all beta-sheet proteins. They exhibit a broad spectrum of interesting biological activities. The secondary structural elements in these toxins include antiparallel double and triple stranded beta-sheets. The three dimensional structures of these toxins reveal an unique asymmetric distribution of the hydrophobic and hydrophilic amino acids. The 3D structures of closely related snake venom toxins such as neurotoxins and cardiotoxin-like basic proteins (CLBP) fail to show similar pattern(s) in the distribution of polar and nonpolar residues. Recently, many novel biological activities have been reported for cardiotoxins. However, to-date, there is no clear structure-function correlation(s) available for snake venom cardiotoxins. The aim of this comprehensive review is to summarize and critically evaluate the progress in research on the structure, dynamics, function and folding aspects of snake venom cardiotoxins.

Protein folding and intermediates

With the exception of the discovery of the rate of formation of the earliest intermediates, there have been no major conceptual leaps in our understanding of protein folding reactions over the past two years. Rather, this period has seen an extension of two established techniques: first, mutational analysis combined with a kinetic definition of the energy landscape of the reaction; and second, the use of hydrogen/deuterium exchange of backbone amide groups combined with NMR. Owing to the application of these methods to a wider range of proteins, it is now possible to draw some general conclusions about the physical processes that direct a protein to its native fold.

Detection of residue contacts in a protein folding intermediate

Protein folding can be described in terms of the development of specific contacts between residues as a highly disordered polypeptide chain converts into the native state. Here we describe an NMR based strategy designed to detect such contacts by observation of nuclear Overhauser effects (NOEs). Experiments with alpha-lactalbumin reveal the existence of extensive NOEs between aromatic and aliphatic protons in the archetypal molten globule formed by this protein at low pH. Analysis of their time development provides direct evidence for near-native compactness of this state. Through a rapid refolding procedure the NOE intensity can be transferred efficiently into the resolved and assigned spectrum of the native state. This demonstrates the viability of using this approach to map out time-averaged interactions between residues in a partially folded protein.

Structural characterisation and comparison of the native and A-states of equine lysozyme

Native state 1H NMR resonance assignments for 125 of the 129 residues of equine lysozyme have enabled measurement of the hydrogen exchange kinetics for over 60 backbone amide and three tryptophan indole hydrogen atoms in the native state. Native holo equine lysozyme hydrogen exchange protection factors are as large as 10(6), the most protected residues being located in elements of secondary structure. High exchange protection in the domain interface correlates with the binding of Ca2+ in this region. Equine lysozyme differs from most non-Ca2+ binding lysozymes in forming a highly populated partially folded state at low pH. The protein in this A-state at pH 2.0 has been found to bind 1-anilino-naphthalene-8-sulphonate with the enhancement of fluorescent intensity and blue shift in the spectral maximum characteristic of molten globules. NMR spectra indicate that the A-state is globally much less ordered than native equine lysozyme but does not contain significant regions of random coil structure. The amides most protected against hydrogen exchange in the A-state (protection factors up to 10(2) at 5 degrees C) correspond to residues of three of the four alpha-helices of the native state; the side-chains of these residues form a hydrophobic cluster that includes five aromatic residues. Circular dichroism and tryptophan fluorescence indicate that these residues are substantially more constrained than similar residues in "classical" molten globules. Taken together, the data suggest a model for the A-state of equine lysozyme in which a more ordered core is surrounded by a less ordered but still compact polypeptide chain.

Fast events in protein folding: relaxation dynamics of secondary and tertiary structure in native apomyoglobin

We report the fast relaxation dynamics of "native" apomyoglobin (pH 5.3) following a 10-ns, laser-induced temperature jump. The structural dynamics are probed using time-resolved infrared spectroscopy. The infrared kinetics monitored within the amide I absorbance of the polypeptide backbone exhibit two distinct relaxation phases which have different spectral signatures and occur on very different time scales (nu = 1633 cm(-1),tau = 48 ns; nu = 1650 cm(-1),tau = 132 micros). We assign these two spectral components to discrete substructures in the protein: helical structure that is solvated (1633 cm(-1)) and native helix that is protected from solvation by interhelix tertiary interactions (1650 cm(-1)). Folding rate coefficients inferred from the observed relaxations at 60 degrees C are k(f)(solvated) = (7 to 20) x 10(6) s(-1) and k(f)(native) = 3.6 x 10(3) s(-1), respectively. The faster rate is interpreted as the intrinsic rate of solvated helix formation, whereas the slower rate is interpreted as the rate of formation of tertiary contacts that determine a native helix. Thus, at 60 degrees C helix formation precedes the formation of tertiary structure by over three orders of magnitude in this protein. Furthermore, the distinct thermodynamics and kinetics observed for the apomyoglobin substructures suggest that they fold independently, or quasi-independently. The observation of inhomogeneous folding for apomyoglobin is remarkable, given the relatively small size and structural simplicity of this protein.

Kinetic role of early intermediates in protein folding

The traditional view that partly folded intermediates are important for directing a protein toward the native state has been challenged by the notion that proteins can intrinsically fold rapidly in a single step if kinetic complications due to slow conformational events are avoided. Intermediates that accumulate within the first few milliseconds of folding are, however, a common observation even for small single-domain proteins. Recent spectroscopic studies, coupled with quantitative kinetic analysis, suggest that folding is facilitated by the rapid formation of compact intermediates with some native-like structural features.

Folding of cytochrome c initiated by submillisecond mixing

Cytochrome c folding was initiated using a new solution mixer that provides a time window which covers over 90% of the burst phase unresolved by conventional stop-flow measurements. Folding was followed by resonance Raman scattering. Kinetic analysis of the high frequency Raman data indicates that a nascent phase occurs within the mixing dead time of 100 microseconds. A significant fraction of the protein was found to be trapped in a misfolded bis-histidine form during the nascent phase at pH 4.5, thereby preventing the protein from folding rapidly and homogeneously. The nascent phase was followed by a haem-ligand exchange phase that populates the native histidine-methionine coordinated form through a thermodynamically controlled equilibrium.

Probing the partly folded states of proteins by limited proteolysis

The folding of a polypeptide chain of a relatively large globular protein into its unique three-dimensional and functionally active structure occurs via folding intermediates. These partly folded states of proteins are difficult to characterize, because they are usually short lived or exist as a distribution of possible conformers. A variety of experimental techniques and approaches have been utilized in recent years in numerous laboratories for characterizing folding intermediates that occur at equilibrium, including spectroscopic techniques, solution X-ray scattering, calorimetry and gel filtration chromatography, as well as genetic methods and theoretical calculations. In this review, we focus on the use of proteolytic enzymes as probes of the structure and dynamics of folding intermediates and we show that this simple biochemical technique can provide useful information, complementing that obtained by other commonly used techniques and approaches. The key result of the proteolysis experiments is that partly folded states (molten globules) of proteins can be sufficiently rigid to prevent extensive proteolysis and appear to maintain significant native-like structure.

Kinetic intermediates in the formation of the cytochrome c molten globule

The relationship between molten globules and transient intermediates in protein folding has been explored by equilibrium and kinetic analysis of the compact acid-denatured A-state of cytochrome c. The chloride-induced formation of the A-state is a complex reaction with structural intermediates resembling those found under native refolding conditions, including a rapidly formed compact state and a subsequent intermediate with interacting N- and C-terminal helices. Together with mutational evidence for specific helix-helix packing interactions, this shows that the A-state is a stable analogue of a late folding intermediate. The L94A mutation blocks all folding steps after the initial collapse and its equilibrium state resembles early kinetic intermediates.

Protein folding monitored at individual residues during a two-dimensional NMR experiment

An approach is described to monitor directly at the level of individual residues the formation of structure during protein folding. A two-dimensional heteronuclear nuclear magnetic resonance (NMR) spectrum was recorded after the rapid initiation of the refolding of a protein labeled with nitrogen-15. The intensities and line shapes of the cross peaks in the spectrum reflected the kinetic time course of the folding events that occurred during the spectral accumulation. The method was used to demonstrate the cooperative nature of the acquisition of the native main chain fold of apo bovine alpha-lactalbumin. The general approach, however, should be applicable to the investigation of a wide range of chemical reactions.

Time-resolved biophysical methods in the study of protein folding

Many of the biophysical techniques developed to characterize native proteins at equilibrium have now been adapted to the structural and thermodynamic characterization of transient intermediate populations during protein folding. Recent advances in these techniques, the use of novel methods of initiating refolding, and a convergence of theoretical and experimental approaches are leading to a detailed understanding of many aspects of the folding process.