Incorporating Synthetic Oligonucleotides via Gene Reassembly (ISOR): a versatile tool for generating targeted libraries

Insertions and Deletions (Indels): A Missing Piece of the Protein Engineering Jigsaw

Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.

Dan Tawfik's Lessons for Protein Engineers about Enzymes Adapting to New Substrates

Natural evolution has been creating new complex systems for billions of years. The process is spontaneous and requires neither intelligence nor moral purpose but is nevertheless difficult to understand. The late Dan Tawfik spent years studying enzymes as they adapted to recognize new substrates. Much of his work focused on gaining fundamental insights, so the practical utility of his experiments may not be obvious even to accomplished protein engineers. Here we focus on two questions fundamental to any directed evolution experiment. Which proteins are the best starting points for such experiments? Which trait(s) of the chosen parental protein should be evolved to achieve the desired outcome? We summarize Tawfik's contributions to our understanding of these problems, to honor his memory and encourage those unfamiliar with his ideas to read his publications.

Pairing of single mutations yields obligate Cre-type site-specific recombinases

Tyrosine site-specific recombinases (SSRs) represent a versatile genome editing tool with considerable therapeutic potential. Recent developments to engineer and evolve SSRs into heterotetramers to improve target site flexibility signified a critical step towards their broad utility in genome editing. However, SSR monomers can form combinations of different homo- and heterotetramers in cells, increasing their off-target potential. Here, we discover that two paired mutations targeting residues implicated in catalysis lead to simple obligate tyrosine SSR systems, where the presence of all distinct subunits to bind as a heterotetramer is obligatory for catalysis. Therefore, only when the paired mutations are applied as single mutations on each recombinase subunit, the engineered SSRs can efficiently recombine the intended target sequence, while the subunits carrying the point mutations expressed in isolation are inactive. We demonstrate the utility of the obligate SSR system to improve recombination specificity of a designer-recombinase for a therapeutic target in human cells. Furthermore, we show that the mutations render the naturally occurring SSRs, Cre and Vika, obligately heteromeric for catalytic proficiency, providing a straight-forward approach to improve their applied properties. These results facilitate the development of safe and effective therapeutic designer-recombinases and advance our mechanistic understanding of SSR catalysis.

Directed Evolution: Methodologies and Applications

Directed evolution aims to expedite the natural evolution process of biological molecules and systems in a test tube through iterative rounds of gene diversifications and library screening/selection. It has become one of the most powerful and widespread tools for engineering improved or novel functions in proteins, metabolic pathways, and even whole genomes. This review describes the commonly used gene diversification strategies, screening/selection methods, and recently developed continuous evolution strategies for directed evolution. Moreover, we highlight some representative applications of directed evolution in engineering nucleic acids, proteins, pathways, genetic circuits, viruses, and whole cells. Finally, we discuss the challenges and future perspectives in directed evolution.

Accessing unexplored regions of sequence space in directed enzyme evolution via insertion/deletion mutagenesis

Insertions and deletions (InDels) are frequently observed in natural protein evolution, yet their potential remains untapped in laboratory evolution. Here we introduce a transposon-based mutagenesis approach (TRIAD) to generate libraries of random variants with short in-frame InDels, and screen TRIAD libraries to evolve a promiscuous arylesterase activity in a phosphotriesterase. The evolution exhibits features that differ from previous point mutagenesis campaigns: while the average activity of TRIAD variants is more compromised, a larger proportion has successfully adapted for the activity. Different functional profiles emerge: (i) both strong and weak trade-off between activities are observed; (ii) trade-off is more severe (20- to 35-fold increased kcat/KM in arylesterase with 60-400-fold decreases in phosphotriesterase activity) and (iii) improvements are present in kcat rather than just in KM, suggesting adaptive solutions. These distinct features make TRIAD an alternative to widely used point mutagenesis, accessing functional innovations and traversing unexplored fitness landscape regions.

A heterodimer of evolved designer-recombinases precisely excises a human genomic DNA locus

