"
Team:BCCS-Bristol/Modelling/GRN/Derivation
From 2010.igem.org
m |
m |
||
| Line 5: | Line 5: | ||
<!--Discussion of how we derived GRN model....--> | <!--Discussion of how we derived GRN model....--> | ||
| - | ==GRN | + | ==GRN Model== |
A Gene Regulatory Network (GRN) model is an established way of representing the biochemical reactions within a bacterial cell [[(7)]]. Representing a network of interacting components using differential equations can be very useful, mainly because one can make analytical and qualitative assessments of a system’s behaviour with very limited information on parameters. This GRN represents a protein called NsrR that represses the production of mRNA encoding GFP, and so indirectly repressing the production of GFP. NsrR can bind to nitrates, effectively destroying its ability to repress the system. | A Gene Regulatory Network (GRN) model is an established way of representing the biochemical reactions within a bacterial cell [[(7)]]. Representing a network of interacting components using differential equations can be very useful, mainly because one can make analytical and qualitative assessments of a system’s behaviour with very limited information on parameters. This GRN represents a protein called NsrR that represses the production of mRNA encoding GFP, and so indirectly repressing the production of GFP. NsrR can bind to nitrates, effectively destroying its ability to repress the system. | ||
| Line 11: | Line 11: | ||
[[Image:GRN.jpg|frameless|center|upright=4|agrEcoli Gene Regulatory Network Equations]] | [[Image:GRN.jpg|frameless|center|upright=4|agrEcoli Gene Regulatory Network Equations]] | ||
| - | + | ==Analysis== | |
The most interesting relationship in this system of equations is the relationship between the concentration of NsrR repressor protein and the concentration of GFP in the system. A useful way of visualising and analysing the system’s dynamics is using a phase plane portrait in 2 dimensions. However, this means eliminating one of the variables. This can be done by making an equilibrium assumption about the amount of mRNA encoding GFP (mGFP for short). To do this, we set Equation 2 equal to 0. Rearranging in terms of [mGFP], this gives: | The most interesting relationship in this system of equations is the relationship between the concentration of NsrR repressor protein and the concentration of GFP in the system. A useful way of visualising and analysing the system’s dynamics is using a phase plane portrait in 2 dimensions. However, this means eliminating one of the variables. This can be done by making an equilibrium assumption about the amount of mRNA encoding GFP (mGFP for short). To do this, we set Equation 2 equal to 0. Rearranging in terms of [mGFP], this gives: | ||
| Line 28: | Line 28: | ||
[[Image:BCCS_GRN_EQ_9.jpg|frameless|center|upright=4|BCCS-Bristol GRN equation 9]] | [[Image:BCCS_GRN_EQ_9.jpg|frameless|center|upright=4|BCCS-Bristol GRN equation 9]] | ||
| + | |||
| + | ==Fixed Point Solution== | ||
Since all the constants in this system of equations are constrained to be real semi-positive numbers (i.e. ≥0), both eigenvalues are always real and negative. This means that the only fixed point on the phase-plane will be a stable node (possibly degenerate or improper). Further, it is impossible for either eigenvalue to be 0 as both breakdown constants of GFP and NsrR are known to be explicitly positive (i.e. >0). Thus it is clear that the system always converges to a single equilibrium point. By finding the nullclines of the equation and working out where they meet, the fixed point can be found: | Since all the constants in this system of equations are constrained to be real semi-positive numbers (i.e. ≥0), both eigenvalues are always real and negative. This means that the only fixed point on the phase-plane will be a stable node (possibly degenerate or improper). Further, it is impossible for either eigenvalue to be 0 as both breakdown constants of GFP and NsrR are known to be explicitly positive (i.e. >0). Thus it is clear that the system always converges to a single equilibrium point. By finding the nullclines of the equation and working out where they meet, the fixed point can be found: | ||
[[Image:BCCS_GRN_EQ_10.jpg|frameless|center|upright=4|BCCS-Bristol GRN equation 10]] | [[Image:BCCS_GRN_EQ_10.jpg|frameless|center|upright=4|BCCS-Bristol GRN equation 10]] | ||
Revision as of 13:59, 24 October 2010
iGEM 2010
Contents |
Derivation
GRN Model
A Gene Regulatory Network (GRN) model is an established way of representing the biochemical reactions within a bacterial cell (7). Representing a network of interacting components using differential equations can be very useful, mainly because one can make analytical and qualitative assessments of a system’s behaviour with very limited information on parameters. This GRN represents a protein called NsrR that represses the production of mRNA encoding GFP, and so indirectly repressing the production of GFP. NsrR can bind to nitrates, effectively destroying its ability to repress the system.
Analysis
The most interesting relationship in this system of equations is the relationship between the concentration of NsrR repressor protein and the concentration of GFP in the system. A useful way of visualising and analysing the system’s dynamics is using a phase plane portrait in 2 dimensions. However, this means eliminating one of the variables. This can be done by making an equilibrium assumption about the amount of mRNA encoding GFP (mGFP for short). To do this, we set Equation 2 equal to 0. Rearranging in terms of [mGFP], this gives:
This can be substituted into Equation 3 to give an expression of [GFP] in terms of only itself and [NsrR]. This is now a system of 2 variables, so with a little further manipulation it can be put into a phase-plane portrait. The system of 2 linear ODEs is:
This can be written in matrix form (Jacobian) to help classify the dynamics of the system. The eigenvalues of the matrix form, M, determine the system’s dynamic states.
The eigenvalues of this matrix are:
Fixed Point Solution
Since all the constants in this system of equations are constrained to be real semi-positive numbers (i.e. ≥0), both eigenvalues are always real and negative. This means that the only fixed point on the phase-plane will be a stable node (possibly degenerate or improper). Further, it is impossible for either eigenvalue to be 0 as both breakdown constants of GFP and NsrR are known to be explicitly positive (i.e. >0). Thus it is clear that the system always converges to a single equilibrium point. By finding the nullclines of the equation and working out where they meet, the fixed point can be found:
