Bifurcations and Dynamics of Cancer Signaling Network Regulated by MicroRNA

MicroRNAs are currently considered as key regulators of a wide variety of biological pathways and regulate many processes of life and obtained more and more attention in recent years. In this paper, we investigate the dynamics of gene network regulated by miR-34a (microRNA) involved in triple negative feedback loop. As we know that the p53 network involve in cancer, How the cancer arise is unclear.We investigate this negative feedback network by usingmathematical model and drive the theoretical results of globally asymptotic stability and provide the sufficient conditions for the periodic oscillation. These results are propitious to understand how p53 network involved in miR-34a induces the cancer.


Introduction
MicroRNAs (miRNAs) are a family of small regulatory RNAs whose function is to regulate the activity and stability of specific mRNA targets through posttranscriptional regulatory mechanism and play a role in repressing translation of mRNA or degrading mRNAs [1][2][3][4][5][6][7][8][9].Recent studies show that microRNAs play a central role in many biological (cellular) processes, including developmental timing, cell proliferation, apoptosis, metabolism, cell differentiation, somitogenesis and tumour genesis [1].In this paper, we focus on miR-34a which may behave as inhibitor depending on biological context.Resecting normal tumor tissues of 25 human hepatocellular carcinoma patients demonstrated an inverse correlation between miR-34a and c-Met-protein [2].MiR-34a is proposed as a marker for the activity of the p53 pathway in chronic lymphocytic leukemia [3] (Table 1).The induction of miR-34a was most pronounced among all differential regulations.Also expression of the primary miR-34a transcript was induced after p53 activation [4].MiR-34a may play a cytoprotective role in cell survival [5].In U87 cell, miR-34a plays a negative role in the regulation of Dll1 target gene through downregulating Dll1 protein but not mRNA [6].
Expression of miR-34a is induced which may mediate the target gene and result in the cells arresting in G1 and G2 phase to repair DNA [7].
In order to understand further the miR-34a involved in the network with p53 and Sirt1, we have planned to model this network with mathematical model.It is well known that the time delay is quite ubiquitous in nature, so we also investigate the relationship between the time delay and the network with miR-34a, p53, and Sirt1.We know that the delay is often caused by a finite signal transmission speed and memory effects, so the time delay can sometimes destabilize the stable unique equilibrium.If the time delay reaches a threshold value, the system will generate the phenomenon with self-oscillation.In nature, the oscillation often occurs in physiological regulatory systems with time delay which can induce the complex behaviors.
In recent years, many scientists deemed that mathematical modeling could be used to investigate the differences at the dynamical level between healthy and pathologic configurations of biological pathways [10].By using the mathematical model, the researchers can detect the key points regulating main properties of biological system and find the methods to solve the different diseases.In this paper, we will investigate the dynamics of the gene network composed of miR-34a, p53, and Sirt1 and reveal how the dynamics of microRNA regulation is affected by time delay associated with translation degradation of mRNA.In Section 2, we give the gene network represented by mathematical model and some theoretical results.In Section 3, we show the numerical analysis of the network.Finally, we summarize the results.
Then, we have the following theoretical results.
The another case we can discuss as above, and ignore it here.

Numerical Analysis
In this section, we present some numerical results of system (1) to verify the analytical predictions obtained in Section 2. Without loss of generality, if we take  11 = 0.11,  12 = 1, which has a positive equilibrium  = (0.9969, 1.1192, 0.8869, 1.8869, 0.5260, 2) which satisfies the conditions indicated in Section 2. When  1 =  2 +  3 = 0,  is asymptotically stable (see Figure 2).When  <  0 = 4.9859, the system (17) is illustrated by the computer simulations (see Figure 2).When  passes through the critical value  0 , the positive equilibrium  loses its stability and the Hopf bifurcation occurs (see Figures 3  and 4).
Figure 3 shows the behaviors and phase portraits of system (17) with  1 =  2 +  3 = 2 <  0 , and the positive equilibrium  is asymptotically stable.
Figure 4 shows the behaviors and phase portraits of system (17) with  1 =  2 +  3 = 8 >  0 and Hopf bifurcation occurs from the positive equilibrium  as showed in Figure 5, we can know that the expression of protein and the oscillation occurs when the time delay overpasses the critical points.Above the critical value of time delay, the expression of small RNA also oscillates periodically and provides the capabilities inside and outside cell to communicate.
In system (1), the time delay approximately represents transportation or diffusion process from nucleus to cytoplasm of mRNA and from cytoplasm to nucleus of protein, respectively.These delays play different roles in the dynamical behaviors of the system (1).

Discussion
From the above discussion, we know that  2 and  3 affect the dynamical behaviors in the form of  2 +  3 due to the cyclic structure of system (1), and we used the sum of  2 and  3 as a bifurcation parameter and then found that the stability depends on the sum of delay.Here we just gave a very simplified model to describe the dynamical behaviors of p53, Sirt1, and miR-34a.In this paper, we have given a detailed analysis of mathematical properties of the system.Our results indicate that the steady state ( 10 ,  20 ,  30 ,  40 ,  50 ,  60 ) is asymptotically stable for all delay values under certain conditions.However, if these conditions are not satisfied, there exists a critical value of the time delay  0 .When  <  0 , the steady state is asymptotically stable; when  >  0 , there are periodical oscillations and the Hopf bifurcation appears.
From the above discussion, we understand how the microRNA regulate the negative feedback loop in cancer signalling network (Figure 1), so we can use these results to explain how the prostate cancer stem cells are regulated by miR-34 [12] and give a clear understanding on the possible reason of prostate cancer.
Finally, it is worth noting that microRNA-mediated regulation has gained recent attention, and computational studies have revealed various regulatory properties unique to microRNAs.These findings will be helpful for our understating of the operating mechanisms and biological implications of microRNA-mediated regulation.They also have great potential for biotechnological and therapeutic applications and synthetic biology.

Conclusions
We analyzed a simple model of the interactions between miR-34a and target protein p53 and Sirt1 and others.Our goal is to explore the oscillatory dynamics and how the microRNAs repress its target protein.Finally, we derive explicit conditions on how the dynamics of a time delay model of the interaction between the microRNA (miR-34a) cluster and p53 and Sirt1 depends on system parameters.Our analysis reveals the complex behavior of the network and there is a limit cycle after a Hopf bifurcation for the time delay parameter and it shows that the analytical results agree with numerical simulations.

Figure 1 :
Figure 1: Schematic diagram of the complex network which illustrates the mechanisms.(a) Original model and (b) abstract model.

Table 1 :
Notation for species concentrations.