Unprecedented scatter plots of calculated versus measured NMR J coupling constants in six densely oxygen functionalized epoxides are found with some B3LYP protocols, an effect attributed to stereoelectronic effects. Hence, 26 other exchange-correlation density functionals (xc DFs) are benchmarked in this work. Very good results are achieved with mPW1PW91 and PBE0 in conjunction with the pcJ-1 basis set (BS) of moderate size. A thorough statistical analysis of 53 relationships between the predicted and observed J datasets is presented. The effects of some xc DFs, including their x and c parts, and BSs on the calculation results are discussed, also in the context of DFT modeling of electron-density distributions. Moreover, related J datasets predicted with 11 different DF methods are considered and compared with the experimental data. Finally, some proposals for further improvement of existing DFs based on the available J (n=1-3) values are briefly outlined, in line with recent results on the DFT electron densities.