Bioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014"

Similar documents
Bioinformatics 3. V18 Kinetic Motifs. Fri, Jan 8, 2016

nutrients growth & division repellants movement

A simple model for the eukaryotic cell cycle. Andrea Ciliberto

Analysis and Simulation of Biological Systems

Network Dynamics and Cell Physiology. John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute

Life Sciences 1a: Section 3B. The cell division cycle Objectives Understand the challenges to producing genetically identical daughter cells

Plant Molecular and Cellular Biology Lecture 8: Mechanisms of Cell Cycle Control and DNA Synthesis Gary Peter

Introduction to Biological Modeling

Cell cycle regulation in the budding yeast

Network Dynamics and Cell Physiology

Lecture 10: Cyclins, cyclin kinases and cell division

Dr. Fred Cross, Rockefeller (KITP Bio Networks 3/26/2003) 1

Network Dynamics and Cell Physiology. John J. Tyson Dept. Biological Sciences Virginia Tech

Modelling Biochemical Reaction Networks. Lecture 17: Modeling the cell cycle, Part I

Basic Synthetic Biology circuits

16 The Cell Cycle. Chapter Outline The Eukaryotic Cell Cycle Regulators of Cell Cycle Progression The Events of M Phase Meiosis and Fertilization

Networks in systems biology

ANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY. Nael H. El-Farra, Adiwinata Gani & Panagiotis D. Christofides

CS-E5880 Modeling biological networks Gene regulatory networks

DNA replication. M (mitosis)

Three different fusions led to three basic ideas: 1) If one fuses a cell in mitosis with a cell in any other stage of the cell cycle, the chromosomes

Enzyme Reactions. Lecture 13: Kinetics II Michaelis-Menten Kinetics. Margaret A. Daugherty Fall v = k 1 [A] E + S ES ES* EP E + P

Michaelis-Menten Kinetics. Lecture 13: Kinetics II. Enzyme Reactions. Margaret A. Daugherty. Fall Substrates bind to the enzyme s active site

Alternating Oscillations and Chaos in a Model of Two Coupled Biochemical Oscillators Driving Successive Phases of the Cell Cycle

Getting In and Out of Mitosis*

Problem Set 2. 1 Competitive and uncompetitive inhibition (12 points) Systems Biology (7.32/7.81J/8.591J)

Introduction to Mathematical Physiology I - Biochemical Reactions

Plant Molecular and Cellular Biology Lecture 10: Plant Cell Cycle Gary Peter

56:198:582 Biological Networks Lecture 11

Biophysical Chemistry 72 (1998)

Study of Tricyclic Cascade Networks using Dynamic Optimization

Overview of the cell cycle

Dynamics of the Cell Cycle: Checkpoints, Sizers, and Timers

MATHEMATICAL MODELING OF MITOSIS

COURSE NUMBER: EH 590R SECTION: 1 SEMESTER: Fall COURSE TITLE: Computational Systems Biology: Modeling Biological Responses

Systems Biology Across Scales: A Personal View XIV. Intra-cellular systems IV: Signal-transduction and networks. Sitabhra Sinha IMSc Chennai

CELL CYCLE AND DIFFERENTIATION

Genetic transcription and regulation

Genetic transcription and regulation

Sic1 as a timer of Clb cyclin waves in the yeast cell cycle design principle of not just an inhibitor

Principles of Synthetic Biology: Midterm Exam

Dynamic Analysis of Interlocked Positive Feedback Loops

12/5/2014. The cell cycle and cell death. The cell cycle: cells duplicate their contents and divide

Regulation of the mammalian cell cycle: a model of the G 1 -to-s transition

Review (V1): Phases of Cell Cycle

Bi 8 Lecture 11. Quantitative aspects of transcription factor binding and gene regulatory circuit design. Ellen Rothenberg 9 February 2016

Multistability in the lactose utilization network of Escherichia coli

Research Article Robust Cell Size Checkpoint from Spatiotemporal Positive Feedback Loop in Fission Yeast

The Cell Cycle/Le Cycle cellulaire SMC6052/BIM6028 IRCM

Numerical analysis of a comprehensive model of M-phase control in Xenopus oocyte extracts and intact embryos

return in class, or Rm B

Analysis of a generic model of eukaryotic cell cycle regulation

Supporting Information. Methods. Equations for four regimes

56:198:582 Biological Networks Lecture 9

Chemical kinetics and catalysis

From quiescence to proliferation: Cdk oscillations drive the mammalian cell cycle

A minimal mathematical model combining several regulatory cycles from the budding yeast cell cycle

