Robust ω-Transaminases by Computational Stabilization of the Subunit Interface

Thermostability and activity improvement in l-threonine aldolase through targeted mutations in V-shaped subunit

l-threonine aldolase (LTA) catalyzes the synthesis of β-hydroxy-α-amino acids, which are important chiral intermediates widely used in the fields of pharmaceuticals and pesticides. However, the limited thermostability of LTA hinders its industrial application. Furthermore, the trade-off between thermostability and activity presents a challenge in the thermostability engineering of this enzyme. This study proposes a strategy to regulate the rigidity of LTA's V-shaped subunit by modifying its opening and hinge regions, distant from the active center, aiming to mitigate the trade-off. With LTA from Bacillus nealsonii as targeted enzyme, a total of 25 residues in these two regions were investigated by directed evolution. Finally, mutant G85A/M207L/A12C was obtained, showing significantly enhanced thermostability with a 20 °C increase in T5060 to 66 °C, and specific activity elevated by 34 % at the optimum temperature. Molecular dynamics simulations showed that the newly formed hydrophobicity and hydrogen bonds improved the thermostability and boosted proton transfer efficiency. This work enhances the thermostability of LTA while preventing the loss of activity. It opens new avenues for the thermostability engineering of other industrially relevant enzymes with active center located at the interface of subunits or domains.

Engineering of a hydroxysteroid dehydrogenase with simultaneous enhancement in activity and thermostability for efficient biosynthesis of ursodeoxycholic acid

Hydroxysteroid dehydrogenases (HSDHs) catalyze the oxidation/reduction of hydroxyl/keto groups of steroids with high regio- or stereoselectivity, playing an essential role in producing optically pure chemicals. In this work, a novel approach was developed to simultaneously improve the stability and activity of 7β-hydroxysteroid dehydrogenase (7β-HSDH) by combining B-factor analysis and computer-aided prediction. Several advantageous mutants were identified, and the most promising variant, S51Y/P202Y, exhibited 2.3-fold improvements in catalytic activity, 3.3-fold in half-life at 40°C, and 4.7-fold in catalytic efficiency (kcat/Km), respectively. Structural modeling analysis showed that the shortened reversible oxidation reaction catalytic distance and the strengthened residue interactions compared to the wild type were attributed to the improved stability and activity of the obtained mutants. To synthesize ursodeoxycholic acid cost-effectively by mutant S51Y/P202Y, a NAD-kinase was employed to facilitate the substitution of nicotinamide adenine dinucleotide phosphate (NADP+) with nicotinamide adenine dinucleotide (NAD+) in the whole-cell catalysis system. The substrate 7-ketolithocholic acid (100 mM) was converted completely in 0.5 h, achieving a space-time yield of 1,887.3 g L-1 d-1. This work provided a general target-oriented strategy for obtaining stable and highly active dehydrogenase for efficient biosynthesis.

IMPORTANCE

Hydroxysteroid dehydrogenases have emerged as indispensable tools in the synthesis of steroids, bile acids, and other steroid derivatives for the pharmaceutical and chemical industries. In this study, a novel approach was developed to simultaneously improve the stability and activity of a hydroxysteroid dehydrogenase by combining B-factor analysis and computer-aided prediction. This semi-rational method was demonstrated to be highly effective for enzyme engineering. In addition, NAD kinase was introduced to convert NAD+ to NADP+ for effective coenzyme regeneration in the whole-cell multienzyme-catalyzed system. This strategy reduces the significant economic costs associated with externally supplemented cofactors in NADP-dependent biosynthetic pathways.

Efficient and easible biocatalysts: Strategies for enzyme improvement. A review

Owing to the environmental friendliness and vast advantages that enzymes offer in the biotechnology and industry fields, biocatalysts are a prolific investigation field. However, the low catalytic activity, stability, and specific selectivity of the enzyme limit the range of the reaction enzymes involved in. A comprehensive understanding of the protein structure and dynamics in terms of molecular details enables us to tackle these limitations effectively and enhance the catalytic activity by enzyme engineering or modifying the supports and solvents. Along with different strategies including computational, enzyme engineering based on DNA recombination, enzyme immobilization, additives, chemical modification, and physicochemical modification approaches can be promising for the wide spread of industrial enzyme usage. This is attributed to the successful application of biocatalysts in industrial and synthetic processes requires a system that exhibits stability, activity, and reusability in a continuous flow process, thereby reducing the production cost. The main goal of this review is to display relevant approaches for improving enzyme characteristics to overcome their industrial application.

Rational redesign of the loop dynamics of carbonyl reductase LfSDR1 to improve the stereoselectivity for asymmetric synthesis of bulky chiral alcohols

Engineering biocatalysts with enhanced stereoselectivity is highly desirable, and active-site loop dynamics play an important role in its regulation. However, knowledge of their precise roles in catalysis and evolution is limited. Here, we used the strategy of Rosetta enzyme design combined molecular dynamic simulations (MDs) to reprogram the landscapes of the key active-site loop dynamics of the carbonyl reductase LfSDR1 to improve stereoselectivity. The key flexible loop in the active site showed the potential to regulate the catalytic properties. A library of virtual variants was produced using the Rosetta design and assessed dynamic effect of the loop with the aid of MDs. A potential candidate was obtained with significant stereoselectivity (ee > 99 %) compared to the wild-type (ee = 42 %) without loss of catalytic activity or thermostability. The molecular basis of the catalytic property enhancement was flanked by MDs, which revealed the role of the G92L mutation in regulating loop dynamics to stabilize the environment of the active site. Finally, a series of the challenge bulky substrate derivatives were assessed using the G92L variant, and all showed improved stereoselectivity ee > 99 %. This study provides novel insights for improving stereoselectivity through rational engineering of the loop dynamics of biocatalysts.

