Phosphate Binding Site Research Articles

Organization and Structural Features of Phosphofructokinase and other Glycolytic Enzymes to Meet their Role in Energy Metabolism

Although glycolytic enzymes can be isolated as cytosolic components, some of them are known to bind to proteins from the cytoplasm and other cellular compartments, such as the mitochondria, the nucleus, or the plasma membrane. In some cases these protein-protein interactions are considered significant for the regulation of energy metabolism, as well as for the modulation of other cell functions. Among glycolytic enzymes, phosphofructokinase (PFK) is thought to play a fundamental role in the control of this pathway, because of the number of metabolic signals that can regulate its complex allosteric behavior. No X-ray structure of a characteristic eukaryotic PFK is available yet. However, sequence data indicated that the eukaryotic enzyme originated by duplication, fusion and divergence of an ancestral prokaryotic gene, such that the duplicated fructose 6-phosphate catalytic site in the C-terminal half became an allosteric site for the activator fructose 2,6-bisphosphate. It has been suggested that both sites are shared in the interface between subunits aligned in an antiparallel orientation. To test the contribution of each terminal domain to these two binding sites, chimeric mammalian PFKs involving exchange of their terminal domains have been found to exhibit affinity properties for fructose 6-phosphate and fructose 2,6-bisphosphate that resembled those of the native isozyme that donated the N-terminal half and the C-terminal half, respectively. Further mutational analyses of muscle PFK led us to investigate the composition of the fructose phosphate binding sites and to gain insight into the structural organization of eukaryotic PFK. Data will also be presented evidencing the relationship between mitochondrial function and the operation of glycolysis, as observed after impairment of mitochondrial DNA transcription.

Serine Palmitoyltransferase, the First Step Enzyme in Sphingolipid Biosynthesis, Is Involved in Nonhost Resistance

An overexpression screen of Nicotiana benthamiana cDNAs identified a gene for the LCB2 subunit of serine palmitoyltransferase (SPT) as a potent inducer of hypersensitive response-like cell death. The pyridoxal 5'-phosphate binding site of NbLCB2 is required for its function as a cell death inducer. NbLCB2 mRNA is accumulated after infection by nonhost pathogen Pseudomonas cichorii. Resistance of N. benthamiana against P. cichorii was compromised by treatment with an SPT inhibitor and in NbLCB2- and NbLCB1-silenced plants. These results suggest that biosynthesis of sphingolipids is necessary for the nonhost resistance of N. benthamiana against P. cichorii.

Characterization of coenzyme binding and selectivity determinants in Mycobacterium tuberculosis flavoprotein reductase A: analysis of Arg199 and Arg200 mutants at the NADP(H) 2′-phosphate binding site

Mycobacterium tuberculosis FprA (flavoprotein reductase A) is an NAD(P)H- and FAD-binding reductase that is structurally/evolutionarily related to adrenodoxin reductase. Structural analysis implicates Arg(199) and Arg(200) in interactions with the NADP(H) 2'-phosphate group. R199A, R200A and R199A/R200A mutants were characterized to explore the roles of these basic residues. All mutations abolished neutral FAD semiquinone stabilization in the NADPH-reduced enzyme, owing to weakened NADPH affinity. Instead, FAD hydroquinone was formed in all mutants, and each displayed substantially enhanced autooxidation rates (20-40-fold) compared with NADPH-reduced WT (wild-type) FprA. Steady-state ferricyanide reduction studies revealed diminished NADPH affinity (higher K(m) values), but lower NADH K(m) values. Despite a lowered k(cat), the R199A/R200A mutant exhibited a 200-fold coenzyme specificity switch towards NADH, although substrate inhibition was observed at high NADH concentrations (K(i)=250 microM). Stopped-flow FAD reduction studies confirmed substantially increased NADPH K(d) values, although the limiting flavin reduction rate constant was similar in all mutants. The R199A mutation abolished electron transfer between hydroquinone FprA and NADP+, while this reaction progressed (via an FADH(2)-NADP+ charge-transfer intermediate) for R200A FprA, albeit more slowly (k(lim)=58.1 s(-1) compared with >300 s(-1)) than in WT. All mutations caused positive shifts in FAD potential (approximately 40-65 mV). Binding of an NADPH analogue (tetrahydro-NADP) induced negative shifts in potential ( approximately 30-40 mV) only for variants with the R200A mutation, indicating distinctive effects of Arg(199)/Arg(200) on coenzyme binding mode and FAD potential. Collectively, these data reveal important roles for the phylogenetically conserved arginines in controlling FprA FAD environment, thermodynamics, coenzyme selectivity and reactivity.

