Trapping the on-pathway folding intermediate of Im7 at equilibrium

The four-helical protein Im7 folds via a rapidly formed on-pathway intermediate (k(UI)=3000 s(-1) at pH 7.0, 10 degrees C) that contains three (helices I, II and IV) of the four native alpha-helices. The relatively slow (k(IN)=300 s(-1)) conversion of this intermediate into the native structure is driven by the folding and docking of the six residue helix III onto the developing hydrophobic core. Here, we describe the structural properties of four Im7* variants designed to trap the protein in the intermediate state by disrupting the stabilising interactions formed between helix III and the rest of the protein structure. In two of these variants (I54A and L53AI54A), hydrophobic residues within helix III have been mutated to alanine, whilst in the other two mutants the sequence encompassing the native helix III was replaced by a glycine linker, three (H3G3) or six (H3G6) residues in length. All four variants were shown to be monomeric, as judged by analytical ultracentrifugation, and highly helical as measured by far-UV CD. In addition, all the variants denature co-operatively and have a stability (DeltaG(UF)) and buried hydrophobic surface area (M(UF)) similar to those of the on-pathway kinetic intermediate. Structural characterisation of these variants using 1-anilino-8-napthalene sulphonic acid (ANS) binding, near-UV CD and 1D (1)H NMR demonstrate further that the trapped intermediate ensemble is highly structured with little exposed hydrophobic surface area. Interestingly, however, the structural properties of the variants I54A and L53AI54A differ in detail from those of H3G3 and H3G6. In particular, the single tryptophan residue, located near the end of helix IV, and distant from helix III, is in a distinct environment in the two sets of mutants as judged by fluorescence, near-UV CD and the sensitivity of tryptophan fluorescence to iodide quenching. Overall, the results confirm previous kinetic analysis that demonstrated the hierarchical folding of Im7 via an on-pathway intermediate, and show that this species is a highly helical ensemble with a well-formed hydrophobic core. By contrast with the native state, however, the intermediate ensemble is flexible enough to change in response to mutation, its structural properties being tailored by residues in the sequence encompassing the native helix III.     

pH-dependent stability and folding kinetics of a protein with an unusual alpha-beta topology: the C-terminal domain of the ribosomal protein L9

The folding kinetics and thermodynamics of the isolated C-terminal domain of the ribosomal protein L9 (CTL9) have been studied as a function of pH. CTL9 is an alpha-beta protein that contains a single beta-sheet with an unusual mixed parallel, anti-parallel topology. The folding is fully reversible and two-state over the entire pH range. Stopped-flow fluorescence and CD experiments yield the same folding rate, and the chevron plots have the characteristic V-shape expected for two-state folding. The values of DeltaG*(H2O) and the m value calculated from the kinetic experiments are in excellent agreement with the equilibrium measurements. The extrapolated initial amplitudes of both the stopped-flow fluorescence and CD measurements show that there is no detectable burst phase intermediate. The domain contains three histidine residues, two of which are largely buried in the native state. They do not participate in salt-bridges or take part in a hydrogen bonded network. NMR measurements reveal that the buried histidine residues have significantly perturbed pK(a) values in the native state. The equilibrium stability and the folding rate are found to be strongly dependent upon their ionization state. There is a linear relationship between the log of the folding rate and DeltaG* (H2O) . The protein is much more stable and folds noticeably faster at pH values above the native state pK(a) values. DeltaG*(H2O) of unfolding increases from 2.90 kcal mol(-1) at pH 5.0 to 6.40 kcal mol(-1) at pH 8.0 while the folding rate increases from 0.60 to 18.7 s(-1). Tanford linkage analysis revealed that the interactions involving the two histidine residues are largely developed in the transition state. The results are compared to other studies of the pH-dependence of folding.     

A change in the apparent m value reveals a populated intermediate under equilibrium conditions in Escherichia coli ribonuclease HI

