Herbicide binding and thermal stability of photosystem II isolated from Thermosynechococcus elongatus

Binding of herbicides to photosystem II inhibits the electron transfer from Q(A) to Q(B) due to competition of herbicides with plastoquinone bound at the Q(B) site. We investigated herbicide binding to monomeric and dimeric photosystem II core complexes (PSIIcc) isolated from Thermosynechococcus elongatus by a combination of different methods (isothermal titration and differential scanning calorimetry, CD spectroscopy and measurements of the oxygen evolution) yielding binding constants, enthalpies and stoichiometries for various herbicides as well as information regarding stabilization/destabilization of the complex. Herbicide binding to detergent-solubilized PSIIcc can be described by a model of single independent binding sites present on this important membrane protein. Interestingly, binding stoichiometries herbicide:PSIIcc are lower than 1:1 and vary depending on the herbicide under study. Strong binding herbicides such as terbutryn stabilize PSIIcc in thermal unfolding experiments and endothermically binding herbicides like ioxynil probably cause large structural changes accompanied with the binding process as shown by differential scanning calorimetry experiments of the unfolding reaction of PSIIcc monomer in the presence of ioxynil. In addition we studied the occupancy of the Q(B) sites with plastoquinone (PQ9) by measuring flash induced fluorescence relaxation yielding a possible explanation for the deviations of herbicide binding from a 1:1 herbicide/binding site model.     

Herbicide binding and thermal stability of photosystem II isolated from Thermosynechococcus elongatus

Binding of herbicides to photosystem II inhibits the electron transfer from Q_A to Q_B due to competition of herbicides with plastoquinone bound at the Q_B site. We investigated herbicide binding to monomeric and dimeric photosystem II core complexes (PSIIcc) isolated from Thermosynechococcus elongatus by a combination of different methods (isothermal titration and differential scanning calorimetry, CD spectroscopy and measurements of the oxygen evolution) yielding binding constants, enthalpies and stoichiometries for various herbicides as well as information regarding stabilization/destabilization of the complex. Herbicide binding to detergent-solubilized PSIIcc can be described by a model of single independent binding sites present on this important membrane protein. Interestingly, binding stoichiometries herbicide:PSIIcc are lower than 1:1 and vary depending on the herbicide under study. Strong binding herbicides such as terbutryn stabilize PSIIcc in thermal unfolding experiments and endothermically binding herbicides like ioxynil probably cause large structural changes accompanied with the binding process as shown by differential scanning calorimetry experiments of the unfolding reaction of PSIIcc monomer in the presence of ioxynil. In addition we studied the occupancy of the Q_B sites with plastoquinone (PQ9) by measuring flash induced fluorescence relaxation yielding a possible explanation for the deviations of herbicide binding from a 1:1 herbicide/binding site model.     

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of DeltaCp of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.     

A calorimetric comparison of the interaction of sodium dodecyl sulfate with cytochrome c and erythrocyte glycoproteins.

The interactions of sodium dodecyl sulfate with cytochrome c and erythrocyte glycoproteins have been studied by the method of titration calorimetry. It was found that the initial addition of sodium dodecyl sulfate to cytochrome c caused an endothermic unfolding of the protein, detectable by circular dichroism (CD). This was followed by the exothermic binding of sodium dodecyl sulfate to the protein, without further CD-detectable conformational changes. In contrast, sodium dodecyl sulfate bound directly to the erythrocyte glycoproteins in an exothermic reaction without any accompanying CD-detectable conformation changes. This indicates that the glycoproteins solubilized in aqueous media have exposed hydrophobic regions which can interact directly with this detergent. The enthalpy changes and stoichiometries of binding are reported.     

Stabilization of proteins by ligand binding: application to drug screening and determination of unfolding energetics

The observed stability of a protein is altered when ligands bind, which results in a shift in the melting temperature (T(m)). Binding to the native state in the absence of binding to the denatured state will necessarily lead to an increase in the T(m), while binding to the unfolded state in the absence of native state binding will decrease the T(m) relative to that of the protein in the absence of ligand. These effects are required by the thermodynamics of reversible folding. However, the relationship between binding affinity and the magnitude of the observed temperature shift is not a simple correlation (i.e., a larger shift in T(m) does not necessarily mean tighter binding) and is complicated by interaction with the denatured state. Using exact simulations, the range of behavior for the dependence of the observed T(m) shift on the energetics of ligand binding is investigated here. Specifically, differential scanning calorimetry (DSC) curves are simulated for protein unfolding in the presence of ligands binding to both the native and denatured states. The results have implications for drug screening and the determination of heat capacity changes for protein unfolding.     

Thermal unfolding and conformational stability of the recombinant domain II of glutamate dehydrogenase from the hyperthermophile Thermotoga maritima

Domain II (residues 189-338, M(r) = 16 222) of glutamate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima was used as a model system to study reversible unfolding thermodynamics of this hyperthermostable enzyme. The protein was produced in large quantities in E.COLI: using a T7 expression system. It was shown that the recombinant domain is monomeric in solution and that it comprises secondary structural elements similar to those observed in the crystal structure of the hexameric enzyme.The recombinant domain is thermostable and undergoes reversible and cooperative thermal unfolding in the pH range 5.90-8.00 with melting temperatures between 75.1 and 68.0 degrees C. Thermal unfolding of the protein was studied using differential scanning calorimetry and circular dichroism spectroscopy. Both methods yielded comparable values. The analysis revealed an unfolding enthalpy at 70 degrees C of 70.2 +/- 4.0 kcal/mol and a DeltaC(p) value of 1.4 +/- 0.3 kcal/mol K. Chemical unfolding of the recombinant domain resulted in m values of 3.36 +/- 0.10 kcal/mol M for unfolding in guanidinium chloride and 1.46 +/- 0.04 kcal/mol M in urea. The thermodynamic parameters for thermal and chemical unfolding equilibria indicate that domain II from T.MARITIMA: glutamate dehydrogenase is a thermostable protein with a DeltaG(max) of 3.70 kcal/mol. However, the thermal and chemical stabilities of the domain are lower than those of the hexameric protein, indicating that interdomain interactions must play a significant role in the stabilization of T. MARITIMA: domain II glutamate dehydrogenase.     

Lectins as probes for assessing the accessibility of N‐linked glycans in relation to the conformational changes of fibronectin

Fibronectin, a ≈450‐kDa protein with 4–9% (w/w) glycosylation, is a key component of extracellular matrices and has a high conformational lability regarding its functions. However, the accessibility and the role of glycosylated moieties associated with the conformational changes of fibronectin are poorly understood. Using lectins as probes, we developed an approach comprising dynamic light scattering, turbidimetry measurements, and isothermal titration calorimetry to assess the accessibility of glycosylated moieties of fibronectin undergoing thermal‐induced conformational changes. Among a set of 14 lectins, fibronectin mainly reacted with mannose‐binding lectins, specifically concanavalin A. When temperature was raised from 25 to 50 °C, fibronectin underwent progressive unfolding, but the conformation of concanavalin A was unaffected. Dynamic light scattering, turbidimetry measurements, and isothermal titration calorimetry showed increased concanavalin A binding to fibronectin during progressive thermal‐induced unfolding of the protein core. Such data suggest that mannosylated residues are progressively exposed as fibronectin unfolds. Because oligosaccharide moieties can be differently exposed to cells, and the cell's responses could be modified physiologically or pathologically, modulation of fibronectin sugar chains could be relevant to its biological functions. Thus, lectins might be useful tools to probe the glycosylation accessibility accompanying changes in protein core folding, for which a better understanding would be of value for biological and biomedical research. Copyright © 2015 John Wiley & Sons, Ltd.     

Mechanism of chaperone function in small heat shock proteins. Two-mode binding of the excited states of T4 lysozyme mutants by alphaA-crystallin

