Design of optimized hypoxia-activated prodrugs using pharmacokinetic/pharmacodynamic modeling

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