Imperial College/Courses/2010/Synthetic Biology/Computer Modelling Practicals/Design: Difference between revisions

From OpenWetWare
Jump to navigationJump to search
No edit summary
No edit summary
Line 11: Line 11:
</div>
</div>


This page includes extra material for the course of synthetic biology.  
This page includes extra material for the course of synthetic biology. '''The material presented in this session is not part of your coursework.
Design of synthetic biological pathways is in general a very complicated affair.  
It is however, useful (very useful) for the rest of the course, especially the mini-iGEM project.
'''
 
Computer-Assisted Design of synthetic biological pathways is in general a very complicated affair. You must, by now, be aware of some of the reasons for this.
* the behaviour depends on the parameters of the system
** there may be many
** we may not know them with enough accuracy - sometimes not at all
** a small change in a parameter may lead to a totally different behaviour (bifurcation)
* initial conditions are also liable to have an influence (the arguments regarding the parameters mostly apply to the initial conditions too)
As you must have seen with the case of the repressilator, 3 genes are enough to generate a pathway with 'interesting' properties.
 
 
The situation is unfortunately worse. Even if there is a combination of parameters that
* some basic properties of the cell  have a significant impact on the effective dynamics of pathways. Take for instance the growth rate:
** it appears in the dilution term of proteins (easy to incorporate into the model)
** but is also affects in a highly nonlinear way the gene copy number
** it affects the concentration of free and bound RNAp and therefore the level of transcription etc..
* some modules in your system may be very hard to model (if at all possible)
** for instance transport of molecules through a membrane and diffusion phenomena can be modelled but it becomes complicated fast
** in a model, errors pile up so much so that after a while the predictive power of your model is negligible.
* your synthetic pathway may 'cross-talk' with natural pathways; since we are not able to model the whole metabolism of the cell this crosstalk effect can not be assessed.




As you must have seen with the case of the repressilator, 3 genes are enough to generate a pathway with 'interesting' properties.
Fortunately, designing simple pathways with predictable properties/functions is possible, even without the extensive use of software.
Fortunately, designing simple pathways with predictable properties/functions is possible, even without the extensive use of software.
This session aims at introducing to you the basic tools and techniques of design.
This session aims at introducing to you the basic tools and techniques of design (without which no computer-assisted design is possible.
'''Please Note:'''
'''Remember: in practice it gets very complicated, very fast...'''
* The material presented in this session is not part of your coursework
 
* It is however, useful (very useful) for the rest of the course, especially the mini-iGEM project




Line 62: Line 80:
|}
|}


'''The following questions constiture the last part of your coursework (Section E):'''
* '''Question 1 :''' Before studying the general properties of the mystery circuit, let us study a simplified version of it. The ODE system indeed contains a lot of symmetries that can be exploited and yield surprising properties. Let us consider a particular choice of initial conditions: the initial conditions of the mRNA terms are all equal and the initial conditions of the protein terms are also all equal.
** '''Q1.1 For Bioengineers only''': It can be shown that for such a choice of initial conditions leads, the mRNA concentrations all remain equal (but vary with time) and so do the protein concentrations.Can you explain briefly why this is the case? '''Note:''' you do not have to submit a full mathematical proof.
** Q1.2: We now call X the mRNA concentration and Y the protein concentration. What system of ODE does X and Y satisfy? Show that this system corresponds to an auto-regulated sytem.
** Q1.3: What is the natural choice for the initial conditions of the system (justify)?
** Q1.4: Simulate the new simplified system for a=10, b=1000 and n=2 for these initial conditions and comment.
** Q1.5: Now simulate the new simplified system for a=200, b=5 and n=2 for these intial conditions and comment.
** Q1.6: For either of these choice of parameters, let the initial conditions of the simplified system vary and comment on your results.
** Q1.7: The simplified ODE system and corresponding auto-regulated behaviour is intersing from a theoretical point of view but has no experimental/practical relevance. Can you explain why?


