Difference between revisions of "Team:William and Mary/Description"

 
(82 intermediate revisions by 4 users not shown)
Line 4: Line 4:
 
{{Team:William_and_Mary/HEADER2}}
 
{{Team:William_and_Mary/HEADER2}}
  
<!--website by John Mitchell, iGEM 2016
 
Bootstrap template by creative-tim-->
 
  
 
<html lang="en">
 
<html lang="en">
 
<cfheader name="X-XSS-Protection" value="1">
 
<cfheader name="X-XSS-Protection" value="1">
  
<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/4.7.0/css/font-awesome.min.css">
+
 
 
          
 
          
         <div class="wrapper">
+
          
+
 
+
<style>
<div>
+
      body {
 +
        font-family: 'PT Sans', sans-serif;
 +
        font-size: 13px;
 +
      }
 +
 
 +
</style>
 +
 
 +
</head>
 +
 
 +
 
 +
<body>
 +
<br></br>
 +
<div style='padding-top: 30px;'></div>
 +
 
 +
<center><img src="https://static.igem.org/mediawiki/2017/3/36/T--William_and_Mary--descript.jpeg" width=380px;/></center>
 +
 
 +
 
 +
<div style='padding-top: 30px;'></div>
 +
 
 +
<div style = 'padding-left: 190px; padding-bottom: 10px;font-size: 25px' ><b>Motivation</b></div>
 +
 
 +
<div style='background: #808080; margin: 0px 190px 20px 190px; height:1px;></div>
 +
 
 +
<div style='padding-top: 0px;'></div>
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 50px;line-height: 25px;' >
 +
One of the main goals of synthetic biology is to create a modular genetic basis for the independent control of circuit behavior properties. Much progress has been made in achieving this aim for properties like gene expression strength (where well-characterized ribosome binding sites (RBSs) can be swapped within a genetic part), circuit architecture (where promoters can be swapped out to introduce connections and feedback architectures), and even gene expression noise (through a combination of the above two modulations). However, in order to move into the next phase of synthetic biology, we need to be able to control the dynamical properties of circuits— we want to move beyond circuits that focus on endpoint, steady-state values and explore the rich variety of dynamical systems. Fundamentally, gaining controlling of dynamical systems implies gaining control of temporal dynamics. 
 +
</div>
 +
 
 +
 
 +
<div style='padding-top: 15px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 50px;line-height: 25px;' >
 +
Currently, there is no good way to control the temporal dynamics of gene expression. Current control strategies require either a rewiring of the circuit architecture to achieve different time-dependent dynamics [1, 2] or a complete circumvention of transcriptional circuitry altogether, relying on post-translational dynamics like phosphorylation [3] or protein-protein interactions [4] to transmit information through a circuit. These approaches are often inaccessible to iGEM teams because they require too drastic an overhaul of existing circuit implementations. To alleviate this issue, and to enable future iGEM teams to create robust dynamical circuits, we created a protein degradation based ‘plug-and-play’ style system that allows modular and predictable control of the gene expression speed of a given circuit without requiring a fundamental redesign of existing circuit architecture.
 +
</div>
 +
 
 +
 
 +
<div style='padding-top: 60px;'></div>
 +
 
 +
 
 +
<div style = 'padding-left: 190px; padding-bottom: 10px;font-size: 25px' ><b>Degradation Based Control of Gene Expression Speed</b></div>
 +
 
 +
 
 +
<div style='background: #808080; margin: 0px 190px 20px 190px; height:1px;></div>
 +
 
 +
<div style='padding-top: 0px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 50px;line-height: 25px;' >The relationship between the degradation rate of a gene's protein product and the response speed of its expression is a well-established concept in the theoretical biology community. Consider the following simple mathematical model presented by Uri Alon [5]. X is an inducible gene which, once activated, produces its protein product x at rate α. If x is degraded at rate γ, then we can use the following differential equation to track the concentration of protein x over time:</div>
 +
 
 +
<figure style='padding-left: 200px;'>
 +
<img src='https://static.igem.org/mediawiki/2017/0/02/T--William_and_Mary--truDescriptionMath1.png' height = "70%" width = "70%"/>
 +
<figcaption><div style='padding-left: 5px;padding-top: 15px; color: #808080; font-size: 14px;'>Figure 1: A simple kinetic model tracking the changing concentration of protein x over time. </div><div style='padding-left: 5px;color: #808080; font-size: 14px;'>Gene circuit component images adapted from Newcastle iGEM 2010; protein images adapted from Cameron and Collins 2014 [7]. </div></figcaption>
 +
