Gene regulation: Dynamical system for proteins
Gene expression with negative feedback
So far we have not really introduced feedback on protein formation; therefore we need a network of genes. The simplest model for negative feedback is represented by the following scheme [cf. Figure 8-76A in Alberts et al. (2015, 6th Edition, p. 515)].
Due to a input signal, the expression of the gene A to the protein A starts. This protein formed is an activating transcription factor for the gene R, which in turn forms a protein R that acts as inhibitory regulator protein for the transcription of gene A.
By using the previous formulas for gene expression, we come to the following set of differential equations for the proteins A and R:
By scaling we can reduce the system of differential equations to the following system of two quantities and that are associated with the scaled concentrations of A and R:
Let us for convenience choose . Then we get:
For an equilibrium we should have
It now turns out that solutions always go in the course of time toward this equilibrium; it is a so-called attracting equilibrium. You see this best in a so-called phase portrait. This is a diagram in which one often draws a number of solution curves at one time and draws the vector field corresponding to the 2-dimensional system of differential equations. At each point in the phase portrait we can namely calculate a vector and draw the vector in the plane. In our example, the phase portrait with one red solution curve that starts in looks as follows: