Doxorubicin (DOXO) is a well-established chemotherapy drug for treatment of different tumors, ranging from breast cancer, melanoma to multiple myeloma (MM). Here, we present a coupled experimental/modeling approach to study DOXO pharmacokinetics in MM cells, investigate its distribution among the extracellular and intracellular compartments during time. Three model candidates are considered and identified. Model selection is performed based on its ability to describe the data both qualitatively and in terms of quantitative indexes. The most parsimonious model consists of a nonlinear structure with a saturation-threshold control of intracellular DOXO efflux by the DOXO bound to the cellular DNA. This structure could explain the hypothesis that MM cells are drug-resistant, likely due to the involvement of P-glycoproteins.The proposed model is able to predict the intracellular (free and bound) DOXO and suggests the presence of a saturation-threshold drug-resistant mechanism.Clinical Relevance- The model can be used to properly understand and guide further experimental setup, e.g., to investigate multiple myeloma cell variability among different cell lines.