New insights into the allosteric mechanism of human hemoglobin from molecular dynamics simulations

Identification of histidine residues that affect the T/R-state conformations of human hemoglobin using constant pH molecular dynamics simulations

Human hemoglobin (Hb) is a tetrameric protein consisting of two α and two β subunits that can adopt a low-affinity T- and high-affinity R-state conformations. Under physiological pH conditions, histidine (His) residues are the main sites for proton binding or release, and their protonation states can affect the T/R-state conformation of Hb. However, it remains unclear which His residues can effectively affect the Hb conformation. Herein, the impact of the 38 His residues of Hb on its T/R-state conformations was evaluated using constant-pH molecular dynamics (CpHMD) simulations at physiological pH while focusing on the His protonation states. Overall, the protonation states of some His residues were found to be correlated with the Hb conformation state. These residues were mainly located in the proximity of the heme (α87 and β92), and at the α1β2 and α2β1 interfaces (α89 and β97). This correlation may be partly explained by how easily hydrogen bonds can be formed, which depends on the protonation states of the His residues. Taken together, these CpHMD-based findings provide new insights into the identification of titratable His residues α87, α89, β92, and β97 that can affect Hb conformational switching under physiological pH conditions.

Slow dynamics measured by phosphorescence lifetime reveals global conformational changes in human adult hemoglobin induced by allosteric effectors

The mechanism underlying allostery in hemoglobin (Hb) is still not completely understood. Various models describing the action of allosteric effectors on Hb function have been published in the literature. It has also been reported that some allosteric effectors-such as chloride ions, inositol hexaphosphate, 2,3-diphospho-glycerate and bezafibrate-considerably lower the oxygen affinity of Hb. In this context, an important question is the extent to which these changes influence the conformational dynamics of the protein. Earlier, we elaborated a challenging method based on phosphorescence quenching, which makes characterizing protein-internal dynamics possible in the ms time range. The experimental technique involves phosphorescence lifetime measurements in thermal equilibrium at varied temperatures from 10 K up to 273 K, based on the signal of Zn-protoporphyrin substituted for the heme in the β-subunits of Hb. The thermal activation of protein dynamics was observed by the enhancement of phosphorescence quenching attributed to O2 diffusion. It was shown that the thermal activation of protein matrix dynamics was clearly distinguishable from the dynamic activation of the aqueous solvent, and was therefore highly specific for the protein. In the present work, the same method was used to study the changes in the parameters of the dynamic activation of human HbA induced by binding allosteric effectors. We interpreted the phenomenon as phase transition between two states. The fitting of this model to lifetime data yielded the change of energy and entropy in the activation process and the quenching rate in the dynamically activated state. The fitted parameters were particularly sensitive to the presence of allosteric effectors and could be interpreted in line with results from earlier experimental studies. The results suggest that allosteric effectors are tightly coupled to the dynamics of the whole protein, and thus underline the importance of global dynamics in the regulation of Hb function.

Molecular Dynamics of Hemoglobin Reveals Structural Alterations and Explains the Interactions Driving Sickle Cell Fibrillation

In sickle cell anemia, deoxyhemoglobin deforms RBCs by forming fibrils inside that disintegrate on oxygenation. We studied 100 ns long all-atom molecular dynamics (MD) for sickle and normal hemoglobin fibril models to understand this process, complemented by multiple 1 μs MD for a single tetramer of sickle and normal hemoglobin in deoxy and oxy states. We find that the presence of hydrophobic residues without a bulky side chain at β-6 in hemoglobin is the reason for the stability of the fibrils. Moreover, the free energy landscapes from MD of hemoglobin starting in the tensed (T) state capture the putative transition from T to relaxed (R) state, associated with oxygen binding. The three conformational wells in the landscapes are characterized by the quaternary changes where one αβ dimer rotates with respect to the other. The conformational changes from the oxygenation of sickle hemoglobin hinder the intermolecular contacts necessary for fibril formation.

Methemoglobin formation in mutant hemoglobin α chains: electron transfer parameters and rates

Hemoglobin-mediated transport of dioxygen (O₂) critically depends on the stability of the reduced (Fe²⁺) form of the heme cofactors. Some protein mutations stabilize the oxidized (Fe³⁺) state (methemoglobin, Hb M), causing methemoglobinemia, and can be lethal above 30%. The majority of the analyses of factors influencing Hb oxidation are retrospective and give insights only for inner-sphere mutations of heme (His58, His87). Herein, we report the first all-atom molecular dynamics simulations on both redox states and calculations of the Marcus electron transfer (ET) parameters for the α chain Hb oxidation and reduction rates for Hb M. The Hb wild-type (WT) and most of the studied α chain variants maintain globin structure except the Hb M Iwate (H87Y). The mutants forming Hb M tend to have lower redox potentials and thus stabilize the oxidized (Fe³⁺) state (in particular, the Hb Miyagi variant with K61E mutation). Solvent reorganization (λ_solv 73-96%) makes major contributions to reorganization free energy, whereas protein reorganization (λ_prot) accounts for 27-30% except for the Miyagi and J-Buda variants (λ_prot ∼4%). Analysis of heme-solvent H-bonding interactions among variants provide insights into the role of Lys61 residue in stabilizing the Fe²⁺ state. Semiclassical Marcus ET theory-based calculations predict experimental k_ET for the Cyt b5-Hb complex and provide insights into relative reduction rates for Hb M in Hb variants. Thus, our methodology provides a rationale for the effect of mutations on the structure, stability, and Hb oxidation reduction rates and has potential for identification of mutations that result in methemoglobinemia.

Computational and drug target analysis of functional single nucleotide polymorphisms associated with Haemoglobin Subunit Beta (HBB) gene

