Design of Optimized Hypoxia-Activated Prodrugs Using Pharmacokinetic/Pharmacodynamic Modeling

Abstract

Hypoxia contributes to resistance of tumors to some cytotoxic drugs and to radiotherapy, but can in principle be exploited with hypoxia-activated prodrugs (HAP). HAP in clinical development fall into two broad groups. Class I HAP (like the benzotriazine N-oxides tirapazamine and SN30000), are activated under relatively mild hypoxia. In contrast, Class II HAP (such as the nitro compounds PR-104A or TH-302) are maximally activated only under extreme hypoxia, but their active metabolites (effectors) diffuse to cells at intermediate O₂ and thus also eliminate moderately hypoxic cells. Here, we use a spatially resolved pharmacokinetic/pharmacodynamic (SR-PK/PD) model to compare these two strategies and to identify the features required in an optimal Class II HAP. The model uses a Green’s function approach to calculate spatial and longitudinal gradients of O₂, prodrug, and effector concentrations, and resulting killing in a digitized 3D tumor microregion to estimate activity as monotherapy and in combination with radiotherapy. An analogous model for a normal tissue with mild hypoxia and short intervessel distances (based on a cremaster muscle microvessel network) was used to estimate tumor selectivity of cell killing. This showed that Class II HAP offer advantages over Class I including higher tumor selectivity and greater freedom to vary prodrug diffusibility and rate of metabolic activation. The model suggests that the largest gains in class II HAP antitumor activity could be realized by optimizing effector stability and prodrug activation rates. We also use the model to show that diffusion of effector into blood vessels is unlikely to materially increase systemic exposure for realistic tumor burdens and effector clearances. However, we show that the tumor selectivity achievable by hypoxia-dependent prodrug activation alone is limited if dose-limiting normal tissues are even mildly hypoxic.

Introduction

Prodrugs that are enzymatically converted to active metabolites (effectors) within tumors are of interest for selective cancer therapy. Three different strategies have been explored for intra-tumor prodrug activation: (1) xenobiotic metabolizing enzymes that are highly expressed in specific tumors (1, 2); (2) exogenous enzymes targeted to tumors using antibody (3), viral (4), or bacterial (5) vectors; (3) endogenous enzymes that reduce prodrugs only under the hypoxic conditions that prevail in tumors. Hypoxia-activated prodrugs (HAP), which provide this third strategy, not only exploit a generic feature that differentiates tumors from most normal tissues but also potentially overcome the resistance of hypoxic cells to radiotherapy and many chemotherapy drugs (6–15). Several HAP have been evaluated clinically (16–19), including ongoing studies with the nitro compounds TH-302 (20, 21), PR-104A [(22–25); administered as the corresponding phosphate ester PR-104] and PR-610 (26). However, many questions remain as to how HAP should be designed for optimal anticancer activity.

In the present study we utilize pharmacokinetic/pharmacodynamic (PK/PD) models to explore relationships between the reaction/diffusion properties of HAP in the tumor microenvironment and their antitumor activity and selectivity. In this context it is useful to distinguish two broad classes of HAP with different PK/PD features (Figure 1). Class I HAP are activated by reduction under relatively mild hypoxia to generate a reactive cytotoxin, which is restricted to the cell in which it is formed. This class is typified by the benzotriazine N-oxides tirapazamine and SN30000 that undergo 1-electron reduction to short-lived cytotoxic radicals (27, 28), which appear not to affect adjoining cells in three-dimensional (3D) cell cultures (29). The O₂ concentration for 50% inhibition of cytotoxicity () in stirred single cell suspensions is approximately 1 μM for both tirapazamine (30, 31) and SN30000 (32). Class II HAP (typified by nitro compounds such as PR-104A and TH-302) require more severe hypoxia for 1-electron reduction, but the resulting radicals are further reduced to a relatively stable effector that can diffuse from the cell of origin (29, 33–35); this diffusion of effector to nearby cells (which may not themselves be capable of prodrug activation) is here referred to as a bystander effect. In the case of PR-104A, a value of ~0.1 μM has been estimated in stirred single cell suspensions (31), and activation of TH-302 also appears to require severe hypoxia (34) as is the case for other nitro compounds (36–38). Given that half-maximal radiosensitization of tumor cells by oxygen requires concentrations of ~4 μM (39, 40), the ability of Class II HAP to overcome radioresistance may depend on bystander effects.

Figure 1

We have suggested (41) that Class II HAP may be preferable to Class I because activation will be confined to the extreme hypoxia found in tumors, thus minimizing toxicity to normal tissues with mild, physiological hypoxia [e.g., retina (42), liver (43–45), esophagus (46), skin (47, 48), and possibly the bone marrow stem cell niche (49) although the oxygenation status of the latter is controversial (50)]. The bystander effect from Class II HAP may contribute to the reported monotherapy activity of PR-104 (33, 51, 52) and TH-302 (53) in preclinical models. However, it is not known under what conditions this theoretical advantage for Class II HAP might be realized, or how this activity might be optimized by prodrug design, or whether diffusion of active metabolites into the tumor microvasculature might contribute to systemic toxicity if this process is too efficient.

Spatially resolved pharmacokinetic/pharmacodynamic (SR-PK/PD) models provide tools for addressing these questions. These models can be used to describe concentration gradients of oxygen, HAP, and their effectors in tumors, using mapped microvascular networks, and to calculate resulting reproductive cell death (clonogenic cell killing). These models include the effects of heterogeneity in inter-capillary distances, vessel diameters, blood flow rates, and vessel oxygen and drug concentrations. We have validated an SR-PK/PD model for Class I HAP by showing that it predicts activity of tirapazamine analogs combined with radiotherapy in human tumor xenografts (32, 54, 55). This modeling clearly demonstrated the need to optimize rates of reductive metabolism such that penetration into hypoxic regions is not compromised by excessive consumption of the prodrug. Recently, we have also reported an SR-PK/PD model for the Class II HAP, PR-104, and used this to estimate that 30–50% of its activity in HCT116 and SiHa xenografts is due to bystander effects both as monotherapy and combined with radiation (35). Here, we use a generalized SR-PK/PD model, in which a HAP is metabolized by an oxygen-inhibited process to a single effector (Figure 2), to ask under what conditions Class II HAP might provide greater tumor activity and selectivity than Class I HAP, and to identify the prodrug features required for optimal antitumor activity. This generalized HAP model makes explicit the diffusion of both the prodrug and effector in the extracellular (interstitial) compartment. We consider two types of effector, which elicit cytotoxicity via irreversible reaction with a target (Case 1, as for an alkylating agent) or by reversible binding to its target (Case 2, as for a receptor ligand).

Figure 2

Methods

The generalized SR-PK/PD model calculates steady-state concentrations of oxygen, HAP, and effector as well as resulting cell killing in digitized 3D tissue microregions using Green’s function methods (35, 54, 56). We used two different tissue microregions that were derived by mapping microvascular anatomy as well as direction and velocity of blood flow in a rat cremaster muscle (56) (“normal” network) and a subcutaneous FaDu tumor xenograft (57) (“tumor” network). The blood vessels are represented by cylindrical segments and vessel walls are treated as part of the tissue space. The model was implemented using a customized version of the Green’s function method written in Visual C++ (Microsoft Visual Studio 2010 Express) (35, 58).

Calculation of oxygenation

Convective transport of oxygen along vessel segments and diffusion into the surrounding tissue (represented as homogeneous medium) was calculated based on previous estimates for blood flow in the networks, blood O₂ content, tissue diffusion, and consumption of O₂ (35, 56). All O₂ transport parameters are summarized in Table 1.

Calculation of pharmacokinetics

Inflow of prodrug to the microvascular networks was defined by its unbound plasma level, using area under the concentration-time curve (AUC) as a time-independent exposure variable compatible with Green’s function formalism. Vessel walls were modeled as offering negligible resistance to radial flux by setting the intravascular resistance constant K (56) to a low value (0.1 s/μm). Extravascular transport was calculated using a (pro)drug transport model (Figure 2) assuming that the tissue consists of two homogeneous (extracellular and intracellular) compartments, with activation of prodrug restricted to the intracellular compartment and diffusion confined to the extracellular compartment. Concentrations in the two compartments were calculated as follows:

