Complexity
Volume 2017 (2017), Article ID 1409865, 20 pages
https://doi.org/10.1155/2017/1409865

Research Article

Bifurcations and Dynamics of the Rb-E2F Pathway Involving miR449

Lingling Li¹ and Jianwei Shen²

¹Department of Basic Courses, Guangdong Polytechnic College, Zhaoqing 526100, China
²Institute of Applied Mathematics, Xuchang University, Xuchang, Henan 461000, China

Correspondence should be addressed to Jianwei Shen; xcjwshen@gmail.com

Received 23 February 2017; Revised 19 April 2017; Accepted 14 May 2017; Published 19 October 2017

Academic Editor: Shanmugam Lakshmanan

Copyright © 2017 Lingling Li and Jianwei Shen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We focused on the gene regulative network involving Rb-E2F pathway and microRNAs (miR449) and studied the influence of time delay on the dynamical behaviors of Rb-E2F pathway by using Hopf bifurcation theory. It is shown that under certain assumptions the steady state of the delay model is asymptotically stable for all delay values; there is a critical value under another set of conditions; the steady state is stable when the time delay is less than the critical value, while the steady state is changed to be unstable when the time delay is greater than the critical value. Thus, Hopf bifurcation appears at the steady state when the delay passes through the critical value. Numerical simulations were presented to illustrate the theoretical results.

1. Introduction

In the last few years, an increasing number of noncoding RNA (or microRNA) have been discovered to play central regulatory roles in gene regulation processes of prokaryotes and eukaryotes [1]. It has become evident that miRNAs regulate a variety of biological processes, and their expression is often deregulated in human malignancy. On one hand, miRNAs play roles in tumorigenesis by modulating oncogenic and tumor suppressor pathways. On the other hand, the expressions of miRNAs can be regulated by several oncogenic or tumor suppressor transcription factors [2]. In this paper, we focus on miR449, which can induce cell senescence and apoptosis and act as a tumor suppressor through regulating Rb/E2F activity [2, 3]. In recent years, large numbers of researches have focused on the mechanisms controlling cellular proliferation associated with human cancer regulated by Rb-E2F pathway experimentally [4–7]. Rb and E2F proteins play important roles in the regulation of cell division, cell growth, and programmed cell death by controlling the expressions of genes involved in these processes, which are best known for their regulation of the cell cycle at the G1/S transition [8]. As the first identified tumor suppressor gene [9], Rb is recognized to play a fundamental role in a signaling pathway that controls cell proliferation [10]. Rb regulates the transcription of genes that are essential for DNA replication and cell cycle progression by binding and inhibiting E2F transcription factors [11]. In the Rb-E2F pathway including negative feedback loops involving miR449, miR449 provides a twofold safety mechanism to avoid excessive E2F-induced proliferation by cell cycle arrest and apoptosis [12]. Mathematical models have been established to explain the nonlinear dynamical behaviors of the Rb-E2F pathway [12–14], which mainly concentrate on the stability and bifurcation of the deterministic systems, but not taking into account the effects of time delay.

Time delay is one of the most important characteristics of gene regulation. In many cases, a gene regulates the expression of another gene by its products (RNAs or proteins). Since it takes time to generate those products and different processes need different amounts of time, time delayed regulation is ubiquitous in cellular processes [15]. Time-delay system is also called system with aftereffect or hereditary system [16, 17], and it is common in mathematical biology, such as population dynamics, the chemostat, neural network, blood cell maturation, transcriptional regulator dynamics, virus dynamics, and genetic network [18–23]. Recent papers have demonstrated that complex dynamical behaviors can arise as a consequence of time delays in biological systems. For example, in gene regulatory networks, time delays may lead to oscillations in protein levels and existing oscillations may become more robust [24, 25]. Oscillatory cellular dynamics, in particular periodic oscillation, plays an important role in maintaining homeostasis of living organisms. In addition, taking into account time delays in models of gene regulatory networks is often essential to capture the whole range of dynamic behaviors. For example, in experiments, a single self-repressed gene has been observed to display oscillatory behavior, but which cannot be deduced by models that ignore time delays. However, this oscillatory behavior is reproduced by a mathematical model including time delays. In addition, theoretical analysis that ignored time delays led to the erroneous conclusion that oscillations were not possible for this single gene [18, 26]. Many kinds of methods have been introduced to infer time delayed gene regulatory network [27–31]. Here, we mainly use the local linearization approach to analyze the nonlinear system.

