The barriers in protein folding

An amino acid code for protein folding

Direct structural information obtained for many proteins supports the following conclusions. The amino acid sequences of proteins can stabilize not only the final native state but also a small set of discrete partially folded native-like intermediates. Intermediates are formed in steps that use as units the cooperative secondary structural elements of the native protein. Earlier intermediates guide the addition of subsequent units in a process of sequential stabilization mediated by native-like tertiary interactions. The resulting stepwise self-assembly process automatically constructs a folding pathway, whether linear or branched. These conclusions are drawn mainly from hydrogen exchange-based methods, which can depict the structure of infinitesimally populated folding intermediates at equilibrium and kinetic intermediates with subsecond lifetimes. Other kinetic studies show that the polypeptide chain enters the folding pathway after an initial free-energy-uphill conformational search. The search culminates by finding a native-like topology that can support forward (native-like) folding in a free-energy-downhill manner. This condition automatically defines an initial transition state, the search for which sets the maximum possible (two-state) folding rate. It also extends the sequential stabilization strategy, which depends on a native-like context, to the first step in the folding process. Thus the native structure naturally generates its own folding pathway. The same amino acid code that translates into the final equilibrium native structure-by virtue of propensities, patterning, secondary structural cueing, and tertiary context-also produces its kinetic accessibility.

Structural basis for the appearance of a molten globule state in chimeric molecules derived from lysozyme and ?-lactalbumin

The problem as to why alpha-lactalbumin, in the absence of Ca(2+), forms a molten globule intermediate, in contrast to its structural homologue lysozyme, has been addressed by the construction of chimeras of human lysozyme in which either the Ca(2+)-binding loop or a part of helix C of bovine alpha-lactalbumin were transplanted. Previously, we have shown that the introduction of both structural elements together in the lysozyme matrix causes the apo form of the resulting chimera to display molten globule behavior during the course of thermal denaturation. In this article, we demonstrate that this molten globule character is not correlated with the Ca(2+)-binding loop. Also, the Del 101 mutant in which Arg101 was deleted to simulate the alpha-lactalbumin conformation of the connecting loop between helix C and helix D, does not show a stable equilibrium intermediate. Rather, the molten globule character of the chimeras has to be related with a specific part of helix C. More particularly, attention is drawn to the four hydrophobic side-chains I93, V96, I99, and L100, the lysozyme counterparts of which are constituted of less bulky valines and alanine. Our observations are discussed in terms of decreased stability of the native form and increased stability of the intermediate molten globule.

Integrity of thermus thermophilus cytochrome c552 synthesized by Escherichia coli cells expressing the host-specific cytochrome c maturation genes, ccmABCDEFGH: biochemical, spectral, and structural characterization of the recombinant protein

We describe the design of Escherichia coli cells that synthesize a structurally perfect, recombinant cytochrome c from the Thermus thermophilus cytochrome c552 gene. Key features are (1) construction of a plasmid-borne, chimeric cycA gene encoding an Escherichia coli-compatible, N-terminal signal sequence (MetLysIleSerIleTyrAlaThrLeu AlaAlaLeuSerLeuAlaLeuProAlaGlyAla) followed by the amino acid sequence of mature Thermus cytochrome c552; and (2) coexpression of the chimeric cycA gene with plasmid-borne, host-specific cytochrome c maturation genes (ccmABCDEFGH). Approximately 1 mg of purified protein is obtained from 1 L of culture medium. The recombinant protein, cytochrome rsC552, and native cytochrome c552 have identical redox potentials and are equally active as electron transfer substrates toward cytochrome ba3, a Thermus heme-copper oxidase. Native and recombinant cytochromes c were compared and found to be identical using circular dichroism, optical absorption, resonance Raman, and 500 MHz 1H-NMR spectroscopies. The 1.7 A resolution X-ray crystallographic structure of the recombinant protein was determined and is indistinguishable from that reported for the native protein (Than, ME, Hof P, Huber R, Bourenkov GP, Bartunik HD, Buse G, Soulimane T, 1997, J Mol Biol 271:629-644). This approach may be generally useful for expression of alien cytochrome c genes in E. coli.

Two-state vs. multistate protein unfolding studied by optical melting and hydrogen exchange

A direct conflict between the stabilization free energy parameters of cytochrome c determined by optical methods and by hydrogen exchange (HX) is quantitatively explained when the partially folded intermediates seen by HX are taken into account. The results support the previous HX measurements of intermediate populations, show how intermediates can elude the standard melting analysis, and illustrate how they confuse the analysis when they are significantly populated within the melting transition region.

Polymer principles of protein calorimetric two-state cooperativity

The experimental calorimetric two-state criterion requires the van't Hoff enthalpy DeltaH(vH) around the folding/unfolding transition midpoint to be equal or very close to the calorimetric enthalpy DeltaH(cal) of the entire transition. We use an analytical model with experimental parameters from chymotrypsin inhibitor 2 to elucidate the relationship among several different van't Hoff enthalpies used in calorimetric analyses. Under reasonable assumptions, the implications of these DeltaH(vH)'s being approximately equal to DeltaH(cal) are equivalent: Enthalpic variations among denatured conformations in real proteins are much narrower than some previous lattice-model estimates, suggesting that the energy landscape theory "folding to glass transition temperature ratio" T(f) /T(g) may exceed 6.0 for real calorimetrically two-state proteins. Several popular three-dimensional lattice protein models, with different numbers of residue types in their alphabets, are found to fall short of the high experimental standard for being calorimetrically two-state. Some models postulate a multiple-conformation native state with substantial pre-denaturational energetic fluctuations well below the unfolding transition temperature, or predict a significant post-denaturational continuous conformational expansion of the denatured ensemble at temperatures well above the transition point, or both. These scenarios either disagree with experiments on protein size and dynamics, or are inconsistent with conventional interpretation of calorimetric data. However, when empirical linear baseline subtractions are employed, the resulting DeltaH(vH)/DeltaH(cal)'s for some models can be increased to values closer to unity, and baseline subtractions are found to correspond roughly to an operational definition of native-state conformational diversity. These results necessitate a re-assessment of theoretical models and experimental interpretations.

