Team:USTC-Software/collaboration

From 2011.igem.org

(Difference between revisions)
 
(9 intermediate revisions not shown)
Line 1: Line 1:
__NOTOC__
__NOTOC__
-
{{:Team:USTC-Software/temp}}
+
{{:Team:USTC-Software/header}}
-
 
+
<html xmlns="http://www.w3.org/1999/xhtml">
<html xmlns="http://www.w3.org/1999/xhtml">
-
 
-
<head>
 
-
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
 
-
<title>Team:USTC-Software - 2011.igem.org/collaboration</title>
 
-
<meta name="description" content="collabration"/>
 
-
<meta name="keywords" content="USTC-Software,iGEM,synthetic biology"/>
 
-
<meta name="author" content="Ao Xu">
 
-
<meta name="copyright" content="Copyright (c) 2011 iGEM USTC Software" />
 
-
 
-
<style type="text/css">
 
-
 
-
body {
 
-
margin: 0;
 
-
padding: 0;
 
-
line-height: 1.5em;
 
-
font-family: Georgia, "Times New Roman", Times, serif;
 
-
font-size: 12px;
 
-
color: #555a4a;
 
-
    background: #a0c73f;
 
-
}
 
-
 
-
a:link, a:visited { color: #0683ab; text-decoration: none; font-weight: normal; }
 
-
a:active, a:hover { color: #0683ab; text-decoration: underline; }
 
-
 
-
img { margin: 0px; padding: 0px; border: none; }
 
-
 
-
#header_wrapper {
 
-
width: 100%;
 
-
height: 265px;
 
-
background: url(https://static.igem.org/mediawiki/2011/0/04/USTC_Software_header.jpg) no-repeat top center;
 
-
}
 
-
 
-
#header {
 
-
position: relative;
 
-
width: 980px;
 
-
height: 265px;
 
-
padding: 0 10px;
 
-
margin: 0 auto;
 
-
background: url(https://static.igem.org/mediawiki/2011/0/04/USTC_Software_header.jpg) no-repeat top center;
 
-
}
 
-
 
-
/* site title */
 
-
 
-
#header  #site_title {
 
-
float: left;
 
-
width: 180px;
 
-
padding: 33px 0;
 
-
}
 
-
 
-
#site_title h1 {
 
-
margin: 0;
 
-
padding: 3px;
 
-
font-family:"Times New Roman", Times, serif;
 
-
color:#000;
 
-
font-size:24px;
 
-
}
 
-
 
-
#site_title h1 a {
 
-
margin: 0px;
 
-
padding: 0px;
 
-
font-size: 30px;
 
-
color: #ffffff;
 
-
font-weight: bold;
 
-
text-decoration: none;
 
-
}
 
-
 
-
#site_title h1 a:hover {
 
-
font-weight: bold;
 
-
text-decoration: none;
 
-
}
 
-
 
-
#site_title h1 a span {
 
-
display: block;
 
-
margin-left: 2px;
 
-
font-size: 14px;
 
-
color: #416902;
 
-
font-weight: bold;
 
-
letter-spacing: 2px;
 
-
}
 
-
 
-
/* end of site title */
 
-
 
-
/* ---------- Navigation ---------- */
 
-
 
-
ul  {
 
-
      list-style-image: none;
 
-
}
 
-
 
-
#nav {
 
-
float: right;
 
-
font-family: Arial, Helvetica, sans-serif;
 
-
font-weight: bold;
 
-
font-size: 12px;
 
-
padding-top: 70px;
 
-
height: 43px;
 
-
}
 
-
 
-
#nav .last {
 
-
padding-right: 2px;
 
-
background-image:url(https://static.igem.org/mediawiki/2011/5/59/USTC_Software_nav_button_sep.jpg);
 
-
background-repeat:no-repeat;
 
-
background-position:right;
 
-
}
 
-
 
-
#nav, #nav ul{
 
-
margin:0;
 
-
padding:0;
 
-
list-style-type:none;
 
-
list-style-position:outside;
 
-
position:relative;
 
-
line-height:1.5em;
 
-
}
 
-
 
-
#nav .on {
 
-
color:#fff !important;
 
-
background-image: url(https://static.igem.org/mediawiki/2011/0/00/USTC_Software_nav_button_bg.jpg) !important;
 
-
background-position: bottom left !important;
 
-
}
 
-
 
