Role of the Subunit Interactions in the Conformational Transitions in Adult Human Hemoglobin: An Explicit Solvent Molecular Dynamics Study

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.

Water dynamics around T0 vs R4 of hemoglobin from local hydrophobicity analysis

The local hydration around tetrameric hemoglobin (Hb) in its T₀ and R₄ conformational substates is analyzed based on molecular dynamics simulations. Analysis of the local hydrophobicity (LH) for all residues at the α₁β₂ and α₂β₁ interfaces, responsible for the quaternary T → R transition, which is encoded in the Monod-Wyman-Changeux model, as well as comparison with earlier computations of the solvent accessible surface area, makes clear that the two quantities measure different aspects of hydration. Local hydrophobicity quantifies the presence and structure of water molecules at the interface, whereas "buried surface" reports on the available space for solvent. For simulations with Hb frozen in its T₀ and R₄ states, the correlation coefficient between LH and buried surface is 0.36 and 0.44, respectively, but it increases considerably if the 95% confidence interval is used. The LH with Hb frozen and flexible changes little for most residues at the interfaces but is significantly altered for a few select ones: Thr41α, Tyr42α, Tyr140α, Trp37β, Glu101β (for T₀) and Thr38α, Tyr42α, Tyr140α (for R₄). The number of water molecules at the interface is found to increase by ∼25% for T₀ → R₄, which is consistent with earlier measurements. Since hydration is found to be essential to protein function, it is clear that hydration also plays an essential role in allostery.

Dynamics of camel and human hemoglobin revealed by molecular simulations

Hemoglobin is one of the most widely studied proteins genetically, biochemically, and structurally. It is an oxygen carrying tetrameric protein that imparts the characteristic red color to blood. Each chain of hemoglobin harbors a heme group embedded in a hydrophobic pocket. Several studies have investigated structural variations present in mammalian hemoglobin and their functional implications. However, camel hemoglobin has not been thoroughly explored, especially from a structural perspective. Importantly, very little is known about how the heme group interacts with hemoglobin under varying conditions of osmolarity and temperature. Several experimental studies have indicated that the tense (T) state is more stable than the relaxed (R) state of hemoglobin under normal physiological conditions. Despite the fact that R state is less stable than the T state, no extensive structural dynamics studies have been performed to investigate global quaternary transitions of R state hemoglobin under normal physiological conditions. To evaluate this, several 500 ns all-atom molecular dynamics simulations were performed to get a deeper understanding of how camel hemoglobin behaves under stress, which it is normally exposed to, when compared to human hemoglobin. Notably, camel hemoglobin was more stable under physiological stress when compared to human hemoglobin. Additionally, when compared to camel hemoglobin, cofactor-binding regions of hemoglobin also exhibited more fluctuations in human hemoglobin under the conditions studied. Several differences were observed between the residues of camel and human hemoglobin that interacted with heme. Importantly, distal residues His58 of α hemoglobin and His63 of β hemoglobin formed more sustained interactions, especially at higher temperatures, in camel hemoglobin. These residues are important for oxygen binding to hemoglobin. Thus, this work provides insights into how camel and human hemoglobin differ in their interactions under stress.

The functional role of the hemoglobin-water interface

The interface between hemoglobin (Hb) and its environment, in particular water, is of great physiological relevance. Here, results from in vitro, in vivo, and computational experiments (molecular dynamics simulations) are summarized and put into perspective. One of the main findings from the computations is that the stability of the deoxy, ligand-free T-state (T₀) can be stabilized relative to the deoxy R-state (R₀) only in sufficiently large simulation boxes for the hydrophobic effect to manifest itself. This effect directly influences protein stability and is operative also under physiological conditions. Furthermore, molecular simulations provide a dynamical interpretation of the Perutz model for Hb function. Results from experiments using higher protein concentrations and realistic cellular environments are also discussed. One of the next great challenges for computational studies, which as we show is likely to be taken up in the near future, is to provide a molecular-level understanding of the dynamics of proteins in such crowded environments.

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.

