Difference between revisions of "Team:USP-Brazil/Model"

 
(6 intermediate revisions by 3 users not shown)
Line 21: Line 21:
 
           <article>
 
           <article>
 
           <div class="mod-content modelo">
 
           <div class="mod-content modelo">
            <h2><a name="model">Model</a></h2>
+
            <h2><a name="model">Model</a></h2>
            <p>As discussed in the Project, we can consider some basic variables of interest in our project:  
+
            <p>As discussed in the <a href="https://2018.igem.org/Team:USP-Brazil/Design">Project,</a> we can consider some basic points of interest in our project:  
Synthetase (S), HSL (H), Receptor (R) and Protein of interest (R),  
+
Synthase (S), HSL (H), Receptor (R) and Protein of interest (R),  
which can be considered as the variables of the systems. Also we can consider the complex formed between the HSL and the Receptor (C).  
+
which can be used as the variables of the systems. Also we can consider the complex formed between the HSL and the Receptor (C). </p>
<strong>NOTE:</strong> For a simplified visualization, we adopt, during the whole discussion of the modelling,  
+
<p> <strong>NOTE:</strong> For a simplified visualization, we adopt, during the whole discussion of the modelling,  
the letters indicated between parentheses as representation of the variables.
+
the letters indicated between parentheses as representation of the variables.
            </p>
+
            </p>
<p>Before starting the model, itself, we need to stablish some assumptions:</p>
+
<p>Before starting the model itself, we need to stablish some assumptions:</p>
<ul>
+
<ul>
<li>The population, during the measurements, is at dynamic equilibria, therefore we can consider it as a constant;</li>
+
<li>The population, during the measurements, is at dynamic equilibrium, therefore we can consider it as a constant;</li>
<li>Transcription, when constitutive, and translation events are very low effected by stochastic effects;</li>
+
<li>Transcription, when constitutive, and translation events are very lowly affected by stochastic effects;</li>
<li>Also, considering the transcription and translation can be coupled in a single term of an equation, due to a  
+
<li>Also, the transcription and translation can be coupled in a single term of an equation, due to a  
linear relation between mRNA and protein (I.e. the translation is not limited by ribosomes in the cell)</li>
+
linear relation between mRNA and protein (therefore, the translation is not limited by ribosomes in the cell)</li>
 
</ul>
 
</ul>
 
<p>We can, then, propose the following reactions:</p>
 
<p>We can, then, propose the following reactions:</p>
Line 43: Line 43:
 
&\rightarrow R \\
 
&\rightarrow R \\
 
R &\rightarrow\\
 
R &\rightarrow\\
H + R &\leftrightarrow C\\
+
H + R &\leftrightharpoons C\\
 
C &\rightarrow\\
 
C &\rightarrow\\
 
C & \rightarrow P\\
 
C & \rightarrow P\\
 
P &\rightarrow \\
 
P &\rightarrow \\
 
\end{align}$$
 
\end{align}$$
 +
<p>Based on this, we can initially construct a graphic view of the model, applying mass action kinects law using now the terms of the equations:</p>
 +
<p>This kind of figure makes it easier to visualize the model and the relations between the classes of the model.
 +
Then, the way to read is the following:</p>
 +
<ul>
 +
<li>Sphere – Variable of the system</li>
 +
<li>Line – Relation between the systems</li>
 +
