AIGR TreeGrafter Function-Inference Stress Test

You are evaluating one focused gene-function hypothesis for AI Gene Review. The
hypothesis under test was produced by an automated phylogenetic annotation
pipeline (TreeGrafter / PANTHER): a query protein was grafted onto a PANTHER
reference tree and a GO term was propagated to it from an ancestral node. Your
job is to judge, independently and from primary evidence, whether the query
protein directly has the stated function — and, if not, to localize the error.

This is not a general gene overview. Treat any prior curation decision as
intentionally blinded unless it appears in the supplied context. Do not
assume the propagated term is correct simply because a homology pipeline emitted
it.

Target Gene

• Organism code: ECOLX
• Taxon: Escherichia coli (NCBITaxon:562)
• Gene directory: mcr-1
• Gene symbol: mcr-1
• UniProt accession: A0A0R6L508

Focus

• Focus type: function_assignment
• Hypothesis slug: function-hypothesis-go-0016776
• Source file: genes/ECOLX/mcr-1/mcr-1-ai-review.yaml
• Source selector: existing_annotations[5].function_hypothesis

Seed Hypothesis (propagated by TreeGrafter/PANTHER)

mcr-1 has phosphotransferase activity, phosphate group as acceptor (GO:0016776).

Term and Decision Context

• Term: phosphotransferase activity, phosphate group as acceptor (GO:0016776)
• Evidence type: IEA
• Original reference: GO_REF:0000118

Reference Context

• GO_REF:0000118
• PMID:29079699

Source Context YAML

term:
  id: GO:0016776
  label: phosphotransferase activity, phosphate group as acceptor
evidence_type: IEA
original_reference_id: GO_REF:0000118

Research Objective

Decide whether mcr-1 directly has the stated function. Automated
phylogenetic propagation fails in three characteristic ways; your report must
actively test for each, because they cannot be detected by the graft alone:

1. Granularity / family-vs-subfamily. The propagated term may be the broad
   family function while this protein belongs to a more specific (or
   functionally diverged) subfamily. Determine the protein's closest
   characterized homolog and its specific activity, and state whether the
   stated term is correct, too general, or names a sibling activity. (Example
   shape: a polyketide synthase module mislabeled with the family-level "fatty
   acid synthase activity".)
2. Pseudo-enzyme / loss of activity. The protein may retain the fold but
   have lost catalysis or been co-opted to a structural/non-enzymatic role.
   Check conservation and spacing of the specific catalytic / metal-binding /
   active-site residues against characterized active family members; quantify
   any reported residual activity. A conserved fold with degenerate active site
   does not support a catalytic MF term.
3. Within-superfamily mis-placement. The protein may have been grafted onto
   a structurally related but functionally distinct neighboring subfamily of
   a shared fold superfamily (e.g. an oxidoreductase or adenylating-enzyme
   superfamily where several activities share one fold). Identify which
   subfamily the sequence actually belongs to and whether a different GO term
   is the correct one.

Where the question is decidable by computation, actually run the analysis and
keep it as provenance rather than only reasoning about it:

• Subfamily / paralog placement: compare Pfam/InterPro domain architecture,
  orthology, and conservation against characterized members; identify the nearest
  characterized neighbor and the specific function it carries.
• Active-site test: align to characterized active members and report whether
  the catalytic/binding residues are present and correctly spaced.
• Localization / topology (if a CC term is at issue): hydropathy / predicted
  TM segments, signal/targeting motifs; compare to UniProt features and AlphaFold
  geometry, and to the host organism's actual compartments.

Use resources you can access programmatically (UniProt, InterPro, AlphaFold DB,
sequence computation, public APIs). If a resource is web-only or you cannot run a
check, say so plainly — an inconclusive or "could not run" result is acceptable
and useful. Never fabricate a result. Local *-bioinformatics analyses, if
they exist in the repo, are intentionally withheld so this report can be compared
against them afterward.

Required Output

Executive Judgment

Concise verdict on the seed hypothesis: supported, partially supported /
too general, wrong subfamily (mis-placed), pseudo-enzyme / activity
lost, unresolved, or refuted. Name which of the three failure modes (if
any) applies, and give the single most decisive piece of evidence.

Independent Family/Function Assignment

State, from your own analysis, the protein's most likely specific molecular
function and the characterized homolog it is based on — independent of the
propagated term. Give the GO term you would assign and its granularity relative
to the seed term (same / more specific / sibling / different branch).

Evidence Matrix

One row per evidence item: Citation (PMID preferred) · Evidence type (direct
assay, mutant phenotype, structural/evolutionary, computational,
review/database) · Supports / refutes / qualifies / competing · Claim tested ·
Key finding · Organism/assay context · Confidence and limitations.

Active-Site / Placement Analysis

Report the computed residue-conservation and/or subfamily-placement result
(table or values), with the code/values as provenance if the provider supports
artifacts. State explicitly whether catalytic residues are intact.

GO Curation Implications

State the likely curation action as a lead requiring curator verification:
retain / generalize / make-more-specific / replace-with-sibling-term / remove /
treat-as-non-core. Avoid "protein binding" as a final recommendation.

Conflicts, Knowledge Gaps, Discriminating Tests

Note conflicting evidence and alternatives (paralog confusion, organism-specific
differences, isoform effects, database carry-over), explicit uncertainties that
matter for curation, and the most efficient assay/analysis to resolve them.

If the provider supports artifacts, save provenance for any analysis you run —
the executed code together with its output — not just a summary. Do not
synthesize a figure implying an analysis you did not run.

AIGR TreeGrafter Function-Inference Stress Test Report: MCR-1 / GO:0016776

Executive Judgment

Verdict: Partially supported / too general (Failure mode 1: Granularity / family-vs-subfamily)

The seed hypothesis that MCR-1 has "phosphotransferase activity, phosphate group as acceptor" (GO:0016776) is not wrong but is insufficiently specific and represents a classic granularity error in phylogenetic propagation. MCR-1 is a well-characterized, catalytically active lipid A phosphoethanolamine (pEtN) transferase (EC 2.7.4.30) that transfers a phosphoethanolamine moiety — not a bare phosphate group — from phosphatidylethanolamine (PE) to the 4ʹ-phosphate of lipid A (zinkle2025mechanisticbasisof pages 7-8, zinkle2025mechanisticbasisof pages 1-2). Because the acceptor in this reaction is indeed a phosphate group on lipid A, GO:0016776 is technically a valid ancestor term in the GO hierarchy. However, the propagated term captures only the broad enzyme commission class (EC 2.7.4.x) and obscures the protein's actual substrate specificity, its phosphoethanolamine-transfer chemistry, and its biological role in colistin resistance. The most decisive piece of evidence is the deep mutational scanning study by Sun & Palzkill (2021) combined with the full-length cryo-EM structure by Zinkle et al. (2025), which together demonstrate that MCR-1 uses Thr285-mediated nucleophilic attack on phosphatidylethanolamine to form a covalent pEtN-enzyme intermediate, followed by transfer of pEtN to the 4ʹ-phosphate of lipid A — a reaction that is more specific than "phosphotransferase activity, phosphate group as acceptor" (sun2021deepmutationalscanning pages 10-11, sun2021deepmutationalscanning pages 3-5, zinkle2025mechanisticbasisof pages 7-8).

No evidence for pseudo-enzyme status (failure mode 2) or within-superfamily mis-placement (failure mode 3) was found. The enzyme is catalytically active with all key residues intact.

Independent Family/Function Assignment

Protein identity: MCR-1 (UniProt A0A0R6L508) is a plasmid-encoded lipid A phosphoethanolamine transferase belonging to the alkaline phosphatase superfamily, classified as a Class I phosphoethanolamine transferase (PET) (anandan2020structureandfunction pages 14-17, samantha2020lipidaphosphoethanolamine pages 8-12, anandan2017structureofa pages 1-2).

Closest characterized homolog: NmEptA (EptA from Neisseria meningitidis), sharing ~36% overall sequence identity with MCR-1, with higher conservation (~50%) in the C-terminal catalytic domain (samantha2020lipidaphosphoethanolamine pages 5-8, sun2021deepmutationalscanning pages 11-13). The full-length crystal structure of NmEptA (Anandan et al. 2017, PNAS) serves as the canonical structural template for this enzyme family (anandan2017structureofa pages 1-2).

Specific molecular function: MCR-1 catalyzes the Zn²⁺-dependent transfer of phosphoethanolamine from phosphatidylethanolamine to the 4ʹ-phosphate group of lipid A, proceeding through a ping-pong mechanism with a covalent pEtN-Thr285 intermediate (zinkle2025mechanisticbasisof pages 7-8, sun2021deepmutationalscanning pages 10-11, thai2023phosphoethanolaminetransferasesas pages 16-18, sun2021deepmutationalscanning pages 3-5).

Recommended GO MF term: The most appropriate term would be a specific GO term for "lipid A phosphoethanolamine transferase activity" or the broader but still more precise "phosphoethanolamine transferase activity". If such a specific GO term does not yet exist, one should be requested. The existing GO:0016776 is a valid ancestor/parent term but should not be the primary annotation.

Granularity relative to seed term: The correct annotation is more specific than GO:0016776 — it is a child/descendant function, not a sibling or different-branch activity.

Evidence Matrix

The following table summarizes the primary evidence evaluated for this assessment:

Table: This table summarizes the key primary and review evidence used to judge whether MCR-1 should carry the broad GO:0016776 annotation. It highlights that MCR-1 is an active lipid A phosphoethanolamine transferase, so the propagated term is best treated as broadly true but too general.

Active-Site / Placement Analysis

MCR-1 retains all catalytic and metal-binding residues characteristic of active Class I lipid A phosphoethanolamine transferases. Deep mutational scanning of 23 active-site positions revealed that 17 positions (75%) are essential for function, with wild-type residues strongly predominating after selection for polymyxin resistance (sun2021deepmutationalscanning pages 10-11, sun2021deepmutationalscanning pages 3-5, sun2021deepmutationalscanning pages 11-13). The core catalytic residues — Thr285 (nucleophile), Glu246, Asp465, His466 (Zn²⁺ coordination) — are identical to their equivalents in NmEptA (Thr280, Glu240, Asp452, His453) (samantha2020lipidaphosphoethanolamine pages 8-12, thai2023phosphoethanolaminetransferasesas pages 14-16, anandan2020structureandfunction pages 14-17). Mutation of Thr285 to alanine completely abolishes catalytic activity while stabilizing the protein, confirming its role as the catalytic nucleophile rather than a structural residue (sun2021deepmutationalscanning pages 11-13, sun2021deepmutationalscanning pages 10-11). Conservative substitutions at metal-binding positions (E246D, D465E) significantly reduce both pEtN transferase activity and polymyxin resistance, demonstrating the stringency of the active-site requirements (sun2021deepmutationalscanning pages 10-11).

The detailed residue-by-residue analysis and subfamily placement comparison is presented below:

Table: This table summarizes residue-level catalytic evidence and subfamily placement for MCR-1. It shows that MCR-1 is an intact, active lipid A phosphoethanolamine transferase in the alkaline phosphatase superfamily and that GO:0016776 is best treated as a broad, too-general annotation.

Key conclusions from active-site analysis:
1. MCR-1 is not a pseudo-enzyme. All catalytic residues are intact and functionally validated by mutagenesis.
2. MCR-1 is correctly placed within the Class I lipid A phosphoethanolamine transferase subfamily, not a neighboring subfamily of the alkaline phosphatase superfamily (which includes enzymes like BcsG that modify cellulose, or EptB/CptA that modify different LPS components) (schumann2024themultifacetedroles pages 4-6, anderson2020theescherichiacoli pages 16-18).
3. MCR-1 shows subfamily-specific site selectivity: it preferentially modifies the 4ʹ-phosphate of lipid A, in contrast to EptA and MCR-9 which prefer the 1-phosphate (schumann2024themultifacetedroles pages 4-6). This further supports a more specific annotation than the broad GO:0016776.

GO Curation Implications

Recommended curation action: Make-more-specific

The current GO:0016776 annotation should be replaced with a more specific term that accurately reflects lipid A phosphoethanolamine transferase activity. The reasoning is as follows:

1. GO:0016776 is not incorrect — it is a valid parent/ancestor term because MCR-1 does transfer a phosphorus-containing group (pEtN) to a phosphate-group acceptor (lipid A phosphate). The enzyme's EC classification (2.7.4.30) falls under the EC 2.7.4.x subclass ("phosphotransferases with a phosphate group as acceptor"), which corresponds to GO:0016776 (anandan2017structureofa pages 1-2).

2. However, GO:0016776 is too general. It does not distinguish MCR-1's pEtN transferase chemistry from unrelated phosphotransferases (e.g., nucleotide kinases, polyphosphate kinases) that also have a phosphate-group acceptor. The term fails to capture the biologically and biochemically defining feature of MCR-1: the transfer of a phosphoethanolamine moiety from a glycerophospholipid donor to lipid A (zinkle2025mechanisticbasisof pages 7-8, thai2023phosphoethanolaminetransferasesas pages 16-18).

3. The ideal annotation would be a term such as "lipid A phosphoethanolamine transferase activity" (if available in GO) or the more general "phosphoethanolamine transferase activity." If no such term exists at the appropriate granularity, a term request to the GO Consortium would be appropriate.

4. This is a characteristic TreeGrafter failure mode 1 (granularity). The PANTHER tree likely captured the broad alkaline-phosphatase-superfamily-level activity at a family node, then propagated it to MCR-1 without resolution to the specific pEtN transferase subfamily.

Conflicts, Knowledge Gaps, and Discriminating Tests

Conflicts and alternatives

• No conflict with enzymatic activity. MCR-1 is unambiguously an active enzyme, ruling out pseudo-enzyme concerns.
• No within-superfamily mis-placement. MCR-1 clearly belongs to the lipid A-modifying PET clade, not to the cellulose-modifying (BcsG), Kdo-modifying (EptB), or core-sugar-modifying (CptA/EptC) PET clades (schumann2024themultifacetedroles pages 4-6, anderson2020theescherichiacoli pages 16-18).
• Paralog confusion risk: The alkaline phosphatase superfamily contains diverse enzymes with the same fold but different activities (phosphatases, sulfatases, phosphodiesterases). TreeGrafter propagation at the superfamily level could erroneously assign any of these activities. In this case, the propagated term (GO:0016776) happens to be a valid ancestor but is too shallow in the hierarchy.

Knowledge gaps

• GO term availability: It is not confirmed from the retrieved evidence whether a specific GO term for "lipid A phosphoethanolamine transferase activity" currently exists in the Gene Ontology. This should be verified against the live GO database.
• Acceptor-site selectivity annotation: MCR-1's preference for the 4ʹ-phosphate over the 1-phosphate of lipid A (schumann2024themultifacetedroles pages 4-6) represents a further level of functional specificity not captured by any standard GO MF term.
• Quantitative kinetic parameters: While the enzyme is clearly active, precise kinetic constants (kcat, KM) for MCR-1 under standardized conditions are not well-established in the retrieved literature, which limits quantitative comparison with other PETs.

Discriminating tests

1. GO term hierarchy check: Query the GO database for child terms of GO:0016776 that specifically describe phosphoethanolamine transferase activity. If such a term exists, it should replace GO:0016776 for MCR-1.
2. InterPro/Pfam annotation review: Verify that MCR-1's InterPro entry maps to the specific PEtN transferase family (e.g., IPR017849/Pfam sulfatase-like domain) rather than a generic alkaline phosphatase superfamily entry, and use this to guide appropriate GO mapping.
3. In vitro kinetic assay with defined substrates: A direct enzymatic assay using purified full-length MCR-1, PE as donor, and Kdo₂-lipid A as acceptor, with product detection by mass spectrometry, would provide the definitive biochemical evidence for curator-verified GO annotation at the most specific level.

References

1. (zinkle2025mechanisticbasisof pages 7-8): Allen P. Zinkle, Mariana Bunoro Batista, Carmen M. Herrera, Satchal K. Erramilli, Brian Kloss, Khuram U. Ashraf, Kamil Nosol, Guozhi Zhang, Rosemary J. Cater, Michael T. Marty, Anthony A. Kossiakoff, M. Stephen Trent, Rie Nygaard, Phillip J. Stansfeld, and Filippo Mancia. Mechanistic basis of antimicrobial resistance mediated by the phosphoethanolamine transferase mcr-1. Nature Communications, Nov 2025. URL: https://doi.org/10.1038/s41467-025-65515-3, doi:10.1038/s41467-025-65515-3. This article has 4 citations and is from a highest quality peer-reviewed journal.

2. (zinkle2025mechanisticbasisof pages 1-2): Allen P. Zinkle, Mariana Bunoro Batista, Carmen M. Herrera, Satchal K. Erramilli, Brian Kloss, Khuram U. Ashraf, Kamil Nosol, Guozhi Zhang, Rosemary J. Cater, Michael T. Marty, Anthony A. Kossiakoff, M. Stephen Trent, Rie Nygaard, Phillip J. Stansfeld, and Filippo Mancia. Mechanistic basis of antimicrobial resistance mediated by the phosphoethanolamine transferase mcr-1. Nature Communications, Nov 2025. URL: https://doi.org/10.1038/s41467-025-65515-3, doi:10.1038/s41467-025-65515-3. This article has 4 citations and is from a highest quality peer-reviewed journal.

3. (sun2021deepmutationalscanning pages 10-11): Zhizeng Sun and Timothy Palzkill. Deep mutational scanning reveals the active-site sequence requirements for the colistin antibiotic resistance enzyme mcr-1. Dec 2021. URL: https://doi.org/10.1128/mbio.02776-21, doi:10.1128/mbio.02776-21. This article has 22 citations and is from a domain leading peer-reviewed journal.

4. (sun2021deepmutationalscanning pages 3-5): Zhizeng Sun and Timothy Palzkill. Deep mutational scanning reveals the active-site sequence requirements for the colistin antibiotic resistance enzyme mcr-1. Dec 2021. URL: https://doi.org/10.1128/mbio.02776-21, doi:10.1128/mbio.02776-21. This article has 22 citations and is from a domain leading peer-reviewed journal.

5. (anandan2020structureandfunction pages 14-17): Anandhi Anandan and Alice Vrielink. Structure and function of lipid a–modifying enzymes. Annals of the New York Academy of Sciences, 1459:19-37, Jan 2020. URL: https://doi.org/10.1111/nyas.14244, doi:10.1111/nyas.14244. This article has 54 citations and is from a peer-reviewed journal.

6. (samantha2020lipidaphosphoethanolamine pages 8-12): Ariela Samantha and Alice Vrielink. Lipid a phosphoethanolamine transferase: regulation, structure and immune response. Journal of Molecular Biology, 432:5184-5196, Aug 2020. URL: https://doi.org/10.1016/j.jmb.2020.04.022, doi:10.1016/j.jmb.2020.04.022. This article has 87 citations and is from a domain leading peer-reviewed journal.

7. (anandan2017structureofa pages 1-2): Anandhi Anandan, Genevieve L. Evans, Karmen Condic-Jurkic, Megan L. O’Mara, Constance M. John, Nancy J. Phillips, Gary A. Jarvis, Siobhan S. Wills, Keith A. Stubbs, Isabel Moraes, Charlene M. Kahler, and Alice Vrielink. Structure of a lipid a phosphoethanolamine transferase suggests how conformational changes govern substrate binding. Proceedings of the National Academy of Sciences, 114:2218-2223, Feb 2017. URL: https://doi.org/10.1073/pnas.1612927114, doi:10.1073/pnas.1612927114. This article has 153 citations and is from a highest quality peer-reviewed journal.

8. (samantha2020lipidaphosphoethanolamine pages 5-8): Ariela Samantha and Alice Vrielink. Lipid a phosphoethanolamine transferase: regulation, structure and immune response. Journal of Molecular Biology, 432:5184-5196, Aug 2020. URL: https://doi.org/10.1016/j.jmb.2020.04.022, doi:10.1016/j.jmb.2020.04.022. This article has 87 citations and is from a domain leading peer-reviewed journal.

9. (sun2021deepmutationalscanning pages 11-13): Zhizeng Sun and Timothy Palzkill. Deep mutational scanning reveals the active-site sequence requirements for the colistin antibiotic resistance enzyme mcr-1. Dec 2021. URL: https://doi.org/10.1128/mbio.02776-21, doi:10.1128/mbio.02776-21. This article has 22 citations and is from a domain leading peer-reviewed journal.

10. (thai2023phosphoethanolaminetransferasesas pages 16-18): Van C. Thai, Keith A. Stubbs, Mitali Sarkar-Tyson, and Charlene M. Kahler. Phosphoethanolamine transferases as drug discovery targets for therapeutic treatment of multi-drug resistant pathogenic gram-negative bacteria. Antibiotics, 12:1382, Aug 2023. URL: https://doi.org/10.3390/antibiotics12091382, doi:10.3390/antibiotics12091382. This article has 13 citations.

11. (sun2021deepmutationalscanning pages 5-7): Zhizeng Sun and Timothy Palzkill. Deep mutational scanning reveals the active-site sequence requirements for the colistin antibiotic resistance enzyme mcr-1. Dec 2021. URL: https://doi.org/10.1128/mbio.02776-21, doi:10.1128/mbio.02776-21. This article has 22 citations and is from a domain leading peer-reviewed journal.

12. (thai2023phosphoethanolaminetransferasesas pages 14-16): Van C. Thai, Keith A. Stubbs, Mitali Sarkar-Tyson, and Charlene M. Kahler. Phosphoethanolamine transferases as drug discovery targets for therapeutic treatment of multi-drug resistant pathogenic gram-negative bacteria. Antibiotics, 12:1382, Aug 2023. URL: https://doi.org/10.3390/antibiotics12091382, doi:10.3390/antibiotics12091382. This article has 13 citations.

13. (thai2023phosphoethanolaminetransferasesas pages 12-14): Van C. Thai, Keith A. Stubbs, Mitali Sarkar-Tyson, and Charlene M. Kahler. Phosphoethanolamine transferases as drug discovery targets for therapeutic treatment of multi-drug resistant pathogenic gram-negative bacteria. Antibiotics, 12:1382, Aug 2023. URL: https://doi.org/10.3390/antibiotics12091382, doi:10.3390/antibiotics12091382. This article has 13 citations.

14. (thai2023phosphoethanolaminetransferasesas pages 7-9): Van C. Thai, Keith A. Stubbs, Mitali Sarkar-Tyson, and Charlene M. Kahler. Phosphoethanolamine transferases as drug discovery targets for therapeutic treatment of multi-drug resistant pathogenic gram-negative bacteria. Antibiotics, 12:1382, Aug 2023. URL: https://doi.org/10.3390/antibiotics12091382, doi:10.3390/antibiotics12091382. This article has 13 citations.

15. (schumann2024themultifacetedroles pages 4-6): Anna Schumann, Ahmed Gaballa, and Martin Wiedmann. The multifaceted roles of phosphoethanolamine-modified lipopolysaccharides: from stress response and virulence to cationic antimicrobial resistance. Microbiology and Molecular Biology Reviews, Dec 2024. URL: https://doi.org/10.1128/mmbr.00193-23, doi:10.1128/mmbr.00193-23. This article has 13 citations and is from a domain leading peer-reviewed journal.

16. (anderson2020theescherichiacoli pages 16-18): Alexander C. Anderson, Alysha J.N. Burnett, Lana Hiscock, Kenneth E. Maly, and Joel T. Weadge. The escherichia coli cellulose synthase subunit g (bcsg) is a zn2+-dependent phosphoethanolamine transferase. Journal of Biological Chemistry, 295:6225-6235, May 2020. URL: https://doi.org/10.1074/jbc.ra119.011668, doi:10.1074/jbc.ra119.011668. This article has 36 citations and is from a domain leading peer-reviewed journal.

17. (zinkle2025mechanisticbasisof pages 3-4): Allen P. Zinkle, Mariana Bunoro Batista, Carmen M. Herrera, Satchal K. Erramilli, Brian Kloss, Khuram U. Ashraf, Kamil Nosol, Guozhi Zhang, Rosemary J. Cater, Michael T. Marty, Anthony A. Kossiakoff, M. Stephen Trent, Rie Nygaard, Phillip J. Stansfeld, and Filippo Mancia. Mechanistic basis of antimicrobial resistance mediated by the phosphoethanolamine transferase mcr-1. Nature Communications, Nov 2025. URL: https://doi.org/10.1038/s41467-025-65515-3, doi:10.1038/s41467-025-65515-3. This article has 4 citations and is from a highest quality peer-reviewed journal.

Artifacts

• Edison artifact artifact-00
• Edison artifact artifact-01