Site-specific recombinases (SSRs) such as the Cre/loxP system are useful genome engineering tools that can be repurposed by altering their DNA-binding specificity. However, SSRs that delete a natural sequence from the human genome have not been reported thus far. Here, we describe the generation of an SSR system that precisely excises a 1.4 kb fragment from the human genome. Through a streamlined process of substrate-linked directed evolution we generated two separate recombinases that, when expressed together, act as a heterodimer to delete a human genomic sequence from chromosome 7. Our data indicates that designer-recombinases can be generated in a manageable timeframe for precision genome editing. A large-scale bioinformatics analysis suggests that around 13% of all human protein-coding genes could be targetable by dual designer-recombinase induced genomic deletion (dDRiGD). We propose that heterospecific designer-recombinases, which work independently of the host DNA repair machinery, represent an efficient and safe alternative to nuclease-based genome editing technologies.

Affinity maturation of antibodies by combinatorial codon mutagenesis versus error-prone PCR

IN VITRO

affinity maturation of therapeutic monoclonal antibodies is commonly applied to achieve desired properties, such as improved binding kinetics and affinity. Currently there are no universally accepted protocols for generation of variegated antibody libraries or selection thereof. Here, we performed affinity maturation using a yeast-based single-chain variable fragment (scFv) expression system to compare two mutagenesis methods: random mutagenesis across the entire V(D)J region by error-prone PCR, and a novel combinatorial mutagenesis process limited to the complementarity-determining regions (CDRs). We applied both methods of mutagenesis to four human antibodies against well-known immuno-oncology target proteins. Detailed sequence analysis showed an even mutational distribution across the entire length of the scFv for the error-prone PCR method and an almost exclusive targeting of the CDRs for the combinatorial method. Though there were distinct mutagenesis profiles for each target antibody and mutagenesis method, we found that both methods improved scFv affinity with similar efficiency. When a subset of the affinity-matured antibodies was expressed as full-length immunoglobulin, the measured affinity constants were mostly comparable to those of the respective scFv, but the full-length antibodies were inferior to their scFv counterparts for one of the targets. Furthermore, we found that improved affinity for the full-length antibody did not always translate into enhanced binding to cell-surface expressed antigen or improved immune checkpoint blocking ability, suggesting that screening with full-length antibody or antigen-binding fragment formats might be advantageous and the subject of a future study.

Golden Mutagenesis: An efficient multi-site-saturation mutagenesis approach by Golden Gate cloning with automated primer design

Site-directed methods for the generation of genetic diversity are essential tools in the field of directed enzyme evolution. The Golden Gate cloning technique has been proven to be an efficient tool for a variety of cloning setups. The utilization of restriction enzymes which cut outside of their recognition domain allows the assembly of multiple gene fragments obtained by PCR amplification without altering the open reading frame of the reconstituted gene. We have developed a protocol, termed Golden Mutagenesis that allows the rapid, straightforward, reliable and inexpensive construction of mutagenesis libraries. One to five amino acid positions within a coding sequence could be altered simultaneously using a protocol which can be performed within one day. To facilitate the implementation of this technique, a software library and web application for automated primer design and for the graphical evaluation of the randomization success based on the sequencing results was developed. This allows facile primer design and application of Golden Mutagenesis also for laboratories, which are not specialized in molecular biology.

Evolution of chalcone isomerase from a noncatalytic ancestor

The emergence of catalysis in a noncatalytic protein scaffold is a rare, unexplored event. Chalcone isomerase (CHI), a key enzyme in plant flavonoid biosynthesis, is presumed to have evolved from a nonenzymatic ancestor related to the widely distributed fatty-acid binding proteins (FAPs) and a plant protein family with no isomerase activity (CHILs). Ancestral inference supported the evolution of CHI from a protein lacking isomerase activity. Further, we identified four alternative founder mutations, i.e., mutations that individually instated activity, including a mutation that is not phylogenetically traceable. Despite strong epistasis in other cases of protein evolution, CHI's laboratory reconstructed mutational trajectory shows weak epistasis. Thus, enantioselective CHI activity could readily emerge despite a catalytically inactive starting point. Accordingly, X-ray crystallography, NMR, and molecular dynamics simulations reveal reshaping of the active site toward a productive substrate-binding mode and repositioning of the catalytic arginine that was inherited from the ancestral fatty-acid binding proteins.

