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?     

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.     

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.     

Amide proton exchange in proteins by EX1 kinetics: studies of the basic pancreatic trypsin inhibitor at variable p2H and temperature

With the use of one-dimensional 1H nuclear magnetic resonance, two-dimensional correlated spectroscopy, and two-dimensional nuclear Overhauser enhancement spectroscopy, the exchange mechanisms for numerous individual amide protons in the basic pancreatic trypsin inhibitor (BPTI) were investigated over a wide range of p2H and temperature. Correlated exchange under an EX1 regime was observed only for the most slowly exchanging protons in the central hydrogen bonds of the antiparallel beta-sheet and only over a narrow range of temperature and p2H, i.e., above ca. 55 degrees C and between p2H 7 and 9, where the opening rates of the structure fluctuations which promote the exchange of these protons are of the order 0.1 min-1. At p2H below 7, the exchange of this most stable group of protons is uncorrelated and is governed by an EX2 mechanism. At p2H above 9, the exchange is also uncorrelated and occurs via either EX2 or EX1 processes promoted by strictly local structure fluctuations. For all other backbone amide protons in BPTI, the exchange was found to be uncorrelated and by an EX2 mechanism under all conditions of p2H and temperature where quantitative measurements could be obtained with the methods used, i.e., for kex approximately less than 5 min-1. From these observations with BPTI it can be concluded that the amide proton exchange in globular proteins is quite generally via EX2 processes, with rare exceptions for measurements with extremely stable protons at high temperature and basic p2H. This emphasizes the need for further development of suitable concepts for the structural interpretation of EX2 amide proton exchange [Wagner, G. (1983) Q. Rev. Biophys. 16, 1-57; Wagner, G., Stassinopoulou, C. I., & Wüthrich, K. (1984) Eur. J. Biochem. 145, 431-436] and for more detailed investigations of the intrinsic exchange rates for solvent-exposed amide protons in the "open" states of a protein [Roder, H., Wagner, G., & Wüthrich, K. (1985) Biochemistry (following paper in this issue)].     

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.     

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.     

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.     

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.     

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.     

Cooperative alpha-helix unfolding in a protein-DNA complex from hydrogen-deuterium exchange

We present experimental evidence for a cooperative unfolding transition of an alpha-helix in the lac repressor headpiece bound to a symmetric variant of the lac operator, as inferred from hydrogen-deuterium (H-D) exchange experiments monitored by NMR spectroscopy. In the EX1 limit, observed exchange rates become pH-independent and exclusively sensitive to local structure fluctuations that expose the amide proton HN to exchange. Close to this regime, we measured decay rates of individual backbone HN signals in D2O, and of their mutual HN-HN NOE by time-resolved two-dimensional (2D) NMR experiments. The data revealed correlated exchange at the center of the lac headpiece recognition helix, Val20-Val23, and suggested that the correlation breaks down at Val24, at the C terminus of the helix. A lower degree of correlation was observed for the exchange of Val9 and Ala10 at the center of helix 1, while no correlation was observed for Val38 and Glu39 at the center of helix 3. We conclude that HN exchange in the recognition helix and, to some extent, in helix 1 is a cooperative event involving the unfolding of these helices, whereas the HN exchange in helix 3 is dominated by random local structure fluctuations.     

Intimate view of a kinetic protein folding intermediate: residue-resolved structure, interactions, stability, folding and unfolding rates, homogeneity

A cytochrome c kinetic folding intermediate was studied by hydrogen exchange (HX) pulse labeling. Advances in the technique and analysis made it possible to define the structured and unstructured regions, equilibrium stability, and kinetic opening and closing rates, all at an amino acid-resolved level. The entire N-terminal and C-terminal helices are formed and docked together at their normal native positions. They fray in both directions from the interaction region, due to a progression in both unfolding and refolding rates, leading to the surprising suggestion that helix propagation may proceed very slowly in the condensed milieu. Several native-like beta turns are formed. Some residues in the segment that will form the native 60s helix are protected but others are not, suggesting energy minimization to some locally non-native conformation in the transient intermediate. All other regions are unprotected, presumably dynamically disordered. The intermediate resembles a partially constructed native state. It is early, on-pathway, and all of the refolding molecules pass through it. These and related results consistently point to distinct, homogeneous, native-like intermediates in a stepwise sequential pathway, guided by the same factors that determine the native structure. Previous pulse labeling efforts have always assumed EX2 exchange during the labeling pulse, often leading to the suggestion of heterogeneous intermediates in alternative parallel pathways. The present work reveals a dominant role for EX1 exchange in the high pH labeling pulse, which will mimic heterogeneous behavior when EX2 exchange is assumed. The general problem of homogeneous versus heterogeneous intermediates and pathways is discussed.     

