Next Article in Journal
Improved Square-Root Cubature Kalman Filtering Algorithm for Nonlinear Systems with Dual Unknown Inputs
Next Article in Special Issue
Bifurcation of Limit Cycles from a Focus-Parabolic-Type Critical Point in Piecewise Smooth Cubic Systems
Previous Article in Journal
Pointwise Estimation of Anisotropic Regression Functions Using Wavelets with Data-Driven Selection Rule
Previous Article in Special Issue
Reduction in Degrees of Freedom for Large-Scale Nonlinear Systems in Computer-Assisted Proofs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New 3D Chaotic Attractor in Gene Regulatory Network

by
Olga Kozlovska
1,†,
Felix Sadyrbaev
2,† and
Inna Samuilik
1,*,†
1
Institute of Applied Mathematics, Riga Technical University, Zunda Embankment 10, LV-1048 Riga, Latvia
2
Institute of Mathematics and Computer Science, University of Latvia, Rainis Boulevard 29, LV-1459 Riga, Latvia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2024, 12(1), 100; https://doi.org/10.3390/math12010100
Submission received: 7 November 2023 / Revised: 8 December 2023 / Accepted: 12 December 2023 / Published: 27 December 2023
(This article belongs to the Special Issue Advances in Nonlinear Dynamics and Chaos: Theory and Application)

Abstract

:
This paper introduces a new 3D chaotic attractor in a gene regulatory network. The proposed model has eighteen parameters. Formulas for characteristic numbers of critical points for three-dimensional systems were considered. We show that the three equilibrium points of the new chaotic 3D system are unstable and deduce that the three-dimensional system exhibits chaotic behavior. The possible outcomes of this 3D model were compared with the results of the Chua circuit. The bifurcation structures of the proposed 3D system are investigated numerically, showing periodic solutions and chaotic solutions. Lyapunov exponents and Kaplan-Yorke dimension are calculated. For calculations, the Wolfram Mathematica is used.
MSC:
34C25; 34D45; 92-08; 92-10

1. Introduction

Deterministic chaos of nonlinear dynamic systems is not chaos, usually understood as complete disorganization and randomness of events; modern ideas about chaos are to some extent closer to the original ancient Greek ones: “chaos” like a boundless disordered mass from which everything that exists arose.
Chaos theory is a branch of science that deals with complex systems and their behavior. Meteorologist Edward Lorenz made a major contribution to chaos theory [1]. During the 1960s, Lorenz developed a computer program that simulated the motion of air masses within the Earth’s atmosphere [2]. The Lorenz model demonstrated a susceptibility to initial conditions [3]. Tiny variations in input data resulted in significant discrepancies in outcomes over time. This reliance on the starting conditions was termed chaos. Lorenz coined a famous feature of chaos called the “butterfly effect”, the notion that a butterfly stirring the air in Hong Kong today can transform storm systems in New York next month [4]. After Lorenz’s studies, many chaotic systems have been presented to the literature Rossler system [5], Chen system [6], LU system [7], Sprott system [8] and others. For a long time, chaos theory was considered a kind of mathematical abstraction that had no confirmation in real conditions. Now it has applications in various scientific disciplines, including physics [9], biology (in the study of uneven heart rate and an uneven number of diseases) [10], meteorology, economics [11], finance [12], geology [13], computer science, engineering, algorithmic trading, politics [14], population dynamics [15], robotics [16], philosophy [17] and mathematics.
Chaos is a complex phenomenon defying easy classification or identification. Although there exists no universally accepted definition for chaos, the solutions of chaotic systems typically exhibit the characteristics [18].
Characteristics of chaos:
  • chaotic behavior is characterized by the presence of an attractor to which all nearby solutions tend to converge over time, provided there is ample time for the process [19];
  • a common feature of chaotic solutions lies in the geometric structure of their attractors. These attractors often exhibit intricate and unconventional shapes, characterized by a twisted and “strange” appearance. This strangeness is indicative of a fractional (fractal) dimension, although it’s important to note that this isn’t always the case [19];
  • sensitivity to initial conditions [18].
Nonlinear ordinary differential equations are the most widespread formalism for modeling genetic regulatory networks [20]. The main contributions of the present study are summarized as follows:
  • formulas for characteristic numbers of critical points for three-dimensional systems were considered;
  • the new chaotic attractor is obtained;
  • the three-dimensional system (7) can have attractors of various kinds;
  • the irregular behavior of solutions near the chaotic attractor is conceivable and may manifest within a narrow parameter range.

2. Materials and Methods

2.1. Chua Circuit