Evolution of cyclohexadienyl dehydratase from an ancestral solute-binding protein

The emergence of enzymes through the neofunctionalization of noncatalytic proteins is ultimately responsible for the extraordinary range of biological catalysts observed in nature. Although the evolution of some enzymes from binding proteins can be inferred by homology, we have a limited understanding of the nature of the biochemical and biophysical adaptations along these evolutionary trajectories and the sequence in which they occurred. Here we reconstructed and characterized evolutionary intermediate states linking an ancestral solute-binding protein to the extant enzyme cyclohexadienyl dehydratase. We show how the intrinsic reactivity of a desolvated general acid was harnessed by a series of mutations radiating from the active site, which optimized enzyme-substrate complementarity and transition-state stabilization and minimized sampling of noncatalytic conformations. Our work reveals the molecular evolutionary processes that underlie the emergence of enzymes de novo, which are notably mirrored by recent examples of computational enzyme design and directed evolution.

Directed evolution of SIRT6 for improved deacylation and glucose homeostasis maintenance

Mammalian SIRT6 is a well-studied histone deacetylase that was recently shown to exhibit high protein deacylation activity enabling the removal of long chain fatty acyl groups from proteins. SIRT6 was shown to play key roles in cellular homeostasis by regulating a variety of cellular processes including DNA repair and glucose metabolism. However, the link between SIRT6 enzymatic activities and its cellular functions is not clear. Here, we utilized a directed enzyme evolution approach to generate SIRT6 mutants with improved deacylation activity. We found that while two mutants show increased deacylation activity at high substrate concentration and improved glucose metabolism they exhibit no improvement and even abolished deacetylation activity on H3K9Ac and H3K56Ac in cells. Our results demonstrate the separation of function between SIRT6 catalytic activities and suggest that SIRT6 deacylation activity in cells is important for glucose metabolism and can be mediated by still unknown acylated cellular proteins.

A generic HTS assay for kinase screening: Validation for the isolation of an engineered malate kinase

An end-point ADP/NAD+ acid/alkali assay procedure, directly applicable to library screening of any type of ATP-utilising/ADP producing enzyme activity, was implemented. Typically, ADP production is coupled to NAD+ co-enzyme formation by the conventional addition of pyruvate kinase and lactate dehydrogenase. Transformation of enzymatically generated NAD+ into a photometrically active alkali derivative product is then achieved through the successive application of acidic/alkali treatment steps. The assay was successfully miniaturized to search for malate kinase activity in a structurally-guided library of LysC aspartate kinase variants comprising 6,700 clones. The screening procedure enabled the isolation of nine positive variants showing novel kinase activity on (L)-malate, the best mutant, LysC V115A:E119S:E434V exhibited strong substrate selectivity for (L)-malate compared to (L)-aspartate with a (kcat/Km)malate/(kcat/Km)aspartate ratio of 86. Double mutants V115A:E119S, V115A:E119C and E119S:E434V were constructed to further probe the origins of stabilising substrate binding energy gains for (L)-malate due to mutation. The introduction of less sterically hindering side-chains in engineered enzymes carrying E119S and V115A mutations increases the effective volume available for substrate binding in the catalytic pocket. Improved binding of the (L)-malate substrate may be assisted by less hindered movement of the Phe184 aromatic side-chain. Additional favourable long-range electostatic effects on binding arising from the E434V surface mutation are conditionally dependent upon the presence of the V115A mutation close to Phe184 in the active-site.

Combinatorial Evolution of DNA with RECODE

