Site-saturation Libraries Research Articles

Directed Evolution of Acoustic Reporter Genes Using High-Throughput Acoustic Screening.

A major challenge in the fields of biological imaging and synthetic biology is noninvasively visualizing the functions of natural and engineered cells inside opaque samples such as living animals. One promising technology that addresses this limitation is ultrasound (US), with its penetration depth of several cm and spatial resolution on the order of 100 μm. Within the past decade, reporter genes for US have been introduced and engineered to link cellular functions to US signals via heterologous expression in commensal bacteria and mammalian cells. These acoustic reporter genes (ARGs) represent a novel class of genetically encoded US contrast agent, and are based on air-filled protein nanostructures called gas vesicles (GVs). Just as the discovery of fluorescent proteins was followed by the improvement and diversification of their optical properties through directed evolution, here we describe the evolution of GVs as acoustic reporters. To accomplish this task, we establish high-throughput, semiautomated acoustic screening of ARGs in bacterial cultures and use it to screen mutant libraries for variants with increased nonlinear US scattering. Starting with scanning site saturation libraries for two homologues of the primary GV structural protein, GvpA/B, two rounds of evolution resulted in GV variants with 5- and 14-fold stronger acoustic signals than the parent proteins. We anticipate that this and similar approaches will help high-throughput protein engineering play as large a role in the development of acoustic biomolecules as it has for their fluorescent counterparts.

Engineering cytochrome P450s for selective alkene to carbonyl oxidation.

The Wacker-Tsuji oxidation is an important aerobic oxidation process to synthesize ethanal from ethene and methyl ketones from 1-alkenes. Current challenges in aerobic alkene oxidation include selective carbonyl product formation beyond methyl ketones. This includes the regioselective oxidation of the terminal carbon atom of 1-alkenes, the regioselective ketone formation with internal alkenes as well as the enantioselective alkene to carbonyl oxidation. Recently, the potential of high-valent metal-oxo species for direct alkene to carbonyl oxidation was explored as carbonyl product formation is frequently reported as a side reaction of alkene epoxidation by cytochrome P450s. It was shown that such promiscuous P450s can be engineered via directed evolution to perform alkene to carbonyl oxidation reactions with high activity and selectivity. Here, we report a protocol to convert promiscuous P450s into efficient and selective enzymes for Wacker-type alkene oxidation. One round of directed evolution is described in detail, which includes the generation and handling of site-saturation libraries, recombinant protein expression, library screening in a 96-well plate format and the rescreening of variants with beneficial mutations. These protocols might be useful to engineer various P450s for selective alkene to carbonyl oxidation, and to engineer enzymes in general.

"NAD-display": Ultrahigh-Throughput in Vitro Screening of NAD(H) Dehydrogenases Using Bead Display and Flow Cytometry.

NAD(H)‐utiliing enzymes have been the subject of directed evolution campaigns to improve their function. To enable access to a larger swath of sequence space, we demonstrate the utility of a cell‐free, ultrahigh‐throughput directed evolution platform for dehydrogenases. Microbeads (1.5 million per sample) carrying both variant DNA and an immobilised analogue of NAD+ were compartmentalised in water‐in‐oil emulsion droplets, together with cell‐free expression mixture and enzyme substrate, resulting in the recording of the phenotype on each bead. The beads’ phenotype could be read out and sorted for on a flow cytometer by using a highly sensitive fluorescent protein‐based sensor of the NAD+:NADH ratio. Integration of this “NAD‐display” approach with our previously described Split & Mix (SpliMLiB) method for generating large site‐saturation libraries allowed straightforward screening of fully balanced site saturation libraries of formate dehydrogenase, with diversities of 2×104. Based on modular design principles of synthetic biology NAD‐display offers access to sophisticated in vitro selections, avoiding complex technology platforms.

“NAD‐display”: Ultrahigh‐Throughput in Vitro Screening of NAD(H) Dehydrogenases Using Bead Display and Flow Cytometry

AbstractNAD(H)‐utiliing enzymes have been the subject of directed evolution campaigns to improve their function. To enable access to a larger swath of sequence space, we demonstrate the utility of a cell‐free, ultrahigh‐throughput directed evolution platform for dehydrogenases. Microbeads (1.5 million per sample) carrying both variant DNA and an immobilised analogue of NAD+ were compartmentalised in water‐in‐oil emulsion droplets, together with cell‐free expression mixture and enzyme substrate, resulting in the recording of the phenotype on each bead. The beads’ phenotype could be read out and sorted for on a flow cytometer by using a highly sensitive fluorescent protein‐based sensor of the NAD+:NADH ratio. Integration of this “NAD‐display” approach with our previously described Split & Mix (SpliMLiB) method for generating large site‐saturation libraries allowed straightforward screening of fully balanced site saturation libraries of formate dehydrogenase, with diversities of 2×104. Based on modular design principles of synthetic biology NAD‐display offers access to sophisticated in vitro selections, avoiding complex technology platforms.

How to engineer glucose oxidase for mediated electron transfer.

Glucose oxidase (GOx) is of high industrial interest for glucose sensing because of its high β-d-glucose specificity. The efficient and specific electrochemical communication between the redox center and electrodes is crucial to ensure accurate glucose determination. The efficiency of the electron transfer rates (ETR) with GOx, together with quinone diamine based mediators, is low and differs even among mediator derivatives. To design optimized enzyme-mediator couples and to describe a mediator binding model, a joint experimental and computational study was performed based on an oxygen-independent GOx variant V7 and two quinone diimine based electron mediators (QDM-1 and QDM-2), which differ in polarity and size, and ferrocenemethanol (FM). A site saturation library at position 414 was screened with all three mediators and yielded four beneficial substitutions Tyr, Met, Leu, and Val. The variants showed increased mediator activity for the more polar QDM-2 with a simultaneously decreased activity for the less polar and smaller QDM-1 and for FM. The variant GOx V7-I414Y exhibited the biggest change for the quinone diimine derivatives compared with V7 (QDM-1: 55.9 U/mg V7, 33.2 U/mg V7-I414Y; QDM-2: 2.7 U/mg V7, 12.9 U/mg V7-I414Y). TheoreticalETR calculated based on the Marcus theory were in good agreement with the experimental results. Molecular docking studies revealed a preferable binding of the two QD mediators directly in the active site, 3.5 Å away from the N5 atom of the flavin adenine dinucleotide (FAD) and in direct vicinity to position 414. In summary, position 414 in the active site was identified to modulate the electron shuttling from the FAD of the GOx to small water-soluble mediators dependent on the polarity and size of residue 414 and on the polarity and size of the mediator. The presented mediator binding model offers a promising possibility for the design of optimized enzyme-mediator couples.

New engineered phenolic biosensors based on the AraC regulatory protein.

Customized transcription factors that control gene expression in response to small molecules can act as endogenous molecular biosensors and are valuable tools for synthetic biology. We previously engineered the Escherichia coli regulatory protein AraC to respond to non-native inducers such as D-arabinose and triacetic acid lactone. Those prior studies involved the construction and screening of individual 4- or 5-site saturation mutagenesis libraries, followed by iterative rounds of positive- and negative fluorescence-activated cell sorting (FACS). Here we describe an improved screening platform and the isolation of several new and potentially useful AraC variants that respond to vanillin and salicylic acid. To increase throughput and reduce total screening time, selection steps were added to the sorting workflow. Two different site-saturation libraries and a random mutagenesis library were pooled together and >108 variants were subjected to iterative FACS and selection in search of variants responding to a panel of compounds. The new phenolic-sensing variants show responses >100-fold over background and are highly specific towards their target compound. The isolation of these variants further demonstrates the potential for engineering the AraC transcriptional regulatory protein for molecular sensing and reporting, and our improved screening system should prove effective in designing similar biosensors.

