Direct determination of protonation states of histidine residues in a 2 A neutron structure of deoxy-human normal adult hemoglobin and implications for the Bohr effect

Histidine in Proteins: pH-Dependent Interplay between π-π, Cation-π, and CH-π Interactions

Histidine (His) stands out as the most versatile natural amino acid due to its side chain's facile propensity to protonate at physiological pH, leading to a transition from aromatic to cationic characteristics and thereby enabling diverse biomolecular interactions. In this study, our objective was to quantify the energetics and geometries of pairwise interactions involving His at varying pH levels. Through quantum chemical calculations, we discovered that His exhibits robust participation in both π-π and cation-π interactions, underscoring its ability to adopt a π or cationic nature, akin to other common residues. Of particular note, we found that the affinity of protonated His for aromatic residues (via cation-π interactions) is greater than the affinity of neutral His for either cationic residues (also via cation-π interactions) or aromatic residues (via π-π interactions). Furthermore, His frequently engages in CH-π interactions, and notably, depending on its protonation state, we found that some instances of hydrogen bonding by His exhibit greater stability than is typical for interamino acid hydrogen bonds. The strength of the pH-dependent pairwise energies of His with aromatic residues is supported by the abundance of pairwise interactions with His of low and high predicted pKa values. Overall, our findings illustrate the contribution of His interactions to protein stability and its potential involvement in conformational changes despite its relatively low abundance in proteins.

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.

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.

Determination of Histidine Protonation States in Proteins by Fast Magic Angle Spinning NMR

Histidine residues play important structural and functional roles in proteins, such as serving as metal-binding ligands, mediating enzyme catalysis, and modulating proton channel activity. Many of these activities are modulated by the ionization state of the imidazole ring. Here we present a fast MAS NMR approach for the determination of protonation and tautomeric states of His at frequencies of 40-62 kHz. The experiments combine 1H detection with selective magnetization inversion techniques and transferred echo double resonance (TEDOR)-based filters, in 2D heteronuclear correlation experiments. We illustrate this approach using microcrystalline assemblies of HIV-1 CACTD-SP1 protein.

In situ remediation of 2,4-dicholrophenoxyacetic acid herbicide using amine-functionalized imidazole coordination complexes

Demand of clean water is always a major concern due to continuous use of toxic pesticides and herbicides to overcome food scarcity. In Asian countries, wide use of ionizable 2,4-D herbicide has worsen problem due to its less binding ability with soil and can easily contaminate drinking water that causes potential risks to humans and environment. The present research focused on synthesis of amino-factionalized coordination complexes using imidazole-based amino benzoic acid ligands for remediation of acidic 2,4-D herbicide. Coordination complexes characterized with FTIR, ¹H-NMR, elemental analysis, thermogravimetric analysis, powder XRD, and BET revealed successful incorporation of functionalized groups with high thermal stability and surface area that make them suitable for adsorption experiments. Batch adsorption experiments conducted at different temperature conditions depicted the spontaneous physisorption process (- ∆G) having endothermic nature (∆H, ∆S). The removal efficiency of the amino-functionalized coordination complex is found to be higher (73%) compared to non-functionalized (35%) and acetic anhydride-functionalized coordination complex (42%). Kinetic studies supported pseudo 2nd-order kinetics with three phases of adsorption depicted by intra-particle diffusion model. Amino-functionalized complexes favored Langmuir isotherm while Freundlich isotherm is best fitted for non-functionalized complexes. The synthesized adsorbents were also proven to be effective for removal of herbicide from irrigated wastewater with average percent removal of 56% for amino functionalized, acetic anhydride functionalized (23%), and non-functionalized (20%).

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M^pro, a cysteine protease, have been determined, facilitating structure-based drug design. M^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41–Cys145, M^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV M^pro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N_δ (HD) and N_ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics.

Room-temperature neutron and X-ray data collection of 3CL Mpro from SARS-CoV-2

The replication of SARS-CoV-2 produces two large polyproteins, pp1a and pp1ab, that are inactive until cleavage by the viral chymotrypsin-like cysteine protease enzyme (3CL Mpro) into a series of smaller functional proteins. At the heart of 3CL Mpro is an unusual catalytic dyad formed by the side chains of His41 and Cys145 and a coordinated water molecule. The catalytic mechanism by which the enzyme operates is still unknown, as crucial information on the protonation states within the active site is unclear. To experimentally determine the protonation states of the catalytic site and of the other residues in the substrate-binding cavity, and to visualize the hydrogen-bonding networks throughout the enzyme, room-temperature neutron and X-ray data were collected from a large H/D-exchanged crystal of ligand-free (apo) 3CL Mpro.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M ^pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M ^pro , a cysteine protease, have been determined, facilitating structure-based drug design. M ^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, M ^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV M ^pro , but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M ^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M ^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α -ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N _δ (HD) and N _ϵ (HE) protonation of His41 and His164, respectively, the α -ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M ^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

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.

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

