High pressure reveals structural determinants for globin hexacoordination: Neuroglobin and myoglobin cases

Redox and spectroscopic properties of mammalian nitrite reductase-like hemoproteins

Besides the canonical pathway of L-arginine oxidation to produce nitric oxide (NO) in vivo, the nitrate-nitrite-NO pathway has been widely accepted as another source for circulating NO in mammals, especially under hypoxia. To date, there have been at least ten heme-containing nitrite reductase-like proteins discovered in mammals with activities mainly identified in vitro, including four globins (hemoglobin, myoglobin, neuroglobin (Ngb), cytoglobin (Cygb)), three mitochondrial respiratory chain enzymes (cytochrome c oxidase, cytochrome bc₁, cytochrome c), and three other heme proteins (endothelial nitric oxide synthase, cytochrome P450 and indoleamine 2,3-dioxygenase 1 (IDO1)). The pathophysiological functions of these proteins are closely related to their redox and spectroscopic properties, as well as their protein structure, although the physiological roles of Ngb, Cygb and IDO1 remain unclear. So far, comprehensive summaries of the redox and spectroscopic properties of these nitrite reductase-like hemoproteins are still lacking. In this review, we have mainly summarized the published data on the application of ultraviolet-visible, electron paramagnetic resonance, circular dichroism and resonance Raman spectroscopies, and X-ray crystallography in studying nitrite reductase-like activity of these 10 proteins, in order to sort out the relationships among enzymatic function, structure and spectroscopic characterization, which might help in understanding their roles in redox biology and medicine.

Structural and (Pseudo-)Enzymatic Properties of Neuroglobin: Its Possible Role in Neuroprotection

Neuroglobin (Ngb), the third member of the globin family, was discovered in human and murine brains in 2000. This monomeric globin is structurally similar to myoglobin (Mb) and hemoglobin (Hb) α and β subunits, but it hosts a bis-histidyl six-coordinated heme-Fe atom. Therefore, the heme-based reactivity of Ngb is modulated by the dissociation of the distal HisE7-heme-Fe bond, which reflects in turn the redox state of the cell. The high Ngb levels (~100–200 μM) present in the retinal ganglion cell layer and in the optic nerve facilitate the O₂ buffer and delivery. In contrast, the very low levels of Ngb (~1 μM) in most tissues and organs support (pseudo-)enzymatic properties including NO/O₂ metabolism, peroxynitrite and free radical scavenging, nitrite, hydroxylamine, hydrogen sulfide reduction, and the nitration of aromatic compounds. Here, structural and (pseudo-)enzymatic properties of Ngb, which are at the root of tissue and organ protection, are reviewed, envisaging a possible role in the protection from neuronal degeneration of the retina and the optic nerve.

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.

Control of distal lysine coordination in a monomeric hemoglobin: A role for heme peripheral interactions

THB1 is a monomeric truncated hemoglobin (TrHb) found in the cytoplasm of the green alga Chlamydomonas reinhardtii. The canonical heme coordination scheme in hemoglobins is a proximal histidine ligand and an open distal site. In THB1, the latter site is occupied by Lys53, which is likely to facilitate Fe(II)/Fe(III) redox cycling but hinders dioxygen binding, two features inherent to the NO dioxygenase activity of the protein. TrHb surveys show that a lysine at a position aligning with Lys53 is an insufficient determinant of coordination, and in this study, we sought to identify factors controlling lysine affinity for the heme iron. We solved the "Lys-off" X-ray structure of THB1, represented by the cyanide adduct of the Fe(III) protein, and hypothesized that interactions that differ between the known "Lys-on" structure and the Lys-off structure participate in the control of Lys53 affinity for the heme iron. We applied an experimental approach (site-directed mutagenesis, heme modification, pH titrations in the Fe(III) and Fe(II) states) and a computational approach (MD simulations in the Fe(II) state) to assess the role of heme propionate-protein interactions, distal helix capping, and the composition of the distal pocket. All THB1 modifications resulted in a weakening of lysine affinity and affected the coupling between Lys53 proton binding and heme redox potential. The results supported the importance of specific heme peripheral interactions for the pH stability of iron coordination and the ability of the protein to undergo redox reactions.

Structural Insights into the Heme Pocket and Oligomeric State of Non-Symbiotic Hemoglobins from Arabidopsis thaliana

Non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana are hexacoordinate heme-proteins that likely have different biological roles, in view of diverse tissue localization, expression pattern, and ligand binding properties. Herein, we expand upon previous biophysical studies on these isoforms, focusing on their oligomeric states and circular dichroism (CD) characteristics. We found that AHb1 exists in solution in a concentration-dependent monomer-dimer equilibrium, while AHb2 is present only as a monomer. The quaternary structure of AHb1 affects its degree of hexacoordination with the formation of the dimer that enhances pentacoordination. Accordingly, the mutant of a conserved residue within the dimeric interface, AHb1-T45A, which is mostly monomeric in solution, has an equilibrium that is shifted toward a hexacoordinate form compared to the wild-type protein. CD studies further support differences in the globin's structure and heme moiety. The Soret CD spectra for AHb2 are opposite in sense to those for AHb1, reflecting different patterns of heme-protein side chain contacts in the two proteins. Moreover, the smaller contribution of the heme to the near-UV CD in AHb2 compared to AHb1 suggests a weaker heme-protein association in AHb2. Our data corroborate the structural diversity of AHb1 and AHb2 and confirm the leghemoglobin-like structural properties of AHb2.

Structural and functional properties of Antarctic fish cytoglobins-1: Cold-reactivity in multi-ligand reactions

While the functions of the recently discovered cytoglobin, ubiquitously expressed in vertebrate tissues, remain uncertain, Antarctic fish provide unparalleled models to study novel protein traits that may arise from cold adaptation. We report here the spectral, ligand-binding and enzymatic properties (peroxynitrite isomerization, nitrite-reductase activity) of cytoglobin-1 from two Antarctic fish, Chaenocephalus aceratus and Dissostichus mawsoni, and present the crystal structure of D. mawsoni cytoglobin-1. The Antarctic cytoglobins-1 display high O₂ affinity, scarcely compatible with an O₂-supply role, a slow rate constant for nitrite-reductase activity, and do not catalyze peroxynitrite isomerization. Compared with mesophilic orthologues, the cold-adapted cytoglobins favor binding of exogenous ligands to the hexa-coordinated bis-histidyl species, a trait related to their higher rate constant for distal-His/heme-Fe dissociation relative to human cytoglobin. At the light of a remarkable 3D-structure conservation, the observed differences in ligand-binding kinetics may reflect Antarctic fish cytoglobin-1 specific features in the dynamics of the heme distal region and of protein matrix cavities, suggesting adaptation to functional requirements posed by the cold environment. Taken together, the biochemical and biophysical data presented suggest that in Antarctic fish, as in humans, cytoglobin-1 unlikely plays a role in O₂ transport, rather it may be involved in processes such as NO detoxification.

Adaptations to environmental change: Globin superfamily evolution in Antarctic fishes

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

Structure of the Human ACP-ISD11 Heterodimer

In recent years, the mammalian mitochondrial protein complex for iron-sulfur cluster assembly has been the focus of important studies. This is partly because of its high degree of relevance in cell metabolism and because mutations of the involved proteins are the cause of several human diseases. Cysteine desulfurase NFS1 is the key enzyme of the complex. At present, it is well-known that the active form of NFS1 is stabilized by the small protein ISD11. In this work, the structure of the human mitochondrial ACP-ISD11 heterodimer was determined at 2.0 Å resolution. ACP-ISD11 forms a cooperative unit stabilized by several ionic interactions, hydrogen bonds, and apolar interactions. The 4'-phosphopantetheine-acyl chain, which is covalently bound to ACP, interacts with several residues of ISD11, modulating together with ACP the foldability of ISD11. Recombinant human ACP-ISD11 was able to interact with the NFS1 desulfurase, thus yielding an active enzyme, and the NFS1/ACP-ISD11 core complex was activated by frataxin and ISCU proteins. Internal motions of ACP-ISD11 were studied by molecular dynamics simulations, showing the persistence of the interactions between both protein chains. The conformation of the dimer is similar to that found in the context of the (NFS1/ACP-ISD11)₂ supercomplex core, which contains the Escherichia coli ACP instead of the human variant. This fact suggests a sequential mechanism for supercomplex consolidation, in which the ACP-ISD11 complex may fold independently and, after that, the NFS1 dimer would be stabilized.

