Team:Michigan/Modeling

From 2013.igem.org

(Difference between revisions)
(7. Optimization)
Line 182: Line 182:
=Future Directions=
=Future Directions=
Data!
Data!
-
<!-- *** What falls between these lines is the Alert Box!  You can remove it from your pages once you have read and understood the alert *** -->
 
-
<html>
 
-
<div id="box" style="width: 700px; margin-left: 137px; padding: 5px; border: 3px solid #000; background-color: #fe2b33;">
 
-
<div id="template" style="text-align: center; font-weight: bold; font-size: large; color: #f6f6f6; padding: 5px;">
 
-
This is a template page. READ THESE INSTRUCTIONS.
 
-
</div>
 
-
<div id="instructions" style="text-align: center; font-weight: normal; font-size: small; color: #f6f6f6; padding: 5px;">
 
-
You are provided with this team page template with which to start the iGEM season.  You may choose to personalize it to fit your team but keep the same "look." Or you may choose to take your team wiki to a different level and design your own wiki.  You can find some examples <a href="https://2008.igem.org/Help:Template/Examples">HERE</a>.
 
-
</div>
 
-
<div id="warning" style="text-align: center; font-weight: bold; font-size: small; color: #f6f6f6; padding: 5px;">
 
-
You <strong>MUST</strong>  have all of the pages listed in the menu below with the names specified.  PLEASE keep all of your pages within your teams namespace. 
 
-
</div>
 
-
</div>
 
-
</html>
 
-
 
-
<!-- *** End of the alert box *** -->
 

Revision as of 01:51, 28 September 2013

Xayona Website Template

Best model.jpg

Contents

Introduction

Switch modeling

• produce both states of switch

• use model to optimize switch

Mass action modeling

Definition

• differential equations

• each equation is a sum of rates

• each rate is proportional to each of its reactants

Benefits

• can minimize assumptions

• already used in chemical rate equations

• standardizes equations:

Mass action.png
    	• x is the species vector
    	• a is an amount of rates changing x
    	• b is a rate constant
    	• c is a number of different species involved
    	• d is a species involved
    	• e is a number of each species involved	
    	• f is a number of species d involved
    	• g is a highest species involved
    	• l is a rate order
    	• m is the highest number of rates involved
    	• n is the number of species

Analytical modeling

• benefits

   • given parameters, it precisely predicts behavior
   • accurately predicts the data needed to determine parameters (unlike numerical analysis)
   • allows avoidance of numerical errors like rounding
   • naturally standardized

ERSESCO Algorithm

The ERSESCO algorithm the Michigan team developed provides a clearly defined method of using symbolic analysis on systems. The algorithm has seven major steps:

1. Equation

Formulate the system into a system of differential equations.

2. Reduction

Apply constraints, boundary conditions, and conservation laws to reduce the system of equations to fewer variables.

3. Solution

Solve the reduced system of equations, if possible.

4. Equilibration

Find the steady states of the system by setting differential equations to zero and solving for critical points.

5. Stabilization

Evaluate the Jacobian matrix of the system of equations at each critical point and find the corresponding eigenvalues. Eigenvalues with a negative real part are stable critical points.

6. Calibration

Solve for the parameters of the system of equations in terms of the system's initial, steady state, time-dependent variables. This step allows the system parameters to be easily characterized using experimental data.

7. Optimization

Optimize the steady state, eigenvalues, parameters, or other desired quantity by equating the derivative of that quantity with respect to each of the system's parameters with zero. Solving these equations gives the optimization conditions for that quantity.

SimBiology

• is a MatLab package provided to iGEM teams

• provides a numerical approximation

• easily creates from a diagram

Model of Recombinase Expression

The Recombinase Expression Model describes the synthesis and degradation of any protein produced at a relatively constant rate. Here, the model is used in particular to describe the synthesis and degradation reaction rates ksyn and kdeg of recombinases in the cell. The model also predicts the recombinase concentration.

Model of the Switch

The Switch Model describes how the switch equilibrates when the concentrations of the recombinases FimE and Hbif approach constant values. This model assumes that each recombinase catalyzes a one way reaction. This model predicts the degree of cooperativity n and m in the binding reaction of the recombinase to DNA for FimE and Hbif, respectively. Furthermore, this model predicts the forward catalysis rate constants, kF and KH, for FimE and Hbif.

Model of Inducible Hbif

The Inducible Hbif Model describes how the switch flips when acted upon by the Lux/HSL-controlled expression of Hbif. In this model, [HSL] represents the concentration of the species HSL:pLux:LuxR, the complex that promotes the expression of Hbif. In this model [HSL] is a function of pLux and LuxR expression levels in the cell, the amount of HSL added, and the fractional occupancies at the complexation equilibrium. It can be assumed that this complexation reaction happens instantaneously with respect to the slow, rate-limiting translation step of Hbif. This model predicts the rate k1 of Lux-induced expression of Hbif, the degradation rate kdeg of Hbif, the forward and reverse rates k2 and k-2 of Hbif catalysis on the switch, and the equilibrium constant K2 for the switch when acted upon by Hbif.

Expression Model

ExpressionModelFigure.png

Equation:

ExpressionModel1.png

Reduction: None

Solution:

ExpressionModel2.png

Equilibration:

ExpressionModel3.png

Stabilization:

ExpressionModel3-1.png

Calibration:

ExpressionModel4.png
ExpressionModel5.png

Optimization:

ExpressionModel6.png

SimBiology:

Expression.jpg

Expression Model Derivation

Switch Model

SwitchModelFigure.png

Equation:

SwitchModel1.png
SwitchModel2.png

Reduction:

SwitchModel3.png

Solution:

SwitchModel3-1.png

Equilibration:

SwitchModel4-1.png

Stabilization:

SwitchModel4-2.png

Calibration:

SwitchModel5-1.png
SwitchModel5-2.png

Optimization:

SwitchModel6.png

SimBiology:

Switch.jpg

Switch Model Derivation

Inducible Hbif Model

InducibleHbifModelFigure.png

Equation:

InducibleHbifModel1.png
InducibleHbifModel2.png

Reduction:

InducibleHbifModel3.png

Solution:

InducibleHbifModel4.png

Equilibration:

InducibleHbifModel5.png

Stabilization:

InducibleHbifModel5-2.png

Calibration:

InducibleHbifModel6-1.png
InducibleHbifModel6-2.png
InducibleHbifModel6-3.png

Optimization of Hbif:

InducibleHbifModel7-1.png

Optimization of ON:

InducibleHbifModel7-2.png

SimBiology:

InducedHbif.jpg

Inducible Hbif Model Derivation

Future Directions

Data!


Home Team Official Team Profile Project Parts Submitted to the Registry Modeling Notebook Safety Attributions