Differences in the probability of cloning specific DNA between primary and amplified libraries: theoretical considerations.

We present a theoretical study of the probability of isolating a particular clone from a DNA library. There are differences in this probability between primary and amplified libraries even if the desired clone represents the same fraction of both libraries. As the result, we must screen severalfold more phage or bacteria in an amplified library than in a primary library.     

Polishing the craft of genetic diversity creation in directed evolution

Genetic diversity creation is a core technology in directed evolution where a high quality mutant library is crucial to its success. Owing to its importance, the technology in genetic diversity creation has seen rapid development over the years and its application has diversified into other fields of scientific research. The advances in molecular cloning and mutagenesis since 2008 were reviewed. Specifically, new cloning techniques were classified based on their principles of complementary overhangs, homologous sequences, overlapping PCR and megaprimers and the advantages, drawbacks and performances of these methods were highlighted. New mutagenesis methods developed for random mutagenesis, focused mutagenesis and DNA recombination were surveyed. The technical requirements of these methods and the mutational spectra were compared and discussed with references to commonly used techniques. The trends of mutant library preparation were summarised. Challenges in genetic diversity creation were discussed with emphases on creating “smart” libraries, controlling the mutagenesis spectrum and specific challenges in each group of mutagenesis methods. An outline of the wider applications of genetic diversity creation includes genome engineering, viral evolution, metagenomics and a study of protein functions. The review ends with an outlook for genetic diversity creation and the prospective developments that can have future impact in this field.     

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.     

Genetic approaches to the study of protein-protein interactions.

This article describes genetic approaches to the study of heterologous protein-protein interactions, focusing on the yeast Saccharomyces cerevisiae as a useful eukaryotic model system. Several methods are described that can be used to search for new interactions, including extragenic suppression, multicopy suppression, synthetic lethality, and transdominant inhibition. Strategies for screening, genetic characterization, and clone identification are described, along with recent examples from the literature. In addition, genetic methods are discussed that can be used to further characterize a newly discovered protein-protein interaction. These include the creation of mutant libraries of a given protein by chemical mutagenesis or polymerase chain reaction, the production of dominant-negative mutants, and strategies for introducing these mutant alleles back into yeast for analysis. Although these genetic methods are quite powerful, they are often just a starting point for further biochemical or cell biological experiments.     

User-friendly algorithms for estimating completeness and diversity in randomized protein-encoding libraries

Directed evolution of proteins depends on the production of molecular diversity by random mutagenesis. While a number of methods have been developed for introducing this diversity, the best ways to sample it are not always clear. Here we used simple statistics to analyse completeness and diversity in randomized libraries generated by oligonucleotide-directed mutagenesis, error-prone polymerase chain reaction (epPCR) and in vitro recombination of highly homologous sequences. For oligonucleotide-directed mutagenesis, we derive equations to estimate how complete a given library is expected to be and also to predict the size of library required to give a fixed probability of being 100% complete. We describe the statistical bases for computer programs which estimate the number of distinct variants represented in epPCR and shuffled libraries, dubbed PEDEL and DRIVeR, respectively. These programs allow the user to calculate (rather than guess) the diversity represented in a given library and also provide empirical guidelines for maximizing this diversity. PEDEL and DRIVeR are available at www.bio.cam.ac.uk/ approximately blackburn/stats.html.     

Unbiased libraries in protein directed evolution

Directed evolution is a powerful approach to study the molecular basis of protein evolution and to engineer proteins for a wide range of applications in synthetic organic chemistry and biotechnology. There are many methods based on random or focused mutagenesis to engineer successfully any protein trait. Focused approaches such as site-directed and saturation mutagenesis have become methods of choice for improving protein activity, selectivity, stability and many other traits because the screening step can be practically handled (bottleneck in directed evolution). Although novel mutagenesis methods based on CRISPR or solid-phase gene synthesis can eliminate bias when creating protein libraries, traditional PCR approaches, although imperfect, remain widely used due to their ease and low cost. One of the most common approaches in focused mutagenesis relies on NNK mutagenesis, however, the primer-based 22c-trick and small-intelligent methods have emerged as key tools for constructing less biased and unbiased libraries when all 20 canonical amino acids are needed for various reasons. In this minireview, we assess studies employing such methods for library creation and their areas of application. We also discuss the advantages and disadvantages of both methods and provide a perspective for creating smarter libraries.     

