Aby zilustrować możliwości aparatury openFLIM-HCA, przedstawiamy trzy przykładowe zastosowania FRET. Pierwsze z nich dotyczy konstruktów FRET zaadaptowanych z modelowych konstruktów FRET opracowanych w laboratorium Stephena Vogela38. Składają się one z serii genetycznie ekspresowalnych konstruktów fluorescencyjnych, w których białko fluorescencyjne donora (mTurquoise) jest połączone z fluoroforem akceptorem (Venus) krótkimi, kontrolowanymi łącznikami o długości 5, 17 i 32 aminokwasów (mTq5V, mTq17V, mTq32V). Różne długości łączników między fluoroforami donora i akceptora skutkują różną wydajnością FRET, a tym samym różnym czasem życia donora. W celu zapewnienia kontroli negatywnej, w konstrukcie z krótkim, 5-aminokwasowym łącznikiem, Venus został zastąpiony białkiem Amber — niefluorescencyjną mutacją YFP (mTq5A), która nie może pełnić roli akceptora FRET38. Zastąpiliśmy białko Cerulean w oryginalnych wektorach pCXV i pC5A poprzez trawienie restrykcyjne wektorów enzymami NheI i BglII oraz ligację fragmentów mTurquoise, trawionych tymi samymi enzymami z wektora pmTurquoise-N1. Obszar łącznika nie został zmodyfikowany podczas tej substytucji.
Rysunek 6 przedstawia przykładowy pomiar FLIM FRET z wykorzystaniem tego modelu systemu wyrażonego w komórkach Cos, w którym na płytce wielodobkowej rozmieszczono po 3 kolumny komórek przetransfekowanych odpowiednio: kontrolą negatywną (kolumny 1-3), konstruktem FRET z linkerem o długości 5 aminokwasów (kolumny 4-6), konstruktem FRET z linkerem o długości 17 aminokwasów (kolumny 7-9) oraz konstruktem FRET z linkerem o długości 32 aminokwasów (kolumny 10-12). Wiersz H zawiera studzienki do szacowania TVB. Rysunek 6(a) pokazuje przykłady automatycznie pozyskanych obrazów czasu życia fluorescencji z typowych pól widzenia w każdej studzience. Widać, że kontrola negatywna (kolumny 1-3) wykazuje najdłuższe czasy życia, a czas życia donora jest najniższy dla konstruktu FRET z najkrótszym linkerem (kolumny 4-6). Rysunek 6(b) pokazuje, jak średni czas życia uśredniony dla wszystkich pól widzenia w każdej studzience zmienia się na całej płytce wielodobkowej. Rysunki 6(c) oraz (d) przedstawiają przykładowe mapy czasu życia z nałożoną intensywnością fluorescencji donora dla pola widzenia odpowiednio mTq5A i mTq17V. Rysunek 6(e) pokazuje profil zaniku intensywności fluorescencji donora uśredniony dla wszystkich komórek obrazowanych w studzience A1, wraz z dopasowaniem do monoeksponencjalnego modelu zaniku oraz IRF (która jest zdominowana przez szerokość bramki GOI). Rysunek 6(f) przedstawia średni czas życia fluorescencji dla każdej kolumny płytki wielodobkowej, uzyskany z pikselowego dopasowania monoeksponencjalnego. Tabelę średniego czasu życia i odchylenia standardowego (STD) można znaleźć w Tabeli 1 poniżej. Pozyskanie danych dla obrazów pokazanych na Rysunku 6 trwało około 160 min. Analiza zajęła 92 s na komputerze 2,6 GHz z 10 rdzeniami i 64 GB pamięci RAM.
Drugi przykładowy test demonstruje zastosowanie FLIM FRET do odczytu oligomeryzacji białka HIV-1 Gag podczas składania wirionów HIV-1 jako sposobu testowania inhibitorów tego procesu25,42. Ekspresja HIV-1 Gag w odpowiednich komórkach gospodarza stanowi system modelowy dla tego etapu cyklu życiowego HIV-1, ponieważ prowadzi do powstawania cząsteczek przypominających wirusy (VLP), które są podobne do niedojrzałych wirionów HIV-1. W tym teście wyeksponowano białko HIV-1 Gag połączone z CFP w komórkach HeLa i porównano działanie różnych inhibitorów poprzez analizę ich wpływu na sygnał FRET wynikający z oligomeryzacji białka Gag. Szczegółowe informacje na temat zastosowanych inhibitorów znajdują się w materiale referencyjnym 42.
Szczegóły przygotowania próbek i pomiarów eksperymentalnych zostały szczegółowo opisane w referencji 25, gdzie wykazaliśmy, że FLIM można wykorzystać do odczytu zarówno heteroFRET pomiędzy agregującymi białkami Gag znakowanymi stochastycznie CFP i YFP, jak i homoFRET w komórkach wykazujących ekspresję wyłącznie białka Gag znakowanego CFP. Chociaż homoFRET zazwyczaj nie powoduje zmiany czasu życia donora, ma to miejsce w szczególnym przypadku CFP, gdzie fluorofor występuje w postaci dwóch izomerów o różnych średnich czasach życia fluorescencji40,41. Odczyt homo-FRET dla CFP charakteryzuje się zaletami w postaci uproszczonego przygotowania próbek w porównaniu do hetero-FRET CFP/YFP oraz większą wydajnością spektralną, co może być istotne w przypadku multipleksowych odczytów FRET.
W naszych poprzednich badaniach zastosowano technikę FLIM FRET w celu oceny oligomeryzacji białka Gag w obecności inhibitora N-mirystylotransferazy (NMT)42. Inhibitor ten zaburza działanie endogennych enzymów odpowiedzialnych za przyłączanie kwasu mirystylowego do białka Gag, co umożliwia wiązanie białek Gag z błoną plazmatyczną i jest warunkiem niezbędnym43 do formowania wirionów HIV lub VLPs. Wykazaliśmy, że odczyty zarówno heteroFRET, jak i homoFRET mogą osiągnąć wartość Z' >0,6 w badaniach zależności dawka-odpowiedź z inhibitorem NMT. Z' jest parametrem stosowanym w przemyśle farmaceutycznym do określania jakości testu; reprezentuje on różnice w wartościach średnich kontroli pozytywnej i negatywnej oraz ich względne rozrzuty, tworząc pojedynczą miarę jakości testu — przy czym wartość Z' > 0,5 jest uznawana za „doskonałą” w przypadku testów farmaceutycznych44. Zauważamy, że tak wysoka wydajność została osiągnięta w próbkach, w których zmiana średniego czasu życia donora wynosiła rzędu 300 ps — co dowodzi statystycznej mocy analizy danych FLIM dla dużej liczby komórek, co pozwala na uzyskanie użytecznych odczytów przy znacznie mniejszych zmianach czasu życia fluorescencji, niż jest to możliwe przy typowej liczbie komórek obrazowanych w manualnych eksperymentach mikroskopowych.
W niniejszym materiale przedstawiamy zestaw danych FLIM FRET z płytki wielodołkowej, w którym porównano cztery związki hamujące oligomeryzację białka Gag HIV (oznaczone jako ICL13, ICL14, ICL15 & ICL16). W tym eksperymencie wykorzystano odczyt homoFRET w komórkach HeLa przejściowo transfekowanych plazmidem Gag-CFP. Zgodnie ze standardowym protokołem opisanym w referencji 32, zastosowano łącznie dziewięć różnych dawek każdego inhibitora, przewidując dwie powtórne studzienki na warunek. Jako kontrolę pozytywną w kolumnie 1 przedstawiono komórki wykazujące ekspresję myr-Gag, zmutowanego białka Gag pozbawionego grupy kwasu mirystynowego, które nie może tworzyć VLP, co symuluje całkowite zahamowanie. Kontrolę negatywną stanowiły komórki wykazujące ekspresję Gag-CFP, którym podano wyłącznie nośnik wykorzystany do wprowadzenia inhibitorów w pozostałych kolumnach (kolumna 11). W każdej studzience pozyskano łącznie osiem pól widzenia (FOV).
Dane dopasowano do modelu pojedynczego rozpadu wykładniczego w skali poszczególnych pikseli, a mapa płyty czasu życia fluorescencji przedstawiająca średni czas życia na studzienkę (uśredniony dla wszystkich pikseli powyżej progu intensywności we wszystkich ośmiu polach widzenia) jest przedstawiona poniżej na Rysunku 7(a). Rysunek 7(b) przedstawia odpowiadający temu wykres czasu życia fluorescencji w funkcji stężenia inhibitora. Trzy z inhibitorów wykazują efekt hamujący (ICL13, ICL15 & ICL16) podczas inkubacji z komórkami wykazującymi ekspresję Gag-CFP. Jeden inhibitor (ICL14) nie wykazuje żadnego efektu przy żadnej badanej dawce. Wartość Z' obliczona dla tej płyty wyniosła 0,51. Wartości uzyskane dla kontroli pozytywnej i negatywnej są przedstawione na Rysunku 7(b) jako punkty w kolorze czarnym przy 100 µM i 0,1 nM.
Krzywe zależności dawka-odpowiedź dla związków ICL13, ICL15 i ICL16 przedstawione na Rysunku 7(b) dopasowano do czteroparametrowego logistycznego modelu regresji nieliniowej z współczynnikiem Hilla ustalonym na 145, przy czym maksymalny i minimalny czas życia ustalono odpowiednio na wartościach uzyskanych z kontroli pozytywnej (kolumna 1) i negatywnej (kolumna 11). Uzyskane wartości IC50 dla tych trzech krzywych wyniosły odpowiednio: 38 nM, 24 nM i 116 nM dla ICL13, ICL15 i ICL16. W celu porównania z wartościami IC50 otrzymanymi za pomocą wewnątrzkomórkowych pomiarów FLIM FRET, określiliśmy IC50 dla hamowania aktywności enzymatycznej rekombinowanej ludzkiej NMT1 w naszym wcześniej opisanym biochemicznym teście enzymatycznym42. Trzy związki, które okazały się aktywne w teście FLIM FRET, są również najaktywniejsze w tym teście (wartości IC50 wynoszące odpowiednio 17 nM, 51 nM i 39 nM dla ICL13, ICL15 i ICL16), podczas gdy związek ICL14, który nie wykazał istotnej odpowiedzi w teście FLIM FRET, był ca. 100-krotnie mniej aktywny wobec NMT1 (IC50 wynoszące 1700 nM). W przypadku związków aktywnych wartości IC50 uzyskane za pomocą tych niezależnych metod analitycznych są niezwykle zbliżone, a niewielkie różnice można prawdopodobnie przypisać odmiennym różnicom w pobieraniu lub metabolizmie związku w komórkach.
Niniejsze niewielkie badanie podkreśla zdolność metody FLIM HCA do rozróżniania potencji różnych związków działających na ten sam wybrany cel. W naszej ocenie jest to pierwsza demonstracja przesiewu z wykorzystaniem zautomatyzowanego FLIM w kontekście wysokiej zawartości danych (high content) w celu wygenerowania wielu krzywych zależności dawka-odpowiedź, co ilustruje potencjał zastosowania FLIM do poszukiwania i/lub charakteryzowania użytecznych związków terapeutycznych.
Trzeci przykładowy eksperyment FLIM FRET dotyczy genetycznie ekspresyjnego biosensora FRET dla białka EPAC (Exchange protein activated by cyclic adenosine monophosphate), które jest czynnikiem wymiany guaniny dla Rap-1, wiążącym specyficznie cAMP. Ten biosensor FRET EPAC (EPAC-SH188) wykorzystuje białko fluorescencyjne mTurquoise2 (mTq2FP) jako fluorofor donora i zawiera dwa białka Venus-FP w celu stworzenia złożonego akceptora46. cAMP jest powszechnym drugim przekaźnikiem i bierze udział w wielu procesach wewnątrzkomórkowych u wielu różnych organizmów. Jest on syntetyzowany przez cyklazę adenylanową z ATP w błonie komórkowej. Tutaj demonstrujemy możliwość odczytu i kwantyfikacji jego aktywności w żywych komórkach HEK293T po dodaniu forskoliny, aktywatora cyklazy adenylanowej, który zwiększa wewnątrzkomórkową produkcję cAMP, a tym samym podnosi jego stężenie w cytoplazmie. Sensor EPAC wykazuje zjawisko FRET w stanie nieaktywnym (niezwiązanym), a czas życia donora wzrasta wraz ze stężeniem cAMP, ponieważ cAMP prowadzi do otwarcia biosensora. Monoeksponencjalny profil zaniku mTq2FP sprawia, że biosensor ten jest lepiej dostosowany do ilościowej analizy czasu życia niż biosensory oparte na CFP 45. Rysunek 8 przedstawia przykład przebiegu czasowego zarejestrowanego za pomocą zautomatyzowanego mikroskopu FLIM z płytką wielodołkową, gdzie wykreślono zmianę średniego czasu życia oraz zmianę frakcji populacji donora uczestniczącego w FRET w czasie po dodaniu 100 µM forskoliny. Dla uproszczenia przyjmujemy, że całkowity sygnał można opisać jako mieszaninę monoeksponencjalnych donorów uczestniczących i nieuczestniczących w FRET, a frakcję populacji aktywowanego biosensora FRET EPAC uzyskano poprzez dopasowanie profilu zaniku dwueksponencjalnego w całej serii czasowej.
Do akwizycji danych cAMP (EPAC) wybrano pięć opóźnień bramki przy szerokości bramki 4,000 ps, z czasem opóźnienia ustawionym na 1,300 ps, 4,300 ps (intensywność szczytowa), 4,919 ps, 6,656 ps, 8,035 ps, 1,0391 ps oraz 16,000 ps. Czas integracji kamery dla każdego opóźnienia wynosił 250 ms. W jednej studience zarejestrowano przebieg czasowy FLIM, wykonując obrazy co 10 s przez 10 min. Punkty danych uzyskane w czasie 120 s i 130 s zostały usunięte, ponieważ dodanie forskoliny spowodowało niewielkie przesunięcie pozycji ogniskowej próbki, a tym samym intensywności fluorescencji rejestrowanej podczas akwizycji FLIM w tych punktach czasowych.
Do analizy danych obrazy wczytano do programu FLIMfit i przeprowadzono segmentację na poziomie pojedynczych komórek. Dane dotyczące czasu życia fluorescencji dopasowano do modelu podwójnego rozpadu wykładniczego z wykorzystaniem IRF oraz stałej wartości t0 uzyskanej z barwnika referencyjnego (75 µM Coumarin 6). Średni czas życia fluorescencji donora uśredniony dla wszystkich komórek przedstawiono na Rysunku 8(a). Rysunek 8(b) przedstawia zmianę frakcji populacji wykazującej FRET w czasie, obliczoną poprzez dopasowanie tych samych danych FLIM do tego samego modelu podwójnego rozpadu wykładniczego

gdzie β to frakcja donora uczestniczącego w FRET. W przypadku punktów czasowych poprzedzających dodanie forskoliny przyjęto istnienie mieszaniny elementów uczestniczących i nieuczestniczących w FRET. W celu wyznaczenia tych czasów życia, do trzech pierwszych punktów czasowych zastosowano dopasowanie dwueksponencjalne, co pozwoliło uzyskać czasy życia τ1 = 3,964 ± 21 ps dla donora nieuczestniczącego w FRET oraz τ2 = 3,022 ± 27 ps dla donora uczestniczącego w FRET. Wartości te zostały ustalone dla późniejszego dopasowania całego zestawu danych w celu uzyskania wartości β dla każdego punktu czasowego.
Podsumowując, zaprezentowaliśmy przykładowe wyniki HCA metodą FLIM dla trzech próbek eksperymentalnych. Pierwszą była płytka 96-dołkowa z komórkami Cos wykazującymi ekspresję konstruktów FRET, składających się z donora mTqFP połączonego z akceptorem Venus FP lub niefluorescencyjnym konstruktem Amber za pomocą łączników peptydowych o różnej długości. Zmiana wydajności FRET, a tym samym czasu życia donora w zależności od długości łącznika, jest wyraźnie widoczna, a odchylenie standardowe średniego czasu życia dla każdego konstruktu wynosiło mniej niż 36 ps. Drugim przykładowym badaniem był „przesiew” czterech potencjalnych inhibitorów mirystoilacji białka Gag HIV, dla których obliczono % inhibicji na podstawie zmienności średniego czasu życia fluorescencji donora w zależności od stężenia inhibitora, zmierzonego w komórkach Hela wykazujących ekspresję białka Gag znakowanego CFP. Odczyt ten opiera się na homoFRET, który zazwyczaj nie powoduje zmiany czasu życia donora, jednak w przypadku CFP następuje ona, ponieważ białko to występuje w wielu izomerach o różnych czasach życia fluorescencji. Może to być użyteczne pod kątem wydajności spektralnej, np. przy multipleksowaniu różnych odczytów FRET25. Trzecim przykładowym badaniem był monitoring w czasie rzeczywistym żywych komórek HEK293T wykazujących ekspresję biosensora FRET dla cAMP, EPAC. Wykazano oczekiwaną odpowiedź po zastosowaniu forskoliny oraz wykorzystaniu programu FLIMfit z modelem podwójnego wykładniczego zaniku, przy czym liczba wykrytych fotonów była współmierna do obrazowania żywych komórek — co wykazaliśmy wcześniej w przypadku biosensora FRET LIBRA dla IP347.

Rycina 1. Schemat szerokopolowego czasowo-bramkowanego FLIM z wykorzystaniem bramkowanego wzmacniacza obrazu optycznego (GOI). Po lewej stronie znajduje się schemat układu doświadczalnego. W prawym górnym rogu schemat impulsu wzbudzenia, położenie obrazów czasowo-bramkowanych oraz dopasowanie monoeksponencjalne do danych. W prawym dolnym rogu rysunek przedstawiający zanik fluorescencji dwóch obiektów ukazanych w układzie doświadczalnym. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2. Schemat szerokopolowego openFLIM-HCA z wykorzystaniem bramkowanego wzmacniacza obrazu optycznego (GOI). Szczegółowy opis komponentów systemu znajduje się w tekście głównym. Aby wyświetlić ten rysunek w większym rozmiarze, kliknij tutaj.

Rysunek 3. Schemat optycznie sekcjonowanego systemu openFLIM-HCA z wykorzystaniem bramkowanego wzmacniacza obrazu optycznego (GOI) z jednostką skanera z tarczą Nipkowa. Szczegółowy opis komponentów systemu znajduje się w tekście głównym. Aby wyświetlić powiększoną wersję tego rysunku, kliknij tutaj.

Rysunek 4. Schemat przepływu danych obrazowych z programu akwizycyjnego Manager do serwera OMERO lub dysku komputera oraz do programu FLIMfit w celu analizy.Przepływ danych rozpoczyna się w lewym górnym rogu, gdzie surowe dane obrazowe są pozyskiwane w oprogramowaniu akwizycyjnym μ-Manager. Elementy wewnątrz przerywanej ramki reprezentują etapy analizy danych FLIM, natomiast elementy na zewnątrz odpowiadają akwizycji i przechowywaniu danych. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 5. Zrzut ekranu interfejsu użytkownika FLIMfit. W lewym górnym rogu lista obecnie załadowanych pól widzenia oraz obraz intensywności fluorescencji dla aktualnie wybranego pola widzenia. W lewym dolnym rogu wybór modelu dopasowania i warunków dopasowania. W prawym górnym rogu zanik fluorescencji dla aktualnego obszaru zainteresowania (niebieskie okręgi), dopasowanie do danych (czerwona linia), IRF (czerwona przerywana linia) oraz rezydua dopasowania poniżej (niebieska linia). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 6. FLIM na płytce wielodołkowej dla modelowych konstruktów FRET mTq5V, mTq17V, mTq32V ekspresowanych w komórkach Cos. (a) przykładowe obrazy czasu życia fluorescencji dla pól widzenia (FOV) w każdej dołku dla różnych konstruktów FRET ekspresowanych w komórkach Cos, zgodnie z opisem w tekście. (b) mapa ciepła średnich czasów życia fluorescencji Venus uśrednionych dla wszystkich komórek w każdym dołku. (c) obraz intensywności mTurquoise Amber nałożony na mapę czasu życia (skala 20 µm). (d) obraz intensywności mTurquoise Venus (17 aminokwasów) nałożony na mapę czasu życia (skala 20 µm). (e) zmierzony profil zaniku intensywności (niebieskie kółka) fluorescencji konstruktu mTq5A (kontrola negatywna) uśredniony dla wszystkich komórek obrazowanych w dołku A1, wraz z dopasowaniem danych do zaniku monoeksponencjalnego (niebieska linia ciągła) oraz IRF (pomarańczowa linia przerywana). (f) wykres średniego czasu życia uśrednionego dla każdej kolumny na płytce wielodołkowej, z słupkami błędów reprezentującymi odchylenie standardowe czasu życia fluorescencji między polami widzenia. Dla (a-d) wartości czasu życia przedstawiono za pomocą fałszywej skali kolorów z wskazanymi granicami w pikosekundach. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 7. Przykładowy screening FLIM inhibitorów Gag na płytce wielodołkowej z użyciem komórek HeLa eksponujących Gag-CFP. (a) Mapa ciepła średniego czasu życia fluorescencji donora na dołek dla płytki wielodołkowej, w której kolumna 1 przedstawia komórki transfekowane myr-Gag-CFP jako kontrolę pozytywną inhibicji, kolumna 11 przedstawia komórki transfekowane Gag-CFP poddane działaniu samego nośnika, a przerywane kolorowe kwadraty wskazują dołki, w których zastosowano jeden z czterech różnych inhibitorów w stężeniach od wysokiego w kolumnie 2 do niskiego w kolumnie 10. Fałszywa skala kolorów reprezentuje średni czas życia fluorescencji w zakresie od 2,200 ps do 3,100 ps. (b) Wykresy czasu życia fluorescencji dla każdego inhibitora z paskami błędów reprezentującymi odchylenie standardowe między powtórzeniami w dołkach. Czarne punkty w położeniach 2 i -4 na osi poziomej to odpowiednio punkty kontroli pozytywnej i negatywnej, uzyskane z kolumn 1 i 11. Linie ciągłe wskazują dopasowanie do danych krzywej odpowiedzi na dawkę dla różnych inhibitorów. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rysunek 8. Przykładowe zastosowanie FLIM do odczytu biosensora FRET EPAC w pomiarze time-lapse z wykorzystaniem komórek HEK293T. Zmiana (a) średniego czasu życia fluorescencji oraz (b) frakcji donora uczestniczącego w FRET w funkcji czasu po dodaniu forskoliny, zaznaczonej zielonym paskiem. Słupki błędów przedstawiają błąd standardowy na komórkę. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.
| na dołek | STD | BŁĄD | | | na dołek | STD | Błąd |
| mT5A | 4008 | 32 | 12 | | mT17V | 3004 | 26 | 10 |
| mT5V | 2717 | 35 | 13 | | mT32V | 3133 | 30 | 11 |
Tabela 1. Średni czas życia i odchylenie standardowe (STD) konstruktów FRET.