In past decades, DNA engineering protocols have led to the rapid development of synthetic biology. To engineer the natural proteins, many directed evolution methods based on molecular biology have been presented for generating genetic diversity or obtaining specific properties. Here, we provide a simple (PCR operation), efficient (larger amount of products), and powerful (multiple point mutations, deletions, insertions, and combinatorial multipoint mutagenesis) RECODE method, which is capable of reediting the target DNA flexibly to restructure regulatory regions and remodel enzymes by using the combined function of the thermostable DNA polymerase and DNA ligase in one pot. RECODE is expected to be an applicable choice to create diverse mutant libraries for rapid evolution and optimization of enzymes and synthetic pathways.

Rational engineering of a native hyperthermostable lactonase into a broad spectrum phosphotriesterase

The redesign of enzyme active sites to alter their function or specificity is a difficult yet appealing challenge. Here we used a structure-based design approach to engineer the lactonase SsoPox from Sulfolobus solfataricus into a phosphotriesterase. The five best variants were characterized and their structure was solved. The most active variant, αsD6 (V27A-Y97W-L228M-W263M) demonstrates a large increase in catalytic efficiencies over the wild-type enzyme, with increases of 2,210-fold, 163-fold, 58-fold, 16-fold against methyl-parathion, malathion, ethyl-paraoxon, and methyl-paraoxon, respectively. Interestingly, the best mutants are also capable of degrading fensulfothion, which is reported to be an inhibitor for the wild-type enzyme, as well as others that are not substrates of the starting template or previously reported W263 mutants. The broad specificity of these engineered variants makes them promising candidates for the bioremediation of organophosphorus compounds. Analysis of their structures reveals that the increase in activity mainly occurs through the destabilization of the active site loop involved in substrate binding, and it has been observed that the level of disorder correlates with the width of the enzyme specificity spectrum. This finding supports the idea that active site conformational flexibility is essential to the acquisition of broader substrate specificity.

Construction of a synthetic metabolic pathway for biosynthesis of the non-natural methionine precursor 2,4-dihydroxybutyric acid

2,4-Dihydroxybutyric acid (DHB) is a molecule with considerable potential as a versatile chemical synthon. Notably, it may serve as a precursor for chemical synthesis of the methionine analogue 2-hydroxy-4-(methylthio)butyrate, thus, targeting a considerable market in animal nutrition. However, no natural metabolic pathway exists for the biosynthesis of DHB. Here we have therefore conceived a three-step metabolic pathway for the synthesis of DHB starting from the natural metabolite malate. The pathway employs previously unreported malate kinase, malate semialdehyde dehydrogenase and malate semialdehyde reductase activities. The kinase and semialdehyde dehydrogenase activities were obtained by rational design based on structural and mechanistic knowledge of candidate enzymes acting on sterically cognate substrates. Malate semialdehyde reductase activity was identified from an initial screening of several natural enzymes, and was further improved by rational design. The pathway was expressed in a minimally engineered Escherichia coli strain and produces 1.8 g l^-1 DHB with a molar yield of 0.15.

A Simple Combinatorial Codon Mutagenesis Method for Targeted Protein Engineering

Directed evolution is a powerful tool for optimizing enzymes, and mutagenesis methods that improve enzyme library quality can significantly expedite the evolution process. Here, we report a simple method for targeted combinatorial codon mutagenesis (CCM). To demonstrate the utility of this method for protein engineering, CCM libraries were constructed for cytochrome P450_BM3, pfu prolyl oligopeptidase, and the flavin-dependent halogenase RebH; 10-26 sites were targeted for codon mutagenesis in each of these enzymes, and libraries with a tunable average of 1-7 codon mutations per gene were generated. Each of these libraries provided improved enzymes for their respective transformations, which highlights the generality, simplicity, and tunability of CCM for targeted protein engineering.

Directed evolution of a soluble human DR3 receptor for the inhibition of TL1A induced cytokine secretion

TNF-like 1A (TL1A) is a cytokine belonging to the TNF superfamily that promotes inflammation in autoimmune diseases. Inhibiting the interaction of TL1A with the endogenous death-domain receptor 3 (DR3) offers a therapeutic approach for treating TL1A-induced autoimmune diseases. Here, we generated improved DR3 variants showing increased TL1A binding affinity and stability using a directed evolution approach. Given the high cysteine content and post-translational modification of DR3, we employed yeast surface display and expression in mammalian cell lines for screening, expression and characterization of improved DR3 variants. A cell-based assay performed with the human TF-1 cell line and CD4⁺ T cells showed that two improved DR3 mutants efficiently inhibited TL1A-induced cell death and secretion of IFN-γ, respectively. These DR3 mutants can be used as drug candidates for the treatment of inflammatory bowel diseases and for other autoimmune diseases, including rheumatic arthritis and asthma.

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.

The evolutionary pathway from a biologically inactive polypeptide sequence to a folded, active structural mimic of DNA

The protein Ocr (overcome classical restriction) from bacteriophage T7 acts as a mimic of DNA and inhibits all Type I restriction/modification (RM) enzymes. Ocr is a homodimer of 116 amino acids and adopts an elongated structure that resembles the shape of a bent 24 bp DNA molecule. Each monomer includes 34 acidic residues and only six basic residues. We have delineated the mimicry of Ocr by focusing on the electrostatic contribution of its negatively charged amino acids using directed evolution of a synthetic form of Ocr, termed pocr, in which all of the 34 acidic residues were substituted for a neutral amino acid. In vivo analyses confirmed that pocr did not display any antirestriction activity. Here, we have subjected the gene encoding pocr to several rounds of directed evolution in which codons for the corresponding acidic residues found in Ocr were specifically re-introduced. An in vivo selection assay was used to detect antirestriction activity after each round of mutation. Our results demonstrate the variation in importance of the acidic residues in regions of Ocr corresponding to different parts of the DNA target which it is mimicking and for the avoidance of deleterious effects on the growth of the host.

Combinatorial Evolution of Enzymes and Synthetic Pathways Using One-Step PCR

DNA engineering is the fundamental motive driving the rapid development of modern biotechnology. Here, we present a versatile evolution method termed "rapidly efficient combinatorial oligonucleotides for directed evolution" (RECODE) for rapidly introducing multiple combinatorial mutations to the target DNA by combined action of a thermostable high-fidelity DNA polymerase and a thermostable DNA Ligase in one reaction system. By applying this method, we rapidly constructed a variant library of the rpoS promoters (with activity of 8-460%), generated a novel heparinase from the highly specific leech hyaluronidase (with more than 30 mutant residues) and optimized the heme biosynthetic pathway by combinatorial evolution of regulatory elements and pathway enzymes (2500 ± 120 mg L(-1) with 20-fold increase). The simple RECODE method enabled researchers the unparalleled ability to efficiently create diverse mutant libraries for rapid evolution and optimization of enzymes and synthetic pathways.

Methods for the directed evolution of proteins

Directed evolution has proved to be an effective strategy for improving or altering the activity of biomolecules for industrial, research and therapeutic applications. The evolution of proteins in the laboratory requires methods for generating genetic diversity and for identifying protein variants with desired properties. This Review describes some of the tools used to diversify genes, as well as informative examples of screening and selection methods that identify or isolate evolved proteins. We highlight recent cases in which directed evolution generated enzymatic activities and substrate specificities not known to exist in nature.

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.

SwiftLib: rapid degenerate-codon-library optimization through dynamic programming

