Team:USTC-Software/collaboration

From 2011.igem.org

(Difference between revisions)
 
(12 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">
 +
<body>
 +
<div id="content_wrapper">
-
<head>
+
<h1>Collaboration</h1>
-
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
+
-
<title>Team:USTC-Software - 2011.igem.org/attribution</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">
+
<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>
-
body {
+
<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>
-
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; }
+
<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>
-
a:active, a:hover { color: #0683ab; text-decoration: underline; }
+
-
img { margin: 0px; padding: 0px; border: none; }
+
<img src="https://static.igem.org/mediawiki/2011/f/f4/USTC_Software_b1.png" width="750" />
-
#header_wrapper {
+
<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>
-
width: 100%;
+
-
height: 265px;
+
-
background: url(https://static.igem.org/mediawiki/2011/0/04/USTC_Software_header.jpg) no-repeat top center;
+
-
}
+
-
#header {
+
<p>From such a mere assembly, our software can generate the network associated with the assembly.</p>
-
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 */
+
<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>
-
#header  #site_title {
+
<img src="https://static.igem.org/mediawiki/2011/9/95/USTC_Software_b2.png" width="750" />
-
float: left;
+
-
width: 180px;
+
-
padding: 33px 0;
+
-
}
+
-
#site_title h1 {
+
<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>
-
margin: 0;
+
-
padding: 3px;
+
-
font-family:"Times New Roman", Times, serif;
+
-
color:#000;
+
-
font-size:24px;
+
-
}
+
-
#site_title h1 a {
+
<p>In order to test the validity of our softare, we tried the time course.</p>
-
margin: 0px;
+
-
padding: 0px;
+
-
font-size: 30px;
+
-
color: #ffffff;
+
-
font-weight: bold;
+
-
text-decoration: none;
+
-
}
+
-
#site_title h1 a:hover {
+
<p>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.</p>
-
font-weight: bold;
+
-
text-decoration: none;
+
-
}
+
-
#site_title h1 a span {
+
<p>At time point 20000s, the theophline is added and at timepoint 40000s AHL both with e-4 moles or 100micro moles.</p>
-
display: block;
+
-
margin-left: 2px;
+
-
font-size: 14px;
+
-
color: #416902;
+
-
font-weight: bold;
+
-
letter-spacing: 2px;
+
-
}
+
-
/* end of site title */
+
<p>All the figures below are in one simulation.</p>
-
/* ---------- Navigation ---------- */
+
<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" />
-
ul  {
+
<p>The parameters we tried are as follows:</p>
-
      list-style-image: none;
+
</html>
-
}
+
-
#nav {
+
    1  NA            6.020e+23 
-
float: right;
+
    2  f              1.000e+00 
-
font-family: Arial, Helvetica, sans-serif;
+
    3  Ve            4.000e-01  # f*0.4
-
font-weight: bold;
+
    4  V              7.000e-16  # f*7e-16
-
font-size: 12px;
+
    5  N              1.000e+00 
-
padding-top: 70px;
+
    6  theo_init      0.000e+00 
-
height: 43px;
+
    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
-
#nav .last {
+
Files: [[File:amitosis.zip]]
-
padding-right: 2px;
+
Mysql input file: amitosis_2.sql function.sql inducer.sql medium.sql
-
background-image:url(https://static.igem.org/mediawiki/2011/5/59/USTC_Software_nav_button_sep.jpg);
+
MoDeL input file: amitosis_2.model
-
background-repeat:no-repeat;
+
MoDel output file: amitosis_2.net
-
background-position:right;
+
SBML output file: amitosis_2.xml
-
}
+
-
#nav, #nav ul{
+
----
-
margin:0;
+
-
padding:0;
+
-
list-style-type:none;
+
-
list-style-position:outside;
+
-
position:relative;
+
-
line-height:1.5em;
+
-
}
+
-
#nav .on {
+
<html>
-
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{
+
<p><strong>OTHER collaborations(not all): bifurcation analysis of the toggle switch</strong></p>
-
    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 {
+
<p>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.</p>
-
    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{
+
<p>We conducted a survey to this phenomenon , trying to figure out the mechanism of it.</p>
-
    float:left;
+
-
    position:relative;
+
-
+
-
}
+
-
#nav ul {
+
<p>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.</p>
-
    position:absolute;
+
-
    width:13em;
+
-
    top:3.8em;
+
-
    display:none;
+
-
z-index: 1;
+
-
}
+
-
#nav li ul a {
+
<p>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.</p>
-
    width:9em;
+
-
    float:left;
+
-
background-color: #2e2e2e;
+
-
background-image: none !important;
+
-
}
+
-
#nav li ul a:hover {
+
<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>
-
background-color: #292929;
+
-
}
+
-
#nav ul ul{
+
<img src="https://static.igem.org/mediawiki/2011/thumb/0/08/USTC_Software_b3.png/800px-USTC_Software_b3.png" width="750" />
-
top:auto;
+
-
}
+
-
#nav li ul ul {
+
<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>
-
    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{
+
<p>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.</p>
-
    display:none;
+
-
    }
+
-
#nav li:hover ul, #nav li li:hover ul, #nav li li li:hover ul, #nav li li li li:hover ul{
+
<p>The task of the inducer is to change the state point of the system, moving it from the orginal region to the oppisite region.</p>
-
    display:block;
+
-
    }
+
-
/*end of nav*/
+
<p>Then we devise a way to move the two nullclines by variating the parameters, trying to separate the two nullclines with less cross sections.</p>
-
/* content_wrapper */
+
<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>
-
+
-
#content {
+
-
        background-color:#a0c73f;
+
-
        width:auto;
+
-
}
+
-
#content_wrapper {
+
<img src="https://static.igem.org/mediawiki/2011/b/be/USTC_Software_b4.png" width="750" />
-
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 {
+
<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>
-
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 {
+
<img src="https://static.igem.org/mediawiki/2011/f/fb/USTC_Software_b5.png" width="750" />
-
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 */
+
<p>So the problem may be the parameter of the promoter strength.</p>
-
/* footer */
+
<p>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.</p>
-
#mfooter {
+
<p>During the spare time, we hang out with our wet friends and talked a lot.</p>
-
width: 980px;
+
-
margin: 0 auto;
+
-
padding: 30px 10px;
+
-
text-align: center;
+
-
color: #333333;
+
-
    background-color: #a0c73f;
+
-
}
+
-
#mfooter a {
+
<p><strong>We are going to collaborate more in the future.</strong></p>
-
color:#000;
+
-
}
+
-
#mfooter .mfooter_menu {
+
<br/><br/>
-
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>
 
-
<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">
 
-
 
-
<h1>Collaboration</h1>
 
</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.