There is overwhelming evidence implicating Haemoglobin Subunit Beta (HBB) protein in the onset of beta thalassaemia. In this study for the first time, we used a combined SNP informatics and computer algorithms such as Neural network, Bayesian network, and Support Vector Machine to identify deleterious non-synonymous Single Nucleotide Polymorphisms (nsSNPs) present in the HBB gene. Our findings highlight three major mutation points (R31G, W38S, and Q128P) within the HBB gene sequence that have significant statistical and computational associations with the onset of beta thalassaemia. The dynamic simulation study revealed that R31G, W38S, and Q128P elicited high structural perturbation and instability, however, the wild type protein was considerably stable. Ten compounds with therapeutic potential against HBB were also predicted by structure-based virtual screening. Interestingly, the instability caused by the mutations was reversed upon binding to a ligand. This study has been able to predict potential deleterious mutants that can be further explored in the understanding of the pathological basis of beta thalassaemia and the design of tailored inhibitors.

Dissimilar flexibility of α and β subunits of human adult hemoglobin influences the protein dynamics and its alteration induced by allosteric effectors

The general question by what mechanism an "effector" molecule and the hemes of hemoglobin interact over widely separated intramolecular distances to change the oxygen affinity has been extensively investigated, and still has remained of central interest. In the present work we were interested in clarifying the general role of the protein matrix and its dynamics in the regulation of human adult hemoglobin (HbA). We used a spectroscopy approach that yields the compressibility (κ) of the protein matrix around the hemes of the subunits in HbA and studied how the binding of heterotropic allosteric effectors modify this parameter. κ is directly related to the variance of volume fluctuation, therefore it characterizes the molecular dynamics of the protein structure. For the experiments the heme groups either in the α or in the β subunits of HbA were replaced by fluorescent Zn-protoporphyrinIX, and series of fluorescence line narrowed spectra were measured at varied pressures. The evaluation of the spectra yielded the compressibility that showed significant dynamic asymmetry between the subunits: κ of the α subunit was 0.17±0.05/GPa, while for the β subunit it was much higher, 0.36±0.07/GPa. The heterotropic effectors, chloride ions, inositol hexaphosphate and bezafibrate did not cause significant changes in κ of the α subunits, while in the β subunits the effectors lead to a significant reduction down to 0.15±0.04/GPa. We relate our results to structural data, to results of recent functional studies and to those of molecular dynamics simulations, and find good agreements. The observed asymmetry in the flexibility suggests a distinct role of the subunits in the regulation of Hb that results in the observed changes of the oxygen binding capability.

Tertiary and quaternary structural basis of oxygen affinity in human hemoglobin as revealed by multiscale simulations

Human hemoglobin (Hb) is a benchmark protein of structural biology that shaped our view of allosterism over 60 years ago, with the introduction of the MWC model based on Perutz structures of the oxy(R) and deoxy(T) states and the more recent Tertiary Two-State model that proposed the existence of individual subunit states -“r” and “t”-, whose structure is yet unknown. Cooperative oxygen binding is essential for Hb function, and despite decades of research there are still open questions related to how tertiary and quaternary changes regulate oxygen affinity. In the present work, we have determined the free energy profiles of oxygen migration and for HisE7 gate opening, with QM/MM calculations of the oxygen binding energy in order to address the influence of tertiary differences in the control of oxygen affinity. Our results show that in the α subunit the low to high affinity transition is achieved by a proximal effect that mostly affects oxygen dissociation and is the driving force of the allosteric transition, while in the β subunit the affinity change results from a complex interplay of proximal and distal effects, including an increase in the HE7 gate opening, that as shown by free energy profiles promotes oxygen uptake.

The time scale of the quaternary structural changes in hemoglobin revealed using the transient grating technique

The quaternary structural transition between the R and T states of human hemoglobin was investigated using the transient grating technique. The results presented herein reveal that the quaternary structural change accompanied by the R-T transition occurs within a few microseconds.

Features of S-nitrosylation based on statistical analysis and molecular dynamics simulation: cysteine acidity, surrounding basicity, steric hindrance and local flexibility

S-Nitrosylation is involved in protein functional regulation and cellular signal transduction. Although intensive efforts have been made, the molecular mechanisms of S-nitrosylation have not yet been fully understood. In this work, we carried out a survey on 213 protein structures with S-nitrosylated cysteine sites and molecular dynamic simulations of hemoglobin as a case study. It was observed that the S-nitrosylated cysteines showed a lower pKa, a higher population of basic residues, a lower population of big-volume residues in the neighborhood, and relatively higher flexibility. The case study of hemoglobin showed that, compared to that in the T-state, Cysβ93 in the R-state hemoglobin possessed the above structural features, in agreement with the previous report that the R-state was more reactive in S-nitrosylation. Moreover, basic residues moved closer to the Cysβ93 in the dep-R-state hemoglobin, while big-volume residues approached the Cysβ93 in the dep-T-state. Using the four characteristics, i.e. cysteine acidity, surrounding basicity, steric hindrance, and local flexibility, a 3-dimensional model of S-nitrosylation was constructed to explain 61.9% of the S-nitrosylated and 58.1% of the non-S-nitrosylated cysteines. Our study suggests that cysteine deprotonation is a prerequisite for protein S-nitrosylation, and these characteristics might be useful in identifying specificity of protein S-nitrosylation.

Structural basis for cooperative oxygen binding and bracelet-assisted assembly of Lumbricus terrestris hemoglobin

