Inference of genomic network dynamics with non-linear ODEs
|
|
- Janis Simpson
- 5 years ago
- Views:
Transcription
1 Inference of genomic network dynamics with non-linear ODEs University of Groningen ernst 11 June 2015 Joint work with Mahdi Mahmoudi, Itai Dattner and Ivan Vujačić
2 Motivation: dynamic genomic processes The amount of RNA of a particular gene depends on the rate of DNA decay (depends on temperature,...) the amount of transcription factor the rate of transcription (depends on other factors,...)
3 Law of Mass Action Given a basic reaction A + B k 1,k 1 C the rate of forward and backward and backward reactions is linearly proportional with concentration A, B and C respectively: d[a] dt = k 1 [C] k 1 [A][B].
4 Enzymes (typically) large proteins catalysts that change the speed of reaction, but not its extent. highly specific: speed up (typically) a single biochemical reaction. are highly regulated themselves by: phosphorylation calcium ATP their own products, etc.
5 Basic problem of enzyme kinetics Let S be a substrate which together with enzyme E results in a product P: S + E P + E. Problem: according to Law of Mass Action rate increases linearly with [S], but it actually saturates. Michaelis and Menten (1913): an intermediate complex C S + E k 1,k 1 C k 2 P + E.
6 Pseudo steady-state approximation Mass action equilibrium: (k 1 + k 2 )[C] = k 1 [S][E] combined with the finite amount of enzyme m: m = [C] + [E], We get so K[C] = [S](m [C]), with K = (k 1 + k 2 )/k 1 [C](K + [S]) = [S]m. The rate of production of product P v P = k 2 [C] = k 2[S] K + [S]
7 E-Coli SOS system 30 genes controlled by protein called LexA. 1 ODE per gene, particular kinetics parameters.
8 Modelling transcription kinetics in a SIM Mean transcription dynamics follow Michaelis-Menten kinetics: µ i (t) = p i (t, η) δ i µ i (t), where p i (t, η) = RNA production term of gene i by protein η, δ i µ i (t) = natural RNA decay of gene i. If gene i only depends on activity η of one protein: η(t) p i (t, η) = α i + β i γ i + η(t) }{{} Activator or p i (t, η) = α i + β i 1 γ i + η(t) }{{} Repressor where α i = basal transcription rate of gene i, β i = maximal reaction-rate for gene i, γ i = half-saturation constant.
9 Cooperativity Enzyme can bind two substrate molecules at different binding sites: which results in the following rate of production of P v P = (k 2K 2 + k 4 [S])m[S] K 1 K 2 + K 2 [S] + [S] 2.
10 Independent binding sites If there are 2 independent binding sites to which enzyme can bind: which results in doubling the rate of production of P v P = k 2 [C] = 2k 2[S] K + [S]
11 ... or more complex systems: fatty acid β-oxidation Differential equations play crucial role in science
12 Statistical model Consider x (t) = f (x(t), t; θ), t [0, T ], Noisy observations of solution x(t; θ, ξ) are available: y(t i ) = x(t i ; θ, ξ) + ε(t i ), i = 1,..., n, where 0 t 1 t n T and with initial values x j (0) = ξ j. Goal: inverse ODE problem 1 Estimate the vector of parameters θ. 2 When initial values ξ are unknown, they are also estimated.
13 Possible approaches 1 Maximum likelihood (Khanin et al. 2007, Xue et al. 2010): PRO: has guaranteed properties of estimator. CON: requires solution of ODE. 2 Reproducing Kernel Hilbert Space (Gonzalez et al. 2014): CON: has no guaranteed properties of estimator. PRO: does not require solution of ODE. 3 Collocation Method (Ramsay et al. 2007): PRO: consistency. PRO: does not require solution of ODE. CON: Iterative procedure with tuning parameters
14 AIM: INFER DYNAMIC PARAMETERS OF SYSTEM Let θ = (θ, ξ), including initial state. Our proposed framework: iteration of 1. x(θ) = argmin x Xm T α,γ (x θ), 2. θ n = argmin θ Θ M n (θ x(θ), Y). Aim Define T α,γ and M n such that: It yields asymptotically efficient estimator. It can handle partially observed systems.
15 Well-posedness in sense of Hadamard Let F : X Y where X, Y are linear normed spaces and solve: F (x) = y, (1) x X, y Y. The problem (1) is well-posed in the sense of Hadamard on (X, Y) if: 1 The solution of (1) exists. 2 It is unique. 3 It is continuous with respect to y. The problem (1) is ill-posed on (X, Y) if it is not well-posed.
16 Solving ODEs Consider for fixed θ. x (t) = f (x(t), t; θ), t [0, T ], Suppressing dependence on θ, define F (x) = x f (x,, θ), then ODE is equivalent to F (x) = 0 d. In principle, this is an ill-posed problem: lack of uniqueness
17 Quasi-solution for ill-posed problems Equation F (x) = y, can be solved on a set S X by minimizing objective functional J (x) = F (x) y 2. Quasi-solution Minimizer of J on S is called quasi-solution, pseudo solution or least squares solution of F (x) = y. Remark: Idea dates back to beginning of 19th century: Gauss, Legendre.
18 J : ODE objective functional Given F (x) = x f (x,, θ), the corresponding ODE objective functional is J (x) = x f (x,, θ) 2 2,w. where di=1 L 2 norm, x 2,w = w T i 0 x i 2(t)dt. with possible weights w = (w 1,..., w d ). Minimizing J = solving ODE
19 What subset S? Finite-dimensional approximation Consider finite-dimensional subspaces X m X (m = 1, 2,...), such that 1 X 1 X m=1 X m is dense in X. X m C 1 [0, T ] linear subspace with basis {h 1,..., h m }, i.e. x X m m x i (t) = β ik h k (t) = βi h(t), k=1 such as cubic splines with increasing knots. Approximate ODE solution The solution x = argmin x Xm x f (x,, θ) 2,w, is quasi-solution of ODE. Remarks: in statistics literature X m s are called sieves.
20 What subset S? Finite-dimensional approximation Approximate ODE solution The solution x = argmin x Xm x f (x,, θ) 2,w, is quasi-solution of ODE. Remarks: in statistics literature X m s are called sieves.
21 Tikhonov regularization Ω is stabilizing functional, which incorporates a priori information on smoothness of solution x. is usually given by norm or semi-norm on X. Tikhonov regularization involves minimization of Tikhonov functional T α (x) = J (x) + αω(x), where α 0 is regularization parameter For ODE we consider norm in Sobolev space {H 2 [0, T ]} d d T Ω(x) = v i {x i (t)} 2 dt. i=1 0
22 Generalized Tikhonov regularization S is similarity functional, which incorporates a priori information on values of x. Generalized Tikhonov regularization involves minimization of generalized Tikhonov functional T α,γ (x) = J (x) + αω(x) + γs(x), where γ 0 is penalty parameter. We consider (for example) minus log-likelihood: d 1 n S(x) = log p(y i (t j ) x i (t j )). i=1 j=1
23 General outline of our method Consider x β = β t h, such that m x i (t) = β ik h k (t) = βi h(t). k=1 1 For fixed θ, GTR solution x = x β is given by 2 For fixed x, obtain solution β = argmin β R dmt α,γ (x β ). θ n = argmin θ Θ M n (θ x(θ), Y). Then iterate steps 1 and 2.
24 Consistency and asymptotic efficiency Theorem (Consistency) Let Assumptions 1-5 from Qi and Zhao (2010) hold. If as n 1 α n 0 2 γ n 0 then θ n θ 0 = o P (1). Theorem (Asymptotic normality & efficiency) Let Assumptions 1-6 from Qi and Zhao (2010) hold. If α n = o(n 2 ) and γ n = o(n 2 ), then θ n is asymptotically normal and efficient.
25 ODE Systems linear in parameters Nonlinear ODE systems that are linear in parameters θ x (t) = g(x(t))θ, where g : R d R d p. Example 1. Lotka-Volterra system { x 1 (t) = θ 1 x 1 (t) θ 2 x 1 (t)x 2 (t), x 2 (t) = θ 3x 2 (t) + θ 4 x 1 (t)x 2 (t), is linear in parameter θ = (θ 1, θ 2, θ 3, θ 4 ) where ( ) x1 (t) x g(x(t)) = 1 x 2 (t) 0 0, 0 0 x 2 (t) x 1 x 2 (t) i.e. (x 1 (t), x 2 (t)) = g(x(t))θ.
26 Example 2. FitzHugh-Nagumo FitzHugh-Nagumo system { x 1 (t) = c{x 1 (t) x 1 (t) 3 /3 + x 2 (t)}, x 2 (t) = 1 c {x 1(t) a + bx 2 (t)}, is linear in parameter θ = (c, 1/c, a/c, b/c) where g(x(t)) = i.e. (x 1 (t), x 2 (t)) = g(x(t))θ. ( x1 (t) x1 3(t)/3 + x ) 2(t) 0 0 0, 0 x 1 (t) 1 x 2 (t)
27 GTR: Window-smoothing Approach We consider X m = piecewise constant functions on m levels. and consider Tikhonov functional (non-dependent on θ!) T 0, (x) = S(x). Obtain x = argmin x X d m T,0 (x) and estimate θ with objective function M n (θ) = T 0 x (t) f ( x(t), θ) q w(t)dt, 1-step GTR estimator!
28 STEP 1. Window smoother for x The window smoother of x is defined as x n (t) = 1 y(t j ), S(t) t j S(t) t S(t). x time x time
29 STEP 2. Objective function for θ Estimate parameters θ and ξ by minimizing M n (ξ, θ) = T 0 x n(t) ξ t 0 g( x n (s)) ds θ 2 dt
30 STEP 2. Objective function for θ Estimate parameters θ and ξ by minimizing M n (ξ, θ) = where T 0 ( T Id Ân = (ξ, θ) x n(t) ξ Â n Ĝ n (t) = Â n = B n = B n t t 0 ) ( ξ θ 0 T 0 T 0 2 g( x n (s)) ds θ dt ) ( T x n(t)dt Ĝn(t) x n (t)dt g( x n (s))ds, Ĝ n (t)dt, T 0 Ĝ n (t) Ĝn(t)dt. ) ( ) ξ. θ
31 STEP 2. Standard least squares form of estimators Minimizing previous quadratic form yields ( ) 1 ( ( ξ, θ) T T Id = Ân  0 x ) n(t)dt n B T n 0 Ĝn(t) x. n (t)dt The estimators of ξ and θ have the form: ξ n = ( TI d θ n = B 1 n T 0 ) Ân B 1 1 T { n  n 0 Id Ĝn(t) { x n (t) ξ n }dt, Ân B 1 n Ĝ n (t) } x n (t)dt,
32 STEP2.... even better! Explicit forms i Ĝ n (t i ) = g( x n (S m )), m=1 Â n (i) = tĝn(t) t i t i (2i 1) 2 g( x n (S i )), Â n = T Ĝn(T ) 2 I (2i 1)g( x n (S i )). 2 i=1 B n = Ĝn(T ) Ân 1 I g( x n (S i )) t 2Ĝn(t) t i 2 t i=1 i I (3i 2 1)g( x n (S i )) g( x n (S i ))) + 3 I 1 (2i 1)Ri g( x n (S i )), 6 2 i=1 i=1 where R i = I m=i+1 g( x n (S m )).
33 Postscript Computational complexity Computing x n : O(nd). Computing ξ n and θ n : O(p 3 + d 3 + ndp 2 ). Consistency: Window-based estimators of θ n and ξ n are n consistent. A reviewer s comment (Statistics and Computing): It s an unusual idea and at first glance it is surprising that it works.
34 A simulation example FitzHugh-Nagumo system, given by the equations, is linear in the parameter FhNdata consist of x 1 (t) = c{x 1(t) x 1 (t) 3 /3 + x 2 (t)}, x 2 (t) = 1 c {x 1(t) a + bx 2 (t)}, θ = (θ 1, θ 2, θ 3, θ 4 ) = (c, 1/c, a/c, b/c). 41 equally spaced observations on the interval [0,20] obtained from the above model with noise levels σ 1 = σ 2 = 0.5 and parameters a = 0.2, b = 0.2 and c = 3.
35 FhNdata from R package CollocInfer x time x time
36 Window estimator provides a sensible initial guess Window, Window + Generalized Profiling, Generalized Profiling u = (1, 1, 1). Estimation of parameters in Fitz-Hugh Nagumo system Method Ini guess a = 0.2 b = 0.2 c = 3 Rel. error Time(sec) Window NO Ramsay Window u u u u u u 6e+99 6e+99 1e+56 6e
37 Application 1. E-Coli SOS system 30 genes controlled by protein called LexA. 1 ODE per gene, particular kinetics parameters.
38 Application 1. SOS system in Ecoli: Model m genes x 1 (t),..., x m (t) and one transcription factor η(t). Model: ẋ k (t) = β k 1 γ k + η(t) + ϕ k δ k x k (t), for k = 1,..., m, t [0, T ]. Gene-dependent kinetics parameters: β k, γ k, δ k and ϕ k.
39 Application 1. SOS system in Ecoli: Data Data set of 14 expression genes. The genes are targets of the master repressor LexA and their expression is studied under UV exposure (40J/m2 ). Abundance of the mrna molecules measured at 0, 5, 10, 20, 40 and 60 minutes. The LexA abundance is unobserved.
40 Application 1. SOS system in Ecoli: Results Reconstructed the activity of the repressor η(t) gene dini. Reconstructed TF LexA Reconstructed gene dini TFA Gene Expression Time Time Confidence intervals for the parameters CI 95% (β) = (0.56, 1.39) CI 95% (δ) = (0.15, 0.27) CI 95% (γ) = (1.71, 2.83), σ = 0.04.
41 Application 1. But not just that... The aim of statistics is not just to fit data, but also compare scientific hypotheses draw scientific conclusions from the models For example, the following scientific conclusions can be drawn: reca has the fastest effective production rate, r = β/(γ + η), uvrb the slowest. reca and umuc are heavily affected by RNA decay. recn, uvrb and umuc has a non-negligible basal production rate α, so they are possibly also transcribed by another TF.
42 Application 2. Auto-regulated system Streptomyces Streptomyces Coelicolor produces antibiotics. CdaR is a crucial regulator. TF CdaR also activates its own transcription.
43 x Application 2. Generalized MM formulation ẋ = β zm γ + z m + ϕ δx ż = ρx τz A B Stability F 1 High abundant Hysteresis loop mrna F 2 γ Bistability γ+ Low abundant mrna Stability z x β = 6 ρ = 0.5 δ = 0.5 τ = 1 ϕ = γ γ C D Value Stability γ = 2.5 x z Value Bistability γ = 7.5 x 1 z 1 x 2 z Time Time
44 Application 2. CdaR data and inference log(aic) A Abundance B mrna, SCO3217 TF, cdar ^ > Hill coefficient Time order A AIC selects Hill coefficient 2. B Data for autoregulated gene SCO3217 (cdar) and obtained smooth functions after estimating parameters. Parameter estimates lie in bistable region, which may provoke bimodal effects at populational level.
45 Discussion Strange that Statistics has not dealt with ODEs until now! Ramsay et al. (2007) has been an important first step. In this talk, We extended it via Generalized Tikhonov Regularization We simplified it via Window Smoothing approach. Genetics may be quite static, genomics is about dynamics. More good work is needed!
A new statistical framework to infer gene regulatory networks with hidden transcription factors
A new statistical framework to infer gene regulatory networks with hidden transcription factors Ivan Vujačić, Javier González, Ernst Wit University of Groningen Johann Bernoulli Institute PEDS II, June
More informationMarkov Chain Monte Carlo Algorithms for Gaussian Processes
Markov Chain Monte Carlo Algorithms for Gaussian Processes Michalis K. Titsias, Neil Lawrence and Magnus Rattray School of Computer Science University of Manchester June 8 Outline Gaussian Processes Sampling
More informationMath 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA Molecular biology and Biochemical kinetics Molecular biology is one of
More informationCS-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 informationBasic Synthetic Biology circuits
Basic Synthetic Biology circuits Note: these practices were obtained from the Computer Modelling Practicals lecture by Vincent Rouilly and Geoff Baldwin at Imperial College s course of Introduction to
More information56:198:582 Biological Networks Lecture 8
56:198:582 Biological Networks Lecture 8 Course organization Two complementary approaches to modeling and understanding biological networks Constraint-based modeling (Palsson) System-wide Metabolism Steady-state
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 informationBioinformatics 3. V18 Kinetic Motifs. Fri, Jan 8, 2016
Bioinformatics 3 V18 Kinetic Motifs Fri, Jan 8, 2016 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
More informationBioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014"
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
More information7.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 information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationGLOBEX Bioinformatics (Summer 2015) Genetic networks and gene expression data
GLOBEX Bioinformatics (Summer 2015) Genetic networks and gene expression data 1 Gene Networks Definition: A gene network is a set of molecular components, such as genes and proteins, and interactions between
More informationSystem Biology - Deterministic & Stochastic Dynamical Systems
System Biology - Deterministic & Stochastic Dynamical Systems System Biology - Deterministic & Stochastic Dynamical Systems 1 The Cell System Biology - Deterministic & Stochastic Dynamical Systems 2 The
More informationCHAPTER : Prokaryotic Genetics
CHAPTER 13.3 13.5: Prokaryotic Genetics 1. Most bacteria are not pathogenic. Identify several important roles they play in the ecosystem and human culture. 2. How do variations arise in bacteria considering
More informationBioinformatics: Network Analysis
Bioinformatics: Network Analysis Kinetics of Gene Regulation COMP 572 (BIOS 572 / BIOE 564) - Fall 2013 Luay Nakhleh, Rice University 1 The simplest model of gene expression involves only two steps: the
More 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 informationBasic modeling approaches for biological systems. Mahesh Bule
Basic modeling approaches for biological systems Mahesh Bule The hierarchy of life from atoms to living organisms Modeling biological processes often requires accounting for action and feedback involving
More informationthe 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 informationHilbert Uniqueness Method and regularity
Hilbert Uniqueness Method and regularity Sylvain Ervedoza 1 Joint work with Enrique Zuazua 2 1 Institut de Mathématiques de Toulouse & CNRS 2 Basque Center for Applied Mathematics Institut Henri Poincaré
More informationA Simple Protein Synthesis Model
A Simple Protein Synthesis Model James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 3, 213 Outline A Simple Protein Synthesis Model
More informationGeneralized Local Regularization for Ill-Posed Problems
Generalized Local Regularization for Ill-Posed Problems Patricia K. Lamm Department of Mathematics Michigan State University AIP29 July 22, 29 Based on joint work with Cara Brooks, Zhewei Dai, and Xiaoyue
More informationIntrinsic Noise in Nonlinear Gene Regulation Inference
Intrinsic Noise in Nonlinear Gene Regulation Inference Chao Du Department of Statistics, University of Virginia Joint Work with Wing H. Wong, Department of Statistics, Stanford University Transcription
More informationAre DNA Transcription Factor Proteins Maxwellian Demons?
Are DNA Transcription Factor Proteins Maxwellian Demons? * physics as a design constraint on molecular evolution. Alexander Grosberg, Longhua Hu, U. Minnesota, NYU U. Minnesota Biophysical Journal 2008
More informationProgram for the rest of the course
Program for the rest of the course 16.4 Enzyme kinetics 17.4 Metabolic Control Analysis 19.4. Exercise session 5 23.4. Metabolic Control Analysis, cont. 24.4 Recap 27.4 Exercise session 6 etabolic Modelling
More informationSupporting Information. Methods. Equations for four regimes
Supporting Information A Methods All analytical expressions were obtained starting from quation 3, the tqssa approximation of the cycle, the derivation of which is discussed in Appendix C. The full mass
More informationLecture 4: Transcription networks basic concepts
Lecture 4: Transcription networks basic concepts - Activators and repressors - Input functions; Logic input functions; Multidimensional input functions - Dynamics and response time 2.1 Introduction The
More information2 Dilution of Proteins Due to Cell Growth
Problem Set 1 1 Transcription and Translation Consider the following set of reactions describing the process of maing a protein out of a gene: G β g G + M M α m M β m M + X X + S 1+ 1 X S 2+ X S X S 2
More informationHybrid Model of gene regulatory networks, the case of the lac-operon
Hybrid Model of gene regulatory networks, the case of the lac-operon Laurent Tournier and Etienne Farcot LMC-IMAG, 51 rue des Mathématiques, 38041 Grenoble Cedex 9, France Laurent.Tournier@imag.fr, Etienne.Farcot@imag.fr
More informationBistability and a differential equation model of the lac operon
Bistability and a differential equation model of the lac operon Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Fall 2016 M. Macauley
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationInference of Gaussian graphical models and ordinary differential equations Vujacic, Ivan
University of Groningen Inference of Gaussian graphical models and ordinary differential equations Vujacic, Ivan IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if
More information1 Periodic stimulations of the Incoherent Feedforward Loop network
1 Periodic stimulations of the Incoherent Feedforward Loop network In this Additional file, we give more details about the mathematical analysis of the periodic activation of the IFFL network by a train
More informationChapter 7: Metabolic Networks
Chapter 7: Metabolic Networks 7.1 Introduction Prof. Yechiam Yemini (YY) Computer Science epartment Columbia University Introduction Metabolic flux analysis Applications Overview 2 1 Introduction 3 Metabolism:
More informationChapter 15 Active Reading Guide Regulation of Gene Expression
Name: AP Biology Mr. Croft Chapter 15 Active Reading Guide Regulation of Gene Expression The overview for Chapter 15 introduces the idea that while all cells of an organism have all genes in the genome,
More information+ regulation. ribosomes
central dogma + regulation rpl DNA tsx rrna trna mrna ribosomes tsl ribosomal proteins structural proteins transporters enzymes srna regulators RNAp DNAp tsx initiation control by transcription factors
More informationStochastic analysis of biochemical reaction networks with absolute concentration robustness
Stochastic analysis of biochemical reaction networks with absolute concentration robustness David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Boston University
More informationSimulation 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 informationBistability in ODE and Boolean network models
Bistability in ODE and Boolean network models Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson)
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 informationLearning parameters in ODEs
Learning parameters in ODEs Application to biological networks Florence d Alché-Buc Joint work with Minh Quach and Nicolas Brunel IBISC FRE 3190 CNRS, Université d Évry-Val d Essonne, France /14 Florence
More informationMCMC Sampling for Bayesian Inference using L1-type Priors
MÜNSTER MCMC Sampling for Bayesian Inference using L1-type Priors (what I do whenever the ill-posedness of EEG/MEG is just not frustrating enough!) AG Imaging Seminar Felix Lucka 26.06.2012 , MÜNSTER Sampling
More informationCollocation approximation of the monodromy operator of periodic, linear DDEs
p. Collocation approximation of the monodromy operator of periodic, linear DDEs Ed Bueler 1, Victoria Averina 2, and Eric Butcher 3 13 July, 2004 SIAM Annual Meeting 2004, Portland 1=Dept. of Math. Sci.,
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationInteractions between proteins, DNA and RNA. The energy, length and time coordinate system to find your way in the cell
Interactions between proteins, DNA and RNA The energy, length and time coordinate system to find your way in the cell Scanning/atomic force microscope (SFM/AFM) mirror laser beam photo diode fluid cell
More informationWritten 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 informationStochastic and deterministic models of biochemical reaction networks
Stochastic and deterministic models of biochemical reaction networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Case Western Reserve November 9th, 215 Goals of the talk
More informationFrom cell biology to Petri nets. Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE
From cell biology to Petri nets Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE Biology = Concentrations Breitling / 2 The simplest chemical reaction A B irreversible,
More 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 informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More informationL3.1: Circuits: Introduction to Transcription Networks. Cellular Design Principles Prof. Jenna Rickus
L3.1: Circuits: Introduction to Transcription Networks Cellular Design Principles Prof. Jenna Rickus In this lecture Cognitive problem of the Cell Introduce transcription networks Key processing network
More information13.4 Gene Regulation and Expression
13.4 Gene Regulation and Expression Lesson Objectives Describe gene regulation in prokaryotes. Explain how most eukaryotic genes are regulated. Relate gene regulation to development in multicellular organisms.
More informationVilla et al. (2005) Structural dynamics of the lac repressor-dna complex revealed by a multiscale simulation. PNAS 102:
Villa et al. (2005) Structural dynamics of the lac repressor-dna complex revealed by a multiscale simulation. PNAS 102: 6783-6788. Background: The lac operon is a cluster of genes in the E. coli genome
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 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 informationTrajectory Smoothing as a Linear Optimal Control Problem
Trajectory Smoothing as a Linear Optimal Control Problem Biswadip Dey & P. S. Krishnaprasad Allerton Conference - 212 Monticello, IL, USA. October 4, 212 Background and Motivation 1 Background and Motivation
More informationRegulation of metabolism
Regulation of metabolism So far in this course we have assumed that the metabolic system is in steady state For the rest of the course, we will abandon this assumption, and look at techniques for analyzing
More informationBioinformatics 2 - Lecture 4
Bioinformatics 2 - Lecture 4 Guido Sanguinetti School of Informatics University of Edinburgh February 14, 2011 Sequences Many data types are ordered, i.e. you can naturally say what is before and what
More informationBioControl - Week 2, Lecture 1
BioControl - Week 2, Lecture Goals of this lecture: CDS tools for bio-molecular systems Study of equilibria Dynamic performance Periodic behaviors Main eample: Negative autoregulation Suggested readings
More informationSystems Biology Across Scales: A Personal View XIV. Intra-cellular systems IV: Signal-transduction and networks. Sitabhra Sinha IMSc Chennai
Systems Biology Across Scales: A Personal View XIV. Intra-cellular systems IV: Signal-transduction and networks Sitabhra Sinha IMSc Chennai Intra-cellular biochemical networks Metabolic networks Nodes:
More informationSparse Approximation and Variable Selection
Sparse Approximation and Variable Selection Lorenzo Rosasco 9.520 Class 07 February 26, 2007 About this class Goal To introduce the problem of variable selection, discuss its connection to sparse approximation
More informationNonparametric Inference In Functional Data
Nonparametric Inference In Functional Data Zuofeng Shang Purdue University Joint work with Guang Cheng from Purdue Univ. An Example Consider the functional linear model: Y = α + where 1 0 X(t)β(t)dt +
More informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
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 informationModel Selection and Geometry
Model Selection and Geometry Pascal Massart Université Paris-Sud, Orsay Leipzig, February Purpose of the talk! Concentration of measure plays a fundamental role in the theory of model selection! Model
More informationFUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models
FUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models Markus W. Covert Stanford University 0 CRC Press Taylor & Francis Group Boca Raton London New York Contents /... Preface, xi
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 informationMultistability in the lactose utilization network of E. coli. Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J.
Multistability in the lactose utilization network of E. coli Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J. Ruby Abrams Motivation Understanding biological switches in the
More informationNetwork motifs in the transcriptional regulation network (of Escherichia coli):
Network motifs in the transcriptional regulation network (of Escherichia coli): Janne.Ravantti@Helsinki.Fi (disclaimer: IANASB) Contents: Transcription Networks (aka. The Very Boring Biology Part ) Network
More informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
More informationBME 5742 Biosystems Modeling and Control
BME 5742 Biosystems Modeling and Control Lecture 24 Unregulated Gene Expression Model Dr. Zvi Roth (FAU) 1 The genetic material inside a cell, encoded in its DNA, governs the response of a cell to various
More informationBiomolecular Feedback Systems
Biomolecular Feedback Systems Domitilla Del Vecchio MIT Richard M. Murray Caltech Version 1.0b, September 14, 2014 c 2014 by Princeton University Press All rights reserved. This is the electronic edition
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationConstrained optimization. Unconstrained optimization. One-dimensional. Multi-dimensional. Newton with equality constraints. Active-set method.
Optimization Unconstrained optimization One-dimensional Multi-dimensional Newton s method Basic Newton Gauss- Newton Quasi- Newton Descent methods Gradient descent Conjugate gradient Constrained optimization
More informationUnbiased Risk Estimation as Parameter Choice Rule for Filter-based Regularization Methods
Unbiased Risk Estimation as Parameter Choice Rule for Filter-based Regularization Methods Frank Werner 1 Statistical Inverse Problems in Biophysics Group Max Planck Institute for Biophysical Chemistry,
More informationDynamics of Quantum Many-Body Systems from Monte Carlo simulations
s of I.N.F.N. and Physics Department, University of Trento - ITALY INT - Seattle - October 24, 2012 Collaborators s of Alessandro Roggero (Trento) Giuseppina Orlandini (Trento) Outline s of : transform
More informationOscillations in Michaelis-Menten Systems
Oscillations in Michaelis-Menten Systems Hwai-Ray Tung Brown University July 18, 2017 Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, 2017 1 / 22 Background A simple
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationComputational methods for continuous time Markov chains with applications to biological processes
Computational methods for continuous time Markov chains with applications to biological processes David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Penn.
More informationMultistability in the lactose utilization network of Escherichia coli
Multistability in the lactose utilization network of Escherichia coli Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J. Ruby Abrams Agenda Motivation Intro to multistability Purpose
More informationState Space and Hidden Markov Models
State Space and Hidden Markov Models Kunsch H.R. State Space and Hidden Markov Models. ETH- Zurich Zurich; Aliaksandr Hubin Oslo 2014 Contents 1. Introduction 2. Markov Chains 3. Hidden Markov and State
More informationRegularity and approximations of generalized equations; applications in optimal control
SWM ORCOS Operations Research and Control Systems Regularity and approximations of generalized equations; applications in optimal control Vladimir M. Veliov (Based on joint works with A. Dontchev, M. Krastanov,
More informationA First Course on Kinetics and Reaction Engineering. Class 9 on Unit 9
A First Course on Kinetics and Reaction Engineering Class 9 on Unit 9 Part I - Chemical Reactions Part II - Chemical Reaction Kinetics Where We re Going A. Rate Expressions - 4. Reaction Rates and Temperature
More informationNonlinear representation, backward SDEs, and application to the Principal-Agent problem
Nonlinear representation, backward SDEs, and application to the Principal-Agent problem Ecole Polytechnique, France April 4, 218 Outline The Principal-Agent problem Formulation 1 The Principal-Agent problem
More informationLecture 15 (10/20/17) Lecture 15 (10/20/17)
Reading: Ch6; 98-203 Ch6; Box 6- Lecture 5 (0/20/7) Problems: Ch6 (text); 8, 9, 0,, 2, 3, 4, 5, 6 Ch6 (study guide-facts); 6, 7, 8, 9, 20, 2 8, 0, 2 Ch6 (study guide-applying); NEXT Reading: Ch6; 207-20
More informationRealization theory for systems biology
Realization theory for systems biology Mihály Petreczky CNRS Ecole Central Lille, France February 3, 2016 Outline Control theory and its role in biology. Realization problem: an abstract formulation. Realization
More informationGene Regulation and Expression
THINK ABOUT IT Think of a library filled with how-to books. Would you ever need to use all of those books at the same time? Of course not. Now picture a tiny bacterium that contains more than 4000 genes.
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationIntroduction to Bayesian methods in inverse problems
Introduction to Bayesian methods in inverse problems Ville Kolehmainen 1 1 Department of Applied Physics, University of Eastern Finland, Kuopio, Finland March 4 2013 Manchester, UK. Contents Introduction
More informationLecture 27. Transition States and Enzyme Catalysis
Lecture 27 Transition States and Enzyme Catalysis Reading for Today: Chapter 15 sections B and C Chapter 16 next two lectures 4/8/16 1 Pop Question 9 Binding data for your thesis protein (YTP), binding
More informationSimple kinetics of enzyme action
Simple kinetics of enzyme action It is established that enzymes form a bound complex to their reactants (i.e. substrates) during the course of their catalysis and prior to the release of products. This
More informationIntroduction to Mathematical Physiology I - Biochemical Reactions
Introduction to Mathematical Physiology I - Biochemical Reactions J. P. Keener Mathematics Department Math Physiology p.1/28 Introduction The Dilemma of Modern Biology The amount of data being collected
More informationStochastic Transcription Elongation via Rule Based Modelling
Stochastic Transcription Elongation via Rule Based Modelling Masahiro HAMANO hamano@jaist.ac.jp SASB, Saint Malo 2015 1 Purpose of This Talk 2 mechano-chemical TE as Rule Based Modelling 3 Table of Contents
More informationEE C128 / ME C134 Final Exam Fall 2014
EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket
More informationCommunication constraints and latency in Networked Control Systems
Communication constraints and latency in Networked Control Systems João P. Hespanha Center for Control Engineering and Computation University of California Santa Barbara In collaboration with Antonio Ortega
More informationFuncICA for time series pattern discovery
FuncICA for time series pattern discovery Nishant Mehta and Alexander Gray Georgia Institute of Technology The problem Given a set of inherently continuous time series (e.g. EEG) Find a set of patterns
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 informationHomogenization for chaotic dynamical systems
Homogenization for chaotic dynamical systems David Kelly Ian Melbourne Department of Mathematics / Renci UNC Chapel Hill Mathematics Institute University of Warwick November 3, 2013 Duke/UNC Probability
More informationBiostatistics-Lecture 16 Model Selection. Ruibin Xi Peking University School of Mathematical Sciences
Biostatistics-Lecture 16 Model Selection Ruibin Xi Peking University School of Mathematical Sciences Motivating example1 Interested in factors related to the life expectancy (50 US states,1969-71 ) Per
More informationENZYME KINETICS. Medical Biochemistry, Lecture 24
ENZYME KINETICS Medical Biochemistry, Lecture 24 Lecture 24, Outline Michaelis-Menten kinetics Interpretations and uses of the Michaelis- Menten equation Enzyme inhibitors: types and kinetics Enzyme Kinetics
More information