Leon O. Chua is a renowned electrical engineer and computer scientist who has made significant contributions to the fields of electronics, circuits, and nonlinear science. He was born on 28 June 1936, in Manila, Philippines. He is a professor in the Electrical Engineering and Computer Sciences department at the University of California, Berkeley. The professor’s scientific projects are related to cellular neural network technologies, bifurcation theory, nonlinear dynamics, and chaos theory. Leon Chua aimed to show that it was possible to create chaos. For this purpose, in 1983 he assembled an electrical circuit [21]. It is a simple oscillatory circuit demonstrating a series of bifurcations and a transition to chaos. In the simple case, Chua’s equations can be written in the following form:
d x d τ = α ( y x f ( x ) ) , d y d τ = x y + z , d z d τ = β y ,
f ( x ) = b x + a b , x 1 , a x , | x | 1 , b x a + b , x 1 ,
where x , y , z -phase variables, α , β , a , b -parameters ( α > 0 , β > 0 , a < 0 , b < 0 ) . Consider
a = 8 7 ; b = 5 7 , α = 9 , β = 14 .
We will replace the piecewise linear function with a cubic polynomial [22].
f ( x ) = 1 16 x 3 7 6 x .
The initial conditions are
x ( 0 ) = 0.1 ; y ( 0 ) = 0.2 ; z ( 0 ) = 1 .
The nullclines are depicted in Figure 1.
The system possesses precisely three critical points, determined and illustrated as the intersection points of the nullclines in Figure 1. There are three critical points at ( 1.633 , 0 , 1.633 ) , ( 0 , 0 , 0 ) and ( 1.633 , 0 , 1.633 ) . Linearization around these points provides us with the characteristic numbers λ . The characteristic numbers are considered in Table 1.
The type of two critical points is a saddle-focus with one-dimensional stable and two-dimensional unstable manifolds. The type of critical point ( 0 , 0 , 0 ) is a saddle-focus with two-dimensional stable and one-dimensional unstable manifolds. The chaotic attractor is depicted in Figure 2 and the graph of solutions is depicted in Figure 3.
In the equilibrium state at the origin, one root is real and positive, indicating the direction in which the initial disturbance will grow. A pair of complex conjugate roots with a negative real part signifies the presence of rotational motion and twisting of the trajectory towards a singular critical point in the plane of rotation. For each of the other two equilibrium states, the motion is unstable in the plane of rotation (the real parts of the complex conjugate pair are greater than zero). At the same time, there exists a direction in which the phase trajectory approaches the critical point.
Lyapunov exponents after 5000 steps were L 1 = 0.3329 , L 2 = 0.0012 , and L 3 = 2.9923 . The largest Lyapunov exponent was greater than 0 as shown in Figure 4.
The Lyapunov exponents characterize the exponential expansion or contraction of phase-space entities, such as one-dimensional lengths, two-dimensional areas, and three-dimensional volumes. The largest Lyapunov exponent denoted as L E 1 characterizes the average rate of divergence over time between two adjacent trajectories separated by a specific distance δ . The sum of the first n Lyapunov exponents describes the divergence or convergence rate of an n-dimensional phase-space volume [23]. The paper [24] states that the presence of at least one positive Lyapunov exponent leads to the divergence of neighboring trajectories, classifying the phase-space motion as “chaotic”, but the positivity of the calculated senior Lyapunov’s exponent is not the criteria of occurrence in the system of chaotic dynamics [25,26]. The computation of the full Lyapunov exponent spectrum is a rather mathematically complicated issue. Computations are performed using Wolfram Mathematica. For Lyapunov exponents calculation the package “lce.m for Mathematica” was used [27].
Meanwhile, according to the Kaplan–Yorke formula,
D K Y = j + 1 | L j + 1 | i = 1 j L i = 2 + L 1 + L 2 | L 3 | = 2.1109
the fractal dimension also further verifies that the new system has a chaotic behavior [28]. For a dissipative dynamical system, the sum of all Lyapunov exponents should equate to a negative value [29]. Only dissipative dynamical systems have attractors [30].

2.2. 3D Gene Regulatory System

Consider the system
d x d t = 1 1 + e μ 1 ( w 11 x + w 12 y + w 13 z θ 1 ) v 1 x , d y d t = 1 1 + e μ 2 ( w 21 x + w 22 y + w 23 z θ 2 ) v 2 y , d z d t = 1 1 + e μ 3 ( w 31 x + w 32 y + w 33 z θ 3 ) v 3 z ,
where μ i , θ i and v i are the parameters, w i j are the coefficients of the so-called regulatory matrix
W = w 11 w 12 w 13 w 21 w 22 w 23 w 31 w 32 w 33 .
The parameters of the GRN have the following biological interpretations:
  • v i —degradation of the i-th gene expression product;
  • w i j —the connection weight or strength of control of gene j on gene i. Positive values of w i j signify activating influences, whereas negative values denote repressing influences;
  • θ i —The impact of external stimuli on gene i is reflected in its ability to modulate the gene’s responsiveness to activating or repressing factors [31].
The nullclines and the critical points for the system are defined by the relations
x = 1 v 1 1 1 + e μ 1 ( w 11 x + w 12 y + w 13 z θ 1 ) , y = 1 v 2 1 1 + e μ 2 ( w 21 x w 22 y + w 23 z θ 2 ) , z = 1 v 3 1 1 + e μ 2 ( w 21 x + w 22 y + w 33 z θ 3 ) .
The sigmoidal function f ( t ) = 1 1 + e μ is used in (7). Sigmoidal functions exhibit a continuous increase from zero to one and possess a solitary inflection point. While various sigmoidal functions exist, the one mentioned above is particularly suitable for analysis and visualization [32]. Such systems were considered in [20,33,34,35].

2.3. Linearized System

The linearized system for any critical point ( x * , y * , z * )
u 1 = v 1 u 1 + μ 1 w 11 g 1 u 1 + μ 1 w 12 g 1 u 2 + μ 1 w 13 g 1 u 3 , u 2 = v 2 u 2 + μ 2 w 21 g 2 u 1 + μ 2 w 22 g 2 u 2 + μ 2 w 23 g 2 u 3 , u 3 = v 3 u 3 + μ 3 w 31 g 3 u 1 + μ 3 w 32 g 3 u 2 + μ 3 w 33 g 3 u 3 ,
where
g 1 = e μ 1 ( w 11 x * + w 12 y * + w 13 z * θ 1 ) [ 1 + e μ 1 ( w 11 x * + w 12 y * + w 13 z * θ 1 ) ] 2 ,
g 2 = e μ 2 ( w 21 x * + w 22 y * + w 23 z * θ 2 ) [ 1 + e μ 2 ( w 21 x * + w 22 y * + w 23 z * θ 2 ) ] 2 ,
g 3 = e μ 3 ( w 31 x * + w 32 y * + w 33 z * θ 3 ) [ 1 + e μ 3 ( w 31 x * + w 32 y * + w 33 z * θ 3 ) ] 2 .
One has
A λ I = μ 1 w 11 g 1 v 1 λ μ 1 w 12 g 1 μ 1 w 13 g 1 μ 2 w 21 g 2 μ 2 w 22 g 2 v 2 λ μ 2 w 23 g 2 μ 3 w 31 g 3 μ 3 w 32 g 3 μ 3 w 33 g 3 v 3 λ
and the characteristic equation is
det | A λ I | = λ 3 + λ 2 ( v 1 v 2 v 3 + μ 1 w 11 g 1 + μ 2 w 22 g 2 + μ 3 w 33 g 3 ) + λ ( g 1 v 3 μ 1 w 11 +
+ μ 2 w 22 g 2 v 3 + g 1 g 2 w 21 μ 1 μ 2 w 12 g 1 g 2 w 11 w 22 μ 1 μ 2 + g 1 g 3 w 31 w 13 μ 1 μ 3
g 1 g 3 w 11 w 33 μ 1 μ 3 + g 2 g 3 w 32 w 23 μ 2 μ 3 g 2 g 3 w 22 w 33 μ 2 μ 3 v 1 ( v 2 + v 3 g 2 w 22 μ 2 g 3 w 33 μ 3 ) +
+ v 2 ( v 3 + g 1 w 11 μ 1 + g 3 w 33 μ 3 ) ) + v 1 ( v 2 ( v 3 + g 3 w 33 μ 3 ) + g 2 μ 2 ( v 3 w 22 + g 3 w 32 w 23 μ 3
g 3 w 22 w 33 μ 3 ) ) + g 1 μ 3 ( v 2 ( v 3 w 11 + g 3 ( w 31 w 13 w 11 w 33 ) μ 3 ) + g 2 μ 2 ( v 3 ( w 21 w 12 w 11 w 22 ) +
+ g 3 ( w 31 w 22 w 13 + w 21 w 32 w 13 + w 31 w 12 w 23 w 11 w 32 w 23 w 21 w 12 w 33 + w 11 w 22 w 33 ) μ 3 ) ) = 0 .
The characteristic equation can be rewritten as
λ 3 + A λ 2 + B λ + C = 0 ,
where
A = ( v 1 + v 2 + v 3 ) + g 1 w 11 μ 1 + g 2 w 22 μ 2 + g 3 w 33 μ 3 ,
B = μ 1 μ 2 w 31 w 13 g 1 g 3 μ 2 μ 3 w 32 w 23 g 2 g 3 + μ 1 μ 2 w 21 w 12 g 1 g 2 ( μ 2 w 22 g 2 v 2 ) ( μ 3 w 33 g 3 v 3 ) ( μ 1 w 11 g 1 v 1 ) ( μ 3 w 33 g 3 v 3 ) ( μ 1 w 11 g 1 v 1 ) ( μ 2 w 22 g 2 v 2 ) ,
C = ( μ 1 w 11 g 1 v 1 ) ( μ 2 w 22 g 2 v 2 ) ( μ 3 w 33 g 3 v 3 ) + μ 1 μ 2 μ 3 w 21 w 32 w 23 g 1 g 2 g 3 + μ 1 μ 2 μ 3 w 31 w 12 w 23 g 1 g 2 g 3 μ 1 μ 3 w 31 w 13 g 1 g 3 ( μ 2 w 22 g 2 v 2 ) μ 2 μ 3 w 32 w 23 g 2 g 3 ( μ 1 w 11 g 1 v 1 ) μ 1 μ 2 w 21 w 12 g 1 g 2 ( μ 3 w 33 g 3 v 3 ) .
Theorem 1.
The vector field ( f 1 ( x , y , z ) , f 2 ( x , y , z ) , f 3 ( x , y , z ) ) , where f 1 , f 2 and f 3 are the right sides of the equations in (7), is directed inward on the boundary of the domain Q 3 : = { ( x , y , z ) : 0 < x < 1 v 1 , 0 < y < 1 v 2 , 0 < z < 1 v 3 } .
Proof of Theorem 1.
Take one of faces of the parallelepiped Q 3 , for example, x = 0 . The vector field there in the x direction is f 1 v 1 , x = f > 0 . Take face x = 1 v 1 . The vector field in the x direction is f 1 v 1 , x = f 1 v 1 , 1 v 1 = f 1 1 < 0 . In both cases, the vector field along the x axis is directed inside Q 3 . Similarly, other faces of Q 3 can be considered. □
Theorem 2.
System (7) has at least one equilibrium (critical point). All equilibria are located in the open box Q 3 : = { ( x , y , z ) : 0 < x < 1 v 1 , 0 < y < 1 v 2 , 0 < z < 1 v 3 } .
This follows from the result of the mapping of a topological ball into itself. The second assertion follows from the fact that nullclines meet and can intersect only in Q 3 .

