Directed evolution for soluble and active periplasmic expression of bovine enterokinase in Escherichia coli

Bovine enterokinase light chain (EK_L) is an industrially useful protease for accurate removal of affinity-purification tags from high-value biopharmaceuticals. However, recombinant expression in Escherichia coli produces insoluble inclusion bodies, requiring solubilisation, refolding, and autocatalytic activation to recover functional enzyme. Error-prone PCR and DNA shuffling of the EK_L gene, T7 promoter, lac operon, ribosome binding site, and pelB leader sequence, yielded 321 unique variants after screening ~ 6500 colonies. The best variants had > 11,000-fold increased total activity in lysates, producing soluble enzyme that no longer needed refolding. Further characterisation identified the factors that improved total activity from an inactive and insoluble starting point. Stability was a major factor, whereby melting temperatures > 48.4 °C enabled good expression at 37 °C. Variants generally did not alter catalytic efficiency as measured by k_cat/K_m, which improved for only one variant. Codon optimisation improved the total activity in lysates produced at 37 °C. However, non-optimised codons and expression at 30 °C gave the highest activity through improved protein quality, with increased k_cat and T_m values. The 321 variants were statistically analysed and mapped to protein structure. Mutations detrimental to total activity and stability clustered around the active site. By contrast, variants with increased total activity tended to combine stabilising mutations that did not disrupt the active site.

Investigation of Active Site Amino Acid Influence on Carbon and Chlorine Isotope Fractionation during Reductive Dechlorination

Reductive dehalogenases (RDases) are corrinoid-dependent enzymes that reductively dehalogenate organohalides in respiratory processes. By comparing isotope effects in biotically-catalyzed reactions to reference experiments with abiotic corrinoid-catalysts, compound-specific isotope analysis (CSIA) has been shown to yield valuable insights into enzyme mechanisms and kinetics, including RDases. Here, we report isotopic fractionation (ε) during biotransformation of chloroform (CF) for carbon (εC = -1.52 ± 0.34‰) and chlorine (εCl = -1.84 ± 0.19‰), corresponding to a ΛC/Cl value of 1.13 ± 0.35. These results are highly suppressed compared to isotope effects observed both during CF biotransformation by another organism with a highly similar RDase (> 95% sequence identity) at the amino acid level, and to those observed during abiotic dehalogenation of CF. Amino acid differences occur at four locations within the two different RDases' active sites, and this study examines whether these differences potentially affect the observed εC, εCl, and ΛC/Cl. Structural protein models approximating the locations of the residues elucidate possible controls on reaction mechanisms and/or substrate binding efficiency. These four locations are not conserved among other chloroalkane reducing RDases with high amino acid similarity (> 90%), suggesting that these locations may be important in determining isotope fractionation within this homologous group of RDases.

Improved thermostability of proteinase K and recognizing the synergistic effect of Rosetta and FoldX approaches

Proteinase K (PRK) is a proteolytic enzyme that has been widely used in industrial applications. However, poor stability has severely limited the uses of PRK. In this work, we used two structure-guided rational design methods, Rosetta and FoldX, to modify PRK thermostability. Fifty-two single amino acid conversion mutants were constructed based on software predictions of residues that could affect protein stability. Experimental characterization revealed that 46% (21 mutants) exhibited enhanced thermostability. The top four variants, D260V, T4Y, S216Q, and S219Q, showed improved half-lives at 69°C by 12.4-, 2.6-, 2.3-, and 2.2-fold that of the parent enzyme, respectively. We also found that selecting mutations predicted by both methods could increase the predictive accuracy over that of either method alone, with 73% of the shared predicted mutations resulting in higher thermostability. In addition to providing promising new variants of PRK in industrial applications, our findings also show that combining these programs may synergistically improve their predictive accuracy.

Improved thermostability of creatinase from Alcaligenes Faecalis through non-biased phylogenetic consensus-guided mutagenesis

Background

Enzymatic quantification of creatinine has become an essential method for clinical evaluation of renal function. Although creatinase (CR) is frequently used for this purpose, its poor thermostability severely limits industrial applications. Herein, we report a novel creatinase from Alcaligenes faecalis (afCR) with higher catalytic activity and lower K_M value, than currently used creatinases. Furthermore, we developed a non-biased phylogenetic consensus method to improve the thermostability of afCR.

Results

