Line 251: | Line 251: | ||
</p> | </p> | ||
<center><div class="card mt-3 mb-3" style="width:500px"> | <center><div class="card mt-3 mb-3" style="width:500px"> | ||
− | <img class="card-img-top" src="" alt="Card image"> | + | <img class="card-img-top" src="https://static.igem.org/mediawiki/2019/7/7f/T--Mingdao--1.png" alt="Card image"> |
<div class="card-body"> | <div class="card-body"> | ||
<p class="card-text"><strong>Equation 7. </strong>Henry's Law, modeling the relationship between the partial pressure of gas and its concentration in a solvent. | <p class="card-text"><strong>Equation 7. </strong>Henry's Law, modeling the relationship between the partial pressure of gas and its concentration in a solvent. |
Revision as of 06:49, 18 October 2019
![Responsive image](https://static.igem.org/mediawiki/2019/d/d1/T--Mingdao--Cair%284%29.png)
Objectives
Factor 1
Factor 2
Final Model
References
Objectives
The objectives of our mathematical model are to:
- Demonstrate and prove the feasibility of our project design by showing how adding human carbonic anhydrase II (abbreviated as CA) affects the consumption of carbon dioxide by the cyanobacteria, synechococcus elongatus, as opposed to natural conditions without CA.
- Predict how the indoor CO2 concentration changes over time, given the number of people in the room, the exact volume of the room, and whether or not CA is utilized.
- Determine the theoretical maximum number of people allowed in a room to maintain the indoor CO2 concentration under a certain level, with our product utilized.
Our model is based on two major assumptions. For the sake of simplicity, we assume that the space under examination is enclosed and airtight, meaning that inflow to and outflow from the room are both eliminated. In addition, regarding the air conditioning, which is a major source of air flow, since internal air is recirculated without interacting with external air with most modern air conditioning systems, it can be assumed that air conditioning does not affect the indoor carbon dioxide level.
OVERVIEW OF MODEL PARAMETERS
Several factors influence the CO2 concentration in an enclosed space (for example: a meeting room), namely the following:
- Additive:
- The number of people in the space.
- The average CO2 exhalation rate per person.
- The average CO2 concentration in exhale.
- Subtractive:
- The natural dissolution rate of CO2 in water (when CA is absent).
- The dissolution rate of CO2 in water when CA is present, expressed as the enzymatic reaction rate of CA.
- The rate of CO2 consumption from the photosynthesis of cyanobacteria.
- The solubility of CO2 in water at various acidities.
- Other:
- The volume of the enclosed space.
- Initial indoor CO2 concentration.
We obtained literature values and models for Factors 1b, 1c, and 2. Factors 1a and 3 are treated as independent variables in the model that can be input accordingly to calculate the indoor CO2 concentration.
Factor 1: Modeling the CO2 produced by human exhalation
In order to come up with a comprehensive mathematical model, we decided to start by modeling the indoor CO2 concentration over time, considering human breathing. The indoor CO2 concentration in this case can be modeled with the following formula (Jung).
![Card image](https://static.igem.org/mediawiki/2019/7/78/T--Mingdao--equa11.png)
Equation 1. The governing equation modeling the concentration of carbon dioxide (C) in an enclosed space at a specific time (t), given the number of people in the room (n), the CO2 concentration in exhale (Cex), the exhalation flow rate (Qex), the volume of the room (Vroom), and the initial indoor CO2 concentration (C0).
In order to better combine this equation with other factors of our model, we rewrote Equation 1 into its differential form and defined it as the instantaneous rate of human production of CO2 by exhalation at time t, vhuman(t).
![Card image](https://static.igem.org/mediawiki/2019/2/2f/T--Mingdao--equa12.png)
Equation 2. The instantaneous change in indoor CO2 concentration (derivative of Equation 1).
Factor 1b: Determining the Exhalation Rate of an Individual (Qex)
Exhalation rate is the amount of air that an individual exhales per unit time. From literature, the exhalation rate of an individual is shown below (Jung). The default value used in our model is calculated by the average of minimum and maximum values. Unit conversions are shown below.
![](https://static.igem.org/mediawiki/2019/e/e3/T--Mingdao--table11.png)
Table 1. Exhalation rate of an individual with unit conversion(Jung)
![](https://static.igem.org/mediawiki/2019/0/0b/T--Mingdao--equa13.png)
Factor 1c: Determining the CO2 Concentration in Exhale (Cex)
The CO2 concentration of an individual's breathing varies highly and is dependent on the individual and their condition. The CO2 concentration in exhale from literature is shown below (Jung). The average of minimum and maximum values is used as the default value in our model. Unit conversions are shown below, with assumptions of an environment with an air pressure of 100 kPa and an air temperature of 25 ℃ (or 297 K).
![](hhttps://static.igem.org/mediawiki/2019/0/06/T--Mingdao--modre.png)
Table 2. CO2 concentration in exhale with unit conversion (Jung)
![](https://static.igem.org/mediawiki/2019/0/08/T--Mingdao--equa14.png)
Factor 2: Modeling the Factors that Reduce Carbon Dioxide
The CO2 fixation process can be illustrated with the following diagram.
![](https://static.igem.org/mediawiki/2019/0/0c/T--Mingdao--modelpic1.png)
Therefore, the process of defining the sequestration of CO2 in our model can be divided into two aspects: determining the natural dissolution rate of CO2 in water (when CA is absent) and determining the increased dissolution rate of CO2 in water (when CA is present). However, the dissolved inorganic carbon (DIC; which implies the combination of H2CO3, HCO3-, CO32-, and all other forms of dissolved carbon dioxide) in water can eventually reach the maximum capacity, which is itself dependent on the pH value of the environment. How rapidly this theoretical maximum is attained is dictated by the photosynthetic rate of cyanobacteria.
Factor 2a: Determining the Natural Dissolution Rate of CO2 in Water
The natural dissolution rate of CO2 in water (vnatural) is calculated using the first-order rate constant equation, as shown in Equation 3, where 0.039 s-1 is the rate constant for the forward reaction of the hydration of carbon dioxide (CO2+H2O->H2CO3). The correlation between the concentration of gaseous carbon dioxide and the natural dissolution rate is illustrated in Figure 1.
![Card image](https://static.igem.org/mediawiki/2019/c/c4/T--Mingdao--equa15.png)
Equation 3. First-order rate constant equation for hydration of carbon dioxide.
![](https://static.igem.org/mediawiki/2019/8/85/T--Mingdao--modelfig1.png)
Figure 1. The linear curve modeling the correlation between the concentration of gaseous carbon dioxide and its natural dissolution rate in water.
Factor 2b: Determining the Enzymatic Reaction Rate of CA
We modeled the enzymatic rate of carbonic anhydrase with Michaelis-Menten kinetics, as shown in Equation 4. vCA represents the enzymatic reaction rate of CA; [CO2(g)] represents the concentration of substrate (in this case, gaseous carbon dioxide); vmax represents the maximum reaction rate possible; and Km is the Michaelis-Menten constant, which is defined as the substrate concentration when the enzymatic reaction rate is half of its maximum (whenv=1/2vmax).
![Card image](https://static.igem.org/mediawiki/2019/2/27/T--Mingdao--equa16.png)
Equation 4. The Michaelis-Menten equation, which models the relationship between the concentration of gaseous carbon dioxide and the enzymatic reaction rate.
![Card image](https://static.igem.org/mediawiki/2019/f/f4/T--Mingdao--equa17.png)
Equation 5. Calculation for vmax, the maximum reaction rate. Kcat denotes the catalyst rate constant of CA and [CA] represents the concentration of the enzyme.
Therefore, by inserting the vmax and Km values into the model, the correlation between CO2 concentration and the reaction rate of CA can be shown. vmax can be calculated with Equation 5.
We obtained the Kcat and Km values for the human carbonic anhydrase II (hCAII) from literature. According to Silverman et al., the Kcat and Km values for hCAII are 1.4×106 s-1 and 9.3×10-3 M, respectively. A [CA] value of 2.883×10-3 M is used based on the enzyme we use. Unit conversions are shown below.
![](https://static.igem.org/mediawiki/2019/b/b1/T--Mingdao--equa18.png)
Using the enzymatic coefficients for the Michaelis-Menten model, we plotted the curve in Figure 2 to model the enzyme kinetics of human carbonic anhydrase II.
![](https://static.igem.org/mediawiki/2019/9/9d/T--Mingdao--modelfig2.png)
Figure 2. The Michaelis-Menten curve modeling the enzyme activities of hCAII at various concentrations of carbon dioxide.
Factor 2c: Determining the Photosynthetic Rate of Cyanobacteria in the Form of CO2 Consumption
After carbon dioxide has fully hydrated and dissolved into water, the last phase in our process of CO2 sequestration is its removal by the photosynthesis of cyanobacteria Synechococcus elongatus PCC7942. The rate at which the dissolved CO2 is consumed by our cyanobacteria can be modeled with Equation 5, which models the relationship between the concentration of dissolved inorganic carbon (DIC) and the rate of carbon-specific growth (Cμ; which is also equivalent to the rate of photosynthesis vphotosynthesis), which is represented in the form of O2 production (Clark). The conversion from the amount of O2 production to CO2 sequestration is then performed using the 1:1 molar ratio between O2 and CO2 obtained from the chemical equation of photosynthesis. The maximum rate of growth μmax and the half saturation constant for Cμ, KG(DIC), specific for the cyanobacteria Synechococcus are then consulted from literature. Calculations for the final function are shown below with the parameters listed in Table 3.
![Card image](https://static.igem.org/mediawiki/2019/2/2c/T--Mingdao--equa19.png)
Equation 6. Relationship between the rate of photosynthesis (vphotosynthesis) and the concentration of dissolved inorganic carbon (DIC) (Clark).
![](https://static.igem.org/mediawiki/2019/6/61/T--Mingdao--modelfig3.png)
Table 3. Parameters used in calculation for Cμ
![](https://static.igem.org/mediawiki/2019/c/c3/T--Mingdao--equa20.png)
Using the converted photosynthetic coefficients that we obtained through literature for our cyanobacteria Synechococcus elongatus, we constructed a graph in Figure 3 showing the relationship between the concentration of dissolved inorganic carbon in water and the rate of CO2 consumption due to photosynthesis.
![](https://static.igem.org/mediawiki/2019/e/e6/T--Mingdao--modelfig4.png)
Figure 3. The Michaelis-Menten curve modeling the photosynthetic activity of Synechococcus elongatus. The CO2 consumption rate at various concentrations of dissolved inorganic carbon (DIC) is shown.
Factor 2d: Determining the Solubility of CO2 in Alkaline Environment
In order to model the maximum dissolvable amount of carbon dioxide in alkaline environments, we first used the Henry's Law (Equation 7), which states that the concentration (C) of a gas in liquid is equal to its partial pressure (Pgas) multiplied by the Henry's law solubility constant (Hcp).
![Card image](https://static.igem.org/mediawiki/2019/7/7f/T--Mingdao--1.png)
Equation 7. Henry's Law, modeling the relationship between the partial pressure of gas and its concentration in a solvent.
From Sander (2015), the Henry's law constant (Hcp) of carbon dioxide at 25 °C with water as a solvent is 3.3 ⨉ 10-4 mol·m-3·Pa-1. Unit conversion is shown below.
Assuming the concentration of CO2 in air is 400 ppm and the air pressure is 1 atm, the concentration of CO2 in water can be calculated as the following:
With the rates of the solubility of CO2 identified both with or without the addition of CA, it is crucial to note that an enzyme does not increase the final yield of products, but instead accelerates its reaction rate. Therefore, it is critical to identify the maximum attainable concentration in liquid. It is also essential to note that the solvent we use in our project is a culture medium for Synechococcus elongatus and is alkaline with a pH value of 11. The effect of pH on solubility is summarized by the following equation, with the maximum concentration from Henry’s Law considered as [CO2 (aq)]:
Equation 8. An extension of the Henry’s Law that determines the maximum solubility of carbon dioxide in water considering the acidity of the reaction environment.
Table 4. Parameters used in calculation for the maximum attainable concentration of carbon dioxide in water Cmax in alkaline environments.Therefore the calculation shows:
Our Final Model
![](https://static.igem.org/mediawiki/2019/f/f1/T--Mingdao--image31.png)
Before combining all the equations together, one thing to note is that our prototype is a closed system as shown above. Therefore, a concentration variable Cdev is defined to denote the concentration of gaseous carbon dioxide in the system. Qin and Qout refer to the inflow and outflow rate of this closed system, while Cin and Cout refer to the inflow and outflow concentrations (Cin is equivalent to the ambient carbon dioxide concentration, Cair). The concentration of gaseous carbon dioxide in the system can be modeled with the differential equation, Equation 10. The rate at which carbon dioxide dissolves in water, vdissolved, is a function of the in-device carbon dioxide concentration (Cdev) and is expressed conditionally, with different definitions depending on whether or not the maximum dissolution capacity is attained, and whether or not CA is added. The definition for vdissolved is shown in Equation 9.
Equation 9. The equation governing the rate at which carbon dioxide dissolves in water, which is a function of Cdev. (From Equations 3, 4, 8)
Equation 10. The equation governing the instantaneous change in in-device gaseous CO2 concentration. (From Equations 6, 9, 10)
Equation 11. The equation governing the instantaneous change in in-water CO2 concentration. (From Equations 6, 9, 10, 11)
Equation 12. The final model that demonstrates the instantaneous change in CO2 concentration in the air, with human exhalation, photosynthesis, and enzyme activity of CA taken into consideration. Vroom signifies the volume of the room modeled while Vwater indicates the in-device volume of the liquid that contains cyanobacteria (From Equations 2, 9, 10)
Our Software(python graphing program)Evaluation of Model
Impact on Experimental Design. In order to testify the model we developed against reality, an experiment was designed to investigate the variables identified through the model. In the experiment, we set the concentration of CO2 and flow rate as the independent variables to test its effect on carbon dioxide sequestration and maintained the volume of cyanobacteria’s biomass and the device as controlled variables to create a simulation on the workings of our device in real life. More details of the experiment are shown in the “Experiment” page.
References
- Clark, D. R., and K. J. Flynn. “The Relationship between the Dissolved Inorganic Carbon Concentration and Growth Rate in Marine Phytoplankton.” Proceedings of the Royal Society of London. Series B: Biological Sciences, vol. 267, no. 1447, 2000, pp. 953–959., doi:10.1098/rspb.2000.1096.
- Collier, Jackie L., et al. “Changes in the Cyanobacterial Photosynthetic Apparatus during Acclimation to Macronutrient Deprivation.” Photosynthesis Research, vol. 42, no. 3, 1994, pp. 173–183., doi:10.1007/bf00018260.
- “Dissolved Carbon Dioxide”, http://ion.chem.usu.edu/~sbialkow/Classes/3650/CO2 Solubility/DissolvedCO2.html.
- Jung, Heejung. “Modeling CO2 Concentrations in Vehicle Cabin.” SAE Technical Paper Series, 2013, doi:10.4271/2013-01-1497.
- Keren, Nir, et al. “A Light-Dependent Mechanism for Massive Accumulation of Manganese in the Photosynthetic BacteriumSynechocystissp. PCC 6803†.” Biochemistry, vol. 41, no. 50, 2002, pp. 15085–15092., doi:10.1021/bi026892s.
- Sander, R. “Compilation of Henrys Law Constants (Version 4.0) for Water as Solvent.” Atmospheric Chemistry and Physics, vol. 15, no. 8, 2015, pp. 4399–4981., doi:10.5194/acp-15-4399-2015.
- Silverman, David N., and Sven Lindskog. “The Catalytic Mechanism of Carbonic Anhydrase: Implications of a Rate-Limiting Protolysis of Water.” Accounts of Chemical Research, vol. 21, no. 1, 1 Jan. 1988, pp. 30–36., doi:10.1021/ar00145a005.
- Turpin, David H., et al. “Predicting The Kinetics Of Dissolved Inorganic Carbon Limited Growth From The Short-Term Kinetics Of Photosynthesis In Synechococcus Leopoliensis (Cyanophyta)1.” Journal of Phycology, vol. 21, no. 3, 2004, pp. 409–418., doi:10.1111/j.0022-3646.1985.00409.x.