A recent molecular dynamics investigation into the stability of hemoglobin concluded that the unliganded protein is only stable in the T state when a solvent box is used in the simulations that is ten times larger than what is usually employed (El Hage et al., 2018). Here, we express three main concerns about that study. In addition, we find that with an order of magnitude more statistics, the reported box size dependence is not reproducible. Overall, no significant effects on the kinetics or thermodynamics of conformational transitions were observed.

What Are We Missing by Not Measuring the Net Charge of Proteins?

The net electrostatic charge (Z) of a folded protein in solution represents a bird's eye view of its surface potentials-including contributions from tightly bound metal, solvent, buffer, and cosolvent ions-and remains one of its most enigmatic properties. Few tools are available to the average biochemist to rapidly and accurately measure Z at pH≠pI. Tools that have been developed more recently seem to go unnoticed. Most scientists are content with this void and estimate the net charge of a protein from its amino acid sequence, using textbook values of pK_a . Thus, Z remains unmeasured for nearly all folded proteins at pH≠pI. When marveling at all that has been learned from accurately measuring the other fundamental property of a protein-its mass-one wonders: what are we missing by not measuring the net charge of folded, solvated proteins? A few big questions immediately emerge in bioinorganic chemistry. When a single electron is transferred to a metalloprotein, does the net charge of the protein change by approximately one elementary unit of charge or does charge regulation dominate, that is, do the pK_a values of most ionizable residues (or just a few residues) adjust in response to (or in concert with) electron transfer? Would the free energy of charge regulation (ΔΔG_z ) account for most of the outer sphere reorganization energy associated with electron transfer? Or would ΔΔG_z contribute more to the redox potential? And what about metal binding itself? When an apo-metalloprotein, bearing minimal net negative charge (e.g., Z=-2.0) binds one or more metal cations, is the net charge abolished or inverted to positive? Or do metalloproteins regulate net charge when coordinating metal ions? The author's group has recently dusted off a relatively obscure tool-the "protein charge ladder"-and used it to begin to answer these basic questions.

Neutron macromolecular crystallography

Neutron diffraction techniques permit direct determination of the hydrogen (H) and deuterium (D) positions in crystal structures of biological macromolecules at resolutions of ∼1.5 and 2.5 Å, respectively. In addition, neutron diffraction data can be collected from a single crystal at room temperature without radiation damage issues. By locating the positions of H/D-atoms, protonation states and water molecule orientations can be determined, leading to a more complete understanding of many biological processes and drug-binding. In the last ca. 5 years, new beamlines have come online at reactor neutron sources, such as BIODIFF at Heinz Maier-Leibnitz Zentrum and IMAGINE at Oak Ridge National Laboratory (ORNL), and at spallation neutron sources, such as MaNDi at ORNL and iBIX at the Japan Proton Accelerator Research Complex. In addition, significant improvements have been made to existing beamlines, such as LADI-III at the Institut Laue-Langevin. The new and improved instrumentations are allowing sub-mm3 crystals to be regularly used for data collection and permitting the study of larger systems (unit-cell edges >100 Å). Owing to this increase in capacity and capability, many more studies have been performed and for a wider range of macromolecules, including enzymes, signalling proteins, transport proteins, sugar-binding proteins, fluorescent proteins, hormones and oligonucleotides; of the 126 structures deposited in the Protein Data Bank, more than half have been released since 2013 (65/126, 52%). Although the overall number is still relatively small, there are a growing number of examples for which neutron macromolecular crystallography has provided the answers to questions that otherwise remained elusive.

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.

Molecular mechanisms of bio-catalysis of heme extraction from hemoglobin

