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.

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.

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

Mechanism-based kinetic models are rigorous tools to analyze enzymatic reactions, but their extension to actual conditions of the biocatalytic synthesis can be difficult. Here, we demonstrate (mechanistic-empirical) hybrid modeling for systematic optimization of the sucrose phosphorylase-catalyzed glycosylation of glycerol from sucrose, to synthesize the cosmetic ingredient α-glucosyl glycerol (GG). The empirical model part was developed to capture nonspecific effects of high sucrose concentrations (up to 1.5 M) on microscopic steps of the enzymatic trans-glycosylation mechanism. Based on verified predictions of the enzyme performance under initial rate conditions (Level 1), the hybrid model was expanded by microscopic terms of the reverse reaction to account for the full-time course of GG synthesis (Level 2). Lastly (Level 3), the application of the hybrid model for comprehensive window-of-operation analysis and constrained optimization of the GG production (~250 g/L) was demonstrated. Using two candidate sucrose phosphorylases (from Leuconostoc mesenteroides and Bifidobacterium adolescentis), we reveal the hybrid model as a powerful tool of "process decision making" to guide rational selection of the best-suited enzyme catalyst. Our study exemplifies a closing of the gap between enzyme kinetic models considered for mechanistic research and applicable in technologically relevant reaction conditions; and it highlights the important benefit thus realizable for biocatalytic process development.

Enzyme-Mediated Lignocellulose Liquefaction Is Highly Substrate-Specific and Influenced by the Substrate Concentration or Rheological Regime.

The high viscosities/yield stresses of lignocellulose slurries makes their industrial processing a significant challenge. However, little is known regarding the degree to which liquefaction and its enzymatic requirements are specific to a substrate's physicochemical and rheological properties. In the work reported here, the substrate- and rheological regime-specificities of liquefaction of various substrates were assessed using real-time in-rheometer viscometry and offline oscillatory rheometry when hydrolyzed by combinations of cellobiohydrolase (Trichoderma reesei Cel7A), endoglucanase (Humicola insolens Cel45A), glycoside hydrolase (GH) family 10 xylanase, and GH family 11 xylanase. In contrast to previous work that has suggested that endoglucanase activity dominates enzymatic liquefaction, all of the enzymes were shown to have at least some liquefaction capacity depending on the substrate and reaction conditions. The contribution of individual enzymes was found to be influenced by the rheological regime; in the concentrated regime, the cellobiohydrolase outperformed the endoglucanase, achieving 2.4-fold higher yield stress reduction over the same timeframe, whereas the endoglucanase performed best in the semi-dilute regime. It was apparent that the significant differences in rheology and liquefaction mechanisms made it difficult to predict the liquefaction capacity of an enzyme or enzyme cocktail at different substrate concentrations.

On the donor substrate dependence of group-transfer reactions by hydrolytic enzymes: Insight from kinetic analysis of sucrose phosphorylase-catalyzed transglycosylation.

Chemical group-transfer reactions by hydrolytic enzymes have considerable importance in biocatalytic synthesis and are exploited broadly in commercial-scale chemical production. Mechanistically, these reactions have in common the involvement of a covalent enzyme intermediate which is formed upon enzyme reaction with the donor substrate and is subsequently intercepted by a suitable acceptor. Here, we studied the glycosylation of glycerol from sucrose by sucrose phosphorylase (SucP) to clarify a peculiar, yet generally important characteristic of this reaction: partitioning between glycosylation of glycerol and hydrolysis depends on the type and the concentration of the donor substrate used (here: sucrose, α-d-glucose 1-phosphate (G1P)). We develop a kinetic framework to analyze the effect and provide evidence that, when G1P is used as donor substrate, hydrolysis occurs not only from the ß-glucosyl-enzyme intermediate (E-Glc), but additionally from a noncovalent complex of E-Glc and substrate which unlike E-Glc is unreactive to glycerol. Depending on the relative rates of hydrolysis of free and substrate-bound E-Glc, inhibition (Leuconostoc mesenteroides SucP) or apparent activation (Bifidobacterium adolescentis SucP) is observed at high donor substrate concentration. At a G1P concentration that excludes the substrate-bound E-Glc, the transfer/hydrolysis ratio changes to a value consistent with reaction exclusively through E-Glc, independent of the donor substrate used. Collectively, these results give explanation for a kinetic behavior of SucP not previously accounted for, provide essential basis for design and optimization of the synthetic reaction, and establish a theoretical framework for the analysis of kinetically analogous group-transfer reactions by hydrolytic enzymes.