pH dependence of the hydrogen exchange in the SH3 domain of alpha-spectrin

Using nuclear magnetic resonance we have measured the hydrogen exchange (HX) in the Src homology region 3 (SH3) domain of alpha-spectrin as a function of pH*. At very acidic pH* values the exchange of most residues appears to occur via global unfolding, although several residues show abnormally large Gibbs energies of exchange, suggesting the presence of some residual structure in the unfolded state. At higher pH* HX occurs mainly via local or partial unfoldings. We have been able to characterize the coupling between the electrostatic interactions in this domain and the conformational fluctuations occurring under native conditions by analyzing the dependence upon pH* of the Gibbs energy of exchange. The SH3 domain seems to be composed of a central core, which requires large structural disruptions to become exposed to the solvent, surrounded by smaller subdomains, which fluctuate independently.     

GuHCl and NaCl-dependent hydrogen exchange in MerP reveals a well-defined core with an unusual exchange pattern

We have analysed hydrogen exchange at amide groups to characterise the energy landscape of the 72 amino acid residue protein MerP. From the guanidine hydrochloride (GuHCl) dependence of exchange in the pre-transitional region we have determined free energy values of exchange (DeltaG(HX)) and corresponding m-values for individual amide protons. Detailed analysis of the exchange patterns indicates that for one set of amide protons there is a weak dependence on denaturant, indicating that the exchange is dominated by local fluctuations. For another set of amide protons a linear, but much stronger, denaturant dependence is observed. Notably, the plots of free energy of exchange versus [GuHCl] for 16 amide protons show pronounced upward curvature, and a close inspection of the structure shows that these residues form a well-defined core in the protein. The hydrogen exchange that was measured at various concentrations of NaCl shows an apparent selective stabilisation of this core. Detailed analysis of this exchange pattern indicates that it may originate from selective destabilisation of the unfolded state by guanidinium ions and/or selective stabilisation of the core in the native state by chloride ions.     

On the nature of conformational openings: native and unfolded-state hydrogen and thiol-disulfide exchange studies of ferric aquomyoglobin.

Native-state amide hydrogen exchange (HX) of proteins in the presence of denaturant has provided valuable details on the structures of equilibrium folding intermediates. Here, we extend HX theory to model thiol group exchange (SX) in single cysteine-containing variants of sperm whale ferric aquomyoglobin. SX is complementary to HX in that it monitors conformational opening events that expose side-chains, rather than the main chain, to solvent. A simple two-process model, consisting of EX2-limited local structural fluctuations and EX1-limited global unfolding, adequately accounts for all HX data. SX is described by the same model except at very low denaturant concentrations and when the bulky labeling reagent 5,5'-dithiobis (2-nitrobenzoic acid) (DTNB) is used. Under these conditions SX can occur by a novel denaturant-dependent process. This anomalous behavior is not observed when the smaller labeling reagent methyl methanethiosulfonate is employed, suggesting that it reflects a denaturant-induced increase in the amplitudes of local structural fluctuations. It also is not seen in heme-free apomyoglobin, which may indicate that local openings are sufficiently large in the absence of denaturant to allow DTNB unhindered access. Differences in SX kinetics obtained using the two labeling reagents provide estimates of the sizes of local opening reactions at different sites in the protein. At all sequence positions examined except for position 73, the same opening event appears to facilitate exchange of both backbone amide and side-chain thiol groups. The C73 thiol group is exposed by a low-energy fluctuation that does not expose its amide group to exchange.     

Energetics of the native energy landscape of a two-domain calcium sensor protein: distinct folding features of the two domains

