W tym filmie przedstawiamy pełny eksperyment śledzenia pojedynczych cząsteczek z wykorzystaniem kropek kwantowych skierowanych przeciwko konkretnemu receptorowi błonowemu. Głównym celem tego eksperymentu jest rozróżnienie różnych typów zachowań dyfuzji molekularnej mierzonych w błonie plazmatycznej żywych komórek. Rzeczywiście, ruchy molekularne występujące w błonie mogą zazwyczaj odbiegać od dyfuzji Browna, będąc na przykład ukierunkowane liniowo lub ograniczone w obrębie nanodomen26-29. Naszym celem jest jednoczesne śledzenie tak wielu receptorów, jak to technicznie możliwe, aby uzyskać obraz różnorodności dynamiki zachodzącej w błonie żywej komórki. Oczekuje się, że w ostateczności pozwoli to na odczytanie mechanizmów regulujących sygnalizację receptorów powierzchni komórki.
1. Hodowla komórkowa
- Przygotować próbkę komórkową: wykorzystać adhezyjne komórki COS-7, które endogennie wykazują ekspresję receptora epidermalnego czynnika wzrostu (EGFR)30. Podczas pracy z żywymi komórkami bez antybiotyków należy upewnić się, że nie występuje utajone, niewykrywalne zanieczyszczenie, stosując przez cały czas przygotowania odpowiednie techniki sterylne.
- Hodować komórki w pożywce kompletnej (DMEM z 10% serum cielęcego, 1% glutaminy, 1% HEPES i 1% pirogronianu sodu, patrz tabela specyficznych odczynników poniżej), w temperaturze 37 °C przy 7% CO2, dbając o to, aby znajdowały się one w fazie wzrostu wykładniczego przed osiągnięciem pełnej konfluencji, przed wysianiem ich do Lab-Tek.
- W dniu poprzedzającym eksperyment wysiać 5 000 komórek/studnię w komorach 8-studzienkowych (Lab-Tek) i inkubować przez noc. Odliczanie komórek zapewnia stałą gęstość komórkową, a tym samym powtarzalny stosunek kropek kwantowych na jedną komórkę.
2. znakowanie komórek
Przygotuj kropki kwantowe z odpowiednią powłoką. Kropki kwantowe są fluorescencyjnymi nanocząstkami zbudowanymi z półprzewodników. Nanocząstki te budzą duże zainteresowanie, ponieważ są bardzo jasne i fotostabilne w porównaniu z klasycznymi sondami fluorescencyjnymi31,32, co pozwala na osiągnięcie odpowiedniego stosunku sygnału do szumu (SNR) w obrazowaniu pojedynczych cząsteczek.
- Przed rozpoczęciem eksperymentu należy przygotować fragmenty Fab przeciwko EGFR z linii komórkowej hybridoma (mAb 108, ATCC HB-9764) poprzez trawienie papainą, zgodnie z wcześniejszym opisem21.
- Należy skoniugować Fab z biotyną zgodnie z instrukcjami producenta (EZ-link sulfo-NHS-LC-biotinylation Kit; Thermo Scientific, Rys. 1A).
- Należy zastosować kropki kwantowe funkcjonalizowane streptawidyną (Invitrogen), emitujące przy 605 nm (optymalna długość fali emisji dla jasności, detekcji i oddzielenia od autofluorescencji komórkowej).
- Należy przygotować roztwór znakujący w pełnej pożywce (patrz tabela specyficznych odczynników), aby nasycić streptawidyny obecne wokół kropek kwantowych, a także zapobiec agregacji i niespecyficznemu wiązaniu do komórek oraz do szkiełka podstawowego.
- Należy przygotować dwa roztwory pośrednie: jeden z kropkami kwantowymi-streptawidyną i drugi z biotynylowanym Fab, każdy w stężeniu 20 nM.
- Należy zmieszać równe objętości tych dwóch roztworów: zmieszać kropki kwantowe z Fab w stosunku 1:1, aby uzyskać końcowe stężenie robocze kompleksu Fab:kropki kwantowe wynoszące 10 nM. Zastosowanie stosunku ekimolarnego sprzyja tworzeniu kompleksów monofunkcjonalnych, składających się z jednej kropki kwantowej i jednego biotynylowanego fragmentu Fab. Sprzyja to znakowaniu monowalentnemu, ograniczając obserwację artefaktów wynikających z sieciowania receptorów33,34.
- Inkubować przez 15 minut w 25 °C, używając wytrząsarki przy 1 200 obrotach na minutę, aby zapobiec agregacji. Mieszanina jest następnie gotowa do użycia na żywych komórkach.
- Inkubować komórki w 100 μl mieszaniny przez 5 minut w 37 °C z 7% CO2.
- Przemyć komórki nieautofluorescencyjną pożywką do obrazowania (bufor HBSS, 1% HEPES, patrz tabela specyficznych odczynników) uprzednio ogrzaną do 37 °C.
- Ostrożnie usunąć nadmiar znaczników, aby zapobiec niepotrzebnemu pogorszeniu SNR przez niezwiązane, poza ogniskową kropki kwantowe, poprzez dokładne przemycie każdej studzienki pożywką do obrazowania, zazwyczaj 5 razy w temperaturze pokojowej, z co najmniej 5-minutową przerwą przed ostatnim płukaniem.
3. Układ optyczny
Układ wideomikroskopowy składa się z czterech głównych części:
- Mikroskop odwrócony z odpowiednim kostkami fluorescencyjnymi (filtry wzbudzenia FF01-457/50, dichroiczny FF495-Di02 oraz emisyjny FF01-617/73, Semrock) oraz obiektywem imersyjnym olejowym 100x o wysokiej aperturze numerycznej (1.3 lub 1.49).
- Lampa rtęciowa 100 W (połączona z mikroskopem za pomocą światłowodu, aby uniknąć zakłóceń w termoregulacji).
- Kamera EMCCD o wysokiej czułości i rozdzielczości 512 x 512 pikseli w celu uzyskania wystarczającego stosunku sygnału do szumu (SNR).
- Inkubator do utrzymywania próbek biologicznych w temperaturze 37 °C podczas eksperymentu.
4. Akwizycja
- Wybierz odizolowane i dobrze rozprzestrzenione komórki. Sprzyja to powstawaniu ładnych, płaskich lamellipodiów, które najlepiej nadają się do obserwacji ruchu planarnego cząsteczek błonowych.
- Oceń stan fizjologiczny komórki na podstawie jej wyglądu zarówno w świetle przechodzącym, jak i fluorescencyjnym: intensywny transport pęcherzykowy, brak oznak nekrozy lub apoptozy, niska autofluorescencja oraz średnio-silne znakowanie kropkami kwantowymi (zazwyczaj do 1 000 kropek kwantowych na komórkę, patrz Ryc. 1B, C).
- Wybierz wysoką gęstość znakowania, co jest krytyczne dla jednoczesnego monitorowania jak największej liczby receptorów. Jest to w rzeczywistości ograniczone zarówno przez parametry fizjologiczne (ekspresja powierzchniowa, dostępność epitopu, ruch lateralny), jak i algorytmiczne (niezachodzące na siebie funkcje rozproszenia punktu (PSF), nawet przy uwzględnieniu rozmycia wynikającego z ruchu, aby umożliwić prawidłową detekcję i ponowne połączenie).
- Dla każdej komórki wykonaj najpierw obraz w polu jasnym (preferencyjnie z użyciem DIC, jeśli jest dostępny), co pozwoli w dalszej kolejności sprawdzić wygląd komórki i granice przestrzenne lamellipodiów.
- Zapisz ten obraz, używając tej samej nazwy co stos wideo (na przykład cell1.tif i cell1.stk), w podfolderze o nazwie dic, ponieważ przy takich konwencjach algorytm może automatycznie odnaleźć pasujący obraz dla każdego stosu.
- Nagraj od 1 do 3 filmów na komórkę, w sposób ciągły, zazwyczaj z częstotliwością 36 ms, co jest najszybszym tempem osiągalnym w pełnej klatce przy użyciu tej kamery. Można jednak rejestrować z wyższymi częstotliwościami, do 1 ms, korzystając z dedykowanych matryc CCD o zwiększonej czułości i/lub mniejszej liczbie pikseli. Technologia transferu klatek zapewnia pomijalne opóźnienie między klatkami.
- Wzmocnienie mnożnika elektronowego powinno być zawsze ustawione na maksimum (tuż poniżej saturacji, jeśli dotyczy), aby umożliwić osiągnięcie czułości na poziomie pojedynczej cząsteczki, z wystarczającym SNR, co najmniej powyżej 20 dB (dla efektywnej detekcji piku1), zazwyczaj około 25-30 dB.
- Zazwyczaj rejestruje się 300 klatek na film, ponieważ trajektorie są rekonstruowane średnio z około 100 klatek, co jest głównie ograniczone przez długie zdarzenia mrugania. Częstotliwość i długość filmu mogą być dostosowane do konkretnego pomiaru, który mógłby wymagać na przykład dłuższych śladów.
5. Analiza MTT
- Wybierz ścieżkę do katalogu zawierającego pliki wideo, aby ocenić dany zestaw danych za pomocą programu Matlab lub Octave.
- Aby rozpocząć w pełni zautomatyzowaną analizę1,35, wpisz polecenie detect_reconnex23 w programie Octave lub MTT23i w programie Matlab. Program ten oznacza „wersję 2.3 z interfejsem użytkownika” i jest dostępny do pobrania wraz z poprzednią wersją 2.2. Najpierw wyświetla on interfejs graficzny z listą wszystkich użytych parametrów, jak pokazano na Rys. 2.
6. Reprezentatywne wyniki
MTT automatycznie analizuje każde nagranie wideo, aby dostarczyć ścieżki wykrytych i oszacowanych celów, uzupełnione o dalsze analizy, takie jak wykrywanie uwięzienia. Pozwala to w efekcie na nałożenie ścieżek na obrazy komórek (Rys. 1C i 3).
Opis MTT
Podstawowa analiza MTT jest przeprowadzana dla każdej klatki, obejmując 3 główne zadania (Ryc. 3):
- Detekcja obecności lub braku celu (kropki kwantowej) w obrębie przesuwającego się podregionu, kolejno centrowanego wokół każdego piksela, gdzie porównuje się dwie hipotezy: obecność sygnału, z PSF modelowanym jako dwuwymiarowy pik Gaussa, lub obecność samego szumu, przy zastosowaniu progu zapewniającego odpowiednio niski poziom fałszywych alarmów, z liczbą mniej niż jednej błędnej detekcji na klatkę. Prowadzi to do uzyskania mapy prawdopodobieństwa detekcji. Każde lokalne maksimum jest traktowane jako potencjalny cel, pod warunkiem, że jego poziom prawdopodobieństwa jest wyższy od wartości minimalnej, określonej zgodnie z pożądanym prawdopodobieństwem fałszywego alarmu (PFA). Limit ten jest ustawiony na poziomie na tyle niskim, aby skutecznie odróżnić sygnały od szumu, minimalizując liczbę błędnych detekcji (domyślnie PFA rzędu 10-6, aby zapewnić mniej niż jeden błąd w obrazach o wymiarach 512 x 512 pikseli), przy jednoczesnym zachowaniu wystarczająco wysokiego prawdopodobieństwa detekcji, osiągającego teoretycznie przewidziane optimum1. Należy zauważyć, że wymóg zastosowania podregionu sprawia, iż brzegi obrazu (3 piksele dla domyślnego okna o wymiarach 7 x 7 pikseli) nie mogą zostać poddane ocenie.
- Estymacja odpowiednich parametrów dla każdego wykrytego celu, takich jak pozycja podpikselowa i intensywność sygnału. Dla każdego wykrytego celu przeprowadza się następnie dopasowanie metodą najmniejszych kwadratów Gaussa-Newtona w celu oszacowania pozycji, szerokości i wysokości wykrytej funkcji Gaussa. Pozwala to w szczególności na wyznaczenie pozycji podpikselowej barwnika (dokładność od 10 do 20 nm dla typowych wartości SNR oraz barwników nieruchomych lub powoli dyfundujących, wzrastająca do ~100 nm dla barwników dyfundujących z prędkością 0,1 μm2/s).
- Rekoneksja nowych celów ze śladami zbudowanymi na podstawie poprzednich klatek. Zbiór nowych celów jest dopasowywany do zbioru poprzednich śladów. W tym celu, aby przypisać każdy cel do śladu (jeśli jest to możliwe), wykorzystuje się wszystkie dostępne informacje statystyczne uzyskane na etapie detekcji – nie tylko pozycję, ale także intensywność, szerokość, mruganie (blinking) oraz statystyki powiązane. Zatem cele nie są przypisywane jedynie do najbliższego śladu: w przypadku krzyżujących się śladów brane będą pod uwagę intensywność, prędkość, szerokość i mruganie. Pozwala to uzyskać statystycznie optymalny wynik rekoneksji. Strategia ta pozwala, w miarę możliwości, uniknąć faworyzowania najbliższych sąsiadów podczas procesu łączenia.
Wykryte piki mogą zostać odrzucone a posteriori, jeśli ich estymacja lub ponowne połączenie okażą się nieudane. Specjalny test obsługuje wykrywanie nowych pików, które inicjowałyby nowe ślady. Test ten wykorzystuje bardziej rygorystyczny PFA (10-7), ponieważ ponowne połączenie piku ze śladem może być interpretowane de facto jako walidacja jego istotności (kryterium to z definicji nie ma zastosowania w przypadku nowych pików).
Analiza trajektorii
W następnym kroku oceniono możliwe przejściowe ograniczenie (transient confinement) za pomocą funkcji odwrotnie proporcjonalnej do dyfuzji lokalnej24-29. Zastosowanie wartości progowej pozwala na zdefiniowanie epizodów ograniczenia lub ich braku. Poprzez iterację tych zdarzeń we wszystkich śladach można zmapować dynamikę błony pod kątem zdarzeń przejściowego ograniczenia/spowolnienia. Można to przedstawić alternatywnie, wykorzystując binarne lub dyskretne wartości tego indeksu ograniczenia.
Domyślnie MTT automatycznie wykonuje te zadania, zapisując 8 parametrów piku w pliku tekstowym: numer klatki, pozycję i oraz j, intensywność sygnału, promień, przesunięcie (offset) oraz mrugnięcie (blink), dla każdej klatki wideo (grupa 7 wierszy) i śladu (kolumna). Te parametry wyjściowe można ponownie załadować do programu Matlab lub Octave za pomocą skryptu fread_data_spt w celu dalszej analizy, np. śladów lub intensywności sygnału, co zilustrowano w skrypcie MTT_example w aneksie.
Dalsze analizy pozwalają na mapowanie śladów w każdej komórce (Ryc. 1C i 3) oraz na przedstawienie rozkładów histogramów dla istotnych parametrów (takich jak intensywności szczytów, SNR lub lokalne wartości dyfuzji). Dla każdego pliku średnia i odchylenie standardowe każdego parametru są zapisywane w pliku tekstowym wraz z obrazem histogramów. Rozkłady logarytmiczne, np. dla przesunięć kwadratowych r2, prowadzą do wyznaczenia średniej geometrycznej. Współczynnik dyfuzji D jest obliczany na podstawie dopasowania liniowego dla pięciu pierwszych punktów krzywej MSD. Wartości te zapewniają przegląd eksperymentu obejmującego na przykład kinetykę reakcji komórkowych lub zabiegi farmaceutyczne/enzymatyczne wpływające na organizację błony. Ponieważ MTT jest kodem otwartoźródłowym, aspekt ten można łatwo dostosować do dowolnych specjalistycznych badań.

Rycina 1. Monitorowanie dynamiki receptorów błonowych za pomocą MTT. (A) Komponenty błony, takie jak EGFR, są znakowane kropkami kwantowymi sprzężonymi z biotynowanymi fragmentami Fab (rysunek schematyczny z przybliżoną skalą każdej cząsteczki). (B) Typowy obraz fluorescencyjny uzyskany z żywej komórki COS-7, z czasem ekspozycji 36 ms, przedstawiający piki ograniczone dyfrakcyjnie, odpowiadające poszczególnym znakowanym receptorom. (C) Wynik analizy MTT przedstawiający zrekonstruowane trajektorie receptorów, nałożone na obraz komórki w jasnym polu.

Rysunek 2. Parametry wejściowe MTT. Uruchomienie MTT23i otwiera graficzny interfejs użytkownika wyświetlający wszystkie parametry wejściowe, nazwy i wartości domyślne, zgodnie z opisem w naszej poprzedniej publikacji1. W algorytmie parametry przestrzenne i czasowe (okna wyszukiwania, promień piku, maksymalna dyfuzja i migotanie) podane są w bezwymiarowych jednostkach standardowych: pikselach i klatkach. Kalibracje mogą być zastosowane a posteriori w celu przeliczenia wyników wyjściowych. Wartości domyślne, odpowiadające Cascade 512BFT z powiększeniem 100x, to rozmiar piksela: 156 nm/pxl oraz opóźnienie klatki: 36 ms/frame.
Badacze powinni zoptymalizować kilka krytycznych parametrów, takich jak maksymalny oczekiwany współczynnik dyfuzji („Diff max”) oraz maksymalny czas zaniku migotania („Toff”). Te dwa limity przestrzenne i czasowe są niemal jedynymi parametrami, które wymagają ponownego rozważenia dla danej warunki eksperymentalnej; pozostałe ustawione są na stabilne wartości domyślne. Na przykład, liczba fałszywych alarmów jest ustawiona bezpośrednio tak, aby zapewnić mniej niż jeden błąd na milion pikseli, a zatem mniej niż jeden na klatkę, co w większości przypadków jest satysfakcjonujące. Wszystkie parametry są zapisywane w pliku tekstowym w folderze wyjściowym, co pozwala użytkownikom na późniejszą weryfikację ustawień wykorzystanych do analizy.

