17512580546 (Talk | contribs) |
|||
(61 intermediate revisions by 7 users not shown) | |||
Line 1: | Line 1: | ||
− | |||
{{NAU-CHINA}} | {{NAU-CHINA}} | ||
{{NAU-CHINA/header}} | {{NAU-CHINA/header}} | ||
Line 6: | Line 5: | ||
<meta charset="utf-8" /> | <meta charset="utf-8" /> | ||
<title>Safety</title> | <title>Safety</title> | ||
− | + | <script type="text/x-mathjax-config"> | |
− | + | MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}}); | |
+ | </script> | ||
+ | <script type="text/javascript" | ||
+ | src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"> | ||
+ | </script> | ||
+ | <script> | ||
+ | $(document).ready(function () { | ||
+ | $('#banner img').attr('src', 'https://static.igem.org/mediawiki/2018/1/1d/T--NAU-China--bannermodeldetails.png') | ||
+ | }); | ||
+ | </script> | ||
<style> | <style> | ||
− | + | .top-title { | |
− | . | + | color: #006030 !important; |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
} | } | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
.main-content h1 { | .main-content h1 { | ||
font-family: 'Avenir Next Condensed',sans-serif; | font-family: 'Avenir Next Condensed',sans-serif; | ||
− | color: # | + | color: #006030 !important; |
− | margin:32px 0 !important; | + | margin: 32px 0 !important; |
− | margin-left:40px !important; | + | margin-left: 40px !important; |
− | font-size: | + | font-size: 40px; |
− | font-weight:bold !important; | + | font-weight: bold !important; |
} | } | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
</style> | </style> | ||
</head> | </head> | ||
<body> | <body> | ||
− | + | <div class="topLine"> | |
− | < | + | <p class="top-title">Model</p> |
− | + | <p class="sec-title">Model Illustration</p> | |
− | + | ||
− | + | ||
− | + | ||
</div> | </div> | ||
− | + | <a href="https://2018.igem.org/Team:NAU-CHINA/Model"> | |
− | + | <img id="icon1" class="guide-icon" src="https://static.igem.org/mediawiki/2018/c/c3/T--NAU-China--overviewgreen.png" /> | |
− | + | </a> | |
+ | <a href="https://2018.igem.org/Team:NAU-CHINA/Auxiliary%20Understanding"> | ||
+ | <img id="icon2" class="guide-icon" src="https://static.igem.org/mediawiki/2018/c/c2/T--NAU-China--auxiliaryunderstandinggreen.png" /> | ||
+ | </a> | ||
+ | <div class="main-content"> | ||
<div class="textblock"> | <div class="textblock"> | ||
− | <h1>Part1 Extracellular signal | + | <p><a href="#Part1 Extracellular signal converting into intracellular signal">Part1 Extracellular signal converting into intracellular signal</a></p> |
+ | <p><a href="#Part2 Internal signal regulating recombinase production">Part2 Internal signal regulating recombinase production</a></p> | ||
+ | <p><a href="#Part3 Recombinase turns on downstream effect expression">Part3 Recombinase turns on downstream effect expression</a></p> | ||
+ | <p><a href="#Part4 Significance of the cell culture environment">Part4 Significance of the cell culture environment</a></p> | ||
+ | <p><a href="#Part5 Recombinase filtering and de-shocking">Part5 Recombinase filtering and de-shocking</a></p> | ||
+ | <p><a href="#Part6: Robustness of the model">Part6 Robustness of the model</a></p> | ||
+ | <h1 id="Part1 Extracellular signal converting into intracellular signal">Part1 Extracellular signal converting into intracellular signal</h1> | ||
<p>Tips: Skip part1 if the gene circuit is directly regulated by intracellular signal.</p> | <p>Tips: Skip part1 if the gene circuit is directly regulated by intracellular signal.</p> | ||
− | <p> | + | <p>SynNotch can convert an extracellular signal into an intracellular signal . Here, we take synNotch-TEV as an example. We use the number of TEV as the intracellular signal after conversion. TEV is an enzyme which is able to release the repressed promoter.</p> |
<p>In this section, 4 reactions are considered:</p> | <p>In this section, 4 reactions are considered:</p> | ||
− | <p>$$ | + | <p>$$G_{synNotch-TEV}\xrightarrow[\quad\quad]{C11}{synNotch-TEV}\ \ \ \ \ \ \ \ \ \ \ (1.1)$$ </p> |
− | <p>(1.1) represents the expression of synNotch which contains the TEV. Although the actual expression progress of the protein involves many chemical reactions from transcription to modification, these intermediate reactions satisfy the | + | |
− | <p>a)The factors affecting these reactions in the system are | + | <p>$$synNotch-TEV+Ts\xrightarrow[\quad\quad]{C12}U+TEV\ \ \ \ \ \ \ \ \ \ \ (1.2)$$</p> |
+ | |||
+ | <p>$$synNotch-TEV\xrightarrow[\quad\quad]{C13}\emptyset\ \ \ \ \ \ \ \ \ \ \ (1.3)$$</p> | ||
+ | |||
+ | <p>$$TEV\xrightarrow[\quad\quad]{C14}\emptyset\ \ \ \ \ \ \ \ \ \ \ (1.4)$$</p> <p>The specific contents are as follows:</p> | ||
+ | <p>$$G_{synNotch-TEV}\xrightarrow[\quad\quad]{C11}{synNotch-TEV}\ \ \ \ \ \ \ \ \ \ \ (1.1)$$</p> | ||
+ | |||
+ | <p>(1.1) represents the expression of synNotch which contains the TEV. Although the actual expression progress of the protein involves many chemical reactions from transcription to modification, these intermediate reactions satisfy the following conditions:</p> | ||
+ | <p>a)The factors affecting these reactions in the system are stable, such as the amount of amino acids, nucleotides, and helicases;</p> | ||
<p>b)The gene is regulated by a constitutive promoter and its expression is not regulated by other products of the system;</p> | <p>b)The gene is regulated by a constitutive promoter and its expression is not regulated by other products of the system;</p> | ||
− | <p>c)The product of | + | <p>c)The product of an intermediate reaction has no effects on the system functions.</p> |
− | <p> | + | <p> |
− | In the same way, the expression and degradation of other proteins involved in our model all suit the conditions of the refinement, so they can be considered as one single chemical reaction.</p> | + | We call these three conditions the post-integration conditions [1], and call these reactions the post-integration reaction. |
+ | In the same way, the expression and degradation of other proteins involved in our model all suit the conditions of the refinement, so they can be considered as one single chemical reaction. | ||
+ | </p> | ||
<p>Obviously (1.1) is a post-integration reaction, and the reaction rate is a constant, for the controlling conditions are unchanged.</p> | <p>Obviously (1.1) is a post-integration reaction, and the reaction rate is a constant, for the controlling conditions are unchanged.</p> | ||
− | <p>$$synNotch | + | <p>$$synNotch-TEV+Ts\xrightarrow[\quad\quad]{C12}U+TEV\ \ \ \ \ \ \ \ \ \ \ (1.2)$$</p> |
− | </p> | + | |
<p>(1.2) represents that the TEV protein is released by the intracellular domain of the sporadic synNotch when the synNotch-TEV protein was stimulated by an external Target signal. The rate of the reaction is calculated by the law of mass action which is in proportion to the number of receptors on the cell membrane and the number of external target signals.</p> | <p>(1.2) represents that the TEV protein is released by the intracellular domain of the sporadic synNotch when the synNotch-TEV protein was stimulated by an external Target signal. The rate of the reaction is calculated by the law of mass action which is in proportion to the number of receptors on the cell membrane and the number of external target signals.</p> | ||
− | <p>$$synNotch\ | + | <p>$$synNotch-TEV\xrightarrow[\quad\quad]{C13}\emptyset\ \ \ \ \ \ \ \ \ \ \ (1.3)$$</p> <p>(1.3) represents the degradation of synNotch membrane proteins linked to TEV proteins The degradation rate is directly proportional to its protein content.</p> |
− | + | <p>$$TEV\xrightarrow[\quad\quad]{C14}\emptyset\ \ \ \ \ \ \ \ \ \ \ (1.4)$$</p> <p>(1.4) represents the degradation of TEV protein.</p> | |
− | <p>$$TEV\ | + | |
− | + | ||
− | + | ||
<p>The degradation rate is directly proportional to its protein content.</p> | <p>The degradation rate is directly proportional to its protein content.</p> | ||
<p>The figures below show the result of simulation using the Gillespie algorithm.</p> | <p>The figures below show the result of simulation using the Gillespie algorithm.</p> | ||
<figure> | <figure> | ||
− | + | <img src="https://static.igem.org/mediawiki/2018/2/25/T--NAU-China--model1.svg"> | |
− | + | <figcaption><b>Fig. 1.1. </b>It shows the number of synNotch-TEV protein changes at different time. We can see that the synNotch-TEV of a single cell will be stable near 20,000 in the absence of a signal (Ts = 0).Because the express rate of synNotch-TEV is certain while the degradation rate of it is directly proportional to its protein content, with the rise of the amount of protein, its express rate will be equal to its degradation rate, making the amount of synNotch-TEV protein keeps stable near 20000.</figcaption> | |
− | + | </figure> | |
<p>It is concluded that the initial value of synNotch-TEV is stable at about 20,000 before the Target signal is received, which is used as the initial condition for subsequent translation.</p> | <p>It is concluded that the initial value of synNotch-TEV is stable at about 20,000 before the Target signal is received, which is used as the initial condition for subsequent translation.</p> | ||
<figure> | <figure> | ||
− | + | <img src="https://static.igem.org/mediawiki/2018/2/21/T--NAU-China--model2.svg"> | |
− | + | <figcaption><b>Fig. 1.2. </b>Relationship between the amount of synNotch-TEV and Ts number level in stable cell lines.</figcaption> | |
− | + | </figure> | |
− | + | <figure> | |
− | + | <img src="https://static.igem.org/mediawiki/2018/3/31/T--NAU-China--model3.svg"> | |
− | + | <figcaption><b>Fig. 1.3. </b> Relationship between the amount of TEV and Ts number level in stable cell lines.</figcaption> | |
− | + | <figcaption><b>Fig. 1.2</b> & <b>Fig. 1.3</b> When stimulated by Target signal, synNotch-TEV is also consumed by the reaction of releasing TEV. Because there are many synNotch-TEV proteins on the membrane, the rate of consumption is higher than the generation and the amount of synNotch-TEV will decrease, which causes the rate of degradation and TEV release to become slower. When the rate of consumption equals to that of generation, the amount of synNotch and the number of intracellular TEV will become stable.</figcaption> | |
− | + | </figure> | |
− | <p><b>Conclusion: </b>When a stabilized cell contacts different | + | <p><b>Conclusion: </b>When a stabilized cell contacts different levels of extracellular signals, synNotch membrane protein will release different amounts of TEV. The amount of TEV protein is in proportion of the strength of target signals.</p> |
− | + | ||
</div> | </div> | ||
</div> | </div> | ||
<div class="main-content"> | <div class="main-content"> | ||
<div class="textblock"> | <div class="textblock"> | ||
− | <h1>Part2 Internal signal regulating recombinase production</h1> | + | <h1 id="Part2 Internal signal regulating recombinase production">Part2 Internal signal regulating recombinase production</h1> |
<p> | <p> | ||
− | + | In this part, intracellular proteins can degrade or destroy the repressors, reducing their number. Different levels of protein have different effects on the relieving repression. Then, using Hill equation, we can calculate the expression rate of the recombinase based on the repressor protein number. At the same time, we use the error-reducing reaction to solve the problem of the extra recombinase. Here, we choose TEV as intracellular protein, and TetR as repressor protein. | |
</p> | </p> | ||
− | |||
<p>In this section, 12 reactions are considered:</p> | <p>In this section, 12 reactions are considered:</p> | ||
− | <p>$$ | + | <p>$$G_{TetR}\xrightarrow[\quad\quad]{C21}TetR\ \ \ \ \ \ \ \ \ \ \ (2.1)$$</p> |
+ | |||
+ | <p>$$TetR\xrightarrow[\quad\quad]{C22}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.2)$$</p> | ||
+ | |||
+ | <p>$$2TetR+tetO\xrightarrow[\quad\quad]{C23}TetOR\ \ \ \ \ \ \ \ \ \ \ (2.3)$$</p> | ||
+ | |||
+ | <p>$$TEV+TetOR\xrightarrow[\quad\quad]{C24}tetO+TEV+U\ \ \ \ \ \ \ \ \ \ \ (2.4)$$</p> | ||
+ | |||
+ | <p>$$TetOR\xrightarrow[\quad\quad]{C25}tetO+U\ \ \ \ \ \ \ \ \ \ \ (2.5)$$</p> | ||
+ | |||
+ | <p>$$G_{rec}\xrightarrow[\quad\quad\quad\quad\quad\quad]{C31-TetOR}rec\ \ \ \ \ \ \ \ \ \ \ (2.6)$$</p> | ||
+ | |||
+ | <p>$$rec\xrightarrow[\quad\quad]{C33}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.7)$$</p> | ||
+ | |||
+ | <p>$$G_{RDF-inhibitor}\xrightarrow[\quad\quad\quad\quad\quad\quad]{C32-TetOR}RDF-inhibitor\ \ \ \ \ \ \ \ \ \ \ (2.8)$$</p> | ||
+ | |||
+ | <p>$$RDF-inhibitor\xrightarrow[\quad\quad]{C34}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.9)$$</p> | ||
+ | |||
+ | <p>$$G_{rec-inhibitor}\xrightarrow[\quad\quad]{C35}rec-inhibitor\ \ \ \ \ \ \ \ \ \ \ (2.10)$$</p> | ||
+ | |||
+ | <p>$$rec-inhibitor\xrightarrow[\quad\quad]{C36}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.11)$$</p> | ||
+ | |||
+ | <p>$$rec+rec-inhibitor\xrightarrow[\quad\quad]{C37}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.12)$$</p> | ||
+ | <p>The specific content is as follows. (As for each chart in this part, we choose the quantity of substance as the vertical axis, and 0.01 hour as the unit of the horizontal axis. Different subplots correspond to different levels of intracellular protein number).</p> | ||
+ | <p>$$G_{TetR}\xrightarrow[\quad\quad]{C21}TetR\ \ \ \ \ \ \ \ \ \ \ (2.1)$$</p> | ||
<p>(2.1) indicates the expression of the TetR protein. Through the analysis of (1.1) in part1, it is easy to conclude that (2.1) is a post-integration reaction, and the rate of the reaction is constant.</p> | <p>(2.1) indicates the expression of the TetR protein. Through the analysis of (1.1) in part1, it is easy to conclude that (2.1) is a post-integration reaction, and the rate of the reaction is constant.</p> | ||
− | <p>$$TetR\ | + | <p>$$TetR\xrightarrow[\quad\quad]{C22}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.2)$$</p> |
− | <p>(2.2) indicates degradation of TetR protein. The degradation rate is directly proportional to its protein content.</p> | + | <p>(2.2) indicates the degradation of TetR protein. The degradation rate is directly proportional to its protein content.</p> |
− | <p>$$ | + | <p>$$2TetR+tetO\xrightarrow[\quad\quad]{C23}TetOR\ \ \ \ \ \ \ \ \ \ \ (2.3)$$</p> |
<p>(2.3) represents the generation of the binary complex TetOR, and the reaction rate is calculated according to the law of mass action in the model.</p> | <p>(2.3) represents the generation of the binary complex TetOR, and the reaction rate is calculated according to the law of mass action in the model.</p> | ||
− | <p>$$TEV+ | + | <p>$$TEV+TetOR\xrightarrow[\quad\quad]{C24}tetO+TEV+U\ \ \ \ \ \ \ \ \ \ \ (2.4)$$</p> |
<p>(2.4) indicates that the TEV is a separate TetOR binary complex. This reaction rate is calculated according to the law of mass action.</p> | <p>(2.4) indicates that the TEV is a separate TetOR binary complex. This reaction rate is calculated according to the law of mass action.</p> | ||
− | <p>$$ | + | <p>$$TetOR\xrightarrow[\quad\quad]{C25}tetO+U\ \ \ \ \ \ \ \ \ \ \ (2.5)$$</p> |
<p>(2.5) indicates that a small amount of TetOR is self-separate. The reaction rate is calculated according to the law of mass action.</p> | <p>(2.5) indicates that a small amount of TetOR is self-separate. The reaction rate is calculated according to the law of mass action.</p> | ||
<p>The figure below shows the calculated results.</p> | <p>The figure below shows the calculated results.</p> | ||
<figure> | <figure> | ||
− | + | <img src="https://static.igem.org/mediawiki/2018/b/b6/T--NAU-China--model4.svg"> | |
− | + | <figcaption><b>Fig. 2.1. </b> After the cells undergo the above reaction, the corresponding number of the repressor is changed at the corresponding signal number (second signal conversion). We can see, when the signal strength is higher than the 9level, TetR is seldom detected. Therefore, we do not need to simulate higher signal number.</figcaption> | |
− | + | </figure> | |
<p>Using the above calculated amount of repressor, the Hill equation is used to determine the expression rate of the recombinase.</p> | <p>Using the above calculated amount of repressor, the Hill equation is used to determine the expression rate of the recombinase.</p> | ||
− | <p>$$ | + | <p>$$G_{rec}\xrightarrow[\quad\quad\quad\quad\quad\quad]{C31-TetOR}rec\ \ \ \ \ \ \ \ \ \ \ (2.6)$$</p> |
<p>(2.6) indicates the generation of rec (rec is the abbreviation of recombinase) , which is also a post-integration reaction, but the expression rate is regulated by the repressor TetR. The reaction rate in the model is calculated using the Hill equation, and the rate is ((N/k) ^n+1) ^-1. Accordingly, N is the number of TetR, and K is the half-maximum of the suppression effect.</p> | <p>(2.6) indicates the generation of rec (rec is the abbreviation of recombinase) , which is also a post-integration reaction, but the expression rate is regulated by the repressor TetR. The reaction rate in the model is calculated using the Hill equation, and the rate is ((N/k) ^n+1) ^-1. Accordingly, N is the number of TetR, and K is the half-maximum of the suppression effect.</p> | ||
− | <p>$$rec\ | + | <p>$$rec\xrightarrow[\quad\quad]{C33}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.7)$$</p> |
− | <p>(2.7) indicates degradation of rec. The degradation rate is directly proportional to rec | + | <p>(2.7) indicates degradation of rec. The degradation rate is directly proportional to rec number.</p> |
<figure> | <figure> | ||
− | + | <img src="https://static.igem.org/mediawiki/2018/f/f0/T--NAU-China--model5.svg"> | |
− | + | <figcaption><b>Fig. 2.2. </b> Different subplots simulate the different number of the repressor calculated from the intracellular signal producing the corresponding expression of rec over time, which is calculated by the Hill equation. When this figure was generated, the response of the minor rec-inhibitor elimination of rec leakage (2.10) - (2.12) was taken into account, making the simulation effect even better. The simulation results without adding this series of noise reduction reactions were shown in <b>Fig. 5.1. </b> </figcaption> | |
− | <p>$$ | + | </figure> |
− | </p> | + | <p>$$G_{RDF-inhibitor}\xrightarrow[\quad\quad\quad\quad\quad\quad]{C32-TetOR}RDF-inhibitor\ \ \ \ \ \ \ \ \ \ \ (2.8)$$</p> |
<p>(2.8) indicates the generation of RDF. Like rec, it is also a post-integration reaction, whose rate is regulated by the repressor TetR. The reaction rate in the model is calculated by the Hill equation, and the rate is ((N/k) ^n+1) ^-1. Similarly, N is the number of TetR, and K is the half-maximum of suppression effect.</p> | <p>(2.8) indicates the generation of RDF. Like rec, it is also a post-integration reaction, whose rate is regulated by the repressor TetR. The reaction rate in the model is calculated by the Hill equation, and the rate is ((N/k) ^n+1) ^-1. Similarly, N is the number of TetR, and K is the half-maximum of suppression effect.</p> | ||
− | <p>$$ | + | <p>$$RDF-inhibitor\xrightarrow[\quad\quad]{C34}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.9)$$</p> |
− | </p> | + | <p>Degradation of RDF-inhibitor is represented by (2.9). The degradation rate is directly proportional to the number of RDF-inhibitor. The roles and meanings of this reaction will be discussed in the next part.</p> |
− | <p>Degradation of RDF-inhibitor is represented by (2.9). The degradation rate is directly proportional to the | + | <p>$$G_{rec-inhibitor}\xrightarrow[\quad\quad]{C35}rec-inhibitor\ \ \ \ \ \ \ \ \ \ \ (2.10)$$</p> |
− | <p> | + | <p>(2.10) indicates the generation of the rec-inhibitor, which is also a post-integration reaction, and the expression rate is constant.</p> |
− | $$ | + | <p>$$rec-inhibitor\xrightarrow[\quad\quad]{C36}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.11)$$</p> |
− | </p> | + | <p>(2.11) indicates the degradation of the rec-inhibitor. The degradation rate is directly proportional to rec-inhibitor number.</p> |
− | + | <p>$$rec+rec-inhibitor\xrightarrow[\quad\quad]{C37}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.12)$$</p> | |
− | <p>$$rec\ | + | <p> (2.12) indicates that the rec-inhibitor combined with the leaked rec. It is used to filter rec. The reaction rate was calculated according to the law of mass action and enables the low expression of rec to be combined to ensure the stability of the system.</p> |
− | + | <figure> | |
− | <p>$$rec+rec\ | + | <img src="https://static.igem.org/mediawiki/2018/b/b8/T--NAU-China--model6.svg"> |
− | + | <figcaption><b>Fig. 2.3. </b> Different sub-figures correspond to different concentrations of TEV, corresponding to different concentrations of TetR. At different concentrations of TetR, the rec concentrations are different, which affect the rec-inhibitor number. At extremely low permeation, rec-inhibitor can almost completely eliminate the rec leakage from the TetR system. </figcaption> | |
− | + | </figure> | |
− | + | <p>It is <b>concluded</b> that internal signal TEV regulates recombinase production. According to (2.4), different concentrations of TEV correspond to different concentrations of TetR. By (2.6), the concentrations of TetR affect rec concentrations. Through the simulation of (2.12), a small amount of rec-inhibitor can almost completely eliminate the rec leakage from the TetR system. This pathway can significantly improves the system stability.</p> | |
− | + | ||
− | <p>It is <b>concluded</b> that internal signal TEV regulates recombinase production. According to (2.4), different concentrations of TEV correspond to different concentrations of TetR. By (2.6), the concentrations of TetR affect rec concentrations. Through the simulation of (2.12), a small amount of rec-inhibitor can almost completely eliminate the rec leakage from the TetR system. This pathway can significantly | + | |
</div> | </div> | ||
</div> | </div> | ||
<div class="main-content"> | <div class="main-content"> | ||
<div class="textblock"> | <div class="textblock"> | ||
− | <h1>Part3 Recombinase turns on downstream effect expression</h1> | + | <h1 id="Part3 Recombinase turns on downstream effect expression">Part3 Recombinase turns on downstream effect expression</h1> |
<p> | <p> | ||
− | + | In this part, we use the above calculated amount of the rec expression at different signal levels. Using rec and RDF, we can switch the state and calculate the expression of the downstream gene. In this section, 5 reactions are considered: | |
− | <p> | + | </p> |
− | $$rec+ | + | <p>$$rec+\overline{G}_{downstream}\xrightarrow[\quad\quad]{C41}rec+G_{downstream}\ \ \ \ \ \ \ \ \ \ \ (3.1)$$</p> |
− | </p> | + | |
+ | <p>$$RDF+G_{dowmstream}\xrightarrow[\quad\quad]{C61}RDF+\overline{G}_{downstream}\ \ \ \ \ \ \ \ \ \ \ (3.2)$$</p> | ||
+ | |||
+ | <p>$$G_{RDF}\xrightarrow[\quad\quad]{C}RDF\ \ \ \ \ \ \ \ \ \ \ (3.3)$$</p> | ||
+ | |||
+ | <p>$$RDF-inhibitor+RDF\xrightarrow[\quad\quad]{C34}\emptyset\ \ \ \ \ \ \ \ \ \ \ \ (3.4)$$</p> | ||
+ | |||
+ | <p>$$G_{gfp}\xrightarrow[\quad\quad]{C}GFP\ \ \ \ \ \ \ \ \ \ \ (3.5)$$</p> | ||
+ | <p>$$rec+\overline{G}_{downstream}\xrightarrow[\quad\quad]{C41}rec+G_{downstream}\ \ \ \ \ \ \ \ \ \ \ (3.1)$$</p> | ||
<p>(3.1) represents that recombinase specifically recognizes the DNA sequence between the two sites and reverses it, so that the downstream gene expression is reversed to an open state. (Downstream genes are GFP and RDF Gene)</p> | <p>(3.1) represents that recombinase specifically recognizes the DNA sequence between the two sites and reverses it, so that the downstream gene expression is reversed to an open state. (Downstream genes are GFP and RDF Gene)</p> | ||
− | <p> | + | <p>$$RDF+G_{dowmstream}\xrightarrow[\quad\quad]{C61}RDF+\overline{G}_{downstream}\ \ \ \ \ \ \ \ \ \ \ (3.2)$$</p> |
− | $$RDF+ | + | |
− | </p> | + | |
<p>(3.2) indicates that the reverse recombination factor specifically recognizes the two sites after the recombinase worked on the DNA, and makes the reversal DNA sequence between the sites restored to the shape.</p> | <p>(3.2) indicates that the reverse recombination factor specifically recognizes the two sites after the recombinase worked on the DNA, and makes the reversal DNA sequence between the sites restored to the shape.</p> | ||
− | <p>$$ | + | <p>$$G_{RDF}\xrightarrow[\quad\quad]{c}RDF\ \ \ \ \ \ \ \ \ \ \ (3.3)$$</p> |
− | </p> | + | <p>(3.3) indicates the expression of the reverse recombination factor, which exists only when the downstream gene is turned into an open state. The expression of the RDF process is not regulated, and also is a post-integration reaction, so the expression rate is constant.</p> |
− | <p>(3.3) indicates the expression of the reverse recombination factor, which exists only when the downstream gene is turned into an open state. The expression of the RDF process is not regulated, and also a post-integration reaction, so the expression rate is constant.</p> | + | <p>$$RDF-inhibitor+RDF\xrightarrow[\quad\quad]{C34}\emptyset\ \ \ \ \ \ \ \ \ \ \ \ (3.4)$$</p> |
− | <p>$$ | + | <p>(3.4) indicates that the RDF-inhibitor is combined with the intracellular RDF when the system requests downstream gene opening. The reaction rate is calculated according to the law of mass action. It makes our system stable in the open state. When it is not requested to open, the increase of RDF is not limited, so that it can function to shut down the system and realize real-time feedback. </p> |
− | </p> | + | |
− | <p>(3.4) indicates that the RDF-inhibitor is combined with the intracellular RDF when the system requests downstream opening. The reaction rate | + | |
− | + | ||
− | + | ||
− | + | ||
<figure> | <figure> | ||
− | + | <img src="https://static.igem.org/mediawiki/2018/7/75/T--NAU-China--model7.svg"> | |
− | + | <figcaption><b>Fig. 3.1. </b> Relationship between the amount of RDF-inhibitor and TEV number level in stable state.</figcaption> | |
− | <p>In this simulation, the off state is the initial state. When the TEV signal | + | </figure> |
− | <p>$$ | + | <figure> |
− | <p>The formula (3.5) indicates the production of the effectors. The downstream gene can be expressed when it is reversed to open. Here, GFP is taken as an example, and the expression is constant after the rectification reaction.</p> | + | <img src="https://static.igem.org/mediawiki/2018/d/d8/T--NAU-China--model8.svg"> |
+ | <figcaption><b>Fig. 3.2. </b> Relationship between the amount of RDF and TEV number level in stable state.</figcaption> | ||
+ | </figure> | ||
+ | <p>In this simulation, the off state is the initial state. When the TEV signal number level is less than 5, there is not enough rec to turn on the downstream gene,and RDF is not expressed; when the TEV signal number level is higher than 5, rec and RDF-inhibitor begin to express simultaneously after de-repression. Then,RDF and GFP are simultaneously expressed after the rec turns on the downstream gene. The combination of RDF-inhibitor and RDF ensures that there is not enough RDF in the cell to shut down the downstream gene, which makes the system stable.</p> | ||
+ | <p>$$G_{gfp}\xrightarrow[\quad\quad]{c}GFP\ \ \ \ \ \ \ \ \ \ \ (3.5)$$</p> | ||
+ | <p>The formula (3.5) indicates the production of the effectors. The downstream gene can be expressed when it is reversed to open state. Here, GFP is taken as an example, and the expression is constant after the rectification reaction.</p> | ||
<p>The figure below shows the opening of the downstream gene. 0 means closed; 1 means open.</p> | <p>The figure below shows the opening of the downstream gene. 0 means closed; 1 means open.</p> | ||
− | + | <p>(a):</p> | |
− | + | <figure> | |
− | + | <img src="https://static.igem.org/mediawiki/2018/3/39/T--NAU-China--model9.svg"> | |
− | + | </figure> | |
<p>(b):</p> | <p>(b):</p> | ||
− | + | <figure> | |
− | + | <img src="https://static.igem.org/mediawiki/2018/0/08/T--NAU-China--model10.svg"> | |
− | + | <figcaption> | |
− | Each subplot in (b) is corresponding to the corresponding subplot in (a). | + | <b>Fig. 3.2. </b> (a) Relationship between TEV number level and the state of the downstream gene in stable state. (b)Expression of GFP in different TEV number levels,<br> |
− | </figcaption></figure> | + | Each subplot in (b) is corresponding to the corresponding subplot in (a). |
− | <p>The results meet our expectation: The system doesn’t express GFP in low TEV | + | </figcaption> |
− | + | </figure> | |
+ | <p>The results meet our expectation: The system doesn’t express GFP in low TEV number level. It begins to express GFP when the concentration level reaches the threshold.</p> | ||
+ | |||
</div> | </div> | ||
</div> | </div> | ||
− | + | <div class="main-content"> | |
<div class="textblock"> | <div class="textblock"> | ||
− | <h1> | + | <h1 id="Part4 Significance of the cell culture environment">Part4 Significance of the cell culture environment</h1> |
− | + | <p>If the normal number of the signal molecules can be added into the culture medium of cells, it will improve the system stability by 10%.</p> | |
− | + | <div class="section"> | |
<h2>First stimulation effect</h2> | <h2>First stimulation effect</h2> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/5/56/T--NAU-China--model11.svg"> |
<figcaption><b>Fig. 4.1. </b> Steady change of cell contact signal, and the amount of synNotch</figcaption> | <figcaption><b>Fig. 4.1. </b> Steady change of cell contact signal, and the amount of synNotch</figcaption> | ||
</figure> | </figure> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/5/5e/T--NAU-China--model12.svg"> |
<figcaption><b>Fig. 4.2. </b>Changes in the amount of intracellular TEV in stable cells exposed to different gradient signals</figcaption> | <figcaption><b>Fig. 4.2. </b>Changes in the amount of intracellular TEV in stable cells exposed to different gradient signals</figcaption> | ||
</figure> | </figure> | ||
− | + | ||
</div> | </div> | ||
− | + | <div class="section"> | |
<h2>Accompanying stimulation effect</h2> | <h2>Accompanying stimulation effect</h2> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/6/ | + | <img src="https://static.igem.org/mediawiki/2018/6/60/T--NAU-China--model13.svg"> |
<figcaption><b>Fig. 4.3. </b> Changes in the amount of synNotch in steady-state cell contact with different gradient signals</figcaption> | <figcaption><b>Fig. 4.3. </b> Changes in the amount of synNotch in steady-state cell contact with different gradient signals</figcaption> | ||
</figure> | </figure> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/4/42/T--NAU-China--model14.svg"> |
<figcaption><b>Fig. 4.4. </b> Changes in the amount of TEV in steady-state cell contact with different gradient signals</figcaption> | <figcaption><b>Fig. 4.4. </b> Changes in the amount of TEV in steady-state cell contact with different gradient signals</figcaption> | ||
</figure> | </figure> | ||
− | <p>It can be seen from the calculation results that the cells need to be stimulated by adding signal molecules of normal | + | <p>It can be seen from the calculation results that the cells need to be stimulated by adding signal molecules of normal number in the culture medium, which can improve the stability of the system by 10%. In the non-irritated culture, the stable amount of synNotch is the maximum amount that cell membrane can carry. In the stimulus culture environment, the stable value of the synNotch is the residual amount of the synNotch at the corresponding signal level. Therefore, in the required environment, the amount of change will not be changed a lot, and the released TEV will not have shock.</p> |
− | + | ||
</div> | </div> | ||
</div> | </div> | ||
Line 263: | Line 256: | ||
<div class="main-content"> | <div class="main-content"> | ||
<div class="textblock"> | <div class="textblock"> | ||
− | <h1> | + | <h1 id="Part5 Recombinase filtering and de-shocking">Part5 Recombinase filtering and de-shocking</h1> |
<div class="section"> | <div class="section"> | ||
<h2>Rec-inhibitor noise reduction</h2> | <h2>Rec-inhibitor noise reduction</h2> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/c/c0/T--NAU-China--model15.svg"> |
<figcaption><b>Fig. 5.1. </b> Amount of rec expression in the case of no rec-inhibitor noise reduction</figcaption> | <figcaption><b>Fig. 5.1. </b> Amount of rec expression in the case of no rec-inhibitor noise reduction</figcaption> | ||
</figure> | </figure> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/c/c8/T--NAU-China--model16.svg"> |
<figcaption><b>Fig. 5.2. </b> Changes in the amount of TEV in steady-state cell contact with different gradient signals</figcaption> | <figcaption><b>Fig. 5.2. </b> Changes in the amount of TEV in steady-state cell contact with different gradient signals</figcaption> | ||
</figure> | </figure> | ||
<p>By comparing the noise reduction, the expression discrimination of the switch signal (rec) can be significantly improved.</p> | <p>By comparing the noise reduction, the expression discrimination of the switch signal (rec) can be significantly improved.</p> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/e/ed/T--NAU-China--model17.svg"> |
<figcaption><b>Fig. 5.3. </b> Switching of downstream genes without rec-inhibitor</figcaption> | <figcaption><b>Fig. 5.3. </b> Switching of downstream genes without rec-inhibitor</figcaption> | ||
</figure> | </figure> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/0/08/T--NAU-China--model18.svg"> |
<figcaption><b>Fig. 5.4. </b> Switching of downstream genes under rec-inhibitor noise reduction</figcaption> | <figcaption><b>Fig. 5.4. </b> Switching of downstream genes under rec-inhibitor noise reduction</figcaption> | ||
</figure> | </figure> | ||
− | + | <p>In general, rec-inhibitor noise reduction doesn’t have significant effect on the total on-off threshold, so our experiment group doesn’t add this reaction. However, for the accuracy of the system we believe that it is worthwhile to reduce noise when time and cost allow.</p> | |
− | + | ||
</div> | </div> | ||
− | + | <div class="section"> | |
<h2>RDF-inhibitor shock remover</h2> | <h2>RDF-inhibitor shock remover</h2> | ||
− | <p>If we do not add RDF-inhibitor, and rely solely on rec and RDF to competitively control downstream genes, we will find that our system may be in a state of constant shock. | + | <p> |
− | The following are the simulated results: | + | If we do not add RDF-inhibitor, and rely solely on rec and RDF to competitively control downstream genes, we will find that our system may be in a state of constant shock. |
− | </p> | + | The following are the simulated results: |
− | + | </p> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <figure> |
+ | <img src="https://static.igem.org/mediawiki/2018/e/ef/T--NAU-China--model19.svg"> | ||
<figcaption><b>Fig. 5.5. </b> Switching of downstream genes without RDF-inhibitor noise reduction. X axis represents whether the downstream gene is active.</figcaption> | <figcaption><b>Fig. 5.5. </b> Switching of downstream genes without RDF-inhibitor noise reduction. X axis represents whether the downstream gene is active.</figcaption> | ||
</figure> | </figure> | ||
− | + | <figure> | |
− | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/6/67/T--NAU-China--model20.svg"> |
<figcaption><b>Fig. 5.6. </b> No RDF-inhibitor noise reduction downstream gene: GFP release</figcaption> | <figcaption><b>Fig. 5.6. </b> No RDF-inhibitor noise reduction downstream gene: GFP release</figcaption> | ||
</figure> | </figure> | ||
− | <p>We can see that | + | <p>We can see that the switching of downstream gene are controlled by signal number. Only when number level is higher than 5, the switch will be open. However, because the expression of RDF is not restricted, it will reset the downstream genes directly, which will turn off the switch. Then the gene will be turned on by the effect of rec again, which will cause periodical shock in the system.</p> |
− | + | <p>When we simulate the reaction, we assume that the open gene is not the same matter as the closed gene. Those two matters converse to each other by reaction (3.1), and (3.2). Because of the shock of effectors’ gene, the expression of GFP is discontinuous.</p> | |
− | + | <p>When the system needs to be turned on, oscillations may occur, which may damage the DNA.</p> | |
− | + | ||
</div> | </div> | ||
</div> | </div> | ||
Line 310: | Line 304: | ||
<div class="main-content"> | <div class="main-content"> | ||
<div class="textblock"> | <div class="textblock"> | ||
− | <h1 | + | <h1 id="Part6: Robustness of the model">Part6 Robustness of the model</h1> |
<p> | <p> | ||
− | + | After adding a lot of disturbances to the previously determined parameters, the calculation results are as follows: | |
</p> | </p> | ||
<figure> | <figure> | ||
− | + | <img src="https://static.igem.org/mediawiki/2018/0/03/T--NAU-China--model21.svg"> | |
− | + | <figcaption><b>Fig. 6.1. </b> The switch diagram after adding a lot of disturbances</figcaption> | |
− | + | </figure> | |
<p>The stability of the algorithm and our path is verified by adding a lot of perturbations on the parameters.</p> | <p>The stability of the algorithm and our path is verified by adding a lot of perturbations on the parameters.</p> | ||
+ | <p><a href="https://static.igem.org/mediawiki/2018/c/cd/T--NAU-China--modelinterpretation.pdf" target="_blank">Full Text Access</a></p> | ||
</div> | </div> | ||
</div> | </div> | ||
+ | |||
+ | |||
+ | |||
+ | </div> | ||
+ | |||
+ | |||
+ | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
</body> | </body> | ||
</html> | </html> | ||
+ | {{NAU-CHINA/footer}} |
Latest revision as of 02:48, 18 October 2018
Model
Model Illustration
Part1 Extracellular signal converting into intracellular signal
Part2 Internal signal regulating recombinase production
Part3 Recombinase turns on downstream effect expression
Part4 Significance of the cell culture environment
Part5 Recombinase filtering and de-shocking
Part1 Extracellular signal converting into intracellular signal
Tips: Skip part1 if the gene circuit is directly regulated by intracellular signal.
SynNotch can convert an extracellular signal into an intracellular signal . Here, we take synNotch-TEV as an example. We use the number of TEV as the intracellular signal after conversion. TEV is an enzyme which is able to release the repressed promoter.
In this section, 4 reactions are considered:
$$G_{synNotch-TEV}\xrightarrow[\quad\quad]{C11}{synNotch-TEV}\ \ \ \ \ \ \ \ \ \ \ (1.1)$$
$$synNotch-TEV+Ts\xrightarrow[\quad\quad]{C12}U+TEV\ \ \ \ \ \ \ \ \ \ \ (1.2)$$
$$synNotch-TEV\xrightarrow[\quad\quad]{C13}\emptyset\ \ \ \ \ \ \ \ \ \ \ (1.3)$$
$$TEV\xrightarrow[\quad\quad]{C14}\emptyset\ \ \ \ \ \ \ \ \ \ \ (1.4)$$
The specific contents are as follows:
$$G_{synNotch-TEV}\xrightarrow[\quad\quad]{C11}{synNotch-TEV}\ \ \ \ \ \ \ \ \ \ \ (1.1)$$
(1.1) represents the expression of synNotch which contains the TEV. Although the actual expression progress of the protein involves many chemical reactions from transcription to modification, these intermediate reactions satisfy the following conditions:
a)The factors affecting these reactions in the system are stable, such as the amount of amino acids, nucleotides, and helicases;
b)The gene is regulated by a constitutive promoter and its expression is not regulated by other products of the system;
c)The product of an intermediate reaction has no effects on the system functions.
We call these three conditions the post-integration conditions [1], and call these reactions the post-integration reaction. In the same way, the expression and degradation of other proteins involved in our model all suit the conditions of the refinement, so they can be considered as one single chemical reaction.
Obviously (1.1) is a post-integration reaction, and the reaction rate is a constant, for the controlling conditions are unchanged.
$$synNotch-TEV+Ts\xrightarrow[\quad\quad]{C12}U+TEV\ \ \ \ \ \ \ \ \ \ \ (1.2)$$
(1.2) represents that the TEV protein is released by the intracellular domain of the sporadic synNotch when the synNotch-TEV protein was stimulated by an external Target signal. The rate of the reaction is calculated by the law of mass action which is in proportion to the number of receptors on the cell membrane and the number of external target signals.
$$synNotch-TEV\xrightarrow[\quad\quad]{C13}\emptyset\ \ \ \ \ \ \ \ \ \ \ (1.3)$$
(1.3) represents the degradation of synNotch membrane proteins linked to TEV proteins The degradation rate is directly proportional to its protein content.
$$TEV\xrightarrow[\quad\quad]{C14}\emptyset\ \ \ \ \ \ \ \ \ \ \ (1.4)$$
(1.4) represents the degradation of TEV protein.
The degradation rate is directly proportional to its protein content.
The figures below show the result of simulation using the Gillespie algorithm.
It is concluded that the initial value of synNotch-TEV is stable at about 20,000 before the Target signal is received, which is used as the initial condition for subsequent translation.
Conclusion: When a stabilized cell contacts different levels of extracellular signals, synNotch membrane protein will release different amounts of TEV. The amount of TEV protein is in proportion of the strength of target signals.
Part2 Internal signal regulating recombinase production
In this part, intracellular proteins can degrade or destroy the repressors, reducing their number. Different levels of protein have different effects on the relieving repression. Then, using Hill equation, we can calculate the expression rate of the recombinase based on the repressor protein number. At the same time, we use the error-reducing reaction to solve the problem of the extra recombinase. Here, we choose TEV as intracellular protein, and TetR as repressor protein.
In this section, 12 reactions are considered:
$$G_{TetR}\xrightarrow[\quad\quad]{C21}TetR\ \ \ \ \ \ \ \ \ \ \ (2.1)$$
$$TetR\xrightarrow[\quad\quad]{C22}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.2)$$
$$2TetR+tetO\xrightarrow[\quad\quad]{C23}TetOR\ \ \ \ \ \ \ \ \ \ \ (2.3)$$
$$TEV+TetOR\xrightarrow[\quad\quad]{C24}tetO+TEV+U\ \ \ \ \ \ \ \ \ \ \ (2.4)$$
$$TetOR\xrightarrow[\quad\quad]{C25}tetO+U\ \ \ \ \ \ \ \ \ \ \ (2.5)$$
$$G_{rec}\xrightarrow[\quad\quad\quad\quad\quad\quad]{C31-TetOR}rec\ \ \ \ \ \ \ \ \ \ \ (2.6)$$
$$rec\xrightarrow[\quad\quad]{C33}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.7)$$
$$G_{RDF-inhibitor}\xrightarrow[\quad\quad\quad\quad\quad\quad]{C32-TetOR}RDF-inhibitor\ \ \ \ \ \ \ \ \ \ \ (2.8)$$
$$RDF-inhibitor\xrightarrow[\quad\quad]{C34}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.9)$$
$$G_{rec-inhibitor}\xrightarrow[\quad\quad]{C35}rec-inhibitor\ \ \ \ \ \ \ \ \ \ \ (2.10)$$
$$rec-inhibitor\xrightarrow[\quad\quad]{C36}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.11)$$
$$rec+rec-inhibitor\xrightarrow[\quad\quad]{C37}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.12)$$
The specific content is as follows. (As for each chart in this part, we choose the quantity of substance as the vertical axis, and 0.01 hour as the unit of the horizontal axis. Different subplots correspond to different levels of intracellular protein number).
$$G_{TetR}\xrightarrow[\quad\quad]{C21}TetR\ \ \ \ \ \ \ \ \ \ \ (2.1)$$
(2.1) indicates the expression of the TetR protein. Through the analysis of (1.1) in part1, it is easy to conclude that (2.1) is a post-integration reaction, and the rate of the reaction is constant.
$$TetR\xrightarrow[\quad\quad]{C22}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.2)$$
(2.2) indicates the degradation of TetR protein. The degradation rate is directly proportional to its protein content.
$$2TetR+tetO\xrightarrow[\quad\quad]{C23}TetOR\ \ \ \ \ \ \ \ \ \ \ (2.3)$$
(2.3) represents the generation of the binary complex TetOR, and the reaction rate is calculated according to the law of mass action in the model.
$$TEV+TetOR\xrightarrow[\quad\quad]{C24}tetO+TEV+U\ \ \ \ \ \ \ \ \ \ \ (2.4)$$
(2.4) indicates that the TEV is a separate TetOR binary complex. This reaction rate is calculated according to the law of mass action.
$$TetOR\xrightarrow[\quad\quad]{C25}tetO+U\ \ \ \ \ \ \ \ \ \ \ (2.5)$$
(2.5) indicates that a small amount of TetOR is self-separate. The reaction rate is calculated according to the law of mass action.
The figure below shows the calculated results.
Using the above calculated amount of repressor, the Hill equation is used to determine the expression rate of the recombinase.
$$G_{rec}\xrightarrow[\quad\quad\quad\quad\quad\quad]{C31-TetOR}rec\ \ \ \ \ \ \ \ \ \ \ (2.6)$$
(2.6) indicates the generation of rec (rec is the abbreviation of recombinase) , which is also a post-integration reaction, but the expression rate is regulated by the repressor TetR. The reaction rate in the model is calculated using the Hill equation, and the rate is ((N/k) ^n+1) ^-1. Accordingly, N is the number of TetR, and K is the half-maximum of the suppression effect.
$$rec\xrightarrow[\quad\quad]{C33}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.7)$$
(2.7) indicates degradation of rec. The degradation rate is directly proportional to rec number.
$$G_{RDF-inhibitor}\xrightarrow[\quad\quad\quad\quad\quad\quad]{C32-TetOR}RDF-inhibitor\ \ \ \ \ \ \ \ \ \ \ (2.8)$$
(2.8) indicates the generation of RDF. Like rec, it is also a post-integration reaction, whose rate is regulated by the repressor TetR. The reaction rate in the model is calculated by the Hill equation, and the rate is ((N/k) ^n+1) ^-1. Similarly, N is the number of TetR, and K is the half-maximum of suppression effect.
$$RDF-inhibitor\xrightarrow[\quad\quad]{C34}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.9)$$
Degradation of RDF-inhibitor is represented by (2.9). The degradation rate is directly proportional to the number of RDF-inhibitor. The roles and meanings of this reaction will be discussed in the next part.
$$G_{rec-inhibitor}\xrightarrow[\quad\quad]{C35}rec-inhibitor\ \ \ \ \ \ \ \ \ \ \ (2.10)$$
(2.10) indicates the generation of the rec-inhibitor, which is also a post-integration reaction, and the expression rate is constant.
$$rec-inhibitor\xrightarrow[\quad\quad]{C36}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.11)$$
(2.11) indicates the degradation of the rec-inhibitor. The degradation rate is directly proportional to rec-inhibitor number.
$$rec+rec-inhibitor\xrightarrow[\quad\quad]{C37}\emptyset\ \ \ \ \ \ \ \ \ \ \ (2.12)$$
(2.12) indicates that the rec-inhibitor combined with the leaked rec. It is used to filter rec. The reaction rate was calculated according to the law of mass action and enables the low expression of rec to be combined to ensure the stability of the system.
It is concluded that internal signal TEV regulates recombinase production. According to (2.4), different concentrations of TEV correspond to different concentrations of TetR. By (2.6), the concentrations of TetR affect rec concentrations. Through the simulation of (2.12), a small amount of rec-inhibitor can almost completely eliminate the rec leakage from the TetR system. This pathway can significantly improves the system stability.
Part3 Recombinase turns on downstream effect expression
In this part, we use the above calculated amount of the rec expression at different signal levels. Using rec and RDF, we can switch the state and calculate the expression of the downstream gene. In this section, 5 reactions are considered:
$$rec+\overline{G}_{downstream}\xrightarrow[\quad\quad]{C41}rec+G_{downstream}\ \ \ \ \ \ \ \ \ \ \ (3.1)$$
$$RDF+G_{dowmstream}\xrightarrow[\quad\quad]{C61}RDF+\overline{G}_{downstream}\ \ \ \ \ \ \ \ \ \ \ (3.2)$$
$$G_{RDF}\xrightarrow[\quad\quad]{C}RDF\ \ \ \ \ \ \ \ \ \ \ (3.3)$$
$$RDF-inhibitor+RDF\xrightarrow[\quad\quad]{C34}\emptyset\ \ \ \ \ \ \ \ \ \ \ \ (3.4)$$
$$G_{gfp}\xrightarrow[\quad\quad]{C}GFP\ \ \ \ \ \ \ \ \ \ \ (3.5)$$
$$rec+\overline{G}_{downstream}\xrightarrow[\quad\quad]{C41}rec+G_{downstream}\ \ \ \ \ \ \ \ \ \ \ (3.1)$$
(3.1) represents that recombinase specifically recognizes the DNA sequence between the two sites and reverses it, so that the downstream gene expression is reversed to an open state. (Downstream genes are GFP and RDF Gene)
$$RDF+G_{dowmstream}\xrightarrow[\quad\quad]{C61}RDF+\overline{G}_{downstream}\ \ \ \ \ \ \ \ \ \ \ (3.2)$$
(3.2) indicates that the reverse recombination factor specifically recognizes the two sites after the recombinase worked on the DNA, and makes the reversal DNA sequence between the sites restored to the shape.
$$G_{RDF}\xrightarrow[\quad\quad]{c}RDF\ \ \ \ \ \ \ \ \ \ \ (3.3)$$
(3.3) indicates the expression of the reverse recombination factor, which exists only when the downstream gene is turned into an open state. The expression of the RDF process is not regulated, and also is a post-integration reaction, so the expression rate is constant.
$$RDF-inhibitor+RDF\xrightarrow[\quad\quad]{C34}\emptyset\ \ \ \ \ \ \ \ \ \ \ \ (3.4)$$
(3.4) indicates that the RDF-inhibitor is combined with the intracellular RDF when the system requests downstream gene opening. The reaction rate is calculated according to the law of mass action. It makes our system stable in the open state. When it is not requested to open, the increase of RDF is not limited, so that it can function to shut down the system and realize real-time feedback.
In this simulation, the off state is the initial state. When the TEV signal number level is less than 5, there is not enough rec to turn on the downstream gene,and RDF is not expressed; when the TEV signal number level is higher than 5, rec and RDF-inhibitor begin to express simultaneously after de-repression. Then,RDF and GFP are simultaneously expressed after the rec turns on the downstream gene. The combination of RDF-inhibitor and RDF ensures that there is not enough RDF in the cell to shut down the downstream gene, which makes the system stable.
$$G_{gfp}\xrightarrow[\quad\quad]{c}GFP\ \ \ \ \ \ \ \ \ \ \ (3.5)$$
The formula (3.5) indicates the production of the effectors. The downstream gene can be expressed when it is reversed to open state. Here, GFP is taken as an example, and the expression is constant after the rectification reaction.
The figure below shows the opening of the downstream gene. 0 means closed; 1 means open.
(a):
(b):
The results meet our expectation: The system doesn’t express GFP in low TEV number level. It begins to express GFP when the concentration level reaches the threshold.
Part4 Significance of the cell culture environment
If the normal number of the signal molecules can be added into the culture medium of cells, it will improve the system stability by 10%.
First stimulation effect
Accompanying stimulation effect
It can be seen from the calculation results that the cells need to be stimulated by adding signal molecules of normal number in the culture medium, which can improve the stability of the system by 10%. In the non-irritated culture, the stable amount of synNotch is the maximum amount that cell membrane can carry. In the stimulus culture environment, the stable value of the synNotch is the residual amount of the synNotch at the corresponding signal level. Therefore, in the required environment, the amount of change will not be changed a lot, and the released TEV will not have shock.
Part5 Recombinase filtering and de-shocking
Rec-inhibitor noise reduction
By comparing the noise reduction, the expression discrimination of the switch signal (rec) can be significantly improved.
In general, rec-inhibitor noise reduction doesn’t have significant effect on the total on-off threshold, so our experiment group doesn’t add this reaction. However, for the accuracy of the system we believe that it is worthwhile to reduce noise when time and cost allow.
RDF-inhibitor shock remover
If we do not add RDF-inhibitor, and rely solely on rec and RDF to competitively control downstream genes, we will find that our system may be in a state of constant shock. The following are the simulated results:
We can see that the switching of downstream gene are controlled by signal number. Only when number level is higher than 5, the switch will be open. However, because the expression of RDF is not restricted, it will reset the downstream genes directly, which will turn off the switch. Then the gene will be turned on by the effect of rec again, which will cause periodical shock in the system.
When we simulate the reaction, we assume that the open gene is not the same matter as the closed gene. Those two matters converse to each other by reaction (3.1), and (3.2). Because of the shock of effectors’ gene, the expression of GFP is discontinuous.
When the system needs to be turned on, oscillations may occur, which may damage the DNA.
Part6 Robustness of the model
After adding a lot of disturbances to the previously determined parameters, the calculation results are as follows:
The stability of the algorithm and our path is verified by adding a lot of perturbations on the parameters.