[go: up one dir, main page]
Next Article in Journal
Multiple Factors Drive Replicating Strand Composition Bias in Bacterial Genomes
Next Article in Special Issue
Combining Protein and Strain Engineering for the Production of Glyco-Engineered Horseradish Peroxidase C1A in Pichia pastoris
Previous Article in Journal
Critical Issues in the Study of Magnesium Transport Systems and Magnesium Deficiency Symptoms in Plants
Previous Article in Special Issue
Effects of Non-Natural Amino Acid Incorporation into the Enzyme Core Region on Enzyme Structure and Function
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 16
Issue 9
10.3390/ijms160923094
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Rational Protein Engineering Guided by Deep Mutational Scanning

by
HyeonSeok Shin
and
Byung-Kwan Cho
*
Department of Biological Sciences and KI for the BioCentury, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 305-701, Korea
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2015, 16(9), 23094-23110; https://doi.org/10.3390/ijms160923094
Submission received: 14 August 2015 / Revised: 4 September 2015 / Accepted: 13 September 2015 / Published: 23 September 2015
(This article belongs to the Special Issue Protein Engineering)
Browse Figures
Figure 1
<p>Schematics of the simplified overview of deep mutational scanning: (<b>a</b>) generation of the initial protein variant library for target protein sequences; (<b>b</b>) screening for protein variants with desired properties; and (<b>c</b>) sequencing and quantification of the mutations under different selection pressures. The asterisks indicate mutations at a specific site and the stacked asterisks indicate enrichment of mutations in specific sites after quantification. For example, mutation counts at different sites are shown with <b>*</b> positions carrying a mutation, <b>**</b> positions carrying two mutations, and <b>****</b> positions carrying four mutations.</p> "> Figure 2
<p>Schematics of the simplified screening systems of the (<b>a</b>) conventional strategy of directed evolution, where iterative assays are performed until a desired phenotype appears; and (<b>b</b>) deep mutational scanning, where the protein variants are screened to a simpler selection pressure. The different phenotypes of the protein variants are shown by gradient of green colored wells. The desirable phenotypes are shown by darker green colored wells and loss of function is shown by white colored wells.</p> "> Figure 3
<p>Methods to rectify sequencing errors: (<b>a</b>) a scheme of how paired end reads with short sequencing reads allow the detection of sequencing errors; (<b>b</b>) schematic showing the concept of the duplex sequencing method; (<b>c</b>) how the consensus sequence is used to remove sequencing errors, adapted by permission from the Macmillan Publishers Ltd: <span class="html-italic">Nature Protocols</span> [<a href="#B61-ijms-16-23094" class="html-bibr">61</a>], copyright 2014. The black bar indicate the target inserts reads for sequencing and the orange and dark blue colored bars at end of the insert reads indicate sequencing adaptors; the yellow and light blue bars indicate the randomized duplex tags; and (<b>d</b>) Hypothetical mapping of the mutation frequency for variant library sequencing. The red line indicates the sequencing error rate of the ampicillin gene used as the cutoff. The asterisks indicate mutations and the orange and blue bars at the ends of the reads indicate the sequencing adaptors.</p> "> Figure 4
<p>(<b>a</b>) A hypothetical mutational map generated to show mutation frequency at each position; Part of the mutational map showing (<b>b</b>) extremely tolerant and critical residues to mutations; (<b>c</b>) tolerant to hydrophobic mutations and (<b>d</b>) tolerant to hydrophilic mutations. The x-axis indicates the protein residues and the y-axis indicates the possible amino acids. The color key represents the mutation frequency at each amino acid. The white color indicates that no mutation was found. The blue color indicates mutation frequency of loss-of-function variants and red color indicates mutation frequency of function-retained variants. The stop codon is indicated by <b>*</b>.</p> ">
Versions Notes

Abstract

:
Sequence–function relationship in a protein is commonly determined by the three-dimensional protein structure followed by various biochemical experiments. However, with the explosive increase in the number of genome sequences, facilitated by recent advances in sequencing technology, the gap between protein sequences available and three-dimensional structures is rapidly widening. A recently developed method termed deep mutational scanning explores the functional phenotype of thousands of mutants via massive sequencing. Coupled with a highly efficient screening system, this approach assesses the phenotypic changes made by the substitution of each amino acid sequence that constitutes a protein. Such an informational resource provides the functional role of each amino acid sequence, thereby providing sufficient rationale for selecting target residues for protein engineering. Here, we discuss the current applications of deep mutational scanning and consider experimental design.

1. Protein Engineering in the Ultrahigh-Throughput Sequencing (uHTS) Era

Protein engineering has extensively been applied to several fields of biotechnology, including the medical sciences and for industrial applications. Since the first report of protein engineering by site-directed mutagenesis in 1982, numerous proteins have been engineered using various techniques developed by technological advances [1]. Advances in protein structure determination techniques have had a great impact on protein behavior prediction and mechanisms, which were integrated with the rational strategy of protein engineering as an informational resource [2,3,4]. Such protein structure data have allowed the developments of modeling methods such as Rosetta and RosettaDesign, which enables the prediction of protein structure for novel as well as engineered proteins [5,6]. Improved computational and modeling power in addition to accumulation of protein structure and mechanism information have introduced a semi-rational method of protein engineering that uses predictive algorithms to preselect potential target sites [7,8]. Developments in high-throughput screening systems have had a great impact on the directed evolution strategy of protein engineering. Developments in application of high-throughput screening methods such as enzyme-linked immunosorbent assay and flow cytometers, such as fluorescence activated cell sorting (FACS), have reduced the time and effort involved in protein engineering [9,10]. Thus far, protein engineering, as a field, has grown along with development of other biotechnologies and it is only natural that recent advances in high-throughput sequencing are applied to protein engineering [11].
Integration of high-throughput sequencing technology with protein engineering involves the coupling of protein variant generation with a high-throughput screening system. Instead of selecting a few selected protein variants, high-throughput sequencing allows the sequencing of millions of protein variants, termed as deep mutational scanning [11]. Depending on the selection pressure, both positive and negative phenotypes by protein mutations can be linked to the sequence space [12]. Such an information load with regard to the protein sequence space has a great impact on the field because it solves some of the problems that occur as a part of current protein engineering strategies. For rational and semi-rational strategies of protein engineering, extensive information on protein structure is required for selecting a target site for engineering. Furthermore, such a selection or pre-selection algorithm shows difficulty in prediction of sites that are distant from binding or active sites [13,14]. On the other hand, directed evolution is another protein engineering strategy comprising two-steps: (1) mutagenesis to generate the mutant library of the protein and (2) screening for the protein variant with the desirable property. For example, phage display has been used to display proteins or peptides on the phage surface and, when followed by affinity captures, allows selection of proteins with desired properties without much prior knowledge on the protein structure or mechanism of function [15,16,17]. Thus, structural or functional information of a protein is not necessary, provided that selection pressure or screening method for desired phenotype is available [18]. However, this approach can only screen for a limited number of engineered products compared to the millions of possible mutations that are contained in the initial mutant library due to the requirement for DNA sequencing. From this perspective, deep mutational scanning provides the required depth for sequencing millions of possible mutations, which, in turn, provides information on the important sites, which is not necessarily when the distance to the known active sites is close. Furthermore, in contrast to the conventional method of directed evolution, where iterative mutagenesis and screening are performed until the desired phenotype is found, the screening system used in deep mutational scanning can be as simple as function on/off that saves time and effort.

