Two-substrate enzyme engineering using small libraries that combine the substrate preferences from two different variant lineages

Improving the range of substrates accepted by enzymes with high catalytic activity remains an important goal for the industrialisation of biocatalysis. Many enzymes catalyse two-substrate reactions which increases the complexity in engineering them for the synthesis of alternative products. Often mutations are found independently that can improve the acceptance of alternatives to each of the two substrates. Ideally, we would be able to combine mutations identified for each of the two alternative substrates, and so reprogramme new enzyme variants that synthesise specific products from their respective two-substrate combinations. However, as we have previously observed for E. coli transketolase, the mutations that improved activity towards aromatic acceptor aldehydes, did not successfully recombine with mutations that switched the donor substrate to pyruvate. This likely results from several active site residues having multiple roles that can affect both of the substrates, as well as structural interactions between the mutations themselves. Here, we have designed small libraries, including both natural and non-natural amino acids, based on the previous mutational sites that impact on acceptance of the two substrates, to achieve up to 630× increases in kcat for the reaction with 3-formylbenzoic acid (3-FBA) and pyruvate. Computational docking was able to determine how the mutations shaped the active site to improve the proximity of the 3-FBA substrate relative to the enamine-TPP intermediate, formed after the initial reaction with pyruvate. This work opens the way for small libraries to rapidly reprogramme enzyme active sites in a plug and play approach to catalyse new combinations of two-substrate reactions.

Selectivity in Enzymatic Phosphorus Recycling from Biopolymers: Isotope Effect, Reactivity Kinetics, and Molecular Docking with Fungal and Plant Phosphatases

Among ubiquitous phosphorus (P) reserves in environmental matrices are ribonucleic acid (RNA) and polyphosphate (polyP), which are, respectively, organic and inorganic P-containing biopolymers. Relevant to P recycling from these biopolymers, much remains unknown about the kinetics and mechanisms of different acid phosphatases (APs) secreted by plants and soil microorganisms. Here we investigated RNA and polyP dephosphorylation by two common APs, a plant purple AP (PAP) from sweet potato and a fungal phytase from Aspergillus niger. Trends of δ¹⁸O values in released orthophosphate during each enzyme-catalyzed reaction in ¹⁸O-water implied a different extent of reactivity. Subsequent enzyme kinetics experiments revealed that A. niger phytase had 10-fold higher maximum rate for polyP dephosphorylation than the sweet potato PAP, whereas the sweet potato PAP dephosphorylated RNA at a 6-fold faster rate than A. niger phytase. Both enzymes had up to 3 orders of magnitude lower reactivity for RNA than for polyP. We determined a combined phosphodiesterase-monoesterase mechanism for RNA and terminal phosphatase mechanism for polyP using high-resolution mass spectrometry and ³¹P nuclear magnetic resonance, respectively. Molecular modeling with eight plant and fungal AP structures predicted substrate binding interactions consistent with the relative reactivity kinetics. Our findings implied a hierarchy in enzymatic P recycling from P-polymers by phosphatases from different biological origins, thereby influencing the relatively longer residence time of RNA versus polyP in environmental matrices. This research further sheds light on engineering strategies to enhance enzymatic recycling of biopolymer-derived P, in addition to advancing environmental predictions of this P recycling by plants and microorganisms.

Characterisation of a hyperthermophilic transketolase from Thermotoga maritima DSM3109 as a biocatalyst for 7-keto-octuronic acid synthesis

Transketolase (TK) is a fundamentally important enzyme in industrial biocatalysis which carries out a stereospecific carbon-carbon bond formation, and is widely used in the synthesis of prochiral ketones. This study describes the biochemical and molecular characterisation of a novel and unusual hyperthermophilic TK from Thermotoga maritima DSM3109 (TKtmar). TKtmar has a low protein sequence homology compared to the already described TKs, with key amino acid residues in the active site highly conserved. TKtmar has a very high optimum temperature (>90 °C) and shows pronounced stability at high temperature (e.g. t1/2 99 and 9.3 h at 50 and 80 °C, respectively) and in presence of organic solvents commonly used in industry (DMSO, acetonitrile and methanol). Substrate screening showed activity towards several monosaccharides and aliphatic aldehydes. In addition, for the first time, TK specificity towards uronic acids was achieved with TKtmar catalysing the efficient conversion of d-galacturonic acid and lithium hydroxypyruvate into 7-keto-octuronic acid, a very rare C8 uronic acid, in high yields (98%, 49 mM).

A 'Split-Gene' Transketolase From the Hyper-Thermophilic Bacterium Carboxydothermus hydrogenoformans: Structure and Biochemical Characterization

A novel transketolase has been reconstituted from two separate polypeptide chains encoded by a ‘split-gene’ identified in the genome of the hyperthermophilic bacterium, Carboxydothermus hydrogenoformans. The reconstituted active α₂β₂ tetrameric enzyme has been biochemically characterized and its activity has been determined using a range of aldehydes including glycolaldehyde, phenylacetaldehyde and cyclohexanecarboxaldehyde as the ketol acceptor and hydroxypyruvate as the donor. This reaction proceeds to near 100% completion due to the release of the product carbon dioxide and can be used for the synthesis of a range of sugars of interest to the pharmaceutical industry. This novel reconstituted transketolase is thermally stable with no loss of activity after incubation for 1 h at 70°C and is stable after 1 h incubation with 50% of the organic solvents methanol, ethanol, isopropanol, DMSO, acetonitrile and acetone. The X-ray structure of the holo reconstituted α₂β₂ tetrameric transketolase has been determined to 1.4 Å resolution. In addition, the structure of an inactive tetrameric β₄ protein has been determined to 1.9 Å resolution. The structure of the active reconstituted α₂β₂ enzyme has been compared to the structures of related enzymes; the E1 component of the pyruvate dehydrogenase complex and D-xylulose-5-phosphate synthase, in an attempt to rationalize differences in structure and substrate specificity between these enzymes. This is the first example of a reconstituted ‘split-gene’ transketolase to be biochemically and structurally characterized allowing its potential for industrial biocatalysis to be evaluated.

Fine-tuning the activity and stability of an evolved enzyme active-site through noncanonical amino-acids

Site-specific saturation mutagenesis within enzyme active sites can radically alter reaction specificity, though often with a trade-off in stability. Extending saturation mutagenesis with a range of noncanonical amino acids (ncAA) potentially increases the ability to improve activity and stability simultaneously. Previously, an Escherichia coli transketolase variant (S385Y/D469T/R520Q) was evolved to accept aromatic aldehydes not converted by wild-type. The aromatic residue Y385 was critical to the new acceptor substrate binding, and so was explored here beyond the natural aromatic residues, to probe side chain structure and electronics effects on enzyme function and stability. A series of five variants introduced decreasing aromatic ring electron density at position 385 in the order para-aminophenylalanine (pAMF), tyrosine (Y), phenylalanine (F), para-cyanophenylalanine (pCNF) and para-nitrophenylalanine (pNTF), and simultaneously modified the hydrogen-bonding potential of the aromatic substituent from accepting to donating. The fine-tuning of residue 385 yielded variants with a 43-fold increase in specific activity for 50 mm 3-HBA and 100% increased k_cat (pCNF), 290% improvement in K_m (pNTF), 240% improvement in k_cat /K_m (pAMF) and decreased substrate inhibition relative to Y. Structural modelling suggested switching of the ring-substituted functional group, from donating to accepting, stabilised a helix-turn (D259-H261) through an intersubunit H-bond with G262, to give a 7.8 °C increase in the thermal transition mid-point, T_m , and improved packing of pAMF. This is one of the first examples in which both catalytic activity and stability are simultaneously improved via site-specific ncAA incorporation into an enzyme active site, and further demonstrates the benefits of expanding designer libraries to include ncAAs.

Oxalyl-CoA Decarboxylase Enables Nucleophilic One-Carbon Extension of Aldehydes to Chiral α-Hydroxy Acids