Thermal Stability of Globins: Implications of Flexibility and Heme Coordination Studied by Molecular Dynamics Simulations

Proteins are sensitive to temperature, and abrupt changes in the normal temperature conditions can have a profound impact on both structure and function, leading to protein unfolding. However, the adaptation of certain organisms to extreme conditions raises questions about the structural features that permit the structure and function of proteins to be preserved under these adverse conditions. To gain insight into the molecular basis of protein thermostability in the globin family, we have examined three representative examples: human neuroglobin, horse heart myoglobin, and Drosophila hemoglobin, which differ in their melting temperatures and coordination states of the heme iron in the absence of external ligands. In order to elucidate the possible mechanisms that govern the thermostability of these proteins, microsecond-scale classical molecular dynamics simulations were performed at different temperatures. Structural fluctuations and essential dynamics were analyzed, indicating that the flexibility of the CD region, which includes the two short C and D helixes and the connecting CD loop, is directly related to the thermostability. We observed that a larger inherent flexibility of the protein produces higher thermostability, probably concentrating the thermal fluctuations observed at high temperature in flexible regions, preventing unfolding. Globally, the results of this work improve our understanding of thermostability in the globin family.

Glyceraldehyde-3-phosphate dehydrogenase is a chaperone that allocates labile heme in cells

Cellular heme is thought to be distributed between a pool of sequestered heme that is tightly bound within hemeproteins and a labile heme pool required for signaling and transfer into proteins. A heme chaperone that can hold and allocate labile heme within cells has long been proposed but never been identified. Here, we show that the glycolytic protein glyceraldehyde-3-phosphate dehydrogenase (GAPDH) fulfills this role by acting as an essential repository and allocator of bioavailable heme to downstream protein targets. We identified a conserved histidine in GAPDH that is needed for its robust heme binding both in vitro and in mammalian cells. Substitution of this histidine, and the consequent decreases in GAPDH heme binding, antagonized heme delivery to both cytosolic and nuclear hemeprotein targets, including inducible nitric-oxide synthase (iNOS) in murine macrophages and the nuclear transcription factor Hap1 in yeast, even though this GAPDH variant caused cellular levels of labile heme to rise dramatically. We conclude that by virtue of its heme-binding property, GAPDH binds and chaperones labile heme to create a heme pool that is bioavailable to downstream proteins. Our finding solves a fundamental question in cell biology and provides a new foundation for exploring heme homeostasis in health and disease.

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.

Structural Plasticity in Globins: Role of Protein Dynamics in Defining Ligand Migration Pathways

Globins are a family of proteins characterized by the presence of the heme prosthetic group and involved in variety of biological functions in the cell. Due to their biological relevance and widespread distribution in all kingdoms of life, intense research efforts have been devoted to disclosing the relationships between structural features, protein dynamics, and function. Particular attention has been paid to the impact of differences in amino acid sequence on the topological features of docking sites and cavities and to the influence of conformational flexibility in facilitating the migration of small ligands through these cavities. Often, tunnels are carved in the interior of globins, and ligand exchange is regulated by gating residues. Understanding the subtle intricacies that relate the differences in sequence with the structural and dynamical features of globins with the ultimate aim of rationalizing the thermodynamics and kinetics of ligand binding continues to be a major challenge in the field. Due to the evolution of computational techniques, significant advances into our understanding of these questions have been made. In this review we focus our attention on the analysis of the ligand migration pathways as well as the function of the structural cavities and tunnels in a series of representative globins, emphasizing the synergy between experimental and theoretical approaches to gain a comprehensive knowledge into the molecular mechanisms of this diverse family of proteins.

Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species

Background

The identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures.

Results

The current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified.

Conclusions

The results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.

Molecular basis of thermal stability in truncated (2/2) hemoglobins

BACKGROUND

