Explore the vast range of titles available.

Myeloid-derived monocyte and macrophages are key cells in the bone that contribute to remodeling and injury repair. However, their temporal polarization status and control of bone-resorbing osteoclasts and bone-forming osteoblasts responses is largely unknown. In this study, we focused on two aspects of monocyte/macrophage dynamics and polarization states over time: 1) the injury-triggered pro- and anti-inflammatory monocytes/macrophages temporal profiles, 2) the contributions of pro- versus anti-inflammatory monocytes/macrophages in coordinating healing response. Bone healing is a complex multicellular dynamic process. While traditional in vitro and in vivo experimentation may capture the behavior of select populations with high resolution, they cannot simultaneously track the behavior of multiple populations. To address this, we have used an integrated coupled ordinary differential equations (ODEs)-based framework describing multiple cellular species to in vivo bone injury data in order to identify and test various hypotheses regarding bone cell populations dynamics. Our approach allowed us to infer several biological insights including, but not limited to,: 1) anti-inflammatory macrophages are key for early osteoclast inhibition and pro-inflammatory macrophage suppression, 2) pro-inflammatory macrophages are involved in osteoclast bone resorptive activity, whereas osteoblasts promote osteoclast differentiation, 3) Pro-inflammatory monocytes/macrophages rise during two expansion waves, which can be explained by the anti-inflammatory macrophages-mediated inhibition phase between the two waves. In addition, we further tested the robustness of the mathematical model by comparing simulation results to an independent experimental dataset. Taken together, this novel comprehensive mathematical framework allowed us to identify biological mechanisms that best recapitulate bone injury data and that explain the coupled cellular population dynamics involved in the process. Furthermore, our hypothesis testing methodology could be used in other contexts to decipher mechanisms in complex multicellular processes.

Spatial agent‐based models are frequently used to investigate the evolution of solid tumours subject to localized cell–cell interactions and microenvironmental heterogeneity. As spatial genomic, transcriptomic and proteomic technologies gain traction, spatial computational models are predicted to become ever more necessary for making sense of complex clinical and experimental data sets, for predicting clinical outcomes, and for optimizing treatment strategies. Here we present a non‐technical step by step guide to developing such a model from first principles. Stressing the importance of tailoring the model structure to that of the biological system, we describe methods of increasing complexity, from the basic Eden growth model up to off‐lattice simulations with diffusible factors. We examine choices that unavoidably arise in model design, such as implementation, parameterization, visualization and reproducibility. Each topic is illustrated with examples drawn from recent research studies and state of the art modelling platforms. We emphasize the benefits of simpler models that aim to match the complexity of the phenomena of interest, rather than that of the entire biological system. Our guide is aimed at both aspiring modellers and other biologists and oncologists who wish to understand the assumptions and limitations of the models on which major cancer studies now so often depend.

The Hybrid Automata Library (HAL) is a Java Library developed for use in mathematical oncology modeling. It is made of simple, efficient, generic components that can be used to model complex spatial systems. HAL's components can broadly be classified into: on- and off-lattice agent containers, finite difference diffusion fields, a GUI building system, and additional tools and utilities for computation and data collection. These components are designed to operate independently and are standardized to make them easy to interface with one another. As a demonstration of how modeling can be simplified using our approach, we have included a complete example of a hybrid model (a spatial model with interacting agent-based and PDE components). HAL is a useful asset for researchers who wish to build performant 1D, 2D and 3D hybrid models in Java, while not starting entirely from scratch. It is available on GitHub at https://github.com/MathOnco/HAL under the MIT License. HAL requires the Java JDK version 1.8 or later to compile and run the source code.

Heterogeneity in strategies for survival and proliferation among the cells that constitute a tumour is a driving force behind the evolution of resistance to cancer therapy. The rules mapping the tumour’s strategy distribution to the fitness of individual strategies can be represented as an evolutionary game. We develop a game assay to measure effective evolutionary games in co-cultures of non-small cell lung cancer cells that are sensitive and resistant to the anaplastic lymphoma kinase inhibitor alectinib. The games are not only quantitatively different between different environments, but targeted therapy and cancer-associated fibroblasts qualitatively switch the type of game being played by the in vitro population from Leader to Deadlock. This observation provides empirical confirmation of a central theoretical postulate of evolutionary game theory in oncology: we can treat not only the player, but also the game. Although we concentrate on measuring games played by cancer cells, the measurement methodology we develop can be used to advance the study of games in other microscopic systems by providing a quantitative description of non-cell-autonomous effects. Growth-rate data from co-cultures of drug-sensitive and drug-resistant cancer cells are quantitatively described using evolutionary game theory, revealing that drug treatment alters the type of game played by the cancer cells from ‘Leader’ to ‘Deadlock’.

Adaptive therapy is a dynamic cancer treatment protocol that updates (or \"adapts\") treatment decisions in anticipation of evolving tumor dynamics. This broad term encompasses many possible dynamic treatment protocols of patient-specific dose modulation or dose timing. Adaptive therapy maintains high levels of tumor burden to benefit from the competitive suppression of treatment-sensitive subpopulations on treatment-resistant subpopulations. This evolution-based approach to cancer treatment has been integrated into several ongoing or planned clinical trials, including treatment of metastatic castrate resistant prostate cancer, ovarian cancer, and BRAF-mutant melanoma. In the previous few decades, experimental and clinical investigation of adaptive therapy has progressed synergistically with mathematical and computational modeling. In this work, we discuss 11 open questions in cancer adaptive therapy mathematical modeling. The questions are split into three sections: 1) the necessary components of mathematical models of adaptive therapy 2) design and validation of dosing protocols, and 3) challenges and opportunities in clinical translation.

Bone-forming osteoblasts and -resorbing osteoclasts control bone injury repair, and myeloid-derived cells such as monocytes and macrophages are known to influence their behavior. However, precisely how these multiple cell types coordinate and regulate each other over time within the bone marrow to restore bone is difficult to dissect using biological approaches. Conversely, mathematical modeling lends itself well to this challenge. Therefore, we generated an ordinary differential equation (ODE) model powered by experimental data (osteoblast, osteoclast, bone volume, pro- and anti-inflammatory myeloid cells) obtained from intra-tibially injured mice. Initial ODE results using only osteoblast/osteoclast populations demonstrated that bone homeostasis could not be recovered after injury, but this issue was resolved upon integration of pro- and anti-inflammatory myeloid population dynamics. Surprisingly, the ODE revealed temporal disconnects between the peak of total bone mineralization/resorption, and osteoblast/osteoclast numbers. Specifically, the model indicated that osteoclast activity must vary greatly (> 17-fold) to return the bone volume to baseline after injury and suggest that osteoblast/osteoclast number alone is insufficient to predict bone the trajectory of bone repair. Importantly, the values of osteoclast activity fall within those published previously. These data underscore the value of mathematical modeling approaches to understand and reveal new insights into complex biological processes.

Multiple myeloma (MM) is an osteolytic malignancy that is incurable due to the emergence of treatment resistant disease. Defining how, when and where myeloma cell intrinsic and extrinsic bone microenvironmental mechanisms cause relapse is challenging with current biological approaches. Here, we report a biology-driven spatiotemporal hybrid agent-based model of the MM-bone microenvironment. Results indicate MM intrinsic mechanisms drive the evolution of treatment resistant disease but that the protective effects of bone microenvironment mediated drug resistance (EMDR) significantly enhances the probability and heterogeneity of resistant clones arising under treatment. Further, the model predicts that targeting of EMDR deepens therapy response by eliminating sensitive clones proximal to stroma and bone, a finding supported by in vivo studies. Altogether, our model allows for the study of MM clonal evolution over time in the bone microenvironment and will be beneficial for optimizing treatment efficacy so as to significantly delay disease relapse. Here, the authors develop a hybrid agent-based model to quantify the contributions of intrinsic cellular mechanisms and bone ecosystem factors to therapy resistance in multiple myeloma. They show that intrinsic mechanisms are essential for resistance, and that the bone microenvironment provides a protective niche that increases the likelihood.

Natural selection in cancer often results in the emergence of increasingly malignant tumor cells that display many if not all of the hallmarks of cancer. One of the most important traits acquired during cancer progression is angiogenesis. Tumor cells capable of secreting pro-angiogenic factors can be seen as cooperators where the improved oxygenation, nutrient delivery and waste disposal resulting from angiogenesis could be seen as a public good. Under this view, the relatively costly secretion of molecular signals required to orchestrate angiogenesis would be undertaken exclusively by cooperating tumor cells but the benefits of angiogenesis would be felt by neighboring tumor cells regardless of their contribution to the process. In this work we detail a mathematical model to better understand how clones capable of secreting pro-angiogenic factors can emerge in a tumor made of non-cooperative tumor cells. Given the importance of the spatial configuration of the tumor in determining the efficacy of the secretion of pro-angiogenic factors as well as the benefits of angiogenesis we have developed a spatial game theoretic approach where interactions and public good diffusion are described by graphs. The results show that structure of the population affects the evolutionary dynamics of the pro-angiogenic clone. Specifically, when the benefit of angiogenesis is represented by sigmoid function with regards to the number of pro-angiogenic clones then the probability of the coexistence of pro-angiogenic and angiogenesis-neutral clones increases. Our results demonstrate that pro-angiogenic clone equilibrates into clusters that appear from surrounding vascular tissues towards the center of tumor. These clusters appear notably less dense after anti-angiogenic therapy.

Intratumor heterogeneity is a feature of cancer that is associated with progression, treatment resistance, and recurrence. However, the mechanisms that allow diverse cancer cell lineages to coexist remain poorly understood. The storage effect is a coexistence mechanism that has been proposed to explain the diversity of a variety of ecological communities, including coral reef fish, plankton, and desert annual plants. Three ingredients are required for there to be a storage effect: (1) temporal variability in the environment, (2) buffered population growth, and (3) species-specific environmental responses. In this article, we argue that these conditions are observed in cancers and that it is likely that the storage effect contributes to intratumor diversity. Data that show the temporal variation within the tumor microenvironment are needed to quantify how cancer cells respond to fluctuations in the tumor microenvironment and what impact this has on interactions among cancer cell types. The presence of a storage effect within a patient’s tumors could have a substantial impact on how we understand and treat cancer.

Acute myeloid leukemia (AML) is a hematologic malignancy originating in the bone marrow and often progressing to extramedullary sites. Despite advances in molecularly targeted therapies and hematopoietic stem cell transplantation, clinical outcomes remain poor. Tyrosine kinase inhibitors (TKIs) provide benefit to a subset of AML patients harboring FLT3-ITD mutations; however, relapse and resistance remain common. These therapeutic failures are driven by both intrinsic properties of leukemic stem cells (LSCs)—a quiescent, self-renewing population—and extrinsic cues from the tumor microenvironment. We previously demonstrated that arteriolar endothelial cells (ECs) produce miR-126, which is transferred to LSCs, promoting quiescence, treatment resistance, and niche retention. During disease progression, TNF-α secreted by expanding blasts suppresses EC miR-126 production. Following TKI administration, blast reduction lowers TNF-ɑ levels, restoring EC miR-126 production, and this miR-126 expression enables LSCs to re-enter quiescence—thereby escaping therapy and facilitating relapse. To explore this dynamic, we developed an agent-based computational model of the AML bone marrow microenvironment, parameterized with in vitro and in vivo data. The model captures vascular niche remodeling and feedback between leukemic populations and endothelial signaling. Simulations reveal that LSC protection mediated by miR-126 can be disrupted by combining TKIs with miRisten, a miR-126 inhibitor. When administered on a defined schedule, this combination dismantles the protective niche and enhances LSC eradication. These findings underscore the therapeutic potential of targeting microenvironmental feedback to overcome resistance and prevent AML relapse.