Experimental studies of protein stability often rely on the determination of an "m value", which describes the denaturant dependence of the free energy change between two states (DeltaG = DeltaG(H2O) - m[denaturant]). Changes in the m value accompanying site specific mutations are usually attributed to structural alterations in the native or unfolded ensemble. Here, we provide an example of significant reduction in the m value resulting from a subtle deviation in two-state behavior not detected by traditional methods. The protein that is studied is a variant of Escherchia coli RNase H in which three residues predicted to be involved in a partially buried salt bridge network were mutated to alanine (R46A, D102A, and D148A). Equilibrium denaturant profiles monitored by both fluorescence and circular dichroism appeared to be cooperative, and a two-state analysis yielded a DeltaG(UN) of approximately -3 kcal/mol with an m value of 1.4 kcal mol(-1) M(-1) (vs 2.3 for RNase H). Analysis of kinetic refolding experiments suggests that the system is actually three-state at equilibrium with an appreciable concentration of an intermediate state under low denaturant concentrations. The stability of the native state determined from a fit of these kinetic data is -6.7 kcal/mol, suggesting that the stability determined by traditional two-state equilibrium analysis is a gross underestimate. The only hint to this loss of two-state behavior was a decrease in the apparent m value, and the presence of the equilibrium intermediate was only identified by a kinetic analysis. Our work serves as a cautionary note; the possibility of a three-state system should be closely addressed before interpreting a change in the m value as a change in the native or unfolded state.     

Destabilization of the Escherichia coli RNase H kinetic intermediate: switching between a two-state and three-state folding mechanism

Escherichia coli RNase H folds through a partially folded kinetic intermediate that mirrors a rarely populated, partially unfolded form detectable by native-state hydrogen exchange under equilibrium conditions. Residue 53 is at the interface of two helices known to be structured in this intermediate. Kinetic refolding studies on mutant proteins varying in size and hydrophobicity at residue 53 support a contribution of hydrophobicity to the stabilities of the kinetic intermediate and the transition state. Packing interactions also play a significant role in the stability of these two states, though they play a much larger role in the native-state stability. One dramatic mutation, I53D, results in the conversion from a three-state to a two-state folding mechanism, which is explained most easily through a simple destabilization of the kinetic intermediate such that it is no longer stable with respect to the unfolded state. These results demonstrate that interactions that stabilize an intermediate can accelerate folding if these same interactions are present in the transition state. Our results are consistent with a hierarchical model of folding, where the intermediate consists of native-like interactions, is on-pathway, and is productive for folding.     

Role of the B helix in early folding events in apomyoglobin: evidence from site-directed mutagenesis for native-like long range interactions

The folding pathways of four mutants in which bulky hydrophobic residues in the B helix of apomyoglobin (ApoMb) are replaced by alanine (I28A, L29A, I30A, and L32A) have been analyzed using equilibrium and kinetic methods employing NMR, CD, fluorescence and mass spectrometry. Hydrogen exchange pulse-labeling followed by mass spectrometry reveals detectable intermediates in the kinetic folding pathways of each of these mutants. Comparison of the quench-flow data analyzed by NMR for the wild-type protein and the mutants showed that the substitutions I28A, L29A and L32A lead to destabilization of the B helix in the burst phase kinetic intermediate, relative to wild-type apomyoglobin. In contrast, the I30A mutation apparently has a slight stabilizing effect on the B helix in the burst phase intermediate; under weak labeling conditions, residues in the C helix region were also relatively stabilized in the mutant compared to the wild-type protein. This suggests that native-like helix B/helix C packing interactions occur in the folding intermediate. The L32A mutant showed significantly lower proton occupancies in the burst phase for several residues in the G helix, specifically F106, I107, E109 and A110, which are in close proximity to L32 in the X-ray structure of myoglobin, providing direct evidence that native-like helix B/helix G contacts are formed in the apomyoglobin burst phase intermediate. The L29A mutation resulted in an increase in burst phase proton occupancies for several residues in the E helix. Since these regions of the B and E helices are not in contact in the native myoglobin structure, these effects suggest the possibility of non-native B/E packing interactions in the kinetic intermediate. The differing effects of these B helix mutations on the apomyoglobin folding process suggests that each side-chain plays a different and important role in forming stable structure in the burst phase intermediate, and points to a role for both native-like and non-native contacts in stabilization of the folding intermediate.     

Folding subdomains of thioredoxin characterized by native-state hydrogen exchange

Native-state hydrogen exchange (HX) studies, used in conjunction with NMR spectroscopy, have been carried out on Escherichia coli thioredoxin (Trx) for characterizing two folding subdomains of the protein. The backbone amide protons of only the slowest-exchanging 24 amino acid residues, of a total of 108 amino acid residues, could be followed at pH 7. The free energy of the opening event that results in an amide hydrogen exchanging with solvent (DeltaG(op)) was determined at each of the 24 amide hydrogen sites. The values of DeltaG(op) for the amide hydrogens belonging to residues in the helices alpha(1), alpha(2), and alpha(4) are consistent with them exchanging with the solvent only when the fully unfolded state is sampled transiently under native conditions. The denaturant-dependences of the values of DeltaG(op) provide very little evidence that the protein samples partially unfolded forms, lower in energy than the unfolded state. The amide hydrogens belonging to the residues in the beta strands, which form the core of the protein, appear to have higher values of DeltaG(op) than amide hydrogens belonging to residues in the helices, suggesting that they might be more stable to exchange. This apparently higher stability to HX of the beta strands might be either because they exchange out their amide hydrogens in a high energy intermediate preceding the globally unfolded state, or, more likely, because they form residual structure in the globally unfolded state. In either case, the central beta strands-beta(3,) beta(2), and beta(4)-would appear to form a cooperatively folding subunit of the protein. The native-state HX methodology has made it possible to characterize the free energy landscape that Trx can sample under equilibrium native conditions.     

Involvement of cysteine residues and domain interactions in the reversible unfolding of lipoxygenase-1

Urea-induced unfolding of lipoxygenase-1 (LOX1) at pH 7.0 was followed by enzyme activity, spectroscopic measurements, and limited proteolysis experiments. Complete unfolding of LOX1 in 9 M urea in the presence of thiol reducing or thiol modifying reagents was observed. The aggregation and oxidative reactions prevented the reversible unfolding of the molecule. The loss of enzyme activity was much earlier than the structural loss of the molecule during the course of unfolding, with the midpoint concentrations being 4.5 and 7.0 M for activity and spectroscopic measurements, respectively. The equilibrium unfolding transition could be adequately fitted to a three-state, two-step model (N left arrow over right arrow I left arrow over right arrow U) and the intermediate fraction was maximally populated at 6.3 M urea. The free energy change (DeltaG(H(2)O)) for the unfolding of native (N) to intermediate (I) was 14.2 +/- 0.28 kcal/mol and for the intermediate to the unfolded state (U) was 11.9 +/- 0.12 kcal/mol. The ANS binding measurements as a function of urea concentration indicated that the maximum binding of ANS was in 6.3 M urea due to the exposure of hydrophobic groups; this intermediate showed significant amount of tertiary structure and retained nearly 60% of secondary structure. The limited proteolysis measurements showed that the initiation of unfolding was from the C-terminal domain. Thus, the stable intermediate observed could be the C-terminal domain unfolded with exposed hydrophobic domain-domain interface. Limited proteolysis experiments during refolding process suggested that the intermediate refolded prior to completely unfolded LOX1. These results confirmed the role of cysteine residues and domain-domain interactions in the reversible unfolding of LOX1. This is the first report of the reversible unfolding of a very large monomeric, multi-domain protein, which also has a prosthetic group.     

Kinetic traps in the folding/unfolding of procaspase-1 CARD domain