To elucidate the mechanism of alphaA-crystallin chaperone function, a detailed thermodynamic analysis of its binding to destabilized, site-directed mutants of T4 lysozyme was carried out. The selected mutants form a ladder of stabilities spanning the 5-10 kcal/mol range of free energy of unfolding. The crystal structures of the majority of the mutants have been previously determined and found to be similar to that of the wild type with no evidence of static local unfolding. Complex formation between alphaA-crystallin and T4 lysozyme was observed directly via the changes in the electron paramagnetic resonance lineshape of a nitroxide introduced at a non-destabilizing, solvent exposed site in T4 lysozyme. AlphaA-crystallin differentially interacts with the mutants, binding the more destabilized ones to a larger extent despite the similar structure of their native states. Our results suggest that the states recognized by alphaA-crystallin are non-native excited states distinct from the unfolded state. Stable complexes are formed when the free energy of binding to alphaA-crystallin is on the order of the free energy associated with the transition from the excited state to the native state. Biphasic binding isotherms reveal two modes of interactions with distinct affinities and stoichiometries. Highly destabilized mutants preferentially bind to the high capacity mode, suggesting conformational preference in the use of each mode. Furthermore, binding can be enhanced by increased temperature and pH, which may be reflecting conformational changes in alphaA-crystallin oligomeric structure.     

Effect of prototypic drugs ibuprofen and warfarin on global chaotropic unfolding of human serum heme-albumin: a fast-field-cycling 1H-NMR relaxometric study

Human serum albumin (HSA) is the most prominent protein in plasma, but it is also found in tissues and secretions throughout the body. The three-domain design of HSA provides a variety of binding sites for many ligands, including heme and drugs. HSA has been used as a model multidomain protein to investigate how interdomain interactions affect the global folding/unfolding process. Here, we report on the reversible chemical denaturation of heme-HSA involving three different conformational states (F, N, and B, occurring at pH 4.0, 7.0, and 9.0, respectively) and on the effect of prototypic drugs ibuprofen and warfarin on thermodynamics of the reversible unfolding process. Chaotropic unfolding of heme-HSA in the F, N, and B conformations is governed by different thermodynamic regimes, with the B form showing an entropic stabilization of the structure that compensates an enthalpic destabilization, and the F form easily unfolding under entropic control. Warfarin and ibuprofen binding stabilizes heme-HSA in both N and B states.     

Cooperative binding of concanavalin A to thymocytes at 4 degrees C and micro-redistribution of concanavalin A receptors.

The mode of binding of 125I-labelled concanavalin A and succinyl-concanavalin A to rat thymocytes at 4 degrees C was investigated. Simultaneously, the free binding sites of the cell-bound lectin molecules were quantified by horseradish peroxidase binding. Concanavalin A showed cooperative binding while succinyl-concanavalin A did not. The number of molecules of concanavalin A bound to the cell surface when it was saturated was twice the number of molecules of succinyl-concanavalin A. We interpret these results as showing that the binding of native concanavalin A to thymocytes at 4 degrees C brings about a cooperative modification of the membrane which leads to appearance of new receptors. Divalent succinyl-concanavalin A has no such effect. Horseradish peroxidase binding to cell-bound lectin was shown to be related to the immobilization of membrane receptors; the more they are immobilized, the more receptor-associated lectin can bind horseradish peroxidase. This allowed us to establish that post-binding events, which we called micro-redistribution, occurred at 4 degrees C when either concanavalin A or succinyl-concanavalin A binds to cells. A cooperative restriction of the micromobility of cell receptors is produced by increasing concentrations of concanavalin A. Succinyl-concanavalin A does not restrict cell receptor mobility at any concentration tested. The results are discussed in terms of cell stimulation and cell agglutination.     

Cooperative binding of inhaled anesthetics and ATP to firefly luciferase

