Orthogonal combinatorial mutagenesis: a codon-level combinatorial mutagenesis method useful for low multiplicity and amino acid-scanning protocols

We describe here a method to generate combinatorial libraries of oligonucleotides mutated at the codon-level, with control of the mutagenesis rate so as to create predictable binomial distributions of mutants. The method allows enrichment of the libraries with single, double or larger multiplicity of amino acid replacements by appropriate choice of the mutagenesis rate, depending on the concentration of synthetic precursors. The method makes use of two sets of deoxynucleoside-phosphoramidites bearing orthogonal protecting groups [4,4′-dimethoxytrityl (DMT) and 9-fluorenylmethoxycarbonyl (Fmoc)] in the 5′ hydroxyl. These phosphoramidites are divergently combined during automated synthesis in such a way that wild-type codons are assembled with commercial DMT-deoxynucleoside-methyl-phosphoramidites while mutant codons are assembled with Fmoc-deoxynucleoside-methyl-phosphoramidites in an NNG/C fashion in a single synthesis column. This method is easily automated and suitable for low mutagenesis rates and large windows, such as those required for directed evolution and alanine scanning. Through the assembly of three oligonucleotide libraries at different mutagenesis rates, followed by cloning at the polylinker region of plasmid pUC18 and sequencing of 129 clones, we concluded that the method performs essentially as intended.     

A computer program for the estimation of protein and nucleic acid sequence diversity in random point mutagenesis libraries

A computer program for the generation and analysis of in silico random point mutagenesis libraries is described. The program operates by mutagenizing an input nucleic acid sequence according to mutation parameters specified by the user for each sequence position and type of point mutation. The program can mimic almost any type of random mutagenesis library, including those produced via error-prone PCR (ep-PCR), mutator Escherichia coli strains, chemical mutagenesis, and doped or random oligonucleotide synthesis. The program analyzes the generated nucleic acid sequences and/or the associated protein library to produce several estimates of library diversity (number of unique sequences, point mutations, and single point mutants) and the rate of saturation of these diversities during experimental screening or selection of clones. This information allows one to select the optimal screen size for a given mutagenesis library, necessary to efficiently obtain a certain coverage of the sequence-space. The program also reports the abundance of each specific protein mutation at each sequence position, which is useful as a measure of the level and type of mutation bias in the library. Alternatively, one can use the program to evaluate the relative merits of preexisting libraries, or to examine various hypothetical mutation schemes to determine the optimal method for creating a library that serves the screen/selection of interest. Simulated libraries of at least 10⁹ sequences are accessible by the numerical algorithm with currently available personal computers; an analytical algorithm is also available which can rapidly calculate a subset of the numerical statistics in libraries of arbitrarily large size. A multi-type double-strand stochastic model of ep-PCR is developed in an appendix to demonstrate the applicability of the algorithm to amplifying mutagenesis procedures. Estimators of DNA polymerase mutation-type-specific error rates are derived using the model. Analyses of an alpha-synuclein ep-PCR library and NNS synthetic oligonucleotide libraries are given as examples.     

Oligonucleotide site-directed mutagenesis in Xenopus egg extracts.

Addition of M13mp18 single-stranded DNA annealed with an oligonucleotide to a Xenopus egg extract results in a rapid and efficient incorporation of the oligonucleotide in a complete double-stranded supercoiled molecule. Both the efficiency of DNA synthesis and the recovery of complete double-stranded molecules are increased relative to the reaction carried out by the classical technique using the E. coli Klenow DNA polymerase, DNA ligase, dNTPs, ATP and ions. Site specific mutagenesis was assayed by reverting a point mutation in the lacz region of M13mp18. The color assay described by Messing and sequencing of the DNA extracted from isolated plaques was used to check for the reversion. A 2 hr incubation of the heteroduplex carrying the mutagenic oligonucleotide in the Klenow-ligase-dNTP mixture allows a recovery of 6% mutant phage after transformation of competent cells with the reaction products. Using the Xenopus egg extract, 83% mutant phage were recovered after the same incubation time, in reactions entirely performed in parallel. The Xenopus extract is stable and contains all components required for the assay, including all ionic and protein factors; thus the only addition is the annealed DNA. Such an eukaryotic system is therefore an attractive alternative to the reconstituted prokaryotic DNA polymerase-DNA ligase system for site specific mutagenesis.     

