Global allostery model of hemoglobin. Modulation of O(2) affinity, cooperativity, and Bohr effect by heterotropic allosteric effectors

Carbon Monoxide Poisoning: From Microbes to Therapeutics

Carbon monoxide (CO) poisoning leads to 50,000-100,000 emergency room visits and 1,500-2,000 deaths each year in the United States alone. Even with treatment, survivors often suffer from long-term cardiac and neurocognitive deficits, highlighting a clear unmet medical need for novel therapeutic strategies that reduce morbidity and mortality associated with CO poisoning. This review examines the prevalence and impact of CO poisoning and pathophysiology in humans and highlights recent advances in therapeutic strategies that accelerate CO clearance and mitigate toxicity. We focus on recent developments of high-affinity molecules that take advantage of the uniquely strong interaction between CO and heme to selectively bind and sequester CO in preclinical models. These scavengers, which employ heme-binding scaffolds ranging from organic small molecules to hemoproteins derived from humans and potentially even microorganisms, show promise as field-deployable antidotes that may rapidly accelerate CO clearance and improve outcomes for survivors of acute CO poisoning.

Mathematical models describing oxygen binding by hemoglobin

Despite the fact that the investigation of the structural and functional properties of hemoglobin dates back more than 150 years, the topic has not lost its relevance today. The most important component of these studies is the development of mathematical models that formalize and generalize the mechanisms determining the cooperative binding of ligands based on data on the structural and functional state of the protein. In this work, we review the mathematical relationships describing oxygen binding by hemoglobin, ranging from the classical Hüfner, Hill, and Adair equations to the Szabo-Karplus and tertiary two-state mathematical models based on the Monod-Wyman-Changeux and Koshland-Némethy-Filmer concepts. The generality of the considered equations as mathematical functions, bearing in their basis a power dependence, is demonstrated. The problems and possible solutions related to approximation of experimental data by the oxygenation equations with correlated fitting parameters are noted. Attention is paid to empirical equations, extended versions of the Hill equation, where the coefficient of cooperation is modulated by Gauss and Lorentz distributions as functions of partial oxygen pressure.

Cooperative oxygen binding in beta-semihemoglobins caused by a chemical modification in the alpha1beta1 interface

A beta-semihemoglobin is an alpha-beta dimer of hemoglobin (Hb) in which the beta-subunit carries heme, while the alpha-subunit is heme-less, in apo form. It is characterised by displaying a high affinity for oxygen, and absence of cooperative binding of oxygen. We have modified chemically the residue beta112Cys (G14), located adjacent to the alpha1beta1 interface, and studied the impact of such a modification on the oligomeric state and oxygenation properties of the derivatives. We also studied the impact of modifying beta93Cys (F9) since its modification was unavoidable. For this, we used N-Ethyl maleimide and iodoacetamide. For the alkylation of beta112Cys (G14) in isolated subunits, we used N-Ethyl maleimide, iodoacetamide, or additionally, 4,4'-Dithiopyridine. Seven native and chemically modified beta-subunit derivatives were prepared and analysed. Only those derivatives treated with iodoacetamide showed oxygenation properties that were indistinguishable from those of native beta-subunits. These derivatives were then converted into their respective semihemoglobin forms, and four additional derivatives were prepared and analysed .in terms of ligation-linked oligomeric state, and oxygenation function, and contrasted against native Hb and unmodified beta-subunits. Strikingly, beta-semiHbs with modifications in beta112Cys showed indications of cooperative oxygen binding in various degrees, which suggested the possibility of assembly of two beta-semiHbs. The derivative modified with 4-Thiopyridine in beta112Cys showed a highly cooperative binding of oxygen (nmax = 1.67). A plausible allosteric scheme that could explain allostery in beta-semiHb system is suggested.

X-ray crystallography and sickle cell disease drug discovery-a tribute to Donald Abraham

X-ray crystallography and structure-based drug discovery have played a major role in the discovery of antisickling agents that target hemoglobin (Hb) for the treatment of sickle cell disease (SCD). Sickle cell disease, the most common inherited hematologic disorder, occurs as a result of a single point mutation of βGlu6 in normal human adult hemoglobin (HbA) to βVal6 in sickle hemoglobin (HbS). The disease is characterized by polymerization of HbS and sickling of red blood cells (RBCs), leading to several secondary pathophysiologies, including but not limited to vaso-occlusion, hemolytic anemia, oxidative stress, inflammation, stroke, pain crisis, and organ damage. Despite the fact that SCD was the first disease to have its molecular basis established, the development of therapies was for a very long time a challenge and took several decades to find therapeutic agents. The determination of the crystal structure of Hb by Max Perutz in the early 60s, and the pioneering X-ray crystallography research by Donald J. Abraham in the early 80s, which resulted in the first structures of Hb in complex with small molecule allosteric effectors of Hb, gave much hope that structure-based drug discovery (SBDD) could be used to accelerate development of antisickling drugs that target the primary pathophysiology of hypoxia-induced HbS polymerization to treat SCD. This article, which is dedicated to Donald J. Abraham, briefly reviews structural biology, X-ray crystallography and structure-based drug discovery from the perspective of Hb. The review also presents the impact of X-ray crystallography in SCD drug development using Hb as a target, emphasizing the major and important contributions by Don Abraham in this field.