C_e and C_i are the extracellular and intracellular concentrations, respectively, of prodrug or effector (denoted by N = P or E), φ_i and φ_e are the intra- and extracellular volume fractions with φ_e = 1 − φ_i, k_tN are first order rate constants for cell transfer between the extracellular and intracellular compartments, D_N is the effective diffusion coefficient of compound N, ▽² is the Laplacian operator, k_metN is the rate constant for drug metabolism and r is the rate of metabolic production with r_P = 0 and rate of effector production equal to the rate of loss of prodrug: r_E = φ_ik_metPC_iP. The rate constant for prodrug activation is O₂-dependent and is given by: where k_metP,max and k_metP are the prodrug activation rate constants under anoxia and at the O₂ concentration [O₂], respectively, and is the O₂ concentration for half-maximum prodrug activation.

Calculation of cell killing

Survival probability (SP) at each point of the tumor microregion was calculated using the PK/PD relationship: where AUC_iE is the intracellular exposure (area under the intracellular concentration-time curve) to the effector and a is a proportionality constant (potency) equal to the reciprocal of the effector AUC to produce a SP of 0.1 as measured by clonogenic cell survival.

We distinguish two cases: case 1, where potency a scales with k_metE, broadly representing irreversible inhibitors such as DNA-alkylating agents for which an increase in reactivity has similar effects on reaction with target and non-target nucleophiles (59, 60) and case 2, where a is independent of k_metE (broadly representing reversible inhibitors). Case 1, where cell killing is proportional to the reactivity of the effector is used throughout this study unless noted otherwise.

To estimate cell SP after radiation treatment, the linear-quadratic model was used as previously (35, 54). The SP from both radiation and prodrug at each point of the tumor microregion was calculated from the sum of log (SP) due to drug and radiation alone. Log cell kill was then calculated as where the SP is averaged over the whole tumor microregion. Prodrug-induced cell kill in addition to radiation was calculated as the difference between overall log cell kill by prodrug + radiation and log cell kill by radiation alone.

Estimation of net transport across the tumor-plasma boundary

The SR-PK/PD model calculates the transport of prodrug and effector across the walls of the 582 vessel segments in the FaDu tumor microvascular network. Summing the 582 values for each compound gives the net transport across the tumor-plasma boundary (in picomoles). This was normalized to the volume of the FaDu tumor microregion.

Compounds

PR-104H (33) and PR-104M (61) were prepared from PR-104A as described and stored in acetonitrile at −80°C. Tetradeuterated internal standards of PR-104H (PR-104H-d4) and PR-104M (PR-104M-d4) were synthesized using the same methods as for the non-deuterated compounds. PR-104H and PR-104H-d4 had a purity of >90% by HPLC while PR-104M and PR-104M-d4 had a purity of 86 and 84%, respectively.

Pharmacokinetic studies in mice

All animal studies were approved by the University of Auckland Animal Ethics Committee. Male NIH-III nude mice (NIH-Lyst^bgFoxn1^nuBtk^xid; 18–20 g body weight), derived from breeding mice purchased from Charles River Laboratories (Wilmington, MA, USA), were dosed i.v. with PR-104H. The dosing solution was prepared fresh for each mouse by dilution of a stock solution of PR-104H in acetonitrile with sterile saline (to <20% acetonitrile) within 3 min of dosing to minimize loss of PR-104H in aqueous solution. At the dose volume of 0.005 ml/g the acetonitrile dose was 755 mg/kg and the vehicle alone did not cause signs of toxicity over the duration of the experiment. Blood was collected by cardiac puncture up to 3 min after cervical dislocation as in (62). Plasma was prepared as described (63) and mixed with three volumes of cold acidified methanol (methanol:ammonium acetate:acetic acid 1000:3.5:0.2, v/w/v) ~6 min after cardiac puncture. Liver tissue was dissected within 3 min of blood sampling, snap-frozen in liquid nitrogen, and stored at −80°C along with plasma samples. Frozen tissue was pulverized at liquid nitrogen temperature using a BioPulverizer™ (BioSpec Products, USA), transferred to tared tubes, weighed, vortex mixed for 30 s with an equal volume of cold phosphate buffered saline and six volumes of cold acidified methanol, and stored at −80°C. For LC-MS/MS analysis [as reported in Ref. (35)], tissue and plasma extracts were centrifuged (13,000 × g, 4°C, 10 min) and supernatants were mixed 1:1 with cold water containing 1 μM PR-104H-d4 and 1 μM PR-104M-d4.

Results

In order to identify optimal properties of HAP, we developed a SR-PK/PD model that captures the key features of both Class I and Class II HAP (Figure 2). This model is a generalized form of our SR-PK/PD model for PR-104 (35), which, unlike earlier SR-PK/PD models (31, 32, 54) accommodates diffusion and reaction of the effector and explicitly considers intra- and extracellular compartments in order to represent the plasma membrane barrier to prodrug uptake and effector efflux. Our generalized model assumes that the effector can be further converted to inactive products within the cell (defined by rate constant k_metE). The model incorporates only one effector (unlike the two alkylating metabolites from PR-104A) and represents an idealized HAP that is not a substrate for O₂-insensitive two-electron-reduction, which is an off-target activation mechanism of some HAP including PR-104A (51). We also ignore any systemic generation of effector outside the modeled tissue region, such as from hepatic metabolism.

Tissue microregions used for SR-PK/PD modeling

We used a digitized 3D tumor microregion (Figure 3A) that was derived by mapping a region of a subcutaneous FaDu tumor xenograft grown in a mouse dorsal window chamber (57). To assess HAP activation under conditions of physiological hypoxia we also used a 3D region of a rat cremaster muscle (56) as an example of a normal tissue with a well-organized and efficient microvascular network (Figure 3B). While the maximum distance to the nearest blood vessel is 50 μm in the cremaster muscle region, 20% of tissue points in the FaDu tumor region are situated further from blood vessels (50–100 μm; Figure 3C). The FaDu tumor region also showed a higher heterogeneity in vessel diameters (Figure 3D), and a lower median blood flow rate (Figure 3E) at the values chosen for total blood inflow into the regions (Q; Table 1). Q per tissue volume was 2.5-fold lower in the FaDu tumor than in the cremaster muscle microregion, which is realistic for a low-perfused hypoxic tumor area. The value for the tumor microregion falls within the range of mean blood flow values measured in human tumors, which show high variability with values that can be lower, similar, or higher than in the tissue of origin (64, 65). A radiobiological hypoxic fraction (<4 μM O₂) of 43% was achieved in the FaDu tumor microregion (Figure 3F) by using O₂ transport parameters previously chosen to model the hypoxic fraction of tumor xenografts (35) (Table 1). The O₂ transport parameters for the cremaster muscle region were as reported (56), with the oxygen consumption rate adjusted to achieve a radiobiological hypoxic fraction of 7.5% (Figure 3F) that is consistent with moderate hypoxia as measured in some normal human tissues using O₂ electrode histography (45). However, the muscle region lacks the severely hypoxic cells (defined as <0.13 μM O₂, i.e., below the value for PR-104A) that comprise 5.0% of the FaDu microregion (Figure 3F).

Figure 3

Modulation of kinetics of prodrug activation

Our previous SR-PK/PD model for Class I HAP showed that high prodrug activation rates limit killing of hypoxic cells by impairing tissue penetration of the prodrug (54). To investigate whether this also applies to HAP with bystander effects and/or activation restricted to more severe hypoxia (low ), we modulated the prodrug activation rate constant for anoxia (k_metP,max) under these conditions. For simulations without a bystander effect, k_tE was set to 0, which traps effector within the intracellular compartment.

Increasing k_metP,max resulted in a decrease of prodrug penetration to the most hypoxic regions, an effect that was less pronounced with a low (Figure 4A) than with a high value (Figure 4D), as expected. Accordingly, killing in the most hypoxic regions (<0.13 μM O₂) was higher and less affected by an increase in k_metP,max for low- HAP (Figures 4B,E; without bystander effects). The observed shift of cell killing to higher O₂ concentrations with increasing k_metP,max reflects increased prodrug activation at intermediate O₂ concentrations. In case of low- HAP this improved complementation of radiation-induced killing (Figure 4B), but in case of high- HAP it resulted in a loss of hypoxic selectivity at a 100-fold higher k_metP,max (Figure 4E). If the effector was allowed to diffuse through the plasma membrane (i.e., the bystander model), killing was increased across the whole tumor microregion, and killing <0.13 μM was less affected by decreasing prodrug penetration with increasing k_metP,max (Figures 4C,F).

Figure 4

Figures 4C,E illustrate the dependence on prodrug activation rate constant for HAP of Class II (low  + bystander) and Class I (high , no bystander) respectively. It should be noted that in the no-bystander simulations the effector potency, a, was set to a ~fivefold lower value (0.0414 μM⁻¹) than for Class II in Table 2. This was done in order to achieve similar overall (average) killing additional to radiation in the tumor microregion for Class I and II HAP under the default conditions (shown in black). Under these conditions, the Class I HAP provides more killing at intermediate O₂ concentrations (~0.2–10 μM), while the Class II HAP provides much higher cell kill under severe hypoxia (<0.2 μM O₂).

Figures 5A,B show overall HAP-mediated killing, averaged across the whole tumor microregion, for the cases demonstrated in Figure 4. In all cases, monotherapy activity increased with increasing prodrug activation rates over the range examined (k_metP,max = 0.001–1 s⁻¹), although more markedly for low- compounds (shown in red). This difference was even more pronounced for killing of radiobiologically hypoxic cells (i.e., killing additional to radiation; Figure 5B), in which case activity of high- compounds (shown in green) decreases at high rates of metabolism reflecting a severe penetration problem (which is not fully offset by allowing bystander diffusion). Class II HAP (dark red) are thus more tolerant of high values of k_metP,max than Class I HAP (light green).

Figure 5

The above analysis suggests that monotherapy activity increases monotonically with k_metP,max up to very high values, but it is important to consider whether this changes selectivity relative to normal tissues. To address this we used the above cremaster muscle microregion (Figure 3), to represent a generic well-perfused normal tissue with mild physiological hypoxia. For the present purpose, we assumed that this “normal tissue” is populated by cells with the same intrinsic sensitivity to effector as cells in the tumor microregion. Cell kill in the normal tissue region showed an even steeper dependence on k_metP,max than the tumor region (Figure 5C). As a consequence, tumor selectivity (assessed using the tumor:normal tissue ratio of log cell kill) fell as rates of bioreductive metabolism increased (Figure 5D). This falling therapeutic ratio reflects the greater compromise to delivery of rapidly metabolized prodrugs in the tumor network because of lower blood flow rates, longer diffusion distances, and more severe hypoxia than in the normal tissue (Figure 3). The model confirmed our expectation that toxicity to physiologically hypoxic normal tissues will be higher for a high- HAP (Figure 5C).

Modulation of prodrug diffusion properties

We next asked whether the activity of class II HAP is sensitive to the diffusibility of the prodrug, as previously shown by the SR-PK/PD model for tirapazamine analogs (54). In the latter study, tissue diffusion was defined by a single parameter, the effective diffusion coefficient as measured by flux through multicellular layer cultures; this represents the weighted averages of the free drug diffusion coefficient within cells, across membranes and in the extracellular matrix, based on treating tissue as homogenous space. In contrast, the present model explicitly considers extra- and intracellular compartments, and overall diffusion is determined by two parameters: the rate constant for transfer between extra- and intracellular compartments, k_tP, and the diffusion coefficient in the extracellular compartment D_P (Figure 2). Since changes in the physicochemical properties of HAP are expected to affect membrane permeability more than paracellular diffusion, k_tP was modulated. The results (Figure 6) show that hypoxic cell killing for Class I HAP is optimal at k_metP,max ~0.1 s⁻¹ [Figure 6A; consistent with the previous SR-PK/PD model for tirapazamine (54)]. Class I HAP also showed an optimum value of k_tP, reflecting the fact that the membrane transfer rate constant controls how much prodrug is available for intracellular activation. In contrast, for Class II HAP killing additional to radiation increased monotonically over the full range of k_metP,max and k_tP (Figure 6B) reflecting the tolerance of high rates of intracellular prodrug activation owing to a low and because effector diffusion offsets compromised prodrug penetration.

Figure 6

Modulation of the effector diffusion range

Next, we investigated the dependence of antitumor activity on effector parameters for Class II HAP. Since tissue transport of the effector is dependent on its membrane transfer (k_tE), paracellular diffusion (D_E) and its stability (k_metE), these parameters were modulated ~10-fold relative to the base model. In order to see an impact of a change in the paracellular diffusion coefficient D_E, k_tE had to be increased 10-fold to 0.1s⁻¹, indicating that effector transport is otherwise limited by membrane transfer. An increase in diffusion coefficient, D_E (Figure 7A) caused increased diffusion of the effector away from the hypoxic regions where it is formed, which decreased killing in hypoxic regions and increased killing in aerobic regions. As a consequence, overall killing in addition to radiation decreased and monotherapy activity increased (Figures 7B,C). A similar effect could be observed when increasing the membrane transfer rate constants k_tE (Figures 7D,E). However, the effect of increased diffusibility (higher D_E or k_tE) on single-agent activity was only minor because the effector was defined to freely permeate the blood vessel walls, and thus lost into the vasculature. The protective effect of this washout was confirmed by simulations with zero vessel permeability to the effector, which showed a much larger increase in HAP monotherapy activity with increasing k_tE (Figures 7G,H). Notably, an increase in effector diffusibility in the presence of vessel permeability elevated the tumor:normal tissue ratio of killing (Figures 7C,F), reflecting that washout of effector into the circulation is more protective in normal tissue (with lower distances to nearest vessel; Figure 3) than in tumor tissue.

Figure 7

The effect of increasing the rate constant for effector reaction k_metE was investigated using the default case where cell killing scales with effector reactivity (case 1) and alternatively where cell killing is independent of reactivity (case 2). The proportionality constant a was set to the values given in Table 2 to achieve equivalent killing of the two effector types at the default k_metE of 0.01 s⁻¹. In case 1, an increase in k_metE produces a proportional increase in the rate of formation of cytotoxic lesions and resulted in higher killing in hypoxic regions (Figure 8A). Killing in aerobic regions (>4 μM O₂) was decreased at high values of k_metE due to associated low effector penetration, and this limited overall activity (Figure 8B). Tumor selectivity decreased with increasing k_metE (Figure 8C) due largely to increasing cell killing in normal tissue.

Figure 8

In case 2 (where cytotoxic lesion production is independent of effector reactivity) decreasing k_metE resulted in higher killing across the whole tumor microregion (Figure 8D), with a large increase in HAP activity, both as a single agent and in addition to radiation, with increasing effector stability (Figure 8E). Killing in the normal tissue microregion was affected to a similar extent, so that tumor selectivity was comparable for different values of k_metE (Figure 8F).

Estimation of the impact of effector washout by blood flow

The above analysis demonstrated that washout of effector by blood flow can, under some conditions, significantly protect perivascular cells from cytotoxicity. This raises the question whether effector washout might have pharmacodynamic effects in downstream tumor regions (redistribution of effector through the tumor vasculature), or might result in significant systemic exposure to the effector. To test the sensitivity of effector washout on different parameters, we used the Class II HAP model to calculate net transport of effector across the tumor/plasma boundary in each vessel segment of the tumor microregion. This showed high sensitivity to the rate constant of effector production (k_metP,max) and effector stability (k_metE) but only moderate sensitivity to effector diffusibility (k_tE and D_E) (Figure 9), because the former parameters substantially affect overall effector concentrations.