The protein folding energy landscape allows a thorough understanding of the protein folding problem which in turn helps in understanding various aspects of biological functions. Characterizing the cooperative unfolding units and the intermediates along the folding funnel of a protein is a challenging task. In this paper, we investigated the native energy landscape of EhCaBP, a calcium sensor, belonging to the same EF-hand superfamily as calmodulin. EhCaBP is a two-domain EF-hand protein consisting of two EF-hands in each domain and binding to four Ca2+ cations. Native-state hydrogen exchange (HX) was used to assess the folding features of the landscape and also to throw light on the structure-folding function paradigm of calcium sensor proteins. HX measurements under the EX2 regime provided the thermodynamic information about the protein folding events under native conditions. HX studies revealed that the unfolding of EhCaBP is not a two-state process. Instead, it proceeds through cooperative units. The C-terminal domain exhibits less denaturant dependence than the N-terminal domain, suggesting that the former is dominated by local fluctuations. It is interesting to note that the N- and C-terminal domains of EhCaBP have distinct folding features. In fact, these observed differences can regulate the domain-dependent target recognition of two-domain Ca2+ sensor proteins.     

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.     

The temperature dependence of the hydrogen exchange in the SH3 domain of alpha-spectrin

The amide hydrogen-deuterium exchange (HX) in the Src homology region 3 (SH3) domain of alpha-spectrin has been measured by nuclear magnetic resonance as a function of temperature between 8 and 46 degrees C. The analysis of the temperature dependence of HX from a statistical thermodynamic point of view has allowed us to estimate the enthalpies and entropies of the conformational processes leading to HX. The results indicate that under native conditions the domain undergoes a wide variety of conformational fluctuations, ranging from local motions, mainly located in loops, turns and chain ends and involving only low enthalpy and entropy, to extensive structural disruptions affecting its core and involving enthalpies and entropies that come fairly close to those observed during global unfolding.     

Molecular mechanisms of pH-driven conformational transitions of proteins: insights from continuum electrostatics calculations of acid unfolding

The acid unfolding of staphylococcal nuclease (SNase) is very cooperative (Whitten and García-Moreno, Biochemistry 2000;39:14292-14304). As many as seven hydrogen ions (H+) are bound preferentially by the acid-unfolded state relative to the native (N) state in the pH range 3.2-3.9. To investigate the mechanism of acid unfolding, structure-based pKa calculations were performed with a variety of continuum electrostatic methods. The calculations reproduced successfully the H+ binding properties of the N state between pH 5 and 9, but they systematically overestimated the number of H+ bound upon acid unfolding. The calculated pKa values of all carboxylic residues in the N state were more depressed than they should be. The discrepancy between the observed and the calculated H+ uptake upon acid unfolding was not improved by using high protein dielectric constants, structures relaxed with molecular dynamics, or other empirical modifications implemented previously by others to maximize agreement between measured and calculated pKa values. This suggests an important role for conformational fluctuations of the backbone as important determinants of pKa values of carboxylic groups. Because no global or subglobal conformational changes have been observed previously for SNase under acidic conditions above the acid-unfolding region, these fluctuations must be local. The acid unfolding of SNase does not seem to involve the disruption of the N state by accruement of intramolecular repulsive interactions, nor the protonation of key ion paired carboxylic residues. It is more consistent with modest contributions from many H+ binding groups, with an important role for local conformational fluctuations in the coupling between H+ binding and the global structural transition.     

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.     

Effects of sucrose on conformational equilibria and fluctuations within the native-state ensemble of proteins

Osmolytes increase the thermodynamic conformational stability of proteins, shifting the equilibrium between native and denatured states to favor the native state. However, their effects on conformational equilibria within native-state ensembles of proteins remain controversial. We investigated the effects of sucrose, a model osmolyte, on conformational equilibria and fluctuations within the native-state ensembles of bovine pancreatic ribonuclease A and S and horse heart cytochrome c. In the presence of sucrose, the far- and near-UV circular dichroism spectra of all three native proteins were slightly altered and indicated that the sugar shifted the native-state ensemble toward species with more ordered, compact conformations, without detectable changes in secondary structural contents. Thermodynamic stability of the proteins, as measured by guanidine HCl-induced unfolding, increased in proportion to sucrose concentration. Native-state hydrogen exchange (HX) studies monitored by infrared spectroscopy showed that addition of 1 M sucrose reduced average HX rate constants at all degrees of exchange of the proteins, for which comparison could be made in the presence and absence of sucrose. Sucrose also increased the exchange-resistant core regions of the proteins. A coupling factor analysis relating the free energy of HX to the free energy of unfolding showed that sucrose had greater effects on large-scale than on small-scale fluctuations. These results indicate that the presence of sucrose shifts the conformational equilibria toward the most compact protein species within native-state ensembles, which can be explained by preferential exclusion of sucrose from the protein surface.     

