Wiele białek spełnia swoją funkcję w dynamicznych kompleksach z innymi cząsteczkami, co jest zapośredniczone przez zmiany konformacyjne i przejściowe asocjacje w szerokim zakresie skal czasowych1,2,3. Powiązane z zewnętrznym źródłem energii (np. ATP), te dynamiczne oddziaływania mogą prowadzić do kierunkowości w cyklu funkcjonalnym i w ostateczności utrzymywać w komórce stan stacjonarny poza równowagą, który jest warunkiem koniecznym dla życia.
Aby w pełni zrozumieć działanie tych maszyn molekularnych, opis statyczny oparty na badaniach strukturalnych jest niewystarczający. Ponadto niezbędna jest znajomość leżącego u podstaw modelu kinetycznego oraz wyznaczenie stałych szybkości kinetycznych. Kilka istniejących metod pozwala badaczom na analizę dynamiki oddziaływań binarnych pomiędzy dwiema wybranymi cząsteczkami, np. rezonans plazmonowy powierzchni, metody relaksacyjne z odczytem spektroskopowym (np. techniki skokowe lub przepływu zatrzymanego), oraz rezonans magnetyczny jądra. Jednakże ich zastosowanie jest w większości przypadków ograniczone do prostych układów dwustanowych (np. (jeden stan związany i jeden niezwiązany) ze względu na uśrednianie właściwe dla eksperymentów w skali makroskopowej. W przypadkach, gdy zaangażowanych jest więcej stanów lub produktów pośrednich, otrzymuje się jedynie złożoną mieszaninę stałych szybkości. Metody pojedynczych cząsteczek, takie jak pęsety optyczne lub magnetyczne bądź dwukolorowy smFRET, t.j. Zastosowanie jednego fluoroforu donora i jednego akceptora przy próbce unieruchomionej na powierzchni pozwala na wyznaczenie stałych szybkości wszystkich obserwowanych zmian konformacyjnych. Jednak w przypadku oddziaływań wpływających na więcej niż jedno miejsce wiązania metody te pozostają ograniczone, a informacje o możliwej korelacji między dwoma (lub większą liczbą) oddziaływaniami będą dostępne jedynie poprzez wnioski pośrednie wyciągnięte z serii eksperymentów.
Wielokolorowy smFRET4,5,6,7,8,9 daje możliwość bezpośredniego badania oddziaływań między tymi komponentami w czasie rzeczywistym i w warunkach zbliżonych do fizjologicznych10. Pozwala to na badanie na przykład wiązania ligandu lub innego białka zależnego od konformacji8,9,11. Ogólne podejście przedstawione w niniejszej pracy polega na znakowaniu badanych białek w określonych pozycjach, przytwierdzeniu jednego białka do powierzchni komory pomiarowej i śledzeniu intensywności fluorescencji w czasie przy użyciu TIRFM typu pryzmatycznego (szczegóły patrz 9,12). Bliskość przestrzenną różnych barwników można następnie określić na podstawie transferu energii między nimi. Strategie znakowania mogą się różnić w zależności od białka (przegląd w 13), a istnieją wytyczne pozwalające uniknąć artefaktów w pomiarach smFRET14.
Ponieważ barwnik donora może przekazywać energię do różnych barwników akceptora w wielokolorowym eksperymencie smFRET, względna pozycja wszystkich barwników nie jest dostępna poprzez wzbudzenie tylko jednego barwnika15,16. Jednak w połączeniu z naprzemiennym wzbudzaniem laserowym (ALEX17, opisane w 18) metoda ta dostarcza pełnych informacji czasowo-przestrzennych z rozdzielczością poniżej sekundy i nanometra.
W zasadzie wysoką rozdzielczość informacji strukturalnych można uzyskać, wykorzystując odległości między barwnikami obliczone na podstawie kombinacji wszystkich intensywności fluorescencji w wielobarwnym eksperymencie smFRET z zastosowaniem ALEX. W niniejszej pracy skupiamy się jednak na identyfikacji i rozdzielaniu stanów, a także na wywodzeniu modeli kinetycznych, w których wielobarwny smFRET jest niezbędny. W przypadku, gdy pożądanym jest „wyłącznie” wyznaczenie struktury metodą triangulacji, można przeprowadzić zestaw prostszych dwubarwnych eksperymentów smFRET o wysokim stosunku sygnału do szumu12,19.
Wykorzystujemy częściową fluorescencję (
) jako wskaźnik transferu energii między dwoma fluoroforami7. Wartość PF jest obliczana na podstawie intensywności fluorescencji analogicznie do wydajności FRET w eksperymencie dwukolorowym:

Gdzie
to natężenie w kanale emisji em po wzbudzeniu kolorem ex, a c to akceptor o najdłuższej długości fali. Kanały detekcji reprezentują tę samą pozycję w komorze z próbką, ale rejestrują różne zakresy widmowe światła fluorescencyjnego. W niniejszym protokole stosuje się te same identyfikatory dla wzbudzenia i emisji (t.j. „niebieski”, „zielony” i „czerwony”).
Z powodu niedociągnięć eksperymentalnych zmierzone intensywności fluorescencji zależą nie tylko od przeniesienia energii, ale także od właściwości fluoroforów i konfiguracji układu. Aby uzyskać rzeczywistą wydajność przeniesienia energii między dwoma fluoroforami, zmierzone intensywności muszą zostać skorygowane. Poniższa procedura opiera się na referencji9. Współczynniki korekcyjne dla pozornego wycieku (lk, tj. detekcji fotonów z jednego fluoroforu w kanale przeznaczonym dla innego barwnika) oraz pozornej wartości gamma (ag, tj. wydajności kwantowej fluorescencji barwnika i wydajności detekcji kanału) są wyznaczane z pojedynczych śladów cząsteczkowych wykazujących zdarzenie wybielenia akceptora.
Przeciek barwnika donora do każdego z możliwych kanałów akceptora oblicza się na podstawie wszystkich punktów danych zarejestrowanych śladów fluorescencji, w których barwnik akceptora uległ wybieleniu, a donor nadal wykazuje fluorescencję (
):

Medianę histogramu wycieku przyjmuje się jako pozorny czynnik wycieku. Po skorygowaniu o wyciek, z tego samego zestawu śladów wyznacza się pozorny czynnik gamma. Oblicza się go poprzez podzielenie zmiany fluorescencji w kanale akceptora przez zmianę fluorescencji w kanale donora po wybieleniu barwnika akceptora:

Gdzie c ponownie oznacza kanał detekcji akceptora o najdłuższej długości fali. Mediana wynikowego rozkładu jest wykorzystywana jako pozorny współczynnik korekcji.
Skorygowane intensywności w każdym kanale uzyskuje się poprzez:

Następnie PF oblicza się zgodnie z wzorem:

Różne populacje można rozdzielić w wielowymiarowej przestrzeni rozpiętej przez PF. Położenie i szerokość każdego stanu są określane poprzez dopasowanie danych do wielowymiarowych funkcji Gaussa. Następna optymalizacja jednego globalnego modelu HMM w oparciu o wszystkie ślady PF zapewnia ilościowy opis zaobserwowanej kinetyki. Wykrywalne są nawet niewielkie zmiany szybkości reakcji.
Ukryte modele Markowa (HMM) umożliwiają wnioskowanie modelu stanów z zestawu zaszumionych przebiegów czasowych. Przyjmuje się, że w dowolnym momencie system znajduje się w jednym z zestawu dyskretnych, ukrytych stanów, a rzeczywista obserwacja (t.j. emisja) jest funkcją probabilistyczną tego stanu ukrytego20. W przypadku danych TIRFM smFRET prawdopodobieństwa emisji bi dla stanu i mogą być modelowane za pomocą ciągłych funkcji gęstości prawdopodobieństwa Gaussa. W regularnie rozmieszczonych dyskretnych punktach czasowych mogą następować przejścia z jednego stanu do drugiego zgodnie z prawdopodobieństwem przejścia, które jest niezmienne w czasie i zależy wyłącznie od stanu bieżącego. Macierz przejść A zawiera te prawdopodobieństwa przejść aij pomiędzy wszystkimi stanami ukrytymi. Rozkład stanu początkowego
określa prawdopodobieństwa specyficzne dla stanu
dla pierwszego punktu czasowego przebiegu. Wykorzystując podejście największej wiarygodności, parametry te można zoptymalizować w celu najlepszego opisania danych za pomocą algorytmów Forward-Backward i Baum-Welcha20,21. Pozwala to uzyskać estymatory największej wiarygodności (MLE). Ostatecznie, sekwencję stanów, która najprawdopodobniej wygenerowała trajektorię obserwacji, można wywnioskować za pomocą algorytmu Viterbi. W przeciwieństwie do innych analiz HMM danych smFRET24,25,26 nie używamy HMM jako zwykłego „wygładzania” danych, lecz wyodrębniamy kinetyczny model stanów z zestawu danych bez konieczności dopasowywania histogramów czasu przebywania27. Analiza HMM jest przeprowadzana za pomocą skryptów wewnętrznych w programie Igor Pro. Implementacja kodu opiera się na referencji21. Na naszej stronie internetowej udostępniamy zestaw oprogramowania i przykładowe dane, aby umożliwić realizację sekcji 5 i 6 niniejszego protokołu (https://www.singlemolecule.uni-freiburg.de/software/3d-fret). Pełne oprogramowanie jest dostępne na żądanie.
Punkty czasowe w danych, dla których PF <-1 lub PF >2 w dowolnym kanale detekcji, otrzymują minimalne prawdopodobieństwo emisji dla wszystkich stanów (10-200). Zapobiega to sztucznym przejściom w tych punktach danych.
Parametry prawdopodobieństw emisji są uzyskiwane z dopasowania histogramu 3D PF do funkcji Gaussa, zgodnie z opisem w kroku 5.7. Parametry te pozostają stałe podczas optymalizacji HMM.
W przedstawionym podejściu początkowy wektor rozkładu stanów oraz macierz przejść są stosowane globalnie do opisu całego zespołu śladów. Są one aktualizowane na podstawie wszystkich N cząsteczek z zestawu danych zgodnie z referencją27.
Parametry początkowe dla rozkładu stanów początkowych są określane na podstawie projekcji 2D histogramu PF (krok 5.3), a prawdopodobieństwa przejścia ustawia się na 0,05, z wyjątkiem prawdopodobieństw pozostania w tym samym stanie, które dobiera się tak, aby prawdopodobieństwo opuszczenia danego stanu było znormalizowane do jedności.
Do wyznaczenia przedziałów ufności (CIs) dla wszystkich stałych szybkości przejść21,22, które służą jako istotne oszacowania ich niepewności, zastosowano metodę profilowania wiarygodności. Aby obliczyć granice przedziału ufności dla konkretnej stałej, prawdopodobieństwo przejścia, które jest przedmiotem analizy, ustala się na wartość inną niż MLE. W ten sposób otrzymuje się model testowy λ'. Test ilorazu wiarygodności (LR) wiarygodności
dla zestawu danych 0 przeprowadza się zgodnie z:

Granicę ufności 95% dla parametru osiąga się, gdy LR przekroczy 3,841, czyli kwantyla 95% rozkładu x2 z jednym stopniem swobody22,23.
Skuteczność tej metody zademonstrowano na przykładzie Hsp90. To powszechnie występujące białko znajduje się u bakterii i eukariotów i stanowi część komórkowej odpowiedzi na stres28. Jest ono obiecującym celem terapeutycznym w leczeniu nowotworów29. Hsp90 jest homodimerem z jedną kieszenią wiążącą nukleotyd w domenie N-końcowej każdej podjednostki30. Może ono przechodzić przejścia pomiędzy co najmniej dwiema globalnie odmiennymi konformacjami: zamkniętą oraz otwartą w obrębie N-końca, przyjmującą kształt litery V19,31,32. Dimericzna natura białka bezpośrednio nasuwa pytanie o interakcje pomiędzy dwoma miejscami wiążącymi nukleotydy w Hsp90.
Poniżej przedstawiamy szczegółowy protokół pozyskiwania i analizy danych w eksperymencie trójkolorowego smFRET dotyczącym nukleotydów oraz białka Hsp90 z drożdży. Analizie poddano zależne od konformacji wiązanie fluorescencyjnie znakowanego AMP-PNP (AMP-PNP*, niehydrolizowalnego analogu ATP). Zastosowanie opisanego postępowania umożliwia badanie wiązania nukleotydów przy jednoczesnej analizie zmian konformacyjnych Hsp90, co pozwala ujawnić kooperatywność pomiędzy dwiema kieszeniami wiążącymi nukleotydy w białku Hsp90.