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.     

Explicit-chain model of native-state hydrogen exchange: implications for event ordering and cooperativity in protein folding

Native-state hydrogen exchange experiments on several proteins have revealed partially unfolded conformations with diverse stabilities. These equilibrium observations have been used to support kinetic arguments that folding proceeds via a sequential "pathway." This interpretative logic is evaluated here by analyzing the relationship between thermodynamic behavior and folding kinetics in a class of simplified lattice protein models. The chain models studied have varying degrees of cooperative interplay (coupling) between local helical conformational preference and favorable nonlocal interactions. When model cooperativity is high, as native conditions are weakened, "isotherms" of free energy of exchange for residues belonging to the same helix merge together before global unfolding. The point of merger depends on the model energetic favorability of the helix. This trend is similar to the corresponding experimental observations. Kinetically, we find that the ordering of helix formation in the very last stage of native core assembly tends to follow the stabilities of their converged isotherms. In a majority (but not all) of folding trajectories, the final assembly of helices that are thermodynamically more stable against exchange precedes that of helices that are less stable against exchange. These model features are in partial agreement with common experimental interpretations. However, the model results also underscore the ensemble nature of the folding process: the kinetics of helix formation is not a discrete, strictly "all-or-none" process as that envisioned by certain non-explicit-chain models. Helices generally undergo many cycles of partial formation and dissolution before their conformations are fixed in the final assembly stage of folding, a kinetic stage that takes up only approximately 2% of the average folding time in the present model; and the ordering of the helices' final assembly in some trajectories can be different from the dominant ordering stipulated by the exchange isotherms.     

Role of native-state structure in rubredoxin native-state hydrogen exchange

Amide exchange rates were measured for Pyrococcus furiosus (Pf) rubredoxin substituted with either Zn(II), Ga(III), or Ge(IV). Base-catalyzed exchange rate constants increase up to 3000-fold per unit charge for the highly protected amides surrounding the active site metal, yielding apparent residue-specific conformational energy decreases of more than 8 kcal/mol in a comparison of the Zn(II)- and Ge(IV)-substituted proteins. However, the exchange kinetics for many of the other amides of the protein are insensitive to these metal substitutions. These differential rates are inversely correlated with the distance between the amide nitrogen and the metal in the X-ray structure, out to a distance of at least 12 A, consistent with an electrostatic potential-dependent shifting of the amide nitrogen pK. This strongly correlated distance dependence is consistent with a nativelike structure for the exchange-competent conformations. The electric field potential within the interior of the rubredoxin structure gives rise to a change of as much as a million-fold in the rate for the exchange-competent state of the individual amide hydrogens. Nevertheless, the strength of these electrostatic interactions in Pf rubredoxin appears to be comparable to those previously reported within other proteins. As a result, contrary to the conventional analysis of hydrogen exchange data, for exchange processes that occur via nonglobal transitions, the residual conformational structure will often modulate the observed rates. Although this necessarily complicates the estimation of the conformational equilibria of these exchange-competent states, this dependence on residual structure can provide insight into the conformation of these transient states.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

EX1 hydrogen exchange and protein folding

Slow amide hydrogen exchange is an increasingly popular tool for investigating structure and function in proteins. The kinetic model for slow hydrogen exchange has two limits, called EX2 and EX1, wherein the thermodynamics and kinetics of protein motions, respectively, are reported by the exchange data. While many laboratories have demonstrated that EX2 exchange can indeed provide accurate results regarding the thermodynamics of protein stability, the potential of EX1 exchange to follow the kinetics of protein unfolding and folding is only beginning to be realized. EX1 hydrogen exchange has advantages over more traditional folding experiments: it provides single-residue resolution, as well as whole-molecule information, the latter of which can be interpreted in terms of the cooperativity of unfolding. However, key questions remain regarding the interpretation of EX1 hydrogen exchange.     

The hydrogen exchange core and protein folding.

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