The iron-containing hemoglobins (Hbs) are essential proteins to serve as oxygen transporters in the blood. Among various kinds of Hbs, the earthworm Hbs are the champions in carrying oxygen due to not only their large size but also the unusually high cooperativity of ligand binding. However, the cooperative oxygen binding mechanisms are still mostly unknown. Here we report the cryo-electron microscopy structure of Lumbricus terrestris Hb in its native, oxygenated state at 9.1 Å resolution, showing remarkable differences from the carbon monoxide-binding X-ray structure. Our structural analysis first indicates that the cooperative ligand binding of L. terrestris Hb requires tertiary and quaternary transitions in the heme pocket and a global subunit movement facilitated by intra-ring and inter-ring contacts. Moreover, the additional sinusoidal bracelet provides the confirmation for the long-standing debate about the additional electron densities absent in the X-ray crystal structure.

Non-site-specific allosteric effect of oxygen on human hemoglobin under high oxygen partial pressure

Protein allostery is essential for vital activities. Allosteric regulation of human hemoglobin (HbA) with two quaternary states T and R has been a paradigm of allosteric structural regulation of proteins. It is widely accepted that oxygen molecules (O₂) act as a “site-specific” homotropic effector, or the successive O₂ binding to the heme brings about the quaternary regulation. However, here we show that the site-specific allosteric effect is not necessarily only a unique mechanism of O₂ allostery. Our simulation results revealed that the solution environment of high O₂ partial pressure enhances the quaternary change from T to R without binding to the heme, suggesting an additional “non-site-specific” allosteric effect of O₂. The latter effect should play a complementary role in the quaternary change by affecting the intersubunit contacts. This analysis must become a milestone in comprehensive understanding of the allosteric regulation of HbA from the molecular point of view.

Collective dynamics underlying allosteric transitions in hemoglobin

Hemoglobin is the prototypic allosteric protein. Still, its molecular allosteric mechanism is not fully understood. To elucidate the mechanism of cooperativity on an atomistic level, we developed a novel computational technique to analyse the coupling of tertiary and quaternary motions. From Molecular Dynamics simulations showing spontaneous quaternary transitions, we separated the transition trajectories into two orthogonal sets of motions: one consisting of intra-chain motions only (referred to as tertiary-only) and one consisting of global inter-chain motions only (referred to as quaternary-only). The two underlying subspaces are orthogonal by construction and their direct sum is the space of full motions. Using Functional Mode Analysis, we were able to identify a collective coordinate within the tertiary-only subspace that is correlated to the most dominant motion within the quaternary-only motions, hence providing direct insight into the allosteric coupling mechanism between tertiary and quaternary conformation changes. This coupling-motion is substantially different from tertiary structure changes between the crystallographic structures of the T- and R-state. We found that hemoglobin's allosteric mechanism of communication between subunits is equally based on hydrogen bonds and steric interactions. In addition, we were able to affect the T-to-R transition rates by choosing different histidine protonation states, thereby providing a possible atomistic explanation for the Bohr effect.

How does hemoglobin generate such diverse functionality of physiological relevance?

The absolute values of the O2-affinities (P50, Klow, and Khigh) of hemoglobin (Hb) are regulated neither by changes in the static T-/R-quaternary and associated tertiary structures nor the ligation states. They are pre-determined and regulated by the extrinsic environmental factors such as pH, buffers, and heterotropic effectors. The effect and role of O2 on Hb are reversibly to drive the structural allosteric equilibrium between the T(deoxy)- and R(oxy)-Hb toward R(oxy)-Hb (the structural allostery). R(oxy)-Hb has a higher O2-affinity (Khigh) relative to that (Klow) of the T(deoxy)-Hb (Khigh>Klow) under any fixed environmental conditions. The apparent O2-affinity of Hb is high, as the globin matrix interferes with the dissociation process of O2, forcing the dissociated O2 geminately to re-bind to the heme Fe. This artificially increases [oxy-Hb] and concomitantly decreases [deoxy-Hb], leading to the apparent increases of the O2-affinity of Hb. The effector-linked high-frequency thermal fluctuations of the globin matrix act as a gating mechanism to modulate such physical, energetic, and kinetic barriers to enhance the dissociation process of O2, resulted in increases in [deoxy-Hb] and concomitant decrease in [oxy-Hb], leading to apparent reductions of the O2-affinity of Hb (the entropic allostery). The heme in Hb is simply a low-affinity O2-trap, the coordination structure of which is not altered by static T-/R-quaternary and associated tertiary structural changes of Hb. Thus, heterotrophic effectors are the signal molecule, which acts as a functional link between these two allosteries and generates the diverse functionality of Hb of physiological relevance. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Small ligand-globin interactions: reviewing lessons derived from computer simulation

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

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.

Unsuspected pathway of the allosteric transition in hemoglobin

Large conformational transitions play an essential role in the function of many proteins, but experiments do not provide the atomic details of the path followed in going from one end structure to the other. For the hemoglobin tetramer, the transition path between the unliganded (T) and tetraoxygenated (R) structures is not known, which limits our understanding of the cooperative mechanism in this classic allosteric system, where both tertiary and quaternary changes are involved. The conjugate peak refinement algorithm is used to compute an unbiased minimum energy path at atomic detail between the two end states. Although the results confirm some of the proposals of Perutz [Perutz MF (1970) Stereochemistry of cooperative effects in haemoglobin. Nature 228:726-734], the subunit motions do not follow the textbook description of a simple rotation of one αβ-dimer relative to the other. Instead, the path consists of two sequential quaternary rotations, each involving different subdomains and axes. The quaternary transitions are preceded and followed by phases of tertiary structural changes. The results explain the recent photodissociation measurements, which suggest that the quaternary transition has a fast (2 μs) as well as a slow (20 μs) component and provide a testable model for single molecule FRET experiments.

Spontaneous quaternary and tertiary T-R transitions of human hemoglobin in molecular dynamics simulation

