Engineering aggregation-resistant proteins by directed evolution

A method was recently described for selecting aggregation-resistant antibody domains. A repertoire of domains displayed on filamentous bacteriophage were heated/cooled and selected for binding to the affinity ligand protein A specific for the folded domains. Here we describe a generalization of the method based on the selection for retained phage infectivity and for the binding of an appended sequence tag, and applicable to any protein displayable in multivalent form on phage.     

Biased mutation-assembling: an efficient method for rapid directed evolution through simultaneous mutation accumulation

We have developed an efficient optimization technique, 'biased mutation-assembling', for improving protein properties such as thermostability. In this strategy, a mutant library is constructed using the overlap extension polymerase chain reaction technique with DNA fragments from wild-type and phenotypically advantageous mutant genes, in which the number of mutations assembled in the wild-type gene is stochastically controlled by the mixing ratio of the mutant DNA fragments to wild-type fragments. A high mixing ratio results in a mutant composition biased to favor multiple-point mutants. We applied this strategy to improve the thermostability of prolyl endopeptidase from Flavobacterium meningosepticum as a case study and found that the proportion of thermostable mutants in a library increased as the mixing ratio was increased. If the proportion of thermostable mutants increases, the screening effort needed to find them should be reduced. Indeed, we isolated a mutant with a 1200-fold longer activity half-life at 60 degrees C than that of wild-type prolyl endopeptidase after screening only 2000 mutants from a library prepared with a high mixing ratio. Our results indicate that an aggressive accumulation of advantageous mutations leads to an increase in the quality of the mutant library and a reduction in the screening effort required to find superior mutants.     

Engineering mammalian cytochrome P450 2B1 by directed evolution for enhanced catalytic tolerance to temperature and dimethyl sulfoxide

The previously laboratory-evolved cytochrome P450 2B1 quadruple mutant V183L/F202L/L209A/S334P (QM), which showed enhanced H(2)O(2)-mediated substrate oxidation, has now been shown to exhibit a >3.0-fold decrease in K(m,HOOH) for 7-ethoxy-4-trifluoromethylcoumarin (7-EFC) O-deethylation compared with the parental enzyme L209A. Subsequently, a streamlined random mutagenesis and a high-throughput screening method were developed using QM to screen and select mutants with enhanced tolerance of catalytic activity to temperature and dimethyl sulfoxide (DMSO). Upon screening >3000 colonies, we identified QM/L295H and QM/K236I/D257N with enhanced catalytic tolerance to temperature and DMSO. QM/L295H exhibited higher activity than QM at a broad range of temperatures (35-55 degrees C) and maintained approximately 1.4-fold higher activity than QM at 45 degrees C for 6 h. In addition, QM/L295H showed a significant increase in T(m,app) compared with L209A. QM/L295H and QM/K236I/D257N exhibited higher activity than QM at a broad range of DMSO concentrations (2.5-15%). Furthermore, QM/K236I/D257N/L295H was constructed by combining QM/K236I/D257N with L295H using site-directed mutagenesis and exhibited a >2-fold higher activity than QM at nearly the entire range of DMSO concentrations. In conclusion, in addition to engineering mammalian cytochromes P450 for enhanced activity, directed evolution can also be used to optimize catalytic tolerance to temperature and organic solvent.     

Digital screening methodology for the directed evolution of transglycosidases

Engineering of glycosidases with efficient transglycosidasesactivity is an alternative to glycosyltransferases or glycosynthasesfor the synthesis of oligosaccharides and glycoconjugates. However,the engineering of transglycosidases by directed evolution methodologiesis hampered by the lack of efficient screening systems for sugar-transferactivity. We report here the development of digital imaging-basedhigh-throughput screening methodology for the directed evolutionof glycosidases into transgalactosidases. Using this methodology,we detected transglycosidase mutants in intact Escherichia colicells by digital imaging monitoring of the activation of non-or low-hydrolytic mutants by an acceptor substrate. We screenedseveral libraries of mutants of β-glycosidase from Thermusthermophilus using this methodology and found variants withup to a 70-fold overall increase in the transglycosidase/hydrolysisactivity ratio. Using natural disaccharide acceptors, thesetransglycosidase mutants were able to synthesise trisaccharides,as a mixture of two regioisomers, with up to 76% yield.     

