Steering directed protein evolution: strategies to manage combinatorial complexity of mutant libraries

Structural bioinformatics studies of glutamate transporters and their AlphaFold2 predicted water-soluble QTY variants and uncovering the natural mutations of L->Q, I->T, F->Y and Q->L, T->I and Y->F

Glutamate transporters play key roles in nervous physiology by modulating excitatory neurotransmitter levels, when malfunctioning, involving in a wide range of neurological and physiological disorders. However, integral transmembrane proteins including the glutamate transporters remain notoriously difficult to study, due to their localization within the cell membrane. Here we present the structural bioinformatics studies of glutamate transporters and their water-soluble variants generated through QTY-code, a protein design strategy based on systematic amino acid substitutions. These include 2 structures determined by X-ray crystallography, cryo-EM, and 6 predicted by AlphaFold2, and their predicted water-soluble QTY variants. In the native structures of glutamate transporters, transmembrane helices contain hydrophobic amino acids such as leucine (L), isoleucine (I), and phenylalanine (F). To design water-soluble variants, these hydrophobic amino acids are systematically replaced by hydrophilic amino acids, namely glutamine (Q), threonine (T) and tyrosine (Y). The QTY variants exhibited water-solubility, with four having identical isoelectric focusing points (pI) and the other four having very similar pI. We present the superposed structures of the native glutamate transporters and their water-soluble QTY variants. The superposed structures displayed remarkable similarity with RMSD 0.528Å-2.456Å, despite significant protein transmembrane sequence differences (41.1%->53.8%). Additionally, we examined the differences of hydrophobicity patches between the native glutamate transporters and their QTY variants. Upon closer inspection, we discovered multiple natural variations of L->Q, I->T, F->Y and Q->L, T->I, Y->F in these transporters. Some of these natural variations were benign and the remaining were reported in specific neurological disorders. We further investigated the characteristics of hydrophobic to hydrophilic substitutions in glutamate transporters, utilizing variant analysis and evolutionary profiling. Our structural bioinformatics studies not only provided insight into the differences between the hydrophobic helices and hydrophilic helices in the glutamate transporters, but they are also expected to stimulate further study of other water-soluble transmembrane proteins.

Combinatorial InVitroFlow-assisted Mutagenesis (CombIMut) yields a 41-fold improved CelA2 cellulase

The combination of diversity generation methods and ultrahigh-throughput screening (uHTS) technologies is key to efficiently explore nature's sequence space and elucidate structure-function relationships of enzymes. Beneficial substitutions often cluster in a few regions and simultaneous amino acid substitutions at multiple positions (e.g., by OmniChange) will likely lead to further improved enzyme variants. An extensive screening effort is required to identify such variants, as the simultaneous randomization of four codons can easily yield over 10⁵ potential enzyme variants. The combination of flow cytometer-based uHTS with cell-free compartmentalization technology using (w/o/w) double emulsions (InVitroFlow), provides analysis capabilities of up to 10⁷ events per hour, thus enabling efficient screening. InVitroFlow is an elegant solution since diversity loss through a transformation of host cells is omitted and emulsion compartments provide a genotype-phenotype linkage through a fluorescence readout. In this work, a multi-site saturation mutagenesis (mSSM) and an OmniChange library with four simultaneously saturated positions in the active site of CelA2 cellulase were screened using InVitroFlow. Screening of over 36 million events, yielded a significantly improved cellulase variant CelA2-M3 (H288F/H524Q) with an 8-fold increase in specific activity compared to the parent CelA2-H288F (83.9 U/mg) and a 41-fold increased specific activity (674.5 U/mg) compared to wildtype CelA2 (16.6 U/mg) for the substrate 4-MUC (4-methylumbelliferyl-β D-cellobioside). This article is protected by copyright. All rights reserved.

Methods for enzyme library creation: Which one will you choose?: A guide for novices and experts to introduce genetic diversity

Enzyme engineering allows to explore sequence diversity in search for new properties. The scientific literature is populated with methods to create enzyme libraries for engineering purposes, however, choosing a suitable method for the creation of mutant libraries can be daunting, in particular for the novices. Here, we address both novices and experts: how can one enter the arena of enzyme library design and what guidelines can advanced users apply to select strategies best suited to their purpose? Section I is dedicated to the novices and presents an overview of established and standard methods for library creation, as well as available commercial solutions. The expert will discover an up-to-date tool to freshen up their repertoire (Section I) and learn of the newest methods that are likely to become a mainstay (Section II). We focus primarily on in vitro methods, presenting the advantages of each method. Our ultimate aim is to offer a selection of methods/strategies that we believe to be most useful to the enzyme engineer, whether a first-timer or a seasoned user.