We present molecular dynamics simulations of unliganded human hemoglobin (Hb) A under physiological conditions, starting from the R, R2, and T state. The simulations were carried out with protonated and deprotonated HC3 histidines His(β)146, and they sum up to a total length of 5.6µs. We observe spontaneous and reproducible T→R quaternary transitions of the Hb tetramer and tertiary transitions of the α and β subunits, as detected from principal component projections, from an RMSD measure, and from rigid body rotation analysis. The simulations reveal a marked asymmetry between the α and β subunits. Using the mutual information as correlation measure, we find that the β subunits are substantially more strongly linked to the quaternary transition than the α subunits. In addition, the tertiary populations of the α and β subunits differ substantially, with the β subunits showing a tendency towards R, and the α subunits showing a tendency towards T. Based on the simulation results, we present a transition pathway for coupled quaternary and tertiary transitions between the R and T conformations of Hb.

As the prototypic allosteric protein, human hemoglobin (Hb) has drawn extensive scientific efforts for many decades. Human Hb exists in two quaternary conformations, the low-affinity (or deoxy) T state, and the high-affinity (or oxy) R state, and the transition between the T and the R state is mainly characterized by a 15° rotation of the α1β1 dimer with respect to the α1β2 dimer. Subsequent binding of molecular oxygen to the four heme groups drives the Hb tetramer from the deoxy T to the oxy R state, rendering the T→R transitions the molecular process that underlies the well-known cooperativity of ligand binding. Despite the wealth of structural information available for Hb, the mechanistic coupling between the quaternary transition and the tertiary transitions in the individual subunits remains poorly understood. We report spontaneous and reproducible T-R transitions of Hb in molecular dynamics simulations, allowing us to study the mechanism underlying the transitions in atomistic detail. We pay special attention to the interplay between the tertiary and quaternary transitions of Hb, as well as to the tertiary t/r populations of the subunits in a particular (T or R) quaternary state. Interestingly, we observe a pronounced asymmetry between the α and β subunits with the β subunits being more strongly linked to the quaternary transitions than the α subunits. The simulations allow us to propose a pathway for coupled quaternary and tertiary transitions between the R and T conformations of Hb.

Effects of active site mutations in haemoglobin I from Lucina pectinata: a molecular dynamic study

Haemoglobin I from Lucina pectinata is a monomeric protein consisting of 142 amino acids. Its active site contains a peculiar arrangement of phenylalanine residues (PheB10, PheCD1 and PheE11) and a distal Gln at position E7. Active site mutations at positions B10, E7 and E11 were performed in deoxy haemoglobin I (HbI), followed by 10 ns molecular dynamic simulations. The results showed that the mutations induced changes in domains far from the active site producing more flexible structures than the native HbI. Distance analyses revealed that the heme pocket amino acids at positions E7 and B10 are extremely sensitive to any heme pocket residue mutation. The high flexibility observed by the E7 position suggests an important role in the ligand binding kinetics in ferrous HbI, while both positions play a major role in the ligand stabilisation processes. Furthermore, our results showed that E11Phe plays a pivotal role in protein stability.

Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes

Conformational transitions between open/closed or free/bound states in proteins possess functional importance. We propose a technique in which the collective modes obtained from an anisotropic network model (ANM) are used in conjunction with a Monte Carlo (MC) simulation approach, to investigate conformational transition pathways and pathway intermediates. The ANM-MC technique is applied to adenylate kinase (AK) and hemoglobin. The iterative method, in which normal modes are continuously updated during the simulation, proves successful in accomplishing the transition between open-closed conformations of AK and tense-relaxed forms of hemoglobin (C(alpha)-root mean square deviations between two end structures of 7.13 A and 3.55 A, respectively). Target conformations are reached by root mean-square deviations of 2.27 A and 1.90 A for AK and hemoglobin, respectively. The intermediate conformations overlap with crystal structures from the AK family within a 3.0-A root mean-square deviation. In the case of hemoglobin, the transition of tense-to-relaxed passes through the relaxed state. In both cases, the lowest-frequency modes are effective during transitions. The targeted Monte Carlo approach is used without the application of collective modes. Both the ANM-MC and targeted Monte Carlo techniques can explore sequences of events in transition pathways with an efficient yet realistic conformational search.

Protein dynamics explain the allosteric behaviors of hemoglobin

Bohr, Hasselbalch, and Krogh discovered homotropic and heterotropic allosteric behaviors of hemoglobin (Hb) in 1903/1904. A chronological description since then of selected principal models of the allosteric mechanism of Hb, such as the Adair scheme, the MWC two-state concerted model, the KNF induced-fit sequential model, the Perutz stereochemical model, the tertiary two-state model, and the global allostery model (an expanded MWC models), is concisely presented, followed by analysis and discussion of their limitations and deficiencies. The determination of X-ray crystallographic structures of deoxy- and ligated-Hb and the structure-based stereochemical model by Perutz are an epoch-making event in this history. However, his assignment of low-affinity deoxy- and high-affinity oxy-quaternary structures of Hb to the T- and R-states, respectively, though apparently reasonable, and as well as his hypothesis that the T-/R-quaternary structural transition regulates the oxygen-affinity, have created confusions and side-tracked studies of Hb on the structure-function relationship. The differences in static molecular structures of Hb between T(deoxy)- and R(oxy)-quaternary states reported in detail by Perutz and others are ligation-linked structural changes, but not related to the control/regulation of the oxygen-affinity. The oxygen-affinity (K(T) and K(R)) of Hb has been shown to be regulated by the heterotropic effector-linked tertiary structural changes without involving the T/R-quaternary changes. However, a recent high-resolution crystallographic analysis of Hb with different oxygen-affinities shows that static molecular structures of Hb determined by crystallography can neither identify the nature of the T(low-affinity) functional state nor decipher the mechanism by which Hb stores free energy in the T(low-affinity) functional state. Molecular dynamics simulations show that fluctuations of helices of oxy-Hb are increased upon de-oxygenation and/or binding 2,3-biphosphoglycerate. These are known to lower the oxygen-affinity of Hb. It is proposed that the coordination mode of the heme Fe with proximal and distal His is modulated by these helical fluctuations, resulting in the modulation of the oxygen-affinity of Hb. Therefore, it is proposed that the oxygen-affinity of Hb is regulated by pentanary (the 5th-order time-dependent or dynamic) tertiary structural changes rather than the T-/R-quaternary structural transitions in Hb. Homotropic and heterotropic allosteric effects of Hb are oxygen- and effector-linked, conformational entropy-driven entropy-enthalpy compensation phenomena and not much to do with static structural changes. The dynamic allostery model, which integrates these observations, provides the structural basis for the global allostery model (an expanded MWC model).

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.

Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics

Recent functional studies reported on human adult hemoglobin (HbA) show that heterotropic effector-linked tertiary structural changes are primarily responsible for modulating the oxygen affinity of hemoglobin. We present the results of 6-ns molecular dynamics simulations performed to gain insights into the dynamical and structural details of these effector-linked tertiary changes. All-atom simulations were carried out on a series of models generated for T- and R-state HbA, and for 2,3-diphosphoglycerate-bound models. Cross-correlation analyses identify both intra- and intersubunit correlated motions that are perturbed by the presence of the effector. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. It is found that HbA in both quaternary states exists as ensembles of tertiary conformations that introduce dynamic heterogeneity in the protein. 2,3-Diphosphoglycerate induces significant perturbations in the fluctuations of both HbA states that translate into the protein visiting different tertiary conformations within each quaternary state. The analysis reveals that the presence of the effector affects the most important components of HbA motions and that heterotropic effectors modify the overall dynamics of the quaternary equilibrium via tertiary changes occurring in regions where conserved functionally significant residues are located, namely in the loop regions between helices C and E, E and F, and F and G, and in concerted helix motions. The changes are not apparent when comparing the available x-ray crystal structures in the presence and absence of effector, but are striking when comparing the respective dynamic tertiary conformations of the R and T tetramers.

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.

Quaternary relaxations in sol-gel encapsulated hemoglobin studied via NIR and UV spectroscopy

In this work, we study the kinetics of the R --> T transition in hemoglobin using a combination of near-infrared and near-ultraviolet spectroscopy. We use a sol-gel encapsulation protocol to decelerate the conformational transitions and to avoid spectral perturbations arising from ligand migration and recombination. We monitor two spectroscopic markers: band III in the near-IR, which is a fine probe of the heme pocket conformation, and the tryptophan band in the near-UV, which probes the formation of the Trpbeta37-Aspalpha94 hydrogen bond, characteristic of the T structure, at the critical alpha1beta2 subunit interface. The time evolution of these two bands is monitored after deoxygenation of encapsulated oxyhemoglobin, obtained by diffusion of a reducing agent into the porous silica matrix. Characteristic spectral shifts are observed: comparison with myoglobin enables us to assign them to quaternary structure relaxations. Band III spectral relaxation is clearly nonexponential, and analysis with the Maximum Entropy Method enables us to identify three processes. On the other hand, near-UV spectral relaxation follows an exponential decay with a time constant closely corresponding to the second process observed in the near IR. Very interestingly, the rates of all processes markedly depend on the viscosity of the co-encapsulated solvent, following a power law. Our results reveal correlations between heme pocket relaxations, induced by the R --> T transition, and structural event(s) occurring at the alpha1beta2 interface and highlight their solvent dependence. The power law viscosity dependence of relaxation rates suggests that the observed protein relaxations are "slaved" to the co-encapsulated solvent. The stepwise character of the quaternary transition is also evidenced.

A DFT study on the relative affinity for oxygen of the alpha and beta subunits of hemoglobin

DFT calculations are carried out on computational models of the active center of the alpha and beta subunits of hemoglobin in both its oxygenated (R) and deoxygenated (T) states. The computational models are defined by the full heme group, including all porphyrin substituents, and the four amino acids closer to it. The role of the protein environment is introduced by freezing the position of the alpha carbon atom of each of the four amino acids to the positions they have in the available PDB structures. Oxygen affinity is then evaluated by computing the energy difference between the optimized structures of the oxygenated and deoxygenated forms of each model. The results indicate a higher affinity of the alpha subunits over the beta ones. Analysis of the computed structures points out to the strength of the hydrogen bond between the distal histidine and the oxygen molecule as a key factor in discriminating the different systems.

Molecular dynamics simulation of deoxy and carboxy murine neuroglobin in water

Globins are respiratory proteins that reversibly bind dioxygen and other small ligands at the iron of a heme prosthetic group. Hemoglobin and myoglobin are the most prominent members of this protein family. Unexpectedly a few years ago a new member was discovered and called neuroglobin (Ngb), being predominantly expressed in the brain. Ngb is a single polypeptide of 151 amino acids and despite the small sequence similarity with other globins, it displays the typical globin fold. Oxygen, nitric oxide, or carbon monoxide can displace the distal histidine which, in ferrous Ngb as well as in ferric Ngb, is bound to the iron, yielding a reversible adduct. Recent crystallographic data on carboxy Ngb show that binding of an exogenous ligand is associated to structural changes involving heme sliding and a topological reorganization of the internal cavities; in particular, the huge internal tunnel that connects the bulk with the active site, peculiar to Ngb, is heavily reorganized. We report the results of extended (90 ns) molecular dynamics simulations in water of ferrous deoxy and carboxy murine neuroglobin, which are both coordinated on the distal site, in the latter case by CO and in the former one by the distal His(64)(E7). The long timescale of the simulations allowed us to characterize the equilibrated protein dynamics and to compare protein structure and dynamical behavior coupled to the binding of an exogenous ligand. We have characterized the heme sliding motion, the topological reorganization of the internal cavities, the dynamics of the distal histidine, and particularly the conformational change of the CD loop, whose flexibility depends ligand binding.

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

