Team:Freiburg/Modeling

From 2012.igem.org




Companel|DNA calculating competitive reactions



PartssymbolT.png



FreigemCompanelDNAlogo.png
During the course of our project we encountered a question that we had never thought about before. We came up with different strategies to synthesize a TAL protein. One of the ideas was to take advantage of the thermodynamic behavior of different DNA strands and synthesize a TAL proteins by means of different primers that anneal at different temperatures to the template strand. We thought it to be possible to elongate a sequence by using this mechanism - simply changing the annealing temperature slightly whilst running a usual PCR. We assumed that the difference in optimal annealing temperatures of different primer sequences would favor one replication product over the others.
Freigemquote.png
The only thing we did not know was if this was significant enough to be used for the assembly of at least ten TAL domains. After unproductive research on this issue we decided to try modeling this scenario ourselves.


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.



Schema4.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

The laws of mass action for the two independent reactions are:

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

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

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


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


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



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.




In case the file does not work properly,
please send us an email to freigem2012(at)googlemail.com





Back to top