Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
|
|
- Gwenda Hutchinson
- 5 years ago
- Views:
Transcription
1 Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA
2 Molecular biology and Biochemical kinetics Molecular biology is one of the most important and rapidly developing areas in the life sciences, and now forms the basis of subjects such as physiology, immunology and genetics. Cellular biology is the study of cells, which make up all living creatures, and which occupy an intermediate level of biological complexity between molecules and multicellular organisms. Biochemical kinetics concerns the concentrations of chemical substances in biological systems as functions of time. Biochemical process are often controlled by enzyme catalysts. Enzymes are proteins that catalyze (i.e. accelerate) chemical reactions. Enzymes are biochemical catalysts. In these reactions, the molecules at the beginning of the process are called substrates, and the enzyme converts these into different molecules, the products. Almost all processes in the cell need enzymes in order to occur at significant rates.
3 (Bio)chemical reaction ma + nb k 1 pc + qd
4 (Bio)chemical reaction ma + nb k 1 pc + qd Reaction rate R = 1 m da dt = 1 db n dt = 1 dc p dt = 1 dd q dt
5 (Bio)chemical reaction ma + nb k 1 pc + qd Reaction rate R = 1 m da dt = 1 db n dt = 1 dc p dt = 1 dd q dt Law of mass action: the reaction rate is proportional to the product of the concentrations of the participating molecules so R = k 1 A m B n, (m + n reaction order)
6 (Bio)chemical reaction ma + nb k 1 pc + qd Reaction rate R = 1 m da dt = 1 db n dt = 1 dc p dt = 1 dd q dt Law of mass action: the reaction rate is proportional to the product of the concentrations of the participating molecules so R = k 1 A m B n, (m + n reaction order) A = mk 1 A m B n, B = nk 1 A m B n, C = pk 1 A m B n, D = qk 1 A m B n.
7 (Bio)chemical reaction ma + nb k 1 pc + qd Reaction rate R = 1 m da dt = 1 db n dt = 1 dc p dt = 1 dd q dt Law of mass action: the reaction rate is proportional to the product of the concentrations of the participating molecules so R = k 1 A m B n, (m + n reaction order) A = mk 1 A m B n, B = nk 1 A m B n, C = pk 1 A m B n, D = qk 1 A m B n. Example: Lotka reaction [Lotka, 1920] A + X k 1 2X, X + Y k 2 2Y, Y k 3 P where A is the grass, X is the rabbits, Y is the foxes, and P is the dead foxes. Assuming that A is a constant (infinite amount of grass). X = k 1 AX k 2 XY, Y = k 2 XY k 3 Y
8 Reversible chemical reaction ma + nb k 1 k2 pc + qd
9 Reversible chemical reaction ma + nb k 1 k2 pc + qd Reaction rate R = 1 m da dt = 1 db n dt = 1 dc p dt = 1 dd q dt
10 Reversible chemical reaction ma + nb k 1 k2 pc + qd Reaction rate R = 1 m R = k 1 A m B n k 2 C p D q, da dt = 1 db n dt = 1 dc p dt = 1 dd q dt A = m(k 1 A m B n k 2 C p D q ), B = n(k 1 A m B n k 2 C p D q ), C = p(k 1 A m B n k 2 C p D q ), D = q(k 1 A m B n k 2 C p D q ).
11 Reversible chemical reaction ma + nb k 1 k2 pc + qd Reaction rate R = 1 m R = k 1 A m B n k 2 C p D q, da dt = 1 db n dt = 1 dc p dt = 1 dd q dt A = m(k 1 A m B n k 2 C p D q ), B = n(k 1 A m B n k 2 C p D q ), C = p(k 1 A m B n k 2 C p D q ), D = q(k 1 A m B n k 2 C p D q ). Example: single molecule reaction A k 1 k2 B A = k 1 A + k 2 B, B = k 1 A k 2 B. Solution: A(t) + B(t) = A 0 + B 0, then A = k 1 A + k 2 (A 0 + B 0 A) = k 2 (A 0 + B 0 ) (k 1 + k 2 )A So lim A(t) = k 2 (A 0 + B 0 ), and lim t k 1 + k B(t) = k 1 (A 0 + B 0 ). 2 t k 1 + k 2 The reaction reaches an equilibrium asymptotically. In equilibrium, A B = k 2 k 1.
12 Autocatalytic reaction A chemical reaction is autocatalytic if the reaction product itself is the catalyst for that reaction. ma + nb k 1 (n + p)b A = mk 1 A m B n, B = pk 1 A m B n ma + nb k 1 k2 (n + p)b A = mk 1 A m B n + mk 2 B n+p, B = pk 1 A m B n pk 2 B n+p. Examples: Lotka reaction [Schnakenberg, 1979] 2X + Y 3X, A Y, X B x = x 2 y x + b, y = x 2 y + a. Brusselator [Prigogine, 1980] A X, 2X + Y 3X, B + X Y + D, X E x = a + x 2 y (b + 1)x, y = bx x 2 y. [Gray-Scott, 1983] U + 2V 3V, V k P, and feeding U (amount F ), removing U and V (amount F ) U = UV 2 + F (1 U), V = UV 2 (F + k)v.
13 Simple autocatalytic reactions A + nx k 1 (n + 1)X k2
14 Simple autocatalytic reactions A + nx k 1 (n + 1)X k2 Mass action kinetics: x(t)=amount of X molecules, a(t)=amount of A molecules dx dt = k 1x n a k 2 x n+1, da dt = k 1x n a + k 2 x n+1.
15 Simple autocatalytic reactions A + nx k 1 (n + 1)X k2 Mass action kinetics: x(t)=amount of X molecules, a(t)=amount of A molecules dx dt = k 1x n a k 2 x n+1, x(t) + a(t)=constant, assumed to be 1. Then da dt = k 1x n a + k 2 x n+1. dx dt = k 1x n (1 x) k 2 x n+1 = k 1 x n (k 1 + k 2 )x n+1.
16 Simple autocatalytic reactions A + nx k 1 (n + 1)X k2 Mass action kinetics: x(t)=amount of X molecules, a(t)=amount of A molecules dx dt = k 1x n a k 2 x n+1, x(t) + a(t)=constant, assumed to be 1. Then Assume k 1 = k 2 = 1. da dt = k 1x n a + k 2 x n+1. dx dt = k 1x n (1 x) k 2 x n+1 = k 1 x n (k 1 + k 2 )x n+1. n = 0: Unimolecular reaction A X (uncatalyzed) dx dt = 1 2x, solution x(t) = (x 0 0.5)e 2t.
17 Solution of autocatalytic equations n = 1: Bimolecular reaction A + X 2X (first order) (logistic equation!) dx dt = x x 0 2x2, solution x(t) = 2x 0 + (1 2x 0 )e t.
18 Solution of autocatalytic equations n = 1: Bimolecular reaction A + X 2X (first order) (logistic equation!) dx dt = x x 0 2x2, solution x(t) = 2x 0 + (1 2x 0 )e t. n = 2: Trimolecular reaction A + 2X dx dt = x2 2x 3, solution t = 2 ln 3X (second order) ( x(1 2x0 ) x 0 (1 2x) ) + 1 x 0 1 x.
19 Solution of autocatalytic equations n = 1: Bimolecular reaction A + X 2X (first order) (logistic equation!) dx dt = x x 0 2x2, solution x(t) = 2x 0 + (1 2x 0 )e t. n = 2: Trimolecular reaction A + 2X dx dt = x2 2x 3, solution t = 2 ln 3X (second order) ( x(1 2x0 ) x 0 (1 2x) ) + 1 x 0 1 x.
20 Solution of autocatalytic equations n = 1: Bimolecular reaction A + X 2X (first order) (logistic equation!) dx dt = x x 0 2x2, solution x(t) = 2x 0 + (1 2x 0 )e t. n = 2: Trimolecular reaction A + 2X dx dt = x2 2x 3, solution t = 2 ln 3X (second order) ( x(1 2x0 ) x 0 (1 2x) ) + 1 x 0 1 x. Red: n = 0; Blue: n = 1; Green: n = 2 k 1 = k 2 = 1, x 0 = 0.05
21 Solution of autocatalytic equations Red: n = 0; Blue: n = 1; Green: n = 2 k 2 = 0, x 0 = 0.01, k 1 are different for three curves and k 1 are chosen so x(5) = 0.5 The solutions for n 1 are all sigmoid function with limits 0 and 1 at t = and t =.
22 Sigmoid functions We define a sigmoid function to be y : (, ) (0, ) such that y (t) > 0, y (t)(t t 0 ) < 0 for all t t 0, lim y(t) = 0. In general any solution of y = yf (y) t with decreasing f (y) is a sigmoid function. In many cases, lim y(t) = N (the t carrying capacity when f (N) = 0), and it is also called a generalized logistic curve. One can choose f (y) to fit the data to an appropriate sigmoid function. Typical sigmoid functions: 1 u(t) = 1 + ae k(t t 0) (t = t 0 is the inflection point where u (t 0 ) = 0) 1 u(t) = (1 + ae k(t t0) ) q (t = t 0 is the inflection point where u (t 0 ) = 0) 0.5 u(t) = 0.5[tanh(t) + 1] = 1 + [ e 2t 2 t ] u(t) = 0.5[erf (t) + 1] = 0.5 π e s2 ds (t + 1) [ π ] u(t) = t 2 + 1, u(t) = arctan(t) + 1 In many other problems, sigmoid functions have range ( 1, 1) instead of (0, 1)
23 Sigmoid functions
24 Hill function Hill function [Hill, 1919] is a function in form f (u) = up h p, for u 0, p > 0. The + up point u = h is the half-saturation point since f (h) = 1/2, the half of the asymptotic limit. When p > 1, u = h is also the inflection point where f (u) = 0. When p =, the limit is the step function f (u) = 0 for u < h and f (u) = 1 for u > h.
25 Hill function Hill function [Hill, 1919] is a function in form f (u) = up h p, for u 0, p > 0. The + up point u = h is the half-saturation point since f (h) = 1/2, the half of the asymptotic limit. When p > 1, u = h is also the inflection point where f (u) = 0. When p =, the limit is the step function f (u) = 0 for u < h and f (u) = 1 for u > h. One can define a generalized Hill function by f : [0, ) [0, 1) such that f (0) = 0, f (u) > 0 for u > 0, and lim f (u) = 1. One can obtain a Hill function by restricting u and shifting a sigmoid function.
26 Hill function Hill function [Hill, 1919] is a function in form f (u) = up h p, for u 0, p > 0. The + up point u = h is the half-saturation point since f (h) = 1/2, the half of the asymptotic limit. When p > 1, u = h is also the inflection point where f (u) = 0. When p =, the limit is the step function f (u) = 0 for u < h and f (u) = 1 for u > h. One can define a generalized Hill function by f : [0, ) [0, 1) such that f (0) = 0, f (u) > 0 for u > 0, and lim f (u) = 1. One can obtain a Hill function by restricting u and shifting a sigmoid function. In ecology, people use the Hill function as the predator functional response. [Holling 1959] p = 1 is the type II, and p > 1 is the type III.
27 Hill function Hill function [Hill, 1919] is a function in form f (u) = up h p, for u 0, p > 0. The + up point u = h is the half-saturation point since f (h) = 1/2, the half of the asymptotic limit. When p > 1, u = h is also the inflection point where f (u) = 0. When p =, the limit is the step function f (u) = 0 for u < h and f (u) = 1 for u > h. One can define a generalized Hill function by f : [0, ) [0, 1) such that f (0) = 0, f (u) > 0 for u > 0, and lim f (u) = 1. One can obtain a Hill function by restricting u and shifting a sigmoid function. In ecology, people use the Hill function as the predator functional response. [Holling 1959] p = 1 is the type II, and p > 1 is the type III. In biochemistry, people use the Hill function as the reaction rate. When p = 1, it is called Michaelis-Menton chemical kinetics. [Menten and Michaelis, 1913]
28 Hill functions
29 Example 1. Consider the chemical reactions A + B k C, B + C k A. 1 Set up a system of 3 differential equations by the law of mass action for A(t), B(t) and C(t). 2 Assume that A(0) = 2, B(0) = 1 and C(0) = 1. Show that A(t) + C(t) = 3 for all t > 0. Reduce the system in part a to a system of two equations for A and B. 3 Show that B(t) = e 3kt so lim B(t) = 0. t 4 What are lim A(t) and lim t t 2. [Schnakenberg, 1979] 2X + Y 3X, A Y, X B x = x 2 y x + b, y = x 2 y + a. 1 Find the equilibrium point, and calculate the Jacobian. 2 Under what condition, the equilibrium point is stable or unstable. 3 Fix a = 1/2 and use b > 0 as a bifurcation parameter. Find a Hopf bifurcation point b 1 where a Hopf bifurcation occurs at the positive equilibrium.
30 Nonlinear reaction rates: Michaelis-Menton In a basic enzyme reaction, an enzyme E binds to a substrate S to form a complex ES, which in turn is converted into a product P and the enzyme E. (Enzyme is a biological catalyst.) E + S k 1 k2 ES k 3 E + P
31 Nonlinear reaction rates: Michaelis-Menton In a basic enzyme reaction, an enzyme E binds to a substrate S to form a complex ES, which in turn is converted into a product P and the enzyme E. (Enzyme is a biological catalyst.) E + S k 1 k2 ES k 3 E + P Let u = [S], v = [E], w = [ES] and z = [P]. Then u = k 1 uv + k 2 w, v = k 1 uv + k 2 w + k 3 w, w = k 1 uv k 2 w k 3 w, z = k 3 w.
32 Nonlinear reaction rates: Michaelis-Menton In a basic enzyme reaction, an enzyme E binds to a substrate S to form a complex ES, which in turn is converted into a product P and the enzyme E. (Enzyme is a biological catalyst.) E + S k 1 k2 ES k 3 E + P Let u = [S], v = [E], w = [ES] and z = [P]. Then u = k 1 uv + k 2 w, v = k 1 uv + k 2 w + k 3 w, w = k 1 uv k 2 w k 3 w, z = k 3 w. 1. v + w is a constant, and assume it is E (Quasi steady state approximation) Assume v = 0 and w = 0. Then u = z E 0 k 1 k 3 u = k 2 + k 3 + k 1 u
33 Nonlinear reaction rates: Michaelis-Menton In a basic enzyme reaction, an enzyme E binds to a substrate S to form a complex ES, which in turn is converted into a product P and the enzyme E. (Enzyme is a biological catalyst.) E + S k 1 k2 ES k 3 E + P Let u = [S], v = [E], w = [ES] and z = [P]. Then u = k 1 uv + k 2 w, v = k 1 uv + k 2 w + k 3 w, w = k 1 uv k 2 w k 3 w, z = k 3 w. 1. v + w is a constant, and assume it is E (Quasi steady state approximation) Assume v = 0 and w = 0. Then u = z E 0 k 1 k 3 u = k 2 + k 3 + k 1 u So in a large network of chemical reactions, one may use S P to replace the three equations above, but the reaction rate is the Michaelis-Menton function.
34 Nonlinear reaction rates: Michaelis-Menton In a basic enzyme reaction, an enzyme E binds to a substrate S to form a complex ES, which in turn is converted into a product P and the enzyme E. (Enzyme is a biological catalyst.) E + S k 1 k2 ES k 3 E + P Let u = [S], v = [E], w = [ES] and z = [P]. Then u = k 1 uv + k 2 w, v = k 1 uv + k 2 w + k 3 w, w = k 1 uv k 2 w k 3 w, z = k 3 w. 1. v + w is a constant, and assume it is E (Quasi steady state approximation) Assume v = 0 and w = 0. Then u = z E 0 k 1 k 3 u = k 2 + k 3 + k 1 u So in a large network of chemical reactions, one may use S P to replace the three equations above, but the reaction rate is the Michaelis-Menton function. There are many more complex enzyme reactions which make chemical reaction network (large system of ODEs satisfying mass action law and other laws).
35 Modeling for the CIMA reaction [Lengyel-Epstein, 1991] Modeling of Turing Structures in the Chlorite-Iodide-Malonic Acid-Starch Reaction System. Science. MA + I 2 IMA + I + H +, CIO 2 + I CIO 2 + (1/2)I 2, CIO 2 + 4I + 4H + CI + 2I 2 + 2H 2 O [CIO 2 ], [I 2 ] and [MA] varying slowly, assumed to be constant Let I = X, CIO 2 = Y and I 2 = A. Then the reaction becomes A k 1 X, X k 2 Y, 4X + Y k 3 P Reaction rates k 1, k 2 are constants, and k 3 is proportional to u = a u 4uv 1 + u 2, v = b(u uv 1 + u 2 ). [X ] [Y ] u + [X ] 2 Here [I ] = u(t), [CIO 2 ] = v(t), a, b > 0. Key parameter: a > 0 (the feeding rate)
36 Activator-inhibitor systems In the law of mass action, a gain is represented by a + term, and a loss is represented by a term. In more general case, if the presence of one chemical u increases the generation of another chemical v, then u is called a activator, and if u decreases the generation of another chemical v, then u is an inhibitor. If u is an activator of v, and v is an inhibitor of u, then in u = f (u, v), v = g(u, v), f v < 0 and g u > 0. [Gierer-Meinhardt, 1972] a = ρ a2 h µaa + ρa, h = ρa 2 µ h h + ρ h. Here a is the activator, and h is the inhibitor. µ a and µ h are decay rates of a and h, ρ a and ρ h are the feeding rates of a and h, and ρ is the reaction rate.
37 More models Water-plant model w = a w wn 2, n = wn 2 mn w(t) water, n(t) plant a > 0: rainfall; w: evaporation; wn 2 : water uptake by plants; wn 2 : plant growth; mn: plant loss [Klausmeier, Science, 1999]
38 More models Water-plant model w = a w wn 2, n = wn 2 mn w(t) water, n(t) plant a > 0: rainfall; w: evaporation; wn 2 : water uptake by plants; wn 2 : plant growth; mn: plant loss [Klausmeier, Science, 1999] Epidemic models S = A αsi AS, I = αsi βi AI, R = βi AR assuming that S(0) + I (0) + R(0) = 1 [Kermack and McKendrick, 1927]
Lecture 27. Transition States and Enzyme Catalysis
Lecture 27 Transition States and Enzyme Catalysis Reading for Today: Chapter 15 sections B and C Chapter 16 next two lectures 4/8/16 1 Pop Question 9 Binding data for your thesis protein (YTP), binding
More informationEnzyme Reactions. Lecture 13: Kinetics II Michaelis-Menten Kinetics. Margaret A. Daugherty Fall v = k 1 [A] E + S ES ES* EP E + P
Lecture 13: Kinetics II Michaelis-Menten Kinetics Margaret A. Daugherty Fall 2003 Enzyme Reactions E + S ES ES* EP E + P E = enzyme ES = enzyme-substrate complex ES* = enzyme/transition state complex EP
More informationMichaelis-Menten Kinetics. Lecture 13: Kinetics II. Enzyme Reactions. Margaret A. Daugherty. Fall Substrates bind to the enzyme s active site
Lecture 13: Kinetics II Michaelis-Menten Kinetics Margaret A. Daugherty Fall 2003 Enzyme Reactions E + S ES ES* EP E + P E = enzyme ES = enzyme-substrate complex ES* = enzyme/transition state complex EP
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 informationHopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd
Hopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd College of William and Mary Williamsburg, Virginia 23187 Mathematical Applications in Ecology and Evolution Workshop
More informationChemical Kinetics. Topic 7
Chemical Kinetics Topic 7 Corrosion of Titanic wrec Casón shipwrec 2Fe(s) + 3/2O 2 (g) + H 2 O --> Fe 2 O 3.H 2 O(s) 2Na(s) + 2H 2 O --> 2NaOH(aq) + H 2 (g) Two examples of the time needed for a chemical
More informationMichaelis-Menton kinetics
Michaelis-Menton kinetics The rate of an enzyme catalyzed reaction in which substrate S is converted into products P depends on the concentration of the enzyme E even though the enzyme does not undergo
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 informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
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 informationTOPIC 6: Chemical kinetics
TOPIC 6: Chemical kinetics Reaction rates Reaction rate laws Integrated reaction rate laws Reaction mechanism Kinetic theories Arrhenius law Catalysis Enzimatic catalysis Fuente: Cedre http://loincognito.-iles.wordpress.com/202/04/titanic-
More information1.2. Introduction to Modeling. P (t) = r P (t) (b) When r > 0 this is the exponential growth equation.
G. NAGY ODE January 9, 2018 1 1.2. Introduction to Modeling Section Objective(s): Review of Exponential Growth. The Logistic Population Model. Competing Species Model. Overview of Mathematical Models.
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 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 information2013 W. H. Freeman and Company. 6 Enzymes
2013 W. H. Freeman and Company 6 Enzymes CHAPTER 6 Enzymes Key topics about enzyme function: Physiological significance of enzymes Origin of catalytic power of enzymes Chemical mechanisms of catalysis
More informationTheory of Ordinary Differential Equations. Stability and Bifurcation I. John A. Burns
Theory of Ordinary Differential Equations Stability and Bifurcation I John A. Burns Center for Optimal Design And Control Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute
More informationChemical Kinetics AP Chemistry Lecture Outline
Chemical Kinetics AP Chemistry Lecture Outline Name: Factors that govern rates of reactions. Generally... (1)...as the concentration of reactants increases, rate (2)...as temperature increases, rate (3)...with
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 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 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 informationUse separation of variables to solve the following differential equations with given initial conditions. y 1 1 y ). y(y 1) = 1
Chapter 11 Differential Equations 11.1 Use separation of variables to solve the following differential equations with given initial conditions. (a) = 2ty, y(0) = 10 (b) = y(1 y), y(0) = 0.5, (Hint: 1 y(y
More informationLecture 11: Enzyme Kinetics, Part I
Biological Chemistry Laboratory Biology 3515/Chemistry 3515 Spring 2018 Lecture 11: Enzyme Kinetics, Part I 13 February 2018 c David P. Goldenberg University of Utah goldenberg@biology.utah.edu Back to
More information1. Introduction to Chemical Kinetics
1. Introduction to Chemical Kinetics objectives of chemical kinetics 1) Determine empirical rate laws H 2 + I 2 2HI How does the concentration of H 2, I 2, and HI change with time? 2) Determine the mechanism
More informationElementary reactions. stoichiometry = mechanism (Cl. + H 2 HCl + H. ) 2 NO 2 ; radioactive decay;
Elementary reactions 1/21 stoichiometry = mechanism (Cl. + H 2 HCl + H. ) monomolecular reactions (decay: N 2 O 4 some isomerisations) 2 NO 2 ; radioactive decay; bimolecular reactions (collision; most
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 informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
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 informationChemical Kinetics. Chapter 13. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
Chemical Kinetics Chapter 13 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Chemical Kinetics Thermodynamics does a reaction take place? Kinetics how fast does
More informationModelling chemical kinetics
Modelling chemical kinetics Nicolas Le Novère, Institute, EMBL-EBI n.lenovere@gmail.com Systems Biology models ODE models Reconstruction of state variable evolution from process descriptions: Processes
More informationRate laws, Reaction Orders. Reaction Order Molecularity. Determining Reaction Order
Rate laws, Reaction Orders The rate or velocity of a chemical reaction is loss of reactant or appearance of product in concentration units, per unit time d[p] = d[s] The rate law for a reaction is of the
More informationLecture 15 (10/20/17) Lecture 15 (10/20/17)
Reading: Ch6; 98-203 Ch6; Box 6- Lecture 5 (0/20/7) Problems: Ch6 (text); 8, 9, 0,, 2, 3, 4, 5, 6 Ch6 (study guide-facts); 6, 7, 8, 9, 20, 2 8, 0, 2 Ch6 (study guide-applying); NEXT Reading: Ch6; 207-20
More informationTheoretical Models for Chemical Kinetics
Theoretical Models for Chemical Kinetics Thus far we have calculated rate laws, rate constants, reaction orders, etc. based on observations of macroscopic properties, but what is happening at the molecular
More informationBiochemistry 3100 Sample Problems Binding proteins, Kinetics & Catalysis
(1) Draw an approximate denaturation curve for a typical blood protein (eg myoglobin) as a function of ph. (2) Myoglobin is a simple, single subunit binding protein that has an oxygen storage function
More informationLokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt
Periodic solutions A periodic solution is a solution (x(t), y(t)) of dx = f(x, y) dt dy = g(x, y) dt such that x(t + T ) = x(t) and y(t + T ) = y(t) for any t, where T is a fixed number which is a period
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationCh 13 Rates of Reaction (Chemical Kinetics)
Ch 13 Rates of Reaction (Chemical Kinetics) Reaction Rates and Kinetics - The reaction rate is how fast reactants are converted to products. - Chemical kinetics is the study of reaction rates. Kinetics
More informationModule 6 : Reaction Kinetics and Dynamics Lecture 28 : Elementary Reactions and Reaction Mechanisms
Module 6 : Reaction Kinetics and Dynamics Lecture 28 : Elementary Reactions and Reaction Mechanisms Objectives In this Lecture you will learn to do the following Define what is an elementary reaction.
More information1.2. Introduction to Modeling
G. NAGY ODE August 30, 2018 1 Section Objective(s): Population Models Unlimited Resources Limited Resources Interacting Species 1.2. Introduction to Modeling 1.2.1. Population Model with Unlimited Resources.
More informationAnalyze the roles of enzymes in biochemical reactions
ENZYMES and METABOLISM Elements: Cell Biology (Enzymes) Estimated Time: 6 7 hours By the end of this course, students will have an understanding of the role of enzymes in biochemical reactions. Vocabulary
More informationCatalysis. v 0 no catalyst v c -- catalyst present. v c. dt with no catalyst) (v c = -d[a]/dt dt with a catalyst)
Catalysis Catalysis provides an additional mechanism by which reactants can be converted to products. The alternative mechanism has a lower activation energy than the reaction in the absence of a catalyst.
More informationWe have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v),
Let us apply the approach presented above to the analysis of population dynamics models. 9. Lotka-Volterra predator-prey model: phase plane analysis. Earlier we introduced the system of equations for prey
More information= dc A dt. The above is a bimolecular elementary reaction. A unimolecular elementary reaction might be HO 2 H + O 2
The above is a bimolecular elementary reaction. A unimolecular elementary reaction might be HO 2 H + O 2 HO 2 just dissociates without any other influence. Rate Laws for Elementary Reactions: 1) A Fragments,
More informationName Student number. UNIVERSITY OF GUELPH CHEM 4540 ENZYMOLOGY Winter 2002 Quiz #1: February 14, 2002, 11:30 13:00 Instructor: Prof R.
UNIVERSITY OF GUELPH CHEM 4540 ENZYMOLOGY Winter 2002 Quiz #1: February 14, 2002, 11:30 13:00 Instructor: Prof R. Merrill Instructions: Time allowed = 90 minutes. Total marks = 30. This quiz represents
More informationC2 Differential Equations : Computational Modeling and Simulation Instructor: Linwei Wang
C2 Differential Equations 4040-849-03: Computational Modeling and Simulation Instructor: Linwei Wang Part IV Dynamic Systems Equilibrium: Stable or Unstable? Equilibrium is a state of a system which does
More informationModeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool
Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationCHAPTER 1: ENZYME KINETICS AND APPLICATIONS
CHAPTER 1: ENZYME KINETICS AND APPLICATIONS EM 1 2012/13 ERT 317 BIOCHEMICAL ENGINEERING Course details Credit hours/units : 4 Contact hours : 3 hr (L), 3 hr (P) and 1 hr (T) per week Evaluations Final
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 information6. Well-Stirred Reactors III
6. Well-Stirred Reactors III Reactors reaction rate or reaction velocity defined for a closed system of uniform pressure, temperature, and composition situation in a real reactor is usually quite different
More informationLecture 4 STEADY STATE KINETICS
Lecture 4 STEADY STATE KINETICS The equations of enzyme kinetics are the conceptual tools that allow us to interpret quantitative measures of enzyme activity. The object of this lecture is to thoroughly
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 informationThe functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z).
Review of Multi-variable calculus: The functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z). The spatial variable represents the environment where the species
More information1.Introduction: 2. The Model. Key words: Prey, Predator, Seasonality, Stability, Bifurcations, Chaos.
Dynamical behavior of a prey predator model with seasonally varying parameters Sunita Gakkhar, BrhamPal Singh, R K Naji Department of Mathematics I I T Roorkee,47667 INDIA Abstract : A dynamic model based
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 informationFrom Friday s material
5.111 Lecture 35 35.1 Kinetics Topic: Catalysis Chapter 13 (Section 13.14-13.15) From Friday s material Le Chatelier's Principle - when a stress is applied to a system in equilibrium, the equilibrium tends
More informationCoupled differential equations
Coupled differential equations Example: dy 1 dy 11 1 1 1 1 1 a y a y b x a y a y b x Consider the case with b1 b d y1 a11 a1 y1 y a1 ay dy y y e y 4/1/18 One way to address this sort of problem, is to
More informationODE Math 3331 (Summer 2014) June 16, 2014
Page 1 of 12 Please go to the next page... Sample Midterm 1 ODE Math 3331 (Summer 2014) June 16, 2014 50 points 1. Find the solution of the following initial-value problem 1. Solution (S.O.V) dt = ty2,
More informationA. One-Substrate Reactions (1) Kinetic concepts
A. One-Substrate Reactions (1) Kinetic concepts (2) Kinetic analysis (a) Briggs-Haldane steady-state treatment (b) Michaelis constant (K m ) (c) Specificity constant (3) Graphical analysis (4) Practical
More informationChem Lecture 4 Enzymes Part 1
Chem 452 - Lecture 4 Enzymes Part 1 Question of the Day: Enzymes are biological catalysts. Based on your general understanding of catalysts, what does this statement imply about enzymes? Introduction Enzymes
More informationAppendix 4. Some Equations for Curve Fitting
Excep for Scientists and Engineers: Numerical Methods by E. Joseph Billo Copyright 0 2007 John Wiley & Sons, Inc. Appendix 4 Some Equations for Curve Fitting This appendix describes a number of equation
More informationChemical Kinetics. Kinetics is the study of how fast chemical reactions occur. There are 4 important factors which affect rates of reactions:
Chemical Kinetics Kinetics is the study of how fast chemical reactions occur. There are 4 important factors which affect rates of reactions: reactant concentration temperature action of catalysts surface
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 informationAfter lectures by. disappearance of reactants or appearance of. measure a reaction rate we monitor the. Reaction Rates (reaction velocities): To
Revised 3/21/2017 After lectures by Dr. Loren Williams (GeorgiaTech) Protein Folding: 1 st order reaction DNA annealing: 2 nd order reaction Reaction Rates (reaction velocities): To measure a reaction
More informationLecture 11: Enzymes: Kinetics [PDF] Reading: Berg, Tymoczko & Stryer, Chapter 8, pp
Lecture 11: Enzymes: Kinetics [PDF] Reading: Berg, Tymoczko & Stryer, Chapter 8, pp. 216-225 Updated on: 2/4/07 at 9:00 pm Key Concepts Kinetics is the study of reaction rates. Study of enzyme kinetics
More informationComputational models: Reaction Diffusion. Matthias Vigelius Week 6
Computational models: Reaction Diffusion Matthias Vigelius Week 6 Reactions (Chemical) reactions Reaction is a transformation from reactants into products: A + 2B D + 3E + k k [ ] Reactions occur with
More informationMath 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations
Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Junping Shi College of William and Mary, USA Equilibrium Model: x n+1 = f (x n ), here f is a nonlinear function
More informationComputer Simulation and Data Analysis in Molecular Biology and Biophysics
Victor Bloomfield Computer Simulation and Data Analysis in Molecular Biology and Biophysics An Introduction Using R Springer Contents Part I The Basics of R 1 Calculating with R 3 1.1 Installing R 3 1.1.1
More informationChemical reaction networks and diffusion
FYTN05 Computer Assignment 2 Chemical reaction networks and diffusion Supervisor: Adriaan Merlevede Office: K336-337, E-mail: adriaan@thep.lu.se 1 Introduction This exercise focuses on understanding and
More informationA population is modeled by the differential equation
Math 2, Winter 2016 Weekly Homework #8 Solutions 9.1.9. A population is modeled by the differential equation dt = 1.2 P 1 P ). 4200 a) For what values of P is the population increasing? P is increasing
More informationFrom cell biology to Petri nets. Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE
From cell biology to Petri nets Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE Biology = Concentrations Breitling / 2 The simplest chemical reaction A B irreversible,
More informationResearch Article Studies on the Existence of Unstable Oscillatory Patterns Bifurcating from Hopf Bifurcations in a Turing Model
Applied Mathematics Volume 214, Article ID 574921, 5 pages http://dx.doi.org/1.1155/214/574921 Research Article Studies on the Existence of Unstable Oscillatory Patterns Bifurcating from Hopf Bifurcations
More informationProblem Set Number 1, j/2.036j MIT (Fall 2014)
Problem Set Number 1, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept., Cambridge, MA 02139) September 9, 2014 Due Fri., September 19, 2014. 1 Get equation from phase line portrait problem
More informationSOME IMPORTANT GRAPH REPRESENTATION
SOME IMPORTNT GRPH REPRESENTTION Use in chemical inetics : () y mx + C () y mx + C (3) y mx UGC POINT: Institute For CSIR, NET, GTE, JM Ph: 96546855, -685555, 6568 (4) y mx Slope Intercept m tan zero (5)
More informationStability and bifurcation in a two species predator-prey model with quintic interactions
Chaotic Modeling and Simulation (CMSIM) 4: 631 635, 2013 Stability and bifurcation in a two species predator-prey model with quintic interactions I. Kusbeyzi Aybar 1 and I. acinliyan 2 1 Department of
More informationLecture 13: Data Analysis for the V versus [S] Experiment and Interpretation of the Michaelis-Menten Parameters
Biological Chemistry Laboratory Biology 3515/Chemistry 3515 Spring 2018 Lecture 13: Data Analysis for the V versus [S] Experiment and Interpretation of the Michaelis-Menten Parameters 20 February 2018
More informationChapter 8. Enzymes: basic concept and kinetics
Chapter 8 Enzymes: basic concept and kinetics Learning objectives: mechanism of enzymatic catalysis Michaelis -Menton Model Inhibition Single Molecule of Enzymatic Reaction Enzymes: catalysis chemical
More informationEnzyme Kinetics. Michaelis-Menten Theory Dehaloperoxidase: Multi-functional Enzyme. NC State University
Enzyme Kinetics Michaelis-Menten Theory Dehaloperoxidase: Multi-functional Enzyme NC State University Michaelis-Menton kinetics The rate of an enzyme catalyzed reaction in which substrate S is converted
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 informationBreakdown of Pattern Formation in Activator-Inhibitor Systems and Unfolding of a Singular Equilibrium
Breakdown of Pattern Formation in Activator-Inhibitor Systems and Unfolding of a Singular Equilibrium Izumi Takagi (Mathematical Institute, Tohoku University) joint work with Kanako Suzuki (Institute for
More information6.1 Enzyme Kinetics: Basic Enzyme Reaction
6. Reaction Kinetics 6.1 Enzyme Kinetics: Basic Enzyme Reaction Biochemical reactions are continually taking place in all living organisms and most of them involve proteins called enzymes, which act as
More informationGlobal Stability Analysis on a Predator-Prey Model with Omnivores
Applied Mathematical Sciences, Vol. 9, 215, no. 36, 1771-1782 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.512 Global Stability Analysis on a Predator-Prey Model with Omnivores Puji Andayani
More informationChemical Kinetics. What quantities do we study regarding chemical reactions? 15 Chemical Kinetics
Chemical Kinetics Chemical kinetics: the study of reaction rate, a quantity conditions affecting it, the molecular events during a chemical reaction (mechanism), and presence of other components (catalysis).
More informationCh 13 Chemical Kinetics. Modified by Dr. Cheng-Yu Lai
Ch 13 Chemical Kinetics Modified by Dr. Cheng-Yu Lai Outline 1. Meaning of reaction rate 2. Reaction rate and concentration 3. Writing a Rate Law 4. Reactant concentration and time 5. Reaction rate and
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 informationChapter 14. Chemical Kinetics
Chapter 14. Chemical Kinetics 14.1 Factors that Affect Reaction Rates The speed at which a chemical reaction occurs is the reaction rate. Chemical kinetics is the study of how fast chemical reactions occur.
More informationSystems of Ordinary Differential Equations
Systems of Ordinary Differential Equations Scott A. McKinley October 22, 2013 In these notes, which replace the material in your textbook, we will learn a modern view of analyzing systems of differential
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 informationDerivation of Reaction-Diffusion Equations
Chapter 1 Derivation of Reaction-Diffusion Equations 1.1 Fick s Law Diffusion mechanism models the movement of many individuals in an environment or media. The individuals can be very small such as basic
More informationChapter 3. Chemistry of Life
Chapter 3 Chemistry of Life Content Objectives Write these down! I will be able to identify: Where living things get energy. How chemical reactions occur. The functions of lipids. The importance of enzymes
More informationENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 7
ENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 7 KINETICS OF ENZYME CATALYSED REACTIONS (CONTD.) So in the last lecture we
More informationChapter 6- An Introduction to Metabolism*
Chapter 6- An Introduction to Metabolism* *Lecture notes are to be used as a study guide only and do not represent the comprehensive information you will need to know for the exams. The Energy of Life
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 informationChapter 14. Enzyme Kinetics
Chapter 4. Enzyme inetics Chemical kinetics Elementary reactions A P (Oerall stoichiometry) (ntermediates) Rate equations aa bb zz P Rate k[a] a [B] b [Z] z k: rate constant The order of the reaction (ab
More informationOn predator-prey models
Predator-prey On models ddd Department of Mathematics College of William and Mary Math 41/CSUMS Talk February 3, 1 Collaborators Sze-Bi Hsu (Tsinghua University, Hsinchu, Taiwan) Junjie Wei (Harbin Institute
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationCHEM-E3205 BIOPROCESS OPTIMIZATION AND SIMULATION
CHEM-E3205 BIOPROCESS OPTIMIZATION AND SIMULATION TERO EERIKÄINEN ROOM D416d tero.eerikainen@aalto.fi COURSE LECTURES AND EXERCISES Week Day Date Time Place Lectures/Execises 37 Mo 12.9.2016 10:15-11:45
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
More informationChapters 8.1 & 8.2 Practice Problems
EXPECTED SKILLS: Chapters 8.1 & 8. Practice Problems Be able to verify that a given function is a solution to a differential equation. Given a description in words of how some quantity changes in time
More informationStability Analysis and Numerical Solution for. the Fractional Order Biochemical Reaction Model
Nonlinear Analysis and Differential Equations, Vol. 4, 16, no. 11, 51-53 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/nade.16.6531 Stability Analysis and Numerical Solution for the Fractional
More informationWhat is an enzyme? Lecture 12: Enzymes & Kinetics I Introduction to Enzymes and Kinetics. Margaret A. Daugherty Fall 2004 KEY FEATURES OF ENZYMES
Lecture 12: Enzymes & Kinetics I Introduction to Enzymes and Kinetics Margaret A. Daugherty Fall 2004 What is an enzyme? General Properties Mostly proteins, but some are actually RNAs Biological catalysts
More information