Multi-level model for the investigation of oncoantigen- driven vaccination effect

BACKGROUND: Cancer stem cell theory suggests that cancers are derived by a population of cells named Cancer Stem Cells (CSCs) that are involved in the growth and in the progression of tumors, and lead to a hierarchical structure characterized by differentiated cell population. This cell heterogeneity affects the choice of cancer therapies, since many current cancer treatments have limited or no impact at all on CSC population, while they reveal a positive effect on the differentiated cell populations. RESULTS: In this paper we investigated the effect of vaccination on a cancer hierarchical structure through a multi-level model representing both population and molecular aspects. The population level is modeled by a system of Ordinary Differential Equations (ODEs) describing the cancer population's dynamics. The molecular level is modeled using the Petri Net (PN) formalism to detail part of the proliferation pathway. Moreover, we propose a new methodology which exploits the temporal behavior derived from the molecular level to parameterize the ODE system modeling populations. Using this multi-level model we studied the ErbB2-driven vaccination effect in breast cancer. CONCLUSIONS: We propose a multi-level model that describes the inter-dependencies between population and genetic levels, and that can be efficiently used to estimate the efficacy of drug and vaccine therapies in cancer models, given the availability of molecular data on the cancer driving force

Engineering simulations for cancer systems biology

Computer simulation can be used to inform in vivo and in vitro experimentation, enabling rapid, low-cost hypothesis generation and directing experimental design in order to test those hypotheses. In this way, in silico models become a scientific instrument for investigation, and so should be developed to high standards, be carefully calibrated and their findings presented in such that they may be reproduced. Here, we outline a framework that supports developing simulations as scientific instruments, and we select cancer systems biology as an exemplar domain, with a particular focus on cellular signalling models. We consider the challenges of lack of data, incomplete knowledge and modelling in the context of a rapidly changing knowledge base. Our framework comprises a process to clearly separate scientific and engineering concerns in model and simulation development, and an argumentation approach to documenting models for rigorous way of recording assumptions and knowledge gaps. We propose interactive, dynamic visualisation tools to enable the biological community to interact with cellular signalling models directly for experimental design. There is a mismatch in scale between these cellular models and tissue structures that are affected by tumours, and bridging this gap requires substantial computational resource. We present concurrent programming as a technology to link scales without losing important details through model simplification. We discuss the value of combining this technology, interactive visualisation, argumentation and model separation to support development of multi-scale models that represent biologically plausible cells arranged in biologically plausible structures that model cell behaviour, interactions and response to therapeutic interventions

Pathway relevance ranking for tumor samples through network-based data integration

The study of cancer, a highly heterogeneous disease with different causes and clinical outcomes, requires a multi-angle approach and the collection of large multi-omics datasets that, ideally, should be analyzed simultaneously. We present a new pathway relevance ranking method that is able to prioritize pathways according to the information contained in any combination of tumor related omics datasets. Key to the method is the conversion of all available data into a single comprehensive network representation containing not only genes but also individual patient samples. Additionally, all data are linked through a network of previously identified molecular interactions. We demonstrate the performance of the new method by applying it to breast and ovarian cancer datasets from The Cancer Genome Atlas. By integrating gene expression, copy number, mutation and methylation data, the method's potential to identify key pathways involved in breast cancer development shared by different molecular subtypes is illustrated. Interestingly, certain pathways were ranked equally important for different subtypes, even when the underlying (epi)-genetic disturbances were diverse. Next to prioritizing universally high-scoring pathways, the pathway ranking method was able to identify subtype-specific pathways. Often the score of a pathway could not be motivated by a single mutation, copy number or methylation alteration, but rather by a combination of genetic and epi-genetic disturbances, stressing the need for a network-based data integration approach. The analysis of ovarian tumors, as a function of survival-based subtypes, demonstrated the method's ability to correctly identify key pathways, irrespective of tumor subtype. A differential analysis of survival-based subtypes revealed several pathways with higher importance for the bad-outcome patient group than for the good-outcome patient group. Many of the pathways exhibiting higher importance for the bad-outcome patient group could be related to ovarian tumor proliferation and survival

Topological Analysis of MAPK Cascade for Kinetic ErbB Signaling