Water Dynamics Around Proteins: T- and R-States of Hemoglobin and Melittin

The water dynamics, as characterized by the local hydrophobicity (LH), is investigated for tetrameric hemoglobin (Hb) and dimeric melittin. For the T₀ to R₀ transition in Hb, it is found that LH provides additional molecular-level insight into the Perutz mechanism, i.e., the breaking and formation of salt bridges at the α₁/β₂ and α₂/β₁ interface is accompanied by changes in LH. For Hb in cubic water boxes with 90 and 120 Å edge length it is observed that following a decrease in LH as a consequence of reduced water density or change of water orientation at the protein/water interface the α/β interfaces are destabilized; this is a hallmark of the Perutz stereochemical model for the T to R transition in Hb. The present work thus provides a dynamical view of the classical structural model relevant to the molecular foundations of Hb function. For dimeric melittin, earlier results by Cheng and Rossky [ Nature 1998, 392, 696-699] are confirmed and interpreted on the basis of LH from simulations in which the protein structure is frozen. For the flexible melittin dimer, the changes in the local hydration can be as much as 30% greater than for the rigid dimer, reflecting the fact that protein and water dynamics are coupled.

Adaptations to environmental change: Globin superfamily evolution in Antarctic fishes

The ancient origins and functional versatility of globins make them ideal subjects for studying physiological adaptation to environmental change. Our goals in this review are to describe the evolution of the vertebrate globin gene superfamily and to explore the structure/function relationships of hemoglobin, myoglobin, neuroglobin and cytoglobin in teleost fishes. We focus on the globins of Antarctic notothenioids, emphasizing their adaptive features as inferred from comparisons with human proteins. We dedicate this review to Guido di Prisco, our co-author, colleague, friend, and husband of C.V. Ever thoughtful, creative, and enthusiastic, Guido spearheaded study of the structure, function, and evolution of the hemoglobins of polar fishes - this review is testimony to his wide-ranging contributions. Throughout his career, Guido inspired younger scientists to embrace polar biological research, and he challenged researchers of all ages to explore evolutionary adaptation in the context of global climate change. Beyond his scientific contributions, we will miss his warmth, his culture, and his great intellect. Guido has left an outstanding legacy, one that will continue to inspire us and our research.

Response to comment on 'Valid molecular dynamics simulations of human hemoglobin require a surprisingly large box size'

We recently reported that molecular dynamics simulations for hemoglobin require a surprisingly large box size to stabilize the T(0) state relative to R(0), as observed in experiments (El Hage et al., 2018). Gapsys and de Groot have commented on this work but do not provide convincing evidence that the conclusions of El Hage et al., 2018 are incorrect. Here we respond to these concerns, argue that our original conclusions remain valid, and raise our own concerns about some of the results reported in the comment by Gapsys and de Groot that require clarification.

Bohr effect and oxygen affinity of carp, eel and human hemoglobin: Quantitative analyses provide rationale for the Root effect

The functional properties of most fish hemoglobins are more complex than those of human hemoglobin. This complexity arises in the form of the Root effect, in which the oxygen affinity of such fish hemoglobins decreases rapidly with pH relative to that of human hemoglobin. Cooperative ligand binding is also diminished below pH ≈ 6.5. The Bohr effect, determined by acid-base titration, has been reported for the Root effect carp and anodic eel hemoglobins. Unlike for mammalian hemoglobins, the Wyman equation for the Bohr effect fails to account quantitatively for these Bohr data. We present a successful quantitative accounting for these data based on evidence for multiple T states in various fish hemoglobins and on their lack of sixhistidine Bohr groups, with pK^oxy > pK^deoxy. On the same bases we also provide a rationale for the higher pH sensitivity of the oxygen affinity of carp compared to human hemoglobin.

Valid molecular dynamics simulations of human hemoglobin require a surprisingly large box size

Recent molecular dynamics (MD) simulations of human hemoglobin (Hb) give results in disagreement with experiment. Although it is known that the unliganded (T[Formula: see text]) and liganded (R[Formula: see text]) tetramers are stable in solution, the published MD simulations of T[Formula: see text] undergo a rapid quaternary transition to an R-like structure. We show that T[Formula: see text] is stable only when the periodic solvent box contains ten times more water molecules than the standard size for such simulations. The results suggest that such a large box is required for the hydrophobic effect, which stabilizes the T[Formula: see text] tetramer, to be manifested. Even in the largest box, T[Formula: see text] is not stable unless His146 is protonated, providing an atomistic validation of the Perutz model. The possibility that extra large boxes are required to obtain meaningful results will have to be considered in evaluating existing and future simulations of a wide range of systems.

Comparative Study of Elastic Network Model and Protein Contact Network for Protein Complexes: The Hemoglobin Case

The overall topology and interfacial interactions play key roles in understanding structural and functional principles of protein complexes. Elastic Network Model (ENM) and Protein Contact Network (PCN) are two widely used methods for high throughput investigation of structures and interactions within protein complexes. In this work, the comparative analysis of ENM and PCN relative to hemoglobin (Hb) was taken as case study. We examine four types of structural and dynamical paradigms, namely, conformational change between different states of Hbs, modular analysis, allosteric mechanisms studies, and interface characterization of an Hb. The comparative study shows that ENM has an advantage in studying dynamical properties and protein-protein interfaces, while PCN is better for describing protein structures quantitatively both from local and from global levels. We suggest that the integration of ENM and PCN would give a potential but powerful tool in structural systems biology.

Induced Fit in Protein Multimerization: The HFBI Case

Hydrophobins, produced by filamentous fungi, are small amphipathic proteins whose biological functions rely on their unique surface-activity properties. Understanding the mechanistic details of the multimerization process is of primary importance to clarify the interfacial activity of hydrophobins. We used free energy calculations to study the role of a flexible β-hairpin in the multimerization process in hydrophobin II from Trichoderma reesei (HFBI). We characterized how the displacement of this β-hairpin controls the stability of the monomers/dimers/tetramers in solution. The regulation of the oligomerization equilibrium of HFBI will necessarily affect its interfacial properties, fundamental for its biological function and for technological applications. Moreover, we propose possible routes for the multimerization process of HFBI in solution. This is the first case where a mechanism by which a flexible loop flanking a rigid patch controls the protein-protein binding equilibrium, already known for proteins with charged binding hot-spots, is described within a hydrophobic patch.

Fungi proliferate by creating a complex hyphal network growing within a wet environment. However, for most fungi to colonize new territories, they must produce spores carried by aerial hyphae and spread them into the air. Aerial structures need to overcome the surface tension of the surrounding water in order to grow into the air. This process requires hydrophobins, a remarkable class of self-associating fungal proteins which lower the surface tension at the air/water interface by creating a thin amphipathic layer. In solution they form multimers in equilibrium with the interfacial layer. Due to their unique surface-activity properties, hydrophobins have been used for a variety of biotechnical applications. We used enhanced sampling molecular dynamics simulations methods to study the multimerization process in solution of a hydrophobin from Trichoderma reesei (HFBI). We clarified the fundamental role of a small flexible region within the HFBI monomer involved in the formation of multimers. A flexible loop flanking a rigid interaction patch is able to fine-tune the interaction energy. This mechanism, already known for charged binding patches, is described here for hydrophobic hot-spots. This result is remarkably important in order to clarify the mechanism of arranging at the interface and enhancing hydrophobin-based technological applications.

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.

Proteomic analysis of the mouse brain after repetitive exposure to hypoxia

