Protein hydrogen exchange mechanism: local fluctuations

Experiments were done to study the dynamic structural motions that determine protein hydrogen exchange (HX) behavior. The replacement of a solvent-exposed lysine residue with glycine (Lys8Gly) in a helix of recombinant cytochrome c does not perturb the native structure, but it entropically potentiates main-chain flexibility and thus can promote local distortional motions and large-scale unfolding. The mutation accelerates amide hydrogen exchange of the mutated residue by about 50-fold, neighboring residues in the same helix by less, and residues elsewhere in the protein not at all, except for Leu98, which registers the change in global stability. The pattern of HX changes shows that the coupled structural distortions that dominate exchange can be several residues in extent, but they expose to exchange only one amide NH at a time. This "local fluctuation" mode of hydrogen exchange may be generally recognized by disparate near-neighbor rates and a low dependence on destabilants (denaturant, temperature, pressure). In contrast, concerted unfolding reactions expose multiple neighboring amide NHs with very similar computed protection factors, and they show marked destabilant sensitivity. In both modes, ionic hydrogen exchange catalysts attack from the bulk solvent without diffusing through the protein matrix.     

Predicting the equilibrium protein folding pathway: Structure-based analysis of staphylococcal nuclease

The equilibrium folding pathway of staphylococcal nucleas (SNase) has been approximated using a statistical thermodynamic formalism that utilizes the high-resolution structure of the native state as a template to generate a large ensemble of partially folded states. Close to 400,000 different states ranging from the native to the completely unfolded states were included in the analysis. The probability of each state was estimated using an empirical structural parametrization of the folding energetics. It is shown that this formalism predicts accurately the stability of the protein, the cooperativity of the folding/unfolding transition observed by differential scanning calorimetry (DSC) or urea denaturation and the thermodynamic parameters for unfolding. More importantly, this formalism provides a quantitative account of the experimental hydrogen exchange protection factors measured under native conditions for SNase. These results suggest that the computer-generated distribution of states approximates well the ensemble of conformations existing in solution. Furthermore, this formalism represents the first model capable of quantitatively predicting within a unified framework the probability distribution of states seen under native conditions and its change upon unfolding. © 1997 Wiley-Liss, Inc.     

Hydrogen exchange study of canine milk lysozyme: stabilization mechanism of the molten globule

The native state (1)H, (15)N resonance assignment of 123 of the 128 nonproline residues of canine milk lysozyme has enabled measurements of the amide hydrogen exchange of over 70 amide hydrogens in the molten globule state. To elucidate the mechanism of protein folding, the molten globule state has been studied as a model of the folding intermediate state. Lysozyme and alpha-lactalbumin are homologous to each other, but their equilibrium unfolding mechanisms differ. Generally, the folding mechanism of lysozyme obeys a two-state model, whereas that of alpha-lactalbumin follows a three-state model. Exceptions to this rule are equine and canine milk lysozymes, which exhibit a partially unfolded state during the equilibrium unfolding; this state resembles the molten globule state of alpha-lactalbumin but with extreme stability. Study of the molten globules of alpha-lactalbumin and equine milk lysozyme showed that the stabilities of their alpha-helices are similar, despite the differences in the thermodynamic stability of their molten globule states. On the other hand, our hydrogen exchange study of the molten globule of canine milk lysozyme showed that the alpha-helices are more stabilized than in alpha-lactalbumin or equine milk lysozyme and that this enhanced stability is caused by the strengthened cooperative interaction between secondary structure elements. Thus, our results underscore the importance of the cooperative interaction in the stability of the molten globule state.     

Equilibrium and kinetic studies on folding of the authentic and recombinant forms of human alpha-lactalbumin by circular dichroism spectroscopy

