Latest revision as of 03:15, 22 October 2019

Document

1. Protein-split Model

Introduction

Proteins are composed of relatively independent subunits or secondary structures. After subtilisin-modified ribonuclease was split and reassembled in 1958 multifarious proteins was split successfully certificated by experiments[1]. However, different split sites produce two different parts, and the assembly effect will be different. Renilla luciferase, wildely used in the reporter systems is split in our experiment, the two parts of the split Renilla luciferase need to be self-assembled under the guidance of the optical protein switch. The split site of Renilla luciferase is choosed from the result of experiments by Paulmurugan R and Gambhir S S[2], which is time-consuming and laborious, so can we predict better sites with calculational methods? We have searched a large number of previous experimental results which hint some kind of regularities (Table.1). The split site is closely related to the secondary structure of the protein, and the inappropriate site is generally located in the loop region of the protein. In addition, the location of the split sites seems to be related to the factors such as Sequence conservation, Length of coil regions, and Center distance. Based on these rules Protein-split Model is established to predict better protein-split positions. We chose nanoLuc protein in the iGEM library to verify our model. The combination of the two nanoLuc parts obtained by the split position of our model is better than the original one verified by our experiment that means our model is convincing(Fig.4, Table.2, Table.3).

Design

We conclude that the protein-split sites may be related to three factors, Secondary structure (factor 1), Length of coil region(factor 2), Center distance(factor 3)(Fig.1), and found 65 split sites (Table.1).

