Then we expanded the scope of investigation to the seconcary coordination shell of a metal center by employing a dual active-center model (shorted as ho-2-MN4/G), which have been demonstrated to disrupt linear relationships in NRR . However, simply introducing a second adjacent metal center just did not bring SCE that can change existing scaling relationships in OER. As shown in Figure As shown in Figure , the adsorption of OH on selected metal centers (inside indigo oval) is much stronger, which makes it thermodynamically unfavorable to be converted into *O of OER because of large ordinate value (∆E*O-∆E*OH). Consequently, OH may stay on these oxygen-philic metal centers, making them potential secondary coordination ligands for the adjacent metal center, similar to the secondary coordination ligands in homogeneous catalysis . Then we selected 4 types of M1N4 (M1= Cd, Sc, Y, and Zn) to host OH based on the relative strength of OH adsorption on each M1N4, and 6 types of M2N4 (M2=Ag, Co, Cu, Ir, Ni, and Rh) as potential active centers for OER. A total of 24 combinations of M1N4-M2N4/G, of which the thermo stability has been verified by a 20-ps ab initio molecular dynamics simulation, were generated and the corresponding adsorption energies of *OH, *O, and *OOH were calculated using the same procedure. In the presence of *OH on M1 sites (denoted as *OH(M1)), the linear scaling of adsorption free energies between *O and *OH (Figure ) on M2 sites is substantially weakened (R 2 decreased from 0.93 in MN4/G to 0.54 in M1N4-M2N4/G).