Directed evolution of green fluorescent protein by a new versatile PCR strategy for site-directed and semi-random mutagenesis

To develop a simple, speedy, economical and widely applicable method for multiple-site mutagenesis, we have substantially modified the Quik-Change™ Site-Directed Mutagenesis Kit protocol (Stratagene, La Jolla, CA). Our new protocol consists of (i) a PCR reaction using an in vitro technique, LDA (ligation-during-amplification), (ii) a DpnI treatment to digest parental DNA and to make megaprimers and (iii) a synthesis of double-stranded plasmid DNA for bacterial transformation. While the Quik Change™ Kit protocol introduces mutations at a single site, requiring two complementary mutagenic oligonucleotides, our new protocol requires only one mutagenic oligonucleotide for a mutation site, and can introduce mutations in a plasmid at multiple sites simultaneously. A targeting efficiency >70% was consistently achieved for multiple-site mutagenesis. Furthermore, the new protocol allows random mutagenesis with degenerative primers, because it does not use two complementary primers. Our mutagenesis strategy was successfully used to alter the fluorescence properties of green fluorescent protein (GFP), creating a new-color GFP mutant, cyan-green fluorescent protein (CGFP). An eminent feature of CGFP is its remarkable stability in a wide pH range (pH 4–12). The use of CGFP would allow us to monitor protein localization quantitatively in acidic organelles in secretory pathways.     

The Fidelity of Template-Directed Oligonucleotide Ligation and the Inevitability of Polymerase Function

The first living systems may have employed template-directed oligonucleotide ligation for replication. The utility of oligonucleotide ligation as a mechanism for the origin and evolution of life is in part dependent on its fidelity. We have devised a method for evaluating ligation fidelity in which ligation substrates are selected from random sequence libraries. The fidelities of chemical and enzymatic ligation are compared under a variety of conditions. While reaction conditions can be found that promote high fidelity copying, departure from these conditions leads to error-prone copying. In particular, ligation reactions with shorter oligonucleotide substrates are less efficient but more faithful. These results support a model for origins in which there was selective pressure for template-directed oligonucleotide ligation to be gradually supplanted by mononucleotide polymerization.     

Thermostable DNA Ligase-Mediated PCR Production of Circular Plasmid (PPCP) and Its Application in Directed Evolution via In situ Error-Prone PCR

Polymerase chain reaction (PCR) is a powerful method to produce linear DNA fragments. Here we describe the Tma thermostable DNA ligase-mediated PCR production of circular plasmid (PPCP) and its application in directed evolution via in situ error-prone PCR. In this thermostable DNA ligase-mediated whole-plasmid amplification method, the resultant DNA nick between the 5′ end of the PCR primer and the extended newly synthesized DNA 3′ end of each PCR cycle is ligated by Tma DNA ligase, resulting in circular plasmid DNA product that can be directly transformed. The template plasmid DNA is eliminated by ‘selection marker swapping’ upon transformation. When performed under an error-prone condition with Taq DNA polymerase, PPCP allows one-step construction of mutagenesis libraries based on in situ error-prone PCR so that random mutations are introduced into the target gene without altering the expression vector plasmid. A significant difference between PPCP and previously published methods is that PPCP allows exponential amplification of circular DNA. We used this method to create random mutagenesis libraries of a xylanase gene and two cellulase genes. Screening of these libraries resulted in mutant proteins with desired properties, demonstrating the usefulness of in situ error-prone PPCP for creating random mutagenesis libraries for directed evolution.     

Novel ceftazidime-resistance β-lactamases generated by a codon-based mutagenesis method and selection

Four known and nine new ceftazidime-resistance β-lactamases were generated by a novel, contaminating codon-based mutagenesis approach. In this method, wild-type codons are spiked with a set of mutant codons during oligonucleotide synthesis, generating random combinatorial libraries of primers that contain few codon replacements per variant. Mutant codons are assembled by tandem addition of a diluted mixture of five Fmoc-dimer amidites to the growing oligo and a mixture of four DMTr-monomer amidites to generate 20 trinucleotides that encode a set of 18 amino acids. Wild-type codons are assembled with conventional chemistry and the whole process takes place in only one synthesis column, making its automation feasible. The random and binomial behavior of this approach was tested in the polylinker region of plasmid pUC19 by the synthesis of three oligonucleotide libraries mutagenized at different rates and cloned as mutagenic cassettes. Additionally, the method was biologically assessed by mutating six contiguous codons that encode amino acids 237–243 (ABL numbering) of the TEM_pUC19 β-lactamase, which is functionally equivalent to the clinically important TEM-1 β-lactamase. The best ceftazidime-recognizing variant was a triple mutant, R164H:E240K: R241A, displaying a 333-fold higher resistance than the wild-type enzyme.     