Efficient biotransformation of isoeugenol to vanillin in recombinant strains of Escherichia coli by using engineered isoeugenol monooxygenase and sol-gel chitosan membrane

Abstract Biotransformation of isoeugenol to vanillin is an attractive approach in industry for vanillin synthesis. However, the biotransformation efficiency is hampered by the low activity of isoeugenol monooxygenase (IEM) and the product inhibition towards IEM. Therefore, it is worth discovering a new IEM and improving its enzymatic properties by protein engineering. Herein, a new IEM720 gene was firstly cloned from microbial metagenome and successfully expressed in E. coli. Secondly, five potentially beneficial amino acid positions were selected by a “three-criteria” in silico system and visual inspection was performed for grouping them into two sets (regions A and B). The variant IEM720-F281Q displayed highest isoeugenol conversion and kcat value after screening the site-saturation libraries of these five positions. Furthermore, the structure-functional relationships of the beneficial substitution on IEM activity were investigated. The final vanillin concentration reached ∼4.5 g/L (∼75% conversion) by combining sol-gel chitosan membrane (for relieve of product inhibition) and E. coli cells harboring IEM720-F281Q.

Site Saturation Mutant Viruses Evolve Neutralizing Antibody Resistance in a Microfluidic Cell Culture System

Viruses constantly evolve to overcome host and environmental pressures during their replication process. These evolutionary adaptations are generally considered to be limited to a particular subset of beneficial mutations. We constructed a series of replication-competent mouse norovirus (MNV-1) site-saturation libraries to test how these viruses evolve in response to a neutralizing antibody. Residues in the A’-B’ (E296) and E’-F’ (V378, A382, D385) loops of the MNV-1 capsid protruding (P) domain were mutated using site directed mutagenesis. Viruses were evolved in the presence and absence of monoclonal neutralizing antibody A6.2 in a droplet microfluidic device such that starting virus populations at low multiplicity of infection could replicate over several passages in the absence of large population effects such as purifying selection. Droplet evolved viruses were subsequently passaged in bulk to determine whether unique combinations of mutations became fixed in the population. Using high depth RNA sequencing, we demonstrate the distribution of all possible amino acid substitutions from initial production in the laboratory, replication in the host, and passaging in the presence or absence of neutralizing antibody. These results demonstrate a method for surveying fitness landscapes that may be broadly applicable to virus-host evolution.

Inversion of Extender Unit Selectivity in the Erythromycin Polyketide Synthase by Acyltransferase Domain Engineering.

Acyltransferase (AT) domains of polyketide synthases (PKSs) select extender units for incorporation into polyketides and dictate large portions of the structures of clinically relevant natural products. Accordingly, there is significant interest in engineering the substrate specificity of PKS ATs in order to site-selectively manipulate polyketide structure. However, previous attempts to engineer ATs have yielded mutant PKSs with relaxed extender unit specificity, rather than an inversion of selectivity from one substrate to another. Here, by directly screening the extender unit selectivity of mutants from active site saturation libraries of an AT from the prototypical PKS, 6-deoxyerythronolide B synthase, a set of single amino acid substitutions was discovered that dramatically impact the selectivity of the PKS with only modest reductions of product yields. One particular substitution (Tyr189Arg) inverted the selectivity of the wild-type PKS from its natural substrate toward a non-natural alkynyl-modified extender unit while maintaining more than twice the activity of the wild-type PKS with its natural substrate. The strategy and mutations described herein form a platform for combinatorial biosynthesis of site-selectively modified polyketide analogues that are modified with non-natural and non-native chemical functionality.

Substitutions at the cofactor phosphate-binding site of a clostridial alcohol dehydrogenase lead to unexpected changes in substrate specificity.

