For the design of an efficient drug delivery system utilizing an ionic liquid (IL) as a carrier, it is prudent to gain molecular/atomistic level insights of a drug with IL in terms of binding and dynamics. In this scenario, the influence of anionic counterpart of imidazolium-based ILs, namely, 1-butyl-3-methyl-imidazolium octyl sulfate [BMIM][OSU] and 1-butyl-3-methyl-imidazolium chloride [BMIM][Cl] in their submicellar region ([IL] = 20 mM) on the model water-soluble anticancer drug doxorubicin hydrochloride (DOX) was probed by employing an arsenal of nuclear magnetic resonance (NMR) approaches. The salient feature of the present study includes the significant interaction of DOX with [BMIM][OSU], whereas the lack of such an interaction with [BMIM][Cl] is gauged by 1H NMR translation self-diffusometry and is further corroborated by 13C chemical shift perturbation. The two-step model was utilized to estimate the bound fraction (pb) and equivalent partition coefficient (K) of DOX with [BMIM][OSU]. A combination of selective and nonselective spin-lattice relaxation rates (R1SEL and R1NS, respectively) enables to gauze the significant interaction of DOX with [BMIM][OSU] over [BMIM][Cl]. Furthermore, 1D transient and truncated driven nuclear Overhauser enhancement (NOE) data analyses in the initial rate limit permits the evaluation of the cross-relaxation efficacy of DOX with the investigated ILs. An Arrhenius-type temperature dependence of the drug's self-diffusion was observed for DOX, DOX-[BMIM][OSU], and DOX-[BMIM][Cl] aqueous mixtures and the corresponding activation energies were evaluated.