Molecular Mechanism of Myoglobin Autoxidation: Insights from Computer Simulations

First-principles study of molecular hydrogen binding to heme in competition with O2, NO and CO

Molecular hydrogen shows antioxidant activity and distinct efficacy towards vascular diseases, but the understanding of this is not yet satisfactory at the atomic level. In this work, we study the binding properties of H2 to the heme group in relation with other diatomic molecules (DMs), including O2, NO and CO, and their displacement reactions, using first-principles calculations. We carry out molecular modeling of the heme group, using iron-porphyrin with the imidazole ligand, i.e., FePIm, and smaller models of Fe(CnHn+2N2)2NH3 with n = 3 and 1, and of molecular complexes of heme-DM and -H. Through analysis of optimized geometries and energetics, it is found that the order of binding strength of DMs or H to the Fe of heme is NO > O2 > CO > H > H2 for FePIm-based systems, while it is H > O2 > NO > CO > H2 for model-based systems. We calculate the activation energies for displacement reactions of H2 and H by other DMs, revealing that the H2 displacements occur spontaneously while the H displacements require a large amount of energy. Finally, our calculations corroborate that the rate constants increase with increasing temperature according to the Arrhenius relation.

Protein dynamics govern the oxyferrous state lifetime of an artificial oxygen transport protein

It has long been known that the alteration of protein side chains that occlude or expose the heme cofactor to water can greatly affect the stability of the oxyferrous heme state. Here, we demonstrate that the rate of dynamically driven water penetration into the core of an artificial oxygen transport protein also correlates with oxyferrous state lifetime by reducing global dynamics, without altering the structure of the active site, via the simple linking of the two monomers in a homodimeric artificial oxygen transport protein using a glycine-rich loop. The tethering of these two helices does not significantly affect the active site structure, pentacoordinate heme-binding affinity, reduction potential, or gaseous ligand affinity. It does, however, significantly reduce the hydration of the protein core, as demonstrated by resonance Raman spectroscopy, backbone amide hydrogen exchange, and pKa shifts in buried histidine side chains. This further destabilizes the charge-buried entatic state and nearly triples the oxyferrous state lifetime. These data are the first direct evidence that dynamically driven water penetration is a rate-limiting step in the oxidation of these complexes. It furthermore demonstrates that structural rigidity that limits water penetration is a critical design feature in metalloenzyme construction and provides an explanation for both the failures and successes of earlier attempts to create oxygen-binding proteins.

First-principles study of mercaptoundecanoic acid molecule adsorption and gas molecule penetration onto silver surface: an insight for corrosion protection

Recently, 11-mercaptoundecanoic acid (MUA) molecule has attracted attention as a promising passivation agent of Ag nanowire (NW) network electrode for corrosion inhibition, but the underneath mechanism has not been elaborated. In this work, we investigate adsorption of MUA molecule on Ag(1 0 0) and Ag(1 1 1) surface, adsorption of air gas molecules of H2O, H2S and O2 on MUA molecular end surface, and their penetrations into the Ag surface using the first-principles calculations. Our calculations reveal that the MUA molecule is strongly bound to the Ag surface with the binding energies ranging from -0.47 to -2.06 eV and the Ag-S bond lengths of 2.68-2.97 Å by Lewis acid-base reaction. Furthermore, we find attractive interactions between the gas molecules and the MUA@Ag complexes upon their adsorptions and calculate activation barriers for their migrations from the outermost end of the complexes to the top of Ag surface. It is found that the penetrations of H2O and H2S are more difficult than the O2 penetration due to their higher activation barriers, while the O2 penetration is still difficult, confirming the corrosion protection of Ag NW network by adsorbing the uniform monolayer of MUA. With these findings, this work can contribute to finding a better passivation agent in the strategy of corrosion protection of Ag NW network electrode.

Antioxidation activity of molecular hydrogen via protoheme catalysis in vivo: an insight from ab initio calculations

Recently, molecular hydrogen has been found to exhibit antioxidation activity through many clinical experiments, but the mechanism has not been fully understandable at atomic level. In this work, we perform systematic ab initio calculations of protoheme-hydrogen complexes to clarify the antioxidation mechanism of molecular hydrogen. We make molecular modeling of iron-protoporphyrin coordinated by imidazole, FeP(Im), and its hydrogen as well as dihydrogen complexes, together with reactive oxygen/nitrogen species (RONS). We carry out structural optimization and Mulliken charge analysis, revealing the two kinds of bonding characteristics between FeP(Im) and H[Formula: see text]: dihydrogen bonding in the end-on asymmetric configuration and Kubas bonding in the side-on symmetric configuration of H[Formula: see text] molecule. The activation barriers for adsorption and dissociation of H[Formula: see text] on and further desorption of H atom from FeP(Im) are found to be below 2.78 eV at most, which is remarkably lower than the H-H bond breaking energy of 4.64 eV in free H[Formula: see text] molecule. We find that the hydrogen bond dissociation energies of FeP(Im)-H[Formula: see text] and -H complexes are lower than those of RONS-H complexes, indicating the decisive role of protoheme as an effective catalyst in RONS antioxidation by molecular hydrogen in vivo.

Heme Protein Binding of Sulfonamide Compounds: A Correlation Study by Spectroscopic, Calorimetric, and Computational Methods

Protein-ligand interaction studies are useful to determine the molecular mechanism of the binding phenomenon, leading to the establishment of the structure-function relationship. Here, we report the binding of well-known antibiotic sulfonamide drugs (sulfamethazine, SMZ; and sulfadiazine, SDZ) with heme protein myoglobin (Mb) using spectroscopic, calorimetric, ζ potential, and computational methods. Formation of a 1:1 complex between the ligand and Mb through well-defined equilibrium was observed. The binding constants obtained between Mb and SMZ/SDZ drugs were on the order of 10⁴ M^-1. SMZ with two additional methyl (-CH₃) substitutions has higher affinity than SDZ. Upon drug binding, a notable loss in the helicity (via circular dichroism) and perturbation of the three-dimensional (3D) protein structure (via infrared and synchronous fluorescence experiments) were observed. The binding also indicated the dominance of non-polyelectrolytic forces between the amino acid residues of the protein and the drugs. The ligand-protein binding distance signified high probability of energy transfer between them. Destabilization of the protein structure upon binding was evident from differential scanning calorimetry results and ζ potential analyses. Molecular docking presented the best probable binding sites of the drugs inside protein pockets. Thus, the present study explores the potential binding characteristics of two sulfonamide drugs (with different substitutions) with myoglobin, correlating the structural and energetic aspects.

An underlying softening mechanism in pale, soft and exudative - Like rabbit meat: The role of reactive oxygen species - Generating systems

This work investigated the role of reactive oxygen species (ROS) - generating systems on the softening of the pale, soft and exudative-like (PSE-like) rabbit meat during aging. PSE-like meat was induced by incubation of post-mortem rabbit Longissimus thoracis et lumborum at 37 °C for 3 h. During aging, PSE-like meat samples had higher values in peroxides value, thiobarbituric acid-reactive substances, metmyoglobin percentage, ferrylmyoglobin content, non-heme iron content, hydroxyl radical content and ROS concentration compared with the normal ones, suggesting that PSE-like incubation could activate lipid-oxidizing system, myoglobin-mediated oxidation system, together with metal-catalyzed oxidation system. Additionally, higher protein carbonyl content was observed in PSE-like meat, along with a significant loss in sulfhydryl group. The results of SDS-PAGE suggested that more serious protein degradation occurred in PSE-like meat. It is plausible that the activated ROS-generating system played an underlying role in the softening texture during the aging period of PSE-like meat.

Common and unique strategies of myoglobin evolution for deep-sea adaptation of diving mammals

Myoglobin (Mb) is highly concentrated in the myocytes of diving mammals such as whales and seals, in comparison with land animals, and its molecular evolution has played a crucial role in their deep-sea adaptation. We previously resurrected ancestral whale Mbs and demonstrated the evolutional strategies for higher solubility under macromolecular crowding conditions. Pinnipeds, such as seals and sea lions, are also expert diving mammals with Mb-rich muscles. In the present study, we resurrected ancestral pinniped Mbs and investigated their biochemical and structural properties. Comparisons between pinniped and whale Mbs revealed the common and distinctive strategies for the deep-sea adaptation. The overall evolution processes, gaining precipitant tolerance and improving thermodynamic stability, were commonly observed. However, the strategies for improving the folding stability differed, and the pinniped Mbs exploited the shielding of hydrophobic surfaces more effectively than the whale Mbs.

Distal lysine (de)coordination in the algal hemoglobin THB1: A combined computer simulation and experimental study

THB1 is a monomeric truncated hemoglobin from the green alga Chlamydomonas reinhardtii. In the absence of exogenous ligands and at neutral pH, the heme group of THB1 is coordinated by two protein residues, Lys53 and His77. THB1 is thought to function as a nitric oxide dioxygenase, and the distal binding of O₂ requires the cleavage of the Fe-Lys53 bond accompanied by protonation and expulsion of the lysine from the heme cavity into the solvent. Nuclear magnetic resonance spectroscopy and crystallographic data have provided dynamic and structural insights of the process, but the details of the mechanism have not been fully elucidated. We applied a combination of computer simulations and site-directed mutagenesis experiments to shed light on this issue. Molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics restrained optimizations were performed to explore the nature of the transition between the decoordinated and lysine-bound states of the ferrous heme in THB1. Lys49 and Arg52, which form ionic interactions with the heme propionates in the X-ray structure of lysine-bound THB1, were observed to assist in maintaining Lys53 inside the protein cavity and play a key role in the transition. Lys49Ala, Arg52Ala and Lys49Ala/Arg52Ala THB1 variants were prepared, and the consequences of the replacements on the Lys (de)coordination equilibrium were characterized experimentally for comparison with computational prediction. The results reinforced the dynamic role of protein-propionate interactions and strongly suggested that cleavage of the Fe-Lys53 bond and ensuing conformational rearrangement is facilitated by protonation of the amino group inside the distal cavity.

An Engineered Glutamate in Biosynthetic Models of Heme-Copper Oxidases Drives Complete Product Selectivity by Tuning the Hydrogen-Bonding Network

Efficiently carrying out the oxygen reduction reaction (ORR) is critical for many applications in biology and chemistry, such as bioenergetics and fuel cells, respectively. In biology, this reaction is carried out by large, transmembrane oxidases such as heme-copper oxidases (HCOs) and cytochrome bd oxidases. Common to these oxidases is the presence of a glutamate residue next to the active site, but its precise role in regulating the oxidase activity remains unclear. To gain insight into its role, we herein report that incorporation of glutamate next to a designed heme-copper center in two biosynthetic models of HCOs improves O2 binding affinity, facilitates protonation of reaction intermediates, and eliminates release of reactive oxygen species. High-resolution crystal structures of the models revealed extended, water-mediated hydrogen-bonding networks involving the glutamate. Electron paramagnetic resonance of the cryoreduced oxy-ferrous centers at cryogenic temperature followed by thermal annealing allowed observation of the key hydroperoxo intermediate that can be attributed to the hydrogen-bonding network. By demonstrating these important roles of glutamate in oxygen reduction biochemistry, this work offers deeper insights into its role in native oxidases, which may guide the design of more efficient artificial ORR enzymes or catalysts for applications such as fuel cells.

A reaction pathway to compound 0 intermediates in oxy-myoglobin through interactions with hydrogen sulfide and His64

Myoglobin (Mb) binds oxygen with high affinity as a low spin singlet complex and thus functions as an oxygen storage protein. Yet, hybrid Density Functional Theory/Molecular Mechanical (DFT/MM) calculations of oxy-Mb models predict that the O₂ bond is much less resistant to breaking in the presence of hydrogen sulfide (H₂S) compared with water. Specifically, a hydrogen atom from H₂S can be transferred to the distal oxygen atom through homolytic cleavage of the S-H bond to form the intermediate Compound (Cpd) 0 structure and a thiyl radical. In the presence of a neutral His64 (Nε protonation, His64-ε) and H₂S, only a metastable Cpd 0 would be formed as the active site is devoid of any additional proton donor to fully break the O₂ bond. In contrast, the calculations predict that the triplet state is significantly favored over the open shell singlet diradical state throughout the entire reaction coordinate in the presence of H₂S and a positively charged His64. Furthermore, a positively charged His64 can readily donate a proton to Cpd 0 to fully break the O₂ bond resulting in a configuration analogous to reported reaction models of a hemoglobin mutant bound to H₂O₂ with H₂S present. Typically, exotic techniques are required to generate Cpd 0 but under the conditions just described the intermediate is readily detected in UV-Vis spectra at room temperature. The effect is observed as a 2 nm red shift of the Soret band from 414 nm to 416 nm (pH 5.0, His64-εδ) and from 416 nm to 418 nm (pH 6.6, His64-ε).

