Kinetic consequences of native state optimization of surface-exposed electrostatic interactions in the Fyn SH3 domain

Optimization of surface exposed charge-charge interactions in the native state has emerged as an effective means to enhance protein stability; but the effect of electrostatic interactions on the kinetics of protein folding is not well understood. To investigate the kinetic consequences of surface charge optimization, we characterized the folding kinetics of a Fyn SH3 domain variant containing five amino acid substitutions that was computationally designed to optimize surface charge-charge interactions. Our results demonstrate that this optimized Fyn SH3 domain is stabilized primarily through an eight-fold acceleration in the folding rate. Analyses of the constituent single amino acid substitutions indicate that the effects of optimization of charge-charge interactions on folding rate are additive. This is in contrast to the trend seen in folded state stability, and suggests that electrostatic interactions are less specific in the transition state compared to the folded state. Simulations of the transition state using a coarse-grained chain model show that native electrostatic contacts are weakly formed, thereby making the transition state conducive to nonspecific, or even nonnative, electrostatic interactions. Because folding from the unfolded state to the folding transition state for small proteins is accompanied by an increase in charge density, nonspecific electrostatic interactions, that is, generic charge density effects can have a significant contribution to the kinetics of protein folding. Thus, the interpretation of the effects of amino acid substitutions at surface charged positions may be complicated and consideration of only native-state interactions may fail to provide an adequate picture.     

Propagation of a single destabilizing mutation throughout the Escherichia coli ribonuclease HI native state

A point mutation (I53A) in the core of Escherichia coli RNase H* is known to destabilize both the native conformation (DeltaG(UN)) and the kinetic intermediate (DeltaG(UI)) by 2 kcal/mole. Here, we have used native-state hydrogen deuterium exchange to ask how this destabilization is propagated throughout the molecule. Stability parameters were obtained for individual residues in I53A and compared with those from the wild-type protein. A destabilization of 2 kcal/mole was observed in residues in the core but was not detected in the periphery of the molecule. These results are consistent with the localized destabilization of the core observed in the early intermediate of the kinetic folding pathway, supporting the resemblance of this kinetic intermediate to the partially unfolded form detected in the native state at equilibrium. A thermodynamic cycle also shows no interaction between Ile 53 and a residue in the periphery. There is, however, an increase in the number of denaturant-independent exchange events in the periphery of I53A, showing that effects of the point mutation are communicated to regions outside the core, although perhaps not through changes in stability. In sum, this work shows that localized regions within a protein can be destabilized independently. Furthermore, it implies a correspondence between the kinetic intermediate and the equilibrium PUF, as the magnitude and localization of the destabilization are the same in both.     

Folding/Unfolding Kinetics of Apomyoglobin

The equilibrium and kinetic folding/unfolding of apomyoglobin (ApoMb) were studied at pH 6.2, 11 °C by recording tryptophan fluorescence. The equilibrium unfolding of ApoMb in the presence of urea was shown to involve accumulation of an intermediate state, which had a higher fluorescence intensity as compared with the native and unfolded states. The folding proceeded through two kinetic phases, a rapid transition from the unfolded to the intermediate state and a slow transition from the intermediate to the native state. The accumulation of the kinetic intermediate state was observed in a wide range of urea concentrations. The intermediate was detected even in the region corresponding to the unfolding limb of the chevron plot. Urea concentration dependence was obtained for the observed folding/unfolding rate. The shape of the dependence was compared with that of two-state proteins characterized by a direct transition from the unfolded to the native state.     

Side chain burial and hydrophobic core packing in protein folding transition states

A critical step in the folding pathway of globular proteins is the formation of a tightly packed hydrophobic core. Several mutational studies have addressed the question of whether tight packing interactions are present during the rate-limiting step of folding. In some of these investigations, substituted side chains have been assumed to form native-like interactions in the transition state when the folding rates of mutant proteins correlate with their native-state stabilities. Alternatively, it has been argued that side chains participate in nonspecific hydrophobic collapse when the folding rates of mutant proteins correlate with side-chain hydrophobicity. In a reanalysis of published data, we have found that folding rates often correlate similarly well, or poorly, with both native-state stability and side-chain hydrophobicity, and it is therefore not possible to select an appropriate transition state model based on these one-parameter correlations. We show that this ambiguity can be resolved using a two-parameter model in which side chain burial and the formation of all other native-like interactions can occur asynchronously. Notably, the model agrees well with experimental data, even for positions where the one-parameter correlations are poor. We find that many side chains experience a previously unrecognized type of transition state environment in which specific, native-like interactions are formed, but hydrophobic burial dominates. Implications of these results to the design and analysis of protein folding studies are discussed.     

