Mathematical Models

Mathematical modeling is an important component in the construction of rationally designed gene networks. The complexity of biological systems makes necessary the use of sophisticated mathematical models to accurately predict network functionality. Models are useful in tuning the individual components of a network to increase system robustness. They also provide a basis for comparison of experimental results, allowing researchers to test the validity of their assumptions.

When dealing with gene activation and repression networks, such as those described in the common [CellularMemory:Biological Designs |Biological Designs] for synthetic cellular memory, one of three mathematical models is often used to model the transcription rate of a gene(s).

The Hill Equation

Three closely related mathematical models can be used to describe binding of two given molecules in a biological system. . All of them can be understood by looking at the definition of the dissociation constant of a given molecular complex. Kd = [A][B]/[AB]. Because [AB] = [Btotal] - [B] you can show that the percentage of unbound B molecules [B]/[Btotal] = 1/(1+X/Kd). If you think of B as a RNA Polymerase binding site and A as a repressor, multiply by the maximal rate of transcription (beta) yields the promoter activity as a function of repressor concentration. (See discussion of Kd below)
The value of K is closely related to the dissociation constant of the inducer/repressor complex. The dissociation constant measures the propensity for this complex to fall apart. Kd = [I][R]/[IR]. This gives the dissociation constant units of concentration. It can be shown mathematically that Kd is equal to [I] at which half of the repressor is bound to the DNA and half is not ([R] = [IR]). Because the amount of repressor bound to the DNA is proportional to the level of transcription, K (being closely related to Kd) is the value at the inducer yields half of the maximum level of transcription.
Michaelis-Menton model doesn't take cooperativity into account.
The Hill equation accounts for cooperativity of binding, but is not a completely accurate description at very low concentrations of repressor/activator.
The Monod, Changeux, and Wymann model is the most accurate model of these three, but is not commonly used because of its increased complexity and the fact that the hill equation accurately describes most situations.