$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Rygorystyczność warunków była początkowo utrzymywana na niskim poziomie, aby zachować cząsteczki wiążące występujące w niskiej liczbie kopii w pierwszych kilku rundach. W pierwszej rundzie zastosowano w szczególności najniższy rygor, aby zachować białka wiążące występujące zazwyczaj jako unikalne sekwencje, co przejawiało się wysokimi stężeniami kortyzolu i DNA oraz brakiem etapów selekcji negatywnej. Sukces tej selekcji aptamerów jest wykazany przez ewolucję pul z niskiego procenta elucji przez cel (runda 2) do wyższej frakcji elucji przez cel, która pozostaje stosunkowo stała w rundach 12-13 (Tabela 1). Taki poziom plateau DNA elucjowanego przez cel jest tradycyjnie punktem końcowym selekcji („wzbogacenia”). Jednak w tej metodzie rygorystyczność selekcji po wzbogaceniu została dodatkowo zwiększona poprzez wydłużenie sondy wychwytującej cDNA (Rycina 1). Wydłużenie tej sondy zwiększa siłę hybrydyzacji pomiędzy cDNA a pulą, co podnosi temperaturę topnienia (Tm) kompleksu i wymaga silniejszej interakcji między celem a sekwencją, aby uwolnić sekwencje wiążące do roztworu. W pierwszej rundzie selekcji zastosowano słabszą interakcję cDNA ze względu na jedną bazę G mniej na końcu 5´ biblioteki. Jest to zgodne z ogólnie niższym rygorem warunków w pierwszej rundzie i nie oczekuje się, aby miało to znaczący wpływ na selekcję, poza umożliwieniem przejścia do następnej rundy większej liczby potencjalnych cząsteczek wiążących, a także sekwencji tła (dehybrydyzacja oparta na termodynamice). Sekwencje tła te były minimalizowane w miarę postępu selekcji w stronę wzbogacenia poprzez stosowanie bardziej rygorystycznych warunków w kolejnych rundach selekcji. Gdy cDNA wydłużono z 7-meru do 8-meru w rundzie 14, procent DNA elucjowanego przez cel spadł z 15,3% do 11,1% przy wzroście Tm z 19,2 °C do 26,7 °C. W rundzie 15 cDNA przedłużono dalej do 9-meru z Tm wynoszącą 31,3 °C, co obniżyło całkowitą ilość elucjowanego materiału do 3,3%.
Dalsze informacje na temat wzbogacenia można uzyskać poprzez sekwencjonowanie kilku pul, zarówno przykładów wzbogaconych, jak i niewzbogaconych, aby śledzić postęp pul w każdej rundzie. Sekwencjonowanie nowej generacji (NGS) stało się potężną metodą sekwencjonowania w procesach selekcji, głównie ze względu na uzyskanie większej liczby odczytów sekwencji (całkowita liczba raportowanych sekwencji) w porównaniu z sekwencjonowaniem metodą Sangera. Wyższa liczba odczytów zapewnia większe pokrycie całkowitej liczby sekwencji w puli, dając pełniejszy obraz różnorodności sekwencji. NGS eliminuje również błąd klonowania32 i pozwala na sekwencjonowanie wielu pul w jednej partii dzięki zastosowaniu kodów kreskowych (bar codes)24,33. NGS wykorzystano do analizy postępu wzbogacenia z trzech oddzielnych pul w tej selekcji próbek (Rycina 2). Jako reprezentatywne punkty selekcji wybrano pule z rundy 6 (lekkie wzbogacenie w porównaniu z rundą 2 według % elucji), rundy 13 (najwyższe wzbogacenie) oraz rundy 15 (najwyższa rygorystyczność). Liczby kopii każdej unikalnej sekwencji przedstawiono jako procent całkowitej liczby odczytów sekwencji dla każdej puli. Sekwencja najwyżej sklasyfikowana (o najwyższej liczbie kopii) stanowiła zaledwie 0,1% wszystkich sekwencji w rundzie 6 (łącznie 33 435 sekwencji, 22 463 unikalne sekwencje). W miarę jak pula staje się maksymalnie wzbogacona w rundzie 13, udział najwyżej sklasyfikowanej sekwencji wzrasta do 15,4% całej puli (łącznie 17 681 sekwencji, 2 987 unikalne sekwencje), co jest wartością 150-krotnie wyższą niż w rundzie 6, pomimo jedynie ok. 5-krotnego wzrostu wzbogacenia obserwowanego w UV. W rundzie 15 (łącznie 28 919 odczytów sekwencji, 2 980 unikalne sekwencje) udział pojedynczej, najwyżej sklasyfikowanej sekwencji skoczył do 44,9% wszystkich sekwencji, mimo że monitorowanie wzbogacenia w UV wykazało, iż ilość elucjowana przez cel była ok. 5x niższa niż w rundzie 13 (Tabela 1). Wzbogacenie jest również widoczne w niższym procencie unikalnych sekwencji (od 67% w rundzie 6 do 10% w rundzie 15), w miarę jak pula ewoluuje z wielu do kilku sekwencji o potencjalnym powinowactwie do celu. Dodatkowo, najwyżej sklasyfikowana sekwencja w rundzie 13 (15-3) była sekwencją o trzeciej najwyższej liczbie kopii w rundzie 15 po wydłużeniu cDNA, a sekwencja sklasyfikowana na drugim miejscu w rundzie 13 stała się sekwencją o najwyższej liczbie kopii w rundzie 15 (15-1):
15-1 GGAATGGATCCACATCCATGGATGGGCAATGCGGGGTGGAGAATGGTTGCCGCACTTCGGCTT
CACTGCAGACTTGACGAAGCTT
15-3 GGAATGGATCCACATCCATGGGAGGGTTGGAAGGGAGGGGCCCGGGGTGGGCCATCGTTCGTT
CACTGCAGACTTGACGAAGCTT
Poprzednie badania analizowały wiązanie tych dwóch sekwencji za pomocą mikrotermeforezy oraz dializy równowagowej24. Wyniki wykazały znaczące wiązanie kortyzolu dla sekwencji 15-1 (Kd = 6.9 ± 2.8 µM w dializie równowagowej; 16.1 ± 0.6 µM w mikrotermeforezie), podczas gdy między kortyzolem a sekwencją 15-3 zaobserwowano minimalne wiązanie. Wskazuje to, że zwiększona długość sondy przechwytującej w 15. rundzie pomogła w zapewnieniu warunków wyższej rygorystyczności w celu wyodrębnienia ligandu o wyższym powinowactwie. Ponadto żadna z sekwencji nie wykazała zauważalnego wiązania z cząsteczką selekcji negatywnej, progesteronem, co dowodzi, że wprowadzenie etapów selekcji negatywnej było korzystne dla procedury.
Fakt, że po zastosowaniu warunków o zwiększonej rygorystyczności w kolejnych rundach selekcji pojawił się ligand o wyższym powinowactwie po maksymalnym wzbogaceniu (określonym tradycyjną metodą analizy procenta DNA elukowanego przez cel), potwierdza zasadność metody wydłużania sondy cDNA do wychwytu po etapie wzbogacenia. Wynik ten pokazuje, że liczba kopii w zsekwencjonowanej puli oraz powinowactwo do celu mogą być powiązane jedynie w określonych warunkach24,34, co stoi w sprzeczności z wcześniejszym doniesieniem sugerującym bezpośrednią korelację35. Podkreśla to również zaletę monitorowania wzbogacenia za pomocą NGS, a nie wyłącznie strategią obliczania % elukowanego, powszechną w większości selekcji, która może naturalnie zmieniać się z rundy na rundę w miarę stosowania warunków o większej rygorystyczności. Niemniej jednak, % elukowanego jest użyteczny do monitorowania wzbogacenia, gdy w kilku rundach stosuje się podobne warunki. Na przykład rundy 6-10 odbywały się w podobnych warunkach, bez istotnej zmiany w stopniu wzbogacenia (Tabela 1). Gdy obserwuje się to przez kilka rund, warunki są prawdopodobnie zbyt rygorystyczne, aby promować wzbogacenie, i badacz powinien zmniejszyć rygorystyczność w następnej rundzie. Przykładowo, stężenie kortyzolu zwiększono z 35 µM do 100 µM w rundach 11-13, a % elukowanego wzrósł z 2,5% w rundzie 10 do 14,5% w rundzie 12. Zatem % elukowanego może służyć jako wytyczna do dostosowywania rygorystyczności w trakcie selekcji, gdy porównuje się podobne warunki między rundami, i może dostarczać informacji uzupełniających do NGS w celu wyznaczenia punktu końcowego selekcji.
Sekwencję 15-1 zastosowano w teście AuNP, aby sprawdzić, czy sekwencja wybrana w sposób mający na celu stworzenie aptameru przełączającego strukturę będzie funkcjonować w teście opierającym się na zmianie strukturalnej po związaniu celu w celu wywołania odpowiedzi. Poprzednie wstępne badania wykazały, że test AuNP z zastosowaniem 15-1 reagował na biomarker stresu, kortyzol, ale nie na dwa inne markery stresu, epinefrynę i norepinefrynę, ani na kwas choliczny, marker choroby wątroby strukturalnie podobny do kortyzolu24. Odpowiedź zaobserwowano w normalnym zakresie wolnego kortyzolu raportowanego w ludzkiej surowicy (~150-600 nM)6, co potwierdziło, że aptamer wybrany pod kątem zmiany strukturalnej po związaniu kortyzolu wywołuje odpowiedź w teście AuNP. W obecnej pracy charakterystyka systemu została rozwinięta poprzez dostosowanie pokrycia (gęstości) 15-1 na powierzchni AuNP. Przy zastosowaniu tego samego zestawu warunków pokrycie DNA zwiększono z 73 D/NP (cząsteczek DNA/AuNP) do 120 D/NP i 200 D/NP (Rysunek 3). Kwas choliczny wywołał minimalną odpowiedź, z wyjątkiem stężenia 10 µM przy 200 D/NP. Odpowiedź testu malała wraz ze wzrostem pokrycia, podobnie jak granice detekcji (LOD). LOD wynosiła 29,5 nM dla 73 D/NP, 145,2 nM dla 120 D/NP i 27,3 µM dla 200 D/NP. Wynik ten jest zgodny z pracą Smith et al., którzy wykazali, że odpowiedź testu AuNP z wykorzystaniem aptameru kokainowego malała, gdy pokrycie aptameru zwiększono z 60 D/NP do 300 D/NP36. Sugeruje to, że zakres detekcji interesującego celu może być dostosowany poprzez optymalizację pokrycia powierzchni DNA w teście AuNP. Mogłoby to być korzystne na przykład przy porównywaniu zakresu wolnego kortyzolu w ludzkiej surowicy (~150-600 nM) i w ludzkiej ślinie (~5-25 nM)6,37. Chociaż płyny fizjologiczne są znacznie bardziej złożone niż ten test buforowy, wyniki stanowią rozszerzenie prac dowodzących zasady działania, wykazując, że warunki AuNP można dostosować w zależności od zastosowania oraz że dalsze prace nad rozszerzeniem medium na płyny biologiczne byłyby interesujące dla środowiska zajmującego się biosensorami.

Rycina 1. Przegląd metody selekcji z przełączaniem struktury. Mała biotynowana sonda przechwytująca cDNA jest unieruchomiona na kuleczkach magnetycznych powlekanych streptawidyną. cDNA jest komplementarne do regionu 5´ biblioteki, co powoduje hybrydyzację biblioteki z kuleczkami. Po dodaniu celu następuje zmiana konformacyjna, która uwalnia sekwencje wiążące z kuleczki, umożliwiając frakcjonowanie przy użyciu magnesu. Nadpływ jest zatrzymywany w celu monitorowania wzbogacenia (% DNA eluowanego przez cel) metodą UV po dializie celu z DNA. Ten krok jest konieczny tylko wtedy, gdy cząsteczka celu absorbuje w tym samym zakresie UV co DNA. Następnie sekwencje są amplifikowane metodą PCR i przygotowywane do kolejnej rundy selekcji. W miarę postępu selekcji stosowana jest selekcja negatywna, w której sekwencje eluowane przez cząsteczkę selekcji negatywnej są odrzucane, a następnie kuleczki z potencjalnymi binderami celu są inkubowane z celem. Długość sondy cDNA jest zwiększana po osiągnięciu maksymalnego wzbogacenia (% eluowanych), aby zwiększyć rygorystyczność selekcji. Rycina ta została przedrukowana za zgodą24. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Runda | [Kortyzol] (µM) | [Progesteron] (µM) | Związane DNA (pmol) | % eluego przez cel | Długość cDNA |
| 1 | 100 | 0 | 357.57 | N/D | 7-mer |
| 2 | 50 | 0 | 145.49 | 1.4 | 7-mer |
| 3 | 50 | 1 | 188.74 | N/D | 7-mer |
| 4 | 50 | 5 | 336.4 | 2.3 | 7-mer |
| 5 | 35 | 5 | 88.05 | N/D | 7-mer |
| 6 | 35 | 10 | 303.89 | 3.1 | 7-mer |
| 7 | 35 | 10 | 255.01 | N/D | 7-mer |
| 8 | 35 | 10 | 236.81 | 2.1 | 7-mer |
| 9 | 35 | 10 | 318.95 | 2.6 | 7-mer |
| 10 | 35 | 10 | 297.18 | 2.5 | 7-mer |
| 11 | 100 | 10 | 147.61 | N/D | 7-mer |
| 12 | 100 | 10 | 154.36 | 14.5 | 7-mer |
| 13 | 100 | 10 | 150.39 | 15.3 | 7-mer |
| 14 | 75 | 20 | 247.71 | 11.1 | 8-mer |
| 15 | 50 | 40 | 409.63 | 3.3 | 9-mer |
Tabela 1. Postęp selekcji i wzbogacenie według % elucji24. Oznaczenie N/A (nie dotyczy) pojawia się przy rundach, w których monitorowanie wzbogacania metodą dializy/UV nie było przeprowadzane w wczesnych etapach selekcji, aby ograniczyć utratę sekwencji o niskiej liczbie kopii. Tabela została przedrukowana za zgodą24.

Rysunek 2. Wzbogacenie selekcji określone przez sekwencjonowanie nowej generacji. (A) Runda 6; (B) Runda 13; (C) Runda 15. Sekwencje zostały uszeregowane według liczby kopii. Udział sekwencji o najwyższej liczbie kopii wzrósł z niskiego procentu całej zsekwencjonowanej puli w rundzie 6 (0,1%) do wysokiego procentu w rundzie 15 (44,9%), w miarę jak pula była wzbogacana o sekwencje wiążące kortyzol. Rysunek przedrukowano za zgodą24. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rycina 3. Odpowiedź testu AuNP przy zmiennym stopniu pokrycia aptamerem 15-1. Gęstość aptameru inkubowanego z AuNP zwiększano z 73 D/NP (czerwony trójkąt) do 120 D/NP (zielony kwadrat) i 200 D/NP (niebieskie koło). Odpowiedzi na kortyzol zaznaczono ciemnymi kolorami, a odpowiedź na kwas cholowy jasnymi kolorami. Odpowiedź testu dla kortyzolu uległa poprawie przy niższych gęstościach, co obniżyło LOD oraz zakres wykrywalności celu. Warunki zapewniające najwyższą odpowiedź na kortyzol (73 D/NP) zostały dalej scharakteryzowane w zakresie liniowym, aby uzyskać więcej punktów w oczekiwanym normalnym zakresie wolnego kortyzolu raportowanym w surowicy ludzkiej (~150-500 nM) i ślinie (5-25 nM)6,37. Minimalna odpowiedź na kwas cholowy przy zastosowaniu tych samych warunków 73 D/NP została szczegółowo opisana w poprzednich pracach24. Wszystkie wykresy przedstawiają średnią ± SEM dla pomiarów powtórzonych dwukrotnie lub trzykrotnie.

Rycina 4. Reprezentatywne wyniki optymalizacji cykli PCR. Od lewej do prawej studzienki reprezentują: 4 cykle, 6 cykli, 8 cykli, 10 cykli, 12 cykli, kontrolę negatywną (brak matrycy DNA), standardową drabinkę DNA 25 bp. Optymalne warunki zaobserwowano dla 8 cykli, gdzie pojedynczy prążek produktu wykazuje wysoką intensywność, bez obecności produktów nadmiernej amplifikacji widocznych przy większej liczbie cykli.

Rycina 5. Optymalizacja czasu inkubacji w teście AuNP. Czas inkubacji w etapie AuNP/DNA przy 73 D/NP został zredukowany z O/N (Rycina 3) do 30 min, a po inkubacji z celem niezwłocznie dodano sól, zamiast czekać 20 min (Rycina 3). (A) Odpowiedź zaobserwowano dla kortyzolu (niebieskie romby), ale nie dla kwasu cholowego (zielone kółka) ani 2MNP (2-metoksynaftalenu; czerwone kwadraty). Wszystkie wykresy przedstawiają średnią ± SEM dla pomiarów wykonanych w dubletach lub tripletach. (B) Od lewej do prawej: próbka ślepa, 10 µM kortyzolu, 10 µM kwasu cholowego, 10 µM 2MNP. Zmiana wizualna dla kortyzolu jest widoczna gołym okiem. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.