$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Wiele białek spełnia swoją funkcję w dynamicznych kompleksach z innymi cząsteczkami, w których pośredniczą zmiany konformacyjne i przejściowe asocjacje w szerokim zakresie skal czasowych1,2,3. W połączeniu z zewnętrznym źródłem energii (np. ATP) te dynamiczne interakcje mogą prowadzić do kierunkowości w cyklu funkcjonalnym i ostatecznie utrzymać nierównowagowy stan ustalony w komórce, co jest warunkiem wstępnym życia.
Aby w pełni zrozumieć te molekularne maszyny, statyczny opis oparty na badaniach strukturalnych nie jest wystarczający. Ponadto niezbędna jest wiedza na temat leżącego u podstaw modelu kinetycznego i określenie stałych szybkości kinetycznej. Kilka istniejących metod pozwala naukowcom badać dynamikę oddziaływań binarnych między dwiema interesującymi cząsteczkami, np. powierzchniowy rezonans plazmonowy, metody relaksacji z odczytem spektroskopowym (np. techniki skoku lub zatrzymanego przepływu) oraz magnetyczny rezonans jądrowy. Jednak ich zastosowanie jest w większości przypadków ograniczone do prostych systemów dwustanowych (np. jeden stan związany i jeden stan niezwiązany) ze względu na uśrednianie nieodłącznie związane z eksperymentami masowymi. W przypadkach, gdy zaangażowanych jest więcej stanów lub produktów pośrednich, dają one tylko złożoną mieszaninę stałych szybkości. Metody jednocząsteczkowe, takie jak pęseta optyczna lub magnetyczna lub dwukolorowy smFRET, tj. jeden fluorofor donorowy i jeden akceptorowy, z próbką unieruchomioną powierzchniowo mogą przywrócić stałe szybkości dla wszystkich obserwowanych zmian konformacyjnych. Jednakże, jeśli chodzi o interakcje wpływające na więcej niż jedno miejsce wiązania, metody te pozostają ograniczone, a informacje na temat możliwej korelacji dwóch (lub więcej) interakcji będą dostępne tylko poprzez pośrednie wnioski z zestawu eksperymentów.
Wielokolorowy smFRET4,5,6,7,8,9 oferuje możliwość badania interakcji między tymi komponentami bezpośrednio, w czasie rzeczywistym i w warunkach zbliżonych do fizjologicznych10. Pozwala to na zbadanie na przykład zależnego od konformacji wiązania liganda lub innego białka8,9,11. Ogólne podejście przedstawione tutaj polega na oznaczaniu interesującego białka (białek) w określonych pozycjach, dołączeniu jednego białka do powierzchni komory pomiarowej i śledzeniu intensywności fluorescencji w czasie na pryzmatycznym TIRFM (szczegóły patrz 9,12). Przestrzenną bliskość 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 (omówione w 13) i istnieją wytyczne dotyczące unikania artefaktów w pomiarach smFRET14.
Ponieważ barwnik donorowy może przekazywać energię do różnych barwników akceptorowych w wielokolorowym eksperymencie smFRET, względna pozycja wszystkich barwników nie jest dostępna na podstawie wzbudzenia tylko jednego barwnika15,16. Ale w połączeniu z naprzemiennym wzbudzeniem laserowym (ALEX17 i omówionym w 18) metoda ta dostarcza wszystkich informacji czasoprzestrzennych w rozdzielczości subsekundowej i subnanometrowej.
Zasadniczo, informacje strukturalne o wysokiej rozdzielczości można uzyskać, korzystając z odległości między barwnikami obliczonymi z kombinacji wszystkich intensywności fluorescencji w wielokolorowym eksperymencie smFRET z ALEX. Tutaj jednak skupiamy się na identyfikacji i separacji stanów oraz ekstrakcji modeli kinetycznych, gdzie wielokolorowy smFRET jest niezbędny. Gdy pożądane jest "tylko" określenie struktury za pomocą triangulacji, można przeprowadzić zestaw prostszych dwukolorowych eksperymentów smFRET o wysokim stosunku sygnału do szumu12,19.
Używamy częściowej fluorescencji (
) jako wskaźnika zastępczego dla transferu energii między dwoma fluoroforami7. PF oblicza się na podstawie intensywności fluorescencji analogicznej do wydajności FRET w eksperymencie dwukolorowym:

Gdzie,
to natężenie w kanale emisyjnym em po wzbudzeniu kolorem ex, a c to akceptor o najdłuższej długości fali. Kanały detekcyjne reprezentują to samo położenie w komorze próbki, ale rejestrują różne zakresy widmowe światła fluorescencyjnego. W tym protokole używany jest ten sam identyfikator wzbudzenia i emisji (tj. "niebieski", "zielony" i "czerwony").
Z powodu niedociągnięć eksperymentalnych, zmierzone intensywności fluorescencji zależą nie tylko od transferu energii, ale także od właściwości fluoroforu i ustawienia. Aby uzyskać rzeczywistą wydajność transferu energii między dwoma fluoroforami, zmierzone intensywności muszą zostać skorygowane. Poniższa procedura jest oparta na reference9. Współczynniki korekcyjne dla pozornego wycieku (lk, tj. detekcja fotonów z fluoroforu w kanale przeznaczonym dla innego barwnika) i pozornego gamma (ag, tj. wydajność kwantowa fluorescencji barwnika i skuteczność wykrywania kanału) uzyskuje się ze śladów pojedynczej cząsteczki, które wykazują zdarzenie bielenia akceptora.
Wyciek barwnika dawcy do każdego możliwego kanału akceptora jest obliczany na podstawie wszystkich punktów danych w zarejestrowanych śladach fluorescencji, w których barwnik akceptorowy wybielił, ale dawca nadal jest fluorescencyjny (
):

Mediana histogramu wycieku jest używana jako współczynnik pozornego wycieku. Po korekcie o nieszczelność pozorny współczynnik gamma jest wyznaczany z tego samego zestawu ścieżek. Oblicza się go, dzieląc zmianę fluorescencji w kanale akceptorowym przez zmianę fluorescencji w kanale donorowym po bieleniu barwnika akceptorowego:

Gdzie c ponownie jest kanałem detekcji dla akceptora o najdłuższej długości fali. Mediana wynikowego rozkładu jest używana jako pozorny współczynnik korygujący.
Poprawione intensywności w każdym kanale są uzyskiwane przez:

PF jest następnie obliczany według:

Różne populacje mogą być oddzielone w wielowymiarowej przestrzeni rozpiętej przez PFs. Położenie i szerokość każdego stanu jest określana przez dopasowanie danych do wielowymiarowych funkcji Gaussa. Późniejsza optymalizacja jednego globalnego HMM w oparciu o wszystkie ślady PF zapewnia ilościowy opis obserwowanej kinetyki. Nawet niewielkie zmiany stawek są wykrywalne.
HMM umożliwiają wnioskowanie o modelu stanu na podstawie zbioru zaszumionych śladów czasu. Uważa się, że system znajduje się w jednym ze zbiorów dyskretnych, ukrytych stanów w danym momencie, a rzeczywista obserwacja (tj. emisja) jest funkcją probabilistyczną tego stanu ukrytego20. W przypadku danych TIRFM smFRET prawdopodobieństwa emisji bi na stan i mogą być modelowane za pomocą ciągłych gaussowskich funkcji gęstości prawdopodobieństwa. W regularnie rozmieszczonych dyskretnych punktach czasowych przejścia z jednego stanu do drugiego mogą zachodzić zgodnie z prawdopodobieństwem przejścia, które jest niezmienne w czasie i zależy tylko od bieżącego stanu. Macierz przejścia A zawiera te prawdopodobieństwa przejścia aij między wszystkimi ukrytymi stanami. Początkowy rozkład stanu
podaje prawdopodobieństwa specyficzne dla stanu
dla pierwszego punktu czasowego śladu czasowego. Korzystając z podejścia maksymalnego prawdopodobieństwa, parametry te można zoptymalizować, aby jak najlepiej opisywać dane za pomocą algorytmów Forward-Backward i Bauma-Welcha20,21. W ten sposób uzyskuje się estymatory maksymalnego prawdopodobieństwa (MLE). Wreszcie, sekwencję stanów, która najprawdopodobniej wytworzyła trajektorię obserwacji, można wywnioskować za pomocą algorytmu Viterbiego. W przeciwieństwie do innych analiz HMM smFRET data24,25,26 nie używamy HMM jako zwykłego "wygładzenia" danych, ale wyodrębniamy model stanu kinetycznego ze zbioru danych bez konieczności dopasowywania histogramów czasu przebywania27. Analiza HMM odbywa się za pomocą własnych skryptów przy użyciu Igor Pro. Implementacja kodu opiera się na reference21. Na naszej stronie internetowej udostępniamy zestaw oprogramowania i przykładowe dane w celu przestrzegania sekcji 5 i 6 niniejszego protokołu (https://www.singlemolecule.uni-freiburg.de/software/3d-fret). Pełne oprogramowanie jest dostępne na życzenie.
Punkty czasowe w danych z PF <-1 lub PF >2 w dowolnym kanale detekcji mają przypisane minimalne prawdopodobieństwo emisji dla wszystkich stanów (10-200). Zapobiega to sztucznym przejściom w tych punktach danych.
Parametry prawdopodobieństwa emisji są uzyskiwane z dopasowania histogramu PF 3D do funkcji Gaussa, jak opisano w kroku 5.7. Parametry te są utrzymywane na stałym poziomie podczas optymalizacji HMM.
W przedstawionym podejściu, wektor rozkładu stanu początkowego i macierz przejścia są używane globalnie do opisania całego zestawu śladów. Są one aktualizowane na podstawie wszystkich cząsteczek N ze zbioru danych zgodnie z reference27.
Parametry początkowe dla początkowego rozkładu stanu są określane na podstawie rzutów 2D histogramu PF (krok 5.3), a prawdopodobieństwa przejścia są ustawiane na 0,05, z wyjątkiem prawdopodobieństwa pozostania w tym samym stanie, które są wybierane tak, że prawdopodobieństwo opuszczenia określonego stanu jest znormalizowane do jedności.
Metoda profilowania prawdopodobieństwa służy do wyznaczania przedziałów ufności (CI) dla wszystkich wskaźników przejścia21,22, które służą jako miarodajne oszacowania ich niepewności. Aby obliczyć granice CI dla określonej stopy, prawdopodobieństwo przejścia odsetek ustala się na wartość inną niż MLE. W ten sposób uzyskuje się model testowy λ'. Test ilorazu prawdopodobieństwa (LR) prawdopodobieństwa
przy zbiorze danych 0 przeprowadza się zgodnie z:

95% pewność graniczna dla parametru jest osiągana, gdy LR przekracza 3,841, 95% kwantyl rozkładu x2 o jednym stopniu swobody22,23.
Moc metody jest demonstrowana za pomocą Hsp90. To obfite białko występuje u bakterii i eukariontów i jest częścią komórkowej odpowiedzi na stres28. Jest to obiecujący cel leku w leczeniu raka29. Hsp90 jest homodimerem z jedną kieszenią wiążącą nukleotyd w domenie N-końcowej każdej podjednostki30. Może przechodzić przejścia między co najmniej dwiema globalnie odrębnymi konformacjami, jedną zamkniętą i jedną N-końcową otwartą, konformacją w kształcie litery V19,31,32. Dimeryczna natura bezpośrednio rodzi pytanie o wzajemne oddziaływanie między dwoma miejscami wiązania nukleotydów w Hsp90.
Poniżej przedstawiamy krok po kroku protokół pozyskiwania i analizy danych z trójkolorowego eksperymentu smFRET na drożdżach Hsp90 i nukleotydach. Analizowane jest zależne od konformacji wiązanie znakowanego fluorescencyjnie AMP-PNP (AMP-PNP*, niehydrolizowanego analogu ATP). Zastosowanie opisanej procedury pozwala na badanie wiązania nukleotydów i jednocześnie zmian konformacyjnych Hsp90, a tym samym ujawnia kooperatywność między dwiema kieszeniami wiążącymi nukleotydy Hsp90.