Aby zweryfikować poprawność proponowanego protokołu, zaprezentowane tutaj eksperymenty PD przeprowadzono z wykorzystaniem biotynilowanego aptameru RNA, zaprojektowanego in silico w celu specyficznego wiązania TDP-4320. Ten RNA wiąże swój białkowy cel z wysokim powinowactwem (Kd = 90 nM)20. W niniejszej pracy ten RNA o sekwencji 5'-CGGUGUUGCU-3' określany jest nazwą „+RNA”. Jako kontrolę negatywną wykorzystano sekwencję odwrotnie komplementarną do +RNA, nazwaną tutaj „-RNA”. Jej sekwencja to 5'-AGCAACACCG-3'. -RNA wykazuje znacznie niższe powinowactwo wiązania do TDP-43 (Kd = 1.5 µM)19. Na potrzeby opisanego tutaj protokołu oligonukleotydy RNA zakupiono w formie sprzężonej z cząsteczką biotyny, aby umożliwić ich wiązanie z kulkami streptawidynowymi. +RNA zakupiono z biotyną-TEG na końcu 3', która zawiera 15-atomowy łącznik z trietylenoglikolu pomiędzy biotyną a grupą fosforanową kwasu nukleinowego; -RNA posiadał natomiast biotynę na końcu 5', sprzężoną z kwasem nukleinowym via łącznik amino-C6. Jeżeli jednak projekt przynęty RNA jest odporny i o ile nie występuje interferencja strukturalna lub chemiczna pomiędzy łącznikiem a RNA, można zastosować inne pozycje koniugacji biotyny oraz inne długości łączników.
Znajomość tożsamości głównego białka, które powinno być związane z sondą +RNA po procesie PD, umożliwiła walidację protokołu poprzez identyfikację TDP-43 w eluate, z wykorzystaniem zarówno spektrometrii mas (MS), jak i western blotu (WB) (Rysunek 1).
Analizę MS przeprowadzono na czterech powtórzeniach PD wykonanych z użyciem +RNA lub -RNA (Rysunek 2). Identyfikacja interaktomów +RNA i -RNA wykracza poza zakres niniejszego protokołu, jednak przedstawiono niektóre wyniki potwierdzające dokładność procedury. Warto zauważyć, że naniesienie istotnie wzbogaconych białek na wykres wulkaniczny (volcano plot) wykazało, że całkowita zawartość białek oraz białka wzbogacone elucyjne z +RNA były znacznie wyższe niż te odzyskane z -RNA (Rysunek 2). Oznacza to, że mimo tej samej długości i struktury (liniowej), +RNA może tworzyć większą liczbę specyficznych oddziaływań, które są utrzymywane aż do etapu elucji w wysokim stężeniu soli. Prawdopodobnie -RNA tworzy natomiast większą liczbę niespecyficznych kontaktów, które są przerywane podczas etapów przemywania. Zgodnie z oczekiwaniami, TDP-43 został zidentyfikowany jako unikalny interaktor +RNA20; średnia kwantyfikacja bez znakowania (LFQ) dla czterech powtórzeń PD wykonanych z +RNA wynosi 31,96 ± 0,56, podczas gdy białko to nie zostało zidentyfikowane wśród interaktorów -RNA. Ponadto, spośród wszystkich unikalnych interaktorów +RNA, TDP-43 okazał się białkiem o największym stopniu wzbogacenia.
W celu dalszej walidacji protokołu wykorzystano opracowany wewnętrznie algorytm catRAPID18,19 do obliczeniowego przewidywania, które inne białka będą wiązać specyficznie +RNA lub -RNA. W szczególności, za pomocą funkcji „interaction propensity” algorytmu catRAPID, zdefiniowanej w naszej poprzedniej pracy27, obliczono wskaźniki interakcji +RNA i -RNA z białkami tworzącymi ludzki proteom. Wśród białek ocenionych z wysoką pewnością przewidziano, że HNRNPH3 będzie wiązać selektywnie +RNA (wskaźnik interakcji z +RNA = 1,01; wskaźnik interakcji z -RNA = 0,21), a PCBP2 będzie oddziaływać specyficznie z -RNA (wskaźnik interakcji z +RNA = -0,5; wskaźnik interakcji z -RNA = 0,31) (Rysunek 3A). Ponadto przewidziano, że białko RBM41 będzie wykazywać niespecyficzne powinowactwo do obu oligonukleotydów RNA (wskaźnik interakcji z +RNA = 0,4; wskaźnik interakcji z -RNA = 0,39) (Rysunek 3A). Analiza MS rzeczywiście potwierdziła obecność HNRNPH3 i PCBP2 odpowiednio w PD dla +RNA i -RNA, natomiast stwierdzono, że RBM41 oddziałuje z obydwoma (Rysunek 3B).
Metodę WB wykorzystano do wykrycia obecności TDP-43 w celu dalszego potwierdzenia wyników oraz podczas optymalizacji protokołu (Rycina 4). W opisanym tutaj postępowaniu pobrano różne próbki na różnych etapach. Próbka wejściowa (IN) składała się z białek całkowitych rozcieńczonych w buforze do lizy. Frakcję przepływową (FT) otrzymano po całonocnej inkubacji białek całkowitych z paciorkami streptawidyny uprzednio powleczonymi biotynilowanym RNA, co reprezentowało frakcję białek, które nie związały się z RNA. Ostatecznie eluat (EL) zawierał wszystkie białka, które specyficznie rozpoznały badane RNA, ponieważ między etapem FT a EL trzy kroki płukania roztworem 150 mM soli i 0,1% triton-X powinny usunąć najsłabsze oddziaływania.
Dla każdej powtórki taką samą ilość (5% v/v) frakcji IN, FT i EL poddano równoległej analizie metodą SDS-PAGE i zabarwiono przeciwciałem anty-TDP-43 (Rycina 4). W przypadku +RNA prążek TDP-43 zaobserwowano w IN oraz w EL, co wskazuje, że białko obecne od początku w całkowitym ekstrakcie białkowym jest zatrzymywane przez +RNA podczas etapów płukania i eluowane dopiero na końcu za pomocą buforu o wysokim stężeniu soli. TDP-43 było również obecne w IN dla -RNA, jednak prążek odpowiadający temu białku jest widoczny także w FT, co wskazuje, że to RNA nie wiąże TDP-43. Brak prążka TDP-43 w EL potwierdza ten wynik.
Podczas optymalizacji protokołu badano elucję białek specyficznie związanych z sekwencjami RNA, wykorzystując zarówno bufor elucyjny zawierający 1 M NaCl (EB1), jak i bufor elucyjny z 2 M NaCl (EB2) (Rysunek 4). Eluaty otrzymane przy użyciu obu buforów EB porównano za pomocą SDS-PAGE, a następnie przeprowadzono blotting z użyciem przeciwciała anty-TDP-43. Uzyskane obrazy przeanalizowano programem ImageJ28, aby określić ilościowo różnice w ilości TDP-43 eluowanego przez oba bufory. Ogólnie nie zaobserwowano istotnych różnic i wyciągnięto wniosek, że w tych testach stężenie soli 1 M jest wystarczające do przerwania nawet najsilniejszych oddziaływań białko-RNA.
Ogólnie rzecz biorąc, przedstawione tutaj wyniki MS i WBs wykazują, że protokół ten jest wydajny w specyficznym wychwytywaniu białek oddziałujących z danym RNA i umożliwia elucję w buforach kompatybilnych z dalszymi analizami.