Hypoxic preconditioning (HPC) is known to have a protective effect against hypoxic damage; however, the precise mechanisms involved remain unknown. In this study, an acute and repetitive hypoxia mouse model, two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) coupled with matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF/TOF-MS), and Western blot experiments were used to identify the differential expression of key proteins in the mouse brain during HPC. Approximately 2100 2D-DIGE spots were observed following gel imaging and spot detection. Significant differences (p < 0.05) in the expression of 66 proteins were observed between the 3× HPC treatment group and the control group, 45 proteins were observed between the 6× HPC treatment group and the control group, and 70 proteins were observed between the 3× HPC treatment group and the 6× HPC group. Consistent results among Western blot, 2D-DIGE and MS methods were observed for the proteins, ATP synthase subunit alpha, malate dehydrogenase, guanine nucleotide-binding protein subunit beta-1 and proteasome subunit alpha type-2. The proteins associated with ATP synthesis and the citric acid cycle were down-regulated, while those linked to glycolysis and oxygen-binding were up-regulated. This proteomic analysis of the mouse brain after HPC furthers understanding of the molecular pathways involved in the protective effect of HPC and these findings provide new insight into the mechanisms of hypoxia and HPC.

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.

Mapping polymerization and allostery of hemoglobin S using point mutations

Hemoglobin is a complex system that undergoes conformational changes in response to oxygen, allosteric effectors, mutations, and environmental changes. Here, we study allostery and polymerization of hemoglobin and its variants by application of two previously described methods: (i) AllosMod for simulating allostery dynamics given two allosterically related input structures and (ii) a machine-learning method for dynamics- and structure-based prediction of the mutation impact on allostery (Weinkam et al. J. Mol. Biol. 2013, 425, 647-661), now applicable to systems with multiple coupled binding sites, such as hemoglobin. First, we predict the relative stabilities of substates and microstates of hemoglobin, which are determined primarily by entropy within our model. Next, we predict the impact of 866 annotated mutations on hemoglobin's oxygen binding equilibrium. We then discuss a subset of 30 mutations that occur in the presence of the sickle cell mutation and whose effects on polymerization have been measured. Seven of these HbS mutations occur in three predicted druggable binding pockets that might be exploited to directly inhibit polymerization; one of these binding pockets is not apparent in the crystal structure, but only in structures generated by AllosMod. For the 30 mutations, we predict that mutation-induced conformational changes within a single tetramer tend not to significantly impact polymerization; instead, these mutations more likely impact polymerization by directly perturbing a polymerization interface. Finally, our analysis of allostery allows us to hypothesize why hemoglobin evolved to have multiple subunits and a persistent low frequency sickle cell mutation.

Subunit disassembly pathway of human hemoglobin revealing the site-specific role of its cysteine residues

Cysteine residues play a unique role in human hemoglobin (Hb) by affecting its cooperative oxygen binding behavior and the stability of its tetrameric structure. However, how these cysteine residues fulfill their biophysical functions from the molecular level is yet unclear. Here we study the subunit disassembly pathway of human hemoglobin using the sulfhydryl reagent, p-hydroxymercuribenzoate (PMB) and investigate the functional roles of cysteine residues in human hemoglobin. We show evidence from the matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry that all three types of cysteine residues, including the surface-exposed βCys93 and the shielded αCys104 and βCys112 are reactive to PMB, resolving an issue long under debate. It is demonstrated that all three types of cysteine residues must be blocked by PMB to accomplish the subunit disassembly, and the PMB-cysteine reactions proceed in a stepwise manner with an order of βCys93, αCys104, and βCys112. The PMB reactions with the three different cysteine residues demonstrate strong site-specificity. The possible influence of PMB-cysteine reactions to the stability of various intersubunit salt bridges has been discussed based on the crystallographic structure of hemoglobin, providing insights in understanding the hemoglobin subunit disassembly pathway and the site-specific functional role of each cysteine residue in hemoglobin.