We applied a non-biased phylogenetic consensus method to identify 59 candidate consensus residues from 24 creatinase family homologs for screening afCR mutants with improved thermostability. Twenty-one amino acids of afCR were selected to mutagenesis and 11 of them exhibited improved thermostability compared to the parent enzyme (afCR-M0). Combination of single-site mutations in sequential screens resulted in a quadruple mutant D17V/T199S/L6P/T251C (M4-2) which showed ~ 1700-fold enhanced half-life at 57 °C and a 4.2 °C higher T₅₀¹⁵ than that of afCR-M0. The mutant retained catalytic activity equivalent to afCR-M0, and thus showed strong promise for application in creatinine detection. Structural homology modeling revealed a wide range of potential molecular interactions associated with individual mutations that contributed to improving afCR thermostability.

Conclusions

Results of this study clearly demonstrated that the non-biased-phylogenetic consensus design for improvement of thermostability in afCR is effective and promising in improving the thermostability of more enzymes.

A Recombinant β-Mannanase from Thermoanaerobacterium aotearoense SCUT27: Biochemical Characterization and Its Thermostability Improvement

β-Mannanase was expressed in Thermoanaerobacterium aotearoense SCUT27 induced by locust bean gum (LBG). The open reading frame encoding a GH26 β-mannanase was identified and encoded a preprotein of 515 amino acids with a putative signal peptide. The enzyme without a signal sequence (Man25) was overexpressed in Escherichia coli with a specific activity of 1286.2 U/mg. Moreover, a facile method for β-mannanase activity screening was established based on agar plates. The optimum temperature for the purified Man25 using LBG as a substrate was 55 °C. The catalytic activity and thermostability of Man25 displayed a strong dependence on calcium ions. Through saturation mutagenesis at the putative Ca²⁺ binding sites in Man25, the best mutant ManM3-3 (D143A) presented improvements in thermostability with 3.6-fold extended half-life at 55 °C compared with that of the wild-type. The results suggest that mutagenesis at metal binding sites could be an efficient approach to increase enzyme thermostability.

Improvement of BsAPA Aspartic Protease Thermostability via Autocatalysis-Resistant Mutation

An aspartic protease gene (Bsapa) was cloned from Bispora sp. MEY-1 and expressed in Pichia pastoris. The recombinant BsAPA showed maximal activity at pH 3.0 and 75 °C and remained stable at 70 °C and below, indicating the thermostable nature of BsAPA. However, heat inactivation still limits the application of BsAPA. To further improve its thermostability, an autocatalysis site (L205-F206) in BsAPA was identified and three mutants (F193W, K204P, and A371V) were generated based on the analysis of the structure neighboring the autocatalysis site. These mutants have improved thermostability, and their half-life at 75 °C increased by 0.5-, 0.2-, and 0.3-fold, respectively. A triple-site mutant (F193W/K204P/A371V) was generated, with 1.5-fold increased half-life at 80 and a 10.7 °C increased T_m, compared with those of the wild-type. These results indicate that autocatalysis of aspartic protease reduces enzyme thermostability. Furthermore, site-directed mutagenesis at regions near the autocatalysis site is an efficient approach to improve aspartic protease thermostability.

Single point mutations reveal amino acid residues important for Chromobacterium violaceum transaminase activity in the production of unnatural amino acids

Unnatural amino acids (UAAs) are chiral amines with high application potential in drug discovery and synthesis of other valuable chemicals. Biocatalysis offers the possibility to synthesise novel optically pure UAAs with different physical and chemical properties. While the biocatalytic potential of transaminases in the synthesis of UAAs has been demonstrated, there is still a need to improve the activity with non-native substrates and to understand which amino acids residues are important for activity with these UAAs. Using a rational design approach, six variants of Chromobacterium violaceum DSM30191 transaminase (CV_TA) carrying a single and one variant carrying two substitutions were generated. Among the variants with a single substitution, CV_Y168F showed a 2 to 2.6-fold increased affinity for 2-oxooctanoic acid (2-OOA) and 3-oxobutyric acid (3-OBA) methyl ester used to synthesise an α- and β-UAA. Analysis of the first half of the transaminase reaction showed no change in the activity with the donor (S)-1-phenylethylamine. The combination of W60C and Y168F substitutions improved the CV_TA affinity for 2-OOA 10-fold compared to the wild type. Other substitutions showed no change, or reduced activity with the tested substrates. Our findings provide structural information on CV_TA and demonstrate the potential of rational design for biosynthesis of UAAs.

Coupled molecular dynamics mediate long- and short-range epistasis between mutations that affect stability and aggregation kinetics

