KI-Retrosynthese: Routen, die sich wirklich umsetzen lassen

Jedes Retrosynthese-Tool liefert einen Baum. Die Frage ist, ob irgendein Ast davon den Kontakt mit dem Abzug übersteht. Eine Route mit Score 1.00 und eine Route, die ein Chemiker tatsächlich fährt, sind zwei verschiedene Dinge – und in dieser Lücke gehen die meisten CASP-Bewertungen schief.

Dieser Artikel begleitet einen echten Lauf von chemreason_retrosynthesis an Celecoxib, zeigt die Rohausgabe und trennt die Vorschläge, die ein Prozesschemiker direkt an die Bench nehmen könnte, von denen, die zuerst eine menschliche Korrektur brauchen.

Was eine KI-Retrosynthese-Engine tatsächlich leistet

In jedem Planer stecken zwei getrennte Aufgaben, und sie versagen unterschiedlich.

• Template-basiert. Reaktionsregeln, die aus einem Korpus wie USPTO oder Reaxys extrahiert wurden, werden auf das Ziel angewendet und rückwärts abgeglichen. Vorhersehbar, chemisch belesen, auditierbar – und strukturell unfähig, eine Diskonnektion vorzuschlagen, die niemand aufgezeichnet hat.
• Template-frei / ML. Ein neuronales Modell lernt die Atom- und Bindungsumlagerung direkt. ChemReason verwendet einen Graph Transformer über Bond-Electron-Matrizen, sodass eine Reaktion eine Änderung der Elektronenbuchhaltung ist statt einer String-Bearbeitung. Das generalisiert über die Template-Bibliothek hinaus und erzeugt gelegentlich etwas mit unmöglicher Valenz.

Die wichtige Konsequenz: Ein hoher Route-Score ist eine Aussage über den *Suchabbruch*, nicht über Ausbeute, Selektivität oder Skalierbarkeit.

Ein echter Lauf: Celecoxib

Ziel-SMILES Cc1ccc(-c2cc(C(F)(F)F)nn2-c2ccc(S(N)(=O)=O)cc2)cc1, 26 Schweratome, max_steps=4, n_proposals=5.

Zunächst das mehrstufige Ergebnis, wörtlich:

n_routes: 1
routes[0]: {depth: 0, score: 1, is_complete: true, steps: [],
            building_blocks: ["Cc1ccc(-c2cc(C(F)(F)F)nn2-c2ccc(S(N)(=O)=O)cc2)cc1"]}

Der Planer lieferte Tiefe 0 mit perfektem Score und ohne Schritte. Das ist kein Fehlschlag. Celecoxib befindet sich im 11.4M-Bausteinkatalog, also lautet die korrekte retrosynthetische Antwort "kaufen". Ein Planer, der stattdessen eine vierstufige Route zu einer käuflichen Verbindung erfunden hätte, wäre schlechter, nicht besser.

Dasselbe passierte bei Paracetamol (CC(=O)Nc1ccc(O)cc1, 11 Atome): n_routes: 1, Tiefe 0, Score 1.0. Beide sind Katalogware.

Das ist die erste praktische Lehre. Prüfen Sie den Katalog, bevor Sie irgendeinen Baum lesen. Wenn Sie dennoch Diskonnektionen für ein käufliches Ziel wollen – aus Kosten-, IP- oder Prozessgründen –, nimmt chemreason_retrosynthesis den Parameter force_decompose=true, und die Einschritt-Vorschläge werden trotzdem zurückgegeben.

Die Einschritt-Vorschläge

Für Celecoxib lieferte das Modell n_valid: 5, n_unique: 5. Alle fünf, unbearbeitet:

Zwei von fünf sind direkt ausführbar. Drei brauchen einen Chemiker, der das Reagenz oder die Oxidationsstufe korrigiert, bevor sie irgendetwas bedeuten.

Bemerkenswert ist, was *fehlt*. Alle fünf Vorschläge diskonnektieren das Sulfonamid an der Peripherie. Keiner diskonnektiert den Pyrazolring – dorthin führt die tatsächliche Herstellroute: ein 1,3-Diketon, kondensiert mit 4-Hydrazinylbenzolsulfonamid. Das Modell fand die einfache Bindung, nicht die strategische. Das ist ein sehr verbreiteter Fehlermodus bei CASP-Tools, und deshalb ist eine Einschritt-Liste eine Anregung für einen Chemiker, kein Plan.

