Global alignment of protein-protein interaction networks by graph - PowerPoint PPT Presentation

Global alignment of protein-protein interaction networks by graph matching methods. Mikhail Zaslavkiy 1 Francis Bach 2 Jean-Philippe Vert 1 1 Mines ParisTech / Institut Curie / INSERM 2 INRIA / Ecole normale superieure de Paris Kyushu University, Department of Informatics, July 14, 2009. JP Vert (ParisTech) Global alignment of PPI networks 1 / 47

Outline 1 Identification of functional orthologs 2 Algorithm for constrained global network alignment 3 Algorithms for balanced global network alignment 4 Experiments 5 Conclusion JP Vert (ParisTech) Global alignment of PPI networks 2 / 47

Functional orthologs Species 1 Species 2 f 1 : MKQALAAADDDDAQ... y 1 : MDDDDALGLLLLA... f 2 :MGDXLLMMAALLLL... y 2 : MHHAAKLLDDAS... ... ... Definition Functional orthologs are pairs of proteins directly inherited from a common ancestor and which play functionally equivalent roles. Our goal Automatic identification of functional orthologs (useful for annotation transfer) JP Vert (ParisTech) Global alignment of PPI networks 4 / 47

Identification of functional orthologs by best-best hit Species 1 Species 2 f 1 : MKQDLARIEQFLDALF... y 1 : MSRLPVLLLLQLLVRGA. . . f 2 : MSKLKIAVSDSCPDCF... y 2 : MELAALCRAGLLLALDA. . . ... ...  y 1 y 2  C = f 1 10 50 C ij -BLAST similarity scores   f 2 27 10 Optimal assignment : f 1 → y 2 , f 2 → y 1 JP Vert (ParisTech) Global alignment of PPI networks 5 / 47

Limitations of sequence comparison-based methods y may be the best hit for f , but f may not be the best hit for y ... ( y 1 , f ) and ( y 2 , f ) may produce very similar blast scores... JP Vert (ParisTech) Global alignment of PPI networks 6 / 47

Clusters of orthologs Many programs produce clusters of orthologous genes from sequence comparison only (COG, KEGG, Inparanoid, ...) Several genes of each species may be in the same cluster How to find functional orthologs within the clusters? JP Vert (ParisTech) Global alignment of PPI networks 7 / 47

Ideas to solve ambiguous functional orthologs Increase the similarity of similarity scores / phylogenetic approaches Comparison of expression profiles across species Functional orthologs tend to have more conserved protein-protein interactions (PPI) across species (Bandyopadhyay et al., 2006) JP Vert (ParisTech) Global alignment of PPI networks 8 / 47

Disambiguation by PPI conservation Idea: If we know that y ∗ and f ∗ are functional orthologs, and there exist interactions f ∗ − f and y ∗ − y 2 . Then the assignment y 2 − f is more likely because it conserves one interaction. JP Vert (ParisTech) Global alignment of PPI networks 9 / 47

Disambiguation by PPI conservation Idea: If we know that y ∗ and f ∗ are functional orthologs, and there exist interactions f ∗ − f and y ∗ − y 2 . Then the assignment y 2 − f is more likely because it conserves one interaction. ??? ??? JP Vert (ParisTech) Global alignment of PPI networks 9 / 47

Disambiguation by PPI conservation Idea: If we know that y ∗ and f ∗ are functional orthologs, and there exist interactions f ∗ − f and y ∗ − y 2 . Then the assignment y 2 − f is more likely because it conserves one interaction. ??? ??? PPI PPI JP Vert (ParisTech) Global alignment of PPI networks 9 / 47

Disambiguation by PPI conservation Idea: If we know that y ∗ and f ∗ are functional orthologs, and there exist interactions f ∗ − f and y ∗ − y 2 . Then the assignment y 2 − f is more likely because it conserves one interaction. PPI PPI JP Vert (ParisTech) Global alignment of PPI networks 9 / 47

Extension to PPI networks JP Vert (ParisTech) Global alignment of PPI networks 10 / 47

Extension to PPI networks matchings JP Vert (ParisTech) Global alignment of PPI networks 10 / 47

Extension to PPI networks 3 conserved interactions JP Vert (ParisTech) Global alignment of PPI networks 10 / 47

Global Network Alignment (GNA) 3 conserved interactions Given two PPI networks and the all-vs-all sequence similarity matrix, find a global matching that maximizes the number of conserved interactions subject to: Constraint GNA: matchings only occur within clusters of orthologs. Balanced GNA: the mean sequence similarity between matched pairs is as large as possible. JP Vert (ParisTech) Global alignment of PPI networks 11 / 47

Complexity of the problems (bad news) Both problems are NP-hard for general graphs and similarity matrix. Therefore we must use algorithms that approximately optimize the criteria, e.g: MRF method (Bandyopadhyay et al., MSB 2006 ) for constrained GNA IsoRank (Singh et al., PNAS 2008 ) for balanced GNA We investigate other algorithms for these problems, borrowing ideas from state-of-the-art graph matching algorithms . JP Vert (ParisTech) Global alignment of PPI networks 12 / 47

Constrained GNA Problem Find matchings within the clusters that maximise the number of conserved interactions JP Vert (ParisTech) Global alignment of PPI networks 14 / 47

Graph of clusters induced by PPI JP Vert (ParisTech) Global alignment of PPI networks 15 / 47

Global optimum Proposition If the graph of clusters generated by the PPI has no cycle, then we can find the optimal matching efficiently with a message passing algorithm. JP Vert (ParisTech) Global alignment of PPI networks 16 / 47

Global optimum by message passing (Similar to Viterbi’s algorithm for HMM) JP Vert (ParisTech) Global alignment of PPI networks 17 / 47

What if the graph of clusters has cycle? The message passing method can not be used... Instead we reformulate the constrained GNA problem as a balanced GNA by setting similarity between proteins in different clusters to −∞ , and use algorithms for balanced GNA. JP Vert (ParisTech) Global alignment of PPI networks 18 / 47

Balanced GNA Given two graphs and a matrix of all-vs-all similarities, find a matching P ∈ P that jointly maximizes: the number of conserved interaction CI ( P ), the mean similarity of matched pairs S ( P ). The trade-off can be found by maximizing over P : P ∈P F ( P ) = (1 − α ) CI ( P ) + α S ( P ) , min where α ∈ [0 , 1] determines the balance between both objectives. JP Vert (ParisTech) Global alignment of PPI networks 20 / 47

Existing methods for balanced GNA P ∈P F ( P ) = (1 − α ) CI ( P ) + α S ( P ) , min When α = 1 this is an optimal assignment problem efficiently solved by the Hungarian algorithm (Kuhn, 1955). When α < 1 this is a general graph matching problem, usually computationally intractable. Existing algorithms include: Exact solution by incomplete enumeration (only for small graphs) Spectral methods (Umeyama, 1986; Singh et al., 2008) Relaxations of the problem into a continuous optimization problem (Almohamad and Duffuaa, 1993; Gold and Rangarajan, 1996). JP Vert (ParisTech) Global alignment of PPI networks 21 / 47

Global alignment of protein-protein interaction networks by graph - PowerPoint PPT Presentation

Global alignment of protein-protein interaction networks by graph matching methods. Mikhail Zaslavkiy 1 Francis Bach 2 Jean-Philippe Vert 1 1 Mines ParisTech / Institut Curie / INSERM 2 INRIA / Ecole normale superieure de Paris Kyushu University,

Global alignment of protein-protein interaction networks by graph matching methods. Mikhail

Hasup Lee, Seungtaek Sun and Ye-Yeong Park ( Group 6 ) Protein-Protein interaction is

Collaboration-based Function Prediction in Protein-Protein Interaction networks Hossein Rahmani

Modeling the dynamics and function of cellular interaction networks Rka Albert Department of

INTERACTION NETWORKS USING DEEP LEARNING Peiying (Colleen) Ruan, PhD, Deep Learning Solution

New prenylchalcones targeting the MDM2-p53 protein-protein interaction: synthesis and evaluation

disease mutations in protein interaction networks and implications for personalized medicine

Atomic Resolution Modeling Protein-Protein interaction network of Large Macromolecular

PEPSI-DOCK A Detailed Data-Driven Protein-Protein Interaction Potential Accelerated By

Pattern Matching in Protein-Protein Interaction Graphs Ga elle Brevier ( Universit e de

Novel Motif Detection Algorithms for Finding Protein-Protein Interaction Sites January Wisniewski

Kpax Protein Structure Alignment Demo: Using Kpax on Linux Dave Ritchie Practical: Homology

Protein-Protein interactions Reducing the complexity Why are protein-protein interactions

Biocuration and rule-based modelling of protein interaction networks in KAMI Sbastien Lgar,

DNA RNA Protein synthesis AMINO ACIDS PROTEIN Protein degradation FUNCTION Some properties

Dynamics of Protein-Protein Interactions: A Probabilistic Model Toward Protein Function Amir

Protein Interaction Prediction: The PIPE and InSiPS Projects Frank Dehne School of Computer

PROTEIN EXPRESSION AND PURIFICATION PROTEIN EXPRESSION AND PURIFICATION Why do we decide to

A Graph Modification Approach for Finding CorePeriphery Structures in Protein Interaction

http://vision.unipv.it/ ANALYSIS OF GEMOMETRICAL AND TOPOLOGICAL ATTITUDE FOR PROTEIN-PROTEIN

Protein Folding Protein Folding Proteins have unique 3-dimensional shapes created by the

Protein Folding Protein Folding Proteins have unique 3-dimensional shapes created by the

N C C C protein sequence but is not fully rigid C C peptide C C bond

DNA What is the major compon onent nt of all cells? s? PROTEIN INS Why would protein