Firefly luciferase is considered a reasonable model of in vivo anesthetic targets despite being destabilized by anesthetics, as reflected by differential scanning calorimetry (DSC). We examined the interaction between two inhaled anesthetics, ATP, luciferase, and temperature, using amide hydrogen exchange, tryptophan fluorescence, and photolabeling in an attempt to examine this apparent discrepancy. In the absence of ATP/Mg2+, halothane and bromoform cause destabilization, as measured by hydrogen exchange, suggesting nonspecific interactions. In the presence of ATP/Mg2+ and at room temperature, the anesthetics produce considerable stabilization with a negative DeltaH, indicating population of a conformer with a specific anesthetic binding site. Stabilizing interactions are lost, however, at unfolding temperatures. We suggest that preferential binding to aggregated forms of luciferase explain the higher temperature destabilization detected with DSC. Our results demonstrate a cooperative binding equilibrium between native ligands and anesthetics, suggesting that similar interactions could underlie actions at biologically relevant targets.     

Synthetic affinity ligands as a novel tool to improve protein stability

Cutinase from Fusarium solani pisi is the model-system for a new approach to assess and enhance protein stability based on the use of synthetic triazine-scaffolded affinity ligands as a novel protein-stabilizing tool. The active site of cutinase is excluded from the main surface regions postulated to be involved in early protein's thermal unfolding events. Hence, these regions are suitable targets for binding complementary affinity ligands with a potential stabilizing effect. A random solid-phase combinatorial library of triazine-bisubstituted molecules was screened for binding cutinase by a rapid fluorescence-based method and affinity chromatography. The best binding substituents were combined with those previously selected by screening a rationally designed library. A second-generation solid-phase biased library was designed and synthesized, following a semi-rational methodology. A dual screening of this library enabled the selection of ligands binding cutinase with higher affinity while retaining its functionality. These compounds were utilized for thermostability assessment with adsorbed cutinase at 60 degrees C and pH 8.0. When bound to different types of ligands, the enzyme showed markedly distinct activity retention profiles, with some synthetic affinity ligands displaying a stabilizing effect on cutinase and others a clearly destabilizing effect, when compared with the free enzyme.     

A differential scanning calorimetric study of the influence of copper and dodecyl trimethyl ammonium bromide on the stability of bovine alpha-lactalbumin

Bovine alpha-lactalbumin (alpha-LA) has been studied by differential scanning calorimetry (DSC), fluorescence spectroscopy and viscometry with various concentrations of Cu2+ and DTAB to elucidate the effect of these ligands on its thermal properties. The DSC profile of dialyzed form of alpha-lactalbumin (m-alpha-LA) contrary to the undialyzed form (holo-form, h-alpha-LA) shows two temperature induced heat absorption peaks. The m-alpha-LA is not a new form of alpha-LA. It contains mixture of the apo (a-alpha-LA) and holo (h-alpha-LA) forms of alpha-LA at low and high temperatures, respectively. Therefore, these two states of alpha-LA (apo and holo) are equilibrating with together after dialyze experiment. The Cu2+ as a metal ion and DTAB as a non metal ion alter the two heat-absorption peaks, in such a manner that, the addition of Cu2+ to the m-alpha-LA increases partial molar heat capacity and enthalpy change values of the h-alpha-LA form at high temperature because the molecular population of the a-alpha-LA form changes into the h-like-alpha-LA. On the contrary, the interaction between the DTAB and the m-alpha-LA increases these thermodynamic values for the a-alpha-LA at low temperature. However, DTAB bound to m-alpha-LA prevents from Ca2+ binding to protein, because there are positive charges repulsion between them. The high temperature peak occurs at the same temperature as the unfolding of the h-alpha-LA, while the low temperature peak lies within the temperature range associated with the unfolding of the a-alpha-LA. The R(s) values of m-alpha-LA, h-alpha-LA and a-alpha-LA forms confirmed the folding and unfolding of the m-alpha-LA during the addition of Cu2+ and DTAB at different concentration, respectively.     

Comparative studies of unfolding and binding of ligands to human serum albumin in the presence of fatty acid: spectroscopic approach

This study was undertaken to investigate the influence of fatty acid binding on the unfolding of HSA and how the fatty acid molecules can influence and/or compete with other ligand molecules bound to the protein. The equilibrium unfolding of fatted and fatty acid free HSA was measured by overlapping of unfolding transition curves monitored by different probes for secondary and tertiary structure and determining changes in free energy of unfolding. Proteins stability was studied by fluorescence spectroscopy, whereas conformational changes were detected by circular dichroism techniques. We have suggested a "molten globule" like intermediate state of HSA at a fairly high concentration of GnHCl (3.2 for fatty acid free and 3.6 for fatted). The free energy of stabilization (DeltaG(D)(H2O)) in the presence of fatty acid was found to be 900 cal mol(-1). We also analyze the effects of fatty acid on binding of ligands using spectroscopic technique and reported the equilibrium constants and free energies obtained from the binding and unfolding experiments.     

Protein stabilization by osmolytes from hyperthermophiles: effect of mannosylglycerate on the thermal unfolding of recombinant nuclease a from Staphylococcus aureus studied by picosecond time-resolved fluorescence and calorimetry

2-O-alpha-Mannosylglycerate, a negatively charged osmolyte widely distributed among (hyper)thermophilic microorganisms, is known to provide notable protection to proteins against thermal denaturation. To study the mechanism responsible for protein stabilization, pico-second time-resolved fluorescence spectroscopy was used to characterize the thermal unfolding of a model protein, Staphylococcus aureus recombinant nuclease A (SNase), in the presence or absence of mannosylglycerate. The fluorescence decay times are signatures of the protein state, and the pre-exponential coefficients are used to evaluate the molar fractions of the folded and unfolded states. Hence, direct determination of equilibrium constants of unfolding from molar fractions was carried out. Van't Hoff plots of the equilibrium constants provided reliable thermodynamic data for SNase unfolding. Differential scanning calorimetry was used to validate this thermodynamic analysis. The presence of 0.5 m potassium mannosylglycerate caused an increase of 7 degrees C in the SNase melting temperature and a 2-fold increase in the unfolding heat capacity. Despite the considerable degree of stabilization rendered by this solute, the nature and population of protein states along unfolding were not altered in the presence of mannosylglycerate, denoting that the unfolding pathway of SNase was unaffected. The stabilization of SNase by mannosylglycerate arises from decreased unfolding entropy up to 65 degrees C and from an enthalpy increase above this temperature. In molecular terms, stabilization is interpreted as resulting from destabilization of the denatured state caused by preferential exclusion of the solute from the protein hydration shell upon unfolding, and stabilization of the native state by specific interactions. The physiological significance of charged solutes in hyperthermophiles is discussed.     

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two‐state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non‐cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier‐less “downhill” folding, as well as for continuous “uphill” unfolding transitions, indicate that gradual non‐cooperative processes may be ubiquitous features on the free energy landscape of protein folding.     

Thermodynamics of intersubunit interactions in cholera toxin upon binding to the oligosaccharide portion of its cell surface receptor, ganglioside GM1

The binding and the energetics of the interaction of cholera toxin with the oligosaccharide portion of ganglioside GM1 (oligo-GM1), the toxin cell surface receptor, have been studied by high-sensitivity isothermal titration calorimetry and differential scanning calorimetry. Previously, we have shown that the association of cholera toxin to ganglioside GM1 enhances the cooperative interactions between subunits in the B-subunit pentamer [Goins, B., & Freire, E. (1988) Biochemistry 27, 2046-2052]. New experiments presented in this paper reveal that the oligosaccharide portion of the receptor is by itself able to enhance the intersubunit cooperative interactions within the B pentamer. This effect is seen in the protein unfolding transition as a shift from independent unfolding of the B promoters toward a cooperative unfolding. To identify the origin of this effect, the binding of cholera toxin to oligo-GM1 has been measured calorimetrically under isothermal conditions. The binding curve at 37 degrees C is sigmoidal, indicating cooperative binding. The binding data can be described in terms of a nearest-neighbor cooperative interaction binding model. In terms of this model, the association of a oligo-GM1 molecule to a B protomer affects the association to adjacent B promoters within the pentameric ring. The measured intrinsic binding enthalpy per protomer is -22 kcal/mol and the cooperative interaction enthalpy -11 kcal/mol. The intrinsic binding constant determined calorimetrically is 1.05 x 10(6) M-1 at 37 degrees C and the cooperative Gibbs free energy equal to -850 cal/mol.(ABSTRACT TRUNCATED AT 250 WORDS)     

Why has porcine VEG protein unusually high stability and suppressed binding ability?

Von Ebner gland protein (VEGP) and odorant-binding protein (OBP) were purified from porcine lingual epithelium and nasal mucosa, respectively. Both VEGP and OBP preparations were homogeneous as indicated by SDS-PAGE, isoelectric focusing, gel-filtration and electrospray mass spectrometry. However, high-sensitivity differential scanning calorimetry (HS-DSC) yielded multiphasic denaturation thermograms for both proteins indicating their conformational heterogeneity. The unfolding transition of VEGP is observed at extremely high temperatures (about 110 degrees C), which is unexpected for a protein with significant structural homology to OBP and other lipocalins. Isothermal titration calorimetry (ITC) did not detect the binding of either aspartame or denatonium saccharide to VEGP nor did it detect binding of 2-isobutyl-3-methoxypyrazine (IBMP) to OBP. Extraction of OBP with mixed organic solvents eliminated the conformational heterogeneity and the protein showed a reversible two-state transition in HS-DSC thereafter. ITC also showed that the extracted OBP was able to bind IBMP. These results imply that tightly bound endogenous ligands increase the thermal stability of OBP and block the binding of other ligands. In contrast to OBP, the extraction of VEGP with organic solvents failed to promote binding or to establish thermal homogeneity, most likely because of the irreversible denaturation of VEGP. Thus, the elucidation of the functional behaviour of VEGP is closely related to the exhaustive purging of its endogenous ligands which otherwise very efficiently mask ligand binding sites of this protein.     

Methionine-121 coordination determines metal specificity in unfolded<Emphasis Type="Italic"> Pseudomonas aeruginosa</Emphasis> azurin

Pseudomonas aeruginosa azurin binds copper so tightly that it remains bound even upon polypeptide unfolding. Copper can be substituted with zinc without change in protein structure, and also in this complex the metal remains bound upon protein unfolding. Previous work has shown that native-state copper ligands Cys112 and His117 are two of at least three metal ligands in the unfolded state. In this study we use isothermal titration calorimetry and spectroscopic methods to test if the native-state ligand Met121 remains a metal ligand upon unfolding. From studies on a point-mutated version of azurin (Met121Ala) and a set of model peptides spanning the copper-binding C-terminal part (including Cys112, His117 and Met121), we conclude that Met121 is a metal ligand in unfolded copper-azurin but not in the case of unfolded zinc-azurin. Combination of unfolding and metal-titration data allow for determination of copper (Cu(II) and Cu(I)) and zinc affinities for folded and unfolded azurin polypeptides, respectively.     

Trans-substitution of the proximal hydrogen bond in myoglobin: II. Energetics, functional consequences, and implications for hemoglobin allostery

The trans-substituted histidine to glycine mutant of sperm whale myoglobin (H93G Mb) is used to study energetics of proximal hydrogen bonding, proximal ligand-heme interactions, and coupling to distal ligand binding. Comparison of mono- and dimethylimidazole structural isomers shows that the hydrogen bond between the proximal ligand and the neighboring Ser92 hydroxyl (position F7) is stabilizing. The range of hydrogen bond stabilities measured here for different distal ligand complexes ranges from -0.7 kcal/mol (monomethylimidazole isomers to MbCO) to -4.1 kcal/mol (dimethylimidazole isomers to MbCN). This range of hydrogen bond stabilities, which is similar to that seen in protein mutagenesis unfolding studies, demonstrates the high sensitivity of the hydrogen bond to modest structural perturbations. The degree to which the 2-methyl group destabilizes proximal ligand binding is found to depend inversely on the total electronic spin. For monomethylimidazole proximal ligands, distal ligand binding weakens the proximal hydrogen bond compared to deoxyMb. Surprisingly, this trend is largely reversed for the dimethylimidazole proximal ligands. These results demonstrate strong coupling between the proximal protein matrix and distal ligand binding. These results provide an explanation for the strong avoidance of hydrogen bonding residues at position F7 in hemoglobin sequences.