2. Deep Mutational Scanning

2.1. Overview

Deep mutational scanning as a concept is quite exciting. It identifies the effect of all possible amino acid changes for each position in a protein and compares the enrichment of mutation to appropriate selection pressure depending on the characteristics of the protein. This is enabled by the delivery of the information between thousands or millions of mutations to protein variants’ phenotype by high-throughput sequencing. Compared to the initial protein variant library, the mutations quantified after selection enable identification of the effect of mutations in certain residues to the functional change of the protein variant. Overall, the elucidation of such vast information on protein sequence space in relation to selection pressure aims to determine the protein sequence and function relationship [11]. In addition, the analysis of the mutations in the protein variants reveals a sequence–structure–function relationship, which is difficult to predict [19,20]. In terms of strategy, deep mutational scanning uses the discovery-based protein engineering methods to determine the protein sequence and function relationship by linking the mutational genotype to the phenotype, which then provides the basis for hypothesis-driven protein engineering.
Because deep mutational scanning is a relatively new concept that differs from conventional protein engineering strategies and experimental designs, this review introduces the basic overview of the methodologies that help determine the experimental design as well as the computational analysis using recent studies. Precautions for data interpretation are also discussed along with perspectives on the current limitations and future developments.
To design a deep mutational scanning experiment for a target protein, three conceptual steps that involve construction of the protein variant library, screening, or selection for functions, and high-throughput sequencing must be considered. First, a library containing the mutated sequence of the target protein is generated, which is used as the input variant library (Figure 1a); Second, the input variant library is expressed and subjected to a screening system that can concordantly link the sequence variant to the functional variant (Figure 1b); Third, the selected protein variants are subjected to high-throughput sequencing and the mutation spectra in different libraries are quantified (Figure 1c). Although Figure 1 illustrates the simplified overview of the processes involved in deep mutational scanning, different methodologies and schemes may be applied. The mutational library can be generated by random mutagenesis, saturated mutagenesis, or DNA synthesis. The selection pressure in deep mutational scanning depends on the protein’s characteristic and the experimental design. It can be applied to assays of binding affinity, enzymatic assays, as well as fluorescence signals [21,22,23,24].
Ijms 16 23094 g001 1024
Figure 1. Schematics of the simplified overview of deep mutational scanning: (a) generation of the initial protein variant library for target protein sequences; (b) screening for protein variants with desired properties; and (c) sequencing and quantification of the mutations under different selection pressures. The asterisks indicate mutations at a specific site and the stacked asterisks indicate enrichment of mutations in specific sites after quantification. For example, mutation counts at different sites are shown with * positions carrying a mutation, ** positions carrying two mutations, and **** positions carrying four mutations.
Figure 1. Schematics of the simplified overview of deep mutational scanning: (a) generation of the initial protein variant library for target protein sequences; (b) screening for protein variants with desired properties; and (c) sequencing and quantification of the mutations under different selection pressures. The asterisks indicate mutations at a specific site and the stacked asterisks indicate enrichment of mutations in specific sites after quantification. For example, mutation counts at different sites are shown with * positions carrying a mutation, ** positions carrying two mutations, and **** positions carrying four mutations.
Ijms 16 23094 g001

2.2. Mutagenesis

The desired sequence space and the size of the protein must be considered for the initial mutant library generation. For instance, a hypothetical protein “hyp-A” constituted by 100 amino acids would require the generation of 1900 different variants for single site mutagenesis of all residues that constitute the protein. To obtain an unbiased mutant library that constitutes all possible sequence space requires synthesis of 1900 oligonucleotides or 3800 primers. The complete coverage of initial library over the sequence space and unbiased concentration of each mutant is important because the initial library is also sequenced as a control [22]. To validate that a mutation site in a protein is important because it has been enriched after selection, the initial library data must show that the mutation was evenly distributed. Thus, a protein variant library generation method such as single-site saturation mutagenesis (S-SSM) is very effective in terms of the completeness and uniformity, as all 20 possibilities of amino acid changes for each residue is designed [23,24,25,26].
Oligonucleotide-directed random mutagenesis is another mutagenesis method frequently used in deep mutational scanning [10,22,27,28,29,30,31,32,33,34,35,36]. This method uses degenerative oligonucleotides that contain the “NNS” or “NNK” codon, where N represents any of the four-nucleotide sequences, S represents either G or C and K represents T or G [10,30]. This method is effective in that all sequence spaces of the target protein can be theoretically covered with one order of oligonucleotide synthesis compared to S-SSM, where 19 oligonucleotide are required [37]. However, there is also a possibility that not all sequence spaces maybe covered in the initial library. As oligonucleotides are randomly synthesized, the possibility of an amino acid being generated from the random combination of NNS or NNK is not equal. There is a chance that the generated NNS sequence may not be evenly distributed and that some of the mutations may be lost in the transformation or expression system. To ensure the complete sequence space of all amino acid variants are covered, the depth required for screening must be carefully calculated. From a practical viewpoint, oligonucleotide-directed random mutagenesis is cost efficient but there is a possibility of some mutations being lost from the initial library. Therefore, while S-SSM is relatively cost-intensive and time consuming, the generated library would be complete and unbiased.
Another method for generating a protein variant library is random mutagenesis by PCR [21,38,39,40,41]. This method can be considered for target proteins or sequences that too are long to use oligonucleotide-directed method, as it is the most efficient in terms of cost and time. However, PCR-based mutagenesis does not cover all possible sequence spaces. Depending on the initial nucleotide of a codon, only 12 mutations out of 19 are possible unless a consecutive mutation is introduced in the codon. Furthermore, random mutagenesis by PCR is known to have some bias in the mutation spectra. Depending on the mutation frequency and the coverage of the mutant library, such a bias may be negligible, as less biased PCR mutagenesis has been developed [42]. For example, Shin et al. used PCR-based random mutagenesis, which showed that the mutation frequency of AT-to-NN was similar to GC-to-NN [21]. Overall, each method for generating the initial mutation library has its advantages and disadvantages and the best fitting method may differ from case to case. Furthermore, although this review followed the sequential steps of experiment, the method to generate the initial mutant library should consider all aspects of the process such as the expression system, selection pressure, and sequencing methods.

2.3. Construction of a Protein Variant Library