Cytochrome c(553), a small heme protein that lacks misligation in its unfolded state, folds with rapid two-state kinetics

Cytochrome c(553) (cyt c(553)) from Desulfovibrio vulgaris is a small helical heme protein that displays apparent two-state equilibrium-unfolding behavior. The covalently attached heme is low-spin, ligated by Met and His residues, in the native state but becomes high-spin upon unfolding at pH 7. Here, we show that in contrast to other c-type heme proteins, where misligations in the unfolded states are prominent, cyt c(553) refolding kinetics at pH 7 proceeds rapidly without detectable intermediates. The extrapolated folding rate constant in water for oxidized cyt c(553) matches exactly that predicted from the cyt c(553) native-state topology: 5300 s(-1 )(experimental) versus 5020 s(-1) (predicted). We therefore conclude that the presence of the oxidized cofactor does not affect the intrinsic formation speed of the cyt c(553 )structural motif.

NMR investigation of ferricytochrome c unfolding: detection of an equilibrium unfolding intermediate and residual structure in the denatured state

Horse ferricytochrome c (cyt c) undergoes exchange of one of its axial heme ligands (Met-80) for one or more non-native ligands under denaturing conditions. We have used (1)H NMR spectroscopy to detect two conformations of paramagnetic cyt c with non-native heme ligation through a range of urea concentrations. One non-native form is an equilibrium unfolding intermediate observed under partially denaturing conditions and is attributed to replacement of Met-80 with one or more Lys side chains. The second non-native form, in which the native Met ligand is replaced by a His, is observed under strongly denaturing conditions. Thermodynamic analysis of these data indicates a relatively small DeltaG (17 kJ/mol) for the transition from native to the Lys-ligated intermediate and a significantly larger DeltaG (47 kJ/mol) for the transition from native to the His-ligated species. Although CD and fluorescence data indicate that the equilibrium unfolding of cyt c is a two-state process, these NMR results implicate an intermediate with His-Lys ligation.

The heme-containing N-fragment (residues 1-56) of cytochrome c is a bis-histidine functional system

The structural and redox properties of a heme-containing fragment (1-56 residues) of cytochrome c have been investigated by spectroscopic (circular dichroism, electronic absorption, and EPR) and voltammetric techniques. The results indicate that the N-fragment lacks ordered secondary structure and has two histidines axially bound to the heme-iron (the native His18 and a misligated His26 or His33). Despite the absence of ordered secondary structure, the peptide chain shields the heme group from solvent, as shown by (i) the pK(a) of protonation of the nonnative histidine ligand (5.18 +/- 0.05), lower than that of the bis-histidine guanidine-unfolded cytochrome c (5.58 +/- 0.05), and (ii) the redox potential, E(o) = 0 +/- 5 mV versus NHE, close to that of bis-histidine cytochrome c mutants but less negative than that of bis-histidine complexes of microperoxidase with short peptides. The electroactive N-fragment may be taken as a "minichrome c" model, with interesting potential for application to biosensor technology; further, the system provides useful information for a deeper understanding of cytochrome c folding and structural/functional organization.

Understanding protein folding via free-energy surfaces from theory and experiment

The ability of protein molecules to fold into their highly structured functional states is one of the most remarkable evolutionary achievements of biology. In recent years, our understanding of the way in which this complex self-assembly process takes place has increased dramatically. Much of the reason for this advance has been the development of energy surfaces (landscapes), which allow the folding reaction to be described and visualized in a meaningful manner. Analysis of these surfaces, derived from the constructive interplay between theory and experiment, has led to the development of a unified mechanism for folding and a recognition of the underlying factors that control the rates and products of the folding process.

Stepwise formation of alpha-helices during cytochrome c folding

Two models have been proposed to describe the folding pathways of proteins. The framework model assumes the initial formation of the secondary structures whereas the hydrophobic collapse model supposes their formation after the collapse of backbone structures. To differentiate between these models for real proteins, we have developed a novel CD spectrometer that enables us to observe the submillisecond time frame of protein folding and have characterized the timing of secondary structure formation in the folding process of cytochrome c (cyt c). We found that approximately 20% of the native helical content was organized in the first phase of folding, which is completed within milliseconds. Furthermore, we suggest the presence of a second intermediate, which has alpha-helical content resembling that of the molten globule state. Our results indicate that many of the alpha-helices are organized after collapse in the folding mechanism of cyt c.

An essential intermediate in the folding of dihydrofolate reductase

The folding of Escherichia coli dihydrofolate reductase was examined at pH 7.8 and 15 degrees C by using stopped-flow fluorescence and absorbance spectroscopies. The formation of a highly fluorescent intermediate occurs with relaxation times ranging between 142 and 343 msec, whereas stopped-flow absorbance spectroscopy using methotrexate binding assays shows a distinct lag phase during these time frames for the native state. The lag in absorbance kinetics and the lack of fast-track folding events indicate that the formation of this ensemble of intermediates is an obligatory step in the folding reaction.

