Team:NCTU Formosa/modeling

From 2011.igem.org

(Difference between revisions)
 
(56 intermediate revisions not shown)
Line 28: Line 28:
#globalWrapper{
#globalWrapper{
-
background-color: #565659;
+
//background-color: #565659;
 +
background-image:url(https://static.igem.org/mediawiki/2011/a/ae/Wallpaper2.jpg);
}
}
Line 179: Line 180:
}
}
 +
<style type = "text/javascript" language = "JavaScript">
 +
function out(id_num){
 +
  var oa = document.getElementById(id_num);
 +
  if(oa.style.display == "none"){
 +
    oa.style.display = "block"; 
 +
  }
 +
  return false;
 +
}
 +
</script>
/*content*/
/*content*/
Line 287: Line 297:
   font-weight: bold;
   font-weight: bold;
}
}
 +
 +
</style>
</style>
 +
<style>
 +
 +
table {
 +
  border: 0;
 +
  font-family: Calibri, Verdana, helvetica, sans-serif;
 +
  font-size:18px;
 +
  background-color: #ede8e2;
 +
  margin: 5px 0 5px 0;
 +
  padding: 5px 0 5px 10px;
 +
  text-align:top;
 +
}
 +
td{vertical-align:top;}
 +
 +
h3 {              font-size: 20px;
 +
                  text-align:right;
 +
                  text-shadow: #C0C0C0  0px 0px 2px;
 +
                  vertical-align: baseline;
 +
            }
 +
div.intro img {   
 +
                        border-color: white ;
 +
border-style:solid;
 +
                        float:left;
 +
                        width:expression(this.width>500?"250px":this.width+"px");
 +
}
 +
 +
 +
 +
 +
</style>
 +
 +
</head>
</head>
Line 298: Line 341:
<ul id="cm-nav">
<ul id="cm-nav">
   <li><a href="https://2011.igem.org/Team:NCTU_Formosa">Home</a></li>
   <li><a href="https://2011.igem.org/Team:NCTU_Formosa">Home</a></li>
-
   <li><a class="arrow no-click">Team </a>
+
   <li><a onClick="out('cm-nav')" class="arrow">Team </a>
       <ul>
       <ul>
         <li><a href="https://2011.igem.org/Team:NCTU_Formosa/members">Members</a></li>
         <li><a href="https://2011.igem.org/Team:NCTU_Formosa/members">Members</a></li>
Line 305: Line 348:
       </ul>
       </ul>
   </li>
   </li>
-
   <li><a class="arrow no-click">Project</a>
+
   <li><a onClick="out('cm-nav')" class="arrow">Project</a>
       <ul class="arrow-pad">
       <ul class="arrow-pad">
         <li><a href="https://2011.igem.org/Team:NCTU_Formosa/introduction">Introduction</a></li>
         <li><a href="https://2011.igem.org/Team:NCTU_Formosa/introduction">Introduction</a></li>
-
         <li><a class="arrow no-click">RNA Thermometer</a>
+
         <li><a onClick="out('cm-nav')" class="arrow">RNA Thermometer</a>
             <ul>
             <ul>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/RNA_design">Design</a></li>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/RNA_design">Design</a></li>
Line 315: Line 358:
             </ul>
             </ul>
         </li>
         </li>
-
         <li><a class="arrow no-click">CI promoter </a>
+
         <li><a onClick="out('cm-nav')" class="arrow">CI promoter </a>
             <ul>
             <ul>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/CI_design">Design</a></li>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/CI_design">Design</a></li>
Line 322: Line 365:
             </ul>
             </ul>
         </li>
         </li>
-
         <li><a class="arrow no-click">Carotenoid synthesis pathway</a>
+
         <li><a onClick="out('cm-nav')" class="arrow">Carotenoid synthesis pathway</a>
             <ul>
             <ul>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/CSP_design">Design</a></li>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/CSP_design">Design</a></li>
Line 328: Line 371:
             </ul>
             </ul>
         </li>
         </li>
-
         <li><a class="arrow no-click">Butanol pathway</a>
