igem Jamboree MIT Nov 4th 2006 Imperial College

Size: px
Start display at page:

Download "igem Jamboree MIT Nov 4th 2006 Imperial College"

Transcription

1 igem Jamboree MIT Nov 4th 2006

2 Engineering a Molecular Predation Oscillator

3 igem Imperial Biomedical Engineers Electrical Engineer Biochemist Biologists Biomedical Engineers Biochemists Dr Mann

4 Project Ideas Bio-Memory Oscillator Bio-Clock Feasibility Originality of Design BioBrick Availability Future Impact

5 Our Definition What is an Oscillator? Device producing a periodic variation in time of a measurable quantity, e.g. amplitude. Engineering Biology

6 The Engineering Approach Specifications Testing/Validation Design Implementation Standard Engineering Certified Practice Modelling

7 The Main Challenges Fluorescence Main challenges of past oscillators: Unstable Noisy Inflexible Repressilator Requirement for a typical Our Specifications: engineering oscillator Sustained Stability: >10 Oscillations periods High SNR: Signal High to Noise Ratio Controllable Flexibility: Controllable Oscillations Amplitude and Frequency Standardized Device for Modular Design Easy Connectivity Easy Connectivity Time (min.) Figure Reference : Michael B. Elowitz & Stanislas Leibler Nature 2000

8 Our Initial Design Ideas Based on Large populations of molecules to reduce influence of noise Oscillations due to population dynamics A well characterized model Molecular Predator - Prey

9 The Lotka-Volterra Model d d X Y ax Prey Growth cxy Predator Growth Prey bxy Killing by Predator dy Predator Death Population Size X :Prey Y : Predator Time

10 Typical LV Simulations Graph of Prey vs. Time Low Frequency High Frequency prey prey Small Amplitude time time prey prey Large Amplitude time time

11 Required Biochemical Properties d A Self promoted Prey Growth expression of A A Prey Killing by of Predator A by B B A d B Expression of B Predator promoted Growth by AB interaction B A Predator Death of B B Population Size A B Time

12 Molecular System Prey Generator Cell Self promoted expression of A A Cell-cell communication Constraint: Chemostat of A by B A Flexibility: B Ratio of populations A A of B B B Predator Generator Cell Expression of B promoted by AB interaction

13 Quorum sensing/quenching Vibrio fischeri Bacillus + plux LuxI aiia BBa_C0062 Forms a complex with to activate plux BioBricks available BBa_C0061 Makes plux BBa_R0062 plux Promoter BBa_C0160 Degrades ptet BBa_F2620 plux ->Pops Receiver

14 Designing the Prey Generator Required Dynamic Self promoted expression of A Useful BioBricks LuxI C0061 tetr F2620 plux Final Construct LuxI tetr plux LuxI

15 Designing the Predator Generator Required Dynamic Expression of B promoted by AB interaction of A by B of B Useful BioBricks plux R0062 C0062 aiia C0060 Natural degradation Final Construct Lactonase Sensing plux aiia Killing

16 System Overview Prey Generator Cell Pool of will oscillate LuxI ptet plux LuxI Predator Generator Cell Lactonase aiia plux

17 Full System set-up reservoir LuxI tetr plux In-flow Cell population Input signal Change in population ratio Prey molecule generator time Wash-out [] Output signal plux aiia Well mixed culture In chemostat time Predator molecule generator

18 Full System set-up reservoir LuxI [] Output signal tetr plux In-flow Prey molecule generator time Wash-out Cell population plux aiia Well mixed culture In chemostat Change in population ratio time Predator molecule generator

19 Path to Our Goal Start! Specifications Design Modelling Implementation Testing Specifications oscillator Based on Lotka-Volterra Based on Quorum Sensing Modelling Lotka-Volterra Modelling Full system

20 Modelling the Full System d Self promoted expression of of by aiia of d Expression of of d aiia Expression of aiia of aiia

21 Modelling the Full System d [ ] Self promoted expression of of by aiia of d [ ] Expression of of d [ aiia ] aiia Expression of aiia of aiia

22 Modelling the Full System Gene Expression d [ ] a a 0 [ ] + [ ] of by aiia of d [ ] [ ][ ] [ ][ ] Production c of c + 0 of d [ aiia ] aiia [ ][ ] [ ][ ] Production c of c0 + aiia of aiia

