Exploiting clinical trial data drastically narrows the window of possible solutions to the problem of clinical adaptation of a multiscale cancer model

A Multidisciplinary Hyper-Modeling Scheme in Personalized In Silico Oncology: Coupling Cell Kinetics with Metabolism, Signaling Networks, and Biomechanics as Plug-In Component Models of a Cancer Digital Twin

The massive amount of human biological, imaging, and clinical data produced by multiple and diverse sources necessitates integrative modeling approaches able to summarize all this information into answers to specific clinical questions. In this paper, we present a hypermodeling scheme able to combine models of diverse cancer aspects regardless of their underlying method or scale. Describing tissue-scale cancer cell proliferation, biomechanical tumor growth, nutrient transport, genomic-scale aberrant cancer cell metabolism, and cell-signaling pathways that regulate the cellular response to therapy, the hypermodel integrates mutation, miRNA expression, imaging, and clinical data. The constituting hypomodels, as well as their orchestration and links, are described. Two specific cancer types, Wilms tumor (nephroblastoma) and non-small cell lung cancer, are addressed as proof-of-concept study cases. Personalized simulations of the actual anatomy of a patient have been conducted. The hypermodel has also been applied to predict tumor control after radiotherapy and the relationship between tumor proliferative activity and response to neoadjuvant chemotherapy. Our innovative hypermodel holds promise as a digital twin-based clinical decision support system and as the core of future in silico trial platforms, although additional retrospective adaptation and validation are necessary.

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.

Studying the regression profiles of cervical tumours during radiotherapy treatment using a patient-specific multiscale model

Apart from offering insight into the biomechanisms involved in cancer, many recent mathematical modeling efforts aspire to the ultimate goal of clinical translation, wherein models are designed to be used in the future as clinical decision support systems in the patient-individualized context. Most significant challenges are the integration of multiscale biodata and the patient-specific model parameterization. A central aim of this study was the design of a clinically-relevant parameterization methodology for a patient-specific computational model of cervical cancer response to radiotherapy treatment with concomitant cisplatin, built around a tumour features-based search of the parameter space. Additionally, a methodological framework for the predictive use of the model was designed, including a scoring method to quantitatively reflect the similarity and bilateral predictive ability of any two tumours in terms of their regression profile. The methodology was applied to the datasets of eight patients. Tumour scenarios in accordance with the available longitudinal data have been determined. Predictive investigations identified three patient cases, anyone of which can be used to predict the volumetric evolution throughout therapy of the tumours of the other two with very good results. Our observations show that the presented approach is promising in quantifiably differentiating tumours with distinct regression profiles.

SpheroidSim-Preliminary evaluation of a new computational tool to predict the influence of cell cycle time and phase fraction on spheroid growth

BACKGROUND

There is a relative paucity of research that integrates materials science and bioengineering with computational simulations to decipher the intricate processes promoting cancer progression. Therefore, a first-generation computational model, SpheroidSim, was developed that includes a biological data set derived from a bioengineered spheroid model to obtain a quantitative description of cell kinetics.

RESULTS

SpheroidSim is a 3D agent-based model simulating the growth of multicellular cancer spheroids. Cell cycle time and phases mathematically motivated the population growth. SpheroidSim simulated the growth dynamics of multiple spheroids by individually defining a collection of specific phenotypic traits and characteristics for each cell. Experimental data derived from a hydrogel-based spheroid model were fit to the predictions providing insight into the influence of cell cycle time (CCT) and cell phase fraction (CPF) on the cell population. A comparison of the number of active cells predicted for each analysis showed that the value and method used to define CCT had a greater effect on the predicted cell population than CPF. The model predictions were similar to the experimental results for the number of cells, with the predicted total number of cells varying by 8% and 12%, respectively, compared to the experimental data.

CONCLUSIONS

SpheroidSim is a first step in developing a biologically based predictive tool capable of revealing fundamental elements in cancer cell physiology. This computational model may be applied to study the effect of the microenvironment on spheroid growth and other cancer cell types that demonstrate a similar multicellular clustering behavior as the population develops. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:1335-1343, 2018.

A Numerical Handling of the Boundary Conditions Imposed by the Skull on an Inhomogeneous Diffusion-Reaction Model of Glioblastoma Invasion Into the Brain: Clinical Validation Aspects

A novel explicit triscale reaction-diffusion numerical model of glioblastoma multiforme tumor growth is presented. The model incorporates the handling of Neumann boundary conditions imposed by the cranium and takes into account both the inhomogeneous nature of human brain and the complexity of the skull geometry. The finite-difference time-domain method is adopted. To demonstrate the workflow of a possible clinical validation procedure, a clinical case/scenario is addressed. A good agreement of the in silico calculated value of the doubling time (ie, the time for tumor volume to double) with the value of the same quantity based on tomographic imaging data has been observed. A theoretical exploration suggests that a rough but still quite informative value of the doubling time may be calculated based on a homogeneous brain model. The model could serve as the main component of a continuous mathematics-based glioblastoma oncosimulator aiming at supporting the clinician in the optimal patient-individualized design of treatment using the patient's multiscale data and experimenting in silico (ie, on the computer).

A Modular Repository-based Infrastructure for Simulation Model Storage and Execution Support in the Context of In Silico Oncology and In Silico Medicine

The plethora of available disease prediction models and the ongoing process of their application into clinical practice - following their clinical validation - have created new needs regarding their efficient handling and exploitation. Consolidation of software implementations, descriptive information, and supportive tools in a single place, offering persistent storage as well as proper management of execution results, is a priority, especially with respect to the needs of large healthcare providers. At the same time, modelers should be able to access these storage facilities under special rights, in order to upgrade and maintain their work. In addition, the end users should be provided with all the necessary interfaces for model execution and effortless result retrieval. We therefore propose a software infrastructure, based on a tool, model and data repository that handles the storage of models and pertinent execution-related data, along with functionalities for execution management, communication with third-party applications, user-friendly interfaces to access and use the infrastructure with minimal effort and basic security features.

Workflow-driven clinical decision support for personalized oncology

BACKGROUND

The adoption in oncology of Clinical Decision Support (CDS) may help clinical users to efficiently deal with the high complexity of the domain, lead to improved patient outcomes, and reduce the current knowledge gap between clinical research and practice. While significant effort has been invested in the implementation of CDS, the uptake in the clinic has been limited. The barriers to adoption have been extensively discussed in the literature. In oncology, current CDS solutions are not able to support the complex decisions required for stratification and personalized treatment of patients and to keep up with the high rate of change in therapeutic options and knowledge.

RESULTS

To address these challenges, we propose a framework enabling efficient implementation of meaningful CDS that incorporates a large variety of clinical knowledge models to bring to the clinic comprehensive solutions leveraging the latest domain knowledge. We use both literature-based models and models built within the p-medicine project using the rich datasets from clinical trials and care provided by the clinical partners. The framework is open to the biomedical community, enabling reuse of deployed models by third-party CDS implementations and supporting collaboration among modelers, CDS implementers, biomedical researchers and clinicians. To increase adoption and cope with the complexity of patient management in oncology, we also support and leverage the clinical processes adhered to by healthcare organizations. We design an architecture that extends the CDS framework with workflow functionality. The clinical models are embedded in the workflow models and executed at the right time, when and where the recommendations are needed in the clinical process.

CONCLUSIONS

In this paper we present our CDS framework developed in p-medicine and the CDS implementation leveraging the framework. To support complex decisions, the framework relies on clinical models that encapsulate relevant clinical knowledge. Next to assisting the decisions, this solution supports by default (through modeling and implementation of workflows) the decision processes as well and exploits the knowledge embedded in those processes.

Numerical simulation of vascular tumour growth under antiangiogenic treatment: addressing the paradigm of single-agent bevacizumab therapy with the use of experimental data

Background

Antiangiogenic agents have been recently added to the oncological armamentarium with bevacizumab probably being the most popular representative in current clinical practice. The elucidation of the mode of action of these agents is a prerequisite for personalized prediction of antiangiogenic treatment response and selection of patients who may benefit from this kind of therapy. To this end, having used as a basis a preexisting continuous vascular tumour growth model which addresses the targeted nature of antiangiogenic treatment, we present a paper characterized by the following three features. First, the integration of a two-compartmental bevacizumab specific pharmacokinetic module into the core of the aforementioned preexisting model. Second, its mathematical modification in order to reproduce the asymptotic behaviour of tumour volume in the theoretical case of a total destruction of tumour neovasculature. Third, the exploitation of a range of published animal datasets pertaining to antitumour efficacy of bevacizumab on various tumour types (breast, lung, head and neck, colon).

Results

Results for both the unperturbed growth and the treatment module reveal qualitative similarities with experimental observations establishing the biologically acceptable behaviour of the model. The dynamics of the untreated tumour has been studied via a parameter analysis, revealing the role of each relevant input parameter to tumour evolution. The combined effect of endogenous proangiogenic and antiangiogenic factors on the angiogenic potential of a tumour is also studied, in order to capture the dynamics of molecular competition between the two key-players of tumoural angiogenesis. The adopted methodology also allows accounting for the newly recognized direct antitumour effect of the specific agent.

Conclusions

Interesting observations have been made, suggesting a potential size-dependent tumour response to different treatment modalities and determining the relative timing of cytotoxic versus antiangiogenic agents administration. Insight into the comparative effectiveness of different antiangiogenic treatment strategies is revealed. The results of a series of in vivo experiments in mice bearing diverse types of tumours (breast, lung, head and neck, colon) and treated with bevacizumab are successfully reproduced, supporting thus the validity of the underlying model.

Reviewers

This article was reviewed by L. Hanin, T. Radivoyevitch and L. Edler.

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.

In silico oncology: exploiting clinical studies to clinically adapt and validate multiscale oncosimulators

This paper presents a brief outline of the notion and the system of oncosimulator in conjunction with a high level description of the basics of its core multiscale model simulating clinical tumor response to treatment. The exemplary case of lung cancer preoperatively treated with a combination of chemotherapeutic agents is considered. The core oncosimulator model is based on a primarily top-down, discrete entity - discrete event multiscale simulation approach. The critical process of clinical adaptation of the model by exploiting sets of multiscale data originating from clinical studies/trials is also outlined. Concrete clinical adaptation results are presented. The adaptation process also conveys important aspects of the planned clinical validation procedure since the same type of multiscale data - although not the same data itself- is to be used for clinical validation. By having exploited actual clinical data in conjunction with plausible literature-based values of certain model parameters, a realistic tumor dynamics behavior has been demonstrated. The latter supports the potential of the specific oncosimulator to serve as a personalized treatment optimizer following an eventually successful completion of the clinical adaptation and validation process.

Computational nanomedicine: modeling of nanoparticle-mediated hyperthermal cancer therapy

Nanoparticle-mediated hyperthermia for cancer therapy is a growing area of cancer nanomedicine because of the potential for localized and targeted destruction of cancer cells. Localized hyperthermal effects are dependent on many factors, including nanoparticle size and shape, excitation wavelength and power, and tissue properties. Computational modeling is an important tool for investigating and optimizing these parameters. In this review, we focus on computational modeling of magnetic and gold nanoparticle-mediated hyperthermia, followed by a discussion of new opportunities and challenges.

Towards in silico oncology: adapting a four dimensional nephroblastoma treatment model to a clinical trial case based on multi-method sensitivity analysis

In the past decades a great progress in cancer research has been made although medical treatment is still widely based on empirically established protocols which have many limitations. Computational models address such limitations by providing insight into the complex biological mechanisms of tumor progression. A set of clinically-oriented, multiscale models of solid tumor dynamics has been developed by the In Silico Oncology Group (ISOG), Institute of Communication and Computer Systems (ICCS)-National Technical University of Athens (NTUA) to study cancer growth and response to treatment. Within this context using certain representative parameter values, tumor growth and response have been modeled under a cancer preoperative chemotherapy protocol in the framework of the SIOP 2001/GPOH clinical trial. A thorough cross-method sensitivity analysis of the model has been performed. Based on the sensitivity analysis results, a reasonable adaptation of the values of the model parameters to a real clinical case of bilateral nephroblastomatosis has been achieved. The analysis presented supports the potential of the model for the study and eventually the future design of personalized treatment schemes and/or schedules using the data obtained from in vitro experiments and clinical studies.

In silico modelling of treatment-induced tumour cell kill: developments and advances

Mathematical and stochastic computer (in silico) models of tumour growth and treatment response of the past and current eras are presented, outlining the aims of the models, model methodology, the key parameters used to describe the tumour system, and treatment modality applied, as well as reported outcomes from simulations. Fractionated radiotherapy, chemotherapy, and combined therapies are reviewed, providing a comprehensive overview of the modelling literature for current modellers and radiobiologists to ignite the interest of other computational scientists and health professionals of the ever evolving and clinically relevant field of tumour modelling.

TumorML: Concept and requirements of an in silico cancer modelling markup language

This paper describes the initial groundwork carried out as part of the European Commission funded Transatlantic Tumor Model Repositories project, to develop a new markup language for computational cancer modelling, TumorML. In this paper we describe the motivations for such a language, arguing that current state-of-the-art biomodelling languages are not suited to the cancer modelling domain. We go on to describe the work that needs to be done to develop TumorML, the conceptual design, and a description of what existing markup languages will be used to compose the language specification.

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

The challenge of modelling cancer presents a major opportunity to improve our ability to reduce mortality from malignant neoplasms, improve treatments and meet the demands associated with the individualization of care needs. This is the central motivation behind the ContraCancrum project. By developing integrated multi-scale cancer models, ContraCancrum is expected to contribute to the advancement of in silico oncology through the optimization of cancer treatment in the patient-individualized context by simulating the response to various therapeutic regimens. The aim of the present paper is to describe a novel paradigm for designing clinically driven multi-scale cancer modelling by bringing together basic science and information technology modules. In addition, the integration of the multi-scale tumour modelling components has led to novel concepts of personalized clinical decision support in the context of predictive oncology, as is also discussed in the paper. Since clinical adaptation is an inelastic prerequisite, a long-term clinical adaptation procedure of the models has been initiated for two tumour types, namely non-small cell lung cancer and glioblastoma multiforme; its current status is briefly summarized.