Very high resolution structure of a trematode hemoglobin displaying a TyrB10-TyrE7 heme distal residue pair and high oxygen affinity

Hypoxic and nitrosative stress conditions modulate expression of myoglobin genes in a carcinogenic hepatobiliary trematode, Clonorchis sinensis

Despite recent evidence suggesting that adult trematodes require oxygen for the generation of bioenergy and eggshells, information on the molecular mechanism by which the parasites acquire oxygen remains largely elusive. In this study, the structural and expressional features of globin genes identified in Clonorchis sinensis, a carcinogenic trematode parasite that invades the hypoxic biliary tracts of mammalian hosts, were investigated to gain insight into the molecules that enable oxygen metabolism. The number of globin paralogs substantially differed among parasitic platyhelminths, ranging from one to five genes, and the C. sinensis genome encoded at least five globin genes. The expression of these Clonorchis genes, named CsMb (CsMb1—CsMb3), CsNgb, and CsGbX, according to their preferential similarity patterns toward respective globin subfamilies, exponentially increased in the worms coinciding with their sexual maturation, after being downregulated in early juveniles compared to those in metacercariae. The CsMb1 protein was detected throughout the parenchymal region of adult worms as well as in excretory-secretory products, whereas the other proteins were localized exclusively in the sexual organs and intrauterine eggs. Stimuli generated by exogenous oxygen, nitric oxide (NO), and nitrite as well as co-incubation with human cholangiocytes variously affected globin gene expression in live C. sinensis adults. Together with the specific histological distributions, these hypoxia-induced patterns may suggest that oxygen molecules transported by CsMb1 from host environments are provided to cells in the parenchyma and intrauterine eggs/sex organs of the worms for energy metabolism and/or, more importantly, eggshell formation by CsMb1 and CsMb3, respectively. Other globin homologs are likely to perform non-respiratory functions. Based on the responsive expression profile against nitrosative stress, an oxygenated form of secreted CsMb1 is suggested to play a pivotal role in parasite survival by scavenging NO generated by host immune cells via its NO dioxygenase activity.

Trematode parasites that invade mammalian tissues have long been believed to produce bioenergy via anaerobic respiration in their definitive hosts. However, recent studies have revealed that these parasites require considerable amounts of oxygen for the generation of hard eggshells during sexual reproduction as well as energy metabolism. Despite these findings, information on the biological mechanisms and relevant molecules responsible for oxygen uptake in the host environment remains largely elusive. Clonorchis sinensis is a carcinogenic trematode parasite that causes clonorchiasis in humans by infecting the bile ducts. Here, we investigated globin genes/proteins in the liver fluke. The genome of C. sinensis encoded at least five globin paralogs (CsMb1, CsMb2, CsMb3, CsNgb, and CsGbX). Temporal expression of these globin genes coincided with the sexual maturation of C. sinensis. Based on the histological localities and induction profiles upon hypoxia, it could be postulated that the oxygen molecules transported by CsMb1 from host environments are provided to cells in the parenchyma and intrauterine eggs/sex organs of the worms by CsMb1 and CsMb3, respectively, for energy metabolism and eggshell formation. Other globin homologs were likely to perform non-respiratory functions. In addition, the oxygenated form of secreted CsMb1 seemed to participate in the scavenging of nitric oxide generated by host immune cells via its nitric oxide dioxygenase activity to increase the survival of the parasite.

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.

A molecule for all seasons: The heme

If life without heme-Fe were at all possible, it would definitely be different. Indeed this complex and versatile iron-porphyrin macrocycle upon binding to different “globins” yields hemeproteins crucial to sustain a variety of vital functions, generally classified, for convenience, in a limited number of functional families. Over-and-above the array of functions briefly outlined below, the spectacular progress in molecular genetics seen over the last 30 years led to the discovery of many hitherto unknown novel hemeproteins in prokaryotes and eukaryotes. Here, we highlight a few basic aspects of the chemistry of the hemeprotein universe, in particular those that are relevant to the control of heme-Fe reactivity and specialization, as sculpted by a variety of interactions with the protein moiety.

Kinetic properties and heme pocket structure of two domains of the polymeric hemoglobin of Artemia in comparison with the native molecule

In this project, we studied some physicochemical properties of two different globin domains of the polymeric hemoglobin of the brine shrimp Artemia salina and compared them with those of the native molecule. Two domains (AsHbC1D1 and AsHbC1D5) were cloned and expressed in BL21(DE3)pLysS strain of Escherichia coli. The recombinant proteins as well as the native hemoglobin (AfHb) were purified from bacteria and frozen Artemia, respectively by standard chromatographic methods and assessed by SDS-PAGE. The heme environment of these proteins was studied by optical spectroscopy and ligand-binding kinetics (e.g. CO association and O2 binding affinity) were measured for the two recombinant proteins and the native hemoglobin. This indicates that the CO association rate for AsHbC1D1 is higher than that of AsHbC1D5 and AfHb, while the calculated P50 value for AsHbC1D1 is lower than that of AsHbC1D5 and AfHb. The geminate and bimolecular rebinding parameters indicate a significant difference between both domains. Moreover, EPR results showed that the heme pocket in AfHb is in a more closed conformation than the heme pocket in myoglobin. Finally, the reduction potential of -0.13V versus the standard hydrogen electrode was determined for AfHb by direct electrochemical measurements. It is about 0.06V higher than the potential of the single domain AsHbC1D5. This work shows that each domain in the hemoglobin of Artemia has different characteristics of ligand binding.

