Schemat strategii powinowactwa do oczyszczania i korelacyjnego profilowania białek w warunkach natywnych z wykorzystaniem niebieskiego barwnika i spektrometrii mas (ABC-MS) przedstawiono na Rysunku 7. Kompleksy białek natywnych wokół białka docelowego są izolowane poprzez oczyszczanie powinowactwa z użyciem przeciwciał przeciwko znacznikowi epitopowemu (w tym przypadku FLAG) oraz elucję konkurencyjną. Kompleksy są rozdzielane za pomocą BN-PAGE, a cała ścieżka żelu jest wycinana na 48 sekcji i przygotowywana do analizy shotgun LC-MS/MS. Ilościowe dane MS są wykorzystywane do wygenerowania profilu migracji dla każdego zidentyfikowanego białka w trakcie rozdzielania niebieskim barwnikiem. Białka, które oddziałują ze sobą, tworząc odrębny kompleks, wykazują podobne profile migracji z nakładającymi się pikami. Gdy białko docelowe bierze udział w więcej niż jednym zgromadzeniu, w jego profilu migracji obserwuje się wiele pików, pod warunkiem że subkompleksy mieszczą się w zakresie rozdzielczości żelu BN-PAGE. Systematyczne porównanie profili migracji można przeprowadzić za pomocą korelacyjnego profilowania białek z wykorzystaniem grupowania hierarchicznego. Dendrogram oraz intensywność pików dla wszystkich frakcji są wizualizowane na mapie ciepła (heat-map), co ułatwia identyfikację oddziałujących białek należących do odrębnych kompleksów białkowych (Rysunek 8).
Strategię tę zastosowano do analizy partnerów oddziałujących z Mta2, kluczową podjednostką kompleksu remodelowania chromatyny NuRD 20. Jak pokazano na Rysunku 9, profil migracji Mta2 wykazał dwa piki o różnych intensywnościach w zakresie od 700 kDa do 1,2 MDa, przy czym pik o niższej masie wykazywał większą obfitość. Inne kluczowe podjednostki NuRD, w tym Mta1/3, Hdac1/2 i Mbd3, wykazały identyczny wzorzec rozdzielenia, choć piki dla niektórych podjednostek, mianowicie Chd4, Gatad2a/b i Rbbp4/7, miały odwrzesny rozkład obfitości (Rysunek 9A, B). W przeciwieństwie do nich, profil Cdk2ap1, czynnika regulacyjnego rekrutującego kompleks NuRD do promotorów genów Wnt 23, wykazał jedynie pik o wyższej masie (Rysunek 9C), podobnie jak Sall4, represor transkrypcyjny, o którym wykazano, że wiąże się z NuRD 20,24,25. Zatem frakcjonowanie białek związanych z Mta2, oczyszczonych metodą powinowactwa, za pomocą blue native PAGE pozwoliło na rozdzielenie dwóch różnych form kompleksu NuRD.
Metoda ABC-MS umożliwia również identyfikację nowych interaktorów przy jednoczesnym przypisaniu ich do konkretnych białek. Poprzez analizę klastrów białek oraz intensywności frakcji przedstawionych na mapie ciepła wykryliśmy silną korelację pomiędzy podjednostkami NuRD a Wdr5, podjednostką regulacyjną kompleksu metylotransferazy MLL 26; Wdr5 wykazywał dwa piki migracji pokrywające się z dwoma pikami NuRD (Rysunek 8), co sugeruje nową interakcję między Wdr5 a NuRD. Potwierdziliśmy tę interakcję za pomocą ko-immunoprecypitacji oraz ko-migracji w chromatografii wykluczania cząsteczek 20.
Wybór metody obliczania metryki odległości w klastrowaniu hierarchicznym wpłynie na kształt klastrów, a w konsekwencji na korelacje. Zalecamy eksperymentowanie z różnymi metrykami, aby uzyskać najlepsze dopasowanie do istniejącej wiedzy na temat oddziaływań białka lub kompleksu, który jest przedmiotem badań. Zazwyczaj najlepsze wyniki osiągaliśmy przy zastosowaniu metryki odległości Manhattan (L1) 20. W literaturze zgłaszano również wykorzystanie korelacji Pearsona lub metryki odległości euklidesowej do identyfikacji kompleksów w oparciu o alternatywne techniki frakcjonowania 10,11,12.
Ilościowe dane MS uzyskane z frakcji Blue Native mogą być również wykorzystane do określenia stechiometrii kompleksów białkowych. Wartość iBAQ (ilościowa analiza oparta na intensywności) zapewnia pomiar względnej obfitości zidentyfikowanych białek 27,28. Wartości iBAQ dla członków kompleksu białkowego we wszystkich frakcjach w obrębie piku profilu migracji są normalizowane do wartości białka przynęci (bait protein) w celu uzyskania ilości względnych. Znormalizowane wartości iBAQ w obrębie piku profilu dla danego członka kompleksu białkowego powinny wykazywać trend poziomy, a wartości tego trendu odzwierciedlają stechiometrię oddziałujących białek w stosunku do białka przynęci. Bardziej szczegółną prezentację obliczeń stechiometrii przedstawiono w 20.

Rysunek 7: Schematyczny przebieg pracy podsumowujący strategię ABC-MS. Rysunek ten został zmodyfikowany na podstawie 20. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 8. Grupowanie hierarchiczne profili migracji BN-PAGE interaktorów Mta2. Białka związane z Mta2 zostały pogrupowane na podstawie podobieństwa ich profili migracji. Przedstawiono jedynie fragment mapy ciepła zawierający kompleks NuRD (obramowany niebieską ramką). Żółta ramka wyróżnia silną korelację Wdr5 z kompleksem NuRD. Adnotowane masy cząsteczkowe oszacowano na podstawie odległości migracji standardów białkowych przeprowadzonych w tym samym żelu. Badanie to zostało pierwotnie opublikowane w Molecular and Cellular Proteomics 20 American Society for Biochemistry and Molecular Biology. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 9. Profil migracji NuRD i białek powiązanych w BN-PAGE. A) Mta2 i podjednostki NuRD wykazujące podobny wzór intensywności. B) Mta2 i podjednostki NuRD o odwrotnym wzorze intensywności. C) Mta2 i Cdk2ap1, który występuje tylko w kompleksie NuRD o wyższej masie cząsteczkowej. Adnotowane masy cząsteczkowe zostały oszacowane na podstawie profilu migracji standardów białkowych. Badanie to zostało pierwotnie opublikowane w Molecular and Cellular Proteomics 20 American Society for Biochemistry and Molecular Biology. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Dane uzupełniające: Przykładowy zestaw danych. Plik proteingroups.txt wygenerowany przez program MaxQuant, zawierający identyfikacje białek i wartości ilościowe z eksperymentu ABC-MS. Badanie to zostało pierwotnie opublikowane w czasopiśmie Molecular and Cellular Proteomics 20 Amerykańskiego Towarzystwa Biochemii i Biologii Molekularnej. Kliknij tutaj, aby pobrać ten plik.