Protein hydrogen exchange: testing current models

To investigate the determinants of protein hydrogen exchange (HX), HX rates of most of the backbone amide hydrogens of Staphylococcal nuclease were measured by NMR methods. A modified analysis was used to improve accuracy for the faster hydrogens. HX rates of both near surface and well buried hydrogens are spread over more than 7 orders of magnitude. These results were compared with previous hypotheses for HX rate determination. Contrary to a common assumption, proximity to the surface of the native protein does not usually produce fast exchange. The slow HX rates for unprotected surface hydrogens are not well explained by local electrostatic field. The ability of buried hydrogens to exchange is not explained by a solvent penetration mechanism. The exchange rates of structurally protected hydrogens are not well predicted by algorithms that depend only on local interactions or only on transient unfolding reactions. These observations identify some of the present difficulties of HX rate prediction and suggest the need for returning to a detailed hydrogen by hydrogen analysis to examine the bases of structure-rate relationships, as described in the companion paper (Skinner et al., Protein Sci 2012;21:996-1005).     

Dynamics of the Tec‐family tyrosine kinase SH3 domains

The Src Homology 3 (SH3) domain is an important regulatory domain found in many signaling proteins. X‐ray crystallography and NMR structures of SH3 domains are generally conserved but other studies indicate that protein flexibility and dynamics are not. We previously reported that based on hydrogen exchange mass spectrometry (HX MS) studies, there is variable flexibility and dynamics among the SH3 domains of the Src‐family tyrosine kinases and related proteins. Here we have extended our studies to the SH3 domains of the Tec family tyrosine kinases (Itk, Btk, Tec, Txk, Bmx). The SH3 domains of members of this family augment the variety in dynamics observed in previous SH3 domains. Txk and Bmx SH3 were found to be highly dynamic in solution by HX MS and Bmx was unstructured by NMR. Itk and Btk SH3 underwent a clear EX1 cooperative unfolding event, which was localized using pepsin digestion and mass spectrometry after hydrogen exchange labeling. The unfolding was localized to peptide regions that had been previously identified in the Src‐family and related protein SH3 domains, yet the kinetics of unfolding were not. Sequence alignment does not provide an easy explanation for the observed dynamics behavior, yet the similarity of location of EX1 unfolding suggests that higher‐order structural properties may play a role. While the exact reason for such dynamics is not clear, such motions can be exploited in intra‐ and intermolecular binding assays of proteins containing the domains.     

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.     

On the Solution Conformation and Dynamics of the HIV-1 Viral Infectivity Factor

Human immunodeficiency virus-1 (HIV-1) has evolved a cunning mechanism to circumvent the antiviral activity of the APOBEC3 family of host cell enzymes. HIV-1 Vif [viral (also called virion) infectivity factor], one of several HIV accessory proteins, targets APOBEC3 proteins for proteasomal degradation and downregulates their expression at the mRNA level. Despite the importance of Vif for HIV-1 infection, there is little conformational data on Vif alone or in complex with other cellular factors due to incompatibilities with many structural techniques and difficulties in producing suitable quantities of the protein for biophysical analysis. As an alternative, we have turned to hydrogen exchange mass spectrometry (HX MS), a conformational analysis method that is well suited for proteins that are difficult to study using X-ray crystallography and/or NMR. HX MS was used to probe the solution conformation of recombinant full-length HIV-1 Vif. Vif specifically interacted with the previously identified binding partner Hck and was able to cause kinase activation, suggesting that the Vif studied by HX MS retained a biochemically competent conformation relevant to Hck interaction. HX MS analysis of Vif alone revealed low deuteration levels in the N-terminal portion, indicating that this region contained structured or otherwise protected elements. In contrast, high deuteration levels in the C-terminal portion of Vif indicated that this region was likely unstructured in the absence of cellular interacting proteins. Several regions within Vif displayed conformational heterogeneity in solution, including the APOBEC3G/F binding site and the HCCH zinc finger. Taken together, these HX MS results provide new insights into the solution conformation of Vif.     

The foldon substructure of staphylococcal nuclease