Protein dynamics viewed by hydrogen exchange

To examine the relationship between protein structural dynamics and measurable hydrogen exchange (HX) data, the detailed exchange behavior of most of the backbone amide hydrogens of Staphylococcal nuclease was compared with that of their neighbors, with their structural environment, and with other information. Results show that H-bonded hydrogens are protected from exchange, with HX rate effectively zero, even when they are directly adjacent to solvent. The transition to exchange competence requires a dynamic structural excursion that removes H-bond protection and allows exposure to solvent HX catalyst. The detailed data often make clear the nature of the dynamic excursion required. These range from whole molecule unfolding, through smaller cooperative unfolding reactions of secondary structural elements, and down to local fluctuations that involve as little as a single peptide group or side chain or water molecule. The particular motion that dominates the exchange of any hydrogen is the one that allows the fastest HX rate. The motion and the rate it produces are determined by surrounding structure and not by nearness to solvent or the strength of the protecting H-bond itself or its acceptor type (main chain, side chain, structurally bound water). Many of these motions occur over time scales that are appropriate for biochemical function.     

Mapping protein energy landscapes with amide hydrogen exchange and mass spectrometry: I. A generalized model for a two-state protein and comparison with experiment

Protein amide hydrogen exchange (HDX) is a convoluted process, whose kinetics is determined by both dynamics of the protein and the intrinsic exchange rate of labile hydrogen atoms fully exposed to solvent. Both processes are influenced by a variety of intrinsic and extrinsic factors. A mathematical formalism initially developed to rationalize exchange kinetics of individual amide hydrogen atoms is now often used to interpret global exchange kinetics (e.g., as measured in HDX MS experiments). One particularly important advantage of HDX MS is direct visualization of various protein states by observing distinct protein ion populations with different levels of isotope labeling under conditions favoring correlated exchange (the so-called EX1 exchange mechanism). However, mildly denaturing conditions often lead to a situation where the overall HDX kinetics cannot be clearly classified as either EX1 or EX2. The goal of this work is to develop a framework for a generalized exchange model that takes into account multiple processes leading to amide hydrogen exchange, and does not require that the exchange proceed strictly via EX1 or EX2 kinetics. To achieve this goal, we use a probabilistic approach that assigns a transition probability and a residual protection to each equilibrium state of the protein. When applied to a small protein chymotrypsin inhibitor 2, the algorithm allows complex HDX patterns observed experimentally to be modeled with remarkably good fidelity. On the basis of the model we are now in a position to begin to extract quantitative dynamic information from convoluted exchange kinetics.     

Internal protein motions, concentrated glycerol, and hydrogen exchange studied in myoglobin

Experiments were carried out to measure the effect of concentrations of glycerol on H-exchange (HX) rates by using myoglobin as a test protein. Concentrated glycerol has only a small slowing effect on the HX kinetics of freely exposed amides, studied in a small molecule model (acetamide). Larger effects occur in structured proteins. The effect of solvent glycerol on different parts of the HX curve of myoglobin was studied by use of a selective "kinetic labeling" approach. Concentrated glycerol exerts an apparently reverse effect on protein H exchange; the faster exchanging "surface" protons are least affected, while the slower and slower amide NH is further slowed by larger and larger factors. These results seem inconsistent with solvent penetration models which generally visualize slower and slower protons as being placed, and undergoing exchange, farther and farther from the solvent-protein interface. On the other hand, the results are as expected for the local unfolding model for protein H exchange since concentrated glycerol is known to stabilize proteins against unfolding. In the local unfolding model, slower exchanging protons are released by way of higher energy and therefore generally larger, unfolding reactions. Larger unfoldings must be more inhibited by the glycerol effect.     

Hydrogen exchange and proteolytic degradation of ribonuclease A. Similarities and distinctions of the kinetic mechanisms

The studies by IR spectroscopy of the temperature dependence of the H-D exchange rate of the RNase A peptide NH atoms permit one to characterize two types of conformation fluctuations, local and global. A comparison with the temperature dependence of the proteolytic degradation rate of RNase A shows that similar in nature fluctuations allow for the H-D exchange of NH atoms and the splitting of peptide bonds of the native protein. In the low temperature region, both processes occur through local fluctuations, by way of the EX2 mechanism, and in the high temperature region, they occur through global fluctuations with the overall denaturation desorganization of the native structure, by way of the EX1 mechanism. The biphasic dependence of the rate of H-D exchange and proteolytic degradation of RNase A on urea concentration is also explained by the combination of local and global fluctuations.     

pH and urea dependence of amide hydrogen-deuterium exchange rates in the beta-trefoil protein hisactophilin.