Detection of a hidden folding intermediate in the focal adhesion target domain: Implications for its function and folding

The focal adhesion target (FAT) domain of focal adhesion kinase has a four-helix bundle structure. Based on a hydrogen exchange-constrained computer simulation study and some indirect experimental results, it has been suggested that a partially unfolded state of the FAT domain with the N-terminal helix unfolded plays an important role in its biological function. Here, using a native-state hydrogen exchange method, we directly detected an intermediate with the N-terminal helix unfolded in a mutant (Y925E) of the FAT domain. In addition, kinetic folding studies on the FAT domain suggest that this intermediate exists on the native side of the rate-limiting transition state for folding. These results provide more direct evidence of the existence of the proposed intermediate and help to understand the folding mechanism of small single domain proteins.     

Conformational heterogeneity of an equilibrium folding intermediate quantified and mapped by scanning mutagenesis

It is challenging to experimentally define an energy landscape for protein folding that comprises multiple partially unfolded states. Experimental results are often ambiguous as to whether a non-native state is conformationally homogeneous. Here, we tested an approach combining systematic mutagenesis and a Br?nsted-like analysis to reveal and quantify conformational heterogeneity of folding intermediate states. Using this method, we resolved an otherwise apparently homogeneous equilibrium folding intermediate of Borrelia burgdorferi OspA into two conformationally distinct species and determined their relative populations. Furthermore, we mapped the structural differences between these intermediate species, which are consistent with the non-native species that we previously proposed based on native-state hydrogen exchange studies. When treated as a single state, the intermediate ensemble exhibited fractional Phi-values for mutations and Hammond-type behaviors that are often observed for folding transition states. We found that a change in relative population of the two species within the intermediate ensemble explains these properties well, suggesting that fractional Phi-values and Hammond-type behaviors exhibited by folding intermediates and transition states may arise more often from conformational heterogeneity than from a single partial structure. Our results are consistent with the presence of multiple minima in a rugged energy landscape predicted from theoretical studies. The method described here provides a promising means to probe a complex folding energy landscape.     

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.     

Relationship between the native-state hydrogen exchange and the folding pathways of barnase

The previous native-state hydrogen exchange experiment with barnase failed to detect any partially unfolded intermediate state which was contrary to the experimental results from kinetic deuterium hydrogen exchange pulse labeling and protein engineering studies. This has been taken to suggest that the native-state hydrogen exchange method cannot be used alone as an analytical tool to study the folding pathways of proteins. Here, we revisited the pulse labeling experiment with barnase and detected no stable folding intermediate. This finding allows a reconciliation of the native-state HX data and the folding pathway of barnase. Along with alternative theoretical interpretations for a curved chevron plot of protein folding, these data suggest that further investigation of the nature of the intermediate of barnase is needed.     

Thermodynamic characterization of an equilibrium folding intermediate of staphylococcal nuclease.

High-sensitivity differential scanning calorimetry and CD spectroscopy have been used to probe the structural stability and measure the folding/unfolding thermodynamics of a Pro117-->Gly variant of staphylococcal nuclease. It is shown that at neutral pH the thermal denaturation of this protein is well accounted for by a 2-state mechanism and that the thermally denatured state is a fully hydrated unfolded polypeptide. At pH 3.5, thermal denaturation results in a compact denatured state in which most, if not all, of the helical structure is missing and the beta subdomain apparently remains largely intact. At pH 3.0, no thermal transition is observed and the molecule exists in the compact denatured state within the 0-100 degrees C temperature interval. At high salt concentration and pH 3.5, the thermal unfolding transition exhibits 2 cooperative peaks in the heat capacity function, the first one corresponding to the transition from the native to the intermediate state and the second one to the transition from the intermediate to the unfolded state. As is the case with other proteins, the enthalpy of the intermediate is higher than that of the unfolded state at low temperatures, indicating that, under those conditions, its stabilization must be of an entropic origin. The folding intermediate has been modeled by structural thermodynamic calculations. Structure-based thermodynamic calculations also predict that the most probable intermediate is one in which the beta subdomain is essentially intact and the rest of the molecule unfolded, in agreement with the experimental data. The structural features of the equilibrium intermediate are similar to those of a kinetic intermediate previously characterized by hydrogen exchange and NMR spectroscopy.     

Detection of a hidden folding intermediate of the third domain of PDZ

The folding pathway of the third domain of PDZ from the synaptic protein PSD-95 was characterized using kinetic and equilibrium methods by monitoring the fluorescence signal from a Trp residue that is incorporated at a near-surface position. Kinetic folding of this domain showed multiple exponential phases, whereas unfolding showed a single exponential phase. The slow kinetic phases were attributed to isomerization of proline residues, since there are five proline residues in this domain. We found that the logarithms of the rate constants for the fast phase of folding and unfolding are linearly dependent on the concentrations of denaturant. The unfolding free energy derived from these rate constants at zero denaturant was close to the value measured using the equilibrium method, suggesting the absence of detectable sub-millisecond folding intermediates. However, native-state hydrogen exchange experiments detected a partially unfolded intermediate under native conditions. It was further confirmed by a protein engineering study. These data suggest that a hidden intermediate exists after the rate-limiting step in the folding of the third domain of PDZ.     

The folding pathway of barnase: the rate-limiting transition state and a hidden intermediate under native conditions

The nature of the rate-limiting transition state at zero denaturant (TS(1)) and whether there are hidden intermediates are the two major unsolved problems in defining the folding pathway of barnase. In earlier studies, it was shown that TS(1) has small phi values throughout the structure of the protein, suggesting that the transition state has either a defined partially folded secondary structure with all side chains significantly exposed or numerous different partially unfolded structures with similar stability. To distinguish the two possibilities, we studied the effect of Gly mutations on the folding rate of barnase to investigate the secondary structure formation in the transition state. Two mutations in the same region of a beta-strand decreased the folding rate by 20- and 50-fold, respectively, suggesting that the secondary structures in this region are dominantly formed in the rate-limiting transition state. We also performed native-state hydrogen exchange experiments on barnase at pD 5.0 and 25 degrees C and identified a partially unfolded state. The structure of the intermediate was investigated using protein engineering and NMR. The results suggest that the intermediate has an omega loop unfolded. This intermediate is more folded than the rate-limiting transition state previously characterized at high denaturant concentrations (TS(2)). Therefore, it exists after TS(2) in folding. Consistent with this conclusion, the intermediate folds with the same rate and denaturant dependence as the wild-type protein, but unfolds faster with less dependence on the denaturant concentration. These and other results in the literature suggest that barnase folds through partially unfolded intermediates that exist after the rate-limiting step. Such folding behavior is similar to those of cytochrome c and Rd-apocyt b(562). Together, we suggest that other small apparently two-state proteins may also fold through hidden intermediates.     

Lack of definable nucleation sites in the rate-limiting transition state of barnase under native conditions.

It has been shown that the burst-phase (submillisecond) intermediate of barnase, if it exists, can be only marginally more stable than the fully unfolded state at pH 6.3 and 25 degrees C. In the study reported here, no stable burst-phase intermediate could be detected, even in the presence of stabilizing salt (0.4 M Na(2)SO(4)). These results suggest that a burst-phase intermediate with even marginal stability does not exist. The absence of such an intermediate in turn suggests the need for re-examination of the rate-limiting transition state (RLTS) under native conditions, which was previously characterized by using a three-state model with a stable intermediate and protein engineering. Surprisingly, mutations throughout the structure of barnase do not significantly affect the folding rate, suggesting a lack of specific favorable interactions among the side-chains in the RLTS. This RLTS is clearly different from that previously characterized under denaturing conditions, indicating that changes take place in the RLTS under native and denaturing conditions. The occurrence of such changes is further supported by the observation that the unfolding rate constants of barnase and its mutants were divergent or convergent as a function of denaturant concentrations. Consistent with changes in the RLTS, a re-analysis of data from native-state hydrogen exchange studies has shown that the logarithm of the unfolding rate constant inflects down under low concentrations of denaturant. Here, we discuss in detail the question of whether changes in the RLTS involve a kinetically silent intermediate that occurs after the initial RLTS.     

Folding of apocytochrome c induced by the interaction with negatively charged lipid micelles proceeds via a collapsed intermediate state

Unfolded apocytochrome c acquires an alpha-helical conformation upon interaction with lipid. Folding kinetic results below and above the lipid's CMC, together with energy transfer measurements of lipid bound states, and salt-induced compact states in solution, show that the folding transition of apocytochrome c from the unfolded state in solution to a lipid-inserted helical conformation proceeds via a collapsed intermediate state (I(C)). This initial compact state is driven by a hydrophobic collapse of the polypeptide chain in the absence of the heme group and may represent a heme-free analogue of an early compact intermediate detected on the folding pathway of cytochrome c in solution. Insertion into the lipid phase occurs via an unfolding step of I(C) through a more extended state associated with the membrane surface (I(S)). While I(C) appears to be as compact as salt-induced compact states in solution with substantial alpha-helix content, the final lipid-inserted state (Hmic) is as compact as the unfolded state in solution at pH 5 and has an alpha-helix content which resembles that of native cytochrome c.     

