Final Report: Evaluation of GO:0004721 (Phosphoprotein Phosphatase Activity) for Pseudomonas putida AceK (Q88EA1)

Executive Judgment

Verdict: Partially supported / too general — Failure mode #1 (granularity / family-vs-subfamily).

The seed hypothesis that AceK (Q88EA1) from Pseudomonas putida KT2440 possesses phosphoprotein phosphatase activity (GO:0004721) is correct in substance but insufficiently specific. AceK is a genuine bifunctional isocitrate dehydrogenase kinase/phosphatase whose phosphatase catalytic residues are fully conserved relative to the experimentally characterized E. coli ortholog (P11071). The single most decisive piece of evidence is the complete conservation of the DXDX(T/V) phosphatase catalytic motif (FYDYDEI, 7/7 positions identical) including the catalytic Asp477→Asp480, which has been shown by mutagenesis and structural studies to be indispensable for the dephosphorylation reaction (PMID: 25272278). However, GO:0004721 is a broad parent term; a more specific child term — GO:0101014 ([isocitrate dehydrogenase (NADP+)] phosphatase activity) — precisely describes AceK's substrate-specific phosphatase function and should replace the propagated annotation, matching the specificity already applied to the kinase half of the enzyme (GO:0008772).

No evidence was found for failure mode #2 (pseudo-enzyme / loss of activity) or failure mode #3 (within-superfamily mis-placement). All catalytic, metal-binding, and ATP-binding residues are intact, and the AlphaFold structural model shows high confidence (pLDDT 89–96) at every active-site position.

Summary

This investigation assessed whether the automated TreeGrafter/PANTHER annotation of phosphoprotein phosphatase activity (GO:0004721) for Pseudomonas putida KT2440 AceK (UniProt Q88EA1) is warranted. AceK is the sole enzyme responsible for reversible phosphorylation of isocitrate dehydrogenase (IDH) in Gram-negative bacteria, toggling the glyoxylate shunt on and off in response to carbon source availability. Its bifunctional kinase/phosphatase mechanism was first characterized in E. coli and represents a landmark discovery in prokaryotic signaling.

Through pairwise sequence alignment against the crystallographically characterized E. coli AceK (P11071; 45.5% identity over 571/578 residues), we confirmed that every residue known to be essential for phosphatase catalysis is conserved in the P. putida ortholog: the catalytic Asp477→Asp480, the APE-motif Glu439→Glu442, the DXDX(T/V) phosphatase motif (100% identity, 7/7 positions), and the ATP-binding lysine K336→K339 required for the unusual ATP/ADP-dependent phosphatase mechanism. AlphaFold structural prediction (AF-Q88EA1-F1, v6) corroborates a well-folded kinase/phosphatase with mean pLDDT of 92.0 and uniformly high confidence at catalytic positions. Domain architecture analysis confirms both proteins share the identical Pfam domains (PF06315 + PF20423), PANTHER subfamily (PTHR39559:SF1), and HAMAP family (MF_00747).

The sole issue is annotation granularity. GO:0004721 is the family-level phosphatase term. A dedicated child term, GO:0101014 ([isocitrate dehydrogenase (NADP+)] phosphatase activity), captures AceK's precise substrate specificity but currently carries zero annotations in QuickGO — a striking asymmetry compared to the kinase counterpart GO:0008772, which has over 2,600 annotations. Additionally, GO:0004722 (protein serine/threonine phosphatase activity) is applicable given that AceK dephosphorylates a specific serine residue on IDH. The recommended curation action is make-more-specific: annotate with GO:0101014 and optionally GO:0004722, replacing the overly broad GO:0004721.

Key Findings

Finding 1: AceK Is a Genuine Bifunctional Kinase/Phosphatase with Conserved Catalytic Residues

Pairwise alignment (EMBOSS Needle, BLOSUM62 matrix, gap penalty 10.0/0.5) of E. coli AceK (P11071, 578 aa) against P. putida AceK (Q88EA1, 571 aa) yielded 45.5% identity (270/594 aligned positions) and 59.8% similarity (355/594), with an alignment score of 1241.5. Both proteins share the same PANTHER subfamily (PTHR39559:SF1), Pfam domain architecture (PF06315 + PF20423), HAMAP family (MF_00747), InterPro family (IPR010452), and dual EC classifications (EC 2.7.11.5 for kinase; EC 3.1.3.- for phosphatase).

Critical catalytic residues mapped from the E. coli crystal structures and mutagenesis data are fully conserved at all positions essential for activity:

Crucially, the two residues that are not conserved (N377→D380 and D403→L406) were experimentally shown to be non-essential for maximal catalytic activity — their mutation in E. coli only affected Km for magnesium, not Vmax (PMID: 11258918). All four essential residue groups (D371, K336, D477, and the DXDX(T/V) motif) are perfectly conserved.

The DXDX(T/V) phosphatase motif, a hallmark of this phosphatase family, spans positions 476–482 in E. coli (FYDYDEI) and is 100% conserved at all seven positions in P. putida:

E. coli:   F473  Y474  D475  Y476  D477  E478  I479
P. putida: F476  Y477  D478  Y479  D480  E481  I482
Match:       =     =     =     =     =     =     =
   7/7 positions CONSERVED (100%)

This conservation provides strong computational evidence that the phosphatase catalytic mechanism — a stepwise addition–elimination reaction using Asp477/Asp480 and a single Mg²⁺ ion (PMID: 25272278) — is retained in P. putida AceK.

Finding 2: GO:0004721 Is Correct but Too General — GO:0101014 Is the Appropriate Specific Term

Systematic traversal of the GO hierarchy via the QuickGO API revealed that GO:0101014 ([isocitrate dehydrogenase (NADP+)] phosphatase activity) is a direct is_a child of GO:0004721 (phosphoprotein phosphatase activity). Its formal definition — "Catalysis of the reaction: [isocitrate dehydrogenase] phosphate + H₂O = [isocitrate dehydrogenase] + phosphate" — precisely describes what AceK does.

A striking annotation asymmetry exists in the GO databases: for AceK's kinase function, TreeGrafter correctly propagates the specific term GO:0008772 ([isocitrate dehydrogenase (NADP+)] kinase activity, 2,635 annotations in QuickGO), but for the phosphatase function it propagates only the broad parent GO:0004721. GO:0101014 currently has zero annotations in QuickGO, suggesting it is systematically underutilized in automated pipelines.

Finding 3: AlphaFold Model Confirms a Well-Folded Active Site with High Confidence

The AlphaFold model (AF-Q88EA1-F1, version 6) for P. putida AceK shows a mean pLDDT of 92.0, with 80.6% of residues scoring above 90 (very high confidence) and 16.6% between 70–90 (confident). At the catalytic residues critical for phosphatase function, pLDDT values are uniformly high:

No low-confidence regions were detected near the active site, ruling out the possibility that structural disorder or a degenerate fold might compromise catalytic function. This provides independent structural evidence against failure mode #2 (pseudo-enzyme).

{{figure:aceK_validation_summary.png|caption=Summary of AceK validation: AlphaFold pLDDT confidence at catalytic residues, phosphatase motif conservation, GO term hierarchy, and overall verdict supporting GO:0101014 as the most appropriate specific annotation.}}

Finding 4: GO:0101014 and GO:0004722 Are Complementary Sibling Terms

Detailed analysis of the GO hierarchy confirmed that GO:0101014 and GO:0004722 are both direct children of GO:0004721 but are not related to each other by any is_a or part_of relationship. They represent independent axes of functional specificity:

• GO:0101014 specifies the substrate identity (isocitrate dehydrogenase)
• GO:0004722 specifies the target residue type (Ser/Thr)

Both are more informative than the seed term GO:0004721 and both are applicable to AceK. The E. coli ortholog P11071 already carries GO:0004722 with IDA evidence. Dual annotation with both terms would provide the most complete and accurate functional description.

Finding 5: ATP-Binding Site Is Conserved — Essential for ATP-Dependent Phosphatase Mechanism

AceK's phosphatase activity is mechanistically unusual: it is strictly dependent on ATP or ADP hydrolysis. This was demonstrated directly by Miller et al. (1996), who showed that "the IDH phosphatase reaction required either ATP or ADP, [but] was not supported by the nonhydrolyzable ATP analogue 5'-adenylyl imidodiphosphate" (PMID: 8702587). This makes ATP-binding site integrity a prerequisite for phosphatase function.

Alignment of the ATP-binding region showed 7/8 positions conserved (88%) between E. coli and P. putida AceK:

The single substitution (I318→V321) is conservative (both hydrophobic aliphatic). The critical lysine K336, whose mutation inhibits enzyme activity, is perfectly conserved as K339. The ATP-dependent phosphatase mechanism is therefore structurally supported.

Independent Family/Function Assignment

Based on this analysis, the most likely specific molecular function of P. putida AceK (Q88EA1) is:

• Primary phosphatase function: [isocitrate dehydrogenase (NADP+)] phosphatase activity (GO:0101014) — a more specific child of the seed term GO:0004721
• Complementary phosphatase function: protein serine/threonine phosphatase activity (GO:0004722) — a sibling of the seed term under GO:0004721
• Kinase function: [isocitrate dehydrogenase (NADP+)] kinase activity (GO:0008772) — the kinase half of this bifunctional enzyme

The nearest characterized homolog is E. coli AceK (P11071), which has been extensively studied by X-ray crystallography (PMID: 20505668), site-directed mutagenesis (PMID: 11258918; PMID: 25272278), enzymology (PMID: 8702587), and metal-binding characterization (PMID: 31235769).

Mechanistic Model / Interpretation

AceK occupies a unique position in prokaryotic signaling: it is a single polypeptide that catalyzes both the phosphorylation (inactivation) and dephosphorylation (reactivation) of isocitrate dehydrogenase (IDH), the gatekeeper enzyme that partitions isocitrate between the TCA cycle and the glyoxylate bypass. The mechanistic model, established from extensive E. coli studies, is:

     AceK (bifunctional)
    ┌─────────────────────────┐
    │                         │
  IDH-active ──────┤  KINASE  (ATP→ADP)      ├──────► IDH-P (inactive)
    │  D371/374, K336/339     │
    │                         │
  IDH-P (inactive)─┤  PHOSPHATASE            ├──────► IDH-active + Pi
    │  D477/480, DXDX(T/V)   │
    │  (ATP/ADP-dependent)    │
    └─────────────────────────┘
              │
    Allosteric regulation:
    AMP → activates phosphatase / inhibits kinase
    Isocitrate → inhibits kinase
    NADPH → inhibits both activities

Both activities share a single active site but rely on distinct catalytic residues: the kinase uses D371 and the canonical protein-kinase catalytic triad, while the phosphatase uses D477 within the DXDX(T/V) motif characteristic of a distinct phosphotransferase family. The phosphatase reaction proceeds via a stepwise addition–elimination mechanism requiring Mg²⁺ and hydrolyzable ATP/ADP (PMID: 25272278; PMID: 22889914). The AMP-binding allosteric site acts as a conformational switch: AMP binding exposes ATP and favors phosphatase activity over kinase activity (PMID: 20505668).

In P. putida KT2440, AceK plays the same regulatory role. The glyoxylate shunt is essential when the bacterium grows on acetate or fatty acids as sole carbon sources, and omics studies have confirmed that isocitrate lyase — the first enzyme of the bypass — is among the most abundant proteins during growth on alternative carbon sources such as butanol (PMID: 26986205). Additionally, studies on ICDH regulation in Pseudomonas fluorescens have demonstrated that ICDH phosphorylation state is dynamically modulated in response to metabolic stress (PMID: 17573472; PMID: 31828449). AceK's regulatory function is therefore physiologically relevant in Pseudomonads.

The conservation of the aceK gene in P. putida, the shared PANTHER subfamily assignment (PTHR39559:SF1), identical Pfam domains, and ≥45% sequence identity with 100% conservation of all known essential catalytic residues collectively indicate that P. putida AceK retains both kinase and phosphatase activities. Structural studies on E. coli AceK have further shown that "the highly stringent AceK binding sites on ICDH are maintained only in Gram-negative bacteria" (PMID: 21870819), and P. putida is a Gram-negative organism — confirming that IDH substrate recognition is expected to be intact.

Genomic context note: In P. putida KT2440, aceK (PP_4565, complement 5184744..5186459) is not in an aceBAK operon as in E. coli. The glyoxylate shunt genes aceA (PP_4116, isocitrate lyase) and aceB (malate synthase) are located elsewhere on the chromosome. This dispersed gene organization is typical for Pseudomonas species and does not affect the functional assignment.

Active-Site / Placement Analysis

Comprehensive Residue Conservation Summary

All 4 essential residue groups are conserved. The 2 non-essential positions that changed (N377, D403) were experimentally shown not to affect maximal activity in E. coli (PMID: 11258918).

Local Alignment Context at Key Regions

Phosphatase motif region:
  E. coli:   GRVVFYDYDEICYMTEV
      |||||||||||...|||
  P. putida: GRVVFYDYDEISFLTEV

Active site region:
  E. coli:   HDRVGRMADTQEFENFV
      .||||||||||||..|.
  P. putida: VDRVGRMADTQEFADFR

Domain Architecture Match

AlphaFold Structural Confidence

• Model: AF-Q88EA1-F1 (AlphaFold v6)
• Global quality: Mean pLDDT = 92.0; 80.6% residues >90 (very high), 16.6% 70–90 (confident)
• Catalytic residue pLDDT: D374 = 91.6, D480 = 89.1, E442 = 95.7, FYDYDEI motif = 89.1–96.3

The AlphaFold model independently confirms the catalytic site is in a well-ordered, confidently predicted conformation consistent with an active enzyme.

Conclusion: All catalytic residues essential for both kinase and phosphatase activity are intact. The ATP-binding site required for ATP-dependent phosphatase activity is highly conserved (88%). The protein is correctly placed in PTHR39559:SF1 with no evidence of pseudo-enzyme status (failure mode #2 excluded) or within-superfamily mis-placement (failure mode #3 excluded).

Evidence Matrix

GO Curation Implications

Recommended curation action: make-more-specific

Key curation observations:

1. Annotation asymmetry: The kinase function receives the specific term GO:0008772 ([isocitrate dehydrogenase (NADP+)] kinase activity) via HAMAP rule (GO_REF:0000120), but the phosphatase function receives only the general GO:0004721 via TreeGrafter (GO_REF:0000118). Both should be at the same level of specificity.

2. GO:0101014 is systematically underutilized: This term exists in the ontology with zero annotations in QuickGO, yet it precisely describes the function of all AceK enzymes. There are 2,635 annotations for the kinase counterpart GO:0008772. This is a clear gap in the automated annotation pipelines.

3. GO:0101014 and GO:0004722 are orthogonal siblings: Both are direct is_a children of GO:0004721 but neither is an ancestor of the other. GO:0101014 specifies the substrate (IDH), while GO:0004722 specifies the residue type (Ser/Thr). Both should be annotated for maximum specificity.

4. Evidence code upgrade: The evidence could be upgraded from IEA to ISS (inferred from sequence or structural similarity) given the 45.5% identity, 100% catalytic residue conservation, and extensive experimental characterization of the E. coli ortholog.

5. Retain, do not remove: The seed term is semantically correct and should not be removed — it should be refined to one or both more specific child terms.

Conflicts, Knowledge Gaps, and Discriminating Tests

Conflicts

No significant conflicts were identified. All evidence — sequence conservation, domain architecture, structural prediction, phylogenetic classification, and physiological context — converges on the conclusion that Q88EA1 is a functional bifunctional AceK with both kinase and phosphatase activities.

Paralog confusion is unlikely: AceK is encoded by a single gene (aceK) in P. putida KT2440, and no paralogous isocitrate dehydrogenase kinase/phosphatases exist in this organism.

Knowledge Gaps

1. No direct experimental data for P. putida AceK exists. All functional inferences are based on homology to the E. coli enzyme. While the conservation is compelling, kinetic parameters (Km, kcat for phosphatase activity) have not been measured for the P. putida enzyme.

2. Allosteric regulation may differ. The AMP-binding allosteric site that switches between kinase and phosphatase modes in E. coli AceK has not been evaluated for conservation in P. putida. Metabolic effector sensitivities could differ between species.

3. ATP/ADP-dependence as an annotation issue. AceK's phosphatase activity is mechanistically unusual — strictly ATP/ADP-dependent. The GO:0004721 definition ("a phosphoprotein + H₂O = a protein + phosphate") does not capture this cofactor requirement. Whether this matters for annotation purposes is a curation policy question.

4. GO:0101014 has zero annotations. Even the experimentally characterized E. coli AceK (P11071) lacks this annotation, suggesting a systematic gap in GO curation for this enzyme family's phosphatase function.

Most Efficient Discriminating Tests

Limitations

1. All functional evidence is inferred by homology. No direct experimental characterization of P. putida AceK phosphatase activity has been published. The conclusion rests on sequence conservation (45.5% identity) and catalytic residue identity relative to the well-characterized E. coli ortholog.

2. Pairwise alignment only. Active-site analysis was performed with pairwise alignment (EMBOSS Needle) rather than a full multiple sequence alignment across the AceK family. A broader MSA would provide stronger phylogenetic signal.

3. AlphaFold pLDDT scores reflect prediction confidence, not experimental validation. High pLDDT at active-site residues supports structural plausibility but does not confirm enzymatic activity.

4. GO hierarchy analysis depends on current GO structure. GO:0101014 may be subject to revision; its current status (zero annotations) should be verified with GO curators before formal adoption.

5. Some API queries were not fully successful due to endpoint limitations during the investigation; domain architecture comparisons relied partly on UniProt cross-references rather than direct InterPro API queries.

Proposed Follow-up Actions

1. Immediate (curator action): Replace GO:0004721 with GO:0101014 for Q88EA1, and consider adding GO:0004722. Upgrade evidence code from IEA to ISS given the strength of the sequence conservation evidence.

2. Short-term (computational): Perform a systematic scan of all AceK orthologs currently annotated with GO:0004721 to identify candidates for upgrade to GO:0101014. The zero-annotation status of this term represents a pipeline-wide gap.

3. Medium-term (experimental): Express recombinant P. putida AceK and confirm phosphatase activity on purified P. putida IDH in vitro, enabling IDA evidence assignment.

4. Long-term (structural): Solve the crystal structure of P. putida AceK to confirm active-site geometry and identify any species-specific allosteric differences that might affect regulation.

Artifacts

• OpenScientist final report
• OpenScientist final report
• OpenScientist aceK validation summary
  
• OpenScientist plot 1