Interprtation abstraite de modles de voies de signalisation - - PowerPoint PPT Presentation

Interprétation abstraite de modèles de voies de signalisation intracellulaire

Jérôme Feret

Joint-work with...

Signalling Pathways

Pathway maps

Differential models

                  

A gap between two worlds

• 1. Databases of proteins interactions in natural language

• 2. ODE-based models

Rule-based approach

• 1. The description level matches with both
• the observation level
• and the intervention level
• f the biologist.

• 2. Model description is very compact.

Formal semantics

Complexity walls

Static analysis of reachable species (I/II)

Static analysis of reachable species (II/II)

• 1. check the consistency of a model [ICCMSE’07]
• 2. compute the properties to allow fast simulation [APLAS’07]
• 3. simplify models,
• 4. compute independent fragments of chemical species [PNAS’09, LICS’10,Chaos’10]

Model reduction

• ments. [PNAS’09,LICS’10,Chaos’10]

Reachability Analysis

• f Rule-based Models

Jérôme Feret

Overview

• 1. Introduction
• 2. Kappa language
• 3. Local views
• 4. Local set of chemical species
• 5. Local rule systems
• 6. Conclusion

A single story

A concurrent story

Overshoot

Overview

• 1. Introduction
• 2. Kappa language
• 3. Local views
• 4. Local set of chemical species
• 5. Local rule systems
• 6. Conclusion

A chemical species

E R R E

E(r!1), R(l!1,r!2), R(r!2,l!3), E(r!3)

A Unbinding/Binding Rule

E R E R

E(r), R(l,r) ← → E(r!1), R(l!1,r)

Internal state

E R E R

R(Y1∼u,l!1), E(r!1) ← → R(Y1∼p,l!1), E(r!1)

Don’t care, Don’t write

R

R

R

R

A contextual rule

R

R

R(Y1∼u,r!_) → R(Y1∼p,r)

Creation/Suppression

R R R

R(r) → R(r!1), R(r!1,l,Y1)

R R R

R(r!1), R(r!1) → R(r)

Early EGF example

egf rules 1

Early EGF example

egf rules 2

Properties of interest

• 1. Show the absence of modelling errors:
• detect dead rules;
• detect overlapping rules;
• detect non exhaustive interactions;
• detect rules with ambiguous molecularity.
• 2. Get idiomatic description of the networks:
• capture causality;
• capture potential interactions;
• capture relationships between site states.

• 3. Allow fast simulation:
• capture accurate approximation of the wake-up relation.

Overview

• 1. Introduction
• 2. Kappa language
• 3. Local views
• 4. Local set of chemical species
• 5. Local rule systems
• 6. Conclusion

Embedding

• Φ is a site-graph morphism:
• i is less specific than Φ(i),
• if there is a link between (i, s) and (i′, s′),

• Φ is an into map (injective):
• Φ(i) = Φ(i′) implies that i = i′.

Set of reachable chemical species

• 1. there is an embedding of the lhs of Rk in the solution c1, . . . , cm;
• 2. the (embedding/rule) produces the solution c′

Inductive definition

• c′

• ∃Rk ∈ R, c1, . . . , cm ∈ X,

• .

• n ∈ N
• .

Local views

E R E R R E

α({R(Y1∼u,l!1), E(r!1)}) = {R(Y1∼u,l!r.E); E(r!l.R)}.

Galois connexion

• Let Local_view be the set of all local views.
• Let α ∈ ℘(Species) → ℘(Local_view) be the function that maps any set
• f chemical species into the set of their local views.
• Let γ ∈ ℘(Local_view) → ℘(Species) be the function that maps any

• The pair (α, γ) forms a Galois connexion:

γ ◦ α

{ }

R

α ◦ γ

{ }

R

;

R R R R

One more question

{ }

R

;

R R R

Abstract rules

#

R R R R E R R R E R R

R

Abstract counterpart to F

• lv′

• ∃Rk ∈ R, lv1, . . . , lvm ∈ X,

• .

• F♯ is monotonic;
• F ◦ γ

Soundness

• 1. (D, ⊆, ∪) and (D♯, ⊑, ⊔) be chain-complete partial orders;
• 2. D −

• 3. F ∈ D → D and F♯ ∈ D♯ → D♯ be monotonic mappings such that:

• 4. x0 ∈ D be an element such that: x0 ⊆ F(x0);

• 1. both lfpx0F and lfpα(x0)F♯ exist,
• 2. lfpx0F ⊆ γ(lfpα(x0)F♯).

Overview

• 1. Introduction
• 2. Kappa language
• 3. Local views
• 4. Local set of chemical species
• 5. Local rule systems
• 6. Conclusion

Concretization

I