Structural and Functional Characterization of the C-Terminal Domain of the Ecdysteroid Phosphate Phosphatase from Bombyx mori Reveals a New Enzymatic Activity

Here, we present the crystal structure of the ecdysone phosphate phosphatase (EPPase) phosphoglycerate mutase (PGM) homology domain, the first structure of a steroid phosphate phosphatase. The structure reveals an alpha/beta-fold common to members of the two histidine (2H)-phosphatase superfamily with strong homology to the Suppressor of T-cell receptor signaling-1 (Sts-1 PGM) protein. The putative EPPase PGM active site contains signature residues shared by 2H-phosphatase enzymes, including a conserved histidine (His80) that acts as a nucleophile during catalysis. The physiological substrate ecdysone 22-phosphate was modeled in a hydrophobic cavity close to the phosphate-binding site. EPPase PGM shows limited substrate specificity with an ability to hydrolyze steroid phosphates, the phospho-tyrosine (pTyr) substrate analogue para-nitrophenylphosphate ( pNPP) and pTyr-containing peptides and proteins. Altogether, our data demonstrate a new protein tyrosine phosphatase (PTP) activity for EPPase. They suggest that EPPase and its closest homologues can be grouped into a distinct subfamily in the large 2H-phosphatase superfamily of proteins.

Structure of l -Xylulose-5-Phosphate 3-Epimerase (UlaE) from the Anaerobic l -Ascorbate Utilization Pathway of Escherichia coli : Identification of a Novel Phosphate Binding Motif within a TIM Barrel Fold

Three catabolic enzymes, UlaD, UlaE, and UlaF, are involved in a pathway leading to fermentation of l-ascorbate under anaerobic conditions. UlaD catalyzes a beta-keto acid decarboxylation reaction to produce L-xylulose-5-phosphate, which undergoes successive epimerization reactions with UlaE (L-xylulose-5-phosphate 3-epimerase) and UlaF (L-ribulose-5-phosphate 4-epimerase), yielding D-xylulose-5-phosphate, an intermediate in the pentose phosphate pathway. We describe here crystallographic studies of UlaE from Escherichia coli O157:H7 that complete the structural characterization of this pathway. UlaE has a triosephosphate isomerase (TIM) barrel fold and forms dimers. The active site is located at the C-terminal ends of the parallel beta-strands. The enzyme binds Zn(2+), which is coordinated by Glu155, Asp185, His211, and Glu251. We identified a phosphate-binding site formed by residues from the beta1/alpha1 loop and alpha3' helix in the N-terminal region. This site differs from the well-characterized phosphate-binding motif found in several TIM barrel superfamilies that is located at strands beta7 and beta8. The intrinsic flexibility of the active site region is reflected by two different conformations of loops forming part of the substrate-binding site. Based on computational docking of the L-xylulose 5-phosphate substrate to UlaE and structural similarities of the active site of this enzyme to the active sites of other epimerases, a metal-dependent epimerization mechanism for UlaE is proposed, and Glu155 and Glu251 are implicated as catalytic residues. Mutation and activity measurements for structurally equivalent residues in related epimerases supported this mechanistic proposal.

Targeting the Ribose and Phosphate Binding Site of p38 Mitogen-Activated Protein (MAP) Kinase: Synthesis and Biological Testing of 2-Alkylsulfanyl-, 4(5)-Aryl-, 5(4)-Heteroaryl-Substituted Imidazoles

Three series of substituted 2-alkylsulfanyl-4-(4-fluorophenyl)imidazoles, 5-pyridinyl-, 1-methyl-5-pyridinyl-, and 5-(2-aminopyridin-4-yl)-imidazoles, were prepared and tested for their ability to inhibit p38 MAP kinase and TNF-alpha release. These compounds were prepared by using different synthetic routes. They were tested by applying a nonradioactive p38 MAP kinase assay and by measurement of TNF-alpha release in human whole blood. Potent inhibitors (IC50 values in the low nanomolar range, as low as 2 nM in the enzyme assay and 37 nM in the human whole blood test) were identified by variation of substituents at the imidazole- C2-thio position as well as at the 2-aminopyridinyl functionality. In contrast to other known kinase inhibitors, these novel imidazole derivatives with the substituents at the imidazole-C2-thio position may interact with the ribose as well as with the phosphate binding site of the p38 MAP kinase.