pH dependence of individual tryptophan N-1 hydrogen exchange rates in lysozyme and its chemically modified derivatives

Nuclear magnetic resonance analyses have been made of the individual hydrogen-deuterium exchange rates of tryptophan indole N-1 hydrogens in native lysozyme and its chemically modified derivatives including lysozyme with an ester cross-linkage between Glu-35 and Trp-108, lysozyme with an internal amide cross-linking between the epsilon-amino group of Lys-13 and the alpha-carboxyl group of Leu-129, and lysozyme with the beta-aspartyl sequence at Asp-101. The pH dependence curves of the exchange rates for Trp-63 and Trp-108 are different from those expected for tryptophan. The pH dependence curve for Trp-108 exchange exhibits the effects from molecular aggregation at pH above 5 and from a transition between the two conformational fluctuations at around pH 4. The exchange rates for tryptophan residues in native lysozyme and modified derivatives are not correlated with the thermodynamic or kinetic parameters in protein denaturation, suggesting that the fluctuations responsible for the exchange are not global ones. The exchange rates for tryptophan residues remote from the modification site are perturbed. Such tryptophan residues are found to be involved in a small but distinct conformational change due to the modification. Therefore, the perturbations of the N-1 hydrogen exchange rates are related to the minor change in local conformation or in conformational strain induced by the chemical modification.     

Restricted backbone conformational and motional flexibilities of loops containing peptidyl-proline bonds dominate the enzyme activity of staphylococcal nuclease

The role of cis-trans isomerizations of peptidyl-proline bonds in the enzyme activity of staphylococcal nuclease (SNase) was examined by mutation of proline residues. The proline-free SNase ([Pro-]SNase), namely, P11A/P31A/P42A/P47T/P56A/P117G-mutant SNase, was adopted for elucidating the correlation between the nuclease activity and the backbone conformational and dynamic states of SNase. The 3D solution structure of [Pro-]SNase has been determined by heteronuclear NMR experiments. Comparing the structure of [Pro-]SNase with the structure of SNase revealed the conformational differences between the two proteins. In the structure of [Pro-]SNase, conformational rearrangements were observed for the loop of residues Ala112-His121 containing a trans Lys116-Gly117 peptide bond and for the C-terminal alpha-helical loop of residues Leu137-Glu142. Mutation of proline at position 117 also caused the conformational rearrangement of the p-loop (Asp77-Leu89), which is remote from the Ala112-His121 loop. The Ala112-His121 loop and p-loop are placed closer to each other in [Pro-]SNase than in SNase. The backbone dynamic features of the omega-loop (Pro42-Pro56) of SNase are different from those of [Pro-]SNase. The backbone of the omega-loop exhibits restricted flexibility with slow conformational exchange motions in SNase, but is highly flexible in [Pro-]SNase. The analysis indicates that the restrained backbone conformation of the Ala112-His121 loop and restricted flexibility of the omega-loop are two dominant factors determining the enzyme activity of SNase. Of the two factors, the former is correlated with the strained cis Lys116-Pro117 peptide bond and the latter is correlated with the cis-trans isomerizations of the His46-Pro47 peptide bond.     

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.     

Probing the high energy states in proteins by proteolysis