A globin domain in a neuronal transmembrane receptor of Caenorhabditis elegans and Ascaris suum: molecular modeling and functional properties

We report the structural and biochemical characterization of GLB-33, a putative neuropeptide receptor that is exclusively expressed in the nervous system of the nematode Caenorhabditis elegans. This unique chimeric protein is composed of a 7-transmembrane domain (7TM), GLB-33 7TM, typical of a G-protein-coupled receptor, and of a globin domain (GD), GLB-33 GD. Comprehensive sequence similarity searches in the genome of the parasitic nematode, Ascaris suum, revealed a chimeric protein that is similar to a Phe-Met-Arg-Phe-amide neuropeptide receptor. The three-dimensional structures of the separate domains of both species and of the full-length proteins were modeled. The 7TM domains of both proteins appeared very similar, but the globin domain of the A. suum receptor surprisingly seemed to lack several helices, suggesting a novel truncated globin fold. The globin domain of C. elegans GLB-33, however, was very similar to a genuine myoglobin-type molecule. Spectroscopic analysis of the recombinant GLB-33 GD showed that the heme is pentacoordinate when ferrous and in the hydroxide-ligated form when ferric, even at neutral pH. Flash-photolysis experiments showed overall fast biphasic CO rebinding kinetics. In its ferrous deoxy form, GLB-33 GD is capable of reversibly binding O2 with a very high affinity and of reducing nitrite to nitric oxide faster than other globins. Collectively, these properties suggest that the globin domain of GLB-33 may serve as a highly sensitive oxygen sensor and/or as a nitrite reductase. Both properties are potentially able to modulate the neuropeptide sensitivity of the neuronal transmembrane receptor.

Peroxidase activity and involvement in the oxidative stress response of roseobacter denitrificans truncated hemoglobin

Roseobacter denitrificans is a member of the widespread marine Roseobacter genus. We report the first characterization of a truncated hemoglobin from R. denitrificans (Rd. trHb) that was purified in the heme-bound form from heterologous expression of the protein in Escherichia coli. Rd. trHb exhibits predominantly alpha-helical secondary structure and absorbs light at 412, 538 and 572 nm. The phylogenetic classification suggests that Rd. trHb falls into group II trHbs, whereas sequence alignments indicate that it shares certain important heme pocket residues with group I trHbs in addition to those of group II trHbs. The resonance Raman spectra indicate that the isolated Rd. trHb contains a ferric heme that is mostly 6-coordinate low-spin and that the heme of the ferrous form displays a mixture of 5- and 6-coordinate states. Two Fe-His stretching modes were detected, notably one at 248 cm^-1, which has been reported in peroxidases and some flavohemoglobins that contain an Fe-His-Asp (or Glu) catalytic triad, but was never reported before in a trHb. We show that Rd. trHb exhibits a significant peroxidase activity with a (k_cat/K_m) value three orders of magnitude higher than that of bovine Hb and only one order lower than that of horseradish peroxidase. This enzymatic activity is pH-dependent with a pK_a value ~6.8. Homology modeling suggests that residues known to be important for interactions with heme-bound ligands in group II trHbs from Mycobacterium tuberculosis and Bacillus subtilis are pointing toward to heme in Rd. trHb. Genomic organization and gene expression profiles imply possible functions for detoxification of reactive oxygen and nitrogen species in vivo. Altogether, Rd. trHb exhibits some distinctive features and appears equipped to help the bacterium to cope with reactive oxygen/nitrogen species and/or to operate redox biochemistry.

Refinement by shifting secondary structure elements improves sequence alignments

Constructing a model of a query protein based on its alignment to a homolog with experimentally determined spatial structure (the template) is still the most reliable approach to structure prediction. Alignment errors are the main bottleneck for homology modeling when the query is distantly related to the template. Alignment methods often misalign secondary structural elements by a few residues. Therefore, better alignment solutions can be found within a limited set of local shifts of secondary structures. We present a refinement method to improve pairwise sequence alignments by evaluating alignment variants generated by local shifts of template-defined secondary structures. Our method SFESA is based on a novel scoring function that combines the profile-based sequence score and the structure score derived from residue contacts in a template. Such a combined score frequently selects a better alignment variant among a set of candidate alignments generated by local shifts and leads to overall increase in alignment accuracy. Evaluation of several benchmarks shows that our refinement method significantly improves alignments made by automatic methods such as PROMALS, HHpred and CNFpred. The web server is available at http://prodata.swmed.edu/sfesa.

Molecular characterization of Clonorchis sinensis secretory myoglobin: delineating its role in anti-oxidative survival

Background

Clonorchiasis is a globally important, neglected food-borne disease caused by Clonorchis sinensis (C. sinensis), and it is highly related to cholangiocarcinoma and hepatocellular carcinoma. Increased molecular evidence has strongly suggested that the adult worm of C. sinensis continuously releases excretory-secretory proteins (ESPs), which play important roles in the parasite-host interactions, to establish successful infection and ensure its own survival. Myoglobin, a hemoprotein, is present in high concentrations in trematodes and ESPs. To further understand the biological function of CsMb and its putative roles in the interactions of C. sinensis with its host, we explored the molecular characterization of CsMb in this paper.