<li>Arrow - indicates the directionality of the system (back of the arrow decreases and front of the arrow grows</li>
 +
<li>Bar (In the line) – Indicates no decreasing of the variable with the line</li>
 +
</ul>
 +
 +
                            <center><img src="https://static.igem.org/mediawiki/2018/d/da/T--USP-Brazil--model_scheme_1.png" style="width: 75%"></center>
 +
 +
<p>The term above the line/on the right of the line represents the parameter that conducts the relation.</p>
 +
<p>The index in the parameter represent in which variable is the parameter assigned. For a summary of the variables and parameters, see the table at the end of the page.</p>
 +
<p>With this data we can finally construct the Mathematical Model that will describe our systems!
 +
<br>
 +
However, due to our design of the project, the protein of interest is going to be used as a relation between YFP and CFP.
 +
This design allows us to do an easier comparison between the systems, having independence of OD measurement and reducing the deviations.
 +
We then stablish the variable ρ as YFP/CFP.</p>
 +
<p>Finally, looking at those representations, we can construct the following differential equations:</p>
 +
$$ \frac{dS}{dt} = \alpha_{S} - \mu_{S} S \\
 +
\frac{dH}{dt} = \alpha_{H} S + \theta C - H( \mu_{H} H + \beta R)\\
 +
\frac{dR}{dt} = \alpha_{R} + \theta C - R( \mu_{R} R + \beta H) \\
 +
\frac{dC}{dt} = \beta R H - C (\mu_{C} + \theta)\\$$
 +
 +
<p>For ρ, there must be a more specifically mathematical work:<br>
 +
So as we define ρ as Yf/Cf, as Yf being YFP e Cf as CFP, we define the following derivation:</p>
 +
 +
$$
 +
\frac{d \rho}{dt} = \frac{Yf'Cf - Cf'Y}{Cf^2} \parallel Cf'  = \frac{dCf}{dt} = 0 (By\; definition)\\
 +
\Rightarrow \frac{d \rho}{dt} = \frac{Yf'Cf - 0'Y}{Cf^2} = \frac{Yf'Cf - 0'Y}{Cf^2} = \frac{Yf'*Cf}{Cf^2}\\
 +
\Rightarrow \frac{d \rho}{dt} = \frac{Yf'}{Cf} = \frac{\alpha_{Yf}C -\mu_{Yf}Y}{Cf} \\
 +
\Rightarrow \frac{d \rho}{dt} = \frac{\alpha_{Yf}C}{Cf} -\mu_{Yf}\rho
 +
$$
 +
 +
<p>Using those general equations, we were able to estimate αYf from each system, establishing the crosstalk!</p>
 +
 +
<table class="olha-a-tabela">
 +
<tr title="I really like this table!"><!-- O comando tr indica a linha -->
 +
  <th>Parameter/Variable Name</th><!-- Cada th é uma coluna do header da tabela -->
 +
  <th>Description</th>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>S</td><!-- cada td é um dos seus itens -->
 +
  <td>System Synthase</td>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>H</td><!-- cada td é um dos seus itens -->
 +
  <td>System HSL (Produced bt the synthase)</td>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>R</td><!-- cada td é um dos seus itens -->
 +
  <td>System Receptor</td>
 +
</tr>
 +
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>C</td><!-- cada td é um dos seus itens -->
 +
  <td>Complex formed by the HSL and Receptor</td>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>Cf</td><!-- cada td é um dos seus itens -->
 +
  <td>CFP used as relative measure</td>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>$$\rho$$</td><!-- cada td é um dos seus itens -->
 +
  <td>Relative Fluorescence (YFP/CFP)</td>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>$$\alpha_{N}$$</td><!-- cada td é um dos seus itens -->
 +
  <td>Production/Expression of the variable, indicated in the index</td>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>$$\mu_{N}$$</td><!-- cada td é um dos seus itens -->
 +
  <td>Degradation of the variable, indicated in the index</td>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>$$\beta$$</td><!-- cada td é um dos seus itens -->
 +
  <td>Formation of the HSL/Receptor Complex</td>
 +
</tr>
 +
 +
<tr><!-- O comando tr indica a linha -->
 +
  <td>$$\theta$$</td><!-- cada td é um dos seus itens -->
 +
  <td>Separation of the Complex HSL/Receptor</td>
 +
</tr>
 +
 +
</table>
 
           </div>
 
           </div>
 
         </article>
 
         </article>

Latest revision as of 03:54, 18 October 2018

Wiki - iGEM Brazil

Model

As discussed in the Project, we can consider some basic points of interest in our project: Synthase (S), HSL (H), Receptor (R) and Protein of interest (R), which can be used as the variables of the systems. Also we can consider the complex formed between the HSL and the Receptor (C).

NOTE: For a simplified visualization, we adopt, during the whole discussion of the modelling, the letters indicated between parentheses as representation of the variables.

Before starting the model itself, we need to stablish some assumptions:

  • The population, during the measurements, is at dynamic equilibrium, therefore we can consider it as a constant;
  • Transcription, when constitutive, and translation events are very lowly affected by stochastic effects;
  • Also, the transcription and translation can be coupled in a single term of an equation, due to a linear relation between mRNA and protein (therefore, the translation is not limited by ribosomes in the cell)

We can, then, propose the following reactions:

$$ \begin{align} &\rightarrow S \\ S &\rightarrow \\ S & \rightarrow H\\ H &\rightarrow \\ &\rightarrow R \\ R &\rightarrow\\ H + R &\leftrightharpoons C\\ C &\rightarrow\\ C & \rightarrow P\\ P &\rightarrow \\ \end{align}$$

Based on this, we can initially construct a graphic view of the model, applying mass action kinects law using now the terms of the equations:

This kind of figure makes it easier to visualize the model and the relations between the classes of the model. Then, the way to read is the following:

  • Sphere – Variable of the system
  • Line – Relation between the systems
  • Arrow - indicates the directionality of the system (back of the arrow decreases and front of the arrow grows
  • Bar (In the line) – Indicates no decreasing of the variable with the line

The term above the line/on the right of the line represents the parameter that conducts the relation.

The index in the parameter represent in which variable is the parameter assigned. For a summary of the variables and parameters, see the table at the end of the page.

With this data we can finally construct the Mathematical Model that will describe our systems!
However, due to our design of the project, the protein of interest is going to be used as a relation between YFP and CFP. This design allows us to do an easier comparison between the systems, having independence of OD measurement and reducing the deviations. We then stablish the variable ρ as YFP/CFP.

Finally, looking at those representations, we can construct the following differential equations:

$$ \frac{dS}{dt} = \alpha_{S} - \mu_{S} S \\ \frac{dH}{dt} = \alpha_{H} S + \theta C - H( \mu_{H} H + \beta R)\\ \frac{dR}{dt} = \alpha_{R} + \theta C - R( \mu_{R} R + \beta H) \\ \frac{dC}{dt} = \beta R H - C (\mu_{C} + \theta)\\$$

For ρ, there must be a more specifically mathematical work:
So as we define ρ as Yf/Cf, as Yf being YFP e Cf as CFP, we define the following derivation:

$$ \frac{d \rho}{dt} = \frac{Yf'Cf - Cf'Y}{Cf^2} \parallel Cf' = \frac{dCf}{dt} = 0 (By\; definition)\\ \Rightarrow \frac{d \rho}{dt} = \frac{Yf'Cf - 0'Y}{Cf^2} = \frac{Yf'Cf - 0'Y}{Cf^2} = \frac{Yf'*Cf}{Cf^2}\\ \Rightarrow \frac{d \rho}{dt} = \frac{Yf'}{Cf} = \frac{\alpha_{Yf}C -\mu_{Yf}Y}{Cf} \\ \Rightarrow \frac{d \rho}{dt} = \frac{\alpha_{Yf}C}{Cf} -\mu_{Yf}\rho $$

Using those general equations, we were able to estimate αYf from each system, establishing the crosstalk!

Parameter/Variable Name Description
S System Synthase
H System HSL (Produced bt the synthase)
R System Receptor
C Complex formed by the HSL and Receptor
Cf CFP used as relative measure
$$\rho$$ Relative Fluorescence (YFP/CFP)
$$\alpha_{N}$$ Production/Expression of the variable, indicated in the index
$$\mu_{N}$$ Degradation of the variable, indicated in the index
$$\beta$$ Formation of the HSL/Receptor Complex
$$\theta$$ Separation of the Complex HSL/Receptor