A comparative NMR study of the polypeptide backbone dynamics of hemoglobin in the deoxy and carbonmonoxy forms

Comparison of the Backbone Dynamics of Dehaloperoxidase-Hemoglobin Isoenzymes

Dehaloperoxidase (DHP) is a multifunctional hemeprotein with a functional switch generally regulated by the chemical class of the substrate. Its two isoforms, DHP-A and DHP-B, differ by only five amino acids and have an almost identical protein fold. However, the catalytic efficiency of DHP-B for oxidation by a peroxidase mechanism ranges from 2- to 6-fold greater than that of DHP-A depending on the conditions. X-ray crystallography has shown that many substrates and ligands have nearly identical binding in the two isoenzymes, suggesting that the difference in catalytic efficiency could be due to differences in the conformational dynamics. We compared the backbone dynamics of the DHP isoenzymes at pH 7 through heteronuclear relaxation dynamics at 11.75, 16.45, and 19.97 T in combination with four 300 ns MD simulations. While the overall dynamics of the isoenzymes are similar, there are specific local differences in functional regions of each protein. In DHP-A, Phe35 undergoes a slow chemical exchange between two conformational states likely coupled to a swinging motion of Tyr34. Moreover, Asn37 undergoes fast chemical exchange in DHP-A. Given that Phe35 and Asn37 are adjacent to Tyr34 and Tyr38, it is possible that their dynamics modulate the formation and migration of the active tyrosyl radicals in DHP-A at pH 7. Another significant difference is that both distal and proximal histidines have a 15-18% smaller S2 value in DHP-B, thus their greater flexibility could account for the higher catalytic activity. The distal histidine grants substrate access to the distal pocket. The greater flexibility of the proximal histidine could also accelerate H2O2 activation at the heme Fe by increased coupling of an amino acid charge relay to stabilize the ferryl Fe(IV) oxidation state in a Poulos-Kraut "push-pull"-type peroxidase mechanism.

Conformational Dynamics of Phytoglobin BvPgb1.2 from Beta vulgaris ssp. vulgaris

Plant hemoglobins, often referred to as phytoglobins, play important roles in abiotic stress tolerance. Several essential small physiological metabolites can be bound to these heme proteins. In addition, phytoglobins can catalyze a range of different oxidative reactions in vivo. These proteins are often oligomeric, but the degree and relevance of subunit interactions are largely unknown. In this study, we delineate which residues are involved in dimer formation of a sugar beet phytoglobin type 1.2 (BvPgb1.2) using NMR relaxation experiments. E. coli cells harboring a phytoglobin expression vector were cultivated in isotope-labeled (²H, ¹³C and ¹⁵N) M9 medium. The triple-labeled protein was purified to homogeneity using two chromatographic steps. Two forms of BvPgb1.2 were examined, the oxy-form and the more stable cyanide-form. Using three-dimensional triple-resonance NMR experiments, sequence-specific assignments for CN-bound BvPgb1.2 were achieved for 137 backbone amide cross-peaks in the ¹H-¹⁵N TROSY spectrum, which amounts to 83% of the total number of 165 expected cross-peaks. A large proportion of the non-assigned residues are located in α-helixes G and H, which are proposed to be involved in protein dimerization. Such knowledge around dimer formation will be instrumental for developing a better understanding of phytoglobins' roles in planta.

Modulation of hemoglobin dynamics by an allosteric effector

Hemoglobin (Hb) is an extensively studied paradigm of proteins that alter their function in response to allosteric effectors. Models of its action have been used as prototypes for structure‐function relationships in many proteins, and models for the molecular basis of its function have been deeply studied and extensively argued. Recent reports suggest that dynamics may play an important role in its function. Relatively little is known about the slow, correlated motions of hemoglobin subunits in various structural states because experimental and computational strategies for their characterization are challenging. Allosteric effectors such as inositol hexaphosphate (IHP) bind to both deoxy‐Hb and HbCO, albeit at different sites, leading to a lowered oxygen affinity. The manner in which these effectors impact oxygen binding is unclear and may involve changes in structure, dynamics or both. Here we use neutron spin echo measurements accompanied by wide‐angle X‐ray scattering to show that binding of IHP to HbCO results in an increase in the rate of coordinated motions of Hb subunits relative to one another with little if any change in large scale structure. This increase of large‐scale dynamics seems to be coupled with a decrease in the average magnitude of higher frequency modes of individual residues. These observations indicate that enhanced dynamic motions contribute to the functional changes induced by IHP and suggest that they may be responsible for the lowered oxygen affinity triggered by these effectors.

