Co-evolution of activity and thermostability of an aldo-keto reductase KmAKR for asymmetric synthesis of statin precursor dichiral diols

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.

Structural-guided design to improve the catalytic performance of aldo-keto reductase KdAKR

Aldo-keto reductases (AKRs) are important biocatalysts that can be used to synthesize chiral pharmaceutical alcohols. In this study, the catalytic activity and stereoselectivity of a NADPH-dependent AKR from Kluyveromyces dobzhanskii (KdAKR) toward t-butyl 6-chloro (5S)-hydroxy-3-oxohexanoate ((5S)-CHOH) were improved by mutating its residues in the loop regions around the substrate-binding pocket. And the thermostability of KdAKR was improved by a consensus sequence method targeted on the flexible regions. The best mutant M6 (Y28A/L58I/I63L/G223P/Y296W/W297H) exhibited a 67-fold higher catalytic efficiency compared to the wild-type (WT) KdAKR, and improved R-selectivity toward (5S)-CHOH (dep value from 47.6% to >99.5%). Moreover, M6 exhibited a 6.3-fold increase in half-life (t1/2 ) at 40°C compared to WT. Under the optimal conditions, M6 completely converted 200 g/L (5S)-CHOH to diastereomeric pure t-butyl 6-chloro-(3R, 5S)-dihydroxyhexanoate ((3R, 5S)-CDHH) within 8.0 h, with a space-time yield of 300.7 g/L/day. Our results deepen the understandings of the structure-function relationship of AKRs, providing a certain guidance for the modification of other AKRs.

Reductive enzymatic dynamic kinetic resolution affording 115 g/L (S)-2-phenylpropanol

Background

Published biocatalytic routes for accessing enantiopure 2-phenylpropanol using oxidoreductases afforded maximal product titers of only 80 mM. Enzyme deactivation was identified as the major limitation and was attributed to adduct formation of the aldehyde substrate with amino acid residues of the reductase.

Results

A single point mutant of Candida tenuis xylose reductase (CtXR D51A) with very high catalytic efficiency (43·10³ s⁻¹ M⁻¹) for (S)-2-phenylpropanal was found. The enzyme showed high enantioselectivity for the (S)-enantiomer but was deactivated by 0.5 mM substrate within 2 h. A whole-cell biocatalyst expressing the engineered reductase and a yeast formate dehydrogenase for NADH-recycling provided substantial stabilization of the reductase. The relatively slow in situ racemization of 2-phenylpropanal and the still limited biocatalyst stability required a subtle adjustment of the substrate-to-catalyst ratio. A value of 3.4 g_substrate/g_{cell-dry-weight} was selected as a suitable compromise between product ee and the conversion ratio. A catalyst loading of 40 g_{cell-dry-weight} was used to convert 1 M racemic 2-phenylpropanal into 843 mM (115 g/L) (S)-phenylpropanol with 93.1% ee.

Conclusion

The current industrial production of profenols mainly relies on hydrolases. The bioreduction route established here represents an alternative method for the production of profenols that is competitive with hydrolase-catalyzed kinetic resolutions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00715-5.

Semirational engineering of an aldo-keto reductase KmAKR for overcoming trade-offs between catalytic activity and thermostability

Enzyme engineering usually generates trade-offs between activity, stability, and selectivity. Herein, we report semirational engineering of an aldo-keto reductase (AKR) KmAKR for simultaneously enhancing its thermostability and catalytic activity. Previously, we constructed KmAKR_M9 (W297H/Y296W/K29H/Y28A/T63M/A30P/T302S/N109K/S196C), which showed outstanding activity towards t-butyl 6-chloro-(3R,5S)-dihydroxyhexanoate ((3R,5S)-CDHH), and t-butyl 6-cyano-(3R,5R)-dihydroxyhexanoate, the key chiral building blocks of rosuvastatin and atorvastatin. Under the guidance of computer-aided design including consensus residues analysis and molecular dynamics (MD) simulations, K164, S182, S232, and Q266 were dug out for their thermostability conferring roles, generating the "best" mutant KmAKR_M13 (W297H/Y296W/K29H/Y28A/T63M/A30P/T302S/N109K/S196C/K164E/S232A/S182H/Q266D). The T_m and T₅₀ ¹⁵ values of KmAKR_M13 were 10.4 and 6.1°C higher than that of KmAKR_M9 , respectively. Moreover, it displayed a significantly elevated organic solvent tolerance over KmAKR_M9 . Structural analysis indicated that stabilization of the α-helixes mainly contributed to thermostability enhancement. Under the optimized conditions, KmAKR_M13 completely asymmetrically reduced 400 g/l t-butyl 6-chloro-(5S)-hydroxy-3-oxohexanoate ((5S)-CHOH) in 8.0 h at a high substrate to catalyst ratio (S/C) of 106.7 g/g, giving diastereomerically pure (3R,5S)-CDHH (>99.5% d.e._P ) with a space-time yield (STY) of 449.2 g/l·d.