3. Results

3.1. 3D Chaotic Attractor

The system (7) with the matrix
W = 0.06 0.825 0.19 1.915 1.76 0.32 3.7 0.128 0.985
and μ 1 = 4 , μ 2 = 5.02 , μ 3 = 8 , v 1 = v 2 = v 3 = 1 and
θ 1 = w 11 + w 12 + w 13 2 ,
θ 2 = w 21 + w 22 + w 23 2 ,
θ 3 = w 31 + w 32 + w 33 2
has three critical points. Linearization around these points provides us with the characteristic numbers λ .
The initial conditions are
x ( 0 ) = 0.592 ; y ( 0 ) = 0.85 ; z ( 0 ) = 2 .
The nullclines are depicted in Figure 5.
The characteristic equation for critical points is (12). The coefficients of characteristic Equation (12) are considered in Table 2.
Characteristic numbers are considered in Table 3.
The type of two critical points is a saddle-focus with one-dimensional stable and two-dimensional unstable manifolds. The type of critical point ( 0.5 , 0.5 , 0.5 ) is a saddle-focus with two-dimensional stable and one-dimensional unstable manifolds. The chaotic attractor is depicted in Figure 6 and Figure 7. The graph of solutions is depicted in Figure 8 and Figure 9.
Lyapunov exponents after 5000 steps were L 1 = 0.0398 , L 2 = 0.0026 and L 3 = 0.7665 . The largest Lyapunov exponent was greater than 0 as shown in Figure 10.

3.2. Bifurcation

The bifurcation theory of differential equation systems, first introduced in the works of Jules Henri Poincare (1854–1912), elucidates qualitative, abrupt alterations in the phase portraits of these systems as their parameters undergo continuous and smooth changes [36]. We take w 23 as a bifurcation parameter (that is, the third element in the second row) in the regulatory matrix (13). Computations are performed using Wolfram Mathematica. Results of calculations for the system (7) with regulatory matrix (13), changing the parameter w 23 . are considered in Table 4.
Lyapunov exponents are considered in Table 5.
Periodic attractors are depicted in Figure 11, Figure 12 and Figure 13. Graphs of solutions are depicted in Figure 14, Figure 15 and Figure 16. The chaotic attractor is depicted in Figure 17 and the graph of solutions is depicted in Figure 18.

4. Discussion

