Whether the nom de guerre is Mathematical Oncology, Computational or Systems Biology, Theoretical Biology, Evolutionary Oncology, Bioinformatics, or simply Basic Science, there is no denying that mathematics continues to play an increasingly prominent role in cancer research. Mathematical Oncology—defined here simply as the use of mathematics in cancer research—complements and overlaps with a number of other fields that rely on mathematics as a core methodology. As a result, Mathematical Oncology has a broad scope, ranging from theoretical studies to clinical trials designed with mathematical models. This Roadmap differentiates Mathematical Oncology from related fields and demonstrates specific areas of focus within this unique field of research. The dominant theme of this Roadmap is the personalization of medicine through mathematics, modelling, and simulation. This is achieved through the use of patient-specific clinical data to: develop individualized screening strategies to detect cancer earlier; make predictions of response to therapy; design adaptive, patient-specific treatment plans to overcome therapy resistance; and establish domain-specific standards to share model predictions and to make models and simulations reproducible. The cover art for this Roadmap was chosen as an apt metaphor for the beautiful, strange, and evolving relationship between mathematics and cancer.

Our nationwide network of BME women faculty collectively argue that racial funding disparity by the National Institutes of Health (NIH) remains the most insidious barrier to success of Black faculty in our profession. We thus refocus attention on this critical barrier and suggest solutions on how it can be dismantled.

Blood travels throughout the body in an extensive network of vessels – arteries, veins and capillaries. This vascular network is not static, but instead dynamically remodels in response to stimuli from cells in the nearby tissue. In particular, the smallest vessels – arterioles, venules and capillaries – can be extended, expanded or pruned, in response to exercise, ischaemic events, pharmacological interventions, or other physiological and pathophysiological events. In this review, we describe the multi-step morphogenic process of angiogenesis – the sprouting of new blood vessels – and the stability of vascular networks in vivo. In particular, we review the known interactions between endothelial cells and the various blood cells and plasma components they convey. We describe progress that has been made in applying computational modelling, quantitative biology and high-throughput experimentation to the angiogenesis process.

This computational study suggests that the rate of VEGF secretion by tumor cells may serve as a biomarker to predict the patient population that is likely to respond to anti-VEGF treatment. Thus, the model predictions have important clinical relevance and may aid clinicians and clinical researchers seeking interpretation of pharmacokinetic and pharmacodynamic observations and optimization of anti-VEGF therapies.

BackgroundVascular endothelial growth factor (VEGF) is a potent regulator of angiogenesis, and its role in cancer biology has been widely studied. Many cancer therapies target angiogenesis, with a focus being on VEGF-mediated signaling such as antibodies to VEGF. However, it is difficult to predict the effects of VEGF-neutralizing agents. We have developed a whole-body model of VEGF kinetics and transport under pathological conditions (in the presence of breast tumor). The model includes two major VEGF isoforms VEGF121 and VEGF165, receptors VEGFR1, VEGFR2 and co-receptors Neuropilin-1 and Neuropilin-2. We have added receptors on parenchymal cells (muscle fibers and tumor cells), and incorporated experimental data for the cell surface density of receptors on the endothelial cells, myocytes, and tumor cells. The model is applied to investigate the action of VEGF-neutralizing agents (called "anti-VEGF") in the treatment of cancer.ResultsThrough a sensitivity study, we examine how model parameters influence the level of free VEGF in the tumor, a measure of the response to VEGF-neutralizing drugs. We investigate the effects of systemic properties such as microvascular permeability and lymphatic flow, and of drug characteristics such as the clearance rate and binding affinity. We predict that increasing microvascular permeability in the tumor above 10-5 cm/s elicits the undesired effect of increasing tumor interstitial VEGF concentration beyond even the baseline level. We also examine the impact of the tumor microenvironment, including receptor expression and internalization, as well as VEGF secretion. We find that following anti-VEGF treatment, the concentration of free VEGF in the tumor can vary between 7 and 233 pM, with a dependence on both the density of VEGF receptors and co-receptors and the rate of neuropilin internalization on tumor cells. Finally, we predict that free VEGF in the tumor is reduced following anti-VEGF treatment when VEGF121 comprises at least 25% of the VEGF secreted by tumor cells.ConclusionsThis study explores the optimal drug characteristics required for an anti-VEGF agent to have a therapeutic effect and the tumor-specific properties that influence the response to therapy. Our model provides a framework for investigating the use of VEGF-neutralizing drugs for personalized medicine treatment strategies.

No abstract

Reprogramming of energy metabolism is a hallmark of cancer that enables the cancer cells to meet the increased energetic requirements due to uncontrolled proliferation. One prominent example is pancreatic ductal adenocarcinoma, an aggressive form of cancer with an overall 5-year survival rate of 5%. The reprogramming mechanism in pancreatic cancer involves deregulated uptake of glucose and glutamine and other opportunistic modes of satisfying energetic demands in a hypoxic and nutrient-poor environment. In the current study, we apply systems biology approaches to enable a better understanding of the dynamics of the distinct metabolic alterations in KRAS-mediated pancreatic cancer, with the goal of impeding early cell proliferation by identifying the optimal metabolic enzymes to target. We have constructed a kinetic model of metabolism represented as a set of ordinary differential equations that describe time evolution of the metabolite concentrations in glycolysis, glutaminolysis, tricarboxylic acid cycle and the pentose phosphate pathway. The model is comprised of 46 metabolites and 53 reactions. The mathematical model is fit to published enzyme knockdown experimental data. We then applied the model to perform in silico enzyme modulations and evaluate the effects on cell proliferation. Our work identifies potential combinations of enzyme knockdown, metabolite inhibition, and extracellular conditions that impede cell proliferation. Excitingly, the model predicts novel targets that can be tested experimentally. Therefore, the model is a tool to predict the effects of inhibiting specific metabolic reactions within pancreatic cancer cells, which is difficult to measure experimentally, as well as test further hypotheses toward targeted therapies.

BackgroundAngiogenesis is important in physiological and pathological conditions, as blood vessels provide nutrients and oxygen needed for tissue growth and survival. Therefore, targeting angiogenesis is a prominent strategy in both tissue engineering and cancer treatment. However, not all of the approaches to promote or inhibit angiogenesis lead to successful outcomes. Angiogenesis-based therapies primarily target pro-angiogenic factors such as vascular endothelial growth factor-A (VEGF) or fibroblast growth factor (FGF) in isolation. However, pre-clinical and clinical evidence shows these therapies often have limited effects. To improve therapeutic strategies, including targeting FGF and VEGF in combination, we need a quantitative understanding of the how the promoters combine to stimulate angiogenesis.ResultsIn this study, we trained and validated a detailed mathematical model to quantitatively characterize the crosstalk of FGF and VEGF intracellular signaling. This signaling is initiated by FGF binding to the FGF receptor 1 (FGFR1) and heparan sulfate glycosaminoglycans (HSGAGs) or VEGF binding to VEGF receptor 2 (VEGFR2) to promote downstream signaling. The model focuses on FGF- and VEGF-induced mitogen-activated protein kinase (MAPK) signaling and phosphorylation of extracellular regulated kinase (ERK), which promotes cell proliferation. We apply the model to predict the dynamics of phosphorylated ERK (pERK) in response to the stimulation by FGF and VEGF individually and in combination. The model predicts that FGF and VEGF have differential effects on pERK. Additionally, since VEGFR2 upregulation has been observed in pathological conditions, we apply the model to investigate the effects of VEGFR2 density and trafficking parameters. The model predictions show that these parameters significantly influence the response to VEGF stimulation.ConclusionsThe model agrees with experimental data and is a framework to synthesize and quantitatively explain experimental studies. Ultimately, the model provides mechanistic insight into FGF and VEGF interactions needed to identify potential targets for pro- or anti-angiogenic therapies.Electronic supplementary materialThe online version of this article (10.1186/s12918-018-0668-5) contains supplementary material, which is available to authorized users.