Solution NMR Spectroscopy for the Study of Enzyme Allostery

Allostery is a ubiquitous biological regulatory process in which distant binding sites within a protein or enzyme are functionally and thermodynamically coupled. Allosteric interactions play essential roles in many enzymological mechanisms, often facilitating formation of enzyme-substrate complexes and/or product release. Thus, elucidating the forces that drive allostery is critical to understanding the complex transformations of biomolecules. Currently, a number of models exist to describe allosteric behavior, taking into account energetics as well as conformational rearrangements and fluctuations. In the following Review, we discuss the use of solution NMR techniques designed to probe allosteric mechanisms in enzymes. NMR spectroscopy is unequaled in its ability to detect structural and dynamical changes in biomolecules, and the case studies presented herein demonstrate the range of insights to be gained from this valuable method. We also provide a detailed technical discussion of several specialized NMR experiments that are ideally suited for the study of enzymatic allostery.

Molecular interactions of hemoglobin with resveratrol: potential protective antioxidant role and metabolic adaptations of the erythrocyte

This article reports the role of resveratrol in the erythrocyte as a result of its interaction with hemoglobin and describes the effect of this interaction on the metabolism, the redox state, and the release of ATP. The drug crosses the erythrocyte membrane and binds to hemoglobin, altering its modulation and the release of ATP. Our study correlates the variation of the phosphorylation balance induced by resveratrol with the change in the intracellular concentration of ATP and with the decrease in ATP release from red blood cell and the consequent paracrine alteration on the vascular epithelium. Molecular docking calculations indicate larger specificity of binding for oxy-hemoglobin that correlates well with the stabilization of the R-quaternary structure and with the functional modulation of resveratrol on the protein. Finally, we locate a putative binding site at the central cavity of hemoglobin and characterize its key interacting residues with the drug. Computational results support the assumption that resveratrol may act as a protector agent against oxidative protein damage by interacting with hemoglobin.

Insight into the allosteric mechanism of Scapharca dimeric hemoglobin

Allosteric regulation is an essential function of many proteins that control a variety of different processes such as catalysis, signal transduction, and gene regulation. Structural rearrangements have historically been considered the main means of communication between different parts of a protein. Recent studies have highlighted the importance, however, of changes in protein flexibility as an effective way to mediate allosteric communication across a protein. Scapharca dimeric hemoglobin (HbI) is the simplest possible allosteric system, with cooperative ligand binding between two identical subunits. Thermodynamic equilibrium studies of the binding of oxygen to HbI have shown that cooperativity is an entropically driven effect. The change in entropy of the system observed upon ligand binding may arise from changes in the protein, the ligand, or the water of the system. The goal of this study is to determine the contribution of the change in entropy of the protein backbone to HbI cooperative binding. Molecular dynamics simulations and nuclear magnetic resonance relaxation techniques have revealed that the fast internal motions of HbI contribute to the cooperative binding to carbon monoxide in two ways: (1) by contributing favorably to the free energy of the system and (2) by participating in the cooperative mechanism at the HbI subunit interface. The internal dynamics of the weakly cooperative HbI mutant, F97Y, were also investigated with the same methods. The changes in backbone NH dynamics observed for F97Y HbI upon ligand binding are not as large as for the wild type, in agreement with the reduced cooperativity observed for this mutant. The results of this study indicate that interface flexibility and backbone conformational entropy of HbI participate in and are important for the cooperative mechanism of carbon monoxide binding.

Differential control of heme reactivity in alpha and beta subunits of hemoglobin: a combined Raman spectroscopic and computational study

