Autocatalytic Mechanism in the Anaerobic Reduction of Metmyoglobin by Sulfide Species

The mechanism of the metal centered reduction of metmyoglobin (MbFeIII) by sulfide species (H2S/HS-) under an argon atmosphere has been studied by a combination of spectroscopic, kinetic, and computational methods. Asymmetric S-shaped time-traces for the formation of MbFeII at varying ratios of excess sulfide were observed at pH 5.3 < pH < 8.0 and 25 °C, suggesting an autocatalytic reaction mechanism. An increased rate at more alkaline pHs points to HS- as relevant reactive species for the reduction. The formation of the sulfanyl radical (HS•) in the slow initial phase was assessed using the spin-trap phenyl N-tert-butyl nitrone. This radical initiates the formation of S-S reactive species as disulfanuidyl/ disulfanudi-idyl radical anions and disulfide (HSSH•-/HSS•2- and HSS-, respectively). The autocatalysis has been ascribed to HSS-, formed after HSSH•-/HSS•2- disproportionation, which behaves as a fast reductant toward the intermediate complex MbFeIII(HS-). We propose a reaction mechanism for the sulfide-mediated reduction of metmyoglobin where only ferric heme iron initiates the oxidation of sulfide species. Beside the chemical interest, this insight into the MbFeIII/sulfide reaction under an argon atmosphere is relevant for the interpretation of biochemical aspects of ectopic myoglobins found on hypoxic tissues toward reactive sulfur species.

Computational Investigation of a Series of Small Molecules as Potential Compounds for Lysyl Hydroxylase-2 (LH2) Inhibition

The catalytic function of lysyl hydroxylase-2 (LH2), a member of the Fe(II)/αKG-dependent oxygenase superfamily, is to catalyze the hydroxylation of lysine to hydroxylysine in collagen, resulting in stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs). Reports show that high amounts of LH2 lead to the accumulation of HLCCs, causing fibrosis and specific types of cancer metastasis. Some members of the Fe(II)/αKG-dependent family have also been reported to have intramolecular O2 tunnels, which aid in transporting one of the required cosubstrates into the active site. While LH2 can be a promising target to combat these diseases, efficacious inhibitors are still lacking. We have used computational simulations to investigate a series of 44 small molecules as lead compounds for LH2 inhibition. Tunneling analyses indicate the existence of several intramolecular tunnels. The lengths of the calculated O2-transporting tunnels in holoenzymes are relatively longer than those in the apoenzyme, suggesting that the ligands may affect the enzyme's structure and possibly block (at least partially) the tunnels. The sequence alignment analysis between LH enzymes from different organisms shows that all of the amino acid residues with the highest occurrence rate in the oxygen tunnels are conserved. Our results suggest that the enolate form of diketone compounds establishes stronger interactions with the Fe(II) in the active site. Branching the enolate compounds with functional groups such as phenyl and pyridinyl enhances the interaction with various residues around the active site. Our results provide information about possible leads for further LH2 inhibition design and development.

Dynamic multiple-scattering treatment of X-ray absorption: Parameterization of a new molecular dynamics force field for myoglobin

We present a detailed analysis of the X-ray absorption near-edge structure (XANES) data on the Fe K-edge of CO Myoglobin based on a combined procedure of Molecular Dynamics (MD) calculations and MXAN (Minuit XANes) data analysis that we call D-MXAN. The ability of performing quantitative XANES data analysis allows us to refine classical force field MD parameters, thus obtaining a reliable tool for the atomic investigation of this important model system for biological macromolecules. The iterative procedure here applied corrects the greatest part of the structural discrepancy between classical MD sampling and experimental determinations. Our procedure, moreover, is able to discriminate between different heme conformational basins visited during the MD simulation, thus demonstrating the necessity of a sampling on the order of tens of nanoseconds, even for an application such X-ray absorption spectroscopy data analysis.

Computational investigation of O2 diffusion through an intra-molecular tunnel in AlkB; influence of polarization on O2 transport

E. Coli AlkB catalyzes the direct dealkylation of various alkylated bases in damaged DNA. The diffusion of molecular oxygen to the active site in AlkB is an essential step for the oxidative dealkylation activity. Despite detailed studies on the stepwise oxidation mechanism of AlkB, there is no conclusive picture of how O2 molecules reach the active site of the protein. Yu et al. (Nature, 439, 879) proposed the existence of an intra-molecular tunnel based on their initial crystal structures of AlkB. We have employed computational simulations to investigate possible migration pathways inside AlkB for O2 molecules. Extensive molecular dynamics (MD) simulations, including explicit ligand sampling and potential of mean force (PMF) calculations, have been performed to provide a microscopic description of the O2 delivery pathway in AlkB. Analysis of intra-molecular tunnels using the CAVER software indicates two possible pathways for O2 to diffuse into the AlkB active site. Explicit ligand sampling simulations suggests that only one of these tunnels provides a viable route. The free energy path for an oxygen molecule to travel along each of these tunnels has been determined with AMBER and AMOEBA. Both PMFs indicate passive transport of O2 from the surface of the protein. However, the inclusion of explicit polarization shows a very large barrier for diffusion of the co-substrate out of the active site, compared with the non-polarizable potential. In addition, our results suggest that the mutation of a conserved residue along the tunnel, Y178, has dramatic effects on the dynamics of AlkB and on the transport of O2 along the tunnel.

Differential Interaction of Myoglobin with Select Fatty Acids of Carbon Chain Lengths C8 to C16

Previous studies have shown that palmitic acid (PAM) and oleic acid (OLE) can bind myoglobin (Mb). How fatty acids (FA) with different carbon chain lengths and sulfate substitution interact with Mb remains uncertain. Indeed, C8:0 and C10:0 fatty acids do not perturb the intensities of the ¹H-NMR MbCN signal intensity at FA:Mb ratios below 2:1. Starting with C12:0, C12:0-C16:0, FA induce a noticeable spectral change. C12:0 and C14:0 FA affect both the 5- and 8-heme methyl signals, whereas the C16:0 FA perturbs only the 8-heme methyl signal. All C12:0-C16:0 saturated FA induce upfield shifts in the -CH₂ peak of different FA in the presence of Mb. Increasing the apparent solubility with a sulfate group substitution enhances the FA interaction of lauric sulfate (LAU 1-SO₄) but not palmitate sulfate acid (PAM 1-SO₄). The detergent (DET) property of FA has no significant contribution. Common positive, neutral, and negative DET at DET:Mb ratio of 1:1 induce no perturbation of the MbCN spectra. The experiment observations establish a basis to investigate the molecular mechanism underlying the FA interaction with Mb.

Insight into the phylogenetic relationship and structural features of vertebrate myoglobin family

Myoglobin (Mb) is studied to clarify the structure-function relationships in protein science. In this work, we report the results of a comparative analysis of amino acid sequences from 298 vertebrate Mbs. Forty-one high conserved residues were identified and seven of them were invariants [E18, G25, F43, V68, L72, H93 (proximal histidine) and H97]. E18 is the only invariant amino acid residue located out of the heme-pocket and Xe-cavities playing a role in interaction between the A and E-helices. A comparative analysis of several parameters related to amino acid composition shows an increase of average mass, accessible surface area and volume per residue from Actinopterygii to Mammalia and Aves. This may be due to an increased number of bulky residues reducing the non-specific cavities volume and thus improving the oxygen flow between the heme site and the outside of the protein. Finally, the phylogenetic analyses of Mb in vertebrates are consistent with an evolution that runs with the diversification of the species, but in which several episodes of gene duplication and lost have occurred, less frequently in the ancestors of great taxons, cartilaginous fishes and non-avian reptiles, most frequently in ray-finned fishes and mammals, and very frequently in birds.

Sulfmyoglobin Conformational Change: A Role in the Decrease of Oxy-Myoglobin Functionality

This work is focused at understanding the interaction of H2S with Myoglobin (Mb), in particular the Sulfmyoglobin (SMb) product, whose physiological role is controversial and not well understood. The scattering curves, Guinier, Kratky, Porod and P(r) plots were analyzed for oxy-Mb and oxy-Hemoglobin I (oxyHbI) in the absence and presence of H2S, using Small and Wide Angle X-ray Scattering (SAXS/WAXS) technique. Three dimensional models were also generated from the SAXS/WAXS data. The results show that SMb formation, produced by oxyMb and H2S interaction, induces a change in the protein conformation where its envelope has a very small cleft and the protein is more flexible, less rigid and compact. Based on the direct relationship between Mb's structural conformation and its functionality, we suggest that the conformational change observed upon SMb formation plays a contribution to the protein decrease in O2 affinity and, therefore, on its functionality.

Benchmarking Rapid TLES Simulations of Gas Diffusion in Proteins: Mapping O2 Migration and Escape in Myoglobin as a Case Study

Standard molecular dynamics (MD) simulations of gas diffusion consume considerable computational time and resources even for small proteins. To combat this, temperature-controlled locally enhanced sampling (TLES) examines multiple diffusion trajectories per simulation by accommodating multiple noninteracting copies of a gas molecule that diffuse independently, while the protein and water molecules experience an average interaction from all copies. Furthermore, gas migration within a protein matrix can be accelerated without altering protein dynamics by increasing the effective temperature of the TLES copies. These features of TLES enable rapid simulations of gas diffusion within a protein matrix at significantly reduced (∼98%) computational cost. However, the results of TLES and standard MD simulations have not been systematically compared, which limits the adoption of the TLES approach. We address this drawback here by benchmarking TLES against standard MD in the simulation of O2 diffusion in myoglobin (Mb) as a case study since this model system has been extensively characterized. We find that 2 ns TLES and 108 ns standard simulations map the same network of diffusion tunnels in Mb and uncover the same docking sites, barriers, and escape portals. We further discuss the influence of simulation time as well as the number of independent simulations on the O2 population density within the diffusion tunnels and on the sampling of Mb's conformational space as revealed by principal component analysis. Overall, our comprehensive benchmarking reveals that TLES is an appropriate and robust tool for the rapid mapping of gas diffusion in proteins when the kinetic data provided by standard MD are not required. Furthermore, TLES provides explicit ligand diffusion pathways, unlike most rapid methods.

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.

Myoglobin oxygen affinity in aquatic and terrestrial birds and mammals

Myoglobin (Mb) is an oxygen binding protein found in vertebrate skeletal muscle, where it facilitates intracellular transport and storage of oxygen. This protein has evolved to suit unique physiological needs in the muscle of diving vertebrates that express Mb at much greater concentrations than their terrestrial counterparts. In this study, we characterized Mb oxygen affinity (P50) from 25 species of aquatic and terrestrial birds and mammals. Among diving species, we tested for correlations between Mb P50 and routine dive duration. Across all species examined, Mb P50 ranged from 2.40 to 4.85 mmHg. The mean P50 of Mb from terrestrial ungulates was 3.72±0.15 mmHg (range 3.70-3.74 mmHg). The P50 of cetaceans was similar to terrestrial ungulates ranging from 3.54 to 3.82 mmHg, with the exception of the melon-headed whale, which had a significantly higher P50 of 4.85 mmHg. Among pinnipeds, the P50 ranged from 3.23 to 3.81 mmHg and showed a trend for higher oxygen affinity in species with longer dive durations. Among diving birds, the P50 ranged from 2.40 to 3.36 mmHg and also showed a trend of higher affinities in species with longer dive durations. In pinnipeds and birds, low Mb P50 was associated with species whose muscles are metabolically active under hypoxic conditions associated with aerobic dives. Given the broad range of potential globin oxygen affinities, Mb P50 from diverse vertebrate species appears constrained within a relatively narrow range. High Mb oxygen affinity within this range may be adaptive for some vertebrates that make prolonged dives.

Non-site-specific allosteric effect of oxygen on human hemoglobin under high oxygen partial pressure

Protein allostery is essential for vital activities. Allosteric regulation of human hemoglobin (HbA) with two quaternary states T and R has been a paradigm of allosteric structural regulation of proteins. It is widely accepted that oxygen molecules (O₂) act as a “site-specific” homotropic effector, or the successive O₂ binding to the heme brings about the quaternary regulation. However, here we show that the site-specific allosteric effect is not necessarily only a unique mechanism of O₂ allostery. Our simulation results revealed that the solution environment of high O₂ partial pressure enhances the quaternary change from T to R without binding to the heme, suggesting an additional “non-site-specific” allosteric effect of O₂. The latter effect should play a complementary role in the quaternary change by affecting the intersubunit contacts. This analysis must become a milestone in comprehensive understanding of the allosteric regulation of HbA from the molecular point of view.

Dynamic void distribution in myoglobin and five mutants

