330 likes | 348 Views
Explore the RB pathway in cancer, its complexities, and key proteins involved. Use BiNoM with CellDesigner to analyze and model pathway interactions extracted from literature for deeper insights.
E N D
Disentangling RB-pathway with BiNoM Laurence Calzone and Andrei Zinovyev Service Bioinformatique, Institut Curie, Paris 26/06/2007 MOCA IHP, Paris
Project on modeling regulation • of cell cycle by • RB/E2F Pathway • Travail de collaboration entre les équipes de : • - bioinformatique d’Emmanuel Barillot • - oncologie moléculaire de François Radvanyi 26/06/2007 MOCA IHP, Paris
CKI (p16) RB-P RB in tumor progression • RB is a cell proliferation regulator: RB is a tumor suppressor that controls cell cycle progression (G0/G1 and G1/S transitions). • The importance of the RB pathway is known in cancer: • RB is targeted and inactivated by viral oncoproteins (E7) • RB pathway is inactivated in most tumoral cells in different ways: - RB gene mutation - Deregulation of the kinases that control its activity Mitogenic Signal Cyc D CDK4/6 E7 RB E2F Progression of Cell Cycle 26/06/2007 MOCA IHP, Paris
CKI (p16) RB-P Construction of RB/E2F pathway from the literature Simplified view Mitogenic Signal Cyc D CDK4/6 E7 RB E2F Progression of Cell Cycle 26/06/2007 MOCA IHP, Paris
RB/E2F pathway: 78 proteins, 169 genes, 208 species, 166 reactions, more than 340 publications… 26/06/2007 MOCA IHP, Paris
Tools CellDesigner (SBML) to describe the interactions between proteins and organize the information extracted from the literature. BiNoM (Cytoscape plug-in) to manipulate the diagrams and query the network 26/06/2007 MOCA IHP, Paris
Comparing our pathway with other pathways METHOD : A subnetwork is extracted from the different databases (including only proteins present in our pathway) REACTOME Twice as many details in our pathway : - Species : 97 in REACTOME, 130 in our diagram - Reactions : 48 in REACTOME, 118 in our diagram - Entities (proteins or group of proteins) : 47 proteins in REACTOME, 74 in our diagram 36 complexes in REACTOME, 98 in our diagram • TRANSPATH - Species : 145 in TRANSPATH, 130 in our diagram - Reactions : 70 in TRANSPATH, 118 in our diagram - Entities (proteins or group of proteins) : 53 proteins in TRANSPATH, 74 in our diagram 28 complexes in TRANSPATH, 98 in our diagram 26/06/2007 MOCA IHP, Paris
Webpage + Publication 26/06/2007 MOCA IHP, Paris
Modeling the pathway Early G1 Late G1 S G2 M G0 26/06/2007 MOCA IHP, Paris
RB/E2F Pathway: paradoxes Learn with the model of the pathway • Loss of Rb with overexpression of p16 • Cf données transcriptome/génome • In the tumors of the bladder: gain/amplification of E2F3 observed with loss of Rb • Cf données transcriptome/génome • -E2F3 is involved as an oncogene but not E2F2 or E2F1 in bladder cancers • -Tumor suppressor and oncogene role of E2F1 Pierce et al., E2F1 has both oncogenic and tumor-suppressive properties in a transgenic model.Mol Cell Biol. 1999 • Sep;19(9):6408-14. • -Oncogenic transformation by Ras requires that RB gene is functionalWilliams et al.,The retinoblastoma protein is required for Ras-induced oncogenic transformation. Mol Cell Biol. 2006 Feb;26(4):1170-82. 26/06/2007 MOCA IHP, Paris
Dividing the pathway in modules • What to do? • Define the sub-networks (by hand and if possible automatically) • Isolate the sub-networks from the rest of the pathway • Establish the links between the different sub-networks • What for? • To read, understand and exploit the pathway • To improve the modules independently from the rest of the pathway • What next? • Modify the links according to the new proteins/genes/interactions added • Refine the subnetworks and model them • The more abstract view can be modeled and already answer questions, propose some predictive phenotypes 26/06/2007 MOCA IHP, Paris
Modular decomposition of the RB/E2F pathway • 16 modules: • RB • E2F Transcription Factors • Activator E2F1 • Repressor E2F4 • Repressor E2F6 • Cyclin-dependent kinases • CycC/CDK3 • CycH/CDK7 • CycD1/CDK4,6 • CycE1/CDK2 • CycA2/CDK2 • CycB1/CDC2 • Cyclin-dependent kinase Inhibitors • p16INK4a, p15INK4b • p21CIP1, p27KIP1 • Phosphorylation of cyclin B kinase • CDC25C • WEE1 • Degradation of cyclin • APC • Apoptosis entry Early G1 Late G1 S G2 M G0 26/06/2007 MOCA IHP, Paris
Some modules: RB RB module RB 26/06/2007 MOCA IHP, Paris
Some modules: E2F1 E2F1 Module E2F1 26/06/2007 MOCA IHP, Paris
Discrete modeling of the pathway( with Denis Thieffry ) 14 nodes (without apoptosis entry module) and influences of one major protein on the others Drawn and computed later in GinSIM (Thieffry) 26/06/2007 MOCA IHP, Paris
Project . • Biological Network Manager • Cytoscape plugin • BiNoM version 1 release date – 1 September 2007The project is open-sourcehttp://bioinfo.curie.fr/projects/binom 26/06/2007 MOCA IHP, Paris
Four BiNoM modules BiNoM Structure analysis BiNoM I/O BiNoM BioPAX Query BiNoM Utilities 26/06/2007 MOCA IHP, Paris
BiNoM I/O CellDesigner CellDesigner User manipulations in BiNoM using Cytoscape graphics BioPAX BioPAX SBML 26/06/2007 MOCA IHP, Paris
BiNoM I/O : Creating and editing BioPAX 26/06/2007 MOCA IHP, Paris
BiNoM I/O : BioPAX Hierarchical pathway structure (useful for MOCA) pathway ‘Pathway structure’ layer Pathway step 26/06/2007 MOCA IHP, Paris
Module depends on the question… (MOCA, 25/06/07) • Divisive (Holme et al, 2003) and agglomerative (Ma et al, 2004) hierarchical clustering based on the shortest path analysis • Eigenpathways (Barret et al, 2006) - Principal component analysis of extreme pathway sampling • Conservation laws (P-invariants) (what basis to chose?) • Exactly solvable modules: • Gorban, Radulescu (2007) “Dynamic and static limitation in multiscale reaction networks”Adv.Chem.Eng (In press and in Arxiv.org) • Radulescu, Zinovyev, Lilienbaum (2007) Model Reduction And Model Comparison For Nf·kb Signaling. Proceedings of FOSBE 2007. • 5) Monotonic modules, timescale separation, stratification of attractors… • Radulescu, Gorban, Vakulenko, Zinovyev. Hierarchies and modules in complex biological systems. Proceedings of ECCS’06, 2006. 26/06/2007 MOCA IHP, Paris
RB-pathway test caseStep 1: Importing CellDesigner File in BiNoM 26/06/2007 MOCA IHP, Paris
In- and Out- digraph layers Bow-tie network structure (Ma et al., 2003) Cyclic part (Strongly connected components) IN OUT Network Output Boundary Conditions IS 26/06/2007 MOCA IHP, Paris
RB-pathway test caseStep 2: Pruning reaction graph 26/06/2007 MOCA IHP, Paris
RB-pathway test caseStep 3: Strongly connected components (Tarjan’s algorithm) 26/06/2007 MOCA IHP, Paris
RB-pathway test caseStep 4: Representing modular structure 26/06/2007 MOCA IHP, Paris
Relevant cycles in chemical reaction networksP.Gleiss, P.Stadler, A.Wagner, D.Fell Adv. Complex Systems (2001) MCB1 MCB3 MCB2 Relevant cycles = union of minimum cyclic bases (MSC) It is unique minimumcyclic decomposition of a graph Vismara’s algorithmis able to deal with families of relevant cycles 26/06/2007 MOCA IHP, Paris
Cycles in reaction networks as functional unitscase 1: mass flow Cycle is an elementary mode (FBA) or T-invariant 26/06/2007 MOCA IHP, Paris
FBA-like path analysis can be also used (Klamt et al., 2006) Cycles in reaction networks as functional unitscase 2: “information” (perturbation) flow 26/06/2007 MOCA IHP, Paris
IN OUT RB-pathway test caseStep 5: Decomposition of Giant Strong Component 15 relevant cycles Agglomerative clustering number of reactions 5 cyclic modules with >50% node intersection Other nodes from IN and OUT part of the graph are connected to cyclic modules 26/06/2007 MOCA IHP, Paris
Modular decomposition 26/06/2007 MOCA IHP, Paris
Modular decomposition (RB-E2F1 module interface) 26/06/2007 MOCA IHP, Paris
What else? • Using semantic information ? • Being a bit more quantitative ? 26/06/2007 MOCA IHP, Paris