Kinetic folding and unfolding of staphylococcal nuclease and its six mutants studied by stopped-flow circular dichroism

Kinetics of refolding and unfolding of staphylococcal nuclease and its six mutants, each carrying single or double amino acid substitutions, are studied by stopped-flow circular dichroism measurements. A transient kinetic intermediate formed within 10 ms after refolding starts possesses a substantial part of the N-domain core β-structure, whereas helices are formed at the later stages. The structure of the kinetic intermediate is less organized than the structure that is known to be formed by a nuclease 1-136 fragment. Only the refolding kinetics are affected by the mutations in all the mutants except two in which the mutations have changed the native structure. From this result and also from the locations of the mutation sites, the major N-terminal domain of the nuclease in the transition state of folding has a structure nearly identical to the native one. On the other hand, the minor C-terminal domain has previously been shown to be still disorganized in the transition state. The effects of the amino acid substitutions on the stability of the native and the transition states are in good agreement with the changes in the hydration free energy, expected for the corresponding amino acid replacements in the unfolded polypeptide. Since side chains of all the mutated residues are not accessible to solvent in the native structure, the result suggests that it is the unfolded state that is mainly affected by the mutations. © 1995 Wiley-Liss, Inc.     

Test for cooperativity in the early kinetic intermediate in lysozyme folding

During the folding of many proteins, collapsed globular states are formed prior to the native structure. The role of these states for the folding process has been widely discussed. Comparison with properties of synthetic homo and heteropolymers had suggested that the initial collapse represented a shift of the ensemble of unfolded conformations to more compact states without major energy barriers. We investigated the folding/unfolding transition of a collapsed state, which transiently populates early in lysozyme folding. This state forms within the dead-time of stopped-flow mixing and it has been shown to be significantly more compact and globular than the denaturant-induced unfolded state. We used the GdmCl-dependence of the dead-time signal change to characterize the unfolding transition of the burst phase intermediate. Fluorescence and far-UV CD give identical unfolding curves, arguing for a cooperative two-state folding/unfolding transition between unfolded and collapsed lysozyme. These results show that collapse leads to a distinct state in the folding process, which is separated from the ensemble of unfolded molecules by a significant energy barrier. NMR, fluorescence and small angle X-ray scattering data further show that some local interactions in unfolded lysozyme exist at denaturant concentrations above the coil-collapse transition. These interactions might play a crucial role in the kinetic partitioning between fast and slow folding pathways.     

The Folding Trajectory of RNase H Is Dominated by Its Topology and Not Local Stability: A Protein Engineering Study of Variants that Fold via Two-State and Three-State Mechanisms

Proteins can sample a variety of partially folded conformations during the transition between the unfolded and native states. Some proteins never significantly populate these high-energy states and fold by an apparently two-state process. However, many proteins populate detectable, partially folded forms during the folding process. The role of such intermediates is a matter of considerable debate. A single amino acid change can convert Escherichia coli ribonuclease H from a three-state folder that populates a kinetic intermediate to one that folds in an apparent two-state fashion. We have compared the folding trajectories of the three-state RNase H and the two-state RNase H, proteins with the same native-state topology but altered regional stability, using a protein engineering approach. Our data suggest that both versions of RNase H fold through a similar trajectory with similar high-energy conformations. Mutations in the core and the periphery of the protein affect similar aspects of folding for both variants, suggesting a common trajectory with folding of the core region followed by the folding of the periphery. Our results suggest that formation of specific partially folded conformations may be a general feature of protein folding that can promote, rather than hinder, efficient folding.     

Monomer topology defines folding speed of heptamer

Small monomeric proteins often fold in apparent two-state processes with folding speeds dictated by their native-state topology. Here we test, for the first time, the influence of monomer topology on the folding speed of an oligomeric protein: the heptameric cochaperonin protein 10 (cpn10), which in the native state has seven beta-barrel subunits noncovalently assembled through beta-strand pairing. Cpn10 is a particularly useful model because equilibrium-unfolding experiments have revealed that the denatured state in urea is that of a nonnative heptamer. Surprisingly, refolding of the nonnative cpn10 heptamer is a simple two-state kinetic process with a folding-rate constant in water (2.1 sec(-1); pH 7.0, 20 degrees C) that is in excellent agreement with the prediction based on the native-state topology of the cpn10 monomer. Thus, the monomers appear to fold as independent units, with a speed that correlates with topology, although the C and N termini are trapped in beta-strand pairing with neighboring subunits. In contrast, refolding of unfolded cpn10 monomers is dominated by a slow association step.     