Changing the cofactor specificity of an enzyme from nicotinamide adenine dinucleotide 2′-phosphate (NADPH) to the more abundant NADH is a common strategy for increasing overall enzyme efficiency in microbial metabolic engineering. The aim of this study was to switch the cofactor specificity of the primary–secondary alcohol dehydrogenase from Clostridium autoethanogenum, a bacterium with considerable promise for the bio-manufacturing of fuels and other petrochemicals, from strictly NADPH-dependent to NADH-dependent. We used insights from a homology model to build a site-saturation library focussed on residue S199, the position deemed most likely to disrupt binding of the 2′-phosphate of NADPH. Although the CaADH(S199X) library did not yield any NADH-dependent enzymes, it did reveal that substitutions at the cofactor phosphate-binding site can cause unanticipated changes in the substrate specificity of the enzyme. Using consensus-guided site-directed mutagenesis, we were able to create an enzyme that was stringently NADH-dependent, albeit with a concomitant reduction in activity. This study highlights the role that distal residues play in substrate specificity and the complexity of enzyme–cofactor interactions.

Enhancing RGI lyase thermostability by targeted single point mutations

Rhamnogalacturonan I lyase (RGI lyase) (EC 4.2.2.-) catalyzes the cleavage of rhamnogalacturonan I in pectins by β-elimination. In this study the thermal stability of a RGI lyase (PL 11) originating from Bacillus licheniformis DSM 13/ATCC14580 was increased by a targeted protein engineering approach involving single amino acid substitution. Nine individual amino acids were selected as targets for site-saturated mutagenesis by the use of a predictive consensus approach in combination with prediction of protein mutant stability changes and B-factor iteration testing. After extensive experimental verification of the thermal stability of the designed mutants versus the original wild-type RGI lyase, several promising single point mutations were obtained, particularly in position Glu434 on the surface of the enzyme protein. The best mutant, Glu434Leu, produced a half-life of 31min at 60°C, corresponding to a 1.6-fold improvement of the thermal stability compared to the original RGI lyase. Gly55Val was the second best mutation with a thermostability half-life increase of 27min at 60°C, and the best mutations following were Glu434Trp, Glu434Phe, and Glu434Tyr, respectively. The data verify the applicability of a combinatorial predictive approach for designing a small site saturation library for improving enzyme thermostability. In addition, new thermostable RGI lyases suitable for enzymatic upgrading of pectinaceous plant biomass materials at elevated temperatures were produced.

A Generic, Whole-Cell–Based Screening Method for Baeyer-Villiger Monooxygenases

Baeyer-Villiger monooxygenases (BVMOs) have been receiving increasing attention as enzymes useful for biocatalytic applications. Industrial requirements call for rapid and extensive redesign of these enzymes. In response to the need for screening large libraries of BVMO mutants, we established a generic screening method that allows screening of Escherichia coli cells expressing active BVMOs in 96-well plate format. For this, we first developed an expression system for production of phenylacetone monooxygenase (PAMO) in the periplasm of E. coli. This allows probing the enzyme for any target substrate while it is also compatible with extracellular coenzyme regeneration. For coenzyme regeneration, we used phosphite dehydrogenase, which forms phosphate upon NADPH recycling. This allowed the use of a chromogenic molybdate-based phosphate determination assay. The screening procedure was supplemented with a detection method for identification of mutant enzymes that act as NADPH oxidases, thereby excluding false-positives. The whole-cell–based screening method was validated by screening site–saturation libraries of PAMO and resulted in the identification of PAMO mutants with altered catalytic properties. This new method can be used for screening libraries of BVMOs for activity with any desired substrate and therefore is a powerful tool for engineering of these enzymes.

Reengineered glucose oxidase for amperometric glucose determination in diabetes analytics

Glucose oxidase is an oxidoreductase exhibiting a high β-d-glucose specificity and high stability which renders glucose oxidase well-suited for applications in diabetes care. Nevertheless, GOx activity is highly oxygen dependent which can lead to inaccuracies in amperometric β-d-glucose determinations. Therefore a directed evolution campaign with two rounds of random mutagenesis (SeSaM followed by epPCR), site saturation mutagenesis studies on individual positions, and one simultaneous site saturation library (OmniChange; 4 positions) was performed. A diabetes care well suited mediator (quinone diimine) was selected and the GOx variant (T30V I94V) served as starting point. For directed GOx evolution a microtiter plate detection system based on the quinone diimine mediator was developed and the well-known ABTS-assay was applied in microtiter plate format to validate oxygen independency of improved GOx variants. Two iterative rounds of random diversity generation and screening yielded to two subsets of amino acid positions which mainly improved activity (A173, A332) and oxygen independency (F414, V560). Simultaneous site saturation of all four positions with a reduced subset of amino acids using the OmniChange method yielded finally variant V7 with a 37-fold decreased oxygen dependency (mediator activity: 7.4U/mg WT, 47.5U/mg V7; oxygen activity: 172.3U/mg WT, 30.1U/mg V7). V7 is still highly β-d-glucose specific, highly active with the quinone diimine mediator and thermal resistance is retained (prerequisite for GOx coating of diabetes test stripes). The latter properties and V7's oxygen insensitivity make V7 a very promising candidate to replace standard GOx in diabetes care applications.

Library construction and evaluation for site saturation mutagenesis

We developed a method for creating and evaluating site-saturation libraries that consistently yields an average of 27.4±3.0 codons of the 32 possible within a pool of 95 transformants. This was verified by sequencing 95 members from 11 independent libraries within the gene encoding alkene reductase OYE 2.6 from Pichia stipitis. Correct PCR primer design as well as a variety of factors that increase transformation efficiency were critical contributors to the method's overall success. We also developed a quantitative analysis of library quality (Q-values) that defines library degeneracy. Q-values can be calculated from standard fluorescence sequencing data (capillary electropherograms) and the degeneracy predicted from an early stage of library construction (pooled plasmids from the initial transformation) closely matched that observed after ca. 1000 library members were sequenced. Based on this experience, we suggest that this analysis can be a useful guide when applying our optimized protocol to new systems, allowing one to focus only on good-quality libraries and reject substandard libraries at an early stage. This advantage is particularly important when lower-throughput screening techniques such as chiral-phase GC must be employed to identify protein variants with desirable properties, e.g., altered stereoselectivities or when multiple codons are targeted for simultaneous randomization.

Comparison of random mutagenesis and semi-rational designed libraries for improved cytochrome P450 BM3-catalyzed hydroxylation of small alkanes

Three semi-rational approaches, combinatorial site-saturation mutagenesis (CSSM) using a reduced amino acid set and two libraries based on C(orbit) and CRAM computational design algorithms targeting up to 10 active site residues, were used to engineer cytochrome P450 BM3 to demethylate dimethyl ether and hydroxylate propane and ethane. These small libraries (343-1028 variants) were all enriched with respect to the fraction functional and maximal activities compared with a random mutagenesis library and individual site-saturation libraries targeting the same residues. Despite high average amino acid substitution levels of 2.6, 5 and 7.5, the CSSM, C(orbit) and CRAM libraries had at least 75% of library members properly folded. Propane- and ethane-hydroxylating P450 BM3 variants were identified using all three mutagenesis approaches, with as few as two amino acid substitutions. The library designed using the CRAM algorithm, which sought to reduce the size of the binding pocket, produced both a higher number of active variants and variants supporting the greatest number of catalytic turnovers. The most active variant E32 supports 16 800 propane turnovers at 36% coupling, which rivals the activity of variants obtained after 10-12 rounds of directed evolution using random and site-saturation mutagenesis. None of the variants in this study achieved the complete re-specialization for propane hydroxylation (including 93% coupling) previously obtained via multiple rounds of mutagenesis and screening. However, these semi-rational approaches allowed for large jumps in sequence space to variants with the desired functions.

Focused Directed Evolution of Pentaerythritol Tetranitrate Reductase by Using Automated Anaerobic Kinetic Screening of Site‐Saturated Libraries

