RXNRECer Enables Fine-grained Enzymatic Function Annotation through Active Learning and Protein Language Models
Authors: Zhenkun Shi, Jun Zhu, Dehang Wang, BoYu Chen, Qianqian Yuan, Zhitao Mao, Fan Wei, Weining Wu, Xiaoping Liao, Hongwu Ma
Organizations: Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, 300308, Tianjin, China. · College of Computer Science and Technology, Harbin Engineering University, Street, Harbin, 10587, State, China.
A key challenge in enzyme annotation is identifying the biochemical reactions catalyzed by proteins. Most existing methods rely on Enzyme Commission (EC) numbers as intermediaries: they first predict an EC number and then retrieve the associated reactions. This indirect strategy introduces ambiguity due to the complex many-to-many mappings among proteins, EC numbers, and reactions, and is further complicated by frequent updates to EC numbers and inconsistencies across databases. To address these challenges, we present RXNRECer, a transformer-based ensemble framework that directly predicts enzyme-catalyzed reactions without relying on EC numbers. It integrates protein language modeling and active learning to capture both high-level sequence semantics and fine-grained transformation patterns. Evaluations on curated cross-validation and temporal test sets demonstrate consistent improvements over six EC-based baselines, with gains of 16.54% in F1 score and 15.43% in accuracy. Beyond accuracy gains, the framework offers clear advantages for downstream applications, including scalable proteome-wide reaction annotation, enhanced specificity in refining generic reaction schemas, systematic annotation of previously uncurated proteins, and reliable identification of enzyme promiscuity. By incorporating large language models, it also provides interpretable rationales for predictions. These capabilities make RXNRECer a robust and versatile solution for EC-free, fine-grained enzyme function prediction, with potential applications across multiple areas of enzyme research and industrial applications.
Figures & tables
Figure 1: Schematic overview of the RXNRECer framework and its active learning workflow. (I) RXNRECer Framework. (a) Sequence-based reaction prediction (RXNRECer-S1). Protein sequences are embedded using a pre-trained language model. The resulting multi-level embeddings are passed through GRU and Transformer layers, followed by a fully connected (FC) layer to produce reaction scores. This EC-independent prediction module is fine-tuned via a reaction-level classification task within an active learning framework. (b) Ensemble integration (RXNRECer-S2). Predictions from EC-based and PLM similarity-based methods are combined with the direct predictor using a voting scheme to improve robustness across diverse reaction types. (c) Interpretation and ranking (RXNRECer-S3). A general-purpose language model (GLM) is prompted to re-rank the candidate reactions and generate concise, biologically meaningful justifications. (II) Active Learning Strategy. (a) Initial model training. The model is initially trained on a small labeled set DLi . (b) Validation. Performance is evaluated on a validation set DV sampled from the unlabeled pool DU . (c–e) Sample selection. An acquisition function S(x) identifies the most informative samples for labeling. (f–g) Fine-tuning. The selected samples are incorporated into the training set DL to refine the model. (h–i) Validation sampling. Additional samples are selected to monitor performance after each fine-tuning iteration.
Figure 2: Comparative evaluation of enzyme reaction prediction methods on two datasets. (I) Results on the 10-fold cross-validation dataset ( ds_rcv ). (II) Results on the independent temporal test dataset ( ds_rcp ). For each panel: (a) Performance of EC-based methods (ECRECer, CatFam, CLEAN, PRIAM, DeepEC, and MSA-via-EC). (b) Performance of PLM-based similarity methods (T5-cosine, T5-euclidean, ESM-cosine, ESM-euclidean, Unirep-cosine, Unirep-euclidean), the end-to-end fine-tuned RXNRECer-S1, and MSA-via-RXN. Results on ds_rcv are averaged over 10 folds with error bars showing standard deviations; results on ds_rcp are from a single evaluation without error bars.
Figure 3: Performance evaluation of ensemble strategies on ds_rcp from three perspectives: (a) macro-averaged Precision (mPrecision), (b) macro-averaged Recall (mRecall), and (c) run time with respect to the number of predictors. Boxplots summarize the distribution of results across different ensemble settings, with individual points representing Recall-boosted (red) and Majority-vote (blue) strategies, and orange dots indicating the average values. The results highlight the trade-offs between precision, recall, and computational cost when adopting different ensemble strategies.
Figure 4: Case study on large-scale reaction prediction for the Fusarium venenatum proteome (FS12832). (I) Coverage Comparison. Heatmap showing the proportion of FS12832 proteins assigned to different functional categories by various methods. (II) Structure Consistency Evaluation. (a) Kernel Density Estimation (KDE) plot comparing TM-score distributions between RXNRECer and CLEAN. (b) Bin-wise TM-score comparison and Wilcoxon test showing RXNRECer achieves higher structural consistency. (III) Non-Enzyme Misannotated as Enzyme. (a) TM-score heatmap showing low structural similarity between CLEAN-predicted proteins and canonical enzymes. (b) Structural alignment confirming CLEAN incorrectly assigned EC 2.4.99.19 to non-enzymes. (IV) RXNRECer vs. CLEAN. (a) TM-score heatmap comparing predictions from CLEAN and RXNRECer. (b) Structural alignment showing RXNRECer correctly identifies enzyme functions, while CLEAN fails.
Figure 5: (I) Representative cases illustrating RXNRECer’s capability to resolve general reaction schemas into substrate-specific enzymatic transformations. (a) refines the generic aldehyde oxidation (RHEA:16185) to a specific reaction acting on acetaldehyde (RHEA:25294); (b) resolves a primary amine oxidation template (RHEA:16153) into a specific transformation involving 2-phenylethylamine (RHEA:22265). (II) Representative cases illustrating RXNRECer’s capability to predict reaction-level functions in reviewed proteins lacking curated catalytic annotations, exemplified by P31852 (a) and P73408 (b).
Figure 6: Representative examples illustrating RXNRECer’s ability to uncover enzyme promiscuity through reaction-level predictions. (I) Predicted CoA-ligase reactions catalyzed by two homologous enzymes (A0A2H5AIX5 and A0A2H5AIY4), both originally annotated with a single 4-coumarate: CoA ligase activity (RHEA:19641). RXNRECer further identifies multiple substrate-specific reactions involving (a) 4-coumarate, (b) ferulate, (c) caffeate, and (d) cinnamate. (II) Schematic of the two-step reaction mechanism for (E)-4-coumarate: CoA ligase. The process proceeds via an AMP-linked intermediate: (a) adenylation of the substrate by ATP to form 4-coumaroyl-AMP, and (b) ligation with CoA to yield 4-coumaroyl-CoA. (III) Structural alignment of A0A2H5AIX5 and A0A2H5AIY4 with O24146 (PDB ID: 5bst), a known rice 4CL1 enzyme. The high structural similarity supports the predicted catalytic versatility across different phenylpropanoid substrates.
Existing reaction models primarily learn molecular transformations, whereas enzy- matic reactions depend jointly on molecular structure and catalytic function. We formulate this problem as learning an enzymatic reaction space linking reactants, products, and Enzyme Commission (EC) annotations. To characterize this space, we introduce VenusRX-Bench, a unified benchmark for forward reaction prediction, single-step retrosynthesis, and EC-number prediction. VenusRX-Bench integrates reactions from multiple biochemical databases with standardized curation, leakage- controlled splits, and consistent evaluation. Benchmarking representative chemical and enzymatic models reveals a clear chemical-to-enzymatic domain gap, driven by limited domain data, catalytic-context dependency, and the difficulty of modeling large biomolecular structures. To bridge this gap, we develop VenusRX, a unified T5-style sequence-to-sequence model for enzymatic reactions. VenusRX jointly learns forward prediction, ret- rosynthesis, and reaction reconstruction, with two-stage training on millions of template-expanded reactions followed by real biochemical reactions. In addition, optional EC conditioning incorporates catalytic context, while Molecule Library- Constrained Decoding improves the generation of complex biomolecules. Across benchmark tasks and challenging generalization splits, VenusRX achieves the best or competitive performance on most evaluated settings over representative chem- ical and enzymatic baselines. Moreover, EC information consistently improves reaction prediction, while learned reaction representations support accurate EC prediction, revealing a bidirectional relationship between reaction structure and catalytic function. Together, VenusRX-Bench and VenusRX provide a unified framework for elucidating and modeling enzymatic reaction space
Yutong Hu, Tianming Huang, Yanbo Zhao +6
Shanghai Jiao Tong University · Shanghai Artificial Intelligence Laboratory · The University of Sydney
Enzyme-reaction retrieval is a fundamental problem in computational biology, underpinning enzyme characterization, reaction mechanism elucidation, and the rational design of metabolic pathways and biocatalysts. As a bidirectional task, it entails both enzyme-to-reaction and reaction-to-enzyme mapping. However, existing approaches suffer from poor generalization across tasks and distributions, with performance highly sensitive to dataset splits and substantial asymmetry between retrieval directions. To address these challenges, we present TIGER, a Text-Informed Generalized Enzyme-Reaction Retrieval framework that leverages protein-to-text generation models to distill textual semantic knowledge from enzyme sequences, providing a generalized representation that bridges enzymes and biochemical reactions. To ensure the quality and reliability of textual semantics, we design a Dynamic Gating Network that adaptively fuses text-derived knowledge with sequence features, enabling more consistent and informative enzyme representations, while a Structure-Shared Feature Projector aligns enzyme and reaction representations within a unified latent space. Extensive experiments demonstrate that, under bidirectional retrieval supervision, TIGER significantly outperforms state-of-the-art baselines across diverse distributions and exhibits strong robustness and transferability across tasks.
Yuhang Zhang, Keyan Ding, Peilin Chen +5
University of Science and Technology of China · 2City University of Hong Kong · 3Zhejiang University
Enzyme function prediction is a hierarchical, knowledge-intensive form of protein function classification. Existing benchmarks expose an anomaly: general LLMs often get the coarse first level right, yet once asked for a complete EC number their accuracy at levels two through four drops to almost zero, while specialized models and tools stay usable. We propose EC-Reason-Bench, a training-free, diagnostic evaluation protocol built to answer two questions: why general LLMs score close to nothing on EC number prediction, and how much of that loss can be recovered without updating a single weight. We break enzyme classification ability into four orthogonal levers that can each be measured on their own: output structure, external knowledge, reasoning structure, and reasoning robustness. We test each lever with an inference-time method against a shared zero-shot baseline reproducing previously reported near-zero performance. Experiments with several strong reasoning LLMs yield four main findings. First, external knowledge is decisive and must precede reasoning: uniformly low closed-book performance rises sharply with open-book access, narrowing model gaps. Second, in closed-book settings, whether cascading and chain-of-thought help or hurt depends on a model's tendency to abstain. Third, once evidence is available the aggregate score of the best LLM setting is indistinguishable from simply voting the EC numbers of the nearest retrieved neighbors; that tie is an artifact of averaging, and it hides a large gain on adversarial evidence set against an equally large loss on multi-functional enzymes. Reasoning over evidence therefore acts as an arbiter of conflicting neighbors rather than as a source of knowledge, and no single-number leaderboard can see it. Fourth, accuracy obeys a law of homology availability.
Linyu Li, Zhi Jin, Yichi Zhang +6
School of Computer Science, Peking University, Beijing, China · College of Computer Science and Technology, Zhejiang University, Hangzhou, China · School of Computer Science, Tibet University, Lhasa, China