Model systems of the FeMo cofactor of nitrogenase have been explored extensively in catalysis to gain insights into their ability for nitrogen fixation that is of vital importance to the human society. Here we investigate the trigonal pyramidal borane-ligand Fe complex by first-principles calculations, and find that the variation of oxidation state of Fe along the reaction path correlates with that of the reverse-dative Fe → B bonding. The redox-flexibility of the reverse-dative Fe → B bonding helps to provide an electron reservoir that buffers and stabilizes the evolution of Fe oxidation state, which is essential for forming the key intermediates of N2 activation. Our work provides insights for understanding and optimizing homogeneous and surface single-atom catalysts with reverse-dative donating ligands for efficient dinitrogen fixation. The extension of this kind of molecular catalytic active center to heterogeneous catalysts with surface single-clusters is also discussed.