A definite similarity was found between the Chua circuits and systems of the form (7), which are used in the mathematical modeling of genetic and neuronal networks. An analysis of critical points and their mutual positions in the Chua system motivated us to study similar configurations in GRN systems. The local analysis of critical points gave hints of what should be arranged in GRN systems to get similar behavior of trajectories as in the Chua system. For this, multiple parameters, built into GRN systems, were useful. The desired behavior can be reached by shifting the nullclines, changing θ parameters, or/and changing the elements of the regulatory matrix, and regulating the form of sigmoidal functions by changing μ parameters. Looking at the data being collected during the study, one may observe some properties that eventually lead to chaotic behavior. As an example, Table 1, Table 2 and Table 3 with the local characteristics of the critical points, should be mentioned. The bifurcation analysis made in the article concerns the varying of one parameter. The changes in the characteristics of critical points and fractal dimensionality of the chaotic attractor, give rise to conclusions and may indicate directions of further numerical experiments. Some questions arise. For instance, the monotone change in the parameter w 23 leads to generally non-monotone changes in characteristics of the critical points and KY fractal dimensionalities of the chaotic attractor. Analyzing the collected data may lead to the formulation of directions for further studies. The questions that should be answered: what is the minimal (or optimal) number of critical points needed for the birth of a strange attractor; what is the minimal configuration (positions, local characteristics) that ensures a transition to chaotic behavior; is the transition to chaotic behavior in the studied system possible without passing through the stable periodic solution stage; could sufficient conditions be formulated to ensure passage to chaotic behavior from a given configuration; what is the description of this configuration; what are the necessary conditions concerning passage to chaotic behavior; is it possible to pass from periodic attractors to chaotic ones through continuous change of KY dimensionality; generally, what is the role of jump changes in a system to go to the strange (chaotic) attractor?
This list can be continued and answering some, or all, of the formulated questions concerning the system (7) (with matrix (13), at least) would indicate the essential progress in understanding the process of formation of a chaotic attractor. So further investigations in this direction promise results, which may be relevant and important for the general theory of dynamical systems.

Author Contributions

Conceptualization, F.S.; methodology, O.K. and I.S.; software, O.K. and I.S.; formal analysis, O.K. and F.S.; investigation, O.K.; resources, I.S.; data curation, I.S. and F.S.; writing—original draft preparation, I.S. and O.K.; writing—review and editing, I.S. and F.S.; visualization, O.K.; supervision, I.S.; project administration, I.S.; funding acquisition, O.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data were created or analyzed in this study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
GRNGene regulatory network