Unlocking the potential of enzyme engineering via rational computational design strategies

Enzymes play a pivotal role in various industries by enabling efficient, eco-friendly, and sustainable chemical processes. However, the low turnover rates and poor substrate selectivity of enzymes limit their large-scale applications. Rational computational enzyme design, facilitated by computational algorithms, offers a more targeted and less labor-intensive approach. There has been notable advancement in employing rational computational protein engineering strategies to overcome these issues, it has not been comprehensively reviewed so far. This article reviews recent developments in rational computational enzyme design, categorizing them into three types: structure-based, sequence-based, and data-driven machine learning computational design. Case studies are presented to demonstrate successful enhancements in catalytic activity, stability, and substrate selectivity. Lastly, the article provides a thorough analysis of these approaches, highlights existing challenges and potential solutions, and offers insights into future development directions.

Lipases for targeted industrial applications, focusing on the development of biotechnologically significant aspects: A comprehensive review of recent trends in protein engineering

Lipases are remarkable biocatalysts, adept at catalyzing the breakdown of diverse compounds into glycerol, fatty acids, and mono- and di-glycerides via hydrolysis. Beyond this, they facilitate esterification, transesterification, alcoholysis, acidolysis, and more, making them versatile in industrial applications. In industrial processes, lipases that exhibit high stability are favored as they can withstand harsh conditions. However, most native lipases are unable to endure adverse conditions, making them unsuitable for industrial use. Protein engineering proves to be a potent technology in the development of lipases that can function effectively under challenging conditions and fulfill criteria for various industrial processes. This review concentrated on new trends in protein engineering to enhance the diversity of lipase genes and employed in silico methods for predicting and comprehensively analyzing target mutations in lipases. Additionally, key molecular factors associated with industrial characteristics of lipases, including thermostability, solvent tolerance, catalytic activity, and substrate preference have been elucidated. The present review delved into how industrial traits can be enhanced through directed evolution (epPCR, gene shuffling), rational design (FRESCO, ASR), combined engineering strategies (i.e. CAST, ISM, and FRISM) as protein engineering methodologies in contexts of biodiesel production, food processing, and applications of detergent, pharmaceutics, and plastic degradation.

Solvent concentration at 50% protein unfolding may reform enzyme stability ranking and process window identification

As water miscible organic co-solvents are often required for enzyme reactions to improve e.g., the solubility of the substrate in the aqueous medium, an enzyme is required which displays high stability in the presence of this co-solvent. Consequently, it is of utmost importance to identify the most suitable enzyme or the appropriate reaction conditions. Until now, the melting temperature is used in general as a measure for stability of enzymes. The experiments here show, that the melting temperature does not correlate to the activity observed in the presence of the solvent. As an alternative parameter, the concentration of the co-solvent at the point of 50% protein unfolding at a specific temperature T in shortc U 50 T is introduced. Analyzing a set of ene reductases,c U 50 T is shown to indicate the concentration of the co-solvent where also the activity of the enzyme drops fastest. Comparing possible rankings of enzymes according to melting temperature andc U 50 T reveals a clearly diverging outcome also depending on the specific solvent used. Additionally, plots ofc U 50 versus temperature enable a fast identification of possible reaction windows to deduce tolerated solvent concentrations and temperature.

Enhancing the Thermal Stability and Enzyme Activity of Ketopantoate Hydroxymethyltransferase through Interface Modification Engineering

Ketopantoate hydroxymethyltransferase (KPHMT) plays a pivotal role in d-pantothenic acid biosynthesis. Most KPHMTs are homodecamers with low thermal stability, posing challenges for protein engineering and limiting output enhancement. Previously, a high-enzyme activity KPHMT mutant (K25A/E189S) from Corynebacterium glutamicum was screened as mother strain (M0). Building upon this strain, our study focused on interface engineering modifications, employing a multifaceted approach including integrating folding-free energy calculation, B-factor analysis, and conserved site analysis. Preliminary screening led to the selection of five mutants in the interface─E106S, E98T, E98N, S247I, and S247D─showing improved thermal stability, culminating in the double-site mutant M8 (M0-E98N/S247D). M8 exhibited a T1/2 value of 288.79 min at 50 °C, showing a 3.29-fold increase compared to M0. Meanwhile, the Tm value of M8 was elevated from 53.2 to 59.6 °C. Investigations of structural and molecular dynamics simulations revealed alterations in surface electrostatic charge distribution and the formation of increased hydrogen bonds between subunits, contributing to enhanced thermal stability. This investigation corroborates the efficacy of interface engineering modifications in bolstering KPHMT stability while showing its potential for positively impacting industrial d-pantothenic acid synthesis.

Rational design of transaminases based on comparative analysis of catalytically active and distance-free modes of the high-energy intermediate state

Bioproduction of chiral amines is limited by low transaminase (TA) activity on nonnatural substrates, leading to the need for protein engineering. To address the challenge of quickly and precisely identifying the engineering targets, a strategy was proposed based on analyzing the mode changes in the high-energy intermediate state (H-state) of the substrate-enzyme complex during catalysis. By substituting the residues with minimal structural changes in catalytically active mode (A-mode) and distance-free mode (F-mode) of the H-state complex with more conserved ones to stabilize it, a TA mutant M5(T295C/L387A/V436A) with 121.9-fold higher activity for synthesizing the (S)-Rivastigmine precursor (S)-1-(3-methoxyphenyl)ethylamine was created. The applicability of this strategy was also validated by engineering another TA for 1.52-fold higher activity and >99% selectivity toward (R)-3-amino-1-butanol biopreparation. The much higher stereoselectivity of the mutant compared with the wild type (28.3%) demonstrated that this strategy is not only advantageous in engineering enzyme activity but also applicable for modulating stereoselectivity.

