$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
GC×GC-TOF MS wzorce wysokiej jakości oliwy z oliwek z pierwszego tłoczenia wykazują około 500 pików 2D powyżej progu stosunku sygnału do szumu (SNR) wynoszącego 100. Taki próg został zdefiniowany przez poprzednie badania nad lotnymi składnikami żywności14,27 jako minimalny względny sygnał powyżej progu w celu uzyskania wiarygodnych widm do analizy porównawczej. Składniki są rozmieszczone w przestrzeni chromatograficznej zgodnie z ich względnym zachowaniem w dwóch wymiarach chromatograficznych, a w szczególności w oparciu o ich lotność/polarność w 1D i lotność w 2D. Tutaj kombinacja kolumn jest polarna × półpolarna (tj. Carbowax 20M × OV1701).
Wzorzec 2D pokazuje wysoki stopień porządku. Względne wzorce retencji dla szeregów i klas homologicznych są pokazane na Rysunek 1A z adnotacjami (grafika dla grup i bąbelki dla pików) dla liniowo nasyconych węglowodorów (), węglowodorów nienasyconych (żółty), liniowo nasyconych aldehydów (niebieski), jednonienasyconych aldehydów (czerwony), wielonienasyconych aldehydów (łosoś), pierwszorzędowych alkoholi (zielony) i krótkołańcuchowych kwasów tłuszczowych (cyjano).
Wykryte piki 2D można następnie zidentyfikować, porównując średnie widmo MS wyodrębnione z całego piku 2D (widmo plamy) lub z największego widma (widmo wierzchołkowe). Rysunek 2 ilustruje wynik wyszukiwania widma wierzchołkowego dla blob 5 i zwraca wysokie dopasowanie podobieństwa (pierwsze 10 trafień) dla (E)-2-hexenal. Eksplorowane bazy danych to bazy danych wstępnie wybrane przez analityka w kroku 8 metody.
Identyfikacja jest weryfikowana przez aktywne indeksowanie retencji. Eksperymentalna wartość IT została obliczona dla pików 2D, tak aby na tym etapie wyszukiwanie w bibliotece nadało priorytet wynikom ze spójnymi wartościami tabelarycznego I T. Okna tolerancji można dostosować w oparciu o doświadczenie analityka, wiarygodność wartości referencyjnej bazy danych zgodnie z fazą stacjonarną oraz zastosowane warunki analityczne. Nowe narzędzia do inteligentnej kalibracji liniowych wskaźników retencji bez eksperymentalnej kalibracji za pomocą n-alkanów zostały niedawno opracowane i omówione w badaniu przeprowadzonym przez Reichenbach et al19.
Zbiór zidentyfikowanych pików 2D (tj. docelowych pików) może zostać wykorzystany do zbudowania szablonu docelowych pików, aby szybko ustalić wiarygodne powiązania między tym samym związkiem na wszystkich chromatogramach próbek. Zbiór docelowych szczytów szablonu jest wizualizowany w Rysunek 1B. Czerwone kółka odpowiadają 196 docelowym związkom, w tym dwóm wzorcom wewnętrznym (IS) połączonym z pikami matrycowymi za pomocą linii łączących. IS są używane do normalizacji odpowiedzi, a linie połączeń pomagają zwizualizować, który z dołączonych IS zostanie przyjęty do normalizacji każdej odpowiedzi szczytu/plamy 2D.
W Rysunek 1B, wypełnione okręgi wskazują pozytywne dopasowania między szczytem szablonu a rzeczywistym wzorcem, podczas gdy puste kółka są dla szczytów szablonu, dla których zgodność nie została zweryfikowana. Dopasowania fałszywie ujemne mogą być ograniczone przez odpowiedni dobór parametrów progowych, widm referencyjnych i funkcji ograniczających13,14,18,19. W przypadku złożonych wzorców z wieloma koelucjami zalecane są funkcje wykrywania pików jonowych, które są oparte na dekonwolucji spektralnej i mogą być prawidłową opcją19. Metadane piku szablonu są pokazane w powiększonym panelu Rysunek 1B dla (E)-2-sześciokątnego.
Specyfika dopasowania szablonu polega na możliwości zastosowania funkcji ograniczających, które ograniczają dodatnią zgodność z tymi kandydującymi pikami, które mieszcząc się w oknie wyszukiwania algorytmu, mają podobieństwo widmowe MS powyżej pewnego progu. W tym przypadku, w kroku 11, progi podobieństwa23 zostały ustalone na 700 zgodnie z wcześniejszymi eksperymentami mającymi na celu zdefiniowanie optymalnych parametrów ograniczających fałszywie ujemne dopasowania14. Podświetlone obszary właściwości piku szablonu w Rysunek 1B pokazuje informacje o referencyjnym ciągu widma MS i funkcji ograniczenia qCLIC (tj. (Match("
Stosując szablon do wszystkich chromatogramów zestawu, można napotkać trudne sytuacje, jak w przypadku częściowego niewspółosiowości wzorców. Może to być spowodowane niespójnościami temperatury pieca, niestabilnością przepływu/ciśnienia gazu nośnego lub ręczną interwencją w system, jak w przypadku substytucji kolumny lub wymiany pętli kapilarnej modulatora14,28. Rysunek 3 pokazuje sytuację częściowego niewspółosiowości między docelowym szablonem a rzeczywistym chromatogramem. Aby zminimalizować niewspółosiowość, interaktywne transformacje szablonów (Rysunek 3, panel sterowania) mogą zmieniać położenie szczytów szablonu w celu lepszego dopasowania. Po zmianie położenia szablon można dopasować w celu nawiązania korespondencji. W tym przykładzie szczyty szablonu (Rysunek 3, krok 12) są poprawnie zgodne z rzeczywistym wzorcem 2D. W przypadku poważnych niewspółosiowości, które nie zostały tutaj omówione, powtórzenie działań dopasuj-przekształcenie-aktualizacja może iteracyjnie dostosować pozycję pików szablonu do rzeczywistego wzorca szczytu12,13,14.
Tutaj docelowe piki (tj. znane anality) dostarczają około 40% wyniku chromatograficznego (196 docelowych pików średnio około 500 wykrywalnych). Pozostałe 60% związków, wraz z informacjami, które przynoszą, nie jest brane pod uwagę w ukierunkowanej analizie. Aby badanie było naprawdę wszechstronne, należy również ustanowić konsekwentne wyrównanie krzyżowe niedocelowych pików 2D. Pierwsza aplikacja, w której dopasowanie szablonu zostało rozszerzone na wszystkie wykrywalne anality, zajmowała się złożonym objętością palonej kawy7. Proces ten jest zautomatyzowany za pomocą oprogramowania (np. Investigator), pokazanego tutaj w krokach 14–15.
W tym procesie, wstępnie wybrane obrazy należące do badanego zestawu próbek (20 próbek) są używane do definiowania wiarygodnych szczytów poprzez krzyżowe dopasowanie wszystkich wzorców obrazów29. Następnie budowany jest chromatogram kompozytowy, na podstawie którego można zidentyfikować wiarygodne piki UT i regiony pików (tj. ślad pików 2D) w tak zwanym szablonie cech17.
Dla analiz uzyskanych przy 70 eV, proces ten określił 144 wiarygodne piki ze złagodzoną niezawodnością29, z których 76 należy do listy docelowych szczytów. Opierając się na tych 144 wiarygodnych pikach, proces wyrównuje wszystkie chromatogramy zgodnie ze średnimi czasami retencji wiarygodnych pików, a następnie łączy je w celu utworzenia chromatogramu kompozytowego. Rysunek 4 pokazuje listę wszystkich próbek oznaczonych zgodnie z regionem produkcji oleju (po lewej) oraz listę wiarygodnych pików/objętości plam w każdej próbce (po prawej).
Szablon nieukierunkowanej cechy składa się z 2pików D z analitów wykrytych w chromatogramie kompozytowym, pokazanych w Rysunek 5A, które są dopasowane przez szablon niezawodnych pików (n = 168 - czerwone kółka dla docelowych pików i zielone kółka dla niedocelowych pików). Widma masowe pików kompozytu, jak również ich czasy retencji, są rejestrowane w szablonie cech, jak pokazano dla octanu (Z)-3-heksenolu w powiększonym obszarze. Obszary pików są pokazane w Rysunek 5B jako grafika w kolorze czerwonym; zamiast tego są one określone przez kontury wszystkich 2pików D wykrytych w chromatogramie złożonym (n = 3578).
Gdy nienadzorowane rozpoznawanie wzorców przez Analizę Głównych Składowych jest stosowane do docelowego rozkładu pików w 20 analizowanych próbkach, oleje sycylijskie i toskańskie grupują się oddzielnie, co sugeruje, że warunki glebowo-klimatyczne i terroir wpływają na względne rozpowszechnienie substancji lotnych. Wyniki są pokazane w Rysunek 6A, a wyniki PCA z niezawodnego rozkładu pików są pokazane w Rysunek 6B. Oba podejścia potwierdzają krzyżowo, że oleje z różnych obszarów geograficznych mają różne, podczas gdy spójne sygnatury chemiczne, niezależnie od tego, czy są to związki docelowe, czy nieukierunkowane, lub jedno i drugie, są mapowane.
Wreszcie, oprogramowanie umożliwia szybkie i skuteczne ponowne wyrównanie wzorców między równoległymi kanałami detekcji. W tym zastosowaniu proponowane jest ponowne wyrównanie dla tandemowych sygnałów jonizacji. Źródło jonów MS multipleksuje się między dwiema energiami jonizacji (tj. 70 i 12 eV) z częstotliwością akwizycji 50 Hz na kanał30. Dwa wynikowe wzorce chromatograficzne są ściśle dopasowane, podczas gdy dane spektralne (tj. sygnatury spektralne i odpowiedzi) dostarczają uzupełniających się informacji o różnych zakresach dynamicznych odpowiedzi26,27. Wyrównane wzorce umożliwiają wyodrębnianie cech (2D pików i regionów pików) z jednoznacznymi identyfikatorami (tj. nazwami chemicznymi dla docelowych pików i unikalną numeracją # dla niedocelowych pików i regionów pików).
Dopasowanie szablonu pozwala na efektywne wyrównanie krzyżowe. W tej sytuacji nie ma dużej niewspółosiowości, ale ograniczenia MS muszą zostać złagodzone, aby umożliwić dopasowania szczytów UT. Z drugiej strony polecane regiony szczytowe UT, które nie mają ograniczeń MS, są natychmiast dopasowywane bez żadnych fałszywie ujemnych dopasowań. Rysunek 5C pokazuje powiększony obszar chromatogramu 12 eV, gdzie dopasowany jest szablon cech zbudowany z danych 70 eV. Wiarygodne piki UT są dodatnio dopasowane ze względu na obniżone ograniczenia qCLIC (np. próg DMF na poziomie 600). Należy zauważyć, że przy 12 eV wykrywa się mniej pików ze względu na ograniczoną fragmentację wywołaną niską energią jonizacji.

Rysunek 1: Dwuwymiarowy wykres konturowy i szablon docelowy. (A) Wykres konturowy frakcji lotnej oliwy z oliwek z pierwszego tłoczenia z Toskanii. Uporządkowane wzory szeregów i klas homologów są wyróżnione różnymi kolorami i liniami: liniowo nasycone węglowodory (linia i kontury 2D), węglowodory nienasycone (żółte), liniowo nasycone aldehydy (niebieskie), aldehydy jednonienasycone (czerwone), aldehydy wielonienasycone (łosoś), pierwszorzędowe alkohole (zielony) i krótkołańcuchowe kwasy tłuszczowe (cyjano). (B) Nałożony na siebie ukierunkowany szablon znanych analitów (czerwone kółka) z liniami łączącymi Standardy Wewnętrzne (IS). Panele pokazują metadane właściwości szczytu/obiektu blob 2D (Decanal) lub właściwości piku szablonu. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Wyszukiwanie Apex MS. Dane wyjściowe wyszukiwania wierzchołka MS dla obiektu blob 5. Lista wpisów bazy danych o najwyższym dopasowaniu podobieństwa i powiązanych metadanych dostępnych w bibliotece. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Wyrównanie szablonu. Przepływ pracy ilustrujący kroki, które umożliwiają ponowne wyrównanie szablonu przez transformację. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Interfejs GC Investigator. Panel badacza ze wszystkimi wybranymi obrazami oznaczonymi zgodnie z regionem produkcji oleju (po lewej) oraz listą wiarygodnych pików/objętości plam w każdej próbce (po prawej). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Szablon ukierunkowany i UT. (A) Wiarygodne wartości szczytowe wynikające z automatycznego przetwarzania w kroku 11; Czerwone kółka odpowiadają znanym analitom, podczas gdy zielone kółka są nieznane. W nałożonym panelu wyświetlane są właściwości obiektu szablonu dla (Z)-3-hexenal. (B) Powiększony obszar, który pokazuje piki UT (czerwone i zielone kółka) oraz obszary pików (czerwona grafika) szablonu UT dopasowane do próbki oleju uzyskanej przy energii jonizacji 70 eV. (C) Matryca UT dopasowana do próbki oleju uzyskanej przy energii jonizacji 12 eV. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Wykresy załadunku PCA. Pokazują one naturalną konformację próbek (oleje z Toskanii i Sycylii) wynikającą z (A) docelowego rozkładu pików lub (B) rozkładu pików UT. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Pliki uzupełniające. Kliknij tutaj, aby pobrać te pliki.