The topology and signature of the regulatory interactions predict the expression pattern of the segment polarity genes in Drosophila m elanogaster

Similar documents
In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

On One Analytic Method of. Constructing Program Controls

Chapter 6: AC Circuits

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Let s treat the problem of the response of a system to an applied external force. Again,

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

2/20/2013. EE 101 Midterm 2 Review

Linear Response Theory: The connection between QFT and experiments

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

FTCS Solution to the Heat Equation

Transcription: Messenger RNA, mrna, is produced and transported to Ribosomes

TSS = SST + SSE An orthogonal partition of the total SS

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys

( ) () we define the interaction representation by the unitary transformation () = ()

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

Cubic Bezier Homotopy Function for Solving Exponential Equations

WiH Wei He

P R = P 0. The system is shown on the next figure:

Relative controllability of nonlinear systems with delays in control

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

National Exams December 2015 NOTES: 04-BS-13, Biology. 3 hours duration

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

Graduate Macroeconomics 2 Problem set 5. - Solutions

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

Political Economy of Institutions and Development: Problem Set 2 Due Date: Thursday, March 15, 2019.

Multi-Fuel and Mixed-Mode IC Engine Combustion Simulation with a Detailed Chemistry Based Progress Variable Library Approach

Robust and Accurate Cancer Classification with Gene Expression Profiling

Robustness Experiments with Two Variance Components

Variants of Pegasos. December 11, 2009

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Chapter 4. Neural Networks Based on Competition

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

Comb Filters. Comb Filters

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Density Matrix Description of NMR BCMB/CHEM 8190

Machine Learning Linear Regression

Department of Economics University of Toronto

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment

On the Boyd- Kuramoto Model : Emergence in a Mathematical Model for Adversarial C2 Systems

Density Matrix Description of NMR BCMB/CHEM 8190

Implementation of Quantized State Systems in MATLAB/Simulink

Switching Regulatory Models of Cellular Stress Response

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)

Solution in semi infinite diffusion couples (error function analysis)

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

Analysis And Evaluation of Econometric Time Series Models: Dynamic Transfer Function Approach

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

Motivations. Opinions, influence networks and centrality. Small deliberative groups. Opinion Dynamics and the Evolution of Influence Networks

Notes on the stability of dynamic systems and the use of Eigen Values.

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

Fall 2010 Graduate Course on Dynamic Learning

Time-interval analysis of β decay. V. Horvat and J. C. Hardy

Lecture 6: Learning for Control (Generalised Linear Regression)

( ) [ ] MAP Decision Rule

5. Stochastic processes (1)

Epistemic Game Theory: Online Appendix

Solutions Problem Set 3 Macro II (14.452)

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

First-order piecewise-linear dynamic circuits

Lecture VI Regression

Multi-Objective Control and Clustering Synchronization in Chaotic Connected Complex Networks*

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

Mean field approximation for PDE-Markov random field models in image analysis

Standard Error of Technical Cost Incorporating Parameter Uncertainty

Sampling Procedure of the Sum of two Binary Markov Process Realizations

Lecture 9: Dynamic Properties

Sliding Mode Controller for Unstable Systems

Chapter 2 Linear dynamic analysis of a structural system

Boolean Modeling of Genetic Regulatory Networks

Stochastic Repair and Replacement with a single repair channel

Learning Objectives. Self Organization Map. Hamming Distance(1/5) Introduction. Hamming Distance(3/5) Hamming Distance(2/5) 15/04/2015

II. Light is a Ray (Geometrical Optics)

Fundamentals of PLLs (I)

Reactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times

Supporting Information: The integrated Global Temperature change Potential (igtp) and relationships between emission metrics

Building Temporal Models for Gesture Recognition

EP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES

Motion in Two Dimensions

Use of Kalman Filtering and Particle Filtering in a Benzene Leachate Transport Model

Hidden Markov Models Following a lecture by Andrew W. Moore Carnegie Mellon University

Born Oppenheimer Approximation and Beyond

Computing Relevance, Similarity: The Vector Space Model

Volatility Interpolation

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

Transient Response in Electric Circuits

Estimation of Cost and. Albert Banal-Estanol

Scattering at an Interface: Oblique Incidence

Dynamic Team Decision Theory

Example: MOSFET Amplifier Distortion

OMXS30 Balance 20% Index Rules

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015

Forecasting customer behaviour in a multi-service financial organisation: a profitability perspective

