Plant hemoglobins: a molecular fossil record for the evolution of oxygen transport

Uncovering the roles of hemoglobins in soybean facing water stress

Water stress drastically hinders crop yield, including soybean - one of the world's most relevant feeding crops - threatening the food security of an ever-growing global population. Hemoglobins (GLBs) are involved in water stress tolerance; however, the role they effectively play in soybean remains underexplored. In this study, in silico and in vivo analyses were performed to identify soybean GLBs, capture their transcriptional profile under water stress, and overexpress promising members to assess how soybean cope with waterlogging. Seven GLBs were found, two GLB1 (non-symbiotic) and five GLB2 (symbiotic or leghemoglobins). Three out of the seven GLBs were differentially expressed in soybean RNA-seq libraries of water stress and were evaluated by real-time PCR. Consistently, GmGLB1-1 and GmGLB1-2 were moderately and highly expressed under waterlogging, respectively. Composite plants with roots overexpressing GmGLB1-1 or GmGLB1-2 (mostly) showed higher transcript abundance of stress-defensive genes involved in anaerobic, nitrogen, carbon, and antioxidant metabolism when subjected to waterlogging. In addition, soybean bearing p35S:GmGLB1-2 had lower H2O2 root content, a reactive oxygen species (ROS), under water excess compared with the control condition. Altogether these results suggest that GmGLB1-2 is a strong candidate for soybean genetic engineering to generate waterlogging-tolerant soybean cultivars.

Three classes of hemoglobins are required for optimal vegetative and reproductive growth of Lotus japonicus: genetic and biochemical characterization of LjGlb2-1

Legumes express two major types of hemoglobins, namely symbiotic (leghemoglobins) and non-symbiotic (phytoglobins), with the latter being categorized into three classes according to phylogeny and biochemistry. Using knockout mutants, we show that all three phytoglobin classes are required for optimal vegetative and reproductive development of Lotus japonicus. The mutants of two class 1 phytoglobins showed different phenotypes: Ljglb1-1 plants were smaller and had relatively more pods, whereas Ljglb1-2 plants had no distinctive vegetative phenotype and produced relatively fewer pods. Non-nodulated plants lacking LjGlb2-1 showed delayed growth and alterations in the leaf metabolome linked to amino acid processing, fermentative and respiratory pathways, and hormonal balance. The leaves of mutant plants accumulated salicylic acid and contained relatively less methyl jasmonic acid, suggesting crosstalk between LjGlb2-1 and the signaling pathways of both hormones. Based on the expression of LjGlb2-1 in leaves, the alterations of flowering and fruiting of nodulated Ljglb2-1 plants, the developmental and biochemical phenotypes of the mutant fed on ammonium nitrate, and the heme coordination and reactivity of the protein toward nitric oxide, we conclude that LjGlb2-1 is not a leghemoglobin but an unusual class 2 phytoglobin. For comparison, we have also characterized a close relative of LjGlb2-1 in Medicago truncatula, MtLb3, and conclude that this is an atypical leghemoglobin.

Kinetic mechanisms for O2 binding to myoglobins and hemoglobins

Antonini and Brunori's 1971 book "Hemoglobin and Myoglobin in Their Reactions with Ligands" was a truly remarkable publication that summarized almost 100 years of research on O₂ binding to these globins. Over the ensuing 50 years, ultra-fast laser photolysis techniques, high-resolution and time resolved X-ray crystallography, molecular dynamics simulations, and libraries of recombinant myoglobin (Mb) and hemoglobin (Hb) variants have provided structural interpretations of O₂ binding to these proteins. The resultant mechanisms provide quantitative descriptions of the stereochemical factors that govern overall affinity, including proximal and distal steric restrictions that affect iron reactivity and favorable positive electrostatic interactions that preferentially stabilize bound O₂. The pathway for O₂ uptake and release by Mb and subunits of Hb has been mapped by screening libraries of site-directed mutants in laser photolysis experiments. O₂ enters mammalian Mb and the α and β subunits of human HbA through a channel created by upward and outward rotation of the distal His at the E7 helical position, is non-covalently captured in the interior of the distal cavity, and then internally forms a bond with the heme Fe(II) atom. O₂ dissociation is governed by disruption of hydrogen bonding interactions with His (E7), breakage of the Fe(II)-O₂ bond, and then competition between rebinding and escape through the E7-gate. The structural features that govern the rates of both the individual steps and overall reactions have been determined and provide the framework for: (1) defining the physiological functions of specific globins and their evolution; (2) understanding the clinical features of hemoglobinopathies; and (3) designing safer and more efficient acellular hemoglobin-based oxygen carriers (HBOCs) for transfusion therapy, organ preservation, and other commercially relevant O₂ transport and storage processes.

Core-modified porphyrins: novel building blocks in chemistry

Various (metallo)porphyrins and related compounds have been intensively investigated by different research groups due to their extremely important role in living organisms along with their versatile applications in technology. The design of novel porphyrinoids by core-modification, or substitution of pyrrole nitrogens, with the elements of other groups of the Periodic Table has been considered as a highly promising methodology for tuning structures and properties of porphyrinoids and thus opening new possible applications for them. Much effort has been given to the modifications of the porphyrin core with elements of the main groups, namely O, S, Se (chalcogens), and the heavier congener of nitrogen, phosphorus. In general, the porphyrin core modification by replacing nitrogens with heteroatoms is a promising and effective strategy for obtaining new compounds with unusual structures and properties (optical, electrochemical, coordinating, etc.) as well as reactivity. These novel molecules can also be employed as promising building or construction blocks in various applications in the nanotechnology area.

Comparison of P- and As-core-modified porphyrins with the parental porphyrin: a computational study

