Analiza interakcji leków (wchłaniania, metabolizmu i efektów) na poziomie pojedynczych komórek jest niezbędna do wykrycia ukrytych lub opornych na leki subpopulacji, a także do zrozumienia wpływu heterogeniczności komórkowej. W niniejszym protokole wykorzystano dwie komplementarne techniki do pomiaru wspomnianych interakcji w pojedynczych komórkach: spektroskopię Ramanowską oraz MS. Spektrometria Ramanowska pozwala na szybką identyfikację komórek pod wpływem leków w oparciu o spektralne biomarkery odpowiedzi na lek. MS jest wykorzystywana do monitorowania wchłaniania i metabolizmu leku w sposób selektywny i półilościowy. Komórki były najpierw przesiewane za pomocą spektroskopii Ramanowskiej, a następnie pobierane indywidualnie do analizy metodą MS.
Porównawcza analiza średniego widma dla każdego stanu (z lekiem i bez leku) przedstawiona jest na Rysunku 2. Uśrednione widma dwóch stanów wyraźnie różnią się w kilku pikach, które zostały wcześniej zidentyfikowane i przypisane do konkretnych związków molekularnych2. W szczególności silne różnice występują w pikach przy 1000 cm- (przypisanych do związków aromatycznych, takich jak fenyloalanina i tyrozyna). Istotność różnic statystycznych powinna zostać oceniona za pomocą dalszych analiz wielowymiarowych.
Zbiór danych wykorzystano następnie do wytrenowania modelu PLS (kroki 4.5-4.8) mającego na celu rozróżnienie dwóch wariantów traktowania komórek (z lekiem: n = 290, bez leku: n = 115). Zdolność prognostyczna w klasyfikacji komórek hodowanych w obecności tamoksyfenu osiągnęła 100% czułości i 72% swoistości w danych testowych (nieznanych dla modelu wytrenowanego z walidacją krzyżową). Czułość jest miarą wyników prawdziwie dodatnich, które zostały poprawnie zidentyfikowane przez model, natomiast swoistość jest miarą rzeczywistych wyników ujemnych zidentyfikowanych przez model. Alternatywne modele, takie jak SVM, LDA i sieci neuronowe, mogą dawać podobne lub lepsze wyniki, choć w niniejszym badaniu nie przeprowadzono kompleksowego porównania.
Na podstawie modelu PLS obliczono wartości VIP, które reprezentują znaczenie długości fal (przesunięć ramanowskich) w rozróżnianiu warunków eksperymentalnych (Rysunek 3). Co istotne, najwyższe piki profilów VIP odpowiadały piki ramanowskim, dla których zaobserwowano wyraźne różnice między dwiema grupami traktowanymi. Potwierdziło to specyficzne różnice molekularne między komórkami traktowanymi a nietraktowanymi. W konsekwencji badacze mogą zidentyfikować potencjalne biomarkery spektralne, które odzwierciedlają odpowiedź pojedynczych komórek na leczenie farmakologiczne. Biomarkery te mogą być poddane dalszym badaniom w celu weryfikacji ich znaczenia biologicznego oraz możliwości generalizacji wyników w różnych warunkach i liniach komórkowych.
System spektrometrii mas pojedynczych żywych komórek (LSC-MS) pozwolił na wykrycie zarówno leku, jak i jego metabolitów w pojedynczych, traktowanych lekiem komórkach HepG2, które wcześniej poddano pomiarom za pomocą spektroskopii Ramanowskiej. Ponadto zastosowano tandemową spektrometrię mas (MS/MS) w celu potwierdzenia struktury obu cząsteczek. Po pozytywnej identyfikacji zmierzono względną zawartość leku oraz jego metabolitów w każdej komórce i porównano je z pikami tła w komórkach nietraktowanych. Zaobserwowano silną zmienność zawartości tamoksyfenu, a zjawisko to było jeszcze bardziej wyraźne w przypadku jego metabolitu, 4-OHT (Rysunek 4). Zbadano również zależność między zawartością tamoksyfenu a jego metabolitami, wykazując istotną dodatnią korelację między nimi (r = 0,54, p = 0,0001, n = 31).

Rycina 1: System do pobierania pojedynczych komórek zamontowany na stoliku mikroskopu. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Uśrednione widmo komórek leczonych lekiem (z tamoksyfenem: n = 295) oraz komórek nieleczonych (bez tamoksyfenu: n = 115). Pików Ramanowskich można zidentyfikować na podstawie literatury. Większość silnych różnic spektralnych jest istotna statystycznie (ANOVA, p ≤ 0.5), jak opisano wcześniej4. Rysunek ten został zmodyfikowany z wcześniejszej publikacji4. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3: Wyniki VIP wyekstrahowane z predykcyjnego modelu PLS. Wyniki VIP odzwierciedlają długości fal, które przyczyniają się do rozróżnienia dwóch klas w modelu. Większość pików odpowiada konkretnym cząsteczkom, które są obserwowane jako biomarkery spektralne efektów działania leków w komórkach poddanych leczeniu. Rysunek został zmodyfikowany na podstawie wcześniejszej publikacji4. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 4: Rozkład obfitości tamoksyfenu i jego metabolitu. Rozkład obfitości tamoksyfenu oraz jego metabolitu, 4-OHT (mierzonej na poziomie pojedynczych komórek) w porównaniu do pików endogennych w komórkach nietraktowanych (kontrola). Rysunek ten został zmodyfikowany z wcześniejszej publikacji4. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.