Computational design of an imine reductase: mechanism-guided stereoselectivity reversion and interface stabilization

Imine reductases (IREDs) are important biocatalysts in the asymmetric synthesis of chiral amines. However, a detailed understanding of the stereocontrol mechanism of IRED remains incomplete, making the design of IRED for producing the desired amine enantiomers challenging. In this study, we investigated the stereoselective catalytic mechanism and designed an (R)-stereoselective IRED from Paenibacillus mucilaginosus (PmIR) using pharmaceutically relevant 2-aryl-substituted pyrrolines as substrates. A putative mechanism for controlling stereoselectivity was proposed based on the crucial role of electrostatic interactions in controlling iminium cation orientation and employed to achieve complete inversion of stereoselectivity in PmIR using computational design. The variant PmIR-Re (Q138M/P140M/Y187E/Q190A/D250M/R251N) exhibited opposite (S)-stereoselectivity, with >96% enantiomeric excess (ee) towards tested 2-aryl-substituted pyrrolines. Computational tools were employed to identify stabilizing mutations at the interface between the two subunits. The variant PmIR-6P (P140A/Q190S/R251N/Q217E/A257R/T277M) showed a nearly 5-fold increase in activity and a 12 °C increase in melting temperature. The PmIR-6P successfully produced (R)-2-(2,5-difluorophenyl)-pyrrolidine, a key chiral pharmaceutical intermediate, at a concentration of 400 mM with an ee exceeding 99%. This study provides insight into the stereocontrol elements of IREDs and demonstrates the potential of computational design for tailored stereoselectivity and thermal stability.

Correlation Analysis of Key Residue Sites between Computational-Aided Design Thermostability d-Amino Acid Oxidase and Ancestral Enzymes

The d-amino acid oxidase (DAAO) from Rhodotorula taiwanensis has proven to have great potential for applications due to its excellent catalytic kinetic parameters. However, its poor thermal stability has limited its performance in biocatalysis. Herein, starting from the variant SHVG of RtwDAAO, this study employed a comprehensive computational design approach for protein stability engineering, resulting in positive substitutions at specific sites (A43S, T45M, C234L, E195Y). The generated variant combination, SHVG/SMLY, exhibited a significant synergistic effect, leading to an extension of the half-life and Tmapp. The ancestral sequence reconstruction revealed the conservation of the variant sites. The association of the variant sites with the highly stable ancestral enzyme was further explored. After determining the contribution of the variant sites to thermal stability, it was applied to other homologous sequences and validated. Molecular dynamics simulations indicated that the increased hydrophobicity of the variant SHVG/SMLY was a key factor for the increased stability, with strengthened intersubunit interactions playing an important role. In addition, the physical properties of the amino acids themselves were identified as crucial factors for thermal stability generality in homologous enzymes, which is important for the rapid acquisition of a series of stable enzymes.

Enhancement of the substrate specificity of D-amino acid oxidase based on tunnel-pocket engineering

D-Amino acid oxidase (DAAO) selectively catalyzes the oxidative deamination of  D-amino acids, making it one of the most promising routes for synthesizing optically pure  L-amino acids, including  L-phosphinothricin ( L-PPT), a chiral herbicide with significant market potential. However, the native DAAOs that have been reported have low activity against unnatural acid substrate  D-PPT. Herein, we designed and screened a DAAO from Rhodotorula taiwanensis (RtwDAAO), and improved its catalytic potential toward  D-PPT through protein engineering. A semirational design approach was employed to create a mutation library based on the tunnel-pocket engineering. After three rounds of iterative saturation mutagenesis, the optimal variant M3rd -SHVG was obtained, exhibiting a >2000-fold increase in relative activity. The kinetic parameters showed that M3rd -SHVG improved the substrate binding affinity and turnover number. This is the optimal parameter reported so far. Further, molecular dynamics simulation revealed that the M3rd -SHVG reshapes the tunnel-pocket and corrects the direction of enzyme-substrate binding, allowing efficiently catalyze unnatural substrates. Our strategy demonstrates that the redesign of tunnel-pockets is effective in improving the activity and kinetic efficiency of DAAO, which provides a valuable reference for enzymatic catalysis. With the M3rd -SHVG as biocatalyst, 500 mM D, L-PPT was completely converted and the yield reached 98%. The results laid the foundation for further industrial production.

Enzymatic Oxy- and Amino-Functionalization in Biocatalytic Cascade Synthesis: Recent Advances and Future Perspectives

Biocatalytic cascades are a powerful tool for building complex molecules containing oxygen and nitrogen functionalities. Moreover, the combination of multiple enzymes in one pot offers the possibility to minimize downstream processing and waste production. In this review, we illustrate various recent efforts in the development of multi-step syntheses involving C-O and C-N bond-forming enzymes to produce high value-added compounds, such as pharmaceuticals and polymer precursors. Both in vitro and in vivo examples are discussed, revealing the respective advantages and drawbacks. The use of engineered enzymes to boost the cascades outcome is also addressed and current co-substrate and cofactor recycling strategies are presented, highlighting the importance of atom economy. Finally, tools to overcome current challenges for multi-enzymatic oxy- and amino-functionalization reactions are discussed, including flow systems with immobilized biocatalysts and cascades in confined nanomaterials.

Enhancing the organic solvent resistance of ω-amine transaminase for enantioselective synthesis of (R)-(+)-1(1-naphthyl)-ethylamine

BACKGROUND

