The research report should be a detailed narrative explaining the function, biological processes, and localization of the gene product. Citations should be given for all claims.

You should prioritize authoritative reviews and primary scientific literature when conducting research. You can supplement
this with annotations you find in gene/protein databases, but these can be outdated or inaccurate.

We are specifically interested in the primary function of the gene - for enzymes, what reaction is catalyzed, and what is the substrate specificity? For transporters, what is the substrate? For structural proteins or adapters, what is the broader structural role? For signaling molecules, what is the role in the pathway.

We are interested in where in or outside the cell the gene product carries out its function.

We are also interested in the signaling or biochemical pathways in which the gene functions. We are less interested in broad pleiotropic effects, except where these elucidate the precise role.

Include evidence where possible. We are interested in both experimental evidence as well as inference from structure, evolution, or bioinformatic analysis. Precise studies should be prioritized over high-throughput, where available.

Research Report: fcs (PP_3356; UniProt Q88HK0) in Pseudomonas putida KT2440

1) Mandatory identity verification (avoid “fcs” ambiguity)

The UniProt-provided target (Q88HK0) is annotated as feruloyl-CoA synthetase (EC 6.2.1.34) encoded by fcs with ordered locus name PP_3356 in Pseudomonas putida KT2440. Primary literature pathway mapping in P. putida KT2440 explicitly labels fcs (PP_3356) at the ferulic acid → feruloyl-CoA step and places it alongside neighboring pathway genes (including vdh (PP_3357) and PP_3358) in the same ferulate/coumarate catabolic branch, matching the intended identity (garciahidalgo2020vanillinproductionin media cfb4b31e, garciahidalgo2020vanillinproductionin pages 10-11).

2) Key concepts, definitions, and current understanding (2023–2024 prioritized)

2.1 Feruloyl-CoA synthetase (Fcs): biochemical role

Feruloyl-CoA synthetase (Fcs) is an acid-thiol ligase / acyl-CoA synthetase in the AMP-binding/ANL enzyme superfamily that catalyzes ATP-dependent CoA thioesterification of hydroxycinnamic acids, including ferulate and p-coumarate, yielding the corresponding CoA thioesters that are required for downstream catabolism and for some engineered biosyntheses (tramontina2023sustainablebiosyntheticpathways pages 4-6, incha2020leveraginghostmetabolism pages 3-4). In P. putida KT2440 specifically, functional work supports that Fcs activates p-coumarate to coumaroyl-CoA (coumaroyl-CoA detected in cells expressing fcs) (incha2023excavatingthegenome pages 86-90, incha2023excavatingthegenomea pages 86-90) and that ferulate is converted to vanillin via Fcs + Ech (ruhl2025productionofvanillin pages 1-3).

Note on reaction stoichiometry: The retrieved excerpts consistently describe ATP-dependent CoA thioesterification but do not provide an explicit balanced equation specifying whether ATP is converted to AMP + PPi vs ADP + Pi for P. putida Fcs (tramontina2023sustainablebiosyntheticpathways pages 4-6, goncalves2022applyingbiochemicaland pages 2-4). Therefore, the canonical EC 6.2.1.34 stoichiometry is not quoted verbatim from the retrieved full-text evidence here.

2.2 Pathway context in KT2440: hydroxycinnamate “peripheral” branch into the β-ketoadipate funnel

A KT2440 pathway map places Fcs (PP_3356) as the entry step for ferulate and p-coumarate catabolism: hydroxycinnamate → CoA-thioester → aromatic aldehyde (vanillin or 4-hydroxybenzaldehyde) → aromatic acid (vanillic acid or 4-hydroxybenzoic acid) → protocatechuate → β-ketoadipate pathway (garciahidalgo2020vanillinproductionin media cfb4b31e). More recent systems/engineering studies describe that native KT2440 converts ferulic acid onward toward protocatechuate and ultimately central metabolism via the β-ketoadipate pathway, and they identify Fcs and Ech as the two key enzymes mediating ferulate-to-vanillin conversion (ruhl2025productionofvanillin pages 1-3).

2.3 Genomic neighborhood and pathway gene set

In P. putida KT2440, fcs (PP_3356) is shown in close association with vdh (PP_3357) and PP_3358 in the ferulate catabolic branch on a gene/pathway map, supporting a local gene neighborhood for hydroxycinnamate utilization (garciahidalgo2020vanillinproductionin media cfb4b31e, garciahidalgo2020vanillinproductionin pages 10-11).

3) Evidence-based functional annotation of Fcs (Q88HK0)

3.1 Supported substrates and substrate scope

Organism-specific evidence (KT2440):
- p-Coumarate → coumaroyl-CoA: Coumaroyl-CoA is directly detected by LC–MS/MS (MRM transition 914.16→407.16) from intracellular extracts of P. putida expressing pBADT-fcs, providing in vivo evidence that Fcs activates p-coumarate to the CoA ester (incha2023excavatingthegenome pages 86-90, incha2023excavatingthegenomea pages 86-90).
- Ferulate entry into the ferulate → vanillin route: Multiple studies place Fcs as the initiating activation step for ferulate metabolism toward vanillin in KT2440 (ruhl2025productionofvanillin pages 1-3, garciahidalgo2020vanillinproductionin media cfb4b31e).

Broader homolog evidence (supporting inference, not KT2440-specific): Metagenome-derived Fcs homologs experimentally accept ferulic, p-coumaric, caffeic, and sinapic acids as substrates, indicating that bacterial Fcs enzymes can exhibit broad hydroxycinnamate specificity (goncalves2022applyingbiochemicaland pages 2-4, goncalves2022applyingbiochemicaland pages 6-9). This supports plausibility of broader substrate scope for Q88HK0 but does not itself prove it in KT2440.

3.2 Functional genomics and phenotypes: why Fcs matters physiologically

RB-TnSeq/fitness analyses and targeted experiments support Fcs as a key determinant of hydroxycinnamate utilization and stress phenotypes:
- A 2024 machine-learning analysis of RB-TnSeq fitness data (ICA) identified a hydroxycinnamate catabolism module (fModule_14) that explicitly groups Fcs:Ech:Vdh and links hydroxycinnamate utilization to assimilation of acetyl-CoA output and glyoxylate shunt dependence (glcB/malate synthase) (borchert2024machinelearninganalysis pages 6-7, borchert2024machinelearninganalysis pages 2-4).
- Accumulation of the activated CoA intermediate can be inhibitory: engineering work reports that coumaroyl-CoA production inhibits growth (increased lag when fcs is expressed with coumarate present) and that coumaroyl-CoA peaks transiently after substrate addition, consistent with toxicity if downstream consumption is limiting (incha2020leveraginghostmetabolism pages 3-4, incha2023excavatingthegenomea pages 8-13).

4) Localization: where Fcs carries out its function

4.1 Intracellular/cytosolic activity (best-supported)

Fcs produces CoA-thioesters that were extracted and quantified from cell pellets/intracellular extracts in KT2440 expressing fcs (incha2023excavatingthegenome pages 86-90, incha2023excavatingthegenomea pages 86-90). This supports the conclusion that the primary catalytic activity relevant to hydroxycinnamate activation is cell-associated (likely cytosolic), consistent with CoA metabolism and canonical bacterial acyl-CoA synthetases.

4.2 Outer membrane vesicles (OMVs): emerging spatial biology (uncertain for Fcs specifically)

A 2020 PNAS study showed P. putida secretes outer membrane vesicles (OMVs) during growth on lignin-rich media; OMV proteomics and functional assays support that OMVs can harbor enzymatically active aromatic catabolic proteins (salvachua2020outermembranevesicles pages 1-2, salvachua2020outermembranevesicles pages 4-6). The study’s β-ketoadipate-pathway heat map/legend lists Fcs, Ech, and Vdh among pathway enzymes detected in extracellular fractions/OMV-related analyses (salvachua2020outermembranevesicles pages 6-7). However, the retrieved excerpts do not provide unambiguous protein-by-protein enrichment values for Fcs in OMVs versus vesicle-free supernatant, so direct OMV localization of Fcs remains unresolved from the currently retrieved text (salvachua2020outermembranevesicles pages 6-7, salvachua2020outermembranevesicles pages 2-4).

5) Regulation and pathway control (recent evidence)

A key 2023 development is experimental evidence that carbon catabolite repression influences the hydroxycinnamate entry branch. In engineered KT2440 derivatives, deleting the global regulator Crc increased utilization rates of p-coumarate and ferulate, and this was attributed to derepression of fcs, ech, and vdh, which were reported to have putative Crc binding sites near their initiation codons (werner2023ligninconversionto pages 2-4). This connects fcs function to a broader regulatory program controlling aromatic assimilation.

6) Current applications and real-world implementations (2023–2024 prioritized)

6.1 Lignin-to-β-ketoadipate funneling (Science Advances, 2023)

Werner et al. (published Sep 2023) engineered P. putida KT2440 to funnel lignin-related aromatics (including hydroxycinnamates) to β-ketoadipic acid. They report high, industrially relevant metrics, including 44.5 ± 1.85 g/L β-ketoadipate at 0.85 ± 0.04 g/L/h in optimized fed-batch and 25 g/L at 0.66 g/L/h from corn stover–derived lignin streams, corresponding to 0.10 g product/g lignin (https://doi.org/10.1126/sciadv.adj0053) (werner2023ligninconversionto pages 2-4, werner2023ligninconversionto pages 1-2). Because fcs/ech/vdh mediate the initial hydroxycinnamate breakdown, pathway derepression via ∆crc is mechanistically tied to improved conversion of ferulate/p-coumarate into the β-ketoadipate funnel (werner2023ligninconversionto pages 2-4).

6.2 Hydrolysate valorization to protocatechuate (Molecules, 2024)

Jin et al. (published Mar 2024) quantified lignin-derived aromatics in corncob hydrolysates, reporting ferulic acid and p-coumaric acid levels exceeding 160/120 mg/L and 530/250 mg/L, respectively, in two hydrolysates. They engineered KT2440 to accumulate protocatechuate by blocking its cleavage and reported PCA titers up to 433.72 mg/L (https://doi.org/10.3390/molecules29071555) (jin2024biologicalvalorizationof pages 1-2, jin2024biologicalvalorizationof pages 2-4). In degradation assays, KT2440 degraded phenolics (including ferulate and p-coumarate) largely within 24 h, with H-type aromatics consumed within 12 h (jin2024biologicalvalorizationof pages 2-4). They also report strong co-substrate effects on p-coumarate consumption (e.g., at 6 h p-CA consumption falling from 73.62% to 41.19% with glucose) (jin2024biologicalvalorizationof pages 4-7). These data are directly relevant to deploying the fcs branch in realistic lignocellulosic streams.

6.3 Functional genomics for strain engineering (mSystems, 2024)

Borchert et al. (published Mar 2024) provide a data-driven approach to improve hydroxycinnamate processing by identifying functional modules and tolerance genes from RB-TnSeq data. Their hydroxycinnamate module (fModule_14) groups Fcs:Ech:Vdh, supporting a systems-level rationale for cofactor/anaplerotic engineering (e.g., glyoxylate shunt via glcB) to improve growth and conversion on ferulate/p-coumarate (https://doi.org/10.1128/msystems.00942-23) (borchert2024machinelearninganalysis pages 6-7, borchert2024machinelearninganalysis pages 2-4).

7) Expert interpretation and authoritative synthesis

Recent synthesis of hydroxycinnamate bioconversion highlights Fcs (EC 6.2.1.34) as the activation gate for CoA-dependent microbial hydroxycinnamate pathways and as a common starting point for producing value-added compounds such as vanillin (Tramontina et al., May 2023, https://doi.org/10.1007/s00253-023-12571-8) (tramontina2023sustainablebiosyntheticpathways pages 4-6). The P. putida KT2440 literature further suggests that a major constraint is not merely enzyme presence but flux balancing and avoidance of toxic intermediate accumulation (e.g., coumaroyl-CoA), motivating engineering strategies that coordinate Fcs with downstream consumption (incha2020leveraginghostmetabolism pages 3-4).

8) Summary evidence map (structured)

The following table consolidates key functional claims, pathway placement, regulation, localization, and application metrics with evidence pointers.

Table: This table summarizes the functional annotation of Pseudomonas putida KT2440 fcs (PP_3356; UniProt Q88HK0), including its enzymatic role, pathway placement, regulation, localization, and engineering relevance. It is useful as a compact evidence map linking gene identity to experimentally supported function and recent applications.

9) Gaps and recommendations for annotation confidence

1. Explicit balanced reaction stoichiometry (AMP + PPi) was not found verbatim in the retrieved KT2440 full-text excerpts; if required for strict annotation, it should be extracted from a primary enzymology paper or curated enzyme database entry that explicitly states the stoichiometry for EC 6.2.1.34 in Pseudomonas Fcs.
2. OMV localization of Fcs is plausible but not definitively established from the retrieved excerpts; confirming whether PP_3356 is enriched in OMVs versus other fractions would require consulting the deposited proteomics datasets and/or the full figure tables from Salvachúa et al. (PXD016114 / MassIVE MSV000084524, MSV000084506) (salvachua2020outermembranevesicles pages 1-2).

Key URLs (with publication dates)

• Werner et al., Sep 2023, Science Advances: https://doi.org/10.1126/sciadv.adj0053 (werner2023ligninconversionto pages 2-4)
• Borchert et al., Mar 2024, mSystems: https://doi.org/10.1128/msystems.00942-23 (borchert2024machinelearninganalysis pages 6-7)
• Jin et al., Mar 2024, Molecules: https://doi.org/10.3390/molecules29071555 (jin2024biologicalvalorizationof pages 1-2)
• Tramontina et al., May 2023, Applied Microbiology and Biotechnology: https://doi.org/10.1007/s00253-023-12571-8 (tramontina2023sustainablebiosyntheticpathways pages 4-6)
• Salvachúa et al., Apr 2020, PNAS: https://doi.org/10.1073/pnas.1921073117 (salvachua2020outermembranevesicles pages 1-2)
• García-Hidalgo et al., Mar 2020, Applied and Environmental Microbiology: https://doi.org/10.1128/AEM.02442-19 (garciahidalgo2020vanillinproductionin pages 10-11)

References

1. (garciahidalgo2020vanillinproductionin media cfb4b31e): Javier García-Hidalgo, Daniel P. Brink, Krithika Ravi, Catherine J. Paul, Gunnar Lidén, and Marie F. Gorwa-Grauslund. Vanillin production in pseudomonas : whole-genome sequencing of pseudomonas sp. strain 9.1 and reannotation of pseudomonas putida cala as a vanillin reductase. Mar 2020. URL: https://doi.org/10.1128/aem.02442-19, doi:10.1128/aem.02442-19. This article has 43 citations and is from a peer-reviewed journal.

2. (garciahidalgo2020vanillinproductionin pages 10-11): Javier García-Hidalgo, Daniel P. Brink, Krithika Ravi, Catherine J. Paul, Gunnar Lidén, and Marie F. Gorwa-Grauslund. Vanillin production in pseudomonas : whole-genome sequencing of pseudomonas sp. strain 9.1 and reannotation of pseudomonas putida cala as a vanillin reductase. Mar 2020. URL: https://doi.org/10.1128/aem.02442-19, doi:10.1128/aem.02442-19. This article has 43 citations and is from a peer-reviewed journal.

3. (tramontina2023sustainablebiosyntheticpathways pages 4-6): Robson Tramontina, Iara Ciancaglini, Ellen K. B. Roman, Micaela G. Chacón, Thamy L. R. Corrêa, Neil Dixon, Timothy D. H. Bugg, and Fabio Marcio Squina. Sustainable biosynthetic pathways to value-added bioproducts from hydroxycinnamic acids. Applied Microbiology and Biotechnology, 107:4165-4185, May 2023. URL: https://doi.org/10.1007/s00253-023-12571-8, doi:10.1007/s00253-023-12571-8. This article has 20 citations and is from a domain leading peer-reviewed journal.

4. (incha2020leveraginghostmetabolism pages 3-4): Matthew R. Incha, Mitchell G. Thompson, Jacquelyn M. Blake-Hedges, Yuzhong Liu, Allison N. Pearson, Matthias Schmidt, Jennifer W. Gin, Christopher J. Petzold, Adam M. Deutschbauer, and Jay D. Keasling. Leveraging host metabolism for bisdemethoxycurcumin production in pseudomonas putida. Jun 2020. URL: https://doi.org/10.1016/j.mec.2019.e00119, doi:10.1016/j.mec.2019.e00119. This article has 76 citations and is from a peer-reviewed journal.

5. (incha2023excavatingthegenome pages 86-90): MR Incha. Excavating the genome mine of pseudomonas putida kt2440. Unknown journal, 2023.

6. (incha2023excavatingthegenomea pages 86-90): MR Incha. Excavating the genome mine of pseudomonas putida kt2440. Unknown journal, 2023.

7. (ruhl2025productionofvanillin pages 1-3): Ilona A. Ruhl, Sean P. Woodworth, Stefan J. Haugen, Hannah M. Alt, Gregg T. Beckham, and Christopher W. Johnson. Production of vanillin from ferulic acid by pseudomonas putida kt2440 using metabolic engineering and in situ product recovery. Microbial Biotechnology, May 2025. URL: https://doi.org/10.1111/1751-7915.70152, doi:10.1111/1751-7915.70152. This article has 12 citations and is from a peer-reviewed journal.

8. (goncalves2022applyingbiochemicaland pages 2-4): Thiago Augusto Gonçalves, Victoria Sodré, Stephanie Nemesio da Silva, Nathalia Vilela, Geizecler Tomazetto, Juscemácia Nascimento Araujo, João Renato C. Muniz, Taícia Pacheco Fill, André Damasio, Wanius Garcia, and Fabio Marcio Squina. Applying biochemical and structural characterization of hydroxycinnamate catabolic enzymes from soil metagenome for lignin valorization strategies. Applied Microbiology and Biotechnology, 106:2503-2516, Mar 2022. URL: https://doi.org/10.1007/s00253-022-11885-3, doi:10.1007/s00253-022-11885-3. This article has 10 citations and is from a domain leading peer-reviewed journal.

9. (goncalves2022applyingbiochemicaland pages 6-9): Thiago Augusto Gonçalves, Victoria Sodré, Stephanie Nemesio da Silva, Nathalia Vilela, Geizecler Tomazetto, Juscemácia Nascimento Araujo, João Renato C. Muniz, Taícia Pacheco Fill, André Damasio, Wanius Garcia, and Fabio Marcio Squina. Applying biochemical and structural characterization of hydroxycinnamate catabolic enzymes from soil metagenome for lignin valorization strategies. Applied Microbiology and Biotechnology, 106:2503-2516, Mar 2022. URL: https://doi.org/10.1007/s00253-022-11885-3, doi:10.1007/s00253-022-11885-3. This article has 10 citations and is from a domain leading peer-reviewed journal.

10. (borchert2024machinelearninganalysis pages 6-7): Andrew J. Borchert, Alissa C. Bleem, Hyun Gyu Lim, Kevin Rychel, Keven D. Dooley, Zoe A. Kellermyer, Tracy L. Hodges, Bernhard O. Palsson, and Gregg T. Beckham. Machine learning analysis of rb-tnseq fitness data predicts functional gene modules in pseudomonas putida kt2440. mSystems, Mar 2024. URL: https://doi.org/10.1128/msystems.00942-23, doi:10.1128/msystems.00942-23. This article has 13 citations and is from a peer-reviewed journal.

11. (borchert2024machinelearninganalysis pages 2-4): Andrew J. Borchert, Alissa C. Bleem, Hyun Gyu Lim, Kevin Rychel, Keven D. Dooley, Zoe A. Kellermyer, Tracy L. Hodges, Bernhard O. Palsson, and Gregg T. Beckham. Machine learning analysis of rb-tnseq fitness data predicts functional gene modules in pseudomonas putida kt2440. mSystems, Mar 2024. URL: https://doi.org/10.1128/msystems.00942-23, doi:10.1128/msystems.00942-23. This article has 13 citations and is from a peer-reviewed journal.

12. (incha2023excavatingthegenomea pages 8-13): MR Incha. Excavating the genome mine of pseudomonas putida kt2440. Unknown journal, 2023.

13. (salvachua2020outermembranevesicles pages 1-2): Davinia Salvachúa, Allison Z. Werner, Isabel Pardo, Martyna Michalska, Brenna A. Black, Bryon S. Donohoe, Stefan J. Haugen, Rui Katahira, Sandra Notonier, Kelsey J. Ramirez, Antonella Amore, Samuel O. Purvine, Erika M. Zink, Paul E. Abraham, Richard J. Giannone, Suresh Poudel, Philip D. Laible, Robert L. Hettich, and Gregg T. Beckham. Outer membrane vesicles catabolize lignin-derived aromatic compounds in pseudomonas putida kt2440. Proceedings of the National Academy of Sciences, 117:9302-9310, Apr 2020. URL: https://doi.org/10.1073/pnas.1921073117, doi:10.1073/pnas.1921073117. This article has 158 citations and is from a highest quality peer-reviewed journal.

14. (salvachua2020outermembranevesicles pages 4-6): Davinia Salvachúa, Allison Z. Werner, Isabel Pardo, Martyna Michalska, Brenna A. Black, Bryon S. Donohoe, Stefan J. Haugen, Rui Katahira, Sandra Notonier, Kelsey J. Ramirez, Antonella Amore, Samuel O. Purvine, Erika M. Zink, Paul E. Abraham, Richard J. Giannone, Suresh Poudel, Philip D. Laible, Robert L. Hettich, and Gregg T. Beckham. Outer membrane vesicles catabolize lignin-derived aromatic compounds in pseudomonas putida kt2440. Proceedings of the National Academy of Sciences, 117:9302-9310, Apr 2020. URL: https://doi.org/10.1073/pnas.1921073117, doi:10.1073/pnas.1921073117. This article has 158 citations and is from a highest quality peer-reviewed journal.

15. (salvachua2020outermembranevesicles pages 6-7): Davinia Salvachúa, Allison Z. Werner, Isabel Pardo, Martyna Michalska, Brenna A. Black, Bryon S. Donohoe, Stefan J. Haugen, Rui Katahira, Sandra Notonier, Kelsey J. Ramirez, Antonella Amore, Samuel O. Purvine, Erika M. Zink, Paul E. Abraham, Richard J. Giannone, Suresh Poudel, Philip D. Laible, Robert L. Hettich, and Gregg T. Beckham. Outer membrane vesicles catabolize lignin-derived aromatic compounds in pseudomonas putida kt2440. Proceedings of the National Academy of Sciences, 117:9302-9310, Apr 2020. URL: https://doi.org/10.1073/pnas.1921073117, doi:10.1073/pnas.1921073117. This article has 158 citations and is from a highest quality peer-reviewed journal.

16. (salvachua2020outermembranevesicles pages 2-4): Davinia Salvachúa, Allison Z. Werner, Isabel Pardo, Martyna Michalska, Brenna A. Black, Bryon S. Donohoe, Stefan J. Haugen, Rui Katahira, Sandra Notonier, Kelsey J. Ramirez, Antonella Amore, Samuel O. Purvine, Erika M. Zink, Paul E. Abraham, Richard J. Giannone, Suresh Poudel, Philip D. Laible, Robert L. Hettich, and Gregg T. Beckham. Outer membrane vesicles catabolize lignin-derived aromatic compounds in pseudomonas putida kt2440. Proceedings of the National Academy of Sciences, 117:9302-9310, Apr 2020. URL: https://doi.org/10.1073/pnas.1921073117, doi:10.1073/pnas.1921073117. This article has 158 citations and is from a highest quality peer-reviewed journal.

17. (werner2023ligninconversionto pages 2-4): Allison Z. Werner, William T. Cordell, Ciaran W. Lahive, Bruno C. Klein, Christine A. Singer, Eric C. D. Tan, Morgan A. Ingraham, Kelsey J. Ramirez, Dong Hyun Kim, Jacob Nedergaard Pedersen, Christopher W. Johnson, Brian F. Pfleger, Gregg T. Beckham, and Davinia Salvachúa. Lignin conversion to β-ketoadipic acid by pseudomonas putida via metabolic engineering and bioprocess development. Science Advances, Sep 2023. URL: https://doi.org/10.1126/sciadv.adj0053, doi:10.1126/sciadv.adj0053. This article has 88 citations and is from a highest quality peer-reviewed journal.

18. (werner2023ligninconversionto pages 1-2): Allison Z. Werner, William T. Cordell, Ciaran W. Lahive, Bruno C. Klein, Christine A. Singer, Eric C. D. Tan, Morgan A. Ingraham, Kelsey J. Ramirez, Dong Hyun Kim, Jacob Nedergaard Pedersen, Christopher W. Johnson, Brian F. Pfleger, Gregg T. Beckham, and Davinia Salvachúa. Lignin conversion to β-ketoadipic acid by pseudomonas putida via metabolic engineering and bioprocess development. Science Advances, Sep 2023. URL: https://doi.org/10.1126/sciadv.adj0053, doi:10.1126/sciadv.adj0053. This article has 88 citations and is from a highest quality peer-reviewed journal.

19. (jin2024biologicalvalorizationof pages 1-2): Xinzhu Jin, Xiaoxia Li, Lihua Zou, Zhaojuan Zheng, and Jia Ouyang. Biological valorization of lignin-derived aromatics in hydrolysate to protocatechuic acid by engineered pseudomonas putida kt2440. Molecules, 29:1555, Mar 2024. URL: https://doi.org/10.3390/molecules29071555, doi:10.3390/molecules29071555. This article has 15 citations.

20. (jin2024biologicalvalorizationof pages 2-4): Xinzhu Jin, Xiaoxia Li, Lihua Zou, Zhaojuan Zheng, and Jia Ouyang. Biological valorization of lignin-derived aromatics in hydrolysate to protocatechuic acid by engineered pseudomonas putida kt2440. Molecules, 29:1555, Mar 2024. URL: https://doi.org/10.3390/molecules29071555, doi:10.3390/molecules29071555. This article has 15 citations.

21. (jin2024biologicalvalorizationof pages 4-7): Xinzhu Jin, Xiaoxia Li, Lihua Zou, Zhaojuan Zheng, and Jia Ouyang. Biological valorization of lignin-derived aromatics in hydrolysate to protocatechuic acid by engineered pseudomonas putida kt2440. Molecules, 29:1555, Mar 2024. URL: https://doi.org/10.3390/molecules29071555, doi:10.3390/molecules29071555. This article has 15 citations.

22. (zhou2025quantitativeanalysisof pages 5-8): Nanqing Zhou, Rebecca A. Wilkes, Xinyu Chen, Kelly P. Teitel, James A. Belgrave, Gregg T. Beckham, Allison Z. Werner, Yanbao Yu, and Ludmilla Aristilde. Quantitative analysis of coupled carbon and energy metabolism for lignin carbon utilization in pseudomonas putida. bioRxiv, Mar 2025. URL: https://doi.org/10.1101/2025.03.24.645021, doi:10.1101/2025.03.24.645021. This article has 3 citations.

23. (incha2023excavatingthegenome pages 8-13): MR Incha. Excavating the genome mine of pseudomonas putida kt2440. Unknown journal, 2023.

Artifacts

• Edison artifact artifact-00
  

Citations

1. ruhl2025productionofvanillin pages 1-3
2. salvachua2020outermembranevesicles pages 6-7
3. werner2023ligninconversionto pages 2-4
4. jin2024biologicalvalorizationof pages 2-4
5. jin2024biologicalvalorizationof pages 4-7
6. tramontina2023sustainablebiosyntheticpathways pages 4-6
7. incha2020leveraginghostmetabolism pages 3-4
8. zhou2025quantitativeanalysisof pages 5-8
9. incha2023excavatingthegenomea pages 86-90
10. salvachua2020outermembranevesicles pages 1-2
11. borchert2024machinelearninganalysis pages 6-7
12. jin2024biologicalvalorizationof pages 1-2
13. garciahidalgo2020vanillinproductionin pages 10-11
14. incha2023excavatingthegenome pages 86-90
15. goncalves2022applyingbiochemicaland pages 2-4
16. goncalves2022applyingbiochemicaland pages 6-9
17. borchert2024machinelearninganalysis pages 2-4
18. incha2023excavatingthegenomea pages 8-13
19. salvachua2020outermembranevesicles pages 4-6
20. salvachua2020outermembranevesicles pages 2-4
21. werner2023ligninconversionto pages 1-2
22. incha2023excavatingthegenome pages 8-13
23. https://doi.org/10.1126/sciadv.adj0053
24. https://doi.org/10.3390/molecules29071555
25. https://doi.org/10.1128/msystems.00942-23
26. https://doi.org/10.1007/s00253-023-12571-8
27. https://doi.org/10.1128/AEM.02442-19
28. https://doi.org/10.1016/j.mec.2019.e00119
29. https://doi.org/10.1007/s00253-022-11885-3
30. https://doi.org/10.1111/1751-7915.70152
31. https://doi.org/10.1101/2025.03.24.645021
32. https://doi.org/10.1073/pnas.1921073117
33. https://doi.org/10.1128/aem.02442-19,
34. https://doi.org/10.1007/s00253-023-12571-8,
35. https://doi.org/10.1016/j.mec.2019.e00119,
36. https://doi.org/10.1111/1751-7915.70152,
37. https://doi.org/10.1007/s00253-022-11885-3,
38. https://doi.org/10.1128/msystems.00942-23,
39. https://doi.org/10.1073/pnas.1921073117,
40. https://doi.org/10.1126/sciadv.adj0053,
41. https://doi.org/10.3390/molecules29071555,
42. https://doi.org/10.1101/2025.03.24.645021,

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: PSEPK
• Taxon: Pseudomonas putida (strain ATCC 47054 / DSM 6125 / CFBP 8728 / NCIMB 11950 / KT2440) (NCBITaxon:160488)
• Gene directory: fcs
• Gene symbol: fcs
• UniProt accession: Q88HK0

Focus

• Focus type: function_assignment
• Hypothesis slug: function-hypothesis-go-0031956
• Source file: genes/PSEPK/fcs/fcs-ai-review.yaml
• Source selector: existing_annotations[2].function_hypothesis

Seed Hypothesis (propagated by TreeGrafter/PANTHER)

fcs has medium-chain fatty acid-CoA ligase activity (GO:0031956).

Term and Decision Context

• Term: medium-chain fatty acid-CoA ligase activity (GO:0031956)
• Evidence type: IEA
• Original reference: GO_REF:0000118

Reference Context

• GO_REF:0000118

Source Context YAML

term:
  id: GO:0031956
  label: medium-chain fatty acid-CoA ligase activity
evidence_type: IEA
original_reference_id: GO_REF:0000118

Research Objective

Decide whether fcs 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: fcs (Q88HK0, Pseudomonas putida KT2440)

Executive Judgment

Verdict: Refuted — Wrong subfamily (mis-placed); Failure mode 3.

The propagated GO term GO:0031956 (medium-chain fatty acid-CoA ligase activity) is incorrect for fcs (Q88HK0) from Pseudomonas putida KT2440. This represents a clear case of within-superfamily mis-placement (failure mode 3): the TreeGrafter/PANTHER pipeline grafted fcs onto a PANTHER node representing medium-chain acyl-CoA synthetases (ACSM subfamily) within the ANL (acyl-CoA synthetase / NRPS adenylation / luciferase) superfamily, when fcs actually belongs to the 4-coumarate:CoA ligase (4CL) / hydroxycinnamate-CoA synthetase subfamily — a functionally distinct, sibling branch of the same adenylate-forming enzyme superfamily.

Most decisive evidence: An fcs deletion mutant in P. putida F1 cannot grow on ferulic acid, caffeic acid, or p-coumaric acid as sole carbon sources, and growth is restored by complementation with fcs in trans (hughes2017pseudomonasputidaf1 pages 3-4). This directly demonstrates that fcs is essential for hydroxycinnamic acid (aromatic) catabolism, not medium-chain fatty acid metabolism. The characterized closest homolog CouL from Rhodococcus jostii RHA1 thioesterifies ferulate, p-coumarate, caffeate, and dihydroferulate but shows no activity with non-aromatic substrates such as vanillate, benzoate, or aliphatic acids (otani2014characterizationofp pages 7-8, otani2014characterizationofp pages 8-9).

Independent Family/Function Assignment

Protein: fcs (Q88HK0, PP_3356), Pseudomonas putida KT2440

Most likely specific molecular function: Feruloyl-CoA synthetase / p-hydroxycinnamoyl-CoA synthetase activity. The enzyme catalyzes the ATP-dependent CoA ligation of hydroxycinnamic acids (ferulic acid, caffeic acid, p-coumaric acid) as the first committed step in the CoA-dependent, non-β-oxidative phenylpropenoid degradation pathway leading to vanillin and vanillate production (jimenez2002genomicanalysisof pages 9-10, ruhl2025productionofvanillin pages 1-3, zhou2020developmentofa pages 6-7).

Reaction catalyzed: Ferulic acid + CoA + ATP → Feruloyl-CoA + AMP + PPi

Recommended GO term: GO:0016207 (4-coumarate-CoA ligase activity) is the most appropriate existing GO molecular function term. This term describes the ATP-dependent CoA ligation of hydroxycinnamic acid substrates, which is precisely the reaction catalyzed by fcs. If a more specific child term for feruloyl-CoA synthetase activity exists or is created, that would be even more precise.

Characterized homolog basis: CouL from Rhodococcus jostii RHA1, which has been biochemically characterized with purified enzyme kinetics showing kcat/Km of ~370 mM⁻¹ s⁻¹ for p-coumarate, with activity on ferulate, caffeate, and dihydroferulate, but no detectable activity on sinapate, vanillate, or benzoate (otani2014characterizationofp pages 7-8, otani2014characterizationofp pages 8-9, otani2014characterizationofp pages 1-2).

Granularity relative to seed term: The correct term (GO:0016207 or a child thereof) is a sibling of GO:0031956 within the broader CoA ligase / acid-thiol ligase activity hierarchy — both are acyl-CoA synthetase activities but act on entirely different substrate classes (aromatic hydroxycinnamic acids vs. aliphatic medium-chain fatty acids). This is not a granularity issue (too general or too specific) but rather a mis-placement onto the wrong neighboring subfamily.

Evidence Matrix

The following table summarizes all evidence items evaluated, showing that every independent line of evidence refutes the assignment of GO:0031956 to fcs.

Table: This table summarizes independent evidence evaluating whether P. putida fcs (Q88HK0) has medium-chain fatty acid-CoA ligase activity. The evidence consistently points instead to hydroxycinnamate/feruloyl-CoA synthetase function and a within-superfamily misannotation.

Active-Site / Placement Analysis

The following analysis compares the structural, catalytic, and SDR-profile features of fcs with characterized 4CL-type (aromatic CoA ligase) and ACSM-type (medium-chain fatty acid CoA ligase) enzymes within the ANL superfamily.

Table: This table compares P. putida fcs with aromatic 4CL-type ligases and medium-chain fatty acid-CoA ligases across catalytic motifs, substrate-pocket architecture, SDRs, and pathway evidence. It is useful for localizing the TreeGrafter error as a within-superfamily misplacement rather than loss of activity.

Summary of Active-Site Analysis

Catalytic residues: The ANL superfamily conserves ten signature motifs (A1–A10) across all members, including the invariant A10 lysine (Lys526 in Nt4CL2) required for adenylation, the P-loop for pyrophosphate coordination, His237 (A4 motif) essential for catalysis, and Thr336 (A5 motif) (li2015structuralbasisfor pages 4-5, li2015structuralbasisfor pages 5-6). These catalytic residues are shared between 4CL-type and ACSM-type enzymes as superfamily-level features. The fcs enzyme is catalytically active as demonstrated by genetic complementation (hughes2017pseudomonasputidaf1 pages 3-4). This rules out failure mode 2 (pseudo-enzyme / loss of activity).

Substrate-binding pocket architecture: The critical distinction between the two subfamilies lies in the substrate-binding pocket. 4CL-type enzymes have broader cavities optimized for flat, planar hydroxycinnamic acid substrates, with residues enabling π-stacking interactions (e.g., Tyr239 equivalent) and specific gate residues controlling ring-substitution tolerance (e.g., Val341 in Nt4CL2) (li2015structuralbasisfor pages 7-8, li2015structuralbasisfor pages 6-7). In contrast, ACSM/medium-chain fatty acid-CoA ligases use hydrophobic tunnel-like pockets that act as "molecular rulers" to select for aliphatic chain length, with bulkier residues along helices α7/α8 and strands β10/β12 that restrict tunnel dimensions (priyadarshan2018fattyacylampligases pages 7-9, khurana2010genomescaleprediction pages 10-13).

SDR-based classification: The Khurana et al. (2010) computational framework using 15 substrate-specificity-determining residues (SDRs) can distinguish 4CL, ACSM, ACSL (long-chain), AcCS (acetyl), and other subfamilies with sensitivity of 0.91–1.0 and specificity of 0.96–1.0 (khurana2010genomescaleprediction pages 6-8, khurana2010genomescaleprediction pages 4-6, khurana2010genomescaleprediction pages 3-4). Key discriminating positions include: position 210 (His in medium-chain vs. Asn in long-chain), position 234 (conserved His in 4CL-type), position 324 (Trp in small-substrate enzymes vs. Gly for larger substrates), and position 301 (khurana2010genomescaleprediction pages 6-8, khurana2010genomescaleprediction pages 10-13). These SDR profiles clearly separate the 4CL aromatic-CoA ligase subfamily from the ACSM fatty acid-CoA ligase subfamily.

Limitation: I was unable to programmatically perform a direct residue-by-residue alignment of Q88HK0 against characterized 4CL and ACSM sequences because the tools available do not include sequence retrieval and alignment capabilities. However, the genomic context, characterized homolog data, and substrate specificity evidence are individually and collectively decisive.

GO Curation Implications

Recommended curation action: Replace with sibling term.

• Remove: GO:0031956 (medium-chain fatty acid-CoA ligase activity) — this is a clear mis-annotation resulting from within-superfamily mis-placement by the TreeGrafter pipeline.
• Replace with: GO:0016207 (4-coumarate-CoA ligase activity) — this correctly captures the ATP-dependent CoA ligation of hydroxycinnamic acid substrates (ferulic acid, p-coumaric acid, caffeic acid) that fcs catalyzes (hughes2017pseudomonasputidaf1 pages 3-4, jimenez2002genomicanalysisof pages 9-10, otani2014characterizationofp pages 7-8).
• Evidence code: The replacement annotation could be supported by IDA (inferred from direct assay) based on the genetic complementation and growth phenotype data in P. putida F1 (hughes2017pseudomonasputidaf1 pages 3-4), combined with ISS (inferred from sequence/structural similarity) to CouL from Rhodococcus jostii RHA1, which has been directly characterized biochemically (otani2014characterizationofp pages 7-8, otani2014characterizationofp pages 8-9).
• Biological process terms such as GO:0019439 (aromatic compound catabolic process) and GO:0009058 (biosynthetic process) related to vanillin biosynthesis from ferulic acid should also be considered for annotation (jimenez2002genomicanalysisof pages 9-10, ruhl2025productionofvanillin pages 1-3, zhou2020developmentofa pages 6-7).

Conflicts, Knowledge Gaps, and Discriminating Tests

Source of the Error

The mis-annotation is a textbook example of the ANL superfamily confounding phylogenetic propagation pipelines. Both 4CL-type hydroxycinnamate-CoA ligases and ACSM-type medium-chain fatty acid-CoA ligases share the AMP-forming/AMP-binding domain (Pfam PF00501), the same two-step adenylation/thioesterification catalytic mechanism, and similar overall folds (priyadarshan2018fattyacylampligases pages 2-4, khurana2010genomescaleprediction pages 4-6). However, their substrate-binding pockets have diverged to accommodate fundamentally different substrate classes (aromatic vs. aliphatic), and their biological roles are unrelated (aromatic catabolism vs. lipid metabolism) (priyadarshan2018fattyacylampligases pages 7-9, li2015structuralbasisfor pages 7-8). TreeGrafter apparently grafted fcs onto the wrong branch of the PANTHER reference tree at the superfamily level rather than the subfamily level.

Key Uncertainties

1. Direct in vitro kinetics on purified fcs (Q88HK0): While the genetic evidence is strong and the characterized homolog CouL provides kinetic data, published purified-enzyme kinetics for the exact Q88HK0 protein were not identified in this search. A 2022 paper by dos Santos et al. (doi:10.1016/j.pep.2022.106109) and a 2025 Methods in Enzymology chapter by Sodré & Squina (doi:10.1016/bs.mie.2025.01.046) describe recombinant feruloyl-CoA synthetase characterization but were unobtainable.
2. SDR alignment not computed: I could not programmatically extract the 15 SDR positions from Q88HK0 and compare them to the 4CL vs. ACSM reference profiles from Khurana et al. (2010). This would be the most direct computational proof of subfamily placement.

Discriminating Tests to Resolve

1. Most efficient resolution: Extract SDR positions from Q88HK0 by aligning to the 1AMU structural template and compare the 15-residue profile against the 4CL, ACSM, and ACSL subfamily PSSMs from Khurana et al. (2010) (khurana2010genomescaleprediction pages 4-6, khurana2010genomescaleprediction pages 14-16). This would computationally confirm subfamily placement in minutes.
2. Direct enzymatic assay: Purify recombinant Q88HK0 and test CoA ligation activity with: (a) hydroxycinnamic acid substrates (ferulic, p-coumaric, caffeic acid), and (b) medium-chain fatty acids (hexanoic, octanoic, decanoic acid). Activity with the former and not the latter would definitively refute GO:0031956.
3. Negative control: Test whether an fcs deletion mutant shows any defect in medium-chain fatty acid utilization. Based on the available evidence, no such defect is expected.

No Conflicting Evidence

No evidence was found supporting the assignment of GO:0031956 to fcs. All nine independent lines of evidence (genetic, transcriptomic, pathway, structural, computational, metabolic engineering) consistently point to hydroxycinnamate-CoA synthetase activity. There is no paralog confusion concern because fcs (PP_3356) is the sole gene in the P. putida KT2440 genome encoding the ferulic acid catabolic pathway's CoA-activating step, and it resides in a well-characterized gene cluster with ech and vdh (jimenez2002genomicanalysisof pages 9-10, darrigo2019analysisofpseudomonas pages 14-16).

References

1. (hughes2017pseudomonasputidaf1 pages 3-4): Jonathan G. Hughes, Xiangsheng Zhang, Juanito V. Parales, Jayna L. Ditty, and Rebecca E. Parales. Pseudomonas putida f1 uses energy taxis to sense hydroxycinnamic acids. Oct 2017. URL: https://doi.org/10.1099/mic.0.000533, doi:10.1099/mic.0.000533. This article has 13 citations and is from a peer-reviewed journal.

2. (otani2014characterizationofp pages 7-8): Hiroshi Otani, Young-Eun Lee, Israël Casabon, and Lindsay D. Eltis. Characterization of p -hydroxycinnamate catabolism in a soil actinobacterium. Journal of Bacteriology, 196:4293-4303, Dec 2014. URL: https://doi.org/10.1128/jb.02247-14, doi:10.1128/jb.02247-14. This article has 74 citations and is from a peer-reviewed journal.

3. (otani2014characterizationofp pages 8-9): Hiroshi Otani, Young-Eun Lee, Israël Casabon, and Lindsay D. Eltis. Characterization of p -hydroxycinnamate catabolism in a soil actinobacterium. Journal of Bacteriology, 196:4293-4303, Dec 2014. URL: https://doi.org/10.1128/jb.02247-14, doi:10.1128/jb.02247-14. This article has 74 citations and is from a peer-reviewed journal.

4. (jimenez2002genomicanalysisof pages 9-10): José Ignacio Jiménez, Baltasar Miñambres, José Luis García, and Eduardo Díaz. Genomic analysis of the aromatic catabolic pathways from pseudomonas putida kt2440. Environmental microbiology, 4 12:824-41, Dec 2002. URL: https://doi.org/10.1046/j.1462-2920.2002.00370.x, doi:10.1046/j.1462-2920.2002.00370.x. This article has 704 citations and is from a domain leading peer-reviewed journal.

5. (ruhl2025productionofvanillin pages 1-3): Ilona A. Ruhl, Sean P. Woodworth, Stefan J. Haugen, Hannah M. Alt, Gregg T. Beckham, and Christopher W. Johnson. Production of vanillin from ferulic acid by pseudomonas putida kt2440 using metabolic engineering and in situ product recovery. Microbial Biotechnology, May 2025. URL: https://doi.org/10.1111/1751-7915.70152, doi:10.1111/1751-7915.70152. This article has 13 citations and is from a peer-reviewed journal.

6. (zhou2020developmentofa pages 6-7): Yueyue Zhou, Lu Lin, Heng Wang, Zhichao Zhang, Jizhong Zhou, and Nianzhi Jiao. Development of a crispr/cas9n-based tool for metabolic engineering of pseudomonas putida for ferulic acid-to-polyhydroxyalkanoate bioconversion. Communications Biology, Mar 2020. URL: https://doi.org/10.1038/s42003-020-0824-5, doi:10.1038/s42003-020-0824-5. This article has 102 citations and is from a peer-reviewed journal.

7. (otani2014characterizationofp pages 1-2): Hiroshi Otani, Young-Eun Lee, Israël Casabon, and Lindsay D. Eltis. Characterization of p -hydroxycinnamate catabolism in a soil actinobacterium. Journal of Bacteriology, 196:4293-4303, Dec 2014. URL: https://doi.org/10.1128/jb.02247-14, doi:10.1128/jb.02247-14. This article has 74 citations and is from a peer-reviewed journal.

8. (li2015structuralbasisfor pages 7-8): Zhi Li and Satish K. Nair. Structural basis for specificity and flexibility in a plant 4-coumarate:coa ligase. Structure, 23:2032-2042, Nov 2015. URL: https://doi.org/10.1016/j.str.2015.08.012, doi:10.1016/j.str.2015.08.012. This article has 93 citations and is from a domain leading peer-reviewed journal.

9. (li2015structuralbasisfor pages 4-5): Zhi Li and Satish K. Nair. Structural basis for specificity and flexibility in a plant 4-coumarate:coa ligase. Structure, 23:2032-2042, Nov 2015. URL: https://doi.org/10.1016/j.str.2015.08.012, doi:10.1016/j.str.2015.08.012. This article has 93 citations and is from a domain leading peer-reviewed journal.

10. (li2015structuralbasisfor pages 6-7): Zhi Li and Satish K. Nair. Structural basis for specificity and flexibility in a plant 4-coumarate:coa ligase. Structure, 23:2032-2042, Nov 2015. URL: https://doi.org/10.1016/j.str.2015.08.012, doi:10.1016/j.str.2015.08.012. This article has 93 citations and is from a domain leading peer-reviewed journal.

11. (li2015structuralbasisfor pages 5-6): Zhi Li and Satish K. Nair. Structural basis for specificity and flexibility in a plant 4-coumarate:coa ligase. Structure, 23:2032-2042, Nov 2015. URL: https://doi.org/10.1016/j.str.2015.08.012, doi:10.1016/j.str.2015.08.012. This article has 93 citations and is from a domain leading peer-reviewed journal.

12. (khurana2010genomescaleprediction pages 6-8): Pankaj Khurana, Rajesh S Gokhale, and Debasisa Mohanty. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles. BMC Bioinformatics, 11:57-57, Jan 2010. URL: https://doi.org/10.1186/1471-2105-11-57, doi:10.1186/1471-2105-11-57. This article has 45 citations and is from a peer-reviewed journal.

13. (khurana2010genomescaleprediction pages 4-6): Pankaj Khurana, Rajesh S Gokhale, and Debasisa Mohanty. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles. BMC Bioinformatics, 11:57-57, Jan 2010. URL: https://doi.org/10.1186/1471-2105-11-57, doi:10.1186/1471-2105-11-57. This article has 45 citations and is from a peer-reviewed journal.

14. (khurana2010genomescaleprediction pages 3-4): Pankaj Khurana, Rajesh S Gokhale, and Debasisa Mohanty. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles. BMC Bioinformatics, 11:57-57, Jan 2010. URL: https://doi.org/10.1186/1471-2105-11-57, doi:10.1186/1471-2105-11-57. This article has 45 citations and is from a peer-reviewed journal.

15. (khurana2010genomescaleprediction pages 10-13): Pankaj Khurana, Rajesh S Gokhale, and Debasisa Mohanty. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles. BMC Bioinformatics, 11:57-57, Jan 2010. URL: https://doi.org/10.1186/1471-2105-11-57, doi:10.1186/1471-2105-11-57. This article has 45 citations and is from a peer-reviewed journal.

16. (khurana2010genomescaleprediction pages 14-16): Pankaj Khurana, Rajesh S Gokhale, and Debasisa Mohanty. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles. BMC Bioinformatics, 11:57-57, Jan 2010. URL: https://doi.org/10.1186/1471-2105-11-57, doi:10.1186/1471-2105-11-57. This article has 45 citations and is from a peer-reviewed journal.

17. (priyadarshan2018fattyacylampligases pages 7-9): Kinatukara Priyadarshan and Rajan Sankaranarayanan. Fatty acyl-amp ligases as mechanistic variants of anl superfamily and molecular determinants dictating substrate specificities. Journal of the Indian Institute of Science, 98:261-272, Jul 2018. URL: https://doi.org/10.1007/s41745-018-0084-2, doi:10.1007/s41745-018-0084-2. This article has 13 citations.

18. (priyadarshan2018fattyacylampligases pages 2-4): Kinatukara Priyadarshan and Rajan Sankaranarayanan. Fatty acyl-amp ligases as mechanistic variants of anl superfamily and molecular determinants dictating substrate specificities. Journal of the Indian Institute of Science, 98:261-272, Jul 2018. URL: https://doi.org/10.1007/s41745-018-0084-2, doi:10.1007/s41745-018-0084-2. This article has 13 citations.

19. (darrigo2019analysisofpseudomonas pages 14-16): Isotta D'Arrigo, João G. R. Cardoso, Maja Rennig, Nikolaus Sonnenschein, Markus J. Herrgård, and Katherine S. Long. Analysis of pseudomonas putida growth on non‐trivial carbon sources using transcriptomics and genome‐scale modelling. Environmental Microbiology Reports, 11:87–97, Nov 2019. URL: https://doi.org/10.1111/1758-2229.12704, doi:10.1111/1758-2229.12704. This article has 38 citations and is from a peer-reviewed journal.

20. (li2015structuralbasisfor pages 1-3): Zhi Li and Satish K. Nair. Structural basis for specificity and flexibility in a plant 4-coumarate:coa ligase. Structure, 23:2032-2042, Nov 2015. URL: https://doi.org/10.1016/j.str.2015.08.012, doi:10.1016/j.str.2015.08.012. This article has 93 citations and is from a domain leading peer-reviewed journal.

21. (li2015structuralbasisfor pages 8-10): Zhi Li and Satish K. Nair. Structural basis for specificity and flexibility in a plant 4-coumarate:coa ligase. Structure, 23:2032-2042, Nov 2015. URL: https://doi.org/10.1016/j.str.2015.08.012, doi:10.1016/j.str.2015.08.012. This article has 93 citations and is from a domain leading peer-reviewed journal.

22. (khurana2010genomescaleprediction pages 8-10): Pankaj Khurana, Rajesh S Gokhale, and Debasisa Mohanty. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles. BMC Bioinformatics, 11:57-57, Jan 2010. URL: https://doi.org/10.1186/1471-2105-11-57, doi:10.1186/1471-2105-11-57. This article has 45 citations and is from a peer-reviewed journal.

23. (khurana2010genomescaleprediction pages 13-14): Pankaj Khurana, Rajesh S Gokhale, and Debasisa Mohanty. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles. BMC Bioinformatics, 11:57-57, Jan 2010. URL: https://doi.org/10.1186/1471-2105-11-57, doi:10.1186/1471-2105-11-57. This article has 45 citations and is from a peer-reviewed journal.

Artifacts

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

Final Report: Evaluation of GO:0031956 (medium-chain fatty acid-CoA ligase activity) for P. putida fcs (Q88HK0)

Executive Judgment

Verdict: REFUTED — Wrong subfamily (within-superfamily mis-placement, failure mode 3)

The seed hypothesis that P. putida KT2440 fcs (Q88HK0) has medium-chain fatty acid-CoA ligase activity (GO:0031956) is refuted with high confidence. The protein is a well-characterized feruloyl-CoA synthetase (EC 6.2.1.34, GO:0050563) that catalyzes the ATP-dependent activation of trans-ferulate and other hydroxycinnamic acids to their CoA thioesters as the first step in aromatic compound catabolism. Both feruloyl-CoA synthetase and medium-chain fatty acid-CoA ligase belong to the AMP-binding/adenylate-forming superfamily (PF00501/IPR000873), but they act on fundamentally different substrates: aromatic phenylpropanoids versus aliphatic C6-C12 fatty acids. The TreeGrafter annotation therefore represents a classic within-superfamily mis-placement: the protein was grafted onto a structurally related but functionally distinct branch of the PANTHER ACYL-COA SYNTHETASE family tree (PTHR43201), and a GO term for fatty acid-CoA ligase activity was incorrectly propagated from an ancestral node.

Most decisive evidence: Direct enzymatic assay of recombinant Fcs protein expressed in E. coli demonstrated feruloyl-CoA synthetase activity, and a chromosomal fcs knockout mutant (KT2440 fcsOmegaKm) completely lost the ability to grow on ferulic acid (PMID: 12764569).

Summary

Pseudomonas putida KT2440 fcs (UniProt Q88HK0) was annotated by TreeGrafter/PANTHER with the GO term GO:0031956 (medium-chain fatty acid-CoA ligase activity) via phylogenetic propagation (IEA, GO_REF:0000118). This investigation tested that annotation against three characteristic failure modes of automated phylogenetic function inference: (1) granularity mismatch (family vs. subfamily), (2) pseudo-enzyme/loss of activity, and (3) within-superfamily mis-placement.

The evidence overwhelmingly identifies within-superfamily mis-placement as the failure mode. The fcs gene product has been directly characterized as a feruloyl-CoA synthetase (EC 6.2.1.34) through heterologous expression in E. coli, enzyme activity assays, and gene knockout experiments in P. putida KT2440 (PMID: 12764569). The protein catalyzes the ATP-dependent activation of ferulic acid (and related hydroxycinnamic acids including caffeic acid and p-coumaric acid) to feruloyl-CoA, the first committed step in the lignin-derived aromatic compound catabolic pathway leading to vanillin and then vanillic acid. This function is captured by GO:0050563 (trans-feruloyl-CoA synthase activity), which is a sibling term to GO:0031956 under the shared parent GO:0016405 (CoA-ligase activity). The two terms diverge at the level of substrate specificity: medium-chain fatty acids (aliphatic C6-C12) versus hydroxycinnamic acids (aromatic phenylpropanoids).

Active-site analysis confirmed that the protein retains the canonical AMP-binding P-loop motif (TSGSTKLPK at positions 228-236) with an intact catalytic lysine (K234), definitively ruling out pseudo-enzyme status. BLAST analysis against SwissProt showed that 4-coumarate-CoA ligase (an aromatic acid-CoA ligase acting on the same hydroxycinnamic acid substrate class) scored higher than FadK (a medium-chain fatty acid-CoA ligase), and CDD classification placed Q88HK0 in the cd05921 subfamily (FCS = feruloyl-CoA synthetase), not in a fatty acid-CoA ligase subfamily. These computational findings are consistent with the biochemical evidence but would not on their own be sufficient — it is the direct assay and knockout data that provide definitive resolution.

Key Findings

Finding 1: fcs is a feruloyl-CoA synthetase, not a medium-chain fatty acid-CoA ligase

The identity of fcs as a feruloyl-CoA synthetase rests on multiple lines of direct experimental evidence from P. putida KT2440, the exact strain specified in the annotation. Jimenez et al. (2002) cloned fcs from P. putida KT2440 and expressed it in E. coli, demonstrating that recombinant strains harboring fcs and the downstream ech gene converted ferulic acid to vanillin — a reaction that proceeds via the feruloyl-CoA intermediate (PMID: 12764569). As stated in that paper: "To confirm the physiological function of these structural genes as feruloyl-CoA synthetase (Fcs), enoyl-CoA hydratase/aldolase (Ech), and vanillin dehydrogenase (Vdh), respectively, they were cloned and expressed in Escherichia coli." The same study constructed a chromosomal fcs knockout (fcsOmegaKm) and showed that the mutant strain was completely unable to grow on ferulic acid as sole carbon source: "The essential involvement of fcs, ech and vdh in the catabolism of ferulic acid in P. putida KT2440 was proven by separately inactivating each gene by insertion of Omega-elements. The corresponding mutant strains KT2440 fcsOmegaKm, KT2440 echOmegaKm, and KT2440 vdhOmegaKm were not able to grow on ferulic acid." This confirms the essential and non-redundant role of fcs in ferulic acid catabolism.

Subsequent studies reinforced this assignment. Graf and Altenbuchner (2014) engineered P. putida KT2440 for industrial vanillin production by overexpressing fcs (feruloyl-CoA synthetase) and ech under a strong tac promoter, treating the feruloyl-CoA synthetase identity as established fact (PMID: 24136472): "The bioconversion was optimized by enhanced chromosomal expression of the structural genes for feruloyl-CoA synthetase (fcs) and enoyl-CoA hydratase/aldolase (ech) by introduction of the strong tac promoter system." Molecular studies in the related species P. putida F1 demonstrated that a mutant lacking fcs was unable to exhibit chemotaxis toward p-coumaric, caffeic, or ferulic acids (PMID: 28954643): "a mutant lacking the gene encoding feruloyl-CoA synthetase (Fcs), which catalyzes the first step in hydroxycinnamic acid degradation, was unable to respond chemotactically toward p-coumaric, caffeic, or ferulic acids."

The substrate specificity of fcs is firmly in the hydroxycinnamic acid class (ferulic acid = 4-hydroxy-3-methoxycinnamic acid, MW 194), which is structurally unrelated to medium-chain fatty acids (C6-C12 saturated aliphatic carboxylic acids). The two substrate classes share only the terminal carboxylate that undergoes adenylation, but differ in the aromatic versus aliphatic backbone, the presence of hydroxyl/methoxy ring substituents, and the alpha,beta-unsaturation of the propanoid side chain.

Finding 2: Active-site residues are intact — not a pseudo-enzyme

Analysis of the Q88HK0 amino acid sequence identified the conserved AMP-binding P-loop motif TSGSTKLPK at positions 228-236. This matches the canonical signature [LIVMFY]xx[STG][STAG]G[ST][TSE][GS]x[PASLIVM]K of active AMP-forming CoA ligases. The catalytic lysine at position 234 (K234) is present and correctly positioned. Additional conserved motifs were identified: PKG at position 323 and the catalytic KLFFF motif at position 348.

Comparison with characterized feruloyl-CoA synthetases (S5M744 from Sphingobium sp. SYK-6 and Q9EY88 from Amycolatopsis sp.) showed conserved P-loop sequences (TSGSTGRPK), confirming that Q88HK0 belongs to the same catalytically active clade. The protein is definitively not a pseudo-enzyme: it retains all the catalytic machinery required for AMP-forming CoA-ligase activity, and direct assays confirm it is enzymatically active. This rules out failure mode 2 (pseudo-enzyme / loss of activity).

P-loop alignment across characterized homologs:

Q88HK0  (P. putida fcs, QUERY):             ...AFAATGPDTIAKFLFTSGSTKLPKAVITTQRMLCA...
S5M744  (Streptomyces FCS, EC 6.2.1.34):    ...PDLPVGLDDVCLLMYTSGSTGRPKGAMLTHGNLTW...
Q9EY88  (Amycolatopsis FCS, EC 6.2.1.34):   ...PDLPVGLDDVCLLMYTSGSTGRPKGAMLTHGNLTW...
                                          *****  **

Finding 3: BLAST-based placement confirms aromatic acid-CoA ligase subfamily affinity

BLAST of Q88HK0 against SwissProt returned hits across the AMP-binding superfamily at 20-25% sequence identity, reflecting the deep divergence within this superfamily. Critically, the sequence similarity ranking favored aromatic acid-CoA ligases over fatty acid-CoA ligases:

4-coumarate-CoA ligase (P41636) acts on the same hydroxycinnamic acid substrate class (coumaric acid, caffeic acid, ferulic acid) as feruloyl-CoA synthetase, and it outscored the medium-chain fatty acid-CoA ligase FadK. This sequence-level signal, combined with the CDD classification of Q88HK0 into subfamily cd05921 (FCS = feruloyl-CoA synthetase), confirms that the protein belongs to the aromatic/hydroxycinnamic acid-CoA ligase branch of the superfamily rather than the fatty acid-CoA ligase branch.

At the ~20-25% identity level, BLAST scores alone would be insufficient for confident function assignment — the score difference between 4CL (77.8 bits) and FadK (76.3 bits) is small. However, the CDD subfamily classification (cd05921) is more informative at this divergence level because it incorporates curated position-specific scoring matrices that capture subfamily-diagnostic residue patterns.

Independent Family/Function Assignment

Mechanistic Model

The fcs gene product operates within a well-characterized three-gene catabolic operon in P. putida KT2440 that degrades ferulic acid (and related hydroxycinnamic acids from lignin breakdown) to central metabolites:

Ferulic acid (4-hydroxy-3-methoxycinnamic acid)
    |
    v  fcs (feruloyl-CoA synthetase, EC 6.2.1.34)
    |  ATP + CoA --> AMP + PPi
    |
Feruloyl-CoA
    |
    v  ech (enoyl-CoA hydratase/aldolase)
    |  Retro-aldol cleavage
    |
Vanillin + Acetyl-CoA
    |
    v  vdh (vanillin dehydrogenase)
    |  NAD+ --> NADH
    |
Vanillic acid
    |
    v  (further catabolism via demethylation --> protocatechuate --> beta-ketoadipate)
    |
TCA cycle

The fcs-catalyzed reaction is the first committed step: ATP-dependent adenylation of ferulic acid's carboxylate to form a feruloyl-AMP intermediate, followed by thioesterification with CoA to yield feruloyl-CoA. This is the same catalytic mechanism used by all AMP-binding superfamily CoA ligases (including fatty acid-CoA ligases), but the substrate-binding pocket of fcs is configured for planar aromatic hydroxycinnamic acids rather than flexible aliphatic fatty acyl chains.

The TreeGrafter mis-assignment arose because the AMP-binding superfamily (Pfam PF00501 / InterPro IPR000873) is extremely large and functionally diverse, encompassing:

• Fatty acid-CoA ligases (long-chain, medium-chain, short-chain)
• 4-coumarate-CoA ligases
• Feruloyl-CoA synthetases
• Benzoate-CoA ligases
• Amino acid adenylating enzymes (NRPS adenylation domains)
• Firefly luciferases
• Many other activities

All share the same core fold and catalytic mechanism (adenylation followed by thioesterification), but differ in their substrate specificity, which is determined by the substrate-binding pocket residues. Phylogenetic grafting onto a tree that does not resolve substrate-specificity-determining residues can easily place an aromatic acid-CoA ligase on a branch annotated with a fatty acid-CoA ligase function.

GO Term Hierarchy Relationship

GO:0003674 (molecular_function)
  └── GO:0003824 (catalytic activity)
      └── GO:0016874 (ligase activity)
  └── GO:0016877 (ligase activity, forming C-S bonds)
      └── GO:0016878 (acid-thiol ligase activity)
  └── GO:0016405 (CoA-ligase activity)
      ├── GO:0031956 (medium-chain fatty acid-CoA ligase) <-- SEED TERM (WRONG)
      │   └── via GO:0015645 --> GO:0140657 (fatty acid CoA-ligase)
      └── GO:0050563 (trans-feruloyl-CoA synthase)        <-- CORRECT TERM

GO:0031956 and GO:0050563 are sibling terms under GO:0016405 (CoA-ligase activity). They share the catalytic mechanism (CoA-thioester bond formation using ATP) but diverge at substrate specificity. The seed term is not merely too general or too specific — it names an entirely different substrate class (aliphatic fatty acids vs. aromatic hydroxycinnamic acids).

Evidence Matrix

Active-Site / Placement Analysis

Active-Site Residue Conservation

The AMP-binding superfamily catalytic machinery was analyzed by searching the Q88HK0 sequence for conserved motifs:

Conclusion: All catalytic residues required for AMP-forming CoA ligase activity are present and correctly positioned. The protein is not a pseudo-enzyme (failure mode 2 is definitively ruled out).

Subfamily Placement Analysis

The TreeGrafter mis-annotation stems from incorrect subfamily placement: Q88HK0 was grafted into PTHR43201:SF32 (2-succinylbenzoate-CoA ligase subfamily), and GO:0031956 (medium-chain fatty acid-CoA ligase) was propagated from an ancestral node. The protein actually belongs to the feruloyl-CoA synthetase subfamily, which CDD correctly identifies as cd05921.

GO Curation Implications

Recommended action: REPLACE with sibling term

The evidence from PMID: 12764569 (direct assay + knockout) is sufficient to support upgrading GO:0050563 from IEA to IDA (Inferred from Direct Assay) or IMP (Inferred from Mutant Phenotype) evidence codes, pending curator review. The protein's function has been characterized in three independent studies spanning nearly two decades.

Failure Mode Analysis Summary

Conflicts, Knowledge Gaps, and Discriminating Tests

Conflicts Identified

1. TreeGrafter vs. UniProt-EC annotation conflict: GO:0031956 (TreeGrafter, IEA) directly contradicts GO:0050563 (UniProt-EC, IEA) on the same protein. The UniProt-EC annotation is correct. Both annotations currently coexist on Q88HK0, creating a contradictory annotation state.

2. PANTHER subfamily mis-assignment: Q88HK0 is placed in PTHR43201:SF32 (2-succinylbenzoate-CoA ligase), which is itself incorrect — the protein is neither a 2-succinylbenzoate-CoA ligase nor a medium-chain fatty acid-CoA ligase. This suggests the PANTHER tree topology in this region may not adequately resolve bacterial aromatic acid-CoA ligases.

3. Cascading biological process error: The TreeGrafter also propagated GO:0006631 (fatty acid metabolic process), which is incorrect. fcs functions in aromatic compound catabolism (hydroxycinnamic acid degradation to vanillic acid via vanillin), not in fatty acid metabolism.

Knowledge Gaps

• Substrate range quantification: While fcs is confirmed to act on ferulic acid, caffeic acid, and p-coumaric acid, detailed Km/Vmax kinetic parameters for each substrate with the P. putida KT2440 enzyme specifically have not been published. Kinetics are available for Amycolatopsis and Streptomyces homologs.

• Structural basis of specificity: No crystal structure of bacterial Fcs enzymes is available. Structural data would definitively show the aromatic binding pocket that distinguishes Fcs from fatty acid-CoA ligases.

• Potential for residual fatty acid activity: It is theoretically possible that fcs has some low-level activity on medium-chain fatty acids, as many AMP-binding superfamily enzymes show some substrate promiscuity. However, no evidence supports this, and the gene knockout phenotype is specific to hydroxycinnamic acid catabolism.

Discriminating Tests

1. Most efficient experimental resolution: A substrate competition assay comparing purified Fcs activity on ferulic acid versus octanoic acid (a C8 medium-chain fatty acid) would unambiguously discriminate between the two GO terms. Based on all available evidence, Fcs would show robust activity on ferulate and negligible activity on octanoate.

2. Computational follow-up (partially completed): BLAST analysis and CDD classification were performed and support the aromatic acid-CoA ligase assignment. A more comprehensive phylogenetic analysis with multiple characterized members of each subfamily would provide stronger computational evidence.

3. Structure prediction comparison: AlphaFold structure comparison with a known medium-chain fatty acid-CoA ligase (e.g., Mycobacterium FadD6) would reveal divergent substrate-binding pocket geometries.

4. TreeGrafter diagnostic: Examining the PANTHER tree to identify why fcs was grafted onto the fatty acid-CoA ligase branch could reveal a systematic issue affecting other aromatic acid-CoA ligases annotated by the same pipeline.

Evidence Base (Literature)

Primary References

1. Jimenez JI, Minambres B, Garcia JL, Diaz E (2002). Functional analyses of genes involved in the metabolism of ferulic acid in Pseudomonas putida KT2440. J Bacteriol. PMID: 12764569

2. The foundational study directly characterizing fcs in P. putida KT2440. Cloned fcs, ech, and vdh; expressed in E. coli to confirm enzymatic activities. Demonstrated that recombinant E. coli expressing fcs and ech converted ferulic acid to vanillin, confirming feruloyl-CoA synthetase activity. Chromosomal fcs knockout (fcsOmegaKm) abolished growth on ferulic acid. This constitutes direct biochemical + genetic evidence in the exact organism and strain.

3. Graf N, Altenbuchner J (2014). Genetic engineering of Pseudomonas putida KT2440 for rapid and high-yield production of vanillin from ferulic acid. Appl Microbiol Biotechnol. PMID: 24136472

4. Applied biotechnology study that engineered P. putida KT2440 for industrial vanillin production by overexpressing fcs and ech under a strong tac promoter. Treats the identity of fcs as feruloyl-CoA synthetase as established, providing independent confirmation from a different laboratory.

5. Ramirez-Morales JE et al. (2017). Pseudomonas putida F1 uses energy taxis to sense hydroxycinnamic acids. Appl Environ Microbiol. PMID: 28954643

6. Demonstrated that a P. putida F1 mutant lacking fcs was unable to respond chemotactically toward p-coumaric, caffeic, or ferulic acids. Confirms that fcs catalyzes the first step in hydroxycinnamic acid degradation in a second P. putida strain, extending functional characterization beyond KT2440.

Limitations and Knowledge Gaps

1. Substrate promiscuity not tested: This analysis did not experimentally test whether Fcs has any residual activity on medium-chain fatty acids. While no published evidence supports such activity, complete exclusion would require a kinetic study with purified enzyme and fatty acid substrates.

2. BLAST limitations at low identity: At 20-25% sequence identity between Q88HK0 and characterized SwissProt entries, BLAST-based placement is supportive but not definitive. The CDD classification and direct biochemical evidence carry more weight.

3. PANTHER tree not examined: The specific tree topology and grafting point that led to the mis-annotation were not examined. This limits understanding of whether the error is systematic or specific to this protein.

4. No crystal structure available: The active-site analysis was performed at the sequence motif level. A structural comparison of substrate-binding pockets would provide stronger evidence for substrate discrimination.

Proposed Follow-up Experiments/Actions

1. Immediate curation action (high priority): Replace GO:0031956 (IEA) with GO:0050563 (IDA/IMP) citing PMID: 12764569. Remove GO:0006631 (fatty acid metabolic process). This is actionable now with high confidence.

2. TreeGrafter pipeline audit (medium priority): Search for other AMP-binding superfamily members annotated by TreeGrafter with GO:0031956 or related fatty acid-CoA ligase terms, and check whether any are actually aromatic acid-CoA ligases. The mis-placement pattern identified here may be systematic.

3. Substrate range characterization (lower priority): If substrate promiscuity is a concern for annotation purposes, purified recombinant fcs could be tested for activity on a panel of fatty acids (C6-C12) alongside ferulic acid. Based on all available evidence, fatty acid activity would be expected to be negligible.

4. Structural analysis (lower priority): AlphaFold structure comparison with characterized fatty acid-CoA ligases versus plant 4-coumarate-CoA ligases could define the structural basis for substrate discrimination and inform improved subfamily classification.

Artifacts

• OpenScientist final report
• OpenScientist final report