Rycina 1: Schemat procedury eksperymentalnej zastosowanej w proponowanym protokole. (A) Biotynilowany oligonukleotyd RNA jest przygotowywany w buforze do lizy w odpowiednim stężeniu. (B) Magnetyczne kulki streptawidynowe są przemywane, blokowane za pomocą tRNA drożdży i obciążane biotynilowanym RNA. (C) Do mieszaniny kulek i RNA dodaje się całkowity ekstrakt białkowy otrzymany z hodowlanych linii komórkowych ssaków. (D) Przeprowadza się wielokrotne przemywanie w celu usunięcia nieswoistych oddziaływań. (E) Swoiste białka oddziałujące są oddzielane od RNA za pomocą roztworu hipertonicznego. (F) Tożsamość interaktorów jest ujawniana za pomocą spektrometrii mas, a wybrane przypadki są walidowane za pomocą western blot. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2Strategia analityczna dla bezetykietowej kwantyfikacji białek metodą MS. (AEluowane białka wytrąca się w zimnym acetonie przez noc. Następnie białka poddaje się denaturacji i przeprowadza trawienie w roztworze. Peptydy proteolityczne są zagęszczane i odsalane.B) Peptydy są analizowane poprzez LC-MS/MS z zastosowaniem „metody shotgun”. (CPrzetwarzanie i analiza danych surowych są wykonywane odpowiednio przy użyciu oprogramowania MaxQuant i Perseus.D) Białka o statystycznie istotnym wzbogaceniu przedstawiono na wykresie wulkanicznym. Kliknij tutaj, aby wyświetlić większą wersję tej ilustracji.