The equilibrium and kinetics of the unfolding and refolding of authentic and recombinant human alpha-lactalbumin, the latter of which had an extra methionine residue at the N-terminus, were studied by circular dichroism spectroscopy, and the results were compared with the results for bovine and goat alpha-lactalbumins obtained in our previous studies. As observed in the bovine and goat proteins, the presence of the extra methionine residue in the recombinant protein remarkably destabilized the native state, and the destabilization was entirely ascribed to an increase in the rate of unfolding. The thermodynamic stability of the native state against the unfolded state was lower, and the thermodynamic stability of the molten globule state against the unfolded state was higher for the human protein than for the other alpha-lactalbumins previously studied. Thus, the population of the molten globule intermediate was higher during the equilibrium unfolding of human alpha-lactalbumin by guanidine hydrochloride. Unlike the molten globule states of the bovine and goat proteins, the human alpha-lactalbumin molten globule showed remarkably more intense circular dichroism ellipticity than the native state in the far-ultraviolet region below 225 nm. During refolding from the unfolded state, human alpha-lactalbumin thus exhibited overshoot kinetics, in which the alpha-helical peptide ellipticity exceeded the native value when the molten globule folding intermediate was formed in the burst phase. The subsequent folding involved reorganization of nonnative secondary structures. It should be noted that the rate constant of the major refolding phase was approximately the same among the three types of alpha-lactalbumin and that the rate constant of unfolding was accelerated 18-600 times in the human protein, and these results interpreted the lower thermodynamic stability of this protein.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.     

Role of arginine 67 in the stabilization of chymotrypsin inhibitor 2: examination of amide proton exchange rates and denaturation thermodynamics of an engineered protein

We have examined the contribution to protein stability of an interaction involving a charged hydrogen bond from an arginyl side chain (Arg67) in the serine proteinase inhibitor chymotrypsin inhibitor 2 (CI-2), by replacing this side chain with an alanyl residue by protein engineering. Using nuclear magnetic resonance spectroscopy (NMR), we have examined the effect of this mutation on the hydrogen-deuterium exchange rates of several backbone amide protons in the native and engineered proteins at 50 degrees C. These exchange rates provide a localized probe at multiple discrete sites throughout the protein and from comparison of native and mutant exchange rates allow calculation of the difference in free energy of exchange (delta delta Gex) resulting from the mutation. The results show that for the majority of amides observed this mutation results in delta delta Gex of ca. 1.7 kcal mol-1 over the whole CI-2 molecule. However, for two relatively exposed amide protons the exchange rates are found to be far less perturbed, implying that local unfolding mechanisms predominate for these protons. Direct measurement of the stability of both proteins to denaturation by guanidinum hydrochloride shows that the interaction contributes 1.4 kcal mol-1 to the stability of the molecule. This value is comparable to those obtained from the NMR exchange measurements and indicates that the exchange processes reflect the differences in stability between the native and mutant proteins.(ABSTRACT TRUNCATED AT 250 WORDS)     

The native state conformational ensemble of the SH3 domain from alpha-spectrin.

The folding/unfolding equilibrium of the alpha-spectrin SH3 domain has been measured by NMR-detected hydrogen/deuterium exchange and by differential scanning calorimetry. Protection factors against exchange have been obtained under native conditions for more than half of the residues in the domain. Most protected residues are located at the beta-strands, the short 3(10) helix, and part of the long RT loop, whereas the loops connecting secondary structure elements show no measurable protection. Apparent stability constants per residue and their corresponding Gibbs energies have been calculated from the exchange experiments. The most stable region of the SH3 domain is defined by the central portions of the beta-strands. The peptide binding region, on the other hand, is composed of a highly stable region (residues 53-57) and a highly unstable region, the loop between residues 34-41 (n-Src loop). All residues in the domain have apparent Gibbs energies lower than the global unfolding Gibbs energy measured by differential scanning calorimetry, indicating that under our experimental conditions the amide exchange of all residues in the SH3 domain occurs primarily via local unfolding reactions. A structure-based thermodynamic analysis has allowed us to predict correctly the thermodynamics of the global unfolding of the domain and to define the ensemble of conformational states that quantitatively accounts for the observed pattern of hydrogen exchange protection. These results demonstrate that under native conditions the SH3 domain needs to be considered as an ensemble of conformations and that the hydrogen exchange data obtained under those conditions cannot be interpreted by a two-state equilibrium. The observation that specific regions of a protein are able to undergo independent local folding/unfolding reactions indicates that under native conditions the scale of cooperative interactions is regional rather than global.     

A statistical mechanical model for hydrogen exchange in globular proteins.