Methods

We expressed CsMb and its mutants in E. coli BL21 and identified its molecular characteristics using bioinformatics analysis and experimental approaches. Reverse transcription PCR analysis was used to measure myoglobin transcripts of C. sinensis with different culture conditions. The peroxidase activity of CsMb was confirmed by spectrophotometry. We co-cultured RAW264.7 cells with recombinant CsMb (rCsMb), and we then measured the production of hydrogen peroxide (H₂O₂) and nitric oxide (NO) in addition to the mRNA levels of inducible nitric oxide synthase (iNOS), Cu-Zn superoxide dismutase (SOD1) and Mn superoxide dismutase (SOD2) in activated RAW264.7 cells.

Results

In the in vitro culture of adult worms, the transcripts of CsMb increased with the increase of oxygen content. Oxidative stress conditions induced by H₂O₂ increased the levels of CsMb transcripts in a dose-dependent manner. Furthermore, CsMb catalyzed oxidation reactions in the presence of H₂O₂, and amino acid 34 of CsMb played an essential role in its reaction with H₂O₂. In addition, CsMb significantly reduced H₂O₂ and NO levels in LPS-activated macrophages, and CsMb downregulated iNOS and SOD expression in activated macrophages.

Conclusion

The present study is the first to investigate the peroxidase activity of CsMb. This investigation suggested that C. sinensis may decrease the redox activation of macrophages by CsMb expression to evade host immune responses. These studies contribute to a better understanding of the role of CsMb in the molecular mechanisms involved in ROS detoxification by C. sinensis.

Adult Opisthorchis felineus major protein fractions deduced from transcripts: comparison with liver flukes Opisthorchis viverrini and Clonorchis sinensis

The epidemiologically important liver flukes Opisthorchis felineus, Opisthorchis viverrini, and Clonorchis sinensis are of interest to health professionals, epidemiologists, pharmacologists, and molecular biologists. Recently the transcriptomes of the latter two species were intensively investigated. However our knowledge on molecular biology of O. felineus is scarce. We report the first results of the O. felineus transcriptome analysis. We isolated and annotated a total of 2560 expressed sequence tag (EST) sequences from adult O. felineus (deposited within the database of expressed sequence tags (dbEST), under accession numbers GenBank: JK624271-JK626790, JK006511-JK006547, JK649790-JK649792). Clustering and analysis resulted in the detection of 267 contigs. Of the protein sequences deduced from these, 82% had homologs in the NCBI (nr) protein database and 63% contained conserved domains, allowing the functions to be interpreted using the Gene Ontology terms. Comprehensive analysis of Opisthorchiidae- and Trematoda-specific substitutions within amino acid sequences deduced for the proteins myoglobin, vitelline precursor protein, cathepsin F, and 28kDa glutathione transferase was carried out. The gene set of the 32 ribosomal proteins for the three Opisthorchiidae species with the addition of available Schistosoma and Fasciola orthologs was created and is provided in the supplementary. The orthologous gene set created was used for inferring phylogeny within the Trematoda with special attention to interrelations within the Opisthorchiidae. The phylogenetic analysis revealed a closer relationship between C. sinensis and O. viverrini and some divergence of O. felineus from either O. viverrini or C. sinensis.

Ligand-rebinding kinetics of 2/2 hemoglobin from the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125

Kinetic studies were performed on ligand rebinding to a cold-adapted globin of the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 (Ph-2/2HbO). This 2/2 hemoglobin displays a rapid spectroscopic phase that is independent of CO concentration, followed by the standard bimolecular recombination. While the geminate recombination usually occurs on a ns timescale, Ph-2/2HbO displays a component of about 1μs that accounts for half of the geminate phase at 8°C, indicative of a relatively slow internal ligand binding. The O2 binding kinetics were measured in competition with CO to allow a short-time exposure of the deoxy hemes to O2 before CO replacement. Indeed Ph-2/2HbO is readily oxidised in the presence of O2, probably due to a superoxide character of the FeO2 bond induced by of a hydrogen-bond donor amino-acid residue. Upon O2 release or iron oxidation a distal residue (probably Tyr) is able to reversibly bind to the heme and as such to compete for binding with an external ligand. The transient hexacoordinated ferrous His-Fe-Tyr conformation after O2 dissociation could initiate the electron transfer from the iron toward its final acceptor, molecular O2 under our conditions. The hexacoordination via the distal Tyr is only partial, indicating a weak interaction between Tyr and the heme under atmospheric pressure. Hydrostatic high pressure enhances the hexacoordination indicating a flexible globin that allows structural changes. The O2 binding affinity for Ph-2/2HbO, poorly affected by the competition with Tyr, is about 1Torr at 8°C, pH7.0, which is compatible for an in vivo O2 binding function; however, this globin is more likely involved in a redox reaction associating diatomic ligands and their derived oxidative species. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Conservation of the three-dimensional structure in non-homologous or unrelated proteins

In this review, we examine examples of conservation of protein structural motifs in unrelated or non-homologous proteins. For this, we have selected three DNA-binding motifs: the histone fold, the helix-turn-helix motif, and the zinc finger, as well as the globin-like fold. We show that indeed similar structures exist in unrelated proteins, strengthening the concept that three-dimensional conservation might be more important than the primary amino acid sequence.

Molecular dynamics free energy calculations to assess the possibility of water existence in protein nonpolar cavities

