Supporting Online Material for
|
|
- Joel Atkins
- 5 years ago
- Views:
Transcription
1 Supporting Online Material for Robust, Tunable Biological Oscillations from Interlinked Positive and Negative Feedback Loops Tony Yu-Chen Tsai, Yoon Sup Choi, Wenzhe Ma, Joseph R. Pomerening, Chao Tang, James E. Ferrell Jr.* *To whom correspondence should be addressed. This PDF file includes: Published 4 July 28, Science 321, 126 (28) DOI: /science SOM Text Fig. S1 Table S1 References
2 Supporting Online Material Robust, tunable biological oscillations from interlinked positive and negative feedback loops Tony Yu-Chen Tsai, Yoon Sup Choi, Wenzhe Ma, Joseph R. Pomerening, Chao Tang, and James E. Ferrell, Jr. Computational modeling For each ordinary differential equation (ODE) model analyzed, we solved the differential equations numerically using a Runge-Kutta method. To construct amplitude vs. frequency curves, we first chose a bifurcation parameter and then identified the range of the parameter over which the system exhibited limit cycle oscillations by locating the Hopf bifurcations at which oscillations were born and extinguished. We then used an iterative algorithm to march along the chosen parameter between the two bifurcations. For each set of parameters we calculated the limit cycle solution of the system, keeping track of the amplitude and frequency of the limit cycle solution. We made use of Mathematica (Wolfram Research, Inc.) and MatCONT (S1), a public toolbox for bifurcation analysis in Matlab (The MathWorks Inc.), for these computations. Below we describe the ODE models, parameter sets, and bifurcation parameters we examined. 1. Negative feedback and positive-plus-negative feedback cell cycle oscillators. This is a nine ODE model of the embryonic mitotic oscillator in Xenopus laevis (S2). Depending upon parameter choice, the model can function as a negative feedback circuit or as a positive-plusnegative feedback circuit. Parameters were chosen so that at a realistic feedback strength of r, the model would reproduce the experimental curves for the steady-state activity of CDK1 as a function of constant cyclin concentration, and for the time course of CDK1 activation when driven by a constant rate of cyclin synthesis (S2, S3). To make it easier to distinguish between complexes (e.g. CDK1-cyclin) and products of individual species (e.g. CDK1 cyclin), for this model we have written all time-dependent species as explicit functions of time. Thus, cdk1cyclin represents the concentration of the CDK1-cyclin complex, whereas cdk1cyclin represents the product of the free CDK1 and cyclin concentrations. dcyclin = k synth " k dest apc act cyclin " ( cdk1 tot " cdk1cyclin " cdk1cyclinyp " cdk1cyclinyptp " cdk1cyclintp )cyclin + k a k d cdk1cyclin 1
3 dcdk1cyclin = ( cdk1 tot " cdk1cyclin " cdk1cyclinyp " cdk1cyclinyptp " cdk1cyclintp )cyclin " k a k d cdk1cyclin " k dest apc act cdk1cyclin " k wee1 wee1 act cdk1cyclin " k wee1basal ( wee1 tot " wee1 act )cdk1cyclin + k cdc25 cdc25 act cdk1cyclinyp + ( cdc25 tot " cdc25 act )cdk1cyclinyp k cdc 25basal dcdk1cyclinyp = k wee1 wee1 act cdk1cyclin + k wee1basal ( wee1 tot " wee1 act )cdk1cyclin " k cdc 25 cdc25 act cdk1cyclinyp " k cdc25basal ( cdc25 tot " cdc25 act )cdk1cyclinyp " k cak cdk1cyclinyp + k pp 2c cdk1cyclinyptp " k dest apc act cdk1cyclinyp dcdk1cyclinyptp = k cak cdk1cyclinyp " k pp 2c cdk1cyclinyptp " k cdc 25 cdc25 act cdk1cyclinyptp " k cdc25basal ( cdc25 tot " cdc25 act )cdk1cyclinyptp + k wee1 wee1 act cdk1cyclintp + k wee1basal ( wee1 tot " wee1 act )cdk1cyclintp " k dest apc act cdk1cyclinyptp dcdk1cylintp = k cdc25 cdc25 act cdk1cyclinyptp + k cdc 25basal ( cdc25 tot " cdc25 act )cdk1cyclinyptp " k wee1 wee1 act cdk1cyclintp " ( wee1 tot " wee1 act )cdk1cyclintp " k dest apc act cdk1cyclintp k wee1basal dcdc25 act " cdk1cyclintp ncdc25 % = k cdc25on $ ncdc ' cdc25 # 25 + cdk1cyclintp ncdc25 tot ( cdc25 act & k cdc 25off cdc25 act dwee1 act dplx1 dapc act ec5 cdc 25 ( ) ( # cdk1cyclintp nwee1 & = "k wee1off % ( wee1 $ nwee1 + cdk1cyclintp nwee1 act + k wee1on wee1 tot " wee1 act ' ec5 wee1 " cdk1cyclintp nplx1 % = k plx1on $ # ec5 nplx1 plx1 + cdk1cyclintp nplx1 ' plx1 tot ( plx1 act & Initial conditions: cyclin[] = cdk1cyclin[] = cdk1cyclinyp[] = " plx1 = k act napc % apcon $ # ec5 napc plx1 + plx1 act napc ' apc tot ( apc act & ( ) ( k plx1off plx1 act ( ) ( k apcoff apc act ( ) 2
4 cdk1cyclinyptp[] = cdk1cyclintp[] = cdc25 act [] = wee1 act [] = wee1 tot plx act [] = apc act [] = Parameters: k dest =.1 k a =.1 k d =.1 k wee1 =.5 k wee1basal = k wee1 /r k cdc 25 =.1 k cdc 25basal = k cdc25 /r cdc2 tot = 23 cdc25 tot 5 wee1 tot 5 apc tot = 5 plx1 tot = 5 nwee1 = 4 ncdc25 = 4 napc = 4 nplx1= 4 ec5 plx1 = 4 ec5 wee1 = 4 ec5 cdc25 = 4 ec5 apc = 4 k cdc 25on.75 k cdc 25off =.2 k apcon k apcoff =.15 k plx1on 3
5 k plxoff =.15 k wee1on =.2 k wee1off.75 k cak =.8 k pp 2c =.8 For the negative feedback model, we took r = 1, which is equivalent to assuming that the activities of Wee1 and Cdc25 are not affected by the activity of CDK1. For the positive-plusnegative feedback model, we took r = 5, 1, or 15 to compare the effects of different feedback strengths. In all cases, the oscillations were found to be born and extinguished through supercritical Hopf bifurcations. Thus a change in the type of bifurcation was not responsible for the change in the shape of the amplitude vs. frequency curves that positive feedback produced. The following five models are negative-feedback-only oscillators: 2. Repressilator. Elowitz and Leibler used a six variable ODE model to guide their construction of the Repressilator, an artificial negative feedback oscillator, in E. coli (S4). The model consists of three mrnas, m1, m2, and m3, which give rise to three proteins, p1, p2, and p3. Each protein inhibits the transcription of the next message (e.g. p1 represses m2). We varied the frequency of the oscillator by varying one of the translation rates by a factor φ. Oscillations were born and extinguished through supercritical Hopf bifurcations. The parameters for the model are as given by Elowitz and Leibler (S4). dm 1 dm 2 dm 3 = "m 1 + # 1+ p 3 n + # = "m 2 + # 1+ p 1 n + # = "m 3 + # 1+ p 2 n + # dp 1 = " ( m 1 # $p 1 ) dp 2 dp 3 = "( m 2 # p 2 ) = "( m 3 # p 3 ) 4
6 Initial conditions: m 1 [] m 2 [] = m 3 [] = p 1 [] = p 2 [] = p 3 [] = Parameters: n = 2 " = 3 " = # =.2 3. Pentilator. If the time scales of transcription and translation are separable, the Repressilator can be reduced to three ODEs. Since a two-ode negative feedback loop cannot generate sustained oscillations, it seemed reasonable that a three ODE system would require careful balancing of kinetic parameters to oscillate, and would therefore oscillate over a limited range of frequencies. We were therefore curious to see if a five-ode analog of the repressilator (the Pentilator ) might have a qualitatively different amplitude/frequency relationship compared to the Repressilator and the other three ODE negative feedback oscillators being examined here. However, as shown in Figure 2, it turned out to behave very similarly to the Repressilator: oscillations were born and extinguished through supercritical Hopf bifurcations, and the operational frequency range of the model was small. dm 1 dm 2 dm 3 dm 4 dm 5 = "m 1 + # 1+ p 5 n + # = "m 2 + # 1+ p 1 n + # = "m 3 + # 1+ p 2 n + # = "m 4 + # 1+ p 3 n + # = "m 5 + # 1+ p 4 n + # dp 1 = " ( m 1 # $p 1 ) dp 2 = "( m 2 # p 2 ) 5
7 dp 3 dp 4 dp 5 = "( m 3 # p 3 ) = "( m 4 # p 4 ) = "( m 5 # p 5 ) Initial conditions: m 1 [] m 2 [] = m 3 [] = m 4 [] = m 5 [] = p 1 [] = p 2 [] = p 3 [] = p 4 [] = p 5 [] = Parameters: n = 2 " = 3 " = # =.2 4. Goodwin oscillator. The Goodwin model is a three-variable ODE system proposed forty years ago (S5), which has been applied to the analysis of circadian oscillations (S6-S1). The parameters used here were taken from Ruoff and co-workers (S8). Because they have emphasized that the frequency of the oscillator is a relatively sensitive function of the degradation rates, we chose the degradation rate constant to be the bifurcation parameter φ. Note that for some of the parameter sets used by Ruoff and co-workers, the Goodwin model produces damped oscillations that approach a stable steady state (or stable spiral) rather than a stable limit cycle (S8). Here we examined only limit cycle oscillations. dx = k 1 1+ z 9 " k 4 x dy = k 2 x " #y dz = k 3y " k 6 z 6
8 Initial conditions: x[] =.348 y[] =.347 z[].735 Parameters: k 1 k 2 k 3 k 4 =.2 k 6 =.1 5. Frzilator. The Frzilator was proposed by Igoshin and co-workers to account for oscillatory swarming motions in Myxococcus (S11). The model consists of three ODEs and assumes that all of the activating and inactivating enzymes are highly saturated by their substrates. This model differs from the other negative feedback models considered here in that oscillations are born and extinguished through subcritical Hopf bifurcations rather than supercritical Hopf bifurcations. This means that there is a non-zero amplitude below which the model cannot generate sustained oscillations. It also means that near the bifurcations, a locally stable steady-state (or spiral) coexists with limit cycle oscillations, and the initial conditions determine which behavior is exhibited by the model. df = " $ 1# f ' $ f ' & ) # d f & ) e %.1+ (1# f )( %.5 + f ( dc = k $ 1# c ' $ c ' c& ) # d c & ) %.5 + (1# c) ( %.5 + c ( de = k $ 1# e ' $ e ' e& ) # d e & ) %.5 + (1# e) ( %.5 + e( Initial conditions: f [] =.53 c[] =.623 e[] =.98 7
9 Parameters: k c = 4 k e = 4 d f d c = 2 d e = 2 6. Metabolator. The Metabolator is a synthetic gene-metabolic oscillator that rewires the glycolytic flux to drive an oscillation between two pools of metabolites in acetate pathway (S12). The design is based on two negative feedback loops. When acetyl phosphate (AcP) level increases, acetyl-coa synthetase (Acs) is activated to drive the flux away from AcP, and phosphate acetyltransferase (Pta) is repressed to decrease the production of AcP. Here we vary the enzyme copy number φ in the model to change the frequency of the oscillations. d[accoa] = V Acs "V pta + V gly " k TCA [AcCoA] d[acp] = V pta "V Ack d[oac " ] = V Ack " k 3 [HOAc] d[hoac] = V AcE " k 3 [HOAc] d[laci] d[aca] d[ pta] ( ) n ( ) + $ " d[laci] n ( ) n ( ) + a 1$ n " d[acs] = # [AcP]/K g,1 1+ [AcP]/K g,1 = a 1# [AcP]/K g,2 1+ [AcP]/K g,2 a = 2 # 1+ [laci]/k g,3 ( ) n + a 2$ " d[pta] where V pta = k 1 [ pta][accoa] K m1 + [AcCoA] V Acs = k 2 [Acs][OAc" ] K m 2 + [OAc " ] 8
10 V Ack = k Ack, f [AcP] " k Ack,r [OAc " ] V AcE ([OAc " ][H + ] " K eq [HOAc]) V out = k 3 ([HOAc] "[HOAc E ]) Initial conditions: AcCoA[] = AcP[] = 2 OAc " [] = HOAc[] =.9 laci[] =.55 Acs[] =.11 pta[] =.844 Parameters: k TCA k 1 = 8 k 2 =.8 K m1 =.6 K m2 =.1 k Ack, f = k Ack,r [H + ] "7.6 K eq "4.5 k 3 =.1 a = a 1 = a 2 = 2 K g,1 = K g,2 K g,3 =.1 d =.6 The following four models are positive-plus-negative feedback oscillators: 7. Meyer and Stryer model of calcium oscillations. This is a three ODE model that accounts for many of the observed characteristics of cytosolic calcium oscillations (S13). The variable x represents cytoplasmic Ca 2+ ; y represents IP 3 ; and z represents ER Ca 2+. The parameter 9
11 R, which represents the fractional activation of a calcium-mobilizing G-protein coupled receptor, is used to vary the oscillator s frequency. dx = c y 3 1 K 1 + y x 2 ( ) z " c 3 2 ( ) + c 2 3 z2 " c 6 (x /c 7 ) c 6 dy = c R x 4 " c 5 y x + K 3 dz = "c y 3 1 K 1 + y x + K 2 x 2 ( ) z + c 3 2 ( ) " c 2 3 z2 Initial conditions: x[] =.1 y[] =.1 z[] = 25 Parameters: c 1 = 6.64 c 2 = 5 c 3 =.313 c 4 c 5 = 2 c 6 =.5 c 7 =.6 K 1 =.1 K 2 =.15 K 3 x + K 2 8. van der Pol oscillator. The van der Pol oscillator model is a one-variable, second order differential equation, which can also be written as a pair of first order differential equations. The model has one parameter, φ, which was used to vary the frequency of oscillations. dx = y dy = "(1# x 2 )y # x 1
12 Initial conditions: x[] =.1 y[] = 9. Fitzhugh-Nagumo oscillator. Parameters for the Fitzhugh-Nagumo relaxation oscillator (S14, S15) were taken from the Vibe website (S16). The coupling parameter φ was varied to change the frequency of the oscillations. dv = v(v "#)(1" v) " w + $ dw = %(v " &w) Initial conditions: v[] = w[] = Parameters: " =.2 # =.112 $ = Cyanobacteria circadian oscillator. The cyanobacterial circadian oscillator has been reconstituted in vitro using three purified proteins, KaiA, KaiB, and KaiC (S17). One of these proteins, KaiC, undergoes a sequence of phosphorylation and dephosphorylation reactions that produces the circadian cycles. Rust and co-workers have modeled these oscillations as a set of three experimentally-constrained ODEs that describe the interconversion of KaiC between a threonine-phosphorylated form (T), a doubly-phosphorylated form (D), and a serinephosphorylated form (S) (S18). The remaining unphosphorylated form of KaiC can be calculated from these three forms and a conservation equation. dt = k UT (S)U + k DT (S)D " k TU (S)T " k TD (S)T dd = k TD(S)T + k SD (S)S " k DT (S)D " k DS (S)D 11
13 ds = k US (S)U + k DS (S)D " k SU (S)S " k SD (S)S where k XY (S) " k XY + k A XY A(S) K 1/ 2 + A(S) A(S) " max{, [KaiA] # 2S} Initial conditions: T[] =.68 D[].36 S[] =.34 Parameters: [KaiA].3 [KaiC] total = U + T + D + S = 3.4 K 1/ 2 =.43 k UT k TU k DS k SU = k TD =.21 =.31 =.11 = k SD = k US = k DT = k A UT =.479 k A TD =.213 k A SD =.56 k A US =.53 k A TU =.798 k A DT =.173 k A DS = " k A SU = #
14 Supplementary Figure S1: Fig. S1. Amplitude vs. frequency curves for randomly-parameterized oscillator models. Each curve corresponds to a different random parameter set. The red curves are from the negativefeedback-only model (Fig. 4A). The same curves are also included in Fig. 4F. The yellow and cyan curves are from the negative-plus-negative feedback model (Fig. 4C). The yellow curves correspond to weaker feedback (k 7 =-1) and the cyan curves to stronger feedback (k 7 =5-6). 13
15 Table S1: Ranges for the random parameter sets Parameter Units Range Distribution k 1 dimensionless 1 Uniform k 2 dimensionless 1 Uniform k 3 dimensionless 1 Uniform k 4 dimensionless 1 Uniform k 5 dimensionless 1 Fixed k 6 dimensionless 1 Uniform k 7 dimensionless Neg Fixed Pos+neg, Neg+neg 1 Uniform Stronger pos+neg, 5 6 Uniform Stronger neg+neg n 1 dimensionless 1 4 Uniform n 2 dimensionless 1 4 Uniform n 3 dimensionless 1 4 Uniform n 4 dimensionless 1 4 Uniform K 1 concentration 4 Uniform K 2 concentration 4 Uniform K 3 concentration 4 Uniform K 4 concentration 4 Uniform The parameters refer to the ordinary differential equation models shown in Figure 4A-C. The seven rate constants (k 1 through k 7 ) were non-dimensionalized by dividing by k 5. 14
16 References and Notes S1. A. Dhooge, W. Govaerts, Y. A. Kuznetwov, ACM Transactions on Mathematical Software 29, 141 (23). S2. J. R. Pomerening, S. Y. Kim, J. E. Ferrell, Jr., Cell 122, 565 (25). S3. J. R. Pomerening, E. D. Sontag, J. E. Ferrell, Jr., Nature Cell Biol. 5, 346 (23). S4. M. B. Elowitz, S. Leibler, Nature 43, 335 (2). S5. B. C. Goodwin, Ed., Oscillatory behavior in enzymatic control processes, vol. 3 (Permagon Press, Oxford UK, 1965), vol. 3, pp S6. P. Ruoff, S. Mohsenzadeh, L. Rensing, Naturwissenschaften 83, 514 (1996). S7. P. Ruoff, M. Vinsjevik, C. Monnerjahn, L. Rensing, J. Theoret. Biol. 29, 29 (21). S8. P. Ruoff, M. Vinsjevik, C. Monnerjahn, L. Rensing, J Biol. Rhythms 14, 469 (1999). S9. P. Ruoff, M. Vinsjevik, S. Mohsenzadeh, L. Rensing, J. Theoret. Biol. 196, 483 (1999). S1. P. Ruoff, L. Rensing, R. Kommedal, S. Mohsenzadeh, Chronobiol. Int. 14, 499 (1997). S11. O. A. Igoshin, A. Goldbeter, D. Kaiser, G. Oster, Proc. Natl. Acad. Sci. U.S.A. 11, 1576 (24). S12. E. Fung et al., Nature 435, 118 (25). S13. T. Meyer, L. Stryer, Proc. Natl. Acad. Sci. U.S.A. 85, 551 (1988). S14. R. FitzHugh, Biophys. J. 1, 445 (1961). S15. J. Nagumo, S. Arimoto, S. Yoshizawa, Proc. IRE 5, 261 (1964). S16. D. S. Eck. Vibe ( 23). S17. M. Nakajima et al., Science 38, 414 (25). S18. M. J. Rust, J. S. Markson, W. S. Lane, D. S. Fisher, E. K. O Shea, Science 318, 89 (27). 15
BE 150 Problem Set #4 Issued: 24 April 2014 Due: May
M. Elowitz and R. M. Murray Spring 2014 CALIFORNIA INSTITUTE OF TECHNOLOGY Biology and Biological Engineering (BBE) BE 150 Problem Set #4 Issued: 24 April 2014 Due: May 1 2014 1. (Repressilator modeling
More informationnutrients growth & division repellants movement
Network Dynamics and Cell Physiology John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute Outline 1. Cell Signaling: Physiology 2. Cell Signaling: Molecular Biology 3. Chemical
More informationCell-free extract characterization
Cell-free extract characterization Neha Bidthanpally, Advisor: Dr. Qiong Yang Department of Biophysics, University of Michigan Abstract A cell-free extract system is an in vitro tool used to study biological
More informationA synthetic oscillatory network of transcriptional regulators
A synthetic oscillatory network of transcriptional regulators Michael B. Elowitz & Stanislas Leibler, Nature, 403, 2000 igem Team Heidelberg 2008 Journal Club Andreas Kühne Introduction Networks of interacting
More informationModelling biological oscillations
Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van
More informationBioinformatics 3. V18 Kinetic Motifs. Fri, Jan 8, 2016
Bioinformatics 3 V18 Kinetic Motifs Fri, Jan 8, 2016 Modelling of Signalling Pathways Curr. Op. Cell Biol. 15 (2003) 221 1) How do the magnitudes of signal output and signal duration depend on the kinetic
More informationBioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014"
Bioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014" Modelling of Signalling Pathways" Curr. Op. Cell Biol. 15 (2003) 221" 1) How do the magnitudes of signal output and signal duration depend on the
More informationNetwork Dynamics and Cell Physiology. John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute
Network Dynamics and Cell Physiology John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute Signal Transduction Network Hanahan & Weinberg (2000) Gene Expression Signal-Response
More informationThe Cdk1-APC/C cell cycle oscillator circuit functions as a time-delayed, ultrasensitive switch
1 Supplementary Information for: The Ck1-APC/C cell cycle oscillator circuit functions as a time-elaye, ultrasensitive switch Qiong Yang*, an James E. Ferrell, Jr.* *Department of Chemical an Systems Biology
More informationA minimal mathematical model combining several regulatory cycles from the budding yeast cell cycle
A minimal mathematical model combining several regulatory cycles from the budding yeast cell cycle K. Sriram, G. Bernot and F. Képès Abstract: A novel topology of regulatory networks abstracted from the
More informationDETC DELAY DIFFERENTIAL EQUATIONS IN THE DYNAMICS OF GENE COPYING
Proceedings of the ASME 2007 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 2007 September 4-7, 2007, Las Vegas, Nevada, USA DETC2007-3424
More informationFUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models
FUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models Markus W. Covert Stanford University 0 CRC Press Taylor & Francis Group Boca Raton London New York Contents /... Preface, xi
More informationANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY. Nael H. El-Farra, Adiwinata Gani & Panagiotis D. Christofides
ANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY Nael H El-Farra, Adiwinata Gani & Panagiotis D Christofides Department of Chemical Engineering University of California, Los Angeles 2003 AIChE
More informationNetwork Dynamics and Cell Physiology. John J. Tyson Dept. Biological Sciences Virginia Tech
Network Dynamics and Cell hysiology John J. Tyson Dept. Biological Sciences Virginia Tech Collaborators Budapest Univ. Techn.. & Econ. Bela Novak Attila Csikasz-Nagy Andrea Ciliberto Virginia Tech Kathy
More informationLIMIT CYCLE OSCILLATORS
MCB 137 EXCITABLE & OSCILLATORY SYSTEMS WINTER 2008 LIMIT CYCLE OSCILLATORS The Fitzhugh-Nagumo Equations The best example of an excitable phenomenon is the firing of a nerve: according to the Hodgkin
More informationdamped oscillation Figure 11.3
Chpater X How to build a Biological oscillator Oscillations are thrilling hearts beat, cells divide every cell cycle, circadian clocks keep, neurons click in trains of regularly spaced spikes. Biological
More informationA Synthetic Oscillatory Network of Transcriptional Regulators
A Synthetic Oscillatory Network of Transcriptional Regulators Michael Elowitz & Stanislas Leibler Nature, 2000 Presented by Khaled A. Rahman Background Genetic Networks Gene X Operator Operator Gene Y
More informationdynamic processes in cells (a systems approach to biology)
dynamic processes in cells (a systems approach to biology) jeremy gunawardena department of systems biology harvard medical school lecture 12 11 october 2016 = vasopressin CCh = carbachol GnRH = gonadotropin
More informationBE 150 Problem Set #2 Issued: 16 Jan 2013 Due: 23 Jan 2013
M. Elowitz and R. M. Murray Winter 2013 CALIFORNIA INSTITUTE OF TECHNOLOGY Biology and Biological Engineering (BBE) BE 150 Problem Set #2 Issued: 16 Jan 2013 Due: 23 Jan 2013 1. (Shaping pulses; based
More informationL'horloge circadienne, le cycle cellulaire et leurs interactions. Madalena CHAVES
L'horloge circadienne, le cycle cellulaire et leurs interactions?? Madalena CHAVES What do all these living organisms have in common? Mus musculus Neurospora crassa Drosophila melanogaster Arabidopsis
More informationAnalysis and Simulation of Biological Systems
Analysis and Simulation of Biological Systems Dr. Carlo Cosentino School of Computer and Biomedical Engineering Department of Experimental and Clinical Medicine Università degli Studi Magna Graecia Catanzaro,
More informationWritten Exam 15 December Course name: Introduction to Systems Biology Course no
Technical University of Denmark Written Exam 15 December 2008 Course name: Introduction to Systems Biology Course no. 27041 Aids allowed: Open book exam Provide your answers and calculations on separate
More informationNIH Public Access Author Manuscript Nat Rev Mol Cell Biol. Author manuscript; available in PMC 2009 December 21.
NIH Public Access Author Manuscript Published in final edited form as: Nat Rev Mol Cell Biol. 2008 December ; 9(12): 981 991. doi:10.1038/nrm2530. Design Principles of Biochemical Oscillators Béla Novak
More informationIntroduction. Dagmar Iber Jörg Stelling. CSB Deterministic, SS 2015, 1.
Introduction Dagmar Iber Jörg Stelling joerg.stelling@bsse.ethz.ch CSB Deterministic, SS 2015, 1 Origins of Systems Biology On this assumption of the passage of blood, made as a basis for argument, and
More informationAlternating Oscillations and Chaos in a Model of Two Coupled Biochemical Oscillators Driving Successive Phases of the Cell Cycle
Alternating Oscillations and Chaos in a Model of Two Coupled Biochemical Oscillators Driving Successive Phases of the Cell Cycle PIERRE-CHARLES ROMOND, a MAURO RUSTICI, b DIDIER GONZE, AND ALBERT GOLDBETER
More informationSupplementary Materials for
www.advances.sciencemag.org/cgi/content/full/1/5/e1500358/dc1 Supplementary Materials for Transplantability of a circadian clock to a noncircadian organism Anna H. Chen, David Lubkowicz, Vivian Yeong,
More informationBioControl - Week 2, Lecture 1
BioControl - Week 2, Lecture Goals of this lecture: CDS tools for bio-molecular systems Study of equilibria Dynamic performance Periodic behaviors Main eample: Negative autoregulation Suggested readings
More informationBiophysical Journal Volume 92 May
Biophysical Journal Volume 92 May 2007 3407 3424 3407 Dynamics of a Minimal Model of Interlocked Positive and Negative Feedback Loops of Transcriptional Regulation by camp-response Element Binding Proteins
More informationSystems biology and complexity research
Systems biology and complexity research Peter Schuster Institut für Theoretische Chemie, Universität Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA Interdisciplinary Challenges for
More informationGene Regulatory Network Identification
Gene Regulatory Network Identification Korkut Uygun and Yinlun Huang Department of Chemical Engineering and Materials Science Wayne State University AIChE Annual National Meeting San Francisco November
More informationThe dynamics of hybrid metabolic-genetic oscillators
The dynamics of hybrid metabolic-genetic oscillators Ed Reznik, Tasso J. Kaper, and Daniel Segrè Citation: Chaos 3, 1313 (13) doi: 1.163/1.4793573 View online: http://dx.doi.org/1.163/1.4793573 View Table
More informationCHEM 515: Chemical Kinetics and Dynamics
Alejandro J. Garza S01163018 Department of Chemistry, Rice University, Houston, TX email: ajg7@rice.edu, ext. 2657 Submitted December 12, 2011 Abstract Spontaneous antispiral wave formation was observed
More informationProceedings of the FOSBE 2007 Stuttgart, Germany, September 9-12, 2007
Proceedings of the FOSBE 2007 Stuttgart, Germany, September 9-12, 2007 A FEEDBACK APPROACH TO BIFURCATION ANALYSIS IN BIOCHEMICAL NETWORKS WITH MANY PARAMETERS Steffen Waldherr 1 and Frank Allgöwer Institute
More informationResearch Article Robust Cell Size Checkpoint from Spatiotemporal Positive Feedback Loop in Fission Yeast
BioMed Research International Volume 213, Article ID 91941, 9 pages http://dx.doi.org/1.1155/213/91941 Research Article Robust Cell Size Checkpoint from Spatiotemporal Positive Feedback Loop in Fission
More informationEngineering of Repressilator
The Utilization of Monte Carlo Techniques to Simulate the Repressilator: A Cyclic Oscillatory System Seetal Erramilli 1,2 and Joel R. Stiles 1,3,4 1 Bioengineering and Bioinformatics Summer Institute,
More informationModelling Biochemical Reaction Networks. Lecture 17: Modeling the cell cycle, Part I
Modelling Biochemical Reaction Networks Lecture 17: Modeling the cell cycle, Part I Marc R. Roussel Department of Chemistry and Biochemistry Recommended reading Fall, Marland, Wagner and Tyson, sections
More informationBasic Synthetic Biology circuits
Basic Synthetic Biology circuits Note: these practices were obtained from the Computer Modelling Practicals lecture by Vincent Rouilly and Geoff Baldwin at Imperial College s course of Introduction to
More informationNo oscillations in the Michaelis-Menten approximation of the dual futile cycle under a sequential and distributive mechanism
No oscillations in the Michaelis-Menten approximation of the dual futile cycle under a sequential and distributive mechanism Luna Bozeman 1 and Adriana Morales 2 1 Department of Mathematical Sciences Clemson
More informationBiomolecular Feedback Systems
Biomolecular Feedback Systems Domitilla Del Vecchio MIT Richard M. Murray Caltech Version 1.0b, September 14, 2014 c 2014 by Princeton University Press All rights reserved. This is the electronic edition
More informationThree different fusions led to three basic ideas: 1) If one fuses a cell in mitosis with a cell in any other stage of the cell cycle, the chromosomes
Section Notes The cell division cycle presents an interesting system to study because growth and division must be carefully coordinated. For many cells it is important that it reaches the correct size
More informationThe role of mrna and protein stability in the function of coupled positive and negative feedback systems in eukaryotic cells.
The role of mrna and protein stability in the function of coupled positive and negative feedback systems in eukaryotic cells. Kristian Moss Bendtsen 1, Mogens H. Jensen 1a, Sandeep Krishna 1,2, Szabolcs
More informationA Model of Periodic Oscillation for Genetic Regulatory Systems
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: FUNDAMENTAL THEORY AND APPLICATIONS, VOL. 49, NO. 10, OCTOBER 2002 1429 A Model of Periodic Oscillation for Genetic Regulatory Systems Luonan Chen, Senior Member,
More informationA simple model for the eukaryotic cell cycle. Andrea Ciliberto
A simple model for the eukaryotic cell cycle Andrea Ciliberto The cell division cycle G1 cell division Start S (DNA Replication) Finish M (mitosis) G2/M G2 Kohn, Mol. Biol. Cell., 1999 How did we get to
More informationChapter 15 Active Reading Guide Regulation of Gene Expression
Name: AP Biology Mr. Croft Chapter 15 Active Reading Guide Regulation of Gene Expression The overview for Chapter 15 introduces the idea that while all cells of an organism have all genes in the genome,
More informationOscillations in Michaelis-Menten Systems
Oscillations in Michaelis-Menten Systems Hwai-Ray Tung Brown University July 18, 2017 Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, 2017 1 / 22 Background A simple
More informationNetwork Dynamics and Cell Physiology
Network Dynamics and Cell Physiology John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute Outline 1. The cell is an information processing system. 2. Information is processed
More informationarxiv: v4 [q-bio.mn] 5 Dec 2017
Describing Function-based Approximations of Biomolecular Systems Abhishek Dey 1, Shaunak Sen 2 Department of Electrical Engineering Indian Institute of Technology Delhi Hauz Khas, New Delhi 110016, INDIA
More informationCS-E5880 Modeling biological networks Gene regulatory networks
CS-E5880 Modeling biological networks Gene regulatory networks Jukka Intosalmi (based on slides by Harri Lähdesmäki) Department of Computer Science Aalto University January 12, 2018 Outline Modeling gene
More informationDNA replication. M (mitosis)
Introduction to cell cycle modeling Attila Csikász-Nagy G1 The cell cycle is the sequence of events whereby a growing cell replicates all its components and divides them more-or-less evenly between two
More informationURL: <
Citation: ngelova, Maia and en Halim, sma () Dynamic model of gene regulation for the lac operon. Journal of Physics: Conference Series, 86 (). ISSN 7-696 Published by: IOP Publishing URL: http://dx.doi.org/.88/7-696/86//7
More informationIntroduction to Biological Modeling
Introduction to Biological odeling Lecture : odeling dynamics Sept. 9, 010 Steve Andrews Brent lab, Basic Sciences Division, FHCRC Last week Why model biology? Example: E. coli chemotaxis Typical modeling
More informationSupporting Information. Methods. Equations for four regimes
Supporting Information A Methods All analytical expressions were obtained starting from quation 3, the tqssa approximation of the cycle, the derivation of which is discussed in Appendix C. The full mass
More informationSupplementary Material for Load Capacity. Improvements in Nucleic Acid Based Systems. Using Partially Open Feedback Control
Supplementary Material for Load Capacity Improvements in Nucleic Acid Based Systems Using Partially Open Feedback Control Vishwesh Kulkarni,, Evgeny Kharisov, Naira Hovakimyan, and Jongmin Kim Institute
More informationDynamical Systems in Biology
Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems
More information7. Well-Stirred Reactors IV
7. Well-Stirred Reactors IV The Belousov-Zhabotinsky Reaction: Models and Experiments Oregonator [based on the FKN mechanism; Field, R. J. & Noyes, R. M.: Oscillations in chemical systems. IV. Limit cycle
More informationStochastic dynamics of small gene regulation networks. Lev Tsimring BioCircuits Institute University of California, San Diego
Stochastic dynamics of small gene regulation networks Lev Tsimring BioCircuits Institute University of California, San Diego Nizhni Novgorod, June, 2011 Central dogma Activator Gene mrna Protein Repressor
More informationIntroduction. ECE/CS/BioEn 6760 Modeling and Analysis of Biological Networks. Adventures in Synthetic Biology. Synthetic Biology.
Introduction ECE/CS/BioEn 6760 Modeling and Analysis of Biological Networks Chris J. Myers Lecture 17: Genetic Circuit Design Electrical engineering principles can be applied to understand the behavior
More informationNetworks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationIntroduction to Mathematical Physiology I - Biochemical Reactions
Introduction to Mathematical Physiology I - Biochemical Reactions J. P. Keener Mathematics Department Math Physiology p.1/28 Introduction The Dilemma of Modern Biology The amount of data being collected
More informationStochastic Models for Circadian Oscillations: Emergence of a Biological Rhythm
Stochastic Models for Circadian Oscillations: Emergence of a Biological Rhythm DIDIER GONZE, JOSÉ HALLOY, ALBERT GOLDBETER Unité de Chronobiologie Théorique, Faculté des Sciences, Université Libre de Bruxelles,
More informationBiophysical Chemistry 72 (1998)
Biophysical Chemistry 72 (1998) 153 167 Link between fertilization-induced Ca 2+ oscillations and relief from metaphase II arrest in mammalian eggs: a model based on calmodulin-dependent kinase II activation
More information56:198:582 Biological Networks Lecture 11
56:198:582 Biological Networks Lecture 11 Network Motifs in Signal Transduction Networks Signal transduction networks Signal transduction networks are composed of interactions between signaling proteins.
More informationIdentification and Modeling of Co-Rhythmic Genes from Micro-array Time Series Data
Proceedings of the 7th World Congress The International Federation of Automatic Control Seoul, Korea, July 6-, 8 Identification and Modeling of Co-Rhythmic Genes from Micro-array Time Series Data Wenxue
More informationTuning an Activator-Repressor Clock Employing Retroactivity
Tuning an Activator-Repressor Clock Employing Retroactivity Alexander Rosenberg, Shridhar Jayanthi and Domitilla Del Vecchio Electrical Engineering Dept., University of Washington, Seattle WA 9895 Electrical
More informationBIOMED Programming & Modeling for BME Final Exam, , Instructor: Ahmet Sacan
BIOMED 201 - Programming & Modeling for BME Final Exam, 2011.12.01, Instructor: Ahmet Sacan Sign the honor code below. No credit will be given for the exam without a signed pledge. I have neither given
More informationIntroduction to Bioinformatics
Systems biology Introduction to Bioinformatics Systems biology: modeling biological p Study of whole biological systems p Wholeness : Organization of dynamic interactions Different behaviour of the individual
More informationPlant Molecular and Cellular Biology Lecture 8: Mechanisms of Cell Cycle Control and DNA Synthesis Gary Peter
Plant Molecular and Cellular Biology Lecture 8: Mechanisms of Cell Cycle Control and DNA Synthesis Gary Peter 9/10/2008 1 Learning Objectives Explain why a cell cycle was selected for during evolution
More informationComputational Methods for Nonlinear Systems
Computational Methods for Nonlinear Systems Cornell Physics 682 / CIS 629 James P. Sethna Christopher R. Myers Computational Methods for Nonlinear Systems Graduate computational science laboratory course
More information2. Mathematical descriptions. (i) the master equation (ii) Langevin theory. 3. Single cell measurements
1. Why stochastic?. Mathematical descriptions (i) the master equation (ii) Langevin theory 3. Single cell measurements 4. Consequences Any chemical reaction is stochastic. k P d φ dp dt = k d P deterministic
More informationQuantitative Analysis of Cellular Metabolic Dissipative, Self-Organized Structures
Int. J. Mol. Sci. 2010, 11, 3540-3599; doi:10.3390/ijms11093540 Review OPEN ACCESS International Journal of Molecular Sciences ISSN 1422-0067 www.mdpi.com/journal/ijms Quantitative Analysis of Cellular
More informationHopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models
Vol.37 (SUComS 6), pp.57-66 http://dx.doi.org/.457/astl.6.37. Hopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models Shuang Guo, Xiuli Li, Jian Yu, Xiu Ren Department of Mathematics, Teacher
More informationCircadian rhythms and molecular noise
CHAOS 16, 026110 2006 Circadian rhythms and molecular noise Didier Gonze and Albert Goldbeter Faculté des Sciences, Université Libre de Bruxelles, Campus Plaine, C.P. 231, B-1050 Brussels, Belgium Received
More informationNotes, Mathematical Cell Biology Course. Leah Edelstein-Keshet
Notes, Mathematical Cell Biology Course Leah Edelstein-Keshet May 2, 2012 ii Leah Edelstein-Keshet Contents 1 Qualitative behaviour of simple ODEs and bifurcations 1 1.0.1 Cubic kinetics.......................
More informationTracking the State of the Hindmarsh-Rose Neuron by Using the Coullet Chaotic System Based on a Single Input
ISSN 1746-7659, England, UK Journal of Information and Computing Science Vol. 11, No., 016, pp.083-09 Tracking the State of the Hindmarsh-Rose Neuron by Using the Coullet Chaotic System Based on a Single
More informationEfficient Leaping Methods for Stochastic Chemical Systems
Efficient Leaping Methods for Stochastic Chemical Systems Ioana Cipcigan Muruhan Rathinam November 18, 28 Abstract. Well stirred chemical reaction systems which involve small numbers of molecules for some
More informationFrom quiescence to proliferation: Cdk oscillations drive the mammalian cell cycle
ORIGINAL RESEARCH ARTICLE published: 02 November 2012 doi: 10.3389/fphys.2012.00413 From quiescence to proliferation: Cdk oscillations drive the mammalian cell cycle Claude Gérard and Albert Goldbeter*
More informationBiomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems
Biomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems Processing information to coordinate activity regulatory & signaling networks
More informationAnalysis of Mode Transitions in Biological Networks
Analysis of Mode Transitions in Biological Networks Nael H. El-Farra, Adiwinata Gani, and Panagiotis D. Christofides Dept. of Chemical Engineering, University of California, Los Angeles, CA 90095 DOI 10.1002/aic.10499
More informationA(γ A D A + γ R D R + γ C R + δ A )
Title: Mechanisms of noise-resistance in genetic oscillators Authors: José M. G. Vilar 1,2, Hao Yuan Kueh 1, Naama Barkai 3, and Stanislas Leibler 1,2 1 Howard Hughes Medical Institute, Departments of
More information1 Periodic stimulations of the Incoherent Feedforward Loop network
1 Periodic stimulations of the Incoherent Feedforward Loop network In this Additional file, we give more details about the mathematical analysis of the periodic activation of the IFFL network by a train
More informationRegulation of metabolism
Regulation of metabolism So far in this course we have assumed that the metabolic system is in steady state For the rest of the course, we will abandon this assumption, and look at techniques for analyzing
More informationDynamic modeling and analysis of cancer cellular network motifs
SUPPLEMENTARY MATERIAL 1: Dynamic modeling and analysis of cancer cellular network motifs Mathieu Cloutier 1 and Edwin Wang 1,2* 1. Computational Chemistry and Bioinformatics Group, Biotechnology Research
More information56:198:582 Biological Networks Lecture 9
56:198:582 Biological Networks Lecture 9 The Feed-Forward Loop Network Motif Subgraphs in random networks We have discussed the simplest network motif, self-regulation, a pattern with one node We now consider
More informationMultistability in the lactose utilization network of Escherichia coli
Multistability in the lactose utilization network of Escherichia coli Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J. Ruby Abrams Agenda Motivation Intro to multistability Purpose
More informationBioinformatics: Network Analysis
Bioinformatics: Network Analysis Kinetics of Gene Regulation COMP 572 (BIOS 572 / BIOE 564) - Fall 2013 Luay Nakhleh, Rice University 1 The simplest model of gene expression involves only two steps: the
More informationPrinciples of Synthetic Biology: Midterm Exam
Principles of Synthetic Biology: Midterm Exam October 28, 2010 1 Conceptual Simple Circuits 1.1 Consider the plots in figure 1. Identify all critical points with an x. Put a circle around the x for each
More informationModelling nucleocytoplasmic transport with application to the intracellular dynamics of the tumor suppressor protein p53
Modelling nucleocytoplasmic transport with application to the intracellular dynamics of the tumor suppressor protein p53 L. Dimitrio Supervisors: R. Natalini (IAC-CNR) and J. Clairambault (INRIA) 5 Septembre
More informationBifurcation Analysis of a Model of Mitotic Control in Frog Eggs
J. theor. Biol. (1998) 195, 69 85 Article No. jt980781 Bifurcation Analysis of a Model of Mitotic Control in Frog Eggs MARK T. BORISUK AND JOHN J. TYSON* Department of Biology, Virginia olytechnic Institute
More informationRegulation and signaling. Overview. Control of gene expression. Cells need to regulate the amounts of different proteins they express, depending on
Regulation and signaling Overview Cells need to regulate the amounts of different proteins they express, depending on cell development (skin vs liver cell) cell stage environmental conditions (food, temperature,
More informationModeling Multiple Steady States in Genetic Regulatory Networks. Khang Tran. problem.
Modeling Multiple Steady States in Genetic Regulatory Networks Khang Tran From networks of simple regulatory elements, scientists have shown some simple circuits such as the circadian oscillator 1 or the
More informationGene Network Science Diagrammatic Cell Language and Visual Cell
Gene Network Science Diagrammatic Cell Language and Visual Cell Mr. Tan Chee Meng Scientific Programmer, System Biology Group, Bioinformatics Institute Overview Introduction Why? Challenges Diagrammatic
More informationA Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationLife Sciences 1a: Section 3B. The cell division cycle Objectives Understand the challenges to producing genetically identical daughter cells
Life Sciences 1a: Section 3B. The cell division cycle Objectives Understand the challenges to producing genetically identical daughter cells Understand how a simple biochemical oscillator can drive the
More informationRegulation of gene expression. Premedical - Biology
Regulation of gene expression Premedical - Biology Regulation of gene expression in prokaryotic cell Operon units system of negative feedback positive and negative regulation in eukaryotic cell - at any
More information56:198:582 Biological Networks Lecture 8
56:198:582 Biological Networks Lecture 8 Course organization Two complementary approaches to modeling and understanding biological networks Constraint-based modeling (Palsson) System-wide Metabolism Steady-state
More informationStochastic Processes around Central Dogma
Stochastic Processes around Central Dogma Hao Ge haoge@pu.edu.cn Beijing International Center for Mathematical Research Biodynamic Optical Imaging Center Peing University, China http://www.bicmr.org/personal/gehao/
More informationLecture 2: Analysis of Biomolecular Circuits
Lecture 2: Analysis of Biomolecular Circuits Richard M. Murray Caltech CDS/BE Goals: Give a short overview of the control techniques applied to biology - uncertainty management - system identification
More information26 Robustness in a model for calcium signal transduction dynamics
26 Robustness in a model for calcium signal transduction dynamics U. Kummer 1,G.Baier 2 and L.F. Olsen 3 1 European Media Laboratory, Villa Bosch, Schloss-Wolfsbrunnenweg 33, 69118 Heidelberg, Germany
More informationSimulating Dark Expressions and Interactions of frq and wc-1 in the Neurospora Circadian Clock
Biophysical Journal Volume 94 February 2008 1221 1232 1221 Simulating Dark Expressions and Interactions of frq and wc-1 in the Neurospora Circadian Clock Christian I. Hong,* y Ingunn W. Jolma,* Jennifer
More informationModeling Cellular Networks
06_4774 12/7/06 3:02 PM Page 151 CHAPTER 6 Modeling Cellular Networks Tae Jun Lee, Dennis Tu, Chee Meng Tan, and Lingchong You 6.1 Introduction Systems-level understanding of cellular dynamics is important
More information