Variations in the stability of NCR ene reductase by rational enzyme loop modulation

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.

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.

Combining Flexible Region Design and Automatic Design to Enhance the Thermal Stability and Catalytic Efficiency of Leucine Dehydrogenase

Leucine dehydrogenase (LeuDH, EC 1.4.1.9) can reversibly catalyze the oxidative deamination of l-leucine and some other specific α-amino acids to form the corresponding α-ketoacids. This reaction has great significance in the field of food additives and the pharmaceutical industry. The LeuDH from Exiguobacterium sibiricum (EsLeuDH) has high catalytic efficiency but limited thermal stability, hindering its widespread industrial application. In this study, a mutant N5F/I12L/A352Y of EsLeuDH (referred to as M2) was developed with enhanced thermal stability and catalytic activity through rational modification. The M2 mutant exhibits a half-life at 60 °C (t1/2(60 °C)) of 975.7 min and a specific activity of 69.6 U mg-1, which is 5.4 and 2.1 times higher than those of EsLeuDH, respectively. This research may facilitate the utilization of EsLeuDH at elevated temperatures, enhancing its potential for industrial applications. The findings offer a practical and efficient approach for optimizing LeuDH and other industrial enzymes.

N-terminal truncation (N-) and directional proton transfer in an old yellow enzyme enables tunable efficient producing (R)- or (S)-citronellal

(R)-Citronellal is a valuable molecule as the precursor for the industrial synthesis of (-)-menthol, one of the worldwide best-selling compounds in the flavors and fragrances field. However, its biocatalytic production, even from the optically pure substrate (E)-citral, is inherently limited by the activity of Old Yellow Enzyme (OYE). Herein, we rationally designed a different approach to increase the activity of OYE in biocatalytic production. The activity of OYE from Corynebacterium glutamicum (CgOYE) is increased, as well as superior thermal stability and pH tolerance via truncating the different lengths of regions at N-terminal of CgOYE. Next, we converted the truncation mutant N31-CgOYE, a protein involved in proton transfer for the asymmetric hydrogenation of CC bonds, into highly (R)- and (S)-stereoselective enzymes using only three mutations. The mixture of racemic (E/Z)-citral is reduced into the (R)-citronellal with ee and conversion up to 99 % by the mutant of CgOYE, overcoming the problem of the reduction for the mixtures of (E/Z)-citral in biocatalytic reaction. The present work provides a general and effective strategy for improving the activity of OYE, in which the partially conserved histidine residues provide "tunable gating" for the enantioselectivity for both the (R)- and (S)-isomerases.

Targeted Mutation of a Non-catalytic Gating Residue Increases the Rate of Pseudomonas aeruginosa d-Arginine Dehydrogenase Catalytic Turnover

Commercial food and l-amino acid industries rely on bioengineered d-amino acid oxidizing enzymes to detect and remove d-amino acid contaminants. However, the bioengineering of enzymes to generate faster biological catalysts has proven difficult as a result of the failure to target specific kinetic steps that limit enzyme turnover, kcat, and the poor understanding of loop dynamics critical for catalysis. Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) oxidizes most d-amino acids and is a good candidate for application in the l-amino acid and food industries. The side chain of the loop L2 E246 residue located at the entrance of the PaDADH active site pocket potentially favors the closed active site conformation and secures the substrate upon binding. This study used site-directed mutagenesis, steady-state, and rapid reaction kinetics to generate the glutamine, glycine, and leucine variants and investigate whether increasing the rate of product release could translate to an increased enzyme turnover rate. Upon E246 mutation to glycine, there was an increased rate of d-arginine turnover kcat from 122 to 500 s-1. Likewise, the kcat values increased 2-fold for the glutamine or leucine variants. Thus, we have engineered a faster biocatalyst for industrial applications by selectively increasing the rate of the PaDADH product release.

Semi-rational protein engineering of a novel ene-reductase from Galdieria sulphuraria for asymmetric reduction of (R)-carvone and ketoisophorone

Asymmetric reduction of (R)-carvone and ketoisophorone by an engineered ene-reductase from Galdieria sulphuraria (GsOYE) combined with glucose dehydrogenase for NADPH regeneration were studied. A semi-rational protein engineering was used to enhance the activity and selectivity of GsOYE. Upon the sequence alignment and molecular docking results, two amino acid residues at positions 66 and 270 were selected as saturation mutation sites. Finally, a single substitution variant of GsOYE-N270A with complete conversion (100%) and diastereoselectivity (de_p >99%) for reduction of (R)-carvone and a double substitution variant GsOYE-Y66P/N270H with improved stereoselectivity for reduction of ketoisophorone were obtained. This article is protected by copyright. All rights reserved.

