Team:Michigan/Modeling

From 2013.igem.org

(Difference between revisions)
(Mass action modeling)
(Mass action modeling)
 
(78 intermediate revisions not shown)
Line 1: Line 1:
{{:Team:Michigan_Template}}
{{:Team:Michigan_Template}}
-
[[File:Best_model.jpg|900px|center]]
 
=Introduction=
=Introduction=
===Switch modeling===
===Switch modeling===
-
• produce both states of switch
+
The goal for our modeling was to create rigorous standard procedures for mathematically representing and analyzing systems.  Using our project as an example, we worked to demonstrate the viability of our approach by calculating the switch positions and optimum parameters.  This standardization should be useful in combination with system specific methods to create better modeling.
-
 
+
-
• use model to optimize switch
+
===Mass action modeling===
===Mass action modeling===
 +
Mass action is a form of modeling that uses systems of homogenous differential and difference equations.  Each equation is a sum of rates and each rate is proportional to the species involved.  The advantage is its straightforward uniform formulation, its minimal assumptions, and its flexibility of complexity.  Mass action itself uses very few assumptions but it can be used to represent systems of varying degrees of simplicity.  The fundamentals of chemical kinetics already champion the use of this modeling in its use of rate equations.  There chemical systems are completely defined by the products, reactants, and reaction rate constants.  This allows chemical modelling to be highly standardized and hence more accessible.  It's actually the uniform nature of mass action that's responsible.  Most importantly for iGEM, all mass action models can actually be put in a standard form as follows:
 +
[[File:mass action.png|670px|center]]
 +
• x is a species
-
Definition
+
• a is an amount of rates changing x
-
differential equations
+
b is a rate constant
-
each equation is a sum of rates
+
c is a number of different species involved
-
each rate is proportional to each of its reactants
+
d is a species involved
-
Benefits
+
• e is a number of each species involved
-
can minimize assumptions
+
f is a number of species d involved
-
already used in chemical rate equations
+
g is a highest species involved
-
standardizes equations and their analysis:
+
l is a rate order
-
<html><iframe src="https://docs.google.com/file/d/0B7Cfp8sKyoyjNE5OdEhsQ3Nrdnc/preview" width="100%" height="480"></iframe></html>
+
-
===Analytical modeling===
+
• m is the highest number of rates involved
-
benefits
+
n is the number of species
-
    • naturally standardized
+
===Analytical modeling===
-
   
+
====Benefits====   
-
    • given parameters, it precisely predicts behavior
+
Given the parameters of the system, analytical (symbolic) modeling precisely predicts a system's behavior. Furthermore, unlike numerical methods, analytical modeling accurately predicts the data needed to determine parameters and avoids the inherent numerical rounding and approximation errors of more crude analyses.  Analytical modeling naturally lends itself to standardization because the results are unique to any given model.  Correct solutions, equilibria, eigenvalues, etc. will always have the same expression regardless of who or what calculated them.  The same can't be said for numerical methods.  With this and the iGEM community in mind, an analytical standard was developed that will allow for the exchange and automation of future modeling.
-
    • accurately predicts the data needed to determine parameters
+
====ERSESCO Algorithm====
 +
The ERSESCO algorithm the Michigan team developed provides a clearly defined method for analyzing mathematical systems. The algorithm has seven major steps:
-
    • allows avoidance of numerical errors like rounding
+
'''1. Equation:'''
 +
Formulate the system into a system of differential equations.
-
• ERSESCO
+
'''2. Reduction:'''
 +
Apply constraints, boundary conditions, and conservation laws to reduce the system of equations to fewer variables.
-
    • equation
+
'''3. Solution:'''
 +
Solve the reduced system of equations, if possible.
-
    • reduction
+
'''4. Equilibration:'''
 +
Find the steady states of the system by setting differential equations to zero and solving for critical points.
-
    • solution
+
'''5. Stabilization:'''
 +
Evaluate the Jacobian matrix of the system of equations at each critical point and find the corresponding eigenvalues. To get stability, find the condition for which all eigenvalues have a negative real part.
-
    • equilibration
+
'''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.
-
    • stabilization
+
'''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.
-
    • calibration
 
-
 
-
    • optimization
 
===SimBiology===
===SimBiology===
-
numerical approximation
+
SimBiology is a useful MATLAB package provided to all iGEM teams. It can be used as a tool to quickly run numerical approximations of biological systems whether or not data for system parameters are known. SimBiology can quickly produce a model from a simple schematic diagram. As a result, the package can quickly provide a qualitative analysis of a system with little work from the user.
-
• diagram / schematic
+
In the following models, our team applied our ERSESCO algorithm to fully analyze the systems at hand. Finally a SimBiology simulation was used as qualitative verification. We present the results from each step of the ERSESCO algorithm, followed by a SimBiology schematic with a graph showing the general behavior of the system, followed by a link to a more detailed derivation of the models.
-
===Model of Recombinase Expression===
 
-
The Recombinase Expression Model describes the synthesis and degradation of any protein in the cell. 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.
+
=Expression Model=
-
===Model of the Switch===
+
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.
-
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===
+
[[File:ExpressionModelFigure.png|center]]
 +
'''Equation:'''
 +
[[File:ExpressionModel1.png|218px|center]]
-
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=
+
'''Reduction: None'''
-
[[File:ExpressionModelFigure.png|center]]
 
-
Equation:
 
-
[[File:ExpressionModel1.png|218px|center]]
 
-
Reduction: None
 
-
Solution:
+
'''Solution:'''
[[File:ExpressionModel2.png|343px|center]]
[[File:ExpressionModel2.png|343px|center]]
-
Equilibration:
+
 
 +
 
 +
'''Equilibration:'''
[[File:ExpressionModel3.png|128px|center]]
[[File:ExpressionModel3.png|128px|center]]
-
Stabilization:
+
 
 +
 
 +
'''Stabilization:'''
[[File:ExpressionModel3-1.png|181px|center]]
[[File:ExpressionModel3-1.png|181px|center]]
-
Calibration:
+
 
 +
 
 +
'''Calibration:'''
[[File:ExpressionModel4.png|269px|center]]
[[File:ExpressionModel4.png|269px|center]]
[[File:ExpressionModel5.png|322px|center]]
[[File:ExpressionModel5.png|322px|center]]
-
Optimization:
+
 
 +
 
 +
'''Optimization:'''
[[File:ExpressionModel6.png|73px|center]]
[[File:ExpressionModel6.png|73px|center]]
-
SimBiology:
+
 
 +
 
 +
'''SimBiology:'''
[[File:Expression.jpg|center]]
[[File:Expression.jpg|center]]
-
[https://docs.google.com/document/d/18ur0oiT7_mHhmd7Nk9vDJtsAGuFecNA1CtgTz9HAjsA/edit?usp=sharing Expression Model Derivation]
+
 
 +
'''Derivation:''' [https://docs.google.com/document/d/18ur0oiT7_mHhmd7Nk9vDJtsAGuFecNA1CtgTz9HAjsA/edit?usp=sharing Expression Model]
=Switch Model=
=Switch Model=
 +
 +
The Switch Model describes how the switch equilibrates when the concentrations of the recombinases FimE and Hbif approach constant values.  The 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.
 +
[[File:SwitchModelFigure.png|center]]
[[File:SwitchModelFigure.png|center]]
-
Equation:
+
 
 +
 
 +
'''Equation:'''
[[File:SwitchModel1.png|395px|center]]
[[File:SwitchModel1.png|395px|center]]
[[File:SwitchModel2.png|280px|center]]
[[File:SwitchModel2.png|280px|center]]
-
Reduction:
+
 
 +
 
 +
'''Reduction:'''
[[File:SwitchModel3.png|460px|center]]
[[File:SwitchModel3.png|460px|center]]
-
Solution:
+
 
 +
 
 +
'''Solution:'''
[[File:SwitchModel3-1.png|882px|center]]
[[File:SwitchModel3-1.png|882px|center]]
-
Equilibration:
+
 
 +
 
 +
'''Equilibration:'''
[[File:SwitchModel4-1.png|313px|center]]
[[File:SwitchModel4-1.png|313px|center]]
-
Stabilization:
+
 
 +
 
 +
'''Stabilization:'''
[[File:SwitchModel4-2.png|607px|center]]
[[File:SwitchModel4-2.png|607px|center]]
-
Calibration:
+
 
 +
 
 +
'''Calibration:'''
[[File:SwitchModel5-1.png|640px|center]]
[[File:SwitchModel5-1.png|640px|center]]
[[File:SwitchModel5-2.png|752px|center]]
[[File:SwitchModel5-2.png|752px|center]]
-
Optimization:
+
 
 +
 
 +
'''Optimization:'''
[[File:SwitchModel6.png|245px|center]]
[[File:SwitchModel6.png|245px|center]]
-
SimBiology:
+
 
 +
 
 +
'''SimBiology:'''
[[File:Switch.jpg|center]]
[[File:Switch.jpg|center]]
-
[https://docs.google.com/document/d/1LewVc-RZQhwpdMyIHFYT9f4rbFrQXvppZOoQjKCZlLM/edit?usp=sharing Switch Model Derivation]
+
 
 +
'''Derivation:''' [https://docs.google.com/document/d/1LewVc-RZQhwpdMyIHFYT9f4rbFrQXvppZOoQjKCZlLM/edit?usp=sharing Switch Model]
=Inducible Hbif Model=
=Inducible Hbif Model=
 +
 +
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.
 +
[[File:InducibleHbifModelFigure.png|center]]
[[File:InducibleHbifModelFigure.png|center]]
-
Equation:
+
 
 +
 
 +
'''Equation:'''
[[File:InducibleHbifModel1.png |355px|center]]
[[File:InducibleHbifModel1.png |355px|center]]
[[File:InducibleHbifModel2.png |242px|center]]
[[File:InducibleHbifModel2.png |242px|center]]
-
Reduction:
+
 
 +
 
 +
'''Reduction:'''
[[File:InducibleHbifModel3.png |460px|center]]
[[File:InducibleHbifModel3.png |460px|center]]
-
Solution:
+
 
 +
 
 +
'''Solution:'''
[[File:InducibleHbifModel4.png |853px|center]]
[[File:InducibleHbifModel4.png |853px|center]]
-
Equilibration:
+
 
 +
 
 +
'''Equilibration:'''
[[File:InducibleHbifModel5.png |204px|center]]
[[File:InducibleHbifModel5.png |204px|center]]
-
Stabilization:
+
 
 +
 
 +
'''Stabilization:'''
[[File:InducibleHbifModel5-2.png |815px|center]]
[[File:InducibleHbifModel5-2.png |815px|center]]
-
Calibration:
+
 
 +
 
 +
'''Calibration:'''
[[File:InducibleHbifModel6-1.png |335px|center]]
[[File:InducibleHbifModel6-1.png |335px|center]]
[[File:InducibleHbifModel6-2.png |395px|center]]
[[File:InducibleHbifModel6-2.png |395px|center]]
[[File:InducibleHbifModel6-3.png |238px|center]]
[[File:InducibleHbifModel6-3.png |238px|center]]
-
Optimization of Hbif:
 
-
[[File:InducibleHbifModel7-1.png |117px|center]]
 
-
Optimization of ON:
 
-
[[File:InducibleHbifModel7-2.png |347px|center]]
 
-
SimBiology:
 
-
[[File:InducedHbif.jpg|center]]
 
-
[https://docs.google.com/document/d/1PNIkP84ksvYOu-HDYZghVVDkhDlMX7Iz9T6sFZnmfr8/edit?usp=sharing Inducible Hbif Model Derivation]
 
-
=Future Directions=
+
'''Optimization of Hbif:'''
-
Data!
+
[[File:InducibleHbifModel7-1.png |117px|center]]
-
<!-- *** 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>
+
'''Optimization of ON:'''
-
<div id="box" style="width: 700px; margin-left: 137px; padding: 5px; border: 3px solid #000; background-color: #fe2b33;">
+
[[File:InducibleHbifModel7-2.png |347px|center]]
-
<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 *** -->
 
 +
'''SimBiology:'''
 +
[[File:InducedHbif.jpg|center]]
-
{| style="color:#1b2c8a;background-color:#0c6;" cellpadding="3" cellspacing="1" border="1" bordercolor="#fff" width="62%" align="center"
+
'''Derivation:''' [https://docs.google.com/document/d/1PNIkP84ksvYOu-HDYZghVVDkhDlMX7Iz9T6sFZnmfr8/edit?usp=sharing Inducible Hbif Model]
-
!align="center"|[[Team:Michigan|Home]]
+
-
!align="center"|[[Team:Michigan/Team|Team]]
+
-
!align="center"|[https://igem.org/Team.cgi?year=2013&team_name=Michigan Official Team Profile]
+
-
!align="center"|[[Team:Michigan/Project|Project]]
+
-
!align="center"|[[Team:Michigan/Parts|Parts Submitted to the Registry]]
+
-
!align="center"|[[Team:Michigan/Modeling|Modeling]]
+
-
!align="center"|[[Team:Michigan/Notebook|Notebook]]
+
-
!align="center"|[[Team:Michigan/Safety|Safety]]
+
-
!align="center"|[[Team:Michigan/Attributions|Attributions]]
+
-
|}
+

Latest revision as of 23:58, 18 October 2013

Xayona Website Template

Contents

Introduction

Switch modeling

The goal for our modeling was to create rigorous standard procedures for mathematically representing and analyzing systems. Using our project as an example, we worked to demonstrate the viability of our approach by calculating the switch positions and optimum parameters. This standardization should be useful in combination with system specific methods to create better modeling.

Mass action modeling

Mass action is a form of modeling that uses systems of homogenous differential and difference equations. Each equation is a sum of rates and each rate is proportional to the species involved. The advantage is its straightforward uniform formulation, its minimal assumptions, and its flexibility of complexity. Mass action itself uses very few assumptions but it can be used to represent systems of varying degrees of simplicity. The fundamentals of chemical kinetics already champion the use of this modeling in its use of rate equations. There chemical systems are completely defined by the products, reactants, and reaction rate constants. This allows chemical modelling to be highly standardized and hence more accessible. It's actually the uniform nature of mass action that's responsible. Most importantly for iGEM, all mass action models can actually be put in a standard form as follows:

Mass action.png

• x is a species

• 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 the parameters of the system, analytical (symbolic) modeling precisely predicts a system's behavior. Furthermore, unlike numerical methods, analytical modeling accurately predicts the data needed to determine parameters and avoids the inherent numerical rounding and approximation errors of more crude analyses. Analytical modeling naturally lends itself to standardization because the results are unique to any given model. Correct solutions, equilibria, eigenvalues, etc. will always have the same expression regardless of who or what calculated them. The same can't be said for numerical methods. With this and the iGEM community in mind, an analytical standard was developed that will allow for the exchange and automation of future modeling.

ERSESCO Algorithm

The ERSESCO algorithm the Michigan team developed provides a clearly defined method for analyzing mathematical 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. To get stability, find the condition for which all eigenvalues have a negative real part.

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

SimBiology is a useful MATLAB package provided to all iGEM teams. It can be used as a tool to quickly run numerical approximations of biological systems whether or not data for system parameters are known. SimBiology can quickly produce a model from a simple schematic diagram. As a result, the package can quickly provide a qualitative analysis of a system with little work from the user.

In the following models, our team applied our ERSESCO algorithm to fully analyze the systems at hand. Finally a SimBiology simulation was used as qualitative verification. We present the results from each step of the ERSESCO algorithm, followed by a SimBiology schematic with a graph showing the general behavior of the system, followed by a link to a more detailed derivation of the models.


Expression Model

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.


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


Derivation: Expression Model

Switch Model

The Switch Model describes how the switch equilibrates when the concentrations of the recombinases FimE and Hbif approach constant values. The 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.


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


Derivation: Switch Model

Inducible Hbif Model

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.


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

Derivation: Inducible Hbif Model