Fluorescence-based screening for engineered aldo-keto reductase KmAKR with improved catalytic performance and extended substrate scope

BACKGROUND

Aldo-keto reductases-catalyzed transformations of ketones to chiral alcohols have become an established biocatalytic process step in the pharmaceutical industry. Previously, we have discovered an aldo-keto reductase (AKR) from Kluyveromyces marxianus that is active to the aliphatic tert-butyl 6-substituted (5R/S)-hydroxy-3-oxohexanoates, but it is inactive to aromatic ketones. In order to acquire an excellent KmAKRmutant for ensuring the simultaneous improvement of activity-thermostability toward tert-butyl 6-cyano-(5R)-hydroxy-3-oxohexanoate ((5R)-1) and broadening the universal application prospects toward more substrates covering both aliphatic and aromatic ketones, a fluorescence-based high-throughput (HT) screening technique was established.

MAIN METHODS AND MAJOR RESULTS

The directed evolution of KmAKR variant M5 (KmAKR-W297H/Y296W/K29H/Y28A/T63M) produced the "best" variant M5-Q213A/T23V. It exhibited enhanced activity-thermostability toward (5R)-1, improved activity toward all 18 test substrates and strict R-stereoselectivity toward 10 substrates in comparison to M5. The enhancement of enzymatic activity and the extension of substrate scope covering aromatic ketones are proposed to be largely attributed to pushing the binding pocket of M5-Q213A/T23V to the enzyme surface, decreasing the steric hindrance at the entrance and enhancing the flexibility of loops surrounding the active center. In addition, combined with 0.94 g dry cell weight (DCW) L^-1 glucose dehydrogenase from Exiguobacterium sibiricum (EsGDH) for NADPH regeneration, 2.81 g DCW L^-1 M5-Q213A/T23V completely converted (5R)-1 of up to 450 g L^-1 at 120 g g^-1 substrates/catalysts (S/C), yielding the corresponding optically pure tert-butyl 6-cyano-(3R,5R)-dihydroxyhexanoate ((3R,5R)-2, > 99.5% d.e._p ) with a space-time yield (STY) of 1.08 kg L^-1 day^-1 .

CONCLUSIONS

A fluorescence-based HT screening system was developed to tailor KmAKR's activity, thermostability and substrate scope. The "best" variant M5-Q213A/T23V holds great potential application for the synthesis of aliphatic/aromatic R-configuration alcohols.

Tailoring an aldo-keto reductase KmAKR for robust thermostability and catalytic efficiency by stepwise evolution and structure-guided consensus engineering

t-Butyl 6-cyano-(3R,5R)-dihydroxyhexanoate ((3R,5R)-2) is an advanced chiral diol intermediate of the cholesterol-lowering drug atorvastatin. KmAKR_M5 (W297H/Y296W/K29H/Y28A/T63M) constructed in our previous work, displayed good biocatalytic performance on (3R,5R)-2. In the present work, stepwise evolution was applied to further enhance the thermostability and activity of KmAKR_M5. For thermostability enhancement, N109 and S196 located far from the active site were picked out by structure-guided consensus engineering, and mutated by site-directed mutagenesis (SDM). For catalytic efficiency improvement, the residues A30 and T302 adjacent to the substrate-binding pocket were subjected to site-saturation mutagenesis (SSM). As a result, the "best" mutant KmAKR_M9 (W297H/Y296W/K29H/Y28A/T63M/A30P/T302S/N109K/S196C) was developed, of which T₅₀¹⁵ and T_m were 5.0 °C and 8.2 °C higher than those of KmAKR_M5. Moreover, compared to KmAKR_M5, KmAKR_M9 displayed a 1.9-fold (846 vs 2436 min) and 6.7-fold (126 vs 972 min) longer half-lives at 40 and 50 °C, respectively. Structural analysis suggested that beneficial mutations introduced additional hydrophobic interactions and hydrogen bonds, contributing rigidification of the flexible loops and the increase of internal forces, hence increasing the thermostability and activity. 5 g DCW (dry cell weight) L^-1KmAKR_M9 completely reduced 350 g L^-1t-butyl 6-cyano-(5R)-hydroxy-3-oxo-hexanoate ((5R)-1), within 3.7 h at 40 °C, yielding optically pure (3R,5R)-2 (d.e._p > 99.5%) with a space-time yield (STY) of 1.82 kg L^-1 d^-1. Hence, KmAKR_M9 is a robust biocatalyst for the synthesis of (3R,5R)-2.