Ogólny schemat metody przedstawiono na Rysunku 1. Podsumowuje on główne etapy opisane w sekcji protokołu: od przygotowania próbki (Rysunek 1A), przez obrazowanie nanostruktur wewnątrzkomórkowych w technice time-lapse (Rysunek 1B), analizę fluktuacji w celu obliczenia serii czasoprzestrzennych funkcji korelacji (Rysunek 1C), aż po dopasowanie danych w celu wyznaczenia średnich właściwości strukturalnych i dynamicznych badanego obiektu (Rysunek 1D).
Kluczowym parametrem jest rozdzielczość czasowa przyjęta do obrazowania interesującego obiektu subkomórkowego. Ta wartość eksperymentalna określi próg czasowy, przy którym mierzony będzie minimalny średni przemieszczenie badanych obiektów. Jednakże preferowanym warunkiem jest ustawienie takiej rozdzielczości czasowej obrazowania, przy której badany obiekt wydaje się „nieruchomy” w obrębie przechwyconej klatki, t. j. wykazuje charakterystyczny rozmiar, który średnio nie ulega deformacji ze względu na prędkość obrazowania. Jest to technicznie możliwe, jeśli badany obiekt jest strukturą subkomórkową lub organellą otoczoną błoną (tak jak w tym przypadku). Typowo struktury subkomórkowe wykazują lokalne współczynniki dyfuzji (D, µm2/s, patrz Tabela 1), które są o kilka rzędów wielkości niższe niż współczynniki dla pojedynczych izolowanych cząsteczek w cytoplazmie (np. GFP21). Walidację można przeprowadzić poprzez sztuczne unieruchomienie badanej organelli (np. za pomocą fiksacji chemicznej). W rzeczywistości warunek ten może służyć jako punkt odniesienia do wyznaczenia rzeczywistego rozmiaru organelli w zastosowanych warunkach eksperymentalnych (np. długość fali wzbudzenia, rozmiar piksela, obiektyw).
W tym przypadku, mimo że w procedurze wykorzystano lizosomy jako organelle testowe, wyniki są ważne niezależnie od badanej struktury. Rysunek 2A przedstawia obraz utrwalonego (tzn. nieruchomego) lizosomu wraz z akwizycją wykonaną na żywych komórkach przy odpowiedniej rozdzielczości czasowej (zazwyczaj poniżej 100 ms/klatkę; np. 65 ms/klatkę w przykładzie na Rysunku 2B) oraz akwizycją wykonaną celowo przy bardzo niskiej rozdzielczości czasowej (np. 10 s/klatkę w przykładzie na Rysunku 2C). Dla każdego warunku rozmiar dyfundującego obiektu wyznacza się w następujący sposób: i) za pomocą narzędzia linii w oprogramowaniu ImageJ wyznacza się profil intensywności plamki; ii) profil intensywności jest wykreślany i interpolowany funkcją Gaussa w celu obliczenia wartości pełnej szerokości połówkowej (FWHM), która z kolei służy jako szacunkowa średnica plamki (Rysunek 2D). Zgodnie z oczekiwaniami, co pokazano na wykresie na Rysunku 2E, akwizycja przy bardzo wysokiej rozdzielczości czasowej (tzn. 65 ms/klatkę) daje średni rozmiar struktury zbliżony do uzyskanego w próbce utrwalonej, zarówno przy użyciu opisanego powyżej standardowego narzędzia, jak i poprzez wyznaczenie punktu przecięcia osi y dla iMSD. Z kolei akwizycja przy niskiej prędkości skutkuje wzrostem pozornego rozmiaru struktury ze względu na naturalną dynamikę struktury podczas obrazowania. W nieoptymalnych warunkach eksperymentalnych wyodrębnione informacje strukturalne/dynamiczne nie odzwierciedlają wiernie właściwości wewnętrznych badanego obiektu.
Po wyborze głównych parametrów eksperymentalnych można przygotować zestawy danych dla docelowych struktur wewnątrzkomórkowych. W przypadku makropinosomów, po 20 min inkubacji komórek z dekstranami 70 kDa, serie czasowe znakowanych struktur wewnątrzkomórkowych pozyskiwano w różnych punktach czasowych po traktowaniu, od 30 min do około 180 min. Co ciekawe, podczas transportu wykryto stopniową zmianę właściwości strukturalnych i dynamicznych makropinosomów (Rysunek 3; rozkłady σ02, Dm, α oraz N dla makropinosomów przedstawiono na wykresach po lewej stronie). Podczas gdy nie stwierdzono wyraźnych zmian w lokalnej dyfuzyjności (Dm) makropinosomów w trakcie transportu, zarówno charakterystyczny rozmiar (σ02), jak i ogólny tryb ruchu (α) zmieniały się w czasie.
Szczególną uwagę zwraca spadek średniego rozmiaru makropinosomów podczas ich transportu (Rycyna 3A, lewa strona), wraz z towarzyszącym temu wzrostem charakteru subdyfuzyjnego ich ruchu (tj. określonego jako spadek wartości α, Rycyna 3C, panel lewy). Dodatkowo, z każdego zarejestrowanego obrazu wyznaczono liczbę makropinosomów: wyniki przedstawione na Rycynie 3D (panel lewy) wyraźnie wykazują wzrost liczby makropinosomów w czasie. Wszystkie te wyniki są zgodne z oczekiwaniami, ponieważ makropinosomy znakowane dekstranem powinny powstawać jako izolowane, duże pęcherzyki otoczone błoną w obrębie błony plazmatycznej (zdolne również do przemieszczania się wzdłuż elementów cytoszkieletu), a następnie stopniowo łączyć się ze szlakiem endo-lizosomalnym, składającym się z dużej populacji mniejszych struktur dyfundujących w sposób losowy.
Zgodnie z wcześniejszymi przewidywaniami, wyniki dla makropinosomów są przeciwstawne do pomiarów wykonanych dla ISG (Rysunek 3, prawa kolumna). Granule insulinowe nie wykazują trendu czasowego parametrów strukturalnych/dynamicznych pochodnych z iMSD (ani ich średniej liczby w komórce) w tym samym oknie czasowym, które obserwowano dla makropinosomów. Co więcej, charakterystyczne wartości σ02, Dm oraz α znacząco różnią się od wartości dla lizosomów, które ponownie posłużyły jako punkt odniesienia. Wynik ten potwierdza przypuszczenie, sformułowane powyżej, że granule są badane w „stanie stacjonarnym”, w którym w dowolnym momencie średnie właściwości strukturalne/dynamiczne całej populacji ISG pozostają niezmienione (tj. są stałe, chyba że zmienią się warunki stanu stacjonarnego, na przykład pod wpływem bodźców zewnętrznych).

