In this combined computational and experimental study, specific chemical interactions affecting the prediction of one-electron and two-electron reduction potentials for anthraquinone derivatives are investigated. For 19 redox reactions in acidic aqueous solution, where AQ is reduced to hydroanthraquinone, density functional theory (DFT) with the polarizable continuum model (PCM) gives a mean absolute deviation (MAD) of 0.037 V for 16 species. DFT(PCM), however, highly overestimates three redox couples with a MAD of 0.194 V, which is almost 5 times that of the remaining 16. These three molecules have ether groups positioned for intramolecular hydrogen bonding that are not balanced with the intermolecular H-bonding of the solvent. This imbalanced description is corrected by quantum mechanics/molecular mechanics (QM/MM) simulations, which include explicit water molecules. The best theoretical estimations result in a good correlation with experiments, V(Theory) = 0.903V(Expt) + 0.007 with an R 2 value of 0.835 and an MAD of 0.033 V. In addition to the aqueous test set, 221 anthraquinone redox couples in aprotic solvent were studied. Five anthraquinone derivatives spanning a range of redox potentials were selected from this library, and their reduction potentials were measured by cyclic voltammetry. DFT(PCM) calculations predict the first reduction potential with high accuracy giving the linear relation, V(Theory) = 0.960V(Expt) − 0.049 with an R 2 value of 0.937 and an MAD of 0.051 V. This approach, however, significantly underestimates the second reduction potential, with an MAD of 0.329 V. It is shown herein that treatment of explicit ion-pair interactions between the anthraquinone derivatives and the cation of the supporting electrolyte is required for the accurate prediction of the second reduction potential. After the correction, V(Theory) = 1.045V(Expt) − 0.088 with an R 2 value 0.910 and an MAD value reduced by more than half to 0.145 V. Finally, molecular design principles are discussed that go beyond simple electron-donating and electron-withdrawing effects to lead to predictable and controllable reduction potentials.