|
|
(161 intermediate revisions not shown) |
Line 1: |
Line 1: |
- | = Introduction =
| + | {{Team:Amsterdam/stylesheet}} |
| + | {{Team:Amsterdam/scripts}} |
| + | {{Team:Amsterdam/Header}} |
| + | {{Team:Amsterdam/Sidebar1}} |
| | | |
- | We’ve established in a very early stage of the design process that we were for growth and degradation speed of plasmid going to use DNA methylation as the molecular mechanism to create our storage mechanism. DNA methylation, usually part of the Type II restriction system, is performed by a group of proteins called the methyltransferases. In choosing DNA-methylation, we were heavily inspired by the DamID technology which has been originated by the Bas van Steensel . By fusing a bacterial methyltransferase (MTase) to a eukaryotic transcription factor which is to be characterized, van Steensel was able to infer the transcription factor binding sites on the eukaryotic genome by reading out where methylation had occurred. We thought we could reverse this idea to create memory unit, at the core constituted by an MTase that would only methylate an especially designated genomic region in reaction to the sensing of a signal. For this, we required an MTase that
| + | <div id="content-area"> |
| + | <div id="sub-menu" class="content-block"> |
| | | |
- | <blockquote> | + | <center> |
- | is not already present in the chassis organism ''E.coli'',
| + | <html> |
| + | <embed src="https://static.igem.org/mediawiki/2012/2/20/PoC1v4.swf" width="500px" height="400px"> |
| + | </html> |
| + | </center> |
| + | == The molecular design in a nutshell == |
| | | |
| + | We established in a very early stage of the design process to use DNA methylation as the molecular mechanism to create our storage mechanism. DNA methylation, which means enriching a specific nucleotide with a methyl group, is performed by a group of proteins called the methyltransferases. In making this design choice, we were heavily influenced by the DamID technology ([[#Steensel1|4]],[[#Steensel2|5]]). Van Steensel was interested in the binding sites in drosophila genomes of various eukaryotic transcription factors. By fusing a bacterial methyltransferase (MTase) to these transcription factors, he was able to infer the transcription factor binding sites by reading out which genomic regions were methylated. |
| | | |
- | has a binding motif of that is not methylated by any other MTases in the system and for which
| + | We thought we could reverse this idea to create a memory unit, at the core constituted by an MTase that would methylate an especially designated genomic region but only in reaction to the sensing of a signal by the cell. The epigenetic status of the bit region could thus either be 1 (methylated or written) ''or'' 0 (ummethylated or unwritten), effectively forming a binary memory unit. |
| | | |
- | | + | Finally, the methylation status of the bit region can be assessed using digestion of the plasmid extracts with the MTase-coupled RE, followed by analyzing the band lengths intensities of a gel electrophoresis of the product. |
- | a restriction enzyme (RE) dependent upon the action of this methyltransferase has been identified.
| + | |
- | </blockquote>
| + | |
- | By finding a way of controlling the activity of the MTase, a genomic region purposed to be methylated, named the ''bit region'', would only be methylated when the cell would be induced to do so by the signal. Finally, the methylation status of the bit regions can then be assessed using digestion with this RE followed by a band length analysis of a gel electrophoresis of the product. This is process is covered more in depth in section {sub:gelinfer}
| + | |
| | | |
| == Choosing a methyltransferase == | | == Choosing a methyltransferase == |
| | | |
- | We wrote a python script to mine the [http://rebase.neb.com/rebase/rebase.html REBASE] database of REs and MTases for a MTase that meets our needs. We found an ideal candidate: M.ScaI, originally from ''Streptomyces caespitosus''. Now, the only remaining design question was how to control the activity of the fusion protein as we would want it to only be active in the presence of a signal. We could either try to control its ''presence'' in the cell by induced transcription or try to directly fuse it to a inducible transcriptional regulator, hoping that the blocked activity of the latter in absence of a signal would also prohibit the MTase from methylating the bit. This design question will be more elaborately discussed in section {chapter_molecular}. | + | To create this system, we first required an MTase that: i) is not already present in the chassis organism ''E.coli'', ii) has a binding motif that is not methylated by any other MTases in the system and iii) for which a restriction enzyme (RE) dependent upon the action of this methyltransferase has been identified. This last condition is always true for any methyltransferase which has been identified to be part of a bacterial [http://en.wikipedia.org/wiki/Restriction_modification_system Restriction/Modification system]. |
| + | We wrote a python script to mine the [http://rebase.neb.com/rebase/rebase.html REBASE] database of REs and MTases for a MTase that meets our needs. We found an ideal candidate: M.ScaI, originally from ''Streptomyces caespitosus''. |
| + | The only remaining design question was how to control the activity of the fusion protein as we would want it to only be active in the presence of the generic signal. This design question will be more elaborately discussed in the next section. |
| | | |
| == Extending the idea to multiple bits == | | == Extending the idea to multiple bits == |
| + | Extending the design idea further, we realized that it should be possible to register and store the presence of multiple signals smultaneously in a single cell. This extension can be accomplished by fusing DNA binding proteins to the methyltransferase and adding their corresponding DNA binding motifs to the bit regions of the DNA. Either traditional bacterial transcriptional regulators could be used for this purpose or the Zinc-Finger Array (ZFA) technology ([[#Fu|1]]). ZFA allows the construction of highly specific artificial protein-DNA interactions ([[#Kaseniit|2]]). |
| | | |
- | Extending the design idea further, we realized that it should be possible to store the presence of different signals – multiple ''bits'' – simultaneously by fusing DNA binding proteins to the methyltransferase and adding their corresponding DNA binding motifs to the bit regions of the DNA. Either traditional bacterial transcriptional regulators could be used for this purpose or the more extendable Zinc-Finger Array (ZFA) technology . The latter allows the construction of highly specific proteins to custom DNA motifs. The large increase in binding affinity of the resulting MTase-ZF fusion protein would have to make sure that the same MTase can be used for different signals with a very low chance of cross-talk between FPs and other signals’ bits.
| + | By fusing the MTase to a Zinc Finger or traditional transcriptional regulator and extending the bit region with the transcriptional regulator’s binding motif, the binding affinity of the fusion protein for the DNA motif could increase about 18-fold ([[#McNamara|3]]) (a C5-cytosine MTase was studied here, not in the class of N4C-MTases to which M.ScaI belongs). This will heavily diminish the possible aspecific cross-talk between FPs and bits assigned to other FPs. |
- | | + | |
- | = Inferring the time of signal onset =
| + | |
- | | + | |
- | Numerous copies of identical plasmids are often present in single cells and plasmids replicate independently of the bacterial chromosome . A plasmid copy number (PCN) has been determined for all plasmids in the Parts Registry, which indicates a likely amount of copies of the plasmid to be present in each cell. Unlike eukaryotes, we discovered that DNA methylation patterns are not copied to the new strand during DNA replication in prokaryotes. This will have to lead to a dilution of the amount of plasmids with a written bit over time, mostly due to cell replication and the ensuing binomial division of the plasmids in the parent cell among the two daughter cells. This seemed like a downside to using a MTase system at first, but it quickly opened our eyes to a very exciting application of this system. Assuming nearly constant proliferation & replication rates for the cells and the within the plasmids, this dilution can be used to infer how long ago the signal of interest was registered by the cell.
| + | |
- | | + | |
- | Using the here presented model, we will examine how this phenomenon can exactly be used to infer the registration time of the signal by the cell. Another goal of this model is to gain insight into how much the various parameters for growth and degradation speed of plasmid in the model affect the ratio <math>F(t) = \frac{\text{written plasmids}(t)}{\text{total plasmids}}</math>. Unknown variables affecting <math>F(t)</math> in a real-life setting would be the time of signal onset, signal duration and signal strength. Additionally the signal response, degradation and proliferation rates of the ''Cellular Logbook'' will influence <math>F(t)</math>, but these should be experimentally determinable before actual deployment and application of our system. This way the can be used to measure a signal and estimate its time of onset. This would allow the inferrance of signal onset <math>s_{\text{on}}</math> from the fraction of methylated over total bits in the population. As this model has the sole purpose of clarification of the desired system behaviour, the units in this model are completely arbitrary.
| + | |
- | | + | |
- | First, let’s model the input signal/compound which is to be registered. Let <math>s_{\text{on}}</math> be the time at which the signal is turned on and <math>s_{\text{off}}</math> be the time at which the signal is turned off. Imagine the to be stationary and positioned along a fluidic stream with the signal to be registered passing the at <math>t = s_{\text{on}}</math> and leaving it at <math>t = s_{\text{off}}</math>. Modelling the signal using the piecewise function <math>S(t)</math> in equation {eqn:signal} seems appropriate. This function has been plotted in Figure {fig:signal}.
| + | |
- | | + | |
- | <math>S(t) = \left\{
| + | |
- | \begin{array}{lr}
| + | |
- | 0 & 0 \le t \le s_{\text{on}} \\
| + | |
- | 1 & s_{\text{on}} \le t \le s_{\text{off}} \\
| + | |
- | 0 & s_{\text{off}} \le t \le \infty
| + | |
- | \end{array}
| + | |
- | \right.
| + | |
- | \label{eqn:signal}
| + | |
- | </math>
| + | |
- | | + | |
- | <blockquote>[[Image:signal.png|image]] {fig:signal} { Plot of the input signal <math>S(t)</math> with <math>s_{\text{on}}</math> at 3 and <math>s_{{\text{off}}}</math> at 4}
| + | |
- | </blockquote>
| + | |
- | A single cell with multiple identical plasmids on which the gene for the methyltransferase and a so called ''bit region'', which is the region especially purposed to be methylated in presence of a signal, have been placed. In this cell the fraction of methylated plasmids <math>F(t)</math> will then be in direct correspondence to how long ago the signal has been registered if the following plausible rules are true:
| + | |
- | | + | |
- | * A well stirred cellular system in which no spatial concentrations differences exist and in which all species concentrations are large enough to be approximated by continuous functions. Such a system is describable by a set of ordinary differential equations (ODE’s).
| + | |
- | * To easy the analysis of the model, a single unique bit with multiple copies per cell is considered here. The system is easily upgradeable to multiple signals however.
| + | |
- | * Assumed is a high response rate ( 3 min) termed <math>\omega</math>, which is the constant rate with which the system responds with methylation to the sensing of a signal.
| + | |
- | * We will assume logistic growth for the memory cell population, with a capacity limit of <math>Ca = 5000</math>.
| + | |
- | * Alternatively, this maximal population number could also be described as the fraction between the cellular proliferation rate (<math>\beta</math>) and degradation rate (<math>\alpha</math>), as the steady state amount of cells in the population will be determined by <math>\frac{\beta}{\alpha}</math> as a solution to the differential equation <math>P'(t) = \beta P - \alpha P</math>.
| + | |
- | * Accumulation of cells due to the division of cells in which te bit has been written is assumed to result in non-written cells; methylation patterns are not copied to the progeny in prokaryotes.
| + | |
- | * <math>P_{0}</math> denote cells that have the single bit set to 0, no write event has taken place in these cells.
| + | |
- | * <math>P_{1}</math> indicate cells in which the bit has been flipped to 1 in response to encountering the signal compound.
| + | |
- | * <math>P_{0+1}</math> indicates the total amount of cells in the sytem, <math>P_{0} + P_{1}</math>
| + | |
- | | + | |
- | From these rules, the following system of ODE’s has been constructed:
| + | |
- | | + | |
- | <math>\frac{dP_{0}}{dt} = k\ P_{0+1}\ (1 - \frac{P_{0+1}}{\text{Ca}}) - \omega\ S(t)\ P_{0} - \alpha\ P_{0} \\
| + | |
- | \frac{dP_{1}}{dt} = \omega\ S(t)\ P_{0} - \alpha\ P_{1}
| + | |
- | </math>
| + | |
- | | + | |
- | <blockquote>
| + | |
- | <table>
| + | |
- | <thead>
| + | |
- | <tr class="header">
| + | |
- | <th align="left">Parameter</th>
| + | |
- | <th align="right">Value</th>
| + | |
- | </tr>
| + | |
- | </thead>
| + | |
- | <tbody>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left">k</td>
| + | |
- | <td align="right">1.5</td>
| + | |
- | </tr>
| + | |
- | <tr class="even">
| + | |
- | <td align="left"><math>\beta</math></td>
| + | |
- | <td align="right">2500</td>
| + | |
- | </tr>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left"><math>\alpha</math></td>
| + | |
- | <td align="right">.5</td>
| + | |
- | </tr>
| + | |
- | <tr class="even">
| + | |
- | <td align="left">Ca</td>
| + | |
- | <td align="right">5000</td>
| + | |
- | </tr>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left"><math>\omega</math></td>
| + | |
- | <td align="right">8</td>
| + | |
- | </tr>
| + | |
- | </tbody>
| + | |
- | </table>
| + | |
- | | + | |
- | { Parameter values for the plasmid methylation model} {tab:populationparms}
| + | |
- | </blockquote>
| + | |
- | Using the parameter values of Table {tab:populationparms} a simulation with a duration of 20 time units is shown in Figure {fig:timelapse}.
| + | |
- | | + | |
- | <blockquote>[[Image:1|image]] { Input signal <math>S(t)</math> with <math>s_{\text{on}}</math> at 3 and <math>s_{{\text{off}}}</math> at 4} {fig:timelapse}
| + | |
- | </blockquote>
| + | |
- | <blockquote>[[Image:2|image]] { Sum of unmethylated and methylated cells, <math>S_{0+1}</math>} {fig:totalcells}
| + | |
- | </blockquote>
| + | |
- | The plasmid population within a is shown to be completely converted to methylated plasmids shortly after . As long as the signal is still present – until <math>s_{\text{off}}</math>, – the bit on all newly copied plasmids will be immediately methylated as the signal is still present. After <math>s_{\text{off}}</math>, the percentage of <math>\frac{\text{methylated} \text{cells}}{\text{total} \text{cells}}</math> will be diluted, due to cell proliferation implicating plasmid division between daughter cells and cell death. The size of the time span with which a hint of methylated cells is still present therefore relates to two factors: positively with the amount of methylated cells at <math>s_{\text{off}}</math> and negatively with the cellular death rate.
| + | |
- | | + | |
- | The bit writing speed termed <math>\omega</math>, or the rate with which the presence of a signal causes the methylation of a bit, is likely to be very high in comparison to the cellular division rate. The underlying required processes such as the diffusion of the signal across the cellular membrane, the signal inducing gene transcription of the FP, mRNA translation and transcription & and FP diffusion over to the bit all occur well within 5-minute time scales. Every single gene on a plasmid will then be expected to be methylated within at most 5 minutes of registering of the signal.
| + | |
- | | + | |
- | The assumption that the population will have reached its steady state level (SS), given by the production rate (<math>\beta</math>) divided by the degradation rate (<math>\alpha</math>), before is a plausible one. The ''Cellular Logbook'' population will have been stored inside of a vesicle consisting of a non-bacterial permeating membrane with a very small volume. Due to this volume and the rapid bacterial division time (<math>\sim</math> 40 mins), the spatial growth limits of the system will then have been reached in a short amount of time. The steady state-level of the bacterial population has already been reached before the signal is registered by the population.
| + | |
- | | + | |
- | To reinforce that the following equation always holds in this model; <math>P_{0} + P_{1} = P_{0+1}</math> the following plot is shown, in which the sum of the methylated and unmethylated cells over time is plotted. This clearly shows the limit, specified by the Capacity limit, which is reached around <math>t=6</math>.
| + | |
- | | + | |
- | == Methylated plasmids over time ==
| + | |
- | | + | |
- | {sub:gelinfer} The monotonically decreasing value of <math>F(t) = \frac{\text{methylated plasmids}}{\text{total plasmids}}</math> can be used to infer , given that the cellular proliferation rate (<math>\beta</math>) and degradation rate (<math>\alpha</math>) are known and constant and the assuming that all bits are methylated during the presence of the signal. Irrespective of the initial amount of plasmids, the population of plasmids within the single cell will have reached a steady state value of <math>\frac{\beta}{\alpha}</math>. After the signal has leaved the medium in which is present, will start to decrease as a function of the degradation rate .
| + | |
- | | + | |
- | <math>\frac{dP_{1}}{dt} = - \alpha\ P_{1}
| + | |
- | </math>
| + | |
- | | + | |
- | Integrating this differential equation, <math>P_{1}</math> will be given by: <math>P_{1}(t) = e^{-\alpha t}
| + | |
- | </math> And multiplying by the steady value <math>\frac{\beta}{\alpha}</math> will yield the amount of methylated plasmids at time <math>t</math>, given that there were <math>\frac{\beta}{\alpha}</math> methylated plasmids at <math>t = 0</math>. <math>P_{1}(t) = \frac{\beta}{\alpha} e^{-\alpha t}
| + | |
- | </math>
| + | |
- | | + | |
- | == In practice ==
| + | |
- | | + | |
- | <blockquote>[[Image:bands|image]] { Gel representations for a range of different <math>F(t)</math> values. Complete methylation of all bits results in a single, bright band at the height of the length of the linearized plasmid. Decreasing the amount of methylated bits shifts the intensity of the top band away to the two bottom bands. These indicate the linearized & digested plasmid.} {fig:gelarray}
| + | |
- | </blockquote>
| + | |
- | In a typical laboratory situation, doing a restriction enzyme assay on the miniprep-extracted plasmid DNA out of followed by gel electrophoresis will be the most convenient way to assess the methylation status of the bits. The relative intensities of the gel bands can then be used to infer , which is exemplified in Figure {fig:gelarray}. Unmethylated bits will result in successfully digested DNA fragments and thus two bands of shorter DNA fragments. Methylated bits will not be cut and will therefore result in one longer band, shown more to the top of the gel. Hence the top band and the two bottom ones are mutually exclusive as they indicate the same (linearized) plasmid DNA to either be cut, resulting in the two bottom bands, or not cut, resulting in the top band.
| + | |
- | | + | |
- | The function for <math>P_{1}(t)</math> is shown in in Figure {fig:readout}. The rate with which the function approaches its limiting value of 0 is directly dependent upon the degradation rate . To get a hands-on feel of the effects on the system an interactive version of this plot, in which , and are controllable via a GUI is included in the attached Mathematica file.
| + | |
- | | + | |
- | <blockquote>[[Image:infer|image]] { Animation of <math>P(t)</math> with its value ranging between 0.01 and 1. Solving for <math>t</math>, the time after which the signal has left/been degraded in the Cell, the solution is shown by the blue dot and the text underneath the plot. The value of <math>P(t)</math> will decrease over time until it becomes 0. An interactive version of this plot has been included in the attached Mathematica file.} {fig:readout}
| + | |
- | </blockquote>
| + | |
- | = Introduction =
| + | |
- | | + | |
- | We’ve established in a very early stage of the design process that we were for growth and degradation speed of plasmid going to use DNA methylation as the molecular mechanism to create our storage mechanism. DNA methylation, usually part of the Type II restriction system, is performed by a group of proteins called the methyltransferases. In choosing DNA-methylation, we were heavily inspired by the DamID technology which has been originated by the Bas van Steensel . By fusing a bacterial methyltransferase (MTase) to a eukaryotic transcription factor which is to be characterized, van Steensel was able to infer the transcription factor binding sites on the eukaryotic genome by reading out where methylation had occurred. We thought we could reverse this idea to create memory unit, at the core constituted by an MTase that would only methylate an especially designated genomic region in reaction to the sensing of a signal. For this, we required an MTase that
| + | |
- | | + | |
- | <blockquote>
| + | |
- | is not already present in the chassis organism ''E.coli'',
| + | |
- | | + | |
- | | + | |
- | has a binding motif of that is not methylated by any other MTases in the system and for which
| + | |
- | | + | |
- | | + | |
- | a restriction enzyme (RE) dependent upon the action of this methyltransferase has been identified.
| + | |
- | </blockquote>
| + | |
- | By finding a way of controlling the activity of the MTase, a genomic region purposed to be methylated, named the ''bit region'', would only be methylated when the cell would be induced to do so by the signal. Finally, the methylation status of the bit regions can then be assessed using digestion with this RE followed by a band length analysis of a gel electrophoresis of the product. This is process is covered more in depth in section {sub:gelinfer}
| + | |
- | | + | |
- | == Choosing a methyltransferase ==
| + | |
- | | + | |
- | We wrote a python script to mine the [http://rebase.neb.com/rebase/rebase.html REBASE] database of REs and MTases for a MTase that meets our needs. We found an ideal candidate: M.ScaI, originally from ''Streptomyces caespitosus''. Now, the only remaining design question was how to control the activity of the fusion protein as we would want it to only be active in the presence of a signal. We could either try to control its ''presence'' in the cell by induced transcription or try to directly fuse it to a inducible transcriptional regulator, hoping that the blocked activity of the latter in absence of a signal would also prohibit the MTase from methylating the bit. This design question will be more elaborately discussed in section {chapter_molecular}.
| + | |
- | | + | |
- | == Extending the idea to multiple bits ==
| + | |
- | | + | |
- | Extending the design idea further, we realized that it should be possible to store the presence of different signals – multiple ''bits'' – simultaneously by fusing DNA binding proteins to the methyltransferase and adding their corresponding DNA binding motifs to the bit regions of the DNA. Either traditional bacterial transcriptional regulators could be used for this purpose or the more extendable Zinc-Finger Array (ZFA) technology . The latter allows the construction of highly specific proteins to custom DNA motifs. The large increase in binding affinity of the resulting MTase-ZF fusion protein would have to make sure that the same MTase can be used for different signals with a very low chance of cross-talk between FPs and other signals’ bits.
| + | |
- | | + | |
- | = Inferring the time of signal onset =
| + | |
- | | + | |
- | Numerous copies of identical plasmids are often present in single cells and plasmids replicate independently of the bacterial chromosome . A plasmid copy number (PCN) has been determined for all plasmids in the Parts Registry, which indicates a likely amount of copies of the plasmid to be present in each cell. Unlike eukaryotes, we discovered that DNA methylation patterns are not copied to the new strand during DNA replication in prokaryotes. This will have to lead to a dilution of the amount of plasmids with a written bit over time, mostly due to cell replication and the ensuing binomial division of the plasmids in the parent cell among the two daughter cells. This seemed like a downside to using a MTase system at first, but it quickly opened our eyes to a very exciting application of this system. Assuming nearly constant proliferation & replication rates for the cells and the within the plasmids, this dilution can be used to infer how long ago the signal of interest was registered by the cell.
| + | |
- | | + | |
- | Using the here presented model, we will examine how this phenomenon can exactly be used to infer the registration time of the signal by the cell. Another goal of this model is to gain insight into how much the various parameters for growth and degradation speed of plasmid in the model affect the ratio <math>F(t) = \frac{\text{written plasmids}(t)}{\text{total plasmids}}</math>. Unknown variables affecting <math>F(t)</math> in a real-life setting would be the time of signal onset, signal duration and signal strength. Additionally the signal response, degradation and proliferation rates of the ''Cellular Logbook'' will influence <math>F(t)</math>, but these should be experimentally determinable before actual deployment and application of our system. This way the can be used to measure a signal and estimate its time of onset. This would allow the inferrance of signal onset <math>s_{\text{on}}</math> from the fraction of methylated over total bits in the population. As this model has the sole purpose of clarification of the desired system behaviour, the units in this model are completely arbitrary.
| + | |
- | | + | |
- | First, let’s model the input signal/compound which is to be registered. Let <math>s_{\text{on}}</math> be the time at which the signal is turned on and <math>s_{\text{off}}</math> be the time at which the signal is turned off. Imagine the to be stationary and positioned along a fluidic stream with the signal to be registered passing the at <math>t = s_{\text{on}}</math> and leaving it at <math>t = s_{\text{off}}</math>. Modelling the signal using the piecewise function <math>S(t)</math> in equation {eqn:signal} seems appropriate. This function has been plotted in Figure {fig:signal}.
| + | |
- | | + | |
- | <math>S(t) = \left\{
| + | |
- | \begin{array}{lr}
| + | |
- | 0 & 0 \le t \le s_{\text{on}} \\
| + | |
- | 1 & s_{\text{on}} \le t \le s_{\text{off}} \\
| + | |
- | 0 & s_{\text{off}} \le t \le \infty
| + | |
- | \end{array}
| + | |
- | \right.
| + | |
- | \label{eqn:signal}
| + | |
- | </math>
| + | |
- | | + | |
- | <blockquote>[[Image:signal.png|image]] {fig:signal} { Plot of the input signal <math>S(t)</math> with <math>s_{\text{on}}</math> at 3 and <math>s_{{\text{off}}}</math> at 4}
| + | |
- | </blockquote>
| + | |
- | A single cell with multiple identical plasmids on which the gene for the methyltransferase and a so called ''bit region'', which is the region especially purposed to be methylated in presence of a signal, have been placed. In this cell the fraction of methylated plasmids <math>F(t)</math> will then be in direct correspondence to how long ago the signal has been registered if the following plausible rules are true:
| + | |
- | | + | |
- | * A well stirred cellular system in which no spatial concentrations differences exist and in which all species concentrations are large enough to be approximated by continuous functions. Such a system is describable by a set of ordinary differential equations (ODE’s).
| + | |
- | * To easy the analysis of the model, a single unique bit with multiple copies per cell is considered here. The system is easily upgradeable to multiple signals however.
| + | |
- | * Assumed is a high response rate ( 3 min) termed <math>\omega</math>, which is the constant rate with which the system responds with methylation to the sensing of a signal.
| + | |
- | * We will assume logistic growth for the memory cell population, with a capacity limit of <math>Ca = 5000</math>.
| + | |
- | * Alternatively, this maximal population number could also be described as the fraction between the cellular proliferation rate (<math>\beta</math>) and degradation rate (<math>\alpha</math>), as the steady state amount of cells in the population will be determined by <math>\frac{\beta}{\alpha}</math> as a solution to the differential equation <math>P'(t) = \beta P - \alpha P</math>.
| + | |
- | * Accumulation of cells due to the division of cells in which te bit has been written is assumed to result in non-written cells; methylation patterns are not copied to the progeny in prokaryotes.
| + | |
- | * <math>P_{0}</math> denote cells that have the single bit set to 0, no write event has taken place in these cells.
| + | |
- | * <math>P_{1}</math> indicate cells in which the bit has been flipped to 1 in response to encountering the signal compound.
| + | |
- | * <math>P_{0+1}</math> indicates the total amount of cells in the sytem, <math>P_{0} + P_{1}</math>
| + | |
- | | + | |
- | From these rules, the following system of ODE’s has been constructed:
| + | |
- | | + | |
- | <math>\frac{dP_{0}}{dt} = k\ P_{0+1}\ (1 - \frac{P_{0+1}}{\text{Ca}}) - \omega\ S(t)\ P_{0} - \alpha\ P_{0} \\
| + | |
- | \frac{dP_{1}}{dt} = \omega\ S(t)\ P_{0} - \alpha\ P_{1}
| + | |
- | </math>
| + | |
- | | + | |
- | <blockquote>
| + | |
- | <table>
| + | |
- | <thead>
| + | |
- | <tr class="header">
| + | |
- | <th align="left">Parameter</th>
| + | |
- | <th align="right">Value</th>
| + | |
- | </tr>
| + | |
- | </thead>
| + | |
- | <tbody>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left">k</td>
| + | |
- | <td align="right">1.5</td>
| + | |
- | </tr>
| + | |
- | <tr class="even">
| + | |
- | <td align="left"><math>\beta</math></td>
| + | |
- | <td align="right">2500</td>
| + | |
- | </tr>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left"><math>\alpha</math></td>
| + | |
- | <td align="right">.5</td>
| + | |
- | </tr>
| + | |
- | <tr class="even">
| + | |
- | <td align="left">Ca</td>
| + | |
- | <td align="right">5000</td>
| + | |
- | </tr>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left"><math>\omega</math></td>
| + | |
- | <td align="right">8</td>
| + | |
- | </tr>
| + | |
- | </tbody>
| + | |
- | </table>
| + | |
- | | + | |
- | { Parameter values for the plasmid methylation model} {tab:populationparms}
| + | |
- | </blockquote>
| + | |
- | Using the parameter values of Table {tab:populationparms} a simulation with a duration of 20 time units is shown in Figure {fig:timelapse}.
| + | |
- | | + | |
- | <blockquote>[[Image:1|image]] { Input signal <math>S(t)</math> with <math>s_{\text{on}}</math> at 3 and <math>s_{{\text{off}}}</math> at 4} {fig:timelapse}
| + | |
- | </blockquote>
| + | |
- | <blockquote>[[Image:2|image]] { Sum of unmethylated and methylated cells, <math>S_{0+1}</math>} {fig:totalcells}
| + | |
- | </blockquote>
| + | |
- | The plasmid population within a is shown to be completely converted to methylated plasmids shortly after . As long as the signal is still present – until <math>s_{\text{off}}</math>, – the bit on all newly copied plasmids will be immediately methylated as the signal is still present. After <math>s_{\text{off}}</math>, the percentage of <math>\frac{\text{methylated} \text{cells}}{\text{total} \text{cells}}</math> will be diluted, due to cell proliferation implicating plasmid division between daughter cells and cell death. The size of the time span with which a hint of methylated cells is still present therefore relates to two factors: positively with the amount of methylated cells at <math>s_{\text{off}}</math> and negatively with the cellular death rate.
| + | |
- | | + | |
- | The bit writing speed termed <math>\omega</math>, or the rate with which the presence of a signal causes the methylation of a bit, is likely to be very high in comparison to the cellular division rate. The underlying required processes such as the diffusion of the signal across the cellular membrane, the signal inducing gene transcription of the FP, mRNA translation and transcription & and FP diffusion over to the bit all occur well within 5-minute time scales. Every single gene on a plasmid will then be expected to be methylated within at most 5 minutes of registering of the signal.
| + | |
- | | + | |
- | The assumption that the population will have reached its steady state level (SS), given by the production rate (<math>\beta</math>) divided by the degradation rate (<math>\alpha</math>), before is a plausible one. The ''Cellular Logbook'' population will have been stored inside of a vesicle consisting of a non-bacterial permeating membrane with a very small volume. Due to this volume and the rapid bacterial division time (<math>\sim</math> 40 mins), the spatial growth limits of the system will then have been reached in a short amount of time. The steady state-level of the bacterial population has already been reached before the signal is registered by the population.
| + | |
- | | + | |
- | To reinforce that the following equation always holds in this model; <math>P_{0} + P_{1} = P_{0+1}</math> the following plot is shown, in which the sum of the methylated and unmethylated cells over time is plotted. This clearly shows the limit, specified by the Capacity limit, which is reached around <math>t=6</math>.
| + | |
- | | + | |
- | == Methylated plasmids over time ==
| + | |
- | | + | |
- | {sub:gelinfer} The monotonically decreasing value of <math>F(t) = \frac{\text{methylated plasmids}}{\text{total plasmids}}</math> can be used to infer , given that the cellular proliferation rate (<math>\beta</math>) and degradation rate (<math>\alpha</math>) are known and constant and the assuming that all bits are methylated during the presence of the signal. Irrespective of the initial amount of plasmids, the population of plasmids within the single cell will have reached a steady state value of <math>\frac{\beta}{\alpha}</math>. After the signal has leaved the medium in which is present, will start to decrease as a function of the degradation rate .
| + | |
- | | + | |
- | <math>\frac{dP_{1}}{dt} = - \alpha\ P_{1}
| + | |
- | </math>
| + | |
- | | + | |
- | Integrating this differential equation, <math>P_{1}</math> will be given by: <math>P_{1}(t) = e^{-\alpha t}
| + | |
- | </math> And multiplying by the steady value <math>\frac{\beta}{\alpha}</math> will yield the amount of methylated plasmids at time <math>t</math>, given that there were <math>\frac{\beta}{\alpha}</math> methylated plasmids at <math>t = 0</math>. <math>P_{1}(t) = \frac{\beta}{\alpha} e^{-\alpha t}
| + | |
- | </math>
| + | |
- | | + | |
- | == In practice ==
| + | |
- | | + | |
- | <blockquote>[[Image:bands|image]] { Gel representations for a range of different <math>F(t)</math> values. Complete methylation of all bits results in a single, bright band at the height of the length of the linearized plasmid. Decreasing the amount of methylated bits shifts the intensity of the top band away to the two bottom bands. These indicate the linearized & digested plasmid.} {fig:gelarray}
| + | |
- | </blockquote>
| + | |
- | In a typical laboratory situation, doing a restriction enzyme assay on the miniprep-extracted plasmid DNA out of followed by gel electrophoresis will be the most convenient way to assess the methylation status of the bits. The relative intensities of the gel bands can then be used to infer , which is exemplified in Figure {fig:gelarray}. Unmethylated bits will result in successfully digested DNA fragments and thus two bands of shorter DNA fragments. Methylated bits will not be cut and will therefore result in one longer band, shown more to the top of the gel. Hence the top band and the two bottom ones are mutually exclusive as they indicate the same (linearized) plasmid DNA to either be cut, resulting in the two bottom bands, or not cut, resulting in the top band.
| + | |
- | | + | |
- | The function for <math>P_{1}(t)</math> is shown in in Figure {fig:readout}. The rate with which the function approaches its limiting value of 0 is directly dependent upon the degradation rate . To get a hands-on feel of the effects on the system an interactive version of this plot, in which , and are controllable via a GUI is included in the attached Mathematica file.
| + | |
- | | + | |
- | <blockquote>[[Image:infer|image]] { Animation of <math>P(t)</math> with its value ranging between 0.01 and 1. Solving for <math>t</math>, the time after which the signal has left/been degraded in the Cell, the solution is shown by the blue dot and the text underneath the plot. The value of <math>P(t)</math> will decrease over time until it becomes 0. An interactive version of this plot has been included in the attached Mathematica file.} {fig:readout}
| + | |
- | </blockquote>
| + | |
- | = Introduction =
| + | |
- | | + | |
- | We’ve established in a very early stage of the design process that we were for growth and degradation speed of plasmid going to use DNA methylation as the molecular mechanism to create our storage mechanism. DNA methylation, usually part of the Type II restriction system, is performed by a group of proteins called the methyltransferases. In choosing DNA-methylation, we were heavily inspired by the DamID technology which has been originated by the Bas van Steensel . By fusing a bacterial methyltransferase (MTase) to a eukaryotic transcription factor which is to be characterized, van Steensel was able to infer the transcription factor binding sites on the eukaryotic genome by reading out where methylation had occurred. We thought we could reverse this idea to create memory unit, at the core constituted by an MTase that would only methylate an especially designated genomic region in reaction to the sensing of a signal. For this, we required an MTase that
| + | |
- | | + | |
- | <blockquote>
| + | |
- | is not already present in the chassis organism ''E.coli'',
| + | |
- | | + | |
- | | + | |
- | has a binding motif of that is not methylated by any other MTases in the system and for which
| + | |
- | | + | |
- | | + | |
- | a restriction enzyme (RE) dependent upon the action of this methyltransferase has been identified.
| + | |
- | </blockquote>
| + | |
- | By finding a way of controlling the activity of the MTase, a genomic region purposed to be methylated, named the ''bit region'', would only be methylated when the cell would be induced to do so by the signal. Finally, the methylation status of the bit regions can then be assessed using digestion with this RE followed by a band length analysis of a gel electrophoresis of the product. This is process is covered more in depth in section {sub:gelinfer}
| + | |
- | | + | |
- | == Choosing a methyltransferase ==
| + | |
- | | + | |
- | We wrote a python script to mine the [http://rebase.neb.com/rebase/rebase.html REBASE] database of REs and MTases for a MTase that meets our needs. We found an ideal candidate: M.ScaI, originally from ''Streptomyces caespitosus''. Now, the only remaining design question was how to control the activity of the fusion protein as we would want it to only be active in the presence of a signal. We could either try to control its ''presence'' in the cell by induced transcription or try to directly fuse it to a inducible transcriptional regulator, hoping that the blocked activity of the latter in absence of a signal would also prohibit the MTase from methylating the bit. This design question will be more elaborately discussed in section {chapter_molecular}.
| + | |
- | | + | |
- | == Extending the idea to multiple bits ==
| + | |
- | | + | |
- | Extending the design idea further, we realized that it should be possible to store the presence of different signals – multiple ''bits'' – simultaneously by fusing DNA binding proteins to the methyltransferase and adding their corresponding DNA binding motifs to the bit regions of the DNA. Either traditional bacterial transcriptional regulators could be used for this purpose or the more extendable Zinc-Finger Array (ZFA) technology . The latter allows the construction of highly specific proteins to custom DNA motifs. The large increase in binding affinity of the resulting MTase-ZF fusion protein would have to make sure that the same MTase can be used for different signals with a very low chance of cross-talk between FPs and other signals’ bits.
| + | |
- | | + | |
- | = Inferring the time of signal onset =
| + | |
- | | + | |
- | Numerous copies of identical plasmids are often present in single cells and plasmids replicate independently of the bacterial chromosome . A plasmid copy number (PCN) has been determined for all plasmids in the Parts Registry, which indicates a likely amount of copies of the plasmid to be present in each cell. Unlike eukaryotes, we discovered that DNA methylation patterns are not copied to the new strand during DNA replication in prokaryotes. This will have to lead to a dilution of the amount of plasmids with a written bit over time, mostly due to cell replication and the ensuing binomial division of the plasmids in the parent cell among the two daughter cells. This seemed like a downside to using a MTase system at first, but it quickly opened our eyes to a very exciting application of this system. Assuming nearly constant proliferation & replication rates for the cells and the within the plasmids, this dilution can be used to infer how long ago the signal of interest was registered by the cell.
| + | |
- | | + | |
- | Using the here presented model, we will examine how this phenomenon can exactly be used to infer the registration time of the signal by the cell. Another goal of this model is to gain insight into how much the various parameters for growth and degradation speed of plasmid in the model affect the ratio <math>F(t) = \frac{\text{written plasmids}(t)}{\text{total plasmids}}</math>. Unknown variables affecting <math>F(t)</math> in a real-life setting would be the time of signal onset, signal duration and signal strength. Additionally the signal response, degradation and proliferation rates of the ''Cellular Logbook'' will influence <math>F(t)</math>, but these should be experimentally determinable before actual deployment and application of our system. This way the can be used to measure a signal and estimate its time of onset. This would allow the inferrance of signal onset <math>s_{\text{on}}</math> from the fraction of methylated over total bits in the population. As this model has the sole purpose of clarification of the desired system behaviour, the units in this model are completely arbitrary.
| + | |
- | | + | |
- | First, let’s model the input signal/compound which is to be registered. Let <math>s_{\text{on}}</math> be the time at which the signal is turned on and <math>s_{\text{off}}</math> be the time at which the signal is turned off. Imagine the to be stationary and positioned along a fluidic stream with the signal to be registered passing the at <math>t = s_{\text{on}}</math> and leaving it at <math>t = s_{\text{off}}</math>. Modelling the signal using the piecewise function <math>S(t)</math> in equation {eqn:signal} seems appropriate. This function has been plotted in Figure {fig:signal}.
| + | |
- | | + | |
- | <math>S(t) = \left\{
| + | |
- | \begin{array}{lr}
| + | |
- | 0 & 0 \le t \le s_{\text{on}} \\
| + | |
- | 1 & s_{\text{on}} \le t \le s_{\text{off}} \\
| + | |
- | 0 & s_{\text{off}} \le t \le \infty
| + | |
- | \end{array}
| + | |
- | \right.
| + | |
- | \label{eqn:signal}
| + | |
- | </math>
| + | |
- | | + | |
- | <blockquote>[[Image:signal.png|image]] {fig:signal} { Plot of the input signal <math>S(t)</math> with <math>s_{\text{on}}</math> at 3 and <math>s_{{\text{off}}}</math> at 4}
| + | |
- | </blockquote>
| + | |
- | A single cell with multiple identical plasmids on which the gene for the methyltransferase and a so called ''bit region'', which is the region especially purposed to be methylated in presence of a signal, have been placed. In this cell the fraction of methylated plasmids <math>F(t)</math> will then be in direct correspondence to how long ago the signal has been registered if the following plausible rules are true:
| + | |
- | | + | |
- | * A well stirred cellular system in which no spatial concentrations differences exist and in which all species concentrations are large enough to be approximated by continuous functions. Such a system is describable by a set of ordinary differential equations (ODE’s).
| + | |
- | * To easy the analysis of the model, a single unique bit with multiple copies per cell is considered here. The system is easily upgradeable to multiple signals however.
| + | |
- | * Assumed is a high response rate ( 3 min) termed <math>\omega</math>, which is the constant rate with which the system responds with methylation to the sensing of a signal.
| + | |
- | * We will assume logistic growth for the memory cell population, with a capacity limit of <math>Ca = 5000</math>.
| + | |
- | * Alternatively, this maximal population number could also be described as the fraction between the cellular proliferation rate (<math>\beta</math>) and degradation rate (<math>\alpha</math>), as the steady state amount of cells in the population will be determined by <math>\frac{\beta}{\alpha}</math> as a solution to the differential equation <math>P'(t) = \beta P - \alpha P</math>.
| + | |
- | * Accumulation of cells due to the division of cells in which te bit has been written is assumed to result in non-written cells; methylation patterns are not copied to the progeny in prokaryotes.
| + | |
- | * <math>P_{0}</math> denote cells that have the single bit set to 0, no write event has taken place in these cells.
| + | |
- | * <math>P_{1}</math> indicate cells in which the bit has been flipped to 1 in response to encountering the signal compound.
| + | |
- | * <math>P_{0+1}</math> indicates the total amount of cells in the sytem, <math>P_{0} + P_{1}</math>
| + | |
- | | + | |
- | From these rules, the following system of ODE’s has been constructed:
| + | |
- | | + | |
- | <math>\frac{dP_{0}}{dt} = k\ P_{0+1}\ (1 - \frac{P_{0+1}}{\text{Ca}}) - \omega\ S(t)\ P_{0} - \alpha\ P_{0} \\
| + | |
- | \frac{dP_{1}}{dt} = \omega\ S(t)\ P_{0} - \alpha\ P_{1}
| + | |
- | </math>
| + | |
- | | + | |
- | <blockquote>
| + | |
- | <table>
| + | |
- | <thead>
| + | |
- | <tr class="header">
| + | |
- | <th align="left">Parameter</th>
| + | |
- | <th align="right">Value</th>
| + | |
- | </tr>
| + | |
- | </thead>
| + | |
- | <tbody>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left">k</td>
| + | |
- | <td align="right">1.5</td>
| + | |
- | </tr>
| + | |
- | <tr class="even">
| + | |
- | <td align="left"><math>\beta</math></td>
| + | |
- | <td align="right">2500</td>
| + | |
- | </tr>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left"><math>\alpha</math></td>
| + | |
- | <td align="right">.5</td>
| + | |
- | </tr>
| + | |
- | <tr class="even">
| + | |
- | <td align="left">Ca</td>
| + | |
- | <td align="right">5000</td>
| + | |
- | </tr>
| + | |
- | <tr class="odd">
| + | |
- | <td align="left"><math>\omega</math></td>
| + | |
- | <td align="right">8</td>
| + | |
- | </tr>
| + | |
- | </tbody>
| + | |
- | </table>
| + | |
| | | |
- | { Parameter values for the plasmid methylation model} {tab:populationparms}
| + | The wide range of available Zinc Fingers could also allows for the usage of the same MTase for multiple signal. This assumes that the site specificity of the Zinc Finger Array is much higher than the site specifity of the MTase, which is very plausible ([[#McNamara|3]]). |
- | </blockquote>
| + | In this setup there would be a very low chance of cross-talk between FPs and other signals’ bits. |
- | Using the parameter values of Table {tab:populationparms} a simulation with a duration of 20 time units is shown in Figure {fig:timelapse}.
| + | |
| | | |
- | <blockquote>[[Image:1|image]] { Input signal <math>S(t)</math> with <math>s_{\text{on}}</math> at 3 and <math>s_{{\text{off}}}</math> at 4} {fig:timelapse}
| + | The ZF-binding sites and methylation motif groups are distributed over the plasmid in a way that the distances between the TF-binding sites are unique: for all binding restriction sites <math>i</math> the intermediate genomic distance with the next bit, $D(i,i+1)$, has to be: |
- | </blockquote>
| + | $$ |
- | <blockquote>[[Image:2|image]] { Sum of unmethylated and methylated cells, <math>S_{0+1}</math>} {fig:totalcells}
| + | D(i,i+1) = \frac{D(i-1,i)}{2} |
- | </blockquote>
| + | $$ |
- | The plasmid population within a is shown to be completely converted to methylated plasmids shortly after . As long as the signal is still present – until <math>s_{\text{off}}</math>, – the bit on all newly copied plasmids will be immediately methylated as the signal is still present. After <math>s_{\text{off}}</math>, the percentage of <math>\frac{\text{methylated} \text{cells}}{\text{total} \text{cells}}</math> will be diluted, due to cell proliferation implicating plasmid division between daughter cells and cell death. The size of the time span with which a hint of methylated cells is still present therefore relates to two factors: positively with the amount of methylated cells at <math>s_{\text{off}}</math> and negatively with the cellular death rate.
| + | By designing the plasmid to adhere to this rule, a unique pattern of fragment lengths between each of the TF binding sites - or bits - will result for each combination of bit statuses. |
| + | One plasmid will then be able to hold |
| + | <math> |
| + | \log_2 (\text{Plasmid BP length}) - \log_2 (256) |
| + | </math> bits |
| + | with the current achievable resolution in gel electrophoresis techniques (resolution down to 300 bp and <math>2^{7} = 256</math>). |
| + | For the typical plasmid base pair length of <math>6 \cdot 10^{4}</math> nucleotides this is <math>8\ \text{bits}</math>. |
| | | |
- | The bit writing speed termed <math>\omega</math>, or the rate with which the presence of a signal causes the methylation of a bit, is likely to be very high in comparison to the cellular division rate. The underlying required processes such as the diffusion of the signal across the cellular membrane, the signal inducing gene transcription of the FP, mRNA translation and transcription & and FP diffusion over to the bit all occur well within 5-minute time scales. Every single gene on a plasmid will then be expected to be methylated within at most 5 minutes of registering of the signal.
| + | [[File:Multiplebits.png|center|thumb|400px|A plasmid storing three bits in which site 0 is always cleaved. Notice the unique distances between bits 1, 2 and 3. The value of each bit can be inferred from the unique band lengths pattern after digestion]] |
| | | |
- | The assumption that the population will have reached its steady state level (SS), given by the production rate (<math>\beta</math>) divided by the degradation rate (<math>\alpha</math>), before is a plausible one. The ''Cellular Logbook'' population will have been stored inside of a vesicle consisting of a non-bacterial permeating membrane with a very small volume. Due to this volume and the rapid bacterial division time (<math>\sim</math> 40 mins), the spatial growth limits of the system will then have been reached in a short amount of time. The steady state-level of the bacterial population has already been reached before the signal is registered by the population. | + | == Extending from multiple to ~64 bits == |
| + | The amount of bits stored on a single plasmid could be enlarged even more. In the previous section all bits used the same MTase and are digested using the same restriction enzyme with the potential to store up to 8 bits. By finding a new MTase-RE combination and engineering the memory plasmid to have bits for this new combination slightly offset from the original MTase-RE binding sites, 8 additional bits will be storable on the plasmid. This requires that during the read-out step, digestions are performed for each RE individually. It is essential that only a single RE-type is for each individual plasmid extract to retain the ability to read out the bits from the unique band lengths. For each RE used, a separate band on the final gel must be used and determination of bit values will proceed in the same way as for single MTase memory modules. A plasmid bp-size of 60k will allow for about 8 different RE’s and thus: |
| | | |
- | To reinforce that the following equation always holds in this model; <math>P_{0} + P_{1} = P_{0+1}</math> the following plot is shown, in which the sum of the methylated and unmethylated cells over time is plotted. This clearly shows the limit, specified by the Capacity limit, which is reached around <math>t=6</math>.
| + | $$ |
| + | 8\ [\text{RE}] \cdot [\frac{8\ \text{bits}}{\text{RE}}] = 64\ \text{bits} |
| + | $$ |
| | | |
- | == Methylated plasmids over time == | + | == Signal strength measurement == |
| + | Binary encoding of a signal's presence is useful, but often in biology gaining insight in the quantity of a signal is even more relevant. By coupling different promoters with different strengths as responses to the same signal and assigning each promoter its own unique bit, the range in which the strength of the signal is located can be measured. |
| | | |
- | {sub:gelinfer} The monotonically decreasing value of <math>F(t) = \frac{\text{methylated plasmids}}{\text{total plasmids}}</math> can be used to infer , given that the cellular proliferation rate (<math>\beta</math>) and degradation rate (<math>\alpha</math>) are known and constant and the assuming that all bits are methylated during the presence of the signal. Irrespective of the initial amount of plasmids, the population of plasmids within the single cell will have reached a steady state value of <math>\frac{\beta}{\alpha}</math>. After the signal has leaved the medium in which is present, will start to decrease as a function of the degradation rate .
| + | == References == |
| + | <span id="Fu"> |
| + | <sup>Fu, Fengli, Jeffry D. Sander, Morgan Maeder, Stacey Thibodeau-Beganny, J. Keith Joung, Drena Dobbs, Leslie Miller, and Daniel F. Voytas. 2009. “Zinc Finger Database (ZiFDB): a repository for information on C2H2 zinc fingers and engineered zinc-finger arrays.” ''Nucleic acids research'' 37 (jan): 279. doi:10.1093/nar/gkn606. [[http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2686427\&tool=pmcentrez\&rendertype=abstract|http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2686427\&tool=pmcentrez\&rendertype=abstract]]</sup> |
| + | </span> |
| | | |
- | <math>\frac{dP_{1}}{dt} = - \alpha\ P_{1} | + | <span id="Kaseniit"> |
- | </math> | + | <sup> |
| + | Kaseniit, K. E., Perli, S. D., & Lu, T. K. (2011). Designing extensible protein-DNA interactions for synthetic biology. 2011 IEEE Biomedical Circuits and Systems Conference (BioCAS), 349–352. doi:10.1109/BioCAS.2011.6107799 |
| + | </sup> |
| + | </span> |
| | | |
- | Integrating this differential equation, <math>P_{1}</math> will be given by: <math>P_{1}(t) = e^{-\alpha t}
| + | <span id="McNamara"> |
- | </math> And multiplying by the steady value <math>\frac{\beta}{\alpha}</math> will yield the amount of methylated plasmids at time <math>t</math>, given that there were <math>\frac{\beta}{\alpha}</math> methylated plasmids at <math>t = 0</math>. <math>P_{1}(t) = \frac{\beta}{\alpha} e^{-\alpha t}
| + | <sup>McNamara, A. R., Hurd, P. J., Smith, A. E. F., & Ford, K. G. (2002). Characterisation of site-biased DNA methyltransferases: specificity, affinity and subsite relationships. Nucleic acids research, 30(17), 3818–30. Retrieved from http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=137423&tool=pmcentrez&rendertype=abstract |
- | </math> | + | </sup> |
| + | </span> |
| | | |
- | == In practice == | + | <span id="Steensel1"> |
| + | <sup> |
| + | van Steensel, B., Delrow, J., & Henikoff, S. (2001). Chromatin profiling using targeted DNA adenine methyltransferase. Nature genetics, 27(3), 304–8. doi:10.1038/85871 |
| + | </span> |
| + | </sup> |
| | | |
- | <blockquote>[[Image:bands|image]] { Gel representations for a range of different <math>F(t)</math> values. Complete methylation of all bits results in a single, bright band at the height of the length of the linearized plasmid. Decreasing the amount of methylated bits shifts the intensity of the top band away to the two bottom bands. These indicate the linearized & digested plasmid.} {fig:gelarray} | + | <span id="Steensel2"> |
- | </blockquote>
| + | <sup> |
- | In a typical laboratory situation, doing a restriction enzyme assay on the miniprep-extracted plasmid DNA out of followed by gel electrophoresis will be the most convenient way to assess the methylation status of the bits. The relative intensities of the gel bands can then be used to infer , which is exemplified in Figure {fig:gelarray}. Unmethylated bits will result in successfully digested DNA fragments and thus two bands of shorter DNA fragments. Methylated bits will not be cut and will therefore result in one longer band, shown more to the top of the gel. Hence the top band and the two bottom ones are mutually exclusive as they indicate the same (linearized) plasmid DNA to either be cut, resulting in the two bottom bands, or not cut, resulting in the top band.
| + | van Steensel, B., & Henikoff, S. (2000). Identification of in vivo DNA targets of chromatin proteins using tethered dam methyltransferase. Nature biotechnology, 18(4), 424–8. doi:10.1038/74487 |
| + | </sup> |
| + | </span> |
| | | |
- | The function for <math>P_{1}(t)</math> is shown in in Figure {fig:readout}. The rate with which the function approaches its limiting value of 0 is directly dependent upon the degradation rate . To get a hands-on feel of the effects on the system an interactive version of this plot, in which , and are controllable via a GUI is included in the attached Mathematica file.
| + | <html></div></html> |
| + | <div id="sub-menu" class="content-block" align="center"> |
| + | [[File:Amsterdam_previous.png|left|link=Team:Amsterdam/project/background/]] |
| + | [[File:Amsterdam_next.png|right|link=Team:Amsterdam/modeling/odemodel]] |
| + | '''Click the arrows to go to the previous or next page''' |
| + | </div> |
| + | </div> |
| | | |
- | <blockquote>[[Image:infer|image]] { Animation of <math>P(t)</math> with its value ranging between 0.01 and 1. Solving for <math>t</math>, the time after which the signal has left/been degraded in the Cell, the solution is shown by the blue dot and the text underneath the plot. The value of <math>P(t)</math> will decrease over time until it becomes 0. An interactive version of this plot has been included in the attached Mathematica file.} {fig:readout}
| + | {{Team:Amsterdam/Foot}} |
- | </blockquote> | + | </html> |