Model
In this page you will find our constructed mathematical model with its corresponding parameters.
![](https://static.igem.wiki/teams/4460/wiki/img/illustrations/rum-igem/component-vector.png)
![](https://static.igem.wiki/teams/4460/wiki/img/illustrations/rum-igem/component-vector-2.png)
Introduction
A mathematical model is known as a tool used to describe the relationship between two or more variables. When referring to modeling, it can be interpreted as simulating reality to formulate ideas precisely and prevent implicit assumptions (Bender, 2000). The bioreactor element was designed according to specific physical parameters and conceptual framework based on the Biodegradation device of the biological circuit. Developing a model contributes to understanding the bioreactor internal process of remediating RDX, by describing behavior and efficiency of substrate consumption into products.
![Anones Lagoon Data as seen from Google Earth](https://static.igem.wiki/teams/4460/wiki/anones-lagoon-contaminated.png)
Figure 1. Current view of contaminated Anones Lagoon from Google Earth
Our Models - Enzyme and Molecule Quantification
For the construction of a mathematical model that could provide insight into the future implementation of the project R-DETOX, it was decided to construct an equation that would result in the time, in months, of detoxification that could be achieved by the Device 2 of our modified genetic circuit. As a software, Microsoft Excel was used since, according to Wu, Ako and Hu in “A Useful Microsoft Excel Add-in Program for Modeling Steady-state Enzyme Kinetics” (2011), Excel is a fast and reliable tool for enzyme kinematic data analysis. It is our project’s purpose to mitigate the contamination in the Anones Lagoon produced by biochemical remains and with this model we could see a clear picture of the quantity of enzyme needed to reduce enough RDX so that the lagoon’s water falls below toxicity levels established by the EPA.
Using the applications Google Earth and Google Maps, which provide detailed information and views about geographical regions worldwide, we could approximate length and width measurements in order to then calculate the estimated volume for the lagoon. We searched for the Anones Lagoon and found measurements for length (848.17 m) and width (313.31 m). For the depth measurement, we were not able to find trustworthy information through research. However, we asked various Vieques residents that knew the area and have previously been involved in investigations of it if they could give us an approximation of the lagoon’s depth. With their answers, it was concluded on average that the lagoon’s depth varied in intervals of 3 through 16 feet (1 through 5 meters) depending on craters in the lagoon’s soil that have been developing over the years. Given this depth approximation, we were able to calculate the lagoon’s volume according to five different depth values: 1 m, 2 m, 3 m, 4 m, and 5 m. Table 1 shows the data approximations for the volume calculation of the Anones Lagoon and Table 2 shows the values for volume calculations in cubic meters and liters according to the five different depth measurements selected.
![Table 1: Anones Lagoon Data](https://static.igem.wiki/teams/4460/wiki/table-1-model.png)
![Table 2: Volume Calculations per depth](https://static.igem.wiki/teams/4460/wiki/table-2-model.png)
![Table 2: Volume Calculations per depth](https://static.igem.wiki/teams/4460/wiki/table-2-model.png)
![Equations for Table 2](https://static.igem.wiki/teams/4460/wiki/eq-table-2.png)
Having estimated the different volume measurements, we were able to calculate the quantity of substrate molecules per liter according to the five different volume values recognized as levels. This body of water has controlled access, therefore, we were not able to get samples directly from the lagoon. We opted for utilizing five possible RDX concentration values to create Models A through E (see Table 3) for each one of them. The set of values selected for the implementation of our model were 1 mg/L, 2 mg/L, 3 mg/L, 4 mg/L, and 5 mg/L. These values were decided to be treated as molecules, hence the proper conversions were made using substrate characteristics (Table 4).
![Table 3: Substrate Molecules per Liter for Models A-E](https://static.igem.wiki/teams/4460/wiki/table-3-model.png)
![Equations for Table 3](https://static.igem.wiki/teams/4460/wiki/eq-table-3.png)
![Table 4: Characteristics](https://static.igem.wiki/teams/4460/wiki/table-4-model.png)
![Equations for Table 4](https://static.igem.wiki/teams/4460/wiki/eq-table-4.png)
The RDX concentration value established by the EPA as the maximum concentration of RDX that can be present in drinking water without causing adverse effects is 0.1 mg/L, as stated by the Agency for Toxic Substances and Disease Registry. This being said, it is important to know how much enzyme is needed to degrade RDX until reaching a concentration value that is considerably acceptable. We calculated the number of substrate (RDX) molecules per liter by multiplying each of the selected substrate concentration levels by the five levels of volume measurements mentioned previously (refer to Tables 5-9).
![Table 5: Model A - Substrate Molecules](https://static.igem.wiki/teams/4460/wiki/table-5-model.png)
![Table 6: Model B - Substrate Molecules](https://static.igem.wiki/teams/4460/wiki/table-6-model.png)
![Table 7: Model C - Substrate Molecules](https://static.igem.wiki/teams/4460/wiki/table-7-model.png)
![Table 8: Model D - Substrate Molecules](https://static.igem.wiki/teams/4460/wiki/table-8-model.png)
![Table 9: Model E - Substrate Molecules](https://static.igem.wiki/teams/4460/wiki/table-9-model.png)
![Equations for Table 5-9](https://static.igem.wiki/teams/4460/wiki/eq-table-5-9.png)
The evaluation of the degradation capacity of a single enzyme molecule was done by computing the quotient of the division of number of substrate molecules by the catalytic constant Kcat (s-1) , or turnover number. The catalytic constant Kcat describes the number of times each enzyme molecule will degrade substrate molecules into products per unit time (seconds). The Kcat value was assumed to be 3.33 s-1 , as determined in the article “The 1.5-Å Structure of XplA-heme, an Unusual Cytochrome P450 Heme Domain That Catalyzes Reductive Biotransformation of Royal Demolition Explosive” (Sabbadin et al., 2009). The experiments held for the article shared our project enzyme, substrate, and aerobic conditions. By completing this division for each of the five concentration levels and volumes, we obtained the amount of substrate molecules that should be degraded by a single enzyme molecule per second. Thus, this result brings insight into how capable our modified genetic circuit will be on achieving the detoxification of the Anones Lagoon.
Considering the magnitude and scope of the project, a rate measured in months for each model was used for further calculations (refer to Tables 5-9). Having this conversion, we were able to establish an approximated range of desired detoxification time in months. For each RDX molecule, enzyme molecules needed to remediate RDX were quantified. In Graphs 1-5, a tendency is noticed that as time progresses, less enzyme is produced because of the fact that less RDX would remain in the lagoon as breakage of RDX rings is occurring. The data showing this behavior is shown in Tables 10-14.
![Model A - Enzyme Molecules per months](https://static.igem.wiki/teams/4460/wiki/table-10-model-a.png)
![Model A - Enzyme Molecules vs time (months)](https://static.igem.wiki/teams/4460/wiki/model-a.png)
![Model B - Enzyme Molecules per months](https://static.igem.wiki/teams/4460/wiki/table-11-model-b.png)
![Model B - Enzyme Molecules vs time (months)](https://static.igem.wiki/teams/4460/wiki/model-b.png)
![Model C - Enzyme Molecules per months](https://static.igem.wiki/teams/4460/wiki/table-12-model-c.png)
![Model C - Enzyme Molecules vs time (months)](https://static.igem.wiki/teams/4460/wiki/model-c.png)
![Model D - Enzyme Molecules per months](https://static.igem.wiki/teams/4460/wiki/table-13-model-d.png)
![Model D - Enzyme Molecules vs time (months)](https://static.igem.wiki/teams/4460/wiki/model-d.png)
![Model E - Enzyme Molecules per months](https://static.igem.wiki/teams/4460/wiki/table-14-model-e.png)
![Model E - Enzyme Molecules vs time (months)](https://static.igem.wiki/teams/4460/wiki/model-e.png)
![Equations for Table 10-14](https://static.igem.wiki/teams/4460/wiki/eq-table-10-14.png)
Conclusion
Creating these five models helped in the understanding of the behavior of our devices. It was observed in our graphs that with the progression of time, the enzymatic production will decrease given that our genetic device will detect that less substrate remains to be degraded. It is shown in our model an approximation of the quantity of enzyme molecules needed to degrade a determined amount of substrate in a period of time. For future implementations, to optimize our model, we would calculate our experiment’s enzymatic efficiency, in other words, the number of enzyme molecules that are produced per liter of culture broth. By having a value for enzymatic efficiency of our degradation device, we will be able to calculate the amount of culture broth liters necessary to perform a complete RDX degradation in a selected period of time. Even though we were not able to obtain experimental results in the laboratory, we were able to make educated assumptions about certain parameters and predict the methodology that would be used to remediate RDX contamination in the Anones Lagoon.
References
Bender, Edward A. An introduction to mathematical modeling. Courier Corporation, 2000.
Google. (n.d.).[Anones Lagoon, Vieques, Puerto Rico]. Retrieved October 5, 2022 from
https://earth.google.com/web/@18.13760326,-65.29697124,5.37465706a,994.512221d,35y,
144.04842575h,44.99997742t,359.99999829r
Google. (n.d.).[Anones Lagoon, Vieques, Puerto Rico]. Retrieved October 5, 2022 from https://www.google.com/maps/place/Laguna+Anones/@18.1367371,-65.2973638,18z/data=!4m5!3m4!1s0x8c04c4c
87c900bfd:0x9ec7b50d6507a9ef!8m2!3d18.1379188!4d-65.2968399
Sabbadin, F., Jackson, R., Haider, K., Tampi, G., Turkenburg, J. P., Hart, S., Bruce, N. C., & Grogan, G. (2009). The 1.5-A structure of XplA-heme, an unusual cytochrome P450 heme domain that catalyzes reductive biotransformation of royal demolition explosive. The Journal of biological chemistry, 284(41), 28467–28475. https://doi.org/10.1074/jbc.M109.031559
Wu, B., R. Ako, and M. Hu. "A useful Microsoft Excel add-in program for modeling steady-state enzyme kinetics." Pharm. Anal. Acta S 11 (2011): 003. https://www.walshmedicalmedia.com/open-access/a-useful-microsoft-excel-add-in-program-for-modeling-steady-state-enzyme-kinetics-2153-2435.S11-003.pdf
© 2022 - Content on this site is licensed under a Creative Commons Attribution 4.0 International license.
The repository used to create this website is available at gitlab.igem.org/2022/rum-uprm.