References

  1. Ramar, R. Design of a New Chaotic System with Sine Function: Dynamic Analysis and Offset Boosting Control. Chaos Theory Appl. 2023, 5, 118–126. [Google Scholar] [CrossRef]
  2. Mansour, M.; Donmez, T.B.; Kutlu, M.C.; Freeman, C. Respiratory Diseases Prediction from a Novel Chaotic System. Chaos Theory Appl. 2023, 5, 20–26. [Google Scholar] [CrossRef]
  3. Keles, Z.; Sonugur, G.; Alcn, M. The Modeling of the Rucklidge Chaotic System with Artificial Neural Networks. Chaos Theory Appl. 2023, 5, 59–64. [Google Scholar] [CrossRef]
  4. Jorgensen, S.E. Encyclopedia of Ecology; Elsevier: Amsterdam, The Netherlands, 2008. [Google Scholar]
  5. Ibrahim, K.M.; Jamal, R.K.; Ali, F.H. Chaotic behaviour of the Rossler model and its analysis by using bifurcations of limit cycles and chaotic attractors. J. Phys. Conf. Ser. 2018, 1003, 012099. [Google Scholar] [CrossRef]
  6. Xiyin, L.; Guoyuan, Q. Mechanical analysis of Chen chaotic system. Chaos Solitons Fractals 2017, 98, 173–177. [Google Scholar] [CrossRef]
  7. Lu, J.; Chen, G.; Guoyuan, Q. A new chaotic attractor coined. Int. J. Bifurc. Chaos 2002, 12, 659–661. [Google Scholar] [CrossRef]
  8. Mota, M.C.; Oliveira., R.D.S.; Qi, G. Dynamic aspects of Sprott BC chaotic system. Discret. Contin. Syst. 2021, 26, 1653–1673. [Google Scholar] [CrossRef]
  9. Kazuyuki, A. Chaos and Its Applications. Procedia IUTAM 2012, 5, 199–203. [Google Scholar] [CrossRef]
  10. Ott, E. Chaos in Dynamical Systems, 2nd ed.; Cambridge University Press: London, UK, 2002. [Google Scholar] [CrossRef]
  11. Feier, S. The Chaos Theory and its Application. J. Phys. Conf. Ser. 2012, 012118, 2021. [Google Scholar]
  12. Klioutchnikov, I.; Sigova, M.; Beizerov, N. Chaos Theory in Finance. Procedia Comput. Sci. 2017, 2017, 368–375. [Google Scholar] [CrossRef]
  13. Biswas, H.R.; Hasan, M.M.; Bala, S.K. Chaos theory and its applications in our real life. Barishal Univ. J. 2018, 5, 123–140. [Google Scholar]
  14. Joan, P.; Dandoy, R. Chaos theory and its application in political science. In Proceedings of the IPSA World Congress, Fukuoka, Japan, 9–13 July 2006. [Google Scholar]
  15. May, M.R. Chaos and the dynamics of biological populations. Nucl. Proc. Suppl. 1987, 2, 225–245. [Google Scholar] [CrossRef]
  16. Zang, X.; Iqbal, S.; Zhu, Y.; Liu, X.; Zhao, J. Applications of Chaotic Dynamics in Robotics. Int. J. Adv. Robot. 2016, 13. [Google Scholar] [CrossRef]
  17. Rana, S.M.S. Bifurcation Analysis and 0–1 Chaos Test of a Discrete T System. Chaos Theory Appl. 2023, 5, 90–104. [Google Scholar] [CrossRef]
  18. Lynch, S.E. Dynamical Systems with Applications Using Mathematica; Springer: Cham, Switzerland, 2017. [Google Scholar]
  19. Hastings, A.; Hom, C.L.; Ellner, S.; Turchin, P.; Godfray, H.C.J. Chaos in Ecology: Is Mother Nature a Strange Attractor? Annu. Rev. Ecol. Syst. 1993, 24, 1–33. [Google Scholar] [CrossRef]
  20. Samuilik, I.; Sadyrbaev, F.; Ogorelova, D. Comparative Analysis of Models of Gene and Neural Networks. Contemp. Math. 2023, 4, 217–229. [Google Scholar] [CrossRef]
  21. Bugajevskij, M.; Ponomarenko, V. Study of Chua Circuit Behavior; Fundamental and Comparative Research; Saratov State University: Saratov, Russia, 1999. [Google Scholar]
  22. Pavlov, A.N. Study of Chua Generator Dynamics Modes; Saratov State University: Saratov, Russia, 2006; pp. 1–13. [Google Scholar]
  23. Syed Ali, M.; Vadivel, R. Decentralized Event-Triggered Exponential Stability for Uncertain Delayed Genetic Regulatory Networks with Markov Jump Parameters and Distributed Delays. Neural Process. Lett. 2018, 47, 1219–1252. [Google Scholar] [CrossRef]
  24. Hoover, W.G.; Tull, C.G.; Posch, H. Negative Lyapunov exponents for dissipative systems. Phys. Lett. A 1998, 131, 211–215. [Google Scholar] [CrossRef]
  25. Magnitskii, N.A. Universal Bifurcation Chaos Theory and Its New Applications. Mathematics 2023, 11, 2536. [Google Scholar] [CrossRef]
  26. Magnitskii, N.A. Bifurcation Theory of Dynamical Chaos. In Chaos Theory; IntechOpen: Rijeka, Croatia, 2018. [Google Scholar]
  27. Software. Available online: www.msandri.it/soft.html (accessed on 21 September 2023).
  28. Wang, J.; Dong, C.; Li, H. A New Variable-Boostable 3D Chaotic System with Hidden and Coexisting Attractors: Dynamical Analysis, Periodic Orbit Coding, Circuit Simulation, and Synchronization. Fractal Fract. 2022, 6, 740. [Google Scholar] [CrossRef]
  29. Sayed, W.S.; Radwan, A.G.; Fahmy, H.A.A. Chaos and Bifurcation in Controllable Jerk-Based chaotic Attractors. In Nonlinear Dynamical Systems with chaotic and Hidden Attractors; Springer: Cham, Switzerland, 2018; pp. 45–70. [Google Scholar]
  30. Nikolov, S.; Nedkova, N. Gyrostat Model Regular And Chaotic Behavior. J. Theor. Appl. Mech. 2015, 45, 15–30. [Google Scholar] [CrossRef]
  31. Vijesh, N.; Kumar, S.; Sreekumar, C.J. Modelling three dimensional gene regulatory networks. WSEAS Trans. Syst. Control 2021, 16, 755–763. [Google Scholar] [CrossRef]
  32. Ogorelova, D.; Sadyrbaev, F.; Samuilik, I. On Targeted Control over Trajectories of Dynamical Systems Arising in Models of Complex Networks. Mathematics 2023, 11, 2206. [Google Scholar] [CrossRef]
  33. Samuilik, I.; Sadyrbaev, F. On trajectories of a system modeling evolution of genetic networks. Math. Biosci. Eng. 2023, 20, 2232–2242. [Google Scholar] [CrossRef]
  34. Vijesh, N.; Chakrabarti, S.; Sreekumar, J. Modeling of gene regulatory networks: A review. J. Biomed. Sci. Eng. 2013, 6, 223–231. [Google Scholar] [CrossRef]
  35. Kozlovska, O.; Sadyrbaev, F. On attractors in systems of ordinary differential equations arising in models of genetic networks. Vibroeng. Procedia 2023, 49, 136–140. [Google Scholar] [CrossRef]
  36. Magnitskii, N.A.; Sidorov, S.V. New Methods for Chaotic Dynamics; World Scientific: Singapore, 2006; pp. 1–384. [Google Scholar] [CrossRef]
Figure 1. The visualization of nullclines for the system (1). ( x red, y green, z blue).
Figure 1. The visualization of nullclines for the system (1). ( x red, y green, z blue).
Mathematics 12 00100 g001
Figure 2. Chua’s double-scroll attractor: Phase portrait for the system (1).
Figure 2. Chua’s double-scroll attractor: Phase portrait for the system (1).
Mathematics 12 00100 g002
Figure 3. Solutions ( x , y , z ) of the system (1) with the initial conditions x ( 0 ) = 0.1 ;   y ( 0 ) = 0.2 ;   z ( 0 ) = 1 .
Figure 3. Solutions ( x , y , z ) of the system (1) with the initial conditions x ( 0 ) = 0.1 ;   y ( 0 ) = 0.2 ;   z ( 0 ) = 1 .
Mathematics 12 00100 g003
Figure 4. Lyapunov exponent spectrum of Chua circuit (1).
Figure 4. Lyapunov exponent spectrum of Chua circuit (1).
Mathematics 12 00100 g004
Figure 5. Nullclines x 1 —red, x 2 —green, x 3 —blue of the system (7) with the regulatory matrix (13).
Figure 5. Nullclines x 1 —red, x 2 —green, x 3 —blue of the system (7) with the regulatory matrix (13).
Mathematics 12 00100 g005
Figure 6. The visualization of chaotic attractor of the system (7) with the regulatory matrix (13).
Figure 6. The visualization of chaotic attractor of the system (7) with the regulatory matrix (13).
Mathematics 12 00100 g006
Figure 7. The visualization of chaotic attractor of the system (7) with the regulatory matrix (13), other view.
Figure 7. The visualization of chaotic attractor of the system (7) with the regulatory matrix (13), other view.
Mathematics 12 00100 g007
Figure 8. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13).
Figure 8. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13).
Mathematics 12 00100 g008
Figure 9. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13).
Figure 9. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13).
Mathematics 12 00100 g009
Figure 10. Lyapunov exponent of the system (7) with the regulatory matrix (13).
Figure 10. Lyapunov exponent of the system (7) with the regulatory matrix (13).
Mathematics 12 00100 g010
Figure 11. The visualization of the periodic attractor, w 23 = 0.569 .
Figure 11. The visualization of the periodic attractor, w 23 = 0.569 .
Mathematics 12 00100 g011
Figure 12. The visualization of the periodic attractor, w 23 = 0.4 .
Figure 12. The visualization of the periodic attractor, w 23 = 0.4 .
Mathematics 12 00100 g012
Figure 13. The visualization of the periodic attractor, w 23 = 8 .
Figure 13. The visualization of the periodic attractor, w 23 = 8 .
Mathematics 12 00100 g013
Figure 14. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13), w 23 = 0.569 .
Figure 14. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13), w 23 = 0.569 .
Mathematics 12 00100 g014
Figure 15. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13), w 23 = 0.4 .
Figure 15. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13), w 23 = 0.4 .
Mathematics 12 00100 g015
Figure 16. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13), w 23 = 8 .
Figure 16. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13), w 23 = 8 .
Mathematics 12 00100 g016
Figure 17. The visualization of the chaotic attractor, w 23 = 0.2401 .
Figure 17. The visualization of the chaotic attractor, w 23 = 0.2401 .
Mathematics 12 00100 g017
Figure 18. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13), w 23 = 0.2401 .
Figure 18. Solutions ( x , y , z ) of the system (7) with the regulatory matrix (13), w 23 = 0.2401 .
Mathematics 12 00100 g018
Table 1. Characteristic numbers.
Table 1. Characteristic numbers.
The Critical Point λ 1 λ 2 λ 3
( 1.633 , 0 , 1.633 ) −4.35630.1782 − 3.1315 i0.1782 + 3.1315 i
( 0 , 0 , 0 ) 2.4730−0.9865 − 2.7734 i−0.9865 − 2.7734 i
( 1.633 , 0 , 1.633 ) −4.35630.1782 − 3.1315 i0.1782 + 3.1315 i
Table 2. Coefficients of characteristic Equation (12).
Table 2. Coefficients of characteristic Equation (12).
The Critical PointABC
( 0.4398 , 0.4545 , 0.0032 ) −0.725248−0.557956−0.787449
( 0.5 , 0.5 , 0.5 ) 1.23880.1959831.59353
( 0.5602 , 0.5455 , 0.9968 ) −0.725248−0.557956−0.787449
Table 3. Characteristic numbers.
Table 3. Characteristic numbers.
The Critical Point λ 1 λ 2 λ 3
( 0.4398 , 0.4545 , 0.0032 ) −0.97800.12638 − 0.888361 i0.12638 + 0.888361 i
( 0.5 , 0.5 , 0.5 ) 1.82477−0.292985 − 0.887376 i−0.292985 + 0.887376 i
( 0.5602 , 0.5455 , 0.9968 ) −0.97800.12638 − 0.888361 i0.12638 + 0.888361 i
Table 4. Results of calculations for the system (7) with regulatory matrix (13), changing the parameter w 23 .
Table 4. Results of calculations for the system (7) with regulatory matrix (13), changing the parameter w 23 .
w 23 x * , y * , z * λ 1 λ 2 λ 3
−0.569(0.2673, 0.6735, 0.00002)−0.9998−0.0051 − 0.6838 i−0.0051 − 0.6838 i
(0.5, 0.5, 0.5)2.0836−0.4224 − 1.2107 i−0.4224 + 1.2107 i
(0.7327, 0.3265, 0.99998)−0.9998−0.0051 − 0.6838 i−0.0051 + 0.6838 i
−0.5(0.3190, 0.6015, 0.0001)−0.99900.0849 − 0.7646 i0.0849 + 0.7646 i
(0.5, 0.5, 0.5)−0.3913−0.3913 − 1.1361−0.3913 + 1.1361
(0.6810, 0.3985, 0.9999)−0.99900.0849 − 0.7646 i0.0849 + 0.7646 i
−0.4(0.3885, 0.5144, 0.0007)−0.99410.1319 − 0.8505 i0.1319 + 0.8505 i
(0.5, 0.5, 0.5)1.9197−0.3404 − 1.0104 i−0.3404 + 1.0104 i
(0.6115, 0.4856, 0.9993)−0.99410.1319 − 0.8505 i0.1319 + 0.8505 i
−0.3(0.4521, 0.4406, 0.0046)−0.96980.1214 − 0.8932 i0.1214 + 0.8932 i
(0.5, 0.5, 0.5)1.7985−0.2799 − 0.8520 i−0.2799 + 0.8520 i
(0.5479, 0.5594, 0.9954)−0.96980.1214 − 0.8932 i0.1214 + 0.8932 i
−0.241(0.4873, 0.4023, 0.0133)−0.92470.1061 − 0.8937 i0.1061 + 0.8937 i
(0.5, 0.5, 0.5)1.7130−0.2371 − 0.7312 i−0.2371 + 0.7312 i
(0.5127, 0.5977, 0.9868)−0.92470.1061 − 0.8937 i0.1061 + 0.8937 i
−0.24(0.4879, 0.4017, 0.0135)−0.92350.1059 − 0.8935 i0.1059 + 0.8935 i
(0.5, 0.5, 0.5)1.7114−0.2363 − 0.7289 i−0.2363 + 0.7289 i
(0.5121, 0.5983, 0.9865)−0.92350.1059 − 0.8935 i0.1059 + 0.8935 i
0(0.5, 0.5, 0.5)−0.70710.9729 − 0.3612 i0.9729 + 0.3612 i
8(0.5, 0.5, 0.5)−3.75502.4969 − 3.1062 i2.4969 + 3.1062 i
Table 5. Lyapunov exponents.
Table 5. Lyapunov exponents.
w 23 LE 1 LE 2 LE 3
0.569 −0.0106−0.0149−0.9997
0.5 −0.00002−0.1529−0.9907
0.4 0.0026−0.2235−0.8975
0.3 0.0309−0.0127−0.7389
0.241 0.0512−0.0062−0.6901
0.24 0.0020−0.0260−0.7109
0−0.0029−0.6131−0.9340
80.0016−0.9139−1.5327
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kozlovska, O.; Sadyrbaev, F.; Samuilik, I. A New 3D Chaotic Attractor in Gene Regulatory Network. Mathematics 2024, 12, 100. https://doi.org/10.3390/math12010100

AMA Style

Kozlovska O, Sadyrbaev F, Samuilik I. A New 3D Chaotic Attractor in Gene Regulatory Network. Mathematics. 2024; 12(1):100. https://doi.org/10.3390/math12010100

Chicago/Turabian Style

Kozlovska, Olga, Felix Sadyrbaev, and Inna Samuilik. 2024. "A New 3D Chaotic Attractor in Gene Regulatory Network" Mathematics 12, no. 1: 100. https://doi.org/10.3390/math12010100

APA Style

Kozlovska, O., Sadyrbaev, F., & Samuilik, I. (2024). A New 3D Chaotic Attractor in Gene Regulatory Network. Mathematics, 12(1), 100. https://doi.org/10.3390/math12010100

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop