# American Institute of Mathematical Sciences

• Previous Article
Sufficient optimality conditions for a class of epidemic problems with control on the boundary
• MBE Home
• This Issue
• Next Article
Newton's method for nonlinear stochastic wave equations driven by one-dimensional Brownian motion
February  2017, 14(1): 249-262. doi: 10.3934/mbe.2017016

## Spatio-temporal models of synthetic genetic oscillators

 School of Mathematics and Statistics, Mathematical Institute, North Haugh, University of St Andrews, St Andrews KY16 9SS, Scotland

* Corresponding author: Mark Chaplain

Received  November 20, 2015 Accepted  April 15, 2016 Published  October 2016

Signal transduction pathways play a major role in many important aspects of cellular function e.g. cell division, apoptosis. One important class of signal transduction pathways is gene regulatory networks (GRNs). In many GRNs, proteins bind to gene sites in the nucleus thereby altering the transcription rate. Such proteins are known as transcription factors. If the binding reduces the transcription rate there is a negative feedback leading to oscillatory behaviour in mRNA and protein levels, both spatially (e.g. by observing fluorescently labelled molecules in single cells) and temporally (e.g. by observing protein/mRNA levels over time). Recent computational modelling has demonstrated that spatial movement of the molecules is a vital component of GRNs and may cause the oscillations. These numerical findings have subsequently been proved rigorously i.e. the diffusion coefficient of the protein/mRNA acts as a bifurcation parameter and gives rise to a Hopf bifurcation. In this paper we first present a model of the canonical GRN (the Hes1 protein) and show the effect of varying the spatial location of gene and protein production sites on the oscillations. We then extend the approach to examine spatio-temporal models of synthetic gene regulatory networks e.g. n-gene repressilators and activator-repressor systems.

Citation: Cicely K. Macnamara, Mark A. J. Chaplain. Spatio-temporal models of synthetic genetic oscillators. Mathematical Biosciences & Engineering, 2017, 14 (1) : 249-262. doi: 10.3934/mbe.2017016
##### References:

show all references

##### References:
Plots showing (a) the hes1 mRNA concentration and (b) the Hes1 protein concentration varying in space and time. The hes1 mRNA gene-site where transcription occurs is located at $x_{m}=0.0$ and the Hes1 protein production sites where translation occurs are located at $x_{p}=\pm0.5$. The plots show sustained oscillations of both hes1 mRNA and Hes1 protein. Baseline parameter set $\mathbb{P}$
Plots showing (a) the total hes1 mRNA concentration and (b) the total Hes1 protein concentration over time. The hes1 mRNA gene-site where transcription occurs is located at $x_{m}=0.0$ and the Hes1 protein production sites where translation occurs are located at $x_{p}=\{\pm0.1, \pm0.3, \pm0.5, \pm0.7, \pm0.9\}$ (see legend). The plots show that if the Hes1 protein production sites are either too close or too far away from the hes1 mRNA gene-site, then oscillations are lost. Baseline parameter set $\mathbb{P}$
Plots showing (a) the total hes1 mRNA concentration and (b) the total Hes1 protein concentration over time. The hes1 mRNA gene-site where transcription occurs is located at $x_{m}=0.0$ and the Hes1 protein production sites where translation occurs are located at $x_{p}=\{\pm0.1, \pm0.3, \pm0.5, \pm0.7, \pm0.9\}$ (see legend). The plots show that if the Hes1 protein production sites are either too close or too far away from the hes1 mRNA gene-site, then oscillations are lost. Parameters are as for the baseline parameter set $\mathbb{P}$, except $\epsilon$ which is taken to be $0.25$
Plots showing (a) the total hes1 mRNA concentration and (b) the total Hes1 protein concentration over time. The hes1 mRNA gene-site where transcription occurs is located at $x_{m}=0.0$. Translation occurs uniformly in the cytoplasm as per equations (3) and (4), with the boundaries beyond which protein production occurs located at $x_{p}=\{0.1, 0.3, 0.5, 0.7, 0.9\}$ (see legend). The plots show that if Hes1 protein production is either too close or too far away from the hes1 mRNA gene-site, then oscillations are lost. Parameters are as for the baseline parameter set $\mathbb{P}$
Plots showing (a) the mRNA concentration and (b) the protein concentration varying in space and time for the three-gene repressilator system. $m_1, p_1$ top row; $m_2, p_2$ middle row; $m_3, p_3$ bottom row. The mRNA gene-sites are in different locations, i.e. $x_{m1}=0.0$, $x_{m2}=\pm0.2$, $x_{m3}=\pm0.4$, and the protein production sites are in different locations, i.e. $x_{p1}=\pm0.5$, $x_{p2}=\pm0.7$ and $x_{p3}=\pm0.9$. The plots show sustained oscillations of each of the three mRNA and protein species. Baseline parameter set $\mathbb{P}$
Plots showing (a) the total mRNA concentration for species 1 and (b) the total protein concentration for species 1 over time for the three-gene repressilator. Case (a) -solid black curve: the gene sites are $x_{m1}=x_{m2}=x_{m3}=0.0$ and protein production sites are $x_{p1}=x_{p2}=x_{p3}=\pm0.5$. Case (b) -blue dotted curve: the gene sites are $x_{m1}=0.0$, $x_{m2}=\pm0.2$, $x_{m2}=\pm0.4$, and protein production sites are $x_{p1}=\pm0.5$, $x_{p2}=\pm0.7$ and $x_{p3}=\pm0.9$. Case (c) -red dashed curve: the gene sites are $x_{m1}=x_{m2}=x_{m3}=0.0$ and protein production sites are $x_{p1}=x_{p2}=x_{p3}=\pm0.9$. Case (d) -cyan dot-dashed curve: the gene sites are $x_{m1}=0.0$, $x_{m2}=\pm0.05$, $x_{m3}=\pm0.1$, and protein production sites are $x_{p1}=\pm0.9$, $x_{m2}=\pm0.95$ and $x_{p3}=\pm1.0$. Baseline parameter set $\mathbb{P}$
Plots showing (a) the total mRNA concentration for species 1 and (b) the total protein concentration for species 1 over time for the three-gene repressilator. For all cases protein production (translation) is constant throughout the cytoplasm $|x|>0.5$. Case (a) -solid black curve: the gene sites are $x_{m1}=x_{m2}=x_{m3}=0.0$. Case (b) -blue dotted curve: the gene sites are $x_{m1}=x_{m2}=x_{m2}=\pm0.2$. Case (c) -red dashed curve: the gene sites are $x_{m1}=0.0$, $x_{m2}=\pm0.2$ and $x_{m3}=\pm0.4$. Case (d) -cyan dot-dashed curve: the gene sites are $x_{m1}=\pm0.1$, $x_{m2}=\pm0.3$ and $x_{m3}=\pm0.5$. Baseline parameter set $\mathbb{P}$
Plots showing (a) the mRNA concentration and (b) the protein concentration varying in space and time for a two-gene activator-repressor system. $m_1, p_1$ top row; $m_2, p_2$ bottom row. The mRNA gene-sites where transcription occurs are located at $x_{m_i}=0.0, i=1, 2$ and the protein production sites where translation occurs are located at $x_{p_i}= \pm0.9, i=1, 2$. The plots show sustained oscillations of mRNA and protein concentrations for both species. Baseline parameter set $\mathbb{P}$, $\beta_m=10.0$
Plots showing (a) the total mRNA concentration for species 1 and (b) the total protein concentration for species 1 over time for the two gene activator-repressor system. The mRNA gene-sites where transcription occurs are located at $x_{m_i}=0.0, i=1, 2$ and the protein production sites where translation occurs are located at $x_{p_i}= \pm0.8, \pm0.88, \pm0.94, \pm1.0, i=1, 2$ [(a) solid line, black, (b) dotted line, blue, (c) dashed line, red, (d) dot-dashed line, cyan, respectively]. The plots show that if the protein production sites are either too close or too far away from the mRNA gene-sites, then sustained oscillations are lost. Baseline parameter set $\mathbb{P}$, $\beta_m=10.0$
 [1] Ö. Uğur, G. W. Weber. Optimization and dynamics of gene-environment networks with intervals. Journal of Industrial & Management Optimization, 2007, 3 (2) : 357-379. doi: 10.3934/jimo.2007.3.357 [2] Kunwen Wen, Lifang Huang, Qiuying Li, Qi Wang, Jianshe Yu. The mean and noise of FPT modulated by promoter architecture in gene networks. Discrete & Continuous Dynamical Systems - S, 2019, 12 (7) : 2177-2194. doi: 10.3934/dcdss.2019140 [3] Qingqing Li, Tianshou Zhou. Interlocked multi-node positive and negative feedback loops facilitate oscillations. Discrete & Continuous Dynamical Systems - B, 2019, 24 (7) : 3139-3155. doi: 10.3934/dcdsb.2018304 [4] William Chad Young, Adrian E. Raftery, Ka Yee Yeung. A posterior probability approach for gene regulatory network inference in genetic perturbation data. Mathematical Biosciences & Engineering, 2016, 13 (6) : 1241-1251. doi: 10.3934/mbe.2016041 [5] Pavol Bokes. Maintaining gene expression levels by positive feedback in burst size in the presence of infinitesimal delay. Discrete & Continuous Dynamical Systems - B, 2019, 24 (10) : 5539-5552. doi: 10.3934/dcdsb.2019070 [6] Roberto Serra, Marco Villani, Alex Graudenzi, Annamaria Colacci, Stuart A. Kauffman. The simulation of gene knock-out in scale-free random Boolean models of genetic networks. Networks & Heterogeneous Media, 2008, 3 (2) : 333-343. doi: 10.3934/nhm.2008.3.333 [7] Baltazar D. Aguda, Ricardo C.H. del Rosario, Michael W.Y. Chan. Oncogene-tumor suppressor gene feedback interactions and their control. Mathematical Biosciences & Engineering, 2015, 12 (6) : 1277-1288. doi: 10.3934/mbe.2015.12.1277 [8] Georg Hetzer, Anotida Madzvamuse, Wenxian Shen. Characterization of turing diffusion-driven instability on evolving domains. Discrete & Continuous Dynamical Systems - A, 2012, 32 (11) : 3975-4000. doi: 10.3934/dcds.2012.32.3975 [9] Lan Qiao, Sining Zheng. Non-simultaneous blow-up for heat equations with positive-negative sources and coupled boundary flux. Communications on Pure & Applied Analysis, 2007, 6 (4) : 1113-1129. doi: 10.3934/cpaa.2007.6.1113 [10] Guy Katriel. Stability of synchronized oscillations in networks of phase-oscillators. Discrete & Continuous Dynamical Systems - B, 2005, 5 (2) : 353-364. doi: 10.3934/dcdsb.2005.5.353 [11] Jian Ren, Feng Jiao, Qiwen Sun, Moxun Tang, Jianshe Yu. The dynamics of gene transcription in random environments. Discrete & Continuous Dynamical Systems - B, 2018, 23 (8) : 3167-3194. doi: 10.3934/dcdsb.2018224 [12] Yong Zhao, Qishao Lu. Periodic oscillations in a class of fuzzy neural networks under impulsive control. Conference Publications, 2011, 2011 (Special) : 1457-1466. doi: 10.3934/proc.2011.2011.1457 [13] Julien Coatléven, Claudio Altafini. A kinetic mechanism inducing oscillations in simple chemical reactions networks. Mathematical Biosciences & Engineering, 2010, 7 (2) : 301-312. doi: 10.3934/mbe.2010.7.301 [14] Elisabeth Logak, Isabelle Passat. An epidemic model with nonlocal diffusion on networks. Networks & Heterogeneous Media, 2016, 11 (4) : 693-719. doi: 10.3934/nhm.2016014 [15] Ying Hao, Fanwen Meng. A new method on gene selection for tissue classification. Journal of Industrial & Management Optimization, 2007, 3 (4) : 739-748. doi: 10.3934/jimo.2007.3.739 [16] Qi Wang, Lifang Huang, Kunwen Wen, Jianshe Yu. The mean and noise of stochastic gene transcription with cell division. Mathematical Biosciences & Engineering, 2018, 15 (5) : 1255-1270. doi: 10.3934/mbe.2018058 [17] Zizi Wang, Zhiming Guo, Huaqin Peng. Dynamical behavior of a new oncolytic virotherapy model based on gene variation. Discrete & Continuous Dynamical Systems - S, 2017, 10 (5) : 1079-1093. doi: 10.3934/dcdss.2017058 [18] Somkid Intep, Desmond J. Higham. Zero, one and two-switch models of gene regulation. Discrete & Continuous Dynamical Systems - B, 2010, 14 (2) : 495-513. doi: 10.3934/dcdsb.2010.14.495 [19] Marek Bodnar. Distributed delays in Hes1 gene expression model. Discrete & Continuous Dynamical Systems - B, 2019, 24 (5) : 2125-2147. doi: 10.3934/dcdsb.2019087 [20] Feng Jiao, Qiwen Sun, Genghong Lin, Jianshe Yu. Distribution profiles in gene transcription activated by the cross-talking pathway. Discrete & Continuous Dynamical Systems - B, 2019, 24 (6) : 2799-2810. doi: 10.3934/dcdsb.2018275

2018 Impact Factor: 1.313