Understanding the molecular mechanism through which proteins are functional at extreme high and low temperatures is one of the key issues in structural biology. To investigate this phenomenon, we have focused on two instructive truncated hemoglobins from Thermobifida fusca (Tf-trHbO) and Mycobacterium tuberculosis (Mt-trHbO); although the two proteins are structurally nearly identical, only the former is stable at high temperatures.

METHODS

We used molecular dynamics simulations at different temperatures as well as thermal melting profile measurements of both wild type proteins and two mutants designed to interchange the amino acid residue, either Pro or Gly, at E3 position.

RESULTS

The results show that the presence of a Pro at the E3 position is able to increase (by 8°) or decrease (by 4°) the melting temperature of Mt-trHbO and Tf-trHbO, respectively. We observed that the ProE3 alters the structure of the CD loop, making it more flexible.

CONCLUSIONS

This gain in flexibility allows the protein to concentrate its fluctuations in this single loop and avoid unfolding. The alternate conformations of the CD loop also favor the formation of more salt-bridge interactions, together augmenting the protein's thermostability.

GENERAL SIGNIFICANCE

These results indicate a clear structural and dynamical role of a key residue for thermal stability in truncated hemoglobins.

QM/MM study of the C-C coupling reaction mechanism of CYP121, an essential cytochrome p450 of Mycobacterium tuberculosis

Among 20 p450s of Mycobacterium tuberculosis (Mt), CYP121 has received an outstanding interest, not only due to its essentiality for bacterial viability but also because it catalyzes an unusual carbon-carbon coupling reaction. Based on the structure of the substrate bound enzyme, several reaction mechanisms were proposed involving first Tyr radical formation, second Tyr radical formation, and C-C coupling. Key and unknown features, being the nature of the species that generate the first and second radicals, and the role played by the protein scaffold each step. In the present work we have used classical and quantum based computer simulation methods to study in detail its reaction mechanism. Our results show that substrate binding promotes formation of the initial oxy complex, Compound I is the responsible for first Tyr radical formation, and that the second Tyr radical is formed subsequently, through a PCET reaction, promoted by the presence of key residue Arg386. The final C-C coupling reaction possibly occurs in bulk solution, thus yielding the product in one oxygen reduction cycle. Our results thus contribute to a better comprehension of MtCYP121 reaction mechanism, with direct implications for inhibitor design, and also contribute to our general understanding of these type of enzymes.

Expression and cell distribution of neuroglobin in the brain tissue after experimental subarachnoid hemorrhage in rats: a pilot study

Neuroglobin (Ngb) is a member of the globin superfamily expressed mainly in the nervous system and retina of vertebrates. Accumulated evidence has clearly demonstrated that Ngb has a neuro-protective role enhancing cell viability under hypoxia and other types of oxidative stress. It was suggested that oxidant stress could play an important role in neuronal injury after subarachnoid hemorrhage (SAH). The present study aims to examine the expression of Ngb in the temporal cortex and its cellular localization after SAH. We used a prechiasmatic cistern model of SAH. Ngb expression was examined at 3, 6, 12, 24, 48, and 72 h after SAH by western blot analysis and real-time polymerase chain reaction (PCR). Immunohistochemistry and immunofluorescence were performed to detect the localization of Ngb. Real-time PCR demonstrated that Ngb mRNA levels increased from 3 h after SAH, peaked at 6 h. Western blot showed Ngb protein levels were significantly increased in SAH groups in the temporal cortex and reached the peak at 24 h after SAH. The immunohistochemical staining demonstrated that Ngb was weakly expressed in the cortex in the control group while the enhanced expression of Ngb could be detected in the SAH groups. In addition, immunofluorescence results revealed that the over-expressed Ngb was located in the neuronal and microglia cell cytoplasm. These findings indicated that Ngb might play an important neuro-protective effect after SAH.

Volumetric properties underlying ligand binding in a monomeric hemoglobin: a high-pressure NMR study