23 Modelling the Full System Gene Expression Enzymatic Reaction d [ ] a a 0 [ ] + [ ] b [ aiia][ ] b + [ ] 0 of d [ ] [ ][ ] [ ][ ] Production c of c + 0 of d aiia ] aiia [ [ ][ ] [ ][ ] Production c of c0 + aiia of aiia

24 Modelling the Full System Gene Expression Enzymatic Reaction d [ ] a a 0 [ ] + [ ] b [ aiia][ ] b + [ ] 0 e[ ] of d [ ] [ ][ ] [ ][ ] Production c of c + 0 d 1 [ ] d aiia ] aiia [ [ ][ ] [ ][ ] Production c of c0 + aiia d [ aiia] 2

25 Full System Simulations Graph of Prey vs. Time Low Frequency High Frequency prey prey Small Amplitude time time prey prey Large Amplitude time time

26 Typical System Behaviours Oscillations with limit cycles Predator Predator No oscillations Prey Prey Prey Prey Time Time

27 Modelling the Full System Population dependent d [ ] Gene Expression a a 0 [ ] + [ ] Enzymatic Reaction b [ aiia][ ] b + [ ] 0 e[ ] of d [ ] d aiia ] aiia [ ][ ] [ ][ ] Production c of c + [ [ ][ ] [ ][ ] 0 Production c of c0 + aiia Wash-out related d 1 [ ] d [ aiia] 2 Constant

28 Path to Our Goal Start! Specifications Design Modelling Implementation Testing Specifications oscillator Specifications test constructs Based on Lotka-Volterra Based on Quorum Sensing Design test constructs Modelling Lotka-Volterra Modelling Full system Modelling Test construct Build test constructs Characterized test constructs

29 Breaking Down the Complexity Prey Generator Predator Generator Prey Sensing Predator Sensing Predator Killing PoPs a a0 PoPs c c0 time b b0

30 Characterization Predator Sensing Test part Predictive model transfer function [GFP] Experimental Data Average with variance and curve fit Experimental data []

31 Characterization Predator Sensing Test part Predictive model transfer function [GFP] Average with variance and curve fitting Experimental data Fitting model to data Parameter extractions []

32 Implementation Registry Catalogue Parts Prey J37034 RS+J37034 Prey with Riboswitch + Assembly process J37023 J37024 AiiA Test Construct J37025 Final predator J37033 J37019 Sensing predator J37031 Sensing predator J

33 On Our Experience Specifications Testing/Validation Design Implementation Modelling

34 Path to Our Goal Start! Specifications Design Modelling Implementation Testing Specifications oscillator Specifications test constructs Based on Lotka-Volterra Based on Quorum Sensing Design test constructs Modelling Lotka-Volterra Modelling Full system Modelling Test construct Modelling characterized sys Build test constructs Build full system Characterized test constructs Characterized Full system Our Goal!

35 Contributions to the Registry Functional Parts Built Sequenced Tested Characterized Documented Final Prey J37015 Sensing Prey T9002 Sensing Predator J37016 Cre/Lox J37027 Intermediate Parts Built Sequenced Built Sequenced J37033 J37034 J37019 J37032 J37023

36 Our Wiki Documentation Communication Organization

37 Thank You Acknowledgements: Prof. Tony Cass Dr. Anna Radomska Dr. Rupert Fray Dr. David Leak Dr. Mauricio Barahona Dr. Danny O'Hare Dr. Geoff Baldwin Susan E. Wryter David Featherbe Ciaran Mckeown James Mansfield Funding: European Commission Deputy Rectors Fund Faculty of Engineering Faculty of Natural Sciences

A synthetic oscillatory network of transcriptional regulators

A synthetic oscillatory network of transcriptional regulators A synthetic oscillatory network of transcriptional regulators Michael B. Elowitz & Stanislas Leibler, Nature, 403, 2000 igem Team Heidelberg 2008 Journal Club Andreas Kühne Introduction Networks of interacting

More information

A Synthetic Oscillatory Network of Transcriptional Regulators

A Synthetic Oscillatory Network of Transcriptional Regulators A Synthetic Oscillatory Network of Transcriptional Regulators Michael Elowitz & Stanislas Leibler Nature, 2000 Presented by Khaled A. Rahman Background Genetic Networks Gene X Operator Operator Gene Y

More information

Engineering of Repressilator

Engineering of Repressilator The Utilization of Monte Carlo Techniques to Simulate the Repressilator: A Cyclic Oscillatory System Seetal Erramilli 1,2 and Joel R. Stiles 1,3,4 1 Bioengineering and Bioinformatics Summer Institute,

More information

Introduction. ECE/CS/BioEn 6760 Modeling and Analysis of Biological Networks. Adventures in Synthetic Biology. Synthetic Biology.

Introduction. ECE/CS/BioEn 6760 Modeling and Analysis of Biological Networks. Adventures in Synthetic Biology. Synthetic Biology. Introduction ECE/CS/BioEn 6760 Modeling and Analysis of Biological Networks Chris J. Myers Lecture 17: Genetic Circuit Design Electrical engineering principles can be applied to understand the behavior

More information

Basic Synthetic Biology circuits

Basic 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 information

Further analysis. Objective Based on the understanding of the LV model from last report, further study the behaviours of the LV equations parameters.

Further analysis. Objective Based on the understanding of the LV model from last report, further study the behaviours of the LV equations parameters. Further analysis Objective Based on the understanding of the LV model from last report, further study the behaviours of the LV equations parameters. Changes have been made Method of finding the mean is

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION doi:10.1038/nature13238 Additional Experimental Results Degradation tag experiments In addition to exploring the effect of variable-length linker (TS repeats) on the phase-shift in module degradation (Supplementary

More information

Non-Linear Models. Non-Linear Models Cont d

Non-Linear Models. Non-Linear Models Cont d Focus on more sophistiated interaction models between systems. These lead to non-linear, rather than linear, DEs; often not soluble exactly in analytical form so use Phase-Plane Analysis. This is a method

More information

Getting Started With The Predator - Prey Model: Nullclines

Getting Started With The Predator - Prey Model: Nullclines Getting Started With The Predator - Prey Model: Nullclines James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline The Predator

More information

2D-Volterra-Lotka Modeling For 2 Species

2D-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 information

Lecture 1. Scott Pauls 1 3/28/07. Dartmouth College. Math 23, Spring Scott Pauls. Administrivia. Today s material.

Lecture 1. Scott Pauls 1 3/28/07. Dartmouth College. Math 23, Spring Scott Pauls. Administrivia. Today s material. Lecture 1 1 1 Department of Mathematics Dartmouth College 3/28/07 Outline Course Overview http://www.math.dartmouth.edu/~m23s07 Matlab Ordinary differential equations Definition An ordinary differential

More information

Ordinary Differential Equations

Ordinary 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 information

Modelling biological oscillations

Modelling 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 information

Lab 5: Nonlinear Systems

Lab 5: Nonlinear Systems Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting

More information

the noisy gene Biology of the Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB)

the noisy gene Biology of the Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB) Biology of the the noisy gene Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB) day III: noisy bacteria - Regulation of noise (B. subtilis) - Intrinsic/Extrinsic