Utilization of an intron located polyadenlyation site resulted in four novel glutamate decarboxylase transcripts

Glutamate decarboxylase (GAD) is the rate-limiting enzyme in the synthesis of gamma-aminobutyric acid (GABA), the most important inhibitory neurotransmitter in central nervous system (CNS). Two homologous forms of GAD encoded by separate genes have been identified in mammalian brain, with molecular weight of 65 kDa (GAD65) and 67 kDa (GAD67). In the present study, four novel GAD67 transcripts produced by alternative splicing and polyadenlyation were cloned from rat testis. These novel GAD67 transcripts were widely expressed in non-neuronal tissues. During rat testis maturation, their expression level showed a time dependent change. These transcripts were predicted to synthesis of GAD proteins truncated of the binding site for pyridoxal phosphate, an essential cofactor, therefore cannot function as a decarboxylase. Thus, post-transcriptional processing mechanism as alternative splicing and polyadenlyation may play a crucial role in regulating rat GAD67 gene expression.

Are histones, tubulin, and actin derived from a common ancestral protein?

Histones and the cytoskeletal components tubulin and actin all act as thermal ratchets, using the energy present in Brownian motion to do work. All three also bind to nucleotides. Here we suggest that histones, tubulin, and actin derive from a common ancestral protein. There is some sequence similarity between histone 2A and the bacterial tubulin homologue FtsZ. Histones and actin also share some sequence similarity in the nucleotides and at phosphate-binding sites. Thus, actin and tubulin may also be related, although this is not obvious from sequence analysis. Indeed, actin and tubulin are closely functionally related and cooperate in many cellular processes. Interestingly, recent advances in nanotechnology suggest that thermal ratchets may be able to impart lifelike properties; thus, the evolution of the ancestral histone, tubulin, and actin thermal ratchet may have been crucial in the development of complexity in living organisms.

Carcinogenic metal compounds: recent insight into molecular and cellular mechanisms

Mechanisms of carcinogenicity are discussed for metals and their compounds, classified as carcinogenic to humans or considered to be carcinogenic to humans: arsenic, antimony, beryllium, cadmium, chromium, cobalt, lead, nickel and vanadium. Physicochemical properties govern uptake, intracellular distribution and binding of metal compounds. Interactions with proteins (e.g., with zinc finger structures) appear to be more relevant for metal carcinogenicity than binding to DNA. In general, metal genotoxicity is caused by indirect mechanisms. In spite of diverse physicochemical properties of metal compounds, three predominant mechanisms emerge: (1) interference with cellular redox regulation and induction of oxidative stress, which may cause oxidative DNA damage or trigger signaling cascades leading to stimulation of cell growth; (2) inhibition of major DNA repair systems resulting in genomic instability and accumulation of critical mutations; (3) deregulation of cell proliferation by induction of signaling pathways or inactivation of growth controls such as tumor suppressor genes. In addition, specific metal compounds exhibit unique mechanisms such as interruption of cell-cell adhesion by cadmium, direct DNA binding of trivalent chromium, and interaction of vanadate with phosphate binding sites of protein phosphatases.

Molecular Recognition of the Protein Phosphatase 1 Glycogen Targeting Subunit by Glycogen Phosphorylase

Disrupting the interaction between glycogen phosphorylase and the glycogen targeting subunit (G(L)) of protein phosphatase 1 is emerging as a novel target for the treatment of type 2 diabetes. To elucidate the molecular basis of binding, we have determined the crystal structure of liver phosphorylase bound to a G(L)-derived peptide. The structure reveals the C terminus of G(L) binding in a hydrophobically collapsed conformation to the allosteric regulator-binding site at the phosphorylase dimer interface. G(L) mimics interactions that are otherwise employed by the activator AMP. Functional studies show that G(L) binds tighter than AMP and confirm that the C-terminal Tyr-Tyr motif is the major determinant for G(L) binding potency. Our study validates the G(L)-phosphorylase interface as a novel target for small molecule interaction.

Functional polyterthiophene-appended uranyl-salophen complex: Electropolymerization and ion-selective response for monohydrogen phosphate

We have synthesized a bis(terthiophene)-appended uranyl-salophen complex, comprising N, N′-bis[4-(5,2′:5′,2″-terthiophen-3′-yl)salicylidene]-1,2-ethanediamine–uranyl complexes (TUS), and used it as a monomer for the electrochemical polymerizations (poly-TUS) on glassy carbon surfaces to prepare functionalized conducting polymer (CP) films. The poly-TUS films prepared from propylene carbonate/0.1 M tetrabutyl ammonium perchlorate (TBAP) on a glassy carbon electrode have both the functionality of ion-to-electron transducers (solid contact) and Lewis-acidic binding sites for a monohydrogen phosphate (MHP) ion-selective electrode (ISE). The CP/poly-TUS sensor showed a linear range between 1.0 × 10 −1 and 1.0 × 10 −4.5 M with a near-Nernstian behavior (−30.4 mV decade −1) at a pH of 8.2. The detection limit of the electrode was 10 −5.0 M and the response time was improved (<10 s) compared to that of conventional ISEs (<20 s). For comparison, a conventional ISE (with an internal aqueous solution) based on a TUS monomer/ o-nitrophenyl octylether ( o-NPOE)/polyvinyl chloride (PVC) liquid membrane with or without tridodecylmethylammonium chloride (TDMACl) as an additive was also constructed and its performance as an MHP-ISE were studied. The superior selectivity and sensitivity of the CP/poly-TUS sensor enabled the direct measurement of MHP in a wide variety of applications.

Structure and Catalytic Mechanism of Eukaryotic Selenocysteine Synthase

In eukaryotes and Archaea, selenocysteine synthase (SecS) converts O-phospho-L-seryl-tRNA [Ser]Sec into selenocysteyl-tRNA [Ser]Sec using selenophosphate as the selenium donor compound. The molecular mechanisms underlying SecS activity are presently unknown. We have delineated a 450-residue core of mouse SecS, which retained full selenocysteyl-tRNA [Ser]Sec synthesis activity, and determined its crystal structure at 1.65 A resolution. SecS exhibits three domains that place it in the fold type I family of pyridoxal phosphate (PLP)-dependent enzymes. Two SecS monomers interact intimately and together build up two identical active sites around PLP in a Schiff-base linkage with lysine 284. Two SecS dimers further associate to form a homotetramer. The N terminus, which mediates tetramer formation, and a large insertion that remodels the active site set SecS aside from other members of the family. The active site insertion contributes to PLP binding and positions a glutamate next to the PLP, where it could repel substrates with a free alpha-carboxyl group, suggesting why SecS does not act on free O-phospho-l-serine. Upon soaking crystals in phosphate buffer, a previously disordered loop within the active site insertion contracted to form a phosphate binding site. Residues that are strictly conserved in SecS orthologs but variant in related enzymes coordinate the phosphate and upon mutation corrupt SecS activity. Modeling suggested that the phosphate loop accommodates the gamma-phosphate moiety of O-phospho-l-seryl-tRNA [Ser]Sec and, after phosphate elimination, binds selenophosphate to initiate attack on the proposed aminoacrylyl-tRNA [Ser]Sec intermediate. Based on these results and on the activity profiles of mechanism-based inhibitors, we offer a detailed reaction mechanism for the enzyme.

Mutagenesis of the phosphate‐binding pocket of KDPG aldolase enhances selectivity for hydrophobic substrates

Narrow substrate specificities often limit the use of enzymes in biocatalysis. To further the development of Escherichia coli 2-keto-3-deoxy-6-phosphogluconate (KDPG) aldolase as a biocatalyst, the molecular determinants of substrate specificity were probed by mutagenesis. Our data demonstrate that S184 is located in the substrate-binding pocket and interacts with the phosphate moiety of KDPG, providing biochemical support for the binding model proposed on the basis of crystallographic data. An analysis of the substrate selectivity of the mutant enzymes indicates that alterations to the phosphate-binding site of KDPG aldolase changes the substrate selectivity. We report mutations that enhance catalysis of aldol cleavage of substrates lacking a phosphate moiety and demonstrate that electrophile reactivity correlates with the hydrophobicity of the substituted side chain. These mutations improve the selectivity for unnatural substrates as compared to KDPG by up to 2000-fold. Furthermore, the S184L KDPG aldolase mutant improves the catalytic efficiency for the synthesis of a precursor for nikkomycin by 40-fold, making it a useful biocatalyst for the preparation of fine chemicals.

The structure/function relationship of a dual-substrate (βα) 8-isomerase

Two structures of phosphoribosyl isomerase A (PriA) from Streptomyces coelicolor, involved in both histidine and tryptophan biosynthesis, were solved at 1.8 Å resolution. A closed conformer was obtained, which represents the first complete structure of PriA, revealing hitherto unnoticed molecular interactions and the occurrence of conformational changes. Inspection of these conformers, including ligand-docking simulations, allowed identification of residues involved in substrate recognition, chemical catalysis and conformational changes. These predictions were validated by mutagenesis and functional analysis. Arg 19 and Ser 81 were shown to play critical roles within the carboxyl and amino phosphate-binding sites, respectively; the catalytic residues Asp 11 and Asp 130 are responsible for both activities; and Thr 166 and Asp 171, which make an unusual contact, are likely to elicit the conformational changes needed for adopting the active site architectures. This represents the first report of the structure/function relationship of this (βα) 8-isomerase.

Stereospecific Proton Transfer by a Mobile Catalyst in Mammalian Fructose-1,6-bisphosphate Aldolase

Class I fructose-1,6-bisphosphate aldolases catalyze the interconversion between the enamine and iminium covalent enzymatic intermediates by stereospecific exchange of the pro(S) proton of the dihydroxyacetone-phosphate C3 carbon, an obligatory reaction step during substrate cleavage. To investigate the mechanism of stereospecific proton exchange, high resolution crystal structures of native and a mutant Lys(146) --> Met aldolase were solved in complex with dihydroxyacetone phosphate. The structural analysis revealed trapping of the enamine intermediate at Lys(229) in native aldolase. Mutation of conserved active site residue Lys(146) to Met drastically decreased activity and enabled trapping of the putative iminium intermediate in the crystal structure showing active site attachment by C-terminal residues 360-363. Attachment positions the conserved C-terminal Tyr(363) hydroxyl within 2.9A of the C3 carbon in the iminium in an orientation consistent with incipient re face proton transfer. We propose a catalytic mechanism by which the mobile C-terminal Tyr(363) is activated by the iminium phosphate via a structurally conserved water molecule to yield a transient phenate, whose developing negative charge is stabilized by a Lys(146) positive charge, and which abstracts the C3 pro(S) proton forming the enamine. An identical C-terminal binding mode observed in the presence of phosphate in the native structure corroborates Tyr(363) interaction with Lys(146) and is consistent with transient C terminus binding in the enamine. The absence of charge stabilization and of a mobile C-terminal catalyst explains the extraordinary stability of enamine intermediates in transaldolases.

Computational Study of Glyceraldehyde-3-phosphate Dehydrogenase of Entamoeba histolytica: Implications for Structure-Based Drug Design

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of the pathogenic protozoa Entamoeba histolytica (Eh) is a major glycolytic enzyme and an attractive drug target since this parasite lacks a functional citric acid cycle and is dependent solely on glycolysis for its energy requirements. The three-dimensional structure of dimeric EhGAPDH in complex with cofactor NAD+ has been generated by homology modeling based on the crystal structure of human liver GAPDH. Our refined model indicates the presence of a parasite specific disulfide bond between two cysteine residues of adjacent monomers in the EhGAPDH dimer, which may be an important target for future drug design. Flexible docking with the substrate glyceraldehyde-3-phosphate (G3P) shows that Cys151, His178, Thr210, and Arg233 are important residues in G3P binding. The inorganic phosphate-binding site of EhGAPDH has been determined by docking study. The binding mode of a natural GAPDH inhibitor, chalepin to EhGAPDH has also been predicted. In search for a better inhibitor for EhGADPH, in silico modification of chalepin has been carried out to form an additional specific polar interaction with Asp194 of EhGAPDH whose equivalent is Leu195 in human GAPDH. In the absence of a crystal structure, our study provides an early insight into the structure of major drug target EhGAPDH, thus, facilitating the inhibitor design.

Common conformational changes in flavodoxins induced by FMN and anion binding: The structure of Helicobacter pylori apoflavodoxin

Flavodoxins, noncovalent complexes between apoflavodoxins and flavin mononucleotide (FMN), are useful models to investigate the mechanism of protein/flavin recognition. In this respect, the only available crystal structure of an apoflavodoxin (that from Anabaena) showed a closed isoalloxazine pocket and the presence of a bound phosphate ion, which posed many questions on the recognition mechanism and on the potential physiological role exerted by phosphate ions. To address these issues we report here the X-ray structure of the apoflavodoxin from the pathogen Helicobacter pylori. The protein naturally lacks one of the conserved aromatic residues that close the isoalloxazine pocket in Anabaena, and the structure has been determined in a medium lacking phosphate. In spite of these significant differences, the isoallozaxine pocket in H. pylori apoflavodoxin appears also closed and a chloride ion is bound at a native-like FMN phosphate site. It seems thus that it is a general characteristic of apoflavodoxins to display closed, non-native, isoalloxazine binding sites together with native-like, rather promiscuous, phosphate binding sites that can bear other available small anions present in solution. In this respect, both binding energy hot spots of the apoflavodoxin/FMN complex are initially unavailable to FMN binding and the specific spot for FMN recognition may depend on the dynamics of the two candidate regions. Molecular dynamics simulations show that the isoalloxazine binding loops are intrinsically flexible at physiological temperatures, thus facilitating the intercalation of the cofactor, and that their mobility is modulated by the anion bound at the phosphate site.

Characterization and Expression Pattern of 1-Aminocyclopropane-1-Carboxylate (ACC) Synthases cDNA from Geranium (Pelargonium x hortorum cv Sincerity)

The enzyme 1-aminocyclopropane-1-caryboxylate (ACC) synthase catalyzes the rate limiting step in the regulation of ethylene biosynthesis and is, primarily, regulated at the transcription level. To understand the role of ACC synthases in stress-induced ethylene bursts in geranium, degenerated oligonucleotide primers were used to amplify a DNA probe from cDNA; and used to screen a cDNA library prepared from an mRNA preparation from mature flowers. Six putative full-length ACC synthase clones were isolated; these clones were divided into two groups on the basis of deduced amino acid sequences and features contained in their mRNA. Clones (e.g. PHS-ACS41, 44, and 45) in group I contain a stop codon (TAG) just before the initiation codon (ATG). Clone PHS-ACS44 from this group contains an “extra” 780 nucleotide sequence at the untranslated 3′-end region (UTR) after the first polyadenylation signal and a short poly(A) tail. The “extra” UTR contains another polyadenylation signal and a poly(A) tail, otherwise the coding sequence is identical to PHS-ACS41. Clone PHS-ACS49 in group II contains stop codon (TAG) twenty one bases upstream from the initiation codon. All these genes contain eleven conserved amino acid residues, the substrate and pyridoxal 5′-phosphate binding site that are characteristic of all ACC synthases. Using RTPCR and gene specific primers, expression studies show that these genes are expressed differentially. PHS-ACS41 and PHSACS44 not only show differences in expression in vivo, but also exhibit differential translational efficiencies in vitro suggesting that the “extra” 780 UTR in vivo may serve as a regulatory element. Our finding of multiple ACC synthase genes and their expression in geranium is consistent with similar observations in other plant species in which ACC synthases are encoded by a multigene family that are differentially regulated by internal and external factors.

Structure–function relationships in a bacterial DING protein

A recombinant DING protein from Pseudomonas fluorescens has been previously shown to have a phosphate-binding site, and to be mitogenic for human cells. Here we report the three-dimensional structure of the protein, confirming a close similarity to the “Venus flytrap” structure seen in other human and bacterial phosphate-binding proteins. Site-directed mutagenesis confirms the role of a key residue involved in phosphate binding, and that the mitogenic activity is not dependent on this property. Deletion of one of the two hinged domains that constitute the Venus flytrap also eliminates phosphate binding whilst enhancing mitogenic activity.

Hemoglobin structure/function and globin-gene evolution in the Arctic fish Liparis tunicatus

The importance of the Arctic, in contributing to the knowledge of the overall ensemble of adaptive processes influencing the evolution of marine organisms, calls for investigations on molecular adaptations in Arctic fish. Unlike the vast majority of Antarctic Notothenioidei, several Arctic species display high hemoglobin multiplicity. The blood of four species, the spotted wolffish of the family Anarhichadidae and three Gadidae, contains three functionally distinct major components. Similar to many Antarctic notothenioids, Arctic Liparis tunicatus (suborder Cottoidei, family Liparidae) has one major hemoglobin (Hb 1) accompanied by a minor component (Hb 2). This paper reports the structural and functional characterisation of Hb 1 of L. tunicatus. This hemoglobin shows low oxygen affinity, and pronounced Bohr and Root effects. The amino-acid sequence of the β chain displays an unusual substitution in NA2 (β2) at the phosphate-binding site, and the replacement of Val E11 (β67) with Ile. Similar to some Antarctic fish Hbs, electron paramagnetic resonance spectra reveal the formation of a ferric penta-coordinated species even at physiological pH. The amino-acid sequences have also been used to gain insight into the evolutionary history of globins of polar fish. L. tunicatus globins appear close to the notothenioid clades as predicted by teleostean phylogenies. Close phylogenetic relationships between Cottoidei and Notothenioidei, together with their life style, seem to be the main factor driving the globin-sequence evolution.