Difference between revisions of "Team:ECUST/Iron Sensing"

Line 66: Line 66:
 
<p>K<sub>i1</sub> is a non-dimensional parameter which repesents the inhibition power, and K<sub>f</sub> is the fixation rate of the Fe-FUR on the FBS. N<sub>pla</sub>1 is the number of pasmids containing the sensor system.  </p>
 
<p>K<sub>i1</sub> is a non-dimensional parameter which repesents the inhibition power, and K<sub>f</sub> is the fixation rate of the Fe-FUR on the FBS. N<sub>pla</sub>1 is the number of pasmids containing the sensor system.  </p>
 
<h2>1.3 Result</h2>
 
<h2>1.3 Result</h2>
 +
<p>We want to know the fittest ki1 for the model to sense the iron and the concentrate of iron.
 +
We make three kinds of fur-box for our sensor system. We want to know which is our best choice. Our experiment result show in the figure 3. </p>
 +
<p>Then we try to change the value of K<sub>i1</sub> to model different strength of promoter with fur-box in our experiment which show in the figure 4. We want our system to make sense in the high level of Fe<sup>2+</sup>, so we choose the fur-2 system. And we finally set the KI1:6.3*10<sup>-5</sup>. </p>
  
  
Line 74: Line 77:
 
</div>
 
</div>
  
<div class="contentbox">
 
<h1 class="box-heading">3.Goals</h1>
 
<p>&nbsp;&nbsp;Our goal of this model is to create a generic quorum sensing model so that:</p>
 
<p>&nbsp;&nbsp;• We can determine the effect of afeR promoter and predict the production of DspB and enterobactin.</p>
 
<p>&nbsp;&nbsp;• We can predict hao long our engineered bacteria would take to remove the biofilm and rust.</p>
 
</div>
 
 
<div class="contentbox">
 
<h1 class="box-heading">4.Materials and Methods</h1>
 
<h2>4.1 HSL Transfer</h2>
 
<p>&nbsp;&nbsp;HSL is produced by iron bacterias and realeased into the water environment. So the first step of our sensing is HSL transfering into our engineered E.coli from the water. And a passive transusion model is used for this process that the transfer rate of HSL can be described as this:</p>
 
<p>&nbsp;&nbsp;• K<sub>HSL,W-C</sub> : transfer coefficient through the membrane (s−1)</p>
 
<p>&nbsp;&nbsp;• We can predict hao long our engineered bacteria would take to remove the biofilm and rust.</p>
 
 
<h2>4.2 AfeR-HSL Complexation</h2>
 
<p>&nbsp;&nbsp;AfeR is produced by engineered E.coli and functions in cell and its concentration is obtained approximating the number of protein per cell, using the E.coli concentration (cell/L) and the Avogadro number.</p>
 
<p>&nbsp;&nbsp;The AfeR-HSL complexation is simply formed that way:</p>
 
<p>&nbsp;&nbsp;Assuming kinetics of AfeR-HSL complexation complexation is fast compared to the rest of the system, we assumed that the free and complexed forms are at equilibrum.</p>
 
<p>&nbsp;&nbsp;• K <sub>eq, AfeR-HSL</sub> : equilibrum constant of the AfeR-HSL complexation (mol/L)</p>
 
 
<h2>4.3 DspB Production</h2>
 
<p>&nbsp;&nbsp;The production of the DspB from the DspB gene includes transcription and translation after activation. In addition, we should also consider its transport and degradation.</p>
 
<h3>4.3.1 DspB Gene Activation</h3>
 
<p>&nbsp;&nbsp;This process is modeled using a Michaelian formalism depending on its activator (AfeR-HSL complexation) concentration. The promoter strength is also taken into account.</p>
 
<p>&nbsp;&nbsp;• DspB <sub>DNA,0/cell</sub> : total number of DspB DNA per cell</p>
 
<p>&nbsp;&nbsp;• DspB <sub>DNA/cell</sub> : number of activated DspB DNA per cell</p>
 
<p>&nbsp;&nbsp;• K <sub>a, AfeR-HSL</sub> : activation constant of the AfeR-HSL complexation (mol/L)</p>
 
<p>&nbsp;&nbsp;• k <sub>p, afeR</sub> : afeR promoter influence</p>
 
 
<h3>4.3.2 DspB Transcription</h3>
 
<p>&nbsp;&nbsp;The DspB transcription depends on the transcription rate of the strain and the length of the DspB gene. The Avogadro number is used to express the transcription velocity in molar concentration in one cell per time unit.</p>
 
<p>&nbsp;&nbsp;• k<sub>transcript</sub> : E.coli transcription rate (nucleotides/s)</p>
 
<p>&nbsp;&nbsp;• RNA polymerase/gene: number of RNA polymerase per gene</p>
 