Resveratrol, a New Allosteric Effector of Hemoglobin, Enhances Oxygen Supply Efficiency and Improves Adaption to Acute Severe Hypoxia

Acute altitude hypoxia represents the cause of multiple adverse consequences. Current treatments are limited by side effects. Recent studies have shown the protective effects of resveratrol (RSV), but the mechanism remains unknown. To address this, the effects of RSV on the structure and function of hemoglobin of adult (HbA) were preliminarily analyzed using surface plasmon resonance (SPR) and oxygen dissociation assays (ODA). Molecular docking was conducted to specifically analyze the binding regions between RSV and HbA. The thermal stability was characterized to further validate the authenticity and effect of binding. Changes in the oxygen supply efficiency of HbA and rat RBCs incubated with RSV were detected ex vivo. The effect of RSV on the anti-hypoxic capacity under acute hypoxic conditions in vivo was evaluated. We found that RSV binds to the heme region of HbA following a concentration gradient and affects the structural stability and rate of oxygen release of HbA. RSV enhances the oxygen supply efficiency of HbA and rat RBCs ex vivo. RSV prolongs the tolerance times of mice suffering from acute asphyxia. By enhancing the oxygen supply efficiency, it alleviates the detrimental effects of acute severe hypoxia. In conclusion, RSV binds to HbA and regulates its conformation, which enhances oxygen supply efficiency and improves adaption to acute severe hypoxia.

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.

Carbon Monoxide Signaling: Examining Its Engagement with Various Molecular Targets in the Context of Binding Affinity, Concentration, and Biologic Response

Carbon monoxide (CO) has been firmly established as an endogenous signaling molecule with a variety of pathophysiological and pharmacological functions, including immunomodulation, organ protection, and circadian clock regulation, among many others. In terms of its molecular mechanism(s) of action, CO is known to bind to a large number of hemoproteins with at least 25 identified targets, including hemoglobin, myoglobin, neuroglobin, cytochrome c oxidase, cytochrome P450, soluble guanylyl cyclase, myeloperoxidase, and some ion channels with dissociation constant values spanning the range of sub-nM to high μM. Although CO's binding affinity with a large number of targets has been extensively studied and firmly established, there is a pressing need to incorporate such binding information into the analysis of CO's biologic response in the context of affinity and dosage. Especially important is to understand the reservoir role of hemoglobin in CO storage, transport, distribution, and transfer. We critically review the literature and inject a sense of quantitative assessment into our analyses of the various relationships among binding affinity, CO concentration, target occupancy level, and anticipated pharmacological actions. We hope that this review presents a picture of the overall landscape of CO's engagement with various targets, stimulates additional research, and helps to move the CO field in the direction of examining individual targets in the context of all of the targets and the concentration of available CO. We believe that such work will help the further understanding of the relationship of CO concentration and its pathophysiological functions and the eventual development of CO-based therapeutics. SIGNIFICANCE STATEMENT: The further development of carbon monoxide (CO) as a therapeutic agent will significantly rely on the understanding of CO's engagement with therapeutically relevant targets of varying affinity. This review critically examines the literature by quantitatively analyzing the intricate relationships among targets, target affinity for CO, CO level, and the affinity state of carboxyhemoglobin and provide a holistic approach to examining the molecular mechanism(s) of action for CO.

Structural origin of cooperativity in human hemoglobin: a view from different roles of α and β subunits in the α2β2 tetramer

This mini-review, mainly based on our resonance Raman studies on the structural origin of cooperative O₂ binding in human adult hemoglobin (HbA), aims to answering why HbA is a tetramer consisting of two α and two β subunits. Here, we focus on the Fe-His bond, the sole coordination bond connecting heme to a globin. The Fe-His stretching frequencies reflect the O₂ affinity and also the magnitude of strain imposed through globin by inter-subunit interactions, which is the origin of cooperativity. Cooperativity was first explained by Monod, Wyman, and Changeux, referred to as the MWC theory, but later explained by the two tertiary states (TTS) theory. Here, we related the higher-order structures of globin observed mainly by vibrational spectroscopy to the MWC theory. It became clear from the recent spectroscopic studies, X-ray crystallographic analysis, and mutagenesis experiments that the Fe-His bonds exhibit different roles between the α and β subunits. The absence of the Fe-His bond in the α subunit in some mutant and artificial Hbs inhibits T to R quaternary structural change upon O₂ binding. However, its absence from the β subunit in mutant and artificial Hbs simply enhances the O₂ affinity of the α subunit. Accordingly, the inter-subunit interactions between α and β subunits are nonsymmetric but substantial for HbA to perform cooperative O₂ binding.

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.