Model-free-based NMR dynamics studies have been undertaken for polypeptide backbone amide N-H bond vectors for both the deoxy and carbonmonoxy forms of chain-specific, isotopically (15N and 2H) labeled tetrameric hemoglobin (Hb) using 15N-relaxation parameters [longitudinal relaxation rate (R1), transverse relaxation rate (R2), and heteronuclear nuclear Overhauser effect (NOE)] measured at two temperatures (29 and 34 degrees C) and two magnetic field strengths (11.7 and 14.1 T). In both deoxy and carbonmonoxy forms of human normal adult hemoglobin (Hb A), the amide N-H bonds of most amino acid residues are rigid on the fast time scale (nanosecond to picosecond), except for the loop regions and certain helix-helix connections. Although rigid in deoxy-Hb A, beta146His has been found to be free from restriction of its backbone motions in the CO form, presumably due to the rupture of its hydrogen bond/salt bridge network. We now have direct dynamics evidence for this structural transition of Hb in solution. While remarkably flexible in the deoxy state, alpha31Arg and beta123Thr, neighbors in the intradimer (alpha1beta1) interface, exhibit stiffening upon CO binding. These findings imply a role for alpha31Arg and beta123Thr in the intradimer communication but contradict the results from X-ray crystallography. We have also found that there is considerable flexibility in the intradimer (alpha1beta1) interface (i.e., B, G, and H helices and the GH corner) and possible involvement of several amino acid residues (e.g., alpha31Arg, beta3Leu, beta41Phe, beta123Thr, and beta146His) in the allosteric pathway. Several amino acid residues at the intradimer interfaces, such as beta109Val, appear to be involved in possible conformational exchange processes. The dynamic picture derived from the present study provides new insights into the traditional description of the stereochemical mechanism for the cooperative oxygenation of Hb A based on X-ray crystallographic results.

A 45-ns molecular dynamics simulation of hemoglobin in water by vectorizing and parallelizing COSMOS90 on the earth simulator: Dynamics of tertiary and quaternary structures

Molecular dynamics (MD) simulations of human adult hemoglobin (HbA) were carried out for 45 ns in water with all degrees of freedom including bond stretching and without any artificial constraints. To perform such large-scale simulations, one of the authors (M.S.) accelerated his own software COSMOS90 on the Earth Simulator by vectorization and parallelization. The dynamical features of HbA were investigated by evaluating root-mean-square deviations from the initial X-ray structure (an oxy T-state hemoglobin with PDB code: 1GZX) and root-mean-square fluctuations around the average structure from the simulation trajectories. The four subunits (alpha(1), alpha(2), beta(1), and beta(2)) of HbA maintained structures close to their respective X-ray structures during the simulations even though no constraints were applied to HbA in the simulations. Dimers alpha(1)beta(1) and alpha(2)beta(2) also maintained structures close to their respective X-ray structures while they moved relative to each other like two stacks of dumbbells. The distance between the two dimers (alpha(1)beta(1) and alpha(2)beta(2)) increased by 2 A (7.4%) in the initial 15 ns and stably fluctuated at the distance with the standard deviation 0.2 A. The relative orientation of the two dimers fluctuated between the initial X-ray angle -100 degrees and about -105 degrees with intervals of a few tens of nanoseconds.

Sulfide-binding hemoglobins: Effects of mutations on active-site flexibility

The dynamics of Hemoglobin I (HbI) from the clam Lucina pectinata, from wild-type sperm whale (SW) myoglobin, and from the L29F/H64Q/V68F triple mutant of SW, both unligated and bound to hydrogen sulfide (H2S), have been studied in molecular dynamics simulations. Features that account for differences in H2S affinity among the three have been examined. Our results verify the existence of an unusual heme rocking motion in unligated HbI that can promote the entrance of large ligands such as H2S. The FQF-mutant partially reproduces the amplitude and relative orientation of the motion of HbI's heme group. Therefore, besides introducing favorable electrostatic interactions with H2S, the three mutations in the distal pocket change the dynamic properties of the heme group. The active-site residues Gln-64(E7), Phe-43(CD1), and His-93(F8) are also shown to be more flexible in unligated HbI than in FQF-mutant and SW. Further contributions to H2S affinity come from differences in hydrogen bonding between the heme propionate groups and nearby amino acid residues.

Dynamics of Allostery in Hemoglobin: Roles of the Penultimate Tyrosine H bonds

The tyrosine residues adjacent to the C termini of the hemoglobin (Hb) subunits, alphaY140 and betaY145, are expected to play important structural roles, because the C termini are the loci of T-state quaternary salt-bridges, and because the tyrosine side-chains bridge the H and F helices via H bonds to the alphaV93 and betaV98 carbonyl groups. These roles have been investigated via measurements of oxygen binding, (1)H NMR spectra, resonance Raman (RR) spectra, and time-resolved resonance Raman (TR(3)) spectra on site mutants in which the Hcdots, three dots, centeredF H bonds are eliminated by replacing the tyrosine residues with phenylalanine. The TR(3) spectra confirm the hypothesis, based on TR(3) studies of wild-type Hb, that the Hcdots, three dots, centeredF H bonds break and then re-form during the sub-microsecond phase of the R-T quaternary transition. The TR(3) spectra support the inference from other mutational studies that the alphabeta dimers act as single dynamic units in this early phase, motions of the E and F helices being coupled tightly across the dimer interface. Formation of T quaternary contacts occurs at about the same rate in the mutants as in HbA. However, these contacts are weakened substantially by the Y/F substitutions. Equilibrium perturbations are apparent also, especially for the alpha-subunits, in which relaxation of the Fe-His bond, strengthening of the Acdots, three dots, centeredE interhelical H bond, and weakening of the "switch" quaternary contact in deoxyHb are all apparent. Structural effects are less marked for the beta-chain Y/F replacement, but the Bohr effect is reduced by 25%, indicating that the salt-bridge and H bond interactions of the adjacent C terminus are loosened. The alpha-chain replacement reduces the Bohr effect much more, consistent with the global perturbations detected by the structure probes.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.