The synthesis of complex molecules from simple, renewable carbon units is the goal of a sustainable economy. Here we explored the biocatalytic potential of the thiamine-diphosphate-dependent (ThDP) oxalyl-CoA decarboxylase (OXC)/2-hydroxyacyl-CoA lyase (HACL) superfamily that naturally catalyzes the shortening of acyl-CoA thioester substrates through the release of the C1 -unit formyl-CoA. We show that the OXC/HACL superfamily contains promiscuous members that can be reversed to perform nucleophilic C1 -extensions of various aldehydes to yield the corresponding 2-hydroxyacyl-CoA thioesters. We improved the catalytic properties of Methylorubrum extorquens OXC by rational enzyme engineering and combined it with two newly described enzymes-a specific oxalyl-CoA synthetase and a 2-hydroxyacyl-CoA thioesterase. This enzymatic cascade enabled continuous conversion of oxalate and aromatic aldehydes into valuable (S)-α-hydroxy acids with enantiomeric excess up to 99 %.

Exploiting correlated molecular-dynamics networks to counteract enzyme activity-stability trade-off

Rigidifying flexible sites is a powerful method to improve enzyme stability. However, if the highly flexible regions form the active site, modifying them risks losing activity due to the activity–stability trade-off. We hypothesized here that regions outside the active site whose dynamics were highly correlated to flexible active sites, would provide good targets for stabilizing mutations. To test this hypothesis, six variants were constructed in the 3M variant of Escherichia coli transketolase. The best variant had a 10.8-fold improved half-life at 55 °C, and increased the T_m and T_agg by 3 °C and 4.3 °C, respectively. The variants even increased the activity, by up to threefold. This study highlights how protein engineering strategies could be potentially improved by considering long-range dynamics.

The directed evolution of enzymes for improved activity or substrate specificity commonly leads to a trade-off in stability. We have identified an activity–stability trade-off and a loss in unfolding cooperativity for a variant (3M) of Escherichia coli transketolase (TK) engineered to accept aromatic substrates. Molecular dynamics simulations of 3M revealed increased flexibility in several interconnected active-site regions that also form part of the dimer interface. Mutating the newly flexible active-site residues to regain stability risked losing the new activity. We hypothesized that stabilizing mutations could be targeted to residues outside of the active site, whose dynamics were correlated with the newly flexible active-site residues. We previously stabilized WT TK by targeting mutations to highly flexible regions. These regions were much less flexible in 3M and would not have been selected a priori as targets using the same strategy based on flexibility alone. However, their dynamics were highly correlated with the newly flexible active-site regions of 3M. Introducing the previous mutations into 3M reestablished the WT level of stability and unfolding cooperativity, giving a 10.8-fold improved half-life at 55 °C, and increased midpoint and aggregation onset temperatures by 3 °C and 4.3 °C, respectively. Even the activity toward aromatic aldehydes increased up to threefold. Molecular dynamics simulations confirmed that the mutations rigidified the active-site via the correlated network. This work provides insights into the impact of rigidifying mutations within highly correlated dynamic networks that could also be useful for developing improved computational protein engineering strategies.

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.

Biotransformation of β-hydroxypyruvate and glycolaldehyde to l-erythrulose by Pichia pastoris strain GS115 overexpressing native transketolase

Transketolase is a proven biocatalytic tool for asymmetric carbon-carbon bond formation, both as a purified enzyme and within bacterial whole-cell biocatalysts. The performance of Pichia pastoris as a host for transketolase whole-cell biocatalysis was investigated using a transketolase-overexpressing strain to catalyze formation of l-erythrulose from β-hydroxypyruvic acid and glycolaldehyde substrates. Pichia pastoris transketolase coding sequence from the locus PAS_chr1-4_0150 was subcloned downstream of the methanol-inducible AOX1 promoter in a plasmid for transformation of strain GS115, generating strain TK150. Whole and disrupted TK150 cells from shake flasks achieved 62% and 65% conversion, respectively, under optimal pH and methanol induction conditions. In a 300 μL reaction, TK150 samples from a 1L fed-batch fermentation achieved a maximum l-erythrulose space time yield (STY) of 46.58 g L^-1 h^-1 , specific activity of 155 U gCDW-1, product yield on substrate (Y_p/s ) of 0.52 mol mol^-1 and product yield on catalyst (Y_p/x ) of 2.23g gCDW-1. We have successfully exploited the rapid growth and high biomass characteristics of Pichia pastoris in whole cell biocatalysis. At high cell density, the engineered TK150 Pichia pastoris strain tolerated high concentrations of substrate and product to achieve high STY of the chiral sugar l-erythrulose. © 2017 The Authors Biotechnology Progress published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers Biotechnol. Prog., 34:99-106, 2018.

Two strategies to engineer flexible loops for improved enzyme thermostability

Flexible sites are potential targets for engineering the stability of enzymes. Nevertheless, the success rate of the rigidifying flexible sites (RFS) strategy is still low due to a limited understanding of how to determine the best mutation candidates. In this study, two parallel strategies were applied to identify mutation candidates within the flexible loops of Escherichia coli transketolase (TK). The first was a “back to consensus mutations” approach, and the second was computational design based on ΔΔG calculations in Rosetta. Forty-nine single variants were generated and characterised experimentally. From these, three single-variants I189H, A282P, D143K were found to be more thermostable than wild-type TK. The combination of A282P with H192P, a variant constructed previously, resulted in the best all-round variant with a 3-fold improved half-life at 60 °C, 5-fold increased specific activity at 65 °C, 1.3-fold improved k_cat and a T_m increased by 5 °C above that of wild type. Based on a statistical analysis of the stability changes for all variants, the qualitative prediction accuracy of the Rosetta program reached 65.3%. Both of the two strategies investigated were useful in guiding mutation candidates to flexible loops, and had the potential to be used for other enzymes.

Droplet millifluidics for kinetic study of transketolase

We present a continuous-flow reactor at the millifluidic scale coupled with an online, non-intrusive spectroscopic monitoring method for determining the kinetic parameters of an enzyme, transketolase (TK) used in biocatalysis for the synthesis of polyols by carboligation. The millifluidic system used is based on droplet flow, a well-established method for kinetic chemical data acquisition. The TK assay is based on the direct quantitative measurement of bicarbonate ions released during the transketolase-catalysed reaction in the presence of hydroxypyruvic acid as the donor, thanks to an irreversible reaction: bicarbonate ions react with phosphoenolpyruvate (PEP) in the presence of PEP carboxylase as the first auxiliary enzyme. The oxaloacetate formed is reduced to malate by NADH in the reaction catalysed by malate dehydrogenase as the second auxiliary enzyme. The extent of oxidation of NADH was measured by spectrophotometry at 340 nm. This system gives a direct, quantitative, generic method to evaluate the TK activity versus different substrates. We demonstrate the accuracy of this strategy to determine the enzymatic kinetic parameters and to study the substrate specificity of a thermostable TK from thermophilic microorganism Geobacillus stearothermophilus, offering promising prospects in biocatalysis. Millifluidic systems are useful in this regard as they can be used to rapidly evaluate the TK activity towards various substrates, and also different sets of conditions, identifying the optimal operating environment while minimizing resource consumption and ensuring high control over the operating conditions.

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.

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.

Microfluidic multi-input reactor for biocatalytic synthesis using transketolase

Biocatalytic synthesis in continuous-flow microreactors is of increasing interest for the production of specialty chemicals. However, the yield of production achievable in these reactors can be limited by the adverse effects of high substrate concentration on the biocatalyst, including inhibition and denaturation. Fed-batch reactors have been developed in order to overcome this problem, but no continuous-flow solution exists. We present the design of a novel multi-input microfluidic reactor, capable of substrate feeding at multiple points, as a first step towards overcoming these problems in a continuous-flow setting. Using the transketolase-(TK) catalysed reaction of lithium hydroxypyruvate (HPA) and glycolaldehyde (GA) to l-erythrulose (ERY), we demonstrate the transposition of a fed-batch substrate feeding strategy to our microfluidic reactor. We obtained a 4.5-fold increase in output concentration and a 5-fold increase in throughput compared with a single input reactor.

CC bond formation using ThDP-dependent lyases

The present review summarizes recent achievements in enzymatic thiamine catalysis during the past three years. With well-established enzymes such as BAL, PDC and TK new reactions have been identified and respective variants were prepared, which enable access to stereoisomeric products. Further we highlight recent progress with 'new' ThDP-dependent enzymes like MenD and PigD, which catalyze the Stetter-like 1,4 addition of aldehydes and YerE, which is the first known ThDP-dependent enzyme accepting ketones as acceptors.