Rysunek 3: Korelacja między przewidywaną skłonnością do interakcji a interakcjami +RNA i -RNA określonymi eksperymentalnie. (A) Wyniki interakcji catRAPID w odniesieniu do HNRNPH3, PCBP2 i RBM41, wskazujące na preferencyjne wiązanie +RNA przez HNRNPH3 i -RNA przez PCBP2, podczas gdy w przypadku RBM41 przewiduje się nieselektywne wiązanie obu sekwencji RNA. (B) Średnie wartości kwantyfikacji bezznakowej (LFQ) wyznaczone za pomocą analizy spektrometrii mas z wyciągów (pull-downs) przeprowadzonych z użyciem +RNA i -RNA. Analiza potwierdza, że HNRNPH3 wiąże wyłącznie +RNA, PCBP2 wiąże wyłącznie -RNA, a RBM41 wiąże oba w równym stopniu. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 4: Walidacja metodą Western blot obecności/braku TDP-43 wśród interaktorów wybranych sekwencji RNA. Membrana WB została poddana działaniu przeciwciała anty-TDP-43. IN = input; FT = flow-through; EL (EB1) = elucja z użyciem buforu elucyjnego 1; EL (EB2) = elucja z użyciem buforu elucyjnego 2; znak „+” oznacza próbki pochodzące z pull-down wykonanego z +RNA; znak „-” oznacza próbki pochodzące z pull-down wykonanego z -RNA; ścieżka 1 zawiera drabinę białkową. TDP-43 zostało wskazane strzałką. WB wskazuje, że TDP-43 znajduje się wśród interaktorów +RNA, ale nie wśród interaktorów -RNA. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
| Nazwa buforu | Skład | |
| 10x bufor do transferu | 250 mM tris, 1.92 M glicyny, 1% SDS, 20% metanolu. Rozcieńczyć 10-krotnie przed użyciem | PD |
| 20X bufor do elektroforezy MES SDS | 1 M MES, 1 M tris, 2% SDS, 20 mM EDTA. Dostosować pH do 7.3. Rozcieńczyć 20-krotnie przed użyciem |
| 4x bufor do nakładania próbek | 0.25 M Tris (zasada), 0.28 M SDS, 40% glicerolu, 20% 2-merkaptoetanolu, 4 mg/ml błękitu bromofenolowego |
| Bufor elucyjny 1 | 20 mM fosforanu pH 7.5, 1 M NaCl, 0.5 mM EDTA, 0.1% Triton X-100, 1 mM DTT (dodać po kwantyfikacji) |
| Bufor elucyjny 2 | 20 mM fosforanu pH 7.5, 2 M NaCl, 0.5 mM EDTA, 0.1% Triton X-100, 1 mM DTT (dodać po kwantyfikacji) |
| Bufor do lizy | 10 mM Tris-HCl pH 7.4, 150 mM NaCl, 0.5 mM EDTA, 0.1% Triton X-100, 1 mM DTT oraz inhibitory proteaz |
| Bufor solny Tris z Tween-20 | 1 M Tris-HCl pH 7.4, 3 M NaCl, 2.0% Tween-20 |
| Bufor piorący 1 | 10 mM Tris-HCl pH 7.4, 150 mM NaCl, 0.5 mM EDTA, 0.1% Triton TM X-100, 1 mM DTT oraz inhibitory proteaz |
| Bufor piorący 2 | 25 mM Hepes pH 8, 150 mM NaCl, 0.5 mM EDTA, 0.1% Triton X-100, 1 mM DTT oraz inhibitory proteaz |
| Bufor A | 0.1% kwasu mrówkowego | MS |
| Bufor B | 60% acetonitrylu, 0.1% kwasu mrówkowego |
| Bufor do denaturacji | 8M mocznika, 50 mM Tris-HCl |
Tabela 1: Bufor PD i MS. Nazwy i skład buforów stosowanych w eksperymentach pull-down (PD) lub w analizie spektrometrią mas (MS).