Team:UC Berkeley/Modeling

From 2008.igem.org

Revision as of 05:54, 25 October 2008 by ChrisBrown (Talk | contribs)


The kinetics of the Lambda phage lysis device was modeled to help our team gain an insight into the behavior of our system. The below equations describe our system


Cdb1.jpg

Where P_H and P_AH represent the mRNA promoter strengths, gamma_mRNA,H and gamma_mRNA,AH are the degradation rates for holin and antiholin respectively and gamma is the protein degradation rate, k_H , k_AH and k_c represent the rate constants for holin and antiholin formation and the coupling rate for the holin-antiholin dimer.

At steady state,

Cdb2.jpg

The system can be simplified into the following transfer function

Cdb3.jpg

Where the system can be divided into three dimensionless parameters which describe the behavior of holin, antiholin and the holin-antiholin dimer

Cdb4.jpg

Using MatLab, the following graph of the system was produced.

Image1.png

The literature indicates that at the time of lysis, cells infected with lambda phage have approximately 1000 holin proteins. Therefore, the critical concentration of holin (Hc) was set at 1000 holin proteins per cell.

The horizontal line at y=1 represents the critical holin concentration needed induce lysis.

Image2.png Image3.png