Clinically driven design of multi-scale cancer models: the ContraCancrum project paradigm

Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics

Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.

Multi-scale network targeting: A holistic systems-biology approach to cancer treatment

The vulnerabilities of cancer at the cellular and, recently, with the introduction of immunotherapy, at the tissue level, have been exploited with variable success. Evaluating the cancer system vulnerabilities at the organismic level through analysis of network topology and network dynamics can potentially predict novel anti-cancer drug targets directed at the macroscopic cancer networks. Theoretical work analyzing the properties and the vulnerabilities of the multi-scale network of cancer needs to go hand-in-hand with experimental research that uncovers the biological nature of the relevant networks and reveals new targetable vulnerabilities. It is our hope that attacking cancer on different spatial scales, in a concerted integrated approach, may present opportunities for novel ways to prevent treatment resistance.

High-throughput binding affinity calculations at extreme scales

Background

Resistance to chemotherapy and molecularly targeted therapies is a major factor in limiting the effectiveness of cancer treatment. In many cases, resistance can be linked to genetic changes in target proteins, either pre-existing or evolutionarily selected during treatment. Key to overcoming this challenge is an understanding of the molecular determinants of drug binding. Using multi-stage pipelines of molecular simulations we can gain insights into the binding free energy and the residence time of a ligand, which can inform both stratified and personal treatment regimes and drug development. To support the scalable, adaptive and automated calculation of the binding free energy on high-performance computing resources, we introduce the High-throughput Binding Affinity Calculator (HTBAC). HTBAC uses a building block approach in order to attain both workflow flexibility and performance.

Results

We demonstrate close to perfect weak scaling to hundreds of concurrent multi-stage binding affinity calculation pipelines. This permits a rapid time-to-solution that is essentially invariant of the calculation protocol, size of candidate ligands and number of ensemble simulations.

Conclusions

As such, HTBAC advances the state of the art of binding affinity calculations and protocols. HTBAC provides the platform to enable scientists to study a wide range of cancer drugs and candidate ligands in order to support personalized clinical decision making based on genome sequencing and drug discovery.

A study on the predictability of acute lymphoblastic leukaemia response to treatment using a hybrid oncosimulator

Efficient use of Virtual Physiological Human (VPH)-type models for personalized treatment response prediction purposes requires a precise model parameterization. In the case where the available personalized data are not sufficient to fully determine the parameter values, an appropriate prediction task may be followed. This study, a hybrid combination of computational optimization and machine learning methods with an already developed mechanistic model called the acute lymphoblastic leukaemia (ALL) Oncosimulator which simulates ALL progression and treatment response is presented. These methods are used in order for the parameters of the model to be estimated for retrospective cases and to be predicted for prospective ones. The parameter value prediction is based on a regression model trained on retrospective cases. The proposed Hybrid ALL Oncosimulator system has been evaluated when predicting the pre-phase treatment outcome in ALL. This has been correctly achieved for a significant percentage of patient cases tested (approx. 70% of patients). Moreover, the system is capable of denying the classification of cases for which the results are not trustworthy enough. In that case, potentially misleading predictions for a number of patients are avoided, while the classification accuracy for the remaining patient cases further increases. The results obtained are particularly encouraging regarding the soundness of the proposed methodologies and their relevance to the process of achieving clinical applicability of the proposed Hybrid ALL Oncosimulator system and VPH models in general.

E-Science technologies in a workflow for personalized medicine using cancer screening as a case study

Objective

We provide an e-Science perspective on the workflow from risk factor discovery and classification of disease to evaluation of personalized intervention programs. As case studies, we use personalized prostate and breast cancer screenings.

Materials and Methods

We describe an e-Science initiative in Sweden, e-Science for Cancer Prevention and Control (eCPC), which supports biomarker discovery and offers decision support for personalized intervention strategies. The generic eCPC contribution is a workflow with 4 nodes applied iteratively, and the concept of e-Science signifies systematic use of tools from the mathematical, statistical, data, and computer sciences.

Results

The eCPC workflow is illustrated through 2 case studies. For prostate cancer, an in-house personalized screening tool, the Stockholm-3 model (S3M), is presented as an alternative to prostate-specific antigen testing alone. S3M is evaluated in a trial setting and plans for rollout in the population are discussed. For breast cancer, new biomarkers based on breast density and molecular profiles are developed and the US multicenter Women Informed to Screen Depending on Measures (WISDOM) trial is referred to for evaluation. While current eCPC data management uses a traditional data warehouse model, we discuss eCPC-developed features of a coherent data integration platform.

Discussion and Conclusion

E-Science tools are a key part of an evidence-based process for personalized medicine. This paper provides a structured workflow from data and models to evaluation of new personalized intervention strategies. The importance of multidisciplinary collaboration is emphasized. Importantly, the generic concepts of the suggested eCPC workflow are transferrable to other disease domains, although each disease will require tailored solutions.

Clinical implications of in silico mathematical modeling for glioblastoma: a critical review

Glioblastoma remains a clinical challenge in spite of years of extensive research. Novel approaches are needed in order to integrate the existing knowledge. This is the potential role of mathematical oncology. This paper reviews mathematical models on glioblastoma from the clinical doctor's point of view, with focus on 3D modeling approaches of radiation response of in vivo glioblastomas based on contemporary imaging techniques. As these models aim to provide a clinically useful tool in the era of personalized medicine, the integration of the latest advances in molecular and imaging science and in clinical practice by the in silico models is crucial for their clinical relevance. Our aim is to indicate areas of GBM research that have not yet been addressed by in silico models and to point out evidence that has come up from in silico experiments, which may be worth considering in the clinic. This review examines how close these models have come in predicting the outcome of treatment protocols and in shaping the future of radiotherapy treatments.

In Silico Oncology: Quantification of the In Vivo Antitumor Efficacy of Cisplatin-Based Doublet Therapy in Non-Small Cell Lung Cancer (NSCLC) through a Multiscale Mechanistic Model

The 5-year survival of non-small cell lung cancer patients can be as low as 1% in advanced stages. For patients with resectable disease, the successful choice of preoperative chemotherapy is critical to eliminate micrometastasis and improve operability. In silico experimentations can suggest the optimal treatment protocol for each patient based on their own multiscale data. A determinant for reliable predictions is the a priori estimation of the drugs’ cytotoxic efficacy on cancer cells for a given treatment. In the present work a mechanistic model of cancer response to treatment is applied for the estimation of a plausible value range of the cell killing efficacy of various cisplatin-based doublet regimens. Among others, the model incorporates the cancer related mechanism of uncontrolled proliferation, population heterogeneity, hypoxia and treatment resistance. The methodology is based on the provision of tumor volumetric data at two time points, before and after or during treatment. It takes into account the effect of tumor microenvironment and cell repopulation on treatment outcome. A thorough sensitivity analysis based on one-factor-at-a-time and latin hypercube sampling/partial rank correlation coefficient approaches has established the volume growth rate and the growth fraction at diagnosis as key features for more accurate estimates. The methodology is applied on the retrospective data of thirteen patients with non-small cell lung cancer who received cisplatin in combination with gemcitabine, vinorelbine or docetaxel in the neoadjuvant context. The selection of model input values has been guided by a comprehensive literature survey on cancer-specific proliferation kinetics. The latin hypercube sampling has been recruited to compensate for patient-specific uncertainties. Concluding, the present work provides a quantitative framework for the estimation of the in-vivo cell-killing ability of various chemotherapies. Correlation studies of such estimates with the molecular profile of patients could serve as a basis for reliable personalized predictions.

Less than 14% of medically treated patients with locally advanced and metastatic non-small cell lung cancer are expected to be alive 5 years after diagnosis. Standard therapeutic strategies include the administration of two drugs in combination, aiming at shrinking the tumor before surgery and improving overall survival. Knowing the sensitivity profile of each patient to different treatment strategies at diagnosis may help choose the most appropriate ones. We develop a methodology for the quantitative estimation of the cytotoxic efficacy of cisplatin-based doublets on cancer cells by applying a simulation model of cancer progression and response. The model incorporates the proliferation cycle, quiescence, differentiation and loss of tumor cells. We evaluate the effect of in vivo microenvironment of real tumors, as expressed by measurable tumor proliferation kinetics, such as how fast the tumor grows, the percentage of cells that are actively dividing, the resistance of stem cells, etc. on treatment outcome so as to derive more accurate estimates. A literature survey guides the selection of values. The methodology is applied to a real clinical dataset of patients. Correlation studies between the derived cytotoxicities and the patients’ molecular profile could lead to predictions of treatment response at the time of diagnosis.

The INTEGRATE project: Delivering solutions for efficient multi-centric clinical research and trials