The expression system for deep mutational scanning is not too different from the conventional protein engineering methods that require a selection system. Similar to any directed evolution strategy in protein engineering, the expression system requires a link between the mutated DNA and the protein variant that enables determination of the corresponding mutation responsible for the changed phenotype [43]. On the contrary, the selection pressure for deep mutational scanning is quite different from that in conventional protein engineering. A directed evolution strategy of protein engineering involves performing repeated screening assays until the desired phenotype is found (Figure 2a). Such repeated assaying is costly, time consuming, and laborious since there is no guarantee that the desired phenotype will appear. On the other hand, deep mutational scanning uses a different strategy in screening the protein variants (Figure 2b). Instead of performing iterative assays to find the protein variants with the desired phenotype, a simple screening pressure such as function on/off or binding affinity is used to divide the selected and unselected protein variants. The massive amount of information from loss of function variants is also analyzed to find sites critical to the mutation or to protein fitness [44].
Recent studies have shown that three types of expression systems are frequently used in deep mutational scanning—the plasmid system, phage or yeast display system, and bacterial or yeast two-hybrid system (Table 1). Among these expression systems, the phage display system is used for the investigation of proteins known for their protein–peptide and protein–DNA binding interactions, where assays using binding affinity are favorable [45]. For example, Fowler et al. used the phage display system to investigate the affinity of binding of the WW domain (named by the two conserved tryptophan residues in the domain) to its peptide ligand [22]. Using beads with the peptide ligands attached, six successive rounds of selection and washing were performed, which allowed the variant WW domain with better affinity to be enriched. The advantage of this display system is that both the protein and the encapsulated DNA are easily accessible. As the protein is displayed in the exterior part of the phage, the use of cytoplasmic phages such as T7 allows efficient delivery of the proteins [46]. Thus, the phage display system is an effective combination that can be used in functional protein binding assays in deep mutational scanning.
Ijms 16 23094 g002 1024
Figure 2. Schematics of the simplified screening systems of the (a) conventional strategy of directed evolution, where iterative assays are performed until a desired phenotype appears; and (b) deep mutational scanning, where the protein variants are screened to a simpler selection pressure. The different phenotypes of the protein variants are shown by gradient of green colored wells. The desirable phenotypes are shown by darker green colored wells and loss of function is shown by white colored wells.
Figure 2. Schematics of the simplified screening systems of the (a) conventional strategy of directed evolution, where iterative assays are performed until a desired phenotype appears; and (b) deep mutational scanning, where the protein variants are screened to a simpler selection pressure. The different phenotypes of the protein variants are shown by gradient of green colored wells. The desirable phenotypes are shown by darker green colored wells and loss of function is shown by white colored wells.
Ijms 16 23094 g002
Table
Table 1. Recent studies using deep mutational scanning.
Table 1. Recent studies using deep mutational scanning.
Mutation Generation Method 1Variant LibrarySequencing Method 2Target Protein 3Reference
ORMPhage displaySolexa/PEPSD95pdz3[10]
ORMBacterial two-hybridIllumina/PEhYAP65[22]
ORMYeast two-hybridIllumina/SEBRCA1[27]
PRMPlasmidIllumina/SEEcFbFP[21]
SMYeast displayIllumina/PEHB80.3[24]
ORMPlasmidIllumina/PEAPH(3′)II[28]
SMPlasmidIllumina/PEBgl3[23]
SMPlasmid454CcdB[26]
ORMPlasmidIllumina/PEPab1[29]
ORMMammalian display vectors454IgG[30]
ORMRibosome display454CDR loops of Fab[47]
ORMPhage displayIllumina/PEhYAP65[48]
1 ORM: Oligonucleotide-directed random mutagenesis; PRM: PCR-based random mutagenesis; SM: Saturated mutagenesis; 2 PE: paired-end; SE: single-end; 3 PDZ domain: post synaptic density protein; hYAP65: human Yes-associated protein 65; BRCA1: breast cancer 1 (early onset); EcFbFP: Escherichia coli flavin mononucleotide binding fluorescent protein; HB80.3: HB80.3 (designed high affinity binding protein); APH(3′)II: Tn5 transposon derived aminoglycoside-3′-phosphotransferase-II; Bgl3: β-glucosidase; CcdB: bacterial toxin protein CcdB; Pab1: poly(A)-binding protein; IgG: immunoglobulin G; CDR: complementary determining region; Fab: fragment antigen-binding region.
The “two-hybrid” system is a screening system used in protein–protein and protein–DNA interactions to activate a downstream gene by binding it to a transcription factor [49,50]. One of the advantages of this system lies in the capability to quantify the downstream reporter genes in vivo. McLaughlin et al. have used the bacterial two-hybrid system to quantitatively link the expression of fluorescent proteins, which enabled FACS to be used as the selection system [10]. Using this system, the authors found nine evolutionarily important sites in PSD05pdz3, which were tested by repeated experiments with different peptide ligands. This study is exemplary in that the discovery-based approach was used to gain information and the gained information was used to successfully apply hypothesis-driven engineering. The two-hybrid system is very effective and efficient in that the reporter system allows rapid quantitative analysis. Although this system requires the preparation of the two-hybrid system for either bacteria or yeast, this system would be a very good candidate for deep mutational scanning in proteins with binding affinity functions.
The plasmid expression system is one of the oldest protein expression systems that is relatively easy to handle [51]. For applications in deep mutational scanning, the plasmid system is favorable for proteins with catalytic activity that require in vivo cell-based screening assays [23,26,28,29]. Unlike the phage display or the two-hybrid method, where the protein’s ligand binding affinity and interaction is measured, proteins with any function can be assayed using the plasmid system. However, besides the advantage of the high degree of freedom in terms of target protein choice, this system does not have any particular advantages in specific screening assays. Thus, the screening system for the protein of choice must be carefully considered for application in deep mutational scanning. For example, the use of a high copy plasmid may cause an overflow of proteins, which would hamper the sensitivity of the screening system; moreover, a saturation point must be measured beforehand.

2.4. Ultra High-Throughput Sequencing (uHTS)