F E

E

Pumping lemma

• We use this rewrite system to enumerate the chemical species of γ(X).
• There are two cases:
• 1. either there is a finite number of rewrite sequences;
• 2. or we encounter cyclic derivations

R.l-r.E ... R.l-r.E

• We only enumerate chemical species that are reached through an acyclic

• It turns out that: if X ∈ Im(α) then each rewrite sequence is the prefix of

Examples

• 1. Make the demo for egf
• 2. Make the demo for fgf
• 3. Make the demo for sfb
• 4. Make the demo for Global invariants

Which information is abstracted away ?

• for EGF cascade (356 chemical species);
• for FGF cascade (79080 chemical species);
• for SFB cascade (≈ 2.1019 chemical species);

• Can we characterize which information is abstracted away ?
• Which is the form of the systems, for which we have no false positive ?
• Do we learn something about the biological systems that we describe ?

Which information is abstracted away ?

• 1. (D, ⊆) be a partial order;
• 2. (D♯, ⊑, ⊔) be chain-complete partial order;
• 3. D −

• 4. F ∈ D → D and F♯ ∈ D♯ → D♯ are monotonic;
• 5. F ◦ γ

• 6. x0, inv ∈ D such that:
• x0 ⊆ F(x0) ⊆ F(inv) ⊆ inv,
• inv = γ(α(inv)),
• and α(F(γ(α(inv)))) = F♯(α(inv));

When is there no false positive ?

• 1. (D, ⊆, ∪) and (D♯, ⊑, ⊔) are chain-complete partial orders;
• 2. (D, ⊆) −

• 3. F : D → D is a monotonic map;
• 4. x0 is a concrete element such that x0 ⊆ F(x0);
• 5. F ◦ γ

• 6. F♯ ◦ α = α ◦ F ◦ γ ◦ α.

• lfpx0F and lfpα(x0)F♯ exist;
• lfpx0F ∈ Im(γ) ⇐

Local set of chemical species

Swapping relation

Swapping closure

• 1. the set X ⊆ Species is local;
• 2. for any tuples (Ci)i∈I, (Dj)j∈J ∈ Species∗ such that:
• (Ci)i∈I ∈ X∗,
• and (Ci)i∈I

Proof (easier implication way)

• X = γ(α(X)),
• (Ci)i∈I ∈ X∗,
• and (Ci)i∈I

Proof (more difficult implication way)

• By def. of α, this is satisfied for any local view in α(X);
• This remains satisfied after unfolding a semi-link with a local view;
• This remains satisfied after binding two semi-links.

Initialization

E F E

E F E

I

Unfolding a semi-link

E F F F E

F E

• pen partial species

Unfolding a semi-link

F

F E

F E

• pen partial species

• SWAP

• Jérôme Feret

Binding two semi-links

F E F E

• pen partial species
• pen partial species

• SWAP

• Jérôme Feret

Bonus

• 1. Each open complex C built with the local views in X is a sub-complex of

• 2. When considering the rewrite system that computes γ(X), any partial

• 3. We have F♯ ◦ α = α ◦ F ◦ γ ◦ α.

Overview

• 1. Introduction
• 2. Kappa language
• 3. Local views
• 4. Local set of chemical species
• 5. Local rule systems
• 6. Conclusion

Outline

• if the set Speciesω of reachable chemical species is close with respect

• then the reachability analysis is exact (i.e. Speciesω = γ(lfpα(Species0)F♯)).

Sufficient conditions

• 1. initial agents are not bound;
• 2. rules are atomic;
• 3. rules are local:
• only agents that interact are tested,
• no cyclic patterns (neither in lhs, nor in rhs);
• 4. binding rules do not interfere i.e. if both:
• A(a∼m,S),B(b∼n,T) → A(a∼m!1,S),B(b∼n!1,T)
• and A(a∼m’,S’),B(b∼n’,T’) → A(a∼m’!1,S’),B(b∼n’!1,T’),

• A(a∼m,S),B(b∼n’,T’) → A(a∼m!1,S),B(b∼n’!1,T’);
• 5. chemical species in γ(α(Speciesω)) are acyclic,

Proof outline

• We consider tuples of complexes in which the same kind of links occur

• We want to swap these links.
• We introduce the history of their computation.
• There are several cases. . .

First case (I/V)

First case (II/V)

First case (III/V)

First case (IV/V)

• are always possible;
• have the same effect on local views;
• commute with the swapping relation

Cn

∼ ,σ

• R,φ
• C′

n R,φ

• Cn+1

∼ ,σ

n+1

First case (V/V)

Second case (I/II)

Second case (II/II)

Sufficient conditions

