Oxygen consumption dynamics in steady-state tumour models

Pig manure degradation and carbon emission: Measuring and modeling combined aerobic-anaerobic transformations

Greenhouse gas emissions from liquid livestock manure storage significantly contribute to global warming. Accurate farm-scale models are essential for predicting these emissions and evaluating manure management strategies, but they rely on multiple parameters describing carbon loss dynamics. Surface respiration may significantly influence carbon loss and methane emission, yet it is not explicitly included in current models. We conducted experiments to measure pig manure surface respiration rate and its effect on organic matter degradation and methane and carbon dioxide emissions. Manure was incubated for 283 days at 10°C or 20°C under aerobic or anaerobic conditions, while measuring methane and carbon dioxide emission. This was followed by anaerobic digestion at 38°C. Surface respiration reduced the organic matter content, and the effect was temperature dependent. Methane emission was not affected by surface respiration, suggesting that substrate availability was not rate-limiting for methanogenesis. Surface respiration rates were 18.1 ± 3.5 g CO2 m-2 day-1 at 10°C and 37.1 ± 13.1 g CO2 m-2 day-1 at 20°C (mean ± standard deviation) and were consistent with microsensor measurements of oxygen consumption in different manure surfaces. Based on these results, temperature- and surface area-dependent respiration was incorporated in the existing anaerobic biodegradation model (ABM). Simulations showed that surface respiration accounts for 29% of carbon losses in a typical pig house and 8% for outdoor storage. Developing and refining algorithms for diverse carbon transformations, such as surface respiration, is crucial for evaluating the potential for methane emission and identification of variables that control emissions at the farm scale.

The 'bIUreactor': An Open-Source 3D Tissue Research Platform

We developed the open-source bIUreactor research platform for studying 3D structured tissues. The versatile and modular platform allows a researcher to generate 3D tissues, culture them with oxygenated perfusion, and provide cyclic loading, all in their own lab (in laboratorium) for an all in cost of $8,000 including 3D printer, printing resin, and electronics. We achieved this by applying a design philosophy that leverages 3D printing, open-source software and hardware, and practical techniques to produce the following: 1. perfusible 3D tissues, 2. a bioreactor chamber for tissue culture, 3. a module for applying cyclic compression, 4. a peristaltic pump for providing oxygenated perfusion to 3D tissues, 5. motor control units, and 6. open-source code for running the control units. By making it widely available for researchers to investigate 3D tissue models and easy for them to use, we intend for the bIUreactor to democratize 3D tissue research, therefore increasing the pace and scale of biomedical research discoveries using 3D tissue models.

A seven-step guide to spatial, agent-based modelling of tumour evolution

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.

Knockout cancer by nano-delivered immunotherapy using perfusion-aided scaffold-based tumor-on-a-chip

Cancer is a multifactorial disease produced by mutations in the oncogenes and tumor suppressor genes, which result in uncontrolled cell proliferation and resistance to cell death. Cancer progresses due to the escape of altered cells from immune monitoring, which is facilitated by the tumor's mutual interaction with its microenvironment. Understanding the mechanisms involved in immune surveillance evasion and the significance of the tumor microenvironment might thus aid in developing improved therapies. Although in vivo models are commonly utilized, they could be better for time, cost, and ethical concerns. As a result, it is critical to replicate an in vivo model and recreate the cellular and tissue-level functionalities. A 3D cell culture, which gives a 3D architecture similar to that found in vivo, is an appropriate model. Furthermore, numerous cell types can be cocultured, establishing cellular interactions between TME and tumor cells. Moreover, microfluidics perfusion can provide precision flow rates, thus simulating tissue/organ function. Immunotherapy can be used with the perfused 3D cell culture technique to help develop successful therapeutics. Immunotherapy employing nano delivery can target the spot and silence the responsible genes, ensuring treatment effectiveness while minimizing adverse effects. This study focuses on the importance of 3D cell culture in understanding the pathophysiology of 3D tumors and TME, the function of TME in drug resistance, tumor progression, and the development of advanced anticancer therapies for high-throughput drug screening.

AngioMT: An in silico platform for digital sensing of oxygen transport through heterogenous microvascular networks

Measuring the capacity of microvascular networks in delivering soluble oxygen and nutrients to its organs is essential in health, disease, and surgical interventions. Here, a finite element method-based oxygen transport program, AngioMT, is designed and validated to predict spatial oxygen distribution and other physiologically relevant transport metrics within both the vascular network and the surrounding tissue. The software processes acquired images of microvascular networks and produces a digital mesh which is used to predict vessel and tissue oxygenation. The image-to-physics translation by AngioMT correlated with results from commercial software, however only AngioMT could provide predictions within the solid tissue in addition to vessel oxygenation. AngioMT predictions were sensitive and positively correlated to spatial heterogeneity and extent of vascularization of 500 different vascular networks formed with variable vasculogenic conditions. The predictions of AngioMT cross-correlate with experimentally-measured oxygen distributions in vivo. The computational power of the software is increased by including calculations of higher order reaction mechanisms, and the program includes defining additional organ and tissue structures for a more physiologically relevant analysis of tissue oxygenation in complex co-cultured systems, or in vivo. AngioMT may serve as a digital performance measuring tool of vascular networks in microcirculation, experimental models of vascularized tissues and organs, and in clinical applications, such as organ transplants.

Hypoxia-Responsive Luminescent CEST MRI Agent for In Vitro and In Vivo Tumor Detection and Imaging

Hypoxia is a feature of most solid tumors and a key determinant of cancer growth and propagation. Sensing hypoxia effectively could lead to more favorable clinical outcomes. Here, we report a molecular antenna-based bimodal probe designed to exploit the complementary advantages of magnetic resonance (MR)- and optical-based imaging. Specifically, we describe the synthesis and evaluation of a dual-action probe (NO₂-Eu) that permits hypoxia-activated chemical exchange saturation transfer (CEST) MR and optical imaging. In CT26 cells, this NO₂-Eu probe not only provides an enhanced CEST MRI signal but also turns "on" the optical signal under hypoxic conditions. Time-dependent in vivo CEST imaging in a hypoxic CT26 tumor xenograft mouse model revealed probe-dependent tumor detection by CEST MRI contrast in the tumor area. We thus suggest that dual-action hypoxia probes, like that reported here, could have a role to play in solid tumor diagnosis and monitoring.

Endocrine disruptor chlorpyrifos promotes migration, invasion, and stemness phenotype in 3D cultures of breast cancer cells and induces a wide range of pathways involved in cancer progression

Organophosphorus chlorpyrifos (CPF) is currently considered an endocrine disruptor (ED), as it can imitate hormone actions both in vitro and in vivo. We recently reported that CPF induces migration and invasion in 2D cultures and changes the expression of key molecular markers involved in epithelial mesenchymal transition in MCF-7 and MDA-MB-231 cell lines. In this study, we investigated whether CPF could behave as a predisposing factor for tumors to become more metastatic and aggressive using 3D culture models. In MCF-7 cells, 0.05 μM CPF induced an increase in the number and size of mammospheres via estrogen receptor alpha (ERα) and c-SRC. Furthermore, 0.05 μM CPF increased the area of spheroids generated from MCF-7 cells, induced invasion using both Matrigel® and type 1 collagen matrices, and increased cell migration capacity via ERα in this 3D model. In turn, 50 μM CPF increased cell migration capacity and invasion using type 1 collagen matrix. In monolayers, CPF increased the phosphorylation and membrane translocation of c-SRC at both concentrations assayed. CPF at 0.05 μM boosted p-AKT, p-GSK-3β and p-P38. While p-AKT rose in a ERα-dependent way, p-GSK-3β was dependent on ERα- and c-SRC, and p-P38 was only dependent on c-SRC. On the other hand, the increase in p-AKT and p-P38 induced by 50 μM CPF was dependent on the c-SRC pathway. We also observed that 0.05 μM CPF increased IGF-1R and IRS-1 expression and that 50 μM CPF induced IGF-1Rβ phosphorylation. In the MDA-MB-231 cell line, 0.05 and 50 μM CPF increased p-c-SRC. Finally, p-AKT and p-GSK-3β were also induced by CPF at 0.05 and 50 μM, and an increase in p-P38 was observed at 50 μM. Taken together, these data provide support for the notion that CPF may represent a risk factor for breast cancer development and progression.