+
         <li><a onClick="out('cm-nav')" class="arrow">Butanol pathway</a>
             <ul>
             <ul>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/BP_design">Design</a></li>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/BP_design">Design</a></li>
Line 334: Line 377:
             </ul>
             </ul>
         </li>
         </li>
-
         <li><a class="arrow no-click">Violacein pathway</a>
+
         <li><a onClick="out('cm-nav')" class="arrow">Violacein pathway</a>
             <ul>
             <ul>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/VP_design">Design</a></li>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/VP_design">Design</a></li>
Line 347: Line 390:
   <li><a href="https://2011.igem.org/Team:NCTU_Formosa/humanpractice">Human Practice</a></li>
   <li><a href="https://2011.igem.org/Team:NCTU_Formosa/humanpractice">Human Practice</a></li>
   <li><a href="https://2011.igem.org/Team:NCTU_Formosa/contributions">Attribution</a></li>
   <li><a href="https://2011.igem.org/Team:NCTU_Formosa/contributions">Attribution</a></li>
-
   <li><a class="arrow no-click">Notebook </a>
+
   <li><a onClick="out('cm-nav')" class="arrow">Notebook </a>
       <ul>
       <ul>
-
           <li><a class="arrow no-click">Protocols</a>
+
           <li><a onClick="out('cm-nav')" class="arrow">Protocols</a>
             <ul>
             <ul>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/protocol">Mutation</a></li>
               <li><a href="https://2011.igem.org/Team:NCTU_Formosa/protocol">Mutation</a></li>
Line 362: Line 405:
<br><br>
<br><br>
<div id="blueBox"><p>Measurement</p></div>
<div id="blueBox"><p>Measurement</p></div>
-
<div id="Box"><h2> Modeling and simulations of high temperature induced device BBa_K098995 – cI promoter & cI repressor
+
<div id="Box"><h2> This new measurement method can calculate the protein expression rate in the different E. coli population density for the protein expression device (=Promoter_RBS device).</h2>
-
</h2>
+
 
 +
<br><b>Summary</b>
 +
<p>In previous studies, the transcriptional strengths of promoter and transcriptional strengths of RBSs was defined as a constant values. But in biological concepts, we know the expression rates of most proteins decreased dramatically while the bacteria at the stationary phase. To overcome this problem, our team developed a new measurement method can calculate the protein expression activity of promoter_RBS device changes with cell density in the culture tube directly.
 +
We provide a simple polynomial equation which can describe linear relationship between the protein expression activity p(s) of promoter_RBS device and cell desity s (the value of OD600).<br>
 +
<p>
 +
p(s) = p<sub>0</sub>+p<sub>1</sub>s
 +
<P>
 +
<br>where p0 denotes zero-order coefficient, p1 denotes first-order coefficient.
 +
This dynamic model captured the experimentally observed differences of green fluorescence intensity for each protein expression circuit. As we expect, promoter-RBS activity decrease linearly related with cell population density. Therefore, if the P0 and P1 parameters of a promoter_RBS device are characterized, we can calculate the protein expression activity of a promoter_RBS device by OD600 value.
 +
Furthermore, this model enable us to rational connect a promote_RBS device with different strength to obtain a target protein expression level in a synthetic genetic circuit.<br></p>
 +
<p><br><b>Motivation</b>
 +
<br>In previous studies to model a protein expression, the transcription rates of promoters and translation rates of ribosome binding sites (RBSs) are defined as a constant value. However, the values are not constant during cell growth. For example, the transcription rates of promoters and the translation rates of RBSs are lager in log phase than bacteria growth in stationary phase. To overcome this problem, we selected four promoters and three RBSs with different regulation strength and constructed 12 protein expression devices which combine promoter, RBS and green fluorescent protein (GFP) in Escherichia coli. The GFP expression levels with time were measured using a flow cytometry, and the experimental data can used to characterize a protein expression rate of a protein expression device which contains a promoters and a RBS. A dynamic model that captured the experimentally observed differences for each protein expression device was developed in this study. Using this method, we can measurement the protein expression rate in the different E. coli population density for the protein expression device.<br></p>
 +