We develop a statistical mechanical theory for the mechanism of hydrogen exchange in globular proteins. Using the HP lattice model, we explore how the solvent accessibilities of chain monomers vary as proteins fluctuate from their stable native conformations. The model explains why hydrogen exchange appears to involve two mechanisms under different conditions of protein stability; (1) a “global unfolding” mechanism by which all protons exchange at a similar rate, approaching that of the denatured protein, and (2) a “stable-state” mechanism by which protons exchange at rates that can differ by many orders of magnitude. There has been some controversy about the stable-state mechanism: does exchange take place inside the protein by solvent penetration, or outside the protein by the local unfolding of a subregion? The present model indicates that the stable-state mechanism of exchange occurs through an ensemble of conformations, some of which may bear very little resemblance to the native structure. Although most fluctuations are small-amplitude motions involving solvent penetration or local unfolding, other fluctuations (the conformational distant relatives) can involve much larger transient excursions to completely different chain folds.     

Role of partial protein unfolding in alcohol‐induced protein aggregation

Proteins aggregate in response to various stresses including changes in solvent conditions. Addition of alcohols has been recently shown to induce aggregation of disease‐related as well as nondisease‐related proteins. Here we probed the biophysical mechanisms underlying alcohol‐induced protein aggregation, in particular the role of partial protein unfolding in aggregation. We have studied aggregation mechanisms due to benzyl alcohol which is used in numerous biochemical and biotechnological applications. We chose cytochrome c as a model protein, for the reason that various optical and structural probes are available to monitor its global and partial unfolding reactions. Benzyl alcohol induced the aggregation of cytochrome c in isothermal conditions and decreased the temperature at which the protein aggregates. However, benzyl alcohol did not perturb the overall native conformation of cytochrome c. Instead, it caused partial unfolding of a local protein region around the methionine residue at position 80. Site‐specific optical probes, two‐dimensional NMR titrations, and hydrogen exchange all support this conclusion. The protein aggregation temperature varied linearly with the melting temperature of the Met80 region. Stabilizing the Met80 region by heme iron reduction drastically decreased protein aggregation, which confirmed that the local unfolding of this region causes protein aggregation. These results indicate that a possible mechanism by which alcohols induce protein aggregation is through partial rather than complete unfolding of native proteins. Proteins 2010. © Wiley‐Liss, Inc.     

NMR hydrogen exchange of the OB-fold protein LysN as a function of denaturant: the most conserved elements of structure are the most stable to unfolding.

The structure of LysN contains an OB-fold motif composed of a structurally conserved five-stranded beta-barrel capped by a poorly conserved alpha-helix between strands beta3 and beta4. Two additional alpha-helices, unique to the LysN structure, flank the N terminus of the OB-fold. The stability of LysN to unfolding has been investigated with NMR native state hydrogen exchange measurements as a function of guanidinium hydrochloride concentration, and equilibrium unfolding transitions monitored by ellipticity at 222 nm and fluorescence at 350 nm. The spectrophotometric measurements suggest an apparent two-state unfolding transition with DeltaGu(0) approximately 6 kcal/mol and m approximately 3 kcal/(molM). By contrast, NMR hydrogen exchange measurements manifest a distribution of DeltaGu(0) and m values which indicate that the protein can undergo subglobal unfolding. The largest DeltaGu(0) values from hydrogen exchange are for residues in the beta-sheet of the protein. These values, which reflect complete unfolding of the protein, are between 3 and 4 kcal/mol higher than those obtained from circular dichroism or fluorescence. This discrepancy may be due to the comparison of NMR hydrogen exchange parameters measured at residue-level resolution, with spectrophotometric parameters that reflect an unresolved super position of unfolding transitions of the alpha-helices and beta-strands. The largest DeltaGu(0) values obtained from hydrogen exchange for the subset of residues in the alpha-helices of the protein, agree with the DeltaGu(0) values obtained from circular dichroism or fluorescence. Based on the hydrogen exchange data, however, the three alpha-helices of LysN are on average 3 kcal/mol less stable than the beta-sheet. Consistent with the subglobal unfolding of LysN evinced by hydrogen exchange, a deletion mutant that lacks the first alpha-helix of the protein retains a cooperatively folded structure. Taken together with previous results on the OB-fold proteins SN and CspA, the present results for LysN suggest that the most conserved elements of structure in the OB-fold motif are the most resistant to denaturation. In all three proteins, stability to denaturation correlates with sequence hydrophobicity.     

