Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

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.     

The apomyoglobin folding pathway revisited: structural heterogeneity in the kinetic burst phase intermediate

Extensive analysis of accurate quench-flow hydrogen exchange results indicates that the burst phase kinetic intermediate in the folding of apomyoglobin (apoMb) from urea is structurally heterogeneous. The structural variability is associated with the partial folding of the E helix during the burst phase (<6.4ms) of the folding process. Analysis of the effects of exchange-out of amide proton labels during the labeling pulse ( approximately pH 10) of the quench-flow process indicates that three of the amide protons in the E helix are in fact largely protected in the burst phase of folding, while the remainder of the E helix has a substantial complement of amide protons that show biphasic kinetics, i.e. are protected partly during the burst phase and partly during the slow phase of folding. The locations of these amide protons can be used to map the sites of structural heterogeneity in the kinetic molten globule. These sites include, besides the E helix, the ends of the A and B helices and part of the C helix. Our results give significant support to the hypothesis that the kinetic molten globule intermediate of apoMb is native-like.     

Population analyses of kinetic partitioning in protein folding

A simple lattice model of protein folding is studied in order to analyze the kinetic partitioning phenomena in the energy landscape perspective. By restricting the area of conformational space, it becomes possible to follow many Monte Carlo trajectories until they reach equilibrium. Alteration of population of trajectories is monitored and the relations between the energy landscape and kinetics are examined. Kinetic partitioning phenomena are categorized into different types in terms of characteristic time constants and partitioning ratio. In a specific partitioning process, refolding proceeds along the parallel pathways; the time constants have a temperature dependence similar to that observed in hen lysozyme. High-energy conformations are classified into groups according to the probability that the trajectories starting from those conformations will reach each energy valley. The partitioning ratio is determined by the way in which the conformational space is organized into these groups.     

Early intermediates in the folding of dihydrofolate reductase from Escherichia coli detected by hydrogen exchange and NMR.

The kinetic folding mechanism for Escherichia coli dihydrofolate reductase postulates two distinct types of transient intermediates. The first forms within 5 ms and has substantial secondary structure but little stability. The second is a set of four species that appear over the course of several hundred milliseconds and have secondary structure, specific tertiary structure, and significant stability (Jennings PA, Finn BE, Jones BE, Matthews CR, 1993, Biochemistry 32:3783-3789). Pulse labeling hydrogen exchange experiments were performed to determine the specific amide hydrogens in alpha-helices and beta-strands that become protected from exchange through the formation of stable hydrogen bonds during this time period. A significant degree of protection was observed for two subsets of the amide hydrogens within the dead time of this experiment (6 ms). The side chains of one subset form a continuous nonpolar strip linking six of the eight strands in the beta-sheet. The other subset corresponds to a nonpolar cluster on the opposite face of the sheet and links three of the strands and two alpha-helices. Taken together, these data demonstrate that the complex strand topology of this eight-stranded sheet can be formed correctly within 6 ms. Measurement of the protection factors at three different folding times (13 ms, 141 ms, and 500 ms) indicates that, of the 13 amide hydrogens displaying significant protection within 6 ms, 8 exhibit an increase in their protection factors from approximately 5 to approximately 50 over this time range; the remaining five exhibit protection factors > 100 at 13 ms. Only approximately half of the population of molecules form this set of stable hydrogen bonds. Thirteen additional hydrogens in the beta-sheet become protected from exchange as the set of native conformers appear, suggesting that the stabilization of this network reflects the global cooperativity of the folding reaction.     

Insights into conformation and dynamics of protein GB1 during folding and unfolding by NMR

Understanding protein stability requires characterization of structural determinants of the folded and unfolded states. Many proteins are capable of populating partially folded states under specific solution conditions. Occasionally, coexistence of the folded and an unfolded state under non- or mildly denaturing conditions can be observed by NMR, allowing us to structurally probe these states under identical conditions. Here we report on a destabilized mutant of the B1 domain of protein G (GB1) whose equilibrium unfolding was systematically investigated. Backbone amide residual dipolar couplings (RDCs), the tryptophan Nepsilon-H resonance and the amide nitrogen transverse relaxation rates (R2s) for varying pH values and different temperatures were measured. The backbone amide RDCs indicate that prior to complete unfolding, two melting hot spots are formed at the turn around T11, L12 and K13 and the N terminus of the helix at A24 and T25. The RDCs for the low pH, thermally unfolded state of GB1 are very small and do not indicate the presence of any native-like structure. Amide nitrogen transverse relaxation rates for GB1 in the folded state at different temperatures exhibit large contributions from exchange processes and the associated dynamics display considerable heterogeneity. Our data provide clear evidence for intermediate conformations and multi-state equilibrium un/folding for this GB1 variant.     

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.     

Time-resolved small-angle X-ray scattering study of the folding dynamics of barnase

Structural changes of barnase during folding were investigated using time-resolved small-angle X-ray scattering (SAXS). The folding of barnase involves a burst-phase intermediate, sometimes designated as the denatured state under physiological conditions, D_phys, and a second hidden intermediate. Equilibrium SAXS measurements showed that the radius of gyration (R_g) of the guanidine unfolded state (U) is 26.9 ± 0.7 Å, which remains largely constant over a wide denaturant concentration range. Time-resolved SAXS measurements showed that the R_g value extrapolated from kinetic R_g data to time zero, R_g,0, is 24.3 ± 0.1 Å, which is smaller than that of U but which is expanded from that of folding intermediates of other proteins with similar chain lengths (19 Å). After the burst-phase change, a single-exponential reduction in R_g² was observed, which corresponds to the formation of the native state for the major component containing the native trans proline isomer. We estimated R_g of the minor component of D_phys containing the non-native cis proline isomer (D_phys,cis) to be 25.7 ± 0.6 Å. Moreover, R_g of the major component of D_phys containing the native proline isomer (D_phys,tra) was estimated as 23.9 ± 0.2 Å based on R_g,0. Consequently, both components of the burst-phase intermediate of barnase (D_phys,tra and D_phys,cis) are still largely expanded. It was inferred that D_phys possesses the N-terminal helix and the center of the β-sheet formed independently and that the formation of the remainder of the protein occurs in the slower phase.     

Interaction of SecB with intermediates along the folding pathway of maltose-binding protein.

SecB, a molecular chaperone involved in protein export in Escherichia coli, displays the remarkable ability to selectively bind many different polypeptide ligands whose only common feature is that of being nonnative. The selectivity is explained in part by a kinetic partitioning between the folding of a polypeptide and its association with SecB. SecB has no affinity for native, stably folded polypeptides but interacts tightly with polypeptides that are nonnative. In order to better understand the nature of the binding, we have examined the interaction of SecB with intermediates along the folding pathway of maltose-binding protein. Taking advantage of forms of maltose-binding protein that are altered in their folding properties, we show that the first intermediate in folding, represented by the collapsed state, binds to SecB, and that the polypeptide remains active as a ligand until it crosses the final energy barrier to attain the native state.     

Structure and heterogeneity of the one- and two-disulfide folding intermediates of tick anticoagulant peptide

Tick anticoagulant peptide (TAP) is a factor Xa-specific inhibitor and is structurally homologous to bovine pancreatic trypsin inhibitor (BPTI). The fully reduced TAP refolds spontaneously to form the native structure under a wide variation of redox buffers. The folding intermediates of TAP consist of at least 22 fractions of one-disulfide, two-disulfide, and three-disulfide scrambled isomers. Three species of well-populated one- and two-disulfide intermediates were isolated and structurally characterized. The predominant one-disulfide species contains TAP-(Cys33—Cys55). Two major two-disulfide isomers were TAP-(Cys33—Cys55, Cys15—Cys39) and TAP-(Cys33—Cys55, Cys5—Cys39). Both Cys33—Cys55 and Cys15—Cys39 are native disulfides of TAP. These three species are structural counterparts of BPTI-(Cys30—Cys51), BPTI-(Cys30—Cys51, Cys14—Cys38), and BPTI-(Cys30—Cys51,Cys5—Cys38), which have been shown to be the major intermediates of BPTI folding. In addition, time-course-trapped folding intermediates of TAP, consisting of about 47% one-disulfide species and 30% two-disulfide species, were collectively digested with thermolysin, and fragmented peptides were analyzed by Edman sequencing and mass spectrometry in order to characterize the disulfide-containing peptides. Among the 15 possible single-disulfide pairings of TAP, 10 (2 native and 8 nonnative) were found as structural components of its one- and two-disulfide folding intermediates. The results demonstrate that the major folding intermediates of TAP bear structural homology to those of BPTI. However, the folding pathway of TAP differs from that of BPTI by (a) a higher degree of heterogeneity of one- and two-disulfide intermediates and (b) the presence of three-disulfide scrambled isomers as folding intermediates. Mechanism(s) that may account for these diversities are proposed and discussed.     

