The scarcity of observed biodegradation reactions constrains the training of product and pathway predictors. We present enviGraph, a framework that generates and curates synthetic reaction data, uses these data to pretrain a graph-to-sequence predictor, and fine-tunes the predictor on observed environmental reactions for product and recursive pathway prediction. Applying inherited environmental transformation rules to structures from PubChem, NORMAN SusDat and EPA DSSTox, followed by curation, produces 776,851 canonical parent-product pairs for pretraining. We compare this dataset with three United States Patent and Trademark Office (USPTO) reaction datasets: USPTO-50K, USPTO-MIT and USPTO-Stereo. We evaluate enviGraph and the sequence-based enviFormer using ten folds split by pathway on the EAWAG Biocatalysis/Biodegradation Database (BBD), regulatory pesticide biodegradation data (Soil) and activated-sludge biodegradation data (Sludge). For both predictors, synthetic pretraining gives the highest BBD mean precision-recall area under the curve (PR-AUC) among the evaluated pretraining datasets for single-generation and multi-generation prediction. enviGraph achieves BBD single-generation PR-AUC of 0.202±0.050 with synthetic pretraining, compared with 0.169±0.065 for the retrospectively selected strongest USPTO condition, USPTO-Stereo (mean ± sample standard deviation across folds). Its exact BBD top-1 and top-5 accuracies with synthetic pretraining are 0.355 and 0.594. The strongest USPTO conditions give higher PR-AUC means than synthetic pretraining in most Soil and Sludge comparisons. enviGraph has higher Soil multi-generation PR-AUC means than enviFormer under each of the four main pretraining datasets, and with synthetic pretraining it also leads enviFormer on BBD single-generation PR-AUC and exact top-1 and top-5 accuracy.
Figures & tables
Figure 1: The integrated enviGraph framework.
Dataset
Compounds
Reactions
Rules
Pathways
BBD
1,399
1,480
290
219
Soil
2,609
2,447
80
317
Sludge
1,070
521
0
184
Table 1: Object counts in the local environmental packages before fold construction.
PSD construction population
Count
Cleaned source records
1,168,683
Unique seed structures
1,090,572
Represented SMIRKS entries, before / after curation
134 / 120
Top-five rule share, before / after curation
56.14% / 23.51%
Table 2: Construction and curation of the PSD synthetic corpus.
BBD
Soil
Sludge
Predictor
Pretraining
Single PR-AUC ↑
Multi PR-AUC ↑
Single PR-AUC ↑
Multi PR-AUC ↑
Single PR-AUC ↑
Multi PR-AUC ↑
enviGraph
PSD-Curated
0.202 ± 0.050
0.110 ± 0.029
0.154 ± 0.037
0.119 ± 0.035
0.017 ± 0.017
0.010 ± 0.004
enviFormer
PSD-Curated
0.175 ± 0.054
0.112 ± 0.037
0.130 ± 0.028
0.088 ± 0.027
0.044 ± 0.045
0.018 ± 0.021
enviGraph
USPTO-50K
0.107 ± 0.027
0.052 ± 0.017
0.103 ± 0.024
0.098 ± 0.046
0.043 ± 0.033
0.019 ± 0.014
enviFormer
USPTO-50K
0.111 ± 0.034
0.051 ± 0.024
0.106 ± 0.019
0.071 ± 0.022
0.060 ± 0.060
0.014 ± 0.010
enviGraph
USPTO-MIT
0.148 ± 0.048
0.068 ± 0.025
0.149 ± 0.038
0.126 ± 0.057
0.112 ± 0.049
0.037 ± 0.015
Table 3: Fine-tuned single-generation and multi-generation PR-AUC.
Figure 2: PR-AUC differences between synthetic and USPTO pretraining.
BBD
Soil
Sludge
Predictor
Pretraining
Top-1 ↑
Top-5 ↑
Top-1 ↑
Top-5 ↑
Top-1 ↑
Top-5 ↑
enviGraph
PSD-Curated
0.355 ± 0.052
0.594 ± 0.077
0.316 ± 0.046
0.563 ± 0.070
0.063 ± 0.035
0.263 ± 0.083
enviFormer
PSD-Curated
0.295 ± 0.057
0.530 ± 0.056
0.268 ± 0.050
0.525 ± 0.063
0.108 ± 0.087
0.344 ± 0.109
enviGraph
USPTO-50K
0.246 ± 0.032
0.461 ± 0.053
0.248 ± 0.052
0.509 ± 0.045
0.125 ± 0.055
0.379 ± 0.139
enviFormer
USPTO-50K
0.229 ± 0.047
0.396 ± 0.069
0.243 ± 0.042
0.475 ± 0.058
0.161 ± 0.099
0.351 ± 0.150
enviGraph
USPTO-MIT
0.282 ± 0.063
0.482 ± 0.059
0.308 ± 0.052
0.580 ± 0.070
0.252 ± 0.052
0.492 ± 0.114
Table 4: Fine-tuned exact top-1 and top-5 accuracy.