The first comparative DFT (B3LYP/6-31G*) study of the Zn-porphyrin and its two derivatives, ZnP(P)4 and ZnP(As)4, is reported. For all three species studied, ZnP, ZnP(P)4 and ZnP(As)4, the singlet was calculated to be the lowest-energy structure and singlet-triplet gap was found to decrease from ca. 41—42 kcal/mol for N to ca. 17—18 kcal/mol for P and to ca. 10 kcal/mol for As. Both ZnP(P)4 and ZnP(As)4 were calculated to attain very pronounced bowl-like shapes. The frontier molecular orbitals (MOs) of the core-modified porphyrins are quite similar to the ZnP frontier MOs. For the HOMO-2 of the core-modified porphyrins due to the ZnP(P)4/ZnP(As)4 bowl-like shapes we might suppose the existence of “internal” electron delocalization inside the ZnP(P)4/ZnP(As)4 “bowls”. Noticeable reduction of the HOMO/LUMO gaps was calculated for ZnP(P)4 and ZnP(As)4, by ca. 1.10 and 1.47 eV, respectively, compared to ZnP. The core-modification of porphyrins by P and especially by As was found to result in significant decrease of the charge on Zn-centers, by ca. 0.61—0.67e for P and by ca. 0.69—0.76e for As. Charges on P- and As-centers were computed to have large positive values, ca. 0.41—0.45e and ca. 0.43—0.47e, for P and As, respectively, compared to significant negative values, ca. −0.65 to −0.66e for N. The porphyrin core-modification by heavier N congeners, P and As, can noticeably modify the structures, electronic, and optical properties of porphyrins, thus affecting their reactivity and potential applications.

Hemoglobins in the legume-Rhizobium symbiosis

Legume nodules have two types of hemoglobins: symbiotic or leghemoglobins (Lbs) and nonsymbiotic or phytoglobins (Glbs). The latter are categorized into three phylogenetic classes differing in heme coordination and O₂ affinity. This review is focused on the roles of Lbs and Glbs in the symbiosis of rhizobia with crop legumes and the model legumes for indeterminate (Medicago truncatula) and determinate (Lotus japonicus) nodulation. Only two hemoglobin functions are well established in nodules: Lbs deliver O₂ to the bacteroids and act as O₂ buffers, preventing nitrogenase inactivation; and Glb1-1 modulates nitric oxide concentration during symbiosis, from the early stage, avoiding the plant's defense response, to nodule senescence. Here, we critically examine early and recent results, update and correct the information on Lbs and Glbs with the latest genome versions, provide novel expression data and identify targets for future research. Crucial unresolved questions include the expression of multiple Lbs in nodules, their presence in the nuclei and in uninfected nodule cells, and, intriguingly, their expression in nonsymbiotic tissues. RNA-sequencing data analysis shows that Lbs are expressed as early as a few hours after inoculation and that their mRNAs are also detectable in roots and pods, which clearly suggests that these heme proteins play additional roles unrelated to nitrogen fixation. Likewise, issues awaiting investigation are the functions of other Glbs in nodules, the spatiotemporal expression profiles of Lbs and Glbs at the mRNA and protein levels, and the molecular mechanisms underlying their regulation during nodule development and in response to stress and hormones.

Plant hemoglobins: a journey from unicellular green algae to vascular plants

Globins (Glbs) are widely distributed in archaea, bacteria and eukaryotes. They can be classified into proteins with 2/2 or 3/3 α-helical folding around the heme cavity. Both types of Glbs occur in green algae, bryophytes and vascular plants. The Glbs of angiosperms have been more intensively studied, and several protein structures have been solved. They can be hexacoordinate or pentacoordinate, depending on whether a histidine is coordinating or not at the sixth position of the iron atom. The 3/3 Glbs of class 1 and the 2/2 Glbs (also called class 3 in plants) are present in all angiosperms, whereas the 3/3 Glbs of class 2 have been only found in early angiosperms and eudicots. The three Glb classes are expected to play different roles. Class 1 Glbs are involved in hypoxia responses and modulate NO concentration, which may explain their roles in plant morphogenesis, hormone signaling, cell fate determination, nutrient deficiency, nitrogen metabolism and plant-microorganism symbioses. Symbiotic Glbs derive from class 1 or class 2 Glbs and transport O₂ in nodules. The physiological roles of class 2 and class 3 Glbs are poorly defined but could involve O₂ and NO transport and/or metabolism, respectively. More research is warranted on these intriguing proteins to determine their non-redundant functions.

Sugar beet hemoglobins: reactions with nitric oxide and nitrite reveal differential roles for nitrogen metabolism

In contrast with human hemoglobin (Hb) in red blood cells, plant Hbs do not transport oxygen, instead research points towards nitrogen metabolism. Using comprehensive and integrated biophysical methods we characterized three sugar beet Hbs: BvHb1.1, BvHb1.2 and BvHb2. Their affinities for oxygen, CO, and hexacoordination were determined. Their role in nitrogen metabolism was studied by assessing their ability to bind NO, to reduce nitrite (NiR, nitrite reductase), and to form nitrate (NOD, NO dioxygenase). Results show that BvHb1.2 has high NOD-like activity, in agreement with the high nitrate levels found in seeds where this protein is expressed. BvHb1.1, on the other side, is equally capable to bind NO as to form nitrate, its main role would be to protect chloroplasts from the deleterious effects of NO. Finally, the ubiquitous, reactive, and versatile BvHb2, able to adopt 'open and closed forms', would be part of metabolic pathways where the balance between oxygen and NO is essential. For all proteins, the NiR activity is relevant only when nitrite is present at high concentrations and both NO and oxygen are absent. The three proteins have distinct intrinsic capabilities to react with NO, oxygen and nitrite; however, it is their concentration which will determine the BvHbs' activity.

Enhancing Statistical Multiple Sequence Alignment and Tree Inference Using Structural Information

For highly divergent sequences, there is often insufficient information to reliably construct alignments and phylogenetic trees. Since protein structure may be strongly conserved despite large divergences in sequence, structural information can be used to help identify homology in such cases.While there exist well-studied models of sequence evolution, structurally informed alignment methods have typically made use of geometric measures of deviation that do not take into account the underlying mutational processes. In order to integrate structural information into sequence-based evolutionary models, we recently developed a stochastic model of structural evolution on a phylogenetic tree and implemented this as the StructAlign plugin for the StatAlign statistical alignment package.In this chapter, we will outline the types of analyses that can be carried out using StructAlign, illustrating how the inclusion of structural information can be used to inform joint estimation of alignments and trees. StructAlign can also be used to infer branch-specific rates of structural evolution, and analysis of an example globin dataset highlights strong variation in the inferred rate across the tree. While structure is more highly conserved within clades, the rate of structural divergence as a function of sequence variation is larger between functionally divergent proteins. Allowing for the rate of structural divergence to vary over the tree results in an improved fit to the empirically observed pairwise RMSD values.