-
#nav a:link, #nav a:active, #nav a:visited{
 
-
    display:block;
 
-
    padding: 14px 35px;
 
-
    color:#a5a5a5;
 
-
    text-decoration:none;
 
-
text-transform: uppercase;
 
-
background-image: url(https://static.igem.org/mediawiki/2011/0/00/USTC_Software_nav_button_bg.jpg);
 
-
background-position: top left;
 
-
}
 
-
 
-
#nav a:hover {
 
-
    color:#fff;
 
-
background-image: url(https://static.igem.org/mediawiki/2011/0/00/USTC_Software_nav_button_bg.jpg);
 
-
background-position: bottom left;
 
-
}
 
-
 
-
#nav li{
 
-
    float:left;
 
-
    position:relative;
 
-
 
-
}
 
-
 
-
#nav ul {
 
-
    position:absolute;
 
-
    width:13em;
 
-
    top:3.8em;
 
-
    display:none;
 
-
z-index: 1;
 
-
}
 
-
 
-
#nav li ul a {
 
-
    width:9em;
 
-
    float:left;
 
-
background-color: #2e2e2e;
 
-
background-image: none !important;
 
-
}
 
-
 
-
#nav li ul a:hover {
 
-
background-color: #292929;
 
-
}
 
-
 
-
#nav ul ul{
 
-
top:auto;
 
-
}
 
-
 
-
#nav li ul ul {
 
-
    left:12em;
 
-
    margin:0px 0 0 10px;
 
-
    }
 
-
 
-
#nav li:hover ul ul, #nav li:hover ul ul ul, #nav li:hover ul ul ul ul{
 
-
    display:none;
 
-
    }
 
-
 
-
#nav li:hover ul, #nav li li:hover ul, #nav li li li:hover ul, #nav li li li li:hover ul{
 
-
    display:block;
 
-
    }
 
-
 
-
/*end of nav*/
 
-
 
-
/* content_wrapper */
 
-
 
-
#content {
 
-
        background-color:#a0c73f;
 
-
        width:auto;
 
-
}
 
-
 
-
#content_wrapper {
 
-
width: 920px;
 
-
margin: 0 auto;
 
-
padding: 30px 40px;
 
-
background: url(https://static.igem.org/mediawiki/2011/b/bc/USTC_Software_content.jpg) repeat-y;
 
-
}
 
-
 
-
#content_wrapper h1 {
 
-
color: #3a5909;
 
-
font-size: 24px;
 
-
padding: 0 0 60px 20px;
 
-
background: url(https://static.igem.org/mediawiki/2011/4/4d/USTC_Software_header_bg.png) no-repeat bottom left;
 
-
}
 
-
 
-
#content_wrapper p {
 
-
font-family: Verdana, Arial;
 
-
font-size: 14px;
 
-
}
 
-
 
-
#content_wrapper_bottom {
 
-
width: 1000px;
 
-
height: 15px;
 
-
margin: 0 auto;
 
-
background: url(https://static.igem.org/mediawiki/2011/a/a5/USTC_Software_content_bottom.jpg) no-repeat;
 
-
    background-color: #a0c73f;
 
-
}
 
-
 
-
/* end of content_wrapper */
 
-
 
-
/* footer */
 
-
 
-
#mfooter {
 
-
width: 980px;
 
-
margin: 0 auto;
 
-
padding: 30px 10px;
 
-
text-align: center;
 
-
color: #333333;
 
-
    background-color: #a0c73f;
 
-
}
 
-
 
-
#mfooter a {
 
-
color:#000;
 
-
}
 
-
 
-
#mfooter .mfooter_menu {
 
-
margin: 0 0 10px 0;
 
-
padding: 0px;
 
-
list-style: none;
 
-
}
 
-
 
-
.mfooter_menu li {
 
-
margin: 0px;
 
-
padding: 0 20px;
 
-
display: inline;
 
-
border-right: 1px solid #000000;
 
-
}
 
-
 
-
.mfooter_menu .last_menu {
 
-
border: none;
 
-
}
 
-
 
-
#footer-box { background-color: transparent; border:none;} /*This is important!*/
 
-
 
-
#totop_box {
 
-
position:fixed;
 
-
width: 77px;
 
-
height: 75px;
 
-
bottom: 21px;
 
-
right: 15px;
 
-
}
 
-
 
-
#totop_box form {
 
-
width: 280px;
 