To search for submolecular foldon units, the spontaneous reversible unfolding and refolding of staphylococcal nuclease under native conditions was studied by a kinetic native-state hydrogen exchange (HX) method. As for other proteins, it appears that staphylococcal nuclease is designed as an assembly of well-integrated foldon units that may define steps in its folding pathway and may regulate some other functional properties. The HX results identify 34 amide hydrogens that exchange with solvent hydrogens under native conditions by way of large transient unfolding reactions. The HX data for each hydrogen measure the equilibrium stability (ΔG_HX) and the kinetic unfolding and refolding rates (k_op and k_cl) of the unfolding reaction that exposes it to exchange. These parameters separate the 34 identified residues into three distinct HX groupings. Two correspond to clearly defined structural units in the native protein, termed the blue and red foldons. The remaining HX grouping contains residues, not well separated by their HX parameters alone, that represent two other distinct structural units in the native protein, termed the green and yellow foldons. Among these four sets, a last unfolding foldon (blue) unfolds with a rate constant of 6 × 10^− 6 s^− 1 and free energy equal to the protein's global stability (10.0 kcal/mol). It represents part of the β-barrel, including mutually H-bonding residues in the β4 and β5 strands, a part of the β3 strand that H-bonds to β5, and residues at the N-terminus of the α2 helix that is capped by β5. A second foldon (green), which unfolds and refolds more rapidly and at slightly lower free energy, includes residues that define the rest of the native α2 helix and its C-terminal cap. A third foldon (yellow) defines the mutually H-bonded β1-β2-β3 meander, completing the native β-barrel, plus an adjacent part of the α1 helix. A final foldon (red) includes residues on remaining segments that are distant in sequence but nearly adjacent in the native protein. Although the structure of the partially unfolded forms closely mimics the native organization, four residues indicate the presence of some nonnative misfolding interactions. Because the unfolding parameters of many other residues are not determined, it seems likely that the concerted foldon units are more extensive than is shown by the 34 residues actually observed.     

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.     

Measuring rapid hydrogen exchange in the homodimeric 36 kDa HIV-1 integrase catalytic core domain

Measurements of rapid hydrogen exchange (HX) of water with protein amide sites contain valuable information on protein structure and function, but current NMR methods for measuring HX rates are limited in their applicability to large protein systems. An alternate method for measuring rapid HX is presented that is well-suited for larger proteins, and we apply the method to the deuterated, homodimeric 36 kDa HIV-1 integrase catalytic core domain (CCD). Using long mixing times for water-amide magnetization exchange at multiple pH values, HX rates spanning more than four orders of magnitude were measured, as well as NOE cross-relaxation rates to nearby exchangeable protons. HX protection factors for the CCD are found to be large (>10(4)) for residues along the dimer interface, but much smaller in many other regions. Notably, the catalytic helix (residues 152-167) exhibits low HX protection at both ends, indicative of fraying at both termini as opposed to just the N-terminal end, as originally thought. Residues in the LEDGF/p75 binding pocket also show marginal stability, with protection factors in the 10-100 range (～1.4-2.7 kcal/mol). Additionally, elevated NOE cross-relaxation rates are identified and, as expected, correspond to proximity of the amide proton to a rapidly exchanging proton, typically from an OH side chain. Indirect NOE transfer between H(2) O and the amide proton of I141, a residue in the partially disordered active site of the enzyme, suggests its proximity to the side chain of S147, an interaction seen in the DNA-bound form for a homologous integrase.     

Simultaneous monitoring of peptide aggregate distributions, structure, and kinetics using amide hydrogen exchange: application to Abeta(1-40) fibrillogenesis

Increasing evidence indicates that soluble aggregates of amyloid beta protein (Abeta) are neurotoxic. However, difficulty in isolating these unstable, dynamic species impedes studies of Abeta and other aggregating peptides and proteins. In this study, hydrogen-deuterium exchange (HX) detected by mass spectrometry (MS) was used to measure Abeta(1-40) aggregate distributions without purification or modification that might alter the aggregate structure or distribution. Different peaks in the mass spectra were assigned to monomer, low molecular weight oligomer, intermediate, and fibril based on HX labeling behavior and complementary assays. After 1 h labeling, the intermediates incorporated approximately ten more deuterons relative to fibrils, indicating a more solvent exposed structure of such intermediates. HX-MS also showed that the intermediate species dissociated much more slowly to monomer than did the very low molecular weight oligomers that were formed at very early times in Abeta aggregation. Atomic force microscopy (AFM) measurements revealed the intermediates were roughly spherical with relatively homogenous diameters of 30-50 nm. Quantitative analysis of the HX mass spectra showed that the amount of intermediate species was correlated with Abeta toxicity patterns reported in a previous study under the same conditions. This study also demonstrates the potential of the HX-MS approach to characterizing complex, multi-component oligomer distributions of aggregating peptides and 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.     

Dynamic motions of free and bound O29 scaffolding protein identified by hydrogen deuterium exchange mass spectrometry