Biocatalysis in high-concentration organic solvents has been applied to produce various industrial products with many advantages. However, using enzymes in organic solvents often suffers from inactivation or decreased catalytic activity and stability. An R-selective ω-amine transaminase from Aspergillus terreus (AtATA) exhibited activity toward 1-acetylnaphthalene. However, AtATA displayed unsatisfactory organic solvent resistance, which is required to enhance the solubility of the hydrophobic substrate 1-acetylnaphthalene. So, improving the tolerance of enzymes in organic solvents is essential.

MAIN METHODS AND RESULTS

The method of regional random mutation combined with combinatorial mutation was used to improve the resistance of AtATA in organic solvents. Enzyme surface areas are structural elements that undergo reversible conformational transitions, thus affecting the stability of the enzyme in organic solvents. Herein, three surface areas containing three loops were selected as potential mutation regions. And the "best" mutant T23I/T200K/P260S (M3) was acquired. In different concentrations of dimethyl sulfoxide (DMSO), the catalytic efficiency (kcat /Km ) toward 1-acetylnaphthalene and the stability (half-life t1/2 ) were higher than the wild-type (WT) of AtATA. The results of decreased Root Mean Square Fluctuation (RMSF) values via 20-ns molecular dynamics (MD) simulations under 15%, 25%, 35%, and 45% DMSO revealed that mutant M3 had lower flexibility, acquiring a more stable protein structure and contributing to its organic solvents stability than WT. Furthermore, M3 was applied to convert 1-acetylnaphthalene for synthesizing (R)-(+)-1(1-naphthyl)-ethylamine ((R)-NEA), which was an intermediate of Cinacalcet Hydrochloride for the treatment of secondary hyperthyroidism and hypercalcemia. Moreover, in a 20-mL scale-up experiment, 10 mM 1-acetylnaphthalene can be converted to (R)-NEA with 85.2% yield and a strict R-stereoselectivity (enantiomeric excess (e.e.) value >99.5%) within 10 h under 25% DMSO.

CONCLUSION

The beneficial mutation sites were identified to tailor AtATA's organic solvents stability via regional random mutation. The "best" mutant T23I/T200K/P260S (M3) holds great potential application for the synthesis of (R)-NEA.

Simultaneously optimizing multiple properties of β-glucosidase Bgl6 using combined (semi-)rational design strategies and investigation of the underlying mechanisms

The performance of β-glucosidase during cellulose saccharification is determined by thermostability, activity and glucose tolerance. However, conflicts between them make it challenging to simultaneously optimize three properties. In this work, such a case was reported using Bgl6-M3 as a starting point. Firstly, four thermostability-enhancing mutations were obtained using computer-aided engineering strategies (mutant M7). Secondly, substrate binding pocket of M7 was reshaped, generating two mutations that increased activity but decreased glucose tolerance (mutant M9). Then a key region lining active site cavity was redesigned, resulting in three mutations that boosted glucose tolerance and activity. Finally, mutant M12 with simultaneously improved thermostability (half-life of 20-fold), activity (k_cat/K_m of 5.6-fold) and glucose tolerance (ΔIC₅₀ of 200 mM) was obtained. Mechanisms for property improvement were elucidated by structural analysis and molecular dynamics simulations. Overall, the strategies used here and new insights into the underlying mechanisms may provide guidance for multi-property engineering of other enzymes.

Turning thermostability of Aspergillus terreus (R)-selective transaminase At-ATA by synthetic shuffling

As versatile and green biocatalysts for the asymmetric amination of ketones, the insufficient thermostability of transaminases always limits its broad application in the pharmaceutical and fine chemical industries. Here, synthetic shuffling technology was used to enhance stability of (R)-selective transaminase from Aspergillus terreus. The results showed that 30 out of 5000 mutants had improved thermostability by color-based screening method, among which mutants with residual enzyme activity higher than 50% at 45 °C for 10 min were selected for further analysis. Especially, the half-inactivation temperature (T₅₀¹⁰), half-life (t_1/2), and melting temperature (T_m) of the best mutant M14 (M280C-H210N-M150C-F115L) were 13.7 °C, 165.8 min, and 13.9 °C higher than that of the wild type (WT), respectively. M14 also exhibited a significant biocatalytic efficiency toward acetophenone and 1-acetylnaphthalene, the yield of which were 265.6% and 117.5% higher than WT, respectively. Based on molecular dynamics simulation, improved catalytic efficiency of M14 could be attributed to its increased hydrogen bonds interaction around the mutation sites. Additionally, the introduction of disulfide bond combined with above mutations has a synergistic effect on the improved protein thermostability.

Improvement of the stability and catalytic efficiency of heparan sulfate N-sulfotransferase for preparing N-sulfated heparosan

The chemo-enzymatic and enzymatic synthesis of heparan sulfate and heparin are considered as an attractive alternative to the extraction of heparin from animal tissues. Sulfation of the hydroxyl group at position 2 of the deacetylated glucosamine is a prerequisite for subsequent enzymatic modifications. In this study, multiple strategies, including truncation mutagenesis based on B-factor values, site-directed mutagenesis guided by multiple sequence alignment, and structural analysis were performed to improve the stability and activity of human N-sulfotransferase. Eventually, a combined variant Mut02 (MBP-hNST-NΔ599-602/S637P/S741P/E839P/L842P/K779N/R782V) was successfully constructed, whose half-life at 37°C and catalytic activity were increased by 105-fold and 1.35-fold, respectively. After efficient overexpression using the Escherichia coli expression system, the variant Mut02 was applied to N-sulfation of the chemically deacetylated heparosan. The N-sulfation content reached around 82.87% which was nearly 1.88-fold higher than that of the wild-type. The variant Mut02 with high stability and catalytic efficiency has great potential for heparin biomanufacturing.

Efficient Exploration of Sequence Space by Sequence-Guided Protein Engineering and Design