Directed evolution on the cold adapted properties of TAB5 alkaline phosphatase

Psychrophilic alkaline phosphatase (AP) from the Antarctic strain TAB5 was subjected to directed evolution in order to identify the key residues steering the enzyme's cold-adapted activity and stability. A round of random mutagenesis and further recombination yielded three thermostable and six thermolabile variants of the TAB5 AP. All of the isolated variants were characterised by their residual activity after heat treatment, Michaelis-Menten kinetics, activation energy and microcalorimetric parameters of unfolding. In addition, they were modelled into the structure of the TAB5 AP. Mutations which affected the cold-adapted properties of the enzyme were all located close to the active site. The destabilised variants H135E and H135E/G149D had 2- and 3-fold higher kcat, respectively, than the wild-type enzyme. Wild-type AP has a complex heat-induced unfolding pattern while the mutated enzymes loose local unfolding transitions and have large shifts of the Tm values. Comparison of the wild-type and mutated TAB5 APs demonstrates that there is a delicate balance between the enzyme activity and stability and that it is possible to improve the activity and thermostability simultaneously as demonstrated in the case of the H135E/G149D variant compared to H135E.     

Combination of computational prescreening and experimental library construction can accelerate enzyme optimization by directed evolution

Chiral compounds can be produced efficiently by using biocatalysts. However, wild-type enzymes often do not meet the requirements of a production process, making optimization by rational design or directed evolution necessary. Here, we studied the lipase-catalyzed hydrolysis of the model substrate 1-(2-naphthyl)ethyl acetate both theoretically and experimentally. We found that a computational equivalent of alanine scanning mutagenesis based on QM/MM methodology can be applied to identify amino acid positions important for the activity of the enzyme. The theoretical results are consistent with concomitant experimental work using complete saturation mutagenesis and high-throughput screening of the target biocatalyst, a lipase from Bacillus subtilis. Both QM/MM-based calculations and molecular biology experiments identify histidine 76 as a residue that strongly affects the catalytic activity. The experiments demonstrate its important influence on enantioselectivity.     

Directed evolution by accumulating tailored mutations: thermostabilization of lactate oxidase with less trade-off with catalytic activity

We assumed that adverse effects posed by introducing multiple mutations could be decomposed into those of each of the component mutations and that the risk could be reduced by the accumulation of mutations that were finely tuned for directed improvement of a specific property. We propose here a directed evolution strategy for improving a specific property with less effect on other ones. This strategy is composed of fine-tuning of mutations and their accumulation by our original mutation-assembling method. In this study, we selected lactate oxidase (LOX) as a model enzyme, because its directed evolution had showed a trade-off between thermostability and catalytic activity. Mutation profiling at each of the sites found by error-prone PCR revealed a strong inverse relationship between the two properties. Thermostable mutations with less effect on catalytic activity were selected at each site and accumulated with ideal combinations by our method. The resultant multiple mutants exhibited 5- to 10-fold superior catalytic activity and comparable thermostability with those created by accumulating thermostable mutations, which were not tuned for catalytic activity. This result demonstrates that the accumulation of fine-tuned mutations is an advantageous approach to reduce the risk of adverse effects posed by accumulating multiple mutations.     

Laboratory evolution of P450 BM3 for mediated electron transfer yielding an activity-improved and reductase-independent variant

One of the main obstacles in employing P450 monooxygenases for preparative chemical syntheses in cell-free systems is their requirement for cofactors such as NAD(P)H. In order to engineer P450 BM3 from Bacillus megaterium for cost-effective process conditions in vitro, a validated medium throughput screening system based on cheap Zn dust as an electron source and Cobalt(III)sepulchrate (Co(III)sep) as a mediator was reported. In the current study, the alternative cofactor system Zn/Co(III)sep was used in a directed evolution experiment to improve the Co(III)sep-mediated electron transfer to P450 BM3. A variant, carrying five mutations (R47F F87A V281G M354S D363H, Table I), P450 BM3 M5 was identified and characterized with respect to its kinetic parameters. P450 BM3 M5 achieved for mediated electron transfer a 2-fold higher k(cat) value and a 3-fold higher catalytic efficiency compared with the starting point mutant P450 BM3 F87A (k(cat): 62 min(-1) compared with 28 min(-1); k(cat)/K(m): 62 microM(-1)min(-1) compared to 19 microM(-1)min(-1)). For obtaining first insights on electron transfer contributions, three reductase-deficient variants were generated. The reductase-deficient mutant of P450 BMP M5 exhibited a catalytic efficiency of 69% and a k(cat) value of 89% of the values obtained for P450 BM3 M5.     

Modulating D-amino acid oxidase substrate specificity: production of an enzyme for analytical determination of all D-amino acids by directed evolution

Recent research on the flavoenzyme D-amino acid oxidase from Rhodotorula gracilis (RgDAAO) has revealed new, intriguing properties of this catalyst and offers novel biotechnological applications. Among them, the reaction of RgDAAO has been exploited in the analytical determination of the D-amino acid content in biological samples. However, because the enzyme does not oxidize acidic D-amino acids, it cannot be used to detect the total amount of D-amino acids. We now present the results obtained using a random mutagenesis approach to produce RgDAAO mutants with a broader substrate specificity. The libraries of RgDAAO mutants were generated by error-prone PCR, expressed in BL21(DE3)pLysS Escherichia coli cells and screened for their ability to oxidize different substrates by means of an activity assay. Five random mutants that have a 'modified' substrate specificity, more useful for the analytical determination of the entire content of D-amino acids than wild-type RgDAAO, have been isolated. With the only exception of Y223 and G199, none of the effective amino acid substitutions lie in segments predicted to interact directly with the bound substrate. The substitutions appear to cluster on the protein surface: it would not have been possible to predict that these substitutions would enhance DAAO activity. We can only conclude that these substitutions synergistically generate small structural changes that affect the dynamics and/or stability of the protein in a way that enhances substrate binding or subsequently catalytic turnover.     

Engineering of Pseudomonas aeruginosa lipase by directed evolution for enhanced amidase activity: mechanistic implication for amide hydrolysis by serine hydrolases

A lipase from Pseudomonas aeruginosa was subjected to directed evolution for increased amidase activity to probe the catalytic mechanism of serine hydrolases for the hydrolysis of amides. Random mutagenesis combined with saturation mutagenesis for all the amino acid residues at the substrate-binding site successfully identified the mutation at the residue 252 next to the catalytic H251 as a hot spot for selectively increasing the amidase activity of the lipase. The saturation mutagenesis targeted for the oxyanion hole (M16 and H83) gave no positive results. The substitutions of Met or Phe for Leu252 significantly increased the amidase activity toward N-(2-naphthyl)oleamide (2), whereas the esterase activity toward structurally similar 2-naphthyl oleate (1) was not affected by the substitution. The triple mutant F207S/A213D/M252F (Sat252) exhibited amidase activity (k(cat)/K(m)) 28-fold higher than that of the wild-type lipase. Kinetic analysis of Sat252 and its parental clone 10F12 revealed that the amidase activity was increased by the increase in the catalytic efficiency (k(cat)). The increase in k(cat) suggested the importance of the leaving group protonation by the catalytic His during the break down of the tetrahedral intermediate in the hydrolysis of amides.     

Construction of rat aldolase C expression plasmid and the hybrid formation between rat aldolase C and human aldolase A or B co-expressed in Escherichia coli

Rat aldolase C cDNA was inserted in an Escherichia coli expressionvector to construct the rat aldolase C expression plasmid, pRAC42.This plasmid produces active rat aldolase C in the transfectedE.coli host cells. The characteristics of the purified enzyme,e.g. mol. wt, electrophoretic mobilities and kinetic parameters,are indistinguishable from those of authentic rat brain aldolaseC. Three different tetrameric hybrid forms, C₃A, C₂A₂ and CA₃,in addition to C₄ and A₄, were found to be produced in the hostcell when E.coli was co-transfected with expression plasmidsfor rat aldolase C and for human aldolase A. Similarly, thehybrid forms, C₃B, C₂B₂ and CB₃, in addition to C₄ and B₄, werealso produced in the cells when co-transfected with the plasmidsfor rat aldolase C and for human aldolase B.     

Alteration of product specificity of Aeropyrum pernix farnesylgeranyl diphosphate synthase (Fgs) by directed evolution

Directed evolution of the C25 farnesylgeranyl diphosphate synthase of Aeropyrum pernix (Fgs) was carried out by error-prone PCR with an in vivo color complementation screen utilizing carotenoid biosynthetic pathway enzymes. Screening yielded 12 evolved clones with C20 geranylgeranyl diphosphate synthase activity which were isolated and characterized in order to understand better the chain elongation mechanism of this enzyme. Analysis of these mutants revealed three different mechanisms of product chain length specificity. Two mutants (A64T and A64V) have a single mutation at the 8th amino acid upstream of a conserved first aspartate-rich motif (FARM), which is involved in the mechanism for chain elongation reaction of all prenyl diphosphate synthases. One mutant (A135T) carries a single mutation at the 7th amino acid upstream of another conserved region (141GQ142), which was recently found to be another important region controlling chain elongation of a type III C20 geranylgeranyl diphosphate synthase and Escherichia coli C15 farnesyl diphosphate synthase. Finally, one mutant carrying four mutations (V84I, H88R, I177 M and M191V) is of interest. Molecular modeling, site-directed mutagenesis and in vitro assays of this mutant suggest that product chain-length distribution can be also controlled by a structural change provoked by a cooperative interaction of amino acids.     

Development of a screening platform for directed evolution using the reef coral fluorescent protein ZsGreen as a solubility reporter

Soluble proteins, with high expression levels, are preferred candidates for structural and functional studies. In cases of low expression, aggregation or inclusion body formation, time-consuming searches for optimal expression or refolding conditions are required. We have developed a high-throughput solubility engineering and screening platform for proteins that are expressed in an insoluble form in Escherichia coli with the aim of obtaining a broad spectrum of best hits with increased solubility in difficult to express target proteins. This process has been developed using error-prone PCR to introduce random base changes in genes of interest. Expression of mutated proteins in fusion with the reef coral fluorescent protein ZsGreen as a solubility marker has enabled the selection of more soluble variants. We have used a colony picker to achieve high-throughput selection of E.coli expressing more soluble target protein-ZsGreen fusions, with increased fluorescence. The whole process enables us to complete one round of mutation, screening and analysis of 20,000 potential soluble clones within approximately 8 weeks. We describe the development of the methods using different model proteins and show one example, the kinase domain from the human EphB2 receptor, as a successful application of the whole platform.     

Directed evolution converts subtilisin E into a functional equivalent of thermitase

We used directed evolution to convert Bacillus subtilis subtilisinE into an enzyme functionally equivalent to its thermophilichomolog thermitase from Thermoactinomyces vulgaris. Five generationsof random mutagenesis, recombination and screening created subtilisinE 5-3H5, whose half-life at 83°C (3.5 min) and temperatureoptimum for activity (T_opt, 76°C) are identical with thoseof thermitase. The T_opt of the evolved enzyme is 17°C higherand its half-life at 65°C is >200 times that of wild-typesubtilisin E. In addition, 5-3H5 is more active towards thehydrolysis of succinyl-Ala-Ala-Pro-Phe-p-nitroanilide than wild-typeat all temperatures from 10 to 90°C. Thermitase differsfrom subtilisin E at 157 amino acid positions. However, onlyeight amino acid substitutions were sufficient to convert subtilisinE into an enzyme equally thermostable. The eight substitutions,which include known stabilizing mutations (N218S, N76D) andalso several not previously reported, are distributed over thesurface of the enzyme. Only two (N218S, N181D) are found inthermitase. Directed evolution provides a powerful tool to unveilmechanisms of thermal adaptation and is an effective and efficientapproach to increasing thermostability without compromisingenzyme activity.     

Concurrent mutations in six amino acids in beta-glucuronidase improve its thermostability

To achieve a thermostable beta-glucuronidase (GUS) and identify key mutation sites, we applied in vitro directed evolution strategy through DNA shuffling and obtained a highly thermostable mutant GUS gene, gus-tr, after four rounds of DNA shuffling and screening. This variant had mutations in 15 nucleic acid sites, resulting in changes in 12 amino acids (AAs). Using gus-tr as the template, we further performed site-directed mutagenesis to reverse the individual mutation to the wild-type protein. We found that six sites (Q493R, T509A, M532T, N550S, G559S and N566S) present in GUS-TR3337, were the key AAs needed to confer its high thermostability. Of these, Q493R and T509A were not reported previously as important residues for thermostability of GUS. Furthermore, all of these six mutations must be present concurrently to confer the high thermostability. We expressed the gus-tr3337 gene and purified the GUS-TR3337 protein that contained the six AA mutations. Compared with the wild-type protein which lost its activity completely after 10 min at 70 degrees C, the mutant GUS-TR3337 protein retained 75% of its activity when heated at 80 degrees C for 10 min. The GUS-TR3337 exhibited high activity even heated at 100 degrees C for 30 min on nitrocellulose filter. The comparison of molecular models of the mutated and wild-type enzyme revealed the relation of protein function and these structural modifications.     

Directed evolution of Tk-subtilisin from a hyperthermophilic archaeon: identification of a single amino acid substitution responsible for low-temperature adaptation

Tk-subtilisin from the hyperthermophilic archaeon Thermococcus kodakaraensis is synthesized in a prepro-form (prepro-Tk-subtilisin), secreted in a pro-form (pro-Tk-subtilisin), and matured to an active form (mat-Tk-subtilisin*; a Ca(2+)-bound active form of matured domain of Tk-subtilisin) upon autoprocessing and degradation of the propeptide [Tk-propeptide; propeptide of Tk-subtilisin (Gly1-Leu69)]. Pro-Tk-subtilisin exhibited halo-forming activity only at 80 degrees C, but not at 70 and 60 degrees C, because Tk-propeptide is not effectively degraded by mat-Tk-subtilisin* and forms an inactive complex with mat-Tk-subtilisin* at <80 degrees C. Random mutagenesis in the entire prepro-Tk-subtilisin gene, followed by screening for mutant proteins with halo-forming activity at 70 and 60 degrees C, allowed us to identify single Gly56 --> Ser mutation in the propeptide region responsible for low-temperature adaptation of pro-Tk-subtilisin. SDS-PAGE analyses and mat-Tk-subtilisin* activity assay of pro-G56S-subtilisin indicated more rapid maturation than pro-Tk-subtilisin. The resultant active form was indistinguishable from mat-Tk-subtilisin* in activity and stability, indicating that Gly56 --> Ser mutation does not seriously affect the folding of the mature domain. However, this mutation greatly destabilized the propeptide, making it unstructured in an isolated form. As a result, Tk-propeptide with Gly56 --> Ser mutation (G56S-propeptide) was more susceptible to proteolytic degradation and less effectively inhibited mat-Tk-subtilisin* activity than Tk-propeptide. These results suggest that pro-G56S-subtilisin is more effectively matured than pro-Tk-subtilisin at lower temperatures, because autoprocessed G56S-propeptide is unstructured upon dissociation from mat-Tk-subtilisin* and is therefore effectively degraded by mat-Tk-subtilisin*.     

Directed evolution of a type I antifreeze protein expressed in Escherichia coli with sodium chloride as selective pressure and its effect on antifreeze tolerance

Both freezing tolerance and NaCI tolerance are improved whenantifreeze proteins are expressed as fusion proteins with twodomains of staphylococcal protein A (SPA) in Escherichia coli.To characterize these properties further we created a randomlymutated expression library in E.coli, based on the winter flounderantifreeze protein HPLC-8 component gene. Low-fidelity PCR productsof this gene were fused to the spa gene encoding two domainsof the SPA. The library was screened for enhanced NaCl toleranceand four clones were selected. The freezing tolerance of eachof the selected clones was enhanced to varying extents. DNAsequencing of the isolated mutants revealed that the amphiphilicproperties of the native antifreeze protein were essentiallyconserved. Furthermore, by studying the primary sequence ofthe randomly mutated clones, in comparison with the degree offreezing tolerance, we have identified clues which help in understandingthe relationship between salt and freezing tolerance.     

Directed evolution of subtilisin E in Bacillus subtilis to enhance total activity in aqueous dimethylformamide

Sequential rounds of error-prone PCR to introduce random mutationsand screenrng of the resultant mutant libraries have been usedto enhance the total catalytic activity of subtilisin E significantlyin a non-natural environment, aqueous dimethylformamide (DMF).Seven DNA substitutions coding for three new amino acid substitutionswere identified in a mutant isolated after two additional generationsof directed evolution carried out on 10M subtilisin E, previously‘evolved’ to increase its specific activity in DMF.A Bacillus subtilis-Escherichia coli shuttle vector was developedin order to increase the size of the mutant library that couldbe established in B.subtilis and the stringency of the screeningprocess was increased to reflect total as well as specific activity.This directed evolution approach has been extremely effectivefor improving enzyme activity in a non-natural environment:the resulting-evolved 13M subtilisin exhibits specific catalyticefficiency towards the hydrolysis of a peptide substrate succlnyl-Ala-Ala-Pro-Phe-p-nitroanilidein 60% DMF solution that is three times that of the parent 10Mand 471 times that of wild type subtilisin E. The total activityof the 13M culture supernatant is enhanced 16-fold over thatof the parent 10M.     

Improvement in the alkaline stability of subtilisin using an efficient random mutagenesis and screening procedure

An efficient random mutagenesis procedure coupled to a replicaplate screen facilitated the isolation of mutant subtilisinsfrom Bacillus amyloliquefaciens that had altered autolytic stabilityunder alkaline conditions. Out of about 4000 clones screened,approximately 70 produced subtilisins with reduced stability(negatives). Two dones produced a more stable subtilisin (positives)and were identified as having a single mutation, either IIe107Valor Lys2l3Arg (the wild-type amino acid is followed by the codonposition and the mutant amino acid). One of the negative mutants,Met50Val, was at a site where other homologous subtilisins containeda Phe. When the Met50Phe mutation was introduced into the B.amyloliquefaciens gene, the mutant subtilisin was more alkalinestable. The double mutant IIe107Val/Lys2l3Arg) was more stablethan the isolated single mutant parents. The triple mutant (Met50Phe/IIel07Val/Lys2l3Arg)was even more stable than IIe107Val/Lys2l3Arg (up to two timesthe autolytic half-time of wild-type at pH 12). These studiesdemonstrate the feasibility for improving the alkaline stabilityof proteins by random mutagenesis and identifying potentialsites where substitutions from homologous proteins can improvealkaline stability.     

The stability and unfolding of an IgG binding protein based upon the B domain of protein A from Staphylococcus aureus probed by tryptophan substitution and fluorescence spectroscopy

The stability and unfolding of an immunoglobulin (Ig) G bindingprotein based upon the B domain of protein A (SpAB) from Staphylococcusaureus were studied by substituting tryptophan residues at strategiclocations within each of the three a-helical regions (al-a3)of the domain. The role of the C-terminal helix, a3, was investigatedby generating two protein constructs, one corresponding to thecomplete SpAB, the other lacking a part of ct3; the Trp substitutionswere made in both one-and two-domain versions of each of theseconstructs. The fluorescence properties of each of the single-tryptophanmutants were studied in the native state and as a function ofguanidine-HCl-mediated unfolding, and their IgG binding activitieswere determined by a competitive enzyme-linked immunosorbentassay. The free energies of folding and of binding to IgG foreach mutant were compared with those for the native domains.The effect of each substitution upon the overall structure andupon the IgG binding interface was modelled by molecular graphicsand energy minimization. These studies indicate that (i) 3 contributesto the overall stability of the domain and to the formationof the IgG binding site in l and 2, and (ii) al unfolds first,followed by 2 and 3 together.