Schnelltriage: Ist das überhaupt herstellbar?

Vollständiges MCTS braucht Minuten an GPU-Zeit. Zum Filtern einer Bibliothek führt chemreason_analyze fünf Einschritt-Vorschläge aus, prüft jede Vorstufe gegen den Katalog und liefert eine Trefferquote zu etwa einem Siebtel der Kosten.

bb_hit_rate ist der Anteil der vorgeschlagenen Vorstufen, die bereits im 11.4M-Katalog stehen. Das ist ein schneller Synthetisierbarkeits-Prior zur Rangfolge einer generierten Bibliothek. Es ist keine Route, und der Ibuprofen-Top-Vorschlag beweist genau das: Die Fragmente sind katalognah, aber die implizierte Friedel-Crafts-Reaktion mit einer alpha-Brom-alpha-hydroxy-Säure ist nicht, wie Ibuprofen hergestellt wird. Trefferquote hoch, Chemie falsch.

Route-Score vs. Route-Umsetzbarkeit

Ein Modell-Score und das Urteil eines Chemikers beantworten unterschiedliche Fragen.

Die letzte Zeile ist nicht theoretisch. Unser Paracetamol-Lauf mit Standardeinstellungen lieferte einen Einschritt-Vorschlag C.CC12CC(Cl)(NC=O)C1(C)C2.C[SH](=O)=O.O – ein bicyclisches Chlor-Formamid ohne sinnvollen Bezug zum Ziel. n_valid: 1 zählt RDKit-parsebar, nicht chemisch sinnvoll. An einem 26-Atom-arzneistoffähnlichen Ziel verhielt sich das Modell gut (5/5 eindeutig und größtenteils interpretierbar); an einem 11-Atom-Katalogartikel erzeugte es Rauschen. Kleine, dicht funktionalisierte Ziele sind die Schwachstelle, und das Tool warnt nicht davor.

Wie man einen Planer nutzt, ohne sich selbst zu täuschen

1. Zuerst den Katalog prüfen. Kommt Tiefe 0 zurück, stoppen. Kaufen ist die Route. 2. Die Einschritt-Liste lesen, bevor man den Baum liest. Der Diskonnektionssatz zeigt, welche Strategie das Modell hat, und ob es die strategische oder nur die bequeme Bindung fand. 3. Jedes Vorstufen-SMILES als Entwurf behandeln. Standardisieren, Valenzen und Oxidationsstufe prüfen, dann entscheiden. 4. `chemreason_analyze` für Bibliotheken, `chemreason_retrosynthesis` für Shortlists nutzen. Die Triage läuft zu etwa einem Siebtel der Kosten einer vollständigen Routensuche pro Molekül – der Unterschied zwischen dem Screening von tausend und zehn Molekülen. 5. Die Route nach den eigenen Randbedingungen bewerten, nicht nach denen des Modells. Bausteinverfügbarkeit im Maßstab, ICH-Q3C-Lösungsmittelklasse, genotoxische Zwischenstufen und ob das eigene Kit die Chemie überhaupt beherrscht.

Was das nicht sagt

• In diesen Läufen wurden keine Bedingungen zurückgegeben. Das Ausgabeschema führt reaction_type, reagents, solvent, temperature und yield_estimate pro Schritt, aber beide Ziele endeten bei Tiefe 0, sodass es keine Schritte zum Befüllen gab. Zitieren Sie keine vorhergesagte Ausbeute, die Sie nicht gesehen haben.
• `n_valid` bedeutet parsebar, nicht plausibel. Chemische Plausibilität ist Ihre eigene Prüfung.
• Der Katalog ist eine Momentaufnahme, nicht Ihre Lieferkette. Ein Bausteintreffer sagt, dass ein Anbieter ihn gelistet hat – nicht, dass er in der benötigten Menge, Reinheit oder Frist liefert.
• Das Modell kennt Ihre Prozessrandbedingungen nicht. Nichts davon berücksichtigt GMP, Gefahrenklassifizierung, PMI oder die regulatorische Geschichte einer Route.
• Retrosynthese ist keine Green-Chemistry-Bewertung. Lösungsmittelauswahl und Verunreinigungs-Überlegungen separat durchführen.
• Das ist Triage. Sie zeigt, wo ein Chemiker seine Aufmerksamkeit einsetzen sollte. Sie ersetzt kein Route-Scouting und ist keine Zulassungsposition.

Häufige Fragen

Was ist KI-Retrosynthese-Planung?

