"

From 2013.igem.org

You do not have permission to edit this page, for the following reason:

You can view and copy the source of this page:

{{:Team:TU-Eindhoven/Template:Base}} {{:Team:TU-Eindhoven/Template:MenuBar}} =Integral Model= {{:Team:TU-Eindhoven/Template:Lead}} The last step in creating a comprehensive model of the contrast mechanism is merging all individual models. [[Team:TU-Eindhoven/ODEModel | Previously]] the model to predict the FNR concentration in a cell was described. This model will be combined with the [[Team:TU-Eindhoven/StochasticModel | decoy sites model]], which was extended to contain the [[Team:TU-Eindhoven/FNRPromotor | specific promoter]] that was used in this research. Here, the properties of the merged model will be discussed. {{:Team:TU-Eindhoven/Template:LeadEnd}} ==Merging the Models== In order to merge the FNR model and the decoy sites model, the exact connection between the two models had to be defined. The FNR model describes the (active) FNR concentration in the bacterial cells. This active FNR species fulfills the role of transcription factor in the decoy sites model. Using this information, the models could simply be combined by merging the ''Active FNR'' species of the FNR model and the ''T'' (transcription factor) of the decoy sites model. The resulting model schematics are depicted in {{:Team:TU-Eindhoven/Template:Figure | id=IntegralModelSchematics }}. {{:Team:TU-Eindhoven/Template:Float | position=right | size=12 }} {{:Team:TU-Eindhoven/Template:Image | filename=IntegralModel.jpg}} {{:Team:TU-Eindhoven/Template:FloatEnd | caption=Schematic of the integral model | id=IntegralModelSchematics }} ==Results== This model was implemented as both a ODE and a stochastic model, various parameters were estimated as described on the [[Team:TU-Eindhoven/FNRPromotor| FNR promoter page]]. One of the most interesting predictions that can be derived from the integral model is the influence of adding decoy sites to the DNA on the response to oxygen. {{:Team:TU-Eindhoven/Template:Figure | id=DecoySitesOxygen}} and {{:Team:TU-Eindhoven/Template:Figure | id=DecoySitesOxygenResponse}} endeavor to visualize this prediction as clearly as possible. {{:Team:TU-Eindhoven/Template:Float | position=left | size=6 }} {{:Team:TU-Eindhoven/Template:Image | filename=DecoySitesOxygen.png}} {{:Team:TU-Eindhoven/Template:FloatEnd | caption=Correlation between the expression rate and the oxygen concentration for different amounts of decoy sites. | id=DecoySitesOxygen }} {{:Team:TU-Eindhoven/Template:Float | position=right | size=6 }} {{:Team:TU-Eindhoven/Template:Image | filename=DecoySitesOxygenResponse.png}} {{:Team:TU-Eindhoven/Template:FloatEnd | caption=Reaction of the expression rate to a change in oxygen concentration. | id=DecoySitesOxygenResponse }} From {{:Team:TU-Eindhoven/Template:Figure | id=DecoySitesOxygenResponse}} it can be concluded that the sharpest response is realized using no decoy sites. Looking at {{:Team:TU-Eindhoven/Template:Figure | id=DecoySitesOxygen}} it can be seen that adding decoy sites decreases the steady state expression rate, which is undesired. For a sharper response a weaker promoter would be needed in combination with stronger decoy sites, of which results are shown in {{:Team:TU-Eindhoven/Template:Figure | id=DecoySitesOxygenSharp}}. {{:Team:TU-Eindhoven/Template:Float | position=left | size=6 }} {{:Team:TU-Eindhoven/Template:Image | filename=DecoySitesOxygenSharp.png}} {{:Team:TU-Eindhoven/Template:FloatEnd | caption=Correlation between the expression rate and the oxygen concentration for different amounts of decoy sites using a weaker promoter (colors have the same meaning as in {{:Team:TU-Eindhoven/Template:Figure | id=DecoySitesOxygen}}). | id=DecoySitesOxygen }} It can be seen that the response is sharper and more switch like, but the problem is that the expression rate decreases significantly. So using the weaker promoter would not allow for generating satisfactory MRI contrast, which is the main purpose of the project. Combining these observations, adding decoy sites to the construct appears to only have negative effects and the idea should be dismissed, based on the modeling results. Another conclusion that can be drawn from the models is that the expression rate needs around one hour to reach its peak, so any observable changes in protein concentration will be delayed by a little more than one hour. ==Sensitivity Analysis== A sensitivity analysis is carried out over the integral model. It aims to see the effect of uncertainties in different parameters on the result curve. Therefore we know which parameters are crucial control points. In the sensitivity analysis, the parameters are added with a random error and put into the model. The method of Latin hypercube sampling is used to generate the error so that the error distribution is uniform over the region. As the result({{:Team:TU-Eindhoven/Template:Figure | id=sensitivityAnalysis}}) shows, the expression rate is most sensitive to the FNR concentration and the expression rate. It is also sensitive to <html>$\alpha_{21}$</html>, <html>$\beta_1$</html> and <html>$\beta_{31}$</html> under anaerobic condition. These parameters can be found in the [[Team:TU-Eindhoven/ODEModel | FNR model ]]. <br> <br> {{:Team:TU-Eindhoven/Template:Float | position=left | size=6 }} {{:Team:TU-Eindhoven/Template:Image | filename=sensitivityAnalysis.png}} {{:Team:TU-Eindhoven/Template:FloatEnd | caption=The Result of Sensitivity Analysis. | id=sensitivityAnalysis }} ==EGFP Estimation== In order to link to the result from [https://2013.igem.org/Team:TU-Eindhoven/AnearobicTesting#Fluorescence_Measurements wetlab], the integral model needs to be adapted for EGFP. The two parameters, the absolute expression rate and the degradation rate of EGFP, are needed for the model. The degradation of EGFP is an exponential function and the half-life time of EGFP is about 26 hours{{:Team:TU-Eindhoven/Template:Ref | id=EGFPdegra | author=Pete Corish and Chris Tyler-Smith | title=Attenuation of green fluorescent protein half-life in mammalian cells | journal=Protein Engineering | edition=12 | pages=1035-1040 | year=1999 }}. The expression rate of EGFP is linear to the [https://2013.igem.org/Team:TU-Eindhoven/StochasticModel#Protein_Expression_Rate promoter occupancy] which is known from the model. The ratio between the absolute expression rate and the relative expression rate will be estimated from the lab data. <!--result of estimation--> '''The source code of all models can be found [[Team:TU-Eindhoven/Code:Models | here]].''' {{:Team:TU-Eindhoven/Template:BaseFooter}} {{:Team:TU-Eindhoven/Template:Sponsors}} {{:Team:TU-Eindhoven/Template:SetTitle | menu=drylab | page=Integral Model }} {{:Team:TU-Eindhoven/Template:UseReferencing}} {{:Team:TU-Eindhoven/Template:UseFigures}} {{:Team:TU-Eindhoven/Template:UseEquations}} {{:Team:TU-Eindhoven/Template:SetHeader | nr=4}}

Return to Team:TU-Eindhoven/IntegralModel.