Identifizierung unbekannter Verunreinigungen mittels LC-MS/MS

Sie haben einen Peak bei 0.19 Flächenprozent, der in der vorherigen Charge nicht da war. ICH Q3A(R2) verlangt, dass alles bei oder über 0.10% identifiziert werden muss, sodass "unbekannt bei RRT 1.12" kein akzeptabler Endzustand ist. Sie haben eine akkurate Masse und ein MS/MS-Spektrum, und keinen Referenzstandard. Das ist der Workflow, der Sie von diesem Spektrum zu einer vertretbaren Strukturhypothese bringt, und die Punkte, an denen er aufhört, zu helfen.

Die Reihenfolge, in der die Arbeit tatsächlich abläuft

1. Formel zuerst. Akkurate Masse plus Isotopenmuster grenzen die Elementzusammensetzung ein. Nichts danach ist vertrauenswürdig, wenn die Formel falsch ist. 2. RDBE und das Delta zur Ausgangssubstanz. Der Sättigungsgrad sagt Ihnen, ob ein Atom hinzugefügt oder ein Ring geschlossen wurde. 3. Fragmentlogik. MS/MS sagt Ihnen, welcher Teil des Moleküls sich verändert hat. 4. Bibliothekssuche. Schnell, wenn sie funktioniert, uninformativ, wenn die Verbindung nicht in der Bibliothek ist. 5. Bestätigung. Referenzstandard, orthogonale Daten oder Synthese. Diesen Schritt übernimmt keine Software.

Schritte 1 bis 4 sind Stunden Arbeit, die sich auf Minuten komprimieren lassen. Schritt 5 dauert weiterhin Wochen.

Das durchgerechnete Beispiel

Ausgangswirkstoff: Propranolol, [M+H]+ 260.1645. Der neue Peak hat [M+H]+ 276.1594 auf einem Q-TOF, kalibriert auf besser als 3 ppm. Das Delta zur Ausgangssubstanz beträgt +15.9949, die klassische Sauerstoffaddition.

Ein Hinweis zur Herkunft vor den Zahlen: Die MS/MS-Peakliste unten ist konstruiert aus publiziertem Fragmentierungsverhalten von Betablockern, nicht an einem Instrument aufgenommen. Jeder zitierte Tool-Output ist ein echter Aufruf gegen die produktiven CovaMS-Tools. Die Peakliste ist die Eingabe; die Ausgaben sind nicht simuliert.

Übergebene Peakliste (m/z, relative Intensität): 276.1594 (100), 258.1489 (25), 199.0754 (40), 183.0804 (10), 157.0648 (15), 116.1070 (80), 98.0964 (30), 74.0600 (20).

Schritt 1: Formel aus akkurater Masse

covams_formula bei [M+H]+, 5 ppm, lieferte drei Kandidaten:

Wortgetreu aus einem covams_formula-Aufruf. Das Tool wendet vor der Bewertung die Senior-Valenzregeln und die Kind-Fiehn-Heuristiken an und rankt dann anhand einer Bayes'schen Kombination aus Massenfehler, einem RDBE-Prior und Element-Verhältnis-Priors, die auf öffentliche Trainingsdaten gefittet wurden. Die zweitplatzierte Formel liegt nur 2.5 ppm entfernt, aber sieben Stickstoffe und ein Schwefel auf 275 Da sind für einen kleinmolekularen Abbaustoff unplausibel, und der Bayes-Term bestraft das um vier Größenordnungen.

Schritt 2: Das Isotopenmuster killt die Alternativen

Akkurate Masse allein unterscheidet nicht zwischen CHNO und CHNOS. Die Isotopenhülle schon. covams_isotope für die beiden führenden Kandidaten, einfach geladen:

Wortgetreu aus zwei covams_isotope-Aufrufen. Das M+1-Verhältnis zählt Kohlenstoffe (rund 1.1% pro Kohlenstoff), und das M+2 bei 5.51% ist die Schwefelsignatur. Liegt Ihr gemessenes M+2 bei etwa 2%, ist die schwefelhaltige Formel ohne weitere Arbeit erledigt. Diese Prüfung kostet nichts und ist das ertragreichste, was Sie mit einem bereits aufgenommenen Full-Scan-Spektrum tun können.

Schritt 3: RDBE sagt Ihnen, welche Art von Veränderung stattgefunden hat

covams_formula berichtet RDBE als 1 + C - (H + Cl + Br + F)/2 + N/2 + P/2. Für C16H21NO3 lieferte es 7.0. Propranolol selbst (C16H21NO2) hat ebenfalls RDBE 7.

Dieser einzelne Vergleich ist aufschlussreich. Die Verunreinigung hat einen Sauerstoff gewonnen und keinen Sättigungsgrad. Also:

• Eine Hydroxylierung (aromatisch oder aliphatisch) passt. RDBE unverändert.
• Ein N-Oxid oder Hydroxylamin passt. RDBE unverändert.
• Ein Keton aus der Oxidation des sekundären Alkohols passt nicht. Dieser Pfad verliert zwei Wasserstoffe, was ein Massendelta von -2.0157 und RDBE 8 ergibt.
• Ein neuer Ring passt ebenfalls nicht. RDBE würde steigen.

Sie haben ein +16-Da-Unbekanntes auf zwei Strukturfamilien eingegrenzt, bevor Sie das MS/MS überhaupt angerührt haben.

Schritt 4: Fragmentannotation, und warum die Trefferzahl die falsche Metrik ist

covams_fragment_annotate gleicht gemessene Peaks gegen kombinatorisch erzeugte Fragmente einer Kandidatenstruktur ab, mit erlaubten gängigen neutralen Verlusten. Gegen 4-Hydroxypropranolol (Ring-Hydroxylierung) bei 10 ppm gelaufen, annotierte es 8 von 8 Peaks, fünf davon mit einer expliziten Fragmentstruktur:

• 116.1070 → CC(C)NCC(O)CO minus H2O, einfache Spaltung, Konfidenz 1.0
• 98.0964 → dieselbe Seitenkette minus 2 H2O, Konfidenz 0.898
• 74.0600 → NCC(O)CO minus H2O, gepaarte Spaltung, Konfidenz 1.0
• 157.0648 → COc1ccc(O)c2ccccc12 minus H2O, Konfidenz 1.0
• 183.0804 → CC(O)COc1ccc(O)c2ccccc12 minus 2 H2O, Konfidenz 0.945

Jetzt der ehrliche Teil. Derselbe Aufruf gegen das N-Hydroxy-Isomer annotierte ebenfalls 8 von 8. Die Coverage unterschied die beiden Strukturen überhaupt nicht. Was sie unterschied, war das Feld *mechanism*:

Methan aus einem Hydroxylamin zu verlieren, um bei 74.0600 zu landen, ist keine Fragmentierung, die jemand vor einem Review-Board verteidigen würde. Die Ring-Hydroxy-Hypothese erklärt dieselben Peaks nur mit Wasserverlusten bei einfachen Bindungsspaltungen. Lesen Sie den Mechanismus und den neutralen Verlust, nicht die Annotationszahl. Ein kombinatorischer Annotator mit großzügigen neutralen Verlusten erklärt fast alles, wenn man ihn lässt.

Der entscheidende Punkt ist der, den die Annotation impliziert: Das intakte Aminopropanol-Ion bei 116.1070 verlangt eine unveränderte Stickstoffseitenkette. Wäre der Sauerstoff an den Stickstoff gegangen, würde dieses Ion auf 132.1019 wandern. Ein Ion, eine Entscheidung.

Schritt 5: Bibliothekssuche, und ihre harte Grenze

covams_analyze_unknown führt Qualitätsbewertung, Processing, Bibliotheksidentifizierung, Klassenvorhersage und Annotation in einem Aufruf durch. Es bewertete dieses Spektrum mit B (spektrale Entropie 2.6201, normalisierte Entropie 0.8734, Precursor vorhanden) und lieferte 10 Kandidaten.

Rang 1 war Homatropin, C16H21NO3, Gesamtscore 0.597, spektrale Ähnlichkeit 0.7023, aus MassBank stammend. 4-Hydroxypropranolol kam auf Rang 6 zurück, Gesamtscore 0.2452, spektrale Ähnlichkeit 0.0709, aus index_formula stammend statt aus einem experimentellen Bibliothekseintrag.

Diese Umkehrung ist die Lektion, kein Bug. Das Bibliotheksranking wird von der Ähnlichkeit zu Spektren dominiert, die in der Bibliothek existieren. Homatropin ist eine Verbindung mit korrekter Formel und einem echten Referenzspektrum. 4-Hydroxypropranolol hatte kein experimentelles Spektrum zum Abgleich, konnte also nur als Kandidat auf Formelebene eintreten und erzielte beim Ähnlichkeitsterm nahezu null. Eine Verbindung ohne Bibliotheksspektrum kann keine Bibliothekssuche gewinnen, egal wie gut sie Ihre Fragmente erklärt. Verunreinigungen und Abbauprodukte sind genau die Verbindungen, die am unwahrscheinlichsten Bibliotheksspektren haben.

Noch ein Vorbehalt: Die Klassenvorhersage lieferte "Hydroxy compounds" mit Wahrscheinlichkeit 1.0 aus einem einzelnen regelbasierten Treffer, während die Laufzusammenfassung predicted_class: "Unknown" meldete. Behandeln Sie die Klassenschicht als groben Hinweis, und lesen Sie das Fehlen eines guten Treffers als Aussage über die Bibliothek, nicht über Ihre Probe.

Was das nicht zeigt

• Es identifiziert nichts. Der Output ist eine gerankte Hypothese. Die Bestätigung braucht einen authentischen Referenzstandard mit passender Retentionszeit und MS/MS, oder NMR an isoliertem Material. Das ist bestenfalls etwa Schymanski-Level 3 (vorläufiger Kandidat).
• Es klärt keine Regiochemie. 4-Hydroxy- gegenüber 5-Hydroxy- gegenüber 7-Hydroxy-Propranolol sind auf dieser Evidenzbasis nicht zu unterscheiden. Chromatographie gegen Standards oder NMR klärt das.
• Es quantifiziert nicht. Flächenprozent ist nicht Massenprozent. Responsefaktoren für ein Abbauprodukt sind meist unbekannt und unterscheiden sich oft um mehr als das Doppelte von der Ausgangssubstanz.
• Es kennt Ihre Matrix nicht. Hilfsstoff-Addukte, In-Source-Fragmente und Dimere sehen alle wie Unbekannte aus. Schließen Sie das aus, bevor Sie mit der Strukturarbeit beginnen.
• Es kann eine nie gesehene Verbindung nicht ranken. Siehe Schritt 5.

Häufige Fragen

Weiterführende Artikel

• Massenbilanz bei Stressbelastung: die Lücke schließen
• Extraktions- und Auslaugstoffe: ein Screening-Workflow
• Studiendesign für Stressbelastung unter ICH Q1A(R2)
• Nitrosamin-Risikobewertung: ein ICH-M7-Workflow

Jeder Aufruf in diesem Artikel läuft kostenlos in CovaSyn, sodass Sie Ihr eigenes Unbekanntes vor Ihrer nächsten OOS-Untersuchung durch covams_formula und covams_fragment_annotate laufen lassen können.