At the moment I am at the European Conference on Mathematical and Theoretical Biology in Krakow. This is a joint conference with the Society for Mathematical Biology and it is very large, with more than 900 participants. This leads to a huge number of parallel sessions and the need to choose very carefully in order to get the most profit from the conference. On one day, for instance, there were two cases with two sessions on immunology occurring simultaneously.

On Tuesday I went to a session on biochemical reaction networks. This included a talk by Gheorghe Craciun with a large expository component which I found enlightening. He raised the question of when a system of ODE with polynomial coefficients can be interpreted as coming from a system of chemical reactions with mass action kinetics. He mentioned a theorem about this and after asking him for details I was able to find a corresponding paper by Hars and Toth. This is in the Colloquia Mathematica Societatis Janos Bolyai, which is a priori not easily accessible. The paper is, however, available as a PDF file on the web page of Janos Toth. A chemical reaction network gives rise to a system of equations of the form where the and are polynomials with positive coefficients. They represent the contributions from reactions where the species with concentration is on the right and left side respectively. The result of Hars and Toth is that any system of this algebraic form can be obtained from a reaction network. It was pointed out by Craciun in his talk that this means that arbitrarily complicated dynamics can be incorporated into systems coming from reaction networks. If we have a system of the form we can replace it by . This changes the system but does not change the orbits of solutions. If, for instance, we start with the Lorenz system with unknowns , and we can simply translate the coordinates so as to move the interesting dynamics into the region where all coordinates are positive and then multiply the result by . This preserves the strange attractor structure. This result may be compared with the Perelson-Wallwork theorem discussed in a previous post.

The construction of the reaction network reproducing the given equations is not very complicated. The main problem is keeping track of the notation. Suppose we start with a system of equations in variables which is polynomial and satisfies the necessary condition already mentioned. The reaction network can be constructed in the following way. (It is not at all unique.) Introduce one species for each . The right hand side of each equation is a sum of terms of the form and one reaction is introduced for each of these terms. To explain what it is suppose without generality that it belongs to the first equation. If hen the reaction transforms the complex to the complex with rate constant . The only species where there is a net production is and so this reaction only contributes to the first equation. Moreover it does so with the desired term. On the other hand if the reaction transforms the complex to the complex with rate constant . The assumption on the system assures that .

July 26, 2011 at 6:05 am |

I liked Gheoge’s idea first, then tried to do some calculations on the Lorenz model, and I think with this premultiplication we often get an equation the solutions of which blow up. Could you find appropriate constants to produce chaotic behaviour in the frst orthant? E. g.

NDSolve[{

x’[t] == x[t] y[t] z[t] (-3 (x[t] – y[t])),

y’[t] == x[t] y[t] z[t] (-x[t] z[t] + 26.5 x[t] – y[t]),

z’[t] == x[t] y[t] z[t] (x[t] y[t] – z[t]),

x[0] == 1, z[0] == 2, y[0] == 5}, {x, y, z}, {t, 0, 200},

MaxSteps -> \[Infinity]];

does not work.

Thanks for the post, János

July 26, 2011 at 11:25 am |

Thanks for the comment. I see now that things are not as easy as I thought. It would be nice if this idea for treating the Lorenz system could be made watertight.

Alan