Rycina 3. Główne etapy analizy MTT. Począwszy od eksperymentalnego stosu obrazów fluorescencyjnych, piki są sekwencyjnie i automatycznie wykrywane, szacowane za pomocą dopasowania Gaussa i ponownie łączone w ramach klatek (pierwszy zakres operacji, górna część schematu blokowego). Ślady mogą być następnie analizowane pod kątem domniemanego uwięzienia, co na koniec prowadzi np. do stworzenia mapy dynamicznej odpowiednich deskryptorów (drugi zakres operacji, dolna część).

Rycina 4. Walencyjność znakowania nie wpływa na MTT. Aby ocenić możliwe obciążenie wprowadzone przez artefakt wielowartościowego znakowania, przeprowadzono analizę MTT w celu śledzenia endogennego EGFR znakowanego za pomocą dwóch różnych schematów, aby wygenerować i przeanalizować mapy trajektorii. (A) Receptory zostały znakowane biotynيلowanymi fragmentami Fab oraz kropkami kwantowymi quantum-dots605-streptawidyną, zgodnie z opisem w protokole. W tym przypadku wielowartościowość kropek kwantowych i streptawidyn może prowadzić do sprzężenia kilku receptorów z jednym barwnikiem. (B) Receptory zostały znakowane fragmentami Fab bezpośrednio sprzężonymi z organicznym barwnikiem Atto647N. W tym przypadku jeden fragment Fab, a zatem jeden receptor, może być sprzężony z więcej niż jednym barwnikiem. (C) Krzywa średniego kwadratu przemieszczenia (MSD) została obliczona dla wszystkich śladów w każdej komórce, znakowanej albo kropkami kwantowymi, albo barwnikami Atto (odpowiednio wykresy lewy i prawy). Współczynniki dyfuzji obliczono za pomocą dopasowania liniowego dla pięciu pierwszych punktów MSD (czerwona przerywana linia). Każdy schemat znakowania doprowadził do podobnych wartości dyfuzji (wykres środkowy). Qdot: quantum-dots605 (n = 5 komórek), Atto: Atto647N (n = 7 komórek), ns: brak istotności (wartość p testu t Studenta > 0,05).