A Computational Model for Oxygen Depletion Hypothesis in FLASH Effect

Experiments have reported low normal tissue toxicities during FLASH irradiation, but the mechanism has not been elaborated. Several hypotheses have been proposed to explain the mechanism. One hypothesis is oxygen depletion. We analyze the time-dependent change of oxygen concentration in the tissue to study the oxygen depletion hypothesis using a computational model. The effects of physical, chemical and physiological parameters on oxygen depletion were explored. The kinetic equation of the model is solved numerically using the finite difference method with rational boundary conditions. Results of oxygen distribution is supported by the experiments of oxygen-sensitivity electrodes and experiments on the expression and distribution of the hypoxia-inducible factors. The analysis of parameters shows that the steady-state oxygen distribution before irradiation is determined by the oxygen consumption rate of the tissue and the microvessel density. The change of oxygen concentration after irradiation has been found to follow a negative exponential function, and the time constant is mainly determined by the microvessel density. The change of oxygen during exposure increases with dose rate and tends to be saturated because of oxygen diffusion. When the dose rate is high enough, the same dose results in the same reduction of oxygen concentration regardless of dose rate. The analysis of the FLASH effect in the brain tissue based on this model does not support the explanation of the oxygen depletion hypothesis. The oxygen depletion hypothesis remains controversial because the oxygen in most normal tissues cannot be depleted to radiation resistance level by FLASH irradiation.

Advances in Engineered Three-Dimensional (3D) Body Articulation Unit Models

Indeed, the body articulation units, commonly referred to as body joints, play significant roles in the musculoskeletal system, enabling body flexibility. Nevertheless, these articulation units suffer from several pathological conditions, such as osteoarthritis (OA), rheumatoid arthritis (RA), ankylosing spondylitis, gout, and psoriatic arthritis. There exist several treatment modalities based on the utilization of anti-inflammatory and analgesic drugs, which can reduce or control the pathophysiological symptoms. Despite the success, these treatment modalities suffer from major shortcomings of enormous cost and poor recovery, limiting their applicability and requiring promising strategies. To address these limitations, several engineering strategies have been emerged as promising solutions in fabricating the body articulation as unit models towards local articulation repair for tissue regeneration and high-throughput screening for drug development. In this article, we present challenges related to the selection of biomaterials (natural and synthetic sources), construction of 3D articulation models (scaffold-free, scaffold-based, and organ-on-a-chip), architectural designs (microfluidics, bioprinting, electrospinning, and biomineralization), and the type of culture conditions (growth factors and active peptides). Then, we emphasize the applicability of these articulation units for emerging biomedical applications of drug screening and tissue repair/regeneration. In conclusion, we put forward the challenges and difficulties for the further clinical application of the in vitro 3D articulation unit models in terms of the long-term high activity of the models.

3D Cell Culture Models as Recapitulators of the Tumor Microenvironment for the Screening of Anti-Cancer Drugs

Today, innovative three-dimensional (3D) cell culture models have been proposed as viable and biomimetic alternatives for initial drug screening, allowing the improvement of the efficiency of drug development. These models are gaining popularity, given their ability to reproduce key aspects of the tumor microenvironment, concerning the 3D tumor architecture as well as the interactions of tumor cells with the extracellular matrix and surrounding non-tumor cells. The development of accurate 3D models may become beneficial to decrease the use of laboratory animals in scientific research, in accordance with the European Union's regulation on the 3R rule (Replacement, Reduction, Refinement). This review focuses on the impact of 3D cell culture models on cancer research, discussing their advantages, limitations, and compatibility with high-throughput screenings and automated systems. An insight is also given on the adequacy of the available readouts for the interpretation of the data obtained from the 3D cell culture models. Importantly, we also emphasize the need for the incorporation of additional and complementary microenvironment elements on the design of 3D cell culture models, towards improved predictive value of drug efficacy.

Modeling of cellular response after FLASH irradiation: a quantitative analysis based on the radiolytic oxygen depletion hypothesis

Purpose.Recent studies suggest ultra-high dose rate (FLASH) irradiation can spare normal tissues from radiotoxicity, while efficiently controlling the tumor, and this is known as the 'FLASH effect'. This study performed theoretical analyses about the impact of radiolytic oxygen depletion (ROD) on the cellular responses after FLASH irradiation.Methods.Monte Carlo simulation was used to model the ROD process, determine the DNA damage, and calculate the amount of oxygen depleted (L_ROD) during FLASH exposure. A mathematical model was applied to analyze oxygen tension (pO₂) distribution in human tissues and the recovery of pO₂after FLASH irradiation. DNA damage and cell survival fractions (SFs) after FLASH irradiation were calculated. The impact of initial cellular pO₂, FLASH pulse number, pulse interval, and radiation quality of the source particles on ROD and subsequent cellular responses were systematically evaluated.Results.The simulated electronL_RODrange was 0.38-0.43μM Gy^-1when pO₂ranged from 7.5 to 160 mmHg. The calculated DNA damage and SFs show that the radioprotective effect is only evident in cells with a low pO₂. Different irradiation setups alter the cellular responses by modifying the pO₂. Single pulse delivery or multi-pulse delivery with pulse intervals shorter than 10-50 ms resulted in fewer DNA damages and higher SFs. Source particles with a low linear energy transfer (LET) have a higher capacity to deplete oxygen, and thus, lead to a more conspicuous radioprotective effect.Conclusions. A systematic analysis of the cellular response following FLASH irradiation was performed to provided suggestions for future FLASH applications. The FLASH radioprotective effect due to ROD may only be observed in cells with a low pO₂. Single pulse delivery or multi-pulse delivery with short pulse intervals are suggested for FLASH irradiation to avoid oxygen tension recovery during pulse intervals. Source particles with low LET are preferred for their conspicuous radioprotective effects.

Microdissected Tissue vs Tissue Slices-A Comparative Study of Tumor Explant Models Cultured On-Chip and Off-Chip

Predicting patient responses to anticancer drugs is a major challenge both at the drug development stage and during cancer treatment. Tumor explant culture platforms (TECPs) preserve the native tissue architecture and are well-suited for drug response assays. However, tissue longevity in these models is relatively low. Several methodologies have been developed to address this issue, although no study has compared their efficacy in a controlled fashion. We investigated the effect of two variables in TECPs, specifically, the tissue size and culture vessel on tissue survival using micro-dissected tumor tissue (MDT) and tissue slices which were cultured in microfluidic chips and plastic well plates. Tumor models were produced from ovarian and prostate cancer cell line xenografts and were matched in terms of the specimen, total volume of tissue, and respective volume of medium in each culture system. We examined morphology, viability, and hypoxia in the various tumor models. Our observations suggest that the viability and proliferative capacity of MDTs were not affected during the time course of the experiments. In contrast, tissue slices had reduced proliferation and showed increased cell death and hypoxia under both culture conditions. Tissue slices cultured in microfluidic devices had a lower degree of hypoxia compared to those in 96-well plates. Globally, our results show that tissue slices have lower survival rates compared to MDTs due to inherent diffusion limitations, and that microfluidic devices may decrease hypoxia in tumor models.

