$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Efektywne tworzenie biblioteki jest kluczowe dla zapewnienia wystarczającego pokrycia przestrzeni mutacyjnej. Aby uniknąć utraty rzadkich wariantów funkcjonalnych, liczba transformantów powinna osiągać od 10× do 100× teoretycznego rozmiaru biblioteki. Ponieważ wydajność transformacji w S. cerevisiae zależy od stężenia komórek i jest najwyższa w komórkach we wczesnej i środkowej fazie logarytmicznej16, zoptymalizowano początkową gęstość komórek przed elektroporacją. Kultury drożdży przesiewano przy początkowych wartościach OD600 w zakresie 1,0–4,0 i inkubowano przez 4 h (n = 2–4 niezależne eksperymenty biologiczne), po czym końcowe wartości OD600 naniesiono na krzywą wzrostu S. cerevisiae, aby potwierdzić, że komórki pozostały w środkowej fazie wzrostu logarytmicznego (Rysunek 5A, Rysunek uzupełniający 1).
Elektroporacja tych kultur z użyciem 4 µg plazmidowego DNA (n = 2 niezależne eksperymenty biologiczne) wykazała, że wydajność transformacji różniła się znacząco w zależności od początkowych warunków OD600 (jednoczynnikowa ANOVA, p = 8,92 × 10−7). W porównaniu z kulturami rozpoczynanymi przy OD600 = 1,0, wszystkie wyższe początkowe gęstości komórek wykazały znacząco większą wydajność transformacji (test wielokrotnych porównań Dunnetta; OD600 = 1,5, skorygowane p = 0,0020; OD600 ≥ 2,0, skorygowane p < 0,0001), przy czym najwyższą wydajność transformacji zaobserwowano przy początkowej wartości OD600 wynoszącej 4,0. Aby wykorzystać naturalną zdolność komórek drożdży do rekombinacji homologicznej, przeprowadziliśmy również elektroporację komórek drożdży z użyciem 4 µg zlinearyzowanego wektora plazmidowego i 12 µg wkładu (tj. biblioteki proteaz) (n = 2 niezależne eksperymenty biologiczne). Wydajność transformacji różniła się znacząco pomiędzy testowanymi początkowymi warunkami OD600 (jednoczynnikowa ANOVA, p = 0,0023). W porównaniu z kulturami rozpoczynanymi przy OD600 = 1,0, wydajności transformacji były znacząco wyższe przy OD600 = 3,0 (skorygowane p = 0,0018), 3,5 (skorygowane p = 0,0062) oraz 4,0 (skorygowane p = 0,0024) (test wielokrotnych porównań Dunnetta). Najwyższą wydajność transformacji zaobserwowano przy umiarkowanej początkowej gęstości subkultury (OD600 = 3,0), co wskazuje na optymalną równowagę między kompetencją komórek a wydajnością rekombinacji w tych warunkach (Rycina 5B). W porównaniu z transformacją samym plazmidem, ko-elektroporacja zlinearyzowanego wektora i wkładu biblioteki proteaz (n = 2 niezależne powtórzenia biologiczne) dała znacząco wyższą wydajność transformacji we wszystkich testowanych początkowych warunkach OD600 (dwuczynnikowa ANOVA dla wydajności transformacji po transformacji log10, p = 9,40 × 10−12). Średnio ko-elektroporacja spowodowała średnio 121-krotny wzrost wydajności transformacji (95% CI, 63–233-krotny) w stosunku do transformacji samym plazmidem (Rycina 5B).
Podczas analizy FACS komórek drożdży zawierających bibliotekę zmutowanej proteazy oraz kasetę selekcji substratu (Ryc. 6A), komórki będą wykazywać albo skróconą kasetę lub proteazę, która nadal rozcina substrat kontroselekcyjny (Q1), wyrażać niefunkcjonalne warianty proteazy (Q2), albo wyrażać warianty proteazy, które skutecznie rozcinają wyłącznie substrat selekcyjny (Q3). Komórki drożdży wyrażające niefunkcjonalne warianty pojawiły się jako dominująca populacja rozciągająca się wzdłuż przekątnej na wykresach PE względem FITC, co odzwierciedla współwystępowanie obu sygnałów (Ryc. 6A, Q2). W przeciwieństwie do nich, funkcjonalne warianty, które rozcięły substrat selekcyjny, oszczędzając natywny substrat (kontroselekcyjny), utworzyły odrębną populację poza przekątną, charakteryzującą się wysoką intensywnością sygnału PE oraz średnią do niskiej intensywności sygnału FITC (Ryc. 6A, Q3). Zastosowanie rygorystycznej strategii bramkowania w celu wyizolowania tej populacji poza przekątną doprowadziło do progresywnego wzbogacania wariantów funkcjonalnych w kolejnych rundach sortowania10. Wzbogacenie to jest widoczne jako coraz wyraźniejsza grupa komórek oddzielona od grupy na przekątnej (Ryc. 6B, od lewej do prawej). Brak wyraźnej populacji poza przekątną może wskazywać na niską jakość biblioteki, niewystarczającą ekspresję lub nieoptymalne warunki barwienia. Jako kontrola negatywna, Ryc. 6C,D przedstawia elektroporowane, indukowane, niebarwione komórki drożdży oraz odpowiadające im histogramy fluorescencji FITC i PE.
W obrębie wzbogaconej populacji poza przekątną można dalej wyodrębnić subpopulacje na podstawie intensywności FITC. Komórki wykazujące średni sygnał FITC mogą być sortowane oddzielnie od tych o niskim sygnale FITC. Sekwencjonowanie tych odrębnych populacji ujawniło unikalne profile mutacji (Tabela 11), co wskazuje, że zróżnicowana wydajność cięcia wiąże się z różnymi wariantami proteazy.

