A main goal of mathematical and computational oncology is to develop quantitative tools to determine the most effective therapies for each individual patient. This involves predicting the right drug to be administered at the right time and at the right dose. Such an approach is known as precision medicine. Mathematical modelling can play an invaluable role in the development of such therapeutic strategies, since it allows for relatively fast, efficient and inexpensive simulations of a large number of treatment schedules in order to find the most effective. This review is a survey of mathematical models that explicitly take into account the spatial architecture of three-dimensional tumours and address tumour development, progression and response to treatments. In particular, we discuss models of epithelial acini, multicellular spheroids, normal and tumour spheroids and organoids, and multicomponent tissues. Our intent is to showcase how these in silico models can be applied to patient-specific data to assess which therapeutic strategies will be the most efficient. We also present the concept of virtual clinical trials that integrate standard-of-care patient data, medical imaging, organ-on-chip experiments and computational models to determine personalized medical treatment strategies.

A novel approach for the selection of step parameters as reaction coordinates in enhanced sampling simulations of DNA is presented. The method uses three atoms per base and does not require coordinate overlays or idealized base pairs. This allowed for a highly efficient implementation of the calculation of all step parameters and their Cartesian derivatives in molecular dynamics simulations. Good correlation between the calculated and actual twist, roll, tilt, shift, and slide parameters is obtained, while the correlation with rise is modest. The method is illustrated by its application to the methylated and unmethylated 5'-CATGTGACGTCACATG-3' double stranded DNA sequence. One-dimensional umbrella simulations indicate that the flexibility of the central CG step is only marginally affected by methylation.

A tumor vasculature that is functionally abnormal results in irregular gradients of metabolites and drugs within the tumor tissue. Recently, significant efforts have been committed to experimentally examine how cellular response to anti-cancer treatments varies based on the environment in which the cells are grown. In vitro studies point to specific conditions in which tumor cells can remain dormant and survive the treatment. In vivo results suggest that cells can escape the effects of drug therapy in tissue regions that are poorly penetrated by the drugs. Better understanding how the tumor microenvironments influence the emergence of drug resistance in both primary and metastatic tumors may improve drug development and the design of more effective therapeutic protocols. This chapter presents a hybrid agent-based model of the growth of tumor micrometastases and explores how microenvironmental factors can contribute to the development of acquired resistance in response to a DNA damaging drug. The specific microenvironments of interest in this work are tumor hypoxic niches and tumor normoxic sanctuaries with poor drug penetration. We aim to quantify how spatial constraints of limited drug transport and quiescent cell survival contribute to the development of drug resistant tumors.

The dynamics of tumor progression is driven by multiple factors, which can be exogenous to the tumor (microenvironment) or intrinsic (genetic, epigenetic or due to intercellular interactions). While tumor heterogeneity has been extensively studied on the level of cell genetic profiles or cellular composition, tumor morphological diversity has not been given as much attention. The limited analysis of tumor morphophenotypes may be attributed to the lack of accurate models, both experimental and computational, capable of capturing changes in tumor morphology with fine levels of spatial detail. Using a three-dimensional, agent-based, lattice-free computational model, we generated a library of multicellular tumor organoids, the experimental analogues of in vivo tumors. By varying three biologically relevant parameters-cell radius, cell division age and cell sensitivity to contact inhibition, we showed that tumor organoids with similar growth dynamics can express distinct morphologies and possess diverse cellular compositions. Taking advantage of the high-resolution of computational modeling, we applied the quantitative measures of compactness and accessible surface area, concepts that originated from the structural biology of proteins. Based on these analyses, we demonstrated that tumor organoids with similar sizes may differ in features associated with drug effectiveness, such as potential exposure to the drug or the extent of drug penetration. Both these characteristics might lead to major differences in tumor organoid's response to therapy. This indicates that therapeutic protocols should not be based solely on tumor size, but take into account additional tumor features, such as their morphology or cellular packing density.

Systemic chemotherapy is one of the main anticancer treatments used for most kinds of clinically diagnosed tumors. However, the efficacy of these drugs can be hampered by the physical attributes of the tumor tissue, such as tortuous vasculature, dense and fibrous extracellular matrix, irregular cellular architecture, tumor metabolic gradients, and non-uniform expression of the cell membrane receptors. This can impede the transport of therapeutic agents to tumor cells in sufficient quantities. In addition, tumor microenvironments undergo dynamic spatio-temporal changes during tumor progression and treatment, which can also obstruct drug efficacy. To examine ways to improve drug delivery on a cell-to-tissue scale (single-cell pharmacology), we developed the microscale pharmacokinetics/pharmacodynamics (microPKPD) modeling framework. Our model is modular and can be adjusted to include only the mathematical equations that are crucial for a biological problem under consideration. This modularity makes the model applicable to a broad range of pharmacological cases. As an illustration, we present two specific applications of the microPKPD methodology that help to identify optimal drug properties. The hypoxia-activated drugs example uses continuous drug concentrations, diffusive-advective transport through the tumor interstitium, and passive transmembrane drug uptake. The targeted therapy example represents drug molecules as discrete particles that move by diffusion and actively bind to cell receptors. The proposed modeling approach takes into account the explicit tumor tissue morphology, its metabolic landscape and/or specific receptor distribution. All these tumor attributes can be assessed from patients' diagnostic biopsies; thus, the proposed methodology can be developed into a tool suitable for personalized medicine, such as neoadjuvant chemotherapy.