Structures of oxygen dissociation intermediates of 400 kDa V2 hemoglobin provide coarse snapshots of the protein allostery

Ever since the historic discovery of the cooperative oxygenation of its multiple subunits, hemoglobin (Hb) has been among the most exhaustively studied allosteric proteins. However, the lack of structural information on the intermediates between oxygenated and deoxygenated forms prevents our detailed understanding of the molecular mechanism of its allostery. It has been difficult to prepare crystals of intact oxy-deoxy intermediates and to individually identify the oxygen saturation for each subunit. However, our recent crystallographic studies have demonstrated that giant Hbs from annelids are suitable for overcoming these problems and can provide abundant information on oxy-deoxy intermediate structures. Here, we report the crystal structures of oxy-deoxy intermediates of a 400 kDa Hb (V2Hb) from the annelid Lamellibrachia satsuma, following up on a series of previous studies of similar giant Hbs. Four intermediate structures had average oxygen saturations of 78%, 69%, 55%, and 26%, as determined by the occupancy refinement of the bound oxygen based on ambient temperature factors. The structures demonstrate that the cooperative oxygen dissociation is weaker, large ternary and quaternary changes are induced at a later stage of the oxygen dissociation process, and the ternary and quaternary changes are smaller with local perturbations. Nonetheless, the overall structural transition seemed to proceed in the manner of the MWC two-state model. Our crystallographic snapshots of the allosteric transition of V2Hb provide important experimental evidence for a more detailed understanding of the allostery of Hbs by extension of the Monod–Wyman–Changeux (MWC) model.

From hemoglobin allostery to hemoglobin-based oxygen carriers

Hemoglobin (Hb) plays its vital role through structural and functional properties evolutionarily optimized to work within red blood cells, i.e., the tetrameric assembly, well-defined oxygen affinity, positive cooperativity, and heterotropic allosteric regulation by protons, chloride and 2,3-diphosphoglycerate. Outside red blood cells, the Hb tetramer dissociates into dimers, which exhibit high oxygen affinity and neither cooperativity nor allosteric regulation. They are prone to extravasate, thus scavenging endothelial NO and causing hypertension, and cause nephrotoxicity. In addition, they are more prone to autoxidation, generating radicals. The need to overcome the adverse effects associated with cell-free Hb has always been a major hurdle in the development of substitutes of allogeneic blood transfusions for all clinical situations where blood is unavailable or cannot be used due to, for example, religious objections. This class of therapeutics, indicated as hemoglobin-based oxygen carriers (HBOCs), is formed by genetically and/or chemically modified Hbs. Many efforts were devoted to the exploitation of the wealth of biochemical and biophysical information available on Hb structure, function, and dynamics to design safe HBOCs, overcoming the negative effects of free plasma Hb. Unfortunately, so far, no HBOC has been approved by FDA and EMA, except for compassionate use. However, the unmet clinical needs that triggered intensive investigations more than fifty years ago are still awaiting an answer. Recently, HBOCs "repositioning" has led to their successful application in organ perfusion fluids.

Hemoglobin allostery and pharmacology

The oxygen demands of the human body require the constant circulation of blood carrying an enormous concentration of hemoglobin (Hb). Oxygen transport depends not only on the amount of Hb, but also on the control over the affinity of the protein for the gas, which can be optimized for the environmental conditions by changes in the concentration of effectors (hydrogen ions, chloride, CO2, and DPG) inside the red cell. Some pathological conditions affecting Hb may benefit from pharmacological interventions to increase or decrease its affinity for oxygen, or otherwise modify its properties, or alter its biosynthesis. Examples of such conditions include sickle cell anemia, thalassemias and inherited hemoglobinopathies. Effective and safe drugs such as voxelotor, bezafibrate and efaproxiral are available that significantly increase or decrease Hb oxygen affinity. Some medical conditions not directly affecting the blood or its oxygen carrying capacity may also be relieved by the manipulation of Hb. For example, the standard treatment of acute cyanide poisoning requires the oxidation of a fraction of the Hb in the bloodstream so that it efficiently scavenges cyanide. Tumors are often extremely hypoxic and therefore strongly resistant to radiotherapy; the sensitivity of cancerous tissue to X-rays may be increased by improved oxygenation through drugs binding Hb. This review attempts to provide a systematic exploration of the pharmacology of Hb, its molecular basis, and its intended and possible uses.

Surface-anchored framework for generating RhD-epitope stealth red blood cells

Rhesus D (RhD) is one of the most important immunogenic antigens on red blood cells (RBCs). However, the supply of RhD-negative blood frequently faces critical shortages in clinical practice, and the positive-to-negative transition of the RhD antigen remains a great challenge. Here, we developed an alternative approach for sheltering the epitopes on RhD-positive RBCs using a surface-anchored framework, which is flexible but can achieve an optimal shield effect with minimal physicochemical influence on the cell. The chemical framework completely obstructed the RhD antigens on the cell surface, and the assessments of both blood transfusion in a mouse model and immunostimulation with human RhD-positive RBCs in a rabbit model confirmed the RhD-epitope stealth characteristics of the engineered RBCs. This work provides an efficient methodology for improving the cell surface for universal blood transfusion and generally indicates the potential of rationally designed cell surface engineering for transfusion and transplantation medicine.

Hemoglobin: Structure, Function and Allostery

This chapter reviews how allosteric (heterotrophic) effectors and natural mutations impact hemoglobin (Hb) primary physiological function of oxygen binding and transport. First, an introduction about the structure of Hb is provided, including the ensemble of tense and relaxed Hb states and the dynamic equilibrium of Hb multistate. This is followed by a brief review of Hb variants with altered Hb structure and oxygen binding properties. Finally, a review of different endogenous and exogenous allosteric effectors of Hb is presented with particular emphasis on the atomic interactions of synthetic ligands with altered allosteric function of Hb that could potentially be harnessed for the treatment of diseases.

Evolutionary and functional insights into the mechanism underlying body-size-related adaptation of mammalian hemoglobin

Hemoglobin (Hb) represents a model protein to study molecular adaptation in vertebrates. Although both affinity and cooperativity of oxygen binding to Hb affect tissue oxygen delivery, only the former was thought to determine molecular adaptations of Hb. Here, we suggest that Hb affinity and cooperativity reflect evolutionary and physiological adaptions that optimized tissue oxygen delivery. To test this hypothesis, we derived the relationship between the Hill coefficient and the relative affinity and conformational changes parameters of the Monod-Wymann-Changeux allosteric model and graphed the 'biophysical Hill landscape' describing this relation. We found that mammalian Hb cooperativity values all reside on a ridge of maximum cooperativity along this landscape that allows for both gross- and fine-tuning of tissue oxygen unloading to meet the distinct metabolic requirements of mammalian tissues for oxygen. Our findings reveal the mechanism underlying body size-related adaptation of mammalian Hb. The generality and implications of our findings are discussed.

Phthalides serve as potent modulators to boost fetal hemoglobin induction therapy for β-hemoglobinopathies

Fetal hemoglobin (HbF) induction therapy has become the most promising strategy for treating β-hemoglobinopathies, including sickle-cell diseases and β-thalassemia. However, subtle but critical structural difference exists between HbF and normal adult hemoglobin (HbA), which inevitably leads to reduced binding of the endogenous modulator 2,3-bisphosphoglycerate (2,3-BPG) to HbF and thus increased oxygen affinity and decreased oxygen transport efficiency of HbF. We combined the oxygen equilibrium experiments, resonance Raman (RR) spectroscopy, and molecular docking modeling, and we discuss 2 phthalides, z-butylidenephthalide and z-ligustilide, that can effectively lower the oxygen affinity of HbF. They adjust it to a level closer to that of HbA and make it a more satisfactory oxygen carrier for adults. From the oxygen equilibrium curve measurements, we show that the 2 phthalides are more effective than 2,3-BPG for modulating HbF. The RR spectra show that phthalides allosterically stabilize the oxygenated HbF in the low oxygen affinity conformation, and the molecular docking modeling reveals that the 2 chosen phthalides interact with HbF via the cleft around the γ1/γ2 interface with a binding strength ∼1.6 times stronger than that of 2,3-BPG. We discuss the implications of z-butylidenephthalide and z-ligustilide in boosting the efficacy of HbF induction therapy to mitigate the clinical severities of β-hemoglobinopathies.

Dynamics of Quaternary Structure Transitions in R-State Carbonmonoxyhemoglobin Unveiled in Time-Resolved X-ray Scattering Patterns Following a Temperature Jump

It is well-known that tetrameric hemoglobin binds ligands cooperatively by undergoing a ligand-induced T → R quaternary structure transition, a structure-function relationship that has long served as a model system for understanding allostery in proteins. However, kinetic studies of the reverse, R → T quaternary structure transition following photolysis of carbonmonoxyhemoglobin (HbCO) reveal complex behavior that may be better explained by the presence of two different R quaternary structures coexisting in thermal equilibrium. Indeed, we report here time-resolved small- and wide-angle X-ray scattering (SAXS/WAXS) patterns of HbCO following a temperature jump that not only provide unambiguous evidence for more than one R state, but also unveil the time scale for interconversion between them. Since the time scale for the photolysis-induced R → T transition is likely different for different R-states, this structural heterogeneity must be accounted for to properly explain the kinetic heterogeneity observed in time-resolved spectroscopic studies following photolysis of HbCO.

Design, Synthesis, and Biological Evaluation of Allosteric Effectors That Enhance CO Release from Carboxyhemoglobin

Carbon monoxide (CO) poisoning causes between 5,000-6,000 deaths per year in the US alone. The development of small molecule allosteric effectors of CO binding to hemoglobin (Hb) represents an important step toward making effective therapies for CO poisoning. To that end, we have found that the synthetic peptide IRL 2500 enhances CO release from COHb in air, but with concomitant hemolytic activity. We describe herein the design, synthesis, and biological evaluation of analogs of IRL 2500 that enhance the release of CO from COHb without hemolysis. These novel structures show improved aqueous solubility and reduced hemolytic activity and could lead the way to the development of small molecule therapeutics for the treatment of CO poisoning.

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.

Engineered synthetic antibodies as probes to quantify the energetic contributions of ligand binding to conformational changes in proteins

Conformational changes in proteins due to ligand binding are ubiquitous in biological processes and are integral to many biological systems. However, it is often challenging to link ligand-induced conformational changes to a resulting biological function because it is difficult to distinguish between the energetic components associated with ligand binding and those due to structural rearrangements. Here, we used a unique approach exploiting conformation-specific and regio-specific synthetic antibodies (sABs) to probe the energetic contributions of ligand binding to conformation changes. Using maltose-binding protein (MBP) as a model system, customized phage-display selections were performed to generate sABs that stabilize MBP in different conformational states, modulating ligand-binding affinity in competitive, allosteric, or peristeric manners. We determined that the binding of a closed conformation-specific sAB (sAB-11M) to MBP in the absence of maltose is entropically driven, providing new insight into designing antibody-stabilized protein interactions. Crystal structures of sABs bound to MBP, together with biophysical data, delineate the basis of free energy differences between different conformational states and confirm the use of the sABs as energy probes for dissecting enthalpic and entropic contributions to conformational transitions. Our work provides a foundation for investigating the energetic contributions of distinct conformational dynamics to specific biological outputs. We anticipate that our approach also may be valuable for analyzing the energy landscapes of regulatory proteins controlling biological responses to environmental changes.

A role of heme side-chains of human hemoglobin in its function revealed by circular dichroism and resonance Raman spectroscopy

Structural changes of heme side-chains of human adult hemoglobin (Hb A) upon ligand (O₂ or CO) dissociation have been studied by circular dichroism (CD) and resonance Raman (RR) spectroscopies. We point out the occurrence of appreciable deformation of heme side-chains like vinyl and propionate groups prior to the out-of-plane displacement of heme iron. Referring to the recent fine resolved crystal structure of Hb A, the deformations of heme side-chains take place only in the β subunits. However, these changes are not observed in the isolated β chain (β₄ homotetramer) and, therefore, are associated with the α-β inter-subunit interactions. For the communications between α and β subunits in Hb A regarding signals of ligand dissociation, possible routes are proposed on the basis of the time-resolved absorption, CD, MCD (magnetic CD), and RR spectroscopies. Our finding of the movements of heme side-chains would serve as one of the clues to solve the cooperative O₂ binding mechanism of Hb A.

Structure and function of haemoglobins

Haemoglobin (Hb) is widely known as the iron-containing protein in blood that is essential for O₂ transport in mammals. Less widely recognised is that erythrocyte Hb belongs to a large family of Hb proteins with members distributed across all three domains of life-bacteria, archaea and eukaryotes. This review, aimed chiefly at researchers new to the field, attempts a broad overview of the diversity, and common features, in Hb structure and function. Topics include structural and functional classification of Hbs; principles of O₂ binding affinity and selectivity between O₂/NO/CO and other small ligands; hexacoordinate (containing bis-imidazole coordinated haem) Hbs; bacterial truncated Hbs; flavohaemoglobins; enzymatic reactions of Hbs with bioactive gases, particularly NO, and protection from nitrosative stress; and, sensor Hbs. A final section sketches the evolution of work on the structural basis for allosteric O₂ binding by mammalian RBC Hb, including the development of newer kinetic models. Where possible, reference to historical works is included, in order to provide context for current advances in Hb research.

Formulation, construction and analysis of kinetic models of metabolism: A review of modelling frameworks

