Design, Synthesis and Anticancer Evaluation of Nitroimidazole Radiosensitisers.

The role of hypoxic tumour cells in resistance to radiotherapy, and in suppression of immune response, continues to endorse tumour hypoxia as a bona fide, yet largely untapped, drug target. Radiotherapy innovations such as stereotactic body radiotherapy herald new opportunities for classical oxygen-mimetic radiosensitisers. Only nimorazole is used clinically as a radiosensitiser, and there is a dearth of new radiosensitisers in development. In this report, we augment previous work to present new nitroimidazole alkylsulfonamides and we document their cytotoxicity and ability to radiosensitise anoxic tumour cells in vitro. We compare radiosensitisation with etanidazole and earlier nitroimidazole sulfonamide analogues and we identify 2-nitroimidazole and 5-nitroimidazole analogues with marked tumour radiosensitisation in ex vivo assays of surviving clonogens and with in vivo tumour growth inhibition.

Protein-metabolite interactomics of carbohydrate metabolism reveal regulation of lactate dehydrogenase.

Metabolic networks are interconnected and influence diverse cellular processes. The protein-metabolite interactions that mediate these networks are frequently low affinity and challenging to systematically discover. We developed mass spectrometry integrated with equilibrium dialysis for the discovery of allostery systematically (MIDAS) to identify such interactions. Analysis of 33 enzymes from human carbohydrate metabolism identified 830 protein-metabolite interactions, including known regulators, substrates, and products as well as previously unreported interactions. We functionally validated a subset of interactions, including the isoform-specific inhibition of lactate dehydrogenase by long-chain acyl-coenzyme A. Cell treatment with fatty acids caused a loss of pyruvate-lactate interconversion dependent on lactate dehydrogenase isoform expression. These protein-metabolite interactions may contribute to the dynamic, tissue-specific metabolic flexibility that enables growth and survival in an ever-changing nutrient environment.

Glyoxylate protects against cyanide toxicity through metabolic modulation.

Although cyanide's biological effects are pleiotropic, its most obvious effects are as a metabolic poison. Cyanide potently inhibits cytochrome c oxidase and potentially other metabolic enzymes, thereby unleashing a cascade of metabolic perturbations that are believed to cause lethality. From systematic screens of human metabolites using a zebrafish model of cyanide toxicity, we have identified the TCA-derived small molecule glyoxylate as a potential cyanide countermeasure. Following cyanide exposure, treatment with glyoxylate in both mammalian and non-mammalian animal models confers resistance to cyanide toxicity with greater efficacy and faster kinetics than known cyanide scavengers. Glyoxylate-mediated cyanide resistance is accompanied by rapid pyruvate consumption without an accompanying increase in lactate concentration. Lactate dehydrogenase is required for this effect which distinguishes the mechanism of glyoxylate rescue as distinct from countermeasures based solely on chemical cyanide scavenging. Our metabolic data together support the hypothesis that glyoxylate confers survival at least in part by reversing the cyanide-induced redox imbalances in the cytosol and mitochondria. The data presented herein represent the identification of a potential cyanide countermeasure operating through a novel mechanism of metabolic modulation.

Tissue Pharmacokinetic Properties and Bystander Potential of Hypoxia-Activated Prodrug CP-506 by Agent-Based Modelling.

Hypoxia-activated prodrugs are bioactivated in oxygen-deficient tumour regions and represent a novel strategy to exploit this pharmacological sanctuary for therapeutic gain. The approach relies on the selective metabolism of the prodrug under pathological hypoxia to generate active metabolites with the potential to diffuse throughout the tumour microenvironment and potentiate cell killing by means of a "bystander effect". In the present study, we investigate the pharmacological properties of the nitrogen mustard prodrug CP-506 in tumour tissues using in silico spatially-resolved pharmacokinetic/pharmacodynamic (SR-PK/PD) modelling. The approach employs a number of experimental model systems to define parameters for the cellular uptake, metabolism and diffusion of both the prodrug and its metabolites. The model predicts rapid uptake of CP-506 to high intracellular concentrations with its long plasma half-life driving tissue diffusion to a penetration depth of 190 µm, deep within hypoxic activating regions. While bioreductive metabolism is restricted to regions of severe pathological hypoxia (<1 µM O2), its active metabolites show substantial bystander potential with release from the cell of origin into the extracellular space. Model predictions of bystander efficiency were validated using spheroid co-cultures, where the clonogenic killing of metabolically defective "target" cells increased with the proportion of metabolically competent "activator" cells. Our simulations predict a striking bystander efficiency at tissue-like densities with the bis-chloro-mustard amine metabolite (CP-506M-Cl2) identified as a major diffusible metabolite. Overall, this study shows that CP-506 has favourable pharmacological properties in tumour tissue and supports its ongoing development for use in the treatment of patients with advanced solid malignancies.