Characterization of the Heme Pocket Structure and Ligand Binding Kinetics of Non-symbiotic Hemoglobins from the Model Legume Lotus japonicus

Plant hemoglobins (Hbs) are found in nodules of legumes and actinorhizal plants but also in non-symbiotic organs of monocots and dicots. Non-symbiotic Hbs (nsHbs) have been classified into two phylogenetic groups. Class 1 nsHbs show an extremely high O2 affinity and are induced by hypoxia and nitric oxide (NO), whereas class 2 nsHbs have moderate O2 affinity and are induced by cold and cytokinins. The functions of nsHbs are still unclear, but some of them rely on the capacity of hemes to bind diatomic ligands and catalyze the NO dioxygenase (NOD) reaction (oxyferrous Hb + NO → ferric Hb + nitrate). Moreover, NO may nitrosylate Cys residues of proteins. It is therefore important to determine the ligand binding properties of the hemes and the role of Cys residues. Here, we have addressed these issues with the two class 1 nsHbs (LjGlb1-1 and LjGlb1-2) and the single class 2 nsHb (LjGlb2) of Lotus japonicus, which is a model legume used to facilitate the transfer of genetic and biochemical information into crops. We have employed carbon monoxide (CO) as a model ligand and resonance Raman, laser flash photolysis, and stopped-flow spectroscopies to unveil major differences in the heme environments and ligand binding kinetics of the three proteins, which suggest non-redundant functions. In the deoxyferrous state, LjGlb1-1 is partially hexacoordinate, whereas LjGlb1-2 shows complete hexacoordination (behaving like class 2 nsHbs) and LjGlb2 is mostly pentacoordinate (unlike other class 2 nsHbs). LjGlb1-1 binds CO very strongly by stabilizing it through hydrogen bonding, but LjGlb1-2 and LjGlb2 show lower CO stabilization. The changes in CO stabilization would explain the different affinities of the three proteins for gaseous ligands. These affinities are determined by the dissociation rates and follow the order LjGlb1-1 > LjGlb1-2 > LjGlb2. Mutations LjGlb1-1 C78S and LjGlb1-2 C79S caused important alterations in protein dynamics and stability, indicating a structural role of those Cys residues, whereas mutation LjGlb1-1 C8S had a smaller effect. The three proteins and their mutant derivatives exhibited similarly high rates of NO consumption, which were due to NOD activity of the hemes and not to nitrosylation of Cys residues.

NO dioxygenase- and peroxidase-like activity of Arabidopsis phytoglobin 3 and its role in Sclerotinia sclerotiorum defense

Phytoglobin 3 appears to be ubiquitous in plants, yet there has been dearth of evidence for their potent physiological functions. Previous crystallographic studies suggest a potential NO dioxygenase like activity of Arabidopsis phytoglobin 3 (AHb3). The present work examined the in vivo function of AHb3 in plant physiology and its role in biotic stress using Arabidopsis- Sclerotinia sclerotorium pathosystem. The gene was found to be ubiquitously expressed in all plant tissues, with moderately increased expression in roots. Its expression was induced upon NO, H₂O₂ and biotic stress. A C-terminal tagged GFP version of the wild type protein revealed its enhanced accumulation in the guard cells. AHb3-GFP was found to be partitioned majorly into the nucleus while residual amounts were present in the cytoplasm. The loss of function AHb3 mutant exhibited reduced root length and fresh weight. AHb3 knockout lines also displayed enhanced susceptibility towards the S. sclerotiorum. Interestingly, these lines displayed enhanced ROS accumulation upon pathogen challenge as suggested by DAB staining. Furthermore, enhanced/decreased NO accumulation in AHb3 knockout/overexpression lines upon treatment with multiple NO donors suggests a potent NO dioxygenase like activity for the protein. Taken together, our data indicate that AHb3 play a crucial role in regulating root length as well as in mediating defense response against S. sclerotiorum, possibly by modulating NO and ROS levels.

Role of Ionic Strength and pH in Modulating Thermodynamic Profiles Associated with CO Escape from Rice Nonsymbiotic Hemoglobin 1

Type 1 nonsymbiotic hemoglobins are found in a wide variety of land plants and exhibit very high affinities for exogenous gaseous ligands. These proteins are presumed to have a role in protecting plant cells from oxidative stress under etiolated/hypoxic conditions through NO dioxygenase activity. In this study we have employed photoacoustic calorimetry, time-resolved absorption spectroscopy, and classical molecular dynamics simulations in order to elucidate thermodynamics, kinetics, and ligand migration pathways upon CO photodissociation from WT and a H73L mutant of type 1 nonsymbiotic hemoglobin from Oryza sativa (rice). We observe a temperature dependence of the resolved thermodynamic parameters for CO photodissociation from CO-rHb1 which we attribute to temperature dependent formation of a network of electrostatic interactions in the vicinity of the heme propionate groups. We also observe slower ligand escape from the protein matrix under mildly acidic conditions in both the WT and H73L mutant (τ = 134 ± 19 and 90 ± 15 ns). Visualization of transient hydrophobic channels within our classical molecular dynamics trajectories allows us to attribute this phenomenon to a change in the ligand migration pathway which occurs upon protonation of the distal His73, His117, and His152. Protonation of these residues may be relevant to the functioning of the protein in vivo given that etiolation/hypoxia can cause a decrease in intracellular pH in plant cells.

Dual targeted poplar ferredoxin NADP(+) oxidoreductase interacts with hemoglobin 1

Previous reports have connected non-symbiotic and truncated hemoglobins (Hbs) to metabolism of nitric oxide (NO), an important signalling molecule involved in wood formation. We have studied the capability of poplar (Populus tremula × tremuloides) Hbs PttHb1 and PttTrHb proteins alone or with a flavin-protein reductase to relieve NO cytotoxicity in living cells. Complementation tests in a Hb-deficient, NO-sensitive yeast (Saccharomyces cerevisiae) Δyhb1 mutant showed that neither PttHb1 nor PttTrHb alone protected cells against NO. To study the ability of Hbs to interact with a reductase, ferredoxin NADP(+) oxidoreductase PtthFNR was characterized by sequencing and proteomics. To date, by far the greatest number of the known dual-targeted plant proteins are directed to chloroplasts and mitochondria. We discovered a novel variant of hFNR that lacks the plastid presequence and resides in cytosol. The coexpression of PttHb1 and PtthFNR partially restored NO resistance of the yeast Δyhb1 mutant, whereas PttTrHb coexpressed with PtthFNR failed to rescue growth. YFP fusion proteins confirmed the interaction between PttHb1 and PtthFNR in plant cells. The structural modelling results indicate that PttHb1 and PtthFNR are able to interact as NO dioxygenase. This is the first report on dual targeting of central plant enzyme FNR to plastids and cytosol.

Structural and Functional Significance of the N- and C-Terminal Appendages in Arabidopsis Truncated Hemoglobin

Plant hemoglobins constitute three distinct groups: symbiotic, nonsymbiotic, and truncated hemoglobins. Structural investigation of symbiotic and nonsymbiotic (class I) hemoglobins revealed the presence of a vertebrate-like 3/3 globin fold in these proteins. In contrast, plant truncated hemoglobins are similar to bacterial truncated hemoglobins with a putative 2/2 α-helical globin fold. While multiple structures have been reported for plant hemoglobins of the first two categories, for plant truncated globins only one structure has been reported of late. Here, we report yet another crystal structure of the truncated hemoglobin from Arabidopsis thaliana (AHb3) with two water molecules in the heme pocket, of which one is distinctly coordinated to the heme iron, unlike the only available crystal structure of AHb3 with a hydroxyl ligand. AHb3 was monomeric in its crystallographic asymmetric unit; however, dimer was evident in the crystallographic symmetry, and the globin indeed existed as a stable dimer in solution. The tertiary structure of the protein exhibited a bacterial-like 2/2 α-helical globin fold with an additional N-terminal α-helical extension and disordered C-termini. To address the role of these extended termini in AHb3, which is yet unknown, N- and C-terminal deletion mutants were created and characterized and molecular dynamics simulations performed. The C-terminal deletion had an insignificant effect on most properties but perturbed the dimeric equilibrium of AHb3 and significantly influenced azide binding kinetics in the ferric state. These results along with the disordered nature of the C-terminus indicated its putative role in intramolecular or intermolecular interactions probably regulating protein-ligand and protein-protein interactions. While the N-terminal deletion did not change the overall globin fold, stability, or ligand binding kinetics, it seemed to have influenced coordination at the heme iron, the hydration status of the active site, and the quaternary structure of AHb3. Evidence indicated that the N-terminus is the predominant factor regulating the quaternary interaction appropriate to physiological requirements, dynamics of the side chains in the heme pocket, and tunnel organization in the protein matrix.

Molecular dynamics simulations indicate that tyrosineB10 limits motions of distal histidine to regulate CO binding in soybean leghemoglobin

Myoglobin (Mb) uses strong electrostatic interaction in its distal heme pocket to regulate ligand binding. The mechanism of regulation of ligand binding in soybean leghemoglobin a (Lba) has been enigmatic and more so due to the absence of gaseous ligand bound atomic resolution three-dimensional structure of the plant globin. While the 20-fold higher oxygen affinity of Lba compared with Mb is required for its dual physiological function, the mechanism by which this high affinity is achieved is only emerging. Extensive mutational analysis combined with kinetic and CO-FT-IR spectroscopic investigation led to the hypothesis that Lba depended on weakened electrostatic interaction between distal HisE7 and bound ligand achieved by invoking B10Tyr, which itself hydrogen bonds with HisE7 thus restricting it in a single conformation detrimental to Mb-like strong electrostatic interaction. Such theory has been re-assessed here using CO-Lba in silico model and molecular dynamics simulation. The investigation supports the presence of at least two major conformations of HisE7 in Lba brought about by imidazole ring flip, one of which makes hydrogen bonds effectively with B10Tyr affecting the former's ability to stabilize bound ligand, while the other does not. However, HisE7 in Lba has limited conformational freedom unlike high frequency of imidazole ring flips observed in Mb and in TyrB10Leu mutant of Lba. Thus, it appears that TyrB10 limits the conformational freedom of distal His in Lba, tuning down ligand dissociation rate constant by reducing the strength of hydrogen bonding to bound ligand, which the freedom of distal His of Mb allows.

Significantly enhanced heme retention ability of myoglobin engineered to mimic the third covalent linkage by nonaxial histidine to heme (vinyl) in synechocystis hemoglobin

Heme proteins, which reversibly bind oxygen and display a particular fold originally identified in myoglobin (Mb), characterize the "hemoglobin (Hb) superfamily." The long known and widely investigated Hb superfamily, however, has been enriched by the discovery and investigation of new classes and members. Truncated Hbs typify such novel classes and exhibit a distinct two-on-two α-helical fold. The truncated Hb from the freshwater cyanobacterium Synechocystis exhibits hexacoordinate heme chemistry and bears an unusual covalent bond between the nonaxial His(117) and a heme porphyrin 2-vinyl atom, which remains tightly associated with the globin unlike any other. It seems to be the most stable Hb known to date, and His(117) is the dominant force holding the heme. Mutations of amino acid residues in the vicinity did not influence this covalent linkage. Introduction of a nonaxial His into sperm whale Mb at the topologically equivalent position and in close proximity to vinyl group significantly increased the heme stability of this prototype globin. Reversed phase chromatography, electrospray ionization-MS, and MALDI-TOF analyses confirmed the presence of covalent linkage in Mb I107H. The Mb mutant with the engineered covalent linkage was stable to denaturants and exhibited ligand binding and auto-oxidation rates similar to the wild type protein. This indeed is a novel finding and provides a new perspective to the evolution of Hbs. The successful attempt at engineering heme stability holds promise for the production of stable Hb-based blood substitute.

Engineering the internal cavity of neuroglobin demonstrates the role of the haem-sliding mechanism

