[edkerk]

Kcats can now come from the EC tree, and the default kcat assignment preset is brenda_ec_tree: custom kcats, the organism's own BRENDA value, the EC tree, then DLKcat. geckopy has the same (CompMetEng/geckopy#67). This is step K4a of the plan in CompMetEng/raven-gecko-parity#134.

BRENDA measurement table

databases/measurements.tsv holds one row per wild-type kcat measurement in BRENDA 2026.1 with an EC number:

• 57,037 rows, with EC number, substrate (name and InChIKey), organism, UniProt accession and sequence id, kcat, temperature, pH and PubMed id;
• kcat in 1/s, the geometric midpoint for a range;
• values more than four orders of magnitude below their EC class, and clusters of values below 0.01 /s that sit three orders of magnitude under the other values of their EC number, are left out as unit errors (423 values).

loadBrendaMeasurements reads it from the BRENDA folder of the model adapter.

EC tree

ecTreeEstimate gives an EC number the median log10 kcat of its measurements, from any organism and substrate, pulled toward its parent EC classes with weight n / (n + k):

value = median of all measurements
for level in (2-digit, 3-digit, 4-digit EC):
    w     = n / (n + k)               # n = measurements at that level, k = 5
    value = w * median(level) + (1 - w) * value

A well-measured EC keeps its own median; a sparsely measured one leans on its class. A reaction with several EC numbers uses the one with most measurements.

kcatsFromECTree writes the values to data/kcats/ec_tree.tsv, with the number of measurements and the deepest EC level used in the detail column. Reactions without an EC number, or whose EC classes have no measurements, get no value. saveKcatSource takes further detail pairs from a detail field.

Preset brenda_ec_tree, now the default

Order	Source
1	custom kcats
2	BRENDA median of measurements of the organism on the reaction's substrate (origin 1, no EC wildcard)
3	EC tree
4	DLKcat, for reactions without an EC number
then	isozymes, standard kcat

It is the default when [kcat_assignment] is left out and in adapterTemplate.toml; gecko_median and gecko_classic stay available, and the full_ecModel tutorial keeps gecko_classic. updateKcatSources gathers ec_tree.tsv and marks it outdated for another k or another BRENDA release.

Why this order

Against measured kcats hidden from the search (enzyme–substrate pairs whose own measurements and those of sequences at ≥ 95% identity are hidden), the EC tree was closer than the GECKO median at origins 2 and 4 and with EC wildcards, by 0.13 to 0.18 log10, and as close at origins 1 and 3. DLKcat was not closer than the BRENDA median at any origin. The organism's own measurement comes first because, in a model, it is often a measurement of that very enzyme. The figures are in CompMetEng/raven-gecko-parity#134.

On the yeast tutorial model before tuning:

Preset	Geometric mean kcat (/s)	Standard kcat	Maximum growth (/h)
gecko_classic	12.7	259	0.082
gecko_median	3.8	259	0.081
brenda_ec_tree	5.3	259	0.082

With brenda_ec_tree, 3,574 reactions take the EC tree, 577 DLKcat and 210 the organism's own BRENDA value.

Tests

ecTreeTests (tc0084–tc0087):

• a worked example whose values are computed by hand, the same as geckopy's;
• the shipped measurement table;
• the file kcatsFromECTree writes;
• the default preset, its order of sources, and updateKcatSources for another k.

ecTreeEstimate gives the same values as geckopy's ECTree on 4,885 EC numbers and combinations from the measurement table.

[github-actions]

Unit test results

61 tests   61 ✅  11s ⏱️
 1 suites   0 💤
 1 files     0 ❌

Results for commit 776e9fd.

♻️ This comment has been updated with latest results.