LMB562

From CSBLwiki

Revision as of 04:49, 12 October 2012

Class Description

This is introductory class to cover topics in synthetic biology by its past (origin), present and future.

The course can be divided into three sessions

1. Fundamental principles of synthetic biology
2. Landmark papers in synthetic biology
3. Current topics and future direction

Additional topics - ethical issues and risks in synthetic biology

1. iGEM
2. DIYBio

Course materials

Fundamental principles

1. Error fetching PMID 16306983: [DrewEndy1]
2. Error fetching PMID 16711359: [SciAm1]
3. Error fetching PMID 18724274: [Embo]

Landmarks

4. Error fetching PMID 10659857: [Toggle]
5. Error fetching PMID 10659856: [Repress]
6. Error fetching PMID 16306980: [Voigt]
7. Error fetching PMID 16306982: [Elo]
8. Error fetching PMID 15858574: [Weiss]
9. Error fetching PMID 16612385: [Keas]
10. Error fetching PMID 16330045: [Chris]

PLoS collections

• PLoS ONE organized Synthetic Biology Collections

Nature WEB focus

• Synthetic systems biology

Science

• Special issue for synthetic biology

PubMed database

• synthetic biology - search results

Google scholar

• synthetic biology - search results

Open questions

1. What is 'synthetic biology'?
2. How is synthetic biology different with traditional genetic engineering?

Practice

library(deSolve)

### Repressilator
### Elowitz & Leibler (2000)
repress = function(t,y,p) {
  dy = rep(0,6);
  dy[1] = -y[1] + p[1]/(1.+y[4]^p[4])+p[2];
  # dy[1] lacI; dy[4] CI
  dy[2] = -y[2] + p[1]/(1.+y[5]^p[4])+p[2];
  # dy[2] tetR; dy[5] LacI
  dy[3] = -y[3] + p[1]/(1.+y[6]^p[4])+p[2];
  # dy[3] cI; dy[6] TetR
  dy[4] = -p[3]*(y[4]-y[3]); # cI
  dy[5] = -p[3]*(y[5]-y[1]); # lacI
  dy[6] = -p[3]*(y[6]-y[2]); # tetR
  return(list(dy))
  # p[1..4] = c(alpha,alpha0,beta,hill)
}
# initial parameters
x0 = 2*runif(6); p = c(50,0,.2,2); t = seq(0,500,by=0.1)
out = lsoda(x0,t,repress,p)
matplot(out[,1],out[,2:4],type="l")

### Genetic toggle switch
### Gardner & Collins (2000)
Toggle=function(t,y,parms) {
  u=y[1]; v=y[2];
  du= -u + parms[1]/(1+v^parms[2]);
  dv= -v + parms[1]/(1+u^parms[3]);
  dY=c(du,dv);
  return(list(dY));
}

x0=c(.2,.1); times=seq(0,50,by=0.2)
out=lsoda(x0,times,Toggle,parms=c(3,2,2))
matplot(out[,1],out[,2:3],type="l",ylim=c(0,3),xlab="time t",ylab="u,v");
# phase portrait
plot(out[,2],out[,3],type="l",xlab="u",ylab="v")