Lecture Network analysis for biological systems
|
|
- Morris Ramsey
- 5 years ago
- Views:
Transcription
1 Lecture Network analysis for biological systems Anja Bråthen Kristoffersen
2 Biological Networks Gene regulatory network: two genes are connected if the expression of one gene modulates expression of another one by either activation or inhibition Protein interaction network: proteins that are connected in physical interactions or metabolic and signaling pathways of the cell Metabolic network: metabolic products and substrates that participate in one reaction Statistical bioinformatics 3
3 What is Gene Regulatory Network? Gene regulatory networks (GRNs) are the on-off switches of a cell operating at the gene level. Two genes are connected if the expression of one gene modulates expression of another one by either activation or inhibition Statistical bioinformatics 4
4 Simplified Representation of Gene Regulatory Network A gene regulatory network can be represented by a directed graph Node represents a gene Directed edge stands for the modulation (regulation) of one node by another: e.g. arrow from gene X to gene Y means gene X affects expression of gene Y Statistical bioinformatics 5
5 Why study Gene Regulatory Network Genes are not independent They regulate each other and act collectively This collective behavior can be observed using microarray Some genes control the response of the cell to changes in the environment by regulating other genes; Potential discovery of triggering mechanism and treatments for disease Statistical bioinformatics 6
6 Network Modeling techniques Boolean network (BN) Bayesian belief network Metabolic network modeling methods Statistical bioinformatics 7
7 Boolean network modeling Boolean: either true or false (1 or 0) Binarization reduces the noise in biological data captures the dynamic behavior in complex systems need a threshold value leads to loss of information Genes are modeled as switch like dynamic elements either on or off Statistical bioinformatics 8
8 Boolean network consist of A set of genes. A set of Boolean functions F = f i (x 1, x 2,, x n ) the function is described with three boolean operators AND / && / & / OR / / / NOT /! / ~ Statistical bioinformatics 9
9 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } x 2 x 1 x 3 Wiring diagram Input (t-1) Truth table Output(t) Statistical bioinformatics 10
10 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } x 2 x 1 x 3 Input (t-1) Output(t) Statistical bioinformatics 11
11 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } x 2 x 1 x 3 Input (t-1) Output(t) Statistical bioinformatics 12
12 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } x 2 x 1 x 3 Input (t-1) Output(t) Statistical bioinformatics 13
13 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } x 2 x 1 x 3 Input (t-1) Output(t) Statistical bioinformatics 14
14 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } Input (t-1) Output(t) x 2 x 1 x 3 Statistical bioinformatics 15
15 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } Input (t-1) Output(t) x 2 x 1 x 3 Statistical bioinformatics 16
16 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } Input (t-1) Output(t) x x 1 x Statistical bioinformatics 17
17 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } Input (t-1) Output(t) x x 1 x Statistical bioinformatics 18
18 Example Graph (G) with 3 genes Given network G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } Input (t-1) Output(t) x 2 x x 3 State transition Statistical bioinformatics 19
19 R code truth table G(V,F), V = {x 1, x 2, x 3 } F = {f 1 = x 2 & x 3, f 2 = x 1, f 3 = x 2 } # f1 state x(t-1) is a vector x with three elements # want to find state x(t) here called y also with three elements y <- rep(na, 3) # allocate space for three elements in y y[1] <- x[2] && x[3] y[2] <- x[1] y[3] <- x[2] # test it by leting x <- c(0,1,0) and run the code over, according to the # truth table should y be 0,0,1. x <- c(0,1,0) y <- rep(na, 3) # allocate space for three elements in y y[1] <- x[2] && x[3] y[2] <- x[1] y[3] <- x[2] Statistical bioinformatics 20
20 Make a function of the truth table R code output <- function(x){ y <- rep(na, 3) y[1] <- x[2] && x[3] y[2] <- x[1] y[3] <- x[2] return(y) } output(c(0,0,1)) [1] Statistical bioinformatics 21
21 R code: Make the output table based on a input table input <- rbind(c(0,0,0), c(0,0,1), c(0,1,0), c(0,1,1), c(1,0,0), c(1,0,1), c(1,1,0), c(1,1,1)) y <- matrix(na, ncol = ncol(input), nrow= nrow(input)) for(i in 1:nrow(y)){ y[i,] <- output(input[i,]) } Statistical bioinformatics 22
22 Search for Boolean functions Chi-square testing-based search Kim H, Lee JK, Park T. (2007). Boolean networks using the chi-square test for inferring large-scale gene regulatory networks. BMC Bioinformatics, 8:37. G1 G2 G3 G4 T T T T T T T T T januar 2014 Statistical bioinformatics 23 T Binary data set from simple network: Four nodes (G1, G2, G3, G4) and 10 time points
23 Observe which gene is on/off the time point before If we want to find a Boolean function, f 4 for node G4, we have to test the independency between G4 at time t and all nodes (G1, G2, G3, G4) at time t-1. We use a chi square distribution and compare expected with observed Statistical bioinformatics 24
24 2 x 2 Contingency table G4 t 0 1 G1 t G4 t 0 1 G2 t Find both the expected and the observed contigency tables G4 t 0 1 G3 t G4 t 0 1 G4 t Statistical bioinformatics 25
25 2 x 2 Contingency table, expected Assume independence between Gi and Gj Only depend on background distribution G1 t G4 t 0 P(G4=0)*P(G1=0) P(G4=0)*P(G1=1) 1 P(G4=1)*P(G1=0) P(G4=1)*P(G1=1) Similar for the other pars of Gi and Gj Statistical bioinformatics 26
26 2 x 2 Contingency table, observed G4 t 0 1 G4 t 0 1 G1 t G3 t G4 t 0 1 G4 t 0 1 G2 t G4 t G1 G2 G3 G4 T T T T T T T T T T Statistical bioinformatics 27
27 2 x 2 Contingency table, observed G1 t G4 t G4 t 0 1 G3 t G4 t 0 1 G4 t 0 1 G2 t G4 t G1 G2 G3 G4 T T T T T T T T T T Statistical bioinformatics 28
28 2 x 2 Contingency table, observed G1 t G4 t G3 t G4 t G2 t G4 t G4 t G4 t G1 G2 G3 G4 T T T T T T T T T T Statistical bioinformatics 29
29 chi square distribution Compare expected with observed Statistical bioinformatics 30
30 Chi-Square Test Result of independence for G4 t node G1 t-1 G2 t-1 G3 t-1 G4 t p-value Statistical bioinformatics 31
31 R code, 2 x 2 Contingency table d1 <- read.table("c:/users/anjab/desktop/infstk/simpledataset4genes10timepoints.txt", header = T, sep = "\t") #I got my dataset read in with the rownames in the first coloumn, #I did not like that so I changed it d2 <- d1[,2:5] rownames(d2) <- d1[,1] #always check that you read it in correctly. head(d2) G1 G2 G3 G4 T T T T T T Statistical bioinformatics 32
32 R code, 2 x 2 contingency table G1 t G4 t 0 a b 1 c d # start with gene 4 and make the 2 x 2 contingency table for gene 1 # find the number a. n <- nrow(d2) posible <- which(d2[2:n,4] == 0) posible [1] a <- sum(length(which(d2[posible, 1] == 0))) a [1] 4 G1 G4 T1 0 1 T2 1 1 T3 1 0 T4 0 0 T5 0 1 T6 1 0 T7 0 0 T8 0 0 T9 0 0 Statistical bioinformatics 33 T10 0 0
33 R code, 2 x 2 contingency table # find the number a, b, c and d. n <- nrow(d2) posible0 <- which(d2[2:n,4] == 0) posible1 <- which(d2[2:n,4] == 1) a <- sum(length(which(d2[posible0, 1] == 0))) b <- sum(length(which(d2[posible0, 1] == 1))) c1 <- sum(length(which(d2[posible1, 1] == 0))) d <- sum(length(which(d2[posible1, 1] == 1))) conttable <- matrix(c(a,b,c1,d), ncol = 2, byrow = T) conttable [,1] [,2] [1,] 4 3 [2,] 2 0 Statistical bioinformatics G1 t G4 t 0 a b 1 c d G1 G4 T1 0 1 T2 1 1 T3 1 0 T4 0 0 T5 0 1 T6 1 0 T7 0 0 T8 0 0 T9 0 0 T10 0 0
34 chisq.test(conttable) > chisq.test(conttable) Pearson's Chi-squared test with Yates' continuity correction data: conttable X-squared = , df = 1, p-value = Warning message: In chisq.test(conttable) : Chi-squared approximation may be incorrect NB, we have a very little dataset, how can we get rid of the warning message?. Look at help(chisq.test) Statistical bioinformatics 35
35 help(chisq.test) You find out that the p-value can be simulated using Monte Carlo simulation. This is for us with a small dataset, a good option. Statistical bioinformatics 36
36 chisq.test(, simulate.p.value = T) > chisq.test(conttable, simulate.p.value = T) Pearson's Chi-squared test with simulated p-value (based on 2000 replicates) data: conttable X-squared = , df = NA, p-value = > chisq.test(conttable, simulate.p.value = T) Pearson's Chi-squared test with simulated p-value (based on 2000 replicates) data: conttable X-squared = , df = NA, p-value = Statistical bioinformatics 37
37 Chi-Square Test Result of independence for all genes G1 t-1 G2 t-1 G3 t-1 G4 t-1 p-value G1 t p-value G2 t p-value G3 t p-value G4 t Statistical bioinformatics 38
38 Make a function that takes two vectors, x (gene at time t) and y (gene at time t-1) chisqres <- function(x,y){ n <- length(x) posible0 <- which(y[2:n] == 0) posible1 <- which(y[2:n] == 1) a <- sum(length(which(x[posible0] == 0))) b <- sum(length(which(x[posible0] == 1))) c1 <- sum(length(which(x[posible1] == 0))) d <- sum(length(which(x[posible1] == 1))) counttable <- matrix(c(a,b,c1,d), ncol = 2, byrow = T) chisq.test(counttable,, simulate.p.value = T)$p.value } 11. januar 2014 Statistical bioinformatics 39
39 Use chisqres() res <- matrix(na, ncol(d2), ncol(d2)) for(i in 1:ncol(d2)){ } for(j in 1:ncol(d2)){ } res[i,j] <- chisqres(d2[,i], d2[,j]) colnames(res) <- paste(colnames(d2), "t-1", sep = "") rownames(res) <- paste(colnames(d2), "t", sep = "") round(res, 3) G1t-1 G2t-1 G3t-1 G4t-1 G1t G2t G3t This only tells us that G1t is dependent on the state G4t-1 had, not which type of dependence. G4t
40 Goes further and calculate tables with three genes: a gene at time t and two genes at time t - 1. Eg. G1t, G3t - 1, G4t - 1 Statistical bioinformatics 41
41 Probabilistic Boolean Network Allow multiple Boolean functions at each node with different probabilities F = {{(f 11, c 11 ),..., (f 1k1, c 1k1 )},..., {(f n1, c n1 ),..., (f nkn, c nkn )}} where k1 is the number of different transition functions for gene 1 the sum of all transition probabilities c 11 to c 1k1 is always 1 Statistical bioinformatics 42
42 Example: Probabilistic Boolean Network Given three genes, and the transition functions: F = { F 1 = {(f 11, c 11 ), (f 12, c 12 )}, F 2 = {(f 21, 1)}, F 3 = {(f 31, c 31 ), (f 32, c 32 )}} Given the truth tabel and probabilities for each transition function There are 2*1*2 different set of transition functions that can be used. They are all listed in tabel K Statistical bioinformatics 43
43 Statistical bioinformatics 44
44 Transition probability matrix P 1 P 2 P 3 P 4 P 1 = c 11 *1* c Statistical bioinformatics
45 R code x1x2x3 f11 f12 f21 f31 f32 State State State State State State State State truthtable1 <- read.table("m:/undervisning/statistical bioinformatics/datasets used/truthtableex.txt", header = T) truthtable <- truthtable1[,2:6] rownames(truthtable) <- substr(truthtable1[,1], 6,8) possiblecomb <- nrow(truthtable) c11 <- 0.6 #probability of function f11 being used c12 <- 1 - c11 c21 <- 1 c31 <- 0.5 c32 <- 1 - c31 A <- matrix(0, nrow = possiblecomb, ncol = possiblecomb) rownames(a) <- rownames(truthtable) colnames(a) <- rownames(truthtable) Statistical bioinfomratics 46
46 possiblemodelstruthtablecolumn <- rbind(c(1,3,4), c(1,3,5), c(2,3,4), c(2,3,5)) probmodels <- c(c11*c21*c31, c11*c21*c32, c12*c21*c31, c12*c21*c32) #probability that each of the possiblemodels are used for(i in 1:nrow(A)){ from <- rownames(a)[i] for (j in 1:nrow(possibleModelsTruthTableColumn)){ modelj <- possiblemodelstruthtablecolumn[j,] to <- paste(truthtable[i,modelj[1]], truthtable[i,modelj[2]], truthtable[i,modelj[3]], sep = "") probmodelsj <- probmodels[j] A[from, to] <- A[from, to] + probmodelsj } } Statistical bioinfomratics 47
47 Synchronous Boolean networks Assume that all genes are updated at the same time This simplification facilitates the analysis of the networks We have until now looked at such simplified networks Statistical bioinformatics 48
48 Asynchronous Boolean networks at each point of time t, only one of the transition functions f i F is chosen at random, and the corresponding Boolean variable is updated. Statistical bioinformatics 49
49 Provides tools for assembling analyzing visualizing Synchronous, asynchronous and probabilistic Boolean networks install.packages("boolnet") library(boolnet) Statistical bioinformatics 50
50 BoolNet, syntaxes targets, factors or targets, factors, probabilities Target is the gene that is effected Factors are those genes effecting it Probabilities occurs when it is more then one transition function Statistical bioinformatics 51
51 BoolNet, syntaxes Example CycD is an input, considered as constant. Translated into a transition rule: CycD, CycD Statistical bioinformatics 52
52 BoolNet, syntaxes Example Rb is expressed if all the genes CycA, CycB, CycD and CycE is absence; it can be expressed in the presence of CycE or CycA if their inhibitory activity is blocked by p27. Translated into a transition rule: First part:! CycA &! CycB &! CycD &! CycE Second part: p27 &! CycB &! CycD Together: Rb, (! CycA &! CycB &! CycD &! CycE) (p27 &! CycB &! CycD) Statistical bioinformatics 53
53 Read a network into R Assume that we have the file cellcycle.txt targets, factors CycD, CycD Rb, (! CycA &! CycB &! CycD &! CycE) (p27 &! CycB &! CycD) E2F, (! Rb &! CycA &! CycB) (p27 &! Rb &! CycB) CycE, (E2F &! Rb) CycA, (E2F &! Rb &! Cdc20 &! (Cdh1 & UbcH10)) (CycA &! Rb &! Cdc20 &! (Cdh1 & UbcH10)) p27, (! CycD &! CycE &! CycA &! CycB) (p27 &! (CycE & CycA) &! CycB &! CycD) Cdc20, CycB Cdh1,(! CycA &! CycB) (Cdc20) (p27 &! CycB) UbcH10,! Cdh1 (Cdh1 & UbcH10 & (Cdc20 CycA CycB)) CycB,! Cdc20 &! Cdh1 Read it into R by: cellcycle <- loadnetwork("cellcycle.txt") cellcycle Statistical bioinformatics 54
54 Reconstruct a network from time series A dataset that are already in BoolNet is the yeasttimeseries: To use this data it has to be binarized Statistical bioinformatics 55
55 Binarization, can be done in many ways. BoolNet support three methods: k-means clustering For each gene, k-means clustering are performed to determine a good separation of groups Edge detector This approach first sorts the measurements for each gene. In the sorted measurements, the algorithm searches for differences of two successive values that satisfy a predefined condition Scan statistic The scan statistic assumes that the measurements for each gene are uniformly and independently distributed. The scan statistic shifts a scanning window across the data and decides for each window position whether there is an unusual accumulation of data points based on an approximated test statistic (see Glaz et al.). 56
56 Reconstruct network. Statistical bioinformatics 57
57 How to read the output Fkh2 = <f(clb1){01}> means Clb1(t) Fkh2(t+1) Sic1 = <f(sic1,clb1){0001}> means Sic1(t) Clb(t) Sic1(t+1) Statistical bioinformatics 58
58 How to ead the output Fkh2 = <f(clb1){01}> means Clb1(t) Fkh2(t+1) Sic1 = <f(sic1,clb1){0001}> means Sic1(t) Clb(t) Sic1(t+1) Statistical bioinformatics 59
59 plotnetworkwiring(net) Statistical bioinformatics 60
60 Creating random networks It is desirable to generate artificial networks To study structural properties of Boolean networks To determine the specific properties of biological networks in comparison to arbitrary networks net <- generaterandomnknetwork(n=10, k=3) Statistical bioinformatics 61
61 Attractors Attractors are stable cycles of states in a Boolean network Attractors in models of gene-regulatory networks are expected to be linked to phenotypes All states that lead to a certain attractor form its basin of attraction Statistical bioinformatics 62
62 Simple attractors occur in synchronous Boolean networks consist of a set of states whose synchronous transitions form a cycle. Complex or loose attractors in asynchronous networks usually more than one possible transition for each state in an asynchronous network a complex attractor is formed by two or more overlapping loops. Steady-state attractors are attractors that consist of only one state. All transitions from this state result in the state itself. Statistical bioinformatics 63
63 Statistical bioinformatics 64
64 Perturbation experiments The generation of perturbed copies of a network is a way to test the robustness of structural properties of the networks to noise and mismeasurements. For example, you could assess the relevance of an attractor by checking whether the same attractor is still found when small random changes are applied to the network. If this is the case, it is less likely that the attractor is an artifact of mismeasurements. perturbednet <- perturbnetwork(cellcycle, perturb="functions", method="bitflip") Statistical bioinformatics 65
65 Generate random networks Generate a random network with as many nodes and edges as your original network Find the attractors in this network Perturbate the random network 1000 times, how many times are the original attractors from the random network found in the perturbated networks Repeate 1000 times 66
66 data(cellcycle) perturbednet <- perturbnetwork(cellcycle, perturb="functions", method="bitflip") Statistical bioinformatics 67
67 Boolean network modeling Positive: Explains the dynamic behavior of living systems efficiently (with possible loops!) Boolean algebra provides a rich set of algorithms already available for supervised learning in binary domain, such as logical analysis of data, and Boolean-based classification algorithms Dichotomization to binary values improves accuracy of classification and simplifies the obtained models by reducing the noise level in experimental data Statistical bioinformatics 68
68 Boolean network modeling Negative: Requires heavy computing times to construct a network structure Needs specific time course data that well capture pathway interactions, but often uncertain whether there are such time points and, even so, whether they were captured well by a time-course experiment Needs a relatively large number of time points Lose quantitative information by dichotomization to binary values Statistical bioinformatics 69
Algorithms and methods of the BoolNet R package
Algorithms and methods of the BoolNet R package Christoph Müssel, Martin Hopfensitz, Hans A. Kestler Abstract This document describes the algorithms and methods that were developed or partially adapted
More informationSynchronous state transition graph
Heike Siebert, FU Berlin, Molecular Networks WS10/11 2-1 Synchronous state transition graph (0, 2) (1, 2) vertex set X (state space) edges (x,f(x)) every state has only one successor attractors are fixed
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 informationState-Feedback Control of Partially-Observed Boolean Dynamical Systems Using RNA-Seq Time Series Data
State-Feedback Control of Partially-Observed Boolean Dynamical Systems Using RNA-Seq Time Series Data Mahdi Imani and Ulisses Braga-Neto Department of Electrical and Computer Engineering Texas A&M University
More information5.3 METABOLIC NETWORKS 193. P (x i P a (x i )) (5.30) i=1
5.3 METABOLIC NETWORKS 193 5.3 Metabolic Networks 5.4 Bayesian Networks Let G = (V, E) be a directed acyclic graph. We assume that the vertices i V (1 i n) represent for example genes and correspond to
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 informationBiological networks CS449 BIOINFORMATICS
CS449 BIOINFORMATICS Biological networks Programming today is a race between software engineers striving to build bigger and better idiot-proof programs, and the Universe trying to produce bigger and better
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 information86 Part 4 SUMMARY INTRODUCTION
86 Part 4 Chapter # AN INTEGRATION OF THE DESCRIPTIONS OF GENE NETWORKS AND THEIR MODELS PRESENTED IN SIGMOID (CELLERATOR) AND GENENET Podkolodny N.L. *1, 2, Podkolodnaya N.N. 1, Miginsky D.S. 1, Poplavsky
More informationComputational Genomics. Systems biology. Putting it together: Data integration using graphical models
02-710 Computational Genomics Systems biology Putting it together: Data integration using graphical models High throughput data So far in this class we discussed several different types of high throughput
More informationarxiv: v1 [q-bio.mn] 7 Nov 2018
Role of self-loop in cell-cycle network of budding yeast Shu-ichi Kinoshita a, Hiroaki S. Yamada b a Department of Mathematical Engineering, Faculty of Engeneering, Musashino University, -- Ariake Koutou-ku,
More informationBoolean networks for modeling and analysis of gene regulation
Boolean networks for modeling and analysis of gene regulation Dao Zhou, Christoph Müssel, Ludwig Lausser, Martin Hopfensitz, Michael Kühl, Hans A. Kestler Ulmer Informatik-Berichte Nr. 2009-10 Oktober
More informationIndex. FOURTH PROOFS n98-book 2009/11/4 page 261
FOURTH PROOFS n98-book 2009/11/4 page 261 Index activity, 10, 13 adaptive control, 215 adaptive controller, 215 adaptive immune system, 20 adaptive intervention strategy, 216 adjacency rule, 111 admissible
More informationA&S 320: Mathematical Modeling in Biology
A&S 320: Mathematical Modeling in Biology David Murrugarra Department of Mathematics, University of Kentucky http://www.ms.uky.edu/~dmu228/as320/ These slides were modified from Matthew Macauley s lecture
More informationLogical modeling of the mammalian cell cycle
Logical modeling of the mammalian cell cycle Pauline Traynard 1,2, Adrien Fauré 3, François Fages 2 and Denis Thieffry 1,2 1 Computational Systems Biology team, Institut de Biologie de l Ecole rmale Supérieure
More informationNetwork Dynamics and Cell Physiology. John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute
Network Dynamics and Cell Physiology John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute Signal Transduction Network Hanahan & Weinberg (2000) Gene Expression Signal-Response
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 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 informationA New Method to Build Gene Regulation Network Based on Fuzzy Hierarchical Clustering Methods
International Academic Institute for Science and Technology International Academic Journal of Science and Engineering Vol. 3, No. 6, 2016, pp. 169-176. ISSN 2454-3896 International Academic Journal of
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 informationRandom Boolean Networks
Random Boolean Networks Boolean network definition The first Boolean networks were proposed by Stuart A. Kauffman in 1969, as random models of genetic regulatory networks (Kauffman 1969, 1993). A Random
More informationnutrients growth & division repellants movement
Network Dynamics and Cell Physiology John J. Tyson Department of Biological Sciences & Virginia Bioinformatics Institute Outline 1. Cell Signaling: Physiology 2. Cell Signaling: Molecular Biology 3. Chemical
More informationInferring Update Sequences in Boolean Gene Regulatory Networks
Inferring Update Sequences in Boolean Gene Regulatory Networks Fabien Tarissan a Camilo La Rota b a Complex System Institute (ISC) & CNRS, Palaiseau, France b Complex System Institute (IXXI), Lyon, France
More informationClassification of Random Boolean Networks
Classification of Random Boolean Networks Carlos Gershenson, School of Cognitive and Computer Sciences University of Sussex Brighton, BN1 9QN, U. K. C.Gershenson@sussex.ac.uk http://www.cogs.sussex.ac.uk/users/carlos
More information6.867 Machine learning, lecture 23 (Jaakkola)
Lecture topics: Markov Random Fields Probabilistic inference Markov Random Fields We will briefly go over undirected graphical models or Markov Random Fields (MRFs) as they will be needed in the context
More informationBoolean models of gene regulatory networks. Matthew Macauley Math 4500: Mathematical Modeling Clemson University Spring 2016
Boolean models of gene regulatory networks Matthew Macauley Math 4500: Mathematical Modeling Clemson University Spring 2016 Gene expression Gene expression is a process that takes gene info and creates
More informationClassification of Random Boolean Networks
in Artificial Life VIII, Standish, Abbass, Bedau (eds)(mit Press) 2002. pp 1 8 1 Classification of Random Boolean Networks Carlos Gershenson, School of Cognitive and Computer Sciences University of Sussex
More informationSig2GRN: A Software Tool Linking Signaling Pathway with Gene Regulatory Network for Dynamic Simulation
Sig2GRN: A Software Tool Linking Signaling Pathway with Gene Regulatory Network for Dynamic Simulation Authors: Fan Zhang, Runsheng Liu and Jie Zheng Presented by: Fan Wu School of Computer Science and
More informationLecture 6: Time-Dependent Behaviour of Digital Circuits
Lecture 6: Time-Dependent Behaviour of Digital Circuits Two rather different quasi-physical models of an inverter gate were discussed in the previous lecture. The first one was a simple delay model. This
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 informationComputational Systems Biology
Computational Systems Biology Vasant Honavar Artificial Intelligence Research Laboratory Bioinformatics and Computational Biology Graduate Program Center for Computational Intelligence, Learning, & Discovery
More informationLab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations
Lab 2 Worksheet Problems Problem : Geometry and Linear Equations Linear algebra is, first and foremost, the study of systems of linear equations. You are going to encounter linear systems frequently in
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 informationFCModeler: Dynamic Graph Display and Fuzzy Modeling of Regulatory and Metabolic Maps
FCModeler: Dynamic Graph Display and Fuzzy Modeling of Regulatory and Metabolic Maps Julie Dickerson 1, Zach Cox 1 and Andy Fulmer 2 1 Iowa State University and 2 Proctor & Gamble. FCModeler Goals Capture
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 informationLearning in Bayesian Networks
Learning in Bayesian Networks Florian Markowetz Max-Planck-Institute for Molecular Genetics Computational Molecular Biology Berlin Berlin: 20.06.2002 1 Overview 1. Bayesian Networks Stochastic Networks
More informationThe Monte Carlo Method: Bayesian Networks
The Method: Bayesian Networks Dieter W. Heermann Methods 2009 Dieter W. Heermann ( Methods)The Method: Bayesian Networks 2009 1 / 18 Outline 1 Bayesian Networks 2 Gene Expression Data 3 Bayesian Networks
More informationOptimal control of Boolean control networks
1 Optimal control of Boolean control networks Ettore Fornasini and Maria Elena Valcher Abstract In this paper we address the optimal control problem for Boolean control networks (BCNs). We first consider
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 informationTHE VINE COPULA METHOD FOR REPRESENTING HIGH DIMENSIONAL DEPENDENT DISTRIBUTIONS: APPLICATION TO CONTINUOUS BELIEF NETS
Proceedings of the 00 Winter Simulation Conference E. Yücesan, C.-H. Chen, J. L. Snowdon, and J. M. Charnes, eds. THE VINE COPULA METHOD FOR REPRESENTING HIGH DIMENSIONAL DEPENDENT DISTRIBUTIONS: APPLICATION
More informationDESIGN OF EXPERIMENTS AND BIOCHEMICAL NETWORK INFERENCE
DESIGN OF EXPERIMENTS AND BIOCHEMICAL NETWORK INFERENCE REINHARD LAUBENBACHER AND BRANDILYN STIGLER Abstract. Design of experiments is a branch of statistics that aims to identify efficient procedures
More informationSelf Similar (Scale Free, Power Law) Networks (I)
Self Similar (Scale Free, Power Law) Networks (I) E6083: lecture 4 Prof. Predrag R. Jelenković Dept. of Electrical Engineering Columbia University, NY 10027, USA {predrag}@ee.columbia.edu February 7, 2007
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 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 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 informationThe Role of Network Science in Biology and Medicine. Tiffany J. Callahan Computational Bioscience Program Hunter/Kahn Labs
The Role of Network Science in Biology and Medicine Tiffany J. Callahan Computational Bioscience Program Hunter/Kahn Labs Network Analysis Working Group 09.28.2017 Network-Enabled Wisdom (NEW) empirically
More informationAccepted Manuscript. Boolean Modeling of Biological Regulatory Networks: A Methodology Tutorial. Assieh Saadatpour, Réka Albert
Accepted Manuscript Boolean Modeling of Biological Regulatory Networks: A Methodology Tutorial Assieh Saadatpour, Réka Albert PII: S1046-2023(12)00277-0 DOI: http://dx.doi.org/10.1016/j.ymeth.2012.10.012
More informationLecture 10: May 27, 2004
Analysis of Gene Expression Data Spring Semester, 2004 Lecture 10: May 27, 2004 Lecturer: Ron Shamir Scribe: Omer Czerniak and Alon Shalita 1 10.1 Genetic Networks 10.1.1 Preface An ultimate goal of a
More informationModels of transcriptional regulation
Models of transcriptional regulation We have already discussed four simple mechanisms of transcriptional regulation, nuclear exclusion nuclear concentration modification of bound activator redirection
More informationApproximate inference for stochastic dynamics in large biological networks
MID-TERM REVIEW Institut Henri Poincaré, Paris 23-24 January 2014 Approximate inference for stochastic dynamics in large biological networks Ludovica Bachschmid Romano Supervisor: Prof. Manfred Opper Artificial
More informationLogic-Based Modeling in Systems Biology
Logic-Based Modeling in Systems Biology Alexander Bockmayr LPNMR 09, Potsdam, 16 September 2009 DFG Research Center Matheon Mathematics for key technologies Outline A.Bockmayr, FU Berlin/Matheon 2 I. Systems
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 informationLecture 9: June 21, 2007
Analysis of Gene Expression Data Spring Semester, 2007 Lecture 9: June 21, 2007 Lecturer: Ron Shamir Scribe: Maria Natanzon and Yevgenia Koren 1 9.1 Genetic networks 9.1.1 Preface An ultimate goal of a
More informationBSc MATHEMATICAL SCIENCE
Overview College of Science Modules Electives May 2018 (2) BSc MATHEMATICAL SCIENCE BSc Mathematical Science Degree 2018 1 College of Science, NUI Galway Fullscreen Next page Overview [60 Credits] [60
More informationSYSTEMS MEDICINE: AN INTEGRATED APPROACH WITH DECISION MAKING PERSPECTIVE. A Dissertation BABAK FARYABI
SYSTEMS MEDICINE: AN INTEGRATED APPROACH WITH DECISION MAKING PERSPECTIVE A Dissertation by BABAK FARYABI Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationDiscovering molecular pathways from protein interaction and ge
Discovering molecular pathways from protein interaction and gene expression data 9-4-2008 Aim To have a mechanism for inferring pathways from gene expression and protein interaction data. Motivation Why
More informationArtificial Neural Networks Examination, March 2004
Artificial Neural Networks Examination, March 2004 Instructions There are SIXTY questions (worth up to 60 marks). The exam mark (maximum 60) will be added to the mark obtained in the laborations (maximum
More informationBayesian Learning. Two Roles for Bayesian Methods. Bayes Theorem. Choosing Hypotheses
Bayesian Learning Two Roles for Bayesian Methods Probabilistic approach to inference. Quantities of interest are governed by prob. dist. and optimal decisions can be made by reasoning about these prob.
More informationhsnim: Hyper Scalable Network Inference Machine for Scale-Free Protein-Protein Interaction Networks Inference
CS 229 Project Report (TR# MSB2010) Submitted 12/10/2010 hsnim: Hyper Scalable Network Inference Machine for Scale-Free Protein-Protein Interaction Networks Inference Muhammad Shoaib Sehgal Computer Science
More informationCS 188: Artificial Intelligence Fall 2008
CS 188: Artificial Intelligence Fall 2008 Lecture 14: Bayes Nets 10/14/2008 Dan Klein UC Berkeley 1 1 Announcements Midterm 10/21! One page note sheet Review sessions Friday and Sunday (similar) OHs on
More informationPackage ENA. February 15, 2013
Package ENA February 15, 2013 Type Package Title Ensemble Network Aggregation Version 1.2-4 Date 2013-02-14 Author Maintainer Depends R (>= 2.13.0), space (>= 0.1), WGCNA (>= 1.20), GeneNet(>= 1.2.5),
More informationUnsupervised machine learning
Chapter 9 Unsupervised machine learning Unsupervised machine learning (a.k.a. cluster analysis) is a set of methods to assign objects into clusters under a predefined distance measure when class labels
More informationLearning and Memory in Neural Networks
Learning and Memory in Neural Networks Guy Billings, Neuroinformatics Doctoral Training Centre, The School of Informatics, The University of Edinburgh, UK. Neural networks consist of computational units
More informationNoisy Attractors and Ergodic Sets in Models. of Genetic Regulatory Networks
Noisy Attractors and Ergodic Sets in Models of Genetic Regulatory Networks Andre S. Ribeiro Institute for Biocomplexity and Informatics, Univ. of Calgary, Canada Department of Physics and Astronomy, Univ.
More informationLecture 8: Temporal programs and the global structure of transcription networks. Chap 5 of Alon. 5.1 Introduction
Lecture 8: Temporal programs and the global structure of transcription networks Chap 5 of Alon 5. Introduction We will see in this chapter that sensory transcription networks are largely made of just four
More informationUsing a Hopfield Network: A Nuts and Bolts Approach
Using a Hopfield Network: A Nuts and Bolts Approach November 4, 2013 Gershon Wolfe, Ph.D. Hopfield Model as Applied to Classification Hopfield network Training the network Updating nodes Sequencing of
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 informationBayesian Networks: Construction, Inference, Learning and Causal Interpretation. Volker Tresp Summer 2014
Bayesian Networks: Construction, Inference, Learning and Causal Interpretation Volker Tresp Summer 2014 1 Introduction So far we were mostly concerned with supervised learning: we predicted one or several
More informationMachine Learning. Neural Networks
Machine Learning Neural Networks Bryan Pardo, Northwestern University, Machine Learning EECS 349 Fall 2007 Biological Analogy Bryan Pardo, Northwestern University, Machine Learning EECS 349 Fall 2007 THE
More informationAnalog Electronics Mimic Genetic Biochemical Reactions in Living Cells
Analog Electronics Mimic Genetic Biochemical Reactions in Living Cells Dr. Ramez Daniel Laboratory of Synthetic Biology & Bioelectronics (LSB 2 ) Biomedical Engineering, Technion May 9, 2016 Cytomorphic
More informationReification of Boolean Logic
526 U1180 neural networks 1 Chapter 1 Reification of Boolean Logic The modern era of neural networks began with the pioneer work of McCulloch and Pitts (1943). McCulloch was a psychiatrist and neuroanatomist;
More information10-810: Advanced Algorithms and Models for Computational Biology. Optimal leaf ordering and classification
10-810: Advanced Algorithms and Models for Computational Biology Optimal leaf ordering and classification Hierarchical clustering As we mentioned, its one of the most popular methods for clustering gene
More information56:198:582 Biological Networks Lecture 10
56:198:582 Biological Networks Lecture 10 Temporal Programs and the Global Structure The single-input module (SIM) network motif The network motifs we have studied so far all had a defined number of nodes.
More informationArtificial Neural Network
Artificial Neural Network Contents 2 What is ANN? Biological Neuron Structure of Neuron Types of Neuron Models of Neuron Analogy with human NN Perceptron OCR Multilayer Neural Network Back propagation
More informationOverview of Research at Bioinformatics Lab
Overview of Research at Bioinformatics Lab Li Liao Develop new algorithms and (statistical) learning methods that help solve biological problems > Capable of incorporating domain knowledge > Effective,
More informationAsynchronous Stochastic Boolean Networks as Gene Network Models
Journal of Computational Biology Journal of Computational Biology: http://mc.manuscriptcentral.com/liebert/jcb Asynchronous Stochastic Boolean Networks as Gene Network Models Journal: Journal of Computational
More informationIllustration of the K2 Algorithm for Learning Bayes Net Structures
Illustration of the K2 Algorithm for Learning Bayes Net Structures Prof. Carolina Ruiz Department of Computer Science, WPI ruiz@cs.wpi.edu http://www.cs.wpi.edu/ ruiz The purpose of this handout is to
More informationarxiv: v1 [cs.sy] 25 Oct 2017
Reconstruct the Logical Network from the Transition Matrix Cailu Wang, Yuegang Tao School of Control Science and Engineering, Hebei University of Technology, Tianjin, 300130, P. R. China arxiv:1710.09681v1
More informationACTA PHYSICA DEBRECINA XLVI, 47 (2012) MODELLING GENE REGULATION WITH BOOLEAN NETWORKS. Abstract
ACTA PHYSICA DEBRECINA XLVI, 47 (2012) MODELLING GENE REGULATION WITH BOOLEAN NETWORKS E. Fenyvesi 1, G. Palla 2 1 University of Debrecen, Department of Experimental Physics, 4032 Debrecen, Egyetem 1,
More informationAN INTRODUCTION TO NEURAL NETWORKS. Scott Kuindersma November 12, 2009
AN INTRODUCTION TO NEURAL NETWORKS Scott Kuindersma November 12, 2009 SUPERVISED LEARNING We are given some training data: We must learn a function If y is discrete, we call it classification If it is
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 informationMarkov Networks.
Markov Networks www.biostat.wisc.edu/~dpage/cs760/ Goals for the lecture you should understand the following concepts Markov network syntax Markov network semantics Potential functions Partition function
More informationComputational methods for predicting protein-protein interactions
Computational methods for predicting protein-protein interactions Tomi Peltola T-61.6070 Special course in bioinformatics I 3.4.2008 Outline Biological background Protein-protein interactions Computational
More informationDifferential Modeling for Cancer Microarray Data
Differential Modeling for Cancer Microarray Data Omar Odibat Department of Computer Science Feb, 01, 2011 1 Outline Introduction Cancer Microarray data Problem Definition Differential analysis Existing
More informationBayesian Networks: Construction, Inference, Learning and Causal Interpretation. Volker Tresp Summer 2016
Bayesian Networks: Construction, Inference, Learning and Causal Interpretation Volker Tresp Summer 2016 1 Introduction So far we were mostly concerned with supervised learning: we predicted one or several
More informationPredicting Protein Functions and Domain Interactions from Protein Interactions
Predicting Protein Functions and Domain Interactions from Protein Interactions Fengzhu Sun, PhD Center for Computational and Experimental Genomics University of Southern California Outline High-throughput
More informationIntroduction to Statistical Inference
Structural Health Monitoring Using Statistical Pattern Recognition Introduction to Statistical Inference Presented by Charles R. Farrar, Ph.D., P.E. Outline Introduce statistical decision making for Structural
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 informationLecture 10: Cyclins, cyclin kinases and cell division
Chem*3560 Lecture 10: Cyclins, cyclin kinases and cell division The eukaryotic cell cycle Actively growing mammalian cells divide roughly every 24 hours, and follow a precise sequence of events know as
More informationPhiladelphia University Faculty of Engineering
Philadelphia University Faculty of Engineering Marking Scheme Exam Paper BSc CE Logic Circuits (630211) Final Exam First semester ate: 03/02/2019 Section 1 Weighting 40% of the module total Lecturer: Coordinator:
More informationANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY. Nael H. El-Farra, Adiwinata Gani & Panagiotis D. Christofides
ANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY Nael H El-Farra, Adiwinata Gani & Panagiotis D Christofides Department of Chemical Engineering University of California, Los Angeles 2003 AIChE
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 informationLogic Regression: Biological Motivation Cyclic Gene Study
Logic Regression Earl F. Glynn March 005 1 Logic Regression 1. Biological Motivation. Background Regression Boolean Algebra 3. Logic Regression Formalities 4. Simple Example in R 5. Take Home Message 1
More informationANAXOMICS METHODOLOGIES - UNDERSTANDING
ANAXOMICS METHODOLOGIES - UNDERSTANDING THE COMPLEXITY OF BIOLOGICAL PROCESSES Raquel Valls, Albert Pujol ǂ, Judith Farrés, Laura Artigas and José Manuel Mas Anaxomics Biotech, c/balmes 89, 08008 Barcelona,
More informationSupplementary Materials
Electronic Supplementary Material (ESI) for Integrative Biology. This journal is The Royal Society of Chemistry 2015 Predicting genetic interactions from Boolean models of biological networks Supplementary
More informationMachine Learning for Data Science (CS4786) Lecture 24
Machine Learning for Data Science (CS4786) Lecture 24 Graphical Models: Approximate Inference Course Webpage : http://www.cs.cornell.edu/courses/cs4786/2016sp/ BELIEF PROPAGATION OR MESSAGE PASSING Each
More informationSupplementary methods
Supplementary methods The rxncon language The network reconstruction was performed using the reaction-contingency (rxncon) language (1). The language is based on the separation of two distinct classes
More informationEfficient Inference in Fully Connected CRFs with Gaussian Edge Potentials
Efficient Inference in Fully Connected CRFs with Gaussian Edge Potentials by Phillip Krahenbuhl and Vladlen Koltun Presented by Adam Stambler Multi-class image segmentation Assign a class label to each
More informationV12 Gene Regulatory Networks, Boolean Networks
V12 Gene Regulatory Networks, Boolean Networks Tue, May 28, 2018 1 Gene Expression Sequence of processes: from DNA to functional proteins transcription transcribed DNA mrna mrna RNA TFs regulation at every
More informationOptimal State Estimation for Boolean Dynamical Systems using a Boolean Kalman Smoother
Optimal State Estimation for Boolean Dynamical Systems using a Boolean Kalman Smoother Mahdi Imani and Ulisses Braga-Neto Department of Electrical and Computer Engineering Texas A&M University College
More information