SynCraft: Reasoning-Driven Molecular Optimization

• SynCraft is a molecular optimization framework that reformulates synthesizability challenges as discrete graph-editing problems.
• It integrates LLM-driven chain-of-thought reasoning with RDKit-based deterministic edits to ensure valid chemical modifications.
• Evaluations on synthesis cliff datasets demonstrate improved synthetic feasibility and structural preservation compared to template-based methods.

SynCraft is a reasoning-driven molecular optimization framework designed to address challenges in molecular synthesizability encountered by generative models in computational chemistry. Unlike prior approaches based on post-hoc filtering or projection into pre-defined templates, SynCraft formulates synthesizability optimization as a discrete graph-editing problem, leveraging LLMs for edit planning and deterministic chemoinformatics toolkits for execution. Central to its design is the navigation of “synthesis cliffs,” regions where minimal atom- or bond-level modifications produce qualitatively enhanced synthetic feasibility while maintaining structural and pharmacophoric integrity (Li et al., 23 Dec 2025).

1. Framework Design and Workflow

SynCraft incorporates a two-component workflow. First, a LLM (Gemini-2.5-Pro) is prompted using Chain-of-Thought reasoning to identify synthetic liabilities within a molecule and propose an explicit, executable edit sequence. The edit commands are JSON-encoded and specify discrete atom-, bond-, and stereochemistry-level operations. Second, a deterministic chemoinformatics toolkit (RDKit) applies these edit operations to the molecular graph, ensuring chemical validity.

Traditional generative approaches issue the top-level prompt “generate a synthesizable analog of this SMILES.” In contrast, SynCraft issues the prompt “given identified synthetic liabilities (and optional biochemical constraints), propose a minimal sequence of atom- and bond-level edit commands.” This paradigm decouples molecular reasoning (LLM) from deterministic graph transformation (RDKit).

A synthesis cliff is operationally defined as a pair $(M_\text{src}, M_\text{tgt})$ where $M_\text{tgt}$ is both highly similar (in ECFP4-Tanimoto space) and synthesizable, whereas $M_\text{src}$ is unsynthesizable. Figure 1 of the reference provides qualitative illustrations of synthesis cliff navigation via edit sequences.

2. Formal Problem Specification

Molecules are modeled as labeled graphs $G = (V, E)$, with atoms $V=\{v_1, \ldots, v_n\}$ carrying element labels $\ell(v_i)\in\{\text{C},\text{N},\ldots\}$ and bonds $E \subset V \times V$ annotated with orders $o(e) \in \{1,2,3,\text{aromatic}\}$. An edit sequence $\mathbf{a} = (a_1, a_2, \ldots, a_T)$, $a_t \in \mathcal{A}$, is drawn from a defined action space $M_\text{tgt}$0.

The synthesizability optimization objective is formalized as: $M_\text{tgt}$1 where $M_\text{tgt}$2 denotes sequential application of $M_\text{tgt}$3 to $M_\text{tgt}$4, $M_\text{tgt}$5 is the ECFP4-Tanimoto similarity, and $M_\text{tgt}$6 is a binary indicator denoting whether SimpRetro can find a valid retrosynthetic route within 30 minutes. The approach is extensible to continuous synthetic accessibility scores $M_\text{tgt}$7.

3. Edit Space and LLM-Driven Editing

3.1 Edit Operations

The edit action set $M_\text{tgt}$8 includes:

• $M_\text{tgt}$9: Remove atom $M_\text{src}$0 and its bonds.
• $M_\text{src}$1: Change atom $M_\text{src}$2’s element.
• $M_\text{src}$3: Introduce new atom $M_\text{src}$4 with element $M_\text{src}$5.
• $M_\text{src}$6: Add a bond.
• $M_\text{src}$7: Remove specified bond.
• $M_\text{src}$8: Adjust bond parameters/aromaticity.
• $M_\text{src}$9 and $G = (V, E)$0: Assign stereochemistry.

These operations are parameterized by atom-map identifiers and, where applicable, bond orders and stereochemical descriptors.

3.2 Prompt Engineering and Inference Pipeline

SynCraft uses retrieval-augmented few-shot prompting. For a given input, it retrieves the top-$G = (V, E)$1 ($G = (V, E)$2) most relevant synthesis cliff pairs $G = (V, E)$3 and formats each example as:

• Source SMILES string
• Chain-of-Thought reasoning statement describing liabilities
• JSON array with the edit sequence

The LLM is instructed to output first an analysis, then the sequence of edits. The core inference pipeline is captured in the following pseudocode:

$\ell(v_i)\in\{\text{C},\text{N},\ldots\}$0

Prompt instructions enforce concise natural-language reasoning followed by edit command output, minimizing syntactic deviations.

4. Optimization, Constraints, and Loss Formulation

Interaction-aware prompting extends SynCraft to incorporate biological constraints derived from structure-based data. A PLIP-derived interaction profile, after docking the ligand into a target protein, is mapped to atom-map indices and described per atom, e.g., $\ell(v_i)\in\{\text{C},\text{N},\ldots\}$1 This is prepended under “[CRITICAL BIOLOGICAL CONSTRAINTS],” explicitly instructing the LLM to preserve or only bioisosterically modify these atoms.

No neural loss is trained; instead, explicit prompting guides the LLM to optimize for

• Structural fidelity (ECFP4-Tanimoto similarity)
• Synthetically feasible output $G = (V, E)$4
• Optional bio-constraints

Ablation results indicate edit-sequence generation avoids “SMILES hallucinations” common with direct sequence generation, yielding higher output validity and similarity at each threshold.

5. Datasets, Benchmarks, and Comparative Evaluation

SynCraft was evaluated using:

• The Synthesis Cliff dataset, containing 3,332 $G = (V, E)$5 pairs from GenBench3D outputs (five generative models), matched to eMolecules neighbors $G = (V, E)$6, $G = (V, E)$7.
• Test sets: 1,025 unsynthesizable Pocket2Mol (P2M) compounds and 1,725 ResGen outputs, both labeled by SimpRetro feasibility.

Baselines included ChemProjector, SynFormer, and ReaSyn (projection-based). The primary metric is the fraction of cases for which a synthesizable analog $G = (V, E)$8 satisfies $G = (V, E)$9 and $V=\{v_1, \ldots, v_n\}$0 for $V=\{v_1, \ldots, v_n\}$1.

Dataset	Similarity $V=\{v_1, \ldots, v_n\}$2	SynCraft (%)	Best Baseline (%)
Pocket2Mol	$V=\{v_1, \ldots, v_n\}$3	42.7	$V=\{v_1, \ldots, v_n\}$430
Pocket2Mol	$V=\{v_1, \ldots, v_n\}$5	28.4	15.4
ResGen	$V=\{v_1, \ldots, v_n\}$6	44.7	37.5
ResGen	$V=\{v_1, \ldots, v_n\}$7	29.1	22.0

Across thresholds, SynCraft consistently outperforms baselines, especially in the medium similarity regime.

6. Medicinal Chemistry Case Studies

6.1 PLK1 Inhibitor Editing (Retrospective)

For input “lig-886” (bearing a 2,5-dimethylpiperazine core causing two chiral centers), SynCraft’s output chain-of-thought was: “Remove methyl groups to avoid stereochemical mixtures.” The executed DEL_BOND and MUTATE_ATOM edits yielded a simple piperazine, matching the human-synthesized IIP0944 lead’s scaffold.

6.2 RIPK1 Candidate Rescue (Prospective)

Case I: For a biaryl linkage between electron-deficient rings, SynCraft suggested inserting an ether via ADD_ATOM (O) and CHANGE_BOND. This enabled a C–O bond formation pathway in retrosynthesis, bypassing problematic cross-coupling, and preserved H-bonding (Vina docking: $V=\{v_1, \ldots, v_n\}$8 kcal/mol).

Case II: Modifying a fused diazepino-purine to a 2,4-disubstituted pyrimidine (scaffold hop), SynCraft’s edits preserved substituent vectors and H-bond acceptor at N1. Docking showed maintained binding efficiency ($V=\{v_1, \ldots, v_n\}$9 kcal/mol).

7. Limitations and Future Research

SynCraft reframes synthesizability as a targeted, minimal graph-editing problem. This approach outperforms template-projection methods, which may introduce large topological distortions. Its modular architecture—separating LLM-driven reasoning from deterministic chemoinformatics execution—guarantees valid chemical outputs while enabling explicit pharmacophore preservation via interaction-aware prompting.

Reported limitations include operational costs (approximately $3–5 per 100 optimizations), limiting suitability for ultra-large library enumeration. Proposed future directions include:

• Integration of a learned, continuous feasibility–fidelity loss
• Automated retrosynthetic feedback loops
• Extension to multi-objective optimization (e.g., ADMET properties)
• Adaptation to broader graph-structured domains

Source code, datasets, and full prompt templates are publicly available at github.com/catalystforyou/SynCraft-Core (Li et al., 23 Dec 2025).