<p><br><b>Methods</b>
 +
<br>Two strategies will be applied in this method:
 +
<br><b>1.Construction of a Promoter-RBS Library and Assay of GFP Expression</b>
 +
<br>We selected four promoters and three RBSs with different regulation strength and constructed 12 protein expression devices which combine promoter, RBS and green fluorescent protein (GFP) in E. coli (Fig.1 &table.1). The GFP expression level with time was measured using a flow cytometry (Fig.3 &Fig.4)<br></p>
 +
<br><div><img src = "https://static.igem.org/mediawiki/2011/b/ba/P-1.jpg" width="450"></div><br>
 +
<b>Figure 1: Combinatorial promoter-RBS architecture reveals constitutive expression.</b>
 +
<br>The reporter gene, green fluorescence protein (GFP) locates downstream of promoter-RBS devices as output. The output expression is controlled by differently combinatorial promoter-RBS devices. Each promoter-RBS device is expected to show the significant expression diversity. The combinatorial library contains 12 sets of promoter-RBS devices. Each set consists of different promoters and RBSs.<br>
 +
<br><div><img src = "https://static.igem.org/mediawiki/2011/0/0c/T-1.jpg" width="450"></div><br>
 +
<b>Table 1:</b> The name of 12 green fluorescence protein devices and their composed biobricks
 +
<br>A (*) refers to construct in two different backbone.<br>
 +
<br><div><img src = "https://static.igem.org/mediawiki/2011/6/65/P-2.jpg" width="450"></div><br>
 +
<b>Figure 2: The green fliorescence intensity changed with time provided the useful hints of the interactiofluorescence intensity per cell was measured using a flow cytometer.</b>
 +
<br>The O.D.ratio and the fluorescence data n function between promoters and RBSs.<br>
 +
<br><div><img src = "https://static.igem.org/mediawiki/2011/f/fa/P-3.png" width="450"></div>
 +
 
 +
<br><b>Figure 3: The time-course expression results were measured with time for 12 sets of GFP expression devices.</b> <br>(a) promoter-RBS devices which contained the same promoter J23105 combined with three different RBSs (B0030, B0032, and B0034) were selected for measuring. (b) promoter-RBS devices which contained the same promoter J23106 combined with three different RBSs (B0030, B0032, and B0034) were selected for measuring. (c) promoter-RBS devices which contain the same promoter J23114 combined with three different RBSs (B0030, B0032, and B0034) were selected for measuring. (d) promoter-RBS devices which contained the same repressible promoter R0040 combined with three different RBSs (B0030, B0032, and B0034) were selected for measuring. Each measuring even was detected every 15 minutes. And all of the data represented the average of three independent measurements. Error bars indicated standard deviations. X-axis indicated the time units, and Y-axis indicated the fluorescence units with different scales. Furthermore, the fluorescent signals changed with time per cell were measured by using a flow cytometer.<br>
 +
<p><br><b>2.Build a Mathematical Model to Characterize Protein Expression Ability of Promoter-RBS Devices.</b>
 +
<br>The logistic growth model is commonly used to describe the bacterium growth curve under the nutrient limited condition as follow [1-2]:<br></p>
 +
<div><img src = "https://static.igem.org/mediawiki/2011/b/b4/P-4.jpg" width="200"></div>  (1)
 +
<p><br>Where s denotes the cell density parameter (optical density at 600 nm (O.D. 600)), smax denotes the maximum value of cell density, and ks denotes growth rate constant. The cell density time course data of the bacteria with different GFP expression devices were used to solve Eq. (1) (Fig. 4). The growth rate constant ks was calculated and listed in Table 2.<br></p>
 +
<div><img src = "https://static.igem.org/mediawiki/2011/5/57/Fig4.jpg" width="700"></div>
 +
<br><b>Figure 4:</b>The bacterial population simulated results of 12 sets of promoter-RBS expression devices. Blue dotted lines were the value of OD600 of individual expression devices. Red curves were the simulated result based on the least squares to estimate parameters.<br>
 +