Unless the native conformation has an unstructured region, proteases cannot effectively digest a protein under native conditions. Digestion must occur from a higher energy form, when at least some part of the protein is exposed to solvent and becomes accessible by proteases. Monitoring the kinetics and denaturant dependence of proteolysis under native conditions yields insight into the mechanism of proteolysis as well as these high-energy conformations. We propose here a generalized approach to exploit proteolysis as a tool to probe high-energy states in proteins. This "native state proteolysis" experiment was carried out on Escherichia coli ribonuclease HI. Mass spectrometry and N-terminal sequencing showed that thermolysin cleaves the peptide bond between Thr92 and Ala93 in an extended loop region of the protein. By comparing the proteolysis rate of the folded protein and a peptidic substrate mimicking the sequence at the cleavage site, the energy required to reach the susceptible state (Delta G(proteolysis)) was determined. From the denaturant dependence of Delta G(proteolysis), we determined that thermolysin digests this protein through a local fluctuation, i.e. localized unfolding with minimal change in solvent assessable surface area. Proteolytic susceptibilities of proteins are discussed based on the finding of this local fluctuation mechanism for proteolysis under native conditions.     

Scope and utility of hydrogen exchange as a tool for mapping landscapes

The ability to determine conformational parameters of protein-folding landscapes is critical for understanding the link between conformation, function, and disease. Monitoring hydrogen exchange (HX) of labile protons at equilibrium enables direct extraction of thermodynamic or kinetic landscape parameters in two limiting extremes. Here, we establish a quantitative framework for relating HX behavior to landscape. We use this framework to demonstrate that the range of predicted global HX behavior for the majority of a set of characterized two-state proteins under near-native conditions does not readily span between both extremes. For most, stability may be quantitatively determined under physiological conditions, with semiquantitative boundaries on kinetics additionally determined using modest experimental perturbations to shift HX behavior. The framework and relationships derived in the simple context of two-state global folding highlight the importance of understanding HX across the entire continuum of behavior, in order to apply HX to map landscapes.     

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.     

Characterization of millisecond time-scale dynamics in the molten globule state of alpha-lactalbumin by NMR

The motional dynamics of the molten globule (MG) state of alpha-lactalbumin have been characterized using (15)N transverse relaxation rates (R2). A modified version of the Carr-Purcell-Meiboom-Gill (CPMG) R2 pulse sequence is proposed in order to overcome the loss of sensitivity that arises from extreme line broadening due to complex dynamics on the millisecond time-scale. Using this pulse sequence, chemical exchange rates were extracted by examining the (15)N transverse relaxation rates as a function of CPMG delay values. The results clearly illustrate that pervasive conformational exchange of 0.2-0.5 ms in the (15)N backbone resonances of the molten globule state of alpha-lactalbumin. The temperature dependence of the conformational exchange rates display standard Arrhenius kinetic behavior between 10 and 30 degrees C. Estimates of the activation energies range from 0.8 to 4. 4 kcal/mol, indicating a low energetic barrier to conformational fluctuations relative to native state proteins. The fluctuations and low energetic barriers may be critical for directing the search for contacts that will result in the transition from the MG state to the native state.     

Hydrogen exchange kinetics of RNase A and the urea:TMAO paradigm

A key paradigm in the biology of adaptation holds that urea affects protein function by increasing the fluctuations of the native state, while trimethylamine N-oxide (TMAO) affects function in the opposite direction by decreasing the normal fluctuations of the native ensemble. Using urea and TMAO separately and together, hydrogen exchange (HX) studies on RNase A at pH* 6.35 were used to investigate the basic tenets of the urea:TMAO paradigm. TMAO (1 M) alone decreases HX rate constants of a select number of sites exchanging from the native ensemble, and low urea alone increases the rate constants of some of the same sites. Addition of TMAO to urea solutions containing RNase A also suppresses HX rate constants. The data show that urea and TMAO independently or in combination affect the dynamics of the native ensemble in opposing ways. The results provide evidence in support of the counteraction aspect of the urea:TMAO paradigm linking structural dynamics with protein function in urea-rich organs and organisms. RNase A is so resistant to urea denaturation at pH* 6.35 that even in the presence of 4.8 M urea, the native ensemble accounts for >99.5% of the protein. An essential test, devised to determine the HX mechanism of exchangeable protons, shows that over the 0-4.8 M urea concentration range nearly 80% of all observed sites convert from EX2 to EX1. The slow exchange sites are all EX1; they do not exhibit global exchange even at urea concentrations (5.8 M) well into the denaturation transition zone, and their energetically distinct activated complexes leading to exchange gives evidence of residual structure. Under these experimental conditions, the use of DeltaG(HX) as a basis for HX analysis of RNase A urea denaturation is invalid.     