Herpesviruses encode an unusual protein-serine/threonine kinase which is nonessential for growth in cultured cells.

We have performed large-scale random oligonucleotide insertion mutagenesis on a 41-kbp genomic segment derived from the unique long (UL) region of the alphaherpesvirus pseudorabies virus (PRV). This procedure has resulted in the generation of a series of PRV strains, each carrying a single gene whose termination of translation is induced by the inserted oligonucleotide. To relate the genes that were involved in the mutagenization to genes previously identified in herpes simplex virus type 1, the prototype alphaherpesvirus, we have performed cross-hybridization studies. In this way, we have mapped the location of the homolog of a gene which was described to have sequence characteristics of a eukaryotic phosphotransferase. We characterized the phenotype of a mutant PRV strain lacking this putative phosphotransferase also the phenotype of a PRV strain lacking, in addition to the UL-encoded putative phosphotransferase, the protein kinase encoded within the unique short region of the virus. To assess the enzymatic activity of the UL region-encoded phosphotransferase, we expressed the gene transiently in a eukaryotic expression system. Immunoprecipitation of the protein followed by kinase assays and phosphoamino acid analyses revealed protein-serine/threonine kinase activity. Implications of sequence divergence of this protein from classical protein-serine/threonine kinases for kinase structure and function are discussed in view of the recent resolution of the structure of the catalytic domain of cyclic AMP-dependent protein kinase.     

Application of Mu in vitro transposition for high-precision mapping of protein–protein interfaces on a yeast two-hybrid platform

High-precision mapping of regions involved in protein–protein interfaces of interacting protein partners is an essential component on a path to understand various cellular functions. Transposon-based systems, particularly those involving in vitro reactions, offer exhaustive insertion mutant libraries and high-throughput platforms for many types of genetic analyses. We present here a genetic strategy to accurately map interacting protein regions at amino acid precision that is based on transposition-assisted construction, sampling, and analysis of a comprehensive insertion mutant library. The methodology integrates random pentapeptide mutagenesis of proteins, yeast two-hybrid screening, and high-resolution genetic footprinting. This straightforward strategy is general, and it provides a rapid and easy means to identify critical contact regions in proteins without the requirement of prior structural knowledge.     

Sequence dependence of cross-hybridization on short oligo microarrays

One of the critical problems in the short oligo microarray technology is how to deal with cross-hybridization that produces spurious data. Little is known about the details of cross-hybridization effect at molecular level. Here, we report a free energy analysis of cross-hybridization on short oligo microarrays using data from a spike-in study. Our analysis revealed that cross-hybridization on the arrays is mostly caused by oligo fragments with a run of 10–16 nt complementary to the probes. Mismatches were estimated to be energetically much more costly in cross-hybridization than that in gene-specific hybridization, implying that the sources of cross-hybridization must be very different between a PM–MM probe pair. Consequently, it is unreliable to use MM probe signal to track cross-hybridizing signal on a corresponding PM probe. Our results also showed that the oligo fragments tend to bind to the 5′ ends of the probes, and are rarely seen at the 3′ ends. These results are useful for microarray design and data analysis.     

Reducing mutational bias in random protein libraries

The success of protein optimization through directed molecular evolution depends to a large extent on the size and quality of the displayed library. Current low-fidelity DNA polymerases that are commonly used during random mutagenesis and recombination in vitro display strong mutational preferences, favoring the substitution of certain nucleotides over others. The result is a biased and reduced functional diversity in the library under selection. In an effort to reduce mutational bias, we combined two different low-fidelity DNA polymerases, Taq and Mutazyme, which have opposite mutational spectra. As a first step, random mutants of the Bacillus thuringiensis cry9Ca1 gene were generated by separate error-prone polymerase chain reactions (PCRs) with each of the two polymerases. Subsequent shuffling by staggered extension process (StEP) of the PCR products resulted in intermediate numbers of AT and GC substitutions, compared to the Taq or Mutazyme error-prone PCR libraries. This strategy should allow generating unbiased libraries or libraries with a specific degree of mutational bias by applying optimal mutagenesis frequencies during error-prone PCR and controlling the concentration of template in the shuffling reaction while taking into account the GC content of the target gene.     