Rysunek 1: Schemat doświadczenia. (A) Komórki wysiano na szalki do hodowli komórkowych przystosowane do aplikacji mikroskopowych 24 h (48 h w przypadku eksperymentów z transfekcją) przed badaniami konfokalnymi. Następnie komórki poddano odpowiedniej obróbce zgodnie z metodą znakowania (patrz protokół), aby zbarwić interesujące organelle cytoplazmatyczne. (B) Typowa akwizycja konfokalna składa się ze stosu obrazów (time-lapse) fragmentu cytoplazmy żywej komórki, opisującego ewolucję czasową dynamiki znakowanych organelli. (C) Film time-lapse jest analizowany za pomocą specjalnie przygotowanego skryptu Matlab, który najpierw oblicza spatiotemporalną funkcję korelacji obrazu i wykonuje dopasowanie Gaussa w celu wykreślenia krzywych iMSD (D) oraz powiązanych z nimi wyekstrahowanych parametrów dopasowania opisujących parametry dynamiki strukturalnej obrazowanych organelli. Skróty: iMSD = imaging-derived mean square displacement; STICS = spatiotemporal image correlation spectroscopy; α = współczynnik dyfuzji anomalnej; Dm = lokalna dyfuzyjność; σ2(τ) = wariancja. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Odpowiednie parametry eksperymentalne. (A) Przykładowy obraz barwionych lizosomów w utrwalonej próbce. Pasek skali = 2 µm. (B) Pierwsza klatka stosu obrazów barwionych lizosomów w żywej komórce, pozyskana przy zastosowaniu odpowiednich parametrów. Rozdzielczość czasowa: 65 ms/frame. Pasek skali = 2 µm. (C) Pierwsza klatka stosu obrazów barwionych lizosomów w żywej komórce, pozyskana przy niskiej prędkości: widoczna jest artefaktyczna deformacja pozornej wielkości lizosomów wynikająca z ruchu organelli podczas obrazowania. Rozdzielczość czasowa: 10 s/frame. Pasek skali = 2 µm. (D) Przykład obliczania wielkości obrazowanych lizosomów w niebieskim ROI na zdjęciach (A), (B) i (C). Profil intensywności wzdłuż niebieskiej linii dopasowano do funkcji Gaussa, aby wyznaczyć FWHM, czyli szacunkową wielkość plamki. Wartości FWHM podano dla każdego dopasowania. (E) Graficzna reprezentacja wartości wielkości uzyskanych poprzez analizę obrazu opisaną w panelu (D) dla wszystkich obrazowanych lizosomów (czarny kwadrat, wartość średnia i odchylenie standardowe), dla lizosomów znajdujących się w niebieskim ROI (niebieski trójkąt) oraz wyznaczonych za pomocą analizy iMSD (czerwone koło). Rysunek pochodzi z 22. Skróty: ROI = region zainteresowania; FWHM = pełna szerokość połówki maksymalnej; iMSD = średni kwadrat przemieszczenia wyznaczony z obrazowania. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Ewolucja czasowa oznakowanych organelli. (A) Wykresy wartości wielkości wyznaczonych metodą iMSD w funkcji postępu czasowego kolejnych akwizycji, przedstawione jako średnia z wartości zmierzonych w akwizycjach wykonanych w 10-minutowym oknie czasowym. Po lewej stronie postępująca redukcja średniej wielkości makropinosomów (czarne okręgi), a po prawej niezmienna w czasie wielkość wydzielniczych ziarnulek insulinowych (czarne kwadraty) w porównaniu z lizosomami, przedstawionymi jako średnia wartość wielkości (grubsza czerwona linia) ± odchylenie standardowe (przerywane czerwone linie). (B) i (C) Postęp czasowy współczynników Dm i α dla makropinosomów (lewa strona) oraz ziarnulek insulinowych (prawa strona) wyznaczonych analizą iMSD. (D) Ewolucja czasowa liczby oznakowanych makropinosomów i ziarnulek insulinowych zmierzonych w pierwszej klatce każdego zarejestrowanego filmu time-lapse. Skróty: iMSD = imaging-derived mean square displacement; α = współczynnik dyfuzji anomalnej; Dm = dyfuzyjność lokalna, N = liczba. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Organellum | Etykietowanie | Linia komórkowa | Rozmiar (nm) | Dm ( × 10-3 μm2/s) | α | N | Ref. |
| wczesny endosom (EE) | CellLight GFP wczesny endosom | HeLa | 395±74 | 3.0±2.4 | 1.02±0.20 | 40 | 10 |
| późny endosom (LE) | CellLight GFP do późnych endosomów | HeLa | 693±102 | 15.4±10.6 | 0.57±0.16 | 58 | 10 |
| Lizosom (LY) | LysoTracker DND-99 | HeLa | 471±76 | 15.3±9.0 | 0.49±0.13 | 143 | 10, 14 |
| Kawawole (CAV) | Kawoleina-EGFP | HeLa | 405±49 | 3.1±1.8 | 1.00±0.22 | 15 | 10 |
| Pęcherzyki powlekane klatryną (CCV) | Transferryna-Alexa 488 | HeLa | 513±62 | 16.2±9.9 | 0.48±0.17 | 33 | 10 |
| granula insulinowa (IG) | C-peptyd-EGFP | INS-1E | 335±56 | 3.0±1.7 | 0.70±0.14 | 107 | 11 |
| wczesny makropinosom (EMCR) | Dekstran fluoresceiny 70 kDa | HeLa | 979±423 | 8.3±9 | 0.79±0.27 | 36 | 10 |
| śródmiędzypoczątkowy makropinosom (IMCR) | Flurosecyna-dekstran 70 kDa | HeLa | 702±180 | 13.7±19.9 | 0.60±0.38 | 29 | 10 |
| późny makropinosom (LMCR) | Fluroresceina-dekstran 70 kDa | HeLa | 592±127 | 5.8±4.7 | 0.39±0.21 | 21 | 10 |
Tabela 1: Strukturalne i dynamiczne parametry wyznaczone za pomocą iMSD. Tabela przedstawia wartości rozmiaru, Dm oraz współczynnika α zmierzone dla różnych organelli, z określeniem strategii znakowania, użytej linii komórkowej oraz liczby przeanalizowanych serii obrazowania. Wartości podano jako średnia ± odchylenie standardowe. Skróty: iMSD = średni kwadrat przemieszczenia wyznaczony z obrazowania; α = współczynnik dyfuzji anomalnej; Dm = lokalna dyfuzyjność, N = liczba; GFP = zielone białko fluorescencyjne; EGFP = ulepszone zielone białko fluorescencyjne.
Plik uzupełniający 1: Szczegóły dotyczące wyprowadzania i analizy śladu iMSD. Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 1: Prosimy kliknąć tutaj, aby pobrać ten plik.