Martina
Sudar
a,
Zvjezdana
Findrik
*a,
Đurđa
Vasić-Rački
a,
Anna
Soler
b and
Pere
Clapés
b
aUniversity of Zagreb, Faculty of Chemical Engineering and Technology, Savska c. 16, HR-10000 Zagreb, Croatia. E-mail: zfindrik@fkit.hr; Fax: +385 1 4597 133; Tel: +385 1 4597 157
bInstitute of Advanced Chemistry of Catalonia, Biotransformation and Bioactive Molecules Group, IQAC-CSIC, Jordi Girona 18-26, 08034 Barcelona, Spain
First published on 11th August 2015
A novel cascade reaction for the production of aldol adduct (3S,4R)-6-[(benzyloxycarbonyl)amino]-5,6-dideoxyhex-2-ulose was studied in this work. The strategy combines three enzymes in one pot: (i) horse liver alcohol dehydrogenase for the oxidation of N-Cbz-3-aminopropanol to the corresponding aldehyde, (ii) NADH oxidase for the regeneration of coenzyme NAD+ and (iii) D-fructose-6-phosphate aldolase from E. coli A129S variant for the aldol addition of dihydroxyacetone to N-Cbz-3-aminopropanal. On the basis of preliminary experiments, optimization of the initial reaction conditions was done using statistical methods, i.e. factorial design of experiments. 79% yield of aldol adduct was achieved in the batch reactor after optimization.
In this work, a novel enzymatic cascade consisting of amino alcohol oxidation to amino aldehyde and aldol addition reaction was developed for the synthesis of the aldol precursor (3, Scheme 1) of D-fagomine. D-Fagomine is a naturally occurring iminosugar which can be found in plants in small quantities.11,12 It is a known inhibitor of glycoprocessing enzymes13–15 which reduces the risk of developing insulin resistance and obesity16 and may promote the adhesion of beneficial bacteria (e.g. Lactobacillus).17 It can be used as dietary supplement or a component of functional food.18 The synthesis of this compound via chemo-enzymatic pathway is a prime example of the efficiency of such method in comparison to purely chemical process.19
![]() | ||
Scheme 1 Reaction scheme of cascade reaction of amino alcohol 1 oxidation to 2 with coenzyme regeneration using NOX and aldol addition of DHA to 2 catalyzed by FSA A129S. |
The strategy devised in this work consisted of the simultaneous use of three enzymes: horse liver alcohol dehydrogenase (HLADH) used for the oxidation of N-Cbz-3-aminopropanol (1) to N-Cbz-3-aminopropanal (2), D-fructose-6-phosphate aldolase (FSA) variant A129S for the aldol addition, and a third enzyme; NADH oxidase (NOX) for coenzyme regeneration (Scheme 1). The product of this cascade reaction is (3S,4R)-6-[(benzyloxycarbonyl)amino]-5,6-dideoxyhex-2-ulose (3, Scheme 1), a precursor of D-fagomine. The aldol addition of dihydroxyacetone (DHA) to 2 catalyzed by FSA to furnish 3, was studied in detail in our previous works.20,21 Since the amino aldehyde 2 is unstable, the purpose of the cascade reaction was to produce it in situ by oxidation of 1. Similar amino alcohol was oxidized by different authors in reactions catalyzed by chloroperoxidase with the addition of tert-butyl hydroperoxide or hydrogen peroxide,22 by laccase/O2/TEMPO23 and alcohol oxidase.23 These systems avoid the coenzyme addition, but all of them are unselective towards the substrate, and result in over-oxidation of the part of amino alcohol to amino acid. We propose the use of horse liver alcohol dehydrogenase, which was successfully tested for oxidation of different amino alcohols,24 and the use of NADH oxidase (NOX) for coenzyme regeneration. NOX is a flavoprotein which oxidizes NADH to NAD+, and was isolated and purified from Lactococcus lactis.25,26 To the best of our knowledge this kind of system (Scheme 1) was not investigated before. Amino alcohol oxidation combined with aldol addition was studied by Mifsud and co-authors23 in one pot – two step consecutive reactions because of the incompatible reaction conditions for both biocatalysts. In the present approach, the reaction was carried out in one-pot with simultaneous addition of three enzymes. This system is a highly promising alternative to the one previously mentioned, since coenzyme regeneration is elegantly resolved by the addition of NADH oxidase, and the problem of thermodynamic equilibrium of oxidation25,27,28 catalyzed by ADH could be overcome. The studied process uses environmentally friendly technology and offers a novel approach to the subject. It spares the enzymes of negative effect of peroxides used in previous works22,23 and of additional chemicals such as TEMPO which pollute the product and complicate its purification. Coenzyme regeneration by NADH oxidase does not add any new compounds to the system, and requires only a minimal amount of oxygen in the solution for the maximum enzyme activity.25
Protocols for purification of 3 and the NMR analysis were identical to those described in the literature.15
Reaction yield of the main product (3) and the by-product (4) were calculated according to eqn (1).
![]() | (1) |
Enzyme operational stability decay rate was described by the kinetics of the first order20 according to eqn (2).
![]() | (2) |
Enzyme operational stability decay rate constants were estimated by non-linear regression methods (simplex and least squares fit) implemented in SCIENTIST software32 using the obtained experimental data. They were estimated by fitting the model to the experimental data – the change in specific activity of enzyme during the experiment. The calculated data were compared with the experimental data, recalculated in the optimization routine and fitted again until a minimal error between experimental and integrated values was achieved. The residual was defined as the sum of squares of the differences between the experimental and calculated data.
Variable name | Unit | Lower limit | Center point | Upper limit |
---|---|---|---|---|
A: NAD+ | mM | 0.10 | 0.55 | 1.00 |
B: FSA | mg mL−1 | 1.00 | 2.00 | 3.00 |
C: HLADH | mg mL−1 | 1.00 | 5.50 | 10.00 |
The experimental design was composed of 18 experiments in total (23 + 1 central point in duplicates). The variable ranges were chosen according to the practical assumptions. If efficient regeneration system is present, coenzyme should not be used in the equimolar concentration in comparison to substrate concentration. That is why the maximum concentration was chosen to be only 1 mM. The lower concentration limit was chosen to be 100-fold lower than the concentration of 1. The concentrations of FSA and HLADH were chosen on the basis of preliminary experiments and previous experience which showed that higher increase of their concentration does not have significant influence on the process outcome at the proposed initial concentrations of 1, DHA, and NOX.
The concentration of 3 and 4 formed after 72 h were measured and used as the responses of the process. It was the purpose of the optimization to find the conditions at which maximum yield of 3 will be achieved, and possibly to minimize the over-oxidation of 1 to 4.
A factorial model is composed of a list of coefficients multiplied by an associated factor (investigated variable). It can be presented by the eqn (3), where the coefficients b are associated with a certain variable, designated as letters A, B and C representing the NAD+, FSA and HLADH concentration. Interactions of these variables are presented as the combination of letters.
Response = b0 + b1A + b2B + b3C + b12AB + b13AC + b23BC + b123ABC | (3) |
The objective was to maximize the yield/concentration of 3 and minimize 4 in the investigated variable range.
Experiment | c NAD+/mM | c DHA/mM | γ NOX/mg mL−1 | Y aldol/% | Y amino acid/% |
---|---|---|---|---|---|
a c amino alcohol = 11.24 mM, γHLADH = 5.07 mg mL−1, 0.015 U, γFSA = 5.07 mg mL−1, 11.71 U. This experiment is not presented in Fig. 2. b c amino alcohol = 12.38 mM, γHLADH = 5.00 mg mL−1, 0.015 U, γFSA = 2.09 mg mL−1, 4.83 U. | |||||
1a | 1.00 | 10.02 | 0.12 | 11.73 | 21.89 |
2b | 0.10 | 25.23 | 0.45 | 32.18 | 0.00 |
3b | 1.00 | 10.00 | 0.45 | 31.09 | 23.56 |
4b | 1.00 | 25.23 | 0.18 | 53.89 | 22.27 |
5b | 1.02 | 25.23 | 0.45 | 75.71 | 15.07 |
![]() | ||
Fig. 2 Aldol adduct formed in the cascade reaction of amino alcohol oxidation with coenzyme regeneration and aldol addition reaction (Vr = 1 mL, 50 mM TEA HCl buffer pH 8.0, 10% v/v acetonitrile, camino alcohol = 12.38 mM, γADH = 5.00 mg mL−1, 0.015 U, γFSA = 2.09 mg mL−1, 4.83 U, initial concentrations of other reactants and enzymes are presented in Table 1; ··■·· – experiment 2, ··□·· – experiment 3, ··○·· – experiment 4, ··●·· – experiment 5). |
Theoretical consideration based on our previous studies conducted with FSA A129S variant20,21 and NOX25 were made on process variables that may have crucial influence on the process outcome; i.e. coenzyme, NOX and DHA concentrations. It was assumed that the efficiency of coenzyme regeneration was very important for the oxido-reduction. The experiments presented in Fig. 2 showed that NAD+ and NOX influence the formation of 3, and by increasing their concentration, higher yield of 3 was achieved. Like oxido-reduction, aldol addition is an equilibrium reaction as well and concentration of DHA higher than equimolar was used to shift the equilibrium towards the aldol adduct. Since the reverse C–C bond-forming processes are often favored by thermodynamic relationships, product concentration might be increased by working at higher substrate concentrations20,34 or by increasing the concentration of one of the reactants.34
With the preliminary experiments presented in Table 2 it can be seen that up to 76% of 3 can be formed (Fig. 3) and a low concentration of 2 is actually present in the reaction media during the reaction. Thus, a positive result of increased DHA concentration can be observed. The results also show that oxidation and aldol addition were, in fact, operating simultaneously and the problem of the instability of 2 was effectively resolved. The synthesized 3 was characterized by 1H and 13C NMR analysis. The results are presented in appendix, Fig. 1 together with the HPLC analysis of the product in comparison to the standard (appendix, Fig. 2).
To eliminate the possibility of side reaction between DHA and HLADH and 3 and ADH, spectrophotometric assays were done. It was found that NADH and DHA induce a low activity of HLADH, but in the absence of NADH, the activity is not present. This is very important because DHA and NAD+ do not react in the HLADH catalyzed reaction. With an efficient NAD+ coenzyme regeneration, there is virtually no NADH in the reaction system, and all coenzyme is present in NAD+ form. 3 did not show any activity with HLADH which is also visible from the batch reactor experiments, i.e. from the fact that concentration of 3 remains constant after achieving maximum reaction yield.
Enzyme | k d/h−1 | t 1/2/h |
---|---|---|
HLADH | 0.008 ± 0.003 | 86.6 |
NOX | 0.066 ± 0.015 | 10.5 |
FSA A129S | 0.078 ± 0.004 | 8.9 |
The t1/2 for FSA is in very good agreement with our previous work,20 which was done in different organic co-solvent. The low organic co-solvent concentration used (up to 10% v/v) did not have a significant influence on aldolase activity.21 Furthermore, NOX is more stable in the reaction system than in the absence of substrate which is also in accordance with our previous findings.25
Run | A: cNAD+/mM | B: γFSA/mg mL−1 | C: γHLADH/mg mL−1 | Response 1 caldol adduct/mM | Response 1 camino acid/mM | Y aldol adduct/% | Y amino acid/% |
---|---|---|---|---|---|---|---|
1 | 0.1 | 1 | 1 | 1.14 | 0.06 | 11.4 | 0.6 |
2 | 1 | 3 | 10 | 7.78 | 0.92 | 77.8 | 9.2 |
3 | 0.1 | 3 | 1 | 0.89 | 0.07 | 8.9 | 0.7 |
4 | 1 | 1 | 1 | 6.78 | 0.15 | 67.8 | 1.5 |
5 | 0.1 | 1 | 10 | 3.14 | 0.05 | 31.4 | 0.5 |
6 | 1 | 3 | 1 | 6.45 | 0.07 | 64.5 | 0.7 |
7 | 0.55 | 2 | 5.5 | 6.06 | 0.11 | 60.6 | 1.1 |
8 | 0.1 | 3 | 10 | 2.84 | 0.01 | 28.4 | 0.1 |
9 | 0.1 | 1 | 10 | 3.37 | 0.05 | 33.7 | 0.5 |
10 | 1 | 3 | 10 | 7.91 | 1.06 | 79.1 | 10.6 |
11 | 0.1 | 3 | 1 | 0.81 | 0.06 | 8.1 | 0.6 |
12 | 1 | 1 | 1 | 6.48 | 0.16 | 64.8 | 1.6 |
13 | 0.1 | 1 | 1 | 0.92 | 0.04 | 9.2 | 0.4 |
14 | 1 | 1 | 10 | 7.35 | 2.37 | 73.5 | 23.7 |
15 | 0.1 | 3 | 10 | 2.91 | 0.05 | 29.1 | 0.5 |
16 | 1 | 3 | 1 | 6.26 | 0.11 | 62.6 | 1.1 |
17 | 1 | 1 | 10 | 6.38 | 1.52 | 63.8 | 15.2 |
18 | 0.55 | 2 | 5.5 | 6.28 | 0.13 | 62.8 | 1.3 |
Statistical analysis37 of the results (Table 4) was done and the effect of each variable or their interaction on the objective function was determined. The calculations were done by the software and the results are presented in Table 5. P values show that concentrations of coenzyme (A) and HLADH (C), as well as interactions between NAD+ and FSA (AB), and NAD+ and HLADH (AC) concentrations are significant process variables in the investigated variable range. They have a significant influence on the equilibrium concentration of 3. Table 5 also implies that all single variables, i.e. NAD+ (A), FSA (B) and HLADH (C) concentrations have significant influence on the formation of 4 in the reaction mixture in the investigated variable range. Additionally, interactions between NAD+ and FSA (AB), and NAD+ and HLADH (AC) concentrations also have significant influence on this optimization objective.
Source | Sum of squares | Degrees of freedom | Mean square | F value | P value |
---|---|---|---|---|---|
Aldol adduct concentration | |||||
Model | 7.79134 | 7 | 1.113049 | 339.6448 | <0.0001 |
A | 6.458416 | 1 | 6.458416 | 1970.774 | <0.0001 |
B | 0.001212 | 1 | 0.001212 | 0.369977 | 0.5581 |
C | 0.891378 | 1 | 0.891378 | 272.0024 | <0.0001 |
AB | 0.027157 | 1 | 0.027157 | 8.286876 | 0.0182 |
AC | 0.385004 | 1 | 0.385004 | 117.4832 | <0.0001 |
BC | 0.012195 | 1 | 0.012195 | 3.72132 | 0.0858 |
ABC | 0.015977 | 1 | 0.015977 | 4.875385 | 0.0546 |
Curvature | 0.427986 | 1 | 0.427986 | 130.5991 | <0.0001 |
Pure error | 0.029494 | 9 | 0.003277 | ||
Cor total | 8.24882 | 17 | |||
![]() |
|||||
Amino acid concentration | |||||
Model | 29.29351 | 7 | 4.184787 | 114.5687 | <0.0001 |
A | 16.57768 | 1 | 16.57768 | 453.8544 | <0.0001 |
B | 0.246378 | 1 | 0.246378 | 6.745198 | 0.0289 |
C | 5.134697 | 1 | 5.134697 | 140.5748 | <0.0001 |
AB | 0.472917 | 1 | 0.472917 | 12.94726 | 0.0058 |
AC | 6.837377 | 1 | 6.837377 | 187.1898 | <0.0001 |
BC | 0.023531 | 1 | 0.023531 | 0.644217 | 0.4429 |
ABC | 0.000922 | 1 | 0.000922 | 0.025252 | 0.8772 |
Curvature | 0.068345 | 1 | 0.068345 | 1.871101 | 0.2045 |
Pure error | 0.328738 | 9 | 0.036526 | ||
Cor total | 29.69059 | 17 |
Multiple regression analysis of experimental data resulted in model eqn (4) and (5), which describe the influence of each investigated variable on the objective functions. Their coefficients are given as coded values.
![]() | (4) |
ln(camino acid) = −1.92465 + 1.017893A − 0.12409B + 0.566497C − 0.17192AB + 0.653709AC − 0.03835BC + 0.007593ABC | (5) |
The model adequacy was checked by ANOVA analysis. The model F values for 3 and 4 concentration objectives were 339.64 and 114.57, respectively. This means that both models are significant. There is only a 0.01% chance that F values this large could occur due to noise in both cases. ANOVA results are summarized in Table 5. The model adequacy was tested by the residual analysis (appendix, Fig. 3), residuals vs. predicted plot (appendix, Fig. 4), predicted vs. actual plot (appendix, Fig. 5) etc. R2 values for 3 and 4 concentrations as objective functions were found to be 0.996 and 0.989, respectively.
![]() | ||
Fig. 4 Graphical interpretations of the model (eqn (4)) for aldol adduct concentration. The influence of coenzyme ((a) cNAD+ = 1.00 mM, (b) cNAD+ = 0.10 mM), HLADH ((c) γADH = 10.00 mg mL−1, (d) γADH = 1.00 mg mL−1) and FSA ((e) γFSA = 3.00 mg mL−1, (f) γFSA = 1.00 mg mL−1) concentration on the aldol adduct concentration. |
![]() | ||
Fig. 5 Graphical interpretations of the model (eqn (5)) for amino acid concentration. The influence of coenzyme ((a) cNAD+ = 1.00 mM, (b) cNAD+ = 0.10 mM), HLADH ((c) γADH = 10.00 mg mL−1, (d) γADH = 1.00 mg mL−1) and FSA ((e) γFSA = 3.00 mg mL−1, (f) γFSA = 1.00 mg mL−1) concentration on the amino acid concentration. |
The model eqn (4), as well as its graphical interpretation presented in Fig. 4, suggests that higher coenzyme (A) and HLADH (C) concentration result in higher concentration of 3. The model eqn (5), as well as its graphical interpretation presented in Fig. 5 show that lower content of 4 would be obtained in the reaction mixture if concentrations of coenzyme (A) and HLADH (C) were decreased. Additionally, higher concentrations of FSA favor the formation of lower concentrations of 4 (Fig. 5c and d) which is due to faster consumption of 2.
It can be seen from the simulations (Fig. 4) how concentrations of HLADH, FSA and NAD+ influence the concentration of 3 (Fig. 4a–d). Fig. 5e–f show that concentration of 4 increases, as the concentration of NAD+ increases. This was suggested by our preliminary experiments. Additionally, concentration of 4 will increase at higher concentration of HLADH (Fig. 5e–f) which implies that over-oxidation of 1 to 4 might be catalyzed by HLADH. This was supported by literature38–40 which reports the ability of many ADHs to catalyze the oxidation of amino aldehyde to amino acid (in the presence of NAD(P)+), especially at higher pH values which favor the oxidation reaction. We conducted additional experiments with HLADH from a different supplier (prepared and purified by prof. Martina Pohl's research group, Research Center Jülich), and it was found that 4 also forms, which further confirms this hypothesis. Considering that this statistical analysis showed that both HLADH and NAD+ concentrations should be increased to produce more 4, it is an additional corroboration to this hypothesis.
Tables 1 and 2 in appendix present the optimal solutions for the separate objective functions calculated by the software. Since similar conditions regarding the concentrations of coenzyme and HLADH favor both the formation of 3 and 4, it was not desirable to join these two objectives into one optimization. This optimization would have been a compromise and would result in lower concentration of 3. Additionally, maximum yield of 3 was 79% (Table 4, exp 10) which is a good result. Hence, the results obtained in this work represent a remarkable achievement in enzymatic cascade oxidation–aldol reactions. For comparison, in a one-pot two-step cascade reaction system with oxidation-C–C bond formation using an analogous N-Cbz-amino alcohol, Mifsud and coauthors23 obtained 19% of aldol adduct and 81% yield of N-Cbz-amino acid in aqueous buffered media.
Optimization of initial conditions of cascade reaction showed that higher concentrations of coenzyme and HLADH favor both the formation of 3 and 4. Maximal yield of 3 achieved in this cascade reaction was 79%. At the same time yield of 4 was 10%. The system might be further improved by finding the appropriate ADH which would not catalyze the over-oxidation of amino aldehyde to amino acid.
c | Molar concentration, mM |
k d | Operational stability decay rate constant, h−1 |
S.A. | Specific activity, U mg−1 |
T | Reaction time, min |
t 1/2 | Half-life time, h |
V.A. | Volume activity, U mL−1 |
V enz | Enzyme volume, mL |
V r | Reactor volume, mL |
Y | Reaction yield, % |
γ | Mass concentration, mg mL−1 |
HLADH | Horse liver alcohol dehydrogenase |
DHA | Dihydroxyacetone |
FSA | D-Fructose-6-phosphate aldolase |
HCl | Hydrochloric acid |
NOX | NADH oxidase |
TEA | Triethanolamine |
TFA | Trifluoroacetic acid |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c5ra14414k |
This journal is © The Royal Society of Chemistry 2015 |