In this paper, firstly, we will introduce Rb-E2F pathway network involving miR449 modeled by time delayed differential equations. Secondly, we will study the dynamical behaviors of the model and derive sufficient conditions of the oscillation by using Hopf bifurcation theory. Particularly, we will also prove that there are periodic solutions under certain conditions. Finally, numerical simulations will be showed to illustrate the theoretical results.

2. A Simple Gene Regulatory Network Mediated by miR449

2.1. A Simple Mathematical Model of Gene Regulation with a Delayed Negative Feedback Loop

A simplified model (Figure 1) of the Rb-E2F pathway begins with growth signals activating CycD. Initially, E2F is bound to and repressed by Rb, a tumor suppressor protein that is dysfunctional in several major cancers. CycD represses the repressor Rb and allows E2F to be turned on. Myc also induces E2F transcription. Subsequently, E2F activates the transcription of CycE, which forms a complex with Cdk2 to further remove Rb repression, establishing a positive feedback loop. E2F also activates its own transcription, constituting another positive feedback loop. An interesting addition to the regulatory mechanism of Rb-E2F network is the recent discovery that miR449 modulates E2F activity. It has been demonstrated that E2F strongly upregulates the expression of miR449. In turn, E2F is inhibited by miR449 through regulating different transcripts. On one hand, miR449 directly affects level of its target transcript Myc and therefore lowers E2F concentration. On the other hand, miR449 directly affects E2F inducer Cdk6 and CycE, thus forming negative feedback loops [12, 14]. Yao et al. [14] provided a mathematical model in the absence of miR449 and indicated that the Rb-E2F pathway acts as a bistable switch to convert signal inputs into all-or-none E2F responses. Yan et al. [12] gave another mathematical model and further investigated the stabilities and bifurcations of E2F, CycE, and miR449 in the participation of miR449. Our main work in this paper is considering the effects of time delays on the dynamic behaviors of the model including miR449. When time delays are taken into account, the time delayed differential equation model of the network including miR449 is described by the following system:

represent the concentrations of E2F, miR449, Myc, CycD-Cdk4/6 complex, CycE-Cdk2 complex, Rb, phosphorylated Rb, and Rb/E2F complex, respectively. And is intensity of growth factor. In the following simulations, all the values of parameters are shown in Table 1 unless specified elsewhere.

Table 1: Parameters for the model.

Rate constant
Value

0.4

0.15

0.15

0.003

0.15

18

0.92

18

0.92

0.25

180

1.4

0.15

0.02

0.6

0.15

1

0.92

0.7

0.92

1

0.5

0.7

0.6

0.15

0.45

0.5

0.03

0.15

1.5

1

0.92

0.35

0.15

1.5

0.7

0.92

0.18

3.6

0.01

180

18

0.92

18

0.92

0.06

18

0.92

18

0.92

18

0.92

18

0.92

3.6

0.01

0.06

180

18

0.92

18

0.92

0.03

Note. The descriptions of the parameters in Table 1 are shown in [13, 15].

Figure 1: Rb-E2F pathway mediated by miR449.

2.2. Oscillation Induced by Time Delay

In this subsection, we consider system (1) with as state variables. The linearized system of (1) at equilibrium point is as follows:

where

Then we can obtain the characteristic equation of (2) at the equilibrium as follows:

where is the identity matrix, and the characteristic equation (5) has the following form:

where the values of are showed in the Appendix.

If we assume that , we will have

where .

(1) If , (7) becomes

where

According to the Routh-Hurwitz criterion, all roots of (8) have negative real parts if and only if all the subdeterminants in the diagonal are positive; that is,

(2) If , considering the transcendental equation (7), clearly is a root of (7) if and only if

Separating the real and imaginary parts of (11), we have

Adding up the squares of both equations of (12), we have

where

Let ; (13) becomes

Denote

Lemma 1. If , (15) has at least one positive root.

Proof. Clearly, , and . Hence, there exists , so that . This completes the proof.

Lemma 2. If , the sufficient condition for (15) has positive roots being and .

Proof. From (16), we have ; suppose the equation has seven real roots and satisfies , and is the local minimum value, if and ; there exists , so that ; this completes the proof.

Lemma 3. If , according to the Routh-Hurwitz criterion, all the roots of (15) have negative real parts if and only if the subdeterminants in the diagonal are positive, that is, , .

Suppose that (15) has positive roots; without loss of generality, we assume that it has eight positive roots, denoted by , respectively. Hence, (13) has eight positive roots, say .

From (12), we can get

where

Define

Let be the root of (7) satisfying .

Lemma 4. Consider the exponential polynomial , where and is polynomial about . As vary, the sum of the orders of the zeros of on the open right half plane can change only if a zero appears on or cross the imaginary axis.
Then, we have the following theoretical results.

Theorem 5. Suppose that conditions (10) are satisfied.

(i)If , then all roots of (7) have negative real parts for all ; thus the steady state of system (2) is absolutely stable.(ii)If or , and , then all the roots of (7) have negative real parts when ; thus the steady state of system (2) is asymptotically stable.(iii)If the condition of (ii) is satisfied, , and , then is a pair of simple purely imaginary roots of (15) and all other roots have negative real parts. Moreover, . Thus, system (2) exhibits the Hopf bifurcation at .

3. Numerical Analysis

In this section, we demonstrate the above theoretical results by numerical method. When we take , and the other parameters are shown in Table 1, system (1) becomes

The system has a positive equilibrium point , . Using Theorem 5, there is a critical value of the time delay . The equilibrium point is stable when (see Figures 2 and 3); the equilibrium point becomes unstable and a Hopf bifurcation occurs when passes through the critical value (see Figure 4). The bifurcation diagrams of system (20) are shown in Figure 5, where the control parameter is the time delay .

Figure 2: An asymptotically stable equilibrium for . Assume initial conditions are , and .

Figure 3: An asymptotically stable equilibrium for . Assume initial conditions are , and .

Figure 4: A periodic solution bifurcated from equilibrium for . Assume initial conditions are , and .

Figure 5: Bifurcation induced by time delay with . Assume initial conditions are , and .

When we take , and the other parameters are shown in Table 1, system (1) becomes

The system has a positive equilibrium point , . Using Theorem 5, there is a critical value of the time delay . The equilibrium point is stable when (see Figures 6 and 7); the equilibrium point becomes unstable and a Hopf bifurcation occurs when passes through the critical value (see Figure 8). The bifurcation diagrams of system (21) are shown in Figure 9, where the control parameter is the time delay .

Figure 6: An asymptotically stable equilibrium for . Assume initial conditions are , and .

Figure 7: An asymptotically stable equilibrium for . Assume initial conditions are , and .

Figure 8: A periodic solution bifurcated from equilibrium for . Assume initial conditions are , and .

Figure 9: Bifurcation induced by time delay with . Assume initial conditions are , and .

4. Discussion

In this paper, we studied the dynamical behaviors of the Rb-E2F pathway including negative feedback loops involving miR449 by using Hopf bifurcation theory. On one hand, we gave the detailed theoretical analysis of system (1) in which there are some dynamical behaviors, and it is shown that, under certain conditions, the equilibrium point is asymptotically stable for all the delay ; however, if these conditions are not met, there is a critical delay ; when , the equilibrium point is asymptotically stable; when , periodic oscillations occur. Thus, Hopf bifurcation appears at the steady state when the delay passes through the critical value . On the other hand, through numerical simulations, we calculated the critical value of , giving the time courses of E2F, CycE, and miR449 both and , drawing the bifurcation diagrams of E2F, CycE, and miR449, respectively, and testifying the correctness of the theory.

Appendix

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This work is supported by National Natural Science Foundation of China (11772291, 11572278, and 11572084) and Innovative Research Team in University of Henan Province (13IRTSTHN019).

References

J. Shen, Z. Liu, W. Zheng, F. Xu, and L. Chen, “Oscillatory dynamics in a simple gene regulatory network mediated by small RNAs,” Physica A. Statistical Mechanics and its Applications, vol. 388, no. 14, pp. 2995–3000, 2009.

X. Yang, M. Feng, X. Jiang et al., “miR-449a and miR-449b are direct transcriptional targets of E2F1 and negatively regulate pRb-E2F1 activity through a feedback loop by targeting CDK6 and CDC25A,” Genes & Development, vol. 23, no. 20, pp. 2388–2393, 2009.

T. Bou Kheir, E. Futoma-Kazmierczak, A. Jacobsen et al., “miR-449 inhibits cell proliferation and is down-regulated in gastric cancer,” Molecular Cancer, vol. 10, article 29, 2011.

E. Khav, Visualing an Rb-E2F cellular switch that controls cell oliferation [Dissertation, thesis], Academic Dissertation, The University of Arizona, 2013.

J. White, E. Stead, R. Faast, S. Conn, P. Cartwright, and S. Dalton, “Developmental activation of the Rb-E2F pathway and establishment of cell cycle-regulated cyclin-dependent kinase activity during embryonic stem cell differentiation,” Molecular Biology of the Cell, vol. 16, no. 4, pp. 2018–2027, 2005.

J. R. Nevins, G. Leone, J. DeGregori, and L. Jakoi, “Role of the Rb/E2F pathway in cell growth control,” Journal of Cellular Physiology, vol. 173, no. 2, pp. 233–236, 1997.

J. R. Nevins, “The Rb/E2F pathway and cancer,” Human Molecular Genetics, vol. 10, no. 7, pp. 699–703, 2001.

L. Hangnoh, Regulation of differentiation-specific genes by the Drosophila RB, E2F, and Myb-interacting proteins complex (dREAM [Dissertation, thesis], Academic Dissertation, The State University of New Jersey, 2011.

D. Hanahan and R. A. Weinberg, “The hallmarks of cancer,” Cell, vol. 100, no. 1, pp. 57–70, 2000.

R. C. Sears and J. R. Nevins, “Signaling networks that link cell proliferation and cell fate,” Journal of Biological Chemistry, vol. 277, no. 14, pp. 11617–11620, 2002.

N. Ghanem, M. G. Andrusiak, D. Svoboda et al., “The Rb/E2F pathway modulates neurogenesis through direct regulation of the Dlx1/Dlx2 bigene cluster,” Journal of Neuroscience, vol. 32, no. 24, pp. 8219–8230, 2012.

F. Yan, H. Liu, J. Hao, and Z. Liu, “Dynamical Behaviors of Rb-E2F Pathway Including Negative Feedback Loops Involving miR449,” PLoS ONE, vol. 7, no. 9, Article ID e43908, 2012.

M. N. Obeyesekere, S. O. Zimmerman, E. S. Tecarro, and G. Auchmuty, “A model of cell cycle behavior dominated by kinetics of a pathway stimulated by growth factors,” Bulletin of Mathematical Biology, vol. 61, no. 5, pp. 917–934, 1999.

G. Yao, T. J. Lee, S. Mori, J. R. Nevins, and L. You, “A bistable Rb-E2F switch underlies the restriction point,” Nature Cell Biology, vol. 10, no. 4, pp. 476–482, 2008.

H. Chen, P. A. Mundra, L. N. Zhao, F. Lin, and J. Zheng, “Highly sensitive inference of time-delayed gene regulation by network deconvolution,” BMC Systems Biology, vol. 8, no. 4, article no. S6, 2014.

J.-P. Richard, “Time-delay systems: an overview of some recent advances and open problems,” Automatica. A Journal of IFAC, the International Federation of Automatic Control, vol. 39, no. 10, pp. 1667–1694, 2003.

C. Y. Ko, M. B. Liu, Z. Song, Z. Qu, and J. N. Weiss, “Multiscale Determinants of Delayed Afterdepolarization Amplitude in Cardiac Tissue,” Biophysical Journal, vol. 112, no. 9, pp. 1949–1961, 2017.

L. Mier-y-Tern-Romero, M. Silber, and V. Hatzimanikatis, “The origins of time-delay in template biopolymerization processes,” PLoS Computational Biology, vol. 6, no. 4, e1000726, 15 pages, 2010.

Q. Zheng, Z. Wang, and J. Shen, “Pattern dynamics of network-organized system with cross-diffusion,” Chinese Physics B, vol. 26, no. 2, p. 020501, 2017.

Q. Zheng and J. Shen, “Dynamics and pattern formation in a cancer network with diffusion,” Communications in Nonlinear Science and Numerical Simulation, vol. 27, no. 1-3, pp. 93–109, 2015.

Y. Xu, Y.-N. Zhu, J. W. Shen, and J. B. Su, “Switch dynamics for stochastic model of genetic toggle switch,” Physica A: Statistical Mechanics and Its Applications, vol. 416, pp. 461–466, 2014.

Q. Zheng and J. Shen, “Pattern formation in the FitzHugh—Nagumo model,” Computers & Mathematics with Applications, vol. 70, no. 5, pp. 1082–1097, 2015.

Q. Zheng and J. Shen, “Turing instability in a gene network with cross-diffusion,” Nonlinear Dynamics. An International Journal of Nonlinear Dynamics and Chaos in Engineering Systems, vol. 78, no. 2, pp. 1301–1310, 2014.

Q. Zheng and J. Shen, “Bifurcations and dynamics of cancer signaling network regulated by MicroRNA,” Discrete Dynamics in Nature and Society, vol. 2013, Article ID 176956, 2013.

G. Orosz, J. Moehlis, and R. M. Murray, “Controlling biological networks by time-delayed signals,” Philosophical Transactions of the Royal Society of London. Series A. Mathematical, Physical and Engineering Sciences, vol. 368, no. 1911, pp. 439–454, 2010.

Z. Song, Z. Qu, and A. Karma, “Stochastic initiation and termination of calcium-mediated triggered activity in cardiac myocytes,” Proceedings of the National Academy of Sciences, vol. 114, no. 3, pp. E270–E279, 2017.

P. Zoppoli, S. Morganella, and M. Ceccarelli, “Timedelay-ARACNE: reverse engineering of gene networks from time-course data by an information theoretic approach,” BMC Bioinformatics, vol. 11, no. 1, article 154, 2010.

Y. Li and A. Ngom, “The max-min high-order dynamic Bayesian network learning for identifying gene regulatory networks from time-series microarray data,” in Proceedings of the 10th Annual IEEE Symposium on Computational Intelligence in Bioinformatics and Computational Biology, CIBCB 2013 - 2013 IEEE Symposium Series on Computational Intelligence, SSCI 2013, pp. 83–90, April 2013.

Y. Cao and P. M. Frank, “Analysis and synthesis of nonlinear time-delay systems via fuzzy control approach,” IEEE Transactions on Fuzzy Systems, vol. 8, no. 2, pp. 200–211, 2000.

Z. Song, C. Y. Ko, M. Nivala, J. N. Weiss, and Z. Qu, “Complex darly and delayed afterdepolarization dynamics caused by voltage-calcium coupling in cardiac myocytes,” Biophysical Journal, vol. 108, no. 8, pp. 261A–262A, 2015.

C. Y. Ko, Z. Song, Z. Qu, and J. N. Weiss, “Multiscale Consequences of Spontaneous Calcium Release on Cardiac Delayed Afterdepolarizations,” Biophysical Journal, vol. 108, no. 2, p. 264a, 2015.