This work describes the development of an automated robotic platform for the rapid screening of enzyme variants generated from directed evolution studies of pentraerythritol tetranitrate (PETN) reductase, a target for industrial biocatalysis. By using a 96-well format, near pure enzyme was recovered and was suitable for high throughput kinetic assays; this enabled rapid screening for improved and new activities from libraries of enzyme variants. Initial characterisation of several single site-saturation libraries targeted at active site residues of PETN reductase, are described. Two mutants (T26S and W102F) were shown to have switched in substrate enantiopreference against substrates (E)-2-aryl-1-nitropropene and α-methyl-trans-cinnamaldehyde, respectively, with an increase in ee (62 % (R) for W102F). In addition, the detection of mutants with weak activity against α,β-unsaturated carboxylic acid substrates showed progress in the expansion of the substrate range of PETN reductase. These methods can readily be adapted for rapid evolution of enzyme variants with other oxidoreductase enzymes.

Thermostabilization of an esterase by alignment-guided focussed directed evolution

Site-saturation libraries of the Pseudomonas fluorescens esterase were created targeting three surface positions to increase its thermostability on the basis of the B-factor iterative test principle. All three positions were saturated simultaneously using our recently developed protocol for the design of 'small, but smart' mutant libraries bearing only consensus-like mutations. Hence, the library size could be significantly reduced while ensuring a high hit rate. Variants could be identified that showed significantly improved stability (8° C higher compared with the wild type) without compromising specific activity. Subsequent iterative saturation mutagenesis gave an esterase mutant with a 9° C increased melting point, but unchanged catalytic properties.

Generation and identification of variants with improved purification yield of Bowman-Birk protease inhibitors carrying protein binding loops

Replacing the chymotrypsin inhibitory loop of soybean Bowman-Birk inhibitor (sBBI) with a VEGF binding peptide (BBI-AV) significantly reduces the overall purification yield when BBI-AV is produced as a fusion protein in a Bacillus subtilis expression system. The low purification yield is primarily due to a higher fraction of molecules with incorrect disulfide bond configurations after production and also after disulfide bond shuffling induced by 2-mercaptoethanol. To improve production yields, site-saturation libraries were generated at 39 out of the 66 amino acid residues of BBI-AV. Initial screens were designed to select for variants with higher trypsin inhibitory activities than the parent after treatment with a reducing agent. Secondary screens were developed to select for variants with the highest purification yields, and to also eliminate any false positives. From the screens, it was found that positively charged substitutions in the exposed hydrophobic patch region (sites 27, 29, 40, 50 & 52) are especially productive. In fact, one substitution, F50R, improves the purification yield to nearly the same level as wild-type sBBI. Productive amino acid substitutions were combined to select for the variant with the best overall yield after purification. Several variants were obtained with higher purification yields than even sBBI. The octuple variants, A13I-S25R-M27A-L29P-S31A-A40H-F50K-V52T and A13I-S25K-M27A-L29R-S31E-A40K-F50Q-V52Q, are particularly productive having greater than a five fold increase in final purification yield over the parent.

Novel mutations in Moloney Murine Leukemia Virus reverse transcriptase increase thermostability through tighter binding to template-primer

In an effort to increase the thermostability of Moloney Murine Leukemia Virus reverse transcriptase (MMLV RT), we screened random and site-saturation libraries for variants that show increased resistance to thermal inactivation. We discovered five mutations E69K, E302R, W313F, L435G and N454K that collectively increase the half-life of MMLV RT at 55°C from less than 5 min to ∼30 min in the presence of template-primer. In addition, these mutations alter the thermal profile by increasing specific activity of the pentuple mutant (M5) over a broad range of cDNA synthesis temperatures (25–70°C). We further show that M5 generates higher cDNA yields and exhibits better RT–PCR performance compared to wild-type RT when used at high temperature to amplify RNA targets containing secondary structure. Finally, we demonstrate that M5 exhibits tighter binding (lower Km) to template-primer, which likely protects against heat inactivation.