Insertions and deletions in protein evolution and engineering

Protein evolution or engineering studies are traditionally focused on amino acid substitutions and the way these contribute to fitness. Meanwhile, the insertion and deletion of amino acids is often overlooked, despite being one of the most common sources of genetic variation. Recent methodological advances and successful engineering stories have demonstrated that the time is ripe for greater emphasis on these mutations and their understudied effects. This review highlights the evolutionary importance and biotechnological relevance of insertions and deletions (indels). We provide a comprehensive overview of approaches that can be employed to include indels in random, (semi)-rational or computational protein engineering pipelines. Furthermore, we discuss the tolerance to indels at the structural level, address how domain indels can link the function of unrelated proteins, and feature studies that illustrate the surprising and intriguing potential of frameshift mutations.

A combination of bioinformatics analysis and rational design strategies to enhance keratinase thermostability for efficient biodegradation of feathers

Keratinase has shown great significance and application potentials in the biodegradation and recycle of keratin waste due to its unique and efficient hydrolysis ability. However, the inherent instability of the enzyme limits its practical utilization. Herein, we obtained a thermostability-enhanced keratinase based on a combination of bioinformatics analysis and rational design strategies for the efficient biodegradation of feathers. A systematical in silico analysis combined with filtering of virtual libraries derived a smart library for experimental validation. Synergistic mutations around the highly flexible loop, the calcium binding site and the non-consensus amino acids generated a dominant mutant which increased the optimal temperature of keratinase from 40 °C to 60 °C, and the half-life at 60 °C was increased from 17.3 min to 66.1 min. The mutant could achieve more than 66% biodegradation of 50 g/L feathers to high-valued keratin product with a major molecular weight of 36 kDa. Collectively, this work provided a promising keratinase variant with enhanced thermostability for efficient conversion of keratin wastes to valuable products. It also generated a general strategy to facilitate enzyme thermostability design which is more targeted and predictable.

Active-site loop variations adjust activity and selectivity of the cumene dioxygenase

Active-site loops play essential roles in various catalytically important enzyme properties like activity, selectivity, and substrate scope. However, their high flexibility and diversity makes them challenging to incorporate into rational enzyme engineering strategies. Here, we report the engineering of hot-spots in loops of the cumene dioxygenase from Pseudomonas fluorescens IP01 with high impact on activity, regio- and enantioselectivity. Libraries based on alanine scan, sequence alignments, and deletions along with a novel insertion approach result in up to 16-fold increases in activity and the formation of novel products and enantiomers. CAVER analysis suggests possible increases in the active pocket volume and formation of new active-site tunnels, suggesting additional degrees of freedom of the substrate in the pocket. The combination of identified hot-spots with the Linker In Loop Insertion approach proves to be a valuable addition to future loop engineering approaches for enhanced biocatalysts.

Assembly of a Rieske non-heme iron oxygenase multicomponent system from Phenylobacterium immobile E DSM 1986 enables pyrazon cis-dihydroxylation in E. coli

Abstract

Phenylobacterium immobile strain E is a soil bacterium with a striking metabolism relying on xenobiotics, such as the herbicide pyrazon, as sole carbon source instead of more bioavailable molecules. Pyrazon is a heterocyclic aromatic compound of environmental concern and its biodegradation pathway has only been reported in P. immobile. The multicomponent pyrazon oxygenase (PPO), a Rieske non-heme iron oxygenase, incorporates molecular oxygen at the 2,3 position of the pyrazon phenyl moiety as first step of degradation, generating a cis-dihydrodiendiol. The aim of this work was to identify the genes encoding for each one of the PPO components and enable their functional assembly in Escherichia coli. P. immobile strain E genome sequencing revealed genes encoding for RO components, such as ferredoxin-, reductase-, α- and β-subunits of an oxygenase. Though, P. immobile E displays three prominent differences with respect to the ROs currently characterized: (1) an operon-like organization for PPO is absent, (2) all the elements are randomly scattered in its DNA, (3) not only one, but 19 different α-subunits are encoded in its genome. Herein, we report the identification of the PPO components involved in pyrazon cis-dihydroxylation in P. immobile, its appropriate assembly, and its functional reconstitution in E. coli. Our results contributes with the essential missing pieces to complete the overall elucidation of the PPO from P. immobile.

Key points

• Phenylobacterium immobile E DSM 1986 harbors the only described pyrazon oxygenase (PPO).

• We elucidated the genes encoding for all PPO components.

• Heterologous expression of PPO enabled pyrazon dihydroxylation in E. coli JW5510.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-021-11129-w.

Old yellow enzymes: structures and structure-guided engineering for stereocomplementary bioreduction

Since the first discovery of old yellow enzyme 1 (OYE1) from Saccharomyces pastorianus in 1932, biocatalytic asymmetric reduction of activated alkenes by OYEs has become a valuable reaction in organic synthesis. To access stereocomplementary C=C-bond bioreduction, the mining of novel OYEs and especially the protein engineering of existing OYEs have been performed, which successfully achieved the stereocomplementary reduction in several cases and further raise the potential of applications. In this review, we analyzed the structures, active sites, and substrate recognition of OYEs, which are the bases for their substrate specificity and stereospecificity. Sequence similarity network of OYEs superfamily was also constructed to investigate the scope of characterized OYEs. The structure-guided engineering to switch the stereoselectivity of OYEs and thus access stereocomplementary bioreduction over the last decade (2009-2020) was then reviewed and discussed, which might give new insights into the mining and engineering of related biocatalysts. KEY POINTS: • The sequence similarity network of OYEs superfamily was constructed and annotated. • The structures and active sites of OYEs from different classes were compared. • "Left/right" binding mode was used to explain the stereopreferences of OYEs. • Structure-guided engineering of OYEs to switch their stereoselectivity was reviewed.

Discovery, Characterisation, Engineering and Applications of Ene Reductases for Industrial Biocatalysis

Recent studies of multiple enzyme families collectively referred to as ene-reductases (ERs) have highlighted potential industrial application of these biocatalysts in the production of fine and speciality chemicals. Processes have been developed whereby ERs contribute to synthetic routes as isolated enzymes, components of multi-enzyme cascades, and more recently in metabolic engineering and synthetic biology programmes using microbial cell factories to support chemicals production. The discovery of ERs from previously untapped sources and the expansion of directed evolution screening programmes, coupled to deeper mechanistic understanding of ER reactions, have driven their use in natural product and chemicals synthesis. Here we review developments, challenges and opportunities for the use of ERs in fine and speciality chemicals manufacture. The ER research field is rapidly expanding and the focus of this review is on developments that have emerged predominantly over the last 4 years.

Computation-Aided Rational Deletion of C-Terminal Region Improved the Stability, Activity, and Expression Level of GH2 β-Glucuronidase

In this study, computation-aided design on the basis of structural analysis was employed to rationally identify a highly dynamic C-terminal region that regulates the stability, expression level, and activity of a GH2 fungal glucuronidase from Aspergillus oryzae Li-3 (PGUS). Then, four mutants with a precisely truncated C-terminal region in different lengths were constructed; among them, mutant D591-604 with a 3.8-fold increase in half-life at 65 °C and a 6.8 kJ/mol increase in Gibbs free energy showed obviously improved kinetic and thermodynamic stability in comparison to PGUS. Mutants D590-604 and D591-604 both showed approximately 2.4-fold increases in the catalytic efficiency k_cat/ K_m and 1.8-fold increases in the expression level. Additionally, the expression level of PGUS was doubled through a C-terminal region swap with bacterial GUS from E. coli (EGUS). Finally, the robust PGUS mutants D590-604 and D591-604 were applied in the preparation of glycyrrhetinic acid with 4.0- and 4.4-fold increases in concentration through glycyrrhizin hydrolysis by a fed-batch process.

Enhanced trypsin thermostability in Pichia pastoris through truncating the flexible region

Background

High thermostability is required for trypsin to have wider industrial applications. Target mutagenesis at flexible regions has been proved to be an efficient protein engineering method to enhance the protein thermostability.

Results

The flexible regions in porcine trypsin were predicted using the methods including molecular dynamic simulation, FlexPred, and FoldUnfold. The amino acids 78–90 was predicted to be the highly flexible region simultaneously by the three methods and hence selected to be the mutation target. We constructed five variants (D3, D5, D7, D9, and D11) by truncating the region. And the variant D9 showed higher thermostability, with a 5 °C increase in T_opt, 5.8 °C rise in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{50}^{10}$$\end{document}T5010, and a 4.5 °C rise in T_m, compared to the wild-type. Moreover, the half-life value of the variant D9 was also found to be dramatically improved by 46 min. Circular dichroism and intrinsic fluorescence indicated that the structures had no significant change between the variant D9 and the wild-type. The surface hydrophobicity of D9 was measured to be lower than that of wild-type, indicating the increased hydrophobic interaction, which could have contributed to the improved thermostability of D9.

Conclusions

