Efficient oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) bifunctional electrocatalysts have been pursued for decades. Meanwhile, single metal atoms embedded in a two-dimensional material substrate (2D-substrate) have emerged as an outstanding catalyst. Herein, we report on computational ORR/OER efficiencies of a series of single atom catalyst systems, with a nitrogen-doped boron phosphide monolayer (N 3 -BP) as the 2D-substrate, and with Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Rh, Pd, Ir, and Pt as the single-atom subject (M). In brief, our density functional theory results show that the overpotentials for ORR/OER are low for CoN 3 -BP, NiN 3 -BP, and PtN 3 -BP, with {η ORR ; η OER } of {0.36; 0.42 V}, {0.29; 0.44 V}, and {0.32; 0.25 V}, respectively. The relevant attributes such as the chemical stability of the 2D-substrate in the ORR/OER environments, immobilization of the single-atom subject on the 2D-substrate, and mechanisms of the ORR/OER activity and the catalyst recovery on the MN 3 -BP catalysts were carefully examined. The key to the comparative study is how the electronic states of the reaction center near the Fermi level of the catalytic system match the frontier orbitals of ORR/OER reaction intermediates. In short, our method predicts the ORR/OER catalytic efficiencies of novel catalysts via a single-atom/2D-substrate design strategy.