Kinetic models are critical to predict the dynamic behaviour of metabolic networks. Mechanistic kinetic models for large networks remain uncommon due to the difficulty of fitting their parameters. Recent modelling frameworks promise new ways to overcome this obstacle while retaining predictive capabilities. In this review, we present an overview of the relevant mathematical frameworks for kinetic formulation, construction and analysis. Starting with kinetic formalisms, we next review statistical methods for parameter inference, as well as recent computational frameworks applied to the construction and analysis of kinetic models. Finally, we discuss opportunities and limitations hindering the development of larger kinetic reconstructions.

Using the MWC model to describe heterotropic interactions in hemoglobin

Hemoglobin is a classical model allosteric protein. Research on hemoglobin parallels the development of key cooperativity and allostery concepts, such as the 'all-or-none' Hill formalism, the stepwise Adair binding formulation and the concerted Monod-Wymann-Changuex (MWC) allosteric model. While it is clear that the MWC model adequately describes the cooperative binding of oxygen to hemoglobin, rationalizing the effects of H+, CO2 or organophosphate ligands on hemoglobin-oxygen saturation using the same model remains controversial. According to the MWC model, allosteric ligands exert their effect on protein function by modulating the quaternary conformational transition of the protein. However, data fitting analysis of hemoglobin oxygen saturation curves in the presence or absence of inhibitory ligands persistently revealed effects on both relative oxygen affinity (c) and conformational changes (L), elementary MWC parameters. The recent realization that data fitting analysis using the traditional MWC model equation may not provide reliable estimates for L and c thus calls for a re-examination of previous data using alternative fitting strategies. In the current manuscript, we present two simple strategies for obtaining reliable estimates for MWC mechanistic parameters of hemoglobin steady-state saturation curves in cases of both evolutionary and physiological variations. Our results suggest that the simple MWC model provides a reasonable description that can also account for heterotropic interactions in hemoglobin. The results, moreover, offer a general roadmap for successful data fitting analysis using the MWC model.

Phthalide Derivatives from Angelica Sinensis Decrease Hemoglobin Oxygen Affinity: A New Allosteric-Modulating Mechanism and Potential Use as 2,3-BPG Functional Substitutes

Angelica sinensis (AS), one of the most versatile herbal medicines remains widely used due to its multi-faceted pharmacologic activities. Besides its traditional use as the blood-nourishing tonic, its anti-hypertensive, anti-cardiovascular, neuroprotective and anti-cancer effects have been reported. Albeit the significant therapeutic effects, how AS exerts such diverse efficacies from the molecular level remains elusive. Here we investigate the influences of AS and four representative phthalide derivatives from AS on the structure and function of hemoglobin (Hb). From the spectroscopy and oxygen equilibrium experiments, we show that AS and the chosen phthalides inhibited the oxygenated Hb from transforming into the high-affinity “relaxed” (R) state, decreasing Hb’s oxygen affinity. It reveals that phthalides cooperate with the endogenous Hb modulator, 2,3-bisphosphoglycerate (2,3-BPG) to synergetically regulate Hb allostery. From the docking modeling, phthalides appear to interact with Hb mainly through its α₁/α₂ interface, likely strengthening four (out of six) Hb “tense” (T) state stabilizing salt-bridges. A new allosteric-modulating mechanism is proposed to rationalize the capacity of phthalides to facilitate Hb oxygen transport, which may be inherently correlated with the therapeutic activities of AS. The potential of phthalides to serve as 2,3-BPG substitutes/supplements and their implications in the systemic biology and preventive medicine are discussed.

Interrelationship among Fe-His Bond Strengths, Oxygen Affinities, and Intersubunit Hydrogen Bonding Changes upon Ligand Binding in the β Subunit of Human Hemoglobin: The Alkaline Bohr Effect

Regulation of the oxygen affinity of human adult hemoglobin (Hb A) at high pH, known as the alkaline Bohr effect, is essential for its physiological function. In this study, structural mechanisms of the alkaline Bohr effect and pH-dependent O₂ affinity changes were investigated via ¹H nuclear magnetic resonance and visible and UV resonance Raman spectra of mutant Hbs, Hb M Iwate (αH87Y) and Hb M Boston (αH58Y). It was found that even though the binding of O₂ to the α subunits is forbidden in the mutant Hbs, the O₂ affinity was higher at alkaline pH than at neutral pH, and concomitantly, the Fe-His stretching frequency of the β subunits was shifted to higher values. Thus, it was confirmed for the β subunits that the stronger the Fe-His bond, the higher the O₂ affinity. It was found in this study that the quaternary structure of α(Fe³⁺)β(Fe²⁺-CO) of the mutant Hb is closer to T than to the ordinary R at neutral pH. The retained Aspβ94-Hisβ146 hydrogen bond makes the extent of proton release smaller upon ligand binding from Hisβ146, known as one of residues contributing to the alkaline Bohr effect. For these T structures, the Aspα94-Trpβ37 hydrogen bond in the hinge region and the Tyrα42-Aspβ99 hydrogen bond in the switch region of the α₁-β₂ interface are maintained but elongated at alkaline pH. Thus, a decrease in tension in the Fe-His bond of the β subunits at alkaline pH causes a substantial increase in the change in global structure upon binding of CO to the β subunit.