Neuroglobin is a member of the globin family involved in neuroprotection; it is primarily expressed in the brain and retina of vertebrates. Neuroglobin belongs to the heterogeneous group of hexacoordinate globins that have evolved in animals, plants and bacteria, endowed with the capability of reversible intramolecular coordination, allowing the binding of small gaseous ligands (O2, NO and CO). In a unique fashion among haemoproteins, ligand-binding events in neuroglobin are dependent on the sliding of the haem itself within a preformed internal cavity, as revealed by the crystal structure of its CO-bound derivative. Point mutants of the neuroglobin internal cavity have been engineered and their functional and structural characterization shows that hindering the haem displacement leads to a decrease in CO affinity, whereas reducing the cavity volume without interfering with haem sliding has negligible functional effects.

Iron bioavailability of maize hemoglobin in a Caco-2 cell culture model

Maize ( Zea mays ) is an important staple crop in many parts of the world but has low iron bioavailability, in part due to its high phytate content. Hemoglobin is a form of iron that is highly bioavailable, and its bioavailability is not inhibited by phytate. It was hypothesized that maize hemoglobin is a highly bioavailable iron source and that biofortification of maize with iron can be accomplished by overexpression of maize globin in the endosperm. Maize was transformed with a gene construct encoding a translational fusion of maize globin and green fluorescent protein under transcriptional control of the maize 27 kDa γ-zein promoter. Iron bioavailability of maize hemoglobin produced in Escherichia coli and of stably transformed seeds expressing the maize globin-GFP fusion was determined using an in vitro Caco-2 cell culture model. Maize flour fortified with maize hemoglobin was found to have iron bioavailability that is not significantly different from that of flour fortified with ferrous sulfate or bovine hemoglobin but is significantly higher than unfortified flour. Transformed maize grain expressing maize globin was found to have iron bioavailability similar to that of untransformed seeds. These results suggest that maize globin produced in E. coli may be an effective iron fortificant, but overexpressing maize globin in maize endosperm may require a different strategy to increase bioavailable iron content in maize.

X-ray crystallographic structural characteristics of Arabidopsis hemoglobin I and their functional implications

Genome of the model dicot flowering plant, Arabidopsis thaliana, a popular tool for understanding molecular biology of plant physiology, encodes all three classes of plant hemoglobins that differ in their sequence, ligand binding and spectral properties. As such these globins are of considerable attention. Crystal structures of few members of plant class I nonsymbiotic hemoglobin have been described earlier. Here we report the crystal structure of Arabidopsis class I hemoglobin (AHb1) to 2.2Ǻ and compare its key features with the structures of similar nonsymbiotic hemoglobin from other species. Crystal structure of AHb1 is homologous to the related members with similar globin fold and heme pocket architecture. The structure is homodimeric in the asymmetric unit with both distal and proximal histidines coordinating to the heme iron atom. Residues lining the dimeric interface are also conserved in AHb1 with the exception of additional electrostatic interaction between H112 and E113 of each subunit and that involving Y119 through two water molecules. In addition, differences in heme pocket non-covalent interactions, a novel Ser residue at F7 position, Xe binding site variability, internal cavity topology differences, CD loop conformation and stability and other such properties might explain kinetic variability in AHb1. Detailed cavity analysis of AHb1 showed the presence of a novel long tunnel connecting the distal pockets of both the monomers. Presence of such tunnel, along with conformational heterogeneity observed in the two chains, might suggest cooperative ligand binding and support its role in NO scavenging. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

CO rebinding kinetics and molecular dynamics simulations highlight dynamic regulation of internal cavities in human cytoglobin

Cytoglobin (Cygb) was recently discovered in the human genome and localized in different tissues. It was suggested to play tissue-specific protective roles, spanning from scavenging of reactive oxygen species in neurons to supplying oxygen to enzymes in fibroblasts. To shed light on the functioning of such versatile machinery, we have studied the processes supporting transport of gaseous heme ligands in Cygb. Carbon monoxide rebinding shows a complex kinetic pattern with several distinct reaction intermediates, reflecting rebinding from temporary docking sites, second order recombination, and formation (and dissociation) of a bis-histidyl heme hexacoordinated reaction intermediate. Ligand exit to the solvent occurs through distinct pathways, some of which exploit temporary docking sites. The remarkable change in energetic barriers, linked to heme bis-histidyl hexacoordination by HisE7, may be responsible for active regulation of the flux of reactants and products to and from the reaction site on the distal side of the heme. A substantial change in both protein dynamics and inner cavities is observed upon transition from the CO-liganded to the pentacoordinated and bis-histidyl hexacoordinated species, which could be exploited as a signalling state. These findings are consistent with the expected versatility of the molecular activity of this protein.

The globins of cyanobacteria and algae

Approximately, 20 years ago, a haemoglobin gene was identified within the genome of the cyanobacterium Nostoc commune. Haemoglobins have now been confirmed in multiple species of photosynthetic microbes beyond N. commune, and the diversity of these proteins has recently come under increased scrutiny. This chapter summarizes the state of knowledge concerning the phylogeny, physiology and chemistry of globins in cyanobacteria and green algae. Sequence information is by far the best developed and the most rapidly expanding aspect of the field. Structural and ligand-binding properties have been described for just a few proteins. Physiological data are available for even fewer. Although activities such as nitric oxide dioxygenation and oxygen scavenging are strong candidates for cellular function, dedicated studies will be required to complete the story on this intriguing and ancient group of proteins.

Nitrite reductase activity of nonsymbiotic hemoglobins from Arabidopsis thaliana

Plant nonsymbiotic hemoglobins possess hexacoordinate heme geometry similar to that of the heme protein neuroglobin. We recently discovered that deoxygenated neuroglobin converts nitrite to nitric oxide (NO), an important signaling molecule involved in many processes in plants. We sought to determine whether Arabidopsis thaliana nonsymbiotic hemoglobins classes 1 and 2 (AHb1 and AHb2, respectively) might function as nitrite reductases. We found that the reaction of nitrite with deoxygenated AHb1 and AHb2 generates NO gas and iron-nitrosyl-hemoglobin species. The bimolecular rate constants for reduction of nitrite to NO are 19.8 ± 3.2 and 4.9 ± 0.2 M(-1) s(-1), respectively, at pH 7.4 and 25 °C. We determined the pH dependence of these bimolecular rate constants and found a linear correlation with the concentration of protons, indicating the requirement for one proton in the reaction. The release of free NO gas during the reaction under anoxic and hypoxic (2% oxygen) conditions was confirmed by chemiluminescence detection. These results demonstrate that deoxygenated AHb1 and AHb2 reduce nitrite to form NO via a mechanism analogous to that observed for hemoglobin, myoglobin, and neuroglobin. Our findings suggest that during severe hypoxia and in the anaerobic plant roots, especially in species submerged in water, nonsymbiotic hemoglobins provide a viable pathway for NO generation via nitrite reduction.

