A Population-level Hybrid Model of Tumour-Immune System Interplay: model construction and analysis.
|
|
- Robyn Golden
- 5 years ago
- Views:
Transcription
1 A Population-level Hybrid Model of Tumour-Immune System Interplay: model construction and analysis. Giulio Caravagna Dipartimento di Informatica, Sistemistica e Comunicazione, Università degli Studi Milano-Bicocca. giulio.caravagna@disco.unimib.it
2 Joint works on T-IS interplay with: A. d Onofrio (European Institute of Oncology); R. Barbuti and P. Milazzo (University of Pisa). References: G.Caravagna, A.d Onofrio, P.Milazzo, R.Barbuti (21) Antitumour Immune Surveillance Through Stochastic Oscillations. Journal of Theoretical Biology, 265 (3), , 21. G.Caravagna, R.Barbuti, A.d Onofrio, (211) Fine-tuning anti-tumor immunotherapies via stochastic simulations. Submitted to BMC Bioinformatics.
3 Step 1: the bio-bakground, as understood by a computer scientist.
4 IS-response triggered by specific (neo)antigens; antigens appear as a consequence of (epi)genetic events; tumours evolve to eventually evade the IS-control;
5 ( [..] naive T cells express neither BLIMP1 (1) nor interleukin-2 (IL-2). Following encounter with antigen, activated T cells secrete high quantities of IL-2 and upregulate the IL-2 receptor (IL-2R). Engagement of IL-2R by IL-2 results in BLIMP1 transcription in early effector cells. BLIMP1 expression increases in late effector cells and represses Il2 expression through an unknown mechanism. [..] ) (1) The B-lymphocyte-induced maturation protein 1 (BLIMP1) guides the fate of effector B and T cells.
6 Immune surveillance: a phenomenon concerning the control, and sometimes the eradication, of a tumour by the IS. (... is a theory that the IS patrols the body not only to recognize and destroy invading pathogens but also host cells that become cancerous. Perhaps potential cancer cells arise frequently throughout life, but the IS usually destroys them as fast as they appear. There is some evidence for this attractive notion. There is also evidence that the IS mounts an attack against established cancers although it often fails. ) Old hypothesis ( 19) recently supported by: new molecular techniques; many epidemiological studies in cancer immunobiology. Take-home lesson: the tumour may be suppressed by the IS. Goal 1: population-level model of T-IS interplay + surveillance.
7 Step 2: a purely deterministic model.
8 Some basic background: the logistic growth function is the ODE dy (1 dx = ry y ) C where: y is the population size, x is time; r defines the growth rate; C is the carrying capacity of the organism.
9 Similarly, the Michaelis-menten kinetics is the ODE modeling the enzyme kinetics of E + S k1 ES k2 E + P with reversible de-comeplexation ES k 1 1 E + S.
10 Building a model: the choice of the abstraction level is crucial. Possible agents: microscopic: molecules, nucleotides, chromosome, proteins,... macroscopic: cells, proteins, tissues,... Consequently involved events: microscopic: transcriptions, complexations, dimerization... macroscopic: apoptosis, coagulation,... This choice is crucial and affects the type of data necessary to analyze the model; the validity of the predictions model.
11 We define a macroscopic population-level model. Agents as populations: people; animals; cells;.. Consequently involved events: people move, get sick, become virus carriers, vaccinate, recover,... cells proliferate, duplicate, cross membranes, die,... Example SIR models;...
12 The deterministic population-level model that we consider appears in: Kirschner and Panetta, Journal of Mathematical Biology 37(3), Why? It has been highly influential (i.e. 14 citations). Cellular populations involved: Tumour cells (T) and Effector cells (E).... plus a protein-level population: Interleukins IL-2 (I). Mathematically, each population is modeled by a density, and its infinitesimal changes are ruled by an ODE. This results in a system of 3 linked ODEs.
13 Tumour cells (T): (i) saturated logistic growth and (ii) cell death induced by the effectors via Michaelis-menten kinetics. T = rt (1 bt ) ate g 2 + T Effector cells (E): (i) linear growth (ct, c antigenicity) models the response, (ii) Michaelis-menten proliferation induced by interleukins and (iii) death as linear apoptosis. E = p 1IE g 1 + I µ 2E + ct Interleukins IL-2 (I): (i) Michaelis-menten proliferation induced by the T-IS interplay and (ii) linear degradation. I = p 2TE g 3 + T µ 3I
14 To analyze a model you need data... Param. Val. Unit Description r.18 days 1 baseline growth rate of the tumor b 1 9 ml 1 carrying capacity of the tumor a 1 ml days 1 baseline strength of the killing rate by effectors c 1 4 days 1 tumor antigenicity g T 1 5 ml 1 5% reduction factor of the killing rate by effectors g E pg l 1 5% reduction factor of IL-stimulated growth rate of effectors g I 1 3 ml 1 5% reduction factor of production rate of ILs p E.1245 days 1 baseline strength of the IL-stimulated growth rate of effectors p I 5 pg days 1 baseline strength of production rate of ILs µ E.3 days 1 inverse of average lifespan of effectors µ I 1 days 1 loss/degradation rate of IL 2... these are fitted on data pertaining to mice.
15 Analysis varying c, i.e. the average number of antigens expressed by a tumour cells.
16 Fact: the equilibrium (,, ) is unstable. (A) bounded growth (i.e. death); (B,C) late-recurrent oscillations (i.e. undetectable miss is different from eradication since tumours eventually evade); (D) dumped oscillations. KP precludes eradication, but small numbers are predicted (T 1).
17 Vexata quæstio: what if we consider (intrinsic) stochastic effects? Step 2: towards hybrid models...
18 Wee need to move from densities to explicit numbers of cells or molecules. We can do that by using the volume V. We estimated volume V [1.16, 3.2]ml for a chimeric mouse: - body weight 2gr (female) up to 4gr grams (male); - blood volume [5.8, 8]ml per 1gr. With V and the Avogadro number we can convert the KP model to a stochastic Gillespie-based model.
19 Is a Gillespie-based model feasible? molecular weight of interleukins is Da; the average number of proteins is huge (#I ); You should know that if A then 1/A. What happens with Gillespie s equation τ = 1 aj (x) log(r 1 ) when #I? A purely stochastic model is computationally intractable. Moreover, when numbers are big stochastic fluctuations are neglectable, and the ODE-approach is not that bad...
20 Deterministic models: densities/concentrations; changes yield differential equations (ODE, DDE, PDE); fluid variations. What is in the middle? Stochastic models: exact numbers; probability distributions; jump processes (piece-wise constant).
21 Deterministic models: densities/concentrations; changes yield differential equations (ODE, DDE, PDE); fluid variations. Hybrid models: exact numbers + densities/concentrations; probability distributions + differential equations; jump processes (piece-wise constant) + fluid variations. Stochastic models: exact numbers; probability distributions; jump processes (piece-wise constant).
22 ([..] we may assume that the dynamics of I (t) is well approximated by a linear ODE with randomly varying coefficients which, however, in the intervals between two consecutive stochastic events are evidently constant [..] ) Hybrid model : with discrete E and T, and continuous I : (i) Reaction-based events (a 4 is time-dependent): T T + 1 a 1 = r 2 T T T 1 a 2 = r 2 bv 1 T 2 T T 1 a 3 = (ate)/(g T V + T ) E E + 1 a 4 (t) = (p E EI (t))/(g1 + I (t)) E E 1 a 5 = µ E E E E + 1 a 6 = ct (ii) interlinked with the ODE for I (t). Next steps to achieve Goal 1: 1. find a way to simulate this model; 2. reproduce all the KP results; 3. investigate whether this model exhibits eradication.
23 Gillespie s SSA states that, if p(r j fires at t + dt x, t) = a j (x)dt then the putative time follows p(τ x, t) = a (x) exp( a (x)τ) where a (x)τ = j t+τ t a j (x)dt. This because the state is constant in [t, t + τ), hence t+τ t t+τ a j (x)dt = a j (x) dt = a j (x)τ. t For time-dependent propensity functions it becomes t+τ t a j (t)dt a function of time yielding a unknown density.
24 E and T constant in any stochastic interval [t i ; t i+1 ); Without immunotherapies I (t) has analytical form t n t < t n+1 I (t) = B n + (I (t n ) B n ) e µ3(t tn) where B n is a constant function of (T n, E n ). t n+1 = t n + τ obtained by generalizing Gillespie: let χ Exp(1) in tn+τ t n a 4 (t)dt + τ 7 a j a j = χ and P(R j ) = 7 i=1 a i i=1 i 4 Transcendental equation no analytical form and costly simulation.
25 Why is it costly: let us assume that we want to find a of a function f, i.e. x. f (x) =? Many techniques, we used the bisection method: find a such that f (a) < ; find b such that f (b) > ; pick the mid-point ( ) b a f 2 and if it is smaller than then assign b a to a, otherwise to b; 2 ( ) b a repeat until f is arbitrary close to. 2
26 The hybrid model is qualitatively equivalent to the KP e+8.8 4e e+8.6 3e+8.5 T 2.5e+8 2e+8 density.4 1.5e+8.3 1e+8.2 5e time (days) eradication time T () = 1 E() = I () = c =.2..however, (left) low-level oscillations yield eradication: T ( 14) =. Probabilistic analysis: what is the empirical PDF of eradication? time-dependent property (i.e. model-checking a path); For PDFs (right) we averaged 1 3 runs; PDFs can not be obtained by KP.
27 PDF of eradication with variable E()/I () and fixed antigenicity density.5.4 density eradication time eradication time T () = 1 I () = E() = 1 i (left) i = 1, 2, 3 (right) i = 4, density eradication time T () = 1 E() = I () = 1 i i = 1, 2,..., 6
28 Part I: (basic) model construction and analysis. Part II: further model modifications and analysis.
29 Immunotherapy: therapeutical-driven modication of the T-IS interplay. Therapy-related patient-dependent issues: 1. How should the IS be stimulated? 2. Which is a good scheduling for a therapy? 3. How much should be the duration of a single session of therapy? Known facts, results of experimental/theoretical studies: fixing such parameters is non-trivial; scheduling plays a major role in the effectiveness of a therapy. Two somewhat contrasting needs arise: avoiding that the therapy-related side-effects kill the patients; no tumour cells at the end of the therapy (i.e eradication). Goal 2: fine-tuning immunotherapies by means of stochastic simulations.
30 We study (i) passive and (ii) active immunotherapies in the hybrid model. (i) Adoptive cellular immunotherapy (ACI) boost directly E i.e. withdraw white blood cells, culture them, re-inject. E = p E IE g E + I µ E E + ct +V σ E (t) (ii) IL-2-based immunotherapy: boost directly I, indirectly E i.e. inject interleukins in the patient. I = p I V TE g I V + T µ I I +σ I (t) We consider piecewise-constant/impulsive and daily/weekly deliveries. d d 1 d k t t s θ θ + A θ 1 θ 1 + A θ k θ k + A t e In both cases the total injected drug is the same.
31 Mathematical modeling of therapies is pretty intuitive: for t s t t e piecewise-constant: σ(t) = k i= d i[ H(t θi ) H(A t + θ i ) ]. impulsive: σ(t) = k i= u iδ(t θ i ). this yields impulsive ODEs (i.e. x(θ + i ) = x(θ i ) + u i ). The simulation algorithm is slightly modified (2 time-dependent react.). Simulation time spans from mins to days. (v.s. few mins for KP). Asymptotic analysis: ODEs predicts eradication only if therapy. i.e. when the therapy stops, tumour restarts growing!
32 IL-2 therapy: {1,...,3 }. T () = 1 5, E() = I () =. 2e+6 T E 2e+6 2 T E e+6 1.5e e e Single-run of piecewise-constant (left) and impulsive (right). (left) t er 16 days, T max 5 T () (right) t er 23 days, T max 2 T () Is the piecewise-constant more efficient?
33 ACI: {1,...,3 }. T () = 1 5, E() = I () =. 14 T E 14 T E Single-run of piecewise-constant (left) and impulsive (right). (left) t er 15 days, T max 1.5 T () (right) t er 15 days, T max 1.5 T () Are the therapies equivalent?
34 Combined ACI & IL-2. T () = 1 5, E() = I () =. 45 T E 45 T E Single-run of impulsive therapies: (left) synch., (right) asynch. deliveries. Θ IL i = i +.5 Θ E i = i i = 1,..., 3 (left) t er 31 days, E max (right) t er 41 days, E max Θ IL i = Θ E i = i i = 1,..., 3 Is the IS response more efficient in the synchronous case?
35 ACI. {1,...,3 }. E() = I () = and bigger initial tumour. 3e+6 T E 1.4e+8 T E 2.5e e e e e e e+7 1e+6 4e+7 5 2e Single-run of impulsive therapies: (left) T () = 1 6, (left) T () = 1 7. In (right) with the same injection rate of (left) T (t) 1 9 and t er. (left) t er 25 days, T max (right) t er 25 days, T max Does T ()/rates affect the final outcome of the T-IS interplay?
36 ACI. {1,...,3 }. T () = 1 5, E() = I () = End of the therapy.6 End of the therapy Empirical PDF of piecewise-constant (left) and impulsive (right) ACIs. Variable rates of injection: multiplicative terms w {5, 2.5, 2, 1.5, 1}. (left), (right) t er < 3 for w 2, t er > 3 for w 1.5 Is [1.5; 2] is a range for w to have t er < 3 (as often desired)?
37 Conclusions A hybrid model enforcing the following hypothesis: (i) daily IL-based: piecewise-constant > impulsive; (ii) daily ACI: piecewise-constant = impulsive; (iii) impulsive ACI: weekly > daily 1 (iv) weekly ACI: impulsive > piecewise-constant; (v) combined impulsive ACI+IL-2: synchronous > asynchronous. What is next? realistic extrinsic noise (e.g. stochastic nonlinear equations); systematic exploration of the parameters space; more complex effects of therapies (e.g. delayed); more efficient simulation strategies; individual-based models. 1 metronomic delivery, recognized to be more effective for other anti-tumor therapies such as anti-angiogenesis therapies and chemotherapies.
Formal Modeling of Biological Systems with Delays
Universita degli Studi di Pisa Dipartimento di Informatica Dottorato di Ricerca in Informatica Ph.D. Thesis Proposal Formal Modeling of Biological Systems with Delays Giulio Caravagna caravagn@di.unipi.it
More informationOn the Interpretation of Delays in Delay Stochastic Simulation of Biological Systems
On the Interpretation of Delays in Delay Stochastic Simulation of Biological Systems Roberto Barbuti Giulio Caravagna Andrea Maggiolo-Schettini Paolo Milazzo Dipartimento di Informatica, Università di
More informationOn the Interpretation of Delays in Delay Stochastic Simulation of Biological Systems
On the Interpretation of Delays in Delay Stochastic Simulation of Biological Systems Roberto Barbuti Giulio Caravagna Andrea Maggiolo-Schettini Paolo Milazzo Dipartimento di Informatica, Università di
More informationIntroduction to some topics in Mathematical Oncology
Introduction to some topics in Mathematical Oncology Franco Flandoli, University of Pisa y, Finance and Physics, Berlin 2014 The field received considerable attentions in the past 10 years One of the plenary
More informationNetworks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationarxiv: v1 [q-bio.to] 21 Apr 2018
arxiv:1805.01009v1 [q-bio.to] 21 Apr 2018 Approximate Analytical Solution of a Cancer Immunotherapy Model by the Application of Differential Transform and Adomian Decomposition Methods Abstract Alireza
More informationThursday. Threshold and Sensitivity Analysis
Thursday Threshold and Sensitivity Analysis SIR Model without Demography ds dt di dt dr dt = βsi (2.1) = βsi γi (2.2) = γi (2.3) With initial conditions S(0) > 0, I(0) > 0, and R(0) = 0. This model can
More informationIntroduction to stochastic multiscale modelling of tumour growth
Introduction to stochastic multiscale modelling of tumour growth Tomás Alarcón ICREA & Centre de Recerca Matemàtica T. Alarcón (ICREA & CRM, Barcelona, Spain) Lecture 1 CIMPA, Santiago de Cuba, June 2016
More informationModeling biological systems with delays in Bio-PEPA
Modeling biological systems with delays in Bio-PEPA Giulio Caravagna Dipartimento di Informatica, Università di Pisa, argo Pontecorvo 3, 56127 Pisa, Italy. caravagn@di.unipi.it Jane Hillston aboratory
More informationZool 3200: Cell Biology Exam 5 4/27/15
Name: Trask Zool 3200: Cell Biology Exam 5 4/27/15 Answer each of the following short answer questions in the space provided, giving explanations when asked to do so. Circle the correct answer or answers
More informationMODELLING TUMOUR-IMMUNITY INTERACTIONS WITH DIFFERENT STIMULATION FUNCTIONS
Int. J. Appl. Math. Comput. Sci., 2003, Vol. 13, No. 3, 307 315 MODELLING TUMOUR-IMMUNITY INTERACTIONS WITH DIFFERENT STIMULATION FUNCTIONS PETAR ZHIVKOV, JACEK WANIEWSKI Institute of Applied Mathematics
More informationStochastic Chemical Kinetics
Stochastic Chemical Kinetics Joseph K Scott November 10, 2011 1 Introduction to Stochastic Chemical Kinetics Consider the reaction I + I D The conventional kinetic model for the concentration of I in a
More informationGene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back
Gene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back June 19, 2007 Motivation & Basics A Stochastic Approach to Gene Expression Application to Experimental
More informationMA 777: Topics in Mathematical Biology
MA 777: Topics in Mathematical Biology David Murrugarra Department of Mathematics, University of Kentucky http://www.math.uky.edu/~dmu228/ma777/ Spring 2018 David Murrugarra (University of Kentucky) Lecture
More informationIntroduction to Population Dynamics
MATH235: Differential Equations with Honors Introduction to Population Dynamics Quote of Population Dynamics Homework Wish I knew what you were looking for. Might have known what you would find. Under
More informationStochastic modelling in Mathematical Biology
Stochastic modelling in Mathematical Biology Daniel Sánchez-Taltavull Centre de Recerca Matemàtica March 4th 2013 D. Sánchez-Taltavull (CRM) Stochastic modelling in Mathematical Biology March 4th 2013
More informationAsynchronous oscillations due to antigenic variation in Malaria Pf
Asynchronous oscillations due to antigenic variation in Malaria Pf Jonathan L. Mitchell and Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX SIAM LS, Pittsburgh, 21 Outline
More information2. The development of revolutionized the of life.
Science 10 Unit 7 Worksheet Chapter 15, Part 1. 1. Briefly describe the three main parts of cell theory: 2. The development of revolutionized the of life. 3. Individual cells need to take in to build and
More informationEfficient Leaping Methods for Stochastic Chemical Systems
Efficient Leaping Methods for Stochastic Chemical Systems Ioana Cipcigan Muruhan Rathinam November 18, 28 Abstract. Well stirred chemical reaction systems which involve small numbers of molecules for some
More informationLecture 4 The stochastic ingredient
Lecture 4 The stochastic ingredient 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 informationMathematical Modeling of Immune Responses to Hepatitis C Virus Infection
East Tennessee State University Digital Commons @ East Tennessee State University Electronic Theses and Dissertations 12-2014 Mathematical Modeling of Immune Responses to Hepatitis C Virus Infection Ivan
More informationStochastic Simulation of Biochemical Reactions
1 / 75 Stochastic Simulation of Biochemical Reactions Jorge Júlvez University of Zaragoza 2 / 75 Outline 1 Biochemical Kinetics 2 Reaction Rate Equation 3 Chemical Master Equation 4 Stochastic Simulation
More informationAnalysis and Simulation of Biological Systems
Analysis and Simulation of Biological Systems Dr. Carlo Cosentino School of Computer and Biomedical Engineering Department of Experimental and Clinical Medicine Università degli Studi Magna Graecia Catanzaro,
More informationHybrid Modeling of Tumor Induced Angiogenesis
University of Iceland Reykjavik, Iceland July 9, 25 Hybrid Modeling of Tumor Induced Angiogenesis L. L. Bonilla, M. Alvaro, V. Capasso, M. Carretero, F. Terragni Gregorio Millán Institute for Fluid Dynamics,
More informationLecture 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 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 informationarxiv: v2 [q-bio.qm] 12 Jan 2017
Approximation and inference methods for stochastic biochemical kinetics - a tutorial review arxiv:1608.06582v2 [q-bio.qm] 12 Jan 2017 David Schnoerr 1,2,3, Guido Sanguinetti 2,3, and Ramon Grima 1,3,*
More informationModelling pulsed immunotherapy of tumour-immune interaction
Accepted Manuscript Modelling pulsed immunotherapy of tumour-immune interaction Jin Yang, Sanyi Tang, Robert A. Cheke PII: S378-47544)34- DOI: http://dx.doi.org/.6/j.matcom.4.9. Reference: MATCOM 4 To
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 informationMODELING THE ABSORPTION, CIRCULATION, AND METABOLISM OF TIRAPAZAMINE
MODELING THE ABSORPTION, CIRCULATION, AND METABOLISM OF TIRAPAZAMINE Olivia Nguyen and Xiaoyu Shi Problem Statement: Tirapazamine (TPZ; 1,2,4-benzotriazin-3-amine 1,4-dioxide) is a hypoxia-activated prodrug
More informationProblem Set 5. 1 Waiting times for chemical reactions (8 points)
Problem Set 5 1 Waiting times for chemical reactions (8 points) In the previous assignment, we saw that for a chemical reaction occurring at rate r, the distribution of waiting times τ between reaction
More informationChemotherapy in heterogeneous cultures of cancer cells with interconversion
Chemotherapy in heterogeneous cultures of cancer cells with interconversion Rui DILãO Institut des Hautes Études Scientifiques 35, route de Chartres 91440 Bures-sur-Yvette (France) Juin 2013 IHES/M/13/20
More informationIntroduction. Dagmar Iber Jörg Stelling. CSB Deterministic, SS 2015, 1.
Introduction Dagmar Iber Jörg Stelling joerg.stelling@bsse.ethz.ch CSB Deterministic, SS 2015, 1 Origins of Systems Biology On this assumption of the passage of blood, made as a basis for argument, and
More informationDynamic-equilibrium solutions of ordinary differential equations and their role in applied problems
Applied Mathematics Letters 21 (2008) 320 325 www.elsevier.com/locate/aml Dynamic-equilibrium solutions of ordinary differential equations and their role in applied problems E. Mamontov Department of Physics,
More informationMITOCW ocw-5_60-s08-lec17_300k
MITOCW ocw-5_60-s08-lec17_300k The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free.
More informationSTOCHASTIC CHEMICAL KINETICS
STOCHASTIC CHEICAL KINETICS Dan Gillespie GillespieDT@mailaps.org Current Support: Caltech (NIGS) Caltech (NIH) University of California at Santa Barbara (NIH) Past Support: Caltech (DARPA/AFOSR, Beckman/BNC))
More informationFormal Methods and Systems Biology: The Calculus of Looping Sequences
Formal Methods and Systems Biology: The Calculus of Looping Sequences Paolo Milazzo Dipartimento di Informatica, Università di Pisa, Italy Verona January 22, 2008 Paolo Milazzo (Università di Pisa) Formal
More informationStochastic Simulation.
Stochastic Simulation. (and Gillespie s algorithm) Alberto Policriti Dipartimento di Matematica e Informatica Istituto di Genomica Applicata A. Policriti Stochastic Simulation 1/20 Quote of the day D.T.
More informationControlling Systemic Inflammation Using NMPC. Using Nonlinear Model Predictive Control with State Estimation
Controlling Systemic Inflammation Using Nonlinear Model Predictive Control with State Estimation Gregory Zitelli, Judy Day July 2013 With generous support from the NSF, Award 1122462 Motivation We re going
More informationReaction time distributions in chemical kinetics: Oscillations and other weird behaviors
Introduction The algorithm Results Summary Reaction time distributions in chemical kinetics: Oscillations and other weird behaviors Ramon Xulvi-Brunet Escuela Politécnica Nacional Outline Introduction
More informationStochastic Modeling of Chemical Reactions
Stochastic Modeling of Chemical Reactions Eric Haseltine Department of Chemical Engineering University of Wisconsin-Madison TWMCC 27 September Texas-Wisconsin Modeling and Control Consortium 22TWMCC Outline
More informationModelling Biochemical Pathways with Stochastic Process Algebra
Modelling Biochemical Pathways with Stochastic Process Algebra Jane Hillston. LFCS, University of Edinburgh 13th April 2007 The PEPA project The PEPA project started in Edinburgh in 1991. The PEPA project
More informationSPA for quantitative analysis: Lecture 6 Modelling Biological Processes
1/ 223 SPA for quantitative analysis: Lecture 6 Modelling Biological Processes Jane Hillston LFCS, School of Informatics The University of Edinburgh Scotland 7th March 2013 Outline 2/ 223 1 Introduction
More informationModelling in Systems Biology
Modelling in Systems Biology Maria Grazia Vigliotti thanks to my students Anton Stefanek, Ahmed Guecioueur Imperial College Formal representation of chemical reactions precise qualitative and quantitative
More informationMathematics, Box F, Brown University, Providence RI 02912, Web site:
April 24, 2012 Jerome L, Stein 1 In their article Dynamics of cancer recurrence, J. Foo and K. Leder (F-L, 2012), were concerned with the timing of cancer recurrence. The cancer cell population consists
More informationSupplementary Figures
Supplementary Figures a x 1 2 1.5 1 a 1 = 1.5 a 1 = 1.0 0.5 a 1 = 1.0 b 0 0 20 40 60 80 100 120 140 160 2 t 1.5 x 2 1 0.5 a 1 = 1.0 a 1 = 1.5 a 1 = 1.0 0 0 20 40 60 80 100 120 140 160 t Supplementary Figure
More informationBioimage Informatics. Lecture 23, Spring Emerging Applications: Molecular Imaging
Bioimage Informatics Lecture 23, Spring 2012 Emerging Applications: Molecular Imaging Lecture 23 April 25, 2012 1 Outline Overview of molecular imaging Molecular imaging modalities Molecular imaging applications
More informationAutomatic Differentiation and Neural Networks
Statistical Machine Learning Notes 7 Automatic Differentiation and Neural Networks Instructor: Justin Domke 1 Introduction The name neural network is sometimes used to refer to many things (e.g. Hopfield
More informationBiology I Level - 2nd Semester Final Review
Biology I Level - 2nd Semester Final Review The 2 nd Semester Final encompasses all material that was discussed during second semester. It s important that you review ALL notes and worksheets from the
More informationStochastic Simulation Methods for Solving Systems with Multi-State Species
Stochastic Simulation Methods for Solving Systems with Multi-State Species Zhen Liu Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of
More informationApproximation of epidemic models by diffusion processes and their statistical inferencedes
Approximation of epidemic models by diffusion processes and their statistical inferencedes Catherine Larédo 1,2 1 UR 341, MaIAGE, INRA, Jouy-en-Josas 2 UMR 7599, LPMA, Université Paris Diderot December
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 informationCybergenetics: Control theory for living cells
Department of Biosystems Science and Engineering, ETH-Zürich Cybergenetics: Control theory for living cells Corentin Briat Joint work with Ankit Gupta and Mustafa Khammash Introduction Overview Cybergenetics:
More informationPrinciples of Synthetic Biology: Midterm Exam
Principles of Synthetic Biology: Midterm Exam October 28, 2010 1 Conceptual Simple Circuits 1.1 Consider the plots in figure 1. Identify all critical points with an x. Put a circle around the x for each
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 information2. 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 informationLecture 7: Simple genetic circuits I
Lecture 7: Simple genetic circuits I Paul C Bressloff (Fall 2018) 7.1 Transcription and translation In Fig. 20 we show the two main stages in the expression of a single gene according to the central dogma.
More informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
More informationWriting Assignment 1: Comparing Harvesting Models
Writing Assignment : Comparing Harvesting Models Paul Allen Fall 28 Abstract We analyze two models for the harvesting of a yogurt culture growing according to the logistic model: constant-rate harvesting
More informationModeling with differential equations
Mathematical Modeling Lia Vas Modeling with differential equations When trying to predict the future value, one follows the following basic idea. Future value = present value + change. From this idea,
More informationCourse plan Academic Year Qualification MSc on Bioinformatics for Health Sciences. Subject name: Computational Systems Biology Code: 30180
Course plan 201-201 Academic Year Qualification MSc on Bioinformatics for Health Sciences 1. Description of the subject Subject name: Code: 30180 Total credits: 5 Workload: 125 hours Year: 1st Term: 3
More informationIntroduction to Cellular Communication *
OpenStax-CNX module: m53235 1 Introduction to Cellular Communication * Steven Telleen This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 4.0 1 Why Cells Communicate
More informationFinal Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger. Project I: Predator-Prey Equations
Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger Project I: Predator-Prey Equations The Lotka-Volterra Predator-Prey Model is given by: du dv = αu βuv = ρβuv
More informationProgrammed Cell Death
Programmed Cell Death Dewajani Purnomosari Department of Histology and Cell Biology Faculty of Medicine Universitas Gadjah Mada d.purnomosari@ugm.ac.id What is apoptosis? a normal component of the development
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationSignal Transduction. Dr. Chaidir, Apt
Signal Transduction Dr. Chaidir, Apt Background Complex unicellular organisms existed on Earth for approximately 2.5 billion years before the first multicellular organisms appeared.this long period for
More information8 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 informationA. Incorrect! The Cell Cycle contains 4 distinct phases: (1) G 1, (2) S Phase, (3) G 2 and (4) M Phase.
Molecular Cell Biology - Problem Drill 21: Cell Cycle and Cell Death Question No. 1 of 10 1. Which of the following statements about the cell cycle is correct? Question #1 (A) The Cell Cycle contains 3
More informationPersistence and Stationary Distributions of Biochemical Reaction Networks
Persistence and Stationary Distributions of Biochemical Reaction Networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Discrete Models in Systems Biology SAMSI December
More informationFigure 1: Ca2+ wave in a Xenopus oocyte following fertilization. Time goes from top left to bottom right. From Fall et al., 2002.
1 Traveling Fronts Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 2 When mature Xenopus oocytes (frog
More informationOn the local kinetics component of biological systems in competition
On the local kinetics component of biological systems in competition A. Ciancio and B.F.F. Flora Abstract. In various previous papers, mathematical models describing the evolution of a tumor in presence
More informationnatural development from this collection of knowledge: it is more reliable to predict the property
1 Chapter 1 Introduction As the basis of all life phenomena, the interaction of biomolecules has been under the scrutiny of scientists and cataloged meticulously [2]. The recent advent of systems biology
More informationHost-Pathogen Interaction. PN Sharma Department of Plant Pathology CSK HPKV, Palampur
Host-Pathogen Interaction PN Sharma Department of Plant Pathology CSK HPKV, Palampur-176062 PATHOGEN DEFENCE IN PLANTS A BIOLOGICAL AND MOLECULAR VIEW Two types of plant resistance response to potential
More informationGRÜNWALD-LETNIKOV SCHEME FOR SYSTEM OF CHRONIC MYELOGENOUS LEUKEMIA FRACTIONAL DIFFERENTIAL EQUATIONS AND ITS OPTIMAL CONTROL OF DRUG TREATMENT
GRÜNWALD-LETNIKOV SHEME FOR SYSTEM OF HRONI MYELOGENOUS LEUKEMIA FRATIONAL DIFFERENTIAL EQUATIONS AND ITS OPTIMAL ONTROL OF DRUG TREATMENT ESMAIL HESAMEDDINI AND MAHIN AZIZI DEPARTMENT OF MATHEMATIAL SIENES,
More informationBiological Pathways Representation by Petri Nets and extension
Biological Pathways Representation by and extensions December 6, 2006 Biological Pathways Representation by and extension 1 The cell Pathways 2 Definitions 3 4 Biological Pathways Representation by and
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationI. Specialization. II. Autonomous signals
Multicellularity Up to this point in the class we have been discussing individual cells, or, at most, populations of individual cells. But some interesting life forms (for example, humans) consist not
More informationSimulating stochastic epidemics
Simulating stochastic epidemics John M. Drake & Pejman Rohani 1 Introduction This course will use the R language programming environment for computer modeling. The purpose of this exercise is to introduce
More informationNumerical Methods in Cancer Models
Numerical Methods in Cancer Models Doron Levy Department of Mathematics and Center for Scientific Computation and Mathematical Modeling (CSCAMM) University of Maryland, College Park Plan Motivation 1 What
More informationCell Cycle Modeling for Budding Yeast with Stochastic Simulation Algorithms
Cell Cycle Modeling for Budding Yeast with Stochastic Simulation Algorithms Tae-Hyuk Ahn 1, Layne T Watson 1,2, Yang Cao 1, Clifford A Shaffer 1, and William T Baumann 3 1 Department of Computer Science,
More informationAn Introduction to Stochastic Simulation
Stephen Gilmore Laboratory for Foundations of Computer Science School of Informatics University of Edinburgh PASTA workshop, London, 29th June 2006 Background The modelling of chemical reactions using
More informationMEDLINE Clinical Laboratory Sciences Journals
Source Type Publication Name ISSN Peer-Reviewed Academic Journal Acta Biochimica et Biophysica Sinica 1672-9145 Y Academic Journal Acta Physiologica 1748-1708 Y Academic Journal Aging Cell 1474-9718 Y
More informationNonlinear pharmacokinetics
5 Nonlinear pharmacokinetics 5 Introduction 33 5 Capacity-limited metabolism 35 53 Estimation of Michaelis Menten parameters(v max andk m ) 37 55 Time to reach a given fraction of steady state 56 Example:
More informationEstimating transition times for a model of language change in an age-structured population
Estimating transition times for a model of language change in an age-structured population W. Garrett Mitchener MitchenerG@cofc.edu http://mitchenerg.people.cofc.edu Abstract: Human languages are stable
More informationUnder the Radar Screen: How Bugs Trick Our Immune Defenses
Under the Radar Screen: How Bugs Trick Our Immune Defenses Session 2: Phagocytosis Marie-Eve Paquet and Gijsbert Grotenbreg Whitehead Institute for Biomedical Research Salmonella Gram negative bacteria
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 informationOnce Upon A Time, There Was A Certain Ludwig
Once Upon A Time, There Was A Certain Ludwig Statistical Mechanics: Ensembles, Distributions, Entropy and Thermostatting Srinivas Mushnoori Chemical & Biochemical Engineering Rutgers, The State University
More informationUCLA: Math 3B Problem set 6 (solutions) Fall, 2018
This week you will practice writing differential equations modelling real world phenomena as well as understanding population models. You will also get practice solving separable differential equations.
More informationThe Calculus of Looping Sequences
The Calculus of Looping Sequences Roberto Barbuti, Giulio Caravagna, Andrea MaggioloSchettini, Paolo Milazzo, Giovanni Pardini Dipartimento di Informatica, Università di Pisa, Italy Bertinoro June 7, 2008
More informationPredator-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 informationCourse Outline. Key Themes CPIB SUMMER SCHOOL 2010: INTRODUCTION TO BIOLOGICAL MODELLING. University of Nottingham. September 2011
CPIB SUMMER SCHOOL 2010: INTRODUCTION TO BIOLOGICAL MODELLING University of Nottingham 5-8 1 Course Outline Aims: To introduce modelling and quantitative approaches to biology To explain where equations
More informationCELL PRACTICE TEST
Name: Date: 1. As a human red blood cell matures, it loses its nucleus. As a result of this loss, a mature red blood cell lacks the ability to (1) take in material from the blood (2) release hormones to
More informationLongtime behavior of stochastically modeled biochemical reaction networks
Longtime behavior of stochastically modeled biochemical reaction networks David F. Anderson Department of Mathematics University of Wisconsin - Madison ASU Math Biology Seminar February 17th, 217 Overview
More informationRoadmap. Sexual Selection. Evolution of Multi-Gene Families Gene Duplication Divergence Concerted Evolution Survey of Gene Families
1 Roadmap Sexual Selection Evolution of Multi-Gene Families Gene Duplication Divergence Concerted Evolution Survey of Gene Families 2 One minute responses Q: How do aphids start producing males in the
More informationgrowth, the replacement of damaged cells, and development from an embryo into an adult. These cell division occur by means of a process of mitosis
The Cellular Basis of Reproduction and Inheritance PowerPoint Lectures for Campbell Biology: Concepts & Connections, Seventh Edition Reece, Taylor, Simon, and Dickey Biology 1408 Dr. Chris Doumen Lecture
More informationExtending the Tools of Chemical Reaction Engineering to the Molecular Scale
Extending the Tools of Chemical Reaction Engineering to the Molecular Scale Multiple-time-scale order reduction for stochastic kinetics James B. Rawlings Department of Chemical and Biological Engineering
More information1. 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 informationOptimal Control of a Cancer Cell Model with Delay
Math. Model. Nat. Phenom. Vol. 5, No. 3, 21, pp. 63-75 DOI: 1.151/mmnp/21535 Optimal Control of a Cancer Cell Model with Delay C. Collins 1, K.R. Fister 2 and M. Williams 3 1 Department of Mathematics,
More informationBiology EOC Review Study Questions
Biology EOC Review Study Questions Microscopes and Characteristics of Life 1. How do you calculate total magnification on a compound light microscope? 2. What is the basic building block of all living
More informationDynamical-Systems Perspective to Stem-cell Biology: Relevance of oscillatory gene expression dynamics and cell-cell interaction
Dynamical-Systems Perspective to Stem-cell Biology: Relevance of oscillatory gene expression dynamics and cell-cell interaction Kunihiko Kaneko Universal Biology Inst., Center for Complex-Systems Biology,
More information