Multiple mutations are typically required to significantly improve protein stability or aggregation kinetics. However, when several substitutions are made in a single protein, the mutations can potentially interact in a nonadditive manner, resulting in epistatic effects, which can hamper protein-engineering strategies to improve thermostability or aggregation kinetics. Here, we have examined the role of protein dynamics in mediating epistasis between pairs of mutations. With Escherichia coli transketolase (TK) as a model, we explored the epistatic interactions between two single variants H192P and A282P, and also between the double-mutant H192P/A282P and two single variants, I365L or G506A. Epistasis was determined for several measures of protein stability, including the following: the free-energy barrier to kinetic inactivation, ∆∆G ^‡; thermal transition midpoint temperatures, T _m; and aggregation onset temperatures, T _agg Nonadditive epistasis was observed between neighboring mutations as expected, but also for distant mutations located in the surface and core regions of different domains. Surprisingly, the epistatic behaviors for each measure of stability were often different for any given pairwise recombination, highlighting that kinetic and thermodynamic stabilities do not always depend on the same structural features. Molecular-dynamics simulations and a pairwise cross-correlation analysis revealed that mutations influence the dynamics of their local environment, but also in some cases the dynamics of regions distant in the structure. This effect was found to mediate epistatic interactions between distant mutations and could therefore be exploited in future protein-engineering strategies.

In vitro evolution of an l-amino acid deaminase active on l-1-naphthylalanine

l-Amino acid deaminase from Proteus myxofaciens (PmaLAAD) is a promising biocatalyst for enantioselective biocatalysis that can be exploited to produce optically pure d-amino acids or α-keto acids.

Computational Identification of Amino-Acid Mutations that Further Improve the Activity of a Chalcone-Flavonone Isomerase from Glycine max

Protein design for improving enzymatic activity remains a challenge in biochemistry, especially to identify target amino-acid sites for mutagenesis and to design beneficial mutations for those sites. Here, we employ a computational approach that combines multiple sequence alignment, positive selection detection, and molecular docking to identify and design beneficial amino-acid mutations that further improve the intramolecular-cyclization activity of a chalcone-flavonone isomerase from Glycine max (GmCHI). By this approach, two GmCHI mutants with higher activities were predicted and verified. The results demonstrate that this approach could determine the beneficial amino-acid mutations for improving the enzymatic activity, and may find more applications in engineering of enzymes.

Structural Analysis of an Evolved Transketolase Reveals Divergent Binding Modes

The S385Y/D469T/R520Q variant of E. coli transketolase was evolved previously with three successive smart libraries, each guided by different structural, bioinformatical or computational methods. Substrate-walking progressively shifted the target acceptor substrate from phosphorylated aldehydes, towards a non-phosphorylated polar aldehyde, a non-polar aliphatic aldehyde, and finally a non-polar aromatic aldehyde. Kinetic evaluations on three benzaldehyde derivatives, suggested that their active-site binding was differentially sensitive to the S385Y mutation. Docking into mutants generated in silico from the wild-type crystal structure was not wholly satisfactory, as errors accumulated with successive mutations, and hampered further smart-library designs. Here we report the crystal structure of the S385Y/D469T/R520Q variant, and molecular docking of three substrates. This now supports our original hypothesis that directed-evolution had generated an evolutionary intermediate with divergent binding modes for the three aromatic aldehydes tested. The new active site contained two binding pockets supporting π-π stacking interactions, sterically separated by the D469T mutation. While 3-formylbenzoic acid (3-FBA) preferred one pocket, and 4-FBA the other, the less well-accepted substrate 3-hydroxybenzaldehyde (3-HBA) was caught in limbo with equal preference for the two pockets. This work highlights the value of obtaining crystal structures of evolved enzyme variants, for continued and reliable use of smart library strategies.

Single active-site mutants are sufficient to enhance serine:pyruvate α-transaminase activity in an ω-transaminase

We have analyzed the natural evolution of transaminase structure and sequence between an α-transaminase serine-pyruvate aminotransferase and an ω-transaminase from Chromobacterium violaceum with < 20% sequence identity, and identified the active-site regions that are least conserved structurally. We also show that these structural changes correlate strongly with transaminase substrate specificity during evolution and therefore might normally be presumed to be essential determinants of substrate specificity. However, key residues are often conserved spatially during evolution and yet originate from within a different region of the sequence via structural reorganizations. In the present study, we also show that α-transaminase-type serine-pyruvate aminotransferase activity can be engineered into the CV2025 ω-transaminase scaffold with any one of many possible single-point mutations at three key positions, without the requirement for significant backbone remodeling, or repositioning of the residue from a different region of sequence. This finding has significant implications for enzyme redesign in which solutions to substrate specificity changes may be found more efficiently than is achieved by engineering in all sequence and structure determinants identified by correlation to substrate specificity.

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.