Plant hemoglobins: important players at the crossroads between oxygen and nitric oxide

Plant hemoglobins constitute a diverse group of hemeproteins and evolutionarily belong to three different classes. Class 1 hemoglobins possess an extremely high affinity to oxygen and their main function consists in scavenging of nitric oxide (NO) at very low oxygen levels. Class 2 hemoglobins have a lower oxygen affinity and they facilitate oxygen supply to developing tissues. Symbiotic hemoglobins in nodules have mostly evolved from class 2 hemoglobins. Class 3 hemoglobins are truncated and represent a clade with a very low similarity to class 1 and 2 hemoglobins. They may regulate oxygen delivery at high O(2) concentrations. Depending on their physical properties, hemoglobins belong either to hexacoordinate non-symbiotic or pentacoordinate symbiotic groups. Plant hemoglobins are plausible targets for improving resistance to multiple stresses.

Oxygen binding to Arabidopsis thaliana AHb2 nonsymbiotic hemoglobin: evidence for a role in oxygen transport

Nonsymbiotic hemoglobins AHb1 and AHb2 discovered in Arabidopsis thaliana are likely to carry out distinct physiological roles, in consideration of their differences in sequence, structure, expression pattern, and tissue localization. Despite a relatively fast autoxidation in the presence of O(2) , we were able to collect O(2) -binding curves for AHb2 in the presence of a reduction enzymatic system. AHb2 binds O(2) noncooperatively with a p50 of 0.021 ± 0.003 Torr, a value consistent with a recently proposed role in O(2) transport. The analysis of the internal cavities derived from the structures sampled in molecular dynamics simulations confirms strong differences with AHb1, proposed to work as a NO deoxygenase in vivo. Overall, our results are consistent with a role for AHb2 as an oxygen carrier, as recently proposed on the basis of experiments on AHb2-overexpressing mutants of A. thaliana.

Structural analysis in nonsymbiotic hemoglobins: what can we learn from inner cavities?

Plants contain three classes of hemoglobins which are not associated with nitrogen fixing bacteria, and have been accordingly termed nonsymbiotic hemoglobins. The function of nonsymbiotic hemoglobins is as yet mostly unknown. A NO dioxygenase activity has been proposed and demonstrated for some of them in vitro. In this context, a sound molecular mechanism that relates the structure with the biological activity is crucial to suggest a given physiological role. Insight into such a mechanism is now facilitated by recent progress made in both experimental and computational techniques. These studies have highlighted a number of key structural features implicated in the function of nonsymbiotic hemoglobins. The bis-histidyl hexacoordination of the heme in both its ferric and ferrous states provides a powerful and general tool to modulate reactivity, protein dynamics, and shape of the cavities. In addition, the specific arrangement of distal cavity residues provides effective protection against autoxidation. Inspection of the static crystal structures available for both liganded and unliganded states seems unsufficient to explain the function of these proteins. Function appears to be intimately linked with protein flexibility, which influences the dynamical behavior of inner cavities, capable of delivering apolar reactants to the reaction site, and removing charged reaction products. In this mini review, we demonstrate how the integration of information derived from experimental assays and computational studies is valuable and can shed light into the linkage between structural plasticity of nonsymbiotic hemoglobins and their biological role.

Contrasting growth and adaptive responses of two oak species to flooding stress: role of non-symbiotic haemoglobin

Soil flooding is an environmental constraint that is increasingly important for forest ecosystems, affecting tree growth and regeneration. As a result, selection pressure will alter forest diversity and distribution by favouring tree species tolerant of soil oxygen deprivation. Sessile and pedunculate oaks are the most abundant oak species and they exhibit a strong differential tolerance to waterlogging. In order to gain some understanding of the mechanisms of tolerance of both species to hypoxia, we undertook the characterization of the physiological, morphological, cellular and molecular responses of both species to flooding stress. Our results indicate that pedunculate oak, the more tolerant species, succeeded in maintaining its growth, water status and photosynthetic activity at a higher level than sessile oak. Furthermore, pedunculate oak developed aerenchyma in its root cortex as well as adventitious roots. The later exhibited a strong accumulation of class1 non-symbiotic haemoglobin localized by in situ hybridization in the protoderm and in some cortical cells. In conclusion, the higher tolerance of pedunculate oak to flooding was associated with an enhanced capacity to maintain photosynthesis and water homeostasis, coupled with the development of adaptive features (aerenchyma, adventitious roots) and with a higher expression of non-symbiotic haemoglobin in the roots.

Phylogenetic relationships of 3/3 and 2/2 hemoglobins in Archaeplastida genomes to bacterial and other eukaryote hemoglobins

Land plants and algae form a supergroup, the Archaeplastida, believed to be monophyletic. We report the results of an analysis of the phylogeny of putative globins in the currently available genomes to bacterial and other eukaryote hemoglobins (Hbs). Archaeplastida genomes have 3/3 and 2/2 Hbs, with the land plant genomes having group 2 2/2 Hbs, except for the unexpected occurrence of two group 1 2/2 Hbs in Ricinus communis. Bayesian analysis shows that plant 3/3 Hbs are related to vertebrate neuroglobins and bacterial flavohemoglobins (FHbs). We sought to define the bacterial groups, whose ancestors shared the precursors of Archaeplastida Hbs, via Bayesian and neighbor-joining analyses based on COBALT alignment of representative sets of bacterial 3/3 FHb-like globins and group 1 and 2 2/2 Hbs with the corresponding Archaeplastida Hbs. The results suggest that the Archaeplastida 3/3 and group 1 2/2 Hbs could have originated from the horizontal gene transfers (HGTs) that accompanied the two generally accepted endosymbioses of a proteobacterium and a cyanobacterium with a eukaryote ancestor. In contrast, the origin of the group 2 2/2 Hbs unexpectedly appears to involve HGT from a bacterium ancestral to Chloroflexi, Deinococcales, Bacilli, and Actinomycetes. Furthermore, although intron positions and phases are mostly conserved among the land plant 3/3 and 2/2 globin genes, introns are absent in the algal 3/3 genes and intron positions and phases are highly variable in their 2/2 genes. Thus, introns are irrelevant to globin evolution in Archaeplastida.

The globin gene family of the cephalochordate amphioxus: implications for chordate globin evolution

BACKGROUND

The lancelet amphioxus (Cephalochordata) is a close relative of vertebrates and thus may enhance our understanding of vertebrate gene and genome evolution. In this context, the globins are one of the best studied models for gene family evolution. Previous biochemical studies have demonstrated the presence of an intracellular globin in notochord tissue and myotome of amphioxus, but the corresponding gene has not yet been identified. Genomic resources of Branchiostoma floridae now facilitate the identification, experimental confirmation and molecular evolutionary analysis of its globin gene repertoire.

RESULTS

We show that B. floridae harbors at least fifteen paralogous globin genes, all of which reveal evidence of gene expression. The protein sequences of twelve globins display the conserved characteristics of a functional globin fold. In phylogenetic analyses, the amphioxus globin BflGb4 forms a common clade with vertebrate neuroglobins, indicating the presence of this nerve globin in cephalochordates. Orthology is corroborated by conserved syntenic linkage of BflGb4 and flanking genes. The kinetics of ligand binding of recombinantly expressed BflGb4 reveals that this globin is hexacoordinated with a high oxygen association rate, thus strongly resembling vertebrate neuroglobin. In addition, possible amphioxus orthologs of the vertebrate globin X lineage and of the myoglobin/cytoglobin/hemoglobin lineage can be identified, including one gene as a candidate for being expressed in notochord tissue. Genomic analyses identify conserved synteny between amphioxus globin-containing regions and the vertebrate β-globin locus, possibly arguing against a late transpositional origin of the β-globin cluster in vertebrates. Some amphioxus globin gene structures exhibit minisatellite-like tandem duplications of intron-exon boundaries ("mirages"), which may serve to explain the creation of novel intron positions within the globin genes.

CONCLUSIONS

The identification of putative orthologs of vertebrate globin variants in the B. floridae genome underlines the importance of cephalochordates for elucidating vertebrate genome evolution. The present study facilitates detailed functional studies of the amphioxus globins in order to trace conserved properties and specific adaptations of respiratory proteins at the base of chordate evolution.

Heme-based sensing by the mammalian circadian protein CLOCK

Heme is emerging as a key player in the synchrony of circadian-coupled transcriptional regulation. Current evidence suggests that levels of circadian-linked transcription are regulated in response to both the availability of intracellular heme and heme-based sensing of carbon monoxide (CO) and possibly nitric oxide (NO). The protein CLOCK is central to the regulation and maintenance of circadian rhythms in mammals. CLOCK comprises two PAS domains, each with a heme binding site. Our studies focus on the functionality of the murine CLOCK PAS-A domain (residues 103-265). We show that CLOCK PAS-A binds iron(III) protoporhyrin IX to form a complex with 1:1 stoichiometry. Optical absorbance and resonance Raman studies reveal that the heme of ferric CLOCK PAS-A is a six-coordinate, low-spin complex whose resonance Raman signature is insensitive to pH over the range of protein stability. Ferrous CLOCK PAS-A is a mixture of five-coordinate, high-spin and six-coordinate, low-spin complexes. Ferrous CLOCK PAS-A forms complexes with CO and NO. Ferric CLOCK PAS-A undergoes reductive nitrosylation in the presence of NO to generate a CLOCK PAS-A-NO, which is a five-coordinate {FeNO}(7) complex. Formation of the highly stable {FeNO}(7) heme complex from either ferrous or ferric heme makes possible the binding of NO at very low concentration, a characteristic of NO sensors. Comparison of the spectroscopic properties and CO-binding kinetics of CLOCK PAS-A with other CO sensor proteins reveals that CLOCK PAS-A exhibits chemical properties consistent with a heme-based gas sensor protein.

Ligand migration and hexacoordination in type 1 non-symbiotic rice hemoglobin

Type 1 non-symbiotic rice hemoglobin (rHb1) shows bis-histidyl heme hexacoordination and is capable of binding diatomic ligands reversibly. The biological function is as yet unclear, but the high oxygen affinity makes it unlikely to be involved in oxygen transport. In order to gain insight into possible physiological roles, we have studied CO rebinding kinetics after laser flash photolysis of rHb1 in solution and encapsulated in silica gel. CO rebinding to wt rHb1 in solution occurs through a fast geminate phase with no sign of rebinding from internal docking sites. Encapsulation in silica gel enhances migration to internal cavities. Site-directed mutagenesis of FB10, a residue known to have a key role in the regulation of hexacoordination and ligand affinity, resulted in substantial effects on the rebinding kinetics, partly inhibiting ligand exit to the solvent, enhancing geminate rebinding and enabling ligand migration within the internal cavities. The mutation of HE7, one of the histidyl residues involved in the hexacoordination, prevents hexacoordination, as expected, but also exposes ligand migration through a complex system of cavities. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Ligand migration and binding in nonsymbiotic hemoglobins of Arabidopsis thaliana

We have studied carbon monoxide (CO) migration and binding in the nonsymbiotic hemoglobins AHb1 and AHb2 of Arabidopsis thaliana using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures. Both proteins have similar amino acid sequences but display pronounced differences in ligand binding properties, at both physiological and cryogenic temperatures. Near neutral pH, the distal HisE7 side chain is close to the heme-bound ligand in the majority of AHb1-CO molecules, as indicated by a low CO stretching frequency at 1921 cm(-1). In this fraction, two CO docking sites can be populated, the primary site B and the secondary site C. When the pH is lowered, a high-frequency stretching band at approximately 1964 cm(-1) grows at the expense of the low-frequency band, indicating that HisE7 protonates and, concomitantly, moves away from the bound ligand. Geminate rebinding barriers are markedly different for the two conformations, and docking site C is not accessible in the low-pH conformation. Rebinding of NO ligands was observed only from site B of AHb1, regardless of conformation. In AHb2, the HisE7 side chain is removed from the bound ligand; rebinding barriers are low, and CO molecules can populate only primary docking site B. These results are interpreted in terms of differences in the active site structures and physiological functions.

Structure and reactivity of hexacoordinate hemoglobins

The heme prosthetic group in hemoglobins is most often attached to the globin through coordination of either one or two histidine side chains. Those proteins with one histidine coordinating the heme iron are called "pentacoordinate" hemoglobins, a group represented by red blood cell hemoglobin and most other oxygen transporters. Those with two histidines are called "hexacoordinate hemoglobins", which have broad representation among eukaryotes. Coordination of the second histidine in hexacoordinate Hbs is reversible, allowing for binding of exogenous ligands like oxygen, carbon monoxide, and nitric oxide. Research over the past several years has produced a fairly detailed picture of the structure and biochemistry of hexacoordinate hemoglobins from several species including neuroglobin and cytoglobin in animals, and the nonsymbiotic hemoglobins in plants. However, a clear understanding of the physiological functions of these proteins remains an elusive goal.

Straight-chain alkyl isocyanides open the distal histidine gate in crystal structures of myoglobin

Crystal structures of methyl, ethyl, propyl, and butyl isocyanide bound to sperm whale myoglobin (Mb) reveal two major conformations. In the in conformer, His(E7) is in a "closed" position, forcing the ligand alkyl chain to point inward. In the out conformer, His(E7) is in an "open" position, allowing the ligand side chain to point outward. A progressive increase in the population of the out conformer is observed with increasing ligand length in P2(1) crystals of native Mb at pH 7.0. This switch from in to out with increasing ligand size also occurs in solution as measured by the decrease in the relative intensity of the low-frequency ( approximately 2075 cm(-1)) versus high-frequency ( approximately 2125 cm(-1)) isocyano bands. In contrast, all four isocyanides in P6 crystals of wild-type recombinant Mb occupy the in conformation. However, mutating either His64 to Ala, creating a "hole" to solvent, or Phe46 to Val, freeing rotation of His64, causes bound butyl isocyanide to point completely outward in P6 crystals. Thus, the unfavorable hindrance caused with crowding a large alkyl side chain into the distal pocket appears to be roughly equal to that for pushing open the His(E7) gate and is easily affected by crystal packing. This structural conclusion supports the "side path" kinetic mechanism for O(2) release, in which the dissociated ligand first moves toward the protein interior and then encounters steric resistance, which is roughly equal to that for escaping to solvent through the His(E7) channel.

Review: correlations between oxygen affinity and sequence classifications of plant hemoglobins

Plants express three phylogenetic classes of hemoglobins (Hb) based on sequence analyses. Class 1 and 2 Hbs are full-length globins with the classical eight helix Mb-like fold, whereas Class 3 plant Hbs resemble the truncated globins found in bacteria. With the exception of the specialized leghemoglobins, the physiological functions of these plant hemoglobins remain unknown. We have reviewed and, in some cases, measured new oxygen binding properties of a large number of Class 1 and 2 plant nonsymbiotic Hbs and leghemoglobins. We found that sequence classification correlates with distinct extents of hexacoordination with the distal histidine and markedly different overall oxygen affinities and association and dissociation rate constants. These results suggest strong selective pressure for the evolution of distinct physiological functions. The leghemoglobins evolved from the Class 2 globins and show no hexacoordination, very high rates of O(2) binding ( approximately 250 muM(-1) s(-1)), moderately high rates of O(2) dissociation ( approximately 5-15 s(-1)), and high oxygen affinity (K(d) or P(50) approximately 50 nM). These properties both facilitate O(2) diffusion to respiring N(2) fixing bacteria and reduce O(2) tension in the root nodules of legumes. The Class 1 plant Hbs show weak hexacoordination (K(HisE7) approximately 2), moderate rates of O(2) binding ( approximately 25 muM(-1) s(-1)), very small rates of O(2) dissociation ( approximately 0.16 s(-1)), and remarkably high O(2) affinities (P(50) approximately 2 nM), suggesting a function involving O(2) and nitric oxide (NO) scavenging. The Class 2 Hbs exhibit strong hexacoordination (K(HisE7) approximately 100), low rates of O(2) binding ( approximately 1 muM(-1) s(-1)), moderately low O(2) dissociation rate constants ( approximately 1 s(-1)), and moderate, Mb-like O(2) affinities (P(50) approximately 340 nM), perhaps suggesting a sensing role for sustained low, micromolar levels of oxygen.

Nonsymbiotic hemoglobins and stress tolerance in plants

Hemoglobins (Hbs) are heme containing proteins found in most organisms including animals, bacteria, and plants. Their structure, size, and function are quite diverse among the different organisms. There are three different types of hemoglobins in plants: symbiotic (sHb), nonsymbiotic (nsHb), and truncated hemoglobins (trHb). The nonsymbiotic hemoglobins are divided into: class 1 hemoglobins (nsHb-1s), which have a very high affinity for oxygen: and class 2 hemoglobins (nsHb-2s), which have lower affinity for oxygen, are similar to the sHbs. nsHb-1s are expressed under hypoxia, osmotic stress, nutrient deprivation, cold stress, rhizobial infection, nitric oxide exposure, and fungal infection. Tolerance to stress is very important for the survival of the plant. Hemoglobins are one of many different strategies that plants have evolved to overcome stress conditions and survive. Hbs also react with NO produced under different stress conditions. Class 1 nsHbs are involved in a metabolic pathway involving NO. Those hemoglobins provide an alternative type of respiration to mitochondrial electron transport under limiting oxygen concentrations. Class 1 nsHbs in hypoxic plants act as part of a soluble, terminal, NO dioxygenase system, yielding nitrate from the reaction of oxyHb with NO. The overall reaction sequence, referred to as the nsHb/NO cycle, consumes NADH and maintains ATP levels via an as yet unknown mechanism. Class 2 nsHbs seem to scavenge NO in a similar fashion as class 1 Hbs and are involved in reducing flowering time in Arabidopsis. nsHbs also show peroxidase-like activity and NO metabolism and possibly protect against nitrosative stress in plant-pathogen interaction and in symbiotic interactions. nsHbs may be involved in other stress conditions such as osmotic, nutrient and cold stress together with NO and the function of nsHbs can be in NO metabolism and signal transduction. However, other possible functions cannot be precluded as Hbs have many different functions in other organisms.