The major transition state in folding need not involve the immobilization of side chains

During protein folding in which few, if any, definable kinetic intermediates are observable, the nature of the transition state is central to understanding the course of the reaction. Current experimental data does not distinguish the relative contributions of side chain immobilization and dehydration phenomena to the major rate-limiting transition state whereas this distinction is central to theoretical models that attempt to simulate the behavior of proteins during folding. Renaturation of the small proteinase inhibitor cystatin under oxidizing versus reducing conditions is the first experimental case in which these processes can be studied independently. Using this example, we show that sidechain immobilization occurs downstream of the major folding transition state. A consequence of this is the existence of states with disordered side chains, which are distinct from kinetic protein folding intermediates and which lie within the folded state free energy well.

Coupled kinetic traps in cytochrome c folding: His-heme misligation and proline isomerization

The effect of His-heme misligation on folding has been investigated for a triple mutant of yeast iso-2 cytochrome c (N26H,H33N,H39K iso-2). The variant contains a single misligating His residue at position 26, a location at which His residues are found in several cytochrome c homologues, including horse, tuna, and yeast iso-1. The amplitude for fast phase folding exhibits a strong initial pH dependence. For GdnHCl unfolded protein at an initial pH<5, the observed refolding at final pH 6 is dominated by a fast phase (tau(2f)=20 ms, alpha(2f)=90 %) that represents folding in the absence of misligation. For unfolded protein at initial pH 6, folding at final pH 6 occurs in a fast phase of reduced amplitude (alpha(2f) approximately 20 %) but the same rate (tau(2f)=20 ms), and in two slower phases (tau(m)=6-8 seconds, alpha(m) approximately 45 %; and tau(1b)=16-20 seconds, alpha(1b) approximately 35 %). Double jump experiments show that the initial pH dependence of the folding amplitudes results from a slow pH-dependent equilibrium between fast and slow folding species present in the unfolded protein. The slow equilibrium arises from coupling of the His protonation equilibrium to His-heme misligation and proline isomerization. Specifically, Pro25 is predominantly in trans in the unligated low-pH unfolded protein, but is constrained in a non-native cis isomerization state by His26-heme misligation near neutral pH. Refolding from the misligated unfolded form proceeds slowly due to the large energetic barrier required for proline isomerization and displacement of the misligated His26-heme ligand.

A new folding intermediate of apomyoglobin from Aplysia limacina: stepwise formation of a molten globule

Apomyoglobin from Aplysia limacina (al-apoMb), despite having only 20 % sequence identity with the more commonly studied mammalian globins (m-apoMbs), properties which result in an increased number of hydrophobic contacts and a loss of most internal salt bridges, shares a number of features of their folding profiles. We show here that it contains an unusually stable core which resists unfolding even at 70 degrees C. The equilibrium intermediate (I(T)) at this high temperature is distinct from the acid unfolded state I(A) which has many properties in common with the acid intermediate observed for the mammalian apoproteins (I(AGH)). It contains a smaller amount of secondary structure (27 % alpha-helical instead of 35 %) and is more highly solvated as evidenced from its fluorescence spectrum (lambda(max)=344 nm instead of 338 nm). Its stability is greatly increased (DeltaDeltaG(w)=-6.75 kcal mol(-1)) in the presence of high salt (2 M KCl), lending support to the view that hydrophobic interactions are responsible for its stability. Kinetic data show classical two-state kinetics between I(A) and the folded state both in the presence and absence of salt. Both I(A) and I(T) can be populated within the dead time of the stopped-flow apparatus, since initiating the refolding reaction from I(T) or I(A) rather than the completely unfolded state does not affect the observed refolding time-course. Our conclusion is that al-apoMb, as other "apo" proteins (including for example alpha-lactalbumin in the absence of Ca(2+)), may be described as "uncoupled" with an unusually high and exploitable tendency to populate partially folded states.

Effects of pH on the kinetic reaction mechanism of myoglobin unfolding studied by time-resolved electrospray ionization mass spectrometry

In most cases, kinetic unfolding reactions of proteins follow a simple one-step mechanism that does not involve any detectable intermediates. One example for a more complicated unfolding reaction is the acid-induced denaturation of holo-myoglobin (hMb). This reaction proceeds through a transient intermediate and can be described by a sequential two-step mechanism (Konermann et al. Biochemistry 1997, 36, 6448-6454). Time-resolved electrospray ionization mass spectrometry (ESI MS) is a new technique for monitoring the kinetics of protein folding and unfolding in solution. Different protein conformations can be distinguished by the different charge state distributions that they generate during ESI. At the same time this technique allows monitoring the loss or binding of noncovalent protein ligands. In this work, time-resolved ESI MS is used to study the dependence of the kinetic unfolding mechanism of hMb on the specific solvent conditions used in the experiment. It is shown that hMb unfolds through a short-lived intermediate only at acidic pH. Under basic conditions no intermediate is observed. These findings are confirmed by the results of optical stopped-flow absorption experiments. This appears to be the first time that a dependence of the kinetic mechanism for protein unfolding on external conditions such as pH has been observed.

Evidence for a requirement for ATP hydrolysis at two distinct steps during a single turnover of the catalytic cycle of human P-glycoprotein

