Three-level hybrid modeling for systematic optimization of biocatalytic synthesis: α-glucosyl glycerol production by enzymatic trans-glycosylation from sucrose

Pushing the boundaries of phosphorylase cascade reaction for cellobiose production II: Model-based multiobjective optimization

The inherent complexity of coupled biocatalytic reactions presents a major challenge for process development with one-pot multienzyme cascade transformations. Kinetic models are powerful engineering tools to guide the optimization of cascade reactions towards a performance suitable for scale up to an actual production. Here, we report kinetic model-based window of operation analysis for cellobiose production (≥100 g/L) from sucrose and glucose by indirect transglycosylation via glucose 1-phosphate as intermediate. The two-step cascade transformation is catalyzed by sucrose and cellobiose phosphorylase in the presence of substoichiometric amounts of phosphate (≤27 mol% of substrate). Kinetic modeling was instrumental to uncover the hidden effect of bulk microviscosity due to high sugar concentrations on decreasing the rate of cellobiose phosphorylase specifically. The mechanistic-empirical hybrid model thus developed gives a comprehensive description of the cascade reaction at industrially relevant substrate conditions. Model simulations serve to unravel opposed relationships between efficient utilization of the enzymes and maximized concentration (or yield) of the product within a given process time, in dependence of the initial concentrations of substrate and phosphate used. Optimum balance of these competing key metrics of process performance is suggested from the model-calculated window of operation and is verified experimentally. The evidence shown highlights the important use of kinetic modeling for the characterization and optimization of cascade reactions in ways that appear to be inaccessible to purely data-driven approaches.

Pushing the boundaries of phosphorylase cascade reaction for cellobiose production I: Kinetic model development

One-pot cascade reactions of coupled disaccharide phosphorylases enable an efficient transglycosylation via intermediary α-d-glucose 1-phosphate (G1P). Such transformations have promising applications in the production of carbohydrate commodities, including the disaccharide cellobiose for food and feed use. Several studies have shown sucrose and cellobiose phosphorylase for cellobiose synthesis from sucrose, but the boundaries on transformation efficiency that result from kinetic and thermodynamic characteristics of the individual enzyme reactions are not known. Here, we assessed in a step-by-step systematic fashion the practical requirements of a kinetic model to describe cellobiose production at industrially relevant substrate concentrations of up to 600 mM sucrose and glucose each. Mechanistic initial-rate models of the two-substrate reactions of sucrose phosphorylase (sucrose + phosphate → G1P + fructose) and cellobiose phosphorylase (G1P + glucose → cellobiose + phosphate) were needed and additionally required expansion by terms of glucose inhibition, in particular a distinctive two-site glucose substrate inhibition of the cellobiose phosphorylase (from Cellulumonas uda). Combined with mass action terms accounting for the approach to equilibrium, the kinetic model gave an excellent fit and a robust prediction of the full reaction time courses for a wide range of enzyme activities as well as substrate concentrations, including the variable substoichiometric concentration of phosphate. The model thus provides the essential engineering tool to disentangle the highly interrelated factors of conversion efficiency in the coupled enzyme reaction; and it establishes the necessary basis of window of operation calculations for targeted optimizations toward different process tasks.