Download
1 / 85
880 likes | 926 Views
Modeling and Analysis Techniques in Systems Biology. CS 6221 Lecture 2 P.S. Thiagarajan. Acknowledgment. Many of the PDF images that appear in the slides to follow are taken from the text book “Systems Biology in Practice” by E. Klipp et.al. The Role of Chemical Reactions.
E N D
Modeling and Analysis Techniques in Systems Biology. CS 6221 Lecture 2 P.S. Thiagarajan
Acknowledgment • Many of the PDF images that appear in the slides to follow are taken from the text book “Systems Biology in Practice” by E. Klipp et.al.
The Role of Chemical Reactions Bio-Chemical reactions Metabolic pathways Signaling pathways Gene regulatory networks A network of Bio-Chemical reactions Interacting Bio-Chemical networks Cell functions
The Role of Chemical Reactions Reaction kinetics Bio-Chemical reactions Metabolic pathways Signaling pathways Gene regulatory networks A network of Bio-Chemical reactions Interacting Bio-Chemical networks Cell functions
Rate Laws • Rate law: • An equation that relates the concentrations of the reactants to the rate. • Differential equations are often used to describe these laws. • Assumption: The reactants participating in the reactions are abundant.
Reaction Kinetics • Kinetics: • Determine reaction rates • Fix reaction law and • determine reaction rate constant • Solve the equation capturing the dynamics. • The reaction rate for a product or reactant in a particular reaction: • the amount (in moles or mass units) per unit time per unit volume that is formed or removed.
Rate Laws Mass action law: The reaction rate is proportional to the probability of collision of the reactants Proportional to the concentration of the reactants to the power of their molecularities.
Mass action law V S1 + S2 P V = k. [S1] [S2] [S1] is the concentration (Moles/ litre) of S1 [S2] is the concentration (Moles/ litre) of S k is the rate constant V, the rate of the reaction
Mass-action Kinetics k1 k3 E + S ES E + P k 2
Assuming mass law kinetics we can write down a system of ordinary differential equations for the 6 species. • But we don’t know how to solve systems of ordinary (non-linear) differential equations even for dimension 4! • We must resort to numerical integration.
Michaelis-Menton Kinetics • Describes the rate of enzyme-mediated reactions in an amalgamated fashion: • Based on mass action law. • Subject to some assumptions • Enzymes • Protein (bio-)catalysts • Catalyst: • A substance that accelerates the rate of a reaction without being used up. • The speed-up can be enormous!
Enzymes • Substrate binds temporarily to the enzyme. • Lowers the activation energy needed for the reaction. • The rate at which an enzyme works is influenced by: • concentration of the substrate • Temperature • beyond a certain point, the protein can get denatured • Its 3 dimensional structure gets disrupted
Enzymes • The rate at which an enzyme works is influenced by: • The presence of inhibitors • molecules that bind to the same site as the substrate (competitive) • prevents the substrate from binding • molecules that bind to some other site of the enzyme but reduces its catalytic power (non-competitive) • pH (the concentration of hydrogen ions in a solution) • affects the 3 dimensional shape
Michaelis-Menton Kinetics k1 k3 E + S ES E + P k 2 • A reversible formation of the Enzyme-Substrate complex ES • Irreversible release of the product P from the enzyme. This is for a single substrate; no backward reaction or negligible if we focus on the initial phase of the reaction.
Michaelis-Menton Kinetics k1 k3 E + S ES E + P k 2 Use mass action law to model each reaction.
This is the rate at which P is being produced. (1) Assumption1: [ES] concentration changes much more slowly than those of [S] and [P](quasi-steady-state) We can then write:
Michaelis-Menton Kinetics (1) (2) (Michaelis constant) Define (3)
Assumption1: [ES] concentration changes much more slowly than those of [S] and [P](quasi-steady-state) Assumption2: The total enzyme concentration does not change with time. [E0] = [E] + [ES] [E0] - initial concentration
Michaelis-Menton Kinetics Vmax is achieved when all of the enzyme (E0) is substrate-bound. (assumption: [S] >> [E0]) at maximum rate, Thus,
Michaelis-Menton Kinetics This is the Michaelis-Menten equation!
Michaelis-Menton Kinetics This is the Michaelis-Menten equation! So what?
Michaelis-Menton Kinetics Consider the case: The KM of an enzyme is therefore the substrate concentration at which the reaction occurs at half of the maximum rate.
Michaelis-Menton Kinetics • KM is an indicator of the affinity that an enzyme has for a given substrate, and hence the stability of the enzyme-substrate complex. • At low [S], it is the availability of substrate that is the limiting factor. • As more substrate is added there is a rapid increase in the initial rate of the reaction.
Curve Plotting • This is not relevant anymore • Good non-linear regression techniques and LARGE amounts of computing power are available.
Reversible form of Michaelis-Menten. Variations E + S ES E + P More complicated equation but similar form.
Variations • Enzymes don’t merely accelerate reactions. • They regulate metabolism: • Their production and degradation adapted to current requirements of the cell. • Enzyme’s effectiveness targeted by inhibitors and activators (effectors).
Variations • Regulatory interactions between an enzyme and an inhibitor are characterized by: • How the enzyme binds the inhibitor I • EI, ESI or both • Which complexes can release the product • ES alone or ESI or both ES and ESI
Competitive Inhibition S and I compete for the binding place High S may out-compete I
Uncompetitive Inhibition Inhibitor binds only to the ES complex. Does not compete but inhibits by binding elsewhere and inhibiting. S can’t out-compete I.
Other forms Inhibitions • Non-competitive inhibition • Mixed inhibition • Partial inhibition
Hill Coefficients • Suppose a dimeric (two identical sub-units linked together) protein has two identical binding sites. • The binding to the first ligand (at the first site) can facilitate binding to the second ligand. • Cooperative binding. • In general, the binding of a ligand to a macromolecule is often enhanced if there are already other ligands present on the same macromolecule • The degree of cooperation is indicated by the Hill coefficient.
Hill Coefficients • A Hill coefficient of 1 indicates completely independent binding. • Independent of whether or not additional ligands are already bound. • A coefficient > 1 indicates cooperative binding. • Oxygen binding to hemoglobin: • Hill coefficient of 2.8 – 3.0
Hill’s equation Hill equation θ - fraction of ligand binding sites filled [L] - ligand concentration KM- ligand concentration producing half occupation (ligand concentration occupying half of the binding sites) n - Hill coefficient, describing cooperativity
Summary • A bio-chemical reaction is governed by a kinetic law. • Mass law, Michalis-Menten, Hill equation,… • Different laws apply under different regimes. • Each law leads to an ODE model of the reaction kinetics. • Often, with an unknown constant of proportionality. (rate constant)
