In recent years, cancer genome sequencing and other high-throughput studies of cancer genomes have generated many notable discoveries. In this review, Novel genomic alteration mechanisms, such as chromothripsis (chromosomal crisis) and kataegis (mutation storms), and their implications for cancer are discussed. Genomic alterations spur cancer genome evolution. Thus, the relationship between cancer clonal evolution and cancer stems cells is commented. The key question in cancer biology concerns how these genomic alterations support cancer development and metastasis in the context of biological functioning. Thus far, efforts such as pathway analysis have improved the understanding of the functional contributions of genetic mutations and DNA copy number variations to cancer development, progression and metastasis. However, the known pathways correspond to a small fraction, plausibly 5-10%, of somatic mutations and genes with an altered copy number. To develop a comprehensive understanding of the function of these genomic alterations in cancer, an integrative network framework is proposed and discussed. Finally, the challenges and the directions of studying cancer omic data using an in
In this paper, we explore the application of ensemble optimal control to derive enhanced strategies for pharmacological cancer treatment, and we tackle the problem of the long-term management of the disease, i.e., when the complete eradication of the tumor is not achievable. In particular, we focus on moving beyond the classical clinical approach of giving the patient the maximal tolerated drug dose (MTD), which does not properly exploit the fight among sensitive and resistant cells for the available resources. Here, we employ a Lotka-Volterra model to describe the competing subpopulations, and we enclose this system within the ensemble control framework. In the first part, we establish general results suitable for application to various cancers. Then, we carry out numerical simulations in the setting of prostate cancer treated with androgen deprivation therapy, yielding a computed policy that is reminiscent of the medical `active surveillance' paradigm. Finally, inspired by the numerical evidence, we propose a variant of the celebrated adaptive therapy (AT), which we call `Off-On' AT.
We investigate a new dynamical system that describes tumor-host interaction. The equation that describes the untreated tumor growth is based on non-extensive statistical mechanics. Recently, this model has been shown to fit successfully exponential, Gompertz, logistic, and power-law tumor growths. We have been able to include as many hallmarks of cancer as possible. We study also the dynamic response of cancer under therapy. Using our model, we can make predictions about the different outcomes when we change the parameters, and/or the initial conditions. We can determine the importance of different factors to influence tumor growth. We discover synergistic therapeutic effects of different treatments and drugs. Cancer is generally untreatable using conventional monotherapy. We consider conventional therapies, oncogene-targeted therapies, tumor-suppressors gene-targeted therapies, immunotherapies, anti-angiogenesis therapies, virotherapy, among others. We need therapies with the potential to target both tumor cells and the tumors' microenvironment. Drugs that target oncogenes and tumor-suppressor genes can be effective in the treatment of some cancers. However, most tumors do reoccur
Cancer development is driven by mutations and selective forces, including the action of the immune system and interspecific competition. When administered to patients, anti-cancer therapies affect the development and dynamics of tumours, possibly with various degrees of resistance due to immunoediting and microenvironment. Tumours are able to express a variety of competing phenotypes with different attributes and thus respond differently to various anti-cancer therapies. In this paper, a mathematical framework incorporating a system of delay differential equations for the immune system activation cycle and an agent-based approach for tumour-immune interaction is presented. The focus is on those metastatic, secondary solid lesions that are still undetected and non-vascularised. By using available experimental data, we analyse the effects of combination therapies on these lesions and investigate the role of mutations on the rates of success of common treatments. Findings show that mutations, growth properties and immunoediting influence therapies' outcomes in nonlinear and complex ways, affecting cancer lesion morphologies, phenotypical compositions and overall proliferation patterns
Embryology has long played a foundational role in shaping our scientific understanding of animal evolution. In recent decades, growing evidence has also highlighted its role in cancer. Despite the indisputable similarities between embryonic development and cancer, there has been limited discussion on the profound embryological implications for the disease. This article explores the understanding of cancer as an embryological and evolutionary phenomenon, offering a fresh perspective on the disease and discussing immediate consequences in the search for therapeutic approaches
The tumor-immune system plays a critical role in colorectal cancer progression. Recent preclinical and clinical studies showed that combination therapy with anti-PD-L1 and cancer vaccines improved treatment response. In this study, we developed a multiscale mathematical model of interactions among tumors, immune cells, and cytokines to investigate tumor evolutionary dynamics under different therapeutic strategies. Additionally, we established a computational framework based on approximate Bayesian computation to generate virtual tumor samples and capture inter-individual heterogeneity in treatment response. The results demonstrated that a multiple low-dose regimen significantly reduced advanced tumor burden compared to baseline treatment in anti-PD-L1 therapy. In contrast, the maximum dose therapy yielded superior tumor growth control in cancer vaccine therapy. Furthermore, cytotoxic T cells were identified as a consistent predictive biomarker both before and after treatment initiation. Notably, the cytotoxic T cells-to-regulatory T cells ratio specifically served as a robust pre-treatment predictive biomarker, offering potential clinical utility for patient stratification and ther
In this paper, a new stochastic framework to determine optimal combination therapies in colon cancer-induced immune response is presented. The dynamics of colon cancer is described through an Itö stochastic process, whose probability density function evolution is governed by the Fokker-Planck equation. An open-loop control optimization problem is proposed to determine the optimal combination therapies. Numerical results with combination therapies comprising of the chemotherapy drug \ind{Doxorubicin} and immunotherapy drug IL-2 validate the proposed framework.
Rapid technological advances have allowed for molecular profiling across multiple omics domains from a single sample for clinical decision making in many diseases, especially cancer. As tumor development and progression are dynamic biological processes involving composite genomic aberrations, key challenges are to effectively assimilate information from these domains to identify genomic signatures and biological entities that are druggable, develop accurate risk prediction profiles for future patients, and identify novel patient subgroups for tailored therapy and monitoring. We propose integrative probabilistic frameworks for high-dimensional multiple-domain cancer data that coherently incorporate dependence within and between domains to accurately detect tumor subtypes, thus providing a catalogue of genomic aberrations associated with cancer taxonomy. We propose an innovative, flexible and scalable Bayesian nonparametric framework for simultaneous clustering of both tumor samples and genomic probes. We describe an efficient variable selection procedure to identify relevant genomic aberrations that can potentially reveal underlying drivers of a disease. Although the work is motivat
In a previous work, we presented a model that integrates cancer cell differentiation and immunotherapy, analysing a particular therapy against cancer stem cells by cytotoxic cell vaccines. As every biological system is exposed to random fluctuations, it is important to study its stochasticity. The influence of demographic and multiplicative noise in the system is carry out on the parameters of reproduction and death in cancer cells. On the other hand, we incorporated fluctuations by adding multiplicative noise. In both cases, we analysed the dynamics for different values of the parameters involved. The final amount of cancer cells decreases for different combinations of these parameters and noise intensity is found.
Evolutionary therapy (ET) aims to steer tumor evolution by adjusting treatment timing and dosing to control rather than eradicate tumor burden. Clinical use requires reliable monitoring of tumor dynamics to inform mathematical models that guide therapy. In cancers such as metastatic castrate-resistant prostate cancer and relapsed platinum-sensitive ovarian cancer, ET models are informed by serial serum biomarkers. For cancers lacking reliable biomarkers, such as metastatic non-small cell lung cancer (NSCLC), radiographic imaging remains the primary method for treatment response assessment, typically using RECIST 1.1 criteria. RECIST, which tracks a few lesions with one-dimensional (1D) measurements and defines progression relative to the nadir, the smallest tumor burden recorded after treatment, was not designed to support ET. It may miss early regrowth, underrepresent tumor burden, and obscure disease trends. Using a virtual NSCLC patient model, we demonstrate that lesion selection and measurement dimensionality strongly affect progression detection. Two-dimensional metrics provide modest improvement, but only 3D volumetric measurements accurately capture both tumor burden and its
In this paper, a new framework for obtaining personalized optimal treatment strategies in colon cancer-induced angiogenesis is presented. The dynamics of colon cancer is given by a Itó stochastic process, which helps in modeling the randomness present in the system. The stochastic dynamics is then represented by the Fokker-Planck (FP) partial differential equation (PDE) that governs the evolution of the associated probability density function. The optimal therapies are obtained using a three step procedure. First, a finite dimensional FP-constrained optimization problem is formulated that takes input individual noisy patient data, and is solved to obtain the unknown parameters corresponding to the individual tumor characteristics. Next, a sensitivity analysis of the optimal parameter set is used to determine the parameters to be controlled, thus, helping in assessing the types of treatment therapies. Finally, a feedback FP control problem is solved to determine the optimal combination therapies. Numerical results with the combination drug, comprising of Bevacizumab and Capecitabine, demonstrate the efficiency of the proposed framework.
Lung cancer is the primary cause of cancer-related mortality, claiming approximately 1.79 million lives globally in 2020, with an estimated 2.21 million new cases diagnosed within the same period. Among these, Non-Small Cell Lung Cancer (NSCLC) is the predominant subtype, characterized by a notably bleak prognosis and low overall survival rate of approximately 25% over five years across all disease stages. However, survival outcomes vary considerably based on the stage at diagnosis and the therapeutic interventions administered. Recent advancements in artificial intelligence (AI) have revolutionized the landscape of lung cancer prognosis. AI-driven methodologies, including machine learning and deep learning algorithms, have shown promise in enhancing survival prediction accuracy by efficiently analyzing complex multi-omics data and integrating diverse clinical variables. By leveraging AI techniques, clinicians can harness comprehensive prognostic insights to tailor personalized treatment strategies, ultimately improving patient outcomes in NSCLC. Overviewing AI-driven data processing can significantly help bolster the understanding and provide better directions for using such syste
Cancer is increasingly perceived as a systems-level, network phenomenon. The major trend of malignant transformation can be described as a two-phase process, where an initial increase of network plasticity is followed by a decrease of plasticity at late stages of tumor development. The fluctuating intensity of stress factors, like hypoxia, inflammation and the either cooperative or hostile interactions of tumor inter-cellular networks, all increase the adaptation potential of cancer cells. This may lead to the bypass of cellular senescence, and to the development of cancer stem cells. We propose that the central tenet of cancer stem cell definition lies exactly in the indefinability of cancer stem cells. Actual properties of cancer stem cells depend on the individual "stress-history" of the given tumor. Cancer stem cells are characterized by an extremely large evolvability (i.e. a capacity to generate heritable phenotypic variation), which corresponds well with the defining hallmarks of cancer stem cells: the possession of the capacity to self-renew and to repeatedly re-build the heterogeneous lineages of cancer cells that comprise a tumor in new environments. Cancer stem cells rep
We present a general computational theory of cancer and its developmental dynamics. The theory is based on a theory of the architecture and function of developmental control networks which guide the formation of multicellular organisms. Cancer networks are special cases of developmental control networks. Cancer results from transformations of normal developmental networks. Our theory generates a natural classification of all possible cancers based on their network architecture. Each cancer network has a unique topology and semantics and developmental dynamics that result in distinct clinical tumor phenotypes. We apply this new theory with a series of proof of concept cases for all the basic cancer types. These cases have been computationally modeled, their behavior simulated and mathematically described using a multicellular systems biology approach. There are fascinating correspondences between the dynamic developmental phenotype of computationally modeled {\em in silico} cancers and natural {\em in vivo} cancers. The theory lays the foundation for a new research paradigm for understanding and investigating cancer. The theory of cancer networks implies that new diagnostic methods
There is a widening recognition that cancer cells are products of complex developmental processes. Carcinogenesis and metastasis formation are increasingly described as systems-level, network phenomena. Here we propose that malignant transformation is a two-phase process, where an initial increase of system plasticity is followed by a decrease of plasticity at late stages of carcinogenesis as a model of cellular learning. We describe the hallmarks of increased system plasticity of early, tumor initiating cells, such as increased noise, entropy, conformational and phenotypic plasticity, physical deformability, cell heterogeneity and network rearrangements. Finally, we argue that the large structural changes of molecular networks during cancer development necessitate a rather different targeting strategy in early and late phase of carcinogenesis. Plastic networks of early phase cancer development need a central hit, while rigid networks of late stage primary tumors or established metastases should be attacked by the network influence strategy, such as by edgetic, multi-target, or allo-network drugs. Cancer stem cells need special diagnosis and targeting, since their dormant and rapid
This paper proposes a general framework for probabilistic certification of cancer therapies. The certification is defined in terms of two key issues which are the tumor contraction and the lower admissible bound on the circulating lymphocytes which is viewed as indicator of the patient health. The certification is viewed as the ability to guarantee with a predefined high probability the success of the therapy over a finite horizon despite of the unavoidable high uncertainties affecting the dynamic model that is used to compute the optimal scheduling of drugs injection. The certification paradigm can be viewed as a tool for tuning the treatment parameters and protocols as well as for getting a rational use of limited or expensive drugs. The proposed framework is illustrated using the specific problem of combined immunotherapy/chemotherapy of cancer.
Prostate cancer can be lethal in advanced stages, for which chemotherapy may become the only viable therapeutic option. While there is no clear clinical management strategy fitting all patients, cytotoxic chemotherapy with docetaxel is currently regarded as the gold standard. However, tumors may regain activity after treatment conclusion and become resistant to docetaxel. This situation calls for new delivery strategies and drug compounds enabling an improved therapeutic outcome. Combination of docetaxel with antiangiogenic therapy has been considered a promising strategy. Bevacizumab is the most common antiangiogenic drug, but clinical studies have not revealed a clear benefit from its combination with docetaxel. Here, we capitalize on our prior work on mathematical modeling of prostate cancer growth subjected to combined cytotoxic and antiangiogenic therapies, and propose an optimal control framework to robustly compute the drug-independent cytotoxic and antiangiogenic effects enabling an optimal therapeutic control of tumor dynamics. We describe the formulation of the optimal control problem, for which we prove the existence of at least a solution and determine the necessary fir
With the advances in artificial intelligence (AI), data-driven algorithms are becoming increasingly popular in the medical domain. However, due to the nonlinear and complex behavior of many of these algorithms, decision-making by such algorithms is not trustworthy for clinicians and is considered a black-box process. Hence, the scientific community has introduced explainable artificial intelligence (XAI) to remedy the problem. This systematic scoping review investigates the application of XAI in breast cancer detection and risk prediction. We conducted a comprehensive search on Scopus, IEEE Explore, PubMed, and Google Scholar (first 50 citations) using a systematic search strategy. The search spanned from January 2017 to July 2023, focusing on peer-reviewed studies implementing XAI methods in breast cancer datasets. Thirty studies met our inclusion criteria and were included in the analysis. The results revealed that SHapley Additive exPlanations (SHAP) is the top model-agnostic XAI technique in breast cancer research in terms of usage, explaining the model prediction results, diagnosis and classification of biomarkers, and prognosis and survival analysis. Additionally, the SHAP mo
Background Precise prediction of cancer types is vital for cancer diagnosis and therapy. Important cancer marker genes can be inferred through predictive model. Several studies have attempted to build machine learning models for this task however none has taken into consideration the effects of tissue of origin that can potentially bias the identification of cancer markers. Results In this paper, we introduced several Convolutional Neural Network (CNN) models that take unstructured gene expression inputs to classify tumor and non-tumor samples into their designated cancer types or as normal. Based on different designs of gene embeddings and convolution schemes, we implemented three CNN models: 1D-CNN, 2D-Vanilla-CNN, and 2D-Hybrid-CNN. The models were trained and tested on combined 10,340 samples of 33 cancer types and 731 matched normal tissues of The Cancer Genome Atlas (TCGA). Our models achieved excellent prediction accuracies (93.9-95.0%) among 34 classes (33 cancers and normal). Furthermore, we interpreted one of the models, known as 1D-CNN model, with a guided saliency technique and identified a total of 2,090 cancer markers (108 per class). The concordance of differential e
Nanorobots are a promising development in targeted drug delivery and the treatment of neurological disorders, with potential for crossing the blood-brain barrier (BBB). These small devices leverage advancements in nanotechnology and bioengineering for precise navigation and targeted payload delivery, particularly for conditions like brain tumors, Alzheimer's disease, and Parkinson's disease. Recent progress in artificial intelligence (AI) and machine learning (ML) has improved the navigation and effectiveness of nanorobots, allowing them to detect and interact with cancer cells through biomarker analysis. This study presents a new reinforcement learning (RL) framework for optimizing nanorobot navigation in complex biological environments, focusing on cancer cell detection by analyzing the concentration gradients of surrounding biomarkers. We utilize a computer simulation model to explore the behavior of nanorobots in a three-dimensional space with cancer cells and biological barriers. The proposed method uses Q-learning to refine movement strategies based on real-time biomarker concentration data, enabling nanorobots to autonomously navigate to cancerous tissues for targeted drug d