Globular proteins contain cavities/voids that play specific roles in controlling protein function. Elongated cavities provide migration channels for the transport of ions and small molecules to the active center of a protein or enzyme. Using Monte Carlo and Molecular Dynamics on fully atomistic protein/water models, a new computational methodology is introduced that takes into account the protein's dynamic structure and maps all the cavities in and on the surface. To demonstrate its utility, the methodology is applied to study cavity structure in myoglobin and five of its mutants. Computed cavity and channel size distributions reveal significant differences relative to the wild type myoglobin. Computer visualization of the channels leading to the heme center indicates restricted ligand access for the mutants consistent with the existing interpretations. The new methodology provides a quantitative measure of cavity structure and distributions and can become a valuable tool for the structural characterization of proteins.

An atomistic view on human hemoglobin carbon monoxide migration processes

A significant amount of work has been devoted to obtaining a detailed atomistic knowledge of the human hemoglobin mechanism. Despite this impressive research, to date, the ligand diffusion processes remain unclear and controversial. Using recently developed computational techniques, PELE, we are capable of addressing the ligand migration processes. First, the methodology was tested on myoglobin's CO migration, and the results were compared with the wealth of theoretical and experimental studies. Then, we explored both hemoglobin tense and relaxed states and identified the differences between the α-and β-subunits. Our results indicate that the proximal site, equivalent to the Xe1 cavity in myoglobin, is never visited. Furthermore, strategically positioned residues alter the diffusion processes within hemoglobin's subunits and suggest that multiple pathways exist, especially diversified in the α-globins. A significant dependency of the ligand dynamics on the tertiary structure is also observed.

Endogenous myoglobin in breast cancer is hypoxia-inducible by alternative transcription and functions to impair mitochondrial activity: a role in tumor suppression?

Recently, immunohistochemical analysis of myoglobin (MB) in human breast cancer specimens has revealed a surprisingly widespread expression of MB in this nonmuscle context. The positive correlation with hypoxia-inducible factor 2α (HIF-2α) and carbonic anhydrase IX suggested that oxygen regulates myoglobin expression in breast carcinomas. Here, we report that MB mRNA and protein levels are robustly induced by prolonged hypoxia in breast cancer cell lines, in part via HIF-1/2-dependent transactivation. The hypoxia-induced MB mRNA originated from a novel alternative transcription start site 6 kb upstream of the ATG codon. MB regulation in normal and tumor tissue may thus be fundamentally different. Functionally, the knockdown of MB in MDA-MB468 breast cancer cells resulted in an unexpected increase of O(2) uptake and elevated activities of mitochondrial enzymes during hypoxia. Silencing of MB transcription attenuated proliferation rates and motility capacities of hypoxic cancer cells and, surprisingly, also fully oxygenated breast cancer cells. Endogenous MB in cancer cells is apparently involved in controlling oxidative cell energy metabolism, contrary to earlier findings on mouse heart, where the targeted disruption of the Mb gene did not effect myocardial energetics and O(2) consumption. This control function of MB seemingly impacts mitochondria and influences cell proliferation and motility, but it does so in ways not directly related to the facilitated diffusion or storage of O(2). Hypothetically, the mitochondrion-impairing role of MB in hypoxic cancer cells is part of a novel tumor-suppressive function.

Time resolved thermodynamics associated with ligand photorelease in heme peroxidases and globins: Open access channels versus gated ligand release

Heme proteins represent a diverse class of biomolecules responsible for an extremely diverse array of physiological functions including electron transport, monooxygenation, ligand transport and storage, cellular signaling, respiration, etc. An intriguing aspect of these proteins is that such functional diversity is accomplished using a single type of heme macrocycle based upon iron protoporphyrin IX. The functional diversity originates from a delicate balance of inter-molecular interactions within the protein matrix together with well choreographed dynamics that modulate the heme electronic structure as well as ligand entry/exit pathways from the bulk solvent to the active site. Of particular interest are the dynamics and energetics associated with the entry/exit of ligands as this process plays a significant role in regulating the rates of heme protein activity. Time-resolved photoacoustic calorimetry (PAC) has emerged as a powerful tool through which to probe the underlying energetics associated with small molecule dissociation and release to the bulk solvent in heme proteins on time scales from tens of nanoseconds to several microseconds. In this review, the results of PAC studies on various classes of heme proteins are summarized highlighting how different protein structures affect the thermodynamics of ligand migration. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.