More information

A New Chassis for Synthetic Biology: Bacteria Without a Cell Wall. L-forms

A New Chassis for Synthetic Biology: Bacteria Without a Cell Wall. L-forms A New Chassis for Synthetic Biology: Bacteria Without a Cell Wall L-forms Pros & Cons of Cell Wall Cell wall Cell membrane Cell membrane DNA ribosomes RNA metabolites Bacterium with cell wall Bacterium

More information

Analog Electronics Mimic Genetic Biochemical Reactions in Living Cells

Analog Electronics Mimic Genetic Biochemical Reactions in Living Cells Analog Electronics Mimic Genetic Biochemical Reactions in Living Cells Dr. Ramez Daniel Laboratory of Synthetic Biology & Bioelectronics (LSB 2 ) Biomedical Engineering, Technion May 9, 2016 Cytomorphic

More information

Is chaos possible in 1d? - yes - no - I don t know. What is the long term behavior for the following system if x(0) = π/2?

Is chaos possible in 1d? - yes - no - I don t know. What is the long term behavior for the following system if x(0) = π/2? Is chaos possible in 1d? - yes - no - I don t know What is the long term behavior for the following system if x(0) = π/2? In the insect outbreak problem, what kind of bifurcation occurs at fixed value

More information

MA 138: Calculus II for the Life Sciences

MA 138: Calculus II for the Life Sciences MA 138: Calculus II for the Life Sciences David Murrugarra Department of Mathematics, University of Kentucky. Spring 2016 David Murrugarra (University of Kentucky) MA 138: Section 11.4.2 Spring 2016 1

More information

RiboTALE: A modular, inducible system for accurate gene expression control

RiboTALE: A modular, inducible system for accurate gene expression control 8 9 10 11 1 1 1 1 1 1 18 19 0 1 RiboTALE: A modular, inducible system for accurate gene expression control Navneet Rai, 1 Aura Ferreiro, Alexander Neckelmann, Amy Soon, Andrew Yao, 1, Justin Siegel, 1,,

More information

CPS 5310: Parameter Estimation. Natasha Sharma, Ph.D.

CPS 5310: Parameter Estimation. Natasha Sharma, Ph.D. Example Suppose our task is to determine the net income for year 2019 based on the net incomes given below Year Net Income 2016 48.3 million 2017 90.4 million 2018 249.9 million Last lecture we tried to

More information

Dynamical Systems in Biology

Dynamical Systems in Biology Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems

More information

First Order Systems of Linear Equations. or ODEs of Arbitrary Order

First Order Systems of Linear Equations. or ODEs of Arbitrary Order First Order Systems of Linear Equations or ODEs of Arbitrary Order Systems of Equations Relate Quantities Examples Predator-Prey Relationships r 0 = r (100 f) f 0 = f (r 50) (Lokta-Volterra Model) Systems

More information

igem 2009 Team Newcastle

igem 2009 Team Newcastle igem 2009 Team Newcastle Introduction Environmental project Heavy metal pollution in soil Cadmium accumulation issue Image: http://www.enst.umd.edu/people/weil/researchprojects.cfm Use of engineered micro-organisms

More information

MA 138 Calculus 2 with Life Science Applications Handout

MA 138 Calculus 2 with Life Science Applications Handout .. MA 138 Calculus 2 with Life Science Applications Handout Alberto Corso alberto.corso@uky.edu Department of Mathematics University of Kentucky February 17, 2017 . Example 4 (Lotka-Volterra Predator-Prey

More information

2. Mathematical descriptions. (i) the master equation (ii) Langevin theory. 3. Single cell measurements

2. Mathematical descriptions. (i) the master equation (ii) Langevin theory. 3. Single cell measurements 1. Why stochastic?. Mathematical descriptions (i) the master equation (ii) Langevin theory 3. Single cell measurements 4. Consequences Any chemical reaction is stochastic. k P d φ dp dt = k d P deterministic

More information

A NUMERICAL STUDY ON PREDATOR PREY MODEL

A NUMERICAL STUDY ON PREDATOR PREY MODEL International Conference Mathematical and Computational Biology 2011 International Journal of Modern Physics: Conference Series Vol. 9 (2012) 347 353 World Scientific Publishing Company DOI: 10.1142/S2010194512005417

More information

Simulation of Chemical Reactions

Simulation of Chemical Reactions Simulation of Chemical Reactions Cameron Finucane & Alejandro Ribeiro Dept. of Electrical and Systems Engineering University of Pennsylvania aribeiro@seas.upenn.edu http://www.seas.upenn.edu/users/~aribeiro/

More information

BIOCHEMICAL OSCILLATIONS VIA STOCHASTIC PROCESSES

BIOCHEMICAL OSCILLATIONS VIA STOCHASTIC PROCESSES BIOCHEMICAL OSCILLATIONS VIA STOCHASTIC PROCESSES EDUARDO S. ZERON 1. Historical Introduction Since the first observation of current oscillations during electrodissolution of an iron wire in nitric acid

More information

Mathematization in Synthetic Biology: Analogies, Templates, and Fictions

Mathematization in Synthetic Biology: Analogies, Templates, and Fictions Mathematization in Synthetic Biology: Analogies, Templates, and Fictions Tarja Knuuttila and Andrea Loettgers T. Knuuttila University of South Carolina, University of Helsinki 901 Sumter St., Byrnes Suite,

More information

Quorum sensing in plantassociated. S. Brook Peterson Parsek Lab UW Microbiology

Quorum sensing in plantassociated. S. Brook Peterson Parsek Lab UW Microbiology Quorum sensing in plantassociated bacteria S. Brook Peterson Parsek Lab UW Microbiology Outline What is quorum sensing? QS in plant associated bacteria What traits are regulated by QS? What benefits does

More information

Lesson 9: Predator-Prey and ode45

Lesson 9: Predator-Prey and ode45 Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and

More information

BIOS 3010: ECOLOGY. Dr Stephen Malcolm. Laboratory 6: Lotka-Volterra, the logistic. equation & Isle Royale

BIOS 3010: ECOLOGY. Dr Stephen Malcolm. Laboratory 6: Lotka-Volterra, the logistic. equation & Isle Royale BIOS 3010: ECOLOGY Dr Stephen Malcolm Laboratory 6: Lotka-Volterra, the logistic equation & Isle Royale This is a computer-based activity using Populus software (P), followed by EcoBeaker analyses of moose

More information

EAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science

EAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Transient Response of a Chemical Reactor Concentration of a substance in a chemical reactor

More information

NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD CHAIN WITH LOTKA-VOLTERRA LINEAR FUNCTIONAL RESPONSE

NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD CHAIN WITH LOTKA-VOLTERRA LINEAR FUNCTIONAL RESPONSE Journal of Sustainability Science and Management Volume 6 Number 1, June 2011: 44-50 ISSN: 1823-8556 Universiti Malaysia Terengganu Publisher NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD

More information

MAS212 Scientific Computing and Simulation

MAS212 Scientific Computing and Simulation MAS212 Scientific Computing and Simulation Dr. Sam Dolan School of Mathematics and Statistics, University of Sheffield Autumn 2017 http://sam-dolan.staff.shef.ac.uk/mas212/ G18 Hicks Building s.dolan@sheffield.ac.uk

More information

Gerardo Zavala. Math 388. Predator-Prey Models

Gerardo Zavala. Math 388. Predator-Prey Models Gerardo Zavala Math 388 Predator-Prey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked

More information

Outline. Calculus for the Life Sciences. What is a Differential Equation? Introduction. Lecture Notes Introduction to Differential Equa

Outline. Calculus for the Life Sciences. What is a Differential Equation? Introduction. Lecture Notes Introduction to Differential Equa Outline Calculus for the Life Sciences Lecture Notes to Differential Equations Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu 1 Department of Mathematics and Statistics Dynamical Systems Group Computational

More information

Cadmium Poisoning Ingestion of cadmium can cause: Cancer Renal failure Bone weakening and fracture Damage to nervous and immune systems Death Danger Level Present in over 60% of hazardous waste sites as

More information

Lecture 1 Modeling in Biology: an introduction

Lecture 1 Modeling in Biology: an introduction Lecture 1 in Biology: an introduction Luca Bortolussi 1 Alberto Policriti 2 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste Via Valerio 12/a, 34100 Trieste. luca@dmi.units.it

More information

7.32/7.81J/8.591J: Systems Biology. Fall Exam #1

7.32/7.81J/8.591J: Systems Biology. Fall Exam #1 7.32/7.81J/8.591J: Systems Biology Fall 2013 Exam #1 Instructions 1) Please do not open exam until instructed to do so. 2) This exam is closed- book and closed- notes. 3) Please do all problems. 4) Use

More information

Lokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt

Lokta-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 information

Motivation and Goals. Modelling with ODEs. Continuous Processes. Ordinary Differential Equations. dy = dt

Motivation and Goals. Modelling with ODEs. Continuous Processes. Ordinary Differential Equations. dy = dt Motivation and Goals Modelling with ODEs 24.10.01 Motivation: Ordinary Differential Equations (ODEs) are very important in all branches of Science and Engineering ODEs form the basis for the simulation

More information

Workshop on Theoretical Ecology and Global Change March 2009

Workshop on Theoretical Ecology and Global Change March 2009 2022-3 Workshop on Theoretical Ecology and Global Change 2-18 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological

More information

APPM 2360 Lab #3: The Predator Prey Model

APPM 2360 Lab #3: The Predator Prey Model APPM 2360 Lab #3: The Predator Prey Model 1 Instructions Labs may be done in groups of 3 or less. One report must be turned in for each group and must be in PDF format. Labs must include each student s:

More information

1 The pendulum equation

1 The pendulum equation Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating

More information

Predator-Prey Population Models

Predator-Prey Population Models 21 Predator-Prey Population Models Tools Used in Lab 21 Hudson Bay Data (Hare- Lynx) Lotka-Volterra Lotka-Volterra with Harvest How can we model the interaction between a species of predators and their

More information

Predator-Prey Population Dynamics

Predator-Prey Population Dynamics Predator-Prey Population Dynamics Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ October 2,

More information

Biomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems

Biomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems Biomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems Processing information to coordinate activity regulatory & signaling networks

More information

1 Introduction. Maria Luisa Guerriero. Jane Hillston

1 Introduction. Maria Luisa Guerriero. Jane Hillston Time-Series Analysis of Circadian Clocks: Quantifying the Effect of Stochasticity on the Distribution of Phase, Period and Amplitude of Circadian Oscillations Maria Luisa Guerriero Centre for Systems Biology

More information

Lecture 4: Importance of Noise and Fluctuations

Lecture 4: Importance of Noise and Fluctuations Lecture 4: Importance of Noise and Fluctuations Jordi Soriano Fradera Dept. Física de la Matèria Condensada, Universitat de Barcelona UB Institute of Complex Systems September 2016 1. Noise in biological

More information

A Discrete Model of Three Species Prey- Predator System

A Discrete Model of Three Species Prey- Predator System ISSN(Online): 39-8753 ISSN (Print): 347-670 (An ISO 397: 007 Certified Organization) Vol. 4, Issue, January 05 A Discrete Model of Three Species Prey- Predator System A.George Maria Selvam, R.Janagaraj

More information

Perfect synthetic biology project

Perfect synthetic biology project Perfect synthetic biology project Something beautiful! Idea Community project Present Evaluate Improve Final product Our problem: Algal Blooms Explosion in the growth of bluegreen algae (cyanobacteria)

More information

Build to understand: synthetic approaches to biology

Build to understand: synthetic approaches to biology REVIEW ARTICLE View Article Online View Journal View Issue Cite this: Integr. Biol., 2016, 8, 394 Received 8th October 2015, Accepted 7th December 2015 DOI: 10.1039/c5ib00252d www.rsc.org/ibiology Build

More information

ANALYTIC SOLUTIONS TO A CLASS OF TWO-DIMENSIONAL LOTKA-VOLTERRA DYNAMICAL SYSTEMS

ANALYTIC SOLUTIONS TO A CLASS OF TWO-DIMENSIONAL LOTKA-VOLTERRA DYNAMICAL SYSTEMS PROCEEDINGS OF THE INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 18 21, 2000, Atlanta, USA pp. 137 142 ANALYTIC SOLUTIONS TO A CLASS OF TWO-DIMENSIONAL LOTKA-VOLTERRA DYNAMICAL

More information

Part 3: Introduction to Master Equation and Complex Initial Conditions in Lattice Microbes

Part 3: Introduction to Master Equation and Complex Initial Conditions in Lattice Microbes Part 3: Introduction to Master Equation Cells: and Complex Initial Conditions in Lattice re cells Microbes en Biophysics, and UC urgh, June 6-8, 2016 rson Joseph R. Peterson and Michael J. Hallock Luthey-Schulten

More information

Section 5. Graphing Systems

Section 5. Graphing Systems Section 5. Graphing Systems 5A. The Phase Plane 5A-1. Find the critical points of each of the following non-linear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A-2. Write each

More information

Computer Simulation and Data Analysis in Molecular Biology and Biophysics

Computer 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 information

Cellular Biophysics SS Prof. Manfred Radmacher

Cellular Biophysics SS Prof. Manfred Radmacher SS 20007 Manfred Radmacher Ch. 12 Systems Biology Let's recall chemotaxis in Dictiostelium synthesis of camp excretion of camp external camp gradient detection cell polarity cell migration 2 Single cells

More information

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

Written Exam 15 December Course name: Introduction to Systems Biology Course no Technical University of Denmark Written Exam 15 December 2008 Course name: Introduction to Systems Biology Course no. 27041 Aids allowed: Open book exam Provide your answers and calculations on separate

More information

BIOL 410 Population and Community Ecology. Predation

BIOL 410 Population and Community Ecology. Predation BIOL 410 Population and Community Ecology Predation Intraguild Predation Occurs when one species not only competes with its heterospecific guild member, but also occasionally preys upon it Species 1 Competitor

More information

Perfect synthetic biology project

Perfect synthetic biology project Perfect synthetic biology project Something beautiful! Idea Community project Present Evaluate Improve Final product Our problem: Algal Blooms Explosion in the growth of bluegreen algae (cyanobacteria)

More information

Lecture 20/Lab 21: Systems of Nonlinear ODEs

Lecture 20/Lab 21: Systems of Nonlinear ODEs Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species

More information

Dynamical Systems and Chaos Part II: Biology Applications. Lecture 10: Coupled Systems. Ilya Potapov Mathematics Department, TUT Room TD325

Dynamical Systems and Chaos Part II: Biology Applications. Lecture 10: Coupled Systems. Ilya Potapov Mathematics Department, TUT Room TD325 Dynamical Systems and Chaos Part II: Biology Applications Lecture 10: Coupled Systems. Ilya Potapov Mathematics Department, TUT Room TD325 Foreword In order to model populations of physical/biological

More information

second order Runge-Kutta time scheme is a good compromise for solving ODEs unstable for oscillators

second order Runge-Kutta time scheme is a good compromise for solving ODEs unstable for oscillators ODE Examples We have seen so far that the second order Runge-Kutta time scheme is a good compromise for solving ODEs with a good precision, without making the calculations too heavy! It is unstable for

More information

Correcting Errors Due to Species Correlations in the Marginal Probability Density Evolution

Correcting Errors Due to Species Correlations in the Marginal Probability Density Evolution Utah State University DigitalCommons@USU All Graduate Theses and Dissertations Graduate Studies 5-2013 Correcting Errors Due to Species Correlations in the Marginal Probability Density Evolution Abiezer

More information

Network Analysis of Biochemical Reactions in Complex Environments

Network 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 information

Alternatives to competition. Lecture 13. Facilitation. Functional types of consumers. Stress Gradient Hypothesis

Alternatives to competition. Lecture 13. Facilitation. Functional types of consumers. Stress Gradient Hypothesis Lecture 13 Finishing Competition and Facilitation Consumer-Resource interactions Predator-prey population dynamics Do predators regulate prey? Lotka-Volterra predator-prey model Predator behavior matters:

More information

Systems Biology in Photosynthesis INTRODUCTION

Systems 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

Journal of the Vol. 36, pp , 2017 Nigerian Mathematical Society c Nigerian Mathematical Society

Journal of the Vol. 36, pp , 2017 Nigerian Mathematical Society c Nigerian Mathematical Society Journal of the Vol. 36, pp. 47-54, 2017 Nigerian Mathematical Society c Nigerian Mathematical Society A CLASS OF GENERALIZATIONS OF THE LOTKA-VOLTERRA PREDATOR-PREY EQUATIONS HAVING EXACTLY SOLUBLE SOLUTIONS

More information

ROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM

ROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 6, Number 1, Winter 1997 ROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM J. CHATTOPADHYAY, E. BERETTA AND F. SOLIMANO ABSTRACT. The present paper deals with

More information

Instabilities In A Reaction Diffusion Model: Spatially Homogeneous And Distributed Systems

Instabilities In A Reaction Diffusion Model: Spatially Homogeneous And Distributed Systems Applied Mathematics E-Notes, 10(010), 136-146 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Instabilities In A Reaction Diffusion Model: Spatially Homogeneous And

More information

MATH 250 Homework 4: Due May 4, 2017

MATH 250 Homework 4: Due May 4, 2017 Due May 4, 17 Answer the following questions to the best of your ability. Solutions should be typed. Any plots or graphs should be included with the question (please include the questions in your solutions).

More information

Phenomenon: Canadian lynx and snowshoe hares

Phenomenon: Canadian lynx and snowshoe hares Outline Outline of Topics Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Phenomenon: Canadian lynx and snowshoe hares All began with

More information

dv dt Predator-Prey Models

dv dt Predator-Prey Models Predator-Prey Models This is a diverse area that includes general models of consumption: Granivores eating seeds Parasitoids Parasite-host interactions Lotka-Voterra model prey and predator: V = victim

More information

Population Dynamics II

Population Dynamics II Population Dynamics II In this class, we shall analyze behavioral patterns of ecosystems, in which more than two species interact with each other. Such systems frequently exhibit chaotic behavior. Chaotic

More information

The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor

The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor , pp. 35-46 http://dx.doi.org/10.14257/ijbsbt.2017.9.3.04 The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor Alemu Geleta Wedajo

More information

Nature Genetics: doi: /ng Supplementary Figure 1. Plasmid diagrams of primary strains used in this work.

Nature Genetics: doi: /ng Supplementary Figure 1. Plasmid diagrams of primary strains used in this work. Supplementary Figure 1 Plasmid diagrams of primary strains used in this work. Refer to Supplementary Table 1 for additional information and full list of strains. Diagrams are not to scale. (a) Strain used

More information

(Refer Slide Time: 00:32)

(Refer Slide Time: 00:32) Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra

More information

CS-E5880 Modeling biological networks Gene regulatory networks

CS-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 information

Predation. Vine snake eating a young iguana, Panama. Vertebrate predators: lions and jaguars

Predation. Vine snake eating a young iguana, Panama. Vertebrate predators: lions and jaguars Predation Vine snake eating a young iguana, Panama Vertebrate predators: lions and jaguars 1 Most predators are insects Parasitoids lay eggs in their hosts, and the larvae consume the host from the inside,

More information

damped oscillation Figure 11.3

damped oscillation Figure 11.3 Chpater X How to build a Biological oscillator Oscillations are thrilling hearts beat, cells divide every cell cycle, circadian clocks keep, neurons click in trains of regularly spaced spikes. Biological

More information

Lotka Volterra Predator-Prey Model with a Predating Scavenger

Lotka Volterra Predator-Prey Model with a Predating Scavenger Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics

More information

Population Ecology. Chapter 44

Population Ecology. Chapter 44 Population Ecology Chapter 44 Stages of Biology O Ecology is the interactions of organisms with other organisms and with their environments O These interactions occur in different hierarchies O The simplest

More information

8 Ecosystem stability

8 Ecosystem stability 8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single

More information

Chapter 2 Biochemical Reaction Networks

Chapter 2 Biochemical Reaction Networks Chapter 2 Biochemical Reaction Networks This chapter is a basic introduction to chemical reactions and chemical species. Different ways of quantifying the abundance of molecules lead to the notions of

More information

The major interactions between two organisms in a mixed culture are : Competition Neutralism Mutualism Commensalism Amensalism Prey-predator

The major interactions between two organisms in a mixed culture are : Competition Neutralism Mutualism Commensalism Amensalism Prey-predator 1 Introduction Major classes of interaction in mixed cultures Simple models describing mixed-cultures interactions Mixed culture in nature Industrial utilization of mixed cultures 2 1 The dynamic of mixed

More information

Biology X, Fall 2010

Biology X, Fall 2010 Biology X, Fall 2010 Signaling games Andreas Witzel CIMS, NYU 2010-11-01 see Skyrms, Signals and Gintis, Game theory evolving Vervet monkeys Vervet monkeys have distinct alarm calls for different predators

More information

Field experiments on competition. Field experiments on competition. Field experiments on competition

Field experiments on competition. Field experiments on competition. Field experiments on competition INTERACTIONS BETWEEN SPECIES Type of interaction species 1 species 2 competition consumer-resource (pred, herb, para) mutualism detritivore-detritus (food is dead) Field experiments on competition Example

More information

Math 128A Spring 2003 Week 12 Solutions

Math 128A Spring 2003 Week 12 Solutions Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems

More information

Interactions. Yuan Gao. Spring Applied Mathematics University of Washington

Interactions. Yuan Gao. Spring Applied Mathematics University of Washington Interactions Yuan Gao Applied Mathematics University of Washington yuangao@uw.edu Spring 2015 1 / 27 Nonlinear System Consider the following coupled ODEs: dx = f (x, y). dt dy = g(x, y). dt In general,

More information

Handbook of Stochastic Methods

Handbook of Stochastic Methods C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical

More information

Graded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo

Graded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo 2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical

More information

1. Population dynamics of rabbits and foxes

1. Population dynamics of rabbits and foxes 1. Population dynamics of rabbits and foxes (a) A simple Lotka Volterra Model We have discussed in detail the Lotka Volterra model for predator-prey relationships dn prey dt = +R prey,o N prey (t) γn prey

More information

Introduction to Bioinformatics

Introduction to Bioinformatics Systems biology Introduction to Bioinformatics Systems biology: modeling biological p Study of whole biological systems p Wholeness : Organization of dynamic interactions Different behaviour of the individual

More information

Mathematical Modelling Lecture 10 Difference Equations

Mathematical Modelling Lecture 10 Difference Equations Lecture 10 Difference Equations phil.hasnip@york.ac.uk Overview of Course Model construction dimensional analysis Experimental input fitting Finding a best answer optimisation Tools for constructing and

More information

MATH 1014 N 3.0 W2015 APPLIED CALCULUS II - SECTION P. Perhaps the most important of all the applications of calculus is to differential equations.

MATH 1014 N 3.0 W2015 APPLIED CALCULUS II - SECTION P. Perhaps the most important of all the applications of calculus is to differential equations. MATH 1014 N 3.0 W2015 APPLIED CALCULUS II - SECTION P Stewart Chapter 9 Differential Equations Perhaps the most important of all the applications of calculus is to differential equations. 9.1 Modeling

More information

THETA-LOGISTIC PREDATOR PREY

THETA-LOGISTIC PREDATOR PREY THETA-LOGISTIC PREDATOR PREY What are the assumptions of this model? 1.) Functional responses are non-linear. Functional response refers to a change in the rate of exploitation of prey by an individual

More information

Summary of the ISBI 2013 Grand Challenge on 3D Deconvolution Microscopy

Summary of the ISBI 2013 Grand Challenge on 3D Deconvolution Microscopy Summary of the ISBI 23 Grand Challenge on 3D Deconvolution Microscopy Cédric Vonesch & Stamatios Lefkimmiatis Biomedical Imaging Group EPFL, Lausanne, Switzerland April 7, 23 Thanks to our sponsors! ISBI

More information

A Producer-Consumer Model With Stoichiometry

A Producer-Consumer Model With Stoichiometry A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann

More information