</figure>
 +
 
 +
<div style='padding-top: 15px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; line-height: 25px;' >We solve the above differential equation to obtain the following function, which gives the protein concentration at a given time t after activation of X: </div>
 +
 
 +
<figure style='padding-left: 200px;'>
 +
<img src='https://static.igem.org/mediawiki/2017/7/71/T--William_and_Mary--truDescriptionMath2.png' height = "30%" width = "30%"/>
 +
<figcaption><div style='padding-left: 5px;padding-top: 15px; color: #808080; font-size: 14px;'>Figure 2: Concentration of protein x as a function of time; the solution to the differential equation in Figure 1. </div></figcaption>
 +
</figure>
 +
 
 +
<div style='padding-top: 15px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; line-height: 25px;' >We are able to graph this function and determine that the steady state concentration for x is equal to production rate over degradation rate. Finally, defining т<sub>1/2</sub> as the time it takes for x to reach half its steady-state concentration, we find the following: </div>
 +
 
 +
<figure style='padding-left: 200px;'>
 +
<img src='https://static.igem.org/mediawiki/2017/6/65/T--William_and_Mary--DescriptionMath3.png' height = "70%" width = "70%"/>
 +
<figcaption><div style='padding-left: 5px;padding-top: 15px; color: #808080; font-size: 14px;'>Figure 3: We find that т<sub>1/2</sub> is a function of γ, the degradation rate for protein x. </div></figcaption>
 +
</figure>
 +
 
 +
<div style='padding-top: 15px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; line-height: 25px;' >The above expression reveals that in this model, it is not only sufficient but also necessary to tune the degradation rate γ in order to control the gene expression speed (represented here as т<sub>1/2</sub>). This property remains true for more complex circuit architectures like feedforward loops, provided that the steady-state concentration is held constant [5]. </div>
 +
 
 +
<div style='padding-top: 15px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 50px;line-height: 25px;' >Thus in order to create a plug-and-play style system to control gene expression speed, we developed and characterized a suite of BioBrick parts which allow for simple, modular and predictable changes to the gene expression speed of arbitrary proteins via protein degradation.</div>
 +
 
 +
 
 +
<div style='padding-top: 60px;'></div>
 +
 
 +
 
 +
<div style = 'padding-left: 190px; padding-bottom: 10px;font-size: 25px' ><b>Orthogonal Degradation Tags</b></div>
 +
 
 +
 
 +
<div style='background: #808080; margin: 0px 190px 20px 190px; height:1px;></div>
 +
 
 +
<div style='padding-top: 0px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 50px;line-height: 25px;' >When we created our gene expression control system, we wanted to make sure that it was both usable across a variety of biological systems and circuits, and easily accessible to other iGEM teams. To this end, we chose to use the <i>Mesoplasma florum</i> Lon (mf-Lon) protease system which was characterized by Gur and Sauer in 2008 [6] and developed into a modular suite of genetic parts by Cameron and Collins in 2014 [7]. This system consists of of a AAA+ protease and its associated protein degradation tags (pdt), which operate in a mechanistically similar manner to the E. Coli endogenous protease ClpXP and its associated ssrA tags. However, unlike ClpXP and ssrA tags, mf-Lon and pdts are completely orthogonal to the endogenous protein degradation systems in E. coli. Using this orthogonal degradation system helps eliminate cross-talk between our system and endogenous E. coli proteins. Further, since there are pdts with a wide range of different affinities, we are able to tune degradation rate, and thus gene expression speed to a wide variety of values. This represents not only a practical advancement in the tuning of circuits, but also an important contribution in our theoretical understanding of biological processes, as to our knowledge there has so far been no direct experimental measurement in the literature of the predicted 1/γ scaling of a gene's speed with its degradation rate. See our <a href='https://2017.igem.org/Team:William_and_Mary/Results' style='text-decoration: underline;'>results</a> page for more details.</div>
 +
 
 +
<div style='padding-top: 15px;'></div>
 +
 
 +
<div style='padding-top: 60px;'></div>
 +
 
 +
 
 +
<div style = 'padding-left: 190px; padding-bottom: 10px;font-size: 25px' ><b>Modeling</b></div>
 +
 
 +
 
 +
<div style='background: #808080; margin: 0px 190px 20px 190px; height:1px;></div>
 +
 
 +