Red blood cell hemolysis in sickle cell disease (SCD) releases free hemoglobin. Extracellular hemoglobin and its degradation products, free heme and iron, are highly toxic due to oxidative stress induction and decrease in nitric oxide availability. We propose an approach that helps to eliminate extracellular hemoglobin toxicity in SCD by employing a bacterial protein system that evolved to extract heme from extracellular hemoglobin. NEAr heme Transporter (NEAT) domains from iron-regulated surface determinant proteins from Staphylococcus aureus specifically bind free heme as well as facilitate its extraction from hemoglobin. We demonstrate that a purified NEAT domain fused with human haptoglobin β-chain is able to remove heme from hemoglobin and reduce heme content and peroxidase activity of hemoglobin. We further use molecular dynamics (MD) simulations to resolve molecular pathway of heme transfer from hemoglobin to NEAT, and to elucidate molecular mechanism of such heme transferring process. Our study is the first of its kind, in which simulations are employed to characterize the process of heme leaving hemoglobin and subsequent rebinding with a NEAT domain. Our MD results highlight important amino acid residues that facilitate heme transfer and will guide further studies for the selection of best NEAT candidate to attenuate free hemoglobin toxicity.

The use of neutron scattering to determine the functional structure of glycoside hydrolase

Neutron diffraction provides different information from X-ray diffraction, because neutrons are scattered by atomic nuclei, whereas X-rays are scattered by electrons. One of the key advantages of neutron crystallography is the ability to visualize hydrogen and deuterium atoms, making it possible to observe the protonation state of amino acid residues, hydrogen bonds, networks of water molecules and proton relay pathways in enzymes. But, because of technical difficulties, less than 100 enzyme structures have been evaluated by neutron crystallography to date. In this review, we discuss the advantages and disadvantages of neutron crystallography as a tool to investigate the functional structure of glycoside hydrolases, with some examples.

Visualizing the Bohr effect in hemoglobin: neutron structure of equine cyanomethemoglobin in the R state and comparison with human deoxyhemoglobin in the T state

Neutron crystallography provides direct visual evidence of the atomic positions of deuterium-exchanged H atoms, enabling the accurate determination of the protonation/deuteration state of hydrated biomolecules. Comparison of two neutron structures of hemoglobins, human deoxyhemoglobin (T state) and equine cyanomethemoglobin (R state), offers a direct observation of histidine residues that are likely to contribute to the Bohr effect. Previous studies have shown that the T-state N-terminal and C-terminal salt bridges appear to have a partial instead of a primary overall contribution. Four conserved histidine residues [αHis72(EF1), αHis103(G10), αHis89(FG1), αHis112(G19) and βHis97(FG4)] can become protonated/deuterated from the R to the T state, while two histidine residues [αHis20(B1) and βHis117(G19)] can lose a proton/deuteron. αHis103(G10), located in the α1:β1 dimer interface, appears to be a Bohr group that undergoes structural changes: in the R state it is singly protonated/deuterated and hydrogen-bonded through a water network to βAsn108(G10) and in the T state it is doubly protonated/deuterated with the network uncoupled. The very long-term H/D exchange of the amide protons identifies regions that are accessible to exchange as well as regions that are impermeable to exchange. The liganded relaxed state (R state) has comparable levels of exchange (17.1% non-exchanged) compared with the deoxy tense state (T state; 11.8% non-exchanged). Interestingly, the regions of non-exchanged protons shift from the tetramer interfaces in the T-state interface (α1:β2 and α2:β1) to the cores of the individual monomers and to the dimer interfaces (α1:β1 and α2:β2) in the R state. The comparison of regions of stability in the two states allows a visualization of the conservation of fold energy necessary for ligand binding and release.

O2 and Water Migration Pathways between the Solvent and Heme Pockets of Hemoglobin with Open and Closed Conformations of the Distal HisE7

Hemoglobin transports O2 by binding the gas at its four hemes. Hydrogen bonding between the distal histidine (HisE7) and heme-bound O2 significantly increases the affinity of human hemoglobin (HbA) for this ligand. HisE7 is also proposed to regulate the release of O2 to the solvent via a transient E7 channel. To reveal the O2 escape routes controlled by HisE7 and to evaluate its role in gating heme access, we compare simulations of O2 diffusion from the distal heme pockets of the T and R states of HbA performed with HisE7 in its open (protonated) and closed (neutral) conformations. Irrespective of HisE7's conformation, we observe the same four or five escape routes leading directly from the α- or β-distal heme pockets to the solvent. Only 21-53% of O2 escapes occur via these routes, with the remainder escaping through routes that encompass multiple internal cavities in HbA. The conformation of the distal HisE7 controls the escape of O2 from the heme by altering the distal pocket architecture in a pH-dependent manner, not by gating the E7 channel. Removal of the HisE7 side chain in the GlyE7 variant exposes the distal pockets to the solvent, and the percentage of O2 escapes to the solvent directly from the α- or β-distal pockets of the mutant increases to 70-88%. In contrast to O2, the dominant water route from the bulk solvent is gated by HisE7 because protonation and opening of this residue dramatically increase the rate of influx of water into the empty distal heme pockets. The occupancy of the distal heme site by a water molecule, which functions as an additional nonprotein barrier to binding of the ligand to the heme, is also controlled by HisE7. Overall, analysis of gas and water diffusion routes in the subunits of HbA and its GlyE7 variant sheds light on the contribution of distal HisE7 in controlling polar and nonpolar ligand movement between the solvent and the hemes.

