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 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.

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.

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.

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).

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.