High-throughput sequencing has generated a large amount of sequence data in different applications. In particular, the sequencing of genomes, transcriptomes, ribosome-interacting RNAs, and protein-interacting DNAs has been widely used in countless studies, regardless of the organism or the application [52,53,54]. To apply high-throughput sequencing to protein engineering, a few aspects of the sequencing process must be considered for appropriate data generation and analysis. Among these, sequencing error is a critical aspect that may hamper the accurate interpretation of the sequencing data. For example, the Illumina Genome Analyzer IIX (GA IIX), which has been used in several studies of deep mutational scanning, is known to have an approximate 0.5%–1% error rate [22,48,55]. The hypothetical hyp-A protein-encoding DNA sequence is 300-bp long, which means that the sequencing result may have three errors that are undistinguishable from the actual mutation. To this end, many methods have been developed to reduce errors in the high-throughput sequencing. For example, Lou et al. developed a circle sequencing method that involves denaturing DNA into single-stranded DNA (ssDNA) followed by circularization [56]. By using random primers and Phi29 polymerase that continuously replicate around the ligated circle, many amplified copies of the reads are generated and sequenced to find the consensus sequence with true mutation [57]. Another method to distinguish sequencing error from the true mutation is to use paired-end sequencing with a short read length, which will result in the formation of an overlapping region between the forward and the reverse paired reads. This overlapping region enables the identification of the true mutation, thus significantly reducing the error rate (Figure 3a). Other methods developed to reduce sequencing error rates use tags or barcodes in the sequencing library construction step to find the consensus of sequenced reads [58,59,60]. For example, the duplex sequencing method is known to enable detection of very low frequency mutations by using randomized barcodes at both ends of a read [61,62]. To detect very low frequency mutations, this method enables the detection of both sequencing and PCR errors that result from the construction of the sequencing library. The sequencing library is constructed using the randomized duplex tags of 11 bp downstream and upstream of the sequence, which allows all the reads to have a unique barcode. The reads are then amplified in both the forward and backward manner, which results in at least two sequence reads with the reversed barcode (Figure 3b). After sequencing, the consensus of each of the sequence reads is obtained for the reads with the same barcode (Figure 3c). By comparing the forward and the backward sequences of the reads, the potential errors generated from PCR are eliminated. Finally, comparing the single strand consensus of the same barcodes, but in different directions, allows the detection of the true mutation. Thus, the use of this method can minimize the error rate to a minimum, which would be very effective for deep mutational scanning experiments that require highly sensitive quantification of the mutation.
In addition to the sequencing library generation methods to reduce sequencing and PCR error rate, extra precaution can be taken to further remove any false positives during the quantification of the mutation. For example, Shin et al. used a plasmid system that had an ampicillin resistance gene [21]. Since the E. coli would not survive in the presence of ampicillin without the resistance gene, mutations that had a lower frequency than the sequencing error rate found in the ampicillin-resistance gene were removed from the gene of interest (GoI) (Figure 3d). In this case, some coverage of the possible sequence space would be lost as a result of the loss of the low frequency mutation. However, the primary interest in deep mutational scanning is directed towards the enriched mutational sites, and a cutoff of the mutation frequency would not affect the major sites. Thus, for mutational scanning system results that are expected to show a high frequency for mutations in iterative cycles, a cutoff for mutation frequencies might help in downstream analysis.
Ijms 16 23094 g003 1024
Figure 3. Methods to rectify sequencing errors: (a) a scheme of how paired end reads with short sequencing reads allow the detection of sequencing errors; (b) schematic showing the concept of the duplex sequencing method; (c) how the consensus sequence is used to remove sequencing errors, adapted by permission from the Macmillan Publishers Ltd: Nature Protocols [61], copyright 2014. The black bar indicate the target inserts reads for sequencing and the orange and dark blue colored bars at end of the insert reads indicate sequencing adaptors; the yellow and light blue bars indicate the randomized duplex tags; and (d) Hypothetical mapping of the mutation frequency for variant library sequencing. The red line indicates the sequencing error rate of the ampicillin gene used as the cutoff. The asterisks indicate mutations and the orange and blue bars at the ends of the reads indicate the sequencing adaptors.
Figure 3. Methods to rectify sequencing errors: (a) a scheme of how paired end reads with short sequencing reads allow the detection of sequencing errors; (b) schematic showing the concept of the duplex sequencing method; (c) how the consensus sequence is used to remove sequencing errors, adapted by permission from the Macmillan Publishers Ltd: Nature Protocols [61], copyright 2014. The black bar indicate the target inserts reads for sequencing and the orange and dark blue colored bars at end of the insert reads indicate sequencing adaptors; the yellow and light blue bars indicate the randomized duplex tags; and (d) Hypothetical mapping of the mutation frequency for variant library sequencing. The red line indicates the sequencing error rate of the ampicillin gene used as the cutoff. The asterisks indicate mutations and the orange and blue bars at the ends of the reads indicate the sequencing adaptors.
Ijms 16 23094 g003
Other aspects of the sequencing process that must be considered are the sequencing depth and library diversity. The sequencing depth can be calculated by considering the length of the target protein and the number of possible variants. For example, hyp-A is 300 bp long and the number of possible variants is 1900. By multiplying the protein length with the number of possible variants, we find that 570,000 is the number of base pairs required for the sequencing depth to be 1×. Thus, sequencing of two variant libraries with 100× depth requires 114 Mb, which is not a problem considering that the current benchtop sequencing platforms can produce gigabytes of data [63]. However, sequencing of very short sequences means that the sequence diversity would be very low, which is known to cause sequencing errors. For an Illumina platform, sequencing is performed by imaging the color produced by the clusters of the platform, and when the diversity is low, the color intensity cannot be properly detected, which causes errors in sequencing [64]. The easiest way to avoid sequencing errors is by increasing the spike-in percentage of the PhiX, which may complicate the calculation of the sequencing output. In addition, an extra sequence containing a series of random N sequences can be attached upstream of the insert reads to solve the low diversity problem. For now, it appears that the best sequencing pipeline for deep mutational scanning is MiSeq (Illumina), which offers numerous cartridges for different read lengths (50, 150 and 250 bp) with paired-end capability and contains a Real Time Analysis (RTA) software known to improve the data quality of low diversity samples [65].

2.5. Data Interpretation

As is the case for all applications of high-throughput sequencing data, the data analysis requires to be coupled to a logically and statistically acceptable workflow to answer the biological question or determine its relevance. For deep mutational scanning, the protein’s sequence and function relationship is the main question towards which the data analysis should be directed. Conceptually, this can be divided into two parts: analysis of the data to generate a mutational map or table and data interpretation of the mutational map. While data analysis involves bioinformatics to determine the abundance of each mutation at all positions, data interpretation is directed at determining the biological implication by calculating the abundance of mutations to different amino acids at different sites.
The nature of deep mutational scanning data is that it contains many mutations. Thus, data analysis is directed at identifying the mutations that result in a codon change in genes coding for different amino acids and determining the abundance of such mutations at all positions for each amino acid. The first part of deep mutational scanning data analysis does not differ from the conventional treatment of sequencing data. After removal of low-quality reads, the data should be mapped to the vector system sequence, the organism genome, and PhiX and only the unmapped reads should be collected. This step essentially removes any possible contaminations. If the sequencing library is constructed with the sequencing error precautions described in Figure 3a,b, the sequencing error can be removed by scripts called Enrich and by several scripts provided by Kennedy et al., respectively [61,66]. After the reads are mapped to the reference sequence, which in this case would be a protein, the mapped reads are extracted by SAMtools, and each read is translated into an amino acid [21,67]. The translated amino acids are then compared to the reference amino acid sequences to define mutations, and the abundance of these mutations is determined. As an extra precautionary measure, in a variant library that is expected to contain only single-site mutations, reads with more than one amino acid should be disregarded. Thus, the basic workflow of the deep mutational scanning data analysis is that the reads need to be translated into the corresponding amino acids first and then compared to the reference to identify and quantify mutations. The downstream analysis for data interpretation is then processed from this mutational table or map, which contains the abundance of each mutation (Figure 4a). It is also notable that there are tools that have been designed specifically for the analysis and visualization of deep mutational scanning data, which are Enrich and dms_tools [66,68].
Ijms 16 23094 g004 1024
Figure 4. (a) A hypothetical mutational map generated to show mutation frequency at each position; Part of the mutational map showing (b) extremely tolerant and critical residues to mutations; (c) tolerant to hydrophobic mutations and (d) tolerant to hydrophilic mutations. The x-axis indicates the protein residues and the y-axis indicates the possible amino acids. The color key represents the mutation frequency at each amino acid. The white color indicates that no mutation was found. The blue color indicates mutation frequency of loss-of-function variants and red color indicates mutation frequency of function-retained variants. The stop codon is indicated by *.
Figure 4. (a) A hypothetical mutational map generated to show mutation frequency at each position; Part of the mutational map showing (b) extremely tolerant and critical residues to mutations; (c) tolerant to hydrophobic mutations and (d) tolerant to hydrophilic mutations. The x-axis indicates the protein residues and the y-axis indicates the possible amino acids. The color key represents the mutation frequency at each amino acid. The white color indicates that no mutation was found. The blue color indicates mutation frequency of loss-of-function variants and red color indicates mutation frequency of function-retained variants. The stop codon is indicated by *.
Ijms 16 23094 g004
Interpretation of deep mutational sequencing data starts from the investigation of sites that have enriched mutations. With this perspective, it is advantageous to generate a colored heatmap from the data in the mutational table because it contains information that is easier to interpret by visualization (Figure 4a). First, the sites that are highly tolerant to mutations and those critical to mutations can be identified (Figure 4b). The sites with mutations indicated in blue color are enriched with every amino acid change, suggesting that these sites are critical to protein function, i.e., these sequences have a structurally or functionally important characteristic. Sites with mutations indicated in red color are bombarded with every amino acid change, suggesting that these sites are tolerant to mutations, i.e., these residues do not affect protein function or that their effects are negligible. This information is important because it identifies structurally important residues that are sequentially and structurally far from the active sites and difficult to predict. Sites with enriched mutations after selection are also quantified by different equations in different studies to show the importance of these residues. The characteristics of the residues can also be identified by investigating the type of amino acid changes. For example, the heatmap in Figure 4c,d is generated with the same dataset as that of the heatmap in Figure 4a (residues 7 to 9 and residues 12 to 13), but is generated with a different order of amino acid changes. From these hypothetical data, it can be inferred that the hydrophobic and hydrophilic characteristic of these residues greatly affects protein function (Figure 4c,d). Residues with a strong hydrophobic characteristic could be a site of dimer formation or a hydrophobic substrate-ligand pocket [21,69]. Such information can be used to define the protein stability in different environment as the function of the protein is closely related to the structural integrity of the protein. The results of the data interpretation are good examples of how the wealth of information generated from deep mutational scanning is used.

