In vivo imaging of model organisms is heavily reliant on fluorescent proteins with high intracellular brightness. Here we describe a practical method for rapid optimization of fluorescent proteins via directed molecular evolution in cultured mammalian cells. Using this method, we were able to perform screening of large gene libraries containing up to 2 × 107 independent random genes of fluorescent proteins expressed in HEK cells, completing one iteration of directed evolution in a course of 8 days. We employed this approach to develop a set of green and near-infrared fluorescent proteins with enhanced intracellular brightness. The developed near-infrared fluorescent proteins demonstrated high performance for fluorescent labeling of neurons in culture and in vivo in model organisms such as Caenorhabditis elegans, Drosophila, zebrafish, and mice. Spectral properties of the optimized near-infrared fluorescent proteins enabled crosstalk-free multicolor imaging in combination with common green and red fluorescent proteins, as well as dual-color near-infrared fluorescence imaging. The described method has a great potential to be adopted by protein engineers due to its simplicity and practicality. We also believe that the new enhanced fluorescent proteins will find wide application for in vivo multicolor imaging of small model organisms.

A prophylactic vaccine is needed to slow the spread of HIV-1 infection. Optimization of the wild-type envelope glycoproteins to create immunogens that can elicit effective neutralizing antibodies is a high priority. Starting with ten genes encoding subtype B HIV-1 gp120 envelope glycoproteins and using in vitro homologous DNA recombination, we created chimeric gp120 variants that were screened for their ability to bind neutralizing monoclonal antibodies. Hundreds of variants were identified with novel antigenic phenotypes that exhibit considerable sequence diversity. Immunization of rabbits with these gp120 variants demonstrated that the majority can induce neutralizing antibodies to HIV-1. One novel variant, called ST-008, induced significantly improved neutralizing antibody responses when assayed against a large panel of primary HIV-1 isolates. Further study of various deletion constructs of ST-008 showed that the enhanced immunogenicity results from a combination of effective DNA priming, an enhanced V3-based response, and an improved response to the constant backbone sequences.

Alteration of antibiotic binding sites through modification of ribosomal RNA (rRNA) is a common form of resistance to ribosome-targeting antibiotics. The rRNA-modifying enzyme Cfr methylates an adenosine nucleotide within the peptidyl transferase center, resulting in the C-8 methylation of A2503 (m8A2503). Acquisition of cfr results in resistance to eight classes of ribosome-targeting antibiotics. Despite the prevalence of this resistance mechanism, it is poorly understood whether and how bacteria modulate Cfr methylation to adapt to antibiotic pressure. Moreover, direct evidence for how m8A2503 alters antibiotic binding sites within the ribosome is lacking. In this study, we performed directed evolution of Cfr under antibiotic selection to generate Cfr variants that confer increased resistance by enhancing methylation of A2503 in cells. Increased rRNA methylation is achieved by improved expression and stability of Cfr through transcriptional and post-transcriptional mechanisms, which may be exploited by pathogens under antibiotic stress as suggested by natural isolates. Using a variant that achieves near-stoichiometric methylation of rRNA, we determined a 2.2 Å cryo-electron microscopy structure of the Cfr-modified ribosome. Our structure reveals the molecular basis for broad resistance to antibiotics and will inform the design of new antibiotics that overcome resistance mediated by Cfr.

Directed evolution has been used for decades to engineer biological systems at or below the organismal level. Above the organismal level, a small number of studies have attempted to artificially select microbial ecosystems, with uneven and generally modest success. Our theoretical understanding of artificial ecosystem selection is limited, particularly for large assemblages of asexual organisms, and we know little about designing efficient methods to direct their evolution. Here, we have developed a flexible modelling framework that allows us to systematically probe any arbitrary selection strategy on any arbitrary set of communities and selected functions. By artificially selecting hundreds of in silico microbial metacommunities under identical conditions, we first show that the main breeding methods used to date, which do not necessarily let communities reach their ecological equilibrium, are outperformed by a simple screen of sufficiently mature communities. We then identify a range of alternative directed evolution strategies that, particularly when applied in combination, are well suited for the top-down engineering of large, diverse and stable microbial consortia. Our results emphasize that directed evolution allows an ecological structure-function landscape to be navigated in search of dynamically stable and ecologically resilient communities with desired quantitative attributes.

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

Directed evolution of orthogonal aminoacyl-tRNA synthetases (AARSs) enables site-specific installation of noncanonical amino acids (ncAAs) into proteins. Traditional evolution techniques typically produce AARSs with greatly reduced activity and selectivity compared to their wild-type counterparts. We designed phage-assisted continuous evolution (PACE) selections to rapidly produce highly active and selective orthogonal AARSs through hundreds of generations of evolution. PACE of a chimeric Methanosarcina spp. pyrrolysyl-tRNA synthetase (PylRS) improved its enzymatic efficiency (kcat/KMtRNA) 45-fold compared to the parent enzyme. Transplantation of the evolved mutations into other PylRS-derived synthetases improved yields of proteins containing noncanonical residues up to 9.7-fold. Simultaneous positive and negative selection PACE over 48 h greatly improved the selectivity of a promiscuous Methanocaldococcus jannaschii tyrosyl-tRNA synthetase variant for site-specific incorporation of p-iodo-L-phenylalanine. These findings offer new AARSs that increase the utility of orthogonal translation systems and establish the capability of PACE to efficiently evolve orthogonal AARSs with high activity and amino acid specificity.

Cells that constitutively diversify their immunoglobulin genes can be used for selection of novel antibodies and for refining existing affinities and specificities. Here, we report an adaptation of the chicken DT40 system wherein its capacity for somatic hypermutation is harnessed to evolve human antibodies expressed as single-chain variable fragments (scFvs). Expression of membrane-anchored scFvs from within the rearranged Igλ locus created self-diversifying scFv libraries from which we could both select scFvs of a desired specificity and evolve both the specificity and affinity of existing scFvs by iterative expansion and selection. From these scFvs, we were able to create fully human IgG antibodies with nanomolar affinities. We further enhanced the functionality of the system by creating a pool of DT40 scFv lines with high levels of mutation driven by the overexpression of a hyperactive variant of activation-induced deaminase. From this library, we successfully isolated scFvs that bound the spliceosome factor CWC15 and the cytokine human IFNγ. Our results demonstrate the flexibility and utility of DT40 for rapid generation of scFv repertoires and efficient selection, evolution and affinity maturation of scFv specificities.

Adenoviruses (Ads) are the most frequently used viruses for oncolytic and gene therapy purposes. Most Ad-based vectors have been generated through rational design. Although this led to significant vector improvements, it is often hampered by an insufficient understanding of Ad's intricate functions and interactions. Here, to evade this issue, we adopted a novel, mutator Ad polymerase-based, 'accelerated-evolution' approach that can serve as general method to generate or optimize adenoviral vectors. First, we site specifically substituted Ad polymerase residues located in either the nucleotide binding pocket or the exonuclease domain. This yielded several polymerase mutants that, while fully supportive of viral replication, increased Ad's intrinsic mutation rate. Mutator activities of these mutants were revealed by performing deep sequencing on pools of replicated viruses. The strongest identified mutators carried replacements of residues implicated in ssDNA binding at the exonuclease active site. Next, we exploited these mutators to generate the genetic diversity required for directed Ad evolution. Using this new forward genetics approach, we isolated viral mutants with improved cytolytic activity. These mutants revealed a common mutation in a splice acceptor site preceding the gene for the adenovirus death protein (ADP). Accordingly, the isolated viruses showed high and untimely expression of ADP, correlating with a severe deregulation of E3 transcript splicing.

Directed evolution is a powerful method for engineering biology in the absence of detailed sequence-function relationships. To enable directed evolution of complex phenotypes encoded by multigene pathways, we require large library sizes for DNA sequences >5-10 kb in length, elimination of genomic hitchhiker mutations, and decoupling of diversification and screening steps. To meet these challenges, we developed Inducible Directed Evolution (IDE), which uses a temperate bacteriophage to package large plasmids and transfer them to naive cells after intracellular mutagenesis. To demonstrate IDE, we evolved a 5-gene pathway from Bacillus licheniformis that accelerates tagatose catabolism in Escherichia coli, resulting in clones with 65% shorter lag times during growth on tagatose after only two rounds of evolution. Next, we evolved a 15.4 kb, 10-gene pathway from Bifidobacterium breve UC2003 that aids E. coli's utilization of melezitose. After three rounds of IDE, we isolated evolved pathways that both reduced lag time by more than 2-fold and enabled 150% higher final optical density. Taken together, this work enhances the capacity and utility of a whole pathway directed evolution approach in E. coli.

Proton-pumping rhodopsins (PPRs) are photoactive retinal-binding proteins that transport ions across biological membranes in response to light. These proteins are interesting for light-harvesting applications in bioenergy production, in optogenetics applications in neuroscience, and as fluorescent sensors of membrane potential. Little is known, however, about how the protein sequence determines the considerable variation in spectral properties of PPRs from different biological niches or how to engineer these properties in a given PPR. Here we report a comprehensive study of amino acid substitutions in the retinal-binding pocket of Gloeobacter violaceus rhodopsin (GR) that tune its spectral properties. Directed evolution generated 70 GR variants with absorption maxima shifted by up to ±80nm, extending the protein's light absorption significantly beyond the range of known natural PPRs. While proton-pumping activity was disrupted in many of the spectrally shifted variants, we identified single tuning mutations that incurred blue and red shifts of 42nm and 22nm, respectively, that did not disrupt proton pumping. Blue-shifting mutations were distributed evenly along the retinal molecule while red-shifting mutations were clustered near the residue K257, which forms a covalent bond with retinal through a Schiff base linkage. Thirty eight of the identified tuning mutations are not found in known microbial rhodopsins. We discovered a subset of red-shifted GRs that exhibit high levels of fluorescence relative to the WT (wild-type) protein.

Low performance of actively targeted nanomedicines required revision of the traditional drug targeting paradigm and stimulated the development of novel phage-programmed, self-navigating drug delivery vehicles. In the proposed smart vehicles, targeting peptides, selected from phage libraries using traditional principles of affinity selection, are substituted for phage proteins discovered through combinatorial avidity selection. Here, we substantiate the potential of combinatorial avidity selection using landscape phage in the discovery of Short Linear Motifs (SLiMs) and their partner domains. We proved an algorithm for analysis of phage populations evolved through multistage screening of landscape phage libraries against the MDA-MB-231 breast cancer cell line. The suggested combinatorial avidity selection model proposes a multistage accumulation of Elementary Binding Units (EBU), or Core Motifs (CorMs), in landscape phage fusion peptides, serving as evolutionary initiators for formation of SLiMs. Combinatorial selection has the potential to harness directed molecular evolution to create novel smart materials with diverse novel, emergent properties.

Pinene, a natural active monoterpene, is widely used as a flavoring agent, perfume, medicine, and biofuel. Although genetically engineered microorganisms have successfully produced pinene, to date, the biological yield of pinene is much lower than that of semiterpenes (isoprene) and sesquiterpenes (farnesene). In addition to the low heterologous expression of geranyl pyrophosphate synthase (GPPS) and pinene synthase (PS), cytotoxicity due to accumulation of the monoterpene also limits the production of pinene in microorganisms. In this study, we attempted to use two strategies to increase the biological yield of pinene. By deleting the random coils of GPPS and PS alone or in combination, a strain with a 335% yield increase was obtained. Additionally, upon computer-guided molecular modeling and docking of GPPS with isopentenyl pyrophosphate (IPP), its substrate, the key sites located within the catalytic pocket for substrate binding, was predicted. After screening, a strain harboring the T273R mutation of GPPS was selected among a batch of mutations of the key sites with a 154% increase in pinene yield.

Using directed evolution and site-directed mutagenesis, we have isolated a highly thermostable variant of Aspergillus niger glucoamylase (GA), designated CR2-1. CR2-1 includes the previously described mutations Asn20Cys and Ala27Cys (forming a new disulfide bond), Ser30Pro, Thr62Ala, Ser119Pro, Gly137Ala, Thr290Ala, His391Tyr and Ser436Pro. In addition, CR2-1 includes several new putative thermostable mutations, Val59Ala, Val88Ile, Ser211Pro, Asp293Ala, Thr390Ser, Tyr402Phe and Glu408Lys, identified by directed evolution. CR2-1 GA has a catalytic efficiency (k(cat)/K(m)) at 35°C and a specific activity at 50°C similar to that of wild-type GA. Irreversible inactivation tests indicated that CR2-1 increases the free energy of thermoinactivation at 80°C by 10 kJ mol(-1) compared with that of wild-type GA. Thus, CR2-1 is more thermostable (by 5 kJ mol(-1) at 80°C) than the most thermostable A. niger GA variant previously described, THS8. In addition, Val59Ala and Glu408Lys were shown to individually increase the thermostability in GA variants by 1 and 2 kJ mol(-1), respectively, at 80°C.

Glycosylation patterns commonly change in cancer, resulting in expression of tumor-associated carbohydrate antigens (TACA). While promising, currently available anti-glycan antibodies are not useful for clinical cancer therapy. Here, we show that potent anti-glycan antibodies can be engineered to acquire cancer therapeutic efficacy. We designed yeast surface display to generate and select for therapeutic antibodies against the TACA SLea (CA19-9) in colon and pancreatic cancers. Elite clones showed increased affinity, better specificity, improved binding of human pancreatic and colon cancer cell lines, and increased complement-dependent therapeutic efficacy. Molecular modeling explained the structural basis for improved antibody functionality at the molecular level. These new tools of directed molecular evolution and selection for effective anti-glycan antibodies, provide insights into the mechanisms of cancer therapy targeting glycosylation, and provide major methodological advances that are likely to open up innovative avenues of research in the field of cancer theranostics.

The direct hydroxylation of benzene to hydroquinone (HQ) under mild reaction conditions is a challenging task for chemical catalysts. Cytochrome P450 (CYP) monooxygenases are known to catalyze the oxidation of a variety of aromatic compounds with atmospheric dioxygen. Protein engineering campaigns led to the identification of novel P450 variants, which yielded improvements in respect to activity, specificity, and stability. An effective screening strategy is crucial for the identification of improved enzymes with desired characteristics in large mutant libraries. Here, we report a first screening system designed for screening of P450 variants capable to produce hydroquinones. The hydroquinone quantification assay is based on the interaction of 4-nitrophenylacetonitrile (NpCN) with hydroquinones under alkaline conditions. In the 96-well plate format, a low detection limit (5 μM) and a broad linear detection range (5 to 250 μM) were obtained. The NpCN assay can be used for the quantification of dihydroxylated aromatic compounds such as hydroquinones, catechols, and benzoquinones. We chose the hydroxylation of pseudocumene by P450 BM3 as a target reaction and screened for improved trimethylhydroquinone (TMHQ) formation. The new P450 BM3 variant AW2 (R47Q, Y51F, I401M, A330P) was identified by screening a saturation mutagenesis library of amino acid position A330 with the NpCN assay. In summary, a 70-fold improved TMHQ formation was achieved with P450 BM3 AW2 when compared to the wild type (WT) and a 1.8-fold improved TMHQ formation compared to the recently reported P450 BM3 M3 (R47S, Y51W, A330F, I401M).

Ancestral sequence reconstruction and resurrection provides useful information for protein engineering, yet its alliance with directed evolution has been little explored. In this study, we have resurrected several ancestral nodes of fungal laccases dating back ∼500 to 250 million years. Unlike modern laccases, the resurrected Mesozoic laccases were readily secreted by yeast, with similar kinetic parameters, a broader stability, and distinct pH activity profiles. The resurrected Agaricomycetes laccase carried 136 ancestral mutations, a molecular testimony to its origin, and it was subjected to directed evolution in order to improve the rate of 1,3-cyclopentanedione oxidation, a β-diketone initiator commonly used in vinyl polymerization reactions.IMPORTANCE The broad variety of biotechnological uses of fungal laccases is beyond doubt (food, textiles, pulp and paper, pharma, biofuels, cosmetics, and bioremediation), and protein engineering (in particular, directed evolution) has become the key driver for adaptation of these enzymes to harsh industrial conditions. Usually, the first requirement for directed laccase evolution is heterologous expression, which presents an important hurdle and often a time-consuming process. In this work, we resurrected a fungal Mesozoic laccase node which showed strikingly high heterologous expression and pH stability. As a proof of concept that the ancestral laccase is a suitable blueprint for engineering, we performed a quick directed evolution campaign geared to the oxidation of the β-diketone 1,3-cyclopentanedione, a poor laccase substrate that is used in the polymerization of vinyl monomers.

Metallohydrolases are ubiquitous in nearly all subclasses of hydrolases, utilizing metal elements to activate a water molecule and facilitate its subsequent dissociation of diverse chemical bonds. However, such a catalytic role of metal ions is rarely found with glycosidases that hydrolyze the glycosidic bonds in sugars. Herein, we design metalloglycosidases by constructing a hydrolytically active Zn-binding site within a barrel-shaped outer membrane protein OmpF. Structure- and mechanism-based redesign and directed evolution have led to the emergence of Zn-dependent glycosidases with catalytic proficiency of 2.8 × 109 and high β-stereoselectivity. Biochemical characterizations suggest that the Zn-binding site constitutes a key catalytic motif along with at least one adjacent acidic residue. This work demonstrates that unprecedented metalloenzymes can be tailor-made, expanding the scope of inorganic reactivities in proteinaceous environments, resetting the structural and functional diversity of metalloenzymes, and providing the potential molecular basis of unidentified metallohydrolases and novel whole-cell biocatalysts.

Microbe-associated molecular patterns (MAMPs) are molecules, or domains within molecules, that are conserved across microbial taxa and can be recognized by a plant or animal immune system. Although MAMP receptors have evolved to recognize conserved epitopes, the MAMPs in some microbial species or strains have diverged sufficiently to render them unrecognizable by some host immune systems. In this study, we carried out in vitro evolution of the Arabidopsis thaliana flagellin receptor FLAGELLIN-SENSING 2 (FLS2) to isolate derivatives that recognize one or more flagellin peptides from bacteria for which the wild-type Arabidopsis FLS2 confers little or no response. A targeted approach generated amino acid variation at FLS2 residues in a region previously implicated in flagellin recognition. The primary screen tested for elevated response to the canonical flagellin peptide from Pseudomonas aeruginosa, flg22. From this pool, we then identified five alleles of FLS2 that confer modest (quantitatively partial) recognition of an Erwinia amylovora flagellin peptide. Use of this Erwinia-based flagellin peptide to stimulate Arabidopsis plants expressing the resulting FLS2 alleles did not lead to a detectable reduction of virulent P. syringae pv. tomato growth. However, combination of two identified mutations into a single allele further increased FLS2-mediated responses to the E. amylovora flagellin peptide. These studies demonstrate the potential to raise the sensitivity of MAMP receptors toward particular targets.

Sulfide-dependent THI4 thiazole synthases could potentially be used to replace plant cysteine-dependent suicide THI4s, whose high protein turnover rates make thiamin synthesis exceptionally energy-expensive. However, sulfide-dependent THI4s are anaerobic or microoxic enzymes and hence unadapted to the aerobic conditions in plants; they are also slow enzymes (kcat < 1 h-1). To improve aerotolerance and activity, we applied continuous directed evolution under aerobic conditions in the yeast OrthoRep system to two sulfide-dependent bacterial THI4s. Seven beneficial single mutations were identified, of which five lie in the active-site cleft predicted by structural modeling and two recapitulate features of naturally aerotolerant THI4s. That single mutations gave substantial improvements suggests that further advance under selection will be possible by stacking mutations. This proof-of-concept study established that the performance of sulfide-dependent THI4s in aerobic conditions is evolvable and, more generally, that yeast OrthoRep provides a plant-like bridge to adapt nonplant enzymes to work better in plants.

Tobacco etch virus protease (TEV) is one of the most widely used proteases in biotechnology because of its exquisite sequence specificity. A limitation, however, is its slow catalytic rate. We developed a generalizable yeast-based platform for directed evolution of protease catalytic properties. Protease activity is read out via proteolytic release of a membrane-anchored transcription factor, and we temporally regulate access to TEV's cleavage substrate using a photosensory LOV domain. By gradually decreasing light exposure time, we enriched faster variants of TEV over multiple rounds of selection. Our TEV-S153N mutant (uTEV1Δ), when incorporated into the calcium integrator FLARE, improved the signal/background ratio by 27-fold, and enabled recording of neuronal activity in culture with 60-s temporal resolution. Given the widespread use of TEV in biotechnology, both our evolved TEV mutants and the directed-evolution platform used to generate them could be beneficial across a wide range of applications.