RetroCast is a comprehensive toolkit for standardizing, scoring, and analyzing multistep retrosynthesis models. It decouples prediction from evaluation, allowing rigorous, apples-to-apples comparison of disparate algorithms on a unified playing field.
The field of retrosynthesis is fragmented.
- Incompatible Outputs: AiZynthFinder outputs bipartite graphs; Retro* outputs precursor maps; DirectMultiStep outputs recursive dictionaries. Comparing them requires writing bespoke parsers for every paper.
- Ad-Hoc Metrics: "Solvability" is often calculated differently across publications, with varying definitions of commercial stock (e.g., using made-to-order libraries vs. actual off-the-shelf compounds).
- Flawed Benchmarks: The standard PaRoutes n5 dataset is heavily skewed (74% of routes are length 3-4), masking performance failures on complex targets. Furthermore, the standard stock definition for PaRoutes creates synthetic "ground truths" that are often physically unobtainable.
RetroCast solves this. It provides a canonical schema, adapters for 10+ models, and a rigorous statistical pipeline to turn retrosynthesis from a qualitative art into a quantitative science.
- Universal Adapters: "Air-gapped" translation layers for AiZynthFinder, *Retro**, DirectMultiStep, SynPlanner, Syntheseus, ASKCOS, RetroChimera, DreamRetro, MultiStepTTL, SynLlama, and *PaRoutes*.
- Canonical Schema: All routes are cast into a strict schema-2
Molecule/Reactiontree. - Rust Core: Every built-in planner adapter, ingest, scoring, and analysis runs in Rust, with canonical SMILES, InChIKeys, and molecular descriptors calculated through RDKit C++.
- Curated Benchmarks: Includes the Reference Series (for algorithm comparison) and Market Series (for practical utility), stratified by route length and topology to eliminate statistical noise.
- Rigorous Statistics: Built-in bootstrapping (95% CI), pairwise tournaments, and probabilistic ranking. No more "Model A is 0.1% better than Model B" without significance testing.
- Reproducibility: Every artifact is tracked via cryptographic manifests (
SHA256).
The Python package includes the native engine and its RDKit C++ libraries; it does not install Python RDKit. It is the right installation when you want both import retrocast and the compatibility CLI:
# Install the Python package as a command-line tool
uv tool install retrocast
# Or add the importable library to a project
uv add retrocastThe standalone Rust executable is published separately on the GitHub Releases page. Download the archive for linux-x86_64, macos-arm64, macos-x86_64, or windows-x86_64; each archive contains retrocast and the RDKit libraries it was built and smoke-tested with. It does not require Python, Conda, Rust, or a separately installed RDKit.
Markdown docs are formatted with oxfmt and kept unwrapped (proseWrap: "never"). Use pnpm docs:fmt to normalize docs formatting or pnpm docs:fmt:check to verify it in review/CI.
Latest Data (Updated Regularly)
For the most up-to-date benchmarks and stocks, use the RetroCast CLI:
# Download a specific benchmark (includes definition + required stock)
retrocast get-data mkt-cnv-160Publication Data (Frozen)
The complete data/ folder as used in the preprint is available at files.ischemist.com/retrocast/publication-data:
# Show available folders and their sizes
curl -fsSL https://files.ischemist.com/retrocast/get-pub-data.sh | bash -s
# Download all benchmark definitions
curl -fsSL https://files.ischemist.com/retrocast/get-pub-data.sh | bash -s -- definitionsyou can verify the integrity of downloaded files against the manifests by running
retrocast verify --allthat command might warn you about missing files---that is expected. Manifests for, say 4-scored, contain hashes of input files from 3-results, and if you downloaded only 4-scored, you will get warnings about missing 3-results files.
a dump of the sqlite db with the stocks, routes, and results loaded into SynthArena can be found in ischemist/syntharena repo.
Have a raw output file from a model? Score it immediately.
# Convert raw AiZynthFinder JSON to RetroCast format
retrocast adapt \
--input raw_predictions.json.gz \
--adapter aizynthfinder \
--output candidates.json.gz
# Align predictions to benchmark-keyed candidates
retrocast collect \
--input candidates.json.gz \
--benchmark data/retrocast/1-benchmarks/definitions/ref-lin-600.json.gz \
--output collected-candidates.json.gz
# Score against a stock file
retrocast score-file \
--benchmark data/retrocast/1-benchmarks/definitions/ref-lin-600.json.gz \
--candidates collected-candidates.json.gz \
--stock data/retrocast/1-benchmarks/stocks/n5-stock.csv.gz \
--output scores.json.gz \
--model-name "My-Experimental-Model"For full-scale benchmarking, RetroCast enforces a structured data lifecycle: Ingest Score Analyze.
model-specific prediction runners were extracted from this repo and now live in project-pandora. use RetroCast here for standardization, scoring, and analysis; use pandora when you need the upstream model-runner scripts.
Inspect project paths:
retrocast configProject mode reads raw model outputs from data/retrocast/2-raw/<model>/<benchmark>/. Pass the adapter on the command line, or put a manifest.json next to the raw results file:
{
"directives": {
"adapter": "directmultistep",
"raw_results_filename": "predictions.json"
}
}Run the pipeline:
# 1. Ingest: Standardize raw outputs from data/retrocast/2-raw/
retrocast ingest --model dms-explorer --dataset ref-lin-600 --adapter directmultistep
# 2. Score: Evaluate against the benchmark's defined stock
retrocast score --model dms-explorer --dataset ref-lin-600
# 3. Analyze: Generate bootstrap statistics
retrocast analyze --model dms-explorer --dataset ref-lin-600The all-in-one command keeps Rust-owned routes in memory between stages and uses one bounded worker pool:
retrocast pipeline \
--raw data/retrocast/2-raw/aizynthfinder-run/mkt-cnv-160 \
--adapter aizynthfinder \
--benchmark data/retrocast/1-benchmarks/definitions/mkt-cnv-160.json.gz \
--stock data/retrocast/1-benchmarks/stocks/buyables-stock.csv.gz \
--output-dir native-results \
--workers 12Output: A schema-v2 analysis report in data/retrocast/5-results/.
RetroCast introduces two new benchmark series derived from PaRoutes, fixing the skew and stock issues of the original dataset. These subsets were selected via seed stability analysis to ensure they are statistically representative of the underlying difficulty distribution.
Target Audience: Algorithm Developers Designed to compare search algorithms (e.g., MCTS vs. Retro* vs. Transformers). Uses the internal PaRoutes stock to isolate search failures from stock availability issues.
| Benchmark | Targets | Description |
|---|---|---|
| ref-lin-600 | 600 | Linear routes stratified by length (100 each for lengths 2–7). |
| ref-cnv-400 | 400 | Convergent routes stratified by length (100 each for lengths 2–5). |
| ref-lng-84 | 84 | All available routes of extreme length (8–10 steps). |
Target Audience: Computational Chemists Designed to assess practical utility. Targets are filtered to be solvable using Buyables, a curated catalog of 300k compounds available for <$100/g.
| Benchmark | Targets | Description |
|---|---|---|
| mkt-lin-500 | 500 | Linear routes solvable with commercial buyables (Stratified). |
| mkt-cnv-160 | 160 | Convergent routes solvable with commercial buyables (Stratified). |
RetroCast is also a library. You can use it to integrate standardization directly into your training or inference loops.
adapt, ingest_candidates, score, and analyze use the same Rust core as the standalone retrocast executable for every built-in adapter. There is no Python/Rust engine selector: Python is an interface to the same implementation.
from retrocast import adapt
from retrocast.adapters import DirectMultiStepAdapter
# Your model's raw output (any supported format)
raw_output = {
"smiles": "CC(=O)Oc1ccccc1C(=O)O",
"children": [...]
}
# Cast to the canonical Route object
adapter = DirectMultiStepAdapter()
candidate = adapt([raw_output], adapter)[0]
route = candidate.route
if route is None:
raise ValueError("Could not adapt route")
print(f"Depth: {route.depth()}")
print(f"Leaves: {[m.smiles for m in route.leaves()]}")RetroCast powers SynthArena, an open-source web platform for visualizing and comparing retrosynthetic routes.
- Compare predictions from any two models side-by-side.
- Visualize ground truth vs. predicted routes with diff overlays.
- Inspect stratified performance metrics interactively.
We distinguish between two fundamental classes of problems in scientific machine learning: quantitative (predicting scalar targets like toxicity or binding affinity) and structural (generating complex objects governed by an underlying grammar). Quantitative problems, analogous to early NLP challenges like sentiment analysis, are often constrained by data scarcity. In contrast, the most transformative AI breakthroughs—from large language models to AlphaFold—have occurred in structural domains.
Mastery of structure is a prerequisite for solving downstream quantitative tasks. Foundation models trained on the structure of language, for instance, now excel at sentiment analysis with little to no task-specific fine-tuning. In organic chemistry, the paramount structural challenge is retrosynthesis: designing a valid synthetic pathway to a molecule of interest. This capability is the key to unlocking critical quantitative problems like predicting synthetic accessibility, a significant bottleneck in drug discovery. Current accessibility heuristics, however, bypass the core structural challenge, relying on learned patterns that correlate with accessibility without ever generating the pathway itself.
A model cannot judge the difficulty of a journey it cannot first articulate.
Achieving structural mastery in retrosynthesis is a long journey—one that requires moving beyond fragmented data formats, inconsistent evaluation methods, and unreliable metrics. Progress demands unified, rigorous infrastructure to standardize outputs, track provenance, and measure improvements with statistical rigor.
RetroCast is that infrastructure.
If you use RetroCast in your research, please cite:
@misc{retrocast,
title = {Procrustean Bed for AI-Driven Retrosynthesis: A Unified Framework for Reproducible Evaluation},
author = {Anton Morgunov and Victor S. Batista},
year = {2025},
eprint = {2512.07079},
archiveprefix = {arXiv},
primaryclass = {cs.LG},
url = {https://arxiv.org/abs/2512.07079}
}MIT License. See LICENSE for details.