-
height: 27px;
 
-
margin: 0;
 
-
padding: 25px 0 0 0;
 
-
}
 
-
 
-
#to_org_box {
 
-
position:fixed;
 
-
width: 77px;
 
-
height: 75px;
 
-
bottom: 21px;
 
-
right: 105px;
 
-
}
 
-
 
-
#to_org_box form {
 
-
width: 280px;
 
-
height: 27px;
 
-
margin: 0;
 
-
padding: 25px 0 0 0;
 
-
}
 
-
/* end of footer*/
 
-
 
-
/* attribution */
 
-
 
-
#sidebar .attribution_list {
 
-
margin: 0;
 
-
padding: 0;
 
-
list-style: none;
 
-
}
 
-
 
-
#sidebar .attribution_list li {
 
-
padding: 0;
 
-
margin: 0;
 
-
}
 
-
 
-
.attribution_list_p{
 
-
font-family: Verdana, Arial;
 
-
font-size: 16px;
 
-
}
 
-
 
-
.attribution_list li a {
 
-
display: block;
 
-
color: #201f1c;
 
-
padding: 5px 0 5px 20px;
 
-
background: url(https://static.igem.org/mediawiki/2011/6/6d/USTC_Software_list.png) center left no-repeat;
 
-
}
 
-
.attribution_list li a:hover {
 
-
color: #537c11;
 
-
text-decoration: none;
 
-
}
 
-
 
-
</style>
 
-
 
-
</head>
 
-
 
-
 
<body>
<body>
-
<div id="header_wrapper">
 
-
 
-
    <div id="header">
 
-
   
 
-
  <div id="site_title">
 
-
            <h1><img src="https://static.igem.org/mediawiki/2011/b/b8/USTC_Software_iGEM_logo.jpg" alt="iGEM" width="173" height="150" /> </br> </br> USTC-Software </h1>
 
-
        </div>
 
-
       
 
-
<!-- navigation -->
 
-
<ul id="nav">
 
-
<li><a href="https://2011.igem.org/Team:USTC-Software">Home</a></li>
 
-
 
-
   
 
-
    <li><a href="https://2011.igem.org/Team:USTC-Software/project">Project</a>
 
-
<ul>
 
-
<li><a href="https://2011.igem.org/Team:USTC-Software/documents">Documents Parser & Models </a></li>
 
-
<li><a href="#">Views</a>
 
-
                                  <ul>
 
-
                                      <li><a href="https://2011.igem.org/Team:USTC-Software/aview">Assembly View</a></li>
 
-
                                      <li><a href="https://2011.igem.org/Team:USTC-Software/bview">Behavior View</a></li>
 
-
                                      <li><a href="https://2011.igem.org/Team:USTC-Software/nview">Network View</a></li>                     
 
-
                                  </ul>
 
-
                                </li>
 
-
 
-
                                <li><a href="#">Technology & Algorithm</a>
 
-
                                  <ul>
 
-
                                    <li><a href="https://2011.igem.org/Team:USTC-Software/parameter">parameter analysis</a></li>
 
-
                                    <li><a href="https://2011.igem.org/Team:USTC-Software/model">MoDeL</a></li>
 
-
                                  </ul>
 
-
                                </li>
 
-
                                   
 
-
                                <li><a href="https://2011.igem.org/Team:USTC-Software/tutorial">Tutorial & Demo</a></li>
 
-
</ul>
 
-
</li>
 
-
 
-
 
-
                <li><a href="https://2011.igem.org/Team:USTC-Software/notebook">Notebook</a></li>
 
-
 
-
<li><a class="on" href="https://2011.igem.org/Team:USTC-Software/team">Team</a>
 
-
                          <ul>
 
-
                            <li><a href="https://2011.igem.org/Team:USTC-Software/members">members</a></li>
 
-
                            <li><a href="https://2011.igem.org/Team:USTC-Software/collaboration">collaboration</a></li>
 
-
                            <li><a href="https://2011.igem.org/Team:USTC-Software/attribution">attribution & contributions</a></li>
 
-
                            <li><a href="https://2011.igem.org/Team:USTC-Software/acknowledgements">acknowledgements</a></li>
 
-
                          </ul>
 
-
                </li>
 
-
 
 
-
  <li><a href="https://2011.igem.org/Team:USTC-Software/human practice">Human practice</a>             
 
-
  <ul>
 
-
    <li><a href="https://2011.igem.org/Team:USTC-Software/safety">safety</a></li>
 
-
    <li><a href="https://2011.igem.org/Team:USTC-Software/meetup">meetup</a></li>
 
-
  </ul>
 
-
</li> 
 
-
 
-
<li><a class="last" href="https://2011.igem.org/Team:USTC-Software/download">Download</a></li>
 
-
           
 
-
</ul>
 
-
<!-- // end of navigation -->
 
-
     
 
-
      <div class="cleaner"></div>
 
-
    </div> <!-- end of header -->
 
-
 
-
</div>
 
-
<!-- end of header_wrapper -->
 
-
 
<div id="content_wrapper">
<div id="content_wrapper">
<h1>Collaboration</h1>
<h1>Collaboration</h1>
-
<p>We had a great time during the summer with our wet team and the Ecoli they raised.</p>
+
<p>We spent a happy summer working with our wet team, USTC-China, through intimate collaboration on how our software tool can help and guide their experimental designs before they are actually constructed in the tube. It cannot be denied that most of our members are lack of biological background but do not worry, they are there available to discuss with us to help us understand their proposals. Specially, we are going to thank Junhui Peng, the modeler of USTC-China, who contributed very much to our cooperation.</p>
-
<p>They helped us a lot with biology background, and we offer necessray modeling and technical support for them. The project they designed is interesting and creative.</p>
+
<p>The main idea of their project is to use quorum sensing mechanism to guide bacteria towards the concentration gradient. The toggle switch plays a central role to run the device by switching the expression of chemotaxis protein cheZ on and off iteratively. When the density of the cell population has reached a certain criteria, the accumulation of AHL molecules exceeds the threshold and thus trigger the toggle switch to change its state from the other one(auotomatic toggle). More information are available <a href="https://2011.igem.org/Team:USTC-China/Project"> here</a>.</p>
-
<p>We discuss a lot with them, expecially with their modeler, Junhui Peng.</p>
+
<p>To test our software tool, we planned to use our rule-based modeling approach to get some insights about their work. It can be divided further to two parts: (1) automatically generate the biochemical network of the quorum sensing subsystem combined with toggle switch as a central part; (2) analyze the phase space of the toggle switch.</p>
-
<p>The main ideal of their project is to use quorum sensing to guide bacteria towards the concentration gradient.The toggle swtich is an important device in their design, which can turn on and off the chemotaxis protein cheZ expression. When the density of the cell population reached certain criteria, there are enough AHL molecules to toggle the swtich to another state(auotomatic toggle).</p>
+
<img src="https://static.igem.org/mediawiki/2011/f/f4/USTC_Software_b1.png" width="750" />
-
<p>For more information about their project, see<a href="https://2011.igem.org/Team:USTC-China/Project"> here</a>.</p>
+
<p>The design-b is our target system to be modeled. In terms of simplicity, we try to ignore their quorum sensing part by setting of parameter of the lasi catalyzing ability low, and add IPTG and AHL deliberately at certain timepoints.</p>
-
 
+
-
<p>We did a lot of collaboration, including using our rulebased modelling approach to automatically generate the network for their design and analyze the phase space of the toggle switch.</p>
+
-
 
+
-
figure1
+
-
 
+
-
<p>We choose their design-b to test our software. In terms of simplicity, we try to ignore their quorum sensing part by setting of parameter of the lasi catalyzing ability low, and add IPTG and AHL deliberately at certain timepoints.</p>
+
<p>From such a mere assembly, our software can generate the network associated with the assembly.</p>
<p>From such a mere assembly, our software can generate the network associated with the assembly.</p>
Line 407: Line 21:
<p>The network <strong>generated</strong> and <strong>visualized</strong> by our software's  rule based modeling part and visualization part is shown below:</p>
<p>The network <strong>generated</strong> and <strong>visualized</strong> by our software's  rule based modeling part and visualization part is shown below:</p>
-
figure2
+
<img src="https://static.igem.org/mediawiki/2011/9/95/USTC_Software_b2.png" width="750" />
<p>There are so many species generated by our approch(however, this is an advantage of our approach, in the future version of our software, we plan to offer visualizing tool to show the structure of the species, ie, the binding state dna chain)</p>
<p>There are so many species generated by our approch(however, this is an advantage of our approach, in the future version of our software, we plan to offer visualizing tool to show the structure of the species, ie, the binding state dna chain)</p>
Line 419: Line 33:
<p>All the figures below are in one simulation.</p>
<p>All the figures below are in one simulation.</p>
-
figure3-9
+
<img src="https://static.igem.org/mediawiki/2011/f/f6/USTC_Software_a1.png" width="750" />
 +
<img src="https://static.igem.org/mediawiki/2011/6/69/USTC_Software_a2.png" width="750" />
 +
<img src="https://static.igem.org/mediawiki/2011/a/ae/USTC_Software_a3.png" width="750" />
 +
<img src="https://static.igem.org/mediawiki/2011/2/21/USTC_Software_a4.png" width="750" />
 +
<img src="https://static.igem.org/mediawiki/2011/2/2e/USTC_Software_a5.png" width="750" />
 +
<img src="https://static.igem.org/mediawiki/2011/4/4f/USTC_Software_a6.png" width="750" />
 +
<img src="https://static.igem.org/mediawiki/2011/c/c0/USTC_Software_a7.png" width="750" />
 +
 
 +
<p>The parameters we tried are as follows:</p>
 +
</html>
 +
 
 +
    1  NA            6.020e+23 
 +
    2  f              1.000e+00 
 +
    3  Ve            4.000e-01  # f*0.4
 +
    4  V              7.000e-16  # f*7e-16
 +
    5  N              1.000e+00 
 +
    6  theo_init      0.000e+00 
 +
    7  ahl_init      0.000e+00 
 +
    8  dna1_init      2.373e-08  # 10/NA/V
 +
    9 dna2_init      2.373e-07  # 100/NA/V
 +
  10  pulse1_start  2.000e+04 
 +
  11  pulse2_start  4.000e+04 
 +
  12  pulse1_conc    1.000e-04 
 +
  13  pulse2_conc    1.000e-04 
 +
  14  transp3_pout  1.000e-01 
 +
  15  transp3_pin    1.000e-01 
 +
  16  transp4_pout  1.000e-01 
 +
  17  transp4_pin    1.000e-01 
 +
  18  rule20_k      5.000e-04 
 +
  19  rule15_k      5.000e-01 
 +
  20  rule17_k      5.000e-01 
 +
  21  rule34_k      1.000e-04 
 +
  22  rule35_k      1.000e-04 
 +
  23  rule24_k      1.155e-02 
 +
  24  rule26_k      5.783e-03 
 +
  25  rule21_k      1.000e+08 
 +
  26  rule1_k        1.000e+07 
 +
  27  rule33_k      2.310e-03 
 +
  28  rule27_k      2.310e-03 
 +
  29  rule22_k      2.000e+01 
 +
  30  rule25_k      1.155e-02 
 +
  31  rule5_k        1.000e+07 
 +
  32  rule32_k      2.310e-03 
 +
  33  rule28_k      2.310e-03 
 +
  34  rule23_k      1.000e-08 
 +
  35  rule29_k      2.310e-03 
 +
  36  rule9_k        1.540e+05 
 +
  37  rule11_k      1.250e+07 
 +
  38  rule30_k      2.310e-03 
 +
  39  rule2_k        1.000e+01 
 +
  40  rule3_k        1.000e+08 
 +
  41  rule31_k      2.310e-03 
 +
  42  rule6_k        1.000e+01 
 +
  43  rule7_k        2.000e+08 
 +
  44  rule10_k      2.000e-01 
 +
  45  rule12_k      1.000e+01 
 +
  46  rule4_k        1.000e-02 
 +
  47  rule18_k      5.000e-04 
 +
  48  rule8_k        1.000e-02 
 +
  49  rule16_k      5.000e-04 
 +
  50  rule13_k      2.000e+10 
 +
  51  rule14_k      4.000e-02 
 +
  52  rule19_k      5.000e-01
 +
 
 +
Files: [[File:amitosis.zip]]
 +
Mysql input file: amitosis_2.sql function.sql inducer.sql medium.sql
 +
MoDeL input file: amitosis_2.model
 +
MoDel output file: amitosis_2.net
 +
SBML output file: amitosis_2.xml
 +
 
 +
----
 +
 
 +
<html>
<p><strong>OTHER collaborations(not all): bifurcation analysis of the toggle switch</strong></p>
<p><strong>OTHER collaborations(not all): bifurcation analysis of the toggle switch</strong></p>
Line 433: Line 119:
<p>We first worked out the nullclines of the toggle switch system, with a symmetric toggle switch equation, we get a figure as below.</p>
<p>We first worked out the nullclines of the toggle switch system, with a symmetric toggle switch equation, we get a figure as below.</p>
-
figure 10
+
<img src="https://static.igem.org/mediawiki/2011/thumb/0/08/USTC_Software_b3.png/800px-USTC_Software_b3.png" width="750" />
<p>There are three cross sections ,or three steady state, two of them are stable steady state and one of them unstable, as indicated by the vector fields.</p>
<p>There are three cross sections ,or three steady state, two of them are stable steady state and one of them unstable, as indicated by the vector fields.</p>
Line 445: Line 131:
<p>Below is a figure of critical parameter values, the bifurcation point.At this critical point, there are only two cross-sections , ie two equlibrium state.</p>
<p>Below is a figure of critical parameter values, the bifurcation point.At this critical point, there are only two cross-sections , ie two equlibrium state.</p>
-
figure11
+
<img src="https://static.igem.org/mediawiki/2011/b/be/USTC_Software_b4.png" width="750" />
<p>If we variate the parameter more, as illustrated as below, there is exactly one equlibrium state, indicating that no matter what the initial condition is, the system will be draged to it’s mono-steady state of one side on and the other side of the toggle switch down.</p>
<p>If we variate the parameter more, as illustrated as below, there is exactly one equlibrium state, indicating that no matter what the initial condition is, the system will be draged to it’s mono-steady state of one side on and the other side of the toggle switch down.</p>
-
figure12
+
<img src="https://static.igem.org/mediawiki/2011/f/fb/USTC_Software_b5.png" width="750" />
<p>So the problem may be the parameter of the promoter strength.</p>
<p>So the problem may be the parameter of the promoter strength.</p>
Line 459: Line 145:
<p><strong>We are going to collaborate more in the future.</strong></p>
<p><strong>We are going to collaborate more in the future.</strong></p>
 +
 +
<br/><br/>
 +
</div>
</div>

Latest revision as of 06:59, 23 October 2011


Team:USTC-Software - 2011.igem.org

Collaboration

We spent a happy summer working with our wet team, USTC-China, through intimate collaboration on how our software tool can help and guide their experimental designs before they are actually constructed in the tube. It cannot be denied that most of our members are lack of biological background but do not worry, they are there available to discuss with us to help us understand their proposals. Specially, we are going to thank Junhui Peng, the modeler of USTC-China, who contributed very much to our cooperation.

The main idea of their project is to use quorum sensing mechanism to guide bacteria towards the concentration gradient. The toggle switch plays a central role to run the device by switching the expression of chemotaxis protein cheZ on and off iteratively. When the density of the cell population has reached a certain criteria, the accumulation of AHL molecules exceeds the threshold and thus trigger the toggle switch to change its state from the other one(auotomatic toggle). More information are available here.

To test our software tool, we planned to use our rule-based modeling approach to get some insights about their work. It can be divided further to two parts: (1) automatically generate the biochemical network of the quorum sensing subsystem combined with toggle switch as a central part; (2) analyze the phase space of the toggle switch.

The design-b is our target system to be modeled. In terms of simplicity, we try to ignore their quorum sensing part by setting of parameter of the lasi catalyzing ability low, and add IPTG and AHL deliberately at certain timepoints.

From such a mere assembly, our software can generate the network associated with the assembly.

The network generated and visualized by our software's rule based modeling part and visualization part is shown below:

There are so many species generated by our approch(however, this is an advantage of our approach, in the future version of our software, we plan to offer visualizing tool to show the structure of the species, ie, the binding state dna chain)

In order to test the validity of our softare, we tried the time course.

We select different interesting species eachtime to see the time course. We worked hard on adjusting the parameters of the system to acquire a good behaviour. And the result is quite satisfactory.

At time point 20000s, the theophline is added and at timepoint 40000s AHL both with e-4 moles or 100micro moles.

All the figures below are in one simulation.

The parameters we tried are as follows:

   1  NA             6.020e+23  
   2  f              1.000e+00  
   3  Ve             4.000e-01   # f*0.4
   4  V              7.000e-16   # f*7e-16
   5  N              1.000e+00  
   6  theo_init      0.000e+00  
   7  ahl_init       0.000e+00  
   8  dna1_init      2.373e-08   # 10/NA/V
   9  dna2_init      2.373e-07   # 100/NA/V
  10  pulse1_start   2.000e+04  
  11  pulse2_start   4.000e+04  
  12  pulse1_conc    1.000e-04  
  13  pulse2_conc    1.000e-04  
  14  transp3_pout   1.000e-01  
  15  transp3_pin    1.000e-01  
  16  transp4_pout   1.000e-01  
  17  transp4_pin    1.000e-01  
  18  rule20_k       5.000e-04  
  19  rule15_k       5.000e-01  
  20  rule17_k       5.000e-01  
  21  rule34_k       1.000e-04  
  22  rule35_k       1.000e-04  
  23  rule24_k       1.155e-02  
  24  rule26_k       5.783e-03  
  25  rule21_k       1.000e+08  
  26  rule1_k        1.000e+07  
  27  rule33_k       2.310e-03  
  28  rule27_k       2.310e-03  
  29  rule22_k       2.000e+01  
  30  rule25_k       1.155e-02  
  31  rule5_k        1.000e+07  
  32  rule32_k       2.310e-03  
  33  rule28_k       2.310e-03  
  34  rule23_k       1.000e-08  
  35  rule29_k       2.310e-03  
  36  rule9_k        1.540e+05  
  37  rule11_k       1.250e+07  
  38  rule30_k       2.310e-03  
  39  rule2_k        1.000e+01  
  40  rule3_k        1.000e+08  
  41  rule31_k       2.310e-03  
  42  rule6_k        1.000e+01  
  43  rule7_k        2.000e+08  
  44  rule10_k       2.000e-01  
  45  rule12_k       1.000e+01  
  46  rule4_k        1.000e-02  
  47  rule18_k       5.000e-04  
  48  rule8_k        1.000e-02  
  49  rule16_k       5.000e-04  
  50  rule13_k       2.000e+10  
  51  rule14_k       4.000e-02  
  52  rule19_k       5.000e-01 

Files: File:Amitosis.zip

Mysql input file: amitosis_2.sql function.sql inducer.sql medium.sql
MoDeL input file: amitosis_2.model
MoDel output file: amitosis_2.net
SBML output file: amitosis_2.xml

OTHER collaborations(not all): bifurcation analysis of the toggle switch

In the simulation of the design-a (see figure above), their modeler Junhui Peng told us that he was confused by a phenomenon: if the promoter P-lasbox strength excess a critical value, no matter how much AHL molecules are produced by lasI (in silicon), the toggle switch just stay in one state and can’t leap to the other.

We conducted a survey to this phenomenon , trying to figure out the mechanism of it.

But the network of the their origial design is a little bit complex to analyze , so we first tried a much simplied model of the toggle switch. We believe that a good work start from the simple thing.

With the help of Professor Haiyan Liu, our team’s mathematical student Hui Li developed an algorithm to analyse the effects of parameters to the phasespace of the toggle swtich. Later, Hui Li speeded up her algorithm.

We first worked out the nullclines of the toggle switch system, with a symmetric toggle switch equation, we get a figure as below.

There are three cross sections ,or three steady state, two of them are stable steady state and one of them unstable, as indicated by the vector fields.

From the vector field we can interpret an important mechanism, after the inducer is added to the system for long enough time, even the induder is removed from the system, the system still stay in the induced state and wouldn’t come back to the orginal one . See the figure above, The state space of this system is divided into two regions, each with an “attracting center”, once the system enters an arbitray state in one region, it will defininitely be “attracted” to the stable steady state of that region.

The task of the inducer is to change the state point of the system, moving it from the orginal region to the oppisite region.

Then we devise a way to move the two nullclines by variating the parameters, trying to separate the two nullclines with less cross sections.

Below is a figure of critical parameter values, the bifurcation point.At this critical point, there are only two cross-sections , ie two equlibrium state.

If we variate the parameter more, as illustrated as below, there is exactly one equlibrium state, indicating that no matter what the initial condition is, the system will be draged to it’s mono-steady state of one side on and the other side of the toggle switch down.

So the problem may be the parameter of the promoter strength.

When Junhui Peng modified the promoter strength to a lower value and simulate, he found that this time with enough amount of the AHL inducer, the system is able to switch between the two steady states.

During the spare time, we hang out with our wet friends and talked a lot.

We are going to collaborate more in the future.