1 / 57

Where Computer Science meets Biology

Systems Biology. Where Computer Science meets Biology. Professor Adriana Compagnoni Computer Science Department. With: Joe Glavy, Svetlana Sukhishvili, Tommy White Amanda DiGuilio (CCBBME) Yifei Bao, Vishakha Sharma, Justin Sousa, and Peter Zafonte (CS). Overview. The Virtual Lab

rich
Download Presentation

Where Computer Science meets Biology

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Systems Biology Where Computer Science meets Biology Professor Adriana Compagnoni Computer Science Department With: Joe Glavy, Svetlana Sukhishvili, Tommy White Amanda DiGuilio (CCBBME) Yifei Bao, Vishakha Sharma, Justin Sousa, and Peter Zafonte (CS)

  2. Overview • The Virtual Lab • Modeling HER2 • Modeling Bio-films • Experimenting with 5 different modeling techniques • Modeling the G-protein Cycle • Incorporating Space

  3. Process model for HER2 Signaling Pathways

  4. HER2

  5. HER2 – Signaling Pathways • HER signaling is initiated when two receptor molecules join together in a process called dimerization, which occurs in response to the presence of a growth factor molecule (ligand). • Dimerization of different members of ErbB family, lead to activation of different kinase mediated intracellular signaling cascades, like • PI3 K/AKT pathway • MAPK pathway • JAK/STAT pathway

  6. Process Algebra Model

  7. Overview of different species

  8. Differential HER2 concentrations a. 100 HER2 are initialized a. 30 HER2 are initialized

  9. Inhibition of Raf activation(activation rate is changed from 27.0 to 1.0E-4 ) Pathways Crosstalk

  10. Computational Modeling of Bio-films for Drug Delivery

  11. 3.2 μm 3.2 μm Hydrogen-bonded multilayer destruction Changing pH fast release of capsule cargo

  12. Computational Modeling For The G-protein Cycle • Using 5 different modeling techniques: • Differential equations • Stochastic Pi-calculus • Stochastic Petri-Nets • Kappa • Bio-Pepa 12

  13. G-Protein Couple Receptors(GPCRs) • G-protein couple receptors (GPCRs) sense molecules outside the cell and activate inside signal transduction pathways. • An estimated 50% of the current pharmaceuticals target GPCRs. 13

  14. Activation cycle of G-proteins by G-protein-coupled receptors 2 1 3 5 4 14

  15. Chemical Reactions and Rates 15

  16. The law of mass action • The law of mass action is a mathematical model that prescribes the evolution of a chemical system in terms of changes of concentrations of the chemical species over time. • In its simplest form, it says that a reaction X+Yk Z has a rate k[X][Y]: the rate is proportional to the concentration of one species ([X]) times the concentration of the other species ([Y]) by the base rate k. 16

  17. Ordinary Differential Equations (ODEs) 17

  18. Process AlgebrasAn alternative ot ODEs Formal languages originally designed to model complex reactive computer systems. Because of the similarities between reactive computer systems and biological systems, process algebra have recently been used to model biological systems. 18

  19. Process Algebras Typically two halves of a communication: Sending and Receiving. !ch(msg) : to send message msg on channel ch. ?ch(msg) : to receive a message msg on channel ch. 19

  20. Process Algebras A = new msg !ch (msg); Ab(msg) B = ?ch(msg); Bb (msg) Ab(msg) = !msg(); A Bb(msg) = ?msg(); B A and B bind A and B dissociate 20

  21. Stochastic Pi-Calculus The stochastic pi-calculus is a process algebra where stochastic rates are imposed on processes, allowing a more accurate description of biological processes. SPiM (stochastic pi machine) is an implementation of the stochastic pi-calculus that can be used to run in-silico simulations that display the change over time in the populations of the different species of the system being modeled.

  22. Step 1 to Step 2 2 1 !bindb Gd G ?bindb Gbg Gd represents alpha Gbg represents the beta-gamma complex 22

  23. Process modeling for G-protein Cycle 23

  24. Process modeling for G-protein Cycle 24

  25. Process modeling for G-protein Cycle 25

  26. Process modeling for G-protein Cycle 26

  27. Process modeling for G-protein Cycle 27

  28. Process modeling for G-protein Cycle 28

  29. SPiM directive plot RL(); abrD(); aT()

  30. Petri Nets Modeling The basic Petri Net is a directed bipartite graph with two kinds of nodes which are either places or transitions and directed arcs which connect nodes. In modeling biological processes, place nodes represent molecular species and transition nodes represent reactions. We use Cell Illustrator to develop our model of the G-protein cycle (www.cellillustrator.com). 30

  31. Petri Nets Modeling NJPLS2010@Stevens

  32. Petri Nets Modeling NJPLS2010@Stevens

  33. Petri Nets Modeling NJPLS2010@Stevens

  34. Petri Nets Modeling NJPLS2010@Stevens

  35. Petri Nets Modeling NJPLS2010@Stevens

  36. Petri Nets Modeling NJPLS2010@Stevens

  37. Kappa Language Modeling Kappa is a formal language for defining agents (typically meant to represent proteins) as sets of sites that constitute abstract resources for interaction. It is used to express rules of interactions between proteins characterized by discrete modification and binding states. The Kappa language modeling platform is Cellucidate.

  38. Kappa Language Modeling R and L are agent names, r is the binding site of R, l is the binding site of L, and 3.32e-09 is the reaction rate. In R(r!1) and L(l!1), 1 is the index of the link that binds R and L at their binding sites. R + L -> RL R(r), L(l)->R(r!1), L(l!1) @ 3.32e-09

  39. Kappa Language Modeling

  40. Bio-PEPA Modeling Bio-PEPA is a process algebra for the modeling and the analysis of biochemical networks. It is a modification of PEPA to deal with some features of biological models, such as stoichiometry and the use of generic kinetic laws.

  41. Bio-PEPA Modeling

  42. Bio-PEPA Modeling

  43. Result (ODEs and Pi-calculus)

  44. Result (Petri Nets and Kappa)

  45. Result (Bio-PEPA)

  46. High Level Notations (motivation) Both ODEs and Petri Nets correspond closely to chemical reactions, and for the average biologist, they are relatively easy to understand. Cellucidate provides a friendly user interface for Kappa that abstracts away from its syntax. SPiM still needs encoding in the stochastic Pi-calculus.

  47. High Level Notation (motivation) In order to hide Pi-calculus communication primitives and enable modeling using only terminology directly obtained from biological processes, we develop a high level notation that can be systematically translated into SPiM programs.

  48. G-protein Cycle 2 bind bind 1 3 activate and dissociate hydrolyze 5 4 dissociate 48

  49. Step 1 to Step 2 2 1 bind(Gd, Gbg, G, 1.0) 49

  50. High Level Notation bind(Gd, Gbg, G, 1.0); bind(R, L, RL, 3.32e-6); activateAnddissociate(G, RL, Ga, Gbg, 1.0e-5); dissociate(RL, R, L, 0.01); hydrolyze(Ga, Gd, 0.11); degrade(R, 4e-4); degrade(RL, 4e-3)

More Related