We have examined the folding and unfolding of the caspase recruitment domain of procaspase-1 (CP1-CARD), a member of the alpha-helical Greek key protein family. The equilibrium folding/unfolding of CP1-CARD is described by a two-state mechanism, and the results show CP1-CARD is marginally stable with a DeltaG(H2O) of 1.1 +/- 0.2 kcal/mole and an m-value of 0.65 +/- 0.06 kcal/mole/M (10 mM Tris-HCl at pH 8.0, 1 mM DTT, 25 degrees C). Consistent with the equilibrium folding data, CP1-CARD is a monomer in solution when examined by size exclusion chromatography. Single-mixing stopped-flow refolding and unfolding studies show that CP1-CARD folds and unfolds rapidly, with no detectable slow phases, and the reactions appear to reach equilibrium within 10 msec. However, double jump kinetic experiments demonstrate the presence of an unfolded-like intermediate during unfolding. The intermediate converts to the fully unfolded conformation with a half-time of 10 sec. Interrupted refolding studies demonstrate the presence of one or more nativelike intermediates during refolding, which convert to the native conformation with a half-time of about 60 sec. Overall, the data show that both unfolding and refolding processes are slow, and the pathways contain kinetically trapped species.     

NMR analysis of the conformational properties of the trapped on-pathway folding intermediate of the bacterial immunity protein Im7

Previous work shows that the transiently populated, on-pathway intermediate in Im7 folding contains three of the four native alpha-helices docked around a core stabilised by native and non-native interactions. To determine the structure and dynamic properties of this species in more detail, we have used protein engineering to trap the intermediate at equilibrium and analysed the resulting proteins using NMR spectroscopy and small angle X-ray scattering. Four variants were created. In L53AI54A, two hydrophobic residues within helix III are truncated, preventing helix III from docking stably onto the developing hydrophobic core. In two other variants, the six residues encompassing the native helix III were replaced with three (H3G3) or six (H3G6) glycine residues. In the fourth variant, YY, two native tyrosine residues (Tyr55 and Tyr56) were re-introduced into H3G6 to examine their role in determining the properties of the intermediate ensemble. All four variants show variable peak intensities and broad peak widths, consistent with these proteins being conformationally dynamic. Chemical shift analyses demonstrated that L53AI54A and YY contain native-like secondary structure in helices I and IV, while helix II is partly formed and helix III is absent. Lack of NOEs and rapid NH exchange for L53AI54A, combined with detailed analysis of the backbone dynamics, indicated that the hydrophobic core of this variant is not uniquely structured, but fluctuates on the NMR timescale. The results demonstrate that though much of the native-like secondary structure of Im7 is present in the variants, their hydrophobic cores remain relatively fluid. The comparison of H3G3/H3G6 and L53AI54A/YY suggests that Tyr55 and/or Tyr56 interact with the three-helix core, leading other residues in this region of the protein to dock with the core as folding progresses. In this respect, the three-helix bundle acts as a template for formation of helix III and the creation of the native fold.     

NMR characterisation of the relationship between frustration and the excited state of Im7

Previous work shows that Im9 folds in a two-state transition while its homologue Im7 folds in a three-state transition via an on-pathway kinetic intermediate state (KIS), with this difference being related to frustration in the structure of Im7. We have used NMR spectroscopy to study conformational dynamics connected to the frustration. A combination of equilibrium peptide N¹H/N²H exchange, model-free analyses of backbone NH relaxation data and relaxation dispersion (RD)-NMR shows that the native state of Im7 is in equilibrium with an intermediate state that is lowly populated [equilibrium intermediate state (EIS)]. Comparison of kinetic and thermodynamic parameters describing the EIS native-state equilibrium obtained by RD-NMR with previously reported parameters describing the KIS native-state equilibrium obtained from stopped-flow fluorescence studies of refolding His-tagged Im7 shows that the KIS and the EIS are the same species. ¹⁵N chemical shifts of the EIS obtained from the RD-NMR analysis show that residues forming helix III in the native state are unstructured in the EIS while other residues experiencing frustration in the native state are in structured regions of the EIS. We show that binding of Im7 and its L53A/I54A variant (which resembles the EIS as shown in previous work) to the cognate partner for Im7, the DNase domain of colicin E7, causes the dynamic processes associated with the frustration to be dampened.     

Thermodynamic and kinetic determinants of Thermotoga maritima cold shock protein stability: a structural and dynamic analysis

