Revision as of 13:00, 17 October 2018

Mathematical Modeling

Introduction

The team decided to use a scaffold to bind together the cellulases that will allow the yeast to degrade cellulose. Experimental data suggests that cellulose gets degraded more effectively this way. Intuitively the scaffold allows the yeast to more efficiently absorb the products of the degradation due to to proximity, but does scaffolding come at a cost? This model aims to simulate a scaffold similar to the one used in the lab mathematically, to see whether it suffers from any disadvantage and whether the choice of cellulases make a difference. For comparison, a system of free enzymes was also modelled in a similar way.

Method

It is assumed that any advantage of scaffolding will also occur for a system of three cellulases. The model follows the degradation of cellulose chains by two enzymes, an exoglucanase and endoglucanase. To simplify the problem further, it is assumed that the cellulase chains are maximally 100 units long and are all completely soluble. This shouldn’t stop the model from working because cellulose is commercially available in soluble form and at short chain length distributions.

Free Enzyme Model

Symbol	Meaning	Unit
C_i	Concentration Of Cellulose chain of length i.	mmol/L
E_CBHI, E_EG2	Concentration Of Enzymes Cellobiohydrolase and Endoglucanase	mmol/L

Figure 1. Table Of Variables For Mechanistic Model

The breakdown of cellulose by two free cellulases is described by balances in figure 2.

*Figure 2. Reaction Balances For Exoglucanases and Endoglucanase*

As described by Levine et al., the complexation step of the cellulase enzyme to the cellulose is so fast compared to the catalysis, that the complex concentration is considered to be in equilibrium.
This allows the for expression of the complexes in terms of the concentrations of the chains, enzymes and equilibrium constants (also called dissociation constants), figure 3.

*Figure 3. Equations for the Concentrations of Complexes from equilibrium assumptions*

It must be kept in mind that the enzymes are complexing with available sites. An endoglucanase binds to the 1-4 glucose bonds, so the concentration available is the concentration of chain i, times the amount of binding sites on chain i. The exoglucanase can only attach to chain ends.
The equilibrium constants are assumed to be the same regardless of the length of the chain. The values for the equilibrium constants were retrieved from the Levine et al. simulation.

Next the concentration balances for each cellulose were derived, to create a system of differential equations. The method for deriving these equations is shown in figure 4 with the chain-length between 3-98 as example.

*Figure 4. Derivation Of the rate balance equation for lengths i from 3 to 99 ( Different equations apply for 1,2,99 and 100)*

The other rate equations, not shown, are derived in exactly the same way. The kcat values for the catalytic steps were also retrieved from Levine’s paper. A summary of all the constants used can be seen in figure 5.

Symbol	Value	Unit
K_eq,EG2	0.053	mmol L^-1
K_eq,CBH1	0.0032	mmol ^L-1
k_EG2,cat	65	s^-1
k_CBH1,cat	9.4	s^-1

Figure 5. Table Of Constants used for free enzyme model.

The concentrations of free enzymes (not in complex) is also a dependent variable, and therefore must be described by an equation for the system to be solvable. This can be done with the general Derivation in figure 6.

*Figure 6. Derivation of concentration of non-complexed enzymes*

Since the complexation is considered instantaneous, initial conditions for the free enzyme are calculated by subtracting the enzymes complexed with the cellulose feed from the enzyme feed.

Scaffolded Enzyme Model

To model the effect of a scaffold, several additional assumptions are made:

The scaffold will benefit from the highest affinity of the constituent enzymes.
The catalytic rate of the scaffold is determined by the ratio of the catalytic rates of the individual enzymes. The smallest one having the highest contribution, to be rate determining.
Upon complexation, The scaffold produces the products of both enzymes in one step. (see figure 7)

*Figure 7. Reaction Balance For the Scaffolded Enzymes*

The cellulose differentials for the scaffold are derived in the same way as the free enzymes. This time, only one enzyme acts on the cellulose, and every reaction leads to the production of cellobiose.

It can be seen the affinity of the Endoglucanase for a particular chain goes up with its length, as evident from the factor (J-1). At times, the affinity of the Endoglucanase is higher than that of the Exoglucanase. According to the assumptions, the scaffold should benefit from the highest affinity ( lowest dissociation constant). To do this, the dissociation constant is allowed vary with the length of the corresponding chain by the equation in figure 8. This equation is a smooth approximation of the maximum of the affinities of both enzymes as a function of chain length. It had to be used because the mathematica solver cannot work with discontinuities.

*Figure 8. Approximation for the equilibrium constant of the scaffold as a function of Chain length (J)*

Figure 9 shows how the dissociation constants of exoglucanase and endoglucanases vary with chain length along with function this model uses to approximate that of the scaffold’s..

Figure 9. Graph showing Effective Dissociation constant of CBH1 (Green), EG2(Yellow), Approximation For Scaffold (Blue). The small interval shows where the approximation becomes very slightly limited at chain length 17, the only place were the approximation can be seen. X-axis = Chain Length(glucose units) , Y-Axis = Affinity(L/mmol)

Results