The 2/2 hemoglobin of the cyanobacterium Synechococcus sp. PCC 7002, GlbN, coordinates the heme iron with two histidines and exists either with a b heme or with a covalently attached heme. The binding of exogenous ligands displaces the distal histidine and induces a conformational rearrangement involving the reorganization of internal void volumes. The formation of passageways within the resulting conformation is thought to facilitate ligand exchange and play a functional role. Here we monitored the perturbation induced by pressure on the ferric bis-histidine and cyanide-bound states of GlbN using (1)H-(15)N HSQC NMR spectroscopy. We inspected the outcome with a statistical analysis of 170 homologous 2/2 hemoglobin sequences. We found that the compression landscape of GlbN, as represented by the variation of an average chemical shift parameter, was highly sensitive to ligand swapping and heme covalent attachment. Stabilization of rare conformers was observed at high pressures and consistent with cavity redistribution upon ligand binding. In all states, the EF loop was found to be exceptionally labile to pressure, suggesting a functional role as a semi-flexible hinge between the adjacent helices. Finally, coevolved clusters presented a common pattern of compensating pressure responses. The high-pressure dissection combined with protein sequence analysis established locations with volumetric signatures relevant to residual communication of 2/2 hemoglobins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2

This study reports a comparative analysis of the topological properties of inner cavities and the intrinsic dynamics of non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. The two proteins belong to the 3/3 globin fold and have a sequence identity of about 60%. However, it is widely assumed that they have distinct physiological roles. In order to investigate the structure-function relationships in these proteins, we have examined the bis-histidyl and ligand-bound hexacoordinated states by atomistic simulations using in silico structural models. The results allow us to identify two main pathways to the distal cavity in the bis-histidyl hexacoordinated proteins. Nevertheless, a larger accessibility to small gaseous molecules is found in AHb2. This effect can be attributed to three factors: the mutation Leu35(AHb1)→Phe32(AHb2), the enhanced flexibility of helix B, and the more favorable energetic profile for ligand migration to the distal cavity. The net effect of these factors would be to facilitate the access of ligands, thus compensating the preference for the fully hexacoordination of AHb2, in contrast to the equilibrium between hexa- and pentacoordinated species in AHb1. On the other hand, binding of the exogenous ligand introduces distinct structural changes in the two proteins. A well-defined tunnel is formed in AHb1, which might be relevant to accomplish the proposed NO detoxification reaction. In contrast, no similar tunnel is found in AHb2, which can be ascribed to the reduced flexibility of helix E imposed by the larger number of salt bridges compared to AHb1. This feature would thus support the storage and transport functions proposed for AHb2. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Role of tertiary structures on the Root effect in fish hemoglobins

Many fish hemoglobins exhibit a marked dependence of oxygen affinity and cooperativity on proton concentration, called Root effect. Both tertiary and quaternary effects have been evoked to explain the allosteric regulation brought about by protons in fish hemoglobins. However, no general rules have emerged so far. We carried out a complementary crystallographic and microspectroscopic characterization of ligand binding to crystals of deoxy-hemoglobin from the Antarctic fish Trematomus bernacchii (HbTb) at pH6.2 and pH8.4. At low pH ligation has negligible structural effects, correlating with low affinity and absence of cooperativity in oxygen binding. At high pH, ligation causes significant changes at the tertiary structural level, while preserving structural markers of the T state. These changes mainly consist in a marked displacement of the position of the switch region CD corner towards an R-like position. The functional data on T-state crystals validate the relevance of the crystallographic observations, revealing that, differently from mammalian Hbs, in HbTb a significant degree of cooperativity in oxygen binding is due to tertiary conformational changes, in the absence of the T-R quaternary transition. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Hydrophobic effect drives oxygen uptake in myoglobin via histidine E7

Since the elucidation of the myoglobin (Mb) structure, a histidine residue on the E helix (His-E7) has been proposed to act as a gate with an open or closed conformation controlling access to the active site. Although it is believed that at low pH, the His-E7 gate is in its open conformation, the full relationship between the His-E7 protonation state, its conformation, and ligand migration in Mb is hotly debated. We used molecular dynamics simulations to first address the effect of His-E7 protonation on its conformation. We observed the expected shift from the closed to the open conformation upon protonation, but more importantly, noted a significant difference between the conformations of the two neutral histidine tautomers. We further computed free energy profiles for oxygen migration in each of the possible His-E7 states as well as in two instructive Mb mutants: Ala-E7 and Trp-E7. Our results show that even in the closed conformation, the His-E7 gate does not create a large barrier to oxygen migration and permits oxygen entry with only a small rotation of the imidazole side chain and movement of the E helix. We identify, instead, a hydrophobic site in the E7 channel that can accommodate an apolar diatomic ligand and enhances ligand uptake particularly in the open His-E7 conformation. This rate enhancement is diminished in the closed conformation. Taken together, our results provide a new conceptual framework for the histidine gate hypothesis.

Biophysical characterisation of neuroglobin of the icefish, a natural knockout for hemoglobin and myoglobin. Comparison with human neuroglobin

The Antarctic icefish Chaenocephalus aceratus lacks the globins common to most vertebrates, hemoglobin and myoglobin, but has retained neuroglobin in the brain. This conserved globin has been cloned, over-expressed and purified. To highlight similarities and differences, the structural features of the neuroglobin of this colourless-blooded fish were compared with those of the well characterised human neuroglobin as well as with the neuroglobin from the retina of the red blooded, hemoglobin and myoglobin-containing, closely related Antarctic notothenioid Dissostichus mawsoni. A detailed structural and functional analysis of the two Antarctic fish neuroglobins was carried out by UV-visible and Resonance Raman spectroscopies, molecular dynamics simulations and laser-flash photolysis. Similar to the human protein, Antarctic fish neuroglobins can reversibly bind oxygen and CO in the Fe(2+) form, and show six-coordination by distal His in the absence of exogenous ligands. A very large and structured internal cavity, with discrete docking sites, was identified in the modelled three-dimensional structures of the Antarctic neuroglobins. Estimate of the free-energy barriers from laser-flash photolysis and Implicit Ligand Sampling showed that the cavities are accessible from the solvent in both proteins.Comparison of structural and functional properties suggests that the two Antarctic fish neuroglobins most likely preserved and possibly improved the function recently proposed for human neuroglobin in ligand multichemistry. Despite subtle differences, the adaptation of Antarctic fish neuroglobins does not seem to parallel the dramatic adaptation of the oxygen carrying globins, hemoglobin and myoglobin, in the same organisms.

Neuroglobin regulates hypoxic response of neuronal cells through Hif-1α- and Nrf2-mediated mechanism

Oxygen sensing in hypoxic neurons has been classically attributed to cytochrome c oxidase and prolyl-4-hydroxylases and involves stabilization of transcription factors, hypoxia-inducible factor-1α (Hif-1α) and nuclear factor erythroid 2-related factor 2 (Nrf2) that mediate survival responses. On the contrary, release of cytochrome c into the cytosol during hypoxic stress triggers apoptosis in neuronal cells. We, here advocate that the redox state of neuroglobin (Ngb) could regulate both Hif-1α and Nrf2 stabilization and cytochrome c release during hypoxia. The hippocampal regions showing higher expression of Ngb were less susceptible to global hypoxia-mediated neurodegeneration. During normoxia, Ngb maintained cytochrome c in the reduced state and prevented its release from mitochondria by using cellular antioxidants. Greater turnover of oxidized cytochrome c and increased utilization of cellular antioxidants during acute hypoxia altered cellular redox status and stabilized Hif-1α and Nrf2 through Ngb-mediated mechanism. Chronic hypoxia, however, resulted in oxidation and degradation of Ngb, accumulation of ferric ions and release of cytochrome c that triggered apoptosis. Administration of N-acetyl-cysteine during hypoxic conditions improved neuronal survival by preventing Ngb oxidation and degradation. Taken together, these results establish a role for Ngb in regulating both the survival and apoptotic mechanisms associated with hypoxia.

Structure and dynamics of Antarctic fish neuroglobin assessed by computer simulations

Neuroglobin (Ngb) is a heme protein, highly conserved along evolution, predominantly found in the nervous system. It is upregulated by hypoxia and ischemia and may have a neuroprotective role under hypoxic stress. Although many other roles have been proposed, the physiological function is still unclear. Antarctic icefishes lack hemoglobin and some species also lack myoglobin, but all have Ngb and thus may help the elucidation of Ngb function. We present the first theoretically derived structure of fish Ngb and describe its behavior using molecular dynamics simulations. Specifically, we sequenced and analyzed Ngbs from a colorless-blooded Antarctic icefish species Chaenocephalus aceratus and a related red-blooded species (Dissostichus mawsoni). Both fish Ngbs are 6-coordinated but have some peculiarities that differentiate them from mammalian counterparts: they have extensions in the N and C termini that can interact with the EF loop, and a gap in the alignment that changes the CD-region structure/dynamics that has been found to play a key role in human neuroglobin. Our results suggest that a single mutation between both fish Ngbs is responsible for significant difference in the behavior of the proteins. The functional role of these characteristics is discussed.