These results showed that the thermostability of variant D9 was increased. The variant D9 could be expected to be a promising tool enzyme for its wider industrial applications. The method of truncating the flexible region used in our study has the potential to be used for enhancing the thermostability of other proteins.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1012-x) contains supplementary material, which is available to authorized users.

Role of an N-terminal extension in stability and catalytic activity of a hyperthermostable α/β hydrolase fold esterase

The carbohydrate esterase family 7 (CE7) enzymes catalyze the deacetylation of acetyl esters of a broad range of alcohols and is unique in its activity towards cephalosporin C. The CE7 fold contains a conserved N-terminal extension that distinguishes it from the canonical α/β hydrolase fold. The hexameric quaternary structure indicates that the N-terminus may affect activity and specificity by controlling access of substrates to the buried active sites via an entrance tunnel. In this context, we characterized the catalytic parameters, conformation and thermal stability of two truncation variants lacking four and ten residues of the N-terminal region of the hyperthermostable Thermotoga maritima CE7 acetyl esterase (TmAcE). The truncations did not affect the secondary structure or the fold but modulated the oligomerization dynamics. A modest increase was observed in substrate specificity for acetylated xylose compared with acetylated glucose. A drastic reduction of ~30-40°C in the optimum temperature for activity of the variants indicated lower thermal stability. The loss of hyperthermostability appears to be an indirect effect associated with an increase in the conformational flexibility of an otherwise rigid neighboring loop containing a catalytic triad residue. The results suggest that the N-terminal extension was evolutionarily selected to preserve the stability of the enzyme.

Structural investigation into the C-terminal extension of the ene-reductase from Ralstonia (Cupriavidus) metallidurans

Ene-reductases (ERs), or Old Yellow Enzymes, catalyze the asymmetric reduction of various activated alkenes. This class of biocatalysts is considered an attractive alternative to current chemical technologies for hydrogenation due to their high selectivity and specificity. Here the X-ray crystal structure of RmER, a "thermophilic"-like ER from Ralstonia (Cupriavidus) metallidurans, is reported. Unlike other members of this class of ERs, RmER is monomeric in solution which we previously related to its atypical elongated C-terminus. A typical dimer interface was however observed in our crystal structure, with the conserved Arg-"finger" forming part of the adjacent monomer's active site and the elongated C-terminus extending into the active site through contacting the "capping" domain. This dimerization also resulted in the loss of one FMN cofactor from each dimer pair. This potential transient dimerization and dissociation of FMN could conceivably explain the rapid rates previously observed when an FMN light-driven cofactor regeneration system was used during catalysis with RmER.

Enantioselective Reduction of Citral Isomers in NCR Ene Reductase: Analysis of an Active-Site Mutant Library

A deeper understanding of the >99 % S-selective reduction of both isomers of citral catalyzed by NCR ene reductase was achieved by active-site mutational studies and docking simulation. Though structurally similar, the E/Z isomers of citral showed a significantly varying selectivity response to introduced mutations. Although it was possible to invert (E)-citral reduction enantioselectivity to ee 46 % (R) by introducing mutation W66A, for (Z)-citral it remained ≥88 % (S) for all single-residue variants. Residue 66 seems to act as a lever for opposite binding modes. This was underlined by a W66A-based double-mutant library that enhanced the (E)-citral derived enantioselectivity to 63 % (R) and significantly lowered the S selectivity for (Z)-citral to 44 % (S). Formation of (R)-citronellal from an (E/Z)-citral mixture is a desire in industrial (-)-menthol synthesis. Our findings pave the way for a rational enzyme engineering solution.

Molecular and catalytic properties of fungal extracellular cellobiose dehydrogenase produced in prokaryotic and eukaryotic expression systems

BACKGROUND

Cellobiose dehydrogenase (CDH) is an extracellular enzyme produced by lignocellulolytic fungi. cdh gene expression is high in cellulose containing media, but relatively low CDH concentrations are found in the supernatant of fungal cultures due to strong binding to cellulose. Therefore, heterologous expression of CDH in Pichia pastoris was employed in the last 15 years, but the obtained enzymes were over glycosylated and had a reduced specific activity.

RESULTS

We compare the well-established CDH expression host P. pastoris with the less frequently used hosts Escherichia coli, Aspergillus niger, and Trichoderma reesei. The study evaluates the produced quantity and protein homogeneity of Corynascus thermophilus CDH in the culture supernatants, the purification, and finally compares the enzymes in regard to cofactor loading, glycosylation, catalytic constants and thermostability.

CONCLUSIONS

