Team:TJU China/Model

<!DOCTYPE >

Dynamic Model of Heavy Metal Detection Biosensor
Minghui Yin,Sherry Dongqi Bao
TianJin University
October 15,2018
1 Introduction
Modeling is a powerful tool in synthetic biology. It provides us with a necessary engineering approach to characterize our pathways quantitatively and predict their performance,thus help us test and modify our design.Through the dynamic model of heavy-metal detection biosensor,we hope to gain insights into the characteristics of our whole circuit's dynamics.
2 Methods
2.1 Analysis of metabolic pathways
Figure 1: Metabolic pathways related to plasmid#1
At the beginning, on the plasmid#1, the promoter $P_{arsR}$ isn't bound with ArsR,thus it is active.ArsR and smURFP are transcribed and translated under the control of the promoters $P_{arsR_{u}}$ and $P_{arsR_{d}}$,with subscript u and d representing upstream and downstream separately.The subscript l of smURFP in the equation means leaky expression without the expression of $As^{3+}$.As ArsR is expressed gradually,it will bind with the promoter $P_{arsR}$ and make it inactive.[1]
On the plasmid#2,the fusion protein of dCas9 and RNAP(RNA polymerase) are produced after transcription and translation,and sgRNA is produced after transcription.
Figure 2: Metabolic pathways related to dCas9/RNAP
dCas9(*RNAP) can bind with its target DNA sequence without cutting, which is at the upstream of the promoter $P_{arsR_{d}}$.Simulataneously,dCas9 can lead RNAP to bind with the promoter $P_{arsR_{d}}$ and enhance the transcription of smURFP.However,because the promoter $P_{arsR_{d}}$ has already bound with ArsR,as a result,RNAP can't bind with the promoter $P_{arsR_{d}}$. can’t bind with the promoter $P_{arsR_{d}}$.
However,at the presence of $As^{3+}$,it can bind with ArsR,then dissociate ArsR and $P_{arsR_{d}}$ , which makes the combination of RNAP and $P_{arsR_{d}}$ possible.
We then take degradation into account:
2.2 Analysis of ODEs
Applying mass action kinetic laws,we obtain the following set of differentiak equations.The several complexes involved:Ars$R^*$$P_{arsR}$,$As^{3+}$,${dCas9}^*$RNAP,${dCas9}^*$RNAP:sgRNA,${dCas9}^*$RNAP:${sgRNA}^*P_{arsR}$, are respectively abbreviated as $cplx_1$,$cplx_2$,$cplx_3$,$cplx_4$,$cplx_5$.
2.3 Simulation
Our simulation is based on two softwares: MATLAB (SimBiology Toolbox) and COPASI.
SimBiology Toolbox provides functions for modeling,simulating and analyzing biochemical pathways by the powerful computing engine of MATLAB.
Figure 3:Reaction map generated from the reaction sets above by SimBiology Toolbox
Figure 4:Simulation of smURFP production as a function of time by MATLAB Through the figure, we can see that the smURFP can gradually increase and reach a steady state after a period in the presence of arsenic ions.
2.4 Sensitivity
A good biosystem should have certain stability towards fluctuations in parameters.A good model should reflect this,and hence a test for robustness can be essential to the model.
Robustness analysis can also pinpoint which reactions/parameters that are important for obtaining a specific biological behavior.A simple measure for sensitivity is to measure the relative change of a system feaure due to a change in a parameter.As for our model,the feature can be the equilibrium concentration of the smURFP(C) for which the sensitivity(S) to a parameter k is:
After analysis, we found that the concentration of smURFP is relatively sensitive to parameters such as ktx3,ktl3,ktx4,kb4,kb6,kd2,kd5, kd6,kd7,kd8,kd11, etc. Among these parameters, except for the parameters that directly affect the production and degradation of smURFP,the rest of them are all related to dCas9-RNAP:sgRNA. It shows that our model reflects the critical role of dCas9-RNAP:sgRNA,which initially confirms our hypothesis:dCas0-RNAP can enhance transcription to increase the concentration of smURFP. However, due to the lack of previous modeling studies on dCas9-RNAP,some kinetic parameters may not be very accurate,and due to time limitation,we have not implemented experiments to measure related parameters,which may lead to some deviations in our model.
The sensitivity of each parameter is shown in the figures below.
(a)sensitivity of ktx1 (b)sensitivity of ktl1
(c)sensitivity of ktx2 (d)sensitivity of ktl2
(e)sensitivity of ktx3 (f)sensitivity of ktl3
(g)sensitivity of ktx4 (h)sensitivity of kb1
(i)sensitivity of kb2 (j)sensitivity of kb3
(a)sensitivity of kb4 (b)sensitivity of kb5
(c)sensitivity of kb6 (d)sensitivity of kd1
(e)sensitivity of kd2 (f)sensitivity of kd3
(g)sensitivity of kd4 (h)sensitivity of kd5
(a)sensitivity of kd6 (b)sensitivity of kd7
(c)sensitivity of kd8 (d)sensitivity of kd9
(e)sensitivity of kd10 (f)sensitivity of kd11
Note:The ordinate axis represents the sensitivity S,and the abscissa axis is the parameter k for which we want to evaluate the sensitivity.
2.5 Application of the model
Since the goal of our project is to increase the sensitivity of biosensors by introducing a complex of dCas9-RNAP and sgRNA, and one of the purposes of our model is to explore whether this complex is effective.So we assume a reasonable and large enough concentration value for this complex. We use the concentration of Glyceraldehyde-3-phosphate dehydrogenase A as the assumed concentration.Glyceraldehyde-3-phosphate dehydrogenase A(gapA) is a crucial enzyme in the glycolytic pathway,and the gene encoding this enzyme is a housekeeping gene in E.coli cells with high expression levels.We find in the literature that the protein mass of gapA is 48645 fg/cell,and its molecular weight is 35492 Da.[4] The amount of abundance of Glyceraldehyde-3-phosphate dehydrogenase A protein per cell can be calculated as follows:
As for the size of E.coli,we found relevant data from the literature,as the figure below shows.[5]
Figure 8:Size of E.coli
The volume of E.coli can be calculated as follows:
Then the concentration of Glyceraldehyde-3-phosphate dehydrogenase A protein in the cell can be determined:
With this concentration,we can get very nice results:
Figure 9:smURFP production with enough dCas9-RNAP:sgRNA
Compared to the diagram without introducing dCas9-RNAP:sgRNA:
Figure 10:smURFP production within a reasonable time frame
Figure 11:smURFP production reached equilibrium but it takes a long time
From these three figures, we can conclude that dCas9-RNAP:sgRNA does have the effect of promoting transcription and increasing fluorescence intensity,thereby increasing sensitivity,as long as its concentration is sufficient.This result enhances the confidence of the experimental group,and they need to try to improve the expression of dCas9-RNAP:sgRNA in E.coli without having to doubt its role.
References
[1] LA Pola-Lopez et al."Novel arsenic biosensor "POLA" obtained by a genetically modified E.coli bioreporter cell" .In:Sensors and Actuators B:Chemical254(2018),pp.1061-1068.
[2] Yves Berset et al."Mechanistic Modeling of Genetic Circults for ArsR Arsenic Regulation".In:ACS synthetic biology 6.5(2017),pp.862-874.
[3] Eyal Karzbrun et al."Coares-grained dynamics of protein synthesis in a cell-free system".In:Phtsical review letters 106.4(2011),p.048104.
[4] Yasushi Ishihama et al."Exponentially modified protein abundance index(emPAI) for estimation of absolute protein amount in proteomics by the number of sequenced peptides per protein".In:Molecular E Cellular Proteomics 4.9(2005),pp.1265-1272.
[5] Nili Crossman,Eliora Z Ron,and Conrad L Woldringh."Changes in cell dimensions during amino acid starvation of Escherichia coli."In:Journal of bacteriology 152.1(1982),pp.35-41.