Nearly 30% of women with early-stage breast cancer develop recurrent disease attributed to resistance to systemic therapy. Prevailing models of chemotherapy failure describe three resistant phenotypes: cells with alterations in transmembrane drug transport, increased detoxification and repair pathways, and alterations leading to failure of apoptosis. Proliferative activity correlates with tumor sensitivity. Cell-cycle status, controlling proliferation, depends on local concentration of oxygen and nutrients. Although physiologic resistance due to diffusion gradients of these substances and drugs is a recognized phenomenon, it has been difficult to quantify its role with any accuracy that can be exploited clinically. We implement a mathematical model of tumor drug response that hypothesizes specific functional relationships linking tumor growth and regression to the underlying phenotype. The model incorporates the effects of local drug, oxygen, and nutrient concentrations within the three-dimensional tumor volume, and includes the experimentally observed resistant phenotypes of individual cells. We conclude that this integrative method, tightly coupling computational modeling with biological data, enhances the value of knowledge gained from current pharmacokinetic measurements, and, further, that such an approach could predict resistance based on specific tumor properties and thus improve treatment outcome. [Cancer Res 2009;69(10):4484-92]