P-glycoprotein (Pgp) is an ATP-dependent hydrophobic natural product anticancer drug efflux pump whose overexpression confers multidrug resistance to tumor cells. The work reported here deals with the elucidation of the energy requirement for substrate interaction with Pgp during the catalytic cycle. We show that the K(d) (412 nM) of the substrate analogue [(125)I]iodoarylazidoprazoin for Pgp is not altered by the presence of the nonhydrolyzable nucleotide 5'-adenylylimididiphosphate and vanadate (K(d) = 403 nM). Though binding of nucleotide per se does not affect interactions with the substrate, ATP hydrolysis results in a dramatic conformational change where the affinity of [(125)I]iodoarylazidoprazoin for Pgp trapped in transition-state conformation (Pgp x ADP x vanadate) is reduced >30-fold. To transform Pgp from this intermediate state of low affinity for substrate to the next catalytic cycle, i.e., a conformation that binds substrate with high affinity, requires conditions that permit ATP hydrolysis. Additionally, there is an inverse correlation (R(2) = 0.96) between 8AzidoADP (or ADP) release and the recovery of substrate binding. These results suggest that the release of nucleotide is necessary for reactivation but not sufficient. The hydrolysis of additional molecule(s) of ATP (or 8AzidoATP) is obligatory for the catalytic cycle to advance to completion. These data are consistent with the observed stoichiometry of two ATP molecules hydrolyzed for the transport of every substrate molecule. Our data demonstrate two distinct roles for ATP hydrolysis in a single turnover of the catalytic cycle of Pgp, one in the transport of substrate and the other in effecting conformational changes to reset the pump for the next catalytic cycle.

Fast kinetics and mechanisms in protein folding

This review describes how kinetic experiments using techniques with dramatically improved time resolution have contributed to understanding mechanisms in protein folding. Optical triggering with nanosecond laser pulses has made it possible to study the fastest-folding proteins as well as fundamental processes in folding for the first time. These include formation of alpha-helices, beta-sheets, and contacts between residues distant in sequence, as well as overall collapse of the polypeptide chain. Improvements in the time resolution of mixing experiments and the use of dynamic nuclear magnetic resonance methods have also allowed kinetic studies of proteins that fold too fast (greater than approximately 10(3) s-1) to be observed by conventional methods. Simple statistical mechanical models have been extremely useful in interpreting the experimental results. One of the surprises is that models originally developed for explaining the fast kinetics of secondary structure formation in isolated peptides are also successful in calculating folding rates of single domain proteins from their native three-dimensional structure.

Antibody-detected folding: kinetics of surface epitope formation are distinct from other folding phases

The rate of macromolecular surface formation in yeast iso-2 cytochrome c and its site-specific mutant, N52I iso-2, has been studied using a monoclonal antibody that recognizes a tertiary epitope including K58 and H39. The results indicate that epitope refolding occurs after fast folding but prior to slow folding, in contrast to horse cytochrome c where surface formation occurs early. The antibody-detected (ad) kinetic phase accompanying epitope formation has k(ad) = 0.2 s(-1) and is approximately 40-fold slower than the fastest detectable event in the folding of yeast iso-2 cytochrome c (k2f approximately 8 s(-1)), but occurs prior to the absorbance- and fluorescence-detected slow folding steps (k1a approximately 0.06 s(-1); k1b approximately 0.09 s(-1)). N5I iso-2 cytochrome c exhibits similar kinetic behavior with respect to epitope formation. A detailed dissection of the mechanistic differences between the folding pathways of horse and yeast cytochromes c identifies possible reasons for the slow surface formation in the latter. Our results suggest that non-native ligation involving H33 or H39 during refolding may slow down the formation of the tertiary epitope in iso-2 cytochrome c. This study illustrates that surface formation can be coupled to early events in protein folding. Thus, the rate of macromolecular surface formation is fine tuned by the residues that make up the surface and the interactions they entertain during refolding.

Redesigning the hydrophobic core of a model beta-sheet protein: destabilizing traps through a threading approach

An off-lattice 46-bead model of a small all-beta protein has been recently criticized for possessing too many traps and long-lived intermediates compared with the folding energy landscape predicted for real proteins and models using the principle of minimal frustration. Using a novel sequence design approach based on threading for finding beneficial mutations for destabilizing traps, we proposed three new sequences for folding in the beta-sheet model. Simulated annealing on these sequences found the global minimum more reliably, indicative of a smoother energy landscape, and simulated thermodynamic variables found evidence for a more cooperative collapse transition, lowering of the collapse temperature, and higher folding temperatures. Folding and unfolding kinetics were acquired by calculating first-passage times, and the new sequences were found to fold significantly faster than the original sequence, with a concomitant lowering of the glass temperature, although none of the sequences have highly stable native structures. The new sequences found here are more representative of real proteins and are good folders in the T(f) > T(g) sense, and they should prove useful in future studies of the details of transition states and the nature of folding intermediates in the context of simplified folding models. These results show that our sequence design approach using threading can improve models possessing glasslike folding dynamics.

Intermolecular aggregations are responsible for the slow kinetics observed in the folding of cytochrome c at neutral pH

Folding of equine cytochrome c at a low protein concentration (26 microM) eliminated a slow kinetic phase (time constant three seconds) that was observed in the previous hydrogen exchange pulse-labeling experiments at pH 6.2 and 10 degrees C. It was demonstrated that this slow folding phase was caused by intermolecular aggregations. Because heterogeneous kinetics is a very general feature in the folding of proteins characterized by pulsed hydrogen exchange coupled with two-dimensional NMR, our experimental results suggest aggregations might also be responsible for the complex folding kinetics of other proteins. This is possible since these experiments were performed at relatively high protein concentrations.

