"
Team:Lethbridge Canada/math
From 2013hs.igem.org
| Line 54: | Line 54: | ||
<div id="math_contents"> | <div id="math_contents"> | ||
| - | <h1>Contents</h1> | + | <h1>Contents:</h1> |
<ul> | <ul> | ||
<li><a href="#the_goal">The Goal</a></li> | <li><a href="#the_goal">The Goal</a></li> | ||
| Line 81: | Line 81: | ||
<div id="the_goal"> | <div id="the_goal"> | ||
| - | <h1>The Goal</h1> | + | <h1>The Goal:</h1> |
<p>By mathematically modeling the maximum rate at which DH5a can produce proteins, the correct ratio of enzyme to protein can be determined. One group is cells are built to maximally express the gene of the oxytocin-neurophysin1 protein. With respect to protein levels, the amount of enzyme required to completely cut each oxytocin-neurophysin1 can be determined. The aim of this is to optimize the expression process to have exactly the amount of enzyme required and not produce an excess or too little. It is important to note that the model determines the theoretical maximum of protein production, not the experimental value.</p> | <p>By mathematically modeling the maximum rate at which DH5a can produce proteins, the correct ratio of enzyme to protein can be determined. One group is cells are built to maximally express the gene of the oxytocin-neurophysin1 protein. With respect to protein levels, the amount of enzyme required to completely cut each oxytocin-neurophysin1 can be determined. The aim of this is to optimize the expression process to have exactly the amount of enzyme required and not produce an excess or too little. It is important to note that the model determines the theoretical maximum of protein production, not the experimental value.</p> | ||
<a href="#header" class="math_back_to_top">Back to top</a> | <a href="#header" class="math_back_to_top">Back to top</a> | ||
| Line 89: | Line 89: | ||
<div id="bacterial_growth"> | <div id="bacterial_growth"> | ||
| - | <h1>Bacterial Growth</h1> | + | <h1>Bacterial Growth:</h1> |
<p>The purpose of including bacterial growth is so the amount of proteins produced per individual cell can be multiplied by the number of DH5a cells present to calculate the total number of proteins produced in that time interval of 1 minute. In the bacterial growth expression, the two variables are represented as such:</p> | <p>The purpose of including bacterial growth is so the amount of proteins produced per individual cell can be multiplied by the number of DH5a cells present to calculate the total number of proteins produced in that time interval of 1 minute. In the bacterial growth expression, the two variables are represented as such:</p> | ||
| Line 106: | Line 106: | ||
<div id="mrna_production"> | <div id="mrna_production"> | ||
| - | <h1>mRNA Production</h1> | + | <h1>mRNA Production:</h1> |
<div id="transcription"> | <div id="transcription"> | ||
| Line 125: | Line 125: | ||
<div id="mrna_decay"> | <div id="mrna_decay"> | ||
| - | <h2>mRNA Decay</h2> | + | <h2>mRNA Decay:</h2> |
<p>The stability of mRNA varies, but it will always decay. Half-lives can rage from seconds to hours.</p> | <p>The stability of mRNA varies, but it will always decay. Half-lives can rage from seconds to hours.</p> | ||
| Line 147: | Line 147: | ||
<div id="transcription_mrna_decay"> | <div id="transcription_mrna_decay"> | ||
| - | <h2>Transcription and mRNA Decay Working Together</h2> | + | <h2>Transcription and mRNA Decay Working Together:</h2> |
<p>When combining the previous two expressions, we can determine the amount of mRNA produced by one cell per minute.</p> | <p>When combining the previous two expressions, we can determine the amount of mRNA produced by one cell per minute.</p> | ||
| Line 161: | Line 161: | ||
<div id="total_mrna"> | <div id="total_mrna"> | ||
| - | <h1>Total mRNA</h1> | + | <h1>Total mRNA:</h1> |
<p>Expanding on the previous expression, we get the amount of total amount of mRNA produced by the total number of cells existing at any given time. This is achieved by multiplying the amount of mRNA produced by a single cell by the amount of cells that currently exist.</p> | <p>Expanding on the previous expression, we get the amount of total amount of mRNA produced by the total number of cells existing at any given time. This is achieved by multiplying the amount of mRNA produced by a single cell by the amount of cells that currently exist.</p> | ||
| Line 174: | Line 174: | ||
<div id="translation"> | <div id="translation"> | ||
| - | <h1>Translation</h1> | + | <h1>Translation:</h1> |
<p>The following expression models the maximum amount of proteins produced per minute.</p> | <p>The following expression models the maximum amount of proteins produced per minute.</p> | ||
| Line 189: | Line 189: | ||
<div id="derive_translation"> | <div id="derive_translation"> | ||
| - | <h2>Deriving the Translational Constant</h2> | + | <h2>Deriving the Translational Constant:</h2> |
<p>By dividing the expression for translation by the expression for transcription, we can determine the rate at which one strand of mRNA can be translated in to one protein.</p> | <p>By dividing the expression for translation by the expression for transcription, we can determine the rate at which one strand of mRNA can be translated in to one protein.</p> | ||
| Line 217: | Line 217: | ||
<div id="protein_production"> | <div id="protein_production"> | ||
| - | <h1>Protein Production</h1> | + | <h1>Protein Production:</h1> |
<p>By taking the expression from Total mRNA section and incorporating the Translational Constant section we can determine the total protein production. By multiplying the amount of mRNA produced by the rate at which mRNA is translated in to protein, we can calculate the amount of protein produced. We do this by finding the area under the curve with respect to time.</p> | <p>By taking the expression from Total mRNA section and incorporating the Translational Constant section we can determine the total protein production. By multiplying the amount of mRNA produced by the rate at which mRNA is translated in to protein, we can calculate the amount of protein produced. We do this by finding the area under the curve with respect to time.</p> | ||
| Line 230: | Line 230: | ||
<div id="equation"> | <div id="equation"> | ||
| - | <h1>Assembling the Equation</h1> | + | <h1>Assembling the Equation:</h1> |
<p>Bacteria don’t double forever. They will reach a maximum population for their environment size. For this reason, to accurately model protein production we need to make a piece-wise equation. We let <i>j</i> equal the amount of time that it takes for DH5a to reach their maximum population size. Before time j we use the expression for start time <i>i</i> up to time <i>j</i>. Because after time <i>j</i>, DH5a no longer double, we replace the time variable with <i>j</i> so the graph from this point will be a horizontal line parallel to the time-axis. Time <i>k</i> is the end time of the protein production period. The area under the curve from <i>i</i> to <i>k</i> is the total amount of protein expressed (<i>n<sub>p</sub></i>).</p> | <p>Bacteria don’t double forever. They will reach a maximum population for their environment size. For this reason, to accurately model protein production we need to make a piece-wise equation. We let <i>j</i> equal the amount of time that it takes for DH5a to reach their maximum population size. Before time j we use the expression for start time <i>i</i> up to time <i>j</i>. Because after time <i>j</i>, DH5a no longer double, we replace the time variable with <i>j</i> so the graph from this point will be a horizontal line parallel to the time-axis. Time <i>k</i> is the end time of the protein production period. The area under the curve from <i>i</i> to <i>k</i> is the total amount of protein expressed (<i>n<sub>p</sub></i>).</p> | ||
| Line 243: | Line 243: | ||
<div id="values"> | <div id="values"> | ||
| - | <h1>Values</h1> | + | <h1>Values:</h1> |
<a href="#header" class="math_back_to_top">Back to top</a> | <a href="#header" class="math_back_to_top">Back to top</a> | ||
</div> | </div> | ||
| Line 250: | Line 250: | ||
<div id="references"> | <div id="references"> | ||
| - | <h1>References</h1> | + | <h1>References:</h1> |
<ul> | <ul> | ||
<li><a href="https://2012hs.igem.org/Team:CIDEB-UANL_Mexico/Math/Overview">https://2012hs.igem.org/Team:CIDEB-UANL_Mexico/Math/Overview</a></li> | <li><a href="https://2012hs.igem.org/Team:CIDEB-UANL_Mexico/Math/Overview">https://2012hs.igem.org/Team:CIDEB-UANL_Mexico/Math/Overview</a></li> | ||
Revision as of 20:34, 12 June 2013
Contents:
The Goal:
By mathematically modeling the maximum rate at which DH5a can produce proteins, the correct ratio of enzyme to protein can be determined. One group is cells are built to maximally express the gene of the oxytocin-neurophysin1 protein. With respect to protein levels, the amount of enzyme required to completely cut each oxytocin-neurophysin1 can be determined. The aim of this is to optimize the expression process to have exactly the amount of enzyme required and not produce an excess or too little. It is important to note that the model determines the theoretical maximum of protein production, not the experimental value.
Back to topBacterial Growth:
The purpose of including bacterial growth is so the amount of proteins produced per individual cell can be multiplied by the number of DH5a cells present to calculate the total number of proteins produced in that time interval of 1 minute. In the bacterial growth expression, the two variables are represented as such:
- bi = initial number of DH5a bacteria present
- t = time (manipulated variable)
We assumed a doubling time of 30 minutes.
Back to topmRNA Production:
Transcription
The expression below calculates the amount of mRNA that is produced by an individual cell in the time period of 1 minute. Credit is given to the 2009 Beijing iGEM team for the average speed of transcription in E. coli.
- 4200nt/min = the maximum rate at which DNA is transcribed in to mRNA in E. coli
- lgene = gene length in (from RBS to the double terminator)
The resultant number is the number of mRNA strands produced per minute in one cell.
Back to topmRNA Decay:
The stability of mRNA varies, but it will always decay. Half-lives can rage from seconds to hours.
- ts = stagnant time of mRNA between transcription and translation
- th = half-life
The above expression will act as a coefficient for how much mRNA will remain after the amount of time that mRNA sits stagnant in the cells between the processes of transcription and translation. However, we have assumed that immediately or very shortly after transcription, the mRNA then undergoes translation. The following logic describes this situation:
Because the exponent of the function is approximately or equal to 0, the entire expression for decay is equal to 1. The numerator of the exponent is essentially 0, so nearly any denominator would give the entire exponent a value of 0 when simplified. Although this logic nullifies the purpose of this term, it is important to note that mRNA decay is existent.
Back to topTranscription and mRNA Decay Working Together:
When combining the previous two expressions, we can determine the amount of mRNA produced by one cell per minute.
Total mRNA:
Expanding on the previous expression, we get the amount of total amount of mRNA produced by the total number of cells existing at any given time. This is achieved by multiplying the amount of mRNA produced by a single cell by the amount of cells that currently exist.
Translation:
The following expression models the maximum amount of proteins produced per minute.
- 2400Aa/min = the rate at which DNA is transcribed in to mRNA in E. coli
- lprotein = length of the protein to be expressed (Aa)
Credit is given to the 2009 Beijing iGEM team for the average speed of translation in E. coli.
Back to topDeriving the Translational Constant:
By dividing the expression for translation by the expression for transcription, we can determine the rate at which one strand of mRNA can be translated in to one protein.
Because there are three nucleotides (also called a codon) of mRNA for every one amino acid, the gene length (nt) over the protein length (Aa) will always be equal to 3.
By simplifying this expression, we get the following term that is the rate of translation from mRNA.
We used the RBS BBa_B0034 which has an efficiency index of 1.0.
Protein Production:
By taking the expression from Total mRNA section and incorporating the Translational Constant section we can determine the total protein production. By multiplying the amount of mRNA produced by the rate at which mRNA is translated in to protein, we can calculate the amount of protein produced. We do this by finding the area under the curve with respect to time.
Assembling the Equation:
Bacteria don’t double forever. They will reach a maximum population for their environment size. For this reason, to accurately model protein production we need to make a piece-wise equation. We let j equal the amount of time that it takes for DH5a to reach their maximum population size. Before time j we use the expression for start time i up to time j. Because after time j, DH5a no longer double, we replace the time variable with j so the graph from this point will be a horizontal line parallel to the time-axis. Time k is the end time of the protein production period. The area under the curve from i to k is the total amount of protein expressed (np).
