$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Efficient library construction is critical to ensure sufficient coverage of the mutational space. To avoid loss of rare functional variants, the number of transformants should reach at least 10× to 100× the theoretical library size. Because transformation efficiency in S. cerevisiae is influenced by cell concentration and is highest in early- to mid-log-phase cells16, we optimized the initial cell density prior to electroporation. Yeast cultures were subcultured at initial OD600 values ranging from 1.0–4.0 and incubated for 4 h (n = 2–4 independent biological experiments), after which the final OD600 values were mapped onto the S. cerevisiae growth curve to confirm that cells remained within the mid-log growth phase (Figure 5A, Supplementary Figure 1).
Electroporation of these cultures with 4 µg plasmid DNA (n = 2 independent biological experiments) showed that transformation efficiency differed significantly across initial OD600 conditions (one-way ANOVA, p = 8.92 × 10−7). Relative to cultures initiated at OD600 = 1.0, all higher starting cell densities produced significantly greater transformation efficiencies (Dunnett's multiple-comparisons test; OD600 = 1.5, adjusted p = 0.0020; OD600 ≥ 2.0, adjusted p < 0.0001), with the highest transformation efficiency observed at an initial OD600 of 4.0. To take advantage of yeast cells' natural ability to perform homologous recombination, we also electroporated 4 µg of the linearized plasmid vector and 12 µg of the insert (i.e., protease library) into yeast cells (n = 2 independent biological experiments). Transformation efficiency differed significantly among the tested initial OD600 conditions (one-way ANOVA, p = 0.0023). Relative to cultures initiated at OD600 = 1.0, transformation efficiencies were significantly higher at OD600 = 3.0 (adjusted p = 0.0018), 3.5 (adjusted p = 0.0062), and 4.0 (adjusted p = 0.0024) (Dunnett's multiple-comparisons test). We observed the highest transformation efficiency at a moderate initial subculture density (OD600 = 3.0), indicating an optimal balance between cell competency and recombination efficiency under these conditions (Figure 5B). Compared with plasmid-only transformation, co-electroporation of a linearized vector and protease library insert (n = 2 independent biological replicates) yielded significantly higher transformation efficiencies across the tested initial OD600 conditions (two-way ANOVA on log10-transformed transformation efficiencies, p = 9.40 × 10−12). On average, co-electroporation produced a mean 121-fold increase in transformation efficiency (95% CI, 63–233-fold) relative to plasmid-only transformation (Figure 5B).
During a FACS run of yeast cells containing the mutant protease library and the substrate selection cassette (Figure 6A), cells will either display a truncated cassette or a protease that still cleaves the counter-selection substrate (Q1), express non-functional protease variants (Q2), or express protease variants that successfully cleave the selection substrate only (Q3). Yeast cells expressing non-functional variants appeared as a dominant population distributed along the diagonal in PE versus FITC plots, reflecting co-retention of both signals (Figure 6A, Q2). By contrast, functional variants that cleaved the selection substrate while sparing the native (counter-selection) substrate formed a distinct off-diagonal population characterized by high-PE and medium-to-low FITC signal intensity (Figure 6A, Q3). Applying a stringent gating strategy to isolate this off-diagonal population resulted in progressive enrichment of functional variants across successive rounds of sorting10. This enrichment is visually observed as an increasingly defined cluster of cells separated from the diagonal cluster (Figure 6B, left to right). If no distinct off-diagonal population is observed, this may indicate low library quality, insufficient expression, or suboptimal staining conditions. As a negative control, Figure 6C,D shows electroporated, induced, unstained yeast cells and the corresponding FITC and PE fluorescence histograms.
Within the enriched off-diagonal population, subpopulations can be further resolved based on FITC intensity. Cells exhibiting a medium FITC signal can be sorted separately from those with a low FITC signal. Sequencing these distinct populations revealed unique mutation profiles (Table 11), indicating that varying cleavage efficiencies are associated with different protease variants.

Figure 1: Vector and Insert Preparation. (A) Generation of potyviral protease libraries via site-saturation mutagenesis. (B) Substrate cassette assembly and integration into yeast display vector (pDD1523). (C) Linearization of yeast display vector (pDD1524) via restriction digestion. Please click here to view a larger version of this figure.

Figure 2: Electroporation of the substrate cassette and the PCR-amplified mutant protease library into competent EBY100 cells. (Dashed box) In vivo circularization of the plasmid containing both protease and substrate cassette. Please click here to view a larger version of this figure.

Figure 3: Induction, staining, and flow cytometric sorting of yeast cells. (A) Induction of protease library and substrate cassette. (B) Fluorescent antibody staining: Staining with anti-HA fluorescein isothiocyanate (FITC) and anti-FLAG phycoerythrin (PE) antibodies. (C) Stained cells are sorted using Fluorescence-Activated Cell Sorting (FACS). (D) Recovery of FACS-sorted cells in selection media. Please click here to view a larger version of this figure.

Figure 4: Plasmid recovery and sequencing of sorted yeast cells. (A) Plating of sorted yeast cells. (B) Plasmid extraction from individual yeast colonies. (C) E. coli transformation of extracted plasmids. (D) Plasmid extraction from individual E. coli colonies. (E) Nanopore sequencing of extracted plasmids. Flask and tube icons adapted from DBCLS and Helicase 11, respectively, via Bioicons (https://bioicons.com) and adapted under the terms of the Creative Commons Attribution 4.0 International License (CC BY 4.0). Please click here to view a larger version of this figure.

Figure 5: Yeast growth curves in different growth media and the influence of initial Optical Density at 600 nm (OD600) of subculture on yeast electroporation efficiency. (A) Growth curves of yeast in 1 × Yeast Extract Peptone Dextrose (YPD), 2 × YPD, Synthetic Dextrose Casamino Acids Agar without antibiotics (SDCAA-), and Synthetic Dextrose Casamino Acids Agar with antibiotics (SDCAA+). (B) Transformation efficiency of yeast electroporated with 4 µg plasmid vector at various subculture initial OD600, and the transformation efficiency using linearized vector (4 µg) and protease library insert (12 µg) at the same initial OD600 values. Data in (A) are presented as OD600 ± SD from two independent biological experiments (n = 2). Data in (B) are presented as mean number of transformants per electroporation on a logarithmic y-axis, with individual biological replicates shown as circles (n = 2). Statistical significance for the plasmid and vector-plus-insert electroporation optimization experiments was determined using one-way ANOVA on log10-transformed electroporation efficiencies followed by Dunnett's multiple-comparisons test using OD600 = 1.0 as the control. ns, not significant; p < 0.01 (**); p < 0.0001 (****). Statistical significance of the difference between plasmid-only and vector-plus-insert transformation efficiencies was determined using two-way ANOVA on log10-transformed transformation efficiencies. p < 0.0001 (****). Please click here to view a larger version of this figure.

Figure 6: Fluorescence-Activated Cell Sorting (FACS)-based enrichment of yeast clones. (A) Schematic representation of the gating strategy. Cells were divided into three regions: Q1, phycoerythrin (PE)-low cells; Q2, PE-high/fluorescein isothiocyanate (FITC)-high cells; and Q3, PE-high cells with low-to-intermediate FITC fluorescence. Cells within Q3 were collected for enrichment. (B) Representative FACS plots of one independent evolutionary campaign showing sequential rounds of sorting with increasing FITC gating stringency. Mid- and low-FITC cells were collected separately for sequencing. Percentages indicate the fraction of live, single yeast cells within the gated population. (C) Representative FACS plot showing electroporated, unstained yeast cells as a negative control. (D) Histograms of FITC-A (left) and PE-A (right) fluorescence intensities corresponding to the data shown in (C). Please click here to view a larger version of this figure.
Table 1: Reaction mixture for 1st-round PCR amplification of the protease library. Composition of the PCR reaction used for the first-round amplification of the protease library, including the final concentrations and amounts of each reagent in a 50 µL reaction. Please click here to download this Table.
Table 2: Cycling conditions for 1st-round PCR amplification of the protease library. The thermal cycling program was used for the first-round PCR amplification of the protease library. Please click here to download this Table.
Table 3: Reaction mixture for 2nd-round PCR amplification of the protease library. Composition of the PCR reaction used for the second-round amplification of the protease library, including the final concentrations and amounts of each reagent in a 50 µL reaction. Please click here to download this Table.
Table 4: Cycling conditions for 2nd-round PCR amplification of the protease library. The thermal cycling program was used for the second-round PCR amplification of the protease library. Please click here to download this Table.
Table 5: Media and Reagents for Yeast Electroporation and Growth. Tabular description of the mass required from each chemical component to formulate specific amounts of 1 × and 2 × Yeast Extract Peptone Dextrose (YPD), electroporation buffer, conditioning buffer, outgrowth media, Synthetic Dextrose Casamino Acids with antibiotics (SDCAA+), Synthetic Galactose Casamino Acids with antibiotics (SGCAA+), Lysogeny Broth (LB) media, LB-chloramphenicol media, Super Optimal Broth (SOB) media, and Super Optimal Broth with Catabolite repression (SOC) media. Instructions for preparing stock solutions, such as 2 M Sorbitol, 1 M dithiothreitol (DTT), and 2 M lithium acetate (LiAc), are included. Additional instructions for preparing SDCAA, LB, and YPD agar plates are provided. Please click here to download this Table.
Table 6: Reaction mixture for annealing and phosphorylation of oligos. Composition of the reaction mixture used for simultaneous annealing and phosphorylation of oligonucleotides for substrate cassette assembly. Please click here to download this Table.
Table 7: Cycling conditions for annealing and phosphorylation of oligos. Temperature program used for oligonucleotide annealing and phosphorylation for substrate cassette assembly. Please click here to download this Table.
Table 8: Reaction mixture for Golden Gate assembly of substrate cassette and vector plasmid. Composition of the Golden Gate assembly reaction used to ligate the substrate cassette into the vector plasmid. Please click here to download this Table.
Table 9: Cycling conditions for Golden Gate assembly of substrate cassette and vector plasmid. Thermal cycling program used for Golden Gate assembly. Please click here to download this Table.
Table 10: Reaction mixture for linearization of vector plasmid. Composition of the restriction digestion reaction used to linearize the vector plasmid prior to library construction. Please click here to download this Table.
Table 11: Nanopore sequencing results of TEVp variants cleaving ENLYFES. Plasmids were extracted from yeast cells after the fourth and final FACS session. Mutations in each clone are shown relative to the wild-type sequence. Please click here to download this Table.
Supplementary Figure 1: Relationship between initial and final EBY100 culture density during subculture. EBY100 cultures were subcultured at varying initial OD600 values, and the final OD600 was measured after incubation under identical growth conditions (30 °C, 220 rpm). Final culture density increased linearly with increasing initial OD600 (R2 = 0.969, p = 5.77 × 10−5). Data points represent means ± SD of (n = 2–4) independent biological replicates.Please click here to download this file.
Supplementary Figure 2: Rational design of the TEVp mutagenesis library. Cartoon representation of TEVp (PDB ID: 1LVB) with the four residues (gold) <5Å away from the P1 residue (purple) targeted for site-saturation mutagenesis.Please click here to download this file.