The cooperativity of burst phase reactions explored

The denaturant-dependence of the major, observable relaxation rates for folding (kobs) of ribonuclease HI from Escherichia coli (RNase H) and phage T4 lysozyme (T4L) reveal that, for both proteins, folding begins with the rapid and transient accumulation of intermediate species in a "burst phase" which precedes the rate-limiting formation of the native state; this is evidenced by a "rollover" in the folding limb of the rate profiles (kobs versus denaturant, or chevron plot). These rate profiles are most simply described by a three-state mechanism (unfolded-to-intermediate-to-native), which implies that the burst phase represents a transition between two distinct thermodynamic states. It is shown here that the equilibrium properties of these burst phase reactions can be equally well modeled by a mechanism involving a continuum of states where the free energy of each state is linearly related to its m-value (the parameter describing the linear relationship between free energy and denaturant). A numerical model is also developed to describe the time evolution of such a system, which exhibits nearly perfect exponential behavior. Both models emphasize how a continuum of states operating under a linear free energy relationship may behave like a two state system. Such a scheme finds experimental justification from an interpretation of recent native state hydrogen exchange data. The analytical model described for a continuum can account for the observed kinetic profiles of several other model proteins. The results, however, appear context specific, suggesting that burst phase reactions are not entirely random and non-specific. The results reported in this study have important implications for the concept of cooperativity in protein folding reactions.

Folding of a model three-helix bundle protein: a thermodynamic and kinetic analysis

The kinetics and thermodynamics of an off-lattice model for a three-helix bundle protein are investigated as a function of a bias gap parameter that determines the energy difference between native and non-native contacts. A simple dihedral potential is used to introduce the tendency to form right-handed helices. For each value of the bias parameter, 100 trajectories of up to one microsecond are performed. Such statistically valid sampling of the kinetics is made possible by the use of the discrete molecular dynamics method with square-well interactions. This permits much faster simulations for off-lattice models than do continuous potentials. It is found that major folding pathways can be defined, although ensembles with considerable structural variation are involved. The large gap models generally fold faster than those with a smaller gap. For the large gap models, the kinetic intermediates are non-obligatory, while both obligatory and non-obligatory intermediates are present for small gap models. Certain large gap intermediates have a two-helix microdomain with one helix extended outward (as in domain-swapped dimers); the small gap intermediates have more diverse structures. The importance of studying the kinetic, as well as the thermodynamics, of folding for an understanding of the mechanism is discussed and the relation between kinetic and equilibrium intermediates is examined. It is found that the behavior of this model system has aspects that encompass both the "new" view and the "old" view of protein folding.

Kinetic and equilibrium intermediate states are different in LYLA1, a chimera of lysozyme and alpha-lactalbumin

For several proteins, a striking resemblance has been observed between the equilibrium partially folded state and the kinetic burst-phase intermediate, observed just after the dead-time in refolding experiments. This has led to the general statement that the conformation of both types of intermediates is similar. We show, at least for one of the proteins investigated here, that, although both states have some common characteristics, they are not identical. LYLA1 is a chimeric protein resulting from the transplantation of the Ca(2+)-binding loop and the adjacent helix C of bovine alpha-lactalbumin into the homologous position (residues 76-102) in human lysozyme. The apo-form of LYLA1 unfolds through a partially folded state, in analogy with the folding behaviour of the structurally homologous alpha-lactalbumin. The folding kinetics of LYLA1 and of its wild-type homologue, human lysozyme, are investigated by means of stopped-flow fluorescence and CD spectroscopy. In the case of human lysozyme, refolding involves parallel pathways as indicated by experiments in the presence of a fluorescent inhibitor. For apo-LYLA1, the burst-phase intermediate is compared with the equilibrium intermediate. At neutral pH, both states correspond, in that an important amount of secondary structure has been established, but the burst-phase intermediate is shown to be significantly less stable than the equilibrium intermediate. At pH 1.85, in the presence of 1.5 M guanidinium hydrochloride (GdnHCl) and at 25 degrees C, the equilibrium partially folded state of LYLA1 is 100% populated. When LYLA1 is rapidly diluted from 6 M GdnHCl to 1.5 M under these conditions, a time-dependent evolution of the fluorescence signal is observed, reflecting the transition from a burst-phase to a different equilibrium intermediate. These results provide strong evidence for the non-identity of both states in this protein.

Direct kinetic evidence for folding via a highly compact, misfolded state

The 2 S seed storage protein, sunflower albumin 8 (SFA-8), contains an unusually high proportion of hydrophobic residues including 16 methionines (some of which may form a surface hydrophobic patch) in a disulfide cross-linked, alpha-helical structure. Circular dichroism and fluorescence spectroscopy show that SFA-8 is highly stable to denaturation by heating or chaotropic agents, the latter resulting in a reversible two-state unfolding transition. The small m(U) (-4.7 M(-1) at 10 degrees C) and DeltaC(p) (-0.95 kcal mol(-1) K(-1)) values indicate that relatively little nonpolar surface of the protein is exposed during unfolding. Commensurate with the unusual distribution of hydrophobic residues, stopped-flow fluorescence data show that the folding pathway of SFA-8 is highly atypical, in that the initial product of the rapid collapse phase of folding is a compact nonnative state (or collection of nonnative states) that must unfold before acquiring the native conformation. The inhibited folding reaction of SFA-8, in which the misfolded state (m(M) = -0.95 M(-1) at 10 degrees C) is more compact than the transition state for folding (m(T) = -2.5 M(-1) at 10 degrees C), provides direct kinetic evidence for the transient misfolding of a protein.

Compactness of the denatured state of a fast-folding protein measured by submillisecond small-angle x-ray scattering

Time-resolved small-angle x-ray scattering was used to measure the radius of gyration of cytochrome c after initiation of folding by a pH jump. Submillisecond time resolution was obtained with a microfabricated diffusional mixer and synchrotron radiation. The results show that the protein first collapses to compact denatured structures before folding very fast to the native state.

The hydrogen exchange core and protein folding

A database of hydrogen-deuterium exchange results has been compiled for proteins for which there are published rates of out-exchange in the native state, protection against exchange during folding, and out-exchange in partially folded forms. The question of whether the slow exchange core is the folding core (Woodward C, 1993, Trends Biochem Sci 18:359-360) is reexamined in a detailed comparison of the specific amide protons (NHs) and the elements of secondary structure on which they are located. For each pulsed exchange or competition experiment, probe NHs are shown explicitly; the large number and broad distribution of probe NHs support the validity of comparing out-exchange with pulsed-exchange/competition experiments. There is a strong tendency for the same elements of secondary structure to carry NHs most protected in the native state, NHs first protected during folding, and NHs most protected in partially folded species. There is not a one-to-one correspondence of individual NHs. Proteins for which there are published data for native state out-exchange and theta values are also reviewed. The elements of secondary structure containing the slowest exchanging NHs in native proteins tend to contain side chains with high theta values or be connected to a turn/loop with high theta values. A definition for a protein core is proposed, and the implications for protein folding are discussed. Apparently, during folding and in the native state, nonlocal interactions between core sequences are favored more than other possible nonlocal interactions. Other studies of partially folded bovine pancreatic trypsin inhibitor (Barbar E, Barany G, Woodward C, 1995, Biochemistry 34:11423-11434; Barber E, Hare M, Daragan V, Barany G, Woodward C, 1998, Biochemistry 37:7822-7833), suggest that developing cores have site-specific energy barriers between microstates, one disordered, and the other(s) more ordered.

Experimental study of the protein folding landscape: unfolding reactions in cytochrome c

Hydrogen exchange results for cytochrome c have been interpreted in terms of transient hydrogen bond-breaking reactions that include large unfolding reactions and small fluctuational distortions. The differential sensitivity of these opening reactions to denaturant, temperature, and protein stability makes it possible to distinguish the different opening reactions and to characterize their structural and thermodynamic parameters. The partially unfolded forms (PUFs) observed are few and discrete, evidently because they are produced by the reversible unfolding of the protein's several intrinsically cooperative secondary structural elements. The PUFs are robust, evidently because the structural elements do not change over a wide range of conditions. The discrete nature of the PUFs and their small number is as expected for classical folding intermediates but not for theoretically derived folding models apparently because the simplified non-protein models usually analyzed in theoretical studies encompass only a single cooperative unit rather than multiple separable units.

Formation of a molten-globule-like state of cytochrome c induced by high concentrations of glycerol

The effect of glycerol on the structure of cytochrome c was investigated by circular dichroism, absorbance and EPR spectroscopy. The results obtained show that an increasing concentration of the organic solvent (70-99.2%, v/v) in aqueous-polyalcohol mixtures converts native cytochrome c into a new, low spin form through a fully reversible, two-state transition. The glycerol-stabilized form (that we call here the G state) retains native-like amounts of alpha-helix structure while rigid tertiary structure and native Fe(III)-Met(80) axial bond are lost. Analysis of data suggests a molten globule character of the G state; support to this view is afforded by the striking similarities between the spectroscopic (and, thus, structural) properties of the G state with those of the acidic molten globule of the protein (A state).

Equilibrium unfolding of a small low-potential cytochrome, cytochrome c553 from Desulfovibrio vulgaris

To understand general aspects of stability and folding of c-type cytochromes, we have studied the folding characteristics of cytochrome c553 from Desulfovibrio vulgaris (Hildenborough). This cytochrome is structurally similar but lacks sequence homology to other heme proteins; moreover, it has an abnormally low reduction potential. Unfolding of oxidized and reduced cytochrome c553 by guanidine hydrochloride (GuHCl) was monitored by circular dichroism (CD) and Soret absorption; the same unfolding curves were obtained with both methods supporting that cytochrome c553 unfolds by an apparent two-state process. Reduced cytochrome c553 is 7(3) kJ/mol more stable than the oxidized form; accordingly, the reduction potential of unfolded cytochrome c553 is 100(20) mV more negative than that of the folded protein. In contrast to many other unfolded cytochrome c proteins, upon unfolding at pH 7.0 both oxidized and reduced heme in cytochrome c553 become high-spin. The lack of heme misligation in unfolded cytochrome c553 implies that its unfolded structure is less constrained than those of cytochromes c with low-spin, misligated hemes.

Folding mechanism of Pseudomonas aeruginosa cytochrome c551: role of electrostatic interactions on the hydrophobic collapse and transition state properties

We report on the folding kinetics of the small 82 residue cytochrome c551from Pseudomonas aeruginosa. The presence of two Trp residues (Trp56 and Trp77) allows the monitoring of fluorescence quenching on refolding in two different regions of the protein. A single His residue (the iron-coordinating His16) permits the study of refolding in the absence of miscoordination events. After identification of the kinetic traps (Pro isomerization and aggregation of denatured protein), overall refolding kinetics is described by two processes: (i) a burstphase collapse (faster than milliseconds) which we show to be a global event leading to a state whose compactness depends on the overall net charge; at the isoeletric pH (4.7), it is maximally compact, while above and below it is more expanded; and (ii) an exponential phase (in the millisecond time range) leading to the native protein via a transition state(s) possibly involving the formation of a specific salt bridge between Lys10 and Glu70, at the contact between the N and C-terminal helices. Comparison with the widely studied horse cytochrome c allows the discussion of similarities and differences in the folding of two proteins which have the same "fold" despite a very low degree of sequence homology (<30 %).

Chaperonin function: folding by forced unfolding

The ability of the GroEL chaperonin to unfold a protein trapped in a misfolded condition was detected and studied by hydrogen exchange. The GroEL-induced unfolding of its substrate protein is only partial, requires the complete chaperonin system, and is accomplished within the 13 seconds required for a single system turnover. The binding of nucleoside triphosphate provides the energy for a single unfolding event; multiple turnovers require adenosine triphosphate hydrolysis. The substrate protein is released on each turnover even if it has not yet refolded to the native state. These results suggest that GroEL helps partly folded but blocked proteins to fold by causing them first to partially unfold. The structure of GroEL seems well suited to generate the nonspecific mechanical stretching force required for forceful protein unfolding.

Cavity defects in the procapsid of bacteriophage P22 and the mechanism of capsid maturation

Bacteriophage P22 belongs to a family of double-stranded DNA viruses that share common morphogenetic features like DNA packaging into a procapsid precursor and maturation. Maturation involves cooperative expansion of the procapsid shell with concomitant lattice stabilization. The expansion is thought to be mediated by movement of two coat protein domains around a hinge. The metastable conformation of subunit within the procapsid lattice is considered to constitute a late folding intermediate. In order to understand the mechanism of expansion it is necessary to characterize the interactions stabilizing procapsid and mature capsid lattices, respectively. We employ pressure dissociation to compare subunit packing within the procapsid and expanded lattice. Procapsid shells contain larger cavities than the expanded shells, presumably due to polypeptide packing defects. These defects contribute to the metastable nature of the procapsid lattice and are cured during expansion. Improved packing contributes to the increased stability of the expanded shell. Comparison of two temperature-sensitive folding (tsf) mutants of coat protein (T294I and W48Q) with wild-type coat revealed that both mutations markedly destabilized the procapsid shell and yet had little effect on relative stability of the monomeric subunit. Thus, the regions affected by these packing defects constitute subunit interfaces of the procapsid shell. The larger activation volume of pressure dissociation observed for both T294I and W48Q indicates that the decreased stability of these particles is due to increase of cavity defects. These defects in the procapsid lattice are cured upon expansion suggesting that the intersubunit contacts affected by tsf mutations are absent or rearranged in the mature shell. The energetics of the in vitro expansion reaction also suggests that entropic stabilization contributes to the large free energy barrier for expansion.

Chaperone-mediated protein folding

The folding of most newly synthesized proteins in the cell requires the interaction of a variety of protein cofactors known as molecular chaperones. These molecules recognize and bind to nascent polypeptide chains and partially folded intermediates of proteins, preventing their aggregation and misfolding. There are several families of chaperones; those most involved in protein folding are the 40-kDa heat shock protein (HSP40; DnaJ), 60-kDa heat shock protein (HSP60; GroEL), and 70-kDa heat shock protein (HSP70; DnaK) families. The availability of high-resolution structures has facilitated a more detailed understanding of the complex chaperone machinery and mechanisms, including the ATP-dependent reaction cycles of the GroEL and HSP70 chaperones. For both of these chaperones, the binding of ATP triggers a critical conformational change leading to release of the bound substrate protein. Whereas the main role of the HSP70/HSP40 chaperone system is to minimize aggregation of newly synthesized proteins, the HSP60 chaperones also facilitate the actual folding process by providing a secluded environment for individual folding molecules and may also promote the unfolding and refolding of misfolded intermediates.

Deciphering the timescales and mechanisms of protein folding using minimal off-lattice models

Considerable insights into the mechanisms and timescales of protein folding have been obtained from detailed studies of minimal off-lattice models. These models are coarse-grained representations of polypeptide chains. Many novel predictions of the mechanisms and timescales of the folding of proteins have been made using simulations of off-lattice models. The concepts derived from these simulations have been used to analyze the recent experiments and simulations of proteins and peptides.

The Greek key protein apo-pseudoazurin folds through an obligate on-pathway intermediate

Folding of the 123 amino acid residue Greek key protein apo-pseudo azurin from Thiosphaera pantotropha has been examined using stopped-flow circular dichroism in 0.5 M Na2SO4 at pH 7.0 and 15 degrees C. The data show that the protein folds from the unfolded state with all eight proline residues in their native isomers (seven trans and one cis) to an intermediate within the dead-time of the stopped-flow mixing (50 ms). The urea dependence of the rates of folding and unfolding of the protein were also determined. The ratio of the folding rate to the unfolding rate (extrapolated into water) is several orders of magnitude too small to account for the equilibrium stability of the protein, consistent with the population of an intermediate. Despite this, the logarithm of the rate of folding versus denaturant concentration is linear. These data can be rationalised by the population of an intermediate under all refolding conditions. Accordingly, kinetic and equilibrium measurements were combined to fit the chevron plot to an on-pathway model (U <==> I <==> N). The fit shows that apo-pseudoazurin rapidly forms a compact species that is stabilised by 25 kJ/mol before folding to the native state at a rate of 2 s-1. Although the data can also be fitted to an off-pathway model (I <==> U <==> N), the resulting kinetic parameters indicate that the protein would have to fold to the native state at a rate of 86,000 s-1 (a time constant of only 12 microseconds). Similarly, models in which this intermediate is bypassed also lead to unreasonably fast refolding rates. Thus, the intermediate populated during the refolding of apo-pseudoazurin appears to be obligate and on the folding pathway. We suggest, based on this study and others, that some intermediates play a critical role in limiting the search to the native state.

Purification and characterization of a plant antimicrobial peptide expressed in Escherichia coli

MiAMP1 is a low-molecular-weight, cysteine-rich, antimicrobial peptide isolated from the nut kernel of Macadamia integrifolia. A DNA sequence encoding MiAMP1 with an additional ATG start codon was cloned into a modified pET vector under the control of the T7 RNA polymerase promoter. The pET vector was cotransformed together with the vector pSB161, which expresses a rare arginine tRNA. The peptide was readily isolated in high yield from the insoluble fraction of the Escherichia coli extract. The purified peptide was shown to have an identical molecular weight to the native peptide by mass spectroscopy indicating that the N-terminal methionine had been cleaved. Analysis by NMR spectroscopy indicated that the refolded recombinant peptide had a similar overall three-dimensional structure to that of the native peptide. The peptide inhibited the growth of phytopathogenic fungi in vitro in a similar manner to the native peptide. To our knowledge, MiAMP1 is the first antimicrobial peptide from plants to be functionally expressed in E. coli. This will permit a detailed structure-function analysis of the peptide and studies of its mode of action on phytopathogens.

Kinetic evidence for an on-pathway intermediate in the folding of cytochrome c

An early folding event of cytochrome c populates a helix-containing intermediate (INC) because of a pH-dependent misligation between the heme iron and nonnative ligands in the unfolded state (U). For folding to proceed, the nonnative ligation error must first be corrected. It is not known whether I is on-pathway, with folding to the native state (N) as in U <-->INC <--> N, or whether the I must first move back through the U and then fold to the N through some alternative path (INC <--> U <--> N). By means of a kinetic test, it is shown here that the cytochrome c I does not first unfold to U. The method used provides an experimental criterion for rejecting the off-pathway I <--> U <--> N option.

Biophysics of the green fluorescent protein

It is almost certainly a truism that interpretation of the fluorescence of a protein matrix-embedded chromophore in terms of the physicochemical character of its environment requires that the tertiary structure of the protein be known to high resolution. This reality derives from the complexity of the photophysics of most fluorescent molecules, complexity that reveals the imperfections of available theory. The accuracy of these dicta is highlighted by the biophysical properties of the green fluorescent protein now being so elegantly elucidated from the application of X-ray crystallography, ultrafast optical spectroscopy, and site-specific mutagenesis. Despite the mass of recent data, however, the physicochemical basis of the green fluorescence cannot be regarded as having been fully defined, nor has the role of protein folding in chromophore formation been solved. In addition, GFP consistently yields surprises typified by the recent experiments of Vanden Bout et al. (1997) on the fluorescence of single molecules mutant (T203F and T2034) GFPs, which showed unique reversible photobleaching and were whimsically termed "blinking molecules" by Moerner (1977). Given the apparent malleability of the GFP sequence and the sensitivity of the chromophore's photophysics to a broad spectrum of physicochemical factors, it is inevitable that additional useful and intriguing biophysical properties will emerge from the study of other mutants. Although on the surface it may seem mundane, determination of the amino acid sequence and tertiary structures of the GFPs from other coelenterates is quite likely to provide very useful insights into the biophysical bases of both protein folding and the green fluorescence per se. Finally, a broader set of spectroscopic techniques need to be applied to the study of GFPs, and future fluorescence examination should include measurements of transient absorption and fluorescence emission anisotropy decays.

A protein folding intermediate of ribonuclease T1 characterized at high resolution by 1D and 2D real-time NMR spectroscopy

The rate-limiting step during the refolding of S54G/P55N ribonuclease T1 is determined by the slow trans-->cis prolyl isomerisation of Pro39. We investigated the refolding of this variant by one-dimensional (1D) and two-dimensional (2D) real-time NMR spectroscopy, initiated by a tenfold dilution from 6 M guanidine hydrochloride at 10 degreesC. Two intermediates could be resolved with the 1D approach. The minor intermediate, which is only present early during refolding, is largely unfolded. The major intermediate, with an incorrect trans Pro39 peptide bond, is highly structured with 33 amide protons showing native chemical shifts and native NOE patterns. They could be assigned in a real-time 2D-NOESY (nuclear Overhauser enhancement spectroscopy) by using a new assignment strategy to generate positive and negative signal intensities for native and non-native NOE cross-peaks, respectively. Surprisingly, amide protons with non-native environments are located not only close to Tyr38-Pro39, but are spread throughout the entire protein, including the C-terminal part of the alpha-helix, beta-strands 3 and 4 and several loop regions. Native secondary and tertiary structure was found for the major intermediate in the N-terminal beta-strands 1 and 2 and the C terminus (connected by the disulfide bonds), the N-terminal part of the alpha-helix, and the loops between beta-strands 4/5 and 5/6. Implications of these native and non-native structure elements of the intermediate for the refolding of S54G/P55N ribonuclease T1 and for cis/trans isomerizations are discussed.