Degenerate codon (DC) libraries efficiently address the experimental library-size limitations of directed evolution by focusing diversity toward the positions and toward the amino acids (AAs) that are most likely to generate hits; however, manually constructing DC libraries is challenging, error prone and time consuming. This paper provides a dynamic programming solution to the task of finding the best DCs while keeping the size of the library beneath some given limit, improving on the existing integer-linear programming formulation. It then extends the algorithm to consider multiple DCs at each position, a heretofore unsolved problem, while adhering to a constraint on the number of primers needed to synthesize the library. In the two library-design problems examined here, the use of multiple DCs produces libraries that very nearly cover the set of desired AAs while still staying within the experimental size limits. Surprisingly, the algorithm is able to find near-perfect libraries where the ratio of amino-acid sequences to nucleic-acid sequences approaches 1; it effectively side-steps the degeneracy of the genetic code. Our algorithm is freely available through our web server and solves most design problems in about a second.

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.

Generating targeted libraries by the combinatorial incorporation of synthetic oligonucleotides during gene shuffling (ISOR)

Protein engineering by directed evolution relies on the use of libraries enriched with beneficial variants. Such libraries should explore large mutational diversities while avoiding high loads of deleterious mutations. Here we describe a simple protocol for incorporating synthetic oligonucleotides that encode designed, site-specific mutations by assembly PCR. This protocol enables a researcher to "hedge the bets," namely, to explore a large number of potentially beneficial mutations in a combinatorial manner such that individual library variants carry a limited number of mutations.

Probabilistic methods in directed evolution: library size, mutation rate, and diversity

Directed evolution has emerged as an important tool for engineering proteins with improved or novel properties. Because of their inherent reliance on randomness, directed evolution protocols are amenable to probabilistic modeling and analysis. This chapter summarizes and reviews in a nonmathematical way some of the probabilistic works related to directed evolution, with particular focus on three of the most widely used methods: saturation mutagenesis, error-prone PCR, and in vitro recombination. The ultimate aim is to provide the reader with practical information to guide the planning and design of directed evolution studies. Importantly, the applications and locations of freely available computational resources to assist with this process are described in detail.

Bioprospecting in the genomic age

The genomic revolution promises great advances in the search for useful biocatalysts. Function-based metagenomic approaches have identified several enzymes with properties that make them useful candidates for a variety of bioprocesses. As DNA sequencing costs continue to decline, the volume of genomic data, along with their corresponding predicted protein sequences, will continue to increase dramatically, necessitating new approaches to leverage this information for gene-based bioprospecting efforts. Additionally, as new functions are discovered and correlated with this sequence information, the knowledge of the often complex relationship between a protein's sequence and function will improve. This in turn will lead to better gene-based bioprospecting approaches and facilitate the tailoring of desired properties through protein engineering projects. In this chapter, we discuss a number of recent advances in bioprospecting within the context of the genomic age.

Protein engineering by random mutagenesis and structure-guided consensus of Geobacillus stearothermophilus Lipase T6 for enhanced stability in methanol

The abilities of enzymes to catalyze reactions in nonnatural environments of organic solvents have opened new opportunities for enzyme-based industrial processes. However, the main drawback of such processes is that most enzymes have a limited stability in polar organic solvents. In this study, we employed protein engineering methods to generate a lipase for enhanced stability in methanol, which is important for biodiesel production. Two protein engineering approaches, random mutagenesis (error-prone PCR) and structure-guided consensus, were applied in parallel on an unexplored lipase gene from Geobacillus stearothermophilus T6. A high-throughput colorimetric screening assay was used to evaluate lipase activity after an incubation period in high methanol concentrations. Both protein engineering approaches were successful in producing variants with elevated half-life values in 70% methanol. The best variant of the random mutagenesis library, Q185L, exhibited 23-fold-improved stability, yet its methanolysis activity was decreased by one-half compared to the wild type. The best variant from the consensus library, H86Y/A269T, exhibited 66-fold-improved stability in methanol along with elevated thermostability (+4.3°C) and a 2-fold-higher fatty acid methyl ester yield from soybean oil. Based on in silico modeling, we suggest that the Q185L substitution facilitates a closed lid conformation that limits access for both the methanol and substrate excess into the active site. The enhanced stability of H86Y/A269T was a result of formation of new hydrogen bonds. These improved characteristics make this variant a potential biocatalyst for biodiesel production.

Engineering V-type nerve agents detoxifying enzymes using computationally focused libraries

VX and its Russian (RVX) and Chinese (CVX) analogues rapidly inactivate acetylcholinesterase and are the most toxic stockpile nerve agents. These organophosphates have a thiol leaving group with a choline-like moiety and are hydrolyzed very slowly by natural enzymes. We used an integrated computational and experimental approach to increase Brevundimonas diminuta phosphotriesterase's (PTE) detoxification rate of V-agents by 5000-fold. Computational models were built of the complex between PTE and V-agents. On the basis of these models, the active site was redesigned to be complementary in shape to VX and RVX and to include favorable electrostatic interactions with their choline-like leaving group. Small libraries based on designed sequences were constructed. The libraries were screened by a direct assay for V-agent detoxification, as our initial studies showed that colorimetric surrogates fail to report the detoxification rates of the actual agents. The experimental results were fed back to improve the computational models. Overall, five rounds of iterating between experiment and model refinement led to variants that hydrolyze the toxic SP isomers of all three V-agents with kcat/KM values of up to 5 × 10(6) M(-1) min(-1) and also efficiently detoxify G-agents. These new catalysts provide the basis for broad spectrum nerve agent detoxification.

Combinatorial engineering of dextransucrase specificity

We used combinatorial engineering to investigate the relationships between structure and linkage specificity of the dextransucrase DSR-S from Leuconostoc mesenteroides NRRL B-512F, and to generate variants with altered specificity. Sequence and structural analysis of glycoside-hydrolase family 70 enzymes led to eight amino acids (D306, F353, N404, W440, D460, H463, T464 and S512) being targeted, randomized by saturation mutagenesis and simultaneously recombined. Screening of two libraries totaling 3.6.10⁴ clones allowed the isolation of a toolbox comprising 81 variants which synthesize high molecular weight α-glucans with different proportions of α(1→3) linkages ranging from 3 to 20 %. Mutant sequence analysis, biochemical characterization and molecular modelling studies revealed the previously unknown role of peptide ⁴⁶⁰DYVHT⁴⁶⁴ in DSR-S linkage specificity. This peptide sequence together with residue S512 contribute to defining +2 subsite topology, which may be critical for the enzyme regiospecificity.

Directed evolution of a soluble human IL-17A receptor for the inhibition of psoriasis plaque formation in a mouse model

Interleukin-17 (IL-17) is a T-cell-derived cytokine that promotes inflammatory pathology in autoimmune diseases. Blocking IL-17A interactions with its endogenous IL-17 receptor (IL-17RA) can constitute an important target for therapeutic intervention. Here, we utilized a directed evolution approach to generate soluble IL-17RA mutants that exhibit increased IL-17A binding affinity and thermostability, relative to the wild-type. Human fibroblast cell-based assay and in vivo analysis in mice indicated that two improved IL-17RA mutants efficiently inhibit the secretion of IL-17A-induced proinflammatory cytokines. Analysis of one of these mutants in a psoriasis mouse model showed its efficacy in promoting the recovery of psoriasis plaques. This mutant can be used as a promising drug candidate for the treatment of psoriasis and may be a therapeutic agent for various other autoimmune diseases.

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.

The evolutionary origins of detoxifying enzymes: the mammalian serum paraoxonases (PONs) relate to bacterial homoserine lactonases

Serum paraoxonases (PONs) are detoxifying lactonases that were first identified in mammals. Three mammalian families are known, PON1, 2, and 3 that reside primarily in the liver. They catalyze essentially the same reaction, lactone hydrolysis, but differ in their substrate specificity. Although some members are highly specific, others have a broad specificity profile. The evolutionary origins and substrate specificities of PONs therefore remain poorly understood. Here, we report a newly identified family of bacterial PONs, and the reconstruction of the ancestor of the three families of mammalian PONs. Both the mammalian ancestor and the characterized bacterial PONX_OCCAL were found to efficiently hydrolyze N-acyl homoserine lactones that mediate quorum sensing in many bacteria, including pathogenic ones. The mammalian PONs may therefore relate to a newly identified family of bacterial, PON-like "quorum-quenching" lactonases. The appearance of PONs in metazoa is likely to relate to innate immunity rather than detoxification. Unlike the bacterial PON, the mammalian ancestor also hydrolyzes, with low efficiency, lactones other than homoserine lactones, thus preceding the detoxifying functions that diverged later in two of the three mammalian families. The bifunctionality of the mammalian ancestor and the trade-off between the quorum-quenching and detoxifying lactonase activities explain the broad and overlapping specificities of some mammalian PONs versus the singular specificity of others.

Integrating the protein and metabolic engineering toolkits for next-generation chemical biosynthesis

Through microbial engineering, biosynthesis has the potential to produce thousands of chemicals used in everyday life. Metabolic engineering and synthetic biology are fields driven by the manipulation of genes, genetic regulatory systems, and enzymatic pathways for developing highly productive microbial strains. Fundamentally, it is the biochemical characteristics of the enzymes themselves that dictate flux through a biosynthetic pathway toward the product of interest. As metabolic engineers target sophisticated secondary metabolites, there has been little recognition of the reduced catalytic activity and increased substrate/product promiscuity of the corresponding enzymes compared to those of central metabolism. Thus, fine-tuning these enzymatic characteristics through protein engineering is paramount for developing high-productivity microbial strains for secondary metabolites. Here, we describe the importance of protein engineering for advancing metabolic engineering of secondary metabolism pathways. This pathway integrated enzyme optimization can enhance the collective toolkit of microbial engineering to shape the future of chemical manufacturing.

What makes a protein fold amenable to functional innovation? Fold polarity and stability trade-offs

Protein evolvability includes two elements--robustness (or neutrality, mutations having no effect) and innovability (mutations readily inducing new functions). How are these two conflicting demands bridged? Does the ability to bridge them relate to the observation that certain folds, such as TIM barrels, accommodate numerous functions, whereas other folds support only one? Here, we hypothesize that the key to innovability is polarity--an active site composed of flexible, loosely packed loops alongside a well-separated, highly ordered scaffold. We show that highly stabilized variants of TEM-1 β-lactamase exhibit selective rigidification of the enzyme's scaffold while the active-site loops maintained their conformational plasticity. Polarity therefore results in stabilizing, compensatory mutations not trading off, but instead promoting the acquisition of new activities. Indeed, computational analysis indicates that in folds that accommodate only one function throughout evolution, for example, dihydrofolate reductase, ≥ 60% of the active-site residues belong to the scaffold. In contrast, folds associated with multiple functions such as the TIM barrel show high scaffold-active-site polarity (~20% of the active site comprises scaffold residues) and >2-fold higher rates of sequence divergence at active-site positions. Our work suggests structural measures of fold polarity that appear to be correlated with innovability, thereby providing new insights regarding protein evolution, design, and engineering.

Evolved stereoselective hydrolases for broad-spectrum G-type nerve agent detoxification

A preferred strategy for preventing nerve agents intoxication is catalytic scavenging by enzymes that hydrolyze them before they reach their targets. Using directed evolution, we simultaneously enhanced the activity of a previously described serum paraoxonase 1 (PON1) variant for hydrolysis of the toxic S(P) isomers of the most threatening G-type nerve agents. The evolved variants show ≤340-fold increased rates and catalytic efficiencies of 0.2-5 × 10(7) M(-1) min(-1). Our selection for prevention of acetylcholinesterase inhibition also resulted in the complete reversion of PON1's stereospecificity, from an enantiomeric ratio (E) < 6.3 × 10(-4) in favor of the R(P) isomer of a cyclosarin analog in wild-type PON1, to E > 2,500 for the S(P) isomer in an evolved variant. Given their ability to hydrolyze G-agents, these evolved variants may serve as broad-range G-agent prophylactics.