The cold shock protein (CSP) from hyperthermophile Thermotoga maritima (TmCSP) is only marginally stable (DeltaG(T(opt)) = 0.3 kcal/mol) at 353 K, the optimum environmental temperature (T(opt)) for T. maritima. In comparison, homologous CSPs from E. coli (DeltaG(T(opt)) = 2.2 kcal/mol) and B. subtilis (DeltaG(T(opt)) = 1.5 kcal/mol) are at least five times more stable at 310 K, the T(opt) for the mesophiles. Yet at the room temperature, TmCSP is more stable (DeltaG(T(R)) = 4.7 kcal/mol) than its homologues (DeltaG(T(R)) = 3.0 kcal/mol for E. coli CSP and DeltaG(T(R)) = 2.1 kcal/mol for B. subtilis CSP). This unique observation suggests that kinetic, rather than thermodynamic, barriers toward unfolding might help TmCSP native structure at high temperatures. Consistently, the unfolding rate of TmCSP is considerably slower than its homologues. High temperature (600 K) complete unfolding molecular dynamics (MD) simulations of TmCSP support our hypothesis and reveal an unfolding scheme unique to TmCSP. For all the studied homologues of TmCSP, the unfolding process first starts at the C-terminal region and N-terminal region unfolds in the end. But for TmCSP, both the terminals resist unfolding for consistently longer simulation times and, in the end, unfold simultaneously. In TmCSP, the C-terminal region is better fortified and has better interactions with the N-terminal region due to the charged residues, R2, E47, E49, H61, K63, and E66, being in spatial vicinity. The electrostatic interactions among these residues are unique to TmCSP. Consistently, the room temperature MD simulations show that TmCSP is more rigid at its N- and C-termini as compared to its homologues from E. coli, B. subtilis, and B. caldolyticus.     

Structural role of the proline residues of the beta-hinge region of p13suc1 as revealed by site-directed mutagenesis and fluorescence studies.

Site-directed mutagenesis, gel filtration, and fluorescence spectroscopy approaches were used to study the molecular hinge mechanism involved in the beta-strand-exchanged dimer formation of the cyclin-dependent protein kinase regulatory subunit p13(suc1) from Schizosaccharomyces pombe. Single and double mutants of residues Pro-90 and Pro-92 (P90V, P92V, and P90V/P92V) were prepared and assayed. Substitution of Pro-90 prevented dimer formation by arm exchange. However, single point mutations did not affect the two-state unfolding transition of wild-type p13(suc1) at equilibrium (i.e., wild type, DeltaG degrees (0,un) = 7.38 +/- 0.35 kcal mol(-1), vs P90V, DeltaG degrees (0,un) = 6.71 +/- 0.18 kcal mol(-1)). On the contrary, the double mutant unfolded with a complex transition, and the reaction was best described by a three-state model (N <==> I <==> U). Resolution of the state-dependent (native vs denatured) intrinsic fluorescence decay amplitudes of p13(suc1) showed that with P90V/P92V these parameters were affected at [GuHCl] significantly less than with wild-type and single mutant proteins. Moreover, with the latter products, fluorescence quenching measurements at 1 M GuHCl revealed linear Stern-Volmer plots with quenching constants typical of tryptophan residues located in a native environment (1.6 M(-1) < K(SV) < 2.3 M(-1)). Dissimilarly, with P90V/P92V a significant deviation from linearity of the Stern-Volmer plot was obtained. Nonlinear least-squares analysis of these data resolved the significant contribution of highly solvent-accessible emitting species (K(SV) = 26 M(-1)) consistent with large exposure of the tryptophan residues. These results are compatible with the existence of an intermediate unfolding state of the double mutation product. Thus, while single residue substitution studies give support to the primary role of Pro-90 in the p13(suc1) dimer formation by domain swapping, double residue substitution studies indicate the important role of the conserved repeat, Pro-x-Pro, for the proper beta-strand spatial organization and stability.     

Apo-azurin folds via an intermediate that resembles the molten-globule