Protein charge ladders reveal that the net charge of ALS-linked superoxide dismutase can be different in sign and magnitude from predicted values

This article utilized "protein charge ladders"-chemical derivatives of proteins with similar structure, but systematically altered net charge-to quantify how missense mutations that cause amyotrophic lateral sclerosis (ALS) affect the net negative charge (Z) of superoxide dismutase-1 (SOD1) as a function of subcellular pH and Zn(2+) stoichiometry. Capillary electrophoresis revealed that the net charge of ALS-variant SOD1 can be different in sign and in magnitude-by up to 7.4 units per dimer at lysosomal pH-than values predicted from standard pKa values of amino acids and formal oxidation states of metal ions. At pH 7.4, the G85R, D90A, and G93R substitutions diminished the net negative charge of dimeric SOD1 by up to +2.29 units more than predicted; E100K lowered net charge by less than predicted. The binding of a single Zn(2+) to mutant SOD1 lowered its net charge by an additional +2.33 ± 0.01 to +3.18 ± 0.02 units, however, each protein regulated net charge when binding a second, third, or fourth Zn(2+) (ΔZ < 0.44 ± 0.07 per additional Zn(2+) ). Both metalated and apo-SOD1 regulated net charge across subcellular pH, without inverting from negative to positive at the theoretical pI. Differential scanning calorimetry, hydrogen-deuterium exchange, and inductively coupled plasma mass spectrometry confirmed that the structure, stability, and metal content of mutant proteins were not significantly affected by lysine acetylation. Measured values of net charge should be used when correlating the biophysical properties of a specific ALS-variant SOD1 protein with its observed aggregation propensity or clinical phenotype.

Looking for Hydrogen Atoms: Neutron Crystallography Provides Novel Insights Into Protein Structure and Function

Neutron crystallography allows direct localization of hydrogen positions in biological macromolecules. Within enzymes, hydrogen atoms play a pivotal role in catalysis. Recent advances in instrumentation and sample preparation have helped to overcome the difficulties of performing neutron diffraction experiments on protein crystals. The application of neutron macromolecular crystallography to a growing number of proteins has yielded novel structural insights. The ability to accurately position water molecules, hydronium ions, and hydrogen atoms within protein structures has helped in the study of low-barrier hydrogen bonds and hydrogen-bonding networks. The determination of protonation states of protein side chains, substrates, and inhibitors in the context of the macromolecule has provided important insights into enzyme chemistry and ligand binding affinities, which can assist in the design of potent therapeutic agents. In this review, we give an overview of the method and highlight advances in knowledge attained through the application of neutron protein crystallography.

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.

Trapping of human hemoglobin by haptoglobin: molecular mechanisms and clinical applications

SIGNIFICANCE

Haptoglobin (Hp) is an abundant plasma protein controlling the fate of hemoglobin (Hb) released from red blood cells after intravascular hemolysis. The complex formed between Hp and Hb is extraordinary strong, and once formed, this protein-protein association can be considered irreversible.

RECENT ADVANCES

A model of the Hp-Hb complex has been generated and the first steps toward understanding the mechanism behind the shielding effects of Hp have been taken. The clinical potential of the complex for modulating inflammatory reactions and for functioning as an Hb-based oxygen carrier have been described.

CRITICAL ISSUES

The three-dimensional structure of the Hp-Hb complex is unknown. Moreover, Hp is not a homogeneous protein. There are two common alleles at the Hp genetic locus denoted Hp1 and Hp2, which when analyzed on the protein levels result in differences between their physiological behavior, particularly in their shielding against Hb-driven oxidative stress. Additional cysteine residues on the α-subunit allow Hp2 to form a variety of native multimers, which influence the biophysical and biological properties of Hp. The multimeric conformations, in turn, also modulate the glycosylation patterns of Hp by steric hindrance.

FUTURE DIRECTIONS