The mechanisms of the proteolytic degradation of native globular proteins. The role of local and global fluctuations of the native structure

Analysis of the proteolytic degradation of the native protein structure carried out by the comparison of the temperature dependence of the hydrogen exchange and proteolytic splitting rates of the hen egg-white lysozyme and human Hb and apoHb. Acceleration of the burst-like (all or none) proteolytic degradation in the high temperature range is provided by the intensification of the global fluctuations with overall unfolding revealed by hydrogen exchange. For Hb and apoHb the rate of burst-like proteolytic degradation and hydrogen exchange weakly depends on temperature in the range, where hydrogen exchange reveals only local fluctuations of the native protein structure. The splitting of the two proteins proceeds by the selfaccelerated burst-like mechanism with the initial rate-limiting single cleavage owing to the local fluctuation of the native structure. The local fluctuations play important role also upon the intracellular burst-like degradation of native proteins.     

Thermal-induced unfolding domains in aldolase identified by amide hydrogen exchange and mass spectrometry.

Amide hydrogen exchange has been measured in short segments of intact rabbit muscle aldolase at temperatures of 14-50 degrees C by the protein fragmentation/mass spectrometry method (Zhang Z, Smith DL, 1993, Protein Sci 2:522-531). Deuterium levels in some segments did not change over the temperature range of the measurements, whereas deuterium levels in other segments increased rapidly with temperature. These results demonstrate that the equilibrium constant for local unfolding, Kunf, of some segments increases with temperature in the low temperature range (14-30 degrees C) of this study. Aldolase begins to lose activity at temperatures above 40 degrees C. In the 40-50 degrees C temperature range, Kunf is greater than 10(-4) in some regions and less than 10(-6) in other regions. This wide range of regional stability in the temperature range where aldolase begins to denature is interpreted in terms of cooperative unfolding/folding domains. Regions of highest stability were located along the hydrophobic subunit binding surface. It is proposed that hydrogen exchange might be used to identify unfolding domains in multidomain proteins whose thermodynamic properties have been determined by differential scanning calorimetry.     

Proteolysis as a probe of thermal unfolding of cytochrome C

Recent hydrogen exchange experiments on native cytochrome c implicate a sequential unfolding pathway in contrast to a simple two-state process. We have studied the heat-induced unfolding of this protein by using spectroscopic measurements to detect changes in conformation and proteolytic enzyme digestion to identify regions of the protein that are labile. Several spectroscopic profiles were monitored: CD at 222 nm, a measurement of secondary structure change in the protein, the absorbance at 280 nm, involving the local environment of Trp 59, and absorbance at 420 nm, the Soret band of the heme. The apparent T_m values for these probes differ, consistent with an unfolding pathway containing intermediates. The limited digestion by proteinase K is consistent with population of an intermediate state in unfolding. We find a single strong region of cleavage at low temperature with retention of structure in each fragment. Proteins 30:435–441, 1998. © 1998 Wiley-Liss, Inc.     

Investigation of apomyoglobin stability depending on urea and temperature at two different pH values

Equilibrium unfolding of apomyoglobin by urea was investigated in the temperature range from 5 to 25 degrees C at two pH values. The thermodynamic parameters of the apomyoglobin native-unfolded state transition were determined. Conformational changes in the protein structure were monitored by tryptophan fluorescence and far UV circular dichroism. Apomyoglobin preserves its native conformation at pH 5.7 and 6.2 in the temperature range used. It was shown that the apomyoglobin stability and its unfolding cooperativity are substantially lower at 5 degrees C than at other temperatures. This fact should be taken in account at the investigation of apomyoglobin.     

Equilibrium hydrogen exchange reveals extensive hydrogen bonded secondary structure in the on-pathway intermediate of Im7

The four-helical immunity protein Im7 folds through an on-pathway intermediate that has a specific, but partially misfolded, hydrophobic core. In order to gain further insight into the structure of this species, we have identified the backbone hydrogen bonds formed in the ensemble by measuring the amide exchange rates (under EX2 conditions) of the wild-type protein and a variant, I72V. In this mutant the intermediate is significantly destabilised relative to the unfolded state (deltadeltaG(ui) = 4.4 kJ/mol) but the native state is only slightly destabilised (deltadeltaG(nu) = 1.8 kJ/mol) at 10 degrees C in 2H2O, pH* 7.0 containing 0.4 M Na2SO4, consistent with the view that this residue forms significant non-native stabilising interactions in the intermediate state. Comparison of the hydrogen exchange rates of the two proteins, therefore, enables the state from which hydrogen exchange occurs to be identified. The data show that amides in helices I, II and IV in both proteins exchange slowly with a free energy similar to that associated with global unfolding, suggesting that these helices form highly protected hydrogen-bonded helical structure in the intermediate. By contrast, amides in helix III exchange rapidly in both proteins. Importantly, the rate of exchange of amides in helix III are slowed substantially in the Im7* variant, I72V, compared with the wild-type protein, whilst other amides exchange more rapidly in the mutant protein, in accord with the kinetics of folding/unfolding measured using chevron analysis. These data demonstrate, therefore, that local fluctuations do not dominate the exchange mechanism and confirm that helix III does not form stable secondary structure in the intermediate. By combining these results with previously obtained Phi-values, we show that the on-pathway folding intermediate of Im7 contains extensive, stable hydrogen-bonded structure in helices I, II and IV, and that this structure is stabilised by both native and non-native interactions involving amino acid side-chains in these helices.     

Aggregation of granulocyte colony stimulating factor under physiological conditions: characterization and thermodynamic inhibition

We have investigated the aggregation of recombinant human granulocyte colony stimulating factor (rhGCSF), a protein that rapidly aggregates and precipitates at pH 6.9 and 37 degrees C. We observed that native monomeric rhGCSF reversibly forms a dimer under physiological conditions and that this dimeric species does not participate in the irreversible aggregation process. Sucrose, a thermodynamic stabilizer, inhibits the aggregation of rhGCSF. We postulate that sucrose acts by reducing the concentration of structurally expanded species, consistent with the hypothesis that preferential exclusion favors most compact species in the native state ensemble. Thermodynamic stability data from unfolding curves and hydrogen-deuterium exchange experimental results support the above hypothesis. Thus, the strategy of stabilizing the native state of the protein under physiological conditions using thermodynamic stabilizers, especially ligands binding with high affinity to the native state, is expected to protect against protein aggregation occurring under such nonperturbing solution conditions.     

Activity-stability relationships in extremophilic enzymes

Psychrophilic, mesophilic, and thermophilic alpha-amylases have been studied as regards their conformational stability, heat inactivation, irreversible unfolding, activation parameters of the reaction, properties of the enzyme in complex with a transition state analog, and structural permeability. These data allowed us to propose an energy landscape for a family of extremophilic enzymes based on the folding funnel model, integrating the main differences in conformational energy, cooperativity of protein unfolding, and temperature dependence of the activity. In particular, the shape of the funnel bottom, which depicts the stability of the native state ensemble, also accounts for the thermodynamic parameters of activation that characterize these extremophilic enzymes, therefore providing a rational basis for stability-activity relationships in protein adaptation to extreme temperatures.     

Mapping the lifetimes of local opening events in a native state protein.

The rate constants for the processes that lead to local opening and closing of the structures around hydrogen bonds in native proteins have been determined for most of the secondary structure hydrogen bonds in the four-helix protein acyl coenzyme A binding protein. In an analysis that combines these results with the energies of activation of the opening processes and the stability of the local structures, three groups of residues in the protein structure have been identified. In one group, the structures around the hydrogen bonds have frequent openings, every 600 to 1,500 s, and long lifetimes in the open state, around 1 s. In another group of local structures, the local opening is a very rare event that takes place only every 15 to 60 h. For these the lifetime in the open state is also around 1 s. The majority of local structures have lifetimes between 2,000 and 20,000 s and relatively short lifetimes of the open state in the range between 30 and 400 ms. Mapping of these groups of amides to the tertiary structure shows that the openings of the local structures are not cooperative at native conditions, and they rarely if ever lead to global unfolding. The results suggest a mechanism of hydrogen exchange by progressive local openings.