A therapeutic vascular conduit to support in vivo cell-secreted therapy

A significant barrier to implementation of cell-based therapies is providing adequate vascularization to provide oxygen and nutrients. Here we describe an approach for cell transplantation termed the Therapeutic Vascular Conduit (TVC), which uses an acellular vessel as a scaffold for a hydrogel sheath containing cells designed to secrete a therapeutic protein. The TVC can be directly anastomosed as a vascular graft. Modeling supports the concept that the TVC allows oxygenated blood to flow in close proximity to the transplanted cells to prevent hypoxia. As a proof-of-principle study, we used erythropoietin (EPO) as a model therapeutic protein. If implanted as an arteriovenous vascular graft, such a construct could serve a dual role as an EPO delivery platform and hemodialysis access for patients with end-stage renal disease. When implanted into nude rats, TVCs containing EPO-secreting fibroblasts were able to increase serum EPO and hemoglobin levels for up to 4 weeks. However, constitutive EPO expression resulted in macrophage infiltration and luminal obstruction of the TVC, thus limiting longer-term efficacy. Follow-up in vitro studies support the hypothesis that EPO also functions to recruit macrophages. The TVC is a promising approach to cell-based therapeutic delivery that has the potential to overcome the oxygenation barrier to large-scale cellular implantation and could thus be used for a myriad of clinical disorders. However, a complete understanding of the biological effects of the selected therapeutic is absolutely essential.

Fabrication of PNIPAm-based thermoresponsive hydrogel microwell arrays for tumor spheroid formation

Complex three-dimensional (3D) cell cultures are being increasingly implemented in biomedical research as they provide important insights into complex cancer biology, and cell-cell and cell-matrix interactions in the tumor microenvironment. However, most methods used today for 3D cell culture are limited by high cost, the need for specialized skills, low throughput and the use of unnatural culture environments. We report the development of a unique biomimetic hydrogel microwell array platform for the generation and stress-free isolation of cancer spheroids. The poly N-isopropylacrylamide-based hydrogel microwell array (PHMA) has thermoresponsive properties allowing for the attachment and growth of cell aggregates/ spheroids at 37 °C, and their easy isolation at room temperature (RT). The reversible phase transition of the microwell arrays at 35 °C was confirmed visually and by differential scanning calorimetry. Swelling/ shrinking studies and EVOS imaging established that the microwell arrays are hydrophilic and swollen at temperatures <35 °C, while they shrink and are hydrophobic at temperatures >35 °C. Spheroid development within the PHMA was optimized for seeding density, incubation time and cell viability. Spheroids of A549, HeLa and MG-63 cancer cell lines, and human lung fibroblast (HLF) cell line generated within the PHMAs had relatively spherical morphology with hypoxic cores. Finally, using MG-63 cell spheroids as representative models, a proof-of-concept drug response study using doxorubicin hydrochloride was conducted. Overall, we demonstrate that the PHMAs are an innovative alternative to currently used 3D cell culture techniques, for the high-throughput generation of cell spheroids for disease modeling and drug screening applications.

3D spheroid models of paediatric SHH medulloblastoma mimic tumour biology, drug response and metastatic dissemination

Studying medulloblastoma, the most common malignant paediatric brain tumour, requires simple yet realistic in vitro models. In this study, we optimised a robust, reliable, three-dimensional (3D) culture method for medulloblastoma able to recapitulate the spatial conformation, cell-cell and cell-matrix interactions that exist in vivo and in patient tumours. We show that, when grown under the same stem cell enriching conditions, SHH subgroup medulloblastoma cell lines established tight, highly reproducible 3D spheroids that could be maintained for weeks in culture and formed pathophysiological oxygen gradients. 3D spheroid culture also increased resistance to standard-of-care chemotherapeutic drugs compared to 2D monolayer culture. We exemplify how this model can enhance in vitro therapeutic screening approaches through dual-inhibitor studies and continual monitoring of drug response. Next, we investigated the initial stages of metastatic dissemination using brain-specific hyaluronan hydrogel matrices. RNA sequencing revealed downregulation of cell cycle genes and upregulation of cell movement genes and key fibronectin interactions in migrating cells. Analyses of these upregulated genes in patients showed that their expression correlated with early relapse and overall poor prognosis. Our 3D spheroid model is a significant improvement over current in vitro techniques, providing the medulloblastoma research community with a well-characterised and functionally relevant culture method.

Determining the parameter space for effective oxygen depletion for FLASH radiation therapy

There has been a recent revival of interest in the FLASH effect, after experiments have shown normal tissue sparing capabilities of ultra-high-dose-rate radiation with no compromise on tumour growth restraint. A model has been developed to investigate the relative importance of a number of fundamental parameters considered to be involved in the oxygen depletion paradigm of induced radioresistance. An example eight-dimensional parameter space demonstrates the conditions under which radiation may induce sufficient depletion of oxygen for a diffusion-limited hypoxic cellular response. Initial results support experimental evidence that FLASH sparing is only achieved for dose rates on the order of tens of Gy/s or higher, for a sufficiently high dose, and only for tissue that is slightly hypoxic at the time of radiation. We show that the FLASH effect is the result of a number of biological, radiochemical and delivery parameters. Also, the threshold dose for a FLASH effect occurring would be more prominent when the parameterisation was optimised to produce the maximum effect. The model provides a framework for further FLASH-related investigation and experimental design. An understanding of the mechanistic interactions producing an optimised FLASH effect is essential for its translation into clinical practice.

Advanced Multi-Dimensional Cellular Models as Emerging Reality to Reproduce In Vitro the Human Body Complexity

A hot topic in biomedical science is the implementation of more predictive in vitro models of human tissues to significantly improve the knowledge of physiological or pathological process, drugs discovery and screening. Bidimensional (2D) culture systems still represent good high-throughput options for basic research. Unfortunately, these systems are not able to recapitulate the in vivo three-dimensional (3D) environment of native tissues, resulting in a poor in vitro–in vivo translation. In addition, intra-species differences limited the use of animal data for predicting human responses, increasing in vivo preclinical failures and ethical concerns. Dealing with these challenges, in vitro 3D technological approaches were recently bioengineered as promising platforms able to closely capture the complexity of in vivo normal/pathological tissues. Potentially, such systems could resemble tissue-specific extracellular matrix (ECM), cell–cell and cell–ECM interactions and specific cell biological responses to mechanical and physical/chemical properties of the matrix. In this context, this review presents the state of the art of the most advanced progresses of the last years. A special attention to the emerging technologies for the development of human 3D disease-relevant and physiological models, varying from cell self-assembly (i.e., multicellular spheroids and organoids) to the use of biomaterials and microfluidic devices has been given.

Modeling the effect of oxygen on the chemical stage of water radiolysis using GPU-based microscopic Monte Carlo simulations, with an application in FLASH radiotherapy