<div style='padding-top: 0px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 50px;line-height: 25px;' >Finally, by constructing a mechanistic model of our degradation system, we were able to rigorously analyze our timeseries datasets using Bayesian Parameter Inference and obtain parameter distributions (Figure 4). We then used this analysis to feed back into our predictive saturation model, which we used to develop the beginnings of a theoretical framework for understanding the relationship between protease-driven speed control and saturation effects on the protease. By combining our characterization measurements with our conceptual insights, future iGEM teams are now enabled to use the pdt system to predictably control the temporal dynamics of their genetic circuits. See our <a href='https://2017.igem.org/Team:William_and_Mary/Model' style='text-decoration: underline;'>modeling</a> page for details
 +
</div>
 +
 
 
<center>
 
<center>
<img src="https://static.igem.org/mediawiki/2017/f/fe/T--William_and_Mary--Description2.jpeg" height = "12%" width = "100%"/>
+
<figure style='padding-left: px;'>
 +
<img src='https://static.igem.org/mediawiki/2017/5/5a/T--William_and_Mary--MCMC-F2.png' width = "45%"/>
 +
<figcaption><div style='padding-left: 20%;padding-right:20%; padding-top: 15px; color: #808080; font-size: 14px;'>Figure 4: MCMC parameter estimation for a simple protein production model correctly estimates the simulated values for both beta and gamma (15 and .03 respectively). In addition, the MCMC identifies a strong positive correlation between beta and gamma. This makes sense as steady-state levels of protein are set by the ratio Beta/Gamma. Therefore as beta goes up we should see a corresponding increase in gamma to best fit the model.
 +
</div></figcaption>
 +
</figure>
 
</center>
 
</center>
 +
 +
<div style='padding-top: 15px;'></div>
 +
 +
<div style='padding-top: 60px;'></div>
 +
 +
<!---References----->
 +
 +
 +
 +
<div style = 'padding-left: 190px; padding-bottom: 10px;font-size: 25px' ><b>References</b></div>
 +
 +
<div style='background: #808080; margin: 0px 190px 20px 190px; height:1px;></div>
 +
 +
<div style='padding-top: 0px;'></div>
 +
 +
<div style='padding-top: 30px;'></div>
 +
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 0px;line-height: 25px;' >[1]  Nitzan Rosenfeld, Michael B Elowitz, and Uri Alon. Negative autoregulation speeds the response times of transcription networks. Journal of molecular biology, 323(5):785–793, 2002.
 
</div>
 
</div>
  
 +
<div style='padding-top: 15px;'></div>
  
  
<div class="section section-we-are-2">
 
<div class="text-area">
 
<div class="container">
 
<div class="row">
 
<p>
 
A fundamental goal of synthetic biology is to create a modular genetic basis for control over circuit behavior properties. While much progress has been made in achieving this objective for properties such as gene expression strength, where well-characterized ribosome binding sites (RBSs) can be conveniently swapped within a genetic part, there is much to be desired in altering gene expression speed. Currently, no such robust mechanism of speed control exists. We intend to provide a means to tune the speed of gene expression in transcriptional circuits, where a pre-characterized genetic part can be inserted into a gene to alter its expression in a predictable way.
 
<p>
 
Based on multiple well-cited claims in literature, a strong relationship can be asserted between gene expression speed and the rate of protein degradation. Using a basic mathematical model of gene expression, one can derive that the speed of gene expression, defined as the time it takes for its protein product to reach half of its steady-state concentration, is a function of the protein’s degradation rate. This reveals that tuning protein degradation rate is essential to controlling gene expression, thus amenable to an approach involving protein degradation tags. Degradation tags are used endogenously to identify misfolded proteins, and different tags have unique protease-binding affinities which confer various degradation rates on the tagged proteins.
 
</p>
 
<p>
 
In 2008, the Sauer lab at MIT reported that Mycoplasma florum’s Lon protease system was orthogonal to the endogenous protein degradation machinery in E. coli. As of now, mf-Lon degradation tags exist only as isolated sequences on the BioBrick Registry. We intend to build a suite of BioBrick parts in the form of [mf-Lon tag] - [Stop Codon] - [Double Terminator] that can be swapped in to directly modulate protein degradation rate just as RBS and promoter sequences can be swapped to modify protein production. This would drastically increase the accessibility of mf-Lon tags and enable other teams to easily amplify their desired protein sequence by simply cloning it into our construct. Should the relationship between gene expression speed and protein degradation rate exist robustly, we will be providing the first modular genetic basis of speed control, fulfilling a core aspiration of synthetic biology to have every gene expression and circuit property accessible at the genetic level.</p>
 