The rapid growth of sequence databases over the past two decades means that protein engineers faced with optimizing a protein for any given task will often have immediate access to a vast number of related protein sequences. These sequences encode information about the evolutionary history of the protein and the underlying sequence requirements to produce folded, stable, and functional protein variants. Methods that can take advantage of this information are an increasingly important part of the protein engineering tool kit. In this Perspective, we discuss the utility of sequence data in protein engineering and design, focusing on recent advances in three main areas: the use of ancestral sequence reconstruction as an engineering tool to generate thermostable and multifunctional proteins, the use of sequence data to guide engineering of multipoint mutants by structure-based computational protein design, and the use of unlabeled sequence data for unsupervised and semisupervised machine learning, allowing the generation of diverse and functional protein sequences in unexplored regions of sequence space. Altogether, these methods enable the rapid exploration of sequence space within regions enriched with functional proteins and therefore have great potential for accelerating the engineering of stable, functional, and diverse proteins for industrial and biomedical applications.

Engineering the Thermostability of the Mono- and Diacylglycerol Lipase SMG1 for the Synthesis of Diacylglycerols

Diacylglycerols (DAGs) display huge application prospectives in food industries. Therefore, new strategies to produce diacylglycerides are needed. Malassezia globose lipase (SMG1) could be used to synthesize DAGs. However, the poor thermostability of SMG1 seriously hampers its application. Herein, a rational design was used to generate a more thermostable SMG1. Compared with the wild type (WT), the M5D mutant (Q34P/A37P/M176V/G177A/M294R/ G28C-P206C), which contains five single-point mutations and one additional disulfide bond, displayed a 14.0 °C increase in the melting temperature (T_m), 5 °C in the optimal temperature, and 1154.3-fold in the half-life (t_1/2) at 55 °C. Meanwhile, the specific activity towards DAGs of the M5D variant was improved by 3.0-fold compared to the WT. Molecular dynamics (MD) simulations revealed that the M5D mutant showed an improved rigid structure. Additionally, the WT and the M5D variants were immobilized and used for the production of DAGs. Compared with the WT, the immobilized M5D-catalyzed esterification showed a 9.1% higher DAG content and a 22.9% increase in residual activity after nine consecutive cycles. This study will pave the way for the industrial application of SMG1.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

Improving Both the Thermostability and Catalytic Efficiency of Phospholipase D from Moritella sp. JT01 through Disulfide Bond Engineering Strategy

Mining of Phospholipase D (PLD) with high activity and stability has attracted strong interest for investigation. A novel PLD from marine Moritella sp. JT01 (MsPLD) was biochemically and structurally characterized in our previous study; however, the short half-life time (t_1/2) under its optimum reaction temperature seriously hampered its further applications. Herein, the disulfide bond engineering strategy was applied to improve its thermostability. Compared with wild-type MsPLD, mutant S148C-T206C/D225C-A328C with the addition of two disulfide bonds exhibited a 3.1-fold t_1/2 at 35 °C and a 5.7 °C increase in melting temperature (T_m). Unexpectedly, its specific activity and catalytic efficiency (k_cat/K_m) also increased by 22.7% and 36.5%, respectively. The enhanced activity might be attributed to an increase in the activation entropy by displacing more water molecules by the transition state. The results of molecular dynamics simulations (MD) revealed that the introduction of double disulfide bonds rigidified the global structure of the mutant, which might cause the enhanced thermostability. Finally, the synthesis capacity of the mutant to synthesize phosphatidic acid (PA) was evaluated. The conversion rate of PA reached about 80% after 6 h reaction with wild-type MsPLD but reached 78% after 2 h with mutant S148C-T206C/D225C-A328C, which significantly reduced the time needed for the reaction to reach equilibrium. The present results pave the way for further application of MsPLD in the food and pharmaceutical industries.

Structure-based interface engineering methodology in designing a thermostable amylose-forming transglucosylase

Many drugs and prebiotics derive their activities from sugar substituents. Due to the prevalence and complexity of these biologically active compounds, enzymatic glycodiversification that facilitates easier access to these compounds can make profound contributions to the pharmaceutical, food, and feed industries. Amylosucrases (ASases) are attractive tools for glycodiversification because of their broad acceptor substrate specificity, but the lack of structural information and their poor thermostability limit their industrial applications. Herein, we reported the crystal structure of ASase from Calidithermus timidus, which displays a homotetrameric quaternary organization not previously observed for other ASases. We employed a workflow composed of five common strategies, including interface engineering, folding energy calculations, consensus sequence, hydrophobic effects enhancement, and B-factor analysis, to enhance the thermostability of C. timidus ASase. As a result, we obtained a quadruple-point mutant M31 ASase with a half-life at 65 °C increased from 22.91 h to 52.93 h, which could facilitate biosynthesis of glucans with a degree of polymerization of more than 20 using sucrose as a substrate at 50 °C. In conclusion, this study provides a structural basis for understanding the multifunctional biocatalyst ASase and presents a powerful methodology to effectively and systematically enhance protein thermostability.

Modeling-Assisted Design of Thermostable Benzaldehyde Lyases from Rhodococcus erythropolis for Continuous Production of α-Hydroxy Ketones