The use of hybrid hemoglobin (Hb), with mesoheme substituted for protoheme, allows separate monitoring of the α or β hemes along the allosteric pathway. Using resonance Raman (rR) spectroscopy in silica gel, which greatly slows protein motions, we have observed that the Fe-histidine stretching frequency, νFeHis, which is a monitor of heme reactivity, evolves between frequencies characteristic of the R and T states, for both α or β chains, prior to the quaternary R-T and T-R shifts. Computation of νFeHis, using QM/MM and the conformational search program PELE, produced remarkable agreement with experiment. Analysis of the PELE structures showed that the νFeHis shifts resulted from heme distortion and, in the α chain, Fe-His bond tilting. These results support the tertiary two-state model of ligand binding (Henry et al., Biophys. Chem. 2002, 98, 149). Experimentally, the νFeHis evolution is faster for β than for α chains, and pump-probe rR spectroscopy in solution reveals an inflection in the νFeHis time course at 3 μs for β but not for α hemes, an interval previously shown to be the first step in the R-T transition. In the α chain νFeHis dropped sharply at 20 μs, the final step in the R-T transition. The time courses are fully consistent with recent computational mapping of the R-T transition via conjugate peak refinement by Karplus and co-workers (Fischer et al., Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 5608). The effector molecule IHP was found to lower νFeHis selectively for α chains within the R state, and a binding site in the α1α2 cleft is suggested.

Solution structure and dynamics of human hemoglobin in the carbonmonoxy form

The solution structure of human adult carbonmonoxy hemoglobin (HbCO A) was refined using stereospecifically assigned methyl groups and residual dipolar couplings based on our previous nuclear magnetic resonance structure. The tertiary structures of individual chains were found to be very similar to the X-ray structures, while the quaternary structures in solution at low salt concentrations resembled the X-ray R structure more than the R2 structure. On the basis of chemical shift perturbation by inositol hexaphosphate (IHP) titration and docking, we identified five possible IHP binding sites in HbCO A. Amide-water proton exchange experiments demonstrated that αThr38 located in the α1β2 interface and several loop regions in both α- and β-chains were dynamic on the subsecond time scale. Side chain methyl dynamics revealed that methyl groups in the α1β2 interface were dynamic, but those in the α1β1 interface were quite rigid on the nanosecond to picosecond and millisecond to microsecond time scales. All the data strongly suggest a dynamic α1β2 interface that allows conformational changes among different forms (like T, R, and R2) easily in solution. Binding of IHP to HbCO A induced small structural and dynamic changes in the α1β2 interface and the regions around the hemes but did not increase the conformational entropy of HbCO A. The binding also caused conformational changes on the millisecond time scale, very likely arising from the relative motion of the α1β1 dimer with respect to the α2β2 dimer. Heterotropic effectors like IHP may change the oxygen affinity of Hb through modulating the relative motion of the two dimers and then further altering the structure of heme binding regions.

Comparative analysis of oxy-hemoglobin and aquomet-hemoglobin by hydrogen/deuterium exchange mass spectrometry

The function of hemoglobin (Hb) as oxygen transporter is mediated by reversible O2 binding to Fe(2+) heme in each of the α and β subunits. X-ray crystallography revealed different subunit arrangements in oxy-Hb and deoxy-Hb. The deoxy state is stabilized by additional contacts, causing a rigidification that results in strong protection against hydrogen/deuterium exchange (HDX). Aquomet-Hb is a dysfunctional degradation product with four water-bound Fe(3+) centers. Heme release from aquomet-Hb is relatively facile, triggering oxidative damage of membrane lipids. Aquomet-Hb crystallizes in virtually the same conformation as oxy-Hb. Hence, it is commonly implied that the solution-phase properties of aquomet-Hb should resemble those of the oxy state. This work compares the structural dynamics of oxy-Hb and aquomet-Hb by HDX mass spectrometry (MS). It is found that the aquomet state exhibits a solution-phase structure that is significantly more dynamic, as manifested by elevated HDX levels. These enhanced dynamics affect the aquomet α and β subunits in a different fashion. The latter undergoes global destabilization, whereas the former shows elevated HDX levels only in the heme binding region. It is proposed that these enhanced dynamics play a role in facilitating heme release from aquomet-Hb. Our findings should be of particular interest to the MS community because oxy-Hb and aquomet-Hb serve as widely used test analytes for probing the relationship between biomolecular structure in solution and in the gas phase. We are not aware of any prior comparative HDX/MS experiments on oxy-Hb and aquomet-Hb.

Coarse-grained and all-atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen-binding protein composed of four subunits (α1, α2, β1, and β2), is a well-known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse-grained (CG) modeling, all-atom simulation, and structural data from X-ray crystallography and wide-angle X-ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot-spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all-atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all-atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling.

Influence of heme post-translational modification and distal ligation on the backbone dynamics of a monomeric hemoglobin