Figure 9

The impact on systemic exposure will depend on the distribution volume (V_D) and clearance (CL) of the effector. To evaluate a specific case, we returned to the SR-PK/PD model for PR-104 (35) because the plasma AUC of its active metabolites (PR-104H and PR-104M) after dosing PR-104 has been determined in three strains of mice (35, 52, 61). To facilitate this analysis, we measured the pharmacokinetics of PR-104H and PR-104M in plasma and liver following i.v. dosing of NIH-III nude mice with synthetic PR-104H (Figure 10). This showed a very short half-life (2.6 min) for PR-104H in plasma, and high PR-104M concentrations in plasma and liver, indicating rapid conversion of PR-104H to PR-104M. The half-lives of PR-104M in plasma and liver (~7 min) were similar to the half-life of PR-104M in culture medium at 37°C [6.6 min (35)], suggesting that CL of PR-104M is mainly due to its chemical instability.

Figure 10

Since washout was most sensitive to k_metP,max (Figure 9), we performed PR-104 SR-PK/PD model simulations using two different values of this parameter, corresponding to rates measured in vitro for HCT116/WT and a cell line with 20-fold higher k_metP,max (HCT116/sPOR#6) (35). This confirmed that washout of metabolites produced in hypoxic regions can increase plasma levels in nearby vessels, an effect that was more pronounced at the higher k_metP,max (white arrows in Figure 11). Using V_D and CL estimates from measured plasma PK following administration of PR-104 or PR-104H (Table 3), we calculated that washout would elevate the free plasma AUC of PR-104H+M in NIH-III nude mice [estimated from the measured PK after i.p. administration of 562 μmol/kg PR-104 (35)] by only 0.8% in case of HCT116/WT tumors, and 3.3% in case of HCT116/sPOR#6 tumors (Table 3). This was calculated assuming all PR-104H+M was released at once and hence represents a maximum addition to the circulating metabolites. This is consistent with published data showing similar plasma concentrations of reduced metabolites in NIH-III nude mice with HCT116/WT and HCT116/sPOR#6 tumors 30 min after i.p. administration of 562 μmol/kg PR-104 (35).

Figure 11

It should be noted that all simulations in this study were performed at blood flows where prodrug delivery and effector extraction were not highly perfusion limited. To demonstrate this (and the effect of blood flow on increasing effector washout into the systemic circulation) we increased the blood flow rate, Q, 2.5-fold, while decreasing inflow pO₂, to maintain a similar hypoxic fraction (Figures S1A,B in Supplementary Material). Under conditions of high prodrug metabolism (k_metP,max of 1 s⁻¹), this increased prodrug extraction from the plasma 1.5-fold and effector washout into the plasma twofold (Figure S1C in Supplementary Material) but overall killing in the microregion increased by only 26% (Figures S1D,E in Supplementary Material) reflecting little increase in perivascular effector concentrations.

Discussion

In this study, SR-PK/PD modeling was utilized to identify strategies for optimization of HAP. The model has a solid biological foundation because it is based on our previous SR-PK/PD model for PR-104 for which parameters have been determined experimentally using three cell lines (35). The latter model was the first to explicitly consider the intra-tumor distribution of bystander metabolites.

Normal tissue toxicity

A novel feature of the present study is the comparison of predictions for a tumor and normal tissue network, the latter based on a mapped microvascular network in a rat cremaster muscle (56) with oxygen consumption rate adjusted to simulate mild hypoxia (Figure 3). This comparison is important because the utility of HAP is likely to be constrained by the requirement that prodrug activation at intermediate O₂ does not cause toxicity in tissues with physiological hypoxia, such as liver (43–45), retina (66), and possibly bone marrow (49). The low tumor:normal tissue ratios of killing determined in our study (Figures 5, 7, and 8) point to limitations on achieving tumor selectivity with HAP. However, it is important to note that our model makes the assumption that intrinsic cellular sensitivity to the released effector (potency parameter a) is the same for cells in the tumor and normal tissue networks. Most effectors of HAP in clinical development are DNA-reactive cytotoxins (13) that provide some tumor-specificity due to the cytokinetics and DNA repair abnormalities (67) of tumor cells relative to normal cells. In addition, release of HAP effectors that target oncogenic pathways in tumor cells has potential for conferring additional tumor selectivity (68). Normal tissue models will need to be revisited when more quantitative information about oxygenation and drug targets in potentially dose-limiting normal tissues is available. Nevertheless, the present results emphasize the importance of using HAP strategies to augment selectivity of effectors that already provide some level of tumor selectivity, rather than as the sole mechanism of tumor targeting.

In addition, normal tissue toxicity due to washout of effectors has been a concern in the development of targeted anticancer prodrugs (69, 70) although to our knowledge no clear link between toxicity and increased effector levels has been established. In the current study systemic toxicity of HAP as a result of release of active metabolites from the tumor appears unlikely; our model suggests that a significant contribution to toxicity would require a combination of extreme assumptions (high tumor burden, intra-tumor activation rates, effector stability, tumor blood flow rates, as well as low systemic CL of effector). However, toxicity in normal tissue sharing a microcirculatory system with the tumor remains a possibility.

Optimal prodrug parameters of class I and class II HAP

This study compared an idealized Class I HAP (high with no-bystander effect, illustrated by tirapazamine and SN30000) and a Class II HAP (low with bystander effect, e.g., PR-104A and TH-302) to assess which of these HAP strategies is better suited for complementation of radiotherapy. Our previous tirapazamine SR-PK/PD model showed that Class I HAP exhibit an optimum k_metP,max, above which hypoxic cell killing decreases because increasing consumption of the prodrug limits its penetration to hypoxic regions (54). This finding could be replicated using the present model, in which intra- and extracellular compartments are made explicit (Figure 5). In contrast, activity of Class II HAP was predicted to increase monotonically with increasing k_metP,max in the investigated range of 0.001–1 s⁻¹; we consider this range to be physiologically and pharmacologically relevant based on measured rate constants (32, 35, 54, 60). In addition, the potential to further increase k_metP,max by drug design would be limited by membrane transport and by the concurrent decrease in tumor selectivity (Figure 5D).

The superiority of Class II over Class I HAP at high k_metP,max can be explained by two factors. Firstly, the lower value minimizes metabolic loss during diffusion into hypoxic target regions (compare Figures 4A,D). Secondly, bystander effects alleviate the impact of poor prodrug penetration due to redistribution of effector between regions of different prodrug exposure. Both of these factors also explain the lower sensitivity of antitumor activity to prodrug diffusibility of class II relative to class I HAP (Figure 6).

Design of optimized class II HAP

Spatially resolved pharmacokinetic/pharmacodynamic modeling represents a valuable tool for the rational design of improved HAP, as has been demonstrated by our previous SR-PK/PD model for Class I HAP (54). The insight that there is an optimum for combinations of D_P and k_metP that can accommodate the competing requirements of efficient prodrug activation and tissue penetration prompted the SR-PK/PD model-guided screening of tirapazamine analogs (71) that identified SN30000 as a prodrug with superior tissue penetration and antitumor activity in combination with radiation (32). The model-based design of class II HAP presents a greater challenge as it requires the specific consideration of the transport parameters of the effector in addition to those of the prodrug. In this study SR-PK/PD modeling was used to identify reaction-diffusion parameters of Class II HAPs that could be optimized. The outcome is summarized in Table 4, which shows that antitumor activity of class II HAP could be enhanced by increasing k_metP,max, although this would need to be balanced against the need to maintain sufficient tumor selectivity.

When looking at the impact of a change in effector metabolism (k_metE) we distinguished two types of effectors: an effector that needs to react irreversibly with its target to elicit cytotoxicity, meaning that cell killing is a function of this chemical reactivity as for alkylating agents (Case 1), and an inhibitor where chemical transformation of the effector (whether spontaneously or by metabolism) gives rise to non-toxic products (Case 2). Modeling showed that effector reactivity should be increased in Case I but decreased in Case 2, in order to improve antitumor activity. In Case 1 the increase in overall killing when increasing reactivity by drug design needs to be balanced against a predicted decrease in tumor selectivity.

The model revealed that bystander effects are limited not only by diffusion and reaction of the effector, but also by effector washout, unless this washout targets downstream regions via a blood-borne bystander effect. The blood vessels act as a sink for effector, a problem that has previously been flagged in the context of drug diffusion into peritoneal tumors from the peritoneal cavity (72). Therefore, the choice of optimal effector diffusibility (defined by k_tE and D_E) depends on the therapeutic context, with low diffusibility generally better for combination with radiation, and high diffusibility preferred for monotherapy activity (Figures 7B,E).

Notably, a greater increase of bystander killing with increasing effector diffusibility has previously been suggested by the finding that bystander effects of dinitrobenzamide mustard prodrugs in multicellular layer co-cultures of NfsB-overexpressing and WT cells correlated with effector lipophilicity (73, 74). This may partially be due to the intimate mixture of activator and target cells in these co-cultures, so that effector washout at the multicellular layer boundaries affects both cell types to the same extent. In contrast, blood vessels in the virtual tumor microregion that act as a sink for effector are generally more distant from hypoxic prodrug-activating regions than from bystander target zones. Therefore, multicellular layer co-cultures have may limited use in assessing antitumor activity for HAP, although they are useful to rank HAP according to the diffusion range of their effectors.

Although an increase in the effector diffusion range provided only a moderate increase in single-agent activity it was found to have a surprisingly large effect in protecting cells in normal tissues where the HAP is activated (Figures 7C,F). Increasing effector diffusibility and prodrug activation rate simultaneously (the latter to compensate for washout) may therefore be a good strategy to improve monotherapy activity.

The suggested HAP model parameters could be modulated in drug design by changing the physicochemical properties that determined these parameters; e.g., one-electron reduction potential in controlling k_metP,max (9, 75, 76) or lipophilicity, pK_a, molecular weight, and number of H-bond donors and acceptors in controlling diffusibility (77). In most cases any change will affect several model parameters, but the value of the SR-PK/PD model lies in part in its ability to accommodate all such changes.

The high prodrug diffusibility and low effector diffusibility that would be ideal for combination with radiotherapy according to the model predictions may not be achievable in all chemical classes of HAP. In case of dinitrobenzamide mustards (PR-104 analogs), the physicochemical properties of prodrug and effector (other than mustard reactivity) will be similar since their structures only differ in one aromatic substituent. The use of fragmenting prodrugs of aliphatic mustards such as TH-302 relaxes this constraint, providing a higher scope for independent optimization of prodrug and effector diffusion properties.

Limitations of the model

The present SR-PK/PD model has two key limitations. Firstly, the model uses a small tumor region (0.12 mm³) that cannot be expected to fully capture the heterogeneity in vascular density and perfusion in tumors (78). This heterogeneity results in macroregional variations in oxygenation and prodrug delivery that may decrease the overall relevance of bystander effects of Class II HAP. For example, macroscopic, well-perfused, and oxygenated areas [observed in some tumor xenografts (35, 79, 80)] would be beyond the reach of the hypoxia-driven bystander effects studied here, unless there is redistribution of diffusible effectors between differently oxygenated tumor regions via the bloodstream (blood-borne bystander effect). Mapped tumor regions of a larger scale are not yet available but could be used to investigate these macroregional effects in the future. Secondly, our steady-state model does not account for time-dependent processes such as cycling hypoxia (81–83) and the kinetics of reoxygenation, cell proliferation, and cell death. Cycling hypoxia is thought to arise from variations in blood flow (84–86), which would affect not only delivery of oxygen but also that of the HAP. Since cycling hypoxia has a dominant periodicity of two to three cycles per hour (82), it may be relevant even for HAP with a residence time in the circulation of few hours and critically important for HAP with long residence times in tumors as is reoxygenation for effectors with long residence times such as AQ4N (87). The implications of variable blood flow and cycling hypoxia should be incorporated in future SR-PK/PD models.

Conclusion

In spite of its limitations, the SR-PK/PD model has identified strategies for optimization of HAP worthy of further investigation. This is the first modeling study that attempts to systematically compare the activity of Class I and Class II HAP and the first that integrates the use of a normal tissue network to assess tumor selectivity. It suggests that class II HAP offer the following advantages relative to class I HAP: they tolerate a wider range of k_metP,max (allowing a wider scope for prodrug design as well as their use in a range of human tumors with different one-electron reductase expression), are less sensitive to prodrug diffusibility (thus making it less critical to optimize this parameter) and have lower toxicity in normal tissues with physiological hypoxia (Figures 4C,F). A possible disadvantage is that complementation of radiation falls off rapidly with increasing effector diffusion away from the hypoxic activating region (Figure 7).

The modeling also showed that the design of optimized Class II HAP is complex due to the need to consider many different processes such as metabolism, membrane transfer, paracellular diffusion, and washout. The spatial and temporal heterogeneity of O₂ and blood flow within tumors (and potentially variations in reductase expression and other determinants of sensitivity) adds further complexity. Therefore it is difficult to design one HAP that fits all tumor types and therapeutic settings. However, important tendencies have been identified by the model, e.g., that antitumor activity may be increased by optimizing effector stability and prodrug activation rates while tumor selectivity may be improved by increasing effector diffusibility.

The present study has implications for targeted anticancer prodrug approaches, as most are limited by spatial heterogeneity in prodrug activation (3, 4, 88). The approach is also potentially applicable to hypoxia-targeting strategies that utilize prodrugs such as gene-directed enzyme prodrug therapy approaches in which the prodrug-activating enzyme is delivered using macrophages that accumulate in hypoxic tumor regions (89), or obligate anaerobic bacteria (90). In approaches with a different spatial distribution of prodrug-activating regions relative to non-activating regions and blood vessels the optimum conditions may differ and could be identified in the future using SR-PK/PD modeling.

Supplementary Material

The Supplementary Material for this article can be found online at http://www.frontiersin.org/Journal/10.3389/fonc.2013.00314/abstract

Impact of blood flow rate on tumor PK/PD of HAP. SR-PK/PD model simulations were performed using the default O₂ transport parameters in Table 1 (Q = 40 nl/min; pO₂ = 40 mm Hg; black) and using a higher blood flow rate and lower inflow pO₂ (Q = 100 nl/min; pO₂ = 29 mm Hg; red). The prodrug activation rate constant k_metP,max was set to a high value of 1 s⁻¹ with all other parameters as in Table 2. (A) O₂ concentration as a function of distance to nearest vessel in the FaDu tumor microregion. (B) The resulting oxygen distribution for the higher blood flow rate case (c.f. Figure 3F). (C) Net amount of prodrug and effector extracted from plasma (− sign) or released into plasma (+ sign). (D) Intracellular prodrug concentration C_iP as a function of O₂ in the FaDu tumor microregion. (E) Resulting killing. Lines indicate average killing in the tumor region (2.47 and 3.04 logs of cell kill for low flow and high flow simulations respectively).

References

• 1

  GraafMMBovenEScheerenHWHaismaHJPinedoHM. Beta-glucuronidase-mediated drug release. Curr Pharm Des (2002) 8:1391–403.10.2174/1381612023394485

  • CrossRef
  • Google Scholar

• 2

  JamiesonDWilsonKPridgeonSMargettsJPEdmondsonRJLeungHYet alNAD (P) H: quinone oxidoreductase 1 and NRH: quinone oxidoreductase 2 activity and expression in bladder and ovarian cancer and lower NRH: quinone oxidoreductase 2 activity associated with an NQO2 exon 3 single-nucleotide polymorphism. Clin Cancer Res (2007) 13:1584.10.1158/1078-0432.CCR-06-1416

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  BagshaweKD. Antibody-directed enzyme prodrug therapy (ADEPT) for cancer. Expert Rev Anticancer Ther (2006) 6:1421–31.10.1586/14737140.6.10.1421

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  DuarteSCarleGFanecaHde LimaMCPierrefite-CarleV. Suicide gene therapy in cancer: where do we stand now?Cancer Lett (2012) 324:160–70.10.1016/j.canlet.2012.05.023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  LehouritisPSpringerCTangneyM. Bacterial-directed enzyme prodrug therapy. J Control Release (2013) 170:120–31.10.1016/j.jconrel.2013.05.005

  • CrossRef
  • Google Scholar

• 6

  BrownJMGiacciaAJ. The unique physiology of solid tumors: opportunities (and problems) for cancer therapy. Cancer Res (1998) 58:1408–16.

  • Pubmed Abstract
  • Google Scholar

• 7

  RauthAMMeloTMisraV. Bioreductive therapies: an overview of drugs and their mechanisms of action. Int J Radiat Oncol Biol Phys (1998) 42:755–62.10.1016/S0360-3016(98)00302-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  StratfordIJWorkmanP. Bioreductive drugs into the next millennium. Anticancer Drug Des (1998) 13:519–28.

  • Google Scholar

• 9

  WardmanP. Electron transfer and oxidative stress as key factors in the design of drugs selectively active in hypoxia. Curr Med Chem (2001) 8:739–61.10.2174/0929867013372959

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  BrownJMWilsonWR. Exploiting tumor hypoxia in cancer treatment. Nat Rev Cancer (2004) 4:437–47.10.1038/nrc1367

  • CrossRef
  • Google Scholar

• 11

  RockwellSDobruckiITKimEYMarrisonSTVuVT. Hypoxia and radiation therapy: past history, ongoing research, and future promise. Curr Mol Med (2009) 9:442–58.10.2174/156652409788167087

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  BrownMHenryS. Kaplan distinguished scientist award lecture 2007. The remarkable yin and yang of tumour hypoxia. Int J Radiat Biol (2010) 86:907–17.10.3109/09553002.2010.492492

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  WilsonWRHayMP. Targeting hypoxia in cancer therapy. Nat Rev Cancer (2011) 11:393–410.10.1038/nrc3064

  • CrossRef
  • Google Scholar

• 14

  WilsonWRHicksKOWangJPruijnFB. Prodrug strategies for targeting tumour hypoxia. In: MelilloG, editor. Hypoxia and Cancer: Biological Implications and Therapeutic Opportunities. New York: Springer (2013). 283 p.

  • Google Scholar

• 15

  HayMPHicksKOWangJ. Hypoxia directed drug strategies to target the tumour microenvironment. In: KoumenisCHammondEGiacciaA, editors. Tumour Microenvironment and Cellular Stress. New York: Springer (2014). 111 p.

  • Google Scholar

• 16

  HafftyBGWilsonLDSonYHChoEIPapacRJFischerDBet alConcurrent chemo-radiotherapy with mitomycin C compared with porfiromycin in squamous cell cancer of the head and neck: final results of a randomized clinical trial. Int J Radiat Oncol Biol Phys (2005) 61:119–28.10.1016/j.ijrobp.2004.07.730

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  PapadopoulosKPGoelSBeeramMWongADesaiKHaigentzMet alA phase 1 open-label, accelerated dose-escalation study of the hypoxia-activated prodrug AQ4N in patients with advanced malignancies. Clin Cancer Res (2008) 14:7110–5.10.1158/1078-0432.CCR-08-0483

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  RischinDPetersLJO’SullivanBGiraltJFisherRYuenKet alTirapazamine, cisplatin, and radiation versus cisplatin and radiation for advanced squamous cell carcinoma of the head and neck (TROG 02.02, HeadSTART): a phase III trial of the Trans-Tasman Radiation Oncology Group. J Clin Oncol (2010) 28:2989–95.10.1200/JCO.2009.27.4449

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  McKeownSRCowenRLWilliamsKJ. Bioreductive drugs: from concept to clinic. Clin Oncol (2007) 19:427–42.10.1016/j.clon.2007.03.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  GanjooKNCranmerLDButrynskiJERushingDAdkinsDOkunoSHet alA phase I study of the safety and pharmacokinetics of the hypoxia-activated prodrug TH-302 in combination with doxorubicin in patients with advanced soft tissue sarcoma. Oncology (2011) 80:50–6.10.1159/000327739

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  BoradMJReddySUronisHSigalDSCohnALSchelmanWRet alRandomized phase II study of the efficacy and safety of gemcitabine + TH-302 (G+T) vs gemcitabine (G) alone in previously untreated patients with advanced pancreatic cancer (Abstract LB-121). AACR Annual Meeting, Chicago (2012).

  • Google Scholar

• 22

  JamesonMBMcKeageMJRamanathanRKRajendranJGuYWilsonWRet alFinal results of a phase Ib trial of PR-104, a pre-prodrug of the bioreductive prodrug PR-104A, in combination with gemcitabine or docetaxel in patients with advanced solid tumors [abstract]. J Clin Oncol (2010) 28:2554.

  • Google Scholar

• 23

  McKeageMJGuYWilsonWRHillAAmiesKMelinkTJet alA phase I trial of PR-104, a pre-prodrug of the bioreductive prodrug PR-104A, given weekly to solid tumour patients. BMC Cancer (2011) 11:432.10.1186/1471-2407-11-432

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  McKeageMJJamesonMBRamanathanRKRajendranJGuYWilsonWRet alPR-104 a bioreductive pre-prodrug combined with gemcitabine or docetaxel in a phase Ib study of patients with advanced solid tumours. BMC Cancer (2012) 12:496.10.1186/1471-2407-12-496

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  Arana YiGYBorthakurGThallPFCovelerALRavandiFJabbourEet alFinal report of phase I/II study of PR104, a hypoxia-activated pro-drug, in relapsed/refractory acute leukemia. J Clin Oncol (2013) 31(Suppl):7074.

  • Google Scholar

• 26

  SmaillJBJaiswalJKAbbattistaMRLuG-LAndersonRFAshoorzadehAet alMechanism of action of the hypoxia-activated irreversible pan-HER inhibitor SN29966. Mol Cancer Ther (2011) 10:A247.10.1158/1535-7163.TARG-11-A247

  • CrossRef
  • Google Scholar

• 27

  YinJGlaserRGatesKS. On the reaction mechanism of tirapazamine reduction chemistry: unimolecular N-OH homolysis, stepwise dehydration, or triazene ring-opening. Chem Res Toxicol (2012) 25:634–45.10.1021/tx200546u

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  ShindeSSMarozAHayMPPattersonAVDennyWAAndersonRF. Characterization of radicals formed following enzymatic reduction of 3-substituted analogues of the hypoxia-selective cytotoxin 3-amino-1,2,4-benzotriazine 1,4-dioxide (tirapazamine). J Am Chem Soc (2010) 132:2591–9.10.1021/ja908689f

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  WilsonWRHicksKOPullenSMFerryDMHelsbyNAPattersonAV. Bystander effects of bioreductive drugs: potential for exploiting pathological tumor hypoxia with dinitrobenzamide mustards. Radiat Res (2007) 167:625–36.10.1667/RR0807.1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  HicksKOSiimBGPruijnFBWilsonWR. Oxygen dependence of the metabolic activation and cytotoxicity of tirapazamine: implications for extravascular transport and activity in tumors. Radiat Res (2004) 161:656–66.10.1667/RR3178

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  HicksKOMyintHPattersonAVPruijnFBSiimBGPatelKet alOxygen dependence and extravascular transport of hypoxia-activated prodrugs: comparison of the dinitrobenzamide mustard PR-104A and tirapazamine. Int J Radiat Oncol Biol Phys (2007) 69:560–71.10.1016/j.ijrobp.2007.05.049

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  HicksKOSiimBGJaiswalJKPruijnFBFraserAMPatelRet alPharmacokinetic/pharmacodynamic modeling identifies SN30000 and SN29751 as tirapazamine analogues with improved tissue penetration and hypoxic cell killing in tumors. Clin Cancer Res (2010) 16:4946–57.10.1158/1078-0432.CCR-10-1439

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  PattersonAVFerryDMEdmundsSJGuYSingletonRSPatelKet alMechanism of action and preclinical antitumor activity of the novel hypoxia-activated DNA crosslinking agent PR-104. Clin Cancer Res (2007) 13:3922–32.10.1158/1078-0432.CCR-07-0478

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  MengFEvansJWBhupathiDBanicaMLanLLorenteGet alMolecular and cellular pharmacology of the hypoxia-activated prodrug TH-302. Mol Cancer Ther (2012) 11:740–51.10.1158/1535-7163.MCT-11-0634

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  FoehrenbacherAPatelKAbbattistaMGuiseCPSecombTWWilsonWRet alThe role of bystander effects in the antitumor activity of the hypoxia-activated prodrug PR-104. Front Oncol (2013) 3:263.10.3389/fonc.2013.00263

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  TaylorYCRauthAM. Oxygen tension, cellular respiration, and redox state as variables influencing the cytotoxicity of the radiosensitizer misonidazole. Radiat Res (1982) 91:104–23.10.2307/3575819

  • CrossRef
  • Google Scholar

• 37

  KochCJ. Unusual oxygen concentration dependence of toxicity of SR-4233, a hypoxic cell toxin. Cancer Res (1993) 53:3992–7.

  • Pubmed Abstract
  • Google Scholar

• 38

  SiimBGAtwellGJWilsonWR. Oxygen dependence of the cytotoxicity and metabolic activation of 4-alkylamino-5-nitroquinoline bioreductive drugs. Br J Cancer (1994) 70:596–603.10.1038/bjc.1994.357

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  WhillansDWHuntJW. A rapid-mixing comparison of the mechanisms of radiosensitization by oxygen and misonidazole in CHO cells. Radiat Res (1982) 90:126–41.10.2307/3575801

  • CrossRef
  • Google Scholar

• 40

  KochCJStobbeCCBumpEA. The effect on the Km for radiosensitization at 0 degree C of thiol depletion by diethylmaleate pretreatment: quantitative differences found using the radiation sensitizing agent misonidazole or oxygen. Radiat Res (1984) 98:141–53.10.2307/3576058

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  DennyWAWilsonWR. Bioreducible mustards: a paradigm for hypoxia-selective prodrugs of diffusible cytotoxins (HPDCs). Cancer Metastasis Rev (1993) 12:135–51.10.1007/BF00689806

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  LeeAEWilsonWR. Hypoxia-dependent retinal toxicity of bioreductive anticancer prodrugs in mice. Toxicol Appl Pharmacol (2000) 163:50–9.10.1006/taap.1999.8834

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  ThurmanRGJiSMatsumuraTLemastersJJ. Is hypoxia involved in the mechanism of alcohol-induced liver injury?Fundam Appl Toxicol (1984) 4:125–33.10.1016/0272-0590(84)90112-X

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  KochCJLordEMShapiroIMClymanRIEvansSM. Imaging hypoxia and blood flow in normal tissues. Adv Exp Med Biol (1997) 428:585–93.10.1007/978-1-4615-5399-1_83

  • CrossRef
  • Google Scholar

• 45

  VaupelPHockelMMayerA. Detection and characterization of tumor hypoxia using pO2 histography. Antioxid Redox Signal (2007) 9:1221–35.10.1089/ars.2007.1628

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  ParliamentMBWiebeLIFrankoAJ. Nitroimidazole adducts as markers for tissue hypoxia: mechanistic studies in aerobic normal tissues and tumour cells. Br J Cancer (1992) 66:1103–8.10.1038/bjc.1992.418

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  StewartFADenekampJRandhawaVS. Skin sensitization by misonidazole: a demonstration of uniform mild hypoxia. Br J Cancer (1982) 45:869–77.10.1038/bjc.1982.139

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  EvansSMSchrlauAEChalianAAZhangPKochCJ. Oxygen levels in normal and previously irradiated human skin as assessed by EF5 binding. J Invest Dermatol (2006) 126:2596–606.10.1038/sj.jid.5700451

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  ParmarKMauchPVergilioJSacksteinRDownJD. Distribution of hematopoietic stem cells in the bone marrow according to regional hypoxia. Proc Natl Acad Sci U S A (2007) 104:5431–6.10.1073/pnas.0701152104

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  Nombela-ArrietaCPivarnikGWinkelBCantyKJHarleyBMahoneyJEet alQuantitative imaging of haematopoietic stem and progenitor cell localization and hypoxic status in the bone marrow microenvironment. Nat Cell Biol (2013) 15:533–43.10.1038/ncb2730

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  GuiseCPAbbattistaMSingletonRSHolfordSDConnollyJDachsGUet alThe bioreductive prodrug PR-104A is activated under aerobic conditions by human aldo-keto reductase 1C3. Cancer Res (2010) 70:1573–84.10.1158/0008-5472.CAN-09-3237

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  HoughtonPJLockRCarolHMortonCLPhelpsDGorlickRet alInitial testing of the hypoxia activated prodrug PR-104 by the Pediatric Preclinical Testing Program. Pediatr Blood Cancer (2011) 57:443–53.10.1002/pbc.22921

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  SunJDLiuQWangJAhluwaliaDFerraroDWangYet alSelective tumor hypoxia targeting by hypoxia-activated prodrug TH-302 inhibits tumor growth in preclinical models of cancer. Clin Cancer Res (2012) 18:758–70.10.1158/1078-0432.CCR-11-1980

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  HicksKOPruijnFBSecombTWHayMPHsuRBrownJMet alUse of three-dimensional tissue cultures to model extravascular transport and predict in vivo activity of hypoxia-targeted anticancer drugs. J Natl Cancer Inst (2006) 98:1118–28.10.1093/jnci/djj306

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  HicksKO. Introducing drug transport early in the design of hypoxia selective anticancer agents using a mathematical modelling approach. In: D’OnofrioACerraiPGandolfiA, editors. New Challenges for Cancer Systems Biomedicine. Milan: Springer (2012). p. 337–54.10.1007/978-88-470-2571-4_18

  • CrossRef
  • Google Scholar

• 56

  SecombTWHsuRParkEYDewhirstMW. Green’s function methods for analysis of oxygen delivery to tissue by microvascular networks. Ann Biomed Eng (2004) 32:1519–29.10.1114/B:ABME.0000049036.08817.44

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  PriesARCornelissenAJSlootAAHinkeldeyMDreherMRHoepfnerMet alStructural adaptation and heterogeneity of normal and tumor microvascular networks. PLoS Comput Biol (2009) 5:e1000394.10.1371/journal.pcbi.1000394

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  SecombTWAlberdingJPHsuRDewhirstMWPriesAR. Angiogenesis: an adaptive dynamic biological patterning problem. PLoS Comput Biol (2013) 9:e1002983.10.1371/journal.pcbi.1002983

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  PanthananickalAHanschCLeoA. Structure-activity relationships in antitumor aniline mustards. J Med Chem (1978) 21:16–26.10.1021/jm00199a004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  HelsbyNAWheelerSJPruijnFBPalmerBDYangSDennyWAet alEffect of nitroreduction on the alkylating reactivity and cytotoxicity of the 2,4-dinitrobenzamide-5-aziridine CB 1954 and the corresponding nitrogen mustard SN 23862: distinct mechanisms of bioreductive activation. Chem Res Toxicol (2003) 16:469–78.10.1021/tx025662b

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  GuYGuiseCPPatelKAbbattistaMRLiJSunXet alReductive metabolism of the dinitrobenzamide mustard anticancer prodrug PR-104 in mice. Cancer Chemother Pharmacol (2011) 67:543–55.10.1007/s00280-010-1354-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  GuYAtwellGJWilsonWR. Metabolism and excretion of the novel bioreductive prodrug PR-104 in mice, rats, dogs and humans. Drug Metab Dispos (2010) 38:498–508.10.1124/dmd.109.030973

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  PatelKLewistonDGuYHicksKOWilsonWR. Analysis of the hypoxia-activated dinitrobenzamide mustard phosphate prodrug PR-104 and its alcohol metabolite PR-104A in plasma and tissues by liquid chromatography-mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci (2007) 856:302–11.10.1016/j.jchromb.2007.06.035

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  VaupelPKallinowskiFOkunieffP. Blood flow, oxygen and nutrient supply, and metabolic microenvironment of human tumors: a review. Cancer Res (1989) 49:6449–65.

  • Pubmed Abstract
  • Google Scholar

• 65

  VaupelP. Tumor microenvironmental physiology and its implications for radiation oncology. Semin Radiat Oncol (2004) 14:198–206.10.1016/j.semradonc.2004.04.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  Wangsa-WirawanNDLinsenmeierRA. Retinal oxygen: fundamental and clinical aspects. Arch Ophthalmol (2003) 121:547–57.10.1001/archopht.121.4.547

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  HelledayTPetermannELundinCHodgsonBSharmaRA. DNA repair pathways as targets for cancer therapy. Nat Rev Cancer (2008) 8:193–204.10.1038/nrc2342

  • CrossRef
  • Google Scholar

• 68

  PattersonAVDonateFJamiesonSThompsonALuGLLeeHHet alPre-clinical characterization of a novel hypoxia-activated irreversible pan-HER inhibitor which induce complete responses in the H1975 xenograft model [abstract]. J Thorac Oncol (2010) 5(6):S249–S50.

  • Google Scholar

• 69

  BagshaweKD. Antibody-directed enzyme/prodrug therapy (ADEPT). Biochem Soc Trans (1990) 18:750–2.

  • Google Scholar

• 70

  BagshaweKDSharmaSKSpringerCJRogersGT. Antibody directed enzyme prodrug therapy (ADEPT). A review of some theoretical, experimental and clinical aspects. Ann Oncol (1994) 5:879–91.

  • Pubmed Abstract
  • Google Scholar

• 71

  HayMPPchalekKPruijnFBHicksKOSiimBGAndersonRFet alHypoxia-selective 3-alkyl 1,2,4-benzotriazine 1,4-dioxides: the influence of hydrogen bond donors on extravascular transport and antitumor activity. J Med Chem (2007) 50:6654–64.10.1021/jm701037w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  DedrickRLFlessnerMF. Pharmacokinetic problems in peritoneal drug administration: tissue penetration and surface exposure. J Natl Cancer Inst (1997) 89:480–7.10.1093/jnci/89.7.480

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  WilsonWRPullenSMHoggAHelsbyNAHicksKODennyWA. Quantitation of bystander effects in nitroreductase suicide gene therapy using three-dimensional cell cultures. Cancer Res (2002) 62:1425–32.

  • Pubmed Abstract
  • Google Scholar

• 74

  AtwellGJYangSPruijnFBPullenSMHoggAPattersonAVet alSynthesis and structure-activity relationships for 2,4-dinitrobenzamide-5 mustards as prodrugs for the Escherichia coli nfsB nitroreductase in gene therapy. J Med Chem (2007) 50:1197–212.10.1021/jm061062o

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  ButlerJHoeyBM. The one-electron reduction potential of several substrates can be related to their reduction rates by cytochrome P-450 reductase. Biochim Biophys Acta (1993) 1161:73–8.10.1016/0167-4838(93)90198-Z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  CenasNAnuseviciusZBironaiteDBachmanovaGIArchakovAIOllingerK. The electron transfer reactions of NADPH: cytochrome P450 reductase with nonphysiological oxidants. Arch Biochem Biophys (1994) 315:400–6.

  • Pubmed Abstract
  • Google Scholar

• 77

  PruijnFBPatelKHayMPWilsonWRHicksKO. Prediction of tumour tissue diffusion coefficients of hypoxia-activated prodrugs from physicochemical parameters. Aust J Chem (2008) 61:687–93.10.1071/CH08240

  • CrossRef
  • Google Scholar

• 78

  GilliesRJSchornackPASecombTWRaghunandN. Causes and effects of heterogeneous perfusion in tumors. Neoplasia (1999) 1:197–207.10.1038/sj.neo.7900037

  • CrossRef
  • Google Scholar

• 79

  Moreno-MerloFNickleeTHedleyDW. Association between tissue hypoxia and elevated non-protein sulphydryl concentrations in human cervical carcinoma xenografts. Br J Cancer (1999) 81:989–93.10.1038/sj.bjc.6690797

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  VukovicVNickleeTHedleyDW. Multiparameter fluorescence mapping of nonprotein sulfhydryl status in relation to blood vessels and hypoxia in cervical carcinoma xenografts 10. Int J Radiat Oncol Biol Phys (2002) 52:837–43.10.1016/S0360-3016(01)02716-X

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  Cardenas-NaviaLIMaceDRichardsonRAWilsonDFShanSDewhirstMW. The pervasive presence of fluctuating oxygenation in tumors. Cancer Res (2008) 68:5812–9.10.1158/0008-5472.CAN-07-6387

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  DewhirstMWCaoYMoellerB. Cycling hypoxia and free radicals regulate angiogenesis and radiotherapy response. Nat Rev Cancer (2008) 8:425–37.10.1038/nrc2397

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  DewhirstMW. Relationships between cycling hypoxia, HIF-1, angiogenesis and oxidative stress. Radiat Res (2009) 172:653–65.10.1667/RR1926.1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  DewhirstMWKimuraHRehmusSWBraunRDPapahadjopoulosDHongKet alMicrovascular studies on the origins of perfusion-limited hypoxia. Br J Cancer Suppl (1996) 27:S247–51.

  • Pubmed Abstract
  • Google Scholar

• 85

  KimuraHBraunRDOngETHsuRSecombTWPapahadjopoulosDet alFluctuations in red cell flux in tumor microvessels can lead to transient hypoxia and reoxygenation in tumor parenchyma. Cancer Res (1996) 56:5522–8.

  • Pubmed Abstract
  • Google Scholar

• 86

  LanzenJBraunRDKlitzmanBBrizelDSecombTWDewhirstMW. Direct demonstration of instabilities in oxygen concentrations within the extravascular compartment of an experimental tumor. Cancer Res (2006) 66:2219–23.10.1158/0008-5472.CAN-03-2958

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87

  PattersonLHMcKeownSR. AQ4N: a new approach to hypoxia-activated cancer chemotherapy. Br J Cancer (2000) 83:1589–93.10.1054/bjoc.2000.1564

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88

  TietzeLFSchmuckK. Prodrugs for targeted tumor therapies: recent developments in ADEPT, GEPT and PMT. Curr Pharm Des (2011) 17:3527–47.10.2174/138161211798194459

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  GriffithsLBinleyKIqballSKanOMaxwellPRatcliffePet alThe macrophage – a novel system to deliver gene therapy to pathological hypoxia. Gene Ther (2000) 7:255–62.10.1038/sj.gt.3301058

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90

  AhnGOBrownM. Targeting tumors with hypoxia-activated cytotoxins. Front Biosci (2007) 12:3483–501.10.2741/2329

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 91

  CarlsonDJKeallPJLooBWJr.ChenZJBrownJM. Hypofractionation results in reduced tumor cell kill compared to conventional fractionation for tumors with regions of hypoxia. Int J Radiat Oncol Biol Phys (2011) 79:1188–95.10.1016/j.ijrobp.2010.10.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92

  PatelKChoySFHicksKOMelinkTJHolfordNHGWilsonWR. A combined pharmacokinetic model for the hypoxia-targeted prodrug PR-104A in humans, dogs, rats and mice predicts species differences in clearance and toxicity. Cancer Chemother Pharmacol (2011) 67:1145–55.10.1007/s00280-010-1412-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar