High-throughput Analysis of in vivo Protein Stability*

Determining the half-life of proteins is critical for an understanding of virtually all cellular processes. Current methods for measuring in vivo protein stability, including large-scale approaches, are limited in their throughput or in their ability to discriminate among small differences in stability. We developed a new method, Stable-seq, which uses a simple genetic selection combined with high-throughput DNA sequencing to assess the in vivo stability of a large number of variants of a protein. The variants are fused to a metabolic enzyme, which here is the yeast Leu2 protein. Plasmids encoding these Leu2 fusion proteins are transformed into yeast, with the resultant fusion proteins accumulating to different levels based on their stability and leading to different doubling times when the yeast are grown in the absence of leucine. Sequencing of an input population of variants of a protein and the population of variants after leucine selection allows the stability of tens of thousands of variants to be scored in parallel. By applying the Stable-seq method to variants of the protein degradation signal Deg1 from the yeast Matα2 protein, we generated a high-resolution map that reveals the effect of ∼30,000 mutations on protein stability. We identified mutations that likely affect stability by changing the activity of the degron, by leading to translation from new start codons, or by affecting N-terminal processing. Stable-seq should be applicable to other organisms via the use of suitable reporter proteins, as well as to the analysis of complex mixtures of fusion proteins.

The regulation of protein stability is critical in order for cells to maintain proper functioning of almost every process. Thus, approaches for measuring the in vivo stability of a protein are essential for the identification of components in proteolytic pathways that affect protein turnover and for an understanding of the consequences of their activities. These approaches include traditional biochemical methods such as the Western blot, in which samples taken from time points after the inhibition of protein expression are fractionated via gel electrophoresis and the relevant protein is visualized with the use of an antibody. Another method allows one to track the degradation rate of newly synthesized proteins by metabolically labeling proteins with a radioisotope and following their radioactivity. A third method fuses a protein to a reporter enzyme like ␤-galactosidase, allowing the steady-state level of a protein to be measured via the enzymatic activity of the reporter enzyme. However, these small-scale methods are limited in the number of samples that they can analyze.
Large-scale methods have been developed that allow the simultaneous quantitation of the in vivo stability of many proteins. For example, Yen et al. fused ϳ8,000 human proteins to green fluorescent protein (GFP) 1 and followed the amount of each protein over time by using fluorescence-activated cell sorting (FACS) (1). To identify proteins, the plasmids encoding the GFP fusions were isolated and PCR products derived from these plasmids were hybridized to a DNA microarray. This method was applied to identify the substrates of a ubiquitin ligase complex (2). However, this method is limited by the number of bins into which protein fusions can be sorted in the FACS analysis and, consequently, how fine changes in stability can be discriminated. Alternatively, quantitative mass spectrometry has been used to analyze the stability of native proteins (3), but this approach often requires costly labeling and extensive data analysis. Moreover, these large-scale methods generally cannot distinguish differences in in vivo stability, which are sometimes significant, that result from small changes in a protein, such as single amino acid substitutions.
We present a method, Stable-seq, for measuring the in vivo stability of large numbers of variants of a protein that combines a simple genetic selection with high-throughput DNA sequencing. Stable-seq is a form of deep mutational scanning (4,5), in which a physical association between each protein variant and the DNA that encodes it allows DNA sequencing to score the frequency of each plasmid in a population. Here, our strategy is to fuse protein variants to a stable biosynthetic enzyme, the stability of which becomes dependent on the stability of the attached variant. Plasmids encoding these fusions are transformed into cells, and the activity of the enzyme is then selected for. The selection results in the enrichment or depletion of each plasmid (5) based on the stability of the fusion protein that it encodes, which in turn determines the cell's growth rate. The frequency of each plasmid in the population before and after selection is determined by DNA sequencing. The ratio of the selected frequency to the input frequency, called the enrichment score, serves as a proxy for in vivo stability (Fig. 1A).
We demonstrate the Stable-seq method via the analysis of a degron, which is a protein degradation signal recognized by the proteolytic machinery (6). We fused the well-characterized degron Deg1 from the yeast Mat␣2 protein (7) to the yeast Leu2 protein, which is necessary for leucine biosynthesis (Fig.  1B). Mat␣2 and Mata1 are transcriptional repressors required in order to specify mating type in Saccharomyces cerevisiae. When Mat␣2 forms a heterodimer with Mata1 in diploid a/␣ cells through its C-terminal tail, it is relatively stable, but it becomes short-lived in haploid ␣ cells through degradation by the ubiquitin proteasome system (8). Deg1, which spans the N-terminal 67 residues of Mat␣2, is recognized by the E3 enzyme Doa10 (9). In the fusion protein, the stability of Leu2 becomes dependent on Deg1. Thus, any mutation in Deg1 that increases its stability results in the presence of more Deg1-Leu2 fusion protein and increased production of leucine. The additional leucine leads to increased growth of yeast, and thus more copies of the Deg1 sequence containing this mutation. Using this approach, we analyzed the effect of ϳ30,000 mutations in Deg1 and identified Deg1 features that affect stability. This approach should work in other organisms with appropriate selections, and it could be scaled up to measure the stability of many different proteins in parallel.