Internal cavities and ligand passageways in human hemoglobin characterized by molecular dynamics simulations

Molecular dynamics simulations of the unliganded T state of human hemoglobin showed the existence of a spontaneous, very wide cavity on the distal side of the alpha subunit. This cavity consists of three tunnels spreading from the vicinity of the iron atom (the ligand binding site) to the surface of the subunit, constituting possible passageways for the entrance of the ligand. A fourth passageway was characterized due to the trajectory of water molecules entering or leaving the heme pocket. Analogous passages were observed in the beta subunits. They all appear and disappear dynamically, although some parts of them are more persistent along the trajectories. The most persistent regions within these tunnels correspond to all the xenon docking sites of human cytoglobin and to some of those of sperm whale and horse heart myoglobins and group I truncated hemoglobins.

Exploring the common dynamics of homologous proteins. Application to the globin family

We present a procedure to explore the global dynamics shared between members of the same protein family. The method allows the comparison of patterns of vibrational motion obtained by Gaussian network model analysis. After the identification of collective coordinates that were conserved during evolution, we quantify the common dynamics within a family. Representative vectors that describe these dynamics are defined using a singular value decomposition approach. As a test case, the globin heme-binding family is considered. The two lowest normal modes are shown to be conserved within this family. Our results encourage the development of models for protein evolution that take into account the conservation of dynamical features.

High frequency dynamics in hemoglobin measured by magnetic relaxation dispersion

The magnetic relaxation dispersion profiles for formate, acetate, and water protons are reported for aqueous solutions of hemoglobin singly and doubly labeled with a nitroxide and mercury(II) ion at cysteines at beta-93. Using two spin labels, one nuclear and one electron spin, a long intramolecular vector is defined between the two beta-93 positions in the protein. The paramagnetic contributions to the observed 1H spin-lattice relaxation rate constant are isolated from the magnetic relaxation dispersion profiles obtained on a dual-magnet apparatus that provides spectral density functions characterizing fluctuations sensed by intermoment dipolar interactions in the time range from the tens of microseconds to approximately 1 ps. Both formate and acetate ions are found to bind specifically within 5 angstroms of the beta-93 spin-label position and the relaxation dispersion has inflection points corresponding to correlation times of 30 ps and 4 ns for both ions. The 4-ns motion is identified with exchange of the anions from the site, whereas the 30-ps correlation time is identified with relative motions of the spin label and the bound anion in the protein environment close to beta-93. The magnetic field dependence of the paramagnetic contributions in both cases is well described by a simple Lorentzian spectral density function; no peaks in the spectral density function are observed. Therefore, the high frequency motions of the protein monitored by the intramolecular vector defined by the electron and nuclear spin are well characterized by a stationary random function of time. Attempts to examine long vector fluctuations by employing electron spin and nuclear spin double-labeling techniques did not yield unambiguous characterization of the high frequency motions of the vector between beta-93 positions on different chains.

Hemoglobin site-mutants reveal dynamical role of interhelical H-bonds in the allosteric pathway: time-resolved UV resonance Raman evidence for intra-dimer coupling

The dynamical effect of eliminating specific tertiary H-bonds in the hemoglobin (Hb) tetramer has been investigated by site-directed mutagenesis and time-resolved absorption and ultraviolet resonance Raman (UVRR) spectroscopy. The Trp alpha 14...Thr alpha 67 and Trp beta 15...Ser beta 72 H-bonds connect the A and E helices in the alpha and beta chains, and are proposed to break in the earliest protein intermediate (Rdeoxy) following photo-deligation of HbCO, along with a second pair of H-bonds involving tyrosine residues. Mutation of the acceptor residues Thr alpha 67 and Ser beta 72 to Val and Ala eliminates the A-E H-bonds, but has been shown to have no significant effect on ligand-binding affinity or cooperativity, or on spectroscopic markers of the T-state quaternary interactions. However, the mutations have profound and unexpected effects on the character of the Rdeoxy intermediate, and on the dynamics of the subsequent steps leading to the T state. Formation of the initial quaternary contact (RT intermediate) is accelerated, by an order of magnitude, but the locking-in of the T state is delayed by a factor of 2. These rate effects are essentially the same for either mutation, or for the double mutation, suggesting that the alpha beta dimer behaves as a mechanically coupled dynamical unit. Further evidence for intra-dimer coupling is provided by the Rdeoxy UVRR spectrum, in which either or both mutations eliminate the tyrosine difference intensity, although only tryptophan H-bonds are directly affected. A possible mechanism for mechanical coupling is outlined, involving transmission of forces through the alpha(1)beta(1) (and alpha(2)beta(2)) interface. The present observations establish that quaternary motions can occur on the approximately 100 ns time-scale. They show also that a full complement of interhelical H-bonds actually slows the initial quaternary motion in Hb, but accelerates the locking in of the T-contacts.

The role of dynamics in allosteric regulation

The biomolecular conformational changes often associated with allostery are, by definition, dynamic processes. Recent publications have disclosed the role of pre-existing equilibria of conformational substates in this process. In addition, the role of dynamics as an entropic carrier of free energy of allostery has been investigated. Recent work thus shows that dynamics is pivotal to allostery, and that it constitutes much more than just the move from the 'T'-state to the 'R'-state. Emerging computational studies have described the actual pathways of allosteric change.