Rysunek 1: Przygotowanie wektora i wkładki. (A) Generowanie bibliotek proteazy potyviralnej poprzez mutagenezę nasycenia miejsc. (B) Montaż kasety substratu i integracja z wektorem do ekspozycji drożdżowej (pDD1523). (C) Linearyzacja wektora do ekspozycji drożdżowej (pDD1524) poprzez trawienie restrykcyjne. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Elektroporacja kasety substratu oraz biblioteki zmutowanej proteazy amplifikowanej metodą PCR do kompetentnych komórek EBY100. (Przerywana ramka) Cyrkularyzacja in vivo plazmidu zawierającego zarówno proteazę, jak i kasetę substratu. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Indukcja, barwienie i sortowanie cytometryczne komórek drożdży. (A) Indukcja biblioteki proteaz i kasety substratu. (B) Barwienie przeciwciałami fluorescencyjnymi: Barwienie przeciwciałami anty-HA znakowanymi izotiocyjanianem fluoresceiny (FITC) oraz anty-FLAG znakowanymi fikoerytryną (PE). (C) Barwione komórki są sortowane przy użyciu fluorescencyjnego sortowania komórek (FACS). (D) Odzyskiwanie komórek posortowanych metodą FACS w pośrodkach selekcyjnych. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Odzyskiwanie plazmidów i sekwencjonowanie posortowanych komórek drożdży. (A) Wysiew posortowanych komórek drożdży. (B) Ekstrakcja plazmidów z poszczególnych kolonii drożdży. (C) Transformacja E. coli wyekstrahowanymi plazmidami. (D) Ekstrakcja plazmidów z poszczególnych kolonii E. coli. (E) Sekwencjonowanie Nanopore wyekstrahowanych plazmidów. Ikony kolb i probówek zaadaptowane odpowiednio z DBCLS i Helicase 11 za pośrednictwem Bioicons (https://bioicons.com) i dostosowane zgodnie z warunkami licencji Creative Commons Attribution 4.0 International License (CC BY 4.0). Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 5: Krzywe wzrostu drożdży w różnych mediach wzrostowych i wpływ początkowej gęstości optycznej przy 600 nm (OD600) subkultury na wydajność elektroporacji drożdży. (A) Krzywe wzrostu drożdży w 1 × Yeast Extract Peptone Dextrose (YPD), 2 × YPD, syntetycznym agarze z dekstrozą i kwasami kazeinowymi bez antybiotyków (SDCAA-) oraz syntetycznym agarze z dekstrozą i kwasami kazeinowymi z antybiotykami (SDCAA+). (B) Wydajność transformacji drożdży elektroporowanych 4 µg wektora plazmidowego przy różnych początkowych wartościach OD600 subkultury oraz wydajność transformacji z użyciem zlinearyzowanego wektora (4 µg) i wstawki z biblioteki proteaz (12 µg) przy tych samych wartościach początkowych OD600. Dane w (A) przedstawiono jako OD600 ± SD z dwóch niezależnych eksperymentów biologicznych (n = 2). Dane w (B) przedstawiono jako średnią liczbę transformantów na jedną elektroporację na logarytmicznej osi y, gdzie poszczególne powtórzenia biologiczne zaznaczono kółkami (n = 2). Istotność statystyczną dla eksperymentów optymalizacji elektroporacji plazmidu oraz wektora z wstawką określono za pomocą jednoczynnikowej analizy wariancji (ANOVA) dla wartości wydajności elektroporacji przekształconych logarytmicznie (log10), a następnie za pomocą testu wielokrotnych porównań Dunnetta, przyjmując OD600 = 1,0 jako grupę kontrolną. ns, brak istotności; p < 0,01 (**); p < 0,0001 (****). Istotność statystyczną różnicy między wydajnością transformacji samego plazmidu a wektora z wstawką określono za pomocą dwuczynnikowej analizy wariancji (ANOVA) dla wartości wydajności transformacji przekształconych logarytmicznie (log10). p < 0,0001 (****). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 6: Wzbogacanie klonów drożdży metodą sortowania komórek aktywowane fluorescencją (FACS). (A) Schematyczna reprezentacja strategii bramkowania. Komórki podzielono na trzy regiony: Q1, komórki z niską fluorescencją fikocyjaniny (PE); Q2, komórki z wysoką fluorescencją PE i wysoką fluorescencją izotiocyjanianu fluoresceiny (FITC); oraz Q3, komórki z wysoką fluorescencją PE i niską do umiarkowanej fluorescencją FITC. Komórki z regionu Q3 zostały zebrane w celu wzbogacenia. (B) Reprezentatywne wykresy FACS z jednej niezależnej kampanii ewolucyjnej, pokazujące sekwencyjne rundy sortowania z coraz większą rygorystycznością bramkowania FITC. Komórki z umiarkowaną i niską fluorescencją FITC zbierano oddzielnie do sekwencjonowania. Procenty wskazują frakcję żywych, pojedynczych komórek drożdży w obrębie bramkowanej populacji. (C) Reprezentatywny wykres FACS pokazujący elektroporowane, niebarwione komórki drożdży jako kontrolę negatywną. (D) Histogramy intensywności fluorescencji FITC-A (lewo) i PE-A (prawo) odpowiadające danym przedstawionym na (C). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Tabela 1: Mieszanina reakcyjna dla pierwszej rundy amplifikacji PCR biblioteki proteaz. Skład reakcji PCR wykorzystywanej do pierwszej rundy amplifikacji biblioteki proteaz, zawierający końcowe stężenia i ilości każdego odczynnika w reakcji o objętości 50 µL. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Warunki cykliczne dla pierwszej rundy amplifikacji PCR biblioteki proteaz.Program cykli termicznych został zastosowany podczas pierwszej rundy amplifikacji PCR biblioteki proteaz. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 3: Mieszanina reakcyjna dla drugiej rundy amplifikacji PCR biblioteki proteaz. Skład reakcji PCR wykorzystywanej do drugiej rundy amplifikacji biblioteki proteaz, zawierający stężenia końcowe i ilości każdego odczynnika w reakcji o objętości 50 µL. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 4: Warunki cykliczności dla drugiej rundy amplifikacji PCR biblioteki proteaz. Program cykli termicznych zastosowano do drugiej rundy amplifikacji PCR biblioteki proteaz. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 5: Pożywki i odczynniki do elektroporacji i hodowli drożdży. Opis tabelaryczny masy wymaganej z każdego składnika chemicznego do przygotowania określonych ilości 1 × i 2 × ekstraktu drożdżowego, peptonu i dekstrozy (YPD), buforu do elektroporacji, buforu kondycjonującego, pożywki do odzysku, syntetycznej dekstrozy z kwasami kazeinowymi i antybiotykami (SDCAA+), syntetycznej galaktozy z kwasami kazeinowymi i antybiotykami (SGCAA+), pożywki Lysogeny Broth (LB), pożywki LB z chloramfenikolem, pożywki Super Optimal Broth (SOB) oraz pożywki Super Optimal Broth z represją kataboliczną (SOC). Dołączono instrukcje przygotowania roztworów zapasowych, takich jak 2 M sorbitol, 1 M ditiotreitol (DTT) i 2 M octan litu (LiAc). Podano dodatkowe instrukcje przygotowania płytek agarowych SDCAA, LB i YPD. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 6: Mieszanina reakcyjna do hybrydyzacji i fosforylacji oligonukleotydów. Skład mieszaniny reakcyjnej stosowanej do jednoczesnej hybrydyzacji i fosforylacji oligonukleotydów w celu montażu kasety substratu. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 7: Warunki cykliczne dla hybrydyzacji i fosforylacji oligonukleotydów. Program temperaturowy stosowany do hybrydyzacji i fosforylacji oligonukleotydów w celu montażu kasety substratowej. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 8: Mieszanina reakcyjna dla metody Golden Gate assembly w celu połączenia kasety substratu i plazmidu wektora. Skład reakcji Golden Gate assembly wykorzystywanej do ligacji kasety substratu do plazmidu wektora. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 9: Warunki cykliczne dla montażu Golden Gate kasety substratu i plazmidu wektora. Program cykli termicznych zastosowany do montażu Golden Gate. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 10: Mieszanina reakcyjna do linearyzacji plazmidu wektora. Skład reakcji trawienia restrykcyjnego stosowanego do linearyzacji plazmidu wektora przed konstrukcją biblioteki. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 11: Wyniki sekwencjonowania nanoporowego wariantów TEVp rozszczepiających ENLYFES. Plazmidy wyekstrahowano z komórek drożdży po czwartej i ostatniej sesji FACS. Mutacje w każdym klonie przedstawiono w odniesieniu do sekwencji typu dzikiego. Kliknij tutaj, aby pobrać tę tabelę.
Rysunek uzupełniający 1: Zależność między początkową a końcową gęstością hodowli EBY100 podczas pasażowania. Hodowle EBY100 pasażowano przy różnych początkowych wartościach OD600, a końcową wartość OD600 mierzono po inkubacji w identycznych warunkach wzrostu (30 °C, 220 rpm). Końcowa gęstość hodowli wzrastała liniowo wraz ze wzrostem początkowej wartości OD600 (R2 = 0.969, p = 5.77 × 10−5). Punkty danych reprezentują średnie ± SD z (n = 2–4) niezależnych powtórzeń biologicznych.Prosimy kliknąć tutaj, aby pobrać ten plik.
Rysunek uzupełniający 2: Racjonalne projektowanie biblioteki mutagenezy TEVp. Schematyczna reprezentacja TEVp (PDB ID: 1LVB) z czterema resztami (złote), znajdującymi się w odległości <5Å od reszty P1 (fioletowej), wybranymi do nasyconej mutagenezy miejsce.Kliknij tutaj, aby pobrać ten plik.