The folding of Pseudomonas aeruginosa apo-azurin was investigated with the intent of identifying putative intermediates. Two apo-mutants were constructed by replacing the main metal-binding ligand C112 with a serine (C112S) and an alanine (C112A). The guanidinium-induced unfolding free energies (DeltaG(U-N)(H2O)) of the C112S and C112A mutants were measured to 36.8 +/- 1 kJ mole(-1) and 26.1 +/- 1 kJ mole(-1), respectively, and the m-value of the transition to 23.5 +/- 0.7 kJ mole(-1) M(-1). The difference in folding free energy (DeltaDeltaG(U-N)(H2O)) is largely attributed to the intramolecular hydrogen bonding properties of the serine Ogamma in the C112S mutant, which is lacking in the C112A structure. Furthermore, only the unfolding rates differ between the two mutants, thus pointing to the energy of the native state as the source of the observed Delta DeltaG(U-N)(H2O). This also indicates that the formation of the hydrogen bonds present in C112S but absent in C112A is a late event in the folding of the apo-protein, thus suggesting that formation of the metal-binding site occurs after the rate-limiting formation of the transition state. In both mutants we also noted a burst-phase intermediate. Because this intermediate was capable of binding 1-anilinonaphtalene-8-sulfonate (ANS), as were an acid-induced species at pH 2.6, we ascribe it molten globule-like status. However, despite the presence of an intermediate, the folding of apo-azurin C112S is well approximated by a two-state kinetic mechanism.     

Investigation of the folding pathway of the TEM-1 β-lactamase

The TEM-1 β-lactamase is a globular protein containing 12 proline residues. The folding mechanism of this enzyme was investigated by kinetic and equilibrium experiments with the help of fluorescence spectroscopy and circular dichroism. The equilibrium denaturation of the protein induced by guanidine hydrochloride occurs in two discrete steps, indicating the existence of a thermodynamically stable intermediate state. Thisstate is 5.2 ± 0.4 kcal/mol less stable than the native conformation and 5.7 ± 0.2 kcal/mol more stable than the fully denaturedprotein. This intermediate state exhibits a high content of native secondary structure elements but is devoid of specific tertiary organization; its relation to the “molten globule” is discussed. Refolding kinetic experimentsrevealed the existence of a transient intermediate conformation between thethermodynamically stable intermediate and the native protein. This transient intermediate appears rapidly during the folding reaction. It exhibits a secondary structure content very similar to that of the native protein and has also recovered a significant amount of tertiary organisation. The final refolding step of the TEM-1 β-lactamase, leading to the native enzyme, is dominated by two major slow kinetic phases which probablyreflect a very complex process kinetically limited by proline cis/transisomerization. © 1995 Wiley-Liss, Inc.     

Increase in the conformational flexibility of beta 2-microglobulin upon copper binding: a possible role for copper in dialysis-related amyloidosis

A key pathological event in dialysis-related amyloidosis is the fibril formation of beta(2)-microglobulin (beta 2-m). Because beta 2-m does not form fibrils in vitro, except under acidic conditions, predisposing factors that may drive fibril formation at physiological pH have been the focus of much attention. One factor that may be implicated is Cu(2+) binding, which destabilizes the native state of beta 2-m and thus stabilizes the amyloid precursor. To address the Cu(2+)-induced destabilization of beta 2-m at the atomic level, we studied changes in the conformational dynamics of beta 2-m upon Cu(2+) binding. Titration of beta 2-m with Cu(2+) monitored by heteronuclear NMR showed that three out of four histidines (His13, His31, and His51) are involved in the binding at pH 7.0. (1)H-(15)N heteronuclear NOE suggested increased backbone dynamics for the residues Val49 to Ser55, implying that the Cu(2+) binding at His51 increased the local dynamics of beta-strand D. Hydrogen/deuterium exchange of amide protons showed increased flexibility of the core residues upon Cu(2+) binding. Taken together, it is likely that Cu(2+) binding increases the pico- to nanosecond fluctuation of the beta-strand D on which His51 exists, which is propagated to the core of the molecule, thus promoting the global and slow fluctuations. This may contribute to the overall destabilization of the molecule, increasing the equilibrium population of the amyloidogenic intermediate.     

