Team:Freiburg/Modeling
From 2012.igem.org
Line 17: | Line 17: | ||
- | + | [[Image:SCHEMA2.png|center|720px|link=]] | |
Revision as of 15:35, 26 October 2012
Companel|DNA calculating competitive reactions
![PartssymbolT.png](/wiki/images/a/a2/PartssymbolT.png)
![FreigemCompanelDNAlogo.png](/wiki/images/9/99/FreigemCompanelDNAlogo.png)
![Freigemquote.png](/wiki/images/2/21/Freigemquote.png)
The result of these efforts is CompanelDNA, an excel application that allows for predicting the equilibrium concentrations of competitive annealing. It can of course be used for every kind of chemical equilibrium where two substances compete for reacting with the same third reactant.
![SCHEMA2.png](/wiki/images/0/05/SCHEMA2.png)
Theoretical background
Although it seems to be a simple scenario, we did not find any preexisting solution to this problem - which came to us as quite a surprise. So we started by looking at the law of mass action behind the annealing reaction. Hereafter we will use the following abbreviations and symbols:
![Freigem1T.png](/wiki/images/1/17/Freigem1T.png)
The laws of mass action for the two independent reactions are:
![Freigem2T.png](/wiki/images/c/ca/Freigem2T.png)
If we want to know what is going to happen to the concentrations we originally put in our reaction/PCR tube, we can express the term in relation to the overall concentrations:
![Freigem3T.png](/wiki/images/1/18/Freigem3T.png)
Bearing in mind that the only unknown variables are (AB) and (CB) (K can be calculated, see below), one would try to solve these two equations with two variables. A number of other equations can furthermore be generated. However, these equations had no exact solutions when we used Mathematica to have them solved. All the proposed solutions are not applicable to a wider combination of parameters. Even when we asked for mathematical advice (thanks to Paul Staab, Munich) we got no further. So we decided to solve graphically instead (thanks to Konrad Schieban, Zurich). In order to do so the equations need to be solved for one variabel:
![Freigem4T.png](/wiki/images/4/43/Freigem4T.png)
These equations can now separately be calculated by varying the variable (AB). When a value of (AB) is found that yields the same value of (CB) in both equations, both conditions are fulfilled. Graphically this would be represented by the intersection of the graphs of both equations. Basically this is what our application does. As a first result it clearly showed that the temperature specificity is not strict enough to be used as a elongation determinant, as we had originally hoped.
How does Companel|DNA work?
When using Companel you need to "bring" your thermodynamical parameters dH° and dS° for each independent annealing reaction. Due to limited resources could no calculator for these numbers be included. However, there are several free calculators online. More information on how to obtain dH° and dS° can be found in the application.
Once this sequence specific information has been entered into the application, the constants K for a series of temperatures are calculated by using the Gibbs function. Subsequently the melting plots for each individual reaction are displayed and give you a first impression of the annealing behavior of your primers.
![FreigemCompanel2.png](/wiki/images/2/24/FreigemCompanel2.png)
Finally you enter concentrations, temperature range and expected product yield. The resulting concentrations of annealing product for each primer are shown graphically and in a list.
![FreigemCompanel1.png](/wiki/images/4/4a/FreigemCompanel1.png)
You just want to have one specific situation at a given temperature calculated? Just use the second part of the application Companel|chem and enter barely the concentrations as well as the two independent equilibrium constants there.
![Constantsonly.png](/wiki/images/f/fd/Constantsonly.png)
What is Companel|DNA good for?
Besides the actual annealing of primers, which it has been intended for, Companel possesses an intrinsic universality when it comes to quantifying competitive processes. For instance a scientist could easily assess the binding priority of two or more (the system can of course be expanded) transcription factors or protein-protein interactions. He or she would just need to measure the two or more separate reactions, which could readily be achieved in an in vitro essay. One could then predict the interaction in vivo - at least thermodynamically. Numerous discoveries during the last years, such as the impact of the circadian clock on metabolism and drug efficacy, have shown that profound research and sometimes even successful treatment can no longer be based on the mere observation of qualitative effects alone. In the challenge to bring both biology and medicine closer to quantitative, predictable and hence more understandable science our application could play a helpful role.
please send us an email to freigem2012(at)googlemail.com