Stability and Bifurcation Analysis of a Diffusive miR-9/Hes1 Network With Time Delay

• Published in: IEEE/ACM transactions on computational biology and bioinformatics Vol. 19; no. 3; pp. 1870 - 1880
• Main Authors: Li, Chengxian, Liu, Haihong, Zhang, Tonghua, Zhang, Yuan
• Format: Journal Article
• Language: English
• Published: United States IEEE 01.05.2022; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Algorithms; Bifurcation; Boundary conditions; Computer Simulation; Delay effects; Delays; Diffusion; Diffusion effects; Eigenvalues; Eigenvectors; Hes1; HES1 protein; Hopf bifurcation; Low level; Mathematical model; MicroRNAs - genetics; miR-9; Models, Biological; mRNA; Oscillators; Proteins; Stability analysis; Steady state; time delay; Time lag; Turing-Hopf bifurcation
• ISSN: 1545-5963; 1557-9964; 1557-9964
• DOI: 10.1109/TCBB.2021.3050006

In this paper, a model of miR-9/Hes1 interaction network involving one time delay and diffusion effect under the Neumann boundary conditions is studied. First of all, the stability of the positive equilibrium and the existence of local Hopf bifurcation and Turing-Hopf bifurcation are investigated by analyzing the associated characteristic equation. Second, a algorithm for determining the direction, stability and period of the corresponding bifurcating periodic solutions is presented. The obtained results suggest that the quiescent progenitors (high steady-state Hes1) can be easily excited into oscillation by time delay whereas the differentiated state (low steady-state Hes1) is basically unaffected, and the integrated effect of delay and diffusion can induce the occurrence of spatially inhomogeneous patterns. More importantly, spatially homogeneous/inhomogeneous periodic solutions can exist simultaneously when the diffusion coefficients of Hes1 mRNA and Hes1 protein are appropriately small, conversely, there is only spatially homogeneous periodic solutions. Intriguingly, both temporal patterns and spatial-temporal patterns show that time delay can prompt Hes1 protein to shift from the high concentration state to the low concentration one ("ON" <inline-formula><tex-math notation="LaTeX">\rightarrow</tex-math> <mml:math><mml:mo>→</mml:mo></mml:math><inline-graphic xlink:href="liu-ieq1-3050006.gif"/> </inline-formula> "OFF"), where Hes1 protein shows low level whereas miR-9 shows high level. Finally, some numerical examples are presented to verify and visualize theoretical results.