<br><div><img src = "https://static.igem.org/mediawiki/2011/1/18/P-5.jpg" width="300"></div>
 +
<br><b>Table 2: </b> The growth constant ks for GFP expression devices
 +
<br>The defined promoter-RBS devices carried out different protein level and comprised different growth constant with their corresponding cell growth rate. We used equation (2) to identify the growth constant of devices via experimental data fitting (figure 4). Each promoter-RBS device corresponded to their growth constant.<br>
 +
 
 +
<p><br>The Ptet promoter (R0040) combined with RBS (B0034) had higher GFP protein expression level and had the smaller ks i.e., slower growth rate. The results indicated that protein expression level of a genetic circuit is negatively correlated to the growth rate of a host bacterium. Bacteria with a low metabolic load have greater growth rate.
 +
The protein expression ability of promoter-RBS devices is indirectly measured by fluorescence intensity of GFP expression devices (Fig. 1 and 3). Since the protein concentration per cell is diluted due to the cell growth and protein degradation, a dynamic model to character the protein expression ability of promoter-RBS devices can be built as follows:<br></p>
 +
<div><img src = "https://static.igem.org/mediawiki/2011/8/81/P-6.jpg" width="300"></div> (2)
 +
<p><br>Where  denotes the protein expression level of a GFP expression device (Fig. 1). p(s) denotes the protein expression ability of a promoter-RBS device, which is the function of cell density parameter s (optical density at 600 nm (O.D. 600)), and rGFP denotes the protein degradation rate of GFP. g(s) denotes the cell growth rate and can be determined from Eq. (1) as:<br></p>
 +
<div><img src = "https://static.igem.org/mediawiki/2011/a/ad/P-7.png" width="300"></div>
 +
(3)
 +
<p><br>The cell growth rate g(s) can be calculated directly by experimental data of O.D. 600 (Fig.4). The negative correlation between cell growth rate g(s) and O.D. 600 indicated that bacteria loss their growth potential gradually when grow in the nutrient limited condition (Fig.5).<br></p>
 +
<div><img src = "https://static.igem.org/mediawiki/2011/6/63/Fig5.jpg" width="700"></div>
 +
<br><b>Figure 5: The growth rate changed with cell density (O.D. 600) in GFP expression devices.</b> Blue dotted lines were the growth rate of individual expression devices with different promoters and RBSs. Red curves were the simulated result based on the least squares estimated parameters in table 2. Figure labels corresponded to table 1. The growth rate showed a significant decreasing current with cells growing. In the saturated stage, cell density reached to the maxima concentration and the cell growth rate almost decreased to 0. All of the data represented the average of three independent measurements. Error bars indicated standard deviations. X-axis indicated the cell density units, and Y-axis indicated the growth rate with the different scales.<br>
 +
<p><br>To calculate the protein expression ability of a promoter-RBS devices p(s) , we assume p(s) is related to cell density parameter s (optical density at 600 nm (O.D. 600)), and a polynomial Eq. (4) is used to trace p(s) as fellow:<br>
 +
p(s) = p<sub>0</sub>+p<sub>1</sub>s &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  (4)
 +
<br>Where p<sub>0</sub> denotes zero-order coefficient, p1 denotes first-order coefficient. The experimental data and trace results were showed in Fig.6. Parameters of our model were determined by the nonlinear regression based on the minimum sum of squired residuals.<b> This dynamic model captured the experimentally observed differences for each genetic circuit (Fig.6).  The fitting results indicate this model is reasonable. A universal mathematical model was found to explain the protein expression activity of Promoter_RBS devices change affected by the cell density (OD600). As we expect, promoter-RBS activity decrease linearly during cell growth (Fig. 7)</b>. Similarly, promoter-RBS activity also decreases with time (Fig. 8). <br></p>
 +
<div><img src = "https://static.igem.org/mediawiki/2011/9/9d/Fig6.jpg" width="700"></div>
 +