Thermodynamics of single peptide bond cleavage in bovine pancreatic trypsin inhibitor (BPTI)

A major goal of this paper was to estimate a dynamic range of equilibrium constant for the opening of a single peptide bond in a model protein, bovine pancreatic trypsin inhibitor (BPTI). Ten mutants of BPTI containing a single Xaa-->Met substitution introduced in different parts of the molecule were expressed in Escherichia coli. The mutants were folded, purified to homogeneity, and cleaved with cyanogen bromide to respective cleaved forms. Conformation of the intact mutants was similar to the wildtype, as judged from their circular dichroism spectra. Substantial conformational changes were observed on the chemical cleavage of three single peptide bonds--Met46-Ser, Met49-Cys, and Met53-Thr--located within the C-terminal helix. Cleavage of those peptide bonds caused a significant destabilization of the molecule, with a drop of the denaturation temperature by 56.4 degrees C to 68 degrees C at pH 4.3. Opening of the remaining seven peptide bonds was related to a 10.8 degrees C to 39.4 degrees C decrease in T(den). Free energies of the opening of 10 single peptide bonds in native mutants (Delta G(op,N)) were estimated from the thermodynamic cycle that links denaturation and cleavage free energies. To calculate those values, we assumed that the free energy of opening of a single peptide bond in the denatured state (Delta G(op,D)) was equal to -2.7 kcal/mole, as reported previously. Calculated Delta G(op,N) values in BPTI were in the range from 0.2 to 10 kcal/mole, which was equivalent to a >1 million-fold difference in equilibrium constants. The values of Delta G(op,N) were the largest for peptide bonds located in the C-terminal helix and significantly lower for peptide bonds in the beta-structure or loop regions. It appears that opening constants for single peptide bonds in various proteins span across 33 orders of magnitude. Typical equilibrium values for a single peptide bond opening in a protein containing secondary structure elements fall into negligibly low values, from 10(-3) to 10(-8), and are efficient to ensure stability against proteolysis.     

GuHCl and NaCl-dependent hydrogen exchange in MerP reveals a well-defined core with an unusual exchange pattern

We have analysed hydrogen exchange at amide groups to characterise the energy landscape of the 72 amino acid residue protein MerP. From the guanidine hydrochloride (GuHCl) dependence of exchange in the pre-transitional region we have determined free energy values of exchange (DeltaG(HX)) and corresponding m-values for individual amide protons. Detailed analysis of the exchange patterns indicates that for one set of amide protons there is a weak dependence on denaturant, indicating that the exchange is dominated by local fluctuations. For another set of amide protons a linear, but much stronger, denaturant dependence is observed. Notably, the plots of free energy of exchange versus [GuHCl] for 16 amide protons show pronounced upward curvature, and a close inspection of the structure shows that these residues form a well-defined core in the protein. The hydrogen exchange that was measured at various concentrations of NaCl shows an apparent selective stabilisation of this core. Detailed analysis of this exchange pattern indicates that it may originate from selective destabilisation of the unfolded state by guanidinium ions and/or selective stabilisation of the core in the native state by chloride ions.     

Identification of native and non-native structure in kinetic folding intermediates of apomyoglobin