In the double-stranded DNA containing bacteriophages, hundreds of copies of capsid protein subunits polymerize to form icosahedral shells, called procapsids, into which the viral genome is subsequently packaged to form infectious virions. High assembly fidelity requires the assistance of scaffolding protein molecules, which interact with the capsid proteins to insure proper geometrical incorporation of subunits into the growing icosahedral lattices. The interactions between the scaffolding and capsid proteins are transient and are subsequently disrupted during DNA packaging. Removal of scaffolding protein is achieved either by proteolysis or alternatively by some form of conformational switch that allows it to dissociate from the capsid. To identify the switch controlling scaffolding protein association and release, hydrogen deuterium exchange was applied to Bacillus subtilis phage Ø29 scaffolding protein gp7 in both free and procapsid-bound forms. The H/D exchange experiments revealed highly dynamic and cooperative opening motions of scaffolding molecules in the N-terminal helix-loop-helix (H-L-H) region. The motions can be promoted by destabilizing the hydrophobic contact between two helices. At low temperature where high energy motions were damped, or in a mutant in which the helices were tethered through the introduction of a disulfide bond, this region displayed restricted cooperative opening motions as demonstrated by a switch in the exchange kinetics from correlated EX1 exchange to uncorrelated EX2 exchange. The cooperative opening rate was increased in the procapsid-bound form, suggesting this region might interact with the capsid protein. Its dynamic nature might play a role in the assembly and release mechanism.     

Structural features of interferon-gamma aggregation revealed by hydrogen exchange

Using hydrogen-deuterium exchange (HX) and electrospray ionization mass spectrometry, we have investigated the stability and structural changes of recombinant human interferon-gamma (IFN-gamma) during aggregation induced by guanidine hydrochloride (GdnHCl) and potassium thiocyanate. First, HX labeling was initiated after the amorphous aggregates were formed to probe the tertiary structure of the aggregated state. Second, labeling was performed at low protein concentrations to assess stability under aggregation prone conditions. In 1 M GdnHCl, the stability of IFN-gamma was greatly reduced and much less protection from HX in solution was observed. Exchange under these conditions was slower in helix C than in the rest of the protein. Aggregates formed in 1 M GdnHCl showed a HX pattern consistent with a partially unfolded state with an intact helix C. Although aggregates formed in 0.3 M KSCN exhibited a HX pattern similar to those formed in GdnHCl, the solution phase HX pattern in 0.3 M KSCN was surprisingly comparable to that of the native state. Varying the aggregation time before performing HX revealed that KSCN first precipitated native protein and then facilitated partial unfolding of the precipitated protein. These results show that helix C, which forms the hydrophobic core of the IFN-gamma dimer, is highly protected from HX under native conditions, is more stable in GdnHCl than the rest of the protein and remains intact in both GdnHCl- and KSCN-induced aggregates. This suggests that native-state HX patterns may presage regions of the protein susceptible to unfolding during aggregation.     

Mapping the energy landscape of repeat proteins using NMR-detected hydrogen exchange

Repeat proteins contain tandem arrays of a simple structural motif. In contrast to globular proteins, repeat proteins are stabilized only by interactions between residues that are relatively close together in the sequence, with no ”long-range” interactions. Our work focuses on the tetratricopeptide repeat (TPR), a 34 amino acid helix-turn-helix motif found in tandem arrays in many natural proteins. Earlier, we reported the design and characterization of a series of consensus TPR (CTPR) proteins, which are built as arrays of multiple tandem copies of a 34 amino acid consensus sequence. Here, we present the results of extensive hydrogen exchange (HX) studies of the folding-unfolding behavior of two CTPR proteins (CTPR2 and CTPR3). We used HX to detect and characterize partially folded species that are populated at low frequency in the nominally folded state. We show that for both proteins the equilibrium folding-unfolding transition is non-two-state, but sequential, with the outermost helices showing a significantly higher probability than inner helices of being unfolded. We show that the experimentally observed unfolding behavior is consistent with the predictions of a simple Ising model, in which individual helices are treated as ”spin-equivalents”. The results that we present have general implications for our understanding of the thermodynamic properties of repeat proteins.     

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.     

Protein structural dynamics at the gas/water interface examined by hydrogen exchange mass spectrometry

Gas/water interfaces (such as air bubbles or foam) are detrimental to the stability of proteins, often causing aggregation. This represents a potential problem for industrial processes, for example, the production and handling of protein drugs. Proteins possess surfactant-like properties, resulting in a high affinity for gas/water interfaces. The tendency of previously buried nonpolar residues to maximize contact with the gas phase can cause significant structural distortion. Most earlier studies in this area employed spectroscopic tools that could only provide limited information. Here we use hydrogen/deuterium exchange (HDX) mass spectrometry (MS) for probing the conformational dynamics of the model protein myoglobin (Mb) in the presence of N₂ bubbles. HDX/MS relies on the principle that unfolded and/or highly dynamic regions undergo faster deuteration than tightly folded segments. In bubble-free solution Mb displays EX2 behavior, reflecting the occurrence of short-lived excursions to partially unfolded conformers. A dramatically different behavior is seen in the presence of N₂ bubbles; EX2 dynamics still take place, but in addition the protein shows EX1 behavior. The latter results from interconversion of the native state with conformers that are globally unfolded and long-lived. These unfolded species likely correspond to Mb that is adsorbed to the surface of gas bubbles. N₂ sparging also induces aggregation. To explain the observed behavior we propose a simple model, that is, “semi-unfolded” ↔ “native” ↔ “globally unfolded” → “aggregated”. This model quantitatively reproduces the experimentally observed kinetics. To the best of our knowledge, the current study marks the first exploration of surface denaturation phenomena by HDX/MS.     