Mutational analysis of the BPTI folding pathway: II. Effects of aromatic-->leucine substitutions on folding kinetics and thermodynamics.

The rates of the individual steps in the disulfide-coupled folding and unfolding of eight BPTI variants, each containing a single aromatic to leucine amino acid replacement, were measured. From this analysis, the contributions of the four phenylalanine and four tyrosine residues to the stabilities of the native protein and the disulfide-bonded folding intermediates were determined. While the substitutions were found to destabilize the native protein by 2 to 7 kcal/mol, they had significantly smaller effects on the intermediates that represent the earlier stages of folding, even when the site of the substitution was located within the ordered regions of the intermediates. These results suggest that stabilizing interactions contribute less to conformational stability in the context of a partially folded intermediate than in a fully folded native protein, perhaps because of decreased cooperativity among the individual interactions. The kinetic analysis also provides new information about the transition states associated with the slowest steps in folding and unfolding, supporting previous suggestions that these transition states are extensively unfolded. Although the substitutions caused large changes in the distribution of folding intermediates and in the rates of some steps in the folding pathway, the kinetically-preferred pathway for all of the variants involved intramolecular disulfide rearrangements, as observed previously for the wild-type protein. These results suggest that the predominance of the rearrangement mechanism reflects conformational constraints present relatively early in the folding pathway.     

Perturbations of the denatured state ensemble: modeling their effects on protein stability and folding kinetics.

By considering the denatured state of a protein as an ensemble of conformations with varying numbers of sequence-specific interactions, the effects on stability, folding kinetics, and aggregation of perturbing these interactions can be predicted from changes in the molecular partition function. From general considerations, the following conclusions are drawn: (1) A perturbation that enhances a native interaction in denatured state conformations always increases the stability of the native state. (2) A perturbation that promotes a non-native interaction in the denatured state always decreases the stability of the native state. (3) A change in the denatured state ensemble can alter the kinetics of aggregation and folding. (4) The loss (or increase) in stability accompanying two mutations, each of which lowers (or raises) the free energy of the denatured state, will be less than the sum of the effects of the single mutations, except in cases where both mutations affect the same set of partially folded conformations. By modeling the denatured state as the ensemble of all non-native conformations of hydrophobic-polar (HP) chains configured on a square lattice, it can be shown that the stabilization obtained from enhancement of native interactions derives in large measure from the avoidance of non-native interactions in the D state. In addition, the kinetic effects of fixing single native contacts in the denatured state or imposing linear gradients in the HH contact probabilities are found, for some sequences, to significantly enhance the efficiency of folding by a simple hydrophobic zippering algorithm. Again, the dominant mechanism appears to be avoidance of non-native interactions. These results suggest stabilization of native interactions and imposition of gradients in the stability of local structure are two plausible mechanisms involving the denatured state that could play a role in the evolution of protein folding and stability.     

Local interactions and the optimization of protein folding

The role of local interactions in protein folding has recently been the subject of some controversy. Here we investigate an extension of Zwanzig's simple and general model of folding in which local and nonlocal interactions are represented by functions of single and multiple conformational degrees of freedom, respectively. The kinetics and thermodynamics of folding are studied for a series of energy functions in which the energy of the native structure is fixed, but the relative contributions of local and nonlocal interactions to this energy are varied over a broad range. For funnel shaped energy landscapes, we find that 1) the rate of folding increases, but the stability of the folded state decreases, as the contribution of local interactions to the energy of the native structure increases, and 2) the amount of native structure in the unfolded state and the transition state vary considerably with the local interaction strength. Simple exponential kinetics and a well-defined free energy barrier separating folded and unfolded states are observed when nonlocal interactions make an appreciable contribution to the energy of the native structure; in such cases a transition state theory type approximation yields reasonably accurate estimates of the folding rate. Bumps in the folding funnel near the native state, which could result from desolvation effects, side chain freezing, or the breaking of nonnative contacts, significantly alter the dependence of the folding rate on the local interaction strength: the rate of folding decreases when the local interaction strength is increased beyond a certain point. A survey of the distribution of strong contacts in the protein structure database suggests that evolutionary optimization has involved both kinetics and thermodynamics: strong contacts are enriched at both very short and very long sequence separations. Proteins 29:282–291, 1997. © 1997 Wiley-Liss, Inc.