An integrated computational analysis of the structure, dynamics, and ligand binding interactions of the human galectin network

Galectins, a family of evolutionarily conserved animal lectins, have been shown to modulate signaling processes leading to inflammation, apoptosis, immunoregulation, and angiogenesis through their ability to interact with poly-N-acetyllactosamine-enriched glycoconjugates. To date 16 human galectin carbohydrate recognition domains have been established by sequence analysis and found to be expressed in several tissues. Given the divergent functions of these lectins, it is of vital importance to understand common and differential features in order to search for specific inhibitors of individual members of the human galectin family. In this work we performed an integrated computational analysis of all individual members of the human galectin family. In the first place, we have built homology-based models for galectin-4 and -12 N-terminus, placental protein 13 (PP13) and PP13-like protein for which no experimental structural information is available. We have then performed classical molecular dynamics simulations of the whole 15 members family in free and ligand-bound states to analyze protein and protein-ligand interaction dynamics. Our results show that all galectins adopt the same fold, and the carbohydrate recognition domains are very similar with structural differences located in specific loops. These differences are reflected in the dynamics characteristics, where mobility differences translate into entropy values which significantly influence their ligand affinity. Thus, ligand selectivity appears to be modulated by subtle differences in the monosaccharide binding sites. Taken together, our results may contribute to the understanding, at a molecular level, of the structural and dynamical determinants that distinguish individual human galectins.

Role of the distal hydrogen-bonding network in regulating oxygen affinity in the truncated hemoglobin III from Campylobacter jejuni

Oxygen affinity in heme-containing proteins is determined by a number of factors, such as the nature and conformation of the distal residues that stabilize the heme bound-oxygen via hydrogen-bonding interactions. The truncated hemoglobin III from Campylobacter jejuni (Ctb) contains three potential hydrogen-bond donors in the distal site: TyrB10, TrpG8, and HisE7. Previous studies suggested that Ctb exhibits an extremely slow oxygen dissociation rate due to an interlaced hydrogen-bonding network involving the three distal residues. Here we have studied the structural and kinetic properties of the G8(WF) mutant of Ctb and employed state-of-the-art computer simulation methods to investigate the properties of the O(2) adduct of the G8(WF) mutant, with respect to those of the wild-type protein and the previously studied E7(HL) and/or B10(YF) mutants. Our data indicate that the unique oxygen binding properties of Ctb are determined by the interplay of hydrogen-bonding interactions between the heme-bound ligand and the surrounding TyrB10, TrpG8, and HisE7 residues.

Human neuroglobin functions as a redox-regulated nitrite reductase

Neuroglobin is a highly conserved hemoprotein of uncertain physiological function that evolved from a common ancestor to hemoglobin and myoglobin. It possesses a six-coordinate heme geometry with proximal and distal histidines directly bound to the heme iron, although coordination of the sixth ligand is reversible. We show that deoxygenated human neuroglobin reacts with nitrite to form nitric oxide (NO). This reaction is regulated by redox-sensitive surface thiols, cysteine 55 and 46, which regulate the fraction of the five-coordinated heme, nitrite binding, and NO formation. Replacement of the distal histidine by leucine or glutamine leads to a stable five-coordinated geometry; these neuroglobin mutants reduce nitrite to NO ∼2000 times faster than the wild type, whereas mutation of either Cys-55 or Cys-46 to alanine stabilizes the six-coordinate structure and slows the reaction. Using lentivirus expression systems, we show that the nitrite reductase activity of neuroglobin inhibits cellular respiration via NO binding to cytochrome c oxidase and confirm that the six-to-five-coordinate status of neuroglobin regulates intracellular hypoxic NO-signaling pathways. These studies suggest that neuroglobin may function as a physiological oxidative stress sensor and a post-translationally redox-regulated nitrite reductase that generates NO under six-to-five-coordinate heme pocket control. We hypothesize that the six-coordinate heme globin superfamily may subserve a function as primordial hypoxic and redox-regulated NO-signaling proteins.