Oxygen plays a critical role in determining the initial DNA damages induced by ionizing radiation. It is important to mechanistically model the oxygen effect in the water radiolysis process. However, due to the computational costs from the many body interaction problem, oxygen is often ignored or treated as a constant continuum radiolysis-scavenger background in the simulations using common microscopic Monte Carlo tools. In this work, we reported our recent progress on the modeling of the chemical stage of the water radiolysis with an explicit consideration of the oxygen effect, based upon our initial development of an open-source graphical processing unit (GPU)-based MC simulation tool, gMicroMC. The inclusion of oxygen mainly reduces the yields of [Formula: see text] and [Formula: see text] chemical radicals, turning them into highly toxic [Formula: see text] and [Formula: see text] species. To demonstrate the practical value of gMicroMC in large scale simulation problems, we applied the oxygen-simulation-enabled gMicroMC to compute the yields of chemical radicals under a high instantaneous dose rate [Formula: see text] to study the oxygen depletion hypothesis in FLASH radiotherapy. A decreased oxygen consumption rate (OCR) was found associated with a reduced initial oxygen concentration level due to reduced probabilities of reactions. With respect to dose rate, for the oxygen concentration of 21% and electron energy of 4.5 [Formula: see text], OCR remained approximately constant (∼0.22 [Formula: see text]) for [Formula: see text]'s of [Formula: see text], [Formula: see text] and reduced to 0.19 [Formula: see text] at [Formula: see text], because the increased dose rate improved the mutual reaction frequencies among radicals, hence reducing their reactions with oxygen. We computed the time evolution of oxygen concentration under the FLASH irradiation setups. At the dose rate of [Formula: see text] and initial oxygen concentrations from 0.01% to 21%, the oxygen is unlikely to be fully depleted with an accumulative dose of 30 Gy, which is a typical dose used in FLASH experiments. The computational efficiency of gMicroMC when considering oxygen molecules in the chemical stage was evaluated through benchmark work to GEANT4-DNA with simulating an equivalent number of radicals. With an initial oxygen concentration of 3% (∼10⁵ molecules), a speedup factor of 1228 was achieved for gMicroMC on a single GPU card when comparing with GEANT4-DNA on a single CPU.

A 3D View of Colorectal Cancer Models in Predicting Therapeutic Responses and Resistance

Although there have been many advances in recent years for the treatment of colorectal cancer (CRC), it still remains the third most common cause of cancer-related deaths worldwide. Many patients with late stage CRC display resistance to multiple different therapeutics. An important aspect in developing effective therapeutics for CRC patients is understanding the interactions that take place in the tumor microenvironment (TME), as it has been shown to contribute to drug resistance in vivo. Much research over the past 100 years has focused on 2D monolayer cultures or in vivo studies, however, the efficacy in translating these to the clinic is very low. More recent studies are turning towards developing an effective 3D model of CRC that is clinically relevant, that can recapitulate the TME in vitro and bridge the gap between 2D cultures and in vivo studies, with the aim of reducing the use of animal models in the future. This review summarises the advantages and limitations of different 3D CRC models. It emphasizes how different 3D models may be optimised to study cellular and extracellular interactions that take place in the TME of CRC in an effort to allow the development of more translatable effective treatment options for patients.

Bioengineering the Oxygen-Deprived Tumor Microenvironment Within a Three-Dimensional Platform for Studying Tumor-Immune Interactions

Oxygen deprivation within tumors is one of the most prevalent causes of resilient cancer cell survival and increased immune evasion in breast cancer (BCa). Current in vitro models do not adequately mimic physiological oxygen levels relevant to breast tissue and its tumor-immune interactions. In this study, we propose an approach to engineer a three-dimensional (3D) model (named 3D engineered oxygen, 3D-O) that supports the growth of BCa cells and generates physio- and pathophysiological oxygen levels to understand the role of oxygen availability in tumor-immune interactions. BCa cells (MDA-MB-231 and MCF-7) were embedded into plasma-derived 3D-O scaffolds that reflected physio- and pathophysiological oxygen levels relevant to the healthy and cancerous breast tissue. BCa cells grown within 3D-O scaffolds were analyzed by flow cytometry, confocal imaging, immunohistochemistry/immunofluorescence for cell proliferation, extracellular matrix protein expression, and alterations in immune evasive outcomes. Exosome secretion from 3D-O scaffolds were evaluated using the NanoSight particle analyzer. Peripheral blood mononuclear cells were incorporated on the top of 3D-O scaffolds and the difference in tumor-infiltrating capabilities as a result of different oxygen content were assessed by flow cytometry and confocal imaging. Lastly, hypoxia and Programmed death-ligand 1 (PD-L1) inhibition were validated as targets to sensitize BCa cells in order to overcome immune evasion. Low oxygen-induced adaptations within 3D-O scaffolds validated known tumor hypoxia characteristics such as reduced BCa cell proliferation, increased extracellular matrix protein expression, increased extracellular vesicle secretion and enhanced immune surface marker expression on BCa cells. We further demonstrated that low oxygen in 3D-O scaffolds significantly influence immune infiltration. CD8+ T cell infiltration was impaired under pathophysiological oxygen levels and we were also able to establish that hypoxia and PD-L1 inhibition re-sensitized BCa cells to cytotoxic CD8+ T cells. Bioengineering the oxygen-deprived BCa tumor microenvironment in our engineered 3D-O physiological and tumorous scaffolds supported known intra-tumoral hypoxia characteristics allowing the study of the role of oxygen availability in tumor-immune interactions. The 3D-O model could serve as a promising platform for the evaluation of immunological events and as a drug-screening platform tool to overcome hypoxia-driven immune evasion.

DEP-Dots for 3D cell culture: low-cost, high-repeatability, effective 3D cell culture in multiple gel systems

It is known that cells grown in 3D are more tolerant to drug treatment than those grown in dispersion, but the mechanism for this is still not clear; cells grown in 3D have opportunities to develop inter-cell communication, but are also closely packed which may impede diffusion. In this study we examine methods for dielectrophoresis-based cell aggregation of both suspension and adherent cell lines, and compare the effect of various drugs on cells grown in 3D and 2D. Comparing viability of pharmacological interventions on 3D cell clusters against both suspension cells and adherent cells grown in monolayer, as well as against a unicellular organism with no propensity for intracellular communication, we suggest that 3D aggregates of adherent cells, compared to suspension cells, show a substantially different drug response to cells grown in monolayer, which increases as the IC₅₀ is approached. Further, a mathematical model of the system for each agent demonstrates that changes to drug response are due to inherent changes in the system of adherent cells from the 2D to 3D state. Finally, differences in the electrophysiological membrane properties of the adherent cell type suggest this parameter plays an important role in the differences found in the 3D drug response.

Trends in Bone Metastasis Modeling

Bone is one of the most common sites for cancer metastasis. Bone tissue is composed by different kinds of cells that coexist in a coordinated balance. Due to the complexity of bone, it is impossible to capture the intricate interactions between cells under either physiological or pathological conditions. Hence, a variety of in vivo and in vitro approaches have been developed. Various models of tumor-bone diseases are routinely used to provide valuable information on the relationship between metastatic cancer cells and the bone tissue. Ideally, when modeling the metastasis of human cancers to bone, models would replicate the intra-tumor heterogeneity, as well as the genetic and phenotypic changes that occur with human cancers; such models would be scalable and reproducible to allow high-throughput investigation. Despite the continuous progress, there is still a lack of solid, amenable, and affordable models that are able to fully recapitulate the biological processes happening in vivo, permitting a correct interpretation of results. In the last decades, researchers have demonstrated that three-dimensional (3D) methods could be an innovative approach that lies between bi-dimensional (2D) models and animal models. Scientific evidence supports that the tumor microenvironment can be better reproduced in a 3D system than a 2D cell culture, and the 3D systems can be scaled up for drug screening in the same way as the 2D systems thanks to the current technologies developed. However, 3D models cannot completely recapitulate the inter- and intra-tumor heterogeneity found in patients. In contrast, ex vivo cultures of fragments of bone preserve key cell-cell and cell-matrix interactions and allow the study of bone cells in their natural 3D environment. Moreover, ex vivo bone organ cultures could be a better model to resemble the human pathogenic metastasis condition and useful tools to predict in vivo response to therapies. The aim of our review is to provide an overview of the current trends in bone metastasis modeling. By showing the existing in vitro and ex vivo systems, we aspire to contribute to broaden the knowledge on bone metastasis models and make these tools more appealing for further translational studies.

Oxygenation Profiles of Human Blood, Cell Culture Medium, and Water for Perfusion of 3D-Bioprinted Tissues using the FABRICA Bioreactor Platform

Persistent and saturated oxygen distribution from perfusion media (i.e., blood, or cell culture media) to cells within cell-dense, metabolically-active biofabricated tissues is required to keep them viable. Improper or poor oxygen supply to cells within the tissue bulk severely limits the tissue culturing potential of many bioreactors. We added an oxygenator module to our modular FABRICA bioreactor in order to provide stable oxygenation to biofabricated tissues during culture. In this proof of concept study of an oxygenated and perfused bioreactor, we characterized the oxygenation of water, cell culture medium, and human blood in the FABRICA as functions of augmenting vacuum (air inlet) pressure, perfusion (volumetric flow) rate, and tubing/oxygenator components. The mean oxygen levels for water and cell culture media were 27.7 ± 2.1% and 27.6 ± 4.1%, respectively. The mean oxygen level for human blood was 197.0 ± 90.0 mmHg, with near-physiologic levels achieved with low-permeability PharMed tubing alone (128.0 ± 14.0 mmHg). Hematologic values pre- and post-oxygenation, respectively were (median ± IQR): Red blood cell: 6.0 ± 0.5 (10⁶/μL) and 6.5 ± 0.4 (10⁶/μL); Hemoglobin: 17.5 ± 1.2 g/dL and 19.2 ± 3.0 g/dL; and Hematocrit: 56.7 ± 2.4% and 61.4 ± 7.5%. The relative stability of the hematologic parameters indicates that blood function and thus blood cell integrity were maintained throughout oxygenation. Already a versatile research tool, the now oxygenated FABRICA provides easy-to-implement, in vivo-like perfusion and stable oxygenation culture conditions in vitro semi-independently of one another, which means the bioreactor has the potential to serve as a platform for investigating the behavior of 3D tissue models (regardless of biofabrication method), performing drug toxicity-testing, and testing pharmaceutical efficacy/safety.

Avascular tumour growth models based on anomalous diffusion

In this study, we model avascular tumour growth in epithelial tissue. This can help us to understand that how an avascular tumour interacts with its microenvironment and what type of physical changes can be observed within the tumour spheroid before angiogenesis. This understanding is likely to assist in the development of better diagnostics, improved therapies, and prognostics. In biological systems, most of the diffusive processes are through cellular membranes which are porous in nature. Due to its porous nature, diffusion in biological systems are heterogeneous. The fractional diffusion equation is well suited to model heterogeneous biological systems, though most of the early studies did not use this fact. They described tumour growth with simple diffusion-based model. We have developed a spherical model based on simple diffusion initially, and then the model is upgraded with fractional diffusion equations to express the anomalous nature of biological system. In this study, two types of fractional models are developed: one of fixed order and the other of variable order. The memory formalism technique is also included in these anomalous diffusion models. These three models are investigated from phenomenological point view by measuring some parameters for characterizing avascular tumour growth over time. Tumour microenvironment is very complex in nature due to several concurrent molecular mechanisms. Diffusion with memory (fixed as well as variable) formation may be an oversimplified technique, and does not reflect the detailed view of the tumour microenvironment. However, it is found that all the models offer realistic and insightful information of the tumour microenvironment at the macroscopic level, and approximate well the physical phenomena. Also, it is observed that the anomalous diffusion based models offer a closer description to clinical facts than the simple model. As the simulation parameters get modified due to different biochemical and biophysical processes, the robustness of the model is determined. It is found that the anomalous diffusion models are moderately sensitive to the parameters.

Close Encounters of the Cell Kind: The Impact of Contact Inhibition on Tumour Growth and Cancer Models

Cancer is a complex phenomenon, and the sheer variation in behaviour across different types renders it difficult to ascertain underlying biological mechanisms. Experimental approaches frequently yield conflicting results for myriad reasons, and mathematical modelling of cancer is a vital tool to explore what we cannot readily measure, and ultimately improve treatment and prognosis. Like experiments, models are underpinned by certain biological assumptions, variation of which can lead to divergent predictions. An outstanding and important question concerns contact inhibition of proliferation (CIP), the observation that proliferation ceases when cells are spatially confined by their neighbours. CIP is a characteristic of many healthy adult tissues, but it remains unclear to which extent it holds in solid tumours, which exhibit regions of hyper-proliferation, and apparent breakdown of CIP. What precisely occurs in tumour tissue remains an open question, which mathematical modelling can help shed light on. In this perspective piece, we explore the implications of different hypotheses and available experimental evidence to elucidate the implications of these scenarios. We also outline how erroneous conclusions about the nature of tumour growth may be arrived at by looking selectively at biological data in isolation, and how this might be circumvented.

Oxygen in the Tumor Microenvironment: Mathematical and Numerical Modeling

There are many reasons to try to achieve a good grasp of the distribution of oxygen in the tumor microenvironment. The lack of oxygen - hypoxia - is a main actor in the evolution of tumors and in their growth and appears to be just as important in tumor invasion and metastasis. Mathematical models of the distribution of oxygen in tumors which are based on reaction-diffusion equations provide partial but qualitatively significant descriptions of the measured oxygen concentrations in the tumor microenvironment, especially when they incorporate important elements of the blood vessel network such as the blood vessel size and spatial distribution and the pulsation of local pressure due to blood circulation. Here, we review our mathematical and numerical approaches to the distribution of oxygen that yield insights both on the role of the distribution of blood vessel density and size and on the fluctuations of blood pressure.

Computational fluid dynamic analysis of bioprinted self-supporting perfused tissue models

Natural tissues are incorporated with vasculature, which is further integrated with a cardiovascular system responsible for driving perfusion of nutrient‐rich oxygenated blood through the vasculature to support cell metabolism within most cell‐dense tissues. Since scaffold‐free biofabricated tissues being developed into clinical implants, research models, and pharmaceutical testing platforms should similarly exhibit perfused tissue‐like structures, we generated a generalizable biofabrication method resulting in self‐supporting perfused (SSuPer) tissue constructs incorporated with perfusible microchannels and integrated with the modular FABRICA perfusion bioreactor. As proof of concept, we perfused an MLO‐A5 osteoblast‐based SSuPer tissue in the FABRICA. Although our resulting SSuPer tissue replicated vascularization and perfusion observed in situ, supported its own weight, and stained positively for mineral using Von Kossa staining, our in vitro results indicated that computational fluid dynamics (CFD) should be used to drive future construct design and flow application before further tissue biofabrication and perfusion. We built a CFD model of the SSuPer tissue integrated in the FABRICA and analyzed flow characteristics (net force, pressure distribution, shear stress, and oxygen distribution) through five SSuPer tissue microchannel patterns in two flow directions and at increasing flow rates. Important flow parameters include flow direction, fully developed flow, and tissue microchannel diameters matched and aligned with bioreactor flow channels. We observed that the SSuPer tissue platform is capable of providing direct perfusion to tissue constructs and proper culture conditions (oxygenation, with controllable shear and flow rates), indicating that our approach can be used to biofabricate tissue representing primary tissues and that we can model the system in silico.

Mechanism-Specific Pharmacodynamics of a Novel Complex-I Inhibitor Quantified by Imaging Reversal of Consumptive Hypoxia with [18F]FAZA PET In Vivo

Tumors lack a well-regulated vascular supply of O₂ and often fail to balance O₂ supply and demand. Net O₂ tension within many tumors may not only depend on O₂ delivery but also depend strongly on O₂ demand. Thus, tumor O₂ consumption rates may influence tumor hypoxia up to true anoxia. Recent reports have shown that many human tumors in vivo depend primarily on oxidative phosphorylation (OxPhos), not glycolysis, for energy generation, providing a driver for consumptive hypoxia and an exploitable vulnerability. In this regard, IACS-010759 is a novel high affinity inhibitor of OxPhos targeting mitochondrial complex-I that has recently completed a Phase-I clinical trial in leukemia. However, in solid tumors, the effective translation of OxPhos inhibitors requires methods to monitor pharmacodynamics in vivo. Herein, ¹⁸F-fluoroazomycin arabinoside ([¹⁸F]FAZA), a 2-nitroimidazole-based hypoxia PET imaging agent, was combined with a rigorous test-retest imaging method for non-invasive quantification of the reversal of consumptive hypoxia in vivo as a mechanism-specific pharmacodynamic (PD) biomarker of target engagement for IACS-010759. Neither cell death nor loss of perfusion could account for the IACS-010759-induced decrease in [¹⁸F]FAZA retention. Notably, in an OxPhos-reliant melanoma tumor, a titration curve using [¹⁸F]FAZA PET retention in vivo yielded an IC₅₀ for IACS-010759 (1.4 mg/kg) equivalent to analysis ex vivo. Pilot [¹⁸F]FAZA PET scans of a patient with grade IV glioblastoma yielded highly reproducible, high-contrast images of hypoxia in vivo as validated by CA-IX and GLUT-1 IHC ex vivo. Thus, [¹⁸F]FAZA PET imaging provided direct evidence for the presence of consumptive hypoxia in vivo, the capacity for targeted reversal of consumptive hypoxia through the inhibition of OxPhos, and a highly-coupled mechanism-specific PD biomarker ready for translation.

Perfluorocarbon nanodroplets can reoxygenate hypoxic tumors in vivo without carbogen breathing

Nanoscale perfluorocarbon (PFC) droplets have enormous potential as clinical theranostic agents. They are biocompatible and are currently used in vivo as contrast agents for a variety of medical imaging modalities, including ultrasound, computed tomography, photoacoustic and 19F-magnetic resonance imaging. PFC nanodroplets can also carry molecular and nanoparticulate drugs and be activated in situ by ultrasound or light for targeted therapy. Recently, there has been renewed interest in using PFC nanodroplets for hypoxic tumor reoxygenation towards radiosensitization based on the high oxygen solubility of PFCs. Previous studies showed that tumor oxygenation using PFC agents only occurs in combination with enhanced oxygen breathing. However, recent studies suggest that PFC agents that accumulate in solid tumors can contribute to radiosensitization, presumably due to tumor reoxygenation without enhanced oxygen breathing. In this study, we quantify the impact of oxygenation due to PFC nanodroplet accumulation in tumors alone in comparison with other reoxygenation methodologies, in particular, carbogen breathing. Methods: Lipid-stabilized, PFC (i.e., perfluorooctyl bromide, CF3(CF2)7Br, PFOB) nanoscale droplets were synthesized and evaluated in xenograft prostate (DU145) tumors in male mice. Biodistribution assessment of the nanodroplets was achieved using a fluorescent lipophilic indocarbocyanine dye label (i.e., DiI dye) on the lipid shell in combination with fluorescence imaging in mice (n≥3 per group). Hypoxia reduction in tumors was measured using PET imaging and a known hypoxia radiotracer, [18F]FAZA (n≥ 3 per group). Results: Lipid-stabilized nanoscale PFOB emulsions (mean diameter of ~250 nm), accumulated in the xenograft prostate tumors in mice 24 hours post-injection. In vivo PET imaging with [18F]FAZA showed that the accumulation of the PFOB nanodroplets in the tumor tissues alone significantly reduced tumor hypoxia, without enhanced oxygen (i.e., carbogen) breathing. This reoxygenation effect was found to be comparable with carbogen breathing alone. Conclusion: Accumulation of nanoscale PFOB agents in solid tumors alone successfully reoxygenated hypoxic tumors to levels comparable with carbogen breathing alone, an established tumor oxygenation method. This study confirms that PFC agents can be used to reoxygenate hypoxic tumors in addition to their current applications as multifunctional theranostic agents.

Oxygen diffusion in ellipsoidal tumour spheroids

Oxygen plays a central role in cellular metabolism, in both healthy and tumour tissue. The presence and concentration of molecular oxygen in tumours has a substantial effect on both radiotherapy response and tumour evolution, and as a result the oxygen micro-environment is an area of intense research interest. Multi-cellular tumour spheroids closely mimic real avascular tumours, and in particular they exhibit physiologically relevant heterogeneous oxygen distribution. This property has made them a vital part of in vitro experimentation. For ideal spheroids, their heterogeneous oxygen distributions can be predicted from theory, allowing determination of cellular oxygen consumption rate (OCR) and anoxic extent. However, experimental tumour spheroids often depart markedly from perfect sphericity. There has been little consideration of this reality. To date, the question of how far an ellipsoid can diverge from perfect sphericity before spherical assumptions break down remains unanswered. In this work, we derive equations governing oxygen distribution (and, more generally, nutrient and drug distribution) in both prolate and oblate tumour ellipsoids, and quantify the theoretical limits of the assumption that the spheroid is a perfect sphere. Results of this analysis yield new methods for quantifying OCR in ellipsoidal spheroids, and how this can be applied to markedly increase experimental throughput and quality.

Highly efficient triazolone/metal ion/polydopamine/MCM-41 sustained release system with pH sensitivity for pesticide delivery

MCM-41 was prepared through the sol-gel method and encapsulated by polydopamine (PDA) before being coordinated with metal ions to form a highly efficient sustained release system (M-PDA-MCM-41) for triazolone delivery. The characterization results confirmed the existence of the coordination bond between the PDA layer and triazolone through the bridge effect from metal ions, which enhanced the interaction between PDA-MCM-41 and triazolone. The adsorption capacity of Fe-PDA-MCM-41 increased up to 173 mg g-1, which was 160% more than that of MCM-41. The sustained release performance of M-PDA-MCM-41 in different pH values was investigated. Under the conditions of pH ≤7, the release speed of triazolone increased with pH decreasing, whereas its release speed in the weak base condition was slower than in the neutral condition. Therefore, the as-synthesized system showed significant pH-sensitivity in the sustained release process, indicating that the sustained release system can be well stored in the neutral or basic environment and activated in the acid environment. Their sustained release curves described by the Korsmeyer-Peppas equation at pH 7 showed the same behaviour, indicating that PDA decoration or metal ion coordination only increases the steric hindrance and the interaction between carrier and triazolone instead of changing the original structure of the pure MCM material in accordance with X-ray diffraction and Brunauer-Emmett-Teller analysis results.

The Relation Between Capillary Transit Times and Hemoglobin Saturation Heterogeneity. Part 1: Theoretical Models

Capillary dysfunction impairs oxygen supply to parenchymal cells and often occurs in Alzheimer's disease, diabetes and aging. Disturbed capillary flow patterns have been shown to limit the efficacy of oxygen extraction and can be quantified using capillary transit time heterogeneity (CTH). However, the transit time of red blood cells (RBCs) through the microvasculature is not a direct measure of their capacity for oxygen delivery. Here we examine the relation between CTH and capillary outflow saturation heterogeneity (COSH), which is the heterogeneity of blood oxygen content at the venous end of capillaries. Models for the evolution of hemoglobin saturation heterogeneity (HSH) in capillary networks were developed and validated using a computational model with moving RBCs. Two representative situations were selected: a Krogh cylinder geometry with heterogeneous hemoglobin saturation (HS) at the inflow, and a parallel array of four capillaries. The heterogeneity of HS after converging capillary bifurcations was found to exponentially decrease with a time scale of 0.15-0.21 s due to diffusive interaction between RBCs. Similarly, the HS difference between parallel capillaries also drops exponentially with a time scale of 0.12-0.19 s. These decay times are substantially smaller than measured RBC transit times and only weakly depend on the distance between microvessels. This work shows that diffusive interaction strongly reduces COSH on a small spatial scale. Therefore, we conclude that CTH influences COSH yet does not determine it. The second part of this study will focus on simulations in microvascular networks from the rodent cerebral cortex. Actual estimates of COSH and CTH will then be given.

Enhanced oxygen permeability in membrane-bottomed concave microwells for the formation of pancreatic islet spheroids

Oxygen availability is a critical factor in regulating cell viability that ultimately contributes to the normal morphogenesis and functionality of human tissues. Among various cell culture platforms, construction of 3D multicellular spheroids based on microwell arrays has been extensively applied to reconstitute in vitro human tissue models due to its precise control of tissue culture conditions as well as simple fabrication processes. However, an adequate supply of oxygen into the spheroidal cellular aggregation still remains one of the main challenges to producing healthy in vitro spheroidal tissue models. Here, we present a novel design for controlling the oxygen distribution in concave microwell arrays. We show that oxygen permeability into the microwell is tightly regulated by varying the poly-dimethylsiloxane (PDMS) bottom thickness of the concave microwells. Moreover, we validate the enhanced performance of the engineered microwell arrays by culturing non-proliferated primary rat pancreatic islet spheroids on varying bottom thickness from 10 μm to 1050 μm. Morphological and functional analyses performed on the pancreatic islet spheroids grown for 14 days prove the long-term stability, enhanced viability, and increased hormone secretion under the sufficient oxygen delivery conditions. We expect our results could provide knowledge on oxygen distribution in 3-dimensional spheroidal cell structures and critical design concept for tissue engineering applications.

STATEMENT OF SIGNIFICANCE

In this study, we present a noble design to control the oxygen distribution in concave microwell arrays for the formation of highly functional pancreatic islet spheroids by engineering the bottom of the microwells. Our new platform significantly enhanced oxygen permeability that turned out to improve cell viability and spheroidal functionality compared to the conventional thick-bottomed 3-D culture system. Therefore, we believe that this could be a promising medical biotechnology platform to further develop high-throughput tissue screening system as well as in vivo-mimicking customised 3-D tissue culture systems.

Hypoxia imaging and radiotherapy: bridging the resolution gap

Oxygen distribution is a major determinant of treatment success in radiotherapy, with well-oxygenated tumour regions responding by up to a factor of three relative to anoxic volumes. Conversely, tumour hypoxia is associated with treatment resistance and negative prognosis. Tumour oxygenation is highly heterogeneous and difficult to measure directly. The recent advent of functional hypoxia imaging modalities such as fluorine-18 fluoromisonidazole positron emission tomography have shown promise in non-invasively determining regions of low oxygen tension. This raises the prospect of selectively increasing dose to hypoxic subvolumes, a concept known as dose painting. Yet while this is a promising approach, oxygen-mediated radioresistance is inherently a multiscale problem, and there are still a number of substantial challenges that must be overcome if hypoxia dose painting is to be successfully implemented. Current imaging modalities are limited by the physics of such systems to have resolutions in the millimetre regime, whereas oxygen distribution varies over a micron scale, and treatment delivery is typically modulated on a centimetre scale. In this review, we examine the mechanistic basis and implications of the radiobiological oxygen effect, the factors influencing microscopic heterogeneity in tumour oxygenation and the consequent challenges in the interpretation of clinical hypoxia imaging (in particular fluorine-18 fluoromisonidazole positron emission tomography). We also discuss dose-painting approaches and outline challenges that must be addressed to improve this treatment paradigm.

Optimizing radiotherapy protocols using computer automata to model tumour cell death as a function of oxygen diffusion processes

The concept of hypofractionation is gaining momentum in radiation oncology centres, enabled by recent advances in radiotherapy apparatus. The gain of efficacy of this innovative treatment must be defined. We present a computer model based on translational murine data for in silico testing and optimization of various radiotherapy protocols with respect to tumour resistance and the microenvironment heterogeneity. This model combines automata approaches with image processing algorithms to simulate the cellular response of tumours exposed to ionizing radiation, modelling the alteration of oxygen permeabilization in blood vessels against repeated doses, and introducing mitotic catastrophe (as opposed to arbitrary delayed cell-death) as a means of modelling radiation-induced cell death. Published data describing cell death in vitro as well as tumour oxygenation in vivo are used to inform parameters. Our model is validated by comparing simulations to in vivo data obtained from the radiation treatment of mice transplanted with human prostate tumours. We then predict the efficacy of untested hypofractionation protocols, hypothesizing that tumour control can be optimized by adjusting daily radiation dosage as a function of the degree of hypoxia in the tumour environment. Further biological refinement of this tool will permit the rapid development of more sophisticated strategies for radiotherapy.

Sulforhodamine Nanothermometer for Multiparametric Fluorescence Lifetime Imaging Microscopy

Live cells function within narrow limits of physiological temperature (T) and O₂ and metabolite concentrations. We have designed a cell-permeable T-sensitive fluorescence lifetime-based nanoprobe based on lipophilic sulforhodamine, which stains 2D and 3D cell models, shows cytoplasmic localization, and has a robust response to T (∼0.037 ns/K). Subsequently, we evaluated the probe and fluorescence lifetime imaging microscopy (FLIM) technique for combined imaging of T and O₂ gradients in metabolically active cells. We found that in adherent 2D culture of HCT116 cells intracellular T and O₂ are close to ambient values. However, in 3D spheroid structures having size >200 μm, T and O₂ gradients become pronounced. These microgradients can be enhanced by treatment with mitochondrial uncouplers or dissipated by drug-induced disaggregation of the spheroids. Thus, we demonstrate the existence of local microgradients of T in 3D cell models and utility of combined imaging of O₂ and T.

Oxygen Partial Pressure Is a Rate-Limiting Parameter for Cell Proliferation in 3D Spheroids Grown in Physioxic Culture Condition

The in situ oxygen partial pressure in normal and tumor tissues is in the range of a few percent. Therefore, when studying cell growth in 3D culture systems, it is essential to consider how the physiological oxygen concentration, rather than the one in the ambient air, influences the proliferation parameters. Here, we investigated the effect of reducing oxygen partial pressure from 21% to 5% on cell proliferation rate and regionalization in a 3D tumor spheroid model. We found that 5% oxygen concentration strongly inhibited spheroid growth, changed the proliferation gradient and reduced the 50% In Depth Proliferation index (IDP50), compared with culture at 21% oxygen. We then modeled the oxygen partial pressure profiles using the experimental data generated by culturing spheroids in physioxic and normoxic conditions. Although hypoxia occurred at similar depth in spheroids grown in the two conditions, oxygen partial pressure was a major rate-limiting factor with a critical effect on cell proliferation rate and regionalization only in spheroids grown in physioxic condition and not in spheroids grown at atmospheric normoxia. Our findings strengthen the need to consider conducting experiment in physioxic conditions (i.e., tissue normoxia) for proper understanding of cancer cell biology and the evaluation of anticancer drugs in 3D culture systems.

The Role of Oxygen in Avascular Tumor Growth

The oxygen status of a tumor has significant clinical implications for treatment prognosis, with well-oxygenated subvolumes responding markedly better to radiotherapy than poorly supplied regions. Oxygen is essential for tumor growth, yet estimation of local oxygen distribution can be difficult to ascertain in situ, due to chaotic patterns of vasculature. It is possible to avoid this confounding influence by using avascular tumor models, such as tumor spheroids, a much better approximation of realistic tumor dynamics than monolayers, where oxygen supply can be described by diffusion alone. Similar to in situ tumours, spheroids exhibit an approximately sigmoidal growth curve, often approximated and fitted by logistic and Gompertzian sigmoid functions. These describe the basic rate of growth well, but do not offer an explicitly mechanistic explanation. This work examines the oxygen dynamics of spheroids and demonstrates that this growth can be derived mechanistically with cellular doubling time and oxygen consumption rate (OCR) being key parameters. The model is fitted to growth curves for a range of cell lines and derived values of OCR are validated using clinical measurement. Finally, we illustrate how changes in OCR due to gemcitabine treatment can be directly inferred using this model.

Estimating oxygen distribution from vasculature in three-dimensional tumour tissue

Regions of tissue which are well oxygenated respond better to radiotherapy than hypoxic regions by up to a factor of three. If these volumes could be accurately estimated, then it might be possible to selectively boost dose to radio-resistant regions, a concept known as dose-painting. While imaging modalities such as 18F-fluoromisonidazole positron emission tomography (PET) allow identification of hypoxic regions, they are intrinsically limited by the physics of such systems to the millimetre domain, whereas tumour oxygenation is known to vary over a micrometre scale. Mathematical modelling of microscopic tumour oxygen distribution therefore has the potential to complement and enhance macroscopic information derived from PET. In this work, we develop a general method of estimating oxygen distribution in three dimensions from a source vessel map. The method is applied analytically to line sources and quasi-linear idealized line source maps, and also applied to full three-dimensional vessel distributions through a kernel method and compared with oxygen distribution in tumour sections. The model outlined is flexible and stable, and can readily be applied to estimating likely microscopic oxygen distribution from any source geometry. We also investigate the problem of reconstructing three-dimensional oxygen maps from histological and confocal two-dimensional sections, concluding that two-dimensional histological sections are generally inadequate representations of the three-dimensional oxygen distribution.

Direct Measurements of Oxygen Gradients in Spheroid Culture System Using Electron Parametric Resonance Oximetry

Advanced in vitro culture from tissues of different origin includes three-dimensional (3D) organoid micro structures that may mimic conditions in vivo. One example of simple 3D culture is spheroids; ball shaped structures typically used as liver and tumour models. Oxygen is critically important in physiological processes, but is difficult to quantify in 3D culture: and the question arises, how small does a spheroid have to be to have minimal micro-environment formation? This question is of particular importance in the growing field of 3D based models for toxicological assessment. Here, we describe a simple non-invasive approach modified for the quantitative measurement and subsequent evaluation of oxygen gradients in spheroids developed from a non-malignant fish cell line (i.e. RTG-2 cells) using Electron Paramagnetic Resonance (EPR) oximetry. Sonication of the paramagnetic probe Lithium phthalocyanine (LiPc) allows for incorporation of probe particulates into spheroid during its formation. Spectra signal strength after incorporation of probe into spheroid indicated that a volume of 20 μl of probe (stock solution: 0.10 mg/mL) is sufficient to provide a strong spectra across a range of spheroid sizes. The addition of non-toxic probes (that do not produce or consume oxygen) report on oxygen diffusion throughout the spheroid as a function of size. We provide evidence supporting the use of this model over a range of initial cell seeding densities and spheroid sizes with the production of oxygen distribution as a function of these parameters. In our spheroid model, lower cell seeding densities (∼2,500 cells/spheroid) and absolute size (118±32 μm) allow control of factors such as pre-existing stresses (e.g. ∼ 2% normoxic/hypoxic interface) for more accurate measurement of treatment response. The applied methodology provides an elegant, widely applicable approach to directly characterize spheroid (and other organoid) cultures in biomedical and toxicological research.

A mechanistic investigation of the oxygen fixation hypothesis and oxygen enhancement ratio

The presence of oxygen in tumours has substantial impact on treatment outcome; relative to anoxic regions, well-oxygenated cells respond better to radiotherapy by a factor 2.5-3. This increased radio-response is known as the oxygen enhancement ratio. The oxygen effect is most commonly explained by the oxygen fixation hypothesis, which postulates that radical-induced DNA damage can be permanently 'fixed' by molecular oxygen, rendering DNA damage irreparable. While this oxygen effect is important in both existing therapy and for future modalities such a radiation dose-painting, the majority of existing mathematical models for oxygen enhancement are empirical rather than based on the underlying physics and radiochemistry. Here we propose a model of oxygen-enhanced damage from physical first principles, investigating factors that might influence the cell kill. This is fitted to a range of experimental oxygen curves from literature and shown to describe them well, yielding a single robust term for oxygen interaction obtained. The model also reveals a small thermal dependency exists but that this is unlikely to be exploitable.

System for measuring oxygen consumption rates of mammalian cells in static culture under hypoxic conditions

Estimating the oxygen consumption rates (OCRs) of mammalian cells in hypoxic environments is essential for designing and developing a three-dimensional (3-D) cell culture system. However, OCR measurements under hypoxic conditions are infrequently reported in the literature. Here, we developed a system for measuring OCRs at low oxygen levels. The system injects nitrogen gas into the environment and measures the oxygen concentration by an optical oxygen microsensor that consumes no oxygen. The developed system was applied to HepG2 cells in static culture. Specifically, we measured the spatial profiles of the local dissolved oxygen concentration in the medium, then estimated the OCRs of the cells. The OCRs, and also the pericellular oxygen concentrations, decreased nonlinearly as the oxygen partial pressure in the environment decreased from 19% to 1%. The OCRs also depended on the culture period and the matrix used for coating the dish surface. Using this system, we can precisely estimate the OCRs of various cell types under environments that mimic 3-D culture conditions, contributing crucial data for an efficient 3-D culture system design.

Mathematical Model of Oxygen Transport in Tuberculosis Granulomas

Pulmonary granulomas--the hallmark of Mycobacterium tuberculosis (MTB) infection--are dense cellular lesions that often feature regions of hypoxia and necrosis, partially due to limited transport of oxygen. Low oxygen in granulomas can impair the host immune response, while MTB are able to adapt and persist in hypoxic environments. Here, we used a physiologically based mathematical model of oxygen diffusion and consumption to calculate oxygen profiles within the granuloma, assuming Michaelis-Menten kinetics. An approximate analytical solution--using a priori and newly estimated parameters from experimental data in a rabbit model of tuberculosis--was able to predict the size of hypoxic and necrotic regions in agreement with experimental results from the animal model. Such quantitative understanding of transport limitations can inform future tuberculosis therapeutic strategies that may include adjunct host-directed therapies that facilitate oxygen and drug delivery for more effective treatment.