Myoglobinopathy is an adult-onset autosomal dominant myopathy with characteristic sarcoplasmic inclusions

Myoglobin, encoded by MB, is a small cytoplasmic globular hemoprotein highly expressed in cardiac myocytes and oxidative skeletal myofibers. Myoglobin binds O_2, facilitates its intracellular transport and serves as a controller of nitric oxide and reactive oxygen species. Here, we identify a recurrent c.292C>T (p.His98Tyr) substitution in MB in fourteen members of six European families suffering from an autosomal dominant progressive myopathy with highly characteristic sarcoplasmic inclusions in skeletal and cardiac muscle. Myoglobinopathy manifests in adulthood with proximal and axial weakness that progresses to involve distal muscles and causes respiratory and cardiac failure. Biochemical characterization reveals that the mutant myoglobin has altered O₂ binding, exhibits a faster heme dissociation rate and has a lower reduction potential compared to wild-type myoglobin. Preliminary studies show that mutant myoglobin may result in elevated superoxide levels at the cellular level. These data define a recognizable muscle disease associated with MB mutation.

Myoglobin is a hemeprotein that reversibly binds oxygen and gives muscle its red color. Here, the authors report a genetic variant in the MB gene that associates with myoglobinopathy, an autosomal dominant progressive myopathy, and altered oxygen binding properties of the mutant protein.

Tracing whale myoglobin evolution by resurrecting ancient proteins

Extant cetaceans, such as sperm whale, acquired the great ability to dive into the ocean depths during the evolution from their terrestrial ancestor that lived about 50 million years ago. Myoglobin (Mb) is highly concentrated in the myocytes of diving animals, in comparison with those of land animals, and is thought to play a crucial role in their adaptation as the molecular aqualung. Here, we resurrected ancestral whale Mbs, which are from the common ancestor between toothed and baleen whales (Basilosaurus), and from a further common quadrupedal ancestor between whale and hippopotamus (Pakicetus). The experimental and theoretical analyses demonstrated that whale Mb adopted two distinguished strategies to increase the protein concentration in vivo along the evolutionary history of deep sea adaptation; gaining precipitant tolerance in the early phase of the evolution, and increase of folding stability in the late phase.

Roles of N- and C-terminal domains in the ligand-binding properties of cytoglobin

Cytoglobin (Cygb) is a member of the hexacoordinated globin protein family and is expressed ubiquitously in rat and human tissues. Although Cygb is reportedly upregulated under hypoxic conditions both in vivo and in vitro, suggesting a physiological function to protect cells under hypoxic/ischemic conditions by scavenging reactive oxygen species or by signal transduction, the mechanisms associated with this function have not been fully elucidated. Recent studies comparing Cygbs among several species suggest that mammalian Cygbs show a distinctly longer C-terminal domain potentially involved in unique physiological functions. In this study, we prepared human Cygb mutants (ΔC, ΔN, and ΔNC) with either one or both terminal domains truncated and investigated the enzymatic functions and structural features by spectroscopic methods. Evaluation of the superoxide-scavenging activity between Cygb variants showed that the ΔC and ΔNC mutants exhibited slightly higher activity involving superoxide scavenging as compared with wild-type Cygb. Subsequent experiments involving ligand titration, flash photolysis, and resonance Raman spectroscopic studies suggested that the truncation of the C- and N-terminal domains resulted in less effective to dissociation constants and binding rates for carbon monoxide, respectively. Furthermore, structural stability was assessed by guanidine hydrochloride and revealed that the C-terminal domain might play a vital role in improving structure, whereas the N-terminal domain did not exert a similar effect. These findings indicated that long terminal domains could be important not only in regulating enzymatic activity but also for structural stability, and that the domains might be relevant to other hypothesized physiological functions for Cygb.

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.

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

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