Whereas E. coli could only express the catalytic dehydrogenase domain of CDH, all eukaryotic hosts could express full length CDH including the cytochrome domain. The CDH produced by T. reesei was most similar to the CDH originally isolated from the fungus C. thermophilus in regard to glycosylation, cofactor loading and catalytic constants. Under the tested experimental conditions the fungal expression hosts produce CDH of superior quality and uniformity compared to P. pastoris.

Applications of protein engineering to members of the old yellow enzyme family

In the 20 years since Massey's initial report in 1995, interest in using alkene reductases to prepare chiral intermediates for synthesis has grown rapidly. While native alkene reductases often show very high stereoselectivities toward favorable substrates, these enzymes have somewhat size-restricted active sites that limit their substrate ranges to small alkenes. In addition, most alkene reductases have the same stereoselectivities, which makes it difficult to access the "other" product enantiomers. Protein engineering strategies have been used to address both of these issues and good progress has been made in several cases. This review summarizes published examples through late 2014 and focuses on studies of six enzymes: Saccharomyces pastorianus OYE 1, tomato OPR1, Zymomonas mobilis NCR, Enterobacter cloacae PB2 PETN reductase, Bacillus subtilis YqjM and Pichia stipitis OYE 2.6.

STRUCTURAL AND FUNCTIONAL CONSEQUENCES OF CIRCULAR PERMUTATION ON THE ACTIVE SITE OF OLD YELLOW ENZYME

Circular permutation of the NADPH-dependent oxidoreductase Old Yellow Enzyme from Saccharomyces pastorianus (OYE1) can significantly enhance the enzyme’s catalytic performance. Termini relocation into four regions of the protein (sectors I–IV) near the active site has proven effective in altering enzyme function. To better understand the structural consequences and rationalize the observed functional gains in these OYE1 variants, we selected representatives from sectors I–III for further characterization by biophysical methods and X-ray crystallography. These investigations not only show trends in enzyme stability and quaternary structure as a function of termini location but also provide a possible explanation for the catalytic gains in our top-performing OYE variant (new N-terminus at residue 303; sector III). Crystallographic analysis indicates that termini relocation into sector III affects the loop β6 region (amino acid positions: 290–310) of OYE1, which forms a lid over the active site. Peptide backbone cleavage greatly enhances local flexibility, effectively converting the loop into a tether and consequently increasing the environmental exposure of the active site. Interestingly, such an active site remodeling does not negatively impact the enzyme’s activity and stereoselectivity; neither does it perturb the conformation of other key active site residues with the exception of Y375. These observations were confirmed in truncation experiments, deleting all residues of the loop β6 region in our OYE variant. Intrigued by the finding that circular permutation leaves most of the key catalytic residues unchanged, we also tested OYE permutants for possible additive or synergistic effects of amino acid substitutions. Distinct functional changes in these OYE variants were detected upon mutations at W116, known in native OYE1 to cause inversion of diastereoselectivity for (S)-carvone reduction. Our findings demonstrate the contribution of loop β6 toward determining the stereoselectivity of OYE1, an important insight for future OYE engineering efforts.

Pichia stipitis OYE 2.6 variants with improved catalytic efficiencies from site-saturation mutagenesis libraries

An earlier directed evolution project using alkene reductase OYE 2.6 from Pichia stipitis yielded 13 active site variants with improved properties toward three homologous Baylis-Hillman adducts. Here, we probed the generality of these improvements by testing the wild-type and all 13 variants against a panel of 16 structurally-diverse electron-deficient alkenes. Several substrates were sterically demanding, and as hoped, creating additional active site volume yielded better conversions for these alkenes. The most impressive improvement was found for 2-butylidenecyclohexanone. The wild-type provided less than 20% conversion after 24h; a triple mutant afforded more than 60% conversion in the same time period. Moreover, even wild-type OYE 2.6 can reduce cyclohexenones with very bulky 4-substituents efficiently.

Computational library design for increasing haloalkane dehalogenase stability

We explored the use of a computational design framework for the stabilization of the haloalkane dehalogenase LinB. Energy calculations, disulfide bond design, molecular dynamics simulations, and rational inspection of mutant structures predicted many stabilizing mutations. Screening of these in small mutant libraries led to the discovery of seventeen point mutations and one disulfide bond that enhanced thermostability. Mutations located in or contacting flexible regions of the protein had a larger stabilizing effect than mutations outside such regions. The combined introduction of twelve stabilizing mutations resulted in a LinB mutant with a 23 °C increase in apparent melting temperature (Tm,app , 72.5 °C) and an over 200-fold longer half-life at 60 °C. The most stable LinB variants also displayed increased compatibility with co-solvents, thus allowing substrate conversion and kinetic resolution at much higher concentrations than with the wild-type enzyme.