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

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.

Theoretical evaluation of the impact of diverse treatment conditions by calculation of the tumor control probability (TCP) of simulated cervical cancer Hyperthermia-Radiotherapy (HT-RT) treatments in-silico

INTRODUCTION

Hyperthermia (HT) induces various cellular biological processes, such as repair impairment and direct HT cell killing. In this context, in-silico biophysical models that translate deviations in the treatment conditions into clinical outcome variations may be used to study the extent of such processes and their influence on combined hyperthermia plus radiotherapy (HT + RT) treatments under varying conditions.

METHODS

An extended linear-quadratic model calibrated for SiHa and HeLa cell lines (cervical cancer) was used to theoretically study the impact of varying HT treatment conditions on radiosensitization and direct HT cell killing effect. Simulated patients were generated to compute the Tumor Control Probability (TCP) under different HT conditions (number of HT sessions, temperature and time interval), which were randomly selected within margins based on reported patient data.

RESULTS

Under the studied conditions, model-based simulations suggested a treatment improvement with a total CEM43 thermal dose of approximately 10 min. Additionally, for a given thermal dose, TCP increased with the number of HT sessions. Furthermore, in the simulations, we showed that the TCP dependence on the temperature/time interval is more correlated with the mean value than with the minimum/maximum value and that comparing the treatment outcome with the mean temperature can be an excellent strategy for studying the time interval effect.

CONCLUSION

The use of thermoradiobiological models allows us to theoretically study the impact of varying thermal conditions on HT + RT treatment outcomes. This approach can be used to optimize HT treatments, design clinical trials, and interpret patient data.

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.

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.

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.

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.

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.

Adaptive Mathematical Model of Tumor Response to Radiotherapy Based on CBCT Data

Mathematical modeling of tumor response to radiotherapy has the potential of enhancing the quality of the treatment plan, which can be even tailored on an individual basis. Lack of extensive in vivo validation has prevented, however, reliable clinical translation of modeling outcomes. Image-guided radiotherapy is a consolidated treatment modality based on computed tomographic (CT) imaging for tumor delineation and volumetric cone beam CT data for periodic checks during treatment. In this study, a macroscopic model of tumor growth and radiation response is proposed, being able to adapt along the treatment course as volumetric tumor data become available. Model parameter learning was based on cone beam CT images in 13 uterine cervical cancer patients, subdivided into three groups (G1, G2, G3) according to tumor type and treatment. Three group-specific parameter sets (PS1, PS2, and PS3) on one general parameter set (PSa) were applied. The corresponding average model fitting errors were 14%, 18%, 13%, and 21%, respectively. The model adaptation testing was performed using volume data of three patients, other than the ones involved in the parameter learning. The extrapolation performance of the general model was improved, while comparable prediction errors were found for the group-specific approach. This suggests that an online parameter tuning can overcome the limitations of a suboptimal patient stratification, which appeared otherwise a critical issue.

Strategy for stochastic dose-rate induced enhanced elimination of malignant tumour without dose escalation

The efficacy of radiation therapy, a primary modality of cancer treatment, depends in general upon the total radiation dose administered to the tumour during the course of therapy. Nevertheless, the delivered radiation also irradiates normal tissues and dose escalation procedure often increases the elimination of normal tissue as well. In this article, we have developed theoretical frameworks under the premise of linear-quadratic-linear (LQL) model using stochastic differential equation and Jensen's inequality for exploring the possibility of attending to the two therapeutic performance objectives in contraposition-increasing the elimination of prostate tumour cells and enhancing the relative sparing of normal tissue in fractionated radiation therapy, within a prescribed limit of total radiation dose. Our study predicts that stochastic temporal modulation in radiation dose-rate appreciably enhances prostate tumour cell elimination, without needing dose escalation in radiation therapy. However, constant higher dose-rate can also enhance the elimination of tumour cells. In this context, we have shown that the sparing of normal tissue with stochastic dose-rate is considerably more than the sparing of normal tissue with the equivalent constant higher dose-rate. Further, by contrasting the stochastic dose-rate effects under LQL and linear-quadratic (LQ) models, we have also shown that the LQ model over-estimates stochastic dose-rate effect in tumour and under-estimates the stochastic dose-rate effect in normal tissue. Our study indicates the possibility of utilizing stochastic modulation of radiation dose-rate for designing enhanced radiation therapy protocol for cancer.

Dendritic cell vaccination for glioblastoma multiforme: review with focus on predictive factors for treatment response

Glioblastoma multiforme (GBM) is the most common and most aggressive type of primary brain cancer. Since median overall survival with multimodal standard therapy is only 15 months, there is a clear need for additional effective and long-lasting treatments. Dendritic cell (DC) vaccination is an experimental immunotherapy being tested in several Phase I and Phase II clinical trials. In these trials, safety and feasibility have been proven, and promising clinical results have been reported. On the other hand, it is becoming clear that not every GBM patient will benefit from this highly personalized treatment. Defining the subgroup of patients likely to respond to DC vaccination will position this option correctly amongst other new GBM treatment modalities, and pave the way to incorporation in standard therapy. This review provides an overview of GBM treatment options and focuses on the currently known prognostic and predictive factors for response to DC vaccination. In this way, it will provide the clinician with the theoretical background to refer patients who might benefit from this treatment.