The results show a sigmoidal production of short-length sugars with time. It can be seen that the scaffold is slightly slower than the system of free enzymes. The degradation of the C100 chains show an exponential decreasing relationship with time. Here, also it can be seen that the free enzymes of are slightly faster than the scaffold. The effect of the difference in kCats of the constituent enzymes was also plotted (fig 12). Decreasing the difference in kCats causes the production of absorbable products to increase in rate, a rightward shift in the function.

*Figure 10. Production Of the soluble sugars By Scaffold (Yelllow), By Free enzyme (Blue). X-axis = Time(hours), Y-Axis= Fraction of glucose monomers in product form*

On figure 12 the effect of varying the difference between the constituent enzyme’s kcats can be seen. As the kCats get closer together, the efficiency goes up. Although for cellulases this effect isn’t very strong, this would be something important to keep in mind in future scaffolding experiments.

Figure 12. Effect of increasing difference in scaffold kcats. As the the difference increases, the graph shifts to the right, efficiency increase. X-axis = Time(hours), Y-Axis= Fraction of glucose monomers in product form

Conclusions

In this model a scaffold was modelled using a system of differential equations. The performance of the scaffold on degrading cellulose chains was compared with that of free enzymes. The effect of varying the kCats between the constituent enzymes was also investigated. The results show that the scaffold performs overall a bit slower than the free enzymes. However, the difference isn't big enough to say that using a scaffold would cause a significant disadvantage in terms of reaction rate. The team could therefore use it and take advantage of the other advantages it offers, such as giving the yeast cell higher proximity to the products. The effect of difference in kCat between the constituent enzymes of the scaffold was also investigated. The results show that decreasing the difference improves the overall performance. This finding does not isignificantly nfluence the choice of enzymes for this particular scaffold but may be useful in future scaffold design.

Main Sources

-Levine, S. E., Fox, J. M., Blanch, H. W. and Clark, D. S. (2010), A mechanistic model of the enzymatic hydrolysis of cellulose. Biotechnol. Bioeng., 107: 37-51. doi:10.1002/bit.22789

-Wen F, Sun J, Zhao H. Yeast Surface Display of Trifunctional Minicellulosomes for Simultaneous Saccharification and Fermentation of Cellulose to Ethanol . Applied and Environmental Microbiology. 2010;76(4):1251-1260. doi:10.1128/AEM.01687-09.

@@ Line 43: / Line 43: @@
       <p>It can be seen the affinity of the Endoglucanase for a particular chain goes up with its length, as evident from the factor (J-1). At times, the affinity of the Endoglucanase is higher than that of the Exoglucanase. According to the assumptions, the scaffold should benefit from the highest affinity ( lowest dissociation constant). To do this, the dissociation constant is allowed vary with the length of the corresponding chain by the equation in figure 8. This equation is a smooth approximation of the maximum of the affinities of both enzymes as a function of chain length. It had to be used because the mathematica solver cannot work with discontinuities.</p>
       <figure><img src="https://static.igem.org/mediawiki/2018/e/ec/T--Groningen--Stochastic-8.png" width="70%"><figcaption><i>Figure 8. Approximation for the equilibrium constant of the scaffold as a function of Chain length (J)</i></figcaption></figure> <p>Figure 9 shows how the dissociation constants of exoglucanase and endoglucanases vary with chain length along with function this model uses to approximate that of the scaffold’s..</p>
-      <figure><img src="https://static.igem.org/mediawiki/2018/e/e8/T--Groningen--Stochastic-3-2.png9.jpg" width="100%"><figcaption><i>Figure 9. Graph showing Effective Dissociation constant of CBH1 (Green), EG2(Yellow), Approximation For Scaffold (Blue). The small interval shows where the approximation becomes very slightly limited at chain length 17, the only place were the approximation can be seen.  X-axis = Time(Hours) , Y-Axis = Affinity(L/mmol)</i></figcaption></figure>
+      <figure><img src="https://static.igem.org/mediawiki/2018/e/e8/T--Groningen--Stochastic-3-2.png9.jpg" width="100%"><figcaption><i>Figure 9. Graph showing Effective Dissociation constant of CBH1 (Green), EG2(Yellow), Approximation For Scaffold (Blue). The small interval shows where the approximation becomes very slightly limited at chain length 17, the only place were the approximation can be seen.  X-axis = Chain Length(glucose units) , Y-Axis = Affinity(L/mmol)</i></figcaption></figure>
       <h4>Results</h4>
       <p>The results show a sigmoidal production of short-length sugars with time. It can be seen that the scaffold is slightly slower than the system of free enzymes. The degradation of the C100 chains show an exponential decreasing relationship with time. Here, also it can be seen that the free enzymes of are slightly faster than the scaffold. The effect of the difference in kCats of the constituent enzymes was also plotted (fig 12).  Decreasing the difference in kCats causes the production of absorbable products to increase in rate, a rightward shift in the function. </p><figure><img src="https://static.igem.org/mediawiki/2018/9/93/T--Groningen--Stochastic-10.jpg" width="100%"><figcaption><i>Figure 10. Production Of the soluble sugars By Scaffold (Yelllow), By Free enzyme (Blue). X-axis = Time(hours), Y-Axis= Fraction of glucose monomers in product form</i></figcaption>

Difference between revisions of "Template:Groningen/Stochastic Modeling"