Definable equilibrium states in the folding of human prion protein

The role of conformational intermediates in the conversion of prion protein from its normal cellular form (PrP(C)) to the disease-associated "scrapie" form (PrP(Sc)) remains unknown. To look for such intermediates in equilibrium conditions, we have examined the unfolding transitions of PrP(C), primarily using the chemical denaturant guanidine hydrochloride (GuHCl). When the protein conformation is assessed by NMR, there is a gradual shift of NMR signals in the regions between residues 125-146 and 186-196. The denaturant dependence of these shifts shows that in aqueous solution the native and locally unfolded conformations are both significantly populated. Following this shift, there is the major unfolding transition to generate a substantially unfolded population. However, analysis of NMR chemical shift and intensity changes shows that there is persistent structure in the molecule well beyond this major cooperative unfolding transition. Residual structure within this state is extensive and encompasses the majority of the secondary structure elements found in the native state of the protein.     

Investigation of the four cooperative unfolding units existing in the MICAL-1 CH domain

The solution structure of human MICAL-1 calpolnin homology (CH) domain is composed of six alpha helices and one 3(10) helix. To study the unfolding of this domain, we carry out native-state hydrogen exchange, intrinsic fluorescence and far-UV circular dichroism experiments. The free energy of unfolding, DeltaG(H2O), is calculated to be 7.11+/-0.58 kcal mol(-1) from GuHCl denaturation at pH 6.5. Four cooperative unfolding units are found using native-state hydrogen exchange experiment. Forty-seven slow-exchange residues can be studied by native-state hydrogen exchange experiments. From the concentration dependence of exchange rates, free energy of amide hydrogen with solvent, DeltaG(HX) and m-value (sensitivity of exposure to denaturant) are obtained, which reveal four cooperative unfolding units. The slowest exchanging protons are distributed throughout the whole hydrophobic core of the protein, which might be the folding core. These results will help us understand the structure of MICAL-1 CH domain more deeply.     

Probing local structural fluctuations in myoglobin by size-dependent thiol-disulfide exchange

All proteins undergo local structural fluctuations (LSFs) or breathing motions. These motions are likely to be important for function but are poorly understood. LSFs were initially defined by amide hydrogen exchange (HX) experiments as opening events, which expose a small number of backbone amides to ¹H/²H exchange, but whose exchange rates are independent of denaturant concentration. Here, we use size-dependent thiol-disulfide exchange (SX) to characterize LSFs in single cysteine-containing variants of myoglobin (Mb). SX complements HX by providing information on motions that disrupt side chain packing interactions. Most importantly, probe reagents of different sizes and chemical properties can be used to characterize the size of structural opening events and the properties of the open state. We use thiosulfonate reagents (126–274 Da) to survey access to Cys residues, which are buried at specific helical packing interfaces in Mb. In each case, the free energy of opening increases linearly with the radius of gyration of the probe reagent. The slope and the intercept are interpreted to yield information on the size of the opening events that expose the buried thiol groups. The slope parameter varies by over 10-fold among Cys positions tested, suggesting that the sizes of breathing motions vary substantially throughout the protein. Our results provide insight to the longstanding question: how rigid or flexible are proteins in their native states?