The objective of the INTEGRATE project (http://www.fp7-integrate.eu/) that has recently concluded successfully was the development of innovative biomedical applications focused on streamlining the execution of clinical research, on enabling multidisciplinary collaboration, on management and large-scale sharing of multi-level heterogeneous datasets, and on the development of new methodologies and of predictive multi-scale models in cancer. In this paper, we present the way the INTEGRATE consortium has approached important challenges such as the integration of multi-scale biomedical data in the context of post-genomic clinical trials, the development of predictive models and the implementation of tools to facilitate the efficient execution of postgenomic multi-centric clinical trials in breast cancer. Furthermore, we provide a number of key "lessons learned" during the process and give directions for further future research and development.

Differentiation resistance through altered retinoblastoma protein function in acute lymphoblastic leukemia: in silico modeling of the deregulations in the G1/S restriction point pathway

BACKGROUND

As in many cancer types, the G1/S restriction point (RP) is deregulated in Acute Lymphoblastic Leukemia (ALL). Hyper-phosphorylated retinoblastoma protein (hyper-pRb) is found in high levels in ALL cells. Nevertheless, the ALL lymphocyte proliferation rate for the average patient is surprisingly low compared to its normal counterpart of the same maturation level. Additionally, as stated in literature, ALL cells possibly reside at or beyond the RP which is located in the late-G1 phase. This state may favor their differentiation resistant phenotype. A major phenomenon contributing to this fact is thought to be the observed limited redundancy in the phosphorylation of retinoblastoma protein (pRb) by the various Cyclin Dependent Kinases (Cdks). The latter may result in partial loss of pRb functions despite hyper-phosphorylation.

RESULTS

To test this hypothesis, an in silico model aiming at simulating the biochemical regulation of the RP in ALL is introduced. By exploiting experimental findings derived from leukemic cells and following a semi-quantitative calibration procedure, the model has been shown to satisfactorily reproduce such a behavior for the RP pathway. At the same time, the calibrated model has been proved to be in agreement with the observed variation in the ALL cell cycle duration.

CONCLUSIONS

The proposed model aims to contribute to a better understanding of the complex phenomena governing the leukemic cell cycle. At the same time it constitutes a significant first step in the creation of a personalized proliferation rate predictor that can be used in the context of multiscale cancer modeling. Such an approach is expected to play an important role in the refinement and the advancement of mechanistic modeling of ALL in the context of the emergent and promising scientific domains of In Silico Oncology and more generally In Silico Medicine.

Integrating Multiscale Modeling with Drug Effects for Cancer Treatment

In this paper, we review multiscale modeling for cancer treatment with the incorporation of drug effects from an applied system's pharmacology perspective. Both the classical pharmacology and systems biology are inherently quantitative; however, systems biology focuses more on networks and multi factorial controls over biological processes rather than on drugs and targets in isolation, whereas systems pharmacology has a strong focus on studying drugs with regard to the pharmacokinetic (PK) and pharmacodynamic (PD) relations accompanying drug interactions with multiscale physiology as well as the prediction of dosage-exposure responses and economic potentials of drugs. Thus, it requires multiscale methods to address the need for integrating models from the molecular levels to the cellular, tissue, and organism levels. It is a common belief that tumorigenesis and tumor growth can be best understood and tackled by employing and integrating a multifaceted approach that includes in vivo and in vitro experiments, in silico models, multiscale tumor modeling, continuous/discrete modeling, agent-based modeling, and multiscale modeling with PK/PD drug effect inputs. We provide an example application of multiscale modeling employing stochastic hybrid system for a colon cancer cell line HCT-116 with the application of Lapatinib drug. It is observed that the simulation results are similar to those observed from the setup of the wet-lab experiments at the Translational Genomics Research Institute.

Semantic biomedical resource discovery: a Natural Language Processing framework

BACKGROUND

A plethora of publicly available biomedical resources do currently exist and are constantly increasing at a fast rate. In parallel, specialized repositories are been developed, indexing numerous clinical and biomedical tools. The main drawback of such repositories is the difficulty in locating appropriate resources for a clinical or biomedical decision task, especially for non-Information Technology expert users. In parallel, although NLP research in the clinical domain has been active since the 1960s, progress in the development of NLP applications has been slow and lags behind progress in the general NLP domain. The aim of the present study is to investigate the use of semantics for biomedical resources annotation with domain specific ontologies and exploit Natural Language Processing methods in empowering the non-Information Technology expert users to efficiently search for biomedical resources using natural language.

METHODS

A Natural Language Processing engine which can "translate" free text into targeted queries, automatically transforming a clinical research question into a request description that contains only terms of ontologies, has been implemented. The implementation is based on information extraction techniques for text in natural language, guided by integrated ontologies. Furthermore, knowledge from robust text mining methods has been incorporated to map descriptions into suitable domain ontologies in order to ensure that the biomedical resources descriptions are domain oriented and enhance the accuracy of services discovery. The framework is freely available as a web application at ( http://calchas.ics.forth.gr/ ).

RESULTS

For our experiments, a range of clinical questions were established based on descriptions of clinical trials from the ClinicalTrials.gov registry as well as recommendations from clinicians. Domain experts manually identified the available tools in a tools repository which are suitable for addressing the clinical questions at hand, either individually or as a set of tools forming a computational pipeline. The results were compared with those obtained from an automated discovery of candidate biomedical tools. For the evaluation of the results, precision and recall measurements were used. Our results indicate that the proposed framework has a high precision and low recall, implying that the system returns essentially more relevant results than irrelevant.

CONCLUSIONS

There are adequate biomedical ontologies already available, sufficiency of existing NLP tools and quality of biomedical annotation systems for the implementation of a biomedical resources discovery framework, based on the semantic annotation of resources and the use on NLP techniques. The results of the present study demonstrate the clinical utility of the application of the proposed framework which aims to bridge the gap between clinical question in natural language and efficient dynamic biomedical resources discovery.

The standardized histogram shift of t2 magnetic resonance image (MRI) signal intensities of nephroblastoma does not predict histopathological diagnostic information

The objective of this study is to assess standardized histograms of signal intensities of T2-weighted magnetic resonance image (MRI) modality before and after preoperative chemotherapy for nephroblastoma (Wilms' tumor). All analyzed patients are enrolled in the International Society of Paediatric Oncology (SIOP) 2001/GPOH trial.1 The question to be answered is whether the comparison of the histograms can add new knowledge by comparing them with the histology of the tumor after preoperative chemotherapy. Twenty-three unilateral nephroblastoma cases were analyzed. All patients were examined by MRI before and after preoperative chemotherapy treatment. T2 modalities of the MRIs were selected, and histogram changes were compared to histopathological data available after surgery. Of the 23 tumors, 22 decreased in volume following chemotherapy (median -57.99%; range 15.65 to -90.82%). The preliminary results suggest that standardized histograms of signal intensities of T2 MRI in nephroblastoma is not predicting histopathological diagnostic information and has no implications for the clinical assessment for further chemotherapy.

Integrated PK-PD and agent-based modeling in oncology

Mathematical modeling has become a valuable tool that strives to complement conventional biomedical research modalities in order to predict experimental outcome, generate new medical hypotheses, and optimize clinical therapies. Two specific approaches, pharmacokinetic-pharmacodynamic (PK-PD) modeling, and agent-based modeling (ABM), have been widely applied in cancer research. While they have made important contributions on their own (e.g., PK-PD in examining chemotherapy drug efficacy and resistance, and ABM in describing and predicting tumor growth and metastasis), only a few groups have started to combine both approaches together in an effort to gain more insights into the details of drug dynamics and the resulting impact on tumor growth. In this review, we focus our discussion on some of the most recent modeling studies building on a combined PK-PD and ABM approach that have generated experimentally testable hypotheses. Some future directions are also discussed.

A survey of MRI-based medical image analysis for brain tumor studies

MRI-based medical image analysis for brain tumor studies is gaining attention in recent times due to an increased need for efficient and objective evaluation of large amounts of data. While the pioneering approaches applying automated methods for the analysis of brain tumor images date back almost two decades, the current methods are becoming more mature and coming closer to routine clinical application. This review aims to provide a comprehensive overview by giving a brief introduction to brain tumors and imaging of brain tumors first. Then, we review the state of the art in segmentation, registration and modeling related to tumor-bearing brain images with a focus on gliomas. The objective in the segmentation is outlining the tumor including its sub-compartments and surrounding tissues, while the main challenge in registration and modeling is the handling of morphological changes caused by the tumor. The qualities of different approaches are discussed with a focus on methods that can be applied on standard clinical imaging protocols. Finally, a critical assessment of the current state is performed and future developments and trends are addressed, giving special attention to recent developments in radiological tumor assessment guidelines.

Patient-specific modelling in drug design, development and selection including its role in clinical decision-making

Genomics has made enormous progress in the twelve years since the publication of the first draft human genome sequence, but it has not yet been translated into the clinic. Despite spiralling development costs, the number of new drug registrations is not increasing. One reason for this lies in the genetic complexity of disease. Most diseases involve dysregulation in pathways that involve many genes, and many (including most cancers) are themselves genetically heterogeneous. Systems biology involves the multi-level simulation of physiology, cell biology and biochemistry using complex computational techniques. We show here using case studies in cancer and HIV how such computational models, and particularly models based on individual patient data, can be used for drug design and development, and in the selection of the appropriate treatment for a given patient in the face of resistance mutations. If these techniques are to be adopted in routine clinical practice, clinicians will need better training in modern approaches to the integrated analysis of large-scale heterogeneous data and multi-scale models, while developers will need to provide much more usable tools. Investment in computational infrastructure is needed so that results can be returned on clinically relevant timescales and data warehouses designed with data protection as well as accessibility in mind.

Simulating radiotherapy effect in high-grade glioma by using diffusive modeling and brain atlases

Applying diffusive models for simulating the spatiotemporal change of concentration of tumour cells is a modern application of predictive oncology. Diffusive models are used for modelling glioblastoma, the most aggressive type of glioma. This paper presents the results of applying a linear quadratic model for simulating the effects of radiotherapy on an advanced diffusive glioma model. This diffusive model takes into consideration the heterogeneous velocity of glioma in gray and white matter and the anisotropic migration of tumor cells, which is facilitated along white fibers. This work uses normal brain atlases for extracting the proportions of white and gray matter and the diffusion tensors used for anisotropy. The paper also presents the results of applying this glioma model on real clinical datasets.

From base pair to bedside: molecular simulation and the translation of genomics to personalized medicine

Despite the promises made that genomic sequencing would transform therapy by introducing a new era of personalized medicine, relatively few tangible breakthroughs have been made. This has led to the recognition that complex interactions at multiple spatial, temporal, and organizational levels may often combine to produce disease. Understanding this complexity requires that existing and future models are used and interpreted within a framework that incorporates knowledge derived from investigations at multiple levels of biological function. It also requires a computational infrastructure capable of dealing with the vast quantities of data generated by genomic approaches. In this review, we discuss the use of molecular modeling to generate quantitative and qualitative insights at the smallest scales of the systems biology hierarchy, how it can play an important role in the development of a systems understanding of disease and in the application of such knowledge to help discover new therapies and target existing ones on a personal level.

High-grade glioma diffusive modeling using statistical tissue information and diffusion tensors extracted from atlases

Glioma, especially glioblastoma, is a leading cause of brain cancer fatality involving highly invasive and neoplastic growth. Diffusive models of glioma growth use variations of the diffusion-reaction equation in order to simulate the invasive patterns of glioma cells by approximating the spatiotemporal change of glioma cell concentration. The most advanced diffusive models take into consideration the heterogeneous velocity of glioma in gray and white matter, by using two different discrete diffusion coefficients in these areas. Moreover, by using diffusion tensor imaging (DTI), they simulate the anisotropic migration of glioma cells, which is facilitated along white fibers, assuming diffusion tensors with different diffusion coefficients along each candidate direction of growth. Our study extends this concept by fully exploiting the proportions of white and gray matter extracted by normal brain atlases, rather than discretizing diffusion coefficients. Moreover, the proportions of white and gray matter, as well as the diffusion tensors, are extracted by the respective atlases; thus, no DTI processing is needed. Finally, we applied this novel glioma growth model on real data and the results indicate that prognostication rates can be improved.

The Virtual Physiological Human

The Virtual Physiological Human is synonymous with a programme in computational biomedicine that aims to develop a framework of methods and technologies to investigate the human body as a whole. It is predicated on the transformational character of information technology, brought to bear on that most crucial of human concerns, our own health and well-being.

Fully automatic segmentation of brain tumor images using support vector machine classification in combination with hierarchical conditional random field regularization

Delineating brain tumor boundaries from magnetic resonance images is an essential task for the analysis of brain cancer. We propose a fully automatic method for brain tissue segmentation, which combines Support Vector Machine classification using multispectral intensities and textures with subsequent hierarchical regularization based on Conditional Random Fields. The CRF regularization introduces spatial constraints to the powerful SVM classification, which assumes voxels to be independent from their neighbors. The approach first separates healthy and tumor tissue before both regions are subclassified into cerebrospinal fluid, white matter, gray matter and necrotic, active, edema region respectively in a novel hierarchical way. The hierarchical approach adds robustness and speed by allowing to apply different levels of regularization at different stages. The method is fast and tailored to standard clinical acquisition protocols. It was assessed on 10 multispectral patient datasets with results outperforming previous methods in terms of segmentation detail and computation times.