Computational model of ductal carcinoma in situ: the effects of contact inhibition on pattern formation

A Multiscale Agent-Based Model of Ductal Carcinoma In Situ

OBJECTIVE

we present a multiscale agent-based model of Ductal Carcinoma in Situ (DCIS) in order to gain a detailed understanding of the cell-scale population dynamics, phenotypic distributions, and the associated interplay of important molecular signaling pathways that are involved in DCIS ductal invasion into the duct cavity (a process we refer to as duct advance rate here).

METHODS

DCIS is modeled mathematically through a hybridized discrete cell-scale model and a continuum molecular scale model, which are explicitly linked through a bidirectional feedback mechanism.

RESULTS

we find that duct advance rates occur in two distinct phases, characterized by an early exponential population expansion, followed by a long-term steady linear phase of population expansion, a result that is consistent with other modeling work. We further found that the rates were influenced most strongly by endocrine and paracrine signaling intensity, as well as by the effects of cell density induced quiescence within the DCIS population.

CONCLUSION

our model analysis identified a complex interplay between phenotypic diversity that may provide a tumor adaptation mechanism to overcome proliferation limiting conditions, allowing for dynamic shifts in phenotypic populations in response to variation in molecular signaling intensity. Further, sensitivity analysis determined DCIS axial advance rates and calcification rates were most sensitive to cell cycle time variation.

SIGNIFICANCE

this model may serve as a useful tool to study the cell-scale dynamics involved in DCIS initiation and intraductal invasion, and may provide insights into promising areas of future experimental research.

A hybrid agent-based model of the developing mammary terminal end bud

Mammary gland ductal elongation is spearheaded by terminal end buds (TEBs), where populations of highly proliferative cells are maintained throughout post-pubertal organogenesis in virgin mice until the mammary fat pad is filled by a mature ductal tree. We have developed a hybrid multiscale agent-based model to study how cellular differentiation pathways, cellular proliferation capacity, and endocrine and paracrine signaling play a role during development of the mammary gland. A simplified cellular phenotypic hierarchy that includes stem, progenitor, and fully differentiated cells within the TEB was implemented. Model analysis finds that mammary gland development was highly sensitive to proliferation events within the TEB, with progenitors likely undergoing 2-3 proliferation cycles before transitioning to a non-proliferative phenotype, and this result is in agreement with our previous experimental work. Endocrine and paracrine signaling were found to provide reliable ductal elongation rate regulation, while variations in the probability a new daughter cell will be of a proliferative phenotype were seen to have minimal effects on ductal elongation rates. Moreover, the distribution of cellular phenotypes within the TEB was highly heterogeneous, demonstrating significant allowable plasticity in possible phenotypic distributions while maintaining biologically relevant growth behavior. Finally, simulation results indicate ductal elongation rates due to cellular proliferation within the TEB may have a greater sensitivity to upstream endocrine signaling than endothelial to stromal paracrine signaling within the TEB. This model provides a useful tool to gain quantitative insights into cellular population dynamics and the effects of endocrine and paracrine signaling within the pubertal terminal end bud.

Mathematical modeling for novel cancer drug discovery and development

INTRODUCTION

Mathematical modeling enables: the in silico classification of cancers, the prediction of disease outcomes, optimization of therapy, identification of promising drug targets and prediction of resistance to anticancer drugs. In silico pre-screened drug targets can be validated by a small number of carefully selected experiments.

AREAS COVERED

This review discusses the basics of mathematical modeling in cancer drug discovery and development. The topics include in silico discovery of novel molecular drug targets, optimization of immunotherapies, personalized medicine and guiding preclinical and clinical trials. Breast cancer has been used to demonstrate the applications of mathematical modeling in cancer diagnostics, the identification of high-risk population, cancer screening strategies, prediction of tumor growth and guiding cancer treatment.

EXPERT OPINION

Mathematical models are the key components of the toolkit used in the fight against cancer. The combinatorial complexity of new drugs discovery is enormous, making systematic drug discovery, by experimentation, alone difficult if not impossible. The biggest challenges include seamless integration of growing data, information and knowledge, and making them available for a multiplicity of analyses. Mathematical models are essential for bringing cancer drug discovery into the era of Omics, Big Data and personalized medicine.

Automated reconstruction algorithm for identification of 3D architectures of cribriform ductal carcinoma in situ

Ductal carcinoma in situ (DCIS) is a pre-invasive carcinoma of the breast that exhibits several distinct morphologies but the link between morphology and patient outcome is not clear. We hypothesize that different mechanisms of growth may still result in similar 2D morphologies, which may look different in 3D. To elucidate the connection between growth and 3D morphology, we reconstruct the 3D architecture of cribriform DCIS from resected patient material. We produce a fully automated algorithm that aligns, segments, and reconstructs 3D architectures from microscopy images of 2D serial sections from human specimens. The alignment algorithm is based on normalized cross correlation, the segmentation algorithm uses histogram equilization, Otsu's thresholding, and morphology techniques to segment the duct and cribra. The reconstruction method combines these images in 3D. We show that two distinct 3D architectures are indeed found in samples whose 2D histological sections are similarly identified as cribriform DCIS. These differences in architecture support the hypothesis that luminal spaces may form due to different mechanisms, either isolated cell death or merging fronds, leading to the different architectures. We find that out of 15 samples, 6 were found to have 'bubble-like' cribra, 6 were found to have 'tube-like' criba and 3 were 'unknown.' We propose that the 3D architectures found, 'bubbles' and 'tubes', account for some of the heterogeneity of the disease and may be prognostic indicators of different patient outcomes.

Biopotential signals of breast cancer versus tumor types and proliferation stages

Clinical studies have shown compelling data of elevated biopotential signals recorded noninvasively from the breasts of women with breast cancer. While these data are compelling and show a strong potential for use in the noninvasive early detection of breast cancer, there remains significant knowledge gaps which must be addressed before this technology can be routinely used for breast cancer detection. A diffusion-drift model is developed to study the spatial and temporal characteristics of the biopotential signals of breast tumors taking into account the morphology and cell division stages. The electric signals of the most common tumor types-papillary, compact, and comedo-are also considered. The largest biopotential signal is observed from the compact tumor, while the smallest signal is observed from the papillary type. The results also show an increase in the time duration of the generated biopotential signals when cancer cells start their transitions at different time instants. The spatial and temporal variations of the biopotential signals are correlated with the tumor pattern which can have important implications for breast cancer detection.

Hybrid models of tumor growth

Cancer is a complex, multiscale process in which genetic mutations occurring at a subcellular level manifest themselves as functional changes at the cellular and tissue scale. The multiscale nature of cancer requires mathematical modeling approaches that can handle multiple intracellular and extracellular factors acting on different time and space scales. Hybrid models provide a way to integrate both discrete and continuous variables that are used to represent individual cells and concentration or density fields, respectively. Each discrete cell can also be equipped with submodels that drive cell behavior in response to microenvironmental cues. Moreover, the individual cells can interact with one another to form and act as an integrated tissue. Hybrid models form part of a larger class of individual-based models that can naturally connect with tumor cell biology and allow for the integration of multiple interacting variables both intrinsically and extrinsically and are therefore perfectly suited to a systems biology approach to tumor growth.

Current trends in mathematical modeling of tumor-microenvironment interactions: a survey of tools and applications

In its simplest description, a tumor is comprised of an expanding population of transformed cells supported by a surrounding microenvironment termed the tumor stroma. The tumor microenvironment has a very complex composition, including multiple types of stromal cells, a dense network of various extracellular matrix (ECM) fibers interpenetrated by the interstitial fluid and gradients of several chemical species that either are dissolved in the fluid or are bound to the ECM structure. In order to study experimentally such complex interactions between multiple players, cancer is dissected and considered at different scales of complexity, such as protein interactions, biochemical pathways, cellular functions or whole organism studies. However, the integration of information acquired from these studies into a common description is as difficult as the disease itself. Computational models of cancer can provide cancer researchers with invaluable tools that are capable of integrating the complexity into organizing principles as well as suggesting testable hypotheses. We will focus in this Minireview on mathematical models in which the whole cell is a main modeling unit. We will present a current stage of such cell-focused mathematical modeling incorporating different stromal components and their interactions with growing tumors, and discuss what modeling approaches can be undertaken to complement the in vivo and in vitro experimentation.