• 1. initial agents are not bound;
• 2. rules are atomic;
• 3. rules are local:
• only agents that interact are tested,
• no cyclic patterns (neither in lhs, nor in rhs);
• 4. binding rules do not interfere i.e. if both:
• A(a∼m,S),B(b∼n,T) → A(a∼m!1,S),B(b∼n!1,T)
• and A(a∼m’,S’),B(b∼n’,T’) → A(a∼m’!1,S’),B(b∼n’!1,T’),

• A(a∼m,S),B(b∼n’,T’) → A(a∼m!1,S),B(b∼n’!1,T’);
• 5. chemical species in γ(α(Speciesω)) are acyclic,

Third case (I/III)

Third case (II/III)

Third case (III/III)

Non local systems

Species0

∆

= R(a∼u) Rules

∆

=          R(a∼u) ↔ R(a∼p) R(a∼u),R(a∼u) → R(a∼u!1),R(a∼u!1) R(a∼p),R(a∼u) → R(a∼p!1),R(a∼p!1) R(a∼p),R(a∼p) → R(a∼p!1),R(a∼p!1)          R(a∼u!1),R(a∼u!1) ∈ Speciesω R(a∼p!1),R(a∼p!1) ∈ Speciesω But R(a∼u!1),R(a∼p!1) ∈ Speciesω.

Non local systems

Species0

∆

= A(a∼u),B(a∼u) Rules

∆

=    A(a∼u),B(a∼u) → A(a∼u!1),B(a∼u!1) A(a∼u!1),B(a∼u!1) → A(a∼p!1),B(a∼u!1) A(a∼u!1),B(a∼u!1) → A(a∼u!1),B(a∼p!1)    A(a∼u!1),B(a∼p!1) ∈ Speciesω A(a∼p!1),B(a∼u!1) ∈ Speciesω But A(a∼p!1),B(a∼p!1) ∈ Speciesω.

Non local systems

Species0

∆

= A(a∼u) Rules

∆

= A(a∼u) ↔ A(a∼p) A(a∼u),A(a∼p) → A(a∼u!1),A(a∼p!1)

• A(a∼u!1),A(a∼p!1) ∈ Speciesω

But A(a∼p!1),A(a∼p!1) ∈ Speciesω.

Non local systems

Species0

∆

= R(a,b) Rules

∆

= { R(a,b),R(a) → R(a,b!1),R(a!1)} R(a,b!2),R(a!2,b!1),R(a!1,b)∈ Speciesω But R(a!1,b!1) ∈ Speciesω.

Overview

• 1. Introduction
• 2. Kappa language
• 3. Local views
• 4. Local set of chemical species
• 5. Local rule systems
• 6. Conclusion

Conclusion

• A scalable static analysis to abstract the reachable chemical species.
• A class of models for which the abstraction is complete.
• Many applications:
• idiomatic description of reachable chemical species;
• dead rule detection;
• rule decontextualization;
• computer-driven kinetic refinement.
• It can also help simulation algorithms:
• wake up/inhibition map (agent-based simulation);
• flat rule system generation (for bounded set of chemical species);
• on the fly flat rule generation (for large/unbounded set)

Formal model reduction

Jérôme Feret

Overview

• 1. Context and motivations
• 2. Handmade ODEs
• 3. Abstract interpretation framework
• 4. Concrete semantics
• 5. Abstract semantics
• 6. Conclusion

Rule-based models

Complexity walls

A breach in the wall(s) ?

Overview

• 1. Context and motivations
• 2. Handmade ODEs

• 3. Abstract interpretation framework
• 4. Concrete semantics
• 5. Abstract semantics
• 6. Conclusion

A simple adapter

A C B

A simple adapter

A C B

A simple adapter

A C B

A simple adapter

A C B

Two subsystems

A C B

Two subsystems

• C

A B B

Two subsystems

• C

A B B

Dependence index

• .

Overview

• 1. Context and motivations
• 2. Handmade ODEs

• 3. Abstract interpretation framework
• 4. Concrete semantics
• 5. Abstract semantics
• 6. Conclusion

A system with a switch

A system with a switch

A system with a switch

Two subsystems

Two subsystems

Two subsystems

Dependence index

Conclusion

Overview

• 1. Context and motivations
• 2. Handmade ODEs
• 3. Abstract interpretation framework

• 4. Concrete semantics
• 5. Abstract semantics
• 6. Conclusion

Differential semantics

• V

• F(ρ)

Overview

• 1. Context and motivations
• 2. Handmade ODEs
• 3. Abstract interpretation framework

• 4. Concrete semantics
• 5. Abstract semantics
• 6. Conclusion

Abstraction

• V♯: a finite set of observables,
• ψ: a mapping from V → R into V♯ → R,
• F♯: a C∞ mapping from V♯ → R+ into V♯ → R;

• ψ is linear with positive coefficients,
• the following diagram commutes:

 

• 

 ψ

Abstraction example

• V

• F(ρ)

• V♯ ∆

• ψ(ρ)

• F♯(ρ♯)

Abstract differential semantics

Abstract differential semantics Proof sketch

Fluid trajectories

t Y(t)

Fluid trajectories

t Y(t) X(t)

Overview

• 1. Context and motivations
• 2. Handmade ODEs
• 3. Abstract interpretation framework
• 4. Concrete semantics
• 5. Abstract semantics
• 6. Conclusion

Differential system

Overview

• 1. Context and motivations
• 2. Handmade ODEs
• 3. Abstract interpretation framework
• 4. Concrete semantics
• 5. Abstract semantics

• 6. Conclusion

Abstract domain

• restrict the exploration;
• drive the intuition (by using the flow of information);
• having efficient way to find suitable abstractions (V♯, ψ)

Contact map

G E R Sh So

Annotated contact map

G E R Sh So

Pattern annotation

G E R Sh So R R

Pattern annotation

G E R Sh So R R

Pattern annotation

G E R Sh So R R

Prefragments and fragments

R R

R R

Are they fragments ?

Are they fragments ?

G So

Are they fragments ?

G So

Are they fragments ?

G So

Are they fragments ?

So G Sh

Are they fragments ?

So G Sh

Are they fragments ?

So G Sh

Are they fragments ?

So G Sh

Are they fragments ?

So G Sh

Are they fragments ?

So G Sh

Are they fragments ?

So G Sh

Are they fragments ?

So G Sh

Are they fragments ?

So G Sh

Are they fragments ?

Are they fragments ? stage 2

R R

Are they fragments ? stage 2

R R

Are they fragments ? stage 2

R R

Are they fragments ? stage 2

R R

Are they fragments ? stage 2

R R

Orthogonal refinement

Overview

• 1. Context and motivations
• 2. Handmade ODEs
• 3. Abstract interpretation framework
• 4. Concrete semantics
• 5. Abstract semantics

• 6. Conclusion

Flow of information

Overview

• 1. Context and motivations
• 2. Handmade ODEs
• 3. Abstract interpretation framework
• 4. Concrete semantics
• 5. Abstract semantics

• 6. Conclusion

Fragments consumption Proper intersection

Sh R Sh R Sh R

u u p

Fragment consumption

Fragment production Proper inter

E R G R G R R G

Fragment production Proper intersection (bis)

E R G R G R R G E R E R

• ments. This is already satisfied thanks to the previous syntactic criteria.

Fragment production

• C′

• SYM(C1, . . . , Cm) · SYM(F);
• therwise, we get no contribution.

Fragment properties

• an annotated contact map satisfies the syntactic criteria,
• fragments are defined by this annotated contact map,
• we know the concentration of fragments;

• we can express the concentration of any connected component occur-

• we can express fragment proper consumption,
• we can express fragment proper production,
• WE HAVE A CONSTRUCTIVE DEFINITION FOR F♯.

Overview

• 1. Context and motivations
• 2. Handmade ODEs
• 3. Abstract interpretation framework
• 4. Concrete semantics
• 5. Abstract semantics
• 6. Conclusion

Experimental results

Related issues Context sensitive approximation

Related issues Approximation of the stochastic semantics

Model reduction of stochastic rules-based models

Jérôme Feret

Joint-work with...

Overview

• 1. Introduction
• 2. Examples of information flow
• 3. Conclusion

ODE fragments

Stochastic fragments ?

Overview

• 1. Introduction
• 2. Examples of information flow
• 3. Conclusion

A model with ubiquitination

?

?

Statistical independence

(n⋆P + n⋆P⋆)(nP⋆ + n⋆P⋆)

• .

• 2.5e-16
• 2e-16
• 1.5e-16
• 1e-16
• 5e-17

Reduced model

Comparison between the two models

• 5e-16
• 4e-16
• 3e-16
• 2e-16
• 1e-16

Coupled semi-reactions

Reduced model

Comparison between the two models

Degree of correlation (in the unreduced model)

Distant control

Reduced model

Comparison between the two models

• 2.5
• 2
• 1.5
• 1
• 0.5

Degree of correlation (in the unreduced model)

• nA⋆

• t

• 2
• 1.5
• 1
• 0.5

Overview

• 1. Introduction
• 2. Examples of information flow
• 3. Conclusion

Hierarchy of semantics

Conclusion

• A framework for reducing stochastic rule-based models.
• We use:

• Algebraic operators combine these abstractions.
• We use backward bisimulations in order to prove statistical invariants,