NMR and protein folding: equilibrium and stopped-flow studies.

NMR studies are now unraveling the structure of intermediates of protein folding using hydrogen-deuterium exchange methodologies. These studies provide information about the time dependence of formation of secondary structure. They require the ability to assign specific resonances in the NMR spectra to specific amide protons of a protein followed by experiments involving competition between folding and exchange reactions. Another approach is to use 19F-substituted amino acids to follow changes in side-chain environment upon folding. Current techniques of molecular biology allow assignments of 19F resonances to specific amino acids by site-directed mutagenesis. It is possible to follow changes and to analyze results from 19F spectra in real time using a stopped-flow device incorporated into the NMR spectrometer.     

Intermediates and the folding of proteins L and G

We use a minimalist protein model, in combination with a sequence design strategy, to determine differences in primary structure for proteins L and G, which are responsible for the two proteins folding through distinctly different folding mechanisms. We find that the folding of proteins L and G are consistent with a nucleation-condensation mechanism, each of which is described as helix-assisted beta-1 and beta-2 hairpin formation, respectively. We determine that the model for protein G exhibits an early intermediate that precedes the rate-limiting barrier of folding, and which draws together misaligned secondary structure elements that are stabilized by hydrophobic core contacts involving the third beta-strand, and presages the later transition state in which the correct strand alignment of these same secondary structure elements is restored. Finally, the validity of the targeted intermediate ensemble for protein G was analyzed by fitting the kinetic data to a two-step first-order reversible reaction, proving that protein G folding involves an on-pathway early intermediate, and should be populated and therefore observable by experiment.     

ALASKA: A Mathematica package for two-state kinetic analysis of protein folding reactions

A Mathematica package (ALASKA) has been developed to simplify the measurement of protein folding kinetics by analysis of ¹H NMR lineshape analysis. This package reads NMR data in ASCII format and can simulate an aromatic ¹ NMR spectrum with or without lineshape broadening from chemical exchange. We describe the analysis of a urea denaturation series of a fast-folding protein, the G46A/G48A variant of monomeric repressor.     

Fast folding of a prototypic polypeptide: the immunoglobulin binding domain of streptococcal protein G.

The folding of the small (56 residues) highly stable B1 immunoglobulin binding domain (GB1) of streptococcal protein G has been investigated by quenched-flow deuterium-hydrogen exchange. This system represents a paradigm for the study of protein folding because it exhibits no complicating features superimposed upon the intrinsic properties of the polypeptide chain. Collapse to a semicompact state exhibiting partial order, reflected in protection factors for ND-NH exchange up to 10-fold higher than that expected for a random coil, occurs within the dead time (< or = 1 ms) of the quenched flow apparatus. This is followed by the formation of the fully native state, as monitored by the fractional proton occupancy of 26 backbone amide groups spread throughout the protein, in a single rapid concerted step with a half-life of 5.2 ms at 5 degrees C.     

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.     

The role of Trp-82 in the folding of intestinal fatty acid binding protein

Multiple phases have been observed during the folding and unfolding of intestinal fatty acid binding protein (WT-IFABP) by stopped-flow fluorescence. Site-directed mutagenesis has been used to examine the role of each of the two tryptophans of this protein in these processes. The unfolding and refolding kinetics of the mutant protein containing only tryptophan 82 (W6Y-IFABP) showed that the tryptophan at this location was critical to the fluorescence signal changes observed throughout the unfolding reaction and early in the refolding reaction. However, the kinetic patterns of the mutant protein containing only tryptophan 6 (W82Y-IFABP) indicated that the tryptophan at this location participated in the fluorescence signal changes observed early in the unfolding reaction and late in the refolding reaction. Together, these data suggest that native-like structure was formed first in the vicinity of tryptophan 82, near the center of the hydrophobic core of this beta-sheet protein, prior to formation of native-like structure in the periphery of the protein.     

A study of intermediates involved in the folding pathway for recombinant human macrophage colony-stimulating factor (M-CSF): evidence for two distinct folding pathways.

The folding pathway for a 150-amino acid recombinant form of the dimeric cytokine human macrophage colony-stimulating factor (M-CSF) has been studied. All 14 cysteine residues in the biologically active homodimer are involved in disulfide linkages. The structural characteristics of folding intermediates blocked with iodoacetamide reveal a rapid formation of a small amount of a non-native dimeric intermediate species followed by a slow progression via both monomeric and dimeric intermediates to the native dimer. The transition from monomer to fully folded dimer is complete within 25 h at room temperature at pH 9.0. The blocked intermediates are stable under conditions of sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and thus represent various dimeric and folded monomeric species of the protein with different numbers of disulfide bridges. Peptide mapping and electrospray ionization mass spectrometry revealed that a folded monomeric species of M-CSF contained three of the four native disulfide bridges, and this folded monomer also showed some biological activity in a cell-based assay. The results presented here strongly suggest that M-CSF can fold via two different pathways, one involving monomeric intermediates and another involving only dimeric intermediates.     

Molecular basis of cooperativity in protein folding. V. Thermodynamic and structural conditions for the stabilization of compact denatured states

The heat-denatured state of proteins has been usually assumed to be a fully hydrated random coil. It is now evident that under certain solvent conditions or after chemical or genetic modifications, the protein molecule may exhibit a hydrophobic core and residual secondary structure after thermal denaturation. This state of the protein has been called the “compact denatured” or “molten globule” state. Recently is has been shown that α-lactalbumin at pH < 5 denatures into a molten globule state upon increasing the temperature (Griko, Y., Freire, E., Privalov, P. L. Biochemistry 33:1889–1899, 1994). This state has a lower heat capacity and a higher enthalpy at low temperatures than the unfolded state. At those temperatures the stabilization of the molten globule state is of an entropic origin since the enthalpy contributes unfavorably to the Gibbs free energy. Since the molten globule is more structured than the unfolded state and, therefore, is expected to have a lower configurational entropy, the net entropic gain must originate primarily from solvent related entropy arising from the hydrophobic effect, and to a lesser extent from protonation or electrostatic effects. In this work, we have examined a large ensemble of partly folded states derived from the native structure of α-lactalbumin in order to identify those states that satisfy the energetic criteria of the molten globule. It was found that only few states satisfied the experimental constraints and that, furthermore, those states were part of the same structural family. In particular, the regions corresponding to the A, B, and C helices were found to be folded, while the β sheet and the D helix were found to be unfolded. At temperatures below 45°C the states exhibiting those structural characteristics are enthalpically higher than the unfolded state in agreement with the experimental data. Interestingly, those states have a heat capacity close to that observed for the acid pH compact denatured state of α-lactalbumin [980 cal (mol.K)^?l]. In addition, the folded regions of these states include those residues found to be highly protected by NMR hydrogen exchange experiments. This work represents an initial attempt to model the structural origin of the thermodynamic properties of partly folded states. The results suggest a number of structural features that are consistent with experimental data. © 1994 Wiley-Liss, Inc.     

Rapid collapse precedes the fast two-state folding of the cold shock protein

The cold shock protein Bc-Csp folds very rapidly in a reaction that is well described by a kinetic two-state mechanism without intermediates. We measured the shortening of six intra-protein distances during folding by F?rster resonance energy transfer (FRET) in combination with stopped-flow experiments. Single tryptophan residues were engineered into the protein as the donors, and single 5-(((acetylamino)ethyl)amino)naphthalene-1-sulfonate (AEDANS) residues were placed as the acceptors at solvent-exposed sites of Bc-Csp. Their R0 value of about 22 A was well suited for following distance changes during the folding of this protein with a high sensitivity. The mutagenesis and the labeling did not alter the refolding kinetics. The changes in energy transfer during folding were monitored by both donor and acceptor emission and reciprocal effects were found. In two cases the donor-acceptor distances were similar in the unfolded and the folded state and, as a consequence, the kinetic changes in energy transfer upon folding were very small. For four donor/acceptor pairs we found that > or =50% of the increase in energy transfer upon folding occurred prior to the rate-limiting step of folding. This reveals that about half of the shortening of the intra-molecular distances upon folding has occurred already before the rate-limiting step and suggests that the fast two-state folding reaction of Bc-Csp is preceded by a very rapid collapse.