$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
1. Przygotowanie mrożono-uwodnionej próbki EM
Ponieważ zespoły 9 białka kapsydu HIV-1 (CA) są stabilne tylko w buforze o wysokiej zawartości soli (1M NaCl), co przyczynia się do silnego szumu tła na obrazach krio-EM, używamy metody szybkiego rozcieńczania i blottingu z tyłu, aby przejściowo zmniejszyć stężenie soli podczas przygotowywania zamrożonej i uwodnionej siatki EM.
- Żar rozładowuje się po stronie węglowej miedzianych siatek R2/1 Quantifoil o oczkach 200 poniżej 25 mA przez 25 sekund.
- Użyj nebulizatora, aby doprowadzić wilgotność w komorze środowiskowej, domowym ręcznym tłoku grawitacyjnym, do 80%.
- W FEI Vitrobot Plunger dewar schłodzić ciekły etan za pomocą ciekłego azotu. Zamontuj dewara do zamrażania wgłębnego na ręcznym tłoku grawitacyjnym
.
- Nałóż 2,5 μl wstępnie zmontowanego roztworu CA na węglową stronę siatki, która jest zamontowana na kleszczach, i załaduj kleszczyki na tłok tak, aby węglowa strona siatki była skierowana od siebie.
- Dodać 3 μl buforu do rozcieńczania o niskiej zawartości soli (100 mM NaCl) na tylną stronę kratki i natychmiast osuszyć tylną stronę siatki kawałkiem bibuły filtracyjnej. Cała tylna powierzchnia siatki powinna być w bliskim kontakcie z bibułą filtracyjną przez około 6 sekund, przed wyjęciem bibuły filtracyjnej. Po wyjęciu bibuły filtracyjnej natychmiast zanurz siatkę w ciekłym etanowym.
- Zdejmij kleszczyki z tłoka i szybko przenieś siatkę do skrzynki do przechowywania siatki.
2. Mikroskopia krioelektronowa zespołów rurowych CA
- Załaduj zamrożoną, uwodnioną siatkę do mikroskopu elektronowego FEI Polara G2 pracującego przy napięciu 200 kv i wyposażonego w kamerę CCD Gatan 4Kx4K.
- W trybie wyszukiwania przy niskiej dawce, przy powiększeniu ~200x i dawce <0,001e-/TΥ2, przesiewaj całą siatkę w poszukiwaniu obszarów z odpowiednim lodem i zapisz pozycje tych obszarów w pliku etapu.
- Przywołaj zapisane pozycje i dalej przesiewaj te obszary w powiększeniu 3,900x w trybie wyszukiwania z niską dawką. Wybierz obszary z jednolitą, cienką warstwą lodu zawierającą dobrze oddzielone, długie rurki nad otworami, do zbierania danych. Zapisz lokalizacje tych obszarów w pliku drugiego etapu.
- Przełącz się w tryb ekspozycji przy powiększeniu 59 000x, wstaw aperturę obiektywu 100 μm i dostosuj stygmatyzm obiektywu i intensywność wiązki dla dawki ̃15 e-/Å 2 na ekspozycję.
- Wróć do trybu wyszukiwania niskiej dawki, przejdź do zapisanej pozycji i zidentyfikuj i wyśrodkuj dobrą rurkę za pomocą kamery CCD. Przełącz się w tryb ostrości, wyreguluj ostrość i ustaw wartość rozmycia, zwykle od 0.5 do 2.5 μm. Przełącz się w tryb ekspozycji, ustaw czas naświetlania 0,3-0,5 sekundy, dla dawki 15 e-/A 2 i zbierz obraz. Obrazy są zbierane na kamerze płytkowej, a filmy powinny być pozostawione na 10 sekund przed wykonaniem ekspozycji.
- Przejdź do następnej zapisanej pozycji i powtórz krok 5, aby zebrać więcej obrazów.
- Filmy są wywoływane w pełnej mocy D19 przez 12 minut i digitalizowane za pomocą skanera Nikon super coolscan 9000 ED o rozmiarze piksela 6,35 μm. Format plików graficznych to TIFF.
3. Indeksowanie spiralne
Obiekt spiralny może być indeksowany przez dwa parametry: porządek Bessela, n, i numer linii warstwy, l. Każda linia warstwy w transformacie Fouriera, scharakteryzowana przez (n, l), odpowiada zbiorowi linii na powierzchni sieci obiektu spiralnego oznaczonego wskaźnikami (h, k), przy użyciu notacji z kraty 2D. Dla dowolnej (h, k) linia warstwy (nhk, lhk) jest liniową kombinacją dwóch podstawowych wektorów (n10, l10) i ( n01, l01), które są wartościami n i l dwóch głównych linii warstwy (1, 0) i (0, 1). l można uzyskać z wysokości linii warstwy mierzonej wzdłuż osi Z w transformacie Fouriera. Wartość n można oszacować za pomocą następującego równania10
2πRr ≈ Jn ≈ 1.1|n|+0.9 ..................... (1)
gdzie Jn to funkcja Bessela, która określa intensywność n-tej linii warstwy, r to promień obiektu spiralnego, a R to promień maksymalnej amplitudy linii warstwy. Numer linii warstwy, l, jest powiązany z n przez regułę wyboru11
l = tn + um .................... (2)
gdzie t i u są stałymi helisy. Dla każdej helisy mogą występować dokładnie jednostki u w dokładnie (lub bardzo blisko) t pełnych obrotów.
- Zapakuj stosunkowo prostą i długą rurkę o jednolitej średnicy za pomocą programu ELAN12 helixboxer i zapisz obraz w formacie MRC.
- Określ odległość powtórzeń śrubowych, c, za pomocą programu opartego na korelacji krzyżowej, takiego jak "imgccf", w pakiecie MRC13.
- Oblicz transformatę Fouriera z nową długością pola, która jest całką z odległości powtórzenia, c.
- Wybierz dwie główne linie warstwy, które definiują dwa podstawowe wektory kraty powierzchni (1,0) i (0,1).
- Zmierz promień rury, r, oraz wysokość i promień dwóch linii warstwy głównej w transformacie Fouriera, odpowiednio l10, R10, l01 R01 (rys. 1).
- Obliczyć n10 i n01 zgodnie z równaniem (1).
- Ponieważ n wartości są tylko szacunkowe, kilka kombinacji n10 i n01 jest testowanych w późniejszych krokach (patrz krok 4.2) w celu znalezienia prawidłowej symetrii śrubowej za pomocą pakietu programu IHRSR.
- Oblicz symetrię opisaną dwiema liczbami rzeczywistymi: obrotem między podjednostkami (Δφ) i wzrostem osiowym (Δz) helisy jednogwiazdowej (n=1). Biorąc pod uwagę wartości l10, l01, n10 i n01, u i t można uzyskać, testując zakres wartości m (na przykład -50<m<50) zgodnie z regułą wyboru. Na koniec Δφ i Δz oblicza się przy użyciu Δφ =360t/u i Δz = c/u.
4. Rekonstrukcja trójwymiarowa
- Segmentacja
cząstek
- Otwórz mikrofotografię zawierającą spiralne cząstki za pomocą programu graficznego Boxer, który jest programem w pakiecie EMAN.
- Pokrój spiralną cząstkę na nakładające się na siebie segmenty. W panelu sterowania Boxera wybierz tryb Helix i ustaw parametry dla pudełka: rozmiar pudełka powinien być większy niż średnica cząstki, a wartość "Olap" powinna wynosić ̃90% rozmiaru pudełka.
- Po kliknięciu lewym przyciskiem myszy na którymś z końców spiralnej cząstki, Boxer automatycznie wygeneruje serię pól cząstek wzdłuż długości helisy.
- Zapisz segmenty w ramkach wraz z ich współrzędnymi.
- Wstępna rekonstrukcja 3D za pomocą programów
IHRSR
- Odwróć kontrast obrazów krio-EM i zastosuj filtrowanie dolnoprzepustowe14 (opcjonalnie) przed przetworzeniem za pomocą iteracyjnej metody rekonstrukcji helisalnej przestrzeni rzeczywistej (IHRSR).
- Otwórz interfejs graficzny programu IHRSR, wpisując "Generator". Zapewnij interfejs graficzny ze wszystkimi informacjami dotyczącymi stosu cząstek w pudełkach, w tym nazwą i ścieżką stosu, liczbą obrazów w stosie, wartościami parametrów symetrii itp. Kliknij przycisk "Zakończ", aby utworzyć skrypt rekonstrukcyjny b25.spi.
- Użyj pełnego lub pustego cylindra jako początkowego odniesienia i pozwól, aby procedura powtarzała się, aż nie nastąpią żadne zmiany w zdefiniowanej symetrii, co zwykle następuje po kilku cyklach. Symetria prawej helisy powinna dawać stabilnie zbieżną rekonstrukcję. Rekonstrukcja wygenerowana w ostatnim cyklu zostanie wykorzystana jako wstępny punkt odniesienia do dalszego udoskonalenia. IHRSR wykonuje rekonstrukcję 3D za pomocą programów SPIDER.
- Rekonstrukcja z iteracyjnym udoskonalaniem15
Rekonstrukcja 3D wygenerowana przez IHRSR jest teraz używana jako początkowe odniesienie do dodatkowych udoskonaleń. Podczas zagęszczenia symetria śrubowa jest ustalana na poziomie Δφ i Δz, które są określane na podstawie procedury IHRSR.
- Określ rozogniskowanie i astygmatyzm obecne na mikrofotografii za pomocą programów CTFFIND3 i CTFTILT16.
- Pomnóż segmenty cząstek przez CTF za pomocą programów SPIDER17. Operacja określana jako FT w pakiecie SPIDER służy do obliczania transformaty Fouriera (FFT) obrazu 2D. Następnie FFT jest mnożony przez wartości CTF, określone w poprzednich krokach, przy użyciu operacji zwanej MU w pakiecie SPIDER. Uzyskana wartość jest następnie przekształcana odwrotnie w celu utworzenia nowego obrazu (skorygowanego przez CTF) w przestrzeni rzeczywistej za pomocą operacji SPIDER, ponownie FT.
- Wykonaj dopasowywanie projekcji, porównując projekcje objętości odniesienia z obrazami skorygowanymi przez CTF, używając wyrównania do wielu odniesień. Zmiana kąta nachylenia poza płaszczyzną jest ograniczona do +/-10° i próbkowana w krokach co 1°. Wprowadź ograniczenia, takie jak wysokie współczynniki korelacji, kąty w płaszczyźnie bliskie 0° lub 180° oraz ograniczone przesunięcia x dla parametrów wyrównania każdego segmentu. W rekonstrukcji należy uwzględnić tylko te segmenty, które spełniają ograniczenia.
- Po każdym iteracyjnym cyklu udoskonalania generowana jest rekonstrukcja 3D przy użyciu projekcji wstecznej i dzielona przez sumę przez CTFˆ2.
- Zastosuj symetrię spiralną, aby wygenerować symetryczną objętość. Iteracyjne udoskonalanie jest przerywane, gdy nie nastąpi dalsza poprawa rozdzielczości nowej rekonstrukcji 3D.
- Pakiet do przetwarzania obrazu SPIDER jest używany do większości etapów udoskonalania. Szereg operacji jest kontrolowany przez skrypty SPIDER, które są tworzonymi przez użytkownika wsadowymi plikami kontrolnymi zawierającymi sekwencje operacji i wartości parametrów. Ostateczna rekonstrukcja jest obliczana za pomocą programów w pakiecie SPIDER.
5. Reprezentatywne wyniki:
Pojedyncza rurka HIV-1 CA A92E (Rys. 1a) została zapakowana, a jej transformacja Fouriera (Rys. 1b) została obliczona dla indeksowania spiralnego. Dla linii warstw (1, 0) i (0, 1), l10=28, l01=37, R10=55, R01=44. Biorąc pod uwagę promień rury wynoszący 211,57A, przybliżyliśmy n10=-12, n01=11 (tutaj ręczność była z góry określona). Przy powtarzalnej odległości 5195,48A symetria rury została określona jako Δz=6,8093A, Δφ=328,88°. Δz i Δφ zostały oczyszczone do 7,1321A i 328,86° za pomocą IHRSR (rys. 2a), a wstępną rekonstrukcję przedstawiono na rys. 2b. Ostateczna rekonstrukcja (rys. 3), po iteracyjnym udoskonaleniu, znacznie poprawiła mapę gęstości w stosunku do początkowego modelu obliczonego za pomocą IHRSR (rys. 2b).

Rysunek 1. Indeksowanie rurki spiralnej HIV-1 CA. (A) Pojedynczy obraz probówki HIV-1 CA A92E. Podziałka liniowa, 30 nm. (B) Transformata Fouriera rurki pokazanej w (A). Indeksy spiralne (n10=-12, n01=11) są oznaczone. Grot strzałki wskazuje linię warstwy w rozdzielczości 23A.

Rysunek 2. Wstępna rekonstrukcja z wykorzystaniem IHRSR. (A) Wyznaczanie symetrii dla każdego cyklu iteracyjnego. Δφ i Δz, zaczynając od wartości początkowych, zbiegają się do wartości stabilnych po 10 iteracyjnych cyklach uszlachetniania. (B) Wstępna mapa gęstości 3D po 10 cyklach iteracyjnych.

Rysunek 3. Mapa gęstości 3D po iteracyjnym doprecyzowaniu. (A-C) Mapa gęstości rur CA jest wyświetlana jako trzy ortogonalne wycinki: równoległe do osi rury i blisko powierzchni (A), prostopadłe do osi rury (B) oraz równoległe do i przez oś rury (C). Podziałka liniowa, 10 nm. (D) Renderowanie powierzchni mapy gęstości 3D wyprofilowanej w czasie 1,8 s, obejmującej 100% objętości.