Organophosphorous ligands in hydrogenase-inspired iron-based catalysts: A DFT study on the energetics of metal protonation as a function of P-atom substitution

How general is the effect of the bulkiness of organic ligands on the basicity of metal-organic catalysts? H2-evolving Mo oxides/sulphides as case studies

Tailoring the activity of an organometallic catalyst usually requires a targeted ligand design. Tuning the ligand bulkiness and tuning the electronic properties are popular approaches, which are somehow interdependent because substituents of different sizes within ligands can determine inter alia the occurrence of different degrees of inductive effects. Ligand basicity, in particular, turned out to be a key property for the modulation of protonation reactions occurring in vacuo at the metals in complexes bearing organophosphorus ligands; however, when the same reactions take place in a polar organic solvent, their energetics becomes dependent on the trade-off between ligand basicity and bulkiness, with the polarity of the solvent playing a key role in this regard [Bancroft et al., Inorg. Chem., 1986, 25, 3675; Rovaletti et al., J. Phys. Org. Chem., 2018, 31, e3748]. In the present contribution, we carried out molecular dynamics and density functional theory calculations on water-soluble Mo-based catalysts for proton reduction, in order to study the energetics of protonation reactions in complexes where the incipient proton binds a catalytically active ligand (i.e., an oxide or a disulphide). We considered complexes either soaked in water or in a vacuum, and featuring N-based ancillary ligands of different bulkiness (i.e. cages constituted either by pyridine or isoquinoline moieties). Our results show that the energetics of protonation events can be affected by ancillary ligand bulkiness even when the metal center does not play the role of the H⁺ acceptor. In vacuo, protonation at the O or S atom in the α position relative to the metal in complexes featuring the bulky isoquinoline-based ligand is more favored by around 10 kcal mol^-1 when compared to the case of the pyridine-based counterparts, a difference that is almost zero when the same reactions occur in water. Such an outcome is rationalized in light of the different electrostatic properties of complexes bearing ancillary ligands of different sizes. The overall picture from theory indicates that such effects of ligand bulkiness can be relevant for the design of green chemistry catalysts that undergo protonation steps in water solutions.

Catalytic H2 evolution/oxidation in [FeFe]-hydrogenase biomimetics: account from DFT on the interplay of related issues and proposed solutions

A DFT overview on selected issues regarding diiron catalysts related to [FeFe]-hydrogenase biomimetic research, with implications for both energy conversion and storage strategies.

Protonation and electrochemical properties of a bisphosphide diiron hexacarbonyl complex bearing amino groups on the phosphide bridge

A bisphosphide-bridged diiron hexacarbonyl complex 3 with NEt2 groups on the phosphide bridge was synthesized to examine a new proton relay system from the NEt2 group to the bridging hydride between the two iron centers. As a precursor of the bridging moiety, peri-Et2NP-PNEt2-bridged naphthylene 5 was synthesized by the reaction of 1,8-dilithionaphthylene with two equivalents of Cl2PNEt2 followed by reductive P-P bond formation by magnesium. The reaction of the diphosphine ligand 5 with Fe2(CO)9 gave the diiron hexacarbonyl complex 3, in which the P-P bond of the ligand was cleaved to form the bisphosphide-bridge. The molecular structure of 3 indicated that the trigonal plane of the NEt2 group was forced to face the Fe-Fe bond to avoid steric congestion with the naphthylene group linking the two phosphide groups. The NEt2 group could be protonated by p-toluenesulfonic acid. Density functional theory (DFT) calculations confirmed that the proton of the N(H)Et2 group adopted a position close to the bridging hydride. The DFT results for the ferrocene analogue 1, in which the 1,8-naphthylene group of 3 was replaced with the 1,1'-ferrocenylene group, also revealed that the most stable orientation of the protonated NHEt2 group was that in the protonated 3. As a result, electrochemical proton reduction reactions using complexes 1 and 3 proceeded with similar catalytic efficiencies. Unfortunately, the catalytic efficiencies (CEs) of these complexes were much lower than those of the complexes with a proton relay system of the terminal hydrogen, indicating that the reactive properties of the bridging hydride in the present proton relay system cannot exceed those of the terminal hydride.

The Challenging in silico Description of Carbon Monoxide Oxidation as Catalyzed by Molybdenum-Copper CO Dehydrogenase

Carbon monoxide (CO) is a highly toxic gas to many living organisms. However, some microorganisms are able to use this molecule as the sole source of carbon and energy. Soil bacteria such as the aerobic Oligotropha carboxidovorans are responsible for the annual removal of about 2x108 tons of CO from the atmosphere. Detoxification through oxidation of CO to CO2 is enabled by the MoCu-dependent CO-dehydrogenase enzyme (MoCu-CODH) which-differently from other enzyme classes with similar function-retains its catalytic activity in the presence of atmospheric O2. In the last few years, targeted advancements have been described in the field of bioengineering and biomimetics, which is functional for future technological exploitation of the catalytic properties of MoCu-CODH and for the reproduction of its reactivity in synthetic complexes. Notably, a growing interest for the quantum chemical investigation of this enzyme has recently also emerged. This mini-review compiles the current knowledge of the MoCu-CODH catalytic cycle, with a specific focus on the outcomes of theoretical studies on this enzyme class. Rather controversial aspects from different theoretical studies will be highlighted, thus illustrating the challenges posed by this system as far as the application of density functional theory and hybrid quantum-classical methods are concerned.