CS-E5880 Modeling biological networks Parameter estimation for biological networks
|
|
- Janis McKinney
- 6 years ago
- Views:
Transcription
1 CS-E5880 Modeling biological networks Parameter estimation for biological networks Jukka Intosalmi Department of Computer Science Aalto University January 16, 2018
2 Outline ODE model calibration as an optimization problem Local gradient based optimization Identifiability and reparameterization ODE model calibration as a statistical estimation problem Practical considerations
3 Parameter estimation Parameter estimation is an important part of ODE modeling Estimation of the reaction rates and initial values Testing if a hypothetical model can produce observed dynamics Model construction in a data-driven manner Experimental data is typically noisy
4 Parameter estimation as an optimization problem (1) Let us consider a one dimensional initial value problem dx dt = f (x, θ), x(0) = x 0, where θ = (θ 1,..., θ d ) is a parameter vector. The system has the solution x(t, θ). The solution can be obtained numerically (recall Euler method) The observed values are denoted by y(t i ), i = 1,..., n where t i are the measurement times.
5 Parameter estimation as an optimization problem (2) Our goal is to minimize the distance between the solution of the ODE model and data In other words, we wish to minimize the objective function L(θ) = n (y(t i ) x(t i, θ)) 2. i=1 For ODE models, closed-form solutions are only rarely available and numerical optimization techniques need to be used. Typically a non-convex optimization problem
6 Optimization techniques Global and local optimization techniques Figure: Illustration of optimization problem (from Fröchlich et al., 2017) In this lecture, we will concentrate on deterministic local optimization Extentions to global optimization through multistart approach
7 Gradient based optimization Search for the nearest minimum of the objective function following the direction of the gradient
8 Gradient descent method 1. Set k = 0 and define initial values for the model parameters θ k 2. Evaluate the objective function L(θ k ) 3. Compute the gradient of the objective function at θ k, L(θ k ) 4. Line search: solve for the r = ˆr minimizing L(θ k r L(θ k )) 5. Set θ k+1 = θ k ˆr L(θ k ) 6. Compare the current objective function value with the previous value; if L(θ k+1 ) > L(θ k ) or if θk+1 θ k < ɛ then the solution is θ k 7. Otherwise, set k := k + 1 and go to step 2
9 Differentiation of the objective function To obtain the gradient, we need to differentiate the objective function L(θ) θ j = n i=1 2(y(t i ) x(t i, θ)) x(t i, θ) θ j. The gradient depends on x(t i,θ) θ j, i = 1,..., n; j = 1,..., d
10 Finite differences approximation for x(t i,θ) θ j The approximation is computed by perturbing the parameters x(t i, θ) x(t i, θ) x(t i, θ + he j ) θ j h where e j is the jth unit vector and h is sufficiently small constant. Requires several numerical solutions for the ODE system No generic way of choosing the constant h In general finite differences approach may result in a poor approximation Better results can be obtained using the sensitivity equations
11 Sensitivity equations We can take a total derivative of the ODE system w.r.t. θ j d dθ j dx dt = d dθ j f (x, θ). By evaluating the total derivative and rearranging, we obtain d dt dx dθ j = f (x, θ) x dx + dθ j f (x, θ) θ j. The sensitivity equation for the parameter θ j can be written in the form ds j f (x, θ) f (x, θ) = s j +, dt x θ j where s j (t) = x(t,θ) θ j are called sensitivities.
12 Solving the sensitivities Analytical formulation for s j (t) = x(t,θ) θ j The sensitivities and the original ODE model form a coupled system and they are solved simultaneously using numerical solvers. The accuracy of the computed sensitivities depends only on the numerical error of the ODE solver!
13 Comparing the finite differences and sensitivity equation approaches Figure: From Raue et al., 2013
14 Example of sensitivity equation construction (1) Let us consider the initial value problem dx dt = α + θ x 1 x + K θ 2x, x(0) = x 0 }{{} f (x,θ) where θ 1 and θ 2 are unknown rate parameters. We have f (x, θ) K = θ 1 x (x + K) 2 θ 2 f (x, θ) = x θ 1 x + K f (x, θ) = x. θ 2
15 Example of sensitivity equation construction (2) The sensitivity equations together with the actual ODE model form the following ODE system dx dt = α + θ x 1 x + K θ 2x ( ) ds 1 dt = K θ 1 (x + K) 2 θ 2 s 1 + x x + K ( ) ds 2 dt = K θ 1 (x + K) 2 θ 2 s 2 x with x(0) = x 0 and s 1 (0) = s 2 (0) = 0.
16 Example of sensitivity equation construction (3) And the system can be solved x s s Time (t) Time (t) Time (t)
17 Sensitivity equation in the general form Let us consider an m-dimensional (x R m ) initial value problem dx dt = f(x, θ), x(0) = x 0, where θ = (θ 1,..., θ d ) is a parameter vector. The sensitivity equations can be expressed in the general form d dx dt dθ = f(x, θ) x dx dθ f(x, θ) +. θ
18 Practical considerations In practice, it is beneficial to carry out the parameter estimation using log-transformed parameters Automated construction of the sensitivity equation Selection of the numerical solver (e.g. Sundials CVODES) Least squares optimization has been used in many applications
19 Global optimization through multistart approach A non-convex optimization problem Multiple local minima Parameter space explored using multiple optimization runs Parallelizable for high-dimensional problems Figure: Illustration of multistart approach (from Raue et al, 2013)
20 Choosing the initial points for multistart (1) Latin hypercube sampling (LHS) Figure: From Raue et al, 2013
21 Choosing the initial points for multistart (2) Figure: From Raue et al, 2013
22 Statistical formulation of the parameter estimation problem Formulating the existing information about the system in statistical terms If we assume normally distributed measurement errors, we can write y(t i ) = x(t i, θ) + ɛ(t i ), where ɛ(t i ) N(0, σ 2 i ). This is equivalent to y(t i ) N(x(t i, θ), σ 2 i ). By assuming independent measurements, we can express the likelihood of the data given the parameters in the form l(θ) = n 1 exp ( (y(t i) x(t i, θ)) 2 ) 2πσi 2 2σi 2. i=1
23 Maximum likelihood estimation The maximum likelihood principle is to choose the value θ = ˆθ which maximizes l(θ). In other words, the maximum likelihood estimate is It is equivalent to solve ˆθ = arg max l(θ) θ ˆθ = arg min L(θ), θ where L(θ) = 2 log(l(θ)) is the negative log-likelihood.
24 Negative log-likelihood (1) n L(θ) = 2 log i=1 1 2πσi 2 exp ( (y(t i) x(t i, θ)) 2 ) 2σ 2 i ( n = log(2πσ2 i ) 1 ( (y(ti ) x(t i, θ)) 2 ) ) 2 σ 2 i=1 i n = (log(2π) + 2 log(σ i ) + (y(t i) x(t i, θ)) 2 ) = i=1 n (2 log(σ i ) + (y(t i) x(t i, θ)) 2 ) + C 1 i=1 σ 2 i σ 2 i
25 Negative log-likelihood (2) If σ i = σ, i = 1,..., n, we have L(θ) = n (2 log(σ) + (y(t i) x(t i, θ)) 2 ) + C 1 i=1 σ 2 = 1 n σ 2 (y(t i ) x(t i, θ)) 2 + C 2 i=1
26 Parameter identifiability In many cases, the parameters are not well determined Infinitely many parameter settings may have equal likelihood The effect of changing one parameter can be compensated by tuning other parameters Figure: From Raue et al, 2010
27 Profile likelihood (1) Parameter identifiability can be studied e.g. using the profile likelihood method. The profile likelihood is defined by PL j (p) = where LL is the log-likelihood. max LL(θ), θ {θ θ j =p} Structural and practical non-identifiability
28 Profile likelihood (2) Figure: From Raue et al, 2010
29 Reparameterization example (1) Modeling metabolic switching in yeast cells. GLUCOSE STATE ethanol production ETHANOL STATE QUIESCENT STATE low glucose level induced transition m g m e m q Figure: Schematic illustration of the model. (Intosalmi et al., unpublished)
30 Reparameterization example (2) dm g dt dm e dt dm q = β 2 m g + β 3 m e dt dg dt = µ 1 m g g γ 1 de dt = µ 1 m g g µ 2 m e e γ 2 γ 3 = µ 1 m g g β 1 1 g + K mg β 2 m g = µ 2 m e e + β 1 1 g + K mg β 3 m e The model output to be linked with data is m g + m e + m q = m.
31 Reparameterization example (3) The initial level of glucose is non-negative Glucose level is a latent variable Normalization w.r.t. the initial glucose level g(0) Formally, Further, we can reparameterize where x 1 = m g, x 4 = x 4 (0) = 1. 1 dg g(0) dt = µ 1 m g g. g(0)γ 1 dx 4 dt = θ 7x 1 x 4, g g(0), θ 7 = µ 1 γ 1, and the initial condition is
32 Reparameterization example (4) In the case of ethanol, we set γ 3 = 1. This results in the expression which has the dimensionless form de dt = µ 1 γ 2 m g g µ 2 m e e dx 5 dt = θ 8x 1 x 4 θ 5 x 2 x 5, where x 5 = e, θ 5 = µ 2, and θ 8 = µ 1 γ 2.
33 Reparameterization example (5) The full system can then be written in the form dx 1 dt = θ 1 1x 1 x 4 θ 2 x 1 θ 4 x 1 x 4 + θ 3 dx 2 dt = θ 1 5x 2 x 5 + θ 2 x 1 θ 6 x 2 x 4 + θ 3 dx 3 dt = θ 4x 1 + θ 6 x 2 dx 4 dt = θ 7x 1 x 4 dx 5 dt = θ 8x 1 x 4 θ 5 x 2 x 5, where (x 1, x 2, x 3, x 4, x 5 ) = (m g, m e, m q, g, e), θ 2 = β 1, θ 3 = K, θ 4 = β 2, and θ 6 = β 3.
34 Reparameterization example (6) Figure: Fitted model (Intosalmi et al., unpublished)
35 T helper cell differentiation example Naive T cells sense the environment and react to different cytokine milieu and differentiate to functionally different T helper subsets Proper regulation and balance between various lymphocyte subsets is central the human immune system For example, so-called T helper 1 cells are involved in defense against intracellular bacteria Molecular mechanisms which control T helper cell differentiation are not at all fully understood A more thorough understanding would allow e.g. better drug design to modulate immune response and help in autoimmune diseases
36 T helper cell differentiation example (2) Different T cell subsets Figure: From (Brusselle et al, 2011)
37 T helper cell differentiation example (3) Immunity Key factors involved Sequential in and Th1 Cell driving Induction T helper by IFN-g 1 differentiation and IL-12 include IFN-γ, T-bet, IL-12 and IL-12Rβ2 The standard/assumed molecular model, until recently, has been the following A B CFigure: From (Schulz et al, 2009) D
38 mrna and protein concentration, respectively. Since naïve cells express T helper cell differentiation example (4) neither IL-12R 2 nor IFN-, only for T-bet we assumed a basal transcript considered had the following form, only differing in the functions f gene de The actualtranscription molecularrate. mechanisms are assumed to follow the following ODE system dtbet dt mrna f 1 Tbet IFN Pr ot,re c Pr ot, Ag Tbet Tbet mrna dtbet dt Pr ot Tbet mrna Tbet Tbet Pr ot dre c dt mrna f Re c Tbet Pr ot, Ag Re c Re c mrna dre c dt Pr ot Re c mrna Re c Re c Pr ot difn dt mrna f IFN Tbet Pr ot,re c Pr ot, Ag IFN IFN m RNA difn dt Pr ot IFN m RNA IFN IFN Pr ot The Figure: rates of From transcription (Schulz are ettaken al, 2009) as saturating, hyperbolic functions Independent activation by multiple signals is described by additive te
39 Sequential Th1 Cell Induction by IFN-g and IL-12 T helper cell differentiation example (5) Figure 3. The Two-Loop Mode (A) Network structure of the literature (B) Schematic representation of t described in Figures 1 and 2: IL-1 T-bet expression and antigen-depen (C and D) Circles represent the e Th cells, isolated from C57BL/6 mice tative example of five independent e shown. Lines represent the best fit of A B One-Loop Model, Can Explain Kinetics of T-Bet, IL-12Rb2, during Th1 Cell Priming loop model. model, integrating the regulatory sion of IL-12Rb2. C D T-bet, IL-12Rb2, and IFN-g mrna over 5 days of Th1 cell differentiatio (C) or the two-loop (D) model. To test the model, mrna time-course measurements have been collected after inducing the Th1 differentiation ODE model is fit to the experimental data using annealing Figure: From (Schulz et al, 2009) r r r flow cytometry in the early phase and in the late phase in WT and Ifngr / cells in the presence and absence of IL-12 and tion of IL-12Rb2 and advance the second wave of T-bet r (Figures 4D and 4E, left). Usi block activation of the calcine pathway, we found that ex inhibition of antigen-depende after 24 hr of TCR stimulation the predicted effect (Figures 4 right). It has been shown previ the induction of IFN-g and IL T-bet could be indirect and m T-bet-dependent repression o r which otherwise would repr (Usui et al., 2006). In the pr T-bet, elevated Stat4 would enh signaling, induce IFN-g, and IL-12Rb2 chain expression. We Stat4 and GATA-3 in the T cell activation and found, in g correlation with T-bet con (Figure S3). In particular, we fin ence in Stat4 and GATA-3 mrn sion when comparing WT and r cient cells in the first 48 hr of s a phase when T-bet and IFN-g are vastly different (Figures 4A Similarly, at hr after th stimulation, comparison of ce in the presence and absence of IL-12 showed that, al expression of T-bet and IL-12Rb2 differed greatly,
40 In Model C, induction of T-bet by IL-12 was added to Mod Sequential Th1 Th1 Cell Cell Induction Induction by IFN-g by IFN-g and and IL-12 IL-12 T helper cell differentiation example (6) Model C Ag (t ) IFN Re c A A B K1B Ag (t ) K2 IFN Pr ot K3 Re c Pr Figure ot Figure 3. The 3. One-Loop One-Loop Mod Tbet Pr ot f Kinetics Kinetics of T- o Re c 4 Immunity, Volume 30 K8 Tbet Pr ot during during Th1 Cell Th1 (A) Network (A) Network struc loop model. loop model. In Model A (one-loop model), T-bet (B) Schematic (B) r In Model is induced D (two-loop by synergistic model), action antigen-dependent of IFN- and repr ant model, model, integratin signals and IL-12R 2 is induced by dependent T-bet: induction of T-bet was added to model described described A: in Figi T-bet T-bet expression sion of sion IL-12Rb2. of Model A (One loop model) Model D (Two - loop model) (C and (C D) andcircle D) C C T-bet, D D T-bet, IL-12Rb2, Ag (t ) IFN Pr ot Ag (t ) IFN Pr ot Re c over 5over days 5 days of T ftbet 2 Pr ot ftbet 2 3 K Ag (t ) K IFN Th cells, Th cells, isolated 1 2 Pr ot K1 Ag (t ) K2 IFN Pr ot K3 Re c Pr ot tative tative example examof Tbet Pr ot Tbet K shown. shown. LinesLine repr fre c 4 Pr ot 4 f Re c 4 K Tbet (C) or(c) theor two-loo the 8 Pr ot K8 Tbet Pr ot K4 Ag (t ) Pr ot ot f Something wrong with thetbet model? 2 another network 3 model Pr Figure: From (Schulz et al, 2009) In Model B, an additional Term was added, describing repression of IL-12R 2 by ant signals (Ag): tion tion of IL-12Rb of the secon w the second r r r r (Figures (Figures 4D an 4 Model B block block activatio pathway, pathway, we Ag (t ) IFN 9 Pr ot ftbet 2 inhibition inhibition of a K1 Ag (t ) K2 IFN Pr ot after after 24 hr24ofht
41 Sequential Th1 Cell Induction by IFN-g and IL-12 T helper cell differentiation example (7) A B C D Figure 3. The Two-Loop Model, Not the One-Loop Model, Can Explain Expression Kinetics of T-Bet, IL-12Rb2, and IFN-g during Th1 Cell Priming (A) Network structure of the literature-based oneloop model. (B) Schematic representation of the two-loop model, integrating the regulatory interactions described in Figures 1 and 2: IL-12-dependent T-bet expression and antigen-dependent repression of IL-12Rb2. (C and D) Circles represent the expression of T-bet, IL-12Rb2, and IFN-g mrna, measured over 5 days of Th1 cell differentiation with naive Th cells, isolated from C57BL/6 mice. A representative example of five independent experiments is shown. Lines represent the best fit of the one-loop (C) or the two-loop (D) model. tion of IL-12Rb2 and advanced onset of the second wave of T-bet expression r r (Figures 4D and 4E, left). Using CsA to block activation of the calcineurin-nfat pathway, we found that experimental inhibition of antigen-dependent signals after 24 hr of TCR stimulation resulted in the predicted effect (Figures 4D and 4E, right). It has been shown previously that the induction of IFN-g and IL-12Rb2 by T-bet could be indirect and mediated by T-bet-dependent repression of GATA-3, r r which otherwise would repress Stat4 (Usui et al., 2006). In the presence of T-bet, elevated Stat4 would enhance IL-12 signaling, induce IFN-g, and increase IL-12Rb2 chain expression. We measured Stat4 and GATA-3 in the course of T cell activation and found, in general, no correlation with T-bet concentrations (Figure S3). In particular, we find no difference in Stat4 and GATA-3 mrna expression when comparing WT and Ifngr-deficient cells in the first 48 hr of stimulation, r r a phase when T-bet and IFN-g expression are vastly different (Figures 4A and 4C). Similarly, at hr after the onset of Figure: From (Schulz et al, 2009) stimulation, comparison of cells cultured flow cytometry in the early phase and in the late phase in WT in the presence and absence of IL-12 showed that, although the /
42 T helper cell differentiation example (8) Immunity The ODE model is predictive Sequential Th1 Cell Induction by IFN-g and IL-12 A B r r Figure 4. The Two-Loop Model Can Simulate and Predict Expression Kinetics under Various Experimental Conditions (A C) The left column shows the expression kinetics of T-bet (A), IL-12Rb2 (B), and IFN-g mrna (C) simulated on the basis of the two-loop model (Figure 3B) under Th1 cell-inducing conditions (black) and in the absence of IFN-g (blue) or IL-12 signaling (red). All parameters were simultaneously fitted to experimental data. As shown in the right column, naive Th cells, isolated from C57BL/6 (black and red) or Ifngr / mice (blue), were stimulated for 5 days with IL-12 (black and blue) or IFN-g, and blocking antibodies to IL-12 (red). T-bet (A), IL-12Rb2 (B), and IFN-g mrna kinetics (C) were quantified. The experiment was performed three times; one representative example is shown. (D and E) As shown in the left column, expression kinetics of T-bet (D) and IL-12Rb2 mrna (E) are compared in silico in cells stimulated with antigen for 24 hr (red) or 48 hr (black) under Th1 cellinducing conditions. As shown in the right column, naive Th cells, isolated from C57BL/6 mice, were cultured under Th1 cellinducing conditions, and cyclosporine A (CsA) was added after 24 hr of stimulation as indicated (red). T-bet (D) and IL-12Rb2 mrna (E) were quantified over 5 days of culture. One representative example out of three independent experiments is shown. r r C expression. Through a mathematical model we could show that these regulatory interactions between T-bet, IFN-g, IL-12Rb2, and the antigen stimulus can explain the major features of the expression kinetics during primary activation, such as two-peaked T-bet expression with a single IFN-g peak and delayed induction of IL-12Rb2. Therefore, we propose that the model described here constitutes the core of the gene network r r controlling Th1 cell differentiation. T-bet imprints the Th1 cell phenotype in that histones D are modified and chromatin is remodeled in regulatory regions on the Ifng gene (Chang and Aune, 2005; Hatton et al., 2006; Shnyreva et al., 2004). Here, we show a quantitative dose dependence of the frequency of IFN-g-producing cells in a recall stimulation on the expression of T-bet in the late phase of primary activation (>48 hr). Importantly, this T-bet dose dependence of Th1 cell priming held true for different stimulation r protocols (e.g., with or without IFN-g signaling), r strongly suggesting that T-bet expression in a specific E time window, rather than other factors, is the decisive event. The underlying mechanism might involve epigenetic regulation rather than direct transcriptional activation, given that the observed dose dependence is maximal between 72 and 120 hr and decreased after 144 hr when the recall response is induced. A temporal constraint on Th1 cell priming has been shown previously in that memory expression of the Ifng gene requires entry into the S phase of the first r r cell cycle (Bird et al., 1998; Richter et al., 1999). However, because this occurs at 20 hr of stimulation, expression through a positive feedback loop of IL-12R and T-bet. whereas the T-bet effect is strongest at 96 hr, S phase entry The expression of T-bet in the late phase, and not in the early seems to be only one of a number of necessary events that allow phase, correlates strongly with the frequency of IFN-g-producing or mediate T-bet action. In agreement with previous reports, we cells in a later recall response, and thus this late IL-12-dependent show that Runx3 and Hlx, which are induced by and cooperate T-bet induction governs the imprinting of the Th1 cell for IFN-g re- with T-bet, were specifically expressed in the late phase of Figure: From (Schulz et al, 2009)
43 Figure S6. Best fit by Model A-D T helper cell differentiation example (9) (A) Network structure of Model A, B, C and D. (B) circles: Expression level of T-bet (top), IL-12R 2 (middle) and IFN- (bottom) mrna was measured over 5 days of Th1 differentiation More using networks naive Th cells, models isolated from C57BL/6 amice. fundamental lines: The best fit of Models question A-D. concerns choosing the correct biological network model Figure: From (Schulz et al, 2009)
44 References Fröchlich et al. (2017) Scalable Inference of Ordinary Differential Equation Models of Biochemical Processes, arxiv preprint arxiv: Raue et al. (2013) Lessons learned from quantitative dynamical modeling in systems biology, PLoS One, 8(9), e Raue et al. (2010) Identifiability and observability analysis for experimental design in nonlinear dynamical models, Chaos, 20, Brusselle GG et al, (2011) New insights into the immunology of chronic obstructive pulmonary disease, Lancet, 378(9795): Schulz et al, (2009) Sequential Polarization and Imprinting of Type 1 T Helper Lymphocytes by Interferon-g and Interleukin-12, Immunity, 30,
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 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 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 informationAssociation studies and regression
Association studies and regression CM226: Machine Learning for Bioinformatics. Fall 2016 Sriram Sankararaman Acknowledgments: Fei Sha, Ameet Talwalkar Association studies and regression 1 / 104 Administration
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 informationWelcome to Class 21!
Welcome to Class 21! Introductory Biochemistry! Lecture 21: Outline and Objectives l Regulation of Gene Expression in Prokaryotes! l transcriptional regulation! l principles! l lac operon! l trp attenuation!
More informationRegulation of Gene Expression
Chapter 18 Regulation of Gene Expression PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley with contributions from
More informationIntroduction 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 informationInferring Transcriptional Regulatory Networks from High-throughput Data
Inferring Transcriptional Regulatory Networks from High-throughput Data Lectures 9 Oct 26, 2011 CSE 527 Computational Biology, Fall 2011 Instructor: Su-In Lee TA: Christopher Miles Monday & Wednesday 12:00-1:20
More informationProkaryotic Gene Expression (Learning Objectives)
Prokaryotic Gene Expression (Learning Objectives) 1. Learn how bacteria respond to changes of metabolites in their environment: short-term and longer-term. 2. Compare and contrast transcriptional control
More informationLecture 1: Supervised Learning
Lecture 1: Supervised Learning Tuo Zhao Schools of ISYE and CSE, Georgia Tech ISYE6740/CSE6740/CS7641: Computational Data Analysis/Machine from Portland, Learning Oregon: pervised learning (Supervised)
More informationCOURSE NUMBER: EH 590R SECTION: 1 SEMESTER: Fall COURSE TITLE: Computational Systems Biology: Modeling Biological Responses
DEPARTMENT: Environmental Health COURSE NUMBER: EH 590R SECTION: 1 SEMESTER: Fall 2017 CREDIT HOURS: 2 COURSE TITLE: Computational Systems Biology: Modeling Biological Responses COURSE LOCATION: TBD PREREQUISITE:
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 informationSupplementary Figure 1. Confirmation of REF52-hE2F1p::4NLS-d4Venus reporter cells and characterization of E2F dynamics. (a) Alignment of E2F dynamics
Supplementary Figure 1. Confirmation of REF52-hE2F1p::4NLS-d4Venus reporter cells and characterization of E2F dynamics. (a) Alignment of E2F dynamics trajectories to endogenous E2F1 mrna expression. Gray
More informationBIOLOGY STANDARDS BASED RUBRIC
BIOLOGY STANDARDS BASED RUBRIC STUDENTS WILL UNDERSTAND THAT THE FUNDAMENTAL PROCESSES OF ALL LIVING THINGS DEPEND ON A VARIETY OF SPECIALIZED CELL STRUCTURES AND CHEMICAL PROCESSES. First Semester Benchmarks:
More informationBi 8 Lecture 11. Quantitative aspects of transcription factor binding and gene regulatory circuit design. Ellen Rothenberg 9 February 2016
Bi 8 Lecture 11 Quantitative aspects of transcription factor binding and gene regulatory circuit design Ellen Rothenberg 9 February 2016 Major take-home messages from λ phage system that apply to many
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 informationPlant Molecular and Cellular Biology Lecture 10: Plant Cell Cycle Gary Peter
Plant Molecular and Cellular Biology Lecture 10: Plant Cell Cycle Gary Peter 9/10/2008 1 Learning Objectives Explain similarities and differences between fungal, mammalian and plant cell cycles Explain
More informationTypes of biological networks. I. Intra-cellurar networks
Types of biological networks I. Intra-cellurar networks 1 Some intra-cellular networks: 1. Metabolic networks 2. Transcriptional regulation networks 3. Cell signalling networks 4. Protein-protein interaction
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 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 information3.B.1 Gene Regulation. Gene regulation results in differential gene expression, leading to cell specialization.
3.B.1 Gene Regulation Gene regulation results in differential gene expression, leading to cell specialization. We will focus on gene regulation in prokaryotes first. Gene regulation accounts for some of
More informationMachine Learning and Computational Statistics, Spring 2017 Homework 2: Lasso Regression
Machine Learning and Computational Statistics, Spring 2017 Homework 2: Lasso Regression Due: Monday, February 13, 2017, at 10pm (Submit via Gradescope) Instructions: Your answers to the questions below,
More informationREVIEW SESSION. Wednesday, September 15 5:30 PM SHANTZ 242 E
REVIEW SESSION Wednesday, September 15 5:30 PM SHANTZ 242 E Gene Regulation Gene Regulation Gene expression can be turned on, turned off, turned up or turned down! For example, as test time approaches,
More informationIntroduction. Gene expression is the combined process of :
1 To know and explain: Regulation of Bacterial Gene Expression Constitutive ( house keeping) vs. Controllable genes OPERON structure and its role in gene regulation Regulation of Eukaryotic Gene Expression
More informationBayesian Inference and the Symbolic Dynamics of Deterministic Chaos. Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2
How Random Bayesian Inference and the Symbolic Dynamics of Deterministic Chaos Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2 1 Center for Complex Systems Research and Department of Physics University
More informationLogistic Regression. William Cohen
Logistic Regression William Cohen 1 Outline Quick review classi5ication, naïve Bayes, perceptrons new result for naïve Bayes Learning as optimization Logistic regression via gradient ascent Over5itting
More informationSimplicity is Complexity in Masquerade. Michael A. Savageau The University of California, Davis July 2004
Simplicity is Complexity in Masquerade Michael A. Savageau The University of California, Davis July 2004 Complexity is Not Simplicity in Masquerade -- E. Yates Simplicity is Complexity in Masquerade One
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 informationFSC-W FSC-H CD4 CD62-L
Supplementary Fig. 1 a SSC-A FSC-A FSC-W FSC-H SSC-W SSC-H CD4 CD62-L b SSC-A FSC-A FSC-W FSC-A FSC-A 7-AAD FSC-A CD4 IL-9 CD4 c SSC-A FSC-A FSC-W FSC-H SSC-W SSC-H 7-AAD KI67 Annexin-V 7-AAD d I L -5
More informationDynamics of the Mixed Feedback Loop Integrated with MicroRNA
The Second International Symposium on Optimization and Systems Biology (OSB 08) Lijiang, China, October 31 November 3, 2008 Copyright 2008 ORSC & APORC, pp. 174 181 Dynamics of the Mixed Feedback Loop
More informationProkaryotic Regulation
Prokaryotic Regulation Control of transcription initiation can be: Positive control increases transcription when activators bind DNA Negative control reduces transcription when repressors bind to DNA regulatory
More informationMathematical model of IL6 signal transduction pathway
Mathematical model of IL6 signal transduction pathway Fnu Steven, Zuyi Huang, Juergen Hahn This document updates the mathematical model of the IL6 signal transduction pathway presented by Singh et al.,
More informationName Period The Control of Gene Expression in Prokaryotes Notes
Bacterial DNA contains genes that encode for many different proteins (enzymes) so that many processes have the ability to occur -not all processes are carried out at any one time -what allows expression
More informationGENE REGULATION AND PROBLEMS OF DEVELOPMENT
GENE REGULATION AND PROBLEMS OF DEVELOPMENT By Surinder Kaur DIET Ropar Surinder_1998@ yahoo.in Mob No 9988530775 GENE REGULATION Gene is a segment of DNA that codes for a unit of function (polypeptide,
More informationWarm-Up. Explain how a secondary messenger is activated, and how this affects gene expression. (LO 3.22)
Warm-Up Explain how a secondary messenger is activated, and how this affects gene expression. (LO 3.22) Yesterday s Picture The first cell on Earth (approx. 3.5 billion years ago) was simple and prokaryotic,
More informationProkaryotic Gene Expression (Learning Objectives)
Prokaryotic Gene Expression (Learning Objectives) 1. Learn how bacteria respond to changes of metabolites in their environment: short-term and longer-term. 2. Compare and contrast transcriptional control
More informationIntroduction. 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 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 informationBig Idea 3: Living systems store, retrieve, transmit and respond to information essential to life processes. Tuesday, December 27, 16
Big Idea 3: Living systems store, retrieve, transmit and respond to information essential to life processes. Enduring understanding 3.B: Expression of genetic information involves cellular and molecular
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 informationLatent Variable models for GWAs
Latent Variable models for GWAs Oliver Stegle Machine Learning and Computational Biology Research Group Max-Planck-Institutes Tübingen, Germany September 2011 O. Stegle Latent variable models for GWAs
More information56:198:582 Biological Networks Lecture 9
56:198:582 Biological Networks Lecture 9 The Feed-Forward Loop Network Motif Subgraphs in random networks We have discussed the simplest network motif, self-regulation, a pattern with one node We now consider
More informationRegulation of Gene Expression
Chapter 18 Regulation of Gene Expression Edited by Shawn Lester PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley
More informationLecture 2 Machine Learning Review
Lecture 2 Machine Learning Review CMSC 35246: Deep Learning Shubhendu Trivedi & Risi Kondor University of Chicago March 29, 2017 Things we will look at today Formal Setup for Supervised Learning Things
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 informationIntroduction to Bioinformatics
CSCI8980: Applied Machine Learning in Computational Biology Introduction to Bioinformatics Rui Kuang Department of Computer Science and Engineering University of Minnesota kuang@cs.umn.edu History of Bioinformatics
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 informationInferring Protein-Signaling Networks
Inferring Protein-Signaling Networks Lectures 14 Nov 14, 2011 CSE 527 Computational Biology, Fall 2011 Instructor: Su-In Lee TA: Christopher Miles Monday & Wednesday 12:00-1:20 Johnson Hall (JHN) 022 1
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 informationSimulation of Gene Regulatory Networks
Simulation of Gene Regulatory Networks Overview I have been assisting Professor Jacques Cohen at Brandeis University to explore and compare the the many available representations and interpretations of
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 informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationInferring biomarkers for Mycobacterium avium subsp. paratuberculosis infection and disease progression in cattle using experimental data
Inferring biomarkers for Mycobacterium avium subsp. paratuberculosis infection and disease progression in cattle using experimental data 1,2,4, Gesham Magombedze, Tinevimbo Shiri, Shigetoshi Eda 4, 5,
More informationEukaryotic Gene Expression
Eukaryotic Gene Expression Lectures 22-23 Several Features Distinguish Eukaryotic Processes From Mechanisms in Bacteria 123 Eukaryotic Gene Expression Several Features Distinguish Eukaryotic Processes
More informationStochastic simulations
Stochastic simulations Application to molecular networks Literature overview Noise in genetic networks Origins How to measure and distinguish between the two types of noise (intrinsic vs extrinsic)? What
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 informationGene Autoregulation via Intronic micrornas and its Functions
Gene Autoregulation via Intronic micrornas and its Functions Carla Bosia Department of Theoretical Physics University of Torino and INFN, Italy cbosia@to.infn.it Annecy-le-Vieux 20-22/10/2010 OUTLINE microrna
More informationComputational Genomics. Reconstructing dynamic regulatory networks in multiple species
02-710 Computational Genomics Reconstructing dynamic regulatory networks in multiple species Methods for reconstructing networks in cells CRH1 SLT2 SLR3 YPS3 YPS1 Amit et al Science 2009 Pe er et al Recomb
More informationInferring Protein-Signaling Networks II
Inferring Protein-Signaling Networks II Lectures 15 Nov 16, 2011 CSE 527 Computational Biology, Fall 2011 Instructor: Su-In Lee TA: Christopher Miles Monday & Wednesday 12:00-1:20 Johnson Hall (JHN) 022
More informationDiscrete Latent Variable Models
Discrete Latent Variable Models Stefano Ermon, Aditya Grover Stanford University Lecture 14 Stefano Ermon, Aditya Grover (AI Lab) Deep Generative Models Lecture 14 1 / 29 Summary Major themes in the course
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationComputational Cell Biology Lecture 4
Computational Cell Biology Lecture 4 Case Study: Basic Modeling in Gene Expression Yang Cao Department of Computer Science DNA Structure and Base Pair Gene Expression Gene is just a small part of DNA.
More informationGenome 541! Unit 4, lecture 2! Transcription factor binding using functional genomics
Genome 541 Unit 4, lecture 2 Transcription factor binding using functional genomics Slides vs chalk talk: I m not sure why you chose a chalk talk over ppt. I prefer the latter no issues with readability
More informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
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 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 informationBioControl - Week 6, Lecture 1
BioControl - Week 6, Lecture 1 Goals of this lecture Large metabolic networks organization Design principles for small genetic modules - Rules based on gene demand - Rules based on error minimization Suggested
More informationPlant Molecular and Cellular Biology Lecture 8: Mechanisms of Cell Cycle Control and DNA Synthesis Gary Peter
Plant Molecular and Cellular Biology Lecture 8: Mechanisms of Cell Cycle Control and DNA Synthesis Gary Peter 9/10/2008 1 Learning Objectives Explain why a cell cycle was selected for during evolution
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 informationRegulation of Transcription in Eukaryotes. Nelson Saibo
Regulation of Transcription in Eukaryotes Nelson Saibo saibo@itqb.unl.pt In eukaryotes gene expression is regulated at different levels 1 - Transcription 2 Post-transcriptional modifications 3 RNA transport
More informationGene regulation I Biochemistry 302. Bob Kelm February 25, 2005
Gene regulation I Biochemistry 302 Bob Kelm February 25, 2005 Principles of gene regulation (cellular versus molecular level) Extracellular signals Chemical (e.g. hormones, growth factors) Environmental
More informationStem Cell Reprogramming
Stem Cell Reprogramming Colin McDonnell 1 Introduction Since the demonstration of adult cell reprogramming by Yamanaka in 2006, there has been immense research interest in modelling and understanding the
More informationChapter 6- An Introduction to Metabolism*
Chapter 6- An Introduction to Metabolism* *Lecture notes are to be used as a study guide only and do not represent the comprehensive information you will need to know for the exams. The Energy of Life
More informationIntroduction to Machine Learning. Lecture 2
Introduction to Machine Learning Lecturer: Eran Halperin Lecture 2 Fall Semester Scribe: Yishay Mansour Some of the material was not presented in class (and is marked with a side line) and is given for
More informationUnderstanding Science Through the Lens of Computation. Richard M. Karp Nov. 3, 2007
Understanding Science Through the Lens of Computation Richard M. Karp Nov. 3, 2007 The Computational Lens Exposes the computational nature of natural processes and provides a language for their description.
More information4. Why not make all enzymes all the time (even if not needed)? Enzyme synthesis uses a lot of energy.
1 C2005/F2401 '10-- Lecture 15 -- Last Edited: 11/02/10 01:58 PM Copyright 2010 Deborah Mowshowitz and Lawrence Chasin Department of Biological Sciences Columbia University New York, NY. Handouts: 15A
More informationRegulation and signaling. Overview. Control of gene expression. Cells need to regulate the amounts of different proteins they express, depending on
Regulation and signaling Overview Cells need to regulate the amounts of different proteins they express, depending on cell development (skin vs liver cell) cell stage environmental conditions (food, temperature,
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 informationUnit 3: Control and regulation Higher Biology
Unit 3: Control and regulation Higher Biology To study the roles that genes play in the control of growth and development of organisms To be able to Give some examples of features which are controlled
More informationA 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 informationLecture 4. f X T, (x t, ) = f X,T (x, t ) f T (t )
LECURE NOES 21 Lecture 4 7. Sufficient statistics Consider the usual statistical setup: the data is X and the paramter is. o gain information about the parameter we study various functions of the data
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationDynamic Causal Modelling for fmri
Dynamic Causal Modelling for fmri André Marreiros Friday 22 nd Oct. 2 SPM fmri course Wellcome Trust Centre for Neuroimaging London Overview Brain connectivity: types & definitions Anatomical connectivity
More informationGenome 541! Unit 4, lecture 3! Genomics assays
Genome 541! Unit 4, lecture 3! Genomics assays Much easier to follow with slides. Good pace.! Having the slides was really helpful clearer to read and easier to follow the trajectory of the lecture.!!
More informationModeling Cellular Networks
06_4774 12/7/06 3:02 PM Page 151 CHAPTER 6 Modeling Cellular Networks Tae Jun Lee, Dennis Tu, Chee Meng Tan, and Lingchong You 6.1 Introduction Systems-level understanding of cellular dynamics is important
More informationInferring Transcriptional Regulatory Networks from Gene Expression Data II
Inferring Transcriptional Regulatory Networks from Gene Expression Data II Lectures 9 Oct 26, 2011 CSE 527 Computational Biology, Fall 2011 Instructor: Su-In Lee TA: Christopher Miles Monday & Wednesday
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 informationControl of Gene Expression
Control of Gene Expression Mechanisms of Gene Control Gene Control in Eukaryotes Master Genes Gene Control In Prokaryotes Epigenetics Gene Expression The overall process by which information flows from
More informationNetwork Biology: Understanding the cell s functional organization. Albert-László Barabási Zoltán N. Oltvai
Network Biology: Understanding the cell s functional organization Albert-László Barabási Zoltán N. Oltvai Outline: Evolutionary origin of scale-free networks Motifs, modules and hierarchical networks Network
More informationMeasuring TF-DNA interactions
Measuring TF-DNA interactions How is Biological Complexity Achieved? Mediated by Transcription Factors (TFs) 2 Regulation of Gene Expression by Transcription Factors TF trans-acting factors TF TF TF TF
More informationAssembling and Testing Physically Based Models of Observed Complex Systems
Assembling and Testing Physically Based Models of Observed Complex Systems with Daniel Rey, Nirag Kadakia, Daniel Breen, Sasha Shirman, Eve Armstrong, Dan Margoliash, Arij Daou, Gert Cauwenberghs, J. Wang,
More informationREGULATION OF GENE EXPRESSION. Bacterial Genetics Lac and Trp Operon
REGULATION OF GENE EXPRESSION Bacterial Genetics Lac and Trp Operon Levels of Metabolic Control The amount of cellular products can be controlled by regulating: Enzyme activity: alters protein function
More informationIntroduction to gradient descent
6-1: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction to gradient descent Derivation and intuitions Hessian 6-2: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction Our
More informationLong noncoding RNAs in T lymphocytes
Review Long noncoding RNAs in T lymphocytes Thomas M. Aune,*,1 Phillip S. Crooke III, and Charles F. Spurlock III* Departments of *Medicine and Mathematics, Vanderbilt University, Nashville, Tennessee,
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 informationIdentifying Signaling Pathways
These slides, excluding third-party material, are licensed under CC BY-NC 4.0 by Anthony Gitter, Mark Craven, Colin Dewey Identifying Signaling Pathways BMI/CS 776 www.biostat.wisc.edu/bmi776/ Spring 2018
More information6.047 / Computational Biology: Genomes, Networks, Evolution Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 6.047 / 6.878 Computational Biology: Genomes, Networks, Evolution Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationSeries 6, May 14th, 2018 (EM Algorithm and Semi-Supervised Learning)
Exercises Introduction to Machine Learning SS 2018 Series 6, May 14th, 2018 (EM Algorithm and Semi-Supervised Learning) LAS Group, Institute for Machine Learning Dept of Computer Science, ETH Zürich Prof
More informationUNIVERSITY OF YORK. BA, BSc, and MSc Degree Examinations Department : BIOLOGY. Title of Exam: Molecular microbiology
Examination Candidate Number: Desk Number: UNIVERSITY OF YORK BA, BSc, and MSc Degree Examinations 2017-8 Department : BIOLOGY Title of Exam: Molecular microbiology Time Allowed: 1 hour 30 minutes Marking
More information