Enantiopure α-hydroxy ketones are important building blocks of active pharmaceutical ingredients (APIs), which can be produced by thiamine-diphosphate-dependent lyases, such as benzaldehyde lyase. Here we report the discovery of a novel thermostable benzaldehyde lyase from Rhodococcus erythropolis R138 (ReBAL). While the overall sequence identity to the only experimentally confirmed benzaldehyde lyase from Pseudomonas fluorescens Biovar I (PfBAL) was only 65 %, comparison of a structural model of ReBAL with the crystal structure of PfBAL revealed only four divergent amino acids in the substrate binding cavity. Based on rational design, we generated two ReBAL variants, which were characterized along with the wild-type enzyme in terms of their substrate spectrum, thermostability and biocatalytic performance in the presence of different co-solvents. We found that the new enzyme variants have a significantly higher thermostability (up to 22 °C increase in T50 ) and a different co-solvent-dependent activity. Using the most stable variant immobilized in packed-bed reactors via the SpyCatcher/SpyTag system, (R)-benzoin was synthesized from benzaldehyde over a period of seven days with a stable space-time-yield of 9.3 mmol ⋅ L-1  ⋅ d-1 . Our work expands the important class of benzaldehyde lyases and therefore contributes to the development of continuous biocatalytic processes for the production of α-hydroxy ketones and APIs.

Insights of conformational dynamics on catalytic activity in the computational stability design of Bacillus subtilis LipA

In protein engineering, the contributions of individual mutations to designed combinatorial mutants are unpredictable. Screening designed mutations that affect enzyme catalytic activity enables evolutions towards efficient activities. Here, Bacillus subtilis LipA (BSLA) was selected as a model protein for thermostabilization designs, and the circular dichroism measurements showed six combinatorial designs with improved stability (from 5.81 °C to 13.61 °C). Based on molecular dynamic simulations, the conformational dynamics of the mutants revealed that mutations alter the populations of conformational states and the increased ensembles of inactive conformations might lead to a reduction in activity. We further demonstrated that the mutations responsible for the reduced enzyme catalytic activity involved a short dynamic correlation path to disturbing the equilibrium conformation of active sites. By removing N82V, which had a close dynamic correlation to the active sites in mutant D3, the redesigned mutant RD3 had an increased activity of 57.6%. By combining computational simulation with experimental verification, this work established that essential sites to counteract the activity-stability trade-off in multipoint combinatorial mutants could be computationally predicted and thus provide a possible strategy by which to indirectly or directly guide protein design.

Enhanced thermal and alkaline stability of L-lysine decarboxylase CadA by combining directed evolution and computation-guided virtual screening

As an important monomer for bio-based nylons PA5X, cadaverine is mainly produced by enzymatic decarboxylation of L-lysine. A key issue with this process is the instability of L-lysine decarboxylase (CadA) during the reaction due to the dissociation of CadA subunits with the accumulation of alkaline cadaverine. In this work, we attempted to improve the thermal and alkaline stability of CadA by combining directed evolution and computation-guided virtual screening. Interestingly, site 477 residue located at the protein surface and not the decamer interface was found as a hotspot in directed evolution. By combinatorial mutagenesis of the positive mutations obtained by directed evolution and virtual screening with the previously reported T88S mutation, K477R/E445Q/T88S/F102V was generated as the best mutant, delivering 37% improvement of cadaverine yield at 50 ºC and pH 8.0. Molecular dynamics simulations suggested the improved rigidity of regional structures, increased number of salt bridges, and enhancement of hydrogen bonds at the multimeric interface as possible origins of the improved stability of the mutant. Using this four-point mutant, 160.7 g/L of cadaverine was produced from 2.0 M Lysine hydrochloride at 50 °C without pH regulation, with a conversion of 78.5%, whereas the wild type produced 143.7 g/L cadaverine, corresponding to 70% conversion. This work shows the combination of directed evolution and virtual screening as an efficient protein engineering strategy.

Microbial Cell Factory of Baccatin III Preparation in Escherichia coli by Increasing DBAT Thermostability and in vivo Acetyl-CoA Supply

Given the rapid development of genome mining in this decade, the substrate channel of paclitaxel might be identified in the near future. A robust microbial cell factory with gene dbat, encoding a key rate-limiting enzyme 10-deacetylbaccatin III-10-O-transferase (DBAT) in paclitaxel biosynthesis to synthesize the precursor baccatin III, will lay out a promising foundation for paclitaxel de novo synthesis. Here, we integrated gene dbat into the wild-type Escherichia coli BW25113 to construct strain BWD01. Yet, it was relatively unstable in baccatin III synthesis. Mutant gene dbat S189V with improved thermostability was screened out from a semi-rational mutation library of DBAT. When it was over-expressed in an engineered strain N05 with improved acetyl-CoA generation, combined with carbon source optimization of fermentation engineering, the production level of baccatin III was significantly increased. Using this combination, integrated strain N05S01 with mutant dbat S189V achieved a 10.50-fold increase in baccatin III production compared with original strain BWD01. Our findings suggest that the combination of protein engineering and metabolic engineering will become a promising strategy for paclitaxel production.

Enhancing Thermostability of Lipase from Pseudomonas alcaligenes for producing L-menthol by the CREATE Strategy

The Lipase from Pseudomonas alcaligenes (PaL) catalyzes the hydrolysis of racemic menthol propionate to produce L-menthol, one of the most important flavoring agents in food, cosmetics and pharmaceuticals industries. However,...

Computational Redesign of an ω-Transaminase from Pseudomonas jessenii for Asymmetric Synthesis of Enantiopure Bulky Amines

ω-Transaminases (ω-TA) are attractive biocatalysts for the production of chiral amines from prochiral ketones via asymmetric synthesis. However, the substrate scope of ω-TAs is usually limited due to steric hindrance at the active site pockets. We explored a protein engineering strategy using computational design to expand the substrate scope of an (S)-selective ω-TA from Pseudomonas jessenii (PjTA-R6) toward the production of bulky amines. PjTA-R6 is attractive for use in applied biocatalysis due to its thermostability, tolerance to organic solvents, and acceptance of high concentrations of isopropylamine as amino donor. PjTA-R6 showed no detectable activity for the synthesis of six bicyclic or bulky amines targeted in this study. Six small libraries composed of 7–18 variants each were separately designed via computational methods and tested in the laboratory for ketone to amine conversion. In each library, the vast majority of the variants displayed the desired activity, and of the 40 different designs, 38 produced the target amine in good yield with >99% enantiomeric excess. This shows that the substrate scope and enantioselectivity of PjTA mutants could be predicted in silico with high accuracy. The single mutant W58G showed the best performance in the synthesis of five structurally similar bulky amines containing the indan and tetralin moieties. The best variant for the other bulky amine, 1-phenylbutylamine, was the triple mutant W58M + F86L + R417L, indicating that Trp58 is a key residue in the large binding pocket for PjTA-R6 redesign. Crystal structures of the two best variants confirmed the computationally predicted structures. The results show that computational design can be an efficient approach to rapidly expand the substrate scope of ω-TAs to produce enantiopure bulky amines.

Biocatalytic Reduction Reactions from a Chemist's Perspective

Reductions play a key role in organic synthesis, producing chiral products with new functionalities. Enzymes can catalyse such reactions with exquisite stereo-, regio- and chemoselectivity, leading the way to alternative shorter classical synthetic routes towards not only high-added-value compounds but also bulk chemicals. In this review we describe the synthetic state-of-the-art and potential of enzymes that catalyse reductions, ranging from carbonyl, enone and aromatic reductions to reductive aminations.

From thiol-subtilisin to omniligase: Design and structure of a broadly applicable peptide ligase

Omniligase-1 is a broadly applicable enzyme for peptide bond formation between an activated acyl donor peptide and a non-protected acyl acceptor peptide. The enzyme is derived from an earlier subtilisin variant called peptiligase by several rounds of protein engineering aimed at increasing synthetic yields and substrate range. To examine the contribution of individual mutations on S/H ratio and substrate scope in peptide synthesis, we selected peptiligase variant M222P/L217H as a starting enzyme and introduced successive mutations. Mutation A225N in the S1′ pocket and F189W of the S2′ pocket increased the synthesis to hydrolysis (S/H) ratio and overall coupling efficiency, whereas the I107V mutation was added to S4 pocket to increase the reaction rate. The final omniligase variants appeared to have a very broad substrate range, coupling more than 250 peptides in a 400-member library of acyl acceptors, as indicated by a high-throughput FRET assay. Crystal structures and computational modelling could rationalize the exceptional properties of omniligase-1 in peptide synthesis

Computational design of noncanonical amino acid-based thioether staples at N/C-terminal domains of multi-modular pullulanase for thermostabilization in enzyme catalysis

Enzyme thermostabilization is considered a critical and often obligatory step in biosynthesis, because thermostability is a significant property of enzymes that can be used to evaluate their feasibility for industrial applications. However, conventional strategies for thermostabilizing enzymes generally introduce non-covalent interactions and/or natural covalent bonds caused by natural amino acid substitutions, and the trade-off between the activity and stability of enzymes remains a challenge. Here, we developed a computationally guided strategy for constructing thioether staples by incorporating noncanonical amino acid (ncAA) into the more flexible N/C-terminal domains of the multi-modular pullulanase from Bacillus thermoleovorans (BtPul) to enhance its thermostability. First, potential thioether staples located in the N/C-terminal domains of BtPul were predicted using RosettaMatch. Next, eight variants involving stable thioether staples were precisely predicted using FoldX and Rosetta ddg_monomer. Six positive variants were obtained, of which T73(O2beY)-171C had a 157% longer half-life at 70 °C and an increase of 7.0 °C in T m, when compared with the wild-type (WT). T73(O2beY)-171C/T126F/A72R exhibited an even more improved thermostability, with a 211% increase in half-life at 70 °C and a 44% enhancement in enzyme activity compared with the WT, which was attributed to further optimization of the local interaction network. This work introduces and validates an efficient strategy for enhancing the thermostability and activity of multi-modular enzymes.

Asymmetric Synthesis of Optically Pure Aliphatic Amines with an Engineered Robust ω-Transaminase

The production of chiral amines by transaminase-catalyzed amination of ketones is an important application of biocatalysis in synthetic chemistry. It requires transaminases that show high enantioselectivity in asymmetric conversion of the ketone precursors. A robust derivative of ω-transaminase from Pseudomonasjessenii (PjTA-R6) that naturally acts on aliphatic substrates was constructed previously by our group. Here, we explore the catalytic potential of this thermostable enzyme for the synthesis of optically pure aliphatic amines and compare it to the well-studied transaminases from Vibrio fluvialis (VfTA) and Chromobacterium violaceum (CvTA). The product yields indicated improved performance of PjTA-R6 over the other transaminases, and in most cases, the optical purity of the produced amine was above 99% enantiomeric excess (e.e.). Structural analysis revealed that the substrate binding poses were influenced and restricted by the switching arginine and that this accounted for differences in substrate specificities. Rosetta docking calculations with external aldimine structures showed a correlation between docking scores and synthetic yields. The results show that PjTA-R6 is a promising biocatalyst for the asymmetric synthesis of aliphatic amines with a product spectrum that can be explained by its structural features.

Construction and Application of PLP Self-sufficient Biocatalysis System for Threonine Aldolase