Calculating complexity of large randomized libraries

Randomized libraries are increasingly popular in protein engineering and other biomedical research fields. Statistics of the libraries are useful to guide and evaluate randomized library construction. Previous works only give the mean of the number of unique sequences in the library, and they can only handle equal molar ratio of the four nucleotides at a small number of mutation sites. We derive formulas to calculate the mean and variance of the number of unique sequences in libraries generated by cassette mutagenesis with mixtures of arbitrary nucleotide ratios. Computer program was developed which utilizes arbitrary numerical precision software package to calculate the statistics of large libraries. The statistics of library with mutations in more than 20 amino acids can be calculated easily. Results show that the nucleotide ratios have significant effects on these statistics. The more skewed the ratio, the larger the library size is needed to obtain the same expected number of unique sequences. The program is freely available at http://graphics.med.yale.edu/cgi-bin/lib_comp.pl.     

Mouse protein arrays from a TH1 cell cDNA library for antibody screening and serum profiling

The mouse is the premier genetic model organism for the study of disease and development. We describe the establishment of a mouse T helper cell type 1 (T(H)1) protein expression library that provides direct access to thousands of recombinant mouse proteins, in particular those associated with immune responses. The advantage of a system based on the combination of large cDNA expression libraries with microarray technology is the direct connection of the DNA sequence information from a particular clone to its recombinant, expressed protein. We have generated a mouse T(H)1 expression cDNA library and used protein arrays of this library to characterize the specificity and cross-reactivity of antibodies. Additionally, we have profiled the autoantibody repertoire in serum of a mouse model for systemic lupus erythematosus on these protein arrays and validated the putative autoantigens on highly sensitive protein microarrays.     

Dissecting the binding energy epitope of a high-affinity variant of human growth hormone: cooperative and additive effects from combining mutations from independently selected phage display mutagenesis libraries

Phage display mutagenesis is a widely used approach to engineering novel protein properties and is especially powerful in probing structure-function relationships in molecular recognition processes. The relative contributions of additive and cooperative binding forces and the influence of conformational diversity in producing a novel protein-protein interface is investigated using as a model an ultra-high-affinity receptor binding variant of human growth hormone (hGHv) that has been previously affinity matured. The modular aspect of how the mutations were grouped in the phage display libraries and combined allowed for a systematic probing of the inherent functional cross-talk between the different secondary structure elements that make up the remodeled hGHv binding surface. We performed an alanine scanning analyses of 35 hGHv residues and determined the kinetics of each variant by surface plasmon resonance (SPR). This analysis showed that there is a significant difference between the additive and cooperative binding forces existing among the selected residues in each library module, and the binding advantage of these residues is maximized over the original wild-type residue when in the context of the other mutations in the library. The degree to which residues in a particular mutagenesis library display binding cooperativity characteristics is generally correlated with the conformational plasticity of the polypeptide chain. Additionally, these cooperativity effects change when the mutations from one library are combined with the mutations from one or several of the other separate libraries. This supports the idea that significant functional cross-talk exists between the combined library modules that can affect the binding energetics of individual residues over a large distance.     

Exploiting recombination in single bacteria to make large phage antibody libraries

The creation of large phage antibody libraries has become an important goal in selecting antibodies against any antigen. Here we describe a method for making libraries so large that the complete diversity cannot be accessed using traditional phage technology. This involves the creation of a primary phage scFv library in a phagemid vector containing two nonhomologous lox sites. Contrary to the current dogma, we found that infecting Cre recombinase-expressing bacteria by such a primary library at a high multiplicity of infection results in the entry of many different phagemid into the cell. Exchange of Vh and Vl genes between such phagemids creates many new V h/Vl combinations, all of which are functional. On the basis of the observed recombination, the library is calculated to have a diversity of 3x1011. A library created using this method was validated by the selection of high affinity antibodies against a large number of different protein antigens.     

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.     

RNA mutagenesis yields highly diverse mRNA libraries for <Emphasis Type="Italic">in vitro</Emphasis>protein evolution

Background  

