Relationship between the native-state hydrogen exchange and folding pathways of a four-helix bundle protein

The hydrogen exchange behavior of a four-helix bundle protein in low concentrations of denaturant reveals some partially unfolded forms that are significantly more stable than the fully unfolded state. Kinetic folding of the protein, however, is apparently two-state in the absence of the accumulation of early folding intermediates. The partially unfolded forms are either as folded as or more folded than the rate-limiting transition state and appear to represent the major intermediates in a folding and unfolding reaction. These results are consistent with the suggestion that partially unfolded intermediates may form after the rate-limiting step for small proteins with apparent two-state folding kinetics.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Thermodynamic and kinetic exploration of the energy landscape of Borrelia burgdorferi OspA by native-state hydrogen exchange

We report a native-state hydrogen-exchange (HX) method to simultaneously obtain both thermodynamic and kinetic information on the formation of multiple excited states in a folding energy landscape. Our method exploits the inherent dispersion and pH dependence of the intrinsic HX rates to cover both the EX2 (thermodynamic) and EX1 (kinetic) regimes. At each concentration of denaturant, HX measurements are performed over a range of pH values. Using this strategy, we dissected Borrelia burgdorferi OspA, a predominantly beta-sheet protein containing a unique single-layer beta-sheet, into five cooperative units and postulated excited states predominantly responsible for HX. More importantly, we determined the interconversion rates between these excited states and the native state. The use of both thermodynamic and kinetic information from native-state HX enabled us to construct a folding landscape of this 28kDa protein, including local minima and maxima, and to discriminate on-pathway and off-pathway intermediates. This method, which we term EX2/EX1 HX, should be a powerful tool for characterizing the complex folding mechanisms exhibited by the majority of proteins.     

The kinetic pathway of folding of barnase

To search for folding intermediates, we have examined the folding and unfolding kinetics of wild-type barnase and four representative mutants under a wide range of conditions that span two-state and multi-state kinetics. The choice of mutants and conditions provided in-built controls for artifacts that might distort the interpretation of kinetics, such as the non-linearity of kinetic and equilibrium data with concentration of denaturant. We measured unfolding rate constants over a complete range of denaturant concentration by using by 1H/2H-exchange kinetics under conditions that favour folding, conventional stopped-flow methods at higher denaturant concentrations and continuous flow. Under conditions that favour multi-state kinetics, plots of the rate constants for unfolding against denaturant concentration fitted quantitatively to the equation for three-state kinetics, with a sigmoid component for a change of rate determining step, as did the refolding kinetics. The position of the transition state on the reaction pathway, as measured by solvent exposure (the Tanford beta value) also moved with denaturant concentration, fitting quantitatively to the same equations with a change of rate determining step. The sigmoid behaviour disappeared under conditions that favoured two-state kinetics. Those data combined with direct structural observations and simulation support a minimal reaction pathway for the folding of barnase that involves two detectable folding intermediates. The first intermediate, I(1), is the denatured state under physiological conditions, D(Phys), which has native-like topology, is lower in energy than the random-flight denatured state U and is suggested by molecular dynamics simulation of unfolding to be on-pathway. The second intermediate, I(2), is high energy, and is proven by the change in rate determining step in the unfolding kinetics to be on-pathway. The change in rate determining step in unfolding with structure or environment reflects the change in partitioning of this intermediate to products or starting materials.     

Stability and folding of the protein complexes of barnase.

Native-like complexes of proteins, formed by the association of two complementary fragments comprising the entire sequence of the protein, can be used to gain insight into the stability and folding of the intact protein. We have studied the structural, thermodynamic and kinetic properties of four barnase complexes, with the cleavage site at different positions of the amino-acid chain (CB36, at position 36; CB56, at position 56; CB68, at position 68; and CB79, at position 79). The four barnase complexes have native-like structure as shown by fluorescence, far-and near-UV CD, size-exclusion chromatography and NMR. The NMR characterization indicated that the structural changes were mainly located in regions close to the cleavage site. The main core of the protein was fully formed and the overall structure was similar to that of intact barnase. The thermal and chemical denaturation showed that all complexes were substantially destabilized. CB56 displayed two denaturation transitions, probably because of the presence of partially folded conformations around the cleavage site. The rate constant for the association/folding of fragments decreased with the decreasing length of the C-terminal fragment. Thus, the larger the fragment (and, consequently, the larger the amount of residual native-like structure), the faster the association. These findings are consistent with the proposed model of barnase folding.     

Prediction of folding rates and transition-state placement from native-state geometry

A variety of experimental and theoretical studies have established that the folding process of monomeric proteins is strongly influenced by the topology of the native state. In particular, folding times have been shown to correlate well with the contact order, a measure of contact locality. Our investigation focuses on identifying additional topologic properties that correlate with experimentally measurable quantities, such as folding rates and transition-state placement, for both two- and three-state folders. The validation against data from 40 experiments shows that a particular topological property that measures the interdependence of contacts, termed cliquishness or clustering coefficient, can account with statistically significant accuracy both for the transition state placement and especially for folding rates. The observed correlations can be further improved by optimally combining the distinct topological information captured by cliquishness and contact order.     