Time-resolved Absorption and UV Resonance Raman Spectra Reveal Stepwise Formation of T Quaternary Contacts in the Allosteric Pathway of Hemoglobin

Hemoglobin undergoes a series of molecular changes on the nanosecond and microsecond time-scale following photodissociation of CO ligands. We have monitored these processes with a combination of transient absorption and resonance Raman (RR) spectroscopy. The latter have been acquired at higher data rates than previously available, thanks to kilohertz Ti:sapphire laser technology, with frequency-quadrupling into the ultraviolet. As a result of improved resolution of the UVRR time-course, a new intermediate has been identified in the pathway from the R (HbCO) to the T (deoxyHb) state. This intermediate is not detected via absorption transients, since the change in heme absorption is insignificant, but its lifetime agrees with a reported magnetic circular dichroism transient, which has been attributed to a quaternary tryptophan interaction. The new UVRR data allow elaboration of the allosteric pathway by establishing that the T-state quaternary contacts are formed in two well-separated steps, with time constants of 2.9 micros and 21 micros, instead of a single 20 micros process. The first step involves the "hinge" region contacts, as monitored by the Trp beta 37...Asp alpha 94 H-bond, while the second involves the "switch" region, as monitored by the Tyr alpha 42...Asp beta 99 H-bond. A working model for the allosteric pathway is presented.

Allosteric Changes in Protein Structure Computed by a Simple Mechanical Model: Hemoglobin T↔R2 Transition

Information on protein dynamics has been usually inferred from spectroscopic studies of parts of the proteins, or indirectly from the comparison of the conformations assumed in the presence of different substrates or ligands. While molecular simulations also provide information on protein dynamics, they usually suffer from incomplete sampling of conformational space, and become prohibitively expensive when exploring the collective dynamics of large macromolecular structures. Here, we explore the dynamics of a well-studied allosteric protein, hemoglobin (Hb), to show that a simple mechanical model based on Gaussian fluctuations of residues can efficiently predict the transition between the tense (T, unliganded) and relaxed (R or R2, O(2) or CO-bound) forms of Hb. The passage from T into R2 is shown to be favored by the global mode of motion, which, in turn is driven by entropic effects. The major difference between the dynamics of the T and R2 forms is the loss of the hinge-bending role of alpha(1)-beta(2) (or alpha(2)-beta(1)) interfacial residues at alpha Phe36-His45 and beta Thr87-Asn102 in the R2 form, which implies a decreased cooperativity in the higher affinity (R2) form of Hb, consistent with many experimental studies. The involvement of the proximal histidine beta His92 in this hinge region suggests that the allosteric propagation of the local structural changes (induced upon O(2) binding) into global ones occur via hinge regions. This is the first demonstration that there is an intrinsic tendency of Hb to undergo T-->R2 transition, induced by purely elastic forces of entropic origin that are uniquely defined for the particular contact topology of the T form.

Molecular dynamics of the FixJ receiver domain: movement of the beta4-alpha4 loop correlates with the in and out flip of Phe101

FixJ is a two-domain response regulator involved in nitrogen fixation in Sinorhizobium meliloti. Recent X-ray characterization of both the native (unphosphorylated) and the active (phosphorylated) states of the protein identify conformational changes of the beta4-alpha4 loop and the conserved residue Phe101 as the key switches in activation. These structures also allowed investigation of the transition between conformations of this two-component regulatory receiver domain by molecular dynamics simulations. The path for the conformational change was studied with a distance constraint directing the system from one state to the other. The simulations provide evidence for a correlation between the conformation of the beta4-alpha4 loop and the orientation of the residue Phe101. A model presenting the sequence of events during the activation/deactivation process is discussed.

NMR investigation of the dynamics of tryptophan side-chains in hemoglobins

NMR relaxation measurements of 15N spin-lattice relaxation rate (R(1)), spin-spin relaxation rate (R(2)), and heteronuclear nuclear Overhauser effect (NOE) have been carried out at 11.7T and 14.1T as a function of temperature for the side-chains of the tryptophan residues of 15N-labeled and/or (2H,15N)-labeled recombinant human normal adult hemoglobin (Hb A) and three recombinant mutant hemoglobins, rHb Kempsey (betaD99N), rHb (alphaY42D/betaD99N), and rHb (alphaV96W), in the carbonmonoxy and the deoxy forms as well as in the presence and in the absence of an allosteric effector, inositol hexaphosphate (IHP). There are three Trp residues (alpha14, beta15, and beta37) in Hb A for each alphabeta dimer. These Trp residues are located in important regions of the Hb molecule, i.e. alpha14Trp and beta15Trp are located in the alpha(1)beta(1) subunit interface and beta37Trp is located in the alpha(1)beta(2) subunit interface. The relaxation experiments show that amino acid substitutions in the alpha(1)beta(2) subunit interface can alter the dynamics of beta37Trp. The transverse relaxation rate (R(2)) for beta37Trp can serve as a marker for the dynamics of the alpha(1)beta(2) subunit interface. The relaxation parameters of deoxy-rHb Kemspey (betaD99N), which is a naturally occurring abnormal human hemoglobin with high oxygen affinity and very low cooperativity, are quite different from those of deoxy-Hb A, even in the presence of IHP. The relaxation parameters for rHb (alphaY42D/betaD99N), which is a compensatory mutant of rHb Kempsey, are more similar to those of Hb A. In addition, TROSY-CPMG experiments have been used to investigate conformational exchange in the Trp residues of Hb A and the three mutant rHbs. Experimental results indicate that the side-chain of beta37Trp is involved in a relatively slow conformational exchange on the micro- to millisecond time-scale under certain experimental conditions. The present results provide new dynamic insights into the structure-function relationship in hemoglobin.