Equilibrium amide hydrogen exchange and protein folding kinetics

The classical Linderstrøm-Lang hydrogen exchange (HX) model is extended to describe the relationship between the HX behaviors (EX1 and EX2) and protein folding kinetics for the amide protons that can only exchange by global unfolding in a three-state system including native (N), intermediate (I), and unfolded (U) states. For these slowly exchanging amide protons, it is shown that the existence of an intermediate (I) has no effect on the HX behavior in an off-pathway three-state system (IUN). On the other hand, in an on-pathway three-state system (UIN), the existence of a stable folding intermediate has profound effect on the HX behavior. It is shown that fast refolding from the unfolded state to the stable intermediate state alone does not guarantee EX2 behavior. The rate of refolding from the intermediate state to the native state also plays a crucial role in determining whether EX1 or EX2 behavior should occur. This is mainly due to the fact that only amide protons in the native state are observed in the hydrogen exchange experiment. These new concepts suggest that caution needs to be taken if one tries to derive the kinetic events of protein folding from equilibrium hydrogen exchange experiments.     

Folding of a de novo designed native-like four-helix bundle protein

The folding of a model native-like dimeric four-helix bundle protein, (alpha(2))(2), was investigated using guanidine hydrochloride, hydrostatic pressure, and low temperature. Unfolding by guanidine hydrochloride followed by circular dichroism and intrinsic fluorescence spectroscopy revealed a highly cooperative transition between the native-like and unfolded states, with free energy of unfolding determined from CD data, DeltaG(unf) = 14.3 +/- 0.8 kcal/mol. However, CD and intrinsic fluorescence data were not superimposable, indicating the presence of an intermediate state during the folding transition. To stabilize the folding intermediate, we used hydrostatic pressure and low temperature. In both cases, dissociation of the dimeric native-like (alpha(2))(2) into folded monomers (alpha(2)) was observed. van't Hoff analysis of the low temperature experiments, assuming a two-state dimer 171-monomer transition, yielded a free energy of dissociation of (alpha(2))(2) of DeltaG(diss) = 11.4 +/- 0.4 kcal/mol, in good agreement with the free energy determined from pressure dissociation experiments (DeltaG(diss) = 10.5 +/- 0.1 kcal/mol). Binding of the hydrophobic fluorescent probe 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid (bis-ANS) to the pressure- and cold-dissociated states of (alpha(2))(2) indicated the existence of molten-globule monomers. In conclusion, we demonstrate that the folding pathway of (alpha(2))(2) can be described by a three-state transition including a monomeric molten globule-like state.     

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.     

Diffusion-collision model study of misfolding in a four-helix bundle protein.

Proteins with complex folding kinetics will be susceptible to misfolding at some stage in the folding process. We simulate this problem by using the diffusion-collision model to study non-native kinetic intermediate misfolding in a four-helix bundle protein. We find a limit on the size of the pairwise hydrophobic area loss in non-native intermediates, such that burying above this limit creates long-lasting non-native kinetic intermediates that would disrupt folding and prevent formation of the native state. Our study of misfolding suggests a method for limiting the production of misfolded kinetic intermediates for helical proteins and could, perhaps, lead to more efficient production of proteins in bulk.     

Residue-wise conformational stability of DLC8 dimer from native-state hydrogen exchange

Dynein light chain (DLC8) is the smallest subunit of the dynein motor complex, which is known to act as a cargo adaptor in intracellular trafficking. The protein exists as a pure dimer at physiological pH and a completely folded monomer below pH 4. Here, we have determined the energy landscape of the dimeric protein using a combination of optical techniques and native-state hydrogen exchange of amide groups, the former giving the global features and the latter yielding the residue level details. The data indicated the presence of intermediates along the equilibrium unfolding transition. The hydrogen exchange data suggested that the molecule has differential stability in its various segments. We deduce from the free energy data that the antiparallel beta-sheets (beta4 and beta5) that form the hydrophobic core of the protein and the alpha2 helix, all of which are highly protected with regard to hydrogen exchange, contribute significantly to the initial step of the protein folding mechanism. Denaturant-dependent hydrogen exchange indicated further that some amides exchange via local fluctuations, whereas there are others which exchange via global unfolding events. Implications of these to cargo adaptability of the dimer are discussed.     

Disulfide crosslinks to probe the structure and flexibility of a designed four-helix bundle protein.

The introduction of disulfide crosslinks is a generally useful method by which to identify regions of a protein that are close together in space. Here we describe the use of disulfide crosslinks to investigate the structure and flexibility of a family of designed 4-helix bundle proteins. The results of these analyses lend support to our working model of the proteins' structure and suggest that the proteins have limited main-chain flexibility.     

Human interleukin 4. The solution structure of a four-helix bundle protein.

Heteronuclear 13C and 15N three-dimensional nuclear magnetic resonance (n.m.r.) techniques have been used to determine the solution structure of human interleukin 4, a four-helix bundle protein. A dynamical simulated annealing protocol was used to calculate an ensemble of structures from an n.m.r. data set of 1735 distance restraints, 101 phi angle restraints and 27 pairs of hydrogen bond restraints. The protein structure has a left-handed up-up-down-down topology for the four helices with the two long overhand loops in the structure being connected by a short section of irregular antiparallel beta-sheet. Analysis of the side-chains in the protein shows a clustering of hydrophobic residues, particularly leucines, in the core of the bundle with the side-chains of charged residues being located on the protein surface. The solution structure has been compared with a recent structure prediction for human interleukin 4 and with crystal structures of other helix bundle proteins.     

Hydrogel biopolymer created from the self-assembly of a designed protein containing a four-helix bundle forming motif

A protein hydrogel system based on the assembly of a four-helix bundle motif was proposed and synthesized by genetic engineering methods. This new polypeptide, named GBH1, consists of identical amphipathic helices of 22 residues in length oriented in opposite fashion to one another at each end of a polypeptide with a total length of 227 amino acids. The middle portion of the polypeptide (residues 79-147) is an unstructured random coil. The region between the amphipathic and unstructured segment is an α-helical stretch (23-78 and 148-204) not possessing a sequence compatible with a coiled-coil conformation, but rather possesses regions that have overwinding of the helix. The thermal unfolding of GBH1 shows more than one inflection point (T(m1) = 30.5 °C, T(m2) = 64.6 °C), indicative of a partially unfolded intermediate and, thus, multiple interactions in the folded state. A qualitative assessment of hydrogel formation with varying pH showed that acidic conditions did not support the gel state, indirectly indicating that the proposed four-helix bundle is the major cross-linking structure and not a leucine zipper motif. Scanning electron microscopy reveals a network of interacting protein molecules forming a spongelike matrix with numerous pores that would be occupied with water molecules.     

Computer modeling and folding of four-helix bundles

In the context of simplified models of globular proteins, the requirements for the unique folding to a four-helix bundle have been addressed through a new Monte Carlo procedure. In particular, the relative importance of secondary versus tertiary interactions in determining the nature of the folded structure is examined. Various cases spanning the extremes where tertiary interactions completely dominate to that where tertiary interactions are negligible have been explored. Not surprisingly, the folding to unique four-helix bundles is found to depend on an adequate balance of the secondary and tertiary interactions. Moreover, because the simplified model is composed of spheres representing α-carbons and side chains, the geometry of the latter being based on small real amino acids, the role played by the side chains, and the problems associated with packing and hard-core repulsions, are considered. Also, possible folding intermediates and their relationship with the experimentally observed molten globule state are explored. From these studies, a general set of rules is extracted which should aid in the further design of more detailed protein models adequate to more fully investigate the protein folding problem. Finally, the relationship between our conclusions and experimental work with specifically designed sequences is briefly discussed. © 1993 Wiley-Liss, Inc.     

Relationship between folding and function in a sequence-specific miniature DNA-binding protein

Previously, we have described a miniature protein-based approach to the design of molecules that bind DNA or protein surfaces with high affinity and specificity. In this approach, the small, well-folded protein avian pancreatic polypeptide acts as a scaffold to present and stabilize an alpha-helical or PPII-helical recognition epitope. The first miniature protein designed in this way, a molecule called p007, presents the alpha-helical recognition epitope found on the bZIP protein GCN4 and binds DNA with nanomolar affinity and exceptional specificity. In this work we use alanine-scanning mutagenesis to explore the contributions of 29 p007 residues to DNA affinity, specificity, and secondary structure. Virtually every residue within the p007 alpha-helix, and most residues within the p007 PPII helix, contribute to both DNA affinity and specificity. These residues include those introduced to make specific and nonspecific DNA contacts, as well as those that complete the miniature protein core. Moreover, there exists a direct correlation between the affinity of a p007 variant for specific DNA and the ability of that variant to select for specific DNA over nonspecific DNA. Although we observe no correlation between alpha-helicity and affinity, we observe a limited correlation between alpha-helicity and sequence specificity that emphasizes the role of coupled binding/folding in the function of p007. Our results imply that formation of a highly evolved set of protein.DNA contacts in the context of a well-packed hydrophobic core, and not the extent of intrinsic alpha-helical structure, is the primary determinant of p007 function.     

Investigation of the structural stability of hUBF HMG box 5 by native-state hydrogen exchange

HMG box 5 of human upstream binding factor (hUBF) consists of three alpha-helices arranged in an L-shape with a hydrophobic core embraced by these helices and stabilized by extensive hydrophobic interactions between nonpolar residues around the core. The GdmCl-induced equilibrium unfolding transition of HMG box 5 of hUBF was monitored by both circular dichroism (CD) and fluorescence spectra. A cooperative two-state unfolding process was observed. The unfolding free energy, DeltaGU(D2O), and the cooperativity of the unfolding reaction, m, are 4.6 +/- 0.16 kcal x mol-1 and 1.62 +/- 0.06 kcal x mol-1 x M-1, respectively. Native-state hydrogen exchange (NHX) experiments under EX2 conditions were performed. NHX results clearly show that the hydrophobic core among the three helices is a slow-exchange core. The three helices would not contribute equally to the stability of the native protein. Helix 3 appears to contribute the least to the stability. The NHX data have also allowed the local, subglobal, and global unfolding structures of hUBF HMG box 5 to be dissected, and common global and subglobal unfolding units were successfully detected.     

Relationship between protein folding kinetics and amino acid properties

The successful prediction of protein-folding rates based on the sequence-predicted secondary structure suggests that the folding rates might be predicted from sequence alone. To pursue this question, we directly predict the folding rates from amino acid sequences, which do not require any information on secondary or tertiary structure. Our work achieves 88% correlation with folding rates determined experimentally for proteins of all folding types and peptide, suggesting that almost all of the information needed to specify a protein's folding kinetics and mechanism is comprised within its amino acid sequence. The influence of residue on folding rate is related to amino acid properties. Hydrophobic character of amino acids may be an important determinant of folding kinetics, whereas other properties, size, flexibility, polarity and isoelectric point, of amino acids have contributed little to the folding rate constant.     

Desolvation shell of hydrogen bonds in folded proteins,protein complexes and folding pathways

A few backbone hydrogen bonds (HBS) in native protein folds are poorly protected from water attack: their desolvation shell contains an inordinately low number of hydrophobic residues. Thus, an approach by solvent-structuring moieties of a binding partner should contribute significantly to enhance their stability. This effect represents an important factor in the site specificity inherent to protein binding, as inferred from a strong correlation between poorly desolvated HBs and binding sites. The desolvation shells were also examined in a dynamic context: except for a few singular under-protected bonds, the size of desolvation shells is preserved along the folding trajectory.