Transcription:

The opology and sgnaure of he regulaory neracons predc he expresson paern of he segmen polary genes n Drosophla m elanogaser Hans Ohmer and Réka Alber Deparmen of Mahemacs Unversy of Mnnesoa

Complex bologcal neworks Bologcal sysems form neworks of neracon gene conrol neworks meabolc neworks neworks of sgnal ransducon and response organ and organsm-level neworks populaon level neworks The funcon of hese sysems canno be undersood n a reducons way or by assumng lnear pahways. I s mporan o uncover he opology of he neracon neworks. sascal descrpon of nework opology sudy he nerplay beween opology and funcon n specfc neworks

Wha s Drosophla melanogaser?

Drosophla segmenaon s governed by a cascade of genes 2:50 h 1:20 h 2:10 h Segmen polary genes: wngless (wg) engraled (en) hedgehog (hh) pached (pc) smoohened (smo) cubus nerrupus (c) sloppy pared (slp)

Goal: relae he paern of he segmen polary genes o her neracons 1. The segmen polary genes nerac va a complex regulaory nework. 2. I s possble o model hs nework based only on he opology and sgnaure of he neracons. 3. The resuls of he model are n good agreemen wh he observed wld ype and muan paerns. 4. The model reveals a remarkable robusness wh respec o perurbaons n he nal condons.

The segmen polary genes are naed by he par-rule genes en hh wg early sages 2:50 h pre-paern 3:00-3:30 h sable paern 4:20-7:20 h 3:30 h c en pc

Sable gene paerns en s expressed n he aneror par of he parasegmen. wg s expressed n he poseror par of he parasegmen. The parasegmenal grooves form beween he wg and en srpes. The pc srpes separae no wo. The c paern s complemenary o ha of en. c en

The expresson of he segmen polary genes s mananed by a regulaory nework Cross-regulaon and feedback loops make hese neworks complex.

Inercellular neracons play a val role Cadgan, Nusse, Genes & Developmen 11, 3286 (1997)

Reconsrucng he opology of he segmen polary nework mrna PROTEIN PROT COMPL represson ranslaon, acvaon, modfcaon cell neghbor cell

The von Dassow model of he segmen polary genes d[ pro] = T d max ν [ mrna] ν Κ + [ mrna] ν [ pro] τ The opology of he nework s somewha dfferen from our reconsrucon. The model nvolves 48 unknown parameers. The soluons leadng o wld ype seady saes are dsrbued homogenously n he bologcally relevan parameer space. von Dassow e al., Naure 406, 188 (2000)

Does he opology of he segmen polary nework deermne s funcon? Boolean model of he nework dynamcs: Transcrps and proens are eher ON (1) or OFF(0). The expresson of a node a mesep s gven by a logcal rule of he expresson of s effecors a me -1. Transcrpon depends on ranscrpon facors; repressors are domnan. Translaon depends on he presence of he ranscrp. Transcrps and proens decay n one sep f no produced.

Rules for ranscrpon and ranslaon en hh pc c + 1 = +1 = +1 = +1 = WG WG ( 1 or + 1 EN CIA no EN = en + 1 WG = EN + 1 CI = c + 1 HH = wg hh + 1 and no and no ) and no CIR EN SLP and no CIR

Rules for pos-ranslaonal processes PH = PTC and ( HH 1 or HH + 1) SMO = no PTC or HH or HH 1 1 + nsananeous 1 CIA CI and ( SMO or hh 1or hh 1) + = CIR + 1 = CI and no SMO and no hh and no hh 1 +ε SMO + 1 +

The expresson of wg, PTC and SLP s auo-regulaed + 1 wg = ( CIA [ wg and and SLP ( CIA and no or SLP CIR ) or 1 PTC + = pc or ( PTC and no HH 1 and no HH + 1) ) and no CIR Eher of he acvaors can couner mrna decay. Free PTC does no decay. ] SLP = SLP + 1 SLP s a source n he segmen polary nework.

Sar he model from an nal sae gvng he expresson of all nodes SLP parasegmen cell wh node ON cell wh node OFF The paerns are assumed o be dencal n each parasegmen. Wld ype nal sae: wg n he las cell of he parasegmen, en/hh n he frs cell of he parasegmen, pc and c complemenary o en

The model reproduces he wld ype seady sae Agrees wh he wld ype paern of he segmen polary genes. nal sae seady sae The ne effec of he neracons s enough o capure he funconng of he nework!

