Parametric sensitivity analysis of oscillatory delay systems with an application to gene regulation.

*(English)*Zbl 1372.92029Summary: A parametric sensitivity analysis for periodic solutions of delay-differential equations is developed. Because phase shifts cause the sensitivity coefficients of a periodic orbit to diverge, we focus on sensitivities of the extrema, from which amplitude sensitivities are computed, and of the period. Delay-differential equations are often used to model gene expression networks. In these models, the parametric sensitivities of a particular genotype define the local geometry of the evolutionary landscape. Thus, sensitivities can be used to investigate directions of gradual evolutionary change. An oscillatory protein synthesis model whose properties are modulated by RNA interference is used as an example. This model consists of a set of coupled delay-differential equations involving three delays. Sensitivity analyses are carried out at several operating points. Comments on the evolutionary implications of the results are offered.

##### MSC:

92C40 | Biochemistry, molecular biology |

##### Keywords:

sensitivity analysis; periodic solutions; delay-differential equations; gene expression; RNA interference##### Software:

DDE-BIFTOOL
PDF
BibTeX
XML
Cite

\textit{B. Ingalls} et al., Bull. Math. Biol. 79, No. 7, 1539--1563 (2017; Zbl 1372.92029)

Full Text:
DOI

##### References:

[1] | Acerenza, L; Sauro, HM; Kacser, H, Control analysis of time-dependent metabolic systems, J Theor Biol, 137, 423-444, (1989) |

[2] | Ahsen, ME; Özbay, H; Niculescu, SI, Analysis of deterministic cyclic gene regulatory network models with delays, Birkhäuser, Cham., (2010) · Zbl 1315.92003 |

[3] | Baek, D; Villén, J; Shin, C; Camargo, FD; Gygi, SP; Bartel, DP, The impact of micrornas on protein output, Nature, 455, 64-71, (2008) |

[4] | Baker CT, Rihan FA (1999) Sensitivity analysis of parameters in modelling with delay-differential equations. Tech. Rep. 349, Manchester Centre for Computational Mathematics |

[5] | Banks H, Robbins D, Sutton KL (2013a) Generalized sensitivity analysis for delay differential equations. In: Control and Optimization with PDE Constraints, Springer, pp 19-44 · Zbl 1279.34076 |

[6] | Banks, HT; Robbins, D; Sutton, KL, Theoretical foundations for traditional and generalized sensitivity functions for nonlinear delay differential equations, Math Biosci Eng, 10, 1301-1333, (2013) · Zbl 1273.34083 |

[7] | Bliss, RD; Painter, PR; Marr, AG, Role of feedback inhibition in stabilizing the classical operon, J Theor Biol, 97, 177-193, (1982) |

[8] | Bocharov, GA; Rihan, FA, Numerical modelling in biosciences using delay differential equations, J Comput Appl Math, 125, 183-199, (2000) · Zbl 0969.65124 |

[9] | Boström K, Wettesten M, Borén J, Bondjers G, Wiklund O, Olofsson SO (1986) Pulse-chase studies of the synthesis and intracellular transport of apolipoprotein B-100 in Hep G2 cells. J Biol Chem 261:13,800-13,806 |

[10] | Buchholtz, F; Schneider, FW, Computer simulation of T3/T7 phage infection using lag times, Biophys Chem, 26, 171-179, (1987) |

[11] | Bueler E, Butcher E (2002) Stability of periodic linear delay-differential equations and the Chebyshev approximation of fundamental solutions, preprint |

[12] | Bure, E; Rozenvasser, E, The study of the sensitivity of oscillatory systems, Autom Remote Control, 35, 1045-1052, (1974) · Zbl 0293.93005 |

[13] | Busenberg, SN; Mahaffy, JM, The effects of dimension and size for a compartmental model of repression, SIAM J Appl Math, 48, 882-903, (1988) · Zbl 0662.92011 |

[14] | Butcher, EA; Ma, H; Bueler, E; Averina, V; Szabo, Z, Stability of linear time-periodic delay-differential equations via Chebyshev polynomials, Int J Numer Methods Eng, 59, 895-922, (2004) · Zbl 1065.70002 |

[15] | Cinquin, O; Demongeot, J, Roles of positive and negative feedback in biological systems, C R Biol, 325, 1085-1095, (2002) |

[16] | Cooke, KL; Grossman, Z, Discrete delay, distributed delay and stability switches, J Math Anal Appl, 86, 592-627, (1982) · Zbl 0492.34064 |

[17] | Cornish-Bowden A, Cárdenas ML (eds) (1990) Control of metabolic processes. Plenum, New York |

[18] | Danø, S; Madsen, MF; Sørensen, PG, Chemical interpretation of oscillatory modes at a Hopf point, Phys Chem Chem Phys, 7, 1674-1679, (2005) |

[19] | Darzacq, X; Shav-Tal, Y; Turris, V; Brody, Y; Shenoy, SM; Phair, RD; Singer, RH, In vivo dynamics of RNA polymerase II transcription, Nat Struct Mol Biol, 14, 796-806, (2007) |

[20] | Dill, H; Linder, B; Fehr, A; Fischer, U, Intronic mir-26b controls neuronal differentiation by repressing its host transcript, ctdsp2, Genes Dev, 26, 25-30, (2012) |

[21] | Driver, RD, Existence and stability of solutions of a delay-differential system, Arch Ration Mech Anal, 10, 401-426, (1962) · Zbl 0105.30401 |

[22] | Ebert, MS; Sharp, PA, Roles for micrornas in conferring robustness to biological processes, Cell, 149, 515-524, (2012) |

[23] | Eden, E; Geva-Zatorsky, N; Issaeva, I; Cohen, A; Dekel, E; Danon, T; Cohen, L; Mayo, A; Alon, U, Proteome half-life dynamics in living human cells, Science, 331, 764-768, (2011) |

[24] | Edmunds LN Jr (1988) Cellular and molecular bases of biological clocks. Springer, New York |

[25] | Engelborghs, K; Luzyanina, T; Roose, D, Numerical bifurcation analysis of delay differential equations using DDE-BIFTOOL, ACM Trans Math Softw, 28, 1-21, (2002) · Zbl 1070.65556 |

[26] | Epstein, IR, Differential delay equations in chemical kinetics: some simple linear model systems, J Chem Phys, 92, 1702-1712, (1990) |

[27] | Ermentrout B (2002) Simulating, analyzing, and animating dynamical systems. SIAM, Philadelphia · Zbl 1003.68738 |

[28] | Falaleeva, M; Stamm, S, Processing of snornas as a new source of regulatory non-coding rnas, BioEssays, 35, 46-54, (2013) |

[29] | Fell, DA, Metabolic control analysis: a survey of its theoretical and experimental development, Biochem J, 286, 313-330, (1992) |

[30] | Feng, J; Sevier, SA; Huang, B; Jia, D; Levine, H, Modeling delayed processes in biological systems, Phys Rev E, 94, 032408, (2016) |

[31] | Ferrell, JE, Tripping the switch fantastic: how a protein kinase can convert graded inputs into switch-like outputs, Trends Biochem Sci, 21, 460-466, (1996) |

[32] | Ferrell, JE; Xiong, W, Bistability in cell signaling: how to make continuous processes discontinuous, and reversible processes irreversible, Chaos, 11, 227-236, (2001) · Zbl 0992.92023 |

[33] | Goldbeter A (1996) Biochemical oscillations and cellular rhythms. Cambridge University Press, Cambridge · Zbl 0837.92009 |

[34] | Goodwin BC (1963) Temporal organization in cells. Academic Press, London |

[35] | Halaney A (1966) Differential Equations. Academic Press, New York, Stability, Oscillations, Time Lags |