Modulation of hemoglobin dynamics by an allosteric effector

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

Experiments on Hemoglobin in Single Crystals and Silica Gels Distinguish among Allosteric Models

Trapping quaternary structures of hemoglobin in single crystals or by encapsulation in silica gels has provided a demanding set of data to test statistical mechanical models of allostery. In this work, we compare the results of those experiments with predictions of the four major allosteric models for hemoglobin: the quaternary two-state model of Monod, Wyman, and Changeux; the tertiary two-state model of Henry et al., which is the simplest extension of the Monod-Wyman-Changeux model to include pre-equilibria of tertiary as well as quaternary conformations; the structure-based model of Szabo and Karplus; and the modification of the latter model by Lee and Karplus. We show that only the tertiary two-state model can provide a near quantitative explanation of the single-crystal and gel experimental results.

From protein structure to function via single crystal optical spectroscopy

The more than 100,000 protein structures determined by X-ray crystallography provide a wealth of information for the characterization of biological processes at the molecular level. However, several crystallographic “artifacts,” including conformational selection, crystallization conditions and radiation damages, may affect the quality and the interpretation of the electron density maps, thus limiting the relevance of structure determinations. Moreover, for most of these structures, no functional data have been obtained in the crystalline state, thus posing serious questions on their validity in infereing protein mechanisms. In order to solve these issues, spectroscopic methods have been applied for the determination of equilibrium and kinetic properties of proteins in the crystalline state. These methods are UV-vis spectrophotometry, spectrofluorimetry, IR, EPR, Raman, and resonance Raman spectroscopy. Some of these approaches have been implemented with on-line instruments at X-ray synchrotron beamlines. Here, we provide an overview of investigations predominantly carried out in our laboratory by single crystal polarized absorption UV-vis microspectrophotometry, the most applied technique for the functional characterization of proteins in the crystalline state. Studies on hemoglobins, pyridoxal 5′-phosphate dependent enzymes and green fluorescent protein in the crystalline state have addressed key biological issues, leading to either straightforward structure-function correlations or limitations to structure-based mechanisms.

Identification of a small molecule that increases hemoglobin oxygen affinity and reduces SS erythrocyte sickling

Small molecules that increase the oxygen affinity of human hemoglobin may reduce sickling of red blood cells in patients with sickle cell disease. We screened 38 700 compounds using small molecule microarrays and identified 427 molecules that bind to hemoglobin. We developed a high-throughput assay for evaluating the ability of the 427 small molecules to modulate the oxygen affinity of hemoglobin. We identified a novel allosteric effector of hemoglobin, di(5-(2,3-dihydro-1,4-benzodioxin-2-yl)-4H-1,2,4-triazol-3-yl)disulfide (TD-1). TD-1 induced a greater increase in oxygen affinity of human hemoglobin in solution and in red blood cells than did 5-hydroxymethyl-2-furfural (5-HMF), N-ethylmaleimide (NEM), or diformamidine disulfide. The three-dimensional structure of hemoglobin complexed with TD-1 revealed that monomeric units of TD-1 bound covalently to β-Cys93 and β-Cys112, as well as noncovalently to the central water cavity of the hemoglobin tetramer. The binding of TD-1 to hemoglobin stabilized the relaxed state (R3-state) of hemoglobin. TD-1 increased the oxygen affinity of sickle hemoglobin and inhibited in vitro hypoxia-induced sickling of red blood cells in patients with sickle cell disease without causing hemolysis. Our study indicates that TD-1 represents a novel lead molecule for the treatment of patients with sickle cell disease.

Experimental basis for a new allosteric model for multisubunit proteins

Monod, Wyman, and Changeux (MWC) explained allostery in multisubunit proteins with a widely applied theoretical model in which binding of small molecules, so-called allosteric effectors, affects reactivity by altering the equilibrium between more reactive (R) and less reactive (T) quaternary structures. In their model, each quaternary structure has a single reactivity. Here, we use silica gels to trap protein conformations and a new kind of laser photolysis experiment to show that hemoglobin, the paradigm of allostery, exhibits two ligand binding phases with the same fast and slow rates in both R and T quaternary structures. Allosteric effectors change the fraction of each phase but not the rates. These surprising results are readily explained by the simplest possible extension of the MWC model to include a preequilibrium between two tertiary conformations that have the same functional properties within each quaternary structure. They also have important implications for the long-standing question of a structural explanation for the difference in hemoglobin oxygen affinity of the two quaternary structures.

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

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

Reverse engineering the cooperative machinery of human hemoglobin

Hemoglobin transports molecular oxygen from the lungs to all human tissues for cellular respiration. Its α₂β₂ tetrameric assembly undergoes cooperative binding and releasing of oxygen for superior efficiency and responsiveness. Over past decades, hundreds of hemoglobin structures were determined under a wide range of conditions for investigation of molecular mechanism of cooperativity. Based on a joint analysis of hemoglobin structures in the Protein Data Bank (Ren, companion article), here I present a reverse engineering approach to elucidate how two subunits within each dimer reciprocate identical motions that achieves intradimer cooperativity, how ligand-induced structural signals from two subunits are integrated to drive quaternary rotation, and how the structural environment at the oxygen binding sites alter their binding affinity. This mechanical model reveals the intricate design that achieves the cooperative mechanism and has previously been masked by inconsistent structural fluctuations. A number of competing theories on hemoglobin cooperativity and broader protein allostery are reconciled and unified.

Structure of fully liganded Hb ζ2β2s trapped in a tense conformation

A variant Hb ζ2β2(s) that is formed from sickle hemoglobin (Hb S; α2β2(s)) by exchanging adult α-globin with embryonic ζ-globin subunits shows promise as a therapeutic agent for sickle-cell disease (SCD). Hb ζ2β2(s) inhibits the polymerization of deoxygenated Hb S in vitro and reverses characteristic features of SCD in vivo in mouse models of the disorder. When compared with either Hb S or with normal human adult Hb A (α2β2), Hb ζ2β2(s) exhibits atypical properties that include a high oxygen affinity, reduced cooperativity, a weak Bohr effect and blunted 2,3-diphosphoglycerate allostery. Here, the 1.95 Å resolution crystal structure of human Hb ζ2β2(s) that was expressed in complex transgenic knockout mice and purified from their erythrocytes is presented. When fully liganded with carbon monoxide, Hb ζ2β2(s) displays a central water cavity, a ζ1-β(s)2 (or ζ2-β(s)1) interface, intersubunit salt-bridge/hydrogen-bond interactions, C-terminal βHis146 salt-bridge interactions, and a β-cleft, that are highly unusual for a relaxed hemoglobin structure and are more typical of a tense conformation. These quaternary tense-like features contrast with the tertiary relaxed-like conformations of the ζ1β(s)1 dimer and the CD and FG corners, as well as the overall structures of the heme cavities. This crystallographic study provides insights into the altered oxygen-transport properties of Hb ζ2β2(s) and, moreover, decouples tertiary- and quaternary-structural events that are critical to Hb ligand binding and allosteric function.

Solution structure and dynamics of human hemoglobin in the carbonmonoxy form

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

Molecular controls of the oxygenation and redox reactions of hemoglobin

SIGNIFICANCE

The broad classes of O(2)-binding proteins known as hemoglobins (Hbs) carry out oxygenation and redox functions that allow organisms with significantly different physiological demands to exist in a wide range of environments. This is aided by allosteric controls that modulate the protein's redox reactions as well as its O(2)-binding functions.

RECENT ADVANCES

The controls of Hb's redox reactions can differ appreciably from the molecular controls for Hb oxygenation and come into play in elegant mechanisms for dealing with nitrosative stress, in the malarial resistance conferred by sickle cell Hb, and in the as-yet unsuccessful designs for safe and effective blood substitutes.

CRITICAL ISSUES

An important basic principle in consideration of Hb's redox reactions is the distinction between kinetic and thermodynamic reaction control. Clarification of these modes of control is critical to gaining an increased understanding of Hb-mediated oxidative processes and oxidative toxicity in vivo.

FUTURE DIRECTIONS

This review addresses emerging concepts and some unresolved questions regarding the interplay between the oxygenation and oxidation reactions of structurally diverse Hbs, both within red blood cells and under acellular conditions. Developing methods that control Hb-mediated oxidative toxicity will be critical to the future development of Hb-based blood substitutes.

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.

Tertiary and quaternary effects in the allosteric regulation of animal hemoglobins

In the last decade, protein allostery has experienced a major resurgence, boosted by the extension of the concept to systems of increasing complexity and by its exploitation for the development of drugs. Expansion of the field into new directions has not diminished the key role of hemoglobin as a test molecule for theory and experimental validation of allosteric models. Indeed, the diffusion of hemoglobins in all kingdoms of life and the variety of functions and of quaternary assemblies based on a common tertiary fold indicate that this superfamily of proteins is ideally suited for investigating the physical and molecular basis of allostery and firmly maintains its role as a main player in the field. This review is an attempt to briefly recollect common and different strategies adopted by metazoan hemoglobins, from monomeric molecules to giant complexes, exploiting homotropic and heterotropic allostery to increase their functional dynamic range. 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.