Protein dynamics and ligand migration interplay as studied by computer simulation

Since proteins are dynamic systems in living organisms, the employment of methodologies contemplating this crucial characteristic results fundamental to allow revealing several aspects of their function. In this work, we present results obtained using classical mechanical atomistic simulation tools applied to understand the connection between protein dynamics and ligand migration. Firstly, we will present a review of the different sampling schemes used in the last years to obtain both ligand migration pathways and the thermodynamic information associated with the process. Secondly, we will focus on representative examples in which the schemes previously presented are employed, concerning the following: i) ligand migration, tunnels, and cavities in myoglobin and neuroglobin; ii) ligand migration in truncated hemoglobin members; iii) NO escape and conformational changes in nitrophorins; iv) ligand selectivity in catalase and hydrogenase; and v) larger ligand migration: the P450 and haloalkane dehalogenase cases. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

The first step of the dioxygenation reaction carried out by tryptophan dioxygenase and indoleamine 2,3-dioxygenase as revealed by quantum mechanical/molecular mechanical studies

Tryptophan dioxygenase (TDO) and indoleamine 2,3-dioxygenase (IDO) are two heme-containing enzymes which catalyze the conversion of L: -tryptophan to N-formylkynurenine (NFK). In mammals, TDO is mostly expressed in liver and is involved in controlling homeostatic serum tryptophan concentrations, whereas IDO is ubiquitous and is involved in modulating immune responses. Previous studies suggested that the first step of the dioxygenase reaction involves the deprotonation of the indoleamine group of the substrate by an evolutionarily conserved distal histidine residue in TDO and the heme-bound dioxygen in IDO. Here, we used classical molecular dynamics and hybrid quantum mechanical/molecular mechanical methods to evaluate the base-catalyzed mechanism. Our data suggest that the deprotonation of the indoleamine group of the substrate by either histidine in TDO or heme-bound dioxygen in IDO is not energetically favorable. Instead, the dioxygenase reaction can be initiated by a direct attack of heme-bound dioxygen on the C(2)=C(3) bond of the indole ring, leading to a protein-stabilized 2,3-alkylperoxide transition state and a ferryl epoxide intermediate, which subsequently recombine to generate NFK. The novel sequential two-step oxygen addition mechanism is fully supported by our recent resonance Raman data that allowed identification of the ferryl intermediate (Lewis-Ballester et al. in Proc Natl Acad Sci USA 106:17371-17376, 2009). The results reveal the subtle differences between the TDO and IDO reactions and highlight the importance of protein matrix in modulating stereoelectronic factors for oxygen activation and the stabilization of both transition and intermediate states.

Effects of point mutations on the structural stability of tuna myoglobins

Structural stabilities of myoglobin (Mb) from several tuna fish species significantly differ from each other, although the amino acid sequence identities are very high (>95%), suggesting that limited number of substitutions greatly affect the stability of Mb. To address this hypothesis, attempts were made to elaborate recombinant tuna Mbs with point mutations on the different residues among fish Mbs. The expression plasmid constructs were based on bigeye tuna Mb cDNA sequence, and the recombinant proteins were expressed as GST-fusion proteins in Escherichia coli. After removal of the GST segment and affinity purification, the stability of five Mb mutants, namely, A49G, T91K, K92Q, V108A, and H112Q, together with the wild type (WT) were measured, taking temperature dependency of alpha-helical content and denaturant (urea and guanidine-HCl) concentration dependency of Soret band absorbance as parameters. As a result, the mutant H112Q showed much higher stability than WT, while the structures of K92Q, T91K and A49G mutants were destabilized. No essential change in helical content was observed for V108A, but the mutant was found to be destabilized easier by the denaturants. These findings suggested that the highly conserved residues among tuna species are responsible for their stability of Mbs, but a few non-conserved residues dramatically give rise to the differences in stability of Mbs among species.