Abstract
The Ni-R active site in [NiFe]-hydrogenase features a bridging hydride between the Ni and Fe, displaced toward the Ni. However, all synthetic Ni-R models reported to date exhibit a hydride displaced toward Fe and display low turnover frequencies for H(2) evolution. Understanding the factors governing the hydride position and activity of Ni-R and biomimetic complexes is crucial for developing efficient hydrogen-evolving catalysts. By utilizing the CCSD theory, DFT, NBO, and QTAIM analysis, we investigated these factors in a Ni-R active-site model (1), and two representative biomimetic complexes, 2* and 3. Our results reveal that the Ni site of 1 inherently prefers a square-planar [S(2)NiSH] configuration with an apically positioned thiolate and that hydride positioning is governed by the strength of [Ni-H-Fe] three-center two-electron bonding, which is modulated by the geometric torsion between the Ni terminal ligands and the bridging thiolates. By modifying the linkers between the Ni terminal ligands and bridging thiolate ligands of 2* and 3, we designed virtual biomimetic complexes (4-10). These complexes exhibit improved hydride nucleophilicity and increased potential for H(2) formation, providing valuable insights into how geometric and electronic factors influence hydride activity and informing the design of more effective biomimetic hydrogenase models.