Optimizing dosage-specific treatments in a multi-Scale model of a tumor growth

The emergence of cell resistance in cancer treatment is a complex phenomenon that emerges from the interplay of processes that occur at different scales. For instance, molecular mechanisms and population-level dynamics such as competition and cell–cell variability have been described as playing a key role in the emergence and evolution of cell resistances. Multi-scale models are a useful tool for studying biology at very different times and spatial scales, as they can integrate different processes occurring at the molecular, cellular, and intercellular levels. In the present work, we use an extended hybrid multi-scale model of 3T3 fibroblast spheroid to perform a deep exploration of the parameter space of effective treatment strategies based on TNF pulses. To explore the parameter space of effective treatments in different scenarios and conditions, we have developed an HPC-optimized model exploration workflow based on EMEWS. We first studied the effect of the cells’ spatial distribution in the values of the treatment parameters by optimizing the supply strategies in 2D monolayers and 3D spheroids of different sizes. We later study the robustness of the effective treatments when heterogeneous populations of cells are considered. We found that our model exploration workflow can find effective treatments in all the studied conditions. Our results show that cells’ spatial geometry and population variability should be considered when optimizing treatment strategies in order to find robust parameter sets.This research has received funding from the Horizon 2020 INFORE Project, GA n° 825070 and the Horizon 2020 PerMedCoE Project, GA n° 951773.Peer ReviewedPostprint (published version

Parallel model exploration for tumor treatment simulations

Abstract Computational systems and methods are often being used in biological research, including the understanding of cancer and the development of treatments. Simulations of tumor growth and its response to different drugs are of particular importance, but also challenging complexity. The main challenges are first to calibrate the simulators so as to reproduce real-world cases, and second, to search for specific values of the parameter space concerning effective drug treatments. In this work, we combine a multi-scale simulator for tumor cell growth and a genetic algorithm (GA) as a heuristic search method for finding good parameter configurations in reasonable time. The two modules are integrated into a single workflow that can be executed in parallel on high performance computing infrastructures. In effect, the GA is used to calibrate the simulator, and then to explore different drug delivery schemes. Among these schemes, we aim to find those that minimize tumor cell size and the probability of emergence of drug resistant cells in the future. Experimental results illustrate the effectiveness and computational efficiency of the approach.This work has received funding from the EU Horizon 2020 RIA program INFORE under grant agreement No 825070Peer ReviewedPostprint (author's final draft

A generalizable data-driven multicellular model of pancreatic ductal adenocarcinoma.

BACKGROUND: Mechanistic models, when combined with pertinent data, can improve our knowledge regarding important molecular and cellular mechanisms found in cancer. These models make the prediction of tissue-level response to drug treatment possible, which can lead to new therapies and improved patient outcomes. Here we present a data-driven multiscale modeling framework to study molecular interactions between cancer, stromal, and immune cells found in the tumor microenvironment. We also develop methods to use molecular data available in The Cancer Genome Atlas to generate sample-specific models of cancer. RESULTS: By combining published models of different cells relevant to pancreatic ductal adenocarcinoma (PDAC), we built an agent-based model of the multicellular pancreatic tumor microenvironment, formally describing cell type-specific molecular interactions and cytokine-mediated cell-cell communications. We used an ensemble-based modeling approach to systematically explore how variations in the tumor microenvironment affect the viability of cancer cells. The results suggest that the autocrine loop involving EGF signaling is a key interaction modulator between pancreatic cancer and stellate cells. EGF is also found to be associated with previously described subtypes of PDAC. Moreover, the model allows a systematic exploration of the effect of possible therapeutic perturbations; our simulations suggest that reducing bFGF secretion by stellate cells will have, on average, a positive impact on cancer apoptosis. CONCLUSIONS: The developed framework allows model-driven hypotheses to be generated regarding therapeutically relevant PDAC states with potential molecular and cellular drivers indicating specific intervention strategies

Metameric representations on optimization of nano particle cancer treatment

In silico evolutionary optimization of cancer treatment based on multiple nano-particle (NP) assisted drug delivery systems was investigated in this study. The use of multiple types of NPs is expected to increase the robustness of the treatment, due to imposing higher complexity on the solution tackling a problem of high complexity, namely the physiology of a tumor. Thus, the utilization of metameric representations in the evolutionary optimization method was examined, along with suitable crossover and mutation operators. An open-source physics-based simulator was utilized, namely PhysiCell, after appropriate modifications, to test the fitness of possible treatments with multiple types of NPs. The possible treatments could be comprised of up to ten types of NPs, simultaneously injected in an area close to the cancerous tumour. Initial results seem to suffer from bloat, namely the best solutions discovered are converging towards the maximum amount of different types of NPs, however, without providing a significant return in fitness when compared with solutions of fewer types of NPs. As the large diversity of NPs will most probably prove to be quite toxic in lab experiments, we opted for methods to reduce the bloat, thus, resolve to therapies with fewer types of NPs. Namely, the bloat control methods studied here were removing types of NPs from the optimization genome as part of the mutation operator and applying parsimony pressure in the replacement operator. By utilizing these techniques, the treatments discovered are composed of fewer types of NPs, while their fitness is not significantly smaller

Optimizing radiation therapy treatments by exploring tumour ecosystem dynamics in-silico

In this contribution, we propose a system-level compartmental population dynamics model of tumour cells that interact with the patient (innate) immune system under the impact of radiation therapy (RT). The resulting in silico - model enables us to analyse the system-level impact of radiation on the tumour ecosystem. The Tumour Control Probability (TCP) was calculated for varying conditions concerning therapy fractionation schemes, radio-sensitivity of tumour sub-clones, tumour population doubling time, repair speed and immunological elimination parameters. The simulations exhibit a therapeutic benefit when applying the initial 3 fractions in an interval of 2 days instead of daily delivered fractions. This effect disappears for fast-growing tumours and in the case of incomplete repair. The results suggest some optimisation potential for combined hyperthermia-radiotherapy. Regarding the sensitivity of the proposed model, cellular repair of radiation-induced damages is a key factor for tumour control. In contrast to this, the radio-sensitivity of immune cells does not influence the TCP as long as the radio-sensitivity is higher than those for tumour cells. The influence of the tumour sub-clone structure is small (if no competition is included). This work demonstrates the usefulness of in silico – modelling for identifying optimisation potentials

Evolutionary computational platform for the automatic discovery of nanocarriers for cancer treatment

We present the EVONANO platform for the evolution of nanomedicines with application to anti-cancer treatments. Our work aims to decrease both the time and cost required to develop nanoparticle designs. EVONANO includes a simulator to grow tumours, extract representative scenarios, and simulate nanoparticle transport through these scenarios in order to predict nanoparticle distribution. The nanoparticle designs are optimised using machine learning to efficiently find the most effective anti-cancer treatments. We demonstrate EVONANO with two examples optimising the properties of nanoparticles and treatment to selectively kill cancer cells over a range of tumour environments. Our platform shows how in silico models that capture both tumour and tissue-scale dynamics can be combined with machine learning to optimise nanomedicine

Computational Modelling of Cancer Systems: From Individual to Collective Cell Behaviour

Debido a su complejidad, el cáncer sigue siendo una de las principales causas de muerte a nivel mundial. La creación de prácticas preventivas adecuadas y terapias innovadoras está limitada por la falta de comprensión de los mecanismos básicos que causan el cáncer. Como tal, se deben desarrollar métodos nuevos y más efectivos que avancen nuestra comprensión del cáncer. En los últimos años, se ha visto un aumento en el uso de modelos computacionales para explicar procesos biológicos que son costosos y difíciles de explorar en entornos experimentales. Estos métodos permiten la traducción de mecanismos biológicos en ecuaciones y suposiciones matemáticas que pueden evaluarse utilizando herramientas informáticas para producir nuevas hipótesis. Además, las tecnologías computacionales se están volviendo más potentes debido a la disponibilidad de datos y la amplia capacidad de procesamiento.El objetivo global de esta tesis es diseñar e implementar modelos computacionales de cáncer, comenzando con comportamientos simples y aislados y progresando hacia fenómenos más complejos. Se abordan tres campos de investigación específicos para lograr este objetivo general: (i) motilidad unicelular, (ii) crecimiento tumoral y (iii) formación de patrones. En el primer objetivo, se presenta un modelo computacional para simular la motilidad celular individual que considera las propiedades mecánicas y químicas del microambiente. Posteriormente, este trabajo fue ampliado para tener en cuenta las interacciones célula-célula y reproducir el crecimiento de estructuras tumorales multicelulares. Por último, todos los eventos biológicos mencionados anteriormente fueron considerados y se añadió la diferenciación celular como el bloque de construcción final de esta tesis para simular la formación de patrones espaciales.Además, esta tesis analiza la relevancia de integrar datos experimentales y métodos computacionales para mejorar la precisión biológica y confirmar los resultados del modelo. En particular, muestra cómo se pueden usar técnicas de calibración y optimización para considerar datos empíricos en el diseño y validación de modelos. Los resultados experimentales cualitativos y cuantitativos, tanto de la literatura como de nuevos experimentos, se reproducen en este artículo para mostrar diferentes enfoques en la integración de datos.En general, esta tesis proporciona un modelo de cómo se pueden utilizar los métodos computacionales para analizar y comprender problemas complejos en la biología del cáncer.Demuestra explícitamente cómo los componentes del modelo pueden representar ciertos aspectos de la biología del cáncer, que pueden mejorarse y reproducirse utilizando datos experimentales. En consecuencia, los comportamientos complejos, como el crecimiento tumoral y la formación de patrones, resultan de la intrincada interacción entre los componentes del modelo.<br /

Examining the efficacy of localised gemcitabine therapy for the treatment of pancreatic cancer using a hybrid agent-based model

The prognosis for pancreatic ductal adenocarcinoma (PDAC) patients has not significantly improved in the past 3 decades, highlighting the need for more effective treatment approaches. Poor patient outcomes and lack of response to therapy can be attributed, in part, to a lack of uptake of perfusion of systemically administered chemotherapeutic drugs into the tumour. Wet-spun alginate fibres loaded with the chemotherapeutic agent gemcitabine have been developed as a potential tool for overcoming the barriers in delivery of systemically administrated drugs to the PDAC tumour microenvironment by delivering high concentrations of drug to the tumour directly over an extended period. While exciting, the practicality, safety, and effectiveness of these devices in a clinical setting requires further investigation. Furthermore, an in-depth assessment of the drug-release rate from these devices needs to be undertaken to determine whether an optimal release profile exists. Using a hybrid computational model (agent-based model and partial differential equation system), we developed a simulation of pancreatic tumour growth and response to treatment with gemcitabine loaded alginate fibres. The model was calibrated using in vitro and in vivo data and simulated using a finite volume method discretisation. We then used the model to compare different intratumoural implantation protocols and gemcitabine-release rates. In our model, the primary driver of pancreatic tumour growth was the rate of tumour cell division. We were able to demonstrate that intratumoural placement of gemcitabine loaded fibres was more effective than peritumoural placement. Additionally, we quantified the efficacy of different release profiles from the implanted fibres that have not yet been tested experimentally. Altogether, the model developed here is a tool that can be used to investigate other drug delivery devices to improve the arsenal of treatments available for PDAC and other difficult-to-treat cancers in the future

High-throughput cancer hypothesis testing with an integrated PhysiCell-EMEWS workflow

Abstract Background Cancer is a complex, multiscale dynamical system, with interactions between tumor cells and non-cancerous host systems. Therapies act on this combined cancer-host system, sometimes with unexpected results. Systematic investigation of mechanistic computational models can augment traditional laboratory and clinical studies, helping identify the factors driving a treatment’s success or failure. However, given the uncertainties regarding the underlying biology, these multiscale computational models can take many potential forms, in addition to encompassing high-dimensional parameter spaces. Therefore, the exploration of these models is computationally challenging. We propose that integrating two existing technologies—one to aid the construction of multiscale agent-based models, the other developed to enhance model exploration and optimization—can provide a computational means for high-throughput hypothesis testing, and eventually, optimization. Results In this paper, we introduce a high throughput computing (HTC) framework that integrates a mechanistic 3-D multicellular simulator (PhysiCell) with an extreme-scale model exploration platform (EMEWS) to investigate high-dimensional parameter spaces. We show early results in applying PhysiCell-EMEWS to 3-D cancer immunotherapy and show insights on therapeutic failure. We describe a generalized PhysiCell-EMEWS workflow for high-throughput cancer hypothesis testing, where hundreds or thousands of mechanistic simulations are compared against data-driven error metrics to perform hypothesis optimization. Conclusions While key notational and computational challenges remain, mechanistic agent-based models and high-throughput model exploration environments can be combined to systematically and rapidly explore key problems in cancer. These high-throughput computational experiments can improve our understanding of the underlying biology, drive future experiments, and ultimately inform clinical practice

Using (and reusing)experimental data in computational models [Presented at ASCB 2017]

These are my slides from my talk on multicellular data sharing and use in computational biology. <div><br></div><div>I presented these at the ASCB-EMBO 2017 Meeting in Subgroup S (http://ascb-embo2017.ascb.org/subgroup-s/)</div><div><br></div><div><br></div><div>{High-throughput cancer hypothesis testing with an integrated PhysiCell-EMEWS workflow [Presented at SuperComputing SC17, Computational Approaches for Cancer Workshop]<br></div