</p>
 
</div>
 
</div>
 
 
 
 
  
+
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 0px;line-height: 25px;' >[2] Shmoolik Mangan and Uri Alon. Structure and function of the feed-forward loop network motif. Proceedings of the National Academy of Sciences, 100(21):11980–11985, 2003.
</div>
+
</div>
</div>
+
<div style='padding-top: 15px;'></div>
<!--
+
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 0px;line-height: 25px;' >
+
[3] Russell M Gordley, Reid E Williams, Caleb J Bashor, Jared E Toettcher, Shude Yan, and Wendell A Lim. Engineering dynamical control of cell fate switching using synthetic phospho- regulons. Proceedings of the National Academy of Sciences, 113(47):13528–13533, 2016.</div>
</div>
+
 
</div>
+
<div style='padding-top: 15px;'></div>
-->
+
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 0px;line-height: 25px;' >
</div>
+
[4] Anselm Levskaya, Orion D Weiner, Wendell A Lim, and Christopher A Voigt. Spatiotemporal control of cell signalling using a light-switchable protein interaction. Nature, 461(7266):997– 1001, 2009.
 +
</div>
 +
<div style='padding-top: 15px;'></div>
 +
 
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 0px;line-height: 25px;' >
 +
[5]  Uri Alon. An introduction to systems biology: design principles of biological circuits. CRC press, 2006
 +
</div>
 +
 
 +
<div style='padding-top: 15px;'></div>
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 0px;line-height: 25px;' >
 +
[6] Eyal Gur and Robert T Sauer. Evolution of the ssra degradation tag in mycoplasma: specificity switch to a different protease. Proceedings of the National Academy of Sciences, 105(42):16113– 16118, 2008.
 +
</div>
 +
 
 +
<div style='padding-top: 15px;'></div>
 +
<div style = 'padding-right: 190px; padding-left: 190px; text-indent: 0px;line-height: 25px;' >
 +
[7] D Ewen Cameron and James J Collins. Tunable protein degradation in bacteria. Nature biotechnology, 32(12):1276–1281, 20
 +
</div>
 +
 
 +
 
 +
 
 +
 
 +
 
 +
</body>
 +
 
 +
 
 
</html>
 
</html>
 +
 +
 +
  
 
{{Team:William_and_Mary/JS_SOURCE2}}
 
{{Team:William_and_Mary/JS_SOURCE2}}
 
{{Team:William_and_Mary/FOOTER_SMALL}}
 
{{Team:William_and_Mary/FOOTER_SMALL}}

Latest revision as of 03:07, 2 November 2017