<br><b>Figure 6:</b> The simulated result of 12 sets of promoter-RBS expression devices. Blue dotted lines were the fluorescence intensities of individual expression devices with different promoters and RBSs. Red curves were the simulated result based on the least squares to estimate parameters in Table 2. Figure labels corresponded to table 1. And all of the data represented the average of three independent measurements. Error bars indicated standard deviations. X-axis indicated the time units, and Y-axis indicated the fluorescence units with different scales. <br>
 +
<br><div><img src = "https://static.igem.org/mediawiki/2011/d/d9/Fig7.jpg" width="700"></div>
 +
<br><b>Figure 7:</b> The activity of different promoter-RBS devices changed with cell density (O.D. 600). We proposed that the activity of promoter-RBS devices will decrease and change with cell growing. Figure labels corresponded to table 1. And all of the data represented the average of three independent measurements. X-axis indicated the cell density units, and Y-axis indicated the device expressed units with different scales.<br>
 +
<br><div><img src = "https://static.igem.org/mediawiki/2011/e/eb/Fig8.jpg" width="700"></div>
 +
<br><b>Figure 8:</b> The activity of different promoter-RBS devices changed with time. We proposed that the activity of promoter-RBS devices will decrease and change with cell growing. Figure labels corresponded to table 1. And all of the data represented the average of three independent measurements. X-axis indicated the cell density units, and Y-axis indicated the device expressed units with different scales.
 +
<p><br>In previous test, we also assumed p(s) is a polynomial as Eq. (5),<br></p>
 +
<div><img src = "https://static.igem.org/mediawiki/2011/f/fd/P-8.jpg" width="300"></div> (5)
 +
<p><br>but the analysis of extra sum of squares was used to compare this alternative. The Eq. (5) did not fit the data significant better than Eq. (4). We concluded that Eq. (5) was overparametrized, and protein expression activity of promoter_RBS can be described simply as Eq. (4).<br></p>  
-
<p> In order to characterize this high temperature induced device <a href=" http://partsregistry.org/Part:BBa_K098995">BBa_K098995</a>, the fluorescence intensity of<a href=" http://partsregistry.org/Part:BBa_K098988">BBa_K098988</a> is measured by the flow cytometry (Figure. 1).
+
<br><b>Reference</b>  
-
<div><img src = https://static.igem.org/mediawiki/2011/e/e4/M-1.1.JPG></div>
+
-
<br><b>Figure1.</b> Part <a href=" http://partsregistry.org/Part:BBa_K098988">BBa_K098988</a> Design. The heat induced device  <a href=" http://partsregistry.org/Part:BBa_K098995">BBa_K098995</a>uses gene <a href=" http://partsregistry.org/Part:BBa_K098997">BBa_K098997</a> coding for cI repressor to inhibit the cI promoter <a href=" http://partsregistry.org/Part:BBa_R0051">BBa_R0051</a>. The activity of cI repressor is decreased by elevating temperature from 30 ℃ to 42 ℃.
+
-
A differential equation is used to calculate protein expression activity of <a href=" http://partsregistry.org/Part:BBa_K098995">BBa_K098995</a>  as follows.<br><br>
+
 +
<br>1. You, L., et al., Programmed population control by cell-cell communication and regulated killing. Nature, 2004. <br>428(6985): p. 868-71.
 +
<br>2. RICHARDS, F.J., A Flexible Growth Function for Empirical Use. Journal of Experimental Botany, 1959. 10(2): p. 290-301.
 +
</p>
 +
</div>
-
<div><img src = "https://static.igem.org/mediawiki/2011/e/e3/M-5.JPG" width="450"></div>
+
<br><br><br>
-
<p> This equation describes the concentration of GFP in <a href=" http://partsregistry.org/Part:BBa_K098988">BBa_K098988</a>change with time (Figure. 1). Alpha-Temp is the protein expression rates corresponding to <a href=" http://partsregistry.org/Part:BBa_K098995">BBa_K098995</a>which is a temperature sensitive expression device. To describe transition during log phase and stationary phase, the alpha-Temp is assumed to zero in stationary phase. Gamma-GFP are decay rates of the GFP proteins. When bacteria divide, the molecular in a bacterium will be dilute. Because bacteria grow faster, the dilution rate d(t) is included in this model and can be calculated from OD ratio of medium (Figure. 2). The values of the kinetic parameters used in the simulation were initially obtained from the literature and experimental data. Data computations were performed with Matlab software. A program was written and used as a subroutine in Matlab for parameter optimization using nonlinear regression (Figure. 3).</p>
+
</body>
-
<div><img src = "https://static.igem.org/mediawiki/2011/e/e8/M-2.jpg" width="450"></div>
+
</html>
-
<br><b>Figure 2. </b> The OD ratio is increased faster in log phase than it in stationary phase. The dilution rate d(t) can be calculated from OD ratio and used in out model.
+
-
<div><img src = "https://static.igem.org/mediawiki/2011/c/c2/M-3.jpg" width="450"></div>
+
-
<br><b>Figure 3. </b>The behavior of high temperature induced device <a href=" http://partsregistry.org/Part:BBa_K098988">BBa_K098988</a> at 25°C, 37 °C and 42°C. Experimental data (dot) and simulated results (line) of the model  suggest this temperature-dependent device can control the expression level of the target protein by the host cell’s incubation. The fitting results indicate our dynamic model can quantitatively assess the protein expression activity of <a href=" http://partsregistry.org/Part:BBa_K098988">BBa_K098988</a>during log phase and stationary phase.
+
-
<p>Using least squares estimation from experimental data, the relative the protein expression activity of <a href=" http://partsregistry.org/Part:BBa_K098988">BBa_K098988</a>  at 25°C, 37 °C and 42°C were estimated (Figure. 4).</p>
+
-
<div><img src = "https://static.igem.org/mediawiki/2011/b/b2/M-4.JPG" width="450"></div>
+
-
<br><b>Figure 4. </b>The relative the protein expression activity of  <a href=" http://partsregistry.org/Part:BBa_K098988">BBa_K098988</a>at 25°C, 37 °C and 42°C estimated using least squares estimation from experimental data. The protein expression activity at 42°C is higher than 25°C, 37 °C
+
-
<p>According to the fitting results (Figure. 3), the dynamic model successfully approximated the behavior of our high-temperature induced system. The model equation presents interesting mathematical properties that can be used to explore how qualitative features of the genetic circuit depend on reaction parameters. This method of dynamic modeling can be used to guide the choice of genetic ‘parts’ for implementation in circuit design in the future.</p>
+
-
<br><b>References </b>
+
-
<p>Alon, U. (2007) An Introduction to Systems Biology: Design Principles of Biological Circuits. Chapman & Hall/CRC.</p>
+

Latest revision as of 15:54, 5 October 2011



Measurement

This new measurement method can calculate the protein expression rate in the different E. coli population density for the protein expression device (=Promoter_RBS device).


Summary

In previous studies, the transcriptional strengths of promoter and transcriptional strengths of RBSs was defined as a constant values. But in biological concepts, we know the expression rates of most proteins decreased dramatically while the bacteria at the stationary phase. To overcome this problem, our team developed a new measurement method can calculate the protein expression activity of promoter_RBS device changes with cell density in the culture tube directly. We provide a simple polynomial equation which can describe linear relationship between the protein expression activity p(s) of promoter_RBS device and cell desity s (the value of OD600).

p(s) = p0+p1s


where p0 denotes zero-order coefficient, p1 denotes first-order coefficient. This dynamic model captured the experimentally observed differences of green fluorescence intensity for each protein expression circuit. As we expect, promoter-RBS activity decrease linearly related with cell population density. Therefore, if the P0 and P1 parameters of a promoter_RBS device are characterized, we can calculate the protein expression activity of a promoter_RBS device by OD600 value. Furthermore, this model enable us to rational connect a promote_RBS device with different strength to obtain a target protein expression level in a synthetic genetic circuit.


Motivation
In previous studies to model a protein expression, the transcription rates of promoters and translation rates of ribosome binding sites (RBSs) are defined as a constant value. However, the values are not constant during cell growth. For example, the transcription rates of promoters and the translation rates of RBSs are lager in log phase than bacteria growth in stationary phase. To overcome this problem, we selected four promoters and three RBSs with different regulation strength and constructed 12 protein expression devices which combine promoter, RBS and green fluorescent protein (GFP) in Escherichia coli. The GFP expression levels with time were measured using a flow cytometry, and the experimental data can used to characterize a protein expression rate of a protein expression device which contains a promoters and a RBS. A dynamic model that captured the experimentally observed differences for each protein expression device was developed in this study. Using this method, we can measurement the protein expression rate in the different E. coli population density for the protein expression device.


Methods
Two strategies will be applied in this method:
1.Construction of a Promoter-RBS Library and Assay of GFP Expression
We selected four promoters and three RBSs with different regulation strength and constructed 12 protein expression devices which combine promoter, RBS and green fluorescent protein (GFP) in E. coli (Fig.1 &table.1). The GFP expression level with time was measured using a flow cytometry (Fig.3 &Fig.4)



Figure 1: Combinatorial promoter-RBS architecture reveals constitutive expression.
The reporter gene, green fluorescence protein (GFP) locates downstream of promoter-RBS devices as output. The output expression is controlled by differently combinatorial promoter-RBS devices. Each promoter-RBS device is expected to show the significant expression diversity. The combinatorial library contains 12 sets of promoter-RBS devices. Each set consists of different promoters and RBSs.


Table 1: The name of 12 green fluorescence protein devices and their composed biobricks
A (*) refers to construct in two different backbone.


Figure 2: The green fliorescence intensity changed with time provided the useful hints of the interactiofluorescence intensity per cell was measured using a flow cytometer.
The O.D.ratio and the fluorescence data n function between promoters and RBSs.


Figure 3: The time-course expression results were measured with time for 12 sets of GFP expression devices.
(a) promoter-RBS devices which contained the same promoter J23105 combined with three different RBSs (B0030, B0032, and B0034) were selected for measuring. (b) promoter-RBS devices which contained the same promoter J23106 combined with three different RBSs (B0030, B0032, and B0034) were selected for measuring. (c) promoter-RBS devices which contain the same promoter J23114 combined with three different RBSs (B0030, B0032, and B0034) were selected for measuring. (d) promoter-RBS devices which contained the same repressible promoter R0040 combined with three different RBSs (B0030, B0032, and B0034) were selected for measuring. Each measuring even was detected every 15 minutes. And all of the data represented the average of three independent measurements. Error bars indicated standard deviations. X-axis indicated the time units, and Y-axis indicated the fluorescence units with different scales. Furthermore, the fluorescent signals changed with time per cell were measured by using a flow cytometer.


2.Build a Mathematical Model to Characterize Protein Expression Ability of Promoter-RBS Devices.
The logistic growth model is commonly used to describe the bacterium growth curve under the nutrient limited condition as follow [1-2]:

(1)


Where s denotes the cell density parameter (optical density at 600 nm (O.D. 600)), smax denotes the maximum value of cell density, and ks denotes growth rate constant. The cell density time course data of the bacteria with different GFP expression devices were used to solve Eq. (1) (Fig. 4). The growth rate constant ks was calculated and listed in Table 2.


Figure 4:The bacterial population simulated results of 12 sets of promoter-RBS expression devices. Blue dotted lines were the value of OD600 of individual expression devices. Red curves were the simulated result based on the least squares to estimate parameters.


Table 2: The growth constant ks for GFP expression devices
The defined promoter-RBS devices carried out different protein level and comprised different growth constant with their corresponding cell growth rate. We used equation (2) to identify the growth constant of devices via experimental data fitting (figure 4). Each promoter-RBS device corresponded to their growth constant.


The Ptet promoter (R0040) combined with RBS (B0034) had higher GFP protein expression level and had the smaller ks i.e., slower growth rate. The results indicated that protein expression level of a genetic circuit is negatively correlated to the growth rate of a host bacterium. Bacteria with a low metabolic load have greater growth rate. The protein expression ability of promoter-RBS devices is indirectly measured by fluorescence intensity of GFP expression devices (Fig. 1 and 3). Since the protein concentration per cell is diluted due to the cell growth and protein degradation, a dynamic model to character the protein expression ability of promoter-RBS devices can be built as follows:

(2)


Where denotes the protein expression level of a GFP expression device (Fig. 1). p(s) denotes the protein expression ability of a promoter-RBS device, which is the function of cell density parameter s (optical density at 600 nm (O.D. 600)), and rGFP denotes the protein degradation rate of GFP. g(s) denotes the cell growth rate and can be determined from Eq. (1) as:

(3)


The cell growth rate g(s) can be calculated directly by experimental data of O.D. 600 (Fig.4). The negative correlation between cell growth rate g(s) and O.D. 600 indicated that bacteria loss their growth potential gradually when grow in the nutrient limited condition (Fig.5).


Figure 5: The growth rate changed with cell density (O.D. 600) in GFP expression devices. Blue dotted lines were the growth rate of individual expression devices with different promoters and RBSs. Red curves were the simulated result based on the least squares estimated parameters in table 2. Figure labels corresponded to table 1. The growth rate showed a significant decreasing current with cells growing. In the saturated stage, cell density reached to the maxima concentration and the cell growth rate almost decreased to 0. All of the data represented the average of three independent measurements. Error bars indicated standard deviations. X-axis indicated the cell density units, and Y-axis indicated the growth rate with the different scales.


To calculate the protein expression ability of a promoter-RBS devices p(s) , we assume p(s) is related to cell density parameter s (optical density at 600 nm (O.D. 600)), and a polynomial Eq. (4) is used to trace p(s) as fellow:
p(s) = p0+p1s             (4)
Where p0 denotes zero-order coefficient, p1 denotes first-order coefficient. The experimental data and trace results were showed in Fig.6. Parameters of our model were determined by the nonlinear regression based on the minimum sum of squired residuals. This dynamic model captured the experimentally observed differences for each genetic circuit (Fig.6). The fitting results indicate this model is reasonable. A universal mathematical model was found to explain the protein expression activity of Promoter_RBS devices change affected by the cell density (OD600). As we expect, promoter-RBS activity decrease linearly during cell growth (Fig. 7). Similarly, promoter-RBS activity also decreases with time (Fig. 8).


Figure 6: The simulated result of 12 sets of promoter-RBS expression devices. Blue dotted lines were the fluorescence intensities of individual expression devices with different promoters and RBSs. Red curves were the simulated result based on the least squares to estimate parameters in Table 2. Figure labels corresponded to table 1. And all of the data represented the average of three independent measurements. Error bars indicated standard deviations. X-axis indicated the time units, and Y-axis indicated the fluorescence units with different scales.


Figure 7: The activity of different promoter-RBS devices changed with cell density (O.D. 600). We proposed that the activity of promoter-RBS devices will decrease and change with cell growing. Figure labels corresponded to table 1. And all of the data represented the average of three independent measurements. X-axis indicated the cell density units, and Y-axis indicated the device expressed units with different scales.


Figure 8: The activity of different promoter-RBS devices changed with time. We proposed that the activity of promoter-RBS devices will decrease and change with cell growing. Figure labels corresponded to table 1. And all of the data represented the average of three independent measurements. X-axis indicated the cell density units, and Y-axis indicated the device expressed units with different scales.


In previous test, we also assumed p(s) is a polynomial as Eq. (5),

(5)


but the analysis of extra sum of squares was used to compare this alternative. The Eq. (5) did not fit the data significant better than Eq. (4). We concluded that Eq. (5) was overparametrized, and protein expression activity of promoter_RBS can be described simply as Eq. (4).


Reference
1. You, L., et al., Programmed population control by cell-cell communication and regulated killing. Nature, 2004.
428(6985): p. 868-71.
2. RICHARDS, F.J., A Flexible Growth Function for Empirical Use. Journal of Experimental Botany, 1959. 10(2): p. 290-301.