EXPERIMENTAL PROCEDURES
Plasmids and Strains-The p416TET off -Deg1-LEU2 plasmid was constructed as follows. The SacI site 5Ј of the GPD1 promoter in p415GPD was changed to an NheI site to replace the GPD1 promoter with the TET off promoter cassette from pCM182 (tTA transcriptional activator, CMV promoter, ADH1 terminator, TetO operator, CYC1 promoter, Multiple Cloning Sites, and CYC1 terminator) (10). The LEU2 and Deg1 (Mat␣2 1-67 ) sequences were sequentially cloned into the Multiple Cloning Sites to generate p416TET off -Deg1-LEU2. A linker (encoding PRRSG) is present between Deg1 and LEU2. Sitedirected mutagenesis was carried out for synonymous codon changes in Deg1 to generate a HindIII site at residue 38 and at residues 42-45 to optimize the annealing temperature for an Illumina sequencing primer ( Fig. 1B and supplemental Table S1). Deg1-LEU2 M1⌬ was generated via site-directed mutagenesis to remove the codon for the first methionine in LEU2. FLAG epitope-tagged versions of LEU2, Deg1-LEU2, and Deg1-LEU2 M1⌬ were cloned into the plasmid with the p416TET off promoter. Point mutant constructs in Deg1 described in Figs. 3C, 5B, and 6B and in supplemental Figs. S2B and S2C were generated via site-directed mutagenesis. All the clones were confirmed by Sanger sequencing. Detailed plasmid maps and sequences are available upon request. Yeast strain BY4741 and isogenic yeast deletion strains are from Open Biosystems.
Construction of the Deg1 Mutant Library-An oligonucleotide encoding amino acids 3-34 of Deg1 was synthesized by Trilink Biosciences (Deg1 3-34 library). In the synthesis, the variable region was doped with 2.088% (0.696% of each non-wild-type nucleotide) to generate on average two nucleotide changes per clone. The oligonucleotide was made double-stranded and then PCR amplified using iProof TM High-Fidelity DNA Polymerase (Bio-Rad) for 15 cycles. Gel- purified PCR fragments were digested with NotI and HindIII and used to replace the wild-type sequence in p416TET off -Deg1-LEU2 M1⌬ to generate a library of ϳ170,000 variants of the Deg1 sequence (input library). The Deg1 N2 library was constructed with an oligonucleotide containing random bases (NNN) at the second codon (Asn) of Deg1, which was cloned into p416TET off -Deg1-LEU2 M1⌬ as described for the Deg1 3-34 library.
Deg1 Library Selection and High-throughput Sequencing-Plasmid DNA of the Deg1 3-34 library was transformed into yeast strain BY4741 via the lithium acetate method to generate 1.1 ϫ 10 6 transformants (11,12). Only moderate transformation efficiency (Ͻ0.5%) was achieved in an effort to avoid the co-transformation of two variant plasmids into a single cell, which was determined to occur Ͻ2% of the time. Approximately 70,000 yeast transformants were plated onto 20 SC-Leu-Ura (13) 15-cm plates and incubated at 30°C for 3 days. Colonies were scraped off the plate and used to prepare DNA via the zymolase method. Briefly, yeast cells were lysed with Qiagen Buffer P1 supplemented with 50 mM DTT and 400 g/ml zymolase 20T at 37°C for 2 h with occasional shaking, and lysing was followed by incubation at Ϫ80°C for 30 min and thawing at 42°C for 1 min. The QIAprep Spin Miniprep Kit (Qiagen) protocol was used to elute plasmid DNA contaminated with yeast genomic DNA. Genomic DNA was digested with Exonuclease I (Affymetrix) and lambda exonuclease (New England Biolabs) and removed by a Zymo DNA Clean & Concentrator TM -5 (Zymo Research) to enrich plasmid DNA in the eluate (selection library). To count the number of plasmids in each pool, we amplified the variable region of Deg1 from input and selection libraries via PCR using iProof TM High-Fidelity DNA Polymerase (Bio-Rad) for 15 cycles, and we sequenced this fragment by MiSeq and HiSeq2000 (Illumina) using the primers listed in supplemental Table S1. The Deg1 N2 library was assayed via same procedures as the Deg1 3-34 library.
Analysis of Sequence Data-The identity of each mutation and its frequency in the input and selected libraries were determined by the Enrich software package (14). We used the E. coli plasmid library of Deg1-LEU2 variants as the input library, because for yeast transformants, even in SC-Ura media (containing 80 mg/l leucine), there was a growth advantage to yeast cells that produced more leucine from a stabilized Deg1. To ensure the quality of sequencing reads, we used paired-end sequencing to read both directions and applied a stringent quality filter. The quality filter included (i) a minimum quality score higher than 20 at every cycle, (ii) fewer than four consecutive mutations, (iii) no ambiguous sequencing bases, and (iv) minimum read counts in the input library of at least 15. The enrichment score (E) of each mutant was calculated as R m /R w , with R m being the frequency of the mutant (m) in the selection library divided by its frequency in the input library, and R w the frequency of the wild type (w) in the selection library divided by its frequency in the input library. This normalization avoids the bias due to different plasmid frequencies in the input population. Separate aliquots of cells from the selection plates were prepared for sequence analysis by either MiSeq or HiSeq2000. Read counts of each variant from the MiSeq and HiSeq2000 runs were highly correlated (R 2 of 0.9999 and 0.9816 for the input and selection libraries, respectively). Therefore, the sets of sequencing data were merged for further analyses. Epistasis scores were calculated as log 2 E m1,2 Ϫ (log 2 E m1 ϩ log 2 E m2 ), where m1,2 represents a double mutant (m1, first single mutant; m2, second single mutant). Epistatic interactions (supplemental Fig. S3) were visualized with a custom script using the D3.js java script library. Other computational and statistical analyses were performed with in-house python scripts and the R statistical package.
Yeast Cell Growth Assay-To determine the range of in vivo protein stability assayed with Deg1-Leu2 fusions, we monitored the growth rates of yeast cells containing variant fusion constructs in liquid culture over 60 h (Synergy H1, BioTek, Winooski, VT). For the spotting assay, yeast cells containing variant Deg1-LEU2 constructs were grown in SC-Ura media, and equal amounts of cells (OD ϭ 0.5) were spotted onto the plates with 5-fold serial dilutions. The plates were incubated at 30°C for 1 to 3 days.
Immunoprecipitation and Western Blotting-Equal amounts of yeast cells containing Deg1-LEU2-FLAG constructs were harvested from an exponentially growing culture in a nonselective media (SC-Ura) and lysed with a bead-beater. Cell extracts were equalized using the Bradford assay, and FLAG-tagged Deg1-Leu2 proteins were immunoprecipitated with FLAG-M2 agarose beads (Sigma A2220). The immunoprecipitates were separated via 10% SDS-PAGE, and a Western blot was visualized using anti-FLAG (Sigma F3165), True-Blot ® (Rockland Immunochemicals Inc., Gilbertsville, PA) to reduce the signal from the IgG heavy chain (which is similar in size to Deg-Leu2-FLAG), and ECL reagents (GE Healthcare RPN2106).