The naure of he neracons can be negraed no he opology The fuure expresson of a node depends on a combnaon of he expresson of oher nodes. hh Inroduce complemenary nodes. Assocae pseudo-nodes o node combnaons. The fuure expresson of nodes depends on he expresson of pseudo-nodes. +1 = EN and no CIR no CIR ( E CR) EN hh +1 = ( ECR) CIR and CIR

The funconal opology reveals acvang srucures me-sep complemen Cycles: wg ( wcacr 2 ) 2 wg ( wscr 2 ) 2 PTC P 2 2H1H 3 Clusers deermne segmen polary en, EN, hh, HH c, CI, CIA, wg, WG, pc, PTC

Wha happens f he nework s perurbed? The mos severe perurbaon s caused by gene muaons. To model a null muaon, we assume ha he mrna s kep OFF, hus he proen canno be ranslaed. The effecs of he muaon propagae hroughou he nework.

wg, en or hh muaons are lehal fnal sae No wg, en and hh srpes, no segmenaon, regardless of nal sae.

pc muaon broadens he srpes fnal sae The wg, en and hh srpes broaden, regardless of nal sae.

c muaon can preserve he prepaern fnal sae for wld ype prepaern The effec of c muaon depends on he nal sae. For wld ype prepaern, he wg, en, hh srpes reman.

Comparson beween hh and c muans predcon for pc pc wg predcon for wg wld ype c muan hh muan Galle e al., Developmen 127, 5509 (2000)

Comparson beween c and pc muans wld ype c muan predcon, c muan pc muan predcon, pc muan en hh Tabaa, Eaon, Kornberg, Genes & Developmen 6, 2635 (1992)

Wha happens f he nal sae s perurbed? Possble number of prepaerns for a sngle node: 16 and complemens Toal number of nework-wde prepaerns: All red nal saes lead o seady saes. N 13 = 16

How many nal saes lead o he wld ype seady sae? 3 3 4 4 4 4 8 8 16 16 16 8 8 16 mnmal prepaern seady sae Toal number of wld-ype nducng prepaerns: 2.4 10 9 = 10 5 N

There are wo oher frequenly occurrng seady saes Broader naon of wg, en or hh leads o broad srpes. Absence of wg leads o a sae wh no segmenaon. 0.98N Galle e al., Developmen 127, 5509 (2000) 10 2 N

The seady saes can be deermned analycally In he sable sae x +1 = x Use he fac ha SLP = SLP2 = and SLP3 = SLP The se of equaons reduces o: 1 = wg1 and no wg2 and no wg wg 2 = wg2 and no wg1 and no wg 3 = wg1 or wg3 wg 4 = wg2 or wg4 wg wg PTC 1 = ( no wg2 and no wg4) or ( PTC1 and no wg1 and no wg3) PTC 2 = ( no wg1 and no wg3) or ( PTC2 and no wg2 and no wg4) PTC = PTC 1 3 4 = 1 0 4 = 4 3 1 Reflecs he assumpon of sably of wg and PTC.

Possble sable paerns wld ype broad lehal double wg dsplaced dsplaced, 2 wg The laer saes have very small probably.

Can he assumpons of he model be changed? 1. Node expresson decays n one sep f no renewed. 1 Assumpon: proens decay n wo seps; EN + = en Same seady saes, only he pah leadng o hem changes. More overall sably. or en 1 2. wg and PTC expresson s easer o manan han ohers. Assumpon: wg and PTC decay f no renewed. c muans have no segmenaon f wg decays; no wld ype seady sae f PTC decays. The model can be modfed o nclude dfferen mescales for ranscrpon, ranslaon, pos-ranscrponal modfcaons, mrna and proen decay.

Can we relax assumpons of mendependence? 3. SLP expresson s a necessary npu. No wld ype seady sae n a SLP muan. 4. The parasegmens reman 4 cells wde. Cell dvson can be ncorporaed no he model.

Conclusons and oulook 1. Boolean modelng confrms ha he funconng of he segmen polary nework s deermned by s opology. 2. The wld ype seady sae s robus o perurbaons n he nal condons. 3. The nework favors seady saes wh srpes even for ecopc prepaerns. 4. Boolean modelng allows for a sysemac deermnaon of seady saes and he nal saes leadng o hem. can be appled o a varey of neworks provdes a verfcaon of he compleeness of he opologcal nformaon