[36] | Hale, JK; Ladeira, LAC, Differentiability with respect to delays, J Differ Equ, 92, 14-26, (1991) · Zbl 0735.34045 |

[37] | Hale JK, Lunel SMV (1993) Introduction to functional differential equations. Springer, New York · Zbl 0787.34002 |

[38] | Heinrich, R; Reder, C, Metabolic control analysis of relaxation processes, J Theor Biol, 151, 343-350, (1991) |

[39] | Heinrich R, Schuster S (1996) The regulation of cellular systems. Chapman & Hall, New York · Zbl 0895.92013 |

[40] | Highkin, HR; Hanson, JB, Possible interaction between light-dark cycles and endogeneous daily rhythms on the growth of tomato plants, Plant Physiol, 29, 301-302, (1954) |

[41] | Hill AV (1910) The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves. J Physiol (Lond) 40:4-7. doi:10.1113/jphysiol.1910.sp001386 (Suppl.) |

[42] | Hillman, WS, Injury of tomato plants by continuous light and unfavorable photoperiodic cycles, Am J Bot, 43, 89-96, (1956) |

[43] | Ingalls, B, Sensitivity analysis: from model parameters to system behaviour, Essays Biochem, 45, 177-193, (2008) |

[44] | Ingalls, BP, Autonomously oscillating biochemical systems: parametric sensitivity of extrema and period, Syst Biol, 1, 62-70, (2004) |

[45] | Ingalls BP (2013) Mathematical modeling in systems biology. MIT Press, Cambridge · Zbl 1312.92003 |

[46] | Ingalls, BP; Sauro, HM, Sensitivity analysis of stoichiometric networks: an extension of metabolic control analysis to non-steady state trajectories, J Theor Biol, 222, 23-36, (2003) |

[47] | Izaurralde, E, Breakers and blockers-mirnas at work, Science, 349, 380-382, (2015) |

[48] | Johnson, CH, Endogenous timekeepers in photosynthetic organisms, Annu Rev Physiol, 63, 695-728, (2001) |

[49] | Kacser, H; Burns, JA, The control of flux, Symp Soc Exp Biol, 27, 65-104, (1973) |

[50] | Kholodenko, BN; Demin, OV; Westerhoff, HV, Control analysis of periodic phenomena in biological systems, J Phys Chem B, 101, 2070-2081, (1997) |

[51] | Klarsfeld, A; Rouyer, F, Effects of Circadian mutations and LD periodicity on the life span of drosophila melanogaster, J Biol Rhythms, 13, 471-478, (1998) |

[52] | Lander ES, Linton LM, Birren B, Nusbaum C et al (2001) Initial sequencing and analysis of the human genome. Nature 409:860-921. doi:10.1038/35057062, errata: Nature 412, 565-566 |

[53] | Lapidot, M; Pilpel, Y, Genome-wide natural antisense transcription: coupling its regulation to its different regulatory mechanisms, EMBO Rep, 7, 1216-1222, (2006) |

[54] | Larter, R, Sensitivity analysis of autonomous oscillators. separation of secular terms and determination of structural stability, J Phys Chem, 87, 3114-3121, (1983) |

[55] | Lenz, SM; Schlöder, JP; Bock, HG, Numerical computation of derivatives in systems of delay differential equations, Math Comput Simul, 96, 124-156, (2014) |

[56] | Lewis, J, Autoinhibition with transcriptional delay: a simple mechanism for the zebrafish somitogenesis oscillator, Curr Biol, 13, 1398-1408, (2003) |

[57] | Long X, Insperger T, Balachandran B (2009) Systems with periodic coefficients and periodically varying delays: semidiscretization-based stability analysis. In: Delay Differential Equations, Springer, pp 131-153 · Zbl 0293.93005 |

[58] | MacDonald, N, Time lag in a model of a biochemical reaction sequence with end product inhibition, J Theor Biol, 67, 549-556, (1977) |

