Oscillations in Michaelis-Menten Systems
|
|
- Virgil Cummings
- 5 years ago
- Views:
Transcription
1 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, / 22
2 Background A simple example of genetic oscillation Gene makes compound (transcription factor) Transcription factor binds to promoter Circadian Clocks Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
3 Background A simple example of genetic oscillation Gene makes compound (transcription factor) Transcription factor binds to promoter Circadian Clocks Mass Action Kinetics Write chemical reaction network as a system of differential equations Michaelis-Menten (MM) Approximation for Biochemical Systems Assume low concentration of intermediates Allows for elimination of variables, reducing differential equations from mass action Michaelis-Menten approximation has oscillations implies original has oscillations Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
4 Michaelis-Menten System (Dual Futile Cycle) Phosphorylation Dephosphorylation S 0 + E k 1 S 0 E k 3 S 1 + E k 4 S 1 E k 6 S 2 + E k2 k5 S 2 + F l 1 S 2 F l 3 S 1 + F l 4 S 1 F l 6 S 0 + F l2 l5 Oscillatory behavior is unknown. Wang and Sontag (2008) showed Michaelis-Menten approximation has no oscillations Bozeman and Morales (REU 2016) showed Michaelis-Menten approximation has no oscillations with more elementary techniques Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
5 Processive vs Distributive Distributive: S 0 + E k 1 S 0 E k 3 S 1 + E k 4 S 1 E k 6 S 2 + E k2 k5 S 2 + F l 1 S 2 F l 3 S 1 + F l 4 S 1 F l 6 S 0 + F l2 l5 Processive: S 0 + E k 1 S 0 E k 7 S 1 E k 6 S 2 + E k2 S 2 + F l 1 S 2 F l 7 S 1 F l 6 S 0 + F l2 Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
6 This Year s Project S 1 + F S 2 + F S 2 F S 1 F S 0 + F S 0 + E S 0 E S 1 E S 2 + E S 1 + E Does it have oscillations? Does its Michaelis-Menten approximation have oscillations? Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
7 Mass Action Kinetics Write chemical reaction network as a system of differential equations Example for rate of [S 0 E] S 0 + E k 1 S 0 E gives Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
8 Mass Action Kinetics Write chemical reaction network as a system of differential equations Example for rate of [S 0 E] S 0 + E k 1 S 0 E gives d[s 0 E] = k 1 [S 0 ][E] Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
9 Mass Action Kinetics Write chemical reaction network as a system of differential equations Example for rate of [S 0 E] S 0 + E k 1 S 0 E gives S 0 + E k 2 S 0 E gives d[s 0 E] = k 1 [S 0 ][E] Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
10 Mass Action Kinetics Write chemical reaction network as a system of differential equations Example for rate of [S 0 E] S 0 + E k 1 S 0 E gives S 0 + E k 2 S 0 E gives d[s 0 E] d[s 0 E] = k 1 [S 0 ][E] = k 2 [S 0 E] Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
11 Mass Action Kinetics Write chemical reaction network as a system of differential equations Example for rate of [S 0 E] S 0 + E k 1 S 0 E gives S 0 + E k 2 S 0 E gives d[s 0 E] = k 1 [S 0 ][E] S 0 + E k 1 k2 S 0 E gives d[s 0 E] = k 2 [S 0 E] Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
12 Mass Action Kinetics Write chemical reaction network as a system of differential equations Example for rate of [S 0 E] S 0 + E k 1 S 0 E gives S 0 + E k 2 S 0 E gives d[s 0 E] = k 1 [S 0 ][E] S 0 + E k 1 k2 S 0 E gives d[s 0 E] = k 2 [S 0 E] d[s 0 E] = k 1 [S 0 ][E] k 2 [S 0 E] Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
13 Corresponding ODEs d[s 0] = l 6[S 1F ] k 1[S 0][E] + k 2[S 0E] d[s 2] = k 6[S 1E] l 1[S 2][F ] + l 2[S 2F ] d[s 1] = k 3[S 0E] k 4[S 1][E] + k 5[S 1E] + l 3[S 2F ] + l 5[S 1F ] l 4[S 1][F ] d[e] = (k 2 + k 3)[S 0E] + (k 5 + k 6)[S 1E] k 1[S 0][E] k 4[S 1][E] d[f ] = (l 2 + l 3)[S 2F ] + (l 5 + l 6)[S 1F ] l 1[S 2][F ] l 4[S 1][F ] d[s 0E] = k 1[S 0][E] (k 2 + k 3)[S 0E] k 7[S 0E] d[s 1E] = k 4[S 1][E] (k 5 + k 6)[S 1E]+k 7[S 0E] d[s 2F ] = l 1[S 2][F ] (l 2 + l 3)[S 2F ] l 7[S 2F ] d[s 1F ] = l 4[S 1][F ] (l 5 + l 6)[S 1F ]+l 7[S 2F ] Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
14 Michaelis-Menten Approximation of ODEs: Conservation Equations Conservation equations refer to conservation of mass. S T = [S 0 ] + [S 1 ] + [S 2 ] + [S 0 E] + [S 1 E] + [S 2 F ] + [S 1 F ] E T = [E] + [S 0 E] + [S 1 E] F T = [F ] + [S 2 F ] + [S 1 F ] S 1 + F S 2 + F S 2 F S 1 F S 0 + F S 0 + E S 0 E S 1 E S 2 + E S 1 + E Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
15 Michaelis-Menten Approximation of ODEs: Conservation Equations Conservation equations refer to conservation of mass. S T = [S 0 ] + [S 1 ] + [S 2 ] + [S 0 E] + [S 1 E] + [S 2 F ] + [S 1 F ] E T = [E] + [S 0 E] + [S 1 E] F T = [F ] + [S 2 F ] + [S 1 F ] d[e] d[s 0 E] d[s 1 E] = (k 2 + k 3 )[S 0 E] + (k 5 + k 6 )[S 1 E] k 1 [S 0 ][E] k 4 [S 1 ][E] = k 1 [S 0 ][E] (k 2 + k 3 )[S 0 E] k 7 [S 0 E] = k 4 [S 1 ][E] (k 5 + k 6 )[S 1 E]+k 7 [S 0 E] Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
16 Michaelis-Menten Approximation of ODEs: Assumption We assume enzyme concentrations are small: E T = εẽt, F T = ε F T, [E] = ε[ẽ], [F ] = ε[ F ], [S 0 E] = ε[ S 0 E], [S 1 E] = ε[ S 1 E], [S 2 F ] = ε[ S 2 F ], [S 1 F ] = ε[ S 1 F ], τ = εt Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
17 Michaelis-Menten Approximation of ODEs: Plug in d[s 0] = l 6[ S 1F ] k 1[S 0][Ẽ] dτ + k2[ S 0E] d[s 2] = k 6[ S 1E] l 1[S 2][ dτ F ] + l 2[ S 2F ] X d[s1] = k 3[S 0E] k 4[S 1][E] + k 5[S 1E] + l 3[S 2F ] + l 5[S 1F ] l 4[S 1][F ] X d[e] = (k 2 + k 3)[S 0E] + (k 5 + k 6)[S 1E] k 1[S 0][E] k 4[S 1][E] X d[f ] = (l 2 + l 3)[S 2F ] + (l 5 + l 6)[S 1F ] l 1[S 2][F ] l 4[S 1][F ] ε d[ S 0E] dτ ε d[ S 1E] dτ ε d[ S 2F ] dτ ε d[ S 1F ] dτ = k 1[S 0][Ẽ] (k2 + k3)[ S 0E] k 7[ S 0E] = k 4[S 1][Ẽ] (k5 + k6)[ S 1E]+k 7[ S 0E] = l 1[S 2][ F ] (l 2 + l 3)[ S 2F ] l 7[ S 2F ] = l 4[S 1][ F ] (l 5 + l 6)[ S 1F ]+l 7[ S 2F ] Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
18 Michaelis-Menten Approximation of ODEs d[s 0 ] dτ d[s 2 ] dτ = (a 1[S 1 ] + a 2 [S 2 ])[ F T ] 1 + c 2 [S 1 ] + d 2 [S 2 ] = (a 4[S 1 ] + a 5 [S 0 ])[ẼT ] 1 + b 1 [S 0 ] + c 1 [S 1 ] [S 1 ] = S T [S 0 ] [S 2 ] a 3 [S 0 ][ẼT ] 1 + b 1 [S 0 ] + c 1 [S 1 ] a 6 [S 2 ][ F T ] 1 + c 2 [S 1 ] + d 2 [S 2 ] b 1 = k 1(k 5 + k 6 + k 7) (k 2 + k 3 + k 7)(k 5 + k 6) k4 c 1 = k 5 + k 6 l4l6 a 1 = l 5 + l 6 l 1l 6l 7 a 2= (l 2 + l 3 + l 7)(l 5 + l 6) k1(k3 + k7) a 3 = (k 2 + k 3 + k 7) d 1 = l 1(l 5 + l 6 + l 7) (l 2 + l 3 + l 7)(l 5 + l 6) l4 c 2 = l 5 + l 6 k4k6 a 4 = k 5 + k 6 k 1k 6k 7 a 5= (k 2 + k 3 + k 7)(k 5 + k 6) l1(l3 + l7) a 6 = (l 2 + l 3 + l 7) Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
19 Dulac s Criterion f ([S 0 ], [S 2 ]) = d[s 0] dτ g([s 0 ], [S 2 ]) = d[s 2] dτ = (a 1[S 1 ] + a 2 [S 2 ])[ F T ] 1 + c 2 [S 1 ] + d 2 [S 2 ] = (a 4[S 1 ] + a 5 [S 0 ])[ẼT ] 1 + b 1 [S 0 ] + c 1 [S 1 ] a 3 [S 0 ][ẼT ] 1 + b 1 [S 0 ] + c 1 [S 1 ] a 6 [S 2 ][ F T ] 1 + c 2 [S 1 ] + d 2 [S 2 ] Theorem (Dulac) If the sign of df d[s 0 ] + dg d[s 2 ] does not change across a simply connected domain in R 2, the system does not exhibit oscillations in the domain. Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
20 Results Theorem (T.) In order for the reduced system to exhibit oscillations, one of the following must be true: (a 5 c 1 a 4 b 1 a 3 c 1 )S T a 3 + a 4 (a 2 c 2 a 1 d 1 a 6 c 2 )S T a 1 + a 6 Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
21 Results Theorem (T.) In order for the reduced system to exhibit oscillations, one of the following must be true: (a 5 c 1 a 4 b 1 a 3 c 1 )S T a 3 + a 4 Corollary (a 2 c 2 a 1 d 1 a 6 c 2 )S T a 1 + a 6 The MM approximation of the distributive model does not have oscillations, as shown earlier by Bozeman & Morales (2016) and Wang & Sontag (2008). Corollary The MM approximation of the processive model does not have oscillations, as immediately implied from work by Conradi et al. (2005) and Conradi & Shiu (2015). Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
22 Results Theorem (T.) In order for the reduced system to exhibit oscillations, one of the following must be true: (a 5 c 1 a 4 b 1 a 3 c 1 )S T a 3 + a 4 Corollary (a 2 c 2 a 1 d 1 a 6 c 2 )S T a 1 + a 6 The MM approximation of the distributive model does not have oscillations, as shown earlier by Bozeman & Morales (2016) and Wang & Sontag (2008). Corollary The MM approximation of the processive model does not have oscillations, as immediately implied from work by Conradi et al. (2005) and Conradi & Shiu (2015). Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
23 Results Theorem (T.) In order for the reduced system to exhibit oscillations, one of the following must be true: (a 5 c 1 a 4 b 1 a 3 c 1 )S T a 3 + a 4 Corollary (a 2 c 2 a 1 d 1 a 6 c 2 )S T a 1 + a 6 The MM approximation of the distributive model does not have oscillations, as shown earlier by Bozeman & Morales (2016) and Wang & Sontag (2008). Corollary The MM approximation of the processive model does not have oscillations, as immediately implied from work by Conradi et al. (2005) and Conradi & Shiu (2015). Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
24 Results Theorem (T.) In order for the reduced system to exhibit oscillations, one of the following must be true: (a 5 c 1 a 4 b 1 a 3 c 1 )S T a 3 + a 4 (a 2 c 2 a 1 d 1 a 6 c 2 )S T a 1 + a 6 Corollary The MM approximation of the mixed-mechanism model does not have oscillations. S 0 + E k 1 S 0 E k 3 S 1 + E k 4 S 1 E k 6 S 2 + E k2 k5 S 2 + F l 1 S 2 F l 7 S 1 F l 6 S 0 + F l2 Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
25 Future Directions: Does MM approx ever oscillate? My guess: Rarely, if ever. Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
26 Future Directions: Does MM approx ever oscillate? My guess: Rarely, if ever. Substitute each parameter with integer between 0 and 10. Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
27 Future Directions: Does MM approx ever oscillate? My guess: Rarely, if ever. Substitute each parameter with integer between 0 and 10. About 95% cannot oscillate by Dulac Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
28 Future Directions: Does MM approx ever oscillate? My guess: Rarely, if ever. Substitute each parameter with integer between 0 and 10. About 95% cannot oscillate by Dulac Of remaining 5%, about 99% reflect pattern below Blue Red Yellow Green Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
29 Future Directions: Does MM approx ever oscillate? My guess: Rarely, if ever. Substitute each parameter with integer between 0 and 10. About 95% cannot oscillate by Dulac Of remaining 5%, about 99% reflect pattern below Blue Red Yellow Green Possible proof from Wang & Sontag s approach (monotone systems theory) Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
30 Future Directions: Does MM approx ever oscillate? Hopf Bifurcations: Equilibrium changes stability type when parameters are changed Implies existence of oscillations generally Routh-Hurwitz Criterion gives necessary conditions for Hopf bifurcation Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
31 Future Directions: Does MM approx ever oscillate? Ex: Selkov model dx = x + ay + x 2 y dy = b ay x 2 y Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
32 Future Directions: Issues with MM Approximation Networks Can Oscillate MM Approx Can Oscillate Distributive Unknown No Processive No No Mixed-Mechanism Yes No Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
33 Future Directions: Issues with MM Approximation Conservation law S T S 0 + S 2 is violated. Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
34 Summary Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
35 Summary Wanted to examine oscillations in systems with processive and distributive elements Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
36 Summary Wanted to examine oscillations in systems with processive and distributive elements Applied Michaelis-Menten approximation to reduce number variables Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
37 Summary Wanted to examine oscillations in systems with processive and distributive elements Applied Michaelis-Menten approximation to reduce number variables Obtained a necessary condition for oscillations Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
38 Summary Wanted to examine oscillations in systems with processive and distributive elements Applied Michaelis-Menten approximation to reduce number variables Obtained a necessary condition for oscillations Discovered mixed-mechanism network has no oscillations, differing from original Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
39 Summary Wanted to examine oscillations in systems with processive and distributive elements Applied Michaelis-Menten approximation to reduce number variables Obtained a necessary condition for oscillations Discovered mixed-mechanism network has no oscillations, differing from original Future Directions Monotone Systems Theory if you think no oscillations Hopf Bifurcations if you think there are oscillations When does Michaelis-Menten Approximation preserve oscillations? Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
40 Thanks for Listening! Special thanks to Advisor: Dr. Anne Shiu Mentors: Nida Obatake, Ola Sobieska, Jonathan Tyler Host: Texas A&M University Funding: National Science Foundation Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, / 22
No 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 informationOscillations in Michaelis-Menten Systems
Oscillations in Michaelis-Menten Systems Hwai-Ray Tung July 20, 2017 Abstract Oscillations play a major role in a number of biological systems, the most notable example in biochemistry being circadian
More informationNetwork Analysis of Biochemical Reactions in Complex Environments
1 Introduction 1 Network Analysis of Biochemical Reactions in Complex Environments Elias August 1 and Mauricio Barahona, Department of Bioengineering, Imperial College London, South Kensington Campus,
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 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 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 informationM. Dechesne R. Sepulchre
Systems and models in chronobiology A delay model for the circadian rhythm M. Dechesne R. Sepulchre Department of Electrical Engineering and Computer Science Monteore Institute University of Liège 24th
More informationMath 480 The Vector Space of Differentiable Functions
Math 480 The Vector Space of Differentiable Functions The vector space of differentiable functions. Let C (R) denote the set of all infinitely differentiable functions f : R R. Then C (R) is a vector space,
More informationHomework Exercises and Solutions for Mathematics Fall 2016
Homework Exercises and Solutions for Mathematics 677 - Fall 216 Remark: Solutions may include maple files or matlab files. Assignment 1: (due Sept. 27, 216 1. In the real world trimolecular reactions are
More informationPersistence and Stationary Distributions of Biochemical Reaction Networks
Persistence and Stationary Distributions of Biochemical Reaction Networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Discrete Models in Systems Biology SAMSI December
More informationAn All-Encompassing Global Convergence Result for Processive Multisite Phosphorylation Systems
An All-Encompassing Global Convergence Result for Processive Multisite Phosphorylation Systems Mitchell Eithun Department of Mathematical Sciences, Ripon College, Ripon, WI 54971 Anne Shiu Department of
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 informationFYTN05/TEK267 Chemical Forces and Self Assembly
FYTN05/TEK267 Chemical Forces and Self Assembly FYTN05/TEK267 Chemical Forces and Self Assembly 1 Michaelis-Menten (10.3,10.4) M-M formalism can be used in many contexts, e.g. gene regulation, protein
More informationMath 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA Molecular biology and Biochemical kinetics Molecular biology is one of
More informationSome special cases
Lecture Notes on Control Systems/D. Ghose/2012 87 11.3.1 Some special cases Routh table is easy to form in most cases, but there could be some cases when we need to do some extra work. Case 1: The first
More informationEnzyme reaction example of Catalysis, simplest form: E + P at end of reaction No consumption of E (ES): enzyme-substrate complex Intermediate
V 41 Enzyme Kinetics Enzyme reaction example of Catalysis, simplest form: k 1 E + S k -1 ES E at beginning and ES k 2 k -2 E + P at end of reaction No consumption of E (ES): enzyme-substrate complex Intermediate
More informationEvery Neural Code Can Be Realized by Convex Sets
Every Neural Code Can Be Realized by Convex Sets Megan K. Franke and Samuel Muthiah July 21, 2017 Abstract Place cells are neurons found in some mammals that fire based on the animal s location in their
More informationCore genetic module: The mixed feedback loop
PHYSICAL REVIEW E 72, 031908 2005 Core genetic module: The mixed feedback loop Paul François and Vincent Hakim Laboratoire de Physique Statistique, * CNRS-UMR 8550, Ecole Normale Supérieure, 24, rue Lhomond
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 informationLecture 10: Proportional, Integral and Derivative Actions
MCE441: Intr. Linear Control Systems Lecture 10: Proportional, Integral and Derivative Actions Stability Concepts BIBO Stability and The Routh-Hurwitz Criterion Dorf, Sections 6.1, 6.2, 7.6 Cleveland State
More informationOscillations and bistability in a model of ERK regulation
Oscillations and bistability in a model of ERK regulation Nida Obatake, Anne Shiu, Xiaoxian Tang, and Angélica Torres 2 arxiv:93.267v [q-bio.mn] 6 Mar 29 Department of Mathematics, Texas A&M University,
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 informationProgram for the rest of the course
Program for the rest of the course 16.4 Enzyme kinetics 17.4 Metabolic Control Analysis 19.4. Exercise session 5 23.4. Metabolic Control Analysis, cont. 24.4 Recap 27.4 Exercise session 6 etabolic Modelling
More informationSingular perturbation theory
Singular perturbation theory Marc R Roussel October 19, 2005 1 Introduction When we apply the steady-state approximation (SSA) in chemical kinetics, we typically argue that some of the intermediates are
More informationChem Lecture 4 Enzymes Part 2
Chem 452 - Lecture 4 Enzymes Part 2 Question of the Day: Is there some easy way to clock how many reactions one enzyme molecule is able to catalyze in an hour? Thermodynamics I think that enzymes are molecules
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 informationCHEM-UA 652: Thermodynamics and Kinetics
CHEM-UA 65: Thermodynamics and Kinetics Notes for Lecture I. THE COMPLEXITY OF MULTI-STEP CHEMICAL REACTIONS It should be clear by now that chemical kinetics is governed by the mathematics of systems of
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 informationEnzyme Kinetics. Jonathan Gent and Douglas Saucedo May 24, 2002
Enzyme Kinetics Jonathan Gent and Douglas Saucedo May 24, 2002 Abstract This paper consists of a mathematical derivation of the Michaelis-Menten equation, which models the rate of reaction of certain enzymatic
More informationProblem Set 2. 1 Competitive and uncompetitive inhibition (12 points) Systems Biology (7.32/7.81J/8.591J)
Problem Set 2 1 Competitive and uncompetitive inhibition (12 points) a. Reversible enzyme inhibitors can bind enzymes reversibly, and slowing down or halting enzymatic reactions. If an inhibitor occupies
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 informationENZYME KINETICS. Medical Biochemistry, Lecture 24
ENZYME KINETICS Medical Biochemistry, Lecture 24 Lecture 24, Outline Michaelis-Menten kinetics Interpretations and uses of the Michaelis- Menten equation Enzyme inhibitors: types and kinetics Enzyme Kinetics
More information( ) ( = ) = ( ) ( ) ( )
( ) Vρ C st s T t 0 wc Ti s T s Q s (8) K T ( s) Q ( s) + Ti ( s) (0) τs+ τs+ V ρ K and τ wc w T (s)g (s)q (s) + G (s)t(s) i G and G are transfer functions and independent of the inputs, Q and T i. Note
More informationA First Course on Kinetics and Reaction Engineering. Class 9 on Unit 9
A First Course on Kinetics and Reaction Engineering Class 9 on Unit 9 Part I - Chemical Reactions Part II - Chemical Reaction Kinetics Where We re Going A. Rate Expressions - 4. Reaction Rates and Temperature
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 informationCHAPTER 6 : LITERATURE REVIEW
CHAPTER 6 : LITERATURE REVIEW Chapter : LITERATURE REVIEW 77 M E A S U R I N G T H E E F F I C I E N C Y O F D E C I S I O N M A K I N G U N I T S A B S T R A C T A n o n l i n e a r ( n o n c o n v e
More informationP E R E N C O - C H R I S T M A S P A R T Y
L E T T I C E L E T T I C E I S A F A M I L Y R U N C O M P A N Y S P A N N I N G T W O G E N E R A T I O N S A N D T H R E E D E C A D E S. B A S E D I N L O N D O N, W E H A V E T H E P E R F E C T R
More informationAlgebraic methods for the study of biochemical reaction networks
Algebraic methods for the study of biochemical reaction networks Alicia Dickenstein Departamento de Matemática, FCEN, Universidad de Buenos Aires, and Instituto de Matemática Luis A. Santaló, UBA CONICET
More informationOutline. Control systems. Lecture-4 Stability. V. Sankaranarayanan. V. Sankaranarayanan Control system
Outline Control systems Lecture-4 Stability V. Sankaranarayanan Outline Outline 1 Outline Outline 1 2 Concept of Stability Zero State Response: The zero-state response is due to the input only; all the
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 informationSystems Biology Across Scales: A Personal View XIV. Intra-cellular systems IV: Signal-transduction and networks. Sitabhra Sinha IMSc Chennai
Systems Biology Across Scales: A Personal View XIV. Intra-cellular systems IV: Signal-transduction and networks Sitabhra Sinha IMSc Chennai Intra-cellular biochemical networks Metabolic networks Nodes:
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 informationReversible reactions
Reversible reactions A reversible enzymic reaction (e.g. the conversion of glucose to fructose, catalysed by glucose isomerase) may be represented by the following scheme where the reaction goes through
More informationElectromagnetic Modeling and Simulation
Electromagnetic Modeling and Simulation Erin Bela and Erik Hortsch Department of Mathematics Research Experiences for Undergraduates April 7, 2011 Bela and Hortsch (OSU) EM REU 2010 1 / 45 Maxwell s Equations
More informationLecture 3 : Bifurcation Analysis
Lecture 3 : Bifurcation Analysis D. Sumpter & S.C. Nicolis October - December 2008 D. Sumpter & S.C. Nicolis General settings 4 basic bifurcations (as long as there is only one unstable mode!) steady state
More informationFeedback-mediated oscillatory coexistence in the chemostat
Feedback-mediated oscillatory coexistence in the chemostat Patrick De Leenheer and Sergei S. Pilyugin Department of Mathematics, University of Florida deleenhe,pilyugin@math.ufl.edu 1 Introduction We study
More informationClassifying Strictly Weakly Integral Modular Categories of Dimension 16p
Classifying Strictly Weakly Integral Modular Categories of Dimension 16p Elena Amparo College of William and Mary July 18, 2017 Elena Amparo (College of William and Mary)Classifying Strictly Weakly Integral
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 informationChapter 14: Chemical Kinetics II. Chem 102 Dr. Eloranta
Chapter 14: Chemical Kinetics II Chem 102 Dr. Eloranta Rate Laws If you are familiar with calculus Experiments would allow you to determine the reaction order and rate constant, but what if you wanted
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 informationMITOCW enzyme_kinetics
MITOCW enzyme_kinetics In beer and wine production, enzymes in yeast aid the conversion of sugar into ethanol. Enzymes are used in cheese-making to degrade proteins in milk, changing their solubility,
More informationHopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators
The Open Acoustics Journal 8 9-3 9 Open Access Hopf ifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators Jianping Cai *a and Jianhe Shen b a Department of
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 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 informationPrevious Class. Today. Cosubstrates (cofactors)
Previous Class Cosubstrates (cofactors) Today Proximity effect Basic equations of Kinetics Steady state kinetics Michaelis Menten equations and parameters Enzyme Kinetics Enzyme kinetics implies characterizing
More informationC a h p a t p e t r e r 6 E z n y z m y e m s
Chapter 6 Enzymes 4. Examples of enzymatic reactions acid-base catalysis: give and take protons covalent catalysis: a transient covalent bond is formed between the enzyme and the substrate metal ion catalysis:
More informationA kinetic mechanism inducing oscillations in simple chemical reactions networks
A kinetic mechanism inducing oscillations in simple chemical reactions networks J. Coatléven and C. Altafini September 1, 9 Abstract It is known that a kinetic reaction network in which one or more secondary
More informationChemical kinetics and catalysis
Chemical kinetics and catalysis Outline Classification of chemical reactions Definition of chemical kinetics Rate of chemical reaction The law of chemical raction rate Collision theory of reactions, transition
More information2 Dilution of Proteins Due to Cell Growth
Problem Set 1 1 Transcription and Translation Consider the following set of reactions describing the process of maing a protein out of a gene: G β g G + M M α m M β m M + X X + S 1+ 1 X S 2+ X S X S 2
More informationSign Conditions for Injectivity of Generalized Polynomial Maps with Applications to Chemical Reaction Networks and Real Algebraic Geometry
Sign Conditions for Injectivity of Generalized Polynomial Maps with Applications to Chemical Reaction Networks and Real Algebraic Geometry Anne Shiu Texas A&M University Solving Polynomial Equations workshop
More informationarxiv: v2 [math.ds] 12 Jul 2011
1 CHEMICAL REACTION SYSTEMS WITH TORIC STEADY STATES MERCEDES PÉREZ MILLÁN, ALICIA DICKENSTEIN, ANNE SHIU, AND CARSTEN CONRADI arxiv:11021590v2 [mathds] 12 Jul 2011 Abstract Mass-action chemical reaction
More informationComputational Methods in Systems and Synthetic Biology
Computational Methods in Systems and Synthetic Biology François Fages EPI Contraintes, Inria Paris-Rocquencourt, France 1 / 62 Need for Abstractions in Systems Biology Models are built in Systems Biology
More information/ / MET Day 000 NC1^ INRTL MNVR I E E PRE SLEEP K PRE SLEEP R E
05//0 5:26:04 09/6/0 (259) 6 7 8 9 20 2 22 2 09/7 0 02 0 000/00 0 02 0 04 05 06 07 08 09 0 2 ay 000 ^ 0 X Y / / / / ( %/ ) 2 /0 2 ( ) ^ 4 / Y/ 2 4 5 6 7 8 9 2 X ^ X % 2 // 09/7/0 (260) ay 000 02 05//0
More information2. At quasi-steady state or equilibrium, the net in-flux of the carrier-substrate complex CS is balanced by the net out-flux of the free carrier C.
Facilitated Transport Instructor: Nam un Wang facilitmcd Process escription In facilitated transport, a carrier molecule C binds to the substrate to form a carrier-substrate complex C at the outer side
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 informationk -2 k 2 K F I-N k 3 k -1 K F U-I k 1 NSP
Symposium on Protein Misfolding Diseases Biochemistry & Molecular Biology website May 1-2, 2007 University of Massachusetts Amherst Tuesday, May 1 1:00-1:30pm Registration 1:30-3:30pm Panel discussion:
More informationDynamics and Mechanism of A Quorum Sensing Network Regulated by Small RNAs in Vibrio Harveyi
Commun. Theor. Phys. 55 (2011) 465 472 Vol. 55, No. 3, March 15, 2011 Dynamics and Mechanism of A Quorum Sensing Network Regulated by Small RNAs in Vibrio Harveyi SHEN Jian-Wei ( å) Institute of Applied
More informationSystems Biology in Photosynthesis INTRODUCTION
1 / 26 Systems Biology in Photosynthesis INTRODUCTION Rainer Machné Institute for Theoretical Chemistry, University of Vienna, Austria PSI - Photon System Instruments, Czech Republic Brno, April, 2011
More information{ } is an asymptotic sequence.
AMS B Perturbation Methods Lecture 3 Copyright by Hongyun Wang, UCSC Recap Iterative method for finding asymptotic series requirement on the iteration formula to make it work Singular perturbation use
More information2. Higher-order Linear ODE s
2. Higher-order Linear ODE s 2A. Second-order Linear ODE s: General Properties 2A-1. On the right below is an abbreviated form of the ODE on the left: (*) y + p()y + q()y = r() Ly = r() ; where L is the
More informationChemistry C2407 Problem Set 2. Solutions to Oxtoby Even Numbered Problems
Chemistry C2407 Problem Set 2 Solutions to Oxtoby Even Numbered Problems Problem 13.32 Reaction Mechanism Elementary Steps: A B + C, A k 1 B + C, C + D k 1 k 3 k 2 E, C + D k 2 E E F Make steady state
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 informationEnzyme Kinetics 2014
V 41 Enzyme Kinetics 2014 Atkins Ch.23, Tinoco 4 th -Ch.8 Enzyme rxn example Catalysis/Mechanism: E + S k -1 ES k 1 ES E is at beginning and k 2 k -2 E + P at end of reaction Catalyst: No consumption of
More informationBiological Pathways Representation by Petri Nets and extension
Biological Pathways Representation by and extensions December 6, 2006 Biological Pathways Representation by and extension 1 The cell Pathways 2 Definitions 3 4 Biological Pathways Representation by and
More informationInference of genomic network dynamics with non-linear ODEs
Inference of genomic network dynamics with non-linear ODEs University of Groningen e.c.wit@rug.nl http://www.math.rug.nl/ ernst 11 June 2015 Joint work with Mahdi Mahmoudi, Itai Dattner and Ivan Vujačić
More informationSystem Biology - Deterministic & Stochastic Dynamical Systems
System Biology - Deterministic & Stochastic Dynamical Systems System Biology - Deterministic & Stochastic Dynamical Systems 1 The Cell System Biology - Deterministic & Stochastic Dynamical Systems 2 The
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 informationMichaelis Menten Kinetics- Identical Independent Binding Sites
Michaelis Menten Kinetics- Identical Independent Binding Sites Dr. M. Vijayalakshmi School of Chemical and Biotechnology SASTRA University Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 8 Table
More informationInteractions between proteins, DNA and RNA. The energy, length and time coordinate system to find your way in the cell
Interactions between proteins, DNA and RNA The energy, length and time coordinate system to find your way in the cell Scanning/atomic force microscope (SFM/AFM) mirror laser beam photo diode fluid cell
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 informationChapter 4: Partial differentiation
Chapter 4: Partial differentiation It is generally the case that derivatives are introduced in terms of functions of a single variable. For example, y = f (x), then dy dx = df dx = f. However, most of
More informationENZYME KINETICS. What happens to S, P, E, ES?
ENZYME KINETICS Go to lecture notes and/or supplementary handouts for the following: 1 Basic observations in enzyme inetics 2 Michaelis-Menten treatment of enzyme inetics 3 Briggs-Haldane treatment of
More informationGlobal Qualitative Analysis for a Ratio-Dependent Predator Prey Model with Delay 1
Journal of Mathematical Analysis and Applications 266, 401 419 (2002 doi:10.1006/jmaa.2001.7751, available online at http://www.idealibrary.com on Global Qualitative Analysis for a Ratio-Dependent Predator
More informationBiophysics 297: Modeling Complex Biological Systems Problem Set 1 Due: 5/12/2004
Biophysics 97: odeling Complex Biological Systems Problem Set 1 Due: 5/1/004 Problem 1: Analyzing a phase space dx Analyze the following equation graphically by plotting against x. ind the fixed points
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 informationSynchronization, Chaos, and the Dynamics of Coupled Oscillators. Supplemental 1. Winter Zachary Adams Undergraduate in Mathematics and Biology
Synchronization, Chaos, and the Dynamics of Coupled Oscillators Supplemental 1 Winter 2017 Zachary Adams Undergraduate in Mathematics and Biology Outline: The shift map is discussed, and a rigorous proof
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationCHEM 251 (4 credits): Description
CHEM 251 (4 credits): Intermediate Reactions of Nucleophiles and Electrophiles (Reactivity 2) Description: An understanding of chemical reactivity, initiated in Reactivity 1, is further developed based
More informationCrime Modeling with Lévy Flights
Crime Modeling with Lévy Flights Jonah Breslau, 1 Sorathan (Tum) Chaturapruek, 2 Daniel Yazdi 3 Mentors: Professor Theodore Kolokolnikov, 4 Professor Scott McCalla 3 1 Pomona College 2 Harvey Mudd College
More informationAn Introduction to Metabolism
An Introduction to Metabolism PREFACE The living cell is a chemical factory with thousands of reactions taking place, many of them simultaneously This chapter is about matter and energy flow during life
More informationA GLOBAL CONVERGENCE RESULT FOR PROCESSIVE MULTISITE PHOSPHORYLATION SYSTEMS
1 A GLOBAL CONVERGENCE RESULT FOR PROCESSIVE MULTISITE PHOSPHORYLATION SYSTEMS CARSTEN CONRADI AND ANNE SHIU Abstract Multisite phosphorylation plays an important role in intracellular signaling There
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 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 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 informationAn Introduction to Metabolism
An Introduction to Metabolism PREFACE The living cell is a chemical factory with thousands of reactions taking place, many of them simultaneously This chapter is about matter and energy flow during life
More informationMIDTERM REVIEW FOR MATH The limit
MIDTERM REVIEW FOR MATH 500 SHUANGLIN SHAO. The limit Define lim n a n = A: For any ε > 0, there exists N N such that for any n N, a n A < ε. The key in this definition is to realize that the choice of
More informationthen the substitution z = ax + by + c reduces this equation to the separable one.
7 Substitutions II Some of the topics in this lecture are optional and will not be tested at the exams. However, for a curious student it should be useful to learn a few extra things about ordinary differential
More informationOscillations in the expression of a self-repressed gene induced by a slow transcriptional dynamics
Oscillations in the expression of a self-repressed gene induced by a slow transcriptional dynamics Pierre-Emmanuel Morant, Quentin Thommen, François Lemaire, Constant Vandermoëre, Benjamin Parent, Marc
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 information