2.6. Limitations and Future Perspectives

The advance of high-throughput sequencing technology and screening assays has led to the development a new strategy of protein engineering called deep mutational scanning. This strategy has succeeded in showing the mutational consequence of different proteins by using different combinations of methods for variant protein generation, screening, and sequencing. However, mutational scanning is not applicable to all proteins because not all proteins are coupled with high-throughput screening assays to investigate their function. Furthermore, analysis of multi-site mutations and their epistasis is difficult because it requires a method that can screen tens of millions of mutants and be applicable only to small-sized proteins because of the sequencing read length.
The availability of a high-throughput screening system is a difficult concern to be addressed. For deep mutational scanning, the screening system requires the DNA to be linked to the phenotype and the screening system needs to be a high-throughput system. Designing such a screening system is difficult for enzymes and proteins for which the interacting substrate is not known, especially if the protein only functions in vivo. Nevertheless, fashioning a screening system is a challenging task but is essential to deep mutational scanning, which is its biggest limitation. A means to overcome this problem would greatly enhance the future prospects of deep mutational scanning.
Analysis of multi-site mutations and their epistasis has been performed in a few studies [10,22,33,48,70]. The results from deep mutational scanning on a WW domain have been used to produce a model for mutational epistasis, which showed a good prediction rate of 70% [22]. Additionally, epistasis in a few mutation sites has been tested for differential ligand binding of PSD95pdz3 [10]. Interestingly, some mutation sites that had a deleterious effect on single-site mutagenesis showed significant epistasis to enhance binding to a different ligand [10]. For large-sized proteins, the sequencing length is insufficient to detect the effect of multi-site mutation in one read. To resolve this issue, tag-directed assembly methods can be applied [71,72,73,74]. These methods require specific barcode for each protein variants to assemble the whole protein sequence with same barcoded reads. However, unlike the random barcode system used to reduce sequencing errors this method require specific barcodes for millions protein variants. This is cost inefficient, as the sequencing library of each protein variant needs to be constructed separately. Another approach to address this problem is to analyze the single-site mutagenesis library first and select a few domains to perform another single-site mutagenesis on the selected mutants. This would require a barcode for the second single-site mutagenesis to distinguish these mutations from those of the first and second round of mutagenesis, and the sequence space will not be highly covered. In addition, the distance between the mutated domains in the first and second rounds needs to be less than the sequencing read length. Although many limitations exist, deep mutational scanning application is the most promising method to study epistasis in different mutation sites.
Overall, deep mutational scanning strategies provide mutational data of a protein by linking mutation to the phenotype. Such mutational data serve as valuable information to understand the proteins’ functional mechanism, which can then be applied to clinical application of antibodies, studies on human genetic diseases, structural protein sciences, and protein engineering [14,23,27,53]. In addition, further developments of this approach are promising because it is based on two rapidly evolving technologies: high-throughput screening and high-throughput sequencing. In this review, we have introduced different methodologies and concepts that can be applied to deep mutational scanning based on recent studies. Although some questions remain regarding the mutational epistasis of large-sized protein targets, deep mutation scanning provides answers to many questions that involve the protein’s sequence and function relationship.

Acknowledgments

The authors acknowledge financial support from the Basic Core Technology Development Program for the Oceans and the Polar Regions (2011-0021053) through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (MISP).

Author Contributions

