Born to sense: biophysical analyses of the oxygen sensing prolyl hydroxylase from the simplest animal Trichoplax adhaerens

Background In humans and other animals, the chronic hypoxic response is mediated by hypoxia inducible transcription factors (HIFs) which regulate the expression of genes that counteract the effects of limiting oxygen. Prolyl hydroxylases (PHDs) act as hypoxia sensors for the HIF system in organisms ranging from humans to the simplest animal Trichoplax adhaerens. Methods We report structural and biochemical studies on the T. adhaerens HIF prolyl hydroxylase (TaPHD) that inform about the evolution of hypoxia sensing in animals. Results High resolution crystal structures (≤1.3 Å) of TaPHD, with and without its HIFα substrate, reveal remarkable conservation of key active site elements between T. adhaerens and human PHDs, which also manifest in kinetic comparisons. Conclusion Conserved structural features of TaPHD and human PHDs include those apparently enabling the slow binding/reaction of oxygen with the active site Fe(II), the formation of a stable 2-oxoglutarate complex, and a stereoelectronically promoted change in conformation of the hydroxylated proline-residue. Comparison of substrate selectivity between the human PHDs and TaPHD provides insights into the selectivity determinants of HIF binding by the PHDs, and into the evolution of the multiple HIFs and PHDs present in higher animals.


Introduction
Hypoxia inducible transcription factors (HIFs) play key roles in maintaining oxygen homeostasis in most, if not all, animals. [1][2][3][4] Extensive work on the HIF mediated response to hypoxia in humans and other higher animals ( Figure 1A and B) has revealed that the levels and transcriptional activity of the α,β-heterodimeric HIF complex are regulated by ferrous iron and 2-oxoglutarate (2OG) dependent oxygenases. [4][5][6] Prolyl-4-hydroxylation of the oxygen dependent degradation domains (ODDs) in HIFα is catalyzed by the HIF prolyl hydroxylases (PHDs or EGLNs), and signals for HIFα degradation via the ubiquitin proteasome system. HIFα prolyl hydroxylation strongly promotes its binding to the von-Hippel Lindau protein (pVHL), a targeting component of an E3 ubiquitin ligase complex 7-9 ( Figure 1C). In some animals, a second 2OG oxygenase, factor inhibiting HIF (FIH), catalyzes asparaginyl-hydroxylation in the C-terminal transcriptional activation domain (CTAD) present in both HIF1α and HIF2α, but not in HIF3α; 10,11 this modification inhibits the interaction between HIF and transcriptional coactivators, which are histone acetyl transferases (CBP/p300). 11 FIH is active at lower oxygen concentrations than the PHDs 12,13 and is likely less important from a fundamental hypoxia-sensing perspective.
The PHDs and FIH belong to structurally distinct 2OG oxygenase subfamilies, with FIH being part of the Jumonji-C (JmjC) 2OG oxygenase subfamily IX, 14 and the PHDs belonging to the PH (prolyl hydroxylase) subfamily VIII. 15 However, unlike the HIFα-PHD-pVHL "hypoxia sensing triad", the FIH-HIFα CTAD "diad" is not universally conserved, being only sporadically present in simple animals, implying it evolved subsequent to the PHD enabled sensing mechanism. 16 Extensive genetic and chemical intervention studies support the physiological importance of the HIFα-PHD-pVHL triad in animals. Indeed, the human PHDs are current therapeutic targets, with PHD inhibitors for the treatment of anemia, via HIF mediated upregulation of erythropoietin, being in late stage clinical trials. 17,18 In humans, there are three PHDs (PHD1-3) and three HIFα paralogues, but only one FIH. Animal studies imply that PHD2, which, unlike PHD1 and 3, contains a highly conserved Myeloid, Nervy, and DEAF-1 (MYND) finger domain, is the most important of the three mammalian PHDs from a physiological perspective. 19,20 Unlike humans and other mammals, there is only one PHD present in (most) lower organisms, which is more similar to PHD2 than PHD1 or PHD3. 16 The human HIF system regulates hundreds of different genes in a context and cell-type dependent manner, with emerging evidence suggesting different roles for HIF1α and HIF2α and the three PHDs. 12,[21][22][23] HIF1α and HIF2α each have two ODDs, the C-terminal and N-terminal (CODD and NODD), respectively, which manifest differential sensitivities to hypoxia. 12,24 In contrast to the extensive studies conducted on HIF1α and HIF2α, the biological role of HIF3α, which contains a single ODD, is less understood, with a role for HIF3α splice variants in the regulation of the hypoxic response being proposed. [25][26][27] The three human PHDs also have different ODD selectivities, with PHD1 and PHD2 accepting both NODD and CODD as substrates; however, PHD3 has a strong preference Figure 1 Overview of the HiF system. Notes: (A) comparison of domain architectures of HiFα and HiFβ in T. adhaerens and humans. in contrast to the multiple PHD/HiFα-isoforms and ODDs present in vertebrates, the prime components of the HiF system in T. adhaerens involve only one PHD and one HiFα, which has a single ODD. 16 arrows indicate assigned prolyl and asparaginyl hydroxylation sites. (B) Domain structures of HiF prolyl hydroxylases (PHDs) in T. adhaerens and humans. Domain acronyms: 2Og dioxygenase domain (oxygenase), MYeloid, Nervy, and DEAF-1 (MYND)-type zinc finger domain (MYND). (C) Outline of the conserved mechanism of the response to chronic hypoxia in animals. The HiF transcription factors are regulated by PHD catalyzed hydroxylation of prolyl residues in the ODD(s) of HiFα under normoxia. recognition of the hydroxylated prolyl residues by pVHl is followed by ubiquitination by the e3 ubiquitin ligase, which tags HiFα for proteasomal degradation. in humans, FiH, which is only sporadically present in non-vertebrate animals, 16 constitutes an additional oxygen dependent regulatory element of HiF activity. FiH catalyzes the hydroxylation of an asparagine residue in the cTaD of HiFα, thus disrupting its interaction with the cBP/p300 transcriptional co-activator, and hindering transcriptional activation. 11 Abbreviations: PHD, Prolyl hydroxylases; HiF, hypoxia-inducible transcription factor; ODD, oxygen-dependent degradation domain; bHlH, basic helix-loop-helix motif; Pas, Per-arnT-sim domain; cTaD, C-terminal transactivation domain; Pac, Pas-associated C-terminal domain; FiH, factor inhibiting HiF.

59
Biophysical analyses of the oxygen sensing PHDs from the simplest animal T. adhaerens for CODD over NODD as a substrate. 21,22 Biochemical studies on the PHDs are supportive of their roles as hypoxic sensors. 21,[28][29][30] An essential requirement for the proposed hypoxia sensing role of the PHDs is that their activity in cells is regulated by oxygen availability. Studies, both in cells and animals, support the proposed hypoxia sensing role of the PHDs, 1-9 raising the question as to whether their hypoxia sensing role is manifested in specific biochemical properties. There is potential for various factors to help enable PHD activity to be limited by oxygen availability in cells, including regulation of localization, concentrations, or activities of the PHDs and HIFα isoforms. Several studies support the proposal that the kinetic properties of the PHDs are directly related to their sensing role. 21,[28][29][30] Strikingly, while PHD2, the most conserved of the human PHDs, binds Fe(II) and 2OG relatively tightly compared to other 2OG oxygenases, forming an unusually stable enzyme.Fe.2OG complex, the PHD2.Fe.2OG.ODD complex reacts unusually slowly with oxygen ( Figure S1). 30,31 Structural analyses on PHD2 have led to the proposal that its slow reaction with oxygen is related to the active site coordination chemistry, with substitution of an Fe(II) bound water by dioxygen being a limiting step in catalysis. [32][33][34][35][36] To investigate the proposal that the apparently special kinetic properties of PHDs are a general feature of their hypoxia sensing roles, we have initiated studies on their evolution. Bioinformatic studies suggest that the HIFα-PHD-pVHL triad only exists in animals, though homologues of the PHDs, but not apparently HIF, are present in earlier organisms. 16,37,38 Roles for PHD homologues in early eukaryotes, ie, hydroxylation of S-phase kinase associated protein 1 (Skp1) in Dictyostelium, 37 and in some bacteria, ie, hydroxylation of the Elongation Factor EF-Tu 38 have been identified. However, while these enzymes also catalyze prolyl C-4 hydroxylation and their overall structures are clearly related to the PHDs, the available evidence is that these enzymes differ in their kinetic properties relative to the human PHDs, [37][38][39] as is the case for the procollagen prolyl-4-hydroxylases. 21,28,[40][41][42][43][44][45][46] The extent to which there is structural conservation of the active sites between the PHDs from different animals has been unclear to date; we have reported experimental evidence for the presence of a functional HIFα-PHD-pVHL system in Trichoplax adhaerens, the simplest known animal ( Figure S2). 16 Interestingly, the genes encoding for TaHIFα and TaPHD are adjacent in the T. adhaerens genome, suggesting linked evolution, and TaPHD activity can be limited by oxygen availability in human cells. 16 Here we report crystallographic and biochemical studies on the single PHD present in T. adhaerens (TaPHD, Uniprot ID: I6QVT6), the results of which reveal striking elements of structural and biochemical conservation between TaPHD and human PHD2 (HsPHD2) active sites, supporting oxygen binding/reaction as a key property in the oxygen sensing capability of the PHDs. The results also provide insights into the evolution of ODD binding by PHDs, and help describe how human PHDs achieve ODD selectivity.

Results
steady-state kinetic studies of TaPHD in comparison to HsPHD2 The HIF system in T. adhaerens is apparently markedly simpler than in humans because it contains only one PHD (TaPHD), and one HIFα isoform with a single ODD (TaODD, Figure 1A). 16 Although the sequences of human HIF1α CODD and NODD (HsHIF1α CODD/NODD), and TaODD differ substantially (Figure 2A), N-terminally truncated TaPHD (aa 64-300, TaPHD 64-300 ) catalyzes prolyl hydroxylation of both the TaODD 25mer and the human HIF1α CODD and HIF1α NODD 19mer substrates, with the latter being hydroxylated less efficiently. 16 To further investigate the substrate preferences of TaPHD 64-300 , its affinities for TaODD 25mer, human HIF sequences, and its co-substrate (2OG) were determined in steady-state kinetic studies and compared to the truncated human PHD isoform 2 (HsPHD2 181-426 ). The kinetic parameters for the catalytic domain of HsPHD2 181-426 have been shown to be similar to those of the full-length protein. 31,59 K m and k cat values ( Figure 2B) for HIF peptide substrates and 2OG were determined under optimized assay conditions described in "Experimental Methods" in Supplementary materials, with the extent of peptide hydroxylation being analyzed by Matrix-Assisted Laser Desorption/Ionization (MALDI) time-of-flight mass spectrometry (MS) ( Figures  2B and S3 -5). Consistent with previous reports, 16 TaPHD 64-300 efficiently catalyzes the hydroxylation of its natural substrate TaODD 25mer and the human HIF1α CODD 19mer peptide. The HsHIF1α NODD 19mer, by contrast, is not turned over above background level by TaPHD 64-300 , at least under the conditions/limits of detection of the assays employed here ( Figure S5). TaPHD 64-300 has a similar affinity to both these substrates, despite their differences in length, as judged by    TaPHD binds HsHiF1α nODD less efficiently than HsHiF1α cODD The observed preference of TaPHD 64-300 for HsHIF1α CODD over HsHIF1α NODD is in agreement with the finding of one ODD in the genome of T. adhaerens, which is closer in sequence to the majority of CODD, compared to NODD sequences ( Figure 2A). 16 While a NODD-like sequence is present in all vertebrates, NODD is apparently not present in some non-vertebrates, eg, Placozoa and Cnidaria, suggesting that the CODD sequence might have evolved prior to that of NODD. 16 Notably, like TaPHD 64-300 , one of the human PHD isoforms, HsPHD3, is highly selective for HsHIF1α/2α CODD over HsHIF1α/2α NODD. 1,21,47 To more directly investigate the extent to which TaPHD 64-300 can bind HsHIF1α NODD and CODD, we performed 1D CLIP HSQC NMR binding studies with HsHIF1α CODD and HsHIF1α NODD peptides labeled with [ 13 C] at the target prolyl residues which undergo hydroxylation during catalysis ( Figure 2C). The [ 13 C]-labeled "reporter" HsHIF1α CODD/ NODD peptides were incubated with Zn(II) and 2OG in the presence and absence of TaPHD 64-300 (1:1 enzyme:ODD, for conditions see "Experimental methods" in Supplementary materials). The decrease of the reporter signal upon binding to the enzyme was measured to investigate the relative binding strength ( Figure 2C). Consistent with the kinetic studies, the NMR results reveal that TaPHD 64-300 strongly discriminates between binding of HsHIF1α CODD and HsHIF1α NODD, as ~40% of the HsHIF1α NODD peptide remained unbound, while HsHIF1α CODD was apparently completely bound to TaPHD  , showing that TaPHD 64-300 binds HsHIF1α CODD better than HsHIF1α NODD.
Monitoring the coupled and uncoupled 2Og turnover of TaPHD 2OG dependent oxygenases catalyze the uncoupled turnover of 2OG to succinate and CO 2 in the absence of their prime substrate to varying extents; 48 in the case of HsPHD2, it has been shown that 2OG decarboxylation is strongly coupled to CODD hydroxylation. 31 To directly compare the extent of uncoupled turnover in TaPHD 64-300 and HsPHD2 181-426 , succinate formation and 2OG depletion were monitored by 1 H-NMR 49 in the presence and absence of their respective natural substrates ( Figure 2D and E). The results demonstrate that the rate of uncoupled 2OG turnover is very low for both TaPHD 64-300 and HsPHD2 181-426 ; 2OG decarboxylation is strongly increased in presence of substrate, implying conservation of (at least some) kinetic properties between HsPHD2 and TaPHD. Overall, these results support the proposal that the strong binding of Fe(II) and 2OG, and the formation of a stable HsPHD2.Fe.2OG complex in the absence of HsHIFα, are conserved properties of the PHDs.

Oxygen dependence of TaPHD
In order for the PHDs to act as hypoxia sensors in cells, their activity must be limited by oxygen availability. This role is proposed to be reflected in the slow reaction of (at least) the HsPHD2.CODD/NODD.Fe.2OG complexes with oxygen. 29,31 It is also proposed to manifest in the unusually high K m value (for 2OG oxygenases) of HsPHD2 for oxygen (K m (O 2 , HsPHD2) >400 µM). 21,28,29,50 In order to study the O 2 dependence of the TaPHD-catalyzed reaction, steady-state kinetics, using reported experimental conditions for HsPHD2; 29 were conducted. TaPHD 64-300 was reacted with its peptidic substrate in sealed glass-vials under varied oxygen concentrations, 29 the initial rates of TaODD 25mer hydroxylation were determined by MALDI-TOF-MS ( Figure 2F and G). No saturation of O 2 was reached for the TaPHD 64-300 -catalyzed reaction with TaODD 25mer at 60% O 2 level, resulting in a K m (O 2 , TaPHD 64-300 ) >400 µM. This result suggests that TaPHD has the potential to react to small changes in atmospheric O 2 partial pressure and could act as an effective oxygen sensor, indicating that the core biochemical features of HsPHD2 in the human oxygen sensing machinery may well be conserved in TaPHD. It has been reported that HsPHD2 reacts relatively slowly with oxygen, compared to other 2OG oxygenases as submit your manuscript | www.dovepress.com Dovepress Dovepress 62 lippl et al determined by stopped flow/UV-visible spectroscopy experiments. 29,31 It was not possible to conduct a similar study on TaPHD for practical reasons (ie, neither TaPHD 64-300 , nor other tested constructs TaPHD 1-300 and TaPHD 79-300 could reach a sufficiently high protein concentration necessary for stopped-flow assays due to aggregation). crystallization of N-terminally truncated T. adhaerens PHD (TaPHD  In order to investigate the extent to which the structural features proposed to be necessary for oxygen sensing are conserved between HsPHD2 and TaPHD, we attempted to obtain a crystal structure of TaPHD in complex with Mn(II), a non-reactive substitute for Fe(II), using a TaPHD construct (aa 64-300, TaPHD 64-300 , MW=29.2 kDa) lacking the N-terminal MYND domain, which could be produced efficiently in Escherichia coli. The shorter construct (TaPHD 64-300 ) yielded crystals that diffracted to 1.2 Å resolution ( Table 1). The structure was solved by molecular replacement, using HsPHD2 (PDB: 2G19), 32 as a search model. To gain insight into the substrate-binding mode of TaPHD, we then co-crystallized TaPHD 64-300 with a 21mer peptide fragment of its substrate, the TaHIFα ODD (TaODD, residues E477-L497), Mn(II) and N-oxalylglycine (NOG), a close 2OG analog. The TaPHD.TaODD.Mn.NOG structure (Table 1) was solved to a resolution of 1.3 Å by molecular replacement, using the structure without TaODD as a search model. The substrate-bound and unbound forms of TaPHD crystallized in the P1 and P2 1 space groups, respectively, reflecting different packing, possibly relating to structural changes in TaPHD induced upon substrate binding.

63
Biophysical analyses of the oxygen sensing PHDs from the simplest animal T. adhaerens Overall structure of TaPHD The structure of TaPHD ( Figure 3A) contains the conserved double-stranded β-helix core-fold (DSBH, or "jelly roll"motif), which is present in all characterized 2OG dependent oxygenases, including the human PHDs. 15 The DSBH is comprised of the major (βI, βVIII, βIII, βVI) and minor (βII, βVII, βIV, βV) β-sheets between which the metal and 2OG binding sites are sandwiched. As with HsPHD2, the major β-sheet core is stabilized by three α-helices (α1-3).

active site of TaPHD
The active site pocket of TaPHD, positioned at one end of the DSBH, deeply embeds the metal ion (Mn, substituting for Fe) and the co-substrate (NOG, substituting for 2OG, Figure 3B and C). The highly buried nature of the 2OG binding site likely contributes to the formation of a stable TaPHD.Fe.2OG complex, as reflected in the observed low rate of uncoupled 2OG turnover by TaPHD (Figures 2D  and 3B).
In the TaPHD.Mn structure without TaODD, the manganese ion is octahedrally coordinated by a conserved triad of residues (H209 TaPHD , D211 TaPHD , and H270 TaPHD ), an acetate ion, and two water molecules. 48,51 Deep within the 2OG binding pocket an additional acetate ion forms a salt bridge with R279 TaPHD , apparently mimicking 2OG C5 carboxylate binding ( Figure 3A).
The 6-coordinate state of the metal is retained in the TaPHD.TaODD complex. Metal-ion binding involves the highly conserved triad of metal-coordinating residues, and an ordered, metal-bound water molecule (W1) that is, in addition to the metal coordination, stabilized via hydrogen bonding with the metal-ligating side chain of D211 TaPHD ( Figures 3C and 6). Electron density corresponding to the 2OG analog inhibitor, NOG, is clearly visible in the TaPHD. TaODD complex structure. NOG coordinates to the metal ion in a bidentate manner and interacts with R279 TaPHD by The metal ion (manganese substituting for iron, purple sphere) is octahedrally coordinated by a triad of residues (H209 TaPHD , D211 TaPHD , and H270 TaPHD ), N-oxalylglycine (nOg), and a water molecule (W1, red sphere). The stable metal-water coordination observed here is conserved in HsPHD2, 32,34 where it is proposed to enable the oxygen sensing ability in HsPHD2. its C5 carboxylate. The 1-carboxylate of NOG is positioned at the coordination site closest to the hydroxylation target, P486 TaODD , making the remaining metal coordination site, occupied by water W1, less accessible, and thus apparently hindering oxygen binding to the metal ion ( Figure 6). 33,34 Crystal structures of the human PHD isoform 2 (HsPHD2) in complex with its substrates HsHIF1α CODD (referred to as HsPHD2.CODD, PDB: 3HQR 34 ) and HsHIF1α NODD (referred to as HsPHD2.NODD, PDB: 5L9V 33 ) have been reported (Figure 4). The overall active site arrangement observed in the TaPHD.TaODD complex is very similar to those observed in the HsPHD2.CODD/NODD complexes (Figures 4-6, 33,34 ). Such conservation includes the octahedral coordination of the metal ion by H313 HsPHD2 , D315 HsPHD2 , and H374 HsPHD2 , NOG, and a metal-bound water molecule. 33,34,52 The observation of a well-ordered metal bound water in both TaPHD and HsPHD2 structures is important from a hypoxia sensing perspective, because displacement of this water from the metal is required for oxygen to bind, and is proposed to be rate limiting in HsPHD2 catalysis. 32,34,53 The buried nature of the metal, the 2OG analog NOG, and the associated ligating water molecule, as observed both in TaPHD and HsPHD2, are likely substantially responsible for the unusual stability of the HsPHD2.Fe.2OG complex, 30 compared to other 2OG oxygenases such as FIH, which is active at lower oxygen concentrations. 12,13,54 Taken together, the observation of conserved active site arrangements in the TaPHD.TaODD and HsPHD2.CODD crystal structures correlates well with the steady-state kinetic parameters determined for 2OG and oxygen, imply-  comparison of ODD binding modes for TaPHD and HsPHD2 The main substrate binding elements in the TaPHD.TaODD complex structure, as in HsPHD2, comprise: the active site containing groove, the β2/β3-finger-loop, and the C-terminal α4-helix ( Figure 3B). The TaODD peptide binds to TaPHD in an extended form; notably, the residues N-terminal to the target proline form a 3 10 -helix. The target proline in TaODD and its N-terminally flanking residues form a YXXLAP motif, which differs by one residue (L->Y) compared to the consensus LXXLAP prolyl-hydroxylation site(s) present in all human, and many, but not all animal, HIFα proteins. 16  [55][56][57]. in (A and D), N-oxalylglycine (nOg) acts as a 2Og analog. note that the Zn (substituting for Fe) in the CrP4H active site is tetrahedrally coordinated, with an acetate binding instead of the 2Og co-substrate. The structures reveal conserved orientations of the "target" proline residues, which in each case, adopt a c-4 endo-conformation. 55
As observed in the HsPHD2.CODD/NODD structures, the substrate proline P486 TaODD pyrrolidine ring adopts a C-4 endo-conformation, with the 4R C-H bond that is cleaved during hydroxylation, being positioned close to the metal center (distance: 4.5 Å, Figures 3C and 6, 33,34 ); indeed, the conformations of the substrate target prolyl rings at the TaPHD and HsPHD2 active sites are nearly identical ( Figure 4A and B). Upon hydroxylation by HsPHD2, the conformation of the pyrrolidine ring in CODD is proposed to switch from C-4 endo to C-4 exo 55 ( Figure 6). When bound to pVHL, the C-4 hydroxylated proline adopts the C-4 exo conformation as favored by stereoelectronic theory (the gauche effect); 55-58 a role for this conformational change in promoting product release from HsPHD2 has been suggested. 34 In the human PHDs, the β2/β3-finger-loop plays a significant role in substrate recognition and in determining NODD/ CODD selectivity. 33,34,59 The β2/β3-finger-loop is partially disordered without its substrate bound, likely reflecting flexibility relative to the enzyme-substrate complex. 33 By contrast, in both the HsPHD2.CODD and HsPHD2.NODD structures, the β2/β3-loop residues fold to enclose the LXXLAP substrate motif at the active site ( Figure 4A and  B). Although, the β2/β3-finger-loop in the TaPHD.TaODD complex also constitutes one of the major substrate interaction sites, it is two residues shorter than in HsPHD2 between N141 TaPHD -V142 TaPHD (Figures S4 and 6). As a consequence, and by contrast to HsPHD2, 33 the β2/β3finger-loop in the TaPHD.TaODD structure is more condensed and stabilized by intramolecular interactions between loop residues, including backbone-backbone H-bonds and salt bridges between side chains ( Figure S3). This structural difference might reflect the fact that there is only a single HIFα isoform with one ODD in T. adhaerens, which contrasts with the situation in humans where there are three HIFα isoform substrates; with both a NODD and a CODD in case of HIF1α and HIF2α (note other HsPHD substrates have also been reported). [60][61][62][63][64][65] It is thus possible that binding of multiple HIF substrates to HsPHD2 33,34 requires a higher degree of flexibility in the β2/ β3-finger-loop, compared to the apparently more rigid β2/ β3-finger-loop in TaPHD.
In both the TaPHD.TaODD and HsPHD2.CODD 34 structures, a salt-bridge between an arginine residue in the C-terminal α4-helix and an aspartate residue in the C-terminus of the respective peptide (R292 TaPHD /D494 TaODD, R396 HsPHD2 / D571 HsHIF1α CODD ) is formed ( Figure S4). However, this saltbridge is not present in the HsPHD2.NODD complex structure, where the HsHIF1α NODD peptide is positioned further from the C-terminal helix and interacts with the C-terminal region of HsPHD2 via hydrophobic contacts. Studies with the TaODD and HsPHD2 variants, including with HsPHD2 P317R HsPHD2 and R396T HsPHD2 , which are selective for CODD or NODD, respectively, support the proposal of a similar binding mode for TaODD and HsHIF1α CODD/NODD to HsPHD2 ( Figure S4).
By contrast with HsHIF1α CODD/NODD peptides, the TaODD adopts a short 3 10 -helix near its Pro-Pro-Pro motif (P489 TaODD -491 TaODD ), which is C-terminal to the target proline P486 TaODD residue ( Figure 4A and B). Notably, the Pro-Pro-Pro motif and the adjacent 3 10 -helix constitute the interaction surface with the C-terminal α4-helix of TaPHD ( Figure 4A and B). In the TaPHD.TaODD complex structure, the TaODD is positioned particularly close to the α4-helix of TaPHD, compared to both HsHIF1α CODD and HsPHD2.NODD, whereas the latter is positioned even more distant from the C-terminus of HsPHD2 ( Figure 4A and B).
In addition to their overall structural conservation, the prolyl-4-hydroxylases adopt strikingly similar substratebinding features (Figures 4-6 and S5). The TaODD, HsHIF1α CODD/NODD, EF-Tu, and the proline-rich (Ser-Pro) 5 substrate backbones are all oriented in the same N→C direction across the active site ( Figure 5A and B). TaODD, HsHIF1α CODD, and HsHIF1α NODD all form a partial helical structure (3 10 -helix), when bound to enzyme ( Figure 4A and B), while isolated HsHIF1α CODD/NODD are predicted to be disordered in solution. 33,34 By contrast, the EF-Tu switch loop retains loop secondary structure upon binding to PPHD, albeit after a large conformational change. 38 Notably, the PHDs all employ a flexible β2/β3-finger-loop (G136 TaPHD -D150 TaPHD , S77 CrP4H -S95 CrP4H , A53 PPHD -D70 PPHD ) that undergoes a major conformational change upon substrate binding and helps to enclose the target proline of the substrate in the active site cleft (Figures 5, 6 and S4). 33,34 In TaPHD, HsPHD2 and PPHD, the C-terminus interacts with the substrate; it is unknown if this is the case for CrP4H, since the shorter peptide substrate used in the CrP4H.(Ser-Pro) 5 structure does not appear to reach the C-terminus of the enzyme. 66 conservation of proline-ring conformation Comparison of the binding modes of the prolyl-4-hydroxylase substrates in the immediate active sites reveals a high degree of conservation ( Figures 5A, B and 6). The substrates all manifest a similar orientation toward the active site residues and the metal center, and the target prolyl residues all adopt the C-4 endo-conformation ( Figure 5A and B). 55 Similarly to the HIF prolyl-4-hydroxylases TaPHD and HsPHD2, PPHD and CrP4H both employ a HXD...H motif for metalbinding in the active site (H143 CrP4H , D145 CrP4H , H227 CrP4H , and H124 PPHD , D126 PPHD, H183 PPHD ). 38,66 In the PPHD crystal structure, the metal ion is octahedrally coordinated and one of the coordination sites is occupied by a water molecule (Figure 6), implying evolutionary conservation in the active site arrangement between the HIF prolyl-4-hydroxylases in animals and the bacterial non-HIF prolyl-4-hydroxylase PPHD. By contrast, in the CrP4H crystal structure, 66 the metal ion is tetrahedrally coordinated; however, in this structure, binding of a water molecule to the metal is apparently blocked by an acetate ion ( Figure 6E), and hence, at least in this detail, it may not be fully representative of the solution structure.

Discussion
It has been proposed that the ferrous iron dependent oxygenases may have evolved as a response to the advent of photosynthetically produced oxygen, possibly from precursors that used relatively bioavailable iron (comparable to eg, zinc) ions, predominantly in non-redox processes. 68,69 The widespread occurrence of 2OG oxygenases in prokaryotes, coupled with their close relationship with two TCA cycle intermediates (2OG and succinate), suggests they may have early origins. The discovery that 2OG oxygenases play key roles in hypoxia sensing in higher animals was therefore of interest from a broad evolutionary perspective. We are interested in exploring the biological distribution of hypoxia sensing mechanisms, as well as 2OG oxygenases and related enzymes. We have found that there is a functional HIFα-PHD-pVHL triad in the simplest animal, T. adhaerens. 16 Although bioinformatics and limited experimental analyses have not provided evidence for HIF transcription factors beyond animals, 16 PHD-like enzymes are much more widely distributed ( Figure S2). 16,37,38 The PHDs are related to the pro-collagen prolyl hydroxylases and are part of a structurally distinct 2OG oxygenase subfamily; the PH VIII subfamily 15 which includes other PHDs of biological interest. In the unicellular eukaryote Dictyostelium discoideum, a cytoplasmic PHD homologue, which is proposed to be involved in hypoxic responses,

69
Biophysical analyses of the oxygen sensing PHDs from the simplest animal T. adhaerens catalyzes hydroxylation of a proline-residue in Skp1, inducing further post-translational modification of Skp1 with a pentasaccharide group that is linked via the hydroxyproline residue alcohol. 37,70 Recent studies have shown that PHD homologues also exist in bacteria, where a PHD present in Pseudomonas species catalyzes hydroxylation of the ribosome associated protein EF-Tu. 38 Central to the hypoxia sensing ability of the PHDs is the capacity for their hydroxylase activity in many cell types to be limited by oxygen availability. 12,29 It is proposed that this is in substantial part due to the slow reaction of oxygen with the active site Fe(II) of the PHDs. 31 By contrast there is evidence that the human PHDs are (at least normally) relatively less sensitive to changes in Fe(II) or 2OG availability, consistent with the proposed specialized role for them as hypoxia sensors. 30 The question then arises as to what extent are these apparently special properties of the human PHDs (especially HsPHD2) conserved in the PHDs/PHD-like enzymes present in animals and beyond. The results presented here provide evidence that the apparently unusual kinetic properties of HsPHD2 are conserved in TaPHD ( Figure S7), likely as a consequence of the conservation of key structural features.
The overall folds of the catalytic domains of TaPHD and the HsPHD2 are very similar, including with respect to the elements involved in substrate binding, in particular the conformationally mobile β2/β3-finger-loop 33 and the C-terminal helix (Figures 5 and 6). Detailed comparison of the PHD-substrate interactions reveals that the active site region interactions between TaPHD and TaODD more closely resemble those between HsPHD2 and HsHIF1α CODD, rather than those of HsPHD2 and HsHIF1α NODD (Figures 4 and 6). Thus, our results support the proposal that the PHD2/CODD-type ODD couple evolved first and is likely the most important of the PHD/ODD interactions in higher animals where there are multiple PHDs and ODDs. 16 Given the likely relatively rigid nature of the Pro-Pro-Pro motif in TaODD, it also appears that TaODD binds to TaPHD in a less flexible manner than does HsHIF1α NODD to HsPHD2, and to a lesser extent, HsHIF1α CODD ( Figure  4). Although solution studies are required to validate this proposal, these differences may further reflect the evolved roles of the HsPHDs (and likely those in other higher animals) in accepting multiple substrates (NODDs/CODDs of HIFα isoforms), and maybe, non-HIF substrates [60][61][62][63][64][65] compared to TaPHD. The conservation in overall fold is also the case with respect to the prokaryotic PHD, PPHD, which acts on the highly abundant globular GTP-utilizing protein EF-Tu, 38 rather than the disordered ODD regions of HIFα ( Figure 5). (A role of PPHD in hypoxia sensing, if any, has yet to be identified).
Crucially, the structural elements which enable the unusually slow reaction of HsPHD2 with oxygen, 31,34 (as well as tight 2OG/Fe(II) binding and low substrate uncoupled turnover), appear to be conserved in TaPHD as observed crystallographically and manifested in the available kinetic studies ( Figures 2 and 6). Since humans and T. adhaerens are at near opposite ends of animal evolution, it seems possible that hypoxia sensing PHD homologues in organisms between them will have similar kinetic properties. Interestingly, TaPHD has a higher affinity for 2OG than HsPHD2 as observed by both K m comparisons and NMR binding studies. The difference in 2OG K m may reflect different 2OG metabolism in T. adhaerens and humans. Further, like HsPHD2, TaPHD forms an unusually stable TaPHD.Fe.2OG complex (half-life >24 h), consistent with both enzymes being specialized to preferentially respond to changes in oxygen, rather than 2OG/Fe(II), availability. 11 Our overall results suggest that the evolution of the HIFα (with a TaODD/CODD type ODD)-PHD2-pVHL triad, may have been an important step in evolution of animals. A key element of this triad is a PHD2type enzyme with appropriate kinetic properties, enabling it to "focus" on hypoxia sensing. The results also suggest that searching for appropriate kinetic and structural properties may be a general approach in helping to identify candidate (hypoxia) sensing enzymes.

Data availability
The TaPHD.Mn and TaPHD.TaODD structures are deposited in the RCSB PBD: entries 6EY1 and 6F0W, respectively.

Publish your work in this journal
Submit your manuscript here: https://www.dovepress.com/hypoxia-journal Hypoxia is an international, peer-reviewed, open access journal that aims to improve understanding of the biological response to hypoxia. The journal will publish original research articles, reviews, methodological advances, clinical studies, and expert opinions that identify developments in the regulation of the physiological and pathological responses to hypoxia and in the therapeutic targeting of hypoxia-responsive pathways. The manuscript management system is completely online and includes a very quick and fair peer-review system, which is all easy to use. Visit http://www.dovepress.com/testimonials.php to read real quotes from published authors.