Motivation
One of the main goals of synthetic biology is to create a modular genetic basis for the independent control of circuit behavior properties. Much progress has been made in achieving this aim for properties like gene expression strength (where well-characterized ribosome binding sites (RBSs) can be swapped within a genetic part), circuit architecture (where promoters can be swapped out to introduce connections and feedback architectures), and even gene expression noise (through a combination of the above two modulations). However, in order to move into the next phase of synthetic biology, we need to be able to control the dynamical properties of circuits— we want to move beyond circuits that focus on endpoint, steady-state values and explore the rich variety of dynamical systems. Fundamentally, gaining controlling of dynamical systems implies gaining control of temporal dynamics.
Currently, there is no good way to control the temporal dynamics of gene expression. Current control strategies require either a rewiring of the circuit architecture to achieve different time-dependent dynamics [1, 2] or a complete circumvention of transcriptional circuitry altogether, relying on post-translational dynamics like phosphorylation [3] or protein-protein interactions [4] to transmit information through a circuit. These approaches are often inaccessible to iGEM teams because they require too drastic an overhaul of existing circuit implementations. To alleviate this issue, and to enable future iGEM teams to create robust dynamical circuits, we created a protein degradation based ‘plug-and-play’ style system that allows modular and predictable control of the gene expression speed of a given circuit without requiring a fundamental redesign of existing circuit architecture.
Degradation Based Control of Gene Expression Speed
The relationship between the degradation rate of a gene's protein product and the response speed of its expression is a well-established concept in the theoretical biology community. Consider the following simple mathematical model presented by Uri Alon [5]. X is an inducible gene which, once activated, produces its protein product x at rate α. If x is degraded at rate γ, then we can use the following differential equation to track the concentration of protein x over time:
Figure 1: A simple kinetic model tracking the changing concentration of protein x over time.
Gene circuit component images adapted from Newcastle iGEM 2010; protein images adapted from Cameron and Collins 2014 [7].
We solve the above differential equation to obtain the following function, which gives the protein concentration at a given time t after activation of X:
Figure 2: Concentration of protein x as a function of time; the solution to the differential equation in Figure 1.
We are able to graph this function and determine that the steady state concentration for x is equal to production rate over degradation rate. Finally, defining т1/2 as the time it takes for x to reach half its steady-state concentration, we find the following:
Figure 3: We find that т1/2 is a function of γ, the degradation rate for protein x.
The above expression reveals that in this model, it is not only sufficient but also necessary to tune the degradation rate γ in order to control the gene expression speed (represented here as т1/2). This property remains true for more complex circuit architectures like feedforward loops, provided that the steady-state concentration is held constant [5].
Thus in order to create a plug-and-play style system to control gene expression speed, we developed and characterized a suite of BioBrick parts which allow for simple, modular and predictable changes to the gene expression speed of arbitrary proteins via protein degradation.
Orthogonal Degradation Tags
When we created our gene expression control system, we wanted to make sure that it was both usable across a variety of biological systems and circuits, and easily accessible to other iGEM teams. To this end, we chose to use the Mesoplasma florum Lon (mf-Lon) protease system which was characterized by Gur and Sauer in 2008 [6] and developed into a modular suite of genetic parts by Cameron and Collins in 2014 [7]. This system consists of of a AAA+ protease and its associated protein degradation tags (pdt), which operate in a mechanistically similar manner to the E. Coli endogenous protease ClpXP and its associated ssrA tags. However, unlike ClpXP and ssrA tags, mf-Lon and pdts are completely orthogonal to the endogenous protein degradation systems in E. coli. Using this orthogonal degradation system helps eliminate cross-talk between our system and endogenous E. coli proteins. Further, since there are pdts with a wide range of different affinities, we are able to tune degradation rate, and thus gene expression speed to a wide variety of values. This represents not only a practical advancement in the tuning of circuits, but also an important contribution in our theoretical understanding of biological processes, as to our knowledge there has so far been no direct experimental measurement in the literature of the predicted 1/γ scaling of a gene's speed with its degradation rate. See our results page for more details.
Modeling
Finally, by constructing a mechanistic model of our degradation system, we were able to rigorously analyze our timeseries datasets using Bayesian Parameter Inference and obtain parameter distributions (Figure 4). We then used this analysis to feed back into our predictive saturation model, which we used to develop the beginnings of a theoretical framework for understanding the relationship between protease-driven speed control and saturation effects on the protease. By combining our characterization measurements with our conceptual insights, future iGEM teams are now enabled to use the pdt system to predictably control the temporal dynamics of their genetic circuits. See our modeling page for details
Figure 4: MCMC parameter estimation for a simple protein production model correctly estimates the simulated values for both beta and gamma (15 and .03 respectively). In addition, the MCMC identifies a strong positive correlation between beta and gamma. This makes sense as steady-state levels of protein are set by the ratio Beta/Gamma. Therefore as beta goes up we should see a corresponding increase in gamma to best fit the model.
References
[1]  Nitzan Rosenfeld, Michael B Elowitz, and Uri Alon. Negative autoregulation speeds the response times of transcription networks. Journal of molecular biology, 323(5):785–793, 2002.
[2] Shmoolik Mangan and Uri Alon. Structure and function of the feed-forward loop network motif. Proceedings of the National Academy of Sciences, 100(21):11980–11985, 2003.
[3] Russell M Gordley, Reid E Williams, Caleb J Bashor, Jared E Toettcher, Shude Yan, and Wendell A Lim. Engineering dynamical control of cell fate switching using synthetic phospho- regulons. Proceedings of the National Academy of Sciences, 113(47):13528–13533, 2016.
[4] Anselm Levskaya, Orion D Weiner, Wendell A Lim, and Christopher A Voigt. Spatiotemporal control of cell signalling using a light-switchable protein interaction. Nature, 461(7266):997– 1001, 2009.
[5]  Uri Alon. An introduction to systems biology: design principles of biological circuits. CRC press, 2006
[6] Eyal Gur and Robert T Sauer. Evolution of the ssra degradation tag in mycoplasma: specificity switch to a different protease. Proceedings of the National Academy of Sciences, 105(42):16113– 16118, 2008.
[7] D Ewen Cameron and James J Collins. Tunable protein degradation in bacteria. Nature biotechnology, 32(12):1276–1281, 20