* '''Question 2 :''' Now we return to the case a=10, b=1000 and n=2. The simplified system studied in the previous question is of course not representative of the overall behaviour of the mystery circuit.
** Q2.1: Let us assume that we can purify one of the proteins so that its initial condition is 1 and the other initial conditions are 0. Run the simulation. What happens?
** Q2.2: Describe the properties of the system (the simulation you have run is representative)
** Q2.3: Can you give a qualitative explanation for the behaviour of the system?
* '''Question 3 :''' Now we return to the case a=200, b=5 and n=2 and we seek to investigate its behaviour
** Q3.1: Again we assume that we can purify one of the proteins so that its initial condition is 1. Run the simulation again. What happens?
** Q3.2: Describe the properties of the system
** Q3.3: Can you give a qualitative explanation for the behaviour of the system?
* '''Question 4:''' The interesting property takes some time to emerge. From an experimental point of view, this is a problem ( the synthetic plasmid is liable to be rejected, the growth medium may run out of nutrients etc...). In iGEM 2007 the Imperial College team investigated playing on the initial conditions of their simple synthetic system so that their system had better proporties. We can do the same here.
** Let us assume that we can purify all the proteins so that you can set up their initial conditions. You want to determine the initial conditions that will make the property emerge fastest. How would you do this?
** '''Please note:''' You are not asked to run all the relevant simulations, just explain how you can solve such a problem.
** The question has at least two distinct components which must be adressed:
*** A practical strategy to browse through the space of possible parameters
*** A practical set of criteria to determine whether the property has emerged
** A few simulations or drawings are always welcome as they usually make explanations simpler...
* '''Question 5:''' The circuit is called the repressilator.Can you explain why?
* '''Hint''': This circuit appeared was in an article published in 2000 by M.Elowitz. Track it if you need help!




Line 102: Line 93:
{| border="1"  cellpadding="10" cellspacing="2" style="width:700px; border:3px green"
{| border="1"  cellpadding="10" cellspacing="2" style="width:700px; border:3px green"
!Model
!Model
!CellDesigner Instructions
!The basic delay
|-
|-
|
|

Revision as of 09:05, 2 February 2010

Complementary Session: Introduction to the Design of Biological circuits


Objectives:

This page includes extra material for the course of synthetic biology. The material presented in this session is not part of your coursework. It is however, useful (very useful) for the rest of the course, especially the mini-iGEM project.

Computer-Assisted Design of synthetic biological pathways is in general a very complicated affair. You must, by now, be aware of some of the reasons for this.

  • the behaviour depends on the parameters of the system
    • there may be many
    • we may not know them with enough accuracy - sometimes not at all
    • a small change in a parameter may lead to a totally different behaviour (bifurcation)
  • initial conditions are also liable to have an influence (the arguments regarding the parameters mostly apply to the initial conditions too)

As you must have seen with the case of the repressilator, 3 genes are enough to generate a pathway with 'interesting' properties.


The situation is unfortunately worse. Even if there is a combination of parameters that

  • some basic properties of the cell have a significant impact on the effective dynamics of pathways. Take for instance the growth rate:
    • it appears in the dilution term of proteins (easy to incorporate into the model)
    • but is also affects in a highly nonlinear way the gene copy number
    • it affects the concentration of free and bound RNAp and therefore the level of transcription etc..
  • some modules in your system may be very hard to model (if at all possible)
    • for instance transport of molecules through a membrane and diffusion phenomena can be modelled but it becomes complicated fast
    • in a model, errors pile up so much so that after a while the predictive power of your model is negligible.
  • your synthetic pathway may 'cross-talk' with natural pathways; since we are not able to model the whole metabolism of the cell this crosstalk effect can not be assessed.


Fortunately, designing simple pathways with predictable properties/functions is possible, even without the extensive use of software. This session aims at introducing to you the basic tools and techniques of design (without which no computer-assisted design is possible. Remember: in practice it gets very complicated, very fast...





Preliminary Simplifications



Download this File, and Open it with CellDesigner.


Model CellDesigner Instructions

It can be shown that after some normalisation the ODE system can be written as:

[math]\displaystyle{ \begin{alignat}{1} \frac{d[mRNA]_{i}}{dt} & = \frac{a}{1+{[Protein]_{j}}^n} - [mRNA]_{i} \\ \frac{d[Protein]_{i}}{dt} & = b[mRNA]_{i} - b[Protein]_{i} \\ \ i=1,2,3; \\ \ j=3,1,2; \\ \end{alignat} }[/math]


Repressilator Genetic Circuit



Negative Auto-Regulation


Model The basic delay

It can be shown that after some normalisation the ODE system can be written as:

[math]\displaystyle{ \begin{alignat}{1} \frac{d[mRNA]}{dt} & = \frac{a}{1+{[Protein]}^n} - [mRNA] \\ \frac{d[Protein]}{dt} & = b[mRNA] - b[Protein] \\ \end{alignat} }[/math]


Repressilator Genetic Circuit

The following questions constiture the last part of your coursework (Section E):