Analysis of Mode Transitions in Biological Networks

Introduction. Dagmar Iber Jörg Stelling. CSB Deterministic, SS 2015, 1.

COMPUTER SIMULATION OF DIFFERENTIAL KINETICS OF MAPK ACTIVATION UPON EGF RECEPTOR OVEREXPRESSION

Biology: Life on Earth

Zool 3200: Cell Biology Exam 5 4/27/15

Lecture 8: Temporal programs and the global structure of transcription networks. Chap 5 of Alon. 5.1 Introduction

Robust Network Topologies for Generating Switch-Like Cellular Responses

A synthetic oscillatory network of transcriptional regulators

Notes, Mathematical Cell Biology Course. Leah Edelstein-Keshet

BioControl - Week 2, Lecture 1

56:198:582 Biological Networks Lecture 10

Biological Pathways Representation by Petri Nets and extension

Src Family Kinases and Receptors: Analysis of Three Activation Mechanisms by Dynamic Systems Modeling

Activation of a receptor. Assembly of the complex

Simple kinetics of enzyme action

L'horloge circadienne, le cycle cellulaire et leurs interactions. Madalena CHAVES

Analog Electronics Mimic Genetic Biochemical Reactions in Living Cells

Enzyme reaction example of Catalysis, simplest form: E + P at end of reaction No consumption of E (ES): enzyme-substrate complex Intermediate

Core genetic module: The mixed feedback loop

CELL CYCLE AND GROWTH REGULATION

Lecture 15 (10/20/17) Lecture 15 (10/20/17)

System Modelling of Mammalian Cell Cycle Regulation Using Multi-Level Hybrid Petri Nets

Key ideas about dynamic models

Proceedings of the FOSBE 2007 Stuttgart, Germany, September 9-12, 2007

Computational models: Reaction Diffusion. Matthias Vigelius Week 6

Cell Cycle Regulation by Chlamydomonas Cyclin-Dependent Protein Kinases

Enzyme Kinetics: The study of reaction rates. For each very short segment dt of the reaction: V k 1 [S]

Building a cellular switch: more lessons from a good egg

1. Introduction to Chemical Kinetics

ENZYME KINETICS. Medical Biochemistry, Lecture 24

Name Student number. UNIVERSITY OF GUELPH CHEM 4540 ENZYMOLOGY Winter 2002 Quiz #1: February 14, 2002, 11:30 13:00 Instructor: Prof R.

NIH Public Access Author Manuscript Nat Rev Mol Cell Biol. Author manuscript; available in PMC 2009 December 21.

The cell cycle entails an ordered series of macromolecular

Computing Algebraic Functions with Biochemical Reaction Networks

Enzymes II. Dr. Mamoun Ahram Summer, 2017

V5 Cell Cycle. In cells with a nucleus (eukaryotes), the cell cycle can be divided in 2 brief periods:

Quantitative Analysis of Cellular Metabolic Dissipative, Self-Organized Structures

Multistability in the lactose utilization network of E. coli. Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J.

Modelling the fission yeast cell cycle Akos Sveiczer, John J. Tyson and Bela Novak Date received (in revised form): 20th October, 2003

Written Exam 15 December Course name: Introduction to Systems Biology Course no

After lectures by. disappearance of reactants or appearance of. measure a reaction rate we monitor the. Reaction Rates (reaction velocities): To

1 Periodic stimulations of the Incoherent Feedforward Loop network

Chapter 15 Active Reading Guide Regulation of Gene Expression

Transcription:

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 " kinetic properties of pathway components?" (2) Can high signal amplification be coupled with fast signaling?" (3) How are signaling pathways designed to ensure that they are safely off in the absence " of stimulation, yet display high signal amplification following receptor activation?" (4) How can different agonists stimulate the same pathway in distinct ways to " elicit a sustained or a transient response, which can have dramatically different consequences?" Bioinformatics 3 WS 13/14" V 20 " 2"

Linear Response" E.g., protein synthesis and degradation (see lecture V10)" S S = signal (e.g., concentration of mrna)" R = response (e.g., concentration of a protein)" R At steady state (which implies S = const):" =>" RSS linearly dependent on S" RSS! 2 1 0 0 1 2 S! k0 = 1, k1 = k2 = 2" Bioinformatics 3 WS 13/14" V 20 " 3"

R phosphorylation/dephosphorylation" S RP forward : R is converted to phosphorylated form RP" backward : RP can be dephosphorylated again to R" S + R => RP" RP => R + T" with Rtot = R + RP" phosphorylated form" T Find steady state for RP: linear until saturation" RPSS! 1 0.1 Output T proportional to RP level:" 0.01 0.01 0.1 1 10 100 S! Rtot = 1, S0 = 1" Bioinformatics 3 WS 13/14" V 20 " 4"