A number of organic synthesis involve threonine aldolase (TA), a pyridoxal phosphate (PLP)-dependent enzyme. Although the addition of exogenous PLP is necessary for the reactions, it increases the cost and complicates the purification of the product. This work constructed a PLP self-sufficient biocatalysis system for TA, which included an improvement of the intracellular PLP level and co-immobilization of TA with PLP. Engineered strain BL-ST was constructed by introducing PLP synthase PdxS/T to Escherichia coli BL21(ED3). The intracellular PLP concentration of the strain increased approximately fivefold to 48.5 μmol/g_DCW. l-TA, from Bacillus nealsonii (BnLTA), was co-expressed in the strain BL-ST with PdxS/T, resulting in the engineered strain BL-BnLTA-ST. Compared with the control strain BL-BnLTA (254.1 U/L), the enzyme activity of the strain BL-BnLTA-ST reached 1518.4 U/L without the addition of exogenous PLP. An efficient co-immobilization system was then designed. The epoxy resin LX-1000HFA wrapped by polyethyleneimine (PEI) acted as a carrier to immobilize the crude enzyme solution of the strain BL-BnLTA-ST mixed with an extra 100 μM of exogenous PLP, resulting in the catalyst HFA^PEI-BnLTA-STPLP 100. HFA^PEI-BnLTA-STPLP 100 exhibited a half-life of approximately 450 h, and the application of the catalyst in the continuous biosynthesis of 3-[4-(methylsulfonyl) phenyl] serine had more than 180 batch reactions (>60%_conv) without the extra addition of exogenous PLP. The excellent compatibility and stability of the system were further confirmed by other TAs. This work introduced a PLP self-sufficient biocatalysis system that can reduce the cost of PLP and contribute to the industrial application of TA. In addition, the system may also be applied in other PLP-dependent enzymes.

Stabilization of ω-transaminase from Pseudomonas fluorescens by immobilization techniques

Transaminases are a class of enzymes with promising applications for the preparation and resolution of a vast diversity of valued amines. Their poor operational stability has fueled many investigations on its stabilization due to their biotechnological relevance. In this work, we screened the stabilization of the tetrameric ω-transaminase from Pseudomonas fluorescens (PfωTA) through both carrier-bound and carrier-free immobilization techniques. The best heterogeneous biocatalyst was the PfωTA immobilized as cross-linked enzyme aggregates (PfωTA-CLEA) which resulted after studying different parameters as the precipitant, additives and glutaraldehyde concentrations. The best conditions for maximum recovered activity (29 %) and maximum thermostability at 60 ºC and 70 ºC (100 % and 71 % residual activity after 1 h, respectively) were achieved by enzyme precipitation with 90% acetone or ethanol, in presence of BSA (100 mg/mL) and employing glutaraldehyde (100 mM) as cross-linker. Studies on different conditions for PfωTA-CLEA preparation yielded a biocatalyst that exhibited 31 and 4.6 times enhanced thermal stability at 60 °C and 70 °C, respectively, compared to its soluble counterpart. The PfωTA-CLEA was successfully used in the bioamination of 4-hydroxybenzaldehyde to 4-hydroxybenzylamine. To the best of our knowledge, this is the first report describing a transaminase cross-linked enzyme aggregates as immobilization strategy to generate a biocatalyst with outstanding thermostability.

Novel Routes in Transformation of Lignocellulosic Biomass to Furan Platform Chemicals: From Pretreatment to Enzyme Catalysis

The constant depletion of fossil fuels along with the increasing need for novel materials, necessitate the development of alternative routes for polymer synthesis. Lignocellulosic biomass, the most abundant carbon source on the planet, can serve as a renewable starting material for the design of environmentally-friendly processes for the synthesis of polyesters, polyamides and other polymers with significant value. The present review provides an overview of the main processes that have been reported throughout the literature for the production of bio-based monomers from lignocellulose, focusing on physicochemical procedures and biocatalysis. An extensive description of all different stages for the production of furans is presented, starting from physicochemical pretreatment of biomass and biocatalytic decomposition to monomeric sugars, coupled with isomerization by enzymes prior to chemical dehydration by acid Lewis catalysts. A summary of all biotransformations of furans carried out by enzymes is also described, focusing on galactose, glyoxal and aryl-alcohol oxidases, monooxygenases and transaminases for the production of oxidized derivatives and amines. The increased interest in these products in polymer chemistry can lead to a redirection of biomass valorization from second generation biofuels to chemical synthesis, by creating novel pathways to produce bio-based polymers.

Concept Study for an Integrated Reactor-Crystallizer Process for the Continuous Biocatalytic Synthesis of (S)-1-(3-Methoxyphenyl)ethylamine

An integrated biocatalysis-crystallization concept was developed for the continuous amine transaminase-catalyzed synthesis of (S)-1-(3-methoxyphenyl)ethylamine, which is a valuable intermediate for the synthesis of rivastigmine, a highly potent drug for the treatment of early stage Alzheimer’s disease. The three-part vessel system developed for this purpose consists of a membrane reactor for the continuous synthesis of the product amine, a saturator vessel for the continuous supply of the amine donor isopropylammonium and the precipitating reagent 3,3-diphenylpropionate and a crystallizer in which the product amine can continuously precipitate as (S)-1-(3-methoxyphenyl)ethylammonium-3,3-diphenylpropionate.

Approaching boiling point stability of an alcohol dehydrogenase through computationally-guided enzyme engineering

Enzyme instability is an important limitation for the investigation and application of enzymes. Therefore, methods to rapidly and effectively improve enzyme stability are highly appealing. In this study we applied a computational method (FRESCO) to guide the engineering of an alcohol dehydrogenase. Of the 177 selected mutations, 25 mutations brought about a significant increase in apparent melting temperature (ΔTm ≥ +3 °C). By combining mutations, a 10-fold mutant was generated with a Tm of 94 °C (+51 °C relative to wild type), almost reaching water's boiling point, and the highest increase with FRESCO to date. The 10-fold mutant's structure was elucidated, which enabled the identification of an activity-impairing mutation. After reverting this mutation, the enzyme showed no loss in activity compared to wild type, while displaying a Tm of 88 °C (+45 °C relative to wild type). This work demonstrates the value of enzyme stabilization through computational library design.