Team:Edinburgh/Cellulases (MATLAB model)

Cellulases (MATLAB model)

The most important part of the biorefinery is the reactor where cellulose is converted to glucose. But accurately predicting how much is converted, using synergy between enzymes is difficult without literature to provide the ordinary differential equations (ODE's) and the kinetic parameters. Therefore this model only looks at the free floating enzyme approach (non-synergy). It is deterministic i.e non random and is set by a series of initial conditions.

Assumptions

The mathematical model is based on the ODE's and kinetic parameters outlined in [http://onlinelibrary.wiley.com/doi/10.1021/bp034316x/full Kadam et al, 2004]. The following are its assumptions and basis:

• Rate equations assume enzyme adsorption follows the Langmuir isotherm model
• Glucose and cellobiose which are the products of cellulose hydrolysis, were assumed to, 'competitively inhibit enzyme hyrolysis' [http://onlinelibrary.wiley.com/doi/10.1021/bp034316x/full Kadam et al, 2004]
• Assume all reactions follow the same temperature dependency Arrhenius relationship. However it should be different for every enzyme component, 'because of their varying degrees of thermostability, with beta-glucocidase being the most thermostable. Hence the assumption is a simplification of reality' [http://onlinelibrary.wiley.com/doi/10.1021/bp034316x/full Kadam et al, 2004]
• Conversion of cellobiose to glucose follows the Michaelis-Menten enzyme kinetic model

Equations

Rate Equations

Cellulose to Cellobiose Reaction with Competitive Glucose, Cellobiose and Xylose Inhibition.

Cellulose to Glucose Reaction with Competitive Glucose, Cellobiose and Xylose Inhibition.

Cellobiose to Glucose Reaction with competitive Glucose, Cellibiose and Xylose Inhibition.

Constants

knr - reaction rate constant for reaction n EnB is the bound concentration for exo and endo-beta-1,4-glucanase for reaction n
Rs - substrate reactivity parameter
S - substrate reactivity at a given time
G2 - concentration of cellobiose
G - concentration of Glucose
X - xylose concentration
KnIG2 - inhibition constant for cellobiose at reaction n
K1IG - inhibition constant for Glucose at reaction n
K1IX - xylose inhibition constant for reaction n
Note: Assuming no xylose inhibition therefore X=0

References

• Kadam KL, Rydholm EC, McMillan JD (2004) [http://onlinelibrary.wiley.com/doi/10.1021/bp034316x/full Development and Validation of a Kinetic Model for Enzymatic Saccharification of Lignocellulosic Biomass]. Biotechnology Progress 20(3): 698–705 (doi: 10.1021/bp034316x).