From 2011.igem.org

Latest revision as of 01:59, 29 October 2011

Modeling of cell-free bisphenol A degradation

In our project we want to construct a cell-free biosensor that can detect bisphenol A (BPA). To be more specific we want to fuse proteins to silica beads that degrade the BPA and detect NAD⁺ that is set free in the reaction by a NAD⁺-dependent DNA ligase (also attached to the beads) in combination with a molecular beacon.

Modeling cell-free processes

The fact that our project aims at the construction of a cell-free biosensor has huge advantages regarding the complexity of the model, as it is possible to simulate our project in a small number of differential equations. There are several reasons why this is possible. First off, we do not need to consider the transcription- and translationrate and other factors influencing the molecular amount of proteins in cells. After extracting our proteins from the E.coli cells we know the exact amount of protein introduced to our cell-free system. Furthermore we can state that the protein amount in our system is constant, as our system does not contain any proteases. Knowing the exact amount of enzymes in our system supplies the modeling with another advantage. We only need to simulate the reactions that are rate-limiting. In case of our model we state that only the BPA concentration is limiting while the enzymes and the NADH are present in surplus quantities. This allows to simplify the first reaction a great deal, as it is not necessary to simulate the electron transport in the Cytochrome P450 monooxygenase system. The third, and by far most important advantage of cell-free modeling is, that we do not need to make a lot of assumptions regarding the system we want to simulate the reaction in. When modeling a reaction in a cell there are a lot of parameters difficult to estimate, for example the contribution of the reactants in the cell or what signaling networks might be triggered by the reaction of interest.

First reaction

The first reaction is the degradation of BPA by the Cytochrome P450 monooxygenase system (the components of this system are cytochrome P450 (<partinfo>K123000</partinfo>), ferredoxin (<partinfo>K123001</partinfo>) and ferredoxin reductase (<partinfo>K525499</partinfo>)). The velocity of the BPA degradation can be calculated by Michaelis Menten enzyme kinetics:

The velocity V₁ was calculated dependent on the enzyme concentration and the k_cat-value of the monooxygenase system which was 3.9 min^-1([http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1317416/ Sasaki, 2005b]).

Second reaction

In the second reaction the NAD⁺-dependent DNA-ligase from E.coli (<partinfo>K525710</partinfo>) is used to transfer the molecular beacon from its closed state (b_cs) to its open state b_os. It can be described by the following reaction scheme:

The rate of a reaction with two substrates is calculated by the following equation (assuming initial product concentration is 0):

where K_iA is the binding affinity of the inhibitor to substrate A, K_mA and K_mB the Michaelis-Menten constants for substrate A respectively B in the specific reaction. So the amount of the molecular beacons open state is described by the equation

For this we assume that the molecular amount of NAD⁺ is equal to the molecular amount of degraded BPA less the amount of beacon in open state. The required electrons for the degradation of one BPA molecule, more precisely for the activation of the oxygen in the cytochrome P450 monooxygenase system, are released by the oxidation of one molecule of the cofactor NADH. The transformation reaction of beacon in closed state to beacon in open state decreases the NAD⁺ amount by one molecule.

The velocity V₂ is calculated by multiplying the k_cat-value of 0.0212 s^-1 found in the [http://onlinelibrary.wiley.com/doi/10.1046/j.1432-1327.2000.01377.x/pdf literature] with the used concentration of NAD⁺. For K_mNAD⁺ and K_mbos we use the values we found in the literature ([http://www.sciencedirect.com/science/article/pii/S0003269702003020 Chen et al., 2002]):

Regarding K_iNAD⁺ we considered AMP and NMN as possible inhibitors ([http://www.sciencedirect.com/science/article/pii/S0003269702003020 Chen et al., 2002]) for the reaction. The literature provides half maximal inhibitory concentrations (IC₅₀-values) of 3.1 mM for AMP and 9.5 µM for NMN. We assume that AMP has no significant influence in our model, as the IC₅₀-value is quite high. The IC₅₀-value of NMN on the other hand is in closer range to concentrations that might be achieved during the performance of our test, so we do consider it in our calculations. To determine K_iNAD⁺ we used the formula

yielding an value of 0.45 µM for NMN.

Implementation of the equations for the first and the second reaction in matlab:

Figure 2 shows the simulation of our model for a timespan of 120 minutes.

Fig. 2: Simulation of extracellular degradation of BPA in combination with the activation of the molecular beacon.

Modeling of intracellular bisphenol A degradation

Regular model

The modeling was done with the software [http://www.berkeleymadonna.com/ Berkeley Madonna] using the [http://en.wikipedia.org/wiki/Runge–Kutta_methods#Common_fourth-order_Runge.E2.80.93Kutta_method common fourth-order Runge-Kutta] method to solve the equations. The model was fitted to the measured data shown above by the function "curve fit" in Berkeley Madonna to calculate the parameters, constants etc.

To model the BPA degradation by E. coli carrying BioBricks for BPA degradation (<partinfo>K123000</partinfo> and <partinfo>K123001</partinfo>) the cell growth has to be described first to calculate a specific BPA degradation rate per cell. The observed growth of E. coli on (our) LB medium was [http://en.wikipedia.org/wiki/Diauxie diauxic] with two different growth phases. Cell growth is a [http://en.wikipedia.org/wiki/First_order_kinetics#First-order_reactions first-order reaction] and is mathematically described as

with the specific growth rate µ and the cell count X. The specific growth rate is dependent on the concentration of the growth limiting substrate (e.g. glucose) and can be described as

with the substrate concentration S, the Monod constant K_S and the maximal specific growth rate µ_max ([http://www.annualreviews.org/doi/abs/10.1146/annurev.mi.03.100149.002103 Monod, 1949]). Because LB medium is a complex medium we cannot measure the substrate concentration so we have to assume an imaginary substrate concentration. Due to the diauxic growth two different substrates with different Monod constants and consumption rates are necessary to model the cell growth. The amount of a substrate S can be modeled as follows

with the specific substrate consumption rate per cell q_S. The whole model for the diauxic growth of E. coli on LB medium with two not measurable (imaginary) substrates looks like:

The specific BPA degradation rate per cell q_D is modeled with an equation like eq. (7). In the beginning of the cultivations, when E. coli growths on the "good" imaginary substrate S₁, no BPA degradation is observed. When this substrate is consumed, the BPA degradation starts. The model for this diauxic behavior is as follows:

Figures 3 - 6 show a comparison between modeled and measured data for cultivations with BPA degrading E. coli. In Tab. 1 the parameters for the model are given, obtained by curve fitting the model to the data.

Model for the fusion protein between FNR, BisdA and BisdB

The cultivations and BPA degradation of E. coli KRX carrying a fusion protein consisting of ferredoxin-NADP⁺ oxidoreductase, ferredoxin_bisd and cytochrome P450_bisd differ from the cultivations with the other BPA degrading BioBricks. Thus, the model for these cultivations has to be adjusted to this behaviour. First of all, no diauxic growth is observed so the growth can be modeled more easily like

The BPA degradation starts when the imaginary substrate is depleted, like observed in the other cultivations with BPA degrading BioBricks. But it seems that the BPA degradation is getting slower with longer cultivation time. So it is modeled with a Monod-like term in which the specific BPA degradation rate is dependent from the BPA concentration:

with the maximal specific BPA degradation rate q_D,max and the constant K_D.

Figures 7 shows a comparison between modeled and measured data for cultivations with BPA degrading fusion protein E. coli. In Tab. 1 the parameters for the model are given, obtained by curve fitting the model to the data.

Comparison between modeled and measured data

Figure 3: Comparison between modeled (lines) and measured (dots) data for cultivations of E. coli KRX carrying BPA degrading BioBrick <partinfo>K525512</partinfo>. The BioBrick <partinfo>K525512</partinfo> (polycistronic bisdAB genes behind medium strong promoter) was cultivated seven times in E. coli KRX in LB + Amp + BPA medium at 30 °C, using 300 mL shaking flasks without baffles with silicon plugs. The BPA concentration (closed dots) and the cell density (open dots) is plotted against the cultivation time.

Figure 4: Comparison between modeled (lines) and measured (dots) data for cultivations of E. coli KRX carrying BPA degrading BioBrick <partinfo>K525517</partinfo>. The BioBrick <partinfo>K525517</partinfo> (fusion protein between BisdA and BisdB behind medium strong promoter) was cultivated five times in E. coli KRX in LB + Amp + BPA medium at 30 °C, using 300 mL shaking flasks without baffles with silicon plugs. The BPA concentration (closed dots) and the cell density (open dots) is plotted against the cultivation time.

Figure 5: Comparison between modeled (lines) and measured (dots) data for cultivations of E. coli KRX carrying BPA degrading BioBrick <partinfo>K525552</partinfo>. The BioBrick <partinfo>K525552</partinfo> (polycistronic fnr : bisdA : bisdB genes behind medium strong promoter) was cultivated four times in E. coli KRX in LB + Amp + BPA medium at 30 °C, using 300 mL shaking flasks without baffles with silicon plugs. The BPA concentration (closed dots) and the cell density (open dots) is plotted against the cultivation time.

Figure 6: Comparison between modeled (lines) and measured (dots) data for cultivations of E. coli KRX carrying BPA degrading BioBrick <partinfo>K525582</partinfo>. The BioBrick <partinfo>K525582</partinfo> (polycistronic fnr : BisdAB (fusion protein) genes behind medium strong promoter) was cultivated four times in E. coli KRX in LB + Amp + BPA medium at 30 °C, using 300 mL shaking flasks without baffles with silicon plugs. The BPA concentration (closed dots) and the cell density (open dots) is plotted against the cultivation time.

Figure 7: Comparison between modeled (lines) and measured (dots) data for cultivations of E. coli KRX carrying BPA degrading BioBrick <partinfo>K525562</partinfo>. The BioBrick <partinfo>K525562</partinfo> (fusion protein between FNR, BisdA and BisdB behind medium strong promoter) was cultivated four times in E. coli KRX in LB + Amp + BPA medium at 30 °C, using 300 mL shaking flasks without baffles with silicon plugs. The BPA concentration (closed dots) and the cell density (open dots) is plotted against the cultivation time.

The specific BPA degradation rate per cell q_D is about 50 % higher when using the fusion protein compared to the polycistronic bisdAB gene. This results in an average 9 hours faster, complete BPA degradation by E. coli carrying <partinfo>K525517</partinfo> compared to <partinfo>K525512</partinfo> as observed during our cultivations. The fusion protein between BisdA and BisdB improves the BPA degradation by E. coli. Introducing a polycistronic ferredoxin-NADP⁺ reductase gene into these systems does not lead to a higher specific BPA degradation rate. The rates are a bit lower, though. The effect that the BisdA | BisdB fusion protein is more efficient than BisdB when expressed polycistronically with BisdA can still be observed in this setup. The cultivations with the expression of the fusion protein FNR | BisdA | BisdB differ from the expression of the other BPA degrading BioBricks. The BPA degradation rate is concentration dependent which is typical for enzymatic reactions but was not observed in the other cultivations. In addition, the growth was faster and not diauxic. The maximal specific BPA degradation rate of this BioBrick is higher than the observed specific BPA degradation rates in the other cultivations. But due to the dependence of the BPA degradation from the BPA concentration the BPA degradation with this BioBrick is not as efficient and not complete. Anyway, the fact that this BioBrick is working is impressive.

References

Sasaki M, Akahira A, Oshiman K, Tsuchido T, Matsumura Y (2005b) Purification of Cytochrome P450 and Ferredoxin, Involved in Bisphenol A Degradation, from Sphingomonas sp. Strain AO1, [http://aem.asm.org/cgi/content/abstract/71/12/8024 Appl Environ Microbiol 71(12):8024-8030].

Georlette D, Jonsson Z, Van Petegem F, Chessa J-P, Van Beeumen J, Huebscher U, Gerday C (2000) A DNA ligase from psychrophile Pseudoalteromonas haloplanktis gives insigth into the adaption of proteins to low temperatures, [http://onlinelibrary.wiley.com/doi/10.1046/j.1432-1327.2000.01377.x/pdf Eur.J.Biochem 267:3502-3512].

Chen X, Hentz N, Hubbard F, Meier T, Sittampalam S, Zhao G (2002) Development of a fluorescence resonance energy transfer assay for measuring the activity of Streptococcus pneumoniae DNA ligase, an enzyme essential for DNA replication, repair and recombination, [http://www.sciencedirect.com/science/article/pii/S0003269702003020 Anal Biochem 309:232-240].

Monod J (1949) The growth of bacterial cultures, Annu Rev Microbiol [http://www.annualreviews.org/doi/abs/10.1146/annurev.mi.03.100149.002103 3:371-394].