Are protein nonpolar cavities filled with water molecules? Although many experimental and theoretical investigations have been done, particularly for the nonpolar cavity of IL-1 beta, the results are still conflicting. To study this problem from the thermodynamic point of view, we calculated hydration free energies of four protein nonpolar cavities by means of the molecular dynamics thermodynamic integration method. In addition to the IL-1 beta cavity (69 A(3)), we selected the three largest nonpolar cavities of AvrPphB (81 A(3)), Trp repressor (87 A(3)), and hemoglobin (108 A(3)) from the structural database, in view of the simulation result from another study that showed larger nonpolar cavities are more likely to be hydrated. The calculations were performed with flexible and rigid protein models. The calculated free energy changes were all positive; hydration of the nonpolar cavities was energetically unfavorable for all four cases. Because hydration of smaller cavities should happen more rarely, we conclude that existing protein nonpolar cavities are not likely to be hydrated. Although a possibility remains for much larger nonpolar cavities, such cases are not found experimentally. We present a hypothesis to explain this: hydrated nonpolar cavities are quite unstable and the conformation could not be maintained.

Rapid comparison of properties on protein surface

The mapping of physicochemical characteristics onto the surface of a protein provides crucial insights into its function and evolution. This information can be further used in the characterization and identification of similarities within protein surface regions. We propose a novel method which quantitatively compares global and local properties on the protein surface. We have tested the method on comparison of electrostatic potentials and hydrophobicity. The method is based on 3D Zernike descriptors, which provides a compact representation of a given property defined on a protein surface. Compactness and rotational invariance of this descriptor enable fast comparison suitable for database searches. The usefulness of this method is exemplified by studying several protein families including globins, thermophilic and mesophilic proteins, and active sites of TIM beta/alpha barrel proteins. In all the cases studied, the descriptor is able to cluster proteins into functionally relevant groups. The proposed approach can also be easily extended to other surface properties. This protein surface-based approach will add a new way of viewing and comparing proteins to conventional methods, which compare proteins in terms of their primary sequence or tertiary structure.

The hemoglobins of the trematodes Fasciola hepatica and Paramphistomum epiclitum: a molecular biological, physico-chemical, kinetic, and vaccination study

The trematode Fasciola hepatica (Fa.he.) is a common parasite of human and livestock. The hemoglobin (Hb) of Fa.he., a potential immunogen, was chosen for characterization in the search for an effective vaccine. Characterization of trematode Hbs show that they are intracellular single-domain globins with the following remarkable features: (1) Fa.he. expresses two Hb isoforms that differ at two amino acid sites (F1: 119Y/123Q; F2: 119F/123L). Both isoforms are monoacetylated at their N-termini; (2) the genes coding for Fa.he. and Paramphistomum epiclitum (Pa.ep.) Hbs are interrupted by two introns at the conserved positions B12.2 and G7.0.; (3) UV/VIS and resonance Raman spectroscopy identify the recombinant Fa.he. HbF2 as a pentacoordinated high-spin ferrous Hb; (4) electron paramagnetic resonance spectroscopy of cyano-met Fa.he. HbF2 proves that the endogenously bound imidazole has no imidazolate character; (5) the major structural determinants of the globin fold are present, they contain a TyrB10/TyrE7 residue pair on the distal side. Although such distal-site pair is a signature for high oxygen affinity, as shown for Pa.ep. Hb, the oxygen-binding rate parameters for Fa.he. Hb are intermediate between those of myoglobin and those of other trematode Hbs; (6) the three-dimensional structure of recombinant Fa.he. HbF2 from this study closely resembles the three-dimensional structure of Pa.ep. determined earlier. The set of distal-site polar interactions observed in Pa.ep. Hb is matched with small but significant structural adjustments; (7) despite the potential immunogenic character of the fluke Hb, vaccination of calves with recombinant Fa.he. HbF2 failed to promote protection against parasitic infection.

Exploring the molecular basis of heme coordination in human neuroglobin

Neuroglobin (Ngb), a recently discovered ancient heme protein, presents the typical globin fold and is around 20% identical to myoglobin (Mb). In contrast with Mb, however, its heme is hexacoordinated (6c). It is expressed in the nervous system and has been the subject of numerous investigations in the last years, but its function is still unclear. The proposed roles include oxygen transport, reactive oxygen species (ROS) detoxification, hypoxia protection, and redox state sensing. All proposed functions require distal histidine dissociation from the heme to yield a reactive iron. With the aim of understanding the 6c to 5c transition, we have performed molecular dynamics simulations for ferrous Ngb in the 6c, 5c, and oxy states. We also computed free energy profiles associated with the transition employing an advanced sampling technique. Finally, we studied the effect of the redox state of CysCD7 and CysD5, which are known to form a disulfide bridge. Our results show that protein oxidation promotes a stabilization of the pentacoordinated species, thus favoring the protein to adopt the more reactive state and supporting the existence of a molecular mechanism whereby O2 would be released under hypoxic conditions, thereby suggesting an O(2) storage function for Ngb. Taken together, our results provide structural information not available experimentally which may shed light on the protein proposed functions, particularly as a redox sensor.

O2 migration pathways are not conserved across proteins of a similar fold