Site-directed mutagenesis has been used to probe the interactions that stabilize the equilibrium and burst phase kinetic intermediates formed by apomyoglobin. Nine bulky hydrophobic residues in the A, E, G and H helices were replaced by alanine, and the effects on protein stability and kinetic folding pathways were determined. Hydrogen exchange pulse-labeling experiments, with NMR detection, were performed for all mutants. All of the alanine substitutions resulted in changes in proton occupancy or an increased rate of hydrogen-deuterium exchange for amides in the immediate vicinity of the mutation. In addition, most mutations affected residues in distant parts of the amino acid sequence, providing insights into the topology of the burst phase intermediate and the interactions that stabilize its structure. Differences between the pH 4 equilibrium molten globule and the kinetic intermediate are evident: the E helix region plays no discernible role in the equilibrium intermediate, but contributes significantly to stabilization of the ensemble of compact intermediates formed during kinetic refolding. Mutations that interfere with docking of the E helix onto the preformed A/B/G/H helix core substantially decrease the folding rate, indicating that docking and folding of the E helix region occurs prior to formation of the apomyoglobin folding transition state. The results of the mutagenesis experiments are consistent with rapid formation of an ensemble of compact burst phase intermediates with an overall native-like topological arrangement of the A, B, E, G, and H helices. However, the experiments also point to disorder in docking of the E helix and to non-native contacts in the kinetic intermediate. In particular, there is evidence for translocation of the H helix by approximately one helical turn towards its N terminus to maximize hydrophobic interactions with helix G. Thus, the burst phase intermediate observed during kinetic refolding of apomyoglobin consists of an ensemble of compact, kinetically trapped states in which the helix docking appears to be topologically correct, but in which there are local non-native interactions that must be resolved before the protein can fold to the native structure.     

Interleukin-1 beta folding between pH 5 and 7: experimental evidence for three-state folding behavior and robust transition state positions late in folding

The thermodynamic stability and folding kinetics of the all beta-sheet protein interleukin-1beta were measured between 0 and 4 M GdmCl concentrations and pH 5-7. Native interleukin-1beta undergoes a 3.5 kcal/mol decrease in thermodynamic stability, Delta, as pH is increased from 5 to 7. The native state parameter m(NU), measuring protein destabilization/[GdmCl], remains constant between pH 5 and 7, indicating that the solvent-exposed surface area difference between the native state and unfolded ensemble is unchanged across this pH range. Similarly, pH changes between 5 and 7 decrease only the thermodynamic stability, DeltaG(H)2(O), and not the m-values, of the kinetic intermediate and transition states. This finding is shown to be consistent with transition state configurations which continue to be the high-energy configurations of the transition state in the face of changing stability conditions. A three-state folding mechanism U right arrow over left arrow I right arrow over left arrow N is shown to be sufficient in characterizing IL-1beta folding under all conditions studied. The m-values of refolding transitions are much larger than the m-values of unfolding transitions, indicating that that the fast, T(2) (U right arrow over left arrow I), and slow, T(1) (I right arrow over left arrow N), transition states are highly similar to the intermediate I and native state N, respectively. Many of the folding properties of interleukin-1beta are shared among other members of the beta-trefoil protein family, although clear differences can exist.     

High-pressure denaturation of apomyoglobin

The pressure denaturation of wild type and mutant apomyoglobin (apoMb) was investigated using a high-pressure, high-resolution nuclear magnetic resonance and high-pressure fluorescence techniques. Wild type apoMb is resistant to pressures up to 80 MPa, and denatures to a high-pressure intermediate, I(p), between 80 and 200 MPa. A further increase of pressure to 500 MPa results in denaturation of the intermediate. The two tryptophans, both in the A helix, remain sequestered from solvent in the high-pressure intermediate, which retains some native NOESY cross peaks in the AGH core as well as between F33 and F43. High-pressure fluorescence shows that the tryptophans remain inaccessible to solvent in the I(p) state. Thus the high-pressure intermediate has some structural properties in common with the apoMb I(2) acid intermediate. The resistance of the AGH core to pressures up to 200 MPa provides further evidence that the intrinsic stability of these alpha-helices is responsible for their presence in a number of equilibrium intermediates as well as in the earliest kinetic folding intermediate. Mutations in the AGH core designed to disrupt packing by burying a charge or increasing the size of a hydrophobic residue significantly perturbed the unfolding of native apoMb to the high-pressure intermediate. The F123W and S108L mutants both unfolded at lower pressures, while retaining some resistance to pressures below 50 MPa. The charge burial mutants, A130K and S108K, are not stable at very low pressures and both denature to the intermediate by 100 MPa, half of the pressure required for wild type apoMb. Thus a similar intermediate state is created independent of the method of perturbation, and mutations have similar effects on native state destabilization for both methods of denaturation. These data suggest that equilibrium intermediates that can be formed through different means are likely to resemble a kinetic intermediate.