A detailed analysis of the influence of Hp glycosylation will be instrumental to generate a deeper understanding of its biological function. Several pathological conditions also modify the glycan compositions allowing Hp to be potentially used as a marker protein for these disorders.

Molecular modeling of the human hemoglobin-haptoglobin complex sheds light on the protective mechanisms of haptoglobin

Hemoglobin (Hb) plays a critical role in human physiological function by transporting O2. Hb is safe and inert within the confinement of the red blood cell but becomes reactive and toxic upon hemolysis. Haptoglobin (Hp) is an acute-phase serum protein that scavenges Hb and the resulting Hb-Hp complex is subjected to CD163-mediated endocytosis by macrophages. The interaction between Hb and Hp is extraordinarily strong and largely irreversible. As the structural details of the human Hb-Hp complex are not yet available, this study reports for the first time on insights of the binding modalities and molecular details of the human Hb-Hp interaction by means of protein-protein docking. Furthermore, residues that are pertinent for complex formation were identified by computational alanine scanning mutagenesis. Results revealed that the surface of the binding interface of Hb-Hp is not flat and protrudes into each binding partner. It was also observed that the secondary structures at the Hb-Hp interface are oriented as coils and α-helices. When dissecting the interface in more detail, it is obvious that several tyrosine residues of Hb, particularly β145Tyr, α42Tyr and α140Tyr, are buried in the complex and protected from further oxidative reactions. Such finding opens up new avenues for the design of Hp mimics which may be used as alternative clinical Hb scavengers.

Hemoglobin redux: combining neutron and X-ray diffraction with mass spectrometry to analyse the quaternary state of oxidized hemoglobins

Improvements in neutron diffraction instrumentation are affording the opportunity to re-examine the structures of vertebrate hemoglobins and to interrogate proton and solvent position changes between the different quaternary states of the protein. For hemoglobins of unknown primary sequence, structural studies of cyanomethemoglobin (CNmetHb) are being used to help to resolve sequence ambiguity in the mass spectra. These studies have also provided additional structural evidence for the involvement of oxidized hemoglobin in the process of erythrocyte senescence. X-ray crystal studies of Tibetan snow leopard CNmetHb have shown that this protein crystallizes in the B state, a structure with a more open dyad, which possibly has relevance to RBC band 3 protein binding and erythrocyte senescence. R-state equine CNmetHb crystal studies elaborate the solvent differences in the switch and hinge region compared with a human deoxyhemoglobin T-state neutron structure. Lastly, comparison of histidine protonation between the T and R state should enumerate the Bohr-effect protons.

Joint X-ray and neutron refinement with phenix.refine

Approximately 85% of the structures deposited in the Protein Data Bank have been solved using X-ray crystallography, making it the leading method for three-dimensional structure determination of macromolecules. One of the limitations of the method is that the typical data quality (resolution) does not allow the direct determination of H-atom positions. Most hydrogen positions can be inferred from the positions of other atoms and therefore can be readily included into the structure model as a priori knowledge. However, this may not be the case in biologically active sites of macromolecules, where the presence and position of hydrogen is crucial to the enzymatic mechanism. This makes the application of neutron crystallography in biology particularly important, as H atoms can be clearly located in experimental neutron scattering density maps. Without exception, when a neutron structure is determined the corresponding X-ray structure is also known, making it possible to derive the complete structure using both data sets. Here, the implementation of crystallographic structure-refinement procedures that include both X-ray and neutron data (separate or jointly) in the PHENIX system is described.

Protonation states of histidine and other key residues in deoxy normal human adult hemoglobin by neutron protein crystallography

The protonation states of the histidine residues key to the function of deoxy (T-state) human hemoglobin have been investigated using neutron protein crystallography. These residues can reversibly bind protons, thereby regulating the oxygen affinity of hemoglobin. By examining the OMIT F(o)-F(c) and 2F(o)-F(c) neutron scattering maps, the protonation states of 35 of the 38 His residues were directly determined. The remaining three residues were found to be disordered. Surprisingly, seven pairs of His residues from equivalent α or β chains, αHis20, αHis50, αHis58, αHis89, βHis63, βHis143 and βHis146, have different protonation states. The protonation of distal His residues in the α(1)β(1) heterodimer and the protonation of αHis103 in both subunits demonstrates that these residues may participate in buffering hydrogen ions and may influence the oxygen binding. The observed protonation states of His residues are compared with their ΔpK(a) between the deoxy and oxy states. Examination of inter-subunit interfaces provided evidence for interactions that are essential for the stability of the deoxy tertiary structure.