Aby skutecznie przeanalizować dane eksperymentalne z testów parowania i wypierania, konieczne jest określenie, w jaki sposób zmiana emisji fluorescencji FAM odpowiada konwersji substratów DNA w produkty. W tym celu należy wyznaczyć odpowiedni zakres intensywności fluorescencji. W przypadku testu parowania emisję fluorescencji 16FA(-), odpowiadającą substratowi ssDNA, porównuje się z emisją 16FA(-) zdenaturowanego z 16AR(+)_40bp, co odpowiada końcowym produktom tej reakcji (Rycina 2A). Odpowiada to maksymalnej wydajności FRET, a tym samym maksymalnemu spadkowi intensywności fluorescencji, jakiego należałoby się spodziewać, gdyby cały substrat ssDNA został przekształcony w produkt dsDNA. W przypadku testu wypierania emisję 16FA(-)_40bp zdenaturowanego z 16AR(+)_40bp, odpowiadającą substratowi, porównuje się z emisją 16FA(-)_40bp, która odpowiada produktowi końcowemu (Rycina 2B). W tym przypadku maksymalny wzrost intensywności fluorescencji FAM odzwierciedla scenariusz, w którym cały substrat dsDNA zostaje przekształcony w produkt ssDNA. Rad51 z S. pombe nie wpływał na emisję fluorescencji FAM ani na wydajność gaszenia FAM przez ROX w obu testach (Rycina 2). Maksymalną wydajność FRET należy mierzyć ponownie przy każdym nowym przygotowaniu oligonukleotydów, ponieważ zależy ona od wydajności znakowania oligonukleotydów.
Reprezentatywne dane dotyczące reakcji parowania i wypierania nici DNA przedstawiono na Rysunku 3. Efekty reakcji spontanicznych pomiędzy substratowymi DNA a fotoblaknięciem były niewielkie w obu testach, co potwierdzają pomijalne zmiany w emisji FAM bez udziału Rad51 w porównaniu do znacznych zmian zaobserwowanych z udziałem Rad51 (Rysunek 3A oraz Rysunek 3B). Na podstawie danych przedstawionych na Rysunku 3A lub Rysunku 3B, zmianę fluorescencji przeliczono odpowiednio na zmianę ilości substratu (FP) lub produktu końcowego (FD), stosując równania opisane w krokach 3.2.2 lub 3.3.2 (Rysunek 3C oraz Rysunek 3D).
Reakcję parowania zasymulowano przy użyciu sekwencyjnego modelu reakcji trójetapowej, składającego się z tworzenia pierwszego trójnićiciowego produktu pośredniego (C1), przejścia pierwszego produktu pośredniego w drugi produkt pośredni (z C1 do C2) oraz uwalniania ssDNA z drugiego produktu pośredniego w celu utworzenia dwóch produktów (D + E) (Rysunek 3E). Aby sprawdzić, czy symulacja z wykorzystaniem sekwencyjnego modelu reakcji trójetapowej jest zgodna z danymi eksperymentalnymi, obliczono rezydua między danymi eksperymentalnymi z testu parowania nici DNA a krzywą teoretyczną uzyskaną z symulacji (Rysunek 3F). Dodatkowo obliczono rezydua między reakcją parowania a krzywą teoretyczną wygenerowaną przy użyciu sekwencyjnego modelu reakcji dwuetapowej (Rysunek 3G). Rezydua dla reakcji parowania i modelu dwuetapowego wykazują systematyczne odchylenie w początkowej fazie, natomiast rezydua dla reakcji parowania i modelu trójetapowego nie wykazują takiego odchylenia. Wskazuje to, że model trójetapowy lepiej nadaje się do symulowania reakcji parowania niż model dwuetapowy.
Aby sprawdzić, czy model trzystopniowy jest zgodny z reakcją wypierania, która wykrywa późny etap wymiany nici DNA, wygenerowano teoretyczną krzywą reakcji wypierania, wykorzystując parametry kinetyczne uzyskane z symulacji reakcji parowania przedstawionej na Ryc. 3C, a następnie porównano ją z danymi eksperymentalnymi reakcji wypierania pokazanymi na Ryc. 3D (Ryc. 3H). Krzywa teoretyczna była zgodna z danymi eksperymentalnymi z testu wypierania. Na podstawie tych wyników wyciągnięto wniosek, że symulacja z wykorzystaniem modelu trzystopniowego pozwala na racjonalną ocenę reakcji wymiany nici DNA pośredniczonej przez białko Rad51.
Reprezentatywne dane z reakcji parowania nici DNA z udziałem Rad51 oraz kompleksu Swi5-Sfr1, białka pomocniczego Rad51, przedstawiono w Rysunek 4AKompleks Swi5-Sfr1 silnie stymulował aktywność parowania Rad51. Podobnie jak w przypadku braku Swi5-Sfr1, w obecności Swi5-Sfr1 reakcja parowania lepiej odpowiadała modelowi trzystopniowemu niż modelowi dwustopniowemu (Rysunek 4B). Poprzez symulację reakcji przy użyciu modelu trzyetapowego obliczono stałe równowagi reakcji dla każdego etapu reakcji z kompleksem Swi5-Sfr1 lub bez niego. Stałe równowagi reakcji wykazały, że kompleks Swi5-Sfr1 nie stymuluje pierwszego etapu reakcji (Rysunek 4C, panel a), w którym powstaje trzyniciowy kompleks pośredni, ale który silnie stymuluje przejście C1-C2 (Rysunek 4C, panel b) oraz uwalnianie ssDNA z kompleksu pośredniego C2 (Rycina 4C, panel c).