COSMIC analysis of the major alpha-helix of barnase during folding

The structures of transition states and intermediates in protein folding may be analysed by protein engineering methods that remove simple interactions that stabilize the folded state. We have now extended the range and reliability of the procedure by using the COSMIC (Combination of Sequential Mutant Interaction Cycles) technique, in which a series of double-mutant cycles is constructed. In each cycle, the side-chains of two amino acid residues that interact in the folded state are mutated separately and together. Kinetic and equilibrium measurements on folding for each cycle show unambiguously whether or not two residues interact during protein folding. A series of such cycles has been constructed to leapfrog along the major alpha-helix of barnase, comprising residues 6 to 18. The helix is found to be intact from its C terminus to residue 12 but begins to unwind towards the N terminus in both the transition state for unfolding and in a folding intermediate.     

Computing folding pathways between RNA secondary structures

Given an RNA sequence and two designated secondary structures A, B, we describe a new algorithm that computes a nearly optimal folding pathway from A to B. The algorithm, RNAtabupath, employs a tabu semi-greedy heuristic, known to be an effective search strategy in combinatorial optimization. Folding pathways, sometimes called routes or trajectories, are computed by RNAtabupath in a fraction of the time required by the barriers program of Vienna RNA Package. We benchmark RNAtabupath with other algorithms to compute low energy folding pathways between experimentally known structures of several conformational switches. The RNApathfinder web server, source code for algorithms to compute and analyze pathways and supplementary data are available at http://bioinformatics.bc.edu/clotelab/RNApathfinder.     

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.     

A kinetic folding intermediate probed by native state hydrogen exchange

Stopped-flow fluorescence studies on the N-terminal domain of rat CD2 (CD2.d1) have demonstrated that folding from the fully denatured state (U) proceeds via the transient accumulation of an apparent intermediate (I) in a so-called burst phase that precedes the rate-limiting transition leading to the native state (N). A previous pH-dependent equilibrium hydrogen exchange (HX) study identified a subset of amides in CD2.d1 which, under EX2 conditions, exchange from N with free energies greater than or equal to the free energy difference between the N and I states calculated from the stopped-flow data. Under EX1 conditions the rates of HX for these amides tend towards an asymptote that matches the global unfolding rate calculated from the stopped-flow data, suggesting that exchange for these amides requires traversing the N-to-I transition state barrier. Exchange for these amides presumably occurs from exchange-competent forms comprising the kinetic burst phase therefore. To explore this idea further, native state HX (NHX) data have been collected for CD2.d1 under EX2 conditions using denaturant concentrations which span either side of the denaturant concentration where, according to the stopped-flow data, the apparent U and I states are iso-energetic. The data fit to a two-component, sub-global (sg)/global (g) NHX mechanism, yielding Delta G and m value parameters (where the m value is a measure of hydrocarbon solvation). Regression analysis demonstrates that the (m(sg), Delta G(sg)) and (m(g), Delta G(g)) values calculated for this subset of amides correspond with those describing the kinetic burst phase transition. This result confirms the ability of the NHX technique to explore the structural and energetic properties of kinetic folding intermediates.     

A Gaussian statistical mechanical model for the equilibrium thermodynamics of barnase folding

Given an all non-hydrogen-atom potential function that implicitly includes solvation effects, it is possible to adjust its parameters to favor the correct native structure for several proteins over decoys produced by ungapped threading. It is also possible to further train it to reproduce the experimental free energy of unfolding in aqueous solution at 298 K for wild-type barnase and 66 mutants. For this, the native state is represented by the crystal structure at a single energy level with a calculated low degeneracy; the denatured state is represented by the extended conformation and a high calculated degeneracy. The same two-state model can be extended to account for the stability of all 67 sequences toward urea denaturation at 298 K by building in a solvation term that depends on urea concentration. With the addition of one more parameter set to give the correct heat capacity of unfolded barnase in solution, it is possible to approximate the experimental thermodynamics of barnase thermal denaturation: melting temperature, width of thermal transition, deltaG, deltaH, deltaS, and deltaCp. This requires a novel sort of statistical mechanical model where the two states each have a Gaussian density of microscopic state distribution as a function of energy.     

Kinetics of unfolding and folding from amide hydrogen exchange in native ubiquitin

Amide hydrogen (NH) exchange is one of the few experimental techniques with the potential for determining the thermodynamics and kinetics of conformational motions at nearly every residue in native proteins. Quantitative interpretation of NH exchange in terms of molecular motions relies on a simple two-state kinetic model: at any given slowly exchanging NH, a closed or exchange-incompetent conformation is in equilibrium with an open or exchange-competent conformation. Previous studies have demonstrated the accuracy of this model in measuring conformational equilibria by comparing exchange data with the thermodynamics of protein unfolding. We report here a test of the accuracy of the model in determining the kinetics of conformational changes in native proteins. The kinetics of folding and unfolding for ubiquitin have been measured by conventional methods and compared with those derived from a comprehensive analysis of the pH dependence of exchange in native ubiquitin. Rate constants for folding and unfolding from these two very different types of experiments show good agreement. The simple model for NH exchange thus appears to be a robust framework for obtaining quantitative information about molecular motions in native proteins.