Factor1: Protein is cut in coil region which is classified as index 1, in the same way, the boundary between coil region and α-helix is classified as index 2, the boundary between coil region and β-strand is classified as index 3，in α-helix or β-strand is classified as index 4. Factor2：If protein-split site is one of the first three types, the factor value was defines by counting the length of the coil region. The length of the coil region is 0 if protein-split site is the index 4 type. Factor3: Center distance was defined as the ratio of the shorter fragment length to half of the full length of the protein(center distence=sublength/(0.5×fulllength).

We can use a typical BP neural network (Dichotomous model) to help us find split sites. If we can find a suitable set of weights{w1, w2}, at the same time, the inputs (I1, I2, I3) = (Factor1, Factor2, Factor3), output(O1, O2)=(1,0), then the protein-split sites selection is actually solved. For any protein, we only need to download its secondary structure from the PDB database, each amino acid residue corresponds to a set of factors. As long as the factors are input into the neural network, a label will be output: T:1:[1, 0] ; F: 0: [0, 1] (Fig. 2). We find a set of appropriate weights {w1, w2} by training the sample to complete our model, and do the sensitivity and specification analysis.

Fig.2 Feedforward neural network(FNN)

Algorithm

See connected protein-split_supplement1.

Database

Input see connected protein-split_supplement2,Output see connected protein-split_supplement3.

Result

Weight matrix:

w1

1	2	3
-0.0306	0.233	0.5924
0.6037	-0.6431	-0.5429
1.3603	-0.6622	0.5412
0.6023	-0.3277	-0.5616

b1

0.7124	-0.464	-0.9178	-0.5768

w2

1	2
0.3485	-0.1632
-0.2581	0.1562
-0.8925	0.9408
-0.2699	-0.1891

b2

0.7532	0.1274

We sample 20 samples as test sets, and the results are as follows.

Fig.3 Confusion matrix

Performance Evaluation Criteria

Conclusion

The model presents to be highly sensitive, but has low specificity which produces high ratios of false positives, meaning that many inappropriate protein-split sites are recognized as appropriate protein-split sites (Fig.3). However, in general, appropriate split sites are much less than inappropriate split sites in a protein with hundreds of amino acid residues, so the low specificity is inevitable. We only need to ensure that very few appropriate sites are not mistaken for inappropriate sites. This model will provide several predicted appropriate sites to instruct experiment, as long as there are one true appropriate sites. Considering our model has a very low false negative, this model is fabulous in practical applications.

Experimental Verification

We chose nanoLuc in iGEM library to split and test our model. Based on the protein-split model, we calculated an optimal split site and predicted that the original split site is not ideal enough (Fig.4). The experimental results showed that the ability of the two parts of nanoLuc(st-1 and sc-1) cut at our predicted split site to self-assemble is much greater than the original split site(st-2 and sc-2)(Table.2 and Table.3).

Fig.4 Predicted split site and the original split site in the sequence of nanoLuc.

Reference

[1]Richards F M. On the enzymic activity of subtilisin-modified ribonuclease[J]. Proceedings of the National Academy of Sciences of the United States of America, 1958, 44(2): 162.
[2]Paulmurugan R, Gambhir S S. Monitoring protein− protein interactions using split synthetic renilla luciferase protein-fragment-assisted complementation[J]. Analytical chemistry, 2003, 75(7): 1584-1589.
[3]Han T, Chen Q, Liu H. Engineered photoactivatable genetic switches based on the bacterium phage T7 RNA polymerase[J]. ACS synthetic biology, 2016, 6(2): 357-366.
[4]Ghosh I, Hamilton A D, Regan L. Antiparallel leucine zipper-directed protein reassembly: application to the green fluorescent protein[J]. Journal of the American Chemical Society, 2000, 122(23): 5658-5659.
[5]Taniuchi H, Anfinsen C B, Sodja A. Nuclease-T: an active derivative of staphylococcal nuclease composed of two noncovalently bonded peptide fragments[J]. Proceedings of the National Academy of Sciences of the United States of America, 1967, 58(3): 1235.
[6]de Prat Gay G, Ruiz-Sanz J, Davis B, et al. The structure of the transition state for the association of two fragments of the barley chymotrypsin inhibitor 2 to generate native-like protein: implications for mechanisms of protein folding[J]. Proceedings of the National Academy of Sciences, 1994, 91(23): 10943-10946.
[7]Bae C, Suchyna T M, Ziegler L, et al. Human PIEZO1 ion channel functions as a split protein[J]. PloS one, 2016, 11(3): e0151289.
[8]Johnsson N, Varshavsky A. Split ubiquitin as a sensor of protein interactions in vivo[J]. Proceedings of the National Academy of Sciences, 1994, 91(22): 10340-10344.
[9]Pelletier J N, Campbell-Valois F X, Michnick S W. Oligomerization domain-directed reassembly of active dihydrofolate reductase from rationally designed fragments[J]. Proceedings of the National Academy of Sciences, 1998, 95(21): 12141-12146.
[10]Galarneau A, Primeau M, Trudeau L E, et al. β-Lactamase protein fragment complementation assays as in vivo and in vitro sensors of protein–protein interactions[J]. Nature biotechnology, 2002, 20(6): 619.
[11]Massoud T F, Paulmurugan R, Gambhir S S. A molecularly engineered split reporter for imaging protein-protein interactions with positron emission tomography[J]. Nature medicine, 2010, 16(8): 921.
[12]Martell J D, Yamagata M, Deerinck T J, et al. A split horseradish peroxidase for the detection of intercellular protein–protein interactions and sensitive visualization of synapses[J]. Nature biotechnology, 2016, 34(7): 774.
[13]Kaihara A, Kawai Y, Sato M, et al. Locating a Protein− Protein Interaction in Living Cells via Split Renilla Luciferase Complementation[J]. Analytical chemistry, 2003, 75(16): 4176-4181.
[14]Paulmurugan R, Gambhir S S. Novel fusion protein approach for efficient high-throughput screening of small molecule–mediating protein-protein interactions in cells and living animals[J]. Cancer research, 2005, 65(16): 7413-7420.
[15]Kim S B, Kanno A, Ozawa T, et al. Nongenomic activity of ligands in the association of androgen receptor with SRC[J]. ACS chemical biology, 2007, 2(7): 484-492.
[16]Luker K E, Smith M C P, Luker G D, et al. Kinetics of regulated protein–protein interactions revealed with firefly luciferase complementation imaging in cells and living animals[J]. Proceedings of the National Academy of Sciences, 2004, 101(33): 12288-12293.
[17]Paulmurugan R, Umezawa Y, Gambhir S S. Noninvasive imaging of protein–protein interactions in living subjects by using reporter protein complementation and reconstitution strategies[J]. Proceedings of the National Academy of Sciences, 2002, 99(24): 15608-15613.
[18]Wehr M C, Laage R, Bolz U, et al. Monitoring regulated protein-protein interactions using split TEV[J]. Nature methods, 2006, 3(12): 985.

2 Microalgae movement model

Introduction

Our project aims to take advantage of the phototaxis characteristics of microalgae to construct a drug carrier which can be directed to the diseased cells under the guidance of light. We first analyzed the movement of microalgae in water. Using Langevin equation to describe the movement of microalgae in water, and then according to the observation of the pause of microalgae movement (movie). We called it "breaststroke", and established the microalgae "breaststroke" model. According to Langevin equation, the viscosity in water can be replaced by the viscosity in blood to imitate the movement in blood. In addition, since red blood cells account for 45% of total volume of blood, the collision between microalgae and red blood cells will change the energy of microalgae. Since the collision between red blood cells and microalgae is random and complex, we assume that the red blood cells and microalgae are spheres. Under this assumption, the energy loss distribution shall conform the three-dimensional trigonometric function, so we established the microalgae movement model. Our model is the first to systematically describe the movement of microalgae in the water and blood. It has great originality and application value.

Ⅰ.The movement of microalgae in water

Design

Parameters

According to Langevin equation, the motion of microalgae in water can be described as（Fig.1）：

By formula transformation, (1) can be simplified as:

Defined：

Fig.1 The movement of microalgae in water

According to the phenomenon that we observed in the experiment: microalgae do not move at a constant speed, it pauses with periodicity during the movement(movie). We assume that the flagellum of microalgae generates a huge thrust, which can make microalgae generate an initial velocity in a very short time. The motion of microalgae in water

Simultaneous equations:

Therefore the motion curve of the microalgae in the water can be determined.

Caculation

(1)Solve : The average velocity of microalgae can be obtained from the experiment, then the distance in one cycle can be obtained:

(2)Solve the expression of :v⁰= 1.96mm/s.

Fig.2 Periodic curve of microalgae movement in water

Ⅱ. The movement of microalgae in the blood

Parameters

Design

In the blood, since the speed of movement of microalgae is related to the light intensity, we can control the light intensity to control the speed of movement of microalgae. Through experiments, we measured the average speeds of microalgae movement under different light conditions in the water（Fig.3）：

Fig.3 The relationship between the microalgae movement speeds and light intensities.

The equation obtained through experimental fitting:

The energy of microalgae in the blood without collision:

Following we will consider the collision in the blood (Fig.4):

Fig.4 Microalgae as a drug carrier in the blood against the flow of blood.

As the figure above, if microalgae swim upstream in the blood, they will be collided by red blood cells when they move in the blood. So, the speed of microalgae movement may change at any time. We find that the percentage of red blood cells in the blood vessels is about 45%. Which means if we take a cross-section of the blood vessels, the percentage of red blood cells is 45%(Fig.5).

Fig.5 cross-section of the blood vessels

Since red blood cells are randomly distributed in blood vessels, we assume that they are also randomly distributed in every small area of every cross section.

Thus, we equivalent the number of red blood cells colliding one microalgae cross section:

The worst case scenario is a head-on collision between the microalgae and the red blood cells. In this case, microalgae lost the most energy. Assume the red blood cell energy is 0 after the collision. By energy conservation, microalgae lost the same energy with red blood cells. While, during the movement of microalgae, the collision angle between microalgae and red blood cells is random. Therefore, we use three-dimensional trigonometric functions to describe the energy loss of microalgae and red blood cells(Fig.6).

Fig.6 The energy loss of microalgae

Through the operation of double integral, we obtained the energy lost of each section:

The energy of microalgae after collision in the blood:

If the microalgae need to swim upstream in the blood, then:

Calculation

(1)Calculate N: By taking the radius of red blood cells and microalgae, the following equation can be obtained:

(2)Calculate Q: The energy they lost in unit time of each blood vessel section:

(3)We can obtain the relationship between the light intensity and different blood flow velocity we set（Fig.7）:

Fig.7 Luminous intensity threshold values that allows microalgae to swim upstream in the blood at different blood flow velocities.

Conclusion

(1) We are the first to describe the microalgae periodic motion phenomenon in water using a precise mathematical expression.

(2) We obtained the minimum light intensity thresholds required to control microalgae to swim as intended at different blood flow velocities. In the future, microalgae may act as a drug carrier. Our model may provide some valuable data and rules for clinical trials.

3 Protein-protein interaction model

Introduction

Protein interaction is widespread in all live units, and its importance is self-evident. In our design, red light can stimulate the combination of N-hrluc-PhyB and C-hrluc-PIF3 and then produce blue light to guide microalgae to move(Fig.1).Therefore, whether our design is reasonable and the experiment successes depends on the degree of combination of N-hrluc-PhyB and C-hrluc-PIF3. Based on chemical thermodynamic formula, the relationship between binding rate and concentration and affinity constant is established in two steps. In addition, as a commonly used light-controlled polymeric protein actuator, this model can not only guide our own experiment, but also provide a quantitative model for the field of optogenetics.

Fig.1 A: The combination of N-hrluc and C-hrluc without the assistance by optical protein switch(PhyB&PIF3). B: The combination of N-hrluc and C-hrluc with the assistance by optical protein switch.

Design

The molecular reaction in the solution can be described that the binding of PhyB and PIF3 stimulated by red light facilitates the binding of the two parts(N-Hrluc, C-hrluc) of the Renilla luciferase. Therefore this molecular reaction in solution can be divided into two stages:

StageⅠ: only the binding of PhyB and PIF3 is considered.

In order to find out the effect of the initial concentration of PhyB and PIF3 on the binding rate of the two proteins, we get a specific equation between the initial concentration of protein and the binding rate of two proteins through dissociation constant Kd.

Parameters

Symbol	Meaning
A	Concentration of PhyB
B	Concentration of PIF3
K_a	Affinity constant
K_d	Dissociation constant
η	Initial concentration of PhyB
A₀	Initial concentration of PhyB
B₀	Initial concentration of PIF3
AB	Concentration of PhyB-PIF3

Hypothesis:

After the binding of two proteins, the solution contains three substances: protein PhyB, protein PIF3 and the combination of them. That is, the material is conserved and there is no protein loss during the process of the reaction.

Calculation

Considering the binding of PhyB and PIF3 in solution：

If the free concentration of PhyB is A, the free concentration of PIF3 is B, and the concentration of PhyB&PIF3 is AB, the binding constant :

The dissociation constant :

Then the expression of the binding rate:

By the (2) (3) (4) simultaneous:

We assume that the initial concentration of PhyB protein is higher than that of PIF3 protein, that is,

Available jointly by (5)(6):

The specific equation between the initial concentration of the two proteins and the binding rate of the two proteins was obtained.
Available from (7):

When the light with the wavelength of 650 nm and a frequency of 20 Hz, we know Kd=500nM [1]（Fig.1）：

Fig.2 Binding Affinity(Kd ) of PIF3&PhyB

Calculation

(1) Similar to the barrel effect, the binding rate depends on the low concentration of protein, and it can be estimated from the diagram that the minimum threshold of the concentration of the two proteins is "1" "0" ^"-6" mol /L when the concentration of the two proteins is equal(Fig.3 and Fig.4).

Fig.3 The binding rate of PhyB and PIF3 related to the concentration.

Fig.4 The trend of binding rate when PhyB and PIF3 are at the same concentration

(2)If the combination rate of A and B is not considered under red light (dark), then

Under the excitation of red light, the value K_d is drastically reduced due to the change of in the configuration of the protein, only whenK_d (dark)/K_d (hv)>10000, the combination rate in the dark can be ignored(Fig.5).

Fig.5 The net increase of the binding rate under

red light eliminated the dark background effect.

StageⅡ: the combination of N-hrluc and C-hrluc.

Under the action of "optical switch"（PhyB&PIF3）, the distance between N-hrluc and C-hrluc becomes closer, which is equivalent to the increasing of concentration. The effect of "optical switch" on protein binding is converted to value by equal.

Parameters

Parameter	Description
B_i,i=1,2	N-hrluc and C-hrluc protein parts respectively
A_i,i=1,2	PhyB and PIF3 respectively
d	Distance between B₁ and B₂

Hypothesis:

(1)The addition of optical switch has no effect on the conformational of hrluc.
(2)After the optical switch bound, the original protein is pulled closer, that is, the distance is shorten.

Calculation

Fig.6 A: The combination of N-hrluc and C-hrluc without the assistance by protein switch(PhyB&PIF3). B: The combination of N-hrluc and C-hrluc with the assistance by protein switch

As shown in the following figure(Fig.6), optical switches A_1,A_2are added to N-hrluc and C-hrluc proteins. Under the irradiation of light, the optical switches A_1,A_2begin to combine. Due to the hypothesis that conformational structure of B1 and B2 does not change, so the affinity constant is unchanged, while the distance of B1 and B2 is reduced, which is equivalent to the increase of concentration.

The concentration of a particle can be obtained as follows:

The volume can be defined by the spacing between proteins, that is,

So:

According to the size of the protein switch(PhyB&PIF3), in the least ideal case, the distance between N-hrluc and C-hrluc is 16 nm. Assuming that the spacing is 15 nm, 10 nm and 5 nm, the equivalent concentrations can be caculated as 3×"1" "0" ^"-3" mol/L, 1×"1" "0" ^"-2" mol/L, 8×"1" "0" ^"-2" mol/L according to formula (12).

Fig.7 N-Luc and C-Luc are catched by the optical protein switch in spatial scale which means a higher binding rate (than being dispersed in solution)

As long as the affinity constant of the two parts of fluorescein enzyme is less than ,it can be considered to be completely binding(Fig.7).

Conclusion

(1) The optical switch proteins(PhyB and PIF3) will bind only when the concentration of optical switch protein reach 10-6mol/L according to our model(Fig.4).

(2) The equivalent concentration of renin luciferase depends on the structure size of optical protein switch(Fig.7).

Reference

[1]Levskaya A, Weiner O D, Lim W A, et al. Spatiotemporal control of cell signalling using a light-switchable protein interaction[J]. Nature, 2009, 461(7266): 997.

@@ Line 1: / Line 1: @@
 {{DUT_China_B/test_style.css}}
 {{Template:DUT_China_B/nav}}
+<!-- <!Doctype html> -->
 <html lang="en">
 <head>
-	<meta charset="UTF-8">
+  <meta charset="UTF-8">
-	<title>Document</title>
+  <title>Document</title>
+  <link rel="stylesheet" type="text/css" href="css/default.css">
-	<style>
+  <style>
-	*{
+  *{
-		margin: 0;
+    margin: 0;
-		padding: 0;
+    padding: 0;
-	}
+  }
-	body{
+  body{
-		background-color: #eee;
+    background-color: FFF;
-	}
+  }
-	.firstimg{
+  .firstimg{
-		width: 100%;
+    width: 100%;
-	}
+  }
-	.firstimg img{
+  .firstimg img{
-		width: 100%
+    width: 100%
-	}
+  }
-	#maintest{
+  #maintest{
-		position: relative;
+    position: relative;
-		float: right;
+    float: right;
-		margin: 0px;
+    margin: 0px;
-		width: 76%;
+    width: 76%;
-		padding: 3%;
+    padding: 3%;
                  font-family: 'Times New Roman' !important;
                  border: solid #D3D3D3 1px;
-		/*border-left: 2px solid #ccc;*/
+    /*border-left: 2px solid #ccc;*/
-		/*clear: both;*/
+    /*clear: both;*/
-	}
+  }
         #maintest p{
         font-size: 20px;
@@ Line 53: / Line 54: @@
        height:100%;
        }
-	/*侧边栏样式*/
+  /*侧边栏样式*/
-	#sides{
+  #sides{
-		width: 22%;
+    width: 22%;
-		padding: 20% 0 10%;
+    padding: 20% 0 10%;
                  background-color:#1C2B42;
-		float: left;
+    float: left;
-	}
+  }
-	#sides a{
+  #sides a{
-		text-decoration: none;
+    text-decoration: none;
-	        color:#FFFFFF !important;
+          color:#FFFFFF !important;
-		font-size: 60px !important；
+    font-size: 60px !important；
                  positon:fixed;
                  top:20px;
-	}
+  }
-	div.table-title {
+  div.table-title {
             display: block;
              margin: auto;
@@ Line 206: / Line 207: @@
-	</style>
+  </style>
 </head>
 <body>
-	<div class="firstimg">
+  <div class="firstimg">
-		<img src="https://static.igem.org/mediawiki/2019/d/dd/T--DUT_China_B--modelbannner1.jpg" alt="parts">
+    <img src="https://static.igem.org/mediawiki/2019/d/dd/T--DUT_China_B--modelbannner1.jpg" alt="parts">
-	</div>
+  </div>
-	<div id="sides">
+  <div id="sides">
 <ul id="menu">
          <li>
-            <a ><font size="4" href="#Inspiration" >Inspiration</font></a>
+         <a href="#Protein-split Model"><font size="5" >Protein-split Model</font></a>
+      <div style="text-align: center; width: 100%; height:40px"></div>
+       </li><br>
+          <li>
+        <a href="#Microalgae movement model"><font size="5" >Microalgae movement model</font></a>
+      <div style="text-align: center; width: 100%; height:40px"></div>
+       </li>
+        <li>
+            <a  href="#Protein-interact model" ><font size="5">Protein-protein interaction model</font></a>
          <div style="text-align: center; width: 100%; height:40px"></div>
          </li>
         <br>
-        <li>
-            <a><font size="4" >Background</font></a>
-      <div style="text-align: center; width: 100%; height:40px"></div>
-       </li>
         <br>
          <li>
-             <a><font size="4">Chlamydomonas reinhardtii</font></a>
+             <a><font size="5"></font></a>
          </li>
@@ Line 231: / Line 237: @@
      </ul>
-	</div>
+  </div>
 <div id="maintest">
+<h1 style="font-family: 'Times New Roman' !important; "><a name="Protein-split Model">1. Protein-split Model</a></h1>
+    <div  class="cart">
+      <h1  style="font-family: 'Times New Roman' !important; "><a name="Inspiration" >Introduction </a><img src="https://static.igem.org/mediawiki/2019/a/ab/T--DUT_China_B--Modelnote.svg" class="icon"> </h1>
+      <p style="font-family: 'Times New Roman' !important;   ">Proteins are composed of relatively independent subunits or secondary structures. After subtilisin-modified ribonuclease was split and reassembled in 1958 multifarious proteins was split successfully certificated by experiments[1]. However, different split sites produce two different parts, and the assembly effect will be different. Renilla luciferase, wildely used in the reporter systems is split in our experiment, the two parts of the split Renilla luciferase need to be self-assembled under the guidance of the optical protein switch. The split site of Renilla luciferase is choosed from the result of experiments by Paulmurugan R and Gambhir S S[2], which is time-consuming and laborious, so can we predict better sites with calculational methods? We have searched a large number of previous experimental results which hint some kind of regularities (Table.1). The split site is closely related to the secondary structure of the protein, and the inappropriate site is generally located in the loop region of the protein. In addition, the location of the split sites seems to be related to the factors such as Sequence conservation, Length of coil regions, and Center distance. Based on these rules Protein-split Model is established to predict better protein-split positions. We chose nanoLuc protein in the iGEM library to verify our model. The combination of the two nanoLuc parts obtained by the split position of our model is better  than the original one verified by our experiment that means our model is convincing(Fig.4, Table.2, Table.3). </p>
+    </div>
+  <!-- 设计 -->
+  <div class="cart">
+      <h1 style="font-family: 'Times New Roman' !important; ">Design<img src="https://static.igem.org/mediawiki/2019/b/b7/T--DUT_China_B--Modeldata.svg" class="icon"> </h1>
+      <p>We conclude that the protein-split sites may be related to three factors, Secondary structure (factor 1), Length of coil region(factor 2), Center distance(factor 3)(Fig.1), and found 65 split sites (Table.1).</P>
+      <div style="text-align: center; width: 100%; ">
+            <img alt="" src="https://static.igem.org/mediawiki/2019/0/07/T--DUT_China_B--modelfigure1.png" style="display: inline-block;width:80%;" />
+            <center> <br> <p style="text-align:center;position:relative;">Factor1: Protein is cut in coil region which is classified as index 1, in the same way, the boundary between coil region and α-helix is classified as index 2, the boundary between coil region and β-strand is classified as index 3，in α-helix or β-strand is classified as index 4.  Factor2：If protein-split site is one of the first three types, the factor value was defines by counting the length of the coil region. The length of the coil region is 0 if protein-split site is the index 4 type.  Factor3:   Center distance was defined as the ratio of the shorter fragment length to half of the full length of the protein(center distence=sublength/(0.5×fulllength).   </p>  </center>
+      </div>
+      <p>We can use a typical BP neural network (Dichotomous model) to help us find split sites. If we can find a suitable set of weights{w1, w2}, at the same time, the inputs (I1, I2, I3) = (Factor1, Factor2, Factor3), output(O1, O2)=(1,0), then the protein-split sites selection is actually solved. For any protein, we only need to download its secondary structure from the PDB database, each amino acid residue corresponds to a set of factors. As long as the factors are input into the neural network, a label will be output: T:1:[1, 0] ; F: 0: [0, 1] (Fig. 2). We find a set of appropriate weights {w1, w2} by training the sample to complete our model, and do the sensitivity and specification analysis.</p>
+      <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/1/1a/T--DUT_China_B--network.png" style="display: inline-block;width:80%;" />
+             <center> <br> <p style="text-align:center;position:relative;">Fig.2 Feedforward neural network(FNN)</p>  </center>
+      </div>
+      <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/9/94/T--DUT_China_B--tabelpage1.png" style="display: inline-block;width:80%;" />
+                <img alt="" src="https://static.igem.org/mediawiki/2019/c/cb/T--DUT_China_B--tabelpage2.png" style="display: inline-block;width:80%;" />
+                <img alt="" src="https://static.igem.org/mediawiki/2019/f/f2/T--DUT_China_B--tabelpage3.png" style="display: inline-block;width:80%;" />
+      </div>
+      <h1 style="font-family: 'Times New Roman' !important; ">Algorithm</h1>
+       <img src="https://static.igem.org/mediawiki/2019/5/53/T--DUT_China_B--Algorithm.jpg">
+      <p style="font-family: 'Times New Roman' !important;">See connected <a href="https://static.igem.org/mediawiki/2019/c/c6/T--DUT_China_B--Protein-split_supplement1.pdf">protein-split_supplement1</a>.</p>
+      <h1 style="font-family: 'Times New Roman' !important; ">Database</h1>
+      <p style="font-family: 'Times New Roman' !important;">Input see connected <a href="https://static.igem.org/mediawiki/2019/3/38/T--DUT_China_B--protein-split_supplement2.xls">protein-split_supplement2</a>,Output see connected <a href="protein-split_supplement3.xlsx">protein-split_supplement3</a>.</p>
-		<div  class="cart">
+      <h1 style="font-family: 'Times New Roman' !important; ">Result</h1>
-               <h1  style="font-family: 'Times New Roman' !important; "><a name="Inspiration" >Introduction </a><img src="https://static.igem.org/mediawiki/2019/9/98/T--DUT_China_B--INSPIRATION.svg" class="icon"> </h1>
+      <h2>Weight matrix:</h2>
+      <h3>w1</h3>
+      <table class="table-fill">
+        <thead>
+        <tr>
+          <th class="text-left">1</th>
+          <th class="text-left">2</th>
+          <th class="text-left">3</th>
+        </tr>
+        </thead>
+        <tbody class="table-hover">
+        <tr>
+          <td class="text-left">-0.0306</td>
+          <td class="text-left">0.233</td>
+          <td class="text-left">0.5924</td>
+        </tr>
+        <tr>
+          <td class="text-left">0.6037</td>
+          <td class="text-left">-0.6431</td>
+          <td class="text-left">-0.5429</td>
+        </tr>
+        <tr>
+          <td class="text-left">1.3603</td>
+          <td class="text-left">-0.6622</td>
+          <td class="text-left">0.5412</td>
+        </tr>
+        <tr>
+          <td class="text-left">0.6023</td>
+          <td class="text-left">-0.3277</td>
+          <td class="text-left">-0.5616</td>
+        </tr>
+        </tbody>
+      </table>
+      <h3>b1</h3>
+      <table class="table-fill">
+        <thead>
+        <tr>
+          <th class="text-left">0.7124</th>
+          <th class="text-left">-0.464</th>
+          <th class="text-left">-0.9178</th>
+          <th class="text-left">-0.5768</th>
+        </tr>
+        </thead>
+      </table>
+      <h3>w2</h3>
+      <table class="table-fill">
+        <thead>
+        <tr>
+          <th class="text-left">1</th>
+          <th class="text-left">2</th>
+        </tr>
+        </thead>
+        <tbody class="table-hover">
+        <tr>
+          <td class="text-left">0.3485</td>
+          <td class="text-left">-0.1632</td>
+        </tr>
+        <tr>
+          <td class="text-left">-0.2581</td>
+          <td class="text-left">0.1562</td>
+        </tr>
+        <tr>
+          <td class="text-left">-0.8925</td>
+          <td class="text-left">0.9408</td>
+        </tr>
+        <tr>
+          <td class="text-left">-0.2699</td>
+          <td class="text-left">-0.1891</td>
+        </tr>
+        </tbody>
+      </table>
+      <h3>b2</h3>
+      <table class="table-fill">
+        <thead>
+        <tr>
+          <th class="text-left">0.7532</th>
+          <th class="text-left">0.1274</th>
+        </tr>
+        </thead>
+      </table>
+      <div style="text-align: center; width: 100%; ">
+     <p>We sample 20 samples as test sets, and the results are as follows.   </p>
+ <img alt="" src="https://static.igem.org/mediawiki/2019/0/0b/T--DUT_China_B--confmat.png" style="display: inline-block;width:80%;" />
+             <center> <br> <p style="text-align:center;position:relative;">Fig.3 Confusion matrix </p>  </center>
+      </div>
+      <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/a/a3/T--DUT_China_B--ROCpara.png" style="display: inline-block;width:80%;" />
+      </div>
+ <h2>Performance Evaluation Criteria</h2>
+  <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/7/78/T--DUT_China_B--Model23.jpg" style="display: inline-block;width:80%;" />
+                <img alt="" src="https://static.igem.org/mediawiki/2019/4/49/T--DUT_China_B--m2formula1.png" style="display: inline-block;width:80%;" />
+                <img alt="" src="https://static.igem.org/mediawiki/2019/e/e5/T--DUT_China_B--m2formula2.png" style="display: inline-block;width:80%;" />
+      </div>
+      <h2>Conclusion</h2><br>
+      <p>The model presents to be highly sensitive, but has low specificity which produces high ratios of false positives, meaning that many inappropriate  protein-split sites are recognized as appropriate  protein-split sites (Fig.3). However, in general, appropriate split  sites are much less than inappropriate split sites in a protein with hundreds of amino acid residues, so the low specificity is inevitable. We only need to ensure that very few appropriate sites are not mistaken  for inappropriate sites. This model will provide several predicted appropriate sites to instruct experiment, as long as there are one true appropriate sites. Considering our model has  a very low false negative, this model is fabulous  in practical applications.
+     </p>
+     <h2>Experimental Verification</h2>
+     <p>We chose nanoLuc in iGEM library to split and test our model. Based on the protein-split model, we calculated an optimal split  site and predicted that the original  split site is not ideal enough  (Fig.4). The experimental results showed that the ability of the two parts of nanoLuc(st-1 and sc-1) cut at our predicted split site to self-assemble is much greater than the original split site(st-2 and sc-2)(Table.2 and Table.3).</p>
+  <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/6/67/T--DUT_China_B--modelfigure2.png" style="display: inline-block;width:80%;" />
+                 <center> <br> <p style="text-align:center;position:relative;">Fig.4 Predicted split site and the original split site in the sequence of nanoLuc.</p>  </center>
+                <img alt="" src="https://static.igem.org/mediawiki/2019/e/e2/T--DUT_China_B--modtabel1.png" style="display: inline-block;width:80%;" />
+                <img alt="" src="https://static.igem.org/mediawiki/2019/b/b6/T--DUT_China_B--modtabel2.png" style="display: inline-block;width:80%;" />
+      </div>
+    </div>
+      <hr>
+<h1 style="font-family: 'Times New Roman' !important; ">Reference</h1>
+  <p>[1]Richards F M. On the enzymic activity of subtilisin-modified ribonuclease[J]. Proceedings of the National Academy of Sciences of the United States of America, 1958, 44(2): 162.<br>
+[2]Paulmurugan R, Gambhir S S. Monitoring protein− protein interactions using split synthetic renilla luciferase protein-fragment-assisted complementation[J]. Analytical chemistry, 2003, 75(7): 1584-1589.<br>
+[3]Han T, Chen Q, Liu H. Engineered photoactivatable genetic switches based on the bacterium phage T7 RNA polymerase[J]. ACS synthetic biology, 2016, 6(2): 357-366.<br>
+[4]Ghosh I, Hamilton A D, Regan L. Antiparallel leucine zipper-directed protein reassembly: application to the green fluorescent protein[J]. Journal of the American Chemical Society, 2000, 122(23): 5658-5659.<br>
+[5]Taniuchi H, Anfinsen C B, Sodja A. Nuclease-T: an active derivative of staphylococcal nuclease composed of two noncovalently bonded peptide fragments[J]. Proceedings of the National Academy of Sciences of the United States of America, 1967, 58(3): 1235.<br>
+[6]de Prat Gay G, Ruiz-Sanz J, Davis B, et al. The structure of the transition state for the association of two fragments of the barley chymotrypsin inhibitor 2 to generate native-like protein: implications for mechanisms of protein folding[J]. Proceedings of the National Academy of Sciences, 1994, 91(23): 10943-10946.<br>
+[7]Bae C, Suchyna T M, Ziegler L, et al. Human PIEZO1 ion channel functions as a split protein[J]. PloS one, 2016, 11(3): e0151289.<br>
+[8]Johnsson N, Varshavsky A. Split ubiquitin as a sensor of protein interactions in vivo[J]. Proceedings of the National Academy of Sciences, 1994, 91(22): 10340-10344.<br>
+[9]Pelletier J N, Campbell-Valois F X, Michnick S W. Oligomerization domain-directed reassembly of active dihydrofolate reductase from rationally designed fragments[J]. Proceedings of the National Academy of Sciences, 1998, 95(21): 12141-12146.<br>
+[10]Galarneau A, Primeau M, Trudeau L E, et al. β-Lactamase protein fragment complementation assays as in vivo and in vitro sensors of protein–protein interactions[J]. Nature biotechnology, 2002, 20(6): 619.<br>
+[11]Massoud T F, Paulmurugan R, Gambhir S S. A molecularly engineered split reporter for imaging protein-protein interactions with positron emission tomography[J]. Nature medicine, 2010, 16(8): 921.<br>
+[12]Martell J D, Yamagata M, Deerinck T J, et al. A split horseradish peroxidase for the detection of intercellular protein–protein interactions and sensitive visualization of synapses[J]. Nature biotechnology, 2016, 34(7): 774.<br>
+[13]Kaihara A, Kawai Y, Sato M, et al. Locating a Protein− Protein Interaction in Living Cells via Split Renilla Luciferase Complementation[J]. Analytical chemistry, 2003, 75(16): 4176-4181.<br>
+[14]Paulmurugan R, Gambhir S S. Novel fusion protein approach for efficient high-throughput screening of small molecule–mediating protein-protein interactions in cells and living animals[J]. Cancer research, 2005, 65(16): 7413-7420.<br>
+[15]Kim S B, Kanno A, Ozawa T, et al. Nongenomic activity of ligands in the association of androgen receptor with SRC[J]. ACS chemical biology, 2007, 2(7): 484-492.<br>
+[16]Luker K E, Smith M C P, Luker G D, et al. Kinetics of regulated protein–protein interactions revealed with firefly luciferase complementation imaging in cells and living animals[J]. Proceedings of the National Academy of Sciences, 2004, 101(33): 12288-12293.<br>
+[17]Paulmurugan R, Umezawa Y, Gambhir S S. Noninvasive imaging of protein–protein interactions in living subjects by using reporter protein complementation and reconstitution strategies[J]. Proceedings of the National Academy of Sciences, 2002, 99(24): 15608-15613.<br>
+[18]Wehr M C, Laage R, Bolz U, et al. Monitoring regulated protein-protein interactions using split TEV[J]. Nature methods, 2006, 3(12): 985.</p>
+<h1 style="font-family: 'Times New Roman' !important; "><a name="Microalgae movement model">2  Microalgae movement model</a></h1>
+    <div  class="cart">
+      <h1  style="font-family: 'Times New Roman' !important; "><a name="Inspiration" >Introduction </a><img src="https://static.igem.org/mediawiki/2019/a/ab/T--DUT_China_B--Modelnote.svg" class="icon"> </h1>
-			<p style="font-family: 'Times New Roman' !important;   ">The combination of N-hrluc-PhyB and C-hrluc-PIF3 stimulated by red light can produce blue light and cause microalgae to move.Therefore, the degree of combination of N-hrluc-PhyB and C-hrluc-PIF3 directly determines whether our design is reasonable and the experiment is successful. Based on the simplest chemical thermodynamic formula, the relationship between binding rate and concentration and affinity constant was established in two steps.In addition, as a commonly used photopolymeric protein actuator in photogenetics, our work can not only guide our own experiments, but also provide a simple quantitative model for this field.</p>
+      <p style="font-family: 'Times New Roman' !important;   ">Our project aims to take advantage of the phototaxis characteristics of microalgae to construct a drug carrier which can be directed to the diseased cells under the guidance of light. We first analyzed the movement of microalgae in water. Using Langevin equation to describe the movement of microalgae in water, and then according to the observation of the pause of microalgae movement (<a href="https://static.igem.org/mediawiki/2019/0/0c/T--DUT_China_B--microalgae_breaststroke.mp4">movie</a>). We called it "breaststroke", and established the microalgae "breaststroke" model. According to Langevin equation, the viscosity in water can be replaced by the viscosity in blood to imitate the movement in blood. In addition, since red blood cells account for 45% of total volume of blood, the collision between microalgae and red blood cells will change the energy of microalgae. Since the collision between red blood cells and microalgae is random and complex, we assume that the red blood cells and microalgae are spheres. Under this assumption, the energy loss distribution shall conform the three-dimensional trigonometric function, so we established the microalgae movement model. Our model is the first to systematically describe the movement of microalgae in the water and blood. It has great originality and application value.</p>
-                </div>
+    </div>
+    <div class="cart">
+                <h2 style="font-family: 'Times New Roman' !important; ">Ⅰ.The movement of microalgae in water</h2>
+               <h2 style="font-family: 'Times New Roman' !important; "> Design<img src="https://static.igem.org/mediawiki/2019/b/b7/T--DUT_China_B--Modeldata.svg" class="icon"> </h2>
+<h3>Parameters</h3>
+         <div style="text-align:center;">     <img src="https://static.igem.org/mediawiki/2019/d/df/T--DUT_China_B--MD3.png"> </div>
+<p>According to Langevin equation, the motion of microalgae in water can be described as（Fig.1）：</p>
+<img src="https://static.igem.org/mediawiki/2019/6/68/T--DUT_China_B--M3F1.jpg">
+    <p>By formula transformation, (1) can be simplified as:</p>
+<img src="https://static.igem.org/mediawiki/2019/b/b2/T--DUT_China_B--M3F2.jpg">
+ <p>Defined：</p>
+<img src="https://static.igem.org/mediawiki/2019/d/d2/T--DUT_China_B--insectmove.png" >
+<p style="text-align:center;">Fig.1 The movement of microalgae in water</p>
+<p>According to the phenomenon that we observed in the experiment: microalgae do not move at a constant speed, it pauses with periodicity during the movement(<a href="https://static.igem.org/mediawiki/2019/0/0c/T--DUT_China_B--microalgae_breaststroke.mp4">movie</a>). We assume that the flagellum of microalgae generates a huge thrust, which can make microalgae generate an initial velocity in a very short time.
+The motion of microalgae in water</p>
+<img src="https://static.igem.org/mediawiki/2019/a/a9/T--DUT_China_B--M3F4.jpg" >
+<p>Simultaneous equations:</p>
+<img src="https://static.igem.org/mediawiki/2019/8/87/T--DUT_China_B--M3F5.jpg"/>
+<div style="text-align:center;"> <img  src="https://static.igem.org/mediawiki/2019/4/4f/T--DUT_China_B--M3F6.jpg "></div>
+<p>Therefore the motion curve of the microalgae in the water can be determined.</p>
+<h3>Caculation</h3>
+<p>(1)Solve : The average velocity  of microalgae can be obtained from the experiment, then the distance in one cycle can be obtained:</p>
+<img src="https://static.igem.org/mediawiki/2019/2/2a/T--DUT_China_B--M3F7.jpg">
+<p>(2)Solve the expression of :<i>v<sup>0</sup></i>= 1.96mm/s.</p>
+<img src="https://static.igem.org/mediawiki/2019/5/5c/T--DUT_China_B--M3F8.jpg">
+<div style="text-align:center;"> <img src="https://static.igem.org/mediawiki/2019/d/df/T--DUT_China_B--M3-2.jpg"></div>
+<p style="text-align:center;">Fig.2 Periodic curve of microalgae movement in water</p>
+<h2> Ⅱ. The movement of microalgae in the blood </h2>
+<h3>Parameters</h3>
+<div style="text-align:center;"> <img src="https://static.igem.org/mediawiki/2019/a/a0/T--DUT_China_B--MD32.png"></div>
+<h3>Design</h3>
+<p>In the blood, since the speed of movement of microalgae is related to the light intensity, we can control the light intensity to control the speed of movement of microalgae. Through experiments, we measured the average speeds of microalgae movement under different light conditions in the water（Fig.3）：</p>
+<div style="text-align:center;"> <img src="https://static.igem.org/mediawiki/2019/4/46/T--DUT_China_B--M3-3.jpg"></div>
+<p style="text-align:center;">Fig.3 The relationship between the microalgae movement speeds and light intensities.</p>
+<p>The equation obtained through experimental fitting:</p>
+<img src="https://static.igem.org/mediawiki/2019/a/a8/T--DUT_China_B--M3F9.jpg">
+<p>The energy of microalgae in the blood without collision:</p>
+<img src="https://static.igem.org/mediawiki/2019/e/ea/T--DUT_China_B--M3F10.jpg">
+<p>Following we will consider the collision in the blood (Fig.4):</p>
+<div style="text-align:center;"> <img src="https://static.igem.org/mediawiki/2019/a/ac/T--DUT_China_B--modelfig4.png"></div>
+<p style="text-align:center;">Fig.4 Microalgae as a drug carrier in the blood against the flow of blood.</p>
+<p>As the figure above, if microalgae swim upstream in the blood, they will be collided by red blood cells when they move in the blood. So, the speed of microalgae movement may change at any time. We find that the percentage of red blood cells in the blood vessels is about 45%. Which means if we take a cross-section of the blood vessels, the percentage of red blood cells is 45%(Fig.5).</p>
+<div style="text-align:center;"> <img src="https://static.igem.org/mediawiki/2019/6/60/T--DUT_China_B--M3-5.jpg"></div>
+<p style="text-align:center;">Fig.5 cross-section of the blood vessels</p>
+<p>Since red blood cells are randomly distributed in blood vessels, we assume that they are also randomly distributed in every small area of every cross section.</p>
+<img  src="https://static.igem.org/mediawiki/2019/7/71/T--DUT_China_B--M3F11.jpg" >
+<p>Thus, we equivalent the number of red blood cells colliding one microalgae cross section:</p>
+<img src="https://static.igem.org/mediawiki/2019/2/2b/T--DUT_China_B--M3F12.jpg">
+<p>The worst case scenario is a head-on collision between the microalgae and the red blood cells. In this case, microalgae lost the most energy. Assume the red blood cell energy is 0 after the collision. By energy conservation, microalgae lost the same energy with red blood cells. While, during the movement of microalgae, the collision angle between microalgae and red blood cells is random. Therefore, we use three-dimensional trigonometric functions to describe the energy loss of microalgae and red blood cells(Fig.6).</p>
+<div style="text-align:center;"> <img src="https://static.igem.org/mediawiki/2019/0/0a/T--DUT_China_B--M3-6.jpg"></div>
+<p style="text-align:center;">Fig.6 The energy loss of microalgae</p>
+<p>Through the operation of double integral, we obtained the energy lost of each section:</p>
+<img src="https://static.igem.org/mediawiki/2019/1/11/T--DUT_China_B--M3F14.jpg">
+<p>The energy of microalgae after collision in the blood:</p>
+<img src="https://static.igem.org/mediawiki/2019/8/89/T--DUT_China_B--M3F15.jpg">
+<p>If the microalgae need to swim upstream in the blood, then:</p>
+<img src="https://static.igem.org/mediawiki/2019/4/4a/T--DUT_China_B--M3F16.jpg" >
+<h3>Calculation</h3>
+<p>(1)Calculate N:
+By taking the radius of red blood cells and microalgae, the following equation can be obtained:</p>
+<img src="https://static.igem.org/mediawiki/2019/a/ad/T--DUT_China_B--M3F17.jpg">
+<p>(2)Calculate Q:
+The energy they lost in unit time of each blood vessel section:</p>
+<img  src="https://static.igem.org/mediawiki/2019/6/63/T--DUT_China_B--M3F18.jpg">
+<p>(3)We can obtain the relationship between the light intensity and different blood flow velocity we set（Fig.7）:</p>
+<img srrc="https://static.igem.org/mediawiki/2019/b/b1/T--DUT_China_B--M3F19.jpg">
+<img src="https://static.igem.org/mediawiki/2019/e/e5/T--DUT_China_B--M3F20.jpg">
+<div style="text-align:center;width:80%;"> <img src="https://static.igem.org/mediawiki/2019/2/29/T--DUT_China_B--M3-7.jpg" ></div>
+<p style="text-align:center;">Fig.7 Luminous intensity threshold values that allows microalgae to swim upstream in the blood at different blood flow velocities.</p>
+<h3>Conclusion</h3>
+<p>(1) We are the first to describe the microalgae periodic motion phenomenon in water using a precise mathematical expression.</p>
+<p>(2) We obtained the minimum light intensity thresholds required to control microalgae to swim as intended at different blood flow velocities. In the future, microalgae may act as a drug carrier. Our model may provide some valuable data and rules for clinical trials.</p></div>
+<h1 style="font-family: 'Times New Roman' !important; "><a name="Protein-interact model">3 Protein-protein interaction model</a></h1>
+    <div  class="cart">
+      <h1  style="font-family: 'Times New Roman' !important; "><a name="Inspiration" >Introduction </a><img src="https://static.igem.org/mediawiki/2019/a/ab/T--DUT_China_B--Modelnote.svg" class="icon"> </h1>
+      <p style="font-family: 'Times New Roman' !important;   ">Protein interaction is widespread in all live units, and its importance is self-evident. In our design, red light can stimulate the combination of N-hrluc-PhyB and C-hrluc-PIF3 and then produce blue light to guide microalgae to move(Fig.1).Therefore, whether our design is reasonable and the experiment successes depends on the degree of combination of N-hrluc-PhyB and C-hrluc-PIF3. Based on chemical thermodynamic formula, the relationship between binding rate and concentration and affinity constant is established in two steps. In addition, as a commonly used light-controlled polymeric protein actuator, this model can not only guide our own experiment, but also provide a quantitative model for the field of optogenetics.</p>
+      <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/5/54/T--DUT_China_B--modelf1.jpg" style="display: inline-block;width:80%;" />
+             <center> <br> <p style="text-align:center;position:relative;">Fig.1  A: The combination of N-hrluc and C-hrluc without the assistance by optical protein switch(PhyB&PIF3).  B: The combination of N-hrluc and C-hrluc with the assistance by optical protein switch.</p>  </center>
+      </div>
+    </div>
-               <div class="cart">
+    <div class="cart">
-               <h1 style="font-family: 'Times New Roman' !important; ">Background <img src="https://static.igem.org/mediawiki/2019/b/b7/T--DUT_China_B--difficultities.svg" class="icon"> </h1>
+              <h1 style="font-family: 'Times New Roman' !important; ">Design<img src="https://static.igem.org/mediawiki/2019/b/b7/T--DUT_China_B--Modeldata.svg" class="icon"> </h1>
-               <p>The molecular reaction in the solution can be divided into two steps:the first step is red light stimulation of PhyB and PIF3 binding, the second step is that the binding of PhyB and PIF3 facilitates the binding of two parts(NHrluc、Chrluc) of the Sea kidney luciferase. The molecular reactions carried out in our model and solution are also divided into two steps. </P>
+               <p>The molecular reaction in the solution can be described that the binding of PhyB and PIF3 stimulated by red light facilitates the binding of the two parts(N-Hrluc, C-hrluc) of the Renilla luciferase. Therefore this molecular reaction in solution can be divided into two stages:</P>
-                <P>The first part: only the binding of PhyB protein and PIF6 protein is considered.</p>
+              <h2>StageⅠ: only the binding of PhyB and PIF3 is considered.</h2>
-		<p>In order to find out the effect of the initial concentration of PhyB protein and PIF3 protein on the binding rate of the two proteins, we considered the binding of only PhyB and PIF3 proteins, and through  of  two kinds of proteins .Then,we get a specific equation between the initial concentration of protein and the binding rate of the two proteins.</p>
+              <p>In order to find out the effect of the initial concentration of PhyB and PIF3 on the binding rate of the two proteins, we get a specific equation between the initial concentration of protein and the binding rate of two proteins through dissociation constant Kd.</p>
-                  <table class="table-fill">
+ <h3>Parameters</h3>
+<table class="table-fill">
       <thead>
       <tr>
@@ Line 290: / Line 548: @@
      </tbody>
       </table>
-		</div>
+   <h3>Hypothesis: </h3><br>
+     <p>After the binding of two proteins, the solution contains three substances: protein PhyB, protein PIF3 and the combination of them. That is, the material is conserved and there is no protein loss during the process of the reaction.</p>
+    <h3>Calculation</h3><br>
-		<div class="cart" >
+                 <p>Considering the binding of PhyB and PIF3 in solution：</p>
-               <h1 style="font-family: 'JosefinSans-Light' !important; ">Chlamydomonas reinhardtii<img src="https://static.igem.org/mediawiki/2019/a/a2/T--DUT_China_B--wei.svg" class="icon"> </h1>
-    <h3>Hypothesis: </h3>
-     <p>(1)&nbsp;After the binding of the two proteins, the solution contains three substances: protein PhyB, protein PIF6 and the
-    combination of the two. That is, the material is conserved and there is no protein loss before and after the reaction.
-     <br>
-    (2) &nbsp; Except for the initial protein concentration, other factors have the same effect on the binding of the two proteins.</p>
-    <h3>Modeling</h3>
-                 <p>Considering the binding of PhyB protein and PIF6 protein in solution, there are/p>
              <div style="text-align: center; width: 100%; ">
-                 <img alt="" src="https://static.igem.org/mediawiki/2019/7/78/T--DUT_China_B--model1.jpg" style="display: inline-block;width:80%;" />
+                 <img alt="" src="https://static.igem.org/mediawiki/2019/a/aa/T--DUT_China_B--modelf11.jpg" style="display: inline-block;width:80%;" />
               <center> <br> <p style="left:45%;position:relative;"> </p>  </center>
                  </div>
-                 <p>If the free concentration of PhyB is A, the free concentration of PIF6 is B, and the concentration of PIF6 is AB, the binding constant is:</P>
+                 <p>If the free concentration of PhyB is A, the free concentration of PIF3 is B, and the concentration of PhyB&PIF3 is AB, the binding constant :</P>
                  <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/b/ba/T--DUT_China_B--model2.jpg" style="display: inline-block;width:80%;" />
               <center> <br> <p style="left:45%;position:relative;"> </p>  </center>
                  </div>
-                 <p>The dissociation constant is</P>
+                 <p>The dissociation constant :</P>
                  <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/c/ca/T--DUT_China_B--model3.jpg" style="display: inline-block;width:80%;" />
               <center> <br> <p style="left:45%;position:relative;"> </p>  </center>
                  </div>
-                <p>Then the expression of the binding rate is</p>
+                <p>Then the expression of the binding rate:</p>
                 <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/e/e0/T--DUT_China_B--model4.png" style="display: inline-block;width:80%;" />
               <center> <br> <p style="left:45%;position:relative;"> </p>  </center>
                  </div>
-                <p>Available jointly by (2) (3) (4):</p>
+                <p>By the (2) (3) (4) simultaneous:</p>
                 <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/1/14/T--DUT_China_B--model5.png" style="display: inline-block;width:80%;" />
               <center> <br> <p style="left:45%;position:relative;"> </p>  </center>
                  </div>
-                <p>Because we assume that the initial concentration of PhyB protein is higher than that of PIF6 protein, that is,</p>
+                <p>We assume that the initial concentration of PhyB protein is higher than that of PIF3 protein, that is,</p>
+                  <div style="text-align: center; width: 100%; ">
+                  <img alt="" src="https://static.igem.org/mediawiki/2019/c/c2/T--DUT_China_B--model6.png" style="display: inline-block;width:80%;" />
+                </div>
                  <p>Available jointly by (5)(6):</p>
                    <div style="text-align: center; width: 100%; ">
-                 <img alt="" src="https://static.igem.org/mediawiki/2019/c/c2/T--DUT_China_B--model6.png" style="display: inline-block;width:80%;" />
+                 <img alt="" src="https://static.igem.org/mediawiki/2019/2/22/T--DUT_China_B--model7.png" style="display: inline-block;width:80%;" />
               <center> <br> <p style="left:45%;position:relative;"> </p>  </center>
                  </div>
-               <p>The specific equation between the initial concentration of the two proteins and the binding rate of the two proteins was obtained.Available from (7)</p>
+              <p>The specific equation between the initial concentration of the two proteins and the binding rate of the two proteins was obtained.<br>
+              Available from (7):
+              </p>
                 <div style="text-align: center; width: 100%; ">
-                 <img alt="" src="https://static.igem.org/mediawiki/2019/2/22/T--DUT_China_B--model7.png" style="display: inline-block;width:80%;" />
+                 <img alt="" src="https://static.igem.org/mediawiki/2019/c/ca/T--DUT_China_B--equation12.png" style="display: inline-block;width:80%;" />
               <center> <br> <p style="left:45%;position:relative;"> </p>  </center>
                  </div>
-                <p>When the light with the wavelength of 650 nm and a frequency of 20 Hz, we know Kd=500nM [1]（fig.1）：</p>
+                <p>When the light with the wavelength of 650 nm and a frequency of 20 Hz, we know Kd=500nM [1]（Fig.1）：</p>
   <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/4/4a/T--DUT_China_B--M1.jpg" style="display: inline-block;width:80%;" />
-              <center> <br> <p style="text-align:center;position:relative;"> Fig.1 Binding Affinity(Kd) of PIF3-PhyB</p>  </center>
+              <center> <br> <p style="text-align:center;position:relative;"> Fig.2 Binding Affinity(Kd ) of PIF3&PhyB</p>  </center>
                  </div>
-                <p>The main results are as follows: (1) similar to the barrel effect, the binding rate depends on the low concentration of protein, and it can be estimated from the diagram that the minimum threshold of the concentration of the two proteins is  mol /L at the same time.</p>
+              <h3>Calculation</h3><br>
+                <p>(1) Similar to the barrel effect, the binding rate depends on the low concentration of protein, and it can be estimated from the diagram that the minimum threshold of the concentration of the two proteins is "1" "0" ^"-6"  mol /L when the concentration of the two proteins is equal(Fig.3 and Fig.4).</p>
                  <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/6/6f/T--DUT_China_B--M2.png" style="display: inline-block;width:80%;" />
-                   <center> <br> <p style="text-align:center;position:relative;"> Fig.2 The binding rate of PhyB and PIF3 is related to the
+                   <center> <br> <p style="text-align:center;position:relative;">Fig.3  The binding rate of PhyB and PIF3 related to the concentration.</p>  </center>
-                      concentration of PhyB and PIF3</p>  </center>
                 </div>
                  <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/f/f3/T--DUT_China_B--M3.png" style="display: inline-block;width:80%;" />
-                   <center> <br> <p style="text-align:center;position:relative;"> Fig.3 The trend of binding rate when PhyB and PIF3 are at the same concentration.</p>  </center>
+                   <center> <br> <p style="text-align:center;position:relative;"> Fig.4 The trend of binding rate when PhyB and PIF3 are at the same concentration </p>  </center>
                  </div>
@@ Line 356: / Line 612: @@
                    <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/0/0c/T--DUT_China_B--model9.jpg" style="display: inline-block;width:80%;" />
-                   <center> <br> <p style="text-align:center;position:relative;"> </p>  </center>
+                   <center> <br> <p style="text-align:center;position:relative;"></p>  </center>
                  </div>
-                 <p>Under the excitation of red light,the value is drastically reduced due to the change of in the configuration of the protein ,only when,The combination rate in the dark can be ignored.</p>
+                 <p>Under the excitation of red light, the value K_d is drastically reduced due to the change of in the configuration of the protein, only whenK_d (dark)/K_d (hv)>10000, the combination rate in the dark can be ignored(Fig.5).</p>
                  <div style="text-align: center; width: 100%; ">
                  <img alt="" src="https://static.igem.org/mediawiki/2019/a/ac/T--DUT_China_B--M4.png" style="display: inline-block;width:80%;" />
-                  <center> <br> <p style="text-align:center;position:relative;">Fig.4 The net increase of the binding rate under red light eliminated the dark background effect. </p>  </center>
+                  <center> <br> <p style="text-align:center;position:relative;">Fig.5 The net increase of the binding rate under</p>
+                 <p style="text-align:center;position:relative;">red light eliminated the dark background effect.</p> </center>
                  </div>
-                   <h3>The second part: the combination of N-hrluc and C-hrluc.</h3>
+                   <h2>StageⅡ: the combination of N-hrluc and C-hrluc.</h2><br>
-<p>Under the action of "optical switch", the distance between N-hrluc and C-hrluc becomes closer, which is equivalent to increasing the concentration. The effect of "optical switch" on protein binding is explored by using the definition of equivalent concentration to replace the "optical switch".</p>
+<p>Under the action of "optical switch"（PhyB&PIF3）, the distance between N-hrluc and C-hrluc becomes closer, which is equivalent to the increasing of concentration. The effect of "optical switch" on protein binding is converted to value by equal.</p>
+<h3>Parameters</h3><br>
+   <table class="table-fill">
      <thead>
       <tr>
-       <th class="text-left">Symbol</th>
+       <th class="text-left">Parameter</th>
-        <th class="text-left">Meaning</th>
+        <th class="text-left">Description</th>
        </tr>
       </thead>
@@ Line 374: / Line 633: @@
       <tr>
       <td class="text-left">B<sub>i</sub>,i=1,2</td>
-      <td class="text-left">Represent N-hrluc and C-hrluc protein parts respectively.</td>
+      <td class="text-left"> N-hrluc and C-hrluc protein parts respectively</td>
        </tr>
       <tr>
        <td class="text-left">A<sub>i</sub>,i=1,2</td>
-       <td class="text-left">Represent PhyB and PIF3 respectively.</td>
+       <td class="text-left">PhyB and PIF3 respectively</td>
       </tr>
        <tr>
@@ Line 386: / Line 645: @@
      </tbody>
       </table>
+       <h3>Hypothesis: </h3>
+    <p>(1)The addition of optical switch has no effect on the conformational of hrluc.
+     <br>
+    (2)After the optical switch bound, the original protein is pulled closer, that is, the distance is shorten.</p>
+   <h3>Calculation</h3><br>
+   <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/1/1e/T--DUT_China_B--modelf6.jpg" style="display: inline-block;width:80%;" />
+                 <center> <br> <p style="text-align:center;position:relative;">Fig.6   A: The combination of N-hrluc and C-hrluc without the assistance by protein switch(PhyB&PIF3).  B: The combination of N-hrluc and C-hrluc with the assistance by protein switch</p>  </center>
+                </div>
+   <p>As shown in the following figure(Fig.6), optical switches A_1,A_2are added to N-hrluc and C-hrluc proteins. Under the irradiation of light, the optical switches A_1,A_2begin to combine. Due to the hypothesis that conformational structure of B1 and B2 does not change, so the affinity constant is unchanged, while the distance of B1 and B2 is reduced, which is equivalent to the increase of concentration.</p>
+ <p>The concentration of a particle can be obtained as follows:</p>
+        <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/4/46/T--DUT_China_B--model10.jpg" style="display: inline-block;width:80%;" />
+                 <center> <br> <p style="text-align:center;position:relative;"> </p>  </center>
+                </div>
+      <p>The volume can be defined by the spacing between proteins, that is,</p>
+        <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/1/16/T--DUT_China_B--model11.jpg" style="display: inline-block;width:80%;" />
+                 <center> <br> <p style="text-align:center;position:relative;"> </p>  </center>
+          </div>
+        <p>So:</p>
+        <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/4/4f/T--DUT_China_B--model12.jpg" style="display: inline-block;width:80%;" />
+                 <center> <br> <p style="text-align:center;position:relative;"> </p>  </center>
+          </div>
+         <p>According to the size of the protein switch(PhyB&PIF3), in the least ideal case, the distance between N-hrluc and C-hrluc is 16 nm. Assuming that the spacing is 15 nm, 10 nm and 5 nm, the equivalent concentrations can be caculated as 3×"1" "0" ^"-3" mol/L, 1×"1" "0" ^"-2" mol/L, 8×"1" "0" ^"-2" mol/L according to formula (12).</p>
+        <div style="text-align: center; width: 100%; ">
+                <img alt="" src="https://static.igem.org/mediawiki/2019/4/4f/T--DUT_China_B--M6.jpg" style="display: inline-block;width:80%;" />
+                 <center> <br> <p style="text-align:center;position:relative;">Fig.7  N-Luc and C-Luc are catched by the optical protein
+switch in spatial scale which means a higher binding rate (than being dispersed in solution)
+</p>  </center>
+          </div>
+<p> As long as the affinity constant of the two parts of fluorescein enzyme is less than  ,it can be considered to be completely binding(Fig.7).</p>
+   <h3>Conclusion</h3>
+   <p>(1) The optical switch proteins(PhyB and PIF3) will bind only when the concentration of optical switch protein reach 10-6mol/L according to our model(Fig.4).</p>
+   <p>(2) The equivalent concentration of renin luciferase depends on the structure size of optical protein switch(Fig.7). </p>
+ </div>
+    <!-- 介绍 -->
-      </div>
        <hr>
+<h1 style="font-family: 'Times New Roman' !important; ">Reference</h1>
+  <p> [1]Levskaya A, Weiner O D, Lim W A, et al. Spatiotemporal control of cell signalling using a light-switchable protein interaction[J]. Nature, 2009, 461(7266): 997. </p>
-	</div>
+  </div>
-	<script type="text/javascript">
+<script type="text/javascript">
-	window.onload=function(){
+  window.onload=function(){
-   		map_height=document.documentElement.clientHeight;//获取页面可见高度
+      map_height=document.documentElement.clientHeight;//获取页面可见高度
-   		document.getElementById("sides").style.height=map_height+"px";
+      document.getElementById("sides").style.height=map_height+"px";
-   	}
+    }
-   	function getScrollTop(){
+    function getScrollTop(){
-		var scrollTop=0;
+    var scrollTop=0;
-		if(document.documentElement.scrollTop){
+    if(document.documentElement.scrollTop){
-			scrollTop=document.documentElement.scrollTop;
+      scrollTop=document.documentElement.scrollTop;
-		}else if(document.body){
+    }else if(document.body){
-			scrollTop=document.body.scrollTop;
+      scrollTop=document.body.scrollTop;
-		}
+    }
-		return scrollTop;
+    return scrollTop;
-		}
+    }
-		document.onscroll = function(){
+    document.onscroll = function(){
-			console.log(getScrollTop())
+      console.log(getScrollTop())
-			if(getScrollTop()>(map_height+170)) {
+      if(getScrollTop()>(map_height+170)) {
-				document.getElementById("sides").style.position="fixed";
+        document.getElementById("sides").style.position="fixed";
-				document.getElementById("sides").style.top="0";
+        document.getElementById("sides").style.top="0";
-				document.getElementById("sides").style.left="0";
+        document.getElementById("sides").style.left="0";
-			}else{
+      }else{
-				document.getElementById("sides").style.position = 'static';
+        document.getElementById("sides").style.position = 'static';
-			}
+      }
-		}
+    }
-	</script>
+  </script>
 </body>
 </html>

Difference between revisions of "Team:DUT China B/Model"

Latest revision as of 03:15, 22 October 2019

1. Protein-split Model

Introduction

Design

Algorithm

Database

Result

Weight matrix:

w1

b1

w2

b2

Performance Evaluation Criteria

Conclusion

Experimental Verification

Reference

2 Microalgae movement model

Introduction

Ⅰ.The movement of microalgae in water

Design

Parameters

Caculation

Ⅱ. The movement of microalgae in the blood

Parameters

Design

Calculation

Conclusion

3 Protein-protein interaction model

Introduction

Design

StageⅠ: only the binding of PhyB and PIF3 is considered.

Parameters

Hypothesis:

Calculation

Calculation

StageⅡ: the combination of N-hrluc and C-hrluc.

Parameters

Hypothesis:

Calculation

Conclusion

Reference