Team:CIDEB-UANL Mexico/math resistance

From 2014hs.igem.org

Revision as of 06:19, 16 June 2014 by Oda.ibarra (Talk | contribs)

iGEM CIDEB 2014 - Project

Resistance Module

The resistance module is based in the use of IrrE gene. IrrE is a protein which is not affected by external factors during its transcription as well as its translation, so we needed to use the established parameters but with the data we obtained from IrrE.


\begin{equation} \large \frac{d\left [ mRNA \right ]}{dt}=\alpha_{1}-d_{1}\left [ mRNA \right ] \end{equation}

We used the parameters for translation and transcription rate from Singapore 2008 iGEM team as well as the speeds at which E. coli carry out transcription and translation assuming a transcription speed of 70nt/s and a translation speed of 40aa/s. So we used them in the equations below with the IrrE gene length (986nt) and protein length (311aa) respectively.


\begin{equation} \large \alpha_{1}=\frac{transcription speed}{gen lenght \cdot(nt)} \alpha_{2}=\frac{translation speed}{protein lenght \cdot(aa)} \end{equation}
\begin{equation} \large \alpha_{1}=\frac{(70)(60)}{986}=4.25\alpha_{2}=\frac{(40)(60)}{311}=7.7 \end{equation}

Then, we used the parameters for degradation rates for proteins and mRNAs from Beijing PKU 2009 iGEM team:


\begin{equation} \large d_{1}=\frac{1}{half-life}+\frac{1}{30}min \end{equation}
\begin{equation} \large d_{2}=\frac{1}{half-life}+\frac{1}{30}min \end{equation}

But since IrrE half-life has not been determined yet, we decided to search about homologous proteins with the same function as IrrE (transcriptional factor) and according to Dibden and Green (2005) the average half-life for transcriptional factors in E. coli is 45 min. They tested FNR proteins (transcriptional factors) through thermo-induciblefnr expression observing that their half-life was 45min average.

For determining the degradation rate of average mRNA we used the information from Selinger’s team (2003). They carried several experiments for finding average mRNA half-life in E. coli. They used mRNAs about 1100nt concluding they have an average half-life of 5min. So with this we found the average mRNA half-life of IrrE was 4.45min.


\begin{equation} \large HL=\frac{1100\cdot(bp)}{5min} \end{equation}

With all these information we could find the degradation rates for both transcription and translation of IrrE.


\begin{equation} \large d_{1}=\frac{1}{4.5}+\frac{1}{30}=0.25d_{2}=\frac{1}{45}+\frac{1}{30}=0.055 \end{equation}

For the simulation we used Simbiology using the previous data in the equations for finding the amount of proteins E. coli would produce at certain time. The following was the result of the simulation:


IMG_0317


But for translation there was another factor we had to consider, the “fpost(p) which were the posttranslational variables affecting the production of the functional protein.


\begin{equation} \large \frac{dpi}{dt}=Ti-Dpmi- f_{post (p)} \end{equation}

We found that IrrE needs to have a positive charge to be functional, accepting one Zn2+ ion from E.coli (Vujicic, 2009). Vujicic’s team developed the structure of IrrE deducing it should have three domains; one specific is the zinc-binding site where the Zn2+ ion binds to make IrrE positive. It becomes positive what makes possible its binding to a substrate forming a substrate complex, but the substrate is unknown. For that reason we only take as a “fpost(p) variable the Zn2+ binding because the data for finding the affinity between IrrE and the unknown substrate was not possible to determine.

As we did not have the affinity between IrrE and Zn2+ we need to find the data from homologous proteins. In fact, we found that proteins with zinc binding domains have an affinity between 0.1 and 0.2 in E. coli (Vorackova, 2012). We used 0.15 as average in the equation for association constant which is defined as the following:


\begin{equation} \large Ka=\frac{\left [ C \right ]}{\left [ S \right ]\left [ E \right ]} \end{equation}

Where “[C]” is the complex formed, “[S]” is the substrate and “[E]” is the enzyme, ligand or ion. Substituting for IrrE it is expressed as below:


\begin{equation} \large Ka=\frac{\left [ C \right ]}{\left [ IrrE \right ]\left [ Zn+ \right ]}0.15=\frac{\left [ C \right ]}{\left [ IrrE \right ]\left [ Zn+ \right ]} \end{equation}

With the association constant we could use it in Simbiology for modelling the functional IrrE production. The results were the following:

IMG_0317


IMG_0317

When we compared both graphs (Graph 1 and Graph 2) we realized that not all the IrrE production was functional. In fact from the total amount produce (about 2400), only about 1700 are functional. So it demonstrates that the rate at which Zn2+ ions binds to IrrE is slower than the rate at which IrrE is produced, leaving a nonfunctional IrrE.


Bibliography

● David Dibden, J. G. (2005). In vivo cycling of the Escherichia coli transcription factor FNR between active and inactive states. Microbiology, 4063-4070.

● Douglas Selinger, R. M. (2003). Global RNA Half-Life Analysis in Escherichia coli Reveals Positional Patterns of Transcript Degradation. Genome Research, 216-223.

● Vorackova Irena, S. S. (2011). Purification of proteins containing zinc finger domains using Immobilized Metal Ion Affinity Chromatography. Protein Expression and Purification, 88-95.


iGEM CIDEB 2014 - Footer