Radiosensitisation of SCCVII tumours and normal tissues in mice by the DNA-dependent protein kinase inhibitor AZD7648.

BACKGROUND AND PURPOSE: Inhibitors of DNA-dependent protein kinase (DNA-PK) are effective radiation sensitisers in preclinical tumours, but little is known about risks of normal tissue radiosensitisation. Here, we evaluate radiosensitisation of head and neck squamous cell carcinoma (HNSCC) cells by DNA-PK inhibitor AZD7648 under oxia and anoxia in vitro, and tumour (SCCVII), oral mucosa and small intestine in mice. MATERIALS AND METHODS: Radiosensitisation of human (UT-SCC-54C) and murine (SCCVII) HNSCC cells by AZD7648 under oxia and anoxia was evaluated by clonogenic assay. Radiosensitisation of SCCVII tumours in C3H mice by oral AZD7648 (75 mg/kg) was determined by ex vivo clonogenic assay 3.5 days post-irradiation, with evaluation of normal tissue surrogate endpoints using 5-ethynyl-2'-deoxyuridine to facilitate detection of regenerating crypts in the ileum and repopulating S-phase cells in the ileum and oral mucosa of the same animals. RESULTS: AZD7648 potently radiosensitised both cell lines, with similar sensitiser enhancement ratios for 10% survival (SER10) under oxia and anoxia. AZD7648 diffused rapidly through multicellular layers, suggesting rapid equilibration between plasma and hypoxic zones in tumours. SCCVII tumours were radiosensitised by AZD7648 (SER10 2.5). AZD7648 also enhanced radiation-induced body weight loss and suppressed regenerating intestinal crypts and repopulating S-phase cells in the ileum and tongue epithelium with SER values similar to SCCVII tumours. CONCLUSION: AZD7648 is a potent radiation sensitiser of both oxic and anoxic tumour cells, but also markedly radiosensitises stem cells in the small intestine and oral mucosa.

Restoring Tumour Selectivity of the Bioreductive Prodrug PR-104 by Developing an Analogue Resistant to Aerobic Metabolism by Human Aldo-Keto Reductase 1C3.

PR-104 is a phosphate ester pre-prodrug that is converted in vivo to its cognate alcohol, PR-104A, a latent alkylator which forms potent cytotoxins upon bioreduction. Hypoxia selectivity results from one-electron nitro reduction of PR-104A, in which cytochrome P450 oxidoreductase (POR) plays an important role. However, PR-104A also undergoes 'off-target' two-electron reduction by human aldo-keto reductase 1C3 (AKR1C3), resulting in activation in oxygenated tissues. AKR1C3 expression in human myeloid progenitor cells probably accounts for the dose-limiting myelotoxicity of PR-104 documented in clinical trials, resulting in human PR-104A plasma exposure levels 3.4- to 9.6-fold lower than can be achieved in murine models. Structure-based design to eliminate AKR1C3 activation thus represents a strategy for restoring the therapeutic window of this class of agent in humans. Here, we identified SN29176, a PR-104A analogue resistant to human AKR1C3 activation. SN29176 retains hypoxia selectivity in vitro with aerobic/hypoxic IC50 ratios of 9 to 145, remains a substrate for POR and triggers γH2AX induction and cell cycle arrest in a comparable manner to PR-104A. SN35141, the soluble phosphate pre-prodrug of SN29176, exhibited superior hypoxic tumour log cell kill (>4.0) to PR-104 (2.5-3.7) in vivo at doses predicted to be achievable in humans. Orthologues of human AKR1C3 from mouse, rat and dog were incapable of reducing PR-104A, thus identifying an underlying cause for the discrepancy in PR-104 tolerance in pre-clinical models versus humans. In contrast, the macaque AKR1C3 gene orthologue was able to metabolise PR-104A, indicating that this species may be suitable for evaluating the toxicokinetics of PR-104 analogues for clinical development. We confirmed that SN29176 was not a substrate for AKR1C3 orthologues across all four pre-clinical species, demonstrating that this prodrug analogue class is suitable for further development. Based on these findings, a prodrug candidate was subsequently identified for clinical trials.