<p>&nbsp;&nbsp;• DNA length (DspB): number of nucleotides on the DspB gene</p>
 
<p>&nbsp;&nbsp;• V <sub>intracell</sub>: volume of a bacterial cell (L)</p>
 
<p>&nbsp;&nbsp;For the convenience of mathematical operation, we merged the ktranscript、RNA polymerase/gene and "V" intracell to a constant.</p>
 
 
<h3>4.3.3 DspB Translation</h3>
 
<p>&nbsp;&nbsp;The DspB translation depends on the translation rate of the strain, the mRNA length and the quantity of mRNA. The translation velocity is expressed in molar concentration in one cell per time unit.</p>
 
<p>&nbsp;&nbsp;• k<sub>translation</sub> : E.coli translation rate (nucleotides/s)</p>
 
<p>&nbsp;&nbsp;• Ribosomes/RNA: number of ribosomes per mRNA</p>
 
<p>&nbsp;&nbsp;• RNA length (DspB): number of nucleotides on the DspB mRNA</p>
 
<p>&nbsp;&nbsp;• [DspB mRNA] : DspB mRNA concentration in one E.coli cell</p>
 
<p>&nbsp;&nbsp;For the convenience of mathematical operation, we merge the ktranslation and Ribosomes/RNA and to a constant.</p>
 
 
<h3>4.3.4 Degradation</h3>
 
<p>&nbsp;&nbsp;Some of the DspB protein and mRNA are degraded. A degradation constant is used to model the degradation velocity.</p>
 
<p>&nbsp;&nbsp;• K<sub>deg,DspB</sub>: DspB degradation constant (s−1)</p>
 
<p>&nbsp;&nbsp;• K<sub>deg,DspB mRNA</sub>: DspB mRNA degradation constant (s−1)</p>
 
 
<h3>4.3.5 DspB Transfer</h3>
 
<p>&nbsp;&nbsp;DspB protein needs to be transferred to the water environment to function. This process is taken into account through a passive transusion model.</p>
 
<p>&nbsp;&nbsp;• K<sub>DspB,C-W</sub> : transfer coefficient through the membrane (s−1)</p>
 
 
 
 
<h2>4.4 Biofilm Removel</h2>
 
<p>&nbsp;&nbsp;The biofilm is removed by the DspB and the process is modeled assuming a Michaelis-Menten kinetics.</p>
 
<p>&nbsp;&nbsp;• kcat,DspB : catalytic constant of the DspB enzyme (s −1)</p>
 
<p>&nbsp;&nbsp;• "KM,D"  : Michaelis constant of the DspB enzyme (mol/L)</p>
 
 
 
<h2>4.5 EntE Production</h2>
 
<p>&nbsp;&nbsp;We treat enterobactin enzymes gene cluster as a whole gene (EntE gene). The production of the enterobactin enzymes from the EntE gene includes transcription and translation after activation. In addition, we should also consider its degradation. Because the enterobactin enzymes function in the cell, we don't need to consider its transport to the water environment.</p>
 
<h3>4.5.1 EntE Gene Activation</h3>
 
<p>&nbsp;&nbsp;This process is modeled using a Michaelian formalism depending on its activator (AfeR-HSL complexation) concentration. The promoter strength is also taken into account.</p>
 
<p>&nbsp;&nbsp;• EntE DNA,0/cell : total number of EntE DNA per cell</p>
 
<p>&nbsp;&nbsp;• EntE DNA/cell : number of activated EntE DNA per cell</p>
 
<p>&nbsp;&nbsp;• K a, AfeR-HSL : activation constant of the AfeR-HSL complexation (mol/L)</p>
 
<p>&nbsp;&nbsp;• k p, afeR : afeR promoter influence</p>
 
 
<h3>4.5.2 EntE Transcription</h3>
 
<p>&nbsp;&nbsp;The EntE transcription depends on the transcription rate of the strain and the length of the EntE gene. The Avogadro number is used to express the transcription velocity in molar concentration in one cell per time unit.</p>
 
<p>&nbsp;&nbsp;• "EntE DNA,/cell"  : number of EntE gene per cell</p>
 
<p>&nbsp;&nbsp;• ktranscript : E.coli transcription rate (nucleotides/s)</p>
 
<p>&nbsp;&nbsp;• RNA polymerase/gene: number of RNA polymerase per gene</p>
 
<p>&nbsp;&nbsp;• DNA length (EntE): number of nucleotides on the EntE gene</p>
 
<p>&nbsp;&nbsp;• Vintracell : volume of a bacterial cell (L)</p>
 
<p>&nbsp;&nbsp;For the convenience of mathematical operation, we merged the ktranscript、RNA polymerase/gene and "V" intracell to a constant.</p>
 
 
 
</div>
 
  
 
</main>
 
</main>

Revision as of 14:20, 17 October 2018

Iron Model

The fur-lac-cecropin AD system

In order to measure the ability of engineered bacteria to kill iron bacteria, we modeled the production of the cecropin AD and the sterilizing ability of the cecropin. We want to model the fur-box fell the concentrate of Fe2+ to drive the LacI to reverse the signal to produce the cecropin AD which play a bactericidal effect. So our model used in the project includes four parts: the iron sensor, the inverter system, the sterilizing system and the chelator system.

The first part focuses on the sensor system to find the Fe2+ needed by our system. And the strength of the best promoter with the fur-box of the three kinds of fur-box was determined. Second, the inverter system our team implemented in the bacteria focused on the cecropin AD produced with time. The third part modeled the sterilization rate of the cecropin AD. The last part integrate those previous parts leading to the cecropin production and sterilization. The model show the time it takes for our bacteria system to work.

1.Sensor model

1.1 Introduction

We first modeled the sensing system using ODEs with the help of experimental results to determine one of our parameters ki1.We had three kinds of fur-box designs. We model our three kinds of fur-box (shown in the figure 2) to find the optimal fur-box and the strength of the promoter with the best kind of fur-box. We finally corrected our model through the experiments. We make this framework like figure 1.

1.2 Methods and materials:

1.2.1 The dynamic simulation of sense iron to FBS:

(1) the iron-FUR complex formation:

$$2\cdot FUR+2\cdot Fe\leftrightarrow Fe_{2}FUR_{2}$$

We think this equation to:

$$FUR+Fe\leftrightarrow FeFUR:K_{FeFUR}$$

We just want to use differential equations more easily. And we can easily divide our [FeFur] by two to get the real complex concentration.

(2) We can easily make the formation (v) and the dissociation (v') speeds:

$$V=K_{FeFUR}\cdot \left[ FUR\right] \cdot \left[ Fe\right]$$

$$V'=d_{ff}\cdot \left[ FeFUR\right] $$

• KFeFUR : Formation constant of FeFur complex (m-1∙s-1)

• dff : FeFUR degradation rate (min-1)

We model the iron input in the bacteria using a linear function of the external iron concentration Ferext with the factor p which is the cell-wall permeability for iron.

$$\dfrac {d\left[ Fe\right] }{dt}=p\cdot Ferext-K_{FeFUR}\cdot \left[ FUR\right] \cdot \left[ Fe\right] +d_{ff}\cdot \left[ FeFUR\right]$$

$$\dfrac {d\left[ FUR\right] }{dt}=FurO-K_{FeFUR}\cdot \left[ FUR\right] \cdot \left[ Fe\right] +d_{ff}\cdot \left[ FeFUR\right]$$

• p : Permeability of cell wall (min-1)

• dff : FeFUR degradation rate (min-1)

We track the free Fe-FUR complex but not those attached to a Fur Binding Sites in our model.

$$\dfrac {d\left[ FeFUR\right] }{dt}=K_{FeFUR}\cdot \left[ FUR\right] \cdot \left[ Fe\right] -d_{ff}\cdot \left[ FeFUR\right] -\dfrac {1}{N_{A}V}.\dfrac {dFBS}{dt}$$

• NA : Avogadro’s constant (mol-1)

• V : Volume of a bacterium (m3)

• FBS : the number of inhibited Fur Binding Sites

We use our Logistic function under its differential form to simulate the inhibition phenomenon. Since it is the Fe-FUR that represses it, the LacI can be expressed as a logistic fuction of the Fe-FUR:

$$\dfrac {dFBS}{d\left[ FeFUR\right] }=\dfrac {K_{i}1}{K_{f}}\cdot FBS\left( \left[ FeFUR\right] \right) \cdot \left( 1-\dfrac {FBS\left( \left[ FeFUR\right] \right) }{N_{pla1}}\right) $$

• Kf : fixation rate of FeFUR (min-1)

• Ki1 : constant repesents the inhibition power (min-1)

• Npla1 : pET28-a plasmid number (nb/cell)

Ki1 is a non-dimensional parameter which repesents the inhibition power, and Kf is the fixation rate of the Fe-FUR on the FBS. Npla1 is the number of pasmids containing the sensor system.

1.3 Result

We want to know the fittest ki1 for the model to sense the iron and the concentrate of iron. We make three kinds of fur-box for our sensor system. We want to know which is our best choice. Our experiment result show in the figure 3.

Then we try to change the value of Ki1 to model different strength of promoter with fur-box in our experiment which show in the figure 4. We want our system to make sense in the high level of Fe2+, so we choose the fur-2 system. And we finally set the KI1:6.3*10-5.