Enzyme: Michaelis-Menten-kinetics" S" kon" E" ES" koff" T" Reaction rate:" Steady state:" Total amount of enzyme is constant:" =>" turnover: " Bioinformatics 3 - WS 12/13" 5"

The MM-equation" Effective turnover according to MM: " Pro:" Cons:" analytical formula for turnover" curve can be easily interpreted: Vmax, KM" enzyme concentration can be ignored" less kinetic information" kon, koff, ET => Vmax, KM" Bioinformatics 3 - WS 12/13" 6"

Sigmoidal Characteristics with MM kinetics" S Same topology as before with Michaelis-Menten kinetics for phosphorylation and dephosphorylation." R RP T this means that S = Rt - RP! KM = R0! Quadratic equation for RP" => sigmoidal characteristics " (threshold behavior)" often found in signalling cascades" RPSS! 10 8 6 4 2 0 0 1 2 3 S! Rt = 10, R0 = RP0 = 1, k1 = k2 = 1" Bioinformatics 3 WS 13/14" V 20 " 7"

Graded Response" 10 RSS! 2 1 RPSS! 1 0.1 RPSS! 8 6 4 0 0 1 2 S! 0.01 0.01 0.1 1 10 100 S! 2 0 0 1 2 3 S Linear, hyperbolic, and sigmoidal characteristic give the same steady state response independent of the previous history" => no hysteresis" BUT: In fast time-dependent scenarios," delay may lead to a modified response " Bioinformatics 3 WS 13/14" V 20 " 8"

Time-dependent Sigmoidal Response" S Direct implementation:" R RP T Parameters: k1 = 1 (mol s) 1, k2 = 1 s 1, R0 = RP0 = 1 mol" Initial conditions: R = 10 mol, RP = 0" Time courses for " S = 1, 1.5, and 2," RP(0) = 0:" RP(t)" equilibrium is reached" faster for " stronger signal" t" Bioinformatics 3 WS 13/14" V 20 " 9"

Adaption - sniffer " Linear response modulated by a second species X" S X R Steady state: Rss independent of S" R changes transiently when S changes, " then goes back to its basal level." found in smell, vision, chemotaxis, " 2 1 S! X! Note: response strength!r depends on rate of change of S." => non-monotonous relation for R(S)" 0 R! 0 1 2 3 4 5 S! k1 = 30, k2 = 40, k3 = k4 = 5" Bioinformatics 3 WS 13/14" V 20 " 10"

Positive Feedback" Feedback via R and EP" => high levels of R will stay" "one-way switch" via bifurcation" Found in processes that are "final":" frog oocyte maturation, apoptosis, " Bioinformatics 3 WS 13/14" V 20 " 11"

Mutual Inhibition - Toggle Switch" Sigmoidal "threshold" in E <=> EP leads to bistable response (hysteresis): toggle switch! Converts continuous external stimulus into " two well defined stable states:" lac operon in bacteria" activation of M-phase promoting factor in frog eggs" Bioinformatics 3 WS 13/14" V 20 " 12"

Negative Feedback" S controls the "demand" for R" => homeostasis" found in biochemical pathways," no transient changes in R for steps in S" (cf. "sniffer")" Bioinformatics 3 WS 13/14" V 20 " 13"

Negative Feedback with Delay" Cyclic activation X => YP => RP => X" => Oscillations (in a range of S)" Proposed mechanism for circadian clocks" Bioinformatics 3 WS 13/14" V 20 " 14"

Circadian Clocks" CK1: casein kinase Rev-erb, ROR: retinoic acidrelated orphan nuclear receptors Cdg: clock-controlled gene(s) Ko & Takahashi Hum Mol Genet 15, R271 (2006) Bioinformatics 3 WS 13/14" V 20 " 15"

Substrate-Depletion Oscillations" R is produced in an autocatalytic reaction from X, finally depleting X " Similar to Lotka-Volterra system (autocatalysis for X, too):" Bioinformatics 3 WS 13/14" V 20 " 16"

The Cell Cycle" DNA separation" (mitosis)" Cell division" (cytokinesis)" DNA replication" cell growth" G1-phase" S-phase" When to take the next step???" Bioinformatics 3 WS 13/14" V 20 " 17"

Cell Cycle Control" Bioinformatics 3 WS 13/14" V 20 " 18"

Cell Cycle Control System" cdc = " "cell division cycle"" Tyson et al, Curr. Op. Cell Biol. 15 (2003) 221" Bioinformatics 3 WS 13/14" V 20 " 19"

Feedback loops control cell cycle" Bioinformatics 3 WS 13/14" V 20 " 20"

G1 => S Toggle Switch" Mutual inhibition between Cdk1-CycB and CKI" (cyclin kinase inhibitor)" Tyson et al, Curr. Op. Cell Biol. 15 (2003) 221" Bioinformatics 3 WS 13/14" V 20 " 21"

Mutual Inhibition"???" Assume: CycB:Cdk1:CKI is stable <=> dissociation is very slow" => same topology" <=> same bistable " behavior (?)" Bioinformatics 3 WS 13/14" V 20 " 22"

Rate Equations: Toggle Switch" R2" A" X" R1" R4" R3" Stoichiometric matrix" "(C)" = catalyst" R1" R2" R3" R4" A" 1" S" (C)" R" 1" 1" (C)" E" (C)" 1" 1" EP" 1" 1" X" 1" Bioinformatics 3 WS 13/14" V 20 " 23"

Rate Equations: G1/S Module" R5" R6" R1" R2" R3" R4" R1" R2" R3" R4" R5" R6" CycB" 1" Cdk1" 1" CycB:Cdk1" 1" 1" (C)" 1" CKI" 1" 1" 1" 1" CKI:P3" 1" 1" CKI:P3" 1" CycB:Cdk1:CKI" 1" -1" Bioinformatics 3 WS 13/14" V 20 " 24"

Comparison: Matrices" R2" A" X" R1" R4" R3" R5" R6" R1" R2" R3" R4" R1" R2" R3" R4" A" 1" S" (C)" R" 1" 1" (C)" E" (C)" 1" 1" EP" 1" 1" X" 1" R1" R2" R3" R4" R5" R6" CycB" 1" Cdk1" 1" CycB:Cdk1" 1" 1" (C)" 1" CKI" 1" 1" 1" 1" CKI:P3" 1" 1" CKI:P3" 1" CycB:Cdk1:CKI" 1" -1" Difference: catalysts vs. substrates" Bioinformatics 3 WS 13/14" V 20 " 25"

Comparison: Equations" R2" A" X" R1" R4" R3" R6" R1" R2" R3" R5" R4" Rename species => same rate equations => same behavior" Bioinformatics 3 WS 13/14" V 20 " 26"

Predicted Behavior: G1 => S" Signal: cell growth = concentration of CycB, Cdk1" Response: activity (concentration) of CycB:Cdk1" Toggle switch:" => above critical cell size CycB:Cdk1 activity will switch on" Tyson et al, Curr. Op. Cell Biol. 15 (2003) 221" Bioinformatics 3 WS 13/14" V 20 " 27"

G2 => M" Toggle switch:" mutual activation between " CycB:Cdk1 and Cdc25 " (phosphatase that activates " the dimer)" mutual inhibition between " CycB:Cdk1 and Wee1 " (kinase that inactivates the dimer)" =>"when the cell grows further during the second gap phase G2, the activity of CycB:Cdk1 will increase by a further step! Tyson et al, Curr. Op. Cell Biol. 15 (2003) 221" Bioinformatics 3 WS 13/14" V 20 " 28"

M => G1" Negative feedback loop" oscillator" i) CycB:Cdk1 activates anaphase " promoting complex (APC)" ii) APC activates Cdc20" iii) Cdc20 degrades CycB" Behavior:" at a critical cell size " CycB:Cdk1 activity increases and decreases again" => at low CycB:Cdk1 level, the G1/S toggle switches off again," => cell cycle completed" Tyson et al, Curr. Op. Cell Biol. 15 (2003) 221" Bioinformatics 3 WS 13/14" V 20 " 29"

Overall Behavior" Cell divides at size 1.46" => daughters start " growing from " size 0.73" => switches to " replication at " size 1.25" G1/S toggle => bistability" M/G1 oscillator" G2/M toggle => bistability" Tyson et al, Curr. Op. Cell Biol. 15 (2003) 221" Bioinformatics 3 WS 13/14" V 20 " 30"

Preventing Cross-Talk" Many enzymes are used " in multiple pathways" => how can different signals cross " the same kinase?" => different temporal signature" (slow vs. transient)" => Dynamic modelling!" Bioinformatics 3 WS 13/14" V 20 " 31"

Summary" Today:! Behavior of cell cycle control circuitry from its modules:" two toggle switches + one oscillator" => map biological system onto motif via" stoichiometric matrices" rate equations" Bioinformatics 3 WS 13/14" V 20 " 32"