KI-Retrosynthese-Planung arbeitet rückwärts von einem Zielmolekül zu käuflichen Ausgangsstoffen. Ein Einschritt-Modell schlägt Vorstufen für eine Diskonnektion vor, und ein Suchalgorithmus verkettet diese Schritte zu vollständigen Routen. ChemReason verwendet für die Einschritt-Vorhersage einen Graph Transformer über Bond-Electron-Matrizen und für die mehrstufige Suche eine Monte-Carlo-Baumsuche über 11.4M kommerzielle Bausteine, mit Routen inklusive Vollständigkeits-Flag und Score.

Was unterscheidet template-basierte von maschinell gelernter Retrosynthese?

Template-basierte Systeme extrahieren Reaktionsregeln aus Reaktionsdatenbanken und wenden sie rückwärts an. Sie sind auditierbar und chemisch konservativ, können aber keine Diskonnektion außerhalb des Regelsatzes vorschlagen. Machine-Learning-Systeme lernen die Bindungs- und Elektronenumlagerung direkt, generalisieren also auf unbekannte Chemie – auf Kosten gelegentlich chemisch unmöglicher Vorstufen. In der Praxis braucht die ML-Ausgabe bei jeder Reagenzstruktur eine chemische Plausibilitätsprüfung durch einen Chemiker.

Warum lieferte das Retrosynthese-Tool keine Schritte für mein Molekül?

Weil das Ziel bereits ein käuflicher Baustein ist. In unserem Lauf lieferten sowohl Celecoxib (26 Schweratome) als auch Paracetamol (11 Schweratome) n_routes: 1 bei depth: 0 mit score: 1 und einer leeren Schrittliste von chemreason_retrosynthesis. Der Planer stoppt, wenn das Ziel im 11.4M-Bausteinkatalog liegt. Werden dennoch Diskonnektionen benötigt, force_decompose=true setzen.

Bedeutet ein hoher Retrosynthese-Score, dass die Route gut ist?

Nein. Der Score zeigt, dass die Suche bei flacher Tiefe bei katalogverfügbarem Material endete. Er enthält keine Information über Ausbeute, Selektivität, Verträglichkeit funktioneller Gruppen, Kosten im Maßstab, Gefahrenklasse oder Durchsatz. In unserem Celecoxib-Lauf waren nur zwei von fünf Einschritt-Vorschlägen direkt unbearbeitet ausführbar; die anderen drei hatten eine nicht standardmäßige Reagenzform oder eine falsche Oxidationsstufe. Die Chemiker-Prüfung ist nicht optional.

Wie kann ich schnell prüfen, ob ein generiertes Molekül synthetisierbar ist?

Vor der vollständigen Routensuche eine schnelle Zugänglichkeitsprüfung nutzen. chemreason_analyze liefert eine Bausteintrefferquote über fünf Einschritt-Vorschläge: Ibuprofen erzielte bb_hit_rate 0.57 und Zugänglichkeit "high", Celecoxib 0.50 und "medium". Das läuft zu etwa einem Siebtel der Kosten der vollständigen Planung, wodurch es für eine generierte Bibliothek nutzbar wird. Als Prior für die Rangfolge behandeln, nicht als Nachweis, dass eine Route existiert.

Findet ein KI-Planer die industrielle Route zu einem bekannten Wirkstoff?

Nicht zuverlässig. Für Celecoxib diskonnektierten alle fünf Vorschläge von chemreason_retrosynthesis das periphere Sulfonamid; keiner schlug den Pyrazolring-Aufbau aus einem 1,3-Diketon und einem Arylhydrazin vor, den die Herstellroute nutzt. Planer sind zu häufigen, spät ansetzenden, einfachen Diskonnektionen tendenziell verzerrt. Nutzen Sie sie, um das eigene Optionsfeld zu erweitern, nicht um strategisches Route-Design zu ersetzen.

Weiterführende Artikel

• Generatives Moleküldesign, das ADMET übersteht
• Applicability Domain: Warum ein Modell sagen können muss "ich weiß es nicht"
• Lösungsmittelauswahl nach ICH Q3C, als Datenfrage
• Was ein Docking-Score von -10 kcal/mol bedeutet – und was nicht

Testen Sie chemreason_analyze zuerst an Ihren eigenen Gerüsten – es ist im kostenlosen CovaSyn-Tarif verfügbar und zeigt innerhalb einer Minute, ob sich der volle Planer lohnt.