All the authors were involved in the discussion and in writing the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dalbadie-McFarland, G.; Cohen, L.W.; Riggs, A.D.; Morin, C.; Itakura, K.; Richards, J.H. Oligonucleotide-directed mutagenesis as a general and powerful method for studies of protein function. Proc. Natl. Acad. Sci. USA 1982, 79, 6409–6413. [Google Scholar] [CrossRef] [PubMed]
  2. Blundell, T.L. Problems and solutions in protein engineering—Towards rational design. Trends Biotechnol. 1994, 12, 145–148. [Google Scholar] [CrossRef]
  3. McManus, S.; Riechmann, L. Use of 2D NMR, protein engineering, and molecular modeling to study the hapten-binding site of an antibody Fv fragment against 2-phenyloxazolone. Biochemistry 1991, 30, 5851–5857. [Google Scholar] [CrossRef] [PubMed]
  4. Hakoshima, T.; Toda, S.; Sugio, S.; Tomita, K.; Nishikawa, S.; Morioka, H.; Fuchimura, K.; Kimura, T.; Uesugi, S.; Ohtsuka, E.; et al. Conformational properties of the guanine-binding site of ribonuclease T1 inferred from the X-ray structure and protein engineering. Protein Eng. 1988, 2, 55–61. [Google Scholar] [CrossRef] [PubMed]
  5. Rohl, C.A.; Strauss, C.E.; Misura, K.M.; Baker, D. Protein structure prediction using rosetta. Methods Enzymol. 2004, 383, 66–93. [Google Scholar] [PubMed]
  6. Liu, Y.; Kuhlman, B. Rosettadesign server for protein design. Nucleic Acids Res. 2006, 34, W235–W238. [Google Scholar] [CrossRef] [PubMed]
  7. Lutz, S. Beyond directed evolution—Semi-rational protein engineering and design. Curr. Opin. Biotechnol. 2010, 21, 734–743. [Google Scholar] [CrossRef] [PubMed]
  8. Gombault, A.; Godin, F.; Sy, D.; Legrand, B.; Chautard, H.; Vallee, B.; Vovelle, F.; Benedetti, H. Molecular basis of the Tfs1/Ira2 interaction: A combined protein engineering and molecular modelling study. J. Mol. Biol. 2007, 374, 604–617. [Google Scholar] [CrossRef] [PubMed]
  9. Engvall, E.; Perlmann, P. Enzyme-linked immunosorbent assay (ELISA). Quantitative assay of immunoglobulin G. Immunochemistry 1971, 8, 871–874. [Google Scholar] [CrossRef]
  10. McLaughlin, R.N., Jr.; Poelwijk, F.J.; Raman, A.; Gosal, W.S.; Ranganathan, R. The spatial architecture of protein function and adaptation. Nature 2012, 491, 138–142. [Google Scholar] [CrossRef] [PubMed]
  11. Araya, C.L.; Fowler, D.M. Deep mutational scanning: Assessing protein function on a massive scale. Trends Biotechnol. 2011, 29, 435–442. [Google Scholar] [CrossRef] [PubMed]
  12. Smith, J.M. Natural selection and the concept of a protein space. Nature 1970, 225, 563–564. [Google Scholar] [CrossRef] [PubMed]
  13. Freeman, A.M.; Mole, B.M.; Silversmith, R.E.; Bourret, R.B. Action at a distance: Amino acid substitutions that affect binding of the phosphorylated chey response regulator and catalysis of dephosphorylation can be far from the chez phosphatase active site. J. Bacteriol. 2011, 193, 4709–4718. [Google Scholar] [CrossRef] [PubMed]
  14. Fowler, D.M.; Fields, S. Deep mutational scanning: A new style of protein science. Nat. Methods 2014, 11, 801–807. [Google Scholar] [CrossRef]
  15. O’Neil, K.T.; Hoess, R.H. Phage display: Protein engineering by directed evolution. Curr. Opin. Struct. Biol. 1995, 5, 443–449. [Google Scholar] [CrossRef]
  16. Fernandez-Gacio, A.; Uguen, M.; Fastrez, J. Phage display as a tool for the directed evolution of enzymes. Trends Biotechnol. 2003, 21, 408–414. [Google Scholar] [CrossRef]
  17. Chen, Y.; Wiesmann, C.; Fuh, G.; Li, B.; Christinger, H.W.; McKay, P.; de Vos, A.M.; Lowman, H.B. Selection and analysis of an optimized anti-VEGF antibody: Crystal structure of an affinity-matured fab in complex with antigen. J. Mol. Biol. 1999, 293, 865–881. [Google Scholar] [CrossRef] [PubMed]
  18. Hibbert, E.G.; Dalby, P.A. Directed evolution strategies for improved enzymatic performance. Microb. Cell Factories 2005, 4. [Google Scholar] [CrossRef] [PubMed]
  19. Shenoy, S.R.; Jayaram, B. Proteins: Sequence to structure and function—Current status. Curr. Protein Pept. Sci. 2010, 11, 498–514. [Google Scholar] [CrossRef] [PubMed]
  20. Sadowski, M.I.; Jones, D.T. The sequence-structure relationship and protein function prediction. Curr. Opin. Struct. Biol. 2009, 19, 357–362. [Google Scholar] [CrossRef] [PubMed]
  21. Shin, H.; Cho, Y.; Choe, D.H.; Jeong, Y.; Cho, S.; Kim, S.C.; Cho, B.K. Exploring the functional residues in a flavin-binding fluorescent protein using deep mutational scanning. PLoS ONE 2014, 9, e97817. [Google Scholar] [CrossRef] [PubMed]
  22. Fowler, D.M.; Araya, C.L.; Fleishman, S.J.; Kellogg, E.H.; Stephany, J.J.; Baker, D.; Fields, S. High-resolution mapping of protein sequence-function relationships. Nat. Methods 2010, 7, 741–746. [Google Scholar] [CrossRef] [PubMed]
  23. Romero, P.A.; Tran, T.M.; Abate, A.R. Dissecting enzyme function with microfluidic-based deep mutational scanning. Proc. Natl. Acad. Sci. USA 2015, 112, 7159–7164. [Google Scholar] [CrossRef] [PubMed]
  24. Whitehead, T.A.; Chevalier, A.; Song, Y.; Dreyfus, C.; Fleishman, S.J.; de Mattos, C.; Myers, C.A.; Kamisetty, H.; Blair, P.; Wilson, I.A.; et al. Optimization of affinity, specificity and function of designed influenza inhibitors using deep sequencing. Nat. Biotechnol. 2012, 30, 543–548. [Google Scholar] [CrossRef] [PubMed]
  25. Zheng, L.; Baumann, U.; Reymond, J.L. An efficient one-step site-directed and site-saturation mutagenesis protocol. Nucleic Acids Res. 2004, 32. [Google Scholar] [CrossRef] [PubMed]
  26. Adkar, B.V.; Tripathi, A.; Sahoo, A.; Bajaj, K.; Goswami, D.; Chakrabarti, P.; Swarnkar, M.K.; Gokhale, R.S.; Varadarajan, R. Protein model discrimination using mutational sensitivity derived from deep sequencing. Structure 2012, 20, 371–381. [Google Scholar] [CrossRef] [PubMed]
  27. Starita, L.M.; Young, D.L.; Islam, M.; Kitzman, J.O.; Gullingsrud, J.; Hause, R.J.; Fowler, D.M.; Parvin, J.D.; Shendure, J.; Fields, S. Massively parallel functional analysis of brca1 ring domain variants. Genetics 2015, 200, 413–422. [Google Scholar] [CrossRef] [PubMed]
  28. Melnikov, A.; Rogov, P.; Wang, L.; Gnirke, A.; Mikkelsen, T.S. Comprehensive mutational scanning of a kinase in vivo reveals substrate-dependent fitness landscapes. Nucleic Acids Res. 2014, 42. [Google Scholar] [CrossRef] [PubMed]
  29. Melamed, D.; Young, D.L.; Gamble, C.E.; Miller, C.R.; Fields, S. Deep mutational scanning of an RRM domain of the Saccharomyces cerevisiae poly(A)-binding protein. RNA 2013, 19, 1537–1551. [Google Scholar] [CrossRef] [PubMed]
  30. Forsyth, C.M.; Juan, V.; Akamatsu, Y.; DuBridge, R.B.; Doan, M.; Ivanov, A.V.; Ma, Z.; Polakoff, D.; Razo, J.; Wilson, K.; et al. Deep mutational scanning of an antibody against epidermal growth factor receptor using mammalian cell display and massively parallel pyrosequencing. MAbs 2013, 5, 523–532. [Google Scholar] [CrossRef] [PubMed]
  31. Bloom, J.D. An experimentally determined evolutionary model dramatically improves phylogenetic fit. Mol. Biol. Evol. 2014, 31, 1956–1978. [Google Scholar] [CrossRef] [PubMed]
  32. Firnberg, E.; Labonte, J.W.; Gray, J.J.; Ostermeier, M. A comprehensive, high-resolution map of a gene’s fitness landscape. Mol. Biol. Evol. 2014, 31, 1581–1592. [Google Scholar] [CrossRef] [PubMed]
  33. Olson, C.A.; Wu, N.C.; Sun, R. A comprehensive biophysical description of pairwise epistasis throughout an entire protein domain. Curr. Biol. 2014, 24, 2643–2651. [Google Scholar] [CrossRef] [PubMed]
  34. Qi, H.; Olson, C.A.; Wu, N.C.; Ke, R.; Loverdo, C.; Chu, V.; Truong, S.; Remenyi, R.; Chen, Z.; Du, Y.; et al. A quantitative high-resolution genetic profile rapidly identifies sequence determinants of hepatitis C viral fitness and drug sensitivity. PLoS Pathog. 2014, 10, e1004064. [Google Scholar] [CrossRef] [PubMed]
  35. Thyagarajan, B.; Bloom, J.D. The inherent mutational tolerance and antigenic evolvability of influenza hemagglutinin. Elife 2014, 3. [Google Scholar] [CrossRef] [PubMed]
  36. Stiffler, M.A.; Hekstra, D.R.; Ranganathan, R. Evolvability as a function of purifying selection in TEM-1 β-lactamase. Cell 2015, 160, 882–892. [Google Scholar] [CrossRef] [PubMed]
  37. Fowler, D.M.; Stephany, J.J.; Fields, S. Measuring the activity of protein variants on a large scale using deep mutational scanning. Nat. Protoc. 2014, 9, 2267–2284. [Google Scholar] [CrossRef] [PubMed]
  38. Al-Mawsawi, L.Q.; Wu, N.C.; Olson, C.A.; Shi, V.C.; Qi, H.; Zheng, X.; Wu, T.T.; Sun, R. High-throughput profiling of point mutations across the HIV-1 genome. Retrovirology 2014, 11. [Google Scholar] [CrossRef] [PubMed]
  39. Wu, N.C.; Young, A.P.; Al-Mawsawi, L.Q.; Olson, C.A.; Feng, J.; Qi, H.; Chen, S.H.; Lu, I.H.; Lin, C.Y.; Chin, R.G.; et al. High-throughput profiling of influenza a virus hemagglutinin gene at single-nucleotide resolution. Sci. Rep. 2014, 4. [Google Scholar] [CrossRef] [PubMed]
  40. Wu, N.C.; Young, A.P.; Al-Mawsawi, L.Q.; Olson, C.A.; Feng, J.; Qi, H.; Luan, H.H.; Li, X.; Wu, T.T.; Sun, R. High-throughput identification of loss-of-function mutations for anti-interferon activity in the influenza a virus ns segment. J. Virol. 2014, 88, 10157–10164. [Google Scholar] [CrossRef] [PubMed]
  41. Wu, N.C.; Olson, C.A.; Du, Y.; Le, S.; Tran, K.; Remenyi, R.; Gong, D.; Al-Mawsawi, L.Q.; Qi, H.; Wu, T.T.; et al. Functional constraint profiling of a viral protein reveals discordance of evolutionary conservation and functionality. PLoS Genet. 2015, 11, e1005310. [Google Scholar] [CrossRef] [PubMed]
  42. McCullum, E.O.; Williams, B.A.; Zhang, J.; Chaput, J.C. Random mutagenesis by error-prone PCR. Methods Mol. Biol. 2010, 634, 103–109. [Google Scholar] [PubMed]
  43. Currin, A.; Swainston, N.; Day, P.J.; Kell, D.B. Synthetic biology for the directed evolution of protein biocatalysts: Navigating sequence space intelligently. Chem. Soc. Rev. 2015, 44, 1172–1239. [Google Scholar] [CrossRef] [PubMed]
  44. Romero, P.A.; Arnold, F.H. Exploring protein fitness landscapes by directed evolution. Nat. Rev. Mol. Cell Biol. 2009, 10, 866–876. [Google Scholar] [CrossRef] [PubMed]
  45. Clackson, T.; Hoogenboom, H.R.; Griffiths, A.D.; Winter, G. Making antibody fragments using phage display libraries. Nature 1991, 352, 624–628. [Google Scholar] [CrossRef] [PubMed]
  46. Dai, M.; Temirov, J.; Pesavento, E.; Kiss, C.; Velappan, N.; Pavlik, P.; Werner, J.H.; Bradbury, A.R. Using T7 phage display to select GFP-based binders. Protein Eng. Des. Sel. 2008, 21, 413–424. [Google Scholar] [CrossRef] [PubMed]
  47. Fujino, Y.; Fujita, R.; Wada, K.; Fujishige, K.; Kanamori, T.; Hunt, L.; Shimizu, Y.; Ueda, T. Robust in vitro affinity maturation strategy based on interface-focused high-throughput mutational scanning. Biochem. Biophys. Res. Commun. 2012, 428, 395–400. [Google Scholar] [CrossRef] [PubMed]
  48. Araya, C.L.; Fowler, D.M.; Chen, W.; Muniez, I.; Kelly, J.W.; Fields, S. A fundamental protein property, thermodynamic stability, revealed solely from large-scale measurements of protein function. Proc. Natl. Acad. Sci. USA 2012, 109, 16858–16863. [Google Scholar] [CrossRef] [PubMed]
  49. Pellis, M.; Muyldermans, S.; Vincke, C. Bacterial two hybrid: A versatile one-step intracellular selection method. Methods Mol. Biol. 2012, 911, 135–150. [Google Scholar] [PubMed]
  50. Evangelista, C.; Lockshon, D.; Fields, S. The yeast two-hybrid system: Prospects for protein linkage maps. Trends Cell Biol. 1996, 6, 196–199. [Google Scholar] [CrossRef]
  51. Monti-Bragadin, C.; Babudri, N.; Samer, L. Expression of the plasmid pKM101—Determined DNA repair system in recA- and lex- strains of Escherichia coli. Mol. Gen. Genet. 1976, 145, 303–306. [Google Scholar] [CrossRef] [PubMed]
  52. Zhou, X.; Ren, L.; Meng, Q.; Li, Y.; Yu, Y.; Yu, J. The next-generation sequencing technology and application. Protein Cell 2010, 1, 520–536. [Google Scholar] [CrossRef] [PubMed]
  53. Tsuchihara, K. Clinical application of next-generation sequencing technologies to achieve cancer precision medicine. Gan To Kagaku Ryoho Cancer Chemother. 2014, 41, 1–6. [Google Scholar]
  54. Shin, H.; Hong, S.J.; Kim, H.; Yoo, C.; Lee, H.; Choi, H.K.; Lee, C.G.; Cho, B.K. Elucidation of the growth delimitation of Dunaliella tertiolecta under nitrogen stress by integrating transcriptome and peptidome analysis. Bioresour. Technol. 2015, 194, 57–66. [Google Scholar] [CrossRef] [PubMed]
  55. Quail, M.A.; Kozarewa, I.; Smith, F.; Scally, A.; Stephens, P.J.; Durbin, R.; Swerdlow, H.; Turner, D.J. A large genome center’s improvements to the illumina sequencing system. Nat. Methods 2008, 5, 1005–1010. [Google Scholar] [CrossRef] [PubMed]
  56. Lou, D.I.; Hussmann, J.A.; McBee, R.M.; Acevedo, A.; Andino, R.; Press, W.H.; Sawyer, S.L. High-throughput DNA sequencing errors are reduced by orders of magnitude using circle sequencing. Proc. Natl. Acad. Sci. USA 2013, 110, 19872–19877. [Google Scholar] [CrossRef] [PubMed]
  57. Dean, F.B.; Nelson, J.R.; Giesler, T.L.; Lasken, R.S. Rapid amplification of plasmid and phage DNA using Phi29 DNA polymerase and multiply-primed rolling circle amplification. Genome Res. 2001, 11, 1095–1099. [Google Scholar] [CrossRef] [PubMed]
  58. Jabara, C.B.; Jones, C.D.; Roach, J.; Anderson, J.A.; Swanstrom, R. Accurate sampling and deep sequencing of the HIV-1 protease gene using a Primer ID. Proc. Natl. Acad. Sci. USA 2011, 108, 20166–20171. [Google Scholar] [CrossRef] [PubMed]
  59. Kinde, I.; Wu, J.; Papadopoulos, N.; Kinzler, K.W.; Vogelstein, B. Detection and quantification of rare mutations with massively parallel sequencing. Proc. Natl. Acad. Sci. USA 2011, 108, 9530–9535. [Google Scholar] [CrossRef] [PubMed]
  60. Gout, J.F.; Thomas, W.K.; Smith, Z.; Okamoto, K.; Lynch, M. Large-scale detection of in vivo transcription errors. Proc. Natl. Acad. Sci. USA 2013, 110, 18584–18589. [Google Scholar] [CrossRef] [PubMed]
  61. Kennedy, S.R.; Schmitt, M.W.; Fox, E.J.; Kohrn, B.F.; Salk, J.J.; Ahn, E.H.; Prindle, M.J.; Kuong, K.J.; Shen, J.C.; Risques, R.A.; et al. Detecting ultralow-frequency mutations by duplex sequencing. Nat. Protoc. 2014, 9, 2586–2606. [Google Scholar] [CrossRef] [PubMed]
  62. Schmitt, M.W.; Kennedy, S.R.; Salk, J.J.; Fox, E.J.; Hiatt, J.B.; Loeb, L.A. Detection of ultra-rare mutations by next-generation sequencing. Proc. Natl. Acad. Sci. USA 2012, 109, 14508–14513. [Google Scholar] [CrossRef] [PubMed]
  63. Loman, N.J.; Misra, R.V.; Dallman, T.J.; Constantinidou, C.; Gharbia, S.E.; Wain, J.; Pallen, M.J. Performance comparison of benchtop high-throughput sequencing platforms. Nat. Biotechnol. 2012, 30, 434–439. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Mitra, A.; Skrzypczak, M.; Ginalski, K.; Rowicka, M. Strategies for achieving high sequencing accuracy for low diversity samples and avoiding sample bleeding using illumina platform. PLoS ONE 2015, 10, e0120520. [Google Scholar] [CrossRef] [PubMed]
  65. Nelson, M.C.; Morrison, H.G.; Benjamino, J.; Grim, S.L.; Graf, J. Analysis, optimization and verification of illumina-generated 16S rRNA gene amplicon surveys. PLoS ONE 2014, 9, e94249. [Google Scholar] [CrossRef] [PubMed]
  66. Fowler, D.M.; Araya, C.L.; Gerard, W.; Fields, S. Enrich: Software for analysis of protein function by enrichment and depletion of variants. Bioinformatics 2011, 27, 3430–3431. [Google Scholar] [CrossRef] [PubMed]
  67. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; 1000 Genome Project Data Processing Subgroup. The sequence alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  68. Bloom, J.D. Software for the analysis and visualization of deep mutational scanning data. BMC Bioinform. 2015, 16. [Google Scholar] [CrossRef] [PubMed]
  69. Kawaguchi, S.; Kuramitsu, S. Thermodynamics and molecular simulation analysis of hydrophobic substrate recognition by aminotransferases. J. Biol. Chem. 1998, 273, 18353–18364. [Google Scholar] [CrossRef] [PubMed]
  70. Bank, C.; Hietpas, R.T.; Jensen, J.D.; Bolon, D.N. A systematic survey of an intragenic epistatic landscape. Mol. Biol. Evol. 2015, 32, 229–238. [Google Scholar] [CrossRef] [PubMed]
  71. Hiatt, J.B.; Patwardhan, R.P.; Turner, E.H.; Lee, C.; Shendure, J. Parallel, tag-directed assembly of locally derived short sequence reads. Nat. Methods 2010, 7, 119–122. [Google Scholar] [CrossRef] [PubMed]
  72. Hong, L.Z.; Hong, S.; Wong, H.T.; Aw, P.P.; Cheng, Y.; Wilm, A.; de Sessions, P.F.; Lim, S.G.; Nagarajan, N.; Hibberd, M.L.; et al. BAsE-Seq: A method for obtaining long viral haplotypes from short sequence reads. Genome Biol. 2014, 15. [Google Scholar] [CrossRef]
  73. Wu, N.C.; de La Cruz, J.; Al-Mawsawi, L.Q.; Olson, C.A.; Qi, H.; Luan, H.H.; Nguyen, N.; Du, Y.; Le, S.; Wu, T.T.; et al. HIV-1 quasispecies delineation by tag linkage deep sequencing. PLoS ONE 2014, 9, e97505. [Google Scholar] [CrossRef] [PubMed]
  74. Borgstrom, E.; Redin, D.; Lundin, S.; Berglund, E.; Andersson, A.F.; Ahmadian, A. Phasing of single DNA molecules by massively parallel barcoding. Nat. Commun. 2015, 6. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Shin, H.; Cho, B.-K. Rational Protein Engineering Guided by Deep Mutational Scanning. Int. J. Mol. Sci. 2015, 16, 23094-23110. https://doi.org/10.3390/ijms160923094

AMA Style

Shin H, Cho B-K. Rational Protein Engineering Guided by Deep Mutational Scanning. International Journal of Molecular Sciences. 2015; 16(9):23094-23110. https://doi.org/10.3390/ijms160923094

Chicago/Turabian Style

Shin, HyeonSeok, and Byung-Kwan Cho. 2015. "Rational Protein Engineering Guided by Deep Mutational Scanning" International Journal of Molecular Sciences 16, no. 9: 23094-23110. https://doi.org/10.3390/ijms160923094

APA Style

Shin, H., & Cho, B.-K. (2015). Rational Protein Engineering Guided by Deep Mutational Scanning. International Journal of Molecular Sciences, 16(9), 23094-23110. https://doi.org/10.3390/ijms160923094

Article Metrics

Back to TopTop