In protein drug development, in vitro molecular optimization or protein maturation can be used to modify protein properties. One basic approach to protein maturation is the introduction of random DNA mutations into the target gene sequence to produce a library of variants that can be screened for the preferred protein properties. Unfortunately, the capability of this approach has been restricted by deficiencies in the methods currently available for random DNA mutagenesis and library generation. Current DNA based methodologies generally suffer from nucleotide substitution bias that preferentially mutate particular base pairs or show significant bias with respect to transitions or transversions. In this report, we describe a novel RNA-based random mutagenesis strategy that utilizes Qβ replicase to manufacture complex mRNA libraries with a mutational spectrum that is close to the ideal.     

Construction and Characterization of Two Bacterial Artificial Chromosome Libraries of Zhikong Scallop, Chlamys farreri Jones et Preston, and Identification of BAC Clones Containing the Genes Involved in Its Innate Immune System

Large-insert bacterial artificial chromosome (BAC) libraries are necessary for advanced genetics and genomics research. To facilitate gene cloning and characterization, genome analysis, and physical mapping of scallop, two BAC libraries were constructed from nuclear DNA of Zhikong scallop, Chlamys farreri Jones et Preston. The libraries were constructed in the BamHI and MboI sites of the vector pECBAC1, respectively. The BamHI library consists of 73,728 clones, and approximately 99% of the clones contain scallop nuclear DNA inserts with an average size of 110 kb, covering 8.0x haploid genome equivalents. Similarly, the MboI library consists of 7680 clones, with an average insert of 145 kb and no insert-empty clones, thus providing a genome coverage of 1.1x. The combined libraries collectively contain a total of 81,408 BAC clones arrayed in 212 384-well microtiter plates, representing 9.1x haploid genome equivalents and having a probability of greater than 99% of discovering at least one positive clone with a single-copy sequence. High-density clone filters prepared from a subset of the two libraries were screened with nine pairs of Overgos designed from the cDNA or DNA sequences of six genes involved in the innate immune system of mollusks. Positive clones were identified for every gene, with an average of 5.3 BAC clones per gene probe. These results suggest that the two scallop BAC libraries provide useful tools for gene cloning, genome physical mapping, and large-scale sequencing in the species.     

Development of a new transformation-competent artificial chromosome (TAC) vector and construction of tomato and rice TAC libraries

Recent research has shown that BIBAC (binary bacterial artificial chromosome) and TAC (transformation-competent artificial chromosome) vector systems are very useful tools for map-based cloning of agronomically important genes in plant species. We have developed a new TAC vector that is suitable for both dicot and monocot transformation. Using this new TAC vector, we constructed large-insert genomic libraries of tomato and rice. The tomato library contains 96,996 clones (28.3-38.5 kb insert size) and has 3.18 haploid genome equivalents. The rice TAC library has 32.7 kb average insert size and has 9.24 haploid genome equivalents. The quality of these two libraries was tested using PCR to verify genome coverage. Individual clones were characterized to confirm insert integrity by Southern analysis, end sequencing and genetic mapping. To investigate the potential application of these TAC libraries in map-based cloning, TAC constructs containing a 45 kb fragment were introduced into the rice genome via Agrobacterium-mediated transformation. Molecular analysis indicates that the 45 kb fragment was successfully transferred into the rice genome. Although rearrangements of the introduced DNA were detected, 50% of regenerated plants contained at least one intact copy of the 45 kb clone and associated vector sequences. These libraries provide us with a valuable resource to rapidly isolate important genes in tomato and rice.     

Structure‐based design of combinatorial mutagenesis libraries

The development of protein variants with improved properties (thermostability, binding affinity, catalytic activity, etc.) has greatly benefited from the application of high‐throughput screens evaluating large, diverse combinatorial libraries. At the same time, since only a very limited portion of sequence space can be experimentally constructed and tested, an attractive possibility is to use computational protein design to focus libraries on a productive portion of the space. We present a general‐purpose method, called “Structure‐based Optimization of Combinatorial Mutagenesis ” (SOCoM ), which can optimize arbitrarily large combinatorial mutagenesis libraries directly based on structural energies of their constituents. SOCoM chooses both positions and substitutions, employing a combinatorial optimization framework based on library‐averaged energy potentials in order to avoid explicitly modeling every variant in every possible library. In case study applications to green fluorescent protein, β‐lactamase, and lipase A, SOCoM optimizes relatively small, focused libraries whose variants achieve energies comparable to or better than previous library design efforts, as well as larger libraries (previously not designable by structure‐based methods) whose variants cover greater diversity while still maintaining substantially better energies than would be achieved by representative random library approaches. By allowing the creation of large‐scale combinatorial libraries based on structural calculations, SOCoM promises to increase the scope of applicability of computational protein design and improve the hit rate of discovering beneficial variants. While designs presented here focus on variant stability (predicted by total energy), SOCoM can readily incorporate other structure‐based assessments, such as the energy gap between alternative conformational or bound states.     

'In vitro evolution' of ligands for HCV-specific serum antibodies

We developed a strategy to improve the properties of ligands selected from phage-displayed random peptide libraries. A site-directed mutagenesis protocol that introduces mutations and extends the size of a target sequence has been set up to generate diversity in a single or in a population of clones. The pool of mutants thus created is screened to identify variants with the desired properties. We refer to this strategy as in vitro evolution' of ligands. Here we report the application of this in vitro evolution protocol to the identification of improved ligands for HCV-specific serum antibodies. A single clone or population of clones were processed to generate a secondary library. Screening of these libraries with sera from HCV-infected patients identified peptides with an enhanced and broadened ability to detect HCV-specific serum antibodies.     

A mathematical model for experimental gene evolution

The purpose of this paper is to determine the optimal mutation rate for random mutagenesis procedures used to make mutant libraries for subsequent screening. When the mutation rate is low, the probability of achieving a rare beneficial mutation is low. When the mutation rate is high, the probability of producing lethal mutations which result in loss of function is also high. We demonstrate that between these two extremes, an optimal mutation rate exists for experimental gene improvement. This rate depends strongly on the number of simultaneous mutations required for a beneficial change of the gene, but only weakly on the number of possible lethal mutations. This model predicts that when mutagenesis is performed at the optimum mutation rate, at least 63% (1--e(-1)) of the cloned genes in a mutant library will be non-functional.     

Efficient isolation and mapping of rad genes of the fungus Coprinus cinereus using chromosome-specific libraries.

We have constructed cosmid libraries from electrophoretically separated chromosomes of the basidiomycete Coprinus cinereus. These libraries greatly facilitate the isolation of genes by complementation of mutant phenotypes and are particularly useful for map-based cloning strategies. From a library constructed from two co-migrating C.cinereus chromosomes, we isolated a clone that complements the C.cinereus rad9-1 mutation. Examination of this clone showed that it complements both the repair and meiotic defects of this mutant. Restriction fragment length polymorphism mapping using a portion of this clone showed that it maps to the rad9 locus. In addition, a single copy of transforming DNA is sufficient to complement the rad9-1 defects. Thus, we believe we have cloned the rad9 gene itself. We also used a chromosome-specific library and backcrossed isolates to rapidly identify a cosmid clone which is tightly linked to the rad11 locus and is therefore a suitable starting point for a chromosome walk. These rapid methods of gene mapping and isolation should be applicable to any organism with separable chromosomes.     

Construction of "small-intelligent" focused mutagenesis libraries using well-designed combinatorial degenerate primers

Site-saturation mutagenesis is a powerful tool for protein optimization due to its efficiency and simplicity. A degenerate codon NNN or NNS (K) is often used to encode the 20 standard amino acids, but this will produce redundant codons and cause uneven distribution of amino acids in the constructed library. Here we present a novel "small-intelligent" strategy to construct mutagenesis libraries that have a minimal gene library size without inherent amino acid biases, stop codons, or rare codons of Escherichia coli by coupling well-designed combinatorial degenerate primers with suitable PCR-based mutagenesis methods. The designed primer mixture contains exactly one codon per amino acid and thus allows the construction of small-intelligent mutagenesis libraries with one gene per protein. In addition, the software tool DC-Analyzer was developed to assist in primer design according to the user-defined randomization scheme for library construction. This small-intelligent strategy was successfully applied to the randomization of halohydrin dehalogenases with one or two randomized sites. With the help of DC-Analyzer, the strategy was proven to be as simple as NNS randomization and could serve as a general tool to efficiently randomize target genes at positions of interest.