Rysunek 1: Projekt eksperymentalny analiz parowania i wypierania nici DNA. Schematy analiz parowania (A) i wypierania (B) nici DNA. Żółte kółka reprezentują monomery Rad51. Zielone kółka z literą „F” i czerwone kółka z literą „R” reprezentują odpowiednio fluoresceinę (FAM) i karboksy-X-rodaminę (ROX). Czarne nici DNA są identyczne pod względem sekwencji i komplementarne do niebieskich nici DNA. Cienkie czarne linie z grotami strzałek wskazują nazwę każdego oligonukleotydu, zgodnie z Tabelą 1. Rysunek został zaadaptowany i zmodyfikowany na podstawie pracy Ito et al.19. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 2: Pomiary maksymalnej wydajności FRET w testach parowania i wypierania. (A) Porównanie widm fluorescencji między substratem ssDNA, 16FA(-), a produktem dsDNA, 16FA(-) sparowanym z 16AR(+)_40bp, w teście parowania. Linie niebieska i czerwona reprezentują widma fluorescencji substratu odpowiednio bez i z dodatkiem Rad51. Linie zielona i fioletowa przedstawiają widma fluorescencji produktu końcowego odpowiednio bez i z dodatkiem Rad51. (B) Porównanie widm fluorescencji między substratem dsDNA, 16FA(-)_40bp sparowanym z 16AR(+)_40bp, a produktem ssDNA, 16FA(-)_40bp, w teście wypierania. Linie niebieska i czerwona reprezentują widma fluorescencji produktu końcowego odpowiednio bez i z dodatkiem Rad51. Linie zielona i fioletowa przedstawiają widma fluorescencji substratu odpowiednio bez i z dodatkiem Rad51. Rysunek ten został zaadaptowany i zmodyfikowany na podstawie pracy Ito et al.19. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3: Reakcje parowania i wypierania nici DNA pośredniczone przez Rad51. (A) Przebieg czasowy znormalizowanej fluorescencji reakcji parowania z Rad51 lub bez niego. (B) Przebieg czasowy znormalizowanej fluorescencji reakcji wypierania z Rad51 lub bez niego. (C) Przebieg czasowy zmiany ilości substratu w reakcji parowania z Rad51. (D) Przebieg czasowy zmiany ilości substratu w reakcji wypierania z Rad51. (E) Schemat sekwencyjnego trzystopniowego modelu reakcji. A i B odpowiadają filamentowi presynaptycznemu i donorowemu dsDNA. C1 odpowiada pierwszemu (niedojrzałemu) tronicowemu intermediatowi. C2 odpowiada drugiemu (dojrzałemu) tronicowemu intermediatowi. D i E odpowiadają heterodupleksowi i ssDNA uwolnionemu z C2. (F i G) Residuwa pomiędzy danymi eksperymentalnymi z testu parowania nici DNA a krzywą teoretyczną uzyskaną w drodze symulacji z użyciem modelu trzystopniowego (F) lub dwustopniowego (G). (H) Czerwone kropki wskazują dane eksperymentalne z reakcji wypierania z Rad51 przedstawionej na panelu D. Niebieska linia wskazuje krzywą teoretyczną produktów końcowych. Krzywa teoretyczna została wygenerowana za pomocą symulacji z wykorzystaniem stałych szybkości reakcji uzyskanych z testu parowania przedstawionego na panelu C. Rysunek zaadaptowany i zmodyfikowany z Ito et al.19. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 4: Swi5-Sfr1 stymuluje drugi i trzeci etap reakcji wymiany nici DNA. (A) Przebieg czasowy zmiany ilości substratu w reakcji parowania z dodatkiem lub bez dodatku Swi5-Sfr1 (S5S1). (B) Reszty między danymi eksperymentalnymi z testu parowania nici DNA z Swi5-Sfr1 a krzywą teoretyczną uzyskaną w drodze symulacji z wykorzystaniem modelu dwuetapowego (linia niebieska) lub trzyetapowego (linia czerwona). (C) Reakcja parowania przedstawiona na Rysunku 4A została zasymulowana za pomocą modelu trzyetapowego z wykorzystaniem programu analitycznego23 (patrz Tabela Materiałów). Stałe równowagi każdego etapu reakcji, K1 (a), K2 (b) i K3 (c), zostały wyznaczone z symulacji. Rysunek został zaadaptowany i zmodyfikowany z pracy Ito et al.19. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
| Oligonukleotydy do testu parowania nici DNA |
| 16FA(-) | 5’-[FAM]-AAATGAACATAAAGTAAATAAGTATAAGGATAATACA
AAATAAGTAAATGAATAAACATAGAAAATAAAGTAAAGGATAT AAA -3’ |
| 16A(-)_40bp | 5’-AAATGAACATAAAGTAAATAAGTATAAGGATAATACAAAA-3’ |
| 16AR(+)_40bp | 5’-TTTTGTATTATCCTTATACTTATTTACTTTATGTTCATTT-[ROX]-3’ |
| Oligonukleotydy do testu wypierania nici DNA |
| 16A(-) | 5’-AAATGAACATAAAGTAAATAAGTATAAGGATAATACAAAATA
AGTAAATGAATAAACATAGAAAATAAAGTAAAGGATAT AAA -3’ |
| 16FA(-)_40bp | 5’-[FAM]-AAATGAACATAAAGTAAATAAGTATAAGGATAATACAAAA-3’ |
| 16AR(+)_40bp | 5’-TTTTGTATTATCCTTATACTTATTTACTTTATGTTCATTT-[ROX]-3’ |
Tabela 1: Lista oligonukleotydów użytych w analizach parowania i wypierania nici DNA. W odpowiednich przypadkach pozycje fluoroforów (fluoresceina amidydowa, FAM; karboksy-x-rodamina, ROX) zostały wskazane w nawiasach kwadratowych.