Human and mouse oligonucleotide-based array CGH

Array-based comparative genomic hybridization is a high resolution method for measuring chromosomal copy number changes. Here we present a validated protocol using in-house spotted oligonucleotide libraries for array comparative genomic hybridization (CGH). This oligo array CGH platform yields reproducible results and is capable of detecting single copy gains, multi-copy amplifications as well as homozygous and heterozygous deletions as small as 100 kb with high resolution. A human oligonucleotide library was printed on amine binding slides. Arrays were hybridized using a hybstation and analysed using BlueFuse feature extraction software, with >95% of spots passing quality control. The protocol allows as little as 300 ng of input DNA and a 90% reduction of Cot-1 DNA without compromising quality. High quality results have also been obtained with DNA from archival tissue. Finally, in addition to human oligo arrays, we have applied the protocol successfully to mouse oligo arrays. We believe that this oligo-based platform using ‘off-the-shelf’ oligo libraries provides an easy accessible alternative to BAC arrays for CGH, which is cost-effective, available at high resolution and easily implemented for any sequenced organism without compromising the quality of the results.     

Combining Rational and Random Strategies in β-Glucosidase Zm-p60.1 Protein Library Construction

Saturation mutagenesis is a cornerstone technique in protein engineering because of its utility (in conjunction with appropriate analytical techniques) for assessing effects of varying residues at selected positions on proteins’ structures and functions. Site-directed mutagenesis with degenerate primers is the simplest and most rapid saturation mutagenesis technique. Thus, it is highly appropriate for assessing whether or not variation at certain sites is permissible, but not necessarily the most time- and cost-effective technique for detailed assessment of variations’ effects. Thus, in the presented study we applied the technique to randomize position W373 in β-glucosidase Zm-p60.1, which is highly conserved among β-glucosidases. Unexpectedly, β-glucosidase activity screening of the generated variants showed that most variants were active, although they generally had significantly lower activity than the wild type enzyme. Further characterization of the library led us to conclude that a carefully selected combination of randomized codon-based saturation mutagenesis and site-directed mutagenesis may be most efficient, particularly when constructing and investigating randomized libraries with high fractions of positive hits.     

TrimerDimer: an oligonucleotide-based saturation mutagenesis approach that removes redundant and stop codons

9-fluorenylmethoxycarbonyl (Fmoc) and 4,4′-dimethoxytrityl (DMTr) are orthogonal hydroxyl protecting groups that have been used in conjunction to assemble oligonucleotide libraries whose variants contain wild-type and mutant codons randomly interspersed throughout a focused DNA region. Fmoc is labile to organic bases and stable to weak acids, whereas DMTr behaves oppositely. Based on these chemical characteristics, we have now devised TrimerDimer, a novel codon-based saturation mutagenesis approach that removes redundant and stop codons during the assembly of degenerate oligonucleotides. In this approach, five DMTr-protected trinucleotide phosphoramidites (dTGG, dATG, dTTT, dTAT and dTGC) and five Fmoc-protected dinucleotide phosphoramidites (dAA, dTT, dAT, dGC and dCG) react simultaneously with a starting oligonucleotide growing on a solid support. The Fmoc group is then removed and the incorporated dimers react with a mixture of three DMTr-protected monomer phosphoramidites (dC, dA and dG) to produce 15 trinucleotides: dCAA, dAAA, dGAA, dCTT, dATT, dGTT, dCAT, dAAT, dGAT, dCGC, dAGC, dGGC, dCCG, dACG and dGCG. After one mutagenic cycle, 20 codons are generated encoding the 20 natural amino acids. TrimerDimer was tested by randomizing the four contiguous codons that encode amino acids L64–G67 of an engineered, nonfluorescent GFP protein. Sequencing of 89 nonfluorescent mutant clones and isolation of two fluorescent mutants confirmed the principle.     

Laboratory-Directed Protein Evolution