Recent advances in computational biology have made it possible to map the complete network and energy profile of gas migration pathways inside proteins. Although networks of O(2) pathways have already been characterized for a small number of proteins, the general properties and locations of these pathways have not been previously compared between proteins. In this study, maps of the O(2) pathways inside 12 monomeric globins were computed. It is found that, despite the conserved tertiary structure fold of the studied globins, the shape and topology of O(2) pathway networks exhibit a large variability between different globins, except when two globins are nearly identical. The locations of the O(2) pathways are, however, found to be correlated with the location of large hydrophobic residues, and a similar correlation is observed in two unrelated protein families: monomeric globins and copper-containing amine oxidases. The results have implications for the evolution of gas pathways in proteins and for protein engineering applications involving modifications of these pathways.

Molecular characterization of two myoglobins of Paragonimus westermani

Myoglobins (Mbs), globin proteins, are present in high concentrations in trematodes. In Paragonimus westermani, 2 cDNAs were found to encode Mbs. The first clone, Pwmyo1, codes a total of 149 amino acids with a calculated mass of 16.6 kDa. The second, Pwmyo2, encodes a 146-amino acid protein with a calculated mass of 16.2 kDa. The predicted secondary structures showed the presence of 8 helices, which is the basic characteristic of Mbs. Sequence alignment revealed a high homology with the other trematode Mbs. The 2 clones contained the characteristic tyrosyl residues at helical positions B10 and distal E7, which are substitutions that have been previously shown to contribute to the high oxygen affinity of Mbs. Polyclonal antibodies against the recombinant Mbs were raised with no cross-reactivity observed. Immunolocalization revealed the proteins to be distributed generally throughout the parenchymal tissues, but absent from the tegument and reproductive organs. The cell mass of the eggs of the worm stained positive to Pwmyo2 but not Pwmyo1, suggesting the stage-specific expression of these Mbs.

Crystallographic evidence for dioxygen interactions with iron proteins

The interaction of dioxygen with iron plays a key role in many important biological processes, such as dioxygen transport in the bloodstream and the reduction of dioxygen by iron in respiration. However, the catalytic mechanisms employed, for example in ligand oxidation, are not fully understood at the current time despite intensive biochemical, spectroscopic and structural studies. This review outlines the structural evidence obtained by X-ray crystallographic methods for the nature of the interactions between dioxygen and the metal in iron-containing proteins. Proteins involved in iron transport or electron transfer are not included.

Ligand Selectivity of Soluble Guanylyl Cyclase

Although soluble guanylyl cyclase (sGC) functions in an environment in which O(2), NO, and CO are potential ligands for its heme moiety, the enzyme displays a high affinity for only its physiological ligand, NO, but has a limited affinity for CO and no affinity for O(2). Recent studies of a truncated version of the sGC beta(1)-subunit containing the heme-binding domain (Boon, E. M., Huang, S H., and Marletta, M. A. (2005) Nat. Chem. Biol., 1, 53-59) showed that introduction of the hydrogen-bonding tyrosine into the distal heme pocket changes the ligand specificity of the heme moiety and results in an oxygen-binding sGC. The hypothesis that the absence of hydrogen-bonding residues in the distal heme pocket is sufficient to provide oxygen discrimination by sGC was put forward. We tested this hypothesis in a context of a complete sGC heterodimer containing both the intact alpha(1)- and beta(1)-subunits. We found that the I145Y substitution in the full-length beta-subunit of the sGC heterodimer did not produce an oxygen-binding enzyme. However, this substitution impeded the association of NO and destabilized the NO.heme complex. The tyrosine in the distal heme pocket also impeded both the binding and dissociation of the CO ligand. We propose that the mechanism of oxygen exclusion by sGC not only involves the lack of hydrogen bonding in the distal heme pocket, but also depends on structural elements from other domains of sGC.

Multiple active site conformers in the carbon monoxide complexes of trematode hemoglobins

Sequence alignment of hemoglobins of the trematodes Paramphistomum epiclitum and Gastrothylax crumenifer with myoglobin suggests the presence of an unusual active site structure in which two tyrosine residues occupy the E7 and B10 helical positions. In the crystal structure of P. epiclitum hemoglobin, such an E7-B10 tyrosine pair at the putative helical positions has been observed, although the E7 Tyr is displaced toward CD region of the polypeptide. Resonance Raman data on both P. epiclitum and G. crumenifer hemoglobins show that interactions of heme-bound ligands with neighboring amino acid residues are unusual. Multiple conformers in the CO complex, termed the C, O, and N conformers, are observed. The conformers are separated by a large difference (approximately 60 cm(-1)) in the frequencies of their Fe-CO stretching modes. In the C conformer the Fe-CO stretching frequency is very high, 539 and 535 cm(-1), for the P. epiclitum and G. crumenifer hemoglobins, respectively. The Fe-CO stretching of the N conformer appears at an unusually low frequency, 479 and 476 cm(-1), respectively, for the two globins. A population of an O conformer is seen in both hemoglobins, at 496 and 492 cm(-1), respectively. The C conformer is stabilized by a strong polar interaction of the CO with the distal B10 tyrosine residue. The O conformer is similar to the ones typically seen in mutant myoglobins in which there are no strong interactions between the CO and residues in the distal pocket. The N conformer possesses an unusual configuration in which a negatively charged group, assigned as the oxygen atom of the B10 Tyr side chain, interacts with the CO. In this conformer, the B10 Tyr assumes an alternative conformation consistent with one of the conformers seen the crystal structure. Implications of the multiple configurations on the ligand kinetics are discussed.

Two distinct heme distal site states define Cerebratulus lacteus mini-hemoglobin oxygen affinity

The nerve tissue hemoglobin of Cerebratulus lacteus (CerHb) is the smallest naturally occurring known hemoglobin. Stabilization of the diatomic bound species (e.g., O(2)) is achieved through a network of hydrogen bonds based on three key residues TyrB10, GlnE7, and ThrE11. The first two residues are typically associated in hemoglobins with enhanced O(2) affinity, related to hydrogen bond stabilization of the heme-bound O(2) resulting in a decrease of the ligand dissociation rates. In contrast to the above observations, the affinity of CerHb for O(2) is only moderate, and the rate of O(2) dissociation is unexpectedly high. To gain insight on the diverse molecular mechanisms controlling ligand affinities, we have analyzed w.t. CerHb and its ThrE11-->Val mutant by means of joint molecular dynamics and quantum mechanics simulation techniques, complementing recent site-directed mutagenesis experiments. Our results suggest that the observed O(2) dissociation rates can only be explained through a dynamic equilibrium between high and low affinity states of the w.t. CerHb heme distal site.

Invertebrate hemoglobins and nitric oxide: how heme pocket structure controls reactivity

Hemoglobins (Hbs), generally defined as 5 or 6 coordinate heme proteins whose primary function is oxygen transport, are now recognized to occur in virtually all phyla of living organisms. Historically, study of their function focused on oxygen as a reversibly bound ligand of the ferrous form of the protein. Other diatomic ligands like carbon monoxide and nitric oxide were considered "non-physiological" but useful probes of structure-function relationships in Hbs. This investigatory landscape changed dramatically in the 1980s when nitric oxide was discovered to activate a heme protein, cyclic guanylate cyclase. Later, its activation was likened to Perutz' description of Hb's allosteric properties being triggered by a ligand-dependent "out-of-plane/into-plane" movement of the heme iron. In 1996, a functional role for nitric oxide in human and mammalian Hbs was demonstrated and since that time, the interest in NO as a physiologically relevant Hb ligand has greatly increased. Concomitantly, non-oxygen binding properties of Hbs have challenged the view that Hbs arose for their oxygen storage and transport properties. In this focused review we discuss some invertebrate Hbs' functionally significant reactions with nitric oxide and how strategic positioning of a few residues in the heme pocket plays an large role in the interplay of diatomic ligands to ferrous and ferric heme iron in these proteins.

pH effects on the binding of oxygen to non-vertebrate monomeric hemoglobins. A linked function model

Monomeric invertabrate hemoglobins with high oxygen affinity usually contain a tyrosine in the distal region of the heme. This feature has stimulated investigations revealing that one of the properties resulting from the presence of the distal tyrosines is a decreased off rate on the binding of oxygen, thus developing the high affinity. Despite that fact that the pK value of the tyrosine group differs significantly from the groups it replaces little attention has been paid to the pH dependence of the binding of oxygen to the high affinity hemoglobins. Such a pH dependence has been reported on two of the monomeric hemoglobins with relatively low oxygen affinity and one monomeric hemoglobin of intermediate affinity. The pH data of these hemoglobins has been analysed with a linked function model involving the hydrogen ion. pK values required for the low-affinity hemoglobins vary from 4.5 to 7.5. When applied to the high-affinity hemoglobins, the linked function model provides reasonable values for the binding parameters. These pK values vary from 3.0 to 9.0.

Tyrosine B10 inhibits stabilization of bound carbon monoxide and oxygen in soybean leghemoglobin

Detailed comparisons of the carbon monoxide FTIR spectra and ligand-binding properties of a library of E7, E11, and B10 mutants indicate significant differences in the role of electrostatic interactions in the distal pockets of wild-type sperm whale myoglobin and soybean leghemoglobin. In myoglobin, strong hydrogen bonds from several closely related conformations of the distal histidine (His(E7)) side chain preferentially stabilize bound oxygen. In leghemoglobin, the imidazole side chain of His(E7) is confined to a single conformation, which only weakly hydrogen bonds to bound ligands. The phenol side chain of Tyr(B10) appears to "fix" the position of His(E7), probably by donating a hydrogen bond to the Ndelta atom of the imidazole side chain. The proximal pocket of leghemoglobin is designed to favor strong coordination bonds between the heme iron and axial ligands. Thus, high oxygen affinity in leghemoglobin is established by a favorable staggered geometry of the proximal histidine. The interaction between His(E7) and Tyr(B10) prevents overstabilization of bound oxygen. If hydrogen bonding from His(E7) were as strong as it is in mammalian myoglobin, the resultant ultrahigh affinity of leghemoglobin would prevent oxygen transport in root nodules.

Thr-E11 regulates O2 affinity in Cerebratulus lacteus mini-hemoglobin

The mini-hemoglobin from Cerebratulus lacteus (CerHb) belongs to a class of globins containing the polar Tyr-B10/Gln-E7 amino acid pair that normally causes low rates of O2 dissociation and ultra-high O2 affinity, which suggest O2 sensing or NO scavenging functions. CerHb, however, has high rates of O2 dissociation (kO2 = 200-600 s(-1)) and moderate O2 affinity (KO2) approximately 1 microm(-1)) as a result of a third polar amino acid in its active site, Thr-E11. When Thr-E11 is replaced by Val, kO2 decreases 1000-fold and KO2 increases 130-fold at pH 7.0, 20 degrees C. The mutation also shifts the stretching frequencies of both heme-bound and photodissociated CO, indicating marked changes of the electrostatic field at the active site. The crystal structure of Thr-E11 --> Val CerHbO2 at 1.70 A resolution is almost identical to that of the wild-type protein (root mean square deviation of 0.12 A). The dramatic functional and spectral effects of the Thr-E11 --> Val mutation are due exclusively to changes in the hydrogen bonding network in the active site. Replacing Thr-E11 with Val "frees" the Tyr-B10 hydroxyl group to rotate toward and donate a strong hydrogen bond to the heme-bound ligand, causing a selective increase in O2 affinity, a decrease of the rate coefficient for O2 dissociation, a 40 cm(-1) decrease in nuCO of heme-bound CO, and an increase in ligand migration toward more remote intermediate sites.

Kinetic modulation in carbonmonoxy derivatives of truncated hemoglobins: the role of distal heme pocket residues and extended apolar tunnel

Truncated hemoglobins (trHbs), are a distinct and newly characterized class of small myoglobin-like proteins that are widely distributed in bacteria, unicellular eukaryotes, and higher plants. Notable and distinctive features associated with trHbs include a hydrogen-bonding network within the distal heme pocket and a long apolar tunnel linking the external solvent to the distal heme pocket. The present work compares the geminate and solvent phase rebinding kinetics from two trHbs, one from the ciliated protozoan Paramecium caudatum (P-trHb) and the other from the green alga Chlamydomonas eugametos (C-trHb). Unusual kinetic patterns are observed including indications of ultrafast (picosecond) geminate rebinding of CO to C-trHb, very fast solvent phase rebinding of CO for both trHbs, time-dependent biphasic CO rebinding kinetics for P-trHb at low CO partial pressures, and for P-trHb, an increase in the geminate yield from a few percent to nearly 100% under high viscosity conditions. Species-specific differences in both the 8-ns photodissociation quantum yield and the rebinding kinetics, point to a pivotal functional role for the E11 residue. The response of the rebinding kinetics to temperature, ligand concentration, and viscosity (glycerol, trehalose) and the viscosity-dependent changes in the resonance Raman spectrum of the liganded photoproduct, together implicate both the apolar tunnel and the static and dynamic properties of the hydrogen-bonding network within the distal heme pocket in generating the unusual kinetic patterns observed for these trHbs.

1H NMR study of the molecular structure and magnetic properties of the active site for the cyanomet complex of O2-avid hemoglobin from the trematode Paramphistomum epiclitum

The solution molecular and electronic structures of the active site in the extremely O2-avid hemoglobin from the trematode Paramphistomum epiclitum have been investigated by 1H NMR on the cyanomet form in order to elucidate the distal hydrogen-bonding to a ligated H-bond acceptor ligand. Comparison of the strengths of dipolar interactions in solution with the alternate crystal structures of methemoglobin establish that the solution structure of wild-type Hb more closely resembles the crystal structure of the recombinant wild-type than the true wild-type met-hemoglobin. The distal Tyr66(E7) is found oriented out of the heme pocket in solution as found in both crystal structures. Analysis of dipolar contacts, dipolar shift and paramagnetic relaxation establishes that the Tyr32(B10) hydrogen proton adopts an orientation that allows it to make a strong H-bond to the bound cyanide. The observation of a significant isotope effect on the heme methyl contact shifts confirms a strong contact between the Tyr32(B10) OH and the ligated cyanide. The quantitative determination of the orientation and anisotropies of the paramagnetic susceptibility tensor reveal that the cyanide is tilted approximately 10 degrees from the heme normal so as to avoid van der Waals overlap with the Tyr32(B10) Oeta. The pattern of heme contact shifts with large low-field shifts for 7-CH3 and 18-CH3 is shown to arise not from the 180 degrees rotation about the alpha-gamma-meso axis, but due to the approximately 45 degrees rotation of the axial His imidazole ring, relative to that in mammalian globins.

A TyrCD1/TrpG8 hydrogen bond network and a TyrB10TyrCD1 covalent link shape the heme distal site of Mycobacterium tuberculosis hemoglobin O

Truncated hemoglobins (Hbs) are small hemoproteins, identified in microorganisms and in some plants, forming a separate cluster within the Hb superfamily. Two distantly related truncated Hbs, trHbN and trHbO, are expressed at different developmental stages in Mycobacterium tuberculosis. Sequence analysis shows that the two proteins share 18% amino acid identities and belong to different groups within the truncated Hb cluster. Although a specific defense role against nitrosative stress has been ascribed to trHbN (expressed during the Mycobacterium stationary phase), no clear functions have been recognized for trHbO, which is expressed throughout the Mycobacterium growth phase. The 2.1-A crystal structure of M. tuberculosis cyano-met trHbO shows that the protein assembles in a compact dodecamer. Six of the dodecamer subunits are characterized by a double conformation for their CD regions and, most notably, by a covalent bond linking the phenolic O atom of TyrB10 to the aromatic ring of TyrCD1, in the heme distal cavity. All 12 subunits display a cyanide ion bound to the heme Fe atom, stabilized by a tight hydrogen-bonded network based on the (globin very rare) TyrCD1 and TrpG8 residues. The small apolar AlaE7 residue leaves room for ligand access to the heme distal site through the conventional "E7 path," as proposed for myoglobin. Different from trHbN, where a 20-A protein matrix tunnel is held to sustain ligand diffusion to an otherwise inaccessible heme distal site, the topologically related region in trHbO hosts two protein matrix cavities.

Distal heme pocket regulation of ligand binding and stability in soybean leghemoglobin

Leghemoglobins facilitate diffusion of oxygen through root tissue to a bacterial terminal oxidase in much the same way that myoglobin transports oxygen from blood to muscle cell mitochondria. Leghemoglobin serves an additional role as an oxygen scavenger to prevent inhibition of nitrogen fixation. For this purpose, the oxygen affinity of soybean leghemoglobin is 20-fold greater than myoglobin, resulting from an 8-fold faster association rate constant combined with a 3-fold slower dissociation rate constant. Although the biochemical mechanism used by myoglobin to bind oxygen has been described in elegant detail, an explanation for the difference in affinity between these two structurally similar proteins is not obvious. The present work demonstrates that, despite their similar structures, leghemoglobin uses methods different from myoglobin to regulate ligand affinity. Oxygen and carbon monoxide binding to a comprehensive set of leghemoglobin distal heme pocket mutant proteins in comparison to their myoglobin counterparts has revealed some of these mechanisms. The "distal histidine" provides a crucial hydrogen bond to stabilize oxygen in myoglobin but has little effect on bound oxygen in leghemoglobin and is retained mainly for reasons of protein stability and prevention of heme loss. Furthermore, soybean leghemoglobin uses an unusual combination of HisE7 and TyrB10 to sustain a weak stabilizing interaction with bound oxygen. Thus, the leghemoglobin distal heme pocket provides a much lower barrier to oxygen association than occurs in myoglobin and oxygen dissociation is regulated from the proximal heme pocket.

[Hemoglobin, from microorganisms to man: a single structural motif, multiple functions]

Haemoglobins from unicellular organisms, plants or animals, share a common structure, which results from the folding, around the heme group, of a polypeptide chain made from 6-8 helices. Nowadays, deciphering the genome of several species allows one to draw the evolutionary tree of this protein going back to 1800 millions of years, at a time when oxygen began to accumulate in the atmosphere. This permits to follow the evolution of the ancestral gene and of its product. It is likely that, only in complex multicellular species, transport and storage of oxygen became the main physiological function of this molecule. In addition, in unicellular organisms and small invertebrates, it is likely that the main function of this protein was to protect the organism from the toxic effect of O2, CO and NO*. The very high oxygen affinity of these molecules, leading them to act rather as a scavenger as an oxygen carrier, supports this hypothesis. Haemoglobins from microorganisms, which may probably be the closest vestiges to the ancestral molecules, are divided into three families. The first one is made from flavohaemoglobins, a group of chimerical proteins carrying a globin domain and an oxido-reduction FAD-dependant domain. The second corresponds to truncated haemoglobins, which are hexacoordinated with very high oxygen-affinity molecules, 20-40 residues shorter than classical haemoglobins. The third group is made from bacterial haemoglobins such as that of Vitreoscilla. Some specific structural arrangements in the region surrounding the heme are cause of their high oxygen affinity. In plants, two types of haemoglobins are present (non-symbiotic and symbiotic), that arose from duplication of an ancestral vegetal gene. Non-symbiotic haemoglobins, which are probably the oldest, are scarcely distributed within tissues having high energetic consumption. Conversely, symbiotic haemoglobins (also named leghaemoglobins) are present at a high concentration (mM) mostly in the rhizomes of legumes, where they are involved in nitrogen metabolism. In some species, haemoglobin was proposed to be an oxygen sensor bringing to the organism information to adjust metabolism or biosynthesis to the oxygen requirement. Elsewhere haemoglobin may act as final electron acceptors in oxido-reduction pathways. Evolution of haemoglobin in invertebrates followed a large variety of scenarios. Some surprising functions as sulphide acquisition in invertebrates living near hydrothermal vents, or a role in the phototrophism of worm need to be mentioned. In invertebrates, the size of haemoglobin varies from monomers to giant molecules associating up to 144 subunits, while in vertebrates it is always a tetramer. In some species, several haemoglobins, with completely different structure and function, may coexist. This demonstrates how hazardous may be to extrapolate the function of a protein from only structural data.

The 109 residue nerve tissue minihemoglobin from Cerebratulus lacteus highlights striking structural plasticity of the alpha-helical globin fold

A very short hemoglobin (CerHb; 109 amino acids) binds O(2) cooperatively in the nerve tissue of the nemertean worm Cerebratulus lacteus to sustain neural activity during anoxia. Sequence analysis suggests that CerHb tertiary structure may be unique among the known globin fold evolutionary variants. The X-ray structure of oxygenated CerHb (R factor 15.3%, at 1.5 A resolution) displays deletion of the globin N-terminal A helix, an extended GH region, a very short H helix, and heme solvent shielding based on specific aromatic residues. The heme-bound O(2) is stabilized by hydrogen bonds to the distal TyrB10-GlnE7 pair. Ligand access to heme may take place through a wide protein matrix tunnel connecting the distal site to a surface cleft located between the E and H helices.