Selectively Targeting Tumor Hypoxia With the Hypoxia-Activated Prodrug CP-506.

Hypoxia-activated prodrugs (HAP) are a promising class of antineoplastic agents that can selectively eliminate hypoxic tumor cells. This study evaluates the hypoxia-selectivity and antitumor activity of CP-506, a DNA alkylating HAP with favorable pharmacologic properties. Stoichiometry of reduction, one-electron affinity, and back-oxidation rate of CP-506 were characterized by fast-reaction radiolytic methods with observed parameters fulfilling requirements for oxygen-sensitive bioactivation. Net reduction, metabolism, and cytotoxicity of CP-506 were maximally inhibited at oxygen concentrations above 1 µmol/L (0.1% O2). CP-506 demonstrated cytotoxicity selectively in hypoxic 2D and 3D cell cultures with normoxic/anoxic IC50 ratios up to 203. Complete resistance to aerobic (two-electron) metabolism by aldo-keto reductase 1C3 was confirmed through gain-of-function studies while retention of hypoxic (one-electron) bioactivation by various diflavin oxidoreductases was also demonstrated. In vivo, the antitumor effects of CP-506 were selective for hypoxic tumor cells and causally related to tumor oxygenation. CP-506 effectively decreased the hypoxic fraction and inhibited growth of a wide range of hypoxic xenografts. A multivariate regression analysis revealed baseline tumor hypoxia and in vitro sensitivity to CP-506 were significantly correlated with treatment response. Our results demonstrate that CP-506 selectively targets hypoxic tumor cells and has broad antitumor activity. Our data indicate that tumor hypoxia and cellular sensitivity to CP-506 are strong determinants of the antitumor effects of CP-506.

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.

Sugar phosphate activation of the stress sensor eIF2B.

The multi-subunit translation initiation factor eIF2B is a control node for protein synthesis. eIF2B activity is canonically modulated through stress-responsive phosphorylation of its substrate eIF2. The eIF2B regulatory subcomplex is evolutionarily related to sugar-metabolizing enzymes, but the biological relevance of this relationship was unknown. To identify natural ligands that might regulate eIF2B, we conduct unbiased binding- and activity-based screens followed by structural studies. We find that sugar phosphates occupy the ancestral catalytic site in the eIF2Bα subunit, promote eIF2B holoenzyme formation and enhance enzymatic activity towards eIF2. A mutant in the eIF2Bα ligand pocket that causes Vanishing White Matter disease fails to engage and is not stimulated by sugar phosphates. These data underscore the importance of allosteric metabolite modulation for proper eIF2B function. We propose that eIF2B evolved to couple nutrient status via sugar phosphate sensing with the rate of protein synthesis, one of the most energetically costly cellular processes.

Identification of small molecule allosteric modulators of 5,10-methylenetetrahydrofolate reductase (MTHFR) by targeting its unique regulatory domain.

The folate and methionine cycles, constituting one-carbon metabolism, are critical pathways for cell survival. Intersecting these two cycles, 5,10-methylenetetrahydrofolate reductase (MTHFR) directs one-carbon units from the folate to methionine cycle, to be exclusively used for methionine and S-adenosylmethionine (AdoMet) synthesis. MTHFR deficiency and upregulation result in diverse disease states, rendering it an attractive drug target. The activity of MTHFR is inhibited by the binding of AdoMet to an allosteric regulatory domain distal to the enzyme's active site, which we have previously identified to constitute a novel fold with a druggable pocket. Here, we screened 162 AdoMet mimetics using differential scanning fluorimetry, and identified 4 compounds that stabilized this regulatory domain. Three compounds were sinefungin analogues, closely related to AdoMet and S-adenosylhomocysteine (AdoHcy). The strongest thermal stabilisation was provided by (S)-SKI-72, a potent inhibitor originally developed for protein arginine methyltransferase 4 (PRMT4). Using surface plasmon resonance, we confirmed that (S)-SKI-72 binds MTHFR via its allosteric domain with nanomolar affinity. Assay of MTHFR activity in the presence of (S)-SKI-72 demonstrates inhibition of purified enzyme with sub-micromolar potency and endogenous MTHFR from HEK293 cell lysate in the low micromolar range, both of which are lower than AdoMet. Nevertheless, unlike AdoMet, (S)-SKI-72 is unable to completely abolish MTHFR activity, even at very high concentrations. Combining binding assays, kinetic characterization and compound docking, this work indicates the regulatory domain of MTHFR can be targeted by small molecules and presents (S)-SKI-72 as an excellent candidate for development of MTHFR inhibitors.

Subcellular Location of Tirapazamine Reduction Dramatically Affects Aerobic but Not Anoxic Cytotoxicity.

Hypoxia is an adverse prognostic feature of solid cancers that may be overcome with hypoxia-activated prodrugs (HAPs). Tirapazamine (TPZ) is a HAP which has undergone extensive clinical evaluation in this context and stimulated development of optimized analogues. However the subcellular localization of the oxidoreductases responsible for mediating TPZ-dependent DNA damage remains unclear. Some studies conclude only nuclear-localized oxidoreductases can give rise to radical-mediated DNA damage and thus cytotoxicity, whereas others identify a broader role for endoplasmic reticulum and cytosolic oxidoreductases, indicating the subcellular location of TPZ radical formation is not a critical requirement for DNA damage. To explore this question in intact cells we engineered MDA-231 breast cancer cells to express the TPZ reductase human NADPH: cytochrome P450 oxidoreductase (POR) harboring various subcellular localization sequences to guide this flavoenzyme to the nucleus, endoplasmic reticulum, cytosol or inner surface of the plasma membrane. We show that all POR variants are functional, with differences in rates of metabolism reflecting enzyme expression levels rather than intracellular TPZ concentration gradients. Under anoxic conditions, POR expression in all subcellular compartments increased the sensitivity of the cells to TPZ, but with a fall in cytotoxicity per unit of metabolism (termed 'metabolic efficiency') when POR is expressed further from the nucleus. However, under aerobic conditions a much larger increase in cytotoxicity was observed when POR was directed to the nucleus, indicating very high metabolic efficiency. Consequently, nuclear metabolism results in collapse of hypoxic selectivity of TPZ, which was further magnified to the point of reversing O2 dependence (oxic > hypoxic sensitivity) by employing a DNA-affinic TPZ analogue. This aerobic hypersensitivity phenotype was partially rescued by cellular copper depletion, suggesting the possible involvement of Fenton-like chemistry in generating short-range effects mediated by the hydroxyl radical. In addition, the data suggest that under aerobic conditions reoxidation strictly limits the TPZ radical diffusion range resulting in site-specific cytotoxicity. Collectively these novel findings challenge the purported role of intra-nuclear reductases in orchestrating the hypoxia selectivity of TPZ.

Crystal structure and interaction studies of human DHTKD1 provide insight into a mitochondrial megacomplex in lysine catabolism.

DHTKD1 is a lesser-studied E1 enzyme among the family of 2-oxoacid de-hydrogenases. In complex with E2 (di-hydro-lipo-amide succinyltransferase, DLST) and E3 (dihydrolipo-amide de-hydrogenase, DLD) components, DHTKD1 is involved in lysine and tryptophan catabolism by catalysing the oxidative de-carboxyl-ation of 2-oxoadipate (2OA) in mitochondria. Here, the 1.9 Šresolution crystal structure of human DHTKD1 is solved in complex with the thi-amine diphosphate co-factor. The structure reveals how the DHTKD1 active site is modelled upon the well characterized homologue 2-oxoglutarate (2OG) de-hydrogenase but engineered specifically to accommodate its preference for the longer substrate of 2OA over 2OG. A 4.7 Šresolution reconstruction of the human DLST catalytic core is also generated by single-particle electron microscopy, revealing a 24-mer cubic scaffold for assembling DHTKD1 and DLD protomers into a megacomplex. It is further demonstrated that missense DHTKD1 variants causing the inborn error of 2-amino-adipic and 2-oxoadipic aciduria impact on the complex formation, either directly by disrupting the interaction with DLST, or indirectly through destabilizing the DHTKD1 protein. This study provides the starting framework for developing DHTKD1 modulators to probe the intricate mitochondrial energy metabolism.

