Recently, 1H, 31P, and 19F NMR magnetic resonance imaging MRI of radiosensitizers has been applied for measuring tumor and tissue oxygenation. In spite of the great importance of such bioreductive drugs, urthermore, surprisingly little detailed computational work on this subject has appeared. Thus, in this work a detailed computational study of BD is presented, calling special attention to the performance of various theoretical methods in reproducing the 13C and 15N, coupling-constant (H-N) data observed in solution. The most sophisticated approach involves density functional based Car–Parrinello molecular dynamics simulations (CPMD) in aqueous solution and averaging chemical shifts over snapshots from the trajectory. In the NMR calculations for these snapshots (performed at the B3LYP level), a small number of discrete water molecules are retained, and the remaining bulk solution effects are included via a polarizable continuum model (PCM). A similarly good accord with experiment is obtained from much less involved, static geometry optimization and NMR computation of pristine 1 employing a PCM approach. Solvent effects on chemical shift are not due to changes in geometric parameters upon solvation, but arise from the direct response of the electronic wave function to the presence of the solvent, which can be represented by discrete molecules and/or the dielectric bulk.