Ligand-induced homo- and hetero-dimer formation of ErbB receptors results in different biological outcomes irrespective of recruitment and activation of similar effector proteins. Earlier experimental research indicated that cells expressing both EGFR (epidermal growth factor receptor) and the ErbB4 receptor (E1/4 cells) induced E1/4 cell-specific B-Raf activation and higher extracellular signal-regulated kinase (ERK) activation, followed by cellular transformation, than cells solely expressing EGFR (E1 cells) in Chinese hamster ovary (CHO) cells. Since our experimental data revealed the presence of positive feedback by ERK on upstream pathways, it was estimated that the cross-talk/feedback pathway structure of the Raf-MEK-ERK cascade might affect ERK activation dynamics in our cell system. To uncover the regulatory mechanism concerning the ERK dynamics, we used topological models and performed parameter estimation for all candidate structures that possessed ERK-mediated positive feedback regulation of Raf. The structure that reliably reproduced a series of experimental data regarding signal amplitude and duration of the signaling molecules was selected as a solution. We found that the pathway structure is characterized by ERK-mediated positive feedback regulation of B-Raf and B-Raf-mediated negative regulation of Raf-1. Steady-state analysis of the estimated structure indicated that the amplitude of Ras activity might critically affect ERK activity through ERK-B-Raf positive feedback coordination with sustained B-Raf activation in E1/4 cells. However, Rap1 that positively regulates B-Raf activity might be less effective concerning ERK and B-Raf activity. Furthermore, we investigated how such Ras activity in E1/4 cells can be regulated by EGFR/ErbB4 heterodimer-mediated signaling. From a sensitivity analysis of the detailed upstream model for Ras activation, we concluded that Ras activation dynamics is dominated by heterodimer-mediated signaling coordination with a large initial speed of dimerization when the concentration of the ErbB4 receptor is considerably high. Such characteristics of the signaling cause the preferential binding of the Grb2-SOS complex to heterodimer-mediated signaling molecules

Phosphoproteomics-Based Modeling Defines the Regulatory Mechanism Underlying Aberrant EGFR Signaling

BACKGROUND: Mutation of the epidermal growth factor receptor (EGFR) results in a discordant cell signaling, leading to the development of various diseases. However, the mechanism underlying the alteration of downstream signaling due to such mutation has not yet been completely understood at the system level. Here, we report a phosphoproteomics-based methodology for characterizing the regulatory mechanism underlying aberrant EGFR signaling using computational network modeling. METHODOLOGY/PRINCIPAL FINDINGS: Our phosphoproteomic analysis of the mutation at tyrosine 992 (Y992), one of the multifunctional docking sites of EGFR, revealed network-wide effects of the mutation on EGF signaling in a time-resolved manner. Computational modeling based on the temporal activation profiles enabled us to not only rediscover already-known protein interactions with Y992 and internalization property of mutated EGFR but also further gain model-driven insights into the effect of cellular content and the regulation of EGFR degradation. Our kinetic model also suggested critical reactions facilitating the reconstruction of the diverse effects of the mutation on phosphoproteome dynamics. CONCLUSIONS/SIGNIFICANCE: Our integrative approach provided a mechanistic description of the disorders of mutated EGFR signaling networks, which could facilitate the development of a systematic strategy toward controlling disease-related cell signaling

Multiscale Modeling of the ErbB Receptor Tyrosine Kinase Signaling Network Through Theory and Experiment

The biochemical processes occurring within a living cell span a spectrum of scales in space and time, ranging from the nano- to the macro-scale. We note that a single cellular process often operates on multiple spatial and temporal scales, and thus it becomes necessary to combine modeling techniques in multiscale approaches, in which different levels of theory are synergized to describe a system at a number of scales or resolutions. In this work we apply a multiscale modeling framework to investigate the molecular regulatory mechanisms governing the activation of the ErbB receptor tyrosine kinases, a family of kinases which are commonly over-expressed or mutated in human cancers, with a focus on the HER3 and HER4 kinases. Our multiscale model of HER3, a kinase which, until recently, has been considered kinase-dead, presents evidence of HER3 catalytic activity and demonstrates that even a weak HER3 signal can be amplified by other cellular signaling mechanisms to induce drug resistance to tyrosine kinase inhibitors in silico. Thus HER3, rather than the commonly-targeted EGFR and HER2 kinases, may represent a superior therapeutic target in specific ErbB-driven cancers. In the second major study, we construct a multiscale model of activity in the HER4 kinase, which has been shown to perform an anti-cancer role in certain tumor cells, by steering the cell toward a program of cellular differentiation and away from a program of uncontrolled proliferation. Our HER4 model, which applies a combined computational and experimental approach, elucidates the molecular mechanisms underlying this HER4-mediated ‘switch’ to the cellular differentiation program, with the ultimate aim of exploiting or modulating the HER4 pathway as a potential therapy in specific ErbB-driven cancers. Furthermore, the model provides structural insights into the effects of several HER4 somatic mutations which have recently been discovered in a subset of cancer patients, and which abrogate the anti-cancer effects of HER4 activity. We have illustrated that multiscale modeling provides a powerful and quantitative platform for investigating the complexity inherent in intracellular signaling pathways and rationalizing the effects of molecular perturbations on downstream signaling events and ultimately, on the cell phenotype