[59] | MacDonald, N, An interference effect of independent delays, IEE Proc D, 134, 38-42, (1987) · Zbl 0604.34042 |

[60] | Meijer, HA; Kong, YW; Lu, WT; Wilczynska, A; Spriggs, RV; Robinson, SW; Godfrey, JD; Willis, AE; Bushell, M, Translational repression and eif4A2 activity are critical for microrna-mediated gene regulation, Science, 340, 82-85, (2013) |

[61] | Mello, CC; Conte, D, Revealing the world of RNA interference, Nature, 431, 338-342, (2004) |

[62] | Mier-y-Terán-Romero, L; Silber, M; Hatzimanikatis, V, The origins of time-delay in template biopolymerization processes, PLoS Comput Biol, 6, e1000726, (2010) |

[63] | Monk, NAM, Oscillatory expression of hes1, p53, and NF-\(κ \)B driven by transcriptional time delays, Curr Biol, 13, 1409-1413, (2003) |

[64] | Ouyang, Y; Andersson, CR; Kondo, T; Golden, SS; Johnson, CH, Resonating Circadian clocks enhance fitness in cyanobacteria, Proc Natl Acad Sci USA, 95, 8660-8664, (1998) |

[65] | Pittendrigh, CS; Minis, DH, Circadian systems: longevity as a function of Circadian resonance in drosophila melanogaster, Proc Natl Acad Sci USA, 69, 1537-1539, (1972) |

[66] | Purcell, O; Savery, NJ; Grierson, CS; Bernardo, M, A comparative analysis of synthetic genetic oscillators, J R Soc Interface, 7, 1503-1524, (2010) |

[67] | Rihan, FA, Sensitivity analysis for dynamic systems with time-lags, J Comput Appl Math, 151, 445-462, (2003) · Zbl 1023.93030 |

[68] | Roussel, CJ; Roussel, MR, Delay-differential equations and the model equivalence problem in chemical kinetics, Phys Can, 57, 114-120, (2001) |

[69] | Roussel, MR, The use of delay differential equations in chemical kinetics, J Phys Chem, 100, 8323-8330, (1996) |

[70] | Roussel, MR; Zhu, R, Validation of an algorithm for delay stochastic simulation of transcription and translation in prokaryotic gene expression, Phys Biol, 3, 274-284, (2006) |

[71] | Rozenwasser E, Yusupov R (1999) Sensitivity of automatic control systems. CRC Press, Boca Raton |

[72] | Saltelli A, Ratto M, Andres T, Campolongo F, Cariboni J, Gatelli D, Saisana M, Tarantola S (2008) Global sensitivity analysis: the primer. Wiley, Chichester · Zbl 1161.00304 |

[73] | Schmiedel, JM; Klemm, SL; Zheng, Y; Sahay, A; Blüthgen, N; Marks, DS; Oudenaarden, A, Microrna control of protein expression noise, Science, 348, 128-132, (2015) |

[74] | Shen, J; Liu, Z; Zheng, W; Xu, F; Chen, L, Oscillatory dynamics in a simple gene regulatory network mediated by small rnas, Phys A, 388, 2995-3000, (2009) |

[75] | Shimoni, Y; Friedlander, G; Hetzroni, G; Niv, G; Altuvia, S; Biham, O; Margalit, H, Regulation of gene expression by small non-coding RNAs: a quantitative view, Mol Syst Biol, 3, 138, (2007) |

[76] | Singh, J; Padgett, RA, Rates of in situ transcription and splicing in large human genes, Nat Struct Mol Biol, 16, 1128-1133, (2009) |

[77] | Smith H (2011) An introduction to delay differential equations with applications to the life sciences, texts in applied mathematics, vol 57. Springer, New York · Zbl 1227.34001 |

[78] | Smolen, P; Baxter, DA; Byrne, JH, Modeling transcriptional control in gene networks-methods, recent results, and future directions, Bull Math Biol, 62, 247-292, (2000) · Zbl 1323.92088 |

[79] | Stokes, A, A Floquet theory for functional differential equation, Proc Natl Acad Sci USA, 48, 1330-1334, (1962) · Zbl 0109.06301 |

[80] | Sweeney BM (1987) Rhythmic phenomena in plants, 2nd edn. Academic Press, San Diego |

[81] | Taylor, SR; Campbell, SA, Approximating chaotic saddles for delay differential equations, Phys Rev E, 75, 046215, (2007) |

[82] | Taylor, SR; Gunawan, R; Petzold, LR; Doyle, FJ, Sensitivity measures for oscillating systems: application to mammalian Circadian gene network, IEEE Trans Automat Control, 53, 177-188, (2008) · Zbl 1366.92018 |

[83] | Tian, T; Burrage, K; Burrage, PM; Carletti, M, Stochastic delay differential equations for genetic regulatory networks, J Comput Appl Math, 205, 696-707, (2007) · Zbl 1112.92029 |

[84] | Tigges, M; Marquez-Lago, TT; Stelling, J; Fussenegger, M, A tunable synthetic Mammalian oscillator, Nature, 457, 309-312, (2009) |

[85] | Tyson, JJ, Classification of instabilities in chemical reaction systems, J Chem Phys, 62, 1010-1015, (1975) · Zbl 0306.65080 |

[86] | Tyson, JJ; Csikasz-Nagy, A; Novak, B, The dynamics of cell cycle regulation, BioEssays, 24, 1095-1109, (2002) |

[87] | Varma A, Morbidelli M, Wu H (2005) Parametric sensitivity in chemical systems. Cambridge University Press, Cambridge |

[88] | Wang, Y; Liu, CL; Storey, JD; Tibshirani, RJ; Herschlag, D; Brown, PO, Precision and functional specificity in mrna decay, Proc Natl Acad Sci USA, 99, 5860-5865, (2002) |

[89] | Wilkins, AK; Tidor, B; White, J; Barton, PI, Sensitivity analysis for oscillating dynamical systems, SIAM J Sci Comput, 31, 2706-2732, (2009) · Zbl 1201.65122 |

[90] | Woelfle, MA; Ouyang, Y; Phanvijhitsiri, K; Johnson, CH, The adaptive value of Circadian clocks: an experimental assessment in cyanobacteria, Curr Biol, 14, 1481-1486, (2004) |

[91] | Yan, X; Hoek, TA; Vale, RD; Tanenbaum, ME, Dynamics of translation of single mrna molecules in vivo, Cell, 165, 976-989, (2016) |

[92] | Yanchuk, S; Perlikowski, P, Delay and periodicity, Phys Rev E, 79, 221, (2009) |

[93] | Zak, DE; Stelling, J; Doyle, FJ, Sensitivity analysis of oscillatory (bio)chemical systems, Comput Chem Eng, 29, 663-673, (2005) |

[94] | Zhang, HM; Kuang, S; Xiong, X; Gao, T; Liu, C; Guo, AY, Transcription factor and microrna co-regulatory loops: important regulatory motifs in biological processes and diseases, Br Bioinform, 16, 45-58, (2015) |

[95] | Zhang, Y; Liu, H; Zhou, J, Oscillatory expression in Escherichia coli mediated by micrornas with transcriptional and translational time delays, IET Syst Biol, 10, 203-209, (2016) |

[96] | Zhdanov, VP, Bistability in gene transcription: interplay of messenger RNA, protein, and nonprotein coding RNA, Biosystems, 95, 75-81, (2009) |

[97] | Zhdanov, VP, Kinetic models of gene expression including non-coding rnas, Phys Rep, 500, 1-42, (2011) |

[98] | Zi, Z, Sensitivity analysis approaches applied to systems biology models, IET Syst Biol, 5, 336-346, (2011) |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.