The cyanobacterium Synechococcus sp. PCC 7002 uses a hemoglobin of the truncated lineage (GlbN) in the detoxification of reactive species generated in the assimilation of nitrate. In view of a sensing or enzymatic role, several states of GlbN are of interest with respect to its structure-activity relationship. Nuclear magnetic resonance spectroscopy was applied to compare the structure and backbone dynamics of six GlbN forms differing in their oxidation state [Fe(II) or Fe(III)], distal ligand to the iron (histidine, carbon monoxide, or cyanide), or heme post-translational modification (b heme or covalently attached heme). Structural properties were assessed with pseudocontact shift calculations. (15)N relaxation data were analyzed by reduced spectral density mapping (picosecond to nanosecond motions) and by inspection of elevated R(2) values (microsecond to millisecond motions). On the picosecond to nanosecond time scale, GlbN exhibited little flexibility and was unresponsive to the differences among the various forms. Regions of slightly higher mobility were the CE turn, the EF loop, and the H-H' kink. In contrast, fluctuations on the microsecond to millisecond time scale depended on the form. Cyanide binding to the ferric state did not enhance motions, whereas reduction to the ferrous bis-histidine state resulted in elevated R(2) values for several amides. This response was attributed, at least in part, to a weakening of the distal histidine coordination. Carbon monoxide binding quenched some of these fluctuations. The results emphasized the role of the distal ligand in dictating backbone flexibility and illustrated the multiple ways in which motions are controlled by the hemoglobin fold.

Heme reactivity is uncoupled from quaternary structure in gel-encapsulated hemoglobin: a resonance Raman spectroscopic study

Encapsulation of hemoglobin (Hb) in silica gel preserves structure and function but greatly slows protein motion, thereby providing access to intermediates along the allosteric pathway that are inaccessible in solution. Resonance Raman (RR) spectroscopy with visible and ultraviolet laser excitation provides probes of heme reactivity and of key tertiary and quaternary contacts. These probes were monitored in gels after deoxygenation of oxyHb and after CO binding to deoxyHb, which initiate conformational change in the R-T and T-R directions, respectively. The spectra establish that quaternary structure change in the gel takes a week or more but that the evolution of heme reactivity, as monitored by the Fe-histidine stretching vibration, ν(FeHis), is completed within two days, and is therefore uncoupled from the quaternary structure. Within each quaternary structure, the evolving ν(FeHis) frequencies span the full range of values between those previously associated with the high- and low-affinity end states, R and T. This result supports the tertiary two-state (TTS) model, in which the Hb subunits can adopt high- and low-affinity tertiary structures, r and t, within each quaternary state. The spectra also reveal different tertiary pathways, involving the breaking and reformation of E and F interhelical contacts in the R-T direction but not the T-R direction. In the latter, tertiary motions are restricted by the T quaternary contacts.

T-quaternary structure of oxy human adult hemoglobin in the presence of two allosteric effectors, L35 and IHP

The cooperative O(2)-binding of hemoglobin (Hb) have been assumed to correlate to change in the quaternary structures of Hb: T(deoxy)- and R(oxy)-quaternary structures, having low and high O(2)-affinities, respectively. Heterotropic allosteric effectors have been shown to interact not only with deoxy- but also oxy-Hbs causing significant reduction in their O(2)-affinities and the modulation of cooperativity. In the presence of two potent effectors, L35 and inositol hexaphosphate (IHP) at pH 6.6, Hb exhibits extremely low O(2)-affinities (K(T)=0.0085mmHg(-1) and K(R)=0.011mmHg(-1)) and thus a very low cooperativity (K(R)/K(T)=1.3 and L(0)=2.4). (1)H-NMR spectra of human adult Hb with these two effectors were examined in order to determine the quaternary state of Hb in solution and to clarify the correlation between the O(2)-affinities and the structural change of Hb caused by the heterotropic effectors. At pH 6.9, (1)H-NMR spectrum of deoxy-Hb in the presence of L35 and IHP showed a marker of the T-quaternary structure (the T-marker) at 14ppm, originated from inter- dimeric α(1)β(2)- (or α(2)β(1)-) hydrogen-bonds, and hyperfine-shifted (hfs) signals around 15-25ppm, caused by high-spin heme-Fe(II)s. Upon addition of O(2), the hfs signals disappeared, reflecting that the heme-Fe(II)s are ligated with O(2), but the T-marker signals still remained, although slightly shifted and broadened, under the partial pressure of O(2) (P(O2)) of 760mmHg. These NMR results accompanying with visible absorption spectroscopy and visible resonance Raman spectroscopy reveal that oxy-Hb in the presence of L35 and IHP below pH 7 takes the ligated T-quaternary structure under the P(O2) of 760mmHg. The L35-concentration dependence of the T-marker in the presence of IHP indicates that there are more than one kind of L35-binding sites in the ligated T-quaternary structure. The stronger binding sites are probably intra-dimeric binding sites between α(1)G- and β(1)G-helices, and the other weaker binding site causes the R→T transition without release of O(2). The fluctuation of the tertiary structure of Hb seems to be caused by both the structural perturbation of α(1)β(1) (or α(2)β(2)) intra-dimeric interface, where the stronger L35-binding sites exist, and by the IHP-binding to the α(1)α(2)- (or β(1)β(2)-) cavity. The tertiary structural fluctuation induced by the allosteric effectors may contribute to the significant reduction of the O(2)-affinity of oxy-Hb, which little depends on the quaternary structures. Therefore, the widely held assumptions of the structure-function correlation of Hb - [the deoxy-state]=[the T-quaternary structure]=[the low O(2)-affinity state] and [the oxy-state]=[the R-quaternary structure]=[the high O(2)-affinity state] and the O(2)-affiny of Hb being regulated by the T/R-quaternary structural transition - are no longer sustainable. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

WAXS studies of the structural diversity of hemoglobin in solution

Specific ligation states of hemoglobin are, when crystallized, capable of taking on multiple quaternary structures. The relationship between these structures, captured in crystal lattices, and hemoglobin structure in solution remains uncertain. Wide-angle X-ray solution scattering (WAXS) is a sensitive probe of protein structure in solution that can distinguish among similar structures and has the potential to contribute to these issues. We used WAXS to assess the relationships among the structures of human and bovine hemoglobins in different liganded forms in solution. WAXS data readily distinguished among the various forms of hemoglobins. WAXS patterns confirm some of the relationships among hemoglobin structures that have been defined through crystallography and NMR and extend others. For instance, methemoglobin A in solution is, as expected, nearly indistinguishable from HbCO A. Interestingly, for bovine hemoglobin, the differences between deoxy-Hb, methemoglobin and HbCO are smaller than the corresponding differences in human hemoglobin. WAXS data were also used to assess the spatial extent of structural fluctuations of various hemoglobins in solution. Dynamics has been implicated in allosteric control of hemoglobin, and increased dynamics has been associated with lowered oxygen affinity. Consistent with that notion, WAXS patterns indicate that deoxy-Hb A exhibits substantially larger structural fluctuations than HbCO A. Comparisons between the observed WAXS patterns and those predicted on the basis of atomic coordinate sets suggest that the structures of Hb in different liganded forms exhibit clear differences from known crystal structures.

Backbone dynamics of deoxy and carbonmonoxy hemoglobin by NMR/SRLS

The slowly relaxing local structure (SRLS) approach, developed for NMR spin relaxation analysis in proteins, is applied herein to amide ¹⁵N relaxation in deoxy and carbonmonoxy hemoglobin. Experimental data including ¹⁵N T₁, T₂ and ¹⁵N-{¹H} NOE, acquired at 11.7 and 14.1 T, and 29 and 34 °C, are analyzed. The restricted local motion of the N-H bond is described in terms of the principal value (S(0)(2)) and orientation (β(D)) of an axial local ordering tensor, S, and the principal values (R(||)(L) and R(⊥)(L)) and orientation (β(O)) of an axial local diffusion tensor, R(L). The parameters c₀² (the potential coefficient in terms of which S(0)(2) is defined), R(||)(L), β(D), and β(O) are determined by data fitting; R(⊥)(L) is set equal to the global motional rate, R(C), found previously to be (5.2-5.8) × 10⁶ 1/s in the temperature range investigated. The principal axis of S is (nearly) parallel to the C(i-1)(α)-C(i)(α) axis; when the two axes are parallel, β(D) = -101.3° (in the frame used). The principal axis of R(L) is (nearly) parallel to the N-H bond; when the two axes are parallel, β(O) = -101.3°. For "rigid" N-H bonds located in secondary structure elements the best-fit parameters are S(0)(2) = 0.88-0.95 (corresponding to local potentials of 8.6-19.9 k(B)T), R(||)(L) = 10⁹-10¹⁰ 1/s, β(D) = -101.3° ± 2.0°, and β(O) = -101.3° ± 4°. For flexible N-H bonds located in loops the best-fit values are S(0)(2) = 0.75-0.80 (corresponding to local potentials of 4.5-5.5 k(B)T), R(||)(L) = (1.0-6.3) × 10⁸ 1/s, β(D) = -101.3° ± 4.0°, and β(O) = -101.3° ± 10°. These results are important in view of their physical clarity, inherent potential for further interpretation, consistency, and new qualitative insights provided (vide infra).

An investigation of the distal histidyl hydrogen bonds in oxyhemoglobin: effects of temperature, pH, and inositol hexaphosphate

On the basis of X-ray crystal structures and electron paramagnetic resonance (EPR) measurements, it has been inferred that the O(2) binding to hemoglobin is stabilized by the hydrogen bonds between the oxygen ligands and the distal histidines. Our previous study by multinuclear nuclear magnetic resonance (NMR) spectroscopy has provided the first direct evidence of such H-bonds in human normal adult oxyhemoglobin (HbO(2) A) in solution. Here, the NMR spectra of uniformly (15)N-labeled recombinant human Hb A (rHb A) and five mutant rHbs in the oxy form have been studied under various experimental conditions of pH and temperature and also in the presence of an organic phosphate, inositol hexaphosphate (IHP). We have found significant effects of pH and temperature on the strength of the H-bond markers, i.e., the cross-peaks for the side chains of the two distal histidyl residues, α58His and β63His, which form H-bonds with the O(2) ligands. At lower pH and/or higher temperature, the side chains of the distal histidines appear to be more mobile, and the exchange with water molecules in the distal heme pockets is faster. These changes in the stability of the H-bonds with pH and temperature are consistent with the changes in the O(2) affinity of Hb as a function of pH and temperature and are clearly illustrated by our NMR experiments. Our NMR results have also confirmed that this H-bond in the β-chain is weaker than that in the α-chain and is more sensitive to changes in pH and temperature. IHP has only a minor effect on these H-bond markers compared to the effects of pH and temperature. These H-bonds are sensitive to mutations in the distal heme pockets but not affected directly by the mutations in the quaternary interfaces, i.e., α(1)β(1) and/or α(1)β(2) subunit interface. These findings provide new insights regarding the roles of temperature, hydrogen ion, and organic phosphate in modulating the structure and function of hemoglobin in solution.

Effector-induced structural fluctuation regulates the ligand affinity of an allosteric protein: binding of inositol hexaphosphate has distinct dynamic consequences for the T and R states of hemoglobin

The present study reports distinct dynamic consequences for the T- and R-states of human normal adult hemoglobin (Hb A) due to the binding of a heterotropic allosteric effector, inositol hexaphosphate (IHP). A nuclear magnetic resonance (NMR) technique based on modified transverse relaxation optimized spectroscopy (TROSY) has been used to investigate the effect of conformational exchange of Hb A in both deoxy and CO forms, in the absence and presence of IHP, at 14.1 and 21.1 T, and at 37 degrees C. Our results show that the majority of the polypeptide backbone amino acid residues of deoxy- and carbonmonoxy-forms of Hb A in the absence of IHP is not mobile on the micros-ms time scale, with the exception of several amino acid residues, that is, beta109Val and beta132Lys in deoxy-Hb A, and alpha40Lys in HbCO A. The mobility of alpha40Lys in HbCO A can be explained by the crystallographic data showing that the H-bond between alpha40Lys and beta146His in deoxy-Hb A is absent in HbCO A. However, the conformational exchange of beta109Val, which is located in the intradimer (alpha 1beta 1 or alpha 2beta 2) interface, is not consistent with the crystallographic observations that show rigid packing at this site. IHP binding appears to rigidify alpha40Lys in HbCO A, but does not significantly affect the flexibility of beta109Val in deoxy-Hb A. In the presence of IHP, several amino acid residues, especially those at the interdimer (alpha 1beta 2 or alpha 2beta 1) interface of HbCO A, exhibit significant conformational exchange. The affected residues include the proximal beta92His in the beta-heme pocket, as well as some other residues located in the flexible joint (betaC helix-alphaFG corner) and switch (alphaC helix-betaFG corner) regions that play an important role in the dimer-dimer rotation of Hb during the oxygenation process. These findings suggest that, upon IHP binding, HbCO A undergoes a conformational fluctuation near the R-state but biased toward the T-state, apparently along the trajectory of its allosteric transition, accompanied by structural fluctuations in the heme pocket of the beta-chain. In contrast, no significant perturbation of the dynamic features on the ms-micros time scale has been observed upon IHP binding to deoxy-Hb A. We propose that the allosteric effector-induced quaternary structural fluctuation may contribute to the reduced ligand affinity of ligated hemoglobin. Conformational exchange mapping of the beta-chain of HbCO A observed at 21.1 T shows significantly increased scatter in the chemical exchange contribution to the transverse relaxation rate ( R ex) values, relative to those at lower fields, due to the enhanced effect of the local chemical shift anisotropy (CSA) fluctuation. A spring-on-scissors model is proposed to interpret the dynamic phenomena induced by the heterotropic effector, IHP.

A quantum-chemical picture of hemoglobin affinity

Understanding the molecular mechanism of hemoglobin cooperativity remains an enduring challenge. Protein forces that control ligand affinity are not directly accessible by experiment. We demonstrate that computational quantum mechanics/molecular mechanics methods can provide reasonable values of ligand binding energies in Hb, and of their dependence on allostery. About 40% of the binding energy differences between the relaxed state and tense state quaternary structures result from strain induced in the heme and its ligands, especially in one of the pyrrole rings. The proximal histidine also contributes significantly, in particular, in the alpha-chains. The remaining energy difference resides in protein contacts, involving residues responsible for locking the quaternary changes. In the alpha-chains, the most important contacts involve the FG corner, at the "hinge" region of the alpha(1)beta(2) quaternary interface. The energy differences are spread more evenly among the beta-chain residues, suggesting greater flexibility for the beta- than for the alpha-chains along the quaternary transition. Despite this chain differentiation, the chains contribute equally to the relaxed substitute state energy difference. Thus, nature has evolved a symmetric response to the quaternary structure change, which is a requirement for maximum cooperativity, via different mechanisms for the two kinds of chains.

Molecular crowding inhibits intramolecular breathing motions in proteins

In aqueous solution some proteins undergo large-scale movements of secondary structures, subunits or domains, referred to as protein "breathing", that define a native-state ensemble of structures. These fluctuations are sensitive to the nature and concentration of solutes and other proteins and are thereby expected to be different in the crowded interior of a cell than in dilute solution. Here we use a combination of wide angle X-ray scattering (WAXS) and computational modeling to derive a quantitative measure of the spatial scale of conformational fluctuations in a protein solution. Concentration-dependent changes in the observed scattering intensities are consistent with a model of structural fluctuations in which secondary structures undergo rigid-body motions relative to one another. This motion increases with decreasing protein concentration or increasing temperature. Analysis of a set of five structurally and functionally diverse proteins reveals a diversity of kinetic behaviors. Proteins with multiple disulfide bonds exhibit little or no increase in breathing in dilute solutions. The spatial extent of structural fluctuations appears highly dependent on both protein structure and concentration and is universally suppressed at very high protein concentrations.

Insights into the solution structure of human deoxyhemoglobin in the absence and presence of an allosteric effector

We present a nuclear magnetic resonance (NMR) study in solution of the structures of human normal hemoglobin (Hb A) in the deoxy or unligated form in the absence and presence of an allosteric effector, inositol hexaphosphate (IHP), using 15N-1H residual dipolar coupling (RDC) measurements. There are several published crystal structures for deoxyhemoglobin A (deoxy-Hb A), and it has been reported that the functional properties of Hb A in single crystals are different from those in solution. Carbonmonoxyhemoglobin A (HbCO A) can also be crystallized in several structures. Our recent RDC studies of HbCO A in the absence and presence of IHP have shown that the solution structure of this Hb molecule is distinctly different from its classical crystal structures (R and R2). To have a better understanding of the structure-function relationship of Hb A under physiological conditions, we need to evaluate its structures in both ligated and unligated states in solution. Here, the intrinsic paramagnetic property of deoxy-Hb A has been exploited for the measurement of RDCs using the magnetic-field dependence of the apparent one-bond 1H-15N J couplings. Our RDC analysis suggests that the quaternary and tertiary structures of deoxy-Hb A in solution differ from its recently determined high-resolution crystal structures. Upon binding of IHP, structural changes in deoxy-Hb A are also observed, and these changes are largely within the alpha1beta1 (or alpha2beta2) dimer itself. These new structural findings allow us to gain a deeper insight into the structure-function relationship of this interesting allosteric protein.