Relationship between the native-state hydrogen exchange and the folding pathways of barnase

The previous native-state hydrogen exchange experiment with barnase failed to detect any partially unfolded intermediate state which was contrary to the experimental results from kinetic deuterium hydrogen exchange pulse labeling and protein engineering studies. This has been taken to suggest that the native-state hydrogen exchange method cannot be used alone as an analytical tool to study the folding pathways of proteins. Here, we revisited the pulse labeling experiment with barnase and detected no stable folding intermediate. This finding allows a reconciliation of the native-state HX data and the folding pathway of barnase. Along with alternative theoretical interpretations for a curved chevron plot of protein folding, these data suggest that further investigation of the nature of the intermediate of barnase is needed.     

Hydrogen-exchange mass spectrometry reveals activation-induced changes in the conformational mobility of p38alpha MAP kinase

Hydrogen-deuterium exchange measurements represent a powerful approach to investigating changes in conformation and conformational mobility in proteins. Here, we examine p38α MAP kinase (MAPK) by hydrogen-exchange (HX) mass spectrometry to determine whether changes in conformational mobility may be induced by kinase phosphorylation and activation. Factors influencing sequence coverage in the HX mass spectrometry experiment, which show that varying sampling depths, instruments, and peptide search strategies yield the highest coverage of exchangeable amides, are examined. Patterns of regional deuteration in p38α are consistent with tertiary structure and similar to deuteration patterns previously determined for extracellular-signal-regulated kinase (ERK) 2, indicating that MAPKs are conserved with respect to the extent of local amide HX. Activation of p38α alters HX in five regions, which are interpreted by comparing X-ray structures of unphosphorylated p38α and X-ray structures of phosphorylated p38γ. Conformational differences account for altered HX within the activation lip, the P + 1 site, and the active site. In contrast, HX alterations are ascribed to activation-induced effects on conformational mobility, within substrate-docking sites (αF-αG, β7-β8), the C-terminal core (αE), and the N-terminal core region (β4-β5, αL16, αC). Activation also decreases HX in a 3-10 helix at the C-terminal extension of p38α. Although this helix in ERK2 forms a dimerization interface that becomes protected from HX upon activation, analytical ultracentrifugation shows that this does not occur in p38α because both unphosphorylated and diphosphorylated forms are monomeric. Finally, HX patterns in monophosphorylated p38α are similar to those in unphosphorylated kinase, indicating that the major activation lip remodeling events occur only after diphosphorylation. Importantly, patterns of activation-induced HX show differences between p38α and ERK2 despite their similarities in overall deuteration, suggesting that although MAPKs are closely related with respect to primary sequence and tertiary structure, they have distinct mechanisms for dynamic control of enzyme function.     

Assessment of differences in the conformational flexibility of hepatitis B virus core‐antigen and e‐antigen by hydrogen deuterium exchange‐mass spectrometry

Hepatitis B virus core-antigen (capsid protein) and e-antigen (an immune regulator) have almost complete sequence identity, yet the dimeric proteins (termed Cp149_d and Cp(−10)149_d, respectively) adopt quite distinct quaternary structures. Here we use hydrogen deuterium exchange-mass spectrometry (HDX-MS) to study their structural properties. We detect many regions that differ substantially in their HDX dynamics. Significantly, whilst all regions in Cp(−10)149_d exchange by EX2-type kinetics, a number of regions in Cp149_d were shown to exhibit a mixture of EX2- and EX1-type kinetics, hinting at conformational heterogeneity in these regions. Comparison of the HDX of the free Cp149_d with that in assembled capsids (Cp149_c) indicated increased resistance to exchange at the C-terminus where the inter-dimer contacts occur. Furthermore, evidence of mixed exchange kinetics were not observed in Cp149_c, implying a reduction in flexibility upon capsid formation. Cp(−10)149_d undergoes a drastic structural change when the intermolecular disulphide bridge is reduced, adopting a Cp149_d-like structure, as evidenced by the detected HDX dynamics being more consistent with Cp149_d in many, albeit not all, regions. These results demonstrate the highly dynamic nature of these similar proteins. To probe the effect of these structural differences on the resulting antigenicity, we investigated binding of the antibody fragment (Fab E1) that is known to bind a conformational epitope on the four-helix bundle. Whilst Fab E1 binds to Cp149_c and Cp149_d, it does not bind non-reduced and reduced Cp(−10)149_d, despite unhindered access to the epitope. These results imply a remarkable sensitivity of this epitope to its structural context.     

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.     

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?