The phosphomimetic mutation S82D in the cancer-associated, FAD-dependent human NADP(H):quinone oxidoreductase 1 (hNQO1) causes a decrease in flavin-adenine dinucleotide-binding affinity and intracellular stability. We test in this work whether the evolutionarily recent neutral mutation R80H in the vicinity of S82 may alter the strong functional effects of S82 phosphorylation through electrostatic interactions. We show using biophysical and bioinformatic analyses that the reverse mutation H80R prevents the effects of S82D phosphorylation on hNQO1 by modulating the local stability. Consistently, in rat NQO1 (rNQO1) which contains R80, the effects of phosphorylation were milder, resembling the behaviour found in hNQO1 when this residue was humanized in rNQO1 (by the R80H mutation). Thus, apparently neutral and evolutionarily divergent mutations may determine the functional response of mammalian orthologues towards phosphorylation.