It is shown, how to construct a system of ordinary differential equations of arbitrary order, which has the periodic attractor and models some genetic network of arbitrary size. The construction is carried out by combining of multiple systems of lower dimensions with known periodic attractors. In our example the six-dimensional system is constructed, using two identical three-dimensional systems, which have stable periodic solutions.
Genetic engineering started 39 years ago at the Miami Winter Symposium, where three independent groups described Agrobacterium tumefaciens-mediated genetic transformation, leading to the production of normal, fertile transgenic plants .
Genetic engineering is a modern area of biotechnology that combines knowledge and techniques from a whole block of related sciences. It is the process of changing an organism’s DNA to introduce new, desirable traits. Many organisms, from bacteria to plants and animals, have been genetically modified for academic, medical, agricultural, and industrial purposes.
Genetic regulatory network (GRN) is responsible for many important functions of living organisms. It exists in any in any cell of any living organism. Huge amounts of experimental data are stored and need systematization and research. One the effective way of treating this is to consider the mathematical models of GRN. If they are appropriate and exhibit sufficient correspondence with experimental information, they (models) can be used for understanding, managing and control over GRN. We focus on models formulated in terms of dynamical systems. Dynamical systems, used for modeling GRN, consist of ordinary differential equations of special structure. Look at the system below.
Consider the general form of writing the -dimensional dynamical system, that is expected to model a genetic regulatory network:
where , and are parameters, and are elements of the regulatory matrix . The parameters of the GRN have the biological interpretations [2,3].
In the system Eq. (1) the sigmoid logistic function is used .
It is to be mentioned, that first system of this form was used to model populations of neurons , . There were attempts to use systems of the form Eq. (1) for the creation of Virtual Network Topology to manage traffic in telecommunication networks . It is recommended also to be acquainted with the articles ,  on the subject.
In this article, we would like to consider the construction of models of GRN with the required properties. Due to the restricted volume of the conference proceedings, we will discuss the creation of models of any dimensionality, considering bounded periodic attractor. Solutions, tending to periodic attractors, reflect rhythmical processes in organisms. Our plan is the following. First, we consider the three-dimensional system of ordinary differential equations of the form Eq. (1) and show, that it has a stable periodic trajectory, which serves as an attractor for other trajectories. Then we construct a greater system, which consists of three-dimensional systems. The periodic attractors in any of three-dimensional systems, when combined in higher dimensional phase space, produce a periodic attractor of a higher dimension.
2. Two-dimensional systems
The two-dimensional system of the form Eq. (1) is defined, if all the parameters are given. Consider 1; 20, 40; 1, –0.3. The coefficients are elements of the matrix and behavior of solutions depends on it.
Consider the system Eq. (1) with regulatory matrix:
The nullclines of the system Eq. (1) together with the stable periodic solution are depicted in Fig. 1 .
Fig. 1The nullclines x1 – red, x2 – green
There is a unique critical point at (0.5; 0.5). The characteristic equation for equilibria (0.5; 0.5) is:
where –7, 62. Solving the equation we have 3.5±7.05. The type of equilibria is an unstable focus.
3. Three-dimensional systems
The three-dimensional system of the form Eq. (1) is defined, if all the parameters are given. Consider 1; 4; 4, 0, 0. The coefficients are elements of regulatory matrix and behavior of solutions depends on it.
Consider the system Eq. (1) with regulatory matrix:
The nullclines of system Eq. (1) together with the stable periodic solution are depicted in Fig. 2. There is a unique equilibria at (0.5; 0.5; 0.5). The central location of the critical point is due to the special choice of the parameters .
Fig. 2The nullclines x1 – red, x2 – green, x3 – blue and the periodic trajectory
The characteristic equation for the critical point (0.5; 0.5; 0.5) is:
where –1, 3, –13. Solving the equation we have 1.102±1.6866; –3.20. The type of critical point is an unstable saddle-focus.
In fact, we have provided an example of the 3D system without stable critical points. More examples of this kind were obtained in the paper .
4. Six-dimensional systems
We are able now to construct systems of any dimensions, combining the above two- and/or three-dimensional systems. Let us construct six-dimensional system with a periodic attractor, combining two identical three-dimensional systems. Consider 6 by 6 block matrix , where the 3 by 3 matrix Eq. (4) is placed at the left-upper and right-lower corners. Other spaces are filled by zeros. The resulting system of the form Eq. (1) describes a GRN with 6 elements and it has the six-dimensional periodic attractor. We can visualize three-dimensional projections of this attractor on subspaces.
The initial conditions for a solution of the six-dimensional system, tending to a six-dimensional periodic attractor, are:
Fig. 3a) The projections of 6D trajectories to 3D subspace(x1,x2,x3) (it is solution of the base three-dimensional system), b) The projections of 6D trajectories to 3D subspace (x1,x2,x4)
Fig. 4a) The projections of 6D trajectories to 3D subspace (x1,x4,x6), b) The projections of 6D trajectories to 3D subspace (x2,x3,x6)
Fig. 5a) The projections of 6D trajectories to 3D subspace(x1,x5,x6), b) The projections of 6D trajectories to 3D subspace (x3,x4,x5)
Periodic attractors can be constructed for higher-dimensional systems of the form Eq. (1), using known periodic trajectories for low-dimensional systems. In our example, we first constructed the three-dimensional system with a regulatory matrix Eq. (4). Then we constructed the six-dimensional square matrix with two blocks on the main diagonal, consisting of the matrix Eq. (4). Other parameters are chosen as for the three-dimensional case. Then we computed the solution of the six-dimensional system with the above initial conditions. This trajectory, projected on , tends to the three-dimensional projection of the attractor. This trajectory, projected on , tends to the three-dimensional projection of the attractor onto the subspace . The initial conditions for these projections, differ. This is why Fig. 4(a) and Fig. 5(a) are not identical.
It is possible, proceeding in this way, to construct attractors for any dimension. In the second section, the two-dimensional system is provided, which has a stable periodic solution. A periodic attractor, therefore, can be constructed for any dimension, combining two and three dimensional blocks in the main regulatory matrix (for example, 4 = 2 + 2, 5 = 3 + 2, 6 = 3 + 3 = 2 + 2 + 2, etc.). In particular, the six-dimensional attractor can be constructed by combining three independent two-dimensional systems.
The situation is different if we combine the six-dimensional system with two three-dimensional systems, which have stable periodic solutions, but with different periods. Recall, that in our example periods for solutions of both three-dimensional components are the same since the three-dimensional systems are identical. The attractor in a resulting six-dimensional system must exist, but it should be a chaotic subset of the phase space, formed by non-periodic six-dimensional trajectories. It is like Lissajous figures, where the periods are incommensurable.
The authors have not disclosed any funding.
I. K. Vasil, “A history of plant biotechnology: from the cell theory of Schleiden and Schwann to biotech crops,” Plant Cell Reports, Vol. 27, No. 9, pp. 1423–1440, Sep. 2008, https://doi.org/10.1007/s00299-008-0571-4
N. Vijesh, S. K. Chakrabarti, and J. Sreekumar, “Modeling of gene regulatory networks: A review,” Journal of Biomedical Science and Engineering, Vol. 6, No. 2, pp. 223–231, 2013, https://doi.org/10.4236/jbise.2013.62a027
I. Samuilik and F. Sadyrbaev, “Modelling three dimensional gene regulatory networks,” WSEAS Transactions on Systems and Control, Vol. 16, pp. 755–763, 2021.
I. Samuilik, “Genetic engineering-construction of a network of four dimensions with a chaotic attractor,” Vibroengineering Procedia, Vol. 44, pp. 66–70, 2022.
H. R. Wilson and J. D. Cowan, “Excitatory and inhibitory interactions in localized populations of model neurons,” Biophysical Journal, Vol. 12, No. 1, pp. 1–24, Jan. 1972, https://doi.org/10.1016/s0006-3495(72)86068-5
V. W. Noonburg, Differential Equations: From Calculus to Dynamical Systems. Providence, Rhode Island: MAA Press, 2019.
Y. Koizumi, T. Miyamura, S.I. Arakawa, E. Oki, K. Shiomoto, and M. Murata, “Adaptive virtual network topology control based on attractor selection,” Journal of Lightwave Technology, Vol. 28, No. 11, pp. 1720–1731, Jun. 2010, https://doi.org/10.1109/jlt.2010.2048412
C. Furusawa and K. Kaneko, “A generic mechanism for adaptive growth rate regulation,” PLoS Computational Biology, Vol. 4, No. 1, p. e3, Jan. 2008, https://doi.org/10.1371/journal.pcbi.0040003
H. de Jong, “Modeling and Simulation of Genetic Regulatory Systems: A Literature Review,” Journal of Computational Biology, Vol. 9, No. 1, pp. 67–103, Jan. 2002, https://doi.org/10.1089/10665270252833208
I. Samuilik, F. Sadyrbaev, and S. Atslega, “Mathematical modeling of nonlinear dynamic systems,” Engineering for Rural Development, Vol. 21, pp. 172–178, 2022.
O. Kozlovska and F. Sadyrbaev, “Models of genetic networks with given properties,” WSEAS Transactions on Computer Research, Vol. 10, pp. 43–49, Apr. 2022, https://doi.org/10.37394/232018.2022.10.6
“Wolfram Mathematica.” https://www.wolfram.com/mathematica/.