Optical coherence tomography complements confocal microscopy for investigation of multicellular tumour spheroids.

Knowledge of optical properties, such as the refractive index (RI), of biological tissues is important in optical imaging, as they influence the distribution and propagation of light in tissue. To accurately study the response of cancerous cells to drugs, optimised imaging protocols are required. This study uses a simple custom-built spectral domain optical coherence tomography (OCT) system to conduct RI measurements of multicellular spheroids, three-dimensional (3D) in-vitro culture systems, of the cell line HCT116. The spheroid RIs are compared to study the effect of growth over time. To improve confocal microscopy imaging protocols, two immersion media (glycerol and ScaleView-A2) matching the spheroid RIs were trialled, with the aim to reduce the RI mismatch between the spheroid and the immersion medium and thus improve imaging depth with confocal microscopy. ScaleView-A2 (n = 1.380) aided in achieving greater depths of imaging of the multicellular spheroids under confocal microscopy. This improvement in imaging depth confirmed the utility of our RI measurements, proving the promising outlook of OCT as a complementary tool to microscopy in cancer research.

Benzotriazine Di-Oxide Prodrugs for Exploiting Hypoxia and Low Extracellular pH in Tumors.

Extracellular acidification is an important feature of tumor microenvironments but has yet to be successfully exploited in cancer therapy. The reversal of the pH gradient across the plasma membrane in cells that regulate intracellular pH (pHi) has potential to drive the selective uptake of weak acids at low extracellular pH (pHe). Here, we investigate the dual targeting of low pHe and hypoxia, another key feature of tumor microenvironments. We prepared eight bioreductive prodrugs based on the benzotriazine di-oxide (BTO) nucleus by appending alkanoic or aminoalkanoic acid sidechains. The BTO acids showed modest selectivity for both low pHe (pH 6.5 versus 7.4, ratios 2 to 5-fold) and anoxia (ratios 2 to 8-fold) in SiHa and FaDu cell cultures. Related neutral BTOs were not selective for acidosis, but had greater cytotoxic potency and hypoxic selectivity than the BTO acids. Investigation of the uptake and metabolism of representative BTO acids confirmed enhanced uptake at low pHe, but lower intracellular concentrations than expected for passive diffusion. Further, the modulation of intracellular reductase activity and competition by the cell-excluded electron acceptor WST-1 suggests that the majority of metabolic reductions of BTO acids occur at the cell surface, compromising the engagement of the resulting free radicals with intracellular targets. Thus, the present study provides support for designing bioreductive prodrugs that exploit pH-dependent partitioning, suggesting, however, that that the approach should be applied to prodrugs with obligate intracellular activation.

Schedule-dependent potentiation of chemotherapy drugs by the hypoxia-activated prodrug SN30000.

Hypoxia-activated prodrugs (HAPs) are hypothesized to improve the therapeutic index of chemotherapy drugs that are ineffective against tumor cells in hypoxic microenvironments. SN30000 (CEN-209) is a benzotriazine di-N-oxide HAP that potentiates radiotherapy in preclinical models, but its combination with chemotherapy has not been explored. Here we apply multiple models (monolayers, multicellular spheroids and tumor xenografts) to identify promising SN30000/chemotherapy combinations (with chemotherapy drugs before, during or after SN30000 exposure). SN30000, unlike doxorubicin, cisplatin, gemcitabine or paclitaxel, was more active against cells in spheroids than monolayers by clonogenic assay. Combinations of SN30000 and chemotherapy drugs in HCT116/GFP and SiHa spheroids demonstrated hypoxia-and schedule-dependent potentiation of gemcitabine or doxorubicin in growth inhibition and clonogenic assays. Co-administration with SN30000 suppressed clearance of gemcitabine in NIH-III mice, likely due to SN30000-induced hypothermia which also modulated extravascular transport of gemcitabine in tumor tissue as assessed from its diffusion through HCT116 multicellular layer cultures. Despite these systemic effects, the same schedules that gave therapeutic synergy in spheroids (SN30000 3 h before or during gemcitabine, but not gemcitabine 3 h before SN30000) enhanced growth delay of HCT116 xenografts without increasing host toxicity. Identification of hypoxic and S-phase cells by immunohistochemistry and flow cytometry established that hypoxic cells initially spared by gemcitabine subsequently reoxygenate and re-enter the cell cycle, and that this repopulation is prevented by SN30000 only when administered with or before gemcitabine. This illustrates the value of spheroids in modeling tumor microenvironment-dependent drug interactions, and the potential of HAPs for overcoming hypoxia-mediated drug resistance.

Development of a Genomic Data Flow Framework: Results of a Survey Administered to NIH-NHGRI IGNITE and eMERGE Consortia Participants.

Precision health's more individualized molecular approach will enrich our understanding of disease etiology and patient outcomes. Universal implementation of precision health will not be feasible, however, until there is much greater automation of processes related to genomic data transmission, transformation, and interpretation. In this paper, we describe a framework for genomic data flow developed by the Clinical Informatics Work Group of the NIH National Human Genome Research Institute (NHGRI) IGNITE Network consortium. We subsequently report the results of a genomic data flow survey administered to sites funded by NIH-NHGRI for large scale genomic medicine implementations. Finally, we discuss insights and challenges identified through these survey results as they relate to both the current and a desirable future state of genomic data flow.

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 O2, 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.

Bystander Effects of Hypoxia-Activated Prodrugs: Agent-Based Modeling Using Three Dimensional Cell Cultures.

Intra-tumor heterogeneity represents a major barrier to anti-cancer therapies. One strategy to minimize this limitation relies on bystander effects via diffusion of cytotoxins from targeted cells. Hypoxia-activated prodrugs (HAPs) have the potential to exploit hypoxia in this way, but robust methods for measuring bystander effects are lacking. The objective of this study is to develop experimental models (monolayer, multilayer, and multicellular spheroid co-cultures) comprising 'activator' cells with high expression of prodrug-activating reductases and reductase-deficient 'target' cells, and to couple these with agent-based models (ABMs) that describe diffusion and reaction of prodrugs and their active metabolites, and killing probability for each cell. HCT116 cells were engineered as activators by overexpressing P450 oxidoreductase (POR) and as targets by knockout of POR, with fluorescent protein and antibiotic resistance markers to enable their quantitation in co-cultures. We investigated two HAPs with very different pharmacology: SN30000 is metabolized to DNA-breaking free radicals under hypoxia, while the dinitrobenzamide PR104A generates DNA-crosslinking nitrogen mustard metabolites. In anoxic spheroid co-cultures, increasing the proportion of activator cells decreased killing of both activators and targets by SN30000. An ABM parameterized by measuring SN30000 cytotoxicity in monolayers and diffusion-reaction in multilayers accurately predicted SN30000 activity in spheroids, demonstrating the lack of bystander effects and that rapid metabolic consumption of SN30000 inhibited prodrug penetration. In contrast, killing of targets by PR104A in anoxic spheroids was markedly increased by activators, demonstrating that a bystander effect more than compensates any penetration limitation. However, the ABM based on the well-studied hydroxylamine and amine metabolites of PR104A did not fit the cell survival data, indicating a need to reassess its cellular pharmacology. Characterization of extracellular metabolites of PR104A in anoxic cultures identified more stable, lipophilic, activated dichloro mustards with greater tissue diffusion distances. Including these metabolites explicitly in the ABM provided a good description of activator and target cell killing by PR104A in spheroids. This study represents the most direct demonstration of a hypoxic bystander effect for PR104A to date, and demonstrates the power of combining mathematical modeling of pharmacokinetics/pharmacodynamics with multicellular culture models to dissect bystander effects of targeted drug carriers.

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.