Team:USTC-China/modeling/gene

From 2014.igem.org

(Difference between revisions)
 
(One intermediate revision not shown)
Line 51: Line 51:
     <div class="large-9 columns" id="content-page">
     <div class="large-9 columns" id="content-page">
       <div class="title"><h1>Gene Network Modeling</h3>
       <div class="title"><h1>Gene Network Modeling</h3>
-
  <div data-magellan-expedition="fixed">
+
    <div data-magellan-expedition="fixed">
-
      <dl class="sub-nav">
+
        <dl class="sub-nav">
-
      <dd data-magellan-arrival="introduction"><a herf="#introduction">Introduction</a></dd>
+
          <dd data-magellan-arrival="introduction"><a href="#introduction">Introduction</a></dd>
-
  <dd data-magellan-arrival="network"><a herf="#network">Network</a></dd>
+
        <dd data-magellan-arrival="network"><a href="#network">Network</a></dd>
-
  <dd data-magellan-arrival="dnainverseenzyme"><a herf="#dnainverseenzyme">DNA Inverse Enzyme</a></dd>
+
        <dd data-magellan-arrival="dnainverseenzyme"><a href="#dnainverseenzyme">DNA Inverse Enzyme</a></dd>
-
  <dd data-magellan-arrival="dgradgrbcircuit"><a herf="#dgradgrbcircuit">DgrA/DgrB</a></dd>
+
        <dd data-magellan-arrival="dgradgrbcircuit"><a href="#dgradgrbcircuit">DgrA/DgrB</a></dd>
-
  <dd data-magellan-arrival="hfiahfsjcircuit"><a herf="#hfiahfsjcircuit">HfiA/HfsJ</a></dd>
+
        <dd data-magellan-arrival="hfiahfsjcircuit"><a href="#hfiahfsjcircuit">HfiA/HfsJ</a></dd>
-
</dl>
+
      </dl>
-
  </div>
+
    </div>
-
 
+
   
-
  </div>
+
    </div>
       <div class="text">
       <div class="text">
Line 71: Line 71:
         <a name="introduction"></a>
         <a name="introduction"></a>
-
<h2 data-magellan-destination="introduction">Introduction<h2>
+
    <h2 data-magellan-destination="introduction">Introduction<h2>
         <p>We've constructed several modeling skeleton in Matlab Simbiology Package, they are:</p>
         <p>We've constructed several modeling skeleton in Matlab Simbiology Package, they are:</p>
Line 82: Line 82:
         </ol>
         </ol>
-
         <a name="networks"></a>
+
         <a name="network"></a>
-
<h2 data-magellan-destination="network">Network<h2>
+
    <h2 data-magellan-destination="network">Network<h2>
         <a name="dnainverseenzyme"></a>
         <a name="dnainverseenzyme"></a>
-
<h3 data-magellan-destination="dnainverseenzyme">DNA Inverse Enzyme<h3>
+
    <h3 data-magellan-destination="dnainverseenzyme">DNA Inverse Enzyme<h3>
         <h4 id="introduction">Introduction</h4>
         <h4 id="introduction">Introduction</h4>
Line 113: Line 113:
         <a name="dgradgrbcircuit"></a>
         <a name="dgradgrbcircuit"></a>
-
<h3 data-magellan-destination="dgradgrbcircuit">DgrA/DgrB Circuit<h3>
+
    <h3 data-magellan-destination="dgradgrbcircuit">DgrA/DgrB Circuit<h3>
Line 136: Line 136:
         <a name="hfiahfsjcircuit"></a>
         <a name="hfiahfsjcircuit"></a>
-
    <h3 data-magellan-destination="hfiahfsjcircuit">HfiA/HfsJ Circuit<h3>
+
      <h3 data-magellan-destination="hfiahfsjcircuit">HfiA/HfsJ Circuit<h3>

Latest revision as of 14:51, 17 October 2014

In this part of modeling work, we mainly center on constructing some crucial circuits in our regulating network, including the light responding upstreams and downstream movement control mechanisms, as well as some middle wares.

The aim of these construction work is to know about the how the messaging delivering system work as a whole, how does the different signals effect each other in the final output level.

Introduction

We've constructed several modeling skeleton in Matlab Simbiology Package, they are:

  1. DgrA/DgrB Circuit in regulating the flagellum movement
  2. HfiA/HfsJ Circuit in regulating the genesis of composing substances holdfast.
  3. The DNA Inverse Circuit which may change the direction in which promoter works.
  4. Light Sensors including Blue, Red and Green.

Network

DNA Inverse Enzyme

Introduction

In this model, our aim is to verify that according to use of a recombinase, we can control the relative expressing level of two different kind of Protein.

This is a small picture illustrating its basic condition:

Equations

$ \frac{d(repressor)}{dt} = \frac{K_{repressor} - Vm_{rep} \cdot DNA^n_{rep}}{Kp_{rep}+DNA^n_{rep}} - K_{bind}\cdot lac\cdot repressor$

$ \frac{d(DNA)}{dt} = -\frac{Vm_{rep}\cdot DNA^n_{rep}}{Kp_{rep}+DNA^n_{rep}}$

$ \frac{d([repressor-DNA])}{dt} = \frac{Vm_{rep}\cdot DNA^n_{rep}}{Kp_{rep}+DNA^n_{rep}}$

$ \frac{d([lac-repressor])}{dt} = K_{bind}\cdot lac\cdot repressor$

$ \frac{d(recombinase)}{dt} = K_{recom}\cdot DNA$

$ \frac{d(PromoterDir1)}{dt} = -\frac{Vm_{cata}\cdot PromoterDir1}{Km_{cata}+PromoterDir1}$

$ \frac{d(PromoterDir2)}{dt} = \frac{Vm_{cata}\cdot PromoterDir1}{Km_{cata}+PromoterDir1}$

$ \frac{d(Protein1)}{dt} = K_{forward}\cdot PromoterDir1$

$ \frac{d(Protein2)}{dt} = K_{backward}\cdot PromoterDir2$

Results

DgrA/DgrB Circuit

Equations

$\frac{d([cdiGMP-DgrA])}{dt} = \frac{(Ka \cdot [cdiGMP] \cdot DgrA - Kf \cdot [cdiGMP-DgrA] \cdot FliL - [Kd_{cdiGMP-DgrA}] \cdot [cdiGMP-DgrA])}{cell}$

$\frac{d([cdiGMP-DgrB])}{dt} = \frac{(Kb \cdot [cdiGMP] \cdot DgrB - [Kd_{cdiGMP-DgrB}] \cdot [cdiGMP-DgrB])}{cell}$

$\frac{d(FliL)}{dt} = \frac{(-Kf \cdot [cdiGMP-DgrA] \cdot FliL + Kg_{FliL})}{cell}$

$\frac{d([cdiGMP-DgrA-FliL])}{dt} = \frac{(Kf \cdot [cdiGMP-DgrA] \cdot FliL)}{cell}$

Results

HfiA/HfsJ Circuit

The binding affinity and cellular concentrations of HfiA and HfsJ are tuned such that this regulatory system is responsive to small changes rather than robust to large changes. This prediction is consistent with a highly responsive and sensitive regulatory system.

So the block is clear. Upstream promoter Expression Rate $\rightarrow$ HfiA $\rightarrow$ HfsJ $\rightarrow$ Catalyzing Rate

Equations

$\frac{d(P_{HfiA})}{dt} = -\frac{Vm_{HfiA} \cdot X^n_{HfiA}/(Kp_{HfiA}+X^n_{HfiA})}{cell}$

$\frac{d([X-P_{HfiA}])}{dt} = \frac{Vm_{HfiA} \cdot X^n_{HfiA}/(Kp_{HfiA}+X^n_{HfiA}) - Kd_{c1 }cdot [X-P_{HfiA}]}{cell}$

$\frac{d(mRNA_{HfiA})}{dt} = \frac{Ktc_{HfiA} \cdot P_{HfiA} - Kd_{mRNA_{HfiA}} \cdot mRNA_{HfiA} }{cell}$

$\frac{d(HfiA)}{dt} = \frac{Ktl_{HfiA} \cdot mRNA_{HfiA} - Kd_{HfiA} \cdot HfiA - Kc \cdot HfiA \cdot HfsJ}{cell}$

$\frac{d(P_{HfsJ})}{dt} = \frac{-Vm_{HfsJ} \cdot P_{HfsJ}^n{_{HfsJ}}/(Kp_{HfsJ}+P_{HfsJ}^{n_{HfsJ}})}{cell}$

$\frac{d([CtrA-P_{HfsJ}])}{dt} = \frac{Vm_{HfsJ} \cdot P_{HfsJ}^n{_{HfsJ}}/(Kp_{HfsJ}+P_{HfsJ}^{n_{HfsJ}}) - Kd_{c2 }cdot [CtrA-P_{HfsJ}]}{cell}$

$\frac{d(mRNA_{HfsJ})}{dt} = \frac{Ktc_{HfsJ} \cdot [CtrA-P_{HfsJ}] - Kd_{mRNA_{HfsJ}} \cdot mRNA_{HfsJ}}{cell}$

$\frac{d(HfsJ)}{dt} = \frac{Ktl_{HfsJ} \cdot mRNA_{HfsJ} - Kd_{HfsJ} \cdot HfsJ - Kc \cdot HfiA \cdot HfsJ - Kf \cdot HfsJ \cdot Sub-Kr \cdot EzCom + Kcat \cdot EzCom}{cell}$

$\frac{d([HfiA-HfsJ])}{dt} = \frac{Kc \cdot HfiA \cdot HfsJ}{cell}$

$\frac{d(Sub)}{dt} = \frac{-Kf \cdot HfsJ \cdot Sub-Kr \cdot EzCom}{cell}$

$\frac{d(Pro)}{dt} = \frac{Kcat \cdot EzCom}{cell}$

$\frac{d(EzCom)}{dt} = \frac{Kf \cdot HfsJ \cdot Sub-Kr \cdot EzCom - Kcat \cdot EzCom}{cell}$

Results

Others

Link to our project files