Team:Northwestern/Modeling

From 2012.igem.org

(Difference between revisions)
Line 115: Line 115:
 +
 +
<table border="1">
 +
<tr>
 +
<td><b>Variable</b></td>
 +
<td><b>Description</b></td>
 +
<td><b>Units</b></td>
 +
</tr>
 +
<tr>
 +
<td>A_mRNA</td>
 +
<td>Messenger RNA for the ClC-ec1 antiporter concentration</td>
 +
<td>uM</td>
 +
</tr>
 +
<tr>
 +
<td>R</td>
 +
<td>gadR transcription factor concentration</td>
 +
<td>uM</td>
 +
</tr>
 +
<tr>
 +
<td>Pgad</td>
 +
<td>Pgad promoter (unactivated) concentration</td>
 +
<td>uM</td>
 +
</tr>
 +
<tr>
 +
<td>Pgad*</td>
 +
<td>Pgad/gadR activated promoter complex concentration</td>
 +
<td>uM</td>
 +
</tr>
 +
<tr>
 +
<td>G_mRNA</td>
 +
<td>GFP messenger RNA concentration</td>
 +
<td>uM</td>
 +
</tr>
 +
<tr>
 +
<td>GFP</td>
 +
<td>GFP protein concentration</td>
 +
<td>uM</td>
 +
</tr>
 +
</br>
 +
<table border="1">
 +
<tr>
 +
<td><b>Value</b></td>
 +
<td><b>Description</b></td>
 +
<td><b>Units</b></td>
 +
</tr>
 +
<tr>
 +
<td>k5</td>
 +
<td>R_mRNA degradation rate coefficient</td>
 +
<td>s-1</td>
 +
</tr>
 +
<tr>
 +
<td>k6</td>
 +
<td>R_mRNA to gadR protein translation rate coefficient</td>
 +
<td>s-1</td>
 +
</tr>
 +
<tr>
 +
<td>k7</td>
 +
<td>gadR protein degradation rate coefficient</td>
 +
<td>s-1</td>
 +
</tr>
 +
<tr>
 +
<td>k8</td>
 +
<td>
 +
Pgad/gadR activation binding rate coefficient</td>
 +
<td>1/M3•s</td>
 +
</tr>
 +
<tr>
 +
<td>k9</td>
 +
<td>Pgad/gadR activation dissociation rate coefficient</td>
 +
<td>s-1</td>
 +
</tr>
 +
<tr>
 +
<td>k11</td>
 +
<td>G_mRNA degradation rate coefficient</td>
 +
<td>s-1</td>
 +
</tr>
 +
<tr>
 +
<td>k12</td>
 +
<td>G_mRNA to GFP protein translation rate coefficient</td>
 +
<td>s-1</td>
 +
</tr>
 +
<tr>
 +
<td>k13</td>
 +
<td>GFP degradation rate coefficient</td>
 +
<td>s-1</td>
 +
</tr>
 +
<tr>
 +
<td>rR</td>
 +
<td>(micromolar) transcription rate of gadR</td>
 +
<td>uM/s</td>
 +
</tr>
 +
<tr>
 +
<td>rL</td>
 +
<td>maximum induced (molar) transcription rate of GFP</td>
 +
<td>uM/s</td>
 +
</tr>
 +
 +
</table>
 +
 +
 +
<p>The equations were solved with MATLAB’s ode15s, an ODE solver for ‘stiff’ functions. A full system plot can be seen below:
 +
 +
<div style="width:800px; margin-left:auto; margin-right:auto">
 +
<img src="https://static.igem.org/mediawiki/2012/5/51/Protein_production_new1.png" style="border:10px solid #333">
 +
</div>
 +
 +
<p>Here, we can see GFP levels (marked black) rising once Pgad becomes activated
 +
This also coincides with rising chloride levels, which continue to increase over time.
 +
 +
 +
 +
 +
 +
<div style="width:800px; margin-left:auto; margin-right:auto">
 +
<img src="https://static.igem.org/mediawiki/2012/1/10/PH_new1.png" style="border:10px solid #333">
 +
</div>
 +
<p>The pH of the cell quickly drops upon the modeled ‘entrance’ into the stomach. Interestingly, the pH of the cell quickly stabilizes to approximately its initial value (y0=7.5).
</div>
</div>
</html>
</html>

Revision as of 01:57, 27 October 2012

Modeling

Overview

The purpose of our modeled system is to tune this system as a specific pH-sensing lysis device for releasing our phytase into the stomach. We aim to characterize the interaction between the H+/Cl- antiporter and the chloride-induced lysis cassette in our E. coli chassis. This model examines the effects of varied promoter strengths for the ClC-ec1 antiporter, as well as the plausibility of utilizing the system for nutritional purposes.

We start with a model cell that has entered the stomach. In this high pH environment, protons begin to leak into the cell. For our purposes, this flux of H+ ions into the cell is modeled by simple diffusion.

The leakage of protons due to the acidic environment begins to disrupt pH homeostasis. However, E. coli has acid resistance mechanisms in order to grow in low pH. One of these mechanisms is the ClC-ec1 antiporter protein, which pumps H+ out of the cell (against its gradient) by utilizing the chloride gradient. The flux of extracellular protons (He) into the cell as well as the production of the antiporter protein are illustrated below:

Each antiporter turnover moves one H+ out while importing two Cl-. We now begin by modeling with first order ordinary differential equations.

Variable Description Units
A_mRNA Messenger RNA for the ClC-ec1 antiporter concentration uM
A Antiporter protein concentration uM
Hi Intracellular proton concentration uM
He Extracellular proton concentration uM
Cli Intracellular chloride concentration uM
Cle Extracellular chloride concentration uM
rA (molar) transcription rate of A_mRNA within cell uM per second

Value Description Units
k1 A_mRNA degradation coefficient s-1
k2 A_mRNA to protein A translation rate coefficient s-1
k3 protein A degradation rate coefficient s-1
k4 antiporter kinetic coefficient /M3•s

The equations that model these species in our system are as follows:

The extracellular proton and chloride rates of change are set to zero since the extracellular volume (the human stomach) is as much as 15 orders of magnitude larger than the individual cell volume. Thus molecular fluxes will have negligible effects on the ionic concentrations.

Our Phytastic E. coli use a chloride-inducible promoter, Pgad, to sense the change in chloride concentration as a way to detect the change in pH. The Pgad operon part, which we obtained from the Chinese University of Hong Kong iGem team, includes a constitutive promoter for a positive regulator gene, gadR. gadR binds upstream of the Pgad promoter, where it can be allosterically activated by chloride ions, which significantly increases transcription. Our Phytastic E. coli have a coding region for lysozymes downstream of the Pgad part so that they can selectively lyse at low pH. The model, however, replaces the lysis cassette with GFP for the purpose of comparing the model to results. The Pgad/GFP system is illustrated below:


Variable Description Units
A_mRNA Messenger RNA for the ClC-ec1 antiporter concentration uM
R gadR transcription factor concentration uM
Pgad Pgad promoter (unactivated) concentration uM
Pgad* Pgad/gadR activated promoter complex concentration uM
G_mRNA GFP messenger RNA concentration uM
GFP GFP protein concentration uM
Value Description Units
k5 R_mRNA degradation rate coefficient s-1
k6 R_mRNA to gadR protein translation rate coefficient s-1
k7 gadR protein degradation rate coefficient s-1
k8 Pgad/gadR activation binding rate coefficient 1/M3•s
k9 Pgad/gadR activation dissociation rate coefficient s-1
k11 G_mRNA degradation rate coefficient s-1
k12 G_mRNA to GFP protein translation rate coefficient s-1
k13 GFP degradation rate coefficient s-1
rR (micromolar) transcription rate of gadR uM/s
rL maximum induced (molar) transcription rate of GFP uM/s

The equations were solved with MATLAB’s ode15s, an ODE solver for ‘stiff’ functions. A full system plot can be seen below:

Here, we can see GFP levels (marked black) rising once Pgad becomes activated This also coincides with rising chloride levels, which continue to increase over time.

The pH of the cell quickly drops upon the modeled ‘entrance’ into the stomach. Interestingly, the pH of the cell quickly stabilizes to approximately its initial value (y0=7.5).