Random Mutagenesis by Insertion of Error-Prone PCR Products to the Chromosome of Bacillus subtilis

Bacillus subtilis is an attractive host for the directed evolution of the enzymes whose substrates cannot be transported across cell membrane. However, the generation of a mutant library in B. subtilis suffers problems of small library size, plasmid instability, and heterozygosity. Here, a large library of random mutant was created by inserting error-prone PCR (epPCR) products to the chromosome of B. subtilis. Specifically, the epPCR product was fused with flanking regions and antibiotic resistant marker using a PCR-based multimerization method, generating insertion construct. The epPCR product was integrated into the chromosome via homologous recombination after the insertion construct was transformed into the supercompetent cells of B. subtilis strain SCK6. The transformation efficiency of the insertion construct was improved through co-expressing homologous recombination-promoting protein NgAgo, raising the number of competent cells, and increasing the length of flanking regions. A library containing 5.31 × 105 random mutants was constructed using per μg insertion construct, which is sufficient for directed evolution. The library generation process was accomplished within 1 day. The effectiveness of this method was confirmed by improving the activity of Methyl Parathion Hydrolase (MPH) toward chlorpyrifos and by enhancing the secretion level of MPH in B. subtilis. Taken together, the present work provides a fast and efficient method to integrate epPCR products into the chromosome of B. subtilis, facilitating directed evolution and expression optimization of target proteins.

Holistic engineering of Cal-A lipase chain-length selectivity identifies triglyceride binding hot-spot

Through the application of a region-focused saturation mutagenesis and randomization approach, protein engineering of the Cal-A enzyme was undertaken with the goal of conferring new triglyceride selectivity. Little is known about the mode of triglyceride binding to Cal-A. Engineering Cal-A thus requires a systemic approach. Targeted and randomized Cal-A libraries were created, recombined using the Golden Gate approach and screened to detect variants able to discriminate between long-chain (olive oil) and short-chain (tributyrin) triglyceride substrates using a high-throughput in vivo method to visualize hydrolytic activity. Discriminative variants were analyzed using an in-house script to identify predominant substitutions. This approach allowed identification of variants that exhibit strong discrimination for the hydrolysis of short-chain triglycerides and others that discriminate towards hydrolysis of long-chain triglycerides. A clear pattern emerged from the discriminative variants, identifying the 217–245 helix-loop-helix motif as being a hot-spot for triglyceride recognition. This was the consequence of introducing the entire mutational load in selected regions, without putting a strain on distal parts of the protein. Our results improve our understanding of the Cal-A lipase mode of action and selectivity. This holistic perspective to protein engineering, where parts of the gene are individually mutated and the impact evaluated in the context of the whole protein, can be applied to any protein scaffold.

High-Throughput Screening Assays for Lipolytic Enzymes

Screening is defined as the identification of hits within a large library of variants of an enzyme or protein with a predefined property. In theory, each variant present in the respective library needs to be assayed; however, to save time and consumables, many screening regimes involve a primary round to identify clones producing active enzymes. Such primary or prescreenings for lipolytic enzyme activity are often carried out on agar plates containing pH indicators or substrates as triolein or tributyrin. Subsequently, high-throughput screening assays are usually performed in microtiter plate (MTP) format using chromogenic or fluorogenic substrates and, if available, automated liquid handling robotics. Here, we describe different assay systems to determine the activity and enantioselectivity of lipases and esterases as well as the synthesis of several substrates. We also report on the construction of a complete site saturation library derived from lipase A of Bacillus subtilis and its testing for detergent tolerance. This approach allows for the identification of amino acids affecting sensitivity or resistance against different detergents.

Novel form of the Michaelis-Menten equation that enables accurate estimation of (kcat/KM)*KI with just two rate measurements; utility in directed evolution

One of applications of directed evolution is to desensitize an enzyme to an inhibitor. kcat,1/KM and KI are three dimensions that when multiplied measure an enzyme's intrinsic capacity for catalysis in the presence of an inhibitor. The ideal values for the individual dimensions depend on substrate and inhibitor concentrations under the conditions of the application. When attempting to optimize those values by directed evolution, (kcat/KM)*KI can be an informative parameter for evaluating libraries of variants, but throughput is limited. We describe a manipulation of the Michaelis-Menten equation for competitive inhibition that isolates (kcat/KM)*KI on one side of the equation. If velocity is measured at constant enzyme and substrate concentrations with two different inhibitor concentrations (one of which can be 0), the data are sufficient to calculate (kcat/KM)*KI with just two rate measurements. The procedure is validated by correlating values obtained by the rapid method with those obtained by substrate saturation kinetics.

Enzyme engineering: A synthetic biology approach for more effective library generation and automated high-throughput screening

The Golden Gate strategy entails the use of type IIS restriction enzymes, which cut outside of their recognition sequence. It enables unrestricted design of unique DNA fragments that can be readily and seamlessly recombined. Successfully employed in other synthetic biology applications, we demonstrate its advantageous use to engineer a biocatalyst. Hot-spots for mutations were individuated in three distinct regions of Candida antarctica lipase A (Cal-A), the biocatalyst chosen as a target to demonstrate the versatility of this recombination method. The three corresponding gene segments were subjected to the most appropriate method of mutagenesis (targeted or random). Their straightforward reassembly allowed combining products of different mutagenesis methods in a single round for rapid production of a series of diverse libraries, thus facilitating directed evolution. Screening to improve discrimination of short-chain versus long-chain fatty acid substrates was aided by development of a general, automated method for visual discrimination of the hydrolysis of varied substrates by whole cells.

ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules

The degree of evolutionary conservation of an amino acid in a protein or a nucleic acid in DNA/RNA reflects a balance between its natural tendency to mutate and the overall need to retain the structural integrity and function of the macromolecule. The ConSurf web server (http://consurf.tau.ac.il), established over 15 years ago, analyses the evolutionary pattern of the amino/nucleic acids of the macromolecule to reveal regions that are important for structure and/or function. Starting from a query sequence or structure, the server automatically collects homologues, infers their multiple sequence alignment and reconstructs a phylogenetic tree that reflects their evolutionary relations. These data are then used, within a probabilistic framework, to estimate the evolutionary rates of each sequence position. Here we introduce several new features into ConSurf, including automatic selection of the best evolutionary model used to infer the rates, the ability to homology-model query proteins, prediction of the secondary structure of query RNA molecules from sequence, the ability to view the biological assembly of a query (in addition to the single chain), mapping of the conservation grades onto 2D RNA models and an advanced view of the phylogenetic tree that enables interactively rerunning ConSurf with the taxa of a sub-tree.

A Robust and Versatile Method of Combinatorial Chemical Synthesis of Gene Libraries via Hierarchical Assembly of Partially Randomized Modules

A major challenge in gene library generation is to guarantee a large functional size and diversity that significantly increases the chances of selecting different functional protein variants. The use of trinucleotides mixtures for controlled randomization results in superior library diversity and offers the ability to specify the type and distribution of the amino acids at each position. Here we describe the generation of a high diversity gene library using tHisF of the hyperthermophile Thermotoga maritima as a scaffold. Combining various rational criteria with contingency, we targeted 26 selected codons of the thisF gene sequence for randomization at a controlled level. We have developed a novel method of creating full-length gene libraries by combinatorial assembly of smaller sub-libraries. Full-length libraries of high diversity can easily be assembled on demand from smaller and much less diverse sub-libraries, which circumvent the notoriously troublesome long-term archivation and repeated proliferation of high diversity ensembles of phages or plasmids. We developed a generally applicable software tool for sequence analysis of mutated gene sequences that provides efficient assistance for analysis of library diversity. Finally, practical utility of the library was demonstrated in principle by assessment of the conformational stability of library members and isolating protein variants with HisF activity from it. Our approach integrates a number of features of nucleic acids synthetic chemistry, biochemistry and molecular genetics to a coherent, flexible and robust method of combinatorial gene synthesis.

The Mutagenesis Assistant Program

Mutagenesis Assistant Program (MAP) is a web-based statistical tool to develop directed evolution strategies by investigating the consequences at the amino acid level of the mutational biases of random mutagenesis methods on any given gene. The latest development of the program, the MAP(2.0)3D server, correlates the generated amino acid substitution patterns of a specific random mutagenesis method to the sequence and structural information of the target protein. The combined information can be used to select an experimental strategy that improves the chances of obtaining functionally efficient and/or stable enzyme variants. Hence, the MAP(2.0)3D server facilitates the "in silico" prescreening of the target gene by predicting the amino acid diversity generated in a random mutagenesis library. Here, we describe the features of MAP(2.0)3D server by analyzing, as an example, the cytochrome P450BM3 monooxygenase (CYP102A1). The MAP(2.0)3D server is available publicly at http://map.jacobs-university.de/map3d.html.

Synthetic biology for the directed evolution of protein biocatalysts: navigating sequence space intelligently

The amino acid sequence of a protein affects both its structure and its function. Thus, the ability to modify the sequence, and hence the structure and activity, of individual proteins in a systematic way, opens up many opportunities, both scientifically and (as we focus on here) for exploitation in biocatalysis. Modern methods of synthetic biology, whereby increasingly large sequences of DNA can be synthesised de novo, allow an unprecedented ability to engineer proteins with novel functions. However, the number of possible proteins is far too large to test individually, so we need means for navigating the 'search space' of possible protein sequences efficiently and reliably in order to find desirable activities and other properties. Enzymologists distinguish binding (Kd) and catalytic (kcat) steps. In a similar way, judicious strategies have blended design (for binding, specificity and active site modelling) with the more empirical methods of classical directed evolution (DE) for improving kcat (where natural evolution rarely seeks the highest values), especially with regard to residues distant from the active site and where the functional linkages underpinning enzyme dynamics are both unknown and hard to predict. Epistasis (where the 'best' amino acid at one site depends on that or those at others) is a notable feature of directed evolution. The aim of this review is to highlight some of the approaches that are being developed to allow us to use directed evolution to improve enzyme properties, often dramatically. We note that directed evolution differs in a number of ways from natural evolution, including in particular the available mechanisms and the likely selection pressures. Thus, we stress the opportunities afforded by techniques that enable one to map sequence to (structure and) activity in silico, as an effective means of modelling and exploring protein landscapes. Because known landscapes may be assessed and reasoned about as a whole, simultaneously, this offers opportunities for protein improvement not readily available to natural evolution on rapid timescales. Intelligent landscape navigation, informed by sequence-activity relationships and coupled to the emerging methods of synthetic biology, offers scope for the development of novel biocatalysts that are both highly active and robust.

Engineering the ligninolytic enzyme consortium

The ligninolytic enzyme consortium is one of the most-efficient oxidative systems found in nature, playing a pivotal role during wood decay and coal formation. Typically formed by high redox-potential oxidoreductases, this array of enzymes can be used within the emerging lignocellulose biorefineries in processes that range from the production of bioenergy to that of biomaterials. To ensure that these versatile enzymes meet industry standards and needs, they have been subjected to directed evolution and hybrid approaches that surpass the limits imposed by nature. This Opinion article analyzes recent achievements in this field, including the incipient groundbreaking research into the evolution of resurrected enzymes, and the engineering of ligninolytic secretomes to create consolidated bioprocessing microbes with synthetic biology applications.

Directed evolution 2.0: improving and deciphering enzyme properties

A KnowVolution: knowledge gaining directed evolution including four phases is proposed in this feature article, which generates improved enzyme variants and molecular understanding.

A high-throughput screening method to reengineer DNA polymerases for random mutagenesis

A screening system for directed evolution of DNA polymerases employing a fluorescent Scorpion probe as a reporter has been developed. The screening system has been validated in a directed evolution experiment of a distributive polymerase from the Y-polymerase family (Dpo4 from Sulfolobus solfataricus) which was improved in elongation efficiency of consecutive mismatches. The engineering campaign yielded improved Dpo4 polymerase variants one of which was successfully benchmarked in a sequence saturation mutagenesis experiment especially with regard to the desirable consecutive transversion mutations ([2.5-fold increase in frequency relative to a reference library prepared with Dpo4 WT). The Scorpion probe screening system enables to reengineer polymerases with low processivity and fidelity, and no secondary activities (i.e. exonuclease activity or strand displacement activity) to match demands in diversity generation for directed protein evolution.

Mutagenic Organized Recombination Process by Homologous IN vivo Grouping (MORPHING) for directed enzyme evolution

Approaches that depend on directed evolution require reliable methods to generate DNA diversity so that mutant libraries can focus on specific target regions. We took advantage of the high frequency of homologous DNA recombination in Saccharomyces cerevisiae to develop a strategy for domain mutagenesis aimed at introducing and in vivo recombining random mutations in defined segments of DNA. Mutagenic Organized Recombination Process by Homologous IN vivo Grouping (MORPHING) is a one-pot random mutagenic method for short protein regions that harnesses the in vivo recombination apparatus of yeast. Using this approach, libraries can be prepared with different mutational loads in DNA segments of less than 30 amino acids so that they can be assembled into the remaining unaltered DNA regions in vivo with high fidelity. As a proof of concept, we present two eukaryotic-ligninolytic enzyme case studies: i) the enhancement of the oxidative stability of a H₂O₂-sensitive versatile peroxidase by independent evolution of three distinct protein segments (Leu28-Gly57, Leu149-Ala174 and Ile199-Leu268); and ii) the heterologous functional expression of an unspecific peroxygenase by exclusive evolution of its native 43-residue signal sequence.

Computational tools for designing smart libraries

Traditional directed evolution experiments are often time-, labor- and cost-intensive because they involve repeated rounds of random mutagenesis and the selection or screening of large mutant libraries. The efficiency of directed evolution experiments can be significantly improved by targeting mutagenesis to a limited number of hot-spot positions and/or selecting a limited set of substitutions. The design of such "smart" libraries can be greatly facilitated by in silico analyses and predictions. Here we provide an overview of computational tools applicable for (a) the identification of hot-spots for engineering enzyme properties, and (b) the evaluation of predicted hot-spots and selection of suitable amino acids for substitutions. The selected tools do not require any specific expertise and can easily be implemented by the wider scientific community.

Computer-Aided Protein Directed Evolution: a Review of Web Servers, Databases and other Computational Tools for Protein Engineering

The combination of computational and directed evolution methods has proven a winning strategy for protein engineering. We refer to this approach as computer-aided protein directed evolution (CAPDE) and the review summarizes the recent developments in this rapidly growing field. We will restrict ourselves to overview the availability, usability and limitations of web servers, databases and other computational tools proposed in the last five years. The goal of this review is to provide concise information about currently available computational resources to assist the design of directed evolution based protein engineering experiment.

Improving the thermostability of alpha-amylase by combinatorial coevolving-site saturation mutagenesis

Background

The generation of focused mutant libraries at hotspot residues is an important strategy in directed protein evolution. Existing methods, such as combinatorial active site testing and residual coupling analysis, depend primarily on the evolutionary conserved information to find the hotspot residues. Hardly any attention has been paid to another important functional and structural determinants, the functionally correlated variation information--coevolution.

Results

In this paper, we suggest a new method, named combinatorial coevolving-site saturation mutagenesis (CCSM), in which the functionally correlated variation sites of proteins are chosen as the hotspot sites to construct focused mutant libraries. The CCSM approach was used to improve the thermal stability of α-amylase from Bacillus subtilis CN7 (Amy7C). The results indicate that the CCSM can identify novel beneficial mutation sites, and enhance the thermal stability of wild-type Amy7C by 8°C ( T5030), which could not be achieved with the ordinarily rational introduction of single or a double point mutation.

Conclusions

Our method is able to produce more thermostable mutant α-amylases with novel beneficial mutations at new sites. It is also verified that the coevolving sites can be used as the hotspots to construct focused mutant libraries in protein engineering. This study throws new light on the active researches of the molecular coevolution.

Mixed oligonucleotides for random mutagenesis: best way of making them

The generation of proteins, especially enzymes, with pre-deliberated, novel properties is a big challenge in the field of protein engineering. This aim, over the years was critically facilitated by newly emerging methods of combinatorial and evolutionary techniques, such as combinatorial gene synthesis followed by functional screening of many structural variants generated in parallel (library). Libraries can be generated by a large number of available methods. Therein the use of mixtures of pre-formed trinucleotide blocks representing codons for the 20 canonical amino acids for oligonucleotide synthesis stands out as allowing fully controlled partial (or total) randomization individually at any number of arbitrarily chosen codon positions of a given gene. This has created substantial demand of fully protected trinucleotide synthons of good reactivity in standard oligonucleotide synthesis. We here review methods for the preparation of oligonucleotide mixtures with a strong focus on codon-specific trinucleotide blocks.

Saccharomyces cerevisiae in directed evolution: An efficient tool to improve enzymes

Over the past 20 years, directed evolution has been seen to be the most reliable approach to protein engineering. Emulating the natural selection algorithm, ad hoc enzymes with novel features can be tailor-made for practical purposes through iterative rounds of random mutagenesis, DNA recombination and screening. Of the heterologous hosts used in laboratory evolution experiments, the budding yeast Saccharomyces cerevisiae has become the best choice to express eukaryotic proteins with improved properties. S. cerevisiae not only allows mutant enzymes to be secreted but also, it permits a wide range of genetic manipulations to be employed, ranging from in vivo cloning to the creation of greater molecular diversity, thanks to its efficient DNA recombination apparatus. Here, we summarize some successful examples of the use of the S. cerevisiae machinery to accelerate artificial evolution, complementing the traditional in vitro methods to generate tailor-made enzymes.

MAP(2.0)3D: a sequence/structure based server for protein engineering

The Mutagenesis Assistant Program (MAP) is a web-based tool to provide statistical analyses of the mutational biases of directed evolution experiments on amino acid substitution patterns. MAP analysis assists protein engineers in the benchmarking of random mutagenesis methods that generate single nucleotide mutations in a codon. Herein, we describe a completely renewed and improved version of the MAP server, the MAP(2.0)3D server, which correlates the generated amino acid substitution patterns to the structural information of the target protein. This correlation aids in the selection of a more suitable random mutagenesis method with specific biases on amino acid substitution patterns. In particular, the new server represents MAP indicators on secondary and tertiary structure and correlates them to specific structural components such as hydrogen bonds, hydrophobic contacts, salt bridges, solvent accessibility, and crystallographic B-factors. Three model proteins (D-amino oxidase, phytase, and N-acetylneuraminic acid aldolase) are used to illustrate the novel capability of the server. MAP(2.0)3D server is available publicly at http://map.jacobs-university.de/map3d.html.

A new and convenient approach for the preparation of β-cyanoethyl protected trinucleotide phosphoramidites

Herein we report a convenient approach for the preparation of fully protected trinucleotide synthons to be used for the synthesis of gene libraries. The trinucleotide synthons bear β-cyanoethyl groups at the phosphate residues, and thus can be used in standard oligonucleotide synthesis without additional steps for deprotection and work-up.

OmniChange: the sequence independent method for simultaneous site-saturation of five codons

Focused mutant library generation methods have been developed to improve mainly "localizable" enzyme properties such as activity and selectivity. Current multi-site saturation methods are restricted by the gene sequence, require subsequent PCR steps and/or additional enzymatic modifications. Here we report, a multiple site saturation mutagenesis method, OmniChange, which simultaneously and efficiently saturates five independent codons. As proof of principle, five chemically cleaved DNA fragments, each carrying one NNK-degenerated codon, were generated and assembled to full gene length in a one-pot-reaction without additional PCR-amplification or use of restriction enzymes or ligases. Sequencing revealed the presence of up to 27 different codons at individual positions, corresponding to 84.4% of the theoretical diversity offered by NNK-degeneration. OmniChange is absolutely sequence independent, does not require a minimal distance between mutated codons and can be accomplished within a day.

High-throughput metabolic engineering: advances in small-molecule screening and selection

Metabolic engineering for the overproduction of high-value small molecules is dependent upon techniques in directed evolution to improve production titers. The majority of small molecules targeted for overproduction are inconspicuous and cannot be readily obtained by screening. We provide a review on the development of high-throughput colorimetric, fluorescent, and growth-coupled screening techniques, enabling inconspicuous small-molecule detection. We first outline constraints on throughput imposed during the standard directed evolution workflow (library construction, transformation, and screening) and establish a screening and selection ladder on the basis of small-molecule assay throughput and sensitivity. An in-depth analysis of demonstrated screening and selection approaches for small-molecule detection is provided. Particular focus is placed on in vivo biosensor-based detection methods that reduce or eliminate in vitro assay manipulations and increase throughput. We conclude by providing our prospectus for the future, focusing on transcription factor-based detection systems as a natural microbial mode of small-molecule detection.

Generation of random mutant libraries with multiple primers in a single reaction

Characterization of multiple sites in a single gene that are important in biological phenotypes is challenging due to the difficulty to generate many mutants representing all or a majority of combinations of mutations in the gene. Using the HIV-1 env and pol genes as templates, four random libraries were generated representing different combinations of mutations introduced by up to 36 mutagenesis primers in a single assay. Over 86% of the clones contained mutations and the mutants tended to have single or fewer mutations in the libraries. When protein size was used as a screening marker, all identified clones contained at least 2 mutations and up to 12 mutations were detected in a single clone. Nearly all mutant clones in each library contained unique mutations, indicating that mutants in the library were generated at random. Closely related mutations which were overlapped by neighboring mutagenesis primers were often introduced in this system. Analysis of the env library showed that some potential N-linked glycosylation sites did not increase the Env molecular mass significantly, suggesting they were not used for glycosylation or only limited carbohydrate moieties were added at these sites. This novel method can serve as a powerful tool to study the biological phenotypes of genes whose functions are determined by multiple sites.

Enhancing stress resistance and production phenotypes through transcriptome engineering

As Saccharomyces cerevisiae is engineered further as a microbial factory for industrially relevant but potentially cytotoxic molecules such as ethanol, issues of cell viability arise that threaten to place a biological limit on output capacity and/or the use of less refined production conditions. Evidence suggests that one naturally evolved mode of survival in deleterious environments involves the complex, multigenic interplay between disparate stress response and homeostasis mechanisms. Rational engineering of such resistance would require a systems-level understanding of cellular behavior that is, in general, not yet available. To circumvent this limitation, we have developed a phenotype discovery approach termed global transcription machinery engineering (gTME) that allows for the generation and selection of nonphysiological traits. We alter gene expression on a genome-wide scale by selecting for dominant mutations in a randomly mutagenized general transcription factor. The gene encoding the mutated transcription factor resides on a plasmid in a strain carrying the unaltered chromosomal allele. Thus, although the dominant mutations may destroy the essential function of the plasmid-borne variant, alteration of the transcriptome with minimal perturbation to normal cellular processes is possible via the presence of the native genomic allele. Achieving a phenotype of interest involves the construction and diversity evaluation of yeast libraries harboring random sequence variants of a chosen transcription factor and the subsequent selection and validation of mutant strains. We describe the rationale and procedures associated with each step in the context of generating strains possessing enhanced ethanol tolerance.

USER friendly DNA recombination (USERec): a simple and flexible near homology-independent method for gene library construction

USER friendly DNA recombination (USERec) is introduced as a near homology-independent method that allows the simultaneous recombination of an unprecedented number of 10 DNA fragments (approximately 40-400 bp) within a day. The large number of fragments and their ease of preparation enables the creation of libraries of much larger genetic diversity (potentially approximately 10(10)-10(11) sequences) than current alternative methods based on DNA truncation (ITCHY, SCRATCHY and SHIPREC) or type IIb restriction enzymes (SISDC). At the same time, the frequency of frameshifts in the recombined library is low (90% of the recombined sequences are in frame). Compared to overlap extension PCR, USERec also requires much reduced crossover sequence constraints (only a 5'-AN(4-8)T-3' motif) and fewer experimental steps. Based on its simplicity and flexibility, and the accessibility of large and high quality recombined DNA libraries, USERec is established as a convenient alternative for the combinatorial assembly of gene fragments (e.g. exon or domain shuffling) and for a number of applications in gene library construction, such as loop grafting and multi-site-directed or random mutagenesis.

Directed evolution of chemotaxis inhibitory protein of Staphylococcus aureus generates biologically functional variants with reduced interaction with human antibodies

Chemotaxis inhibitory protein of Staphylococcus aureus (CHIPS) is a protein that binds and blocks the C5a receptor (C5aR) and formylated peptide receptor, thereby inhibiting the immune cell recruitment associated with inflammation. If CHIPS was less reactive with existing human antibodies, it would be a promising anti-inflammatory drug candidate. Therefore, we applied directed evolution and computational/rational design to the CHIPS gene in order to generate new CHIPS variants displaying lower interaction with human IgG, yet retaining biological function. The optimization was performed in four rounds: one round of random mutagenesis to add diversity into the CHIPS gene and three rounds of DNA recombination by Fragment INduced Diversity (FIND). Every round was screened by phage selection and/or ELISA for decreased interaction with human IgG and retained C5aR binding. The mean binding of human anti-CHIPS IgG decreased with every round of evolution. For further optimization, new amino acid substitutions were introduced by rational design, based on the mutations identified during directed evolution. Finally, seven CHIPS variants with low interaction with human IgG and retained C5aR blocking capacity could be identified.

Directed enzyme evolution: climbing fitness peaks one amino acid at a time

Directed evolution can generate a remarkable range of new enzyme properties. Alternate substrate specificities and reaction selectivities are readily accessible in enzymes from families that are naturally functionally diverse. Activities on new substrates can be obtained by improving variants with broadened specificities or by step-wise evolution through a sequence of more and more challenging substrates. Evolution of highly specific enzymes has been demonstrated, even with positive selection alone. It is apparent that many solutions exist for any given problem, and there are often many paths that lead uphill, one step at a time.

In silico modelling of directed evolution: Implications for experimental design and stepwise evolution

We model the process of directed evolution (DE) in silico using genetic algorithms. Making use of the NK fitness landscape model, we analyse the effects of mutation rate, crossover and selection pressure on the performance of DE. A range of values of K, the epistatic interaction of the landscape, are considered, and high- and low-throughput modes of evolution are compared. Our findings suggest that for runs of or around ten generations' duration-as is typical in DE-there is little difference between the way in which DE needs to be configured in the high- and low-throughput regimes, nor across different degrees of landscape epistasis. In all cases, a high selection pressure (but not an extreme one) combined with a moderately high mutation rate works best, while crossover provides some benefit but only on the less rugged landscapes. These genetic algorithms were also compared with a "model-based approach" from the literature, which uses sequential fixing of the problem parameters based on fitting a linear model. Overall, we find that purely evolutionary techniques fare better than do model-based approaches across all but the smoothest landscapes.

GLUE-IT and PEDEL-AA: new programmes for analyzing protein diversity in randomized libraries

There are many methods for introducing random mutations into nucleic acid sequences. Previously, we described a suite of programmes for estimating the completeness and diversity of randomized DNA libraries generated by a number of these protocols. Our programmes suggested some empirical guidelines for library design; however, no information was provided regarding library diversity at the protein (rather than DNA) level. We have now updated our web server, enabling analysis of translated libraries constructed by site-saturation mutagenesis and error-prone PCR (epPCR). We introduce GLUE-Including Translation (GLUE-IT), which finds the expected amino acid completeness of libraries in which up to six codons have been independently varied (according to any user-specified randomization scheme). We provide two tools for assisting with experimental design: CodonCalculator, for assessing amino acids corresponding to given randomized codons; and AA-Calculator, for finding degenerate codons that encode user-specified sets of amino acids. We also present PEDEL-AA, which calculates amino acid statistics for libraries generated by epPCR. Input includes the parent sequence, overall mutation rate, library size, indel rates and a nucleotide mutation matrix. Output includes amino acid completeness and diversity statistics, and the number and length distribution of sequences truncated by premature termination codons. The web interfaces are available at http://guinevere.otago.ac.nz/stats.html.

Advances in laboratory evolution of enzymes

We address recent developments in the area of laboratory, or directed evolution, with a focus on enzymes and on new methodologies of generic potential. We survey three main areas: (i) library making techniques, including the application of computational and rational methods for library design; (ii) screening and selection techniques, including recent applications of enzyme screening by FACS (fluorescence activated cell sorter); (iii) new approaches for performing directed evolution, and in particular, the application of 'neutral drifts' (libraries generated by rounds of mutation and selection for the enzyme's original function) and of consensus mutations to generate highly evolvable starting points for directed evolution.

Enzyme optimization: moving from blind evolution to statistical exploration of sequence-function space

Directed evolution is a powerful tool for the creation of commercially useful enzymes, particularly those approaches that are based on in vitro recombination methods, such as DNA shuffling. Although these types of search algorithms are extraordinarily efficient compared with purely random methods, they do not explicitly represent or interrogate the genotype-phenotype relationship and are essentially blind in nature. Recently, however, researchers have begun to apply multivariate statistical techniques to model protein sequence-function relationships and guide the evolutionary process by rapidly identifying beneficial diversity for recombination. In conjunction with state-of-the-art library generation methods, the statistical approach to sequence optimization is now being used routinely to create enzymes efficiently for industrial applications.