Strategy for Measuring the Stability of a Protein by Fusing
It to a Reporter Protein-We first sought to confirm that fusion of Deg1 to Leu2 generates a fusion protein with a rapid turnover. We tested the growth rate of cells containing Deg1-Leu2 variants via a spotting assay in which 5-fold serial dilutions of a yeast culture were plated on control (-uracil) or selection (-leucine -uracil) plates ( Fig. 2A). Yeast expressing Leu2 with no Deg1 sequence grew well under selection, but A, spotting assay of Deg1-Leu2 variants with 5-fold serial dilutions. Growth on the -Ura plate, which requires only the presence of the URA3 transformation marker, serves as the spotting control, and growth on the -Leu -Ura plate selects for stable versions of Leu2. B, Western blot analysis of C-terminally FLAG-tagged Deg1-Leu2 variants in DOA ϩ and doa10⌬ cells. The full-length Deg1-Leu2 and Deg1-Leu2 M1⌬ proteins are unstable in DOA ϩ cells, but Deg1-Leu2 produces a Leu2-sized band. Both Deg1-Leu2 and Deg1-Leu2 M1⌬ produce a full-length band in doa10⌬ cells, as well as a smaller band that runs between Deg1-Leu2 and Leu2 that is likely due to cleavage of the full-length protein.
the expression of Deg1-Leu2 resulted in only moderately reduced growth under selection. We surmised that this modest reduction was due to translation of the fusion downstream from the initiator codon. Although Deg1 contains no methionine except for the initiator, translation could begin from the next in-frame methionine, the start codon of Leu2, especially given the poor context (15) for the Deg1 initiator in this fusion construct. Deletion of the first methionine of Leu2 (Deg1-Leu2 M1⌬ ) nearly eliminated the growth of yeast under selection ( Fig. 2A), indicating that Leu2 had been made unstable by its fusion to Deg1.
The use of the alternative start codon likely occurred by leaky scanning of the ribosome (15). Initiation codon selection by the eukaryotic ribosome is often determined by the context surrounding the first AUG of the open reading frame. In S. cerevisiae, a 5Ј-untranslated region rich with A's, especially an A at position -3, is highly favored (15). When this context is not favorable, the next AUG in a better context may be used as an alternative start codon. In the case of Deg1-Leu2, an unfavorable GCGGCCGC precedes the first AUG. The use of the LEU2 AUG was confirmed by Western blot analysis of FLAG epitope-tagged Deg1-Leu2 variants (supplemental Fig. S1A).
In DOA ϩ cells expressing Deg1-Leu2, only a band the size of Leu2 was apparent, whereas these cells expressing Deg1-Leu2 M1⌬ showed no detectable Leu2 band (Fig. 2B), consistent with the inability of these cells to grow in selection media (supplemental Fig. S1A). In doa10⌬ cells, in which the degron is not targeted for degradation, the Deg1-LEU2 plasmid produced a protein consistent with translation starting from the initiator methionine of Deg1 and another protein the size of Leu2; however, the Deg1-LEU2 M1⌬ plasmid produced the larger species but no Leu2-sized protein (Fig. 2B). The degradation of the Deg1-Leu2 fusion was dependent on both an E3 (Doa10) and an E2 (Ubc7), like other Deg1 fusion proteins (8,9) (supplemental Fig. S1B and Fig. 2B). Given the failure of Deg1-LEU2 M1⌬ to provide sufficient Leu2 function, stable variants of Deg1 should result in the production of more leucine, faster growth of yeast on selection media, and thus more copies of the plasmid encoding these variants. We demonstrated that even highly stable variants do not saturate the assay (supplemental Fig. S1C).
Stable-seq Analysis of a Deg1 Mutant Library-To apply the Stable-seq method to many Deg1 variants simultaneously, we used a doped oligonucleotide to mutate residues 3 to 34 of Deg1 fused to Leu2 M1⌬ (Fig. 1B), generating a library of ϳ170,000 Deg1 variant plasmids. Yeast cells transformed with this library formed colonies of similar size if no selection for leucine was imposed, but the transformants produced far fewer colonies, of different sizes, on a selection plate (Fig. 3A). We interpret the leucine selection results as support for the rationale that colony size is dependent on the amount of Deg1-Leu2 M1⌬ fusion protein, which in turn is determined by the stability of Deg1. We harvested the cells from selection plates, isolated plasmids, and sequenced the DNA encod-ing Deg1. Comparing the frequency of each variant in the selected yeast to that in the input plasmid library allowed us to assay the effect of mutations on the stability of ϳ30,000 variants of Deg1 (supplemental Table S2). Mutations present in the input library but not present after selection might be the result of extreme instability, or they may have been lost due to the limited number of colonies sampled after selection. Fig. 3B shows the log 2 enrichment scores of the single mutations observed, covering 71% of all the possible single mutations. Of the 13 previously identified mutations (7) for which we had DNA sequence data, 10 had positive log 2 enrichment scores (Ͼ1.5). These scores indicate that the mutations increased in frequency after selection, in accordance with their behavior in the spotting assay (Fig. 3C). The Deg1 residues in which these previously identified mutations occur in several cases could also be mutated to other amino acids with similar or greater enrichment scores (e.g. F18D, F18N, S20N, I22G, L29D, and I32D). In addition to mutations in the previously identified residues, we also found novel mutations with high enrichment scores, including D8R, S21Q, and K27S (Fig. 3B).
In order to compare the scores obtained via Stable-seq with those determined via a previous approach, we examined the growth under leucine selection of cells expressing Deg1-Leu2 that contained one of 13 characterized mutations in Deg1 that lead to stabilization, which had been identified based on the ␤-galactosidase activity of a Deg1-␤-galactosidase fusion (7). These mutations resulted in better growth on the selection plate for cells carrying 10 of these variants (Fig. 3C). The lack of correlation for the other three variants (K19Q, S20P, and I32S) may be due to the different behavior of the Deg1-Leu2 versus the Deg1-␤-galactosidase fusion protein in different expression systems and assays.
Unlike the previous small-scale study (7), we also identified mutations with negative log 2 enrichment scores. These less stable mutants included more than 60% of all the single mutants observed (supplemental Fig. S2A). For example, C33F showed greater instability than wild-type Deg1 (supplemental Fig. S2B), but the degradation was still dependent on Doa10 (supplemental Fig. S2C). These variant Deg1 sequences may be better recognized by the E3 ligase Doa10 than the wild-type Deg1.
Alternative Start Codons and Their Effect on Protein Stability-Changes to methionine at 10 positions, between residues 14 and 32, showed a strong stabilizing effect (Fig.  4A). The sequencing data for several mutants containing a new methionine were also confirmed by a spotting assay (Fig.  4B). It is likely that the new methionines serve as alternative start codons via a leaky scanning mechanism (15). The effect on the stability of truncated Deg1 variants due to alternative start codons correlated well with the results of a previous deletion study (7). Support for this interpretation also comes from 33 double mutants that combine a stop codon with a new methionine C-terminal to this stop (Fig. 4C), which likely initiate or reinitiate translation at the new methionine. Of these double mutants, 21 had a positive log 2 enrichment score. The location of the new methionine, or the distance between the stop codon and the new methionine, did not correlate with the enrichment score (data not shown), but the double mutants with a stop codon at residue 17 followed by a methionine strongly stabilized Deg1, suggesting that the underlying mechanism is more complex.
N-terminal Processing Effects on in vivo Protein Stability-N-terminal acetylation is the major post-translational modification in eukaryotes, with more than 50% of the proteins in S. cerevisiae undergoing this modification (16). At the N terminus of a protein, the initiator methionine is removed if the second residue is small enough (e.g. Ala, Ser, Thr, Val, or Cys) to be accessed by methionine aminopeptidases, and the exposed second residues are acetylated by the NatA complex (16,17). Initiator methionines followed by a larger residue are often acetylated by other N-terminal acetyltransferases (e.g. NatB or NatC) depending on the property of the second residue (16,17). N-terminal acetylation has been proposed as another type of degradation signal based on work using Deg1 fusion proteins as model substrates (18). This concept was also explored with physiological substrates (19,20), which further expanded the functions of the N-end rule pathway (21,22).
To determine whether the Stable-seq method could detect changes that affect N-terminal processing, we generated another library that had a random nucleotide sequence specifying only the second residue (Asn) of Deg1 (Deg1 N2 library), and we subjected the library to the same assay as the Deg1 3-34 library. We analyzed all 20 amino acids, as well as the stop codon, and found that the strongest increases on stability were due to the mutations N2K and N2R (Figs. 5A and 5B), consistent with a lack of N-terminal acetylation when the second residue is basic (18).
Acetylation of the wild-type Deg1 occurs via the action of NatB (18). Nat3 is the catalytic subunit of the NatB N-terminal acetyltransferase complex that acts on Mat␣2. Deletion of the NAT3 gene resulted in stabilization of the Deg1-Leu2 M1⌬ fusion (Fig. 5C), as observed for another Deg1 fusion (18). This stabilization is consistent with the role of  (7). The mutations identified by Johnson et al. (7) resulted in increases in stability, measured by ␤-galactosidase or pulse-chase assays, as shown (␤-gal). The log 2 E values are enrichment scores calculated from DNA sequence data: variant frequencies after leucine selection were divided by frequencies in the input library and then normalized to the wild-type ratio of frequencies.
N-terminal acetylation in protein degradation as proposed by Hwang et al. (18).
Epistatic Effects Observed in Double Mutants-In addition to yielding single mutants for analysis, the Stable-seq method also generated data for Ͼ17,000 variants (58% of the total variants observed) that contain two mutations. Double mutants can be examined for epistasis, in which the interaction between two mutations causes the double mutant to behave unexpectedly given the behavior of the two constituent single mutations. We used a multiplicative predictive model in which epistasis scores were calculated by subtracting the sum of the log2 enrichment scores of each single mutant (predicted stability) from the log2 enrichment score of the double mutant (observed stability) (Fig. 6A). Based on this model, positive epistasis indicates that the double mutant displayed more stability than predicted, and negative epistasis the opposite. We examined a few double mutants that showed positive epistasis via a spotting assay and validated these unexpected increases in stability, including cases in which neither single mutation alone resulted in an increase in stability (Fig. 6B).
We further analyzed the double mutants with the most positive epistasis scores (highest 1%) and the most negative epistasis scores (lowest 1%) (supplemental Fig. S3). Positive epistasis in the highest 1% was distributed among many different mutations, with no single mutation accounting for a large fraction of the total number of epistatic interactions. However, much of the negative epistasis in the lowest 1% could be accounted for by a small number of mutations. In most of these cases, the strong stabilizing effect of one mutation in this small grouping was not further increased by the presence of many different second stabilizing mutations. Some of the mutations in the small grouping are changes to a methionine, and others occur in the set of 13 previously identified mutations (7).
Clustering of Read Counts-We analyzed the effect of the proposed mechanisms of stabilization for all the single mutations with positive log 2 enrichment scores. These mutations were clustered into five groups, with the median enrichment scores and ranges of each shown by a boxplot (Fig. 7A). Mutations that generate a new methionine or that had been previously identified in the study by Johnson et al. (7) resulted in the highest median enrichment scores. Although there are relatively few mutations in these two groups, they account for more than half of the sequence reads of enriched mutations from the selection library (Fig. 7B). DISCUSSION Here, we provide a method, Stable-seq, that uses highthroughput DNA sequencing to assess in vivo protein stability. We show that Stable-seq can identify key features of a protein domain that affect stability, including mutations that affect the full-length domain, alternative start codons that likely truncate the domain, and mutations that appear to affect N-terminal acetylation. The strong correlation between the scores generated via DNA sequencing and the stabilities of Deg1-␤galactosidase variants determined via ␤-galactosidase assay or pulse-chase analysis (7) indicates that the high-throughput Stable-seq assay is measuring stabilities in a useful range. Moreover, Stable-seq does not require the use of multiple time-points to calculate protein stability as other methods do, and its use of DNA sequencing reactions to compare input and selected populations allows a fine-grained discrimination of protein stability.
Changes of internal residues to a methionine can serve as an alternative start codon. Through this process, the protein becomes truncated at the N terminus, which for Deg1 results in protein stabilization, because the degron is no longer functional. This proposed mechanism is supported by double mutants that contain an upstream termination codon followed by a new methionine. We also found that the identity of the second residue affects stability. This mechanism is likely mediated by whether or not N-terminal acetylation occurs, as the Deg1-Leu2 M1⌬ fusion is stabilized in the nat3 mutant, which does not carry out NatB-mediated acetylation. Both of these mechanisms may interact together in complex ways, as each new methionine that serves as an alternative start codon is coupled to a new residue in the second position.
We identified double mutants that showed much greater or lesser stability than would be expected based on the behavior of the constituent single mutants. A mutation to a methionine generally showed negative epistasis when it combined with another mutation at an upstream location (supplemental Fig.  S3, bottom panel), indicating that the double mutant was less stabilized than expected. These results support the idea that the new methionine serves as an alternative start codon, because the N-terminally truncated Deg1 that initiates from the new methionine would not contain the upstream stabilizing mutation.
Stable-seq is based on the assumption that the transcription and translation of variants are the same, and thus the level of the nutritional marker is dependent solely on the stability of the variants. However, it is possible that other factors will influence the function of the metabolic enzyme (e.g. Leu2); for example, mutations could change the folding of the enzyme or protein-protein interaction. When the method is applied to assay diverse proteins simultaneously, the presence of factors such as different protein localization signals will make additional controls necessary in order for the resultant stabilities to be validated.
Stable-seq could be adapted to analyze the stability of other degrons and other proteins. The method could be scaled up to handle large complements of proteins simultaneously if libraries of Leu2 fusions with random genomic or cDNA inserts or a collection of defined open reading frames were assayed. Stable-seq should also be amenable to other model organisms or to tissue culture cells if appropriate selection markers such as proteins that confer drug resistance are used.