Amide hydrogen/deuterium exchange rates were measured as a function of pH and urea for 37 slowly exchanging amides in the beta-trefoil protein hisactophilin. The rank order of exchange rates is generally maintained under different solution conditions, and trends in the pH and urea dependence of exchange rates are correlated with the rank order of exchange rates. The observed trends are consistent with the expected behavior for exchange of different amides via global and/or local unfolding. Analysis of the pH dependence of exchange in terms of rate constants for structural opening and closing reveals a wide range of rates in different parts of the hisactophilin structure. The slowest exchanging amides have the slowest opening and closing rates. Many of the slowest exchanging amides are located in trefoil 2, but there are also some slow exchanging amides in trefoils 1 and 3. Slow exchangers tend to be near the interface between the beta-barrel and the beta-hairpin triplet portions of this single-domain structure. The pattern of exchange behaviour in hisactophilin is similar to that observed previously in interleukin-1 beta, indicating that exchange properties may be conserved among beta-trefoil proteins. Comparisons of opening and closing rates in hisactophilin with rates obtained for other proteins reveal clear trends for opening rates; however, trends in closing rates are less apparent, perhaps due to inaccuracies in the values used for intrinsic exchange rates in the data fitting. On the basis of the pH and urea dependence of exchange rates and optical measurements of stability and folding, EX2 is the main exchange mechanism in hisactophilin, but there is also evidence for varying levels of EX1 exchange at low and high pH and high urea concentrations. Equilibrium intermediates in which subglobal portions of structure are cooperatively disrupted are not apparent from analysis of the urea dependence of exchange rates. There is, however, a strong correlation between the Gibbs free energy of opening and the denaturant dependence of opening for all amides, which suggests exchange from a continuum of states with different levels of structure. Intermediates are not very prominent either in equilibrium exchange experiments or in quenched-flow kinetic studies; hence, hisactophilin may not form partially folded states as readily as IL-1 beta and other beta-trefoil proteins.     

The reactivity of the thiol groups of calf thymus deoxyribonucleohistone

The reactivities of the two cysteine thiol groups of calf thymus F3 histone were investigated using 5,5'-dithiobis-[2- nitrobenzoic acid], (DTNB). In isolated histone, both thiol groups were available for reaction. However, analysis of reaction profiles of native deoxyribonucleohistone, (DNH), in various solvent conditions, together with gel electrophoresis studies of DNH modified with DTNB, showed that only one of the thiol groups is normally modified by the reagent. If NaCl is present (above 1.OM) the other thiol group can also be modified. The reactivities of both groups were largely independent of the degree of DNH supercoiling and of the binding of F3 to the DNA.     

Ellman's reagent: 5,5'-dithiobis(2-nitrobenzoic acid)--a reexamination.

Accurate determination of thiol groups by means of Ellman's reagent [5,5′-dithiobis(2-nitrobenzoic acid), DTNB] has been limited by uncertainty about the molar absorption coefficient of the dianion of the product, 2-nitro-5-thiobenzoic acid (TNB). A procedure is described for the purification of TNB by reduction of commercial DTNB followed by gel chromatography and crystallization. Pure DTNB is prepared by reoxidation of purified TNB followed by gel chromatography and crystallization. The molar absorption coefficient of the dianion of TNB is 14,150 at 412 nm in dilute aqueous salt solutions. This value was confirmed independently by reduction of purified DTNB with cysteine. Titration of sulphydryl groups with DTNB can be done at pH 7.27 where the thiol group of TNB is 99.8% in the intensely-colored conjugate base form while the hydroxide-promoted hydrolysis of DTNB is minimal.     

Interactions between aspartate animotransferase protomers observed during holoenzyme reconstitution..

Aspartate aminotransferase from pig heart cytosol consists of 2 identical protomers. Car?ymethylation of the easily accessible cysteines 45 and 82 of the holoenzyme produces a modified protein which however is still fully active. If the apoenzyme obtained from the car?ymethylated holoenzyme is subjected to the thiol specific reagent DTNB, another group per protomer can be titrated. There is some evidence that the group modified is the thiol 390 which is titrated in “syncatalytic” conditions, i.e. 2 mM ketoglutarate and 70 mM glutamate. The full reactivation of the car?ymethylated apoenzyme is obtained on binding of one PLP at each active site of the dimer. During the reactivation the number of thiols 390 accessible to DTNB decreases from 2 to 0 when PLP is added. The variation of the number of reactive thiol residues is not linear with coenzyme binding. Apo-holoenzyme car?ymethylated hybrid does not react with DTNB. Thus, a strong interaction between the 2 protomers is demonstrated.