Theoretical simulation of oxygen transport to tumors by three-dimensional networks of microvessels

Spatially-resolved pharmacokinetic/pharmacodynamic modelling of bystander effects of a nitrochloromethylbenzindoline hypoxia-activated prodrug

PURPOSE

Hypoxia-activated prodrugs (HAPs) have the potential for eliminating chemo- and radiation-resistant hypoxic tumour cells, but their activity is often compromised by limited penetration into hypoxic zones. Nitrochloromethylbenzindoline (nitroCBI) HAPs are reduced in hypoxic cells to highly cytotoxic DNA minor groove alkylating aminoCBI metabolites. In this study, we investigate whether a lead nitroCBI, SN30548, generates a significant bystander effect through the diffusion of its aminoCBI metabolite and whether this compensates for any diffusion limitations of the prodrug in tumour tissue.

METHODS

Metabolism and uptake of the nitroCBI in oxic and anoxic cells, and diffusion through multicellular layer cultures, was characterised by LC-MS/MS. To quantify bystander effects, clonogenic cell killing of HCT116 cells was assessed in multicellular spheroid co-cultures comprising cells transfected with cytochrome P450 oxidoreductase (POR) or E. coli nitroreductase NfsA. Spatially-resolved pharmacokinetic/pharmacodynamic (PK/PD) models, parameterised by the above measurements, were developed for spheroids and tumours using agent-based and Green's function modelling, respectively.

RESULTS

NitroCBI was reduced to aminoCBI by POR under anoxia and by NfsA under oxia, and was the only significant cytotoxic metabolite in both cases. In spheroid co-cultures comprising 30% NfsA-expressing cells, non-metabolising cells were as sensitive as the NfsA cells, demonstrating a marked bystander effect. Agent-based PK/PD models provided good prediction of cytotoxicity in spheroids, while use of the same parameters in a Green's function model for a tumour microregion demonstrated that local diffusion of aminoCBI overcomes the penetration limitation of the prodrug.

CONCLUSIONS

The nitroCBI HAP SN30548 generates a highly efficient bystander effect through local diffusion of its active metabolite in tumour tissue.

Effects of impaired microvascular flow regulation on metabolism-perfusion matching and organ function

Impaired tissue oxygen delivery is a major cause of organ damage and failure in critically ill patients, which can occur even when systemic parameters, including cardiac output and arterial hemoglobin saturation, are close to normal. This review addresses oxygen transport mechanisms at the microcirculatory scale, and how hypoxia may occur in spite of adequate convective oxygen supply. The structure of the microcirculation is intrinsically heterogeneous, with wide variations in vessel diameters and flow pathway lengths, and consequently also in blood flow rates and oxygen levels. The dynamic processes of structural adaptation and flow regulation continually adjust microvessel diameters to compensate for heterogeneity, redistributing flow according to metabolic needs to ensure adequate tissue oxygenation. A key role in flow regulation is played by conducted responses, which are generated and propagated by endothelial cells and signal upstream arterioles to dilate in response to local hypoxia. Several pathophysiological conditions can impair local flow regulation, causing hypoxia and tissue damage leading to organ failure. Therapeutic measures targeted to systemic parameters may not address or may even worsen tissue oxygenation at the microvascular level. Restoration of tissue oxygenation in critically ill patients may depend on restoration of endothelial cell function, including conducted responses.

Dual-emissive, oxygen-sensing boron nanoparticles quantify oxygen consumption rate in breast cancer cells

SIGNIFICANCE

Decreasing the oxygen consumption rate (OCR) of tumor cells is a powerful method for ameliorating tumor hypoxia. However, quantifying the change in OCR is challenging in complex experimental systems.

AIM

We present a method for quantifying the OCR of two tumor cell lines using oxygen-sensitive dual-emissive boron nanoparticles (BNPs). We hypothesize that our BNP results are equivalent to the standard Seahorse assay.

APPROACH

We quantified the spectral emissions of the BNP and accounted for external oxygen diffusion to quantify OCR over 24 h. The BNP-computed OCR of two breast cancer cell lines, E0771 and 4T07, were compared with their respective Seahorse assays. Both cell lines were also irradiated to quantify radiation-induced changes in the OCR.

RESULTS

Using a Bland-Altman analysis, our BNPs OCR was equivalent to the standard Seahorse assay. Moreover, in an additional experiment in which we irradiated the cells at their 50% survival fraction, the BNPs were sensitive enough to quantify 24% reduction in OCR after irradiation.

CONCLUSIONS

Our results conclude that the BNPs are a viable alternative to the Seahorse assay for quantifying the OCR in cells. The Bland-Altman analysis showed that these two methods result in equivalent OCR measurements. Future studies will extend the OCR measurements to complex systems including 3D cultures and in vivo models, in which OCR measurements cannot currently be made.

Morphometric, Hemodynamic, and Biomechanical Factors Influencing Blood Flow and Oxygen Concentration in the Human Lamina Cribrosa

Purpose

We developed a combined biomechanical and hemodynamic model of the human eye to estimate blood flow and oxygen concentration within the lamina cribrosa (LC) and rank the factors that influence LC oxygen concentration.

Methods

We generated 5000 finite-element eye models with detailed microcapillary networks of the LC and computed the oxygen concentration of the lamina retinal ganglion cell axons. For each model, we varied the intraocular pressure (IOP) from 10 mm Hg to 55 mm Hg in 5-mm Hg increments, the cerebrospinal fluid pressure (13 ± 2 mm Hg), cup depth (0.2 ± 0.1 mm), scleral stiffness (±20% of the mean values), LC stiffness (0.41 ± 0.2 MPa), LC radius (1.2 ± 0.12 mm), average LC pore size (5400 ± 2400 µm²), the microcapillary arrangement (radial, isotropic, or circumferential), and perfusion pressure (50 ± 9 mm Hg). Blood flow was assumed to originate from the LC periphery and drain via the central retinal vein. Finally, we performed linear regressions to rank the influence of each factor on the LC tissue oxygen concentration.

Results

LC radius and perfusion pressure were the most important factors in influencing the oxygen concentration of the LC. IOP was another important parameter, and eyes with higher IOP had higher compressive strain and slightly lower oxygen concentration. In general, superior–inferior regions of the LC had significantly lower oxygen concentration than the nasal–temporal regions, resulting in an hourglass pattern of oxygen deficiency.

Conclusions

To the best of our knowledge, this study is the first to implement a comprehensive hemodynamical model of the eye that accounts for the biomechanical forces and morphological parameters of the LC. The results provide further insight into the possible relationship of biomechanical and vascular pathways leading to ischemia-induced optic neuropathy.

A Systems Pharmacology Model for Drug Delivery to Solid Tumors by Antibody-Drug Conjugates: Implications for Bystander Effects

Antibody-drug conjugates (ADCs) are cancer drugs composed of a humanized antibody linked to a cytotoxic payload, allowing preferential release of payload in cancer cells expressing the antibody-targeted antigen. Here, a systems pharmacology model is used to simulate ADC transport from blood to tumor tissue and ADC uptake by tumor cells. The model includes effects of spatial gradients in drug concentration in a three-dimensional network of tumor blood vessels with realistic geometry and accounts for diffusion of ADC in the tumor extracellular space, binding to antigen, internalization, intracellular processing, and payload efflux from cells. Cells that process an internalized ADC-antigen complex may release payload that can be taken up by other "bystander" cells. Such bystander effects are included in the model. The model is used to simulate conditions in previous experiments, showing good agreement with experimental results. Simulations are used to analyze the relationship of bystander effects to payload properties and single-dose administrations. The model indicates that exposure of payload to cells distant from vessels is sensitive to the free payload diffusivity in the extracellular space. When antigen expression is heterogeneous, the model indicates that the amount of payload accumulating in non-antigen-expressing cells increases linearly with dose but depends only weakly on the percentage of antigen-expressing cells. The model provides an integrated mechanistic framework for understanding the effects of spatial gradients on drug distribution using ADCs and for designing ADCs to achieve more effective payload distribution in solid tumors, thereby increasing the therapeutic index of the ADC.

An approach for vascular tumor growth based on a hybrid embedded/homogenized treatment of the vasculature within a multiphase porous medium model

The aim of this work is to develop a novel computational approach to facilitate the modeling of angiogenesis during tumor growth. The preexisting vasculature is modeled as a 1D inclusion and embedded into the 3D tissue through a suitable coupling method, which allows for nonmatching meshes in 1D and 3D domain. The neovasculature, which is formed during angiogenesis, is represented in a homogenized way as a phase in our multiphase porous medium system. This splitting of models is motivated by the highly complex morphology, physiology, and flow patterns in the neovasculature, which are challenging and computationally expensive to resolve with a discrete, 1D angiogenesis and blood flow model. Moreover, it is questionable if a discrete representation generates any useful additional insight. By contrast, our model may be classified as a hybrid vascular multiphase tumor growth model in the sense that a discrete, 1D representation of the preexisting vasculature is coupled with a continuum model describing angiogenesis. It is based on an originally avascular model which has been derived via the thermodynamically constrained averaging theory. The new model enables us to study mass transport from the preexisting vasculature into the neovasculature and tumor tissue. We show by means of several illustrative examples that it is indeed capable of reproducing important aspects of vascular tumor growth phenomenologically.

Implementation of the Chick Chorioallantoic Membrane (CAM) Model in Radiation Biology and Experimental Radiation Oncology Research

Radiotherapy (RT) is part of standard cancer treatment. Innovations in treatment planning and increased precision in dose delivery have significantly improved the therapeutic gain of radiotherapy but are reaching their limits due to biologic constraints. Thus, a better understanding of the complex local and systemic responses to RT and of the biological mechanisms causing treatment success or failure is required if we aim to define novel targets for biological therapy optimization. Moreover, optimal treatment schedules and prognostic biomarkers have to be defined for assigning patients to the best treatment option. The complexity of the tumor environment and of the radiation response requires extensive in vivo experiments for the validation of such treatments. So far in vivo investigations have mostly been performed in time- and cost-intensive murine models. Here we propose the implementation of the chick chorioallantoic membrane (CAM) model as a fast, cost-efficient model for semi high-throughput preclinical in vivo screening of the modulation of the radiation effects by molecularly targeted drugs. This review provides a comprehensive overview on the application spectrum, advantages and limitations of the CAM assay and summarizes current knowledge of its applicability for cancer research with special focus on research in radiation biology and experimental radiation oncology.

An Intratumor Pharmacokinetic/Pharmacodynamic Model for the Hypoxia-Activated Prodrug Evofosfamide (TH-302): Monotherapy Activity is Not Dependent on a Bystander Effect

Tumor hypoxia contributes to resistance to anticancer therapies. Hypoxia-activated prodrugs (HAPs) selectively target hypoxic cells and their activity can extend to well-oxygenated areas of tumors via diffusion of active metabolites. This type of bystander effect has been suggested to be responsible for the single agent activity of the clinical-stage HAP evofosfamide (TH-302) but direct evidence is lacking. To dissect the contribution of bystander effects to TH-302 activity, we implemented a Green's function pharmacokinetic (PK) model to simulate the spatial distribution of O₂, TH-302 and its cytotoxic metabolites, bromo-isophosphoramide mustard (Br-IPM) and its dichloro derivative isophosphoramide mustard (IPM), in two digitized tumor microvascular networks. The model was parameterized from literature and experimentally, including measurement of diffusion coefficients of TH-302 and its metabolites in multicellular layer cultures. The latter studies demonstrate that Br-IPM and IPM cannot diffuse significantly from the cells in which they are generated, although evidence was obtained for diffusion of the hydroxylamine metabolite of TH-302. The spatially resolved PK model was linked to a pharmacodynamic (PD) model that describes cell killing probability at each point in the tumor microregion as a function of Br-IPM and IPM exposure. The resulting PK/PD model accurately predicted previously reported monotherapy activity of TH-302 in H460 tumors, without invoking a bystander effect, demonstrating that the notable single agent activity of TH-302 in tumors can be accounted for by significant bioreductive activation of TH-302 even in oxic regions, driven by the high plasma concentrations achievable with this well-tolerated prodrug.

Differential Impact of Flavonoids on Redox Modulation, Bioenergetics, and Cell Signaling in Normal and Tumor Cells: A Comprehensive Review

SIGNIFICANCE

Flavonoids can interact with multiple molecular targets to elicit their cellular effects, leading to changes in signal transduction, gene expression, and/or metabolism, which can, subsequently, affect the entire cell and organism. Immortalized cell lines, derived from tumors, are routinely employed as a surrogate for mechanistic studies, with the results extrapolated to tissues in vivo. Recent Advances: We review the activities of selected flavonoids on cultured tumor cells derived from various tissues in comparison to corresponding primary cells or tissues in vivo, mainly using quercetin and flavanols (epicatechin and (-)-epigallocatechin gallate) as exemplars. Several studies have indicated that flavonoids could retard cancer progression in vivo in animal models as well as in tumor cell models.

CRITICAL ISSUES

Extrapolation from in vitro and animal models to humans is not straightforward given both the extensive conjugation and complex microbiota-dependent metabolism of flavonoids after consumption, as well as the heterogeneous metabolism of different tumors.

FUTURE DIRECTIONS

Comparison of data from studies on primary cells or in vivo are essential not only to validate results obtained from cultured cell models, but also to highlight whether any differences may be further exploited in the clinical setting for chemoprevention. Tumor cell models can provide a useful mechanistic tool to study the effects of flavonoids, provided that the limitations of each model are understood and taken into account in interpretation of the data.

An agent-based model for drug-radiation interactions in the tumour microenvironment: Hypoxia-activated prodrug SN30000 in multicellular tumour spheroids

Multicellular tumour spheroids capture many characteristics of human tumour microenvironments, including hypoxia, and represent an experimentally tractable in vitro model for studying interactions between radiotherapy and anticancer drugs. However, interpreting spheroid data is challenging because of limited ability to observe cell fate within spheroids dynamically. To overcome this limitation, we have developed a hybrid continuum/agent-based model (ABM) for HCT116 tumour spheroids, parameterised using experimental models (monolayers and multilayers) in which reaction and diffusion can be measured directly. In the ABM, cell fate is simulated as a function of local oxygen, glucose and drug concentrations, determined by solving diffusion equations and intracellular reactions. The model is lattice-based, with cells occupying discrete locations on a 3D grid embedded within a coarser grid that encompasses the culture medium; separate solvers are employed for each grid. The generated concentration fields account for depletion in the medium and specify concentration-time profiles within the spheroid. Cell growth and survival are determined by intracellular oxygen and glucose concentrations, the latter based on direct measurement of glucose diffusion/reaction (in multilayers) for the first time. The ABM reproduces known features of spheroids including overall growth rate, its oxygen and glucose dependence, peripheral cell proliferation, central hypoxia and necrosis. We extended the ABM to describe in detail the hypoxia-dependent interaction between ionising radiation and a hypoxia-activated prodrug (SN30000), again using experimentally determined parameters; the model accurately simulated clonogenic cell killing in spheroids, while inclusion of reversible cell cycle delay was required to account for the marked spheroid growth delay after combined radiation and SN30000. This ABM of spheroid growth and response exemplifies the utility of integrating computational and experimental tools for investigating radiation/drug interactions, and highlights the critical importance of understanding oxygen, glucose and drug concentration gradients in interpreting activity of therapeutic agents in spheroid models.

Studies in 3D cultures, notably multicellular tumour spheroids that mimic many features of solid tumours, have great potential for speeding up anticancer drug development. However the increased complexity of 3D cultures makes interpretation of experiments more difficult. We have developed a hybrid continuum/agent-based mathematical model, validated by experiments, to aid interpretation of spheroid experiments in developing drugs designed to eliminate radiation-resistant hypoxic cells. This model includes key features of the tumour microenvironment including oxygen and glucose transport and regions of hypoxia where the cells are resistant to radiation, but sensitive to hypoxia-activated prodrugs such as SN30000. This enables us to predict the growth and cell response in untreated spheroids and compare the results to spheroids treated with radiation and SN30000. We demonstrate good prediction of cellular responses in spheroids treated with radiation and SN30000 and good agreement with spheroid regrowth after treatment when additional effects of cellular growth delay are added. This demonstrates that the modelling approach has potential to improve interpretation of experimental investigations of drug and radiation combinations.

Computational modelling of drug delivery to solid tumour: Understanding the interplay between chemotherapeutics and biological system for optimised delivery systems

Drug delivery to solid tumour involves multiple physiological, biochemical and biophysical processes taking place across a wide range of length and time scales. The therapeutic efficacy of anticancer drugs is influenced by the complex interplays among the intrinsic properties of tumours, biophysical aspects of drug transport and cellular uptake. Mathematical and computational modelling allows for a well-controlled study on the individual and combined effects of a wide range of parameters on drug transport and therapeutic efficacy, which would not be possible or economically viable through experimental means. A wide spectrum of mathematical models has been developed for the simulation of drug transport and delivery in solid tumours, including PK/PD-based compartmental models, microscopic and macroscopic transport models, and molecular dynamics drug loading and release models. These models have been used as a tool to identify the limiting factors and for optimal design of efficient drug delivery systems. This article gives an overview of the currently available computational models for drug transport in solid tumours, together with their applications to novel drug delivery systems, such as nanoparticle-mediated drug delivery and convection-enhanced delivery.

Roles of PFKFB3 in cancer

The understanding of 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 (PFK-2/FBPase 3, PFKFB3) has advanced considerably since its initial identification in human macrophages in the mid-1990s. As a vital regulator of glycolysis, accumulating studies have suggested that PFKFB3 is associated with many aspects of cancer, including carcinogenesis, cancer cell proliferation, vessel aggressiveness, drug resistance and tumor microenvironment. In this review, we summarize current knowledge of PFKFB3 regulation by several signal pathways and its function in cancer development in different cell types in cancer tissues. Ubiquitous PFKFB3 has emerged as a potential target for anti-neoplastic therapy.

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.

A computational study of cancer hyperthermia based on vascular magnetic nanoconstructs

The application of hyperthermia to cancer treatment is studied using a novel model arising from the fundamental principles of flow, mass and heat transport in biological tissues. The model is defined at the scale of the tumour microenvironment and an advanced computational scheme called the embedded multiscale method is adopted to solve the governing equations. More precisely, this approach involves modelling capillaries as one-dimensional channels carrying flow, and special mathematical operators are used to model their interaction with the surrounding tissue. The proposed computational scheme is used to analyse hyperthermic treatment of cancer based on systemically injected vascular magnetic nanoconstructs carrying super-paramagnetic iron oxide nanoparticles. An alternating magnetic field is used to excite the nanoconstructs and generate localized heat within the tissue. The proposed model is particularly adequate for this application, since it has a unique capability of incorporating microvasculature configurations based on physiological data combined with coupled capillary flow, interstitial filtration and heat transfer. A virtual tumour model is initialized and the spatio-temporal distribution of nanoconstructs in the vascular network is analysed. In particular, for a reference iron oxide concentration, temperature maps of several different hypothesized treatments are generated in the virtual tumour model. The observations of the current study might in future guide the design of more efficient treatments for cancer hyperthermia.

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.

Modelling mass and heat transfer in nano-based cancer hyperthermia

We derive a sophisticated mathematical model for coupled heat and mass transport in the tumour microenvironment and we apply it to study nanoparticle delivery and hyperthermic treatment of cancer. The model has the unique ability of combining the following features: (i) realistic vasculature; (ii) coupled capillary and interstitial flow; (iii) coupled capillary and interstitial mass transfer applied to nanoparticles; and (iv) coupled capillary and interstitial heat transfer, which are the fundamental mechanisms governing nano-based hyperthermic treatment. This is an improvement with respect to previous modelling approaches, where the effect of blood perfusion on heat transfer is modelled in a spatially averaged form. We analyse the time evolution and the spatial distribution of particles and temperature in a tumour mass treated with superparamagnetic nanoparticles excited by an alternating magnetic field. By means of numerical experiments, we synthesize scaling laws that illustrate how nano-based hyperthermia depends on tumour size and vascularity. In particular, we identify two distinct mechanisms that regulate the distribution of particle and temperature, which are characterized by perfusion and diffusion, respectively.

A computational model of drug delivery through microcirculation to compare different tumor treatments

Starting from the fundamental laws of filtration and transport in biological tissues, we develop a computational model to capture the interplay between blood perfusion, fluid exchange with the interstitial volume, mass transport in the capillary bed, through the capillary walls and into the surrounding tissue. These phenomena are accounted at the microscale level, where capillaries and interstitial volume are viewed as two separate regions. The capillaries are described as a network of vessels carrying blood flow. We apply the model to study drug delivery to tumors. The model can be adapted to compare various treatment options. In particular, we consider delivery using drug bolus injection and nanoparticle injection into the blood stream. The computational approach is suitable for a systematic quantification of the treatment performance, enabling the analysis of interstitial drug concentration levels, metabolization rates and cell surviving fractions. Our study suggests that for the treatment based on bolus injection, the drug dose is not optimally delivered to the tumor interstitial volume. Using nanoparticles as intermediate drug carriers overrides the shortcomings of the previous delivery approach. This work shows that the proposed theoretical and computational framework represents a promising tool to compare the efficacy of different cancer treatments.

The Role of Bystander Effects in the Antitumor Activity of the Hypoxia-Activated Prodrug PR-104

Activation of prodrugs in tumors (e.g., by bioreduction in hypoxic zones) has the potential to generate active metabolites that can diffuse within the tumor microenvironment. Such “bystander effects” may offset spatial heterogeneity in prodrug activation but the relative importance of this effect is not understood. Here, we quantify the contribution of bystander effects to antitumor activity for the first time, by developing a spatially resolved pharmacokinetic/pharmacodynamic (SR-PK/PD) model for PR-104, a phosphate ester pre-prodrug that is converted systemically to the hypoxia-activated prodrug PR-104A. Using Green’s function methods we calculated concentrations of oxygen, PR-104A and its active metabolites, and resultant cell killing, at each point of a mapped three-dimensional tumor microregion. Model parameters were determined in vitro, using single cell suspensions to determine relationships between PR-104A metabolism and clonogenic cell killing, and multicellular layer (MCL) cultures to measure tissue diffusion coefficients. LC-MS/MS detection of active metabolites in the extracellular medium following exposure of anoxic single cell suspensions and MCLs to PR-104A confirmed that metabolites can diffuse out of cells and through a tissue-like environment. The SR-PK/PD model estimated that bystander effects contribute 30 and 50% of PR-104 activity in SiHa and HCT116 tumors, respectively. Testing the model by modulating PR-104A-activating reductases and hypoxia in tumor xenografts showed overall clonogenic killing broadly consistent with model predictions. Overall, our data suggest that bystander effects are important in PR-104 antitumor activity, although their reach may be limited by macroregional heterogeneity in hypoxia and reductase expression in tumors. The reported computational and experimental techniques are broadly applicable to all targeted anticancer prodrugs and could be used to identify strategies for rational prodrug optimization.

Part III. Molecular changes induced by high nitric oxide adaptation in human breast cancer cell line BT-20 (BT-20-HNO): a switch from aerobic to anaerobic metabolism

Nutrient deprivation and reactive oxygen species (ROS) play an important role in breast cancer mitochondrial adaptation. Adaptations to these conditions allow cells to survive in the stressful microenvironment of the tumor bed. This study is directed at defining the consequences of High Nitric Oxide (HNO) exposure to mitochondria in human breast cancer cells. The breast cancer cell line BT-20 (parent) was adapted to HNO as previously reported, resulting in the BT-20-HNO cell line. Both cell lines were analyzed by a variety of methods including MTT, LDH leakage assay, DNA sequencing, and Western blot analysis. The LDH assay and the gene chip data showed that BT-20-HNO was more prone to use the glycolytic pathway than the parent cell line. The BT-20-HNO cells were also more resistant to the apoptotic inducing agent salinomycin, which suggests that p53 may be mutated in these cells. Polymerase chain reaction (PCR) followed by DNA sequencing of the p53 gene showed that it was, in fact, mutated at the DNA-binding site (L194F). Western blot analysis showed that p53 was significantly upregulated in these cells. These results suggest that free radicals, such as nitric oxide (NO), pressure human breast tumor cells to acquire an aggressive phenotype and resistance to apoptosis. These data collectively provide a mechanism by which the dysregulation of ROS in the mitochondria of breast cancer cells can result in DNA damage.

Uptake rate of cationic mitochondrial inhibitor MKT-077 determines cellular oxygen consumption change in carcinoma cells

OBJECTIVE

Since tumor radiation response is oxygen-dependent, radiosensitivity can be enhanced by increasing tumor oxygenation. Theoretically, inhibiting cellular oxygen consumption is the most efficient way to increase oxygen levels. The cationic, rhodacyanine dye-analog MKT-077 inhibits mitochondrial respiration and could be an effective metabolic inhibitor. However, the relationship between cellular MKT-077 uptake and metabolic inhibition is unknown. We hypothesized that rat and human mammary carcinoma cells would take up MKT-077, causing a decrease in oxygen metabolism related to drug uptake.

METHODS

R3230Ac rat breast adenocarcinoma cells were exposed to MKT-077. Cellular MKT-077 concentration was quantified using spectroscopy, and oxygen consumption was measured using polarographic electrodes. MKT-077 uptake kinetics were modeled by accounting for uptake due to both the concentration and potential gradients across the plasma and mitochondrial membranes. These kinetic parameters were used to model the relationship between MKT-077 uptake and metabolic inhibition. MKT-077-induced changes in oxygen consumption were also characterized in MDA-MB231 human breast carcinoma cells.

RESULTS

Cells took up MKT-077 with a time constant of ∼1 hr, and modeling showed that over 90% of intracellular MKT-077 was bound or sequestered, likely by the mitochondria. The uptake resulted in a rapid decrease in oxygen consumption, with a time constant of ∼30 minutes. Surprisingly the change in oxygen consumption was proportional to uptake rate, not cellular concentration. MKT-077 proved a potent metabolic inhibitor, with dose-dependent decreases of 45-73% (p = 0.003).

CONCLUSIONS

MKT-077 caused an uptake rate-dependent decrease in cellular metabolism, suggesting potential efficacy for increasing tumor oxygen levels and radiosensitivity in vivo.

19F MRI oximetry: simulation of perfluorocarbon distribution impact

In (19)F MRI oximetry, a method used to image tumour hypoxia, perfluorocarbons serve as oxygenation markers. The goal of this study is to evaluate the impact of perfluorocarbon distribution and concentration in (19)F MRI oximetry through a computer simulation. The simulation studies the correspondence between (19)F measured (pO(FNMR)(2)) and actual tissue oxygen tension (pO(2)) for several tissue perfluorocarbon distributions. For this, a Krogh tissue model is implemented which incorporates the presence of perfluorocarbons in blood and tissue. That is, in tissue the perfluorocarbons are distributed homogeneously according to Gaussian diffusion profiles, or the perfluorocarbons are concentrated in the capillary wall. Using these distributions, the oxygen tension in the simulation volume is calculated. The simulated mean oxygen tension is then compared with pO(FNMR)(2), the (19)F MRI-based measure of pO(2) and with pO(0)(2), pO(2) in the absence of perfluorocarbons. The agreement between pO(FNMR)(2) and actual pO(2) is influenced by vascular density and perfluorocarbon distribution. The presence of perfluorocarbons generally gives rise to a pO(2) increase in tissue. This effect is enhanced when perfluorocarbons are also present in blood. Only the homogeneous perfluorocarbon distribution in tissue with no perfluorocarbons in blood guarantees small deviations of pO(FNMR)(2) from pO(2). Hence, perfluorocarbon distribution in tissue and blood has a serious impact on the reliability of (19)F MRI-based measures of oxygen tension. In addition, the presence of perfluorocarbons influences the actual oxygen tension. This finding may be of great importance for further development of (19)F MRI oximetry.

Emerging metabolic targets in cancer therapy

Cancer cells are different from normal cells in their metabolic properties. Normal cells mostly rely on mitochondrial oxidative phosphorylation to produce energy. In contrast, cancer cells depend mostly on glycolysis, the aerobic breakdown of glucose into ATP. This altered energy dependency is known as the "Warburg effect" and is a hallmark of cancer cells. In recent years, investigating the metabolic changes within cancer cells has been a rapidly growing area. Emerging evidence shows that oncogenes that drive the cancer-promoting signals also drive the altered metabolism. Although the exact mechanisms underlying the Warburg effect are unclear, the existing evidence suggests that increased glycolysis plays an important role in support malignant behavior of cancer cells. A thorough understanding of the unique metabolism of cancer cells will help to design of more effective drugs targeting metabolic pathways, which will greatly impact the capacity to effectively treat cancer patients. Here we provide an overview of the current understanding of the Warburg effect upon tumor cell growth and survival, and discussion on the potential metabolic targets for cancer therapy.

An oxygen-consuming phantom simulating perfused tissue to explore oxygen dynamics and (19)F MRI oximetry

OBJECTIVE

This study presents a reproducible phantom which mimics oxygen-consuming tissue and can be used for the validation of (19)F MRI oximetry.

MATERIALS AND METHODS

The phantom consists of a haemodialysis filter of which the outer compartment is filled with a gelatin matrix containing viable yeast cells. Perfluorocarbon emulsions can be added to the gelatin matrix to simulate sequestered perfluorocarbons. A blood-substituting perfluorocarbon fluid is pumped through the lumen of the fibres in the filter. (19)F relaxometry MRI is performed with a fast 2D Look-Locker imaging sequence on a clinical 3T scanner.

RESULTS

Acute and perfusion-related hypoxia were simulated and imaged spatially and temporally using the phantom.

CONCLUSIONS

The presented experimental setup can be used to simulate oxygen consumption by somatic cells in vivo and for validating computational biophysical models of hypoxia, as measured with (19)F MRI oximetry.

Oxygen consumption can regulate the growth of tumors, a new perspective on the Warburg effect

Background

The unique metabolism of tumors was described many years ago by Otto Warburg, who identified tumor cells with increased glycolysis and decreased mitochondrial activity. However, “aerobic glycolysis” generates fewer ATP per glucose molecule than mitochondrial oxidative phosphorylation, so in terms of energy production, it is unclear how increasing a less efficient process provides tumors with a growth advantage.

Methods/Findings

We carried out a screen for loss of genetic elements in pancreatic tumor cells that accelerated their growth as tumors, and identified mitochondrial ribosomal protein L28 (MRPL28). Knockdown of MRPL28 in these cells decreased mitochondrial activity, and increased glycolysis, but paradoxically, decreased cellular growth in vitro. Following Warburg's observations, this mutation causes decreased mitochondrial function, compensatory increase in glycolysis and accelerated growth in vivo. Likewise, knockdown of either mitochondrial ribosomal protein L12 (MRPL12) or cytochrome oxidase had a similar effect. Conversely, expression of the mitochondrial uncoupling protein 1 (UCP1) increased oxygen consumption and decreased tumor growth. Finally, treatment of tumor bearing animals with dichloroacetate (DCA) increased pyruvate consumption in the mitochondria, increased total oxygen consumption, increased tumor hypoxia and slowed tumor growth.

Conclusions

We interpret these findings to show that non-oncogenic genetic changes that alter mitochondrial metabolism can regulate tumor growth through modulation of the consumption of oxygen, which appears to be a rate limiting substrate for tumor proliferation.

Dynamic contrast-enhanced magnetic resonance imaging as a predictor of clinical outcome in canine spontaneous soft tissue sarcomas treated with thermoradiotherapy

PURPOSE

This study tests whether dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) parameters obtained from canine patients with soft tissue sarcomas, treated with hyperthermia and radiotherapy, are predictive of therapeutic outcome.

EXPERIMENTAL DESIGN

Thirty-seven dogs with soft tissue sarcomas had DCE-MRI done before and following the first hyperthermia. Signal enhancement for tumor and reference muscle were fitted empirically, yielding a washin/washout rate for the contrast agent and tumor area under the signal enhancement curve (AUC) calculated from 0 to 60 seconds, 90 seconds, and the time of maximal enhancement in the reference muscle. These parameters were then compared with local tumor control, metastasis-free survival, and overall survival.

RESULTS

Pretherapy rate of contrast agent washout was positively predictive of improved overall and metastasis-free survival with hazard ratio of 0.67 (P = 0.015) and 0.68 (P = 0.012), respectively. After the first hyperthermia washin rate, AUC60, AUC90, and AUCt-max were predictive of improved overall and metastasis-free survival with hazard ratio ranging from 0.46 to 0.53 (P < 0.002) and 0.44 to 0.55 (P < 0.004), respectively. DCE-MRI parameters were compared with extracellular pH and (31)P MR spectroscopy results (previously published) in the same patients showing a correlation. This suggested that an increase in perfusion after therapy was effective in eliminating excess acid from the tumor.

CONCLUSIONS

This study shows that DCE-MRI has utility predicting overall and metastasis-free survival in canine patients with soft tissue sarcomas. To our knowledge, this is the first time that DCE-MRI parameters are predictive of clinical outcome for soft tissue sarcomas.

MUC1 oncoprotein promotes autophagy in a survival response to glucose deprivation

Tumor cells survive under conditions of nutrient deprivation by mechanisms that are not fully understood. The MUC1 oncoprotein is aberrantly overexpressed by most human carcinomas and blocks oxidative stress-induced death. The present studies show that MUC1 inhibits the induction of necrosis in response to the deprivation of glucose. MUC1 suppressed glucose deprivation-induced increases in reactive oxygen species (ROS) and thereby depletion of ATP and cell death. Cells respond to oxidative stress and energy depletion with the induction of autophagy. Our results demonstrate that MUC1 blocks depletion of ATP and sustains growth of glucose-deprived cells by a mechanism sensitive to the autophagy inhibitor, 3-methyladenine. Silencing expression of ATG7, a protein essential for the formation of autophagic vacuoles, also attenuated the MUC1-sustained increases in ATP and growth in response to glucose deprivation. Moreover, we found that MUC1 stimulates AMPK activation and thereby promotes lysosomal turnover of LC3-II, a marker of starvation-induced autophagic activity. These results indicate that MUC1 suppresses glucose deprivation-induced increases in ROS and thereby promotes ATP production and survival. The findings also indicate that the overexpression of MUC1 as found in human cancers could provide a survival advantage in microenvironments with low glucose levels.

Noninvasive multimodality imaging of the tumor microenvironment: registered dynamic magnetic resonance imaging and positron emission tomography studies of a preclinical tumor model of tumor hypoxia

In vivo knowledge of the spatial distribution of viable, necrotic, and hypoxic areas can provide prognostic information about the risk of developing metastases and regional radiation sensitivity and may be used potentially for localized dose escalation in radiation treatment. In this study, multimodality in vivo magnetic resonance imaging (MRI) and positron emission tomography (PET) imaging using stereotactic fiduciary markers in the Dunning R3327-AT prostate tumor were performed, focusing on the relationship between dynamic contrast-enhanced (DCE) MRI using Magnevist (Gd-DTPA) and dynamic (18)F-fluoromisonidazole ((18)F-Fmiso) PET. The noninvasive measurements were verified using tumor tissue sections stained for hematoxylin/eosin and pimonidazole. To further validate the relationship between (18)F-Fmiso and pimonidazole uptake, (18)F digital autoradiography was performed on a selected tumor and compared with the corresponding pimonidazole-stained slices. The comparison of Akep values (kep = rate constant of movement of Gd-DTPA between the interstitial space and plasma and A = amplitude in the two-compartment model (Hoffmann U, Brix G, Knopp MV, Hess T and Lorenz WJ (1995). Magn Reson Med 33, 506-514) derived from DCE-MRI studies and from early (18)F-Fmiso uptake PET studies showed that tumor vasculature is a major determinant of early (18)F-Fmiso uptake. A negative correlation between the spatial map of Akep and the slope map of late (last 1 hour of the dynamic PET scan) (18)F-Fmiso uptake was observed. The relationships between DCE-MRI and hematoxylin/eosin slices and between (18)F-Fmiso PET and pimonidazole slices confirm the validity of MRI/PET measurements to image the tumor microenvironment and to identify regions of tumor necrosis, hypoxia, and well-perfused tissue.

Her2/neu signaling blockade improves tumor oxygenation in a multifactorial fashion in Her2/neu+ tumors

PURPOSE

Tumor hypoxia reduces the efficacy of radiation and chemotherapy as well as altering gene expression that promotes cell survival and metastasis. The growth factor receptor, Her2/neu, is overexpressed in 25-30% of breast tumors. Tumors that are Her2(+) may have an altered state of oxygenation, relative to Her2(-) tumors, due to differences in tumor growth rate and angiogenesis.

METHODS

Her2 blockade was accomplished using an antibody to the receptor (trastuzumab; Herceptin). This study examined the effects of Her2 blockade on tumor angiogenesis, vascular architecture, and hypoxia in Her2(+) and Her2(-) MCF7 xenograft tumors.

RESULTS

Treatment with trastuzumab in Her2(+) tumors significantly improved tumor oxygenation, increased microvessel density, and improved vascular architecture compared with the control-treated Her2(+) tumors. The Her2(+) xenografts treated with trastuzumab also demonstrated decreased proliferation indices when compared with control-treated xenografts. These results indicate that Her2 blockade can improve tumor oxygenation by decreasing oxygen consumption (reducing tumor cell proliferation and inducing necrosis) and increasing oxygen delivery (vascular density and architecture).

CONCLUSIONS

These results support the use of trastuzumab as an adjunct in the treatment of breast tumors with chemotherapy or radiotherapy, as improvements in tumor oxygenation should translate into improved treatment response.

Elevated CAIX Expression is Associated with an Increased Risk of Distant Failure in Early-Stage Cervical Cancer

Tumor hypoxia is associated with adverse outcome in many malignancies. The goal of this study was to determine if elevated expression of carbonic anhydrase IX (CAIX), a biomarker of hypoxia, predicts for recurrence in early-stage cervical cancer. The charts of all patients with early-stage cervical cancer, primarily FIGO IB, treated by radical hysterectomy at our institution from 1988-2001 were reviewed. Adequate pathologic specimens from patients who recurred or who had at least three years follow-up and remained disease-free were stained for CAIX. An immunohistochemical score (IHC) was generated from the extent/intensity of staining. Outcome, as measured by freedom from recurrence (FFR), distant metastases (FFDM) and local recurrence (FFLR), was analyzed as a function of age, IHC, lymph node status (LN) and histology. Forty-two relapsing patients and 76 non-relapsing patients were evaluated. In univariate analysis, +LN, though not IHC or histology, was a significant predictor of any recurrence. Both +LN and higher IHC were associated with decreased FFDM but not FFLR. Patients with both +LN and elevated IHC more frequently exhibited distant metastases as first site of failure (5-year FFDM 50%) than patients with only +LN, elevated IHC or neither feature (70, 85 and 95%, respectively, p = 0.0004). In multivariable analysis, only +LN was significantly associated with poorer FFDM (hazard ratio 4.6, p = 0.0015) though there was a strong trend with elevated CAIX expression (p = 0.069). Elevated CAIX expression is associated with more frequent distant metastases in early-stage cervical cancer, suggesting that patients with this characteristic may benefit from more aggressive treatment.

Elevated CAIX Expression is Associated with an Increased Risk of Distant Failure in Early-Stage Cervical Cancer

Tumor hypoxia is associated with adverse outcome in many malignancies. The goal of this study was to determine if elevated expression of carbonic anhydrase IX (CAIX), a biomarker of hypoxia, predicts for recurrence in early-stage cervical cancer. The charts of all patients with early-stage cervical cancer, primarily FIGO IB, treated by radical hysterectomy at our institution from 1988–2001 were reviewed. Adequate pathologic specimens from patients who recurred or who had at least three years follow-up and remained disease-free were stained for CAIX. An immunohistochemical score (IHC) was generated from the extent/intensity of staining. Outcome, as measured by freedom from recurrence (FFR), distant metastases (FFDM) and local recurrence (FFLR), was analyzed as a function of age, IHC, lymph node status (LN) and histology. Forty-two relapsing patients and 76 non-relapsing patients were evaluated. In univariate analysis, +LN, though not IHC or histology, was a significant predictor of any recurrence. Both +LN and higher IHC were associated with decreased FFDM but not FFLR. Patients with both +LN and elevated IHC more frequently exhibited distant metastases as first site of failure (5-year FFDM 50%) than patients with only +LN, elevated IHC or neither feature (70, 85 and 95%, respectively, p = 0.0004). In multivariable analysis, only +LN was significantly associated with poorer FFDM (hazard ratio 4.6, p = 0.0015) though there was a strong trend with elevated CAIX expression (p = 0.069). Elevated CAIX expression is associated with more frequent distant metastases in early-stage cervical cancer, suggesting that patients with this characteristic may benefit from more aggressive treatment.

Oxygen Dependence and Extravascular Transport of Hypoxia-Activated Prodrugs: Comparison of the Dinitrobenzamide Mustard PR-104A and Tirapazamine

PURPOSE

To compare oxygen dependence and tissue transport properties of a new hypoxia-activated prodrug, PR-104A, with tirapazamine, and to evaluate the implications for antitumor activity when combined with radiotherapy.

METHODS AND MATERIALS

Oxygen dependence of cytotoxicity was measured by clonogenic assay in SiHa cell suspensions. Tissue transport parameters were determined using SiHa multicellular layers. Spatially resolved pharmacokinetic (PK) and pharmacodynamic (PD) models were developed to predict cell killing in SiHa tumors and tested by clonogenic assay 18 h after treatment with the corresponding phosphate ester, PR-104.

RESULTS

The K-value (oxygen concentration to halve cytotoxic potency) of PR-104A was 0.126 +/- 0.021 microM (10-fold lower than tirapazamine at 1.30 +/- 0.28 microM). The diffusion coefficient of PR-104A in multicellular layers (4.42 +/- 0.15 x 10(-7) cm2 s(-1)) was lower than that of tirapazamine (1.30 +/- 0.05 x 10(-6) cm2 s(-1)) but PK modeling predicted better penetration to hypoxic cells in tumors because of its slower metabolism. The tirapazamine PK/PD model successfully predicted the measured activity in combination with single-dose radiation against SiHa tumors, and the PR-104A model underpredicted the activity, which was greater for PR-104 than for tirapazamine (at equivalent host toxicity) both with radiation and as a single agent.

CONCLUSION

PR-104/PR-104A has different PK/PD properties from tirapazamine and superior activity with single-dose radiotherapy against SiHa xenografts. We have inferred that PR-104A is better able to kill cells at intermediate partial pressure of oxygen in tumors than implied by the PK/PD model, most likely because of a bystander effect resulting from diffusion of its activated metabolites from severely hypoxic zones.

Use of three-dimensional tissue cultures to model extravascular transport and predict in vivo activity of hypoxia-targeted anticancer drugs

BACKGROUND

Because of the inefficient vasculature of solid tumors, anticancer drugs must penetrate relatively long distances through the extravascular compartment. The requirement for such diffusion may limit their activity, especially that of hypoxia-targeted drugs. We tested whether a three-dimensional pharmacokinetic/pharmacodynamic (PK/PD) model based on a representative mapped tumor microvascular network could predict the therapeutic activity of anticancer drugs in mouse xenograft tumors.

METHODS

Diffusion coefficients of the hypoxia-activated anticancer drug tirapazamine (TPZ) and of 15 TPZ analogs were estimated by measuring their transport through HT29 colon cancer multicellular layers (MCLs). Anoxic cytotoxic potency (by clonogenic assay) and metabolism of the TPZ analogs were measured in HT29 cell suspensions, and their plasma pharmacokinetics was measured in CD-1 nude mice. This information was used to create a spatially resolved PK/PD model for the tumor microvascular network. Model predictions were compared with actual hypoxic cell kill as measured by clonogenic assays on HT29 xenograft tumors 18 hours after treatment with each TPZ analog.

RESULTS

Modeling TPZ transport in the tumor microvascular network showed substantial drug depletion in the most hypoxic regions, with predicted maximum cell kill of only 3 logs, compared with more than 10 logs if there were no transport impediment. A large range of tissue diffusion coefficients (0.027 x 10(-6)-1.87 x 10(-6) cm2/s) was observed for the TPZ analogs. There was a strong correlation between model-predicted and measured hypoxic cell kill (R2 = 0.89) but a poor correlation when the model did not include extravascular transport (R2 = 0.32).

CONCLUSIONS

Extravascular transport in tumors, and its consequences for tumor cell killing, can be predicted by measuring drug penetration through MCLs in vitro and modeling pharmacokinetics at each position in three-dimensional microvascular networks.

Direct Demonstration of Instabilities in Oxygen Concentrations within the Extravascular Compartment of an Experimental Tumor

To test the hypothesis that temporal variations in microvessel red cell flux cause unstable oxygen levels in tumor interstitium, extravascular oxygenation of R3230Ac mammary tumors grown in skin-fold window chambers was measured using recessed tip polarographic microelectrodes. Red cell fluxes in microvessels surrounding pO2 measurement locations were measured using fluorescently labeled red cells. Temporal pO2 instability was observed in all experiments. Median pO2 was inversely related to radial distance from microvessels. Transient fluctuations above and below 10 mm Hg were consistently seen, except in one experiment near the oxygen diffusion distance limit (140 microm) where pO2 fluctuations were <2 mm Hg and median pO2 was <5 mm Hg. Vascular stasis was not seen in these experiments. These results show that fluctuations in red cell flux, as opposed to vascular stasis, can cause temporal variations in pO2 that extend from perivascular regions to the maximum oxygen diffusion distance.

Assessment of the tumor microenvironment in cervix cancer using dynamic contrast enhanced CT, interstitial fluid pressure and oxygen measurements

PURPOSE

Interstitial fluid pressure (IFP) and oxygen (pO(2)) measurements are prognostic factors in cervical cancer. The purpose of this study was to determine the relationship between IFP and oxygenation and parameters derived from dynamic contrast-enhanced computed tomography (DCE-CT).

METHODS AND MATERIALS

Dynamic contrast-enhanced computed tomography was performed in 32 patients with cervical cancer before radiation therapy. Images were acquired during intravenous contrast injection at 1 per s for 120 s and 1 per 15 s for 60 s. DCE-CT was analyzed using CT Perfusion 3 software (GE Medical Systems) to derive tumor blood flow (BF), permeability surface area product, blood volume, and mean transit time. Further analysis was performed to obtain relative peak enhancement, residual enhancement at 3 min after contrast injection (RE), time to peak and initial slope. Nodal status and tumor size were assessed with MRI. From in vivo IFP (n = 31) and pO(2) (n = 31) tumor measurements median pO(2) (mO(2)), percentage measurements less than 5 mm Hg (HP5) and mean IFP values were calculated.

RESULTS

There was a positive correlation between BF and mO(2) (r = 0.47, p = 0.007) and between RE and HP5 (r = 0.39, p = 0.03). There was no correlation between IFP and DCE-CT parameters.

CONCLUSION

There is a moderately positive, correlation between tumor oxygenation and BF as well as RE and HP5. Further study is required to determine if DCE-CT parameters are useful predictors of tumor behavior in cervical cancer.

Hyperspectral imaging of hemoglobin saturation in tumor microvasculature and tumor hypoxia development

Tumor hypoxia has been shown to have prognostic value in clinical trials involving radiation, chemotherapy, and surgery. Tumor oxygenation studies at microvascular levels can provide understanding of oxygen transport on scales at which oxygen transfer to tissue occurs. To fully grasp the significance of blood oxygen delivery and hypoxia at microvascular levels during tumor growth and angiogenesis, the spatial and temporal relationship of the data must be preserved and mapped. Using tumors grown in window chamber models, hyperspectral imaging can provide serial spatial maps of blood oxygenation in terms of hemoglobin saturation at the microvascular level. We describe our application of hyperspectral imaging for in vivo microvascular tumor oxygen transport studies using red fluorescent protein (RFP) to identify all tumor cells, and hypoxia-driven green fluorescent protein (GFP) to identify the hypoxic fraction. 4T1 mouse mammary carcinoma cells, stably transfected with both reporter genes, are grown in dorsal skin-fold window chambers. Hyperspectral imaging is used to create image maps of hemoglobin saturation, and classify image pixels where RFP alone is present (tumor cells), or both RFP and GFP are present (hypoxic tumor cells). In this work, in vivo calibration of the imaging system is described and in vivo results are shown.

The relationship between hypoxia and angiogenesis

Recent studies have generated a large amount of data supporting the hypothesis that hypoxia drives tumor angiogenesis. The relationship between the two is often considered a matter of supply and demand: ineffectively-vascularized tumor tissue becomes hypoxic, stimulating neoangiogenesis to improve the influx of oxygen, thereby diminishing the angiogenic drive. Although this paradigm is logically pleasing, much of what is known about tumor biology argues against such a straightforward relationship. In fact, some preclinical data convincingly shows that tumor hypoxia and angiogenesis do not always go hand in hand. It is important to begin to explore means of reconciling these discrepancies. Although poor oxygenation is a strong stimulus for tumor angiogenesis, (1) the pathogenesis of tumor hypoxia is much more complicated than the supply-demand paradigm lets on and (2) hypoxia is not necessarily sufficient or necessary for neovascularization to occur. These subtleties may help to explain why so much data disagrees with the current hypoxia-angiogenesis model and may begin to build a better understanding of the role hypoxia plays in tumor vascularization. This article will review what is known about hypoxia and angiogenesis in nononcological processes and will apply these lessons to tumor biology to more deeply describe their relationship.

Why do cancers have high aerobic glycolysis?

If carcinogenesis occurs by somatic evolution, then common components of the cancer phenotype result from active selection and must, therefore, confer a significant growth advantage. A near-universal property of primary and metastatic cancers is upregulation of glycolysis, resulting in increased glucose consumption, which can be observed with clinical tumour imaging. We propose that persistent metabolism of glucose to lactate even in aerobic conditions is an adaptation to intermittent hypoxia in pre-malignant lesions. However, upregulation of glycolysis leads to microenvironmental acidosis requiring evolution to phenotypes resistant to acid-induced cell toxicity. Subsequent cell populations with upregulated glycolysis and acid resistance have a powerful growth advantage, which promotes unconstrained proliferation and invasion.

Tumor-dependent kinetics of partial pressure of oxygen fluctuations during air and oxygen breathing

The primary purpose of this study was to examine the kinetics of partial pressure of oxygen (pO2) fluctuations in fibrosarcoma (FSA) and 9L tumors under air and O2 breathing conditions. The overall hypothesis was that key factors relating to oxygen tension fluctuations would vary between the two tumor types and as a function of the oxygen content of the breathing gas. To assist in the interpretation of the temporal data, spatial pO2 distributions were measured in 10 FSA and 8 9L tumors transplanted into the subcutis of the hind leg of Nembutal-anesthetized (50 mg/kg) Fischer 344 rats. Recessed-tip oxygen microelectrodes were inserted into the tumor, and linear pO2 measurements were recorded in 50-microm steps along a 3-mm path, and blood pressure was simultaneously measured via femoral arterial access. Additionally, pO2 was measured at a single location for 90 to 120 minutes in FSA (n=11) or 9L tumors (n=12). Rats were switched from air to 100% O2 breathing after 45 minutes. Temporal pO2 records were evaluated for their potential radiobiological significance by assessing the number of times they crossed a 10-mm-Hg threshold. In addition, the data were subjected to Fourier analysis for air and O2 breathing. FSA and 9L tumors had spatial median pO2 measurements of 4 and 1 mm Hg, respectively. 9L had more low pO2 measurements < or =2.5 mm Hg than did FSA, whereas between 2.5 and 10 mm Hg this pattern was reversed. Pimonidazole staining patterns in FSA and 9L tumors supported these results. Temporal pO2 instability was observed in all experiments during air and O2 breathing. Threshold analyses indicated that the 10 mm Hg threshold was crossed 2 to 5 times per hour, independent of tumor type. However, the magnitude of 9L pO2 fluctuations was approximately eight times greater than FSA fluctuations, as assessed with Fourier transform analysis (Wilcoxon, P < 0.005). O2 breathing significantly increased median pO2 in FSA from 3 to 8 mm Hg (P < 0.005) and caused a significant increase in frequency and magnitude of pO2 fluctuations. One hundred percent O2 breathing had no effect on 9L tumor pO2, and it decreased the magnitude of pO2 fluctuations with borderline significance. These results show that these two tumors differ significantly with respect to spatial and temporal oxygenation conditions under air and O2 breathing. Fluctuations of pO2 of the type reported herein are predicted to significantly affect radiotherapy response and could be a source for genetic instability, increased angiogenesis, and metastases.

The relationship between the tumor physiologic microenvironment and angiogenesis

This article examines the pathophysiology of tumors, with an emphasis on how these features influence angiogenesis in tumors.

Physiological noise in murine solid tumours using T2*-weighted gradient-echo imaging: a marker of tumour acute hypoxia?

T2*-weighted gradient-echo magnetic resonance imaging (T2*-weighted GRE MRI) was used to investigate spontaneous fluctuations in tumour vasculature non-invasively. FSa fibrosarcomas, implanted intramuscularly (i.m.) in the legs of mice, were imaged at 4.7 T, over a 30 min or 1 h sampling period. On a voxel-by-voxel basis, time courses of signal intensity were analysed using a power spectrum density (PSD) analysis to isolate voxels for which signal changes did not originate from Gaussian white noise or linear drift. Under baseline conditions, the tumours exhibited spontaneous signal fluctuations showing spatial and temporal heterogeneity over the tumour. Statistically significant fluctuations occurred at frequencies ranging from 1 cycle/3 min to 1 cycle/h. The fluctuations were independent of the scanner instabilities. Two categories of signal fluctuations were reported: (i) true fluctuations (TFV), i.e., sequential signal increase and decrease, and (ii) profound drop in signal intensity with no apparent signal recovery (SDV). No temporal correlation between tumour and contralateral muscle fluctuations was observed. Furthermore, treatments aimed at decreasing perfusion-limited hypoxia, such as carbogen combined with nicotinamide and flunarizine, decreased the incidence of tumour T2*-weighted GRE fluctuations. We also tracked dynamic changes in T2* using multiple GRE imaging. Fluctuations of T2* were observed; however, fluctuation maps using PSD analysis could not be generated reliably. An echo-time dependency of the signal fluctuations was observed, which is typical to physiological noise. Finally, at the end of T2*-weighted GRE MRI acquisition, a dynamic contrast-enhanced MRI was performed to characterize the microenvironment in which tumour signal fluctuations occurred in terms of vessel functionality, vascularity and microvascular permeability. Our data showed that TFV were predominantly located in regions with functional vessels, whereas SDV occurred in regions with no contrast enhancement as the result of vessel functional impairment. Furthermore, transient fluctuations appeared to occur preferentially in neoangiogenic hyperpermeable vessels. The present study suggests that spontaneous T2*-weighted GRE fluctuations are very likely to be related to the spontaneous fluctuations in blood flow and oxygenation associated with the pathophysiology of acute hypoxia in tumours. The disadvantage of the T2*-weighted GRE MRI technique is the complexity of signal interpretation with regard to pO2 changes. Compared to established techniques such as intravital microscopy or histological assessments, the major advantage of the MRI technique lies in its capacity to provide simultaneously both temporal and detailed spatial information on spontaneous fluctuations throughout the tumour.

Oxygen dependence of the metabolic activation and cytotoxicity of tirapazamine: implications for extravascular transport and activity in tumors

The hypoxic cytotoxin tirapazamine (TPZ) is currently in phase III clinical trial and appears to have clinical activity. One hypothesis as to why TPZ has been used more successfully in the clinic than most other bioreductive drugs is that its unusual O(2) dependence allows killing of radioresistant cells at "intermediate" O(2) concentrations. We have determined the O(2) dependence of the metabolism of TPZ to its reduction product SR 4317, and its cytotoxicity, in stirred suspensions of HT29 colon carcinoma cells while monitoring O(2) in solution with an Oxylite trade mark probe. The O(2) dependence of the cytotoxicity of TPZ is entirely accounted for by its inhibition of the metabolism of TPZ, with a K(O(2)) value (O(2) concentration for 50% inhibition) of 1.21 +/- 0.09 (SEM) microM. We used this experimental O(2) dependence to extend a recent (Hicks et al., Cancer Res. 63, 5970-5977, 2003) pharmacokinetic/pharmacodynamic model for the cytotoxicity of TPZ in anoxic HT29 multicellular layers to model cell killing in tumors. The model indicates that the O(2) dependence of killing by TPZ complements that of radiation well during fractionated radiotherapy. It predicts that lowering K(O(2)) would decrease killing in radioresistant cells at intermediate O(2) concentrations, while higher K(O(2)) values would exacerbate metabolic consumption of TPZ and thus further impede its penetration into hypoxic regions. Raising K(O(2)) would also increase metabolic activation at physiological O(2) concentrations, thereby compromising hypoxic selectivity. We conclude that the K(O(2)) value of TPZ is indeed close to the optimum for a bioreductive drug of this class (i.e. one that kills only cells in which it is reduced).