Use of nonionic surfactants for improvement of terpene production in Saccharomyces cerevisiae

To facilitate enzyme and pathway engineering, a selection was developed for improved sesquiterpene titers in Saccharomyces cerevisiae. α-Bisabolene, a candidate advanced biofuel, was found to protect yeast against the disruptive action of nonionic surfactants such as Tween 20 (T20). An experiment employing competition between two strains of yeast, one of which makes twice as much bisabolene as the other, demonstrated that growth in the presence of T20 provided sufficient selective pressure to enrich the high-titer strain to form 97% of the population. Following this, various methods were used to mutagenize the bisabolene synthase (BIS) coding sequence, coupled with selection by subculturing in the presence of T20. Mutagenesis targeting the BIS active site did not yield an improvement in bisabolene titers, although mutants were found which made a mixture of α-bisabolene and β-farnesene, another candidate biofuel. Based on evidence that the 3' end of the BIS mRNA may be unstable in yeast, we randomly recoded the last 20 amino acids of the enzyme and, following selection in T20, found a variant which increased specific production of bisabolene by more than 30%. Since T20 could enrich a mixed population, efficiently removing strains that produced little or no bisabolene, we investigated whether it could also be applied to sustain high product titers in a monoculture for an extended period. Cultures grown in the presence of T20 for 14 days produced bisabolene at titers up to 4-fold higher than cultures grown with an overlay of dodecane, used to sequester the terpene product, and 20-fold higher than cultures grown without dodecane.

Computational tools for designing smart libraries

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

Changing the specificity of α-amino acid ester hydrolase toward para-hydroxyl cephalosporins synthesis by site-directed saturation mutagenesis

α-Amino acid ester hydrolases (AEHs) catalyze the synthesis of β-lactam antibiotics containing an α-amino group with decreased activity toward antibiotics with a p-hydroxyl group. The AEH gene from Xanthomonas rubrillineans was cloned and expressed in Escherichia coli. Based on the crystal structure of the AEH and cefprozil complex, 13 residues not directly involved in substrate recognition were mutated individually. The resulting ~1,300 mutants were screened for activity using cefprozil as a model product based on spectrophotometric assay in a 96-well format. Mutants with improved cefprozil synthetic activity revealed the particular importance of positions 87, 131 and 175 for specificity. The mutant V131S with the highest initial rates of synthesis toward three p-hydroxyl cephalosporins showed 23 %, 17 % and 64 % increase in maximum product accumulation of cefadroxil, cefprozil and cefatrizine, respectively.

Directed evolution to re-adapt a co-evolved network within an enzyme

We have previously used targeted active-site saturation mutagenesis to identify a number of transketolase single mutants that improved activity towards either glycolaldehyde (GA), or the non-natural substrate propionaldehyde (PA). Here, all attempts to recombine the singles into double mutants led to unexpected losses of specific activity towards both substrates. A typical trade-off occurred between soluble expression levels and specific activity for all single mutants, but many double mutants decreased both properties more severely suggesting a critical loss of protein stability or native folding. Statistical coupling analysis (SCA) of a large multiple sequence alignment revealed a network of nine co-evolved residues that affected all but one double mutant. Such networks maintain important functional properties such as activity, specificity, folding, stability, and solubility and may be rapidly disrupted by introducing one or more non-naturally occurring mutations. To identify variants of this network that would accept and improve upon our best D469 mutants for activity towards PA, we created a library of random single, double and triple mutants across seven of the co-evolved residues, combining our D469 variants with only naturally occurring mutations at the remaining sites. A triple mutant cluster at D469, E498 and R520 was found to behave synergistically for the specific activity towards PA. Protein expression was severely reduced by E498D and improved by R520Q, yet variants containing both mutations led to improved specific activity and enzyme expression, but with loss of solubility and the formation of inclusion bodies. D469S and R520Q combined synergistically to improve k_cat 20-fold for PA, more than for any previous transketolase mutant. R520Q also doubled the specific activity of the previously identified D469T to create our most active transketolase mutant to date. Our results show that recombining active-site mutants obtained by saturation mutagenesis can rapidly destabilise critical networks of co-evolved residues, whereas beneficial single mutants can be retained and improved upon by randomly recombining them with natural variants at other positions in the network.