Systematic approaches to directed evolution of proteins have been documented since the 1970s. The ability to recruit new protein functions arises from the considerable substrate ambiguity of many proteins. The substrate ambiguity of a protein can be interpreted as the evolutionary potential that allows a protein to acquire new specificities through mutation or to regain function via mutations that differ from the original protein sequence. All organisms have evolutionarily exploited this substrate ambiguity. When exploited in a laboratory under controlled mutagenesis and selection, it enables a protein to “evolve” in desired directions. One of the most effective strategies in directed protein evolution is to gradually accumulate mutations, either sequentially or by recombination, while applying selective pressure. This is typically achieved by the generation of libraries of mutants followed by efficient screening of these libraries for targeted functions and subsequent repetition of the process using improved mutants from the previous screening. Here we review some of the successful strategies in creating protein diversity and the more recent progress in directed protein evolution in a wide range of scientific disciplines and its impacts in chemical, pharmaceutical, and agricultural sciences.     

Empirical Establishment of Oligonucleotide Probe Design Criteria

Criteria for the design of gene-specific and group-specific oligonucleotide probes were established experimentally via an oligonucleotide array that contained perfect match (PM) and mismatch probes (50-mers and 70-mers) based upon four genes. The effects of probe-target identity, continuous stretch, mismatch position, and hybridization free energy on specificity were tested. Little hybridization was observed at a probe-target identity of ≤85% for both 50-mer and 70-mer probes. PM signal intensities (33 to 48%) were detected at a probe-target identity of 94% for 50-mer oligonucleotides and 43 to 55% for 70-mer probes at a probe-target identity of 96%. When the effects of sequence identity and continuous stretch were considered independently, a stretch probe (>15 bases) contributed an additional 9% of the PM signal intensity compared to a nonstretch probe (≤15 bases) at the same identity level. Cross-hybridization increased as the length of continuous stretch increased. A 35-base stretch for 50-mer probes or a 50-base stretch for 70-mer probes had approximately 55% of the PM signal. Little cross-hybridization was observed for probes with a minimal binding free energy greater than −30 kcal/mol for 50-mer probes or −40 kcal/mol for 70-mer probes. Based on the experimental results, a set of criteria are suggested for the design of gene-specific and group-specific oligonucleotide probes, and the experimentally established criteria should provide valuable information for new software and algorithms for microarray-based studies.     

Fitness Loss and Library Size Determination in Saturation Mutagenesis

Saturation mutagenesis is a widely used directed evolution technique, in which a large number of protein variants, each having random amino acids in certain predetermined positions, are screened in order to discover high-fitness variants among them. Several metrics for determining the library size (the number of variants screened) have been suggested in the literature, but none of them incorporates the actual fitness of the variants discovered in the experiment. We present the results of an extensive simulation study, which is based on probabilistic models for protein fitness landscape, and which investigates how the result of a saturation mutagenesis experiment – the fitness of the best variant discovered – varies as a function of the library size. In particular, we study the loss of fitness in the experiment: the difference between the fitness of the best variant discovered, and the fitness of the best variant in variant space. Our results are that the existing criteria for determining the library size are conservative, so smaller libraries are often satisfactory. Reducing the library size can save labor, time, and expenses in the laboratory.     

A general method of site-specific mutagenesis using a modification of the Thermus aquaticus polymerase chain reaction

A specific mutagenic change in the cDNA of human protein S was introduced by a modification of the polymerase chain reaction that permits the introduction of a mutation at any position in a double-stranded DNA molecule. The method employed four synthetic oligonucleotide primers. One oligonucleotide contained a single-base mismatch to direct the mutagenesis; the other three oligonucleotides were designed to allow selective amplification of the mutated sequence with Thermus aquaticus polymerase. The mutagenized cDNA was cloned into a plasmid vector and transformed into Escherichia coli RR1 cells for characterization. The desired cytosine to guanine change in the target cDNA was confirmed by the predicted appearance of an AluI restriction site and by dideoxynucleotide sequencing. No other sequence changes were detected within the amplified region. This method of site-specific mutagenesis can be applied to any linear double-stranded DNA large enough for primer annealing and obviates specialized cloning vectors, DNA constructs, and selection techniques. It has the advantage over a recently published PCR technique (R. Higuchi, B. Krummel, and R. Saki (1988) Nucleic Acids Res. 16, 7351-7367) in requiring no diafiltration to remove primers between steps and in requiring only a single mutagenic oligonucleotide to be synthesized for each mutant construct made after the initial one.