Assessing directed evolution methods for the generation of biosynthetic enzymes with potential in drug biosynthesis

To address the synthesis of increasingly structurally diverse small-molecule drugs, methods for the generation of efficient and selective biological catalysts are becoming increasingly important. 'Directed evolution' is an umbrella term referring to a variety of methods for improving or altering the function of enzymes using a nature-inspired twofold strategy of mutagenesis followed by selection. This article provides an objective assessment of the effectiveness of directed evolution campaigns in generating enzymes with improved catalytic parameters for new substrates from the last decade, excluding studies that aimed to select for only improved physical properties and those that lack kinetic characterization. An analysis of the trends of methodologies and their success rates from 81 qualifying examples in the literature reveals the average fold improvement for k (cat) (or V (max)), K (m) and k (cat)/K (m) to be 366-, 12- and 2548-fold, respectively, whereas the median fold improvements are 5.4, 3 and 15.6. Further analysis by enzyme class, library-generation methodology and screening methodology explores relationships between successful campaigns and the methodologies employed.

Switch of substrate specificity of hyperthermophilic acylaminoacyl peptidase by combination of protein and solvent engineering

The inherent evolvability of promiscuous enzymes endows them with great potential to be artificially evolved for novel functions. Previously, we succeeded in transforming a promiscuous acylaminoacyl peptidase (apAAP) from the hyperthermophilic archaeon Aeropyrum pernix K1 into a specific carboxylesterase by making a single mutation. In order to fulfill the urgent requirement of thermostable lipolytic enzymes, in this paper we describe how the substrate preference of apAAP can be further changed from p-nitrophenyl caprylate (pNP-C8) to p-nitrophenyl laurate (pNP-C12) by protein and solvent engineering. After one round of directed evolution and subsequent saturation mutagenesis at selected residues in the active site, three variants with enhanced activity towards pNP-C12 were identified. Additionally, a combined mutant W474V/F488G/R526V/T560W was generated, which had the highest catalytic efficiency (k (cat)/K (m)) for pNP-C12, about 71-fold higher than the wild type. Its activity was further increased by solvent engineering, resulting in an activity enhancement of 280-fold compared with the wild type in the presence of 30% DMSO. The structural basis for the improved activity was studied by substrate docking and molecular dynamics simulation. It was revealed that W474V and F488G mutations caused a significant change in the geometry of the active center, which may facilitate binding and subsequent hydrolysis of bulky substrates. In conclusion, the combination of protein and solvent engineering may be an effective approach to improve the activities of promiscuous enzymes and could be used to create naturally rare hyperthermophilic enzymes.

Iterative saturation mutagenesis accelerates laboratory evolution of enzyme stereoselectivity: rigorous comparison with traditional methods

Efficacy in laboratory evolution of enzymes is currently a pressing issue, making comparative studies of different methods and strategies mandatory. Recent reports indicate that iterative saturation mutagenesis (ISM) provides a means to accelerate directed evolution of stereoselectivity and thermostability, but statistically meaningful comparisons with other methods have not been documented to date. In the present study, the efficacy of ISM has been rigorously tested by applying it to the previously most systematically studied enzyme in directed evolution, the lipase from Pseudomonas aeruginosa as a catalyst in the stereoselective hydrolytic kinetic resolution of a chiral ester. Upon screening only 10,000 transformants, unprecedented enantioselectivity was achieved (E = 594). ISM proves to be considerably more efficient than all previous systematic efforts utilizing error-prone polymerase chain reaction at different mutation rates, saturation mutagenesis at hot spots, and/or DNA shuffling, pronounced positive epistatic effects being the underlying reason.

Exploring protein fitness landscapes by directed evolution

Directed evolution circumvents our profound ignorance of how a protein's sequence encodes its function by using iterative rounds of random mutation and artificial selection to discover new and useful proteins. Proteins can be tuned to adapt to new functions or environments by simple adaptive walks involving small numbers of mutations. Directed evolution studies have shown how rapidly some proteins can evolve under strong selection pressures and, because the entire 'fossil record' of evolutionary intermediates is available for detailed study, they have provided new insight into the relationship between sequence and function. Directed evolution has also shown how mutations that are functionally neutral can set the stage for further adaptation.